HET-E and HET-D belong to a new subfamily of WD40 proteins involved in vegetative incompatibility specificity 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.

The Bcvic1 and Bcvic2 vegetative incompatibility genes in Botrytis cinerea encode proteins with domain architectures involved in allorecognition in other filamentous fungi

Vegetative incompatibility is a fungal allorecognition system characterised by the inability of genetically distinct conspecific fungal strains to form a viable heterokaryon and is controlled by multiple polymorphic loci termed vic (vegetative incompatibility) or het (heterokaryon incompatibility). We have genetically identified and characterised the first vic locus in the economically important, plant-pathogenic, necrotrophic fungus Botrytis cinerea. A bulked segregant approach coupled with whole genome Illumina sequencing of near-isogenic lines of B. cinerea was used to map a vic locus to a 60-kb region of the genome. Within that locus, we identified two adjacent, highly polymorphic open reading frames, Bcvic1 and Bcvic2, which encode predicted proteins that contain domain architectures implicated in vegetative incompatibility in other filamentous fungi. Bcvic1 encodes a predicted protein containing a putative serine esterase domain, a NACHT family of NTPases domain, and several Ankyrin repeats. Bcvic2 encodes a putative syntaxin protein containing a SNARE domain; such proteins typically function in vesicular transport. Deletion of Bcvic1 and Bcvic2 individually had no effect on vegetative incompatibility. However, deletion of the region containing both Bcvic1 and Bcvic2 resulted in mutant lines that were severely restricted in growth and showed loss of vegetative incompatibility. Complementation of these mutants by ectopic expression restored the growth and vegetative incompatibility phenotype, indicating that Bcvic1 and Bcvic2 are controlling vegetative incompatibility at this vic locus.

Genome sequence and annotation of Periconia digitata a hopeful biocontrol agent of phytopathogenic oomycetes

The Periconia fungal genus belongs to the phylum Ascomycota, order Pleosporales, family Periconiaceae. Periconia are found in many habitats, but little is known about their ecology. Several species from this genus produce bioactive molecules. Periconia digitata extracts were shown to be deadly active against the pine wilt nematode. Furthermore, P. digitata was shown to inhibit the plant pathogenic oomycete Phytophthora parasitica. Because P. digitata has great potential as a biocontrol agent and high quality genomic resources are still lacking in the Periconiaceae family, we generated long-read genomic data for P. digitata. Using PacBio Hifi sequencing technology, we obtained a highly-contiguous genome assembled in 13 chromosomes and totaling ca. 39 Mb. In addition, we produced a reference transcriptome, based on 12 different culture conditions, and proteomic data to support the genome annotation. Besides representing a new reference genome within the Periconiaceae, this work will contribute to our better understanding of the Eukaryotic tree of life and opens new possibilities in terms of biotechnological applications.

Microbial gasdermins: More than a billion years of pyroptotic-like cell death

In the recent past, the concept of immunity has been extended to eukaryotic and prokaryotic microorganisms, like fungi and bacteria. The latest findings have drawn remarkable evolutionary parallels between metazoan and microbial defense-related genes, unveiling a growing number of shared transkingdom components of immune systems. One such component is the gasdermin family of pore-forming proteins - executioners of a highly inflammatory immune cell death program in mammals, termed pyroptosis. Pyroptotic cell death limits the spread of intracellular pathogens by eliminating infected cells and coordinates the broader inflammatory response to infection. The microbial gasdermins have similarly been implicated in defense-related cell death reactions in fungi, bacteria and archaea. Moreover, the discovery of the molecular regulators of gasdermin cytotoxicity in fungi and bacteria, has established additional evolutionary links to mammalian pyroptotic pathways. Here, we focus on the gasdermin proteins in microorganisms and their role in organismal defense and provide perspective on this remarkable case study in comparative immunology.

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.

Bacterial NLR-related proteins protect against phage

Bacteria use a wide range of immune pathways to counter phage infection. A subset of these genes shares homology with components of eukaryotic immune systems, suggesting that eukaryotes horizontally acquired certain innate immune genes from bacteria. Here, we show that proteins containing a NACHT module, the central feature of the animal nucleotide-binding domain and leucine-rich repeat containing gene family (NLRs), are found in bacteria and defend against phages. NACHT proteins are widespread in bacteria, provide immunity against both DNA and RNA phages, and display the characteristic C-terminal sensor, central NACHT, and N-terminal effector modules. Some bacterial NACHT proteins have domain architectures similar to the human NLRs that are critical components of inflammasomes. Human disease-associated NLR mutations that cause stimulus-independent activation of the inflammasome also activate bacterial NACHT proteins, supporting a shared signaling mechanism. This work establishes that NACHT module-containing proteins are ancient mediators of innate immunity across the tree of life.

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.

Genome sequences of 24 Aspergillus niger sensu stricto strains to study strain diversity, heterokaryon compatibility, and sexual reproduction

