MicroRNA-like RNA Functions Are Required for the Biosynthesis of Active Compounds in the Medicinal Fungus Sanghuangporus vaninii

ABSTRACT miRNA-like RNAs (milRNAs) have been recognized as sequence-specific regulators of posttranscriptional regulation of gene expression in eukaryotes. However, the functions of hundreds of fungal milRNAs in the biosynthesis of metabolic components are obscure. Sanghuangporus produces diverse bioactive compounds and is widely used in Asian countries. Here, genes encoding two Dicers, four Argonautes, and four RdRPs were identified and characterized in Sanghuangporus vanini. Due to the lack of an efficient gene manipulation system, the efficacy of spray-induced gene silencing (SIGS) was determined in S. vanini, which showed efficient double-stranded RNA (dsRNA) uptake and gene silencing efficiency. SIGS-mediated gene knockdown showed that SVRDRP-3, SVRDRP-4, SVDICER-1, and SVDICER-2 were critical for mycelial biomass, flavonoid, triterpenoid, and polysaccharide production. Illumina deep sequencing was performed to characterize the milRNAs from S. vanini mycelium and fruiting body. A total of 31 milRNAs were identified, out of which, SvmilR10, SvmilR17, and SvmilR33 were Svrdrp-4- and Svdicer-1-dependent milRNAs. Importantly, SIGS-mediated overexpression of SvmilR10 and SvmilR33 resulted in significant changes in the yields of flavonoids, triterpenoids, and polysaccharides. Further analysis showed that these milRNA target genes encoding the retrotransposon-derived protein PEG1 and histone-lysine N-methyltransferase were potentially downregulated in the milRNA overexpressing strain. Our results revealed that S. vanini has high external dsRNA and small RNA uptake efficiency and that milRNAs may play crucial regulatory roles in the biosynthesis of bioactive compounds. IMPORTANCE Fungi can take up environmental RNA that can silence fungal genes with RNA interference, which prompts the development of SIGS. Efficient dsRNA and milRNA uptake in S. vanini, successful dsRNA-targeted gene block, and the increase in intracellular miRNA abundance showed that SIGS technology is an effective and powerful tool for the functional dissection of fungal genes and millRNAs. We found that the RdRP, Dicer, and Argonaute genes are critical for mycelial biomass and bioactive compound production. Our study also demonstrated that overexpressed SVRDRP-4- and SVDICER-1-dependent milRNAs (SvmilR10 and SvmilR33) led to significant changes in the yields of the three active compounds. This study not only provides the first report on SIGS-based gene and milRNA function exploration, but also provides a theoretical platform for exploration of the functions of milRNAs involved in biosynthesis of metabolic compounds in fungi.

domain, a Dicer_dimer domain, and two RNase III domains; Dicer-2 also includes a PAZ domain (Fig. 1B). AGO-1, AGO-3, and AGO-4 contain four conserved domains, an ArgoN domain, an ArgoL1 domain, a PAZ domain, and a Piwi domain, whereas no ArgoN domain was searched in AGO-2 (Fig. 1C).
Phylogenetic analysis showed that the number of RdRPs varied among different fungi; S. vaninii and S. baumii possessed four RdRPs, while three RdRPs were detected in N. crassa and A. flavus (Fig. 1D). RdRP-3 and RdRP-4 in Sanghuangporus species formed supported monophyletic clades that are sister to each other, suggesting distinct roles among these proteins in different fungi. The number of Dicers varied from two (S. vaninii and S. baumii) to four (A. flavus) (Fig. 1E). In the Dicer-1 group, S. vaninii and S. baumii formed separate clades and shared 56% sequence similarity with N. crassa and A. flavus. The phylogenetic tree constructed on the basis of the AGO protein sequences showed that AGO-2, AGO-3, and AGO-4 from S. vaninii and S. baumii were resolved into a distinct group (Fig. 1F).
To further explore the potential roles of milRNAs in S. vaninii, the expression profiles of four Svrdrps, two Svdicers, and four Svagos in solid PDA (20 days), liquid PDS (10 days), and fruiting bodies (3 years old) were assessed with qPCR. The results showed that all 10 genes were actively transcribed at all investigated developmental stages (Fig. 2). In particular, the expression trends of all four Svrdrps were highly consistent in each stage and showed decreased expression in solid PDA samples compared with that in both liquid PDS and fruiting body samples ( Fig. 2A). The highest gene expression levels of the two Svdicers were observed in the fruiting body, and Svdicer-2 expression decreased from the liquid PDS sample to the solid PDA sample, Silencing of AGO, Dicer, and RdRP genes with dsRNA uptake. To determine whether dsRNA was taken up by S. vaninii, the efficiency of dsRNA uptake was tested with fluorescein-labeled Svago-1 dsRNA (Svago-1-dsRNA). Mycelium from both solid PDA and liquid PDS medium were treated with fluorescein-labeled Svago-1-dsRNA to determine the fungal dsRNA uptake efficiency at different stages. Fluorescence signals were observed after the removal of labeled dsRNA from the fungal cells, and the results showed that mycelia from both solid PDA and liquid PDS medium showed very efficient dsRNA uptake (Fig. 3A). We also found that the expression levels of Svago-1 and Svago-1 were significantly downregulated after treatment with Svago-1 dsRNA (Fig. 3B). Remarkably, dsRNA-targeted silencing was more efficiently taken up in liquid PDS than in solid PDA; thus, the fungus cultured in liquid PDS medium was used for gene silencing with dsRNA uptake in the following study. Moreover, the silencing efficiency of Svago-1 was confirmed by using semiquantitative RT-PCR (Fig. 3C).
Combining the above results, silencing of the AGO, Dicer and RdRP genes was performed with dsRNA uptake. qPCR results showed that all genes were found to be significantly downregulated (2-to 5-fold decrease in transcript levels) in the dsRNA-treated fungus (Fig. 3D). Thus, the application of dsRNA uptake to silence the targeted genes could be an effective strategy to explore gene function in S. vaninii.
Mycelial biomass and bioactive compound production from different genesilenced fungi. To further explore the possible roles of the AGO, Dicer and RdRP genes in the production of bioactive compounds, yields of flavonoids, triterpenoids, and polysaccharides were detected. First, we determined the fungal biomass by measuring the mycelial dry weight. Compared with NC (gfp-dsRNA treatment), Svrdrp-3-dsRNA and Svrdrp-4-dsRNA showed considerable decreases, whereas Svdicer-1-dsRNA, Svdicer-2-dsRNA, and Svago-1-dsRNA revealed significant increases in mycelial biomass (Fig. 4A). Flavonoid production of Svdicer-1-dsRNA, Svdicer-2-dsRNA, and Svago-1-dsRNA was significantly higher than that in NC (Fig. 4B). However, compared with NC, a reduction was observed in flavonoid production in Svrdrp-1-, Svrdrp-3-, Svrdrp-4-, or Svago-4-silenced mycelium. For triterpenoid production, Svrdrp-1-dsRNA, Svrdrp-3-dsRNA, Svrdrp-4-dsRNA, Svdicer-1-dsRNA, and Svago-4-dsRNA were weaker but stronger in Svdicer-2-dsRNA and Svago-1-dsRNA than NC (Fig. 4C). Polysaccharide contents increased significantly by silencing Svdicer-1 but decreased significantly after Svrdrp-1, Svrdrp-3, Svrdrp-4, or Svago-4 was silenced in S. vaninii (Fig. 4D). Overall, Svrdrp-3, Svrdrp-4, Svdicer-1, and Svdicer-2 were critical for mycelial biomass and the production of the three bioactive compounds. Analysis of milRNA characteristics in S. vaninii. To investigate the milRNA characteristics, we analyzed the milRNAs generated in S. vaninii from liquid PDS medium-cultured fungus (MY) and fruiting body (FB). A total of 31 milRNAs were obtained, and six milRNAs were expressed exclusively in MY, whereas twenty-three were expressed in both MY and FB ( Fig. 5A and Table S1). The first base at the 59 terminus of milRNAs 21 to 23 bp in length had a strong preference for "C," in contrast to milRNAs with lengths of 24 to 29 bp (Fig. 5B). All milRNA abundance was normalized with TPM normalization, and an absolute value of log 2 ratio $ 1 (P , 0.05) was used as the threshold to determine significant differences in milRNA expression. A total of 15 milRNAs were found to have significantly different expression between FB and MY: compared with their expression in MY, four milRNAs were upregulated while the remainder were downregulated in FB (Fig. 5C). To validate the high-sequencing data, we analyzed the expression profiles of all differentially expressed milRNAs with qPCR. The expression profiles of all 15 milRNAs were consistent with the sequencing data and showed significantly differential expression between FB and MY (Fig. 5D).
To understand the biological functions triggered by milRNAs, the target genes of these milRNAs were identified. A total of 213 unique coding genes were predicted to be targets of these milRNAs, and a total of 212 milRNA-target pairs were identified ( Fig. 5E and Table S2). Moreover, analysis of interaction networks suggested that most Fungal miRNA Roles on Metabolic Component Biosynthesis Microbiology Spectrum of these milRNAs regulated several genes, and some milRNAs, such as SvmilR9, SvmilR16, SvmilR24, and SvmilR46, had connections with more than 20 genes. RdRP-and Dicer-dependent biogenesis of milRNAs involved in mycelial biomass and bioactive compound production. As in our study above, Svrdrp-3, Svrdrp-4, Svdicer-1, and Svdicer-2 were critical for mycelial biomass and the production of the three bioactive compounds. To identify these gene-dependent milRNAs, the expression characteristics of five milRNAs with high abundance (more than 100 counts in both FB and MY) were examined from these gene-silenced fungi. Silencing Svrdrp-4 or Svdicer-1 resulted in a reduction in the expression of SvmilR10, SvmilR17, and SvmilR33, whereas SvmilR33 showed increased expression of Svrdrp-3-dsRNA compared with that in NC (Fig. 6A). Thus, SvmilR10, SvmilR17 and SvmilR33 were regarded as Svrdrp-4and Svdicer-1-dependent milRNAs.
To investigate whether synthetic milRNA was taken up by S. vaninii, Cy3-labeled miRNA was added to the liquid PDS culture of mycelium. Fluorescence signals were clearly observed at 48 h (Fig. 6B). Most importantly, the milRNA expression levels of SvmilR10, SvmilR17 and SvmilR33 were considerably increased (more than 15-fold increase in transcript levels) in S. vaninii treated with different agomir miRNAs (Fig. 6C). Thus, S. vaninii took up external miRNAs efficiently.
Next, we investigated the mycelial biomass and bioactive compound production from these agomir miRNA-treated fungi. Compared with the control (water treatment) and negative-control (NC) treatment, agomir-17 showed a significant reduction in fungal biomass. The yields of both flavonoids and polysaccharides were increased, while the production of triterpenoids was decreased in the fungi treated with agomir-10 ( Fig. 6D). Overexpressed SvmilR33 in the fungus showed the opposite change for the Fungal miRNA Roles on Metabolic Component Biosynthesis Microbiology Spectrum triterpenoid and polysaccharide contents, in which polysaccharide yield was decreased while triterpenoid yield was increased. Subsequently, the expression levels of predicted targets of SvmilR10, SvmilR17 and SvmilR33 were studied in agomir miRNA-treated fungi. The results showed that the gene encoding the retrotransposon-derived protein PEG10 was significantly decreased in agomir-10-treated fungus, and the expression level of the histone-lysine N-methyltransferase gene was downregulated in the fungus after agomir-10 treatment (Fig. 6E). These results indicate a role for Svrdrp-4-and Svdicer-1-dependent milRNAs in the regulation of mycelial biomass and bioactive compound production by gene silencing mediation.

