Evolutionary strata on young mating-type chromosomes despite the lack of sexual antagonism

Consequences of partially recessive deleterious genetic variation for the evolution of inversions suppressing recombination between sex chromosomes1

The evolution of suppressed recombination between sex chromosomes is widely hypothesized to be driven by sexually antagonistic selection (SA), where tighter linkage between the sex-determining gene(s) and nearby SA loci is favored when it couples male-beneficial alleles to the proto-Y chromosome, and female-beneficial alleles to the proto-X. Although difficult to test empirically, the SA selection hypothesis overshadows several alternatives, including an incomplete but often-repeated "sheltering" hypothesis which suggests that expansion of the sex-linked region (SLR) reduces the homozygous expression of deleterious mutations at selected loci. Here, we use population genetic models to evaluate the consequences of partially recessive deleterious mutational variation for the evolution of otherwise neutral chromosomal inversions expanding the SLR on proto-Y chromosomes. Both autosomal and SLR-expanding inversions face a race against time: lightly-loaded inversions are initially beneficial, but eventually become deleterious as they accumulate new mutations, after which their chances of fixing become negligible. In contrast, initially unloaded inversions eventually become neutral as their deleterious load reaches the same equilibrium as non-inverted haplotypes. Despite the differences in inheritance and indirect selection, SLR-expanding inversions exhibit similar evolutionary dynamics to autosomal inversions over many biologically plausible parameter conditions. Differences emerge when the population average mutation load is quite high; in this case large autosomal inversions that are lucky enough to be mutation-free can rise to intermediate to high frequencies where selection in homozygotes becomes important (Y-linked inversions never appear as homozygous karyotypes); conditions requiring either high mutation rates, highly recessive deleterious mutations, weak selection, or a combination thereof.

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.

Repeated co-option of HMG-box genes for sex determination in brown algae and animals

In many eukaryotes, genetic sex determination is not governed by XX/XY or ZW/ZZ systems but by a specialized region on the poorly studied U (female) or V (male) sex chromosomes. Previous studies have hinted at the existence of a dominant male-sex factor on the V chromosome in brown algae, a group of multicellular eukaryotes distantly related to animals and plants. The nature of this factor has remained elusive. Here, we demonstrate that an HMG-box gene acts as the male-determining factor in brown algae, mirroring the role HMG-box genes play in sex determination in animals. Over a billion-year evolutionary timeline, these lineages have independently co-opted the HMG box for male determination, representing a paradigm for evolution's ability to recurrently use the same genetic "toolkit" to accomplish similar tasks.

Karyotype Diversification and Chromosome Rearrangements in Squamate Reptiles

Karyotype diversification represents an important, yet poorly understood, driver of evolution. Squamate reptiles are characterized by a high taxonomic diversity which is reflected at the karyotype level in terms of general structure, chromosome number and morphology, and insurgence of differentiated simple or multiple-sex-chromosome systems with either male or female heterogamety. The potential of squamate reptiles as unique model organisms in evolutionary cytogenetics has been recognised in recent years in several studies, which have provided novel insights into the chromosome evolutionary dynamics of different taxonomic groups. Here, we review and summarize the resulting complex, but promising, general picture from a systematic perspective, mapping some of the main squamate karyological characteristics onto their phylogenetic relationships. We highlight how all the major categories of balanced chromosome rearrangements contributed to the karyotype evolution in different taxonomic groups. We show that distinct karyotype evolutionary trends may occur, and coexist, with different frequencies in different clades. Finally, in light of the known squamate chromosome diversity and recent research advances, we discuss traditional and novel hypotheses on karyotype evolution and propose a scenario of circular karyotype 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.

Social antagonism facilitates supergene expansion in ants

Antagonistic selection has long been considered a major driver of the formation and expansion of sex chromosomes. For example, sexually antagonistic variation on an autosome can select for suppressed recombination between that autosome and the sex chromosome, leading to a neo-sex chromosome. Autosomal supergenes, chromosomal regions containing tightly linked variants affecting the same complex trait, share similarities with sex chromosomes, raising the possibility that sex chromosome evolution models can explain the evolution of genome structure and recombination in other contexts. We tested this premise in a Formica ant species, wherein we identified four supergene haplotypes on chromosome 3 underlying colony social organization and sex ratio. We discovered a novel rearranged supergene variant (9r) on chromosome 9 underlying queen miniaturization. The 9r is in strong linkage disequilibrium with one chromosome 3 haplotype (P2) found in multi-queen (polygyne) colonies. We suggest that queen miniaturization is strongly disfavored in the single-queen (monogyne) background and is thus socially antagonistic. As such, divergent selection experienced by ants living in alternative social "environments" (monogyne and polygyne) may have contributed to the emergence of a genetic polymorphism on chromosome 9 and associated queen-size dimorphism. Consequently, an ancestral polygyne-associated haplotype may have expanded to include the polymorphism on chromosome 9, resulting in a larger region of suppressed recombination spanning two chromosomes. This process is analogous to the formation of neo-sex chromosomes and consistent with models of expanding regions of suppressed recombination. We propose that miniaturized queens, 16%-20% smaller than queens without 9r, could be incipient intraspecific social parasites.

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.

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.

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.

On the evolution of variation in sexual reproduction through the prism of eukaryotic microbes

Almost all eukaryotes undergo sexual reproduction to generate diversity and select for fitness in their population pools. Interestingly, the systems by which sex is defined are highly diverse and can even differ between evolutionarily closely related species. While the most commonly known form of sex determination involves males and females in animals, eukaryotic microbes can have as many as thousands of different mating types for the same species. Furthermore, some species have found alternatives to sexual reproduction and prefer to grow clonally and yet undergo infrequent facultative sexual reproduction. These organisms are mainly invertebrates and microbes, but several examples are also present among vertebrates suggesting that alternative modes of sexual reproduction evolved multiple times throughout evolution. In this review, we summarize the sex-determination modes and variants of sexual reproduction found across the eukaryotic tree of life and suggest that eukaryotic microbes provide unique opportunities to study these processes in detail. We propose that understanding variations in modes of sexual reproduction can serve as a foundation to study the evolution of sex and why and how it evolved in the first place.

A divergent RWP-RK transcription factor determines mating type in heterothallic Closterium

The Closterium peracerosum-strigosum-littorale complex (Closterium, Zygnematophyceae) has an isogamous mating system. Members of the Zygnematophyceae are the closest relatives to extant land plants and are distantly related to chlorophytic models, for which a genetic basis of mating type (MT) determination has been reported. We thus investigated MT determination in Closterium. We sequenced genomes representing the two MTs, mt+ and mt-, in Closterium and identified CpMinus1, a gene linked to the mt- phenotype. We analyzed its function using reverse genetics methods. CpMinus1 encodes a divergent RWP-RK domain-containing-like transcription factor and is specifically expressed during gamete differentiation. Introduction of CpMinus1 into an mt+ strain was sufficient to convert it to a phenotypically mt- strain, while CpMinus1-knockout mt- strains were phenotypically mt+. We propose that CpMinus1 is the major MT determinant that acts by evoking the mt- phenotype and suppressing the mt+ phenotype in heterothallic Closterium. CpMinus1 likely evolved independently in the Zygnematophyceae lineage, which lost an egg-sperm anisogamous mating system. mt- specific regions possibly constitute an MT locus flanked by common sequences that undergo some recombination.

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.

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.

Sex in protists: A new perspective on the reproduction mechanisms of trypanosomatids

The Protist kingdom individuals are the most ancestral representatives of eukaryotes. They have inhabited Earth since ancient times and are currently found in the most diverse environments presenting a great heterogeneity of life forms. The unicellular and multicellular algae, photosynthetic and heterotrophic organisms, as well as free-living and pathogenic protozoa represents the protist group. The evolution of sex is directly associated with the origin of eukaryotes being protists the earliest protagonists of sexual reproduction on earth. In eukaryotes, the recombination through genetic exchange is a ubiquitous mechanism that can be stimulated by DNA damage. Scientific evidences support the hypothesis that reactive oxygen species (ROS) induced DNA damage can promote sexual recombination in eukaryotes which might have been a decisive factor for the origin of sex. The fact that some recombination enzymes also participate in meiotic sex in modern eukaryotes reinforces the idea that sexual reproduction emerged as consequence of specific mechanisms to cope with mutations and alterations in genetic material. In this review we will discuss about origin of sex and different strategies of evolve sexual reproduction in some protists such that cause human diseases like malaria, toxoplasmosis, sleeping sickness, Chagas disease, and leishmaniasis.

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'.

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.

Consequences of partially recessive deleterious genetic variation for the evolution of inversions suppressing recombination between sex chromosomes

The evolution of suppressed recombination between sex chromosomes is widely hypothesized to be driven by sexually antagonistic selection (SA), where tighter linkage between the sex-determining gene(s) and nearby SA loci is favored when it couples male-beneficial alleles to the proto-Y chromosome, and female-beneficial alleles to the proto-X. Despite limited empirical evidence, the SA selection hypothesis overshadows several alternatives, including an incomplete but often-repeated "sheltering hypothesis" that suggests that expansion of the sex-linked region (SLR) reduces homozygous expression of partially recessive deleterious mutations at selected loci. Here, we use population genetic models to evaluate the consequences of deleterious mutational variation for the evolution of neutral chromosomal inversions expanding the SLR on proto-Y chromosomes. We find that SLR-expanding inversions face a race against time: lightly loaded inversions are initially beneficial, but eventually become deleterious as they accumulate new mutations, and must fix before this window of opportunity closes. The outcome of this race is strongly influenced by inversion size, the mutation rate, and the dominance coefficient of deleterious mutations. Yet, small inversions have elevated fixation probabilities relative to neutral expectations for biologically plausible parameter values. Our results demonstrate that deleterious genetic variation can plausibly drive recombination suppression in small steps and would be most consistent with empirical patterns of small evolutionary strata or gradual recombination arrest.

Modelling the consequences of the dikaryotic life cycle of mushroom-forming fungi on genomic conflict

Generally, sexual organisms contain two haploid genomes, one from each parent, united in a single diploid nucleus of the zygote which links their fate during growth. A fascinating exception to this is Basidiomycete fungi, where the two haploid genomes remain separate in a dikaryon, retaining the option to fertilize subsequent monokaryons encountered. How the ensuing nuclear competition influences the balance of selection within and between individuals is largely unexplored. We test the consequences of the dikaryotic life cycle for mating success and mycelium-level fitness components. We assume a trade-off between mating fitness at the level of the haploid nucleus and fitness of the fungal mycelium. We show that the maintenance of fertilization potential by dikaryons leads to a higher proportion of fertilized monokaryons, but that the ensuing intradikaryon selection for increased nuclear mating fitness leads to reduced mycelium fitness relative to a diploid life cycle. However, this fitness reduction is lower compared to a hypothetical life cycle where dikaryons can also exchange nuclei. Prohibition of fusion between dikaryons therefore reduces the level of nuclear parasitism. The number of loci influencing fitness is an important determinant of the degree to which average mycelium-level fitness is reduced. The results of this study crucially hinge upon a trade-off between nucleus and mycelium-level fitness. We discuss the evidence for this assumption and the implications of an alternative that there is a positive relationship between nucleus and mycelium-level fitness.

Convergence Analysis of Rust Fungi and Anther Smuts Reveals Their Common Molecular Adaptation to a Phytoparasitic Lifestyle

Convergent evolution between distantly related taxa often mirrors adaptation to similar environments. Rust fungi and anther smuts, which belong to different classes in Pucciniomycotina, have independently evolved a phytoparasitic lifestyle, representing an example of convergent evolution in the fungal kingdom. To investigate their adaptations and the genetic bases underlying their phytoparasitic lifestyles, we performed genome-wide convergence analysis of amino acid substitutions, evolutionary rates, and gene gains and losses. Convergent substitutions were detected in ATPeV0D and RP-S27Ae, two genes important for the generation of turgor pressure and ribosomal biosynthesis, respectively. A total of 51 positively selected genes were identified, including eight genes associated with translation and three genes related to the secretion pathway. In addition, rust fungi and anther smuts contained more proteins associated with oligopeptide transporters and vacuolar proteases than did other fungi. For rust fungi and anther smuts, these forms of convergence suggest four adaptive mechanisms for a phytoparasitic lifestyle: 1) reducing the metabolic demand for hyphal growth and penetration at the pre-penetration stage, 2) maintaining the efficiency of protein synthesis during colonization, 3) ensuring the normal secretion of rapidly evolving secreted proteins, and 4) improving the capacity for oligopeptide metabolism. Our results are the first to shed light on the genetic convergence mechanisms and molecular adaptation underlying phytoparasitic lifestyles in fungi.

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.

Y recombination arrest and degeneration in the absence of sexual dimorphism

