Genome-wide identification of microRNAs in larch and stage-specific modulation of 11 conserved microRNAs and their targets during somatic embryogenesis

MicroRNAs (miRNAs) are emerging as essential regulators of biological processes. Somatic embryogenesis is one of the most important techniques for gymnosperm-breeding programs, but there is little understanding of its underlying mechanism. To investigate the roles of miRNAs during somatic embryogenesis in larch, we constructed a small RNA library from somatic embryos. High-throughput sequencing of the library identified 83 conserved miRNAs from 35 families, 16 novel miRNAs, and 14 plausible miRNA candidates, with a high proportion specific to larch or gymnosperms. qRT-PCR analysis demonstrated that both the conserved and novel or candidate miRNAs were expressed in larch. Several miRNA precursor sequences were obtained via RACE. We predicted 110 target genes using bioinformatics, and validated 9 of them by 5' RACE. 11 conserved miRNA families including 17 miRNAs with critical functions in plant development and six target mRNAs were detected by qRT-PCR in the larch SE. Stage-specific expression of miRNAs and their targets indicate their possible modulation on SE of larch: miR171a/b might exert function on PEMs, while miR171c acts in the induction process of larch SE; miR397 and miR398 mainly involved in modulation of PEM propagation and transition to single embryo; miR162 and miR168 exert their regulatory function during total SE process, especially during stages 5-8; miR156, miR159, miR160, miR166, miR167, and miR390 might play regulatory roles during cotyledonary embryo development. These findings indicate that larch and possibly other gymnosperms have complex mechanisms of gene regulation involving specific and common miRNAs operating post-transcriptionally during embryogenesis.