Polymorphisms in Arsenic(+III Oxidation State) Methyltransferase (AS3MT) Predict Gene Expression of AS3MT as Well as Arsenic Metabolism
Karin Engström1, Marie Vahter2, Simona Jurkovic Mlakar3, Gabriela Concha4, Barbro Nermell2, Rubhana Raqib5, Alejandro Cardozo6, Karin Broberg1
1 Department of Laboratory Medicine, Section of Occupational and Environmental Medicine, Lund University, Lund, Sweden, 2 Institute of Environmental Medicine, Section for Metals and Health, Karolinska Institutet, Stockholm, Sweden, 3 University Medical Centre Ljubljana, Clinical Institute of Clinical Chemistry and Biochemistry, Ljubljana, Slovenia, 4 National Food Administration, Toxicology Division, Uppsala, Sweden, 5 International Centre for Diarrhoeal Disease Research, Bangladesh, Dhaka, Bangladesh, 6 San Antonio de Los Cobres Hospital, San Antonio de Los Cobres, Argentina
Environ Health Perspect 119:182-188 (2011). http://dx.doi.org/10.1289/ehp.1002471 [online 24 September 2010]
Background: Arsenic (As) occurs as monomethylarsonic acid (MMA) and dimethylarsinic acid (DMA) in humans, and the methylation pattern demonstrates large interindividual differences. The fraction of urinary MMA is a marker for susceptibility to As-related diseases.
Objectives: We evaluated the impact of polymorphisms in five methyltransferase genes on As metabolism in two populations, one in South America and one in Southeast Asia. The methyltransferase genes were arsenic(+III oxidation state) methyltransferase (AS3MT), DNA-methyltransferase 1a and 3b (DNMT1a and DNMT3b, respectively), phosphatidylethanolamine N-methyltransferase (PEMT), and betaine-homocysteine methyltransferase (BHMT). AS3MT expression was analyzed in peripheral blood.
Methods: Subjects were women exposed to As in drinking water in the Argentinean Andes [n = 172; median total urinary As (U-As), 200 µg/L] and in rural Bangladesh (n = 361; U-As, 100 µg/L; all in early pregnancy). Urinary As metabolites were measured by high-pressure liquid chromatography/inductively coupled plasma mass spectrometry. Polymorphisms (n = 22) were genotyped with Sequenom, and AS3MT expression was measured by quantitative real-time polymerase chain reaction using TaqMan expression assays.
Results: Six AS3MT polymorphisms were significantly associated with As metabolite patterns in both populations (p ≤ 0.01). The most frequent AS3MT haplotype in Bangladesh was associated with a higher percentage of MMA (%MMA), and the most frequent haplotype in Argentina was associated with a lower %MMA and a higher percentage of DMA. Four polymorphisms in the DNMT genes were associated with metabolite patterns in Bangladesh. Noncoding AS3MT polymorphisms affected gene expression of AS3MT in peripheral blood, demonstrating that one functional impact of AS3MT polymorphisms may be altered levels of gene expression.
Conclusions: Polymorphisms in AS3MT significantly predicted As metabolism across these two very different populations, suggesting that AS3MT may have an impact on As metabolite patterns in populations worldwide.
Key words: arsenic, AS3MT, BHMT, DNMT1a, DNMT3b, gene expression, methylation, one-carbon metabolism, PEMT, polymorphism
Address correspondence to K. Engström, Department of Laboratory Medicine, Section of Occupational and Environmental Medicine, Lund University, 221 85 Lund, Sweden. Telephone: 46-0-46-173148. Fax: 46-0-46-143702. E-mail: Karin.email@example.com
This research was supported by grants from the Swedish Council for Working Life and Social Research, the Karolinska Institutet, the Kungliga Fysiografiska Sällskapet (Royal Physiographic Society in Lund), and the European Union Sixth Framework Programme for Research and Technology Development (“PHIME” contract FOOD-CT-2006-016253). The present study is part of the Maternal and Infant Nutrition Interventions in Matlab (MINIMat) study, funded by UNICEF; Swedish International Development Cooperation Agency (Sida); U.K. Medical Research Council; Swedish Research Council; Department for International Development; International Centre for Diarrhoeal Disease Research, Bangladesh; Global Health Research Fund, Japan; Child Health and Nutrition Research Initiative; Uppsala University; and U.S. Agency for International Development.
The authors declare they have no actual or potential competing financial interests.
Received 21 May 2010; accepted 24 September 2010; online 24 September 2010.
Millions of people around the world are exposed to high concentrations of arsenic (As) in drinking water, and these exposures are associated with a number of adverse health effects (International Agency for Research on Cancer 2004; National Research Council 2001). There is a marked individual variability in As metabolism, as measured by the distribution of inorganic As (iAs) and the main methylated metabolites monomethylarsonic acid (MMA) and dimethylarsinic acid (DMA) in urine (Vahter 2002). Several studies have shown an association between an increased fraction of MMA in urine, probably reflecting the highly toxic MMAIII in the tissues, and an increased risk of various As-related adverse health effects (Chen et al. 2005; Chung et al. 2009; Lindberg et al. 2008; Tseng 2007). As is metabolized by a series of reduction and methylation reactions. S-Adenosyl methionine (SAM) is the main methyl donor, implying that As methylation is dependent on factors influencing one-carbon metabolism (Vahter 2002). As(+III oxidation state) methyltransferase (AS3MT) methylates As (Lin et al. 2002). We previously reported an unusual metabolism of As in a mainly indigenous population in the Argentinean Andes, with a lower percentage of MMA and a higher percentage of DMA compared with most other populations (Vahter et al. 1995). This is thought to be partly due to their high frequencies of AS3MT alleles associated with a lower percentage of MMA (%MMA) and higher percentage of DMA (%DMA) (Schläwicke Engström et al. 2007, 2009). The importance of AS3MT polymorphisms for As metabolism has also been also emphasized in other studies, both for noncoding variants (Agusa et al. 2009; Meza et al. 2005; Valenzuela et al. 2009) and for the exonic variant Met287Thr (Hernandez et al. 2008; Lindberg et al. 2007b; Valenzuela et al. 2009). However, there are few data on influential genetic variants in Bangladesh, where about 50 million people drink water with elevated As concentrations and where the inhabitants have shown an efficient methylation of As to DMA, despite prevalent malnutrition, which is likely to impair one-carbon metabolism (Li et al. 2008). Functional data for AS3MT polymorphisms are scarce, apart from a few studies on the Met287Thr variant (Drobna et al. 2004; Wood et al. 2006). Polymorphisms in the promoter and intronic regions may affect the expression of the AS3MT gene, which in turn is likely to affect the metabolite pattern.
The methylation of As may partly be an AS3MT-independent process, because methylated As species have been found in As3mt-knockout mice, albeit a lower fraction than found in As3mt wild-type mice (Drobna et al. 2009). This may be explained by an As-methylating capacity of other methyltransferases. Such candidates could be the DNA methyltransferases (DNMT1a, DNMT3b). DNMTs and AS3MT use SAM as a methyl donor, and there are several parallel characteristics of As methylation and DNA methylation (Vahter 2007).
Several different enzymes are involved in the one-carbon metabolism that generates SAM. For example, the phosphatidylethanolamine N-methyltransferase (PEMT) and betaine-homocysteine methyltransferase (BHMT) may influence As metabolism either by regulating the efficiency of the one carbon metabolism or by direct methylation of As.
In this study, we evaluated the impact of polymorphisms in the above-mentioned methyltransferases (AS3MT, DNMT1a, DNMT3b, PEMT, and BHMT) on As metabolism in two different female populations: in the Andean part of northern Argentina and in rural Bangladesh. Both populations are known to have a wide range of exposure to As in drinking water (Concha et al. 2006; Schläwicke Engström et al. 2007; Vahter et al. 2006), and both populations show efficient methylation of As (Li et al. 2008; Vahter et al. 1995). Furthermore, we investigated the association between AS3MT polymorphisms and AS3MT gene expression in peripheral blood (Argentina), and the association between AS3MT gene expression and As metabolite pattern.
Materials and Methods
Study areas and populations. The study participants in Argentina were mainly women from San Antonio de los Cobres (SAC), a village in the northern Argentinean Andes. All subjects in SAC consumed drinking water from the same source; this water contained about 200 µg As/L, with small variations over time (Concha et al. 2006). To achieve a wider range of exposure, we also included in the study women from three nearby villages: Tolar Grande (3.5 µg As/L in drinking water; n = 21), Olacapato (12 µg As/L; n = 9), and Salar de Pocitos (73 µg As/L; n = 5). The study subjects in these three villages were very similar to the study subjects from SAC regarding ethnicity and socioeconomic factors. In total, urine and blood samples were collected from 176 women (not first-degree relatives) in 2008 (172 of which had both metabolite and genotyping data). The samples did not overlap with the women from the same region included in our earlier studies (Schläwicke Engström et al. 2007, 2009).
The Bangladeshi women studied were from a longitudinal mother–child study of As-related health effects in Matlab, a rural area 53 km southeast of Dhaka. This study group is part of a large, population-based, randomized food and micronutrient supplementation trial in pregnancy (Maternal and Infant Nutrition Interventions in Matlab MINIMat). The study group and sampling procedures were described in detail by Li et al. (2008) and Vahter et al. (2006). Urine and blood were sampled in early pregnancy, urine during gestational week 8, and blood during gestational week 14, on average. From the 2,119 women who were enrolled during 2002, 500 women were randomly chosen for the present study. In total, 440 of the 500 women remained in the trial to receive micronutrient supplementation (no supplementation had been received before the sampling of urine), and 408 of those women provided a blood sample. Of these, urine for gestational week 8 was available for 361 individuals.
The Ethical Review Committee of the International Centre for Diarrhoeal Disease Research, Bangladesh; the Health Ministry of Salta (Argentina); and the Ethics Committee of the Karolinska Institutet (Sweden) approved this study. Oral and written informed consent regarding sample collection was obtained from all participants in Argentina and Bangladesh. The results of the present study were evaluated blind for the individuals’ identity, and none of the participants will receive feedback regarding their own results.
Assessment of As exposure. Exposure to iAs was assessed by the total concentration of As metabolites in urine (U-As; sum of iAs, MMA, and DMA). Speciation of As metabolites in urine was performed using high-pressure liquid chromatography (HPLC) hyphenated with hydride generation and inductively coupled plasma mass spectrometry (ICPMS) (Agilent 1100 series system, Agilent 7500ce; Agilent Technologies, Tokyo, Japan, and Waldbronn, Germany), employing adequate quality control (Lindberg et al. 2007a; Schläwicke Engström et al. 2007). As concentrations were adjusted to the mean specific gravity for each population (1.020 g/mL in Argentina and 1.012 g/mL in Bangladesh), measured by a hand refractometer (Atago, Tokyo, Japan).
Genotyping. DNA was isolated from peripheral blood samples from both groups of subjects, and genotyping was performed by Swegene’s DNA Facility at Malmö University Hospital (Malmö, Sweden) using Sequenom technology (San Diego, CA, USA). Genotyped polymorphisms (22 polymorphisms in five methyltransferases) are shown in Supplemental Material, Table 1 (doi:10.1289/ehp.1002471).
Gene expression. We analyzed expression of the AS3MT gene (which covers 11 exons) by quantitative real-time polymerase chain reaction (qPCR) assays in peripheral blood samples from participants in Argentina. The rationale for the gene expression analyses was to evaluate whether the AS3MT polymorphisms affected gene expression, and if their possible effect on gene expression was due to alternative splicing of AS3MT. To evaluate potential alternative splicing, we ran qPCR assays for four different gene segments. If there was an alternative splicing of the gene, we assumed that the polymorphism that was closest to the amplified region was more likely to have an effect on the splicing of the nearby segments than those situated farther from the amplified region. The assays covered segments near the polymorphisms rs3740400 [assay 1: Taqman assay Hs00960536_g1, exons 1–2; polymorphism numbers from the National Center for Biotechnology Information (NCBI 2010b)], rs3740393 (assay 2: Hs00960529_m1, exons 5–6), and rs3740390 (assay 3: Hs00960532_g1, exons 8–9; all from Applied Biosystems, Foster City, CA, USA) and polymorphism rs1046778 [assay 4: 3´ untranslated region, base pairs 37,166–37,264 in AY817668 (NCBI 2010a)], using Applied Biosystems’ Assays-by-Design. We also evaluated whether AS3MT polymorphisms were associated with the expression of cyclin M2 (CNNM2; Hs00929648_m1, exons 1–2), which is located 16,500 base pairs from the AS3MT 3´ end, with rs1046778 as the closest polymorphism.
For each gene expression assay, we tested associations with the genotype of the closest polymorphism in subsets of the total population that included approximately equal numbers of women with each of the three possible genotypes for that polymorphism [n per assay = 27–29; matched on a group level for age, weight, body mass index (BMI), U-As, and %MMA]. We repeated the analysis of assay 4 on a larger number of subjects with samples available, matched for age, weight, BMI, and U-As (n = 55). Peripheral blood was collected in PAX tubes (PreAnalytiX GmbH, Hombrechtikon, Switzerland), and all samples were frozen and stored at –20°C, after no more than 24 hr at room temperature after sampling. RNA was extracted with the PAXgene Blood RNA Kit (PreAnalytiX) and stored at –80°C. RNA concentration and purity were evaluated on a NanoDrop spectrophotometer (Thermo Scientific, Wilmington, DE, USA; for all samples, 260/280 nm mean optical density ratio = 2.3; range, 1.7–3.3), and RNA integrity (RIN) was evaluated on a Bioanalyzer 2100 (Agilent, Santa Clara, CA, USA; for 17 samples, mean RIN = 8.5; range, 7.6–9.1). For reverse transcription, we used the High Capacity RNA-to-cDNA Kit (Applied Biosystems). Preamplification for AS3MT was performed before qPCR using TaqMan PreAmp MasterMix (Applied Biosystems). AS3MT expression was normalized to hypoxanthine phosphoribosyltransferase 1 (HPRT1) and TATA-binding protein (TBP). HPRT1 and TBP were quantified using SYBR Green PCR Master Mix (Applied Biosystems) in nonpreamplified cDNA. The qPCR (20 µL reaction) was run under standard conditions on an ABI7900 instrument (Applied Biosystems) and performed in triplicate with a negative control included for each assay. Samples with an SD < 0.2 for the threshold cycle (CT) value were accepted. CT values of each sample were normalized to the mean CT value of each reference gene (ΔCT = AS3MT CT – reference gene CT). The ΔCT values were converted to a linear form using the term 2–ΔCT. The coefficients of variance for the mean CT values (2–CT) for all samples were 22% for HPRT1 and 32% for TBP, and we found no difference in expression between individuals with high and low U-As. We reran assay 4 for the 55 individuals to evaluate repeatability. The Pearson product-moment correlation coefficient (r) for the two runs was 0.99.
Statistical analyses. The study groups from Argentina and Bangladesh were analyzed separately. Deviations from Hardy-Weinberg equilibrium were tested using chi-square analysis. We performed linkage disequilibrium (LD) analysis using Haploview (Barrett et al. 2005). Haplotypes were inferred using PHASE software (Stephens and Donnelly 2003).
In consideration of normally distributed residuals, the percentage of iAs (%iAs), %MMA, and total U-As were natural log (ln)-transformed for both study populations in all analyses. Associations between genotypes or haplotypes (independent variables) and urinary As metabolites (dependent variables: %iAs, %MMA, and %DMA) were analyzed by one-way analysis of variance. For each metabolite (dependent variables), we used a multivariable regression model (using the general linear model) that included genotype or haplotype (as categorical variables) and ln(U-As) (as a continuous and independent variable) to derive p-values for associations of genotypes or haplotypes with each metabolite. Additionally, for each metabolite, we used a multivariable regression model that included genotype or haplotype, ln(U-As), and adjustments for other variables. For each metabolite, all potentially influential independent variables were tested in univariate analyses, and those with p < 0.2 were included in the latter multivariable regression model. The potentially influential variables tested for the Argentinean study group were body weight, BMI, age, parity (all modeled as continuous variables), and coca use (yes/no), whereas those for the Bangladeshi study group were gestational week, body weight, BMI, age, and parity (all modeled as continuous variables).
Each polymorphism was modeled as a categorical variable (zero, one, or two alleles). When the frequency of a homozygous genotype was very low (one to three individuals), the group was pooled with the heterozygotes. Haplotypes were also modeled as categorical variables (zero, one, or two copies of the haplotype). The genotype or haplotype associated with a low %MMA and a high %DMA was used as the reference genotype, because the aim was to discover which genotypes are associated with high fractions of MMA, a risk factor for As-related health effects. Observed geometric or arithmetic mean percent metabolite values are reported according to genotype for polymorphisms that were significantly associated with one or more metabolites (p < 0.05) in models adjusted for ln(U-As) only or in models adjusted for ln(U-As) and other variables. Observed geometric or arithmetic mean percent metabolite values are also reported for the most common AS3MT haplotypes. We used models for hypothesis testing only, and no model-based estimates are reported.
We evaluated associations of AS3MT genotype with AS3MT gene expression using the nonparametric Kruskal-Wallis test. The association between gene expression in assay 4 and As the metabolite pattern in urine was evaluated using Spearman’s rank correlation coefficient (rs).
All statistical analyses were performed using SPSS (version 17; SPSS Inc., Chicago, IL, USA). We considered p < 0.05 to be statistically significant.
General characteristics. All polymorphisms were in Hardy-Weinberg equilibrium except for DNMT1a rs2228611 in Argentina (χ2 = 4.9, p = 0.03; in Bangladesh, χ2 = 0.4, p 0.52).
Descriptive data for both study groups are presented in Table 1. No smokers were included in the study groups. Compared with women in Bangladesh, women in Argentina had U-As about twice as high (median, 200 µg/L; Bangladesh median, 100), lower %MMA (7.9% vs. 10%), and higher %DMA (80 vs. 75%). The U-As range was wide in both populations (10–1,200 in Argentina; 5–1,200 in Bangladesh).
AS3MT genotype analysis. The seven noncoding AS3MT polymorphisms had completely different allelic frequencies in the two study groups [see Supplemental Material, Table 1 (doi:10.1289/ehp.1002471)]. In Argentina, all AS3MT polymorphisms, except for rs11191439 (Met287Thr), showed strong LD (R2 between 0.82 and 0.98; see Supplemental Material, Figure 1). In Bangladesh, we detected two LD clusters: one cluster with rs3740400, rs10748835, and rs1046778 (cluster 1) and one cluster with rs3740393, rs3740390, and rs11191453 (cluster # 2) (see Supplemental Material, Figure 1). The polymorphism rs11191439 showed no LD with any of the other AS3MT polymorphisms evaluated in either study population.
Box plots depicting relative AS3MT gene expression according to genotype in women from Argentina; expression for the reference genotype (associated with low %MMA and high %DMA) is set to 1. (A) Relative gene expression based on assay 1 (covering exons 1–2) according to rs3740400 genotype. (B) Relative gene expression based on assay 4 (covering part of the 3´ untranslated region) according to rs1046778 genotype. Numbers below the boxes indicate the number of individuals with that genotype. The bottom and top of each box indicate 25th and 75th percentiles, respectively; the line inside the box is the median; whiskers represent the smallest and the largest values that are not outliers. Circles and asterisks indicate outliers; squares indicate outliers outside the range of the y-axis, which had a gene expression of 25 (CT genotype) and 23 (TT genotype).
In Argentina, all AS3MT polymorphisms, except for the nonsynonymous rs11191439, were significantly associated with the metabolite pattern, and for each polymorphism the allele associated with lower %iAs, lower %MMA, and higher %DMA was the most common (Table 2). In Bangladesh, all AS3MT polymorphisms except rs7085104 were significantly associated with the metabolite pattern. However, patterns of associations varied between the two LD clusters: Polymorphisms in cluster 1 were significantly associated with %MMA but not %iAs or %DMA, whereas polymorphisms in cluster 2 were significantly associated with %iAs and %DMA but not %MMA. In most cases, the alleles associated with lower %MMA or lower %iAs and higher %DMA were the less frequent ones in Bangladesh.
In both populations, participants possessing one or two copies of haplotype 1 had higher %MMA than did those who did not carry this haplotype (Table 3). Haplotype 1 was the most common haplotype in Bangladesh (52%). Haplotype 2 was associated with lower %iAs, lower %MMA, and higher %DMA in both populations (although %MMA was not statistically significant in Bangladesh) and was more frequent in Argentina (70%) than in Bangladesh (16.5%). Haplotype 3 was significantly associated with lower %MMA, but this haplotype was present almost exclusively in Bangladesh (16.5%). Haplotype 4 was associated with higher %iAs and lower %DMA in both populations, although this was statistically significant only in Bangladesh. This was the only haplotype containing the C-allele of rs11191439, and we found no individuals with two copies of this haplotype in either population.
AS3MT expression analysis. Heterozygous and homozygous genotypes that included alleles associated with high %MMA and low %DMA among the women in Argentina were associated with higher AS3MT expression than were corresponding reference genotypes associated with low %MMA and high %DMA. Associations were statistically significant for rs3740400 based on assay 1 (observed gene expressions for the AA and CA genotypes were 240% and 220% higher, respectively, than for the CC genotype; p = 0.009) and for rs1046778 based on assay 4 (observed gene expressions for the TT and CT genotypes were 69% and 36%, respectively, higher than for the CC genotype; p = 0.008) (Figure 1); the p-values for the other assays were 0.18 (for rs3740393 based on assay 2) and 0.077 (for rs3740390 based on assay 3). We ran both endogenous control genes for all assays. The results were very similar when using HPRT1 or TBP as endogenous controls, but the HPRT1 gene demonstrated less variability (the coefficient of variance was 22% for HPRT1 vs. 32% for TBP). We found no association between the AS3MT rs1046778 polymorphism and expression of the downstream gene CNNM2 (p = 0.82). Carriers of each genotype demonstrated a similar expression pattern for all AS3MT gene segments analyzed, and we found no evidence that a certain genotype caused alternative splicing. We repeated the gene expression analysis based on assay 4 on a larger set of individuals (n = 55) and observed the same pattern: relative to the rs1046778 CC genotype, AS3MT expression was 175% and 80% higher among those with rs1046778 TT and CT genotypes, respectively (p < 0.001; Figure 1). Among these 55 individuals, expression of AS3MT (assay 4) was positively (but not significantly) associated with %iAs (rs = 0.16, p = 0.24) and %MMA (rs = 0.18, p = 0.18) and inversely associated with %DMA (rs = –0.28, p = 0.090).
DNMT1a and DNMT3b. The minor allele frequencies for the five polymorphisms analyzed in DNMT1a were 4–40% in Argentina and 2–44% in Bangladesh, whereas the minor allele frequencies for the three polymorphisms analyzed in DNMT3b were 5–6% in Argentina and 23–44% in Bangladesh [see Supplemental Material, Table 1 (doi:10.1289/ehp.1002471)]. The LD pattern for the DNMT genes was similar in Bangladesh and Argentina. The rs10854076 and rs2228612 polymorphisms in DNMT1a showed strong LD in both populations (R2 = 0.99 in Argentina, 0.84 in Bangladesh), whereas the R2 values between the other polymorphisms in DNMT1a were rather low. All DNMT3b polymorphisms showed LD in both study populations, with R2 values ranging from 0.39 to 0.79 in Bangladesh and from 0.88 to 1 in Argentina.
We found no significant associations between DNMT polymorphisms and As metabolism among the Argentinean women. However, four DNMT polymorphisms were significantly associated with metabolite patterns among women in Bangladesh (Table 2): The DNMT1a polymorphisms rs2228612 and rs10845076 (in complete LD) predicted %iAs (p = 0.028, highest for the common rs2228612 AA variant), whereas the DNMT3b polymorphisms rs6087990 and rs2424913 predicted %iAs (p = 0.030 and 0.044, highest for the common rs6087990 CC and rs2424913 TT genotypes, respectively) and %DMA (p-values = 0.013 and 0.021, respectively, with lower mean values for the CC and TT genotypes). We observed no significant associations with metabolite patterns in either population for DNMT1a rs2228611, rs16999593, or rs7253062 or for DNMT3b rs2424932 (data not shown).
PEMT and BHMT. For PEMT, rs1531100 and rs4244598 showed LD for both populations (R2 = 0.70 in Bangladesh and 1 in Argentina), but rs2278952 and rs897453 showed no LD with any other PEMT polymorphism evaluated. The two BHMT polymorphisms showed no LD. Polymorphisms in PEMT or BHMT were not statistically significant predictors of As metabolite patterns.
This study demonstrated strong associations between AS3MT polymorphisms and As metabolism in two populations, one from the Argentinean Andes and one from Bangladesh, both of which are heterogeneous in genetic background and in living conditions. The finding of two LD clusters with different patterns of metabolite fractions suggests the presence of several influential polymorphisms along the gene. Noncoding polymorphisms in AS3MT were significantly associated with altered gene expression. In addition, four polymorphisms in DNMT1a and DNMT3b were associated with As metabolites among the pregnant women from Bangladesh.
One strength of this study is that both populations had wide ranges of both total U-As and As metabolites. Also, the study subjects within each population were rather homogeneous regarding social status and food habits. However, the populations differed from each other in a number of ways, not only genetically but also regarding food habits, nutritional status, lifestyle factors, and pregnancy (all women from Bangladesh were in early pregnancy). This strengthens the hypothesis that the impact of AS3MT on As metabolite pattern should be valid for people worldwide and not strongly modified by environmental or biological factors.
We are aware that relatively few individuals were included in the gene expression analyses. Nevertheless, genotype was highly predictive of the gene expression of AS3MT. This was apparent despite the fact that gene expression was measured in peripheral blood, which has a lower expression of AS3MT than in the main As-methylating tissues, that is, liver and kidney (Su et al. 2004). We used blood as a proxy for liver because sampling liver tissues is a difficult and invasive procedure. It is, however, important to note that extrapolation of gene expression from blood to liver is difficult because of tissue-specific differences in transcription factors and other regulating elements, as well as differences in metabolite distributions between blood and liver, that may affect the result.
Six polymorphisms were associated with As metabolite pattern in subjects from both Argentina and Bangladesh. Several AS3MT polymorphisms were associated with %iAs, which we did not detect in our earlier studies in Argentina (Schläwicke Engström et al. 2007, 2009), possibly due to a lower power in those studies (which included about 100 fewer individuals). We confirmed associations between six AS3MT polymorphisms and As metabolism observed in our earlier studies in another group of individuals from the same region in Argentina (Schläwicke Engström et al. 2007, 2009). Additionally, two polymorphisms not investigated in any of our earlier studies (rs11191453 and rs1046778) significantly predicted the metabolite pattern. Because many polymorphisms showed strong LD and because the significance levels are dependent on the genotype frequencies, it is not possible to attribute causal effects to specific variants. However, in Bangladesh we observed associations with %MMA mainly for polymorphisms in cluster 1 (rs3740400, rs10748835, and rs1046778), which are situated in the 5´ and the 3´ end of the gene, whereas we observed associations with %iAs and %DMA mainly for polymorphisms in cluster 2 (rs3740393, rs3740390, and rs11191453), which are situated in the middle and the 3´ end of the gene. We confirmed previous findings of a higher %MMA and lower %DMA in association with the Thr allele in Met287Thr (rs111191439) (Hernandez et al. 2008; Lindberg et al. 2007b; Valenzuela et al. 2009). The Thr allele in Met287Thr has also been associated with a marginally higher risk (p = 0.055) of premalignant skin lesions (Valenzuela et al. 2009). This polymorphism had rather low frequencies in the populations studied, which were from Chile (Hernandez et al. 2008), Eastern Europe (Lindberg et al. 2007b), and Mexico (Valenzuela et al. 2009), where the allele frequencies were around 5–10%.
The fact that the women from Bangladesh were in early pregnancy may have had an impact on the association between polymorphisms in methyltransferases and As metabolite patterns. Pregnancy is known to induce methylation reactions required for tissue growth and development, as well as As methylation (Concha et al. 1998; Szyf 2009), which may have obscured associations with polymorphisms in methyltransferase genes among the women in Bangladesh. However, to determine the effect of pregnancy on the associations between genotype and As metabolite pattern, pregnant and nonpregnant women from the same population should be compared, or a study with repeated sampling during pregnancy could be performed.
The frequencies of the AS3MT polymorphisms associated with As metabolism deviated markedly between the Argentinean and Bangladeshi populations. In the Argentinean population, the alleles and haplotypes associated with more efficient methylation were very frequent, which is interesting considering the fact that the indigenous people have been living in this As-contaminated area for many generations (Figueroa et al. 1992), suggesting a selection for genotypes associated with efficient As metabolism. In Bangladesh, where most of the tube wells have been installed in the last 20–40 years, the alleles associated with more efficient methylation were the least common; nevertheless, allele frequencies were rather high (~ 45%) for the polymorphisms in LD cluster 1. Still, the general evidence of associations between AS3MT genotype on the As metabolite fraction was similar in both populations. The present study, as well as other studies (Agusa et al. 2009; Hernandez et al. 2008; Lindberg et al. 2007b; Meza et al. 2005), strongly support a causal association of AS3MT genotype on As metabolism worldwide.
A higher %MMA is linked to several adverse health effects (Hsueh et al. 1997; Lindberg et al. 2008; Tseng 2007; Yu et al. 2000). Individuals in Bangladesh who had an MMA level > 8% had a significantly higher risk (p < 0.01) for skin lesions compared with the reference group (where %MMA was < 8) (Lindberg et al. 2008). As demonstrated here, AS3MT genotypes significantly predict the fraction of metabolites and may explain part of the individual differences in the fraction of MMA detected and, in turn, explain differences in susceptibility to skin lesions. Thus, AS3MT genotype is probably a susceptibility factor for As toxicity, and the allelic distribution of AS3MT polymorphisms should be taken into account in future risk evaluations of As in drinking water and food.
In the present study, we demonstrate for the first time that noncoding AS3MT polymorphisms were associated with AS3MT gene expression levels. The common homozygotes in the Argentinean population were associated with lower %MMA and higher %DMA and were also associated with lower AS3MT expression relative to other genotypes. We could not determine whether associations were causal effects due to a single AS3MT polymorphism or a specific haplotype, because all polymorphisms, spread throughout the gene, showed strong LD in Argentina. However, polymorphisms in the 5´ part of the gene (assay 1, 5´ of rs3740400) and 3´ untranslated region (assay 4, including rs1046778), which are likely to involve regulatory sequences, were significantly associated with gene expression. In addition, expression was also associated, although not significantly, with metabolite patterns in urine. We saw no indications that the genotype influenced the splicing of the gene, because we observed a similar pattern of gene expression for all four gene segments analyzed.
Wood et al. (2006) evaluated the potential impact of polymorphisms on AS3MT protein levels in vitro, where the transient expression of the wild-type AS3MT sequence and the variant allozyme of Met287Thr were evaluated in COS-1 monkey kidney cells. Endogenous activity of untransfected COS-1 cells was negligible. The 287Thr allozyme, which has been associated with a higher %MMA and lower %DMA, had higher in vitro protein levels than did those of the wild type. In the present study, very few Argentinean subjects possessed the Thr allele, so it was not possible to evaluate whether this polymorphism was associated with AS3MT gene expression. However, for the polymorphisms that we included in the gene expression analyses, the alleles associated with higher %MMA and lower %DMA were associated with higher gene expression, similar to the association between the 287Thr allele, linked to higher %MMA and lower %DMA, and higher in vitro protein levels. Gene expression and protein levels are not always correlated because posttranscriptional regulation processes may be involved. Further studies are needed to evaluate the impact of polymorphisms on the protein level.
We observed associations with DNMT polymorphisms in Bangladesh only, and these results need to be confirmed. DNMT1a rs16999593 was previously associated with lower %DMA in Argentina (Schläwicke Engström et al. 2009), but this association was not evident in the present study population. However, some of the DNMT polymorphisms had a very low allele frequency in Argentina.
Polymorphisms in AS3MT significantly predicted As metabolism in these two very different populations, which suggests that an effect of AS3MT on As metabolite pattern may be present in populations worldwide. AS3MT polymorphisms were also associated with AS3MT expression in peripheral blood, such that polymorphisms associated with efficient methylation were associated with reduced expression of the AS3MT gene. Evidence of an effect of the polymorphisms is of importance for risk assessment of As, because a low methylation efficiency is associated with increased risk for As toxicity.
- Agusa T, Iwata H, Fujihara J, Kunito T, Takeshita H, Minh TB, et al. 2009. Genetic polymorphisms in AS3MT and arsenic metabolism in residents of the Red River Delta, Vietnam. Toxicol Appl Pharmacol 236:131–141.
- Barrett JC, Fry B, Maller J, Daly MJ. 2005. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21:263–265.
- Chen CJ, Hsu LI, Wang CH, Shih WL, Hsu YH, Tseng MP, et al. 2005. Biomarkers of exposure, effect, and susceptibility of arsenic-induced health hazards in Taiwan. Toxicol Appl Pharmacol 206:198–206.
- Chung WH, Sung BH, Kim SS, Rhim H, Kuh HJ. 2009. Synergistic interaction between tetra-arsenic oxide and paclitaxel in human cancer cells in vitro. Int J Oncol 34:1669–1679.
- Concha G, Nermell B, Vahter M.. 2006. Spatial and temporal variations in arsenic exposure via drinking-water in northern Argentina. J Health Popul Nutr 24:317–326.
- Concha G, Vogler G, Lezcano D, Nermell B, Vahter M.. 1998. Exposure to inorganic arsenic metabolites during early human development. Toxicol Sci 44:185–190.
- Drobna Z, Naranmandura H, Kubachka KM, Edwards BC, Herbin-Davis K, Styblo M, et al. 2009. Disruption of the arsenic (+3 oxidation state) methyltransferase gene in the mouse alters the phenotype for methylation of arsenic and affects distribution and retention of orally administered arsenate. Chem Res Toxicol 22:1713–1720.
- Drobna Z, Waters SB, Walton FS, LeCluyse EL, Thomas DJ, Styblo M. 2004. Interindividual variation in the metabolism of arsenic in cultured primary human hepatocytes. Toxicol Appl Pharmacol 201:166–177.
- Figueroa LT, Razmilic BB, González MU. 1992. Corporal distribution of arsenic in mummied bodies owned to an arsenical habitat. In: International Seminar Proceedings: Arsenic in the Environment and Its Incidence on Health (Sancha AM, ed). Santiago, Chile:Universidad de Chile, Facultad de Ciencias Fisicas y Matematicas. pp. 77–82.
- Hernandez A, Xamena N, Surralles J, Sekaran C, Tokunaga H, Quinteros D, et al. 2008. Role of the Met(287)Thr polymorphism in the AS3MT gene on the metabolic arsenic profile. Mutat Res 637:80–92.
- Hsueh YM, Chiou HY, Huang YL, Wu WL, Huang CC, Yang MH, et al. 1997. Serum beta-carotene level, arsenic methylation capability, and incidence of skin cancer. Cancer Epidemiol Biomarkers Prev 6:589–596.
- International Agency for Research on Cancer 2004. Some drinking-water disinfectants and contaminants, including arsenic. IARC Monog Eval Carcinog Risk Hum 84:1–477.
- Li L, Ekstrom EC, Goessler W, Lonnerdal B, Nermell B, Yunus M, et al. 2008. Nutritional status has marginal influence on the metabolism of inorganic arsenic in pregnant Bangladeshi women. Environ Health Perspect 116:315–321.
- Lin S, Shi Q, Nix FB, Styblo M, Beck MA, Herbin-Davis KM, et al. 2002. A novel S-adenosyl-L-methionine:arsenic(III) methyltransferase from rat liver cytosol. J Biol Chem 277:10795–10803.
- Lindberg AL, Goessler W, Grander M, Nermell B, Vahter M. 2007. . Evaluation of the three most commonly used analytical methods for determination of inorganic arsenic and its metabolites in urine. Toxicol Lett 168:310–318.
- Lindberg AL, Kumar R, Goessler W, Thirumaran R, Gurzau E, Koppova K, et al. 2007. . Metabolism of low-dose inorganic arsenic in a central European population: influence of sex and genetic polymorphisms. Environ Health Perspect 115:1081–1086.
- Lindberg AL, Rahman M, Persson LA, Vahter M. 2008. 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 230:9–16.
- Meza MM, Yu L, Rodriguez YY, Guild M, Thompson D, Gandolfi AJ, et al. 2005. Developmentally restricted genetic determinants of human arsenic metabolism: association between urinary methylated arsenic and CYT19 polymorphisms in children. Environ Health Perspect 113:775–781.
- NCBI (National Center for Biotechnology Information) 2010a. Nucleotide Database. Available: http://www.ncbi.nlm.nih.gov/sites/entrez?db=nucleotide [accessed 11 May 2010]
- NCBI (National Center for Biotechnology Information) 2010b. Single Nucleotide Polymorphism Database. Available: http://www.ncbi.nlm.nih.gov/snp/ [accessed 6 January 2011]
- National Research Council 2001. Arsenic in Drinking Water: 2001 Update. Washington, DC:National Academy Press.
- Schläwicke Engström K, Broberg K, Concha G, Nermell B, Warholm M, Vahter M.. 2007. Genetic polymorphisms influencing arsenic metabolism: evidence from Argentina. Environ Health Perspect 115:599–605.
- Schläwicke Engström K, Nermell B, Concha G, Stromberg U, Vahter M, Broberg K.. 2009. Arsenic metabolism is influenced by polymorphisms in genes involved in one-carbon metabolism and reduction reactions. Mutat Res 667:4–14.
- Stephens M, Donnelly P.. 2003. A comparison of Bayesian methods for haplotype reconstruction from population genotype data. Am J Hum Genet 73:1162–1169.
- Su AI, Wiltshire T, Batalov S, Lapp H, Ching KA, Block D, et al. 2004. A gene atlas of the mouse and human protein-encoding transcriptomes. Proc Natl Acad Sci USA 101:6062–6067.
- Szyf M.. 2009. The early life environment and the epigenome. Biochim Biophys Acta 1790:878–885.
- Tseng CH. 2007. Arsenic methylation, urinary arsenic metabolites and human diseases: current perspective. J Environ Sci Health C Environ Carcinog Ecotoxicol Rev 25:1–22.
- Vahter M.. 2002. Mechanisms of arsenic biotransformation. Toxicology 181–182:211–217.
- Vahter ME. 2007. Interactions between arsenic-induced toxicity and nutrition in early life. J Nutr 137:2798–2804.
- Vahter M, Concha G, Nermell B, Nilsson R, Dulout F, Natarajan AT. 1995. A unique metabolism of inorganic arsenic in native Andean women. Eur J Pharmacol 293:455–462.
- Vahter ME, Li L, Nermell B, Rahman A, El Arifeen S, Rahman M, et al. 2006. Arsenic exposure in pregnancy: a population-based study in Matlab, Bangladesh. J Health Popul Nutr 24:236–245.
- Valenzuela OL, Drobná Z, Hernández-Castellanos E, Sánchez-Peña LC, Garcia-Vargas GG, Borja-Aburto VH, et al. 2009. Association of AS3MT polymorphisms and the risk of premalignant arsenic skin lesions. Toxicol Appl Pharmacol 239(2):200–207.
- Wood TC, Salavagionne OE, Mukherjee B, Wang L, Klumpp AF, Thomae BA, et al. 2006. Human arsenic methyltransferase (AS3MT) pharmacogenetics: gene resequencing and functional genomics studies. J Biol Chem 281:7364–7373.
- Yu RC, Hsu KH, Chen CJ, Froines JR. 2000. Arsenic methylation capacity and skin cancer. Cancer Epidemiol Biomarkers Prev 9:1259–1262.
Come visit with EHP’s Science Editor, Jane Schroeder, at the 2016 ISEE Meeting in Rome, Italy, from 1–4 September 2016. This is a great opportunity to meet a member of our team, and to learn more about the journal. A number of our Associate Editors will be in attendance as well. Jane is also co-hosting a half-day writing and publishing workshop on 4 September, following the conference.
EHP is pleased to present the abstracts from the 28th annual meeting of the International Society for Environmental Epidemiology (ISEE), held in Rome, Italy, 1–4 September 2016, and hosted by the Department of Epidemiology Lazio Regional Health Service, ASL Roma 1, and the Italian Epidemiological Association. The focus of this year’s conference is current and future challenges in exposure assessment, study design, and data analyses.
Featured Children’s Health
Birgit Claus Henn, Adrienne S. Ettinger, Marianne R. Hopkins, Rebecca Jim, Chitra Amarasiriwardena, David C. Christiani, Brent A. Coull, David C. Bellinger, and Robert O. Wright
Diane Gilbert-Diamond, Jennifer A. Emond, Emily R. Baker, Susan A. Korrick, and Margaret R. Karagas
Benjamin B. Green, Margaret R. Karagas, Tracy Punshon, Brian P. Jackson, David J. Robbins, E. Andres Houseman, and Carmen J. Marsit