Differential DNA methylation in purified human blood cells: implications for cell lineage and studies on disease susceptibility

Animal and plant protein intake during infancy and childhood DNA methylation: a meta-analysis in the NutriPROGRAM consortium

BACKGROUND

Higher early-life animal protein intake is associated with a higher childhood obesity risk compared to plant protein intake. Differential DNA methylation may represent an underlying mechanism.

METHODS

We analysed associations of infant animal and plant protein intakes with DNA methylation in early (2-6 years, N = 579) and late (7̄-12 years, N = 604) childhood in two studies. Study-specific robust linear regression models adjusted for relevant confounders were run, and then meta-analysed using a fixed-effects model. We also performed sex-stratified meta-analyses. Follow-up analyses included pathway analysis and eQTM look-up.

RESULTS

Infant animal protein intake was not associated with DNA methylation in early childhood, but was associated with late-childhood DNA methylation at cg21300373 (P = 4.27 × 10¯8, MARCHF1) and cg10633363 (P = 1.09 × 10¯7, HOXB9) after FDR correction. Infant plant protein intake was associated with early-childhood DNA methylation at cg25973293 (P = 2.26 × 10-7, C1orf159) and cg15407373 (P = 2.13 × 10-7, MBP) after FDR correction. There was no overlap between the findings from the animal and plant protein analyses. We did not find enriched functional pathways at either time point using CpGs associated with animal and plant protein. These CpGs were not previously associated with childhood gene expression. Sex-stratified meta-analyses showed sex-specific DNA methylation associations for both animal and plant protein intake.

CONCLUSION

Infant animal protein intake was associated with DNA methylation at two CpGs in late childhood. Infant plant protein intake was associated with DNA methylation in early childhood at two CpGs. A potential mediating role of DNA methylation at these CpGs between infant protein intake and health outcomes requires further investigation.

Development of an epigenetic clock resistant to changes in immune cell composition

Epigenetic clocks are age predictors that use machine-learning models trained on DNA CpG methylation values to predict chronological or biological age. Increases in predicted epigenetic age relative to chronological age (epigenetic age acceleration) are connected to aging-associated pathologies, and changes in epigenetic age are linked to canonical aging hallmarks. However, epigenetic clocks rely on training data from bulk tissues whose cellular composition changes with age. Here, we found that human naive CD8+ T cells, which decrease in frequency during aging, exhibit an epigenetic age 15-20 years younger than effector memory CD8+ T cells from the same individual. Importantly, homogenous naive T cells isolated from individuals of different ages show a progressive increase in epigenetic age, indicating that current epigenetic clocks measure two independent variables, aging and immune cell composition. To isolate the age-associated cell intrinsic changes, we created an epigenetic clock, the IntrinClock, that did not change among 10 immune cell types tested. IntrinClock shows a robust predicted epigenetic age increase in a model of replicative senescence in vitro and age reversal during OSKM-mediated reprogramming.

An assessment of compositional methods for the analysis of DNA methylation-based deconvolution estimates

DNA methylation (DNAm)-based deconvolution estimates contain relative data, forming a composition, that standard methods (testing directly on cell proportions) are ill-suited to handle. In this study we examined the performance of an alternative method, analysis of compositions of microbiomes (ANCOM), for the analysis of DNAm-based deconvolution estimates. We performed two different simulation studies comparing ANCOM to a standard approach (two sample t-test performed directly on cell proportions) and analyzed a real-world data from the Women's Health Initiative to evaluate the applicability of ANCOM to DNAm-based deconvolution estimates. Our findings indicate that ANCOM can effectively account for the compositional nature of DNAm-based deconvolution estimates. ANCOM adequately controls the false discovery rate while maintaining statistical power comparable to that of standard methods.

Maternal Adverse Childhood Experiences and Biological Aging During Pregnancy and in Newborns

Importance

Adverse childhood experiences (ACEs), potentially traumatic experiences occurring before the age of 18 years, are associated with epigenetic aging later in life and may be transmitted across generations.

Objective

To test evidence of the transmission of biological embedding of life experience across generations by analyzing maternal ACEs and epigenetic clocks measured in mothers during pregnancy and in their children at birth.

Design, Setting, and Participants

For this cross-sectional study, data from the Accessible Resource for Integrated Epigenomic Studies (ARIES) substudy of the Avon Longitudinal Study of Parents and Children (ALSPAC) were analyzed. The ALSPAC study recruited 14 541 women who gave birth in the Avon Health District in the UK between April 1, 1991, and December 31, 1992. The ARIES substudy comprised 1018 mother-offspring dyads based on the availability of DNA samples profiled in 2014. Epigenetic age was estimated using DNA methylation-based epigenetic clocks (including Horvath, Hannum, GrimAge, PhenoAge, and DunedinPACE) in mothers during pregnancy and the Knight and Bohlin cord blood epigenetic clocks in newborns. Analyses were performed between October 1, 2022, and November 30, 2023.

Exposures

A composite measure of maternal ACEs was the primary exposure in both maternal and offspring models; as a secondary analysis, individual ACEs were measured separately. The Edinburgh Postnatal Depression Scale (EPDS) was used to investigate depression during pregnancy as an exposure.

Main Outcomes and Measures

Changes in epigenetic age acceleration (EAA) were investigated as the primary outcome in maternal models during pregnancy. Changes in epigenetic gestational age acceleration (GAA) were the primary outcome in offspring analyses. Linear regression analyses were used to determine the association between maternal ACEs and both outcomes.

Results

This study included 883 mother-child dyads. The mean (SD) maternal age at delivery was 29.8 (4.3) years. Pregnant women with higher ACE scores exhibited higher GrimAge EAA (β, 0.22 [95% CI, 0.12 to 0.33] years; P < .001). Maternal ACEs were not associated with GAA in newborns using P < .05 as a cutoff to determine statistical significance. Depression was associated with higher GrimAge EAA (β, 0.06 [95% CI, 0.02 to 0.10] years; P = .01) in mothers during pregnancy, but not in newborns, and did not mediate the association between ACEs and EAA.

Conclusions and Relevance

The findings of this study suggest that maternal ACEs may be associated with epigenetic aging later in life, including during pregnancy, supporting a role for maternal ACEs in offspring development and health later in life.

Biological aging of different blood cell types

Biological age (BA) captures detrimental age-related changes. The best-known and most-used BA indicators include DNA methylation-based epigenetic clocks and telomere length (TL). The most common biological sample material for epidemiological aging studies, whole blood, is composed of different cell types. We aimed to compare differences in BAs between blood cell types and assessed the BA indicators' cell type-specific associations with chronological age (CA). An analysis of DNA methylation-based BA indicators, including TL, methylation level at cg16867657 in ELOVL2, as well as the Hannum, Horvath, DNAmPhenoAge, and DunedinPACE epigenetic clocks, was performed on 428 biological samples of 12 blood cell types. BA values were different in the majority of the pairwise comparisons between cell types, as well as in comparison to whole blood (p < 0.05). DNAmPhenoAge showed the largest cell type differences, up to 44.5 years and DNA methylation-based TL showed the lowest differences. T cells generally had the "youngest" BA values, with differences across subsets, whereas monocytes had the "oldest" values. All BA indicators, except DunedinPACE, strongly correlated with CA within a cell type. Some differences such as DNAmPhenoAge-difference between naïve CD4 + T cells and monocytes were constant regardless of the blood donor's CA (range 20-80 years), while for DunedinPACE they were not. In conclusion, DNA methylation-based indicators of BA exhibit cell type-specific characteristics. Our results have implications for understanding the molecular mechanisms underlying epigenetic clocks and underscore the importance of considering cell composition when utilizing them as indicators for the success of aging interventions.

Epigenetic age oscillates during the day

Since their introduction, epigenetic clocks have been extensively used in aging, human disease, and rejuvenation studies. In this article, we report an intriguing pattern: epigenetic age predictions display a 24-h periodicity. We tested a circadian blood sample collection using 17 epigenetic clocks addressing different aspects of aging. Thirteen clocks exhibited significant oscillations with the youngest and oldest age estimates around midnight and noon, respectively. In addition, daily oscillations were consistent with the changes of epigenetic age across different times of day observed in an independant populational dataset. While these oscillations can in part be attributed to variations in white blood cell type composition, cell count correction methods might not fully resolve the issue. Furthermore, some epigenetic clocks exhibited 24-h periodicity even in the purified fraction of neutrophils pointing at plausible contributions of intracellular epigenomic oscillations. Evidence for circadian variation in epigenetic clocks emphasizes the importance of the time-of-day for obtaining accurate estimates of epigenetic age.

Socioeconomic Status, Lifestyle, and DNA Methylation Age Among Racially and Ethnically Diverse Adults: NIMHD Social Epigenomics Program

Importance

Variation in DNA methylation at specific loci estimates biological age, which is associated with morbidity, mortality, and social experiences. Aging estimates known as epigenetic clocks, including the Dunedin Pace of Aging Calculated From the Epigenome (DunedinPACE), were trained on data predominately from individuals of European ancestry; however, limited research has explored DunedinPACE in underrepresented populations experiencing health disparities.

Objective

To investigate associations of neighborhood and individual sociobehavioral factors with biological aging in a racially and ethnically diverse population.

Design, Setting, and Participants

