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.     

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.     

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 role of immune cell composition in epigenetic age acceleration: Insights from cell-type deconvolution applied to human blood epigenetic clocks

Aging is a significant risk factor for various human disorders, and DNA methylation clocks have emerged as powerful tools for estimating biological age and predicting health-related outcomes. Methylation data from blood DNA has been a focus of more recently developed DNA methylation clocks. However, the impact of immune cell composition on epigenetic age acceleration (EAA) remains unclear as only some clocks incorporate partial cell type composition information when analyzing EAA. We investigated associations of 12 immune cell types measured by cell-type deconvolution with EAA predicted by six widely-used DNA methylation clocks in data from >10,000 blood samples. We observed significant associations of immune cell composition with EAA for all six clocks tested. Across the clocks, nine or more of the 12 cell types tested exhibited significant associations with EAA. Higher memory lymphocyte subtype proportions were associated with increased EAA, and naïve lymphocyte subtypes were associated with decreased EAA. To demonstrate the potential confounding of EAA by immune cell composition, we applied EAA in rheumatoid arthritis. Our research maps immune cell type contributions to EAA in human blood and offers opportunities to adjust for immune cell composition in EAA studies to a significantly more granular level. Understanding associations of EAA with immune profiles has implications for the interpretation of epigenetic age and its relevance in aging and disease research. Our detailed map of immune cell type contributions serves as a resource for studies utilizing epigenetic clocks across diverse research fields, including aging-related diseases, precision medicine, and therapeutic interventions.     

Epigenome-wide profiling identifies significant differences in DNA methylation between matched-pairs of T- and B-lymphocytes from healthy individuals

Multiple reports now describe changes to the DNA methylome in rheumatoid arthritis and in many cases have analyzed methylation in mixed cell populations from whole blood. However, these approaches may preclude the identification of cell type-specific methylation, which may subsequently bias identification of disease-specific changes. To address this possibility, we conducted genome-wide DNA methylation profiling using HumanMethylation450 BeadChips to identify differences within matched pairs of T-lymphocytes and B-lymphocytes isolated from the peripheral blood of 10 healthy females. Array data were processed and differential methylation identified using NIMBL software. Validation of array data was performed by bisulfite pyrosequencing. Genome-wide DNA methylation was initially determined by analysis of LINE-1 sequences and was higher in B-lymphocytes than matched T-lymphocytes (69.8% vs. 65.2%, P ≤ 0.01). Pairwise analysis identified 679 CpGs, representing 250 genes, which were differentially methylated between T-lymphocytes and B-lymphocytes. The majority of sites (76.6%) were hypermethylated in B-lymphocytes. Pyrosequencing of selected candidates confirmed the array data in all cases. Hierarchical clustering revealed perfect segregation of samples into two distinct clusters based on cell type. Differentially methylated genes showed enrichment for biological functions/pathways associated with leukocytes and T-lymphocytes. Our work for the first time shows that T-lymphocytes and B-lymphocytes possess intrinsic differences in DNA methylation within a restricted set of functionally related genes. These data provide a foundation for investigating DNA methylation in diseases in which these cell types play important and distinct roles.     

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.     

Blood-based profiles of DNA methylation predict the underlying distribution of cell types

The potential influence of underlying differences in relative leukocyte distributions in studies involving blood-based profiling of DNA methylation is well recognized and has prompted development of a set of statistical methods for inferring changes in the distribution of white blood cells using DNA methylation signatures. However, the extent to which this methodology can accurately predict cell-type proportions based on blood-derived DNA methylation data in a large-scale epigenome-wide association study (EWAS) has yet to be examined. We used publicly available data deposited in the Gene Expression Omnibus (GEO) database (accession number GSE37008), which consisted of both blood-derived epigenome-wide DNA methylation data assayed using the Illumina Infinium HumanMethylation27 BeadArray and complete blood cell (CBC) counts among a community cohort of 94 non-diseased individuals. Constrained projection (CP) was used to obtain predictions of the proportions of lymphocytes, monocytes and granulocytes for each of the study samples based on their DNA methylation signatures. Our findings demonstrated high consistency between the average CBC-derived and predicted percentage of monocytes and lymphocytes (17.9% and 17.6% for monocytes and 82.1% and 81.4% for lymphocytes), with root mean squared error (rMSE) of 5% and 6%, for monocytes and lymphocytes, respectively. Similarly, there was moderate-high correlation between the CP-predicted and CBC-derived percentages of monocytes and lymphocytes (0.60 and 0.61, respectively), and these results were robust to the number of leukocyte differentially methylated regions (L-DMRs) used for CP prediction. These results serve as further validation of the CP approach and highlight the promise of this technique for EWAS where DNA methylation is profiled using whole-blood genomic DNA.     

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.     

