PLoS OneplosplosonePLoS ONE1932-6203Public Library of ScienceSan Francisco, USA212979373031146PONE-D-10-0160310.1371/journal.pone.0016506Research ArticleBiologyComputational BiologyGenomicsGenome Analysis ToolsGenome-Wide Association StudiesDevelopmental BiologyMorphogenesisBirth DefectsGeneticsHuman GeneticsGenome-Wide Association StudiesEpigeneticsGenome-Wide Association StudiesGenomicsGenome Analysis ToolsGenome-Wide Association StudiesGenomic MedicineMedicineCardiovascularCongenital Heart DiseaseObstetrics and GynecologyPregnancyPregnancy ComplicationsMaternal Genome-Wide DNA Methylation Patterns and Congenital Heart DefectsMaternal Methylation and Congenital Heart DefectsChowdhuryShimul*EricksonStephen W.MacLeodStewart L.ClevesMario A.HuPingKarimMohammad A.HobbsCharlotte A.Department of Pediatrics, College of Medicine, University of Arkansas for Medical Sciences, Arkansas Children's Hospital Research Institute, Little Rock, Arkansas, United States of AmericaTjwaMarcEditorUniversity of Frankfurt - University Hospital Frankfurt, Germany* E-mail: schowdhury@uams.edu

Conceived and designed the experiments: SC SWE SLM MAC CAH. Performed the experiments: SC. Analyzed the data: SWE SC PH MAC. Contributed reagents/materials/analysis tools: CAH, MAC. Wrote the paper: SC SWE SLM MAC CAH. Assisted in study conception and design and interpretation of results: SC SWE SLM MAC MAK CAH.

2011241201161e16506392010312011Chowdhury et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

The majority of congenital heart defects (CHDs) are thought to result from the interaction between multiple genetic, epigenetic, environmental, and lifestyle factors. Epigenetic mechanisms are attractive targets in the study of complex diseases because they may be altered by environmental factors and dietary interventions. We conducted a population based, case-control study of genome-wide maternal DNA methylation to determine if alterations in gene-specific methylation were associated with CHDs. Using the Illumina Infinium Human Methylation27 BeadChip, we assessed maternal gene-specific methylation in over 27,000 CpG sites from DNA isolated from peripheral blood lymphocytes. Our study sample included 180 mothers with non-syndromic CHD-affected pregnancies (cases) and 187 mothers with unaffected pregnancies (controls). Using a multi-factorial statistical model, we observed differential methylation between cases and controls at multiple CpG sites, although no CpG site reached the most stringent level of genome-wide statistical significance. The majority of differentially methylated CpG sites were hypermethylated in cases and located within CpG islands. Gene Set Enrichment Analysis (GSEA) revealed that the genes of interest were enriched in multiple biological processes involved in fetal development. Associations with canonical pathways previously shown to be involved in fetal organogenesis were also observed. We present preliminary evidence that alterations in maternal DNA methylation may be associated with CHDs. Our results suggest that further studies involving maternal epigenetic patterns and CHDs are warranted. Multiple candidate processes and pathways for future study have been identified.

Introduction

Birth defects are the leading cause of infant mortality, and congenital heart defects (CHDs) are among the most fatal of all birth defects [1]. Multiple genes have been implicated in CHD development [2], but for the majority of infants diagnosed with a CHD, an established causative gene or teratogenic agent cannot be identified [3], [4]. Identification of risk factors for CHDs are further complicated by the fact that both the maternal and fetal genetic susceptibilities may affect the intrauterine environment during gestation, and both may contribute to the development of CHDs [5], [6]. The complex nature of non-syndromic altered cardiogenesis presents a significant challenge to investigators interested in deciphering the etiology of CHDs.

Maternal folate supplementation has been shown to reduce the risk of CHDs [7]. Folate plays a key role in multiple cellular processes that are in increased demand during pregnancy, including amino acid and nucleotide synthesis, as well as DNA methylation [8]. Previous studies suggest that alterations in DNA methylation may contribute to the development of birth defects [9], [10], [11]. Multiple factors that have been shown to affect DNA methylation patterns including diet, genotype, and environmental exposures [10], [12], [13] have also been recognized as maternal factors implicated in the modification of fetal phenotypes [14], [15], [16]. Maternal genes and environmental exposures may modify the fetus through direct interactions or through alterations in the intrauterine environment during development [17], [18]. Although alterations in maternal one-carbon plasma metabolites indicative of a cellular hypomethylation status have been associated with an increased risk of CHDs [19], [20], a comprehensive study assessing maternal genomic methylation patterns has not yet been conducted.

To understand the mechanisms whereby gene-environment interactions lead to complex disease, epigenetic mechanisms must be considered [21]. DNA methylation is the best characterized epigenetic mechanism. It involves the covalent addition of a methyl group to the cytosine base within the context of CpG dinucleotides. CpG methylation is involved in gene silencing, genomic imprinting, chromosomal stability and protection against parasitic repetitive elements in various cells and tissues [22], [23]. Tightly controlled DNA methylation is essential in early fetal development, and in regulating genomic programming [24], [25]. DNA methylation patterns are tissue-specific and vary depending on cell type [26]. However, DNA methylation patterns of DNA isolated from peripheral blood cells have been shown to be a potential marker of exposure and disease [14], [27].

Through the use of genome-wide DNA methylation array technology, we sought to determine if alterations in maternal gene-specific DNA methylation were associated with CHDs. We then further assessed candidate genes that were differentially methylated between mothers who had CHD-affected pregnancies and control mothers, and identified potential biological processes and pathways of biological relevance to CHDs.

Results

Genome-wide gene-specific DNA methylation was measured in 367 participants (180 cases and 187 controls). Various lifestyle characteristics were analyzed in our sample population to determine if these differed between cases and controls (Table 1). The majority of cases (61.7%) and controls (60.4%) were less than 30 years old. The study population consisted mostly of Caucasian women. Smoking was the only selected covariate that was significantly more prevalent in cases (30.0%) than in controls (18.7%; P = 0.015).

