Dicer-Like Proteins Regulate Sexual Development via the Biogenesis of Perithecium-Specific MicroRNAs in a Plant Pathogenic Fungus Fusarium graminearum

Trans-Kingdom RNA Dialogues: miRNA and milRNA Networks as Biotechnological Tools for Sustainable Crop Defense and Pathogen Control

MicroRNAs (miRNAs) are a class of non-coding RNAs approximately 20-24 nucleotides in length, which play a crucial role during gene regulation in plant-pathogen interaction. They negatively regulate the expression of target genes, primarily at the transcriptional or post-transcriptional level, through complementary base pairing with target gene sequences. Recent studies reveal that during pathogen infection, miRNAs produced by plants and miRNA-like RNAs (milRNAs) produced by fungi can regulate the expression of endogenous genes in their respective organisms and undergo trans-kingdom transfer. They can thereby negatively regulate the expression of target genes in recipient cells. These findings provide novel perspectives for deepening our understanding of the regulatory mechanisms underlying plant-pathogen interactions. Here, we summarize and discuss the roles of miRNAs and milRNAs in mediating plant-pathogen interactions via multiple pathways, providing new insights into the functions of these RNAs and their modes of action. Collectively, these insights lay a theoretical foundation for the targeted management of crop diseases.

Global discovery, expression pattern, and regulatory role of miRNA-like RNAs in Ascosphaera apis infecting the Asian honeybee larvae

Ascosphaera apis, a specialized fungal pathogen, causes lethal infection in honeybee larvae. miRNA-like small RNAs (milRNAs) are fungal small non-coding RNAs similar to miRNAs, which have been shown to regulate fungal hyphal growth, spore formation, and pathogenesis. Based on the transcriptome data, differentially expressed miRNA-like RNAs (DEmilRNAs) in A. apis infecting the Apis cerana cerana worker 4-, 5-, and 6-day-old larvae (Aa-4, Aa-5, and Aa-6) were screened and subjected to trend analysis, followed by target prediction and annotation as well as investigation of regulatory networks, with a focus on sub-networks relative to MAPK signaling pathway, glycerolipid metabolism, superoxide dismutase, and enzymes related to chitin synthesis and degradation. A total of 606 milRNAs, with a length distribution ranging from 18 nt to 25 nt, were identified. The first nucleotide of these milRNAs presented a bias toward U, and the bias patterns across bases of milRNAs were similar in the aforementioned three groups. There were 253 milRNAs, of which 68 up-and 54 down-regulated milRNAs shared by these groups. Additionally, the expression and sequences of three milRNAs were validated by stem-loop RT-PCR and Sanger sequencing. Trend analysis indicated that 79 DEmilRNAs were classified into three significant profiles (Profile4, Profile6, and Profile7). Target mRNAs of DEmilRNAs in these three significant profiles were engaged in 42 GO terms such as localization, antioxidant activity, and nucleoid. These targets were also involved in 120 KEGG pathways including lysine biosynthesis, pyruvate metabolism, and biosynthesis of antibiotics. Further investigation suggested that DEmilRNA-targeted mRNAs were associated with the MAPK signaling pathway, glycerolipid metabolism, superoxide dismutase, and enzymes related to chitin synthesis and degradation. Moreover, the binding relationships between aap-milR10516-x and ChsD as well as between aap-milR-2478-y and mkh1 were confirmed utilizing a combination of dual-luciferase reporter gene assay and RT-qPCR. Our data not only provide new insights into the A. apis proliferation and invasion, but also lay a basis for illustrating the DEmilRNA-modulated mechanisms underlying the A. apis infection.

An Antisense Long Non-Coding RNA, LncRsn, Is Involved in Sexual Reproduction and Full Virulence in Fusarium graminearum

Fusarium head blight (FHB), primarily caused by Fusarium graminearum, is a devastating crop disease that leads to significant declines in wheat yield and quality worldwide. Long non-coding RNAs (lncRNAs) are found to play significant functions in various biological processes, but their regulatory functions in the sexual reproduction and pathogenicity of F. graminearum have not been studied extensively. This study identified an antisense lncRNA, named lncRsn, located in the transcription initiation site region between the 5'-flanking gene FgSna and the 3'-flanking gene FgPta. A deletion mutant of lncRsn (ΔlncRsn) was constructed through homologous recombination. ΔlncRsn exhibited huge reductions in pathogen and sexual reproduction. Additionally, the deletion of lncRsn disrupted the biosynthesis of deoxynivalenol (DON) and impaired the formation of infection structures. RT-qPCR analysis reveals that lncRsn may negatively regulate the transcription of the target gene FgSna. This study found that lncRsn plays an important role in sexual and asexual reproduction, pathogenicity, virulence, osmotic stress, and cell wall integrity (CWI) in F. graminearum. Further characterization of pathogenesis-related genes and the reaction between lncRsn and protein-coding genes will aid in developing novel approaches for controlling F. graminearum diseases.

Harnessing RNA interference for the control of Fusarium species: A critical review

Fusarium fungi are a pervasive threat to global agricultural productivity. They cause a spectrum of plant diseases that result in significant yield losses and threaten food safety by producing mycotoxins that are harmful to human and animal health. In recent years, the exploitation of the RNA interference (RNAi) mechanism has emerged as a promising avenue for the control of Fusarium-induced diseases, providing both a mechanistic understanding of Fusarium gene function and a potential strategy for environmentally sustainable disease management. However, despite significant progress in elucidating the presence and function of the RNAi pathway in different Fusarium species, a comprehensive understanding of its individual protein components and underlying silencing mechanisms remains elusive. Accordingly, while a considerable number of RNAi-based approaches to Fusarium control have been developed and many reports of RNAi applications in Fusarium control under laboratory conditions have been published, the applicability of this knowledge in agronomic settings remains an open question, and few convincing data on RNAi-based disease control under field conditions have been published. This review aims to consolidate the current knowledge on the role of RNAi in Fusarium disease control by evaluating current research and highlighting important avenues for future investigation.

RNAi-Based Approaches to Control Mycotoxin Producers: Challenges and Perspectives

