Comparative Small RNA and Degradome Sequencing Provides Insights into Antagonistic Interactions in the Biocontrol Fungus Clonostachys rosea.

Necrotrophic mycoparasitism is an intricate process involving recognition, physical mycelial contact, and killing of host fungi (mycohosts). During such interactions, mycoparasites undergo a complex developmental process involving massive regulatory changes of gene expression to produce a range of chemical compounds and proteins that contribute to the parasitism of the mycohosts. Small RNAs (sRNAs) are vital components of posttranscriptional gene regulation, although their role in gene expression regulation during mycoparasitisms remain understudied. Here, we investigated the role of sRNA-mediated gene regulation in mycoparasitism by performing sRNA and degradome tag sequencing of the mycoparasitic fungus Clonostachys rosea interacting with the plant-pathogenic mycohosts Botrytis cinerea and Fusarium graminearum at two time points. The majority of differentially expressed sRNAs were downregulated during the interactions with the mycohosts compared to a C. rosea self-interaction control, thus allowing desuppression (upregulation) of mycohost-responsive genes. Degradome analysis showed a positive correlation between high degradome counts and antisense sRNA mapping and led to the identification of 201 sRNA-mediated potential gene targets for 282 differentially expressed sRNAs. Analysis of sRNA potential gene targets revealed that the regulation of genes coding for membrane proteins was a common response against both mycohosts. The regulation of genes involved in oxidative stress tolerance and cellular metabolic and biosynthetic processes was exclusive against F. graminearum, highlighting common and mycohost-specific gene regulation of C. rosea. By combining these results with transcriptome data collected during a previous study, we expand the understanding of the role of sRNA in regulating interspecific fungal interactions and mycoparasitism. IMPORTANCE Small RNAs (sRNAs) are emerging as key players in pathogenic and mutualistic fungus-plant interactions; however, their role in fungus-fungus interactions remains elusive. In this study, we employed the necrotrophic mycoparasite Clonostachys rosea and the plant-pathogenic mycohosts Botrytis cinerea and Fusarium graminearum and investigated the sRNA-mediated gene regulation in mycoparasitic interactions. The combined approach of sRNA and degradome tag sequencing identified 201 sRNA-mediated putative gene targets for 282 differentially expressed sRNAs, highlighting the role of sRNA-mediated regulation of mycoparasitism in C. rosea. We also identified 36 known and 13 novel microRNAs (miRNAs) and their potential gene targets at the endogenous level and at a cross-species level in B. cinerea and F. graminearum, indicating a role of cross-species RNA interference (RNAi) in mycoparasitism, representing a novel mechanism in biocontrol interactions. Furthermore, we showed that C. rosea adapts its transcriptional response, and thereby its interaction mechanisms, based on the interaction stages and identity of the mycohost.

Role of Dicer-Dependent RNA Interference in Regulating Mycoparasitic Interactions.

Dicer-like proteins (DCLs) play a vital role in RNA interference (RNAi), by cleaving RNA filament into small RNAs. Although DCL-mediated RNAi can regulate interspecific communication between pathogenic/mutualistic organisms and their hosts, its role in mycoparasitic interactions is yet to be investigated. In this study, we deleted dcl genes in the mycoparasitic fungus Clonostachys rosea and characterize the functions of DCL-dependent RNAi in mycoparasitism. Deletion of dcl2 resulted in a mutant with reduced secondary metabolite production, antagonism toward the plant-pathogenic fungus Botrytis cinerea, and reduced ability to control Fusarium foot rot disease on wheat, caused by Fusarium graminearum. Transcriptome sequencing of the in vitro interaction between the C. rosea Δdcl2 strain and B. cinerea or F. graminearum identified the downregulation of genes coding for transcription factors, membrane transporters, hydrolytic enzymes, and secondary metabolites biosynthesis enzymes putatively involved in antagonistic interactions, in comparison with the C. rosea wild-type interaction. A total of 61 putative novel microRNA-like RNAs (milRNAs) were identified in C. rosea, and 11 were downregulated in the Δdcl2 mutant. In addition to putative endogenous gene targets, these milRNAs were predicted to target B. cinerea and F. graminearum virulence factor genes, which showed an increased expression during interaction with the Δdcl2 mutant incapable of producing the targeting milRNAs. In summary, this study constitutes the first step in elucidating the role of RNAi in mycoparasitic interactions, with important implications for biological control of plant diseases, and poses the base for future studies focusing on the role of cross-species RNAi regulating mycoparasitic interactions. IMPORTANCE Small RNAs mediated RNA interference (RNAi) known to regulate several biological processes. Dicer-like endoribonucleases (DCLs) play a vital role in the RNAi pathway by generating sRNAs. In this study, we investigated a role of DCL-mediated RNAi in interference interactions between mycoparasitic fungus Clonostachys rosea and the two fungal pathogens Botrytis cinerea and Fusarium graminearum (here called mycohosts). We found that the dcl mutants were not able to produce 11 sRNAs predicted to finetune the regulatory network of genes known to be involved in production of hydrolytic enzymes, antifungal compounds, and membrane transporters needed for antagonistic action of C. rosea. We also found C. rosea sRNAs putatively targeting known virulence factors in the mycohosts, indicating RNAi-mediated cross-species communication. Our study expanded the understanding of underlying mechanisms of cross-species communication during interference interactions and poses a base for future works studying the role of DCL-based cross-species RNAi in fungal interactions.

