The transcriptional response to nonself in the fungus Podospora anserina

het-B allorecognition in Podospora anserina is determined by pseudo-allelic interaction of genes encoding a HET and lectin fold domain protein and a PII-like protein

Filamentous fungi display allorecognition genes that trigger regulated cell death (RCD) when strains of unlike genotype fuse. Podospora anserina is one of several model species for the study of this allorecognition process termed heterokaryon or vegetative incompatibility. Incompatibility restricts transmission of mycoviruses between isolates. In P. anserina, genetic analyses have identified nine incompatibility loci, termed het loci. Here we set out to clone the genes controlling het-B incompatibility. het-B displays two incompatible alleles, het-B1 and het-B2. We find that the het-B locus encompasses two adjacent genes, Bh and Bp that exist as highly divergent allelic variants (Bh1/Bh2 and Bp1/Bp2) in the incompatible haplotypes. Bh encodes a protein with an N-terminal HET domain, a cell death inducing domain bearing homology to Toll/interleukin-1 receptor (TIR) domains and a C-terminal domain with a predicted lectin fold. The Bp product is homologous to PII-like proteins, a family of small trimeric proteins acting as sensors of adenine nucleotides in bacteria. We show that although the het-B system appears genetically allelic, incompatibility is in fact determined by the non-allelic Bh1/Bp2 interaction while the reciprocal Bh2/Bp1 interaction plays no role in incompatibility. The highly divergent C-terminal lectin fold domain of BH determines recognition specificity. Population studies and genome analyses indicate that het-B is under balancing selection with trans-species polymorphism, highlighting the evolutionary significance of the two incompatible haplotypes. In addition to emphasizing anew the central role of TIR-like HET domains in fungal RCD, this study identifies novel players in fungal allorecognition and completes the characterization of the entire het gene set in that species.

Emergence of the fungal immune system

Investigation of fungal biology has been frequently motivated by the fact that many fungal species are important plant and animal pathogens. Such efforts have contributed significantly toward our understanding of fungal pathogenic lifestyles (virulence factors and strategies) and the interplay with host immune systems. In parallel, work on fungal allorecognition systems leading to the characterization of fungal regulated cell death determinants and pathways, has been instrumental for the emergent concept of fungal immunity. The uncovered evolutionary trans-kingdom parallels between fungal regulated cell death pathways and innate immune systems incite us to reflect further on the concept of a fungal immune system. Here, I briefly review key findings that have shaped the fungal immunity paradigm, providing a perspective on what I consider its most glaring knowledge gaps. Undertaking to fill such gaps would establish firmly the fungal immune system inside the broader field of comparative immunology.

Allorecognition genes drive reproductive isolation in Podospora anserina

Allorecognition, the capacity to discriminate self from conspecific non-self, is a ubiquitous organismal feature typically governed by genes evolving under balancing selection. Here, we show that in the fungus Podospora anserina, allorecognition loci controlling vegetative incompatibility (het genes), define two reproductively isolated groups through pleiotropic effects on sexual compatibility. These two groups emerge from the antagonistic interactions of the unlinked loci het-r (encoding a NOD-like receptor) and het-v (encoding a methyltransferase and an MLKL/HeLo domain protein). Using a combination of genetic and ecological data, supported by simulations, we provide a concrete and molecularly defined example whereby the origin and coexistence of reproductively isolated groups in sympatry is driven by pleiotropic genes under balancing selection.

Transcriptomic and Ultrastructural Analyses of Pyricularia Oryzae Treated With Fungicidal Peptaibol Analogs of Trichoderma Trichogin

Eco-friendly analogs of Trichogin GA IV, a short peptaibol produced by Trichoderma longibrachiatum, were assayed against Pyricularia oryzae, the causal agent of rice blast disease. In vitro and in vivo screenings allowed us to identify six peptides able to reduce by about 70% rice blast symptoms. One of the most active peptides was selected for further studies. Microscopy analyses highlighted that the treated fungal spores could not germinate and the fluorescein-labeled peptide localized on the spore cell wall and in the agglutinated cytoplasm. Transcriptomic analysis was carried out on P. oryzae mycelium 3 h after the peptide treatment. We identified 1,410 differentially expressed genes, two-thirds of which upregulated. Among these, we found genes involved in oxidative stress response, detoxification, autophagic cell death, cell wall biogenesis, degradation and remodeling, melanin and fatty acid biosynthesis, and ion efflux transporters. Molecular data suggest that the trichogin analogs cause cell wall and membrane damages and induce autophagic cell death. Ultrastructure observations on treated conidia and hyphae confirmed the molecular data. In conclusion, these selected peptides seem to be promising alternative molecules for developing effective bio-pesticides able to control rice blast disease.