Current theory proposes that degenerated sex chromosomes-such as the mammalian Y-evolve through three steps: (i) recombination arrest, linking male-beneficial alleles to the Y chromosome; (ii) Y degeneration, resulting from the inefficacy of selection in the absence of recombination; and (iii) dosage compensation, correcting the resulting low expression of X-linked genes in males. We investigate a model of sex chromosome evolution that incorporates the coevolution of cis and trans regulators of gene expression. We show that the early emergence of dosage compensation favors the maintenance of Y-linked inversions by creating sex-antagonistic regulatory effects. This is followed by degeneration of these nonrecombining inversions caused by regulatory divergence between the X and Y chromosomes. In contrast to current theory, the whole process occurs without any selective pressure related to sexual dimorphism.

Obligate sexual reproduction of a homothallic fungus closely related to the Cryptococcus pathogenic species complex

<title>eLife digest</title>.

Fungi are enigmatic organisms that flourish in soil, on decaying plants, or during infection of animals or plants. Growing in myriad forms, from single-celled yeast to multicellular molds and mushrooms, fungi have also evolved a variety of strategies to reproduce. Normally, fungi reproduce in one of two ways: either they reproduce asexually, with one individual producing a new individual identical to itself, or they reproduce sexually, with two individuals of different 'mating types' contributing to produce a new individual. However, individuals of some species exhibit 'homothallism' or self-fertility: these individuals can produce reproductive cells that are universally compatible, and therefore can reproduce sexually with themselves or with any other cell in the population.

Homothallism has evolved multiple times throughout the fungal kingdom, suggesting it confers advantage when population numbers are low or mates are hard to find. Yet some homothallic fungi been overlooked compared to heterothallic species, whose mating types have been well characterised. Understanding the genetic basis of homothallism and how it evolved in different species can provide insights into pathogenic species that cause fungal disease.

With that in mind, Passer, Clancey et al. explored the genetic basis of homothallism in Cryptococcus depauperatus, a close relative of C. neoformans, a species that causes fungal infections in humans. A combination of genetic sequencing techniques and experiments were applied to analyse, compare, and manipulate C. depauperatus' genome to see how this species evolved self-fertility.

Passer, Clancey et al. showed that C. depauperatus evolved the ability to reproduce sexually by itself via a unique evolutionary pathway. The result is a form of homothallism never reported in fungi before. C. depauperatus lost some of the genes that control mating in other species of fungi, and acquired genes from the opposing mating types of a heterothallic ancestor to become self-fertile.

Passer, Clancey et al. also found that, unlike other Cryptococcus species that switch between asexual and sexual reproduction, C. depauperatus grows only as long, branching filaments called hyphae, a sexual form. The species reproduces sexually with itself throughout its life cycle and is unable to produce a yeast (asexual) form, in contrast to other closely related species.

This work offers new insights into how different modes of sexual reproduction have evolved in fungi. It also provides another interesting case of how genome plasticity and evolutionary pressures can produce similar outcomes, homothallism, via different evolutionary paths. Lastly, assembling the complete genome of C. depauperatus will foster comparative studies between pathogenic and non-pathogenic Cryptococcus species.

Draft genome sequences of strains CBS6241 and CBS6242 of the basidiomycetous yeast Filobasidium floriforme

The Tremellomycetes are a species-rich group within the basidiomycete fungi; however, most analyses of this group to date have focused on pathogenic Cryptococcus species within the order Tremellales. Recent genome-assisted studies of other Tremellomycetes have identified interesting features with respect to biotechnological applications as well as the evolution of genes involved in mating and sexual development. Here, we report genome sequences of two strains of Filobasidium floriforme, a species from the order Filobasidiales, which branches basally to the Tremellales, Trichosporonales, and Holtermanniales. The assembled genomes of strains CBS6241 and CBS6242 are 27.4 Mb and 26.4 Mb in size, respectively, with 8314 and 7695 predicted protein-coding genes. Overall sequence identity at nucleic acid level between the strains is 97%. Among the predicted genes are pheromone precursor and pheromone receptor genes as well as two genes encoding homedomain (HD) transcription factors, which are predicted to be part of the mating type (MAT) locus. Sequence analysis indicates that CBS6241 and CBS6242 carry different alleles for both the pheromone/receptor genes as well as the HD transcription factors. Orthology inference identified 1482 orthogroups exclusively found in F. floriforme, some of which were involved in carbohydrate transport and metabolism. Subsequent CAZyme repertoire characterization identified 267 and 247 enzymes for CBS6241 and CBS6242, respectively, the second highest number of CAZymes among the analyzed Tremellomycete species. In addition, F. floriforme contains five CAZymes absent in other species and several plant-cell-wall degrading CAZymes with the highest copy number in Tremellomycota, indicating the biotechnological potential of this species.

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.

A neutral model for the loss of recombination on sex chromosomes

The loss of recombination between sex chromosomes has occurred repeatedly throughout nature, with important implications for their subsequent evolution. Explanations for this remarkable convergence have generally invoked only adaptive processes (e.g. sexually antagonistic selection); however, there is still little evidence for these hypotheses. Here we propose a model in which recombination on sex chromosomes is lost due to the neutral accumulation of sequence divergence adjacent to (and thus, in linkage disequilibrium with) the sex determiner. Importantly, we include in our model the fact that sequence divergence, in any form, reduces the probability of recombination between any two sequences. Using simulations, we show that, under certain conditions, a region of suppressed recombination arises and expands outwards from the sex-determining locus, under purely neutral processes. Further, we show that the rate and pattern of recombination loss are sensitive to the pre-existing recombination landscape of the genome and to sex differences in recombination rates, with patterns consistent with evolutionary strata emerging under some conditions. We discuss the applicability of these results to natural systems. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)'.

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.

The evolving species concepts used for yeasts: from phenotypes and genomes to speciation networks

Here we review how evolving species concepts have been applied to understand yeast diversity. Initially, a phenotypic species concept was utilized taking into consideration morphological aspects of colonies and cells, and growth profiles. Later the biological species concept was added, which applied data from mating experiments. Biophysical measurements of DNA similarity between isolates were an early measure that became more broadly applied with the advent of sequencing technology, leading to a sequence-based species concept using comparisons of parts of the ribosomal DNA. At present phylogenetic species concepts that employ sequence data of rDNA and other genes are universally applied in fungal taxonomy, including yeasts, because various studies revealed a relatively good correlation between the biological species concept and sequence divergence. The application of genome information is becoming increasingly common, and we strongly recommend the use of complete, rather than draft genomes to improve our understanding of species and their genome and genetic dynamics. Complete genomes allow in-depth comparisons on the evolvability of genomes and, consequently, of the species to which they belong. Hybridization seems a relatively common phenomenon and has been observed in all major fungal lineages that contain yeasts. Note that hybrids may greatly differ in their post-hybridization development. Future in-depth studies, initially using some model species or complexes may shift the traditional species concept as isolated clusters of genetically compatible isolates to a cohesive speciation network in which such clusters are interconnected by genetic processes, such as hybridization.

