Genesis of a fungal non-self recognition repertoire

Reconstructing NOD-like receptor alleles with high internal conservation in Podospora anserina using long-read sequencing

NOD-like receptors (NLRs) are intracellular immune receptors that detect pathogen-associated cues and trigger defense mechanisms, including regulated cell death. In filamentous fungi, some NLRs mediate heterokaryon incompatibility, a self/non-self recognition process that prevents the vegetative fusion of genetically distinct individuals, reducing the risk of parasitism. The het-d and het-e NLRs in Podospora anserina are highly polymorphic incompatibility genes (het genes) whose products recognize different alleles of the het-c gene via a sensor domain composed of WD40 repeats. These repeats display unusually high sequence identity maintained by concerted evolution. However, some sites within individual repeats are hypervariable and under diversifying selection. Despite extensive genetic studies, inconsistencies in the reported WD40 domain sequence have hindered functional and evolutionary analyses. Here we demonstrate that the WD40 domain can be accurately reconstructed from long-read sequencing (Oxford Nanopore and PacBio) data, but not from Illumina-based assemblies. Functional alleles are usually formed by 11 highly conserved repeats, with different repeat combinations underlying the same phenotypic het-d and het-e incompatibility reactions. Protein structure models suggest that their WD40 domain folds into two 7-blade β-propellers composed of the highly conserved repeats, as well as three cryptic divergent repeats at the C-terminus. We additionally show that one particular het-e allele does not have an incompatibility reaction with common het-c alleles, despite being 11-repeats long. Our findings provide a robust foundation for future research into the molecular mechanisms and evolutionary dynamics of het NLRs, while also highlighting both the fragility and the flexibility of β-propellers as immune sensor domains.

What are the 100 most cited fungal genera?

The global diversity of fungi has been estimated between 2 to 11 million species, of which only about 155 000 have been named. Most fungi are invisible to the unaided eye, but they represent a major component of biodiversity on our planet, and play essential ecological roles, supporting life as we know it. Although approximately 20 000 fungal genera are presently recognised, the ecology of most remains undetermined. Despite all this diversity, the mycological community actively researches some fungal genera more commonly than others. This poses an interesting question: why have some fungal genera impacted mycology and related fields more than others? To address this issue, we conducted a bibliometric analysis to identify the top 100 most cited fungal genera. A thorough database search of the Web of Science, Google Scholar, and PubMed was performed to establish which genera are most cited. The most cited 10 genera are Saccharomyces, Candida, Aspergillus, Fusarium, Penicillium, Trichoderma, Botrytis, Pichia, Cryptococcus and Alternaria. Case studies are presented for the 100 most cited genera with general background, notes on their ecology and economic significance and important research advances. This paper provides a historic overview of scientific research of these genera and the prospect for further research. Citation: Bhunjun CS, Chen YJ, Phukhamsakda C, Boekhout T, Groenewald JZ, McKenzie EHC, Francisco EC, Frisvad JC, Groenewald M, Hurdeal VG, Luangsa-ard J, Perrone G, Visagie CM, Bai FY, Błaszkowski J, Braun U, de Souza FA, de Queiroz MB, Dutta AK, Gonkhom D, Goto BT, Guarnaccia V, Hagen F, Houbraken J, Lachance MA, Li JJ, Luo KY, Magurno F, Mongkolsamrit S, Robert V, Roy N, Tibpromma S, Wanasinghe DN, Wang DQ, Wei DP, Zhao CL, Aiphuk W, Ajayi-Oyetunde O, Arantes TD, Araujo JC, Begerow D, Bakhshi M, Barbosa RN, Behrens FH, Bensch K, Bezerra JDP, Bilański P, Bradley CA, Bubner B, Burgess TI, Buyck B, Čadež N, Cai L, Calaça FJS, Campbell LJ, Chaverri P, Chen YY, Chethana KWT, Coetzee B, Costa MM, Chen Q, Custódio FA, Dai YC, Damm U, de Azevedo Santiago ALCM, De Miccolis Angelini RM, Dijksterhuis J, Dissanayake AJ, Doilom M, Dong W, Alvarez-Duarte E, Fischer M, Gajanayake AJ, Gené J, Gomdola D, Gomes AAM, Hausner G, He MQ, Hou L, Iturrieta-González I, Jami F, Jankowiak R, Jayawardena RS, Kandemir H, Kiss L, Kobmoo N, Kowalski T, Landi L, Lin CG, Liu JK, Liu XB, Loizides M, Luangharn T, Maharachchikumbura SSN, Makhathini Mkhwanazi GJ, Manawasinghe IS, Marin-Felix Y, McTaggart AR, Moreau PA, Morozova OV, Mostert L, Osiewacz HD, Pem D, Phookamsak R, Pollastro S, Pordel A, Poyntner C, Phillips AJL, Phonemany M, Promputtha I, Rathnayaka AR, Rodrigues AM, Romanazzi G, Rothmann L, Salgado-Salazar C, Sandoval-Denis M, Saupe SJ, Scholler M, Scott P, Shivas RG, Silar P, Souza-Motta CM, Silva-Filho AGS, Spies CFJ, Stchigel AM, Sterflinger K, Summerbell RC, Svetasheva TY, Takamatsu S, Theelen B, Theodoro RC, Thines M, Thongklang N, Torres R, Turchetti B, van den Brule T, Wang XW, Wartchow F, Welti S, Wijesinghe SN, Wu F, Xu R, Yang ZL, Yilmaz N, Yurkov A, Zhao L, Zhao RL, Zhou N, Hyde KD, Crous PW (2024). What are the 100 most cited fungal genera? Studies in Mycology 108: 1-411. doi: 10.3114/sim.2024.108.01.

Plants interfere with non-self recognition of a phytopathogenic fungus via proline accumulation to facilitate mycovirus transmission

Non-self recognition is a fundamental aspect of life, serving as a crucial mechanism for mitigating proliferation of molecular parasites within fungal populations. However, studies investigating the potential interference of plants with fungal non-self recognition mechanisms are limited. Here, we demonstrate a pronounced increase in the efficiency of horizontal mycovirus transmission between vegetatively incompatible Sclerotinia sclerotiorum strains in planta as compared to in vitro. This increased efficiency is associated with elevated proline concentration in plants following S. sclerotiorum infection. This surge in proline levels attenuates the non-self recognition reaction among fungi by inhibition of cell death, thereby facilitating mycovirus transmission. Furthermore, our field experiments reveal that the combined deployment of hypovirulent S. sclerotiorum strains harboring hypovirulence-associated mycoviruses (HAVs) together with exogenous proline confers substantial protection to oilseed rape plants against virulent S. sclerotiorum. This unprecedented discovery illuminates a novel pathway by which plants can counteract S. sclerotiorum infection, leveraging the weakening of fungal non-self recognition and promotion of HAVs spread. These promising insights provide an avenue to explore for developing innovative biological control strategies aimed at mitigating fungal diseases in plants by enhancing the efficacy of horizontal HAV transmission.

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.

