Transcriptomic and functional analyses on a Botrytis cinerea multidrug-resistant (MDR) strain provides new insights into the potential molecular mechanisms of MDR and fitness.

Botrytis cinerea is a notorious pathogen causing pre- and post-harvest spoilage in many economically important crops. Excessive application of site-specific fungicides to control the pathogen has led to the selection of strains possessing target site alterations associated with resistance to these fungicides and/or strains overexpressing efflux transporters associated with multidrug resistance (MDR). MDR in B. cinerea has been correlated with the overexpression of atrB and mfsM2, encoding an ATP-binding cassette (ABC) and a major facilitator superfamily (MFS) transporter, respectively. However, it remains unknown whether other transporters may also contribute to the MDR phenotype. In the current study, the transcriptome of a B. cinerea multidrug-resistant (MDR) field strain was analysed upon exposure to the fungicide fludioxonil, and compared to the B05.10 reference strain. The transcriptome of this field strain displayed significant differences as compared to B05.10, including genes involved in sugar membrane transport, toxin production and virulence. Among the induced genes in the field strain, even before exposure to fludioxonil, were several putatively encoding ABC and MFS transmembrane transporters. Overexpression of a highly induced MFS transporter gene in the B05.10 strain led to an increased tolerance to the fungicides fluopyram and boscalid, indicating an involvement in efflux transport of these compounds. Overall, the data from this study give insights towards better understanding the molecular mechanisms involved in MDR and fitness cost, contributing to the development of more efficient control strategies against this pathogen.

Transcriptomic analysis identifies candidate genes for Aphanomyces root rot disease resistance in pea.

BACKGROUND: Aphanomyces euteiches is a soil-borne oomycete that causes root rot in pea and other legume species. Symptoms of Aphanomyces root rot (ARR) include root discoloration and wilting, leading to significant yield losses in pea production. Resistance to ARR is known to be polygenic but the roles of single genes in the pea immune response are still poorly understood. This study uses transcriptomics to elucidate the immune response of two pea genotypes varying in their levels of resistance to A. euteiches. RESULTS: In this study, we inoculated roots of the pea (P. sativum L.) genotypes 'Linnea' (susceptible) and 'PI180693' (resistant) with two different A. euteiches strains varying in levels of virulence. The roots were harvested at 6 h post-inoculation (hpi), 20 hpi and 48 hpi, followed by differential gene expression analysis. Our results showed a time- and genotype-dependent immune response towards A. euteiches infection, involving several WRKY and MYB-like transcription factors, along with genes associated with jasmonic acid (JA) and abscisic acid (ABA) signaling. By cross-referencing with genes segregating with partial resistance to ARR, we identified 39 candidate disease resistance genes at the later stage of infection. Among the genes solely upregulated in the resistant genotype 'PI180693', Psat7g091800.1 was polymorphic between the pea genotypes and encoded a Leucine-rich repeat receptor-like kinase reminiscent of the Arabidopsis thaliana FLAGELLIN-SENSITIVE 2 receptor. CONCLUSIONS: This study provides new insights into the gene expression dynamics controlling the immune response of resistant and susceptible pea genotypes to A. euteiches infection. We present a set of 39 candidate disease resistance genes for ARR in pea, including the putative immune receptor Psat7g091800.1, for future functional validation.

Verticillium longisporum phospholipase VlsPLA2 is a virulence factor that targets host nuclei and modulates plant immunity.

Phospholipase A2 (PLA2 ) is a lipolytic enzyme that hydrolyses phospholipids in the cell membrane. In the present study, we investigated the role of secreted PLA2 (VlsPLA2 ) in Verticillium longisporum, a fungal phytopathogen that mostly infects plants belonging to the Brassicaceae family, causing severe annual yield loss worldwide. Expression of the VlsPLA2 gene, which encodes active PLA2 , is highly induced during the interaction of the fungus with the host plant Brassica napus. Heterologous expression of VlsPLA2 in Nicotiana benthamiana resulted in increased synthesis of certain phospholipids compared to plants in which enzymatically inactive PLA2 was expressed (VlsPLA2 ΔCD ). Moreover, VlsPLA2 suppresses the hypersensitive response triggered by the Cf4/Avr4 complex, thereby suppressing the chitin-induced reactive oxygen species burst. VlsPLA2 -overexpressing V. longisporum strains showed increased virulence in Arabidopsis plants, and transcriptomic analysis of this fungal strain revealed that the induction of the gene contributed to increased virulence. VlsPLA2 was initially localized to the host nucleus and then translocated to the chloroplasts at later time points. In addition, VlsPLA2 bound to the vesicle-associated membrane protein A (VAMPA) and was transported to the nuclear membrane. In the nucleus, VlsPLA2 caused major alterations in the expression levels of genes encoding transcription factors and subtilisin-like proteases, which play a role in plant immunity. In conclusion, our study showed that VlsPLA2 acts as a virulence factor, possibly by hydrolysing host nuclear envelope phospholipids, which, through a signal transduction cascade, may suppress basal plant immune responses.