10.1371/journal.pone.0016506.t001Selected characteristics of cases and controls.
Cases (n = 180) n (%)Controls (n = 187) n (%)p-value1
Age (y)
<30111 (61.7)113 (60.4)
≥3069 (38.3)74 (39.6)0.831
Race
Caucasian136 (75.6)148 (79.1)
African American32 (17.8)22 (11.8)
Others12 (6.7)17 (9.1)0.205
Smoker
Missing0 (0.0)1 (0.5)
No126 (70.0)151 (80.7)
Yes54 (30.0)35 (18.7)0.015
Alcohol drinker
Missing1 (0.6)2 (1.1)
No93 (51.7)85 (45.4)
Yes86 (47.8)100 (53.5)0.294
Vitamin supplementation
Missing0 (0.0)1 (0.5)
No99 (55.0)111 (59.4)
Yes81 (45.0)75 (40.1)0.400
Maternal education
Missing5 (2.8)0 (0.0)
High school or lower83 (46.1)79 (42.2)
College education or higher92 (51.1)108 (57.8)0.342
Household income
Missing13 (7.2)12 (6.4)
Less than $10,00034 (18.9)23 (12.3)
$10,000–$30,00060 (33.3)55 (29.4)
$30,000–$50,00039 (21.7)51 (27.3)
More than $50,00034 (18.9)46 (24.6)0.137
BMI (kg/m2)
Missing8 (4.4)9 (4.8)
Underweight (<18.5)3 (1.7)6 (3.2)
Normal (18.5–24.9)57 (31.7)78 (41.7)
Overweight (25.0–29.9)44 (24.4)34 (18.2)
Obese (≥30.0)68 (37.8)60 (32.1)0.118

Fisher's exact test for categorical variables.

To illustrate the representative distribution of methylation levels across more than 27,000 CpG sites, Figure 1A displays histograms of β-values for four randomly selected study subjects. In each sample, the distribution is bimodal with a high peak of hypomethylated CpG sites and a low peak of hypermethylated CpG sites. These individuals, therefore, have a large proportion of CpG sites characterized by low methylation, a smaller proportion of sites characterized by high methylation, and a certain proportion of sites in between the two peaks with moderate levels of methylation. The distributions displayed are typical of all subjects. The observed distribution of β-values was expected as the array design is biased to interrogate the promoter regions of genes within CpG islands.

10.1371/journal.pone.0016506.g001Distribution of β-values for selected subjects and CpG sites.

A) The distribution of β values for 27,249 CpG sites is displayed for four randomly selected subjects. The distribution for the four subjects is representative of the distribution for the entire sample population. B) The distribution of β values at selected CpG sites across all 367 study subjects. The four CpG sites were selected to illustrate typical methylation patterns observed in our sample population. Each CpG site displays skewed distributions, with a small amount of variation between individuals. Three of the four panels show substantially skewed distributions, most dramatically for the sites massed nearest 0 or 1. We did not observe drastic shifts in DNA methylation patterns in our sample population at individual CpG sites.

Figure 1B displays histograms of β-values at selected CpG sites across all study subjects. The four CpG sites were selected to illustrate typical methylation patterns observed in our sample population. Each CpG site has a typical methylation level, with a small amount of variation between individuals. Three of the four panels in Figure 1B show substantially skewed distributions, most dramatically for the data massed nearest 0 or 1. We did not observe drastic shifts in DNA methylation patterns in our sample population at individual CpG sites.

To determine if the associations between DNA methylation and cardiac defects were related to specific cardiac phenotypes, an analysis stratified by cardiac phenotypes was conducted (data not shown). Cardiac phenotypes included conotruncal, septal, and obstructive defects. Cases that had more than one cardiac defect type were included in all relevant strata. Among cases, 47.2% had atrial or ventricular septal defects (n = 85), 35.4% had right- and left-sided obstructive defects (n = 63), and 16.7% had conotruncal defects (n = 30). The 14 cases that did not fall within one of the three described main defect groups were excluded from the stratified analysis due to a small sample size. Analyses with sub-phenotype indicator variables included in the regression model were conducted and no statistically significant findings of sub-phenotype-specific methylation patterns were observed. Because gene-specific methylation patterns did not deviate significantly between cardiac phenotypes, cases were combined for subsequent analyses to maximize power to test our study hypotheses.

Figure 2 displays a quantile-quantile (Q-Q) plot of –log10(p-values) for the 27,249 tests of association between gene-specific methylation and disease status. The observed quantiles are consistently higher than their expected values under the null hypothesis of no disease association, providing evidence of site-specific disease association for a large number of CpG sites.

10.1371/journal.pone.0016506.g002Quantile-Quantile (Q-Q) plot of observed versus expected p-values.

A quantile-quantile (Q-Q) plot of –log10(p-values) for the tests of association between gene-specific methylation and disease status. The observed quantiles are consistently higher than their expected values under the null hypothesis of no disease association, providing evidence of site-specific disease association for a large number of CpG sites.

We identified 425 CpG sites that were differentially methylated between cases and controls with P<0.005, although an FDR threshold of 0.05 was not reached for any single CpG site (Table S1). The top 50 CpG sites ranked by statistical significance are listed in Table 2. The 425 CpG sites of interest encompassed 415 genes. The vast majority of these CpG sites were located within CpG islands (386 CpG sites, 90.8%). Given that the array platform contains 76% of sites within CpG islands, this result was statistically significant (P<0.001). Additionally, we observed that case mothers were hypermethylated in the majority of loci (379 sites, 89.2%) compared to hypomethylated loci (46 sites, 10.8%). Finally, the list of differentially methylated genes includes four genes previous shown to be imprinted (Table 3) as well as 14 miRNA sites (Table 4).

