Characterizing the hypomethylated DNA methylation profile of nucleated red blood cells from cord blood

Placental and Immune Cell DNA Methylation Reference Panel for Bulk Tissue Cell Composition Estimation in Epidemiological Studies

To distinguish DNA methylation (DNAm) from cell proportion changes in whole placental tissue research, we developed a robust cell type-specific DNAm reference to estimate cell composition. We collated newly collected and existing cell type DNAm profiles quantified via Illumina EPIC or 450k microarrays. To estimate cell composition, we deconvoluted whole placental samples (n=36) with robust partial correlation based on the top 50 hyper- and hypomethylated sites per cell type. To test deconvolution performance, we evaluated RMSE in predicting principal component one of DNAm variation in 204 external placental samples. We analyzed DNAm profiles (n=368,435 sites) from 12 cell types: cytotrophoblasts (n=18), endothelial cells (n=19), Hofbauer cells (n=26), stromal cells (n=21), syncytiotrophoblasts (n=4), six lymphocyte types (n=36), and nucleated red blood cells (n=11). Median cell composition was consistent with placental biology: 60.4% syncytiotrophoblast, 17.1% stromal, 8.8% endothelial, 4.5% cytotrophoblast, 3.9% Hofbauer, 1.7% nucleated red blood cells, and 1.2% neutrophils. Our expanded reference outperformed an existing reference in predicting DNAm variation (15.4% variance explained, IQR=21.61) with cell composition estimates (RMSE:10.51 vs. 11.43, p-value<0.001). This cell type reference can robustly estimate cell composition from whole placental DNAm data to detect important cell types, reveal biological mechanisms, and improve casual inference.

Epigenomic signature of major congenital heart defects in newborns with Down syndrome

BACKGROUND

Congenital heart defects (CHDs) affect approximately half of individuals with Down syndrome (DS), but the molecular reasons for incomplete penetrance are unknown. Previous studies have largely focused on identifying genetic risk factors associated with CHDs in individuals with DS, but comprehensive studies of the contribution of epigenetic marks are lacking. We aimed to identify and characterize DNA methylation differences from newborn dried blood spots (NDBS) of DS individuals with major CHDs compared to DS individuals without CHDs.

METHODS

We used the Illumina EPIC array and whole-genome bisulfite sequencing (WGBS) to quantitate DNA methylation for 86 NDBS samples from the California Biobank Program: (1) 45 DS-CHD (27 female, 18 male) and (2) 41 DS non-CHD (27 female, 14 male). We analyzed global CpG methylation and identified differentially methylated regions (DMRs) in DS-CHD versus DS non-CHD comparisons (both sex-combined and sex-stratified) corrected for sex, age of blood collection, and cell-type proportions. CHD DMRs were analyzed for enrichment in CpG and genic contexts, chromatin states, and histone modifications by genomic coordinates and for gene ontology enrichment by gene mapping. DMRs were also tested in a replication dataset and compared to methylation levels in DS versus typical development (TD) WGBS NDBS samples.

RESULTS

We found global CpG hypomethylation in DS-CHD males compared to DS non-CHD males, which was attributable to elevated levels of nucleated red blood cells and not seen in females. At a regional level, we identified 58, 341, and 3938 CHD-associated DMRs in the Sex Combined, Females Only, and Males Only groups, respectively, and used machine learning algorithms to select 19 Males Only loci that could distinguish CHD from non-CHD. DMRs in all comparisons were enriched for gene exons, CpG islands, and bivalent chromatin and mapped to genes enriched for terms related to cardiac and immune functions. Lastly, a greater percentage of CHD-associated DMRs than background regions were differentially methylated in DS versus TD samples.

CONCLUSIONS

A sex-specific signature of DNA methylation was detected in NDBS of DS-CHD compared to DS non-CHD individuals. This supports the hypothesis that epigenetics can reflect the variability of phenotypes in DS, particularly CHDs.

Prenatal exposures to per- and polyfluoroalkyl substances and epigenetic aging in umbilical cord blood: The Healthy Start study

BACKGROUND

Per- and polyfluoroalkyl substances (PFAS) are ubiquitous, environmentally persistent chemicals, and prenatal exposures have been associated with adverse child health outcomes. Prenatal PFAS exposure may lead to epigenetic age acceleration (EAA), defined as the discrepancy between an individual's chronologic and epigenetic or biological age.