A metabolomic study of vegetative incompatibility in Cryphonectria parasitica

Vegetative incompatibility (VI) is a form of non-self allorecognition in filamentous fungi that restricts conspecific hyphal fusion and the formation of heterokaryons. In the chestnut pathogenic fungus, Cryphonectria parasitica, VI is controlled by six vic loci and has been of particular interest because it impedes the spread of hypoviruses and thus biocontrol strategies. We use nuclear magnetic resonance and high-resolution mass spectrometry to characterize alterations in the metabolome of C. parasitica over an eight-day time course of vic3 incompatibility. Our findings support transcriptomic data that indicated remodeling of secondary metabolite profiles occurs during vic3 -associated VI. VI-associated secondary metabolites include novel forms of calbistrin, decumbenone B, a sulfoxygenated farnesyl S-cysteine analog, lysophosphatidylcholines, and an as-yet unidentified group of lipid disaccharides. The farnesyl S-cysteine analog is structurally similar to pheromones predicted to be produced during VI and is here named 'crypheromonin'. Mass features associated with C. parasitica secondary metabolites skyrin, rugulosin and cryphonectric acid were also detected but were not VI specific. Partitioning of VI-associated secondary metabolites was observed, with crypheromonins and most calbistrins accumulating in the growth medium over time, whereas lysophosphatidylcholines, lipid disaccharide-associated mass features and other calbistrin-associated mass features peaked at distinct time points in the mycelium. Secondary metabolite biosynthetic gene clusters and potential biological roles associated with the detected secondary metabolites are discussed.

Transcriptome analysis implicates secondary metabolite production, redox reactions, and programmed cell death during allorecognition in Cryphonectria parasitica

The underlying molecular mechanisms of programmed cell death associated with fungal allorecognition, a form of innate immunity, remain largely unknown. In this study, transcriptome analysis was used to infer mechanisms activated during barrage formation in vic3-incompatible strains of Cryphonectria parasitica, the chestnut blight fungus. Pronounced differential expression occurred in barraging strains of genes involved in mating pheromone (mf2-1, mf2-2), secondary metabolite production, detoxification (including oxidative stress), apoptosis-related, RNA interference, and HET-domain genes. Evidence for secondary metabolite production and reactive oxygen species (ROS) accumulation is supported through UPLC-HRMS analysis and cytological staining, respectively. Differential expression of mating-related genes and HET-domain genes was further examined by RT-qPCR of incompatible interactions involving each of the six vegetative incompatibility (vic) loci in C. parasitica and revealed distinct recognition process networks. We infer that vegetative incompatibility in C. parasitica activates defence reactions that involve secondary metabolism, resulting in increased toxicity of the extra- and intracellular environment. Accumulation of ROS (and other potential toxins) may result in detoxification failure and activation of apoptosis, sporulation, and the expression of associated pheromone genes. The incompatible reaction leaves abundant traces of a process-specific metabolome as conidiation is initiated.

A RID-like putative cytosine methyltransferase homologue controls sexual development in the fungus Podospora anserina