Mycotoxin contamination of food and feed is a worldwide problem that needs to be addressed with highly efficient and biologically safe techniques. RNA interference (RNAi) is a natural mechanism playing an important role in different processes in eukaryotes, including the regulation of gene expression, maintenance of genome stability, protection against viruses and others. Recently, RNAi-based techniques have been widely applied for the purposes of food safety and management of plant diseases, including those caused by mycotoxin-producing fungi. In this review, we summarize the current state-of-the-art RNAi-based approaches for reducing the aggressiveness of key toxigenic fungal pathogens and mycotoxin contamination of grain and its products. The ways of improving RNAi efficiency for plant protection and future perspectives of this technique, including progress in methods of double-stranded RNA production and its delivery to the target cells, are also discussed.

Challenges and Opportunities Arising from Host-Botrytis cinerea Interactions to Outline Novel and Sustainable Control Strategies: The Key Role of RNA Interference

The necrotrophic plant pathogenic fungus Botrytis cinerea (Pers., 1794), the causative agent of gray mold disease, causes significant losses in agricultural production. Control of this fungal pathogen is quite difficult due to its wide host range and environmental persistence. Currently, the management of the disease is still mainly based on chemicals, which can have harmful effects not only on the environment and on human health but also because they favor the development of strains resistant to fungicides. The flexibility and plasticity of B. cinerea in challenging plant defense mechanisms and its ability to evolve strategies to escape chemicals require the development of new control strategies for successful disease management. In this review, some aspects of the host-pathogen interactions from which novel and sustainable control strategies could be developed (e.g., signaling pathways, molecules involved in plant immune mechanisms, hormones, post-transcriptional gene silencing) were analyzed. New biotechnological tools based on the use of RNA interference (RNAi) are emerging in the crop protection scenario as versatile, sustainable, effective, and environmentally friendly alternatives to the use of chemicals. RNAi-based fungicides are expected to be approved soon, although they will face several challenges before reaching the market.

RNA-dependent RNA polymerases regulate ascospore discharge through the exonic-sRNA-mediated RNAi pathway

Ascospores, forcibly released into the air from perithecia, are the primary inoculum for Fusarium head blight. In Fusarium graminearum, the biological functions of four RNA-dependent RNA polymerases (RdRPs) (Fgrdrp1-4) have been reported, but their regulatory mechanisms are poorly understood and the function of Fgrdrp5 is still unknown. In this study, we found that in addition to Fgrdrp1 and Fgrdrp2, Fgrdrp5 also plays an important role in ascospore discharge, and they all participate in the generation of turgor pressure in a polyol-dependent manner. Moreover, these three genes all affect the maturation of ascospores. Deep sequencing and co-analysis of small RNA and mRNA certified that Fgrdrp1, Fgrdrp2, and Fgrdrp5 partly share their functions in the biogenesis and accumulation of exonic small interference RNA (ex-siRNA), and these three RdRPs negatively regulate the expression levels of ex-siRNA corresponding genes, including certain genes associated with ascospore development or discharge. Furthermore, the differentially expressed genes of deletion mutants, those involved in lipid and sugar metabolism or transport as well as sexual development-related transcription factors, may also contribute to the defects in ascospore maturation or ascospore discharge. In conclusion, our study suggested that the components of the dicer-dependent ex-siRNA-mediated RNA interference pathway include at least Fgrdrp1, Fgrdrp2, and Fgrdrp5.

IMPORTANCE

We found that in addition to Fgrdrp1 and Fgrdrp2, Fgrdrp5 also plays important roles in ascospore maturation and ascospore discharge of Fusarium graminearum. These three RNA-dependent RNA polymerases participate in the biogenesis and accumulation of exonic small interference RNA and then regulate ascospore discharge.

Advances in Understanding Fusarium graminearum: Genes Involved in the Regulation of Sexual Development, Pathogenesis, and Deoxynivalenol Biosynthesis

The wheat head blight disease caused by Fusarium graminearum is a major concern for food security and the health of both humans and animals. As a pathogenic microorganism, F. graminearum produces virulence factors during infection to increase pathogenicity, including various macromolecular and small molecular compounds. Among these virulence factors, secreted proteins and deoxynivalenol (DON) are important weapons for the expansion and colonization of F. graminearum. Besides the presence of virulence factors, sexual reproduction is also crucial for the infection process of F. graminearum and is indispensable for the emergence and spread of wheat head blight. Over the last ten years, there have been notable breakthroughs in researching the virulence factors and sexual reproduction of F. graminearum. This review aims to analyze the research progress of sexual reproduction, secreted proteins, and DON of F. graminearum, emphasizing the regulation of sexual reproduction and DON synthesis. We also discuss the application of new gene engineering technologies in the prevention and control of wheat head blight.

Genome-wide identification of Sclerotinia sclerotiorum small RNAs and their endogenous targets

BACKGROUND

Several phytopathogens produce small non-coding RNAs of approximately 18-30 nucleotides (nt) which post-transcriptionally regulate gene expression. Commonly called small RNAs (sRNAs), these small molecules were also reported to be present in the necrotrophic pathogen Sclerotinia sclerotiorum. S. sclerotiorum causes diseases in more than 400 plant species, including the important oilseed crop Brassica napus. sRNAs can further be classified as microRNAs (miRNAs) and short interfering RNAs (siRNAs). Certain miRNAs can activate loci that produce further sRNAs; these secondary sRNA-producing loci are called 'phased siRNA' (PHAS) loci and have only been described in plants. To date, very few studies have characterized sRNAs and their endogenous targets in S. sclerotiorum.

RESULTS

