Differentiating mechanisms of toxicity using global gene expression analysis in Saccharomyces cerevisiae

Genotoxicity Assessment of Haloacetaldehyde Disinfection Byproducts via a Simplified Yeast-Based Toxicogenomics Assay

Haloacetaldehydes (HALs) represent the third-largest category of disinfection byproducts (DBPs) in drinking water in terms of weight. As a subset of unregulated DBPs, only a few HALs have undergone assessment, yielding limited information regarding their genotoxicity mechanisms. Herein, we developed a simplified yeast-based toxicogenomics assay to evaluate the genotoxicity of five specific HALs. This assay recorded the protein expression profiles of eight Saccharomyces cerevisiae strains fused with green fluorescent protein, including all known DNA damage and repair pathways. High-resolution real-time pathway activation data and protein expression profiles in conjunction with clustering analysis revealed that the five HALs induced various DNA damage and repair pathways. Among these, chloroacetaldehyde and trichloroacetaldehyde were found to be positively associated with genotoxicity, while dichloroacetaldehyde, bromoacetaldehyde, and tribromoacetaldehyde displayed negative associations. The protein effect level index, which are molecular end points derived from a toxicogenomics assay, exhibited a statistically significant positive correlation with the results of traditional genotoxicity assays, such as the comet assay (rp = 0.830 and p < 0.001) and SOS/umu assay (rp = 0.786 and p = 0.004). This yeast-based toxicogenomics assay, which employs a minimal set of gene biomarkers, can be used for mechanistic genotoxicity screening and assessment of HALs and other chemical compounds. These results contribute to bridging the knowledge gap regarding the molecular mechanisms underlying the genotoxicity of HALs and enable the categorization of HALs based on their distinct DNA damage and repair mechanisms.

Stability of Rad51 recombinase and persistence of Rad51 DNA repair foci depends on post-translational modifiers, ubiquitin and SUMO

The DNA double-strand breaks are particularly deleterious, especially when an error-free repair pathway is unavailable, enforcing the error-prone recombination pathways to repair the lesion. Cells can resume the cell cycle but at the expense of decreased viability due to genome rearrangements. One of the major players involved in recombinational repair of DNA damage is Rad51 recombinase, a protein responsible for presynaptic complex formation. We previously showed that an increased level of this protein promotes the usage of illegitimate recombination. Here we show that the level of Rad51 is regulated via the ubiquitin-dependent proteolytic pathway. The ubiquitination of Rad51 depends on multiple E3 enzymes, including SUMO-targeted ubiquitin ligases. We also demonstrate that Rad51 can be modified by both ubiquitin and SUMO. Moreover, its modification with ubiquitin may lead to opposite effects: degradation dependent on Rad6, Rad18, Slx8, Dia2, and the anaphase-promoting complex, or stabilization dependent on Rsp5. We also show that post-translational modifications with SUMO and ubiquitin affect Rad51's ability to form and disassemble DNA repair foci, respectively, influencing cell cycle progression and cell viability in genotoxic stress conditions. Our data suggest the existence of a complex E3 ligases network that regulates Rad51 recombinase's turnover, its molecular activity, and access to DNA, limiting it to the proportions optimal for the actual cell cycle stage and growth conditions, e.g., stress. Dysregulation of this network would result in a drop in cell viability due to uncontrolled genome rearrangement in the yeast cells. In mammals would promote the development of genetic diseases and cancer.

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.

Endosomal trafficking and DNA damage checkpoint kinases dictate survival to replication stress by regulating amino acid uptake and protein synthesis

Atg6^{Beclin 1} mediates autophagy and endosomal trafficking. We investigated how Atg6 influences replication stress. Combining genetic, genomic, metabolomic, and proteomic approaches, we found that the Vps34-Vps15-Atg6^{Beclin 1}-Vps38^UVRAG-phosphatydilinositol-3 phosphate (PtdIns(3)P) axis sensitizes cells to replication stress by favoring the degradation of plasma membrane amino acid (AA) transporters via endosomal trafficking and ESCRT proteins, while the PtdIns(3)P phosphatases Ymr1 and Inp53 promote survival to replication stress by reversing this process. An impaired AA uptake triggers activation of Gcn2, which attenuates protein synthesis by phosphorylating eIF2α. Mec1^Atr-Rad53^Chk1/Chk2 activation during replication stress further hinders translation efficiency by counteracting eIF2α dephosphorylation through Glc7^PP1. AA shortage-induced hyperphosphorylation of eIF2α inhibits the synthesis of 65 stress response proteins, thus resulting in cell sensitization to replication stress, while TORC1 promotes cell survival. Our findings reveal an integrated network mediated by endosomal trafficking, translational control pathways, and checkpoint kinases linking AA availability to the response to replication stress.

Genetic interaction between F-box motif encoding YDR131C and retrograde signaling-related RTG1 regulates the stress response and apoptosis in Saccharomyces cerevisiae

The retrograde signaling pathway is well conserved from yeast to humans, which regulates cell adaptation during stress conditions and prevents cell death. One of its components, RTG1 encoded Rtg1p in association with Rtg3p communicates between mitochondria, nucleus, and peroxisome during stress for adaptation, by regulation of transcription. The F-box motif protein encoded by YDR131C  constitutes a part of SCF ^Ydr131c -E3 ligase complex, with unknown function; however, it is known that retrograde signaling is modulated by the E3 ligase complex. This study reports epistasis interaction between YDR131C and RTG1, which regulates cell growth, response to genotoxic stress, decreased apoptosis, resistance to petite mutation, and cell wall integrity. The cells of ydr131cΔrtg1Δ genetic background exhibits growth rate improvement however, sensitivity to hydroxyurea, itraconazole antifungal agent and synthetic indoloquinazoline-based alkaloid (8-fluorotryptanthrin, RK64), which disrupts the cell wall integrity in Saccharomyces cerevisiae. The epistatic interaction between YDR131C and RTG1 indicates a link between protein degradation and retrograde signaling pathways.

Host Chromatin Regulators Required for Aggregatibacter actinomycetemcomitans Cytolethal Distending Toxin Activity in Saccharomyces cerevisiae Model