DISCUSSION
In fungi, the RdRP family generates dsRNA from certain single-stranded RNAs or from the target mRNA, activating or amplifying the silencing response, respectively. Different RDRP proteins are involved in diverse functional mechanisms; for example, RdRP1 is  (27). In this study, four RdRPs were identified from S. vaninii and S. baumii, which was more than most filamentous fungi. Fungi exhibit diversity in the number of RdRPs, suggesting that milRNA-produced mechanisms are distinct among different fungi (28,29). In addition, RdRP-3, RdRP-4, AGO-2, AGO-3, and AGO-4 from S. vaninii and S. baumii shared low similarity with the other fungal RdRPs and AGOs, suggesting that Sanghuangporus species may possess specific RNAi mechanisms.
Recently, impressive numbers of studies on SIGS indicated that RNAi is an effective control strategy against fungal pathogens and confirmed that fungi could take up environmental RNA (dsRNA and sRNA), which can then silence fungal genes through environmental RNA interference (16,30). Undoubtedly, SIGS provides new strategies to study the gene function of some fungi that lack efficient gene manipulation systems. A recent study showed that Sanghuangporus species were binuclear fungi, which was limited for functional gene analysis with homologous recombination by DNA-mediated transformation and CRISPR/Cas9-based gene editing (9). Thus, RNA uptake efficiency was examined in S. vaninii because SIGS was dependent on the efficiency of fungal RNA uptake (31). We observed efficient dsRNA uptake in S. vaninii from both solid PDA and liquid PDS medium, but the silencing efficiency of fungus from PDA was less efficient than that from PDS, suggesting that differences in culture environment and cell types could affect fungal RNA uptake efficiency, which was consistent with an earlier report (16). In Fusarium graminearum, exogenous application of dsRNA triggering SIGS of CYP51 was more efficient than host-induced gene silencing, suggesting that fungi seem to respond more efficiently to dsRNA than to siRNA (32). In this study, although milRNA abundances increased more than 15-fold, only weak characteristic changes occurred in the fungus, indicating low gene silencing efficiency by sRNA uptake. Therefore, exogenous application of dsRNA is a powerful genetic tool used for gene silencing. In fact, other gene editing technologies should be further explored in the fungus because of the time-consuming, high-cost, and low-efficiency stability of SIGSbased functional gene analysis.
In fungi, RdRP produces dsRNA molecules from single-stranded RNA precursors, and then Dicer cleaves the dsRNA and generates sRNAs (33). Fungal biomass and all three bioactive compound yields changed significantly in the S. vaninii treated with Svrdrp-3-, Svrdrp-4-, Svdicer-1-, or Svdicer-2-targeted dsRNAs, suggesting that Svrdrp-3-, Svrdrp-4-, Svdicer-1-, or Svdicer-2-regulated milRNAs were involved in fungal development and secondary metabolite production. The AGO protein forms a complex (RISC, RNA-induced silencing complex) to block gene transcription by targeting complementary specific miRNA sequences as a guide (34). MiRNAs were loaded into AGO4 in the plant cytoplasm and then mediated DNA methylation in the nucleus to target gene translational repression (35). In this study, silencing Svago-4 resulted in a decrease in the yield of all three bioactive compounds in S. vaninii, indicating that miRNAs regulated the production of these compounds within SVAGO-4 to form RISC. In G. lucidum, F. graminearum and C. militaris, 134, 143 and 38milRNAs were expressed in the sexual developmental stage, and only 92, 49, and 19 milRNAs were obtained from the asexual development stage (22,29,36). A total of 31 milRNAs were obtained in S. vaninii, twenty-nine milRNAs were detected in the asexual stage (MY), and 25 milRNAs were expressed in the sexual developmental stage (FB), suggesting that milRNAs were needed for asexual development rather than for sexual development in S. vaninii, which is contrary to these fungi. Although thirty-one milRNAs were obtained in S. vaninii, we have not been able to find any of the milRNAs to be conserved, which is consistent with G. lucidum, F. graminearum and C. militaris. Species-specific miRNAs are considered to be miRNAs that have evolved recently; thus, milRNAs may have evolved independently and recently in fungi (37).
Various studies have shown the important aspect of miRNAs in regulating secondary metabolite biosynthesis in different plant species (35). For example, miR828a and miR396b target anthocyanin regulatory protein-C1 and kaempferol 3-O-beta-d-galactosyltransferase genes in Rauvolfia serpentine (35,38). In Arabidopsis, miR5021 targets genes encoding the GcpE protein and terpene cyclase involved in terpenoid biosynthesis, miR858a targets genes encoding the R2R3-MYB transcription factor involved in regulating flavonoid biosynthesis, and miR156 regulates anthocyanin accumulation by targeting the SPL9 gene (39)(40)(41). Here, we found that SVRDRP-3, SVRDRP-4, SVDICER-1, and SVDICER-2 were critical for bioactive compound production, and SvmilR10, SvmilR17, and SvmilR33 were regarded as Svrdrp-4-and Svdicer-1-dependent milRNAs. In the SvmilR10-overexpressing strain, flavonoid and polysaccharide contents were increased but triterpenoid production was decreased. One SvmilR10 target, the retrotransposon-derived protein PEG10 gene, was critical for normal cell homeostasis and downregulated in the strain after Agomir-10 uptake, suggesting that SvmilR10 could control bioactive compound production with gene regulation (42,43). Furthermore, overexpression of SvmilR33 in the fungus resulted in decreased polysaccharide contents but increased triterpenoid yield, suggesting that different molecular mechanisms were used by milRNAs while producing different bioactive compounds.
In summary, core proteins of the miRNA-based silencing mechanism (RdRPs, Dicers and AGOs) and milRNAs were analyzed and characterized from the medicinal fungus S. vaninii. Importantly, milRNA roles in fungal biomass and bioactive compound production were analyzed by gene silencing with external dsRNA and sRNA uptake. These results suggested that milRNAs may be involved in the biosynthesis of bioactive compounds in S. vaninii.