10.1371/journal.pone.0016506.t002Top 50 differentially methylated genes ranked by statistical significance.
SYMBOLGene nameGene_IDp-value1q-value2Chr3CpG IslandDirection4
CTSL2cathepsin L215151.6E-050.133789Yes+
RREB1ras responsive element binding protein 162391.6E-050.133786Yes+
BSNbassoon (presynaptic cytomatrix protein)89272.2E-050.133783Yes+
HSPC268hypothetical protein HSPC2681547914.6E-050.170867Yes+
NEUROG1neurogenin 147625.5E-050.170865Yes+
MAP4K5mitogen-activated protein kinase kinase kinase kinase 5111837.3E-050.1708614Yes+
FTMTferritin mitochondrial940337.4E-050.170865Yes-
CPEcarboxypeptidase E13638.8E-050.170864Yes+
EXT2exostosin 221320.00010.1708611Yes+
UBE3Aubiquitin protein ligase E3A73370.000130.1708615Yes+
DNAJB6DnaJ (Hsp40) homolog, subfamily B, member 6100490.000140.170867Yes+
UQCRHubiquinol-cytochrome c reductase hinge protein73880.000140.170861Yes+
ORF1-FL49chromosome 5 open reading frame 32844180.000160.170865Yes+
H2AFY2H2A histone family, member Y2555060.000180.1708610Yes+
PKN2protein kinase N255860.000220.170861Yes+
STRN3striatin, calmodulin binding protein 3299660.000220.1708614Yes+
C8orf35chromosome 8 open reading frame 35551740.000220.170868Yes+
NDUFA7NADH dehydrogenase (ubiquinone) 1 alpha subcomplex, 747010.000230.1708619Yes-
SLC36A4solute carrier family 36 (proton/amino acid symporter), member 41201030.000230.1708611Yes+
C6orf72chromosome 6 open reading frame 721162540.000230.170866Yes+
C6orf199chromosome 6 open reading frame 1992212640.000250.170866Yes+
CCDC74Acoiled-coil domain containing 74A905570.000260.170862Yes+
TERF1telomeric repeat binding factor (NIMA-interacting) 170130.000280.170868Yes+
SNX3sorting nexin 387240.000280.170866Yes+
HNRPH3heterogeneous nuclear ribonucleoprotein H3 (2H9)31890.00030.1708610Yes+
CGNL1cingulin-like 1849520.000310.1708615Yes-
ITM2Bintegral membrane protein 2B94450.000310.1708613Yes+
ZNF544zinc finger protein 544273000.000320.1708619Yes+
MGC14376hypothetical protein MGC14376849810.000350.1708617Yes+
ICAM3intercellular adhesion molecule 333850.000360.1708619No+
RRP22Ras-related protein 22106330.000360.1708622Yes+
AGTPBP1ATP/GTP binding protein 1232870.000370.170869Yes+
WNT5Awingless-type MMTV integration site family, member 5A74740.000370.170863Yes+
GMDSGDP-mannose 4,6-dehydratase27620.000390.170866Yes+
XPR1xenotropic and polytropic retrovirus receptor 192130.00040.170861Yes+
GDF3growth differentiation factor 395730.000410.1708612Yes+
RASEFRAS and EF-hand domain containing1581580.000410.170869Yes+
KCTD4potassium channel tetramerisation domain containing 43866180.000410.1708613No-
RFXAPregulatory factor X-associated protein59940.000410.1708613Yes+
KCNA4potassium voltage-gated channel, shaker-related subfamily, member 437390.000420.1708611Yes+
EGFRepidermal growth factor receptor19560.000430.170867Yes+
LGR6leucine-rich repeat-containing G protein-coupled receptor 6593520.000450.170861Yes+
MAPK13mitogen-activated protein kinase 1356030.000470.170866Yes+
SSFA2sperm specific antigen 267440.000490.170862Yes+
PHCAphytoceramidase, alkaline553310.000490.1708611Yes+
TBC1D1TBC1 (tre-2/USP6, BUB2, cdc16) domain family, member 1232160.000490.170864Yes+
LRRC3Bleucine rich repeat containing 3B1161350.000490.170863Yes+
SHCBP1SHC SH2-domain binding protein 1798010.00050.1708616Yes+
PMS2PMS2 postmeiotic segregation increased 2 (S. cerevisiae)53950.000530.170867Yes+
ZNF304zinc finger protein 304573430.000530.1708619Yes+

p-values were calculated by multiple linear regression and randomization testing (see Materials and Methods).

q-values were calculated using the qvalue package in the R statistical programming environment under default settings (see Materials and Methods).

Chr =  chromosome.

Direction of methylation in cases.

10.1371/journal.pone.0016506.t003Differentially methylated imprinted genes.
SYMBOLGene functionp-value1q-value2Direction of methylation in casesChromosomeCpG Island
UBE3AImprinted0.00012620.1708606+15Yes
ZNF264Imprinted0.00082940.1743975+19Yes
SNURFImprinted0.00395260.1822456-15Yes
INSImprinted0.00439080.1822456-11No

p-values were calculated by multiple linear regression and randomization testing (see Materials and Methods).

q-values were calculated using the qvalue package in the R statistical programming environment under default settings (see Materials and Methods).

10.1371/journal.pone.0016506.t004Differentially methylated miRNAs.
GenemiRNAp-value1q-value2Direction of methylation in casesChromosomeCpG Island
HOXB4HSA-MIR-10A;0.00433840.1822456+17Yes
BTG4HSA-MIR-34B;HSA-MIR-34C;0.00066630.1708606+11Yes
CCDC55HSA-MIR-423;0.00295420.1822456-17Yes
SACM1LHSA-MIR-565;0.00080430.1743975+3Yes
HSPC268HSA-MIR-594;0.00004630.1708606+7Yes
HSPC268HSA-MIR-594;0.00457670.1822456+7Yes
WDR37HS_104;0.00333060.1822456+10Yes
DUSP4HS_128;0.00340430.1822456+8Yes
C10orf4HS_138;0.00261440.1822456+10Yes
PIGOHS_140;0.00445430.1822456+9Yes
PPME1HS_176;0.00155820.1822456+11Yes
DUS2LHS_253;HS_253;0.00413220.1822456-16Yes
CDC7HS_255;0.00493220.1822456+1Yes
FBXO28HS_78;0.00150550.1822456+1Yes

p-values were calculated by multiple linear regression model and randomization testing (see Materials and Methods).

q-values were calculated using the qvalue package in the R statistical programming environment under default settings (see Materials and Methods).

Further analysis of the differentially methylated genes was conducted for potential biological significance. Gene Set Enrichment Analysis (GSEA) was performed to determine if the differentially methylated genes were significantly enriched in relevant biological processes and pathways of biological relevance to fetal heart development.

The 415 genes of interest were input into the GSEA program, and the top ten gene set overlaps were generated. Figure 3A summarizes the results of the most significant overlaps for the biological process analysis. Within the top ten most significant gene set overlaps, multiple functional categories associated with development were strongly enriched within our differentially methylated gene list. Specifically, the biological process categories of nucleic acid metabolism, signal transduction, anatomical structure development, multicellular organismal development, and system development have potential direct functional implications in the development of CHDs. The inclusion of these specific processes within the top ten most significant gene set overlaps provided evidence of strong biological relevance of the differentially methylated genes to fetal development. GSEA was also conducted to determine if significant canonical pathway overlaps could be established from the differentially methylated gene list (Figure 3B). Pathway analysis revealed functional overlaps that have previously been shown to be involved in embryonic heart development. These pathways included cytokine-cytokine receptor interaction, the regulation of actin cytoskeleton, tight junction, calcium signaling pathway, smooth muscle contraction and the Wnt signaling pathway. The GSEA results provided further evidence that the genes that were found to be differentially methylated appear to have biological relevance to CHDs. The results summaries generated from the GSEA, which includes the specific differentially methylated genes involved within each biological process and canonical pathway overlap, is included in Table S2 and Table S3 respectively.

10.1371/journal.pone.0016506.g003Gene Set Enrichment Analysis (GSEA) of biological processes and canonical pathways for differentially methylated genes.

A) The top ten significantly enriched biological processes are displayed with their corresponding –log10(p-value). Multiple processes appear to be directly related to fetal heart development. B) The top ten significantly enriched canonical pathways are displayed with the corresponding −log10(p-value). Multiple pathways had previously been implicated in CHDs.

Discussion

To our knowledge, we have performed the first population based maternal case-control study to provide evidence of association between maternal gene-specific DNA methylation and CHDs. Further analysis of the differentially methylated genes revealed functionally relevant enrichment in biological processes and pathways previously shown to be involved in fetal development.

No individual CpG site achieved genome-wide statistical significance. A combination of low variation in site-specific methylation between individuals and relatively modest sample size may account for the inability to reach statistical genome-wide significance. Thus, caution in interpreting the results of differentially methylated genes of interest must be noted. The advent of genome-wide association studies (GWAS) has led to an influx of studies structured to detect genetic variants associated with various human diseases, but limited studies have been conducted using genome-wide DNA methylation platforms. Issues regarding adjustments for multiple testing, statistically significant thresholds, and biologically relevant alteration levels, have not been established for high-throughput DNA methylation studies. The increased demand for genome wide epigenetic studies will lead to clearer definition of these issues in the future.

In the current study, gene-specific methylation was measured in DNA isolated from peripheral blood lymphocytes. We observed slight alterations in methylation between cases and controls. Factors such as genetic variation, metabolites involved in one-carbon metabolism, and environmental exposures are possible mechanisms of differential methylation between cases and controls. The interpretation of results from DNA methylation studies must always be evaluated with caution as epigenetic marks may vary depending on cell type and can be influenced by various factors. Previous studies have used Illumina Infinium technology to assess DNA methylation in peripheral blood lymphocytes and have identified CpG sites associated with disease [27], [28]. The accumulation of studies that have assessed DNA methylation patterns in lymphocytes in multiple diseases and exposures provides evidence that methylation patterns in peripheral blood cells may be useful in identifying biomarkers of various diseases.

It has been established that various maternal factors may cause altered fetal growth and development [29], but the impact of maternal genetics and epigenetics in the development of birth defects is not well understood. Maternal genes may exert effects through direct interactions with the fetus [30], [31] or through changes to the intrauterine environment. Maternal blood provides nutrients to the developing fetus and could provide an avenue of signal exchange. Multiple studies have been conducted to determine the impact of maternal genetic effects and the risk of birth defects [32], [33], [34]. Although these studies encompass a relatively small number of genes with respect to the genome, they suggest that maternal genetic effects, independent of fetal genotype, may influence the risk of birth defects. Additionally, alterations in gene expression in maternal blood are associated with adverse pregnancy outcomes [35], [36] providing further evidence that maternal genetic influence may impact the fetus. Epigenetics are considered to be the link between environmental exposures and their effects on the genome [37] and may explain how maternal genes may exert effects on the developing fetus.

Genes previously implicated in CHDs were identified within our differentially methylated genes. For example, multiple genes involved in the mitogen-activated protein kinases (MAPK) pathway are included in the list of differentially methylated genes. The MAPK pathway is involved in multiple cellular processes including transcription, proliferation, migration, survival, and differentiation [38]. It has been postulated that this family of genes may play a role in CHDs [39]. Additional genes implicated in heart development that were differentially methylated in our study include, but are not exclusive to, EGFR [40], GATA4 [41], and Wnt5a [42]. The inclusion of multiple miRNA targets as candidate genes is intriguing as the diverse role of miRNAs in cardiac function has received increased interest [43]. Although genetic studies have been conducted on these genes of interest, our results indicate that epigenetics should be considered as well.

Alterations in maternal DNA methylation may also indirectly affect the fetus through changes in the maternal environment. Changes in the maternal environment can result in direct changes in gene expression in the developing fetus [44], and if certain harmful exposures occur during organogenesis, various structural anomalies may occur [13]. Alterations in maternal methylation in genes involved in the maintenance of the intrauterine environment may lead to an increased susceptibility to teratogenic agents. Genes that were found to be differentially methylated, such as GPX3, provide an example of a gene involved in the sustainment of a normal maternal environment. The GPX3 gene is involved in antioxidant potential in the body. Altered expression in GPX3 may result in a lower anti-oxidant potential, which may potentially lead to increased oxidative stress for the fetus.

It is though that hypomethylation in the parental germline can influence genomic instability, and that this instability may subsequently cause genomic instability in the progeny [45]. Although it is unclear how maternal lymphocyte methylation patterns correlate with germ cells or the fertilized oocyte, maternal methylation patterns in peripheral blood DNA have recently been shown to correlate with infant methylation patterns in certain loci [46], [47]. The phenomenon of passing intact methylation patterns to future generations, termed transgenerational epigenetic inheritance, has been used to explain non-Mendelian inheritance patterns in complex diseases. Although speculation regarding whether epigenetic marks are transmitted intact from parent to offspring exists [48], the vast majority of the evidence is found in animal models [49].

The importance of epigenetic factors in the development of complex disease is now unquestioned. Embryonic heart development requires a series of highly complex, coordinated and rapid morphogenesis processes. Specifically for CHDs, epigenetic factors may contribute to cardiomyocyte differentiation and chamber morphogenesis [50], and it has been proposed that epigenetic mechanisms are necessary to coordinate genetic programs in heart development [51]. The complex series of events involved in heart development indicates that multiple genes and pathways are responsible for CHDs [52]. Although many studies have been conducted to investigate the developmental genetics of CHDs, limited studies are available regarding DNA methylation within these genes and pathways. Our analysis of biological processes revealed overlaps involving system, anatomical structure, and organismal development. The listed processes are involved in fetal development and provide candidate processes for additional epigenetic studies. Nucleic acid metabolism and signal transduction are also of importance to development as the increase demand for DNA synthesis and properly functioning signaling cascades are important during pregnancy [8], [53]. Canonical pathways enriched within our genes of interest have previously been implicated in heart development and some of the pathways appear to be interconnected. For example, alterations in intracellular calcium levels can affect smooth muscle contraction downstream [54]. Other significant pathways that have previously been implicated in heart development include ribosomal proteins [55], the actin cytoskeleton [52], Wnt signaling [56], and cytokine-cytokine receptor interaction [57]. Future epigenetic studies should be conducted within the candidate pathways and processes described.

