Skip to main content

Transcriptome-wide analysis of changes in the fetal placenta associated with prenatal arsenic exposure in the New Hampshire Birth Cohort Study



Increasing evidence suggests that prenatal exposure to arsenic, even at common environmental levels, adversely affects child health. These adverse effects include impaired fetal growth, which can carry serious health implications lifelong. However, the mechanisms by which arsenic affects fetal health and development remain unclear.


We addressed this question using a group of 46 pregnant women selected from the New Hampshire Birth Cohort Study (NHBCS), a US cohort exposed to low-to-moderate arsenic levels in drinking water through the use of unregulated private wells. Prenatal arsenic exposure was assessed using maternal urine samples taken at mid-gestation. Samples of the fetal portion of the placenta were taken from the base of the umbilical cord insertion at the time of delivery, stored in RNAlater and frozen. We used RNA sequencing to analyze changes in global gene expression in the fetal placenta associated with in utero arsenic exposure, adjusting for maternal age. Gene set enrichment analysis and enrichment mapping were then used to identify biological processes represented by the differentially expressed genes. Since our previous analyses have identified considerable sex differences in placental gene expression associated with arsenic exposure, we analyzed male and female samples separately.


At FDR < 0.05, no genes were differentially expressed in female placenta, while 606 genes were differentially expressed in males. Genes showing the most significant associations with arsenic exposure in females were LEMD1 and UPK3B (fold changes 2.51 and 2.48), and in males, FIBIN and RANBP3L (fold changes 0.14 and 0.15). In gene set enrichment analyses, at FDR < 0.05, a total of 211 gene sets were enriched with differentially expressed genes in female placenta, and 154 in male placenta. In female but not male placenta, 103 of these gene sets were also associated with reduced birth weight.


Our results reveal multiple biological functions in the fetal placenta that are potentially affected by increased arsenic exposure, a subset of which is sex-dependent. Further, our data suggest that in female infants, the mechanisms underlying the arsenic-induced reduction of birth weight may involve activation of stress response pathways.

Peer Review reports


An abundance of epidemiological studies have linked prenatal arsenic exposure to a range of adverse infant health outcomes, including spontaneous abortion and neonatal mortality; reduced birth weight; and increased risk of infections in infanthood (reviewed in [1]). Although the effects of high-level exposure are more severe, reports suggest that even levels close to the World Health Organization’s recommended drinking water limit of 10 μg/L, which are estimated to affect millions of people worldwide [2], can be harmful. The New Hampshire Birth Cohort Study (NHBCS) is an ongoing cohort study that was initiated in 2009 to elucidate the effects of such common levels of prenatal exposure to arsenic on maternal and infant health [3]. Participants in the study are pregnant women who use an unregulated private water supply in a US region with low to moderate groundwater arsenic levels. Analyses of the NHBCS have already revealed associations of arsenic exposure with fetal growth measures including birth weight [4]; childhood infections [5]; and various physiological and molecular changes in both the cord blood and placenta [6,7,8,9,10].

Birth weight, a measure of fetal growth, can be an important indicator of risk of both childhood and adult disease conditions including neurocognitive disorders, diabetes, hypertension, and renal disease [11]. Several epidemiological studies have found decreased birth weight associated with greater in utero exposure to arsenic [12, 13], and this has been observed even at relatively low levels of exposure [4, 14]. However, thus far, the mechanism by which arsenic impairs fetal growth is unclear. We have addressed this question by analyzing how arsenic exposure is associated with gene expression in the fetal placenta. The fetal placenta plays a central role in the regulation of fetal growth, controlling the flow of nutrients and oxygen, producing essential hormones, and acting as a protective barrier. Moreover, it accumulates arsenic to up to three times the level in maternal blood [15]. Thus, it is likely that the effects of arsenic on the fetal placenta play a major role in the etiology of low birth weight caused by prenatal exposure. We used RNA sequencing (RNA-seq) as an unbiased, transcriptome-wide approach to identify genes whose expression in the fetal placenta is related to arsenic exposure. Further, using bioinformatic analyses, we identified biological processes related to arsenic exposure, and to birth weight, and used these data to identify potential mechanisms through which arsenic impacts fetal placental function to reduce fetal growth and birth weight. In this initial study, we focused on a group of 46 infants with the lowest and highest levels of prenatal arsenic exposure among a subcohort of the NHBCS. Based on the sexual dimorphism observed in previous analyses of the NHBCS [16, 17], we analyzed the placentas of male and female infants separately.

Materials and methods

The New Hampshire birth cohort study (NHBCS)

The study cohort was selected from a subcohort of 312 mother-child pairs who were enrolled in the New Hampshire Birth Cohort Study (NHBCS) [3] between February 2012 and September 2013. Participants were English-speaking, mentally competent women between 18 and 45 years of age, pregnant with a singleton infant, whose home water supply was from a private, unregulated well at their home. Demographic data, pregnancy history and outcome, and lifestyle factor information were collected using prenatal and delivery records and questionnaires. All subjects provided written informed consent in accordance with the requirements of the Institutional Review Board of Dartmouth College.

Study cohort

