Fluctuating selection facilitates the discovery of broadly effective but difficult to reach adaptive outcomes in yeast.

Evolutionary compromises are thought to be common under fluctuating selection because the mutations that best enable adaptation to one environmental context can often be detrimental to others. Yet, prior experimental work has shown that generalists can sometimes perform as well as specialists in their own environments. Here we use a highly replicated evolutionary experiment (N = 448 asexual lineages of the brewer's yeast) to show that even though fluctuation between two environmental conditions often induces evolutionary compromises (at least early on), it can also help reveal difficult to reach adaptive outcomes that ultimately improve performance in both environments. Specifically, we begin by showing that yeast adaptation to chemical stress can involve fitness trade-offs with stress-free environments and that, accordingly, lineages that are repeatedly exposed to occasional stress tend to respond by trading performance for breadth of adaptation. We then show that on rare occasions, fluctuating selection leads to the evolution of no-cost generalists that can even outcompete constant selection specialists in their own environments. We propose that the discovery of these broader and more effective adaptive outcomes under fluctuating selection could be partially facilitated by changes in the adaptive landscape that result from having to deal with fitness trade-offs across different environmental conditions. Overall, our findings indicate that reconciling the short- and long-term evolutionary consequences of fluctuating selection could significantly improve our understanding of the evolution of specialization and generalism.

A connection between the ribosome and two S. pombe tRNA modification mutants subject to rapid tRNA decay.

tRNA modifications are crucial in all organisms to ensure tRNA folding and stability, and accurate translation. In both the yeast Saccharomyces cerevisiae and the evolutionarily distant yeast Schizosaccharomyces pombe, mutants lacking certain tRNA body modifications (outside the anticodon loop) are temperature sensitive due to rapid tRNA decay (RTD) of a subset of hypomodified tRNAs. Here we show that for each of two S. pombe mutants subject to RTD, mutations in ribosomal protein genes suppress the temperature sensitivity without altering tRNA levels. Prior work showed that S. pombe trm8Δ mutants, lacking 7-methylguanosine, were temperature sensitive due to RTD, and that one class of suppressors had mutations in the general amino acid control (GAAC) pathway, which was activated concomitant with RTD, resulting in further tRNA loss. We now find that another class of S. pombe trm8Δ suppressors have mutations in rpl genes, encoding 60S subunit proteins, and that suppression occurs with minimal restoration of tRNA levels and reduced GAAC activation. Furthermore, trm8Δ suppression extends to other mutations in the large or small ribosomal subunit. We also find that S. pombe tan1Δ mutants, lacking 4-acetylcytidine, are temperature sensitive due to RTD, that one class of suppressors have rpl mutations, associated with minimal restoration of tRNA levels, and that suppression extends to other rpl and rps mutations. However, although S. pombe tan1Δ temperature sensitivity is associated with some GAAC activation, suppression by an rpl mutation only modestly inhibits GAAC activation. We propose a model in which ribosomal protein mutations result in reduced ribosome concentrations, leading to both reduced ribosome collisions and a reduced requirement for tRNA, with these effects having different relative importance in trm8Δ and tan1Δ mutants. This model is consistent with our results in S. cerevisiae trm8Δ trm4Δ mutants, known to undergo RTD, fueling speculation that this model applies across eukaryotes.

Divergence in the Saccharomyces Species' Heat Shock Response Is Indicative of Their Thermal Tolerance.

The Saccharomyces species have diverged in their thermal growth profile. Both Saccharomyces cerevisiae and Saccharomyces paradoxus grow at temperatures well above the maximum growth temperature of Saccharomyces kudriavzevii and Saccharomyces uvarum but grow more poorly at lower temperatures. In response to thermal shifts, organisms activate a stress response that includes heat shock proteins involved in protein homeostasis and acquisition of thermal tolerance. To determine whether Saccharomyces species have diverged in their response to temperature, we measured changes in gene expression in response to a 12 °C increase or decrease in temperature for four Saccharomyces species and their six pairwise hybrids. To ensure coverage of subtelomeric gene families, we sequenced, assembled, and annotated a complete S. uvarum genome. In response to heat, the cryophilic species showed a stronger stress response than the thermophilic species, and the hybrids showed a mixture of parental responses that depended on the time point. After an initial strong response indicative of high thermal stress, hybrids with a thermophilic parent resolved their heat shock response to become similar to their thermophilic parent. Within the hybrids, only a small number of temperature-responsive genes showed consistent differences between alleles from the thermophilic and cryophilic species. Our results show that divergence in the heat shock response is mainly a consequence of a strain's thermal tolerance, suggesting that cellular factors that signal heat stress or resolve heat-induced changes are relevant to thermal divergence in the Saccharomyces species.

A connection between the ribosome and two S. pombe tRNA modification mutants subject to rapid tRNA decay.

tRNA modifications are crucial in all organisms to ensure tRNA folding and stability, and accurate translation in the ribosome. In both the yeast Saccharomyces cerevisiae and the evolutionarily distant yeast Schizosaccharomyces pombe, mutants lacking certain tRNA body modifications (outside the anticodon loop) are temperature sensitive due to rapid tRNA decay (RTD) of a subset of hypomodified tRNAs. Here we show that for each of two S. pombe mutants subject to RTD, mutations in ribosomal protein genes suppress the temperature sensitivity without altering tRNA levels. Prior work showed that S. pombe trm8Δ mutants, lacking 7-methylguanosine, were temperature sensitive due to RTD and that one class of suppressors had mutations in the general amino acid control (GAAC) pathway, which was activated concomitant with RTD, resulting in further tRNA loss. We now find that another class of S. pombe trm8Δ suppressors have mutations in rpl genes, encoding 60S subunit proteins, and that suppression occurs with minimal restoration of tRNA levels and reduced GAAC activation. Furthermore, trm8Δ suppression extends to other mutations in the large or small ribosomal subunit. We also find that S. pombe tan1Δ mutants, lacking 4-acetylcytidine, are temperature sensitive due to RTD, that one class of suppressors have rpl mutations, associated with minimal restoration of tRNA levels, and that suppression extends to other rpl and rps mutations. However, although S. pombe tan1Δ temperature sensitivity is associated with some GAAC activation, suppression by an rpl mutation does not significantly inhibit GAAC activation. These results suggest that ribosomal protein mutations suppress the temperature sensitivity of S. pombe trm8Δ and tan1Δ mutants due to reduced ribosome concentrations, leading to both a reduced requirement for tRNA, and reduced ribosome collisions and GAAC activation. Results with S. cerevisiae trm8Δ trm4Δ mutants are consistent with this model, and fuel speculation that similar results will apply across eukaryotes.

Divergence in the Saccharomyces species' heat shock response is indicative of their thermal tolerance.

The Saccharomyces species have diverged in their thermal growth profile. Both S. cerevisiae and S. paradoxus grow at temperatures well above the maximum growth temperature of S. kudriavzevii and S. uvarum, but grow more poorly at lower temperatures. In response to thermal shifts, organisms activate a stress response that includes heat shock proteins involved in protein homeostasis and acquisition of thermal tolerance. To determine whether Saccharomyces species have diverged in their response to temperature we measured changes in gene expression in response to a 12°C increase or decrease in temperature for four Saccharomyces species and their six pairwise hybrids. To ensure coverage of subtelomeric gene families we sequenced, assembled and annotated a complete S. uvarum genome. All the strains exhibited a stronger response to heat than cold treatment. In response to heat, the cryophilic species showed a stronger response than the thermophilic species. The hybrids showed a mixture of parental stress responses depending on the time point. After the initial response, hybrids with a thermophilic parent were more similar to S. cerevisiae and S. paradoxus, and the S. cerevisiae × S. paradoxus hybrid showed the weakest heat shock response. Within the hybrids a small subset of temperature responsive genes showed species specific responses but most were also hybrid specific. Our results show that divergence in the heat shock response is indicative of a strain's thermal tolerance, suggesting that cellular factors that signal heat stress or resolve heat induced changes are relevant to thermal divergence in the Saccharomyces species.

Experimental evolution of Saccharomyces uvarum at high temperature yields elevation of maximal growth temperature and loss of the mitochondrial genome.

An organism's upper thermal tolerance is a major driver of its ecology and is a complex, polygenic trait. Given the wide variance in this critical phenotype across the tree of life, it is quite striking that this trait has not proven very evolutionarily labile in experimental evolution studies of microbes. In stark contrast to recent studies, William Henry Dallinger in the 1880s reported increasing the upper thermal limit of microbes he experimentally evolved by >40°C using a very gradual temperature ramping strategy. Using a selection scheme inspired by Dallinger, we sought to increase the upper thermal limit of Saccharomyces uvarum . This species has a maximum growth temperature of 34-35°C, considerably lower than S. cerevisiae . After 136 passages on solid plates at progressively higher temperatures, we recovered a clone that can grow at 36°C, a gain of ~1.5°C. Additionally, the evolved clone lost its mitochondrial genome and cannot respire. In contrast, an induced rho 0 derivative of the ancestor shows a decrease in thermotolerance. Also, incubation of the ancestor at 34°C for 5 days increased the frequency of petite mutants drastically compared to 22°C, supporting the notion that mutation pressure rather than selection drove loss of mtDNA in the evolved clone. These results demonstrate that S. uvarum 's upper thermal limit can be elevated slightly via experimental evolution and corroborate past observations in S. cerevisiae that high temperature selection schemes can inadvertently lead to production of the potentially undesirable respiratory incompetent phenotype in yeasts.

Macroevolutionary diversity of traits and genomes in the model yeast genus Saccharomyces.

Species is the fundamental unit to quantify biodiversity. In recent years, the model yeast Saccharomyces cerevisiae has seen an increased number of studies related to its geographical distribution, population structure, and phenotypic diversity. However, seven additional species from the same genus have been less thoroughly studied, which has limited our understanding of the macroevolutionary events leading to the diversification of this genus over the last 20 million years. Here, we show the geographies, hosts, substrates, and phylogenetic relationships for approximately 1,800 Saccharomyces strains, covering the complete genus with unprecedented breadth and depth. We generated and analyzed complete genome sequences of 163 strains and phenotyped 128 phylogenetically diverse strains. This dataset provides insights about genetic and phenotypic diversity within and between species and populations, quantifies reticulation and incomplete lineage sorting, and demonstrates how gene flow and selection have affected traits, such as galactose metabolism. These findings elevate the genus Saccharomyces as a model to understand biodiversity and evolution in microbial eukaryotes.

Off the deep end: What can deep learning do for the gene expression field?

After a COVID-related hiatus, the fifth biennial symposium on Evolution and Core Processes in Gene Regulation met at the Stowers Institute in Kansas City, Missouri July 21 to 24, 2022. This symposium, sponsored by the American Society for Biochemistry and Molecular Biology (ASBMB), featured experts in gene regulation and evolutionary biology. Topic areas covered enhancer evolution, the cis-regulatory code, and regulatory variation, with an overall focus on bringing the power of deep learning (DL) to decipher DNA sequence information. DL is a machine learning method that uses neural networks to learn complex rules that make predictions about diverse types of data. When DL models are trained to predict genomic data from DNA sequence information, their high prediction accuracy allows the identification of impactful genetic variants within and across species. In addition, the learned sequence rules can be extracted from the model and provide important clues about the mechanistic underpinnings of the cis-regulatory code.

Using colony size to measure fitness in Saccharomyces cerevisiae.

Competitive fitness assays in liquid culture have been a mainstay for characterizing experimental evolution of microbial populations. Growth of microbial strains has also been extensively characterized by colony size and could serve as a useful alternative if translated to per generation measurements of relative fitness. To examine fitness based on colony size, we established a relationship between cell number and colony size for strains of Saccharomyces cerevisiae robotically pinned onto solid agar plates in a high-density format. This was used to measure growth rates and estimate relative fitness differences between evolved strains and their ancestors. After controlling for edge effects through both normalization and agar-trimming, we found that colony size is a sensitive measure of fitness, capable of detecting 1% differences. While fitnesses determined from liquid and solid mediums were not equivalent, our results demonstrate that colony size provides a sensitive means of measuring fitness that is particularly well suited to measurements across many environments.

The evolutionary patterns of barley pericentromeric chromosome regions, as shaped by linkage disequilibrium and domestication.

The distribution of recombination events along large cereal chromosomes is uneven and is generally restricted to gene-rich telomeric ends. To understand how the lack of recombination affects diversity in the large pericentromeric regions, we analysed deep exome capture data from a final panel of 815 Hordeum vulgare (barley) cultivars, landraces and wild barleys, sampled from across their eco-geographical ranges. We defined and compared variant data across the pericentromeric and non-pericentromeric regions, observing a clear partitioning of diversity both within and between chromosomes and germplasm groups. Dramatically reduced diversity was found in the pericentromeres of both cultivars and landraces when compared with wild barley. We observed a mixture of completely and partially differentiated single-nucleotide polymorphisms (SNPs) between domesticated and wild gene pools, suggesting that domesticated gene pools were derived from multiple wild ancestors. Patterns of genome-wide linkage disequilibrium, haplotype block size and number, and variant frequency within blocks showed clear contrasts among individual chromosomes and between cultivars and wild barleys. Although most cultivar chromosomes shared a single major pericentromeric haplotype, chromosome 7H clearly differentiated the two-row and six-row types associated with different geographical origins. Within the pericentromeric regions we identified 22 387 non-synonymous SNPs, 92 of which were fixed for alternative alleles in cultivar versus wild accessions. Surprisingly, only 29 SNPs found exclusively in the cultivars were predicted to be 'highly deleterious'. Overall, our data reveal an unconventional pericentromeric genetic landscape among distinct barley gene pools, with different evolutionary processes driving domestication and diversification.

Sleep behavior and chronotype before and throughout pregnancy.

STUDY OBJECTIVE: To compare sleep behavior before and during pregnancy. METHODS: In this prospective cohort study, healthy women were followed from pre-pregnancy until delivery. At pre-pregnancy and each trimester, participants completed validated questionnaires of chronotype and sleep quality and timing, including the Munich ChronoType Questionnaire, Epworth Sleepiness Scale, and Pittsburgh Sleep Quality Index. The primary outcomes were sleep period start and end times, sleep duration, sleep midpoint, and social jetlag, compared between pre-pregnancy and each trimester. Wrist actigraphy was used to measure the same outcomes in a subset of participants. RESULTS: Eighty-six women were included in analysis of questionnaires. Of these, 37 provided complete actigraphy data. Questionnaire and actigraphy data indicate that participants had less social jetlag during pregnancy than before pregnancy. Sleep period start times were earlier on both work and free days in the first and second trimesters than pre-pregnancy, and returned to pre-pregnancy times by the third trimester. Actigraphy data revealed that, compared to pre-pregnancy, participants had longer sleep periods in all trimesters on work days and in the first trimester on free days. Sleep surveys revealed that participants had poorer sleep quality in the first and third trimesters and more sleepiness in the first trimester than pre-pregnancy. CONCLUSION: The first trimester of pregnancy is characterized by earlier sleep period start time, longer sleep duration, and poorer sleep quality than pre-pregnancy. Sleep quality temporarily improves in the second trimester, and sleep period start time returns to pre-pregnancy time by the third trimester. STUDY RATIONALE: Multiple parameters of sleep have been studied in the context of pregnancy and pregnancy outcomes, but rarely in comparison to pre-pregnancy or longitudinally through pregnancy. STUDY IMPACT: Actigraphy and questionnaire data reveal sleep timing and quality change throughout pregnancy. These data on sleep changes in healthy pregnancy can be used as a baseline to identify sleep-related risk factors throughout pregnancy.

The structure and diversity of strain-level variation in vaginal bacteria.

The vaginal microbiome plays an important role in human health and species of vaginal bacteria have been associated with reproductive disease. Strain-level variation is also thought to be important, but the diversity, structure and evolutionary history of vaginal strains is not as well characterized. We developed and validated an approach to measure strain variation from metagenomic data based on SNPs within the core genomes for six species of vaginal bacteria: Gardnerella vaginalis, Lactobacillus crispatus, Lactobacillus iners, Lactobacillus jensenii, Lactobacillus gasseri and Atopobium vaginae. Despite inhabiting the same environment, strain diversity and structure varies across species. All species except L. iners are characterized by multiple distinct groups of strains. Even so, strain diversity is lower in the Lactobacillus species, consistent with a more recent colonization of the human vaginal microbiome. Both strain diversity and the frequency of multi-strain samples is related to species-level diversity of the microbiome in which they occur, suggesting similar ecological factors influencing diversity within the vaginal niche. We conclude that the structure of strain-level variation provides both the motivation and means of testing whether strain-level differences contribute to the function and health consequences of the vaginal microbiome.

High-throughput analysis of adaptation using barcoded strains of Saccharomyces cerevisiae.

BACKGROUND: Experimental evolution of microbes can be used to empirically address a wide range of questions about evolution and is increasingly employed to study complex phenomena ranging from genetic evolution to evolutionary rescue. Regardless of experimental aims, fitness assays are a central component of this type of research, and low-throughput often limits the scope and complexity of experimental evolution studies. We created an experimental evolution system in Saccharomyces cerevisiae that utilizes genetic barcoding to overcome this challenge. RESULTS: We first confirm that barcode insertions do not alter fitness and that barcode sequencing can be used to efficiently detect fitness differences via pooled competition-based fitness assays. Next, we examine the effects of ploidy, chemical stress, and population bottleneck size on the evolutionary dynamics and fitness gains (adaptation) in a total of 76 experimentally evolving, asexual populations by conducting 1,216 fitness assays and analyzing 532 longitudinal-evolutionary samples collected from the evolving populations. In our analysis of these data we describe the strengths of this experimental evolution system and explore sources of error in our measurements of fitness and evolutionary dynamics. CONCLUSIONS: Our experimental treatments generated distinct fitness effects and evolutionary dynamics, respectively quantified via multiplexed fitness assays and barcode lineage tracking. These findings demonstrate the utility of this new resource for designing and improving high-throughput studies of experimental evolution. The approach described here provides a framework for future studies employing experimental designs that require high-throughput multiplexed fitness measurements.

Associations between the vaginal microbiome and Candida colonization in women of reproductive age.

BACKGROUND: The composition of bacteria within the vaginal microbiome has garnered a lot of recent attention and has been associated with reproductive health and disease. Despite the common occurrence of yeast (primarily Candida) within the vaginal microbiome, there is still an incomplete picture of relationships between yeast and bacteria (especially lactobacilli), as well as how such associations are governed. Such relationships could be important to a more holistic understanding of the vaginal microbiome and its connection to reproductive health. OBJECTIVE: The objective of the study was to perform molecular characterization of clinical specimens to define associations between vaginal bacteria (especially Lactobacillus species) and Candida colonization. In vitro studies were conducted to test the 2 most common dominant Lactobacillus species (Lactobacillus crispatus and Lactobacillus iners) in their ability to inhibit Candida growth and to examine the basis for such inhibition. STUDY DESIGN: A nested cross-sectional study of reproductive-age women from the Contraceptive CHOICE Project was conducted. Vaginal swabs from 299 women were selected to balance race and bacterial vaginosis status, resulting in a similar representation of black and white women in each of the 3 Nugent score categories (normal [0-3], intermediate [4-6], and bacterial vaginosis [7-10]). Sequencing of the 16S ribosomal gene (V4 region) was used to determine the dominant Lactobacillus species present (primarily Lactobacillus iners and Lactobacillus crispatus), defined as >50% of the community. Subjects without dominance by a single Lactobacillus species were classified as Diverse. A Candida-specific quantitative polymerase chain reaction targeting the internally transcribed spacer 1 was validated using vaginal samples collected from a second cohort of women and used to assess Candida colonization. Two hundred fifty-five nonpregnant women with sufficient bacterial biomass for analysis were included in the final analysis. Generalized linear models were used to evaluate associations between Lactobacillus dominance, sociodemographic and risk characteristics, and vaginal Candida colonization. In separate in vitro studies, the potential of cell-free supernatants from Lactobacillus crispatus and Lactobacillus iners cultures to inhibit Candida growth was evaluated. RESULTS: Forty-two women (16%) were vaginally colonized with Candida. Microbiomes characterized as Diverse (38%), Lactobacillus iners-dominant (39%), and Lactobacillus crispatus-dominant (20%) were the most common. The microbiome, race, and Candida colonization co-varied with a higher prevalence of Candida among black women and Lactobacillus iners-dominant communities compared with white women and Lactobacillus crispatus-dominant communities. Lactobacillus iners-dominant communities were more likely to harbor Candida than Lactobacillus crispatus-dominant communities (odds ratio, 2.85, 95% confidence interval, 1.03-7.21; Fisher exact test, P = .048). In vitro, Lactobacillus crispatus produced greater concentrations of lactic acid and exhibited significantly more pH-dependent growth inhibition of Candida albicans, suggesting a potential mechanism for the clinical observations. CONCLUSION: In nonpregnant women, Lactobacillus iners-dominant communities were significantly more likely to harbor Candida than Lactobacillus crispatus-dominant communities, suggesting that Lactobacillus species have different relationships with Candida. In vitro experiments indicate that Lactobacillus crispatus may impede Candida colonization more effectively than Lactobacillus iners through a greater production of lactic acid.

Riding the Rhythm of Melatonin Through Pregnancy to Deliver on Time.

Pregnancy is influenced by the circadian ("circa" or approximately; diem or day) system, which coordinates physiology and behavior with predictable daily changes in the environment such as light/dark cycles. For example, most species deliver around a particular time of day. In mammals, circadian rhythms are controlled by the master circadian pacemaker, the suprachiasmatic nucleus. One key way that the suprachiasmatic nucleus coordinates circadian rhythms throughout the body is by regulating production of the sleep-promoting hormone melatonin. Serum melatonin concentration, which peaks at night and is suppressed during the day, is one of the best biological indicators of circadian timing. Circadian misalignment causes maternal disturbances in the temporal organization of many physiological processes including melatonin synthesis, and these disturbances of the circadian system have been linked to an increased risk for pregnancy complications. Here, we review evidence that melatonin helps regulate the maternal and fetal circadian systems and the timing of birth. Finally, we discuss the potential for melatonin-based therapeutic strategies to alleviate poor pregnancy outcomes such as preeclampsia and preterm birth.

The Fate of Deleterious Variants in a Barley Genomic Prediction Population.

Targeted identification and purging of deleterious genetic variants has been proposed as a novel approach to animal and plant breeding. This strategy is motivated, in part, by the observation that demographic events and strong selection associated with cultivated species pose a "cost of domestication." This includes an increase in the proportion of genetic variants that are likely to reduce fitness. Recent advances in DNA resequencing and sequence constraint-based approaches to predict the functional impact of a mutation permit the identification of putatively deleterious SNPs (dSNPs) on a genome-wide scale. Using exome capture resequencing of 21 barley lines, we identified 3855 dSNPs among 497,754 total SNPs. We generated whole-genome resequencing data of Hordeum murinum ssp. glaucum as a phylogenetic outgroup to polarize SNPs as ancestral vs. derived. We also observed a higher proportion of dSNPs per synonymous SNPs (sSNPs) in low-recombination regions of the genome. Using 5215 progeny from a genomic prediction experiment, we examined the fate of dSNPs over three breeding cycles. Adjusting for initial frequency, derived alleles at dSNPs reduced in frequency or were lost more often than other classes of SNPs. The highest-yielding lines in the experiment, as chosen by standard genomic prediction approaches, carried fewer homozygous dSNPs than randomly sampled lines from the same progeny cycle. In the final cycle of the experiment, progeny selected by genomic prediction had a mean of 5.6% fewer homozygous dSNPs relative to randomly chosen progeny from the same cycle.

Multiple Changes Underlie Allelic Divergence of CUP2 Between Saccharomyces Species.

Under the model of micromutationism, phenotypic divergence between species is caused by accumulation of many small-effect changes. While mapping the causal changes to single nucleotide resolution could be difficult for diverged species, genetic dissection via chimeric constructs allows us to evaluate whether a large-effect gene is composed of many small-effect nucleotide changes. In a previously described non-complementation screen, we found an allele difference of CUP2, a copper-binding transcription factor, underlies divergence in copper resistance between Saccharomyces cerevisiae and S. uvarum Here, we tested whether the allele effect of CUP2 was caused by multiple nucleotide changes. By analyzing chimeric constructs containing four separate regions in the CUP2 gene, including its distal promoter, proximal promoter, DNA binding domain and transcriptional activation domain, we found that all four regions of the S. cerevisiae allele conferred copper resistance, with the proximal promoter showing the largest effect, and that both additive and epistatic effects are likely involved. These findings support a model of multiple changes underlying evolution and suggest an important role of both protein coding and cis-regulatory changes in evolution.

A polyploid admixed origin of beer yeasts derived from European and Asian wine populations.

Strains of Saccharomyces cerevisiae used to make beer, bread, and wine are genetically and phenotypically distinct from wild populations associated with trees. The origins of these domesticated populations are not always clear; human-associated migration and admixture with wild populations have had a strong impact on S. cerevisiae population structure. We examined the population genetic history of beer strains and found that ale strains and the S. cerevisiae portion of allotetraploid lager strains were derived from admixture between populations closely related to European grape wine strains and Asian rice wine strains. Similar to both lager and baking strains, ale strains are polyploid, providing them with a passive means of remaining isolated from other populations and providing us with a living relic of their ancestral hybridization. To reconstruct their polyploid origin, we phased the genomes of two ale strains and found ale haplotypes to both be recombinants between European and Asian alleles and to also contain novel alleles derived from extinct or as yet uncharacterized populations. We conclude that modern beer strains are the product of a historical melting pot of fermentation technology.

Mitochondria-encoded genes contribute to evolution of heat and cold tolerance in yeast.

Genetic analysis of phenotypic differences between species is typically limited to interfertile species. Here, we conducted a genome-wide noncomplementation screen to identify genes that contribute to a major difference in thermal growth profile between two reproductively isolated yeast species, Saccharomyces cerevisiae and Saccharomyces uvarum. The screen identified only a single nuclear-encoded gene with a moderate effect on heat tolerance, but, in contrast, revealed a large effect of mitochondrial DNA (mitotype) on both heat and cold tolerance. Recombinant mitotypes indicate that multiple genes contribute to thermal divergence, and we show that protein divergence in COX1 affects both heat and cold tolerance. Our results point to the yeast mitochondrial genome as an evolutionary hotspot for thermal divergence.

Mitochondrial DNA and temperature tolerance in lager yeasts.

A growing body of research suggests that the mitochondrial genome (mtDNA) is important for temperature adaptation. In the yeast genus Saccharomyces, species have diverged in temperature tolerance, driving their use in high- or low-temperature fermentations. Here, we experimentally test the role of mtDNA in temperature tolerance in synthetic and industrial hybrids (Saccharomyces cerevisiae × Saccharomyces eubayanus or Saccharomyces pastorianus), which cold-brew lager beer. We find that the relative temperature tolerances of hybrids correspond to the parent donating mtDNA, allowing us to modulate lager strain temperature preferences. The strong influence of mitotype on the temperature tolerance of otherwise identical hybrid strains provides support for the mitochondrial climactic adaptation hypothesis in yeasts and demonstrates how mitotype has influenced the world's most commonly fermented beverage.