|
|
|
|
Vol. 15, Issue 9, 4321-4336, September 2004
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||



* Institute of Toxicology, Swiss Federal Institute of Technology ETH, CH-8603 Schwerzenbach, Switzerland;
Pharma Division, F. Hoffmann-La Roche Ltd., CH-4070 Basel, Switzerland; and
Ordway Research Institute, Albany, NY 12208-03479
Submitted May 7, 2004;
Revised June 9, 2004;
Accepted June 14, 2004
Monitoring Editor: Keith Yamamoto
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
AFB1 is a mutagen in Saccharomyces cerevisiae (Sengstag et al., 1996
), Escherichia coli, rainbow trout, mice, rat and human cells (reviewed in Smela et al., 2001
), and a recombinagen in yeast and in human cells (Stettler and Sengstag, 2001
). In yeast, AFB1 can induce mitotic, homologous recombination resulting in heteroallelic gene conversion and translocations (Sengstag et al., 1996
). After yeast cells are exposed to low doses of AFB1 in the expected range of human exposure, there is a strong stimulation of recombination but not mutation (unpublished data). In human lymphoblastoid cell line TK6, AFB1 exposure increases heteroallelic recombination at the thymidine kinase locus resulting in loss of heterozygosity (Stettler and Sengstag, 2001
). Thus, understanding the molecular basis for the recombinogenicity of AFB1 in yeast may help understand the potent carcinogenicity of AFB1 compared with toxins with similar mutagenicity.
The remarkable recombinogenicity of AFB1 may result from a combination of factors. First, specific AFB1-DNA adducts may enzymatically or spontaneously convert to DNA double-strand breaks, thus directly initiating recombination. The N7 adduct 8,9-dihydro-8-(N7-guanyl)-9-hydroxyaflatoxin B1 is the major product in vitro (Essigmann et al., 1977
) and in vivo (Lin et al., 1977
; Croy et al., 1978
). The positively charged imidazole ring of the principal DNA adduct promotes depurination, giving rise to an apurinic (AP) site, which can further yield single-strand breaks by
-elimination (Friedberg et al., 1995
). Clusters of these single-strand breaks could yield double-strand breaks. Alternatively, mildly alkaline conditions can subsequently result in the formation of a chemically and biologically stable foramidopyrimidine derivative (AFB1-FAPY), which represents a significant product in vivo (Croy and Wogan, 1981
). AP sites can be removed by the base excision repair (BER) pathway, and the AFB1-N7-guanine adducts can be removed by the nucleotide excision repair (NER; Leadon et al., 1981
). The AFB1-FAPY adduct, however, is a nonrepairable, persistent lesion (Martin and Garner, 1977
) that interferes with DNA replication. Such interference could indirectly stimulate recombination (Friedberg et al., 1995
) and generate DNA double-strand breaks. However, chromosomal fragments have not been detected by pulse-field electrophoresis after yeast cells were exposed to AFB1 (unpublished data).
Alternatively, exposure to AFB1 could also elicit a stress response in yeast that stimulates more recombination than mutation. We thus investigated the global cellular response to a 4-h. exposure to AFB1. DNA microarrays have been used successfully in yeast to investigate the global transcriptional response after exposure to saline (Posas et al., 2000
), methyl methanesulfonate (MMS; Jelinsky and Samson, 1999
, Gasch et al., 2000
), and ionizing radiation (Gasch et al., 2001
). The current mRNA expression analysis shows that a large fraction of the AFB1-induced genes is involved in maintenance of DNA integrity. Because the majority of the transcriptionally upregulated DNA repair genes belong to the NER or recombinational repair (RR) pathway, we exposed the respective rad1 and rad51 repair mutants to AFB1 and measured translocation frequencies. To strengthen the correlation between AFB1-associated recombination and RAD51 induction, we measured AFB1-associated recombination in mec1 checkpoint mutants, defective in the DNA damage inducibility of RAD51, and in mec1 mutants expressing higher basal levels of RAD51. Our data suggest that AFB1 upregulates a recombinational repair pathway that involves RAD51 and RAD1.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Yeast strains contain two overlapping his3 fragments on chromosomes II and IV and were derived from YB109 (Fasullo and Dave, 1994
). Translocation frequencies were determined by selecting for His+ recombinants that are generated by mitotic recombination between the his3 fragments. (Fasullo and Davis, 1987
). YMK2181 (MATa/MAT
, ura3-52/ura3-52, his3-
200/his3-
200, ade2-101/ade2-101, trp1-
1/TRP1, gal3-/gal3-, leu2-3112/leu2-3112, GAL1::his3-
5'/GAL1::his3-
5', trp1::his3-
3'/trp1::his3-
3', leu2-
3', leu2-
5', kanMX4, HOcs), and YB110 (MATa/MAT
, ade2-101/ade2-101, ura3-52/ura3-52, his3-
200/his3-
200, trp1-
1/trp1-
1, leu2/LEU2, GAL1::his3-
5'/GAL1, trp1::his3
3'/trp1-
1, LYS2/lys2-801; Fasullo and Dave, 1994
) have been previously used to measure DNA damage-associated translocation. YB150 (rad1) and YB195 (rad51) are identical to YB110 (Rad+) except for the rad51 and rad1 disruptions, respectively. We replaced the ade2-101 allele in YB109 and with ade2-n (YB318) and the ade2-101 allele in YA102 with ade2-a (YB336) by two-step gene replacement using the plasmid pKH9 (Huang and Symington, 1994
).
mec1 checkpoint mutants that measure AFB1-associated translocations contain either mec1-21 or the mec1 null mutation. The original MAT
mec1-21 (YA16) strain is derived from W303 (Sanchez et al., 1996
). We backcrossed YA16 10 times with strains in the S288c background (YB163 and FY251 [Dong and Fasullo, 2003
] and YB336) to generate meiotic segregants YB316 (MAT
ura3-52 his3-
200, trp1-
1, ade2-a, mec1-21) and YB314 (MAT
ura3-52 his3-
200, trp1
-1, ADE2, mec1-21) by tetrad dissections. YB318 was crossed with YB314 to generate the meiotic segregant YB319 (MATa-inc, ura3-52, his3-
200, ade2-n, trp1-
1, leu2, lys2, GAL1::his3-
5', trp1::his3-
3', mec1-21). YB325 (MATa/MAT
, ade2-a/ade2-n, ura3-52/ura3-52, his3-
200/his3-
200, trp1-
1/trp1-
1, leu2/LEU2, GAL1::his3-
5'/GAL1, trp1::his3-
3'/trp1-
1, lys3801/lys2-801, mec1-21/mec1-21) was then used to measure translocations and heteroallelic recombination in the mec1-21 background. To measure translocations in the mec1 null mutant, we first introduced the sml1::kanMX allele in YB318 and YB315 by PCR-mediated gene replacement (Goldstein and McCusker, 1999
) to make YB320 and YB317, respectively, because lethality conferred by mec1 deletions is suppressed by sml1 mutations (Zhao et al., 1998
). YB323 (MATa/MAT
, ade2-a/ade2-n, ura3-52/ura3-52, his3-
200/his3-
200, trp1-
1/trp1-
1, leu2/LEU2, GAL1::his3-
5'/GAL1, sml1::kanMX/sml1::kanMX, trp1::his3
3'/trp1-
1, lys2801/lys2801) was then derived by a diploid cross of YB320 and YB317. The mec1
::TRP1 allele (Zhao et al., 2000
) was introduced into YB320 and YB317 to make YB321 and YB322, respectively. YB324 (MATa/MAT
, ade2-a/ade2-n, ura3-52/ura3-52, his3-
200/his3-
200, trp1-
1/trp1-
1, leu2/LEU2, GAL1::his3-
5'/GAL1, sml1::kanMX/sml1::kanMX, trp1::his3
3'/trp1-
1, lys2801/lys2801, mec1
::TRP1/mec1
::TRP1) was then derived by a diploid cross. To overexpress RAD51 in the mec1 mutants, pR51.3 (Leu+), containing RAD51 on a 2 µ plasmid, was introduced into YB325 (Sung and Stratton, 1996
).
The 2 µ URA3 plasmids pMK637 (this work) or pSB229 (Eugster et al., 1992
), containing hCYP1A2+hOR and hCYP1A1+hOR cDNAs, respectively, or the LEU2 plasmid pCS512 (Sengstag et al., 1996
), containing hCYP1A1+hOR cDNAs, were first introduced into yeast strains by DNA transformation to metabolically activate the AFB1 and benzo-(a)-pyrene-7,8-dihydrodiol (BaP-DHD; Klebe et al., 1983
). The 2 µ URA3 plasmid pCS316, containing the hCYP1A1+hOR cDNA in the opposite orientation as in pSB229 (Eugster et al., 1992
), was introduced into YB110, YB324, and YB335 to measure AFB1-associated translocations in mec1 checkpoint mutants. pMK637 was introduced into the strain YMK2181 to measure AFB1 related-changes in gene expression using the oligonucleotide arrays. pCS512 was introduced into YB150 and the plasmid pSB229 was introduced into YB195 and YB110 to measure chromosomal translocation frequency and drug killing after exposure to ethyl methanesulfonate (EMS), AFB1, and BaP-DHD.
Exposure of Yeast Strains to DNA-damaging Agents
In brief, exponentially growing yeast cells were collected by centrifugation and resuspended in 0.1 M sodium phosphate buffer (pH 7.5); the final cell density was 4 x 108 cells/ml. To measure the stimulation of recombination, 1 ml of the cells in 0.1 M sodium phosphate buffer (pH 7.5) was exposed to chemicals for 4 h. at 30°C in a rotary shaker. The cells were then pelleted in a clinical centrifuge, washed, and diluted in supplemented minimal medium. To measure the net frequencies of recombination, the spontaneous frequencies were subtracted from the DNA damageassociated frequency. To measure AFB1-associated changes in gene expression, 2 ml of cells in 0.1 M sodium phosphate buffer was exposed to 25 µM AFB1 for 4 h at 30°C in a rotary shaker. Cells were then centrifuged and resuspended in the appropriate buffers to extract nucleic acids.
Preparation of Nucleic Acids for Oligonucleotide Arrays and Hybridization
After AFB1 exposure, cells were washed once, resuspended in 0.5 ml RLT buffer (Qiagen GmbH, Hilden, Germany) supplemented with 1% mercaptoethanol (Riedel-deHaën, Hannover, Germany) and transferred to a glass tube. Acid-washed glass beads (Ø 0.450.55 mm, Merck, Darmstadt, Germany) were added up to the meniscus and the cells were disrupted by heavy vortexing three times for 3 min. After addition of 3.3 ml RLT buffer, the lysate was recovered with a glass capillary. Total RNA was isolated using the RNeasy Midi Kit (Qiagen AG, Basel, Switzerland) according to the manufacturer's protocol. RNA quality was assessed on an agarose gel. Poly(A)+ RNA was amplified and biotin-labeled as follows. Starting with 20 µg total RNA, double-stranded cDNA was constructed using the GibcoBRL Superscript choice system (Life Technologies AG, Basel, Switzerland) and a T7-(T)24 primer to introduce a T7 promoter. Double-stranded cDNA was purified by three successive phenol:chloroform:isoamyl alcohol extractions and a subsequent alcohol precipitation. Phase-Lock Gel (5 Prime to 3 Prime, Boulder, CO) was used for all organic extractions to increase recovery. Using
0.20.5 µg cDNA as a template, a biotin-labeled riboprobe was synthesized with the help of the T7 Megascrip system (Ambion, Austin, TX) and two biotin-labeled nucleotides (Bio-11-CTP and Bio-16-UTP, Enzo Diagnostics, Farmingdale, NY), which replaced one third of the provided CTP and UTP. The 6-h in vitro transcription reaction yielded
50 µg cRNA, which was purified by RNA affinity resin (RNeasy spin columns, Qiagen). An aliquot was separated on a 0.8% agarose gel to check sample integrity. Subsequently, 40 µg of the transcript were used to hybridize a set of four commercially available oligonucleotide expression arrays (GeneChip Ye6100 arrays, Affymetrix, Santa Clara, CA) comprising a total of more than 260,000 oligonucleotides complementary to 6218 yeast open reading frames (ORFs). The biotinylated cRNA samples were fragmented to increase hybridization efficiency and specificity and to reduce potential problems caused by nucleic acid secondary structure (Wodicka et al., 1997
). Chip hybridization, washing, and staining with a streptavidin-phycoerythrin conjugate were performed using Affymetrix instrumentation according to the company's recommended protocols. The arrays were read at 7.5 µm with a confocal scanner (Molecular Dynamics, Sunnyvale, CA) and analyzed with GENECHIP software, version 3.0. A threshold of 20 arbitrary fluorescence units was assigned to any gene with a calculated expression level below 20, because discrimination of mRNA levels in this low range could not be performed. Chip hybridization and mRNA quality were verified with controls on the arrays consisting of 3', middle, and 5' regions of housekeeping genes (actin, SPT15, SRB4) and marker oligonucleotides at the corners, edges, and in the middle of the array (Wodicka et al., 1997
; unpublished data).
Statistical Analysis of the AFB1/DMSO Data Set
mRNA levels were expressed as the average difference of hybridization signals, measured as fluorescence intensity, between perfect match and central-mismatch oligonucleotide probe sets (Wodicka et al., 1997
), and supplemented with an absent/present call generated by the Affymetrix software. Data from different chips were normalized using the parameter of total chip signal. We calculated the mean of the average differences of two chips each of AFB1 (AFB1+) and solvent (AFB1) exposed cells. Only ORFs deviating <40% of this mean value (purity
0.6) were used for further analysis; 5630 ORFs fulfilled this criterion. The data sets were then imported into a MS Excel spreadsheet for further calculations and logical operations.
Preparation of RNA for Quantitative PCR Analysis
RNA was extracted from control cells, and cells were exposed to AFB1 (Shirra et al., 2001
). RNA quality was assessed on a 0.8% agarose gel. DNaseI (0.05 U/ml, BD Biosciences, San Diego, CA) was added to ensure that no DNA was present in the extraction and after digestion at 37°C for 30 min, was inactivated in 1 mM EDTA (pH. 8.0). After extraction in phenol:chloroform:isoamyl alcohol (25:24:1, pH 4.5) and chloroform extraction, the aqueous layer was precipitated in 0.2 M NaOAc, 70% EtOH. The RNA pellet was then resuspended in TrisEDTA. One milligram of RNA was used for the reverse transcription reaction (first-strand cDNA synthesis), using a protocol described in the reverse transcription system kit (Promega, Madison, WI). cDNA was measured in a iCycler (Bio-Rad, Richmond, CA) by quantitative PCR (QPCR) using the IQ Green SYBR supermix kit (Bio-Rad). Cycle conditions included denaturation at 95°C, followed by 35 cycles of 95°C denaturation, 57°C reannealing, and 72°C reaction; a 95°C denaturation step; and a 55°C reannealing step. Rad51 cDNA was measured using oligos 5'-CAACTTGGGCGACCACTT G-3' and 5'-AAAGGCTGGCCGACCAAT-3'. Act1 cDNA was measured using oligos 5'-CCACCAATCCAGACGGAGACT-3' and 5'-GCCGAAAGAATG CAAAAG GA-3'. Rad1 cDNA was measured using 5'-CTAATTGTGCCTCATCGACCAA-3' and 5'-GGATGCCAATAAACCGTCAGTATC-3'.
Measurements of DNA Damageassociated Recombination Frequencies in Checkpoint and rad Mutants and in Wild Type
We measured the frequency AFB1, EMS, and BaP-DHDassociated translocations and drug toxicity in the rad mutants, YB195pSB229 (rad51) and YB150pCS512 (rad1); checkpoint mutants, YB324pCS316 (sml1, mec1) and YB325pCS316 (mec1); and the Rad+ proficient strain YB110pCS316, as previously described (Sengstag et al., 1996
). YB195, YB150, YB324, and YB325 transformants were grown in YM His-Ade-Trp-Lys and YB110 transformants in YM His-Ade-Trp. After exposure to chemical agents, cells were resuspended to a density of 8 x 108 cells/ml, 100250 µl was plated directly on YM Ade-Ura-Trp-Leu-Lys to select for His+ recombinants, and the appropriate dilution was plated on YPD to measure viability. Selection plates were incubated at 30°C, and the colonies were counted after 7 days.
Chemicals
Benzo-(a)-pyrene-7,8-dihydrodiol (BaP-DHD; Midwest Research Institute, Kansas City, MO) and aflatoxin B1 (AFB1, Fluka, Buchs, Switzerland) were dissolved in DMSO. Ethyl-methane-sulfonate (EMS) was obtained from Eastman Kodak (Rochester, NY). DNA modifying enzymes were obtained from New England Biolabs, Inc. (Beverly, MA), 5-fluoroorotic acid (FOA) from Toronto Research Chemicals Inc. (Toronto, Ontario, Canada) and zymolyase was purchased from Seikagaku Corp. (Tokyo, Japan).
| RESULTS |
|---|
|
|
|---|
|
Of 6218 ORFs, 183 (2.9%) showed at least a threefold change in transcript levels after AFB1 exposure. One hundred seventeen were upregulated; the strongest induction was 35.1-fold (YEL068C; Table 1). Of the 66 genes that were downregulated, the strongest repression was 14.1-fold (YDR306C; Table 2). Most of the products of the responsive genes are located in the nucleus (79%; Figure 2). To gain an overview of the transcriptional response to AFB1, the ORFs were assigned to functional categories according to the MIPS database (Munich Information Center for Protein Sequences; Mewes et al., 1997
). The category of cell growth, cell division, and DNA synthesis contains 30 AFB1-responsive genes, the most number of any category. The second most numerous is the metabolism category, which contains 28 AFB1-responsive genes. However, comparing the percentage of responsive genes in each respective category, the category of cell rescue, defense, cell death, and aging contains the largest percentage of AFB1-responsive genes (4.8%; Figure 2). A more detailed view is provided by the analysis of the subcategories (Figure 3). Thus, there are genes in several functional categories whose expression is either induced or repressed after AFB1 exposure.
|
|
|
|
Because the genotoxicity of AFB1 may result from AFB1-induced DNA damage, we identified AFB1-inducible DNA repair genes. Of 109 genes involved in DNA damage repair according to the Yeast Protein Database (Payne and Garrels, 1997
), 15 (14%) were upregulated (RAD51, CDC2(POL3), DST1, RAD3, RSP5, RFA1, RAD16, MLH1, MMS21, DIN7, MET18, HPR5, RFA2, MSH6, RAD1) and 3 were repressed (DDR48, SIR4, DNL4) at least twofold. Furthermore, of the 16 genes assigned to specific repair pathways, 11 genes (69%) function in recombinational repair, 7 genes (44%) in nucleotide excision repair, and 4 genes function in both pathways. Only 2 of the repressed genes are not in either pathway but function in nonhomologous end joining (NHEJ; Table 3). Analysis of cell cycle periodicity of the repair genes showed that changes in expression levels after AFB1 exposure are not simply caused by changes in cell cycle progression (Keller-Seitz, 2001
). Furthermore, 44 genes exhibiting at least a twofold or greater change in expression are involved in damage signaling, stress response, or cell cycle control (Table 4). Thus, AFB1 exposure induces DNA repair genes in NER, MMR, and recombinational DNA repair.
|
|
AFB1-inducible repair genes that function in multiple recombination pathways include RAD51 (7.7-fold increased) and RAD1 (2.1-fold increased). We confirmed that RAD51 RNA increases after AFB1 exposure by QPCR, using actin RNA as a control (Figure 4). YB110 (pCS316) was exposed to 25 µm AFB1 for 4 h, and RNA was extracted for QPCR. The amount of RAD51 RNA increased fivefold in YB110 cells treated with AFB1, whereas actin RNA did not significantly increase (Figure 4). However, we found that the amount of RAD1 mRNA increased less than twofold (unpublished data). The differences between the QPCR results and the microarray results are likely due to the greater sensitivity of the microarrays (Etienne et al., 2004
).
|
Although we did not prove that the induced expression of the NER genes and RAD51 is necessary for AFB1-associated recombination, we did use rad1 (YB150 pCS512) and rad51 mutant (YB195 pSB229) yeast strains to determine whether RAD1 or RAD51 function in AFB1-associated recombination and lethality (Figure 5). We measured the frequencies of AFB1-associated translocations by selecting for His+ recombinants as previously described (Fasullo and Davis, 1987
). Besides AFB1, the carcinogen benzo-(a)-pyrene-7,8-dihydrodiol (BaP-DHD) and the mutagen ethyl methanesulfonate (EMS) were also tested. Although the viability of both mutant strains is slightly decreased after exposure to AFB1, the rad51 strain is hypersensitive to EMS and the rad1 strain is hypersensitive to BaP-DHD (Figure 5). Compared with the wild-type (Rad+) strain, the frequency of AFB1-associated translocations was slightly decreased in both the rad51 and rad1 mutants, whereas the frequency of EMS-associated recombination increased in the rad51 strain and the frequency of BaP-DHDassociated translocations was significantly higher in the rad1 mutant (Figure 5). At EMS concentrations greater than 40 mM, EMS-associated recombination could not be measured because of the extreme EMS toxicity. Because the frequencies of spontaneous recombination are (1.4 ± 0.4) x 107 and (3.8 ± 1.8) x 108 in the wild-type and rad1 strains, respectively, and low compared with the DNA damageassociated frequencies, the DNA damageassociated frequencies are similar to the net recombination frequencies. In the rad51 mutant the spontaneous frequency (avg.) was (1.3 ± 0.2) x 106, and thus the net recombination frequencies (avg.) for the highest level of AFB1-associated, EMS-associated, and BaP-DHDassociated translocations were 3.4 x 106, 96 x 106, 7.7 x 106, respectively. The highest net AFB1-associated frequency is still
25% lower in the rad51 diploid than in wild-type, whereas the highest net BaP-DHDassociated frequency was about threefold higher in the rad51 diploid than in wild type. These results indicate that RAD51 and RAD1 function in AFB1-associated recombination, whereas RAD51 and RAD1 suppress EMS-associated and BaP-DHDassociated recombination, respectively.
|
To further understand the correlation between RAD51 expression and AFB1-associated recombination, we measured translocation frequencies in mec1 checkpoint strains. mec1 checkpoint mutants are deficient in RAD51 induction after MMS and x-ray exposure (Gasch et al., 2001
). We found that the mec1-21 (YB325) and mec1
::TRP1 (YB324) mutants cannot induce RAD51 levels after exposure to 25 µM AFB1 (Figure 4). The translocation frequency increased 26-fold after wild type (YB110 pCS316) was exposed to 25 µM AFB1, consistent with previous studies (Sengstag et al., 1996
). We found no significant increase in translocation frequencies after both the mec1 deletion mutant (YB324 pCS316) and the mec1-21 mutant (YB325 pCS316) were exposed to 25 µM AFB1 (Table 5). However, because the sml1 null mutant exhibited a decrease in AFB1-associated recombination, the decrease in AFB1-associated recombination in the mec1 null mutant could be partially conferred by the sml1 mutation. To increase Rad51 in the mec1-21 strain (YB325), we introduced the 2 µ LEU2 plasmid (pR51.3; Sung and Stratton, 1996
) containing RAD51 expressed from a strong constitutive PGK (phosphoglycerol kinase) promoter by selecting for Leu+ transformants. By QPCR, we found that the basal level of RAD51 RNA before and after AFB1 exposure in the Leu+ transformants was the same and more than a thousand fold higher than the basal level of RAD51 RNA in mec1-21 (YB325), whereas there was no change in ACT1 RNA levels (unpublished data). In YB325 (pR51.3) cells, translocation frequencies increased ninefold after AFB1 exposure. RAD51 overexpression in mec1-21 also increased lethality after AFB1 exposure, suggesting that other detrimental recombination events may also be generated. These data indicate that an increase in RAD51 expression can enhance AFB1-associated recombination in a strain deficient in the DNA damage inducibility of RAD51.
|
| DISCUSSION |
|---|
|
|
|---|
The gene expression patterns reveal that AFB1-induced transcripts are not evenly distributed among yeast chromosomes (Figure 1). The differences ranged from
3-fold induction of 10% of ORFs on chromosome V to only 0.7% on chromosome XV. On the level of individual ORFs, a subset of genes that exhibit significant change in expression are closely linked; these include RAD3 (YER171W, 4.6-fold increased), which is located between ADK2 (YER170W, 5.1-fold increased) and BRR2 (YER172C, 2.2-fold increased), and RAD51 (YER095W, 7.7-fold increased), that is located close to UBP9 (YER098W, 7.4-fold increased) and SWI4 (YER111C, 5.1-fold increased). The linked ORFs are on different DNA strands (Goffeau et al., 1996
), suggesting that changes in chromosome structure may alter expression of multiple genes. Hence, our data suggest that some DNA damage responsive genes might be organized as clusters in coregulated chromosomal regions.
The complex response pattern caused by AFB1 reflects the broad range of toxic effects in the cell; however, the pattern of gene expression after exposure to AFB1 does not reflect a general stress response or a general response to DNA damaging agents, such as alkylating agents (Jelinsky and Samson, 1999
). Five stress response genes, DDR48, PAI3, YML131W, YKL100C, and YNL116W, that are upregulated after MMS exposure (Jelinsky and Samson, 1999
) and saline stress conditions (Posas et al., 2000
) are not induced after exposure to AFB1. The only gene assigned to the general stress response that was induced with MMS (17.8-fold) and AFB1 (2.6-fold) was the DNA damage-inducible gene DDI1. It is unlikely that the AFB1 solvent DMSO triggered the general stress response, because the only difference between the gene expression patterns after exposure to DMSO and H2O was the DMSO-dependent induction of YML131W. In addition, RNR1, which is generally induced after exposure alkylating agents and UV (Elledge and Davis, 1990
), was not induced after exposure to AFB1. We therefore suggest that the gene expression pattern of the cells exposed to AFB1 did not result from a general stress response but from specific AFB1-induced toxicity, such as AFB1-induced DNA damage.
The gene expression patterns may provide further insights into the recombinogencity of AFB1. Many of the responsive genes are directly or indirectly involved in recombination (reviewed in Aguilera et al., 2000
; Nicholson et al., 2000
). Several are also induced after diploid cells are exposed to
rays; these include RAD51, SRS2, RFA1, RFA2, and MSH6 (Mercier et al., 2001
). Among the AFB1-inducible DNA repair genes, RAD51 (7.7-fold increased) exhibited the strongest induction. The proteins encoded by the AFB1-inducible genes RFA1 and RFA2 are subunits of the replication factor A (RPA) and promote Rad51-stimulated DNA pairing and strand exchange in vitro (Sung, 1994
) by removing secondary DNA structures (Sung and Robberson, 1995
). Interestingly, numerous genes of the NER pathway were induced after AFB1 exposure; these included RAD1, RAD3, RAD16, and MET18. RAD1 functions in several RAD51-independent mitotic recombination events (Davies et al., 1995
; Saparbaev et al., 1996
; Paques and Haber, 1999
; Aguilera et al., 2000
; Haber, 2000
; Nicholson et al., 2000
) as well as in the spontaneous generation of homology-directed translocations (Fasullo et al., 1998
).
Several genes that are upregulated after AFB1 exposure also function to decrease particular recombination events. For example, the Hpr5/Srs2 helicase is suggested to function as an antirecombinase preventing excessive and aberrant RAD51-mediated recombination events (Klein, 2000
). In addition, some genes of the MMR pathway, including MLH1, MLH3, and MSH6, were induced, and Mlh1 and Mlh3 can both reduce recombination between repeated sequences containing mismatches (Nicholson et al., 2000
). Although the induction of these genes may seem contradictory to the notion that AFB1 stimulates recombination, HPR5/SRS2 is also upregulated in meiosis during which higher levels of heteroallelic and ectopic recombination occur. Mitotic, heteroallelic recombination is not decreased in mismatch repair mutants (Saparbaev et al., 1996
), and msh2 mutants do not exhibit decreased mitotic recombination between his3 fragments positioned on different chromosomes (unpublished data). Yeast Mlh1 interacts with Sgs1, a protein encoded by a human BLM homologue (Foury, 1997
), and may be involved in some aspect of general recombination (Pedrazzi et al., 2001
). We also speculate that the induction of the mismatch repair proteins may contribute to the weak mutagenicity of AFB1. Thus, the upregulation of MSH6 and HPR5 is consistent with the genotoxic properties of AFB1.
Besides genes encoding DNA repair functions, genes involved in damage signaling, stress response, and cell cycle progression were also upregulated after AFB1 exposure. Overexpression of SPO12, which was strongly induced (8.5-fold) after AFB1 exposure, is thought to reduce cyclin-dependent kinase activity and trigger exit from mitosis (Shah et al., 2001
). The observation that 9 (TOR1, TOR2, CTK1, MAD1, PTK1, PCL10, VPS34, PHO85, GAT1) of 43 genes displaying at least twofold change in expression have functions in TOR (target of rapamycin) signaling suggests a role of this pathway in response to AFB1 toxicity. Genes involved in TOR signaling, including VPS34 (phosphatidylinositol 3-kinase), are involved in regulatory mechanisms modulating protein synthesis and degradation and are important for promoting growth (Keith and Schreiber, 1995
; Thomas and Hall, 1997
; Dennis et al., 1999
). The significance of these genes may be further elucidated when specific cell cycle checkpoints are identified that are triggered by AFB1 exposure.
The downregulation of some genes may also function to increase the recombinogenicity of AFB1-induced DNA lesions. For example, the gene encoding protein kinase C (PKC1) is downregulated (3.9-fold decreased); null pkc1 mutants are inviable and arrest during S phase, whereas other mutations in PKC1 confer a hyper-recombinogenic phenotype (Huang and Symington, 1994
). Mutations in two other AFB1 downregulated signaling genes, MCM3 and SPT6, also confer hyper-recombinogenic phenotypes (Aguilera et al., 2000
). These data provide further evidence that the downregulation as well as the upregulation of specific genes may contribute to the recombinogenic cellular response to AFB1.
The mechanism by which the AFB1-induced changes in gene expression increase recombination is unknown. However, these changes may aid in identifying genes that contribute to the genotoxicity of AFB1, compared with other DNA-damaging agents. For example, AFB1-associated recombination depends on the function of several AFB1-inducible genes, such as RAD1 and RAD51. We demonstrated that higher levels of RAD51 message correlates with higher frequencies of AFB1-associated translocations in checkpoint mutants deficient in RAD51 induction. In mammalian cells overexpression of RAD51 also increases the frequency of chromosomal rearrangements and translocations (Richardson et al., 2004
). Thus, an increase in RAD51 levels in particular mammalian or yeast cells may increase recombination. Considering that mutations in upstream regulatory regions of other DNA damage-inducible genes, such as RAD54, do not confer a decrease in either radiation resistance or recombination (Cole and Mortimer, 1989
), further experiments are necessary to understand whether RAD51 induction per se is required for AFB1-associated translocations.
We had previously observed that rad51 diploid mutants exhibit 30- and 10-fold higher frequencies of translocations after x-ray and UV exposure, respectively (Fasullo et al., 2001
), whereas, rad1 mutants exhibit decreased frequencies of x-rayassociated translocations (Fasullo et al., 1998
), compared with wild type. Higher frequencies of x-rayassociated translocations were detected in rad51 mutants even when survival was low (Fasullo et al., 2001
). Thus, it is unlikely that AFB1-induced lethality caused the recombination defect. Most His+ recombinants generated in the rad51 mutant contained nonreciprocal translocations; whereas the majority of translocations identified after AFB1 and UV exposure in wild-type strains are reciprocal translocations. Nonreciprocal translocations may occur by break-induced replication (BIR) when a DNA polymerase replicates past a single-strand nick or when chromosomal fragments are inherited in subsequent divisions (Fasullo et al., 1998
). Considering that we are unable to detect chromosomal fragments after AFB1 exposure, we speculate that AFB1 lesions do not trigger replication fork collapse.
RAD1 and RAD51 play a different function in EMS or BaP-DHDassociated recombination (Figure 6). Although DNA damage generated by alkylating agents, such as EMS, is mainly repaired by BER (Friedberg et al., 1995
), DNA damage that results from agents that form bulky adducts, such as BaP-DHD, is mainly repaired by NER (Hess et al., 1997
). rad1 mutants are more sensitive to AFB1 than to EMS and are extremely sensitive to BaP-DHD; this suggests that NER is not the main pathway or is redundant in the repair of AFB1 lesions. We speculate that AFB1-induced lesions may require RAD1 to either initiate or process a recombination intermediate; RAD1-dependent recombination pathways have been extensively demonstrated by different groups (for review see Aguilera et al., 2000
; Haber, 2000
; Sung et al., 2000
) and may participate in crossing-over (Symington et al., 2000
). Interestingly, BaP-DHD is more recombinogenic in rad1 mutants than in wild-type strains, suggesting that more recombination events occur when the bulky adduct is not excised.
|
Similarly, rad51 mutants are more sensitive to EMS than to AFB1, indicating that RAD51 is involved in repair of AFB1 lesions, but more important in the repair of EMS lesions. The stimulation of recombination by particular alkylating agents, such as EMS, is dependent on cell division, suggesting that EMS-associated recombination occurs after the DNA polymerase encounters the unrepaired lesion (Galli and Schiestl, 1999
; Aguilera et al., 2000
). The hypersensitivity of the rad51 strain to EMS likely results from the rad51 defect in double-strand break repair; we speculate that double-strand breaks could be generated during BER if a DNA polymerase transverses a single-strand nick or gap and could stimulate a BIR mechanism. Additional experiments would be necessary to demonstrate that BER is a mechanism for generating more EMS-induced translocations in rad51 mutants.
DNA damage and the subsequent repair are thought to account largely for the carcinogenicity of AFB1. Recombinogenicity of a genotoxin could be pivotal in carcinogenesis as demonstrated by negative results in several mutagenesis tests (Schiestl, 1989
; Galli and Schiestl, 1995
, 1998
). Results shown here help elucidate mechanisms by which changes in gene expression contribute to the genotoxicity of a compound. It will be interesting to investigate whether AFB1 also changes the gene expression of orthologous genes in mammalian cells.
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
| Footnotes |
|---|
Corresponding author. E-mail address: mfasullo{at}ordwayresearch.org.
| REFERENCES |
|---|
|
|
|---|
Bailey, G.S. ((1994). ). Role of aflatoxin-DNA adducts in the cancer process. In: The Toxicology of Aflatoxins: Human Health, Veterinary, and Agricultural Significance, ed. D.L. Eaton and J.D. Groopmann, San Diego: Academic Press, 137148.
Buss, P., Caviezel, M., and Lutz, W.K. ((1990). ). Linear dose-response relationship for DNA adducts in rat liver from chronic exposure to aflatoxin B1. Carcinogenesis 11, , 21332135.
Cole, G.M., and Mortimer, R.K. ((1989). ). Failure to induce a DNA repair gene, RAD54, in Saccharomyces cerevisiae does not affect DNA repair or recombination phenotypes. Mol. Cell. Biol. 9, , 33143322.
Croy, R.G., Essigmann, J.M., Reinhold, V.N., and Wogan, G.N. ((1978). ). Identification of the principal aflatoxin B1-DNA adduct formed in vivo in rat liver. Proc. Natl. Acad. Sci. USA 75, , 17451749.
Croy, R.G., and Wogan, G.N. ((1981). ). Quantitative comparison of covalent aflatoxin-DNA adducts formed in rat and mouse livers and kidneys. J. Natl. Cancer Inst. 66, , 761768.
Davies, A.A., Friedberg, E.C., Tomkinson, A.E., Wood, R.D., and West, S.C. ((1995). ). Role of the Rad1 and Rad10 proteins in nucleotide excision repair and recombination. J. Biol. Chem. 270, , 2463824641.
Dennis, P.B., Fumagalli, S., and Thomas, G. ((1999). ). Target of rapamycin (TOR): balancing the opposing forces of protein synthesis and degradation. Curr. Opin. Genet. Dev. 9, , 4954.[CrossRef][Medline]
Dong, Z., and Fasullo, M. ((2003). ). Multiple recombination pathways for sister chromatid exchange in Saccharomyces cerevisiae: role of RAD1 and the RAD52 epistasis group genes. Nucleic Acids Res. 10, , 25762585.
Eaton, D.L., and Gallagher, E.P. ((1994). ). Mechanisms of aflatoxin carcinogenesis. Annu. Rev. Pharmacol. Toxicol. 34, , 135172.[CrossRef][Medline]
Elledge, S.J., and Davis, R.W. ((1990). ). Two genes differentially regulated in the cell cycle and by DNA damaging agents encode alternative regulatory subunits of ribonucleotide reductase. Genes Dev. 4, , 740751.
Essigmann, J.M., Croy, R.G., Nadzan, A.M., Busby, W.F.J., Reinhold, V.N., Buchi, G., and Wogan, G.N. ((1977). ). Structural identification of the major DNA adduct formed by aflatoxin B1 in vitro. Proc. Natl. Acad. Sci. USA 74, , 18701874.
Etienne, W., Meyer, M.H., Peppers, J., and Meyer, R.A. ((2004). ). Comparison of mRNA gene expression by RT-PCR and DNA microarray. Biotechniques 36, , 618626.[Medline]
Eugster, H.P., Bärtsch, S., Würgler, F.E., and Sengstag, C. ((1992). ). Functional co-expression of human oxidoreductase and cytochrome P450 1A1 in Saccharomyces cerevisiae results in increased EROD activity. Biochem. Biophys. Res. Commun. 185, , 641647.[CrossRef][Medline]
Fasullo, M., and Dave, P. ((1994). ). Mating type regulates the radiation-associated stimulation of reciprocal translocation events in Saccharomyces cerevisiae. Mol. Gen. Genet. 243, , 6370.[CrossRef][Medline]
Fasullo, M.T., and Davis, R.W. ((1987). ). Recombination substrates designed to study recombination between unique and repetitive sequences in vivo. Proc. Natl. Acad. Sci. USA 84, , 62156219.
Fasullo, M.T., Bennett, T., AhChing, P., and Koudelik, J. ((1998). ). The Saccharomyces cerevisiae RAD9 checkpoint reduces the DNA damage-associated stimulation of directed reciprocal translocations. Mol. Cell. Biol. 18, , 11902000.
Fasullo, M.T., Giallanza, P., Bennett, T., Cera, C., and Dong, Z. ((2001). ). Saccharomyces cerevisiae rad51 mutants are defective in DNA damage-stimulated sister chromatid exchange but exhibit increase rates of homology-directed translocations. Genetics 158, , 959972.
Foury, F. ((1997). ). Human genetic diseases: a cross-talk between man and yeast. Gene 195, , 110.[CrossRef][Medline]
Friedberg, E.C., Walker, G.C., and Siede, W. ((1995). ). DNA Repair and Mutagenesis, Washington, DC: ASM Press.
Galli, A., and Schiestl, R.H. ((1995). ). Salmonella test positive and negative carcinogens show different effects on intrachromosomal recombination in G2 cell cycle arrested yeast cells. Carcinogenesis 16, , 659663.
Galli, A., and Schiestl, R.H. ((1998). ). Effect of Salmonella assay negative and positive carcinogens on intrachromosomal recombination in S-phase arrested yeast cells. Mutat. Res. 419, , 5368.[Medline]
Galli, A., and Schiestl, R.H. ((1999). ). Cell division transforms mutagenic lesions into deletion-recombinagenic lesions in yeast cells. Mutat. Res. 429, , 1326.[Medline]
Gasch, A., Spellman, P., Kao, C., Carmel-Harel, O., Eisen, M., Storz, G., Botstein, D., and Brown, P.O. ((2000). ). Genomic expression programs in the response of yeast cells to environmental changes. Mol. Biol. Cell 11, , 42414257.
Gasch, A., Huang, M., Metzner, S., Botstein, D., Elledge, S., and Brown, P. ((2001). ). Genomic expression responses to DNA-damaging agents and the regulatory role of the yeast ATR homolog Mec1p. Mol. Biol. Cell 12, , 29873003.
Goffeau, A. et al. ((1996). ). Life with 6000 genes. Science 274, , 563547.
Goldstein, A.L., and McCusker, J.H. ((1999). ). Three new dominant drug resistance cassettes for gene disruption in Saccharomyces cerevisiae. Yeast 15, , 15411553.[CrossRef][Medline]
Haber, J.E. ((2000). ). Recombination: a frank view of exchanges and vice versa. Curr. Opin. Cell. Biol. 12, , 286292.[CrossRef][Medline]
Hess, M.T., Gunz, D., Luneva, N., Geacintov, N.E., and Naegeli, H. ((1997). ). Base pair conformation-dependent excision of benzo[a]pyrene diol epoxideguanine adducts by human nucleotide excision repair enzymes. Mol. Cell. Biol. 17, , 70697076.[Abstract]
Hsu, I.C., Metcalf, R.A., Sun, T., Welsh, J.A., Wang, N.J., and Harris, C.C. ((1991). ). Mutational hotspot in the p53 gene in human hepatocellular carcinomas. Nature 350, , 427428.[CrossRef][Medline]
Huang, K., and Symington, L. ((1994). ). Mutation of the gene encoding protein kinase C1 stimulates mitotic recombination in Saccharomyces cerevisiae. Mol. Cell. Biol. 14, , 60396045.
Jelinsky, S.A., and Samson, L.D. ((1999). ). Global response of Saccharomyces cerevisiae to an alkylating agent. Proc. Natl. Acad. Sci. USA 96, , 14861491.
Keith, C.T., and Schreiber, S.L. ((1995). ). PIK-related kinases: DNA repair, recombination, and cell cycle checkpoints. Science 270, , 5051.
Keller-Seitz, M. ((2001). ). PhD thesis no. 14321. Swiss Federal Institute of Technology (ETH).
Klebe, R.J., Harriss, J.V., Sharp, Z.D., and Douglas, M.G. ((1983). ). A general method for polyethylene-glycol-induced genetic transformation of bacteria and yeast. Gene 25, , 333341.[CrossRef][Medline]
Klein, H.L. ((2000). ). A radical solution to death. Nat. Genet. 25, , 132134.[CrossRef][Medline]
Leadon, S.A., Tyrrell, R.M., and Cerutti, P.A. ((1981). ). Excision repair of aflatoxin B1-DNA adducts in human fibroblasts. Cancer Res. 41, , 51255129.[Medline]
Lin, J.-K., Miller, J.A., and Miller, E.C. ((1977). ). 2,3-Dihydro-2-(guan-7-yl)-3-hydroxy-aflatoxin B1, a major acid hydrolysis product of aflatoxin B1-DNA or -ribosomal RNA adducts formed in hepatic microsome-mediated reactions in rat liver in vivo. Cancer Res. 37, , 44304438.
Martin, C.N., and Garner, R.C. ((1977). ). Aflatoxin B1-oxide generated by chemical or enzymatic oxidation of aflatoxin B1 causes guanine substitution in nucleic acids. Nature 267, , 863865.[CrossRef][Medline]
Mercier, G., Denis, Y., Marc, P., Picard, L., and Dutriex, M. ((2001). ). Transcriptional induction of repair genes during slowing of replication in irradiated Saccharomyces cerevisiae. Mutat. Res. 487, , 157172.[Medline]
Mewes, H.W., Albermann, K., Heumann, K., Liebl, S., and Pfeiffer, F. ((1997). ). MIPS: a database for protein sequences, homology data and yeast genome information. Nucleic Acids Res. 25, , 2830.
Nicholson, A., Hendrix, M., Jinks-Robertson, S., and Crouse, G.F. ((2000). ). Regulation of mitotic homeologous recombination in yeast. Functions of mismatch repair and nucleotide excision repair genes. Genetics 154, , 133146.
Otteneder, M., and Lutz, W.K. ((1999). ). Correlation of DNA adduct levels with tumor incidence: carcinogenic potency of DNA adducts. Mutat. Res. 424, , 237247.[Medline]
Paques, F., and Haber, J.E. ((1999). ). Multiple pathways of recombination induced by double-strand breaks in Saccharomyces cerevisiae. Microbiol. Mol. Biol. Res. 63, , 349404.
Payne, W.E., and Garrels, J.I. ((1997). ). Yeast Protein database (YPD): a database for the complete proteome of Saccharomyces cerevisiae. Nucleic Acids Res. 25, , 5762.
Pedrazzi, G., et al. ((2001). ). Direct association of Bloom's syndrome gene product with the human mismatch repair protein. MLH1. Nucleic Acids Res. 29, , 43784386.
Posas, F., Chambers, J.R., Heyman, J.A., Hoeffler, J.P., de Nadal, E., and Arino, J. ((2000). ). The transcriptional response of yeast to saline stress. J. Biol. Chem. 275, , 1724917255.