Whole-Genome Sequencing and Intraspecific Analysis of the Yeast Species Lachancea quebecensis

Genomic factors shape carbon and nitrogen metabolic niche breadth across Saccharomycotina yeasts

Organisms exhibit extensive variation in ecological niche breadth, from very narrow (specialists) to very broad (generalists). Two general paradigms have been proposed to explain this variation: (i) trade-offs between performance efficiency and breadth and (ii) the joint influence of extrinsic (environmental) and intrinsic (genomic) factors. We assembled genomic, metabolic, and ecological data from nearly all known species of the ancient fungal subphylum Saccharomycotina (1154 yeast strains from 1051 species), grown in 24 different environmental conditions, to examine niche breadth evolution. We found that large differences in the breadth of carbon utilization traits between yeasts stem from intrinsic differences in genes encoding specific metabolic pathways, but we found limited evidence for trade-offs. These comprehensive data argue that intrinsic factors shape niche breadth variation in microbes.

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.

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

sppIDer: A Species Identification Tool to Investigate Hybrid Genomes with High-Throughput Sequencing

The genomics era has expanded our knowledge about the diversity of the living world, yet harnessing high-throughput sequencing data to investigate alternative evolutionary trajectories, such as hybridization, is still challenging. Here we present sppIDer, a pipeline for the characterization of interspecies hybrids and pure species, that illuminates the complete composition of genomes. sppIDer maps short-read sequencing data to a combination genome built from reference genomes of several species of interest and assesses the genomic contribution and relative ploidy of each parental species, producing a series of colorful graphical outputs ready for publication. As a proof-of-concept, we use the genus Saccharomyces to detect and visualize both interspecies hybrids and pure strains, even with missing parental reference genomes. Through simulation, we show that sppIDer is robust to variable reference genome qualities and performs well with low-coverage data. We further demonstrate the power of this approach in plants, animals, and other fungi. sppIDer is robust to many different inputs and provides visually intuitive insight into genome composition that enables the rapid identification of species and their interspecies hybrids. sppIDer exists as a Docker image, which is a reusable, reproducible, transparent, and simple-to-run package that automates the pipeline and installation of the required dependencies (https://github.com/GLBRC/sppIDer; last accessed September 6, 2018).

André Goffeau's imprinting on second generation yeast "genomologists"

All authors of the present paper have worked in labs that participated to the sequencing effort of the Saccharomyces cerevisiae reference genome, and we owe to this the fact that we have all chosen to work on genomics of yeasts. S. cerevisiae has been a popular model species for genetics since the 20th century as well as being a model for general eukaryotic cellular processes. Although it has also been used empirically in fermentation for millennia, there was until recently, a lack of knowledge about the natural and evolutionary history of this yeast. The achievement of the international effort to sequence its genome was the foundation for understanding many eukaryotic biological processes but also represented the first step towards the study of the genome and ecological diversity of yeast populations worldwide. We will describe recent advances in yeast comparative and population genomics that find their origins in the S. cerevisiae genome project initiated and pursued by André Goffeau.

Incorporating comparative genomics into the design-test-learn cycle of microbial strain engineering

Engineering microbes with new properties is an important goal in industrial engineering, to establish biological factories for production of biofuels, commodity chemicals and pharmaceutics. But engineering microbes to produce new compounds with high yield remains a major challenge toward economically viable production. Incorporating several modern approaches, including synthetic and systems biology, metabolic modeling and regulatory rewiring, has proven to significantly advance industrial strain engineering. This review highlights how comparative genomics can also facilitate strain engineering, by identifying novel genes and pathways, regulatory mechanisms and genetic background effects for engineering. We discuss how incorporating comparative genomics into the design-test-learn cycle of strain engineering can provide novel information that complements other engineering strategies.

PHENOS: a high-throughput and flexible tool for microorganism growth phenotyping on solid media

Background

Microbial arrays, with a large number of different strains on a single plate printed with robotic precision, underpin an increasing number of genetic and genomic approaches. These include Synthetic Genetic Array analysis, high-throughput Quantitative Trait Loci (QTL) analysis and 2-hybrid techniques. Measuring the growth of individual colonies within these arrays is an essential part of many of these techniques but is useful for any work with arrays. Measurement is typically done using intermittent imagery fed into complex image analysis software, which is not especially accurate and is challenging to use effectively. We have developed a simple and fast alternative technique that uses a pinning robot and a commonplace microplate reader to continuously measure the thickness of colonies growing on solid agar, complemented by a technique for normalizing the amount of cells initially printed to each spot of the array in the first place. We have developed software to automate the process of combining multiple sets of readings, subtracting agar absorbance, and visualizing colony thickness changes in a number of informative ways.

Results

The “PHENOS” pipeline (PHENotyping On Solid media), optimized for Saccharomyces yeasts, produces highly reproducible growth curves and is particularly sensitive to low-level growth. We have empirically determined a formula to estimate colony cell count from an absorbance measurement, and shown this to be comparable with estimates from measurements in liquid. We have also validated the technique by reproducing the results of an earlier QTL study done with conventional liquid phenotyping, and found PHENOS to be considerably more sensitive.

Conclusions

“PHENOS” is a cost effective and reliable high-throughput technique for quantifying growth of yeast arrays, and is likely to be equally very useful for a range of other types of microbial arrays. A detailed guide to the pipeline and software is provided with the installation files at https://github.com/gact/phenos.

Electronic supplementary material

The online version of this article (10.1186/s12866-017-1143-y) contains supplementary material, which is available to authorized users.

Genome Diversity and Evolution in the Budding Yeasts (Saccharomycotina)

Considerable progress in our understanding of yeast genomes and their evolution has been made over the last decade with the sequencing, analysis, and comparisons of numerous species, strains, or isolates of diverse origins. The role played by yeasts in natural environments as well as in artificial manufactures, combined with the importance of some species as model experimental systems sustained this effort. At the same time, their enormous evolutionary diversity (there are yeast species in every subphylum of Dikarya) sparked curiosity but necessitated further efforts to obtain appropriate reference genomes. Today, yeast genomes have been very informative about basic mechanisms of evolution, speciation, hybridization, domestication, as well as about the molecular machineries underlying them. They are also irreplaceable to investigate in detail the complex relationship between genotypes and phenotypes with both theoretical and practical implications. This review examines these questions at two distinct levels offered by the broad evolutionary range of yeasts: inside the best-studied Saccharomyces species complex, and across the entire and diversified subphylum of Saccharomycotina. While obviously revealing evolutionary histories at different scales, data converge to a remarkably coherent picture in which one can estimate the relative importance of intrinsic genome dynamics, including gene birth and loss, vs. horizontal genetic accidents in the making of populations. The facility with which novel yeast genomes can now be studied, combined with the already numerous available reference genomes, offer privileged perspectives to further examine these fundamental biological questions using yeasts both as eukaryotic models and as fungi of practical importance.

Dissection of quantitative traits by bulk segregant mapping in a protoploid yeast species

Since more than a decade ago, Saccharomyces cerevisiae has been used as a model to dissect complex traits, revealing the genetic basis of a large number of traits in fine detail. However, to have a more global view of the genetic architecture of traits across species, the examination of the molecular basis of phenotypes within non-conventional species would undoubtedly be valuable. In this respect, the Saccharomycotina yeasts represent ideal and potential non-model organisms. Here we sought to assess the feasibility of genetic mapping by bulk segregant analysis in the protoploid Lachancea kluyveri (formerly S. kluyveri) yeast species, a distantly related species to S. cerevisiae For this purpose, we designed a fluorescent mating-type marker, compatible with any mating-competent strains representative of this species, to rapidly create a large population of haploid segregants (>10(5) cells). Quantitative trait loci can be mapped by selecting and sequencing an enriched pool of progeny with extreme phenotypic values. As a test bed, we applied this strategy and mapped the causal loci underlying halotolerance phenotypes in L. kluyveri Overall, this study demonstrates that bulk segregant mapping is a powerful way for investigating the genetic basis of natural variations in non-model yeast organisms and more precisely in L. kluyveri.

Reconstruction of ancestral chromosome architecture and gene repertoire reveals principles of genome evolution in a model yeast genus

Reconstructing genome history is complex but necessary to reveal quantitative principles governing genome evolution. Such reconstruction requires recapitulating into a single evolutionary framework the evolution of genome architecture and gene repertoire. Here, we reconstructed the genome history of the genus Lachancea that appeared to cover a continuous evolutionary range from closely related to more diverged yeast species. Our approach integrated the generation of a high-quality genome data set; the development of AnChro, a new algorithm for reconstructing ancestral genome architecture; and a comprehensive analysis of gene repertoire evolution. We found that the ancestral genome of the genus Lachancea contained eight chromosomes and about 5173 protein-coding genes. Moreover, we characterized 24 horizontal gene transfers and 159 putative gene creation events that punctuated species diversification. We retraced all chromosomal rearrangements, including gene losses, gene duplications, chromosomal inversions and translocations at single gene resolution. Gene duplications outnumbered losses and balanced rearrangements with 1503, 929, and 423 events, respectively. Gene content variations between extant species are mainly driven by differential gene losses, while gene duplications remained globally constant in all lineages. Remarkably, we discovered that balanced chromosomal rearrangements could be responsible for up to 14% of all gene losses by disrupting genes at their breakpoints. Finally, we found that nonsynonymous substitutions reached fixation at a coordinated pace with chromosomal inversions, translocations, and duplications, but not deletions. Overall, we provide a granular view of genome evolution within an entire eukaryotic genus, linking gene content, chromosome rearrangements, and protein divergence into a single evolutionary framework.

Large-Scale Survey of Intraspecific Fitness and Cell Morphology Variation in a Protoploid Yeast Species

It is now clear that the exploration of the genetic and phenotypic diversity of nonmodel species greatly improves our knowledge in biology. In this context, we recently launched a population genomic analysis of the protoploid yeast Lachancea kluyveri (formerly Saccharomyces kluyveri), highlighting a broad genetic diversity (π = 17 × 10⁻³) compared to the yeast model organism, S. cerevisiae (π = 4 × 10⁻³). Here, we sought to generate a comprehensive view of the phenotypic diversity in this species. In total, 27 natural L. kluyveri isolates were subjected to trait profiling using the following independent approaches: (i) analyzing growth in 55 growth conditions and (ii) investigating 501 morphological changes at the cellular level. Despite higher genetic diversity, the fitness variance observed in L. kluyveri is lower than that in S. cerevisiae. However, morphological features show an opposite trend. In addition, there is no correlation between the origins (ecological or geographical) of the isolate and the phenotypic patterns, demonstrating that trait variation follows neither population history nor source environment in L. kluyveri. Finally, pairwise comparisons between growth rate correlation and genetic diversity show a clear decrease in phenotypic variability linked to genome variation increase, whereas no such a trend was identified for morphological changes. Overall, this study reveals for the first time the phenotypic diversity of a distantly related species to S. cerevisiae. Given its genetic properties, L. kluyveri might be useful in further linkage mapping analyses of complex traits, and could ultimately provide a better insight into the evolution of the genotype–phenotype relationship across yeast species.