Cytolethal distending toxin (CDT) is a bacterial genotoxin that causes host cell cycle arrest and death. We previously employed a Saccharomyces cerevisiae model with inducible expression of the CDT catalytic subunit from Aggregatibacter actinomycetemcomitans, AaCdtB, and showed that a wide variety of host factors play a role in facilitating the activity of CdtB. Our observation that a yeast H2B mutant defective in chromatin condensation was partially resistant to CdtB implies that chromatin structure may affect CDT function. In this study, we identified host chromatin regulatory genes required for CdtB cytotoxicity. We found that the deletion of HTZ1 or certain subunits of SWR, INO80, and SIR complexes increased cellular resistance to CdtB. We hypothesized that CdtB may interact with Htz1 or the chromatin, but immunoprecipitation experiments failed to detect physical interaction between CdtB and Htz1 or the chromatin. However, we observed reduced nuclear localization of CdtB in several mutants, suggesting that impaired nuclear translocation may, at least partly, explain the mechanisms of CdtB resistance. In addition, mutations in chromatin regulatory genes induce changes in the global gene expression profile, and these may indirectly affect CdtB toxicity. Our results suggest that decreased expression of endoplasmic reticulum (ER)-Golgi transport-related genes that may be involved in CdtB transport and/or increased expression of DNA repair genes may contribute to CdtB resistance. These results suggest that the functions of chromatin regulators may contribute to the activity of CDT in host cells.

Microbial bioassays in environmental toxicity testing

Accidental spills and the misuse of chemicals may lead to current and legacy environmental contamination and pose concerns over possible (eco)toxicological secondary effects and risks toward non-target microbes and higher eukaryotes, including humans, in ecosystems. In the last decades, scientists and regulators have faced requests to thoroughly screen, prioritize and predict the possible deleterious effects of the huge numbers of existing and emerging xenobiotics, wastewaters and environmental samples on biological systems. In this context, it has become necessary to develop and validate (eco)toxicity bioassays based on microorganisms (e.g., bacteria, microalga, yeast, filamentous fungi, protozoa) as test-organisms whose data should be meaningful for environmental (micro)organisms that may be exposed to contaminated environments. These generally simple, fast and cost-effective bioassays may be preliminary and complementary to the more complex and long-term whole-organism animal-based traditional ecotoxicity tests. With the goal of highlighting the potential offered by microbial-based bioassays as non-animal alternatives in (eco)toxicity testing, the present chapter provides an overview of the current state-of-the art in the development and use of microbial toxicity bioassays through the examination of relatively recent examples with a diverse range of toxicity endpoints. It goes into the (eco)toxicological relevance of these bioassays, ranging from the more traditional microalga- and bacterial-based assays already accepted at regulatory level and commercially available to the more innovative microbial transcriptional profiling and gene expression bioassays, including some examples of biosensors.

Toxicological response of the model fungus Saccharomyces cerevisiae to different concentrations of commercial graphene nanoplatelets

Graphene nanomaterials have attracted a great interest during the last years for different applications, but their possible impact on different biological systems remains unclear. Here, an assessment to understand the toxicity of commercial polycarboxylate functionalized graphene nanoplatelets (GN) on the unicellular fungal model Saccharomyces cerevisiae was performed. While cell proliferation was not negatively affected even in the presence of 800 mg L-1 of the nanomaterial for 24 hours, oxidative stress was induced at a lower concentration (160 mg L-1), after short exposure periods (2 and 4 hours). No DNA damage was observed under a comet assay analysis under the studied conditions. In addition, to pinpoint the molecular mechanisms behind the early oxidative damage induced by GN and to identify possible toxicity pathways, the transcriptome of S. cerevisiae exposed to 160 and 800 mg L-1 of GN was studied. Both GN concentrations induced expression changes in a common group of genes (337), many of them related to the fungal response to reduce the nanoparticles toxicity and to maintain cell homeostasis. Also, a high number of genes were only differentially expressed in the GN800 condition (3254), indicating that high GN concentrations can induce severe changes in the physiological state of the yeast.

SWI/SNF chromatin remodelling complex contributes to clearance of cytoplasmic protein aggregates and regulates unfolded protein response in Saccharomyces cerevisiae

Chromatin remodelling complexes are multi-subunit assemblies, each containing a catalytic ATPase and translocase that is capable of mobilizing nucleosomes to alter the chromatin structure. SWI/SNF remodelling complexes with higher DNA translocation efficiency evict histones or slide the nucleosomes away from each other making DNA accessible for transcription and repair machinery. Chromatin remodelling at the promoter of stress-responsive genes by SWI/SNF becomes necessary during the heat and proteotoxic stress. While the involvement of SWI/SNF in transcription of stress-responsive genes has been studied extensively, the regulation of proteostasis by SWI/SNF is not well understood. This study demonstrates critical functions of SWI/SNF in response to cadmium-induced proteotoxic stress. Deletion of either ATPase-translocase subunit of SWI/SNF complex (Swi2/Snf2) or a regulatory subunit Swi3 abrogates the clearance of cadmium-induced protein aggregates. Our results suggest that Snf2 and Swi3 regulate the protein folding in endoplasmic reticulum (ER) that reduces the chances of forming unfolded protein aggregates under the proteotoxic stress of cadmium. The Ire1-mediated unfolded protein response (UPR) maintains ER homeostasis by upregulating the expression of chaperones and ER-associated degradation (ERAD) components. We found that Snf2 maintains normal oxidative environment essential for Ire1 activity. Deletion of SNF2 reduced the Ire1 activity and UPR, indicating involvement of Snf2 in Ire1-mediated ER proteostasis. Together, these findings suggest that SWI/SNF complex regulates ER homeostasis and protein folding crucial for tolerating proteotoxic stress.

Coordination of the Cell Wall Integrity and High-Osmolarity Glycerol Pathways in Response to Ethanol Stress in Saccharomyces cerevisiae

