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

Expanding the genetic toolbox of Rhodotorula toruloides by identification and validation of six novel promoters induced or repressed under nitrogen starvation

BACKGROUND

The non-conventional yeast Rhodotorula toruloides is an emerging host organism in biotechnology by merit of its natural capacity to accumulate high levels of carotenoids and intracellular storage lipids from a variety of carbon sources. While the number of genetic engineering strategies that employ R. toruloides is increasing, the lack of genetic tools available for modification of this yeast is still limiting strain development. For instance, several strong, constitutive R. toruloides promoters have been characterized, but to date, only five inducible promoters have been identified. Although nitrogen-limited cultivation conditions are commonly used to induce lipid accumulation in this yeast, no promoters regulated by nitrogen starvation have been described for R. toruloides.

RESULTS

In this study, we used a combination of genomics and transcriptomics methods to identify novel R. toruloides promoter sequences that are either inducible or repressible by nitrogen starvation. RNA sequencing was used to assess gene expression in the recently isolated strain R. toruloides BOT-A2 during exponential growth and during nitrogen starvation, when cultivated with either glucose or xylose as the carbon source. The genome of BOT-A2 was sequenced using a combination of long- and short-read sequencing and annotated with support of the RNAseq data. Differential expression analysis was used to identify genes with a |log2 fold change|≥ 2 when comparing their expression during nitrogen depletion to that during exponential growth. The promoter regions from 16 of these genes were evaluated for their ability to drive the expression of a fluorescent reporter gene. Three promoters that were clearly upregulated under nitrogen starvation and three that were downregulated were selected and further characterized. One promoter, derived from gene RTBOTA2_003877, was found to function like an on-off switch, as it was only upregulated under full nitrogen depletion and downregulated in the presence of the nitrogen source.

CONCLUSIONS

Six new R. toruloides promoters that were either upregulated or downregulated under nitrogen-starvation were identified. These substantially contribute to the available promoters when engineering this organism and are foreseen to be particularly useful for future engineering strategies requiring specific regulation of target genes in accordance with nitrogen availability.

De Novo Genome Sequencing and Assembly of Rhodosporidium toruloides Strain Δdao1e

Rhodosporidium toruloides is an important industrial strain with high oil production. Here, we report the de novo sequencing of R. toruloides strain Δdao1e, which was derived from R. toruloides strain ATCC 10657, by combining sequencing results from the Illumina and PacBio Sequel platforms.

Induction of resistance mechanisms in Rhodotorula toruloides for growth in sugarcane hydrolysate with high inhibitor content

The oleaginous yeast Rhodotorula toruloides is a potential lipid producer for biodiesel production. However, this yeast shows growth inhibition due to harmful compounds when cultivated in hemicellulose hydrolysate. Here, we present a comparative analysis of colony selection and heterologous adaptive laboratory enhancement (ALE) strategies for obtaining robust strains. We implemented these ALE strategies for R. toruloides in a culture medium containing sugarcane hemicellulose hydrolysate. Our comparison study showed that the strain obtained with heterogeneous ALE strategy (Rth) reached a µ_max of 55% higher than the parental strain. It also exhibited higher biomass production (6.51 g/l) and lipid content (60%). ALE with colony selection strategy (Rtc) had a fitness gain in terms of shortening of the lag phase (9 h) when compared to Rth and parental strain (11.67, 12.33 h, respectively). When cultivated in Eucalyptus urograndis hemicellulose hydrolysate, the Rth strain achieved a high lipid content, 64%. Kinetics studies showed a strong effect of acetic acid as a repressor of xylose consumption during R. toruloides cultivation.Key points• Distinct adaptive laboratory strategies resulted in strains with different physiologies.• Heterologous adaptive laboratory enhancement provided the best results (fitness gain of 55% in µmax).• The Rth strain achieved a lipid content of 64.3% during cultivation in eucalyptus hemicellulose hydrolysate.

Understanding and exploiting the fatty acid desaturation system in Rhodotorula toruloides

BACKGROUND

Rhodotorula toruloides is a robust producer of triacylglycerol owing to its fast growth rate and strong metabolic flux under conditions of high cell density fermentation. However, the molecular basis of fatty acid biosynthesis, desaturation and regulation remains elusive.

RESULTS

We present the molecular characterization of four fatty acid desaturase (FAD) genes in R. toruloides. Biosynthesis of oleic acid (OA) and palmitoleic acid (POA) was conferred by a single-copy ∆9 Fad (Ole1) as targeted deletion of which abolished the biosynthesis of all unsaturated fatty acids. Conversion of OA to linoleic acid (LA) and α-linolenic acid (ALA) was predominantly catalyzed by the bifunctional ∆12/∆15 Fad2. FAD4 was found to encode a trifunctional ∆9/∆12/∆15 FAD, playing important roles in lipid and biomass production as well as stress resistance. Furthermore, an abundantly transcribed OLE1-related gene, OLE2 encoding a 149-aa protein, was shown to regulate Ole1 regioselectivity. Like other fungi, the transcription of FAD genes was controlled by nitrogen levels and fatty acids in the medium. A conserved DNA motif, (T/C)(G/A)TTGCAGA(T/C)CCCAG, was demonstrated to mediate the transcription of OLE1 by POA/OA. The applications of these FAD genes were illustrated by engineering high-level production of OA and γ-linolenic acid (GLA).

CONCLUSION

Our work has gained novel insights on the transcriptional regulation of FAD genes, evolution of FAD enzymes and their roles in UFA biosynthesis, membrane stress resistance and, cell mass and total fatty acid production. Our findings should illuminate fatty acid metabolic engineering in R. toruloides and beyond.

Diversity and phylogeny of basidiomycetous yeasts from plant leaves and soil: Proposal of two new orders, three new families, eight new genera and one hundred and seven new species

Nearly 500 basidiomycetous yeast species were accepted in the latest edition of The Yeasts: A Taxonomic Study published in 2011. However, this number presents only the tip of the iceberg of yeast species diversity in nature. Possibly more than 99 % of yeast species, as is true for many groups of fungi, are yet unknown and await discovery. Over the past two decades nearly 200 unidentified isolates were obtained during a series of environmental surveys of yeasts in phyllosphere and soils, mainly from China. Among these isolates, 107 new species were identified based on the phylogenetic analyses of nuclear ribosomal DNA (rDNA) [D1/D2 domains of the large subunit (LSU), the small subunit (SSU), and the internal transcribed spacer region including the 5.8S rDNA (ITS)] and protein-coding genes [both subunits of DNA polymerase II (RPB1 and RPB2), the translation elongation factor 1-α (TEF1) and the mitochondrial gene cytochrome b (CYTB)], and physiological comparisons. Forty-six of these belong to 16 genera in the Tremellomycetes (Agaricomycotina). The other 61 are distributed in 26 genera in the Pucciniomycotina. Here we circumscribe eight new genera, three new families and two new orders based on the multi-locus phylogenetic analyses combined with the clustering optimisation analysis and the predicted similarity thresholds for yeasts and filamentous fungal delimitation at genus and higher ranks. Additionally, as a result of these analyses, three new combinations are proposed and 66 taxa are validated.

The Engineering Potential of Rhodosporidium toruloides as a Workhorse for Biotechnological Applications

Moving our society towards a bioeconomy requires efficient and sustainable microbial production of chemicals and fuels. Rhodotorula (Rhodosporidium) toruloides is a yeast that naturally synthesizes substantial amounts of specialty chemicals and has been recently engineered to (i) enhance its natural production of lipids and carotenoids, and (ii) produce novel industrially relevant compounds. The use of R. toruloides by companies and research groups has exponentially increased in recent years as a result of recent improvements in genetic engineering techniques and the availability of multiomics information on its genome and metabolism. This review focuses on recent engineering approaches in R. toruloides for bioproduction and explores its potential as a biotechnological chassis.

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.

Genome sequence of the oleaginous yeast Rhodotorula toruloides strain CGMCC 2.1609

Most eukaryotic oleaginous species are yeasts and among them the basidiomycete red yeast, Rhodotorula (Rhodosporidium) toruloides (Pucciniomycotina) is known to produce high quantities of lipids when grown in nitrogen-limiting media, and has potential for biodiesel production. The genome of the CGMCC 2.1609 strain of this oleaginous red yeast was sequenced using a hybrid of Roche 454 and Illumina technology generating 13 × coverage. The de novo assembly was carried out using MIRA and scaffolded using MAQ and BAMBUS. The sequencing and assembly resulted in 365 scaffolds with total genome size of 33.4 Mb. The complete genome sequence of this strain was deposited in GenBank and the accession number is LKER00000000. The annotation is available on Figshare (doi:10.6084/m9.figshare.4754251).

