Yeast genomic expression studies using DNA microarrays

Commensal Yeast Promotes Salmonella Typhimurium Virulence

Enteric pathogens engage in complex interactions with the host and the resident microbiota to establish gut colonization. Although mechanistic interactions between enteric pathogens and bacterial commensals have been extensively studied, whether and how commensal fungi affect pathogenesis of enteric infections remains largely unknown. Here we show that colonization with the common human gut commensal fungus Candida albicans worsened infections with the enteric pathogen Salmonella enterica serovar Typhimurium. Presence of C. albicans in the mouse gut increased Salmonella cecum colonization and systemic dissemination. We investigated the underlying mechanism and found that Salmonella binds to C. albicans via Type 1 fimbriae and uses its Type 3 Secretion System (T3SS) to deliver effector proteins into C. albicans . A specific effector, SopB, was sufficient to manipulate C. albicans metabolism, triggering increased arginine biosynthesis in C. albicans and the release of millimolar amounts of arginine into the extracellular environment. The released arginine, in turn, induced T3SS expression in Salmonella , increasing its invasion of epithelial cells. C. albicans deficient in arginine production was unable to increase Salmonella virulence in vitro or in vivo . In addition to modulating pathogen invasion, arginine also directly influenced the host response to infection. Arginine-producing C. albicans dampened the inflammatory response during Salmonella infection, whereas C. albicans deficient in arginine production did not. Arginine supplementation in the absence of C. albicans increased the systemic spread of Salmonella and decreased the inflammatory response, phenocopying the presence of C. albicans . In summary, we identified C. albicans colonization as a susceptibility factor for disseminated Salmonella infection, and arginine as a central metabolite in the cross-kingdom interaction between fungi, bacteria, and host.

Premature aging in aneuploid yeast is caused in part by aneuploidy-induced defects in Ribosome Quality Control

Premature aging is a hallmark of Down syndrome, caused by trisomy of human chromosome 21, but the reason is unclear and difficult to study in humans. We used an aneuploid model in wild yeast to show that chromosome amplification disrupts nutrient-induced cell-cycle arrest, quiescence entry, and healthy aging, across genetic backgrounds and amplified chromosomes. We discovered that these defects are due in part to aneuploidy-induced dysfunction in Ribosome Quality Control (RQC). Compared to euploids, aneuploids entering quiescence display aberrant ribosome profiles, accumulate RQC intermediates, and harbor an increased load of protein aggregates. Although they have normal proteasome capacity, aneuploids show signs of ubiquitin dysregulation, which impacts cyclin abundance to disrupt arrest. Remarkably, inducing ribosome stalling in euploids produces similar aberrations, while up-regulating limiting RQC subunits or proteins in ubiquitin metabolism alleviates many of the aneuploid defects. Our results provide implications for other aneuploidy disorders including Down syndrome.

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.

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.

Understanding the differences in 2G ethanol fermentative scales through omics data integration

In this work, we evaluated the fermentative performance and metabolism modifications of a second generation (2G) industrial yeast by comparing an industrial condition during laboratory and industrial scale fermentations. Fermentations were done using industrial lignocellulosic hydrolysate and a synthetic medium containing inhibitors and analyses were carried out through transcriptomics and proteomics of these experimental conditions. We found that fermentation profiles were very similar, but there was an increase in xylose consumption rate during fermentations using synthetic medium when compared to lignocellulosic hydrolysate, likely due to the presence of unknown growth inhibitors contained in the hydrolysate. We also evaluated the bacterial community composition of the industrial fermentation setting and found that the presence of homofermentative and heterofermentative bacteria did not significantly change the performance of yeast fermentation. In parallel, temporal differentially expressed genes (tDEG) showed differences in gene expression profiles between compared conditions, including heat shocks and the presence of up-regulated genes from the TCA cycle during anaerobic xylose fermentation. Thus, we indicate HMF as a possible electron acceptor in this rapid respiratory process performed by yeast, in addition to demonstrating the importance of culture medium for the performance of yeast within industrial fermentation processes, highlighting the uniquenesses according to scales.

Differential activation mechanisms of two isoforms of Gcr1 transcription factor generated from spliced and un-spliced transcripts in Saccharomyces cerevisiae

Gcr1, an important transcription factor for glycolytic genes in Saccharomyces cerevisiae, was recently revealed to have two isoforms, Gcr1^U and Gcr1^S, produced from un-spliced and spliced transcripts, respectively. In this study, by generating strains expressing only Gcr1^U or Gcr1^S using the CRISPR/Cas9 system, we elucidate differential activation mechanisms of these two isoforms. The Gcr1^U monomer forms an active complex with its coactivator Gcr2 homodimer, whereas Gcr1^S acts as a homodimer without Gcr2. The USS domain, 55 residues at the N-terminus existing only in Gcr1^U, inhibits dimerization of Gcr1^U and even acts in trans to inhibit Gcr1^S dimerization. The Gcr1^S monomer inhibits the metabolic switch from fermentation to respiration by directly binding to the ALD4 promoter, which can be restored by overexpression of the ALD4 gene, encoding a mitochondrial aldehyde dehydrogenase required for ethanol utilization. Gcr1^U and Gcr1^S regulate almost the same target genes, but show unique activities depending on growth phase, suggesting that these isoforms play differential roles through separate activation mechanisms depending on environmental conditions.

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.

Comparison of RNA isolation methods on RNA-Seq: implications for differential expression and meta-analyses

BACKGROUND

The increasing number of transcriptomic datasets has allowed for meta-analyses, which can be valuable due to their increased statistical power. However, meta-analyses can be confounded by so-called "batch effects," where technical variation across different batches of RNA-seq experiments can clearly produce spurious signals of differential expression and reduce our power to detect true differences. While batch effects can sometimes be accounted for, albeit with caveats, a better strategy is to understand their sources to better avoid them. In this study, we examined the effects of RNA isolation method as a possible source of batch effects in RNA-seq design.

RESULTS

Based on the different chemistries of "classic" hot phenol extraction of RNA compared to common commercial RNA isolation kits, we hypothesized that specific mRNAs may be preferentially extracted depending upon method, which could masquerade as differential expression in downstream RNA-seq analyses. We tested this hypothesis using the Saccharomyces cerevisiae heat shock response as a well-validated environmental response. Comparing technical replicates that only differed in RNA isolation method, we found over one thousand transcripts that appeared "differentially" expressed when comparing hot phenol extraction with the two kits. Strikingly, transcripts with higher abundance in the phenol-extracted samples were enriched for membrane proteins, suggesting that indeed the chemistry of hot phenol extraction better solubilizes those species of mRNA.

CONCLUSIONS

Within a self-contained experimental batch (e.g. control versus treatment), the method of RNA isolation had little effect on the ability to identify differentially expressed transcripts. However, we suggest that researchers performing meta-analyses across different experimental batches strongly consider the RNA isolation methods for each experiment.

Wolbachia pipientis grows in Saccharomyces cerevisiae evoking early death of the host and deregulation of mitochondrial metabolism

Wolbachia sp. has colonized over 70% of insect species, successfully manipulating host fertility, protein expression, lifespan, and metabolism. Understanding and engineering the biochemistry and physiology of Wolbachia holds great promise for insect vector-borne disease eradication. Wolbachia is cultured in cell lines, which have long duplication times and are difficult to manipulate and study. The yeast strain Saccharomyces cerevisiae W303 was used successfully as an artificial host for Wolbachia wAlbB. As compared to controls, infected yeast lost viability early, probably as a result of an abnormally high mitochondrial oxidative phosphorylation activity observed at late stages of growth. No respiratory chain proteins from Wolbachia were detected, while several Wolbachia F1 F0 -ATPase subunits were revealed. After 5 days outside the cell, Wolbachia remained fully infective against insect cells.

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.

Genotype-by-Environment-by-Environment Interactions in the Saccharomyces cerevisiae Transcriptomic Response to Alcohols and Anaerobiosis

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.

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.

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

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

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

GSF2 deletion increases lactic acid production by alleviating glucose repression in Saccharomyces cerevisiae

