- Open Access
Genome-scale hypomethylation in the cord blood DNAs associated with early onset preeclampsia
Clinical Epigeneticsvolume 7, Article number: 21 (2015)
Preeclampsia is one of the leading causes of fetal and maternal morbidity and mortality worldwide. Preterm babies of mothers with early onset preeclampsia (EOPE) are at higher risks for various diseases later on in life, including cardiovascular diseases. We hypothesized that genome-wide epigenetic alterations occur in cord blood DNAs in association with EOPE and conducted a case control study to compare the genome-scale methylome differences in cord blood DNAs between 12 EOPE-associated and 8 normal births.
Bioinformatics analysis of methylation data from the Infinium HumanMethylation450 BeadChip shows a genome-scale hypomethylation pattern in EOPE, with 51,486 hypomethylated CpG sites and 12,563 hypermethylated sites (adjusted P <0.05). A similar trend also exists in the proximal promoters (TSS200) associated with protein-coding genes. Using summary statistics on the CpG sites in TSS200 regions, promoters of 643 and 389 genes are hypomethylated and hypermethylated, respectively. Promoter-based differential methylation (DM) analysis reveals that genes in the farnesoid X receptor and liver X receptor (FXR/LXR) pathway are enriched, indicating dysfunction of lipid metabolism in cord blood cells. Additional biological functional alterations involve inflammation, cell growth, and hematological system development. A two-way ANOVA analysis among coupled cord blood and amniotic membrane samples shows that a group of genes involved in inflammation, lipid metabolism, and proliferation are persistently differentially methylated in both tissues, including IL12B, FAS, PIK31, and IGF1.
These findings provide, for the first time, evidence of prominent genome-scale DNA methylation modifications in cord blood DNAs associated with EOPE. They may suggest a connection between inflammation and lipid dysregulation in EOPE-associated newborns and a higher risk of cardiovascular diseases later in adulthood.
Preeclampsia is a pregnancy-specific disease characterized by hypertension (blood pressure >140/90 mmHg) and proteinuria (>0.3 g/day). It affects 5% to 8% of all pregnancies  and is one of the leading causes of fetal/neonatal and maternal morbidity and mortality . Preeclampsia may develop at any time after 20 weeks of pregnancy, and depending on the time of the onset of symptoms (earlier or later than 34 weeks), it can be categorized as early onset preeclampsia (EOPE) or late onset preeclampsia (LOPE). Although EOPE and LOPE share some clinical symptoms, they have distinct etiologies, biochemical markers, and maternal and fetal outcomes . EOPE is considered a more severe, mostly fetal disorder, and it is typically associated with inadequate placental implantation and subsequent placental dysfunction, intrauterine growth restriction, low birth weight, and adverse maternal and fetal outcomes, including perinatal death. By contrast, LOPE is considered a maternal disorder often associated with a normal placenta, larger placental volume, normal fetal growth, and favorable maternal and fetal outcomes. In this study, we focus on EOPE given its severity.
The etiology of EOPE is unknown. It is hypothesized that the process is initiated by the shallow implantation of the placenta into the uterine wall . During normal pregnancy, extravillous trophoblasts remodel the spiral arteries in the decidua lining during the formation of the maternal placenta. However, in EOPE, such remodeling is incomplete, leading to poor placental development and function . The incompletely transformed spiral arteries are unable to provide an adequate blood supply to the placenta, resulting in ischemia and hypoxic conditions, which further lead to inflammatory stress, apoptosis, and other phenomena in the trophoblasts [6-8]. In mothers, the pathological manifestations of these conditions include hypertension, thrombocytopenia, hemolysis, seizures, acute atherosis, renal failure, and proteinuria [9-11]. Whereas in the fetus, the consequences include intrauterine growth restriction and preterm delivery with associated issues such as neurodevelopmental disorders  and long-term adult onset disorders including diabetes, congestive heart failure, and hypertension .
Increasing evidence shows that epigenetic processes, especially DNA methylation, play important roles in preeclampsia . While the majority of epigenetic studies have been focusing on detecting differentially methylated regions in the maternal blood [15,16] or placental tissues [17-20], only a few studies have investigated the alteration of local or global methylation and epigenetic factors in neonatal cord blood cells associated with preeclampsia. At the individual gene level, IGF2 but not GNAS has been found to be significantly hypomethylated in cord blood cells associated with maternal preeclampsia . At the global level, controversy exists regarding the plausible relationship between DNA methylation and preterm birth, a major consequence of maternal EOPE. While some authors have proposed that preeclampsia might epigenetically program placental tissue but not umbilical cord blood , others have found that preterm birth is associated with overall lower LINE-1 methylation levels in cord blood, which could be a surrogate of global DNA methylation levels . However, to our knowledge, there has been no epigenome-scale methylation study on cord blood DNAs in association with EOPE. In this study, we have for the first time used the Infinium Human Methylation 450 K BeadChip (450 K) to quantitatively analyze methylation on over 485,000 CpG sites in the cord blood DNAs from EOPE mothers (cases) and full-term controls. Our bioinformatics and systems biology analyses have revealed a global hypomethylation pattern in DNA extracted from cord blood in infants born at <34 weeks due to EOPE and have found many changes related to immune response and cardiovascular dysfunction. The discovery could shed light on the connection to adult-onset diseases in the newborns from EOPE mothers.
Global CpG hypomethylation in cord blood DNAs associated with maternal EOPE
We compared the genome-scale methylation profile of 20 cord blood DNAs associated with maternal EOPE cases and normal term delivery controls (12 cases vs. 8 controls). A summary of physiological and clinical factors associated with the maternal samples is tabulated in Table 1, including mother’s age, gestational age, gravida, parity, baby gender, and birth weight. We did not observe statistical differences in maternal age, gravida, parity, and sex of babies between cases and controls. As a consequence of EOPE, the patients had to deliver the babies prematurely to terminate the potentially fatal EOPE condition. Thus, both gestational age and mean weight of the babies were significantly different between EOPE and full-term controls (31.5 weeks versus 39.3 weeks, P = 2.28e − 8 and 1,350 g versus P = 1.17e − 9 and 3,420 g, respectively).
We first performed Principal Component Analysis (PCA) with the M-values of the 384,632 pre-processed methylation data points (Figure 1A). The cases were separable from the controls with the first three principal components, indicating that the cord blood cell methylation data provide strong signatures to the maternal preeclampsia condition. To systematically determine how the physiological and clinical factors contributed to variance in the data, we performed analysis of variance (ANOVA) and computed the average F-statistic for each factor (Figure 1B). Given that gestational age and baby weight differences are the direct consequences of EOPE and thus they are highly correlated to EOPE, we use the term “preeclampsia” to represent the overall phenotype of these three factors. ANOVA analysis shows that preeclampsia condition is the highest source of variance, followed by parity and gravida. To determine the best statistical model for the differential methylation (DM) test, we used Akaike’s information criterion (AIC) to compare the effects of physiological and clinical factors on the linear model. The result shows that the model with just preeclampsia had the best score and was selected more often than the model containing parity and/or gravida. We therefore conclude that maternal preeclampsia is the primary source for the differences in methylation among samples.
The DM analysis using a moderated t-test identified a total of 64,049 significant CpG sites after the Benjamini-Hochberg (BH) multiple hypothesis testing correction with a threshold of P <0.05. A global hypomethylation pattern was observed, with 51,486 hypomethylated sites and 12,563 hypermethylated sites (Figure 1C). To validate the results of the Illumina 450 K array, we performed pyrosequencing on seven selected CpG methylation sites (cg26478401 from FAS, cg2607695 from PI3KR1, cg06111286 from IL12B, cg06652329 from IGF1, cg08198483 from miR-24-2, cg11671363 and cg08537847 from miR-145). A coefficient of correlation (R 2) as high as 0.9387 (P < 2.2e − 6) was obtained by comparing the methylation percentages determined by pyrosequencing and the β-values determined by Illumina 450 K (Figure 2A). For individual CpG sites, the directions of fold change in the 450 K compared to those of the pyrosequencing were consistent (Figure 2B, C, D, E, F, G, H), and the differences in pyrosequencing results between the cases and controls were significant. These results show the robustness and validity of the 450 K data.
Genome-scale patterns of differentially methylated CpG sites in CpG islands and surrounding areas
To examine the methylation profile across genomic locations, we counted the number of DM sites, hypermethylated sites, and hypomethylated sites in different genomic categories relative to CpG islands (Figure 3A). The CpG-related categories are based on the annotation given by Illumina. CpG islands are defined as regions with a GC content of 50% or greater, length greater than 200 bp and with a ratio of observed CpGs to the expected number of CpGs greater than 0.6 . CpG shores are defined as 0 to 2 kb upstream (north) or downstream (south) of the CpG islands and CpG shelves as 2 to 4 kb upstream or downstream of the CpG islands. The open sea regions refer to regions not located in any of the abovementioned categories relating to CpG islands.
Significant differences between observed frequencies and expected frequencies exist in all categories based on CpG-related annotations as stated in the previous paragraph (Figure 3A). The most striking result is the significant (about fourfold) reduction in the percentages of the overall DM and hypermethylated and hypomethylated sites in the CpG islands compared to the background percentage of probes in CpG islands (for example, the percentage of all annotated probes in CpG islands). The hypomethylation pattern is also evident in the heatmap results of CpG islands with unsupervised hierarchical clustering, where the case and control groups are clearly separated (Figure 4A). Conversely, there is significantly more DM, including both hypermethylation and hypomethylation in the open sea region compared to the background percentage (Figure 3A). The majority of CpG islands are associated with gene promoters whereas the open seas are often associated with gene body and intergenic regions. These results show that methylation changes in non-promoter regions of cord blood DNAs are mostly associated with maternal EOPE (Figure 3B).
As examples, we selected the most statistically significant 20 hypomethylated and 20 hypermethylated CpG sites with an average β-value difference greater than 0.2 and present their dot plots in Additional file 1: Figure S1 and Additional file 2: Figure S2. The top three hypermethylated loci were cg03110921 (associated with DLG5), cg20968678 (associated with TMEM169), and cg14087413 (associated with TMEM8B). The top three hypomethylated loci were cg18188328 (associated with JARID2), cg02522196 (not associated with genes), and cg02291365 (not associated with genes). The dot plots have consistent differential methylation patterns in these CpG sites between preeclampsia cases and controls.
Genome-scale patterns of differentially methylated CpG sites in association with genes
According to the Illumina 450 K annotation, we categorized the CpG sites in association with genes into seven groups, namely the TSS1500 (200 to 1,500 bp upstream of transcription starting sites), TSS200 (the transcription starting sites to 200 bp upstream interval), 5′ UTR, first exon, gene body, 3′ UTR, and intergenic regions. Significant differences between observed and background frequencies are seen in all the seven CpG categories in association with genes (Figure 3B). As expected from the annotations by CpG island and surrounding areas, the most striking results are the significantly higher percentages of overall DM, hypermethylated and hypomethylated sites in the gene body and intergenic regions, compared to the background percentage of probes located in these regions. Similar results are also observed in the 3′ UTR. These results highlight the significance of non-promoter-associated differential methylation in EOPE. Consistent with the observations in the CpG islands, we find that there is a significantly lower percentage of DM and hypermethylated and hypomethylated sites in the TSS200 and first exon regions compared to their background probe percentages. The hypomethylation pattern in the TSS200 region is also evident in the heatmap results of CpG islands with unsupervised hierarchical clustering, where the case and control groups are clearly separated (Figure 4B). On the other hand, the TSS1500 regions have a significantly lower percentage of hypomethylated sites but about the same percentage of hypermethylated sites compared to background percentages.
Promoter-associated network and pathway analysis
Given that methylation of CpG islands and promoter regions have received the most attention in determining how epigenetic changes affect gene expression [25,26], we focused the rest of the analysis on the TSS200 region. For each gene, we took the sum of the Beta values and performed a DM test on the summed Beta values, using a moderated t-test. As a result, 643 genes with hypomethylated promoters and 389 genes with hypermethylated promoters were identified.
We then performed an Ingenuity Pathway Analysis (IPA) on the genes with significant DM at TSS200 to identify pathways and networks that were altered in EOPE (Figure 5A, C, D, Tables 2 and 3). The top five networks were as follows: 1) inflammatory response, cellular function and maintenance, and respiratory disease; 2) cellular development, cellular growth and proliferation, and hematological system development and function; 3) cell-to-cell signaling and interaction, hematological system development and function, and immune cell trafficking; 4) molecular transport, cellular movement, and hematological system development and function; and 5) cellular development, cellular growth and proliferation, and hematological system development and function. Many genes previously shown to be important in EOPE were found in these top networks, such as PI3K1, FAS, IGF1, and IL12B. In addition, several microRNAs appeared in the top five networks, such as mir-10a, mir-181c, mir-21, mir-101-1, and mir-29a. Interestingly, the top canonical pathways were as follows: (1) FXR/RXR activation, (2) granulocyte adhesion and diapedesis, (3) altered T-cell and B-cell signaling in rheumatoid arthritis, (4) hepatic cholestasis, and (5) agranulocyte adhesion and diapedesis.
We next used the software package EpiMod to validate the network findings from IPA. Previously, EpiMod was used to identify DM hotspots in the Illumina Infinium 27 K platform, by integrating methylation data with protein interaction networks using the methylation of promoter-associated CpG sites . We have extended the package to include promoter-associated probes from Illumina 450 K. Similarly to the IPA analysis, we used the sum of Beta values within the TSS200 as the measure for DM, in order to perform the EpiMod analysis. Many key genes in the EpiMod networks were also hubs in the IPA network analysis, and the networks themselves were associated with many of the same functions. We show one exemplary module from EpiMods with the hotspot genes IL12B and receptors IL12RB1/B2 (Figure 5B). Similar to the best ranked network in IPA (Figure 5A), this epi-module shows a clear enrichment of genes involved in the inflammatory response. Consistently, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis on the genes found in this epi-module shows that “Cytokine-cytokine receptor interaction” is the most significant pathway (FDR = 2.5e − 21).
Differential methylation in the non-coding genomic regions
Since more than 95% of the genome does not encode proteins, we next examined the CpG sites in non-coding RNAs. We first analyzed microRNAs, given that they are relatively better understood compared to other non-coding RNAs (ncRNAs). The microRNAs present in the top ten IPA networks were grouped based on their methylation status (hypo- vs. hypermethylation), and their targets were predicted with the TargetScan program. This list of target genes was then subjected to KEGG analysis using the DAVID online functional annotation tool . Pathways with a BH-adjusted P value less than 0.05 were considered significant. The pathways that were enriched in the hypomethylated and hypermethylated microRNA target gene sets were similar, with common signaling pathways such as the “ErbB signaling pathway” and the “neurotrophin signaling pathway” (Table 4). Other common biological processes such as “endocytosis,” as well as diseases such as “Chronic myeloid leukemia” were also enriched.
We next predicted the function of various long intergenic non-coding RNAs (lincRNAs) that had DM in TSS200 regions. In order to represent the protein-coding genes with significant DM at the pathway level, the summed β-values of the TSS200 regions of DM protein-coding genes were used to generate a pathway-based deregulation score matrix (PDS matrix). This matrix is a measure of how much the activities of various pathways deviate from normal. The summed β-values of the TSS200 regions of DM lincRNAs were then correlated with the scores in the PDS matrix using the Spearman correlation coefficient (ρ). Only lincRNA-pathway pairs that satisfied the cutoff absolute (ρ) >0.90 and P <0.05 are considered significant. A total of 96 KEGG and BIOCARTA pathways stand out, many of which are related to the immune response and inflammation [see Additional file 3: Table S2].
Tissue differences and genes of interest
We previously reported a global hypermethylation coupled with hypomethylation in promoters from the chorioamniotic membrane of EOPE patients . Some of the patients who participated in that study gave cord blood samples for this study. To determine the contribution of tissue specificity vs. disease (EOPE) specificity, we compared the subsets of methylation data from the coupled chorioamniotic membrane and cord blood samples. An ANOVA model containing terms for preeclampsia and the tissue type was used to gain insight on their contributions. The resulting F-statistics show that tissue type rather than EOPE is the dominant factor responsible for the variance among the methylation data [see Additional file 4: Figure S3].
We further analyzed the DM sites in both chorioamniotic membranes and the cord blood, with moderated t-tests subject to BH multiple hypothesis testing (MHT). The analysis revealed that 10,510 CpG sites were considered DM, and among them, 10,355 CpG sites were consistently hypermethylated or hypomethylated in both tissues. These sites were mapped to the TSS200 regions of 579 genes. Multiple genes from the IPA and EpiMod networks were found among these 579 genes, including the highly connected genes CD86, IL19, OSCAR, TRAF3IP2, FAS, PIK3R1, IFNG, IL1RN, and ISG20. These results suggest that a small number of genes involved in inflammation and lipid metabolism are epigenetically modified similarly in both membrane tissues and cord blood of fetal origin.
Comparison with other cord blood methylation studies of preeclampsia
Compared to the number of methylation studies in placenta and maternal peripheral blood tissues related to EOPE, there are surprisingly few reports on umbilical cord blood DNAs. This might be partially due to the fact that EOPE is a maternal disease. However, given the well-established link between preterm birth and adulthood onset diseases such as hypertension, cardiovascular diseases, and obesity, it is plausible that some long-term consequences are carried through the changes made by “epigenetics markers” at or even before birth. In this study, we describe in depth the significant genome-scale hypomethylation observed in cord blood DNAs associated with maternal EOPE, based on the Illumina 450 K. Our study confirms and expands on the very few earlier methylation studies of cord blood samples from preterm deliveries (or more specifically, preterm deliveries due to EOPE), which had either lower throughput or lacked gene level resolution for DNA methylation. Burris et al. used LINE-1 methylation levels as a surrogate measurement of the global methylation and found that after adjustment for other confounding factors, a lower cord blood LINE-1 level methylation was associated with a shorter gestation (−0.45 weeks, (95% CI −0.83, −0.08)) and higher odds of preterm birth (OR 4.55 (95% CI 1.18, 17.5)) . On the other hand, He et al. conducted a pilot study to examine DNA methylation levels of IGF2 and GNAS in cord blood associated with preeclampsia. They found that the average methylation level of IGF2 but not GNAS was significantly lower in the preeclampsia samples . We used a much higher throughput method with the Illumina 450 K and also found that IGF2 has one CpG site with a significantly lower methylation level in the gene body (cg01668279, FDR adjusted P = 0.0429). However, unlike the other study, we found that seven out of nine CpG sites associated with GNAS have significantly lower methylation in EOPE, in regions ranging from the TSS1500 to the gene body. The overall agreement between our study and the majority of other studies suggest that EOPE has broad effects resulting in abnormal DNA methylation even in newborn cord blood cells.
Potential complications in the study due to gestational age and birth weight
In our current experimental design of comparing EOPE vs. full-term delivery, earlier gestational age and lower mean birth weights are naturally coupled with EOPE. We considered EOPE as the major contributing factor of DM, because the majority of the best-fit models have EOPE as the only effective predictor. Also, the inclusion of highly correlated factors in the statistical model, such as gestational age and birth weight that are coupled with EOPE, can inflate the standard error  and potentially mask the causative effect of EOPE. Admittedly, these two factors may confound the observed differences due to EOPE. Hence, we compared our results with other studies on the effect of gestational age on placental or cord blood DNA methylation pattern without preeclampsia [30-32]. Studies using cord blood show little evidence of changes in methylation levels within the third trimester [33,34]. In one study, Michels et al. found less than 1% of effect comparing cord blood methylation of preterm babies to full-term babies . In another study, Nomura et al. found no significant effect on methylation due to gestational age, in both placenta and cord blood (P = 0.31 and P = 0.82, respectively) . In the placenta, the DNA methylation pattern does change as gestation progresses. However, the major changes of methylation occur from the first to the second trimesters and considerably fewer changes happen after the second trimester . Moreover, within the third trimester, the placental methylation patterns from babies delivered earlier (around 34 weeks) and those of full-term babies are clustered together and they cannot be distinguished. This suggests that the gestational ages within the last trimester have much less effect on the placental DNA methylation. Novakovic et al. also stated that the variation among samples in the third trimester was primarily due to environmentally induced or stochastic reasons. Altogether, these studies suggest that very small or minimal differences in the observed methylation patterns are due to gestational age in our study. This is confirmed by Cruickshank et al. who found that only 1,555 DM sites from neonatal blood spots were statistically significantly associated with gestational age , whereas we found 64,049 DM sites associated with EOPE after MHT correction. Further, we compared the DM sites in both studies and found that only 319 sites are commonly hypomethylated and 139 sites are commonly hypermethylated [see Additional file 5: Table S3]. In summary, we conclude that gestational age and birth weight have insignificant effect compared to the disease of EOPE.
Methylation changes related to inflammation
Evidence of changes in CpG methylation related to the inflammatory process is strong in our data. In the top-ranked IPA network (Figure 4A) and EpiMod module (Figure 4B), we see that a large number of DM genes (for example, IL12B, IL19, IL20, CCL2) are connected to the nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) complex, although NF-κB itself is not differentially methylated. Although it is well known that increased inflammatory stress exists in preeclampsia placentas [35,36] and maternal tissues , little is known about immunological responses in newborns. Some studies have shown that Interleukin-6 (IL6) was significantly elevated in the cord blood from neonates born small with intrauterine growth retardation. However, no direct measurements of IL6 were made from those born to preeclamptic mothers . In another study, the authors found evidence of higher NKp30 levels in NK cells, but not in the monocytes of cord blood from preeclamptic pregnancies, suggesting that inflammatory responses may be cell-type specific . Since we used the extracted genomic DNA from cord blood cells without the purification of a subpopulation, the subpopulation-specific inflammation response in EOPE could not be studied. However, we do know that the DNA methylation changes occur in many genes involved in inflammation. Meanwhile, we also observe that some methylation differences are very interesting. For example, the hyper-methylation of IL12B is validated by our pyrosequencing, although the Beta value only decreases by 0.1 from 0.3 in cases (P < 0.05) (Figure 2E). Corresponding to this finding, others reported that preeclampsia patients had significantly less IL-12 in villous trophoblasts , as well as significantly less IL-12 secretion from the decidua . Moreover, a positive correlation between IL-12 and fetal birth weight in the severe preeclamptic group was observed . All the information seems to support the hypothesis that IL12B is repressed at the gene expression level in the cord blood cells. It will be very interesting to validate the IL12B gene expression in fresh cord blood tissue.
Methylation changes related to potential cardiovascular problems
Interestingly, we found that FXR/RXR activation and hepatic cholestasis are among the top enriched pathways in the IPA analysis of TSS200 DM regions. An important target of FXR/RXR is fatty acid synthase, FAS , which was consistently hypermethylated in both cord blood and the placental chorioamniotic membrane. We propose two hypotheses that could explain the methylation changes occurring in these pathways. One hypothesis arises from the fact that fetal liver is the main source of hemopoietic cells during the majority of the gestational period ; thus, the DM related to FXR/RXR activation and cholestasis in the cord blood cells likely reflects the dysfunction of liver tissues of newborns associated with maternal EOPE. The nuclear receptors farnesoid X receptor (FXR) and liver X receptor (LXR) form heterodimers with the retinoid X receptors (RXRs) to regulate cholesterol and lipid metabolism, and the primary tissue responsible for cholesterol metabolism is the liver, where cholesterol biosynthesis occurs and cholesterol low-density lipoprotein (LDL) is taken up from the plasma by the LDL receptor. Another hypothesis is that the DM observed in the FXR/RXR pathway from the cord blood cells may be attributed to circulating macrophages within the cord blood cells. LXR/RXR dimers in macrophages activate ABC transporters, such as ABCA1 and ABCG1, which efflux cholesterol from macrophages and prevent toxicity due to cholesterol accumulation . In adults, macrophages overloaded with cholesterol (or foam cells), are hallmarks of atherosclerosis . Regardless of the mechanism leading to the DM among genes related to FXR/RXR activation in cord blood cells, the FXR/RXR pathway is known to play an important role in maintaining cholesterol homeostasis. Abnormal methylation in this pathway might be indicative of long-term risks for coronary heart disease due to lipid metabolism dysfunction. Other processes, such as the inflammatory stress mentioned earlier, might also mediate or contribute to higher risks of cardiovascular diseases later in life for babies born to preeclamptic mothers.
The samples in this study were obtained from the pre-processed genomic DNA of cord blood cells stored in the biorepository. As no cord blood serum was stored, we could not obtain inflammatory molecule signals or the lipid profiles in the cord blood. We plan to follow up with another human study to validate the candidate biomarker findings in the protein-coding genes and noncoding RNAs at the methylation level and the expression level.
To our knowledge, this study is the first to address the genome-scale methylation changes in the cord blood cells associated with EOPE patients. The newborn cord blood cells are thought to originate from the fetal liver; thus, much of the methylation changes involve signatures of liver functions. We have found significant genome-scale DNA methylation modifications that may suggest a link between inflammation and lipid dysregulation in newborns. These findings may suggest a link between maternal EOPE and newborns with an increased lifetime risk of cardiovascular diseases.
We utilized samples from the RMATRIX Hawaii Biorepository (HiBR), which has obtained its own IRB approval. From this repository, a total of 21 cord blood samples were initially identified. Maternal cases had EOPE at <34 weeks gestation and controls were full-term pregnancies without preeclampsia. To ensure the quality of the cord blood samples, maternal patients with known medical or obstetrical disorders associated with preeclampsia were excluded, including diabetes mellitus or gestational diabetes mellitus, hypertension, renal disease, hyperthyroidism, and systemic lupus erythematous. Only singleton pregnancies were included.
Genomic DNA from whole cord blood cells was extracted by the HiBR laboratory. Clotted cord blood was centrifuged through a Clotspin Basket (Qiagen) at 2,000 × g for 5 min to disperse the clot. Three volumes of Red Blood Cell (RBC) Lysis Solution (Qiagen) were added to one volume of dispersed cord blood. The samples were briefly vortexed, shaken for 15 min at room temperature then centrifuged as above to pellet leukocytes and clot particulates. The pellet was resuspended in Cell Lysis Solution (Qiagen) containing Proteinase K (Qiagen) at 20 mg/ml and then incubated at 55°C overnight. The Gentra Systems Autopure LS (Qiagen) was used to extract the DNA from the sample using the Compromised Cell Lysate Protocol. The purified DNA solution was then incubated in a 65°C water bath for 1 h and shaken at room temperature overnight.
Genome-scale methylome profiling using the HumanMethylation450 BeadChip
Methylation profiling was done using the Illumina Infinium methylation assay and the HumanMethylation450 BeadChip. First, the cord blood DNA samples were pre-screened for the suitability of the Methylation450 assay using Nanodrop to check the purity and the PicoGreen assay to check the dsDNA concentration. The EZ DNA Methylation kit was used for the bisulfite conversion (Zymo Research), and 100 to 250 ng of bisulfite-converted DNA was hybridized onto the Infinium Human Methylation450 BeadChip following the Illumina Infinium HD Methylation protocol.
Primer design was carried out using the PyroMark Assay Design 2.0 software. One of the primers was biotinylated to enable capture by Streptavidin Sepharose. Bisulfite-treated DNA was amplified followed by pyrosequencing using the PyroMark PCR kit (Qiagen) on a PyroMark Q24 instrument using PyroMark Gold Q24 Reagents, according to the manufacturer’s protocol. The primers for the seven CpG sites selected for validation are listed in Additional file 6: Table S1.
Methylation data preprocessing
Raw data preprocessing was performed using the R package “minfi” . The background signal was subtracted from the raw intensity data, and any probes with a detection P value greater than 0.05 were eliminated. The intensity value for each probe was transformed first to methylation Beta values and then to M-values (the logit transform of Beta values) and normalized using the SWAN algorithm . Probes with known SNPs based on the latest annotation by Illumina were removed from downstream analysis due to the concern that they could confound intensity changes during later analysis steps. Probes for CpG sites on sex chromosomes were removed to eliminate the effect of the genders of newborns.
Site-level differential methylation (DM) analysis
Before DM analysis, samples were examined for the differences using the primary physiological or clinical parameters. As a result, 12 cases (one case was eliminated due to the outlier birth weight) and eight controls were determined suitable for DM. The AIC was used to select the primary physiological or clinical parameters in a linear model that best explains the differences among the methylation data. The model that optimized the AIC score most frequently was the model with “preeclampsia” as the only contributing variable. The importance of this term was verified using ANOVA analysis, where the preeclampsia factor was found to have the largest value for its average F-statistic. A moderated t-test from the “limma” package was used to determine the statistical significance of difference in methylation for each CpG site between cases and controls. The resulting P value was then adjusted for MHT using the BH correction. CpG sites with an adjusted P value less than 0.05 were considered significant. Correlation heatmaps, PCA plots, volcano plots, and bar graphs, and so on were generated using R 3.0.2.
Ingenuity pathway analysis (IPA)
Since the promoter region plays a critical role in transcriptional regulation, the sum of the Beta values of CpG sites in the TSS200 region was used as a proxy for gene methylation levels. The “limma” package’s moderated t-test was used to determine the statistical significance of methylation differences between cases and controls for each gene. Similar to the site-level analysis, P values were corrected for MHT using the BH correction, and adjusted P values less than 0.05 were considered significant. The log2-fold change in methylation values was calculated, and values for significantly differentially methylated genes were used as inputs to IPA to generate gene interaction networks and canonical pathways that are significantly associated with EOPE.
Similar to the IPA analysis, the sums of Beta values for sites in TSS200 regions were used as inputs for the functions in the EpiMod package . The EpiMod package was originally developed to analyze Illumina 27 K data. We modified the code to analyze 450 K data by replacing the methylation values of each gene with the sum of TSS200, similar to the newer package Functional Epigenetic Modules (FEM) from the same developers of EpiMods .
Non-coding RNA methylation analysis
Two types of ncRNA were analyzed: microRNA and lincRNA. TargetScan was used to predict mRNA targets for the microRNAs that were found in the top ten IPA networks. The microRNA targets were then subjected to GO and KEGG analysis to determine the enriched gene sets among the predicted targets. For lincRNA analysis, the CpG sites were associated with annotated lincRNAs from the UCSC browser, and lincRNAs were tested for differential methylation analysis in the TSS200 regions, similar to the approach used for the protein-coding genes described earlier (see the “Methods” section for IPA analysis). To predict the function of differentially methylated lincRNAs, the “pathifier” package was used to create a PDS matrix , using the DM protein-coding RefSeq genes per their sum of β-values in the TSS200 region. Spearman’s rho (ρ) correlation coefficient was used to calculate the correlation between the lincRNA methylation and the pathways built from DM RefSeq genes in the samples. Pathway-lincRNA pairs with |ρ| > 0.90 were kept for functional analysis.
Akaike’s information criterion
Benjamini-Hochberg (BH) correction, analysis of variance
differential methylation or differentially methylated
early onset preeclampsia
RMATRIX Hawaii Biorepository
ingenuity pathway analysis
late onset preeclampsia
multiple hypothesis testing
Srinivas SK, Edlow AG, Neff PM, Sammel MD, Andrela CM, Elovitz MA. Rethinking IUGR in preeclampsia: dependent or independent of maternal hypertension? J Perinatol. 2009;29(10):680–4.
Say L, Chou D, Gemmill A, Tunçalp Ö, Moller A-B, Daniels J, et al. Global causes of maternal death: a WHO systematic analysis. Lancet Global Health. 2014;2(6):e323–33.
Raymond D, Peterson E. A critical review of early-onset and late-onset preeclampsia. Obstet Gynecol Surv. 2011;66(8):497–506.
Jung J-J, Tiwari A, Inamdar SM, Thomas CP, Goel A, Choudhury A. Secretion of soluble vascular endothelial growth factor receptor 1 (sVEGFR1/sFlt1) requires Arf1, Arf6, and Rab11 GTPases. PLoS One. 2012;7(9):e44572.
Pijnenborg R, Vercruysse L, Hanssens M. The uterine spiral arteries in human pregnancy: facts and controversies. Placenta. 2006;27(9–10):939–58.
Allaire AD, Ballenger KA, Wells SR, McMahon MJ, Lessey BA. Placental apoptosis in preeclampsia. Obstet Gynecol. 2000;96(2):271–6.
Leung DN, Smith SC, To KF, Sahota DS, Baker PN. Increased placental apoptosis in pregnancies complicated by preeclampsia. Am J Obstet Gynecol. 2001;184(6):1249–50.
Tal R. The role of hypoxia and hypoxia-inducible factor-1alpha in preeclampsia pathogenesis. Biol Reprod. 2012;87(6):134.
Boos CJ, Lip GYH. Is hypertension an inflammatory process? Curr Pharm Des. 2006;12(13):1623–35.
Nizet V, Johnson RS. Interdependence of hypoxic and innate immune responses. Nat Rev Immunol. 2009;9(9):609–17.
Palazón A, Aragonés J, Morales-Kastresana A, de Landázuri MO, Melero I. Molecular pathways: hypoxia response in immune cells fighting or promoting cancer. Clin Cancer Res. 2012;18(5):1207–13.
Batalle D, Eixarch E, Figueras F, Muñoz-Moreno E, Bargallo N, Illa M, et al. Altered small-world topology of structural brain networks in infants with intrauterine growth restriction and its association with later neurodevelopmental outcome. Neuroimage. 2012;60(2):1352–66.
Cosmi E, Fanelli T, Visentin S, Trevisanuto D, Zanardo V. Consequences in infants that were intrauterine growth restricted. J Pregnancy. 2011;2011:364381.
Choudhury M, Friedman JE. Epigenetics and microRNAs in preeclampsia. Clin Exp Hypertens. 2012;34(5):334–41.
Mousa AA, Archer KJ, Cappello R, Estrada-Gutierrez G, Isaacs CR, Strauss 3rd JF, et al. DNA methylation is altered in maternal blood vessels of women with preeclampsia. Reprod Sci. 2012;19(12):1332–42.
Anderson CM, Ralph JL, Wright ML, Linggi B, Ohm JE. DNA methylation as a biomarker for preeclampsia. Biol Res Nurs. 2014;16(4):409–20.
Blair JD, Yuen RKC, Lim BK, McFadden DE. Dadelszen Pv, Robinson WP: Widespread DNA hypomethylation at gene enhancer regions in placentas associated with early-onset preeclampsia. Mol Hum Reprod. 2013;19(10):697–708.
Bourque DK, Avila L, Penaherrera M, Von Dadelszen P, Robinson WP. Decreased placental methylation at the H19/IGF2 imprinting control region is associated with normotensive intrauterine growth restriction but not preeclampsia. Placenta. 2010;31(3):197–202.
Kulkarni A, Chavan-Gautam P, Mehendale S, Yadav H, Joshi S. Global DNA methylation patterns in placenta and its association with maternal hypertension in pre-eclampsia. DNA Cell Biol. 2011;30(2):79–84.
Yuen RKC, Penaherrera MS, von Dadelszen P, McFadden DE, Robinson WP. DNA methylation profiling of human placentas reveals promoter hypomethylation of multiple genes in early-onset preeclampsia. Eur J Hum Genet. 2010;18(9):1006–12.
Hart SN, Therneau TM, Zhang Y, Poland GA, Kocher J-P. Calculating sample size estimates for RNA sequencing data. J Comput Biol. 2013;20(12):970–8.
Nomura Y, Lambertini L, Rialdi A, Lee M, Mystal EY, Grabie M, et al. Global methylation in the placenta and umbilical cord blood from pregnancies with maternal gestational diabetes, preeclampsia, and obesity. Reprod Sci. 2014;21(1):131–7.
Burris HH, Rifas-Shiman SL, Baccarelli A, Tarantini L, Boeke CE, Kleinman K, et al. Associations of LINE-1 DNA methylation with preterm birth in a prospective cohort study. J Dev Orig Health Dis. 2012;3(3):173–81.
Gardiner-Garden M, Frommer M. CpG islands in vertebrate genomes. J Mol Biol. 1987;196(2):261–82.
Zöchbauer-Müller S, Fong KM, Virmani AK, Geradts J, Gazdar AF, Minna JD. Aberrant promoter methylation of multiple genes in non-small cell lung cancers. Cancer Res. 2001;61(1):249–55.
West J, Beck S, Wang X, Teschendorff AE. An integrative network algorithm identifies age-associated differential methylation interactome hotspots targeting stem-cell differentiation pathways. Sci Rep. 2013;3:1630.
da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
Ching T, Song MA, Tiirikainen M, Molnar J, Berry M, Towner D, et al. Genome-wide hypermethylation coupled with promoter hypomethylation in the chorioamniotic membranes of early onset pre-eclampsia. Mol Hum Reprod. 2014;20(9):885–904.
Miles J, Shevlin M. Applying regression and correlation: a guide for students and researchers. Sage; 2001:131–148
Bachmayer N, Rafik Hamad R, Liszka L, Bremme K, Sverremark-Ekstrom E. Aberrant uterine natural killer (NK)-cell expression and altered placental and serum levels of the NK-cell promoting cytokine interleukin-12 in pre-eclampsia. Am J Reprod Immunol. 2006;56(5–6):292–301.
Novakovic B, Yuen RK, Gordon L, Penaherrera MS, Sharkey A, Moffett A, et al. Evidence for widespread changes in promoter methylation profile in human placenta in response to increasing gestational age and environmental/stochastic factors. BMC Genomics. 2011;12:529.
Cruickshank MN, Oshlack A, Theda C, Davis PG, Martino D, Sheehan P, et al. Analysis of epigenetic changes in survivors of preterm birth reveals the effect of gestational age and evidence for a long term legacy. Genome Med. 2013;5(10):96.
Michels KB, Harris HR, Barault L. Birthweight, maternal weight trajectories and global DNA methylation of LINE-1 repetitive elements. PLoS One. 2011;6(9):e25254.
Nomura Y, Lambertini L, Rialdi A, Lee M, Mystal EY, Grabie M, et al. Global methylation in the placenta and umbilical cord blood from pregnancies with maternal gestational diabetes, preeclampsia, and obesity. Reprod Sci. 2013;21(1):131–7. doi:10.1177/1933719113492206.
Sikkema JM, Van Rijn BB, Franx A, Bruinse HW, De Roos R, Stroes ESG, et al. Placental superoxide is increased in pre-eclampsia. Placenta. 2001;22(4):304–8.
Dechend R, Viedt C, Müller DN, Ugele B, Brandes RP, Wallukat G, et al. AT1 receptor agonistic antibodies from preeclamptic patients stimulate NADPH oxidase. Circulation. 2003;107(12):1632–9.
Redman CW, Sacks GP, Sargent IL. Preeclampsia: an excessive maternal inflammatory response to pregnancy. Am J Obstet Gynecol. 1999;180(2 Pt 1):499–506.
Lausten-Thomsen U, Olsen M, Greisen G, Schmiegelow K. Inflammatory markers in umbilical cord blood from small-for-gestational-age newborns. Fetal Pediatr Pathol. 2014;33(2):114–8.
Sohlberg E, Saghafian-Hedengren S, Bachmayer N, Hamad RR, Bremme K, Holmlund U. Pre-eclampsia affects cord blood NK cell expression of activation receptors and serum cytokine levels but not CB monocyte characteristics. Am J Reprod Immunol. 2014;71(2):178–88.
Cemgil Arikan D, Aral M, Coskun A, Ozer A. Plasma IL-4, IL-8, IL-12, interferon-gamma and CRP levels in pregnant women with preeclampsia, and their relation with severity of disease and fetal birth weight. J Matern Fetal Neonatal Med. 2012;25(9):1569–73.
Spann NJ, Garmire LX, McDonald JG, Myers DS, Milne SB, Shibata N, et al. Regulated accumulation of desmosterol integrates macrophage lipid metabolism and inflammatory responses. Cell. 2012;151(1):138–52.
Tavassoli M. Embryonic and fetal hemopoiesis: an overview. Blood Cells. 1991;17(2):269–81. discussion 282–266.
Kennedy MA, Barrera GC, Nakamura K, Baldan A, Tarr P, Fishbein MC, et al. ABCG1 has a critical role in mediating cholesterol efflux to HDL and preventing cellular lipid accumulation. Cell Metab. 2005;1(2):121–31.
Glass CK, Witztum JL. Atherosclerosis. The road ahead. Cell. 2001;104(4):503–16.
Aryee MJ, Jaffe AE, Corrada-Bravo H, Ladd-Acosta C, Feinberg AP, Hansen KD, et al. Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics. 2014;30(10):1363–9.
Maksimovic J, Gordon L, Oshlack A. SWAN: Subset-quantile Within Array Normalization for Illumina Infinium HumanMethylation450 BeadChips. Genome Biol. 2012;13(6):R44.
Jiao Y, Widschwendter M, Teschendorff AE. A systems-level integrative framework for genome-wide DNA methylation and gene expression data identifies differential gene expression modules under epigenetic control. Bioinformatics. 2014;30(16):2360–6.
Drier Y, Sheffer M, Domany E. Pathway-based personalized analysis of cancer. Proc Natl Acad Sci U S A. 2013;110(16):6388–93.
This study was partially supported by RMATRIX award U54MD007584 and RCMI BRIDGES award G12 MD007601 from the National Institute on Minority Health and Health Disparities, National Institutes of Health (NIH) and research funding provided by the NIGMS P20 COBRE GM103457, NIEHS K01 ES025434-01, and Hawaiian Community Foundation to L.X. Garmire. M. Berry is supported by NIH grants RO1 DK047320 and G12 MD007601. The Genomics Shared Resource (GSR) at UHCC is supported by the P30-CA071789.
We thank the University of Hawaii Biorepository team, including Dr. Abby Collier, Dr. Timothy Dye, Dr. Joshua Astern, Will Chen, and others for providing samples and resources. The University of Hawaii Biorepository is supported by the National Institute on Minority Health and Health Disparities U54MD007584, G12MD007601, and the National Institute of General Medical Sciences P20GM103466, from the National Institutes of Health.
The authors declare that they have no competing interests.
LXG designed and supervised the data analysis. MB and DT initiated the project and oversaw the sample collection. MT assisted with the sample selection and laboratory procedures. JH and TC analyzed the data and selected the CpG sites and samples for validation. MAS carried out bisulfite conversion for the samples to be utilized in the 450 K array and quality check for the array. MAS performed the 450 K assays and the preliminary QC analysis. JM designed for primers of pyrosequencing and performed the bisulfite conversion and pyrosequencing for the technical validation. JH, TC, and LXG wrote the manuscript. All authors have read, revised, and approved the manuscript.
Travers Ching and James Ha contributed equally to this work.
Top 20 hypomethylated CpG sites. Dot plots of statistically significant top 20 hypomethylated CpG sites with Beta value difference greater than 0.2.
Top 20 hypermethylated CpG sites. Dot plots of statistically significant top 20 hypermethylated CpG sites.
LincRNA pathway analysis. KEGG and BIOCARTA pathways with high correlations (ρ > 0.90) and significant P values (P < 0.05) for lincRNAs.
ANOVA plot of two factor analysis. Average F-statistics plot of tissue and disease factors using the summed TSS200 methylation Beta values in cord blood and paired chorioamniotic membrane samples.
CpG sites related to preeclampsia and gestational age. A list of CpG from a comparison of the differentially methylated CpG sites in preeclampsia with the results of Cruickshank et al.  on gestational age.
Pyrosequencing primers. A list of the primer sequences used to validate important CpG sites.