Shared forces of sex chromosome evolution in haploid-mating and diploid-mating organisms: Microbotryum violaceum and other model organisms

Early research on anther-smut disease: A fuller view of science?

The anther-smut host-pathogen system has provided extensive insights into the evolutionary ecology of disease resistance, transmission modes, host shifts, pathogen specialization, and disease evolution in metapopulations. It also has led to unexpected insights into sex ratio distorters, sex chromosome evolution, and transposable elements in fungi. In addition, anther-smut disease played a major role in Linnaeus' germ theory and the correspondence on parasitic castration between Darwin and Becker, one of the first female botanists. Here, we explicitly highlight some of the realities in the process of science, using an unusual autobiographical approach to describe how we came to collaborate on this system in the 1980s. Using perspectives from our different career stages, we present a surprising narrative that could not be deduced from merely reading the published papers. While our work was grounded in previous ecological and evolutionary theory, it was the product as much of empirical failures and intellectual roadblocks, as the result of a progressive scientific method. Our experiences illustrate not only the "human dimension of science" but more importantly show that linear sequences of hypothesis testing do not necessarily lead to new study systems and new ideas. We suggest there is a need to re-evaluate the scientific method in ecology and evolution, especially where the challenge is to engage in a productive dialog between natural history and theory.

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.

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.

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.

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.

Meiotic recombination in the offspring of Microbotryum hybrids and its impact on pathogenicity

BACKGROUND

Hybridization is a central mechanism in evolution, producing new species or introducing important genetic variation into existing species. In plant-pathogenic fungi, adaptation and specialization to exploit a host species are key determinants of evolutionary success. Here, we performed experimental crosses between the two pathogenic Microbotryum species, M. lychnidis-dioicae and M. silenes-acaulis that are specialized to different hosts. The resulting offspring were analyzed on phenotypic and genomic levels to describe genomic characteristics of hybrid offspring and genetic factors likely involved in host-specialization.

RESULTS

Genomic analyses of interspecific fungal hybrids revealed that individuals were most viable if the majority of loci were inherited from one species. Interestingly, species-specific loci were strictly controlled by the species' origin of the mating type locus. Moreover we detected signs of crossing over and chromosome duplications in the genomes of the analyzed hybrids. In Microbotryum, mitochondrial DNA was found to be uniparentally inherited from the a2 mating type. Genome comparison revealed that most gene families are shared and the majority of genes are conserved between the two species, indicating very similar biological features, including infection and pathogenicity processes. Moreover, we detected 211 candidate genes that were retained under host-driven selection of backcrossed lines. These genes and might therefore either play a crucial role in host specialization or be linked to genes that are essential for specialization.

CONCLUSION

The combination of genome analyses with experimental selection and hybridization is a promising way to investigate host-pathogen interactions. This study manifests genetic factors of host specialization that are required for successful biotrophic infection of the post-zygotic stage, but also demonstrates the strong influence of intra-genomic conflicts or instabilities on the viability of hybrids in the haploid host-independent stage.

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.

Comparative whole-genome analysis reveals artificial selection effects on Ustilago esculenta genome

Ustilago esculenta, infects Zizania latifolia, and induced host stem swollen to be a popular vegetable called Jiaobai in China. It is the long-standing artificial selection that maximizes the occurrence of favourable Jiaobai, and thus maintaining the plant-fungi interaction and modulating the fungus evolving from plant pathogen to entophyte. In this study, whole genome of U. esculenta was sequenced and transcriptomes of the fungi and its host were analysed. The 20.2 Mb U. esculenta draft genome of 6,654 predicted genes including mating, primary metabolism, secreted proteins, shared a high similarity to related Smut fungi. But U. esculenta prefers RNA silencing not repeat-induced point in defence and has more introns per gene, indicating relatively slow evolution rate. The fungus also lacks some genes in amino acid biosynthesis pathway which were filled by up-regulated host genes and developed distinct amino acid response mechanism to balance the infection-resistance interaction. Besides, U. esculenta lost some surface sensors, important virulence factors and host range-related effectors to maintain the economic endophytic life. The elucidation of the U. esculenta genomic information as well as expression profiles can not only contribute to more comprehensive insights into the molecular mechanism underlying artificial selection but also into smut fungi-host interactions.

Multiple convergent supergene evolution events in mating-type chromosomes

Convergent adaptation provides unique insights into the predictability of evolution and ultimately into processes of biological diversification. Supergenes (beneficial gene linkage) are striking examples of adaptation, but little is known about their prevalence or evolution. A recent study on anther-smut fungi documented supergene formation by rearrangements linking two key mating-type loci, controlling pre- and post-mating compatibility. Here further high-quality genome assemblies reveal four additional independent cases of chromosomal rearrangements leading to regions of suppressed recombination linking these mating-type loci in closely related species. Such convergent transitions in genomic architecture of mating-type determination indicate strong selection favoring linkage of mating-type loci into cosegregating supergenes. We find independent evolutionary strata (stepwise recombination suppression) in several species, with extensive rearrangements, gene losses, and transposable element accumulation. We thus show remarkable convergence in mating-type chromosome evolution, recurrent supergene formation, and repeated evolution of similar phenotypes through different genomic changes.

Supergenes result from beneficial linkage and recombination suppression between two or more genes. Giraud and colleagues use whole genome sequencing data to show convergent evolution of supergenes on mating-type chromosomes in multiple closely-related fungal lineages.

Massive Expansion of Gypsy-Like Retrotransposons in Microbotryum Fungi

Transposable elements (TEs) are selfish, autonomously replicating DNA sequences that constitute a major component of eukaryotic genomes and contribute to genome evolution through their movement and amplification. Many fungal genomes, including the anther-smut fungi in the basidiomycete genus Microbotryum, have genome defense mechanisms, such as repeat-induced point mutation (RIP), which hypermutate repetitive DNA and limit TE activity. Little is known about how hypermutation affects the tempo of TE activity and their sequence evolution. Here we report the identification of a massive burst-like expansion of Gypsy-like retrotransposons in a strain of Microbotryum. This TE expansion evidently occurred in the face of RIP-like hypermutation activity. By examining the fitness of individual TE insertion variants, we found that RIP-like mutations impair TE fitness and limit proliferation. Our results provide evidence for a punctuated pattern of TE expansion in a fungal genome, similar to that observed in animals and plants. While targeted hypermutation is often thought of as an effective protection against mobile element activity, our findings suggest that active TEs can persist and undergo selection while they proliferate in genomes that have RIP-like defenses.

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

Sex chromosomes can display successive steps of recombination suppression known as "evolutionary strata," which are thought to result from the successive linkage of sexually antagonistic genes to sex-determining genes. However, there is little evidence to support this explanation. Here we investigate whether evolutionary strata can evolve without sexual antagonism using fungi that display suppressed recombination extending beyond loci determining mating compatibility despite lack of male/female roles associated with their mating types. By comparing full-length chromosome assemblies from five anther-smut fungi with or without recombination suppression in their mating-type chromosomes, we inferred the ancestral gene order and derived chromosomal arrangements in this group. This approach shed light on the chromosomal fusion underlying the linkage of mating-type loci in fungi and provided evidence for multiple clearly resolved evolutionary strata over a range of ages (0.9-2.1 million years) in mating-type chromosomes. Several evolutionary strata did not include genes involved in mating-type determination. The existence of strata devoid of mating-type genes, despite the lack of sexual antagonism, calls for a unified theory of sex-related chromosome evolution, incorporating, for example, the influence of partially linked deleterious mutations and the maintenance of neutral rearrangement polymorphism due to balancing selection on sexes and mating types.

Fungal Sex: The Basidiomycota

Fungi of the Basidiomycota, representing major pathogen lineages and mushroom-forming species, exhibit diverse means to achieve sexual reproduction, with particularly varied mechanisms to determine compatibilities of haploid mating partners. For species that require mating between distinct genotypes, discrimination is usually based on both the reciprocal exchange of diffusible mating pheromones, rather than sexes, and the interactions of homeodomain protein signals after cell fusion. Both compatibility factors must be heterozygous in the product of mating, and genetic linkage relationships of the mating pheromone/receptor and homeodomain genes largely determine the complex patterns of mating-type variation. Independent segregation of the two compatibility factors can create four haploid mating genotypes from meiosis, referred to as tetrapolarity. This condition is thought to be ancestral to the basidiomycetes. Alternatively, cosegregation by linkage of the two mating factors, or in some cases the absence of the pheromone-based discrimination, yields only two mating types from meiosis, referred to as bipolarity. Several species are now known to have large and highly rearranged chromosomal regions linked to mating-type genes. At the population level, polymorphism of the mating-type genes is an exceptional aspect of some basidiomycete fungi, where selection under outcrossing for rare, intercompatible allelic variants is thought to be responsible for numbers of mating types that may reach several thousand. Advances in genome sequencing and assembly are yielding new insights by comparative approaches among and within basidiomycete species, with the promise to resolve the evolutionary origins and dynamics of mating compatibility genetics in this major eukaryotic lineage.

Transcriptional analysis of mating and pre-infection stages of the anther smut, Microbotryum lychnidis-dioicae

Microbotryum lychnidis-dioicae is an obligate biotrophic parasite of the wildflower species Silene latifolia. This dikaryotic fungus, commonly known as an anther smut, requires that haploid, yeast-like sporidia of opposite mating types fuse and differentiate into dikaryotic hyphae that penetrate host tissue as part of the fungal life cycle. Mating occurs under conditions of cool temperatures and limited nutrients. Further development requires host cues or chemical mimics, including a variety of lipids, e.g. phytols. To identify global changes in transcription associated with developmental shifts, RNA-Seq was conducted at several in vitro stages of fungal propagation, i.e. haploid cells grown independently on rich and nutrient-limited media, mated cells on nutrient-limited media as well as a time course of such mated cells exposed to phytol. Comparison of haploid cells grown under rich and nutrient-limited conditions identified classes of genes probably associated with general nutrient availability, including components of the RNAi machinery. Some gene enrichment patterns comparing the nutrient-limited and mated transcriptomes suggested gene expression changes associated with the mating programme (e.g. homeodomain binding proteins, secreted proteins, proteins unique to M. lychnidis-dioicae¸ multicopper oxidases and RhoGEFs). Analysis for phytol treatment compared with mated cells alone allowed identification of genes likely to be involved in the dikaryotic switch (e.g. oligopeptide transporters). Gene categories of particular note in all three conditions included those in the major facilitator superfamily, proteins containing PFAM domains of the secretory lipase family as well as proteins predicted to be secreted, many of which have the hallmarks of fungal effectors with potential roles in pathogenicity.

