Skip to main content

Genetic variation in Glutathione S-Transferase Omega-1, Arsenic Methyltransferase and Methylene-tetrahydrofolate Reductase, arsenic exposure and bladder cancer: a case–control study



Ingestion of groundwater with high concentrations of inorganic arsenic has been linked to adverse health outcomes, including bladder cancer, however studies have not consistently observed any elevation in risk at lower concentrations. Genetic variability in the metabolism and clearance of arsenic is an important consideration in any investigation of its potential health risks. Therefore, we examined the association between genes thought to play a role in the metabolism of arsenic and bladder cancer.


Single nucleotide polymorphisms (SNPs) in GSTO-1, As3MT and MTHFR were genotyped using DNA from 219 bladder cancer cases and 273 controls participating in a case–control study in Southeastern Michigan and exposed to low to moderate (<50 μg/L) levels of arsenic in their drinking water. A time-weighted measure of arsenic exposure was constructed using measures from household water samples combined with past residential history, geocoded and merged with archived arsenic data predicted from multiple resources.


While no single SNP in As3MT was significantly associated with bladder cancer overall, several SNPs were associated with bladder cancer among those exposed to higher arsenic levels. Individuals with one or more copies of the C allele in rs11191439 (the Met287Thr polymorphism) had an elevated risk of bladder cancer (OR = 1.17; 95% CI = 1.04-1.32 per 1 μg/L increase in average exposure). However, no association was observed between average arsenic exposure and bladder cancer among TT homozygotes in the same SNP. Bladder cancer cases were also 60% less likely to be homozygotes for the A allele in rs1476413 in MTHFR compared to controls (OR = 0.40; 95% CI = 0.18-0.88).


Variation in As3MT and MTHFR is associated with bladder cancer among those exposed to relatively low concentrations of inorganic arsenic. Further investigation is warranted to confirm these findings.

Peer Review reports


Cancer of the urinary bladder (or bladder) is the seventh most common malignancy worldwide among men and the 17th among women with approximately 386,000 new cases diagnosed each year [1, 2]. In the United States (U.S.), bladder cancer is the fourth most common cancer diagnosed among men and the eleventh among women with a projected 73,510 new cases to be diagnosed in 2012 [3]. A vast majority of bladder tumors arise in the urothelium and the most common histology is transitional cell carcinoma (TCC).

Ingestion of inorganic arsenic primarily from groundwater supplies has been linked to a number of adverse health conditions including both benign and malignant skin lesions, peripheral vascular disease and cancers of the kidney, lung and bladder [4, 5]. The data to support these relationships were first gathered from investigations in areas of chronic and endemic arseniasis, however, in the past decade there has been an increasing amount of attention devoted to study of the health consequences of low to moderate concentrations of arsenic consumed from contaminated household drinking water (<100 μg/L). While studies in Argentina, Finland and the U.S. indicate an increased risk of bladder cancer among smokers exposed to higher levels of arsenic, they do not report any increase in risk overall associated with concentrations in this range [69]. In a recent report, our findings from a population-based case–control study conducted in Southeastern Michigan are in agreement with other U.S. studies indicating no significant association overall between arsenic exposure from drinking water and bladder cancer risk [10]. However the limitations inherent in the design of these studies, particularly the difficulty in estimating arsenic exposure over an extended period of time to account for the long-induction latency period for cancer, introduces the possibility of exposure misclassification which may explain the null findings. Results have not been uniformly negative, as studies in Northeastern Taiwan, Chile and others in Argentina have reported increases in bladder cancer risk and mortality associated with arsenic concentrations as low as 10-50 μg/L depending upon the length of exposure [1115].

Crucial to any investigation of an environmental exposure such as arsenic and potential health risks is the consideration of interindividual variability in its metabolism and clearance. The primary pathway for the metabolism of inorganic arsenic is methylation to monomethylated (monomethylarsonic acid [MMAV+, monomethylarsonous acid [MMAIII+) and eventually dimethylated (dimethylarsinic acid [DMAV+, dimethylarsinous acid [DMAIII+) forms [16]. In humans, the process is incomplete and all forms are excreted in urine in varying proportions, with arsenite and monomethylated forms considered the most bioactive and linked to enhanced risk of arsenic-induced health hazards including bladder cancer [1720]. The typical urine profile in populations exposed to inorganic arsenic contains 10-30% inorganic arsenic, 10-20% MMA and 60-80% DMA [16]. The interindividual variation in the relative proportion of inorganic arsenic and its methylated metabolites excreted in urine has been attributed to differences in gender, age, pregnancy, smoking behavior and diet [21]. Moreover, genetics have long been suspected to play an important role in arsenic metabolism [21, 22]. Polymorphisms in a number of genes have been associated with variation in methylation capacity with Glutathione S-Transferase Omega-1 (GSTO-1), Methylene-tetrahydrofolate Reductase (MTHFR) and Arsenic (+3) Methyltransferase (As3MT) as three of the more widely investigated genes [2333].

Therefore, we set out to examine the association between specific polymorphisms in each of these genes with bladder cancer in a pilot case–control study. This study is one of the first to investigate the interaction between genes involved in arsenic metabolism and lifetime arsenic exposure on risk of bladder cancer conducted in a residentially stable population in Southeast Michigan.


A population-based case–control study was conducted in eleven counties in Southeastern Michigan based on evidence of moderately elevated (<100 μg/L) arsenic concentrations in groundwater supplies. It has been estimated that approximately 8% of the total population in the counties from which our participants were selected are exposed to arsenic concentrations which exceed the Maximum Contaminant Level (MCL) of 10 μg/L set by the U.S. Environmental Protection Agency [34]. A detailed summary of study methods has been previously reported [10]. Eligible participants were residents of Genesee, Huron, Ingham, Jackson, Lapeer, Livingston, Oakland, Sanilac, Shiawasee, Tuscola and Washtenaw counties for at least 5 years prior to diagnosis (bladder cancer cases) and interview (controls). Bladder cancer cases diagnosed between January 1st, 2000 and December 31st, 2004 were identified from the Michigan Cancer Surveillance Program, the state’s cancer registry, within the Division of Vital Records, Michigan Department of Community Health. Cases were aged 21 to 80 years at date of diagnosis, histopathologically confirmed, with no history of invasive cancer prior to bladder cancer diagnosis. Of the 1,634 potentially eligible cases, 25% (n = 411) completed all phases of the initial study protocol. The most common reasons for non-participation were 1.) refusing permission to registry for initial release of contact information (41%); and 2.) death prior to registry contact (22%). Disease-free controls were identified through random-digit dialing and frequency matched to cases by age (± 5 years), gender and race. Among the 2,132 eligible controls identified, 27% (n = 566) completed all phases of the initial protocol.

All participants completed a computer-assisted telephone interview conducted by trained interviewers from the Michigan Public Health Institute (MPHI). The interview collected information on demographics, water and other fluid consumption both at home and away from home, a dietary questionnaire, smoking, medical and occupational history. Patterns of fluid consumption over the participant’s adult lifetime were reported, including a detailed characterization of any major changes in fluid consumption. After the initial interview, research team members scheduled a visit to the home of each participant to collect complete residential and occupational histories. Participants were asked to provide information on duration of residence/employment, drinking water source at each location, change in water source and a street address (or nearest cross streets). Each residence in the study area was geocoded and assigned a specific coordinate in ArcGIS (v.9.0; ESRI, Redlands, CA) with no geocoding of residences outside the study area. During the in-home interview, household water and biospecimens were also collected.

A genetic component to the larger case–control study was initiated in 2005, at the approximate midpoint in subject recruitment. The protocol was amended to request that participants provide a saliva specimen during their scheduled in-home visit. A total of 219 cases and 273 controls, representing ~98% of subjects recruited during this time, submitted saliva samples for genetic studies. Genomic DNA was isolated from these specimens using the Puregene DNA Purification buccal cell kit (Gentra Systems, Minneapolis, MN). All protocols were reviewed and approved by the University of Michigan Health Sciences Institutional Review Board and all participants provided written informed consent.

Measure of lifetime arsenic exposure

Our approach for estimating average lifetime arsenic exposure and analysis of household water samples has been previously described [10, 35]. Briefly, an estimation of arsenic concentrations from past residences on public well water supplies relied on data that were extracted from the Michigan state arsenic database and via the contact of specific municipality within the study area where the participant reported a prior residency. Public surface water supplies from past residences were assigned the average measured level of surface water in the study population (0.30 μg/L). For past private wells inside the study area, a geostatistical model was developed for predicting arsenic concentrations [36]. For public supplies or private wells outside of the study area, arsenic concentrations were estimated using city averages from the U.S. Geological Survey (USGS) database [37]. Estimates of average lifetime residential exposure to arsenic were constructed independent of case–control status, combining the measured concentrations from household water samples collected from the subject’s current residence and data gathered from past residences with the assumption that arsenic concentrations were relatively stable within a particular well. This assumption is supported by our own internal investigation [38] and in the literature [39]. The measure was then weighted by duration of residence with an average (in μg/L) calculated across the adult lifetime [10]. The household water samples were analyzed for arsenic concentration using an inductively coupled plasma mass spectrometry unit (ICP-MS, Model 7500c, Agilent Technologies, Santa Clara, CA) at the University of Michigan.

SNP selection and genotyping methods

A total of twelve single nucleotide polymorphisms (SNPs) in GSTO-1 (rs4925, rs2282326), As3MT (rs1046778, rs11191438, rs11191439, rs3740400, rs7085104) and MTHFR (rs1801131, rs1801133, rs6541003, rs9651118, rs4846048, rs1476413) were genotyped using Taqman SNP genotyping assays (Applied Biosystems, Foster City, CA). These SNPs were selected based upon previous citation/functional relevance in the literature, or a minor allele frequency of ≥ 10% among CEPH (Centre d’Etude du Polymorphisme Humain) Utah European ancestry samples (CEU) in HapMap and an ability to tag major haplotype blocks within each gene in the same samples [40]. We used the ABI PRISM 7900HT Sequence Detection System and the SDS software (Applied Biosystems, version 2.1) to distinguish SNP alleles. PCR primers were purchased from Invitrogen Life Technologies (Carlsbad, California) and are available upon request. Each PCR reaction contained 25μL total volume (2.5μL of 10X PCR Buffer (Invitrogen Life Technologies, Inc.), 0.5 μL 50 mM MgCl2, 0.5μL 10 mM dNTPs, 2.5 μL each of the two PCR primers at 5 μM concentration, 2 μL of template DNA at 20 ng/μL, 0.125 μL of PlatinumTaq Polymerase (Invitrogen), and 14.375 μL of H20. PCR products were cleaned using MinElute PCR purification kit (Germantown, Massachusetts) and sequenced using an ABI prism 3100 Genetic Analyzer using Big Dye Terminator v1.1 chemistries (Applied Biosystems, Foster City, California). The average genotyping call rate was 97.5%, with call rates ≥ 94.6% for each SNP. Undetermined samples from an assay were sequenced if DNA was available, resulting in a final genotyping call rate of 99.5%, with call rates ≥ 98.8% for each SNP. A random subset of samples (2%) were also duplicated and verified by direct sequencing without discrepancy.

Statistical methods

Statistical analyses were performed using Statistical Analysis Software (SAS Institute Inc. v. 9.1, Cary, N.C.). All tests were evaluated using a two-sided hypothesis test with statistical significance interpreted as having a p-value less than or equal to 0.05. Differences in age, gender, race, smoking history, average arsenic exposure and the distribution of genotypes for each SNP between bladder cancer cases and controls were tested using either the Wilcoxon rank-sum test or a simple chi-square test. For each SNP, we calculated the minor allele frequency among controls (Table 1).

Table 1 Single nucleotide polymorphisms in Glutathione S-Transferase Omega, Arsenic Methyltransferase and 5,10 Methylene-tetrahydrofolate Reductase

The observed genotype distribution for each SNP among controls was tested for consistency with HWE expectations using Pearson’s chi-square test. Linkage disequilibrium (LD) between selected SNPs within each gene was evaluated using the PLINK software package ( )[41] Unconditional logistic regression was used to determine whether genotypes within each SNP were associated with bladder cancer under the log-additive, dominant and recessive models. A post-hoc test comparing homozygotes for the minor allele and heterozygotes each to homozygotes for the major allele (the referent) was also performed. Multivariable models simultaneously adjusted for the effect of each SNP, age (continuous), and smoking history (ever/never). Because adjustment for race and gender did not alter the odds ratios, they were not included in the final models. The joint effect of each SNP and arsenic concentration was evaluated by including an interaction term in the regression models which combined genotype (dominant model) and average arsenic concentration measured continuously. Additional analyses stratified by genotype were performed, calculating odds ratios for bladder cancer associated with arsenic measured both continuously and according to quartiles in the distribution of exposure among controls. The p-values presented in the tables correspond to the β coefficient for the interaction term. Potential haplotypes for As3MT were constructed from genotype data and differences in haplotype frequency (excluding rare haplotypes [<2%]) between cases and controls were evaluated using the Haplo.stat software package The same package was used to estimate the association between haplotype (under a log-additive model) and bladder cancer, adjusting for age and smoking history [42].


Relevant characteristics of the 219 bladder cancer cases and 273 controls participating in the pilot investigation are summarized in Table 2. By design, cases and controls were similar with respect to age, race and gender. Approximately 75% of participants were men, 95% were white and the mean age at diagnosis (among cases) and at interview (among controls) was 66 years (standard deviation ± 9.5 years). Bladder cancer cases were more likely to report that they had either formerly smoked or were current smokers at the time of diagnosis compared to controls at time of interview (78% versus 53%; p < 0.0001). Furthermore, the estimated mean number of pack-years among smokers was higher among cases than controls. As has been reported in the larger case–control study, there was no significant difference between genotyped cases and controls in lifetime average exposure to arsenic (p = 0.45) with an estimated median concentration of 1.36 μg/L (interquartile range = 3.33 μg/L) in the study population.

Table 2 Characteristics of Study Participants

The genotype distribution observed for each of the tested SNPs among controls was consistent with HWE. Table 3 summarizes the findings of our post-hoc testing of associations between genotype for each SNP and bladder cancer risk adjusted for age and smoking history, however this adjustment produced little difference from the crude estimates of risk. We did not observe any significant associations between SNPs in GSTO-1, As3MT and MTHFR and bladder cancer risk overall with one exception. Bladder cancer cases were 60% less likely to be homozygotes for the A allele in rs1476413 in MTHFR compared to controls (OR = 0.40; 95% CI = 0.18-0.88).

Table 3 Odds of bladder cancer by genotype in genes of interest (Adjusted Odds Ratios and 95% Confidence Intervals)

We then stratified the study population according to genotype, examining the differences in the distribution between time-weighted average (TWA) of arsenic exposure (divided into quartiles and measured continuously) between cases and controls (Tables 4 and 5). There was evidence for interaction between several of the SNPs in As3MT and arsenic exposure. Among persons with one or more copies of the C allele in rs11191439 (Met287Thr), higher average arsenic exposure was associated with an increase in the odds of bladder cancer (OR = 1.17; 95% CI = 1.04-1.32 per 1 μg/L increase). However, no association was observed between average arsenic exposure and bladder cancer among individuals who were TT homozygotes for the same SNP. This same polymorphism was also associated with an increase in the odds of bladder cancer among all study subjects, however the estimate was not statistically significant (OR = 2.75; 95%CI = 0.47-16.77 among CC homozygotes). Interactions with SNPs rs3740400, rs11191438 and rs7085104 and arsenic exposure were also observed (Table 4). These SNPs were in strong LD with one another (r2 values ranging from 0.83-0.95) ( Additional file 1: Table S1). And while there was nominal evidence for interaction with the rs4925 SNP in GSTO-1 and arsenic exposure, we not find evidence for interaction with selected SNPs in MTHFR and arsenic (Table 5). Haplotype-based association analyses for the As3MT gene were consistent with the single SNP results (Table 6). No evidence for an association between haplotypes in the As3MT gene and bladder cancer was detected when using the entire sample of participants. There was, however, an association between the As3MT haplotype which included the Met287Thr polymorphism (compared to the most common or referent haplotype) and bladder cancer for participants exposed to arsenic concentrations higher than 3.72 μg/L (the uppermost quartile) (OR = 3.90; 95% CI = 1.39-10.9), but not among those exposed to lower concentrations.

Table 4 Odds of bladder cancer associated with As3MT SNPs and arsenic from drinking water (in quartiles†)
Table 5 Odds of bladder cancer associated with GSTO-1 and MTHFR SNPs and arsenic from drinking water†
Table 6 Common haplotypes in the As3MT gene and bladder cancer by average level of arsenic from drinking water


In 2006, the U.S. Environmental Protection Agency lowered the MCL for arsenic concentrations in drinking water provided by public water supplies in the U.S. from 50 μg/L to 10 μg/L, in response to growing concern and mounting evidence that inorganic arsenic was hazardous to human health at relatively low levels [4, 43]. Results from studies conducted across several regions in the U.S. generally do not indicate these levels elevate risk of bladder cancer, however some studies report elevated risks among certain subgroups in the population (i.e. smokers) in conjunction with low-level arsenic exposure[6, 810]. Findings from the current investigation suggest that possessing one or more copies of the Met287Thr polymorphism in the As3MT gene (rs11191439), coupled with higher average arsenic exposure is associated with an increase in bladder cancer risk. Our study population had levels of arsenic primarily below the current MCL, therefore the findings pertain to potential genetic susceptibility to effects of arsenic at relatively low-levels.

While our study is the first to report an interaction between As3MT (10q24.32) and arsenic exposure on bladder cancer risk, these findings are consistent with evidence gleaned from several other studies demonstrating that polymorphisms in As3MT are involved in the methylation of inorganic arsenic [2429, 33, 44]. A study of urinary arsenic species in a Mexican population exposed to arsenic at concentrations averaging 5 μg/L to 43 μg/L was the first to report that three polymorphic sites in As3MT were associated with a higher ratio of DMAV to MMAV, with a single polymorphism (rs11191453) also associated with a higher ratio of arsenite to MMAV in children, but not adults [33]. Subsequent studies have pointed to the Met287Thr polymorphism as being associated with enhanced As3MT enzymatic activity [45] and a higher percentage of MMA to DMA in urine samples among individuals carrying at least one copy of the variant (C) allele [24, 29, 44, 46]. Furthermore, the findings of Lindberg et al., suggest that the influence of the Met387Thr polymorphism is greater among men than women explaining approximately 20% of the variation in methylation ability in men, but just 4% in women [47]. In a recent report, no association was observed between any of eight tested SNPs in As3MT and bladder cancer overall and no clear evidence of effect modification by A3sMT genotype. However, despite the study’s limited statistical power to address interaction between arsenic exposure and genotype, a number of the odds ratios reported among participants exposed to higher arsenic concentrations were elevated, but with wide confidence intervals [48].

Gomez-Rubio and colleagues caution against attributing enhanced arsenic methylation to any single SNP in As3MT, or for that matter any one gene, as they identified a large (347 kb) region in strong LD including four genes in addition to As3MT [49]. We observed strong pairwise LD between three of four As3MT polymorphisms interacting with arsenic exposure ( Additional file 1: Table S1). However, the Met287Thr polymorphism did not strongly co-segregate with the other tested SNPs. Our results suggest that homozygotes for the more common allele in each of these SNPs (rs3740400, rs11191438 and rs7085104) had a lower risk of bladder cancer in conjunction with arsenic exposure, however minor allele carriers were at an increased risk. Prior investigations of rs3740400, or those observed in LD with this particular SNP, suggest an altered DMA to MMA ratio associated with the minor allele [28, 29, 50]. However, the observed racial and ethnic variation in allele frequency and patterns of LD, coupled with the worldwide geographic variability in not only exposure to arsenic but other risk factors, make the comparison between previous findings and our own difficult and an area deserving further investigation.

There was some evidence, albeit weak, to support an interaction between the rs4925 SNP in GSTO-1 and arsenic exposure in this investigation. Specifically, there was evidence of a positive association between arsenic exposure and bladder cancer among those homozygous wildtype for the GSTO-1 polymorphism that was not present for those with at least one variant allele. GSTO-1, also on chromosome 10 (10q25.1), has demonstrated the ability to catalyze the reduction of MMAV to MMAIII, considered the rate-limiting step in arsenic biotransformation [51]. A strong pattern of LD has been demonstrated in European Americans among the more common (MAF > 10%) polymorphisms identified [32]. The Ala140Asp (rs4925) polymorphism in GSTO-1 has been the most widely-studied, linked to risk of pre-cancerous skin lesions and cancer in arsenic-exposed populations [26, 52, 53], however studies have not confirmed any association between specific variants in GSTO-1 with the relative proportion of arsenic and its metabolites [47, 52, 54].

A greater proportion of controls compared to cases were homozygotes for a single polymorphism in MTHFR (rs1476413). The function, if any, of this particular SNP is unknown, however other polymorphisms in the MTHFR gene have been linked to a number of cancers [55]. Findings with respect to bladder cancer have been mixed [23, 31, 5659]. Methylene-tetrahydrofolate reductase is an enzyme which plays a notable role in the metabolism of folate, and also considered important in arsenic methylation [60]. The 677 C → T (Ala222Val) and 1298 A → C (Glu429Ala) polymorphisms in MTHFR have been investigated relative to their influence on arsenic methylation capacity suggesting that homozygotes for the variant T allele in MTHFR 677 excrete a smaller proportion of DMA to inorganic arsenic in their urine and the opposite effect among homozygotes for the variant C allele in MTHFR 1298 compared to homozygotes for the wild-type for each SNP [23, 47]. We did not find an association with either of these polymorphisms in our investigation.

Consideration of potential sources of bias is important in interpreting the findings from this investigation and a description of potential bias in the larger case–control study has been previously published [10]. The participation rate was relatively low (~25%) introducing the possibility of selection bias. To address this, we compared known characteristics of participating cases and potentially eligible, non-participating cases and found no difference between the two groups with respect to race, gender and urban versus rural residency. However, because participating cases were more likely to be diagnosed with less invasive disease and at a younger age compared to non-participating cases, our findings might not be generalizable to older patients diagnosed with advanced bladder cancer. As nearly one-quarter of potentially eligible cases died prior to registry contact, survival bias is a consideration. Deceased cases were more likely than participating cases to reside in rural areas versus urban areas at time of diagnosis (10% versus 6%). As rural areas tend to have generally higher and more variable levels of arsenic exposure, it is possible that arsenic exposure may be associated with more aggressive disease and earlier death producing a bias for the overall association between arsenic exposure and bladder cancer toward the null. However, published findings to support this notion are inconsistent [6163]. The reasons for non-participation among eligible controls were unavailable; however an evaluation of participating controls suggested that they were geographically representative of the general population of the study area and by design similar in distribution to the cases with respect to age, gender and race. A more detailed description of this evaluation and its findings has been previously reported [10].

Lastly, because the genetic component did not begin until the midpoint in recruitment, genetic material was collected on just 492 (50.4%) of all study participants. A comparison of characteristics between participants with and without DNA found no difference between these groups by gender or smoking history (data not shown). A higher proportion of participants with DNA were white (95%) compared to those without DNA (89%) and exposed to higher arsenic concentrations (median concentration among those with DNA 1.39 μg/L compared to 1.00 μg/L among those without). However, the differences were similar between cases and controls and were a reflection of the overall composition of the study participants recruited in the latter half of the overall investigation. Nearly all (~98%) subjects recruited submitted a saliva specimen during their scheduled in-home visit.

Misclassification of lifetime arsenic exposure is a potential source of bias and may have been introduced from either inaccurate recall of residential history by participants or lack of information on groundwater arsenic concentrations for past residences on private wells. However, one of the strengths of the study was the residential stability of our population. On average, 80% of the person-years used to construct the measure of average arsenic exposure were spent in southeastern Michigan, with just 5% of reported residential histories among subjects on private wells outside of the study area and inadequate for geocoding. We do not anticipate that these sources of misclassification would differ systematically by genotype or disease status and would therefore tend to bias our results toward the null. Lastly, the study was of limited power to detect modest differences between genotype and bladder cancer risk particularly when the “at-risk” allele frequency was less common. And in the case of Met287Thr, imprecision in the estimate of risk associated with higher arsenic exposure due to the small number of exposed cases and controls. The number of participants with arsenic levels greater than the current MCL was prohibitively small to examine gene-environment interaction at arsenic concentrations greater than 10 μg/L.


Results from this pilot investigation suggest that variation in As3MT and MTHFR is associated with bladder cancer among those exposed to relatively low arsenic concentrations. The Met287Thr polymorphism in As3MT, in particular, was associated with an increased risk of bladder cancer among participants exposed to higher levels of arsenic exposure. Interactions between several other SNPs in As3MT and arsenic exposure were also observed; however LD in the region makes it difficult to attribute findings to any single SNP. The underlying mechanism may also involve other genes in close proximity to As3MT. Future studies in larger populations are clearly warranted to confirm these findings as well as to broaden the exploration of genes to their potential association with arsenic in bladder cancer. If validated, these findings highlight the importance of gene-environment interaction in the evaluation of arsenic-mediated health risks in populations exposed to even low levels of arsenic.





Arsenic (3+) Methyltransferase


Aspartic Acid


Dimethylated arsenic


Deoxyribonucleic acid


Glutamic Acid


Glutathione S-Transferase Omega 1


Hardy-Weinberg Equilibrium


Linkage Disequilibrium


Minor allele frequency




Magnesium Choride


Michigan Public Health Institute


Methylene-tetrahydrofolate Reductase


Maximum Contaminant Level


Monomethylated arsenic


Odds Ratio


Polymerase chain reaction


Single Nucleotide Polymorphism


Transitional Cell Carcinoma




Time weighted average


United States






Micrograms per liter.


  1. Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D: Global cancer statistics. CA Cancer J Clin. 2011, 61: 69-90. 10.3322/caac.20107.

    Article  Google Scholar 

  2. Parkin DM, Bray F, Ferlay J, Pisani P: Global cancer statistics, 2002. CA Cancer J Clin. 2005, 55: 74-108. 10.3322/canjclin.55.2.74.

    Article  Google Scholar 

  3. Siegel R, Naishadham D, Jemal A: Cancer statistics, 2012. CA Cancer J Clin. 2012, 62: 10-29. 10.3322/caac.20138.

    Article  Google Scholar 

  4. National Research Council: Arsenic in Drinking Water. 2001, National Academy of Sciences Press, Washington D.C

    Google Scholar 

  5. IARC: monographs on the evaluation of carcinogenic risks to humans. [84]. 2004, International Agency for Research on Cancer, Lyon, France

    Google Scholar 

  6. Bates MN, Smith AH, Cantor KP: Case–control study of bladder cancer and arsenic in drinking water. Am J Epidemiol. 1995, 141: 523-530.

    CAS  Google Scholar 

  7. Kurttio P, Pukkala E, Kahelin H, Auvinen A, Pekkanen J: Arsenic concentrations in well water and risk of bladder and kidney cancer in Finland. Environ Health Perspect. 1999, 107: 705-710. 10.1289/ehp.99107705.

    Article  CAS  Google Scholar 

  8. Karagas MR, Tosteson TD, Morris JS, Demidenko E, Mott LA, Heaney J, Schned A: Incidence of transitional cell carcinoma of the bladder and arsenic exposure in New Hampshire. Canc Causes Contr. 2004, 15: 465-472.

    Article  Google Scholar 

  9. Steinmaus C, Yuan Y, Bates MN, Smith AH: Case–control study of bladder cancer and drinking water arsenic in the western United States. Am J Epidemiol. 2003, 158: 1193-1201. 10.1093/aje/kwg281.

    Article  Google Scholar 

  10. Meliker JR, Slotnick MJ, AvRuskin GA, Schottenfeld D, Jacquez GM, Wilson ML, Goovaerts P, Franzblau A, Nriagu JO: Lifetime exposure to arsenic in drinking water and bladder cancer: a population-based case–control study in Michigan, USA. Canc Causes Contr. 2010, 21: 745-757. 10.1007/s10552-010-9503-z.

    Article  Google Scholar 

  11. Hopenhayn-Rich C, Biggs ML, Fuchs A, Bergoglio R, Tello EE, Nicolli H, Smith AH: Bladder cancer mortality associated with arsenic in drinking water in Argentina. Epidemiology. 1996, 7: 117-124. 10.1097/00001648-199603000-00003.

    Article  CAS  Google Scholar 

  12. Smith AH, Goycolea M, Haque R, Biggs ML: Marked increase in bladder and lung cancer mortality in a region of Northern Chile due to arsenic in drinking water. Am J Epidemiol. 1998, 147: 660-669. 10.1093/oxfordjournals.aje.a009507.

    Article  CAS  Google Scholar 

  13. Chiou HY, Chiou ST, Hsu YH, Chou YL, Tseng CH, Wei ML, Chen CJ: Incidence of transitional cell carcinoma and arsenic in drinking water: a follow-up study of 8,102 residents in an arseniasis-endemic area in northeastern Taiwan. Am J Epidemiol. 2001, 153: 411-418. 10.1093/aje/153.5.411.

    Article  CAS  Google Scholar 

  14. Marshall G, Ferreccio C, Yuan Y, Bates MN, Steinmaus C, Selvin S, Liaw J, Smith A: H: Fifty-year study of lung and bladder cancer mortality in Chile related to arsenic in drinking water. J Natl Canc Inst. 2007, 99: 920-928. 10.1093/jnci/djm004.

    Article  CAS  Google Scholar 

  15. Chen CL, Chiou HY, Hsu LI, Hsueh YM, Wu MM, Wang YH, Chen CJ: Arsenic in drinking water and risk of urinary tract cancer: a follow-up study from northeastern Taiwan. Canc Epidemiol Biomarkers Prev. 2010, 19: 101-110. 10.1158/1055-9965.EPI-09-0333.

    Article  CAS  Google Scholar 

  16. Vahter M: Methylation of inorganic arsenic in different mammalian species and population groups. Sci Prog. 1999, 82 (Pt 1): 69-88.

    CAS  Google Scholar 

  17. Chen YC, Su HJ, Guo YL, Hsueh YM, Smith TJ, Ryan LM, Lee MS, Christiani DC: Arsenic methylation and bladder cancer risk in Taiwan. Canc Causes Contr. 2003, 14: 303-310. 10.1023/A:1023905900171.

    Article  Google Scholar 

  18. Steinmaus C, Bates MN, Yuan Y, Kalman D, Atallah R, Rey OA, Biggs ML, Hopenhayn C, Moore LE, Hoang BK, Smith AH: Arsenic methylation and bladder cancer risk in case–control studies in Argentina and the United States. J Occup Environ Med. 2006, 48: 478-488. 10.1097/01.jom.0000200982.28276.70.

    Article  CAS  Google Scholar 

  19. Huang YK, Huang YL, Hsueh YM, Yang MH, Wu MM, Chen SY, Hsu LI, Chen CJ: Arsenic exposure, urinary arsenic speciation, and the incidence of urothelial carcinoma: a twelve-year follow-up study. Canc Causes Contr. 2008, 19: 829-839. 10.1007/s10552-008-9146-5.

    Article  Google Scholar 

  20. Hernandez A, Marcos R: Genetic variations associated with interindividual sensitivity in the response to arsenic exposure. Pharmacogenomics. 2008, 9: 1113-1132. 10.2217/14622416.9.8.1113.

    Article  CAS  Google Scholar 

  21. Tseng CH: A review on environmental factors regulating arsenic methylation in humans. Toxicol Appl Pharmacol. 2009, 235: 338-350. 10.1016/j.taap.2008.12.016.

    Article  CAS  Google Scholar 

  22. Chung JS, Kalman DA, Moore LE, Kosnett MJ, Arroyo AP, Beeris M, Mazumder DN, Hernandez AL, Smith AH: Family correlations of arsenic methylation patterns in children and parents exposed to high concentrations of arsenic in drinking water. Environ Health Perspect. 2002, 110: 729-733.

    Article  CAS  Google Scholar 

  23. Steinmaus C, Moore LE, Shipp M, Kalman D, Rey OA, Biggs ML, Hopenhayn C, Bates MN, Zheng S, Wiencke JK, Smith AH: Genetic polymorphisms in MTHFR 677 and 1298, GSTM1 and T1, and metabolism of arsenic. J Toxicol Environ Health. 2007, 70: 159-170. 10.1080/15287390600755240.

    Article  CAS  Google Scholar 

  24. Hernandez A, Xamena N, Surralles J, Sekaran C, Tokunaga H, Quinteros D, Creus A, Marcos R: Role of the Met(287)Thr polymorphism in the AS3MT gene on the metabolic arsenic profile. Mutat Res. 2008, 637: 80-92. 10.1016/j.mrfmmm.2007.07.004.

    Article  CAS  Google Scholar 

  25. Marnell LL, Garcia-Vargas GG, Chowdhury UK, Zakharyan RA, Walsh B, Avram MD, Kopplin MJ, Cebrian ME, Silbergeld EK, Aposhian HV: Polymorphisms in the human monomethylarsonic acid (MMA V) reductase/hGSTO1 gene and changes in urinary arsenic profiles. Chem Res Toxicol. 2003, 16: 1507-1513. 10.1021/tx034149a.

    Article  CAS  Google Scholar 

  26. Chung CJ, Hsueh YM, Bai CH, Huang YK, Huang YL, Yang MH, Chen CJ: Polymorphisms in arsenic metabolism genes, urinary arsenic methylation profile and cancer. Canc Causes Contr. 2009, 20: 1653-1661. 10.1007/s10552-009-9413-0.

    Article  Google Scholar 

  27. Schlawicke EK, Nermell B, Concha G, Stromberg U, Vahter M, Broberg K: Arsenic metabolism is influenced by polymorphisms in genes involved in one-carbon metabolism and reduction reactions. Mutat Res. 2009, 667: 4-14. 10.1016/j.mrfmmm.2008.07.003.

    Article  Google Scholar 

  28. Schlawicke EK, Broberg K, Concha G, Nermell B, Warholm M, Vahter M: Genetic polymorphisms influencing arsenic metabolism: evidence from Argentina. Environ Health Perspect. 2007, 115: 599-605. 10.1289/ehp.9734.

    Article  Google Scholar 

  29. Agusa T, Fujihara J, Takeshita H, Iwata H: Individual Variations in Inorganic Arsenic Metabolism Associated with AS3MT Genetic Polymorphisms. Int J Mol Sci. 2011, 12: 2351-2382. 10.3390/ijms12042351.

    Article  CAS  Google Scholar 

  30. Agusa T, Iwata H, Fujihara J, Kunito T, Takeshita H, Minh TB, Trang PT, Viet PH, Tanabe S: Genetic polymorphisms in glutathione S-transferase (GST) superfamily and arsenic metabolism in residents of the Red River Delta, Vietnam. Toxicol Appl Pharmacol. 2010, 242: 352-362. 10.1016/j.taap.2009.11.007.

    Article  CAS  Google Scholar 

  31. Karagas MR, Park S, Nelson HH, Andrew AS, Mott L, Schned A, Kelsey KT: Methylenetetrahydrofolate reductase (MTHFR) variants and bladder cancer: a population-based case–control study. Int J Hyg Environ Health. 2005, 208: 321-327. 10.1016/j.ijheh.2005.04.005.

    Article  CAS  Google Scholar 

  32. Yu L, Kalla K, Guthrie E, Vidrine A, Klimecki WT: Genetic variation in genes associated with arsenic metabolism: glutathione S-transferase omega 1–1 and purine nucleoside phosphorylase polymorphisms in European and indigenous Americans. Environ Health Perspect. 2003, 111: 1421-1427. 10.1289/ehp.6420.

    Article  CAS  Google Scholar 

  33. Meza MM, Yu L, Rodriguez YY, Guild M, Thompson D, Gandolfi AJ, Klimecki WT: Developmentally restricted genetic determinants of human arsenic metabolism: association between urinary methylated arsenic and CYT19 polymorphisms in children. Environ Health Perspect. 2005, 113: 775-781. 10.1289/ehp.7780.

    Article  CAS  Google Scholar 

  34. Frost FJ, Muller T, Petersen HV, Thomson B, Tollestrup K: Identifying US populations for the study of health effects related to drinking water arsenic. J Expo Anal Environ Epidemiol. 2003, 13: 231-239. 10.1038/sj.jea.7500275.

    Article  CAS  Google Scholar 

  35. Meliker JR, Slotnick MJ, AvRuskin GA, Kaufmann A, Fedewa SA, Goovaerts P, Jacquez GJ, Nriagu JO: Individual lifetime exposure to inorganic arsenic using a space-time information system. Int Arch Occup Environ Health. 2007, 80: 184-197. 10.1007/s00420-006-0119-2.

    Article  CAS  Google Scholar 

  36. Meliker JR, AvRuskin GA, Slotnick MJ, Goovaerts P, Schottenfeld D, Jacquez GM, Nriagu JO: Validity of spatial models of arsenic concentrations in private well water. Environ Res. 2008, 106: 42-50. 10.1016/j.envres.2007.09.001.

    Article  CAS  Google Scholar 

  37. Arsenic in Groundwater of the United States: USGS. 2000, Accessed March 30, 2000. [http//]

    Google Scholar 

  38. Slotnick MJ, Meliker JR, Nriagu JO: Effects of time and point-of-use devices on arsenic levels in Southeastern Michigan drinking water, USA. Sci Total Environ. 2006, 369: 42-50. 10.1016/j.scitotenv.2006.04.021.

    Article  CAS  Google Scholar 

  39. Steinmaus CM, Yuan Y, Smith AH: The temporal stability of arsenic concentrations in well water in western Nevada. Environ Res. 2005, 99: 164-168. 10.1016/j.envres.2004.10.003.

    Article  CAS  Google Scholar 

  40. Gibbs RA, Belmont JW, Hardenbol P, Willis TD, Yu FL, Yang HM, Chang LY, Huang W, Liu B, Shen Y, Tam PKH, Tsui LC, Waye MMY, Wong JTF, Zeng CQ, Zhang QR, Chee MS, Galver LM, Kruglyak S, Murray SS, Oliphant AR, Montpetit A, Hudson TJ, Chagnon F, Ferretti V, Leboeuf M, Phillips MS, Verner A, Kwok PY, Duan SH: The International HapMap Project. Nature. 2003, 426: 789-796. 10.1038/nature02168.

    Article  CAS  Google Scholar 

  41. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81: 559-575. 10.1086/519795.

    Article  CAS  Google Scholar 

  42. Schaid DJ, Rowland CM, Tines DE, Jacobson RM, Poland GA: Score tests for association between traits and haplotypes when linkage phase is ambiguous. Am J Hum Genet. 2002, 70: 425-434. 10.1086/338688.

    Article  Google Scholar 

  43. Lubin JH, Beane-Freeman LE, Cantor KP: Inorganic arsenic in drinking water: an evolving public health concern. J Natl Canc Inst. 2007, 99: 906-907. 10.1093/jnci/djm012.

    Article  CAS  Google Scholar 

  44. Hernandez A, Xamena N, Sekaran C, Tokunaga H, Sampayo-Reyes A, Quinteros D, Creus A, Marcos R: High arsenic metabolic efficiency in AS3MT287Thr allele carriers. Pharmacogenet Genomics. 2008, 18: 349-355. 10.1097/FPC.0b013e3282f7f46b.

    Article  CAS  Google Scholar 

  45. Wood TC, Salavagionne OE, Mukherjee B, Wang L, Klumpp AF, Thomae BA, Eckloff BW, Schaid DJ, Wieben ED, Weinshilboum RM: Human arsenic methyltransferase (AS3MT) pharmacogenetics: gene resequencing and functional genomics studies. J Biol Chem. 2006, 281: 7364-7373. 10.1074/jbc.M512227200.

    Article  CAS  Google Scholar 

  46. Valenzuela OL, Drobna Z, Hernandez-Castellanos E, Sanchez-Pena LC, Garcia-Vargas GG, Borja-Aburto VH, Styblo M, Del Razo LM: Association of 6. Toxicol Appl Pharmacol. 2009, 239: 200-207. 10.1016/j.taap.2009.06.007.

    Article  CAS  Google Scholar 

  47. Lindberg AL, Kumar R, Goessler W, Thirumaran R, Gurzau E, Koppova K, Rudnai P, Leonardi G, Fletcher T, Vahter M: Metabolism of low-dose inorganic arsenic in a central European population: influence of sex and genetic polymorphisms. Environ Health Perspect. 2007, 115: 1081-1086. 10.1289/ehp.10026.

    Article  CAS  Google Scholar 

  48. Lesseur C, Gilbert-Diamond D, Andrew AS, Ekstrom M, Li Z, Kelsey KT, Marsit CJ, Karagas MR: A case–control study of polymorphisms in xenobiotic and arsenic metabolism genes and arsenic-related bladder cancer in New Hampshire. Toxicol Lett. 2012, 210: 100-106. 10.1016/j.toxlet.2012.01.015.

    Article  CAS  Google Scholar 

  49. Gomez-Rubio P, Meza-Montenegro MM, Cantu-Soto E, Klimecki WT: Genetic association between intronic variants in AS3MT and arsenic methylation efficiency is focused on a large linkage disequilibrium cluster in chromosome 10. J Appl Toxicol. 2010, 30: 260-270.

    CAS  Google Scholar 

  50. Engstrom K, Vahter M, Mlakar SJ, Concha G, Nermell B, Raqib R, Cardozo A, Broberg K: Polymorphisms in arsenic (+III oxidation state) methyltransferase (AS3MT) predict gene expression of AS3MT as well as arsenic metabolism. Environ Health Perspect. 2011, 119: 182-188. 10.1289/ehp.119-a182.

    Article  Google Scholar 

  51. Zakharyan RA, Sampayo-Reyes A, Healy SM, Tsaprailis G, Board PG, Liebler DC, Aposhian HV: Human monomethylarsonic acid (MMA(V)) reductase is a member of the glutathione-S-transferase superfamily. Chem Res Toxicol. 2001, 14: 1051-1057. 10.1021/tx010052h.

    Article  CAS  Google Scholar 

  52. Ahsan H, Chen Y, Kibriya MG, Slavkovich V, Parvez F, Jasmine F, Gamble MV, Graziano JH: Arsenic metabolism, genetic susceptibility, and risk of premalignant skin lesions in Bangladesh. Canc Epidemiol Biomark Prev. 2007, 16: 1270-1278. 10.1158/1055-9965.EPI-06-0676.

    Article  CAS  Google Scholar 

  53. Marahatta SB, Punyarit P, Bhudisawasdi V, Paupairoj A, Wongkham S, Petmitr S: Polymorphism of glutathione S-transferase omega gene and risk of cancer. Cancer Lett. 2006, 236: 276-281. 10.1016/j.canlet.2005.05.020.

    Article  CAS  Google Scholar 

  54. De CS, Ghosh P, Sarma N, Majumdar P, Sau TJ, Basu S, Roychoudhury S, Ray K, Giri AK: Genetic variants associated with arsenic susceptibility: study of purine nucleoside phosphorylase, arsenic (+3) methyltransferase, and glutathione S-transferase omega genes. Environ Health Perspect. 2008, 116: 501-505.

    Google Scholar 

  55. Dong LM, Potter JD, White E, Ulrich CM, Cardon LR, Peters U: Genetic susceptibility to cancer: the role of polymorphisms in candidate genes. JAMA. 2008, 299: 2423-2436. 10.1001/jama.299.20.2423.

    Article  CAS  Google Scholar 

  56. Moore LE, Malats N, Rothman N, Real FX, Kogevinas M, Karami S, Garcia-Closas R, Silverman D, Chanock S, Welch R, Tardon A, Serra C, Carrato A, Dosemeci M, Garcia-Closas M: Polymorphisms in one-carbon metabolism and trans-sulfuration pathway genes and susceptibility to bladder cancer. Int J Cancer. 2007, 120: 2452-2458. 10.1002/ijc.22565.

    Article  CAS  Google Scholar 

  57. Safarinejad MR, Shafiei N, Safarinejad S: Genetic susceptibility of methylenetetrahydrofolate reductase (MTHFR) gene C677T, A1298C, and G1793A polymorphisms with risk for bladder transitional cell carcinoma in men. Med Oncol. 2011, S1: S398-S412.

    Article  Google Scholar 

  58. Chung CJ, Pu YS, Su CT, Chen HW, Huang YK, Shiue HS, Hsueh YM: Polymorphisms in one-carbon metabolism pathway genes, urinary arsenic profile, and urothelial carcinoma. Canc Causes Contr. 2010, 21: 1605-1613. 10.1007/s10552-010-9589-3.

    Article  Google Scholar 

  59. Lin J, Spitz MR, Wang Y, Schabath MB, Gorlov IP, Hernandez LM, Pillow PC, Grossman HB, Wu X: Polymorphisms of folate metabolic genes and susceptibility to bladder cancer: a case–control study. Carcinogenesis. 2004, 25: 1639-1647. 10.1093/carcin/bgh175.

    Article  CAS  Google Scholar 

  60. Vahter M: Genetic polymorphism in the biotransformation of inorganic arsenic and its role in toxicity. Toxicol Lett. 2000, 112–113: 209-217.

    Article  Google Scholar 

  61. Moore LE, Smith AH, Eng C, DeVries S, Kalman D, Bhargava V, Chew K, Ferreccio C, Rey OA, Hopenhayn C, Biggs ML, Bates MN, Waldman FM: P53 alterations in bladder tumors from arsenic and tobacco exposed patients. Carcinogenesis. 2003, 24: 1785-1791. 10.1093/carcin/bgg136.

    Article  CAS  Google Scholar 

  62. Kwong RC, Karagas MR, Kelsey KT, Mason RA, Tanyos SA, Schned AR, Marsit CJ, Andrew AS: Arsenic exposure predicts bladder cancer survival in a US population. World J Urol. 2010, 28: 487-492. 10.1007/s00345-009-0477-y.

    Article  CAS  Google Scholar 

  63. Chen CH, Chiou HY, Hsueh YM, Chen CJ, Yu HJ, Pu YS: Clinicopathological characteristics and survival outcome of arsenic related bladder cancer in Taiwan. J Urol. 2009, 181: 547-552.

    Article  CAS  Google Scholar 

Download references


We would like to thank Dr. James Montie, former chair of the Department of Urology, for his support of this project, the University of Michigan Medical School and the University of Michigan Comprehensive Cancer Center. This work was supported by the National Cancer Institute grants at the National Institutes of Health R01 CA9600210; K07 CA127214; K07 CA127214-02 S1.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Jennifer L Beebe-Dimmer.

Additional information

Competing interests

The authors declare they have no competing interests.

Authors’ contributions

JBD participated in the design and conception of the study, selection of genes, and collection of saliva specimens (with JN and JM), provided guidance during analyses and was primarily responsible for manuscript preparation. PI, SM and KA were responsible for the isolation of DNA from specimens, genotyping and sample quality control procedures. GK and EML were responsible for the conduct of the statistical analyses. AGS, JN, JM, DS, KA participated in the design and conception of the study as well as manuscript preparation and editing. KA supervised all of the genotyping efforts in her lab. All authors have read and approved of final manuscript.

Electronic supplementary material


Additional file 1 :Table S1. Pairwise linkage disequilibrium analysis of 5 selected single nucleotide polymorphisms in As3MT among participants of University of Michigan Bladder Cancer Study. (DOC 32 KB)

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Beebe-Dimmer, J.L., Iyer, P.T., Nriagu, J.O. et al. Genetic variation in Glutathione S-Transferase Omega-1, Arsenic Methyltransferase and Methylene-tetrahydrofolate Reductase, arsenic exposure and bladder cancer: a case–control study. Environ Health 11, 43 (2012).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: