BackgroundMicroRNAs (miRNAs) regulate gene expression by mediating gene silencing at transcriptional and post-transcriptional levels in higher plants. miRNAs and related target genes have been widely studied in model plants such as Arabidopsis and rice; however, the number of identified miRNAs in soybean (Glycine max) is limited, and global identification of the related miRNA targets has not been reported in previous research.ResultsIn our study, a small RNA library and a degradome library were constructed from developing soybean seeds for deep sequencing. We identified 26 new miRNAs in soybean by bioinformatic analysis and further confirmed their expression by stem-loop RT-PCR. The miRNA star sequences of 38 known miRNAs and 8 new miRNAs were also discovered, providing additional evidence for the existence of miRNAs. Through degradome sequencing, 145 and 25 genes were identified as targets of annotated miRNAs and new miRNAs, respectively. GO analysis indicated that many of the identified miRNA targets may function in soybean seed development. Additionally, a soybean homolog of Arabidopsis SUPPRESSOR OF GENE SLIENCING 3 (AtSGS3) was detected as a target of the newly identified miRNA Soy_25, suggesting the presence of feedback control of miRNA biogenesis.ConclusionsWe have identified large numbers of miRNAs and their related target genes through deep sequencing of a small RNA library and a degradome library. Our study provides more information about the regulatory network of miRNAs in soybean and advances our understanding of miRNA functions during seed development.

Cytosine methylation is an important mechanism for dynamical regulation of gene expression and transposable element (TE) mobility during plant developmental processes. Here, we identified the transcription start sites of genes using high-throughput sequencing and then analyzed the DNA methylation status in soybean roots, stems, leaves, and cotyledons of developing seeds at single-base resolution. Profiling of DNA methylation in different organs revealed 2162 differentially methylated regions among organs, and a portion of hypomethylated regions were correlated with high expression of neighboring genes. Because of the different distribution of class I TEs (retrotransposons) and class II TEs (DNA transposons), the promoters of the lowest-expressed genes showed higher levels of CG and CHG methylation but a lower level of CHH methylation. We further found that the CHH methylation level of class II TEs was higher than class I TEs, possibly due to the presence of more smRNAs in class II TEs. In cotyledons of developing seeds, smRNA abundance was roughly positively correlated with hypermethylated regions but negatively related to hypomethylated regions. These studies provide significant insights into the complicated interplays among DNA methylation, smRNA abundance, TE distribution, and gene expression in soybean.