Resurgence of Less-Studied Smut Fungi as Models of Phytopathogenesis in the Omics Age

The smut fungi form a large, diverse, and nonmonophyletic group of plant pathogens that have long served as both important pests of human agriculture and, also, as fertile organisms of scientific investigation. As modern techniques of molecular genetic analysis became available, many previously studied species that proved refractive to these techniques fell by the wayside and were neglected. Now, as the advent of rapid and affordable next-generation sequencing provides genomic and transcriptomic resources for even these "forgotten" fungi, several species are making a comeback and retaking prominent places in phytopathogenic research. In this review, we highlight several of these smut fungi, with special emphasis on Microbotryum lychnidis-dioicae, an anther smut whose molecular genetic tools have finally begun to catch up with its historical importance in classical genetics and now provide mechanistic insights for ecological studies, evolution of host-pathogen interaction, and investigations of emerging infectious disease.

Reliable transformation system for Microbotryum lychnidis-dioicae informed by genome and transcriptome project

Microbotryum lychnidis-dioicae is a member of a species complex infecting host plants in the Caryophyllaceae. It is used as a model system in many areas of research, but attempts to make this organism tractable for reverse genetic approaches have not been fruitful. Here, we exploited the recently obtained genome sequence and transcriptome analysis to inform our design of constructs for use in Agrobacterium-mediated transformation techniques currently available for other fungi. Reproducible transformation was demonstrated at the genomic, transcriptional and functional levels. Moreover, these initial proof-of-principle experiments provide evidence that supports the findings from initial global transcriptome analysis regarding expression from the respective promoters under different growth conditions of the fungus. The technique thus provides for the first time the ability to stably introduce transgenes and over-express target M. lychnidis-dioicae genes.

Deciphering evolutionary strata on plant sex chromosomes and fungal mating-type chromosomes through compositional segmentation

Sex chromosomes have evolved from a pair of homologous autosomes which differentiated into sex determination systems, such as XY or ZW system, as a consequence of successive recombination suppression between the gametologous chromosomes. Identifying the regions of recombination suppression, namely, the "evolutionary strata", is central to understanding the history and dynamics of sex chromosome evolution. Evolution of sex chromosomes as a consequence of serial recombination suppressions is well-studied for mammals and birds, but not for plants, although 48 dioecious plants have already been reported. Only two plants Silene latifolia and papaya have been studied until now for the presence of evolutionary strata on their X chromosomes, made possible by the sequencing of sex-linked genes on both the X and Y chromosomes, which is a requirement of all current methods that determine stratum structure based on the comparison of gametologous sex chromosomes. To circumvent this limitation and detect strata even if only the sequence of sex chromosome in the homogametic sex (i.e. X or Z chromosome) is available, we have developed an integrated segmentation and clustering method. In application to gene sequences on the papaya X chromosome and protein-coding sequences on the S. latifolia X chromosome, our method could decipher all known evolutionary strata, as reported by previous studies. Our method, after validating on known strata on the papaya and S. latifolia X chromosome, was applied to the chromosome 19 of Populus trichocarpa, an incipient sex chromosome, deciphering two, yet unknown, evolutionary strata. In addition, we applied this approach to the recently sequenced sex chromosome V of the brown alga Ectocarpus sp. that has a haploid sex determination system (UV system) recovering the sex determining and pseudoautosomal regions, and then to the mating-type chromosomes of an anther-smut fungus Microbotryum lychnidis-dioicae predicting five strata in the non-recombining region of both the chromosomes.

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

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

Size Does Matter: Staging of Silene latifolia Floral Buds for Transcriptome Studies

Dioecious plants in the Caryophyllaceae family are susceptible to infection by members of the anthericolous smut fungi. In our studies of the Silene latifolia/Microbotryum lychnidis-dioicae pathosystem, we were interested in characterizing the plant-pathogen interaction at the molecular level before and during teliosporogenesis. This takes place during floral bud development, and we hoped to capture the interaction by Illumina Next-Gen RNA-Sequencing. Using previous literature that documented the stages of the floral buds for S. latifolia, we examined the floral buds from plants grown and infected under growth chamber conditions, using the disserting microscope to determine the stage of floral buds based on the morphology. We compiled the information and determined the size of floral buds that correspond to the desired stages of development for tissue collection, for the purpose of RNA-sequencing. This offers a practical approach for researchers who require a large number of floral buds/tissue categorized by stages of development, ascertaining whether infected/uninfected buds are at comparable stages of development and whether this also holds true for male vs. female buds. We also document our experience in infecting the plants and some of the unusual morphologies we observed after infection.

Sex and parasites: genomic and transcriptomic analysis of Microbotryum lychnidis-dioicae, the biotrophic and plant-castrating anther smut fungus

Background

The genus Microbotryum includes plant pathogenic fungi afflicting a wide variety of hosts with anther smut disease. Microbotryum lychnidis-dioicae infects Silene latifolia and replaces host pollen with fungal spores, exhibiting biotrophy and necrosis associated with altering plant development.

Results

