Recombination suppression and evolutionary strata around mating-type loci in fungi: documenting patterns and understanding evolutionary and mechanistic causes

Discovery of New Genomic Configuration of Mating-Type Loci in the Largest Lineage of Lichen-Forming Fungi

The genetic architecture of mating-type loci in lichen-forming fungi has been characterized in very few taxa. Despite the limited data, and in contrast to all other major fungal lineages, arrangements that have both mating-type alleles in a single haploid genome have been hypothesized to be absent from the largest lineage of lichen-forming fungi, the Lecanoromycetes. We report the discovery of both mating-type alleles from the haploid genomes of three species within this group. Our results demonstrate that Lecanoromycetes are not an outlier among Ascomycetes.

Why do sex chromosomes progressively lose recombination?

Progressive recombination loss is a common feature of sex chromosomes. Yet, the evolutionary drivers of this phenomenon remain a mystery. For decades, differences in trait optima between sexes (sexual antagonism) have been the favoured hypothesis, but convincing evidence is lacking. Recent years have seen a surge of alternative hypotheses to explain progressive extensions and maintenance of recombination suppression: neutral accumulation of sequence divergence, selection of nonrecombining fragments with fewer deleterious mutations than average, sheltering of recessive deleterious mutations by linkage to heterozygous alleles, early evolution of dosage compensation, and constraints on recombination restoration. Here, we explain these recent hypotheses and dissect their assumptions, mechanisms, and predictions. We also review empirical studies that have brought support to the various hypotheses.

Absence of chromosome axis protein recruitment prevents meiotic recombination chromosome-wide in the budding yeast Lachancea kluyveri

Meiotic recombination shows broad variations across species and along chromosomes and is often suppressed at and around genomic regions determining sexual compatibility such as mating type loci in fungi. Here, we show that the absence of Spo11-DSBs and meiotic recombination on Lakl0C-left, the chromosome arm containing the sex locus of the Lachancea kluyveri budding yeast, results from the absence of recruitment of the two chromosome axis proteins Red1 and Hop1, essential for proper Spo11-DSBs formation. Furthermore, cytological observation of spread pachytene meiotic chromosomes reveals that Lakl0C-left does not undergo synapsis. However, we show that the behavior of Lakl0C-left is independent of its particularly early replication timing and is not accompanied by any peculiar chromosome structure as detectable by Hi-C in this yet poorly studied yeast. Finally, we observed an accumulation of heterozygous mutations on Lakl0C-left and a sexual dimorphism of the haploid meiotic offspring, supporting a direct effect of this absence of meiotic recombination on L. kluyveri genome evolution and fitness. Because suppression of meiotic recombination on sex chromosomes is widely observed across eukaryotes, the mechanism for recombination suppression described here may apply to other species, with the potential to impact sex chromosome evolution.

Genome biology and evolution of mating-type loci in four cereal rust fungi

Permanent heterozygous loci, such as sex- or mating-compatibility regions, often display suppression of recombination and signals of genomic degeneration. In Basidiomycota, two distinct loci confer mating compatibility. These loci encode homeodomain (HD) transcription factors and pheromone receptor (Pra)-ligand allele pairs. To date, an analysis of genome level mating-type (MAT) loci is lacking for obligate biotrophic basidiomycetes in the Pucciniales, an order containing serious agricultural plant pathogens. Here, we focus on four species of Puccinia that infect oat and wheat, including P. coronata f. sp. avenae, P. graminis f. sp. tritici, P. triticina and P. striiformis f. sp. tritici. MAT loci are located on two separate chromosomes supporting previous hypotheses of a tetrapolar mating compatibility system in the Pucciniales. The HD genes are multiallelic in all four species while the PR locus appears biallelic, except for P. graminis f. sp. tritici, which potentially has multiple alleles. HD loci are largely conserved in their macrosynteny, both within and between species, without strong signals of recombination suppression. Regions proximal to the PR locus, however, displayed signs of recombination suppression and genomic degeneration in the three species with a biallelic PR locus. Our observations support a link between recombination suppression, genomic degeneration, and allele diversity of MAT loci that is consistent with recent mathematical modelling and simulations. Finally, we confirm that MAT genes are expressed during the asexual infection cycle, and we propose that this may support regulating nuclear maintenance and pairing during infection and spore formation. Our study provides insights into the evolution of MAT loci of key pathogenic Puccinia species. Understanding mating compatibility can help predict possible combinations of nuclear pairs, generated by sexual reproduction or somatic recombination, and the potential evolution of new virulent isolates of these important plant pathogens.

The interplay of local adaptation and gene flow may lead to the formation of supergenes