We used Illumina sequencing to characterize sRNAs from fungal mycelial mats of S. sclerotiorum spread over B. napus leaves. In total, eight sRNA libraries were prepared from in vitro, 12 h post-inoculation (HPI), and 24 HPI mycelial mat samples. Cluster analysis identified 354 abundant sRNA clusters with reads of more than 100 Reads Per Million (RPM). Differential expression analysis revealed upregulation of 34 and 57 loci at 12 and 24 HPI, respectively, in comparison to in vitro samples. Among these, 25 loci were commonly upregulated. Altogether, 343 endogenous targets were identified from the major RNAs of 25 loci. Almost 88% of these targets were annotated as repeat element genes, while the remaining targets were non-repeat element genes. Fungal degradome reads confirmed cleavage of two transposable elements by one upregulated sRNA. Altogether, 24 milRNA loci were predicted with both mature and milRNA* (star) sequences; these are both criteria associated previously with experimentally verified miRNAs. Degradome sequencing data confirmed the cleavage of 14 targets. These targets were related to repeat element genes, phosphate acetyltransferases, RNA-binding factor, and exchange factor. A PHAS gene prediction tool identified 26 possible phased interfering loci with 147 phasiRNAs from the S. sclerotiorum genome, suggesting this pathogen might produce sRNAs that function similarly to miRNAs in higher eukaryotes.

CONCLUSIONS

Our results provide new insights into sRNA populations and add a new resource for the study of sRNAs in S. sclerotiorum.

milR20 negatively regulates the development of fruit bodies in Pleurotus cornucopiae

The mechanism underlying the development of fruit bodies in edible mushroom is a widely studied topic. In this study, the role of milRNAs in the development of fruit bodies of Pleurotus cornucopiae was studied by comparative analyses of the mRNAs and milRNAs at different stages of development. The genes that play a crucial role in the expression and function of milRNAs were identified and subsequently expressed and silenced at different stages of development. The total number of differentially expressed genes (DEGs) and differentially expressed milRNAs (DEMs) at different stages of development was determined to be 7,934 and 20, respectively. Comparison of the DEGs and DEMs across the different development stages revealed that DEMs and its target DEGs involved in the mitogen-activated protein kinase (MAPK) signaling pathway, protein processing in endoplasmic reticulum, endocytosis, aminoacyl-tRNA biosynthesis, RNA transport, and other metabolism pathways, which may play important roles in the development of the fruit bodies of P. cornucopiae. The function of milR20, which targeted pheromone A receptor g8971 and was involved in the MAPK signaling pathway, was further verified by overexpression and silencing in P. cornucopiae. The results demonstrated that the overexpression of milR20 reduced the growth rate of mycelia and prolonged the development of the fruit bodies, while milR20 silencing had an opposite effect. These findings indicated that milR20 plays a negative role in the development of P. cornucopiae. This study provides novel insights into the molecular mechanism underlying the development of fruit bodies in P. cornucopiae.

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

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.

Trichoderma atroviride hyphal regeneration and conidiation depend on cell-signaling processes regulated by a microRNA-like RNA

The ability to respond to injury is essential for the survival of an organism and involves analogous mechanisms in animals and plants. Such mechanisms integrate coordinated genetic and metabolic reprogramming events requiring regulation by small RNAs for adequate healing of the wounded area. We have previously reported that the response to injury of the filamentous fungus Trichoderma atroviride involves molecular mechanisms closely resembling those of plants and animals that lead to the formation of new hyphae (regeneration) and the development of asexual reproduction structures (conidiophores). However, the involvement of microRNAs in this process has not been investigated in fungi. In this work, we explore the participation of microRNA-like RNAs (milRNAs) molecules by sequencing messenger and small RNAs during the injury response of the WT strain and RNAi mutants. We found that Dcr2 appears to play an important role in hyphal regeneration and is required to produce the majority of sRNAs in T. atroviride. We also determined that the three main milRNAs produced via Dcr2 are induced during the damage-triggered developmental process. Importantly, elimination of a single milRNA phenocopied the main defects observed in the dcr2 mutant. Our results demonstrate the essential role of milRNAs in hyphal regeneration and asexual development by post-transcriptionally regulating cellular signalling processes involving phosphorylation events. These observations allow us to conclude that fungi, like plants and animals, in response to damage activate fine-tuning regulatory mechanisms.

FoQDE2-dependent milRNA promotes Fusarium oxysporum f. sp. cubense virulence by silencing a glycosyl hydrolase coding gene expression

MicroRNAs (miRNAs) are small non-coding RNAs that regulate protein-coding gene expression primarily found in plants and animals. Fungi produce microRNA-like RNAs (milRNAs) that are structurally similar to miRNAs and functionally important in various biological processes. The fungus Fusarium oxysporum f. sp. cubense (Foc) is the causal agent of Banana Fusarium vascular wilt that threatens global banana production. It remains uncharacterized about the biosynthesis and functions of milRNAs in Foc. In this study, we investigated the biological function of milRNAs contributing to Foc pathogenesis. Within 24 hours post infecting the host, the Argonaute coding gene FoQDE2, and two Dicer coding genes FoDCL1 and FoDCL2, all of which are involved in milRNA biosynthesis, were significantly induced. FoQDE2 deletion mutant exhibited decreased virulence, suggesting the involvement of milRNA biosynthesis in the Foc pathogenesis. By small RNA sequencing, we identified 364 small RNA-producing loci in the Foc genome, 25 of which were significantly down-regulated in the FoQDE2 deletion mutant, from which milR-87 was verified as a FoQDE2-depedent milRNA based on qRT-PCR and Northern blot analysis. Compared to the wild-type, the deletion mutant of milR-87 was significantly reduced in virulence, while overexpression of milR-87 enhanced disease severity, confirming that milR-87 is crucial for Foc virulence in the infection process. We furthermore identified FOIG_15013 (a glycosyl hydrolase-coding gene) as the direct target of milR-87 based on the expression of FOIG_15013-GFP fusion protein. The FOIG_15013 deletion mutant displayed similar phenotypes as the overexpression of milR-87, with a dramatic increase in the growth, conidiation and virulence. Transient expression of FOIG_15013 in Nicotiana benthamiana leaves activates the host defense responses. Collectively, this study documents the involvement of milRNAs in the manifestation of the devastating fungal disease in banana, and demonstrates the importance of milRNAs in the pathogenesis and other biological processes. Further analyses of the biosynthesis and expression regulation of fungal milRNAs may offer a novel strategy to combat devastating fungal diseases.

