Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

Open Access Research article

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

Rebecca C Fry12, Michael S DeMott1, Joseph P Cosgrove13, Thomas J Begley4, Leona D Samson12 and Peter C Dedon12*

  • * Corresponding author: Peter C Dedon pcdedon@mit.edu

  • † Equal contributors

Author Affiliations

1 Department of Biological Engineering, Massachusetts Institute of Technology, NE47-277, 77 Massachusetts Avenue, Cambridge, MA 02139, USA

2 Center for Environmental Health Sciences, Massachusetts Institute of Technology, 56-235, 77 Massachusetts Avenue, Cambridge, MA 02139, USA

3 Novartis Institute for Biomedical Research, 250 Massachusetts Ave., 3C-141, Cambridge, MA 02139, USA

4 University of Albany, Gen*NY*Sis Center, Room 211, 1 Discovery Drive, Rensselear, NY 12144-3456, USA

For all author emails, please log on.

BMC Genomics 2006, 7:313  doi:10.1186/1471-2164-7-313


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/7/313


Received:4 August 2006
Accepted:12 December 2006
Published:12 December 2006

© 2006 Fry et al; licensee 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.

Abstract

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 γ1I, 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.

Background

Exposure to DNA damaging agents can cause mutation and cell death and may ultimately lead to disease. Protection from this damage is provided by a host of DNA repair and cell cycle checkpoint proteins that collectively represent numerous pathways to help in the recovery response [2]. In humans, there are approximately 150 DNA repair and cell cycle checkpoint proteins that serve to ensure the repair of damage caused to DNA [3,4] and most of these have functional homologues in S. cerevisiae. Recently, studies have shown that proteins with additional cellular functions beyond DNA repair and cell cycle regulation may ameliorate the toxic effects of agents that cause DNA damage [5-8].

Genome-wide phenotyping studies to identify genes involved in cellular recovery after exposure to DNA alkylating agents, such as methyl methane sulfonate (MMS), indicate that a vast array of cellular processes are required for the recovery of S. cerevisiae. Transcriptional profiling demonstrated that up to 30% of the S. cerevisiae ~6000 genes respond upon exposure to MMS in a time-, agent- and dose-dependent manner [9,10]. Given this vast transcriptional response to damaging agents like MMS and the knowledge that, in addition to DNA, these agents can also damage proteins, RNA and lipids, we set out to identify the transcriptional response specifically caused by damage to DNA.

Here we compared the responses of S. cerevisiae upon exposure to γ-radiation, a non-selective oxidant that attacks DNA (base and sugar), lipids, carbohydrates, proteins and small metabolites in cells, with the response of S. cerevisiae to a panel of enediyne antibiotics that are known to damage DNA and not other cellular molecules (calicheamicin γ1I, esperamicin A1 and neocarzinostatin; structures shown in Figure 1) [11,12]. The enediyne family is a structurally diverse group of DNA-cleaving molecules that undergo reductive activation, presumably by glutathione in vivo, to form a diradical intermediate that binds with high affinity (109 M-1) in the minor groove of DNA and abstract hydrogen atoms from deoxyribose. This "radiomimetic" damage leads to the formation of well-defined proportions of single- and double-stranded DNA lesions unique to each type of enediyne [11,12] the proportions of which are noted in Table 1. With no damage to DNA bases, the enediynes specifically oxidize the deoxyribose moiety to produce either direct strand breaks with various sugar residues or phosphate groups attached to the 3'- or 5'-ends of the breaks, or various unstable oxidized abasic sites [11].

Table 1. Proportions of Single- and Double-Strand Damage Produced by γ-Radiation and Enediynes

thumbnailFigure 1. Chemical structures and relative ratios of double-strand to single-strand lesions generated. Three enediynes were used for treatment in S. cerevisiae including calicheamicin γ1I, esperamicin A1 and neocarzinostatin. Each of these results in varying ratios of double-strand (DS) to single-strand (SS) lesions in DNA as outlined in Table 1.

