Genome-wide Mining for Micrornas and Their Targets in Betula Luminifera Using High-Throughput Sequencing and Degradome Analyses
Junhong Zhang,Menghui Huang,Junsheng Liang,Ying Pan,Longjun Cheng,Jun Wu,Zaikang Tong
DOI: https://doi.org/10.1007/s11295-016-1047-2
2016-01-01
Tree Genetics & Genomes
Abstract:microRNAs (miRNAs) play vital roles in plants regulating a panoply of biological processes, such as development, hormone signaling, and stress response, by inhibiting target genes at the post-transcriptional level. However, the roles of miRNAs in Betula luminifera remain elusive. To mine for B. luminifera miRNAs and targets, we used a deep sequencing approach to analyze the sRNAs and degradome sequencing of mixed samples, including roots, stems, and leaves. A total of 114 known miRNAs or miRNA*s from 44 families, and 24 novel miRNAs and 17 miRNA*s plus 15 plausible miRNA candidates were identified, of which 36 known miRNAs, 29 miRNA*s, and all novel miRNAs had precursor sequences. Additionally, 49 targets for 19 known miRNA families and seven miRNA*s, and seven targets for novel miRNAs were identified using a high-throughput degradome-sequencing approach. The conserved miRNA targets were mainly transcription factors, whereas the miRNA* targets were mainly protein-coding genes, with preferential propensity to functional enzymes. A Gene Ontology analysis showed that the predicted targets were classified into 62 biological processes, 20 cellular components, and 28 molecular functions, respectively. We found two different targets for miR396a* and miR396c*, and the target changed when the miR156d precursor cleavage site was shifted toward the 5′-end by two nucleotides, indicating the diverse regulatory roles of MIRNA genes. Furthermore, three targets identified by degradome sequencing were validated further through 5′ rapid amplification of cDNA ends. The expression patterns of the randomly selected miRNAs varied among different tissues. miR164 expression was induced under nitrogen starvation, with tissue-specific expression patterns, and was negatively correlated with the NAC1 target gene in roots and leaves, but not in stems. This study is a transcriptome-based analysis of miRNAs and the degradome, providing useful information to explore the functions of miRNAs and their target genes in B. luminifera.