Regulatory networks revealed by transcriptional profiling of damaged Saccharomyces cerevisiae cells: Rpn4 links base excision repair with proteasomes

Functional Genomics and a New Era in Radiation Biology and Oncology

Ionizing radiation is a ubiquitous stress to which all life is continuously exposed, and thus complex mechanisms have evolved to regulate cellular responses to radiation, including cell cycle arrest, DNA repair, and programmed cell death. Changes in gene expression shape part of the response to radiation, and have historically provided insight into the underlying mechanisms of that response. However, the advent of microarrays, which can measure expression of all the genes in a cell simultaneously, has transformed the study of gene expression, and is beginning to have an impact on both basic mechanistic and clinical studies. This article provides an overview of concepts in gene expression and microarray technology, and highlights their impacts on the study of radiation biology.

The sphingolipid long-chain base-Pkh1/2-Ypk1/2 signaling pathway regulates eisosome assembly and turnover

Eisosomes are recently described fungal structures that play roles in the organization of the plasma membrane and endocytosis. Their major protein components are Pil1 and Lsp1, and previous studies showed that these proteins are phosphorylated by the sphingolipid long-chain base-activated Pkh1 and Pkh2 protein kinases in vitro. We show that Pkh1 and Pkh2 phosphorylate Pil1 and Lsp1 in vivo to produce species B, and that heat stress, which activates Pkh1 and Pkh2, generates a more highly phosphorylated species, C. Cells with low Pkh activity lack species B and C and contain abnormally organized eisosomes. To verify that Pil1 phosphorylation is essential for correct eisosome organization, phosphorylated serine and threonine residues were identified and changed to alanines. A variant Pil1 protein lacking five phosphorylation sites did not form eisosomes during log phase growth, indicating that phosphorylation is critical for eisosome organization. We also found that eisosomes are dynamic structures and disassemble when the Ypk protein kinases, which are activated by the sphingolipid-Pkh signaling pathway, are inactivated or when the sphingolipid signal is pharmacologically blocked with myriocin. We conclude that eisosome formation and turnover are regulated by the sphingolipid-Pkh1/2-Ypk1/2 signaling pathway. These data and previous data showing that endocytosis is regulated by the sphingolipid-Pkh1/2-Ypk1/2 signaling pathway suggest that Pkh1 and -2 respond to changes in membrane sphingolipids and transmit this information to eisosomes via Pil1 phosphorylation. Eisosomes then control endocytosis to align the composition and function of the plasma membrane to match demand.

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.

Genome-wide analysis identifies MYND-domain protein Mub1 as an essential factor for Rpn4 ubiquitylation

The proteasome homeostasis in Saccharomyces cerevisiae is regulated by a negative feedback circuit in which the Rpn4 transcription factor upregulates the proteasome genes and is rapidly degraded by the proteasome. Previous work has identified Ubr2 and Rad6 as the cognate E3 and E2 enzymes for Rpn4 ubiquitylation. However, our recent attempts to ubiquitylate Rpn4 using purified Ubr2 and Rad6 proteins in a reconstitution system have been unsuccessful, suggesting that an additional factor is required for Rpn4 ubiquitylation. Here, we screened the entire collection of the single-gene-deletion yeast mutants generated by the Saccharomyces Genome Deletion Project and identified the mub1Delta mutant defective in ubiquitin-dependent degradation of Rpn4. An in vitro reconstitution ubiquitylation assay confirms that Mub1 is the missing factor for Rpn4 ubiquitylation. We further show that Mub1 directly interacts with Ubr2 and Rpn4. The MYND domain of Mub1 may play an important role in Rpn4 ubiquitylation. Interestingly, Mub1 itself is a short-lived protein and its degradation is dependent on the Ubr2/Rad6 ubiquitin ligase. Together, these data suggest that Mub1 and Ubr2 cooperate to transfer ubiquitin to Rpn4 from Rad6 and that Mub1 may switch from a partner to a substrate of the Ubr2/Rad6 ubiquitin ligase.

Protein damage by reactive electrophiles: targets and consequences

It has been 60 years since the Millers first described the covalent binding of carcinogens to tissue proteins. Protein covalent binding was gradually overshadowed by the emergence of DNA adduct formation as the dominant paradigm in chemical carcinogenesis but re-emerged in the early 1970s as a critical mechanism of drug and chemical toxicity. Technology limitations hampered the characterization of protein adducts until the emergence of mass spectrometry-based proteomics in the late 1990s. The time since then has seen rapid progress in the characterization of the protein targets of electrophiles and the consequences of protein damage. Recent integration of novel affinity chemistries for electrophile probes, shotgun proteomics methods, and systems modeling tools has led to the identification of hundreds of protein targets of electrophiles in mammalian systems. The technology now exists to map the targets of damage to critical components of signaling pathways and metabolic networks and to understand mechanisms of damage at a systems level. The implementation of sensitive, specific analyses for protein adducts from both xenobiotic-derived and endogenous electrophiles offers a means to link protein damage to clinically relevant health effects of both chemical exposures and disease processes.

Transient transcriptional responses to stress are generated by opposing effects of mRNA production and degradation

The state of the transcriptome reflects a balance between mRNA production and degradation. Yet how these two regulatory arms interact in shaping the kinetics of the transcriptome in response to environmental changes is not known. We subjected yeast to two stresses, one that induces a fast and transient response, and another that triggers a slow enduring response. We then used microarrays following transcriptional arrest to measure genome-wide decay profiles under each condition. We found condition-specific changes in mRNA decay rates and coordination between mRNA production and degradation. In the transient response, most induced genes were surprisingly destabilized, whereas repressed genes were somewhat stabilized, exhibiting counteraction between production and degradation. This strategy can reconcile high steady-state level with short response time among induced genes. In contrast, the stress that induces the slow response displays the more expected behavior, whereby most induced genes are stabilized, and repressed genes are destabilized. Our results show genome-wide interplay between mRNA production and degradation, and that alternative modes of such interplay determine the kinetics of the transcriptome in response to stress.

Trm9-catalyzed tRNA modifications link translation to the DNA damage response

Transcriptional and posttranslational signals are known mechanisms that promote efficient responses to DNA damage. We have identified Saccharomyces cerevisiae tRNA methyltransferase 9 (Trm9) as an enzyme that prevents cell death via translational enhancement of DNA damage response proteins. Trm9 methylates the uridine wobble base of tRNAARG(UCU) and tRNAGLU(UUC). We used computational and molecular approaches to predict that Trm9 enhances the translation of some transcripts overrepresented with specific arginine and glutamic acid codons. We found that translation elongation factor 3 (YEF3) and the ribonucleotide reductase (RNR1 and RNR3) large subunits are overrepresented with specific arginine and glutamic acid codons, and we demonstrated that Trm9 significantly enhances Yef3, Rnr1, and Rnr3 protein levels. Additionally, we identified 425 genes, which included YEF3, RNR1, and RNR3, with a unique codon usage pattern linked to Trm9. We propose that Trm9-specific tRNA modifications enhance codon-specific translation elongation and promote increased levels of key damage response proteins.

Transcriptional networks in S. cerevisiae linked to an accumulation of base excision repair intermediates