Expression profiling experiments have been performed using some of these agents independently [13,14]. However, there are no studies that directly compare the responses of these enediynes to each other and to that of γ-radiation under identical conditions. We establish four key findings: (i) under conditions of similar cell survival, exposure to non-selective γ-radiation results in more extensive reprogramming of S. cerevisiae transcription than does exposure to the DNA-selective enediynes; (ii) in response to DNA-strand breaks induced by both the non-selective and selective treatments, S. cerevisiae induces genes involved in DNA repair and the general stress response and represses genes encoding for ribosomal biogenesis; (iii) a considerable fraction of the response upon exposure to γ-radiation was not mimicked by treatment with the enediynes; and (iv) the "DNA damage signature" as described by Gasch et al. [1] in S. cerevisiae is associated with single-strand breaks in DNA. These results have implications for our understanding of the cellular responses arising from damage to specific components of the cell.

Results and discussion

Genome-wide responses of S. cerevisiae upon exposure to enediynes and γ-radiation are partially overlapping

The expression patterns of yeast genes were monitored in S. cerevisiae exposed to the enediynes calicheamicin γ1I, esperamicin A1 and neocarzinostatin, and to γ-radiation. Care was taken in these studies to control the conditions of the treatments, including cytotoxicity (doses producing ~20–25% lethality; [see Additional files 3 and 4]), exposure time (15 min for all chemicals; for γ-radiation, cells were incubated for 13 min at ambient temperature following the ~2 min irradiation), and exposure temperature (ambient temperature for all exposures). The reproducibility of the dose-lethality relationships of the agents is indicated by the highly similar concentrations of calicheamicin and neocarzinostatin producing 20–25% lethality in the DBY747 strain employed for these studies and in another strain, BY4741: 1–1.4 nM vs. 0.8–1 nM for calicheamicin, and 3.3–4.2 nM vs. 2–3 nM for neocarzinostatin, respectively [see Additional files 3 and 4]. Statistically significant transcriptional changes were determined for treated cells versus untreated (water or methanol vehicle for irradiated and enediyne-treated cells, respectively) using local pooled error analysis with a correction for false discovery rate (Table 2). The stringent selection criteria for gene expression change (see Methods) identified changes in 225 ORFs (~4% of the genome) across all experiments. Hierarchical cluster analysis displays the similarities and differences between the expression profiles of the enediyne and γ-radiation treated S. cerevisiae (Figure 2A). The highest similarities in expression responses were seen for the three enediyne treatments, which showed more similarity to each other than to γ-radiation as visualized in the dendrogram in the hierarchical cluster (Figure 2A). The difference between the response of the enediynes and γ-radiation was also identified using principal component analysis (PCA) with the highest similarity found between the responses of the enediynes, which clustered separately from γ-radiation (Figure 2B). PCA showed that the largest percentage of variability within the dataset (58.2%) was attributable to treatment effect (Figure 2B). Although S. cerevisiae treated with the enediynes did show highly similar responses in genome-wide modulation relative to γ-radiation, there were minor differences within the responses induced by neocarzinostatin, calicheamicin and esperamicin. These differences influence the presence of the subdendrogram generated via hierarchical cluster analysis within the enediynes (Figure 2A), as well as the principal component analysis where neocarzinostatin response differs moderately from that of calicheamicin and esperamicin (Figure 2B). Although these minor differences exist within the response to enediyne treatment, the overall responses were highly similar to each other and strikingly different to that of γ-radiation.

Additional File 1. Differential enrichment of GO categories identified in the comparison of yeast responses to enediynes and γ-radiation. The table lists the gene ontology categories identified as significantly enriched within the four sectors identified in the comparison of the enediynes response to that of γ-radiation.

Format: PDF Size: 57KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional File 2. GO categories identified as associated with single-strand or double-strand damage. The table is a list of the gene ontology categories enriched in the gene sets that were identified as associated with single-strand or double-strand damage.