OBJECTIVES

We estimated associations of maternal serum PFAS concentrations with EAA in umbilical cord blood DNA methylation using linear regression, and a multivariable exposure-response function of the PFAS mixture using Bayesian kernel machine regression.

METHODS

Five PFAS were quantified in maternal serum (median: 27 weeks of gestation) among 577 mother-infant dyads from a prospective cohort. Cord blood DNA methylation data were assessed with the Illumina HumanMethylation450 array. EAA was calculated as the residuals from regressing gestational age on epigenetic age, calculated using a cord-blood specific epigenetic clock. Linear regression tested for associations between each maternal PFAS concentration with EAA. Bayesian kernel machine regression with hierarchical selection estimated an exposure-response function for the PFAS mixture.

RESULTS

In single pollutant models we observed an inverse relationship between perfluorodecanoate (PFDA) and EAA (-0.148 weeks per log-unit increase, 95% CI: -0.283, -0.013). Mixture analysis with hierarchical selection between perfluoroalkyl carboxylates and sulfonates indicated the carboxylates had the highest group posterior inclusion probability (PIP), or relative importance. Within this group, PFDA had the highest conditional PIP. Univariate predictor-response functions indicated PFDA and perfluorononanoate were inversely associated with EAA, while perfluorohexane sulfonate had a positive association with EAA.

CONCLUSIONS

Maternal mid-pregnancy serum concentrations of PFDA were negatively associated with EAA in cord blood, suggesting a pathway by which prenatal PFAS exposures may affect infant development. No significant associations were observed with other PFAS. Mixture models suggested opposite directions of association between perfluoroalkyl sulfonates and carboxylates. Future studies are needed to determine the importance of neonatal EAA for later child health outcomes.

Epigenomic signature of major congenital heart defects in newborns with Down syndrome

Background

Congenital heart defects (CHDs) affect approximately half of individuals with Down syndrome (DS) but the molecular reasons for incomplete penetrance are unknown. Previous studies have largely focused on identifying genetic risk factors associated with CHDs in individuals with DS, but comprehensive studies of the contribution of epigenetic marks are lacking. We aimed to identify and characterize DNA methylation differences from newborn dried blood spots (NDBS) of DS individuals with major CHDs compared to DS individuals without CHDs.

Methods

We used the Illumina EPIC array and whole-genome bisulfite sequencing (WGBS) to quantitate DNA methylation for 86 NDBS samples from the California Biobank Program: 1) 45 DS-CHD (27 female, 18 male) and 2) 41 DS non-CHD (27 female, 14 male). We analyzed global CpG methylation and identified differentially methylated regions (DMRs) in DS-CHD vs DS non-CHD comparisons (both sex-combined and sex-stratified) corrected for sex, age of blood collection, and cell type proportions. CHD DMRs were analyzed for enrichment in CpG and genic contexts, chromatin states, and histone modifications by genomic coordinates and for gene ontology enrichment by gene mapping. DMRs were also tested in a replication dataset and compared to methylation levels in DS vs typical development (TD) WGBS NDBS samples.

Results

We found global CpG hypomethylation in DS-CHD males compared to DS non-CHD males, which was attributable to elevated levels of nucleated red blood cells and not seen in females. At a regional level, we identified 58, 341, and 3,938 CHD-associated DMRs in the Sex Combined, Females Only, and Males Only groups, respectively, and used machine learning algorithms to select 19 Males Only loci that could distinguish CHD from non-CHD. DMRs in all comparisons were enriched for gene exons, CpG islands, and bivalent chromatin and mapped to genes enriched for terms related to cardiac and immune functions. Lastly, a greater percentage of CHD-associated DMRs than background regions were differentially methylated in DS vs TD samples.

Conclusions

A sex-specific signature of DNA methylation was detected in NDBS of DS-CHD compared to DS non-CHD individuals. This supports the hypothesis that epigenetics can reflect the variability of phenotypes in DS, particularly CHDs.

Epigenome-wide association study of bronchopulmonary dysplasia in preterm infants: results from the discovery-BPD program

BACKGROUND

