Search tips
Search criteria 


Logo of plosgenPLoS GeneticsSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)View this Article
PLoS Genet. 2012 February; 8(2): e1002522.
Published online 2012 February 23. doi:  10.1371/journal.pgen.1002522
PMCID: PMC3285587

Genome-Wide Association Study Identifies Chromosome 10q24.32 Variants Associated with Arsenic Metabolism and Toxicity Phenotypes in Bangladesh

Mark I. McCarthy, Editor


Arsenic contamination of drinking water is a major public health issue in many countries, increasing risk for a wide array of diseases, including cancer. There is inter-individual variation in arsenic metabolism efficiency and susceptibility to arsenic toxicity; however, the basis of this variation is not well understood. Here, we have performed the first genome-wide association study (GWAS) of arsenic-related metabolism and toxicity phenotypes to improve our understanding of the mechanisms by which arsenic affects health. Using data on urinary arsenic metabolite concentrations and approximately 300,000 genome-wide single nucleotide polymorphisms (SNPs) for 1,313 arsenic-exposed Bangladeshi individuals, we identified genome-wide significant association signals (P<5×10−8) for percentages of both monomethylarsonic acid (MMA) and dimethylarsinic acid (DMA) near the AS3MT gene (arsenite methyltransferase; 10q24.32), with five genetic variants showing independent associations. In a follow-up analysis of 1,085 individuals with arsenic-induced premalignant skin lesions (the classical sign of arsenic toxicity) and 1,794 controls, we show that one of these five variants (rs9527) is also associated with skin lesion risk (P = 0.0005). Using a subset of individuals with prospectively measured arsenic (n = 769), we show that rs9527 interacts with arsenic to influence incident skin lesion risk (P = 0.01). Expression quantitative trait locus (eQTL) analyses of genome-wide expression data from 950 individual's lymphocyte RNA suggest that several of our lead SNPs represent cis-eQTLs for AS3MT (P = 10−12) and neighboring gene C10orf32 (P = 10−44), which are involved in C10orf32-AS3MT read-through transcription. This is the largest and most comprehensive genomic investigation of arsenic metabolism and toxicity to date, the only GWAS of any arsenic-related trait, and the first study to implicate 10q24.32 variants in both arsenic metabolism and arsenical skin lesion risk. The observed patterns of associations suggest that MMA% and DMA% have distinct genetic determinants and support the hypothesis that DMA is the less toxic of these two methylated arsenic species. These results have potential translational implications for the prevention and treatment of arsenic-associated toxicities worldwide.

Author Summary

Exposure to arsenic through drinking water is a serious public health issue in many countries, including Bangladesh and the United States. Although there is substantial inter-individual variation in arsenic metabolism and toxicity, the biological basis of this variation is not well understood. Here, we have conducted the first genome-wide association study of arsenic-related traits within a unique population cohort of arsenic-exposed Bangladeshi individuals. Using data on 1,313 well-characterized individuals, we identify multiple association signals for urinary arsenic metabolite concentrations in the 10q24.32 regions, near the AS3MT (arsenite methyltransferase) gene. In a subsequent analysis of >2,000 individuals, we show for the first time that variants that influence arsenic metabolism can also influence risk for arsenical skin lesions (the classical sign of arsenic toxicity) through interaction with arsenic exposure. Using array-based genome-wide gene expression data, we show that several of our lead genetic variants are associated with expression of AS3MT and neighboring gene C10orf32, providing a potential mechanism by which 10q24.32 variants influence arsenic metabolism and toxicity. Knowledge of variation in this region and associated biological processes could be used to develop intervention and pharmacological strategies aimed at preventing large numbers of arsenic-related deaths in arsenic-exposed populations.


Over 100 million individuals worldwide are exposed to arsenic through drinking water, including approximately 56 million people in Bangladesh [1] and 13 million in the United States [2]. Arsenic is a class I human carcinogen, and chronic exposure to high levels of arsenic (>300 µg/L) is associated with substantial increased risk for a wide array of diseases including cancers of the lung [3], bladder [4], liver [5], skin [6], and kidney [7], [8], as well as neurological [9], [10] and cardiovascular [11] diseases. Emerging evidence suggests that arsenic may have adverse effects on health even at concentrations as low as 10–50 µg/L, as recent studies in Bangladesh have observed dose-response relationships with mortality [12], [13] and arsenical skin lesion risk [14] in populations with low to moderate arsenic exposure over many years. Arsenical skin lesions are a classical sign of arsenic toxicity, an indicator of susceptibility to arsenic-related disease, and a precursor to arsenic-induced skin cancers [6]. Once individuals are chronically exposed to arsenic, risk for arsenic-related diseases and mortality remains high for several decades even after cessation of exposure [15], [16].

Consumed arsenic enters the blood as AsV and AsIII, known collectively as inorganic arsenic (iAs). Once consumed, iAs is methylated using S-Adenosyl methionine (SAM) as the methyl donor, producing monomethylarsonic acid (MMA) and then dimethylarsinic acid (DMA). MMA is believed to be the more toxic of these metabolites, with the DMA/MMA ratio showing an inverse association with arsenic toxicity in several studies [17][20] and DMA being more readily excreted in urine and expelled from the body. Arsenic metabolite concentrations are often expressed as percentages of all arsenic species present in urine (i.e., iAs%, MMA%, DMA%) or as ratios that reflect methylation efficiency (e.g., DMA%/MMA%, MMA%/iAs%).