Format: PDF Size: 48KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional File 3. Cytotoxicity dose-response curves for S. cerevisiae strain DBY747 treated with enediynes. The figure plots the cytotoxicity dose-response curve for exposure of strain DBY747 to enediynes.

Format: PDF Size: 157KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional File 4. Cytotoxicity dose-response curves for S. cerevisiae strain BY4741 treated with enediynes. The figure plots the cytotoxicity dose-response curve for exposure of strain BY4741 to enediynes.

Format: PDF Size: 125KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Table 2. Gene expression changes in response to enediyne or γ-radiation treatmenta

thumbnailFigure 2. Hierarchical cluster and principal component analysis of damage-induced gene expression ratios. WT S. cerevisiae was exposed to damage induced by enediynes: calicheamicin γ1I (cal), esperamicin A1 (esp) and neocarzinostatin (ncs) or γ-radiation (gamma). Significantly modulated transcripts in any one experiment (see Methods) were identified resulting in 225 ORFs and clustered hierarchically (A) or analyzed for principal components (B).

Common and unique responses of S. cerevisiae upon exposure to enediynes or γ-radiation

As shown in Figure 3, comparisons of the gene expression responses of the collective enediynes to γ-radiation identify four primary sectors: (I) those that were similarly induced upon exposure to enediynes and γ-radiation; (II) those that were induced upon exposure to γ-radiation but repressed upon exposure to enediynes; (III) those that were similarly repressed upon exposure to enediynes and γ-radiation; and (IV) those that were repressed upon exposure to γ-radiation and induced by enediynes (Figure 3). Gene Ontology enrichment analysis identified biological processes significantly enriched within each of these four sectors (Figure 3 and [see 1]). In general, the absolute fold change of the enediyne treatment was lower than fold changes observed with γ-radiation. We hypothesize that this difference is based on the specificity of the damage generated by the enediynes versus the nonspecific nature of the cellular damage caused by exposure to γ-radiation.

thumbnailFigure 3. Scatter plot of relative transcriptional modulation in enediyne or γ-irradiated S. cerevisiae. The average expression ratio for enediynes was compared to that of γ-irradiated cells. Four primary sectors are identified (I-IV) demonstrating common or unique responses. Subsets of biological processes enriched (see Methods) are shown within each of the sectors [for complete list see 1].

Common responses were identified to treatment with both enediynes and γ-radiation representing 55% of the modulated transcripts (124 of 225 ORFs). Genes induced by both enediynes and γ-radiation (Sector I) were represented by 13 Gene Ontology categories, 7 of which are involved in DNA repair [see 1]. The most significant class enriched was that of Response to Stress (p < 3.7 × 10-5) and includes MSN2, MSN4, DUN1, and RAD51. Genes repressed by both enediynes and γ-radiation (Sector III) were represented by 31 Gene Ontology categories. Among these 31 categories, significant enrichment was seen in the Ribosomal Large Subunit Biogenesis (p < 0.006) [see 1]. Interestingly, both the induction of stress response genes and repression of ribosomal protein synthesis genes are two of the hallmarks of the Environmental Stress Response (ESR) in S. cerevisiae. The ESR was identified by Brown and colleagues where a commonly responding set of genes (~900) was modulated upon exposure to a wide variety of environmental stressors including temperature shock, hydrogen peroxide and hyper- and hypo-osmotic shock [15].

Despite the common response to treatment with enediynes and γ-radiation described above, a substantial fraction (45%) of the genomic response induced by γ-radiation was repressed by the enediynes and vise versa suggesting that much of the response to the enediynes is not "radiomimetic." Specifically 101 transcripts of the total 225 modulated transcripts demonstrated opposite directionality between γ-radiation and the enediynes. Genes that were induced by γ-radiation but repressed by enediynes (Sector II) were enriched for 19 Gene Ontology categories [see 1]. The most significantly enriched class was Mutagenesis (p < 2 × 10-6) and includes POL30 and CDC8. The final class included transcripts repressed upon exposure to γ-radiation and induced by enediynes (Sector IV). A significantly enriched Gene Ontology class includes Chromatin Assembly/Disassembly (p < 0.0007) including HTB2, HHT1, and HHT2 [see 1].