This cohort study, part of the Multiethnic Cohort study conducted from May 1993 to September 1996 to examine racial and ethnic disparities in chronic diseases, integrated biospecimen and self-reported data collected between April 2004 and November 2005 from healthy Hawaii residents aged 45 to 76 years. These participants self-identified as of Japanese American, Native Hawaiian, or White racial and ethnic background. Data were analyzed from January 2022 to May 2024.

Main Outcomes and Measures

DNA methylation data were generated from monocytes enriched from cryopreserved lymphocytes and used to derive DunedinPACE scores from November 2017 to June 2021. Neighborhood social economic status (NSES) was estimated from 1990 US Census Bureau data to include factors such as educational level, occupation, and income. Individual-level factors analyzed included educational level, body mass index (BMI), physical activity (PA), and diet quality measured by the Healthy Eating Index (HEI). Linear regression analysis of DunedinPACE scores was used to examine their associations with NSES and sociobehavioral variables.

Results

A total of 376 participants were included (113 [30.1%] Japanese American, 144 [38.3%] Native Hawaiian, and 119 [31.6%] White; 189 [50.3%] were female). Mean (SE) age was 57.81 (0.38) years. Overall, mean (SE) DunedinPACE scores were significantly higher among females than among males (1.28 [0.01] vs 1.25 [0.01]; P = .005); correlated negatively with NSES (R = -0.09; P = .08), HEI (R = -0.11; P = .03), and educational attainment (R = -0.15; P = .003) and positively with BMI (R = 0.31; P < .001); and varied by race and ethnicity. Native Hawaiian participants exhibited a higher mean (SE) DunedinPACE score (1.31 [0.01]) compared with Japanese American (1.25 [0.01]; P < .001) or White (1.22 [0.01]; P < .001) participants. Controlling for age, sex, HEI, BMI, and NSES, linear regression analyses revealed a negative association between educational level and DunedinPACE score among Japanese American (β, -0.005 [95% CI, -0.013 to 0.002]; P = .03) and Native Hawaiian (β, -0.003 [95% CI, -0.011 to 0.005]; P = .08) participants, yet this association was positive among White participants (β, 0.007; 95% CI, -0.001 to 0.015; P = .09). Moderate to vigorous PA was associated with lower DunedinPACE scores only among Native Hawaiian participants (β, -0.006; 95% CI, -0.011 to -0.001; P = .005), independent of NSES.

Conclusions and Relevance

In this study of a racially and ethnically diverse sample of 376 adults, low NSES was associated with a higher rate of biological aging measured by DunedinPACE score, yet individual-level factors such as educational level and physical activity affected this association, which varied by race and ethnicity. These findings support sociobehavioral interventions in addressing health inequities.

Maximizing Insights from Longitudinal Epigenetic Age Data: Simulations, Applications, and Practical Guidance

Background

Epigenetic Age (EA) is an age estimate, developed using DNA methylation (DNAm) states of selected CpG sites across the genome. Although EA and chronological age are highly correlated, EA may not increase uniformly with time. Departures, known as epigenetic age acceleration (EAA), are common and have been linked to various traits and future disease risk. Limited by available data, most studies investigating these relationships have been cross-sectional - using a single EA measurement. However, the recent growth in longitudinal DNAm studies has led to analyses of associations with EA over time. These studies differ in (i) their choice of model; (ii) the primary outcome (EA vs. EAA); and (iii) in their use of chronological age or age-independent time variables to account for the temporal dynamic. We evaluated the robustness of each approach using simulations and tested our results in two real-world examples, using biological sex and birthweight as predictors of longitudinal EA.

Results

Our simulations showed most accurate effect sizes in a linear mixed model or generalized estimating equation, using chronological age as the time variable. The use of EA versus EAA as an outcome did not strongly impact estimates. Applying the optimal model in real-world data uncovered an accelerated EA rate in males and an advanced EA that decelerates over time in children with higher birthweight.

Conclusion

Our results can serve as a guide for forthcoming longitudinal EA studies, aiding in methodological decisions that may determine whether an association is accurately estimated, overestimated, or potentially overlooked.

Neural-net-based cell deconvolution from DNA methylation reveals tumor microenvironment associated with cancer prognosis

DNA methylation is a pivotal epigenetic modification that defines cellular identity. While cell deconvolution utilizing this information is considered useful for clinical practice, current methods for deconvolution are limited in their accuracy and resolution. In this study, we collected DNA methylation data from 945 human samples derived from various tissues and tumor-infiltrating immune cells and trained a neural network model with them. The model, termed MEnet, predicted abundance of cell population together with the detailed immune cell status from bulk DNA methylation data, and showed consistency to those of flow cytometry and histochemistry. MEnet was superior to the existing methods in the accuracy, speed, and detectable cell diversity, and could be applicable for peripheral blood, tumors, cell-free DNA, and formalin-fixed paraffin-embedded sections. Furthermore, by applying MEnet to 72 intrahepatic cholangiocarcinoma samples, we identified immune cell profiles associated with cancer prognosis. We believe that cell deconvolution by MEnet has the potential for use in clinical settings.

Epigenome-wide analysis of frailty: Results from two European twin cohorts

Epigenetics plays an important role in the aging process, but it is unclear whether epigenetic factors also influence frailty, an age-related state of physiological decline. In this study, we performed a meta-analysis of epigenome-wide association studies in four samples drawn from the Swedish Adoption/Twin Study of Aging (SATSA) and the Longitudinal Study of Aging Danish Twins (LSADT) to explore the association between DNA methylation and frailty. Frailty was defined using the frailty index (FI), and DNA methylation levels were measured in whole blood using Illumina's Infinium HumanMethylation450K and MethylationEPIC arrays. In the meta-analysis consisting of a total of 829 participants, we identified 589 CpG sites that were statistically significantly associated with either the continuous or categorical FI (false discovery rate <0.05). Many of these CpGs have previously been associated with age and age-related diseases. The identified sites were also largely directionally consistent in a longitudinal analysis using mixed-effects models in SATSA, where the participants were followed up to a maximum of 20 years. Moreover, we identified three differentially methylated regions within the MGRN1, MIR596, and TAPBP genes that have been linked to neuronal aging, tumor growth, and immune functions. Furthermore, our meta-analysis results replicated 34 of the 77 previously reported frailty-associated CpGs at p < 0.05. In conclusion, our findings demonstrate robust associations between frailty and DNA methylation levels in 589 novel CpGs, previously unidentified for frailty, and strengthen the role of neuronal/brain pathways in frailty.

Severe induction of aberrant DNA methylation by nodular gastritis in adults

BACKGROUND

Nodular gastritis (NG) is characterized by marked antral lymphoid follicle formation, and is a strong risk factor for diffuse-type gastric cancer in adults. However, it is unknown whether aberrant DNA methylation, which is induced by atrophic gastritis (AG) and is a risk for gastric cancer, is induced by NG. Here, we analyzed methylation induction by NG.

METHODS

Gastric mucosal samples were obtained from non-cancerous antral tissues of 16 NG and 20 AG patients with gastric cancer and 5 NG and 6 AG patients without, all age- and gender-matched. Genome-wide methylation analysis and expression analysis were conducted by a BeadChip array and RNA-sequencing, respectively.

RESULTS

Clustering analysis of non-cancerous antral tissues of NG and AG patients with gastric cancer was conducted using methylation levels of 585 promoter CpG islands (CGIs) of methylation-resistant genes, and a large fraction of NG samples formed a cluster with strong methylation induction. Promoter CGIs of CDH1 and DAPK1 tumor-suppressor genes were more methylated in NG than in AG. Notably, methylation levels of these genes were also higher in the antrum of NG patients without cancer. Genes related to lymphoid follicle formation, such as CXCL13/CXCR5 and CXCL12/CXCR4, had higher expression in NG, and genes involved in DNA demethylation TET2 and IDH1, had only half the expression in NG.

CONCLUSIONS

Severe aberrant methylation, involving multiple tumor-suppressor genes, was induced in the gastric antrum and body of patients with NG, in accordance with their high gastric cancer risk.

Demographic characteristics and epigenetic biological aging among post-9/11 veterans: Associations of DunedinPACE with sex, race, and age

Measures of epigenetic aging derived from DNA methylation (DNAm) have enabled the assessment of biological aging in new populations and cohorts. In the present study, we used an epigenetic measure of aging, DunedinPACE, to examine rates of aging across demographic groups in a sample of 2,309 United States military veterans from the VISN 6 MIRECC's Post-Deployment Mental Health Study. As assessed by DunedinPACE, female veterans were aging faster than male veterans (β = 0.39, 95 % CI [0.29, 0.48], p < .001), non-Hispanic Black veterans were aging faster than non-Hispanic White veterans (β = 0.58, 95 % CI [0.50, 0.66], p < .001), and older veterans were biologically aging faster than younger veterans (β = 0.21, 95 % CI [0.18, 0.25], p < .001). In secondary analyses, these differences in rates of aging were not explained by a variety of biopsychosocial covariates. In addition, the percentage of European genetic admixture in non-Hispanic Black veterans was not associated with DunedinPACE. Our findings suggest that female and non-Hispanic Black veterans are at greater risk of accelerated aging among post-9/11 veterans. Interventions that slow aging might provide relatively greater benefit among veterans comprising these at-risk groups.

Brain cell-type shifts in Alzheimer's disease, autism, and schizophrenia interrogated using methylomics and genetics

