- Research article
- Open Access
- Open Peer Review
GSTT2 promoter polymorphisms and colorectal cancer risk
BMC Cancervolume 7, Article number: 16 (2007)
Glutathione S-transferases are a group of enzymes that participate in detoxification and defense mechanisms against toxic carcinogens and other compounds. These enzymes play an important role in human carcinogenesis. In the present study, we sought to determine whether GSTT2 promoter single nucleotide polymorphisms (SNPs) are associated with colorectal cancer risk.
A total of 436 colorectal cancer patients and 568 healthy controls were genotyped for three GSTT2 promoter SNPs (-537G>A, -277T>C and -158G>A), using real-time TaqMan assay and direct sequencing. An electrophoretic mobility shift assay (EMSA) was performed to determine the effects of polymorphisms on protein binding to the GSTT2 promoter.
The -537A allele (-537G/A or A/A) was significantly associated with colorectal cancer risk (OR = 1.373, p = 0.025), while the -158A allele (-158G/A or A/A) was involved in protection against colorectal cancer (OR = 0.539, p = 0.032). Haplotype 2 (-537A, -277T, -158G) was significantly associated with colorectal cancer risk (OR = 1.386, p = 0.021), while haplotype 4 (-537G, -277C, -158A) protected against colorectal cancer (OR = 0.539, p = 0.032). EMSA data revealed lower promoter binding activity in the -537A allele than its -537G counterpart.
Our results collectively suggest that SNPs and haplotypes of the GSTT2 promoter region are associated with colorectal cancer risk in the Korean population.
Colorectal cancer (CRC), a predominant cause of tumor-related death in Western nations, is becoming more prevalent in Asian countries such as Korea [1, 2]. Genetic and environmental factors are modulators of the carcinogenesis process . Studies on the contribution of genetic polymorphisms to cancer development have examined numerous genes, including oncogenes, tumor suppressors, DNA repair genes and those encoding Phase I and Phase II enzymes . Phase II enzymes attach to additional substrates in an attempt to detoxify the activated metabolite in preparation for final breakdown or excretion. Glutathione-S-transferases (GSTs) are a family of Phase II detoxification enzymes that protect cellular macromolecules by catalyzing the conjugation of glutathione (GSH) to a wide variety of endogenous and exogenous electrophilic compounds . GSTs are divided into at least six classes (alpha, mu, pi, theta, omega and zeta), each of which consists of one or more isoforms . GST polymorphisms are associated with bladder, colorectal, gastric, prostate and various other human cancers [5, 7–11].
The theta class consists of two genes, GSTT1 and GSTT2, located at 22q11.2, and separated by about 50 kb. GSTT1 and GSTT2 share 55% amino acid sequence identity . The GSTT1-null genotype is associated with increased CRC risk [7, 13, 14]. A common characteristic of the theta class is their affinity for the organic hydroperoxide species; for example, cumene hydroperoxide is a good substrate for GSTT2 . Thus, GSTT2 activity is important for the protection of cells against toxic products of oxygen and lipid peroxidation , which represents a major source of endogenous DNA damage in humans that contributes significantly to cancer and other genetic diseases . Pool-Zobel et al.  observed that GSTT2 involved in defense against oxidative stress in primary human colon cells is up-regulated upon incubation with butyrate. A comprehensive study of Phase I and Phase II metabolism gene polymorphisms revealed no significant contribution of GSTT2 polymorphism to CRC risk . However, distinct effects of GSST2 promoter haplotypes on gene expression were observed with a luciferase reporter assay . Genetic alterations may induce abnormal gene expression. To date, no case-control analysis of GSTT2 promoter polymorphisms has been performed in human cancers.
In the present study, we sought to determine whether four GSTT2 promoter polymorphisms, -537G>A, -277T>C, -158G>A, and -129T>C, and their haplotypes are associated with CRC risk. We additionally employed EMSA to determine whether GSTT2 promoter polymorphisms affect the binding of putative transcription factors.
Samples and DNA extraction
Four hundred and thirty-six colorectal cancer samples were collected from the Seoul National University Hospital and the National Cancer Center Hospital, Korea. Clinical characteristics, including age, sex, tumor location, and TNM stage, were additionally analyzed for the association study with GSTT2 genotyping. DNA was extracted from normal colorectal tissues. Fresh cancer tissues were stored at -70°C in a liquid nitrogen tank. Normal population samples were obtained from 568 healthy Korean individuals as controls, and DNA was isolated from blood samples. All colorectal cancer samples were collected from Korean patients enrolled between 1990 and 2003 at Seoul National University and National Cancer Center. Normal controls were selected from cancer-free samples enrolled from the Cancer Cohort Study Branch of the National Cancer Center. The mean ages of colorectal cancer patients at diagnosis and normal subjects were 58.8 years (58.8 ± 12.9 years) and 52.1 years (52.1 ± 10.5 years), respectively.
Total genomic DNA was extracted with TRIzol reagent, according to the manufacturer's instructions (Invitrogen, CA, USA). This study was approved by the institutional review board of the National Cancer Center and informed consent was obtained from all subjects prior to inclusion in the study.
The -537G>A (rs140186) and -277T>C (rs9624369) polymorphisms were screened using real-time TaqMan assay, and -158G>A and -129T>C polymorphisms were screened by direct sequencing. For the TaqMan assay, both PCR primers and MGB TaqMan probes were designed using the commercial Assay by Design service (Applied Biosystems, Foster City, CA, USA). One allelic probe was labeled with FAM dye, and the other with VIC. PCR was performed using the TaqMan Universal Master MIX without UNG in a 7900 HT fast real-time system (Applied Biosystems, Foster City, CA, USA). Data files were assessed using SDS2.1.1 software (Applied Biosystems, Foster City, CA, USA). We used BigDye Terminator and an ABI PRISM 3730 DNA analyzer (Applied Biosystems, Foster City, CA, USA) for direct sequencing analysis. Amplification for direct sequencing was performed in a final volume of 15 μl containing 10 ng genomic DNA, 10 pmol of each primer (forward: 5'-CCACTGGGTGAAACTCTAAG-3', reverse: 5'-GGACACCAGGTCAAGAAAC-3'), 0.25 mM each dNTP, 0.5 U of Taq polymerase and reaction buffer (Genecraft Ltd, Munster, Germany). Reactions were performed in a programmable thermal cycler (MWG Biotech AG, Eberberg, Germany) under the following conditions: denaturation for 5 min at 94°C, followed by 40 cycles at 94°C for 30 s, 55°C for 30 s, 72°C for 1 min, and a final extension of 10 min at 72°C.
Microsatellite instability (MSI) analysis
Two microsatellite markers (BAT-25 and BAT-26) were used to assess MSI status. DHPLC (denaturing high-performance liquid chromatography) and the capillary-based method were employed . Cancer samples were classified as MSI when at least one MSI was evident in two markers. MSI data were obtained from 200 previously reported samples  and 236 patients analyzed in this study.
Electrophoretic mobility shift assay (EMSA)
EMSA was performed with the Gel Shift Assay System (Promega Corporation, Madison, WI, USA). Briefly, the following consensus oligonucleotide pairs corresponding to the GSTT2 promoter sequence were synthesized (Bioneer, Seoul, South Korea) (bold letters specify polymorphisms): -537G (including the G allele) forward 5'-TAAGATCCCTTTTAGGGGATCCCATTCGCTC-3' and reverse 5'-GAGCGAATGGGATCCCCTAAAAGGGATCTTA-3', -537A (including the A allele) forward 5'-TAAGATCCCTTTTAGAGGATCCCATTCGCTC-3', and reverse 5'-GAGCGAATGGGATCCTCTAAAAGGGATCTTA-3'; -158G (including the G allele) forward 5'-AGGAACCGAAGGGGCGAGGCGGGTCCGGGGG-3' and reverse 5'-TCCTTGGCTTCCCCGCTCCGCCCAGGCCCCC-3', -158A (including the A allele) forward 5'-AGGAACCGAAGGGGCAAGGCGGGTCCGGGGG-3', and reverse 5'-TCCTTGGTTCCCCGTTCCGCCCAGGCCCCC-3'. Primer pairs were annealed and labeled with [γ-32P]ATP (Amersham Biosciences, Buckinghamshire, UK). Binding reactions were performed using HeLa nuclear extracts (Promega Corporation, Madison, WI, USA), according to the manufacturer's instructions. The 32P-labeled probe (1 μl) was incubated with 5 mg of HeLa nuclear extracts for 1 h at room temperature. Following the binding reaction, DNA-protein complexes were resolved by electrophoresis on a 4% non-denaturing acrylamide gel, which was subsequently transferred to 3 M blotting paper, dried, and exposed to film.
GSTT2 genotyping results were analyzed for categorical variables using the χ2 test, and for continuous variables using student's t-test (SPSS 12.0, Chicago, IL). We examined Lewontin's D' (|D'|) and the linkage disequilibrium coefficient, r 2, between all pairs of biallelic loci [21, 22]. Haplotypes of each individual were inferred with the algorithm developed by Stephens et al. [21, 22], which uses a Bayesian approach incorporating prior expectations of haplotypic structure based on population genetics and the coalescent theory (PHASE version 2.0). The genotype and haplotype-specific risks were estimated as odds ratios (OR) with associated 95% confidence intervals (CI) by logistic regression (SPSS 12.0, Chicago, IL). The observed genotype frequencies were analyzed with a chi-square test to determine whether they were in Hardy-Weinberg equilibrium (HWE). A p value of less than 0.05 was considered statistically significant.
Results and discussion
We initially selected three single nucleotide polymorphisms (SNPs) within the GSTT2 promoter region (-537G>A, -277T>C and -129T>C) highlighted in an earlier study  and further identified a novel -158G>A polymorphism. The rare -129T>C polymorphism with less than 1% minor allele frequency (MAF) was excluded, and three GSTT2 promoter polymorphisms (-537G>A, -277T>C and -158G>A) were analyzed for CRC risk. The MAF values of the three SNPs were 0.42 (-537A), 0.07 (-277C) and 0.03 (-158A) (Table 1). In control samples, the genotype distribution did not deviate from HWE (Table 1). The GSTT2 genotype results were not influenced by other potential CRC risk factors, such as age, sex, location, and TNM stage (Table 2). Moreover, the MSI status was not associated with GSTT2 polymorphisms (Table 2). Genotypes were analyzed for association with CRC risk using logistic regression models adjusted for age and sex (Table 1). The frequency of the -537A allele (G/A or A/A) was significantly increased in CRC patients compared to healthy controls (OR = 1.373, p = 0.025, Table 1). In contrast, the -158A allele (G/A or A/A) was more frequently observed in healthy controls than CRC patients (OR = 0.539, p = 0.032, Table 1). Moreover, the -277T/C polymorphism was not related to CRC risk (p = 0.125).
Linkage disequilibrium (LD) was tested for all SNP pairs. Lewontin's D' values of the pairs were 0.965 (-537G>A:-277T>C), 1.000 (-537G>A:-158G>A) and 1.000 (-277T>C:-158G>A), and the r 2 values were 0.043, 0.017, and 0.385, respectively. Good level of LD (D' > 0.965) was observed between all SNP pairs. Four haplotypes with more than 2% frequency were selected for haplotype association analysis. Haplotype frequencies were 0.511 (HT1), 0.422 (HT2), 0.033 (HT3), and 0.032 (HT4) (Table 3). Haplotypes were analyzed for association with CRC risk using logistic regression models adjusted for age and sex (Table 3). Haplotype 2 (HT2; -537A, -277T, -158G) was correlated with significantly increased CRC risk, compared to non-HT2 (OR = 1.386, p = 0.021, Table 3). In contrast, haplotype 4 (HT4; -537G, -277C, -158A) was associated with protection against CRC, compared with non-HT4 (OR = 0.539, p = 0.032, Table 3). No significant relationship with CRC risk was evident for haplotypes 1 (HT1; -537G, -277T, -158G) (p = 0.336) and 3 (HT3; -537G, -277C, -158G) (p = 0.984). The -537A polymorphism positively associated with CRC risk was only observed in HT2, while the -158A polymorphism with a protective role against cancer was present in HT4. Individual SNP analysis of -537G>A and -158G>A is almost consistent with HT2 and HT4, respectively. This may be due to limited individual SNP numbers.
After Bonferroni correction (the threshold of significance was 0.017, 3 polymorphisms were analyzed), the associated p values did not retain significance. However, since SNPs within the same gene displayed good level of LD, test statistics for the 3 polymorphisms were not independent, and the significance of association with CRC risk is noteworthy. However, further functional evidence is required to confirm our results that GSTT2 promoter polymorphisms are associated with CRC risk.
Accordingly, to determine whether the two polymorphisms affect the binding of putative transcription factors, EMSA was performed using allele-specific consensus oligonucleotide probes. Upon incubation of radiolabeled oligonucleotides specific for -537G and -537A in the presence of HeLa nuclear extracts, the former oligonucleotide displayed stronger band intensity than the latter (Lanes 2 and 8) (Figure 1). To confirm the presence of the DNA-protein complex, competition assays were performed with increasing amounts of unlabeled oligonucleotides (10-, 50-, 100-fold excess). Band densities decreased with increasing concentrations of unlabeled specific competitor (Lanes 3–5 and Lanes 9–11). In the presence of unlabeled non-specific probes, such as a 100-fold excess of the -537A competitor, the -537G oligonucleotide remained unaffected (lane 6). However, the -537A oligonucleotide probe could not bind the transcription factor (lane 12), and displayed no band with a 100-fold excess of -537G probe. Therefore, we conclude that the -537G oligonucleotide has more specific DNA-protein binding capacity than the -537A probe. EMSA was performed to investigate the binding between the -537G and A alleles in HeLa nuclear extracts. The -537A allele displayed weak transcription factor binding activity compared to the -537G allele. The -158G/A polymorphism did not display differences in binding activity in EMSA. Our results strongly indicate that the -537G/A polymorphism of GSTT2 specifically influences transcription factor binding activity, leading to a decrease or increase in GSTT2 expression.
GSTT2 exhibits high glutathione peroxidase activity with cumene hydroperoxide as a substrate . In human colon cancer cells, butyrate and flavonoids that contribute to detoxification of dietary carcinogens induce upregulation of GSTT2 to protect against toxic products of oxygen and lipid peroxidation [18, 23]. This gene may have an important role in carcinogenesis and sensitivity of tumors against oxidation stress. GSTT2 promoter polymorphisms significantly reduce luciferase activity in two human cell lines (HEK239t and TE671) . SNPs or mutations within the promoter region reportedly affect transcription activity and gene expression . Therefore, SNPs and haplotypes of the GSTT2 promoter region result in distinct GSTT2 activities, and the -537A allele may be associated with cancer development arising from oxidation stress.
In summary, case-control analysis of CRC patients reveals that GSTT2 polymorphisms are associated with CRC risk. Individual SNP analyses and haplotype results were consistent. EMSA results revealed lower promoter binding activity of the -537A allele, compared to the -537G allele. To our knowledge, this is the first report that GSTT2 promoter polymorphisms and their haplotypes are associated with colorectal cancer risk. Our data collectively suggest that GSTT2 promoter polymorphisms are involved in CRC development, and it would thus be beneficial to include GSTT2 promoter SNPs when screening for relationships between GST families and human cancers.
Shin HR, Jung KW, Won YJ, Park JG: 2002 Annual report of the Korea Central Cancer Registry: Based on registered data from 139 hospitals. Cancer Res Treat. 2004, 36: 103-114.
Kim IJ, Shim Y, Kang HC, Park JH, Ku JL, Park HW, Park HR, Lim SB, Jeong ST, Kim WH, Park JG: Robust microsatellite instability (MSI) analysis by denaturing high-performance liquid chromatography (DHPLC). J Hum Genet. 2003, 48 (10): 525-530. 10.1007/s10038-003-0070-y.
Sheweita SA, Tilmisany AK: Cancer and phase II drug-metabolizing enzymes. Curr Drug Metab. 2003, 4 (1): 45-58. 10.2174/1389200033336919.
Cartwright RA, Glashan RW, Rogers HJ, Ahmad RA, Barham-Hall D, Higgins E, Kahn MA: Role of N-acetyltransferase phenotypes in bladder carcinogenesis: a pharmacogenetic epidemiological approach to bladder cancer. Lancet. 1982, 2 (8303): 842-845. 10.1016/S0140-6736(82)90810-8.
Townsend D, Tew K: Cancer drugs, genetic variation and the glutathione-S-transferase gene family. Am J Pharmacogenomics. 2003, 3 (3): 157-172. 10.2165/00129785-200303030-00002.
Hayes JD, Flanagan JU, Jowsey IR: Glutathione transferases. Annu Rev Pharmacol Toxicol. 2005, 45: 51-88. 10.1146/annurev.pharmtox.45.120403.095857.
Butler WJ, Ryan P, Roberts-Thomson IC: Metabolic genotypes and risk for colorectal cancer. J Gastroenterol Hepatol. 2001, 16 (6): 631-635. 10.1046/j.1440-1746.2001.02501.x.
Peters WH, Wormskamp NG, Thies E: Expression of glutathione S-transferases in normal gastric mucosa and in gastric tumors. Carcinogenesis. 1990, 11 (9): 1593-1596. 10.1093/carcin/11.9.1593.
Nakazato H, Suzuki K, Matsui H, Koike H, Okugi H, Ohtake N, Takei T, Nakata S, Hasumi M, Ito K, Kurokawa K, Yamanaka H: Association of genetic polymorphisms of glutathione-S-transferase genes (GSTM1, GSTT1 and GSTP1) with familial prostate cancer risk in a Japanese population. Anticancer Res. 2003, 23 (3C): 2897-2902.
Lee KM, Park SK, Kim SU, Doll MA, Yoo KY, Ahn SH, Noh DY, Hirvonen A, Hein DW, Kang D: N-acetyltransferase (NAT1, NAT2) and glutathione S-transferase (GSTM1, GSTT1) polymorphisms in breast cancer. Cancer Lett. 2003, 196 (2): 179-186. 10.1016/S0304-3835(03)00311-2.
Balendiran GK, Dabur R, Fraser D: The role of glutathione in cancer. Cell Biochem Funct. 2004, 22 (6): 343-352. 10.1002/cbf.1149.
Tan KL, Webb GC, Baker RT, Board PG: Molecular cloning of a cDNA and chromosomal localization of a human theta-class glutathione S-transferase gene (GSTT2) to chromosome 22. Genomics. 1995, 25 (2): 381-387. 10.1016/0888-7543(95)80037-M.
Deakin M, Elder J, Hendrickse C, Peckham D, Baldwin D, Pantin C, Wild N, Leopard P, Bell DA, Jones P, Duncan H, Brannigan K, Alldersea J, Fryer AA, Strange RC: Glutathione S-transferase GSTT1 genotypes and susceptibility to cancer: studies of interactions with GSTM1 in lung, oral, gastric and colorectal cancers. Carcinogenesis. 1996, 17 (4): 881-884. 10.1093/carcin/17.4.881.
Zhang H, Ahmadi A, Arbman G, Zdolsek J, Carstensen J, Nordenskjold B, Soderkvist P, Sun XF: Glutathione S-transferase T1 and M1 genotypes in normal mucosa, transitional mucosa and colorectal adenocarcinoma. Int J Cancer. 1999, 84 (2): 135-138. 10.1002/(SICI)1097-0215(19990420)84:2<135::AID-IJC7>3.0.CO;2-C.
Hussey AJ, Hayes JD: Characterization of a human class-Theta glutathione S-transferase with activity towards 1-menaphthyl sulphate. Biochem J. 1992, 286 (Pt 3): 929-935.
Tan KL: Purification and characterization of a recombinant human Theta-class glutathione transferase (GSTT2-2). Biochem J. 1996, 315 (Pt 3): 727-732.
Marnett LJ: Oxy radicals, lipid peroxidation and DNA damage. Toxicology. 2002, 181–182: 219-22. 10.1016/S0300-483X(02)00448-1.
Pool-Zobel BL, Selvaraju V, Sauer J, Kautenburger T, Kiefer J, Richter KK, Soom M, Wolfl S: Butyrate may enhance toxicological defence in primary, adenoma and tumor human colon cells by favourably modulating expression of glutathione S-transferases genes, an approach in nutrigenomics. Carcinogenesis. 2005, 26 (6): 1064-1076. 10.1093/carcin/bgi059.
Landi S, Gemignani F, Moreno V, Gioia-Patricola L, Chabrier A, Guino E, Navarro M, de Oca J, Capella G, Canzian F, Bellvitge Colorectal Cancer Study Group: A comprehensive analysis of phase I and phase II metabolism gene polymorphisms and risk of colorectal cancer. Pharmacogenet Genomics. 2005, 15 (8): 535-546.
Guy CA, Hoogendoorn B, Smith SK, Coleman S, O'Donovan MC, Buckland PR: Promoter polymorphisms in glutathione-S-transferase genes affect transcription. Pharmacogenetics. 2004, 14 (1): 45-51. 10.1097/00008571-200401000-00005.
Stephens M, Smith NJ, Donnelly P: A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001, 68 (4): 978-989. 10.1086/319501.
Stephens M, Donnelly P: A comparison of bayesian methods for haplotype reconstruction from population genotype data. Am J Hum Genet. 2003, 73 (5): 1162-1169. 10.1086/379378.
Veeriah S, Kautenburger T, Habermann N, Sauer J, Dietrich H, Will F, Pool-Zobel BL: Apple flavonoids inhibit growth of HT29 human colon cancer cells and modulate expression of genes involved in the biotransformation of xenobiotics. Mol Carcinog. 2006, 45 (3): 164-174. 10.1002/mc.20158.
Coleman SL, Buckland PR, Hoogendoorn B, Guy C, Smith K, O'Donovan MC: Experimental analysis of the annotation of promoters in the public database. Hum Mol Genet. 2002, 11 (16): 1817-1821. 10.1093/hmg/11.16.1817.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/7/16/prepub
This work was supported by a research grant from the National Cancer Center, Korea and the BK21 project for Medicine, Dentistry and Pharmacy.
The author(s) declare that they have no competing interests.
SGJ, IJK, and HCK participated in the design of the study, and drafted and wrote the manuscript. SGJ, IJK, HCK, and HWP participated in the statistical analysis. SAA and HJY participated in the production of genotype data. SGJ and KK carried out the EMSA, and HRS and JSL participated in the acquisition and interpretation of data. JGP participated in the design and coordination of the study. All authors read and approved the final manuscript.
Sang-Geun Jang, Il-Jin Kim, Hio Chung Kang contributed equally to this work.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.