Exploring a diverse world of effector domains and amyloid signaling motifs in fungal NLR proteins

NLR proteins are intracellular receptors constituting a conserved component of the innate immune system of cellular organisms. In fungi, NLRs are characterized by high diversity of architectures and presence of amyloid signaling. Here, we explore the diverse world of effector and signaling domains of fungal NLRs using state-of-the-art bioinformatic methods including MMseqs2 for fast clustering, probabilistic context-free grammars for sequence analysis, and AlphaFold2 deep neural networks for structure prediction. In addition to substantially improving the overall annotation, especially in basidiomycetes, the study identifies novel domains and reveals the structural similarity of MLKL-related HeLo- and Goodbye-like domains forming the most abundant superfamily of fungal NLR effectors. Moreover, compared to previous studies, we found several times more amyloid motif instances, including novel families, and validated aggregating and prion-forming properties of the most abundant of them in vitro and in vivo. Also, through an extensive in silico search, the NLR-associated amyloid signaling was identified in basidiomycetes. The emerging picture highlights similarities and differences in the NLR architectures and amyloid signaling in ascomycetes, basidiomycetes and other branches of life.

Comparative genomic analysis reveals contraction of gene families with putative roles in pathogenesis in the fungal boxwood pathogens Calonectria henricotiae and C. pseudonaviculata

Background

Boxwood blight disease caused by Calonectria henricotiae and C. pseudonaviculata is of ecological and economic significance in cultivated and native ecosystems worldwide. Prior research has focused on understanding the population genetic and genomic diversity of C. henricotiae and C. pseudonaviculata, but gene family evolution in the context of host adaptation, plant pathogenesis, and trophic lifestyle is poorly understood. This study applied bioinformatic and phylogenetic methods to examine gene family evolution in C. henricotiae, C. pseudonaviculata and 22 related fungi in the Nectriaceae that vary in pathogenic and saprobic (apathogenic) lifestyles.

Results

A total of 19,750 gene families were identified in the 24 genomes, of which 422 were rapidly evolving. Among the six Calonectria species, C. henricotiae and C. pseudonaviculata were the only species to experience high levels of rapid contraction of pathogenesis-related gene families (89% and 78%, respectively). In contrast, saprobic species Calonectria multiphialidica and C. naviculata, two of the closest known relatives of C. henricotiae and C. pseudonaviculata, showed rapid expansion of pathogenesis-related gene families.

Conclusions

Our results provide novel insight into gene family evolution within C. henricotiae and C. pseudonaviculata and suggest gene family contraction may have contributed to limited host-range expansion of these pathogens within the plant family Buxaceae.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02035-4.

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.

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.

Identification of NLR-associated Amyloid Signaling Motifs in Bacterial Genomes

In filamentous fungi, amyloid signaling sequences allow Nod-like receptors (NLRs) to activate downstream cell-death inducing proteins with HeLo and HeLo-like (HELL) domains and amyloid RHIM and RHIM-related motifs control immune defense pathways in mammals and flies. Herein, we show bioinformatically that analogous amyloid signaling motifs exist in bacteria. These short motifs are found at the N terminus of NLRs and at the C terminus of proteins with a domain we term BELL. The corresponding NLR and BELL proteins are encoded by adjacent genes. We identify 10 families of such bacterial amyloid signaling sequences (BASS), one of which (BASS3) is homologous to RHIM and a fungal amyloid motif termed PP. BASS motifs occur nearly exclusively in bacteria forming multicellular structures (mainly in Actinobacteria and Cyanobacteria). We analyze experimentally a subset of seven of these motifs (from the most common BASS1 family and the RHIM-related BASS3 family) and find that these sequences form fibrils in vitro. Using a fungal in vivo model, we show that all tested BASS-motifs form prions and that the NLR-side motifs seed prion-formation of the corresponding BELL-side motif. We find that BASS3 motifs show partial prion cross-seeding with mammalian RHIM and fungal PP-motifs and that proline mutations on key positions of the BASS3 core motif, conserved in RHIM and PP-motifs, abolish prion formation. This work expands the paradigm of prion amyloid signaling to multicellular prokaryotes and suggests a long-term evolutionary conservation of these motifs from bacteria, to fungi and animals.

Amyloid Signaling in Filamentous Fungi and Bacteria

Amyloids are implicated in many protein misfolding diseases. Amyloid folds, however, also display a range of functional roles particularly in the microbial world. The templating ability of these folds endows them with specific properties allowing their self-propagation and protein-to-protein transmission in vivo. This property, the prion principle, is exploited by specific signaling pathways that use transmission of the amyloid fold as a way to convey information from a receptor to an effector protein. I describe here amyloid signaling pathways involving fungal nucleotide binding and oligomerization domain (NOD)-like receptors that were found to control nonself recognition and programmed cell death processes. Studies on these fungal amyloid signaling motifs stem from the characterization of the fungal [Het-s] prion protein and have led to the identification in fungi but also in multicellular bacteria of several distinct families of signaling motifs, one of which is related to RHIM [receptor-interacting protein (RIP) homotypic interaction motif], an amyloid motif regulating mammalian necroptosis.

WD40 Repeat Proteins: Signalling Scaffold with Diverse Functions

The WD40 domain is one of the most abundant and interacting domains in the eukaryotic genome. In proteins the WD domain folds into a β-propeller structure, providing a platform for the interaction and assembly of several proteins into a signalosome. WD40 repeats containing proteins, in lower eukaryotes, are mainly involved in growth, cell cycle, development and virulence, while in higher organisms, they play an important role in diverse cellular functions like signal transduction, cell cycle control, intracellular transport, chromatin remodelling, cytoskeletal organization, apoptosis, development, transcriptional regulation, immune responses. To play the regulatory role in various processes, they act as a scaffold for protein-protein or protein-DNA interaction. So far, no WD40 domain has been identified with intrinsic enzymatic activity. Several WD40 domain-containing proteins have been recently characterized in prokaryotes as well. The review summarizes the vast array of functions performed by different WD40 domain containing proteins, their domain organization and functional conservation during the course of evolution.

Divergence of Funneliformis mosseae populations over 20 years of laboratory cultivation, as revealed by vegetative incompatibility and molecular analysis

