Evolutionary restoration of fertility in an interspecies hybrid yeast, by whole-genome duplication after a failed mating-type switch

The promise and pitfalls of synteny in phylogenomics

Reconstructing the tree of life remains a central goal in biology. Early methods, which relied on small numbers of morphological or genetic characters, often yielded conflicting evolutionary histories, undermining confidence in the results. Investigations based on phylogenomics, which use hundreds to thousands of loci for phylogenetic inquiry, have provided a clearer picture of life's history, but certain branches remain problematic. To resolve difficult nodes on the tree of life, 2 recent studies tested the utility of synteny, the conserved collinearity of orthologous genetic loci in 2 or more organisms, for phylogenetics. Synteny exhibits compelling phylogenomic potential while also raising new challenges. This Essay identifies and discusses specific opportunities and challenges that bear on the value of synteny data and other rare genomic changes for phylogenomic studies. Synteny-based analyses of highly contiguous genome assemblies mark a new chapter in the phylogenomic era and the quest to reconstruct the tree of life.

Microtubule search-and-capture model evaluates the effect of chromosomal volume conservation on spindle assembly during mitosis

Variation in the chromosome numbers can arise from the erroneous mitosis or fusion and fission of chromosomes. While the mitotic errors lead to an increase or decrease in the overall chromosomal substance in the daughter cells, fission and fusion keep this conserved. Variations in chromosome numbers are assumed to be a crucial driver of speciation. For example, the members of the muntjac species are known to have very different karyotypes with the chromosome numbers varying from 2n=70+3B in the brown brocket deer to 2n=46 in the Chinese muntjac and 2n=6/7 in the Indian muntjac. The chromosomal content in the nucleus of these closely related mammals is roughly the same and various chromosome fusion and fission pathways have been suggested as the evolution process of these karyotypes. Similar trends can also be found in lepidoptera and yeast species which show a wide variation of chromosome numbers. The effect of chromosome number variation on the spindle assembly time and accuracy is still not properly addressed. We computationally investigate the effect of conservation of the total chromosomal substance on the spindle assembly during prometaphase. Our results suggest that chromosomal fusion pathways aid the microtubule-driven search and capture of the kinetochore in cells with monocentric chromosomes. We further report a comparative analysis of the site and percentage of amphitelic captures, dependence on cell shape, and position of the kinetochore in respect to chromosomal volume partitioning.

Cellular Responses and Targets in Food Spoilage Yeasts Exposed to Antifungal Prenylated Isoflavonoids

Prenylated isoflavonoids are phytochemicals with promising antifungal properties. Recently, it was shown that glabridin and wighteone disrupted the plasma membrane (PM) of the food spoilage yeast Zygosaccharomyces parabailii in distinct ways, which led us to investigate further their modes of action (MoA). Transcriptomic profiling with Z. parabailii showed that genes encoding transmembrane ATPase transporters, including Yor1, and genes homologous to the pleiotropic drug resistance (PDR) subfamily in Saccharomyces cerevisiae were upregulated in response to both compounds. Gene functions involved in fatty acid and lipid metabolism, proteostasis, and DNA replication processes were overrepresented among genes upregulated by glabridin and/or wighteone. Chemogenomic analysis using the genome-wide deletant collection for S. cerevisiae further suggested an important role for PM lipids and PM proteins. Deletants of gene functions involved in biosynthesis of very-long-chain fatty acids (constituents of PM sphingolipids) and ergosterol were hypersensitive to both compounds. Using lipid biosynthesis inhibitors, we corroborated roles for sphingolipids and ergosterol in prenylated isoflavonoid action. The PM ABC transporter Yor1 and Lem3-dependent flippases conferred sensitivity and resistance, respectively, to the compounds, suggesting an important role for PM phospholipid asymmetry in their MoAs. Impaired tryptophan availability, likely linked to perturbation of the PM tryptophan permease Tat2, was evident in response to glabridin. Finally, substantial evidence highlighted a role of the endoplasmic reticulum (ER) in cellular responses to wighteone, including gene functions associated with ER membrane stress or with phospholipid biosynthesis, the primary lipid of the ER membrane. IMPORTANCE Preservatives, such as sorbic acid and benzoic acid, inhibit the growth of undesirable yeast and molds in foods. Unfortunately, preservative tolerance and resistance in food spoilage yeast, such as Zygosaccharomyces parabailii, is a growing challenge in the food industry, which can compromise food safety and increase food waste. Prenylated isoflavonoids are the main defense phytochemicals in the Fabaceae family. Glabridin and wighteone belong to this group of compounds and have shown potent antifungal activity against food spoilage yeasts. The present study demonstrated the mode of action of these compounds against food spoilage yeasts by using advanced molecular tools. Overall, the cellular actions of these two prenylated isoflavonoids share similarities (at the level of the plasma membrane) but also differences. Tryptophan import was specifically affected by glabridin, whereas endoplasmic reticulum membrane stress was specifically induced by wighteone. Understanding the mode of action of these novel antifungal agents is essential for their application in food preservation.

Extreme genome diversity and cryptic speciation in a harmful algal-bloom-forming eukaryote

Harmful algal blooms of the toxic haptophyte Prymnesium parvum are a recurrent problem in many inland and estuarine waters around the world. Strains of P. parvum vary in the toxins they produce and in other physiological traits associated with harmful algal blooms, but the genetic basis for this variation is unknown. To investigate genome diversity in this morphospecies, we generated genome assemblies for 15 phylogenetically and geographically diverse strains of P. parvum, including Hi-C guided, near-chromosome-level assemblies for two strains. Comparative analysis revealed considerable DNA content variation between strains, ranging from 115 to 845 Mbp. Strains included haploids, diploids, and polyploids, but not all differences in DNA content were due to variation in genome copy number. Haploid genome size between strains of different chemotypes differed by as much as 243 Mbp. Syntenic and phylogenetic analyses indicate that UTEX 2797, a common laboratory strain from Texas, is a hybrid that retains two phylogenetically distinct haplotypes. Investigation of gene families variably present across the strains identified several functional categories associated with metabolic and genome size variation in P. parvum, including genes for the biosynthesis of toxic metabolites and proliferation of transposable elements. Together, our results indicate that P. parvum comprises multiple cryptic species. These genomes provide a robust phylogenetic and genomic framework for investigations into the eco-physiological consequences of the intra- and inter-specific genetic variation present in P. parvum and demonstrate the need for similar resources for other harmful algal-bloom-forming morphospecies.

Improving CRISPR-Cas9 mediated genome integration in interspecific hybrid yeasts