Variation in mate-recognition pheromones of the fungal genus Microbotryum

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

Evolution of Mating Systems in Basidiomycetes and the Genetic Architecture Underlying Mating-Type Determination in the Yeast Leucosporidium scottii

In most fungi, sexual reproduction is bipolar; that is, two alternate sets of genes at a single mating-type (MAT) locus determine two mating types. However, in the Basidiomycota, a unique (tetrapolar) reproductive system emerged in which sexual identity is governed by two unlinked MAT loci, each of which controls independent mechanisms of self/nonself recognition. Tetrapolar-to-bipolar transitions have occurred on multiple occasions in the Basidiomycota, resulting, for example, from linkage of the two MAT loci into a single inheritable unit. Nevertheless, owing to the scarcity of molecular data regarding tetrapolar systems in the earliest-branching lineage of the Basidiomycota (subphylum Pucciniomycotina), it is presently unclear if the last common ancestor was tetrapolar or bipolar. Here, we address this question, by investigating the mating system of the Pucciniomycotina yeast Leucosporidium scottii. Using whole-genome sequencing and chromoblot analysis, we discovered that sexual reproduction is governed by two physically unlinked gene clusters: a multiallelic homeodomain (HD) locus and a pheromone/receptor (P/R) locus that is biallelic, thereby dismissing the existence of a third P/R allele as proposed earlier. Allele distribution of both MAT genes in natural populations showed that the two loci were in strong linkage disequilibrium, but independent assortment of MAT alleles was observed in the meiotic progeny of a test cross. The sexual cycle produces fertile progeny with similar proportions of the four mating types, but approximately 2/3 of the progeny was found to be nonhaploid. Our study adds to others in reinforcing tetrapolarity as the ancestral state of all basidiomycetes.

Genome sequencing provides insight into the reproductive biology, nutritional mode and ploidy of the fern pathogen Mixia osmundae

Mixia osmundae (Basidiomycota, Pucciniomycotina) represents a monotypic class containing an unusual fern pathogen with incompletely understood biology. We sequenced and analyzed the genome of M. osmundae, focusing on genes that may provide some insight into its mode of pathogenicity and reproductive biology. Mixia osmundae has the smallest plant pathogenic basidiomycete genome sequenced to date, at 13.6 Mb, with very few repeats, high gene density, and relatively few significant gene family gains. The genome shows that the yeast state of M. osmundae is haploid and the lack of segregation of mating genes suggests that the spores produced on Osmunda spp. fronds are probably asexual. However, our finding of a complete complement of mating and meiosis genes suggests the capacity to undergo sexual reproduction. Analyses of carbohydrate active enzymes suggest that this fungus is a biotroph with the ability to break down several plant cell wall components. Analyses of publicly available sequence data show that other Mixia members may exist on other plant hosts and with a broader distribution than previously known.

Molecular characterization of KU70 and KU80 homologues and exploitation of a KU70-deficient mutant for improving gene deletion frequency in Rhodosporidium toruloides

Background

Rhodosporidium toruloides is a β-carotenoid accumulating, oleaginous yeast that has great biotechnological potential. The lack of reliable and efficient genetic manipulation tools have been a major hurdle blocking its adoption as a biotechnology platform.

Results

We report for the first time the development of a highly efficient targeted gene deletion method in R. toruloides ATCC 10657 via Agrobacterium tumefaciens-mediated transformation. To further improve targeting frequency, the KU70 and KU80 homologs in R. toruloides were isolated and characterized in detail. A KU70-deficient mutant (∆ku70e) generated with the hygromycin selection cassette removed by the Cre-loxP recombination system showed a dramatically improved targeted gene deletion frequency, with over 90% of the transformants being true knockouts when homology sequence length of at least 1 kb was used. Successful gene targeting could be made with homologous flanking sequences as short as 100 bp in the ∆ku70e strain. KU70 deficiency did not perturb cell growth although an elevated sensitivity to DNA mutagenic agents was observed. Compared to the other well-known oleaginous yeast, Yarrowia lipolytica, R. toruloides KU70/KU80 genes contain much higher density of introns and are the most GC-rich KU70/KU80 genes reported.