Mating-type distribution within a phylogenetic tree, heterokaryon compatibility, and subsequent diploid formation were studied in 24 Aspergillus niger sensu stricto strains. The genomes of the 24 strains were sequenced and analyzed revealing an average of 6.1 ± 2.0 variants/kb between Aspergillus niger sensu stricto strains. The genome sequences were used together with available genome data to generate a phylogenetic tree revealing 3 distinct clades within Aspergillus niger sensu stricto. The phylogenetic tree revealed that both MAT1-1 and MAT1-2 mating types were present in each of the 3 clades. The phylogenetic differences were used to select for strains to analyze heterokaryon compatibility. Conidial color markers (fwnA and brnA) and auxotrophic markers (pyrG and nicB) were introduced via CRISPR/Cas9-based genome editing in a selection of strains. Twenty-three parasexual crosses using 11 different strains were performed. Only a single parasexual cross between genetically highly similar strains resulted in a successful formation of heterokaryotic mycelium and subsequent diploid formation, indicating widespread heterokaryon incompatibility as well as multiple active heterokaryon incompatibility systems between Aspergillus niger sensu stricto strains. The 2 vegetatively compatible strains were of 2 different mating types and a stable diploid was isolated from this heterokaryon. Sclerotium formation was induced on agar media containing Triton X-100; however, the sclerotia remained sterile and no ascospores were observed. Nevertheless, this is the first report of a diploid Aspergillus niger sensu stricto strain with 2 different mating types, which offers the unique possibility to screen for conditions that might lead to ascospore formation in A. niger.

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.

Characterization of Vegetative Incompatibility in Morchella importuna and Location of the Related-Genes by Bulk Segregant Analysis

Vegetative incompatibility (VI) is a widespread phenomenon developed in Morchella importuna, a species of ascomycete fungus that is cultivated on a rapidly expanding scale in China. Understanding the genetic bases of this nonself-recognition phenomenon is beneficial for resolving some problems that are associated with the production of this highly prized edible fungus, such as crossbreeding, strain classification, and pathogen transmission. VI is genetically controlled by het genes, organized in two different systems, namely allelic and nonallelic. These het genes have been well characterized in Podospora anserina and Neurospora crassa. In this work, putative het-homologs were identified in the genome of M. importuna, but their low allelic polymorphism in different vegetative compatibility groups (VCGs) suggested that VI in this fungus might not be regulated by these het genes. The progeny derived from vegetative compatible parents became a VCG, while the single-ascospore strains from vegetative incompatible parents were divided into four VCGs, and the interaction between the inter-group strains led to the formation of two types of barrages, viz., thin dark line and raised aggregate of hyphae. The Bulk Segregant Analysis confirmed that the genes mimpvic32 and mimpvic33 were linked to VI reactions in M. importuna; nevertheless, the formation of barrages also occurred between the pairs carrying the same allele of these two genes. In sum, the VI control system in M. importuna was complicated, and there were more other allelic or non-allelic VI-related genes.

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.

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.

Conflict, Competition, and Cooperation Regulate Social Interactions in Filamentous Fungi

Social cooperation impacts the development and survival of species. In higher taxa, kin recognition occurs via visual, chemical, or tactile cues that dictate cooperative versus competitive interactions. In microbes, the outcome of cooperative versus competitive interactions is conferred by identity at allorecognition loci, so-called kind recognition. In syncytial filamentous fungi, the acquisition of multicellularity is associated with somatic cell fusion within and between colonies. However, such intraspecific cooperation entails risks, as fusion can transmit deleterious genotypes or infectious components that reduce fitness, or give rise to cheaters that can exploit communal goods without contributing to their production. Allorecognition mechanisms in syncytial fungi regulate somatic cell fusion by operating precontact during chemotropic interactions, during cell adherence, and postfusion by triggering programmed cell death reactions. Alleles at fungal allorecognition loci are highly polymorphic, fall into distinct haplogroups, and show evolutionary signatures of balancing selection, similar to allorecognition loci across the tree of life.

Integrative Activity of Mating Loci, Environmentally Responsive Genes, and Secondary Metabolism Pathways during Sexual Development of Chaetomium globosum

Fungal diversity has amazed evolutionary biologists for decades. One societally important aspect of this diversity manifests in traits that enable pathogenicity. The opportunistic pathogen Chaetomium globosum is well adapted to a high-humidity environment and produces numerous secondary metabolites that defend it from predation. Many of these chemicals can threaten human health. Understanding the phases of the C. globosum life cycle in which these products are made enables better control and even utilization of this fungus. Among its intriguing traits is that it both is self-fertile and lacks any means of propagule-based asexual reproduction. By profiling genome-wide gene expression across the process of sexual reproduction in C. globosum and comparing it to genome-wide gene expression in the model filamentous fungus N. crassa and other closely related fungi, we revealed associations among mating-type genes, sexual developmental genes, sexual incompatibility regulators, environmentally responsive genes, and secondary metabolic pathways.