Upon exposure to DNA damaging agents, Saccharomyces cerevisiae respond by activating a massive transcriptional program that reflects the fact that "DNA damaging" agents also damage other cellular macromolecules. To identify the transcriptional response that is specific to DNA damage, we have modulated the first two enzymes in the base excision repair (BER) pathway generating yeast strains with varied levels of the repair intermediates, abasic sites or strand breaks. We show that the number of abasic sites is significantly increased when the 3-methyladenine DNA glycosylase (Mag): AP endonuclease (Apn1) ratio is increased and that spontaneous frame shift mutation is considerably elevated when either Mag, or Mag plus Apn1, expression is elevated. Expression profiling identified 633 ORFs with differential expression associated with BER modulation. Analysis of transcriptional networks associated with the accumulation of DNA repair intermediates identifies an enrichment for numerous biological processes. Moreover, most of the BER-activated transcriptional response was independent of the classical yeast environmental stress response (ESR). This study highlights that DNA damage in the form of abasic sites or strand breaks resulting from BER modulation is a trigger for substantial genome-wide change and that this response is largely ESR-independent. Taken together, these results suggest that a branch point exists in the current model for DNA damage-signaled transcription in S. cerevisiae.

The protein degradation response of Saccharomyces cerevisiae to classical DNA-damaging agents

Genome wide experiments indicate both proteasome- and vacuole-mediated protein degradation modulate sensitivity to classical DNA-damaging agents. Here, we show that global protein degradation is significantly increased upon methyl methanesulfonate (MMS) exposure. In addition, global protein degradation is similarly increased upon exposure to 4-nitroquinoline-N-oxide (4NQO) and UV and, to a lesser extent, tert-butyl hydroperoxide. The proteasomal inhibitor MG132 decreases both MMS-induced and 4NQO-induced protein degradation, while addition of the vacuolar inhibitor phenylmethanesulfonyl fluoride does not. The addition of both inhibitors grossly inhibits cell growth upon MMS exposure over and above the growth inhibition induced by MMS alone. The MMS-induced protein degradation response remains unchanged in several ubiquitin-proteasome and vacuolar mutants, presumably because these mutants are not totally deficient in either essential pathway. Furthermore, MMS-induced protein degradation is independent of Mec1, Mag1, Rad23, and Rad6, suggesting that the protein degradation response is not transduced through the classical Mec1 DNA damage response pathway or through repair intermediates generated by the base excision, nucleotide excision, or postreplication-DNA repair pathways. These results identify the regulation of protein degradation as an important factor in the recovery of cells from toxicity induced by classical DNA-damaging agents.

The putative roles of the ubiquitin/proteasome pathway in resistance to anticancer therapy

The ubiquitin/proteasome (UP) pathway plays a significant role in many important biological functions and alterations in this pathway have been shown to contribute to the pathology of many human diseases, including cancer. Proteasome inhibition has been well established as a rational strategy for the treatment of multiple myeloma and is currently under investigation for the treatment of other haematological malignancies and solid tumours. Recent evidence suggests that proteasome inhibition may also sensitise tumour cells to the actions of both conventional chemotherapy and radiotherapy, suggesting that this pathway may modify clinical response to anticancer therapy. However, conflicting evidence exists as to the roles of the UP pathway in resistance to treatment. This review endeavours to discuss such roles.

Toxicogenomic analysis incorporating operon-transcriptional coupling and toxicant concentration-expression response: analysis of MX-treated Salmonella

Background

Deficiencies in microarray technology cause unwanted variation in the hybridization signal, obscuring the true measurements of intracellular transcript levels. Here we describe a general method that can improve microarray analysis of toxicant-exposed cells that uses the intrinsic power of transcriptional coupling and toxicant concentration-expression response data. To illustrate this approach, we characterized changes in global gene expression induced in Salmonella typhimurium TA100 by 3-chloro-4-(dichloromethyl)-5-hydroxy-2(5H)-furanone (MX), the primary mutagen in chlorinated drinking water. We used the co-expression of genes within an operon and the monotonic increases or decreases in gene expression relative to increasing toxicant concentration to augment our identification of differentially expressed genes beyond Bayesian-t analysis.

Results

Operon analysis increased the number of altered genes by 95% from the list identified by a Bayesian t-test of control to the highest concentration of MX. Monotonic analysis added 46% more genes. A functional analysis of the resulting 448 differentially expressed genes yielded functional changes beyond what would be expected from only the mutagenic properties of MX. In addition to gene-expression changes in DNA-damage response, MX induced changes in expression of genes involved in membrane transport and porphyrin metabolism, among other biological processes. The disruption of porphyrin metabolism might be attributable to the structural similarity of MX, which is a chlorinated furanone, to ligands indigenous to the porphyrin metabolism pathway. Interestingly, our results indicate that the lexA regulon in Salmonella, which partially mediates the response to DNA damage, may contain only 60% of the genes present in this regulon in E. coli. In addition, nanH was found to be highly induced by MX and contains a putative lexA regulatory motif in its regulatory region, suggesting that it may be regulated by lexA.

Conclusion

Operon and monotonic analyses improved the determination of differentially expressed genes beyond that of Bayesian-t analysis, showing that MX alters cellular metabolism involving pathways other than DNA damage. Because co-expression of similarly functioning genes also occurs in eukaryotes, this method has general applicability for improving analysis of toxicogenomic data.

Genomic evolution of the proteasome system among hemiascomycetous yeasts

Components of the proteasome-ubiquitin pathway are highly conserved throughout eukaryotic organisms. In S. cerevisiae, the expression of proteasomal genes is subject to concerted control by a transcriptional regulator, Rpn4p, interacting with a highly conserved cis-regulatory element, PACE, located in the upstream regions of these genes. Taking advantage of sequence data accumulated from 15 Hemiascomycetes, we performed an in silico study to address the problem of how this system might have evolved among these species. We found that in all these species the Rpn4p homologues are well conserved in terms of sequence and characteristic domain features. The "PACE patterns" turned out to be nearly identical among the Saccharomyces "sensu stricto" species, whereas in the evolutionary more distant species the putatively functional cis-regulatory motifs revealed deviations from the "canonical" PACE nonamere sequence in one or two nucleotides. Our findings suggest that during evolution of the Hemiascomycetes such slightly divergent ancestral motifs have converged into a unique PACE element for the majority of the proteasomal genes within the most recent species of this class. Likewise, the Rpn4 factors within the most recent species of this class show a higher degree of similarity in sequence than their ancestral counterparts. By contrast, we did not detect PACE-like motifs among the proteasomal genes in other eukaryotes, such as S. pombe, several filamentous fungi, A. thaliana, or humans, leaving the interesting question which type of concerted regulation of the proteasome system has developed in species other than the Hemiascomycetes.

Top1- and Top2-mediated topological transitions at replication forks ensure fork progression and stability and prevent DNA damage checkpoint activation

DNA topoisomerases solve topological problems during chromosome metabolism. We investigated where and when Top1 and Top2 are recruited on replicating chromosomes and how their inactivation affects fork integrity and DNA damage checkpoint activation. We show that, in the context of replicating chromatin, Top1 and Top2 act within a 600-base-pair (bp) region spanning the moving forks. Top2 exhibits additional S-phase clusters at specific intergenic loci, mostly containing promoters. TOP1 ablation does not affect fork progression and stability and does not cause activation of the Rad53 checkpoint kinase. top2 mutants accumulate sister chromatid junctions in S phase without affecting fork progression and activate Rad53 at the M-G1 transition. top1 top2 double mutants exhibit fork block and processing and phosphorylation of Rad53 and gamma H2A in S phase. The exonuclease Exo1 influences fork processing and DNA damage checkpoint activation in top1 top2 mutants. Our data are consistent with a coordinated action of Top1 and Top2 in counteracting the accumulation of torsional stress and sister chromatid entanglement at replication forks, thus preventing the diffusion of topological changes along large chromosomal regions. A failure in resolving fork-related topological constrains during S phase may therefore result in abnormal chromosome transitions, DNA damage checkpoint activation, and chromosome breakage during segregation.

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.