We determined the haploid genome sequence for M. lychnidis-dioicae and analyzed whole transcriptome data from plant infections and other stages of the fungal lifecycle, revealing the inventory and expression level of genes that facilitate pathogenic growth. Compared to related fungi, an expanded number of major facilitator superfamily transporters and secretory lipases were detected; lipase gene expression was found to be altered by exposure to lipid compounds, which signaled a switch to dikaryotic, pathogenic growth. In addition, while enzymes to digest cellulose, xylan, xyloglucan, and highly substituted forms of pectin were absent, along with depletion of peroxidases and superoxide dismutases that protect the fungus from oxidative stress, the repertoire of glycosyltransferases and of enzymes that could manipulate host development has expanded. A total of 14 % of the genome was categorized as repetitive sequences. Transposable elements have accumulated in mating-type chromosomal regions and were also associated across the genome with gene clusters of small secreted proteins, which may mediate host interactions.

Conclusions

The unique absence of enzyme classes for plant cell wall degradation and maintenance of enzymes that break down components of pollen tubes and flowers provides a striking example of biotrophic host adaptation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1660-8) contains supplementary material, which is available to authorized users.

Chaos of Rearrangements in the Mating-Type Chromosomes of the Anther-Smut Fungus Microbotryum lychnidis-dioicae

Sex chromosomes in plants and animals and fungal mating-type chromosomes often show exceptional genome features, with extensive suppression of homologous recombination and cytological differentiation between members of the diploid chromosome pair. Despite strong interest in the genetics of these chromosomes, their large regions of suppressed recombination often are enriched in transposable elements and therefore can be challenging to assemble. Here we show that the latest improvements of the PacBio sequencing yield assembly of the whole genome of the anther-smut fungus, Microbotryum lychnidis-dioicae (the pathogenic fungus causing anther-smut disease of Silene latifolia), into finished chromosomes or chromosome arms, even for the repeat-rich mating-type chromosomes and centromeres. Suppressed recombination of the mating-type chromosomes is revealed to span nearly 90% of their lengths, with extreme levels of rearrangements, transposable element accumulation, and differentiation between the two mating types. We observed no correlation between allelic divergence and physical position in the nonrecombining regions of the mating-type chromosomes. This may result from gene conversion or from rearrangements of ancient evolutionary strata, i.e., successive steps of suppressed recombination. Centromeres were found to be composed mainly of copia-like transposable elements and to possess specific minisatellite repeats identical between the different chromosomes. We also identified subtelomeric motifs. In addition, extensive signs of degeneration were detected in the nonrecombining regions in the form of transposable element accumulation and of hundreds of gene losses on each mating-type chromosome. Furthermore, our study highlights the potential of the latest breakthrough PacBio chemistry to resolve complex genome architectures.

Degeneration of the nonrecombining regions in the mating-type chromosomes of the anther-smut fungi

Dimorphic mating-type chromosomes in fungi are excellent models for understanding the genomic consequences of recombination suppression. Their suppressed recombination and reduced effective population size are expected to limit the efficacy of natural selection, leading to genomic degeneration. Our aim was to identify the sequences of the mating-type chromosomes (a1 and a2) of the anther-smut fungi and to investigate degeneration in their nonrecombining regions. We used the haploid a1 Microbotryum lychnidis-dioicae reference genome sequence. The a1 and a2 mating-type chromosomes were both isolated electrophoretically and sequenced. Integration with restriction-digest optical maps identified regions of recombination and nonrecombination in the mating-type chromosomes. Genome sequence data were also obtained for 12 other Microbotryum species. We found strong evidence of degeneration across the genus in the nonrecombining regions of the mating-type chromosomes, with significantly higher rates of nonsynonymous substitution (dN/dS) than in nonmating-type chromosomes or in recombining regions of the mating-type chromosomes. The nonrecombining regions of the mating-type chromosomes also showed high transposable element content, weak gene expression, and gene losses. The levels of degeneration did not differ between the a1 and a2 mating-type chromosomes, consistent with the lack of homogametic/heterogametic asymmetry between them, and contrasting with X/Y or Z/W sex chromosomes.

The evolution of sex chromosomes in organisms with separate haploid sexes

The evolution of dimorphic sex chromosomes is driven largely by the evolution of reduced recombination and the subsequent accumulation of deleterious mutations. Although these processes are increasingly well understood in diploid organisms, the evolution of dimorphic sex chromosomes in haploid organisms (U/V) has been virtually unstudied theoretically. We analyze a model to investigate the evolution of linkage between fitness loci and the sex-determining region in U/V species. In a second step, we test how prone nonrecombining regions are to degeneration due to accumulation of deleterious mutations. Our modeling predicts that the decay of recombination on the sex chromosomes and the addition of strata via fusions will be just as much a part of the evolution of haploid sex chromosomes as in diploid sex chromosome systems. Reduced recombination is broadly favored, as long as there is some fitness difference between haploid males and females. The degeneration of the sex-determining region due to the accumulation of deleterious mutations is expected to be slower in haploid organisms because of the absence of masking. Nevertheless, balancing selection often drives greater differentiation between the U/V sex chromosomes than in X/Y and Z/W systems. We summarize empirical evidence for haploid sex chromosome evolution and discuss our predictions in light of these findings.

Recent gene-capture on the UV sex chromosomes of the moss Ceratodon purpureus