Few neuropsychiatric disorders have replicable biomarkers, prompting high-resolution and large-scale molecular studies. However, we still lack consensus on a more foundational question: whether quantitative shifts in cell types-the functional unit of life-contribute to neuropsychiatric disorders. Leveraging advances in human brain single-cell methylomics, we deconvolve seven major cell types using bulk DNA methylation profiling across 1270 postmortem brains, including from individuals diagnosed with Alzheimer's disease, schizophrenia, and autism. We observe and replicate cell-type compositional shifts for Alzheimer's disease (endothelial cell loss), autism (increased microglia), and schizophrenia (decreased oligodendrocytes), and find age- and sex-related changes. Multiple layers of evidence indicate that endothelial cell loss contributes to Alzheimer's disease, with comparable effect size to APOE genotype among older people. Genome-wide association identified five genetic loci related to cell-type composition, involving plausible genes for the neurovascular unit (P2RX5 and TRPV3) and excitatory neurons (DPY30 and MEMO1). These results implicate specific cell-type shifts in the pathophysiology of neuropsychiatric disorders.

Current Challenges of Methylation-Based Liquid Biopsies in Cancer Diagnostics

In current clinical practice, effective cancer testing and screening paradigms are limited to specific types of cancer, exhibiting varying efficiency, acceptance, and adherence. Cell-free DNA (cfDNA) methylation profiling holds promise in providing information about the presence of malignity regardless of its type and location while leveraging blood-based liquid biopsies as a method to obtain analytical samples. However, technical difficulties, costs and challenges resulting from biological variations, tumor heterogeneity, and exogenous factors persist. This method exploits the mechanisms behind cfDNA release but faces issues like fragmentation, low concentrations, and high background noise. This review explores cfDNA methylation's origins, means of detection, and profiling for cancer diagnostics. The critical evaluation of currently available multi-cancer early detection methods (MCEDs) as well as tests targeting single genes, emphasizing their potential and limits to refine strategies for early cancer detection, are explained. The current methodology limitations, workflows, comparisons of clinically approved liquid biopsy-based methylation tests for cancer, their utilization in companion diagnostics as well as the biological limitations of the epigenetics approach are discussed, aiming to help healthcare providers as well as researchers to orient themselves in this increasingly complex and evolving field of diagnostics.

Validating reference-based algorithms to determine cell-type heterogeneity in ovarian cancer DNA methylation studies

Information about cell composition in tissue samples is crucial for biomarker discovery and prognosis. Specifically, cancer tissue samples present challenges in deconvolution studies due to mutations and genetic rearrangements. Here, we optimized a robust, DNA methylation-based protocol, to be used for deconvolution of ovarian cancer samples. We compared several state-of-the-art methods (HEpiDISH, MethylCIBERSORT and ARIC) and validated the proposed protocol in an in-silico mixture and in an external dataset containing samples from ovarian cancer patients and controls. The deconvolution protocol we eventually implemented is based on MethylCIBERSORT. Comparing deconvolution methods, we paid close attention to the role of a reference panel. We postulate that a possibly high number of samples (in our case: 247) should be used when building a reference panel to ensure robustness and to compensate for biological and technical variation between samples. Subsequently, we tested the performance of the validated protocol in our own study cohort, consisting of 72 patients with malignant and benign ovarian disease as well as in five external cohorts. In conclusion, we refined and validated a reference-based algorithm to determine cell type composition of ovarian cancer tissue samples to be used in cancer biology studies in larger cohorts.

DNA Methylation-Derived Immune Cell Proportions and Cancer Risk, Including Lung Cancer, in Black Participants

Prior cohort studies assessing cancer risk based on immune cell subtype profiles have predominantly focused on White populations. This limitation obscures vital insights into how cancer risk varies across race. Immune cell subtype proportions were estimated using deconvolution based on leukocyte DNA methylation markers from blood samples collected at baseline on participants without cancer in the Atherosclerosis Risk in Communities (ARIC) Study. Over a mean of 17.5 years of follow-up, 668 incident cancers were diagnosed in 2,467 Black participants. Cox proportional hazards regression was used to examine immune cell subtype proportions and overall cancer incidence and site-specific incidence (lung, breast, and prostate cancers). Higher T regulatory cell proportions were associated with statistically significantly higher lung cancer risk (hazard ratio = 1.22, 95% confidence interval = 1.06-1.41 per percent increase). Increased memory B cell proportions were associated with significantly higher risk of prostate cancer (1.17, 1.04-1.33) and all cancers (1.13, 1.05-1.22). Increased CD8+ naïve cell proportions were associated with significantly lower risk of all cancers in participants ≥55 years (0.91, 0.83-0.98). Other immune cell subtypes did not display statistically significant associations with cancer risk. These results in Black participants align closely with prior findings in largely White populations. Findings from this study could help identify those at high cancer risk and outline risk stratifying to target patients for cancer screening, prevention, and other interventions. Further studies should assess these relationships in other cancer types, better elucidate the interplay of B cells in cancer risk, and identify biomarkers for personalized risk stratification.

Associations in cell type-specific hydroxymethylation and transcriptional alterations of pediatric central nervous system tumors

Although intratumoral heterogeneity has been established in pediatric central nervous system tumors, epigenomic alterations at the cell type level have largely remained unresolved. To identify cell type-specific alterations to cytosine modifications in pediatric central nervous system tumors, we utilize a multi-omic approach that integrated bulk DNA cytosine modification data (methylation and hydroxymethylation) with both bulk and single-cell RNA-sequencing data. We demonstrate a large reduction in the scope of significantly differentially modified cytosines in tumors when accounting for tumor cell type composition. In the progenitor-like cell types of tumors, we identify a preponderance differential Cytosine-phosphate-Guanine site hydroxymethylation rather than methylation. Genes with differential hydroxymethylation, like histone deacetylase 4 and insulin-like growth factor 1 receptor, are associated with cell type-specific changes in gene expression in tumors. Our results highlight the importance of epigenomic alterations in the progenitor-like cell types and its role in cell type-specific transcriptional regulation in pediatric central nervous system tumors.

Green space exposure and blood DNA methylation at birth and in childhood - A multi-cohort study

Green space exposure has been associated with improved mental, physical and general health. However, the underlying biological mechanisms remain largely unknown. The aim of this study was to investigate the association between green space exposure and cord and child blood DNA methylation. Data from eight European birth cohorts with a total of 2,988 newborns and 1,849 children were used. Two indicators of residential green space exposure were assessed: (i) surrounding greenness (satellite-based Normalized Difference Vegetation Index (NDVI) in buffers of 100 m and 300 m) and (ii) proximity to green space (having a green space ≥ 5,000 m2 within a distance of 300 m). For these indicators we assessed two exposure windows: (i) pregnancy, and (ii) the period from pregnancy to child blood DNA methylation assessment, named as cumulative exposure. DNA methylation was measured with the Illumina 450K or EPIC arrays. To identify differentially methylated positions (DMPs) we fitted robust linear regression models between pregnancy green space exposure and cord blood DNA methylation and between cumulative green space exposure and child blood DNA methylation. Two sensitivity analyses were conducted: (i) without adjusting for cellular composition, and (ii) adjusting for air pollution. Cohort results were combined through fixed-effect inverse variance weighted meta-analyses. Differentially methylated regions (DMRs) were identified from meta-analysed results using the Enmix-combp and DMRcate methods. There was no statistical evidence of pregnancy or cumulative exposures associating with any DMP (False Discovery Rate, FDR, p-value < 0.05). However, surrounding greenness exposure was inversely associated with four DMRs (three in cord blood and one in child blood) annotated to ADAMTS2, KCNQ1DN, SLC6A12 and SDK1 genes. Results did not change substantially in the sensitivity analyses. Overall, we found little evidence of the association between green space exposure and blood DNA methylation. Although we identified associations between surrounding greenness exposure with four DMRs, these findings require replication.

Donor whole blood DNA methylation is not a strong predictor of acute graft versus host disease in unrelated donor allogeneic haematopoietic cell transplantation

Allogeneic hematopoietic cell transplantation (HCT) is used to treat many blood-based disorders and malignancies, however it can also result in serious adverse events, such as the development of acute graft-versus-host disease (aGVHD). This study aimed to develop a donor-specific epigenetic classifier to reduce incidence of aGVHD by improving donor selection. Genome-wide DNA methylation was assessed in a discovery cohort of 288 HCT donors selected based on recipient aGVHD outcome; this cohort consisted of 144 cases with aGVHD grades III-IV and 144 controls with no aGVHD. We applied a machine learning algorithm to identify CpG sites predictive of aGVHD. Receiver operating characteristic (ROC) curve analysis of these sites resulted in a classifier with an encouraging area under the ROC curve (AUC) of 0.91. To test this classifier, we used an independent validation cohort (n = 288) selected using the same criteria as the discovery cohort. Attempts to validate the classifier failed with the AUC falling to 0.51. These results indicate that donor DNA methylation may not be a suitable predictor of aGVHD in an HCT setting involving unrelated donors, despite the initial promising results in the discovery cohort. Our work highlights the importance of independent validation of machine learning classifiers, particularly when developing classifiers intended for clinical use.

The transition from genomics to phenomics in personalized population health

