Evidence for balancing selection operating at the het-c heterokaryon incompatibility locus in a group of filamentous fungi

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.

Candida albicans and Candida glabrata: global priority pathogens

SUMMARYA significant increase in the incidence of Candida-mediated infections has been observed in the last decade, mainly due to rising numbers of susceptible individuals. Recently, the World Health Organization published its first fungal pathogen priority list, with Candida species listed in medium, high, and critical priority categories. This review is a synthesis of information and recent advances in our understanding of two of these species-Candida albicans and Candida glabrata. Of these, C. albicans is the most common cause of candidemia around the world and is categorized as a critical priority pathogen. C. glabrata is considered a high-priority pathogen and has become an increasingly important cause of candidemia in recent years. It is now the second most common causative agent of candidemia in many geographical regions. Despite their differences and phylogenetic divergence, they are successful as pathogens and commensals of humans. Both species can cause a broad variety of infections, ranging from superficial to potentially lethal systemic infections. While they share similarities in certain infection strategies, including tissue adhesion and invasion, they differ significantly in key aspects of their biology, interaction with immune cells, host damage strategies, and metabolic adaptations. Here we provide insights on key aspects of their biology, epidemiology, commensal and pathogenic lifestyles, interactions with the immune system, and antifungal resistance.

The Narrow Footprint of Ancient Balancing Selection Revealed by Heterokaryon Incompatibility Genes in Aspergillus fumigatus

In fungi, fusion between individuals leads to localized cell death, a phenomenon termed heterokaryon incompatibility. Generally, the genes responsible for this incompatibility are observed to be under balancing selection resulting from negative frequency-dependent selection. Here, we assess this phenomenon in Aspergillus fumigatus, a human pathogenic fungus with a very low level of linkage disequilibrium as well as an extremely high crossover rate. Using complementation of auxotrophic mutations as an assay for hyphal compatibility, we screened sexual progeny for compatibility to identify genes involved in this process, called het genes. In total, 5/148 (3.4%) offspring were compatible with a parent and 166/2,142 (7.7%) sibling pairs were compatible, consistent with several segregating incompatibility loci. Genetic mapping identified five loci, four of which could be fine mapped to individual genes, of which we tested three through heterologous expression, confirming their causal relationship. Consistent with long-term balancing selection, trans-species polymorphisms were apparent across several sister species, as well as equal allele frequencies within A. fumigatus. Surprisingly, a sliding window genome-wide population-level analysis of an independent dataset did not show increased Tajima's D near these loci, in contrast to what is often found surrounding loci under balancing selection. Using available de novo assemblies, we show that these balanced polymorphisms are restricted to several hundred base pairs flanking the coding sequence. In addition to identifying the first het genes in an Aspergillus species, this work highlights the interaction of long-term balancing selection with rapid linkage disequilibrium decay.

SERS probing of fungal HET-s fibrils formed at neutral and acidic pH conditions

Advances in precision medical diagnostics require accurate and sensitive characterization of pathogens. In particular, health conditions associated with protein misfolding require an identification of proteinaceous amyloid fibrils or their precursors. These pathogenic entities express specific molecular structures, which require ultra-sensitive, molecular-level detection methods. A potentially transformative technique termed nanoplasmonics employs electro-optical phenomena in the vicinity of specially engineered metal nanostructures. A signature application of nanoplasmonics exploits enhancement of inelastic scattering of light in specific locations near metallic nanostructures, known as surface-enhanced Raman scattering (SERS). We applied SERS complemented with confocal microscopy imaging for ultra-sensitive, non-invasive, and label-free characterization of the fungal prion HET-s (218-289) as a model for β-sheet rich amyloid structures. This characterization employed Au-coated dielectric supports as plasmonic substrates. After confirming the formation of HET-s fibrils at both pH 7.5 and 2.8 using negative staining transmission electron microscopy, we subjected the fibril-containing solutions to multimodal analysis using confocal microscopy and SERS. The SERS spectral fingerprints from all HET-s samples expressed vibrational markers for β-structure, unstructured backbone, and aromatic side-chains. However, relative intensities of major SERS bands were pronouncedly different for the two pH levels. We have analyzed potential origins of the most pronounced SERS bands and proposed hypothetical mechanistic models that could explain the observed SERS fingerprints from HET-s fibrils grown at pH 7.5 and 2.8.

Manipulating multi-level selection in a fungal entomopathogen reveals social conflicts and a method for improving biocontrol traits

Changes in parasite virulence are commonly expected to lead to trade-offs in other life history traits that can affect fitness. Understanding these trade-offs is particularly important if we want to manipulate the virulence of microbial biological control agents. Theoretically, selection across different spatial scales, i.e. between- and within-hosts, shapes these trade-offs. However, trade-offs are also dependent on parasite biology. Despite their applied importance the evolution of virulence in fungal parasites is poorly understood: virulence can be unstable in culture and commonly fails to increase in simple passage experiments. We hypothesized that manipulating selection intensity at different scales would reveal virulence trade-offs in a fungal pathogen of aphids, Akanthomyces muscarius. Starting with a genetically diverse stock we selected for speed of kill, parasite yield or infectivity by manipulating competition within and between hosts and between-populations of hosts over 7 rounds of infection. We characterized ancestral and evolved lineages by whole genome sequencing and by measuring virulence, growth rate, sporulation and fitness. While several lineages showed increases in virulence, we saw none of the trade-offs commonly found in obligately-killing parasites. Phenotypically similar lineages within treatments often shared multiple single-nucleotide variants, indicating strong convergent evolution. The most dramatic phenotypic changes were in timing of sporulation and spore production in vitro. We found that early sporulation led to reduced competitive fitness but could increase yield of spores on media, a trade-off characteristic of social conflict. Notably, the selection regime with strongest between-population competition and lowest genetic diversity produced the most consistent shift to early sporulation, as predicted by social evolution theory. Multi-level selection therefore revealed social interactions novel to fungi and showed that these biocontrol agents have the genomic flexibility to improve multiple traits-virulence and spore production-that are often in conflict in other parasites.

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.

Emergence of the fungal immune system

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

Allorecognition genes drive reproductive isolation in Podospora anserina

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

Large-scale fungal strain sequencing unravels the molecular diversity in mating loci maintained by long-term balancing selection

Balancing selection, an evolutionary force that retains genetic diversity, has been detected in multiple genes and organisms, such as the sexual mating loci in fungi. However, to quantify the strength of balancing selection and define the mating-related genes require a large number of strains. In tetrapolar basidiomycete fungi, sexual type is determined by two unlinked loci, MATA and MATB. Genes in both loci define mating type identity, control successful mating and completion of the life cycle. These loci are usually highly diverse. Previous studies have speculated, based on culture crosses, that species of the non-model genus Trichaptum (Hymenochaetales, Basidiomycota) possess a tetrapolar mating system, with multiple alleles. Here, we sequenced a hundred and eighty strains of three Trichaptum species. We characterized the chromosomal location of MATA and MATB, the molecular structure of MAT regions and their allelic richness. The sequencing effort was sufficient to molecularly characterize multiple MAT alleles segregating before the speciation event of Trichaptum species. Analyses suggested that long-term balancing selection has generated trans-species polymorphisms. Mating sequences were classified in different allelic classes based on an amino acid identity (AAI) threshold supported by phylogenetics. 17,550 mating types were predicted based on the allelic classes. In vitro crosses allowed us to support the degree of allelic divergence needed for successful mating. Even with the high amount of divergence, key amino acids in functional domains are conserved. We conclude that the genetic diversity of mating loci in Trichaptum is due to long-term balancing selection, with limited recombination and duplication activity. The large number of sequenced strains highlighted the importance of sequencing multiple individuals from different species to detect the mating-related genes, the mechanisms generating diversity and the evolutionary forces maintaining them.

Cytoplasmic Mixing, Not Nuclear Coexistence, Can Explain Somatic Incompatibility in Basidiomycetes

Nonself recognition leading to somatic incompatibility (SI) is commonly used by mycologists to distinguish fungal individuals. Despite this, the process remains poorly understood in basidiomycetes as all current models of SI are based on genetic and molecular research in ascomycete fungi. Ascomycete fungi are mainly found in a monokaryotic stage, with a single type of haploid nuclei, and only briefly during mating do two genomes coexist in heterokaryotic cells. The sister phylum, Basidiomycota, differs in several relevant aspects. Basidiomycete fungi have an extended heterokaryotic stage, and SI is generally observed between heterokaryons instead of between homokaryons. Additionally, considerable nuclear migration occurs during a basidiomycete mating reaction, introducing a nucleus into a resident homokaryon with cytoplasmic mixing limited to the fused or neighboring cells. To accommodate these differences, we describe a basidiomycete model for nonself recognition using post-translational modification, based on a reader-writer system as found in other organisms. This post-translational modification combined with nuclear migration allows for the coexistence of two genomes in one individual while maintaining nonself recognition during all life stages. Somewhat surprisingly, this model predicts localized cell death during mating, which is consistent with previous observations but differs from the general assumptions of basidiomycete mating. This model will help guide future research into the mechanisms behind basidiomycete nonself recognition.

The Predicted Mannosyltransferase GT69-2 Antagonizes RFW-1 To Regulate Cell Fusion in Neurospora crassa