For this study, we initially aimed to select 12 infants of each sex with the highest prenatal arsenic exposure within our NHBCS subcohort, and the 12 infants of each sex with the lowest exposure, based on maternal urinary arsenic levels excluding arsenobetaine (U-As, as detailed under “Urine sample collection and arsenic measurement” below). Power analysis, performed using the RNASeqPower R package [18], indicated adequate statistical power (80%) to detect a 2.1-fold change in gene expression using this sample size. One male sample was later found to not have an available urine sample, and one female sample (with low arsenic exposure) was found to be an outlier by principal component analysis of gene expression, as described below. After this exclusion, median average maternal gestational U-As levels were as follows: male high exposure group: 21.04 μg/L (IQR 25.55 μg/L); male low exposure group: 0.67 μg/L (IQR 0.31 μg/L); female high exposure group: 13.81 μg/L (IQR 8.87 μg/L); female low exposure group: 0.74 μg/L (IQR 0.21 μg/L). Other demographic details are provided in Table 1.

Table 1 Demographic details of the study cohort

Arsenic exposure assessment

For assessment of prenatal arsenic exposure in this study, we used total maternal urinary arsenic, excluding arsenobetaine (U-As), measured from spot urine samples at mid-gestation. An advantage of using this measure of exposure is that, due to the relatively high levels of arsenic in urine, we were able to measure species levels individually and remove arsenobetaine, which is considered non-toxic [19]. To minimize variation in arsenic intake from household water during gestation, participants were selected for the NHBCS who had not changed residence since their last menstrual period, and were not planning to move residence before delivery. Household water arsenic was previously shown to be strongly associated with maternal urinary arsenic in the NHBCS [3], consistent with it being a major source of exposure.

Urine sample collection and arsenic measurement

Details of sample collection and arsenic measurement have been described previously [6]. Briefly, maternal spot urine samples were collected at approximately 24–28 weeks of gestation, into acid-washed containers containing diammonium diethyldithiocarbamate to stabilize trivalent methyl arsenic species [20], and frozen at − 80 °C until analysis (within 24 h of collection). Samples were analyzed for levels of five individual arsenic species: arsenite (AsIII), arsenate (AsV), dimethylarsinic acid (DMAV), monomethylarsonic acid (MMAV), and arsenobetaine, using high-performance liquid chromatography inductively coupled plasma mass spectrometry (ICP-MS). The detection limits were: 0.15 μg/L for AsIII, 0.10 μg/L for AsV, 0.14 μg/L for MMA, and 0.11 μg/L for DMA. Numbers of samples below the detection limits were as follows: AsIII: 29, AsV: 0, DMAV: 30, and MMAV: 22. We calculated the sum of AsIII, AsV, DMAV and MMAV, with arsenobetaine excluded, as it is considered non-toxic and is not metabolized, and used the resulting value, denoted U-As, as a measure of arsenic exposure. Proportions of different arsenic species did not vary significantly between high and low groups, except in the case of AsV, which was slightly higher among low exposure than high exposure males (P = 0.032, see Additional file 1).

Placenta biopsy and gene profiling

Placental biopsies were taken by medical staff in the delivery room. Biopsies measured approximately 1 cm deep and 1–2 cm across, and were taken at the base of the umbilical cord insertion to minimize heterogeneity. Any maternal decidua was removed. Samples were immediately placed in tubes containing RNAlater (Life Technologies), and then frozen at − 80 °C within 24 h. Total RNA was extracted using the RNA/DNA extraction kit (Norgen Biotek, Thorold, ON), quantified using a NanoDrop spectrophotometer, and stored at − 80 °C. RNA quality was determined using an Agilent 2100 Bioanalyzer, and global gene expression analysis was performed by RNA sequencing (RNA-seq) at the Oncogenomics Core Facility at the University of Miami. Specifically, TruSeq Stranded Total RNA-seq Library Prep kits (Illumina) were used to convert total RNA to cDNA libraries, which were then sequenced using the Genome Analyzer IIx system (Illumina). The male and female RNA-seq assays were performed separately; however, identical procedures were used for both assays.

Bioinformatic analysis

A summary of the analysis pathway is provided in Fig. 1. Specifically, sequencing fragments were aligned to the human reference genome (GRCh38) using the STAR algorithm [21], and reads mapping to genomic features (transcribed RNA) were counted using featureCounts [22]. Differential expression was assessed using DESeq2 [23]. Specifically, samples were dichotomized into high and low arsenic exposure groups. A linear model: raw read counts ~ groups + maternal age at enrollment + batch was then fitted to each gene, and Wald significance tests and their corresponding p-values were used to assess differences in expression related to arsenic exposure, adjusted for maternal age at enrollment and assay batch [24]. For analyses of differences in expression related to birth weight, birth weight (in kilograms) was considered as a continuous variable. Gestational age, enrollment age and assay batch were also included as covariates in the linear models for association tests with birth weight. The adjustments (excluding assay batch) were based on previous analyses of the NHBCS, in which a range of variables (maternal age, maternal smoking status (never, former, current), maternal education level, infant birth weight, infant sex and gestational age) were assessed as potential confounders using a series of linear regression models [25]. In our data, we found the above variables to be associated with U-As/birth weight and gene expression. To account for multiple comparisons, we computed the False Discovery Rate (FDR) [26].