Saccharomyces pastorianus is not a classical taxon, it is an interspecific hybrid resulting from the cross of Saccharomyces cerevisiae and Saccharomyces eubayanus. Exhibiting heterosis for phenotypic traits such as wort α-oligosaccharide consumption and fermentation at low temperature, it has been domesticated to become the main workhorse of the brewing industry. Although CRISPR-Cas9 has been shown to be functional in S. pastorianus, repair of CRISPR- induced double strand break is unpredictable and preferentially uses the homoeologous chromosome as template, preventing targeted introduction of the desired repair construct. Here, we demonstrate that lager hybrids can be edited with near 100% efficiency at carefully selected landing sites on the chimeric SeScCHRIII. The landing sites were systematically selected and evaluated for (i) absence of loss of heterozygosity upon CRISPR-editing, (ii) efficiency of the gRNA, and (iii) absence of effect on strain physiology. Successful examples of highly efficient single and double gene integration illustrated that genome editing can be applied in interspecies hybrids, paving the way to a new impulse to lager yeast strain development. DATA AVAILABILITY: Data underlying graphs and figures found in this manuscript are deposited at the 4TU research dat center (https://data.4tu.nl/info/en/) and available through the doi: 10.4121/21648329.

Evolution of yeast hybrids by aborted meiosis

Sterile hybrids are broadly considered evolutionary dead-ends because of their faulty sexual reproduction. While sterility in obligate sexual organisms is a clear constraint in perpetuating the species, some facultative sexual microbes such as yeasts can propagate asexually and maintain genome plasticity. Moreover, incomplete meiotic pathways in yeasts represent alternative routes to the standard meiosis that generates genetic combinations in the population and fuel adaptation. Here, we review how aborting meiosis promotes genome-wide allele shuffling in sterile Saccharomyces hybrids and describe approaches to identify evolved clones in a cell population. We further discuss possible implications of this process in generating phenotypic novelty and report cases of abortive meiosis across yeast species.

Evolution and molecular bases of reproductive isolation

The most challenging problem in speciation research is disentangling the relative strength and order in which different reproductive barriers evolve. Here, we review recent developments in the study of reproductive isolation in yeasts. With over a thousand genome-sequenced isolates readily available for testing the viability, sterility, and fitness of both intraspecies and interspecies hybrid crosses, Saccharomyces yeasts are an ideal model to study such fundamental questions. Our survey demonstrates that, while chromosomal-level mutations are widespread at the intraspecific level, anti-recombination-driven chromosome missegregation is the primary reproductive barrier between species. Finally, despite their strength, all of these postzygotic barriers can be resolved through the asexual life history of hybrids.

Genome doubling enabled the expansion of yeast vesicle traffic pathways

Vesicle budding and fusion in eukaryotes depend on a suite of protein types, such as Arfs, Rabs, coats and SNAREs. Distinct paralogs of these proteins act at distinct intracellular locations, suggesting a link between gene duplication and the expansion of vesicle traffic pathways. Genome doubling, a common source of paralogous genes in fungi, provides an ideal setting in which to explore this link. Here we trace the fates of paralog doublets derived from the 100-Ma-old hybridization event that gave rise to the whole genome duplication clade of budding yeast. We find that paralog doublets involved in specific vesicle traffic functions and pathways are convergently retained across the entire clade. Vesicle coats and adaptors involved in secretory and early-endocytic pathways are retained as doublets, at rates several-fold higher than expected by chance. Proteins involved in later endocytic steps and intra-Golgi traffic, including the entire set of multi-subunit and coiled-coil tethers, have reverted to singletons. These patterns demonstrate that selection has acted to expand and diversify the yeast vesicle traffic apparatus, across species and time.

Mating-Type Switching in Budding Yeasts, from Flip/Flop Inversion to Cassette Mechanisms

Mating-type switching is a natural but unusual genetic control process that regulates cell identity in ascomycete yeasts. It involves physically replacing one small piece of genomic DNA by another, resulting in replacement of the master regulatory genes in the mating pathway and hence a switch of cell type and mating behavior. In this review, we concentrate on recent progress that has been made on understanding the origins and evolution of mating-type switching systems in budding yeasts (subphylum Saccharomycotina). Because of the unusual nature and the complexity of the mechanism in Saccharomyces cerevisiae, mating-type switching was assumed until recently to have originated only once or twice during yeast evolution. However, comparative genomics analysis now shows that switching mechanisms arose many times independently-at least 11 times in budding yeasts and once in fission yeasts-a dramatic example of convergent evolution. Most of these lineages switch mating types by a flip/flop mechanism that inverts a section of a chromosome and is simpler than the well-characterized 3-locus cassette mechanism (MAT/HML/HMR) used by S. cerevisiae. Mating-type switching (secondary homothallism) is one of the two possible mechanisms by which a yeast species can become self-fertile. The other mechanism (primary homothallism) has also emerged independently in multiple evolutionary lineages of budding yeasts, indicating that homothallism has been favored strongly by natural selection. Recent work shows that HO endonuclease, which makes the double-strand DNA break that initiates switching at the S. cerevisiae MAT locus, evolved from an unusual mobile genetic element that originally targeted a glycolytic gene, FBA1.

Multiple Hybridization Events Punctuate the Evolutionary Trajectory of Malassezia furfur

Malassezia species are important fungal skin commensals and are part of the normal microbiota of humans and other animals. However, under certain circumstances these fungi can also display a pathogenic behavior. For example, Malassezia furfur is a common commensal of human skin and yet is often responsible for skin disorders but also systemic infections. Comparative genomics analysis of M. furfur revealed that some isolates have a hybrid origin, similar to several other recently described hybrid fungal pathogens. Because hybrid species exhibit genomic plasticity that can impact phenotypes, we sought to elucidate the genomic evolution and phenotypic characteristics of M. furfur hybrids in comparison to their parental lineages. To this end, we performed a comparative genomics analysis between hybrid strains and their presumptive parental lineages and assessed phenotypic characteristics. Our results provide evidence that at least two distinct hybridization events occurred between the same parental lineages and that the parental strains may have originally been hybrids themselves. Analysis of the mating-type locus reveals that M. furfur has a pseudobipolar mating system and provides evidence that after sexual liaisons of mating compatible cells, hybridization involved cell-cell fusion leading to a diploid/aneuploid state. This study provides new insights into the evolutionary trajectory of M. furfur and contributes with valuable genomic resources for future pathogenicity studies. IMPORTANCEMalassezia furfur is a common commensal member of human/animal microbiota that is also associated with several pathogenic states. Recent studies report involvement of Malassezia species in Crohn's disease, a type of inflammatory bowel disease, pancreatic cancer progression, and exacerbation of cystic fibrosis. A recent genomics analysis of M. furfur revealed the existence of hybrid isolates and identified their putative parental lineages. In this study, we explored the genomic and phenotypic features of these hybrids in comparison to their putative parental lineages. Our results revealed the existence of a pseudobipolar mating system in this species and showed evidence for the occurrence of multiple hybridization events in the evolutionary trajectory of M. furfur. These findings significantly advance our understanding of the evolution of this commensal microbe and are relevant for future studies exploring the role of hybridization in the adaptation to new niches or environments, including the emergence of pathogenicity.

Genome analysis of five recently described species of the CUG-Ser clade uncovers Candida theae as a new hybrid lineage with pathogenic potential in the Candida parapsilosis species complex

Candida parapsilosis species complex comprises three important pathogenic species: Candida parapsilosis sensu stricto, Candida orthopsilosis and Candida metapsilosis. The majority of C. orthopsilosis and all C. metapsilosis isolates sequenced thus far are hybrids, and most of the parental lineages remain unidentified. This led to the hypothesis that hybrids with pathogenic potential were formed by the hybridization of non-pathogenic lineages that thrive in the environment. In a search for the missing hybrid parentals, and aiming to get a better understanding of the evolution of the species complex, we sequenced, assembled and analysed the genome of five close relatives isolated from the environment: Candida jiufengensis, Candida pseudojiufengensis, Candida oxycetoniae, Candida margitis and Candida theae. We found that the linear conformation of mitochondrial genomes in Candida species emerged multiple times independently. Furthermore, our analyses discarded the possible involvement of these species in the mentioned hybridizations, but identified C. theae as an additional hybrid in the species complex. Importantly, C. theae was recently associated with a case of infection, and we also uncovered the hybrid nature of this clinical isolate. Altogether, our results reinforce the hypothesis that hybridization is widespread among Candida species, and potentially contributes to the emergence of lineages with opportunistic pathogenic behaviour.

Molecular Tools for Leveraging the Potential of the Acid-Tolerant Yeast Zygosaccharomyces bailii as Cell Factory

Microorganisms offer a tremendous potential as cell factories, and they are indeed been used by humans since the previous centuries for biotransformations. Among them, yeasts combine the advantage of a unicellular state with a eukaryotic organization. Moreover, in the era of biorefineries, their biodiversity can offer solutions to specific process constraints. Zygosaccharomyces bailii, an ascomycete budding yeast, is widely known for its peculiar tolerance to different stresses, among which are organic acids. Moreover, the recent reclassification of the species, including diverse hybrids, is further expanding both fundamental and applied interests. It is therefore reasonable that despite the possibility to apply with this yeast some of the molecular tools and protocols routinely used to manipulate Saccharomyces cerevisiae, adjustments and optimizations are necessary. Here we describe in detail the methods for determining chromosome number, size, and aneuploidy, transformation, classical target gene disruption or gene integration, and designing of episomal expression plasmids helpful for engineering the yeast Z. bailii .

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

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

Trends in yeast diversity discovery

Yeasts, usually defined as unicellular fungi, occur in various fungal lineages. Hence, they are not a taxonomic unit, but rather represent a fungal lifestyle shared by several unrelated lineages. Although the discovery of new yeast species occurs at an increasing speed, at the current rate it will likely take hundreds of years, if ever, before they will all be documented. Many parts of the earth, including many threatened habitats, remain unsampled for yeasts and many others are only superficially studied. Cold habitats, such as glaciers, are home to a specific community of cold-adapted yeasts, and, hence, there is some urgency to study such environments at locations where they might disappear soon due to anthropogenic climate change. The same is true for yeast communities in various natural forests that are impacted by deforestation and forest conversion. Many countries of the so-called Global South have not been sampled for yeasts, despite their economic promise. However, extensive research activity in Asia, especially China, has yielded many taxonomic novelties. Comparative genomics studies have demonstrated the presence of yeast species with a hybrid origin, many of them isolated from clinical or industrial environments. DNA-metabarcoding studies have demonstrated the prevalence, and in some cases dominance, of yeast species in soils and marine waters worldwide, including some surprising distributions, such as the unexpected and likely common presence of Malassezia yeasts in marine habitats.

Biodiversity of Zygosaccharomyces species in food systems

Zygosaccharomyces species are among the most problematic food spoilage yeasts. The two most infamous species are Zygosaccharomyces balii and Zygosaccharomyces rouxii, although they may also take a positive role during the production of some fermented foods. DNA sequence based yeast identification aided by freely available reference databases of barcoding DNA sequences has boosted the description rate of novel yeast species in the last two decades. The genus Zygosaccharomyces has been considerably expanded as well. Especially the number of the extremely osmotolerant Zygosaccharomyces species, related to Z. rouxii and regularly found in high-sugar foods, has enlarged. A brief account of recent developments in the taxonomy and biodiversity of this important food associated genus is given in this review.

Genomic and Phenotypic Analysis of COVID-19-Associated Pulmonary Aspergillosis Isolates of Aspergillus fumigatus

The ongoing global pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for coronavirus disease 2019 (COVID-19), first described in Wuhan, China. A subset of COVID-19 patients has been reported to have acquired secondary infections by microbial pathogens, such as opportunistic fungal pathogens from the genus Aspergillus. To gain insight into COVID-19-associated pulmonary aspergillosis (CAPA), we analyzed the genomes and characterized the phenotypic profiles of four CAPA isolates of Aspergillus fumigatus obtained from patients treated in the area of North Rhine-Westphalia, Germany. By examining the mutational spectrum of single nucleotide polymorphisms, insertion-deletion polymorphisms, and copy number variants among 206 genes known to modulate A. fumigatus virulence, we found that CAPA isolate genomes do not exhibit significant differences from the genome of the Af293 reference strain. By examining a number of factors, including virulence in an invertebrate moth model, growth in the presence of osmotic, cell wall, and oxidative stressors, secondary metabolite biosynthesis, and the MIC of antifungal drugs, we found that CAPA isolates were generally, but not always, similar to A. fumigatus reference strains Af293 and CEA17. Notably, CAPA isolate D had more putative loss-of-function mutations in genes known to increase virulence when deleted. Moreover, CAPA isolate D was significantly more virulent than the other three CAPA isolates and the A. fumigatus reference strains Af293 and CEA17, but similarly virulent to two other clinical strains of A. fumigatus. These findings expand our understanding of the genomic and phenotypic characteristics of isolates that cause CAPA. IMPORTANCE The global pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the etiological agent of coronavirus disease 2019 (COVID-19), has already killed millions of people. COVID-19 patient outcome can be further complicated by secondary infections, such as COVID-19-associated pulmonary aspergillosis (CAPA). CAPA is caused by Aspergillus fungal pathogens, but there is little information about the genomic and phenotypic characteristics of CAPA isolates. We conducted genome sequencing and extensive phenotyping of four CAPA isolates of Aspergillus fumigatus from Germany. We found that CAPA isolates were often, but not always, similar to other reference strains of A. fumigatus across 206 genetic determinants of infection-relevant phenotypes, including virulence. For example, CAPA isolate D was more virulent than other CAPA isolates and reference strains in an invertebrate model of fungal disease, but similarly virulent to two other clinical strains. These results expand our understanding of COVID-19-associated pulmonary aspergillosis.

Genomic diversity, chromosomal rearrangements, and interspecies hybridization in the Ogataea polymorpha species complex

The methylotrophic yeast Ogataea polymorpha has long been a useful system for recombinant protein production, as well as a model system for methanol metabolism, peroxisome biogenesis, thermotolerance, and nitrate assimilation. It has more recently become an important model for the evolution of mating-type switching. Here, we present a population genomics analysis of 47 isolates within the O. polymorpha species complex, including representatives of the species O. polymorpha, Ogataea parapolymorpha, Ogataea haglerorum, and Ogataea angusta. We found low levels of nucleotide sequence diversity within the O. polymorpha species complex and identified chromosomal rearrangements both within and between species. In addition, we found that one isolate is an interspecies hybrid between O. polymorpha and O. parapolymorpha and present evidence for loss of heterozygosity following hybridization.

Interspecific hybridization as a driver of fungal evolution and adaptation

Cross-species gene transfer is often associated with bacteria, which have evolved several mechanisms that facilitate horizontal DNA exchange. However, the increased availability of whole-genome sequences has revealed that fungal species also exchange DNA, leading to intertwined lineages, blurred species boundaries or even novel species. In contrast to prokaryotes, fungal DNA exchange originates from interspecific hybridization, where two genomes are merged into a single, often highly unstable, polyploid genome that evolves rapidly into stabler derivatives. The resulting hybrids can display novel combinations of genetic and phenotypic variation that enhance fitness and allow colonization of new niches. Interspecific hybridization led to the emergence of important pathogens of humans and plants (for example, various Candida and 'powdery mildew' species, respectively) and industrially important yeasts, such as Saccharomyces hybrids that are important in the production of cold-fermented lagers or cold-cellared Belgian ales. In this Review, we discuss the genetic processes and evolutionary implications of fungal interspecific hybridization and highlight some of the best-studied examples. In addition, we explain how hybrids can be used to study molecular mechanisms underlying evolution, adaptation and speciation, and serve as a route towards development of new variants for industrial applications.

The Interspecific Fungal Hybrid Verticillium longisporum Displays Subgenome-Specific Gene Expression

Hybridization is an important evolutionary mechanism that can enable organisms to adapt to environmental challenges. It has previously been shown that the fungal allodiploid species Verticillium longisporum, the causal agent of verticillium stem striping in rapeseed, originated from at least three independent hybridization events between two haploid Verticillium species. To reveal the impact of genome duplication as a consequence of hybridization, we studied the genome and transcriptome dynamics upon two independent V. longisporum hybridization events, represented by the hybrid lineages “A1/D1” and “A1/D3.” We show that V. longisporum genomes are characterized by extensive chromosomal rearrangements, including between parental chromosomal sets. V. longisporum hybrids display signs of evolutionary dynamics that are typically associated with the aftermath of allodiploidization, such as haploidization and more relaxed gene evolution. The expression patterns of the two subgenomes within the two hybrid lineages are more similar than those of the shared A1 parent between the two lineages, showing that the expression patterns of the parental genomes homogenized within a lineage. However, as genes that display differential parental expression in planta do not typically display the same pattern in vitro, we conclude that subgenome-specific responses occur in both lineages. Overall, our study uncovers genomic and transcriptomic plasticity during the evolution of the filamentous fungal hybrid V. longisporum and illustrates its adaptive potential.

Network motif-based analysis of regulatory patterns in paralogous gene pairs

Current high-throughput experimental techniques make it feasible to infer gene regulatory interactions at the whole-genome level with reasonably good accuracy. Such experimentally inferred regulatory networks have become available for a number of simpler model organisms such as S. cerevisiae, and others. The availability of such networks provides an opportunity to compare gene regulatory processes at the whole genome level, and in particular, to assess similarity of regulatory interactions for homologous gene pairs either from the same or from different species. We present here a new technique for analyzing the regulatory interaction neighborhoods of paralogous gene pairs. Our central focus is the analysis of S. cerevisiae gene interaction graphs, which are of particular interest due to the ancestral whole-genome duplication (WGD) that allows to distinguish between paralogous transcription factors that are traceable to this duplication event and other paralogues. Similar analysis is also applied to E. coli and C. elegans networks. We compare paralogous gene pairs according to the presence and size of bi-fan arrays, classically associated in the literature with gene duplication, within other network motifs. We further extend this framework beyond transcription factor comparison to obtain topology-based similarity metrics based on the overlap of interaction neighborhoods applicable to most genes in a given organism. We observe that our network divergence metrics show considerably larger similarity between paralogues, especially those traceable to WGD. This is the case for both yeast and C. elegans, but not for E. coli regulatory network. While there is no obvious cross-species link between metrics, different classes of paralogues show notable differences in interaction overlap, with traceable duplications tending toward higher overlap compared to genes with shared protein families. Our findings indicate that divergence in paralogous interaction networks reflects a shared genetic origin, and that our approach may be useful for investigating structural similarity in the interaction networks of paralogous genes.

Extraordinary diversity of telomeres, telomerase RNAs and their template regions in Saccharomycetaceae

Telomerase RNA (TR) carries the template for synthesis of telomere DNA and provides a scaffold for telomerase assembly. Fungal TRs are long and have been compared to higher eukaryotes, where they show considerable diversity within phylogenetically close groups. TRs of several Saccharomycetaceae were recently identified, however, many of these remained uncharacterised in the template region. Here we show that this is mainly due to high variability in telomere sequence. We predicted the telomere sequences using Tandem Repeats Finder and then we identified corresponding putative template regions in TR candidates. Remarkably long telomere units and the corresponding putative TRs were found in Tetrapisispora species. Notably, variable lengths of the annealing sequence of the template region (1–10 nt) were found. Consequently, species with the same telomere sequence may not harbour identical TR templates. Thus, TR sequence alone can be used to predict a template region and telomere sequence, but not to determine these exactly. A conserved feature of telomere sequences, tracts of adjacent Gs, led us to test the propensity of individual telomere sequences to form G4. The results show highly diverse values of G4-propensity, indicating the lack of ubiquitous conservation of this feature across Saccharomycetaceae.

The neutral rate of whole-genome duplication varies among yeast species and their hybrids

Hybridization and polyploidization are powerful mechanisms of speciation. Hybrid speciation often coincides with whole-genome duplication (WGD) in eukaryotes. This suggests that WGD may allow hybrids to thrive by increasing fitness, restoring fertility and/or increasing access to adaptive mutations. Alternatively, it has been suggested that hybridization itself may trigger WGD. Testing these models requires quantifying the rate of WGD in hybrids without the confounding effect of natural selection. Here we show, by measuring the spontaneous rate of WGD of more than 1300 yeast crosses evolved under relaxed selection, that some genotypes or combinations of genotypes are more prone to WGD, including some hybrids between closely related species. We also find that higher WGD rate correlates with higher genomic instability and that WGD increases fertility and genetic variability. These results provide evidence that hybridization itself can promote WGD, which in turn facilitates the evolution of hybrids.

The revenge of Zygosaccharomyces yeasts in food biotechnology and applied microbiology

Non-conventional yeasts refer to a huge and still poorly explored group of species alternative to the well-known model organism Saccharomyces cerevisiae. Among them, Zygosaccharomyces rouxii and the sister species Zygosaccharomyces bailii are infamous for spoiling food and beverages even in presence of several food preservatives. On the other hand, their capability to cope with a wide range of process conditions makes these yeasts very attractive factories (the so-called "ZygoFactories") for bio-converting substrates poorly permissive for the growth of other species. In balsamic vinegar Z. rouxii is the main yeast responsible for converting highly concentrated sugars into ethanol, with a preference for fructose over glucose (a trait called fructophily). Z. rouxii has also attracted much attention for the ability to release important flavor compounds, such as fusel alcohols and the derivatives of 4-hydroxyfuranone, which markedly contribute to fragrant and smoky aroma in soy sauce. While Z. rouxii was successfully proposed in brewing for producing low ethanol beer, Z. bailii is promising for lactic acid and bioethanol production. Recently, several research efforts exploited omics tools to pinpoint the genetic bases of distinctive traits in "ZygoFactories", like fructophily, tolerance to high concentrations of sugars, lactic acid and salt. Here, I provided an overview of Zygosaccharomyces industrially relevant phenotypes and summarized the most recent findings in disclosing their genetic bases. I suggest that the increasing number of genomes available for Z. rouxii and other Zygosaccharomyces relatives, combined with recently developed genetic engineering toolkits, will boost the applications of these yeasts in biotechnology and applied microbiology.

The Role of Structural Variation in Adaptation and Evolution of Yeast and Other Fungi

Mutations in DNA can be limited to one or a few nucleotides, or encompass larger deletions, insertions, duplications, inversions and translocations that span long stretches of DNA or even full chromosomes. These so-called structural variations (SVs) can alter the gene copy number, modify open reading frames, change regulatory sequences or chromatin structure and thus result in major phenotypic changes. As some of the best-known examples of SV are linked to severe genetic disorders, this type of mutation has traditionally been regarded as negative and of little importance for adaptive evolution. However, the advent of genomic technologies uncovered the ubiquity of SVs even in healthy organisms. Moreover, experimental evolution studies suggest that SV is an important driver of evolution and adaptation to new environments. Here, we provide an overview of the causes and consequences of SV and their role in adaptation, with specific emphasis on fungi since these have proven to be excellent models to study SV.

Identification of a novel gene controlling homothallism in the yeast Kazachstania naganishii isolated in Japan

A novel gene controlling homothallic life cycle was identified in the yeast Kazachstania naganishii isolated in Japan. This gene was isolated by means of complementing a mutation, mti1, which had led to heterothallism from original homothallism in the yeast. The configuration of original mutation in MTI1 gene revealed that a truncated product is formed due to occurrence of a stop codon by a nucleotide insertion. When the gene was disrupted with a marker, the disruptant spore clone was haploid and stably heterothallic. Disfunction of the gene caused inability to self-diploidize due to defect of mating-type interconversion. The gene MTI1 (for Mating Type Interconversion) is a weak homolog of the Saccharomyces cerevisiae VID22/ENV11, which has been reported to function in vacuolar protein processing. K. naganishii has a gene representing significant homology with the HO gene of S. cerevisiae on chromosome V, which has not been clarified to be involved in regulation of life cycle in K. naganishii. The MTI1 gene defined in this study is located on K. naganishii chromosome IV and does not represent significant homology to the above ScHO-like gene and any other genes concerning life cycles of yeasts. From the viewpoint of gene evolution, it is extremely interesting that the MTI1 gene is a new type of gene controlling homothallism in addition to an HO-type gene, leading to discovery of an unknown mechanism regulating life cycles in yeasts.

Domesticating a food spoilage yeast into an organic acid-tolerant metabolic engineering host: Lactic acid production by engineered Zygosaccharomyces bailii

Lactic acid represents an important class of commodity chemicals, which can be produced by microbial cell factories. However, due to the toxicity of lactic acid at lower pH, microbial production requires the usage of neutralizing agents to maintain neutral pH. Zygosaccharomyces bailii, a food spoilage yeast, can grow under the presence of organic acids used as food preservatives. This unique trait of the yeast might be useful for producing lactic acid. With the goal of domesticating the organic acid-tolerant yeast as a metabolic engineering host, seven Z. bailii strains were screened in a minimal medium with 10 g/L of acetic, or 60 g/L of lactic acid at pH 3. The Z. bailii NRRL Y7239 strain was selected as the most robust strain to be engineered for lactic acid production. By applying a PAN-ARS-based CRISPR-Cas9 system consisting of a transfer RNA promoter and NAT selection, we demonstrated the targeted deletion of ADE2 and site-specific integration of Rhizopus oryzae ldhA coding for lactate dehydrogenase into the PDC1 locus. The resulting pdc1::ldhA strain produced 35 g/L of lactic acid without ethanol production. This study demonstrates the feasibility of the CRISPR-Cas9 system in Z. bailii, which can be applied for a fundamental study of the species.

Hybridization and the origin of new yeast lineages

Hybrids originate from the mating of two diverged organisms, resulting in novel lineages that have chimeric genomes. Hybrids may exhibit unique phenotypic traits that are not necessarily intermediate between those present in the progenitors. These unique traits may enable them to thrive in new environments. Many hybrid lineages have been discovered among yeasts in the Saccharomycotina, of which many have industrial or clinical relevance, but this might reflect a bias toward investigating species with relevance to humans. Hybridization has also been proposed to be at the root of the whole-genome duplication in the lineage leading to Saccharomyces cerevisiae. Thus, hybridization seems to have played a prominent role in the evolution of Saccharomycotina yeasts, although it is still unclear how common this evolutionary process has been during the evolution of this and other fungal clades. Similarly, the evolutionary aftermath of hybridization, including implications at the genomic, transcriptional, physiological or ecological levels, remains poorly understood. In this review, I survey recent findings from genomic analysis of yeast hybrids of industrial or clinical relevance, and discuss the evolutionary implications of genomic hybridization for the origin of new lineages, including when such hybridization results in a whole-genome duplication.

Pathogenic Allodiploid Hybrids of Aspergillus Fungi

Interspecific hybridization substantially alters genotypes and phenotypes and can give rise to new lineages. Hybrid isolates that differ from their parental species in infection-relevant traits have been observed in several human-pathogenic yeasts and plant-pathogenic filamentous fungi but have yet to be found in human-pathogenic filamentous fungi. We discovered 6 clinical isolates from patients with aspergillosis originally identified as Aspergillus nidulans (section Nidulantes) that are actually allodiploid hybrids formed by the fusion of Aspergillus spinulosporus with an unknown close relative of Aspergillus quadrilineatus, both in section Nidulantes. Evolutionary genomic analyses revealed that these isolates belong to Aspergillus latus, an allodiploid hybrid species. Characterization of diverse infection-relevant traits further showed that A. latus hybrid isolates are genomically and phenotypically heterogeneous but also differ from A. nidulans, A. spinulosporus, and A. quadrilineatus. These results suggest that allodiploid hybridization contributes to the genomic and phenotypic diversity of filamentous fungal pathogens of humans.

Diversity and Postzygotic Evolution of the Mitochondrial Genome in Hybrids of Saccharomyces Species Isolated by Double Sterility Barrier

Eukaryotic species are reproductively isolated by sterility barriers that prevent interspecies fertilization (prezygotic sterility barrier) or the fertilization results in infertile offspring (postzygotic sterility barrier). The Saccharomyces species are isolated by postzygotic sterility barriers. Their allodiploid hybrids form no viable gametes (ascospores) and the viable ascospores of the allotetraploids cannot fertilize (conjugate). Our previous work revealed that this mechanism of reproductive isolation differs from those operating in plants and animals and we designated it double sterility barrier (the failure of homeologous chromosomes to pair and the repression of mating by mating-type heterozygosity). Other studies implicated nucleo-mitochondrial incompatibilities in the sterility of the Saccharomyces hybrids, a mechanism assumed to play a central role in the reproductive isolation of animal species. In this project the mitochondrial genomes of 50 cevarum (S. cerevisiae × S. uvarum) hybrids were analyzed. 62% had S. cerevisiae mitotypes, 4% had S. uvarum mitotypes, and 34% had recombinant mitotypes. All but one hybrid formed viable spores indicating that they had genomes larger than allodiploid. Most of these spores were sterile (no sporulation in the clone of vegetative descendants; a feature characteristic of allodiploids). But regardless of their mitotypes, most hybrids could also form fertile alloaneuploid spore clones at low frequencies upon the loss of the MAT-carrying chromosome of the S. uvarum subgenome during meiosis. Hence, the cevarum alloploid nuclear genome is compatible with both parental mitochondrial genomes as well as with their recombinants, and the sterility of the hybrids is maintained by the double sterility barrier (determined in the nuclear genome) rather than by nucleo-mitochondrial incompatibilities. During allotetraploid sporulation both the nuclear and the mitochondrial genomes of the hybrids could segregate but no correlation was observed between the sterility or the fertility of the spore clones and their mitotypes. Nucleo-mitochondrial incompatibility was manifested as respiration deficiency in certain meiotic segregants. As respiration is required for meiosis-sporulation but not for fertilization (conjugation), these segregants were deficient only in sporulation. Thus, the nucleo-mitochondrial incompatibility affects the sexual processes only indirectly through the inactivation of respiration and causes only partial sterility in certain segregant spore clones.

MAT heterozygosity and the second sterility barrier in the reproductive isolation of Saccharomyces species

The genetic analysis of large numbers of Saccharomyces cerevisiae × S. uvarum ("cevarum") and S. kudriavzevii × S. uvarum ("kudvarum") hybrids in our previous studies revealed that these species are isolated by a postzygotic double-sterility barrier. We proposed a model in which the first barrier is due to the abruption of the meiotic process by the failure of the chromosomes of the subgenomes to pair (and recombine) in meiosis and the second barrier is assumed to be the result of the suppression of mating by allospecific MAT heterozygosity. While the former is analogous to the major mechanism of postzygotic reproductive isolation in plants and animals, the latter seems to be Saccharomyces specific. To bolster the assumed involvement of MAT in the second sterility barrier, we produced synthetic alloploid two-species cevarum and kudvarum hybrids with homo- and heterothallic backgrounds as well as three-species S. cerevisiae × S. kudvarum × S. uvarum ("cekudvarum") hybrids by mass-mating and examined their MAT loci using species- and cassette-specific primer pairs. We found that the allospecific MAT heterozygosity repressed MAT switching and mating in the hybrids and in the viable but sterile spores produced by the cevarum hybrids that had increased (allotetraploid) genomes. The loss of heterozygosity by meiotic malsegregation of MAT-carrying chromosomes in the latter hybrids broke down the sterility barrier. The resulting spores nullisomic for the S. uvarum chromosome produced vegetative cells capable of MAT switching and conjugation, opening the way for GARMe (Genome Autoreduction in Meiosis), the process that leads to chimeric genomes.

The yeast mating-type switching endonuclease HO is a domesticated member of an unorthodox homing genetic element family

The mating-type switching endonuclease HO plays a central role in the natural life cycle of Saccharomyces cerevisiae, but its evolutionary origin is unknown. HO is a recent addition to yeast genomes, present in only a few genera close to Saccharomyces. Here we show that HO is structurally and phylogenetically related to a family of unorthodox homing genetic elements found in Torulaspora and Lachancea yeasts. These WHO elements home into the aldolase gene FBA1, replacing its 3' end each time they integrate. They resemble inteins but they operate by a different mechanism that does not require protein splicing. We show that a WHO protein cleaves Torulaspora delbrueckii FBA1 efficiently and in an allele-specific manner, leading to DNA repair by gene conversion or NHEJ. The DNA rearrangement steps during WHO element homing are very similar to those during mating-type switching, and indicate that HO is a domesticated WHO-like element.

In the same way as a sperm from a male and an egg from a female join together to form an embryo in most animals, yeast cells have two sexes that coordinate how they reproduce. These are called “mating types” and, rather than male or female, an individual yeast cell can either be mating type “a” or “alpha”. Every yeast cell contains the genes for both mating types, and each cell’s mating type is determined by which of those genes it has active.

Only one mating type gene can be ‘on’ at a time, but some yeast species can swap mating type on demand by switching the corresponding genes ‘on’ or ‘off’. This switch is unusual. Rather than simply activate one of the genes it already has, the yeast cell keeps an inactive version of each mating type gene tucked away, makes a copy of the gene it wants to be active and pastes that copy into a different location in its genome. To do all of this yeast need another gene called HO. This gene codes for an enzyme that cuts the DNA at the location of the active mating type gene. This makes an opening that allows the cell to replace the ‘a’ gene with the ‘alpha’ gene, or vice versa. This system allows yeast cells to continue mating even if all the cells in a colony start off as the same mating type. But, cutting into the DNA is risky, and can damage the health of the cell. So, why did yeast cells evolve a system that could cause them harm?

To find out where the HO gene came from, Coughlan et al. searched through all the available genomes from yeast species for other genes with similar sequences and identified a cluster which they nicknamed “weird HO” genes, or WHO genes for short. Testing these genes revealed that they also code for enzymes that make cuts in the yeast genome, but the way the cell repairs the cuts is different. The WHO genes are jumping genes. When the enzyme encoded by a WHO gene makes a cut in the genome, the yeast cell copies the gene into the gap, allowing the gene to ‘jump’ from one part of the genome to another. It is possible that this was the starting point for the evolution of the HO gene. Changes to a WHO gene could have allowed it to cut into the mating type region of the yeast genome, giving the yeast an opportunity to ‘domesticate’ it. Over time, the yeast cell stopped the WHO gene from jumping into the gap and started using the cut to change its mating type.

Understanding how cells adapt genes for different purposes is a key question in evolutionary biology. There are many other examples of domesticated jumping genes in other organisms, including in the human immune system. Understanding the evolution of HO not only sheds light on how yeast mating type switching evolved, but on how other species might harness and adapt their genes.

Genome structure reveals the diversity of mating mechanisms in Saccharomyces cerevisiae x Saccharomyces kudriavzevii hybrids, and the genomic instability that promotes phenotypic diversity

Interspecific hybridization has played an important role in the evolution of eukaryotic organisms by favouring genetic interchange between divergent lineages to generate new phenotypic diversity involved in the adaptation to new environments. This way, hybridization between Saccharomyces species, involving the fusion between their metabolic capabilities, is a recurrent adaptive strategy in industrial environments. In the present study, whole-genome sequences of natural hybrids between Saccharomyces cerevisiae and Saccharomyces kudriavzevii were obtained to unveil the mechanisms involved in the origin and evolution of hybrids, as well as the ecological and geographic contexts in which spontaneous hybridization and hybrid persistence take place. Although Saccharomyces species can mate using different mechanisms, we concluded that rare-mating is the most commonly used, but other mechanisms were also observed in specific hybrids. The preponderance of rare-mating was confirmed by performing artificial hybridization experiments. The mechanism used to mate determines the genomic structure of the hybrid and its final evolutionary outcome. The evolution and adaptability of the hybrids are triggered by genomic instability, resulting in a wide diversity of genomic rearrangements. Some of these rearrangements could be adaptive under the stressful conditions of the industrial environment.

Investigation of Genetic Relationships Between Hanseniaspora Species Found in Grape Musts Revealed Interspecific Hybrids With Dynamic Genome Structures

Hanseniaspora, a predominant yeast genus of grape musts, includes sister species recently reported as fast evolving. The aim of this study was to investigate the genetic relationships between the four most closely related species, at the population level. A multi-locus sequence typing strategy based on five markers was applied on 107 strains, confirming the clear delineation of species H. uvarum, H. opuntiae, H. guilliermondii, and H. pseudoguilliermondii. Huge variations were observed in the level of intraspecific nucleotide diversity, and differences in heterozygosity between species indicate different life styles. No clear population structure was detected based on geographical or substrate origins. Instead, H. guilliermondii strains clustered into two distinct groups, which may reflect a recent step toward speciation. Interspecific hybrids were detected between H. opuntiae and H. pseudoguilliermondii. Their characterization using flow cytometry, karyotypes and genome sequencing showed different genome structures in different ploidy contexts: allodiploids, allotriploids, and allotetraploids. Subculturing of an allotriploid strain revealed chromosome loss equivalent to one chromosome set, followed by an auto-diploidization event, whereas another auto-diploidized tetraploid showed a segmental duplication. Altogether, these results suggest that Hanseniaspora genomes are not only fast evolving but also highly dynamic.

Origin of Lactose Fermentation in Kluyveromyces lactis by Interspecies Transfer of a Neo-functionalized Gene Cluster during Domestication

Humans have used yeasts to make cheese and kefir for millennia, but the ability to ferment the milk sugar lactose is found in only a few yeast species, of which the foremost is Kluyveromyces lactis [1]. Two genes, LAC12 (lactose permease) and LAC4 (lactase), are sufficient for lactose uptake and hydrolysis to glucose and galactose [2]. Here, we show that these genes have a complex evolutionary history in the genus Kluyveromyces that is likely the result of human activity during domestication. We show that the ancestral Lac12 was bifunctional, able to import both lactose and cellobiose into the cell. These disaccharides were then hydrolyzed by Lac4 in the case of lactose or Cel2 in the case of cellobiose. A second cellobiose transporter, Cel1, was also present ancestrally. In the K. lactis lineage, the ancestral LAC12 and LAC4 were lost and a separate upheaval in the sister species K. marxianus resulted in loss of CEL1 and quadruplication of LAC12. One of these LAC12 genes became neofunctionalized to encode an efficient lactose transporter capable of supporting fermentation, specifically in dairy strains of K. marxianus, where it formed a LAC4-LAC12-CEL2 gene cluster, although another remained a cellobiose transporter. Then, the ability to ferment lactose was acquired very recently by K. lactis var. lactis by introgression of LAC12 and LAC4 on a 15-kb subtelomeric region from a dairy strain of K. marxianus. The genomic history of the LAC genes shows that strong selective pressures were imposed on yeasts by early dairy farmers.

Genome expansion by allopolyploidization in the fungal strain Coniochaeta 2T2.1 and its exceptional lignocellulolytic machinery

BACKGROUND

Particular species of the genus Coniochaeta (Sordariomycetes) exhibit great potential for bioabatement of furanic compounds and have been identified as an underexplored source of novel lignocellulolytic enzymes, especially Coniochaeta ligniaria. However, there is a lack of information about their genomic features and metabolic capabilities. Here, we report the first in-depth genome/transcriptome survey of a Coniochaeta species (strain 2T2.1).

RESULTS

The genome of Coniochaeta sp. strain 2T2.1 has a size of 74.53 Mbp and contains 24,735 protein-encoding genes. Interestingly, we detected a genome expansion event, resulting ~ 98% of the assembly being duplicated with 91.9% average nucleotide identity between the duplicated regions. The lack of gene loss, as well as the high divergence and strong genome-wide signatures of purifying selection between copies indicates that this is likely a recent duplication, which arose through hybridization between two related Coniochaeta-like species (allopolyploidization). Phylogenomic analysis revealed that 2T2.1 is related Coniochaeta sp. PMI546 and Lecythophora sp. AK0013, which both occur endophytically. Based on carbohydrate-active enzyme (CAZy) annotation, we observed that even after in silico removal of its duplicated content, the 2T2.1 genome contains exceptional lignocellulolytic machinery. Moreover, transcriptomic data reveal the overexpression of proteins affiliated to CAZy families GH11, GH10 (endoxylanases), CE5, CE1 (xylan esterases), GH62, GH51 (α-l-arabinofuranosidases), GH12, GH7 (cellulases), and AA9 (lytic polysaccharide monoxygenases) when the fungus was grown on wheat straw compared with glucose as the sole carbon source.

CONCLUSIONS

We provide data that suggest that a recent hybridization between the genomes of related species may have given rise to Coniochaeta sp. 2T2.1. Moreover, our results reveal that the degradation of arabinoxylan, xyloglucan and cellulose are key metabolic processes in strain 2T2.1 growing on wheat straw. Different genes for key lignocellulolytic enzymes were identified, which can be starting points for production, characterization and/or supplementation of enzyme cocktails used in saccharification of agricultural residues. Our findings represent first steps that enable a better understanding of the reticulate evolution and "eco-enzymology" of lignocellulolytic Coniochaeta species.

Evolutionary and functional patterns of shared gene neighbourhood in fungi

Gene clusters comprise genomically co-localized and potentially co-regulated genes that tend to be conserved across species. In eukaryotes, multiple examples of metabolic gene clusters are known, particularly among fungi and plants. However, little is known about how gene clustering patterns vary among taxa or with respect to functional roles. Furthermore, mechanisms of the formation, maintenance and evolution of gene clusters remain unknown. We surveyed 341 fungal genomes to discover gene clusters shared by different species, independently of their functions. We inferred 12,120 cluster families, which comprised roughly one third of the gene space and were enriched in genes associated with diverse cellular functions. Additionally, most clusters did not encode transcription factors, suggesting that they are regulated distally. We used phylogenomics to characterize the evolutionary history of these clusters. We found that most clusters originated once and were transmitted vertically, coupled to differential loss. However, convergent evolution-that is, independent appearance of the same cluster-was more prevalent than anticipated. Finally, horizontal gene transfer of entire clusters was somewhat restricted, with the exception of those associated with secondary metabolism. Altogether, our results provide insights on the evolution of gene clustering as well as a broad catalogue of evolutionarily conserved gene clusters whose function remains to be elucidated.

Spontaneous whole-genome duplication restores fertility in interspecific hybrids

Interspecies hybrids often show some advantages over parents but also frequently suffer from reduced fertility, which can sometimes be overcome through sexual reproduction that sorts out genetic incompatibilities. Sex is however inefficient due to the low viability or fertility of hybrid offspring and thus limits their evolutionary potential. Mitotic cell division could be an alternative to fertility recovery in species such as fungi that can also propagate asexually. Here, to test this, we evolve in parallel and under relaxed selection more than 600 diploid yeast inter-specific hybrids that span from 100,000 to 15 M years of divergence. We find that hybrids can recover fertility spontaneously and rapidly through whole-genome duplication. These events occur in both hybrids between young and well-established species. Our results show that the instability of ploidy in hybrid is an accessible path to spontaneous fertility recovery.

Hybridization across species can lead to offspring with reduced fertility. Here, the authors experimentally evolve yeast and show that whole-genome duplication during asexual reproduction can restore fertility in hybrids over a relatively short evolutionary timespan.

Multiple Reinventions of Mating-type Switching during Budding Yeast Evolution

Cell type in budding yeasts is determined by the genotype at the mating-type (MAT) locus, but yeast species differ widely in their mating compatibility systems and life cycles. Among sexual yeasts, heterothallic species are those in which haploid strains fall into two distinct and stable mating types (MATa and MATα), whereas homothallic species are those that can switch mating types or that appear not to have distinct mating types [1, 2]. The evolutionary history of these mating compatibility systems is uncertain, particularly regarding the number and direction of transitions between homothallism and heterothallism, and regarding whether the process of mating-type switching had a single origin [3, 4, 5]. Here, we inferred the mating compatibility systems of 332 budding yeast species from their genome sequences. By reference to a robust phylogenomic tree [6], we detected evolutionary transitions between heterothallism and homothallism, and among different forms of homothallism. We find that mating-type switching has arisen independently at least 11 times during yeast evolution and that transitions from heterothallism to homothallism greatly outnumber transitions in the opposite direction (31 versus 3). Although the 3-locus MAT-HML-HMR mechanism of mating-type switching as seen in Saccharomyces cerevisiae had a single evolutionary origin in budding yeasts, simpler “flip/flop” mechanisms of switching evolved separately in at least 10 other groups of yeasts. These results point to the adaptive value of homothallism and mating-type switching to unicellular fungi.

•

Mating-type switching by flip-flopping arose at least 10 separate times

•

Mating-type switching by copy-and-paste arose only once in budding yeasts

•

Transitions toward homothallism outnumber transitions away from homothallism

•

Heterothallic species can become homothallic by DNA introgression

The extremely halotolerant black yeast Hortaea werneckii - a model for intraspecific hybridization in clonal fungi

The polymorphic black yeast Hortaea werneckii (Capnodiales, Ascomycota) is extremely halotolerant (growth from 0 to 30% [w/v] NaCl) and has been extensively studied as a model for halotolerance in Eukaryotes for over two decades. Its most frequent sources are hypersaline environments and adjacent sea-water habitats in temperate, subtropical and tropical climates. Although typically saprobic, H. werneckii can also act as a commensal coloniser on human skin, causing tinea nigra on hands and soles. Here, we report that addition of NaCl to culture media expands the growth range of H. werneckii to 37 °C, which explains its colonisation of human skin, with its increased salinity. The morphological and physiological plasticity/ versatility of H. werneckii indicate that a species complex might be involved. This was investigated in this polyphasic taxonomic analysis based on the global diversity of H. werneckii strains collected from hypersaline environments, and from humans and animals. Analysis of D1/D2domains of 28S and internal transcribed spacer rDNA revealed 10 and 17 genotypes, respectively, that were not always compliant. The genotypes have global distributions. Human and environmental strains with the same genotypes are intermingled. Due to the limited number of phylogenetically informative characters in the ribosomal DNA dataset, the partial genes encoding for β-tubulin (BTB) and mini-chromosome maintenance protein (MCM7) were also sequenced. The use of these genes was hampered by ambiguous sequences obtained by Sanger sequencing, as a consequence of the diploid and highly heterozygous genome of many H. werneckii strains. Analysis of the BTB and MCM7 genes showed that in some cases two copies of the gene from the same genome are positioned in distant phylogenetic clusters of the intraspecific gene tree. Analysis of whole-genome sequences of selected H. werneckii strains generally confirmed the phylogenetic distances estimated on the basis of ribosomal genes, but also showed substantial reticulation within the phylogenetic history of the strains. This is in line with the hypothesis that the diploid genomes of H. werneckii were formed by hybridizations, which have sometimes occurred between relatively divergent strains.

Zygosaccharomyces pseudobailii, another yeast interspecies hybrid that regained fertility by damaging one of its MAT loci

Interspecies hybridization is an important evolutionary mechanism in yeasts. The genus Zygosaccharomyces in particular contains numerous hybrid strains and/or species. Here, we investigated the genome of Zygosaccharomyces strain MT15, an isolate from Maotai-flavor Chinese liquor fermentation. We found that it is an interspecies hybrid and identified it as Zygosaccharomyces pseudobailii. The Z. bailii species complex consists of three species: Z. bailii, which is not a hybrid and whose 10 Mb genome is designated 'A', and two hybrid species Z. parabailii ('AB' genome, 20 Mb) and Z. pseudobailii ('AC' genome, 20 Mb). The A, B and C subgenomes are all approximately 7%-10% different from one another in nucleotide sequence, and are derived from three different parental species. Despite being hybrids, Z. pseudobailii and Z. parabailii are capable of mating and sporulating. We previously showed that Z. parabailii regained fertility when one copy of its MAT locus became broken into two parts, causing the allodiploid hybrid to behave as a haploid gamete. In Z. pseudobailii, we find that a very similar process occurred after hybridization, when a deletion of 1.5 kb inactivated one of the two copies of its MAT locus. The half-sibling species Z. parabailii and Z. pseudobailii therefore went through remarkably parallel but independent steps to regain fertility after they were formed by separate interspecies hybridizations.

Interplay of Chimeric Mating-Type Loci Impairs Fertility Rescue and Accounts for Intra-Strain Variability in Zygosaccharomyces rouxii Interspecies Hybrid ATCC42981

The pre-whole genome duplication (WGD) Zygosaccharomyces clade comprises several allodiploid strain/species with industrially interesting traits. The salt-tolerant yeast ATCC42981 is a sterile and allodiploid strain which contains two subgenomes, one of them resembling the haploid parental species Z. rouxii. Recently, different mating-type-like (MTL) loci repertoires were reported for ATCC42981 and the Japanese strain JCM22060, which are considered two stocks of the same strain. MTL reconstruction by direct sequencing approach is challenging due to gene redundancy, structure complexities, and allodiploid nature of ATCC42981. Here, DBG2OLC and MaSuRCA hybrid de novo assemblies of ONT and Illumina reads were combined with in vitro long PCR to definitively solve these incongruences. ATCC42981 exhibits several chimeric MTL loci resulting from reciprocal translocation between parental haplotypes and retains two MATa/MATα expression loci, in contrast to MATα in JCM22060. Consistently to these reconstructions, JCM22060, but not ATCC42981, undergoes mating and meiosis. To ascertain whether the damage of one allele at the MAT locus regains the complete sexual cycle in ATCC42981, we removed the MATα expressed locus by gene deletion. The resulting MATa/- hemizygous mutants did not show any evidence of sporulation, as well as of self- and out-crossing fertility, probably because incomplete silencing at the chimeric HMLα cassette masks the loss of heterozygosity at the MAT locus. We also found that MATα deletion switched off a2 transcription, an activator of a-specific genes in pre-WGD species. These findings suggest that regulatory scheme of cell identity needs to be further investigated in Z. rouxii protoploid yeast.

Physiological Genomics of the Highly Weak-Acid-Tolerant Food Spoilage Yeasts of Zygosaccharomyces bailii sensu lato

Zygosaccharomyces bailii and two closely related species, Z. parabailii and Z. pseudobailii ("Z. bailii species complex", "Z. bailii sensu lato" or simply "Z. bailii (s.l.)"), are frequently implicated in the spoilage of acidified preserved foods and beverages due to their tolerance to very high concentrations of weak acids used as food preservatives. The recent sequencing and annotation of these species' genomes have clarified their genomic organization and phylogenetic relationship, which includes events of interspecies hybridization. Mechanistic insights into their adaptation and tolerance to weak acids (e.g., acetic and lactic acids) are also being revealed. Moreover, the potential of Z. bailii (s.l.) to be used in industrial biotechnological processes as interesting cell factories for the production of organic acids, reduction of the ethanol content, increase of alcoholic beverages aroma complexity, as well as of genetic source for increasing weak acid resistance in yeast, is currently being considered. This chapter includes taxonomical, ecological, physiological, and biochemical aspects of Z. bailii (s.l.). The focus is on the exploitation of physiological genomics approaches that are providing the indispensable holistic knowledge to guide the effective design of strategies to overcome food spoilage or the rational exploitation of these yeasts as promising cell factories.

Tempo and Mode of Genome Evolution in the Budding Yeast Subphylum

Budding yeasts (subphylum Saccharomycotina) are found in every biome and are as genetically diverse as plants or animals. To understand budding yeast evolution, we analyzed the genomes of 332 yeast species, including 220 newly sequenced ones, which represent nearly one-third of all known budding yeast diversity. Here, we establish a robust genus-level phylogeny comprising 12 major clades, infer the timescale of diversification from the Devonian period to the present, quantify horizontal gene transfer (HGT), and reconstruct the evolution of 45 metabolic traits and the metabolic toolkit of the budding yeast common ancestor (BYCA). We infer that BYCA was metabolically complex and chronicle the tempo and mode of genomic and phenotypic evolution across the subphylum, which is characterized by very low HGT levels and widespread losses of traits and the genes that control them. More generally, our results argue that reductive evolution is a major mode of evolutionary diversification.

TERribly Difficult: Searching for Telomerase RNAs in Saccharomycetes

The telomerase RNA in yeasts is large, usually >1000 nt, and contains functional elements that have been extensively studied experimentally in several disparate species. Nevertheless, they are very difficult to detect by homology-based methods and so far have escaped annotation in the majority of the genomes of Saccharomycotina. This is a consequence of sequences that evolve rapidly at nucleotide level, are subject to large variations in size, and are highly plastic with respect to their secondary structures. Here, we report on a survey that was aimed at closing this gap in RNA annotation. Despite considerable efforts and the combination of a variety of different methods, it was only partially successful. While 27 new telomerase RNAs were identified, we had to restrict our efforts to the subgroup Saccharomycetacea because even this narrow subgroup was diverse enough to require different search models for different phylogenetic subgroups. More distant branches of the Saccharomycotina remain without annotated telomerase RNA.

Genomic evidence for intraspecific hybridization in a clonal and extremely halotolerant yeast

Background

The black yeast Hortaea werneckii (Dothideomycetes, Ascomycota) is one of the most extremely halotolerant fungi, capable of growth at NaCl concentrations close to saturation. Although dothideomycetous fungi are typically haploid, the reference H. werneckii strain has a diploid genome consisting of two subgenomes with a high level of heterozygosity.

Results

In order to explain the origin of the H. werneckii diploid genome we here report the genome sequencing of eleven strains isolated from different habitats and geographic locations. Comparison of nine diploid and two haploid strains showed that the reference genome was likely formed by hybridization between two haploids and not by endoreduplication as suggested previously. Results also support additional hybridization events in the evolutionary history of investigated strains, however exchange of genetic material in the species otherwise appears to be rare. Possible links between such unusual reproduction and the extremotolerance of H. werneckii remain to be investigated.

Conclusions

H. werneckii appears to be able to form persistent haploid as well as diploid strains, is capable of occasional hybridization between relatively heterozygous haploids, but is otherwise limited to clonal reproduction. The reported data and the first identification of haploid H. werneckii strains establish this species as a good model for studying the effects of ploidy and hybridization in an extremotolerant system unperturbed by frequent genetic recombination.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4751-5) contains supplementary material, which is available to authorized users.