Certain methodological limitations of this study should be considered. The blood obtained to isolate genomic DNA was collected after pregnancy. We have previously described the rationale for the use of case-control study designs in human birth defects research studies [20]. We did not measure DNA methylation changes during gestation and neonatal development, which represents a period of high epigenomic susceptibility [58]. Contrasting reports exist regarding gene-specific methylation changes over time [59], [60], and it is unclear how maternal methylation patterns may change during pregnancy. Another limitation of our study is the lack of infant methylation data to determine if differential methylation is observed between infants with CHDs and controls, and if these methylation patterns are associated with maternal DNA methylation patterns. Previous studies have shown global DNA hypomethylation in the proband for neural tube defects [61], [62], but studies in gene-specific and global methylation in a family-based design should be conducted for other congenital anomalies including CHDs. The design of the HumanMethylation27 BeadChip allows the assessment of over 27,000 CpG sites, which represents 0.01–0.18% of the single-copy genome [48]. The assay is focused primarily on one or two CpG sites within gene promoter regions. Thus, most of the neighboring CpG sites are not studied. Although our study has revealed a strong bias in DNA alterations occurring in CpG islands, the possibly exists for DNA methylation alterations occurring outside CpG islands and in CpG sites not included on the assay. Although the Illumina Infinium HumanMethylation27 platform provides a high-throughput and a reliable examination of a large number of genes [63], studies in birth defects investigating a larger number of CpG sites should be conducted. Our study provides additional evidence that maternal epigenetics and genetics should be included in the investigation of birth defects.

Alterations in maternal DNA methylation converge with evidence from previous studies that indicate folate-dependent genetic and metabolic susceptibilities increase the risk of CHDs. The study has notable strengths including analyzing genome-wide gene specific DNA methylation in a large population of mothers. Extensive statistical modeling was used to adjust for any potential experimental confounders. We also adjusted for various lifestyle factors that have been shown to affect DNA methylation. The data presented here strongly suggest that maternal epigenetics may be a determinant of CHD risk. The non-invasive nature of the determination of DNA methylation in lymphocytes isolated from peripheral blood draw makes it an attractive candidate for disease risk assessment. Identification of genetic and epigenetic profiles associated with an increased risk of CHDs may add new dimensions to preconception risk assessment as well as further elucidate mechanisms involved in CHD development. Because DNA methylation may be influenced by diet or drugs, epigenetics may provide an avenue for future therapeutic intervention and prevention of CHDs.

Materials and MethodsStudy Population

The National Birth Defects Research and Prevention Study (NBDPS) is an ongoing multi-site population-based case-control study investigating the etiology of 30 non-syndromic birth defects [64]. The NBDPS is the largest case-control study of birth defects ever conducted in the US. The study population and eligibility criteria for the NBDPS have been previously outlined [20]. For the current study, case women were Arkansas residents who participated in the NBDPS and who delivered a singleton live birth with a non-syndromic CHD. Cases for whom the pregnancy was also affected by a known single gene disorder, chromosomal abnormality, or syndrome were excluded. All diagnostic tests on cardiac NBDPS case infants were reviewed by a pediatric cardiologist to ensure uniform criteria were used for diagnoses. To be eligible for the study described, case infants were required to have at least one cardiac lesion, which included conotruncal, septal, and/or obstructive lesions. Using a classification system developed for the NBDPS, which incorporated three dimensions of cardiac phenotype, cardiac complexity, and extracardiac anomalies [65], case women for the current study included those who carried fetuses with simple or association CHDs. Control mothers were Arkansas residents who had a singleton live birth without birth defects during the same period as the cases, who participated in the NBDPS, and were randomly selected from birth certificate data or hospital discharge logs. Case and control mothers spoke either English or Spanish.

From eligible Arkansas NBDPS participants, we assessed DNA methylation in 180 cases and 187 controls for whom frozen blood samples were available. Blood samples were collected by research nurses during home visits in Arkansas that occurred between 2001 and 2008. Samples were collected by venipuncture after pregnancy. All participants signed informed consent approved by the UAMS Institutional Review Board.

Covariates

Information regarding selected lifestyle factors was obtained from in-home interviews conducted by research nurses using a Block Abbreviated Food-Frequency Questionnaire [66] and from a NBDPS-structured computer-assisted telephone interview. The information on covariates reflects data up to one month prior to the visit. Variables considered as covariates in this study due to their potential effects on DNA methylation status were age, race, body mass index (BMI), and the use of multivitamins, cigarettes, and alcohol.

Sample preparation and processing

Fasting blood samples were collected in EDTA-Vacutainer tubes and immediately chilled on ice before they were centrifuged at 4000×g for 10 min at 4°C. DNA was isolated from the stored frozen blood samples according the manufacturer's protocol for the PureGene DNA isolation kit (Gentra Systems, Minneapolis, MN).

DNA quantification and bisulfite conversion

Genomic DNA was quantified via the RNase P 20x assay (Applied Biosystems, Forest City, CA). In a total reaction volume of 25 µl, 2 µl of genomic DNA was used for absolute quantification for the RNase P assay on the ABI 7900HT Real Time PCR System, according to the manufacturer's protocol. After genomic DNA quantification, 500 ng of genomic DNA underwent bisulfite modification utilizing the EZ DNA Methylation-Direct Kit (Zymo Research, Orange, CA). The bisulfite-converted DNA was resuspended in 12 µl TE buffer and stored at -80°C until the samples were ready for analysis.

Infinium Assay for Methylation

DNA methylation analysis was conducted using the Illumina Infinium Human Methylation27 BeadChip. Each HumanMethylation27 BeadChip consists of 12 arrays, and up to 4 bead chips were processed simultaneously. The assay allows the interrogation of over 27,000 CpG sites located within the proximal promoter regions of over 14,000 consensus coding sequences (CCDS) genes throughout the genome [63]. In addition, the assay includes 110 miRNA promoters and also includes imprinted genes.

The standard protocol provided by Illumina was used for DNA methylation analysis. Briefly, 4 µl of bisulfite converted DNA was isothermally amplified at 37°C overnight. The amplified DNA product was fragmented by an endpoint enzymatic process. Fragmented DNA was precipitated, resuspended, and applied to an Infinium Human Methylation27 BeadChip and hybridized overnight. During hybridization, the amplified and fragmented DNA samples anneal to specific oligomers that are covalently linked to over 27,000 different bead types. Each bead type corresponds to the nucleotide identity and thus the methylation status at a bisulfite-converted cytosine in a specific CpG site. The bead chips were then subjected to a single-base extension reaction using the hybridized DNA as a template incorporating fluorescently labeled nucleotides of two different colors, each corresponding to the cytosine (methylated) or uracil (unmethylated) identity of the bisulfite-converted nucleotide at a specific CpG site. The fluorescently stained chip was imaged by the Illumina BeadArray Reader. Illumina's Genome Studio program was used to analyze BeadArray data to assign site-specific DNA methylation β-values to each CpG site. The proportion of methylation (β) for each subject at each CpG site was computed by first subtracting the background signal intensity of negative controls from both the methylated and unmethylated signals and then taking the ratio of the methylated signal intensity to the sum of both methylated and unmethylated signals. Thus, the β-value is a continuous variable ranging between 0 and 1.