A proteasome for all occasions

In the ubiquitin-proteasome system, substrates fated for destruction first acquire covalent modification by ubiquitin, and are subsequently destroyed by the proteasome. Traditionally, 26S proteasomes have been seen as largely uniform in their composition and functional capacity. Accordingly, cells can control proteasome abundance via transcriptional pathways that mediate concerted regulation of all known proteasome genes. However, recent evidence suggests that the proteasome is also subject to subunit-specific modes of regulation, which serve to alter proteasome function and may generate ensembles of compositionally distinct proteasomes. These modes of proteasome regulation provide varied means to adapt protein degradation pathways to changing conditions in the cell.

Biting the hand that feeds: Rpn4-dependent feedback regulation of proteasome function

The 26S proteasome of eukaryotic cells mediates ubiquitin-dependent as well as ubiquitin-independent degradation of proteins in many regulatory processes as well as in protein quality control. The proteasome itself is a dynamic complex with varying compositions and interaction partners. Studies in Saccharomyces cerevisiae have revealed that expression of proteasome subunit genes is coordinately controlled by the Rpn4 transcriptional activator. The cellular level of Rpn4 itself is subject to a complex regulation, which, aside of a transcriptional control of its gene, intriguingly involves ubiquitin-dependent as well as ubiquitin-independent control of its stability by the proteasome. A novel study by Ju et al. [D. Ju, H. Yu, X. Wang, Y. Xie, Ubiquitin-mediated degradation of Rpn4 is controlled by a phosphorylation-dependent ubiquitylation signal, Biochim. Biophys. Acta (in press), doi:10.1016/j.bbamcr.2007.04.012] now revealed another level of complexity by showing that phosphorylation of a specific serine residue in Rpn4 is required for its efficient targeting by the Ubr2 ubiquitin ligase.

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

Exposure to DNA-damaging agents can activate cell cycle checkpoint and DNA repair processes to ensure genetic integrity. Such exposures also can affect the transcription of many genes required for these processes. In the budding yeast Saccharomyces cerevisiae, changes of global gene expression as a result of a DNA-damaging agent were previously identified by using DNA chip technology. DNA microarray analysis is a powerful tool for identifying genes whose expressions are changed in response to environmental changes. Transcriptional levels, however, do not necessarily reflect cellular protein levels. Green fluorescent protein (GFP) has been widely used as a reporter of gene expression and subcellular protein localization. We have used 4156 yeast strains expressing full-length, chromosome-tagged GFP fusion proteins to monitor changes of protein levels in response to the DNA-damaging agent, methyl methanesulphonate (MMS). Through flow cytometry, we identified 157 proteins whose levels were increased at least three-fold following treatment with MMS. Of 157 responsible genes, transcriptions of 57 were previously not known to be induced by MMS. Immunoblot experiments with tandem affinity-tagged yeast strains under the same experimental conditions confirmed these newly found proteins as inducible. These results suggest, therefore, that the 57 protein expressions are regulated by different mechanisms, such as post-translational modifications, and not by transcriptional regulation.

Expression of the phospholipid scramblase (PLSCR) gene family during the acute phase response

Phospholipid scramblase 1 (PLSCR1) is a member of PLSCR gene family that has been implicated in multiple cellular processes including movement of phospholipids, gene regulation, immuno-activation, and cell proliferation/apoptosis. In the present study, we identified PLSCR1 as a positive intracellular acute phase protein that is upregulated by LPS in liver, heart, and adipose tissue, but not skeletal muscle. LPS administration resulted in a marked increase in PLSCR1 mRNA and protein levels in the liver. This stimulation occurred rapidly (within 2 h), and was very sensitive to LPS (half-maximal response at 0.1 microg/mouse). Moreover, two other APR-inducers, zymosan and turpentine, also produced significant increases in PLSCR1 mRNA and protein levels, indicating that PLSCR1 was stimulated in a number of models of the APR. To determine signaling pathways by which LPS stimulated PLSCR1, we examined the effect of proinflammatory cytokines in vitro and in vivo. TNFalpha, IL-1beta, and IL-6 all stimulated PLSCR1 in cultured Hep B3 hepatocytes, whereas only TNFalpha stimulated PLSCR1 in cultured 3T3-L1 adipocytes, suggesting cell type-specific effects of cytokines. Furthermore, the LPS-stimulated increase in liver PLSCR1 mRNA was greatly attenuated by 80% in TNFalpha and IL-1beta receptor null mice as compared to wild-type controls. In contrast, PLSCR1 levels in adipose tissue were induced to a similar extent in TNFalpha and IL-1beta receptor null mice and controls. These results indicate that maximal stimulation of PLSCR1 by LPS in liver required TNFalpha and/or IL-1beta, whereas the stimulation of PLSCR1 in adipose tissue is not dependent on TNFalpha and/or IL-1beta. These data provide evidence that PLSCR1 is a positive intracellular acute phase protein with a tissue-specific mechanism for up-regulation.

Ubiquitin-mediated degradation of Rpn4 is controlled by a phosphorylation-dependent ubiquitylation signal

A ubiquitylation signal of a protein substrate is defined as a short primary sequence or a structural feature recognized by a specific E3. Our previous work has mapped the ubiquitylation signal of Rpn4, the transcription activator for the Saccharomyces cerevisiae proteasome genes, to an N-terminal acidic domain (NAD) consisting of amino acids 211-229. However, the molecular mechanism by which Ubr2, the cognate E3, recognizes NAD remains unclear. Here we show that phosphorylation of either Ser-214 or Ser-220 enhances the binding of NAD to Ubr2. However, phosphorylation of Ser-220 but not Ser-214 plays a predominant role in Rpn4 ubiquitylation and degradation. Interestingly, NAD does not constitute the major Ubr2-binding site of Rpn4 even though it serves as the ubiquitylation signal essential for Rpn4 degradation. By contrast, the stable binding with Ubr2 conferred by other domains of Rpn4 is not required for Rpn4 degradation. Our results indicate that ubiquitin-mediated degradation of Rpn4 is controlled by a phosphorylation-dependent ubiquitylation signal. This study also suggests that binding to E3 may be only a part of the function of a ubiquitylation signal.

Quantification of N-(deoxyguanosin-8-yl)-4-aminobiphenyl adducts in human lymphoblastoid TK6 cells dosed with N-hydroxy-4-acetylaminobiphenyl and their relationship to mutation, toxicity, and gene expression profiling