Early life socioeconomic factors and genomic DNA methylation in mid-life

Epigenetic modifications may be one mechanism linking early life factors, including parental socioeconomic status (SES), to adult onset disease risk. However, SES influences on DNA methylation patterns remain largely unknown. In a US birth cohort of women, we examined whether indicators of early life and adult SES were associated with white blood cell methylation of repetitive elements (Sat2, Alu and LINE-1) in adulthood. Low family income at birth was associated with higher Sat2 methylation (β = 19.7, 95% CI: 0.4, 39.0 for lowest vs. highest income quartile) and single parent family was associated with higher Alu methylation (β = 23.5, 95% CI: 2.6, 44.4), after adjusting for other early life factors. Lower adult education was associated with lower Sat2 methylation (β = -16.7, 95% CI: -29.0, -4.5). There were no associations between early life SES and LINE-1 methylation. Overall, our preliminary results suggest possible influences of SES across the life-course on genomic DNA methylation in adult women. However, these preliminary associations need to be replicated in larger prospective studies.     

Prenatal exposure to polycyclic aromatic hydrocarbons could increase the risk of low birth weight by affecting the DNA methylation states in a Chinese cohort

Polycyclic aromatic hydrocarbons (PAHs), as a kind of endocrine disruptors, can enter the fetus body cross the placental barrier from prenatal PAHs exposure to cause adverse birth outcomes. However, it is controversial association between prenatal PAHs exposure and low birth weight (LBW) of their infants. So the present study aimed to estimate the effects of prenatal PAHs exposure during the pregnancy on the risk of LBW in a Chinese cohort through modifying the DNA methylation states. A longitudinal prospective study with 407 pregnant women was established from May to October 2019. The prenatal PAHs exposure during the pregnancy was assessed using the internal dose such as the PAHs metabolites and PAH-DNA adducts in the umbilical cord blood. The methylation levels of genomic DNA and growth-related genes (IGF1 and IGF2) were assessed, while the expressions of these genes were both determined by RT-PCR and Elisa methods. The growth outcomes and relevant Z-scores were recorded at birth. The correlations between the DNA methylation status and concentrations of PAHs, expression levels of growth-related genes and body weight/WAZ were investigated as the measures. According to the PAH-DNA adducts, the subjects were divided into two groups: PAHs-exposed group (PAH-DNA adducts>0, n = 55) and non-exposed group (PAH-DNA adducts = 0, n = 352). Compared with the non-exposed group, it displayed marked decreased birth weight, and increased concentrations of PAHs and DNA methylation levels of the global genomic, IGF1 and IGF2 with their lower expressions in the PAHs-exposed group. These hypermethylation (global genomic, CpG14 and CpG15 of IGF1, and CpG14 of IGF2) were positively correlated with the contents of PAHs in the umbilical cord blood, and negatively correlated with the growth outcomes and their expressions. Totally, prenatal PAHs exposures may contribute to an increased risk of LBW of their infants by modulating the DNA methylation states of genomic DNA and growth-related genes (IGF1 and IGF2) in the umbilical cord blood, which could provide the prenatal prevention of PAHs exposure from possible environmental media except from the occupation and tobacco usage to ensure the health of their infants.     

DNA methylation arrays as surrogate measures of cell mixture distribution

ABSTRACT: BACKGROUND: There has been a long-standing need in biomedical research for a method that quantifies the normally mixed composition of leukocytes beyond what is possible by simple histological or flow cytometric assessments. The latter is restricted by the labile nature of protein epitopes, requirements for cell processing, and timely cell analysis. In a diverse array of diseases and following numerous immune-toxic exposures, leukocyte composition will critically inform the underlying immuno-biology to most chronic medical conditions. Emerging research demonstrates that DNA methylation is responsible for cellular differentiation, and when measured in whole peripheral blood, serves to distinguish cancer cases from controls. RESULTS: Here we present a method, similar to regression calibration, for inferring changes in the distribution of white blood cells between different subpopulations (e.g. cases and controls) using DNA methylation signatures, in combination with a previously obtained external validation set consisting of signatures from purified leukocyte samples. We validate the fundamental idea in a cell mixture reconstruction experiment, then demonstrate our method on DNA methylation data sets from several studies, including data from a Head and Neck Squamous Cell Carcinoma (HNSCC) study and an ovarian cancer study. Our method produces results consistent with prior biological findings, thereby validating the approach. CONCLUSIONS: Our method, in combination with an appropriate external validation set, promises new opportunities for large-scale immunological studies of both disease states and noxious exposures.     

LINE-1 DNA methylation is inversely correlated with cord plasma homocysteine in man: A preliminary study

Folic acid supplementation during pregnancy has known beneficial effects. It reduces risk of neural tube defects and low birth weight. Folate and other one-carbon intermediates might secure these clinical effects via DNA methylation. However, most data on the effects of folate on the epigenome is derived from animal or in vitro models. We examined the relationship between cord blood methylation and maternal folic acid intake, cord blood folate and homocysteine using data from 24 pregnant women. Genome-wide methylation was determined by the level of methylation of LINE-1 repeats using Pyrosequencing. We show that cord plasma homocysteine (p = 0.001, r = -0.688), but not serum folate or maternal folic acid intake, is inverse correlated with LINE-1 methylation. This remained significant after correction for potential confounders (p = 0.004). These data indicate that levels of folate-associated intermediates in cord blood during late pregnancy have significant consequences for the fetal epigenome.     

Quantitative,high-resolution epigenetic profiling of CpG loci identifies associations with cord blood plasma homocysteine and birth weight in humans

Supplementation with folic acid during pregnancy is known to reduce the risk of neural tube defects and low birth weight. It is thought that folate and other one-carbon intermediates might secure these clinical effects via DNA methylation. We examined the effects of folate on the human methylome using quantitative interrogation of 27,578 CpG loci associated with 14,496 genes at single-nucleotide resolution across 12 fetal cord blood samples. Consistent with previous studies, the majority of CpG dinucleotides located within CpG islands exhibited hypomethylation while those outside CpG islands showed mid-high methylation. However, for the first time in human samples, unbiased analysis of methylation across samples revealed a significant correlation of methylation patterns with plasma homocysteine, LINE-1 methylation and birth weight centile. Additionally, CpG methylation significantly correlated with either birth weight or LINE-1 methylation were predominantly located in CpG islands. These data indicate that levels of folate-associated intermediates in cord blood reflect their influence and consequences for the fetal epigenome and potentially on pregnancy outcome. In these cases, their influence might be exerted during late gestation or reflect those present during the peri-conceptual period.Key words: cord blood, birth weight, folic acid, homocysteine, BeadArray, hierarchical clustering, Illumina     

Quantitative,high-resolution epigenetic profiling of CpG loci identifies associations with cord blood plasma homocysteine and birth weight in humans

Supplementation with folic acid during pregnancy is known to reduce the risk of neural tube defects and low birth weight. It is thought that folate and other one-carbon intermediates might secure these clinical effects via DNA methylation. We examined the effects of folate on the human methylome using quantitative interrogation of 27,578 CpG loci associated with 14,496 genes at single-nucleotide resolution across 12 fetal cord blood samples. Consistent with previous studies, the majority of CpG dinucleotides located within CpG islands exhibited hypo-methylation while those outside CpG islands showed mid-high methylation. However, for the first time in human samples, unbiased analysis of methylation across samples revealed a significant correlation of methylation patterns with plasma homocysteine, LINE-1 methylation and birth weight centile. Additionally, CpG methylation significantly correlated with either birth weight or LINE-1 methylation were predominantly located in CpG islands. These data indicate that levels of folate-associated intermediates in cord blood reflect their influence and consequences for the fetal epigenome and potentially on pregnancy outcome. In these cases, their influence might be exerted during late gestation or reflect those present during the peri-conceptual period.