Genetic architecture of ethanol-responsive transcriptome variation in Saccharomyces cerevisiae strains

The fitness trade-off between growth and stress resistance determines the phenotypic landscape

BACKGROUND

A central challenge in biology is to discover a principle that determines individual phenotypic differences within a species. The growth rate is particularly important for a unicellular organism, and the growth rate under a certain condition is negatively associated with that of another condition, termed fitness trade-off. Therefore, there should exist a common molecular mechanism that regulates multiple growth rates under various conditions, but most studies so far have focused on discovering those genes associated with growth rates under a specific condition.

RESULTS

In this study, we found that there exists a recurrent gene expression signature whose expression levels are related to the fitness trade-off between growth preference and stress resistance across various yeast strains and multiple conditions. We further found that the genomic variation of stress-response, ribosomal, and cell cycle regulators are potential causal genes that determine the sensitivity between growth and survival. Intriguingly, we further observed that the same principle holds for human cells using anticancer drug sensitivities across multiple cancer cell lines.

CONCLUSIONS

Together, we suggest that the fitness trade-off is an evolutionary trait that determines individual growth phenotype within a species. By using this trait, we can possibly overcome anticancer drug resistance in cancer cells.

Contribution of the yeast bi-chaperone system in the restoration of the RNA helicase Ded1 and translational activity under severe ethanol stress

Preexposure to mild stress often improves cellular tolerance to subsequent severe stress. Severe ethanol stress (10% v/v) causes persistent and pronounced translation repression in Saccharomyces cerevisiae. However, it remains unclear whether preexposure to mild stress can mitigate translation repression in yeast cells under severe ethanol stress. We found that the translational activity of yeast cells pretreated with 6% (v/v) ethanol was initially significantly repressed under subsequent 10% ethanol but was then gradually restored even under severe ethanol stress. We also found that 10% ethanol caused the aggregation of Ded1, which plays a key role in translation initiation as a DEAD-box RNA helicase. Pretreatment with 6% ethanol led to the gradual disaggregation of Ded1 under subsequent 10% ethanol treatment in wild-type cells but not in fes1Δhsp104Δ cells, which are deficient in Hsp104 with significantly reduced capacity for Hsp70. Hsp104 and Hsp70 are key components of the bi-chaperone system that play a role in yeast protein quality control. fes1Δhsp104Δ cells did not restore translational activity under 10% ethanol, even after pretreatment with 6% ethanol. These results indicate that the regeneration of Ded1 through the bi-chaperone system leads to the gradual restoration of translational activity under continuous severe stress. This study provides new insights into the acquired tolerance of yeast cells to severe ethanol stress and the resilience of their translational activity.

Integrative Analysis of the Ethanol Tolerance of Saccharomyces cerevisiae

Ethanol (EtOH) alters many cellular processes in yeast. An integrated view of different EtOH-tolerant phenotypes and their long noncoding RNAs (lncRNAs) is not yet available. Here, large-scale data integration showed the core EtOH-responsive pathways, lncRNAs, and triggers of higher (HT) and lower (LT) EtOH-tolerant phenotypes. LncRNAs act in a strain-specific manner in the EtOH stress response. Network and omics analyses revealed that cells prepare for stress relief by favoring activation of life-essential systems. Therefore, longevity, peroxisomal, energy, lipid, and RNA/protein metabolisms are the core processes that drive EtOH tolerance. By integrating omics, network analysis, and several other experiments, we showed how the HT and LT phenotypes may arise: (1) the divergence occurs after cell signaling reaches the longevity and peroxisomal pathways, with CTA1 and ROS playing key roles; (2) signals reaching essential ribosomal and RNA pathways via SUI2 enhance the divergence; (3) specific lipid metabolism pathways also act on phenotype-specific profiles; (4) HTs take greater advantage of degradation and membraneless structures to cope with EtOH stress; and (5) our EtOH stress-buffering model suggests that diauxic shift drives EtOH buffering through an energy burst, mainly in HTs. Finally, critical genes, pathways, and the first models including lncRNAs to describe nuances of EtOH tolerance are reported here.

Adaptability of wine yeast to ethanol-induced protein denaturation

This year marks the 200th anniversary of the birth of Dr Louis Pasteur (1822-1895), who revealed that alcoholic fermentation is performed by yeast cells. Subsequently, details of the mechanisms of alcoholic fermentation and glycolysis in yeast cells have been elucidated. However, the mechanisms underlying the high tolerance and adaptability of yeast cells to ethanol are not yet fully understood. This review presents the response and adaptability of yeast cells to ethanol-induced protein denaturation. Herein, we describe the adverse effects of severe ethanol stress on intracellular proteins and the responses of yeast cells. Furthermore, recent findings on the acquired resistance of wine yeast cells to severe ethanol stress that causes protein denaturation are discussed, not only under laboratory conditions, but also during the fermentation process at 15°C to mimic the vinification process of white wine.

Functional mining of novel terpene synthases from metagenomes

BACKGROUND

Terpenes are one of the most diverse and abundant classes of natural biomolecules, collectively enabling a variety of therapeutic, energy, and cosmetic applications. Recent genomics investigations have predicted a large untapped reservoir of bacterial terpene synthases residing in the genomes of uncultivated organisms living in the soil, indicating a vast array of putative terpenoids waiting to be discovered.

RESULTS

We aimed to develop a high-throughput functional metagenomic screening system for identifying novel terpene synthases from bacterial metagenomes by relieving the toxicity of terpene biosynthesis precursors to the Escherichia coli host. The precursor toxicity was achieved using an inducible operon encoding the prenyl pyrophosphate synthetic pathway and supplementation of the mevalonate precursor. Host strain and screening procedures were finely optimized to minimize false positives arising from spontaneous mutations, which avoid the precursor toxicity. Our functional metagenomic screening of human fecal metagenomes yielded a novel β-farnesene synthase, which does not show amino acid sequence similarity to known β-farnesene synthases. Engineered S. cerevisiae expressing the screened β-farnesene synthase produced 120 mg/L β-farnesene from glucose (2.86 mg/g glucose) with a productivity of 0.721 g/L∙h.

CONCLUSIONS

A unique functional metagenomic screening procedure was established for screening terpene synthases from metagenomic libraries. This research proves the potential of functional metagenomics as a sequence-independent avenue for isolating targeted enzymes from uncultivated organisms in various environmental habitats.

Complex genetics cause and constrain fungal persistence in different parts of the mammalian body

Determining how genetic polymorphisms enable certain fungi to persist in mammalian hosts can improve understanding of opportunistic fungal pathogenesis, a source of substantial human morbidity and mortality. We examined the genetic basis of fungal persistence in mice using a cross between a clinical isolate and the lab reference strain of the budding yeast Saccharomyces cerevisiae. Employing chromosomally encoded DNA barcodes, we tracked the relative abundances of 822 genotyped, haploid segregants in multiple organs over time and performed linkage mapping of their persistence in hosts. Detected loci showed a mix of general and antagonistically pleiotropic effects across organs. General loci showed similar effects across all organs, while antagonistically pleiotropic loci showed contrasting effects in the brain vs the kidneys, liver, and spleen. Persistence in an organ required both generally beneficial alleles and organ-appropriate pleiotropic alleles. This genetic architecture resulted in many segregants persisting in the brain or in nonbrain organs, but few segregants persisting in all organs. These results show complex combinations of genetic polymorphisms collectively cause and constrain fungal persistence in different parts of the mammalian body.