Bronchopulmonary dysplasia (BPD) is a lung disease in premature infants caused by therapeutic oxygen supplemental and characterized by impaired pulmonary development which persists into later life. While advances in neonatal care have improved survival rates of premature infants, cases of BPD have been increasing with limited therapeutic options for prevention and treatment. This study was designed to explore the relationship between gestational age (GA), birth weight, and estimated blood cell-type composition in premature infants and to elucidate early epigenetic biomarkers associated with BPD.

METHODS

Cord blood DNA from preterm neonates that went on to develop BPD (n = 14) or not (non-BPD, n = 93) was applied to Illumina 450 K methylation arrays. Blood cell-type compositions were estimated using DNA methylation profiles. Multivariable robust regression analysis elucidated CpGs associated with BPD risk. cDNA microarray analysis of cord blood RNA identified differentially expressed genes in neonates who later developed BPD.

RESULTS

The development of BPD and the need for oxygen supplementation were strongly associated with GA (BPD, p < 1.0E-04; O2 supplementation, p < 1.0E-09) and birth weight (BPD, p < 1.0E-02; O2 supplementation, p < 1.0E-07). The estimated nucleated red blood cell (NRBC) percent was negatively associated with birth weight and GA, positively associated with hypomethylation of the tobacco smoke exposure biomarker cg05575921, and high-NRBC blood samples displayed a hypomethylation profile. Epigenome-wide association study (EWAS) identified 38 (Bonferroni) and 275 (false discovery rate 1%) differentially methylated CpGs associated with BPD. BPD-associated CpGs in cord blood were enriched for lung maturation and hematopoiesis pathways. Stochastic epigenetic mutation burden at birth was significantly elevated among those who developed BPD (adjusted p = 0.02). Transcriptome changes in cord blood cells reflected cell cycle, development, and pulmonary disorder events in BPD.

CONCLUSIONS

While results must be interpreted with caution because of the small size of this study, NRBC content strongly impacted DNA methylation profiles in preterm cord blood and EWAS analysis revealed potential insights into biological pathways involved in BPD pathogenesis.

Immune cell type and DNA methylation vary with reproductive status in women: possible pathways for costs of reproduction

Background

Consistent with evolutionarily theorized costs of reproduction (CoR), reproductive history in women is associated with life expectancy and susceptibility to certain cancers, autoimmune disorders and metabolic disease. Immunological changes originating during reproduction may help explain some of these relationships.

Methodology

To explore the potential role of the immune system in female CoR, we characterized leukocyte composition and regulatory processes using DNA methylation (DNAm) in a cross-sectional cohort of young (20–22 years old) women differing in reproductive status.

Results

Compared to nulliparity, pregnancy was characterized by differential methylation at 828 sites, 96% of which were hypomethylated and enriched for genes associated with T-cell activation, innate immunity, pre-eclampsia and neoplasia. Breastfeeding was associated with differential methylation at 1107 sites (71% hypermethylated), enriched for genes involved in metabolism, immune self-recognition and neurogenesis. There were no significant differences in DNAm between nulliparous and parous women. However, compared to nullipara, pregnant women had lower proportions of B, CD4T, CD8T and natural killer (NK) cells, and higher proportions of granulocytes and monocytes. Monocyte counts were lower and NK counts higher among breastfeeding women, and remained so among parous women.

Implications

Our findings point to widespread differences in DNAm during pregnancy and lactation. These effects appear largely transient, but may accumulate with gravidity become detectable as women age. Nulliparous and parous women differed in leukocyte composition, consistent with more persistent effects of reproduction on cell type. These findings support transient (leukocyte DNAm) and persistent (cell composition) changes associated with reproduction in women, illuminating potential pathways contributing to CoR.

Lay Summary: Evolutionary theory and epidemiology support costs of reproduction (CoR) to women’s health that may involve changes in immune function. We report differences in immune cell composition and gene regulation during pregnancy and breastfeeding. While many of these differences appear transient, immune cell composition may remain, suggesting mechanisms for female CoR.

Enhanced cell deconvolution of peripheral blood using DNA methylation for high-resolution immune profiling

