Introduction
The environment of the developing child is an important determinant of disease susceptibility in later life.Reference Gluckman 1 This phenomenon may be mediated through epigenetic changes, such as alterations of DNA methylation, which alter developmental programming.Reference Attig, Gabory and Junien 2 – Reference Waterland and Michels 5 Lately, it has become increasingly clear that toxic chemicals can also affect DNA methylation.Reference Baccarelli and Bollati 6 – Reference Perera and Herbstman 8
Arsenic, which frequently occurs at elevated concentrations in drinking water worldwide, is a potent carcinogen, a general toxicant and an endocrine disrupter. 9 , 10 Besides the accumulating evidence for associations between in utero arsenic exposure and impaired fetal and child health and development,Reference Hamadani, Tofail and Nermell 11 – Reference Saha, Engstrom and Hamadani 14 both human and experimental animal studies indicate higher risk for cancer, cardiovascular effects and increased mortality later in childhood or adulthood when exposure starts prenatally or shortly after birth.Reference Liaw, Marshall and Yuan 15 – Reference Yorifuji, Tsuda and Grandjean 20 These results suggest persistent long-term effects of early-life changes, possibly through alterations of DNA methylation at the 5-methylcytosine position. Studies on epigenetic effects of arsenic mainly have concerned changes in DNA methylation in relation to cancer development, which is closely related to DNA methylation instability.Reference Bhattacharjee, Banerjee and Giri 21 – Reference Jones and Baylin 25 In particular, arsenic-related changes in gene methylation status have been proposed to silence tumor-suppressor genes, possibly eading to long-term changes in the activity of genes controlling cell transformation.Reference Bhattacharjee, Banerjee and Giri 21 These results suggest that arsenic has epigenetic effects related to its carcinogenicity. Genome-wide DNA methylation has been assessed in a small group (n=16) of adults.Reference Smeester, Rager and Bailey 26 Studies on the effects of arsenic exposure on DNA methylation in newborns are few and have generally examined non-gene-related markers of DNA methylation.Reference Intarasunanont, Navasumrit and Waraprasit 27 – Reference Pilsner, Hall and Liu 29 A recent study by Koestler et al.Reference Koestler, Avissar-Whiting, Houseman, Karagas and Marsit 30 investigated the association between low maternal arsenic exposure measured in late pregnancy and genome-wide effects on 5-methylcytosine methylation in cord blood. Therefore, the aim of the present study was to evaluate cord blood DNA methylation of gene-related CpG sites, covering a large part of the whole epigenome, in relation to a large range of maternal arsenic exposure in both early and late pregnancy. Considering the increasing evidence for sexual dimorphism both for programming trajectories and in response to the environmental insults,Reference Gabory, Attig and Junien 31 including arsenic,Reference Pilsner, Hall and Liu 29 , Reference Nohara, Baba and Murai 32 we evaluated the results in female and male newborns separately.
Methods
Study area and design
The study area, Matlab, is a rural area southeast of Dhaka, where the International Centre for Diarrhoeal Disease Research, Bangladesh (ICDDR,B), runs a health and demographic surveillance system. In a population-based screening in 2002–2003, we found a total range of <1 to 3640 µg/l, with 70% of the tube wells in Matlab exceeding the WHO-recommended limit of arsenic in drinking water of 10 µg/l.Reference Rahman, Vahter and Wahed 33 In an effort to decrease adverse health effect because of polluted drinking water, most families (>95%) have had a well drilled in their yard past few decades ago, or obtain water from community-owned wells.Reference Winston, Escamilla and Perez-Heydrich 34 Unfortunately, a considerable fraction of those wells were later found to contain arsenic, because of the hydrogeological conditions, and this has been related to impaired fetal growth and increased infant morbidity and mortality. Other toxic exposures are low. This is a very pristine rural area of Bangladesh, with very little environmental pollution, for example, essentially no cars.
The study was nested in a large randomized population-based food and micronutrient supplementation trial (MINIMat) in pregnancy.Reference Persson, Arifeen and Ekström 35 Pregnancy was identified by urine test, usually in gestational week (GW) 8. Women who tested positive were advised to visit the health facility for confirmation of pregnancy by ultrasound. The current study is based on 127 women who delivered at any of the health centers before early afternoon (5:00 am–2:30 pm, owing to logistics in transferring samples to the laboratory) from May 2003 to June 2004. According to the questionnaire, all women were non-smokers, while 73 (58%) of the women chewed betel leaves (6 out of 73 women reported betel chewing with tobacco) during pregnancy.
Birth weight was measured mostly within 24 h of delivery (all within 72 h), using electronic scales. Socioeconomic status was based on family assets using data from the health demographic surveillance system database.Reference Saha, Frongillo and Alam 36 Participants gave written informed consent and the study was approved by the Research Review Committee and Ethical Review Committee of ICDDR,B and the Regional Ethics Committee at the Karolinska Institutet, Sweden.
Arsenic exposure
Arsenic exposure was based on the sum concentration of arsenic metabolites [inorganic arsenic, methylarsonic acid (MMA) and dimethylarsinic acid] in spot urine samples collected, as previously described.Reference Vahter, Li and Nermell 37 This is an established biomarker that reflects exposure to inorganic arsenic from all sources. The half-life of arsenic in the body is 3–4 days, which is why spot urine reflects the ongoing exposure. However, because the exposure in the area is mainly from drinking water, with some contribution from rice, the main staple food, the exposure is fairly stable over time.Reference Gardner, Nermell and Kippler 38 Because of the local custom that young women, especially first-time mothers, move to their parents’ home for delivery, we evaluated arsenic exposure in both early and late pregnancy, using urine samples from 127 women that were collected between GW 5 and 14 (median GW 8) and between GW 26 and 36 (median GW 30), respectively. The relative amounts of the different metabolites in urine were used as markers of methylation efficiency.Reference Vahter 39 The reason is that we recently found that the efficiency of arsenic methylation increases very early in pregnancy.Reference Gardner, Nermell and Kippler 38 In particular, the percentage of MMA, which has been associated with increased risk of toxic effects in adults,Reference Li, Engström and Vahter 40 , Reference Pierce, Kibriya and Tong 41 decreased very early in pregnancy. The arsenic metabolite concentrations were measured using hydride generation atomic absorption spectrophotometry with a detection limit of 1.3±0.27 µg/l. All arsenic concentrations were adjusted for variation in urine dilution by specific gravity to the average value, 1.012 g/ml.
DNA isolation and epigenetic analysis
Cord blood mononuclear cells were separated by Ficoll (Pharmacia-Upjohn) density gradient centrifugation. The epigenetic analysis has been described in detail elsewhere.Reference Engström, Hossain and Lauss 42 In brief, 1 µg DNA (50 ng/µl) was bisulfite-treated (Zymo Research). The DNA samples were randomized for gender and arsenic exposure on two 96-well plates for analysis with the Infinium HumanMethylation450K BeadChip (lllumina). All chips were from the same batch. We also included four controls in duplicate, with the duplicate positioned on different 96-well plates and different BeadChips. Three of the control samples were DNA extracted from blood and one was demethylated DNA (Zymo Research). The SCIBLU facility (Lund, Sweden) used 200 ng bisulphite-treated DNA for hybridization to the BeadChip. Unsupervised hierarchial cluster analysis showed that all duplicates of the control samples clustered together.
For all samples, the number of detected CpGs (P-value=0.01) were at least >98%. All 1500 probes (the top 500 selected for all children, top 500 selected for boys and top 500 selected for girls) performed very well with a detection P-value of at least 0.01 (the absolute majority P<0.0001) among 80% of the samples. In an ongoing study on DNA samples from whole blood we are evaluating the concordance between results from the HumanMethylation450K BeadChip and pyrosequencing. We found a correlation of 0.85–0.92 (P>0.0001) for four different CpG sites tested for both methods.
Measurement of LINE1 methylation, as a marker for global methylation status, was performed on bisulphite-treated cord blood DNA from a subgroup of 80 newborns (median of maternal arsenic in urine at GW 8 was 66 µg/l, range 3–740 µg/l): 40 boys and 40 girls. Measurement of four CpG sites in the LINE1 sequence was analyzed by the use of pyrosequencing, as described previously.Reference Hossain, Vahter, Concha and Broberg 43
Statistical analysis
Details of the statistical analyses were reported previously.Reference Kippler, Engström and Mlakar 7 Missing beta values (0.1%) were imputed using k-nearest neighbor imputation (k=10). Principal component analysis (PCA), using the R package swamp, was used to capture the major directions of variation in the data and influencing factors. For each of the top principal components, we fitted a univariate linear model with each of the sample annotations as regressors. The log10 P-values from the models’ F-statistics were plotted as a heat map. The results showed that the analysis plate (two 96-well plates) was associated with methylation levels in the first (P=10−6) and third components (P<10−10). We removed the plate influence at each CpG site by using the residuals from the initial linear regression model of methylation with the analysis plate as regressor. The residuals of the linear model added to the total mean before correction became the new data for each CpG site. In this way, the methylation levels of the sites were not influenced by the analysis plate.
We first evaluated whether the maternal arsenic exposure (urinary arsenic) overall was associated with genome-wide DNA methylation by performing 482,421 separate linear regression models, one for each CpG site, with methylation as the dependent variable and arsenic concentration as the independent variable. The range of arsenic was rescaled to vary between 0 and 1 to minimize the potential effects of skewed data. We tested whether the slope was statistically significant in all 482,421 models. Had there been no effect of arsenic on any CpG site, then the slopes would result from mere sampling error. In this case, the P-values would be distributed uniformly over the 0–1 range.
The associations of CpG site methylation fractions with urinary arsenic were then determined using Spearman correlation. Resulting P-values were corrected for multiple testing by the Benjamini–Hochberg method to obtain false discovery rates.
Some of the associations between DNA methylation and maternal urinary arsenic were non-linear, but after log2 transformation of the arsenic concentrations the associations showed linear patterns, and linear regression analyses were applied. We did not observe any typical genotype-related clustering of DNA methylation when visually inspecting the associations between arsenic and DNA methylation. In multivariable-adjusted analysis, we controlled for the maternal characteristics age, body mass index (BMI) (at GW 8), SES, as well as GW at urine collection, gestational age at delivery and fetal sex (except for analyses stratified for sex). In sensitivity analysis, we additionally adjusted for birth weight, betel chewing including with tobacco and cadmium.
For the pathway analysis, we used the Ingenuity Pathway Analysis Tool (Ingenuity H Systems). We selected the top 500 gene-annotated CpG sites associated with arsenic in maternal urine in boys and girls separately (from Illumina). Differentially methylated genes are mapped to genetic networks/diseases available in the Ingenuity Pathway Knowledge Base and then ranked by score.
Results
Background data
The studied women were on average 25.3 years old and with a BMI of 20.4 kg/m2 (Table 1). The gestational age at birth was 39 weeks on average: 12 births (six boys and six girls) being premature (range 34.4–36.7 weeks at birth). Boys had slightly higher birth weight (2838 g) than girls (2724 g). The median maternal urinary arsenic concentrations were 66 μg/l in GW 8 and 89 μg/l in GW 30, with wide variations. Maternal urinary arsenic concentration did not differ by sex of the child. Women in the study sample were to a larger extent nulliparous, had higher SES and had children with somewhat higher birth weight, compared with the total cohort of 1729 pregnant women studied for arsenic exposure.
Table 1 Characteristics of the 127 mother–child pairs and all other women enrolled in the MINIMat trialFootnote a
BMI, body mass index; GW, gestational week; MMA, methylarsonic acid; DMA, dimethylarsinic acid.
a Enrollment in the MINIMat trial from October 2002 to October 2003.
b Adjusted to average specific gravity of 1.012 g/ml.
c In MINIMat, urinary arsenic was measured at GW 8 (n=1729) and GW 30 (n=1031).
d Percent of total metabolite concentration in urine. Urinary arsenic refers to inorganic arsenic and its methylated metabolites in urine.
Maternal arsenic exposure and genome-wide DNA methylation
The PCA showed that maternal arsenic exposure in early gestation (GW 8) was associated with cord blood DNA methylation in the eighth component and, to a lesser extent, in the fourth component (Supplementary Figure S1). Arsenic exposure measured in late gestation (GW 30) showed weaker association with DNA methylation than that in GW 8.
We first evaluated whether maternal arsenic exposure in early gestation was associated with newborn genome-wide DNA methylation by analyzing all CpG sites in cord blood v. arsenic concentrations in maternal urine in GW 8 (482,421 separate models). Low P-values were more frequent than expected from a uniform distribution (Kolmogorov–Smirnov test, P-value<10–15, Supplementary Figure S2), and there were marked differences between sexes; significant effects of arsenic were considerably more pronounced in boys compared with girls (Fig. 1). The proportion of small P-values was greater than expected in boys (P-value<0.0001), but not in girls (P-value=0.99).
Fig. 1 Distribution of the P-values for the coefficient associated with maternal urinary arsenic concentrations in early gestation from linear regression analysis of CpG methylation (n=482,421 sites) in cord blood. (a) Cord blood from girls, (b) Cord blood from boys.
We selected the top 500 CpG loci with the lowest P-values based on their association with urinary arsenic levels in GW 8 and found hypomethylation with increasing arsenic exposure for 277 (55%) sites, and hypermethylation for 223 (45%) sites (all infants, test if more hypomethylation than expected gave P=0.018). There were major sex differences: among boys 372 (74%) of the top 500 sites showed hypomethylation with increasing arsenic exposure, compared with only 207 (41%) among girls (difference between boys and girls P-value<0.001, Fisher’s exact test). We found no single overlap between the top 500 arsenic-related sites in boys and those in girls. The 20 CpG sites with the strongest association with urinary arsenic in GW 8 in boys and girls are listed in Table 2. Generally, stronger associations were found for boys compared with girls, as estimated from correlation coefficients and P-values; three CpG sites in boys (cg15255455, cg13659051 and cg17646418), but none in girls, were significantly associated with arsenic exposure after adjustment for multiple comparisons (Table 2). There were stronger associations between DNA methylation and maternal arsenic exposure measured in early than in late gestation, but the directions were the same. Examples of the associations are given in Fig. 2, which indicate that the effects of arsenic on DNA methylation start at low exposure levels, well below 100 µg/l.
Figure 2 Scatterplot of the fractions (0 to 1) of methylation of the CpG sites in cord blood and the concentrations of arsenic metabolites in maternal urine in early pregnancy (gestational week 8) for: (a) cg13659051 and (c) cg02975107, among top five CpG sites in boys; and (b) cg06411879 and (d) cg09606015, among top five CpG sites in girls. The associations are indicated by Lowess lines.
Table 2 Top 20 CpG sites with the strongest correlations (r S ) with maternal urinary arsenic concentrations in gestational week (GW) 8 by child sex, as well as regression analysis (GW 8)Footnote a
Chr, chromosome number; FDR, false discovery rate; CI, confidence intervals
In the regression analysis, the arsenic exposure variable was log2-transformed.
a Correlations in GW 30 are also shown.
b P-value adjusted for FDR.
c Adjusted for mother age, body mass index, gestational age at birth, socioeconomic status, exact gestational weeks at urine collection.
For the CpG sites in Table 2, only 1 (cg22614624) of the 40 sites contained an single-nucleotide polymorphism (SNP) within ⩽10 bp from the query site (Supplementary Table S2). This SNP has a minor allele frequency of 0.349, based on very little data and none from Asian populations (www.ncbi.nlm.nih.gov/projects/SNP). We also analyzed the data without the sex chromosomes, but this did not influence the results. The reported P-values were still significant as well as the order of the top SNPs reported.
In the subsequent linear regression analyses (log2-transformed urinary arsenic concentrations in GW 8), the associations between CpG methylation and arsenic exposure in early gestation were very robust when taking other potentially influential covariates into account (Table 2). Apart from GW of urine collection, few covariates were even close to significant in the models and their inclusion did not markedly change the estimates for arsenic. In a sensitivity analysis, we additionally adjusted for birth weight, betel chewing and maternal cadmium in blood, but this only marginally changed the effect estimates. Among the 20 top sites in boys, the effect of arsenic on the degree of CpG methylation ranged between 0.3 and 2.7% for a doubling of arsenic concentrations in maternal urine (µg/l) in early pregnancy. Arsenic caused in boys mainly lower methylation (all 12 sites with >1% change involved hypomethylation). The corresponding range for girls was 0.3–5.4% (of the seven sites with >1% change, two showed lower and five higher methylation).
We performed analysis of LINE1 in a subgroup of 40 boys and 40 girls. There were no significant differences in LINE1 methylation (expressed as mean percentage of methylation of the four CpG sites) between the sexes (mean 68.0% in boys and 68.3% in girls, P=0.75, ANOVA). Furthermore, there were no significant associations between urinary arsenic in GW 8 and LINE1 methylation in all newborns (P=0.86, adjusted for the same covariates as in the analysis above for individual CpG sites P=0.59, ANCOVA), or in boys or girls separately (boys: P=0.73, adjusted P=0.59; girls: P=0.46, adjusted P=0.45).
We additionally stratified the regression analyses by the percentage of MMA in urine (GW 8), as a measure of maternal arsenic methylation efficiency, mothers with relatively high percentage of MMA were considered to have poorer methylation efficiency (Supplementary Table S3). For several of the top 20 sites, the arsenic-related effect on CpG methylation was markedly stronger in newborns to mothers with poor arsenic methylation.
Arsenic exposure in early gestation and epigenetic changes related to disease pathways
We performed pathway analysis of the top 500 gene-annotated CpG sites for which methylation was affected by arsenic in boys and girls. In both boys and girls, there was overrepresentation of networks related to embryonic, organic and cellular development (Table 3). In contrast, cancer-related networks were found to be overrepresented in boys (146 genes, Supplementary Table S4), but not in girls. Of the 146 genes, 67% were hypomethylated with increasing arsenic exposure, the remaining were hypermethylated. When searching for related diseases that were overrepresented in boys, cancer was the first one listed (P-value=4.4×10–5) followed by reproductive system diseases, hematological diseases, hereditary disorder, and skeletal and muscular disorders. In girls, inflammatory disease ranked the highest (P-value=2.7×10–5), followed by neurological diseases, skeletal and muscular disorders, endocrine system disorders and gastrointestinal diseases.
Table 3 Pathway analysisFootnote a , stratified for sex, between maternal urinary arsenic concentrations in early gestation CpG methylation in cord blood
a The input data to the Ingenuity Pathway Analysis software was the top annotated genes from the top 500 scoring CpG sites in cord blood correlated to arsenic in maternal blood in boys and girls, respectively.
We performed linear regression analysis of the top 10 CpG sites identified among the 146 genes associated with cancer in boys (and in girls for comparison). The results (Table 4) show that the associations in boys were robust in the models, and that methylation in these sites was not associated with arsenic in girls. Three CpG sites (cg13659051, cg24042517 and cg20415517) associated with cancer-related genes were also among the top sites in Table 2: LRRC26, HOXB9 and BRSK2.
Table 4 Top 10 cancer-related CpG sites in boys (and girls for comparison) with the strongest correlations (r S ) with maternal arsenic urinary concentrations in gestational week 8
Chr, chromosome number; FDR, false discovery rate; CI, confidence intervals.
In the regression analysis arsenic exposure variable was log2-transformed.
a P-value adjusted for FDR.
b Adjusted for mother’s age, body mass index, gestational age at birth, socioeconomic status, exact gestational weeks at urine collection.
Discussion
This first study of sex-differential genome-wide cord blood DNA methylation in relation to arsenic exposure provides strong evidence that maternal arsenic exposure in early gestation is associated with modified DNA methylation in the newborn child. For many of the CpG sites, arsenic exposure changed the DNA methylation by several percent, and the changes started at low exposure levels, below 50 µg/l in the urine. At the genome-wide DNA methylation level, the associations with arsenic exposure were much more evident in boys than in girls. In particular, arsenic seemed to induce hypomethylation in boys; that is, in 74 of the top 500 CpG sites methylation decreased with increasing arsenic exposure. In a subgroup of the newborns, we evaluated whether the decrease in methylation with increasing arsenic exposure in boys, reflects a general genomic feature, such as interference of arsenic with the methyl donor pool. However, we did not find any overall effect of arsenic on global methylation in either of the sexes. Furthermore, we did not observe multiple probes affected for individual genes but rather strong effects on individual CpG sites in relation to arsenic. This observation, together with the data for global methylation, suggests that arsenic affects accessible areas of the chromatin with, certain structural features, or impacts on specific transcription factors or other chromatin-associated proteins. The pathway analysis showed enrichment of DNA methylation changes in cancer-related genes in boys, but not in girls. This is interesting considering the reported sex differences in cancer development in mice following prenatal arsenic exposure, with later-life occurrence of liver carcinoma and adrenal cortical adenoma in male, but not female, offspring.Reference Waalkes, Liu and Diwan 44 A further observation is noteworthy in relation to our data that the follow-up of the male liver DNA showed significant reduction in methylation of GC-rich regions.Reference Xie, Liu and Benbrahim-Tallaa 19 In contrast, mice that were continuously exposed also after birth showed much less sex differences in cancer development,Reference Tokar, Diwan, Ward, Delker and Waalkes 45 which is in agreement with the observed marked increase in liver cancer mortality in both boys and girls exposed to arsenic through drinking water in northern Chile as young children, but not in those exposed before birth.Reference Liaw, Marshall and Yuan 15 It should be noted that we have identified arsenic-related impaired fetal growth in boys,Reference Kippler, Wagatsuma and Rahman 46 although postnatal growth seems to be more affected in girls.Reference Saha, Engstrom and Hamadani 14 Here we identified different processes involved in development, as the top pathways for arsenic-related DNA methylation change in both boys and girls. One can speculate that arsenic might interfere with already existing sex-specific gene expression early in lifeReference Gabory, Attig and Junien 31 as seen in mice exposed to arsenic,Reference Waalkes, Liu and Chen 47 or alternatively it causes sex-specific DNA methylation for genes that should not differ between the sexes. Thus, it is essential to follow-up on the epigenetic effects of arsenic as the children grow.
The much stronger associations with arsenic exposure in early compared with late gestation may reflect interference with the de novo DNA methylation patterns established in early fetal development.Reference Waterland and Michels 5 The rapid upregulation of the second step in arsenic methylation, resulting in a decrease in the percentage of the most toxic metabolite MMA very early in pregnancy,Reference Gardner, Nermell and Kippler 38 may further contribute to the observed weaker epigenetic associations with late gestation exposure levels. Indeed, our data also suggest that a low capacity for methylation of arsenic, resulting in a higher fraction of urinary MMA, a known susceptibility factor for arsenic toxicity,Reference Li, Engström and Vahter 40 , Reference Pierce, Kibriya and Tong 41 was associated with stronger epigenetic effects. However, the results in relation to methylation capacity need to be followed up, as they were based on a small number of children in each comparison group. Koestler et al.Reference Koestler, Avissar-Whiting, Houseman, Karagas and Marsit 30 also found the effects of prenatal exposure to arsenic on cord blood methylation; however, in contrast to our data, they found increasing methylation with increasing arsenic exposure for most top sites. The results are difficult to compare with ours as the women in their study had markedly lower exposure (median maternal urinary arsenic was 4 v. 66 µg/l in our study) measured late in pregnancy, and moreover they did not stratify their results for sex.
The top sites associated with arsenic might represent possible new candidates for arsenic epigenome toxicity. In boys, three CpG sites associated with different genes (PLIN5, LRRC26 and RPS6KA), not previously linked to arsenic, were significantly hypomethylated in relation to arsenic exposure. Perilipin 5 (PLIN5) is a lipid metabolism-related protein; overexpression of PLIN5 in rat skeletal muscle promoted oxidative gene expression and intramyocellular lipid content.Reference Bosma, Sparks and Hooiveld 48 Heart-specific overexpression of PLIN5 in mice was associated with increased triglyceride concentrations and strong increased expression of oxidative-induced genes via NF-E2-related factor 2 antioxidative pathway.Reference Wang, Sreenivasan and Gong 49 LRRC26 is a protein associated with calcium-activated potassium channels, so-called BK channels. It is normally expressed not only in the fetal brain and thymus, but also in prostate cancer where its expression actually reduced tumor growth and metastasis.Reference Liu, Xiang and Zhang 50 RPS6KA is a type of serine/threonine kinase, implicated in controlling cell growth and differentiation and it has been identified as a risk marker for colon and rectal cancer.Reference Slattery, Lundgreen, Herrick and Wolff 51 In addition, HOXB9 and BRKS2 are worth mentioning as they were among the cancer-related genes and among the top 20 genes for boys. HOXB9, a class I homebox gene, is overexpressed in breast cancer and it induces an epithelial-to-mesenchymal transition, a key factor in metastasis.Reference Chiba, Comaills and Shiotani 52 BRSK2 is an AMP-activated protein kinase and a multifunctional regulator of cell-cycle progression. Further studies are needed to examine the DNA methylation of these newly identified genes related to arsenic exposure later in life and in relation to cancer risk. Moreover, as we measure changes in DNA methylation in blood, they need to be followed up in relation to potential hematopoetic effects.
We have previously reported effects of maternal cadmium exposure on 5-methylcytosine methylation in the same cohort of newborns.Reference Kippler, Engström and Mlakar 7 However, the cadmium-associated sites were not overlapping with those identified here and adjustment for cadmium did not influence the result. Betel chewing was not influential in this study on arsenic, which might be owing to the fact that very few of the mothers mixed betel leaves with tobacco. Further adjustments for SES or other potential influential factors did not alter the associations; however, we cannot completely exclude that there might be residual confounding. A potential caveat of this study is that we were not able to sort cells in the blood samples during the field studies. We therefore measured DNA methylation in cord blood mononuclear cells, which are a mixture of different cell types with partly different methylation patterns. Thus, we cannot exclude that a potential cell-specific effect of arsenic might blur associations between DNA methylation and arsenic exposure. However, in the recent study by Koestler et al.,Reference Koestler, Avissar-Whiting, Houseman, Karagas and Marsit 30 white blood cell distributions explained only a small proportion of the variability in patterns of cord blood DNA methylation associated with maternal arsenic exposure (3% for total arsenic in urine). Furthermore, some results from the Ingenuity Pathway Analysis (IPA) analysis, although plausible, should be carefully interpreted as a recent paper describes that certain diseases and pathways, for example, cancer and developmental-related pathways, tend to be overrepresented in the IPA.Reference Harper, Peters and Gamble 53
In conclusion, we show a clear relationship between maternal exposure to arsenic, especially in early pregnancy, and epigenetic effects early in life, particularly in boys. How that relates to the observed health effects of arsenic in early and later life remains to be elucidated. Our results emphasize the need to evaluate epigenetic effects of fetal environmental factors by sex and critical windows of exposure.
Acknowledgments
None.
Financial Support
This research was supported by the Swedish Council for Working Life and Social Research (K.B., co-PI), the Swedish Research Council (M.V., PI), the Swedish International Development Cooperation Agency (SIDA; M.V., PI) and the Karolinska Institutet (M.V.).
Conflicts of Interest
None.
Ethical Standards
The authors assert that all procedures contributing to this work comply with the ethical standards of the Helsinki Declaration of 1975, as revised in 2008, and has been approved by the institutional committees Research Review Committee and Ethical Review Committee of ICDDR,B in Bangladesh and the Regional Ethics Committee at the Karolinska Institutet in Sweden.
Open access