The fungus Fusarium oxysporum f. sp. cubense (Foc) is the causal agent of Banana Fusarium vascular wilt that threatens global banana production. However, knowledge about pathogenesis of Foc is limited. In particular, pathogenic regulatory mechanism of the microRNA like small RNAs (milRNAs) found in Foc is unknown. Here, we found that FoQDE2, an Argonaute coding gene, and two Dicer coding genes FoDCL1 and FoDCL2, which are involved in milRNA biosynthesis, are significantly induced during the early infection stage of Foc. The results suggested that the milRNAs biosynthesis mediated by these genes may play an active role in the infection process of Foc. Based on this assumption, we subsequently found a FoQDE2-dependent milRNA (milR-87) and identified its target gene. Functional analysis showed that FoQDE2, milR-87 and its target gene were involved in the pathogenicity of Foc in different degree. The studies help us gain insight into the pathogenesis with FoQDE2, milR-87, and its target gene as central axis in Foc. The identified pathogenicity-involved milRNA provides an active target for developing novel and efficient biocontrol agents against Banana Fusarium wilt.

lncRsp1, a long noncoding RNA, influences Fgsp1 expression and sexual reproduction in Fusarium graminearum

Long noncoding RNAs (lncRNAs) are crucial regulators of gene expression in many biological processes, but their biological functions remain largely unknown, especially in fungi. Fusarium graminearum is an important pathogen that causes the destructive disease Fusarium head blight (FHB) or head scab disease on wheat and barley. In our previous RNA sequencing (RNA-Seq) study, we discovered that lncRsp1 is an lncRNA that is located +99 bp upstream of a putative sugar transporter gene, Fgsp1, with the same transcription direction. Functional studies revealed that ΔlncRsp1 and ΔFgsp1 were normal in growth and conidiation but had defects in ascospore discharge and virulence on wheat coleoptiles. Moreover, lncRsp1 and Fgsp1 were shown to negatively regulate the expression of several deoxynivalenol (DON) biosynthesis genes, TRI4, TRI5, TRI6, and TRI13, as well as DON production. Further analysis showed that the overexpression of lncRsp1 enhanced the ability of ascospore release and increased the mRNA expression level of the Fgsp1 gene, while lncRsp1-silenced strains reduced ascospore discharge and inhibited Fgsp1 expression during the sexual reproduction stage. In addition, the lncRsp1 complementary strains lncRsp1-LC-1 and lncRsp1-LC-2 restored ascospore discharge to the level of the wild-type strain PH-1. Taken together, our results reveal the distinct and specific functions of lncRsp1 and Fgsp1 in F. graminearum and principally demonstrate that lncRsp1 can affect the release of ascospores by regulating the expression of Fgsp1.

A Novel DCL2-Dependent Micro-Like RNA Vm-PC-3p-92107_6 Affects Pathogenicity by Regulating the Expression of Vm-VPS10 in Valsa mali

Dicer proteins are mainly responsible for generating small RNAs (sRNAs), which are involved in gene silencing in most eukaryotes. In previous research, two DCL proteins in Valsa mali, the pathogenic fungus causing apple tree Valsa canker, were found associated with both the pathogenicity and generation of sRNAs. In this study, the differential expression of small interfering RNAs (siRNAs) and miRNA-like RNAs (milRNAs) was analyzed based on the deep sequencing of the wild type and Vm-DCL2 mutant, respectively. Overall, the generation of 40 siRNAs and 18 milRNAs was evidently associated with Vm-DCL2. The target genes of milRNAs were then identified using degradome sequencing; according to the prediction results, most candidate targets are related to pathogenicity. Further, expression of Vm-PC-3p-92107_6 was confirmed in the wild type but not in the Vm-DCL2 mutant. Moreover, the pathogenicity of Vm-PC-3p-92107_6 deletion mutants (ΔVm-PC-3p-92107_6) and the over-expression transformants (Vm-PC-3p-92107_6-OE) was significantly increased and decreased, respectively. Based on those degradome results, vacuolar protein sorting 10 (Vm-VPS10) was identified as the target of Vm-PC-3p-92107_6. Co-expression analysis in tobacco leaves further confirmed that Vm-PC-3p-92107_6 could suppress the expression of Vm-VPS10. Meanwhile, the expression levels of Vm-PC-3p-92107_6 and Vm-VPS10 displayed divergent trends in ΔVm-PC-3p-92107_6 and Vm-PC-3p-92107_6-OE, respectively. Perhaps most importantly, ΔVm-VPS10 featured a significant reduction in pathogenicity. Taken together, our results indicate that a DCL2-dependent milRNA Vm-PC-3p-92107_6 plays roles in pathogenicity by regulating the expression of Vm-VPS10. This study lays a foundation for the comprehensive analysis of pathogenic mechanisms of V. mali and deepens our understanding of the generation and function of fungal sRNA.

Small RNAs in Plant Immunity and Virulence of Filamentous Pathogens

Gene silencing guided by small RNAs governs a broad range of cellular processes in eukaryotes. Small RNAs are important components of plant immunity because they contribute to pathogen-triggered transcription reprogramming and directly target pathogen RNAs. Recent research suggests that silencing of pathogen genes by plant small RNAs occurs not only during viral infection but also in nonviral pathogens through a process termed host-induced gene silencing, which involves trans-species small RNA trafficking. Similarly, small RNAs are also produced by eukaryotic pathogens and regulate virulence. This review summarizes the small RNA pathways in both plants and filamentous pathogens, including fungi and oomycetes, and discusses their role in host-pathogen interactions. We highlight secondarysmall interfering RNAs of plants as regulators of immune receptor gene expression and executors of host-induced gene silencing in invading pathogens. The current status and prospects of trans-species gene silencing at the host-pathogen interface are discussed.

RNAi-Related Dicer and Argonaute Proteins Play Critical Roles for Meiocyte Formation, Chromosome-Axes Lengths and Crossover Patterning in the Fungus Sordaria macrospora

RNA interference (RNAi) is a cellular process involving small RNAs that target and regulate complementary RNA transcripts. This phenomenon has well-characterized roles in regulating gene and transposon expression. In addition, Dicer and Argonaute proteins, which are key players of RNAi, also have functions unrelated to gene repression. We show here that in the filamentous Ascomycete Sordaria macrospora, genes encoding the two Dicer (Dcl1 and Dcl2) and the two Argonaute (Sms2 and Qde2) proteins are dispensable for vegetative growth. However, we identified roles for all four proteins in the sexual cycle. Dcl1 and Sms2 are essential for timely and successful ascus/meiocyte formation. During meiosis per se, Dcl1, Dcl2, and Qde2 modulate, more or less severely, chromosome axis length and crossover numbers, patterning and interference. Additionally, Sms2 is necessary both for correct synaptonemal complex formation and loading of the pro-crossover E3 ligase-protein Hei10. Moreover, meiocyte formation, and thus meiotic induction, is completely blocked in the dcl1 dcl2 and dcl1 sms2 null double mutants. These results indicate complex roles of the RNAi machinery during major steps of the meiotic process with newly uncovered roles for chromosomes-axis length modulation and crossover patterning regulation.

A novel antisense long non-coding RNA participates in asexual and sexual reproduction by regulating the expression of GzmetE in Fusarium graminearum

Fusarium graminearum is an important worldwide pathogen that causes Fusarium head blight in wheat, barley, maize and other grains. LncRNAs play important roles in many biological processes, but little is known about their functions and mechanisms in filamentous fungi. Here, we report that a natural antisense RNA, GzmetE-AS, is transcribed from the opposite strand of GzmetE. GzmetE encodes a homoserine O-acetyltransferase, which is important for sexual development and plant infection. The expression of GzmetE-AS was increased significantly during the conidiation stage, while GzmetE was upregulated in the late stage of sexual reproduction. Overexpression of GzmetE-AS inhibited the transcription of GzmetE. In contrast, the expression of GzmetE was significantly increased in GzmetE-AS transcription termination strain GzmetE-AS-T. Furthermore, GzmetE-AS-T produced more perithecia and facilitated the ascospore discharge, resembling the phenotype of GzmetE overexpressing strains. However, overexpression of GzmetE-AS in ∆dcl1/2 strain cannot inhibit the expression of GzmetE, and the GzmetE nat-siRNA is also significantly reduced in ∆dcl1/2 mutant. Taken together, we have identified a novel antisense lncRNA GzmetE-AS, which is involved in asexual and sexual reproduction by regulating its antisense gene GzmetE through RNAi pathway. Our findings reveal that the lncRNA plays critical roles in the development of F. graminearum.

The Evolutionary Significance of RNAi in the Fungal Kingdom

RNA interference (RNAi) was discovered at the end of last millennium, changing the way scientists understood regulation of gene expression. Within the following two decades, a variety of different RNAi mechanisms were found in eukaryotes, reflecting the evolutive diversity that RNAi entails. The essential silencing mechanism consists of an RNase III enzyme called Dicer that cleaves double-stranded RNA (dsRNA) generating small interfering RNAs (siRNAs), a hallmark of RNAi. These siRNAs are loaded into the RNA-induced silencing complex (RISC) triggering the cleavage of complementary messenger RNAs by the Argonaute protein, the main component of the complex. Consequently, the expression of target genes is silenced. This mechanism has been thoroughly studied in fungi due to their proximity to the animal phylum and the conservation of the RNAi mechanism from lower to higher eukaryotes. However, the role and even the presence of RNAi differ across the fungal kingdom, as it has evolved adapting to the particularities and needs of each species. Fungi have exploited RNAi to regulate a variety of cell activities as different as defense against exogenous and potentially harmful DNA, genome integrity, development, drug tolerance, or virulence. This pathway has offered versatility to fungi through evolution, favoring the enormous diversity this kingdom comprises.

Characterization of microRNA-like RNAs associated with sclerotial development in Sclerotinia sclerotiorum

Sclerotinia sclerotiorum is a model necrotrophic pathogen causing great economic losses worldwide. Sclerotia are dormant structures that play significant biological and ecological roles in the life and disease cycles of S. sclerotiorum and other species of sclerotia-forming fungi. microRNA-like RNAs (milRNAs) as non-coding small RNAs play regulatory roles in fungal development and pathogenicity. Therefore, milRNAs associated with sclerotial development in S. sclerotiorum were investigated in this study. A total of 275 milRNAs with induced expression during sclerotia development were identified, in which 51 were differentially expressed. The target genes of all milRNAs were predicted. The putative functions of the targets regulated by milRNAs were annotated by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. The expression levels of six selected milRNAs that coordinated with their corresponding targets were validated by qRT-PCR. Among these six milRNAs, Ssc-milR-240 was potentially associated with sclerotial development by epigenetic regulation of its target histone acetyltransferase. This study will facilitate the better understanding of the milRNA regulation associated with sclerotial development in S. sclerotiorum and even other sclerotia-forming fungi. This work will provide novel insights into the molecular regulations of fungal morphogenesis and the candidate targets of milRNAs used for the sustainable management of plant diseases caused by S. sclerotiorum.

Chilling Stress Triggers VvAgo1-Mediated miRNA-Like RNA Biogenesis in Volvariella volvacea

In Volvariella volvacea, an important species of edible mushroom, cryogenic autolysis is a typical phenomenon that occurs during abnormal metabolism. Analysis of gene expression profiling and qPCR showed that chilling stress (CS) significantly and continuously upregulated only one type of Argonaute in V. volvacea, i.e., VvAgo1. Structural and evolutionary analysis revealed that VvAgo1 belongs to the Ago-like family, and its evolution has involved gene duplication, subsequent gene loss, and purifying selection. Analysis of its interaction network and expression suggested that CS triggers VvAgo1-mediated miRNA-like RNA (milRNA) biogenesis in V. volvacea V23 but not in VH3 (a composite mutant strain from V23 with improved CS resistance). Small RNA sequencing and qPCR analysis confirmed that CS triggered the increased milRNA expression in V23 and not in VH3. The predicted target genes of the increased milRNAs were enriched in several pathways, such as signal transduction and ubiquitination. Heatmap analysis showed that CS altered the expression profile of milRNAs with their target genes related to signal transduction and ubiquitination in V23. Combined analysis of transcriptome and proteome data confirmed that most of the target genes of the increased milRNAs were not translated into proteins. Our observations indicate that CS might trigger VvAgo1-mediated RNAi to facilitate the cryogenic autolysis of V. volvacea.