Improving lactic acid (LA) tolerance is important for cost-effective microbial production of LA under acidic fermentation conditions. Previously, we generated LA-tolerant D-LA-producing S. cerevisiae strain JHY5310 by laboratory adaptive evolution of JHY5210. In this study, we performed whole genome sequencing of JHY5310, identifying four loss-of-function mutations in GSF2, SYN8, STM1, and SIF2 genes, which are responsible for the LA tolerance of JHY5310. Among the mutations, a nonsense mutation in GSF2 was identified as the major contributor to the improved LA tolerance and LA production in JHY5310. Deletion of GSF2 in the parental strain JHY5210 significantly improved glucose uptake and D-LA production levels, while derepressing glucose-repressed genes including genes involved in the respiratory pathway. Therefore, more efficient generation of ATP and NAD⁺ via respiration might rescue the growth defects of the LA-producing strain, where ATP depletion through extensive export of lactate and proton is one of major reasons for the impaired growth. Accordingly, alleviation of glucose repression by deleting MIG1 or HXK2 in JHY5210 also improved D-LA production. GSF2 deletion could be applied to various bioprocesses where increasing biomass yield or respiratory flux is desirable.

Leveraging Genetic-Background Effects in Saccharomyces cerevisiae To Improve Lignocellulosic Hydrolysate Tolerance

A major obstacle to sustainable lignocellulosic biofuel production is microbe inhibition by the combinatorial stresses in pretreated plant hydrolysate. Chemical biomass pretreatment releases a suite of toxins that interact with other stressors, including high osmolarity and temperature, which together can have poorly understood synergistic effects on cells. Improving tolerance in industrial strains has been hindered, in part because the mechanisms of tolerance reported in the literature often fail to recapitulate in other strain backgrounds. Here, we explored and then exploited variations in stress tolerance, toxin-induced transcriptomic responses, and fitness effects of gene overexpression in different Saccharomyces cerevisiae (yeast) strains to identify genes and processes linked to tolerance of hydrolysate stressors. Using six different S. cerevisiae strains that together maximized phenotypic and genetic diversity, first we explored transcriptomic differences between resistant and sensitive strains to identify common and strain-specific responses. This comparative analysis implicated primary cellular targets of hydrolysate toxins, secondary effects of defective defense strategies, and mechanisms of tolerance. Dissecting the responses to individual hydrolysate components across strains pointed to synergistic interactions between osmolarity, pH, hydrolysate toxins, and nutrient composition. By characterizing the effects of high-copy gene overexpression in three different strains, we revealed the breadth of the background-specific effects of gene fitness contributions in synthetic hydrolysate. Our approach identified new genes for engineering improved stress tolerance in diverse strains while illuminating the effects of genetic background on molecular mechanisms.

IMPORTANCE

Recent studies on natural variation within Saccharomyces cerevisiae have uncovered substantial phenotypic diversity. Here, we took advantage of this diversity, using it as a tool to infer the effects of combinatorial stress found in lignocellulosic hydrolysate. By comparing sensitive and tolerant strains, we implicated primary cellular targets of hydrolysate toxins and elucidated the physiological states of cells when exposed to this stress. We also explored the strain-specific effects of gene overexpression to further identify strain-specific responses to hydrolysate stresses and to identify genes that improve hydrolysate tolerance independent of strain background. This study underscores the importance of studying multiple strains to understand the effects of hydrolysate stress and provides a method to find genes that improve tolerance across strain backgrounds.

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.

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.

CK2-dependent inhibitory phosphorylation is relieved by Ppt1 phosphatase for the ethanol stress-specific activation of Hsf1 in Saccharomyces cerevisiae

Ethanol, the major fermentation product of Saccharomyces cerevisiae, has long been known as an inducer of heat shock response, but the underlying mechanisms by which ethanol activates heat shock transcription factor (HSF) are not well understood. We demonstrate that CK2-dependent phosphorylation on S608 is an ethanol stress-specific repression mechanism of Hsf1, which does not affect the basal or heat-induced activity of Hsf1. This repression is relieved by dephosphorylation by Ppt1 which directly interacts with Hsf1 via its tetratricopeptide repeat (TPR) domain. In response to ethanol stress, PPT1 deletion and CK2 overexpression exert synergistic inhibitory effects on Hsf1 activation, whereas Hsf1(S608A) mutant shows enhanced activation. Therefore, regulation of the Hsf1 S608 phosphorylation status by reciprocal actions of CK2 and Ppt1 might play an important role to determine Hsf1 sensitivity towards ethanol stress.

Natural variation in the yeast glucose-signaling network reveals a new role for the Mig3p transcription factor

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.

Cellular memory of acquired stress resistance in Saccharomyces cerevisiae

Cellular memory of past experiences has been observed in several organisms and across a variety of experiences, including bacteria "remembering" prior nutritional status and amoeba "learning" to anticipate future environmental conditions. Here, we show that Saccharomyces cerevisiae maintains a multifaceted memory of prior stress exposure. We previously demonstrated that yeast cells exposed to a mild dose of salt acquire subsequent tolerance to severe doses of H(2)O(2). We set out to characterize the retention of acquired tolerance and in the process uncovered two distinct aspects of cellular memory. First, we found that H(2)O(2) resistance persisted for four to five generations after cells were removed from the prior salt treatment and was transmitted to daughter cells that never directly experienced the pretreatment. Maintenance of this memory did not require nascent protein synthesis after the initial salt pretreatment, but rather required long-lived cytosolic catalase Ctt1p that was synthesized during salt exposure and then distributed to daughter cells during subsequent cell divisions. In addition to and separable from the memory of H(2)O(2) resistance, these cells also displayed a faster gene-expression response to subsequent stress at >1000 genes, representing transcriptional memory. The faster gene-expression response requires the nuclear pore component Nup42p and serves an important function by facilitating faster reacquisition of H(2)O(2) tolerance after a second cycle of salt exposure. Memory of prior stress exposure likely provides a significant advantage to microbial populations living in ever-changing environments.

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.

Dynamic changes in nucleosome occupancy are not predictive of gene expression dynamics but are linked to transcription and chromatin regulators

The response to stressful stimuli requires rapid, precise, and dynamic gene expression changes that must be coordinated across the genome. To gain insight into the temporal ordering of genome reorganization, we investigated dynamic relationships between changing nucleosome occupancy, transcription factor binding, and gene expression in Saccharomyces cerevisiae yeast responding to oxidative stress. We applied deep sequencing to nucleosomal DNA at six time points before and after hydrogen peroxide treatment and revealed many distinct dynamic patterns of nucleosome gain and loss. The timing of nucleosome repositioning was not predictive of the dynamics of downstream gene expression change but instead was linked to nucleosome position relative to transcription start sites and specific cis-regulatory elements. We measured genome-wide binding of the stress-activated transcription factor Msn2p over time and found that Msn2p binds different loci with different dynamics. Nucleosome eviction from Msn2p binding sites was common across the genome; however, we show that, contrary to expectation, nucleosome loss occurred after Msn2p binding and in fact required Msn2p. This negates the prevailing model that nucleosomes obscuring Msn2p sites regulate DNA access and must be lost before Msn2p can bind DNA. Together, these results highlight the complexities of stress-dependent chromatin changes and their effects on gene expression.

Comparative genomics of xylose-fermenting fungi for enhanced biofuel production

Cellulosic biomass is an abundant and underused substrate for biofuel production. The inability of many microbes to metabolize the pentose sugars abundant within hemicellulose creates specific challenges for microbial biofuel production from cellulosic material. Although engineered strains of Saccharomyces cerevisiae can use the pentose xylose, the fermentative capacity pales in comparison with glucose, limiting the economic feasibility of industrial fermentations. To better understand xylose utilization for subsequent microbial engineering, we sequenced the genomes of two xylose-fermenting, beetle-associated fungi, Spathaspora passalidarum and Candida tenuis. To identify genes involved in xylose metabolism, we applied a comparative genomic approach across 14 Ascomycete genomes, mapping phenotypes and genotypes onto the fungal phylogeny, and measured genomic expression across five Hemiascomycete species with different xylose-consumption phenotypes. This approach implicated many genes and processes involved in xylose assimilation. Several of these genes significantly improved xylose utilization when engineered into S. cerevisiae, demonstrating the power of comparative methods in rapidly identifying genes for biomass conversion while reflecting on fungal ecology.

A dynamic model of proteome changes reveals new roles for transcript alteration in yeast

By characterizing dynamic changes in yeast protein abundance following osmotic shock, this study shows that the correlation between protein and mRNA differs for transcripts that increase versus decrease in abundance, and reveals physiological reasons for these differences.

The correlation between protein and mRNA change is very high at transcripts that increase in abundance, but negligible at reduced transcripts following NaCl shock.

Modeling and experimental data suggest that reducing levels of high-abundance transcripts helps to direct translational machinery to newly made transcripts.

The transient burst of transcript increase serves to accelerate changes in protein abundance.

Post-transcriptional regulation of protein abundance is pervasive, although most of the variance in protein change is explained by changes in mRNA abundance.