The origins and maintenance of the rich fungal diversity have been longstanding issues in evolutionary biology. To investigate how differences in expression regulation contribute to divergences in development and ecology among closely related species, transcriptomes were compared between Chaetomium globosum, a homothallic pathogenic fungus thriving in highly humid ecologies, and Neurospora crassa, a heterothallic postfire saprotroph. Gene expression was quantified in perithecia at nine distinct morphological stages during nearly synchronous sexual development. Unlike N. crassa, expression of all mating loci in C. globosum was highly correlated. Key regulators of the initiation of sexual development in response to light stimuli—including orthologs of N. crassasub-1, sub-1-dependent gene NCU00309, and asl-1—showed regulatory dynamics matching between C. globosum and N. crassa. Among 24 secondary metabolism gene clusters in C. globosum, 11—including the cochliodones biosynthesis cluster—exhibited highly coordinated expression across perithecial development. C. globosum exhibited coordinately upregulated expression of histidine kinases in hyperosmotic response pathways—consistent with gene expression responses to high humidity we identified in fellow pathogen Fusarium graminearum. Bayesian networks indicated that gene interactions during sexual development have diverged in concert with the capacities both to reproduce asexually and to live a self-compatible versus self-incompatible life cycle, shifting the hierarchical roles of genes associated with conidiation and heterokaryon incompatibility in N. crassa and C. globosum. This divergence supports an evolutionary history of loss of conidiation due to unfavorable combinations of heterokaryon incompatibility in homothallic species.

Inter-strain expression of sequence-diverse HET domain genes severely inhibits growth of Aspergillus oryzae

In the Pezizomycotina (filamentous ascomycete) species, genes that encode proteins with an HET domain (Pfam: PF06985) are reportedly involved in heterokaryon incompatibility (HI) in which cell death or growth defects are induced after fusion of cells that are genetically incompatible owing to diversities in their nucleotide sequence. HET domain genes are commonly found in Pezizomycotina genomes and are functionally characterized in only a few species. Here, we compared 44 HET domain genes between an incompatible strain pair of Aspergillus oryzae RIB40 and RIB128 and performed inter-strain expression of 37 sequence-diverse genes for mimicking HI. Four HET domain genes were identified to cause severe growth inhibition in a strain- or sequence-specific manner. Furthermore, SNPs responsible for the inhibition of cell growth were identified. This study provides an important insight into the physiological significance of sequence diversity of HET domain genes and their potential functions in HI of A. oryzae.

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.

Drivers of genetic diversity in secondary metabolic gene clusters within a fungal species

Filamentous fungi produce a diverse array of secondary metabolites (SMs) critical for defense, virulence, and communication. The metabolic pathways that produce SMs are found in contiguous gene clusters in fungal genomes, an atypical arrangement for metabolic pathways in other eukaryotes. Comparative studies of filamentous fungal species have shown that SM gene clusters are often either highly divergent or uniquely present in one or a handful of species, hampering efforts to determine the genetic basis and evolutionary drivers of SM gene cluster divergence. Here, we examined SM variation in 66 cosmopolitan strains of a single species, the opportunistic human pathogen Aspergillus fumigatus. Investigation of genome-wide within-species variation revealed 5 general types of variation in SM gene clusters: nonfunctional gene polymorphisms; gene gain and loss polymorphisms; whole cluster gain and loss polymorphisms; allelic polymorphisms, in which different alleles corresponded to distinct, nonhomologous clusters; and location polymorphisms, in which a cluster was found to differ in its genomic location across strains. These polymorphisms affect the function of representative A. fumigatus SM gene clusters, such as those involved in the production of gliotoxin, fumigaclavine, and helvolic acid as well as the function of clusters with undefined products. In addition to enabling the identification of polymorphisms, the detection of which requires extensive genome-wide synteny conservation (e.g., mobile gene clusters and nonhomologous cluster alleles), our approach also implicated multiple underlying genetic drivers, including point mutations, recombination, and genomic deletion and insertion events as well as horizontal gene transfer from distant fungi. Finally, most of the variants that we uncover within A. fumigatus have been previously hypothesized to contribute to SM gene cluster diversity across entire fungal classes and phyla. We suggest that the drivers of genetic diversity operating within a fungal species shown here are sufficient to explain SM cluster macroevolutionary patterns.

All organisms produce metabolites, which are small molecules important for growth, reproduction, and other essential functions. Some organisms, including fungi, plants, and bacteria, make specialized forms of metabolites known as “secondary” metabolites that are ecologically important and improve their producers’ chances of survival and reproduction. In fungi, the genes in pathways that synthesize secondary metabolites are typically located next to each other in the genome and organized in contiguous gene clusters. These gene clusters, along with the metabolites they produce, are highly distinct, even between otherwise similar fungi, and it is often difficult to reconstruct how these differences evolved. To understand how secondary metabolic pathways evolve in fungi, we compared secondary metabolic gene clusters in 66 strains of one species of filamentous fungus, the human pathogen Aspergillus fumigatus. We show that these gene clusters vary extensively within this species, and describe the genetic processes that cause these differences. We identify 5 types of variants: single nucleotide changes, gene and gene cluster gain and loss, different gene clusters at the same genomic position, and mobile gene clusters that “jump” around the genome. These results provide a road map to the types and frequencies of genomic changes underlying the extensive diversity of fungal secondary metabolites.