RNA silencing proteins and small RNAs in oomycete plant pathogens and biocontrol agents.

Introduction: Oomycetes cause several damaging diseases of plants and animals, and some species also act as biocontrol agents on insects, fungi, and other oomycetes. RNA silencing is increasingly being shown to play a role in the pathogenicity of Phytophthora species, either through trans-boundary movement of small RNAs (sRNAs) or through expression regulation of infection promoting effectors. Methods: To gain a wider understanding of RNA silencing in oomycete species with more diverse hosts, we mined genome assemblies for Dicer-like (DCL), Argonaute (AGO), and RNA dependent RNA polymerase (RDRP) proteins from Phytophthora plurivora, Ph. cactorum, Ph. colocasiae, Pythium oligandrum, Py. periplocum, and Lagenidium giganteum. Moreover, we sequenced small RNAs from the mycelium stage in each of these species. Results and discussion: Each of the species possessed a single DCL protein, but they differed in the number and sequence of AGOs and RDRPs. SRNAs of 21nt, 25nt, and 26nt were prevalent in all oomycetes analyzed, but the relative abundance and 5' base preference of these classes differed markedly between genera. Most sRNAs mapped to transposons and other repeats, signifying that the major role for RNA silencing in oomycetes is to limit the expansion of these elements. We also found that sRNAs may act to regulate the expression of duplicated genes. Other sRNAs mapped to several gene families, and this number was higher in Pythium spp., suggesting a role of RNA silencing in regulating gene expression. Genes for most effector classes were the source of sRNAs of variable size, but some gene families showed a preference for specific classes of sRNAs, such as 25/26 nt sRNAs targeting RxLR effector genes in Phytophthora species. Novel miRNA-like RNAs (milRNAs) were discovered in all species, and two were predicted to target transcripts for RxLR effectors in Ph. plurivora and Ph. cactorum, indicating a putative role in regulating infection. Moreover, milRNAs from the biocontrol Pythium species had matches in the predicted transcriptome of Phytophthora infestans and Botrytis cinerea, and L. giganteum milRNAs matched candidate genes in the mosquito Aedes aegypti. This suggests that trans-boundary RNA silencing may have a role in the biocontrol action of these oomycetes.

Insights into the ecological generalist lifestyle of Clonostachys fungi through analysis of their predicted secretomes.

Introduction: The fungal secretome comprise diverse proteins that are involved in various aspects of fungal lifestyles, including adaptation to ecological niches and environmental interactions. The aim of this study was to investigate the composition and activity of fungal secretomes in mycoparasitic and beneficial fungal-plant interactions. Methods: We used six Clonostachys spp. that exhibit saprotrophic, mycotrophic and plant endophytic lifestyles. Genome-wide analyses was performed to investigate the composition, diversity, evolution and gene expression of Clonostachys secretomes in relation to their potential role in mycoparasitic and endophytic lifestyles. Results and discussion: Our analyses showed that the predicted secretomes of the analyzed species comprised between 7 and 8% of the respective proteomes. Mining of transcriptome data collected during previous studies showed that 18% of the genes encoding predicted secreted proteins were upregulated during the interactions with the mycohosts Fusarium graminearum and Helminthosporium solani. Functional annotation of the predicted secretomes revealed that the most represented protease family was subclass S8A (11-14% of the total), which include members that are shown to be involved in the response to nematodes and mycohosts. Conversely, the most numerous lipases and carbohydrate-active enzyme (CAZyme) groups appeared to be potentially involved in eliciting defense responses in the plants. For example, analysis of gene family evolution identified nine CAZyme orthogroups evolving for gene gains (p ≤ 0.05), predicted to be involved in hemicellulose degradation, potentially producing plant defense-inducing oligomers. Moreover, 8-10% of the secretomes was composed of cysteine-enriched proteins, including hydrophobins, important for root colonization. Effectors were more numerous, comprising 35-37% of the secretomes, where certain members belonged to seven orthogroups evolving for gene gains and were induced during the C. rosea response to F. graminearum or H. solani. Furthermore, the considered Clonostachys spp. possessed high numbers of proteins containing Common in Fungal Extracellular Membranes (CFEM) modules, known for their role in fungal virulence. Overall, this study improves our understanding of Clonostachys spp. adaptation to diverse ecological niches and establishes a basis for future investigation aiming at sustainable biocontrol of plant diseases.

Evaluation of pea genotype PI180693 partial resistance towards aphanomyces root rot in commercial pea breeding.

The cultivation of vining pea (Pisum sativum) faces a major constraint with root rot diseases, caused by a complex of soil-borne pathogens including the oomycetes Aphanomyces euteiches and Phytophtora pisi. Disease resistant commercial varieties are lacking but the landrace PI180693 is used as a source of partial resistance in ongoing pea breeding programs. In this study, the level of resistance and their interaction with A. euteiches virulence levels of six new back-crossed pea breeding lines, deriving from the cross between the susceptible commercial cultivar Linnea and PI180693, were evaluated for their resistance towards aphanomyces root rot in growth chamber and green house tests. Resistance towards mixed infections by A. euteiches and P. pisi and commercial production traits were evaluated in field trials. In growth chamber trials, pathogen virulence levels had a significant effect on plant resistance, as resistance was more consistent against A. euteiches strains exhibiting high or intermediate virulence compared with lowly virulent strains. In fact, line Z1701-1 showed to be significantly more resistant than both parents when inoculated with a lowly virulent strain. In two separate field trials in 2020, all six breeding lines performed equally well as the resistant parent PI180693 at sites only containing A. euteiches, as there were no differences in disease index. In mixed infections, PI180693 exhibited significantly lower disease index scores than Linnea. However, breeding lines displayed higher disease index scores compared with PI180693, indicating higher susceptibility towards P. pisi. Data on seedling emergence from the same field trials suggested that PI180693 was particularly sensitive towards seed decay/damping off disease caused by P. pisi. Furthermore, the breeding lines performed equally well as Linnea in traits important for green pea production, again emphasizing the commercial potential. In summary, we show that the resistance from PI180693 interacts with virulence levels of the pathogen A. euteiches and is less effective towards root rot caused by P. pisi. Our results show the potential use of combining PI180693 partial resistance against aphanomyces root rot with commercially favorable breeding traits in commercial breeding programs.

Computational Analysis of HTS Data and Its Application in Plant Pathology.

High-throughput sequencing is a basic tool of biological research, and it is extensively used in plant pathology projects. Here, we describe how to handle data coming from a variety of sequencing experiments, focusing on the analysis of Illumina reads. We describe how to perform genome assembly and annotation with DNA reads, correctly analyze RNA-seq data to discover differentially expressed genes, handle amplicon sequencing data from microbial communities, and utilize small RNA sequencing data to predict miRNA sequences and their putative targets.

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.

Cellulose amendment promotes P solubilization by Penicillium aculeatum in non-sterilized soil.

Successful application of microbial biofertilizers, such as phosphorus (P) solubilizing fungi to agroecosystems, is constrained from the lack of knowledge about their ecology; for example in terms of how they respond to an external input of carbon (C) to get established in the soil. In two soil incubation experiments we examined the performance of the P solubilizing fungus Penicillium aculeatum in non-sterile and semi-sterile (γ-irradiated) soil with different C and P sources. Results from the first experiment with C sources showed that starch and cellulose generally improved P solubilization by P. aculeatum measured as water extractable P (Pwep), though only significantly in non-sterile soil. This coincided with an increased population density of P. aculeatum measured with a hygromycin B resistant strain of this fungus. Soil respiration used to measure soil microbial activity was overall much higher in treatments with C compounds than without C in both non-sterile and semi-sterile soil. However, soil respiration was highest with cellulose in semi-sterile soil, especially in combination with P. aculeatum. Hence, for the second experiment with P sources (tricalcium phosphate (TCP) and sewage sludge ash) cellulose was used as a C source for P. aculeatum growth in all treatments. Main results showed that P. aculeatum in combination with cellulose soil amendment increased soil Pwep independent of soil sterilization and P source treatments. Soil resin P (Pres) and microbial P (Pmic), which represents stocks of potentially plant available P, were also affected from P. aculeatum inoculation. Increased soil Pres from TCP and sewage sludge ash was observed with P. aculeatum independent of soil type. On the other hand soil Pmic was higher after P. aculeatum inoculation only in semi-sterile soil. Population density of P. aculeatum measured with qPCR was maintained or increased in non-sterile and semi-sterile soil, respectively, compared to the original inoculum load of P. aculeatum. In conclusion, our results underline the importance of C source addition for P. aculeatum if used as a biofertilizer. For this, cellulose seems to be a promising option promoting P. aculeatum growth and P solubilization also in non-sterilized soil.

