Cytoplasmic male sterility (CMS) plays an important role in the production of soybean hybrid seeds. MicroRNAs (miRNAs) are a class of non-coding endogenous ~21 nt small RNAs that play crucial roles in flower and pollen development by targeting genes in plants. Here, two small RNA libraries and two degradome libraries were constructed from the flower buds of the soybean CMS line NJCMS1A and its restorer (Rf) line NJCMS1C. Following high-throughput sequencing, 558 known miRNAs, 103 novel miRNAs on the other arm of known pre-miRNAs, 10 novel miRNAs, and a number of base-edited miRNAs were identified. Among the identified miRNAs, 76 differentially expressed miRNAs were discovered with greater than two-fold changes between NJCMS1A and NJCMS1C. By degradome analysis, a total of 466 distinct transcripts targeted by 200 miRNAs and 122 distinct transcripts targeted by 307 base-edited miRNAs were detected. Further integrated analysis of transcriptome and small RNA found some miRNAs and their targets’ expression patterns showing a negative correlation, such as miR156b-GmSPL and miR4413b-GmPPR. Previous reports showed that these targets might be related to flower bud development, suggesting that miRNAs might act as regulators of soybean CMS fertility. These findings may provide a better understanding of the miRNA-mediated regulatory networks in CMS mechanisms of soybean.

Exploration of miRNA-mediated fertility regulation network of cytoplasmic male sterility during flower bud development in soybean