Supergenes are genetic architectures resulting in the segregation of alternative combinations of alleles underlying complex phenotypes. The co-segregation of alleles at linked loci is often facilitated by polymorphic chromosomal rearrangements suppressing recombination locally. Supergenes are involved in many complex polymorphisms, including sexual, colour or behavioural polymorphisms in numerous plants, fungi, mammals, fish, and insects. Despite a long history of empirical and theoretical research, the formation of supergenes remains poorly understood. Here, using a two-island population genetic model, we explore how gene flow and the evolution of overdominant chromosomal inversions may jointly lead to the formation of supergenes. We show that the evolution of inversions in differentiated populations, both under disruptive selection, leads to an increase in frequency of poorly adapted, immigrant haplotypes. Indeed, rare allelic combinations, such as immigrant haplotypes, are more frequently reshuffled by recombination than common allelic combinations, and therefore benefit from the recombination suppression generated by inversions. When an inversion capturing a locally adapted haplotype spreads but is associated with a fitness cost hampering its fixation (e.g. a recessive mutation load), the maintenance of a non-inverted haplotype in the population is enhanced; under certain conditions, the immigrant haplotype persists alongside the inverted local haplotype, while the standard local haplotype disappears. This establishes a stable, local polymorphism with two non-recombining haplotypes encoding alternative adaptive strategies, that is, a supergene. These results bring new light to the importance of local adaptation, overdominance, and gene flow in the formation of supergenes and inversion polymorphisms in general.

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.

Comparative genomics of Cryptococcus and Kwoniella reveals pathogenesis evolution and contrasting karyotype dynamics via intercentromeric recombination or chromosome fusion

A large-scale comparative genomic analysis was conducted for the global human fungal pathogens within the Cryptococcus genus, compared to non-pathogenic Cryptococcus species, and related species from the sister genus Kwoniella. Chromosome-level genome assemblies were generated for multiple species of both genera, resulting in a dataset encompassing virtually all of their known diversity. Although Cryptococcus and Kwoniella have comparable genome sizes (about 19.2 and 22.9 Mb) and similar gene content, hinting at pre-adaptive pathogenic potential, our analysis found evidence in pathogenic Cryptococcus species of specific examples of gene gain (via horizontal gene transfer) and gene loss, which might represent evolutionary signatures of pathogenic development. Genome analysis also revealed a significant variation in chromosome number and structure between the two genera. By combining synteny analysis and experimental centromere validation, we found that most Cryptococcus species have 14 chromosomes, whereas most Kwoniella species have fewer (11, 8, 5 or even as few as 3). Reduced chromosome number in Kwoniella is associated with formation of giant chromosomes (up to 18 Mb) through repeated chromosome fusion events, each marked by a pericentric inversion and centromere loss. While similar chromosome inversion-fusion patterns were observed in all Kwoniella species with fewer than 14 chromosomes, no such pattern was detected in Cryptococcus. Instead, Cryptococcus species with less than 14 chromosomes, underwent chromosome reductions primarily through rearrangements associated with the loss of repeat-rich centromeres. Additionally, Cryptococcus genomes exhibited frequent interchromosomal translocations, including intercentromeric recombination facilitated by transposons shared between centromeres. Taken together, our findings advance our understanding of genomic changes possibly associated with pathogenicity in Cryptococcus and provide a foundation to elucidate mechanisms of centromere loss and chromosome fusion driving distinct karyotypes in closely related fungal species, including prominent global human pathogens.

SM, Lee F, Lockhart SR, Luethy P, Martin I, Kwon-Chung KJ, Meyer W, Nguyen MH, Ostrosky-Zeichner L, Palavecino E, Pancholi P, Pappas PG, Procop GW, Redhead SA, Rhoads DD, Riedel S, Stevens B, Sullivan KO, Vergidis P, Roilides E, Seyedmousavi A, Tao L, Vicente VA, Vitale RG, Wang QM, Wengenack NL, Westblade L, Wiederhold N, White L, Wojewoda CM, Zhang SX. A conceptual framework for nomenclatural stability and validity of medically important fungi: a proposed global consensus guideline for fungal name changes supported by ABP, ASM, CLSI, ECMM, ESCMID-EFISG, EUCAST-AFST, FDLC, IDSA, ISHAM, MMSA, and MSGERC

The rapid pace of name changes of medically important fungi is creating challenges for clinical laboratories and clinicians involved in patient care. We describe two sources of name change which have different drivers, at the species versus the genus level. Some suggestions are made here to reduce the number of name changes. We urge taxonomists to provide diagnostic markers of taxonomic novelties. Given the instability of phylogenetic trees due to variable taxon sampling, we advocate to maintain genera at the largest possible size. Reporting of identified species in complexes or series should where possible comprise both the name of the overarching species and that of the molecular sibling, often cryptic species. Because the use of different names for the same species will be unavoidable for many years to come, an open access online database of the names of all medically important fungi, with proper nomenclatural designation and synonymy, is essential. We further recommend that while taxonomic discovery continues, the adaptation of new name changes by clinical laboratories and clinicians be reviewed routinely by a standing committee for validation and stability over time, with reference to an open access database, wherein reasons for changes are listed in a transparent way.

Dynamics of transposable element accumulation in the non-recombining regions of mating-type chromosomes in anther-smut fungi

In the absence of recombination, the number of transposable elements (TEs) increases due to less efficient selection, but the dynamics of such TE accumulations are not well characterized. Leveraging a dataset of 21 independent events of recombination cessation of different ages in mating-type chromosomes of Microbotryum fungi, we show that TEs rapidly accumulated in regions lacking recombination, but that TE content reached a plateau at ca. 50% of occupied base pairs by 1.5 million years following recombination suppression. The same TE superfamilies have expanded in independently evolved non-recombining regions, in particular rolling-circle replication elements (Helitrons). Long-terminal repeat (LTR) retrotransposons of the Copia and Ty3 superfamilies also expanded, through transposition bursts (distinguished from gene conversion based on LTR divergence), with both non-recombining regions and autosomes affected, suggesting that non-recombining regions constitute TE reservoirs. This study improves our knowledge of genome evolution by showing that TEs can accumulate through bursts, following non-linear decelerating dynamics.

The fate of sex chromosomes during the evolution of monoicy from dioicy in liverworts

Liverworts comprise one of six primary land plant lineages, with the predicted origin of extant liverwort diversity dating to the Silurian. The ancestral liverwort has been inferred to have been dioicous (unisexual) with chromosomal sex determination in which the U chromosome of females and the V chromosome of males were dimorphic with an extensive non-recombining region. In liverworts, sex is determined by a U chromosomal "feminizer" gene that promotes female development, and in its absence, male development ensues. Monoicy (bisexuality) has independently evolved multiple times within liverworts. Here, we explore the evolution of monoicy, focusing on the monoicous species Ricciocarpos natans, and propose that the evolution of monoicy in R. natans involved the appearance of an aneuploid spore that possessed both U and V chromosomes. Chromosomal rearrangements involving the U chromosome resulted in distribution of essential U chromosome genes, including the feminizer, to several autosomal locations. By contrast, we infer that the ancestral V chromosome was inherited largely intact, probably because it carries numerous dispersed "motility" genes distributed across the chromosome. The genetic networks for sex differentiation in R. natans appear largely unchanged except that the feminizer is developmentally regulated, allowing for temporally separated differentiation of female and male reproductive organs on a single plant. A survey of other monoicous liverworts suggests that similar genomic rearrangements may have occurred repeatedly in lineages transitioning to monoicy from dioicy. These data provide a foundation for understanding how genetic networks controlling sex determination can be subtly rewired to produce profound changes in sexual systems.

Frequent transitions in mating-type locus chromosomal organization in Malassezia and early steps in sexual reproduction

Fungi in the basidiomycete genus Malassezia are the most prevalent eukaryotic microbes resident on the skin of human and other warm-blooded animals and have been implicated in skin diseases and systemic disorders. Analysis of Malassezia genomes revealed that key adaptations to the skin microenvironment have a direct genomic basis, and the identification of mating/meiotic genes suggests a capacity to reproduce sexually, even though no sexual cycle has yet been observed. In contrast to other bipolar or tetrapolar basidiomycetes that have either two linked mating-type-determining (MAT) loci or two MAT loci on separate chromosomes, in Malassezia species studied thus far the two MAT loci are arranged in a pseudobipolar configuration (linked on the same chromosome but capable of recombining). By generating additional chromosome-level genome assemblies, and an improved Malassezia phylogeny, we infer that the pseudobipolar arrangement was the ancestral state of this group and revealed six independent transitions to tetrapolarity, seemingly driven by centromere fission or translocations in centromere-flanking regions. Additionally, in an approach to uncover a sexual cycle, Malassezia furfur strains were engineered to express different MAT alleles in the same cell. The resulting strains produce hyphae reminiscent of early steps in sexual development and display upregulation of genes associated with sexual development as well as others encoding lipases and a protease potentially relevant for pathogenesis of the fungus. Our study reveals a previously unseen genomic relocation of mating-type loci in fungi and provides insight toward the identification of a sexual cycle in Malassezia, with possible implications for pathogenicity.

Structure Analysis of the MatA Locus of Sexual Compatibility in the Edible Mushroom Pleurotus ostreatus

The edible oyster mushroom Pleurotus ostreatus is one of the most cultivated species worldwide. Morphogenesis associated with the maturation of fruit bodies is controlled by two unlinked loci of sexual compatibility matA and matB with multiple alleles (tetrapolar system of sexual compatibility). Quantitative analysis of the alleles of mating compatibility loci in 17 natural isolates collected in the Moscow region was performed in mon-mon (monokaryons-monokaryon) and di-mon (dikaryon-monokaryon) crossings. Four monokaryotic testers strains which were heteroallelic at both mating type loci were obtained for each of the five natural mushroom isolates by using original technique of sterile spore prints on Petri dishes and mon-mon crossing. Twelve natural isolates were crossed via di-mon mating with the four monokaryotic testers M-38. Genetic analysis of the alleles of sexual compatibility loci in 17 natural isolates revealed multiple alleles at both loci: at least ten alleles at matA locus and eight alleles at matB locus. Structural organization analysis of the matA locus was performed in silico for homokaryotic strains PC9 and PC15 based on the whole-genome sequencing data available at DOE Joint Genome Institute. The matA locus has an extremely divergent structure: there are one copy of the homeodomain gene hd1 and one copy of the hd2 gene in the PC9 strain, whereas the matA locus of the PC15 strain is composed by two copies of hd1.1 and hd1.2 genes (class HD1 homeodomain proteins) and one copy of hd2 gene (class HD2 proteins). Comprehensive analysis of amino acid sequences of HD1 and HD2 homeodomain proteins demonstrated that the proteins have a globular structure with the nuclear localization and contain a variable N-terminus and a more conserved DNA-binding domain with a specific conserved motif  WFXNXR in the third ɑ-helix. The results suggest that multiple alleles of the matA locus of sexual compatibility in basidiomycete fungi is achieved due to both different copy number of the coding hd genes within the locus and the variability of the coding gene sequences.

The Sordariomycetes: an expanding resource with Big Data for mining in evolutionary genomics and transcriptomics

Advances in genomics and transcriptomics accompanying the rapid accumulation of omics data have provided new tools that have transformed and expanded the traditional concepts of model fungi. Evolutionary genomics and transcriptomics have flourished with the use of classical and newer fungal models that facilitate the study of diverse topics encompassing fungal biology and development. Technological advances have also created the opportunity to obtain and mine large datasets. One such continuously growing dataset is that of the Sordariomycetes, which exhibit a richness of species, ecological diversity, economic importance, and a profound research history on amenable models. Currently, 3,574 species of this class have been sequenced, comprising nearly one-third of the available ascomycete genomes. Among these genomes, multiple representatives of the model genera Fusarium, Neurospora, and Trichoderma are present. In this review, we examine recently published studies and data on the Sordariomycetes that have contributed novel insights to the field of fungal evolution via integrative analyses of the genetic, pathogenic, and other biological characteristics of the fungi. Some of these studies applied ancestral state analysis of gene expression among divergent lineages to infer regulatory network models, identify key genetic elements in fungal sexual development, and investigate the regulation of conidial germination and secondary metabolism. Such multispecies investigations address challenges in the study of fungal evolutionary genomics derived from studies that are often based on limited model genomes and that primarily focus on the aspects of biology driven by knowledge drawn from a few model species. Rapidly accumulating information and expanding capabilities for systems biological analysis of Big Data are setting the stage for the expansion of the concept of model systems from unitary taxonomic species/genera to inclusive clusters of well-studied models that can facilitate both the in-depth study of specific lineages and also investigation of trait diversity across lineages. The Sordariomycetes class, in particular, offers abundant omics data and a large and active global research community. As such, the Sordariomycetes can form a core omics clade, providing a blueprint for the expansion of our knowledge of evolution at the genomic scale in the exciting era of Big Data and artificial intelligence, and serving as a reference for the future analysis of different taxonomic levels within the fungal kingdom.

Frequent transitions in mating-type locus chromosomal organization in Malassezia and early steps in sexual reproduction

Fungi in the basidiomycete genus Malassezia are the most prevalent eukaryotic microbes resident on the skin of human and other warm-blooded animals and have been implicated in skin diseases and systemic disorders. Analysis of Malassezia genomes revealed that key adaptations to the skin microenvironment have a direct genomic basis, and the identification of mating/meiotic genes suggests a capacity to reproduce sexually, even though no sexual cycle has yet been observed. In contrast to other bipolar or tetrapolar basidiomycetes that have either two linked mating-type-determining ( MAT ) loci or two MAT loci on separate chromosomes, in Malassezia species studied thus far the two MAT loci are arranged in a pseudobipolar configuration (linked on the same chromosome but capable of recombining). By incorporating newly generated chromosome-level genome assemblies, and an improved Malassezia phylogeny, we infer that the pseudobipolar arrangement was the ancestral state of this group and revealed six independent transitions to tetrapolarity, seemingly driven by centromere fission or translocations in centromere- flanking regions. Additionally, in an approach to uncover a sexual cycle, Malassezia furfur strains were engineered to express different MAT alleles in the same cell. The resulting strains produce hyphae reminiscent of early steps in sexual development and display upregulation of genes associated with sexual development as well as others encoding lipases and a protease potentially relevant for pathogenesis of the fungus. Our study reveals a previously unseen genomic relocation of mating-type loci in fungi and provides insight towards the discovery of a sexual cycle in Malassezia , with possible implications for pathogenicity.

Significance Statement

Malassezia , the dominant fungal group of the mammalian skin microbiome, is associated with numerous skin disorders. Sexual development and yeast-to-hyphae transitions, governed by genes at two mating-type ( MAT ) loci, are thought to be important for fungal pathogenicity. However, Malassezia sexual reproduction has never been observed. Here, we used chromosome-level assemblies and comparative genomics to uncover unforeseen transitions in MAT loci organization within Malassezia , possibly related with fragility of centromeric-associated regions. Additionally, by expressing different MAT alleles in the same cell, we show that Malassezia can undergo hyphal development and this phenotype is associated with increased expression of key mating genes along with other genes known to be virulence factors, providing a possible connection between hyphal development, sexual reproduction, and pathogenicity.

The Unique Homothallic Mating-Type Loci of the Fungal Tree Pathogens Chrysoporthe syzygiicola and Chrysoporthe zambiensis from Africa

Chrysoporthe syzygiicola and C. zambiensis are ascomycete tree pathogens first described from Zambia, causing stem canker on Syzygium guineense and Eucalyptus grandis, respectively. The taxonomic descriptions of these two species were based on their anamorphic states, as no sexual states are known. The main purpose of this work was to use whole genome sequences to identify and define the mating-type (MAT1) loci of these two species. The unique MAT1 loci for C. zambiensis and C. syzygiicola consist of the MAT1-1-1, MAT1-1-2, and MAT1-2-1 genes, but the MAT1-1-3 gene is absent. Genes canonically associated with opposite mating types were present at the single mating-type locus, suggesting that C. zambiensis and C. syzygiicola have homothallic mating systems.