Sex chromosomes evolve from ordinary autosomes through the expansion and subsequent degeneration of a region of suppressed recombination that is inherited through one sex. Here we investigate the relative timing of these processes in the UV sex chromosomes of the moss Ceratodon purpureus using molecular population genetic analyses of eight newly discovered sex-linked loci. In this system, recombination is suppressed on both the female-transmitted (U) sex chromosome and the male-transmitted (V) chromosome. Genes on both chromosomes therefore should show the deleterious effects of suppressed recombination and sex-limited transmission, while purifying selection should maintain homologs of genes essential for both sexes on both sex chromosomes. Based on analyses of eight sex-linked loci, we show that the nonrecombining portions of the U and V chromosomes expanded in at least two events (~0.6-1.3 MYA and ~2.8-3.5 MYA), after the divergence of C. purpureus from its dioecious sister species, Trichodon cylindricus and Cheilothela chloropus. Both U- and V-linked copies showed reduced nucleotide diversity and limited population structure, compared to autosomal loci, suggesting that the sex chromosomes experienced more recent selective sweeps that the autosomes. Collectively these results highlight the dynamic nature of gene composition and molecular evolution on nonrecombining portions of the U and V sex chromosomes.

Extensive divergence between mating-type chromosomes of the anther-smut fungus

Genomic regions that determine mating compatibility are subject to distinct evolutionary forces that can lead to a cessation of meiotic recombination and the accumulation of structural changes between members of the homologous chromosome pair. The relatively recent discovery of dimorphic mating-type chromosomes in fungi can aid the understanding of sex chromosome evolution that is common to dioecious plants and animals. For the anther-smut fungus, Microbotryum lychnidis-dioicae (= M. violaceum isolated from Silene latifolia), the extent of recombination cessation on the dimorphic mating-type chromosomes has been conflictingly reported. Comparison of restriction digest optical maps for the two mating-type chromosomes shows that divergence extends over 90% of the chromosome lengths, flanked at either end by two pseudoautosomal regions. Evidence to support the expansion of recombination cessation in stages from the mating-type locus toward the pseudoautosomal regions was not found, but evidence of such expansion could be obscured by ongoing processes that affect genome structure. This study encourages the comparison of forces that may drive large-scale recombination suppression in fungi and other eukaryotes characterized by dimorphic chromosome pairs associated with sexual life cycles.

Linkage to the mating-type locus across the genus Microbotryum: insights into nonrecombining chromosomes

Parallels have been drawn between the evolution of nonrecombining regions in fungal mating-type chromosomes and animal and plant sex chromosomes, particularly regarding the stages of recombination cessation forming evolutionary strata of allelic divergence. Currently, evidence and explanations for recombination cessation in fungi are sparse, and the presence of evolutionary strata has been examined in a minimal number of fungal taxa. Here, the basidiomycete genus Microbotryum was used to determine the history of recombination cessation for loci on the mating-type chromosomes. Ancestry of linkage with mating type for 13 loci was assessed across 20 species by a phylogenetic method. No locus was found to exhibit trans-specific polymorphism for alternate alleles as old as the mating pheromone receptor, indicating that ages of linkage to mating type varied among the loci. The ordering of loci in the ancestry of linkage to mating type does not agree with their previously proposed assignments to evolutionary strata. This study suggests that processes capable of influencing divergence between alternate alleles may act at loci in the nonrecombining regions (e.g., gene conversion) and encourages further work to dissect the evolutionary processes acting upon genomic regions that determine mating compatibility.

Contrasted patterns of molecular evolution in dominant and recessive self-incompatibility haplotypes in Arabidopsis

Self-incompatibility has been considered by geneticists a model system for reproductive biology and balancing selection, but our understanding of the genetic basis and evolution of this molecular lock-and-key system has remained limited by the extreme level of sequence divergence among haplotypes, resulting in a lack of appropriate genomic sequences. In this study, we report and analyze the full sequence of eleven distinct haplotypes of the self-incompatibility locus (S-locus) in two closely related Arabidopsis species, obtained from individual BAC libraries. We use this extensive dataset to highlight sharply contrasted patterns of molecular evolution of each of the two genes controlling self-incompatibility themselves, as well as of the genomic region surrounding them. We find strong collinearity of the flanking regions among haplotypes on each side of the S-locus together with high levels of sequence similarity. In contrast, the S-locus region itself shows spectacularly deep gene genealogies, high variability in size and gene organization, as well as complete absence of sequence similarity in intergenic sequences and striking accumulation of transposable elements. Of particular interest, we demonstrate that dominant and recessive S-haplotypes experience sharply contrasted patterns of molecular evolution. Indeed, dominant haplotypes exhibit larger size and a much higher density of transposable elements, being matched only by that in the centromere. Overall, these properties highlight that the S-locus presents many striking similarities with other regions involved in the determination of mating-types, such as sex chromosomes in animals or in plants, or the mating-type locus in fungi and green algae.

Self-incompatibility is a common genetic system preventing selfing through recognition and rejection of self-pollen in hermaphroditic flowering plants. In the Brassicaceae family, this system is controlled by a single genomic region, called the S-locus, where many distinct specificities segregate in natural populations. In this study, we obtained genomic sequences comprising the S-locus in two closely related Brassicaceae species, Arabidopsis lyrata and A. halleri, and analyzed their diversity and patterns of molecular evolution. We report compelling evidence that the S-locus presents many similar properties with other genomic regions involved in the determination of mating-types in mammals, insects, plants, or fungi. In particular, in spite of their diversity, these genomic regions all show absence of similarity in intergenic sequences, large depth of genealogies, highly divergent organization, and accumulation of transposable elements. Moreover, some of these features were found to vary according to dominance of the S-locus specificities, suggesting that dominance/recessivity interactions are key drivers of the evolution of this genomic region.

DNA polymorphism in recombining and non-recombining mating-type-specific loci of the smut fungus Microbotryum

The population-genetic processes leading to the genetic degeneration of non-recombining regions have mainly been studied in animal and plant sex chromosomes. Here, we report population genetic analysis of the processes in the non-recombining mating-type-specific regions of the smut fungus Microbotryum violaceum. M. violaceum has A1 and A2 mating types, determined by mating-type-specific 'sex chromosomes' that contain 1-2 Mb long non-recombining regions. If genetic degeneration were occurring, then one would expect reduced DNA polymorphism in the non-recombining regions of this fungus. The analysis of DNA diversity among 19 M. violaceum strains, collected across Europe from Silene latifolia flowers, revealed that (i) DNA polymorphism is relatively low in all 20 studied loci (π∼0.15%), (ii) it is not significantly different between the two mating-type-specific chromosomes nor between the non-recombining and recombining regions, (iii) there is substantial population structure in M. violaceum populations, which resembles that of its host species, S. latifolia, and (iv) there is significant linkage disequilibrium, suggesting that widespread selfing in this species results in a reduction of the effective recombination rate across the genome. We hypothesise that selfing-related reduction of recombination across the M. violaceum genome negates the difference in the level of DNA polymorphism between the recombining and non-recombining regions, and may possibly lead to similar levels of genetic degeneration in the mating-type-specific regions of the non-recombining 'sex chromosomes' and elsewhere in the genome.

Consequences of reproductive mode on genome evolution in fungi

An organism's reproductive mode is believed to be a major factor driving its genome evolution. In theory, sexual inbreeding and asexuality are associated with lower effective recombination levels and smaller effective population sizes than sexual outbreeding, giving rise to reduced selection efficiency and genetic hitchhiking. This, in turn, is predicted to result in the accumulation of deleterious mutations and other genomic changes, for example the accumulation of repetitive elements. Empirical data from plants and animals supporting/refuting these theories are sparse and have yielded few conclusive results. A growing body of data from the fungal kingdom, wherein reproductive behavior varies extensively within and among taxonomic groups, has provided new insights into the role of mating systems (e.g., homothallism, heterothallism, pseudohomothallism) and asexuality, on genome evolution. Herein, we briefly review the theoretical relationships between reproductive mode and genome evolution and give examples of empirical data on the topic derived to date from plants and animals. We subsequently focus on the available data from fungi, which suggest that reproductive mode alters the rates and patterns of genome evolution in these organisms, e.g., protein evolution, mutation rate, codon usage, frequency of genome rearrangements and repetitive elements, and variation in chromosome size.

The evolution of sex: a perspective from the fungal kingdom

Sex is shrouded in mystery. Not only does it preferentially occur in the dark for both fungi and many animals, but evolutionary biologists continue to debate its benefits given costs in light of its pervasive nature. Experimental studies of the benefits and costs of sexual reproduction with fungi as model systems have begun to provide evidence that the balance between sexual and asexual reproduction shifts in response to selective pressures. Given their unique evolutionary history as opisthokonts, along with metazoans, fungi serve as exceptional models for the evolution of sex and sex-determining regions of the genome (the mating type locus) and for transitions that commonly occur between outcrossing/self-sterile and inbreeding/self-fertile modes of reproduction. We review here the state of the understanding of sex and its evolution in the fungal kingdom and also areas where the field has contributed and will continue to contribute to illuminating general principles and paradigms of sexual reproduction.

Evolutionary strata in a small mating-type-specific region of the smut fungus Microbotryum violaceum

DNA sequence analysis and genetic mapping of loci from mating-type-specific chromosomes of the smut fungus Microbotryum violaceum demonstrated that the nonrecombining mating-type-specific region in this species comprises approximately 25% ( approximately 1 Mb) of the chromosome length. Divergence between homologous mating-type-linked genes in this region varies between 0 and 8.6%, resembling the evolutionary strata of vertebrate and plant sex chromosomes.

Silene as a model system in ecology and evolution

The genus Silene, studied by Darwin, Mendel and other early scientists, is re-emerging as a system for studying interrelated questions in ecology, evolution and developmental biology. These questions include sex chromosome evolution, epigenetic control of sex expression, genomic conflict and speciation. Its well-studied interactions with the pathogen Microbotryum has made Silene a model for the evolution and dynamics of disease in natural systems, and its interactions with herbivores have increased our understanding of multi-trophic ecological processes and the evolution of invasiveness. Molecular tools are now providing new approaches to many of these classical yet unresolved problems, and new progress is being made through combining phylogenetic, genomic and molecular evolutionary studies with ecological and phenotypic data.

Hybrid sterility and inviability in the parasitic fungal species complex Microbotryum

Microbotryum violaceum, the anther-smut fungus, forms a complex of sibling species which specialize on different plants. Previous studies have shown the presence of partial ecological isolation and F1 inviability, but did not detect assortative mating apart from a high selfing rate. We investigated other post-mating barriers and show that F1 hybrid sterility, the inability of gametes to mate, increased gradually with the increasing genetic distance between the parents. F2 hybrids showed a reduced ability to infect the plants that was also correlated with the genetic distance. The host on which the F2 hybrids were passaged caused a selection for alleles derived from the pathogen species originally isolated from that host, but this effect was not detectable for the most closely related species. The post-mating barriers thus remain weak among the closest species pairs, suggesting that premating barriers are sufficient to initiate divergence in this system.