A TALE-inspired computational screen for proteins that contain approximate tandem repeats

TAL (transcription activator-like) effectors (TALEs) are bacterial proteins that are secreted from bacteria to plant cells to act as transcriptional activators. TALEs and related proteins (RipTALs, BurrH, MOrTL1 and MOrTL2) contain approximate tandem repeats that differ in conserved positions that define specificity. Using PERL, we screened ~47 million protein sequences for TALE-like architecture characterized by approximate tandem repeats (between 30 and 43 amino acids in length) and sequence variability in conserved positions, without requiring sequence similarity to TALEs. Candidate proteins were scored according to their propensity for nuclear localization, secondary structure, repeat sequence complexity, as well as covariation and predicted structural proximity of variable residues. Biological context was tentatively inferred from co-occurrence of other domains and interactome predictions. Approximate repeats with TALE-like features that merit experimental characterization were found in a protein of chestnut blight fungus, a eukaryotic plant pathogen.

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.

Exploring the genomic diversity of black yeasts and relatives (Chaetothyriales, Ascomycota)

The order Chaetothyriales (Pezizomycotina, Ascomycetes) harbours obligatorily melanised fungi and includes numerous etiologic agents of chromoblastomycosis, phaeohyphomycosis and other diseases of vertebrate hosts. Diseases range from mild cutaneous to fatal cerebral or disseminated infections and affect humans and cold-blooded animals globally. In addition, Chaetothyriales comprise species with aquatic, rock-inhabiting, ant-associated, and mycoparasitic life-styles, as well as species that tolerate toxic compounds, suggesting a high degree of versatile extremotolerance. To understand their biology and divergent niche occupation, we sequenced and annotated a set of 23 genomes of main the human opportunists within the Chaetothyriales as well as related environmental species. Our analyses included fungi with diverse life-styles, namely opportunistic pathogens and closely related saprobes, to identify genomic adaptations related to pathogenesis. Furthermore, ecological preferences of Chaetothyriales were analysed, in conjuncture with the order-level phylogeny based on conserved ribosomal genes. General characteristics, phylogenomic relationships, transposable elements, sex-related genes, protein family evolution, genes related to protein degradation (MEROPS), carbohydrate-active enzymes (CAZymes), melanin synthesis and secondary metabolism were investigated and compared between species. Genome assemblies varied from 25.81 Mb (Capronia coronata) to 43.03 Mb (Cladophialophora immunda). The bantiana-clade contained the highest number of predicted genes (12 817 on average) as well as larger genomes. We found a low content of mobile elements, with DNA transposons from Tc1/Mariner superfamily being the most abundant across analysed species. Additionally, we identified a reduction of carbohydrate degrading enzymes, specifically many of the Glycosyl Hydrolase (GH) class, while most of the Pectin Lyase (PL) genes were lost in etiological agents of chromoblastomycosis and phaeohyphomycosis. An expansion was found in protein degrading peptidase enzyme families S12 (serine-type D-Ala-D-Ala carboxypeptidases) and M38 (isoaspartyl dipeptidases). Based on genomic information, a wide range of abilities of melanin biosynthesis was revealed; genes related to metabolically distinct DHN, DOPA and pyomelanin pathways were identified. The MAT (MAting Type) locus and other sex-related genes were recognized in all 23 black fungi. Members of the asexual genera Fonsecaea and Cladophialophora appear to be heterothallic with a single copy of either MAT-1-1 or MAT-1-2 in each individual. All Capronia species are homothallic as both MAT1-1 and MAT1-2 genes were found in each single genome. The genomic synteny of the MAT-locus flanking genes (SLA2-APN2-COX13) is not conserved in black fungi as is commonly observed in Eurotiomycetes, indicating a unique genomic context for MAT in those species. The heterokaryon (het) genes expansion associated with the low selective pressure at the MAT-locus suggests that a parasexual cycle may play an important role in generating diversity among those fungi.

Suppressor of fusion, a Fusarium oxysporum homolog of Ndt80, is required for nutrient-dependent regulation of anastomosis

Heterokaryon formation is an essential step in asexual recombination in Fusarium oxysporum. Filamentous fungi have an elaborate nonself recognition machinery to prevent formation and proliferation of heterokaryotic cells, called heterokaryon incompatibility (HI). In F. oxysporum the regulation of this machinery is not well understood. In Neurospora crassa, Vib-1, a putative transcription factor of the p53-like Ndt80 family of transcription factors, has been identified as global regulator of HI. In this study we investigated the role of the F. oxysporum homolog of Vib-1, called Suf, in vegetative hyphal and conidial anastomosis tube (CAT) fusion and HI. We identified a novel function for an Ndt80 homolog as a nutrient-dependent regulator of anastomosis. Strains carrying the SUF deletion mutation display a hyper-fusion phenotype during vegetative growth as well as germling development. In addition, conidial paring of incompatible SUF deletion strains led to more heterokaryon formation, which is independent of suppression of HI. Our data provides further proof for the divergence in the functions of different members Ndt80 family. We propose that Ndt80 homologs mediate responses to nutrient quality and quantity, with specific responses varying between species.

A WDR Gene Is a Conserved Member of a Chitin Synthase Gene Cluster and Influences the Cell Wall in Aspergillus nidulans

WD40 repeat (WDR) proteins are pleiotropic molecular hubs. We identify a WDR gene that is a conserved genomic neighbor of a chitin synthase gene in Ascomycetes. The WDR gene is unique to fungi and plants, and was called Fungal Plant WD (FPWD). FPWD is within a cell wall metabolism gene cluster in the Ascomycetes (Pezizomycotina) comprising chsD, a Chs activator and a GH17 glucanase. The FPWD, AN1556.2 locus was deleted in Aspergillus nidulans strain SAA.111 by gene replacement and only heterokaryon transformants were obtained. The re-annotation of Aspergilli genomes shows that AN1556.2 consists of two tightly linked separate genes, i.e., the WDR gene and a putative beta-flanking gene of unknown function. The WDR and the beta-flanking genes are conserved genomic neighbors localized within a recently identified metabolic cell wall gene cluster in genomes of Aspergilli. The heterokaryons displayed increased susceptibility to drugs affecting the cell wall, and their phenotypes, observed by optical, confocal, scanning electron and atomic force microscopy, suggest cell wall alterations. Quantitative real-time PCR shows altered expression of some cell wall-related genes. The possible implications on cell wall biosynthesis are discussed.

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.

Genome-Based Selection and Characterization of Fusarium circinatum-Specific Sequences

Fusarium circinatum is an important pathogen of pine trees and its management in the commercial forestry environment relies largely on early detection, particularly in seedling nurseries. The fact that the entire genome of this pathogen is available opens new avenues for the development of diagnostic tools for this fungus. In this study we identified open reading frames (ORFs) unique to F. circinatum and determined that they were specific to the pathogen. The ORF identification process involved bioinformatics-based screening of all the putative F. circinatum ORFs against public databases. This was followed by functional characterization of ORFs found to be unique to F. circinatum. We used PCR- and hybridization-based approaches to confirm the presence of selected unique genes in different strains of F. circinatum and their absence from other Fusarium species for which genome sequence data are not yet available. These included species that are closely related to F. circinatum as well as those that are commonly encountered in the forestry environment. Thirty-six ORFs were identified as potentially unique to F. circinatum. Nineteen of these encode proteins with known domains while the other 17 encode proteins of unknown function. The results of our PCR analyses and hybridization assays showed that three of the selected genes were present in all of the strains of F. circinatum tested and absent from the other Fusarium species screened. These data thus indicate that the selected genes are common and unique to F. circinatum. These genes thus could be good candidates for use in rapid, in-the-field diagnostic assays specific to F. circinatum. Our study further demonstrates how genome sequence information can be mined for the identification of new diagnostic markers for the detection of plant pathogens.

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.

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.

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.

Diversification of a protein kinase cascade: IME-2 is involved in nonself recognition and programmed cell death in Neurospora crassa

Kinase cascades and the modification of proteins by phosphorylation are major mechanisms for cell signaling and communication, and evolution of these signaling pathways can contribute to new developmental or environmental response pathways. The Saccharomyces cerevisiae kinase Ime2 has been well characterized for its role in meiosis. However, recent studies have revealed alternative functions for Ime2 in both S. cerevisiae and other fungi. In the filamentous fungus Neurospora crassa, the IME2 homolog (ime-2) is not required for meiosis. Here we determine that ime-2 interacts genetically with a transcription factor vib-1 during nonself recognition and programmed cell death (PCD). Mutations in vib-1 (Δvib-1) suppress PCD due to nonself recognition events; however, a Δvib-1 Δime-2 mutant restored wild-type levels of cell death. A role for ime-2 in the post-translational processing and localization of a mitochondrial matrix protein was identified, which may implicate mitochondria in N. crassa nonself recognition and PCD. Further, Δvib-1 strains do not produce extracellular proteases, but protease secretion reverted to near wild-type levels in a Δvib-1 Δime-2 strain. Mass spectrometry analysis revealed that the VIB-1 protein is phosphorylated at several sites, including a site that matches the IME-2 consensus. The genetic and biochemical data for ime-2 and vib-1 indicate that IME-2 is a negative regulator of VIB-1 and suggest parallel negative regulation by IME-2 of a cell death pathway in N. crassa that functions in concert with the VIB-1 cell death pathway. Thus, IME2 kinase function has evolved following the divergence of S. cerevisiae and N. crassa and provides insight into the evolution of kinases and their regulatory targets.

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.

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.

Diverse interactions mediate asymmetric incompatibility by the het-6 supergene complex in Neurospora crassa