Fig. 1
figure 1

RNA-seq data analysis. Schematic diagram summarizing the analyses described in this report. DE; differentially expressed, GSEA; gene set enrichment analysis, EMap; enrichment map

To further evaluate the potential molecular and functional effects of arsenic exposure, we performed gene set enrichment analysis (GSEA) [27]. This involved identifying gene sets (representing biological processes) enriched with differentially expressed genes. The “canonical pathways” (CP) gene set collection from the Molecular Signatures Database (MSigDB), comprising 1330 gene sets, was used. Gene sets were ranked by normalized enrichment score (NES), a linear measurement of the degree to which a gene set is overrepresented at the top or bottom of a list of genes ranked according to their differential expression [27]. Enrichment Maps [28] were then generated to visualize the major biological themes of the GSEA results, by grouping the identified gene sets into clusters based on common genes.

RNA-seq quality control data are provided in Additional file 2. An average of 58.2 million and 68.0 million reads were sequenced from female and male placenta samples respectively. Of these, 96.55 and 98.06% respectively were aligned to the reference genome, and 49.21 and 46.09% respectively were assigned to genes (mapped to transcribed RNA). Principal component analysis (PCA) of the gene expression data revealed one female sample (with low arsenic exposure) to be an outlier (Additional file 3A); this sample was therefore excluded from subsequent analyses. PCA did not reveal any outliers among the male samples (Additional file 3B), but showed that the male and female samples clustered separately (Additional file 3C).


The RNA-seq data were analyzed as summarized in Fig. 1. First, we identified genes that were differentially expressed (DE) between high- and low-U-As groups (Additional file 4). At P < 0.05, 1748 genes were DE in female placenta, and 2438 genes in male placenta. The top DE genes were LEMDI and UPK3B in females (fold changes 2.51 and 2.48), and FIBIN and RANBP3L in males (fold changes 0.14 and 0.15). A total of 458 DE genes at P < 0.05 were common to males and females (data not shown). After FDR adjustment for multiple comparisons, no genes were DE in females at FDR < 0.05, while 606 genes were differentially expressed in males (Additional file 4). To validate these results, we performed qPCR for four of the top DE genes in females using new tissue samples from the same placentas, and found that the results were largely in keeping with our RNA-seq data (Additional file 5).

Next, we performed gene set enrichment analysis (GSEA) to identify gene sets enriched with DE genes in our high versus low U-As groups. Such gene sets represent biological processes associated with arsenic exposure. At a significance level of FDR < 0.05, in females, 87 gene sets were upregulated, and 124 gene sets were downregulated, while in males, 7 gene sets were upregulated, and 146 gene sets were downregulated (Additional file 6A-D). Table 2 shows the top five gene sets in each group, based on normalized enrichment scores (NES). In females (Table 2A), the top upregulated gene sets were related to cellular respiration and the ubiquitin proteasome system, while downregulated gene sets involved protein synthesis and the extracellular matrix. In males (Table 2B), the top upregulated gene sets were related to the unfolded protein response and endocytosis/recycling of cell surface proteins, while downregulated gene sets involved the extracellular matrix, as in females, and smooth muscle contraction.

Table 2 Top gene sets associated with arsenic exposure

We then compared the male and female GSEA results (Additional file 6E). Six gene sets were upregulated and 61 gene sets were downregulated at high U-As in both male and female placentas.

To better characterize these results, we generated enrichment maps (EMaps) [28]. These maps group the gene sets identified by GSEA into clusters, to identify the major biological themes associated with arsenic exposure in female and/or male fetal placentas (Additional files 7 and 8). Comparison of the clusters in each EMap revealed both similarities and differences between the sexes. For example, in both male and female placentas, gene set clusters related to the extracellular matrix (ECM) and muscle contraction were downregulated, and a cluster related to the unfolded protein response was upregulated. In contrast, several gene set clusters were unique to females, in particular a large upregulated cluster related to proteasomal degradation, while male placentas showed downregulation of a large gene set cluster related to a range of cell processes/components including G protein coupled receptor (GPCR) signaling, and T cell receptor activation, only a subset of which was also downregulated in females. In addition, clusters related to autoimmune responses/inflammation, and interferon gamma signaling, were largely unique to male placentas. These data suggest that multiple cellular functions in the fetal placenta may be altered by increased arsenic exposure at the common environmental levels encountered in our study, and that a subset of these potential effects are fetal sex-dependent.

Previous studies have indicated that higher prenatal arsenic exposure may decrease birth weight [4, 12,13,14]; however, the underlying mechanisms are unknown. Thus, we next used our RNA-seq data to identify biological processes that may mediate adverse effects of arsenic on fetal growth, by comparing the gene sets identified by GSEA at high versus low U-As with those that were associated with birth weight (as a continuous variable). In female placenta, at FDR < 0.05, 183 gene sets were associated with birth weight (Additional file 9A), and 103 of these were also associated with U-As (Additional file 10A). Of these, 51 gene sets were upregulated and 16 downregulated at high U-As and low birth weight, representing a total of 66 biological processes potentially mediating female birth weight - arsenic exposure associations. The top 10 of these included downregulation of the core matrisome (ECM components), and upregulation of electron transport chain and proteasome components (Table 3). An EMap was generated of the 66 gene sets, revealing 7 gene-set clusters (Fig. 2). The upregulated clusters were related to the ubiquitin proteasome pathway, cellular respiration, mRNA synthesis/processing, and protein glycosylation; while the downregulated clusters involved ECM components and interactors, endothelial G-protein-coupled receptors, and growth factors involved in wound healing and angiogenesis.