DNA methyltransferases are ubiquitous enzymes conserved in bacteria, plants and opisthokonta. These enzymes, which methylate cytosines, are involved in numerous biological processes, notably development. In mammals and higher plants, methylation patterns established and maintained by the cytosine DNA methyltransferases (DMTs) are essential to zygotic development. In fungi, some members of an extensively conserved fungal-specific DNA methyltransferase class are both mediators of the Repeat Induced Point mutation (RIP) genome defense system and key players of sexual reproduction. Yet, no DNA methyltransferase activity of these purified RID (RIP deficient) proteins could be detected in vitro. These observations led us to explore how RID-like DNA methyltransferase encoding genes would play a role during sexual development of fungi showing very little genomic DNA methylation, if any. To do so, we used the model ascomycete fungus Podospora anserina. We identified the PaRid gene, encoding a RID-like DNA methyltransferase and constructed knocked-out ΔPaRid defective mutants. Crosses involving P. anserina ΔPaRid mutants are sterile. Our results show that, although gametes are readily formed and fertilization occurs in a ΔPaRid background, sexual development is blocked just before the individualization of the dikaryotic cells leading to meiocytes. Complementation of ΔPaRid mutants with ectopic alleles of PaRid, including GFP-tagged, point-mutated and chimeric alleles, demonstrated that the catalytic motif of the putative PaRid methyltransferase is essential to ensure proper sexual development and that the expression of PaRid is spatially and temporally restricted. A transcriptomic analysis performed on mutant crosses revealed an overlap of the PaRid-controlled genetic network with the well-known mating-types gene developmental pathway common to an important group of fungi, the Pezizomycotina.

Population genomics of Fusarium graminearum reveals signatures of divergent evolution within a major cereal pathogen

The cereal pathogen Fusarium graminearum is the primary cause of Fusarium head blight (FHB) and a significant threat to food safety and crop production. To elucidate population structure and identify genomic targets of selection within major FHB pathogen populations in North America we sequenced the genomes of 60 diverse F. graminearum isolates. We also assembled the first pan-genome for F. graminearum to clarify population-level differences in gene content potentially contributing to pathogen diversity. Bayesian and phylogenomic analyses revealed genetic structure associated with isolates that produce the novel NX-2 mycotoxin, suggesting a North American population that has remained genetically distinct from other endemic and introduced cereal-infecting populations. Genome scans uncovered distinct signatures of selection within populations, focused in high diversity, frequently recombining regions. These patterns suggested selection for genomic divergence at the trichothecene toxin gene cluster and thirteen additional regions containing genes potentially involved in pathogen specialization. Gene content differences further distinguished populations, in that 121 genes showed population-specific patterns of conservation. Genes that differentiated populations had predicted functions related to pathogenesis, secondary metabolism and antagonistic interactions, though a subset had unique roles in temperature and light sensitivity. Our results indicated that F. graminearum populations are distinguished by dozens of genes with signatures of selection and an array of dispensable accessory genes, suggesting that FHB pathogen populations may be equipped with different traits to exploit the agroecosystem. These findings provide insights into the evolutionary processes and genomic features contributing to population divergence in plant pathogens, and highlight candidate genes for future functional studies of pathogen specialization across evolutionarily and ecologically diverse fungi.

Overlapping Podospora anserina Transcriptional Responses to Bacterial and Fungal Non Self Indicate a Multilayered Innate Immune Response

Recognition and response to non self is essential to development and survival of all organisms. It can occur between individuals of the same species or between different organisms. Fungi are established models for conspecific non self recognition in the form of vegetative incompatibility (VI), a genetically controlled process initiating a programmed cell death (PCD) leading to the rejection of a fusion cell between genetically different isolates of the same species. In Podospora anserina VI is controlled by members of the hnwd gene family encoding for proteins analogous to NOD Like Receptors (NLR) immune receptors in eukaryotes. It was hypothesized that the hnwd controlled VI reaction was derived from the fungal innate immune response. Here we analyze the P. anserina transcriptional responses to two bacterial species, Serratia fonticola to which P. anserina survives and S. marcescens to which P. anserina succumbs, and compare these to the transcriptional response induced under VI conditions. Transcriptional responses to both bacteria largely overlap, however the number of genes regulated and magnitude of regulation is more important when P. anserina survives. Transcriptional responses to bacteria also overlap with the VI reaction for both up or down regulated gene sets. Genes up regulated tend to be clustered in the genome, and display limited phylogenetic distribution. In all three responses we observed genes related to autophagy to be up-regulated. Autophagy contributes to the fungal survival in all three conditions. Genes encoding for secondary metabolites and histidine kinase signaling are also up regulated in all three conditions. Transcriptional responses also display differences. Genes involved in response to oxidative stress, or encoding small secreted proteins are essentially expressed in response to bacteria, while genes encoding NLR proteins are expressed during VI. Most functions encoded in response to bacteria favor survival of the fungus while most functions up regulated during VI would lead to cell death. These differences are discussed in the frame of a multilayered response to non self in fungi.