Arbuscular mycorrhizal fungi (AMF) are widespread, important plant symbionts. They absorb and translocate mineral nutrients from the soil to host plants through an extensive extraradical mycelium, consisting of indefinitely large networks of nonseptate, multinucleated hyphae which may be interconnected by hyphal fusions (anastomoses). This work investigated whether different lineages of the same isolate may lose the ability to establish successful anastomoses, becoming vegetatively incompatible, when grown separately. The occurrence of hyphal incompatibility among five lineages of Funneliformis mosseae, originated from the same ancestor isolate and grown in vivo for more than 20 years in different European locations, was assessed by systematic detection of anastomosis frequency and cytological studies. Anastomosis frequencies ranged from 60 to 80% within the same lineage and from 17 to 44% among different lineages. The consistent detection of protoplasm continuity and nuclei in perfect fusions showed active protoplasm flow both within and between lineages. In pairings between different lineages, post-fusion incompatible reactions occurred in 6-48% of hyphal contacts and pre-fusion incompatibility in 2-17%. Molecular fingerprinting profiles showed genetic divergence among lineages, with overall Jaccard similarity indices ranging from 0.85 to 0.95. Here, phenotypic divergence among the five F. mosseae lineages was demonstrated by the reduction of their ability to form anastomosis and the detection of high levels of vegetative incompatibility. Our data suggest that potential genetic divergence may occur in AMF over only 20 years and represent the basis for detailed studies on the relationship between genes regulating anastomosis formation and hyphal compatibility in AMF.

Balancing selection at nonself recognition loci in the chestnut blight fungus, Cryphonectria parasitica, demonstrated by trans-species polymorphisms, positive selection, and even allele frequencies

Balancing selection has been inferred in diverse organisms for nonself recognition genes, including those involved in immunity, mating compatibility, and vegetative incompatibility. Although selective forces maintaining polymorphisms are known for genes involved in immunity and mating, mechanisms of balancing selection for vegetative incompatibility genes in fungi are being debated. We hypothesized that allorecognition and its consequent inhibition of virus transmission contribute to the maintenance of polymorphisms in vegetative incompatibility loci (vic) in the chestnut blight fungus, Cryphonectria parasitica. Balancing selection was demonstrated at two loci, vic2 and vic6, by trans-species polymorphisms in C. parasitica, C. radicalis, and C. japonica and signatures of positive selection in gene sequences. In addition, more than half (31 of 54) of allele frequency estimates at six vic loci in nine field populations of C. parasitica from Asia and the eastern US were not significantly different from 0.5, as expected at equilibrium for two alleles per locus under balancing selection. At three vic loci, deviations from 0.5 were predicted based on the effects of heteroallelism on virus transmission. Twenty-five of 27 allele frequency estimates were greater than or equal to 0.5 for the allele that confers significantly stronger inhibition of virus transmission at three loci with asymmetric transmission. These results are consistent with the allorecognition hypothesis that vegetative incompatibility genes are under selection because of their role in reducing infection by viruses.

Regulated Forms of Cell Death in Fungi

Cell death occurs in all domains of life. While some cells die in an uncontrolled way due to exposure to external cues, other cells die in a regulated manner as part of a genetically encoded developmental program. Like other eukaryotic species, fungi undergo programmed cell death (PCD) in response to various triggers. For example, exposure to external stress conditions can activate PCD pathways in fungi. Calcium redistribution between the extracellular space, the cytoplasm and intracellular storage organelles appears to be pivotal for this kind of cell death. PCD is also part of the fungal life cycle, in which it occurs during sexual and asexual reproduction, aging, and as part of development associated with infection in phytopathogenic fungi. Additionally, a fungal non-self-recognition mechanism termed heterokaryon incompatibility (HI) also involves PCD. Some of the molecular players mediating PCD during HI show remarkable similarities to major constituents involved in innate immunity in metazoans and plants. In this review we discuss recent research on fungal PCD mechanisms in comparison to more characterized mechanisms in metazoans. We highlight the role of PCD in fungi in response to exogenic compounds, fungal development and non-self-recognition processes and discuss identified intracellular signaling pathways and molecules that regulate fungal PCD.

Molecular Mechanisms Regulating Cell Fusion and Heterokaryon Formation in Filamentous Fungi

For the majority of fungal species, the somatic body of an individual is a network of interconnected cells sharing a common cytoplasm and organelles. This syncytial organization contributes to an efficient distribution of resources, energy, and biochemical signals. Cell fusion is a fundamental process for fungal development, colony establishment, and habitat exploitation and can occur between hyphal cells of an individual colony or between colonies of genetically distinct individuals. One outcome of cell fusion is the establishment of a stable heterokaryon, culminating in benefits for each individual via shared resources or being of critical importance for the sexual or parasexual cycle of many fungal species. However, a second outcome of cell fusion between genetically distinct strains is formation of unstable heterokaryons and the induction of a programmed cell death reaction in the heterokaryotic cells. This reaction of nonself rejection, which is termed heterokaryon (or vegetative) incompatibility, is widespread in the fungal kingdom and acts as a defense mechanism against genome exploitation and mycoparasitism. Here, we review the currently identified molecular players involved in the process of somatic cell fusion and its regulation in filamentous fungi. Thereafter, we summarize the knowledge of the molecular determinants and mechanism of heterokaryon incompatibility and place this phenomenon in the broader context of biotropic interactions and immunity.

Whole genome sequencing and comparative genomics of closely related Fusarium Head Blight fungi: Fusarium graminearum, F. meridionale and F. asiaticum

BACKGROUND

The Fusarium graminearum species complex is composed of many distinct fungal species that cause several diseases in economically important crops, including Fusarium Head Blight of wheat. Despite being closely related, these species and individuals within species have distinct phenotypic differences in toxin production and pathogenicity, with some isolates reported as non-pathogenic on certain hosts. In this report, we compare genomes and gene content of six new isolates from the species complex, including the first available genomes of F. asiaticum and F. meridionale, with four other genomes reported in previous studies.

RESULTS

A comparison of genome structure and gene content revealed a 93-99% overlap across all ten genomes. We identified more than 700 k base pairs (kb) of single nucleotide polymorphisms (SNPs), insertions, and deletions (indels) within common regions of the genome, which validated the species and genetic populations reported within species. We constructed a non-redundant pan gene list containing 15,297 genes from the ten genomes and among them 1827 genes or 12% were absent in at least one genome. These genes were co-localized in telomeric regions and select regions within chromosomes with a corresponding increase in SNPs and indels. Many are also predicted to encode for proteins involved in secondary metabolism and other functions associated with disease. Genes that were common between isolates contained high levels of nucleotide variation and may be pseudogenes, allelic, or under diversifying selection.

CONCLUSIONS