Filamentous fungi undergo somatic cell fusion to create a syncytial, interconnected hyphal network which confers a fitness benefit during colony establishment. However, barriers to somatic cell fusion between genetically different cells have evolved that reduce invasion by parasites or exploitation by maladapted genetic entities (cheaters). Here, we identified a predicted mannosyltransferase, glycosyltransferase family 69 protein (GT69-2) that was required for somatic cell fusion in Neurospora crassa Cells lacking GT69-2 prematurely ceased chemotropic signaling and failed to complete cell wall dissolution and membrane merger in pairings with wild-type cells or between Δgt69-2 cells (self fusion). However, loss-of-function mutations in the linked regulator of cell fusion and cell wall remodeling-1 (rfw-1) locus suppressed the self-cell-fusion defects of Δgt69-2 cells, although Δgt69-2 Δrfw-1 double mutants still failed to undergo fusion with wild-type cells. Both GT69-2 and RFW-1 localized to the Golgi apparatus. Genetic analyses indicated that RFW-1 negatively regulates cell wall remodeling-dependent processes, including cell wall dissolution during cell fusion, separation of conidia during asexual sporulation, and conidial germination. GT69-2 acts as an antagonizer to relieve or prevent negative functions on cell fusion by RFW-1. In Neurospora species and N. crassa populations, alleles of gt69-2 were highly polymorphic and fell into two discrete haplogroups. In all isolates within haplogroup I, rfw-1 was conserved and linked to gt69-2 All isolates within haplogroup II lacked rfw-1. These data indicated that gt69-2/rfw-1 are under balancing selection and provide new mechanisms regulating cell wall remodeling during cell fusion and conidial separation.IMPORTANCE Cell wall remodeling is a dynamic process that balances cell wall integrity versus cell wall dissolution. In filamentous fungi, cell wall dissolution is required for somatic cell fusion and conidial separation during asexual sporulation. In the filamentous fungus Neurospora crassa, allorecognition checkpoints regulate the cell fusion process between genetically different cells. Our study revealed two linked loci with transspecies polymorphisms and under coevolution, rfw-1 and gt69-2, which form a coordinated system to regulate cell wall remodeling during somatic cell fusion, conidial separation, and asexual spore germination. RFW-1 acts as a negative regulator of these three processes, while GT69-2 functions antagonistically to RFW-1. Our findings provide new insight into the mechanisms involved in regulation of fungal cell wall remodeling during growth and development.

An atlas of CNV maps in cattle, goat and sheep

Copy number variation (CNV) is the most prevalent type of genetic structural variation that has been recognized as an important source of phenotypic variation in humans, animals and plants. However, the mechanisms underlying the evolution of CNVs and their function in natural or artificial selection remain unknown. Here, we generated CNV region (CNVR) datasets which were diverged or shared among cattle, goat, and sheep, including 886 individuals from 171 diverse populations. Using 9 environmental factors for genome-wide association study (GWAS), we identified a series of candidate CNVRs, including genes relating to immunity, tick resistance, multi-drug resistance, and muscle development. The number of CNVRs shared between species is significantly higher than expected (P<0.00001), and these CNVRs may be more persist than the single nucleotide polymorphisms (SNPs) shared between species. We also identified genomic regions under long-term balancing selection and uncovered the potential diversity of the selected CNVRs close to the important functional genes. This study provides the evidence that balancing selection might be more common in mammals than previously considered, and might play an important role in the daily activities of these ruminant species.

Creating novel specificities in a fungal nonself recognition system by single step homologous recombination events

In many organisms, two component systems have evolved to discriminate self from nonself. While the molecular function of the two components has been elucidated in several systems, the evolutionary events leading to the large number of different specificities for self-nonself recognition found in most systems remain obscure. We have investigated the variation within a multiallelic nonself recognition system in the phytopathogenic basidiomycete Ustilago maydis by means of sequence analysis and functional studies. The multiallelic b mating type locus of U. maydis ensures outbreeding during sexual development. Nonself recognition is specified by the two homeodomain proteins, bE and bW, encoded by the b locus. While bE-bW combinations from the same allele do not dimerize, bE and bW proteins originating from different alleles form a heterodimeric complex that functions as master regulator for sexual and pathogenic development. We show that novel specificities of the b mating type locus have arisen by single homologous recombination events between distinct b alleles that lead to a simultaneous exchange of subdomains involved in dimerization in both bE and bW, altering the specificity of both proteins in a single step.

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.

Proteomic investigation of interhyphal interactions between strains of Agaricus bisporus

Hyphae of filamentous fungi undergo polar extension, bifurcation and hyphal fusion to form reticulating networks of mycelia. Hyphal fusion or anastomosis, a ubiquitous process among filamentous fungi, is a vital strategy for how fungi expand over their substrate and interact with or recognise self- and non-self hyphae of neighbouring mycelia in their environment. Morphological and genetic characterisation of anastomosis has been studied in many model fungal species, but little is known of the direct proteomic response of two interacting fungal isolates. Agaricus bisporus, the most widely cultivated edible mushroom crop worldwide, was used as an in vitro model to profile the proteomes of interacting cultures. The globally cultivated strain (A15) was paired with two distinct strains; a commercial hybrid strain and a wild isolate strain. Each co-culture presented a different interaction ranging from complete vegetative compatibility (self), lack of interactions, and antagonistic interactions. These incompatible strains are the focus of research into disease-resistance in commercial crops as the spread of intracellular pathogens, namely mycoviruses, is limited by the lack of interhyphal anastomosis. Unique proteomic responses were detected between all co-cultures. An array of cell wall modifying enzymes, plus fungal growth and morphogenesis proteins were found in significantly (P < 0.05) altered abundances. Nitrogen metabolism dominated in the intracellular proteome, with evidence of nitrogen starvation between competing, non-compatible cultures. Changes in key enzymes of A. bisporus morphogenesis were observed, particularly via increased abundance of glucanosyltransferase in competing interactions and certain chitinases in vegetative compatible interactions only. Carbohydrate-active enzyme arsenals are expanded in antagonistic interactions in A. bisporus. Pathways involved in carbohydrate metabolism and genetic information processing were higher in interacting cultures, most notably during self-recognition. New insights into the differential response of interacting strains of A. bisporus will enhance our understanding of potential barriers to viral transmission through vegetative incompatibility. Our results suggest that a differential proteomic response occurs between A. bisporus at strain-level and findings from this work may guide future proteomic investigation of fungal anastomosis.

Programmed Cell Death in Neurospora crassa Is Controlled by the Allorecognition Determinant rcd-1

Nonself recognition following cell fusion between genetically distinct individuals of the same species in filamentous fungi often results in a programmed cell death (PCD) reaction, where the heterokaryotic fusion cell is compartmentalized and rapidly killed. The allorecognition process plays a key role as a defense mechanism that restricts genome exploitation, resource plundering, and the spread of deleterious senescence plasmids and mycoviruses. Although a number of incompatibility systems have been described that function in mature hyphae, less is known about the PCD pathways in asexual spores, which represent the main infectious unit in various human and plant fungal pathogens. Here, we report the identification of regulator of cell death-1 (rcd-1), a novel allorecognition gene, controlling PCD in germinating asexual spores of Neurospora crassa; rcd-1 is one of the most polymorphic genes in the genomes of wild N. crassa isolates. The coexpression of two antagonistic rcd-1-1 and rcd-1-2 alleles was necessary and sufficient to trigger cell death in fused germlings and in hyphae. Based on analysis of wild populations of N. crassa and N. discreta, rcd-1 alleles appeared to be under balancing selection and associated with trans-species polymorphisms. We shed light on genomic rearrangements that could have led to the emergence of the incompatibility system in Neurospora and show that rcd-1 belongs to a much larger gene family in fungi. Overall, our work contributes toward a better understanding of allorecognition and PCD in an underexplored developmental stage of filamentous fungi.

Bioinformatic Analysis of the Type VI Secretion System and Its Potential Toxins in the Acinetobacter Genus

Several Acinetobacter strains are important nosocomial pathogens, with Acinetobacter baumannii as the species of greatest concern worldwide due to its multi-drug resistance and recent appearance of hyper-virulent strains in the clinical setting. Acinetobacter colonization of the environment and the host is associated with a multitude of factors which remain poorly characterized. Among them, the secretion systems (SS) encoded by Acinetobacter species confer adaptive advantages depending on the niche occupied. Different SS have been characterized in this group of microorganisms, including T6SS used by several Acinetobacter species to outcompete other bacteria and in some A. baumannii strains for Galleria mellonella colonization. Therefore, to better understand the distribution of the T6SS in this genus we carried out an in-depth comparative genomic analysis of the T6SS in 191 sequenced strains. To this end, we analyzed the gene content, sequence similarity, synteny and operon structure of each T6SS loci. The presence of a single conserved T6SS-main cluster (T6SS-1), with two different genetic organizations, was detected in the genomes of several ecologically diverse species. Furthermore, a second main cluster (T6SS-2) was detected in a subgroup of 3 species of environmental origin. Detailed analysis also showed an impressive genetic versatility in T6SS-associated islands, carrying VgrG, PAAR and putative toxin-encoding genes. This in silico study represents the first detailed intra-species comparative analysis of T6SS-associated genes in the Acinetobacter genus, that should contribute to the future experimental characterization of T6SS proteins and effectors.

Long-term balancing selection drives evolution of immunity genes in Capsella

Genetic drift is expected to remove polymorphism from populations over long periods of time, with the rate of polymorphism loss being accelerated when species experience strong reductions in population size. Adaptive forces that maintain genetic variation in populations, or balancing selection, might counteract this process. To understand the extent to which natural selection can drive the retention of genetic diversity, we document genomic variability after two parallel species-wide bottlenecks in the genus Capsella. We find that ancestral variation preferentially persists at immunity related loci, and that the same collection of alleles has been maintained in different lineages that have been separated for several million years. By reconstructing the evolution of the disease-related locus MLO2b, we find that divergence between ancient haplotypes can be obscured by referenced based re-sequencing methods, and that trans-specific alleles can encode substantially diverged protein sequences. Our data point to long-term balancing selection as an important factor shaping the genetics of immune systems in plants and as the predominant driver of genomic variability after a population bottleneck.

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

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

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.

Long-term balancing selection contributes to adaptation in Arabidopsis and its relatives

BACKGROUND

In contrast to positive selection, which reduces genetic variation by fixing beneficial alleles, balancing selection maintains genetic variation within a population or species and plays crucial roles in adaptation in diverse organisms. However, which genes, genome-wide, are under balancing selection and the extent to which these genes are involved in adaptation are largely unknown.

RESULTS

We performed a genome-wide scan for genes under balancing selection across two plant species, Arabidopsis thaliana and its relative Capsella rubella, which diverged about 8 million generations ago. Among hundreds of genes with shared coding-region polymorphisms, we find evidence for long-term balancing selection in five genes: AT1G35220, AT2G16570, AT4G29360, AT5G38460, and AT5G44000. These genes are involved in the response to biotic and abiotic stress and other fundamental biochemical processes. More intriguingly, for these genes, we detected significant ecological diversification between the two haplotype groups, suggesting that balancing selection has been very important for adaptation.

CONCLUSIONS

Our results indicate that beyond the well-known S-locus genes and resistance genes, many loci are under balancing selection. These genes are mostly correlated with resistance to stress or other fundamental functions and likely play a more important role in adaptation to diverse habitats than previously thought.

Balancing selection on allorecognition genes in the colonial ascidian Botryllus schlosseri

Allorecognition is the capability of an organism to recognize its own or related tissues. The colonial ascidian Botryllus schlosseri, which comprises five genetically distinct and divergent species (Clades A-E), contains two adjacent genes that control allorecognition: fuhc^sec and fuhc^tm. These genes have been characterized extensively in Clade A and are highly polymorphic. Using alleles from 10 populations across the range of Clade A, we investigated the type and strength of selection maintaining this variation. Both fuhc genes exhibit higher within-population variation and lower population differentiation measures (F_ST) than neutral loci. The fuhc genes contain a substantial number of codons with >95% posterior probability of d_N/d_S > 1. fuhc^sec and fuhc^tm also have polymorphisms shared between Clade A and Clade E that were present prior to speciation (trans-species polymorphisms). These results provide robust evidence that the fuhc genes are evolving under balancing selection.

Characterization of the Sclerotinia sclerotiorum cell wall proteome

We used a proteomic analysis to identify cell wall proteins released from Sclerotinia sclerotiorum hyphal and sclerotial cell walls via a trifluoromethanesulfonic acid (TFMS) digestion. Cell walls from hyphae grown in Vogel's glucose medium (a synthetic medium lacking plant materials), from hyphae grown in potato dextrose broth and from sclerotia produced on potato dextrose agar were used in the analysis. Under the conditions used, TFMS digests the glycosidic linkages in the cell walls to release intact cell wall proteins. The analysis identified 24 glycosylphosphatidylinositol (GPI)-anchored cell wall proteins and 30 non-GPI-anchored cell wall proteins. We found that the cell walls contained an array of cell wall biosynthetic enzymes similar to those found in the cell walls of other fungi. When comparing the proteins in hyphal cell walls grown in potato dextrose broth with those in hyphal cell walls grown in the absence of plant material, it was found that a core group of cell wall biosynthetic proteins and some proteins associated with pathogenicity (secreted cellulases, pectin lyases, glucosidases and proteases) were expressed in both types of hyphae. The hyphae grown in potato dextrose broth contained a number of additional proteins (laccases, oxalate decarboxylase, peroxidase, polysaccharide deacetylase and several proteins unique to Sclerotinia and Botrytis) that might facilitate growth on a plant host. A comparison of the proteins in the sclerotial cell wall with the proteins in the hyphal cell wall demonstrated that sclerotia formation is not marked by a major shift in the composition of cell wall protein. We found that the S. sclerotiorum cell walls contained 11 cell wall proteins that were encoded only in Sclerotinia and Botrytis genomes.

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.

Large-scale genotyping of highly polymorphic loci by next-generation sequencing: how to overcome the challenges to reliably genotype individuals?

Studying the different roles of adaptive genes is still a challenge in evolutionary ecology and requires reliable genotyping of large numbers of individuals. Next-generation sequencing (NGS) techniques enable such large-scale sequencing, but stringent data processing is required. Here, we develop an easy to use methodology to process amplicon-based NGS data and we apply this methodology to reliably genotype four major histocompatibility complex (MHC) loci belonging to MHC class I and II of Alpine marmots (Marmota marmota). Our post-processing methodology allowed us to increase the number of retained reads. The quality of genotype assignment was further assessed using three independent validation procedures. A total of 3069 high-quality MHC genotypes were obtained at four MHC loci for 863 Alpine marmots with a genotype assignment error rate estimated as 0.21%. The proposed methodology could be applied to any genetic system and any organism, except when extensive copy-number variation occurs (that is, genes with a variable number of copies in the genotype of an individual). Our results highlight the potential of amplicon-based NGS techniques combined with adequate post-processing to obtain the large-scale highly reliable genotypes needed to understand the evolution of highly polymorphic functional genes.

Nuclear interactions in a heterokaryon: insight from the model Neurospora tetrasperma

A heterokaryon is a tissue type composed of cells containing genetically different nuclei. Although heterokaryosis is commonly found in nature, an understanding of the evolutionary implications of this phenomenon is largely lacking. Here, we use the filamentous ascomycete Neurospora tetrasperma to study the interplay between nuclei in heterokaryons across vegetative and sexual developmental stages. This fungus harbours nuclei of two opposite mating types (mat A and mat a) in the same cell and is thereby self-fertile. We used pyrosequencing of mat-linked SNPs of three heterokaryons to demonstrate that the nuclear ratio is consistently biased for mat A-nuclei during mycelial growth (mean mat A/mat a ratio 87%), but evens out during sexual development (ratio ranging from 40 to 57%). Furthermore, we investigated the association between nuclear ratio and expression of alleles of mat-linked genes and found that expression is coregulated to obtain a tissue-specific bias in expression ratio: during mycelial extension, we found a strong bias in expression for mat A-linked genes, that was independent of nuclear ratio, whereas at the sexual stage we found an expression bias for genes of the mat a nuclei. Taken together, our data indicate that nuclei cooperate to optimize the fitness of the heterokaryon, via both altering their nuclear ratios and coregulation genes expressed in the different nuclei.

Advantageous diversity maintained by balancing selection in humans

Most human polymorphisms are neutral or slightly deleterious, but some genetic variation is advantageous and maintained in populations by balancing selection. Considered a rarity and overlooked for years, balanced polymorphisms have recently received renewed attention with several lines of evidence showing their relevance in human evolution. From theoretical work on its role in adaptation to empirical studies that identify its targets, recent developments have showed that balancing selection is more prevalent than previously thought. Here we review these developments and discuss their implications in our understanding of the influence of balancing selection in human evolution. We also review existing evidence on the biological functions that benefit most from advantageous diversity, and the functional consequences of these variants. Overall, we argue that balancing selection must be considered an important selective force in human evolution.

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.

Hyperparasites influence population structure of the chestnut blight pathogen, Cryphonectria parasitica

Vegetative compatibility (VC) is commonly used to characterize structure and diversity in fungal populations. In the chestnut blight fungus, Cryphonectria parasitica, high VC diversity is hypothesized to be responsible for the failure of hyperparasitic mycoviruses to spread through pathogen populations in North America. To test this hypothesis, we assessed VC diversity at three recovering sites in Michigan where mycoviruses had invaded and compared them with four epidemic population sites where mycoviruses were absent. VC diversity was assessed for samples collected in 1996 and 2009, which allowed us to determine how C. parasitica populations changed with time. Twelve VC types were found in 1996 while 29 were found in 2009; 75% of types overlapped between the sample dates. Sites where mycoviruses were present had unique VC structures with the exception of the recovering population site at County Line where the main VC group was also detected at two epidemic sites. With one exception, epidemic sites contained more VC groups and displayed higher population level diversity than recovering sites. Mating-type analyses of blight populations revealed that two of three recovering populations were significantly skewed for MAT2 suggesting asexual reproduction, while epidemic sites with a long history of blight infection had ratios near 50:50 suggesting sexual reproduction. We propose that selection in the largely asexual C. parasitica populations at two recovering sites favors the most-fit fungal genotype by mycovirus combination and results in reduced diversity relative to the sexually reproducing pathogen populations at epidemic sites.

Extensive trans-specific polymorphism at the mating type locus of the root decay fungus Heterobasidion

Incompatibility systems in which individuals bearing identical alleles reject each other favor the maintenance of a diversity of alleles. Mushroom mating type loci (MAT) encode for dozens or hundreds of incompatibility alleles whose loss from the population is greatly restricted through negative frequency selection, leading to a system of alleles with highly divergent sequences. Here, we use DNA sequences of homeodomain (HD) encoding genes at the MAT locus of five closely related species of the root rot basidiomycete Heterobasidion annosum sensu lato to show that the extended coalescence time of MAT alleles greatly predates speciation in the group, contrasting loci outside of MAT that show allele divergences largely consistent with the species phylogeny with those of MAT, which show rampant trans-species polymorphism. We observe a roughly 6-fold greater genealogical depth and polymorphism of MAT compared with non-MAT that argues for the maintenance of balanced polymorphism for a minimum duration of 24 My based on a molecular-clock calibrated species phylogeny. As with other basidiomycete HD genes, balancing selection appears to be concentrated at the specificity-determining region in the N-terminus of the protein based on identification of codons under selection and the absence of recombination within the region. However, the elevated polymorphism extends into the nonspecificity determining regions as well as a neighboring non-MAT gene, the mitochondrial intermediate peptidase (MIP). In doing so, increased divergence should decrease recombination among alleles and as a by-product create incompatibilities in the functional domains not involved in allele recognition but in regulating sexual development.