Modern health care faces several serious challenges, including an ageing population and its inherent burden of chronic diseases, rising costs and marginal quality metrics. By assessing and optimizing the health trajectory of each individual using a data-driven personalized approach that reflects their genetics, behaviour and environment, we can start to address these challenges. This assessment includes longitudinal phenome measures, such as the blood proteome and metabolome, gut microbiome composition and function, and lifestyle and behaviour through wearables and questionnaires. Here, we review ongoing large-scale genomics and longitudinal phenomics efforts and the powerful insights they provide into wellness. We describe our vision for the transformation of the current health care from disease-oriented to data-driven, wellness-oriented and personalized population health.

Computational deconvolution of DNA methylation data from mixed DNA samples

In this review, we provide a comprehensive overview of the different computational tools that have been published for the deconvolution of bulk DNA methylation (DNAm) data. Here, deconvolution refers to the estimation of cell-type proportions that constitute a mixed sample. The paper reviews and compares 25 deconvolution methods (supervised, unsupervised or hybrid) developed between 2012 and 2023 and compares the strengths and limitations of each approach. Moreover, in this study, we describe the impact of the platform used for the generation of methylation data (including microarrays and sequencing), the applied data pre-processing steps and the used reference dataset on the deconvolution performance. Next to reference-based methods, we also examine methods that require only partial reference datasets or require no reference set at all. In this review, we provide guidelines for the use of specific methods dependent on the DNA methylation data type and data availability.

A subtype of schizophrenia patients with altered methylation level of genes related to immune cell activity

BACKGROUND

Epigenetic changes are plausible molecular sources of clinical heterogeneity in schizophrenia. A subgroup of schizophrenia patients with elevated inflammatory or immune-dysregulation has been reported by previous studies. However, little is known about epigenetic changes in genes related to immune activation in never-treated first-episode patients with schizophrenia (FES) and its consistency with that in treated long-term ill (LTS) patients.

METHODS

In this study, epigenome-wide profiling with a DNA methylation array was applied using blood samples of both FES and LTS patients, as well as their corresponding healthy controls. Non-negative matrix factorization (NMF) and k -means clustering were performed to parse heterogeneity of schizophrenia, and the consistency of subtyping results from two cohorts. was tested.

RESULTS

This study identified a subtype of patients in FES participants (47.5%) that exhibited widespread methylation level alterations of genes enriched in immune cell activity and a significantly higher proportion of neutrophils. This clustering of FES patients was validated in LTS patients, with high correspondence in epigenetic and clinical features across two cohorts.

CONCLUSIONS

In summary, this study demonstrated a subtype of schizophrenia patients across both FES and LTS cohorts, defined by widespread alterations in methylation profile of genes related to immune function and distinguishing clinical features. This finding illustrates the promise of novel treatment strategies targeting immune dysregulation for a subpopulation of schizophrenia patients.

The Epigenetic Evolution of Glioma Is Determined by the IDH1 Mutation Status and Treatment Regimen

Tumor adaptation or selection is thought to underlie therapy resistance in glioma. To investigate longitudinal epigenetic evolution of gliomas in response to therapeutic pressure, we performed an epigenomic analysis of 132 matched initial and recurrent tumors from patients with IDH-wildtype (IDHwt) and IDH-mutant (IDHmut) glioma. IDHwt gliomas showed a stable epigenome over time with relatively low levels of global methylation. The epigenome of IDHmut gliomas showed initial high levels of genome-wide DNA methylation that was progressively reduced to levels similar to those of IDHwt tumors. Integration of epigenomics, gene expression, and functional genomics identified HOXD13 as a master regulator of IDHmut astrocytoma evolution. Furthermore, relapse of IDHmut tumors was accompanied by histologic progression that was associated with survival, as validated in an independent cohort. Finally, the initial cell composition of the tumor microenvironment varied between IDHwt and IDHmut tumors and changed differentially following treatment, suggesting increased neoangiogenesis and T-cell infiltration upon treatment of IDHmut gliomas. This study provides one of the largest cohorts of paired longitudinal glioma samples with epigenomic, transcriptomic, and genomic profiling and suggests that treatment of IDHmut glioma is associated with epigenomic evolution toward an IDHwt-like phenotype.

SIGNIFICANCE

Standard treatments are related to loss of DNA methylation in IDHmut glioma, resulting in epigenetic activation of genes associated with tumor progression and alterations in the microenvironment that resemble treatment-naïve IDHwt glioma.

Genome-wide DNA methylation analysis of cannabis use disorder in a veteran cohort enriched for posttraumatic stress disorder

Cannabis use has been increasing over the past decade, not only in the general US population, but particularly among military veterans. With this rise in use has come a concomitant increase in cannabis use disorder (CUD) among veterans. Here, we performed an epigenome-wide association study for lifetime CUD in an Iraq/Afghanistan era veteran cohort enriched for posttraumatic stress disorder (PTSD) comprising 2,310 total subjects (1,109 non-Hispanic black and 1,201 non-Hispanic white). We also investigated CUD interactions with current PTSD status and examined potential indirect effects of DNA methylation (DNAm) on the relationship between CUD and psychiatric diagnoses. Four CpGs were associated with lifetime CUD, even after controlling for the effects of current smoking (AHRR cg05575921, LINC00299 cg23079012, VWA7 cg22112841, and FAM70A cg08760398). Importantly, cg05575921, a CpG strongly linked to smoking, remained associated with lifetime CUD even when restricting the analysis to veterans who reported never smoking cigarettes. Moreover, CUD interacted with current PTSD to affect cg05575921 and cg23079012 such that those with both CUD and PTSD displayed significantly lower DNAm compared to the other groups. Finally, we provide preliminary evidence that AHRR cg05575921 helps explain the association between CUD and any psychiatric diagnoses, specifically mood disorders.

Associations of prenatal one-carbon metabolism nutrients and metals with epigenetic aging biomarkers at birth and in childhood in a US cohort

Epigenetic gestational age acceleration (EGAA) at birth and epigenetic age acceleration (EAA) in childhood may be biomarkers of the intrauterine environment. We investigated the extent to which first-trimester folate, B12, 5 essential, and 7 non-essential metals in maternal circulation are associated with EGAA and EAA in early life. Bohlin EGAA and Horvath pan-tissue and skin and blood EAA were calculated using DNA methylation measured in cord blood (N=351) and mid-childhood blood (N=326; median age = 7.7 years) in the Project Viva pre-birth cohort. A one standard deviation increase in individual essential metals (copper, manganese, and zinc) was associated with 0.94-1.2 weeks lower Horvath EAA at birth, and patterns of exposures identified by exploratory factor analysis suggested that a common source of essential metals was associated with Horvath EAA. We also observed evidence nonlinear associations of zinc with Bohlin EGAA, magnesium and lead with Horvath EAA, and cesium with skin and blood EAA at birth. Overall, associations at birth did not persist in mid-childhood; however, arsenic was associated with greater EAA at birth and in childhood. Prenatal metals, including essential metals and arsenic, are associated with epigenetic aging in early life, which might be associated with future health.

Blood-based DNA methylation and exposure risk scores predict PTSD with high accuracy in military and civilian cohorts

Background

Incorporating genomic data into risk prediction has become an increasingly useful approach for rapid identification of individuals most at risk for complex disorders such as PTSD. Our goal was to develop and validate Methylation Risk Scores (MRS) using machine learning to distinguish individuals who have PTSD from those who do not.

Methods

Elastic Net was used to develop three risk score models using a discovery dataset (n = 1226; 314 cases, 912 controls) comprised of 5 diverse cohorts with available blood-derived DNA methylation (DNAm) measured on the Illumina Epic BeadChip. The first risk score, exposure and methylation risk score (eMRS) used cumulative and childhood trauma exposure and DNAm variables; the second, methylation-only risk score (MoRS) was based solely on DNAm data; the third, methylation-only risk scores with adjusted exposure variables (MoRSAE) utilized DNAm data adjusted for the two exposure variables. The potential of these risk scores to predict future PTSD based on pre-deployment data was also assessed. External validation of risk scores was conducted in four independent cohorts.

Results

The eMRS model showed the highest accuracy (92%), precision (91%), recall (87%), and f1-score (89%) in classifying PTSD using 3730 features. While still highly accurate, the MoRS (accuracy = 89%) using 3728 features and MoRSAE (accuracy = 84%) using 4150 features showed a decline in classification power. eMRS significantly predicted PTSD in one of the four independent cohorts, the BEAR cohort (beta = 0.6839, p-0.003), but not in the remaining three cohorts. Pre-deployment risk scores from all models (eMRS, beta = 1.92; MoRS, beta = 1.99 and MoRSAE, beta = 1.77) displayed a significant (p < 0.001) predictive power for post-deployment PTSD.

Conclusion

Results, especially those from the eMRS, reinforce earlier findings that methylation and trauma are interconnected and can be leveraged to increase the correct classification of those with vs. without PTSD. Moreover, our models can potentially be a valuable tool in predicting the future risk of developing PTSD. As more data become available, including additional molecular, environmental, and psychosocial factors in these scores may enhance their accuracy in predicting the condition and, relatedly, improve their performance in independent cohorts.

Maternal educational attainment in pregnancy and epigenome-wide DNA methylation changes in the offspring from birth until adolescence