Ploidy Variation in Kluyveromyces marxianus Separates Dairy and Non-dairy Isolates

Kluyveromyces marxianus is traditionally associated with fermented dairy products, but can also be isolated from diverse non-dairy environments. Because of thermotolerance, rapid growth and other traits, many different strains are being developed for food and industrial applications but there is, as yet, little understanding of the genetic diversity or population genetics of this species. K. marxianus shows a high level of phenotypic variation but the only phenotype that has been clearly linked to a genetic polymorphism is lactose utilisation, which is controlled by variation in the LAC12 gene. The genomes of several strains have been sequenced in recent years and, in this study, we sequenced a further nine strains from different origins. Analysis of the Single Nucleotide Polymorphisms (SNPs) in 14 strains was carried out to examine genome structure and genetic diversity. SNP diversity in K. marxianus is relatively high, with up to 3% DNA sequence divergence between alleles. It was found that the isolates include haploid, diploid, and triploid strains, as shown by both SNP analysis and flow cytometry. Diploids and triploids contain long genomic tracts showing loss of heterozygosity (LOH). All six isolates from dairy environments were diploid or triploid, whereas 6 out 7 isolates from non-dairy environment were haploid. This also correlated with the presence of functional LAC12 alleles only in dairy haplotypes. The diploids were hybrids between a non-dairy and a dairy haplotype, whereas triploids included three copies of a dairy haplotype.