During fermentation, a high ethanol concentration is a major stress that influences the vitality and viability of yeast cells, which in turn leads to a termination of the fermentation process. In this study, we show that the BCK1 and SLT2 genes encoding mitogen-activated protein kinase kinase kinase (MAPKKK) and mitogen-activated protein kinase (MAPK) of the cell wall integrity (CWI) pathway, respectively, are essential for ethanol tolerance, suggesting that the CWI pathway is involved in the response to ethanol-induced cell wall stress. Upon ethanol exposure, the CWI pathway induces the expression of specific cell wall-remodeling genes, including FKS2, CRH1, and PIR3 (encoding β-1,3-glucan synthase, chitin transglycosylase, and O-glycosylated cell wall protein, respectively), which eventually leads to the remodeling of the cell wall structure. Our results revealed that in response to ethanol stress, the high-osmolarity glycerol (HOG) pathway plays a collaborative role with the CWI pathway in inducing cell wall remodeling via the upregulation of specific cell wall biosynthesis genes such as the CRH1 gene. Furthermore, the substantial expression of CWI-responsive genes is also triggered by external hyperosmolarity, suggesting that the adaptive changes in the cell wall are crucial for protecting yeast cells against not only cell wall stress but also osmotic stress. On the other hand, the cell wall stress-inducing agent calcofluor white has no effect on promoting the expression of GPD1, a major target gene of the HOG pathway. Collectively, these findings suggest that during ethanol stress, the CWI and HOG pathways collaboratively regulate the transcription of specific cell wall biosynthesis genes, thereby leading to adaptive changes in the cell wall.IMPORTANCE The budding yeast Saccharomyces cerevisiae has been widely used in industrial fermentations, including the production of alcoholic beverages and bioethanol. During fermentation, an increased ethanol concentration is the main stress that affects yeast metabolism and inhibits ethanol production. This work presents evidence that in response to ethanol stress, both CWI and HOG pathways cooperate to control the expression of cell wall-remodeling genes in order to build the adaptive strength of the cell wall. These findings will contribute to a better understanding of the molecular mechanisms underlying adaptive responses and tolerance of yeast to ethanol stress, which is essential for successful engineering of yeast strains for improved ethanol tolerance.

Transcriptional responses to DNA damage

In response to the threat of DNA damage, cells exhibit a dramatic and multi-factorial response spanning from transcriptional changes to protein modifications, collectively known as the DNA damage response (DDR). Here, we review the literature surrounding the transcriptional response to DNA damage. We review differences in observed transcriptional responses as a function of cell cycle stage and emphasize the importance of experimental design in these transcriptional response studies. We additionally consider topics including structural challenges in the transcriptional response to DNA damage as well as the connection between transcription and protein abundance.

TGMI: an efficient algorithm for identifying pathway regulators through evaluation of triple-gene mutual interaction

Despite their important roles, the regulators for most metabolic pathways and biological processes remain elusive. Presently, the methods for identifying metabolic pathway and biological process regulators are intensively sought after. We developed a novel algorithm called triple-gene mutual interaction (TGMI) for identifying these regulators using high-throughput gene expression data. It first calculated the regulatory interactions among triple gene blocks (two pathway genes and one transcription factor (TF)), using conditional mutual information, and then identifies significantly interacted triple genes using a newly identified novel mutual interaction measure (MIM), which was substantiated to reflect strengths of regulatory interactions within each triple gene block. The TGMI calculated the MIM for each triple gene block and then examined its statistical significance using bootstrap. Finally, the frequencies of all TFs present in all significantly interacted triple gene blocks were calculated and ranked. We showed that the TFs with higher frequencies were usually genuine pathway regulators upon evaluating multiple pathways in plants, animals and yeast. Comparison of TGMI with several other algorithms demonstrated its higher accuracy. Therefore, TGMI will be a valuable tool that can help biologists to identify regulators of metabolic pathways and biological processes from the exploded high-throughput gene expression data in public repositories.

Yeast-based genotoxicity tests for assessing DNA alterations and DNA stress responses: a 40-year overview

By damaging DNA molecules, genotoxicants cause genetic mutations and also increase human susceptibility to cancers and genetic diseases. Over the past four decades, several assays have been developed in the budding yeast Saccharomyces cerevisiae to screen potential genotoxic substances and provide alternatives to animal-based genotoxicity tests. These yeast-based genotoxicity tests are either DNA alteration-based or DNA stress-response reporter-based. The former, which came first, were developed from the genetic studies conducted on various types of DNA alterations in yeast cells. Despite their limited throughput capabilities, some of these tests have been used as short-term genotoxicity tests in addition to bacteria- or mammalian cell-based tests. In contrast, the latter tests are based on the emergent transcriptional induction of DNA repair-related genes via activation of the DNA damage checkpoint kinase cascade triggered by DNA damage. Some of these reporter assays have been linked to DNA damage-responsive promoters to assess chemical carcinogenicity and ecotoxicity in environmental samples. Yeast-mediated genotoxicity tests are being continuously improved by increasing the permeability of yeast cell walls, by the ectopic expression of mammalian cytochrome P450 systems, by the use of DNA repair-deficient host strains, and by integrating them into high-throughput formats or microfluidic devices. Notably, yeast-based reporter assays linked with the newer toxicogenomic approaches are becoming powerful short-term genotoxicity tests for large numbers of compounds. These tests can also be used to detect polluted environmental samples, and as effective screening tools during anticancer drug development.

Unique molecular mechanisms for maintenance and alteration of genetic information in the budding yeast Saccharomyces cerevisiae