Adaptive regulation of virulence genes by microRNA-like RNAs in Valsa mali

MicroRNAs play important roles in the regulation of gene expression in plants and animals. However, little information is known about the action mechanism and function of fungal microRNA-like RNAs (milRNAs). In this study, combining deep sequencing, molecular and histological assays, milRNAs and their targets in the phytopathogenic fungus Valsa mali were isolated and identified. A critical milRNA, Vm-milR16, was identified to adaptively regulate the expression of virulence genes. Fourteen isolated milRNAs showed high expression abundance. Based on the assessment of a pathogenicity function of these milRNAs, Vm-milR16 was found to be a critical milRNA in V. mali by regulating sucrose non-fermenting 1 (VmSNF1), 4,5-DOPA dioxygenase extradiol (VmDODA), and a hypothetical protein (VmHy1). During V. mali infection, Vm-milR16 is downregulated, while its targets are upregulated. Overexpression of Vm-milR16, but not mutated Vm-milR16, significantly reduces the expression of targets and virulence of V. mali. Furthermore, deletion of VmSNF1, VmDODA and VmHy1 significantly reduce virulence of V. mali. All three targets seem to be essential for oxidative stress response and VmSNF1 is required for expression of pectinase genes during V. mali-host interaction. Our results demonstrate Vm-milRNAs contributing to the infection of V. mali on apple trees by adaptively regulating virulence genes.

Identification and characterization of Dicer-like genes in leaf rust pathogen (Puccinia triticina) of wheat

Puccinia triticina (P. triticina) is one of the most devastating fungal pathogens of wheat which causes significant annual yield loss to the crop. Understanding the gene regulatory mechanism of the biotrophic pathogen is one of the important aspects of host-pathogen interaction studies. Dicer-like genes are considered as important mediators of RNAi-based gene regulation. In this study, we report the presence of three Dicer-like genes (Pt-DCL1, Pt-DCL2, Pt-DCL3) in P. triticina genome identified through computational and biological analyses. Quantitative real-time PCR studies revealed an increase in the expression of these genes in germinating spore stages. Heterologous expression combined with mass spectrometry analysis of Pt-DCL2 confirmed the presence of a canonical Dicer-like gene in P. triticina. Phylogenetic analysis of the Pt-DCLs with the Dicer-like proteins from other organisms showed a distinct cluster of rust pathogens from the order Pucciniales. The results indicated a species-specific duplication of Dicer-like genes within the wheat rust pathogens. This study, for the first time, reports the presence of Dicer-dependent RNAi pathway in P. triticina that may play a role in gene regulatory mechanism of the pathogen during its development. Our study serves as a vital source of information for further RNAi-based molecular studies for better understanding and management of the wheat leaf rust disease.

RNA-Spray-Mediated Silencing of Fusarium graminearum AGO and DCL Genes Improve Barley Disease Resistance

Over the last decade, several studies have revealed the enormous potential of RNA-silencing strategies as a potential alternative to conventional pesticides for plant protection. We have previously shown that targeted gene silencing mediated by an in planta expression of non-coding inhibitory double-stranded RNAs (dsRNAs) can protect host plants against various diseases with unprecedented efficiency. In addition to the generation of RNA-silencing (RNAi) signals in planta, plants can be protected from pathogens, and pests by spray-applied RNA-based biopesticides. Despite the striking efficiency of RNA-silencing-based technologies holds for agriculture, the molecular mechanisms underlying spray-induced gene silencing (SIGS) strategies are virtually unresolved, a requirement for successful future application in the field. Based on our previous work, we predict that the molecular mechanism of SIGS is controlled by the fungal-silencing machinery. In this study, we used SIGS to compare the silencing efficiencies of computationally-designed vs. manually-designed dsRNA constructs targeting ARGONAUTE and DICER genes of Fusarium graminearum (Fg). We found that targeting key components of the fungal RNAi machinery via SIGS could protect barley leaves from Fg infection and that the manual design of dsRNAs resulted in higher gene-silencing efficiencies than the tool-based design. Moreover, our results indicate the possibility of cross-kingdom RNA silencing in the Fg-barley interaction, a phenomenon in which sRNAs operate as effector molecules to induce gene silencing between species from different kingdoms, such as a plant host and their interacting pathogens.

Facilitative and synergistic interactions between fungal and plant viruses

Plants and fungi are closely associated through parasitic or symbiotic relationships in which bidirectional exchanges of cellular contents occur. Recently, a plant virus was shown to be transmitted from a plant to a fungus, but it is unknown whether fungal viruses can also cross host barriers and spread to plants. In this study, we investigated the infectivity of Cryphonectria hypovirus 1 (CHV1, family Hypoviridae), a capsidless, positive-sense (+), single-stranded RNA (ssRNA) fungal virus in a model plant, Nicotiana tabacum CHV1 replicated in mechanically inoculated leaves but did not spread systemically, but coinoculation with an unrelated plant (+)ssRNA virus, tobacco mosaic virus (TMV, family Virgaviridae), or other plant RNA viruses, enabled CHV1 to systemically infect the plant. Likewise, CHV1 systemically infected transgenic plants expressing the TMV movement protein, and coinfection with TMV further enhanced CHV1 accumulation in these plants. Conversely, CHV1 infection increased TMV accumulation when TMV was introduced into a plant pathogenic fungus, Fusarium graminearum In the in planta F. graminearum inoculation experiment, we demonstrated that TMV infection of either the plant or the fungus enabled the horizontal transfer of CHV1 from the fungus to the plant, whereas CHV1 infection enhanced fungal acquisition of TMV. Our results demonstrate two-way facilitative interactions between the plant and fungal viruses that promote cross-kingdom virus infections and suggest the presence of plant-fungal-mediated routes for dissemination of fungal and plant viruses in nature.