Table 3 Top 10 gene sets enriched at high arsenic exposure and reduced birth weight in female fetal placenta
Fig. 2
figure 2

Gene sets differentially related to arsenic exposure and birth weight in female fetal placenta. Enrichment map of gene sets that are associated with increased arsenic exposure and reduced birth weight in female placenta (with a significance level of FDR < 0.05). Clusters thus represent biological themes that may underlie female fetal growth inhibition by arsenic. Red circles; upregulated gene sets, blue circles; downregulated gene sets, large yellow circles; gene set clusters, green lines indicate overlapping genes between gene sets. Clusters were labeled based on the most common themes of the gene sets they comprise. Singleton gene sets are omitted (see Additional file 10A for a complete list of gene sets)

In male fetal placenta, 147 gene sets were associated with birth weight at FDR < 0.05 (Additional file 9B), and 100 of these were also negatively associated with U-As (Additional file 10B). However, no gene sets were positively associated with U-As, and none showed differential associations with U-As and birth weight. Therefore, in contrast to females, this analysis did not reveal any biological processes that might mediate decreases in birth weight related to arsenic exposure in male fetal placenta.

To explore the sexual dimorphism we observed in the biological processes linking arsenic exposure and birth weight, we compared the gene sets associated with birth weight in male and female placentas (Additional file 9C). Interestingly, at FDR < 0.05, 22 gene sets were associated with birth weight in both sexes, but 16 of these showed opposite directions of association, namely negative in males but positive in females. Among these 16 gene sets, 14 were also negatively associated with U-As in both sexes, and comprised all of the downregulated gene-set clusters potentially linking increased arsenic exposure with lower female infant birth weight shown in the EMap in Fig. 2. These clusters were related to ECM components and interactors, endothelial G-protein-coupled receptors, and wound healing and angiogenesis.


Our RNA-seq analysis identified numerous genes whose expression in fetal placenta appeared to associate with arsenic exposure in a sex-dependent manner, although these associations did not withstand adjustment for multiple comparisons in females. Further, GSEA and enrichment mapping revealed both common and sex-specific biological processes, represented by gene sets and gene-set clusters, that were associated with arsenic exposure. The most significant processes common to both sexes included upregulation of the unfolded protein response (UPR). The UPR is activated by endoplasmic reticulum (ER) stress, an accumulation of unfolded proteins in the ER. Arsenic has been shown in human and mouse cell culture studies to activate the UPR [29, 30], likely by interfering with oxidative protein folding [31]. Also common to both sexes was downregulation of ECM components, endothelial GPCRs, and wound healing and angiogenesis. These results are somewhat reminiscent of a study in mice that found reduced expression of ECM genes, and disruption of arterial ECM, in the heart and lungs following chronic arsenic exposure [32]. These gene sets were also inversely associated with birth weight in female placentas, suggesting that reduced expression of ECM components, endothelial GPCRs, and factors involved in wound healing and angiogenesis may contribute to lower female infant birth weight caused by arsenic exposure. However, interestingly, we found that, in male placentas, these gene sets were positively associated with birth weight. One explanation for this is that these responses are not alone sufficient to affect fetal growth, but that, in the females in our analysis, they have occurred in parallel with other functional changes that directly reduce fetal growth.

Such functional changes may be represented by four gene set clusters that were upregulated with higher arsenic exposure and lower birth weight in female, but not male, placentas. These clusters represent the proteasome, cellular respiration, mRNA synthesis/processing, and protein glycosylation. Some of these pathways may be activated in response to arsenic-induced oxidative stress. Numerous studies have found arsenic to induce oxidative stress and reactive oxygen species production (reviewed in [33, 34]). To counteract the effects of oxidative stress, cells activate the antioxidant response pathway. This involves activation of the NFE2L2 (Nrf2) transcription factor, leading to transcription of a suite of genes that have antioxidant response elements (AREs) in their promoters. Notably, these include many subunits of the 26S proteasome, which is upregulated to remove oxidatively damaged proteins [35]. NFE2L2 has also been shown to stimulate mitochondrial biogenesis [36], which may explain the increased expression of electron transport chain components. Increased transcription and synthesis of NFE2L2 targets may also explain the observed upregulation of mRNA synthesis and processing, and protein glycosylation pathways. Increased protein glycosylation may also be a feedback response to impaired protein folding in the ER in response to arsenic, as mentioned above.