The high-fidelity transmission of genetic information is crucial for the survival of organisms, the cells of which have the ability to protect DNA against endogenous and environmental agents, including reactive oxygen species (ROS), ionizing radiation, and various chemical compounds. The basis of protection mechanisms has been evolutionarily conserved from yeast to humans; however, each organism often has a specialized mode of regulation that uses different sets of machineries, particularly in lower eukaryotes. The divergence of molecular mechanisms among related organisms has provided insights into the evolution of cellular machineries to a higher architecture. Uncommon characteristics of machineries may also contribute to the development of new applications such as drugs with novel mechanisms of action. In contrast to the cellular properties for maintaining genetic information, living organisms, particularly microbes, inevitably undergo genetic alterations in order to adapt to environmental conditions. The maintenance and alteration of genetic information may be inextricably linked to each other. In this review, we describe recent findings on the unconventional molecular mechanisms of DNA damage response and DNA double-strand break (DSB) repair in the budding yeast Saccharomyces cerevisiae. We also introduce our previous research on genetic and phenotypic instabilities observed in a clonal population of clinically-derived S. cerevisiae. The molecular mechanisms of this case were associated with mutations to generate tyrosine-inserting tRNA-Tyr ochre suppressors and the position effects of mutation frequencies among eight tRNA-Tyr loci dispersed in the genome. Phenotypic variations among different strain backgrounds have also been observed by another type of nonsense suppressor, the aberrant form of the translation termination factor. Nonsense suppressors are considered to be responsible for the genome-wide translational readthrough of termination codons, including natural nonsense codons. The nonsense suppressor-mediated acquisition of phenotypic variations may be advantageous for adaptation to environmental conditions and survival during evolution.

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

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

Crystal Structure of Saccharomyces cerevisiae ECM4, a Xi-Class Glutathione Transferase that Reacts with Glutathionyl-(hydro)quinones

Glutathionyl-hydroquinone reductases (GHRs) belong to the recently characterized Xi-class of glutathione transferases (GSTXs) according to unique structural properties and are present in all but animal kingdoms. The GHR ScECM4 from the yeast Saccharomyces cerevisiae has been studied since 1997 when it was found to be potentially involved in cell-wall biosynthesis. Up to now and in spite of biological studies made on this enzyme, its physiological role remains challenging. The work here reports its crystallographic study. In addition to exhibiting the general GSTX structural features, ScECM4 shows extensions including a huge loop which contributes to the quaternary assembly. These structural extensions are probably specific to Saccharomycetaceae. Soaking of ScECM4 crystals with GS-menadione results in a structure where glutathione forms a mixed disulfide bond with the cysteine 46. Solution studies confirm that ScECM4 has reductase activity for GS-menadione in presence of glutathione. Moreover, the high resolution structures allowed us to propose new roles of conserved residues of the active site to assist the cysteine 46 during the catalytic act.

DNA Damage Response Checkpoint Activation Drives KP1019 Dependent Pre-Anaphase Cell Cycle Delay in S. cerevisiae

Careful regulation of the cell cycle is required for proper replication, cell division, and DNA repair. DNA damage–including that induced by many anticancer drugs–results in cell cycle delay or arrest, which can allow time for repair of DNA lesions. Although its molecular mechanism of action remains a matter of debate, the anticancer ruthenium complex KP1019 has been shown to bind DNA in biophysical assays and to damage DNA of colorectal and ovarian cancer cells in vitro. KP1019 has also been shown to induce mutations and induce cell cycle arrest in Saccharomyces cerevisiae, suggesting that budding yeast can serve as an appropriate model for characterizing the cellular response to the drug. Here we use a transcriptomic approach to verify that KP1019 induces the DNA damage response (DDR) and find that KP1019 dependent expression of HUG1 requires the Dun1 checkpoint; both consistent with KP1019 DDR in budding yeast. We observe a robust KP1019 dependent delay in cell cycle progression as measured by increase in large budded cells, 2C DNA content, and accumulation of Pds1 which functions to inhibit anaphase. Importantly, we also find that deletion of RAD9, a gene required for the DDR, blocks drug-dependent changes in cell cycle progression, thereby establishing a causal link between the DDR and phenotypes induced by KP1019. Interestingly, yeast treated with KP1019 not only delay in G2/M, but also exhibit abnormal nuclear position, wherein the nucleus spans the bud neck. This morphology correlates with short, misaligned spindles and is dependent on the dynein heavy chain gene DYN1. We find that KP1019 creates an environment where cells respond to DNA damage through nuclear (transcriptional changes) and cytoplasmic (motor protein activity) events.

Cross-Talk between Carbon Metabolism and the DNA Damage Response in S. cerevisiae

Yeast cells with DNA damage avoid respiration, presumably because products of oxidative metabolism can be harmful to DNA. We show that DNA damage inhibits the activity of the Snf1 (AMP-activated) protein kinase (AMPK), which activates expression of genes required for respiration. Glucose and DNA damage upregulate SUMOylation of Snf1, catalyzed by the SUMO E3 ligase Mms21, which inhibits SNF1 activity. The DNA damage checkpoint kinases Mec1/ATR and Tel1/ATM, as well as the nutrient-sensing protein kinase A (PKA), regulate Mms21 activity toward Snf1. Mec1 and Tel1 are required for two SNF1-regulated processes-glucose sensing and ADH2 gene expression-even without exogenous genotoxic stress. Our results imply that inhibition of Snf1 by SUMOylation is a mechanism by which cells lower their respiration in response to DNA damage. This raises the possibility that activation of DNA damage checkpoint mechanisms could contribute to aerobic fermentation (Warburg effect), a hallmark of cancer cells.

Bicluster Sampled Coherence Metric (BSCM) provides an accurate environmental context for phenotype predictions

Background

Biclustering is a popular method for identifying under which experimental conditions biological signatures are co-expressed. However, the general biclustering problem is NP-hard, offering room to focus algorithms on specific biological tasks. We hypothesize that conditional co-regulation of genes is a key factor in determining cell phenotype and that accurately segregating conditions in biclusters will improve such predictions. Thus, we developed a bicluster sampled coherence metric (BSCM) for determining which conditions and signals should be included in a bicluster.

Results

Our BSCM calculates condition and cluster size specific p-values, and we incorporated these into the popular integrated biclustering algorithm cMonkey. We demonstrate that incorporation of our new algorithm significantly improves bicluster co-regulation scores (p-value = 0.009) and GO annotation scores (p-value = 0.004). Additionally, we used a bicluster based signal to predict whether a given experimental condition will result in yeast peroxisome induction. Using the new algorithm, the classifier accuracy improves from 41.9% to 76.1% correct.

Conclusions