Stepwise recombination suppression around the mating-type locus in an ascomycete fungus with self-fertile spores

Recombination is often suppressed at sex-determining loci in plants and animals, and at self-incompatibility or mating-type loci in plants and fungi. In fungal ascomycetes, recombination suppression around the mating-type locus is associated with pseudo-homothallism, i.e. the production of self-fertile dikaryotic sexual spores carrying the two opposite mating types. This has been well studied in two species complexes from different families of Sordariales: Podospora anserina and Neurospora tetrasperma. However, it is unclear whether this intriguing association holds in other species. We show here that Schizothecium tetrasporum, a fungus from a third family in the order Sordariales, also produces mostly self-fertile dikaryotic spores carrying the two opposite mating types. This was due to a high frequency of second meiotic division segregation at the mating-type locus, indicating the occurrence of a single and systematic crossing-over event between the mating-type locus and the centromere, as in P. anserina. The mating-type locus has the typical Sordariales organization, plus a MAT1-1-1 pseudogene in the MAT1-2 haplotype. High-quality genome assemblies of opposite mating types and segregation analyses revealed a suppression of recombination in a region of 1.47 Mb around the mating-type locus. We detected three evolutionary strata, indicating a stepwise extension of recombination suppression. The three strata displayed no rearrangement or transposable element accumulation but gene losses and gene disruptions were present, and precisely at the strata margins. Our findings indicate a convergent evolution of self-fertile dikaryotic sexual spores across multiple ascomycete fungi. The particular pattern of meiotic segregation at the mating-type locus was associated with recombination suppression around this locus, that had extended stepwise. This association between pseudo-homothallism and recombination suppression across lineages and the presence of gene disruption at the strata limits are consistent with a recently proposed mechanism of sheltering deleterious alleles to explain stepwise recombination suppression.

Lessons from the meiotic recombination landscape of the ZMM deficient budding yeast Lachancea waltii

Meiotic recombination is a driving force for genome evolution, deeply characterized in a few model species, notably in the budding yeast Saccharomyces cerevisiae. Interestingly, Zip2, Zip3, Zip4, Spo16, Msh4, and Msh5, members of the so-called ZMM pathway that implements the interfering meiotic crossover pathway in S. cerevisiae, have been lost in Lachancea yeast species after the divergence of Lachancea kluyveri from the rest of the clade. In this context, after investigating meiosis in L. kluyveri, we determined the meiotic recombination landscape of Lachancea waltii. Attempts to generate diploid strains with fully hybrid genomes invariably resulted in strains with frequent whole-chromosome aneuploidy and multiple extended regions of loss of heterozygosity (LOH), which mechanistic origin is so far unclear. Despite the lack of multiple ZMM pro-crossover factors in L. waltii, numbers of crossovers and noncrossovers per meiosis were higher than in L. kluyveri but lower than in S. cerevisiae, for comparable genome sizes. Similar to L. kluyveri but opposite to S. cerevisiae, L. waltii exhibits an elevated frequency of zero-crossover bivalents. Lengths of gene conversion tracts for both crossovers and non-crossovers in L. waltii were comparable to those observed in S. cerevisiae and shorter than in L. kluyveri despite the lack of Mlh2, a factor limiting conversion tract size in S. cerevisiae. L. waltii recombination hotspots were not shared with either S. cerevisiae or L. kluyveri, showing that meiotic recombination hotspots can evolve at a rather limited evolutionary scale within budding yeasts. Finally, L. waltii crossover interference was reduced relative to S. cerevisiae, with interference being detected only in the 25 kb distance range. Detection of positive inference only at short distance scales in the absence of multiple ZMM factors required for interference-sensitive crossovers in other systems likely reflects interference between early recombination precursors such as DSBs.

Onset and stepwise extensions of recombination suppression are common in mating-type chromosomes of Microbotryum anther-smut fungi

Sex chromosomes and mating-type chromosomes can display large genomic regions without recombination. Recombination suppression often extended stepwise with time away from the sex- or mating-type-determining genes, generating evolutionary strata of differentiation between alternative sex or mating-type chromosomes. In anther-smut fungi of the Microbotryum genus, recombination suppression evolved repeatedly, linking the two mating-type loci and extended multiple times in regions distal to the mating-type genes. Here, we obtained high-quality genome assemblies of alternative mating types for four Microbotryum fungi. We found an additional event of independent chromosomal rearrangements bringing the two mating-type loci on the same chromosome followed by recombination suppression linking them. We also found, in a new clade analysed here, that recombination suppression between the two mating-type loci occurred in several steps, with first an ancestral recombination suppression between one of the mating-type locus and its centromere; later, completion of recombination suppression up to the second mating-type locus occurred independently in three species. The estimated dates of recombination suppression between the mating-type loci ranged from 0.15 to 3.58 million years ago. In total, this makes at least nine independent events of linkage between the mating-type loci across the Microbotryum genus. Several mating-type locus linkage events occurred through the same types of chromosomal rearrangements, where similar chromosome fissions at centromeres represent convergence in the genomic changes leading to the phenotypic convergence. These findings further highlight Microbotryum fungi as excellent models to study the evolution of recombination suppression.