DNA methylation microarrays can be employed to interrogate cell-type composition in complex tissues. Here, we expand reference-based deconvolution of blood DNA methylation to include 12 leukocyte subtypes (neutrophils, eosinophils, basophils, monocytes, naïve and memory B cells, naïve and memory CD4 + and CD8 + T cells, natural killer, and T regulatory cells). Including derived variables, our method provides 56 immune profile variables. The IDOL (IDentifying Optimal Libraries) algorithm was used to identify libraries for deconvolution of DNA methylation data for current and previous platforms. The accuracy of deconvolution estimates obtained using our enhanced libraries was validated using artificial mixtures and whole-blood DNA methylation with known cellular composition from flow cytometry. We applied our libraries to deconvolve cancer, aging, and autoimmune disease datasets. In conclusion, these libraries enable a detailed representation of immune-cell profiles in blood using only DNA and facilitate a standardized, thorough investigation of immune profiles in human health and disease.

Cell Types in Environmental Epigenetic Studies: Biological and Epidemiological Frameworks

PURPOSE OF REVIEW

This article introduces the roles of perinatal DNA methylation in human health and disease, highlights the challenges of tissue and cellular heterogeneity to studying DNA methylation, summarizes approaches to overcome these challenges, and offers recommendations in conducting research in environmental epigenetics.

RECENT FINDINGS

Epigenetic modifications are essential for human development and are labile to environmental influences, especially during gestation. Epigenetic dysregulation is also a hallmark of multiple diseases. Environmental epigenetic studies routinely measure DNA methylation in readily available tissues. However, tissues and cell types exhibit specific epigenetic patterning and heterogeneity between samples complicates epigenetic studies. Failure to account for cell-type heterogeneity limits identification of biological mechanisms and biases study results. Tissue-level epigenetic measures represent a convolution of epigenetic signals from individual cell types. Tissue-specific epigenetics is an evolving field and the use of disease-affected target, surrogate, or multiple tissues has inherent trade-offs and affects inference. Likewise, experimental and bioinformatic approaches to accommodate cell-type heterogeneity have varying assumptions and inherent trade-offs that affect inference. The relationships between exposure, disease, tissue-level DNA methylation, cell type-specific DNA methylation, and cell-type heterogeneity must be carefully considered in study design and analysis. Causal diagrams can inform study design and analytic strategies. Properly addressing cell-type heterogeneity limits sources of potential bias, avoids misinterpretation of study results, and allows investigators to distinguish shifts in cell-type proportions from direct changes to cellular epigenetic programming, both of which provide insights into environmental disease etiology and aid development of novel methods for prevention and treatment.

Global DNA methylation in placental tissues from pregnant with preeclampsia: A systematic review and pathway analysis

Pre-eclampsia (PE) is the major cause of fetal and maternal mortality and can be classified according to gestational age of onset into early-onset (EOPE, <34 weeks of gestation) and late- (LOPE, ≥34 weeks of gestation). DNA methylation (DNAm) may help to understand the abnormal placentation in PE. Therefore, we performed a systematic review to assess the role of global DNAm on pathophysiology of PE, focused on fetal and maternal tissues of placenta from pregnant with PE, including EOPE and LOPE. We searched the databases EMBASE, Medline/PubMed, Cochrane Central Register of Controlled Trials, Scopus, Lilacs, Scielo and Google Scholar, and followed the MOOSE guidelines. Moreover, we performed pathway analysis with the overlapping genes from the included studies. Twelve out of 24 included studies in the qualitative analysis considered the classification into EOPE and LOPE. We did not found heterogeneity in the criteria used for diagnosis of PE, and a few studies evaluated whether confounding factors would influence placental DNAm. Fourteen out of 24 included studies showed hypomethylation in placental tissue from pregnant with PE compared to controls. The differences in DNAm are specific to genes or differentially methylated regions, and more evident in EOPE and preterm PE compared to controls, rather than LOPE and term PE. The overlapping genes from included studies revealed pathways relevant to pathophysiology of PE. Our findings highlighted the heterogeneous results of the included studies, mainly focused on North America and China. Replication studies in different populations should use the same placental tissues, techniques to assess DNAm and pipelines for bioinformatic analysis.

The tissue-specific aspect of genome-wide DNA methylation in newborn and placental tissues: implications for epigenetic epidemiologic studies

