- Open Access
Aberrant methylation of Pax3 gene and neural tube defects in association with exposure to polycyclic aromatic hydrocarbons
Clinical Epigeneticsvolume 11, Article number: 13 (2019)
Neural tube defects (NTDs) are common and severe congenital malformations. Pax3 is an essential gene for neural tube closure in mice but it is unknown whether altered expression or methylation of PAX3 contributes to human NTDs. We examined the potential role of hypermethylation of Pax3 in the development of NTDs by analyzing human NTD cases and a mouse model in which NTDs were induced by benzo[a]pyrene (BaP), a widely studied polycyclic aromatic hydrocarbon (PAH).
We extracted methylation information of PAX3 in neural tissues from array data of ten NTD cases and eight non-malformed controls. A validation study was then performed in a larger independent population comprising 73 NTD cases and 29 controls. Finally, we examined methylation patterns and expression of Pax3 in neural tissues from mouse embryos of dams exposed to BaP or BaP and vitamin E.
Seven CpG sites in PAX3 were hypermethylated in NTD fetuses as compared to controls in the array data. In the validation phase, significantly higher methylation levels in the body region of PAX3 were observed in NTD cases than in controls (P = 0.003). And mean methylation intensity in the body region of PAX3 in fetal neural tissues was positively correlated with median concentrations of PAH in maternal serum. In the mouse model, BaP-induced NTDs were associated with hypermethylation of specific CpG sites within both the promoter and body region of Pax3. Supplementation with vitamin E via chow decreased the rate of NTDs, partly recovered the repressed total antioxidant capacity in mouse embryos exposed to BaP, and this was accompanied by the normalization of Pax3 methylation level and gene expression.
Hypermethylation of Pax3 may play a role in the development of NTDs; DNA methylation aberration may be caused by exposure to BaP, with possible involvement of oxidative stress.
Neural tube defects (NTDs) arise from a failed or disordered closure of the neural tube during embryogenesis. The occurrence of NTDs is around 0.5–2/1000 pregnancies worldwide . Fetuses affected with an NTD are often stillborn (e.g., anencepahlics), and most surviving infants suffer from life-long disabilities. The etiology of NTDs is complex, involving both genetic and non-genetic factors . Over 300 genes have been identified to be involved in the regulation of neural tube closure in mouse NTD mutants [3, 4]; however, only a few of these genes have successfully been validated in human NTDs. In addition to potential additive effects of multiple risk alleles, increasing attention has also focused on the potential for epigenetic alterations to contribute to the occurrence of NTDs by mediating the interplay of fetal genetics and environmental factors .
Epigenetic modification can cause changes in gene expression that are not directly related to the DNA sequence itself, of which DNA methylation is one of the best understood epigenetic mechanisms . During early development, a tight regulation of genome-wide erasure of epigenetic footprints with resetting of the methylation signature takes place, making developing fetuses particularly susceptible to epigenetic dysregulation as a consequence of environmental exposure [7, 8]. Abnormal genome-wide methylation during embryogenesis has been linked to developmental abnormalities at birth, including NTDs. Recent studies have demonstrated that global DNA hypomethylation, evaluated using LINE-1 methylation as an indicator in human fetuses, was associated with an increased risk for NTDs . Aberrations of methylation at specific genes are reported to be involved in NTDs, including imprinted genes [10, 11], DNA repair genes , planar-cell polarity genes [13, 14], and HOX genes .
Pax3 is a paired-homeodomain-containing transcription factor essential for promoting neural crest induction, maintenance, migration, and differentiation . Previous studies in mice have shown that Pax3 function is required for neural tube closure. Several alleles of Pax3 cause NTDs in mice and homozygous Splotch (Sp2H) embryos develop NTDs with 100% penetrance . However, the role of PAX3 in human NTDs remains unclear. Recessive and dominant mutations in PAX3 in humans are known to cause Waardenburg syndrome, an autosomal dominant condition that affects neural crest-derived structures and also includes spina bifida as part of its phenotypic spectrum . A 5-bp deletion in exon 5 of the PAX3 gene was reported in a patient with spina bifida . In contrast, the results from a case-control study including 74 infants with spina bifida and 87 non-malformed controls indicated that variants in PAX3 were not strong risk factors for human spina bifida . Therefore, exploring mechanisms other than coding sequence variants, such as methylation modification, in PAX3 may provide novel insight into the etiology of human NTDs.
We hypothesized that aberrant DNA methylation of Pax3 plays a role in the formation of NTDs. To test this hypothesis, we first compared the difference in methylation levels of CpG sites within PAX3 using genomic methylation array data with DNA from neural tissues of NTD cases and non-malformed controls. The methylation status of the CpG sites that were found to differ was then validated in a larger NTD case-control population. We tested whether there was any correlation between CpG site methylation levels in fetal neural tissues and maternal serum concentrations of polycyclic aromatic hydrocarbons (PAHs), a class of ubiquitous environmental pollutants that have been shown to be associated with the risk for NTDs in epidemiological studies . Finally, we assessed the methylation level and gene expression of Pax3 in neural tissues from mouse embryos exposed in utero to benzo[a]pyrene (BaP), a widely studied PAH that induces NTDs . In addition, we assayed markers of oxidative stress in fetal mice to further explore the possible mechanisms by which BaP might affect methylation regulation.
Methylation of PAX3 gene in genomic microarray
A detailed description on the genome-wide methylation results obtained from Infinium HumanMethylation450 BeadChip (HM450K), using DNA isolated from neural tissues from ten NTD cases and eight unrelated non-malformed controls can be found elsewhere . In brief, out of 485,199 CpG sites across the entire genome, 23,294 (4.8%) were differentially methylated between cases and controls. Of the differentially methylated CpG sites, 12,383 (53.2%) were significantly hypermethylated and 10,911 (46.8%) were significantly hypomethylated in NTD cases when compared to controls. The characteristics of the NTD cases and controls in phase 1 are presented in Additional file 1: Table S1.
In the PAX3 gene, the focus of this study, a total of 54 CpG sites were extracted from the HM450K array data, of which 47 CpG sites (87.0%) were found to be hypermethylated. And 7 out of the 47 CpG sites exhibited statistically greater methylation in NTD fetuses than in controls (Additional file 2: Table S2). Analysis of the genomic location of the seven significantly hypermethylated CpG sites showed that one CpG was located at TSS1500, one in the 5′UTR, and the remaining five CpGs were within the body of the gene (Fig. 1a).
Validation of PAX3 methylation in an independent cohort of NTD cases and controls
In order to validate these findings, 73 NTD cases and 29 controls were used to examine the differentially methylated regions of PAX3 identified in phase 1 using Sequenom EpiTYPER. The demographic characteristics of the subjects are described in Table 1. Three DNA amplicons were developed that cover 35 CpGs (Fig. 1a). In the present study, TSS1500, TSS200, 5′UTR, and the 1st exon were defined as the promoter region of PAX3 gene. As shown in Fig. 2a, no significant difference in overall mean methylation levels in the promoter region of PAX3 was observed between NTD cases and controls. However, NTD samples exhibited significantly higher DNA methylation levels (13.9 ± 10.1%) than control samples (7.8 ± 4.2%) in the gene body region (P = 0.003) (Fig. 2b). Furthermore, within the gene body, all CpG sites except CpG_12 showed significantly higher levels of methylation in the NTD samples compared to controls (Fig. 2e). After Benjanmini-Hochberg correction, methylation levels in five of the nine CpG sites were still significantly higher in NTD cases.
The relationship between hypermethylation of PAX3 and the risk of NTDs was examined. Methylation levels were categorized according to the median methylation level of the controls. As shown in Table 2, a higher level of methylation in the gene body region was associated with 6.24-fold increased risk for NTDs (95% CI 1.30–29.97).
Correlation between PAH concentrations in maternal serum and PAX3 methylation in fetal neural tissues
A previous study from our team found that higher PAH concentrations in maternal serum were associated with an increased risk for NTDs . We therefore conducted a correlation analysis between PAH concentrations in maternal serum and PAX3 methylation levels in fetal neural tissues (N = 51 mother-fetus pairs) to further our understanding of the potential relationship of PAH exposure, gene methylation, and NTDs risk. Notably, a significant positive correlation was found between mean methylation levels in the PAX3 gene body and median concentrations of high-molecular-weight PAHs in maternal serum (r = 0.310. P = 0.027) (Additional file 3: Table S3).
Disturbed methylation of Pax3 in mouse embryos exposed to BaP
We further investigated the finding of a correlation between maternal serum PAH concentrations and PAX3 hypermethylation in humans by utilizing a BaP-induced NTD mouse model. Three amplicons developed according to the mouse genome were used (Fig. 1b). The body region of Pax3 showed a trend toward higher levels of methylation in the BaP-treated group compared with controls (Fig. 3a), but this did not reach statistical significance. Furthermore, in terms of specific CpG sites, two significantly hypermethylated CpG sites within the promoter and two within the body region of Pax3 were respectively detected in the BaP-exposed mouse fetuses when compared to the non-exposed fetuses (Fig. 3b–d).
The role of oxidative stress in BaP-related disturbed Pax3 methylation and expression
Since BaP exposure was hypothesized to favor the generation of reactive oxygen species (ROS) in the embryos, we inferred that oxidative stress may get involved in the aberrant Pax3 methylation induced by BaP. As shown in Fig. 4, compared to the control group, the levels of total antioxidant capacity (TAC) and the activity of superoxide dismutases (SODs) were decreased. Real-time PCR and whole mount in situ hybridization showed that BaP exposure significantly decreased the transcription of Pax3 (P = 0.008) (Fig. 5). After co-supplemented with vitamin E, a commonly used antioxidant, the repressed levels of TAC and SODs were partly restored (Fig. 4). Coincidentally, the hypermethylation of Pax3 in specific CpG sites induced by BaP treatment was recovered and Pax3 expression was also normalized by vitamin E supplementation (Figs. 3 and 5), which suggested a causal effect of oxidative stress on both Pax3 methylation and expression.
Pax3 is a key gene, encoding a transcription factor required for neural tube closure . In this study, we evaluated the methylation status of PAX3 in a two-phase design study to examine whether there is a potential role of PAX3 methylation in the development of human NTDs. Differentially hypermethylated CpG sites were found in the promoter and gene body region within PAX3 in the neural tissues of NTD cases. Association analysis showed that a higher methylation level in the gene body region was associated with an elevated risk for NTDs. Moreover, PAH concentrations in maternal serum, which were known to be associated with increased risk of NTDs, were positively correlated to methylation levels at several CpG sites. In the BaP-induced NTD mouse model, hypermethylation of Pax3 gene and suppressed gene expression were observed in embryos with BaP treatment, along with reduced TAC level. We previously showed that vitamin E supplementation could decrease the rate of NTDs and alleviate oxidative stress in BaP exposed embryos . Here, we further demonstrated that vitamin E mitigated the shifts in Pax3 methylation and gene expression.
Studies on the role of the Pax3 in NTDs have mostly focused on gene depletion or loss of function mutants . In mice, mutations in Pax3 give rise to the Splotch (Sp) phenotype, which includes exencephaly, spina bifida, and neural crest abnormalities in homozygous mutant embryos [25, 26]. The human PAX3 gene exhibits high homology with mouse Pax3. Mutations within the human PAX3 gene have been found in Waardenburg syndromes, a condition which is occasionally associated with NTDs . A previous study identified two spina bifida patients who had small interstitial chromosomal deletions involving PAX3 . Exon sequencing of PAX3 in 114 cases with spina bifida also identified two common variants; however, without unaffected individuals, the influence of these variants on the risk of spina bifida could not be determined . To the contrary, the screening results from 74 spina bifida cases and 87 control infants suggested that variant in PAX3 was not a major contributor to the overall burden of NTDs at population level . Recently, evidence from animal studies suggested that aberrant methylation of Pax3 was involved in the development of NTDs induced by hyperglycemia [29, 30]. However, to date, no study investigating the role of PAX3 methylation in the etiology of human NTDs has been reported. In our study, we analyzed the methylation status of PAX3 gene in neural tissues from human fetuses, and our results showed that the methylation level in PAX3 was significantly higher in NTD cases than in non-malformed controls.
Unexpectedly, in our study, compared with promoter region, the body region of PAX3 was more consistently hypermethylated in NTD cases and the differentially methylated region was mainly located in intron 4. Our results imply that a methylation change in the body region of PAX3 may be an epigenetic component of human NTDs, which is reminiscent with findings for the HOXB7 gene in a myelomeningocele (spina bifida) case-control study . In the latter study, a genome-wide methylation assay found three CpGs in HOXB7 gene body to be hypomethylated in myelomeningocele patients when compared to controls, and these were further verified in a larger population, using the Sequenom EpiTYPER platform. Studies on GRHL3 and SOX18 also reported that altered methylation within the gene body was associated with the risk for NTDs [14, 31].
It is generally accepted that DNA methylation changes in the promoter and gene body regions may have differential effects on gene expression. DNA methylation in promoter regions is usually negatively associated with gene expression, whereas in the gene body, the impact is not consistent, with both positive and negative impacts having been reported [32, 33]. In the present study, a negative correlation between gene body DNA methylation and transcription level of Pax3 was observed in mouse embryos, which was in line with the finding in hyperglycemia-induced NTDs . The underlying mechanisms of gene body methylation in regulating gene expression have not been well understood. Recently, increasing evidence demonstrating the role of DNA methylation in alternative splicing regulation has been reported, which is essential for providing tissue-specific features for some genes . DNA methylation is normally more abundant in exons compared to the flanking introns, which is a marker for distinguishing exons from introns. Previous studies proposed that the change of exon methylation levels would affect the recognition of exons while splicing . It is thus reasonable to hypothesize that the methylation status of intron might also be crucial for alternative splicing. Further, cis-acting element has been recognized in intron, which could regulate gene expression . Within intron 4 of PAX3, a number of transcription factors were predicted to bind at several CpG sites examined in the present study, with high predictive values (predictive value > 8, JASPAR, Fig. 1c), which might be responsible for the observed repressed gene expression. However, more research is necessary to examine this point in detail.
PAX3 hypermethylation by itself is not likely to be a sole cause for NTDs but rather be part of a complex combination of environmental and epigenetic risk factors. Our previous epidemiological studies have suggested that maternal exposure to PAHs was associated with an elevated risk of NTDs in the offspring [21, 24]. Cohort studies reported that prenatal PAH exposure was associated with lower global methylation and hypermethylation of interferon γ in umbilical cord white blood cells [36, 37]. Several toxicological studies have also suggested that BaP exposure could disrupt DNA methylation status [38, 39]. Importantly, in the present study, we found that differentially methylated CpG sites in PAX3 gene in fetal neural tissues were positively correlated with PAH concentrations in maternal serum. In line with these findings in human subjects, our mouse experiment showed that the methylation level of Pax3 was elevated after BaP treatment, indicating that the presence of BaP impacts Pax3 methylation regulation. These findings support the hypothesis that hypermethylation of Pax3 gene is involved in abnormal closure of neural tube secondary to PAH exposure.
In considering how BaP exposure could affect methylation, accumulating evidence has demonstrated that oxidative stress in response to various environmental insults or maternal dietary factors is responsible for aberrant DNA methylation [40,41,42,43]. In our study, hypermethylation of Pax3 and repressed gene expression, along with reduced TAC level, were observed in BaP-treated mouse embryos. Co-administration with vitamin E could rescue NTDs induced by BaP, partly normalized the TAC level, and attenuate the hypermethylation of Pax3 and the repressed gene expression. All these findings support the postulation that the oxidative stress may be causally involved in aberrant DNA methylation in the BaP-treated group, which might be the underlying mechanism for the development of NTDs caused by BaP. Consistent with this idea, a previous investigation revealed that epigallocatechin gallate, which is the major polyphenol in green tea and is believed to act as an antioxidant, could block hypermethylation of several neural tube closure essential genes induced by maternal diabetes, including Pax3 . Although it is not precisely known, previous investigations proposed that oxidative stress could influence the synthesis of S-adenosylmethionine , the expression of DNA methyltransferase , and the activity of the ten eleven translocation enzymes , which are all critical for DNA methylation. On the other hand, epigenetic regulation may also contribute to impairment of antioxidant gene expression . In this regard, the ROS production and epigenetic regulation establish an interconnected cycle, which would amplify external factors toward the progression of pathological disorders.
DNA methylation is considered tissue-specific and therefore appropriate biological samples for methylation study are crucial. One advantage of our study is that neural tissues from fetuses were used for assessment of methylation, which makes our methylation study more relevant to the outcome, as compared to those that used blood DNA as the surrogate . The limitation is that we did not investigate the expression of PAX3 gene in our human subjects, as fresh neural tissue samples of terminated NTD cases for RNA assay are extremely difficult to collect. Although our mouse experiment may provide a clue, additional studies are needed to correlate the methylation levels of the PAX3 gene body with PAX3 gene or protein expression values. Another limitation is that more case mothers reported folic acid supplementation than control mothers, which might cause the differences in methylation levels between the two groups. However, when analyses were performed by folic acid supplementation status in the case group, no differences in the methylation level of PAX3 were found between those with or those without folic acid use (Additional file 4: Figure S1), suggesting that folic acid has no impact on methylation level in the present study.
In conclusion, we found that hypermethylation of Pax3 and the downregulation of the gene may be important events in the development of NTDs during embryogenesis. A relationship exists between higher maternal serum PAH concentrations and PAX3 hypermethylation in fetal neural tissues, which is further supported by our mouse experiments. In addition, oxidative stress may be involved in the process of environmental exposure and methylation modification. Our study provides novel evidence on the interaction between genetic and environmental factors in the etiology of NTDs.
The human subjects were recruited from rural counties in Shanxi province of northern China (Pinding, Zezhou, Xiyang, Shouyang, and Taigu), where NTD prevalence is among the highest in the world . As described in more detail in our previous report , NTD cases were terminated fetuses affected by an NTD; controls were terminated non-malformed fetuses. Information on maternal sociodemographic characteristics, lifestyle, and folic acid supplementation was collected through in-person interviews and by viewing medical records. Maternal venous blood samples were collected at pregnancy termination. Fetal spinal cord and brain tissues were collected by autopsy performed by experienced pathologists. All samples were stored at − 80 °C until assay. The study protocol was approved by the Institutional Review Board of Peking University, and written informed consent was obtained from all participating women.
ICR mice of 8–9 weeks old were used, as described previously . Briefly, BaP (Sigma, USA) dissolved in corn oil was intraperitoneally administered into pregnant mice from E7 for four consecutive days (250–350 mg/kg). Mice in the vitamin E co-exposure group were fed with chow supplemented with the water-soluble (±)-α-tocopherol succinate form of vitamin E (Sigma, USA) beginning from E0.5 (0.125%, w/w) and treated with BaP from E7 (250 mg/kg). On E10.5, embryos were collected and carefully inspected for NTDs. All experimental procedures were approved by the Institutional Animal Care and Use Committee of Peking University (certificate no. LA2013-36).
DNA methylation level analysis
DNA from neural tissues of human subjects was extracted with QIAamp DNA Mini Kit (QIAGEN, Hilden, Germany) and E10.5 mouse DNA was extracted from frozen neural tissues using TIANamp Genomic DNA Kit (TIANGEN Biotech, Beijing, China). Then, 500 ng of genomic DNA from each sample was bisulfite-treated with EZ DNA methylation kit (Zymo Research, CA, USA). The bisulfite conversion reaction was performed in duplicate for each sample to minimize potential bias caused by variable conversion efficiency, and pooled bisulfite-treated DNA was used for subsequent analysis. Human methylation study was performed in two phases. In phase 1, Infinium HumanMethylation450 BeadChip assay (450K, Illumina, San Diego, CA, USA) was used for genomic methylation assay of bisulfite-treated DNA from ten NTD cases and eight non-malformed controls, which has been described in detail in our previous study . Methylation data of PAX3 gene was extracted from the array data and used for validation in the next phase.
In phase 2, 73 NTD cases and 29 non-malformed control fetuses were included. The locations of the amplicons to target the aberrant CpG regions are shown in Fig. 1a. Bisulphite DNA was amplified by PCR and primers for the PAX3 gene were designed using the online tool Epidesigner (www.epidesigner.com). Primer sequences are listed in Additional file 5: Table S4. After reverse transcription, fragmentation, and analysis on a mass spectrometer (Sequenom, Inc., San Diego, USA), EpiTYPER Analyzer software was used for translating mass signal patterns into quantitative DNA methylation levels of different CpG sites.
Methylation analysis of mouse neural tissues was assayed with the same methods performed in phase 2 of the human methylation study. The regions of gene sequence analyzed were the same as those in human PAX3 gene and three amplicons were designed, with two in the promoter region and one in the body region, as shown in Fig. 1b. See Additional file 6: Table S5 for PCR primer sequences.
PAH concentrations in maternal blood were determined with an Agilent 7890A-5975C gas chromatograph and mass spectrometer equipped with a HP-5MS capillary column (30 m × 0.25 mm × 0.25 μm), as described previously . In the present study, 51 mother-fetus pairs were available for PAH-methylation correlation analysis.
Oxidative stress assessment
The antioxidant capacity was determined by the reduced ferric reducing antioxidant power (FRAP) assay according to the manufacturer’s instructions (Beyotime Institute of Biotechnology, China) for determining the level of TAC of mouse embryos. The TAC aims to measure both small molecule and protein antioxidants, including polyphenols, flavonoids, vitamins, and enzymes like glutathione peroxidase and superoxide dismutase. Briefly, neural tissues from E10.5 embryo samples were homogenized in PBS. Freshly prepared FRAP reagent was warmed to 37 °C before use. Then, 5 μl of the diluted sample was added to 180 μl of the FRAP reagent. The absorbance of the mixture was measured at 593 nm using a Synergy 2 Multi-Mode Microplate Reader (BioTek, USA) after incubation for 4 min. The TAC of each sample was calculated from the standard curve constructed using FeSO4 solution, and the results were expressed as mmol/g protein, adjusted according to the protein concentration of the samples. All samples were loaded in duplicate. Protein concentration was determined using the BCA protein assay kit (Beyotime Institute of Biotechnology, China).
SOD activity was determined using the total superoxide dismutase assay kit with WST-8 (Beyotime Institute of Biotechnology, China). Briefly, 20 μl of sample (E10.5 embryo samples were homogenized in lysis buffer) was mixed with 160 μl of WST-8/enzyme working solution. Then, 20 μl of reaction triggering working solution was added. After incubation at 37 °C for 30 min, the absorbance was determined at 450 nm using a Synergy 2 Multi-Mode Microplate Reader (BioTek, USA). The SOD activity was expressed as U/mg protein.
RNA isolation and real-time PCR
RNA was isolated from E10.5 embryos neural tissues using Trizol (Invitrogen); genomic DNA was removed by DNase I digestion (DNA-free, Ambion) and then reverse-transcribed using random hexamers (Superscript VILO cDNA synthesis kit). The abundance of mRNA of Pax3 were analyzed using real-time PCR (iTaq™ Universal SYBR Green Supermix, BioRad) on a 7500 Fast Real Time PCR system (Applied Biosystems), with each sample analyzed in triplicate. Primers are listed in Additional file 7: Table S6. Relative quantification of gene expression level was normalized according to the Gapdh gene expression.
Whole mount in situ hybridization
Whole mount in situ hybridization on E9.5 control embryos and embryos treated with BaP was performed according to the procedures described by Yun et al. . Pax3 probe was cloned by real-time PCR into pGEM-T (Promega) and used to generate digoxygenin-labeled cRNA probes by reverse transcription using T7 RNA polymerase (Roche). For detection, anti-digoxigenin-AP antibody (1:2000, Roche) in 1% sheep serum was used and incubated overnight. Color detection was carried out using NBT/BCIP developing solution (Roche) in NTMT. After color development, embryos were imaged with a DFC490 camera (Leica), and then embedded. Sections of 40 μm thickness were obtained using a vibratome.
In the human subject study, differences in proportions of population characteristics between NTD cases and controls were examined with Pearson’s χ2 test. Independent t test was performed to evaluate the difference in methylation of CpGs between NTDs and controls in phase 1, and adjusted for multiple testing with the Benjamini-Hochberg false discovery rate (FDR) methods. In phase 2, Shapiro-Wilk test was used to examine the distribution of methylation values of NTD cases and controls, and independent samples t test was used to identify CpG sites that were differentially methylated between cases and controls. Odds ratio (OR) was calculated by logistic regression to evaluate the association between higher methylation levels of PAX3 with the risk of NTDs, adjusting for maternal education, occupation, periconceptional folic acid use, fetal sex, and gestational age. Correlation between differentially methylated CpG sites in neural tissues and PAH concentrations in maternal serum was analyzed with Spearman’s correlation analysis. In the mouse study, data on methylation level of Pax3 gene, oxidative stress markers, and the abundance of mRNA were presented as mean ± SE (SD). A one-way analysis of variance (ANOVA) followed by LSD (equal variances assumed) or Dunnett’s T3 (equal variances not assumed) was used for testing the differences between groups. A two-tail P value of < 0.05 was considered statistically significant. Statistical analyses were conducted using SPSS 18.0.
False discovery rate
Ferric reducing antioxidant power
Neural tube defects
Polycyclic aromatic hydrocarbon
Reactive oxygen species
Total antioxidant capacity
Mitchell LE. Epidemiology of neural tube defects. Am J Med Genet C Semin Med Genet. 2005;135c(1):88–94.
Greene ND, Stanier P, Copp AJ. Genetics of human neural tube defects. Hum Mol Genet. 2009;18(R2):R113–29.
Harris MJ, Juriloff DM. Mouse mutants with neural tube closure defects and their role in understanding human neural tube defects. Birth Defects Res A Clin Mol Teratol. 2007;79(3):187–210.
Wilde JJ, Petersen JR, Niswander L. Genetic, epigenetic, and environmental contributions to neural tube closure. Annu Rev Genet. 2014;48:583–611.
Greene ND, Stanier P, Moore GE. The emerging role of epigenetic mechanisms in the etiology of neural tube defects. Epigenetics. 2011;6(7):875–83.
Weinhold B. Epigenetics: the science of change. Environ Health Perspect. 2006;114(3):A160–7.
Morgan HD, Santos F, Green K, Dean W, Reik W. Epigenetic reprogramming in mammals. Hum Mol Genet. 2005;14 Spec No 1:R47–58.
Kim JK, Samaranayake M, Pradhan S. Epigenetic mechanisms in mammals. Cell Mol Life Sci. 2009;66(4):596–612.
Wang L, Wang F, Guan J, Le J, Wu L, Zou J, et al. Relation between hypomethylation of long interspersed nucleotide elements and risk of neural tube defects. Am J Clin Nutr. 2010;91(5):1359–67.
Wu L, Wang L, Shangguan S, Chang S, Wang Z, Lu X, et al. Altered methylation of IGF2 DMR0 is associated with neural tube defects. Mol Cell Biochem. 2013;380(1–2):33–42.
Bai B, Zhang Q, Liu X, Miao C, Shangguan S, Bao Y, et al. Different epigenetic alterations are associated with abnormal IGF2/Igf2 upregulation in neural tube defects. PLoS One. 2014;9(11):e113308.
Liu Z, Wang Z, Li Y, Ouyang S, Chang H, Zhang T, et al. Association of genomic instability, and the methylation status of imprinted genes and mismatch-repair genes, with neural tube defects. Eur J Hum Genet. 2012;20(5):516–20.
Shangguan S, Wang L, Chang S, Lu X, Wang Z, Wu L, et al. DNA methylation aberrations rather than polymorphisms of FZD3 gene increase the risk of spina bifida in a high-risk region for neural tube defects. Birth Defects Res A Clin Mol Teratol. 2015;103(1):37–44.
Tian T, Wang L, Shen Y, Zhang B, Finnell RH, Ren A. Hypomethylation of GRHL3 gene is associated with the occurrence of neural tube defects. Epigenomics. 2018;10(7):891–901.
Rochtus A, Izzi B, Vangeel E, Louwette S, Wittevrongel C, Lambrechts D, et al. DNA methylation analysis of Homeobox genes implicates HOXB7 hypomethylation as risk factor for neural tube defects. Epigenetics. 2015;10(1):92–101.
Phelan SA, Ito M, Loeken MR. Neural tube defects in embryos of diabetic mice: role of the Pax-3 gene and apoptosis. Diabetes. 1997;46(7):1189–97.
Greene ND, Massa V, Copp AJ. Understanding the causes and prevention of neural tube defects: insights from the splotch mouse model. Birth Defects Res A Clin Mol Teratol. 2009;85(4):322–30.
Pingault V, Ente D, Dastot-Le Moal F, Goossens M, Marlin S, Bondurand N. Review and update of mutations causing Waardenburg syndrome. Hum Mutat. 2010;31(4):391–406.
Hol FA, Hamel BC, Geurds MP, Mullaart RA, Barr FG, Macina RA, et al. A frameshift mutation in the gene for PAX3 in a girl with spina bifida and mild signs of Waardenburg syndrome. J Med Genet. 1995;32(1):52–6.
Lu W, Zhu H, Wen S, Laurent C, Shaw GM, Lammer EJ, et al. Screening for novel PAX3 polymorphisms and risks of spina bifida. Birth Defects Res A Clin Mol Teratol. 2007;79(1):45–9.
Ren A, Qiu X, Jin L, Ma J, Li Z, Zhang L, et al. Association of selected persistent organic pollutants in the placenta with the risk of neural tube defects. Proc Natl Acad Sci U S A. 2011;108(31):12770–5.
Lin S, Ren A, Wang L, Huang Y, Wang Y, Wang C, et al. Oxidative stress and apoptosis in benzo[a]pyrene-induced neural tube defects. Free Radic Biol Med. 2018;116:149–58.
Wang L, Lin S, Zhang J, Tian T, Jin L, Ren A. Fetal DNA hypermethylation in tight junction pathway is associated with neural tube defects: a genome-wide DNA methylation analysis. Epigenetics. 2017;12(2):157–65.
Wang B, Jin L, Ren A, Yuan Y, Liu J, Li Z, et al. Levels of polycyclic aromatic hydrocarbons in maternal serum and risk of neural tube defects in offspring. Environ Sci Technol. 2015;49(1):588–96.
Fleming A, Copp AJ. Embryonic folate metabolism and mouse neural tube defects. Science. 1998;280(5372):2107–9.
Fleming A, Copp AJ. A genetic risk factor for mouse neural tube defects: defining the embryonic basis. Hum Mol Genet 2000;9(4):575–581.
Nye JS, Balkin N, Lucas H, Knepper PA, McLone DG, Charrow J. Myelomeningocele and Waardenburg syndrome (type 3) in patients with interstitial deletions of 2q35 and the PAX3 gene: possible digenic inheritance of a neural tube defect. Am J Med Genet. 1998;75(4):401–8.
Agopian AJ, Bhalla AD, Boerwinkle E, Finnell RH, Grove ML, Hixson JE, et al. Exon sequencing of PAX3 and T (brachyury) in cases with spina bifida. Birth Defects Res A Clin Mol Teratol. 2013;97(9):597–601.
Wei D, Loeken MR. Increased DNA methyltransferase 3b (Dnmt3b)-mediated CpG island methylation stimulated by oxidative stress inhibits expression of a gene required for neural tube and neural crest development in diabetic pregnancy. Diabetes. 2014;63(10):3512–22.
Zhong J, Xu C, Reece EA, Yang P. The green tea polyphenol EGCG alleviates maternal diabetes-induced neural tube defects by inhibiting DNA hypermethylation. Am J Obstet Gynecol. 2016;215(3):368. e1-.e10
Rochtus A, Winand R, Laenen G, Vangeel E, Izzi B, Wittevrongel C, et al. Methylome analysis for spina bifida shows SOX18 hypomethylation as a risk factor with evidence for a complex (epi)genetic interplay to affect neural tube development. Clin Epigenetics. 2016;8:108.
Bell AC, Felsenfeld G. Methylation of a CTCF-dependent boundary controls imprinted expression of the Igf2 gene. Nature. 2000;405(6785):482–5.
Zhang X, Yazaki J, Sundaresan A, Cokus S, Chan SW, Chen H, et al. Genome-wide high-resolution mapping and functional analysis of DNA methylation in arabidopsis. Cell. 2006;126(6):1189–201.
Lev Maor G, Yearim A, Ast G. The alternative role of DNA methylation in splicing regulation. Trends Genet. 2015;31(5):274–80.
Yearim A, Gelfman S, Shayevitch R, Melcer S, Glaich O, Mallm JP, et al. HP1 is involved in regulating the global impact of DNA methylation on alternative splicing. Cell Rep. 2015;10(7):1122–34.
Tang WY, Levin L, Talaska G, Cheung YY, Herbstman J, Tang D, et al. Maternal exposure to polycyclic aromatic hydrocarbons and 5’-CpG methylation of interferon-gamma in cord white blood cells. Environ Health Perspect. 2012;120(8):1195–200.
Herbstman JB, Tang D, Zhu D, Qu L, Sjodin A, Li Z, et al. Prenatal exposure to polycyclic aromatic hydrocarbons, benzo[a]pyrene-DNA adducts, and genomic DNA methylation in cord blood. Environ Health Perspect. 2012;120(5):733–8.
Fang X, Dong W, Thornton C, Willett KL. Benzo[a]pyrene effects on glycine N-methyltransferase mRNA expression and enzyme activity in Fundulus heteroclitus embryos. Aquat Toxicol. 2010;98(2):130–8.
Sadikovic B, Rodenhiser DI. Benzopyrene exposure disrupts DNA methylation and growth dynamics in breast cancer cells. Toxicol Appl Pharmacol. 2006;216(3):458–68.
Dolinoy DC, Huang D, Jirtle RL. Maternal nutrient supplementation counteracts bisphenol A-induced DNA hypomethylation in early development. Proc Natl Acad Sci U S A. 2007;104(32):13056–61.
Bernal AJ, Dolinoy DC, Huang D, Skaar DA, Weinhouse C, Jirtle RL. Adaptive radiation-induced epigenetic alterations mitigated by antioxidants. FASEB J. 2013;27(2):665–71.
Franco R, Schoneveld O, Georgakilas AG, Panayiotidis MI. Oxidative stress, DNA methylation and carcinogenesis. Cancer Lett. 2008;266(1):6–11.
Angelini F, Pagano F, Bordin A, Milan M, Chimenti I. The impact of environmental factors in influencing epigenetics related to oxidative states in the cardiovascular system. Oxid Med Cell Longev. 2017;2017:2712751.
Hitchler MJ, Domann FE. An epigenetic perspective on the free radical theory of development. Free Radic Biol Med. 2007;43(7):1023–36.
Weitzman SA, Turk PW, Milkowski DH, Kozlowski K. Free radical adducts induce alterations in DNA cytosine methylation. Proc Natl Acad Sci U S A. 1994;91(4):1261–4.
Chia N, Wang L, Lu X, Senut MC, Brenner C, Ruden DM. Hypothesis: environmental regulation of 5-hydroxymethylcytosine by oxidative stress. Epigenetics. 2011;6(7):853–6.
Cyr AR, Hitchler MJ, Domann FE. Regulation of SOD2 in cancer by histone modifications and CpG methylation: closing the loop between redox biology and epigenetics. Antioxid Redox Signal. 2013;18(15):1946–55.
Liu J, Zhang L, Li Z, Jin L, Zhang Y, Ye R, et al. Prevalence and trend of neural tube defects in five counties in Shanxi province of northern China, 2000 to 2014. Birth Defects Res A Clin Mol Teratol. 2016;106(4):267–74.
Pai YJ, Leung KY, Savery D, Hutchin T, Prunty H, Heales S, et al. Glycine decarboxylase deficiency causes neural tube defects and features of non-ketotic hyperglycinemia in mice. Nat Commun. 2015;6:6388.
Technical assistance from Kit-Yi Leung, Yanhui Xiao, and Wenli Ni is gratefully acknowledged. All participants contributing to the original study subjects recruit are acknowledged.
This work was supported by the National Key Research and Development Program, Ministry of Science and Technology, the Peoples’ Republic of China (Grant No. 2016YFC1000501), and grants from the National Natural Science Foundation of China (Grant No. 31371523 and 81773441), the Ministry of Education of China (Grant No. 20130001110064), and the UK Medical Research Council (N003713, K022741). NG is supported by Great Ormond Street Hospital Children’s Charity and by the National Institute for Health Research Biomedical Research Centre at Great Ormond Street Hospital for Children NHS Foundation Trust and University College London.
Availability of data and materials
The data sets supporting the results of this article are included within the article and its additional files.
Ethics approval and consent to participate
The study protocol was approved by the Institutional Review Board of Peking University, and written informed consent was obtained from all participating women. The animal experiments were approved by the Institutional Animal Care and Use Committee of Peking University (certificate no. LA2013-36).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Characteristics of NTD cases and controls in phase 1 for methylation assay. (DOCX 18 kb)
Table S2. Methylation of PAX3 gene in phase 1 using the HumanMethylation450 BeadChip assay. (DOCX 24 kb)
Table S3. Correlation analysis of DNA methylation of PAX3 in fetal neural tissues and PAH concentrations in maternal serum (N = 51 mother-fetus pairs). (DOCX 22 kb)
Figure S1. PAX3 methylation pattern assayed by Sequenom EpiTYPER in NTD cases with and without folic acid supplementation. (DOCX 148 kb)
Table S4. The PCR primer sequences in Sequenom EpiTYPER sequencing. (DOCX 15 kb)
Table S5. The PCR primer sequences for mouse methylation analysis. (DOCX 15 kb)
Table S6. Sequences of primers for real-time PCR in mouse study. (DOCX 15 kb)