Ancient trans-specific polymorphism at pheromone receptor genes in basidiomycetes

In the majority of sexual organisms, reproduction occurs almost exclusively through the combination of distinct and alternate forms, called sexes or mating types. In some fungi, there can be dozens to hundreds of alternate alleles that determine compatible mating types. Such extensive polymorphism is expected to be maintained by balancing selection, and in extreme cases may give rise to trans-specific polymorphism. Here, we analyzed sequences of two pheromone receptors in the Microbotryum fungal species complex (Basidiomycota), which has only two alternate mating types. Several lines of evidence strongly suggest that the pheromone receptors are two allelic sequences acting to determine the alternate A1 and A2 mating types required for mating in Microbotryum. Phylogenetic trees of pheromone receptors in the Microbotryum species complex indicated a trans-specific polymorphism: the Microbotryum sequences from a given mating type were all more similar to the pheromone receptors of distantly related classes of fungi than to the alternate pheromone receptor in the Microbotryum species. A phylogenetic tree built using other known pheromone receptors from basidiomycetes showed that trans-specific polymorphism is widespread. The pheromone receptor alleles from Microbotryum appeared as the oldest, being at least 370 million years old. This represents the oldest known trans-specific polymorphism known in any organism so far, which may be due to the existence of sex chromosomes, obligate sexuality, mitochondrial inheritance linked to the mating type, and a highly selfing mating system in Microbotryum.

Identification of mating type genes in the bipolar basidiomycetous yeast Rhodosporidium toruloides: first insight into the MAT locus structure of the Sporidiobolales

Rhodosporidium toruloides is a heterothallic, bipolar, red yeast that belongs to the Sporidiobolales, an order within a major lineage of basidiomycetes, the Pucciniomycotina. In contrast to other basidiomycetes, considerably less is known about the nature of the mating type (MAT) loci that control sexual reproduction in this lineage. Three genes (RHA1, RHA2, and RHA3) encoding precursors of the MAT A1 pheromone (rhodotorucine A) were previously identified and formed the basis for a genome walking approach that led to the identification of additional MAT genes in complementary mating strains of R. toruloides. Two mating type-specific alleles encoding a p21-activated kinase (PAK; Ste20 homolog) were found between the RHA2 and RHA3 genes, and identification in MAT A2 strains of a gene encoding a presumptive pheromone precursor enabled prediction of the structure of rhodotorucine a. In addition, a putative pheromone receptor gene (STE3 homolog) was identified upstream of RHA1. Analyses of genomic data from two closely related species, Sporobolomyces roseus and Sporidiobolus salmonicolor, identified syntenic regions that contain homologs of all the above-mentioned genes. Notably, six novel pheromone precursor genes were uncovered, which encoded, similarly to the RHA genes, multiple tandem copies of the peptide moiety. This suggests that this structure, which is unique among fungal lipopeptide pheromones, seems to be prevalent in red yeasts. Species comparisons provided evidence for a large, multigenic MAT locus structure in the Sporidiobolales, but no putative homeodomain transcription factor genes (which are present in all basidiomycetous MAT loci characterized thus far) could be found in any of the three species in the vicinity of the MAT genes identified.

Expressed sequences tags of the anther smut fungus, Microbotryum violaceum, identify mating and pathogenicity genes

Background

The basidiomycete fungus Microbotryum violaceum is responsible for the anther-smut disease in many plants of the Caryophyllaceae family and is a model in genetics and evolutionary biology. Infection is initiated by dikaryotic hyphae produced after the conjugation of two haploid sporidia of opposite mating type. This study describes M. violaceum ESTs corresponding to nuclear genes expressed during conjugation and early hyphal production.

Results

A normalized cDNA library generated 24,128 sequences, which were assembled into 7,765 unique genes; 25.2% of them displayed significant similarity to annotated proteins from other organisms, 74.3% a weak similarity to the same set of known proteins, and 0.5% were orphans. We identified putative pheromone receptors and genes that in other fungi are involved in the mating process. We also identified many sequences similar to genes known to be involved in pathogenicity in other fungi. The M. violaceum EST database, MICROBASE, is available on the Web and provides access to the sequences, assembled contigs, annotations and programs to compare similarities against MICROBASE.

Conclusion

This study provides a basis for cloning the mating type locus, for further investigation of pathogenicity genes in the anther smut fungi, and for comparative genomics.

A linkage map reveals a complex basis for segregation distortion in an interpopulation cross in the moss Ceratodon purpureus

We report the construction of a linkage map for the moss Ceratodon purpureus (n = 13), based on a cross between geographically distant populations, and provide the first experimental confirmation of maternal chloroplast inheritance in bryophytes. From a mapping population of 288 recombinant haploid gametophytes, genotyped at 121 polymorphic AFLP loci, three gene-based nuclear loci, one chloroplast marker, and sex, we resolved 15 linkage groups resulting in a map length of approximately 730 cM. We estimate that the map covers more than three-quarters of the C. purpureus genome. Approximately 35% of the loci were sex linked, not including those in recombining pseudoautosomal regions. Nearly 45% of the loci exhibited significant segregation distortion (alpha = 0.05). Several pairs of unlinked distorted loci showed significant deviations from multiplicative genotypic frequencies, suggesting that distortion arises from genetic interactions among loci. The distorted autosomal loci all exhibited an excess of the maternal allele, suggesting that these interactions may involve nuclear-cytoplasmic factors. The sex ratio of the progeny was significantly male biased, and the pattern of nonrandom associations among loci indicates that this results from interactions between the sex chromosomes. These results suggest that even in interpopulation crosses, multiple mechanisms act to influence segregation ratios.