Integrated multi-omics reveals common properties underlying stress granule and P-body formation

Non-membrane-bound compartments such as P-bodies (PBs) and stress granules (SGs) play important roles in the regulation of gene expression following environmental stresses. We have systematically and quantitatively determined the protein and mRNA composition of PBs and SGs formed before and after nutrient stress. We find that high molecular weight (HMW) complexes exist prior to glucose depletion that we propose may act as seeds for further condensation of proteins forming mature PBs and SGs. We identify an enrichment of proteins with low complexity and RNA binding domains, as well as long, structured mRNAs that are poorly translated following nutrient stress. Many proteins and mRNAs are shared between PBs and SGs including several multivalent RNA binding proteins that promote condensate interactions during liquid-liquid phase separation. We uncover numerous common protein and RNA components across PBs and SGs that support a complex interaction profile during the maturation of these biological condensates. These interaction networks represent a tuneable response to stress, highlighting previously unrecognized condensate heterogeneity. These studies therefore provide an integrated and quantitative understanding of the dynamic nature of key biological condensates.

Trans-acting genetic variation affects the expression of adjacent genes

Gene expression differences among individuals are shaped by trans-acting expression quantitative trait loci (eQTLs). Most trans-eQTLs map to hotspot locations that influence many genes. The molecular mechanisms perturbed by hotspots are often assumed to involve "vertical" cascades of effects in pathways that can ultimately affect the expression of thousands of genes. Here, we report that trans-eQTLs can affect the expression of adjacent genes via "horizontal" mechanisms that extend along a chromosome. Genes affected by trans-eQTL hotspots in the yeast Saccharomyces cerevisiae were more likely to be located next to each other than expected by chance. These paired hotspot effects tended to occur at adjacent genes that also show coexpression in response to genetic and environmental perturbations, suggesting shared mechanisms. Physical proximity and shared chromatin state, in addition to regulation of adjacent genes by similar transcription factors, were independently associated with paired hotspot effects among adjacent genes. Paired effects of trans-eQTLs can occur at neighboring genes even when these genes do not share a common function. This phenomenon could result in unexpected connections between regulatory genetic variation and phenotypes.

Mutational sources of trans-regulatory variation affecting gene expression in Saccharomyces cerevisiae

Heritable variation in a gene’s expression arises from mutations impacting cis- and trans-acting components of its regulatory network. Here, we investigate how trans-regulatory mutations are distributed within the genome and within a gene regulatory network by identifying and characterizing 69 mutations with trans-regulatory effects on expression of the same focal gene in Saccharomyces cerevisiae. Relative to 1766 mutations without effects on expression of this focal gene, we found that these trans-regulatory mutations were enriched in coding sequences of transcription factors previously predicted to regulate expression of the focal gene. However, over 90% of the trans-regulatory mutations identified mapped to other types of genes involved in diverse biological processes including chromatin state, metabolism, and signal transduction. These data show how genetic changes in diverse types of genes can impact a gene’s expression in trans, revealing properties of trans-regulatory mutations that provide the raw material for trans-regulatory variation segregating within natural populations.

Acquired Resistance to Severe Ethanol Stress in Saccharomyces cerevisiae Protein Quality Control

Acute severe ethanol stress (10% [vol/vol]) damages proteins and causes the intracellular accumulation of insoluble proteins in Saccharomyces cerevisiae On the other hand, a pretreatment with mild stress increases tolerance to subsequent severe stress, which is called acquired stress resistance. It currently remains unclear whether the accumulation of insoluble proteins under severe ethanol stress may be mitigated by increasing protein quality control (PQC) activity in cells pretreated with mild stress. In the present study, we examined the induction of resistance to severe ethanol stress in PQC and confirmed that a pretreatment with 6% (vol/vol) ethanol or mild thermal stress at 37°C significantly reduced insoluble protein levels and the aggregation of Lsg1, which is prone to denaturation and aggregation by stress, in yeast cells under 10% (vol/vol) ethanol stress. The induction of this stress resistance required the new synthesis of proteins; the expression of proteins comprising the bichaperone system (Hsp104, Ssa3, and Fes1), Sis1, and Hsp42 was upregulated during the pretreatment and maintained under subsequent severe ethanol stress. Since the pretreated cells of deficient mutants in the bichaperone system (fes1Δ hsp104Δ and ssa2Δ ssa3Δ ssa4Δ) failed to sufficiently reduce insoluble protein levels and Lsg1 aggregation, the enhanced activity of the bichaperone system appears to be important for the induction of adequate stress resistance. In contrast, the importance of proteasomes and aggregases (Btn2 and Hsp42) in the induction of stress resistance has not been confirmed. These results provide further insights into the PQC activity of yeast cells under severe ethanol stress, including the brewing process.IMPORTANCE Although the budding yeast S. cerevisiae, which is used in the production of alcoholic beverages and bioethanol, is highly tolerant of ethanol, high concentrations of ethanol are also stressful to the yeast and cause various adverse effects, including protein denaturation. A pretreatment with mild stress improves the ethanol tolerance of yeast cells; however, it currently remains unclear whether it increases PQC activity and reduces the levels of denatured proteins. In the present study, we found that a pretreatment with mild ethanol upregulated the expression of proteins involved in PQC and mitigated the accumulation of insoluble proteins, even under severe ethanol stress. These results provide novel insights into ethanol tolerance and the adaptive capacity of yeast. They may also contribute to research on the physiology of yeast cells during the brewing process, in which the concentration of ethanol gradually increases.

Physical, genetic and functional interactions between the eisosome protein Pil1 and the MBOAT O-acyltransferase Gup1

The Saccharomyces cerevisiae MBOAT O-acyltransferase Gup1 is involved in many processes, including cell wall and membrane composition and integrity, and acetic acid-induced cell death. Gup1 was previously shown to interact physically with the mitochondrial membrane VDAC (Voltage-Dependent Anion Channel) protein Por1 and the ammonium transceptor Mep2. By co-immunoprecipitation, the eisosome core component Pil1 was identified as a novel physical interaction partner of Gup1. The expression of PIL1 and Pil1 protein levels were found to be unaffected by GUP1 deletion. In ∆gup1 cells, Pil1 was distributed in dots (likely representing eisosomes) in the membrane, identically to wt cells. However, ∆gup1 cells presented 50% less Pil1-GFP dots/eisosomes, suggesting that Gup1 is important for eisosome formation. The two proteins also interact genetically in the maintenance of cell wall integrity, and during arsenite and acetic acid exposure. We show that Δgup1 Δpil1 cells take up more arsenite than wt and are extremely sensitive to arsenite and to acetic acid treatments. The latter causes a severe apoptotic wt-like cell death phenotype, epistatically reverting the ∆gup1 necrotic type of death. Gup1 and Pil1 are thus physically, genetically and functionally connected.

Simultaneous quantification of mRNA and protein in single cells reveals post-transcriptional effects of genetic variation

Trans-acting DNA variants may specifically affect mRNA or protein levels of genes located throughout the genome. However, prior work compared trans-acting loci mapped in separate studies, many of which had limited statistical power. Here, we developed a CRISPR-based system for simultaneous quantification of mRNA and protein of a given gene via dual fluorescent reporters in single, live cells of the yeast Saccharomyces cerevisiae. In large populations of recombinant cells from a cross between two genetically divergent strains, we mapped 86 trans-acting loci affecting the expression of ten genes. Less than 20% of these loci had concordant effects on mRNA and protein of the same gene. Most loci influenced protein but not mRNA of a given gene. One locus harbored a premature stop variant in the YAK1 kinase gene that had specific effects on protein or mRNA of dozens of genes. These results demonstrate complex, post-transcriptional genetic effects on gene expression.

Systematic analysis of nuclear gene function in respiratory growth and expression of the mitochondrial genome in S. cerevisiae

The production of metabolic energy in form of ATP by oxidative phosphorylation depends on the coordinated action of hundreds of nuclear-encoded mitochondrial proteins and a handful of proteins encoded by the mitochondrial genome (mtDNA). We used the yeast Saccharomyces cerevisiae as a model system to systematically identify the genes contributing to this process. Integration of genome-wide high-throughput growth assays with previously published large data sets allowed us to define with high confidence a set of 254 nuclear genes that are indispensable for respiratory growth. Next, we induced loss of mtDNA in the yeast deletion collection by growth on ethidium bromide-containing medium and identified twelve genes that are essential for viability in the absence of mtDNA (i.e. petite-negative). Replenishment of mtDNA by cytoduction showed that respiratory-deficient phenotypes are highly variable in many yeast mutants. Using a mitochondrial genome carrying a selectable marker, ARG8^m, we screened for mutants that are specifically defective in maintenance of mtDNA and mitochondrial protein synthesis. We found that up to 176 nuclear genes are required for expression of mitochondria-encoded proteins during fermentative growth. Taken together, our data provide a comprehensive picture of the molecular processes that are required for respiratory metabolism in a simple eukaryotic cell.

A tps1Δ persister-like state in Saccharomyces cerevisiae is regulated by MKT1

Trehalose metabolism in yeast has been linked to a variety of phenotypes, including heat resistance, desiccation tolerance, carbon-source utilization, and sporulation. The relationships among the several phenotypes of mutants unable to synthesize trehalose are not understood, even though the pathway is highly conserved. One of these phenotypes is that tps1Δ strains cannot reportedly grow on media containing glucose or fructose, even when another carbon source they can use (e.g. galactose) is present. Here we corroborate the recent observation that a small fraction of yeast tps1Δ cells do grow on glucose, unlike the majority of the population. This is not due to a genetic alteration, but instead resembles the persister phenotype documented in many microorganisms and cancer cells undergoing lethal stress. We extend these observations to show that this phenomenon is glucose-specific, as it does not occur on another highly fermented carbon source, fructose. We further demonstrate that this phenomenon appears to be related to mitochondrial complex III function, but unrelated to inorganic phosphate levels in the cell, as had previously been suggested. Finally, we found that this phenomenon is specific to S288C-derived strains, and is the consequence of a variant in the MKT1 gene.

Dynamic genetic architecture of yeast response to environmental perturbation shed light on origin of cryptic genetic variation

Cryptic genetic variation could arise from, for example, Gene-by-Gene (G-by-G) or Gene-by-Environment (G-by-E) interactions. The underlying molecular mechanisms and how they influence allelic effects and the genetic variance of complex traits is largely unclear. Here, we empirically explored the role of environmentally influenced epistasis on the suppression and release of cryptic variation by reanalysing a dataset of 4,390 haploid yeast segregants phenotyped on 20 different media. The focus was on 130 epistatic loci, each contributing to segregant growth in at least one environment and that together explained most (69–100%) of the narrow sense heritability of growth in the individual environments. We revealed that the epistatic growth network reorganised upon environmental changes to alter the estimated marginal (additive) effects of the individual loci, how multi-locus interactions contributed to individual segregant growth and the level of expressed genetic variance in growth. The estimated additive effects varied most across environments for loci that were highly interactive network hubs in some environments but had few or no interactors in other environments, resulting in changes in total genetic variance across environments. This environmentally dependent epistasis was thus an important mechanism for the suppression and release of cryptic variation in this population. Our findings increase the understanding of the complex genetic mechanisms leading to cryptic variation in populations, providing a basis for future studies on the genetic maintenance of trait robustness and development of genetic models for studying and predicting selection responses for quantitative traits in breeding and evolution.

Many biological traits are polygenic, with complex interplay between underlying genes and the surrounding environment. As a result, individuals with the same allele might have distinctive phenotypes due to differences in the polygenic background and/or the environment. Such differences often create additional genetic variation that is highly relevant to quantitative and evolutionary genetics by limiting our ability to accurately predict the phenotypes in medical or agricultural applications and providing opportunities for long term evolution. Previously, yeast growth regulating genes were found to be organised in large interacting networks. Here, we found that these networks were reorganised upon environmental changes, and that this resulted in altered effect sizes of individual genes, and how the whole network contributed to growth and the level of total genetic variance, providing a basis for future studies on the genetic maintenance of trait robustness and development of genetic models for studying and predicting selection responses for quantitative traits.

Mkt1 is required for RNAi-mediated silencing and establishment of heterochromatin in fission yeast

Constitutive domains of repressive heterochromatin are maintained within the fission yeast genome through self-reinforcing mechanisms involving histone methylation and small RNAs. Non-coding RNAs generated from heterochromatic regions are processed into small RNAs by the RNA interference pathway, and are subject to silencing through both transcriptional and post-transcriptional mechanisms. While the pathways involved in maintenance of the repressive heterochromatin state are reasonably well understood, less is known about the requirements for its establishment. Here, we describe a novel role for the post-transcriptional regulatory factor Mkt1 in establishment of heterochromatin at pericentromeres in fission yeast. Loss of Mkt1 does not affect maintenance of existing heterochromatin, but does affect its recovery following depletion, as well as de novo establishment of heterochromatin on a mini-chromosome. Pathway dissection revealed that Mkt1 is required for RNAi-mediated post-transcriptional silencing, downstream of small RNA production. Mkt1 physically associates with pericentromeric transcripts, and is additionally required for maintenance of silencing and heterochromatin at centromeres when transcriptional silencing is impaired. Our findings provide new insight into the mechanism of RNAi-mediated post-transcriptional silencing in fission yeast, and unveil an important role for post-transcriptional silencing in establishment of heterochromatin that is dispensable when full transcriptional silencing is imposed.

Transcriptional and epigenetic control of regulated cell death in yeast

Unicellular organisms like yeast can undergo controlled demise in a manner that is partly reminiscent of mammalian cell death. This is true at the levels of both mechanistic and functional conservation. Yeast offers the combination of unparalleled genetic amenability and a comparatively simple biology to understand both the regulation and evolution of cell death. In this minireview, we address the capacity of the nucleus as a regulatory hub during yeast regulated cell death (RCD), which is becoming an increasingly central question in yeast RCD research. In particular, we explore and critically discuss the available data on stressors and signals that specifically impinge on the nucleus. Moreover, we also analyze the current knowledge on nuclear factors as well as on transcriptional control and epigenetic events that orchestrate yeast RCD. Altogether we conclude that the functional significance of the nucleus for yeast RCD in undisputable, but that further exploration beyond correlative work is necessary to disentangle the role of nuclear events in the regulatory network.

DNA variants affecting the expression of numerous genes in trans have diverse mechanisms of action and evolutionary histories