Heterokaryon incompatibility (HI) in filamentous fungi is a form of nonself recognition that operates during the vegetative phase of the life cycle. One HI gene complex in Neurospora crassa, the het-6 locus, comprises two incompatibility genes, het-6 and un-24, each having two allelic variants, Oak Ridge (OR) and Panama (PA). The un-24 gene also encodes the large subunit of ribonucleotide reductase while het-6 appears to be a member of a repetitive gene family with no other known function aside from HI. These two genes are in severe linkage disequilibrium such that only un-24(OR)het-6(OR) and un-24(PA)het-6(PA) haplotypes occur in nature. In this study we unravel several genetic interactions that govern the HI functions of this gene complex. We use novel un-24(PA)het-6(OR) strains and het-6 deletion strains to demonstrate that nonallelic interactions occur between un-24 and het-6 and reveal an allelic incompatibility interaction between the OR and PA forms of un-24 that is asymmetrically enhanced by the presence of het-6(OR) or het-6(PA). We also show how two allelic forms of vib-1, a suppressor of het-c- and mat-associated incompatibility, differentially act as recessive suppressors of HI associated with nonallelic interactions between un-24(PA) and het-6(OR). In contrast, vib-1 is a dominant suppressor of HI associated with allelic differences at un-24 and a dominant partial suppressor of the un-24(OR) and het-6(PA) nonallelic interaction. The range of suppressor activities is largely explained by an interesting differential effect on het-6(OR) and het-6(PA) transcript levels by VIB-1.

The role of macroautophagy in development of filamentous fungi

Autophagy (macroautophagy) is a bulk degradative pathway by which cytoplasmic components are delivered to the vacuole for recycling. This process is conserved from yeast to human, where it is implicated in cancer and neurodegenerative diseases. During the last decade, many ATG genes involved in autophagy have been identified, initially in Saccharomyces cerevisiae. This review summarizes the knowledge on the molecular mechanisms of autophagy using yeast as model system. Although many of the core components involved in autophagy are conserved from yeast to human, there are, nevertheless, significant differences between these organisms, for example, during autophagy initiation. Autophagy also plays an essential role in filamentous fungi especially during differentiation. Remarkably, in these species autophagy may reflect features of both yeast and mammals. This is exemplified by the finding that filamentous fungi lack the S. cerevisiae clade-specific Atg31 protein, but contain Atg101, which is absent in this clade. A reappraisal of genome data further suggests that, similar to yeast and mammals, filamentous fungi probably also contain two distinct phosphatidylinositol 3-kinase complexes. This review also summarizes the state of knowledge on the role of autophagy in filamentous fungi during differentiation, such as pathogenic development, programmed cell death during heteroincompatibility, and spore formation.

Evolution and diversity of a fungal self/nonself recognition locus

BACKGROUND

Self/nonself discrimination is an essential feature for pathogen recognition and graft rejection and is a ubiquitous phenomenon in many organisms. Filamentous fungi, such as Neurospora crassa, provide a model for analyses of population genetics/evolution of self/nonself recognition loci due to their haploid nature, small genomes and excellent genetic/genomic resources. In N. crassa, nonself discrimination during vegetative growth is determined by 11 heterokaryon incompatibility (het) loci. Cell fusion between strains that differ in allelic specificity at any of these het loci triggers a rapid programmed cell death response.

METHODOLOGY/PRINCIPAL FINDINGS

In this study, we evaluated the evolution, population genetics and selective mechanisms operating at a nonself recognition complex consisting of two closely linked loci, het-c (NCU03493) and pin-c (NCU03494). The genomic position of pin-c next to het-c is unique to Neurospora/Sordaria species, and originated by gene duplication after divergence from other species within the Sordariaceae. The het-c pin-c alleles in N. crassa are in severe linkage disequilibrium and consist of three haplotypes, het-c1/pin-c1, het-c2/pin-c2 and het-c3/pin-c3, which are equally frequent in population samples and exhibit trans-species polymorphisms. The absence of recombinant haplotypes is correlated with divergence of the het-c/pin-c intergenic sequence. Tests for positive and balancing selection at het-c and pin-c support the conclusion that both of these loci are under non-neutral balancing selection; other regions of both genes appear to be under positive selection. Our data show that the het-c2/pin-c2 haplotype emerged by a recombination event between the het-c1/pin-c1 and het-c3/pin-c3 approximately 3-12 million years ago.

CONCLUSIONS/SIGNIFICANCE

These results support models by which loci that confer nonself discrimination form by the association of polymorphic genes with genes containing HET domains. Distinct allele classes can emerge by recombination and positive selection and are subsequently maintained by balancing selection and divergence of intergenic sequence resulting in recombination blocks between haplotypes.

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.

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.

Comparative genomic analyses of the human fungal pathogens Coccidioides and their relatives