The Loci of repeated evolution: a catalog of genetic hotspots of phenotypic variation

What is the nature of the genetic changes underlying phenotypic evolution? We have catalogued 1008 alleles described in the literature that cause phenotypic differences among animals, plants, and yeasts. Surprisingly, evolution of similar traits in distinct lineages often involves mutations in the same gene ("gene reuse"). This compilation yields three important qualitative implications about repeated evolution. First, the apparent evolution of similar traits by gene reuse can be traced back to two alternatives, either several independent causative mutations or a single original mutational event followed by sorting processes. Second, hotspots of evolution-defined as the repeated occurrence of de novo mutations at orthologous loci and causing similar phenotypic variation-are omnipresent in the literature with more than 100 examples covering various levels of analysis, including numerous gain-of-function events. Finally, several alleles of large effect have been shown to result from the aggregation of multiple small-effect mutations at the same hotspot locus, thus reconciling micromutationist theories of adaptation with the empirical observation of large-effect variants. Although data heterogeneity and experimental biases prevented us from extracting quantitative trends, our synthesis highlights the existence of genetic paths of least resistance leading to viable evolutionary change.

The fester locus in Botryllus schlosseri experiences selection

Background

Allorecognition, the ability of an organism to distinguish self from non-self, occurs throughout the entire tree of life. Despite the prevalence and importance of allorecognition systems, the genetic basis of allorecognition has rarely been characterized outside the well-known MHC (Major Histocompatibility Complex) in vertebrates and SI (Self-Incompatibility) in plants. Where loci have been identified, their evolutionary history is an open question. We have previously identified the genes involved in self/non-self recognition in the colonial ascidian Botryllus schlosseri, and we can now begin to investigate their evolution. In B. schlosseri, colonies sharing 1 or more alleles of a gene called FuHC (Fusion Histocompatibility) will fuse. Protein products of a locus called fester, located ~300 kb from FuHC, have been shown to play multiple roles in the histocompatibility reaction, as activating and/or inhibitory receptors. We test whether the proteins encoded by this locus are evolving neutrally or are experiencing balancing, directional, or purifying selection.

Results

Nearly all of the variation in the fester locus resides within populations. The 13 housekeeping genes (12 nuclear genes and mitochondrial cytochrome oxidase I) have substantially more structure among populations within groups and among groups than fester. All polymorphism statistics (Tajima's D, Fu and Li's D* and F*) are significantly negative for the East Coast A-type alleles, and Fu and Li's F* statistic is significantly negative for the West Coast A-type alleles. These results are likely due to selection rather than demography, given that 10 of the housekeeping loci have no populations with significant values for any of the polymorphism statistics. The majority of codons in the fester proteins have ω values < 1, but 15–27 codons have > 95% posterior probability of ω values > 1.

Conclusion

Fester proteins are evolving non-neutrally. The polymorphism statistics are consistent with either purifying selection or directional selection. The ω statistics show that the majority of the protein is experiencing purifying selection (ω < 1), but that 15–27 codons are undergoing either balancing or directional selection: ω > 1 is compatible with either scenario. The distribution of variation within and among populations points towards balancing selection and away from directional selection. While these data do not provide unambiguous support for a specific type of selection, they contribute to our evolutionary understanding of a critical biological process by determining the forces that affect loci involved in allorecognition.

Recent and ancient signature of balancing selection around the S-locus in Arabidopsis halleri and A. lyrata

Balancing selection can maintain different alleles over long evolutionary times. Beyond this direct effect on the molecular targets of selection, balancing selection is also expected to increase neutral polymorphism in linked genome regions, in inverse proportion to their genetic map distances from the selected sites. The genes controlling plant self-incompatibility are subject to one of the strongest forms of balancing selection, and they show clear signatures of balancing selection. The genome region containing those genes (the S-locus) is generally described as nonrecombining, and the physical size of the region with low recombination has recently been established in a few species. However, the size of the region showing the indirect footprints of selection due to linkage to the S-locus is only roughly known. Here, we improved estimates of this region by surveying synonymous polymorphism and estimating recombination rates at 12 flanking region loci at known physical distances from the S-locus region boundary, in two closely related self-incompatible plants Arabidopsis halleri and A. lyrata. In addition to studying more loci than previous studies and using known physical distances, we simulated an explicit demographic scenario for the divergence between the two species, to evaluate the extent of the genomic region whose diversity departs significantly from neutral expectations. At the closest flanking loci, we detected signatures of both recent and ancient indirect effects of selection on the S-locus flanking genes, finding ancestral polymorphisms shared by both species, as well as an excess of derived mutations private to either species. However, these effects are detected only in a physically small region, suggesting that recombination in the flanking regions is sufficient to quickly break up linkage disequilibrium with the S-locus. Our approach may be useful for distinguishing cases of ancient versus recently evolved balancing selection in other systems.

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.

