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.

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.

Next generation biofuels including longer-chain alcohols such as butanol are attractive as renewable, high-energy fuels. A barrier to microbial production of butanols is the increased toxicity compared to ethanol; however, the cellular targets and microbial defense mechanisms remain poorly understood, especially under anaerobic conditions used frequently in industry. Here we took a comparative approach to understand the response of Saccharomyces cerevisiae to 1-butanol, isobutanol, or ethanol, across three genetic backgrounds of varying tolerance in aerobic and anaerobic conditions. We find that strains have different growth properties and alcohol tolerances with and without oxygen availability, as well as unique and common responses to each of the three alcohols. Our results provide evidence for strain-by-alcohol-by-oxygen interactions that moderate how cells respond to alcohol stress.

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.

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 Y1S2P3T4S5P6S7 heptad repeats and sequential phosphorylations of Ser7, Ser5 and Ser2 occur universally at Pol II-transcribed genes. Phosphorylation of Thr4, however, appears to selectively modulate transcription of specific classes of genes. Here, we identify ten 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 place phospho-Thr4 marks to regulate expression of targeted genes.

The Crabtree effect, in which fermentative metabolism is preferred at the expense of respiration, is a hallmark of budding yeast's glucose response and a model for the Warburg effect in human tumors. While the glucose-responsive transcriptional repressors Mig1p and Mig2p play well-characterized roles in the Crabtree effect, little function for the related Mig3p transcription factor has been uncovered, despite numerous investigations of laboratory yeast strains. Here we studied a wild isolate of Saccharomyces cerevisiae to uncover a critical role for Mig3p that has been lost in S288c-derived laboratory strains. We found that Mig3p affects the expression of hundreds of glucose-responsive genes in the oak strain YPS163, both during growth under standard conditions and upon ethanol treatment. Our results suggest that Mig3p may act as a multifunctional activator/repressor that plays separate roles under standard vs. stress conditions and that this function has been largely lost in the lab strains. Population analysis suggests that the lab strain and several wild strains harbor mutations that diminish Mig3p function. Thus, by expanding our attention to multiple genetic backgrounds, we have uncovered an important missing link in a key metabolic response.

In nature, stressful environments often occur in combination or close succession, and thus the ability to prepare for impending stress likely provides a significant fitness advantage. Organisms exposed to a mild dose of stress can become tolerant to what would otherwise be a lethal dose of subsequent stress; however, the mechanism of this acquired stress tolerance is poorly understood. To explore this, we exposed the yeast gene-deletion libraries, which interrogate all essential and non-essential genes, to successive stress treatments and identified genes necessary for acquiring subsequent stress resistance. Cells were exposed to one of three different mild stress pretreatments (salt, DTT, or heat shock) and then challenged with a severe dose of hydrogen peroxide (H(2)O(2)). Surprisingly, there was little overlap in the genes required for acquisition of H(2)O(2) tolerance after different mild-stress pretreatments, revealing distinct mechanisms of surviving H(2)O(2) in each case. Integrative network analysis of these results with respect to protein-protein interactions, synthetic-genetic interactions, and functional annotations identified many processes not previously linked to H(2)O(2) tolerance. We tested and present several models that explain the lack of overlap in genes required for H(2)O(2) tolerance after each of the three pretreatments. Together, this work shows that acquired tolerance to the same severe stress occurs by different mechanisms depending on prior cellular experiences, underscoring the context-dependent nature of stress tolerance.

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.

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.

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.

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.

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.

Organisms simplify the orchestration of gene expression by coregulating genes whose products function together in the cell. Many proteins serve different roles depending on the demands of the organism, and therefore the corresponding genes are often coexpressed with different groups of genes under different situations. This poses a challenge in analyzing whole-genome expression data, because many genes will be similarly expressed to multiple, distinct groups of genes. Because most commonly used analytical methods cannot appropriately represent these relationships, the connections between conditionally coregulated genes are often missed.

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.

Lignocellulosic plant material is a viable source of biomass to produce alternative energy including ethanol and other biofuels. However, several factors—including toxic by products 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.

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.

HIV-1 replication requires myriad interactions between cellular proteins and the viral unspliced RNA. These interactions are important in archetypal RNA processes such as transcription and translation as well as for more specialized functions including alternative splicing and packaging of unspliced genomic RNA into virions. We present here a hybridization capture strategy for purification of unspliced full-length HIV RNA-protein complexes preserved in vivo by formaldehyde crosslinking, and coupled with mass spectrometry to identify HIV RNA-protein interactors in HIV-1 infected cells. One hundred eighty-nine proteins were identified to interact with unspliced HIV RNA including Rev and Gag/Gag-Pol, 24 host proteins previously shown to bind segments of HIV RNA, and over 90 proteins previously shown to impact HIV replication. Further analysis using siRNA knockdown techniques against several of these proteins revealed significant changes to HIV expression. These results demonstrate the utility of the approach for the discovery of host proteins involved in HIV replication. Additionally, because this strategy only requires availability of 30 nucleotides of the HIV-RNA for hybridization with a capture oligonucleotide, it is readily applicable to any HIV system of interest regardless of cell type, HIV-1 virus strain, or experimental perturbation.

The heterochromatin protein 1 (HP1) family is involved in various functions with maintenance of chromatin structure. During murine somatic cell reprogramming, we find that early depletion of HP1γ reduces the generation of induced pluripotent stem cells, while late depletion enhances the process, with a concomitant change from a centromeric to nucleoplasmic localization and elongation-associated histone H3.3 enrichment. Depletion of heterochromatin anchoring protein SENP7 increased reprogramming efficiency to a similar extent as HP1γ, indicating the importance of HP1γ release from chromatin for pluripotency acquisition. HP1γ interacted with OCT4 and DPPA4 in HP1α and HP1β knockouts and in H3K9 methylation depleted H3K9M embryonic stem cell (ESC) lines. HP1α and HP1γ complexes in ESCs differed in association with histones, the histone chaperone CAF1 complex, and specific components of chromatin-modifying complexes such as DPY30, implying distinct functional contributions. Taken together, our results reveal the complex contribution of the HP1 proteins to pluripotency.

A major goal in evolutionary biology is to understand how adaptive evolution has influenced natural variation, but identifying loci subject to positive selection has been a challenge. Here we present the adaptive loss of a pair of paralogous genes in specific Saccharomyces cerevisiae subpopulations. We mapped natural variation in freeze-thaw tolerance to two water transporters, AQY1 and AQY2, previously implicated in freeze-thaw survival. However, whereas freeze-thaw-tolerant strains harbor functional aquaporin genes, the set of sensitive strains lost aquaporin function at least 6 independent times. Several genomic signatures at AQY1 and/or AQY2 reveal low variation surrounding these loci within strains of the same haplotype, but high variation between strain groups. This is consistent with recent adaptive loss of aquaporins in subgroups of strains, leading to incipient balancing selection. We show that, although aquaporins are critical for surviving freeze-thaw stress, loss of both genes provides a major fitness advantage on high-sugar substrates common to many strains' natural niche. Strikingly, strains with non-functional alleles have also lost the ancestral requirement for aquaporins during spore formation. Thus, the antagonistic effect of aquaporin function-providing an advantage in freeze-thaw tolerance but a fitness defect for growth in high-sugar environments-contributes to the maintenance of both functional and nonfunctional alleles in S. cerevisiae. This work also shows that gene loss through multiple missense and nonsense mutations, hallmarks of pseudogenization presumed to emerge after loss of constraint, can arise through positive selection.

Yeast responding to stress activate a large gene expression program called the Environmental Stress Response that consists of approximately 600 repressed genes and approximately 300 induced genes. Numerous factors are implicated in regulating subsets of Environmental Stress Response genes; however, a complete picture of Environmental Stress Response regulation remains unclear. We investigated the role of the histone deacetylase Rpd3p, previously linked to the upstream regions of many Environmental Stress Response genes, in producing Environmental Stress Response gene expression changes in response to stress.