Comparative genomics highlights the importance of drug efflux transporters during evolution of mycoparasitism in Clonostachys subgenus Bionectria (Fungi, Ascomycota, Hypocreales).

Various strains of the mycoparasitic fungal species Clonostachys rosea are used commercially as biological control agents for the control of fungal plant diseases in agricultural crop production. Further improvements of the use and efficacy of C. rosea in biocontrol require a mechanistic understanding of the factors that determines the outcome of the interaction between C. rosea and plant pathogenic fungi. Here, we determined the genome sequences of 11 Clonostachys strains, representing five species in Clonostachys subgenus Bionectria, and performed a comparative genomic analysis with the aim to identify gene families evolving under selection for gene gains or losses. Several gene families predicted to encode proteins involved in biosynthesis of secondary metabolites, including polyketide synthases, nonribosomal peptide syntethases and cytochrome P450s, evolved under selection for gene gains (p ≤ .05) in the Bionectria subgenus lineage. This was accompanied with gene copy number increases (p ≤ .05) in ATP-binding cassette (ABC) transporters and major facilitator superfamily (MFS) transporters predicted to contribute to drug efflux. Most Clonostachys species were also characterized by high numbers of auxiliary activity (AA) family 9 lytic polysaccharide monooxygenases, AA3 glucose-methanol-choline oxidoreductases and additional carbohydrate-active enzyme gene families with putative activity (or binding) towards xylan and rhamnose/pectin substrates. Particular features of the C. rosea genome included expansions (p ≤ .05) of the ABC-B4 multidrug resistance transporters, the ABC-C5 multidrug resistance-related transporters and the 2.A.1.3 drug:H + antiporter-2 MFS drug resistance transporters. The ABC-G1 pleiotropic drug resistance transporter gene abcG6 in C. rosea was induced (p ≤ .009) by exposure to the antifungal Fusarium mycotoxin zearalenone (1121-fold) and various fungicides. Deletion of abcG6 resulted in mutants with reduced (p < .001) growth rates on media containing the fungicides boscalid, fenhexamid and iprodione. Our results emphasize the role of biosynthesis of, and protection against, secondary metabolites in Clonostachys subgenus Bionectria.

Natural variation of root lesion nematode antagonism in the biocontrol fungus Clonostachys rosea and identification of biocontrol factors through genome-wide association mapping.