Identification of a Novel Nematotoxic Protein by Challenging the Model Mushroom Coprinopsis cinerea with a Fungivorous Nematode

The dung of herbivores, the natural habitat of the model mushroom Coprinopsis cinerea, is a nutrient-rich but also very competitive environment for a saprophytic fungus. We showed previously that C. cinerea expresses constitutive, tissue-specific armories against antagonists such as animal predators and bacterial competitors. In order to dissect the inducible armories against such antagonists, we sequenced the poly(A)-positive transcriptome of C. cinerea vegetative mycelium upon challenge with fungivorous and bacterivorous nematodes, Gram-negative and Gram-positive bacteria and mechanical damage. As a response to the fungivorous nematode Aphelenchus avenae, C. cinerea was found to specifically induce the transcription of several genes encoding previously characterized nematotoxic lectins. In addition, a previously not characterized gene encoding a cytoplasmic protein with several predicted Ricin B-fold domains, was found to be strongly upregulated under this condition. Functional analysis of the recombinant protein revealed a high toxicity toward the bacterivorous nematode Caenorhabditis elegans. Challenge of the mycelium with A. avenae also lead to the induction of several genes encoding putative antibacterial proteins. Some of these genes were also induced upon challenge of the mycelium with the bacteria Escherichia coli and Bacillus subtilis. These results suggest that fungi have the ability to induce specific innate defense responses similar to plants and animals.

Signal transduction by a fungal NOD-like receptor based on propagation of a prion amyloid fold

In the fungus Podospora anserina, the [Het-s] prion induces programmed cell death by activating the HET-S pore-forming protein. The HET-s β-solenoid prion fold serves as a template for converting the HET-S prion-forming domain into the same fold. This conversion, in turn, activates the HET-S pore-forming domain. The gene immediately adjacent to het-S encodes NWD2, a Nod-like receptor (NLR) with an N-terminal motif similar to the elementary repeat unit of the β-solenoid fold. NLRs are immune receptors controlling cell death and host defense processes in animals, plants and fungi. We have proposed that, analogously to [Het-s], NWD2 can activate the HET-S pore-forming protein by converting its prion-forming region into the β-solenoid fold. Here, we analyze the ability of NWD2 to induce formation of the β-solenoid prion fold. We show that artificial NWD2 variants induce formation of the [Het-s] prion, specifically in presence of their cognate ligands. The N-terminal motif is responsible for this prion induction, and mutations predicted to affect the β-solenoid fold abolish templating activity. In vitro, the N-terminal motif assembles into infectious prion amyloids that display a structure resembling the β-solenoid fold. In vivo, the assembled form of the NWD2 N-terminal region activates the HET-S pore-forming protein. This study documenting the role of the β-solenoid fold in fungal NLR function further highlights the general importance of amyloid and prion-like signaling in immunity-related cell fate pathways.

The fungus Podospora anserina uses a prion amyloid fold as a signal transduction device between a Nod-like receptor and a downstream cell death execution protein.

Although amyloids are best known as protein aggregates that are responsible for fatal neurodegenerative diseases, amyloid structures can also fulfill functional roles in cells. In particular, the controlled formation of amyloid structures appears to be involved in different signaling processes in the context of programmed cell death and host defense. The [Het-s] prion of the filamentous fungus Podospora anserina is a model system in which the 3-D structure of the prion form has been solved. The [Het-s] prion works as an activation switch for a second protein termed HET-S. HET-S is a pore-forming protein that is activated when the [Het-s] prion causes its C-terminal domain to adopt an amyloid-like fold. The protein encoded by the gene adjacent to het-S is a Nod-like receptor (NLR) called NWD2. NLRs are immune receptors that control host defense and cell death processes in plants, animals, and fungi. We show that NWD2 can template the formation of the [Het-s] prion fold in a ligand-controlled manner. NWD2 has an N-terminal motif homologous to the HET-S/s prion-forming region; we find that this region is both necessary and sufficient for its prion-inducing activity, and our functional and structural approaches reveal that the N-terminal region of NWD2 adopts a fold closely related to that of the HET-S/s prion. This study illustrates how the controlled formation of a prion amyloid fold can be used in a signaling process whereby a Nod-like receptor protein activates a downstream cell death execution domain.

