DNA methylation of cord blood cell types: Applications for mixed cell birth studies

Epigenome-wide association studies of disease widely use DNA methylation measured in blood as a surrogate tissue. Cell proportions can vary between people and confound associations of exposure or outcome. An adequate reference panel for estimating cell proportions from adult whole blood for DNA methylation studies is available, but an analogous cord blood cell reference panel is not yet available. Cord blood has unique cell types and the epigenetic signatures of standard cell types may not be consistent throughout the life course. Using magnetic bead sorting, we isolated cord blood cell types (nucleated red blood cells, granulocytes, monocytes, natural killer cells, B cells, CD4⁺T cells, and CD8⁺T cells) from 17 live births at Johns Hopkins Hospital. We confirmed enrichment of the cell types using fluorescence assisted cell sorting and ran DNA from the separated cell types on the Illumina Infinium HumanMethylation450 BeadChip array. After filtering, the final analysis was on 104 samples at 429,794 probes. We compared cell type specific signatures in cord to each other and methylation at 49.2% of CpG sites on the array differed by cell type (F-test P < 10^?8). Differences between nucleated red blood cells and the remainder of the cell types were most pronounced (36.9% of CpG sites at P < 10^?8) and 99.5% of these sites were hypomethylated relative to the other cell types. We also compared the mean-centered sorted cord profiles to the available adult reference panel and observed high correlation between the overlapping cell types for granulocytes and monocytes (both r=0.74), and poor correlation for CD8⁺T cells and NK cells (both r=0.08). We further provide an algorithm for estimating cell proportions in cord blood using the newly developed cord reference panel, which estimates biologically plausible cell proportions in whole cord blood samples.     

Cord blood buffy coat DNA methylation is comparable to whole cord blood methylation

Cord blood DNA methylation is associated with numerous health outcomes and environmental exposures. Whole cord blood DNA reflects all nucleated blood cell types, while centrifuging whole blood separates red blood cells, generating a white blood cell buffy coat. Both sample types are used in DNA methylation studies. Cell types have unique methylation patterns and processing can impact cell distributions, which may influence comparability. We evaluated differences in cell composition and DNA methylation between cord blood buffy coat and whole cord blood samples. Cord blood DNA methylation was measured with the Infinium EPIC BeadChip (Illumina) in eight individuals, each contributing buffy coat and whole blood samples. We analyzed principal components (PC) of methylation, performed hierarchical clustering, and computed correlations of mean-centered methylation between pairs. We conducted moderated t-tests on single sites and estimated cell composition. DNA methylation PCs were associated with individual (P_PC1 = 1.4 × 10^?9; P_PC2 = 2.9 × 10^?5; P_PC3 = 3.8 × 10^-5; P_PC4 = 4.2 × 10^-6; P_PC5 = 9.9 × 10^-13, P_PC6 = 1.3 × 10^?11) and not with sample type (P_PC1-6>0.7). Samples hierarchically clustered by individual. Pearson correlations of mean-centered methylation between paired samples ranged from r = 0.66 to r = 0.87. No individual site significantly differed between buffy coat and whole cord blood when adjusting for multiple comparisons (five sites had unadjusted P<10^?5). Estimated cell type proportions did not differ by sample type (P = 0.46), and estimated proportions were highly correlated between paired samples (r = 0.99). Differences in methylation and cell composition between buffy coat and whole cord blood are much lower than inter-individual variation, demonstrating that both sample preparation types can be analytically combined and compared.     

Deciphering the heterogeneity in DNA methylation patterns during stem cell differentiation and reprogramming

Background

Human induced pluripotent stem cells (iPSCs) have a wide range of applications throughout the fields of basic research, disease modeling and drug screening. Epigenetic instable iPSCs with aberrant DNA methylation may divide and differentiate into cancer cells. Unfortunately, little effort has been taken to compare the epigenetic variation in iPSCs with that in differentiated cells. Here, we developed an analytical procedure to decipher the DNA methylation heterogeneity of mixed cells and further exploited it to quantitatively assess the DNA methylation variation in the methylomes of adipose-derived stem cells (ADS), mature adipocytes differentiated from ADS cells (ADS-adipose) and iPSCs reprogrammed from ADS cells (ADS-iPSCs).