Biological control is a promising approach to reduce plant diseases caused by nematodes to ensure high productivity in agricultural production. Large-scale analyses of genetic variation in fungal species used for biocontrol can generate knowledge regarding interaction mechanisms that can improve efficacy of biocontrol applications. In this study, we performed a genome-wide association study (GWAS) for in vitro antagonism against the root lesion nematode Pratylenchus penetrans in 53 previously genome re-sequenced strains of the biocontrol fungus Clonostachys rosea. Nematode mortality in C. rosea potato dextrose broth (PDB) culture filtrates was highly variable and showed continuous variation (p < .001) between strains, indicating a polygenic inheritance. Twenty-one strains produced culture filtrates with higher (p ≤ .05) nematode mortality compared with the PDB control treatment, while ten strains lowered (p ≤ .05) the mortality. The difference in in vitro antagonism against P. penetrans correlated with antagonism against the soybean cyst nematode Heterodera glycines, indicating lack of host specificity in C. rosea. An empirical Bayesian multiple hypothesis testing approach identified 279 single nucleotide polymorphism markers significantly (local false sign rate < 10-10) associated with the trait. Genes present in the genomic regions associated with nematicidal activity included several membrane transporters, a chitinase and genes encoding proteins predicted to biosynthesize secondary metabolites. Gene deletion strains of the predicted nonribosomal peptide synthetase genes nps4 and nps5 were generated and showed increased (p ≤ .001) fungal growth and conidiation rates compared to the wild type. Deletion strains also exhibited reduced (p < .001) nematicidal activity and reduced (p ≤ .05) biocontrol efficacy against nematode root disease and against fusarium foot rot on wheat. In summary, we show that the GWAS approach can be used to identify biocontrol factors in C. rosea, specifically the putative nonribosomal peptide synthetases NPS4 and NPS5.

Deciphering common and specific transcriptional immune responses in pea towards the oomycete pathogens Aphanomyces euteiches and Phytophthora pisi.

BACKGROUND: Root rot caused by Aphanomyces euteiches is one of the most destructive pea diseases while a distantly related species P. pisi has been recently described as the agent of pea and faba bean root rot. These two oomycete pathogens with different pathogenicity factor repertories have both evolved specific mechanisms to infect pea. However, little is known about the genes and mechanisms of defence against these pathogens in pea. In the present study, the transcriptomic response of pea to these two pathogens was investigated at two time points during early phase of infection using a Medicago truncatula microarray. RESULTS: Of the 37,976 genes analysed, 574 and 817 were differentially expressed in response to A. euteiches at 6 hpi and 20 hpi, respectively, while 544 and 611 genes were differentially regulated against P. pisi at 6 hpi and 20 hpi, respectively. Differentially expressed genes associated with plant immunity responses were involved in cell wall reinforcement, hormonal signalling and phenylpropanoid metabolism. Activation of cell wall modification, regulation of jasmonic acid biosynthesis and induction of ethylene signalling pathway were among the common transcriptional responses to both of these oomycetes. However, induction of chalcone synthesis and the auxin pathway were specific transcriptional changes against A. euteiches. CONCLUSIONS: Our results demonstrate a global view of differentially expressed pea genes during compatible interactions with P. pisi and A. euteiches at an early phase of infection. The results suggest that distinct signalling pathways are triggered in pea by these two pathogens that lead to common and specific immune mechanisms in response to these two oomycetes. The generated knowledge may eventually be used in breeding pea varieties with resistance against root rot disease.

Biocontrol agents efficiently inhibit sporulation of Botrytis aclada on necrotic leaf tips but spread to adjacent living tissue is not prevented.

Ulocladium atrum (isolates 385 and 302) consistently inhibited Botrytis aclada sporulation on dead onion leaf pieces under constant moist conditions and with an interrupted wetness period of 9 h. Clonostachys rosea (isolate 201) was as effective as U. atrum under constant moist conditions, but was ineffective if exposed to a drying period. No sporulation of B. aclada was observed 8 and 12 days after inoculation in the presence of U. atrum 302. C. rosea 201 significantly reduced B. aclada sporulation 8 days, but not 12 days after inoculation. When U. atrum 302 or C. rosea 201 was applied 1 day prior to B. aclada the antagonistic effect was higher compared to when the antagonists were applied on the same day. C. rosea 201 and U. atrum 302 did not obstruct the growth of B. aclada from necrotic onion leaf tips into living tissue, when artificially induced necrotic leaf tips were infested with B. aclada 24 h prior to antagonists. Three days after antagonist application, no symptoms could be observed on the healthy leaf tissue, nor was there sporulation on the necrotic leaf tip. However, B. aclada was immunologically detected 2 cm below the inoculation site. We conclude that under constant moist conditions the antagonists C. rosea 201 and U. atrum 302 cannot stop the progress of B. aclada from necrotic into fresh leaf tissue.