Diversity and variability of NOD-like receptors in fungi

Nucleotide-binding oligomerization domain (NOD)-like receptors (NLRs) are intracellular receptors that control innate immunity and other biotic interactions in animals and plants. NLRs have been characterized in plant and animal lineages, but in fungi, this gene family has not been systematically described. There is however previous indications of the involvement of NLR-like genes in nonself recognition and programmed cell death in fungi. We have analyzed 198 fungal genomes for the presence of NLRs and have annotated a total of 5,616 NLR candidates. We describe their phylogenetic distribution, domain organization, and evolution. Fungal NLRs are characterized by a great diversity of domain organizations, suggesting frequently occurring combinatorial assortments of different effector, NOD and repeat domains. The repeat domains are of the WD, ANK, and TPR type; no LRR motifs were found. As previously documented for WD-repeat domains of fungal NLRs, TPR, and ANK repeats evolve under positive selection and show highly conserved repeats and repeat length polymorphism, suggesting the possibility of concerted evolution of these repeats. We identify novel effector domains not previously found associated with NLRs, whereas others are related to effector domains of plant or animals NLRs. In particular, we show that the HET domain found in fungal NLRs may be related to Toll/interleukin-1 receptor domains found in animal and plant immune receptors. This description of fungal NLR repertoires reveals both similarities and differences with plant and animals NLR collections, highlights the importance of domain reassortment and repeat evolution and provides a novel entry point to explore the evolution of NLRs in eukaryotes.

Gene expression associated with intersterility in Heterobasidion

Intersterility (IS) is thought to prevent mating compatibility between homokaryons that belong to different species. Although IS in Heterobasidion is regulated by the genes located at the IS loci, it is not yet known how the IS genes influence sexual compatibility and heterokaryon formation. To increase our understanding of the molecular events underlying IS, we studied mRNA abundance changes during IS compatible and incompatible interactions over time. The clustering of the transcripts into expression profiles, followed by the application of Gene Ontology (GO) enrichment pathway analysis of each of the clusters, allowed inference of biological processes participating in IS. These analyses identified events involved in mating and sexual development (i.e., linked with IS compatibility), which included processes associated with cell-cell adhesion and recognition, cell cycle control and signal transduction. We also identified events potentially involved in overriding mating between individuals belonging to different species (i.e., linked with IS incompatibility), which included reactive oxygen species (ROS) production, responses to stress (especially to oxidative stress), signal transduction and metabolic biosynthesis. Our findings thus enabled detection and characterization of gene expression changes associated with IS in Heterobasidion, as well as identification of important processes and pathways associated with this phenomenon. Overall, the results of this study increase current knowledge regarding the molecular mechanisms underpinning IS in Heterobasidion and allowed for the establishment of a vital baseline for further studies.

Contribution of specific residues of the β-solenoid fold to HET-s prion function, amyloid structure and stability

The [Het-s] prion of the fungus Podospora anserina represents a good model system for studying the structure-function relationship in amyloid proteins because a high resolution solid-state NMR structure of the amyloid prion form of the HET-s prion forming domain (PFD) is available. The HET-s PFD adopts a specific β-solenoid fold with two rungs of β-strands delimiting a triangular hydrophobic core. A C-terminal loop folds back onto the rigid core region and forms a more dynamic semi-hydrophobic pocket extending the hydrophobic core. Herein, an alanine scanning mutagenesis of the HET-s PFD was conducted. Different structural elements identified in the prion fold such as the triangular hydrophobic core, the salt bridges, the asparagines ladders and the C-terminal loop were altered and the effect of these mutations on prion function, fibril structure and stability was assayed. Prion activity and structure were found to be very robust; only a few key mutations were able to corrupt structure and function. While some mutations strongly destabilize the fold, many substitutions in fact increase stability of the fold. This increase in structural stability did not influence prion formation propensity in vivo. However, if an Ala replacement did alter the structure of the core or did influence the shape of the denaturation curve, the corresponding variant showed a decreased prion efficacy. It is also the finding that in addition to the structural elements of the rigid core region, the aromatic residues in the C-terminal semi-hydrophobic pocket are critical for prion propagation. Mutations in the latter region either positively or negatively affected prion formation. We thus identify a region that modulates prion formation although it is not part of the rigid cross-β core, an observation that might be relevant to other amyloid models.