Local adaptation in fungi

In this review, I explore the pervasive but underappreciated role of local adaptation in fungi. It has been difficult historically to study local adaptation in fungi because of the limited understanding of fungal species and their traits, but new hope has been offered with technological advances in sequencing. The filamentous nature of fungi invalidates some assumptions made in evolution because of their ability to exist as multinucleate entities with genetically different nuclei sharing the same cytoplasm. Many insights on local adaptation have come from studying fungi, and much of the empirical evidence gathered about local adaptation in the context of host-pathogen interactions comes from studying fungal virulence genes, drug resistance, and environmental adaptation. Together, these insights paint a picture of the variety of processes involved in fungal local adaptation and their connections to the unusual cell biology of Fungi (multinucleate, filamentous habit), but there is much that remains unknown, with major gaps in our knowledge of fungal species, their phenotypes, and the ways by which they adapt to local conditions.

Homokaryotic High-Quality Genome Assembly of Medicinal Fungi Wolfiporia hoelen Reveals Auto-Regulation and High-Temperature Adaption of Probable Two-Speed Genome

Sclerotia of Wolfiporia hoelen are one of the most important traditional Chinese medicines and are commonly used in China, Japan, Korea, and other Asian countries. In the present study, we presented the first high-quality homokaryotic genome of W. hoelen with 14 chromosomes which was evaluated with assembly index, telomere position detection, and whole-genome collinearity. A 64.44 Mb genome was assembled with a Contig N50 length of 3.76 Mb. The imbalanced distribution of transposons and chromosome characters revealed the probable two-speed genome of W. hoelen. High consistency between methylation and transposon conserved the genome stability. The expansion of the gene family about signal transduction and nutritional transport has intimate relationships with sclerotial formation. Up-regulation of expression for distinctive decomposition enzymes, ROS clearance genes, biosynthesis of unsaturated fatty acids, and change of the cell wall components maintained high-speed growth of mycelia that may be the high-temperature adaption strategy of W. hoelen. Further, the analysis of mating-control genes demonstrated that HD3 probably had no function on mating recognition, with the HD protein in a distant genetic with known species. Overall, the high-quality genome of W. hoelen provided crucial information for genome structure and stability, high-temperature adaption, and sexual and asexual process.

Association mapping of colour variation in a butterfly provides evidence that a supergene locks together a cluster of adaptive loci

Supergenes are genetic architectures associated with discrete and concerted variation in multiple traits. It has long been suggested that supergenes control these complex polymorphisms by suppressing recombination between sets of coadapted genes. However, because recombination suppression hinders the dissociation of the individual effects of genes within supergenes, there is still little evidence that supergenes evolve by tightening linkage between coadapted genes. Here, combining a landmark-free phenotyping algorithm with multivariate genome-wide association studies, we dissected the genetic basis of wing pattern variation in the butterfly Heliconius numata. We show that the supergene controlling the striking wing pattern polymorphism displayed by this species contains several independent loci associated with different features of wing patterns. The three chromosomal inversions of this supergene suppress recombination between these loci, supporting the hypothesis that they may have evolved because they captured beneficial combinations of alleles. Some of these loci are, however, associated with colour variations only in a subset of morphs where the phenotype is controlled by derived inversion forms, indicating that they were recruited after the formation of the inversions. Our study shows that supergenes and clusters of adaptive loci in general may form via the evolution of chromosomal rearrangements suppressing recombination between co-adapted loci but also via the subsequent recruitment of linked adaptive mutations. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

Organization and Unconventional Integration of the Mating-Type Loci in Morchella Species

True morels (Morchella spp.) are a group of delicious fungi in high demand worldwide, and some species of morels have been successfully cultivated in recent years. To better understand the sexual reproductive mechanisms of these fungi, we characterized the structure of the mating-type loci from ten morel species, and seven of them were obtained using long-range PCR amplification. Among the studied species, eight were heterothallic, two were homothallic, and four types of composition were observed in the MAT loci. In three of the five black morel species, the MAT1-1-1, MAT1-1-10, and MAT1-1-11 genes were in the MAT1-1 idiomorph, and only the MAT1-2-1 gene was in the MAT1-2 idiomorph, while an integration event occurred in the other two species and resulted in the importation of the MAT1-1-11 gene into the MAT1-2 idiomorph and survival as a truncated fragment in the MAT1-1 idiomorph. However, the MAT1-1-11 gene was not available in the four yellow morels and one blushing morel species. M. rufobrunnea, a representative species of the earliest diverging branch of true morels, along with another yellow morel Mes-15, were confirmed to be homothallic, and the MAT1-1-1, MAT1-1-10, and MAT1-2-1 genes were arranged in a tandem array. Therefore, we hypothesized that homothallism should be the ancestral reproductive state in Morchella. RT-PCR analyses revealed that four mating genes could be constitutively expressed, while the MAT1-1-10 gene underwent alternative splicing to produce different splice variants.