MATERIALS AND METHODS
Fungal strain and culture conditions. The fungus used in this study was S. vaninii Kangneng, which has been used for large-scale industrial production in Jiangsu KONEN Biological Engineering Co. Ltd. (Anhui, China) (8,26). PDA plates (20% potato, 2% dextrose, and 1.5% agar, wt/vol) were used for the fungal culture at 28°C in the dark. A total of 250 mL flasks containing 100 mL of liquid PDS (20% potato, 2% mulberry tree sawdust) medium inoculated with the mycelia were used for the collection of mycelia from liquid medium. For fruiting body production, 10-day-cultured mycelia from liquid medium were inoculated in short log sections as described in our previous study (8).
Phylogenetic analyses of RdRPs, Dicers, and AGOs. RdRPs, Dicers, and AGOs from S. vaninii (SV), Sanghuangporus baumii (SB); Neurospora crassa (NC); Aspergillus flavus (AF); and Arabidopsis thaliana (AT) were used for phylogenetic analysis to further determine these conserved RNAi proteins (29). The BLASTp algorithm, underpinned by the Pfam and CDD databases, was used for searches of conserved protein domains or motifs. The amino acid sequences of RdRP, Dicer, and AGO proteins from different organisms were obtained from GenBank according to our previous study (29). Alignment sequences were used for phylogenetic analysis by MEGA7 (http://www.megasoftware.net) using a neighbor-joining method.
RNA isolation, quantitative PCR (qPCR) and stem-loop real-time PCR. According to our previous study, the production of bioactive compounds from S. vaninii was different at different stages; in this study, milRNA functions in the biosynthesis of active compounds were explored. Thus, fungal samples grown at the stage with high bioactive compound production were collected with different cultivation methods (solid PDA, liquid PDS and fruiting bodies) (26). Different fungal samples collected from solid PDA (20 days), liquid PDS (10 days), and fruiting bodies (3 years old) were used for RNA extraction with TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's instructions. RNase-free DNase I was used for DNA genome removal, and the RNA concentration was evaluated. Then, cDNA was synthesized with a gDNA remover kit (Toyobo, Osaka, Japan) according to the manufacturer's instructions. qPCR amplification was performed with a SYBR green kit (TaKaRa) and Bio-Rad CFX96TM system (Hercules, CA, United States). Transcripts of different protein-coding genes were normalized to the control gene. The stem-loop real-time PCR method was performed to quantitate milRNA expression. In brief, a stem-loop RT primer was used to reverse-transcribe mature milRNAs to complementary DNAs (cDNAs). The 20 mL reverse-transcription reactions contained 500 ng of total RNA, 50 nM each individual stem-loop RT primer, 0.5 mM dNTP (TaKaRa, Dalian, China), 10 U M-MLV reverse transcriptase (TaKaRa, Dalian, China), and 4 U RNase inhibitor. qPCR amplification was performed as mentioned above. Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) was used as the control. The relative gene expression levels were calculated with the 2 2DDCt method. All experiments were performed with six replicates. Details of the primers used in this study are given in Table S3.
In vitro synthesis of dsRNA and agomir milRNA. To silence different core components of the RNAi machinery in S. vaninii, dsRNAs targeting RNA-dependent RNA polymerases, dicers, and argonautes were synthesized with T7 RNA polymerase (T7 RiboMAX Express RNAi system kit, Promega, Madison, WI, USA) based on established protocols (44). Briefly, ;300 bp DNA fragments, introduced at both the 59 and 39 ends with the T7 promoter sequence, were produced via PCR (PCR) from cDNA. The DNA fragments were used for in vitro transcription after confirmation by DNA sequencing, and then the final dsRNA was dissolved in nuclease-free water. The green fluorescent protein (GFP) gene was used as the negative control. Fluorescein-labeled dsRNA was synthesized with a Fluorescein RNA Labeling Mix kit following the manufacturer's instructions (Sigma, St. Louis, MO, USA) as previously described (45). Agomir miRNA, chemically modified and cholesterylated, was a stable miRNA mimic and synthesized by Shanghai GenePharma (GenePharma, Shanghai, China).
Fungal uptake of dsRNA or agomir miRNA in vitro. Conidia from the fungus grown on PDA medium were counted and diluted to 5 Â 10 6 spores/mL. Fifty mL of conidial suspension was transferred into 20 mL of liquid PDS or solid PDA medium cultured at 28°C in the dark for 10 days. Then, 100 mL of dsRNA (100 ng/mL) or agomir miRNA (50 ng/mL) was added to the liquid PDS medium or evenly spread over the surface of the solid PDA medium by brushing lightly with a soft brush. Synthesized dsRNA targeting gfp (gfp-dsRNA) was used as the negative control (NC). Before confocal microscopy and qPCR examination, dsRNA or agomir miRNA outside the fungal cells was removed by washing four times with sterilized water. Fungal mycelium was added to external gfp-dsRNA or agomir miRNA (micrON agomir NC number 24), and the mixtures were immediately used for dsRNA or agomir miRNA removal. Then, total RNA was extracted for qPCR to ensure clean external RNA removal.
For confocal microscopy examination, fungal mycelium was collected from PDS or PDA medium after 48 h of culture with the addition of fluorescent dsRNA or Cy3-labeled agomir miRNA. The fluorescent signal was analyzed using a confocal microscope (Zeiss LSM 880, Oberkochen, Germany).
Determination of active compounds. The yields of the active ingredients (polysaccharides, flavonoids, and triterpenoid) among fungi treated with dsRNA or agomir miRNA were determined as described in our previous report (26). In brief, samples were dried and ground into powder, and 50 mg of powder was used for the extraction and determination of polysaccharide, flavonoid, and triterpenoid contents by measuring the absorbance at 620 nm, 384 nm, and 551 nm, respectively.
sRNA library construction and sequencing. The mycelium (MY: harvested after 10 days of growth in the PDS medium) and fruiting bodies (FB: 3 years old), currently widely used methods for the Sanghuangporus collection, of S. vaninii were used to construct small RNA and cDNA libraries. For fungal sRNA library construction, fresh fruiting bodies were pulverized to powder as previously described (26). Samples were ground into a powder using liquid nitrogen, and RNAiso Plus reagent (TaKaRa, Shiga, Japan) was used for total RNA extraction. A NanoDrop ND-2000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, United States) and an Agilent 2100 Bioanalyzer (Agilent, United States) were used for RNA concentration evaluation after genomic DNA removal with RNase-free DNase I. To maximize target coverage, equal amounts of total RNA from each of the three replicates of samples were pooled for library construction. RNA (15 to 30 nt) was collected from total RNA on a 15% denaturing polyacrylamide gel, and then 59 and 39 of these RNAs were ligated to specific adaptors. cDNA libraries were sequenced with the Illumina HiSeq 2000 platform (BGI, Shenzhen, China) after reverse transcription.
Data analysis of small RNAs was performed as described in a previous study (29). Briefly, clean reads were obtained from raw reads by filtering out poor-quality reads after adaptor sequence removal. Adaptortrimmed unique sequences were aligned to the S. vaninii genome (https://www.ncbi.nlm.nih.gov/nuccore/ VVIT00000000.1/). Known noncoding RNAs from mapped reads, such as tRNA, snoRNA, and rRNA, were removed, and the remaining reads were aligned with mature miRNA sequences in the miRbase database to identify known milRNAs. According to the characteristics of the precursor of miRNA (a typical hairpin structure), the mideep 2 (https://www.mdc-berlin.de/content/mirdeep2-documentation) method was used to predict new miRNAs, and the secondary structure was successfully predicted. The abundance of milRNAs was Fungal miRNA Roles on Metabolic Component Biosynthesis Microbiology Spectrum normalized according to transcripts per million (TPM) normalization, and milRNA prediction was performed with PsRobot (https://tools4mirs.org/software/target_prediction/psrobot/) and TargetFinder (http://plantgrn .noble.org/psRNATarget/) according to plantlike target interactions and their targets with default parameters (35,46). Statistical analysis. Data were analyzed by one-way analysis of variance (ANOVA) and Tukey's HSD test. SPSS v23.0 software for Windows (SPSS Inc., Chicago, IL, USA) was used for statistical analysis. All results are expressed as the mean 6 standard error of the mean (SD). P , 0.05 was considered significant.
Data availability. Data are available in a public, open access repository. All sequence data from this study has been submitted to Sequence Read Archive (https://www.ncbi.nlm.nih.gov/sra) and can be accessed with accessions PRJNA798153, SRR17649062, and SRR17649063.

SUPPLEMENTAL MATERIAL
Supplemental material is available online only. SUPPLEMENTAL FILE 1, XLSX file, 0.02 MB.