Conclusions

The KU70-deficient mutant generated herein was effective in improving gene deletion frequency and allowed shorter homology sequences to be used for gene targeting. It retained the key oleaginous and fast growing features of R. toruloides. The strain should facilitate both fundamental and applied studies in this important yeast, with the approaches taken here likely to be applicable in other species in subphylum Pucciniomycotina.

Biogenesis of the Saccharomyces cerevisiae pheromone a-factor, from yeast mating to human disease

The mating pheromone a-factor secreted by Saccharomyces cerevisiae is a farnesylated and carboxylmethylated peptide and is unusually hydrophobic compared to other extracellular signaling molecules. Mature a-factor is derived from a precursor with a C-terminal CAAX motif that directs a series of posttranslational reactions, including prenylation, endoproteolysis, and carboxylmethylation. Historically, a-factor has served as a valuable model for the discovery and functional analysis of CAAX-processing enzymes. In this review, we discuss the three modules comprising the a-factor biogenesis pathway: (i) the C-terminal CAAX-processing steps carried out by Ram1/Ram2, Ste24 or Rce1, and Ste14; (ii) two sequential N-terminal cleavage steps, mediated by Ste24 and Axl1; and (iii) export by a nonclassical mechanism, mediated by the ATP binding cassette (ABC) transporter Ste6. The small size and hydrophobicity of a-factor present both challenges and advantages for biochemical analysis, as discussed here. The enzymes involved in a-factor biogenesis are conserved from yeasts to mammals. Notably, studies of the zinc metalloprotease Ste24 in S. cerevisiae led to the discovery of its mammalian homolog ZMPSTE24, which cleaves the prenylated C-terminal tail of the nuclear scaffold protein lamin A. Mutations that alter ZMPSTE24 processing of lamin A in humans cause the premature-aging disease progeria and related progeroid disorders. Intriguingly, recent evidence suggests that the entire a-factor pathway, including all three biogenesis modules, may be used to produce a prenylated, secreted signaling molecule involved in germ cell migration in Drosophila. Thus, additional prenylated signaling molecules resembling a-factor, with as-yet-unknown roles in metazoan biology, may await discovery.

Genome sequence of the oleaginous red yeast Rhodosporidium toruloides MTCC 457

We report the de novo assembled 20.05-Mb draft genome of the red yeast Rhodosporidium toruloides MTCC 457, predicted to encode 5,993 proteins, 4 rRNAs, and 125 tRNAs. Proteins known to be unique to oleaginous fungi are present among the predicted proteins. The genome sequence will be valuable for molecular genetic analysis and manipulation of lipid accumulation in this yeast and for developing it as a potential host for biofuel production.

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.

Discovery of a modified tetrapolar sexual cycle in Cryptococcus amylolentus and the evolution of MAT in the Cryptococcus species complex

Sexual reproduction in fungi is governed by a specialized genomic region called the mating-type locus (MAT). The human fungal pathogenic and basidiomycetous yeast Cryptococcus neoformans has evolved a bipolar mating system (a, α) in which the MAT locus is unusually large (>100 kb) and encodes >20 genes including homeodomain (HD) and pheromone/receptor (P/R) genes. To understand how this unique bipolar mating system evolved, we investigated MAT in the closely related species Tsuchiyaea wingfieldii and Cryptococcus amylolentus and discovered two physically unlinked loci encoding the HD and P/R genes. Interestingly, the HD (B) locus sex-specific region is restricted (∼2 kb) and encodes two linked and divergently oriented homeodomain genes in contrast to the solo HD genes (SXI1α, SXI2a) of C. neoformans and Cryptococcus gattii. The P/R (A) locus contains the pheromone and pheromone receptor genes but has expanded considerably compared to other outgroup species (Cryptococcus heveanensis) and is linked to many of the genes also found in the MAT locus of the pathogenic Cryptococcus species. Our discovery of a heterothallic sexual cycle for C. amylolentus allowed us to establish the biological roles of the sex-determining regions. Matings between two strains of opposite mating-types (A1B1×A2B2) produced dikaryotic hyphae with fused clamp connections, basidia, and basidiospores. Genotyping progeny using markers linked and unlinked to MAT revealed that meiosis and uniparental mitochondrial inheritance occur during the sexual cycle of C. amylolentus. The sexual cycle is tetrapolar and produces fertile progeny of four mating-types (A1B1, A1B2, A2B1, and A2B2), but a high proportion of progeny are infertile, and fertility is biased towards one parental mating-type (A1B1). Our studies reveal insights into the plasticity and transitions in both mechanisms of sex determination (bipolar versus tetrapolar) and sexual reproduction (outcrossing versus inbreeding) with implications for similar evolutionary transitions and processes in fungi, plants, and animals.