Small RNA Bidirectional Crosstalk During the Interaction Between Wheat and Zymoseptoria tritici

Cross-kingdom RNA interference (RNAi) has been shown to play important roles during plant-pathogen interactions, and both plants and pathogens can use small RNAs (sRNAs) to silence genes in each other. This bidirectional cross-kingdom RNAi was still unexplored in the wheat-Zymoseptoria tritici pathosystem. Here, we performed a detailed analysis of the sRNA bidirectional crosstalk between wheat and Z. tritici. Using a combination of small RNA sequencing (sRNA-seq) and microRNA sequencing (mRNA-seq), we were able to identify known and novel sRNAs and study their expression and their action on putative targets in both wheat and Z. tritici. We predicted the target genes of all the sRNAs in either wheat or Z. tritici transcriptome and used degradome analysis to validate the cleavage of these gene transcripts. We could not find any clear evidence of a cross-kingdom RNAi acting by mRNA cleavage in this pathosystem. We also found that the fungal sRNA enrichment was lower in planta than during in vitro growth, probably due to the lower expression of the only Dicer gene of the fungus during plant infection. Our results support the recent finding that Z. tritici sRNAs cannot play important roles during wheat infection. However, we also found that the fungal infection induced wheat sRNAs regulating the expression of specific wheat genes, including auxin-related genes, as an immune response. These results indicate a role of sRNAs in the regulation of wheat defenses during Z. tritici infection. Our findings contribute to improve our understanding of the interactions between wheat and Z. tritici.

RNA-Spray-Mediated Silencing of Fusarium graminearum AGO and DCL Genes Improve Barley Disease Resistance

Over the last decade, several studies have revealed the enormous potential of RNA-silencing strategies as a potential alternative to conventional pesticides for plant protection. We have previously shown that targeted gene silencing mediated by an in planta expression of non-coding inhibitory double-stranded RNAs (dsRNAs) can protect host plants against various diseases with unprecedented efficiency. In addition to the generation of RNA-silencing (RNAi) signals in planta, plants can be protected from pathogens, and pests by spray-applied RNA-based biopesticides. Despite the striking efficiency of RNA-silencing-based technologies holds for agriculture, the molecular mechanisms underlying spray-induced gene silencing (SIGS) strategies are virtually unresolved, a requirement for successful future application in the field. Based on our previous work, we predict that the molecular mechanism of SIGS is controlled by the fungal-silencing machinery. In this study, we used SIGS to compare the silencing efficiencies of computationally-designed vs. manually-designed dsRNA constructs targeting ARGONAUTE and DICER genes of Fusarium graminearum (Fg). We found that targeting key components of the fungal RNAi machinery via SIGS could protect barley leaves from Fg infection and that the manual design of dsRNAs resulted in higher gene-silencing efficiencies than the tool-based design. Moreover, our results indicate the possibility of cross-kingdom RNA silencing in the Fg-barley interaction, a phenomenon in which sRNAs operate as effector molecules to induce gene silencing between species from different kingdoms, such as a plant host and their interacting pathogens.

Different Components of the RNA Interference Machinery Are Required for Conidiation, Ascosporogenesis, Virulence, Deoxynivalenol Production, and Fungal Inhibition by Exogenous Double-Stranded RNA in the Head Blight Pathogen Fusarium graminearum

In filamentous fungi, gene silencing through RNA interference (RNAi) shapes many biological processes, including pathogenicity. We explored the requirement of key components of fungal RNAi machineries, including DICER-like 1 and 2 (FgDCL1, FgDCL2), ARGONAUTE 1 and 2 (FgAGO1, FgAGO2), AGO-interacting protein FgQIP (QDE2-interacting protein), RecQ helicase (FgQDE3), and four RNA-dependent RNA polymerases (FgRdRP1, FgRdRP2, FgRdRP3, FgRdRP4), in the ascomycete mycotoxin-producing fungal pathogen Fusarium graminearum (Fg) for sexual and asexual multiplication, pathogenicity, and its sensitivity to double-stranded (ds)RNA. We corroborate and extend earlier findings that conidiation, ascosporogenesis, and Fusarium head blight (FHB) symptom development require an operable RNAi machinery. The involvement of RNAi in conidiation is dependent on environmental conditions as it is detectable only under low light (<2 μmol m⁻² s⁻¹). Although both DCLs and AGOs partially share their functions, the sexual ascosporogenesis is mediated primarily by FgDCL1 and FgAGO2, while FgDCL2 and FgAGO1 contribute to asexual conidia formation and germination. FgDCL1 and FgAGO2 also account for pathogenesis as their knockout (KO) results in reduced FHB development. Apart from KO mutants Δdcl2 and Δago1, mutants Δrdrp2, Δrdrp3, Δrdrp4, Δqde3, and Δqip are strongly compromised for conidiation, while KO mutations in all RdPRs, QDE3, and QIP strongly affect ascosporogenesis. Analysis of trichothecenes mycotoxins in wheat kernels showed that the relative amount of deoxynivalenol (DON), calculated as [DON] per amount of fungal genomic DNA was reduced in all spikes infected with RNAi mutants, suggesting the possibility that the fungal RNAi pathways affect Fg’s DON production. Moreover, silencing of fungal genes by exogenous target gene-specific double-stranded RNA (dsRNA) (spray-induced gene silencing, SIGS) is dependent on DCLs, AGOs, and QIP, but not on QDE3. Together these data show that in F. graminearum, different key components of the RNAi machinery are crucial in different steps of fungal development and pathogenicity.

sRNA Profiling Combined With Gene Function Analysis Reveals a Lack of Evidence for Cross-Kingdom RNAi in the Wheat - Zymoseptoria tritici Pathosystem

Cross-kingdom small RNA (sRNA) silencing has recently emerged as a mechanism facilitating fungal colonization and disease development. Here we characterized RNAi pathways in Zymoseptoria tritici, a major fungal pathogen of wheat, and assessed their contribution to pathogenesis. Computational analysis of fungal sRNA and host mRNA sequencing datasets was used to define the global sRNA populations in Z. tritici and predict their mRNA targets in wheat. 389 in planta-induced sRNA loci were identified. sRNAs generated from some of these loci were predicted to target wheat mRNAs including those potentially involved in pathogen defense. However, molecular approaches failed to validate targeting of selected wheat mRNAs by fungal sRNAs. Mutant strains of Z. tritici carrying deletions of genes encoding key components of RNAi such as Dicer-like (DCL) and Argonaute (AGO) proteins were generated, and virulence bioassays suggested that these are dispensable for full infection of wheat. Nonetheless, our results did suggest the existence of non-canonical DCL-independent pathway(s) for sRNA biogenesis in Z. tritici. dsRNA targeting essential fungal genes applied in vitro or generated from an RNA virus vector in planta in a procedure known as HIGS (Host-Induced Gene Silencing) was ineffective in preventing Z. tritici growth or disease. We also demonstrated that Z. tritici is incapable of dsRNA uptake. Collectively, our data suggest that RNAi approaches for gene function analyses in this fungal species and potentially also as a control measure may not be as effective as has been demonstrated for some other plant pathogenic fungi.

Nonproteinaceous effectors: the terra incognita of plant-fungal interactions

Molecular plant-fungal interaction studies have mainly focused on small secreted protein effectors. However, accumulating evidence shows that numerous fungal secondary metabolites are produced at all stages of plant colonization, especially during early asymptomatic/biotrophic phases. The discovery of fungal small RNAs targeting plant transcripts has expanded the fungal repertoire of nonproteinaceous effectors even further. The challenge now is to develop specific functional methods to fully understand the biological roles of these effectors. Studies on fungal extracellular vesicles are also needed because they could be the universal carriers of all kinds of fungal effectors. With this review, we aim to stimulate the nonproteinaceous effector research field to move from descriptive to functional studies, which should bring a paradigm shift in plant-fungal interactions.

Symptomatic plant viroid infections in phytopathogenic fungi

Viroids are pathogenic agents that have a small, circular noncoding RNA genome. They have been found only in plant species; therefore, their infectivity and pathogenicity in other organisms remain largely unexplored. In this study, we investigate whether plant viroids can replicate and induce symptoms in filamentous fungi. Seven plant viroids representing viroid groups that replicate in either the nucleus or chloroplast of plant cells were inoculated to three plant pathogenic fungi, Cryphonectria parasitica, Valsa mali, and Fusarium graminearum By transfection of fungal spheroplasts with viroid RNA transcripts, each of the three, hop stunt viroid (HSVd), iresine 1 viroid, and avocado sunblotch viroid, can stably replicate in at least one of those fungi. The viroids are horizontally transmitted through hyphal anastomosis and vertically through conidia. HSVd infection severely debilitates the growth of V. mali but not that of the other two fungi, while in F. graminearum and C. parasitica, with deletion of dicer-like genes, the primary components of the RNA-silencing pathway, HSVd accumulation increases. We further demonstrate that HSVd can be bidirectionally transferred between F. graminearum and plants during infection. The viroids also efficiently infect fungi and induce disease symptoms when the viroid RNAs are exogenously applied to the fungal mycelia. These findings enhance our understanding of viroid replication, host range, and pathogenicity, and of their potential spread to other organisms in nature.

Degradation of Fungal MicroRNAs Triggered by Short Tandem Target Mimics Is via the Small-RNA-Degrading Nuclease

MicroRNAs (miRNAs) have been recognized as sequence-specific regulators of the genome, transcriptome, and proteome in eukaryotes. However, the functions and working mechanisms of hundreds of fungal miRNA-like (miR-like) RNAs are obscure. Here, we report that a short tandem target mimic (STTM) triggered the degradation of several fungal miR-like RNAs in two different fungal species, Metarhizium robertsii and Aspergillus flavus, and that small-RNA-degrading nucleases (SDNs) were indispensable for such degradation. STTMs were most effective when the fungal polymerase II (Pol II) promoter was used for their expression, while the Pol III promoter was less effective. The length of the STTM spacer, approximately 48 to 96 nucleotides, and the number of miR-like RNA binding sites, from 2 to 4 copies, showed no significant difference in the degradation of miR-like RNAs. STTMs modulated the miR-like RNA expression levels in at least two different fungal species, which further impacted fungal asexual growth and sporulation. Further analysis showed that the degraded miR-like RNAs in STTM mutants led to the upregulation of potential target genes involved in fungal development and conidial production, which result in different phenotypes in these mutants. The STTM technology developed in this study is an effective and powerful tool for the functional dissection of fungal miR-like RNAs.IMPORTANCE The development and application of STTM technology to block miR-like RNAs in M. robertsii and A. flavus may allow for efficient generation of miR-like RNA mutants in various fungi, providing a powerful tool for functional genomics of small RNA molecules in fungi.

milR4 and milR16 Mediated Fruiting Body Development in the Medicinal Fungus Cordyceps militaris

Cordyceps militaris readily performs sexual reproduction, thus providing a remarkably rich model for understanding the processes involved in sexual development. It could regulate expression of human genes by diet-derived miRNA-like RNAs (milRNAs). However, the study of miRNAs in C. militaris has been limited. In the present study, genes encoding Dicers, Argonautes, and RNA-dependent RNA polymerases were identified. Illumina deep sequencing was performed to characterize the milRNAs in C. militaris at asexual and sexual development stages. Total 38 milRNAs were identified and five milRNAs were validated by northern blot and qRT-PCR, out of which, 19 were specific for sexual development. Importantly, the fungi could not form fruiting bodies after disruption of milR4, while the perithecium was formed in advance after over-expression of milR4. Abnormal pale yellow fruiting body primordium, covered with abnormal primordium, was formed in the strain with miR16 disruption. Although no milR4 or milR16 target genes were identified, differential expression of many different genes involved in mycelium growth and sexual development (mating process, mating signaling, and fruiting body development) among these mutants were found. Overall, milRNAs play vital roles in sexual development in C. militaris.