While most Ascomycetes tend to associate principally with plants, the dimorphic fungi Coccidioides immitis and Coccidioides posadasii are primary pathogens of immunocompetent mammals, including humans. Infection results from environmental exposure to Coccidiodies, which is believed to grow as a soil saprophyte in arid deserts. To investigate hypotheses about the life history and evolution of Coccidioides, the genomes of several Onygenales, including C. immitis and C. posadasii; a close, nonpathogenic relative, Uncinocarpus reesii; and a more diverged pathogenic fungus, Histoplasma capsulatum, were sequenced and compared with those of 13 more distantly related Ascomycetes. This analysis identified increases and decreases in gene family size associated with a host/substrate shift from plants to animals in the Onygenales. In addition, comparison among Onygenales genomes revealed evolutionary changes in Coccidioides that may underlie its infectious phenotype, the identification of which may facilitate improved treatment and prevention of coccidioidomycosis. Overall, the results suggest that Coccidioides species are not soil saprophytes, but that they have evolved to remain associated with their dead animal hosts in soil, and that Coccidioides metabolism genes, membrane-related proteins, and putatively antigenic compounds have evolved in response to interaction with an animal host.

Nonallelic interactions between het-c and a polymorphic locus, pin-c, are essential for nonself recognition and programmed cell death in Neurospora crassa

Nonself recognition in filamentous fungi is conferred by genetic differences at het (heterokaryon incompatibility) loci. When individuals that differ in het specificity undergo hyphal fusion, the heterokaryon undergoes a programmed cell death reaction or is highly unstable. In Neurospora crassa, three allelic specificities at the het-c locus are conferred by a highly polymorphic domain. This domain shows trans-species polymorphisms indicative of balancing selection, consistent with the role of het loci in nonself recognition. We determined that a locus closely linked to het-c, called pin-c (partner for incompatibility with het-c) was required for het-c nonself recognition and heterokaryon incompatibility (HI). The pin-c alleles in isolates that differ in het-c specificity were extremely polymorphic. Heterokaryon and transformation tests showed that nonself recognition was mediated by synergistic nonallelic interactions between het-c and pin-c, while allelic interactions at het-c increased the severity of the HI phenotype. The pin-c locus encodes a protein containing a HET domain; predicted proteins containing HET domains are frequent in filamentous ascomycete genomes. These data suggest that nonallelic interactions may be important in nonself recognition in filamentous fungi and that proteins containing a HET domain may be a key factor in these interactions.

Identification of the het-r vegetative incompatibility gene of Podospora anserina as a member of the fast evolving HNWD gene family

In fungi, vegetative incompatibility is a conspecific non-self recognition mechanism that restricts formation of viable heterokaryons when incompatible alleles of specific het loci interact. In Podospora anserina, three non-allelic incompatibility systems have been genetically defined involving interactions between het-c and het-d, het-c and het-e, het-r and het-v. het-d and het-e are paralogues belonging to the HNWD gene family that encode proteins of the STAND class. HET-D and HET-E proteins comprise an N-terminal HET effector domain, a central GTP binding site and a C-terminal WD repeat domain constituted of tandem repeats of highly conserved WD40 repeat units that define the specificity of alleles during incompatibility. The WD40 repeat units of the members of this HNWD family are undergoing concerted evolution. By combining genetic analysis and gain of function experiments, we demonstrate that an additional member of this family, HNWD2, corresponds to the het-r non-allelic incompatibility gene. As for het-d and het-e, allele specificity at the het-r locus is determined by the WD repeat domain. Natural isolates show allelic variation for het-r.

Comparative and functional characterization of intragenic tandem repeats in 10 Aspergillus genomes

Intragenic tandem repeats (ITRs) are consecutive repeats of three or more nucleotides found in coding regions. ITRs are the underlying cause of several human genetic diseases and have been associated with phenotypic variation, including pathogenesis, in several clades of the tree of life. We have examined the evolution and functional role of ITRs in 10 genomes spanning the fungal genus Aspergillus, a clade of relevance to medicine, agriculture, and industry. We identified several hundred ITRs in each of the species examined. ITR content varied extensively between species, with an average 79% of ITRs unique to a given species. For the fraction of conserved ITR regions, sequence comparisons within species and between close relatives revealed that they were highly variable. ITR-containing proteins were evolutionarily less conserved, compositionally distinct, and overrepresented for domains associated with cell-surface localization and function relative to the rest of the proteome. Furthermore, ITRs were preferentially found in proteins involved in transcription, cellular communication, and cell-type differentiation but were underrepresented in proteins involved in metabolism and energy. Importantly, although ITRs were evolutionarily labile, their functional associations appeared. To be remarkably conserved across eukaryotes. Fungal ITRs likely participate in a variety of developmental processes and cell-surface-associated functions, suggesting that their contribution to fungal lifestyle and evolution may be more general than previously assumed.

The fungus-specific HET domain mediates programmed cell death in Podospora anserina

Vegetative incompatibility is a programmed cell death reaction that occurs when fungal cells of unlike genotypes fuse. Genes defining vegetative incompatibility (het genes) are highly polymorphic, and most if not all incompatibility systems include a protein partner bearing the fungus-specific domain termed the HET domain. The nonallelic het-C/het-E incompatibility system is the best-characterized incompatibility system in Podospora anserina. Cell death is triggered by interaction of specific alleles of het-C, encoding a glycolipid transfer protein, and het-E, encoding a HET domain and a WD repeat domain involved in recognition. We show here that overexpression of the isolated HET domain from het-E results in cell death. This cell death is characterized by induction of autophagy, increased vacuolization, septation, and production of lipid droplets, which are hallmarks of cell death by incompatibility. In addition, the HET domain lethality is suppressed by the same mutations as vegetative incompatibility, but not by the inactivation of het-C. These results establish the HET domain as the mediator of cell death by incompatibility and lead to a modular conception of incompatibility systems whereby recognition is ensured by the variable regions of incompatibility proteins and cell death is triggered by the HET domain.