There is considerable inter-individual variation in arsenic metabolism, as some individuals are able to methylate, and thus excrete, arsenic more efficiently than others [21], [22]. Similarly, because high inter-individual variability in toxicity is observed among individuals with similar levels of exposure to arsenic, genetic susceptibility factors for arsenical skin lesions are believed to exist [23].

In light of the enormous global health impact of arsenic exposure and the remarkable inter-individual variability in arsenic metabolism and toxicity, we performed the first genome-wide association study (GWAS) of common arsenic-related phenotypes. We identified multiple genetic variants in the 10q24.32 region near AS3MT (arsenite methyltransferase, previously known as CYT19) that show robust associations with urinary concentrations of arsenic metabolites, risk for arsenical skin lesions, and local gene expression, including transcript levels of AS3MT.


GWAS of arsenic metabolites

We assessed genome-wide associations for the three arsenic metabolites measured in urine (iAs%, MMA%, and DMA%) using high-quality data on 259,597 single-nucleotide polymorphisms (SNPs) from 1,313 individuals randomly selected from a large population-based cohort of Bangladeshi individuals exposed to a wide range of arsenic concentrations through drinking water. Associations were assessed using mixed linear models [24] to account for existence of related individuals in our sample (Figure S1). The strongest association signals, genome-wide, for both DMA% and MMA% were in the 10q24.32 region (P<5×10−8) (Figures S2 and S3), which contains the AS3MT gene and substantial LD spanning ~1 Mb (Figure S4).

For DMA%, the strongest 10q24.32 association was for rs9527 (P = 2.7×10−9; Figure 1). After conditioning on rs9527, a strong residual association signal remained (rs11191527; P = 8.0×10−8), the strength of which was weaker without adjustment for rs9527 (P = 2.3×10−5) due to mild LD between these SNPs (D′ = 0.26, r2 = 0.03 in our data; D′ = 0.27; r2 = 0.03 in HapMap GIH). After conditioning on both SNPs, there was very little evidence of additional association in the region. Analyses of imputed and measured genotypes produced the same two association signals, but with imputed SNPs rs3740394 and rs17115073 showing slightly stronger association than rs9527 and rs11191527, respectively (Figure S5).

Figure 1
Multiple variants in the 10q24.32 region show independent associations with DMA% (n = 1,313).

The strongest association observed for MMA% was rs4919694 (P = 2.9×10−8) (Figure 2). After conditioning on rs4919694, residual association was still observed (rs4290163; P = 7.0×10−5). This association is much weaker without adjustment for rs4919694 (P = 0.03) due to LD between rs4919694 and rs4290163 (D′ = 0.80, r2 = 0.09 in our data; D′ = 0.80, r2 = 0.04 in HapMap GIH). Aftern conditioning on both SNPs, residual association was observed for rs11191659 (P = 0.0009), a SNP in moderate LD with rs9527, the top SNP from the %DMA analysis (D′ = 0.66, r2 = 0.23 in our data; D′ = 0.82, r2 = 0.30 in HapMap GIH). Conditioning on all three SNPs eliminated the 10q24.32 association signal. Imputation of unobserved genotypes in the region did not reveal associations stronger than those observed for the measured genotypes (Figure S6). Multivariate models for %DMA and %MMA including all five of the above-mentioned SNPs are described in Table 1. Outside of the 10q24.32 region, there was no genome-wide significant (P<5×10−8) association signal for DMA% or MMA%.

Figure 2
Multiple variants in the 10q24.32 region show independent associations with MMA% (n = 1,313).
Table 1
Multivariate associations between arsenic metabolites and genotyped SNPs in the 10q24.32 region showing the strongest univariate associations with DMA% and MMA% (n = 1,313).

The 10q24.32 association results for the DMA%/MMA% ratio (the “secondary methylation index”, log-transformed), were very similar to the MMA% results, as these phenotypes were strongly correlated (r = −0.84; Table S1). Associations for 10q24.32 SNPs with iAs% and MMA%/iAs% (the “primary methylation index” (PMI) log-transformed) were much weaker than for DMA% and MMA%; The strongest association in the 10q24.32 region observed for iAs% was rs9527 (P = 0.0009) and no association of P<0.001 was observed for log(PMI). In genome-wide analyses of iAs% and PMI, no SNP reached genome-wide significance (Figure S7).

Association of SNPs with skin lesions risk

Because variants influencing arsenic metabolism may alter susceptibility to arsenic toxicity, we investigated the roles of metabolite-associated SNPs in arsenic-induced premalignant skin lesions, the hallmark of chronic arsenic toxicity. For our five lead SNPs, we tested association with skin lesion status among 1,085 skin lesion cases and 1,794 population controls, using the ROADTRIPS method that was developed for case-control association testing in the presence of cryptic relatedness [25]. The rs9527 allele associated with decreased DMA% (A) was associated with increased skin lesion risk (P = 0.0005), consistent with the hypothesis that DMA is less toxic than MMA (Table 2). rs11191659 showed suggestive association (P = 0.02), also consistent with this hypothesis.

Table 2
Association between the 10q24.32 genotyped variants and arsenical skin lesion risk and SNP-arsenic interaction estimates.

To confirm that these associations with skin lesions were due to gene-arsenic interaction, we tested the interaction between rs9527 and arsenic exposure using a subset of 69 incident skin lesion cases and 700 controls with prospectively-measured arsenic exposure (measured in both water and urine at baseline, prior to skin lesion incidence and arsenic mitigation efforts [26]). We found SNP-arsenic interaction for both rs9527 (multiplicative interaction P = 0.01; additive interaction P = 0.004) and rs11191659 (multiplicative interaction P = 0.02; additive interaction P = 0.001), where water arsenic exposure showed stronger association with skin lesions in the presence of the risk allele (Table 2 and Table S2). There were no significant main effects for either of these SNPs in the context of models that included SNP-arsenic interaction terms.

For the subset of individuals with available genotype, arsenic metabolite, and skin lesion data (82 cases, 1211 controls), DMA% showed evidence of partial mediation of the association between rs9527 and skin lesions (accounting for 13% of the observed association).

eQTL analyses of the 10q24.32 region

To investigate the role of our lead SNPs in gene regulation, used genome-wide expression data derived from lymphocyte RNA obtained at baseline for 950 participants (Illumina HumanHT-12 array) and examined SNP-expression associations for all 30 genes in the 10q24.32 LD region (Table S3). Several of our lead SNPs showed association with AS3MT expression at P<5×10−5 (rs4919694, rs9527, rs4290163). However, after examining associations for all SNPs in this region, C10orf32 intronic SNP rs7096169 showed the strongest association with AS3MT expression (P = 8×10−12; Figure 3 and Figure S8), and conditioning on rs7096169 eliminated the eQTL signal. rs7096169 was not one of our lead SNPs, but it was associated with DMA% (P = 0.001; MMA% P = 0.28). Interestingly, the rs9527 risk allele (A) was associated with decreased C10orf32 expression (P = 2.6×10−41; Figure S8), the strongest eQTL signal for C10orf32 expression in the region (Figure 3) and the strongest genome-wide eQTL effect for rs9527 (Figure S9). C10orf32 is ~4 kb upstream of AS3MT, and these genes are involved in C10orf32-AS3MT read-through transcription, producing a transcript that is a candidate for nonsense-mediated mRNA decay. Thus, it is possible that the eQTL signal observed for C10orf32 represents a regulatory mechanism that influences read-through transcript production. After conditioning on rs9527, the residual eQTL signal was best represented by rs11083790 (P = 10−5). Conditioning on both SNPs eliminated the eQTL signal. Interestingly, C10orf32 expression was also associated with arsenic exposure (measured as total arsenic in urine, collected at the same time as blood; P = 0.001), while AS3MT expression was not (P = 0.37). None of our lead SNPs modified the association between arsenic exposure and C10orf32 expression.

Figure 3
Variants in the 10q24.32 region are associated with transcript levels of AS3MT and C10orf32.

Our lead SNPs were also associated with USMG5 expression (Table S3), a gene ~500 kb downstream of AS3MT, but these associations appear to be due to moderate LD with downstream variants showing very strong association with USMG5 expression (e.g., rs12220267; P = 10−210; Figure S8).


The role of AS3MT in arsenic metabolism has been described [27], and several prior studies have evaluated associations between candidate AS3MT variants arsenic-related traits in Bangladesh and elsewhere [28][34]. A recent review [35] highlighted two AS3MT SNPs, rs11191439 (Met287Thr) and rs3740393 (intronic), as being consistently related to arsenic metabolism across diverse populations. The most recent and comprehensive Bangladeshi study of AS3MT SNPs [28] reported three association signals for arsenic metabolites, best represented by HapMap3 SNPs rs1046778 (for MMA%), rs11191439 (DMA% and iAs%), and rs3740390 (DMA% and iAs%), a proxy for rs3740393 (r2 = 0.91). After imputation, we were able to replicate rs11191439 (DMA% P = 4.2×10−6; MMA% P = 5.8×10−7) and rs1046778 (MMA% P = 8.9×10−7; DMA% P = 0.0002), which were strongly correlated with lead SNPs rs4919694 (r2 = 0.69) and rs4290163 (r2 = 0.63), respectively. After conditioning on our lead SNPs, these associations were no longer significant. The evidence for rs3740390 was less convincing (DMA% P = 0.54; MMA% P = 0.007), as this SNP was not strongly correlated with any of our lead SNPs (Figure S10). We identified two novel 10q24.32 association signals, represented by rs9527 and rs11191527, which were not strongly correlated with any previously-reported SNP (Figure S10). These SNPs were likely missed in prior studies due to limited coverage of the SNPs in this region.

The identities of the functional variants in this region remain unclear. rs9527 lies in the 5′ UTR of C10orf32, a transcription factor binding region (GATA-1 and TAL1 (SC-12984)) and a DNase hypersensitivity site. If causal, rs9527 could also exert its effects through regulation of AS3MT-C10orf32 read-through transcription. However, the LD block represented by rs9527 includes transcription factor binding site SNP rs12416687 and miRNA SNPs rs11191401, rs12573077, rs7904252, and rs9527. Detailed information on potential functional variants from HapMap3 (GIH) for each of the 5 SNPs identified is contained in Tables S4, S5, S6, S7, S8. However, genetic variation in this population has not been comprehensively characterized (especially rare variation), and the underlying functional variants may not be present in HapMap3. It is also possible that the underlying causal variants have implications for surrounding genes. For example, rs4919694 and rs11191527 are intronic SNPs within the CNNM2 gene, which is involved in magnesium reabsorption by the kidney [36]. It is possible that magnesium and iAs interact [37], [38], influencing the amount of free arsenic available for methylation.

To our knowledge, this study is the largest genetic association study of arsenic metabolites to date, the only GWAS of arsenic-related traits, the first study to implicate 10q24.32 SNPs in both arsenic metabolism and arsenical skin lesion risk, and one of the earliest GWAS conducted in the developing country setting. Our results suggest that MMA% and DMA% have distinct genetic determinants and highlight the importance of conditional analyses, as LD among alleles with opposing effects can mask associations in univariate analyses. The associations observed in this study are likely due to the effects of unmeasured, potentially rare variants in LD with the measured SNPs and/or substantial allelic heterogeneity, whereby multiple 10q24.32 variants influence arsenic metabolism.

Considering the substantial LD in this region [39], the variation in allele frequencies and LD patterns among the various arsenic-exposed populations under study [40], and the apparent allelic heterogeneity with respect to arsenic metabolism, future DNA sequencing studies are needed to help identify causal variants in the 10q24.32 region. Identifying these variants will help clarify the links between the association signals observed for %DMA, %MMA, and AS3MT/C10orf32 expression. These association signals appear largely independent in our dataset, but perhaps there are underlying causal variants that influence all of these phenotypes. Developing a better understanding the effects of functional variation related to AS3MT will also provide a more nuanced understanding of the biology of arsenic methylation, which can in turn help us better understand how variation in methylation efficiency affects health. Finally, knowledge of this causal variation and the methylation processes that they influence could potentially be exploited for intervention strategies that aim to prevent large numbers of deaths arsenic-exposed populations, by defining susceptibility subgroups and exploiting the biological processes uncovered by genomics for developing pharmacological treatments.

Materials and Methods

Study descriptions

The DNA samples genotyped in this study were obtained at baseline recruitment from individuals participating in one of the following studies: The Health Effects of Arsenic Longitudinal Study (HEALS) [41] or the Bangladesh Vitamin E and Selenium Trial (BEST) [42]. GWAS analyses of arsenic metabolites were conducted using urinary arsenic metabolite and SNP data on 1,313 individuals randomly selected from the HEALS study. Analyses of skin lesion data were conducted using genotype data from 1,085 skin lesion cases and 1,794 controls drawn from both studies, including the 1,313 HEALS individuals with metabolite data. Skin lesion cases included individuals with keratosis, melanosis, and leukomelanosis. Gene expression analyses were based on lymphocyte RNA extracted at baseline recruitment for the first 950 BEST participants. A summary of these overlapping sets of samples is provided in Figure S11.

The Health Effects of Arsenic Longitudinal Study (HEALS [41]) is a prospective investigation of health outcomes associated with arsenic exposure through drinking water in a cohort of adults in Araihazar, Bangladesh, a rural area east of the capital city, Dhaka. Between October 2000 and May 2002, we recruited healthy married individuals (age 18–75 years) who were residents of the study area for at least five years and primarily consumed drinking water from a local well. We enumerated 65,876 individuals residing in Araihazar, from which we identified a sampling frame of 14,828 eligible residents. Of these 14,828 individuals, 11,746 men and women were enrolled. During 2006–2008, additional recruitment of 8,287 participants from the same underlying source population expanded the cohort size to over 20,000 individuals. All 5,966 wells in the study area were tested for arsenic using graphite furnace atomic absorption spectrometry and individuals reported the primary well from which they drank. At baseline, trained study physicians, blinded to the arsenic measurements, conducted in-person interviews and clinical evaluations and collected spot urine and blood samples from participants in their homes using structured protocols. Similar in-person follow-up interviews were conducted biennially for the entire cohort during the following periods: follow-up 1 during September 2002 to May 2004, follow-up 2 during June 2004 to August 2006, and follow-up 3 during January 2007 to February 2009. At baseline and each follow-up interview, a structured protocol was used to ascertain skin lesions by the study physicians, who had undergone training for the detection and diagnosis of skin lesions [43]. The study protocol was approved by the Institutional Review Boards of The University of Chicago, Columbia University, and the Bangladesh Medical Research Council. Informed consent was obtained from all participants.

The Bangladesh Vitamin E and Selenium Trial (BEST) is a 2×2 factorial randomized chemoprevention trial evaluating the long-term effects of vitamin E and selenium supplementation on non-melanoma skin cancer (NMSC) risk. BEST participants are residents of Araihazar (the same geographic area as HEALS participants with 132 overlapping participants), Matlab, and surrounding areas. BEST uses many of the same study protocols as does HEALS, especially arsenic exposure assessment and biospecimen collection protocols. All participants were required to have existing arsenic-related skin lesions to be eligible. A total of 7,000 individuals have been randomized to one of the four treatment arms: vitamin E only (100 IU/day), L-selenomethionine only (200 µg/day), both vitamin E and selenium, and placebo. Participants have been actively followed for 6 years and systematic ascertainment of histopathologically-confirmed NMSC has been conducted (including BCC and SCC). For all participants, biological samples, including all fractions of blood including DNA and RNA, urine, toenails, and tumor samples have been collected at baseline, along with clinical and covariate data, creating a biological and data repository that is available for research purposes. The study protocol was approved by the Ethical Review Committee of International Center for Diarrheal Disease Research, Bangladesh, the Bangladesh Medical Research Council, and the Institutional Review Boards of The University of Chicago and Columbia University. Informed consent was obtained from all participants.

In each study, urinary arsenic was measured using graphite furnace atomic absorption spectrometry in a single laboratory [44]. Urinary creatinine was measured by a colorimetric diagnostics kit (Sigma, St Louis, MO, USA). Total urinary arsenic concentration was divided by creatinine to obtain creatinine-adjusted total arsenic concentration (µg/g creatinine) [45]. Urinary arsenic metabolites (arsenobetaine, arsenocholine, arsenite, arsenate, monomethylarseonous acid, and dimethylarsenic acid) were distinguished as described by Ahsan et al. [17], using a high-performance liquid chromatography method for separation of arsenic metabolites, followed by detection using inductively coupled plasma-mass spectrometry with dynamic reaction cell. The percentage of iAs, MMA and DMA in total arsenic was calculated after subtracting asenobetaine and arsenocholine (i.e., nontoxic organic arsenic from dietary sources). Because these metabolites lie on the same biological pathway and are expressed as a percentage of arsenic species, their values show substantial correlation (Table S1).

Genotyping and quality control

For BEST samples, DNA extraction was carried out from the whole blood using the QIAamp 96 DNA Blood Kit (cat # 51161) from Qiagen, Valencia, USA. For HEALS samples, DNA was extracted from clot blood using Flexigene DNA kit (Cat # 51204) from Qiagen. Concentration and quality of all extracted DNA were checked by Nanodrop 1000. As starting material, 250 ng of DNA was used on the Illumina Infinium HD SNP array according to Illumina's protocol. Samples were processed on HumanCytoSNP-12 v2.1 chips with 299,140 markers and read on the BeadArray Reader. Image data was processed in BeadStudio software to generate genotype calls.

Prior to genotype QC, our genotype data consisted of 2,920 samples typed for 299,140 SNPs. First, we removed DNA samples with very poor call rates (<90%; n = 8) and SNPs that were poorly called (<90%) or monomorphic (n = 39,276). Individuals with gender mismatches were removed (n = 10), as were technical replicate DNA samples run to assure high genotyping accuracy (n = 21). No individuals had outlying autosomal heterozygosity or inbreeding values. After inspecting distributions of SNP and samples call rates, we excluded samples with call rates <97% (n = 2) and SNPs with call rates <95% (n = 103). SNPs with HWE p-values<10−7 were excluded (n = 164). This QC resulted in 2,879 individuals with high-quality genotype data for 259,597 SNPs. All QC was performed using PLINK [46].

Gene expression

RNA was extracted from mononuclear cells preserved in buffer RLT, stored at −86°C using RNeasy Micro Kit (cat# 74004) from Qiagen, Valencia, USA. Concentration and quality of all extracted RNA were checked on Nanodrop 1000. cRNA synthesis was done from 250 ng of RNA using Illumina TotalPrep 96 RNA Amplification kit. As recommended by Illumina we used 750 ng of cRNA on HumanHT-12-v4 for gene expression. The chip contains a total of 47,231 probes covering 31,335 genes.

Statistical analyses and software tools

Pair-wise kinship coefficients were estimated using PLINK [46] and their distribution is shown in Figure S1. To assess population structure that was unrelated to the relative pairs present in our dataset, we removed one individual from each related pair (kinship coefficient >0.05) and assessed population structure in this dataset of 403 individuals using principal components analysis as implemented in EIGENSTRAT [47]. We found very little evidence of population stratification (Figure S12), with the eigenvalues from the first ten principle components being between 1.123 and 1.184. All SNP association tests for urinary metabolites were conducted using a mixed model that accounted for cryptic relatedness as implemented in EMMAX [24] (rather than principle components), adjusting for water arsenic, sex, and age. All regional association plots were generated using LocusZoom [48]. Association testing for skin lesion status was conducting using PLINK [46] and the ROADTRIPS [25] software developed for case-control association testing in samples with unknown population and pedigree structure. We conducted local imputation for the 10q24.32 region using MACH, the GIH reference panel, and imputation parameters suggested by the developers [49]. The estimated genotype and allele error rates were 0.034 and 0.017, respectively. LD structure in the 10q24.32 region was visualized using Haploview [50]. Information on the potential functional consequences of SNPs in the 10q24.32 regions was obtained using the NIEHS's SNPinfo Web Server [51]. Interaction analyses was conducted using only HEALS incident cases (n = 69) and controls (n = 701). For BEST participants and some HEALS participants arsenic exposure (based on water and urine) was not measured prior to arsenic mitigation efforts [26], so the measured exposure status for these individuals is not likely to reflect long-term arsenic exposure status. Interactions were tested using the SAS 9.2 PROC MIXED procedure, using the “bn” matrix derived using EMMAX. To assess mediation of the association between SNPs and skin lesions, we used the “proportion explained” (PE) equation for odds ratios (PEOR = (ORxy−ORxy|m)/(ORxy−1) where x is an exposure, y is a binary outcome, and m is a potential mediating factor [52]). Genome-wide eQTL analysis for our five lead SNPs was performed using the significance of microarray method as implemented in BRB Array Tools. Promising eQTL effects were then examined using EMMAX as described above, treating the expression values as a quantitative trait. In a similar fashion, arsenic exposure was tested for association with expression traits of interest and for interaction with SNPs in relation to expression traits using PROC MIXED.

Supporting Information

Figure S1

Distribution of all pair-wise kinship coefficients among 2,879 Bangladeshi individuals with measured genome-wide SNP data. Kinship values are truncated at 0.05. The observed clusters of observations centered at 0.5 0.25, and 0.125 represent full siblings or parent-offspring, half siblings, and, and first cousins pairs, respectively. One individual from each pair of twins or duplicate samples (kinship coefficient of 1.0) was removed from the analysis dataset.


Figure S2

GWAS Results for DMA% (including Manhattan plot, QQ plot, and the strongest associated SNPs).


Figure S3

GWAS Results for MMA% (including Manhattan plot, QQ plot, and the strongest associated SNPs).


Figure S4

Summary of linkage disequilibrium (LD) in the 10q24.32 region. LD data for all SNPs in the GIH HapMap3 panel are shown above (A) and the data for SNPs typed in this study are shown below (B). Dark red squares represent a D′ value near 1 and white squares represent a D′ value near zero.


Figure S5

DMA% associations results for imputed and genotyped SNPs in the 10q24.32 region (n = 1,313). P-values were generated using mixed-models adjusted for age, sex, and water arsenic concentration. The strongest associated SNP is labeled in each panel. The top panel shows the overall association results. The second panel shows P-values from models that are adjusted for rs3740394. The third panel shows P-values from models adjusted for both rs3740394 and rs17115073.


Figure S6

MMA% associations results for imputed and genotyped SNPs in the 10q24.32 region (n = 1,313). P-values were generated using mixed-models adjusted for age, sex, and water arsenic concentration. The strongest associated SNP is labeled in each panel. The top panel shows the overall association results. The second panel shows P-values from models that are adjusted for rs4919694. The third panel shows P-values from models adjusted for both rs4919694 and rs4290163.


Figure S7

Q-Q plots for genome-wide association scans of the iAs% and the primary methylation index (PMI = MMA%/iAs%). Results are based on 1,310 samples and 259,597 SNPs. The EMMAX model is adjusted for sex and water arsenic.


Figure S8

Cis-eQTL signals in the 10q24.32 region. Expression values for AS3MT, C10orf32, and USMG5 are shown by the minor allele count for rs7096160, rs9527, and rs12220267, respectively. Mean expression values are shown as dotted lines.


Figure S9

Association for our 5 lead SNPs with genome-wide transcript levels. Genome-wide eQTL analysis was performed using the Significance of Microarray method as implemented in BRB Array Tools. The first and second most strongly associated transcripts (in red) are C10orf32 and USMG5, respectively. Both are in the 10q24.32 region.


Figure S10

Linkage Disequilibrium between lead SNPs and previously reported variants. Our lead SNPs are shown in boxes and the variants representing the signals previously reported in a Bangladeshi study (Engstrom et al. [30]) are shown in circles.


Figure S11

An overview of the participants and samples used in this work.


Figure S12

Scatter plot of the first two principle components for 403 unrelated study participants (no pair-wise kinship value>0.05).


Table S1

Pair-wise correlations among the arsenic-related urinary phenotypes examined in this study (n = 1,333).


Table S2

Regression models for interaction between arsenic exposure and rs9257 in relation to skin lesion risk (69 skin lesion cases, 701 controls).


Table S3

P-values from association test for our 5 lead 10q24.32 SNPs and expression values for all genes in the 10q24.32 LD region (n = 950 individuals).


Table S4

Functional information for SNPs in LD with rs9527.


Table S5

Functional information for SNPs in LD with rs11191527.


Table S6

Functional information for SNPs in LD with rs4919694.


Table S7

Functional information for SNPs in LD with rs3740394.


Table S8

Functional information for SNPs in LD with rs11191659.



The authors would like to sincerely thank all study participants and research staff that have contributed to this research.


The authors have declared that no competing interests exist.

This work was supported by NIH Grants P42ES010349, R01CA102484, R01CA107431, and P30CA014599 and by Department of Defense grant W81XWH-10-1-0499. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


1. Smith AH, Lingas EO, Rahman M. Contamination of drinking-water by arsenic in Bangladesh: a public health emergency. Bull World Health Organ. 2000;78:1093–1103. [PubMed]
2. United States Environmental Protection Agency. Fact Sheet: drinking water standard for arsenic. 2009. EPI815-R-00-015.
3. Celik I, Gallicchio L, Boyd K, Lam TK, Matanoski G, et al. Arsenic in drinking water and lung cancer: a systematic review. Environ Res. 2008;108:48–55. [PubMed]
4. Mink PJ, Alexander DD, Barraj LM, Kelsh MA, Tsuji JS. Low-level arsenic exposure in drinking water and bladder cancer: a review and meta-analysis. Regul Toxicol Pharmacol. 2008;52:299–310. [PubMed]
5. Liu J, Waalkes MP. Liver is a target of arsenic carcinogenesis. Toxicol Sci. 2008;105:24–32. [PMC free article] [PubMed]
6. Yu HS, Liao WT, Chai CY. Arsenic carcinogenesis in the skin. J Biomed Sci. 2006;13:657–666. [PubMed]
7. Chen CJ, Chen CW, Wu MM, Kuo TL. Cancer potential in liver, lung, bladder and kidney due to ingested inorganic arsenic in drinking water. Br J Cancer. 1992;66:888–892. [PMC free article] [PubMed]
8. Yuan Y, Marshall G, Ferreccio C, Steinmaus C, Liaw J, et al. Kidney cancer mortality: fifty-year latency patterns related to arsenic exposure. Epidemiology. 2010;21:103–108. [PubMed]
9. Brinkel J, Khan MH, Kraemer A. A systematic review of arsenic exposure and its social and mental health effects with special reference to Bangladesh. Int J Environ Res Public Health. 2009;6:1609–1619. [PMC free article] [PubMed]
10. Vahidnia A, van der Voet GB, de Wolff FA. Arsenic neurotoxicity–a review. Hum Exp Toxicol. 2007;26:823–832. [PubMed]
11. States JC, Srivastava S, Chen Y, Barchowsky A. Arsenic and cardiovascular disease. Toxicol Sci. 2009;107:312–323. [PMC free article] [PubMed]
12. Argos M, Kalra T, Rathouz PJ, Chen Y, Pierce B, et al. Arsenic exposure from drinking water, and all-cause and chronic-disease mortalities in Bangladesh (HEALS): a prospective cohort study. Lancet. 2010;376:252–258. [PMC free article] [PubMed]
13. Sohel N, Persson LA, Rahman M, Streatfield PK, Yunus M, et al. Arsenic in drinking water and adult mortality: a population-based cohort study in rural Bangladesh. Epidemiology. 2009;20:824–830. [PubMed]
14. Argos M, Kalra T, Pierce BL, Chen Y, Parvez F, et al. A prospective study of arsenic exposure from drinking water and incidence of skin lesions in Bangladesh. Am J Epidemiol. 2011;174:185–194. [PMC free article] [PubMed]
15. Yang CY, Chiu HF, Chang CC, Ho SC, Wu TN. Bladder cancer mortality reduction after installation of a tap-water supply system in an arsenious-endemic area in southwestern Taiwan. Environ Res. 2005;98:127–132. [PubMed]
16. Chang CC, Ho SC, Tsai SS, Yang CY. Ischemic heart disease mortality reduction in an arseniasis-endemic area in southwestern Taiwan after a switch in the tap-water supply system. J Toxicol Environ Health A. 2004;67:1353–1361. [PubMed]
17. Ahsan H, Chen Y, Kibriya MG, Slavkovich V, Parvez F, et al. Arsenic metabolism, genetic susceptibility, and risk of premalignant skin lesions in Bangladesh. Cancer Epidemiol Biomarkers Prev. 2007;16:1270–1278. [PubMed]
18. Kile ML, Hoffman E, Rodrigues EG, Breton CV, Quamruzzaman Q, et al. A pathway-based analysis of urinary arsenic metabolites and skin lesions. Am J Epidemiol. 2011;173:778–786. [PMC free article] [PubMed]
19. Lindberg AL, Rahman M, Persson LA, Vahter M. The risk of arsenic induced skin lesions in Bangladeshi men and women is affected by arsenic metabolism and the age at first exposure. Toxicol Appl Pharmacol. 2008;230:9–16. [PubMed]
20. Valenzuela OL, Borja-Aburto VH, Garcia-Vargas GG, Cruz-Gonzalez MB, Garcia-Montalvo EA, et al. Urinary trivalent methylated arsenic species in a population chronically exposed to inorganic arsenic. Environ Health Perspect. 2005;113:250–254. [PMC free article] [PubMed]
21. Drobna Z, Waters SB, Walton FS, LeCluyse EL, Thomas DJ, et al. Interindividual variation in the metabolism of arsenic in cultured primary human hepatocytes. Toxicol Appl Pharmacol. 2004;201:166–177. [PubMed]
22. Chung JS, Kalman DA, Moore LE, Kosnett MJ, Arroyo AP, et al. 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. [PMC free article] [PubMed]
23. Hernandez A, Marcos R. Genetic variations associated with interindividual sensitivity in the response to arsenic exposure. Pharmacogenomics. 2008;9:1113–1132. [PubMed]
24. Kang HM, Sul JH, Service SK, Zaitlen NA, Kong SY, et al. Variance component model to account for sample structure in genome-wide association studies. Nat Genet. 2010;42:348–354. [PMC free article] [PubMed]
25. Thornton T, McPeek MS. ROADTRIPS: case-control association testing with partially or completely unknown population and pedigree structure. Am J Hum Genet. 2010;86:172–184. [PubMed]
26. Chen Y, van Geen A, Graziano JH, Pfaff A, Madajewicz M, et al. Reduction in urinary arsenic levels in response to arsenic mitigation efforts in Araihazar, Bangladesh. Environ Health Perspect. 2007;115:917–923. [PMC free article] [PubMed]
27. Song X, Geng Z, Li X, Hu X, Bian N, et al. New insights into the mechanism of arsenite methylation with the recombinant human arsenic (+3) methyltransferase (hAS3MT). Biochimie. 2010;92:1397–1406. [PubMed]
28. Engstrom K, Vahter M, Mlakar SJ, Concha G, Nermell B, et al. 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. [PMC free article] [PubMed]
29. 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. [PMC free article] [PubMed]
30. Schlawicke Engstrom K, Nermell B, Concha G, Stromberg U, Vahter M, et al. Arsenic metabolism is influenced by polymorphisms in genes involved in one-carbon metabolism and reduction reactions. Mutat Res. 2009;667:4–14. [PubMed]
31. Chung CJ, Hsueh YM, Bai CH, Huang YK, Huang YL, et al. Polymorphisms in arsenic metabolism genes, urinary arsenic methylation profile and cancer. Cancer Causes Control. 2009;20:1653–1661. [PubMed]
32. Agusa T, Iwata H, Fujihara J, Kunito T, Takeshita H, et al. Genetic polymorphisms in AS3MT and arsenic metabolism in residents of the Red River Delta, Vietnam. Toxicol Appl Pharmacol. 2009;236:131–141. [PubMed]
33. Hwang YH, Chen YH, Su YN, Hsu CC, Yuan TH. Genetic polymorphism of As3MT and delayed urinary DMA excretion after organic arsenic intake from oyster ingestion. J Environ Monit. 2010;12:1247–1254. [PubMed]
34. Sampayo-Reyes A, Hernandez A, El-Yamani N, Lopez-Campos C, Mayet-Machado E, et al. Arsenic induces DNA damage in environmentally exposed Mexican children and adults. Influence of GSTO1 and AS3MT polymorphisms. Toxicol Sci. 2010;117:63–71. [PubMed]
35. 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. [PMC free article] [PubMed]
36. Stuiver M, Lainez S, Will C, Terryn S, Gunzel D, et al. CNNM2, encoding a basolateral protein required for renal Mg2+ handling, is mutated in dominant hypomagnesemia. Am J Hum Genet. 2011;88:333–343. [PubMed]
37. Stachowicz M, Hiemstra T, van Riemsdijk WH. Multi-competitive interaction of As(III) and As(V) oxyanions with Ca(2+), Mg(2+), PO(3−)(4), and CO(2−)(3) ions on goethite. J Colloid Interface Sci. 2008;320:400–414. [PubMed]
38. Srivastava D, Subramanian RB, Madamwar D, Flora SJ. Protective effects of selenium, calcium, and magnesium against arsenic-induced oxidative stress in male rats. Arh Hig Rada Toksikol. 2010;61:153–159. [PubMed]
39. Wood TC, Salavagionne OE, Mukherjee B, Wang L, Klumpp AF, et al. Human arsenic methyltransferase (AS3MT) pharmacogenetics: gene resequencing and functional genomics studies. J Biol Chem. 2006;281:7364–7373. [PubMed]
40. Fujihara J, Yasuda T, Kato H, Yuasa I, Panduro A, et al. Genetic variants associated with arsenic metabolism within human arsenic (+3 oxidation state) methyltransferase show wide variation across multiple populations. Arch Toxicol. 2011;85:119–125. [PubMed]
41. Ahsan H, Chen Y, Parvez F, Argos M, Hussain AI, et al. Health Effects of Arsenic Longitudinal Study (HEALS): description of a multidisciplinary epidemiologic investigation. J Expo Sci Environ Epidemiol. 2006;16:191–205. [PubMed]
42. Verret WJ, Chen Y, Ahmed A, Islam T, Parvez F, et al. A randomized, double-blind placebo-controlled trial evaluating the effects of vitamin E and selenium on arsenic-induced skin lesions in Bangladesh. J Occup Environ Med. 2005;47:1026–1035. [PubMed]
43. Ahsan H, Chen Y, Parvez F, Zablotska L, Argos M, et al. Arsenic exposure from drinking water and risk of premalignant skin lesions in Bangladesh: baseline results from the Health Effects of Arsenic Longitudinal Study. Am J Epidemiol. 2006;163:1138–1148. [PubMed]
44. Nixon DE, Mussmann GV, Eckdahl SJ, Moyer TP. Total arsenic in urine: palladium-persulfate vs nickel as a matrix modifier for graphite furnace atomic absorption spectrophotometry. Clin Chem. 1991;37:1575–1579. [PubMed]
45. Nermell B, Lindberg AL, Rahman M, Berglund M, Persson LA, et al. Urinary arsenic concentration adjustment factors and malnutrition. Environ Res. 2008;106:212–218. [PubMed]
46. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–575. [PubMed]
47. Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, et al. Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006;38:904–909. [PubMed]
48. Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, et al. LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics. 2010;26:2336–2337. [PMC free article] [PubMed]
49. Li Y, Willer CJ, Ding J, Scheet P, Abecasis GR. MaCH: using sequence and genotype data to estimate haplotypes and unobserved genotypes. Genet Epidemiol. 2010;34:816–834. [PMC free article] [PubMed]
50. Barrett JC. Haploview: Visualization and analysis of SNP genotype data. Cold Spring Harb Protoc. 2009;2009:pdb ip71. [PubMed]
51. Xu Z, Taylor JA. SNPinfo: integrating GWAS and candidate gene information into functional SNP selection for genetic association studies. Nucleic Acids Res. 2009;37:W600–605. [PMC free article] [PubMed]
52. Hafeman DM. “Proportion explained”: a causal interpretation for standard measures of indirect effect? Am J Epidemiol. 2009;170:1443–1448. [PubMed]

Articles from PLoS Genetics are provided here courtesy of Public Library of Science