Results

We observed that the degree of DNA methylation variation varies across distinct genomic regions with promoter and 5’UTR regions exhibiting low methylation variation and Satellite showing high methylation variation. Compared with differentiated cells, ADS-iPSCs possess globally decreased methylation variation, in particular in repetitive elements. Interestingly, DNA methylation variation decreases in promoter regions during differentiation but increases during reprogramming. Methylation variation in promoter regions is negatively correlated with gene expression. In addition, genes showing a bipolar methylation pattern, with both completely methylated and completely unmethylated reads, are related to the carbohydrate metabolic process, cellular development, cellular growth, proliferation, etc.

Conclusions

This study delivers a way to detect cell-subset specific methylation genes in a mixed cell population and provides a better understanding of methylation dynamics during stem cell differentiation and reprogramming.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-978) contains supplementary material, which is available to authorized users.     

Validation of a DNA methylation reference panel for the estimation of nucleated cells types in cord blood

Cord blood is widely used as surrogate tissue in epigenome-wide association studies of prenatal conditions. Cell type composition variation across samples can be an important confounder of epigenome-wide association studies in blood that constitute a mixture of cells. We evaluated a newly developed cord blood reference panel to impute cell type composition from DNA methylation levels, including nucleated red blood cells (nRBCs). We estimated cell type composition from 154 unique cord blood samples with available DNA methylation data as well as direct measurements of nucleated cell types. We observed high correlations between the estimated and measured composition for nRBCs (r = 0.92, R² = 0.85), lymphocytes (r = 0.77, R² = 0.58), and granulocytes (r = 0.72, R² = 0.52), and a moderate correlation for monocytes (r = 0.51, R² = 0.25) as well as relatively low root mean square errors from the residuals ranging from 1.4 to 5.4%. These results validate the use of the cord blood reference panel and highlight its utility and limitations for epidemiological studies.     

Effect of prenatal arsenic exposure on DNA methylation and leukocyte subpopulations in cord blood

Prenatal arsenic exposure is associated with increased risk of disease in adulthood. This has led to considerable interest in arsenic’s ability to disrupt fetal programming. Many studies report that arsenic exposure alters DNA methylation in whole blood but these studies did not adjust for cell mixture. In this study, we examined the relationship between arsenic in maternal drinking water collected ≤ 16 weeks gestational age and DNA methylation in cord blood (n = 44) adjusting for leukocyte-tagged differentially methylated regions. DNA methylation was quantified using the Infinium HumanMethylation 450 BeadChip array. Recursively partitioned mixture modeling examined the relationship between arsenic and methylation at 473,844 CpG sites. Median arsenic concentration in water was 12 µg/L (range < 1- 510 µg/L). Log₁₀ arsenic was associated with altered DNA methylation across the epigenome (P = 0.002); however, adjusting for leukocyte distributions attenuated this association (P = 0.013). We also observed that arsenic had a strong effect on the distribution of leukocytes in cord blood. In adjusted models, every log₁₀ increase in maternal drinking water arsenic exposure was estimated to increase CD8+ T cells by 7.4% (P = 0.0004) and decrease in CD4+ T cells by 9.2% (P = 0.0002). These results show that prenatal exposure to arsenic had an exposure-dependent effect on specific T cell subpopulations in cord blood and altered DNA methylation in cord blood. Future research is needed to determine if these small changes in DNA methylation alter gene expression or are associated with adverse health effects.     

Effect of prenatal arsenic exposure on DNA methylation and leukocyte subpopulations in cord blood