Maternal educational attainment (MEA) shapes offspring health through multiple potential pathways. Differential DNA methylation may provide a mechanistic understanding of these long-term associations. We aimed to quantify the associations of MEA with offspring DNA methylation levels at birth, in childhood and in adolescence. Using 37 studies from high-income countries, we performed meta-analysis of epigenome-wide association studies (EWAS) to quantify the associations of completed years of MEA at the time of pregnancy with offspring DNA methylation levels at birth (n = 9 881), in childhood (n = 2 017), and adolescence (n = 2 740), adjusting for relevant covariates. MEA was found to be associated with DNA methylation at 473 cytosine-phosphate-guanine sites at birth, one in childhood, and four in adolescence. We observed enrichment for findings from previous EWAS on maternal folate, vitamin-B12 concentrations, maternal smoking, and pre-pregnancy BMI. The associations were directionally consistent with MEA being inversely associated with behaviours including smoking and BMI. Our findings form a bridge between socio-economic factors and biology and highlight potential pathways underlying effects of maternal education. The results broaden our understanding of bio-social associations linked to differential DNA methylation in multiple early stages of life. The data generated also offers an important resource to help a more precise understanding of the social determinants of health.

Twin study: genotype-dependent epigenetic factors affecting free thyroxine levels in the normal range

Aim: To explore the clinical application of DNA methylation affecting thyroid function, we evaluated the association of DNA methylation with free thyroxine (FT4) and TSH measurements in monozygotic twins. Materials & methods: Discordant pairs for FT4 or TSH levels were examined for the relationship between the within-pair difference of each measurement and the DNA methylation levels using epigenome-wide association studies. The contribution of polymorphisms to the methylation sensitivity was also examined. Results: We found two CpG sites significantly associated with FT4 levels, and also some CpG sites showing significant differences in their methylation levels within FT4-discordant pairs depending on the polymorphism in EPHB2. Conclusion: The FT4 level may be associated with a combination of methylation and polymorphisms in the EPHB2 gene.

Maternal prenatal social experiences and offspring epigenetic age acceleration from birth to mid-childhood

PURPOSE

Investigate associations of maternal social experiences with offspring epigenetic age acceleration (EAA) from birth through mid-childhood among 205 mother-offspring dyads of minoritized racial and ethnic groups.

METHODS

We used linear regression to examine associations of maternal experiences of racial bias or discrimination (0 = none, 1-2 = intermediate, or 3+ = high), social support (tertile 1 = low, 2 = intermediate, 3 = high), and socioeconomic status index (tertile 1 = low, 2 = intermediate, 3 = high) during the prenatal period with offspring EAA according to Horvath's Pan-Tissue, Horvath's Skin and Blood, and Intrinsic EAA clocks at birth, 3 years, and 7 years.

RESULTS

In comparison to children of women who did not experience any racial bias or discrimination, those whose mothers reported highest levels of racial bias or discrimination had lower Pan-Tissue clock EAA in early (-0.50 years; 90% CI: -0.91, -0.09) and mid-childhood (-0.75 years; -1.41, -0.08). We observed similar associations for the Skin and Blood clock and Intrinsic EAA. Maternal experiences of discrimination were not associated with Pan-Tissue EAA at birth. Neither maternal social support nor socioeconomic status predicted offspring EAA.

CONCLUSIONS

Children whose mothers experienced higher racial bias or discrimination exhibited slower EAA. Future studies are warranted to confirm these findings and establish associations of early-life EAA with long-term health outcomes.

Comprehensive evaluation of smoking exposures and their interactions on DNA methylation

BACKGROUND

Smoking impacts DNA methylation, but data are lacking on smoking-related differential methylation by sex or dietary intake, recent smoking cessation (<1 year), persistence of differential methylation from in utero smoking exposure, and effects of environmental tobacco smoke (ETS).

METHODS

We meta-analysed data from up to 15,014 adults across 5 cohorts with DNA methylation measured in blood using Illumina's EPIC array for current smoking (2560 exposed), quit < 1 year (500 exposed), in utero (286 exposed), and ETS exposure (676 exposed). We also evaluated the interaction of current smoking with sex or diet (fibre, folate, and vitamin C).

FINDINGS

Using false discovery rate (FDR < 0.05), 65,857 CpGs were differentially methylated in relation to current smoking, 4025 with recent quitting, 594 with in utero exposure, and 6 with ETS. Most current smoking CpGs attenuated within a year of quitting. CpGs related to in utero exposure in adults were enriched for those previously observed in newborns. Differential methylation by current smoking at 4-71 CpGs may be modified by sex or dietary intake. Nearly half (35-50%) of differentially methylated CpGs on the 450 K array were associated with blood gene expression. Current smoking and in utero smoking CpGs implicated 3049 and 1067 druggable targets, including chemotherapy drugs.

INTERPRETATION

Many smoking-related methylation sites were identified with Illumina's EPIC array. Most signals revert to levels observed in never smokers within a year of cessation. Many in utero smoking CpGs persist into adulthood. Smoking-related druggable targets may provide insights into cancer treatment response and shared mechanisms across smoking-related diseases.

FUNDING

Intramural Research Program of the National Institutes of Health, Norwegian Ministry of Health and Care Services and the Ministry of Education and Research, Chief Scientist Office of the Scottish Government Health Directorates and the Scottish Funding Council, Medical Research Council UK and the Wellcome Trust.

Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects

Bulk-tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, and context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from the blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell-type proportions, we demonstrate that cell-type iQTLs could be considered as proxies for cell-type-specific QTL effects, particularly for the most abundant cell type in the tissue. The interpretation of age iQTLs, however, warrants caution because the moderation effect of age on the genotype and molecular phenotype association could be mediated by changes in cell-type composition. Finally, we show that cell-type iQTLs contribute to cell-type-specific enrichment of diseases that, in combination with additional functional data, could guide future functional studies. Overall, this study highlights the use of iQTLs to gain insights into the context specificity of regulatory effects.

Matched analysis of detailed peripheral blood and tumor immune microenvironment profiles in bladder cancer

Background: Bladder cancer and therapy responses hinge on immune profiles in the tumor microenvironment (TME) and blood, yet studies linking tumor-infiltrating immune cells to peripheral immune profiles are limited. Methods: DNA methylation cytometry quantified TME and matched peripheral blood immune cell proportions. With tumor immune profile data as the input, subjects were grouped by immune infiltration status and consensus clustering. Results: Immune hot and cold groups had different immune compositions in the TME but not in circulating blood. Two clusters of patients identified with consensus clustering had different immune compositions not only in the TME but also in blood. Conclusion: Detailed immune profiling via methylation cytometry reveals the significance of understanding tumor and systemic immune relationships in cancer patients.

DNA methylation at birth and lateral ventricular volume in childhood: a neuroimaging epigenetics study

BACKGROUND

Lateral ventricular volume (LVV) enlargement has been repeatedly linked to schizophrenia; yet, what biological factors shape LVV during early development remain unclear. DNA methylation (DNAm), an essential process for neurodevelopment that is altered in schizophrenia, is a key molecular system of interest.

METHODS

In this study, we conducted the first epigenome-wide association study of neonatal DNAm in cord blood with LVV in childhood (measured using T1-weighted brain scans at 10 years), based on data from a large population-based birth cohort, the Generation R Study (N = 840). Employing both probe-level and methylation profile score (MPS) approaches, we further examined whether epigenetic modifications identified at birth in cord blood are: (a) also observed cross-sectionally in childhood using peripheral blood DNAm at age of 10 years (Generation R, N = 370) and (b) prospectively associated with LVV measured in young adulthood in an all-male sample from the Avon Longitudinal Study of Parents and Children (ALSPAC, N = 114).

RESULTS

At birth, DNAm levels at four CpGs (annotated to potassium channel tetramerization domain containing 3, KCTD3; SHH signaling and ciliogenesis regulator, SDCCAG8; glutaredoxin, GLRX) prospectively associated with childhood LVV after genome-wide correction; these genes have been implicated in brain development and psychiatric traits including schizophrenia. An MPS capturing a broader epigenetic profile of LVV - but not individual top hits - showed significant cross-sectional associations with LVV in childhood in Generation R and prospectively associated with LVV in early adulthood within ALSPAC.

CONCLUSIONS

This study finds suggestive evidence that DNAm at birth prospectively associates with LVV at different life stages, albeit with small effect sizes. The prediction of MPS on LVV in a childhood sample and an independent male adult sample further underscores the stability and reproducibility of DNAm as a potential marker for LVV. Future studies with larger samples and comparable time points across development are needed to further elucidate how DNAm associates with this clinically relevant brain structure and risk for neuropsychiatric disorders, and what factors explain the identified DNAm profile of LVV at birth.

Prenatal BRCA1 epimutations contribute significantly to triple-negative breast cancer development

BACKGROUND

Normal cell BRCA1 epimutations have been associated with increased risk of triple-negative breast cancer (TNBC). However, the fraction of TNBCs that may have BRCA1 epimutations as their underlying cause is unknown. Neither are the time of occurrence and the potential inheritance patterns of BRCA1 epimutations established.

METHODS

To address these questions, we analyzed BRCA1 methylation status in breast cancer tissue and matched white blood cells (WBC) from 408 patients with 411 primary breast cancers, including 66 TNBCs, applying a highly sensitive sequencing assay, allowing allele-resolved methylation assessment. Furthermore, to assess the time of origin and the characteristics of normal cell BRCA1 methylation, we analyzed umbilical cord blood of 1260 newborn girls and 200 newborn boys. Finally, we assessed BRCA1 methylation status among 575 mothers and 531 fathers of girls with (n = 102) and without (n = 473) BRCA1 methylation.

RESULTS

