| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Molecular Biology, Pathobiology and Genetics |
Laboratory of Computational Biology and Risk Analysis, National Institute of Environmental Health Sciences, NIH, Research Triangle Park, North Carolina
Request for reprints: Douglas A. Bell, Environmental Genomics Section, C3-03, P.O. Box 12233, National Institute of Environmental Health Sciences, NIH, Research Triangle Park, NC 27709. Phone: 919-541-7686; Fax: 919-541-4634. E-mail: bell1{at}niehs.nih.gov.
| Abstract |
|---|
|
|
|---|
Key Words: SNPs gene expression glioma ChIP GSTM3 promoter
| Introduction |
|---|
|
|
|---|
GSTM3 is expressed in testis, brain, lung, and lymphocytes as well as in many tumor types (NCI-60; refs. 1719). A recently described GSTM3 3-bp deletion polymorphism (GSTM3*A/*B) has been shown to have a frequency of 0.16 in the Caucasian population (20, 21). This polymorphism was associated with an increased risk of cutaneous basal cell carcinomas and laryngeal carcinomas and interacted with other GST variants in modifying astrocytoma risk but had no effect on pharyngeal carcinoma risk (2024). Anttila et al. (17) observed variation in GSTM3 expression in lung tissue and proposed a relationship among smoking exposure, GSTM3 expression level, and GSTM1-null genotype.
To explore the mechanism for the Anttila et al. observations of GSTM3 expression variation, we identified human cell lines with low expression GSTM3 phenotypes in a Centre d'Etude du Polymorphisme Humain (CEPH) family and sequenced the GSTM3 promoters of these individuals. Using expression analyses (quantitative real-time PCR and Western blot), luciferase reporter assays and allele-specific chromatin immunoprecipitation (ChIP), we have shown that a GSTM3 promoter polymorphism causes a dramatic reduction in GSTM3 transcription and protein expression in the lymphoblast and glioma cell lines we investigated.
| Materials and Methods |
|---|
|
|
|---|
DNA Extraction, Genotyping, and Sequencing. For genotyping, we used a community-based sample of healthy, unrelated (19-53 years old) volunteers of African and European descent from Durham and Chapel Hill, NC (protocol 86E-0037), and placental DNA samples from anonymous full-term, normal, unrelated maternity patients with Asian ancestry were provided by Dr. L.L. Hsieh (Chang Gung University, Taiwan, China). DNA was isolated from all samples with standard proteinase K digestion followed by phenol-chloroform extraction and ethanol precipitation. Genotyping of the promoter single nucleotide polymorphism (SNP; rs1332018) was carried out using a Taqman allelic discrimination assay (Applied Biosystems, Foster City, CA). The primers used were forward: 5'-GCGGAGGCGGAGTCA-3' and reverse: 5'-TCGAGGACTAGGGAAACTGTG-3'. The probes used were VIC: 5'-CCCCCTTATGTAGGGTATA-3' and FAM: 5'-CCCTTATGTCGGGTATA-3'. The reaction mix included 20 ng genomic DNA, 1x primer (0.2 µmol/L)/probe (0.9 µmol/L) mix, and 1 x Taqman Universal Master Mix (Applied Biosystems) in a final volume of 10 µL. Thermocycling consisted of a 95°C, 10-minute hold for AmpliTaq Gold activation followed by 40 cycles of 95°C, 15 seconds and 60°C, 60 seconds. VIC and FAM probe signals were recorded post-PCR using an ABI PRISM 7700 instrument (Applied Biosystems). GSTM3 promoter sequencing in lymphoblasts was done using an ABI PRISM BigDye Terminator Cycle Sequencing kit and an ABI PRISM 377 DNA sequencer (Applied Biosystems). Sequences were compared with the GSTM3 consensus (Genbank accession no. AF043105) using the PileUp program of GCG Wisconsin package version 10.3 (Accelrys, San Diego, CA).
Measurement of GSTM3 Expression Using Real-time PCR. Total RNA was isolated from cells using a RNeasy Mini kit (Qiagen, Inc., Valencia, CA) treated with DNase and measured for A260/A280 nm ratios using a NanoDrop spectrometer (NanoDrop, Rockland, ME). cDNA was generated from 1 µg total RNA using the SuperScript First-Strand Synthesis System (Invitrogen). The cDNA reaction solution was diluted to a final concentration of 1.33 ng/µL. Real-time PCR constituents were 22.5 µL DNA, 1x Taqman Universal Master Mix, and 1x primer/probe mix (Hs00356079, Applied Biosystems) in a 50 µL final volume. Real-time PCR amplification was done using the same conditions as for genotyping, except for the addition of an initial 50°C, 2-minute hold for optimal AmpErase UNG activity. Signals were recorded during PCR using an ABI PRISM 7700 instrument. All gene expression results were normalized to HPRT expression. Gene expression was calculated using the
Ct method (Applied Biosystems).
Measurement of GSTM3 Protein Level Using Western Blot. Cells were harvested and incubated (30 minutes on ice) in lysis buffer. After centrifugation (13,000 rpm, 15 minutes), proteins in the supernatants were quantified (Bio-Rad Laboratories, Hercules, CA). Cell lysates (10 µg) were denatured (80°C, 15 minutes), loaded onto an 18% Tris-glycine gel (Invitrogen), run on a XCell Surelock Mini-Cell (Invitrogen), and transferred onto a 0.2 µm nitrocellulose membrane (Invitrogen). The membrane was rinsed (PBS with 0.05% Tween), blocked (1% bovine serum albumin in PBS), rinsed again, and incubated in primary rabbit anti-human GSTM3 antibody (1:2,000 dilution) generously provided by Dr. John D. Hayes (Biomedical Research Centre, University of Dundee, Dundee, United Kingdom). Protein was detected by enhanced chemiluminescence (Amersham Pharmacia Biotech, Little Chalfont, United Kingdom) autoradiography. The membrane was stripped and similarly blotted with ß-actin antibody.
Preparation of pGL3 Vector Constructs. A fragment of the GSTM3 5' untranslated region, 948 to +3 nucleotides (relative to codon 1 start, Genbank accession no. AF043105), was generated by PCR amplification using genomic DNA from GM07354 cells (783AA/63AA). PCR was done using Expand High-Fidelity PCR System (Roche Molecular Diagnostics Corp., Indianapolis, IN) according to the manufacturer's instructions. The primers used were forward: 5'-AAGAAGCTCGAGGTTGGACAGATGAATGAGTG-3' and reverse: 5'-AAGAAGAGATCTCATGGTGACGGGCTTCCGAG-3'. The PCR product was purified and inserted into the firefly luciferase reporter plasmid, pGL3-enhancer vector (Promega, Madison, WI), resulting in the reporter plasmid. In pGL3, the four possible haplotypes for the A783G and A63C sites were generated using the QuickChange Site-Directed Mutagenesis kit (Stratagene, La Jolla, CA). The primers used were forward: 5'-CCGCCCCCTTATGTgGGGTATAAAGCCCCTC-3' and reverse: 5'-GGAGGGGCTTTATACCCcACATAAGGGGGCG-3' for A783G site and forward: 5'-GATTAATTTCCCTAAATGCcTGAGCATTAGG-3' and reverse: 5'-CCTAATGCTCAgGCATTTAGGGAAATTAATC-3' for A63C site. Sequences were verified using the dye terminator cycle sequencing method on an ABI PRISM.
Transient Transfection. The human glioma cell line SF-268 was plated in six-well plates at a density of 2 x 105 cells per well and incubated at 37°C in RPMI 1640 containing 15% fetal bovine serum and 1% penicillin-streptomycin (Invitrogen) for 24 hours. After washing the cells with PBS, the cells were cotransfected (37°C, 4 hours) with 3 µg pGL3-enhancer construct, 0.3 µg pRL-TK (for normalization of luciferase expression), and 6 µL Transfectam (Promega) in 500 µL RPMI 1640 without fetal bovine serum and then incubated (37°C, 48 hours) with fresh medium. The pRL-TK, pGL-3 control, and pGL-3 empty vectors were used to establish transfection conditions. Cells were washed with PBS twice and harvested after incubation (2 hours, room temperature) in 500 µL passive lysis buffer (Promega) with agitation. Luciferase activity was determined using 20 µL cell lysate aliquots. The lysates were mixed with 100 µL luciferase firefly assay reagent and read on a model TD-20/20 luminometer (Turner Designs, Sunnyvale, CA) every 0.5 second for 25 seconds in a 96-well plate. Luciferase activity was normalized to the intensity of pRL-TK (Promega) and expressed as relative light units. All luciferase assays were carried out in triplicate.
ChIP Assay. We have used an approach similar to Knight et al. (25), combining a ChIP protocol in heterozygous 63AC (SNB-19) cells with Taqman quantitative genotyping. SNB-19 cells (3 x 107) were used in 10 independent experiments carried out with a ChIP procedure modified from Upstate Biotechnology (Lake Placid, NY). Cells were fixed using 1% formaldehyde and incubated (37°C, 10 minutes). Fixation was stopped with 125 mmol/L glycine (37°C, 5 minutes). Cells were washed, pelleted (4°C, 2,000 rpm, 4 minutes), and lysed with 200 µL SDS lysis buffer (Upstate Biotechnology) containing the protease inhibitors. The lysate was sonicated (20 seconds, six times) using a Sonicator 3000 (Misonix, Farmingdale, NY) at 90% of output level, resulting in an average of 500-bp DNA length, determined by 1% agarose Tris-borate EDTA gel electrophoresis. Samples were centrifuged (13,000 rpm, 10 minutes) and the supernatant was diluted 10-fold in ChIP dilution buffer (Upstate Biotechnology) to a 2 mL final volume. For each sample, 20 µL supernatant was kept as input/starting material. Immunoprecipitations were done by incubating (4°C, 30 minutes) 2 mL diluted supernatant in 80 µL salmon sperm DNA/protein A agarose 50% slurry with agitation and then by centrifugation (2,000 rpm, 3 minutes). The supernatant was incubated (4°C, overnight) in 3 µg RNA pol II antibody (Santa Cruz Biotechnology, Santa Cruz, CA). A no-antibody immunoprecipitation served as a negative control. Following precipitation, protein-DNA complexes were washed and eluted from the antibody by adding elution buffer (Upstate Biotechnology). Protein-DNA cross-links were reversed in 0.2 mol/L NaCl by heating (65°C, 4 hours). The protein was digested in 40 µg/mL proteinase K (45°C, 1 hour, Qiagen), 0.01 EDTA, and 0.40 mol/L Tris-HCl (pH 6.5). DNA was extracted with phenol-chloroform, precipitated with ice-cold 100% ethanol, and resuspended in 50 µL Tris-EDTA. DNA concentrations were estimated using a NanoDrop method according to the manufacturer's directions.
Quantitative Allelic Discrimination Assay. Quantitative allelic discrimination of the A63C SNP was done using an ABI PRISM 7700 with the same primer and probe set used for genotyping. A standard curve was created using 63AA and 63CC genomic DNA mixtures at varying percentages (0%, 10%, 20%, ..., 100%). PCR constituents and thermocycling were the same as for A63C genotyping, except that 11.25 µL immunoprecipitated DNA (estimated to be
10 ng) in a 25 µL final volume and 36 cycles were used. VIC and FAM signals were recorded post-PCR using an ABI PRISM 7700.
Statistical Analysis. The results for gene expression, transient transfection, and allele-specific ChIP assays are given as mean ± SD of at least triplicate measurements. The significance of the results was determined using the Student's t test at
= 0.05.
| Results and Discussion |
|---|
|
|
|---|
We screened GSTM1 genotype and GSTM3 expression in lymphoblast cell lines from a CEPH family (CEPH1345) containing parents heterozygous for the GSTM1*1/*0 null genotype and children with various combinations of the GSTM1 alleles. These individuals were all homozygous for the GSTM3*A allele (the GSTM3*B allele is a 3-nucleotide deletion in intron 6). No relationship between GSTM1 genotype and GSTM3 expression was observed1; however, several children displayed GSTM3 low expression phenotypes (Fig. 1). GSTM3 promoter sequencing in these children revealed two SNPs, A783G and A63C, relative to the start of codon 1. The individuals homozygous for the 63CC genotype had 8-fold lower expression compared with those with the 63AA genotype (P < 0.003). Mean GSTM3 expression in the two parents and two other children, all 63AC heterozygotes, was 1.6-fold lower than the 63AA children, suggesting a recessive mode of action for this SNP. GM07350 displays somewhat lower expression than the other three heterozygous cell lines, possibly due to the presence of the 783 G allele. However, this was not consistent across all cell lines.
|
50 kb) of GSTM3 to the GSTM1*0/*1 deletion polymorphism, the GSTM3 63C allele can occur on the same chromosome with either GSTM1*1 or GSTM1*0 allele (all four haplotypes are present), suggesting that these polymorphisms are not strongly in linkage.
|
|
900 bp of the GSTM3 promoter (containing the 783A/63A haplotype) and cloned this fragment into a pGL3-enhancer vector upstream of the luciferase reporter. After transient transfection into SF-268 cells, strong luciferase expression was observed for the 783A/63A vector compared with an empty pGL3-enhancer vector (Fig. 3). Using in vitro mutagenesis, we generated luciferase reporter constructs containing the four possible 783A and 63A haplotypes and measured luciferase expression activity in SF-268 cells (Fig. 3). Relative to the 783A/63A allele, both 63C containing alleles reduced expression. The 783A/63C allele reduced luciferase expression
4-fold from 9.4 ± 1.36 to 2.36 ± 0.3 (P < 0.0004) and the 783G/63C allele was reduced 2.7-fold. The 783G polymorphism also reduced luciferase expression but to a lesser extent (1.6-fold lower in the 783G/63A haplotype relative to the 783A/63A haplotype; P < 0.02). These observations indicate that the A63C site is more important for promoter activity than the A783G site. We hypothesized that the A63C site, which is 4-nucleotide upstream of the TATA box, could alter the binding of RNA pol II and examined this possibility using ChIP.
|
We validated the quantitative genotyping component by measuring the VIC (63A) and FAM (63C) probe signals in cell line DNA with the three possible A63C genotypes (63A/A, 63C/C, and 63A/C) and in genomic DNA mixtures of varying 63A/A and 63C/C DNA percentages (0%, 10%, 20%, ..., 100%). The DNA mixture experiment showed a linear relationship (r2 = 0.93) between allele composition and VIC/FAM ratios (Fig. 4A). VIC/FAM ratios for the homozygous 63A/A and 63C/C DNA samples were 2.82 ± 0.41 and 0.28 ± 0.03 (Fig. 4B), respectively, and the heterozygous 63A/C cell line was 0.93 ± 0.02. Furthermore, the VIC/FAM ratio observed from a heterozygous 63A/C cell line, GM01032A (0.93 ± 0.02), was not significantly different from the 50%:50% allele mixture experiment (1.17 ± 0.04). VIC and FAM signals were linear down to at least 1 ng genomic DNA per reaction (data not shown). This low detection limit and low variability indicated that this method would permit quantitative allelic discrimination in small quantities of immunoprecipitated DNA and would allow multiple analyses from the same RNA pol IIprecipitated sample.
|
To determine the relative binding of RNA polymerase to the 63A and 63C promoter alleles, we quantified the presence of each of the alleles in the DNA recovered from the ChIP procedure. Figure 4C displays the mean VIC/FAM ratios obtained from quantitative genotyping (in triplicate) of the unprecipitated cell lysate reference sample (left) and the DNA recovered from 10 parallel RNA pol II ChIP experiments (mean value is plotted). The mean VIC/FAM ratio for immunoprecipitated DNA was 2.5 ± 0.3 and this difference (from reference value, 1.04 ± 0.02) was highly significant (P < 0.0001). Using the mixed DNA samples as a standard curve, we estimated that the precipitated DNA contained at least 90% of A allele and 10% of C allele (compared with the 90% A:10% C DNA mixture; VIC/FAM ratio = 2.14). SNB-19 cell extracts mixed with protein A agarose beads with no antibody present showed baseline VIC and FAM signal levels (data not shown), indicating an absence of nonspecifically precipitated GSTM3 promoter DNA. These allele-specific ChIP results are consistent with both gene expression (Figs. 1 and 2) and luciferase reporter assays results (Fig. 3). It is interesting to note that the magnitude of the allele difference for the allele-specific ChIP assays and the mRNA expression assays are similar (
9-fold difference), although results for the luciferase reporter assays showed a smaller difference (
4-fold).
Allele-specific ChIP is a powerful technique for assessing functional differences between promoter alleles. This method provides in vivo relative determination of transcription factor loading on an endogenous gene promoter under normal regulatory control in a chromosomal environment. This approach may reflect the "true" biological differences between promoter alleles in vivo and as such could prove to be more biologically relevant than the usual allele comparisons using dual luciferase reporter gene constructs. A previous study by Knight et al. (25) determined relative RNA pol II binding to tumor necrosis factor
promoter by combining ChIP with primer extension/matrix-assisted laser desorption/ionization time-of-flight mass spectrometrybased genotyping (Sequenom, San Diego, CA). In the present study, we modified the Knight et al. approach by combining RNA pol II ChIP with allele quantification using a 5' nuclease-based allelic discrimination assay (Taqman). As there are many thousands of laboratories with instruments capable of using the 5' nuclease chemistry, this approach should be more accessible than the matrix-assisted laser desorption/ionization time-of-flight mass spectrometry method.
Biological Impact of the A63C Polymorphism in the GSTM3 Promoter. GSTM3 is the ancestral gene in the GSTµ gene family (2), and several biological functions have been shown. Like other GSTµ proteins, GSTM3 catalyzes the conjugation of numerous toxic or carcinogenic compounds to glutathione, and low levels of GSTM3 expression in the cell might lead to greater amounts of cellular DNA damage. It has also been suggested that GSTM3 expression could affect sensitivity to chemotherapeutic agents. In addition, GSTM3 was identified recently as a cytosolic glutathione-dependent prostaglandin E2 synthase in human brain (34), suggesting an endo genous role for GSTM3 protein levels in regulation of prostaglandin signaling pathways. Because prostaglandin E2 plays an important role in cell proliferation signaling, GSTM3 expression variation could have an impact on growth or cell proliferation during tumorigenesis in the brain or perhaps in other GSTM3-expressing tissues.
In summary, we have shown that the GSTM3 A63C SNP has a major impact on regulation of GSTM3 expression in vivo. GSTM3 has several potentially important biological functions, and the effect and prevalence of a SNP in multiple ethnic populations suggests the possibility of a significant role in carcinogenesis or tumor biology. To test this hypothesis, we have initiated studies examining the relationship between this GSTM3 SNP and the risk of brain, lung, and testicular cancer.
| Acknowledgments |
|---|
We thank Dr. L.L. Hsieh for providing DNA samples from Taiwanese maternity patients, Prof. Julian Knight (University of Oxford, Oxford, United Kingdom) for kindly sharing his ChIP protocol, Dr. John Hayes for sharing his GSTM3 antibody, and the National Cancer Institute Developmental Therapeutics Program and Susan Holbeck (National Cancer Institute) for tumor cell lines and helpful advice.
| Footnotes |
|---|
Received 8/ 9/04. Revised 10/ 6/04. Accepted 10/26/04.
| References |
|---|
|
|
|---|
-class glutathione S-transferase and D-dopachrome tautomerase gene complex. Biochem J 1998;334:61723.
glutathione S-transferase T1-1 can catalyse the activation of dichloromethane, a liver and lung carcinogen in the mouse. Comparison of the tissue distribution of GST T1-1 with that of classes
, µ and
GST in human. Biochem J 1997;326:83746.
This article has been cited by other articles:
![]() |
J. A. Schwartzbaum, A. Ahlbom, S. Lonn, M. Warholm, A. Rannug, A. Auvinen, H. C. Christensen, R. Henriksson, C. Johansen, C. Lindholm, et al. An International Case-Control Study of Glutathione Transferase and Functionally Related Polymorphisms and Risk of Primary Adult Brain Tumors Cancer Epidemiol. Biomarkers Prev., March 1, 2007; 16(3): 559 - 565. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |