The response to single-gene duplication implicates translation as a key vulnerability in aneuploid yeast

Aneuploidy produces myriad consequences in health and disease, yet models of the deleterious effects of chromosome amplification are still widely debated. To distinguish the molecular determinants of aneuploidy stress, we measured the effects of duplicating individual genes in cells with varying chromosome duplications, in wild-type cells and cells sensitized to aneuploidy by deletion of RNA-binding protein Ssd1. We identified gene duplications that are nearly neutral in wild-type euploid cells but significantly deleterious in euploids lacking SSD1 or SSD1+ aneuploid cells with different chromosome duplications. Several of the most deleterious genes are linked to translation; in contrast, duplication of other translational regulators, including eI5Fa Hyp2, benefit ssd1Δ aneuploids over controls. Using modeling of aneuploid growth defects, we propose that the deleterious effects of aneuploidy emerge from an interaction between the cumulative burden of many amplified genes on a chromosome and a subset of duplicated genes that become toxic in that context. Our results suggest that the mechanism behind their toxicity is linked to a key vulnerability in translation in aneuploid cells. These findings provide a perspective on the dual impact of individual genes and overall genomic burden, offering new avenues for understanding aneuploidy and its cellular consequences.

Comparative modeling reveals the molecular determinants of aneuploidy fitness cost in a wild yeast model

Although implicated as deleterious in many organisms, aneuploidy can underlie rapid phenotypic evolution. However, aneuploidy will only be maintained if the benefit outweighs the cost, which remains incompletely understood. To quantify this cost and the molecular determinants behind it, we generated a panel of chromosome duplications in Saccharomyces cerevisiae and applied comparative modeling and molecular validation to understand aneuploidy toxicity. We show that 74-94% of the variance in aneuploid strains' growth rates is explained by the additive cost of genes on each chromosome, measured for single-gene duplications using a genomic library, along with the deleterious contribution of snoRNAs and beneficial effects of tRNAs. Machine learning to identify properties of detrimental gene duplicates provided no support for the balance hypothesis of aneuploidy toxicity and instead identified gene length as the best predictor of toxicity. Our results present a generalized framework for the cost of aneuploidy with implications for disease biology and evolution.

Integrative annotation scores of variants for impact on RNA binding protein activities

MOTIVATION

The ENCODE project generated a large collection of eCLIP-seq RNA binding protein (RBP) profiling data with accompanying RNA-seq transcriptomes of shRNA knockdown of RBPs. These data could have utility in understanding the functional impact of genetic variants, however their potential has not been fully exploited. We implement INCA (Integrative annotation scores of variants for impact on RBP activities) as a multi-step genetic variant scoring approach that leverages the ENCODE RBP data together with ClinVar and integrates multiple computational approaches to aggregate evidence.

RESULTS

INCA evaluates variant impacts on RBP activities by leveraging genotypic differences in cell lines used for eCLIP-seq. We show that INCA provides critical specificity, beyond generic scoring for RBP binding disruption, for candidate variants and their linkage-disequilibrium partners. As a result, it can, on average, augment scoring of 46.2% of the candidate variants beyond generic scoring for RBP binding disruption and aid in variant prioritization for follow-up analysis.

AVAILABILITY AND IMPLEMENTATION

INCA is implemented in R and is available at https://github.com/keleslab/INCA.

Gene-by-environment interactions influence the fitness cost of gene copy-number variation in yeast

Variation in gene copy number can alter gene expression and influence downstream phenotypes; thus copy-number variation provides a route for rapid evolution if the benefits outweigh the cost. We recently showed that genetic background significantly influences how yeast cells respond to gene overexpression, revealing that the fitness costs of copy-number variation can vary substantially with genetic background in a common-garden environment. But the interplay between copy-number variation tolerance and environment remains unexplored on a genomic scale. Here, we measured the tolerance to gene overexpression in four genetically distinct Saccharomyces cerevisiae strains grown under sodium chloride stress. Overexpressed genes that are commonly deleterious during sodium chloride stress recapitulated those commonly deleterious under standard conditions. However, sodium chloride stress uncovered novel differences in strain responses to gene overexpression. West African strain NCYC3290 and North American oak isolate YPS128 are more sensitive to sodium chloride stress than vineyard BC187 and laboratory strain BY4743. Consistently, NCYC3290 and YPS128 showed the greatest sensitivities to overexpression of specific genes. Although most genes were deleterious, hundreds were beneficial when overexpressed-remarkably, most of these effects were strain specific. Few beneficial genes were shared between the sodium chloride-sensitive isolates, implicating mechanistic differences behind their sodium chloride sensitivity. Transcriptomic analysis suggested underlying vulnerabilities and tolerances across strains, and pointed to natural copy-number variation of a sodium export pump that likely contributes to strain-specific responses to overexpression of other genes. Our results reveal extensive strain-by-environment interactions in the response to gene copy-number variation, raising important implications for the accessibility of copy-number variation-dependent evolutionary routes under times of stress.

Advances in S. cerevisiae Engineering for Xylose Fermentation and Biofuel Production: Balancing Growth, Metabolism, and Defense

Genetically engineering microorganisms to produce chemicals has changed the industrialized world. The budding yeast Saccharomyces cerevisiae is frequently used in industry due to its genetic tractability and unique metabolic capabilities. S. cerevisiae has been engineered to produce novel compounds from diverse sugars found in lignocellulosic biomass, including pentose sugars, like xylose, not recognized by the organism. Engineering high flux toward novel compounds has proved to be more challenging than anticipated since simply introducing pathway components is often not enough. Several studies show that the rewiring of upstream signaling is required to direct products toward pathways of interest, but doing so can diminish stress tolerance, which is important in industrial conditions. As an example of these challenges, we reviewed S. cerevisiae engineering efforts, enabling anaerobic xylose fermentation as a model system and showcasing the regulatory interplay's controlling growth, metabolism, and stress defense. Enabling xylose fermentation in S. cerevisiae requires the introduction of several key metabolic enzymes but also regulatory rewiring of three signaling pathways at the intersection of the growth and stress defense responses: the RAS/PKA, Snf1, and high osmolarity glycerol (HOG) pathways. The current studies reviewed here suggest the modulation of global signaling pathways should be adopted into biorefinery microbial engineering pipelines to increase efficient product yields.

Gene-by-environment interactions influence the fitness cost of gene copy-number variation in yeast

Variation in gene copy number can alter gene expression and influence downstream phenotypes; thus copy-number variation (CNV) provides a route for rapid evolution if the benefits outweigh the cost. We recently showed that genetic background significantly influences how yeast cells respond to gene over-expression (OE), revealing that the fitness costs of CNV can vary substantially with genetic background in a common-garden environment. But the interplay between CNV tolerance and environment remains unexplored on a genomic scale. Here we measured the tolerance to gene OE in four genetically distinct Saccharomyces cerevisiae strains grown under sodium chloride (NaCl) stress. OE genes that are commonly deleterious during NaCl stress recapitulated those commonly deleterious under standard conditions. However, NaCl stress uncovered novel differences in strain responses to gene OE. West African strain NCYC3290 and North American oak isolate YPS128 are more sensitive to NaCl stress than vineyard BC187 and laboratory strain BY4743. Consistently, NCYC3290 and YPS128 showed the greatest sensitivities to gene OE. Although most genes were deleterious, hundreds were beneficial when overexpressed - remarkably, most of these effects were strain specific. Few beneficial genes were shared between the NaCl-sensitive isolates, implicating mechanistic differences behind their NaCl sensitivity. Transcriptomic analysis suggested underlying vulnerabilities and tolerances across strains, and pointed to natural CNV of a sodium export pump that likely contributes to strain-specific responses to OE of other genes. Our results reveal extensive strain-by-environment interaction in the response to gene CNV, raising important implications for the accessibility of CNV-dependent evolutionary routes under times of stress.

PKA regulatory subunit Bcy1 couples growth, lipid metabolism, and fermentation during anaerobic xylose growth in Saccharomyces cerevisiae