We found concordant tumor and mosaic WBC BRCA1 epimutations in 10 out of 66 patients with TNBC and in four out of six patients with estrogen receptor (ER)-low expression (< 10%) tumors (combined: 14 out of 72; 19.4%; 95% CI 11.1-30.5). In contrast, we found concordant WBC and tumor methylation in only three out of 220 patients with 221 ER ≥ 10% tumors and zero out of 114 patients with 116 HER2-positive tumors. Intraindividually, BRCA1 epimutations affected the same allele in normal and tumor cells. Assessing BRCA1 methylation in umbilical WBCs from girls, we found mosaic, predominantly monoallelic BRCA1 epimutations, with qualitative features similar to those in adults, in 113/1260 (9.0%) of individuals, but no correlation to BRCA1 methylation status either in mothers or fathers. A significantly lower fraction of newborn boys carried BRCA1 methylation (9/200; 4.5%) as compared to girls (p = 0.038). Similarly, WBC BRCA1 methylation was found less common among fathers (16/531; 3.0%), as compared to mothers (46/575; 8.0%; p = 0.0003).

CONCLUSIONS

Our findings suggest prenatal BRCA1 epimutations might be the underlying cause of around 20% of TNBC and low-ER expression breast cancers. Such constitutional mosaic BRCA1 methylation likely arise through gender-related mechanisms in utero, independent of Mendelian inheritance.

Identification of intergenerational epigenetic inheritance by whole genome DNA methylation analysis in trios

Genome-wide association studies have identified thousands of loci associated with common diseases and traits. However, a large fraction of heritability remains unexplained. Epigenetic modifications, such as the observed in DNA methylation have been proposed as a mechanism of intergenerational inheritance. To investigate the potential contribution of DNA methylation to the missing heritability, we analysed the methylomes of four healthy trios (two parents and one offspring) using whole genome bisulphite sequencing. Of the 1.5 million CpGs (19%) with over 20% variability between parents in at least one family and compatible with a Mendelian inheritance pattern, only 3488 CpGs (0.2%) lacked correlation with any SNP in the genome, marking them as potential sites for intergenerational epigenetic inheritance. These markers were distributed genome-wide, with some preference to be located in promoters. They displayed a bimodal distribution, being either fully methylated or unmethylated, and were often found at the boundaries of genomic regions with high/low GC content. This analysis provides a starting point for future investigations into the missing heritability of simple and complex traits.

Uncertainty quantification of reference-based cellular deconvolution algorithms

The majority of epigenetic epidemiology studies to date have generated genome-wide profiles from bulk tissues (e.g., whole blood) however these are vulnerable to confounding from variation in cellular composition. Proxies for cellular composition can be mathematically derived from the bulk tissue profiles using a deconvolution algorithm; however, there is no method to assess the validity of these estimates for a dataset where the true cellular proportions are unknown. In this study, we describe, validate and characterize a sample level accuracy metric for derived cellular heterogeneity variables. The CETYGO score captures the deviation between a sample's DNA methylation profile and its expected profile given the estimated cellular proportions and cell type reference profiles. We demonstrate that the CETYGO score consistently distinguishes inaccurate and incomplete deconvolutions when applied to reconstructed whole blood profiles. By applying our novel metric to >6,300 empirical whole blood profiles, we find that estimating accurate cellular composition is influenced by both technical and biological variation. In particular, we show that when using a common reference panel for whole blood, less accurate estimates are generated for females, neonates, older individuals and smokers. Our results highlight the utility of a metric to assess the accuracy of cellular deconvolution, and describe how it can enhance studies of DNA methylation that are reliant on statistical proxies for cellular heterogeneity. To facilitate incorporating our methodology into existing pipelines, we have made it freely available as an R package (https://github.com/ds420/CETYGO).

DNA methylation differences within INS, PTPN22 and IL2RA promoters in lymphocyte subsets in children with type 1 diabetes and controls

We elucidated the effect of four known T1D-susceptibility associated single nucleotide polymorphism (SNP) markers in three genes (rs12722495 and rs2104286 in IL2RA, rs689 in INS and rs2476601 in PTPN22) on CpG site methylation of their proximal promoters in different lymphocyte subsets using pyrosequencing. The study cohort comprised 25 children with newly diagnosed T1D and 25 matched healthy controls. The rs689 SNP was associated with methylation at four CpG sites in INS promoter: -234, -206, -102 and -69. At all four CpG sites, the susceptibility genotype AA was associated with a higher methylation level compared to the other genotypes. We also found an association between rs12722495 and methylation at CpG sites -373 and -356 in IL2RA promoter in B cells, where the risk genotype AA was associated with lower methylation level compared to the AG genotype. The other SNPs analyzed did not demonstrate significant associations with CpG site methylation in the examined genes. Additionally, we compared the methylation between children with T1D and controls, and found statistically significant methylation differences at CpG -135 in INS in CD8+ T cells (p = 0.034), where T1D patients had a slightly higher methylation compared to controls (87.3 ± 7.2 vs. 78.8 ± 8.9). At the other CpG sites analyzed, the methylation was similar. Our results not only confirm the association between INS methylation and rs689 discovered in earlier studies but also report this association in sorted immune cells. We also report an association between rs12722495 and IL2RA promoter methylation in B cells. These results suggest that at least part of the genetic effect of rs689 and rs12722495 on T1D pathogenesis may be conveyed by DNA methylation.

Association of DNA methylation signatures with premature ageing and cardiovascular death in patients with end-stage kidney disease: a pilot epigenome-wide association study

Patients with end-stage kidney disease (ESKD) display features of premature aging. There is strong evidence that changes in DNA methylation (DNAm) contribute to age-related pathologies; however, little is known about their association with premature aging and cardiovascular mortality in patients with ESKD. We assayed genome-wide DNAm in a pilot case-control study of 60 hemodialysis patients with (n=30, cases) and without (n=30, controls) a fatal cardiovascular event. DNAm was profiled on the Illumina EPIC BeadChip. Four established DNAm clocks (i.e., Horvath-, Hannum-, Pheno-, and GrimAge) were used to estimate epigenetic age (DNAmAge). Epigenetic age acceleration (EAA) was derived as the residuals of regressing DNAmAge on chronological age (chroAge), and its association with cardiovascular death was examined using multivariable conditional logistic regression. An epigenome-wide association study (EWAS) was performed to identify differentially methylated CpGs associated with cardiovascular death. All clocks performed well at predicting chroAge (correlation between DNAmAges and chroAge of r=0.76-0.89), with GrimAge showing the largest deviation from chroAge (a mean of +21.3 years). There was no significant association of EAAs with cardiovascular death. In the EWAS, a CpG (cg22305782) in the FBXL19 gene had the strongest association with cardiovascular death with significantly lower DNAm in cases vs. controls (P_FDR=2.0x10^-6). FBXL19 is involved in cell apoptosis, inflammation, and adipogenesis. Overall, we observed more accelerated aging in patients with ESKD, although there was no significant association of EAAs with cardiovascular death. EWAS suggests a potential novel DNAm biomarker for premature cardiovascular mortality in ESKD.

A novel approach to risk exposure and epigenetics-the use of multidimensional context to gain insights into the early origins of cardiometabolic and neurocognitive health

BACKGROUND

Each mother-child dyad represents a unique combination of genetic and environmental factors. This constellation of variables impacts the expression of countless genes. Numerous studies have uncovered changes in DNA methylation (DNAm), a form of epigenetic regulation, in offspring related to maternal risk factors. How these changes work together to link maternal-child risks to childhood cardiometabolic and neurocognitive traits remains unknown. This question is a key research priority as such traits predispose to future non-communicable diseases (NCDs). We propose viewing risk and the genome through a multidimensional lens to identify common DNAm patterns shared among diverse risk profiles.

METHODS

We identified multifactorial Maternal Risk Profiles (MRPs) generated from population-based data (n = 15,454, Avon Longitudinal Study of Parents and Children (ALSPAC)). Using cord blood HumanMethylation450 BeadChip data, we identified genome-wide patterns of DNAm that co-vary with these MRPs. We tested the prospective relation of these DNAm patterns (n = 914) to future outcomes using decision tree analysis. We then tested the reproducibility of these patterns in (1) DNAm data at age 7 and 17 years within the same cohort (n = 973 and 974, respectively) and (2) cord DNAm in an independent cohort, the Generation R Study (n = 686).

RESULTS

We identified twenty MRP-related DNAm patterns at birth in ALSPAC. Four were prospectively related to cardiometabolic and/or neurocognitive childhood outcomes. These patterns were replicated in DNAm data from blood collected at later ages. Three of these patterns were externally validated in cord DNAm data in Generation R. Compared to previous literature, DNAm patterns exhibited novel spatial distribution across the genome that intersects with chromatin functional and tissue-specific signatures.

CONCLUSIONS

To our knowledge, we are the first to leverage multifactorial population-wide data to detect patterns of variability in DNAm. This context-based approach decreases biases stemming from overreliance on specific samples or variables. We discovered molecular patterns demonstrating prospective and replicable relations to complex traits. Moreover, results suggest that patterns harbour a genome-wide organisation specific to chromatin regulation and target tissues. These preliminary findings warrant further investigation to better reflect the reality of human context in molecular studies of NCDs.

Increased expression of individual genes in whole blood is associated with late-stage lung cancer at and close to diagnosis