Natural microenvironments change rapidly, and living creatures must respond quickly and efficiently to thrive within this flux. At all cellular levels—signaling, transcription, translation, metabolism, cell growth, and division—the response is dynamic and coordinated. Some aspects of this response, such as dynamic changes of the transcriptome, are well understood. But other aspects, like the response of the proteome, have remained obscured primarily because of previous limitations in technology. Without coordinated time-course data, it has remained impossible to correctly characterize the correlations and dependencies between these two essential levels of cell biology.

This work presents an extended picture of the coordinated response of the transcriptome and proteome as cells respond to an abrupt environmental change. To assay proteomic dynamics, we developed a strategy for large-scale, multiplexed quantitation using isobaric tags and high mass accuracy mass spectrometry. This sensitive yet efficient platform allows for the expedient collection of quantitative time-course proteomic data at six time points, sufficiently reproducible to permit meaningful interpretation of variation across biological replicates. Time-course transcriptome data were generated from paired biological samples, allowing us to examine the relationships between changes in mRNA and protein for each gene in terms of direction and intensity, as well as the characteristics of the temporal profiles for each gene.

It was immediately obvious that a single measure of correlation across the entire data set was a meaningless metric. We therefore analyzed relationships between mRNA and protein for different subsets of data. In response to osmotic shock, hundreds of transcripts are highly induced, and their temporal pattern reveals a transient peak of maximal induction, which resolves into a new elevated level as cells acclimate (Figure 2). For this group of genes, there is extremely high correlation between peak mRNA change and protein change (R²∼0.8). But the dynamics of the molecules differ: while mRNA levels transiently overshoot their final levels, proteins gradually rise in abundance toward their new, elevated state. We observed, however, that a measure of efficiency connects the two profiles. The time it takes for a protein to acclimate to its new state correlates with the magnitude of the excess mRNA induction. Thus, the cell imparts an urgency to protein induction by transiently producing excess transcript.

The most surprising result, however, involves transcripts that decrease in abundance. In response to osmotic shock, the cell transiently reduces over 600 transcripts, many of which are among the most highly expressed in unstressed cells. But protein levels for these genes remain, for the most part, almost completely unchanged. The stark absence of protein repression is independent of basal protein abundance, independent of reported protein half-lives, reproducible across biological replicates, and validated by quantitative western blots. Furthermore, since we do detect a handful of proteins whose abundance is significantly reduced, our technology is capable of identifying protein loss. Thus, we conclude that transcript reduction serves another purpose besides reducing protein levels.

To explore alternate interpretations of the consequence of transcriptional repression, we devised a mass-action kinetic model, which describes protein changes based on mRNA dynamics in the context of transient changes in the rates of cell division. The model successfully recapitulated the observed data, allowing us to alter modeling parameters to test various hypotheses.

In response to osmotic shock, overall rates of translation temporarily decrease and cell growth transiently arrests before resuming at a slower rate. We reasoned that mRNA reduction might lower the rate of new protein synthesis, but that retarded production is balanced by reduced cell division. We explored both aspects of this logic with our model.

As expected, removing cell division from our model led to a calculated decrease of protein levels, indicating that reduced growth is necessary for maintaining protein levels. However, when we computationally held mRNA levels stable and calculated protein levels in the absence of mRNA repression, we did not find the expected increase in protein abundance.

We then considered the possibility that one function of the regulated repression of these highly abundant transcripts was to liberate proteins essential for translation, such as ribosomes or translation initiation factors. To explore this, we examined a mutant lacking the Dot6p/Tod6p transcriptional repressors, which fails to properly repress ∼250 genes in response to osmotic shock. In the wild type, the mRNA for a Dot6p/Tod6p target (ARX1) decreased seven-fold, and the remaining transcript was generally unassociated with poly-ribosomes. In the mutant, however, the mRNA levels were reduced only two-fold, while the remaining transcript continued to bind ribosomes. Therefore, failure to reduce transcript levels led to a persistent association with poly-ribosomes, thereby consuming translational machinery.

Our hypothesis is, therefore, that widespread changes in the transcriptome promote efficient translation of new proteins. Transcript increase serves to increase abundance of the encoded proteins, while reduction of some of the most abundant and highly translated mRNAs supports this project by liberating translational capacity. While it is not clear what factors are the limiting elements, it is clear that a full picture of cellular biology requires exploring the dynamics of the cellular response.