Genome walking in eukaryotes

Genome walking is a molecular procedure for the direct identification of nucleotide sequences from purified genomes. The only requirement is the availability of a known nucleotide sequence from which to start. Several genome walking methods have been developed in the last 20 years, with continuous improvements added to the first basic strategies, including the recent coupling with next generation sequencing technologies. This review focuses on the use of genome walking strategies in several aspects of the study of eukaryotic genomes. In a first part, the analysis of the numerous strategies available is reported. The technical aspects involved in genome walking are particularly intriguing, also because they represent the synthesis of the talent, the fantasy and the intelligence of several scientists. Applications in which genome walking can be employed are systematically examined in the second part of the review, showing the large potentiality of this technique, including not only the simple identification of nucleotide sequences but also the analysis of large collections of mutants obtained from the insertion of DNA of viral origin, transposons and transfer DNA (T-DNA) constructs. The enormous amount of data obtained indicates that genome walking, with its large range of applicability, multiplicity of strategies and recent developments, will continue to have much to offer for the rapid identification of unknown sequences in several fields of genomic research.

Is sex necessary?

Fungal sexual reproductive modes have markedly high diversity and plasticity, and asexual species have been hypothesized to arise frequently from sexual fungal species. A recent study on the red yeasts provides further support for the notion that sexual ancestors may give rise to shorter-lived asexual species. However, presumed asexual species may also be cryptically sexual, as revealed by other recent studies.

See research article: http://www.biomedcentral.com/1471-2148/11/249

Evidence for maintenance of sex determinants but not of sexual stages in red yeasts, a group of early diverged basidiomycetes

Background

The red yeasts are an early diverged group of basidiomycetes comprising sexual and asexual species. Sexuality is based on two compatible mating types and sexual identity is determined by MAT loci that encode homeodomain transcription factors, peptide pheromones and their receptors. The objective of the present study was to investigate the presence and integrity of MAT genes throughout the phylogenetic diversity of red yeasts belonging to the order Sporidiobolales.

Results

We surveyed 18 sexual heterothallic and self-fertile species and 16 asexual species. Functional pheromone receptor homologues (STE3.A1 and STE3.A2) were found in multiple isolates of most of the sexual and asexual species. For each of the two mating types, sequence comparisons with whole-genome data indicated that synteny tended to be conserved along the pheromone receptor region. For the homeodomain transcription factor, likelihood methods suggested that diversifying selection acting on the self/non-self recognition region promotes diversity in sexual species, while rapid evolution seems to be due to relaxed selection in asexual strains.

Conclusions

The majority of both sexual and asexual species of red yeasts have functional pheromone receptors and homeodomain homologues. This and the frequent existence of asexual strains within sexual species, makes the separation between sexual and asexual species imprecise. Events of loss of sexuality seem to be recent and frequent, but not uniformly distributed within the Sporidiobolales. Loss of sex could promote speciation by fostering the emergence of asexual lineages from an ancestral sexual stock, but does not seem to contribute to the generation of exclusively asexual lineages that persist for a long time.

Having sex, yes, but with whom? Inferences from fungi on the evolution of anisogamy and mating types