Lung cancer (LC) mortality rates are still increasing globally. As survival is linked to stage, there is a need to identify markers for earlier LC diagnosis and individualized treatment. The whole blood transcriptome of LC patients represents a source of potential LC biomarkers. We compared expression of > 60,000 genes in whole blood specimens taken from LC cases at diagnosis (n = 128) and controls (n = 62) using genome-wide RNA sequencing, and identified 14 candidate genes associated with LC. High expression of ANXA3, ARG1 and HP was strongly associated with lower survival in late-stage LC cases (hazard ratios (HRs) = 2.81, 2.16 and 2.54, respectively). We validated these markers in two independent population-based studies with pre-diagnostic whole blood specimens taken up to eight years prior to LC diagnosis (n = 163 cases, 184 matched controls). ANXA3 and ARG1 expression was strongly associated with LC in these specimens, especially with late-stage LC within two years of diagnosis (odds ratios (ORs) = 3.47 and 5.00, respectively). Additionally, blood CD4 T cells, NK cells and neutrophils were associated with LC at diagnosis and improved LC discriminative ability beyond candidate genes. Our results indicate that in whole blood, increased expression levels of ANXA3, ARG1 and HP are diagnostic and prognostic markers of late-stage LC.

Alterations to DNA methylation patterns induced by chemotherapy treatment are associated with negative impacts on the olfactory pathway

BACKGROUND

Exposure to cytotoxic chemotherapy treatment may alter DNA methylation (DNAm) in breast cancer patients.

METHODS

We performed DNAm analysis in 125 breast cancer patients with blood drawn before and after chemotherapy, using the Illumina MethylationEPIC array. DNAm changes of 588,798 individual CpGs (including 41,207 promoter regions) were evaluated using linear regression models adjusted for monocyte proportion. Gene set enrichment analyses (GSEA) were conducted to identify key Gene Ontology (GO) biological processes or Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways associated with chemotherapy. Results were validated in a separate cohort of breast cancer patients who were treated (n = 1273) and not treated (n = 872) by chemotherapy (1808 blood, 337 saliva).

RESULTS

A total of 141 differentially methylated CpGs and 11 promoters were significantly associated with chemotherapy after multiple testing corrections in both the paired sample and single time point analyses. GSEA of promoter regions (pre-ranked by test statistics) identified six suppressed biological processes (p < 4.67e-8) related to sensory perception and detection of chemical stimuli, including smell perception (GO:0007606, GO:0007608, GO:0009593, GO:0050906, GO:0050907, and GO:0050911). The same six biological processes were significantly suppressed in the validation dataset (p < 9.02e-14). The KEGG pathway olfactory transduction (hsa04740) was also found to be significantly suppressed (ppaired-samples = 1.72e-9, psingle-timepoint-blood = 2.03e-15 and psingle-timepoint-saliva = 7.52e-56).

CONCLUSION

The enrichment of imprinted genes within biological processes and pathways suggests a biological mechanism by which chemotherapy could affect the perception of smell.

Targeted DNA Methylation Analysis Facilitates Leukocyte Counts in Dried Blood Samples

BACKGROUND

Cell-type specific DNA methylation (DNAm) can be employed to determine the numbers of leukocyte subsets in blood. In contrast to conventional methods for leukocyte counts, which are based on cellular morphology or surface marker protein expression, the cellular deconvolution based on DNAm levels is applicable for frozen or dried blood. Here, we further enhanced targeted DNAm assays for leukocyte counts in clinical application.

METHODS

DNAm profiles of 40 different studies were compiled to identify CG dinucleotides (CpGs) with cell-type specific DNAm using a computational framework, CimpleG. DNAm levels at these CpGs were then measured with digital droplet PCR in venous blood from 160 healthy donors and 150 patients with various hematological disorders. Deconvolution was further validated with venous blood (n = 75) and capillary blood (n = 31) that was dried on Whatman paper or on Mitra microsampling devices.

RESULTS

In venous blood, automated cell counting or flow cytometry correlated well with epigenetic estimates of relative leukocyte counts for granulocytes (r = 0.95), lymphocytes (r = 0.97), monocytes (r = 0.82), CD4 T cells (r = 0.84), CD8 T cells (r = 0.94), B cells (r = 0.96), and NK cells (r = 0.72). Similar correlations and precisions were achieved for dried blood samples. Spike-in with a reference plasmid enabled accurate epigenetic estimation of absolute leukocyte counts from dried blood samples, correlating with conventional venous (r = 0.86) and capillary (r = 0.80) blood measurements.

CONCLUSIONS

The advanced selection of cell-type specific CpGs and utilization of digital droplet PCR analysis provided accurate epigenetic blood counts. Analysis of dried blood facilitates self-sampling with a finger prick, thereby enabling easier accessibility to testing.

Host fecal DNA specific methylation signatures mark gut dysbiosis and inflammation in children affected by autism spectrum disorder

The gut-brain axis involves several bidirectional pathway communications including microbiome, bacterial metabolites, neurotransmitters as well as immune system and is perturbed both in brain and in gastrointestinal disorders. Consistently, microbiota-gut-brain axis has been found altered in autism spectrum disorder (ASD). We reasoned that such alterations occurring in ASD may impact both on methylation signatures of human host fecal DNA (HFD) and possibly on the types of human cells shed in the stools from intestinal tract giving origin to HFD. To test this hypothesis, we have performed whole genome methylation analysis of HFD from an age-restricted cohort of young children with ASD (N = 8) and healthy controls (N = 7). In the same cohort we have previously investigated the fecal microbiota composition and here we refined such analysis and searched for eventual associations with data derived from HFD methylome analysis. Our results showed that specific epigenetic signatures in human fecal DNA, especially at genes related to inflammation, associated with the disease. By applying methylation-based deconvolution algorithm, we found that the HFD derived mainly from immune cells and the relative abundance of those differed between patients and controls. Consistently, most of differentially methylated regions fitted with genes involved in inflammatory response. Interestingly, using Horvath epigenetic clock, we found that ASD affected children showed both epigenetic and microbiota age accelerated. We believe that the present unprecedented approach may be useful for the identification of the ASD associated HFD epigenetic signatures and may be potentially extended to other brain disorders and intestinal inflammatory diseases.

The Origin of Highly Elevated Cell-Free DNA in Healthy Individuals and Patients with Pancreatic, Colorectal, Lung, or Ovarian Cancer

Cell-free DNA (cfDNA) concentrations from patients with cancer are often elevated compared with those of healthy controls, but the sources of this extra cfDNA have never been determined. To address this issue, we assessed cfDNA methylation patterns in 178 patients with cancers of the colon, pancreas, lung, or ovary and 64 patients without cancer. Eighty-three of these individuals had cfDNA concentrations much greater than those generally observed in healthy subjects. The major contributor of cfDNA in all samples was leukocytes, accounting for ∼76% of cfDNA, with neutrophils predominating. This was true regardless of whether the samples were derived from patients with cancer or the total plasma cfDNA concentration. High levels of cfDNA observed in patients with cancer did not come from either neoplastic cells or surrounding normal epithelial cells from the tumor's tissue of origin. These data suggest that cancers may have a systemic effect on cell turnover or DNA clearance.

SIGNIFICANCE

The origin of excess cfDNA in patients with cancer is unknown. Using cfDNA methylation patterns, we determined that neither the tumor nor the surrounding normal tissue contributes this excess cfDNA-rather it comes from leukocytes. This finding suggests that cancers have a systemic impact on cell turnover or DNA clearance. See related commentary by Thierry and Pisareva, p. 2122. This article is featured in Selected Articles from This Issue, p. 2109.

Mouse methylation profiles for leukocyte cell types, and estimation of leukocyte fractions in inflamed gastrointestinal DNA samples

Precise analysis of tissue DNA and RNA samples is often hampered by contaminating non-target cells whose amounts are highly variable. DNA methylation profiles are specific to cell types, and can be utilized for assessment of the fraction of such contaminating non-target cells. Here, we aimed 1) to identify methylation profiles specific to multiple types of mouse leukocytes, and 2) to estimate the fraction of leukocytes infiltrating inflamed tissues using DNA samples. First, genome-wide DNA methylation analysis was conducted for three myeloid-lineage cells and four lymphoid-lineage cells isolated by fluorescence-activated cell sorting after magnetic-activated cell sorting from leukocytes in the spleen. Clustering analysis using CpG sites within enhancers separated the three myeloid-lineage cells and four lymphoid-lineage cells while that using promoter CpG islands (TSS200CGIs) did not. Among the 266,108 CpG sites analyzed, one CpG site was specifically hypermethylated (β value ≥ 0.7) in B cells, and four, seven, 183, and 34 CpG sites were specifically hypomethylated (β value < 0.2) in CD4+ T cells, CD8+ T cells, B cells, and NK cells, respectively. Importantly, cell type-specific hypomethylated CpG sites were located at genes involved in cell type-specific biological functions. Then, marker CpG sites to estimate the leukocyte fraction in a tissue with leukocyte infiltration were selected, and an estimation algorithm was established. The fractions of infiltrating leukocytes were estimated to be 1.6-12.4% in the stomach (n = 10) with Helicobacter pylori-induced inflammation and 1.5-4.3% in the colon with dextran sulfate sodium-induced colitis (n = 4), and the fractions were highly correlated with those estimated histologically using Cd45-stained tissue sections [R = 0.811 (p = 0.004)]. These results showed that mouse methylation profiles at CpG sites within enhancers reflected leukocyte cell lineages, and the use of marker CpG sites successfully estimated the leukocyte fraction in inflamed gastric and colon tissues.