We demonstrate that the proposed BSCM helps determine which signals ought to be co-clustered, resulting in more accurately assigned bicluster membership. Furthermore, we show that BSCM can be extended to more accurately detect under which experimental conditions the genes are co-clustered. Features derived from this more accurate analysis of conditional regulation results in a dramatic improvement in the ability to predict a cellular phenotype in yeast. The latest cMonkey is available for download at https://github.com/baliga-lab/cmonkey2. The experimental data and source code featured in this paper is available http://AitchisonLab.com/BSCM. BSCM has been incorporated in the official cMonkey release.

Quantitative proteomics of the yeast Hsp70/Hsp90 interactomes during DNA damage reveal chaperone-dependent regulation of ribonucleotide reductase

The highly conserved molecular chaperones Hsp90 and Hsp70 are indispensible for folding and maturation of a significant fraction of the proteome, including many proteins involved in signal transduction and stress response. To examine the dynamics of chaperone-client interactions after DNA damage, we applied quantitative affinity-purification mass spectrometry (AP-MS) proteomics to characterize interactomes of the yeast Hsp70 isoform Ssa1 and Hsp90 isoform Hsp82 before and after exposure to methyl methanesulfonate. Of 256 proteins identified and quantified via (16)O(/18)O labeling and LC-MS/MS, 142 are novel Hsp70/90 interactors. Nearly all interactions remained unchanged or decreased after DNA damage, but 5 proteins increased interactions with Ssa1 and/or Hsp82, including the ribonucleotide reductase (RNR) subunit Rnr4. Inhibiting Hsp70 or 90 chaperone activity destabilized Rnr4 in yeast and its vertebrate homolog hRMM2 in breast cancer cells. In turn, pre-treatment of cancer cells with chaperone inhibitors sensitized cells to the RNR inhibitor gemcitabine, suggesting a novel chemotherapy strategy. All MS data have been deposited in the ProteomeXchange with identifier PXD001284.

BIOLOGICAL SIGNIFICANCE

This study provides the dynamic interactome of the yeast Hsp70 and Hsp90 under DNA damage which suggest key roles for the chaperones in a variety of signaling cascades. Importantly, the cancer drug target ribonucleotide reductase was shown to be a client of Hsp70 and Hsp90 in both yeast and breast cancer cells. As such, this study highlights the potential of a novel cancer therapeutic strategy that exploits the synergy of chaperone and ribonucleotide reductase inhibitors.

Sky1 regulates the expression of sulfur metabolism genes in response to cisplatin

Cisplatin is commonly used in cancer therapy and yeast cells are also sensitive to this compound. We present a transcriptome analysis discriminating between RNA changes induced by cisplatin treatment, which are dependent on or independent of SKY1 function – a gene whose deletion increases resistance to the drug. Gene expression changes produced by addition of cisplatin to W303 and W303-Δsky1 cells were recorded using DNA microarrays. The data, validated by quantitative PCR, revealed 122 differentially expressed genes: 69 upregulated and 53 downregulated. Among the upregulated genes, those related to sulfur metabolism were over-represented and partially dependent on Sky1. Deletions of MET4 or other genes encoding co-regulators of the expression of sulfur-metabolism-related genes, with the exception of MET28, did not modify the cisplatin sensitivity of yeast cells. One of the genes with the highest cisplatin-induced upregulation was SEO1, encoding a putative permease of sulfur compounds. We also measured the platinum, sulfur and glutathione content in W303, W303-Δsky1 and W303-Δseo1 cells after cisplatin treatment, and integration of the data suggested that these transcriptional changes might represent a cellular response that allowed chelation of cisplatin with sulfur-containing amino acids and also helped DNA repair by stimulating purine biosynthesis. The transcription pattern of stimulation of sulfur-containing amino acids and purine synthesis decreased, or even disappeared, in the W303-Δsky1 strain.

Dissection of DNA damage responses using multiconditional genetic interaction maps

To protect the genome, cells have evolved a diverse set of pathways designed to sense, signal, and repair multiple types of DNA damage. To assess the degree of coordination and crosstalk among these pathways, we systematically mapped changes in the cell's genetic network across a panel of different DNA-damaging agents, resulting in ~1,800,000 differential measurements. Each agent was associated with a distinct interaction pattern, which, unlike single-mutant phenotypes or gene expression data, has high statistical power to pinpoint the specific repair mechanisms at work. The agent-specific networks revealed roles for the histone acetyltranferase Rtt109 in the mutagenic bypass of DNA lesions and the neddylation machinery in cell-cycle regulation and genome stability, while the network induced by multiple agents implicates Irc21, an uncharacterized protein, in checkpoint control and DNA repair. Our multiconditional genetic interaction map provides a unique resource that identifies agent-specific and general DNA damage response pathways.

Identification of the molecular mechanisms underlying the cytotoxic action of a potent platinum metallointercalator

Platinum-based DNA metallointercalators are structurally different from the covalent DNA binders such as cisplatin and its derivatives but have potent in vitro activity in cancer cell lines. However, limited understanding of their molecular mechanisms of cytotoxic action greatly hinders their further development as anticancer agents. In this study, a lead platinum-based metallointercalator, [(5,6-dimethyl-1,10-phenanthroline) (1S,2S-diaminocyclohexane)platinum(II)](2+) (56MESS) was found to be 163-fold more active than cisplatin in a cisplatin-resistant cancer cell line. By using transcriptomics in a eukaryotic model organism, yeast Saccharomyces cerevisiae, we identified 93 genes that changed their expressions significantly upon exposure of 56MESS in comparison to untreated controls (p ≤ 0.05). Bioinformatic analysis of these genes demonstrated that iron and copper metabolism, sulfur-containing amino acids and stress response were involved in the cytotoxicity of 56MESS. Follow-up experiments showed that the iron and copper concentrations were much lower in 56MESS-treated cells compared to controls as measured by inductively coupled plasma optical emission spectrometry. Deletion mutants of the key genes in the iron and copper metabolism pathway and glutathione synthesis were sensitive to 56MESS. Taken together, the study demonstrated that the cytotoxic action of 56MESS is mediated by its ability to disrupt iron and copper metabolism, suppress the biosynthesis of sulfur-containing amino acids and attenuate cellular defence capacity. As these mechanisms are in clear contrast to the DNA binding mechanism for cisplatin and its derivative, 56MESS may be able to overcome cisplatin-resistant cancers. These findings have provided basis to further develop the platinum-based metallointercalators as anticancer agents.