Organisms have evolved elaborate physiological pathways that regulate growth, proliferation, metabolism, and stress response. These pathways must be properly coordinated to elicit the appropriate response to an ever-changing environment. While individual pathways have been well studied in a variety of model systems, there remains much to uncover about how pathways are integrated to produce systemic changes in a cell, especially in dynamic conditions. We previously showed that deletion of Protein Kinase A (PKA) regulatory subunit BCY1 can decouple growth and metabolism in Saccharomyces cerevisiae engineered for anaerobic xylose fermentation, allowing for robust fermentation in the absence of division. This provides an opportunity to understand how PKA signaling normally coordinates these processes. Here, we integrated transcriptomic, lipidomic, and phospho-proteomic responses upon a glucose to xylose shift across a series of strains with different genetic mutations promoting either coupled or decoupled xylose-dependent growth and metabolism. Together, results suggested that defects in lipid homeostasis limit growth in the bcy1Δ strain despite robust metabolism. To further understand this mechanism, we performed adaptive laboratory evolutions to re-evolve coupled growth and metabolism in the bcy1Δ parental strain. The evolved strain harbored mutations in PKA subunit TPK1 and lipid regulator OPI1, among other genes, and evolved changes in lipid profiles and gene expression. Deletion of the evolved opi1 gene partially reverted the strain's phenotype to the bcy1Δ parent, with reduced growth and robust xylose fermentation. We suggest several models for how cells coordinate growth, metabolism, and other responses in budding yeast and how restructuring these processes enables anaerobic xylose utilization.

The IV International Symposium on Fungal Stress and the XIII International Fungal Biology Conference

For the first time, the International Symposium on Fungal Stress was joined by the XIII International Fungal Biology Conference. The International Symposium on Fungal Stress (ISFUS), always held in Brazil, is now in its fourth edition, as an event of recognized quality in the international community of mycological research. The event held in São José dos Campos, SP, Brazil, in September 2022, featured 33 renowned speakers from 12 countries, including: Austria, Brazil, France, Germany, Ghana, Hungary, México, Pakistan, Spain, Slovenia, USA, and UK. In addition to the scientific contribution of the event in bringing together national and international researchers and their work in a strategic area, it helps maintain and strengthen international cooperation for scientific development in Brazil.

Modeling single-cell phenotypes links yeast stress acclimation to transcriptional repression and pre-stress cellular states

Stress defense and cell growth are inversely related in bulk culture analyses; however, these studies miss substantial cell-to-cell heterogeneity, thus obscuring true phenotypic relationships. Here, we devised a microfluidics system to characterize multiple phenotypes in single yeast cells over time before, during, and after salt stress. The system measured cell and colony size, growth rate, and cell-cycle phase along with nuclear trans-localization of two transcription factors: stress-activated Msn2 that regulates defense genes and Dot6 that represses ribosome biogenesis genes during an active stress response. By tracking cells dynamically, we discovered unexpected discordance between Msn2 and Dot6 behavior that revealed subpopulations of cells with distinct growth properties. Surprisingly, post-stress growth recovery was positively corelated with activation of the Dot6 repressor. In contrast, cells lacking Dot6 displayed slower growth acclimation, even though they grow normally in the absence of stress. We show that wild-type cells with a larger Dot6 response display faster production of Msn2-regulated Ctt1 protein, separable from the contribution of Msn2. These results are consistent with the model that transcriptional repression during acute stress in yeast provides a protective response, likely by redirecting translational capacity to induced transcripts.

Comparative chemical genomic profiling across plant-based hydrolysate toxins reveals widespread antagonism in fitness contributions

The budding yeast Saccharomyces cerevisiae has been used extensively in fermentative industrial processes, including biofuel production from sustainable plant-based hydrolysates. Myriad toxins and stressors found in hydrolysates inhibit microbial metabolism and product formation. Overcoming these stresses requires mitigation strategies that include strain engineering. To identify shared and divergent mechanisms of toxicity and to implicate gene targets for genetic engineering, we used a chemical genomic approach to study fitness effects across a library of S. cerevisiae deletion mutants cultured anaerobically in dozens of individual compounds found in different types of hydrolysates. Relationships in chemical genomic profiles identified classes of toxins that provoked similar cellular responses, spanning inhibitor relationships that were not expected from chemical classification. Our results also revealed widespread antagonistic effects across inhibitors, such that the same gene deletions were beneficial for surviving some toxins but detrimental for others. This work presents a rich dataset relating gene function to chemical compounds, which both expands our understanding of plant-based hydrolysates and provides a useful resource to identify engineering targets.

Comparative functional genomics identifies an iron-limited bottleneck in a Saccharomyces cerevisiae strain with a cytosolic-localized isobutanol pathway

Metabolic engineering strategies have been successfully implemented to improve the production of isobutanol, a next-generation biofuel, in Saccharomyces cerevisiae. Here, we explore how two of these strategies, pathway re-localization and redox cofactor-balancing, affect the performance and physiology of isobutanol producing strains. We equipped yeast with isobutanol cassettes which had either a mitochondrial or cytosolic localized isobutanol pathway and used either a redox-imbalanced (NADPH-dependent) or redox-balanced (NADH-dependent) ketol-acid reductoisomerase enzyme. We then conducted transcriptomic, proteomic and metabolomic analyses to elucidate molecular differences between the engineered strains. Pathway localization had a large effect on isobutanol production with the strain expressing the mitochondrial-localized enzymes producing 3.8-fold more isobutanol than strains expressing the cytosolic enzymes. Cofactor-balancing did not improve isobutanol titers and instead the strain with the redox-imbalanced pathway produced 1.5-fold more isobutanol than the balanced version, albeit at low overall pathway flux. Functional genomic analyses suggested that the poor performances of the cytosolic pathway strains were in part due to a shortage in cytosolic Fe-S clusters, which are required cofactors for the dihydroxyacid dehydratase enzyme. We then demonstrated that this cofactor limitation may be partially recovered by disrupting iron homeostasis with a fra2 mutation, thereby increasing cellular iron levels. The resulting isobutanol titer of the fra2 null strain harboring a cytosolic-localized isobutanol pathway outperformed the strain with the mitochondrial-localized pathway by 1.3-fold, demonstrating that both localizations can support flux to isobutanol.

Breadth and Specificity in Pleiotropic Protein Kinase A Activity and Environmental Responses

Protein Kinase A (PKA) is an essential kinase that is conserved across eukaryotes and plays fundamental roles in a wide range of organismal processes, including growth control, learning and memory, cardiovascular health, and development. PKA mediates these responses through the direct phosphorylation of hundreds of proteins-however, which proteins are phosphorylated can vary widely across cell types and environmental cues, even within the same organism. A major question is how cells enact specificity and precision in PKA activity to mount the proper response, especially during environmental changes in which only a subset of PKA-controlled processes must respond. Research over the years has uncovered multiple strategies that cells use to modulate PKA activity and specificity. This review highlights recent advances in our understanding of PKA signaling control including subcellular targeting, phase separation, feedback control, and standing waves of allosteric regulation. We discuss how the complex inputs and outputs to the PKA network simultaneously pose challenges and solutions in signaling integration and insulation. PKA serves as a model for how the same regulatory factors can serve broad pleiotropic functions but maintain specificity in localized control.

Crabtree/Warburg-like aerobic xylose fermentation by engineered Saccharomyces cerevisiae

Bottlenecks in the efficient conversion of xylose into cost-effective biofuels have limited the widespread use of plant lignocellulose as a renewable feedstock. The yeast Saccharomyces cerevisiae ferments glucose into ethanol with such high metabolic flux that it ferments high concentrations of glucose aerobically, a trait called the Crabtree/Warburg Effect. In contrast to glucose, most engineered S. cerevisiae strains do not ferment xylose at economically viable rates and yields, and they require respiration to achieve sufficient xylose metabolic flux and energy return for growth aerobically. Here, we evolved respiration-deficient S. cerevisiae strains that can grow on and ferment xylose to ethanol aerobically, a trait analogous to the Crabtree/Warburg Effect for glucose. Through genome sequence comparisons and directed engineering, we determined that duplications of genes encoding engineered xylose metabolism enzymes, as well as TKL1, a gene encoding a transketolase in the pentose phosphate pathway, were the causative genetic changes for the evolved phenotype. Reengineered duplications of these enzymes, in combination with deletion mutations in HOG1, ISU1, GRE3, and IRA2, increased the rates of aerobic and anaerobic xylose fermentation. Importantly, we found that these genetic modifications function in another genetic background and increase the rate and yield of xylose-to-ethanol conversion in industrially relevant switchgrass hydrolysate, indicating that these specific genetic modifications may enable the sustainable production of industrial biofuels from yeast. We propose a model for how key regulatory mutations prime yeast for aerobic xylose fermentation by lowering the threshold for overflow metabolism, allowing mutations to increase xylose flux and to redirect it into fermentation products.

Natural variation in the consequences of gene overexpression and its implications for evolutionary trajectories

Copy number variation through gene or chromosome amplification provides a route for rapid phenotypic variation and supports the long-term evolution of gene functions. Although the evolutionary importance of copy-number variation is known, little is understood about how genetic background influences its tolerance. Here, we measured fitness costs of over 4000 overexpressed genes in 15 Saccharomyces cerevisiae strains representing different lineages, to explore natural variation in tolerating gene overexpression (OE). Strain-specific effects dominated the fitness costs of gene OE. We report global differences in the consequences of gene OE, independent of the amplified gene, as well as gene-specific effects that were dependent on the genetic background. Natural variation in the response to gene OE could be explained by several models, including strain-specific physiological differences, resource limitations, and regulatory sensitivities. This work provides new insight on how genetic background influences tolerance to gene amplification and the evolutionary trajectories accessible to different backgrounds.

The role of stress-activated RNA-protein granules in surviving adversity

Severe environmental stress can trigger a plethora of physiological changes and, in the process, significant cytoplasmic reorganization. Stress-activated RNA-protein granules have been implicated in this cellular overhaul by sequestering pre-existing mRNAs and influencing their fates during and after stress acclimation. While the composition and dynamics of stress-activated granule formation has been well studied, their function and impact on RNA-cargo has remained murky. Several recent studies challenge the view that these granules degrade and silence mRNAs present at the onset of stress and instead suggest new roles for these structures in mRNA storage, transit, and inheritance. Here we discuss recent evidence for revised models of stress-activated granule functions and the role of these granules in stress survival and recovery.

Genetic variation in aneuploidy prevalence and tolerance across Saccharomyces cerevisiae lineages

Individuals carrying an aberrant number of chromosomes can vary widely in their expression of aneuploidy phenotypes. A major unanswered question is the degree to which an individual's genetic makeup influences its tolerance of karyotypic imbalance. Here we investigated within-species variation in aneuploidy prevalence and tolerance, using Saccharomyces cerevisiae as a model for eukaryotic biology. We analyzed genotypic and phenotypic variation recently published for over 1,000 S. cerevisiae strains spanning dozens of genetically defined clades and ecological associations. Our results show that the prevalence of chromosome gain and loss varies by clade and can be better explained by differences in genetic background than ecology. The relationships between lineages with high aneuploidy frequencies suggest that increased aneuploidy prevalence emerged multiple times in S. cerevisiae evolution. Separate from aneuploidy prevalence, analyzing growth phenotypes revealed that some genetic backgrounds-such as the European Wine lineage-show fitness costs in aneuploids compared to euploids, whereas other clades with high aneuploidy frequencies show little evidence of major deleterious effects. Our analysis confirms that chromosome gain can produce phenotypic benefits, which could influence evolutionary trajectories. These results have important implications for understanding genetic variation in aneuploidy prevalence in health, disease, and evolution.

CRISpy-Pop: A Web Tool for Designing CRISPR/Cas9-Driven Genetic Modifications in Diverse Populations

CRISPR/Cas9 is a powerful tool for editing genomes, but design decisions are generally made with respect to a single reference genome. With population genomic data becoming available for an increasing number of model organisms, researchers are interested in manipulating multiple strains and lines. CRISpy-pop is a web application that generates and filters guide RNA sequences for CRISPR/Cas9 genome editing for diverse yeast and bacterial strains. The current implementation designs and predicts the activity of guide RNAs against more than 1000 Saccharomyces cerevisiae genomes, including 167 strains frequently used in bioenergy research. Zymomonas mobilis, an increasingly popular bacterial bioenergy research model, is also supported. CRISpy-pop is available as a web application (https://CRISpy-pop.glbrc.org/) with an intuitive graphical user interface. CRISpy-pop also cross-references the human genome to allow users to avoid the selection of guide RNAs with potential biosafety concerns. Additionally, CRISpy-pop predicts the strain coverage of each guide RNA within the supported strain sets, which aids in functional population genetic studies. Finally, we validate how CRISpy-pop can accurately predict the activity of guide RNAs across strains using population genomic data.

Phosphoproteome Response to Dithiothreitol Reveals Unique Versus Shared Features of Saccharomyces cerevisiae Stress Responses

To cope with sudden changes in the external environment, the budding yeast Saccharomyces cerevisiae orchestrates a multifaceted response that spans many levels of physiology. Several studies have interrogated the transcriptome response to endoplasmic reticulum (ER) stress and the role of regulators such as the Ire1 kinase and Hac1 transcription factors. However, less is known about responses to ER stress at other levels of physiology. Here, we used quantitative phosphoproteomics and computational network inference to uncover the yeast phosphoproteome response to the reducing agent dithiothreitol (DTT) and the upstream signaling network that controls it. We profiled wild-type cells and mutants lacking IRE1 or MAPK kinases MKK1 and MKK2, before and at various times after DTT treatment. In addition to revealing downstream targets of these kinases, our inference approach predicted new regulators in the DTT response, including cell-cycle regulator Cdc28 and osmotic-response kinase Rck2, which we validated computationally. Our results also revealed similarities and surprising differences in responses to different stress conditions, especially in the response of protein kinase A targets. These results have implications for the breadth of signaling programs that can give rise to common stress response signatures.

The genetic basis of aneuploidy tolerance in wild yeast

Aneuploidy is highly detrimental during development yet common in cancers and pathogenic fungi – what gives rise to differences in aneuploidy tolerance remains unclear. We previously showed that wild isolates of Saccharomyces cerevisiae tolerate chromosome amplification while laboratory strains used as a model for aneuploid syndromes do not. Here, we mapped the genetic basis to Ssd1, an RNA-binding translational regulator that is functional in wild aneuploids but defective in laboratory strain W303. Loss of SSD1 recapitulates myriad aneuploidy signatures previously taken as eukaryotic responses. We show that aneuploidy tolerance is enabled via a role for Ssd1 in mitochondrial physiology, including binding and regulating nuclear-encoded mitochondrial mRNAs, coupled with a role in mitigating proteostasis stress. Recapitulating ssd1Δ defects with combinatorial drug treatment selectively blocked proliferation of wild-type aneuploids compared to euploids. Our work adds to elegant studies in the sensitized laboratory strain to present a mechanistic understanding of eukaryotic aneuploidy tolerance.

PKA and HOG signaling contribute separable roles to anaerobic xylose fermentation in yeast engineered for biofuel production

Lignocellulosic biomass offers a sustainable source for biofuel production that does not compete with food-based cropping systems. Importantly, two critical bottlenecks prevent economic adoption: many industrially relevant microorganisms cannot ferment pentose sugars prevalent in lignocellulosic medium, leaving a significant amount of carbon unutilized. Furthermore, chemical biomass pretreatment required to release fermentable sugars generates a variety of toxins, which inhibit microbial growth and metabolism, specifically limiting pentose utilization in engineered strains. Here we dissected genetic determinants of anaerobic xylose fermentation and stress tolerance in chemically pretreated corn stover biomass, called hydrolysate. We previously revealed that loss-of-function mutations in the stress-responsive MAP kinase HOG1 and negative regulator of the RAS/Protein Kinase A (PKA) pathway, IRA2, enhances anaerobic xylose fermentation. However, these mutations likely reduce cells' ability to tolerate the toxins present in lignocellulosic hydrolysate, making the strain especially vulnerable to it. We tested the contributions of Hog1 and PKA signaling via IRA2 or PKA negative regulatory subunit BCY1 to metabolism, growth, and stress tolerance in corn stover hydrolysate and laboratory medium with mixed sugars. We found mutations causing upregulated PKA activity increase growth rate and glucose consumption in various media but do not have a specific impact on xylose fermentation. In contrast, mutation of HOG1 specifically increased xylose usage. We hypothesized improving stress tolerance would enhance the rate of xylose consumption in hydrolysate. Surprisingly, increasing stress tolerance did not augment xylose fermentation in lignocellulosic medium in this strain background, suggesting other mechanisms besides cellular stress limit this strain's ability for anaerobic xylose fermentation in hydrolysate.

Rewired cellular signaling coordinates sugar and hypoxic responses for anaerobic xylose fermentation in yeast

Microbes can be metabolically engineered to produce biofuels and biochemicals, but rerouting metabolic flux toward products is a major hurdle without a systems-level understanding of how cellular flux is controlled. To understand flux rerouting, we investigated a panel of Saccharomyces cerevisiae strains with progressive improvements in anaerobic fermentation of xylose, a sugar abundant in sustainable plant biomass used for biofuel production. We combined comparative transcriptomics, proteomics, and phosphoproteomics with network analysis to understand the physiology of improved anaerobic xylose fermentation. Our results show that upstream regulatory changes produce a suite of physiological effects that collectively impact the phenotype. Evolved strains show an unusual co-activation of Protein Kinase A (PKA) and Snf1, thus combining responses seen during feast on glucose and famine on non-preferred sugars. Surprisingly, these regulatory changes were required to mount the hypoxic response when cells were grown on xylose, revealing a previously unknown connection between sugar source and anaerobic response. Network analysis identified several downstream transcription factors that play a significant, but on their own minor, role in anaerobic xylose fermentation, consistent with the combinatorial effects of small-impact changes. We also discovered that different routes of PKA activation produce distinct phenotypes: deletion of the RAS/PKA inhibitor IRA2 promotes xylose growth and metabolism, whereas deletion of PKA inhibitor BCY1 decouples growth from metabolism to enable robust fermentation without division. Comparing phosphoproteomic changes across ira2Δ and bcy1Δ strains implicated regulatory changes linked to xylose-dependent growth versus metabolism. Together, our results present a picture of the metabolic logic behind anaerobic xylose flux and suggest that widespread cellular remodeling, rather than individual metabolic changes, is an important goal for metabolic engineering.

Noncanonical CTD kinases regulate RNA polymerase II in a gene-class-specific manner

Phosphorylation of the carboxyl-terminal domain (CTD) of the largest subunit of RNA polymerase II (Pol II) governs stage-specific interactions with different cellular machines. The CTD consists of Y₁S₂P₃T₄S₅P₆S₇ heptad repeats, and sequential phosphorylations of Ser7, Ser5 and Ser2 occur universally across Pol II-transcribed genes. Phosphorylation of Thr4, however, appears to selectively modulate transcription of specific classes of genes. Here, we identify 10 new Thr4 kinases from different kinase structural groups. Irreversible chemical inhibition of the most active Thr4 kinase, Hrr25, reveals a novel role for this kinase in transcription termination of specific class of noncoding snoRNA genes. Genome-wide profiles of Hrr25 reveal a selective enrichment at 3ʹ regions of noncoding genes that display termination defects. Importantly, phospho-Thr4 marks placed by Hrr25 are recognized by Rtt103, a key component of the termination machinery. Our results suggest that these uncommon CTD kinases selectively place phospho-Thr4 marks to regulate expression of targeted genes.

Natural Variation in the Multidrug Efflux Pump SGE1 Underlies Ionic Liquid Tolerance in Yeast

Imidazolium ionic liquids (IILs) have a range of biotechnological applications, including as pretreatment solvents that extract cellulose from plant biomass for microbial fermentation into sustainable bioenergy. However, residual levels of IILs, such as 1-ethyl-3-methylimidazolium chloride ([C2C1im]Cl), are toxic to biofuel-producing microbes, including the yeast Saccharomyces cerevisiae. S. cerevisiae strains isolated from diverse ecological niches differ in genomic sequence and in phenotypes potentially beneficial for industrial applications, including tolerance to inhibitory compounds present in hydrolyzed plant feedstocks. We evaluated >100 genome-sequenced S. cerevisiae strains for tolerance to [C2C1im]Cl and identified one strain with exceptional tolerance. By screening a library of genomic DNA fragments from the [C2C1im]Cl-tolerant strain for improved IIL tolerance, we identified SGE1, which encodes a plasma membrane multidrug efflux pump, and a previously uncharacterized gene that we named ionic liquid tolerance 1 (ILT1), which encodes a predicted membrane protein. Analyses of SGE1 sequences from our panel of S. cerevisiae strains together with growth phenotypes implicated two single nucleotide polymorphisms (SNPs) that associated with IIL tolerance and sensitivity. We confirmed these phenotypic effects by transferring the SGE1 SNPs into a [C2C1im]Cl-sensitive yeast strain using CRISPR/Cas9 genome editing. Further studies indicated that these SNPs affect Sge1 protein stability and cell surface localization, influencing the amount of toxic IILs that cells can pump out of the cytoplasm. Our results highlight the general potential for discovering useful biotechnological functions from untapped natural sequence variation and provide functional insight into emergent SGE1 alleles with reduced capacities to protect against IIL toxicity.

Decoupling Yeast Cell Division and Stress Defense Implicates mRNA Repression in Translational Reallocation during Stress

Stress tolerance and rapid growth are often competing interests in cells. Upon severe environmental stress, many organisms activate defense systems concurrent with growth arrest. There has been debate as to whether aspects of the stress-activated transcriptome are regulated by stress or an indirect byproduct of reduced proliferation. For example, stressed Saccharomyces cerevisiae cells mount a common gene expression program called the environmental stress response (ESR) [1] comprised of ∼300 induced (iESR) transcripts involved in stress defense and ∼600 reduced (rESR) mRNAs encoding ribosomal proteins (RPs) and ribosome biogenesis factors (RiBi) important for division. Because ESR activation also correlates with reduced growth rate in nutrient-restricted chemostats and prolonged G1 in slow-growing mutants, an alternate proposal is that the ESR is simply a consequence of reduced division [2-5]. A major challenge is that past studies did not separate effects of division arrest and stress defense; thus, the true responsiveness of the ESR-and the purpose of stress-dependent rESR repression in particular-remains unclear. Here, we decoupled cell division from the stress response by following transcriptome, proteome, and polysome changes in arrested cells responding to acute stress. We show that the ESR cannot be explained by changes in growth rate or cell-cycle phase during stress acclimation. Instead, failure to repress rESR transcripts reduces polysome association of induced transcripts, delaying production of their proteins. Our results suggest that stressed cells alleviate competition for translation factors by removing mRNAs and ribosomes from the translating pool, directing translational capacity toward induced transcripts to accelerate protein production.

Network inference reveals novel connections in pathways regulating growth and defense in the yeast salt response

Cells respond to stressful conditions by coordinating a complex, multi-faceted response that spans many levels of physiology. Much of the response is coordinated by changes in protein phosphorylation. Although the regulators of transcriptome changes during stress are well characterized in Saccharomyces cerevisiae, the upstream regulatory network controlling protein phosphorylation is less well dissected. Here, we developed a computational approach to infer the signaling network that regulates phosphorylation changes in response to salt stress. We developed an approach to link predicted regulators to groups of likely co-regulated phospho-peptides responding to stress, thereby creating new edges in a background protein interaction network. We then use integer linear programming (ILP) to integrate wild type and mutant phospho-proteomic data and predict the network controlling stress-activated phospho-proteomic changes. The network we inferred predicted new regulatory connections between stress-activated and growth-regulating pathways and suggested mechanisms coordinating metabolism, cell-cycle progression, and growth during stress. We confirmed several network predictions with co-immunoprecipitations coupled with mass-spectrometry protein identification and mutant phospho-proteomic analysis. Results show that the cAMP-phosphodiesterase Pde2 physically interacts with many stress-regulated transcription factors targeted by PKA, and that reduced phosphorylation of those factors during stress requires the Rck2 kinase that we show physically interacts with Pde2. Together, our work shows how a high-quality computational network model can facilitate discovery of new pathway interactions during osmotic stress.

Cells sense and respond to stressful environments by utilizing complex signaling networks that integrate diverse signals to coordinate a multi-faceted physiological response. Much of this response is controlled by post-translational protein phosphorylation. Although many regulators that mediate changes in protein phosphorylation are known, how these regulators inter-connect in a single regulatory network that can transmit cellular signals is not known. It is also unclear how regulators that promote growth and regulators that activate the stress response interconnect to reorganize resource allocation during stress. Here, we developed an integrated experimental and computational workflow to infer the signaling network that regulates phosphorylation changes during osmotic stress in the budding yeast Saccharomyces cerevisiae. The workflow integrates data measuring protein phosphorylation changes in response to osmotic stress with known physical interactions between yeast proteins from large-scale datasets, along with other information about how regulators recognize their targets. The resulting network suggested new signaling connections between regulators and pathways, including those involved in regulating growth and defense, and predicted new regulators involved in stress defense. Our work highlights the power of using network inference to deliver new insight on how cells coordinate a diverse adaptive strategy to stress.

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

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

Genetic background effects in quantitative genetics: gene-by-system interactions

Proper cell function depends on networks of proteins that interact physically and functionally to carry out physiological processes. Thus, it seems logical that the impact of sequence variation in one protein could be significantly influenced by genetic variants at other loci in a genome. Nonetheless, the importance of such genetic interactions, known as epistasis, in explaining phenotypic variation remains a matter of debate in genetics. Recent work from our lab revealed that genes implicated from an association study of toxin tolerance in Saccharomyces cerevisiae show extensive interactions with the genetic background: most implicated genes, regardless of allele, are important for toxin tolerance in only one of two tested strains. The prevalence of background effects in our study adds to other reports of widespread genetic-background interactions in model organisms. We suggest that these effects represent many-way interactions with myriad features of the cellular system that vary across classes of individuals. Such gene-by-system interactions may influence diverse traits and require new modeling approaches to accurately represent genotype-phenotype relationships across individuals.

Genome-wide association across Saccharomyces cerevisiae strains reveals substantial variation in underlying gene requirements for toxin tolerance

Cellulosic plant biomass is a promising sustainable resource for generating alternative biofuels and biochemicals with microbial factories. But a remaining bottleneck is engineering microbes that are tolerant of toxins generated during biomass processing, because mechanisms of toxin defense are only beginning to emerge. Here, we exploited natural diversity in 165 Saccharomyces cerevisiae strains isolated from diverse geographical and ecological niches, to identify mechanisms of hydrolysate-toxin tolerance. We performed genome-wide association (GWA) analysis to identify genetic variants underlying toxin tolerance, and gene knockouts and allele-swap experiments to validate the involvement of implicated genes. In the process of this work, we uncovered a surprising difference in genetic architecture depending on strain background: in all but one case, knockout of implicated genes had a significant effect on toxin tolerance in one strain, but no significant effect in another strain. In fact, whether or not the gene was involved in tolerance in each strain background had a bigger contribution to strain-specific variation than allelic differences. Our results suggest a major difference in the underlying network of causal genes in different strains, suggesting that mechanisms of hydrolysate tolerance are very dependent on the genetic background. These results could have significant implications for interpreting GWA results and raise important considerations for engineering strategies for industrial strain improvement.

Understanding the genetic architecture of complex traits is important for elucidating the genotype-phenotype relationship. Many studies have sought genetic variants that underlie phenotypic variation across individuals, both to implicate causal variants and to inform on architecture. Here we used genome-wide association analysis to identify genes and processes involved in tolerance of toxins found in plant-biomass hydrolysate, an important substrate for sustainable biofuel production. We found substantial variation in whether or not individual genes were important for tolerance across genetic backgrounds. Whether or not a gene was important in a given strain background explained more variation than the alleleic differences in the gene. These results suggest substantial variation in gene contributions, and perhaps underlying mechanisms, of toxin tolerance.

The Substrates of Nonsense-Mediated mRNA Decay in Caenorhabditis elegans

Nonsense-mediated mRNA decay (NMD) is a conserved pathway that strongly influences eukaryotic gene expression. Inactivating or inhibiting NMD affects the abundance of a substantial fraction of the transcriptome in numerous species. Transcripts whose abundance is altered in NMD-deficient cells may represent either direct substrates of NMD or indirect effects of inhibiting NMD. We present a genome-wide investigation of the direct substrates of NMD in Caenorhabditis elegans Our goals were (i) to identify mRNA substrates of NMD and (ii) to distinguish those mRNAs from others whose abundance is indirectly influenced by the absence of NMD. We previously demonstrated that Upf1p/SMG-2, the central effector of NMD in all studied eukaryotes, preferentially associates with mRNAs that contain premature translation termination codons. We used this preferential association to distinguish direct from indirect effects by coupling immunopurification of Upf1/SMG-2 with high-throughput mRNA sequencing of NMD-deficient mutants and NMD-proficient controls. We identify 680 substrates of NMD, 171 of which contain novel spliced forms that (i) include sequences of annotated introns and (ii) have not been previously documented in the C. elegans transcriptome. NMD degrades unproductively spliced mRNAs with sufficient efficiency in NMD-proficient strains that such mRNAs were not previously known. Two classes of genes are enriched among the identified NMD substrates: (i) mRNAs of expressed pseudogenes and (ii) mRNAs of gene families whose gene number has recently expanded in the C. elegans genome. Our results identify novel NMD substrates and provide a context for understanding NMD's role in normal gene expression and genome evolution.

Proteoform Suite: Software for Constructing, Quantifying, and Visualizing Proteoform Families

We present an open-source, interactive program named Proteoform Suite that uses proteoform mass and intensity measurements from complex biological samples to identify and quantify proteoforms. It constructs families of proteoforms derived from the same gene, assesses proteoform function using gene ontology (GO) analysis, and enables visualization of quantified proteoform families and their changes. It is applied here to reveal systemic proteoform variations in the yeast response to salt stress.

Single-cell RNA sequencing reveals intrinsic and extrinsic regulatory heterogeneity in yeast responding to stress

From bacteria to humans, individual cells within isogenic populations can show significant variation in stress tolerance, but the nature of this heterogeneity is not clear. To investigate this, we used single-cell RNA sequencing to quantify transcript heterogeneity in single Saccharomyces cerevisiae cells treated with and without salt stress to explore population variation and identify cellular covariates that influence the stress-responsive transcriptome. Leveraging the extensive knowledge of yeast transcriptional regulation, we uncovered significant regulatory variation in individual yeast cells, both before and after stress. We also discovered that a subset of cells appears to decouple expression of ribosomal protein genes from the environmental stress response in a manner partly correlated with the cell cycle but unrelated to the yeast ultradian metabolic cycle. Live-cell imaging of cells expressing pairs of fluorescent regulators, including the transcription factor Msn2 with Dot6, Sfp1, or MAP kinase Hog1, revealed both coordinated and decoupled nucleocytoplasmic shuttling. Together with transcriptomic analysis, our results suggest that cells maintain a cellular filter against decoupled bursts of transcription factor activation but mount a stress response upon coordinated regulation, even in a subset of unstressed cells.

Genetically identical cells growing in the same environment can vary in their cellular state and behavior. Such heterogeneity may explain why some cells in an isogenic population can survive sudden severe environmental stress whereas other cells succumb. Cell-to-cell variation in gene expression has been linked to variable stress survival, but how and why transcript levels vary across the transcriptome in single cells is only beginning to emerge. Here, we used single-cell RNA sequencing (scRNA-seq) to measure cell-to-cell heterogeneity in the transcriptome of budding yeast (Saccharomyces cerevisiae). We find surprising patterns of variation across known sets of transcription factor targets, indicating that cells vary in their transcriptome profile both before and after stress exposure. scRNA-seq analysis combined with live-cell imaging of transcription factor activation dynamics revealed some cells in which the stress response was coordinately activated and other cells in which the traditional response was decoupled, suggesting unrecognized regulatory nuances that expand our understanding of stress response and survival.

Multiplexed Sequence-Specific Capture of Chromatin and Mass Spectrometric Discovery of Associated Proteins

Comprehensive understanding of a gene's expression and regulation at the molecular level requires identification of all proteins interacting with the gene. HyCCAPP (Hybridization Capture of Chromatin Associated Proteins for Proteomics) is an approach that uses single-stranded DNA oligonucleotides to capture specific genomic sequences in cross-linked chromatin fragments and identify associated proteins by mass spectrometry. Previous studies have shown HyCCAPP to provide useful information on protein-DNA interactions, revealing the proteins associated with the GAL1-10 region in yeast. We present here a multiplexed version of HyCCAPP. Utilizing a toehold-mediated capture/release strategy, HyCCAPP is targeted to multiple genomic loci in parallel, and the protein binders at each locus are eluted in a programmable and selective fashion. Multiplexed HyCCAPP was applied to four genes (25S rDNA, ARX1, CTT1, and RPL30) in S. cerevisiae under normal and stressed conditions. Capture and release efficiencies and specificities were comparable to those obtained without multiplexing. Using mass spectrometry-based bottom-up proteomics, hundreds of proteins were discovered at each locus in each condition. Statistical analysis revealed 34-88 enriched proteins in each gene capture. Many of these proteins had expected functions, including DNA-related and ribosome biogenesis-associated activities. Multiplexed HyCCAPP provides a useful strategy for the identification of proteins interacting with specific chromatin regions.

Directed Evolution Reveals Unexpected Epistatic Interactions That Alter Metabolic Regulation and Enable Anaerobic Xylose Use by Saccharomyces cerevisiae

The inability of native Saccharomyces cerevisiae to convert xylose from plant biomass into biofuels remains a major challenge for the production of renewable bioenergy. Despite extensive knowledge of the regulatory networks controlling carbon metabolism in yeast, little is known about how to reprogram S. cerevisiae to ferment xylose at rates comparable to glucose. Here we combined genome sequencing, proteomic profiling, and metabolomic analyses to identify and characterize the responsible mutations in a series of evolved strains capable of metabolizing xylose aerobically or anaerobically. We report that rapid xylose conversion by engineered and evolved S. cerevisiae strains depends upon epistatic interactions among genes encoding a xylose reductase (GRE3), a component of MAP Kinase (MAPK) signaling (HOG1), a regulator of Protein Kinase A (PKA) signaling (IRA2), and a scaffolding protein for mitochondrial iron-sulfur (Fe-S) cluster biogenesis (ISU1). Interestingly, the mutation in IRA2 only impacted anaerobic xylose consumption and required the loss of ISU1 function, indicating a previously unknown connection between PKA signaling, Fe-S cluster biogenesis, and anaerobiosis. Proteomic and metabolomic comparisons revealed that the xylose-metabolizing mutant strains exhibit altered metabolic pathways relative to the parental strain when grown in xylose. Further analyses revealed that interacting mutations in HOG1 and ISU1 unexpectedly elevated mitochondrial respiratory proteins and enabled rapid aerobic respiration of xylose and other non-fermentable carbon substrates. Our findings suggest a surprising connection between Fe-S cluster biogenesis and signaling that facilitates aerobic respiration and anaerobic fermentation of xylose, underscoring how much remains unknown about the eukaryotic signaling systems that regulate carbon metabolism.

The yeast Saccharomyces cerevisiae is being genetically engineered to produce renewable biofuels from sustainable plant material. Efficient biofuel production from plant material requires conversion of the complex suite of sugars found in plant material, including the five-carbon sugar xylose. Because it does not efficiently metabolize xylose, S. cerevisiae has been engineered with a minimal set of genes that should overcome this problem; however, additional genetic changes are required for optimal fermentative conversion of xylose into biofuel. Despite extensive knowledge of the regulatory networks controlling glucose metabolism, less is known about the regulation of xylose metabolism and how to rewire these networks for effective biofuel production. Here we report genetic mutations that enabled the conversion of xylose into bioethanol by a previously ineffective yeast strain. By comparing altered protein and metabolite abundance within yeast cells containing these mutations, we determined that the mutations synergistically alter metabolic pathways to improve the rate of xylose conversion. One change in a gene with well-characterized aerobic mitochondrial functions was found to play an unexpected role in anaerobic conversion of xylose into ethanol. The results of this work will allow others to rapidly generate yeast strains for the conversion of xylose into biofuels and other products.

HyCCAPP as a tool to characterize promoter DNA-protein interactions in Saccharomyces cerevisiae

Currently available methods for interrogating DNA-protein interactions at individual genomic loci have significant limitations, and make it difficult to work with unmodified cells or examine single-copy regions without specific antibodies. In this study, we describe a physiological application of the Hybridization Capture of Chromatin-Associated Proteins for Proteomics (HyCCAPP) methodology we have developed. Both novel and known locus-specific DNA-protein interactions were identified at the ENO2 and GAL1 promoter regions of Saccharomyces cerevisiae, and revealed subgroups of proteins present in significantly different levels at the loci in cells grown on glucose versus galactose as the carbon source. Results were validated using chromatin immunoprecipitation. Overall, our analysis demonstrates that HyCCAPP is an effective and flexible technology that does not require specific antibodies nor prior knowledge of locally occurring DNA-protein interactions and can now be used to identify changes in protein interactions at target regions in the genome in response to physiological challenges.

Further support for aneuploidy tolerance in wild yeast and effects of dosage compensation on gene copy-number evolution

In our prior work by Hose et al., we performed a genome-sequencing survey and reported that aneuploidy was frequently observed in wild strains of S. cerevisiae. We also profiled transcriptome abundance in naturally aneuploid isolates compared to isogenic euploid controls and found that 10–30% of amplified genes, depending on the strain and affected chromosome, show lower-than-expected expression compared to gene copy number. In Hose et al., we argued that this gene group is enriched for genes subject to one or more modes of dosage compensation, where mRNA abundance is decreased in response to higher dosage of that gene. A recent manuscript by Torres et al. refutes our prior work. Here, we provide a response to Torres et al., along with additional analysis and controls to support our original conclusions. We maintain that aneuploidy is well tolerated in the wild strains of S. cerevisiae that we studied and that the group of genes enriched for those subject to dosage compensation show unique evolutionary signatures.

DOI:http://dx.doi.org/10.7554/eLife.14409.001

Cells package their DNA into structures called chromosomes. Sometimes when a cell divides, it fails to allocate the right number of chromosomes to each new cell and so they end up with too many or too few chromosomes. The extra copies of the genes on an additional chromosome can be harmful to the cells, because the levels of the proteins encoded by those genes may rise abnormally.

Some organisms counteract the harmful effect of having additional chromosomes through a process called dosage compensation. Proteins are produced using genetic information via two steps: first a gene’s DNA sequence is copied into a molecule of RNA, which is then translated into a protein. Dosage compensation can inactivate single genes or whole chromosomes via various means to ensure that the levels of RNA expressed remain normal, even in the presence of extra genes.

In 2015, researchers from the University of Wisconsin-Madison reported that dosage compensation occurs in wild strains of budding yeast and effectively protects against the harmful effects of having extra chromosomes. However, these findings conflicted with earlier studies of laboratory strains of this yeast, and earlier in 2016, other researchers re-analysed the previous study’s data and challenged its findings.

Now, Gasch et al. – who conducted the work reported in 2015 – provide additional controls and computational experiments that support their original analysis. The latest analysis confirmed that the genes identified in the first study are indeed commonly duplicated in wild yeast populations, yet the expression of these genes remains controlled. This is consistent with a model of dosage compensation, for at least some of duplicated genes.

Gasch et al. believe that part of the difference in interpretation of the data relates to perspective. The challenging researchers tested to see if there was a mechanism of dosage compensation that acted across entire chromosomes, which is known to occur in the case of sex chromosomes in mammals. Gasch et al. on the other hand took a different approach and looked to identify effects at the level of individual genes.

Together, the analyses show that, while there is no evidence for a widespread mechanism, the expression of a select set of genes in wild yeast is consistent with gene-specific dosage compensation. Future work will now undoubtedly test the mechanisms behind the gene-specific effects, and explore why wild yeast strains are more tolerant to extra chromosomes than laboratory strains.

DOI:http://dx.doi.org/10.7554/eLife.14409.002

The Power of Natural Variation for Model Organism Biology

Genetic background effects have long been recognized and, in some cases studied, but they are often viewed as a nuisance by molecular biologists. We suggest that genetic variation currently represents a critical frontier for molecular studies. Human genetics has seen a surge of interest in genetic variation and its contributions to disease, but insights into disease mechanisms are difficult since information about gene function is lacking. By contrast, model organism genetics has excelled at revealing molecular mechanisms of cellular processes, but often de-emphasizes genetic variation and its functional consequences. We argue that model organism biology would benefit from incorporating natural variation, both to capture how well laboratory lines exemplify the species they represent and to inform on molecular processes and their variability. Such a synthesis would also greatly expand the relevance of model systems for studies of complex trait variation, including disease.

A unique ecological niche fosters hybridization of oak-tree and vineyard isolates of Saccharomyces cerevisiae

Differential adaptation to distinct niches can restrict gene flow and promote population differentiation within a species. However, in some cases the distinction between niches can collapse, forming a hybrid niche with features of both environments. We previously reported that distinctions between vineyards and oak soil present an ecological barrier that restricts gene flow between lineages of Saccharomyces cerevisiae. Vineyard isolates are tolerant to stresses associated with grapes while North American oak strains are particularly tolerant to freeze-thaw cycles. Here, we report the isolation of S. cerevisiae strains from Wisconsin cherry trees, which display features common to vineyards (e.g. high sugar concentrations) and frequent freeze-thaw cycles. Genome sequencing revealed that the isolated strains are highly heterozygous and represent recent hybrids of the oak × vineyard lineages. We found that the hybrid strains are phenotypically similar to vineyard strains for some traits, but are more similar to oak strains for other traits. The cherry strains were exceptionally good at growing in cherry juice, raising the possibility that they have adapted to this niche. We performed transcriptome profiling in cherry, oak and vineyard strains and show that the cherry-tree hybrids display vineyard-like or oak-like expression, depending on the gene sets, and in some cases, the expression patterns linked back to shared stress tolerances. Allele-specific expression in these natural hybrids suggested concerted cis-regulatory evolution at sets of functionally regulated genes. Our results raise the possibility that hybridization of the two lineages provides a genetic solution to the thriving in this unique niche.

Comparative genomics of Saccharomyces cerevisiae natural isolates for bioenergy production

Lignocellulosic plant material is a viable source of biomass to produce alternative energy including ethanol and other biofuels. However, several factors—including toxic byproducts from biomass pretreatment and poor fermentation of xylose and other pentose sugars—currently limit the efficiency of microbial biofuel production. To begin to understand the genetic basis of desirable traits, we characterized three strains of Saccharomyces cerevisiae with robust growth in a pretreated lignocellulosic hydrolysate or tolerance to stress conditions relevant to industrial biofuel production, through genome and transcriptome sequencing analysis. All stress resistant strains were highly mosaic, suggesting that genetic admixture may contribute to novel allele combinations underlying these phenotypes. Strain-specific gene sets not found in the lab strain were functionally linked to the tolerances of particular strains. Furthermore, genes with signatures of evolutionary selection were enriched for functional categories important for stress resistance and included stress-responsive signaling factors. Comparison of the strains’ transcriptomic responses to heat and ethanol treatment—two stresses relevant to industrial bioethanol production—pointed to physiological processes that were related to particular stress resistance profiles. Many of the genotype-by-environment expression responses occurred at targets of transcription factors with signatures of positive selection, suggesting that these strains have undergone positive selection for stress tolerance. Our results generate new insights into potential mechanisms of tolerance to stresses relevant to biofuel production, including ethanol and heat, present a backdrop for further engineering, and provide glimpses into the natural variation of stress tolerance in wild yeast strains.

Exploiting the yeast stress-activated signaling network to inform on stress biology and disease signaling

Healthy cells utilize intricate systems to monitor their environment and mount robust responses in the event of cellular stress. Whether stress arises from external insults or defects due to mutation and disease, cells must be able to respond precisely to mount the appropriate defenses. Multi-faceted stress responses are generally coupled with arrest of growth and cell-cycle progression, which both limits the transmission of damaged materials and serves to reallocate limited cellular resources toward defense. Therefore, stress defense versus rapid growth represent competing interests in the cell. How eukaryotic cells set the balance between defense versus proliferation, and in particular knowledge of the regulatory networks that control this decision, are poorly understood. In this perspective, we expand upon our recent work inferring the stress-activated signaling network in budding yeast, which captures pathways controlling stress defense and regulators of growth and cell-cycle progression. We highlight similarities between the yeast and mammalian stress responses and explore how stress-activated signaling networks in yeast can inform on signaling defects in human cancers.

Dosage compensation can buffer copy-number variation in wild yeast

Aneuploidy is linked to myriad diseases but also facilitates organismal evolution. It remains unclear how cells overcome the deleterious effects of aneuploidy until new phenotypes evolve. Although laboratory strains are extremely sensitive to aneuploidy, we show here that aneuploidy is common in wild yeast isolates, which show lower-than-expected expression at many amplified genes. We generated diploid strain panels in which cells carried two, three, or four copies of the affected chromosomes, to show that gene-dosage compensation functions at >30% of amplified genes. Genes subject to dosage compensation are under higher expression constraint in wild populations—but they show elevated rates of gene amplification, suggesting that copy-number variation is buffered at these genes. We find that aneuploidy provides a clear ecological advantage to oak strain YPS1009, by amplifying a causal gene that escapes dosage compensation. Our work presents a model in which dosage compensation buffers gene amplification through aneuploidy to provide a natural, but likely transient, route to rapid phenotypic evolution.

DOI:http://dx.doi.org/10.7554/eLife.05462.001

Evolution is driven by changes to the genes and other genetic information found in the DNA of an organism. These changes might, for example, alter the physical characteristics of the organism, or change how efficiently crucial tasks are carried out inside cells. Whatever the change, if it makes it easier for the organism to survive and reproduce, it is more likely to be passed on to future generations.

DNA is organized inside cells in structures called chromosomes. Most of the cells in animals, plants, and fungi contain two copies of each chromosome. However, sometimes mistakes happen during cell division and extra copies of a chromosome—and hence the genes contained within it—may end up in a cell. These extra copies of genes might help to speed up the rate at which a species evolves, as the ‘spare’ copies are free to adapt to new roles. However, having extra copies of genes can also often be harmful, and in humans can cause genetic disorders such as Down syndrome.

In the laboratory, chromosomes are commonly studied in a species of yeast called Saccharomyces cerevisiae. This species consists of several groups—or strains—that are genetically distinct from each other. Over the years, breeding the yeast for experiments has created laboratory strains that have lost some of the characteristics seen in wild strains. Earlier studies suggested that these cells fail to grow properly if they contain extra copies of chromosomes. Now, Hose et al. have studied nearly 50 wild strains of Saccharomyces cerevisiae. In these, extra copies of chromosomes are commonplace, and seemingly have no detrimental effect on growth. Instead, Hose et al. found that cells with too many copies of a gene use many of those genes less often than would be expected. This process is known as ‘dosage compensation’. This dosage compensation has not been observed in laboratory strains, in part because the extra gene copies make them sickly and hard to study.

Together, the results provide examples of how dosage compensation could help new traits to evolve in a species by reducing the negative effects of duplicated genes. This knowledge may have broad application, from suggesting methods to alleviate human disorders to implicating new ways to engineer useful traits in yeast and other microbes.

DOI:http://dx.doi.org/10.7554/eLife.05462.002

Ecological and Genetic Barriers Differentiate Natural Populations of Saccharomyces cerevisiae

How populations that inhabit the same geographical area become genetically differentiated is not clear. To investigate this, we characterized phenotypic and genetic differences between two populations of Saccharomyces cerevisiae that in some cases inhabit the same environment but show relatively little gene flow. We profiled stress sensitivity in a group of vineyard isolates and a group of oak-soil strains and found several niche-related phenotypes that distinguish the populations. We performed bulk-segregant mapping on two of the distinguishing traits: The vineyard-specific ability to grow in grape juice and oak-specific tolerance to the cell wall damaging drug Congo red. To implicate causal genes, we also performed a chemical genomic screen in the lab-strain deletion collection and identified many important genes that fell under quantitative trait loci peaks. One gene important for growth in grape juice and identified by both the mapping and the screen was SSU1, a sulfite-nitrite pump implicated in wine fermentations. The beneficial allele is generated by a known translocation that we reasoned may also serve as a genetic barrier. We found that the translocation is prevalent in vineyard strains, but absent in oak strains, and presents a postzygotic barrier to spore viability. Furthermore, the translocation was associated with a fitness cost to the rapid growth rate seen in oak-soil strains. Our results reveal the translocation as a dual-function locus that enforces ecological differentiation while producing a genetic barrier to gene flow in these sympatric populations.

Pathway connectivity and signaling coordination in the yeast stress-activated signaling network

Stressed cells coordinate a multi-faceted response spanning many levels of physiology. Yet knowledge of the complete stress-activated regulatory network as well as design principles for signal integration remains incomplete. We developed an experimental and computational approach to integrate available protein interaction data with gene fitness contributions, mutant transcriptome profiles, and phospho-proteome changes in cells responding to salt stress, to infer the salt-responsive signaling network in yeast. The inferred subnetwork presented many novel predictions by implicating new regulators, uncovering unrecognized crosstalk between known pathways, and pointing to previously unknown ‘hubs’ of signal integration. We exploited these predictions to show that Cdc14 phosphatase is a central hub in the network and that modification of RNA polymerase II coordinates induction of stress-defense genes with reduction of growth-related transcripts. We find that the orthologous human network is enriched for cancer-causing genes, underscoring the importance of the subnetwork's predictions in understanding stress biology.

Competence for chemical reprogramming of sexual fate correlates with an intersexual molecular signature in Caenorhabditis elegans

In multicellular organisms, genetic programs guide cells to adopt cell fates as tissues are formed during development, maintained in adults, and repaired after injury. Here we explore how a small molecule in the environment can switch a genetic program from one fate to another. Wild-type Caenorhabditis elegans XX adult hermaphrodites make oocytes continuously, but certain mutant XX adults make sperm instead in an otherwise hermaphrodite soma. Thus, puf-8; lip-1 XX adults make only sperm, but they can be switched from sperm to oocyte production by treatment with a small-molecule MEK inhibitor. To ask whether this chemical reprogramming is common, we tested six XX sperm-only mutants, but found only one other capable of cell fate switching, fbf-1; lip-1. Therefore, reprogramming competence relies on genotype, with only certain mutants capable of responding to the MEK inhibitor with a cell fate change. To gain insight into the molecular basis of competence for chemical reprogramming, we compared polyadenylated transcriptomes of competent and noncompetent XX sperm-only mutants in the absence of the MEK inhibitor and hence in the absence of cell fate reprogramming. Despite their cellular production of sperm, competent mutants were enriched for oogenic messenger RNAs relative to mutants lacking competence for chemical reprogramming. In addition, competent mutants expressed the oocyte-specific protein RME-2, whereas those lacking competence did not. Therefore, mutants competent for reprogramming possess an intersexual molecular profile at both RNA and protein levels. We suggest that this intersexual molecular signature is diagnostic of an intermediate network state that poises the germline tissue for changing its cellular fate in response to environmental cues.

Discovery of Chromatin-Associated Proteins via Sequence-Specific Capture and Mass Spectrometric Protein Identification in Saccharomyces cerevisiae

DNA–protein interactions play critical roles in the control of genome expression and other fundamental processes. An essential element in understanding how these systems function is to identify their molecular components. We present here a novel strategy, Hybridization Capture of Chromatin Associated Proteins for Proteomics (HyCCAPP), to identify proteins that are interacting with any given region of the genome. This technology identifies and quantifies the proteins that are specifically interacting with a genomic region of interest by sequence-specific hybridization capture of the target region from in vivo cross-linked chromatin, followed by mass spectrometric identification and quantification of associated proteins. We demonstrate the utility of HyCCAPP by identifying proteins associated with three multicopy and one single-copy loci in yeast. In each case, a locus-specific pattern of target-associated proteins was revealed. The binding of previously unknown proteins was confirmed by ChIP in 11 of 17 cases. The identification of many previously known proteins at each locus provides strong support for the ability of HyCCAPP to correctly identify DNA-associated proteins in a sequence-specific manner, while the discovery of previously unknown proteins provides new biological insights into transcriptional and regulatory processes at the target locus.

NeuCode labels for relative protein quantification

We describe a synthesis strategy for the preparation of lysine isotopologues that differ in mass by as little as 6 mDa. We demonstrate that incorporation of these molecules into the proteomes of actively growing cells does not affect cellular proliferation, and we discuss how to use the embedded mass signatures (neutron encoding (NeuCode)) for multiplexed proteome quantification by means of high-resolution mass spectrometry. NeuCode SILAC amalgamates the quantitative accuracy of SILAC with the multiplexing of isobaric tags and, in doing so, offers up new opportunities for biological investigation. We applied NeuCode SILAC to examine the relationship between transcript and protein levels in yeast cells responding to environmental stress. Finally, we monitored the time-resolved responses of five signaling mutants in a single 18-plex experiment.

Integrated module and gene-specific regulatory inference implicates upstream signaling networks

Regulatory networks that control gene expression are important in diverse biological contexts including stress response and development. Each gene's regulatory program is determined by module-level regulation (e.g. co-regulation via the same signaling system), as well as gene-specific determinants that can fine-tune expression. We present a novel approach, Modular regulatory network learning with per gene information (MERLIN), that infers regulatory programs for individual genes while probabilistically constraining these programs to reveal module-level organization of regulatory networks. Using edge-, regulator- and module-based comparisons of simulated networks of known ground truth, we find MERLIN reconstructs regulatory programs of individual genes as well or better than existing approaches of network reconstruction, while additionally identifying modular organization of the regulatory networks. We use MERLIN to dissect global transcriptional behavior in two biological contexts: yeast stress response and human embryonic stem cell differentiation. Regulatory modules inferred by MERLIN capture co-regulatory relationships between signaling proteins and downstream transcription factors thereby revealing the upstream signaling systems controlling transcriptional responses. The inferred networks are enriched for regulators with genetic or physical interactions, supporting the inference, and identify modules of functionally related genes bound by the same transcriptional regulators. Our method combines the strengths of per-gene and per-module methods to reveal new insights into transcriptional regulation in stress and development.

Defining flexible vs. inherent promoter architectures: the importance of dynamics and environmental considerations

The degree to which nucleosome positioning regulates transcription is an ongoing debate. To address this question, we recently followed dynamic changes in nucleosome occupancy, transcription factor binding and gene expression in yeast cells responding to oxidative stress. Integrating across these dynamic processes revealed new insights into the functions of nucleosome reorganization. Here, we used our data to address the extent to which upstream promoter architecture is a static feature inherent to specific genes vs. a dynamic platform that changes across conditions. Our results argue that, while some aspects of promoter architecture are fixed across environments, the level to which promoters are "open" or "covered" by nucleosomes depends on the conditions investigated.

An Aspergillus nidulans bZIP response pathway hardwired for defensive secondary metabolism operates through aflR

The eukaryotic bZIP transcription factors are critical players in organismal response to environmental challenges. In fungi, the production of secondary metabolites (SMs) is hypothesized as one of the responses to environmental insults, e.g. attack by fungivorous insects, yet little data to support this hypothesis exists. Here we establish a mechanism of bZIP regulation of SMs through RsmA, a recently discovered YAP-like bZIP protein. RsmA greatly increases SM production by binding to two sites in the Aspergillus nidulans AflR promoter region, a C6 transcription factor known for activating production of the carcinogenic and anti-predation SM, sterigmatocystin. Deletion of aflR in an overexpression rsmA (OE:rsmA) background not only eliminates sterigmatocystin production but also significantly reduces asperthecin synthesis. Furthermore, the fungivore, Folsomia candida, exhibited a distinct preference for feeding on wild type rather than an OE:rsmA strain. RsmA may thus have a critical function in mediating direct chemical resistance against predation. Taken together, these results suggest RsmA represents a bZIP pathway hardwired for defensive SM production.