The advantage of sex has been among the most debated issues in biology. Surprisingly, the question of why sexual reproduction generally requires the combination of distinct gamete classes, such as small and large gametes, or gametes with different mating types, has been much less investigated. Why do systems with alternative gamete classes (i.e. systems with either anisogamy or mating types or both) appear even though they restrict the probability of finding a compatible mating partner? Why does the number of gamete classes vary from zero to thousands, with most often only two classes? We review here the hypotheses proposed to explain the origin, maintenance, number, and loss of gamete classes. We argue that fungi represent highly suitable models to help resolve issues related to the evolution of distinct gamete classes, because the number of mating types vary from zero to thousands across taxa, anisogamy is present or not, and because there are frequent transitions between these conditions. We review the nature and number of gamete classes in fungi, and we attempt to draw inferences from these data on the evolutionary forces responsible for their appearance, loss or maintenance, and number.

A deviation from the bipolar-tetrapolar mating paradigm in an early diverged basidiomycete

In fungi, sexual identity is determined by specialized genomic regions called MAT loci which are the equivalent to sex chromosomes in some animals and plants. Usually, only two sexes or mating types exist, which are determined by two alternate sets of genes (or alleles) at the MAT locus (bipolar system). However, in the phylum Basidiomycota, a unique tetrapolar system emerged in which four different mating types are generated per meiosis. This occurs because two functionally distinct molecular recognition systems, each encoded by one MAT region, constrain the selection of sexual partners. Heterozygosity at both MAT regions is a pre-requisite for mating in both bipolar and tetrapolar basidiomycetes. Tetrapolar mating behaviour results from the absence of genetic linkage between the two regions bringing forth up to thousands of mating types. The subphylum Pucciniomycotina, an early diverged lineage of basidiomycetes encompassing important plant pathogens such as the rusts and saprobes like Rhodosporidium and Sporidiobolus, has been so far poorly explored concerning the content and organization of MAT loci. Here we show that the red yeast Sporidiobolus salmonicolor has a mating system unlike any previously described because occasional disruptions of the genetic cohesion of the bipolar MAT locus originate new mating types. We confirmed that mating is normally bipolar and that heterozygosity at both MAT regions is required for mating. However, a laboratory cross showed that meiotic recombination may occur within the bipolar MAT locus, explaining tetrapolar features like increased allele number and evolution rates of some MAT genes. This pseudo-bipolar system deviates from the classical bipolar–tetrapolar paradigm and, to our knowledge, has never been observed before. We propose a model for MAT evolution in the Basidiomycota in which the pseudo-bipolar system may represent a hitherto unforeseen gradual form of transition from an ancestral tetrapolar system to bipolarity.

Sexual reproduction in fungi is regulated by genomic regions called MAT loci that determine sexual identity in a manner comparable to that of sex chromosomes in animal and plants. In most fungi, sexual reproduction is bipolar, i.e., two alternate and distinct sets of genes at the MAT locus determine two mating types (the equivalent to sexes). In the Basidiomycota, which is the fungal lineage that includes the mushrooms, a unique (tetrapolar) sexual reproduction system evolved in which the mating type is determined by two functionally independent and genetically unlinked regions. In addition, the tetrapolar system functions with multiple alleles for at least one the two classes of MAT genes. This potentially generates thousands of mating types, as observed for many mushroom species. Here we report on the molecular characterization of the mating system in the basidiomycetous red yeast Sporidiobolus salmonicolor, which belongs to the Pucciniomycotina, the earliest derived lineage of the Basidiomycota that has remained virtually unexplored with respect to gene content and structure of MAT loci. Our results revealed for the first time a mating system that is neither tetrapolar nor strictly bipolar, and we propose a model for the evolution of basidiomycete MAT loci that accommodates this novel finding.

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.

Orchestration of sexual reproduction and virulence by the fungal mating-type locus

The mating-type locus (MAT) orchestrates sexual reproduction in fungi. Sexual reproduction is related not only to fitness of an organism, but also correlated with virulence in certain pathogens. In the dandruff-associated fungus Malassesia globosa, although the sexual cycle remains to be discovered, whole genome analysis has led to the hypothesis that mating may occur on host skin. Furthermore, the MAT locus of M. globosa and U. hordei provides evidence that transitions between tetrapolar and bipolar systems have independently occurred. These results, together with studies recapitulating the ancestral tetrapolar mating system in Cryptococcus and the structure of MAT in related smut fungi, have furthered understanding on transitions between different mating systems and the evolution of MAT in the Basidiomycota.