Our analyses also identified some biological processes that were associated with arsenic exposure exclusively in male placentas. Among the most significant of these were upregulation of genes related to transferrin endocytosis, and downregulation of those related to autoimmune responses/inflammation, and interferon gamma signaling. Transferrin is an iron transport protein that is expressed in both the cytotrophoblast and syncytiotrophoblast of the fetal placenta [37,38,39]. In support of our findings, a previous study reported increased transferrin expression in the syncytiotrophoblast of pregnancies complicated by maternal drug abuse, gestational diabetes or pregnancy-induced hypertension, suggesting that this may represent a response to intrauterine stress [37]. Placental expression of the transferrin receptor was also increased in conditions of iron deficiency [40]. Interferon gamma (IFNG) is a proinflammatory cytokine that is produced by various immune cells, including natural killer and CD4+ T helper 1 (Th1) cells [41], and plays an central role in the development of autoimmunity [42]. IFNG and its receptors are also expressed by the trophoblast of the fetal placenta; however this tissue shows a dampened response to IFNG [41], and the fetal/neonatal immune system tends to be tolerogenic and Th2-biased, with minimal IFNG production, to avoid responses to maternal alloantigens [43, 44]. Our results appear to suggest that, in male fetal placenta, arsenic exposure increases this anti-inflammatory bias. This finding aligns with a previous study of the NHBCS, which found increased numbers of Th2-type cells in cord blood at high arsenic exposure, although sex differences were not observed [10]. Studies of adults chronically exposed to arsenic have also shown immune effects, including reduced expression of IFNG [45]. Arsenic exposure, including prenatal exposure [5, 46], has been associated with increased susceptibility to various infections, and some studies have recorded sex differences in such associations, with males tending to be more susceptible (reviewed in [47]).

Limitations of this study include the small sample size, the large range of U-As values within the “high” and “low”-arsenic groups, and the variation in these ranges between male and female placenta samples. In addition, RNA sequencing of the male and female samples were performed separately. Therefore, although the samples were treated and analyzed in an identical manner, it may be that a subset of the differences observed is due to experimental variation, rather than sex. A further limitation is that, in this study, we did not examine levels of other metal toxicants in maternal urine, so cannot rule out the possibility of a confounding effect by other elements. Another potential confounding factor is that the participants are exposed to arsenic from different sources, e.g. rice, seafood, and drinking water, which contains different proportions of arsenic species and therefore may have different effects. The relative proportion of AsV in males was slightly higher in the low arsenic group (Additional file 1), which may have influenced some of the noted associations. Models were adjusted for RNA-seq batch, maternal age, and in the birth weight analysis, gestational age. Additional potential confounders, i.e., maternal smoking status, and maternal education level were not found to significantly associate with U-As in our previous analyses of the NHBCS [25]. Another important point to acknowledge is that we are using gene expression after delivery to draw conclusions about prenatal gene expression. There are likely to be numerous changes occurring in the fetal placenta within this period, and therefore caution must be taken in the interpretation of our findings.

In this study, we chose to use maternal urinary arsenic (excluding arsenobetaine), measured at mid-gestation, as our measure of prenatal arsenic exposure. To minimize variation in arsenic exposure between mid-gestation and delivery, participants were selected who did not plan to move residence between mid-gestation and delivery. Thus, the household water supply, a major component of overall arsenic exposure [3], remained constant. An advantage of using urinary arsenic as the exposure measure was the ability to measure individual species, and thus remove arsenobetaine from our calculation, which is considered non-toxic [19]. In contrast, placental arsenic is present in much lower levels, and therefore, when it is used for exposure assessment, arsenobetaine cannot be accounted for. However, clearly, mid-gestational U-As has a number of disadvantages: primarily the time difference between urine sampling and placental expression analysis, and the fact that only a single sample was taken, as well as the increased metabolism of arsenic during pregnancy [48]. These factors must be acknowledged as important limitations of the current study. Future studies should include similar studies using placental arsenic for assessment of prenatal exposure.


Our results suggest that common levels of arsenic exposure are associated with multiple changes in the human placental transcriptome, a subset of which was sex-specific. Further, we identified some potential sex-dependent mechanisms for the known adverse effects of arsenic on birth weight. Overall, our findings offer insights into potential mechanisms through which prenatal arsenic exposure may impact the fetal placenta in a sex-dependent manner to affect fetal health and development, which may be further explored in future studies.

Availability of data and materials

The dataset supporting the conclusions of this article is included within the article and its additional files.



Enrichment map


False discovery rate


Gene set enrichment analysis


New Hampshire Birth Cohort Study