Gene expression profiles that are anchored to phenotypic endpoints may lead to the identification of signatures that predict mutagenicity or carcinogenicity. The study presented here describes the analysis of DNA adducts in the human TK6 lymphoblastoid cell line after exposure to N-hydroxy-4-aminobiphenyl, a mutagenic metabolite of 4-aminobiphenyl. A validated nano-LC microelectrospray mass spectrometry assay is reported for the detection and quantification of N-(deoxyguanosin-8-yl)-4-aminobiphenyl (dG-C8-ABP), the principal DNA adduct of 4-aminobiphenyl. Limits of quantification, based on a signal-to-noise ratio of 10:1, are determined to correspond to approximately 27 fg of dG-C8-ABP injected on-column. The assay has been used to measure the steady-state levels of the adduct in the human TK6 lymphoblastoid cell line as a function of dose (0.5, 1.0, and 10.0 microM) and time (2, 6, and 27 h) after exposure to N-hydroxy-4-aminobiphenyl. The levels of dG-C8-ABP adducts in the cells, ranging from 18 to 500 adducts in 10(9) nucleotides, were then correlated to cell toxicity, induced mutation at the TK (thymidine kinase) and HPRT loci, and gene expression profiling through microarray analysis. Cell cultures were evaluated for toxicity by growth curve extrapolation, mutation assays were performed on the HPRT and TK loci, and gene expression profiles were generated by analyses using microarray technology. In the mutation assay analysis, as the toxicant concentration increased, there was an increase in mutation fraction, indicating a direct correlation to metabolite dosing level and mutations occurring at these two loci. Statistical analysis of the gene expression data determined that a total of 2250 genes exhibited statistically significant changes in expression after treatment with N-OH-AABP (P < 0.05). Among the genes identified, 2245 were up-regulated, whereas 5 genes that had functions in cell survival and cell growth and, hence, could be indicators of toxicity, were down-regulated relative to controls. The results demonstrate the value of anchoring gene expression patterns to phenotypic markers, such as DNA adduct levels, toxicity, and mutagenicity.

Genome-wide transcriptional plasticity underlies cellular adaptation to novel challenge

By recruiting the essential HIS3 gene to the GAL regulatory system and switching to a repressing glucose medium, we confronted yeast cells with a novel challenge they had not encountered before along their history in evolution.

Adaptation to this challenge involved a global transcriptional response of a sizeable fraction of the genome, which relaxed on the time scale of the population adaptation, of order of 10 generations.

For a large fraction of the responding genes there is no simple biological interpretation, connecting them to the specific cellular demands imposed by the novel challenge.

Strikingly, repeating the experiment did not reproduce similar transcription patterns neither in the transient phase nor in the adapted state in glucose.

These results suggest that physiological selection operates on the new metabolic configurations generated by the non-specific large scale transcriptional response to eventually stabilize an adaptive state.

Cells adjust their transcriptional state to accommodate environmental and genetic perturbations. Some common perturbations, such as changes in nutrient composition, elicit well-characterized transcriptional responses that can be understood by simple engineering-like design principles as satisfying specific demands imposed by the perturbation. However, cells also have the ability to adapt to novel and unforeseen challenges. This ability is central in realizing the evolvability potential of cells as they respond to dramatic genetic or environmental changes along evolution. Little is known about the mechanisms underlying such adaptations to novel challenges; in particular, the role of the transcriptional regulatory network in such adaptations has not been characterized. Genome-wide measurements have revealed that, in many cases, perturbations lead to a global transcriptional response involving a sizeable fraction of the genome (Gasch et al, 2000; Jelinsky et al, 2000; Causton et al, 2001; Ideker et al, 2001; Lai et al, 2005). Such global behavior suggests that general collective properties of the genetic network, rather than specific pre-designed pathways, determine an important part of the transcriptional response. It is not known however what fraction of genes within such massive transcriptional responses is essential to the specific cellular demands. It is also unknown whether the non-pre-designed part of the response can have a functional role in adaptation to novel challenges.

To study these questions, we confronted yeast cells with a novel challenge they had not encountered before along their history in evolution. A strain of the yeast Saccharomyces cerevisiae was engineered to recruit the gene HIS3, an essential enzyme from the histidine biosynthesis pathway (Hinnebusch, 1992), to the GAL regulatory system, responsible for galactose utilization (Stolovicki et al, 2006). The GAL system is known to be strongly repressed when the cells are exposed to glucose. Therefore, upon switching to a medium containing glucose and lacking histidine, the GAL system and with it HIS3 are highly repressed immediately following the switch and the cells encounter a severe challenge. We have recently shown that a cell population carrying this rewired genome can adapt to grow competitively in a chemostat in a medium containing pure glucose (Stolovicki et al, 2006). This adaptation occurred on a timescale of ∼10 generations; applying a stronger environmental pressure in the form of a competitive inhibitor to HIS3 (3AT) resulted in a similar adaptation albeit with a longer timescale. Figure 1 shows the dynamics of the population's cell density (blue lines, measured by OD) following a medium switch from galactose to glucose in the chemostat without (A) and with (B) 3AT. The experiments revealed that adaptation occurs on physiological timescales (much shorter than required by spontaneous random mutations), but the mechanisms underlying this adaptation have remained unclear (Stolovicki et al, 2006).

Yeast cells had not encountered recruitment of HIS3 to the GAL system along their evolutionary history, and their genome could not possibly have been selected to specifically address glucose repression of HIS3. This experiment, therefore, provides a unique opportunity to characterize the spontaneous transcriptional response during adaptation to a novel challenge and to assess the functional role of the regulatory system in this adaptation. We used DNA microarrays to measure the genome-wide expression levels at time points along the adaptation process, with and without 3AT. These measurements revealed that a sizeable fraction of the genome responded by induction or repression to the switch into glucose. Superimposed on the OD traces, Figure 1 shows the results of a clustering analysis of the expression of genes as measured by the arrays along time in the experiments. This analysis revealed two dominant clusters, each containing hundreds of genes in each experiment, which responded to the medium switch to glucose by a strong transient induction or repression followed by relaxation to steady state on the timescale of the adaptation process, ∼ 10 generations. The two clusters in each experiment show similar but opposite dynamics.

A detailed analysis of the gene content in the two clusters revealed that only a small portion of the response was induced by a change in carbon source (15% overlap between the corresponding clusters in the two experiments, with and without 3AT). Moreover, it revealed a very low overlap with the universal stress response observed for a wide range of environmental stresses (Gasch et al, 2000; Causton et al, 2001) and with the typical response to amino-acid starvation (Natarajan et al, 2001). Additionally, all known specific responses to stress in the literature are characterized by transient induction or repression with relaxation to steady state within a generation time (Gasch et al, 2000; Koerkamp et al, 2002; Wu et al, 2004), whereas in our experiments relaxation of the transcriptional response occurs over many generations. Taken together, these results show that the transcriptional response observed here is neither a metabolic response to the change in carbon source nor is it a standard response to stress or amino-acid starvation. This raises the possibility that it is a spontaneous collective response that is largely composed of genes that do not have a specific function. This possibility was tested directly by repeating the experiment with different populations and comparing their responses. This procedure revealed reproducible adaptation dynamics and steady states in terms of population density, but showed significantly different transcriptional transient responses and steady states for the two repeated experiments. Thus, a significant portion of the genes that changed their expression during the adaptation process do not have a well-defined and reproducible function in the challenging environment.

The application of a stronger environmental pressure in the form of 3AT had a dramatic effect on the global characteristics of the transcriptional response: it induced a markedly higher correlation among the hundreds of responding genes. Figure 3A compares the array data in color code for the two experiments. It is seen that the emergent pattern of transcription exhibits a higher degree of order by the introduction of high external pressure. Observation of the transcriptional patterns for specific metabolic pathways illustrates the different contributions to the correlated dynamics (Figure 3B–D). A general energetic module such as glycolysis exhibited similar patterns of induction and relaxation in experiments with and without 3AT (Figure 3B). However, in general, we found that more than one-third of the known metabolic modules (30 out of 88 modules described in KEGG) exhibited high expression correlation among their genes when the environmental pressure was high but not when it was low. As an example, Figure 3C shows the histidine biosynthesis pathway and Figure 3D the purine pathway. Note the highly ordered trajectories in the lower panels (with 3AT) compared to the disordered ones in the upper panels (no 3AT). This order extends also between genes belonging to different and even distant metabolic modules. It indicates that a global transcriptional regulatory mechanism is in operation, rather than a local specific one. Surprisingly, genes belonging to the same metabolic pathway exhibited simultaneous positively and negatively correlated dynamics. Thus, an important conclusion of this work is that the global transcriptional response to a novel challenge cannot be explained by a simple cellular or metabolic logic. This is to be expected if the response had not been specifically selected in evolution and was not pre-designed for the challenge.

Our data clearly reveal that the massive transcriptional response underlies the adaptation process to a novel challenge. The novelty of the challenge presented to the cells excludes the possibility that this response has been specifically selected toward this challenge. Thus, transcriptional regulation has dynamic properties resulting in a general massive nonspecific response to a novel perturbation. Such a response in turn allows for metabolic rearrangements, which by feeding back on transcription lead to adaptation of the cells to the unforeseen situation. The drastic change in the expression state of the cell opens multiple new metabolic pathways. Physiological selection works then on these multiple metabolic pathways to stabilize an adaptive state that causes relaxation of the perturbed expression pattern. This scenario, involving the creation of a library of possibilities and physiological selection over this library, is compatible with our understanding of a broad class of biological systems, placing the cellular metabolic/regulatory networks on the same footing as the neural or the immune systems (Gerhart and Kirschner, 1997).

Cells adjust their transcriptional state to accommodate environmental and genetic perturbations. An open question is to what extent transcriptional response to perturbations has been specifically selected along evolution. To test the possibility that transcriptional reprogramming does not need to be ‘pre-designed' to lead to an adaptive metabolic state on physiological timescales, we confronted yeast cells with a novel challenge they had not previously encountered. We rewired the genome by recruiting an essential gene, HIS3, from the histidine biosynthesis pathway to a foreign regulatory system, the GAL network responsible for galactose utilization. Switching medium to glucose in a chemostat caused repression of the essential gene and presented the cells with a severe challenge to which they adapted over approximately 10 generations. Using genome-wide expression arrays, we show here that a global transcriptional reprogramming (>1200 genes) underlies the adaptation. A large fraction of the responding genes is nonreproducible in repeated experiments. These results show that a nonspecific transcriptional response reflecting the natural plasticity of the regulatory network supports adaptation of cells to novel challenges.

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.

Functional genomics of the yeast DNA-damage response

Two high-throughput studies of the DNA-damage response in yeast reveal new regulatory pathways and genes involved.

High-throughput approaches are beginning to have an impact on many areas of yeast biology. Two recent studies, using different experimental platforms, provide insight into new pathways involved in the response of yeast to DNA damage.

A novel yeast-based reporter assay system for the sensitive detection of genotoxic agents mediated by a DNA damage-inducible LexA-GAL4 protein

Yeast-based genotoxicity testing systems can sensitively detect DNA damaging agents in the environment. We have developed a novel "indirect" reporter assay system based on a recombinant yeast containing both a sensor and a reporter plasmid. The sensor plasmid contains a gene encoding the artificial transcription factor of the Escherichia coli LexA DNA binding domain fused to the transcriptional activation domain of yeast Gal4p, which is regulated by the DNA damage-inducible RNR2 promoter. The reporter plasmid contains the E. coli lacZ gene with the LexA binding site in the 5'-upstream region, allowing transcriptional activation by the induced LexA-GAL4 protein. The activity of DNA damage-dependent beta-galactosidase (beta-gal) in the "indirect" reporter assay system was compared with that of a current yeast-based "direct" reporter system. The "indirect" system exhibited 1.5- to 5-fold greater beta-gal activity upon induction by alkylating agents or camptothecin. To increase the sensitivity of the new reporter system further, several deletion yeast strains were tested, and enhanced induction of reporter activity was observed in DNA repair-deficient mag1Delta cells. The "indirect" 96-well microtiter plate assay system is a potentially inexpensive and sensitive method for detecting genotoxic activities in a wide range of compounds, and in polluted environmental samples.

Toxicogenomics: a pivotal piece in the puzzle of toxicological research

Toxicogenomics, resulting from the merge of conventional toxicology with functional genomics, being the scientific field studying the complex interactions between the cellular genome, toxic agents in the environment, organ dysfunction and disease state. When an organism is exposed to a toxic agent the cells respond by altering the pattern of gene expression. Genes are transcribed into mRNA, which in turn is translated into proteins that serve in a variety of cellular functions. Toxicogenomics through microarray technology, offers large-scale detection and quantification of mRNA transcripts, related to alterations in mRNA stability or gene regulation. This may prove advantageous in toxicological research. In the present review, the applications of toxicogenomics, especially to mechanistic and predictive toxicology are reported. The limitations arising from the use of this technology are also discussed. Additionally, a brief report of other approaches, using other -omic technologies (proteomics and metabonomics) that overcome limitations and give global information related to toxicity, is included.

The DNA-damage signature in Saccharomyces cerevisiae is associated with single-strand breaks in DNA

Background

Upon exposure to agents that damage DNA, Saccharomyces cerevisiae undergo widespread reprogramming of gene expression. Such a vast response may be due not only to damage to DNA but also damage to proteins, RNA, and lipids. Here the transcriptional response of S. cerevisiae specifically induced by DNA damage was discerned by exposing S. cerevisiae to a panel of three "radiomimetic" enediyne antibiotics (calicheamicin γ₁^I, esperamicin A1 and neocarzinostatin) that bind specifically to DNA and generate varying proportions of single- and double-strand DNA breaks. The genome-wide responses were compared to those induced by the non-selective oxidant γ-radiation.

Results

Given well-controlled exposures that resulted in similar and minimal cell death (~20–25%) across all conditions, the extent of gene expression modulation was markedly different depending on treatment with the enediynes or γ-radiation. Exposure to γ-radiation resulted in more extensive transcriptional changes classified both by the number of genes modulated and the magnitude of change. Common biological responses were identified between the enediynes and γ-radiation, with the induction of DNA repair and stress response genes, and the repression of ribosomal biogenesis genes. Despite these common responses, a fraction of the response induced by gamma radiation was repressed by the enediynes and vise versa, suggesting that the enediyne response is not entirely "radiomimetic." Regression analysis identified 55 transcripts with gene expression induction associated both with double- or single-strand break formation. The S. cerevisiae "DNA damage signature" genes as defined by Gasch et al. [1] were enriched among regulated transcripts associated with single-strand breaks, while genes involved in cell cycle regulation were associated with double-strand breaks.

Conclusion

Dissection of the transcriptional response in yeast that is specifically signaled by DNA strand breaks has identified that single-strand breaks provide the signal for activation of transcripts encoding proteins involved in the DNA damage signature in S. cerevisiae, and double-strand breaks signal changes in cell cycle regulation genes.

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.

Development of RNR3- and RAD54-GUS reporters for testing genotoxicity in Saccharomyces cerevisiae

S. cerevisiae RNR3 and RAD54 gene transcription becomes strongly activated upon DNA damage. This property was used to construct yeast strains in which DNA damage can be monitored by a very sensitive fluorogenic assay in a convenient 96-well microtiter plate format. These strains carried stably integrated fusions of RNR3 or RAD54 promoters to the E. coli beta-glucuronidase GUS gene. GUS activity was measured by fluorogenic detection, a method that greatly increases the precision and sensitivity of the assay. Detection levels were similar to those of real-time quantitative PCR methods and close to the limits of biological response. The two reporters differed in terms of fold-induction, activation kinetics, sensitivity and specificity upon exposure to a variety of genotoxic compounds. While RNR3-GUS showed the fastest response, RAD54-GUS showed the highest sensitivity: similar to previous reported sensitivities for bacterial and eukaryotic genotoxic detection systems. These reporter strains may complement current genotoxicity tests, but they also have the advantages of higher flexibility, requirement for shorter incubation times, and the capability of being fully automated. In addition, the intrinsic features of the system facilitate its easy improvement by genetic manipulating the yeast strain or by introducing mammalian metabolizing enzymes.