The transcriptome and proteome change dynamically as cells respond to environmental stress; however, prior proteomic studies reported poor correlation between mRNA and protein, rendering their relationships unclear. To address this, we combined high mass accuracy mass spectrometry with isobaric tagging to quantify dynamic changes in ∼2500 Saccharomyces cerevisiae proteins, in biological triplicate and with paired mRNA samples, as cells acclimated to high osmolarity. Surprisingly, while transcript induction correlated extremely well with protein increase, transcript reduction produced little to no change in the corresponding proteins. We constructed a mathematical model of dynamic protein changes and propose that the lack of protein reduction is explained by cell-division arrest, while transcript reduction supports redistribution of translational machinery. Furthermore, the transient ‘burst' of mRNA induction after stress serves to accelerate change in the corresponding protein levels. We identified several classes of post-transcriptional regulation, but show that most of the variance in protein changes is explained by mRNA. Our results present a picture of the coordinated physiological responses at the levels of mRNA, protein, protein-synthetic capacity, and cellular growth.

Exploiting natural variation in Saccharomyces cerevisiae to identify genes for increased ethanol resistance

Ethanol production from lignocellulosic biomass holds promise as an alternative fuel. However, industrial stresses, including ethanol stress, limit microbial fermentation and thus prevent cost competitiveness with fossil fuels. To identify novel engineering targets for increased ethanol tolerance, we took advantage of natural diversity in wild Saccharomyces cerevisiae strains. We previously showed that an S288c-derived lab strain cannot acquire higher ethanol tolerance after a mild ethanol pretreatment, which is distinct from other stresses. Here, we measured acquired ethanol tolerance in a large panel of wild strains and show that most strains can acquire higher tolerance after pretreatment. We exploited this major phenotypic difference to address the mechanism of acquired ethanol tolerance, by comparing the global gene expression response to 5% ethanol in S288c and two wild strains. Hundreds of genes showed variation in ethanol-dependent gene expression across strains. Computational analysis identified several transcription factor modules and known coregulated genes as differentially expressed, implicating genetic variation in the ethanol signaling pathway. We used this information to identify genes required for acquisition of ethanol tolerance in wild strains, including new genes and processes not previously linked to ethanol tolerance, and four genes that increase ethanol tolerance when overexpressed. Our approach shows that comparative genomics across natural isolates can quickly identify genes for industrial engineering while expanding our understanding of natural diversity.

Toxicogenomics using yeast DNA microarrays

Development of genomics and bioinformatics enable us to analyze the global gene expression profiles of cells by DNA microarray. Changes in gene expression patterns indicate changes in its physiological conditions. Following the exposure of an organism or cell to toxic chemicals or other environmental stresses, the global genetic responses can be expeditiously and easily analyzed. Baker's yeast, Saccharomyces cerevisiae, is one of the most studied and useful model eukaryotes. The biggest advantage of yeast genomics is the available functional information for each gene and a considerable number of data are accumulating in the field of toxicity assessment using yeast DNA microarray. In this review, we discuss the toxicogenomics of metal ions, alcohols and aldehydes, and other chemicals.

Transient genotype-by-environment interactions following environmental shock provide a source of expression variation for essential genes

Understanding complex genotype-by-environment interactions (GEIs) is crucial for understanding phenotypic variation. An important factor often overlooked in GEI studies is time. We measured the contribution of GEIs to expression variation in four nonlaboratory Saccharomyces cerevisiae strains responding dynamically to a 25 degrees -37 degrees heat shock. GEI was a major force explaining expression variation, affecting 55% of the genes analyzed. Importantly, almost half of these expression patterns showed GEI influence only during the transition between environments, but not in acclimated cells. This class reveals a genotype-by-environment-by-time interaction that affected expression of a large fraction of yeast genes. Strikingly, although transcripts subject to persistent GEI effects were enriched for nonessential genes with upstream TATA elements, those displaying transient GEIs were enriched for essential genes regardless of TATA regulation. Genes subject to persistent GEI influences showed relaxed constraint on acclimated gene expression compared to the average yeast gene, whereas genes restricted to transient GEIs did not. We propose that transient GEI during the transition between environments provides a previously unappreciated source of expression variation, particularly for essential genes.

The histone deacetylase Rpd3p is required for transient changes in genomic expression in response to stress

