Identification of milRNAs and their target genes in Ganoderma lucidum by high-throughput sequencing and degradome analysis

• Published in: Fungal genetics and biology 2020-03, Vol.136, p.103313-103313, Article 103313
• Main Authors: Shao, Junjie, Wang, Liqiang, Liu, Yang, Qi, Qianru, Wang, Bin, Lu, Shanfa, Liu, Chang
• Format: Article
• Language: English
• Subjects: Degradome sequencing; Ganoderma lucidum; milRNAs; Small RNA sequencing; Stem-loop RT-qPCR
• ISSN: 1087-1845; 1096-0937
• DOI: 10.1016/j.fgb.2019.103313

•For the first time, Identied the milRNAs and their targets experimentally in G. lucidum.•Identified many miRNA features that may be unique to basidiomycetes.•The milRNAs’ functions are developmental stages-dependent.•Different mechanism of regulation for miRNA might be present.•Conserved miRNA biosynthesis and processing machineries are present in basdiomycetes. MicroRNAs (miRNAs in animals and plants or milRNAs in fungi) are endogenous noncoding RNAs that can regulate gene expression. However, little information is known about milRNAs and their target genes in Ganoderma lucidum. Here, we systematically predicted and characterised the milRNAs and their target genes across the three developmental stages of G. lucidum. A total of 168 unique milRNAs were predicted using a small RNA sequencing method. For them, 1612 target sequences corresponding to 1311 unique genes were predicted by degradome sequencing. We selected 42 predicted milRNAs and performed RT-PCR amplification and Sanger sequencing of the products. Five products were found to have sequences similar to those predicted, confirming the presence of milRNAs in G. lucidum, and demonstrating the difficulty in their validation. Among the 168 milRNAs, 111 were found to be significantly differentially expressed across the three developmental stages (q ≤ 0.05). The expression levels of 12 milRNAs were measured by stem-loop quantitative real-time polymerase chain reaction. Eight of them were in line with the sequencing results (r ≥ 0.9, p ≤ 0.05). These 12 milRNAs and their target genes form 16 milRNA-target gene pairs. The expression profiles of 8 of these 16 miRNA-target pairs were negatively correlated, according to real-time quantitative analysis, whereas the other eight pairs were positively correlated. Furthermore, the results of functional enrichment analysis showed that the target genes of milRNAs mapped to the Gene Ontology terms ‘GTP binding’ and ‘FAD binding’ were enriched in specific developmental stages. These target genes were related to the biosynthesis of triterpenes and polysaccharides and lignin degradation pathway in G. lucidum. In summary, this study indicates that milRNAs may play crucial regulatory roles in various biological processes of G. lucidum and open up new avenues for research on milRNAs’ biosyntheses and functions in basidiomycetes.