Whole-genome sequencing and genome regions of special interest: Lessons from major histocompatibility complex, sex determination, and plant self-incompatibility

Whole‐genome sequencing of non‐model organisms is now widely accessible and has allowed a range of questions in the field of molecular ecology to be investigated with greater power. However, some genomic regions that are of high biological interest remain problematic for assembly and data‐handling. Three such regions are the major histocompatibility complex (MHC), sex‐determining regions (SDRs) and the plant self‐incompatibility locus (S‐locus). Using these as examples, we illustrate the challenges of both assembling and resequencing these highly polymorphic regions and how bioinformatic and technological developments are enabling new approaches to their study. Mapping short‐read sequences against multiple alternative references improves genotyping comprehensiveness at the S‐locus thereby contributing to more accurate assessments of allelic frequencies. Long‐read sequencing, producing reads of several tens to hundreds of kilobase pairs in length, facilitates the assembly of such regions as single sequences can span the multiple duplicated gene copies of the MHC region, and sequence through repetitive stretches and translocations in SDRs and S‐locus haplotypes. These advances are adding value to short‐read genome resequencing approaches by allowing, for example, more accurate haplotype phasing across longer regions. Finally, we assessed further technical improvements, such as nanopore adaptive sequencing and bioinformatic tools using pangenomes, which have the potential to further expand our knowledge of a number of genomic regions that remain challenging to study with classical resequencing approaches.

Pervasive Phenotypic Impact of a Large Nonrecombining Introgressed Region in Yeast

To explore the origin of the diversity observed in natural populations, many studies have investigated the relationship between genotype and phenotype. In yeast species, especially in Saccharomyces cerevisiae, these studies are mainly conducted using recombinant offspring derived from two genetically diverse isolates, allowing to define the phenotypic effect of genetic variants. However, large genomic variants such as interspecies introgressions are usually overlooked even if they are known to modify the genotype–phenotype relationship. To have a better insight into the overall phenotypic impact of introgressions, we took advantage of the presence of a 1-Mb introgressed region, which lacks recombination and contains the mating-type determinant in the Lachancea kluyveri budding yeast. By performing linkage mapping analyses in this species, we identified a total of 89 loci affecting growth fitness in a large number of conditions and 2,187 loci affecting gene expression mostly grouped into two major hotspots, one being the introgressed region carrying the mating-type locus. Because of the absence of recombination, our results highlight the presence of a sexual dimorphism in a budding yeast for the first time. Overall, by describing the phenotype–genotype relationship in the Lachancea kluyveri species, we expanded our knowledge on how genetic characteristics of large introgression events can affect the phenotypic landscape.

The Genomic Architecture and Evolutionary Fates of Supergenes

Supergenes are genomic regions containing sets of tightly linked loci that control multi-trait phenotypic polymorphisms under balancing selection. Recent advances in genomics have uncovered significant variation in both the genomic architecture as well as the mode of origin of supergenes across diverse organismal systems. Although the role of genomic architecture for the origin of supergenes has been much discussed, differences in the genomic architecture also subsequently affect the evolutionary trajectory of supergenes and the rate of degeneration of supergene haplotypes. In this review, we synthesize recent genomic work and historical models of supergene evolution, highlighting how the genomic architecture of supergenes affects their evolutionary fate. We discuss how recent findings on classic supergenes involved in governing ant colony social form, mimicry in butterflies, and heterostyly in flowering plants relate to theoretical expectations. Furthermore, we use forward simulations to demonstrate that differences in genomic architecture affect the degeneration of supergenes. Finally, we discuss implications of the evolution of supergene haplotypes for the long-term fate of balanced polymorphisms governed by supergenes.

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

Genomic regions determining sexual compatibility often display recombination suppression, as occurs in sex chromosomes, plant self-incompatibility loci and fungal mating-type loci. Regions lacking recombination can extend beyond the genes determining sexes or mating types, by several successive steps of recombination suppression. Here we review the evidence for recombination suppression around mating-type loci in fungi, sometimes encompassing vast regions of the mating-type chromosomes. The suppression of recombination at mating-type loci in fungi has long been recognized and maintains the multiallelic combinations required for correct compatibility determination. We review more recent evidence for expansions of recombination suppression beyond mating-type genes in fungi ('evolutionary strata'), which have been little studied and may be more pervasive than commonly thought. We discuss testable hypotheses for the ultimate (evolutionary) and proximate (mechanistic) causes for such expansions of recombination suppression, including (1) antagonistic selection, (2) association of additional functions to mating-type, such as uniparental mitochondria inheritance, (3) accumulation in the margin of nonrecombining regions of various factors, including deleterious mutations or transposable elements resulting from relaxed selection, or neutral rearrangements resulting from genetic drift. The study of recombination suppression in fungi could thus contribute to our understanding of recombination suppression expansion across a broader range of organisms.

When Sex Chromosomes Recombine Only in the Heterogametic Sex: Heterochiasmy and Heterogamety in Hyla Tree Frogs

Sex chromosomes are classically predicted to stop recombining in the heterogametic sex, thereby enforcing linkage between sex-determining (SD) and sex-antagonistic (SA) genes. With the same rationale, a pre-existing sex asymmetry in recombination is expected to affect the evolution of heterogamety, for example, a low rate of male recombination might favor transitions to XY systems, by generating immediate linkage between SD and SA genes. Furthermore, the accumulation of deleterious mutations on nonrecombining Y chromosomes should favor XY-to-XY transitions (which discard the decayed Y), but disfavor XY-to-ZW transitions (which fix the decayed Y as an autosome). Like many anuran amphibians, Hyla tree frogs have been shown to display drastic heterochiasmy (males only recombine at chromosome tips) and are typically XY, which seems to fit the above expectations. Instead, here we demonstrate that two species, H. sarda and H. savignyi, share a common ZW system since at least 11 Ma. Surprisingly, the typical pattern of restricted male recombination has been maintained since then, despite female heterogamety. Hence, sex chromosomes recombine freely in ZW females, not in ZZ males. This suggests that heterochiasmy does not constrain heterogamety (and vice versa), and that the role of SA genes in the evolution of sex chromosomes might have been overemphasized.