The genomic resources we have contributed will be useful for the identification of genes that contribute to the phenotypic variation and niche specialization that have been reported among members of the F. graminearum species complex.

Evolution of Protein Domain Repeats in Metazoa

Repeats are ubiquitous elements of proteins and they play important roles for cellular function and during evolution. Repeats are, however, also notoriously difficult to capture computationally and large scale studies so far had difficulties in linking genetic causes, structural properties and evolutionary trajectories of protein repeats. Here we apply recently developed methods for repeat detection and analysis to a large dataset comprising over hundred metazoan genomes. We find that repeats in larger protein families experience generally very few insertions or deletions (indels) of repeat units but there is also a significant fraction of noteworthy volatile outliers with very high indel rates. Analysis of structural data indicates that repeats with an open structure and independently folding units are more volatile and more likely to be intrinsically disordered. Such disordered repeats are also significantly enriched in sites with a high functional potential such as linear motifs. Furthermore, the most volatile repeats have a high sequence similarity between their units. Since many volatile repeats also show signs of recombination, we conclude they are often shaped by concerted evolution. Intriguingly, many of these conserved yet volatile repeats are involved in host-pathogen interactions where they might foster fast but subtle adaptation in biological arms races. KEY WORDS: protein evolution, domain rearrangements, protein repeats, concerted evolution.

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.

Insights into Adaptations to a Near-Obligate Nematode Endoparasitic Lifestyle from the Finished Genome of Drechmeria coniospora

Nematophagous fungi employ three distinct predatory strategies: nematode trapping, parasitism of females and eggs, and endoparasitism. While endoparasites play key roles in controlling nematode populations in nature, their application for integrated pest management is hindered by the limited understanding of their biology. We present a comparative analysis of a high quality finished genome assembly of Drechmeria coniospora, a model endoparasitic nematophagous fungus, integrated with a transcriptomic study. Adaptation of D. coniospora to its almost completely obligate endoparasitic lifestyle led to the simplification of many orthologous gene families involved in the saprophytic trophic mode, while maintaining orthologs of most known fungal pathogen-host interaction proteins, stress response circuits and putative effectors of the small secreted protein type. The need to adhere to and penetrate the host cuticle led to a selective radiation of surface proteins and hydrolytic enzymes. Although the endoparasite has a simplified secondary metabolome, it produces a novel peptaibiotic family that shows antibacterial, antifungal and nematicidal activities. Our analyses emphasize the basic malleability of the D. coniospora genome: loss of genes advantageous for the saprophytic lifestyle; modulation of elements that its cohort species utilize for entomopathogenesis; and expansion of protein families necessary for the nematode endoparasitic lifestyle.

Whole genome resequencing of Botrytis cinerea isolates identifies high levels of standing diversity

How standing genetic variation within a pathogen contributes to diversity in host/pathogen interactions is poorly understood, partly because most studied pathogens are host-specific, clonally reproducing organisms which complicates genetic analysis. In contrast, Botrytis cinerea is a sexually reproducing, true haploid ascomycete that can infect a wide range of diverse plant hosts. While previous work had shown significant genomic variation between two isolates, we proceeded to assess the level and frequency of standing variation in a population of B. cinerea. To begin measuring standing genetic variation in B. cinerea, we re-sequenced the genomes of 13 different isolates and aligned them to the previously sequenced T4 reference genome. In addition one of these isolates was resequenced from four independently repeated cultures. A high level of genetic diversity was found within the 13 isolates. Within this variation, we could identify clusters of genes with major effect polymorphisms, i.e., polymorphisms that lead to a predicted functional knockout, that surrounded genes involved in controlling vegetative incompatibility. The genotype at these loci was able to partially predict the interaction of these isolates in vegetative fusion assays showing that these loci control vegetative incompatibility. This suggests that the vegetative incompatibility loci within B. cinerea are associated with regions of increased genetic diversity. The genome re-sequencing of four clones from the one isolate (Grape) that had been independently propagated over 10 years showed no detectable spontaneous mutation. This suggests that B. cinerea does not display an elevated spontaneous mutation rate. Future work will allow us to test if, and how, this diversity may be contributing to the pathogen's broad host range.

Sphingolipid transfer proteins defined by the GLTP-fold

Glycolipid transfer proteins (GLTPs) originally were identified as small (~24 kDa), soluble, amphitropic proteins that specifically accelerate the intermembrane transfer of glycolipids. GLTPs and related homologs now are known to adopt a unique, helically dominated, two-layer 'sandwich' architecture defined as the GLTP-fold that provides the structural underpinning for the eukaryotic GLTP superfamily. Recent advances now provide exquisite insights into structural features responsible for lipid headgroup selectivity as well as the adaptability of the hydrophobic compartment for accommodating hydrocarbon chains of differing length and unsaturation. A new understanding of the structural versatility and evolutionary premium placed on the GLTP motif has emerged. Human GLTP-motifs have evolved to function not only as glucosylceramide binding/transferring domains for phosphoinositol 4-phosphate adaptor protein-2 during glycosphingolipid biosynthesis but also as selective binding/transfer proteins for ceramide-1-phosphate. The latter, known as ceramide-1-phosphate transfer protein, recently has been shown to form GLTP-fold while critically regulating Group-IV cytoplasmic phospholipase A2 activity and pro-inflammatory eicosanoid production.

Identification of Allorecognition Loci in Neurospora crassa by Genomics and Evolutionary Approaches

Understanding the genetic and molecular bases of the ability to distinguish self from nonself (allorecognition) and mechanisms underlying evolution of allorecognition systems is an important endeavor for understanding cases where it becomes dysfunctional, such as in autoimmune disorders. In filamentous fungi, allorecognition can result in vegetative or heterokaryon incompatibility, which is a type of programmed cell death that occurs following fusion of genetically different cells. Allorecognition is genetically controlled by het loci, with coexpression of any combination of incompatible alleles triggering vegetative incompatibility. Herein, we identified, characterized, and inferred the evolutionary history of candidate het loci in the filamentous fungus Neurospora crassa. As characterized het loci encode proteins carrying an HET domain, we annotated HET domain genes in 25 isolates from a natural population along with the N. crassa reference genome using resequencing data. Because allorecognition systems can be affected by frequency-dependent selection favoring rare alleles (i.e., balancing selection), we mined resequencing data for HET domain loci whose alleles displayed elevated levels of variability, excess of intermediate frequency alleles, and deep gene genealogies. From these analyses, 34 HET domain loci were identified as likely to be under balancing selection. Using transformation, incompatibility assays and genetic analyses, we determined that one of these candidates functioned as a het locus (het-e). The het-e locus has three divergent allelic groups that showed signatures of positive selection, intra- and intergroup recombination, and trans-species polymorphism. Our findings represent a compelling case of balancing selection functioning on multiple alleles across multiple loci potentially involved in allorecognition.

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 expansion shapes genome architecture in the human pathogen Lichtheimia corymbifera: an evolutionary genomics analysis in the ancient terrestrial mucorales (Mucoromycotina)