DNA variants that alter gene expression contribute to variation in many phenotypic traits. In particular, trans-acting variants, which are often located on different chromosomes from the genes they affect, are an important source of heritable gene expression variation. However, our knowledge about the identity and mechanism of causal trans-acting variants remains limited. Here, we developed a fine-mapping strategy called CRISPR-Swap and dissected three expression quantitative trait locus (eQTL) hotspots known to alter the expression of numerous genes in trans in the yeast Saccharomyces cerevisiae. Causal variants were identified by engineering recombinant alleles and quantifying the effects of these alleles on the expression of a green fluorescent protein-tagged gene affected by the given locus in trans. We validated the effect of each variant on the expression of multiple genes by RNA-sequencing. The three variants differed in their molecular mechanism, the type of genes they reside in, and their distribution in natural populations. While a missense leucine-to-serine variant at position 63 in the transcription factor Oaf1 (L63S) was almost exclusively present in the reference laboratory strain, the two other variants were frequent among S. cerevisiae isolates. A causal missense variant in the glucose receptor Rgt2 (V539I) occurred at a poorly conserved amino acid residue and its effect was strongly dependent on the concentration of glucose in the culture medium. A noncoding variant in the conserved fatty acid regulated (FAR) element of the OLE1 promoter influenced the expression of the fatty acid desaturase Ole1 in cis and, by modulating the level of this essential enzyme, other genes in trans. The OAF1 and OLE1 variants showed a non-additive genetic interaction, and affected cellular lipid metabolism. These results demonstrate that the molecular basis of trans-regulatory variation is diverse, highlighting the challenges in predicting which natural genetic variants affect gene expression.

Differences in the DNA sequence of individual genomes contribute to differences in many traits, such as appearance, physiology, and the risk for common diseases. An important group of these DNA variants influences how individual genes across the genome are turned on or off. In this paper, we describe a strategy for identifying such “trans-acting” variants in different strains of baker’s yeast. We used this strategy to reveal three single DNA base changes that each influences the expression of dozens of genes. These three DNA variants were very different from each other. Two of them changed the protein sequence, one in a transcription factor and the other in a sugar sensor. The third changed the expression of an enzyme, a change that in turn caused other genes to alter their expression. One variant existed in only a few yeast isolates, while the other two existed in many isolates collected from around the world. This diversity of DNA variants that influence the expression of many other genes illustrates how difficult it is to predict which DNA variants in an individual’s genome will have effects on the organism.

Membrane Fluidity of Saccharomyces cerevisiae from Huangjiu (Chinese Rice Wine) Is Variably Regulated by OLE1 To Offset the Disruptive Effect of Ethanol

An evolution and resequencing strategy was used to research the genetic basis of Saccharomyces cerevisiae BR20 (with 18 vol% ethanol tolerance) and the evolved strain F23 (with 25 vol% ethanol tolerance). Whole-genome sequencing and RNA sequencing (RNA-seq) indicated that the enhanced ethanol tolerance under 10 vol% ethanol could be attributed to amino acid metabolism, whereas 18 vol% ethanol tolerance was due to fatty acid metabolism. Ultrastructural analysis indicated that F23 exhibited better membrane integrity than did BR20 under ethanol stress. At low concentrations (<5 vol%), the partition of ethanol into the membrane increased the membrane fluidity, which had little effect on cell growth. However, the toxic effects of medium and high ethanol concentrations (5 to 20 vol%) tended to decrease the membrane fluidity. Under high ethanol stress (>10 vol%), the highly tolerant strain was able to maintain a relatively constant fluidity by increasing the content of unsaturated fatty acid (UFA), whereas less-tolerant strains show a continuous decrease in fluidity and UFA content. OLE1, which was identified as the only gene with a differential single-nucleotide polymorphism (SNP) mutation site related to fatty acid metabolism, was significantly changed in response to ethanol. The role of OLE1 in membrane fluidity was positively validated in its overexpressed transformants. Therefore, OLE1 lowered the rate of decline in membrane fluidity and thus enabled the yeast to better fight the deleterious effects of ethanol.IMPORTANCE Yeasts with superior ethanol tolerance are desirable for winemakers and wine industries. In our previous work, strain F23 was evolved with superior ethanol tolerance and fermentation activity to improve the flavor profiles of Chinese rice wine. Therefore, exploring the genomic variations and ethanol tolerance mechanism of strain F23 could contribute to an understanding of its effect on the flavor characteristics in the resulting Chinese rice wine. The cellular membrane plays a vital role in the ethanol tolerance of yeasts; however, how the membrane is regulated to fight the toxic effect of ethanol remains to be elucidated. This study suggests that the membrane fluidity is variably regulated by OLE1 to offset the disruptive effect of ethanol. Current work will help develop more ethanol-tolerant yeast strains for wine industries and contribute to a deep understanding of its high flavor-producing ability.

Multimodal Microorganism Development: Integrating Top-Down Biological Engineering with Bottom-Up Rational Design

Biological engineering has unprecedented potential to solve society's most pressing challenges. Engineering approaches must consider complex technical, economic, and social factors. This requires methods that confer gene/pathway-level functionality and organism-level robustness in rapid and cost-effective ways. This article compares foundational engineering approaches - bottom-up, gene-targeted engineering, and top-down, whole-genome engineering - and identifies significant complementarity between them. Cases drawn from engineering Saccharomyces cerevisiae exemplify the synergy of a combined approach. Indeed, multimodal engineering streamlines strain development by leveraging the complementarity of whole-genome and gene-targeted engineering to overcome the gap in design knowledge that restricts rational design. As biological engineers target more complex systems, this dual-track approach is poised to become an increasingly important tool to realize the promise of synthetic biology.

Mechanisms of Yeast Adaptation to Wine Fermentations

Cells face genetic and/or environmental changes in order to outlast and proliferate. Characterization of changes after stress at different "omics" levels is crucial to understand the adaptation of yeast to changing conditions. Wine fermentation is a stressful situation which yeast cells have to cope with. Genome-wide analyses extend our cellular physiology knowledge by pointing out the mechanisms that contribute to sense the stress caused by these perturbations (temperature, ethanol, sulfites, nitrogen, etc.) and related signaling pathways. The model organism, Saccharomyces cerevisiae, was studied in response to industrial stresses and changes at different cellular levels (transcriptomic, proteomic, and metabolomics), which were followed statically and/or dynamically in the short and long terms. This chapter focuses on the response of yeast cells to the diverse stress situations that occur during wine fermentations, which induce perturbations, including nutritional changes, ethanol stress, temperature stress, oxidative stress, etc.

Independent Mechanisms for Acquired Salt Tolerance versus Growth Resumption Induced by Mild Ethanol Pretreatment in Saccharomyces cerevisiae

Microbes in nature frequently experience “boom or bust” cycles of environmental stress. Thus, microbes that can anticipate the onset of stress would have an advantage. One way that microbes anticipate future stress is through acquired stress resistance, where cells exposed to a mild dose of one stress gain the ability to survive an otherwise lethal dose of a subsequent stress. In the budding yeast Saccharomyces cerevisiae, certain stressors can cross protect against high salt concentrations, though the mechanisms governing this acquired stress resistance are not well understood. In this study, we took advantage of wild yeast strains to understand the mechanism underlying ethanol-induced cross protection against high salt concentrations. We found that mild ethanol stress allows cells to resume growth on high salt, which involves a novel role for a well-studied salt transporter. Overall, this discovery highlights how leveraging natural variation can provide new insights into well-studied stress defense mechanisms.