Epigenetic programming is essential for lineage differentiation, embryogenesis and placentation in early pregnancy. In epigenetic association studies, DNA methylation is often examined in DNA derived from white blood cells, although its validity to other tissues of interest remains questionable. Therefore, we investigated the tissue specificity of epigenome-wide DNA methylation in newborn and placental tissues. Umbilical cord white blood cells (UC-WBC, n = 25), umbilical cord blood mononuclear cells (UC-MNC, n = 10), human umbilical vein endothelial cells (HUVEC, n = 25) and placental tissue (n = 25) were obtained from 36 uncomplicated pregnancies. Genome-wide DNA methylation was measured by the Illumina HumanMethylation450K BeadChip. Using UC-WBC as a reference tissue, we identified 3595 HUVEC tissue-specific differentially methylated regions (tDMRs) and 11,938 placental tDMRs. Functional enrichment analysis showed that HUVEC and placental tDMRs were involved in embryogenesis, vascular development and regulation of gene expression. No tDMRs were identified in UC-MNC. In conclusion, the extensive amount of genome-wide HUVEC and placental tDMRs underlines the relevance of tissue-specific approaches in future epigenetic association studies, or the use of validated representative tissues for a certain disease of interest, if available. To this purpose, we herewith provide a relevant dataset of paired, tissue-specific, genome-wide methylation measurements in newborn tissues.

Impact of ex vivo Sample Handling on DNA Methylation Profiles in Human Cord Blood and Neonatal Dried Blood Spots

The profiling of DNA methylation modifications in peripheral blood has significant potential to determine risk factors for human disease. Little is known concerning the sensitivity of DNA methylation profiles to ex vivo sample handling. Here, we studied typical conditions prior to sample storage associated with cord blood samples obtained from clinical investigations using reduced representation bisulfite sequencing. We examined both whole blood collected shortly after birth and dried blood spots, a potentially important source of neonatal blood for investigation of the DNA methylome and the Developmental Origins of Health and Disease in human cohorts because they are routinely collected during clinical care. Samples were matched across different time conditions, as they were from the same cord blood samples obtained from the same individuals. Maintaining whole blood ex vivo up to 24 h (4°C) or dried blood spots up to 7 days (room temp.) had little effect on DNA methylation profiles. Minimal differences were detected between cord blood immediately frozen and dried blood spots. Our results indicate that DNA methylation profiles are resilient to ex vivo sample handling conditions prior to storage. These data will help guide future human studies focused toward determination of DNA methylation modifications in whole blood.

Smoking-associated AHRR demethylation in cord blood DNA: impact of CD235a+ nucleated red blood cells

BACKGROUND

Numerous studies have demonstrated that DNA methylation levels in the aryl hydrocarbon receptor repressor (AHRR) gene measured in cord blood are significantly associated with prenatal tobacco smoke exposure and can be used as a fetal exposure biomarker. The mechanism driving this demethylation has not been determined and it is unclear if all cord blood cell types are impacted. Nucleated red blood cells (nRBCs/CD235a+ cells) are developmentally immature RBCs that display genome-wide hypomethylation and are observed at increased frequency in the cord blood of smoking mothers. We tested if AHRR methylation levels in CD235a+ nRBCs or nRBC counts influenced AHRR methylation in whole cord blood.

METHODS

Cord blood was collected from smoking (n = 34) and nonsmoking (n = 19) mothers and DNA was prepared from whole cord blood, isolated CD235a+ nRBCs, and CD14+ monocytes. AHRR methylation in cord blood DNA was measured using Illumina 850K arrays (cg05575921, chr5:373378). Pyrosequencing was used to compare methylation levels among cord blood, CD235a+, and CD14+ cells. We measured nRBC percentages using conventional complete blood counts and estimated percent nRBCs by a deconvolution model.

RESULTS

Methylation levels in AHRR were significantly lower in nRBCs relative to whole cord blood and CD14+ monocytes. While AHRR methylation levels in the cell types were significantly correlated across all subjects, methylation values at the chr5:373378 CpG averaged 14.6% lower in nRBCs (range 0.4 to 24.8%; p = 3.8E-13) relative to CD14+, with nonsmokers showing a significantly greater hypomethylation (- 4.1%, p = 1.8E-02). Methylation level at the AHRR chr5:373378 CpG was strongly associated with self-reported smoking in both CD14+ monocytes (t test p = 5.7E-09) and nRBCs (p = 4.8E-08), as well as cotinine levels (regression p = 1.1E-07 and p = 3.6E-04, respectively). For subjects with whole blood 850K data, robust linear regression models adjusting for estimated cell type composition, either including nRBCs counts or estimates, modestly increased the association between smoking and cg05575921 methylation.

