Title: Hypermethylation of miRNA Genes During Nodule Development.

DNA methylation has recently emerged as a powerful regulatory mechanism controlling the expression of key regulators of various developmental processes, including nodulation. However, the functional role of DNA methylation in regulating the expression of microRNA (miRNA) genes during the formation and development of nitrogen-fixing nodules remains largely unknown. In this study, we profiled DNA methylation patterns of miRNA genes during nodule formation, development, and early senescence stages in soybean () through the analysis of methylC-seq data. Absolute DNA methylation levels in the CG, CHH, and CHH sequence contexts over the promoter and primary transcript regions of miRNA genes were significantly higher in the nodules compared with the corresponding root tissues at these three distinct nodule developmental stages. We identified a total of 82 differentially methylated miRNAs in the nodules compared with roots. Differential DNA methylation of these 82 miRNAs was detected only in the promoter (69), primary transcript region (3), and both in the promoter and primary transcript regions (10). The large majority of these differentially methylated miRNAs were hypermethylated in nodules compared with the corresponding root tissues and were found mainly in the CHH context and showed stage-specific methylation patterns. Differentially methylated regions in the promoters of 25 miRNAs overlapped with transposable elements, a finding that may explain the vulnerability of miRNAs to DNA methylation changes during nodule development. Gene expression analysis of a set of promoter-differentially methylated miRNAs pointed to a negative association between DNA methylation and miRNA expression. Gene Ontology and pathways analyses indicate that changes in DNA methylation of miRNA genes are reprogrammed and contribute to nodule development through indirect regulation of genes involved in cellular processes and pathways with well-established roles in nodulation.