Evolutionary stasis of the pseudoautosomal boundary in strepsirrhine primates

Sex chromosomes are typically comprised of a non-recombining region and a recombining pseudoautosomal region. Accurately quantifying the relative size of these regions is critical for sex-chromosome biology both from a functional and evolutionary perspective. The evolution of the pseudoautosomal boundary (PAB) is well documented in haplorrhines (apes and monkeys) but not in strepsirrhines (lemurs and lorises). Here, we studied the PAB of seven species representing the main strepsirrhine lineages by sequencing a male and a female genome in each species and using sex differences in coverage to identify the PAB. We found that during primate evolution, the PAB has remained unchanged in strepsirrhines whereas several recombination suppression events moved the PAB and shortened the pseudoautosomal region in haplorrhines. Strepsirrhines are well known to have much lower sexual dimorphism than haplorrhines. We suggest that mutations with antagonistic effects between males and females have driven recombination suppression and PAB evolution in haplorrhines.

Higher Gene Flow in Sex-Related Chromosomes than in Autosomes during Fungal Divergence

Nonrecombining sex chromosomes are widely found to be more differentiated than autosomes among closely related species, due to smaller effective population size and/or to a disproportionally large-X effect in reproductive isolation. Although fungal mating-type chromosomes can also display large nonrecombining regions, their levels of differentiation compared with autosomes have been little studied. Anther-smut fungi from the Microbotryum genus are castrating pathogens of Caryophyllaceae plants with largely nonrecombining mating-type chromosomes. Using whole genome sequences of 40 fungal strains, we quantified genetic differentiation among strains isolated from the geographically overlapping North American species and subspecies of Silene virginica and S. caroliniana. We inferred that gene flow likely occurred at the early stages of divergence and then completely stopped. We identified large autosomal genomic regions with chromosomal inversions, with higher genetic divergence than the rest of the genomes and highly enriched in selective sweeps, supporting a role of rearrangements in preventing gene flow in genomic regions involved in ecological divergence. Unexpectedly, the nonrecombining mating-type chromosomes showed lower divergence than autosomes due to higher gene flow, which may be promoted by adaptive introgressions of less degenerated mating-type chromosomes. The fact that both mating-type chromosomes are always heterozygous and nonrecombining may explain such patterns that oppose to those found for XY or ZW sex chromosomes. The specific features of mating-type chromosomes may also apply to the UV sex chromosomes determining sexes at the haploid stage in algae and bryophytes and may help test general hypotheses on the evolutionary specificities of sex-related chromosomes.

Mating-type locus rearrangements and shifts in thallism states in Citrus-associated Phyllosticta species

Currently, eight Phyllosticta species are known to be associated with several Citrus hosts, incorporating diverse lifestyles: while some of them are endophytic (P. capitalensis and P. citribraziliensis), others are pathogenic (P. citriasiana, P. citricarpa, P. citrichinaensis and P. paracitricarpa). Sexual reproduction plays a key role in the interaction between these Phyllosticta species and their Citrus hosts, especially for the spread and persistence of the pathogenic species in the environment. Given this, differences in sexual reproduction strategies could be related to the differences in lifestyles. To evaluate this hypothesis, we characterized the mating-type loci of six Citrus-associated Phyllosticta species from whole genome assemblies. Mating-type genes in the Citrus-associated Phyllosticta species are highly variable in their sequence content, but the genomic locations and organization of the mating-type loci are conserved. Phyllosticta citriasiana, P. citribraziliensis, P. citricarpa and P. paracitricarpa are heterothallic, while P. capitalensis and P. citrichinaensis are homothallic. In addition, the P. citrichinaensis MAT1-2 idiomorph occurs in a separate location from the mating-type locus. Ancestral state reconstruction suggests that homothallism is the ancestral thallism state in Phyllosticta, with a shift to heterothallism in Phyllosticta species that are pathogenic to Citrus. Moreover, the homothallic strategies of P. capitalensis and P. citrichinaensis result from independent evolutionary events, as P. capitalensis locus likely represents the ancestral state, and P. citrichinaensis homothallism has risen through a reversion in a heterothallic ancestor and underwent remodelling events. As the pathogenic species P. citriasiana, P. citricarpa and P. paracitricarpa are heterothallic and incapable of selfing, disease management practices focused in preventing the occurrence of sexual reproduction could assist in the control of Citrus Black Spot and Citrus Tan Spot diseases. This study emphasizes the importance of studying Citrus-Phyllosticta interactions under evolutionary and genomic perspectives, as these approaches can provide valuable information about the association between Phyllosticta species and their hosts, and also serve as guidance for the improvement of disease management practices.

Genomic architecture and evolutionary antagonism drive allelic expression bias in the social supergene of red fire ants

Supergene regions maintain alleles of multiple genes in tight linkage through suppressed recombination. Despite their importance in determining complex phenotypes, our empirical understanding of early supergene evolution is limited. Here we focus on the young ‘social’ supergene of fire ants, a powerful system for disentangling the effects of evolutionary antagonism and suppressed recombination. We hypothesize that gene degeneration and social antagonism shaped the evolution of the fire ant supergene, resulting in distinct patterns of gene expression. We test these ideas by identifying allelic differences between supergene variants, characterizing allelic expression across populations, castes and body parts, and contrasting allelic expression biases with differences in expression between social forms. We find strong signatures of gene degeneration and gene-specific dosage compensation. On this background, a small portion of the genes has the signature of adaptive responses to evolutionary antagonism between social forms.

Red fire ants (Solenopsis invicta) are native to South America, but the species has spread to North America, Australia and New Zealand where it can be an invasive pest. A reason for this species’ invasiveness types of colonies : one with a single egg-laying queen and another with several queens. However, it is not possible to simply add more queens to a colony with one queen. Instead, the number of queens in a colony is controlled genetically, by a chromosome known as the ‘social chromosome’.