ELECTRONIC SUPPLEMENTARY MATERIAL

The online version of this article (doi:10.1007/s12154-011-0070-x) contains supplementary material, which is available to authorized users.

S-Glutathionyl-(chloro)hydroquinone reductases: a novel class of glutathione transferases

Sphingobium chlorophenolicum completely mineralizes PCP (pentachlorophenol). Two GSTs (glutathione transferases), PcpC and PcpF, are involved in the degradation. PcpC uses GSH to reduce TeCH (tetrachloro-p-hydroquinone) to TriCH (trichloro-p-hydroquinone) and then to DiCH (dichloro-p-hydroquinone) during PCP degradation. However, oxidatively damaged PcpC produces GS-TriCH (S-glutathionyl-TriCH) and GS-DiCH (S-glutathionyl-TriCH) conjugates. PcpF converts the conjugates into TriCH and DiCH, re-entering the degradation pathway. PcpF was further characterized in the present study. It catalysed GSH-dependent reduction of GS-TriCH via a Ping Pong mechanism. First, PcpF reacted with GS-TriCH to release TriCH and formed disulfide bond between its Cys53 residue and the GS moiety. Then, a GSH came in to regenerate PcpF and release GS-SG. A TBLASTN search revealed that PcpF homologues were widely distributed in bacteria, halobacteria (archaea), fungi and plants, and they belonged to ECM4 (extracellular mutant 4) group COG0435 in the conserved domain database. Phylogenetic analysis grouped PcpF and homologues into a distinct group, separated from Omega class GSTs. The two groups shared conserved amino acid residues, for GSH binding, but had different residues for the binding of the second substrate. Several recombinant PcpF homologues and two human Omega class GSTs were produced in Escherichia coli and purified. They had zero or low activities for transferring GSH to standard substrates, but all had reasonable activities for GSH-dependent reduction of disulfide bond (thiol transfer), dehydroascorbate and dimethylarsinate. All the tested PcpF homologues reduced GS-TriCH, but the two Omega class GSTs did not. Thus PcpF homologues were tentatively named S-glutathionyl-(chloro)hydroquinone reductases for catalysing the GSH-dependent reduction of GS-TriCH.

Graph ranking for exploratory gene data analysis

Background

Microarray technology has made it possible to simultaneously monitor the expression levels of thousands of genes in a single experiment. However, the large number of genes greatly increases the challenges of analyzing, comprehending and interpreting the resulting mass of data. Selecting a subset of important genes is inevitable to address the challenge. Gene selection has been investigated extensively over the last decade. Most selection procedures, however, are not sufficient for accurate inference of underlying biology, because biological significance does not necessarily have to be statistically significant. Additional biological knowledge needs to be integrated into the gene selection procedure.

Results

We propose a general framework for gene ranking. We construct a bipartite graph from the Gene Ontology (GO) and gene expression data. The graph describes the relationship between genes and their associated molecular functions. Under a species condition, edge weights of the graph are assigned to be gene expression level. Such a graph provides a mathematical means to represent both species-independent and species-dependent biological information. We also develop a new ranking algorithm to analyze the weighted graph via a kernelized spatial depth (KSD) approach. Consequently, the importance of gene and molecular function can be simultaneously ranked by a real-valued measure, KSD, which incorporates the global and local structure of the graph. Over-expressed and under-regulated genes also can be separately ranked.

Conclusion

The gene-function bigraph integrates molecular function annotations into gene expression data. The relevance of genes is described in the graph (through a common function). The proposed method provides an exploratory framework for gene data analysis.

Carbohydrate utilization and the lager yeast transcriptome during brewery fermentation

The fermentable carbohydrate composition of wort and the manner in which it is utilized by yeast during brewery fermentation have a direct influence on fermentation efficiency and quality of the final product. In this study the response of a brewing yeast strain to changes in wort fermentable carbohydrate concentration and composition during full-scale (3275 hl) brewery fermentation was investigated by measuring transcriptome changes with the aid of oligonucleotide-based DNA arrays. Up to 74% of the detectable genes showed a significant (p</=0.01) differential expression pattern during fermentation and the majority of these genes showed transient or prolonged peaks in expression following the exhaustion of the monosaccharides from the wort. Transcriptional activity of many genes was consistent with their known responses to glucose de/repression under laboratory conditions, despite the presence of di- and trisaccharide sugars in the wort. In a number of cases the transcriptional response of genes was not consistent with their known responses to glucose, suggesting a degree of complexity during brewery fermentation which cannot be replicated in small-scale wort fermentations or in laboratory experiments involving defined media.

Improving prediction of chemical carcinogenicity by considering multiple mechanisms and applying toxicogenomic approaches

While scientific knowledge of the potential health significance of chemical exposures has grown, experimental methods for predicting the carcinogenicity of environmental agents have not been substantially updated in the last two decades. Current methodologies focus first on identifying genotoxicants under the premise that agents capable of directly damaging DNA are most likely to be carcinogenic to humans. Emphasis on the distinction between genotoxic and non-genotoxic carcinogens is also motivated by assumed implications for the dose-response curve; it is purported that genotoxicants would lack a threshold in the low dose region, in contrast to non-genotoxic agents. However, for the vast majority of carcinogens, little if any empirical data exist to clarify the nature of the cancer dose-response relationship at low doses in the exposed human population. Recent advances in scientific understanding of cancer biology-and increased appreciation of the multiple impacts of carcinogens on this disease process-support the view that environmental chemicals can act through multiple toxicity pathways, modes and/or mechanisms of action to induce cancer and other adverse health outcomes. Moreover, the relationship between dose and a particular outcome in an individual could take multiple forms depending on genetic background, target tissue, internal dose and other factors besides mechanisms or modes of action; inter-individual variability and susceptibility in response are, in turn, key determinants of the population dose-response curve. New bioanalytical approaches (e.g., transcriptomics, proteomics, and metabolomics) applied in human, animal and in vitro studies could better characterize a wider array of hazard traits and improve the ability to predict the potential carcinogenicity of chemicals.