Experimental Design and Quality Control

A total of one batch of two bead chips (24 samples) and eight batches of four bead chips were assayed. Cases and controls were randomly placed within the bead chips. Replicate samples were included in each batch for quality control and normalization measures. All replicate samples achieved a correlation coefficient of greater than 0.985 (data not shown). Using Illumina's Genome Studio, background adjusted β-values and assay control probe information was used for initial quality control. Subjects were required to pass three initial quality control criteria in order to be included in further analysis: 1) Subjects achieved a 95% CpG site call rate, which equals 26,250 sites called at an α<0.05 significance level; 2) Background signal for each subject was under 1000 units, determined from the controls dashboard in Genome Studio; and 3) Clear separation was observed in the non-polymorphic controls, also determined from the controls dashboard. Samples that did not meet these criteria were removed from the analysis. Individual CpG sites were examined using the detection p-value metric provided by Illumina in which the signal generated from each CpG site is compared to negative controls. A threshold of p<0.05 was used as a cutoff. CpG sites that did not reach this threshold were eliminated from the analysis. Using this criterion, 329 of 27,578 CpG sites (1.2%) had missing or invalid β-values in ten or more samples; these sites were dropped from subsequent analysis, leaving 27,249 CpG sites among 367 subjects.

Principal Component Analysis

Principal component analysis (PCA) was performed on the β-value data matrix of 367 subjects (rows) and 27,249 CpG sites (columns). Data were logit-transformed prior to PCA and subsequent association testing. As described above, 329 CpG sites with missing or invalid β-values in ten or more samples were removed from the data; before performing PCA, any remaining missing β-values were imputed using k nearest neighbor averaging [67].

The first 20 principal components (PCs) explain 62.8% of overall variance in the logit-transformed β-value matrix, and no single PC among the remaining 347 PCs explains more than 0.34% (Figure S1). Regression analysis (not shown) revealed that these 20 PCs are variously associated with experimental factors that had been previously shown to be potential confounders under the Infinium Methylation platform [27], as well as with several phenotypic and lifestyle factors, including case/control status, age, race, and body mass index (BMI), and the use of multivitamins, cigarettes, and alcohol. Regression analysis of the top 20 PCs therefore informed the selection of covariates to be included in the CpG-site-specific regression model, which is described in the following section.

Determination of differentially methylated sites between cases and controls

After initial quality control and PCA of the β value data, a sample of 367 subjects (180 cases and 187 controls) was tested for association between disease status (i.e. case/control) and maternal gene-specific methylation. In this analysis, each CpG site was tested for association by regressing its logit-transformed β value on case/control status using multiple linear regression models. Additional covariates included these models were experimental batch (as a factor), bisulfite-conversion (BSC) efficiency, age, BMI, race (as a factor), and indicator variables of perinatal vitamin usage, alcohol drinking, and tobacco smoking. BSC efficiency for each subject has two channels, red and green. The red and green channels are highly correlated but are not sufficiently to be collinear, thus both channels were included in the models (Figure S2).

Under standard theory a t-statistic, derived from the least-squares estimate of the case/control parameter and its standard error, was compared with the theoretical t-distribution to test the null hypothesis of no association between disease status and site-specific methylation, , where  = 1,…,27,249 indicates the CpG site. Because of the potential for spurious association due to chip bias, however, a randomization testing approach was used to evaluate statistical significance. For each pseudo-dataset, disease status was randomized in a two-stage process for each chip. First, disease status was randomly permuted among the subjects on that chip, and then with 50% probability, disease status was swapped between case and control for all subjects on the chip.

In each pseudo-dataset, therefore, each subject was equally likely to be assigned case or control status, which means that, is known to be (synthetically) true. However, because the randomization is performed separately within chips, the structure of the experimental design, and therefore the chance for spurious disease association due to chip effect, is retained. The empirical null distributions resulting from this randomization approach had heavier tails than the theoretical null distribution, leading to less significant p-values. The randomization approach is therefore a conservative approach and properly controls type-I error in the presence of chip-to-chip variation of β-values.

To properly account for the large number of statistical tests being performed, false discovery rate (FDR) q-values [68] were computed for each CpG site. In our context, FDR is the proportion of detected disease associations (under a given hypothesis testing procedure) that are false, while the q-value is defined as the minimum FDR at which a given test can be considered significant. The q-value is therefore the FDR equivalent of the p-value and was computed using the qvalue package in the R statistical programming environment, under default settings [69].

Biological Processes and Pathway Analysis

Based on the results from our association testing, a list of differentially methylated CpG sites with p-values less than 0.005 was generated. The CpG sites of interest were to mapped to their corresponding genes and the list of genes were then tested for potential overlaps with biological processes and pathways using the Gene Set Enrichment Analysis (GSEA) [70]. The software is designed to evaluate microarrays at the level of gene sets. Gene sets are defined as groups of genes that share common biological function, chromosomal location or regulation. The gene set database used in the analysis is established based on prior biological knowledge.

Supporting Information

Variance explained by first 20 (of 367) principal components. The proportion of variance explained by each of the first 20 of 367 principal components, from a PCA of the 367-by-27,249 matrix of logit-transformed methylation β-values. Combined, the first 20 PCs explain 62.8% of overall variance, with no single PC among the remaining 347 explaining more than 0.34%.

(PDF)

Click here for additional data file.

Bisulfite conversion efficiency metrics scatterplot. Scatterplot of the red and green channels of bisulfite conversion efficiency for the 367 samples used in association testing, overlaid with a locally weighted scatterplot smoothing (LOWESS) curve.

(PDF)

Click here for additional data file.

List of differentially methylated CpG sites at P<0.005.

(DOC)

Click here for additional data file.

Gene Set Enrichment Analysis for Biological Processes.

(XML)

Click here for additional data file.

Gene Set Enrichment Analysis for Canonical Pathways.

(XML)

Click here for additional data file.

We acknowledge the generous participation of the numerous Arkansas families that made this research study possible.

Competing Interests: The authors have declared that no competing interests exist.

Funding: This research is supported by grants from the National Institute for Child Health and Development (5-R01-HD039054-08), the Centers for Disease Control and Prevention (3-U50DD613236-10W1), the Arkansas Biosciences Institute (the major research component of the Tobacco Settlement Proceeds Act of 2000), and the University of Arkansas for Medical Sciences College of Medicine Children's University Medical Group Grant Program. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References HoyertDL 2003 Mortality associated with birth defects: influence of successive disease classification revisions. Birth Defects Res A Clin Mol Teratol 67 651 655 14703789 WolfMBassonCT 2010 The molecular genetics of congenital heart disease: a review of recent developments. Curr Opin Cardiol JenkinsKJCorreaAFeinsteinJABottoLBrittAE 2007 Noninherited risk factors and congenital cardiovascular defects: current knowledge: a scientific statement from the American Heart Association Council on Cardiovascular Disease in the Young: endorsed by the American Academy of Pediatrics. Circulation 115 2995 3014 17519397 PierpontMEBassonCTBensonDWJrGelbBDGigliaTM 2007 Genetic basis for congenital heart defects: current knowledge: a scientific statement from the American Heart Association Congenital Cardiac Defects Committee, Council on Cardiovascular Disease in the Young: endorsed by the American Academy of Pediatrics. Circulation 115 3015 3038 17519398 HobbsCAClevesMAKarimMAZhaoWMacLeodSL 2010 Maternal folate-related gene environment interactions and congenital heart defects. Obstet Gynecol 116 316 322 20664391 MitchellLEWeinbergCR 2005 Evaluation of offspring and maternal genetic effects on disease risk using a family-based approach: the "pent" design. Am J Epidemiol 162 676 685 16093287 ShawGMO'MalleyCDWassermanCRTolarovaMMLammerEJ 1995 Maternal periconceptional use of multivitamins and reduced risk for conotruncal heart defects and limb deficiencies among offspring. Am J Med Genet 59 536 545 8585581 OommenAMGriffinJBSarathGZempleniJ 2005 Roles for nutrients in epigenetic events. J Nutr Biochem 16 74 77 15681164 BlomHJShawGMdenHMFinnellRH 2006 Neural tube defects and folate: case far from closed. Nat Rev Neurosci 7 724 731 16924261 LiDPickellLLiuYWuQCohnJS 2005 Maternal methylenetetrahydrofolate reductase deficiency and low dietary folate lead to adverse reproductive outcomes and congenital heart defects in mice. Am J Clin Nutr 82 188 195 16002818 OkanoMBellDWHaberDALiE 1999 DNA methyltransferases Dnmt3a and Dnmt3b are essential for de novo methylation and mammalian development. Cell 99 247 257 10555141 WaterlandRAJirtleRL 2003 Transposable elements: targets for early nutritional effects on epigenetic gene regulation. Mol Cell Biol 23 5293 5300 12861015 OzanneSEConstanciaM 2007 Mechanisms of disease: the developmental origins of disease and the role of the epigenotype. Nat Clin Pract Endocrinol Metab 3 539 546 17581623 BaccarelliAWrightROBollatiVTarantiniLLitonjuaAA 2009 Rapid DNA Methylation Changes after Exposure to Traffic Particles. Am J Respir Crit Care Med CandiloroILDobrovicA 2009 Detection of MGMT promoter methylation in normal individuals is strongly associated with the T allele of the rs16906252 MGMT promoter single nucleotide polymorphism. Cancer Prev Res (Phila Pa) 2 862 867 CooneyCADaveAAWolffGL 2002 Maternal methyl supplements in mice affect epigenetic variation and DNA methylation of offspring. J Nutr 132 2393S 2400S 12163699 FurnessDLFenechMFKhongYTRomeroRDekkerGA 2008 One-carbon metabolism enzyme polymorphisms and uteroplacental insufficiency. Am J Obstet Gynecol 199 276 278 18771981 MacLennanNKJamesSJMelnykSPirooziAJerniganS 2004 Uteroplacental insufficiency alters DNA methylation, one-carbon metabolism, and histone acetylation in IUGR rats. Physiol Genomics 18 43 50 15084713 van DrielLMdeJRHelbingWAvan ZelstBDOttenkampJ 2008 Maternal global methylation status and risk of congenital heart diseases. Obstet Gynecol 112 277 283 18669723 HobbsCAClevesMAMelnykSZhaoWJamesSJ 2005 Congenital heart defects and abnormal maternal biomarkers of methionine and homocysteine metabolism. Am J Clin Nutr 81 147 153 15640474 LiuLLiYTollefsbolTO 2008 Gene-environment interactions and epigenetic basis of human diseases. Curr Issues Mol Biol 10 25 36 18525104 JonesPABaylinSB 2002 The fundamental role of epigenetic events in cancer. Nat Rev Genet 3 415 428 12042769 GronbaekKHotherCJonesPA 2007 Epigenetic changes in cancer. APMIS 115 1039 1059 18042143 EhrlichM 2003 Expression of various genes is controlled by DNA methylation during mammalian development. J Cell Biochem 88 899 910 12616529 ShamesDSMinnaJDGazdarAF 2007 DNA methylation in health, disease, and cancer. Curr Mol Med 7 85 102 17311535 ChalitchagornKShuangshotiSHourpaiNKongruttanachokNTangkijvanichP 2004 Distinctive pattern of LINE-1 methylation level in normal tissues and the association with carcinogenesis. Oncogene 23 8841 8846 15480421 TeschendorffAEMenonUGentry-MaharajARamusSJGaytherSA 2009 An epigenetic signature in peripheral blood predicts active ovarian cancer. PLoS One 4 e8274 20019873 BellCGTeschendorffAERakyanVKMaxwellAPBeckS 2010 Genome-wide DNA methylation analysis for diabetic nephropathy in type 1 diabetes mellitus. BMC Med Genomics 3 33 20687937 WuGBazerFWCuddTAMeiningerCJSpencerTE 2004 Maternal nutrition and fetal development. J Nutr 134 2169 2172 15333699 PoplikerMShatzAAviviAUllrichASchlessingerJ 1987 Onset of endogenous synthesis of epidermal growth factor in neonatal mice. Dev Biol 119 38 44 3539666 LetterioJJGeiserAGKulkarniABRocheNSSpornMB 1994 Maternal rescue of transforming growth factor-beta 1 null mice. Science 264 1936 1938 8009224 LiDRozenR 2006 Maternal folate deficiency affects proliferation, but not apoptosis, in embryonic mouse heart. J Nutr 136 1774 1778 16772436 JugessurAShiMGjessingHKLieRTWilcoxAJ 2010 Maternal genes and facial clefts in offspring: a comprehensive search for genetic associations in two population-based cleft studies from Scandinavia. PLoS ONE 5 e11493 20634891 DoolinMTBarbauxSMcDonnellMHoessKWhiteheadAS 2002 Maternal genetic effects, exerted by genes involved in homocysteine remethylation, influence the risk of spina bifida. Am J Hum Genet 71 1222 1226 12375236 SunCJZhangLZhangWY 2009 Gene expression profiling of maternal blood in early onset severe preeclampsia: identification of novel biomarkers. J Perinat Med 37 609 616 19681734 EnquobahrieDAWilliamsMAQiuCMuhieSYSlentz-KeslerK 2009 Early pregnancy peripheral blood gene expression and risk of preterm delivery: a nested case control study. BMC Pregnancy Childbirth 9 56 20003277 JaenischRBirdA 2003 Epigenetic regulation of gene expression: how the genome integrates intrinsic and environmental signals. Nat Genet 33 Suppl 245 254 12610534 RoyJKaziMHedinUThybergJ 2001 Phenotypic modulation of arterial smooth muscle cells is associated with prolonged activation of ERK1/2. Differentiation 67 50 58 11270123 BakerJE 2004 Oxidative stress and adaptation of the infant heart to hypoxia and ischemia. Antioxid Redox Signal 6 423 429 15025944 JoziasseICvan de SmagtJJSmithKBakkersJSieswerdaGJ 2008 Genes in congenital heart disease: atrioventricular valve formation. Basic Res Cardiol 103 216 227 18392768 GargVKathiriyaISBarnesRSchlutermanMKKingIN 2003 GATA4 mutations cause human congenital heart defects and reveal an interaction with TBX5. Nature 424 443 447 12845333 SchleiffarthJRPersonADMartinsenBJSukovichDJNeumannA 2007 Wnt5a is required for cardiac outflow tract septation in mice. Pediatr Res 61 386 391 17515859 vanREOlsonEN 2007 MicroRNAs: powerful new regulators of heart disease and provocative therapeutic targets. J Clin Invest 117 2369 2376 17786230 NyirendaMJDeanSLyonsVChapmanKESecklJR 2006 Prenatal programming of hepatocyte nuclear factor 4alpha in the rat: A key mechanism in the ‘foetal origins of hyperglycaemia’? Diabetologia 49 1412 1420 16570165 FilkowskiJNIlnytskyyYTammingaJKoturbashIGolubovA 2010 Hypomethylation and genome instability in the germline of exposed parents and their progeny is associated with altered miRNA expression. Carcinogenesis 31 1110 1115 19959559 MirabelloLSavageSAKordeLGadallaSMGreeneMH 2010 LINE-1 methylation is inherited in familial testicular cancer kindreds. BMC Med Genet 11 77 20478068 KileMLBaccarelliATarantiniLHoffmanEWrightRO 2010 Correlation of global and gene-specific DNA methylation in maternal-infant pairs. PLoS One 5 e13730 21060777 FeinbergAP 2010 Genome-scale approaches to the epigenetics of common human disease. Virchows Arch 456 13 21 19844740 MorganDKWhitelawE 2008 The case for transgenerational epigenetic inheritance in humans. Mamm Genome 19 394 397 18663528 SrivastavaD 2006 Making or breaking the heart: from lineage determination to morphogenesis. Cell 126 1037 1048 16990131 BruneauBG 2008 The developmental genetics of congenital heart disease. Nature 451 943 948 18288184 ZhuHCabreraRMWlodarczykBJBozinovDWangD 2007 Differentially expressed genes in embryonic cardiac tissues of mice lacking Folr1 gene activity. BMC Dev Biol 7 128 18028541 BottiniECosmiENicotraMSanteusanioGLaTM 1998 The genetics of signal transduction and the feto-maternal relationship. A study of cytosolic low molecular weight phosphotyrosine phosphatase. Dis Markers 14 143 150 10427472 KajimotoHHashimotoKBonnetSNHaromyAHarryG 2007 Oxygen activates the Rho/Rho-kinase pathway and induces RhoB and ROCK-1 expression in human and rabbit ductus arteriosus by increasing mitochondria-derived reactive oxygen species: a newly recognized mechanism for sustaining ductal constriction. Circulation 115 1777 1788 17353442 KaynakBvonHAMebusSSeelowDHennigS 2003 Genome-wide array analysis of normal and malformed human hearts. Circulation 107 2467 2474 12742993 PhillipsHMRheeHJMurdochJNHildrethVPeatJD 2007 Disruption of planar cell polarity signaling results in congenital heart defects and cardiomyopathy attributable to early cardiomyocyte disorganization. Circ Res 101 137 145 17556662 QingMSchumacherKHeiseRWoltjeMVazquez-JimenezJF 2003 Intramyocardial synthesis of pro- and anti-inflammatory cytokines in infants with congenital cardiac defects. J Am Coll Cardiol 41 2266 2274 12821258 DolinoyDCHuangDJirtleRL 2007 Maternal nutrient supplementation counteracts bisphenol A-induced DNA hypomethylation in early development. Proc Natl Acad Sci U S A 104 13056 13061 17670942 BjornssonHTSigurdssonMIFallinMDIrizarryRAAspelundT 2008 Intra-individual change over time in DNA methylation with familial clustering. JAMA 299 2877 2883 18577732 EckhardtFLewinJCorteseRRakyanVKAttwoodJ 2006 DNA methylation profiling of human chromosomes 6, 20 and 22. Nat Genet 38 1378 1385 17072317 WangLWangFGuanJLeJWuL 2010 Relation between hypomethylation of long interspersed nucleotide elements and risk of neural tube defects. Am J Clin Nutr 91 1359 1367 20164316 ChenXGuoJLeiYZouJLuX 2010 Global DNA hypomethylation is associated with NTD-affected pregnancy: A case-control study. Birth Defects Res A Clin Mol Teratol 88 575 581 20641100 BibikovaMLeJBarnesBSaedinia-MelynkSZhouLShenR Genome-wide DNA methylation profiling using Infinium assay. Epigenomics 1[1], 177-200. 2009 YoonPWRasmussenSALynbergMCMooreCAAnderkaM 2001 The National Birth Defects Prevention Study. Public Health Rep 116 Suppl 1 32 40 11889273 BottoLDLinAERiehle-ColarussoTMalikSCorreaA 2007 Seeking causes: Classifying and evaluating congenital heart defects in etiologic studies. Birth Defects Res A Clin Mol Teratol 79 714 727 17729292 BlockGHartmanAMNaughtonD 1990 A reduced dietary questionnaire: development and validation. Epidemiology 1 58 64 2081241 TroyanskayaOCantorMSherlockGBrownPHastieT 2001 Missing value estimation methods for DNA microarrays. Bioinformatics 17 520 525 11395428 StoreyJDTaylorJESigmundD Strong control, conservative point estimation, and simultaneous conservative consistency of false discovery rates: A unified approach. Journal of the Royal Statistical Society, Series B 66 187 205 2004 R Development Core Team R: A language and environment for statistical computing. R Foundation for Statistical Computing 2008 SubramanianATamayoPMoothaVKMukherjeeSEbertBL 2005 Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A 102 15545 15550 16199517