A Verticillium longisporum pleiotropic drug transporter determines tolerance to the plant host ß-pinene monoterpene.

Terpenes constitute a major part of secondary metabolites secreted by plants in the rhizosphere. However, their specific functions in fungal-plant interactions have not been investigated thoroughly. In this study we investigated the role of monoterpenes in interactions between oilseed rape (Brassica napus) and the soilborne pathogen Verticillium longisporum. We identified seven monoterpenes produced by B. napus, and production of α-pinene, ß-pinene, 3-carene, and camphene was significantly increased upon fungal infection. Among them, ß-pinene was chosen for further analysis. Transcriptome analysis of V. longisporum on exposure to ß-pinene resulted in identification of two highly expressed pleotropic drug transporters paralog genes named VlAbcG1a and VlAbcG1b. Overexpression of VlAbcG1a in Saccharomyces cerevisiae increased tolerance to ß-pinene, while deletion of the VlAbcG1a homologous gene in Verticillium dahliae resulted in mutants with increased sensitivity to certain monoterpenes. Furthermore, the VlAbcG1a overexpression   strain displayed an increased tolerance to ß-pinene and increased virulence in tomato plants. Data from this study give new insights into the roles of terpenes in plant-fungal pathogen interactions and the mechanisms fungi deploy to cope with the toxicity of these secondary metabolites.

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.

Spray-Induced Gene Silencing as a Potential Tool to Control Potato Late Blight Disease.

Phytophthora infestans causes late blight disease on potato and tomato and is currently controlled by resistant cultivars or intensive fungicide spraying. Here, we investigated an alternative means for late blight control by spraying potato leaves with double-stranded RNAs (dsRNA) that target the P. infestans genes essential for infection. First, we showed that the sporangia of P. infestans expressing green fluorescent protein (GFP) can take up in vitro synthesized dsRNAs homologous to GFP directly from their surroundings, including leaves, which led to the reduced relative expression of GFP. We further demonstrate the potential of spray-induced gene silencing (SIGS) in controlling potato late blight disease by targeting developmentally important genes in P. infestans such as guanine-nucleotide binding protein ß-subunit (PiGPB1), haustorial membrane protein (PiHmp1), cutinase (PiCut3), and endo-1,3(4)-ß-glucanase (PiEndo3). Our results demonstrate that SIGS can potentially be used to mitigate potato late blight; however, the degree of disease control is dependent on the selection of the target genes.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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.

Fungal X-Intrinsic Protein Aquaporin from Trichoderma atroviride: Structural and Functional Considerations.

The major intrinsic protein (MIP) superfamily is a key part of the fungal transmembrane transport network. It facilitates the transport of water and low molecular weight solutes across biomembranes. The fungal uncharacterized X-Intrinsic Protein (XIP) subfamily includes the full protein diversity of MIP. Their biological functions still remain fully hypothetical. The aim of this study is still to deepen the diversity and the structure of the XIP subfamily in light of the MIP counterparts-the aquaporins (AQPs) and aquaglyceroporins (AQGPs)-and to describe for the first time their function in the development, biomass accumulation, and mycoparasitic aptitudes of the fungal bioagent Trichoderma atroviride. The fungus-XIP clade, with one member (TriatXIP), is one of the three clades of MIPs that make up the diversity of T. atroviride MIPs, along with the AQPs (three members) and the AQGPs (three members). TriatXIP resembles those of strict aquaporins, predicting water diffusion and possibly other small polar solutes due to particularly wider ar/R constriction with a Lysine substitution at the LE2 position. The XIP loss of function in ∆TriatXIP mutants slightly delays biomass accumulation but does not impact mycoparasitic activities. ∆TriatMIP forms colonies similar to wild type; however, the hyphae are slightly thinner and colonies produce rare chlamydospores in PDA and specific media, most of which are relatively small and exhibit abnormal morphologies. To better understand the molecular causes of these deviant phenotypes, a wide-metabolic survey of the ∆TriatXIPs demonstrates that the delayed growth kinetic, correlated to a decrease in respiration rate, is caused by perturbations in the pentose phosphate pathway. Furthermore, the null expression of the XIP gene strongly impacts the expression of four expressed MIP-encoding genes of T. atroviride, a plausible compensating effect which safeguards the physiological integrity and life cycle of the fungus. This paper offers an overview of the fungal XIP family in the biocontrol agent T. atroviride which will be useful for further functional analysis of this particular MIP subfamily in vegetative growth and the environmental stress response in fungi. Ultimately, these findings have implications for the ecophysiology of Trichoderma spp. in natural, agronomic, and industrial systems.