Regression analysis identified gene induction that associates with double-strand or single-strand damage

Using the various proportions of single- and double-strand deoxyribose oxidation data shown in Table 1, we performed regression analysis of gene expression responses to treatment with each of the enediynes and γ-radiation. The results identified 55 ORFs with transcriptional induction that associates either with double-strand or single-strand damage (Figure 4A and 4B, respectively). As shown in Table 3, positive association was identified with 29 ORFs demonstrating elevated expression in response to double-strand lesions and 27 ORFS with elevated expression in response to single-strand lesions (Table 3). While these correlations are significant (r2 > 0.6) for direct strand breaks, the correlation increases substantially when total deoxyribose oxidation is considered (r2 > 0.8; Table 3). The distinction between direct strand breaks and total deoxyribose oxidation involves the formation of oxidized abasic sites by γ-radiation and the enediynes. These lesions consist of the 2'-deoxyribonolactone and 2-deoxypentose-4-ulose derived from 1'- and 4'-oxidation of deoxyribose in DNA, respectively, and are present in more than 75% of calicheamicin- and neocarzinostatin-induced double-strand lesions DNA damage [12]. The abasic sites are not taken into account in most studies of "strand break" formation by oxidizing agents (e.g., see ref. [16]. Indeed, the 1:3 ratio of direct double- to single-strand breaks observed by Hammersten and coworkers for calicheamicin-induced damage in cells [16] is similar to the 1:1 to 1:2 ratio for direct breaks that they and others observed in plasmid DNA in vitro [12,16-18]. Hydrolysis of the abasic sites dramatically shifts the ratio of lesions toward double-strand breaks, as shown in Table 1. The abasic sites are DNA lesions that cause significant distortion of the DNA helix, so they may indeed be recognized by proteins involved in the detection of double-strand DNA lesions, as suggested by the stronger correlation observed for total deoxyribose oxidation than for direct strand breaks (Table 3), an hypothesis that warrants further study.

Table 3. Elevated expression of ORFs associated with single- and double-strand damage caused by exposure of S. cerevisiae to γ-radiation and enediynesa

thumbnailFigure 4. Trend line graph of ORFS with transcriptional modulation associated with single-strand or double-strand DNA breaks. ORFs with associated transcriptional modulation in response to treatment with enediynes or γ-radiation were identified using regression analysis (see Methods). A) ORFs with modulation associated with double-strand breaks and B) ORFs with modulation associated with single-strand breaks [for complete list of ORFs see 2].

Interestingly, ORFs with regulatory association with single-strand breaks show an enrichment of gene ontology categories including Response to Stress (HSP26, YCL033C, DUN1, RAD51, HYR1, ASF1, AHP1, CRS5), and DNA replication (POL1, RFA2) [see 2]. Within this single-strand break associated gene set there is an enrichment of genes belonging to the "DNA damage signature" in S. cerevisiae including DUN1, RAD51, RNR2, and YBR070C. The DNA damage signature, identified by Brown and colleagues, includes nine genes: the DNA damage repair genes RAD51 and RAD54, ribonucleotide reductase subunits RNR2 and RNR4, DNA damage activated kinase DUN1, and uncharacterized genes YER004W and YBR070C [1]. The "DNA damage" signature was identified by comparing expression programs in S. cerevisiae that were modulated by a range of environmental stresses including heat shock, oxidative stress, reductive stress, osmotic shock and amino acid starvation to those elicited specifically by DNA damage induced by the alkylating agent methylmethane sulfonate (MMS) and ionizing radiation (γ-radiation) [1,15]. Although it is not surprising that here we identify the DNA damage signature gene set as associated with γ-radiation (as this gene set was characterized in response to γ-radiation), it is noteworthy that there is a specificity of the activation of the DNA damage signature gene set to single-strand damage.

The induction of the DNA damage signature gene set was muted in cells deficient for the ATR homolog, mec1, as well as cells deficient for dun1, which indicates these genes are downstream and targets for the Mec1 pathway [1]. In contrast to the single-strand break associated gene set, genes associated with double-strand damage are enriched for gene ontology categories that include cell cycle regulation [see 2]. Interestingly, no DNA repair or stress responses categories were significantly enriched within this gene set. These results suggest that there is differential genomic modulation in response to single-strand or double-strand DNA damage in S. cerevisiae (Figure 5).

thumbnailFigure 5. Model of specific responses in S. cerevisiae to double- or single-strand breaks.

Comparisons with published studies of enediyne-induced changes in gene expression

Expression profiling experiments have been performed using calicheamicin [13] and neocarzinostatin [14]. Several major experimental differences make it impossible to directly compare the data generated here with those of the published studies. The calicheamicin studies performed by Schaus et al. employed the YPH500α strain of S. cerevisiae and significantly higher (8–80 nM) drug concentrations [13], while the neocarzinostatin studies of Watanabe et al. [14] utilized strain BY4741 and the holoantibiotic form of the drug (DNA-cleaving chromophore bound to 17,000 Da apo-protein; see ref. [11]). In addition to the absence of publicly available expression data sets for these studies, neither reported any measure of enediyne-induced cytotoxicity or other biological response that could be used to normalize the resulting transcriptional data.

Conclusion

Genome-wide expression profiling in S. cerevisiae has resulted in the understanding that a robust response is mounted upon exposure to agents that damage DNA. Here we discern the components of the transcriptional response of S. cerevisiae that are specifically due to particular types of DNA damage, namely that of single- or double-strand DNA damage arising from deoxyribose oxidation. These studies were conducted with a panel of three DNA-selective enediyne antibiotics (calicheamicin γ1I, esperamicin A1 and neocarzinostatin) that produce different proportions of double- and single-strand deoxyribose damage in DNA, with genome-wide responses compared to those induced by the non-selective γ-radiation. We find at doses producing similar toxicity, exposure to non-selective γ-radiation results in more extensive reprogramming of the S. cerevisiae transcriptome than exposure to each of the three enediynes. The extensive response to γ-radiation may reflect the non-specific nature of the oxidative attack on DNA (base and sugar), lipids, carbohydrates, proteins and small metabolites in cells. A striking finding is that only a modest fraction of the response upon exposure to γ-radiation was mimicked by treatment with the enediynes. We find that, in response to DNA strand breaks induced by both the non-selective and selective agents, yeast induce genes involved in DNA repair and stress response and repress genes encoding for ribosomal biogenesis. We also identify that the DNA damage signature in S. cerevisiae is more closely associated with single-strand breaks in DNA, than with double-strand breaks.

Methods

Yeast strains, culture, and reagents

Saccharomyces cerevisiae strain DBY747 (MATa, his3Δ1, leu2Δ112, ura3Δ52, trp1Δ289a gals, Can1, CUPr) was used in this study. Enediynes were used as methanolic stocks and obtained as follows: Calicheamicin was obtained from Wyeth Research (Cambridge, MA), neocarzinostatin was obtained from Kayaku Co. Ltd. (Tokyo, Japan; no longer available; currently produced by Sigma Chemical Company, St. Louis, MO), esperamicin A1 was obtained from Bristol Myers Squibb (Wallingford, CT; no longer available). The neocarzinostatin chromophore was isolated by methanol extraction and the concentration determined as described elsewhere [19]. Cells were grown and maintained in YPD (10 g yeast extract, 20 g peptone, 20 g dextrose, 20 g agar per liter). As shown in [3], samples of the DBY747 WT strain were grown to mid-log phase and treated with doses of agents determined to produce ~20–25% lethality following a 15 min exposure: calicheamicin (1 nM), neocarzinostatin (4 nM), esperamicin A1 (3 nM), and γ-radiation (210 Gy delivered at 105 Gy/min in a 60Co source with 13 min post-irradiation incubation at ambient temperature). Similar studies were performed with the BY4741 strain, with similar results [see 4]. All exposures were performed at ambient temperature under identical conditions of light and temperature, and cells were snap frozen in liquid nitrogen immediately after the 15 min exposure period.

RNA preparation/cRNA synthesis

Total RNA was extracted from pelleted cells using a hot phenol protocol [20] that has been used successfully in previous expression profiling studies [5,9,10], including one with calicheamicin [14]. AE buffered phenol (50 mM sodium acetate pH 5.3, 10 mM EDTA, pH 8.0) was utilized. RNA was precipitated from the aqueous phase of the phenol extraction. Labeled cRNA was generated as follows. Total RNA was converted into single-stranded cDNA using a modified oligo(dT) primer with a 5' T7 RNA polymerase promoter sequence and reverse transcriptase (SuperScript II RT, Gibco). Double-stranded cDNA was generated using DNA polymerase and DNA ligase (Invitrogen Life Technologies) and purified. Biotin-labeled cRNA was generated using in vitro transcription with T7 RNA polymerase (ENZO BioArray HighYield RNA transcript Labeling Kit, Affymetrix, CA). Reactions were carried out for 5 h at 37°C and cRNA was purified using RNeasy spin columns (Qiagen, CA). cRNA was quantified at UV260 and 15 μg of RNA was fragmented randomly using (200 mM Tris-Acetate, 500 mM potassium acetate, 150 mM magnesium acetate) at 94°C for 35 min.

GeneChip® hybridizations and image analysis

Hybridizations were performed as follows. Fragmented cRNA was hybridized to GeneChip® (YG-S98, Affymetrix, CA) at a concentration of 0.05 μg/μl in 200 μl of Affymetrix buffer (100 mM MES, 1 M NaCl, 20 mM EDTA, 0.01% Tween 20) with GeneChip® eukaryotic hybridization controls (GeneChip® Eukaryotic Hybridization Controls Kit, Affymetrix, CA) in the presence of 0.1 mg/ml herring sperm DNA and 0.5 mg/ml acetylated BSA at 40°C for 16 h with constant rotation. Arrays were rinsed after hybridization with 200 μl of stringent wash buffer (100 mM MES, 0.1 M NaCl, 0.01% Tween 20) followed by a non-stringent wash (6XSSPE, 0.01% Tween 20). 20XSSPE had the following composition (3 M NaCl, 0.2 M NaH2PO4, 0.02 M EDTA). Staining was done with 2 μg/ml streptavidin-phycoerytherin and 1 mg/ml acetylated BSA in 6xSSPE-T. Arrays were scanned using a HP G2500A GeneArray scanner.

Data pre-processing and differential gene testing

Untreated and treated samples were analyzed in biological triplicate on YG-S98 arrays. Quantile normalization was carried out using the Robust Multichip Average (RMA) algorithm [21]. Transcripts that were absent across all experiments were identified using Absent/Present calls (Microarray Suite 5.0) and filtered for transcripts that were not expressed in any experiment. Differential gene expression was calculated using the dual filtering criteria of ≥1.5 fold change and statistical significance determined through the Local Pooled Error test (LPE) [22] with an adjustment for false discovery rate calculation of p value of ≤0.10 [23]. All microarray data have been submitted to the Gene Expression Omnibus Database (Series, GSE5301).

Gene Ontology Enrichment Analysis

Statistical evaluation of co-regulated groups of genes was carried out through Gene Ontology Enrichment Analysis within the Functional Specification Database (funspec.med.utoronto.ca) [24]. Co-regulated yeast ORFs were classified according to the Gene Ontology Biological Process and the hypergeometric distribution (p < 0.01) used to assess enrichment of a particular gene category. The hypergeometric distribution asses for each gene ontology category, the probability (p-value) of observing such an overlap by chance is calculated as:

<a onClick="popup('http://www.biomedcentral.com/1471-2164/7/313/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2164/7/313/mathml/M1">View MathML</a>

where G is the size of the genome, C is the number of genes in the genome having that attribute, n is the size of the query cluster, of which k are known to possess the attribute.

Regression analysis

To identify genes with expression associated with double- or single-strand damage, all modulated ORFs across experimental conditions (225 total) were analyzed for association with induced expression ratio relative to percent single or double-strand lesion formation. Linear regression was performed on expression ratios relative to percent of direct single-strand or double-strand breaks with R2>0.6 and percent of total deoxyribose oxidation events with R2>0.8 selected for associated gene sets (Figure 4).

Authors' contributions

RF carried out the analysis of all genomics experiments and drafted the manuscript. MD drafted the manuscript and participated in the analysis of the genomics experiments. JC carried out the experimental treatments of S. cerevisiae and the hybridizations to GeneChips. TB assisted in the experimental treatments and gene expression profiling experiments. LS participated in the oversight of the study and helped to draft the manuscript. PD conceived of the study and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.

Acknowledgements

The authors wish to thank Ms Mary Montesalzo and Dr. Izle Skare at the University of Massachusetts at Lowell for expert assistance with the γ-irradiations. This work was supported by NIH grant CA103146 (PCD). Transcriptional profiling studies were performed in the Genomics and Bioinformatics Core of the MIT Center for Environmental Health Sciences supported by NIEHS grant ES002109. Support for RF and TB was provided by U19-ES11399 and CA55042.