Like many other animals, red fire ants are diploid: their cells have two copies of each chromosome, which can carry two different versions of each gene. The social chromosome is no different, and it comes in two variants, SB and Sb. Each ant can therefore have either two SB chromosomes, leading to a colony with a single queen; or one SB chromosome and one Sb chromosome, leading to a colony with multiple queens. Ants with two copies of the Sb variant die when they are young, so the Sb version is inherited in a similar way to how the Y chromosome is passed on in humans. However, the social chromosome in red fire ants appeared less than one million years ago, making it much younger than the human Y chromosome, which is 180 million years old. This makes the social chromosome a good candidate for examining the early evolution of special chromosome variants that are only inherited.

How differences between the SB and the Sb chromosomes are evolving is an open question, however. Perhaps each version of the social chromosome has been optimised through natural selection to one colony type. Another suggestion is that the Sb chromosome has degenerated over time because its genes cannot be ‘reshuffled’ as they would be on normal chromosomes.

Martinez-Ruiz et al. compared genetic variants on the SB and Sb chromosomes, along with their expression in different types of ant colonies. The analysis showed that the Sb variant is in fact breaking down because of the lack of gene shuffling. This loss is compensated by intact copies of the same genes found on the SB variant, which explains why ants with the Sb variant can only survive if they also carry the SB version. Only a handful of genes on the social chromosomes appear to have been optimised by natural selection. Therefore Martinez-Ruiz et al. concluded the differences between the two chromosomes that lead to different colony types are collateral effects of Sb’s inability to reshuffle its genes.

This work reveals how a special chromosome similar to the Y chromosome in humans evolved. It also shows how multiple complex evolutionary forces can shape a species’ genetic makeup and social forms.

Towards yeast taxogenomics: lessons from novel species descriptions based on complete genome sequences

In recent years, ‘multi-omic’ sciences have affected all aspects of fundamental and applied biological research. Yeast taxonomists, though somewhat timidly, have begun to incorporate complete genomic sequences into the description of novel taxa, taking advantage of these powerful data to calculate more reliable genetic distances, construct more robust phylogenies, correlate genotype with phenotype and even reveal cryptic sexual behaviors. However, the use of genomic data in formal yeast species descriptions is far from widespread. The present review examines published examples of genome-based species descriptions of yeasts, highlights relevant bioinformatic approaches, provides recommendations for new users and discusses some of the challenges facing the genome-based systematics of yeasts.

Genome assembly of the basket willow, Salix viminalis, reveals earliest stages of sex chromosome expansion

BACKGROUND

Sex chromosomes have evolved independently multiple times in eukaryotes and are therefore considered a prime example of convergent genome evolution. Sex chromosomes are known to emerge after recombination is halted between a homologous pair of chromosomes, and this leads to a range of non-adaptive modifications causing gradual degeneration and gene loss on the sex-limited chromosome. However, the proximal causes of recombination suppression and the pace at which degeneration subsequently occurs remain unclear.

RESULTS

Here, we use long- and short-read single-molecule sequencing approaches to assemble and annotate a draft genome of the basket willow, Salix viminalis, a species with a female heterogametic system at the earliest stages of sex chromosome emergence. Our single-molecule approach allowed us to phase the emerging Z and W haplotypes in a female, and we detected very low levels of Z/W single-nucleotide divergence in the non-recombining region. Linked-read sequencing of the same female and an additional male (ZZ) revealed the presence of two evolutionary strata supported by both divergence between the Z and W haplotypes and by haplotype phylogenetic trees. Gene order is still largely conserved between the Z and W homologs, although the W-linked region contains genes involved in cytokinin signaling regulation that are not syntenic with the Z homolog. Furthermore, we find no support across multiple lines of evidence for inversions, which have long been assumed to halt recombination between the sex chromosomes.

CONCLUSIONS

Our data suggest that selection against recombination is a more gradual process at the earliest stages of sex chromosome formation than would be expected from an inversion and may result instead from the accumulation of transposable elements. Our results present a cohesive understanding of the earliest genomic consequences of recombination suppression as well as valuable insights into the initial stages of sex chromosome formation and regulation of sex differentiation.

Genome assembly of the basket willow, Salix viminalis, reveals earliest stages of sex chromosome expansion

BACKGROUND

Sex chromosomes have evolved independently multiple times in eukaryotes and are therefore considered a prime example of convergent genome evolution. Sex chromosomes are known to emerge after recombination is halted between a homologous pair of chromosomes, and this leads to a range of non-adaptive modifications causing gradual degeneration and gene loss on the sex-limited chromosome. However, the proximal causes of recombination suppression and the pace at which degeneration subsequently occurs remain unclear.

RESULTS

Here, we use long- and short-read single-molecule sequencing approaches to assemble and annotate a draft genome of the basket willow, Salix viminalis, a species with a female heterogametic system at the earliest stages of sex chromosome emergence. Our single-molecule approach allowed us to phase the emerging Z and W haplotypes in a female, and we detected very low levels of Z/W single-nucleotide divergence in the non-recombining region. Linked-read sequencing of the same female and an additional male (ZZ) revealed the presence of two evolutionary strata supported by both divergence between the Z and W haplotypes and by haplotype phylogenetic trees. Gene order is still largely conserved between the Z and W homologs, although the W-linked region contains genes involved in cytokinin signaling regulation that are not syntenic with the Z homolog. Furthermore, we find no support across multiple lines of evidence for inversions, which have long been assumed to halt recombination between the sex chromosomes.

CONCLUSIONS

Our data suggest that selection against recombination is a more gradual process at the earliest stages of sex chromosome formation than would be expected from an inversion and may result instead from the accumulation of transposable elements. Our results present a cohesive understanding of the earliest genomic consequences of recombination suppression as well as valuable insights into the initial stages of sex chromosome formation and regulation of sex differentiation.

Sex Chromosome Evolution: So Many Exceptions to the Rules

Genomic analysis of many nonmodel species has uncovered an incredible diversity of sex chromosome systems, making it possible to empirically test the rich body of evolutionary theory that describes each stage of sex chromosome evolution. Classic theory predicts that sex chromosomes originate from a pair of homologous autosomes and recombination between them is suppressed via inversions to resolve sexual conflict. The resulting degradation of the Y chromosome gene content creates the need for dosage compensation in the heterogametic sex. Sex chromosome theory also implies a linear process, starting from sex chromosome origin and progressing to heteromorphism. Despite many convergent genomic patterns exhibited by independently evolved sex chromosome systems, and many case studies supporting these theoretical predictions, emerging data provide numerous interesting exceptions to these long-standing theories, and suggest that the remarkable diversity of sex chromosomes is matched by a similar diversity in their evolution. For example, it is clear that sex chromosome pairs are not always derived from homologous autosomes. In addition, both the cause and the mechanism of recombination suppression between sex chromosome pairs remain unclear, and it may be that the spread of recombination suppression is a more gradual process than previously thought. It is also clear that dosage compensation can be achieved in many ways, and displays a range of efficacy in different systems. Finally, the remarkable turnover of sex chromosomes in many systems, as well as variation in the rate of sex chromosome divergence, suggest that assumptions about the inevitable linearity of sex chromosome evolution are not always empirically supported, and the drivers of the birth-death cycle of sex chromosome evolution remain to be elucidated. Here, we concentrate on how the diversity in sex chromosomes across taxa highlights an equal diversity in each stage of sex chromosome evolution.