Functional characterization of the AGL1 aegerolysin in the mycoparasitic fungus Trichoderma atroviride reveals a role in conidiation and antagonism.

Aegerolysins are small secreted pore-forming proteins that are found in both prokaryotes and eukaryotes. The role of aegerolysins in sporulation, fruit body formation, and in lysis of cellular membrane is suggested in fungi. The aim of the present study was to characterize the biological function of the aegerolysin gene agl1 in the mycoparasitic fungus Trichoderma atroviride, used for biological control of plant diseases. Gene expression analysis showed higher expression of agl1 during conidiation and during growth in medium supplemented with cell wall material from the plant pathogenic fungus Rhizoctonia solani as the sole carbon source. Expression of agl1 was supressed under iron-limiting condition, while agl1 transcript was not detected during T. atroviride interactions with the prey fungi Botrytis cinerea or R. solani. Phenotypic analysis of agl1 deletion strains (Δagl1) showed reduced conidiation compared to T. atroviride wild type, thus suggesting the involvement of AGL1 in conidiation. Furthermore, the Δagl1 strains display reduced antagonism towards B. cinerea and R. solani based on a secretion assay, although no difference was detected during direct interactions. These data demonstrate the role of AGL1 in conidiation and antagonism in the mycoparasitic fungus T. atroviride.

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.

Combining transcriptomics and genetic linkage based information to identify candidate genes associated with Heterobasidion-resistance in Norway spruce.

The Heterobasidion annosum s.l species complex comprises the most damaging forest pathogens to Norway spruce. We revisited previously identified Quantitative Trait Loci (QTLs) related to Heterobasidion-resistance in Norway spruce to identify candidate genes associated with these QTLs. We identified 329 candidate genes associated with the resistance QTLs using a gene-based composite map for Pinaceae. To evaluate the transcriptional responses of these candidate genes to H. parviporum, we inoculated Norway spruce plants and sequenced the transcriptome of the interaction at 3 and 7 days post inoculation. Out of 298 expressed candidate genes 124 were differentially expressed between inoculation and wounding control treatment. Interestingly, PaNAC04 and two of its paralogs in the subgroup III-3 of the NAC family transcription factors were found to be associated with one of the QTLs and was also highly induced in response to H. parviporum. These genes are possibly involved in the regulation of biosynthesis of flavonoid compounds. Furthermore, several of the differentially expressed candidate genes were associated with the phenylpropanoid pathway including a phenylalanine ammonia-lyase, a cinnamoyl-CoA reductase, a caffeoyl-CoA O-methyltransferase and a PgMYB11-like transcription factor gene. Combining transcriptome and genetic linkage analyses can help identifying candidate genes for functional studies and molecular breeding in non-model species.

LysM Proteins Regulate Fungal Development and Contribute to Hyphal Protection and Biocontrol Traits in Clonostachys rosea.

