Global protein expression profiling of budding yeast in response to DNA damage

Saccharomyces cerevisiae employs complex regulation strategies to tolerate low pH stress during ethanol production

Background

Industrial bioethanol production may involve a low pH environment caused by inorganic acids, improving the tolerance of Saccharomyces cerevisiae to a low pH environment is of industrial importance to increase ethanol yield, control bacterial contamination, and reduce production cost. In our previous study, acid tolerance of a diploid industrial Saccharomyces cerevisiae strain KF-7 was chronically acclimatized by continuous ethanol fermentation under gradually increasing low-pH stress conditions. Two haploid strains B3 and C3 having excellent low pH tolerance were derived through the sporulation of an isolated mutant. Diploid strain BC3 was obtained by mating these two haploids. In this study, B3, C3, BC3, and the original strain KF-7 were subjected to comparison transcriptome analysis to investigate the molecular mechanism of the enhanced phenotype.

Result

The comparison transcriptome analysis results suggested that the upregulated vitamin B1 and B6 biosynthesis contributed to the low pH tolerance. Amino acid metabolism, DNA repairment, and general stress response might also alleviate low pH stress.

Conclusion

Saccharomyces cerevisiae seems to employ complex regulation strategies to tolerate low pH during ethanol production. The findings provide guides for the construction of low pH-tolerant industrial strains that can be used in industrial fermentation processes.

Quantitative Characterisation of Low Abundant Yeast Mitochondrial Proteins Reveals Compensation for Haplo-Insufficiency in Different Environments

The quantification of low abundant membrane-binding proteins such as transcriptional factors and chaperones has proven difficult, even with the most sophisticated analytical technologies. Here, we exploit and optimise the non-invasive Fluorescence Correlation Spectroscopy (FCS) for the quantitation of low abundance proteins, and as proof of principle, we choose two interacting proteins involved in the fission of mitochondria in yeast, Fis1p and Mdv1p. In Saccharomyces cerevisiae, the recruitment of Fis1p and Mdv1p to mitochondria is essential for the scission of the organelles and the retention of functional mitochondrial structures in the cell. We use FCS in single GFP-labelled live yeast cells to quantify the protein abundance in homozygote and heterozygote cells and to investigate the impact of the environments on protein copy number, bound/unbound protein state and mobility kinetics. Both proteins were observed to localise predominantly at mitochondrial structures, with the Mdv1p bound state increasing significantly in a strictly respiratory environment. Moreover, a compensatory mechanism that controls Fis1p abundance upon deletion of one allele was observed in Fis1p but not in Mdv1p, suggesting differential regulation of Fis1p and Mdv1p protein expression.

Transcriptional Profiling of the Candida albicans Response to the DNA Damage Agent Methyl Methanesulfonate

The infection of a mammalian host by the pathogenic fungus Candida albicans involves fungal resistance to reactive oxygen species (ROS)—induced DNA damage stress generated by the defending macrophages or neutrophils. Thus, the DNA damage response in C. albicans may contribute to its pathogenicity. Uncovering the transcriptional changes triggered by the DNA damage—inducing agent MMS in many model organisms has enhanced the understanding of their DNA damage response processes. However, the transcriptional regulation triggered by MMS remains unclear in C. albicans. Here, we explored the global transcription profile in response to MMS in C. albicans and identified 306 defined genes whose transcription was significantly affected by MMS. Only a few MMS-responsive genes, such as MGT1, DDR48, MAG1, and RAD7, showed potential roles in DNA repair. GO term analysis revealed that a large number of induced genes were involved in antioxidation responses, and some downregulated genes were involved in nucleosome packing and IMP biosynthesis. Nevertheless, phenotypic assays revealed that MMS-induced antioxidation gene CAP1 and glutathione metabolism genes GST2 and GST3 showed no direct roles in MMS resistance. Furthermore, the altered transcription of several MMS—responsive genes exhibited RAD53—related regulation. Intriguingly, the transcription profile in response to MMS in C. albicans shared a limited similarity with the pattern in S. cerevisiae, including COX17, PRI2, and MGT1. Overall, C. albicans cells exhibit global transcriptional changes to the DNA damage agent MMS; these findings improve our understanding of this pathogen’s DNA damage response pathways.

Telomere Interacting Proteins and TERRA Regulation

Telomere shortening rates inversely correlate with life expectancy and hence it is critical to understand how telomere shortening is regulated. Recently, the telomeric non-coding RNA, TERRA has been implicated in the regulation of replicative senescence. To better understand how TERRA is regulated we employed a proteomics approach to look for potential RNA regulators that associate with telomeric sequences. Based on the results, we have identified proteins that may regulate TERRA in both a positive and negative manner, depending on the state of the telomere. In this mini-review, we discuss and speculate about these data to expand our understanding of TERRA and telomere interactors with respect to telomere shortening dynamics.

Molecular modeling, docking and dynamics analysis of lipid droplet associated enzyme Ypr147cp from Saccharomyces cerevisiae

Ypr147cp of Saccharomyces cerevisiae was localized to lipid droplets. The recombinant Ypr147cp showed both triacylglycerol lipase and ester hydrolase activities. Knock out of YPR147C led to accumulation of TAG in ypr147cΔ when compared to wild type (WT). Transmission electron microscopic analysis of ypr147cΔ cells show increased lipid bodies. Moreover, the lipid profiling confirmed the accumulation of fatty acids derived from neutral and phospholipids in ypr147cΔ cells. Sequence analysis of Ypr147cp show the presence of an a/b hydrolase domain with the conserved GXSXG lipase motif. The YPR147c homology model was built and the modeled protein was analysed using RMSD and root mean square fluctuation (RMSF) for a 100 ns simulation trajectory. Docking the acetate, butyrate and palmitate ligands with the model confirmed covalent binding of ligands with the Ser207 of the GXSXG motif. Thus, Ypr147cp is a lipid droplet associated triacylglycerol lipase having short chain ester hydrolyzing capacity.

