The signature of balancing selection: fungal mating compatibility gene evolution

Evolution of unexpected diversity in a putative mating type locus and its correlation with genome variability reveals likely asexuality in the model mycorrhizal fungus Rhizophagus irregularis

BACKGROUND

Arbuscular mycorrhizal fungi (AMF) form mutualistic partnerships with approximately 80% of plant species. AMF, and their diversity, play a fundamental role in plant growth, driving plant diversity, and global carbon cycles. Knowing whether AMF are sexual or asexual has fundamental consequences for how they can be used in agricultural applications. Evidence for and against sexuality in the model AMF, Rhizophagus irregularis, has been proposed. The discovery of a putative mating-type locus (MAT locus) in R. irregularis, and the previously suggested recombination among nuclei of a dikaryon R. irregularis isolate, potentially suggested sexuality. Unless undergoing frequent sexual reproduction, evolution of MAT-locus diversity is expected to be very low. Additionally, in sexual species, MAT-locus evolution is decoupled from the evolution of arbitrary genome-wide loci.

RESULTS

We studied MAT-locus diversity of R. irregularis. This was then compared to diversification in a phosphate transporter gene (PTG), that is not involved in sex, and to genome-wide divergence, defined by 47,378 single nucleotide polymorphisms. Strikingly, we found unexpectedly high MAT-locus diversity indicating that either it is not involved in sex, or that AMF are highly active in sex. However, a strongly congruent evolutionary history of the MAT-locus, PTG and genome-wide arbitrary loci allows us to reject both the hypothesis that the MAT-locus is involved in mating and that the R. irregularis lineage is sexual.

CONCLUSION

Our finding shapes the approach to developing more effective AMF strains and is highly informative as it suggests that introduced strains applied in agriculture will not exchange DNA with native populations.

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.

Unique footprints of balancing selection in bovine genome

Balancing selection is the process of selection that preserves various alleles within a population. Studying the areas undergoing balancing selection is essential, because it preserves genetic diversity in a population. Finding genes that exhibit signs of balancing selection during the domestication of cattle is the goal of this study. To identify regions where polymorphism has persisted in the cattle population for millions of years, we examined the genome of cattle. In this study, we used bovine SNP 50 k data to conduct a detailed genome-wide assessment of selection signatures for balancing selection. We have included the genotyped data from 427 animals, including five taurines, two crossbreds, and eight Indian cattle breeds. For this study, we employed Tajima's D approach to identify signature regions undergoing balancing selection. Using the NCBI database, PANTHER 17.0, and CattleQTL database, the annotation was carried out after finding the relevant areas under balancing selection. The number of genomic regions undergoing balancing selection in Ayrshire, Brown-Swiss, Frieswal, Gir, Guernsey, Hariana, Holstein Friesian, Jersey, Kankrej, Nelore, Ongole, Red Sindhi, Sahiwal, Tharparkar, and Vrindavani was 11, 13, 13, 19, 18, 11, 17, 14, 14, 12, 10, 12, 13, 13, and 11, respectively. We have observed multiple immune system-related genes going through balancing selection, including KIT, NFATC2, GBP4, LRRC32, SYT7, RAG1, RAG2, LOC513659, and ZBTB17. In our study, we found that the majority of the immune-related genes and a few genes associated with growth, reproduction, production, and adaptation are undergoing balancing selection.

Sex Without Sexes: Can the Cost of Finding a Mate Explain Diversity in Fungal Mating Systems?

Eukaryotes have evolved myriad ways of uniting gametes during sexual reproduction. A repeated pattern is the convergent evolution of a mating system with the fusion of larger gametes with smaller gametes (anisogamy) from that of fusion between morphologically identical gametes (isogamy). In anisogamous species, sexes are defined as individuals that produce only one gamete type. Although sexes abound throughout Eukarya, in fungi there are no biological sexes, because even in anisogamous species, individuals are hermaphroditic and produce both gamete types. For this reason, the term mating types is preferred over sexes, and, thus defined, only individuals of differing mating types can mate (homoallelic incompatibility). In anisogamous fungal species, there is scant evidence that there are more than two mating types, and this may be linked to genetic constraints, such as the use of mating types to determine the inheritance of cytoplasmic genomes. However, the mushroom fungi (Agaricomycetes) stand out as having both large numbers of mating types within a species, which will allow nearly all individuals to be compatible with each other, and reciprocal exchange of nuclei during mating, which will avoid cytoplasmic mixing and cyto-nuclear conflicts. Although the limitation of mating types to two in most fungi is consistent with the cyto-nuclear conflicts model, there are many facets of the Agaricomycete life cycle that also suggest they will demand a high outbreeding efficiency. Specifically, they are mostly obligately sexual and outcrossing, inhabit complex competitive niches, and display broadcast spore dispersal. Subsequently, the Agaricomycete individual pays a high cost to being choosy when encountering a mate. Here, I discuss the costs of mate finding and choice and demonstrate how most fungi have multiple ways of reducing these costs, which can explain why mating types are mostly limited to two per species. Nevertheless, it is perplexing that fungi have not evolved multiple mating types on more occasions nor evolved sexes. The few exceptions to these rules suggest that it is dictated by both molecular and evolutionary constraints.

Population structure and mixed reproductive strategies in Bipolaris maydis from single and multiple corn cultivars in Fujian Province, China

Bipolaris maydis is the pathogenic microorganism of southern corn leaf blight, a persistent biotic constraint responsible for substantial yield losses of corn worldwide. In the present study, 96 isolates from six representative fields growing single and multiple sweet corn cultivars in Pingnan, Fuqing, and Jian'ou in Fujian Province, which are characterized by different geographical characteristics and cropping patterns, were genetically analyzed using inter-simple sequence repeat (ISSR) markers to assess the impact of geographical origins and corn cultivars on B. maydis population differentiation. B. maydis isolates originated from diverse regions possessed higher genetic variety than those from single and multiple sweet corn cultivars. Phylogenetic analysis showed that the isolates from single and multiple sweet corn cultivars were randomly grouped into different clusters, with those from the same location tending to form clusters. A greater genetic differentiation among different geographical populations than between those from single and multiple sweet corn cultivars was observed by pairwise comparison. Hierarchical analysis indicated that among-population variation was higher when comparatively analyzed B. maydis populations from different locations than in those from single and multiple sweet corn cultivars. In conclusion, these results suggest that geographical origin acts a more considerable role in genetic differentiation of B. maydis than corn cultivar. Two divided genetic clusters were detected in the B. maydis populations from single and multiple sweet corn cultivars at the three locations in Fujian Province, with major genetic variation being derived within populations. The high haplotypic diversity and expected mating type ratio of 1:1 in combination with significant linkage disequilibrium suggested that a mixed reproductive strategy occurs in the B. maydis population in Fujian Province. This study will enrich the information on the role that geographical origins and corn cultivars play in the population structure of the pathogen as well as the reproductive strategies in B. maydis population in Fujian Province.

Introgression between highly divergent fungal sister species

To understand how species evolve and adapt to changing environments, it is important to study gene flow and introgression due to their influence on speciation and radiation events. Here, we apply a novel experimental system for investigating these mechanisms using natural populations. The system is based on two fungal sister species with morphological and ecological similarities occurring in overlapping habitats. We examined introgression between these species by conducting whole genome sequencing of individuals from populations in North America and Europe. We assessed genome-wide nucleotide divergence and performed crossing experiments to study reproductive barriers. We further used ABBA-BABA statistics together with a network analysis to investigate introgression, and conducted demographic modelling to gain insight into divergence times and introgression events. The results revealed that the species are highly divergent and incompatible in vitro. Despite this, small regions of introgression were scattered throughout the genomes and one introgression event likely involves a ghost population (extant or extinct). This study demonstrates that introgression can be found among divergent species and that population histories can be studied without collections of all the populations involved. Moreover, the experimental system is shown to be a useful tool for research on reproductive isolation in natural populations.

Genomic determination of breeding systems and trans-specific evolution of HD MAT genes in suilloid fungi

Studying the signatures of evolution can help to understand genetic processes. Here, we demonstrate how the existence of balancing selection can be used to identify the breeding systems of fungi from genomic data. The breeding systems of fungi are controlled by self-incompatibility loci that determine mating types between potential mating partners, resulting in strong balancing selection at the loci. Within the fungal phylum Basidiomycota, two such self-incompatibility loci, namely HD MAT locus and P/R MAT locus, control mating types of gametes. Loss of function at one or both MAT loci results in different breeding systems and relaxes the MAT locus from balancing selection. By investigating the signatures of balancing selection at MAT loci, one can infer a species' breeding system without culture-based studies. Nevertheless, the extreme sequence divergence among MAT alleles imposes challenges for retrieving full variants from both alleles when using the conventional read-mapping method. Therefore, we employed a combination of read-mapping and local de novo assembly to construct haplotypes of HD MAT alleles from genomes in suilloid fungi (genera Suillus and Rhizopogon). Genealogy and pairwise divergence of HD MAT alleles showed that the origins of mating types predate the split between these two closely related genera. High sequence divergence, trans-specific polymorphism, and the deeply diverging genealogy confirm the long-term functionality and multiallelic status of HD MAT locus in suilloid fungi. This work highlights a genomics approach to studying breeding systems regardless of the culturability of organisms based on the interplay between evolution and genetics.

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.

Genetic Differentiation and Mixed Reproductive Strategies in the Northern Corn Leaf Blight Pathogen Setosphaeria turcica From Sweet Corn in Fujian Province, China

The northern corn leaf blight (NCLB) pathogen Setosphaeria turcica (Luttrell) Leonard and Suggs is one of the main biotic constraints on sweet corn (Zea mays L.) yield and quality in Fujian Province, China. Currently, however, there is comparatively little information available regarding the distribution of mating types, population genetics, and reproductive strategies of this pathogen in Fujian. In this study, we investigated the distribution of mating types and population genetics of 117 isolates of S. turcica collected from seven of the main sweet corn-growing regions in Fujian Province, based on multiple polymerase chain reaction analyses using two mating type-specific primer pairs and 11 inter-simple sequence repeat markers. Furthermore, we examined the mode of reproduction of Fujian S. turcica populations. Both MAT1-1 and MAT1-2 mating types were detected throughout all seven sampling locations. The majority of MAT1-2 isolates were detected from Dongyou, Jian’ou, Pingnan, Songxi, and Longyan, whereas a large proportion of the detected MAT1-1 isolates were among those collected from Dongfeng and Nanjing. Furthermore, we detected five shared multi-locus haplotypes among S. turcica isolates from Dongyou, Jian’ou, Pingnan, Nanjing, and Songxi, whereas no shared haplotypes were observed between the Dongfeng (or Longyan) population and these five populations. Pairwise comparisons of the indices Φ_PT and Nm, and population structure and principal coordinate analyses indicated genetic differentiation between both the regional and the mating type populations of S. turcica in Fujian. The skewed mating type ratio associated with low a haplotypic diversity and evident linkage disequilibrium reveals a mixed reproductive strategy for S. turcica populations in Fujian Province. The findings of this study advance our current understanding of the genetic diversity, population structure, and reproductive strategies of S. turcica populations infecting sweet corn in Fujian Province, and will potentially contribute to further resistance breeding efforts.

Population Genomics of the Maize Pathogen Ustilago maydis: Demographic History and Role of Virulence Clusters in Adaptation

The tight interaction between pathogens and their hosts results in reciprocal selective forces that impact the genetic diversity of the interacting species. The footprints of this selection differ between pathosystems because of distinct life-history traits, demographic histories, or genome architectures. Here, we studied the genome-wide patterns of genetic diversity of 22 isolates of the causative agent of the corn smut disease, Ustilago maydis, originating from five locations in Mexico, the presumed center of origin of this species. In this species, many genes encoding secreted effector proteins reside in so-called virulence clusters in the genome, an arrangement that is so far not found in other filamentous plant pathogens. Using a combination of population genomic statistical analyses, we assessed the geographical, historical, and genome-wide variation of genetic diversity in this fungal pathogen. We report evidence of two partially admixed subpopulations that are only loosely associated with geographic origin. Using the multiple sequentially Markov coalescent model, we inferred the demographic history of the two pathogen subpopulations over the last 0.5 Myr. We show that both populations experienced a recent strong bottleneck starting around 10,000 years ago, coinciding with the assumed time of maize domestication. Although the genome average genetic diversity is low compared with other fungal pathogens, we estimated that the rate of nonsynonymous adaptive substitutions is three times higher in genes located within virulence clusters compared with nonclustered genes, including nonclustered effector genes. These results highlight the role that these singular genomic regions play in the evolution of this pathogen.

Population genomics and evolution of a fungal pathogen after releasing exotic strains to control insect pests for 20 years

Entomopathogenic fungi are one of the key regulators of insect populations in nature. Some species such as Beauveria bassiana with a wide host range have been developed as promising alternatives to chemical insecticides for the biocontrol of insect pests. However, the long-term persistence of the released strains, the effect on non-target hosts and local fungal populations remains elusive, but they are considerable concerns with respect to environmental safety. Here we report the temporal features of the Beauveria population genomics and evolution over 20 years after releasing exotic strains to control pine caterpillar pests. We found that the isolates within the biocontrol site were mostly of clonal origins. The released strains could persist in the environment for a long time but with low recovery rates. Similar to the reoccurrence of host jumping by local isolates, the infection of non-target insects by the released strains was evident to endemically occur in association with host seasonality. No obvious dilution effect on local population structure was evident by the releases. However, the population was largely replaced by genetically divergent isolates once per decade but evolved with a pattern of balancing selection and towards expansion through adaptation, non-random outcrossing and isolate migration. This study not only unveils the real-time features of entomopathogenic fungal population genomics and evolution but also provides added values to alleviate the concerns of environmental safety regarding the biocontrol application of mycoinsecticides.

Characterization of Natural Isolates of Bipolaris maydis Associated with Mating Types, Genetic Diversity, and Pathogenicity in Fujian Province, China

Due to the natural destructiveness and persistence of the southern corn leaf blight (SCLB) fungus Bipolaris maydis (Nisikado et Miyake) Shoem, the characterization of B. maydis field isolates is essential to guide the rational distribution of resistant materials in corn-growing regions. In the present study, 102 field isolates collected from seven locations covering the entire region of Fujian Province, China, were assessed for mating type distribution, genetic diversity, and pathogenicity toward local sweet corn cultivars. Mating type detection via polymerase chain reaction indicated that 36.3 and 63.7% of isolates were MAT1-1 and MAT1-2, respectively; more than 80% of these isolates were confirmed using cross assays with known mating type isolates. Thirteen intersimple sequence repeat (ISSR) markers within and among two mating type populations revealed a high level of DNA polymorphism for all combined isolates and between MAT1-1 and MAT1-2 populations. The MAT1-2 population was more diverse based on DNA polymorphism than the MAT1-1 population. The value of G_ST was 0.0070, ranging from 0.0399 to 0.3044 based on analysis of combined isolates and individual regional populations, respectively, suggesting the presence of genetic differentiation in the two mating type populations from different locations. Pathogenicity assays revealed that both MAT1-1 and MAT1-2 populations were pathogenic to all 11 local sweet corn cultivars tested in this study. The potential of sexual reproduction, existence of genetic diversity in the two mating type populations, and pathogenicity suggest that B. maydis populations have independently clonally adapted under natural field conditions during corn cultivation.

Invasion and Extinction Dynamics of Mating Types Under Facultative Sexual Reproduction

In sexually reproducing isogamous species, syngamy between gametes is generally not indiscriminate, but rather restricted to occurring between complementary self-incompatible mating types. A longstanding question regards the evolutionary pressures that control the number of mating types observed in natural populations, which ranges from two to many thousands. Here, we describe a population genetic null model of this reproductive system, and derive expressions for the stationary probability distribution of the number of mating types, the establishment probability of a newly arising mating type, and the mean time to extinction of a resident type. Our results yield that the average rate of sexual reproduction in a population correlates positively with the expected number of mating types observed. We further show that the low number of mating types predicted in the rare-sex regime is primarily driven by low invasion probabilities of new mating type alleles, with established resident alleles being very stable over long evolutionary periods. Moreover, our model naturally exhibits varying selection strength dependent on the number of resident mating types. This results in higher extinction and lower invasion rates for an increasing number of residents.

Genetic and Genomic Analyses Reveal Boundaries between Species Closely Related to Cryptococcus Pathogens

Speciation is a central mechanism of biological diversification. While speciation is well studied in plants and animals, in comparison, relatively little is known about speciation in fungi. One fungal model is the Cryptococcus genus, which is best known for the pathogenic Cryptococcus neoformans/Cryptococcus gattii species complex that causes >200,000 new human infections annually. Elucidation of how these species evolved into important human-pathogenic species remains challenging and can be advanced by studying the most closely related nonpathogenic species, Cryptococcus amylolentus and Tsuchiyaea wingfieldii However, these species have only four known isolates, and available data were insufficient to determine species boundaries within this group. By analyzing full-length chromosome assemblies, we reappraised the phylogenetic relationships of the four available strains, confirmed the genetic separation of C. amylolentus and T. wingfieldii (now Cryptococcus wingfieldii), and revealed an additional cryptic species, for which the name Cryptococcus floricola is proposed. The genomes of the three species are ∼6% divergent and exhibit significant chromosomal rearrangements, including inversions and a reciprocal translocation that involved intercentromeric ectopic recombination, which together likely impose significant barriers to genetic exchange. Using genetic crosses, we show that while C. wingfieldii cannot interbreed with any of the other strains, C. floricola can still undergo sexual reproduction with C. amylolentus However, most of the resulting spores were inviable or sterile or showed reduced recombination during meiosis, indicating that intrinsic postzygotic barriers had been established. Our study and genomic data will foster additional studies addressing fungal speciation and transitions between nonpathogenic and pathogenic Cryptococcus lineages.IMPORTANCE The evolutionary drivers of speciation are critical to our understanding of how new pathogens arise from nonpathogenic lineages and adapt to new environments. Here we focus on the Cryptococcus amylolentus species complex, a nonpathogenic fungal lineage closely related to the human-pathogenic Cryptococcus neoformans/Cryptococcus gattii complex. Using genetic and genomic analyses, we reexamined the species boundaries of four available isolates within the C. amylolentus complex and revealed three genetically isolated species. Their genomes are ∼6% divergent and exhibit chromosome rearrangements, including translocations and small-scale inversions. Although two of the species (C. amylolentus and newly described C. floricola) were still able to interbreed, the resulting hybrid progeny were usually inviable or sterile, indicating that barriers to reproduction had already been established. These results advance our understanding of speciation in fungi and highlight the power of genomics in assisting our ability to correctly identify and discriminate fungal species.

Estimating the number of sexual events per generation in a facultatively sexual haploid population

In populations of facultatively sexual organisms, the proportion of sexually produced offspring contributed to each generation is a critical determinant of their evolutionary potential. However, estimating this parameter in natural populations has proved difficult. Here we develop a population genetic model for estimating the number of sexual events occurring per generation for facultatively sexual haploids possessing a biallelic mating-type locus (e.g., Chlamydomonas, ascomycete fungi). Our model treats the population as two subpopulations possessing opposite mating-type alleles, which exchange genes only when a sexual event takes place. Where mating types are equally abundant, we show that, for a neutral genetic marker, genetic differentiation between mating-type subpopulations is a simple function of the effective population size, the frequency of sexual reproduction, and the recombination fraction between the genetic marker and the mating-type locus. We employ simulations to examine the effects of linkage of markers to the mating-type locus, inequality of mating-type frequencies, mutation rate, and selection on this relationship. Finally, we apply our model to estimate the number of sexual reproduction events per generation in populations of four species of facultatively sexual ascomycete fungi, which have been jointly scored for mating type and a range of polymorphic molecular markers. Relative estimates are in line with expectations based on the known reproductive biology of these species.

Evidence for Sexual Reproduction: Identification, Frequency, and Spatial Distribution of Venturia effusa (Pecan Scab) Mating Type Idiomorphs

Venturia effusa (syn. Fusicladium effusum), causal agent of pecan scab, is the most prevalent pathogen of pecan (Carya illinoinensis), causing severe yield losses in the southeastern United States. V. effusa is currently known only by its asexual (conidial) stage. However, the degree and distribution of genetic diversity observed within and among populations of V. effusa are typical of a sexually reproducing fungal pathogen, and comparable with other dothideomycetes with a known sexual stage, including the closely related apple scab pathogen, V. inaequalis. Using the mating type (MAT) idiomorphs from V. inaequalis, we identified a single MAT gene, MAT1-1-1, in a draft genome of V. effusa. The MAT1-1-1 locus is flanked by two conserved genes encoding a DNA lyase (APN2) and a hypothetical protein. The MAT locus spanning the flanking genes was amplified and sequenced from a subset of 14 isolates, of which 7 contained MAT1-1-1 and the remaining samples contained MAT1-2-1. A multiplex polymerase chain reaction screen was developed to amplify MAT1-1-1, MAT1-2-1, and a conserved reference gene encoding β-tubulin, and used to screen 784 monoconidial isolates of V. effusa collected from 11 populations of pecan across the southeastern United States. A hierarchical sampling protocol representing region, orchard, and tree allowed for analysis of MAT structure at different spatial scales. Analysis of this collection revealed the frequency of the MAT idiomorphs is in a 1:1 equilibrium of MAT1-1:MAT1-2. The apparent equilibrium of the MAT idiomorphs provides impetus for a renewed effort to search for the sexual stage of V. effusa. [Formula: see text] Copyright © 2018 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license .

Evolutionary Pathways for the Generation of New Self-Incompatibility Haplotypes in a Nonself-Recognition System

Self-incompatibility (SI) is a genetically based recognition system that functions to prevent self-fertilization and mating among related plants. An enduring puzzle in SI is how the high diversity observed in nature arises and is maintained. Based on the underlying recognition mechanism, SI can be classified into two main groups: self-recognition (SR) and nonself-recognition (NSR). Most work has focused on diversification within SR systems despite expected differences between the two groups in the evolutionary pathways and outcomes of diversification. Here, we use a deterministic population genetic model and stochastic simulations to investigate how novel S-haplotypes evolve in a gametophytic NSR [SRNase/S Locus F-box (SLF)] SI system. For this model, the pathways for diversification involve either the maintenance or breakdown of SI and can vary in the order of mutations of the female (SRNase) and male (SLF) components. We show analytically that diversification can occur with high inbreeding depression and self-pollination, but this varies with evolutionary pathway and level of completeness (which determines the number of potential mating partners in the population), and, in general, is more likely for lower haplotype number. The conditions for diversification are broader in stochastic simulations of finite population size. However, the number of haplotypes observed under high inbreeding and moderate-to-high self-pollination is less than that commonly observed in nature. Diversification was observed through pathways that maintain SI as well as through self-compatible intermediates. Yet the lifespan of diversified haplotypes was sensitive to their level of completeness. By examining diversification in a NSR SI system, this model extends our understanding of the evolution and maintenance of haplotype diversity observed in a recognition system common in flowering plants.

Negative Frequency-Dependent Selection Is Frequently Confounding

Persistent genetic variation within populations presents an evolutionary problem, as natural selection and genetic drift tend to erode genetic diversity. Models of balancing selection were developed to account for the maintenance of genetic variation observed in natural populations. Negative frequency-dependent selection is a powerful type of balancing selection that maintains many natural polymorphisms, but it is also commonly misinterpreted. This review aims to clarify the processes underlying negative frequency-dependent selection, describe classes of polymorphisms that can and cannot result from these processes, and discuss the empirical data needed to accurately identify processes that generate or maintain diversity in nature. Finally, the importance of accurately describing the processes affecting genetic diversity within populations as it relates to research progress is considered.

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.

Global population structure and adaptive evolution of aflatoxin-producing fungi

Aflatoxins produced by several species in Aspergillus section Flavi are a significant problem in agriculture and a continuous threat to human health. To provide insights into the biology and global population structure of species in section Flavi, a total of 1,304 isolates were sampled across six species (A. flavus, A. parasiticus, A. nomius, A. caelatus, A. tamarii, and A. alliaceus) from single fields in major peanut‐growing regions in Georgia (USA), Australia, Argentina, India, and Benin (Africa). We inferred maximum‐likelihood phylogenies for six loci, both combined and separately, including two aflatoxin cluster regions (aflM/alfN and aflW/aflX) and four noncluster regions (amdS, trpC, mfs and MAT), to examine population structure and history. We also employed principal component and STRUCTURE analysis to identify genetic clusters and their associations with six different categories (geography, species, precipitation, temperature, aflatoxin chemotype profile, and mating type). Overall, seven distinct genetic clusters were inferred, some of which were more strongly structured by G chemotype diversity than geography. Populations of A. flavus S in Benin were genetically distinct from all other section Flavi species for the loci examined, which suggests genetic isolation. Evidence of trans‐speciation within two noncluster regions, whereby A. flavus S_BG strains from Australia share haplotypes with either A. flavus or A. parasiticus, was observed. Finally, while clay soil and precipitation may influence species richness in Aspergillus section Flavi, other region‐specific environmental and genetic parameters must also be considered.

Heterogeneous gene duplications can be adaptive because they permanently associate overdominant alleles

Gene duplications are widespread in genomes, but their role in contemporary adaptation is not fully understood. Although mostly deleterious, homogeneous duplications that associate identical repeats of a locus often increase the quantity of protein produced, which can be selected in certain environments. However, another type exists: heterogeneous gene duplications, which permanently associate two (or more) alleles of a single locus on the same chromosome. They are far less studied, as only few examples of contemporary heterogeneous duplications are known. Haldane proposed in 1954 that they could be adaptive in situations of heterozygote advantage, or overdominance, but this hypothesis was never tested. To assess its validity, we took advantage of the well-known model of insecticide resistance in mosquitoes. We used experimental evolution to estimate the fitnesses associated with homozygous and heterozygous genotypes in different selection regimes. It first showed that balanced antagonist selective pressures frequently induce overdominance, generating stable polymorphic equilibriums. The frequency of equilibrium moreover depends on the magnitude of two antagonistic selective pressures, the survival advantage conferred by the resistant allele versus the selective costs it induces. We then showed that heterogeneous duplications are selected over single-copy alleles in such contexts. They allow the fixation of the heterozygote phenotype, providing an alternative and stable intermediate fitness trade-off. By allowing the rapid fixation of divergent alleles, this immediate advantage could contribute to the rarity of overdominance. More importantly, it also creates new material for long-term genetic innovation, making a crucial but underestimated contribution to the evolution of new genes and gene families.

Comparative Genomics of the Ectomycorrhizal Sister Species Rhizopogon vinicolor and Rhizopogon vesiculosus (Basidiomycota: Boletales) Reveals a Divergence of the Mating Type B Locus

Divergence of breeding system plays an important role in fungal speciation. Ectomycorrhizal fungi, however, pose a challenge for the study of reproductive biology because most cannot be mated under laboratory conditions. To overcome this barrier, we sequenced the draft genomes of the ectomycorrhizal sister species Rhizopogon vinicolor Smith and Zeller and R. vesiculosus Smith and Zeller (Basidiomycota, Boletales)-the first genomes available for Basidiomycota truffles-and characterized gene content and organization surrounding their mating type loci. Both species possess a pair of homeodomain transcription factor homologs at the mating type A-locus as well as pheromone receptor and pheromone precursor homologs at the mating type B-locus. Comparison of Rhizopogon genomes with genomes from Boletales, Agaricales, and Polyporales revealed synteny of the A-locus region within Boletales, but several genomic rearrangements across orders. Our findings suggest correlation between gene content at the B-locus region and breeding system in Boletales with tetrapolar species possessing more diverse gene content than bipolar species. Rhizopogon vinicolor possesses a greater number of B-locus pheromone receptor and precursor genes than R. vesiculosus, as well as a pair of isoprenyl cysteine methyltransferase genes flanking the B-locus compared to a single copy in R. vesiculosus Examination of dikaryotic single nucleotide polymorphisms within genomes revealed greater heterozygosity in R. vinicolor, consistent with increased rates of outcrossing. Both species possess the components of a heterothallic breeding system with R. vinicolor possessing a B-locus region structure consistent with tetrapolar Boletales and R. vesiculosus possessing a B-locus region structure intermediate between bipolar and tetrapolar Boletales.

Unholy marriages and eternal triangles: how competition in the mushroom life cycle can lead to genomic conflict

In the vast majority of sexual life cycles, fusion between single-celled gametes is directly followed by nuclear fusion, leading to a diploid zygote and a lifelong commitment between two haploid genomes. Mushroom-forming basidiomycetes differ in two key respects. First, the multicellular haploid mating partners are fertilized in their entirety, each cell being a gamete that simultaneously can behave as a female, i.e. contributing the cytoplasm to a zygote by accepting nuclei, and a male gamete, i.e. only donating nuclei to the zygote. Second, after gamete union, the two haploid genomes remain separate so that the main vegetative stage, the dikaryon, has two haploid nuclei per cell. Only when the dikaryon produces mushrooms, do the nuclei fuse to enter a short diploid stage, immediately followed by meiosis and haploid spore formation. So in basidiomycetes, gamete fusion and genome mixing (sex) are separated in time. The 'living apart together' of nuclei in the dikaryon maintains some autonomy for nuclei to engage in a relationship with a different nucleus. We show that competition among the two nuclei of the dikaryon for such 'extramarital affairs' may lead to genomic conflict by favouring genes beneficial at the level of the nucleus, but deleterious at that of the dikaryon.This article is part of the themed issue 'Weird sex: the underappreciated diversity of sexual reproduction'.

Fruiting Body Formation in Volvariella volvacea Can Occur Independently of Its MAT-A-Controlled Bipolar Mating System, Enabling Homothallic and Heterothallic Life Cycles

Volvariella volvacea is an important crop in Southeast Asia, but erratic fruiting presents a serious challenge for its production and breeding. Efforts to explain inconsistent fruiting have been complicated by the multinucleate nature, typical lack of clamp connections, and an incompletely identified sexual reproductive system. In this study, we addressed the life cycle of V. volvacea using whole genome sequencing, cloning of MAT loci, karyotyping of spores, and fruiting assays. Microscopy analysis of spores had previously indicated the possible coexistence of heterothallic and homothallic life cycles. Our analysis of the MAT loci showed that only MAT-A, and not MAT-B, controlled heterokaryotization. Thus, the heterothallic life cycle was bipolar. Karyotyping of single spore isolates (SSIs) using molecular markers supported the existence of heterokaryotic spores. However, most SSIs were clearly not heterokaryotic, yet contained structural variation (SV) markers relating to both alleles of both parents. Heterokaryons from crossed, self-sterile homokaryons could produce fruiting bodies, agreeing with bipolar heterothallism. Meanwhile, some SSIs with two different MAT-A loci also produced fruiting bodies, which supported secondary homothallism. Next, SSIs that clearly contained only one MAT-A locus (homothallism) were also able to fruit, demonstrating that self-fertile SSIs were not, per definition, secondary homothallic, and that a third life cycle or genetic mechanism must exist. Finally, recombination between SV markers was normal, yet 10 out of 24 SV markers showed 1:2 or 1:3 distributions in the spores, and large numbers of SSIs contained doubled SV markers. This indicated selfish genes, and possibly partial aneuploidy.

Introgression maintains the genetic integrity of the mating-type determining chromosome of the fungus Neurospora tetrasperma

Genome evolution is driven by a complex interplay of factors, including selection, recombination, and introgression. The regions determining sexual identity are particularly dynamic parts of eukaryotic genomes that are prone to molecular degeneration associated with suppressed recombination. In the fungus Neurospora tetrasperma, it has been proposed that this molecular degeneration is counteracted by the introgression of nondegenerated DNA from closely related species. In this study, we used comparative and population genomic analyses of 92 genomes from eight phylogenetically and reproductively isolated lineages of N. tetrasperma, and its three closest relatives, to investigate the factors shaping the evolutionary history of the genomes.We found that suppressed recombination extends across at least 6 Mbp (∼ 63%) of the mating-type (mat) chromosome in N. tetrasperma and is associated with decreased genetic diversity, which is likely the result primarily of selection at linked sites. Furthermore, analyses of molecular evolution revealed an increased mutational load in this region, relative to recombining regions. However, comparative genomic and phylogenetic analyses indicate that the mat chromosomes are temporarily regenerated via introgression from sister species; six of eight lineages show introgression into one of their mat chromosomes, with multiple Neurospora species acting as donors. The introgressed tracts have been fixed within lineages, suggesting that they confer an adaptive advantage in natural populations, and our analyses support the presence of selective sweeps in at least one lineage. Thus, these data strongly support the previously hypothesized role of introgression as a mechanism for the maintenance of mating-type determining chromosomal regions.

Population genetic analysis reveals cryptic sex in the phytopathogenic fungus Alternaria alternata

Reproductive mode can impact population genetic dynamics and evolutionary landscape of plant pathogens as well as on disease epidemiology and management. In this study, we monitored the spatial dynamics and mating type idiomorphs in ~700 Alternaria alternata isolates sampled from the main potato production areas in China to infer the mating system of potato early blight. Consistent with the expectation of asexual species, identical genotypes were recovered from different locations separated by hundreds of kilometers of geographic distance and spanned across many years. However, high genotype diversity, equal MAT1-1 and MAT1-2 frequencies within and among populations, no genetic differentiation and phylogenetic association between two mating types, combined with random association amongst neutral markers in some field populations, suggested that sexual reproduction may also play an important role in the epidemics and evolution of the pathogen in at least half of the populations assayed despite the fact that no teleomorphs have been observed yet naturally or artificially. Our results indicated that A. alternata may adopt an epidemic mode of reproduction by combining many cycles of asexual propagation with fewer cycles of sexual reproduction, facilitating its adaptation to changing environments and making the disease management on potato fields even more difficult.

Variation in mate-recognition pheromones of the fungal genus Microbotryum

Mate recognition is an essential life-cycle stage that exhibits strong conservation in function, whereas diversification of mating signals can contribute directly to the integrity of species boundaries through assortative mating. Fungi are simple models, where compatibility is based on the recognition of pheromone peptides by corresponding receptor proteins, but clear patterns of diversification have not emerged from the species examined, which are few compared with mate signaling studies in plant and animal systems. In this study, candidate loci from Microbotryum species were used to characterize putative pheromones that were synthesized and found to be functional across multiple species in triggering a mating response in vitro. There is no significant correlation between the strength of a species' response and its genetic distance from the pheromone sequence source genome. Instead, evidence suggests that species may be strong or weak responders, influenced by environmental conditions or developmental differences. Gene sequence comparisons reveals very strong purifying selection on the a₁ pheromone peptide and corresponding receptor, but significantly less purifying selection on the a₂ pheromone peptide that corresponds with more variation across species in the receptor. This represents an exceptional case of a reciprocally interacting mate-recognition system in which the two mating types are under different levels of purifying selection.

Diversification and distinctive structural features of S-RNase alleles in the genus Solanum

The multigenic and multiallelic S-locus in plants is responsible for the gametophytic self-incompatibility system, which is important to prevent the detrimental effects of self-fertilization and inbreeding depression. Several studies have discussed the importance of punctual mutations, recombination, and natural selection in the generation of allelic diversity in the S-locus. However, there has been no wide-ranging study correlating the molecular evolution and structural aspects of the corresponding proteins in Solanum. Therefore, we evaluated the molecular evolution of one gene in this locus and generated a statistically well-supported phylogenetic tree, as well as evidence of positive selection, helping us to understand the diversification of S alleles in Solanum. The three-dimensional structures of some of the proteins corresponding to the major clusters of the phylogenetic tree were constructed and subsequently submitted to molecular dynamics to stabilize the folding and obtain the native structure. The positively selected amino acid residues were predominantly located in the hyper variable regions and on the surface of the protein, which appears to be fundamental for allele specificity. One of the positively selected residues was identified adjacent to a conserved strand that is crucial for enzymatic catalysis. Additionally, we have shown significant differences in the electrostatic potential among the predicted molecular surfaces in S-RNases. The structural results indicate that local changes in the three-dimensional structure are present in some regions of the molecule, although the general structure seems to be conserved. No previous study has described such structural variations in S-RNases.

Nucleotide variability at its limit? Insights into the number and evolutionary dynamics of the sex-determining specificities of the honey bee Apis mellifera

Deciphering the evolutionary processes driving nucleotide variation in multiallelic genes is limited by the number of genetic systems in which such genes occur. The complementary sex determiner (csd) gene in the honey bee Apis mellifera is an informative example for studying allelic diversity and the underlying evolutionary forces in a well-described model of balancing selection. Acting as the primary signal of sex determination, diploid individuals heterozygous for csd develop into females, whereas csd homozygotes are diploid males that have zero fitness. Examining 77 of the functional heterozygous csd allele pairs, we established a combinatorical criteria that provide insights into the minimum number of amino acid differences among those pairs. Given a data set of 244 csd sequences, we show that the total number of csd alleles found in A. mellifera ranges from 53 (locally) to 87 (worldwide), which is much higher than was previously reported (20). Using a coupon-collector model, we extrapolate the presence of in total 116–145 csd alleles worldwide. The hypervariable region (HVR) is of particular importance in determining csd allele specificity, and we provide for this region evidence of high evolutionary rate for length differences exceeding those of microsatellites. The proportion of amino acids driven by positive selection and the rate of nonsynonymous substitutions in the HVR-flanking regions reach values close to 1 but differ with respect to the HVR length. Using a model of csd coalescence, we identified the high originating rate of csd specificities as a major evolutionary force, leading to an origin of a novel csd allele every 400,000 years. The csd polymorphism frequencies in natural populations indicate an excess of new mutations, whereas signs of ancestral transspecies polymorphism can still be detected. This study provides a comprehensive view of the enormous diversity and the evolutionary forces shaping a multiallelic gene.

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.

Evolution of uni- and bifactorial sexual compatibility systems in fungi

Mating systems, that is, whether organisms give rise to progeny by selfing, inbreeding or outcrossing, strongly affect important ecological and evolutionary processes. Large variations in mating systems exist in fungi, allowing the study of their origin and consequences. In fungi, sexual incompatibility is determined by molecular recognition mechanisms, controlled by a single mating-type locus in most unifactorial fungi. In Basidiomycete fungi, however, which include rusts, smuts and mushrooms, a system has evolved in which incompatibility is controlled by two unlinked loci. This bifactorial system probably evolved from a unifactorial system. Multiple independent transitions back to a unifactorial system occurred. It is still unclear what force drove evolution and maintenance of these contrasting inheritance patterns that determine mating compatibility. Here, we give an overview of the evolutionary factors that might have driven the evolution of bifactoriality from a unifactorial system and the transitions back to unifactoriality. Bifactoriality most likely evolved for selfing avoidance. Subsequently, multiallelism at mating-type loci evolved through negative frequency-dependent selection by increasing the chance to find a compatible mate. Unifactoriality then evolved back in some species, possibly because either selfing was favoured or for increasing the chance to find a compatible mate in species with few alleles. Owing to the existence of closely related unifactorial and bifactorial species and the increasing knowledge of the genetic systems of the different mechanisms, Basidiomycetes provide an excellent model for studying the different forces that shape breeding systems.

Mutualism and asexual reproduction influence recognition genes in a fungal symbiont

Mutualism between microbes and insects is common and alignment of the reproductive interests of microbial symbionts with this lifestyle typically involves clonal reproduction and vertical transmission by insect partners. Here the Amylostereum fungus-Sirex woodwasp mutualism was used to consider whether their prolonged association and predominance of asexuality have affected the mating system of the fungal partner. Nucleotide information for the pheromone receptor gene rab1, as well as the translation elongation factor 1α gene and ribosomal RNA internal transcribed spacer region were utilized. The identification of rab1 alleles in Amylostereum chailletii and Amylostereum areolatum populations revealed that this gene is more polymorphic than the other two regions, although the diversity of all three regions was lower than what has been observed in free-living Agaricomycetes. Our data suggest that suppressed recombination might be implicated in the diversification of rab1, while no evidence of balancing selection was detected. We also detected positive selection at only two codons, suggesting that purifying selection is important for the evolution of rab1. The symbiotic relationship with their insect partners has therefore influenced the diversity of this gene and influenced the manner in which selection drives and maintains this diversity in A. areolatum and A. chailletii.

Sexual selection in fungi

The significance of sexual selection, the component of natural selection associated with variation in mating success, is well established for the evolution of animals and plants, but not for the evolution of fungi. Even though fungi do not have separate sexes, most filamentous fungi mate in a hermaphroditic fashion, with distinct sex roles, that is, investment in large gametes (female role) and fertilization by other small gametes (male role). Fungi compete to fertilize, analogous to 'male-male' competition, whereas they can be selective when being fertilized, analogous to female choice. Mating types, which determine genetic compatibility among fungal gametes, are important for sexual selection in two respects. First, genes at the mating-type loci regulate different aspects of mating and thus can be subject to sexual selection. Second, for sexual selection, not only the two sexes (or sex roles) but also the mating types can form the classes, the members of which compete for access to members of the other class. This is significant if mating-type gene products are costly, thus signalling genetic quality according to Zahavi's handicap principle. We propose that sexual selection explains various fungal characteristics such as the observed high redundancy of pheromones at the B mating-type locus of Agaricomycotina, the occurrence of multiple types of spores in Ascomycotina or the strong pheromone signalling in yeasts. Furthermore, we argue that fungi are good model systems to experimentally study fundamental aspects of sexual selection, due to their fast generation times and high diversity of life cycles and mating 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.

Insight into trade-off between wood decay and parasitism from the genome of a fungal forest pathogen

Parasitism and saprotrophic wood decay are two fungal strategies fundamental for succession and nutrient cycling in forest ecosystems. An opportunity to assess the trade-off between these strategies is provided by the forest pathogen and wood decayer Heterobasidion annosum sensu lato. We report the annotated genome sequence and transcript profiling, as well as the quantitative trait loci mapping, of one member of the species complex: H. irregulare. Quantitative trait loci critical for pathogenicity, and rich in transposable elements, orphan and secreted genes, were identified. A wide range of cellulose-degrading enzymes are expressed during wood decay. By contrast, pathogenic interaction between H. irregulare and pine engages fewer carbohydrate-active enzymes, but involves an increase in pectinolytic enzymes, transcription modules for oxidative stress and secondary metabolite production. Our results show a trade-off in terms of constrained carbohydrate decomposition and membrane transport capacity during interaction with living hosts. Our findings establish that saprotrophic wood decay and necrotrophic parasitism involve two distinct, yet overlapping, processes.

Evaluation of the potential for sexual reproduction in field populations of Cercospora beticola from USA

Cercospora leaf spot, caused by the hemibiotrophic fungal pathogen Cercospora beticola, is the most economically damaging foliar disease of sugarbeet worldwide. Although most C. beticola populations display characteristics reminiscent of sexual recombination, no teleomorph has been described. To assess whether populations in northern United States have characteristics consistent with sexual reproduction, 1024 isolates collected over a 3-y period were analyzed for frequency and distribution of mating type genes. After clone correction, an approximately equal distribution of mating types was found for each sampling year. Mating type frequency was also assessed in individual lesions. Lesions always consisted of isolates with a single mating type and microsatellite haplotype, but both mating types and up to five microsatellite haplotypes could be found on an individual leaf. The MAT1-1-1 and MAT1-2-1 genes were sequenced from 28 MAT1-1 and 28 MAT1-2 isolates, respectively. Three MAT1-1-1 nucleotide haplotypes were identified that encoded a single amino acid sequence. For MAT1-2-1, five nucleotide haplotypes were identified that encoded four protein variants. MAT1-1-1 and MAT1-2-1 gene expression analyses were conducted on plants inoculated with either or both mating types. MAT1-1-1 expression remained low, but MAT1-2-1 spiked during late stages of colonization. A segment of the MAT1-2-1 coding sequence was also found in MAT1-1 isolates. Taken together, these results suggest that C. beticola has the potential for sexual reproduction.

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.

Interspecific sex in grass smuts and the genetic diversity of their pheromone-receptor system

The grass smuts comprise a speciose group of biotrophic plant parasites, so-called Ustilaginaceae, which are specifically adapted to hosts of sweet grasses, the Poaceae family. Mating takes a central role in their life cycle, as it initiates parasitism by a morphological and physiological transition from saprobic yeast cells to pathogenic filaments. As in other fungi, sexual identity is determined by specific genomic regions encoding allelic variants of a pheromone-receptor (PR) system and heterodimerising transcription factors. Both operate in a biphasic mating process that starts with PR-triggered recognition, directed growth of conjugation hyphae, and plasmogamy of compatible mating partners. So far, studies on the PR system of grass smuts revealed diverse interspecific compatibility and mating type determination. However, many questions concerning the specificity and evolutionary origin of the PR system remain unanswered. Combining comparative genetics and biological approaches, we report on the specificity of the PR system and its genetic diversity in 10 species spanning about 100 million years of mating type evolution. We show that three highly syntenic PR alleles are prevalent among members of the Ustilaginaceae, favouring a triallelic determination as the plesiomorphic characteristic of this group. Furthermore, the analysis of PR loci revealed increased genetic diversity of single PR locus genes compared to genes of flanking regions. Performing interspecies sex tests, we detected a high potential for hybridisation that is directly linked to pheromone signalling as known from intraspecies sex. Although the PR system seems to be optimised for intraspecific compatibility, the observed functional plasticity of the PR system increases the potential for interspecific sex, which might allow the hybrid-based genesis of newly combined host specificities.

Evolution of the S-locus region in Arabidopsis relatives

The S locus, a single polymorphic locus, is responsible for self-incompatibility (SI) in the Brassicaceae family and many related plant families. Despite its importance, our knowledge of S-locus evolution is largely restricted to the causal genes encoding the S-locus receptor kinase (SRK) receptor and S-locus cysteine-rich protein (SCR) ligand of the SI system. Here, we present high-quality sequences of the genomic region of six S-locus haplotypes: Arabidopsis (Arabidopsis thaliana; one haplotype), Arabidopsis lyrata (four haplotypes), and Capsella rubella (one haplotype). We compared these with reference S-locus haplotypes of the self-compatible Arabidopsis and its SI congener A. lyrata. We subsequently reconstructed the likely genomic organization of the S locus in the most recent common ancestor of Arabidopsis and Capsella. As previously reported, the two SI-determining genes, SCR and SRK, showed a pattern of coevolution. In addition, consistent with previous studies, we found that duplication, gene conversion, and positive selection have been important factors in the evolution of these two genes and appear to contribute to the generation of new recognition specificities. Intriguingly, the inactive pseudo-S-locus haplotype in the self-compatible species C. rubella is likely to be an old S-locus haplotype that only very recently became fixed when C. rubella split off from its SI ancestor, Capsella grandiflora.

Comparative genomics of the mating-type loci of the mushroom Flammulina velutipes reveals widespread synteny and recent inversions

BACKGROUND

Mating-type loci of mushroom fungi contain master regulatory genes that control recognition between compatible nuclei, maintenance of compatible nuclei as heterokaryons, and fruiting body development. Regions near mating-type loci in fungi often show adapted recombination, facilitating the generation of novel mating types and reducing the production of self-compatible mating types. Compared to other fungi, mushroom fungi have complex mating-type systems, showing both loci with redundant function (subloci) and subloci with many alleles. The genomic organization of mating-type loci has been solved in very few mushroom species, which complicates proper interpretation of mating-type evolution and use of those genes in breeding programs.

METHODOLOGY/PRINCIPAL FINDINGS

We report a complete genetic structure of the mating-type loci from the tetrapolar, edible mushroom Flammulina velutipes mating type A3B3. Two matB3 subloci, matB3a that contains a unique pheromone and matB3b, were mapped 177 Kb apart on scaffold 1. The matA locus of F. velutipes contains three homeodomain genes distributed over 73 Kb distant matA3a and matA3b subloci. The conserved matA region in Agaricales approaches 350 Kb and contains conserved recombination hotspots showing major rearrangements in F. velutipes and Schizophyllum commune. Important evolutionary differences were indicated; separation of the matA subloci in F. velutipes was diverged from the Coprinopsis cinerea arrangement via two large inversions whereas separation in S. commune emerged through transposition of gene clusters.

CONCLUSIONS/SIGNIFICANCE

In our study we determined that the Agaricales have very large scale synteny at matA (∼350 Kb) and that this synteny is maintained even when parts of this region are separated through chromosomal rearrangements. Four conserved recombination hotspots allow reshuffling of large fragments of this region. Next to this, it was revealed that large distance subloci can exist in matB as well. Finally, the genes that were linked to specific mating types will serve as molecular markers in breeding.

Origin and diversification dynamics of self-incompatibility haplotypes

Self-incompatibility (SI) is a genetic system found in some hermaphrodite plants. Recognition of pollen by pistils expressing cognate specificities at two linked genes leads to rejection of self pollen and pollen from close relatives, i.e., to avoidance of self-fertilization and inbred matings, and thus increased outcrossing. These genes generally have many alleles, yet the conditions allowing the evolution of new alleles remain mysterious. Evolutionary changes are clearly necessary in both genes, since any mutation affecting only one of them would result in a nonfunctional self-compatible haplotype. Here, we study diversification at the S-locus (i.e., a stable increase in the total number of SI haplotypes in the population, through the incorporation of new SI haplotypes), both deterministically (by investigating analytically the fate of mutations in an infinite population) and by simulations of finite populations. We show that the conditions allowing diversification are far less stringent in finite populations with recurrent mutations of the pollen and pistil genes, suggesting that diversification is possible in a panmictic population. We find that new SI haplotypes emerge fastest in populations with few SI haplotypes, and we discuss some implications for empirical data on S-alleles. However, allele numbers in our simulations never reach values as high as observed in plants whose SI systems have been studied, and we suggest extensions of our models that may reconcile the theory and data.