High natural prevalence of a fungal prion

Prions are infectious proteins that cause fatal diseases in mammals. Prions have also been found in fungi, but studies on their role in nature are scarce. The proposed biological function of fungal prions is debated and varies from detrimental to benign or even beneficial. [Het-s] is a prion of the fungus Podospora anserina. The het-s locus exists as two antagonistic alleles that constitute an allorecognition system: the het-s allele encoding the protein variant capable of prion formation and the het-S allele encoding a protein variant that cannot form a prion. We document here that het-s alleles, capable of prion formation, are nearly twice as frequent as het-S alleles in a natural population of 112 individuals. Then, we report a 92% prevalence of [Het-s] prion infection among the het-s isolates and find evidence of the role of the [Het-s]/het-S allorecognition system on the incidence of infection by a deleterious senescence plasmid. We explain the het-s/het-S allele ratios by the existence of two selective forces operating at different levels. We propose that during the somatic stage, the role of [Het-s]/HET-S in allorecognition leads to frequency-dependent selection for which an equilibrated frequency would be optimal. However, in the sexual cycle, the [Het-s] prion causes meiotic drive favoring the het-s allele. Our findings indicate that [Het-s] is a selected and, therefore, widespread prion whose activity as selfish genetic element is counteracted by balancing selection for allorecognition polymorphism.

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.

Parasitism and maintenance of diversity in a fungal vegetative incompatibility system: the role of selection by deleterious cytoplasmic elements

In fungi, horizontal transmission of deleterious cytoplasmic elements is reduced by the vegetative incompatibility system. This self/non-self recognition system may select for greater diversity of fungal incompatibility phenotypes in a frequency-dependent manner but the link between the diversity of fungal phenotypes and the virulence of cytoplasmic parasites has been poorly studied. We used an epidemiological model to show that even when transmission between incompatibility types is permitted, parasite pressure can lead to high levels of polymorphism for vegetative incompatibility systems. Moreover, high levels of polymorphism in host populations can select for less virulent cytoplasmic parasites. This feedback mechanism between parasite virulence and vegetative incompatibility system polymorphism of host populations may account for the general avirulence of most known mycoviruses. Furthermore, this mechanism provides a new perspective on the particular ecology and evolution of the host/parasite interactions acting between fungi and their cytoplasmic parasites.

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.

A 454 multiplex sequencing method for rapid and reliable genotyping of highly polymorphic genes in large-scale studies

Background

High-throughput sequencing technologies offer new perspectives for biomedical, agronomical and evolutionary research. Promising progresses now concern the application of these technologies to large-scale studies of genetic variation. Such studies require the genotyping of high numbers of samples. This is theoretically possible using 454 pyrosequencing, which generates billions of base pairs of sequence data. However several challenges arise: first in the attribution of each read produced to its original sample, and second, in bioinformatic analyses to distinguish true from artifactual sequence variation. This pilot study proposes a new application for the 454 GS FLX platform, allowing the individual genotyping of thousands of samples in one run. A probabilistic model has been developed to demonstrate the reliability of this method.

Results

DNA amplicons from 1,710 rodent samples were individually barcoded using a combination of tags located in forward and reverse primers. Amplicons consisted in 222 bp fragments corresponding to DRB exon 2, a highly polymorphic gene in mammals. A total of 221,789 reads were obtained, of which 153,349 were finally assigned to original samples. Rules based on a probabilistic model and a four-step procedure, were developed to validate sequences and provide a confidence level for each genotype. The method gave promising results, with the genotyping of DRB exon 2 sequences for 1,407 samples from 24 different rodent species and the sequencing of 392 variants in one half of a 454 run. Using replicates, we estimated that the reproducibility of genotyping reached 95%.

Conclusions

This new approach is a promising alternative to classical methods involving electrophoresis-based techniques for variant separation and cloning-sequencing for sequence determination. The 454 system is less costly and time consuming and may enhance the reliability of genotypes obtained when high numbers of samples are studied. It opens up new perspectives for the study of evolutionary and functional genetics of highly polymorphic genes like major histocompatibility complex genes in vertebrates or loci regulating self-compatibility in plants. Important applications in biomedical research will include the detection of individual variation in disease susceptibility. Similarly, agronomy will benefit from this approach, through the study of genes implicated in productivity or disease susceptibility traits.