Incorporation of a unified protein abundance dataset into the Saccharomyces genome database

The identification and accurate quantitation of protein abundance has been a major objective of proteomics research. Abundance studies have the potential to provide users with data that can be used to gain a deeper understanding of protein function and regulation and can also help identify cellular pathways and modules that operate under various environmental stress conditions. One of the central missions of the Saccharomyces Genome Database (SGD; https://www.yeastgenome.org) is to work with researchers to identify and incorporate datasets of interest to the wider scientific community, thereby enabling hypothesis-driven research. A large number of studies have detailed efforts to generate proteome-wide abundance data, but deeper analyses of these data have been hampered by the inability to compare results between studies. Recently, a unified protein abundance dataset was generated through the evaluation of more than 20 abundance datasets, which were normalized and converted to common measurement units, in this case molecules per cell. We have incorporated these normalized protein abundance data and associated metadata into the SGD database, as well as the SGD YeastMine data warehouse, resulting in the addition of 56 487 values for untreated cells grown in either rich or defined media and 28 335 values for cells treated with environmental stressors. Abundance data for protein-coding genes are displayed in a sortable, filterable table on Protein pages, available through Locus Summary pages. A median abundance value was incorporated, and a median absolute deviation was calculated for each protein-coding gene and incorporated into SGD. These values are displayed in the Protein section of the Locus Summary page. The inclusion of these data has enhanced the quality and quantity of protein experimental information presented at SGD and provides opportunities for researchers to access and utilize the data to further their research.

Definition of the Minimal Contents for the Molecular Simulation of the Yeast Cytoplasm

The cytoplasm is a densely packed environment filled with macromolecules with hindered diffusion. Molecular simulation of the diffusion of biomolecules under such macromolecular crowding conditions requires the definition of a simulation cell with a cytoplasmic-like composition. This has been previously done for prokaryote cells (E. coli) but not for eukaryote cells such as yeast as a model organism. Yeast proteomics datasets vary widely in terms of cell growth conditions, the technique used to determine protein composition, the reported relative abundance of proteins, and the units in which abundances are reported. We determined that the gene ontology profiles of the most abundant proteins across these datasets are similar, but their abundances vary greatly. To overcome this problem, we chose five mass spectrometry proteomics datasets that fulfilled the following criteria: high internal consistency, consistency with published experimental data, and freedom from GFP-tagging artifacts. Using these datasets, the contents of a simulation cell containing a single 80S ribosome were defined, such that the macromolecular density and the mass ratio of ribosomal-to-cytoplasmic proteins were consistent with experiment and chosen datasets. Finally, multiple tRNAs were added, consistent with their experimentally-determined number in the yeast cell. The resulting composition can be readily used in molecular simulations representative of yeast cytoplasmic macromolecular crowding conditions to characterize a variety of phenomena, such as protein diffusion, protein-protein interactions and biological processes such as protein translation.

Mutations in the PCNA DNA Polymerase Clamp of Saccharomyces cerevisiae Reveal Complexities of the Cell Cycle and Ploidy on Heterochromatin Assembly

In Saccharomyces cerevisiae, transcriptional silencing at HML and HMR maintains mating-type identity. The repressive chromatin structure at these loci is replicated every cell cycle and must be re-established quickly to prevent transcription of the genes at these loci. Mutations in a component of the replisome, the proliferating cell nuclear antigen (PCNA), encoded by POL30, cause a loss of transcriptional silencing at HMR We used an assay that captures transient losses of silencing at HML and HMR to perform extended genetic analyses of the pol30-6, pol30-8, and pol30-79 alleles. All three alleles destabilized silencing only transiently and only in cycling cells. Whereas pol30-8 caused loss of silencing by disrupting the function of Chromatin Assembly Factor 1, pol30-6 and pol30-79 acted through a separate genetic pathway, but one still dependent on histone chaperones. Surprisingly, the silencing-loss phenotypes of pol30-6 and pol30-79 depended on ploidy, but not on POL30 dosage or mating-type identity. Separately from silencing loss, the pol30-6 and pol30-79 alleles also displayed high levels of mitotic recombination in diploids. These results established that histone trafficking involving PCNA at replication forks is crucial to the maintenance of chromatin state and genome stability during DNA replication. They also raised the possibility that increased ploidy may protect chromatin states when the replisome is perturbed.

Lsm12 Mediates Deubiquitination of DNA Polymerase η To Help Saccharomyces cerevisiae Resist Oxidative Stress

In Saccharomyces cerevisiae, the Y family DNA polymerase η (Polη) regulates genome stability in response to different forms of environmental stress by translesion DNA synthesis. To elucidate the role of Polη in oxidative stress-induced DNA damage, we deleted or overexpressed the corresponding gene RAD30 and used transcriptome analysis to screen the potential genes associated with RAD30 to respond to DNA damage. Under 2 mM H₂O₂ treatment, the deletion of RAD30 resulted in a 2.2-fold decrease in survival and a 2.8-fold increase in DNA damage, whereas overexpression of RAD30 increased survival and decreased DNA damage by 1.2- and 1.4-fold, respectively, compared with the wild-type strain. Transcriptome and phenotypic analyses identified Lsm12 as a main factor involved in oxidative stress-induced DNA damage. Deleting LSM12 caused growth defects, while its overexpression enhanced cell growth under 2 mM H₂O₂ treatment. This effect was due to the physical interaction of Lsm12 with the UBZ domain of Polη to enhance Polη deubiquitination through Ubp3 and consequently promote Polη recruitment. Overall, these findings demonstrate that Lsm12 is a novel regulator mediating Polη deubiquitination to promote its recruitment under oxidative stress. Furthermore, this study provides a potential strategy to maintain the genome stability of industrial strains during fermentation.IMPORTANCE Polη was shown to be critical for cell growth in the yeast Saccharomyces cerevisiae, and deletion of its corresponding gene RAD30 caused a severe growth defect under exposure to oxidative stress with 2 mM H₂O₂ Furthermore, we found that Lsm12 physically interacts with Polη and promotes Polη deubiquitination and recruitment. Overall, these findings indicate Lsm12 is a novel regulator mediating Polη deubiquitination that regulates its recruitment in response to DNA damage induced by oxidative stress.

Intracellular mechanism by which genotoxic stress activates yeast SAPK Mpk1

Stress-activated MAP kinases (SAPKs) respond to a wide variety of stressors. In most cases, the pathways through which specific stress signals are transmitted to the SAPKs are not known. The Saccharomyces cerevisiae SAPK Mpk1 (Slt2) is a well-characterized component of the cell-wall integrity (CWI) signaling pathway, which responds to physical and chemical challenges to the cell wall. However, Mpk1 is also activated in response to genotoxic stress through an unknown pathway. We show that, in contrast to cell-wall stress, the pathway for Mpk1 activation by genotoxic stress does not involve the stimulation of the MAP kinase kinases (MEKs) that function immediately upstream of Mpk1. Instead, DNA damage activates Mpk1 through induction of proteasomal degradation of Msg5, the dual-specificity protein phosphatase principally responsible for maintaining Mpk1 in a low-activity state in the absence of stress. Blocking Msg5 degradation in response to genotoxic stress prevented Mpk1 activation. This work raises the possibility that other Mpk1-activating stressors act intracellularly at different points along the canonical Mpk1 activation pathway.

Unification of Protein Abundance Datasets Yields a Quantitative Saccharomyces cerevisiae Proteome

Protein activity is the ultimate arbiter of function in most cellular pathways, and protein concentration is fundamentally connected to protein action. While the proteome of yeast has been subjected to the most comprehensive analysis of any eukaryote, existing datasets are difficult to compare, and there is no consensus abundance value for each protein. We evaluated 21 quantitative analyses of the S. cerevisiae proteome, normalizing and converting all measurements of protein abundance into the intuitive measurement of absolute molecules per cell. We estimate the cellular abundance of 92% of the proteins in the yeast proteome and assess the variation in each abundance measurement. Using our protein abundance dataset, we find that a global response to diverse environmental stresses is not detected at the level of protein abundance, we find that protein tags have only a modest effect on protein abundance, and we identify proteins that are differentially regulated at the mRNA abundance, mRNA translation, and protein abundance levels.

Cell cycle-dependent positive and negative functions of Fun30 chromatin remodeler in DNA damage response

The evolutionally conserved Fun30 chromatin remodeler in Saccharomyces cerevisiae has been shown to contribute to cellular resistance to genotoxic stress inflicted by camptothecin (CPT), methyl methanesulfonate (MMS) and hydroxyurea (HU). Fun30 aids in extensive DNA resection of DNA double stranded break (DSB) ends, which is thought to underlie its role in CPT-resistance. How Fun30 promotes MMS- or HU-resistance has not been resolved. Interestingly, we have recently found Fun30 to also play a negative role in cellular tolerance to MMS and HU in the absence of the Rad5-dependent DNA damage tolerance pathway. In this report, we show that Fun30 acts to down regulate Rad9-dependent DNA damage checkpoint triggered by CPT or MMS, but does not affect Rad9-independent intra-S phase replication checkpoint induced by MMS or HU. These results support the notion that Fun30 contributes to cellular response to DSBs by preventing excessive DNA damage checkpoint activation in addition to its role in facilitating DNA end resection. On the other hand, we present evidence suggesting that Fun30's negative function in MMS- and HU-tolerance in the absence of Rad5 is not related to its regulation of checkpoint activity. Moreover, we find Fun30 to be cell cycle regulated with its abundance peaking in G2/M phase of the cell cycle. Importantly, we demonstrate that artificially restricting Fun30 expression to G2/M does not affect its positive or negative function in genotoxin-resistance, but confining Fun30 to S phase abolishes its functions. These results indicate that both positive and negative functions of Fun30 in DNA damage response occur mainly in G2/M phase.

Orphan proteins of unknown function in the mitochondrial intermembrane space proteome: New pathways and metabolic cross-talk

The mitochondrial intermembrane space (IMS) is involved in protein transport, lipid homeostasis and metal ion exchange, while further acting in signalling pathways such as apoptosis. Regulation of these processes involves protein modifications, as well as stress-induced import or release of proteins and other signalling molecules. Even though the IMS is the smallest sub-compartment of mitochondria, its redox state seems to be tightly regulated. However, the way in which this compartment participates in the cross-talk between the multiple organelles and the cytosol is far from understood. Here we focus on newly identified IMS proteins that may represent future challenges in mitochondrial research. We present an overview of the import pathways, the recently discovered new components of the IMS proteome and how these relate to key aspects of cell signalling and progress made in stem cell and cancer research.

•

A brief overview of the classic mitochondrial import pathways is featured

•

Recent studies assigning a number of new proteins to the mitochondrial IMS are discussed

•

Analysis of the expanded IMS proteomes can provide insights into organelle cross-talk and signalling pathways

Combined Gene Expression and RNAi Screening to Identify Alkylation Damage Survival Pathways from Fly to Human

Alkylating agents are a key component of cancer chemotherapy. Several cellular mechanisms are known to be important for its survival, particularly DNA repair and xenobiotic detoxification, yet genomic screens indicate that additional cellular components may be involved. Elucidating these components has value in either identifying key processes that can be modulated to improve chemotherapeutic efficacy or may be altered in some cancers to confer chemoresistance. We therefore set out to reevaluate our prior Drosophila RNAi screening data by comparison to gene expression arrays in order to determine if we could identify any novel processes in alkylation damage survival. We noted a consistent conservation of alkylation survival pathways across platforms and species when the analysis was conducted on a pathway/process level rather than at an individual gene level. Better results were obtained when combining gene lists from two datasets (RNAi screen plus microarray) prior to analysis. In addition to previously identified DNA damage responses (p53 signaling and Nucleotide Excision Repair), DNA-mRNA-protein metabolism (transcription/translation) and proteasome machinery, we also noted a highly conserved cross-species requirement for NRF2, glutathione (GSH)-mediated drug detoxification and Endoplasmic Reticulum stress (ER stress)/Unfolded Protein Responses (UPR) in cells exposed to alkylation. The requirement for GSH, NRF2 and UPR in alkylation survival was validated by metabolomics, protein studies and functional cell assays. From this we conclude that RNAi/gene expression fusion is a valid strategy to rapidly identify key processes that may be extendable to other contexts beyond damage survival.

Evidence that COG0325 proteins are involved in PLP homeostasis

Pyridoxal 5'-phosphate (PLP) is an essential cofactor for nearly 60 Escherichia coli enzymes but is a highly reactive molecule that is toxic in its free form. How PLP levels are regulated and how PLP is delivered to target enzymes are still open questions. The COG0325 protein family belongs to the fold-type III class of PLP enzymes and binds PLP but has no known biochemical activity although it occurs in all kingdoms of life. Various pleiotropic phenotypes of the E. coli COG0325 (yggS) mutant have been reported, some of which were reproduced and extended in this study. Comparative genomic, genetic and metabolic analyses suggest that these phenotypes reflect an imbalance in PLP homeostasis. The E. coli yggS mutant accumulates the PLP precursor pyridoxine 5'-phosphate (PNP) and is sensitive to an excess of pyridoxine but not of pyridoxal. The pyridoxine toxicity phenotype is complemented by the expression of eukaryotic yggS orthologs. It is also suppressed by the presence of amino acids, specifically isoleucine, threonine and leucine, suggesting the PLP-dependent enzyme transaminase B (IlvE) is affected. These genetic results lay a foundation for future biochemical studies of the role of COG0325 proteins in PLP homeostasis.

Identification and functional evaluation of the reductases and dehydrogenases from Saccharomyces cerevisiae involved in vanillin resistance

Background

Vanillin, a type of phenolic released during the pre-treatment of lignocellulosic materials, is toxic to microorganisms and therefore its presence inhibits the fermentation. The vanillin can be reduced to vanillyl alcohol, which is much less toxic, by the ethanol producer Saccharomyces cerevisiae. The reducing capacity of S. cerevisiae and its vanillin resistance are strongly correlated. However, the specific enzymes and their contribution to the vanillin reduction are not extensively studied. In our previous work, an evolved vanillin-resistant strain showed an increased vanillin reduction capacity compared with its parent strain. The transcriptome analysis suggested the reductases and dehydrogenases of this vanillin resistant strain were up-regulated. Using this as a starting point, 11 significantly regulated reductases and dehydrogenases were selected in the present work for further study. The roles of these reductases and dehydrogenases in the vanillin tolerance and detoxification abilities of S. cerevisiae are described.

Results

Among the candidate genes, the overexpression of the alcohol dehydrogenase gene ADH6, acetaldehyde dehydrogenase gene ALD6, glucose-6-phosphate 1-dehydrogenase gene ZWF1, NADH-dependent aldehyde reductase gene YNL134C, and aldo-keto reductase gene YJR096W increased 177, 25, 6, 15, and 18 % of the strain μ_max in the medium containing 1 g L⁻¹ vanillin. The in vitro detected vanillin reductase activities of strain overexpressing ADH6, YNL134C and YJR096W were notably higher than control. The vanillin specific reduction rate increased by 8 times in ADH6 overexpressed strain but not in YNL134C and YJR096W overexpressed strain. This suggested that the enzymes encoded by YNL134C and YJR096W might prefer other substrate and/or could not show their effects on vanillin on the high background of Adh6p in vivo. Overexpressing ALD6 and ZWF1 mainly increased the [NADPH]/[NADP⁺] and [GSH]/[GSSG] ratios but not the vanillin reductase activities. Their contribution to strain growth and vanillin reduction were balancing the redox state of strain when vanillin was presented.

Conclusions

Beside the reported Adh6p, the enzymes encoded by YNL134C and YJR096W were proved to have vanillin reduction activity in present study. While ALD6 and ZWF1 did not directly reduce vanillin to vanillyl alcohol, their contribution to vanillin resistance primarily depended on the enhancement of the reducing equivalent supply.

Electronic supplementary material

The online version of this article (doi:10.1186/s12896-016-0264-y) contains supplementary material, which is available to authorized users.

High-throughput fluorescence microscopic analysis of protein abundance and localization in budding yeast

Proteins directly carry out and regulate cellular functions. As a result, changes in protein levels within a cell directly influence cellular processes. Similarly, it is intuitive that the intracellular localization of proteins is a key component of their functionality. Optimal activity is achieved by a combination of protein concentration, co-compartmentalization with substrates, co-factors and regulators and sequestration from deleterious locales. The proteome within a cell is highly dynamic and changes in response to different environmental conditions. High-throughput microscopic analysis in the budding yeast Saccharomyces cerevisiae has afforded proteome-wide views of protein organization in living cells, and of how protein abundance and location is regulated and remodeled in response to stress.

Integrated RNA- and protein profiling of fermentation and respiration in diploid budding yeast provides insight into nutrient control of cell growth and development

Diploid budding yeast undergoes rapid mitosis when it ferments glucose, and in the presence of a non-fermentable carbon source and the absence of a nitrogen source it triggers sporulation. Rich medium with acetate is a commonly used pre-sporulation medium, but our understanding of the molecular events underlying the acetate-driven transition from mitosis to meiosis is still incomplete. We identified 263 proteins for which mRNA and protein synthesis are linked or uncoupled in fermenting and respiring cells. Using motif predictions, interaction data and RNA profiling we find among them 28 likely targets for Ume6, a subunit of the conserved Rpd3/Sin3 histone deacetylase-complex regulating genes involved in metabolism, stress response and meiosis. Finally, we identify 14 genes for which both RNA and proteins are detected exclusively in respiring cells but not in fermenting cells in our sample set, including CSM4, SPR1, SPS4 and RIM4, which were thought to be meiosis-specific. Our work reveals intertwined transcriptional and post-transcriptional control mechanisms acting when a MATa/α strain responds to nutritional signals, and provides molecular clues how the carbon source primes yeast cells for entering meiosis.

BIOLOGICAL SIGNIFICANCE

Our integrated genomics study provides insight into the interplay between the transcriptome and the proteome in diploid yeast cells undergoing vegetative growth in the presence of glucose (fermentation) or acetate (respiration). Furthermore, it reveals novel target genes involved in these processes for Ume6, the DNA binding subunit of the conserved histone deacetylase Rpd3 and the co-repressor Sin3. We have combined data from an RNA profiling experiment using tiling arrays that cover the entire yeast genome, and a large-scale protein detection analysis based on mass spectrometry in diploid MATa/α cells. This distinguishes our study from most others in the field-which investigate haploid yeast strains-because only diploid cells can undergo meiotic development in the simultaneous absence of a non-fermentable carbon source and nitrogen. Indeed, we report molecular clues how respiration of acetate might prime diploid cells for efficient spore formation, a phenomenon that is well known but poorly understood.

Mass spectrometry-based quantification of the cellular response to methyl methanesulfonate treatment in human cells

Faithful transmission of genetic material is essential for cell viability and organism health. The occurrence of DNA damage, due to either spontaneous events or environmental agents, threatens the integrity of the genome. The consequences of these insults, if allowed to perpetuate and accumulate over time, are mutations that can lead to the development of diseases such as cancer. Alkylation is a relevant DNA lesion produced endogenously as well as by exogenous agents including certain chemotherapeutics. We sought to better understand the cellular response to this form of DNA damage using mass spectrometry-based proteomics. For this purpose, we performed sub-cellular fractionation to monitor the effect of methyl methanesulfonate (MMS) treatment on protein localization to chromatin. The levels of over 500 proteins were increased in the chromatin-enriched nuclear lysate including histone chaperones. Levels of ubiquitin and subunits of the proteasome were also increased within this fraction, suggesting that ubiquitin-mediated degradation by the proteasome has an important role in the chromatin response to MMS treatment. Finally, the levels of some proteins were decreased within the chromatin-enriched lysate including components of the nuclear pore complex. Our spatial proteomics data demonstrate that many proteins that influence chromatin organization are regulated in response to MMS treatment, presumably to open the DNA to allow access by other DNA damage response proteins. To gain further insight into the cellular response to MMS-induced DNA damage, we also performed phosphorylation enrichment on total cell lysates to identify proteins regulated via post-translational modification. Phosphoproteomic analysis demonstrated that many nuclear phosphorylation events were decreased in response to MMS treatment. This reflected changes in protein kinase and/or phosphatase activity in response to DNA damage rather than changes in total protein abundance. Using these two mass spectrometry-based approaches, we have identified a novel set of MMS-responsive proteins that will expand our understanding of DNA damage signaling.

Genome-wide single-cell-level screen for protein abundance and localization changes in response to DNA damage in S. cerevisiae

An effective response to DNA damaging agents involves modulating numerous facets of cellular homeostasis in addition to DNA repair and cell-cycle checkpoint pathways. Fluorescence microscopy-based imaging offers the opportunity to simultaneously interrogate changes in both protein level and subcellular localization in response to DNA damaging agents at the single-cell level. We report here results from screening the yeast Green Fluorescent Protein (GFP)-fusion library to investigate global cellular protein reorganization on exposure to the alkylating agent methyl methanesulfonate (MMS). Broad groups of induced, repressed, nucleus- and cytoplasm-enriched proteins were identified. Gene Ontology and interactome analyses revealed the underlying cellular processes. Transcription factor (TF) analysis identified principal regulators of the response, and targets of all major stress-responsive TFs were enriched amongst the induced proteins. An unexpected partitioning of biological function according to the number of TFs targeting individual genes was revealed. Finally, differential modulation of ribosomal proteins depending on methyl methanesulfonate dose was shown to correlate with cell growth and with the translocation of the Sfp1 TF. We conclude that cellular responses can navigate different routes according to the extent of damage, relying on both expression and localization changes of specific proteins.

Towards systematic discovery of signaling networks in budding yeast filamentous growth stress response using interventional phosphorylation data

Reversible phosphorylation is one of the major mechanisms of signal transduction, and signaling networks are critical regulators of cell growth and development. However, few of these networks have been delineated completely. Towards this end, quantitative phosphoproteomics is emerging as a useful tool enabling large-scale determination of relative phosphorylation levels. However, phosphoproteomics differs from classical proteomics by a more extensive sampling limitation due to the limited number of detectable sites per protein. Here, we propose a comprehensive quantitative analysis pipeline customized for phosphoproteome data from interventional experiments for identifying key proteins in specific pathways, discovering the protein-protein interactions and inferring the signaling network. We also made an effort to partially compensate for the missing value problem, a chronic issue for proteomics studies. The dataset used for this study was generated using SILAC (Stable Isotope Labeling with Amino acids in Cell culture) technique with interventional experiments (kinase-dead mutations). The major components of the pipeline include phosphopeptide meta-analysis, correlation network analysis and causal relationship discovery. We have successfully applied our pipeline to interventional experiments identifying phosphorylation events underlying the transition to a filamentous growth form in Saccharomyces cerevisiae. We identified 5 high-confidence proteins from meta-analysis, and 19 hub proteins from correlation analysis (Pbi2p and Hsp42p were identified by both analyses). All these proteins are involved in stress responses. Nine of them have direct or indirect evidence of involvement in filamentous growth. In addition, we tested four of our predicted proteins, Nth1p, Pbi2p, Pdr12p and Rcn2p, by interventional phenotypic experiments and all of them present differential invasive growth, providing prospective validation of our approach. This comprehensive pipeline presents a systematic way for discovering signaling networks using interventional phosphoproteome data and can suggest candidate proteins for further investigation. We anticipate the methodology to be applicable as well to other interventional studies via different experimental platforms.

Transcriptional regulation of fermentative and respiratory metabolism in Saccharomyces cerevisiae industrial bakers' strains

Bakers' yeast-producing companies grow cells under respiratory conditions, at a very high growth rate. Some desirable properties of bakers' yeast may be altered if fermentation rather than respiration occurs during biomass production. That is why differences in gene expression patterns that take place when industrial bakers' yeasts are grown under fermentative, rather than respiratory conditions, were examined. Macroarray analysis of V1 strain indicated changes in gene expression similar to those already described in laboratory Saccharomyces cerevisiae strains: repression of most genes related to respiration and oxidative metabolism and derepression of genes related to ribosome biogenesis and stress resistance in fermentation. Under respiratory conditions, genes related to the glyoxylate and Krebs cycles, respiration, gluconeogenesis, and energy production are activated. DOG21 strain, a partly catabolite-derepressed mutant derived from V1, displayed gene expression patterns quite similar to those of V1, although lower levels of gene expression and changes in fewer number of genes as compared to V1 were both detected in all cases. However, under fermentative conditions, DOG21 mutant significantly increased the expression of SNF1 -controlled genes and other genes involved in stress resistance, whereas the expression of the HXK2 gene, involved in catabolite repression, was considerably reduced, according to the pleiotropic stress-resistant phenotype of this mutant. These results also seemed to suggest that stress-resistant genes control desirable bakers' yeast qualities.

DNA damage during G2 phase does not affect cell cycle progression of the green alga Scenedesmus quadricauda

DNA damage is a threat to genomic integrity in all living organisms. Plants and green algae are particularly susceptible to DNA damage especially that caused by UV light, due to their light dependency for photosynthesis. For survival of a plant, and other eukaryotic cells, it is essential for an organism to continuously check the integrity of its genetic material and, when damaged, to repair it immediately. Cells therefore utilize a DNA damage response pathway that is responsible for sensing, reacting to and repairing damaged DNA. We have studied the effect of 5-fluorodeoxyuridine, zeocin, caffeine and combinations of these on the cell cycle of the green alga Scenedesmus quadricauda. The cells delayed S phase and underwent a permanent G2 phase block if DNA metabolism was affected prior to S phase; the G2 phase block imposed by zeocin was partially abolished by caffeine. No cell cycle block was observed if the treatment with zeocin occurred in G2 phase and the cells divided normally. CDKA and CDKB kinases regulate mitosis in S. quadricauda; their kinase activities were inhibited by Wee1. CDKA, CDKB protein levels were stabilized in the presence of zeocin. In contrast, the protein level of Wee1 was unaffected by DNA perturbing treatments. Wee1 therefore does not appear to be involved in the DNA damage response in S. quadricauda. Our results imply a specific reaction to DNA damage in S. quadricauda, with no cell cycle arrest, after experiencing DNA damage during G2 phase.

A genetically encoded probe for the identification of proteins that form sulfenic acid in response to H2O2 in Saccharomyces cerevisiae

It is widely known that reactive oxygen species (ROS), such as hydrogen peroxide, play important roles in cellular signaling and initiation of oxidative stress responses via thiol modifications. Identification of the targets of these modifications will provide a better understanding of the relationship between ROS and human diseases, such as cancer and atherosclerosis. Sulfenic acid is the principle product of a reaction between hydrogen peroxide and a reactive protein cysteine. This reversible post-translational modification plays an important role in enzyme active sites, signaling transduction via disulfide bond formation, as well as an intermediate to overoxidation products during oxidative stress. By re-engineering the C-terminal cysteine rich domain (cCRD) of the Yap1 transcription factor, we were able to create a genetically encoded probe for the general detection and identification of proteins that form sulfenic acid in vivo. The Yap1-cCRD probe has been used previously in the identification of proteins that form sulfenic acid in Escherichia coli. Here we demonstrate the successful use of the Yap1-cCRD probe in the identification of proteins that form sulfenic acid in response to hydrogen peroxide in Saccharomyces cerevisiae.

Phenotypes associated with Saccharomyces cerevisiae Hug1 protein, a putative negative regulator of dNTP Levels, reveal similarities and differences with sequence-related Dif1

Saccharomyces cerevisiae Hugl is a small protein of unknown function that is highly inducible following replication stress and DNA damage. Its deletion suppresses the lethality of deletion of checkpoint kinase Mecl. Although DNA damage responses were largely normal in the HUG1 deletion mutant, we found enhanced resistance towards heat in logarithmic phase. In response to simultaneous carbon and replication stress, overall growth delay and less pseudohyphal filament formation were evident. These novel phenotypes are shared with deletion mutants of the negative regulators of ribonucleotide reductase, Difl and Smll. Microarray analysis showed the influence of Hugl on the expression of a large number of transcripts, including stress-related transcripts. Elevated dNTP levels in hugl Δ cells may result in a stress response reflected by the observed phenotypes and transcript profiles. However, in contrast to a deletion of structurally related Difl, Rnr2-Rnr4 subcellular localization is not grossly altered in a Hugl deletion mutant. Thus, although Hugl appears to be derived from the Rnr2-Rnr4 binding region of Difl, its mechanism of action must be independent of determining the localization of Rnr2-Rnr4.

Autophagy-dependent regulation of the DNA damage response protein ribonucleotide reductase 1

Protein synthesis and degradation are posttranscriptional pathways used by cells to regulate protein levels. We have developed a systems biology approach to identify targets of posttranscriptional regulation and we have employed this system in Saccharomyces cerevisiae to study the DNA damage response. We present evidence that 50% to 75% of the transcripts induced by alkylation damage are regulated posttranscriptionally. Significantly, we demonstrate that two transcriptionally-induced DNA damage response genes, RNR1 and RNR4, fail to show soluble protein level increases after DNA damage. To determine one of the associated mechanisms of posttranscriptional regulation, we tracked ribonucleotide reductase 1 (Rnr1) protein levels during the DNA damage response. We show that RNR1 is actively translated after damage and that a large fraction of the corresponding Rnr1 protein is packaged into a membrane-bound structure and transported to the vacuole for degradation, with these last two steps dependent on autophagy proteins. We found that inhibition of target of rapamycin (TOR) signaling and subsequent induction of autophagy promoted an increase in targeting of Rnr1 to the vacuole and a decrease in soluble Rnr1 protein levels. In addition, we demonstrate that defects in autophagy result in an increase in soluble Rnr1 protein levels and a DNA damage phenotype. Our results highlight roles for autophagy and TOR signaling in regulating a specific protein and demonstrate the importance of these pathways in optimizing the DNA damage response.

The impact of multifunctional genes on "guilt by association" analysis

Many previous studies have shown that by using variants of "guilt-by-association", gene function predictions can be made with very high statistical confidence. In these studies, it is assumed that the "associations" in the data (e.g., protein interaction partners) of a gene are necessary in establishing "guilt". In this paper we show that multifunctionality, rather than association, is a primary driver of gene function prediction. We first show that knowledge of the degree of multifunctionality alone can produce astonishingly strong performance when used as a predictor of gene function. We then demonstrate how multifunctionality is encoded in gene interaction data (such as protein interactions and coexpression networks) and how this can feed forward into gene function prediction algorithms. We find that high-quality gene function predictions can be made using data that possesses no information on which gene interacts with which. By examining a wide range of networks from mouse, human and yeast, as well as multiple prediction methods and evaluation metrics, we provide evidence that this problem is pervasive and does not reflect the failings of any particular algorithm or data type. We propose computational controls that can be used to provide more meaningful control when estimating gene function prediction performance. We suggest that this source of bias due to multifunctionality is important to control for, with widespread implications for the interpretation of genomics studies.

Genomewide expression profile analysis of the Candida glabrata Pdr1 regulon

The ABC transporters Candida glabrata Cdr1 (CgCdr1), CgPdh1, and CgSnq2 are known to mediate azole resistance in the pathogenic fungus C. glabrata. Activating mutations in CgPDR1, a zinc cluster transcription factor, result in constitutive upregulation of these ABC transporter genes but to various degrees. We examined the genomewide gene expression profiles of two matched azole-susceptible and -resistant C. glabrata clinical isolate pairs. Of the differentially expressed genes identified in the gene expression profiles for these two matched pairs, there were 28 genes commonly upregulated with CgCDR1 in both isolate sets including YOR1, LCB5, RTA1, POG1, HFD1, and several members of the FLO gene family of flocculation genes. We then sequenced CgPDR1 from each susceptible and resistant isolate and found two novel activating mutations that conferred increased resistance when they were expressed in a common background strain in which CgPDR1 had been disrupted. Microarray analysis comparing these reengineered strains to their respective parent strains identified a set of commonly differentially expressed genes, including CgCDR1, YOR1, and YIM1, as well as genes uniquely regulated by specific mutations. Our results demonstrate that while CgPdr1 activates a broad repertoire of genes, specific activating mutations result in the activation of discrete subsets of this repertoire.

A network of conserved damage survival pathways revealed by a genomic RNAi screen

Damage initiates a pleiotropic cellular response aimed at cellular survival when appropriate. To identify genes required for damage survival, we used a cell-based RNAi screen against the Drosophila genome and the alkylating agent methyl methanesulphonate (MMS). Similar studies performed in other model organisms report that damage response may involve pleiotropic cellular processes other than the central DNA repair components, yet an intuitive systems level view of the cellular components required for damage survival, their interrelationship, and contextual importance has been lacking. Further, by comparing data from different model organisms, identification of conserved and presumably core survival components should be forthcoming. We identified 307 genes, representing 13 signaling, metabolic, or enzymatic pathways, affecting cellular survival of MMS-induced damage. As expected, the majority of these pathways are involved in DNA repair; however, several pathways with more diverse biological functions were also identified, including the TOR pathway, transcription, translation, proteasome, glutathione synthesis, ATP synthesis, and Notch signaling, and these were equally important in damage survival. Comparison with genomic screen data from Saccharomyces cerevisiae revealed no overlap enrichment of individual genes between the species, but a conservation of the pathways. To demonstrate the functional conservation of pathways, five were tested in Drosophila and mouse cells, with each pathway responding to alkylation damage in both species. Using the protein interactome, a significant level of connectivity was observed between Drosophila MMS survival proteins, suggesting a higher order relationship. This connectivity was dramatically improved by incorporating the components of the 13 identified pathways within the network. Grouping proteins into "pathway nodes" qualitatively improved the interactome organization, revealing a highly organized "MMS survival network." We conclude that identification of pathways can facilitate comparative biology analysis when direct gene/orthologue comparisons fail. A biologically intuitive, highly interconnected MMS survival network was revealed after we incorporated pathway data in our interactome analysis.

Histone chaperone specificity in Rtt109 activation

Rtt109 is a histone acetyltransferase that requires a histone chaperone for the acetylation of histone 3 at lysine 56 (H3K56). Rtt109 forms a complex with the chaperone Vps75 in vivo and is implicated in DNA replication and repair. Here we show that both Rtt109 and Vps75 bind histones with high affinity, but only the complex is efficient for catalysis. The C-terminal acidic domain of Vps75 contributes to activation of Rtt109 and is necessary for in vivo functionality of Vps75, but it is not required for interaction with either Rtt109 or histones. We demonstrate that Vps75 is a structural homolog of yeast Nap1 by solving its crystal structure. Nap1 and Vps75 interact with histones and Rtt109 with comparable affinities. However, only Vps75 stimulates Rtt109 enzymatic activity. Our data highlight the functional specificity of Vps75 in Rtt109 activation.

Identification and dissection of a complex DNA repair sensitivity phenotype in Baker's yeast

Complex traits typically involve the contribution of multiple gene variants. In this study, we took advantage of a high-density genotyping analysis of the BY (S288c) and RM strains of Saccharomyces cerevisiae and of 123 derived spore progeny to identify the genetic loci that underlie a complex DNA repair sensitivity phenotype. This was accomplished by screening hybrid yeast progeny for sensitivity to a variety of DNA damaging agents. Both the BY and RM strains are resistant to the ultraviolet light-mimetic agent 4-nitroquinoline 1-oxide (4-NQO); however, hybrid progeny from a BYxRM cross displayed varying sensitivities to the drug. We mapped a major quantitative trait locus (QTL), RAD5, and identified the exact polymorphism within this locus responsible for 4-NQO sensitivity. By using a backcrossing strategy along with array-assisted bulk segregant analysis, we identified one other locus, MKT1, and a QTL on Chromosome VII that also link to the hybrid 4-NQO-sensitive phenotype but confer more minor effects. This work suggests an additive model for sensitivity to 4-NQO and provides a strategy for mapping both major and minor QTL that confer background-specific phenotypes. It also provides tools for understanding the effect of genetic background on sensitivity to genotoxic agents.

Regulation of Saccharomyces cerevisiae DNA polymerase eta transcript and protein

RAD30-encoded DNA polymerase eta functions as a translesion polymerase that can bypass the most frequent types of UV-induced pyrimidine photoproducts in an error-free manner. Although its transcript is UV-inducible in Saccharomyces cerevisiae, Rad30 (studied as a Rad30-Myc fusion) is a stable protein whose levels do not fluctuate following UV treatment or during cell cycle progression. Rad30 protein is subject to monoubiquitination whose level is upregulated in G1 and downregulated during S-phase reentry. This downregulation is accelerated in UV-treated cells. A missense mutation (L577Q) of the ubiquitin binding domain (UBZ) confers a reduced degree of ubiquitination outside of G1 and a complete failure to stably interact with ubiquitinated substrates. This mutation confers a phenotype resembling a complete RAD30 deletion, thus attesting to the significance of the UBZ motif for polymerase eta function in vivo.