Adaptive Response and Tolerance to Acetic Acid in Saccharomyces cerevisiae and Zygosaccharomyces bailii: A Physiological Genomics Perspective

Acetic acid is an important microbial growth inhibitor in the food industry; it is used as a preservative in foods and beverages and is produced during normal yeast metabolism in biotechnological processes. Acetic acid is also a major inhibitory compound present in lignocellulosic hydrolysates affecting the use of this promising carbon source for sustainable bioprocesses. Although the molecular mechanisms underlying Saccharomyces cerevisiae response and adaptation to acetic acid have been studied for years, only recently they have been examined in more detail in Zygosaccharomyces bailii. However, due to its remarkable tolerance to acetic acid and other weak acids this yeast species is a major threat in the spoilage of acidic foods and beverages and considered as an interesting alternative cell factory in Biotechnology. This review paper emphasizes genome-wide strategies that are providing global insights into the molecular targets, signaling pathways and mechanisms behind S. cerevisiae and Z. bailii tolerance to acetic acid, and extends this information to other weak acids whenever relevant. Such comprehensive perspective and the knowledge gathered in these two yeast species allowed the identification of candidate molecular targets, either for the design of effective strategies to overcome yeast spoilage in acidic foods and beverages, or for the rational genome engineering to construct more robust industrial strains. Examples of successful applications are provided.