Lichtheimia species are the second most important cause of mucormycosis in Europe. To provide broader insights into the molecular basis of the pathogenicity-associated traits of the basal Mucorales, we report the full genome sequence of L. corymbifera and compared it to the genome of Rhizopus oryzae, the most common cause of mucormycosis worldwide. The genome assembly encompasses 33.6 MB and 12,379 protein-coding genes. This study reveals four major differences of the L. corymbifera genome to R. oryzae: (i) the presence of an highly elevated number of gene duplications which are unlike R. oryzae not due to whole genome duplication (WGD), (ii) despite the relatively high incidence of introns, alternative splicing (AS) is not frequently observed for the generation of paralogs and in response to stress, (iii) the content of repetitive elements is strikingly low (<5%), (iv) L. corymbifera is typically haploid. Novel virulence factors were identified which may be involved in the regulation of the adaptation to iron-limitation, e.g. LCor01340.1 encoding a putative siderophore transporter and LCor00410.1 involved in the siderophore metabolism. Genes encoding the transcription factors LCor08192.1 and LCor01236.1, which are similar to GATA type regulators and to calcineurin regulated CRZ1, respectively, indicating an involvement of the calcineurin pathway in the adaption to iron limitation. Genes encoding MADS-box transcription factors are elevated up to 11 copies compared to the 1-4 copies usually found in other fungi. More findings are: (i) lower content of tRNAs, but unique codons in L. corymbifera, (ii) Over 25% of the proteins are apparently specific for L. corymbifera. (iii) L. corymbifera contains only 2/3 of the proteases (known to be essential virulence factors) in comparison to R. oryzae. On the other hand, the number of secreted proteases, however, is roughly twice as high as in R. oryzae.

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.

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.

Hypovirus molecular biology: from Koch's postulates to host self-recognition genes that restrict virus transmission

The idea that viruses can be used to control fungal diseases has been a driving force in mycovirus research since the earliest days. Viruses in the family Hypoviridae associated with reduced virulence (hypovirulence) of the chestnut blight fungus, Cryphonectria parasitica, have held a prominent place in this research. This has been due in part to the severity of the chestnut blight epidemics in North America and Europe and early reports of hypovirulence-mediated mitigation of disease in European forests and successful application for control of chestnut blight in chestnut orchards. A more recent contributing factor has been the development of a hypovirus/C. parasitica experimental system that has overcome many of the challenges associated with mycovirus research, stemming primarily from the exclusive intracellular lifestyle shared by all mycoviruses. This chapter will focus on hypovirus molecular biology with an emphasis on the development of the hypovirus/C. parasitica experimental system and its contributions to fundamental and practical advances in mycovirology and the broader understanding of virus-host interactions and fungal pathogenesis.

The transcriptional response to nonself in the fungus Podospora anserina

In fungi, heterokaryon incompatibility is a nonself recognition process occurring when filaments of different isolates of the same species fuse. Compatibility is controlled by so-called het loci and fusion of strains of unlike het genotype triggers a complex incompatibility reaction that leads to the death of the fusion cell. Herein, we analyze the transcriptional changes during the incompatibility reaction in Podospora anserina. The incompatibility response was found to be associated with a massive transcriptional reprogramming: 2231 genes were up-regulated by a factor 2 or more during incompatibility. In turn, 2441 genes were down-regulated. HET, NACHT, and HeLo domains previously found to be involved in the control of heterokaryon incompatibility were enriched in the up-regulated gene set. In addition, incompatibility was characterized by an up-regulation of proteolytic and other hydrolytic activities, of secondary metabolism clusters and toxins and effector-like proteins. The up-regulated set was found to be enriched for proteins lacking orthologs in other species and chromosomal distribution of the up-regulated genes was uneven with up-regulated genes residing preferentially in genomic islands and on chromosomes IV and V. There was a significant overlap between regulated genes during incompatibility in P. anserina and Neurospora crassa, indicating similarities in the incompatibility responses in these two species. Globally, this study illustrates that the expression changes occurring during cell fusion incompatibility in P. anserina are in several aspects reminiscent of those described in host-pathogen or symbiotic interactions in other fungal species.

MITEs in the promoters of effector genes allow prediction of novel virulence genes in Fusarium oxysporum

BACKGROUND

The plant-pathogenic fungus Fusarium oxysporum f.sp.lycopersici (Fol) has accessory, lineage-specific (LS) chromosomes that can be transferred horizontally between strains. A single LS chromosome in the Fol4287 reference strain harbors all known Fol effector genes. Transfer of this pathogenicity chromosome confers virulence to a previously non-pathogenic recipient strain. We hypothesize that expression and evolution of effector genes is influenced by their genomic context.

RESULTS

To gain a better understanding of the genomic context of the effector genes, we manually curated the annotated genes on the pathogenicity chromosome and identified and classified transposable elements. Both retro- and DNA transposons are present with no particular overrepresented class. Retrotransposons appear evenly distributed over the chromosome, while DNA transposons tend to concentrate in large chromosomal subregions. In general, genes on the pathogenicity chromosome are dispersed within the repeat landscape. Effector genes are present within subregions enriched for DNA transposons. A miniature Impala (mimp) is always present in their promoters. Although promoter deletion studies of two effector gene loci did not reveal a direct function of the mimp for gene expression, we were able to use proximity to a mimp as a criterion to identify new effector gene candidates. Through xylem sap proteomics we confirmed that several of these candidates encode proteins secreted during plant infection.

CONCLUSIONS

Effector genes in Fol reside in characteristic subregions on a pathogenicity chromosome. Their genomic context allowed us to develop a method for the successful identification of novel effector genes. Since our approach is not based on effector gene similarity, but on unique genomic features, it can easily be extended to identify effector genes in Fo strains with different host specificities.

Mycovirus transmission revisited by in situ pairings of vegetatively incompatible isolates of Cryphonectria parasitica

In disease ecology, parasite transmission is a key parameter important at both epidemiological and evolutionary scales. Mycoviruses can be transmitted both horizontally and vertically. Their horizontal transmission is strongly restricted by the host vegetative compatibility system, which controls the outcome of somatic fusion in fungi, and by the same way, may limit mycovirus transmission. However, most of current knowledge and predictive capabilities regarding these host/pathogen systems are derived from studies pairing fungal mycelia on artificial medium. An original bioassay method, using infected mycelia as well as asexual spores, had been developed to assess in situ transmission of Cryphonectria Hypovirus-1 (CHV1), a mycovirus of Cryphonectria parasitica that causes chestnut blight. For every pair of different vegetative compatibility types tested, rates of CHV1 transmission were always superior in situ than in vitro. This study supports the hypothesis that the natural ability of CHV1 to migrate within a fungal population composed of different vegetative compatible types may have been underestimated by in vitro essays. This result offers opportunities for a biological control of fungal diseases with mycoviruses.

Evidence for selection on a chordate histocompatibility locus

Allorecognition is the ability of an organism to differentiate self or close relatives from unrelated individuals. The best known applications of allorecognition are the prevention of inbreeding in hermaphroditic species (e.g., the self-incompatibility [SI] systems in plants), the vertebrate immune response to foreign antigens mediated by MHC loci, and somatic fusion, where two genetically independent individuals physically join to become a chimera. In the few model systems where the loci governing allorecognition outcomes have been identified, the corresponding proteins have exhibited exceptional polymorphism. But information about the evolution of this polymorphism outside MHC is limited. We address this subject in the ascidian Botryllus schlosseri, where allorecognition outcomes are determined by a single locus, called FuHC (Fusion/HistoCompatibility). Molecular variation in FuHC is distributed almost entirely within populations, with very little evidence for differentiation among different populations. Mutation plays a larger role than recombination in the creation of FuHC polymorphism. A selection statistic, neutrality tests, and distribution of variation within and among different populations all provide evidence for selection acting on FuHC, but are not in agreement as to whether the selection is balancing or directional.

Evolutionary genetics of the hydroid allodeterminant alr2

We surveyed genetic variation in alr2, an allodeterminant of the colonial hydroid Hydractinia symbiolongicarpus. We generated cDNA from a sample of 239 Hydractinia colonies collected at Lighthouse Point, Connecticut, and identified 473 alr2 alleles, 198 of which were unique. Rarefaction analysis suggested that the sample was near saturation. Most alleles were rare, with 86% occurring at frequencies of 1% or less. Alleles were highly variable, diverging on average by 18% of the amino acids in a predicted extracellular domain of the molecule. Analysis of 152 full-length alleles confirmed the existence of two structural types, defined by exons 4-8 of the gene. Several residues of the predicted immunoglobulin superfamily-like domains display signatures of positive selection. We also identified 77 unique alr2 pseudogene sequences from 85 colonies. Twenty-seven of these sequences matched expressed alr2 sequences from other colonies. This observation is consistent with pseudogenes contributing to alr2 diversification through sequence donation. A more limited collection of animals was made from a distant, relict population of H. symbiolongicarpus. Sixty percent of the unique sequences identified in this sample were found to match sequences from the Lighthouse Point population. The large number of alr2 alleles, their degree of divergence, the predominance of rare alleles in the population, their persistence over broad spatial and temporal scales, and the signatures of positive selection in multiple residues of the putative recognition domain paint a consistent picture of negative-frequency-dependent selection operating in this system. The genetic diversity observed at alr2 is comparable to that of the most highly polymorphic genetic systems known to date.

Genomic clustering and homology between HET-S and the NWD2 STAND protein in various fungal genomes

Background

Prions are infectious proteins propagating as self-perpetuating amyloid polymers. The [Het-s] prion of Podospora anserina is involved in a cell death process associated with non-self recognition. The prion forming domain (PFD) of HET-s adopts a β-solenoid amyloid structure characterized by the two fold repetition of an elementary triangular motif. [Het-s] induces cell death when interacting with HET-S, an allelic variant of HET-s. When templated by [Het-s], HET-S undergoes a trans-conformation, relocates to the cell membrane and induces toxicity.

Methodology/Principal Findings

Here, comparing HET-s homologs from different species, we devise a consensus for the HET-s elementary triangular motif. We use this motif to screen genomic databases and find a match to the N-terminus of NWD2, a STAND protein, encoded by the gene immediately adjacent to het-S. STAND proteins are signal transducing ATPases which undergo ligand-induced oligomerisation. Homology modelling predicts that the NWD2 N-terminal region adopts a HET-s-like fold. We propose that upon NWD2 oligomerisation, these N-terminal extensions adopt the β-solenoid fold and template HET-S to adopt the amyloid fold and trigger toxicity. We extend this model to a putative prion, the σ infectious element in Nectria haematococca, because the s locus controlling propagation of σ also encodes a STAND protein and displays analogous features. Comparative genomic analyses indicate evolutionary conservation of these STAND/prion-like gene pairs, identify a number of novel prion candidates and define, in addition to the HET-s PFD motif, two distinct, novel putative PFD-like motifs.

Conclusions/Significance

We suggest the existence, in the fungal kingdom, of a widespread and evolutionarily conserved mode of signal transduction based on the transmission of an amyloid-fold from a NOD-like STAND receptor protein to an effector protein.

Creation and maintenance of variation in allorecognition Loci: molecular analysis in various model systems

Allorecognition is the ability of an organism to differentiate self or close relatives from unrelated conspecifics. Effective allorecognition systems are critical to the survival of organisms; they prevent inbreeding and facilitate fusions between close relatives. Where the loci governing allorecognition outcomes have been identified, the corresponding proteins often exhibit exceptional polymorphism. Two important questions about this polymorphism remain unresolved: how is it created, and how is it maintained. Because the genetic bases of several allorecognition systems have now been identified, including alr1 and alr2 in Hydractinia, fusion histocompatibility in Botryllus, the het (vic) loci in fungi, tgrB1 and tgrC1 in Dictyostelium, and self-incompatibility (SI) loci in several plant families, we are now poised to achieve a clearer understanding of how these loci evolve. In this review, we summarize what is currently known about the evolution of allorecognition loci, highlight open questions, and suggest future directions.

Self/nonself recognition in Tuber melanosporum is not mediated by a heterokaryon incompatibility system

Vegetative incompatibility is a widespread phenomenon in filamentous ascomycetes, which limits formation of viable heterokaryons. Whether this phenomenon plays a role in maintaining the homokaryotic state of the hyphae during the vegetative growth of Tuber spp. Gene expression, polymorphism analysis as well as targeted in vitro experiments allowed us to test whether a heterokaryon incompatibility (HI) system operates in Tuber melanosporum. HI is controlled by different genetic systems, often involving HET domain genes and their partners whose interaction can trigger a cell death reaction. Putative homologues to HI-related genes previously characterized in Neurospora crassa and Podospora anserina were identified in the T. melanosporum genome. However, only two HET domain genes were found. In many other ascomycetes HET domains have been found within different genes including some members of the NWD (NACHT and WD-repeat associated domains) gene family of P. anserina. More than 50 NWD homologues were found in T. melanosporum but none of these contain a HET domain. All these T. melanosporum paralogs showed a conserved gene organization similar to the microexon genes only recently characterized in Schistosoma mansoni. Expression data of the annotated HI-like genes along with low allelic polymorphism suggest that they have cellular functions unrelated to HI. Moreover, morphological analyses did not provide evidence for HI reactions between pairs of genetically different T. melanosporum strains. Thus, the maintenance of the genetic integrity during the vegetative growth of this species likely depends on mechanisms that act before hyphal fusion.

Molecular characterization of vegetative incompatibility genes that restrict hypovirus transmission in the chestnut blight fungus Cryphonectria parasitica

Genetic nonself recognition systems such as vegetative incompatibility operate in many filamentous fungi to regulate hyphal fusion between genetically dissimilar individuals and to restrict the spread of virulence-attenuating mycoviruses that have potential for biological control of pathogenic fungi. We report here the use of a comparative genomics approach to identify seven candidate polymorphic genes associated with four vegetative incompatibility (vic) loci of the chestnut blight fungus Cryphonectria parasitica. Disruption of candidate alleles in one of two strains that were heteroallelic at vic2, vic6, or vic7 resulted in enhanced virus transmission, but did not prevent barrage formation associated with mycelial incompatibility. Detailed characterization of the vic6 locus revealed the involvement of nonallelic interactions between two tightly linked genes in barrage formation, heterokaryon formation, and asymmetric, gene-specific influences on virus transmission. The combined results establish molecular identities of genes associated with four C. parasitica vic loci and provide insights into how these recognition factors interact to trigger incompatibility and restrict virus transmission.

Endophytic life strategies decoded by genome and transcriptome analyses of the mutualistic root symbiont Piriformospora indica

Recent sequencing projects have provided deep insight into fungal lifestyle-associated genomic adaptations. Here we report on the 25 Mb genome of the mutualistic root symbiont Piriformospora indica (Sebacinales, Basidiomycota) and provide a global characterization of fungal transcriptional responses associated with the colonization of living and dead barley roots. Extensive comparative analysis of the P. indica genome with other Basidiomycota and Ascomycota fungi that have diverse lifestyle strategies identified features typically associated with both, biotrophism and saprotrophism. The tightly controlled expression of the lifestyle-associated gene sets during the onset of the symbiosis, revealed by microarray analysis, argues for a biphasic root colonization strategy of P. indica. This is supported by a cytological study that shows an early biotrophic growth followed by a cell death-associated phase. About 10% of the fungal genes induced during the biotrophic colonization encoded putative small secreted proteins (SSP), including several lectin-like proteins and members of a P. indica-specific gene family (DELD) with a conserved novel seven-amino acids motif at the C-terminus. Similar to effectors found in other filamentous organisms, the occurrence of the DELDs correlated with the presence of transposable elements in gene-poor repeat-rich regions of the genome. This is the first in depth genomic study describing a mutualistic symbiont with a biphasic lifestyle. Our findings provide a significant advance in understanding development of biotrophic plant symbionts and suggest a series of incremental shifts along the continuum from saprotrophy towards biotrophy in the evolution of mycorrhizal association from decomposer fungi.

Gene expression associated with vegetative incompatibility in Amylostereum areolatum

In filamentous fungi, vegetative compatibility among individuals of the same species is determined by the genes encoded at the heterokaryon incompatibility (het) loci. The hyphae of genetically similar individuals that share the same allelic specificities at their het loci are able to fuse and intermingle, while different allelic specificities at the het loci result in cell death of the interacting hyphae. In this study, suppression subtractive hybridization (SSH) followed by pyrosequencing and quantitative reverse transcription PCR were used to identify genes that are selectively expressed when vegetatively incompatible individuals of Amylostereum areolatum interact. The SSH library contained genes associated with various cellular processes, including cell-cell adhesion, stress and defence responses, as well as cell death. Some of the transcripts encoded proteins that were previously implicated in the stress and defence responses associated with vegetative incompatibility. Other transcripts encoded proteins known to be associated with programmed cell death, but have not previously been linked with vegetative incompatibility. Results of this study have considerably increased our knowledge of the processes underlying vegetative incompatibility in Basidiomycetes in general and A. areolatum in particular.

Conformational folding and stability of the HET-C2 glycolipid transfer protein fold: does a molten globule-like state regulate activity?

The glycolipid transfer protein (GLTP) superfamily is defined by the human GLTP fold that represents a novel motif for lipid binding and transfer and for reversible interaction with membranes, i.e., peripheral amphitropic proteins. Despite limited sequence homology with human GLTP, we recently showed that HET-C2 GLTP of Podospora anserina is organized conformationally as a GLTP fold. Currently, insights into the folding stability and conformational states that regulate GLTP fold activity are almost nonexistent. To gain such insights into the disulfide-less GLTP fold, we investigated the effect of a change in pH on the fungal HET-C2 GLTP fold by taking advantage of its two tryptophans and four tyrosines (compared to three tryptophans and 10 tyrosines in human GLTP). pH-induced conformational alterations were determined by changes in (i) intrinsic tryptophan fluorescence (intensity, emission wavelength maximum, and anisotropy), (ii) circular dichroism over the near-UV and far-UV ranges, including thermal stability profiles of the derivatized molar ellipticity at 222 nm, (iii) fluorescence properties of 1-anilinonaphthalene-8-sulfonic acid, and (iv) glycolipid intermembrane transfer activity monitored by Förster resonance energy transfer. Analyses of our recently determined crystallographic structure of HET-C2 (1.9 Å) allowed identification of side chain electrostatic interactions that contribute to HET-C2 GLTP fold stability and can be altered by a change in pH. Side chain interactions include numerous salt bridges and interchain cation-π interactions, but not intramolecular disulfide bridges. Histidine residues are especially important for stabilizing the local positioning of the two tryptophan residues and the conformation of adjacent chains. Induction of a low-pH-induced, molten globule-like state inhibited glycolipid intermembrane transfer by the HET-C2 GLTP fold.

WD-repeat instability and diversification of the Podospora anserina hnwd non-self recognition gene family

Background

Genes involved in non-self recognition and host defence are typically capable of rapid diversification and exploit specialized genetic mechanism to that end. Fungi display a non-self recognition phenomenon termed heterokaryon incompatibility that operates when cells of unlike genotype fuse and leads to the cell death of the fusion cell. In the fungus Podospora anserina, three genes controlling this allorecognition process het-d, het-e and het-r are paralogs belonging to the same hnwd gene family. HNWD proteins are STAND proteins (signal transduction NTPase with multiple domains) that display a WD-repeat domain controlling recognition specificity. Based on genomic sequence analysis of different P. anserina isolates, it was established that repeat regions of all members of the gene family are extremely polymorphic and undergoing concerted evolution arguing for frequent recombination within and between family members.

Results

Herein, we directly analyzed the genetic instability and diversification of this allorecognition gene family. We have constituted a collection of 143 spontaneous mutants of the het-R (HNWD2) and het-E (hnwd5) genes with altered recognition specificities. The vast majority of the mutants present rearrangements in the repeat arrays with deletions, duplications and other modifications as well as creation of novel repeat unit variants.

Conclusions

We investigate the extreme genetic instability of these genes and provide a direct illustration of the diversification strategy of this eukaryotic allorecognition gene family.

De novo assembly of a 40 Mb eukaryotic genome from short sequence reads: Sordaria macrospora, a model organism for fungal morphogenesis

Filamentous fungi are of great importance in ecology, agriculture, medicine, and biotechnology. Thus, it is not surprising that genomes for more than 100 filamentous fungi have been sequenced, most of them by Sanger sequencing. While next-generation sequencing techniques have revolutionized genome resequencing, e.g. for strain comparisons, genetic mapping, or transcriptome and ChIP analyses, de novo assembly of eukaryotic genomes still presents significant hurdles, because of their large size and stretches of repetitive sequences. Filamentous fungi contain few repetitive regions in their 30-90 Mb genomes and thus are suitable candidates to test de novo genome assembly from short sequence reads. Here, we present a high-quality draft sequence of the Sordaria macrospora genome that was obtained by a combination of Illumina/Solexa and Roche/454 sequencing. Paired-end Solexa sequencing of genomic DNA to 85-fold coverage and an additional 10-fold coverage by single-end 454 sequencing resulted in approximately 4 Gb of DNA sequence. Reads were assembled to a 40 Mb draft version (N50 of 117 kb) with the Velvet assembler. Comparative analysis with Neurospora genomes increased the N50 to 498 kb. The S. macrospora genome contains even fewer repeat regions than its closest sequenced relative, Neurospora crassa. Comparison with genomes of other fungi showed that S. macrospora, a model organism for morphogenesis and meiosis, harbors duplications of several genes involved in self/nonself-recognition. Furthermore, S. macrospora contains more polyketide biosynthesis genes than N. crassa. Phylogenetic analyses suggest that some of these genes may have been acquired by horizontal gene transfer from a distantly related ascomycete group. Our study shows that, for typical filamentous fungi, de novo assembly of genomes from short sequence reads alone is feasible, that a mixture of Solexa and 454 sequencing substantially improves the assembly, and that the resulting data can be used for comparative studies to address basic questions of fungal biology.

Structural determination and tryptophan fluorescence of heterokaryon incompatibility C2 protein (HET-C2), a fungal glycolipid transfer protein (GLTP), provide novel insights into glycolipid specificity and membrane interaction by the GLTP fold

HET-C2 is a fungal protein that transfers glycosphingolipids between membranes and has limited sequence homology with human glycolipid transfer protein (GLTP). The human GLTP fold is unique among lipid binding/transfer proteins, defining the GLTP superfamily. Herein, GLTP fold formation by HET-C2, its glycolipid transfer specificity, and the functional role(s) of its two Trp residues have been investigated. X-ray diffraction (1.9 A) revealed a GLTP fold with all key sugar headgroup recognition residues (Asp(66), Asn(70), Lys(73), Trp(109), and His(147)) conserved and properly oriented for glycolipid binding. Far-UV CD showed secondary structure dominated by alpha-helices and a cooperative thermal unfolding transition of 49 degrees C, features consistent with a GLTP fold. Environmentally induced optical activity of Trp/Tyr/Phe (2:4:12) detected by near-UV CD was unaffected by membranes containing glycolipid but was slightly altered by membranes lacking glycolipid. Trp fluorescence was maximal at approximately 355 nm and accessible to aqueous quenchers, indicating free exposure to the aqueous milieu and consistent with surface localization of the two Trps. Interaction with membranes lacking glycolipid triggered significant decreases in Trp emission intensity but lesser than decreases induced by membranes containing glycolipid. Binding of glycolipid (confirmed by electrospray injection mass spectrometry) resulted in a blue-shifted emission wavelength maximum (approximately 6 nm) permitting determination of binding affinities. The unique positioning of Trp(208) at the HET-C2 C terminus revealed membrane-induced conformational changes that precede glycolipid uptake, whereas key differences in residues of the sugar headgroup recognition center accounted for altered glycolipid specificity and suggested evolutionary adaptation for the simpler glycosphingolipid compositions of filamentous fungi.

Fungal incompatibility: evolutionary origin in pathogen defense?

In fungi, cell fusion between genetically unlike individuals triggers a cell death reaction known as the incompatibility reaction. In Podospora anserina, the genes controlling this process belong to a gene family encoding STAND proteins with an N-terminal cell death effector domain, a central NACHT domain and a C-terminal WD-repeat domain. These incompatibility genes are extremely polymorphic, subject to positive Darwinian selection and display a remarkable genetic plasticity allowing for constant diversification of the WD-repeat domain responsible for recognition of non-self. Remarkably, the architecture of these proteins is related to pathogen-recognition receptors ensuring innate immunity in plants and animals. Here, we hypothesize that these P. anserina incompatibility genes could be components of a yet-unidentified innate immune system of fungi. As already proposed in the case of plant hybrid necrosis or graft rejection in mammals, incompatibility could be a by-product of pathogen-driven divergence in host defense genes.

Monitoring autophagy in the filamentous fungus Podospora anserina

Autophagy has been monitored in the filamentous fungus Podospora anserina using electron, light, and fluorescence microscopy. In this organism autophagy can be induced either by starvation or rapamycin treatment or by het gene incompatibility. Incompatible HET products signal a cell death reaction referred to as cell death by incompatibility. In het-R het-V strain bearing the two incompatible het-R and het-V genes, cell death is induced by a simple shift in growth temperature, as incompatibility is thermosensitive. In this strain large autophagosomes are formed as revealed by electron microscopy or using the GFP-PaATG8 marker. This strain constitutes an alternative model to study autophagy. Analysis of the three autophagy mutants, DeltaPaATG1, DeltaPaATG8, and DeltapspA, reveals that autophagy is essential for aerial hyphae and female organ differentiation and involved in spore germination. During the incompatibility reaction, autophagy might protect cells from cell death as suggested by accelerated cell death observed in autophagy mutants.