Sex determination in Chlamydomonas

The sex-determination system of the unicellular green alga, Chlamydomonas reinhardtii, is governed by genes in the mating-type (MT) locus and entails additional genes located in autosomes. Gene expression is initiated by nitrogen starvation, and cells differentiate into plus or minus gametes within 6h. Reviewed is our current understanding of gametic differentiation and fertilization, initiation of zygote development, and the uniparental inheritance of organelle genomes.

Common sex-linked deleterious alleles in a plant parasitic fungus alter infection success but show no pleiotropic advantage

Microbotryum violaceum is a fungus that causes the sterilizing anther smut disease in Caryophyllaceae. Its diploid teliospores normally produce equal proportions of haploid sporidia of its two mating types. However natural populations contain high frequencies of individuals producing sporidia of only one mating type ('biased strains'). This mating type-ratio bias is caused by deleterious alleles at haploid phase ('haplo-lethals') linked to the mating type locus that can be transmitted only by intra-tetrad selfing. We used experimental inoculations to test some of the hypotheses proposed to explain the maintenance of haplo-lethals. We found a disadvantage of biased strains in infection ability and high intra-tetrad mating rates. Biased strains had no higher competitive ability nor shorter latency and their higher spore production per flower appeared insufficient to compensate their disadvantages. These findings were only consistent with the hypothesis that haplo-lethals are maintained under a metapopulation structure because of high intra-tetrad selfing rates, founder effects and selection at the population level.

Chromosomal sex-determining regions in animals, plants and fungi

The independent evolution of sex chromosomes in many eukaryotic species raises questions about the evolutionary forces that drive their formation. Recent advances in our understanding of these genomic structures in mammals in parallel with alternate models such as the monotremes, fish, dioecious plants, and fungi support the idea of a remarkable convergence in structure to form large, non-recombining regions with discrete evolutionary strata. The discovery that evolutionary events similar to those that have transpired in humans have also occurred during the formation of sex chromosomes in organisms as divergent as the plant Silene, the fungus Cryptococcus and the fish medaka highlights the importance of future studies in these systems. Such investigation will broaden our knowledge of the evolution and plasticity of these ubiquitous genomic features underlying sexual dimorphism and reproduction.

Blocked recombination along the mating-type chromosomes of Neurospora tetrasperma involves both structural heterozygosity and autosomal genes

The Neurospora tetrasperma mating-type chromosomes have been shown to be structurally heterozygous by reciprocal introgression of these chromosomes between N. tetrasperma and N. crassa. This structural heterozygosity correlates with both a previously described recombination block and cytologically visible unpaired chromosomes at pachytene. Genes on the autosomes are also implicated in blocking recombination.

Repeat-induced point mutation and the population structure of transposable elements in Microbotryum violaceum

Repeat-induced point mutation (RIP) is a genome defense in fungi that hypermutates repetitive DNA and is suggested to limit the accumulation of transposable elements. The genome of Microbotryum violaceum has a high density of transposable elements compared to other fungi, but there is also evidence of RIP activity. This is the first report of RIP in a basidiomycete and was obtained by sequencing multiple copies of the integrase gene of a copia-type transposable element and the helicase gene of a Helitron-type element. In M. violaceum, the targets for RIP mutations are the cytosine residues of TCG trinucleotide combinations. Although RIP is a linkage-dependent process that tends to increase the variation among repetitive sequences, a chromosome-specific substructuring was observed in the transposable element population. The observed chromosome-specific patterns are not consistent with RIP, but rather suggest an effect of gene conversion, which is also a linkage-dependent process but results in a homogenization of repeated sequences. Particular sequences were found more widely distributed within the genome than expected by chance and may reflect the recently active variants. Therefore, sequence variation of transposable elements in M. violaceum appears to be driven by selection for transposition ability in combination with the context-specific forces of the RIP and gene conversion.

Repetitive DNA in the automictic fungus Microbotryum violaceum

The small genomes of fungi are expected to have little repetitive content other than rDNA genes. Moreover, among asexual or highly selfing lineages, the diversity of repetitive elements is also expected to be very low. However, in the automictic fungus Microbotryum violaceum, a very large proportion of random DNA fragments from the autosomes and the fungal sex chromosomes are repetitive in nature, either as retrotransposon or helicase sequences. Among the retrotransposon sequences, examples were found from each major kind of elements, including copia, gypsy, and non-LTR sequences. The most numerous were copia-like elements, which are believed to be rare in fungi, particularly among basidiomycetes. The many helicase sequences appear to belong to the recently discovered Helitron type of transposable elements. Also, sequences that could not be identified as a known type of gene were also very repetitive within the database of random fragments from M. violaceum. The differentiated pair of fungal sex chromosomes and suppression of recombination may be the major forces determining the highly repetitive content in the small genome of M. violaceum.