Chromatin-immunoprecipitation and computational analysis implicate Rpd3p as an important co-factor in the network of genes regulating the yeast environmental stress response.

Background

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.

Results

We found that the Rpd3-Large complex is required for proper expression of both induced and repressed Environmental Stress Response genes under multiple stress conditions. Cells lacking RPD3 or the Rpd3-Large subunit PHO23 had a major defect in Environmental Stress Response initiation, particularly during the transient phase of expression immediately after stress exposure. Chromatin-immunoprecipitation showed a direct role for Rpd3-Large at representative genes; however, there were different effects on nucleosome occupancy and histone deacetylation at different promoters. Computational analysis implicated regulators that may act with Rpd3p at Environmental Stress Response genes. We provide genetic and biochemical evidence that Rpd3p is required for binding and action of the stress-activated transcription factor Msn2p, although the contribution of these factors differs for different genes.

Conclusions

Our results implicate Rpd3p as an important co-factor in the Environmental Stress Response regulatory network, and suggest the importance of histone modification in producing transient changes in gene expression triggered by stress.

Variations in stress sensitivity and genomic expression in diverse S. cerevisiae isolates

Interactions between an organism and its environment can significantly influence phenotypic evolution. A first step toward understanding this process is to characterize phenotypic diversity within and between populations. We explored the phenotypic variation in stress sensitivity and genomic expression in a large panel of Saccharomyces strains collected from diverse environments. We measured the sensitivity of 52 strains to 14 environmental conditions, compared genomic expression in 18 strains, and identified gene copy-number variations in six of these isolates. Our results demonstrate a large degree of phenotypic variation in stress sensitivity and gene expression. Analysis of these datasets reveals relationships between strains from similar niches, suggests common and unique features of yeast habitats, and implicates genes whose variable expression is linked to stress resistance. Using a simple metric to suggest cases of selection, we found that strains collected from oak exudates are phenotypically more similar than expected based on their genetic diversity, while sake and vineyard isolates display more diverse phenotypes than expected under a neutral model. We also show that the laboratory strain S288c is phenotypically distinct from all of the other strains studied here, in terms of stress sensitivity, gene expression, Ty copy number, mitochondrial content, and gene-dosage control. These results highlight the value of understanding the genetic basis of phenotypic variation and raise caution about using laboratory strains for comparative genomics.

Much attention has been given to the ways in which organisms evolve new phenotypes and the influence of the environment on this process. A major focus of study is defining the genetic basis for phenotypes important for organismal fitness. As a first step toward this goal, we surveyed phenotypic variation in diverse yeast strains collected from different environments by characterizing variations in stress resistance and genomic expression. We uncovered many phenotypic differences across yeast strains, both in stress tolerance and gene expression. The similarities and differences of the strains analyzed uncovered phenotypes shared by strains that live in similar environments, suggesting common features of yeast niches as well as mechanisms that different strains use to thrive in those conditions. We provide evidence that some characteristics of strains isolated from oak tree soil have been selected for, perhaps because of the shared selective pressures imposed by their environment. One theme emerging from our studies is that the laboratory strain of yeast, long used as a model for yeast physiology and basic biology, is aberrant compared to all other strains. This result raises caution about making general conclusions about yeast biology based on a single strain with a specific genetic makeup.

Stress-activated genomic expression changes serve a preparative role for impending stress in yeast

Yeast cells respond to stress by mediating condition-specific gene expression changes and by mounting a common response to many stresses, called the environmental stress response (ESR). Giaever et al. previously revealed poor correlation between genes whose expression changes in response to acute stress and genes required to survive that stress, raising question about the role of stress-activated gene expression. Here we show that gene expression changes triggered by a single dose of stress are not required to survive that stimulus but rather serve a protective role against future stress. We characterized the increased resistance to severe stress in yeast preexposed to mild stress. This acquired stress resistance is dependent on protein synthesis during mild-stress treatment and requires the "general-stress" transcription factors Msn2p and/or Msn4p that regulate induction of many ESR genes. However, neither protein synthesis nor Msn2/4p is required for basal tolerance of a single dose of stress, despite the substantial expression changes triggered by each condition. Using microarrays, we show that Msn2p and Msn4p play nonredundant and condition-specific roles in gene-expression regulation, arguing against a generic general-stress function. This work highlights the importance of condition-specific responses in acquired stress resistance and provides new insights into the role of the ESR.

Antioxidant protection of resveratrol and catechin in Saccharomyces cerevisiae

Moderate consumption of red wine reduces the risk of heart disease and extends lifespan, but the relative contribution of wine polyphenols to these effects is unclear. In this work, the capacity of resveratrol and catechin to protect the eukaryotic microorganism Saccharomyces cerevisiae against oxidative stress caused by different agents, hydrogen peroxide, carbon tetrachloride, and cadmium, was evaluated. Under all stress conditions, both polyphenols increased tolerance, although their protection was more evident under peroxide exposure. By using mutant strains deficient in specific antioxidant defense systems (superoxide dismutases, catalase, or glutathione), it was observed that increased H2O2 tolerance produced by both polyphenols was associated with catalase, as well as the rise in survival rates caused by resveratrol under CCl4. The acquisition of tolerance was correlated with a reduction in lipid peroxidation, indicating that the antioxidant property of resveratrol and catechin involves protection against membrane oxidation.

Transcription-associated mutagenesis in yeast is directly proportional to the level of gene expression and influenced by the direction of DNA replication

A high level of transcription has been associated with elevated spontaneous mutation and recombination rates in eukaryotic organisms. To determine whether the transcription level is directly correlated with the degree of genomic instability, we have developed a tetracycline-regulated LYS2 reporter system to modulate the transcription level over a broad range in Saccharomyces cerevisiae. We find that spontaneous mutation rate is directly proportional to the transcription level, suggesting that movement of RNA polymerase through the target initiates a mutagenic process(es). Using this system, we also investigated two hypotheses that have been proposed to explain transcription-associated mutagenesis (TAM): (1) transcription impairs replication fork progression in a directional manner and (2) DNA lesions accumulate under high-transcription conditions. The effect of replication fork progression was probed by comparing the mutational rates and spectra in yeast strains with the reporter gene placed in two different orientations near a well-characterized replication origin. The effect of endogenous DNA damage accumulation was investigated by studying TAM in strains defective in nucleotide excision repair or in lesion bypass by the translesion polymerase Polzeta. Our results suggest that both replication orientation and endogenous lesion accumulation play significant roles in TAM, particularly in terms of mutation spectra.

Sampling Saccharomyces cerevisiae cells by rapid filtration improves the yield of mRNAs

To optimize the recovery of mRNAs extracted from yeast, different methods for sampling the yeast cells have been compared. For Saccharomyces cerevisiae strains growing on gluconeogenic carbon sources (derepressed cells) rapid filtration allowed much higher yields (3-10 fold) than centrifugation at room temperature or at 4 degrees C. Recovery of total RNA was similar with the different procedures. For S. cerevisiae growing on glucose, filtration caused a 2-4 fold improvement on the mRNA yields obtained from cells sampled by centrifugation. It was also observed that, when derepressed cells of S. cerevisiae W303-1A were collected by filtration and flash-frozen, part of the 25S and 18S rRNAs (up to 50%) was recovered in an unprocessed 32S or 33S form.

yMGV: a cross-species expression data mining tool

The yeast Microarray Global Viewer (yMGV @ http://transcriptome.ens.fr/ymgv) was created 3 years ago as a database that houses a collection of Saccharomyces cerevisiae and Schizosaccharo myces pombe microarray data sets published in 82 different articles. yMGV couples data mining tools with a user-friendly web interface so that, with a few mouse clicks, one can identify the conditions that affect the expression of a gene or list of genes regulated in a set of experiments. One of the major new features we present here is a set of tools that allows for inter-organism comparisons. This should enable the fission yeast community to take advantage of the large amount of available information on budding yeast transcriptome. New tools and ongoing developments are also presented here.

Standardization of protocols in cDNA microarray analysis

Systematic variations can occur at various steps of a cDNA microarray experiment and affect the measurement of gene expression levels. Accepted standards integrated into every cDNA microarray analysis can assess these variabilities and aid the interpretation of cDNA microarray experiments from different sources. A universally applicable approach to evaluate parameters such as input and output ratios, signal linearity, hybridization specificity and consistency across an array, as well as normalization strategies, is the utilization of exogenous control genes as spike-in and negative controls. We suggest that the use of such control sets, together with a sufficient number of experimental repeats, in-depth statistical analysis and thorough data validation should be made mandatory for the publication of cDNA microarray data.