Total maternal urinary arsenic concentration (excluding arsenobetaine)


  1. Farzan SF, Karagas MR, Chen Y. In utero and early life arsenic exposure in relation to long-term health and disease. Toxicol Appl Pharmacol. 2013;272:384–90.

    Article  CAS  Google Scholar 

  2. Nordstrom DK. Public health. Worldwide occurrences of arsenic in ground water. Science. 2002;296:2143–5.

    Article  CAS  Google Scholar 

  3. Gilbert-Diamond D, Cottingham KL, Gruber JF, Punshon T, Sayarath V, Gandolfi AJ, Baker ER, Jackson BP, Folt CL, Karagas MR. Rice consumption contributes to arsenic exposure in US women. Proc Natl Acad Sci U S A. 2011;108:20656–60.

    Article  CAS  Google Scholar 

  4. Gilbert-Diamond D, Emond JA, Baker ER, Korrick SA, Karagas MR. Relation between in utero arsenic exposure and birth outcomes in a cohort of mothers and their newborns from New Hampshire. Environ Health Perspect. 2016;124:1299–307.

    Article  CAS  Google Scholar 

  5. Farzan SF, Li Z, Korrick SA, Spiegelman D, Enelow R, Nadeau K, Baker E, Karagas MR. Infant infections and Respiratory symptoms in relation to in utero arsenic exposure in a U.S. cohort. Environ Health Perspect. 2016;124:840–7.

    Article  CAS  Google Scholar 

  6. Koestler DC, Avissar-Whiting M, Houseman EA, Karagas MR, Marsit CJ. Differential DNA methylation in umbilical cord blood of infants exposed to low levels of arsenic in utero. Environ Health Perspect. 2013;121:971–7.

    Article  Google Scholar 

  7. Farzan SF, Brickley EB, Li Z, Gilbert-Diamond D, Gossai A, Chen Y, Howe CG, Palys T, Karagas MR. Maternal and infant inflammatory markers in relation to prenatal arsenic exposure in a U.S. pregnancy cohort. Environ Res. 2017;156:426–33.

    Article  CAS  Google Scholar 

  8. Gossai A, Lesseur C, Farzan S, Marsit C, Karagas MR, Gilbert-Diamond D. Association between maternal urinary arsenic species and infant cord blood leptin levels in a New Hampshire pregnancy cohort. Environ Res. 2015;136:180–6.

    Article  CAS  Google Scholar 

  9. Green BB, Karagas MR, Punshon T, Jackson BP, Robbins DJ, Houseman EA, Marsit CJ. Epigenome-wide assessment of DNA methylation in the placenta and arsenic exposure in the New Hampshire birth cohort study (USA). Environ Health Perspect. 2016;124:1253–60.

    Article  CAS  Google Scholar 

  10. Nadeau KC, Li Z, Farzan S, Koestler D, Robbins D, Fei DL, Malipatlolla M, Maecker H, Enelow R, Korrick S, Karagas MR. In utero arsenic exposure and fetal immune repertoire in a US pregnancy cohort. Clin Immunol. 2014;155:188–97.

    Article  CAS  Google Scholar 

  11. Malin GL, Morris RK, Riley RD, Teune MJ, Khan KS. When is birthweight at term (>/=37 weeks’ gestation) abnormally low? A systematic review and meta-analysis of the prognostic and predictive ability of current birthweight standards for childhood and adult outcomes. BJOG. 2015;122:634–42.

    Article  CAS  Google Scholar 

  12. Huyck KL, Kile ML, Mahiuddin G, Quamruzzaman Q, Rahman M, Breton CV, Dobson CB, Frelich J, Hoffman E, Yousuf J, et al. Maternal arsenic exposure associated with low birth weight in Bangladesh. J Occup Environ Med. 2007;49:1097–104.

    Article  CAS  Google Scholar 

  13. Llanos MN, Ronco AM. Fetal growth restriction is related to placental levels of cadmium, lead and arsenic but not with antioxidant activities. Reprod Toxicol. 2009;27:88–92.

    Article  CAS  Google Scholar 

  14. Rahman A, Vahter M, Smith AH, Nermell B, Yunus M, El Arifeen S, Persson LA, Ekstrom EC. Arsenic exposure during pregnancy and size at birth: a prospective cohort study in Bangladesh. Am J Epidemiol. 2009;169:304–12.

    Article  Google Scholar 

  15. Concha G, Vogler G, Lezcano D, Nermell B, Vahter M. Exposure to inorganic arsenic metabolites during early human development. Toxicol Sci. 1998;44:185–90.

    Article  CAS  Google Scholar 

  16. Winterbottom EF, Fei DL, Koestler DC, Giambelli C, Wika E, Capobianco AJ, Lee E, Marsit CJ, Karagas MR, Robbins DJ. GLI3 links environmental arsenic exposure and human fetal growth. EBioMedicine. 2015;2:536–43.

    Article  Google Scholar 

  17. Winterbottom EF, Koestler DC, Fei DL, Wika E, Capobianco AJ, Marsit CJ, Karagas MR, Robbins DJ. The aquaglyceroporin AQP9 contributes to the sex-specific effects of in utero arsenic exposure on placental gene expression. Environ Health. 2017;16:59.

    Article  CAS  Google Scholar 

  18. Hart SN, Therneau TM, Zhang Y, Poland GA, Kocher JP. Calculating sample size estimates for RNA sequencing data. J Comput Biol. 2013;20:970–8.

    Article  CAS  Google Scholar 

  19. Navas-Acien A, Francesconi KA, Silbergeld EK, Guallar E. Seafood intake and urine concentrations of total arsenic, dimethylarsinate and arsenobetaine in the US population. Environ Res. 2011;111:110–8.

    Article  CAS  Google Scholar 

  20. Wang Z, Zhou J, Lu X, Gong Z, Le XC. Arsenic speciation in urine from acute promyelocytic leukemia patients undergoing arsenic trioxide treatment. Chem Res Toxicol. 2004;17:95–103.

    Article  CAS  Google Scholar 

  21. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, Batut P, Chaisson M, Gingeras TR. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21.

    Article  CAS  Google Scholar 

  22. Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30:923–30.

    Article  CAS  Google Scholar 

  23. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.

    Article  CAS  Google Scholar 

  24. Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:3.

    Article  Google Scholar 

  25. Fei DL, Koestler DC, Li Z, Giambelli C, Sanchez-Mejias A, Gosse JA, Marsit CJ, Karagas MR, Robbins DJ. Association between in utero arsenic exposure, placental gene expression, and infant birth weight: a US birth cohort study. Environ Health. 2013;12:58.

    Article  CAS  Google Scholar 

  26. Benjamini Y, Hochberg Y. Controlling the false discovery rate - a practical and powerful approach to multiple testing. J Royal Stat Soc Series B Methodol. 1995;57:289–300.

    Google Scholar 

  27. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545–50.

    Article  CAS  Google Scholar 

  28. Merico D, Isserlin R, Stueker O, Emili A, Bader GD. Enrichment map: a network-based method for gene-set enrichment visualization and interpretation. PLoS One. 2010;5:e13984.

    Article  CAS  Google Scholar 

  29. Weng CY, Chiou SY, Wang L, Kou MC, Wang YJ, Wu MJ. Arsenic trioxide induces unfolded protein response in vascular endothelial cells. Arch Toxicol. 2014;88:213–26.

    Article  CAS  Google Scholar 

  30. Srivastava RK, Li C, Chaudhary SC, Ballestas ME, Elmets CA, Robbins DJ, Matalon S, Deshane JS, Afaq F, Bickers DR, Athar M. Unfolded protein response (UPR) signaling regulates arsenic trioxide-mediated macrophage innate immune function disruption. Toxicol Appl Pharmacol. 2013;272:879–87.

    Article  CAS  Google Scholar 

  31. Ramadan D, Rancy PC, Nagarkar RP, Schneider JP, Thorpe C. Arsenic(III) species inhibit oxidative protein folding in vitro. Biochemistry. 2009;48:424–32.

    Article  CAS  Google Scholar 

  32. Hays AM, Lantz RC, Rodgers LS, Sollome JJ, Vaillancourt RR, Andrew AS, Hamilton JW, Camenisch TD. Arsenic-induced decreases in the vascular matrix. Toxicol Pathol. 2008;36:805–17.

    Article  CAS  Google Scholar 

  33. Jomova K, Valko M. Advances in metal-induced oxidative stress and human disease. Toxicology. 2011;283:65–87.

    Article  CAS  Google Scholar 

  34. Sinha D, Biswas J, Bishayee A. Nrf2-mediated redox signaling in arsenic carcinogenesis: a review. Arch Toxicol. 2013;87:383–96.

    Article  CAS  Google Scholar 

  35. Livneh I, Cohen-Kaplan V, Cohen-Rosenzweig C, Avni N, Ciechanover A. The life cycle of the 26S proteasome: from birth, through regulation and function, and onto its death. Cell Res. 2016;26:869–85.

    Article  CAS  Google Scholar 

  36. Dinkova-Kostova AT, Abramov AY. The emerging role of Nrf2 in mitochondrial function. Free Radic Biol Med. 2015;88:179–88.

    Article  CAS  Google Scholar 

  37. Kralova A, Svetlikova M, Madar J, Ulcova-Gallova Z, Bukovsky A, Peknicova J. Differential transferrin expression in placentae from normal and abnormal pregnancies: a pilot study. Reprod Biol Endocrinol. 2008;6:27.

    Article  CAS  Google Scholar 

  38. Morris Buus R, Boockfor FR. Transferrin expression by placental trophoblastic cells. Placenta. 2004;25:45–52.

    Article  CAS  Google Scholar 

  39. Verrijt CE, Kroos MJ, Verhoeven AJ, van Eijk HG, van Dijk JP. Transferrin in cultured human term cytotrophoblast cells: synthesis and heterogeneity. Mol Cell Biochem. 1997;173:177–81.

    Article  CAS  Google Scholar 

  40. Gambling L, Danzeisen R, Fosset C, Andersen HS, Dunford S, Srai SK, HJ MC. Iron and copper interactions in development and the effect on pregnancy outcome. J Nutr. 2003;133:1554S–6S.

    Article  CAS  Google Scholar 

  41. Murphy SP, Tayade C, Ashkar AA, Hatta K, Zhang J, Croy BA. Interferon gamma in successful pregnancies. Biol Reprod. 2009;80:848–59.

    Article  CAS  Google Scholar 

  42. Pollard KM, Cauvi DM, Toomey CB, Morris KV, Kono DH. Interferon-gamma and systemic autoimmunity. Discov Med. 2013;16:123–31.

    Google Scholar 

  43. Adkins B. Development of neonatal Th1/Th2 function. Int Rev Immunol. 2000;19:157–71.

    Article  CAS  Google Scholar 

  44. Garcia AM, Fadel SA, Cao S, Sarzotti M. T cell immunity in neonates. Immunol Res. 2000;22:177–90.

    Article  CAS  Google Scholar 

  45. Biswas R, Ghosh P, Banerjee N, Das JK, Sau T, Banerjee A, Roy S, Ganguly S, Chatterjee M, Mukherjee A, Giri AK. Analysis of T-cell proliferation and cytokine secretion in the individuals exposed to arsenic. Hum Exp Toxicol. 2008;27:381–6.

    Article  CAS  Google Scholar 

  46. Farzan SF, Korrick S, Li Z, Enelow R, Gandolfi AJ, Madan J, Nadeau K, Karagas MR. In utero arsenic exposure and infant infection in a United States cohort: a prospective study. Environ Res. 2013;126:24–30.

    Article  CAS  Google Scholar 

  47. Ferrario D, Gribaldo L, Hartung T. Arsenic exposure and Immunotoxicity: a review including the possible influence of age and sex. Curr Environ Health Rep. 2016;3:1–12.

    Article  CAS  Google Scholar 

  48. Gardner RM, Nermell B, Kippler M, Grander M, Li L, Ekstrom EC, Rahman A, Lonnerdal B, Hoque AM, Vahter M. Arsenic methylation efficiency increases during the first trimester of pregnancy independent of folate status. Reprod Toxicol. 2011;31:210–8.

    Article  CAS  Google Scholar 