Prenatal arsenic exposure is associated with increased risk of disease in adulthood. This has led to considerable interest in arsenic’s ability to disrupt fetal programming. Many studies report that arsenic exposure alters DNA methylation in whole blood but these studies did not adjust for cell mixture. In this study, we examined the relationship between arsenic in maternal drinking water collected ≤ 16 weeks gestational age and DNA methylation in cord blood (n = 44) adjusting for leukocyte-tagged differentially methylated regions. DNA methylation was quantified using the Infinium HumanMethylation 450 BeadChip array. Recursively partitioned mixture modeling examined the relationship between arsenic and methylation at 473,844 CpG sites. Median arsenic concentration in water was 12 µg/L (range < 1- 510 µg/L). Log₁₀ arsenic was associated with altered DNA methylation across the epigenome (P = 0.002); however, adjusting for leukocyte distributions attenuated this association (P = 0.013). We also observed that arsenic had a strong effect on the distribution of leukocytes in cord blood. In adjusted models, every log₁₀ increase in maternal drinking water arsenic exposure was estimated to increase CD8+ T cells by 7.4% (P = 0.0004) and decrease in CD4+ T cells by 9.2% (P = 0.0002). These results show that prenatal exposure to arsenic had an exposure-dependent effect on specific T cell subpopulations in cord blood and altered DNA methylation in cord blood. Future research is needed to determine if these small changes in DNA methylation alter gene expression or are associated with adverse health effects.     

Comparison of DNA methylation measured by Illumina 450K and EPIC BeadChips in blood of newborns and 14-year-old children

Analysis of DNA methylation helps to understand the effects of environmental exposures as well as the role of epigenetics in human health. Illumina, Inc. recently replaced the HumanMethylation450 BeadChip (450K) with the EPIC BeadChip, which nearly doubles the measured CpG sites to >850,000. Although the new chip uses the same underlying technology, it is important to establish if data between the two platforms are comparable within cohorts and for meta-analyses. DNA methylation was assessed by 450K and EPIC using whole blood from newborn (n = 109) and 14-year-old (n = 86) participants of the Center for the Health Assessment of Mothers and Children of Salinas. The overall per-sample correlations were very high (r >0.99), although many individual CpG sites, especially those with low variance of methylation, had lower correlations (median r = 0.24). There was also a small subset of CpGs with large mean methylation β-value differences between platforms, in both the newborn and 14-year datasets. However, estimates of cell type proportion prediction by 450K and EPIC were highly correlated at both ages. Finally, differentially methylated positions between boys and girls replicated very well by both platforms in newborns and older children. These findings are encouraging for application of combined data from EPIC and 450K platforms for birth cohorts and other population studies. These data in children corroborate recent comparisons of the two BeadChips in adults and in cancer cell lines. However, researchers should be cautious when characterizing individual CpG sites and consider independent methods for validation of significant hits.     

Novel region discovery method for Infinium 450K DNA methylation data reveals changes associated with aging in muscle and neuronal pathways

We describe a methodology for detecting differentially methylated regions (DMRs) and variably methylated regions (VMRs), in data from Infinium 450K arrays that are very widely used in epigenetic studies. Region detection is more specific than single CpG analysis as it increases the extent of common findings between studies, and is more powerful as it reduces the multiple testing problem inherent in epigenetic whole‐genome association studies (EWAS). In addition, results driven by single erroneous probes are removed. We have used multiple publicly available Infinium 450K data sets to generate a consensus list of DMRs for age, supporting the hypothesis that aging is associated with specific epigenetic modifications. The consensus aging DMRs are significantly enriched for muscle biogenesis pathways. We find a massive increase in VMRs with age and in regions of the genome associated with open chromatin and neurotransmission. Old age VMRs are significantly enriched for neurotransmission pathways. EWAS studies should investigate the role of this interindividual variation in DNA methylation, in the age‐associated diseases of sarcopenia and dementia.     

Genome-wide DNA methylation associations with spontaneous preterm birth in US blacks: findings in maternal and cord blood samples

Preterm birth (PTB) affects one in six Black babies in the United States. Epigenetics is believed to play a role in PTB; however, only a limited number of epigenetic studies of PTB have been reported, most of which have focused on cord blood DNA methylation (DNAm) and/or were conducted in white populations. Here we conducted, by far, the largest epigenome-wide DNAm analysis in 300 Black women who delivered early spontaneous preterm (sPTB, n = 150) or full-term babies (n = 150) and replicated the findings in an independent set of Black mother-newborn pairs from the Boston Birth Cohort. DNAm in maternal blood and/or cord blood was measured using the Illumina HumanMethylation450 BeadChip. We identified 45 DNAm loci in maternal blood associated with early sPTB, with a false discovery rate (FDR) <5%. Replication analyses confirmed sPTB associations for cg03915055 and cg06804705, located in the promoter regions of the CYTIP and LINC00114 genes, respectively. Both loci had comparable associations with early sPTB and early medically-indicated PTB, but attenuated associations with late sPTB. These associations could not be explained by cell composition, gestational complications, and/or nearby maternal genetic variants. Analyses in the newborns of the 110 Black women showed that cord blood methylation levels at both loci had no associations with PTB. The findings from this study underscore the role of maternal DNAm in PTB risk, and provide a set of maternal loci that may serve as biomarkers for PTB. Longitudinal studies are needed to clarify temporal relationships between maternal DNAm and PTB risk.     

Age-related epigenome-wide DNA methylation and hydroxymethylation in longitudinal mouse blood

DNA methylation at cytosine-phosphate-guanine (CpG) dinucleotides changes as a function of age in humans and animal models, a process that may contribute to chronic disease development. Recent studies have investigated the role of an oxidized form of DNA methylation – 5-hydroxymethylcytosine (5hmC) – in the epigenome, but its contribution to age-related DNA methylation remains unclear. We tested the hypothesis that 5hmC changes with age, but in a direction opposite to 5-methylcytosine (5mC), potentially playing a distinct role in aging. To characterize epigenetic aging, genome-wide 5mC and 5hmC were measured in longitudinal blood samples (2, 4, and 10 months of age) from isogenic mice using two sequencing methods – enhanced reduced representation bisulfite sequencing and hydroxymethylated DNA immunoprecipitation sequencing. Examining the epigenome by age, we identified 28,196 unique differentially methylated CpGs (DMCs) and 8,613 differentially hydroxymethylated regions (DHMRs). Mouse blood showed a general pattern of epigenome-wide hypermethylation and hypo-hydroxymethylation with age. Comparing age-related DMCs and DHMRs, 1,854 annotated genes showed both differential 5mC and 5hmC, including one gene – Nfic – at five CpGs in the same 250 bp chromosomal region. At this region, 5mC and 5hmC levels both decreased with age. Reflecting these age-related epigenetic changes, Nfic RNA expression in blood decreased with age, suggesting that age-related regulation of this gene may be driven by 5hmC, not canonical DNA methylation. Combined, our genome-wide results show age-related differential 5mC and 5hmC, as well as some evidence that changes in 5hmC may drive age-related DNA methylation and gene expression.     

Gestational intake of methyl donors and global LINE-1 DNA methylation in maternal and cord blood: Prospective results from a folate-replete population

Maternal diet affects offspring DNA methylation in animal models, but evidence from humans is limited. We investigated the extent to which gestational intake of methyl donor nutrients affects global DNA methylation in maternal and umbilical cord blood. Among mother-infant pairs in Project Viva, a folate-replete US population, we estimated maternal intakes of vitamin B12, betaine, choline, folate, cadmium, zinc and iron periconceptionally and during the second trimester. We examined associations of these nutrients with DNA methylation, measured as %5-methyl cytosines (%5mC) in Long Interspersed Nuclear Element-1 (LINE-1), in first trimester (n = 830) and second trimester (n = 671) maternal blood and in cord blood at delivery (n = 516). Cord blood methylation was higher for male than female infants {mean [standard deviation (SD)] 84.8 [0.6] vs. 84.4 [0.7]%}. In the multivariable-adjusted model, maternal intake of methyl donor nutrients periconceptionally and during the second trimester of pregnancy was not positively associated with first trimester, second trimester or cord blood LINE-1 methylation. Periconceptional betaine intake was inversely associated with cord blood methylation [regression coefficient = -0.08% (95% confidence interval (CI): -0.14,-0.01)] but this association was attenuated after adjustment for dietary cadmium, which itself was directly associated with first trimester methylation and inversely associated with cord blood methylation. We also found an inverse association between periconceptional choline [-0.10%, 95% CI: -0.17,-0.03 for each SD (~63 mg/d)] and cord blood methylation in males only. In this folate-replete population, we did not find positive associations between intake of methyl donor nutrients during pregnancy and DNA methylation overall, but among males, higher early pregnancy intakes of choline were associated with lower cord blood methylation.     

Differential DNA methylation identified in the blood and retina of AMD patients

Age-related macular degeneration (AMD) is a major cause of blindness in the western world. While genetic studies have linked both common and rare variants in genes involved in regulation of the complement system to increased risk of development of AMD, environmental factors, such as smoking and nutrition, can also significantly affect the risk of developing the disease and the rate of disease progression. Since epigenetics has been implicated in mediating, in part, the disease risk associated with some environmental factors, we investigated a possible epigenetic contribution to AMD. We performed genome-wide DNA methylation profiling of blood from AMD patients and controls. No differential methylation site reached genome-wide significance; however, when epigenetic changes in and around known GWAS-defined AMD risk loci were explored, we found small but significant DNA methylation differences in the blood of neovascular AMD patients near age-related maculopathy susceptibility 2 (ARMS2), a top-ranked GWAS locus preferentially associated with neovascular AMD. The methylation level of one of the CpG sites significantly correlated with the genotype of the risk SNP rs10490924, suggesting a possible epigenetic mechanism of risk. Integrating genome-wide DNA methylation analysis of retina samples with and without AMD together with blood samples, we further identified a consistent, replicable change in DNA methylation in the promoter region of protease serine 50 (PRSS50). These methylation changes may identify sites in novel genes that are susceptible to non-genetic factors known to contribute to AMD development and progression.     

An epigenome-wide association meta-analysis of prenatal maternal stress in neonates: A model approach for replication

Prenatal maternal stress exposure has been associated with neonatal differential DNA methylation. However, the available evidence in humans is largely based on candidate gene methylation studies, where only a few CpG sites were evaluated. The aim of this study was to examine the association between prenatal exposure to maternal stress and offspring genome-wide cord blood methylation using different methods. First, we conducted a meta-analysis and follow-up pathway analyses. Second, we used novel region discovery methods [i.e., differentially methylated regions (DMRs) analyses]. To this end, we used data from two independent population-based studies, the Generation R Study (n = 912) and the Avon Longitudinal Study of Parents and Children (ALSPAC, n = 828), to (i) measure genome-wide DNA methylation in cord blood and (ii) extract a prenatal maternal stress composite. The meta-analysis (n_total = 1,740) revealed no epigenome-wide (meta P <1.00e-07) associations of prenatal maternal stress exposure with neonatal differential DNA methylation. Follow-up analyses of the top hits derived from our epigenome-wide meta-analysis (meta P <1.00e-04) indicated an over-representation of the methyltransferase activity pathway. We identified no Bonferroni-corrected (P <1.00e-06) DMRs associated with prenatal maternal stress exposure. Combining data from two independent population-based samples in an epigenome-wide meta-analysis, the current study indicates that there are no large effects of prenatal maternal stress exposure on neonatal DNA methylation. Such replication efforts are essential in the search for robust associations, whether derived from candidate gene methylation or epigenome-wide studies.     

Blood glutathione redox status and global methylation of peripheral blood mononuclear cell DNA in Bangladeshi adults

Oxidative stress and DNA methylation are metabolically linked through the relationship between one-carbon metabolism and the transsulfuration pathway, but possible modulating effects of oxidative stress on DNA methylation have not been extensively studied in humans. Enzymes involved in DNA methylation, including DNA methyltransferases and histone deacetylases, may show altered activity under oxidized cellular conditions. Additionally, in vitro studies suggest that glutathione (GSH) depletion leads to global DNA hypomethylation, possibly through the depletion of S-adenosylmethionine (SAM). We tested the hypothesis that a more oxidized blood GSH redox status is associated with decreased global peripheral blood mononuclear cell (PBMC) DNA methylation in a sample of Bangladeshi adults. Global PBMC DNA methylation and whole blood GSH, glutathione disulfide (GSSG), and SAM concentrations were measured in 320 adults. DNA methylation was measured by using the [³H]-methyl incorporation assay; values are inversely related to global DNA methylation. Whole blood GSH redox status (E_h) was calculated using the Nernst equation. We found that a more oxidized blood GSH E_h was associated with decreased global DNA methylation (B ± SE, 271 ± 103, p = 0.009). Blood SAM and blood GSH were associated with global DNA methylation, but these relationships did not achieve statistical significance. Our findings support the hypothesis that a more oxidized blood GSH redox status is associated with decreased global methylation of PBMC DNA. Furthermore, blood SAM does not appear to mediate this association. Future research should explore mechanisms through which cellular redox might influence global DNA methylation.     

Type 2 diabetes mellitus in relation to global LINE-1 DNA methylation in peripheral blood: A cohort study

In the last years, epigenetic processes have emerged as a promising area of complex diseases research. DNA methylation measured in Long Interspersed Nucleotide Element 1 (LINE-1) sequences has been considered a surrogate marker for global genome methylation. New findings have suggested the potential involvement of epigenetic mechanisms in Type 2 diabetes (T2DM) as a crucial interface between the effects of genetic predisposition and environmental influences. Our study evaluated whether global DNA methylation predicted increased risk from T2DM or other carbohydrate metabolism disorders in a cohort study. We used a prospective cohort intervention study and a control group. We collected phenotypic, anthropometric, biochemical, and nutritional information from all subjects. Global LINE-1 DNA methylation was quantified by pyrosequencing technology. Subjects that did not improve their carbohydrate metabolism status showed lower levels of global LINE-1 DNA methylation (63.9 ± 1.7 vs. 64.7 ± 2.4) and they practiced less intense physical activity (5.8% vs. 21.5%). Logistic regression analyses showed a significant association between LINE-1 DNA methylation and metabolic status after adjustment for sex, age, BMI, and physical activity. Our study showed that lower LINE-1 DNA methylation levels were associated with a higher risk metabolic status worsening, independent of other classic risk factors. This finding highlights the potential role for epigenetic biomarkers as predictors of T2DM risk or other related metabolic disorders.     

Type 2 diabetes mellitus in relation to global LINE-1 DNA methylation in peripheral blood: A cohort study

In the last years, epigenetic processes have emerged as a promising area of complex diseases research. DNA methylation measured in Long Interspersed Nucleotide Element 1 (LINE-1) sequences has been considered a surrogate marker for global genome methylation. New findings have suggested the potential involvement of epigenetic mechanisms in Type 2 diabetes (T2DM) as a crucial interface between the effects of genetic predisposition and environmental influences. Our study evaluated whether global DNA methylation predicted increased risk from T2DM or other carbohydrate metabolism disorders in a cohort study. We used a prospective cohort intervention study and a control group. We collected phenotypic, anthropometric, biochemical, and nutritional information from all subjects. Global LINE-1 DNA methylation was quantified by pyrosequencing technology. Subjects that did not improve their carbohydrate metabolism status showed lower levels of global LINE-1 DNA methylation (63.9 ± 1.7 vs. 64.7 ± 2.4) and they practiced less intense physical activity (5.8% vs. 21.5%). Logistic regression analyses showed a significant association between LINE-1 DNA methylation and metabolic status after adjustment for sex, age, BMI, and physical activity. Our study showed that lower LINE-1 DNA methylation levels were associated with a higher risk metabolic status worsening, independent of other classic risk factors. This finding highlights the potential role for epigenetic biomarkers as predictors of T2DM risk or other related metabolic disorders.