References

  1. 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.

    Mol Biol Cell 2001, 12(10):2987-3003. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Friedberg EC, Walker GC, Siede W: DNA Repair and Mutagenesis. Washington, D.C. , ASM Press; 1995.

  3. Wood RD, Mitchell M, Lindahl T: Human DNA repair genes, 2005.

    Mutat Res 2005, 577(1-2):275-283. PubMed Abstract | Publisher Full Text OpenURL

  4. Wood RD, Mitchell M, Sgouros J, Lindahl T: Human DNA repair genes.

    Science 2001, 291(5507):1284-1289. PubMed Abstract | Publisher Full Text OpenURL

  5. Begley TJ, Rosenbach AS, Ideker T, Samson LD: Damage recovery pathways in Saccharomyces cerevisiae revealed by genomic phenotyping and interactome mapping.

    Mol Cancer Res 2002, 1(2):103-112. PubMed Abstract | Publisher Full Text OpenURL

  6. Begley TJ, Rosenbach AS, Ideker T, Samson LD: Hot spots for modulating toxicity identified by genomic phenotyping and localization mapping.

    Mol Cell 2004, 16(1):117-125. PubMed Abstract | Publisher Full Text OpenURL

  7. 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.

    Proc Natl Acad Sci U S A 2002, 99(16):10605-10610. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Birrell GW, Giaever G, Chu AM, Davis RW, Brown JM: A genome-wide screen in Saccharomyces cerevisiae for genes affecting UV radiation sensitivity.

    Proc Natl Acad Sci U S A 2001, 98(22):12608-12613. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. 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.

    Mol Cell Biol 2000, 20(21):8157-8167. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Jelinsky SA, Samson LD: Global response of Saccharomyces cerevisiae to an alkylating agent.

    Proc Natl Acad Sci U S A 1999, 96(4):1486-1491. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. Dedon PC, Goldberg IH: Free-radical mechanisms involved in the formation of sequence-dependent bistranded DNA lesions by the antitumor antibiotics bleomycin, neocarzinostatin, and calicheamicin.

    Chemical Research in Toxicology 1992, 5(3):311-332. PubMed Abstract | Publisher Full Text OpenURL

  12. Dedon PC, Salzberg AA, Xu J: Exclusive production of bistranded DNA damage by calicheamicin.

    Biochemistry 1993, 32(April 13):3617-3622. PubMed Abstract | Publisher Full Text OpenURL

  13. Schaus SE, Cavalieri D, Myers AG: Gene transcription analysis of Saccharomyces cerevisiae exposed to neocarzinostatin protein-chromophore complex reveals evidence of DNA damage, a potential mechanism of resistance, and consequences of prolonged exposure.

    Proc Natl Acad Sci U S A 2001, 98(20):11075-11080. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Watanabe CM, Supekova L, Schultz PG: Transcriptional effects of the potent enediyne anti-cancer agent Calicheamicin gamma(I)(1).

    Chem Biol 2002, 9(2):245-251. PubMed Abstract | Publisher Full Text OpenURL

  15. 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.

    Mol Biol Cell 2000, 11(12):4241-4257. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Elmroth K, Nygren J, Martensson S, Ismail IH, Hammarsten O: Cleavage of cellular DNA by calicheamicin gamma1.

    DNA Repair (Amst) 2003, 2(4):363-374. PubMed Abstract | Publisher Full Text OpenURL

  17. Drak J, Iwasawa N, Danishefsky S, Crothers DM: The carbohydrate domain of calicheamicin g1I determines its sequence specificity for DNA cleavage.

    Proc Natl Acad Sci USA 1991, 88:7464-7468. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Walker S, Landovitz R, Ding WD, Ellestad GE, Kahne D: Cleavage behavior of calicheamicin g1 and calicheamicin T.

    Proc Natl Acad Sci USA 1992, 89(May):4608-4612. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Xi Z, Mao QK, Goldberg IH: Mechanistic studies on the base-catalyzed transformation of neocarzinostatin chromophore: roles of bulged DNA.

    Biochemistry 1999, 38(14):4342-4354. PubMed Abstract | Publisher Full Text OpenURL

  20. Boffa LC, Bolognesi C: Methylating agents: their target amino acids in nuclear proteins.

    Carcinogenesis 1985, 6(9):1399-1401. PubMed Abstract OpenURL

  21. Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP: Summaries of Affymetrix GeneChip probe level data.

    Nucleic Acids Res 2003, 31(4):e15. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Jain N, Thatte J, Braciale T, Ley K, O'Connell M, Lee JK: Local-pooled-error test for identifying differentially expressed genes with a small number of replicated microarrays.

    Bioinformatics 2003, 19(15):1945-1951. PubMed Abstract | Publisher Full Text OpenURL

  23. Benjamini Y Hochberg, Y.: Controlling the false discovery rate: A practical and powerful approach to multiple testing.

    J R Stat Soc Ser B Methodol 1995, 57:289-300. OpenURL

  24. Robinson MD, Grigull J, Mohammad N, Hughes TR: FunSpec: a web-based cluster interpreter for yeast.

    BMC Bioinformatics 2002, 3:35. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  25. Fulford J, Nikjoo H, Goodhead DT, O'Neill P: Yields of SSB and DSB induced in DNA by Al(K) ultrasoft X-rays and alpha-particles: comparison of experimental and simulated yields.

    Int J Radiat Biol 2001, 77(10):1053-1066. PubMed Abstract | Publisher Full Text OpenURL

  26. Milligan JR, Aguilera JA, Paglinawan RA, Ward JF, Limoli CL: DNA strand break yields after post-high LET irradiation incubation with endonuclease-III and evidence for hydroxyl radical clustering.

    Int J Radiat Biol 2001, 77(2):155-164. PubMed Abstract | Publisher Full Text OpenURL

  27. Yu L Golik, J., Harrison, R., and Dedon, P.: The deoxyfucose-anthranilate of esperamicin A1 confers intercalative DNA binding and causes a switch in the chemistry of bistranded DNA lesions.

    Journal of the American Chemical Society 1994, 116:9733-9738. Publisher Full Text OpenURL

  28. Dedon PC, Goldberg IH: Influence of thiol structure on neocarzinostatin activation and expression of DNA damage.

    Biochemistry 1992, 31(7):1909-1917. PubMed Abstract | Publisher Full Text OpenURL