Download references


We thank Dr. Sion Williams, Yoslayma Cardentey, and other staff at the Oncogenomics Core Facility at the University of Miami for their technical expertise. Further, we are grateful to the members of the Robbins, Capobianco, Marsit, and Karagas laboratories for helpful discussions and technical assistance during the course of this work. We thank all participants of the New Hampshire Birth Cohort Study.


This publication was made possible by U.S. Environmental Protection Agency (US EPA) grant RD-83544201, and National Institute of Environmental Health Sciences (NIEHS) grant P01 ES022832. Its contents are solely the responsibility of the grantee and do not necessarily represent the official views of the US EPA. Further, the US EPA does not endorse the purchase of any commercial products or services mentioned in the publication.

Author information

Authors and Affiliations



EFW and DJR designed research; EFW performed qPCR experiments; YB, XS, LW, XC, and EFW analyzed data; EFW, AJC, CJM, MRK, and DJR wrote the paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Emily F. Winterbottom.

Ethics declarations

Ethics approval and consent to participate

All research involving human participants was approved by The Committee for the Protection of Human Subjects at Dartmouth College. Written informed consent was required for enrollment.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1.

Relative proportions of arsenic species in high and low arsenic groups. Upper and lower ends of boxes indicate the 25th and 75th percentiles, respectively, and black band represents the median. Error bars represent minimum and maximum values, excluding outliers, which are depicted as dots. P values are based on a Wilcoxon signed rank test.

Additional file 2.

RNA-seq quality control data for A) female, and B) male placenta samples.

Additional file 3.

Principal component analyses. A) female placenta samples; B) male placenta samples; C) male and female samples.

Additional file 4.

Differentially expressed genes at high versus low U-As adjusted for maternal age at enrollment and batch. A) female placenta, B) male placenta.

Additional file 5.

qPCR validation of selected differentially expressed genes. RNA was extracted from repeat samples of the same placentas used for the RNA-seq analysis, and qPCR was performed using TaqMan probes designed against a subset of the top 10 differentially expressed genes in female placentas. Error bars show SEM. * P < 0.05.

Additional file 6.

Gene sets enriched with differentially expressed genes at high versus low U-As at FDR < 0.05, in A-B) female and C-D) male fetal placenta, and E) comparison.

Additional file 7.

Gene set enrichment analysis of arsenic-exposed female fetal placenta. Enrichment map showing gene sets enriched with differentially expressed genes at high versus low arsenic exposure (U-As levels) in female fetal placenta. Clustered gene sets with a significance level of FDR < 0.05 are shown. The “canonical pathways” gene set collection from MSigDB was used. Red circles; upregulated gene sets, blue circles; downregulated gene sets, large yellow circles; gene set clusters, green lines indicate overlapping genes between gene sets; words that appear most frequently in the gene set titles are shown. Singleton gene sets are omitted (see Additional file 6A-B for a complete list of gene sets).

Additional file 8.

Gene set enrichment analysis of arsenic-exposed male fetal placenta. Enrichment map showing gene sets enriched with differentially expressed genes at high versus low arsenic exposure (U-As levels) in male fetal placenta. Clustered gene sets with a significance level of FDR < 0.05 are shown. The “canonical pathways” gene set collection from MSigDB was used. Red circles; upregulated gene sets, blue circles; downregulated gene sets, large yellow circles; gene set clusters, green lines indicate overlapping genes between gene sets; words that appear most frequently in the gene set titles are shown. Singleton gene sets are omitted (see Additional file 6C-D for a complete list of gene sets).

Additional file 9.

Gene sets enriched with differentially expressed genes at high versus low birth weight at FDR < 0.05, in A) female and B) male fetal placenta, and C) comparison.

Additional file 10.

Gene sets enriched in both birth weight and U-As analyses in A) female and B) male fetal placenta at FDR < 0.05. Gene sets showing opposite directions of association in 9A) are asterisked.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Winterbottom, E.F., Ban, Y., Sun, X. et al. Transcriptome-wide analysis of changes in the fetal placenta associated with prenatal arsenic exposure in the New Hampshire Birth Cohort Study. Environ Health 18, 100 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: