Research | Open | Published:
In utero exposure to maternal smoking is associated with DNA methylation alterations and reduced neuronal content in the developing fetal brain
Epigenetics & Chromatinvolume 10, Article number: 4 (2017)
Intrauterine exposure to maternal smoking is linked to impaired executive function and behavioral problems in the offspring. Maternal smoking is associated with reduced fetal brain growth and smaller volume of cortical gray matter in childhood, indicating that prenatal exposure to tobacco may impact cortical development and manifest as behavioral problems. Cellular development is mediated by changes in epigenetic modifications such as DNA methylation, which can be affected by exposure to tobacco.
In this study, we sought to ascertain how maternal smoking during pregnancy affects global DNA methylation profiles of the developing dorsolateral prefrontal cortex (DLPFC) during the second trimester of gestation. When DLPFC methylation profiles (assayed via Illumina, HM450) of smoking-exposed and unexposed fetuses were compared, no differentially methylated regions (DMRs) passed the false discovery correction (FDR ≤ 0.05). However, the most significant DMRs were hypomethylated CpG Islands within the promoter regions of GNA15 and SDHAP3 of smoking-exposed fetuses. Interestingly, the developmental up-regulation of SDHAP3 mRNA was delayed in smoking-exposed fetuses. Interaction analysis between gestational age and smoking exposure identified significant DMRs annotated to SYCE3, C21orf56/LSS, SPAG1 and RNU12/POLDIP3 that passed FDR. Furthermore, utilizing established methods to estimate cell proportions by DNA methylation, we found that exposed DLPFC samples contained a lower proportion of neurons in samples from fetuses exposed to maternal smoking. We also show through in vitro experiments that nicotine impedes the differentiation of neurons independent of cell death.
We found evidence that intrauterine smoking exposure alters the developmental patterning of DNA methylation and gene expression and is associated with reduced mature neuronal content, effects that are likely driven by nicotine.
Numerous studies have established that maternal smoking during pregnancy is associated with impaired executive function and behavioral problems in the offspring [1–3]. Maternal smoking is associated with altered fetal brain development  and reduced volumes of cortical gray matter in childhood , indicating that exposure to tobacco smoke constituents in utero may impact brain development and subsequently result in neurodevelopmental abnormalities. Offspring exposed to smoking after birth does not exhibit the same adverse trajectories [6, 7], suggesting biologically mediated mechanisms during gestation. Cigarette smoke is a highly complex mixture of more than 5000 chemicals of which approximately 100 are known to be hazardous . Linking specific compound(s) with defined phenotypes has proven difficult. Indirect biological mechanisms caused by cigarette constituents other than nicotine have been proposed, including hypoxia/ischemia and DNA damage . Exposure to nicotine prenatally has a direct impact on brain development. In rodents, prenatal exposure to nicotine is reported to induce abnormal dendritic morphology and reduced synapse density in the cerebral cortex and nucleus accumbens . Additionally, prenatal nicotine exposure during primate brain development up-regulates nicotinic acetylcholine receptors (nAChR), causes cell death, and alters cell size and neurite outgrowth in a regionally dependent manner . Furthermore, nicotine replacement therapy has been suggested to increase the risk for behavioral impairments (for review see [12–17]).
Prenatal exposure to environmental factors such as alcohol  and industrial chemicals (lead, methylmercury, PCBs, reviewed in [19, 20]) often manifests as neurodevelopmental disorders. Epigenetic modifications such as DNA methylation regulate gene activity necessary for cell differentiation . Exposure to tobacco smoke can induce alterations in epigenetic patterning that are associated with a wide spectrum of human diseases including cardiovascular, pulmonary, neurobehavioral disorders and cancer [22–28]. Maternal smoking during pregnancy alters DNA methylation in the blood of newborns  and can cause DNA methylation changes that persist into childhood . In relation to neurological function, differences in DNA methylation have been reported between offspring of smokers and non-smokers in the promoters of catechol-O-methyltransferase (COMT) and monoamine oxidase A (MAOA), genes thought to be involved in nicotine dependence and other neurobehavioral disorders [31, 32]. Further, an increase in the DNA methylation of the brain-derived neurotrophic factor-6 (BDNF-6) promoter/5′UTR has been found in adolescents exposed to maternal smoking during pregnancy . To our knowledge, no studies have directly examined the epigenetic changes of the developing human brain exposed in utero to maternal cigarette smoking. Here we interrogate DNA methylation patterns in the developing cortex of human fetuses exposed to maternal smoking on a genome scale.
Fetuses were from second-trimester elective saline abortions performed for non-medical reasons. Fetal sample groups, exposed (N = 14) and unexposed (N = 10) to maternal smoking, were balanced for sex and gestational age (weeks since the first day of last normal menstrual period) (Table 1). All mothers of the exposed group were active smokers prior to and during pregnancy, whereas no mothers of the unexposed group were active smokers prior to or during pregnancy (Additional file 4: Table S1). Alcohol is a well-described teratogen that affects neuroanatomical development [18, 34]. No mother reported any alcohol abuse or dependence prior to or during pregnancy. A higher proportion of mothers reported consuming some (“any”) alcohol during pregnancy in the unexposed (60%) then exposed (29%) groups, although this difference was not significant (p value = 0.12, Chi-Square, Additional file 4: Table S1). All samples were identified as Caucasians.
Sample dissection and processing
Upon delivery, the products of conception were refrigerated, and within hours, they were moved to a −80 °C freezer. For examination, they were placed at −20 °C overnight. Working quickly over dry ice, the brain was removed without thawing. The cortical plate was sampled in the region that becomes the DLPFC in order to obtain post-migratory NeuN-immunoreactive (NeuN+) neurons, which normally become numerous in cortical layers 4–6 between 14 and 20 weeks gestational age, and in layers 2 and 3 between 20 and 24 weeks (Additional file 1: Figure S1) . This region was chosen because it is involved in decision-making and working memory, and its function is compromised in neurodevelopmental and psychiatric conditions, including autism spectrum disorder (ASD). It is readily identified and accessible in second-trimester fetal brain. During the second trimester, the cerebral hemispheric wall in the frontal region grows from a thickness of ~2 mm at 12 weeks to ~6–8 mm at 18 weeks and ~18 mm at 26 weeks, with cortical plate thickness of ~0.5–1, ~1.5, and ~2 mm, respectively [36–38]. We obtained tissue from the cortical plate from frozen fetal brains by scraping the dorsal prefrontal region of the left hemisphere to a depth of approximately 0.5 mm for the youngest fetuses, where there was no gross demarcation between plate and subplate. In the older fetuses, we were guided by a change in color at the junction of the cortical plate and subplate at approximately the predicted depth. These sample specimens for DNA methylation and gene expression assays were stored at −80 °C for further processing.
Human induced Neuronal Precursor Cells (hiNPC)
All hiNPC lines were derived as previously described . To match the in vivo data generated from postmortem studies, hiNPC lines (NSB553-3-C, NSB2607-4-1 and NSB690-2-1) used in this study were derived from three Caucasian males, and for full details of the donors of the fibroblasts and validation of the hiPSC and NPC lines, please see . Cell culture; NPCs were maintained at high density, grown on growth factor-reduced Matrigel (BD Biosciences)-coated plates in NPC media (Dulbecco’s Modified Eagle Medium/Ham’s F12 Nutrient Mixture (ThermoFisher Scientific), 1× N2, 1× B27-RA (ThermoFisher Scientific) and 20 ng/ml−1 FGF2 and split 1:3 every week with Accutase (Millipore, Billerica, MA, USA). Neural differentiation; NPCs were dissociated with Accutase and plated at 2.0 × 105 cells per cm2 in NPC media onto growth factor-reduced Matrigel-coated plates. For neuronal differentiation, medium was changed to neural differentiation medium (DMEM/F12, 1× N2, 1× B27-RA, 20 ng/ml−1 BDNF (Peprotech), 20 ng/ml−1 GDNF (Peprotech), 1 mM dibutyryl-cyclic AMP (Sigma), 200 nM ascorbic acid (Sigma) 1–2 days later. NPC-derived neurons were differentiated for 3 and 6 weeks before being assayed.
Nicotine (N0267-100MG, Sigma) was diluted at three different concentrations [100 nM (low), 10 μM (med) and 1 mM (high)] in neuronal media and added every 2nd day with a complete media change. Control wells were treated with equal volume of vehicle (ethanol) added to neuronal media.
The cell impermeant nuclei dye TO-PRO3® (ThermoFisher Scientific, T3605) was added at week 3 of differentiation. Three plates, each containing triplicates of a hiNPC line, were imaged with an Odyssey® infrared imaging system (LI-COR). TO-PRO3® fluorescence intensity was normalized to control (vehicle treated) wells.
At 3 and 6 weeks of differentiation, cells were washed once with 1× PBS and then fixed in 4% paraformaldehyde (Electron Microscopy Services) for 15 min. Following 3 washes with 1× PBS, cells were then blocked and permeabilized with 1% v/v BSA Fraction V (BSA, ThermoFisher Scientific) with 0.3% v/v Triton-X 100 (T-100X, Sigma). Primary antibodies (Rb-Ki67, 1:500, Abcam, ab15580 and Ms-TUJ1, 1:1000, Covance, MMS-435P) were added overnight in 1%BSA/0.5%T-100X. Appropriate secondary antibodies (AlexaFluor Dk secondaries, Ms-680 and Rb-800) were incubated for 2.5 h in 1%BSA/0.5%T-100X. Following 3 washes with 1× PBS, plates were imaged with an Odyssey® infrared imaging system (LI-COR). Fluorescence intensity was normalized to control wells. Statistical differences between nicotine-treated and vehicle-treated controls were determined by Student’s ttest using R Language 3.03 .
Illumina Infinium Human Methylation BeadChip sample processing
DNA from fetal brains and hiNPCs were isolated and bisulfite converted (Zymo Research), and CpG methylation was determined using Illumina Infinium Human Methylation BeadChip microarrays (HM450), as described previously .
DNA methylation data preprocessing
The analyses were performed using R Language 3.03  an environment for statistical computing and Bioconductor 2.13 . Raw data files (.idat) were processed by minfi package . All samples displayed a mean probe-wise detection call for the 485,512 array probes <0.0005. The data were normalized, background subtracted and further normalized by SWAN . M values were used in feature selection models. Beta values (logistic transformed M values) were used for sample sex determination and DNA methylation reporting. Probes mapping to multiple locations (N = 19,834), Infinium type I probes with a SNP at the interrogated CpG (N = 13,708) and probes mapping to the X- and Y-chromosomes were removed from analysis (N = 11,648), as described , leaving 452,930 analyzable probes.
DNA methylation analysis
Differentially methylated probes (DMPs) display a mean difference in DNA methylation of at least 20%, corresponding to a methylation difference detectable by the HM450 with 99% confidence . DMPs were mapped to refSeq gene annotations and analyzed using Ingenuity Pathways Analysis (IPA) software (Ingenuity Systems, www.ingenuity.com). Differentially methylated regions (DMR) were found using the bumphunter algorithm applied to DNA methylation M values . Specifically, for each CpG site, we estimate the difference between the M values for the exposed and unexposed adjusting for gestational age, sex and sample chip assignment. An interaction term was included between smoking exposure and gestational age for interaction DMR analysis. The methylation difference estimates are smoothed based on the predefined CpG clusters where the maximal gap between neighboring CpG sites is 500 bp, while the largest cluster size is set to 1500 bp. The smoothed regional methylation difference estimates were obtained using a predefined threshold to identify the putative DMRs, with associated significance levels obtained empirically based on 1000 permutations. Cell-proportion estimates were performed using the methods described in Jaffe et al. . Briefly, publicly available HM450 data from ESC-derived NPC (H9) , adult cortical NeuN+ and NeuN− cells  were quantile normalized together  and 227 unique probes that separated the 3 cell types were used in a nonlinear mixed modeling  to estimate the proportion of each of the 3 cell types within our HM450 fetal dataset. Cell-proportion estimates were also generated for publicly available HM450 data from dissected postnatal DLPFC aged 4, 6 and 10 months, produced by the BrainSpan Consortium .
Gene expression analysis
Total RNA was isolated from the same 24 fetal samples used for DNA methylation analysis (ToTALLY RNA™ Total RNA Isolation Kit, Ambion). Fetal mRNA was analyzed using Nanostring nCounter Elements technology. Gene expression analysis of fetal DLPFC samples was performed for the 2 most significant smoking-DMRs (SDHAP3 and GNA15) and 3 of the 5 genes annotated to the most significant interaction DMRs (C21orf56, POLDIP3 and SYCE3). Housekeeping gene selection: We used the Nanostring nCounter Elements technology and selected 4 housekeeping genes for expression normalization. Previously, Penna et al.  investigated the stability of a panel of housekeeping genes for mRNA normalization in human postmortem brain samples. Additionally, Madden et al.  described a subset of ubiquitously expressed transcripts ideal for using as housekeeping genes within brain tissue. We selected 4 housekeeping genes, 3 of which were identified by both Penna et al. and Madden et al. (GAPDH, YWHAZ and CYC) and SDHA, identified by the former group and that we have previously used successfully in Nanostring interrogation of rat mRNA . Negative control subtraction and normalization to housekeeping genes was performed using the nSolver Analysis Software. Sample fold-changes (FC) were calculated relative to sample FS5777 (one of the 24 fetal samples analyzed chosen at random) gene expression levels for each gene independently. Any expression values of 3 standard deviations from the group mean were deemed outliers and removed from the analysis. No more than 1 result for any assay was removed.
Fetal samples exposed to maternal smoking were matched to unexposed fetal samples by age and sex when available (Table 1). Fetal brain weight and total weight were highly correlated (R 2 = 0.94, Additional file 1: Figure S1a), and within either early or late second-trimester samples, no significant difference in brain weight was observed between exposure groups (p value = 0.3 and 0.8, respectively, Student’s t test) (Additional file 1: Figure S1b, c). The cortical plate was sampled from the presumptive DLPFC in an effort to obtain post-migratory NeuN+ neurons, which normally become numerous in cortical layers 4–6 between 14 and 20 weeks gestational age, and in layers 2 and 3 between 20 and 24 weeks .
Maternal smoking-associated differential DNA methylation in the fetal cortex
To examine DMRs associated with maternal smoking exposure in the fetal cortex, we performed DNA methylation microarray analysis among exposure groups (“Methods” section). Across the 452930 CpG sites examined, no DMRs passed multiple testing correction (family wise error rate, fwer, cutoff ≤0.05, Fig. 1a). This was likely due to the small number of fetal brains available for analysis. Notably, smoking exposure DMRs with the highest point-wise significance was found within the gene promoters of SDHAP3 and G protein subunit Alpha 15 (GNA15) (Fig. 1b, c). Both DMRs were hypomethylated in smoking exposed (Fig. 1b, c). Gene expression analysis of SDHAP3 and GNA15 (Fig. 1d, e) did not reveal exposure-related differences. However, we observed up-regulation of mRNA between early and late second-trimester samples that were restricted to smoking-exposed fetuses for expression of SDHAP3 and larger for exposed than for unexposed for expression of GNA15 (p value = 0.005 vs 0.02, Fig. 1e). These observations led us to hypothesize that maternal smoking exposure has temporal effects on gene expression and DNA methylation during fetal cortical development.
Developmental interaction with maternal smoking exposure
To explore the effects of smoking exposure on fetal development, we performed an interaction DMR analysis between gestational age and smoking exposure (“Methods” section). Four significant DMRs (fwer p value ≤ 0.05) located within the promoters of synaptonemal complex central element protein 3 (SYCE3), chromosome 21 open reading frame 56 (C21orf56/LSS), sperm-associated antigen 1 (SPAG1) and RNA, U12 Small Nuclear (RNU12/POLDIP3) (Fig. 2a) were identified. Within the promoter region of C21orf56, exposed fetal cortices at 22.4 weeks (late) exhibited higher total DNA methylation levels than samples examined at 16.6 weeks (early). Conversely, unexposed fetal cortices displayed higher DNA methylation levels in early samples compared to late (Fig. 2b). Additionally, we found an intragenic region of SPAG1 that was conservatively hypomethylated in all samples except a subset of late second-trimester exposed fetal cortices that were hypermethylated (Additional file 2: Figure S2). Nanostring expression analysis of C21orf56 demonstrated no difference in mRNA levels between exposure/developmental groups (Fig. 2c).
No discernable difference in group-wise DNA methylation patterns by exposure/development was evident in the promoter regions of RUN12/POLDIP3 or SYCE3, regions that had been identified using the bumphunter algorithm as significant (Additional file 3: Figure S3a, b). Gene expression analysis of SYCE3 revealed low gene expression in early second-trimester smoking exposed (p value = 0.02) (Additional file 3: Figure S3c), analogous to the SDHAP3 mRNA results indicative of developmental delay in smoking exposed. No difference in POLDIP3 mRNA was found between exposure/development (Additional file 3: Figure S3d).
Global DNA methylation of fetal cortex exposed to maternal smoking
We examined global DNA methylation changes during cortical development in response to maternal smoking by calculating differentially methylated probes (DMPs) between early and late second-trimester samples in each exposure group (methods). We identified 574/371 hyper/hypomethylated DMPs (Additional file 4: Table S2) in unexposed samples and 399/178 hyper/hypomethylated DMPs (Additional file 4: Table S3) in exposed samples (Fig. 3a). Unsupervised hierarchical clustering of both hypomethylated and hypermethylated DMPs separated samples by gestational age (Fig. 3a). More DMPs were found between early and late second-trimester unexposed samples than between smoking-exposed samples. Notably, a higher proportion of these DMPs were hypermethylated in the unexposed (64%) compared to exposed samples (44%). Hypomethylation is a feature of pluripotent/multipotent stem cells and as cells differentiate, the acquisition of DNA methylation restricts cell lineage and drives cell specification . We postulated that reduced hypermethylated DMPs in smoking-exposed samples could reflect alterations or delay in cell-type specification occurring in the developing cortices in response to smoking exposure. Gene ontology analysis implicated canonical pathways of “Role of NFAT in Cardiac Hypertrophy,” “Th2 Pathway” and “Th1 and Th2 Activation Pathway” associated with genes annotated to DMPs between early and late second-trimester unexposed fetal samples (Additional file 4: Table S4). NFAT is a transcription factor that mediates axon growth in developing neurons (reviewed in ). Conversely, canonical pathways associated with genes annotated to DMPs between early and late second-trimester exposed fetal samples implicated “T Helper Cell Differentiation,” “Epithelial Adherens Junction Signaling” and “Factors Promoting Cardiogenesis in Vertebrates” (Additional file 4: Table S4).
Our cortical plate sectioning of the presumptive DLPFC aimed to enrich for post-migratory NeuN+ neurons (Fig. 3b). Other investigators have established cell deconvolution algorithms with remarkable accuracy in estimating NeuN+ proportions of whole brain tissues using DNA methylation profiles [51, 60]. Using the DNA methylation profiles generated from the fetal cortices, we were able to estimate the cell proportions (CP estimates) of NeuN+, NeuN− and neural precursor cells (NPCs) within the fetal DLPFC sections (“Methods” section). CP estimates revealed our fetal cortical sections contained a high proportion of NPCs (mean = 16%) compared to CP estimates of NPCs of cortical sections (postnatal 4–10 months) produced by the BrainSpan consortium (mean = 7%) [50, 53, 57] (Fig. 3c). These NPCs are presumably undergoing maturation, and thus, our sections provide a rare window into the effects of maternal smoking exposure on neuro-cellular development.
CP estimates of fetal DLPFC revealed significantly fewer NeuN+ in smoking-exposed fetal DLPFC (p value = 0.04, Fig. 4). Any decrease in CP estimate will be balanced by an increase in the proportion of another cellular population, and indeed, we observed a higher proportion of NeuN− in smoking-exposed fetal samples (p value = 0.08, Fig. 4). No difference in NPC proportions was observed between smoking exposure groups (Fig. 4). The results indicate that tobacco exposure is associated with a reduction in NeuN+ in the developing DLPFC. Using ReFACTor  method that adjusts for the possible effect of variation in cell-type proportions, we found no difference in methylation associated with smoking exposure or developmental stage interaction. The methylation change associated with smoking exposure reflects differences in rate of cell differentiation during development as observed by the estimated CP changes. To investigate the mechanism of nicotine exposure on neuron development, we turned to an in vitro model of development.
In vitro modeling of neurodevelopment in response to nicotine exposure
To model the effects of nicotine on human neuronal development, we exposed human induced neural precursor cells (hiNPCs) to nicotine over a 6-week period of neuronal differentiation. Following 3 weeks of nicotine exposure, differentiating hiNPCs exhibited a dosage-dependent increase in TO-PRO®-3 fluorescence (live cells are impermeable to the dye, whereas the dye penetrates compromised membranes characteristic of dead cells) (p value <0.05) and an increase in cell proliferation (Ki67) that was significant at high dose (p value = 0.01, Fig. 5). We observed a decrease in the staining of βIII Tubulin (TUJ1), a marker for mature neurons, and notably, the most significant decrease was found at 0.1 uM (low) nicotine exposure (p value = 0.0003, Fig. 5). No significant differences in cell death or proliferation were detectable by 6 weeks; however, TUJ1 showed a dosage-dependent increase that was significant at 1000 μM nicotine (p value = 0.03, Fig. 5). These data demonstrate that nicotine elicits dosage-dependent effects on neuronal maturation. Interestingly, the biologically meaningful nicotine exposure (0.1 μM) restricts early neuron development despite exhibiting a lower toxicity (cell death) than higher exposures. Further, the differences in TUJ1 staining between 3 and 6 weeks indicate that less differentiated hiNPCs (3 weeks) are more susceptible to the inhibiting effects of nicotine on neuron development (Fig. 5).
In this study, we profiled genome-scale DNA methylation of fetal brain development in response to exposure to maternal smoking in utero. Although limited in number, the fetal brain samples have well-characterized maternal health history and exposure data, thus providing a rare opportunity to investigate the impact of nicotine exposure on early human brain development. Notably, our DNA methylation profiling was performed on whole tissue sections from the developing DLPFC that consist of a mixture of neuronal and non-neuronal cell types. Established methods enable the isolation of neuronal nuclei ; however, due to the fragility of fetal neuronal nuclei, this technique cannot be applied. Estimating cell proportions using DNA methylation profiles revealed our fetal DLPFC sections contained approximately 16% NPCs that are presumably undergoing maturation, providing a rare window into human neuro-cellular development in response to maternal smoking exposure.
Embryogenesis is a stage of rapid neurological transformation and growth in which epigenomic landscapes undergo dramatic change [63, 64]. In the developing DLPFC, genome-wide DNA methylation changes that distinguish early and late second-trimester samples were clearly reduced in fetal cortex of smoking exposed indicating alterations in cell-type differentiation. Indeed, our most significant DMRs were identified by interaction analysis between gestational age and maternal smoking exposure. We found a delay in the up-regulation of expression of SYCE3, a gene that is conserved among mammals and whose loss leads to a block in synapsis initiation resulting in meiotic arrest . We also identified an interaction DMR in the promoter of C21orf56 (also known as SPATC1L). C21orf56 is a spermatogenesis and centriole-associated 1-like gene found on chromosome 21 of little-known function. Although we did not observe developmental/exposure dependent changes in gene expression, we observed a gain in promoter methylation in late second-trimester exposed samples that could reduce the transcriptional potential of C21orf56 later in development.
Genes such as GNA15 and SDHAP3 that contained maternal smoking-associated DMRs displayed a developmental delay in mRNA up-regulation in smoking exposed. Notably, SDHAP3 is a subunit of the succinate dehydrogenase complex located within the mitochondrial membrane and functions in electron transport chain transfer of electrons to coenzyme Q . It has been reported that mutations within succinate dehydrogenase subunits actually increase levels of oxidative stress . Intriguingly, this same DMR was recently found hypermethylated in the cerebellum of patients diagnosed with ASD  and differentially methylated in the DLPFC of patients diagnosed with schizophrenia (SCZ) within 3 independent studies . Furthermore, in a separate report, GNA15 was found to be differentially methylated in the PFC of ASD patients . Both ASD and SCZ probably have prenatal origins [71–73]. Taken together, these results reveal a potential link between maternal smoking-associated DNA methylation perturbation and potential increase risk for neurodevelopmental abnormalities. Notably, GNA15 is transcriptionally modifiable by acute doses of nicotine in neuroblastoma cell lines , indicating nicotine as a potential causative agent.
Cell deconvolution algorithms have shown remarkable accuracy in estimating NeuN+ proportions from DNA methylation profiles from whole brain tissue [51, 60]. CP estimates within our fetal DLPFC revealed a smoking exposure-associated reduction in NeuN+ cells supporting previous observations of reduced gray matter in the cortex of smoking-exposed children . The adverse effects of maternal smoking on fetal development are well described; however, it was estimated that 30% of smokers attempting to quit smoking use cessation aids that contain nicotine . Nicotine is a well-studied substance in tobacco and has been shown to induce oxidative stress in rodent [76, 77] and human neurons . Exposure of hiNPCs to 100 nM nicotine resulted in the lowest amount of toxicity but the greatest suppression of neuronal differentiation (B3-Tubulin). These results recapitulate the reduction in the estimated proportion of NeuN+ cells we observed in human samples and implicate nicotine as a causative agent in impeding neuronal development. These data provide direct evidence from primary tissue of in utero exposure to teratogenic agents as found in cigarettes—warranting further investigations of the in utero environment on fetal development and how it impacts offspring health and disease risk through the lifespan.
In summary, we have found evidence that intrauterine smoking exposure alters the developmental patterning of DNA methylation and gene expression and is associated with reduced mature neuronal content, effects that are likely driven by nicotine through mechanisms independent of cell death.
autism spectrum disorder
brain-derived neurotrophic factor-6
chromosome 21 open reading frame 56
dorsolateral prefrontal cortex
differentially methylated positions
differentially methylated region
false discovery rate
- GNA15 :
G protein subunit alpha 15
human induced pluripotent stem cell
Illumina Infinium Human Methylation Microarray 450K platform
monoamine oxidase A
nicotinic acetylcholine receptors
- NeuN+ :
neural progenitor cell
RNA, U12 Small Nuclear
sperm-associated antigen 1
Cornelius MD, et al. Effects of prenatal cigarette smoke exposure on neurobehavioral outcomes in 10-year-old children of adolescent mothers. Neurotoxicol Teratol. 2011;33(1):137–44.
Huijbregts SC, et al. Hot and cool forms of inhibitory control and externalizing behavior in children of mothers who smoked during pregnancy: an exploratory study. J Abnorm Child Psychol. 2008;36(3):323–33.
Robinson M, et al. Smoking cessation in pregnancy and the risk of child behavioural problems: a longitudinal prospective cohort study. J Epidemiol Community Health. 2010;64(7):622–9.
Roza SJ, et al. Effects of maternal smoking in pregnancy on prenatal brain development. The Generation R Study. Eur J Neurosci. 2007;25(3):611–7.
El Marroun H, et al. Prenatal tobacco exposure and brain morphology: a prospective study in young children. Neuropsychopharmacology. 2014;39(4):792–800.
Eskenazi B, Bergmann JJ. Passive and active maternal smoking during pregnancy, as measured by serum cotinine, and postnatal smoke exposure. I. Effects on physical growth at age 5 years. Am J Epidemiol. 1995;142(9 Suppl):S10-8.
Wakschlag LS, et al. Maternal smoking during pregnancy and severe antisocial behavior in offspring: a review. Am J Public Health. 2002;92(6):966–74.
Talhout R, et al. Hazardous compounds in tobacco smoke. Int J Environ Res Public Health. 2011;8(2):613–28.
Li Y, et al. Perinatal nicotine exposure increases vulnerability of hypoxic–ischemic brain injury in neonatal rats: role of angiotensin II receptors. Stroke J Cereb Circ. 2012;43(9):2483–90.
Mychasiuk R, et al. Long-term alterations to dendritic morphology and spine density associated with prenatal exposure to nicotine. Brain Res. 2013;1499:53–60.
Slotkin TA, et al. Effects of prenatal nicotine exposure on primate brain development and attempted amelioration with supplemental choline or vitamin C: neurotransmitter receptors, cell signaling and cell development biomarkers in fetal brain regions of rhesus monkeys. Neuropsychopharmacology. 2005;30(1):129–44.
Bruin JE, Gerstein HC, Holloway AC. Long-term consequences of fetal and neonatal nicotine exposure: a critical review. Toxicol Sci Off J Soc Toxicol. 2010;116(2):364–74.
Dwyer JB, McQuown SC, Leslie FM. The dynamic effects of nicotine on the developing brain. Pharmacol Ther. 2009;122(2):125–39.
Button TM, Maughan B, McGuffin P. The relationship of maternal smoking to psychological problems in the offspring. Early Human Dev. 2007;83(11):727–32.
Lim R, Sobey CG. Maternal nicotine exposure and fetal programming of vascular oxidative stress in adult offspring. Br J Pharmacol. 2011;164(5):1397–9.
Mahar I, et al. Developmental hippocampal neuroplasticity in a model of nicotine replacement therapy during pregnancy and breastfeeding. PLoS ONE. 2012;7(5):e37219.
Pauly JR, Slotkin TA. Maternal tobacco smoking, nicotine replacement and neurobehavioural development. Acta Paediatr. 2008;97(10):1331–7.
Jones KL, Smith DW. Recognition of the fetal alcohol syndrome in early infancy. Lancet. 1973;302(7836):999–1001.
Grandjean P, Landrigan PJ. Developmental neurotoxicity of industrial chemicals. Lancet. 2006;368(9553):2167–78.
Julvez J, Grandjean P. Neurodevelopmental toxicity risks due to occupational exposure to industrial chemicals during pregnancy. Ind Health. 2009;47(5):459–68.
Gidekel S, Bergman Y. A unique developmental pattern of Oct-3/4 DNA methylation is controlled by a cis-demodification element. J Biol Chem. 2002;277(37):34521–30.
Adcock IM, et al. Epigenetic regulation of airway inflammation. Curr Opin Immunol. 2007;19(6):694–700.
Furniss CS, et al. Line region hypomethylation is associated with lifestyle and differs by human papillomavirus status in head and neck squamous cell carcinomas. Cancer Epidemiol Biomarkers Prev. 2008;17(4):966–71.
Kanai Y. Alterations of DNA methylation and clinicopathological diversity of human cancers. Pathol Int. 2008;58(9):544–58.
Lawrence J, et al. Foetal nicotine exposure causes PKCepsilon gene repression by promoter methylation in rat hearts. Cardiovasc Res. 2011;89(1):89–97.
Prescott SL, Clifton V. Asthma and pregnancy: emerging evidence of epigenetic interactions in utero. Curr Opin Allergy Clin Immunol. 2009;9(5):417–26.
Satta R, et al. Nicotine decreases DNA methyltransferase 1 expression and glutamic acid decarboxylase 67 promoter methylation in GABAergic interneurons. Proc Natl Acad Sci USA. 2008;105(42):16356–61.
Zochbauer-Muller S, et al. Aberrant methylation of multiple genes in the upper aerodigestive tract epithelium of heavy smokers. Int J Cancer. 2003;107(4):612–6.
Joubert BR, et al. DNA methylation in newborns and maternal smoking in pregnancy: genome-wide consortium meta-analysis. Am J Hum Genet. 2016;98(4):680–96.
Ladd-Acosta C, et al. Presence of an epigenetic signature of prenatal cigarette smoke exposure in childhood. Environ Res. 2016;144(Pt A):139–48.
Philibert RA, et al. MAOA methylation is associated with nicotine and alcohol dependence in women. Am J Med Genet B Neuropsychiatr Genet. 2008;147B(5):565–70.
Xu Q, et al. Determination of methylated CpG sites in the promoter region of catechol-O-methyltransferase (COMT) and their involvement in the etiology of tobacco smoking. Front Psychiatry. 2010;1:16.
Toledo-Rodriguez M, et al. Maternal smoking during pregnancy is associated with epigenetic modifications of the Brain-derived neurotrophic factor-6 exon in adolescent offspring. Am J Med Genet B Neuropsychiatr Genet. 2010;153B(7):1350–4.
Roebuck TM, Mattson SN, Riley EP. A review of the neuroanatomical findings in children with fetal alcohol syndrome or prenatal exposure to alcohol. Alcohol Clin Exp Res. 1998;22(2):339–44.
Sarnat HB, Nochlin D, Born DE. Neuronal nuclear antigen (NeuN): a marker of neuronal maturation in early human fetal nervous system. Brain Dev. 1998;20(2):88–94.
Kostovic I, et al. Laminar organization of the human fetal cerebrum revealed by histochemical markers and magnetic resonance imaging. Cereb Cortex. 2002;12(5):536–44.
Rados M, Judas M, Kostovic I. In vitro MRI of brain development. Eur J Radiol. 2006;57(2):187–98.
Shepard TH, et al. Organ weight standards for human fetuses. Pediatr Pathol. 1988;8(5):513–24.
Lee IS, et al. Characterization of molecular and cellular phenotypes associated with a heterozygous CNTNAP2 deletion using patient-derived hiPSC neural cells. Npj Schizophr. 2015;1:15019.
Topol A, et al. Dysregulation of miRNA-9 in a subset of schizophrenia patient-derived neural progenitor cells. Cell Rep. 2016;15(5):1024–36.
R Development Core Team. R: a language and environment for statistical computing. R Foundation for Statistical Computing. 2014; http://www.R-project.org/.
Bibikova M, et al. High-throughput DNA methylation profiling using universal bead arrays. Genome Res. 2006;16(3):383–93.
Gentleman RC, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5(10):R80.
Aryee MJ, et al. Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics. 2014;30(10):1363–9.
Maksimovic J, Gordon L, Oshlack A. SWAN: subset-quantile within array normalization for illumina infinium HumanMethylation450 BeadChips. Genome Biol. 2012;13(6):R44.
Naeem H, et al. Reducing the risk of false discovery enabling identification of biologically significant genome-wide methylation status using the HumanMethylation450 array. BMC Genom. 2014;15:51.
Bibikova M, et al. High density DNA methylation array with single CpG site resolution. Genomics. 2011;98(4):288–95.
Jaffe AE, et al. Bump hunting to identify differentially methylated regions in epigenetic epidemiology studies. Int J Epidemiol. 2012;41(1):200–9.
Jaffe AE, et al. Developmental regulation of human cortex transcription and its clinical relevance at single base resolution. Nat Neurosci. 2015;18(1):154–61.
Kim M, et al. Dynamic changes in DNA methylation and hydroxymethylation when hES cells undergo differentiation toward a neuronal lineage. Hum Mol Genet. 2014;23(3):657–67.
Guintivano J, Aryee MJ, Kaminsky ZA. A cell epigenotype specific model for the correction of brain cellular heterogeneity bias and its application to age, brain region and major depression. Epigenetics. 2013;8(3):290–302.
Houseman EA, et al. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinform. 2012;13:86.
BrainSpan: Atlas of the developing human brain [Internet]. Funded by ARRA Awards 1RC2MH089921-01, R.M.-., and 1RC2MH089929-01. © 2011. http://developinghumanbrain.org.
Penna I, et al. Selection of candidate housekeeping genes for normalization in human postmortem brain samples. Int J Mol Sci. 2011;12(9):5461–70.
Velculescu VE, et al. Analysis of human transcriptomes. Nat Genet. 1999;23(4):387–8.
Haghighi F, et al. Neuronal DNA methylation profiling of blast-related traumatic brain injury. J Neurotrauma. 2015;32(16):1200–9.
Sasaki T, et al. Type of feeding during infancy and later development of schizophrenia. Schizophr Res. 2000;42:79–82.
Mohn F, et al. Lineage-specific polycomb targets and de novo DNA methylation define restriction and potential of neuronal progenitors. Mol Cell. 2008;30(6):755–66.
Nguyen T, Di Giovanni S. NFAT signaling in neural development and axon growth. Int J Dev Neurosci. 2008;26(2):141–5.
Montano CM, et al. Measuring cell-type specific differential methylation in human brain tissue. Genome Biol. 2013;14(8):R94.
Rahmani E, et al. Sparse PCA corrects for cell type heterogeneity in epigenome-wide association studies. Nat Methods. 2016;13(5):443–5.
Jiang Y, et al. Isolation of neuronal chromatin from brain tissue. BMC Neurosci. 2008;9:42.
Lister R, et al. Global epigenomic reconfiguration during mammalian brain development. Science. 2013;341(6146):1237905.
Numata S, et al. DNA methylation signatures in development and aging of the human prefrontal cortex. Am J Hum Genet. 2012;90(2):260–72.
Schramm S, et al. A novel mouse synaptonemal complex protein is essential for loading of central element proteins, recombination, and fertility. PLoS Genet. 2011;7(5):e1002088.
Yankovskaya V, et al. Architecture of succinate dehydrogenase and reactive oxygen species generation. Science. 2003;299(5607):700–4.
Slane BG, et al. Mutation of succinate dehydrogenase subunit C results in increased O2.–, oxidative stress, and genomic instability. Cancer Res. 2006;66(15):7615–20.
Ladd-Acosta C, et al. Common DNA methylation alterations in multiple brain regions in autism. Mol Psychiatry. 2014;19(8):862–71.
Wockner LF, et al. Brain-specific epigenetic markers of schizophrenia. Transl Psychiatry. 2015;5:e680.
Nardone S, et al. DNA methylation analysis of the autistic brain reveals multiple dysregulated biological pathways. Transl Psychiatry. 2014;4:e433.
Christensen J, et al. Prenatal valproate exposure and risk of autism spectrum disorders and childhood autism. JAMA. 2013;309(16):1696–703.
Pidsley R, et al. Methylomic profiling of human brain tissue supports a neurodevelopmental origin for schizophrenia. Genome Biol. 2014;15(10):483.
Jaffe AE, et al. Mapping DNA methylation across development, genotype and schizophrenia in the human frontal cortex. Nat Neurosci. 2016;19(1):40–7.
Wang J, et al. Genome-wide expression analysis reveals diverse effects of acute nicotine exposure on neuronal function-related genes and pathways. Front Psychiatry. 2011;2:5.
Census U.D.o.C.B. Current population survey, January 2007: Tobacco use supplement file: Tech. Doc. CPS-07. Washington, DC: US Department of Commerce; 2007.
Bhagwat SV, et al. Preferential effects of nicotine and 4-(N-methyl-N-nitrosamine)-1-(3-pyridyl)-1-butanone on mitochondrial glutathione S-transferase A4-4 induction and increased oxidative stress in the rat brain. Biochem Pharmacol. 1998;56(7):831–9.
Helen A, et al. Antioxidant effect of onion oil (Allium cepa. Linn) on the damages induced by nicotine in rats as compared to alpha-tocopherol. Toxicol Lett. 2000;116(1–2):61–8.
Guan ZZ, Yu WF, Nordberg A. Dual effects of nicotine on oxidative stress and neuroprotection in PC12 cells. Neurochem Int. 2003;43(3):243–9.
ZC, NM, SC, MM, YG and FH designed and performed DNA methylation and gene expression experiments and analysis. ZC, FH, BH, KB and SMH designed and performed cell culture studies. AD, GR, AS and ITI designed and collected mothers’ history and fetal brain samples. AD performed fetal brain dissections. ZC, FH, AD, YG, KB, MM and BH wrote and edited the manuscript. All authors read and approved the final manuscript.
Please see title page for author affiliations and contact information.
This work was supported in part through the computational resources and staff expertise provided by Scientific Computing at the Icahn School of Medicine at Mount Sinai.
The authors declare that they have no competing interests.
Availability of data and materials
All raw Illumina HM450 DNA methylation data used in this study have been deposited in GEO (GSE90871).
Consent for publication
All authors have read and consent to the publication of this research article.
Ethical approval and consent to participate
This study was approved by the Institutional Review Board at the Icahn School of Medicine at Mount Sinai. Informed consent was obtained by a psychologist with no involvement in the mothers’ clinical care, who also interviewed them to obtain a medical history, family history and history of environmental exposures.
The Haghighi Laboratory is supported by the National Institute of Health (NIH) Grant R01MH094774. ZC is supported by a NIDA T32 training grant in Drug Abuse Research from the NIH, USA. Kristen J. Brennand is a New York Stem Cell Foundation—Robertson Investigator. The Brennand Laboratory is supported by the Brain and Behavior Research Foundation, NIH Grants R01 MH101454 and R01 MH106056, and the New York Stem Cell Foundation. Research reported in this paper was supported by the Office of Research Infrastructure of the National Institutes of Health under award number S10OD018522. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.