Sheltering of deleterious mutations explains the stepwise extension of recombination suppression on sex chromosomes and other supergenes

Many organisms have sex chromosomes with large nonrecombining regions that have expanded stepwise, generating "evolutionary strata" of differentiation. The reasons for this remain poorly understood, but the principal hypotheses proposed to date are based on antagonistic selection due to differences between sexes. However, it has proved difficult to obtain empirical evidence of a role for sexually antagonistic selection in extending recombination suppression, and antagonistic selection has been shown to be unlikely to account for the evolutionary strata observed on fungal mating-type chromosomes. We show here, by mathematical modeling and stochastic simulation, that recombination suppression on sex chromosomes and around supergenes can expand under a wide range of parameter values simply because it shelters recessive deleterious mutations, which are ubiquitous in genomes. Permanently heterozygous alleles, such as the male-determining allele in XY systems, protect linked chromosomal inversions against the expression of their recessive mutation load, leading to the successive accumulation of inversions around these alleles without antagonistic selection. Similar results were obtained with models assuming recombination-suppressing mechanisms other than chromosomal inversions and for supergenes other than sex chromosomes, including those without XY-like asymmetry, such as fungal mating-type chromosomes. However, inversions capturing a permanently heterozygous allele were found to be less likely to spread when the mutation load segregating in populations was lower (e.g., under large effective population sizes or low mutation rates). This may explain why sex chromosomes remain homomorphic in some organisms but are highly divergent in others. Here, we model a simple and testable hypothesis explaining the stepwise extensions of recombination suppression on sex chromosomes, mating-type chromosomes, and supergenes in general.

Population Genomics Reveals Population Structure and Mating-Type Loci in Marssonina brunnea

Marssonina brunnea is an important fungal pathogen of poplar trees. We collected 32 M. brunnea f.sp. multigermtubi (MbMu) and three M. brunnea f.sp. monogermtubi (MbMo) isolates from four poplar species in three Chinese regions and performed genome resequencing. An annotation of SNPs of MbMu indicated that the SNPs potentially have a functional effect on 69.2% of the predicted genes. Using the SNP dataset of nonredundant isolates, a structure and principal component analysis revealed that MbMu and MbMo belong to two genetically distinct populations. By contrast, subpopulation structures could not be found among MbMu isolates. A neighbor-net analysis and a homoplasy index test provided evidence of recombination among MbMu isolates. The short distance (109–174 bp) of linkage disequilibrium half-decay supported the presence of a high level of recombination in the MbMu population. The genetic architectures of the MAT loci of MbMu and MbMo were revealed by searching genome assemblies or by homology-based cloning, and a BLAST search verified each isolate carrying one of the two opposite MAT loci. This study revealed that the MbMu population contains a wide range of functional variants, shows high-frequency recombination, and exhibits heterothallic mating systems, indicating high evolutionary potential and a resultant threat to poplar plantations.

Tempo of degeneration across independently evolved non-recombining regions

Recombination is beneficial over the long term, allowing more effective selection. Despite long-term advantages of recombination, local recombination suppression can evolve and lead to genomic degeneration, in particular on sex chromosomes. Here, we investigated the tempo of degeneration in nonrecombining regions, that is, the function curve for the accumulation of deleterious mutations over time, leveraging on 22 independent events of recombination suppression identified on mating-type chromosomes of anther-smut fungi, including newly identified ones. Using previously available and newly generated high-quality genome assemblies of alternative mating types of 13 Microbotryum species, we estimated degeneration levels in terms of accumulation of nonoptimal codons and nonsynonymous substitutions in nonrecombining regions. We found a reduced frequency of optimal codons in the nonrecombining regions compared with autosomes, that was not due to less frequent GC-biased gene conversion or lower ancestral expression levels compared with recombining regions. The frequency of optimal codons rapidly decreased following recombination suppression and reached an asymptote after ca. 3 Ma. The strength of purifying selection remained virtually constant at d_N/d_S = 0.55, that is, at an intermediate level between purifying selection and neutral evolution. Accordingly, nonsynonymous differences between mating-type chromosomes increased linearly with stratum age, at a rate of 0.015 per My. We thus develop a method for disentangling effects of reduced selection efficacy from GC-biased gene conversion in the evolution of codon usage and we quantify the tempo of degeneration in nonrecombining regions, which is important for our knowledge on genomic evolution and on the maintenance of regions without recombination.

Insights on life cycle and cell identity regulatory circuits for unlocking genetic improvement in Zygosaccharomyces and Kluyveromyces yeasts

Evolution has provided a vast diversity of yeasts that play fundamental roles in nature and society. This diversity is not limited to genotypically homogeneous species with natural interspecies hybrids and allodiploids that blur species boundaries frequently isolated. Thus, life cycle and the nature of breeding systems have profound effects on genome variation, shaping heterozygosity, genotype diversity and ploidy level. The apparent enrichment of hybrids in industry-related environments suggests that hybridization provides an adaptive route against stressors and creates interest in developing new hybrids for biotechnological uses. For example, in the Saccharomyces genus where regulatory circuits controlling cell identity, mating competence and meiosis commitment have been extensively studied, this body of knowledge is being used to combine interesting traits into synthetic F1 hybrids, to bypass F1 hybrid sterility and to dissect complex phenotypes by bulk segregant analysis. Although these aspects are less known in other industrially promising yeasts, advances in whole-genome sequencing and analysis are changing this and new insights are being gained, especially in the food-associated genera Zygosaccharomyces and Kluyveromyces. We discuss this new knowledge and highlight how deciphering cell identity circuits in these lineages will contribute significantly to identify the genetic determinants underpinning complex phenotypes and open new avenues for breeding programmes.