Transcriptional Response to Lactic Acid Stress in the Hybrid Yeast Zygosaccharomyces parabailii

Lactic acid has a wide range of applications starting from its undissociated form, and its production using cell factories requires stress-tolerant microbial hosts. The interspecies hybrid yeast Zygosaccharomyces parabailii has great potential to be exploited as a novel host for lactic acid production, due to high organic acid tolerance at low pH and a fermentative metabolism with a high growth rate. Here we used mRNA sequencing (RNA-seq) to analyze Z. parabailii's transcriptional response to lactic acid added exogenously, and we explore the biological mechanisms involved in tolerance. Z. parabailii contains two homeologous copies of most genes. Under lactic acid stress, the two genes in each homeolog pair tend to diverge in expression to a significantly greater extent than under control conditions, indicating that stress tolerance is facilitated by interactions between the two gene sets in the hybrid. Lactic acid induces downregulation of genes related to cell wall and plasma membrane functions, possibly altering the rate of diffusion of lactic acid into cells. Genes related to iron transport and redox processes were upregulated, suggesting an important role for respiratory functions and oxidative stress defense. We found differences in the expression profiles of genes putatively regulated by Haa1 and Aft1/Aft2, previously described as lactic acid responsive in Saccharomyces cerevisiae. Furthermore, formate dehydrogenase (FDH) genes form a lactic acid-responsive gene family that has been specifically amplified in Z. parabailii in comparison to other closely related species. Our study provides a useful starting point for the engineering of Z. parabailii as a host for lactic acid production.