The conserved ATPase Get3/Arr4 modulates the activity of membrane-associated proteins in Saccharomyces cerevisiae

The regulation of cellular membrane dynamics is crucial for maintaining proper cell growth and division. The Cdc48-Npl4-Ufd1 complex is required for several regulated membrane-associated processes as part of the ubiquitin-proteasome system, including ER-associated degradation and the control of lipid composition in yeast. In this study we report the results of a genetic screen in Saccharomyces cerevisiae for extragenic suppressors of a temperature-sensitive npl4 allele and the subsequent analysis of one suppressor, GET3/ARR4. The GET3 gene encodes an ATPase with homology to the regulatory component of the bacterial arsenic pump. Mutants of GET3 rescue several phenotypes of the npl4 mutant and transcription of GET3 is coregulated with the proteasome, illustrating a functional relationship between GET3 and NPL4 in the ubiquitin-proteasome system. We have further found that Get3 biochemically interacts with the trans-membrane domain proteins Get1/Mdm39 and Get2/Rmd7 and that Deltaget3 is able to suppress phenotypes of get1 and get2 mutants, including sporulation defects. In combination, our characterization of GET3 genetic and biochemical interactions with NPL4, GET1, and GET2 implicates Get3 in multiple membrane-dependent pathways.

Widespread, but Non-identical, Association of Proteasomal 19 and 20 S Proteins with Yeast Chromatin

It has recently become clear that various aspects of nucleic acid metabolism and the ubiquitin-proteasome pathway intersect in several direct and important ways. To begin to assess the scope of some of these activities in the yeast Saccharomyces cerevisiae, we assessed the physical and functional association of proteasomal proteins from both the 20 S core and 19 S regulatory particles with approximately 6400 yeast genes. Genome-wide chromatin immunoprecipitation analyses revealed that proteasome substituents are associated with the majority of yeast genes. Many of these associations correlated strongly with expression levels and the presence of RNA polymerase II. Although the data support the presence of the intact 26 S proteasome on most genes, several hundred yeast genes were cross-linked to either the 20 or 19 S complex but not both, consistent with some degree of independent function for the proteasomal subcomplexes.

Non-random clustering of stress-related genes during evolution of the S. cerevisiae genome

BACKGROUND

Coordinately regulated genes often physically cluster in eukaryotic genomes, for reasons that remain unclear.

RESULTS

Here we provide evidence that many S. cerevisiae genes induced by starvation and other stresses reside in non-random clusters, where transcription of these genes is repressed in the absence of stress. Most genes essential for growth or for rapid, post-transcriptional responses to stress in cycling cells map between these gene clusters. Genes that are transcriptionally induced by stresses include a large fraction of rapidly evolving paralogues of duplicated genes that arose during an ancient whole genome duplication event. Many of these rapidly evolving paralogues have acquired new or more specialized functions that are less essential for growth. The slowly evolving paralogues of these genes are less likely to be transcriptionally repressed in the absence of stress, and are frequently essential for growth or for rapid stress responses that may require constitutive expression of these genes in cycling cells.

CONCLUSION

Our findings suggest that a fundamental organizing principle during evolution of the S. cerevisiae genome has been clustering of starvation and other stress-induced genes in chromosome regions that are transcriptionally repressed in the absence of stress, from which most genes essential for growth or rapid stress responses have been excluded. Chromatin-mediated repression of many stress-induced genes may have evolved since the whole genome duplication in parallel with functions for proteins encoded by these genes that are incompatible with growth. These functions likely provide fitness effects that escape detection in assays of reproductive capacity routinely employed to assess evolutionary fitness, or to identify genes that confer stress-resistance in cycling cells.

DNA microarray analyses reveal a post-irradiation differential time-dependent gene expression profile in yeast cells exposed to X-rays and γ-rays

Ionizing radiation (IR) is the most enigmatic of genotoxic stress inducers in our environment that has been around from the eons of time. IR is generally considered harmful, and has been the subject of numerous studies, mostly looking at the DNA damaging effects in cells and the repair mechanisms therein. Moreover, few studies have focused on large-scale identification of cellular responses to IR, and to this end, we describe here an initial study on the transcriptional responses of the unicellular genome model, yeast (Saccharomyces cerevisiae strain S288C), by cDNA microarray. The effect of two different IR, X-rays, and gamma (gamma)-rays, was investigated by irradiating the yeast cells cultured in YPD medium with 50 Gy doses of X- and gamma-rays, followed by resuspension of the cells in YPD for time-course experiments. The samples were collected for microarray analysis at 20, 40, and 80 min after irradiation. Microarray analysis revealed a time-course transcriptional profile of changed gene expressions. Up-regulated genes belonged to the functional categories mainly related to cell cycle and DNA processing, cell rescue defense and virulence, protein and cell fate, and metabolism (X- and gamma-rays). Similarly, for X- and gamma-rays, the down-regulated genes belonged to mostly transcription and protein synthesis, cell cycle and DNA processing, control of cellular organization, cell fate, and C-compound and carbohydrate metabolism categories, respectively. This study provides for the first time a snapshot of the genome-wide mRNA expression profiles in X- and gamma-ray post-irradiated yeast cells and comparatively interprets/discusses the changed gene functional categories as effects of these two radiations vis-à-vis their energy levels.

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.

A systems approach to mapping DNA damage response pathways

Failure of cells to respond to DNA damage is a primary event associated with mutagenesis and environmental toxicity. To map the transcriptional network controlling the damage response, we measured genomewide binding locations for 30 damage-related transcription factors (TFs) after exposure of yeast to methyl-methanesulfonate (MMS). The resulting 5272 TF-target interactions revealed extensive changes in the pattern of promoter binding and identified damage-specific binding motifs. As systematic functional validation, we identified interactions for which the target changed expression in wild-type cells in response to MMS but was nonresponsive in cells lacking the TF. Validated interactions were assembled into causal pathway models that provide global hypotheses of how signaling, transcription, and phenotype are integrated after damage.

Parallel induction of ATM-dependent pro- and antiapoptotic signals in response to ionizing radiation in murine lymphoid tissue

The ATM protein kinase, functionally missing in patients with the human genetic disorder ataxia-telangiectasia, is a master regulator of the cellular network induced by DNA double-strand breaks. The ATM gene is also frequently mutated in sporadic cancers of lymphoid origin. Here, we applied a functional genomics approach that combined gene expression profiling and computational promoter analysis to obtain global dissection of the transcriptional response to ionizing radiation in murine lymphoid tissue. Cluster analysis revealed a prominent pattern characterizing dozens of genes whose response to irradiation was Atm-dependent. Computational analysis identified significant enrichment of the binding site signatures of NF-kappaB and p53 among promoters of these genes, pointing to the major role of these two transcription factors in mediating the Atm-dependent transcriptional response in the irradiated lymphoid tissue. Examination of the response showed that pro- and antiapoptotic signals were simultaneously induced, with the proapoptotic pathway mediated by p53 targets, and the prosurvival pathway by NF-kappaB targets. These findings further elucidate the molecular network induced by IR, point to novel putative NF-kappaB targets, and suggest a mechanistic model for cellular balancing between pro- and antiapoptotic signals induced by IR in lymphoid tissues, which has implications for cancer management. The emerging model suggests that restoring the p53-mediated apoptotic arm while blocking the NF-kappaB-mediated prosurvival arm could effectively increase the radiosensitivity of lymphoid tumors.

A stress regulatory network for co‐ordinated activation of proteasome expression mediated by yeast heat shock transcription factor

Heat shock transcription factor (HSF) mediates the transcriptional response of eukaryotic cells to heat, infection and inflammation, pharmacological agents, and other stresses. Although genes encoding heat shock proteins (HSPs) are the best characterized targets of HSF, recent genome-wide localization of Saccharomyces cerevisiae HSF revealed novel HSF targets involved in a wide range of cellular functions. One such target, the RPN4 gene, encodes a transcription factor that directly activates expression of a number of genes encoding proteasome subunits. Here we demonstrate that HSF co-ordinates a feed-forward gene regulatory circuit for RPN4 activation. We show that HSF activates expression of PDR3, encoding a multidrug resistance (MDR) transcription factor that also directly activates RPN4 gene expression. We demonstrate that the HSF binding site (HSE) in the RPN4 promoter is primarily responsible for heat- or methyl methanesulphonate induction of RPN4, with a minor contribution of Pdr3 binding sites (PDREs), while a Yap1 binding site (YRE) is responsible for RPN4 induction in response to oxidative stress. Furthermore, heat-induced expression of Rpn4 protein leads to expression of Rpn4 targets at later stages of heat stress, providing a temporal controlling mechanism for proteasome synthesis upon stress conditions that could result in irreversibly damaged proteins. In addition, the overlapping transcriptional regulatory networks involving HSF, Yap1 and Pdr3 suggest a close linkage between stress responses and pleiotropic drug resistance.

Unique role for the UbL-UbA protein Ddi1 in turnover of SCFUfo1 complexes

SCF complexes are E3 ubiquitin-protein ligases that mediate degradation of regulatory and signaling proteins and control G1/S cell cycle progression by degradation of G1 cyclins and the cyclin-dependent kinase inhibitor, Sic1. Interchangeable F-box proteins bind the core SCF components; each recruits a specific subset of substrates for ubiquitylation. The F-box proteins themselves are rapidly turned over by autoubiquitylation, allowing rapid recycling of SCF complexes. Here we report a role for the UbL-UbA protein Ddi1 in the turnover of the F-box protein, Ufo1. Ufo1 is unique among F-box proteins in having a domain comprising multiple ubiquitin-interacting motifs (UIMs) that mediate its turnover. Deleting the UIMs leads to stabilization of Ufo1 and to cell cycle arrest at G1/S of cells with long buds resembling skp1 mutants. Cells accumulate substrates of other F-box proteins, indicating that the SCF pathway of substrate ubiquitylation is inhibited. Ufo1 interacts with Ddi1 via its UIMs, and Deltaddi1 cells arrest when full-length UFO1 is overexpressed. These results imply a role for the UIMs in turnover of SCF(Ufo1) complexes that is dependent on Ddi1, a novel activity for an UbL-UbA protein.

Fructose-1,6-bisphosphatase mediates cellular responses to DNA damage and aging in Saccharomyces cerevisiae

Response to DNA damage, lack of nutrients and other stress conditions is an essential property of living systems. The coordinate response includes DNA damage repair, activation of alternate biochemical pathways, adjustment of cellular proliferation and cell cycle progression as well as drastic measures like cellular suicide which prevents proliferation of severely damaged cells. Investigating the transcriptional response of Saccharomyces cerevisiae to low doses of the alkylating agent methylmethane sulfonate (MMS) we observed induction of genes involved in glucose metabolism. RT-PCR analysis showed that the expression of the key enzyme in gluconeogenesis fructose-1,6-bisphosphatase (FBP1) was clearly up-regulated by MMS in glucose-rich medium. Interestingly, deletion of FBP1 led to reduced sensitivity to MMS, but not to other DNA-damaging agents, such as 4-NQO or phleomycin. Reintroduction of FBP1 in the knockout restored the wild-type phenotype while overexpression increased MMS sensitivity of wild-type, shortened life span and increased induction of RNR2 after treatment with MMS. Deletion of FBP1 reduced production of reactive oxygen species (ROS) in response to MMS treatment and in untreated aged cells, and increased the amount of cells able to propagate and to form colonies, but had no influence on the genotoxic effect of MMS. Our results indicate that FBP1 influences the connection between DNA damage, aging and oxidative stress through either direct signalling or an intricate adaptation in energy metabolism.

Analysis of cellular responses to aflatoxin B(1) in yeast expressing human cytochrome P450 1A2 using cDNA microarrays

Aflatoxin B1 (AFB(1)) is a potent human hepatotoxin and hepatocarcinogen produced by the mold Aspergillus flavus. In human, AFB(1) is bioactivated by cytochrome P450 (CYP450) enzymes, primarily CYP1A2, to the genotoxic epoxide that forms N(7)-guanine DNA adducts. To characterize the transcriptional responses to genotoxic insults from AFB(1), a strain of Saccharomyces cerevisiae engineered to express human CYP1A2 was exposed to doses of AFB(1) that resulted in minimal lethality, but substantial genotoxicity. Flow cytometric analysis demonstrated a dose and time dependent S phase delay under the same treatment conditions, indicating a checkpoint response to DNA damage. Replicate cDNA microarray analyses of AFB(1) treated cells showed that about 200 genes were significantly affected by the exposure. The genes activated by AFB(1)-treatment included RAD51, DUN1 and other members of the DNA damage response signature reported in a previous study with methylmethane sulfonate and ionizing radiation [A.P. Gasch, M. Huang, S. Metzner, D. Botstein, S.J. Elledge, P.O. Brown, Genomic expression responses to DNA-damaging agents and the regulatory role of the yeast ATR homolog Mec1p, Mol. Biol. Cell 12 (2001) 2987-3003]. However, unlike previous studies using highly cytotoxic doses, environmental stress response genes [A.P. Gasch, P.T. Spellman, C.M. Kao, O. Carmel-Harel, M.B. Eisen, G. Storz, D. Botstein, P.O. Brown, Genomic expression programs in the response of yeast cells to environmental changes, Mol. Biol. Cell 11 (2000) 4241-4257] were largely unaffected by our dosing regimen. About half of the transcripts affected are also known to be cell cycle regulated. The most strongly repressed transcripts were those encoding the histone genes and a group of genes that are cell cycle regulated and peak in M phase and early G1. These include most of the known daughter-specific genes. The rapid and coordinated repression of histones and M/G1-specific transcripts cannot be explained by cell cycle arrest, and suggested that there are additional signaling pathways that directly repress these genes in cells under genotoxic stress.

Genome-wide responses to DNA-damaging agents

Genome-wide studies of mRNA regulation and phenotypic responses have shown that eukaryotic cells mount a robust and multifaceted response upon exposure to DNA-damaging agents. The integration of theses studies over frameworks provided by protein-protein interactions, protein-DNA interactions, and subcellular localization information have led to the identification of networked responses to damage. Taken together, these studies illustrate that cellular protection from DNA and other macromolecular damage involves an intricate network of proteins involved in many different cellular functions, some of them expected (e.g., DNA repair and cell cycle checkpoints) but many of them unexpected (e.g., protein trafficking and degradation). This review highlights many of the studies that detail genome-wide responses to DNA-damaging agents and examines how these datasets have been used to build a systems view of cellular responses to damage.

Assessment and integration of publicly available SAGE, cDNA microarray, and oligonucleotide microarray expression data for global coexpression analyses

Large amounts of gene expression data from several different technologies are becoming available to the scientific community. A common practice is to use these data to calculate global gene coexpression for validation or integration of other "omic" data. To assess the utility of publicly available datasets for this purpose we have analyzed Homo sapiens data from 1202 cDNA microarray experiments, 242 SAGE libraries, and 667 Affymetrix oligonucleotide microarray experiments. The three datasets compared demonstrate significant but low levels of global concordance (rc<0.11). Assessment against Gene Ontology (GO) revealed that all three platforms identify more coexpressed gene pairs with common biological processes than expected by chance. As the Pearson correlation for a gene pair increased it was more likely to be confirmed by GO. The Affymetrix dataset performed best individually with gene pairs of correlation 0.9-1.0 confirmed by GO in 74% of cases. However, in all cases, gene pairs confirmed by multiple platforms were more likely to be confirmed by GO. We show that combining results from different expression platforms increases reliability of coexpression. A comparison with other recently published coexpression studies found similar results in terms of performance against GO but with each method producing distinctly different gene pair lists.

Functional Division of Substrate Processing Cofactors of the Ubiquitin-Selective Cdc48 Chaperone

Ubiquitin-dependent protein degradation usually involves escort factors that target ubiquitylated substrates to the proteasome. A central element in a major escort pathway is Cdc48, a chaperone-like AAA ATPase that collects ubiquitylated substrates via alternative substrate-recruiting cofactors. Cdc48 also associates with Ufd2, an E4 multiubiquitylation enzyme that adds further ubiquitin moieties to preformed ubiquitin conjugates to promote degradation. Here, we show that E4 can be counteracted in vivo by two distinct mechanisms. First, Ufd3, a WD40 repeat protein, directly competes with Ufd2, because both factors utilize the same docking site on Cdc48. Second, Cdc48 also binds Otu1, a deubiquitylation enzyme, which disassembles multiubiquitin chains. Notably, Cdc48 can bind Otu1 and Ufd3 simultaneously, making a cooperation of both inhibitory mechanisms possible. We propose that the balance between the distinct substrate-processing cofactors may determine whether a substrate is multiubiquitylated and routed to the proteasome for degradation or deubiquitylated and/or released for other purposes.

A haploid-specific transcriptional response to irradiation in Saccharomyces cerevisiae

Eukaryotic cells respond to DNA damage by arresting the cell cycle and modulating gene expression to ensure efficient DNA repair. We used global transcriptome analysis to investigate the role of ploidy and mating-type in inducing the response to damage in various Saccharomyces cerevisiae strains. We observed a response to DNA damage specific to haploid strains that seemed to be controlled by chromatin regulatory proteins. Consistent with these microarray data, we found that mating-type factors controlled the chromatin-dependent silencing of a reporter gene. Both these analyses demonstrate the existence of an irradiation-specific response in strains (haploid or diploid) with only one mating-type factor. This response depends on the activities of Hdf1 and Sir2. Overall, our results suggest the existence of a new regulation pathway dependent on mating-type factors, chromatin structure remodeling, Sir2 and Hdf1 and independent of Mec1 kinase.

Identifying the conserved network of cis-regulatory sites of a eukaryotic genome

A major focus of genome research has been to decipher the cis-regulatory code that governs complex transcriptional regulation. We report a computational approach for identifying conserved regulatory motifs of an organism directly from whole genome sequences of several related species without reliance on additional information. We first construct phylogenetic profiles for each promoter, then use a BLAST-like algorithm to efficiently search through the entire profile space of all of the promoters in the genome to identify conserved motifs and the promoters that contain them. Statistical significance is estimated by modified Karlin-Altschul statistics. We applied this approach to the analysis of 3,524 Saccharomyces cerevisiae promoters and identified a highly organized regulatory network involving 3,315 promoters and 296 motifs. This network includes nearly all of the currently known motifs and covers >90% of known transcription factor binding sites. Most of the predicted coregulated gene clusters in the network have additional supporting evidence. Theoretical analysis suggests that our algorithm should be applicable to much larger genomes, such as the human genome, without reaching its statistical limitation.

Proteasome-associated proteins: regulation of a proteolytic machine

The proteasome is a compartmentalized, ATP-dependent protease composed of more than 30 subunits that recognizes and degrades polyubiquitinated substrates. Despite its physiological importance, many aspects of the proteasome's structural organization and regulation remain poorly understood. In addition to the proteins that form the proteasome holocomplex, there is increasing evidence that proteasomal function is affected by a wide variety of associating proteins. A group of ubiquitin-binding proteins assist in delivery of substrates to the proteasome, whereas proteasome-associated ubiquitin ligases and deubiquitinating enzymes may alter the dynamics of ubiquitin chains already associated with the proteasome. Some proteins appear to influence the overall stability of the complex, and still others have the capacity to activate or inhibit the hydrolytic activity of the core particle. The increasing number of interacting proteins identified suggests that proteasomes, as they exist in the cell, are larger and more diverse in composition than previously assumed. Thus, the study of proteasome-associated proteins will lead to new perspectives on the dynamics of this uniquely complex proteolytic machine.

Gene expression responses to DNA damage are altered in human aging and in Werner Syndrome

The accumulation of DNA damage and mutations is considered a major cause of cancer and aging. While it is known that DNA damage can affect changes in gene expression, transcriptional regulation after DNA damage is poorly understood. We characterized the expression of 6912 genes in human primary fibroblasts after exposure to three different kinds of cellular stress that introduces DNA damage: 4-nitroquinoline-1-oxide (4NQO), gamma-irradiation, or UV-irradiation. Each type of stress elicited damage specific gene expression changes of up to 10-fold. A total of 85 genes had similar changes in expression of 3-40-fold after all three kinds of stress. We examined transcription in cells from young and old individuals and from patients with Werner syndrome (WS), a segmental progeroid condition with a high incidence of cancer, and found various age-associated transcriptional changes depending upon the type of cellular stress. Compared to young individuals, both WS and old individuals had similarly aberrant transcriptional responses to gamma- and UV-irradiation, suggesting a role for Werner protein in stress-induced gene expression. Our results suggest that aberrant DNA damage-induced gene regulation may contribute to the aging process and the premature aging in WS.

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

Genotoxic stress triggers a variety of biological responses including the transcriptional activation of genes regulating DNA repair, cell survival and cell death. Genomic approaches, which monitor gene expressions across large numbers of genes, can serve as a powerful tool for exploring mechanisms of toxicity. Here, using five different agents, we investigated whether the analysis of genome-wide expression profiles in Saccharomyces cerevisiae could provide insights into mechanisms of genotoxicity versus cytotoxicity. To differentiate the genotoxic stress-associated expression signatures from that of a general cytotoxic stress, we compared gene expression profiles following the treatment with DNA-reactive (cisplatin, MMS, bleomycin) and DNA non-reactive (ethanol and sodium chloride) compounds. Although each of the tested chemicals produced a distinct gene expression profile, we were able to identify a gene expression signature consisting of a relatively small number of biologically relevant genes capable of differentiating genotoxic and cytotoxic stress. The gene set includes such upregulated genes as HUG1, ECM4 and previously uncharacterized gene, YLR297W in the genotoxic and GAP1, CGR1 in the cytotoxic group. Our results indicate the potential of gene expression profile analysis for elucidating mechanism of action of genotoxic agents.