Natural variation of heterokaryon incompatibility gene het-c in Podospora anserina reveals diversifying selection

In filamentous fungi, allorecognition takes the form of heterokaryon incompatibility, a cell death reaction triggered when genetically distinct hyphae fuse. Heterokaryon incompatibility is controlled by specific loci termed het-loci. In this article, we analyzed the natural variation in one such fungal allorecognition determinant, the het-c heterokaryon incompatibility locus of the filamentous ascomycete Podospora anserina. The het-c locus determines an allogenic incompatibility reaction together with two unlinked loci termed het-d and het-e. Each het-c allele is incompatible with a specific subset of the het-d and het-e alleles. We analyzed variability at the het-c locus in a population of 110 individuals, and in additional isolates from various localities. We identified a total of 11 het-c alleles, which define 7 distinct incompatibility specificity classes in combination with the known het-d and het-e alleles. We found that the het-c allorecognition gene of P. anserina is under diversifying selection. We find a highly unequal allele distribution of het-c in the population, which contrasts with the more balanced distribution of functional groups of het-c based on their allorecognition function. One explanation for the observed het-c diversity in the population is its function in allorecognition. However, alleles that are most efficient in allorecognition are rare. An alternative and not exclusive explanation for the observed diversity is that het-c is involved in pathogen recognition. In Arabidopsis thaliana, a homolog of het-c is a pathogen effector target, supporting this hypothesis. We hypothesize that the het-c diversity in P. anserina results from both its functions in pathogen-defense, and allorecognition.

New insights into mycoviruses and exploration for the biological control of crop fungal diseases

Mycoviruses are viruses that infect fungi. A growing number of novel mycoviruses have expanded our knowledge of virology, particularly in taxonomy, ecology, and evolution. Recent progress in the study of mycoviruses has comprehensively improved our understanding of the properties of mycoviruses and has strengthened our confidence to explore hypovirulence-associated mycoviruses that control crop diseases. In this review, the advantages of using hypovirulence-associated mycoviruses to control crop diseases are discussed, and, as an example, the potential for Sclerotinia sclerotiorum hypovirulence-associated DNA virus 1 (SsHADV-1) to control the stem rot of rapeseed (Brassica napus) is also introduced. Fungal vegetative incompatibility is likely to be the key factor that limits the wide utilization of mycoviruses to control crop diseases; however, there are suggested strategies for resolving this problem.

Distribution and evolution of het gene homologs in the basidiomycota

In filamentous fungi a system known as somatic incompatibility (SI) governs self/non-self recognition. SI is controlled by a regulatory signaling network involving proteins encoded at the het (heterokaryon incompatible) loci. Despite the wide occurrence of SI, the molecular identity and structure of only a small number of het genes and their products have been characterized in the model fungi Neurospora crassa and Podospora anserina. Our aim was to identify and study the distribution and evolution of putative het gene homologs in the Basidiomycota. For this purpose we used the information available for the model fungi to identify homologs of het genes in other fungi, especially the Basidiomycota. Putative het-c, het-c2 and un-24 homologs, as well as sequences containing the NACHT, HET or WD40 domains present in the het-e, het-r, het-6 and het-d genes were identified in certain members of the Ascomycota and Basidiomycota. The widespread phylogenetic distribution of certain het genes may reflect the fact that the encoded proteins are involved in fundamental cellular processes other than SI. Although homologs of het-S were previously known only from the Sordariomycetes (Ascomycota), we also identified a putative homolog of this gene in Gymnopus luxurians (Basidiomycota, class Agaricomycetes). Furthermore, with the exception of un-24, all of the putative het genes identified occurred mostly in a multi-copy fashion, some with lineage and species-specific expansions. Overall our results indicated that gene duplication followed by gene loss and/or gene family expansion, as well as multiple events of domain fusion and shuffling played an important role in the evolution of het gene homologs of Basidiomycota and other filamentous fungi.