IMPORTANCE Hybrid yeasts are important in biotechnology because of their tolerance to harsh industrial conditions. The molecular mechanisms of tolerance can be studied by analyzing differential gene expression under conditions of interest and relating gene expression patterns to protein functions. However, hybrid organisms present a challenge to the standard use of mRNA sequencing (RNA-seq) to study transcriptional responses to stress, because their genomes contain two similar copies of almost every gene. Here we used stringent mapping methods and a high-quality genome sequence to study the transcriptional response to lactic acid stress in Zygosaccharomyces parabailii ATCC 60483, a natural interspecies hybrid yeast that contains two complete subgenomes that are approximately 7% divergent in sequence. Beyond the insights we gained into lactic acid tolerance in this study, the methods we developed will be broadly applicable to other yeast hybrid strains.

Mechanism for Restoration of Fertility in Hybrid Zygosaccharomyces rouxii Generated by Interspecies Hybridization

The mechanism of whole-genome duplication (WGD) in yeast has been intensively studied because it has a large impact on yeast evolution. WGD has shaped the genomic architecture of modern Saccharomyces cerevisiae; however, the mechanism for restoring fertility after interspecies hybridization, which would be involved in the process of WGD, has not been thoroughly elucidated. In this study, we obtained a draft genome sequence of the salt-tolerant yeast Zygosaccharomyces rouxii NBRC110957 and revealed that it is a hybrid lineage of Z. rouxii (allodiploid) with two subgenomes equivalent to NBRC1876. Because this allodiploid yeast can mate with other allodiploid strains and form spores, it can be a good model of restoring fertility after interspecies hybridization. We observed that NBRC110957 and NBRC1876 contain six mating-type-like (MTL) loci. There are no large deletions or deleterious mutations in MTL loci, except for several-base-pair deletions in the X region in certain MTL loci. We also assigned only one mating-type (MAT) locus that exclusively determines mating types from six MTL loci. These results suggest that it is possible to recover mating competence regardless of whether cells lose one MAT locus through random gene loss by mitotically dividing after interspecies hybridization. Moreover, we propose that perturbation of gene expression and substantial breakdown of MAT heterozygosity caused by chromosomal rearrangement at MTL loci play roles in restoring the mating competence of allodiploids. This scenario can provide a mechanism for restoring fertility after interspecies hybridization that is compatible with random gene loss models and suggests genomic plasticity during WGD in yeast.IMPORTANCE A whole-genome duplication occurred in an ancestor of the baker's yeast Saccharomyces cerevisiae The origins of this complex and multifaceted process, which requires intra- or interspecies hybridization followed by dysfunction of one mating-type (MAT) locus to regain mating competence, has not been thoroughly elucidated. In this study, we provide a mechanism for regaining fertility in an interspecies hybrid, Zygosaccharomyces rouxii The draft genome sequence analysis and mating test showed that the Z. rouxii strain used in this study is an intact interspecies hybrid, suggesting that it is possible to recover fertility regardless of whether cells lose one MAT locus.