CONCLUSIONS

Prenatal smoke exposure was highly significantly associated with AHRR methylation in cord blood, CD14+ monocytes, and CD235a+ nRBCs. AHRR methylation levels in nRBCs and nRBC counts had minimal effect on cord blood methylation measurements. However, regression models using estimated nRBCs or actual nRBC counts outperformed those lacking these covariates.

Locus-specific DNA methylation prediction in cord blood and placenta

DNA methylation is known to be responsive to prenatal exposures, which may be a part of the mechanism linking early developmental exposures to future chronic diseases. Many studies use blood to measure DNA methylation, yet we know that DNA methylation is tissue specific. Placenta is central to fetal growth and development, but it is rarely feasible to collect this tissue in large epidemiological studies; on the other hand, cord blood samples are more accessible. In this study, based on paired samples of both placenta and cord blood tissues from 169 individuals, we investigated the methylation concordance between placenta and cord blood. We then employed a machine-learning-based model to predict locus-specific DNA methylation levels in placenta using DNA methylation levels in cord blood. We found that methylation correlation between placenta and cord blood is lower than other tissue pairs, consistent with existing observations that placenta methylation has a distinct pattern. Nonetheless, there are still a number of CpG sites showing robust association between the two tissues. We built prediction models for placenta methylation based on cord blood data and documented a subset of 1,012 CpG sites with high correlation between measured and predicted placenta methylation levels. The resulting list of CpG sites and prediction models could help to reveal the loci where internal or external influences may affect DNA methylation in both placenta and cord blood, and provide a reference data to predict the effects on placenta in future study even when the tissue is not available in an epidemiological study.

Maternal swimming pool exposure during pregnancy in relation to birth outcomes and cord blood DNA methylation among private well users

Swimming in pools during pregnancy may expose the fetus to water disinfection by-products (DBP). As yet, our understanding of the impacts on DBPs on the fetus is uncertain. Individuals with public water systems are typically exposed to DBPs through drinking, showering and bathing, whereas among those on private water systems, swimming in pools may be the primary exposure source. We analyzed the effects of maternal swimming on birth outcomes and cord blood epigenetic changes in the New Hampshire Birth Cohort Study, a cohort of pregnant women with households on private water systems. Information about swimming in pools during pregnancy was obtained from 1033 women via questionnaires. Swimming pool use and duration were modeled using linear regression with newborn weight, length, and head circumference (z-scores) and genome wide cord blood DNA methylation as the outcomes and with adjustment for potential confounders. Overall 19.7% of women reported swimming in a pool during pregnancy. Among swimmers, duration of swimming was inversely related to head circumference (-0.02 z-score per 10% increase in duration, P = 0.004). No associations were observed with birth weight, length or DNA methylation modifications. Our findings suggest swimming pool exposure may impact the developing fetus although longer-term studies are needed.

An optimized library for reference-based deconvolution of whole-blood biospecimens assayed using the Illumina HumanMethylationEPIC BeadArray

Genome-wide methylation arrays are powerful tools for assessing cell composition of complex mixtures. We compare three approaches to select reference libraries for deconvoluting neutrophil, monocyte, B-lymphocyte, natural killer, and CD4+ and CD8+ T-cell fractions based on blood-derived DNA methylation signatures assayed using the Illumina HumanMethylationEPIC array. The IDOL algorithm identifies a library of 450 CpGs, resulting in an average R2 = 99.2 across cell types when applied to EPIC methylation data collected on artificial mixtures constructed from the above cell types. Of the 450 CpGs, 69% are unique to EPIC. This library has the potential to reduce unintended technical differences across array platforms.

Cell-type deconvolution from DNA methylation: a review of recent applications