The anti-tumour agent, cisplatin, and its clinically ineffective isomer, transplatin, produce unique gene expression profiles in human cells

Cisplatin is a DNA-damaging anti-cancer agent that is widely used to treat a range of tumour types. Despite its clinical success, cisplatin treatment is still associated with a number of dose-limiting toxic side effects. The purpose of this study was to clarify the molecular events that are important in the anti-tumour activity of cisplatin, using gene expression profiling techniques. Currently, our incomplete understanding of this drug's mechanism of action hinders the development of more efficient and less harmful cisplatin-based chemotherapeutics. In this study the effect of cisplatin on gene expression in human foreskin fibroblasts has been investigated using human 19K oligonucleotide microarrays. In addition its clinically inactive isomer, transplatin, was also tested. Dualfluor microarray experiments comparing treated and untreated cells were performed in quadruplicate. Cisplatin treatment was shown to significantly up- or down-regulate a consistent subset of genes. Many of these genes responded similarly to treatment with transplatin, the therapeutically inactive isomer of cisplatin. However, a smaller proportion of these transcripts underwent differential expression changes in response to the two isomers. Some of these genes may constitute part of the DNA damage response induced by cisplatin that is critical for its anti-tumour activity. Ultimately, the identification of gene expression responses unique to clinically active compounds, like cisplatin, could thus greatly benefit the design and development of improved chemotherapeutics.

A systems approach to delineate functions of paralogous transcription factors: role of the Yap family in the DNA damage response

Duplication of genes encoding transcription factors plays an essential role in driving phenotypic variation. Because regulation can occur at multiple levels, it is often difficult to discern how each duplicated factor achieves its regulatory specificity. In these cases, a "systems approach" may distinguish the role of each factor by integrating complementary large-scale measurements of the regulatory network. To explore such an approach, we integrate growth phenotypes, promoter binding profiles, and gene expression patterns to model the DNA damage response network controlled by the Yeast-specific AP-1 (YAP) family of transcription factors. This analysis reveals that YAP regulatory specificity is achieved by at least three mechanisms: (i) divergence of DNA-binding sequences into two subfamilies; (ii) condition-specific combinatorial regulation by multiple Yap factors; and (iii) interactions of Yap 1, 4, and 6 with chromatin remodeling proteins. Additional microarray experiments establish that Yap 4 and 6 regulate gene expression through interactions with the histone deacetylase, Hda1. The data further highlight differences among Yap paralogs in terms of their regulatory mode of action (activation vs. repression). This study suggests how other large TF families might be disentangled in the future.

Identification of genes affecting hydrogen sulfide formation in Saccharomyces cerevisiae

A screen of the Saccharomyces cerevisiae deletion strain set was performed to identify genes affecting hydrogen sulfide (H(2)S) production. Mutants were screened using two assays: colony color on BiGGY agar, which detects the basal level of sulfite reductase activity, and production of H(2)S in a synthetic juice medium using lead acetate detection of free sulfide in the headspace. A total of 88 mutants produced darker colony colors than the parental strain, and 4 produced colonies significantly lighter in color. There was no correlation between the appearance of a dark colony color on BiGGY agar and H(2)S production in synthetic juice media. Sixteen null mutations were identified as leading to the production of increased levels of H(2)S in synthetic juice using the headspace analysis assay. All 16 mutants also produced H(2)S in actual juices. Five of these genes encode proteins involved in sulfur containing amino acid or precursor biosynthesis and are directly associated with the sulfate assimilation pathway. The remaining genes encode proteins involved in a variety of cellular activities, including cell membrane integrity, cell energy regulation and balance, or other metabolic functions. The levels of hydrogen sulfide production of each of the 16 strains varied in response to nutritional conditions. In most cases, creation of multiple deletions of the 16 mutations in the same strain did not lead to a further increase in H(2)S production, instead often resulting in decreased levels.

Specific transcriptional responses induced by 8-methoxypsoralen and UVA in yeast

Treatment of eukaryotic cells with 8-methoxypsoralen plus UVA irradiation (8-MOP/UVA) induces pyrimidine monoadducts and interstrand crosslinks and initiates a cascade of events leading to cytotoxic, mutagenic and carcinogenic responses. Transcriptional activation plays an important part in these responses. Our previous study in Saccharomyces cerevisiae showed that the repair of these lesions involves the transient formation of DNA double-strand breaks and the enhanced expression of landmark DNA damage response genes such as RAD51, RNR2 and DUN1, as well as the Mec1/Rad53 kinase signaling cascade. We have now used DNA microarrays to examine genome-wide transcriptional changes produced after induction of 8-MOP/UVA photolesions. We found that 128 genes were strongly induced and 29 genes strongly repressed. Modifications in gene expression concern numerous biological processes. Compared to other genotoxic treatments, c. 42% of the response genes were specific to 8-MOP/UVA treatment. In addition to common DNA damage response genes and genes induced by environmental stresses, a large fraction of 8-MOP/UVA response genes correspond to membrane-related functions.

Application of toxicogenomics to genetic toxicology risk assessment

Based on the assumption that compounds having similar toxic modes of action induce specific gene expression changes, the toxicity of unknown compounds can be predicted after comparison of their molecular fingerprints with those obtained with compounds of known toxicity. These predictive models will therefore rely on the characterization of marker genes. Toxicogenomics (TGX) also provides mechanistic insight into the mode of toxicity, and can therefore be used as an adjunct to the standard battery of genotoxicity tests. Promising results, highlighting the ability of TGX to differentiate genotoxic from non-genotoxic carcinogens, as well as DNA-reactive from non-DNA reactive genotoxins, have been reported. Additional data suggested the possibility of ranking genotoxins according to the nature of their interactions with DNA. This new approach could contribute to the improvement of risk assessment. TGX could be applied as a follow-up testing strategy in case of positive in vitro genotoxicity findings, and could contribute to improve our ability to identify the molecular mechanism of action and to possibly better assess dose-response curves. TGX has been found to be less sensitive than the standard genotoxicity end-points, probably because it measures the whole cell population response, when compared with standard tests designed to detect rare events in a small number of cells. Further validation will be needed (1) to better link the profiles obtained with TGX to the established genotoxicity end-points, (2) to improve the gene annotation tools, and (3) to standardise study design and data analysis and to better evaluate the impact of variability between platforms and bioinformatics approaches.

Toxicogenomics: overview and potential applications for the study of non-covalent DNA interacting chemicals

Non-covalent DNA interacting agents, DNA-groove binding chemicals and DNA intercalators, are generally considered less cytotoxic than agents producing covalent DNA adducts and other DNA damage. Although the impact of non-covalent compound-DNA interactions on convoluted molecular and biochemical pathways is not well characterized, the most prominent effects include DNA conformational and related structural perturbations, interference with normal DNA protein interactions, such as topoisomerases, as well as effects on mitochondrial DNA and function. The cellular responses to such perturbations would be expected to include changes in transcription of a variety of genes. The emerging field of toxicogenomics seeks to exploit gene responses to define expression profiling signatures for various types of drugs and toxicants, and to provide mechanistic insight into their cellular effects. There are a variety of examples whereby different classes of genotoxicants and non-genotoxic agents can be distinguished by gene expression profiling using functional genomics approaches, which survey global transcriptional responses. In this review, we will discuss the promises and precautions in the use of functional genomics approaches to characterize stress agents including non-covalent DNA interacting agents.

Early transcriptional response of wine yeast after rehydration: osmotic shock and metabolic activation

The inoculation of active dry wine yeast (ADWY) is one of the most common practices in winemaking. We have used DNA microarray technology to examine the genetic expression patterns of a commercial ADWY strain after rehydration. After rehydration of ADWY for 30 min, a further hour in water after rehydration did not lead to any relevant changes in global gene expression. Expression changes in rehydrated cells upon incubation in a sorbitol solution at the same osmotic pressure as in complete must were rather limited, whereas the presence of fermentable carbon sources or the complete medium (synthetic must) produced very similar transcriptional responses. The main responses were the activation of some genes of the fermentation pathway and of the nonoxidative branch of the pentose phosphate pathway, and the induction of a huge cluster of genes related to ribosomal biogenesis and protein synthesis. The presence of cycloheximide in fermentable medium produced a similar but stronger transcriptional response. Whereas the viabilities of rehydrated cells incubated for 1 h in these different media were similar, yeast vitality, which represents the fermentative capacity of the yeast, showed a positive correlation with the availability of a fermentable carbon source.

Analyzing the dose-dependence of the Saccharomyces cerevisiae global transcriptional response to methyl methanesulfonate and ionizing radiation

BACKGROUND

One of the most crucial tasks for a cell to ensure its long term survival is preserving the integrity of its genetic heritage via maintenance of DNA structure and sequence. While the DNA damage response in the yeast Saccharomyces cerevisiae, a model eukaryotic organism, has been extensively studied, much remains to be elucidated about how the organism senses and responds to different types and doses of DNA damage. We have measured the global transcriptional response of S. cerevisiae to multiple doses of two representative DNA damaging agents, methyl methanesulfonate (MMS) and gamma radiation.

RESULTS

Hierarchical clustering of genes with a statistically significant change in transcription illustrated the differences in the cellular responses to MMS and gamma radiation. Overall, MMS produced a larger transcriptional response than gamma radiation, and many of the genes modulated in response to MMS are involved in protein and translational regulation. Several clusters of coregulated genes whose responses varied with DNA damaging agent dose were identified. Perhaps the most interesting cluster contained four genes exhibiting biphasic induction in response to MMS dose. All of the genes (DUN1, RNR2, RNR4, and HUG1) are involved in the Mec1p kinase pathway known to respond to MMS, presumably due to stalled DNA replication forks. The biphasic responses of these genes suggest that the pathway is induced at lower levels as MMS dose increases. The genes in this cluster with a threefold or greater transcriptional response to gamma radiation all showed an increased induction with increasing gamma radiation dosage.

CONCLUSION

Analyzing genome-wide transcriptional changes to multiple doses of external stresses enabled the identification of cellular responses that are modulated by magnitude of the stress, providing insights into how a cell deals with genotoxicity.

Gene expression profile analysis: an emerging approach to investigate mechanisms of genotoxicity

The response to stress triggers transcriptional activation of genes involved in cell survival and/or cell death. Thus, the monitoring of gene expression levels in large gene sets or whole genomes in response to various agents (toxicogenomics) has been proposed as a tool for investigating mechanisms of toxicity. Although standard in vitro genetic toxicity testing provides relatively simple and accurate hazard detection, interpretation of positive findings, i.e., in vitro chromosome aberrations, in terms of relevant risk to humans is difficult, due to the limited insight into the underlying mechanisms. Therefore, the development of experimental approaches capable of differentiating a wide range of genotoxic mechanisms is expected to significantly improve risk assessment. The goal of this review is to summarize current developments in toxicogenomic analysis of genotoxic stress, and to provide a perspective on the application of gene expression profile analysis in genetic toxicology.

Annotation of unknown yeast ORFs by correlation analysis of microarray data and extensive literature searches

Changes in the expression of genes were used to elucidate the metabolic pathways and regulatory mechanisms that respond to environmental or genetic modifications. Results from previously published chemostat datasets were merged with novel data generated in the present study. ORFs displaying significant changes in expression that correlated with those of other ORFs were analysed using GO mapping tools and supplemented by literature information. The strategy developed was used to propose annotations for ORFs of unknown function. The following ORFs were assigned functions as a result of this study: YMR090w, YGL157w, YGR243w, YLR327c, YER121w, YFR017c, YGR067c, YKL187c, YGR236c (SPG1), YMR107w (SPG4), YMR206w, YER067w, YJL103c, YNL175C (NOP13) YJL200C, YDL070C (FMP16) and YGR173W.