Evolutionary Genomics of Sex-Related Chromosomes at the Base of the Green Lineage

Although sex is now accepted as a ubiquitous and ancestral feature of eukaryotes, direct observation of sex is still lacking in most unicellular eukaryotic lineages. Evidence of sex is frequently indirect and inferred from the identification of genes involved in meiosis from whole genome data and/or the detection of recombination signatures from genetic diversity in natural populations. In haploid unicellular eukaryotes, sex-related chromosomes are named mating-type (MTs) chromosomes and generally carry large genomic regions where recombination is suppressed. These regions have been characterized in Fungi and Chlorophyta and determine gamete compatibility and fusion. Two candidate MT+ and MT− alleles, spanning 450–650 kb, have recently been described in Ostreococcus tauri, a marine phytoplanktonic alga from the Mamiellophyceae class, an early diverging branch in the green lineage. Here, we investigate the architecture and evolution of these candidate MT+ and MT− alleles. We analyzed the phylogenetic profile and GC content of MT gene families in eight different genomes whose divergence has been previously estimated at up to 640 Myr, and found evidence that the divergence of the two MT alleles predates speciation in the Ostreococcus genus. Phylogenetic profiles of MT trans-specific polymorphisms in gametologs disclosed candidate MTs in two additional species, and possibly a third. These Mamiellales MT candidates are likely to be the oldest mating-type loci described to date, which makes them fascinating models to investigate the evolutionary mechanisms of haploid sex determination in eukaryotes.

Size Variation of the Nonrecombining Region on the Mating-Type Chromosomes in the Fungal Podospora anserina Species Complex

Sex chromosomes often carry large nonrecombining regions that can extend progressively over time, generating evolutionary strata of sequence divergence. However, some sex chromosomes display an incomplete suppression of recombination. Large genomic regions without recombination and evolutionary strata have also been documented around fungal mating-type loci, but have been studied in only a few fungal systems. In the model fungus Podospora anserina (Ascomycota, Sordariomycetes), the reference S strain lacks recombination across a 0.8-Mb region around the mating-type locus. The lack of recombination in this region ensures that nuclei of opposite mating types are packaged into a single ascospore (pseudohomothallic lifecycle). We found evidence for a lack of recombination around the mating-type locus in the genomes of ten P. anserina strains and six closely related pseudohomothallic Podospora species. Importantly, the size of the nonrecombining region differed between strains and species, as indicated by the heterozygosity levels around the mating-type locus and experimental selfing. The nonrecombining region is probably labile and polymorphic, differing in size and precise location within and between species, resulting in occasional, but infrequent, recombination at a given base pair. This view is also supported by the low divergence between mating types, and the lack of strong linkage disequilibrium, chromosomal rearrangements, transspecific polymorphism and genomic degeneration. We found a pattern suggestive of evolutionary strata in P. pseudocomata. The observed heterozygosity levels indicate low but nonnull outcrossing rates in nature in these pseudohomothallic fungi. This study adds to our understanding of mating-type chromosome evolution and its relationship to mating systems.

Mating-Type Locus Organization and Mating-Type Chromosome Differentiation in the Bipolar Edible Button Mushroom Agaricus bisporus

In heterothallic basidiomycete fungi, sexual compatibility is restricted by mating types, typically controlled by two loci: PR, encoding pheromone precursors and pheromone receptors, and HD, encoding two types of homeodomain transcription factors. We analysed the single mating-type locus of the commercial button mushroom variety, Agaricus bisporus var. bisporus, and of the related variety burnettii. We identified the location of the mating-type locus using genetic map and genome information, corresponding to the HD locus, the PR locus having lost its mating-type role. We found the mip1 and β-fg genes flanking the HD genes as in several Agaricomycetes, two copies of the β-fg gene, an additional HD2 copy in the reference genome of A. bisporus var. bisporus and an additional HD1 copy in the reference genome of A. bisporus var. burnettii. We detected a 140 kb-long inversion between mating types in an A. bisporus var. burnettii heterokaryon, trapping the HD genes, the mip1 gene and fragments of additional genes. The two varieties had islands of transposable elements at the mating-type locus, spanning 35 kb in the A. bisporus var. burnettii reference genome. Linkage analyses showed a region with low recombination in the mating-type locus region in the A. bisporus var. burnettii variety. We found high differentiation between β-fg alleles in both varieties, indicating an ancient event of recombination suppression, followed more recently by a suppression of recombination at the mip1 gene through the inversion in A. bisporus var. burnettii and a suppression of recombination across whole chromosomes in A. bisporus var. bisporus, constituting stepwise recombination suppression as in many other mating-type chromosomes and sex chromosomes.