Recent advances in cell-type deconvolution approaches are adding to our understanding of the biology underlying disease development and progression. DNA methylation (DNAm) can be used as a biomarker of cell types, and through deconvolution approaches, to infer underlying cell type proportions. Cell-type deconvolution algorithms have two main categories: reference-based and reference-free. Reference-based algorithms are supervised methods that determine the underlying composition of cell types within a sample by leveraging differentially methylated regions (DMRs) specific to cell type, identified from DNAm measures of purified cell populations. Reference-free algorithms are unsupervised methods for use when cell-type specific DMRs are not available, allowing scientists to estimate putative cellular proportions or control for potential confounding from cell type. Reference-based deconvolution is typically applied to blood samples and has potentiated our understanding of the relation between immune profiles and disease by allowing estimation of immune cell proportions from archival DNA. Bioinformatic analyses using DNAm to infer immune cell proportions, part of a new field known as Immunomethylomics, provides a new direction for consideration in epigenome wide association studies (EWAS).

Genome-wide DNA methylation drives human embryonic stem cell erythropoiesis by remodeling gene expression dynamics

AIM

To investigate the role of DNA methylation during erythrocyte production by human embryonic stem cells (hESCs).

METHODS

We employed an erythroid differentiation model from hESCs, and then tracked the genome-wide DNA methylation maps and gene expression patterns through an Infinium HumanMethylation450K BeadChip and an Ilumina Human HT-12 v4 Expression Beadchip, respectively.

RESULTS

A negative correlation between DNA methylation and gene expression was substantially enriched during the later differentiation stage and was present in both the promoter and the gene body. Moreover, erythropoietic genes with differentially methylated CpG sites that were primarily enriched in nonisland regions were upregulated, and demethylation of their gene bodies was associated with the presence of enhancers and DNase I hypersensitive sites. Finally, the components of JAK-STAT-NF-κB signaling were DNA hypomethylated and upregulated, which targets the key genes for erythropoiesis.

CONCLUSION

Erythroid lineage commitment by hESCs requires genome-wide DNA methylation modifications to remodel gene expression dynamics.

Cord blood hematopoietic cells from preterm infants display altered DNA methylation patterns

Background

Premature infants are highly vulnerable to infection. This is partly attributable to the preterm immune system, which differs from that of the term neonate in cell composition and function. Multiple studies have found differential DNA methylation (DNAm) between preterm and term infants’ cord blood; however, interpretation of these studies is limited by the confounding factor of blood cell composition. This study evaluates the epigenetic impact of preterm birth in isolated hematopoietic cell populations, reducing the concern of cell composition differences.

Methods

Genome-wide DNAm was measured using the Illumina 450K array in T cells, monocytes, granulocytes, and nucleated red blood cells (nRBCs) isolated from cord blood of 5 term and 5 preterm (<31 weeks gestational age) newborns. DNAm of hematopoietic cells was compared globally across the 450K array and through site-specific linear modeling.

Results

Nucleated red blood cells (nRBCs) showed the most extensive changes in DNAm, with 9258 differentially methylated (DM) sites (FDR < 5%, |Δβ| > 0.10) discovered between preterm and term infants compared to the <1000 prematurity-DM sites identified in white blood cell populations. The direction of DNAm change with gestational age at these prematurity-DM sites followed known patterns of hematopoietic differentiation, suggesting that term hematopoietic cell populations are more epigenetically mature than their preterm counterparts. Consistent shifts in DNAm between preterm and term cells were observed at 25 CpG sites, with many of these sites located in genes involved in growth and proliferation, hematopoietic lineage commitment, and the cytoskeleton. DNAm in preterm and term hematopoietic cells conformed to previously identified DNAm signatures of fetal liver and bone marrow, respectively.

Conclusions

This study presents the first genome-wide mapping of epigenetic differences in hematopoietic cells across the late gestational period. DNAm differences in hematopoietic cells between term and <31 weeks were consistent with the hematopoietic origin of these cells during ontogeny, reflecting an important role of DNAm in their regulation. Due to the limited sample size and the high coincidence of prematurity and multiple births, the relationship between cause of preterm birth and DNAm could not be evaluated. These findings highlight gene regulatory mechanisms at both cell-specific and systemic levels that may be involved in fetal immune system maturation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-017-0339-1) contains supplementary material, which is available to authorized users.