Lysin motif (LysM) modules are approximately 50 amino acids long and bind to peptidoglycan, chitin and its derivatives. Certain LysM proteins in plant pathogenic and entomopathogenic fungi are shown to scavenge chitin oligosaccharides and thereby dampen host defense reactions. Other LysM proteins can protect the fungal cell wall against hydrolytic enzymes. In this study, we investigated the biological function of LysM proteins in the mycoparasitic fungus Clonostachys rosea. The C. rosea genome contained three genes coding for LysM-containing proteins and gene expression analysis revealed that lysm1 and lysm2 were induced during mycoparasitic interaction with Fusarium graminearum and during colonization of wheat roots. Lysm1 was suppressed in germinating conidia, while lysm2 was induced during growth in chitin or peptidoglycan-containing medium. Deletion of lysm1 and lysm2 resulted in mutants with increased levels of conidiation and conidial germination, but reduced ability to control plant diseases caused by F. graminearum and Botrytis cinerea. The Δlysm2 strain showed a distinct, accelerated mycelial disintegration phenotype accompanied by reduced biomass production and hyphal protection against hydrolytic enzymes including chitinases, suggesting a role of LYSM2 in hyphal protection against chitinases. The Δlysm2 and Δlysm1Δlysm2 strains displayed reduced ability to colonize wheat roots, while only Δlysm1Δlysm2 failed to suppress expression of the wheat defense response genes PR1 and PR4. Based on our data, we propose a role of LYSM1 as a regulator of fungal development and of LYSM2 in cell wall protection against endogenous hydrolytic enzymes, while both are required to suppress plant defense responses. Our findings expand the understanding of the role of LysM proteins in fungal-fungal interactions and biocontrol.

Deletion of the Nonribosomal Peptide Synthetase Gene nps1 in the Fungus Clonostachys rosea Attenuates Antagonism and Biocontrol of Plant Pathogenic Fusarium and Nematodes.

Secondary metabolites produced by biological control agents may influence the outcome of their interactions with plant pathogenic microorganisms and plants. In the present study, we investigated the role of the nonribosomal peptide synthetase gene nps1 expressed by the biocontrol fungus Clonostachys rosea. A gene expression analysis showed that nps1 was induced during confrontations with the plant pathogenic fungus Botrytis cinerea. Gene deletion strains of nps1 displayed increased growth rates and conidiation. However, the nematicidal activity of culture filtrates from C. rosea Δnps1 strains was significantly weaker than that from wild-type filtrates (P ≤ 0.001); after 24 h of incubation with culture filtrates from nps1 deletion strains, only 13 to 33% of a mixed community of nematodes were dead compared with 42% of nematodes incubated with wild-type culture filtrates. The Δnps1 strains also showed reduced biocontrol efficacy during pot experiments, thus failing to protect wheat seedlings from foot rot disease caused by the plant pathogenic fungus Fusarium graminearum. Furthermore, C. rosea Δnps1 strains were not able to reduce populations of plant-parasitic nematodes in soil or in roots of wheat as efficiently as the wild-type strain. Both C. rosea wild-type and Δnps1 strains increased the dry shoot weight and shoot length of wheat by 20 and 13%, respectively. We showed that NPS1, a putative nonribosomal peptide synthetase encoded by nps1, is a biocontrol factor, presumably by producing a hitherto unknown nonribosomal peptide compound with antifungal and nematicidal properties that contributes to the biocontrol properties of C. rosea.

Out in the Cold: Identification of Genomic Regions Associated With Cold Tolerance in the Biocontrol Fungus Clonostachys rosea Through Genome-Wide Association Mapping.

There is an increasing importance for using biocontrol agents in combating plant diseases sustainably and in the long term. As large scale genomic sequencing becomes economically viable, the impact of single nucleotide polymorphisms (SNPs) on biocontrol-associated phenotypes can be easily studied across entire genomes of fungal populations. Here, we improved a previously reported genome assembly of the biocontrol fungus Clonostachys rosea strain IK726 using the PacBio sequencing platform, which resulted in a total genome size of 70.7 Mbp and 21,246 predicted genes. We further performed whole-genome re-sequencing of 52 additional C. rosea strains isolated globally using Illumina sequencing technology, in order to perform genome-wide association studies in conditions relevant for biocontrol activity. One such condition is the ability to grow at lower temperatures commonly encountered in cryic or frigid soils in temperate regions, as these will be prevalent for protecting growing crops in temperate climates. Growth rates at 10°C on potato dextrose agar of the 53 sequenced strains of C. rosea were measured and ranged between 0.066 and 0.413 mm/day. Performing a genome wide association study, a total of 1,478 SNP markers were significantly associated with the trait and located in 227 scaffolds, within or close to (< 1000 bp distance) 265 different genes. The predicted gene products included several chaperone proteins, membrane transporters, lipases, and proteins involved in chitin metabolism with possible roles in cold tolerance. The data reported in this study provides a foundation for future investigations into the genetic basis for cold tolerance in fungi, with important implications for biocontrol.