Differential Gene Expression between Fungal Mating Types Is Associated with Sequence Degeneration

Degenerative mutations in non-recombining regions, such as in sex chromosomes, may lead to differential expression between alleles if mutations occur stochastically in one or the other allele. Reduced allelic expression due to degeneration has indeed been suggested to occur in various sex-chromosome systems. However, whether an association occurs between specific signatures of degeneration and differential expression between alleles has not been extensively tested, and sexual antagonism can also cause differential expression on sex chromosomes. The anther-smut fungus Microbotryum lychnidis-dioicae is ideal for testing associations between specific degenerative signatures and differential expression because 1) there are multiple evolutionary strata on the mating-type chromosomes, reflecting successive recombination suppression linked to mating-type loci; 2) separate haploid cultures of opposite mating types help identify differential expression between alleles; and 3) there is no sexual antagonism as a confounding factor accounting for differential expression. We found that differentially expressed genes were enriched in the four oldest evolutionary strata compared with other genomic compartments, and that, within compartments, several signatures of sequence degeneration were greater for differentially expressed than non-differentially expressed genes. Two particular degenerative signatures were significantly associated with lower expression levels within differentially expressed allele pairs: upstream insertion of transposable elements and mutations truncating the protein length. Other degenerative mutations associated with differential expression included nonsynonymous substitutions and altered intron or GC content. The association between differential expression and allele degeneration is relevant for a broad range of taxa where mating compatibility or sex is determined by genes located in large regions where recombination is suppressed.

Fitness differences suppress the number of mating types in evolving isogamous species

Sexual reproduction is not always synonymous with the existence of two morphologically different sexes; isogamous species produce sex cells of equal size, typically falling into multiple distinct self-incompatible classes, termed mating types. A long-standing open question in evolutionary biology is: what governs the number of these mating types across species? Simple theoretical arguments imply an advantage to rare types, suggesting the number of types should grow consistently; however, empirical observations are very different. While some isogamous species exhibit thousands of mating types, such species are exceedingly rare, and most have fewer than 10. In this paper, we present a mathematical analysis to quantify the role of fitness variation-characterized by different mortality rates-in determining the number of mating types emerging in simple evolutionary models. We predict that the number of mating types decreases as the variance of mortality increases.

The Evolution of Sexual Reproduction and the Mating-Type Locus: Links to Pathogenesis of Cryptococcus Human Pathogenic Fungi

Cryptococcus species utilize a variety of sexual reproduction mechanisms, which generate genetic diversity, purge deleterious mutations, and contribute to their ability to occupy myriad environmental niches and exhibit a range of pathogenic potential. The bisexual and unisexual cycles of pathogenic Cryptococcus species are stimulated by properties associated with their environmental niches and proceed through well-characterized signaling pathways and corresponding morphological changes. Genes governing mating are encoded by the mating-type (MAT) loci and influence pathogenesis, population dynamics, and lineage divergence in Cryptococcus. MAT has undergone significant evolutionary changes within the Cryptococcus genus, including transition from the ancestral tetrapolar state in nonpathogenic species to a bipolar mating system in pathogenic species, as well as several internal reconfigurations. Owing to the variety of established sexual reproduction mechanisms and the robust characterization of the evolution of mating and MAT in this genus, Cryptococcus species provide key insights into the evolution of sexual reproduction.

Cause and Effectors: Whole-Genome Comparisons Reveal Shared but Rapidly Evolving Effector Sets among Host-Specific Plant-Castrating Fungi

Plant pathogens use molecular weapons to successfully infect their hosts, secreting a large portfolio of various proteins and enzymes. Different plant species are often parasitized by host-specific pathogens; however, it is still unclear whether the molecular basis of such host specialization involves species-specific weapons or different variants of the same weapons. We therefore compared the genes encoding secreted proteins in three plant-castrating pathogens parasitizing different host plants, producing their spores in plant anthers by replacing pollen. We validated our predictions for secretion signals for some genes and checked that our predicted secreted proteins were often highly expressed during plant infection. While we found few species-specific secreted proteins, numerous genes encoding secreted proteins showed signs of rapid evolution and of natural selection. Our study thus found that most changes among closely related host-specific pathogens involved rapid adaptive changes in shared molecular weapons rather than innovations for new weapons.

Plant pathogens utilize a portfolio of secreted effectors to successfully infect and manipulate their hosts. It is, however, still unclear whether changes in secretomes leading to host specialization involve mostly effector gene gains/losses or changes in their sequences. To test these hypotheses, we compared the secretomes of three host-specific castrating anther smut fungi (Microbotryum), two being sister species. To address within-species evolution, which might involve coevolution and local adaptation, we compared the secretomes of strains from differentiated populations. We experimentally validated a subset of signal peptides. Secretomes ranged from 321 to 445 predicted secreted proteins (SPs), including a few species-specific proteins (42 to 75), and limited copy number variation, i.e., little gene family expansion or reduction. Between 52% and 68% of the SPs did not match any Pfam domain, a percentage that reached 80% for the small secreted proteins, indicating rapid evolution. In comparison to background genes, we indeed found SPs to be more differentiated among species and strains, more often under positive selection, and highly expressed in planta; repeat-induced point mutations (RIPs) had no role in effector diversification, as SPs were not closer to transposable elements than background genes and were not more RIP affected. Our study thus identified both conserved core proteins, likely required for the pathogenic life cycle of all Microbotryum species, and proteins that were species specific or evolving under positive selection; these proteins may be involved in host specialization and/or coevolution. Most changes among closely related host-specific pathogens, however, involved rapid changes in sequences rather than gene gains/losses.

Consequences of recombination for the evolution of the mating type locus in Chlamydomonas reinhardtii

