Analyzing the dose-dependence of the Saccharomyces cerevisiae global transcriptional response to methyl methanesulfonate and ionizing radiation
© Benton et al; licensee BioMed Central Ltd. 2006
Received: 14 April 2006
Accepted: 01 December 2006
Published: 01 December 2006
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.
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.
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.
DNA is the conduit for the transmission of genetic information between generations of a species . The linear sequence of DNA is labile and can be altered by both intracellular and extracellular sources in a variety of ways. For example, within a cell errors by replication machinery can lead to point mutations while the physical stresses of replication can fracture strands of DNA. External sources can mutate DNA by alkylating single strands, causing interstrand crosslinks, or even breaking strands apart, resulting in the loss of some bases upon ligation [2–4].
One of a cell's most crucial tasks is maintaining the integrity of genetic information from one generation to the next . Although some mutations are beneficial in the long term by allowing for adaptation to changing environments via natural selection, most mutations are detrimental and can lead to consequences as serious as cell death or transformation; the correlation between mutagenic activity and carcinogenic activity is approximately 90% . This high correlation has increased interest in identifying causes of DNA damage and the effects of this damage on the cellular environment with the hopes of both preventing cancer and developing better chemotherapeutic agents for the treatment of existing cancerous cells [7–10].
Saccharomyces cerevisiae is commonly used as a model organism for studying the DNA damage response because its genome is well-annotated and the cells are easily cultured and genetically manipulated. Furthermore, because yeasts are eukaryotic, the gene regulation and biochemical pathways involved in responding to DNA damage are expected to be similar to, although perhaps simpler than, those of higher organisms. Indeed, the basic cellular responses to DNA damage have been shown to be similar in yeast and humans [11, 12]. Many DNA damage sensory and repair mechanisms in humans have been identified from their homology to similar mechanisms in yeast .
The study of the DNA damage response in yeast has uncovered a broad variety of mechanisms used by the cells to sense and repair the damage. When cells sense a defect in their DNA, several signal transduction pathways are induced which in turn activate various checkpoints that cause cell-cycle arrest at the G1, S, or G2 phases until the defect can be repaired . These signal transduction pathways regulate genome-wide expression patterns, causing hundreds of genes to be induced or repressed in direct response to the perceived damage . Many of the induced genes are also upregulated in response to cell cycle arrest in stationary phase or in response to general cell stresses, indicating that DNA damage response pathways in yeast overlap with other cell cycle pathways [16, 17]. Once cell division is arrested, repairs can be made in a variety of ways depending on the type of damage which has occurred. Single-strand lesions can be repaired by base excision repair, nucleotide excision repair, or mismatch repair using the non-altered strand as a template, while double-strand breaks are repaired by homologous recombination or non-homologous end joining [18–22].
To better define the cellular pathways involved in the recognition and repair of DNA damage, it is important to identify the global set of genes induced or repressed in response to the damage. These genes have been identified in a number of ways. Traditionally, pools of mutants were screened for sensitivity to various DNA damaging agents, enabling the proposal of epistasis groups based on common responses [23, 24]. In the last 10 years, the proliferation of spotted and oligonucleotide microarrays have enabled genome-wide profiling in response to multiple DNA damaging agents [15, 25–27]. These types of studies identified hundreds of genes not previously known to be modulated in response to DNA damage. Recently the two approaches have been combined by culture of gene deletion mutant strains tagged with unique identification sequences in the presence of a variety of DNA damaging agents, followed by PCR amplification and microarray hybridization of the identification sequences to assess the relative fitness level of each strain as a result of the exposure [28–31]. With each new study, genes not previously known to be regulated by DNA damage are identified, suggesting that there are additional genes and pathways whose control is crucial in DNA damage response yet to be distinguished.
In this study, we identified genes involved in the transcriptional response to DNA damage based on variation in expression upon exposure to two representative DNA damaging agents, methylmethane sulfonate (MMS) and ionizing gamma radiation (γ-ray), over three orders of magnitude of agent dose. MMS is a DNA alkylating agent known to methylate guanine residues at the N-7 position [30, 32]. While this methylation is often tolerable to the cells, MMS can also methylate guanine at the O-6 position or adenine at the N-3 site [30, 32]. These latter lesions inhibit DNA synthesis and must be repaired for replication to properly occur [30, 32]. Ionizing radiation induces multiple types of DNA damage, but double strand breaks are considered to be the most lethal . MMS generally produces single-strand lesions which inhibit DNA replication at the S phase checkpoint, while double-strand breaks caused by γ-ray inhibit the G2/M checkpoint . Using spotted microarrays, Gasch et al. performed studies where S. cerevisiae cells were cultured in the presence of 0.02% MMS and the transcriptional response was assessed at multiple timepoints up to 120 min . Jelinsky et al. performed a similar work using oligonucleotide arrays to monitor the transcriptional response of S. cerevisiae to 0.1% MMS at three timepoints up to 60 min . Both groups identified clusters of genes co-regulated temporally in response to MMS.
Collectively, studies of the S. cerevisiae transcriptional response to MMS have reported results for multiple MMS doses. Jelinsky et al. reported that 693 genes showed a significant transcriptional modulation in response to 0.1% MMS for one hr , while Gasch et al. reported that more than 750 genes showed a significant change in transcription levels when exposed to 0.02% MMS for two hr . More recently, using oligonucleotide arrays, Caba et al. identified 1,377 genes with a significant transcriptional response to 0.12% MMS for two hr . These results suggest a dose-dependent element to the transcriptional response of S. cerevisiae to MMS, although yeast strains and culture conditions varied in these studies. Jelinsky et al. did examine multiple doses and reported a total of 1,863 genes responding to at least one of three doses; however, these doses varied by less than one order of magnitude .
Similarly, Jelinsky et al. observed 248 genes with a significant transcriptional response to 300 Gy of gamma irradiation . Gasch et al. reported over 500 genes whose expression changed as a result of exposure to 170 Gy of gamma irradiation . De Sanctis et al. used spotted microarrays to compare the global transcriptional responses of wild-type and rad53 strains to ionizing radiation administered in doses of 0.8 Gy, 8 Gy, 80 Gy, 300 Gy, and 400 Gy . They observed 72 genes with a 2-fold or more change in transcription level in response to 80 Gy (with 56 genes being induced and 16 genes being repressed) and 86 genes with a 2-fold or more change in transcription level in response to 400 Gy. Recently, Mercier et al. used spotted microarrays to identify approximately 1400 genes with a significant change in transcription level after exposure to 200 Gy of ionizing radiation . In this work, we used oligonucleotide microarrays to quantify the transcriptional response of S. cerevisiae cells to multiple doses, varying over three orders of magnitude, of MMS and γ-ray to show how the global response varies with dose and to identify sets of genes that respond similarly to dose changes.
Results and Discussion
Hierarchical clustering of array data
Global transcriptional response to MMS and γ-ray
Figure 1 shows an overview of the global transcriptional responses to both MMS and γ-ray, obtained by simultaneous clustering of genes and treatments. Clustering of treatments shows that all MMS experiments grouped together, as did the γ-ray-treated samples. Furthermore, within the cluster for a specific type of damage, the highest doses tested were clustered together (0.01% and 0.1% MMS and 10 Gy and 100 Gy for γ-ray). Clearly, exposure to MMS and γ-ray evokes a widespread transcriptional response in S. cerevisiae cells. Many of the clustered genes exhibited a transcriptional response to both MMS and γ-ray, although, it is clear that the response often varied with dose, as expected (see Additional Files 1, 2, 3, 4).
For further clarity, a Venn diagram (Figure 2) is presented to illustrate the number of genes modulated by at least 3-fold to each DNA damaging agent (2-A), each dose of MMS (2-B), and each dose of γ-ray (2-C).
Clusters responsive to MMS and γ-ray
In addition to merely identifying genes that were modulated in response to MMS and/or γ-ray, we sought to identify clusters of genes that exhibited similar or opposite dose dependent responses to both MMS and γ-ray. Clustering diagrams of four clusters (A-D) highlighted in Figure 1 are shown in Figure 3.
Cluster 3-A contains 3 genes, all of which have been shown to play a role in mother-daughter cell separation and all of which were down-regulated in response to all doses of MMS and γ-ray. One of these genes is CTS1, a chitinase required for cell separation after mitosis, whose transcription is induced during the late M and early G1 stages of the cell cycle [36, 37]. Microarray analysis of S. cerevisiae exposed to 5-fluorocytosine, an inhibitor of DNA synthesis, also identified down-regulation of CTS1 . DSE2 is involved in the separation of daughter cells from mother cells by degrading the cell wall from the daughter side [39, 40]. ASH1 localizes in daughter cells during mitosis and inhibits the transcription of HO, preventing the daughter cells from changing mating types . The fact that the transcription factor ASH1 was downregulated after DNA damage exposure correlates well with observations of Workman et al. that Ash1p binds more genes in the absence of MMS than when MMS is present . The co-regulation of these 3 genes implies that the cell may be actively trying to stop transmission of damaged DNA to daughter cells by inhibiting cell division.
Cluster 3-B contains 2 genes repressed in response to MMS and induced by γ-ray exposure. FAS1 is involved in the synthesis of fatty acids [43, 44]. A third gene with similar behavior but excluded from this cluster due to higher than acceptable FDR values, HMG1, catalyzes the rate-limiting step in sterol biosynthesis . Gasch et al. observed the repression of genes responsible for ergosterol synthesis in response to MMS , while De Sanctis et al. noted many of these same genes were induced in response to low doses of ionizing radiation, approximately 80 Gy or less .
Cluster 3-C contains 2 genes induced in response to MMS and repressed in response to γ-ray. One of the genes in this cluster, YDL016C, is a hypothetical ORF. The other gene, CHA4, is a transcription factor for CHA1. In the presence of serine or threonine as the sole nitrogen source, CHA1 is actively transcribed, but its expression is not detectable when cells are grown in the presence of other nitrogen sources such as ammonium . A recent study by Yang et al. suggests CHA4 may play a role in cell-cycle regulation  as one of four putative interaction partners of the two forkhead transcription factors, FKH1 and FKH2, which are involved in cell-cycle regulation and mitotic exit and seem to play a role in yeast dimorphism . This is particularly interesting since 5 of the genes clustered in Figure 1 (ASH1, BEM3, CDC24, TEC1, and UBA4) play a role in regulating pseudohyphal growth in S. cerevisiae.
The cluster highlighted in Figure 3-D is composed of 54 genes whose transcription is induced upon exposure to all doses of MMS and γ-ray used in this study. Several of the genes in this cluster are known to be involved in cellular responses to stress or DNA damage. GAD1 has been shown to mediate resistance to oxidative stress induced by hydrogen peroxide and diamide , and NTG1 has also been identified as crucial in repairing oxidative stress in DNA . RAD16 is part of the nuclear excision repair pathway and has been shown to be involved in double-strand break repair in the fission yeast, Schizosaccharomyces pombe, suggesting an explanation for its induction in response to γ-ray . NEJ1 is involved in the regulation of non-homologous end joining repair in response to double-strand breaks . A potential reason for NEJ1 induction by MMS may involve creation of single-strand breaks following excision repair of methylated bases. When the replication fork confronts these breaks, double-strand breaks can be caused by the forces of replication , which could elicit a partial transcriptional response in genes responsible for repairing double-strand breaks. This could also address induction of RAD54, known to be crucial in homologous recombination for the repair of double-strand breaks , in response to MMS. The genes shown in this cluster suggest there is some overlap of DNA-damage repair pathways in sensing and responding to damage induced by MMS and γ-ray, as members of the base excision repair pathway, nucleotide excision repair pathway, non-homologous end joining pathway, and homologous recombination pathway were all induced.
To gain additional insight into the coregulation of Cluster 3-D, the MEME system was used to identify highly conserved motifs within the regions 500 base pairs upstream from all of the genes in the cluster . Thirty seven of the genes contained a 16 bp region with an E-value of 0.05 (see Additional file 5). To our knowledge, this sequence (T-T- [CT]- [CA]- [CT]-C- [CT]- [CA]-T- [TC]- [TC]-T- [TC]-C- [TC]- [TA]) has not been shown to be a transcription factor binding site in S. cerevisiae [56, 57]. Although these 37 genes are involved in a wide variety of biological processes, it is possible that this conserved motif comprises or contains a site to which a transcription factor induced by DNA-damage can bind and simultaneously regulate several genes induced by both MMS and γ-ray, but clearly further experiments are required to prove this notion. A comparison of this putative transcription factor binding site with the JASPAR database showed closest homology with a mouse pancreatic development transcriptional repressor (81.78%)  and a human nuclear receptor ligand activated transcription factor (67.67%) .
Genes modulated by only one dose of γ-ray
During analysis of the dose dependent transcriptional response to MMS and γ-ray, we sought to identify genes whose transcription was only modulated in response to a single dose of damaging agent. While all of the genes found to respond to only one dose of MMS responded only to the highest dose, 44 genes responded to only one of the γ-ray doses tested and interestingly some responded only to the lowest or to the intermediate dose (Figure 4).
Expression of 5 genes was induced by a single γ-ray dose, with perhaps the most interesting being UGA1. UGA1 is a γ-aminobutyrate transaminase that along with GAD1, discussed earlier, helps regulate the tolerance of cells to oxidative stresses . Its induction at the low dose of 10 Gy indicates the sensitivity of this particular stress response in S. cerevisiae. Many more genes, 39, were repressed at only one dose of γ-ray. At the lowest dose tested, 16 genes were repressed, including 3 involved in protein synthesis regulation (TFS1, BSD2, and PPH21). This down-regulation is consistent with features of the environmental stress response in yeast . Another interesting gene in this cluster is ECO1, which is required for cohesion of sister chromatids during DNA replication, a necessity for repair of double-strand breaks caused by replication . The fact that 44 genes were induced or repressed at only one dose of MMS or γ-ray illustrates one benefit of analyzing genome-wide transcriptional changes to external stresses at multiple doses.
The MEME system was used to identify conserved motifs within the clusters shown in Figure 4. All 16 of the genes induced by 1 Gy of γ-ray conserved a similar 12 bp sequence ([CG]-T-T-T-T-T-T-T- [TC]- [TC]- [TA]-T) with an E-value of 0.08. This motif contains a reported binding site of Azf1p, TTTTTCTT [56, 57]. Azf1p is a transcription factor shown to regulate cell cycle genes in response to DNA damage caused by MMS . While Azf1p is not known to regulate all of the genes in this cluster, it is possible that Azf1p is involved in regulating these genes in response to DNA damage caused by γ-ray, but further experimentation is necessary to test this relationship. A comparison of this putative transcription factor binding site with the JASPAR database showed closest homology with a human (79.17%)  and a rat (80.45%) forkhead transcription factor .
Global transcriptional response to either MMS or γ-ray
To identify genes that were coregulated at all doses of either MMS or γ-ray, the genes were reclustered such that MMS and γ-ray exposure were considered separately (Figure 5) (see Additional Files 3, 4, and 6). All of the genes in each cluster in Figure 5 had at least one transcriptional fold change of 3 or more. One of the most visually evident differences between Figure 5-A and 5-B is the relative number of genes induced or repressed by MMS and γ-ray. Of the 714 genes modulated by MMS (5-A), 581 (81%) are induced (when the dose yielding the highest modulation is considered). Of the 482 genes induced or repressed by γ-ray, 228 (47%) are induced. This suggests that cells modulate entire cellular processes differently based on the type of damage occurring. For example, MMS causes an increased repression of cellular transport and metabolism genes, while γ-ray causes an increased repression of protein regulation and translation genes (see Additional file 6).
Comparing genes with roles in DNA replication/repair and cell-cycle progression and how they were differentially modulated by MMS and γ-ray provides insight into cell responses to these agents. In the cluster shown in Figure 1, a total of 85 genes fell into the DNA replication/repair and cell-cycle progression category, 4 of which are nucleotide excision repair pathway genes. In the MMS cluster (Figure 5-A), a total of 63 genes fell into this category, and all 4 nucleotide excision repair pathway genes are present. In the γ-ray cluster (Figure 5-B), 59 genes fell into the DNA replication/repair and cell-cycle progression category, but only one of the nucleotide excision repair genes, RAD16, was significantly modulated. This variation suggests that, as would be expected, the nucleotide excision repair pathway is much more active in response to MMS than γ-ray.
An interesting subset of the genes involved in cellular metabolic processes includes 8 genes involved in ergosterol synthesis – ERG28, ERG2, ERG11, ERG6, ERG25, ERG13, ERG1, and MVD1. These genes all appear in the MMS and γ-ray cluster (Figure 1) and all 8 also appear in the MMS cluster (Figure 5-A), but none are present in the γ-ray cluster (Figure 5-B), suggesting a link between the MMS response and ergosterol synthesis. Gasch et al. previously noted four of these eight genes (along with several other ergosterol synthesis genes) were down-regulated in response to 0.02% MMS . They attributed this regulation to the fact that many of these genes are regulated by Hap1p and suggested that Hap1p could be methylated by MMS. Jelinsky and Samson reported the repression of ERG3, a sterol desaturase, and SUR4, a sterol isomerase, upon exposure to 0.01% MMS .
Considering the role of DAP1 offers additional insight into the link between ergosterol synthesis and the MMS response. Hand et al. found that dap1 strains were extremely sensitive to MMS, but still underwent cell cycle arrest similar to wild type cells . They proposed that DAP1 is involved in the response to MMS, but not in triggering cell-cycle arrest. DAP1 is also required for normal ergosterol levels in S. cerevisiae; however dap1 strains did not possess increased membrane permeability despite lower levels of ergosterol . This demonstrated that the MMS sensitivity of dap1 was not due to a compromised membrane allowing more MMS inside the cell, but rather some other pathway interactions within the cell. The down-regulation of ergosterol synthesis genes in response to MMS may allow cells to divert energy from normal ergosterol synthesis levels to fighting the effects of MMS. An alternative explanation could be that the lack of cell division due to cell-cycle arrest has lessened the need for ergosterol. Neither of these explanations, however, explains why these genes are down regulated in response to MMS but not γ-ray.
Clusters with dose-dependent responses to MMS or γ-ray
To illustrate the dose-dependent response of certain genes to either MMS or γ-ray, six smaller clusters from Figure 5 in which one of three specific dose-dependent responses was observed: monotonically increasing, monotonically decreasing, or biphasic (where peak induction or repression occurred at the intermediate dose) were identified and are presented in Figure 6.
Cluster 6-A contains 171 genes which exhibit monotonically increasing expression in response to MMS: the genes are induced at all doses of MMS tested and the level of induction increases with MMS dose. Four of the genes in this cluster are directly involved in sensing or repairing DNA damage – MAG1, HSM3, YIM1, and DDR48. HSM3 is involved in DNA mismatch repair in slowly dividing cells [64, 65], YIM1 is believed to be involved in the cellular response to DNA-damaging agents , and DDR48 is induced in response to heat shock or DNA-lesion producing treatments . Additionally, ten other genes (LCB5, GLR1, GPX2, YAP1, GRE2, SSA4, PRE1, PRE3, UBI4, and YNR064C) involved in stress response were also monotonically induced by escalating MMS concentrations. YAP1 is particularly intriguing because it is considered to be the primary regulator of the oxidative stress response , and it serves as a transcription factor for multiple genes and other transcription factors during the MMS response . PRE1 has recently been shown to be modulated in response to MMS by the transcription factor Swi6p . Another interesting feature of cluster A is that 26 of the genes are involved in protein ubiquitination and/or catabolism. MMS induction of protein degradation has been observed in yeast and higher organisms [27, 69] and may be related to the degradation of methylated proteins.
Cluster 6-B contains 23 genes which exhibit monotonically decreasing repression with increasing MMS concentration. The largest subset within this cluster is a group of 6 genes (PTR2, RPS7A, RPS22B, MRP10, KAP123, and NOG2) involved in protein synthesis or transport. PTR2 has been shown to be down-regulated as part of the stress response initiated by MMS exposure . This correlates well with the observation stated earlier: in response to environmental stress, many genes involved in protein synthesis and metabolism in S. cerevisiae are repressed .
The MEME system was used to identify conserved motifs within cluster 6-B. Thirteen of the genes in this cluster conserved a 14 bp region (T-G-C-G-A-T-G-A-G- [CA]- [TA]- [GT]-A-G) with an E-value of 1.2 E-10. Within this motif is a highly conserved 9 bp region proposed to be a binding site for Cha4p (TGCGATGAG) . However, to date Cha4p has not been shown to regulate any of the genes in this cluster , and CHA4 expression is upregulated at all MMS doses tested in this study (Figure 3-C). A comparison of this putative transcription factor binding site with the JASPAR database showed closest homology with a mouse pancreatic development transcriptional repressor (72.57%) and a human zinc finger transcription factor involved in the differentiation response to Ras in human carcinomas (67.67%) .
Cluster 6-C, a group of four genes with a biphasic induction in response to MMS treatment, is very intriguing. All of the genes in this cluster are involved in the Mec1p kinase pathway known to respond to DNA damage caused by MMS. Ribonucleotide reductase is an enzyme that consists of four subunits, two large (RNR1 and RNR3) and two small (RNR2 and RNR4), and catalyzes the rate-limiting step of dNTP synthesis [72, 73]. An increase in dNTP levels is often associated with DNA damage and is considered necessary for DNA repair and resumption of normal transcription . When a cell senses genotoxicity or a stalled replication fork, the signal is transduced through the kinase checkpoint Mec1p, leading to the activation of the kinases Rad53p and Dun1p . Dun1p plays a role in the phosphorylation and eventual degradation of Sml1p, an inhibitor of RNR gene transcription . When Sml1p is destroyed, RNR genes are transcribed and cellular dNTP levels rise, leading to a resumption of DNA synthesis. HUG1, also transcribed in a kinase checkpoint dependent manner, is downstream from DUN1 and thought to play a role in the recovery from the transcriptional response to DNA damage and cell-cycle arrest . The biphasic induction of DUN1, RNR2, RNR4, and HUG1 was unexpected and may suggest that the Mec1p pathway is eventually scaled back or even shut off at high doses of MMS. RNR1 and RNR3 were also biphasically induced in response to MMS but were excluded from clustering analyses because of prohibitive FDR analysis results.
There are several possibilities for the decreased induction of genes involved in the Mec1p pathway at the highest MMS doses. First, there may be alternative cellular responses to different MMS doses. Gasch et al. reported that mec1 mutants show an induction of RNR2, RNR4, and DUN1 in response to DNA damage, although the response is muted, indicating there may be other mechanisms to derepress the RNR genes . More recently, Dubacq and coworkers described a Snf1p kinase directed pathway activated in response to three DNA-damaging agents that cause stalled replication forks independently of and in parallel to the MEC1 pathway . They described snf1 mutants that were sensitive to hydroxyurea, MMS (0.02% and 0.04%), and cadmium ions but not to UV radiation, γ-ray, hydrogen peroxide, camptothecin, and phleomycin. Based on these results, the activity of the Snf1p kinase mediated pathway in response to higher doses of MMS warrants further study.
A second possibility addressing the biphasic response is that MMS may methylate a member of the Mec1p pathway, leading to its inactivation or destruction. In our experiments, RAD53 showed a monotonically increasing induction in response to MMS treatments above 0.001%, although it was excluded from our clustering analysis as a result of FDR filtering. This fact in conjunction with the cluster shown in Figure 6-C suggests the possibility of a failure in the Mec1p pathway between Rad53 and Dun1. This could have been caused by methylation of a protein with susceptibility to MMS at high doses that was degraded in response to its alkylation by MMS effectively creating a roadblock in the DNA-damage response.
A final possibility, though perhaps less likely, it that the RNR enzyme itself is methylated by MMS, leading to its inability to catalyze dNTP synthesis. Although MMS has been shown to increase dNTP levels at low doses , it is possible that RNR is alkylated at higher doses since RNR has an active cysteine residue that is involved in the catalysis of dNTP synthesis . Once a minimum MMS concentration threshold is passed, RNR could be rendered incapable of catalyzing dNTP synthesis and dNTP levels in the cell would drop. This drop in dNTP concentration could be sensed by the cell and mistakenly interpreted as a completion of DNA synthesis, thus inactivating the Mec1p kinase pathway.
Cluster 6-D contains 5 genes that are induced in a monotonically increasing manner in response to γ-ray. PRM5 has been shown to be induced in response to compromised cell wall integrity . The responses of RNR4, RNR2, and HUG1 shown in this cluster are very interesting. All of these genes were induced at all doses of γ-ray, with the highest induction coming at the highest dose. RNR1 and RNR3 are also achieve maximum induction by γ-ray at the highest dose, but were excluded from this clustering analysis by statistical filtering. This behavior offers an interesting contrast to the biphasic response these genes exhibited upon MMS exposure.
The final two clusters shown in Figure 6 (6-E and 6-F) both show a biphasic response to γ-ray. The three genes in 6-E are repressed, while the five genes in cluster 6-F are biphasically induced by γ-ray. An interesting member of this cluster is RAD1, which functions in the repair of double strand breaks that lack overlapping end sequences .
Approximately 25% of genes responding to either MMS or γ-ray were classified as biological process unknown (Figure 7-A). This number is reasonable considering 34% of S. cerevisiae predicted ORFs are currently unverified . Within this category, 67% of the modulated genes were upregulated.
A second group of genes in Figure 7-A, approximately 9%, was composed of genes involved in DNA synthesis and repair or cell-cycle progression. Sixty-five percent of these genes were upregulated in response to DNA damage. Of the 85 genes classified as DNA synthesis and repair or cell-cycle progression genes, 24 are specifically involved in DNA damage repair, including MAG1 which was highly induced in response to MMS and is known to be a member of the base excision repair pathway involved in the repair of N7-alkylguanines, among the deleterious alkylations caused by MMS [83–87]. RAD16 and RAD7, which encode two parts of a heterotrimeric complex that removes damaged oligonucleotides during nucleotide excision repair , were also induced in response to MMS. RNR2 and RNR4 were highly induced in response to γ-ray; these genes have been shown to be transcriptionally induced in response to hydroxyurea and ionizing radiation . A second subset of the DNA synthesis and repair or cell-cycle progression genes includes 9 genes (6 induced and 3 repressed) that function in cell-cycle checkpoints, such as CDC53, HRT1, and PTK2, all involved in the G1/S transition of the cell cycle. This regulation presumably halts cell cycle progression, preventing the synthesis of damaged DNA.
A third major category of genes whose expression changes upon MMS or γ-ray exposure includes those involved in RNA regulation and transcription. Of the 100 genes in this category, 57 were induced and 43 were repressed. Many of these genes are involved in mRNA catabolism, export from the nucleus, or polyadenylation. For example, HRP1 and CLP1, both induced in response to MMS, are subunits of cleavage factor I, required for the cleavage and polyadenylation of 3' mRNA ends . Other RNA regulation and transcription genes modulated in response to MMS and γ-ray include those involved in the positive and negative RNA polymerase II promoter regulators, spliceosome activity, and ribosome synthesis. The modulation of these genes correlates well with the observation that the transcriptional response of genes involved in cell growth-related processes are regulated by the environmental stress response .
Another category of genes modulated by MMS and γ-ray includes genes involved in protein regulation and translation, accounting for 21% of the genes in Figure 7-A. Seventy-three percent of the genes in this category were induced. In response to environmental stress, expression many genes involved in protein synthesis and metabolism in S. cerevisiae changes, presumably to allow the cell to conserve energy .
Finally, a group of 41 genes involved in the response to various general cell stresses are modulated upon exposure to MMS and/or γ-ray. 38 of these genes were induced by more than 3 fold upon exposure to DNA damage. This group includes genes that respond to drugs, osmotic stress, various metal ions, and oxidative stress, indicating the regulation of many genes involved in general stress response pathways.
Figure 7-B, which shows the functional classification of genes with a 3 fold or greater transcriptional response to MMS only, shows little variation from panel A in the fraction of genes composing each functional category. Changes from panel A to panel C, which shows genes with at least a 3 fold transcriptional response to γ-ray only, were mildly more apparent. Additionally, the percentage of genes involved in DNA replication/repair and cell-cycle progression increased slightly while the fraction of genes involved in protein metabolism and translation and the fraction involved in general cellular metabolism decreased slightly. This change may be due to fewer cellular proteins being damaged by γ-ray than by MMS. This is supported by a comparison of panels B and C, which perhaps indicates that MMS elicits a broader response of genes involved in protein regulation and translation than does γ-ray, presumably from methylation and subsequent degradation of proteins following MMS exposure, leading to a higher rate of protein turnover in the cell.
The global transcriptional response of S. cerevisiae as a function of dose of DNA-damaging agent was measured by exposing S. cerevisiae cells to multiple concentrations, varying over 3 orders of magnitude, of two representative DNA-damaging agents, MMS and γ-ray, and hybridizing their cRNA to oligonucleotide microarrays containing over 6,300 hypothetical and known S. cerevisiae ORFs. Hierarchical clustering of genes with a statistically significant change in transcription enabled the identification of several clusters of coregulated genes whose responses varied with DNA damaging agent dose. Our study showed many genes involved in sensing and repairing DNA damage are regulated by both MMS and γ-ray exposure and suggests interconnection between the pathways that respond to various types of DNA damage. There were differences in the responses as well. For example, 8 genes involved in the ergosterol synthesis pathway were down-regulated in response to MMS, but not significantly modulated by γ-ray exposure, verifying earlier observations [26, 27]. Furthermore, we identified 44 genes which responded to only one dose of γ-ray, highlighting one benefit of measuring the transcriptional response to external stresses over multiple doses. MEME analysis identified conserved motifs upstream of multiple genes in these clusters, potentially signifying common regulatory elements.
Further analysis revealed 6 clusters of genes with a dose-dependent response to either MMS or γ-ray. Clusters were identified based on their mode of response, either monotonically increasing, monotonically decreasing, or biphasic (where peak modulation occurred at the intermediate MMS or γ-ray dose). Analyzing genome-wide transcriptional changes to multiple doses of external stresses enables the identification of cellular responses that are modulated by extent of damage, providing new insights into how a cell deals with genotoxicity. Perhaps the most interesting example of this was a cluster that contained 4 genes (DUN1, RNR2, RNR4, and HUG1) that are members of the Mec1p kinase pathway known to respond to MMS induced damage, presumably due to stalled DNA replication forks. These genes showed a biphasic induction in response to MMS treatment, while many of the same genes, RNR2, RNR4, and HUG1, were monotonically increasingly induced by γ-ray.
Strains and culture conditions
Three independent experiments were performed for each exposure condition. S. cerevisiae strain SPY810 (W303 MAT a ura3-52 his3::hisG leu2::hisG) was grown in rich medium (YPD) . Cultures were inoculated overnight at 30°C and diluted 50 fold with YPD the next morning. Cells were returned to 30°C with shaking at approximately 250 rpm until mid log-phase growth was achieved.
DNA damage induction
A culture of mid log-phase yeast cells (OD600 ~ 0.5) was divided into four equal samples. Fresh MMS (Acros Organics, 99% pure, Fisher Scientific, Pittsburgh, PA, USA) was added to three of the samples to final concentrations of 0.001% (v/v), 0.01%, and 0.1%. These doses allowed for greater than 90% cell survival after one hour of exposure to all doses, as confirmed by flow cytometry quantifying fractions of cells excluding propidium iodide (see Additional file 7). The fourth sample had no MMS added but was otherwise handled identically to the MMS-treated cultures. After addition of MMS, cells were returned to 30°C incubation with light shaking for 60 min.
A second culture of mid log phase yeast cells was divided into four equal samples. Using a Cs137 source (7 Gy/min), three of the aliquots were exposed at room temperature in YPD to ionizing radiation to final doses of 1 Gy, 10 Gy, and 100 Gy. These doses allowed for greater than 90% cell survival one hr following exposure, as confirmed by flow cytometry quantifying fractions of cells excluding propidium iodide (see Additional file 7). The fourth sample was not irradiated but was otherwise handled identically to the irradiated cultures. After irradiation, cells were returned to 30°C incubation with light shaking for 60 min.
For each of the culture conditions previously described, total RNA was recovered from approximately 5 × 107 S. cerevisiae cells using the RNeasy Mini Kit (Qiagen, Chatsworth, CA, USA). The manufacturer's protocol for enzymatic lysis was followed using 250 U of zymolyase (Associates of Cape Cod, Seikagaku America, Falmouth, MA, USA) per recovery column to generate spheroplasts. Once isolated from the spheroplasts, RNA integrity was examined using OD260/280 absorption ratio and agarose gel electrophoresis with ethidium bromide staining. RNA was stored in RNase free TE at -80°C.
Microarray probe preparation
Microarray methods were modified from . Total RNA was converted to double-stranded cDNA using the Superscript Double-Stranded cDNA synthesis kit from Invitrogen (Carlsbad, CA, USA) and a Proligo oligo dT primer (Sigma-Aldrich, St. Louis, MO, USA) containing a T7 RNA polymerase promoter region (5'-GGCCAGTGAATTGTAATACGACTCACTATAGGGAGGCGG-T24-3'). To synthesize the first strand of cDNA, 5 μg denatured total RNA was combined with 1× first strand buffer, 10 mM DTT, 500 μM (ea.) dNTPs, 5 pmol dT primer, and 600 U of SuperScript II reverse transcriptase, then incubated at 42°C for 60 min.
To synthesize double-stranded cDNA, the RNA-cDNA hybrid formed in the previous step was combined with 1× second strand buffer, 200 μM (ea.) dNTPs, 0.07 U/μL DNA ligase, 0.267 U/μL DNA polymerase I, and 0.0133 U/μL RNase H and incubated at 16°C for 2 hr. After 2 hr, the ends were polished with T4 DNA polymerase for 5 min. The reaction was stopped by adding 0.5 M EDTA. An equal volume of 25:24:1 phenol:chloroform:isoamyl alcohol (Fisher Scientific, Pittsburgh, PA, USA) was added and the solution was vortexed. The aqueous phase was removed and double-stranded cDNA was precipitated in ethanol.
The double-stranded cDNA was converted to cRNA by in vitro transcription using the Ambion MEGAscript T7 IVT kit (Austin, TX, USA). Double-stranded cDNA was combined with 7.5 mM ATP and GTP, 5.625 mM CTP and UTP, and 1.875 mM biotinylated UTP and biotinylated CTP (Perkin-Elmer, Wellesley, MA, USA). After mixing, the solution was incubated at 37°C for 5 hr. The cRNA was recovered using a RNeasy Mini Kit (Qiagen, Valencia, CA, USA). After elution of cRNA, the total volume of cRNA solution was 40 μL; 10 μL of 5× fragmentation buffer (Affymetrix, Santa Clara, CA, USA) was added and the solution was incubated at 95°C for 35 min. Fragmented cRNA was stored at -80°C until microarray hybridization.
This study used NimbleGen™ (Madison, WI, USA) maskless photolithographic, oligonucleotide microarrays . The arrays were designed based on the complete S. cerevisiae genome sequence and annotation of 6,322 verified and hypothetical open reading frames (ORFs) obtained from the Saccharomyces Genome Database. A high-density oligonucleotide array was designed for this study containing 15 probe pairs per ORF, where each pair contains a perfect match (PM) oligonucleotide and a mismatch (MM) oligonucleotide. The PM sequence is a 24-mer oligonucleotide probe from the S. cerevisiae genome. MM probes were identical except for substitutions at positions 6 and 12 from the 5' end. To generate MM sequences, bases were systematically substituted by inversion as follows: A → T, C → G, G → C, and T → A. On each array, genes were represented by an identical number of oligos to minimize bias.
Four microarrays were hybridized at a time using the NimbleGen™ hybriwheel. Three arrays contained samples exposed to one of three doses of DNA-damaging agent (either MMS or γ-ray) and one was a control array with an untreated sample. Pre-hybridization solution consisting of 0.1 μg/μL herring sperm DNA, 0.5 μg/μL acetylated BSA, and 1× MES hybridization buffer was heated to 95°C for 5 min and then cooled to 45°C for 5 min. The solution was centrifuged at 12,000 g for 10 min. 400 μL of pre-hybridization solution was added to each array, followed by incubation at 45°C for 15 min. After the incubation, the pre-hybridization solution was removed from the arrays.
Hybridization solution consisting of 10 μg of cRNA, 0.33 pmole CPK6 oligos (NimbleGen™ Systems), 0.1 μg/μL herring sperm DNA, 0.5 μg/μL acetylated BSA, and 1× MES hybridization buffer was heated to 95°C for 5 min and then cooled to 45°C for 5 min. The solution was then centrifuged at 12,000 g for 10 min. 300 μL of hybridization solution was added to each microarray. The array was incubated for 16 hours at 45°C and 5 rpm on a hybriwheel.
After hybridization, the hybridization solution was removed and the arrays were rinsed 3× with 350 μL nonstringent wash buffer (6× SSPE, 0.01% Tween 20) (1× SSPE is 0.18 M NaCl, 10 mM NaH2PO4, and 1 mM EDTA [pH 7.7]) at room temperature. After the final rinse, arrays were incubated with 350 μL stringent wash buffer (0.1 M MES buffer [pH 6.5], 26 mM NaCl, 0.01% Tween 20) at 45°C for 30 min, during which the stringent wash buffer was replaced every 5 min. After the 30 min incubation, the stringent wash buffer was removed and 350 μL of 1× stain solution consisting of 1× stain buffer, 2 mg BSA, and Cy3 (Amersham, Piscataway, NJ, USA) was applied to each array. After the arrays were incubated with stain solution for 25 min at room temperature, they were rinsed with nonstringent wash buffer before being dipped in 1× NimbleGen™ final wash buffer for 30 seconds. Finally, the microarrays were blown dry with high pressure argon.
Microarray data analysis
Hybridized microarrays were scanned using a Gene Pix 4000B scanner (Axon Instruments, Molecular Devices Corporation, Sunnyvale, CA, USA) and signal intensity values for each probe were extracted using NimbleScan software. For each probe pair, the MM signal was subtracted from the corresponding PM signal to reduce the effect of nonspecific hybridization signal when estimating transcript abundance. It has been noted that that this may be unnecessary  but results of a large number of 'spike-in' experiments suggest that subtraction of nonspecific signal is preferable and reduces false negative and false positive rates . The data from the 15 probe pairs specific for each ORF were averaged, excluding any probe pairs where the signal was more than 2 standard deviations from the average of that probe set. For purposes of normalization of the array intensity data, array-specific scaling factors were calculated by assuming a constant mean signal intensity of 1,000 signal units for each array. The signal for each ORF was obtained by multiplying the raw signal intensity by the array-specific scaling factor. The ratio of signal for each ORF from an MMS or γ-ray treated sample was calculated relative to the signal from a paired untreated control sample. Ratios of signal intensities were log2 transformed for hierarchical clustering. All microarray data are available for public download through the Gene Expression Omnibus at http://www.ncbi.nlm.nih.gov/projects/geo/ (Series GSE6018).
Expression changes were estimated from three replicates of each experimental treatment for all 6,322 ORFs represented on the microarrays. False discovery rate (FDR) analysis as described by Storey  was used to filter gene responses for false positives. An FDR cutoff of 0.05 was applied, suggesting 1 false positive in every 20 genes considered to be differentially expressed. Storey's method is appropriate for comparing expression values on genome wide studies, since it allows one to control false positives without being so conservative that a large number of false negatives will be introduced. Genes that met this statistical significance criteria in at least one experimental treatment and changed by at least 3-fold in one experimental treatment were selected for hierarchical clustering (average-linkage clustering using uncentered Pearson correlation coefficients in Cluster ). All files generated from Cluster were viewed using Java Treeview.
Groups of co-regulated genes were identified by clustering patterns of changes in transcript levels across experimental conditions. Changes were measured for each gene as ratios of signal between experimental and control samples. This identified transcripts that show similar relative levels of change compared to the control, but does not capture differences in absolute RNA abundance levels between genes. An alternative clustering of genes based on the actual signal intensity values, either Perfect Match minus Mismatch or Perfect Match only, could be used to identify genes with similar expression levels across conditions.
MEME system motif analysis
The MEME system was used to search for highly conserved regions of DNA upstream from clusters of coregulated genes . DNA sequences containing the 500 bp upstream from each gene represented in a cluster were simultaneously analyzed with the MEME parameters of minimum width of 5 bp, maximum width of 20 bp, and maximum number of motifs to find of 6. Only enriched regions with E-values of 0.10 or less are reported in the text.
The position-specific probability matrices of conserved motifs identified with MEME were compared to known motifs in the JASPAR database . This comparison used a modified Needleman-Wunsch algorithm as previously described . The comparison yielded a raw score (where the maximum score is twice the width of the shorter motif) and a percentage of maximum score achieved. The highest score for human and rat or mouse genes are reported in the text.
This work was supported by an NSF CAREER award (BES-0238680) and a Cancer Research and Prevention Foundation grant to SPP. MGB is supported by an NHGRI training grant to the Genomic Sciences Training Program (5T32HG002760). Microarrays were generously provided by NimbleGen™. The authors would like to thank Sandra Splinter BonDurant for helpful discussions on RNA isolation, cDNA synthesis, and hybridization techniques. John Luecke provided expertise in the NimbleGen™ hybridization protocols and chip preparation. Audrey Gasch contributed valuable guidance on microarray data analysis. George Bell provided valuable guidance on FDR analysis.
- Hershey AD, Chase M: Independent functions of viral protein and nucleic acid in growth of bacteriophage. The Journal of general physiology. 1952, 36 (1): 39-56. 10.1085/jgp.36.1.39.PubMedPubMed CentralView ArticleGoogle Scholar
- Grzesiuk E: The role of mutation frequency decline and SOS repair systems in methyl methanesulfonate mutagenesis. Acta Biochimica Polonica. 1998, POLAND , 45 (2): 523-533.Google Scholar
- Nohmi T, Masumura K: Molecular nature of intrachromosomal deletions and base substitutions induced by environmental mutagens. Environmental and molecular mutagenesis. 2005, United States , 45 (2-3): 150-161. 10.1002/em.20110.Google Scholar
- Sankaranarayanan K, Wassom JS: Ionizing radiation and genetic risks XIV. Potential research directions in the post-genome era based on knowledge of repair of radiation-induced DNA double-strand breaks in mammalian somatic cells and the origin of deletions associated with human genomic disorders. Mutation research. 2005, Netherlands , 578 (1-2): 333-370.Google Scholar
- Friedberg EC, Walker GC, Siede W: DNA Repair and Mutagenesis. 1995, Washington, DC , ASM PressGoogle Scholar
- Brusick DJ: In vitro mutagenesis assays as predictors of chemical carcinogenesis in mammals. Clinical toxicology. 1977, UNITED STATES , 10 (1): 79-109.Google Scholar
- Madhusudan S, Middleton MR: The emerging role of DNA repair proteins as predictive, prognostic and therapeutic targets in cancer. Cancer treatment reviews. 2005, England , 31 (8): 603-617. 10.1016/j.ctrv.2005.09.006.Google Scholar
- Marriott SJ, Semmes OJ: Impact of HTLV-I Tax on cell cycle progression and the cellular DNA damage repair response. Oncogene. 2005, England , 24 (39): 5986-5995. 10.1038/sj.onc.1208976.Google Scholar
- Mullighan CG, Flotho C, Downing JR: Genomic assessment of pediatric acute leukemia. Cancer journal (Sudbury, Mass). 2005, United States , 11 (4): 268-282.Google Scholar
- Seve P, Dumontet C: Chemoresistance in non-small cell lung cancer. CurrMedChemAnticancer Agents. 2005, Netherlands , 5 (1): 73-88. 10.2174/1568011053352604.Google Scholar
- Critchlow SE, Jackson SP: DNA end-joining: from yeast to man. Trends in biochemical sciences. 1998, ENGLAND , 23 (10): 394-398. 10.1016/S0968-0004(98)01284-5.Google Scholar
- Weinert T: Yeast checkpoint controls and relevance to cancer. Cancer surveys. 1997, UNITED STATES , 29: 109-132.Google Scholar
- Nickoloff JA, Hoekstra MF: DNA Damage and Repair: Contemporary Cancer Research. 1998, Totowa, NJ , Humana PressGoogle Scholar
- Elledge SJ: Cell cycle checkpoints: preventing an identity crisis. Science. 1996, UNITED STATES , 274 (5293): 1664-1672. 10.1126/science.274.5293.1664.Google Scholar
- Jelinsky SA, Estep P, Church GM, Samson LD: Regulatory networks revealed by transcriptional profiling of damaged Saccharomyces cerevisiae cells: Rpn4 links base excision repair with proteasomes. Molecular and cellular biology. 2000, UNITED STATES , 20 (21): 8157-8167. 10.1128/MCB.20.21.8157-8167.2000.Google Scholar
- Fry RC, Begley TJ, Samson LD: Genome-wide responses to DNA-damaging agents. Annual Review of Microbiology. 2005, United States , 59: 357-377. 10.1146/annurev.micro.59.031805.133658.Google Scholar
- Gasch AP, Spellman PT, Kao CM, Carmel-Harel O, Eisen MB, Storz G, Botstein D, Brown PO: Genomic expression programs in the response of yeast cells to environmental changes. Molecular biology of the cell. 2000, UNITED STATES , 11 (12): 4241-4257.Google Scholar
- Aylon Y, Kupiec M: DSB repair: the yeast paradigm. DNA Repair (Amst). 2004, Netherlands , 3 (8-9): 797-815. 10.1016/j.dnarep.2004.04.013.Google Scholar
- Doetsch PW, Morey NJ, Swanson RL, Jinks-Robertson S: Yeast base excision repair: interconnections and networks. Progress in nucleic acid research and molecular biology. 2001, United States , 68: 29-39.Google Scholar
- Hoeijmakers JH: Nucleotide excision repair. II: From yeast to mammals. Trends in genetics : TIG. 1993, ENGLAND , 9 (6): 211-217. 10.1016/0168-9525(93)90121-W.Google Scholar
- Marti TM, Kunz C, Fleck O: DNA mismatch repair and mutation avoidance pathways. Journal of cellular physiology. 2002, United States , Wiley-Liss, Inc, 191 (1): 28-41. 10.1002/jcp.10077.Google Scholar
- Pastwa E, Blasiak J: Non-homologous DNA end joining. Acta Biochimica Polonica. 2003, Poland , 50 (4): 891-908.Google Scholar
- Prakash L, Prakash S: Isolation and characterization of MMS-sensitive mutants of Saccharomyces cerevisiae. Genetics. 1977, UNITED STATES , 86 (1): 33-55.Google Scholar
- Snow R: Mutants of yeast sensitive to ultraviolet light. Journal of Bacteriology. 1967, UNITED STATES , 94 (3): 571-575.Google Scholar
- Caba E, Dickinson DA, Warnes GR, Aubrecht J: Differentiating mechanisms of toxicity using global gene expression analysis in Saccharomyces cerevisiae. Mutation research. 2005, Netherlands , 575 (1-2): 34-46.Google Scholar
- Gasch AP, Huang M, Metzner S, Botstein D, Elledge SJ, Brown PO: Genomic expression responses to DNA-damaging agents and the regulatory role of the yeast ATR homolog Mec1p. Molecular biology of the cell. 2001, United States , 12 (10): 2987-3003.Google Scholar
- Jelinsky SA, Samson LD: Global response of Saccharomyces cerevisiae to an alkylating agent. Proceedings of the National Academy of Sciences of the United States of America. 1999, UNITED STATES , 96 (4): 1486-1491. 10.1073/pnas.96.4.1486.Google Scholar
- Giaever G, Chu AM, Ni L, Connelly C, Riles L, Veronneau S, Dow S, Lucau-Danila A, Anderson K, Andre B, Arkin AP, Astromoff A, El-Bakkoury M, Bangham R, Benito R, Brachat S, Campanaro S, Curtiss M, Davis K, Deutschbauer A, Entian KD, Flaherty P, Foury F, Garfinkel DJ, Gerstein M, Gotte D, Guldener U, Hegemann JH, Hempel S, Herman Z, Jaramillo DF, Kelly DE, Kelly SL, Kotter P, LaBonte D, Lamb DC, Lan N, Liang H, Liao H, Liu L, Luo C, Lussier M, Mao R, Menard P, Ooi SL, Revuelta JL, Roberts CJ, Rose M, Ross-Macdonald P, Scherens B, Schimmack G, Shafer B, Shoemaker DD, Sookhai-Mahadeo S, Storms RK, Strathern JN, Valle G, Voet M, Volckaert G, Wang CY, Ward TR, Wilhelmy J, Winzeler EA, Yang Y, Yen G, Youngman E, Yu K, Bussey H, Boeke JD, Snyder M, Philippsen P, Davis RW, Johnston M: Functional profiling of the Saccharomyces cerevisiae genome. Nature. 2002, England , 418 (6896): 387-391. 10.1038/nature00935.Google Scholar
- Hanway D, Chin JK, Xia G, Oshiro G, Winzeler EA, Romesberg FE: Previously uncharacterized genes in the UV- and MMS-induced DNA damage response in yeast. Proceedings of the National Academy of Sciences of the United States of America. 2002, United States , 99 (16): 10605-10610. 10.1073/pnas.152264899.Google Scholar
- Lee W, St Onge RP, Proctor M, Flaherty P, Jordan MI, Arkin AP, Davis RW, Nislow C, Giaever G: Genome-Wide Requirements for Resistance to Functionally Distinct DNA-Damaging Agents. PLoS Genet. 2005, 1 (2): e24-10.1371/journal.pgen.0010024.PubMedPubMed CentralView ArticleGoogle Scholar
- Winzeler EA, Shoemaker DD, Astromoff A, Liang H, Anderson K, Andre B, Bangham R, Benito R, Boeke JD, Bussey H, Chu AM, Connelly C, Davis K, Dietrich F, Dow SW, El Bakkoury M, Foury F, Friend SH, Gentalen E, Giaever G, Hegemann JH, Jones T, Laub M, Liao H, Liebundguth N, Lockhart DJ, Lucau-Danila A, Lussier M, M'Rabet N, Menard P, Mittmann M, Pai C, Rebischung C, Revuelta JL, Riles L, Roberts CJ, Ross-MacDonald P, Scherens B, Snyder M, Sookhai-Mahadeo S, Storms RK, Veronneau S, Voet M, Volckaert G, Ward TR, Wysocki R, Yen GS, Yu K, Zimmermann K, Philippsen P, Johnston M, Davis RW: Functional characterization of the S. cerevisiae genome by gene deletion and parallel analysis. Science. 1999, UNITED STATES , 285 (5429): 901-906. 10.1126/science.285.5429.901.Google Scholar
- Chang M, Bellaoui M, Boone C, Brown GW: A genome-wide screen for methyl methanesulfonate-sensitive mutants reveals genes required for S phase progression in the presence of DNA damage. Proceedings of the National Academy of Sciences of the United States of America. 2002, United States , 99 (26): 16934-16939. 10.1073/pnas.262669299.Google Scholar
- Friedl AA, Beisker W, Hahn K, Eckardt-Schupp F, Kellerer AM: Application of pulsed field gel electrophoresis to determine gamma-ray-induced double-strand breaks in yeast chromosomal molecules. International journal of radiation biology. 1993, ENGLAND , 63 (2): 173-181.Google Scholar
- Mercier G, Berthault N, Touleimat N, Kepes F, Fourel G, Gilson E, Dutreix M: A haploid-specific transcriptional response to irradiation in Saccharomyces cerevisiae. Nucleic acids research. 2005, England , 33 (20): 6635-6643. 10.1093/nar/gki959.Google Scholar
- De Sanctis V, Bertozzi C, Costanzo G, Di Mauro E, Negri R: Cell cycle arrest determines the intensity of the global transcriptional response of Saccharomyces cerevisiae to ionizing radiation. Radiation research. 2001, United States , 156 (4): 379-387. 10.1667/0033-7587(2001)156[0379:CCADTI]2.0.CO;2.Google Scholar
- Kuranda MJ, Robbins PW: Chitinase is required for cell separation during growth of Saccharomyces cerevisiae. Journal of Biological Chemistry. 1991, UNITED STATES , 266 (29): 19758-19767.Google Scholar
- O'Conallain C, Doolin MT, Taggart C, Thornton F, Butler G: Regulated nuclear localisation of the yeast transcription factor Ace2p controls expression of chitinase (CTS1) in Saccharomyces cerevisiae. Molecular & general genetics : MGG. 1999, GERMANY , 262 (2): 275-282. 10.1007/s004380051084.Google Scholar
- Zhang L, Zhang Y, Zhou Y, Zhao Y, Zhou Y, Cheng J: Expression profiling of the response of Saccharomyces cerevisiae to 5-fluorocytosine using a DNA microarray. International journal of antimicrobial agents. 2002, Netherlands , Elsevier Science B.V. and International Society of Chemotherapy, 20 (6): 444-450. 10.1016/S0924-8579(02)00201-7.Google Scholar
- Tadi D, Hasan RN, Bussereau F, Boy-Marcotte E, Jacquet M: Selection of genes repressed by cAMP that are induced by nutritional limitation in Saccharomyces cerevisiae. Yeast (Chichester, West Sussex). 1999, ENGLAND , John Wiley & Sons, Ltd, 15 (16): 1733-1745.Google Scholar
- Ufano S, Pablo ME, Calzada A, del Rey F, Vazquez de Aldana CR: Swm1p subunit of the APC/cyclosome is required for activation of the daughter-specific gene expression program mediated by Ace2p during growth at high temperature in Saccharomyces cerevisiae. Journal of cell science. 2004, England , 117 (Pt 4): 545-557. 10.1242/jcs.00880.Google Scholar
- Bobola N, Jansen RP, Shin TH, Nasmyth K: Asymmetric accumulation of Ash1p in postanaphase nuclei depends on a myosin and restricts yeast mating-type switching to mother cells. Cell. 1996, UNITED STATES , 84 (5): 699-709. 10.1016/S0092-8674(00)81048-X.Google Scholar
- Workman CT, Mak HC, McCuine S, Tagne JB, Agarwal M, Ozier O, Begley TJ, Samson LD, Ideker T: A systems approach to mapping DNA damage response pathways. Science. 2006, 312 (5776): 1054-1059. 10.1126/science.1122088.PubMedPubMed CentralView ArticleGoogle Scholar
- Schweizer M, Roberts LM, Holtke HJ, Takabayashi K, Hollerer E, Hoffmann B, Muller G, Kottig H, Schweizer E: The pentafunctional FAS1 gene of yeast: its nucleotide sequence and order of the catalytic domains. Molecular & general genetics : MGG. 1986, GERMANY, WEST , 203 (3): 479-486. 10.1007/BF00422073.Google Scholar
- Stukey JE, McDonough VM, Martin CE: The OLE1 gene of Saccharomyces cerevisiae encodes the delta 9 fatty acid desaturase and can be functionally replaced by the rat stearoyl-CoA desaturase gene. Journal of Biological Chemistry. 1990, UNITED STATES , 265 (33): 20144-20149.Google Scholar
- Parks LW, Smith SJ, Crowley JH: Biochemical and physiological effects of sterol alterations in yeast--a review. Lipids. 1995, 30 (3): 227-230. 10.1007/BF02537825.PubMedView ArticleGoogle Scholar
- Petersen JG, Kielland-Brandt MC, Nilsson-Tillgren T, Bornaes C, Holmberg S: Molecular genetics of serine and threonine catabolism in Saccharomyces cerevisiae. Genetics. 1988, UNITED STATES , 119 (3): 527-534.Google Scholar
- Yang YL, Suen J, Brynildsen MP, Galbraith SJ, Liao JC: Inferring yeast cell cycle regulators and interactions using transcription factor activities. BMC genomics [computer file]. 2005, England , 6 (1): 90-10.1186/1471-2164-6-90.Google Scholar
- Zhu G, Spellman PT, Volpe T, Brown PO, Botstein D, Davis TN, Futcher B: Two yeast forkhead genes regulate the cell cycle and pseudohyphal growth. Nature. 2000, ENGLAND , 406 (6791): 90-94. 10.1038/35017581.Google Scholar
- Coleman ST, Fang TK, Rovinsky SA, Turano FJ, Moye-Rowley WS: Expression of a glutamate decarboxylase homologue is required for normal oxidative stress tolerance in Saccharomyces cerevisiae. Journal of Biological Chemistry. 2001, UNITED STATES , 276 (1): 244-250. 10.1074/jbc.M007103200.Google Scholar
- Alseth I, Eide L, Pirovano M, Rognes T, Seeberg E, Bjoras M: The Saccharomyces cerevisiae homologues of endonuclease III from Escherichia coli, Ntg1 and Ntg2, are both required for efficient repair of spontaneous and induced oxidative DNA damage in yeast. Molecular and cellular biology. 1999, UNITED STATES , 19 (5): 3779-3787.Google Scholar
- Prudden J, Evans JS, Hussey SP, Deans B, O'Neill P, Thacker J, Humphrey T: Pathway utilization in response to a site-specific DNA double-strand break in fission yeast. The EMBO journal. 2003, England , 22 (6): 1419-1430. 10.1093/emboj/cdg119.Google Scholar
- Ooi SL, Shoemaker DD, Boeke JD: A DNA microarray-based genetic screen for nonhomologous end-joining mutants in Saccharomyces cerevisiae. Science. 2001, United States , 294 (5551): 2552-2556. 10.1126/science.1065672.Google Scholar
- Lundin C, North M, Erixon K, Walters K, Jenssen D, Goldman AS, Helleday T: Methyl methanesulfonate (MMS) produces heat-labile DNA damage but no detectable in vivo DNA double-strand breaks. Nucleic acids research. 2005, England , 33 (12): 3799-3811. 10.1093/nar/gki681.Google Scholar
- Symington LS: Role of RAD52 epistasis group genes in homologous recombination and double-strand break repair. Microbiology and molecular biology reviews : MMBR. 2002, United States , 66 (4): 630-70, table of contents. 10.1128/MMBR.66.4.630-670.2002.Google Scholar
- Bailey TL, Elkan C: Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol. 1994, 2: 28-36.PubMedGoogle Scholar
- Harbison CT, Gordon DB, Lee TI, Rinaldi NJ, Macisaac KD, Danford TW, Hannett NM, Tagne JB, Reynolds DB, Yoo J, Jennings EG, Zeitlinger J, Pokholok DK, Kellis M, Rolfe PA, Takusagawa KT, Lander ES, Gifford DK, Fraenkel E, Young RA: Transcriptional regulatory code of a eukaryotic genome. Nature. 2004, 431 (7004): 99-104. 10.1038/nature02800.PubMedPubMed CentralView ArticleGoogle Scholar
- MacIsaac KD, Wang T, Gordon DB, Gifford DK, Stormo GD, Fraenkel E: An improved map of conserved regulatory sites for Saccharomyces cerevisiae. BMC Bioinformatics. 2006, 7: 113-10.1186/1471-2105-7-113.PubMedPubMed CentralView ArticleGoogle Scholar
- Smith SB, Ee HC, Conners JR, German MS: Paired-homeodomain transcription factor PAX4 acts as a transcriptional repressor in early pancreatic development. Mol Cell Biol. 1999, 19 (12): 8272-8280.PubMedPubMed CentralView ArticleGoogle Scholar
- Okuno M, Arimoto E, Ikenobu Y, Nishihara T, Imagawa M: Dual DNA-binding specificity of peroxisome-proliferator-activated receptor gamma controlled by heterodimer formation with retinoid X receptor alpha. Biochem J. 2001, 353 (Pt 2): 193-198. 10.1042/0264-6021:3530193.PubMedPubMed CentralView ArticleGoogle Scholar
- Sjogren C, Nasmyth K: Sister chromatid cohesion is required for postreplicative double-strand break repair in Saccharomyces cerevisiae. Current biology : CB. 2001, England , 11 (12): 991-995. 10.1016/S0960-9822(01)00271-8.Google Scholar
- Overdier DG, Ye H, Peterson RS, Clevidence DE, Costa RH: The winged helix transcriptional activator HFH-3 is expressed in the distal tubules of embryonic and adult mouse kidney. J Biol Chem. 1997, 272 (21): 13725-13730. 10.1074/jbc.272.21.13725.PubMedView ArticleGoogle Scholar
- Overdier DG, Porcella A, Costa RH: The DNA-binding specificity of the hepatocyte nuclear factor 3/forkhead domain is influenced by amino-acid residues adjacent to the recognition helix. Mol Cell Biol. 1994, 14 (4): 2755-2766.PubMedPubMed CentralView ArticleGoogle Scholar
- Hand RA, Jia N, Bard M, Craven RJ: Saccharomyces cerevisiae Dap1p, a novel DNA damage response protein related to the mammalian membrane-associated progesterone receptor. EukaryotCell. 2003, United States , 2 (2): 306-317. 10.1128/EC.2.2.306-317.2003.Google Scholar
- Fedorova IV, Gracheva LM, Kovaltzova SV, Evstuhina TA, Alekseev SY, Korolev VG: The yeast HSM3 gene acts in one of the mismatch repair pathways. Genetics. 1998, 148 (3): 963-973.PubMedPubMed CentralGoogle Scholar
- Fedorova IV, Kovaltzova SV, Korolev VG: The yeast HSM3 gene is involved in DNA mismatch repair in slowly dividing cells. Genetics. 2000, 154 (1): 495-496.PubMedPubMed CentralGoogle Scholar
- Birrell GW, Brown JA, Wu HI, Giaever G, Chu AM, Davis RW, Brown JM: Transcriptional response of Saccharomyces cerevisiae to DNA-damaging agents does not identify the genes that protect against these agents. Proc Natl Acad Sci U S A. 2002, 99 (13): 8778-8783. 10.1073/pnas.132275199.PubMedPubMed CentralView ArticleGoogle Scholar
- Treger JM, McEntee K: Structure of the DNA damage-inducible gene DDR48 and evidence for its role in mutagenesis in Saccharomyces cerevisiae. Molecular and cellular biology. 1990, UNITED STATES , 10 (6): 3174-3184.Google Scholar
- Rodrigues-Pousada CA, Nevitt T, Menezes R, Azevedo D, Pereira J, Amaral C: Yeast activator proteins and stress response: an overview. FEBS letters. 2004, Netherlands , 567 (1): 80-85. 10.1016/j.febslet.2004.03.119.Google Scholar
- van Laar T, van der Eb AJ, Terleth C: A role for Rad23 proteins in 26S proteasome-dependent protein degradation?. Mutation research. 2002, Netherlands , 499 (1): 53-61.Google Scholar
- Lee TI, Rinaldi NJ, Robert F, Odom DT, Bar-Joseph Z, Gerber GK, Hannett NM, Harbison CT, Thompson CM, Simon I, Zeitlinger J, Jennings EG, Murray HL, Gordon DB, Ren B, Wyrick JJ, Tagne JB, Volkert TL, Fraenkel E, Gifford DK, Young RA: Transcriptional regulatory networks in Saccharomyces cerevisiae. Science. 2002, 298 (5594): 799-804. 10.1126/science.1075090.PubMedView ArticleGoogle Scholar
- Thiagalingam A, De Bustros A, Borges M, Jasti R, Compton D, Diamond L, Mabry M, Ball DW, Baylin SB, Nelkin BD: RREB-1, a novel zinc finger protein, is involved in the differentiation response to Ras in human medullary thyroid carcinomas. Mol Cell Biol. 1996, 16 (10): 5335-5345.PubMedPubMed CentralView ArticleGoogle Scholar
- Elledge SJ, Davis RW: Two genes differentially regulated in the cell cycle and by DNA-damaging agents encode alternative regulatory subunits of ribonucleotide reductase. Genes & development. 1990, UNITED STATES , 4 (5): 740-751.Google Scholar
- Huang M, Elledge SJ: Identification of RNR4, encoding a second essential small subunit of ribonucleotide reductase in Saccharomyces cerevisiae. Molecular and cellular biology. 1997, UNITED STATES , 17 (10): 6105-6113.Google Scholar
- Chabes A, Georgieva B, Domkin V, Zhao X, Rothstein R, Thelander L: Survival of DNA damage in yeast directly depends on increased dNTP levels allowed by relaxed feedback inhibition of ribonucleotide reductase. Cell. 2003, United States , 112 (3): 391-401. 10.1016/S0092-8674(03)00075-8.Google Scholar
- Zhou Z, Elledge SJ: DUN1 encodes a protein kinase that controls the DNA damage response in yeast. Cell. 1993, UNITED STATES , 75 (6): 1119-1127. 10.1016/0092-8674(93)90321-G.Google Scholar
- Zhao X, Rothstein R: The Dun1 checkpoint kinase phosphorylates and regulates the ribonucleotide reductase inhibitor Sml1. Proceedings of the National Academy of Sciences of the United States of America. 2002, United States , 99 (6): 3746-3751. 10.1073/pnas.062502299.Google Scholar
- Basrai MA, Velculescu VE, Kinzler KW, Hieter P: NORF5/HUG1 is a component of the MEC1-mediated checkpoint response to DNA damage and replication arrest in Saccharomyces cerevisiae. Molecular and cellular biology. 1999, UNITED STATES , 19 (10): 7041-7049.Google Scholar
- Dubacq C, Chevalier A, Mann C: The protein kinase Snf1 is required for tolerance to the ribonucleotide reductase inhibitor hydroxyurea. Molecular and cellular biology. 2004, United States , 24 (6): 2560-2572. 10.1128/MCB.24.6.2560-2572.2004.Google Scholar
- Jung US, Levin DE: Genome-wide analysis of gene expression regulated by the yeast cell wall integrity signalling pathway. Mol Microbiol. 1999, 34 (5): 1049-1057. 10.1046/j.1365-2958.1999.01667.x.PubMedView ArticleGoogle Scholar
- Ma JL, Kim EM, Haber JE, Lee SE: Yeast Mre11 and Rad1 proteins define a Ku-independent mechanism to repair double-strand breaks lacking overlapping end sequences. Mol Cell Biol. 2003, 23 (23): 8820-8828. 10.1128/MCB.23.23.8820-8828.2003.PubMedPubMed CentralView ArticleGoogle Scholar
- Gene Ontology C: Creating the gene ontology resource: design and implementation. Genome research. 2001, United States , 11 (8): 1425-1433. 10.1101/gr.180801.Google Scholar
- Hong EL, Balakrishnan R, Christie KR, Costanzo MC, Dwight SS, Engel SR, Fisk DG, Hirschman JE, Livestone MS, Nash R, Park J, Oughtred R, Skrzypek M, Starr B, Theesfeld CL, Andrada R, Binkley G, Dong Q, Lane C, Hitz B, Miyasato S, Schroeder M, Sethuraman A, Weng S, Dolinski K, Botstein D, Cherry JM: Saccharomyces Genome Database. 2006, 2006:Google Scholar
- Chen J, Derfler B, Maskati A, Samson L: Cloning a eukaryotic DNA glycosylase repair gene by the suppression of a DNA repair defect in Escherichia coli. Proc Natl Acad Sci U S A. 1989, 86 (20): 7961-7965. 10.1073/pnas.86.20.7961.PubMedPubMed CentralView ArticleGoogle Scholar
- Chen J, Derfler B, Samson L: Saccharomyces cerevisiae 3-methyladenine DNA glycosylase has homology to the AlkA glycosylase of E. coli and is induced in response to DNA alkylation damage. Embo J. 1990, 9 (13): 4569-4575.PubMedPubMed CentralGoogle Scholar
- Chen J, Samson L: Induction of S.cerevisiae MAG 3-methyladenine DNA glycosylase transcript levels in response to DNA damage. Nucleic Acids Res. 1991, 19 (23): 6427-6432.PubMedPubMed CentralView ArticleGoogle Scholar
- Kunz BA, Henson ES, Karthikeyan R, Kuschak T, McQueen SA, Scott CA, Xiao W: Defects in base excision repair combined with elevated intracellular dCTP levels dramatically reduce mutation induction in yeast by ethyl methanesulfonate and N-methyl-N'-nitro-N-nitrosoguanidine. Environmental and molecular mutagenesis. 1998, UNITED STATES , 32 (2): 173-178. 10.1002/(SICI)1098-2280(1998)32:2<173::AID-EM13>3.0.CO;2-M.Google Scholar
- McHugh PJ, Gill RD, Waters R, Hartley JA: Excision repair of nitrogen mustard-DNA adducts in Saccharomyces cerevisiae. Nucleic acids research. 1999, ENGLAND , 27 (16): 3259-3266. 10.1093/nar/27.16.3259.Google Scholar
- Yu S, Owen-Hughes T, Friedberg EC, Waters R, Reed SH: The yeast Rad7/Rad16/Abf1 complex generates superhelical torsion in DNA that is required for nucleotide excision repair. DNA Repair (Amst). 2004, Netherlands , Elsevier B.V, 3 (3): 277-287. 10.1016/j.dnarep.2003.11.004.Google Scholar
- Mulder KW, Winkler GS, Timmers HT: DNA damage and replication stress induced transcription of RNR genes is dependent on the Ccr4-Not complex. Nucleic acids research. 2005, England , 33 (19): 6384-6392. 10.1093/nar/gki938.Google Scholar
- Kessler MM, Henry MF, Shen E, Zhao J, Gross S, Silver PA, Moore CL: Hrp1, a sequence-specific RNA-binding protein that shuttles between the nucleus and the cytoplasm, is required for mRNA 3'-end formation in yeast. Genes & development. 1997, UNITED STATES , 11 (19): 2545-2556.Google Scholar
- Sherman F: Getting started with yeast. Methods Enzymol. 2002, 350: 3-41.PubMedView ArticleGoogle Scholar
- Nuwaysir EF, Huang W, Albert TJ, Singh J, Nuwaysir K, Pitas A, Richmond T, Gorski T, Berg JP, Ballin J, McCormick M, Norton J, Pollock T, Sumwalt T, Butcher L, Porter D, Molla M, Hall C, Blattner F, Sussman MR, Wallace RL, Cerrina F, Green RD: Gene expression analysis using oligonucleotide arrays produced by maskless photolithography. Genome research. 2002, United States , 12 (11): 1749-1755. 10.1101/gr.362402.Google Scholar
- Zhou Y, Abagyan R: Match-only integral distribution (MOID) algorithm for high-density oligonucleotide array analysis. BMC Bioinformatics. 2002, 3: 3-10.1186/1471-2105-3-3.PubMedPubMed CentralView ArticleGoogle Scholar
- Choe SE, Boutros M, Michelson AM, Church GM, Halfon MS: Preferred analysis methods for Affymetrix GeneChips revealed by a wholly defined control dataset. Genome Biol. 2005, 6 (2): R16-10.1186/gb-2005-6-2-r16.PubMedPubMed CentralView ArticleGoogle Scholar
- Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci U S A. 2003, 100 (16): 9440-9445. 10.1073/pnas.1530509100.PubMedPubMed CentralView ArticleGoogle Scholar
- Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proceedings of the National Academy of Sciences of the United States of America. 1998, UNITED STATES , 95 (25): 14863-14868. 10.1073/pnas.95.25.14863.Google Scholar
- Sandelin A, Alkema W, Engstrom P, Wasserman WW, Lenhard B: JASPAR: an open-access database for eukaryotic transcription factor binding profiles. Nucleic Acids Res. 2004, 32 (Database issue): D91-4. 10.1093/nar/gkh012.PubMedPubMed CentralView ArticleGoogle Scholar
- Sandelin A, Hoglund A, Lenhard B, Wasserman WW: Integrated analysis of yeast regulatory sequences for biologically linked clusters of genes. Funct Integr Genomics. 2003, 3 (3): 125-134. 10.1007/s10142-003-0086-6.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.