Genome-Scale Methylation Analysis Identifies Immune Profiles and Age Acceleration Associations with Bladder Cancer Outcomes

BACKGROUND

Immune profiles have been associated with bladder cancer outcomes and may have clinical applications for prognosis. However, associations of detailed immune cell subtypes with patient outcomes remain underexplored and may contribute crucial prognostic information for better managing bladder cancer recurrence and survival.

METHODS

Bladder cancer case peripheral blood DNA methylation was measured using the Illumina HumanMethylationEPIC array. Extended cell-type deconvolution quantified 12 immune cell-type proportions, including memory, naïve T and B cells, and granulocyte subtypes. DNA methylation clocks determined biological age. Cox proportional hazards models tested associations of immune cell profiles and age acceleration with bladder cancer outcomes. The partDSA algorithm discriminated 10-year overall survival groups from clinical variables and immune cell profiles, and a semi-supervised recursively partitioned mixture model (SS-RPMM) with DNA methylation data was applied to identify a classifier for 10-year overall survival.

RESULTS

Higher CD8T memory cell proportions were associated with better overall survival [HR = 0.95, 95% confidence interval (CI) = 0.93-0.98], while higher neutrophil-to-lymphocyte ratio (HR = 1.36, 95% CI = 1.23-1.50), CD8T naïve (HR = 1.21, 95% CI = 1.04-1.41), neutrophil (HR = 1.04, 95% CI = 1.03-1.06) proportions, and age acceleration (HR = 1.06, 95% CI = 1.03-1.08) were associated with worse overall survival in patient with bladder cancer. partDSA and SS-RPMM classified five groups of subjects with significant differences in overall survival.

CONCLUSIONS

We identified associations between immune cell subtypes and age acceleration with bladder cancer outcomes.

IMPACT

The findings of this study suggest that bladder cancer outcomes are associated with specific methylation-derived immune cell-type proportions and age acceleration, and these factors could be potential prognostic biomarkers.

Genome-wide analysis of DNA methylation and its relationship with serum homocysteine levels in patients with hypertension

BACKGROUND

Homocysteine (Hcy) is an independent risk factor for cardiovascular diseases, and elevated plasma Hcy levels could aggravate vascular injury in hypertension. Hyperhomocysteinemia can change the methylation status of global DNA and specific genes. In the present study, we aim to examine the comprehensive influence of Hcy levels on DNA methylation status in patients with hypertension.

METHODS

Epigenome-wide methylation profiles of the peripheral leukocyte DNA of 218 patients with hypertension were analyzed using the Illumina Infinium Methylation EPIC BeadChip. Differentially methylated positions (DMPs) associated with serum Hcy levels were identified by mixed linear regression with the adjustment of potential confounders. Gene Ontology analysis and Kyoto Encyclopedia of Genes and Genomes pathway analysis were conducted to determine the potential functions of the identified DMPs. The association between the methylation level of DMPs and carotid-femoral pulse wave velocity (Cf-PWV) was also analyzed.

RESULTS

Five DMPs at cg13169662, cg03179312, cg21976560, cg25262698, and cg09433843 showed significant association with serum Hcy levels (false discovery rate-corrected P  < 0.05). An additional six CpG sites met the threshold for suggestive significance ( P  < 1 × 10 -6 ), among which three DMPs (cg25781123, cg26463106, and cg06679221) were annotated to THUMPD3 . Furthermore, the methylation levels of cg13169662 and cg25262698 (RPRD1A) were significantly associated with Cf-PWV.

CONCLUSION

Our results suggest that Hcy could induce DNA methylation alteration in patients with hypertension. Further functional research is warranted to elucidate the concrete role of DMPs in hypertension.

Genome-wide DNA methylation landscape of four Chinese populations and epigenetic variation linked to Tibetan high-altitude adaptation

DNA methylation (DNAm) is one of the major epigenetic mechanisms in humans and is important in diverse cellular processes. The variation of DNAm in the human population is related to both genetic and environmental factors. However, the DNAm profiles have not been investigated in the Chinese population of diverse ethnicities. Here, we performed double-strand bisulfite sequencing (DSBS) for 32 Chinese individuals representing four major ethnic groups including Han Chinese, Tibetan, Zhuang, and Mongolian. We identified a total of 604,649 SNPs and quantified DNAm at more than 14 million CpGs in the population. We found global DNAm-based epigenetic structure is different from the genetic structure of the population, and ethnic difference only partially explains the variation of DNAm. Surprisingly, non-ethnic-specific DNAm variations showed stronger correlation with the global genetic divergence than these ethnic-specific DNAm. Differentially methylated regions (DMRs) among these ethnic groups were found around genes in diverse biological processes. Especially, these DMR-genes between Tibetan and non-Tibetans were enriched around high-altitude genes including EPAS1 and EGLN1, suggesting DNAm alteration plays an important role in high-altitude adaptation. Our results provide the first batch of epigenetic maps for Chinese populations and the first evidence of the association of epigenetic changes with Tibetans' high-altitude adaptation.

Analysis of DNA methylation at birth and in childhood reveals changes associated with season of birth and latitude

BACKGROUND

Seasonal variations in environmental exposures at birth or during gestation are associated with numerous adult traits and health outcomes later in life. Whether DNA methylation (DNAm) plays a role in the molecular mechanisms underlying the associations between birth season and lifelong phenotypes remains unclear.

METHODS

We carried out epigenome-wide meta-analyses within the Pregnancy And Childhood Epigenetic Consortium to identify associations of DNAm with birth season, both at differentially methylated probes (DMPs) and regions (DMRs). Associations were examined at two time points: at birth (21 cohorts, N = 9358) and in children aged 1-11 years (12 cohorts, N = 3610). We conducted meta-analyses to assess the impact of latitude on birth season-specific associations at both time points.

RESULTS

We identified associations between birth season and DNAm (False Discovery Rate-adjusted p values < 0.05) at two CpGs at birth (winter-born) and four in the childhood (summer-born) analyses when compared to children born in autumn. Furthermore, we identified twenty-six differentially methylated regions (DMR) at birth (winter-born: 8, spring-born: 15, summer-born: 3) and thirty-two in childhood (winter-born: 12, spring and summer: 10 each) meta-analyses with few overlapping DMRs between the birth seasons or the two time points. The DMRs were associated with genes of known functions in tumorigenesis, psychiatric/neurological disorders, inflammation, or immunity, amongst others. Latitude-stratified meta-analyses [higher (≥ 50°N), lower (< 50°N, northern hemisphere only)] revealed differences in associations between birth season and DNAm by birth latitude. DMR analysis implicated genes with previously reported links to schizophrenia (LAX1), skin disorders (PSORS1C, LTB4R), and airway inflammation including asthma (LTB4R), present only at birth in the higher latitudes (≥ 50°N).

CONCLUSIONS

In this large epigenome-wide meta-analysis study, we provide evidence for (i) associations between DNAm and season of birth that are unique for the seasons of the year (temporal effect) and (ii) latitude-dependent variations in the seasonal associations (spatial effect). DNAm could play a role in the molecular mechanisms underlying the effect of birth season on adult health outcomes.

A meta-analysis of immune-cell fractions at high resolution reveals novel associations with common phenotypes and health outcomes

BACKGROUND

Changes in cell-type composition of tissues are associated with a wide range of diseases and environmental risk factors and may be causally implicated in disease development and progression. However, these shifts in cell-type fractions are often of a low magnitude, or involve similar cell subtypes, making their reliable identification challenging. DNA methylation profiling in a tissue like blood is a promising approach to discover shifts in cell-type abundance, yet studies have only been performed at a relatively low cellular resolution and in isolation, limiting their power to detect shifts in tissue composition.

METHODS

Here we derive a DNA methylation reference matrix for 12 immune-cell types in human blood and extensively validate it with flow-cytometric count data and in whole-genome bisulfite sequencing data of sorted cells. Using this reference matrix, we perform a directional Stouffer and fixed effects meta-analysis comprising 23,053 blood samples from 22 different cohorts, to comprehensively map associations between the 12 immune-cell fractions and common phenotypes. In a separate cohort of 4386 blood samples, we assess associations between immune-cell fractions and health outcomes.

RESULTS

Our meta-analysis reveals many associations of cell-type fractions with age, sex, smoking and obesity, many of which we validate with single-cell RNA sequencing. We discover that naïve and regulatory T-cell subsets are higher in women compared to men, while the reverse is true for monocyte, natural killer, basophil, and eosinophil fractions. Decreased natural killer counts associated with smoking, obesity, and stress levels, while an increased count correlates with exercise and sleep. Analysis of health outcomes revealed that increased naïve CD4 + T-cell and N-cell fractions associated with a reduced risk of all-cause mortality independently of all major epidemiological risk factors and baseline co-morbidity. A machine learning predictor built only with immune-cell fractions achieved a C-index value for all-cause mortality of 0.69 (95%CI 0.67-0.72), which increased to 0.83 (0.80-0.86) upon inclusion of epidemiological risk factors and baseline co-morbidity.

CONCLUSIONS

This work contributes an extensively validated high-resolution DNAm reference matrix for blood, which is made freely available, and uses it to generate a comprehensive map of associations between immune-cell fractions and common phenotypes, including health outcomes.