Recombination suppression in sex chromosomes and mating type loci can lead to degeneration as a result of reduced selection efficacy and Muller's ratchet effects. However, genetic exchange in the form of noncrossover gene conversions may still take place within crossover-suppressed regions. Recent work has found evidence that gene conversion may explain the low degrees of allelic differentiation in the dimorphic mating-type locus (MT) of the isogamous alga Chlamydomonas reinhardtii. However, no one has tested whether gene conversion is sufficient to avoid the degeneration of functional sequence within MT. Here, we calculate degree of linkage disequilibrium (LD) across MT as a proxy for recombination rate and investigate its relationship to patterns of population genetic variation and the efficacy of selection in the region. We find that degree of LD predicts selection efficacy across MT, and that purifying selection is stronger in shared genes than in MT-limited genes to the point of being equivalent to that of autosomal genes. We argue that while crossover suppression is needed in the mating-type loci of many isogamous systems, these loci are less likely to experience selection to differentiate further. Thus, recombination can act in these regions and prevent degeneration caused by Hill-Robertson effects.

Convergent evolution of linked mating-type loci in basidiomycete fungi

Sexual development is a key evolutionary innovation of eukaryotes. In many species, mating involves interaction between compatible mating partners that can undergo cell and nuclear fusion and subsequent steps of development including meiosis. Mating compatibility in fungi is governed by the mating type (MAT) loci. In basidiomycetes, the ancestral state is hypothesized to be tetrapolar, with two genetically unlinked MAT loci containing homeodomain transcription factor genes (HD locus) and pheromone and pheromone receptor genes (P/R locus), respectively. Alleles at both loci must differ between mating partners for completion of sexual development. However, there are also basidiomycetes with bipolar mating systems, which can arise through genomic linkage of the HD and P/R loci. In the order Tremellales, bipolarity is found only in the pathogenic Cryptococcus species. Here, we describe the analysis of MAT loci from 24 species of the Trichosporonales, a sister order to the Tremellales. In all of the species analyzed, the MAT loci are fused and a single HD gene is present in each mating type, similar to the organization in the pathogenic Cryptococci. However, the HD and P/R allele combinations in the Trichosporonales are different from those in the pathogenic Cryptococci. This and the existence of tetrapolar species in the Tremellales suggest that fusion of the HD and P/R loci occurred independently in the Trichosporonales and pathogenic Cryptococci, supporting the hypothesis of convergent evolution towards fused MAT regions, similar to previous findings in other fungal groups. Unlike the fused MAT loci in several other basidiomycete lineages though, the gene content and gene order within the fused MAT loci are highly conserved in the Trichosporonales, and there is no apparent suppression of recombination extending from the MAT loci to adjacent chromosomal regions, suggesting different mechanisms for the evolution of physically linked MAT loci in these groups.

Evolution of asymmetric gamete signaling and suppressed recombination at the mating type locus

The two partners required for sexual reproduction are rarely the same. This pattern extends to species which lack sexual dimorphism yet possess self-incompatible gametes determined at mating-type regions of suppressed recombination, likely precursors of sex chromosomes. Here we investigate the role of cellular signaling in the evolution of mating-types. We develop a model of ligand-receptor dynamics, and identify factors that determine the capacity of cells to send and receive signals. The model specifies conditions favoring the evolution of gametes producing ligand and receptor asymmetrically and shows how these are affected by recombination. When the recombination rate evolves, the conditions favoring asymmetric signaling also favor tight linkage of ligand and receptor loci in distinct linkage groups. These results suggest that selection for asymmetric gamete signaling could be the first step in the evolution of non-recombinant mating-type loci, paving the road for the evolution of anisogamy and sexes.

Sexual reproduction, from birds to bees, relies on distinct classes of sex cells, known as gametes, fusing together. Most single cell organisms, rather than producing eggs and sperm, have similar sized gametes that fall into distinct ‘mating types’. However, only sex cells belonging to different mating types can fuse together and sexually reproduce.

At first glance, it seems illogical that cells from the same mating type cannot reproduce with each other, as this restricts eligible partners within a population and makes finding a mate more difficult. Yet the typical pattern amongst single cell organisms is still two distinct classes of sex cells, just as in birds and bees. How did the obsession with mating between two different types become favored during evolution?

One possibility is that cells with different mating types can recognize and communicate with each other more easily. Cells communicate by releasing proteins known as ligands, which bind to specific receptors found on the cell’s surface. Using mathematical modelling, Hadjivasiliou and Pomiankowski showed that natural selection typically favors ‘asymmetric’ signaling, whereby cells evolve to either produce receptor A with ligand B, or have the reverse pattern and produce receptor B with ligand A. These asymmetric mutants are favored because they avoid producing ligands that clog or activate the receptors on their own surface. As a result, different types of cells are better at recognizing each other and mating more quickly.

When cells sexually reproduce they exchange genetic material with each other to produce offspring with a combination of genes that differ to their own. However, if the genes coding for ligand and receptor pairs were constantly being ‘swapped’, this could lead to new combinations, and a loss of asymmetric signaling. Hadjivasiliou and Pomiankowski showed that for asymmetric signaling to evolve, natural selection favors the genes encoding these non-compatible ligand and receptor pairs to be closely linked within the genome. This ensures that the mis-matching ligand and receptor are inherited together, preventing cells from producing pairs which can bind to themselves.

This study provides an original way to address an evolutionary question which has long puzzled biologists. These findings raise further questions about how gametes evolved to become the sperm and egg, and how factors such as signaling between cells can determine the sex of more complex organisms, such as ourselves.

Understanding Adaptation, Coevolution, Host Specialization, and Mating System in Castrating Anther-Smut Fungi by Combining Population and Comparative Genomics

Anther-smut fungi provide a powerful system to study host-pathogen specialization and coevolution, with hundreds of Microbotryum species specialized on diverse Caryophyllaceae plants, castrating their hosts through manipulation of the hosts' reproductive organs to facilitate disease transmission. Microbotryum fungi have exceptional genomic characteristics, including dimorphic mating-type chromosomes, that make this genus anexcellent model for studying the evolution of mating systems and their influence on population genetics structure and adaptive potential. Important insights into adaptation, coevolution, host specialization, and mating system evolution have been gained using anther-smut fungi, with new insights made possible by the recent advent of genomic approaches. We illustrate with Microbotryum case studies how using a combination of comparative genomics, population genomics, and transcriptomics approaches enables the integration of different evolutionary perspectives across different timescales. We also highlight current challenges and suggest future studies that will contribute to advancing our understanding of the mechanisms underlying adaptive processes in populations of fungal pathogens.