Genesis of a fungal non-self recognition repertoire

Conspecific allorecognition, the ability for an organism to discriminate its own cells from those of another individual of the same species, has been developed by many organisms. Allorecognition specificities are determined by highly polymorphic genes. The processes by which this extreme polymorphism is generated remain largely unknown. Fungi are able to form heterokaryons by fusion of somatic cells, and somatic non self-recognition is controlled by heterokaryon incompatibility loci (het loci). Herein, we have analyzed the evolutionary features of the het-d and het-e fungal allorecognition genes. In these het genes, allorecognition specificity is determined by a polymorphic WD-repeat domain. We found that het-d and het-e belong to a large gene family with 10 members that all share the WD-repeat domain and show that repeats of all members of the family undergo concerted evolution. It follows that repeat units are constantly exchanged both within and between members of the gene family. As a consequence, high mutation supply in the repeat domain is ensured due to the high total copy number of repeats. We then show that in each repeat four residues located at the protein/protein interaction surface of the WD-repeat domain are under positive diversifying selection. Diversification of het-d and het-e is thus ensured by high mutation supply, followed by reshuffling of the repeats and positive selection for favourable variants. We also propose that RIP, a fungal specific hypermutation process acting specifically on repeated sequences might further enhance mutation supply. The combination of these evolutionary mechanisms constitutes an original process for generating extensive polymorphism at loci that require rapid diversification.

Sexual and vegetative compatibility genes in the aspergilli

Gene flow within populations can occur by sexual and/or parasexual means. Analyses of experimental and in silico work are presented relevant to possible gene flow within the aspergilli. First, the discovery of mating-type (MAT) genes within certain species of Aspergillus is described. The implications for self-fertility, sexuality in supposedly asexual species and possible uses as phylogenetic markers are discussed. Second, the results of data mining for heterokaryon incompatibility (het) and programmed cell death (PCD) related genes in the genomes of two heterokaryon incompatible isolates of the asexual species Aspergillus niger are reported. Het-genes regulate the formation of anastomoses and heterokaryons, may protect resources and prevent the spread of infectious genetic elements. Depending on the het locus involved, hetero-allelism is not tolerated and fusion of genetically different individuals leads to growth inhibition or cell death. The high natural level of heterokaryon incompatibility in A. niger blocks parasexual analysis of the het-genes involved, but in silico experiments in the sequenced genomes allow us to identify putative het-genes. Homologous sequences to known het- and PCD-genes were compared between different sexual and asexual species including different Aspergillus species, Sordariales and the yeast Saccharomyces cerevisiae. Both het- and PCD-genes were well conserved in A. niger. However some point mutations and other small differences between the het-genes in the two A. niger isolates examined may hint to functions in heterokaryon incompatibility reactions.

Cell death by incompatibility in the fungus Podospora

Filamentous fungi are naturally able of somatic fusions. When cells of unlike genotype at specific het loci fuse, non-self recognition operates in the fusion cell and a cell death reaction termed cell death by incompatibility is triggered. In Podospora anserina cell death by incompatibility is characterized by a dramatic vacuolar enlargement, induction of autophagy and cell lysis. Autophagy contributes neither to vacuolar morphological changes nor to cell death but rather protects cells against death. Autophagy could be involved in selective elimination of pro-death signals. Vacuole collapse and cytoplasm acidification might be the cause of cell death by incompatibility.

The homologue of het-c of Neurospora crassa lacks vegetative compatibility function in Fusarium proliferatum

For two fungal strains to be vegetatively compatible and capable of forming a stable vegetative heterokaryon they must carry matching alleles at a series of loci variously termed het or vic genes. Cloned het/vic genes from Neurospora crassa and Podospora anserina have no obvious functional similarity and have various cellular functions. Our objective was to identify the homologue of the Neurospora het-c gene in Fusarium proliferatum and to determine if this gene has a vegetative compatibility function in this economically important and widely dispersed fungal pathogen. In F. proliferatum and five other closely related Fusarium species we found a few differences in the DNA sequence, but the changes were silent and did not alter the amino acid sequence of the resulting protein. Deleting the gene altered sexual fertility as the female parent, but it did not alter male fertility or existing vegetative compatibility interactions. Replacement of the allele-specific portion of the coding sequence with the sequence of an alternate allele in N. crassa did not result in a vegetative incompatibility response in transformed strains of F. proliferatum. Thus, the fphch gene in Fusarium appears unlikely to have the vegetative compatibility function associated with its homologue in N. crassa. These results suggest that the vegetative compatibility phenotype may result from convergent evolution. Thus, the genes involved in this process may need to be identified at the species level or at the level of a group of species and could prove to be attractive targets for the development of antifungal agents.