All living organisms must recognize and respond to various environmental stresses throughout their lifetime. In natural environments, cells frequently encounter fluctuating concentrations of different stressors that can occur in combination or sequentially. Thus, the ability to anticipate an impending stress is likely ecologically relevant. One possible mechanism for anticipating future stress is acquired stress resistance, where cells preexposed to a mild sublethal dose of stress gain the ability to survive an otherwise lethal dose of stress. We have been leveraging wild strains of Saccharomyces cerevisiae to investigate natural variation in the yeast ethanol stress response and its role in acquired stress resistance. Here, we report that a wild vineyard isolate possesses ethanol-induced cross protection against severe concentrations of salt. Because this phenotype correlates with ethanol-dependent induction of the ENA genes, which encode sodium efflux pumps already associated with salt resistance, we hypothesized that variation in ENA expression was responsible for differences in acquired salt tolerance across strains. Surprisingly, we found that the ENA genes were completely dispensable for ethanol-induced survival of high salt concentrations in the wild vineyard strain. Instead, the ENA genes were necessary for the ability to resume growth on high concentrations of salt following a mild ethanol pretreatment. Surprisingly, this growth acclimation phenotype was also shared by the lab yeast strain despite lack of ENA induction under this condition. This study underscores that cross protection can affect both viability and growth through distinct mechanisms, both of which likely confer fitness effects that are ecologically relevant.

IMPORTANCE Microbes in nature frequently experience “boom or bust” cycles of environmental stress. Thus, microbes that can anticipate the onset of stress would have an advantage. One way that microbes anticipate future stress is through acquired stress resistance, where cells exposed to a mild dose of one stress gain the ability to survive an otherwise lethal dose of a subsequent stress. In the budding yeast Saccharomyces cerevisiae, certain stressors can cross protect against high salt concentrations, though the mechanisms governing this acquired stress resistance are not well understood. In this study, we took advantage of wild yeast strains to understand the mechanism underlying ethanol-induced cross protection against high salt concentrations. We found that mild ethanol stress allows cells to resume growth on high salt, which involves a novel role for a well-studied salt transporter. Overall, this discovery highlights how leveraging natural variation can provide new insights into well-studied stress defense mechanisms.

Genetics of trans-regulatory variation in gene expression

Heritable variation in gene expression forms a crucial bridge between genomic variation and the biology of many traits. However, most expression quantitative trait loci (eQTLs) remain unidentified. We mapped eQTLs by transcriptome sequencing in 1012 yeast segregants. The resulting eQTLs accounted for over 70% of the heritability of mRNA levels, allowing comprehensive dissection of regulatory variation. Most genes had multiple eQTLs. Most expression variation arose from trans-acting eQTLs distant from their target genes. Nearly all trans-eQTLs clustered at 102 hotspot locations, some of which influenced the expression of thousands of genes. Fine-mapped hotspot regions were enriched for transcription factor genes. While most genes had a local eQTL, most of these had no detectable effects on the expression of other genes in trans. Hundreds of non-additive genetic interactions accounted for small fractions of expression variation. These results reveal the complexity of genetic influences on transcriptome variation in unprecedented depth and detail.

Linkage mapping of yeast cross protection connects gene expression variation to a higher-order organismal trait

Gene expression variation is extensive in nature, and is hypothesized to play a major role in shaping phenotypic diversity. However, connecting differences in gene expression across individuals to higher-order organismal traits is not trivial. In many cases, gene expression variation may be evolutionarily neutral, and in other cases expression variation may only affect phenotype under specific conditions. To understand connections between gene expression variation and stress defense phenotypes, we have been leveraging extensive natural variation in the gene expression response to acute ethanol in laboratory and wild Saccharomyces cerevisiae strains. Previous work found that the genetic architecture underlying these expression differences included dozens of “hotspot” loci that affected many transcripts in trans. In the present study, we provide new evidence that one of these expression QTL hotspot loci affects natural variation in one particular stress defense phenotype—ethanol-induced cross protection against severe doses of H₂O₂. A major causative polymorphism is in the heme-activated transcription factor Hap1p, which we show directly impacts cross protection, but not the basal H₂O₂ resistance of unstressed cells. This provides further support that distinct cellular mechanisms underlie basal and acquired stress resistance. We also show that Hap1p-dependent cross protection relies on novel regulation of cytosolic catalase T (Ctt1p) during ethanol stress in a wild oak strain. Because ethanol accumulation precedes aerobic respiration and accompanying reactive oxygen species formation, wild strains with the ability to anticipate impending oxidative stress would likely be at an advantage. This study highlights how strategically chosen traits that better correlate with gene expression changes can improve our power to identify novel connections between gene expression variation and higher-order organismal phenotypes.

A major goal in genetics is to understand how individuals with different genetic makeups respond to their environment. Understanding these “gene-environment interactions” is important for the development of personalized medicine. For example, gene-environment interactions can explain why some people are more sensitive to certain drugs or are more likely to get certain cancers. While the underlying causes of gene-environment interactions are unclear, one possibility is that differences in gene expression across individuals are responsible. In this study, we examined that possibility using baker’s yeast as a model. We were interested in a phenomenon called acquired stress resistance, where cells exposed to a mild dose of one stress can become resistant to an otherwise lethal dose of severe stress. This response is observed in diverse organisms ranging from bacteria to humans, though the specific mechanisms governing acquisition of higher stress resistance are poorly understood. To understand the differences between yeast strains with and without the ability to acquire further stress resistance, we employed genetic mapping. We found that part of the variation in acquired stress resistance was due to sequence differences in a key regulatory protein, thus providing new insight into how different individuals respond to acute environmental change.

Transcriptome profiling of Issatchenkia orientalis under ethanol stress

Issatchenkia orientalis, a non-Saccharomyces yeast that can resist a wide variety of environmental stresses, has potential use in winemaking and bioethanol production. Little is known about gene expression or the physiology of I. orientalis under ethanol stress. In this study, high-throughput RNA sequencing was used to investigate the transcriptome profile of I. orientalis in response to ethanol. 502 gene transcripts were differentially expressed, of which 451 were more abundant, and 51 less abundant, in cells subjected to 4 h of ethanol stress (10% v/v). Annotation and statistical analyses suggest that multiple genes involved in ergosterol biosynthesis, trehalose metabolism, and stress response are differentially expressed under these conditions. The up-regulation of molecular chaperones HSP90 and HSP70, and also genes associated with the ubiquitin–proteasome proteolytic pathway suggests that ethanol stress may cause aggregation of misfolded proteins. Finally, ethanol stress in I. orientalis appears to have a nitrogen starvation effect, and many genes involved in nutrient uptake were up-regulated.

Membrane fluidification by ethanol stress activates unfolded protein response in yeasts

The toxic effect of ethanol is one of the most important handicaps for many biotechnological applications of yeasts, such as bioethanol production. Elucidation of ethanol stress response will help to improve yeast performance in biotechnological processes. In the yeast Saccharomyces cerevisiae, ethanol stress has been recently described as an activator of the unfolded protein response (UPR), a conserved intracellular signalling pathway that regulates the transcription of ER homoeostasis‐related genes. However, the signal and activation mechanism has not yet been unravelled. Here, we studied UPR's activation after ethanol stress and observed the upregulation of the key target genes, like INO1, involved in lipid metabolism. We found that inositol content influenced UPR activation after ethanol stress and we observed significant changes in lipid composition, which correlate with a major membrane fluidity alteration by this amphipathic molecule. Then, we explored the hypothesis that membrane fluidity changes cause UPR activation upon ethanol stress by studying UPR response against fluidification or rigidification agents and by studying a mutant, erg2, with altered membrane fluidity. The results suggest that the membrane fluidification effects of ethanol and other agents are the signal for UPR activation, a mechanism that has been proposed in higher eukaryotes.

Gene-gene and gene-environment interactions in complex traits in yeast

One of the fundamental questions in biology is how the genotype regulates the phenotype. An increasing number of studies indicate that, in most cases, the effect of a genetic locus on the phenotype is context-dependent, i.e. it is influenced by the genetic background and the environment in which the phenotype is measured. Still, the majority of the studies, in both model organisms and humans, that map the genetic regulation of phenotypic variation in complex traits primarily identify additive loci with independent effects. This does not reflect an absence of the contribution of genetic interactions to phenotypic variation, but instead is a consequence of the technical limitations in mapping gene-gene interactions (GGI) and gene-environment interactions (GEI). Yeast, with its detailed molecular understanding, diverse population genomics and ease of genetic manipulation, is a unique and powerful resource to study the contributions of GGI and GEI in the regulation of phenotypic variation. Here we review studies in yeast that have identified GGI and GEI that regulate phenotypic variation, and discuss the contribution of these findings in explaining missing heritability of complex traits, and how observations from these GGI and GEI studies enhance our understanding of the mechanisms underlying genetic robustness and adaptability that shape the architecture of the genotype-phenotype map.

Different contributions of local- and distant-regulatory changes to transcriptome divergence between stickleback ecotypes

Differential gene expression can play an important role in phenotypic evolution and divergent adaptation. Although differential gene expression can be caused by both local- and distant-regulatory changes, we know little about their relative contribution to transcriptome evolution in natural populations. Here, we conducted expression quantitative trait loci (eQTL) analysis to investigate the genetic architecture underlying transcriptome divergence between marine and stream ecotypes of threespine sticklebacks (Gasterosteus aculeatus). We identified both local and distant eQTLs, some of which constitute hotspots, regions with a disproportionate number of significant eQTLs relative to the genomic background. The majority of local eQTLs including those in the hotspots caused expression changes consistent with the direction of transcriptomic divergence between ecotypes. Genome scan analysis showed that many local eQTLs overlapped with genomic regions of high differentiation. In contrast, nearly half of the distant eQTLs including those in the hotspots caused opposite expression changes, and few overlapped with regions of high differentiation, indicating that distant eQTLs may act as a constraint of transcriptome evolution. Finally, a comparison between two salinity conditions revealed that nearly half of eQTL hotspots were environment specific, suggesting that analysis of genetic architecture in multiple conditions is essential for predicting response to selection.

Natural gene expression variation studies in yeast

The rise of sequence information across different yeast species and strains is driving an increasing number of studies in the emerging field of genomics to associate polymorphic variants, mRNA abundance and phenotypic differences between individuals. Here, we gathered evidence from recent studies covering several layers that define the genotype-phenotype gap, such as mRNA abundance, allele-specific expression and translation efficiency to demonstrate how genetic variants co-evolve and define an individual's genome. Moreover, we exposed several antecedents where inter- and intra-specific studies led to opposite conclusions, probably owing to genetic divergence. Future studies in this area will benefit from the access to a massive array of well-annotated genomes and new sequencing technologies, which will allow the fine breakdown of the complex layers that delineate the genotype-phenotype map. Copyright © 2016 John Wiley & Sons, Ltd.

An integrative analysis of transcriptomic response of ethanol tolerant strains to ethanol in Saccharomyces cerevisiae

The accumulation of ethanol is one of the main environmental stresses that Saccharomyces cerevisiae cells are exposed to in industrial alcoholic beverage and bioethanol production processes. Despite the known impacts of ethanol, the molecular mechanisms underlying ethanol tolerance are still not fully understood. Novel gene targets leading to ethanol tolerance were previously identified via a network approach and the investigations of the deletions of these genes resulted in the improved ethanol tolerance of pmt7Δ/pmt7Δ and yhl042wΔ/yhl042wΔ strains. In the present study, an integrative system based approach was used to investigate the global transcriptional changes in these two ethanol tolerant strains in response to ethanol and hence to elucidate the mechanisms leading to the observed tolerant phenotypes. In addition to strain specific biological processes, a number of common and already reported biological processes were found to be affected in the reference and both ethanol tolerant strains. However, the integrative analysis of the transcriptome with the transcriptional regulatory network and the ethanol tolerance network revealed that each ethanol tolerant strain had a specific organization of the transcriptomic response. Transcription factors around which most important changes occur were determined and active subnetworks in response to ethanol and functional clusters were identified in all strains.

Fitness Trade-Offs Lead to Suppressor Tolerance in Yeast

Genetic variation among individuals within a population provides the raw material for phenotypic diversity upon which natural selection operates. Some given variants can act on multiple standing genomic variations simultaneously and release previously inaccessible phenotypes, leading to increased adaptive potential upon challenging environments. Previously, we identified such a variant related to a tRNA nonsense suppressor in yeast. When introduced into other genetic backgrounds, the suppressor led to an increased population phenotypic variance on various culture conditions, conferring background and environment specific selective advantages. Nonetheless, most isolates are intolerant to the suppressor on rich media due to a severe fitness cost. Here, we found that the tolerance to suppressor is related to a surprising level of fitness outburst, showing a trade-off effect to accommodate the cost of carrying the suppressor. To dissect the genetic basis of such trade-offs, we crossed strains with contrasting tolerance levels on rich media, and analyzed the fitness distribution patterns in the offspring. Combining quantitative tetrad analysis and bulk segregant analysis, we identified two genes, namely MKT1 and RGA1, involved in suppressor tolerance. We showed that alleles from the tolerant parent for both genes conferred a significant gain of fitness, which increased the suppressor tolerance. Our results present a detailed dissection of suppressor tolerance in yeast and provide insights into the molecular basis of trade-offs between fitness and evolutionary potential.

Experimental Evolution of Gene Expression and Plasticity in Alternative Selective Regimes

Little is known of how gene expression and its plasticity evolves as populations adapt to different environmental regimes. Expression is expected to evolve adaptively in all populations but only those populations experiencing environmental heterogeneity are expected to show adaptive evolution of plasticity. We measured the transcriptome in a cadmium-enriched diet and a salt-enriched diet for experimental populations of Drosophila melanogaster that evolved for ~130 generations in one of four selective regimes: two constant regimes maintained in either cadmium or salt diets and two heterogeneous regimes that varied either temporally or spatially between the two diets. For populations evolving in constant regimes, we find a strong signature of counter-gradient evolution; the evolved expression differences between populations adapted to alternative diets is opposite to the plastic response of the ancestral population that is naïve to both diets. Based on expression patterns in the ancestral populations, we identify a set of genes for which we predict selection in heterogeneous regimes to result in increases in plasticity and we find the expected pattern. In contrast, a set of genes where we predicted reduced plasticity did not follow expectation. Nonetheless, both gene sets showed a pattern consistent with adaptive expression evolution in heterogeneous regimes, highlighting the difference between observing “optimal” plasticity and improvements in environment-specific expression. Looking across all genes, there is evidence in all regimes of differences in biased allele expression across environments (“allelic plasticity”) and this is more common among genes with plasticity in total expression.

Different developmental environments change how genes are expressed and what phenotypes are produced. Here we examine how the responsiveness of gene expression to different environments (“expression plasticity”) evolves in populations adapted to constant environments or heterogeneous ones (temporal or spatial heterogeneity) using experimental populations of D. melanogaster. We find the plastic response of the ancestral population that is naïve to both environments is generally opposed by the evolved differences between populations adapted to alternative environments. Populations that live in heterogeneous environments show evidence of adaptive expression evolution in genes predicted to evolve changes in plasticity.

Variation in Indole-3-Acetic Acid Production by Wild Saccharomyces cerevisiae and S. paradoxus Strains from Diverse Ecological Sources and Its Effect on Growth

Phytohormone indole-3-acetic acid (IAA) is the most common naturally occurring and most thoroughly studied plant growth regulator. Microbial synthesis of IAA has long been known. Microbial IAA biosynthesis has been proposed as possibly occurring through multiple pathways, as has been proven in plants. However, the biosynthetic pathways of IAA and the ecological roles of IAA in yeast have not been widely studied. In this study, we investigated the variation in IAA production and its effect on the growth of Saccharomyces cerevisiae and its closest relative Saccharomyces paradoxus yeasts from diverse ecological sources. We found that almost all Saccharomyces yeasts produced IAA when cultured in medium supplemented with the primary precursor of IAA, L-tryptophan (L-Trp). However, when cultured in medium without L-Trp, IAA production was only detected in three strains. Furthermore, exogenous added IAA exerted stimulatory and inhibitory effects on yeast growth. Interestingly, a negative correlation was observed between the amount of IAA production in the yeast cultures and the IAA inhibition ratio of their growth.

Genomic saturation mutagenesis and polygenic analysis identify novel yeast genes affecting ethyl acetate production, a non-selectable polygenic trait

Isolation of mutants in populations of microorganisms has been a valuable tool in experimental genetics for decades. The main disadvantage, however, is the inability of isolating mutants in non-selectable polygenic traits. Most traits of organisms, however, are non-selectable and polygenic, including industrially important properties of microorganisms. The advent of powerful technologies for polygenic analysis of complex traits has allowed simultaneous identification of multiple causative mutations among many thousands of irrelevant mutations. We now show that this also applies to haploid strains of which the genome has been loaded with induced mutations so as to affect as many non-selectable, polygenic traits as possible. We have introduced about 900 mutations into single haploid yeast strains using multiple rounds of EMS mutagenesis, while maintaining the mating capacity required for genetic mapping. We screened the strains for defects in flavor production, an important non-selectable, polygenic trait in yeast alcoholic beverage production. A haploid strain with multiple induced mutations showing reduced ethyl acetate production in semi-anaerobic fermentation, was selected and the underlying quantitative trait loci (QTLs) were mapped using pooled-segregant whole-genome sequence analysis after crossing with an unrelated haploid strain. Reciprocal hemizygosity analysis and allele exchange identified PMA1 and CEM1 as causative mutant alleles and TPS1 as a causative genetic background allele. The case of CEM1 revealed that relevant mutations without observable effect in the haploid strain with multiple induced mutations (in this case due to defective mitochondria) can be identified by polygenic analysis as long as the mutations have an effect in part of the segregants (in this case those that regained fully functional mitochondria). Our results show that genomic saturation mutagenesis combined with complex trait polygenic analysis could be used successfully to identify causative alleles underlying many non-selectable, polygenic traits in small collections of haploid strains with multiple induced mutations.

Using global transcription machinery engineering (gTME) to improve ethanol tolerance of Zymomonas mobilis

BACKGROUND

With the increasing global crude oil crisis and resulting environmental concerns, the production of biofuels from renewable resources has become increasingly important. One of the major challenges faced during the process of biofuel production is the low tolerance of the microbial host towards increasing biofuel concentrations.

RESULTS

Here, we demonstrate that the ethanol tolerance of Zymomonas mobilis can be greatly enhanced through the random mutagenesis of global transcription factor RpoD protein, (σ(70)). Using an enrichment screening, four mutants with elevated ethanol tolerance were isolated from error-prone PCR libraries. All mutants showed significant growth improvement in the presence of ethanol stress when compared to the control strain. After an ethanol (9 %) stress exposure lasting 22 h, the rate of glucose consumption was approximately 1.77, 1.78 and 1.39 g L(-1) h(-1) in the best ethanol-tolerant strain ZM4-mrpoD4, its rebuilt mutant strain ZM4-imrpoD and the control strain, respectively. Our results indicated that both ZM4-mrpoD4 and ZM4-imrpoD consumed glucose at a faster rate after the initial 9 % (v/v) ethanol stress, as nearly 0.64 % of the initial glucose remained after 54 h incubation versus approximately 5.43 % for the control strain. At 9 % ethanol stress, the net ethanol productions by ZM4-mrpoD4 and ZM4-imrpoD during the 30-54 h were 13.0-14.1 g/l versus only 6.6-7.7 g/l for the control strain. The pyruvate decarboxylase activity of ZM4-mrpoD4 was 62.23 and 68.42 U/g at 24 and 48 h, respectively, which were 2.6 and 1.6 times higher than the control strain. After 24 and 48 h of 9 % ethanol stress, the alcohol dehydrogenase activities of ZM4-mrpoD4 were also augmented, showing an approximate 1.4 and 1.3 times increase, respectively, when compared to the control strain. Subsequent quantitative real-time PCR analysis under these stress conditions revealed that the relative expression of pdc in cultured (6 and 24 h) ZM4-mrpoD4 increased by 9.0- and 12.7-fold when compared to control strain.

CONCLUSIONS

Collectively, these results demonstrate that the RpoD mutation can enhance ethanol tolerance in Z. mobilis. Our results also suggested that RpoD may play an important role in resisting high ethanol concentration in Z. mobilis and manipulating RpoD via global transcription machinery engineering (gTME) can provide an alternative and useful approach for strain improvement for complex phenotypes.

How do yeast cells become tolerant to high ethanol concentrations?

The brewer's yeast Saccharomyces cerevisiae displays a much higher ethanol tolerance compared to most other organisms, and it is therefore commonly used for the industrial production of bioethanol and alcoholic beverages. However, the genetic determinants underlying this yeast's exceptional ethanol tolerance have proven difficult to elucidate. In this perspective, we discuss how different types of experiments have contributed to our understanding of the toxic effects of ethanol and the mechanisms and complex genetics underlying ethanol tolerance. In a second part, we summarize the different routes and challenges involved in obtaining superior industrial yeasts with improved ethanol tolerance.

Adaptation to High Ethanol Reveals Complex Evolutionary Pathways

Tolerance to high levels of ethanol is an ecologically and industrially relevant phenotype of microbes, but the molecular mechanisms underlying this complex trait remain largely unknown. Here, we use long-term experimental evolution of isogenic yeast populations of different initial ploidy to study adaptation to increasing levels of ethanol. Whole-genome sequencing of more than 30 evolved populations and over 100 adapted clones isolated throughout this two-year evolution experiment revealed how a complex interplay of de novo single nucleotide mutations, copy number variation, ploidy changes, mutator phenotypes, and clonal interference led to a significant increase in ethanol tolerance. Although the specific mutations differ between different evolved lineages, application of a novel computational pipeline, PheNetic, revealed that many mutations target functional modules involved in stress response, cell cycle regulation, DNA repair and respiration. Measuring the fitness effects of selected mutations introduced in non-evolved ethanol-sensitive cells revealed several adaptive mutations that had previously not been implicated in ethanol tolerance, including mutations in PRT1, VPS70 and MEX67. Interestingly, variation in VPS70 was recently identified as a QTL for ethanol tolerance in an industrial bio-ethanol strain. Taken together, our results show how, in contrast to adaptation to some other stresses, adaptation to a continuous complex and severe stress involves interplay of different evolutionary mechanisms. In addition, our study reveals functional modules involved in ethanol resistance and identifies several mutations that could help to improve the ethanol tolerance of industrial yeasts.

Organisms can evolve resistance to specific stress factors, which allows them to thrive in environments where non-adapted organisms fail to grow. However, the molecular mechanisms that underlie adaptation to complex stress factors that interfere with basic cellular processes are poorly understood. In this study, we reveal how yeast populations adapt to high ethanol concentrations, an ecologically and industrially relevant stress that is still poorly understood. We exposed six independent populations of genetically identical yeast cells to gradually increasing ethanol levels, and we monitored the changes in their DNA sequence over a two-year period. Together with novel computational analyses, we could identify the mutational dynamics and molecular mechanisms underlying increased ethanol resistance. Our results show how adaptation to high ethanol is complex and can be reached through different mutational pathways. Together, our study offers a detailed picture of how populations adapt to a complex continuous stress and identifies several mutations that increase ethanol resistance, which opens new routes to obtain superior biofuel yeast strains.

Population perspectives on functional genomic variation in yeast

Advances in high-throughput sequencing have facilitated large-scale surveys of genomic variation in the budding yeast,Saccharomyces cerevisiae These surveys have revealed extensive sequence variation between yeast strains. However, much less is known about how such variation influences the amount and nature of variation for functional genomic traits within and between yeast lineages. We review population-level studies of functional genomic variation, with a particular focus on how population functional genomic approaches can provide insights into both genome function and the evolutionary process. Although variation in functional genomics phenotypes is pervasive, our understanding of the consequences of this variation, either in physiological or evolutionary terms, is still rudimentary and thus motivates increased attention to appropriate null models. To date, much of the focus of population functional genomic studies has been on gene expression variation, but other functional genomic data types are just as likely to reveal important insights at the population level, suggesting a pressing need for more studies that go beyond transcription. Finally, we discuss how a population functional genomic perspective can be a powerful approach for developing a mechanistic understanding of the processes that link genomic variation to organismal phenotypes through gene networks.

Clade- and species-specific features of genome evolution in the Saccharomycetaceae

Many aspects of the genomes of yeast species in the family Saccharomycetaceae have been well conserved during evolution. They have similar genome sizes, genome contents, and extensive collinearity of gene order along chromosomes. Gene functions can often be inferred reliably by using information from Saccharomyces cerevisiae. Beyond this conservative picture however, there are many instances where a species or a clade diverges substantially from the S. cerevisiae paradigm—for example, by the amplification of a gene family, or by the absence of a biochemical pathway or a protein complex. Here, we review clade-specific features, focusing on genomes sequenced in our laboratory from the post-WGD genera Naumovozyma, Kazachstania and Tetrapisispora, and from the non-WGD species Torulaspora delbrueckii. Examples include the loss of the pathway for histidine synthesis in the cockroach-associated species Tetrapisispora blattae; the presence of a large telomeric GAL gene cluster in To. delbrueckii; losses of the dynein and dynactin complexes in several independent yeast lineages; fragmentation of the MAT locus and loss of the HO gene in Kazachstania africana; and the patchy phylogenetic distribution of RNAi pathway components.

The authors review species-specific evolutionary attributes of yeast genomes.

Temporal expression profiling identifies pathways mediating effect of causal variant on phenotype

Even with identification of multiple causal genetic variants for common human diseases, understanding the molecular processes mediating the causal variants’ effect on the disease remains a challenge. This understanding is crucial for the development of therapeutic strategies to prevent and treat disease. While static profiling of gene expression is primarily used to get insights into the biological bases of diseases, it makes differentiating the causative from the correlative effects difficult, as the dynamics of the underlying biological processes are not monitored. Using yeast as a model, we studied genome-wide gene expression dynamics in the presence of a causal variant as the sole genetic determinant, and performed allele-specific functional validation to delineate the causal effects of the genetic variant on the phenotype. Here, we characterized the precise genetic effects of a functional MKT1 allelic variant in sporulation efficiency variation. A mathematical model describing meiotic landmark events and conditional activation of MKT1 expression during sporulation specified an early meiotic role of this variant. By analyzing the early meiotic genome-wide transcriptional response, we demonstrate an MKT1-dependent role of novel modulators, namely, RTG1/3, regulators of mitochondrial retrograde signaling, and DAL82, regulator of nitrogen starvation, in additively effecting sporulation efficiency. In the presence of functional MKT1 allele, better respiration during early sporulation was observed, which was dependent on the mitochondrial retrograde regulator, RTG3. Furthermore, our approach showed that MKT1 contributes to sporulation independent of Puf3, an RNA-binding protein that steady-state transcription profiling studies have suggested to mediate MKT1-pleiotropic effects during mitotic growth. These results uncover interesting regulatory links between meiosis and mitochondrial retrograde signaling. In this study, we highlight the advantage of analyzing allele-specific transcriptional dynamics of mediating genes. Applications in higher eukaryotes can be valuable for inferring causal molecular pathways underlying complex dynamic processes, such as development, physiology and disease progression.

The causal path from a genetic variant to a complex phenotype such as disease progression is often not known. Studying gene expression variation is one approach to identify the mediating genes, however, it is difficult to distinguish causative from correlative genes. This becomes a challenge especially when studying developmental and physiological traits, since they involve dynamic processes contributing to the variation and only single static expression profiling is performed. As a proof of concept, we addressed this challenge here in yeast, by studying genome-wide gene expression in the presence of the causative polymorphism of MKT1 as the sole genetic variant, during the time phase when it contributes to sporulation efficiency variation. Our analysis during early sporulation identified mitochondrial retrograde signaling and nitrogen starvation as novel regulators, acting additively to regulate sporulation efficiency. Furthermore, we showed that PUF3, a known interactor of MKT1 had an independent role in sporulation. Our results highlight the role of differential mitochondrial signaling for efficient meiosis, providing insights into the factors regulating infertility. In addition, our study has implications for characterizing the molecular effects of causal genetic variants on dynamic biological processes during development and disease progression.

Identification of loci that cause phenotypic variation in diverse species with the reciprocal hemizygosity test

The reciprocal hemizygosity test is a straightforward genetic test that can positively identify genes that have evolved to contribute to a phenotypic difference between strains or between species. The test involves a comparison between hybrids that are genetically identical throughout the genome except at the test locus, which is rendered hemizygous for alternative alleles from the two parental strains. If the two reciprocal hemizygotes display different phenotypes, then the two parental alleles must have evolved. New methods for targeted mutagenesis will allow application of the reciprocal hemizygosity test in many organisms. This review discusses the principles, advantages, and limitations of the test.