The correlation of methylation levels measured using Illumina 450K and EPIC BeadChips in blood samples

Critical evaluation of the reliability of DNA methylation probes on the Illumina MethylationEPIC v1.0 BeadChip microarrays

DNA methylation (DNAm) plays a crucial role in a number of complex diseases. However, the reliability of DNAm levels measured using Illumina arrays varies across different probes. Previous research primarily assessed probe reliability by comparing duplicate samples between the 450k-450k or 450k-EPIC platforms, with limited investigations on Illumina EPIC v1.0 arrays. We conducted a comprehensive assessment of the EPIC v1.0 array probe reliability using 69 blood DNA samples, each measured twice, generated by the Alzheimer's Disease Neuroimaging Initiative study. We observed higher reliability in probes with average methylation beta values of 0.2 to 0.8, and lower reliability in type I probes or those within the promoter and CpG island regions. Importantly, we found that probe reliability has significant implications in the analyses of Epigenome-wide Association Studies (EWAS). Higher reliability is associated with more consistent effect sizes in different studies, the identification of differentially methylated regions (DMRs) and methylation quantitative trait locus (mQTLs), and significant correlations with downstream gene expression. Moreover, blood DNAm measurements obtained from probes with higher reliability are more likely to show concordance with brain DNAm measurements. Our findings, which provide crucial reliability information for probes on the EPIC v1.0 array, will serve as a valuable resource for future DNAm studies.

Associations between epigenetic aging and diabetes mellitus in a Swedish longitudinal study

Diabetes mellitus type 2 (T2D) is associated with accelerated biological aging and the increased risk of onset of other age-related diseases. Epigenetic changes in DNA methylation levels have been found to serve as reliable biomarkers for biological aging. This study explores the relationship between various epigenetic biomarkers of aging and diabetes risk using longitudinal data. Data from the Swedish Adoption/Twin Study of Aging (SATSA) was collected from 1984 to 2014 and included 536 individuals with at least one epigenetic measurement. The following epigenetic biomarkers of aging were employed: DNAm PAI-1, DNAmTL, DunedinPACE, PCHorvath1, PCHorvath2, PCHannum, PCPhenoAge, and PCGrimAge. Firstly, longitudinal analysis of biomarker trajectories was done. Secondly, linear correlations between the biomarkers and time to diabetes were studied within individuals developing diabetes. Thirdly, Cox proportional hazards (PH) models were used to assess the associations between these biomarkers and time of diabetes diagnosis, with adjustments for chronological age, sex, education, smoking, blood glucose, and BMI. The longitudinal trajectories of the biomarkers revealed differences between individuals with and without diabetes. Smoothened average curves for DunedinPACE and DNAm PAI-1 were higher for individuals with diabetes around the age 60-70, compared to controls. Likewise, DunedinPACE and DNAm PAI-1 were higher closer to diabetes onset. However, no significant associations were found between the epigenetic biomarkers of aging and risk of diabetes in Cox PH models. Our findings suggest the potential value of developing epigenetic biomarkers specifically tailored to T2D, should we wish to model and explore the potential for predicting the disease.

Placental PFAS concentrations are associated with perturbations of placental DNA methylation at loci with important roles on cardiometabolic health

The placenta is crucial for fetal development, is affected by PFAS toxicity, and evidence is accumulating that gestational PFAS perturb the epigenetic activity of the placenta. Gestational PFAS exposure is can adversely affect offspring, yet individual and cumulative impacts of PFAS on the placental epigenome remain underexplored. Here, we conducted an epigenome-wide association study (EWAS) to examine the relationships between placental PFAS levels and DNA methylation in a cohort of mother-infant dyads in Arkansas. We measured 17 PFAS in human placental tissues and quantified placental DNA methylation levels via the Illumina EPIC Microarray. We tested for differential DNA methylation with individual PFAS, and with mixtures of multiple PFAS. Our results demonstrated that numerous epigenetic loci were perturbed by PFAS, with PFHxS exhibiting the most abundant effects. Mixture analyses suggested cumulative effects of PFOA and PFOS, while PFHxS may act more independently. We additionally explored whether sex-specific effects may be present and concluded that future large studies should explicitly test for sex-specific effects. The genes that are annotated to our PFAS-associated epigenetic loci are primarily involved in growth processes and cardiometabolic health, while some genes are involved in neurodevelopment. These findings shed light on how prenatal PFAS exposures affect birth outcomes and children's health, emphasizing the importance of understanding PFAS mechanisms in the in-utero environment.

Placental DNA methylation signatures of prenatal air pollution exposure and potential effects on birth outcomes: an analysis of three prospective cohorts

BACKGROUND

Pregnancy air pollution exposure (PAPE) has been linked to a wide range of adverse birth and childhood outcomes, but there is a paucity of data on its influence on the placental epigenome, which can regulate the programming of physiological functions and affect child development. This study aimed to investigate the association between prenatal air pollutant exposure concentrations and changes in placental DNA methylation patterns, and to explore the potential windows of susceptibility and sex-specific alterations.

METHODS

This multi-site study used three prospective population-based mother-child cohorts: EDEN, PELAGIE, and SEPAGES, originating from four French geographical regions (Nancy, Poitiers, Brittany, and Grenoble). Pregnant women were included between 2003 and 2006 for EDEN and PELAGIE, and between 2014 and 2017 for SEPAGES. The main eligibility criteria were: being older than 18 years, having a singleton pregnancy, and living and planning to deliver in one of the maternity clinics in one of the study areas. A total of 1539 mother-child pairs were analysed, measuring placental DNA methylation using Illumina BeadChips. We used validated spatiotemporally resolved models to estimate PM2·5, PM10, and NO2 exposure over each trimester of pregnancy at the maternal residential address. We conducted a pooled adjusted epigenome-wide association study to identify differentially methylated 5'-C-phosphate-G-3' (CpG) sites and regions (assessed using the Infinium HumanMethylationEPIC BeadChip array, n=871), including sex-specific and sex-linked alterations, and independently validated our results (assessed using the Infinium HumanMethylation450 BeadChip array, n=668).

FINDINGS

We identified four CpGs and 28 regions associated with PAPE in the total population, 469 CpGs and 87 regions in male infants, and 150 CpGs and 66 regions in female infants. We validated 35% of the CpGs available. More than 30% of the identified CpGs were related to one (or more) birth outcome and most significant alterations were enriched for neural development, immunity, and metabolism related genes. The 28 regions identified for both sexes overlapped with imprinted genes (four genes), and were associated with neurodevelopment (nine genes), immune system (seven genes), and metabolism (five genes). Most associations were observed for the third trimester for female infants (134 of 150 CpGs), and throughout pregnancy (281 of 469 CpGs) and the first trimester (237 of 469 CpGs) for male infants.

INTERPRETATION

These findings highlight the molecular pathways through which PAPE might affect child health in a widespread and sex-specific manner, identifying the genes involved in the major physiological functions of a developing child. Further studies are needed to elucidate whether these epigenetic changes persist and affect health later in life.

FUNDING

French Agency for National Research, Fondation pour la Recherche Médicale, Fondation de France, and the Plan Cancer.

Toward a more holistic approach to the study of exposures and child outcomes

Aim: The current work was designed to demonstrate the application of the exposome framework in examining associations between exposures and children's long-term neurodevelopmental and behavioral outcomes. Methods: Longitudinal data were collected from birth through age 6 from 402 preterm infants. Three statistical methods were utilized to demonstrate the exposome framework: exposome-wide association study, cumulative exposure and machine learning models, with and without epigenetic data. Results: Each statistical approach answered a distinct research question regarding the impact of exposures on longitudinal child outcomes. Findings highlight associations between exposures, epigenetics and executive function. Conclusion: Findings demonstrate how an exposome-based approach can be utilized to understand relationships between internal (e.g., DNA methylation) and external (e.g., prenatal risk) exposures and long-term developmental outcomes in preterm children.

Epigenetic associations in HPA axis genes related to bronchopulmonary dysplasia and antenatal steroids

BACKGROUND

Bronchopulmonary dysplasia (BPD), a common morbidity among very preterm infants, is associated with chronic disease and neurodevelopmental impairments. A hypothesized mechanism for these outcomes lies in altered glucocorticoid (GC) activity. We hypothesized that BPD and its treatments may result in epigenetic differences in the hypothalamic-pituitary-adrenal (HPA) axis, which is modulated by GC, and could be ascertained using an established GC risk score and DNA methylation (DNAm) of HPA axis genes.

METHODS

DNAm was quantified from buccal tissue (ECHO-NOVI) and from neonatal blood spots (ELGAN ECHO) via the EPIC microarray. Prenatal maternal characteristics, pregnancy complication, and neonatal medical complication data were collected from medical record review and maternal interviews.

RESULTS

The GC score was not associated with steroid exposure or BPD. However, six HPA genes involved in stress response regulation demonstrated differential methylation with antenatal steroid exposure; two CpGs within FKBP5 and POMC were differentially methylated with BPD severity. These findings were sex-specific in both cohorts; males had greater magnitude of differential methylation within these genes.

CONCLUSIONS

These findings suggest that BPD severity and antenatal steroids are associated with DNAm at some HPA genes in very preterm infants and the effects appear to be sex-, tissue-, and age-specific.

IMPACT

This study addresses bronchopulmonary dysplasia (BPD), an important health outcome among preterm neonates, and interrogates a commonly studied pathway, the hypothalamic-pituitary-adrenal (HPA) axis. The combination of BPD, the HPA axis, and epigenetic markers has not been previously reported. In this study, we found that BPD itself was not associated with epigenetic responses in the HPA axis in infants born very preterm; however, antenatal treatment with steroids was associated with epigenetic responses.

Former smoking associated with epigenetic modifications in human granulosa cells among women undergoing assisted reproduction

Smoking exposure during adulthood can disrupt oocyte development in women, contributing to infertility and possibly adverse birth outcomes. Some of these effects may be reflected in epigenome profiles in granulosa cells (GCs) in human follicular fluid. We compared the epigenetic modifications throughout the genome in GCs from women who were former (N = 15) versus never smokers (N = 44) undergoing assisted reproductive technologies (ART). This study included 59 women undergoing ART. Smoking history including time since quitting was determined by questionnaire. GCs were collected during oocyte retrieval and DNA methylation (DNAm) levels were profiled using the Infinium MethylationEPIC BeadChip. We performed an epigenome-wide association study with robust linear models, regressing DNAm level at individual loci on smoking status, adjusting for age, ovarian stimulation protocol, and three surrogate variables. We performed differentially methylated regions (DMRs) analysis and over-representation analysis of the identified CpGs and corresponding gene set. 81 CpGs were differentially methylated among former smokers compared to never smokers (FDR < 0.05). We identified 2 significant DMRs (KCNQ1 and RHBDD2). The former smoking-associated genes were enriched in oxytocin signaling, adrenergic signaling in cardiomyocytes, platelet activation, axon guidance, and chemokine signaling pathway. These epigenetic variations have been associated with inflammatory responses, reproductive outcomes, cancer development, neurodevelopmental disorder, and cardiometabolic health. Secondarily, we examined the relationships between time since quitting and DNAm at significant CpGs. We observed three CpGs in negative associations with the length of quitting smoking (p < 0.05), which were cg04254052 (KCNIP1), cg22875371 (OGDHL), and cg27289628 (LOC148145), while one in positive association, which was cg13487862 (PLXNB1). As a pilot study, we demonstrated epigenetic modifications associated with former smoking in GCs. The study is informative to potential biological pathways underlying the documented association between smoking and female infertility and biomarker discovery for smoking-associated reproductive outcomes.

Epigenome-wide association study identifies neonatal DNA methylation associated with two-year attention problems in children born very preterm

Prior research has identified epigenetic predictors of attention problems in school-aged children but has not yet investigated these in young children, or children at elevated risk of attention problems due to preterm birth. The current study evaluated epigenome-wide associations between neonatal DNA methylation and attention problems at age 2 years in children born very preterm. Participants included 441 children from the Neonatal Neurobehavior and Outcomes in Very Preterm Infants (NOVI) Study, a multi-site study of infants born < 30 weeks gestational age. DNA methylation was measured from buccal swabs collected at NICU discharge using the Illumina MethylationEPIC Bead Array. Attention problems were assessed at 2 years of adjusted age using the attention problems subscale of the Child Behavior Checklist (CBCL). After adjustment for multiple testing, DNA methylation at 33 CpG sites was associated with child attention problems. Differentially methylated CpG sites were located in genes previously linked to physical and mental health, including several genes associated with ADHD in prior epigenome-wide and genome-wide association studies. Several CpG sites were located in genes previously linked to exposure to prenatal risk factors in the NOVI sample. Neonatal epigenetics measured at NICU discharge could be useful in identifying preterm children at risk for long-term attention problems and related psychiatric disorders, who could benefit from early prevention and intervention efforts.

Biological stability of DNA methylation measurements over varying intervals of time and in the presence of acute stress

Identifying factors that influence the stability of DNA methylation measurements across biological replicates is of critical importance in basic and clinical research. Using a within-person between-group experimental design (n = 31, number of observations = 192), we report the stability of biological replicates over a variety of unique temporal scenarios, both in the absence and presence of acute psychosocial stress, and between individuals who have experienced early life adversity (ELA) and non-exposed individuals. We found that varying time intervals, acute stress, and ELA exposure influenced the stability of repeated DNA methylation measurements. In the absence of acute stress, probes were less stable as time passed; however, stress exerted a stabilizing influence on probes over longer time intervals. Compared to non-exposed individuals, ELA-exposed individuals had significantly lower probe stability immediately following acute stress. Additionally, we found that across all scenarios, probes used in most epigenetic-based algorithms for estimating epigenetic age or immune cell proportions had average or below-average stability, except for the Principal Component and DunedinPACE epigenetic ageing clocks, which were enriched for more stable probes. Finally, using highly stable probes in the absence of stress, we identified multiple probes that were hypomethylated in the presence of acute stress, regardless of ELA status. Two hypomethylated probes are located near the transcription start site of the glutathione-disulfide reductase gene (GSR), which has previously been shown to be an integral part of the stress response to environmental toxins. We discuss implications for future studies concerning the reliability and reproducibility of DNA methylation measurements.Abbreviations: DNAm - DNA methylation, CpG - 5'-cytosine-phosphate-guanine-3,' ICC - Interclass correlation coefficient, ELA - Early-life adversity, PBMCs - Peripheral blood mononuclear cells, mQTL - Methylation quantitative trait loci, TSS - Transcription start site, GSR - Glutathione-disulfide reductase gene, TSST - Trier social stress test, PC - Principal component.

A consistent pattern of slide effects in Illumina DNA methylation BeadChip array data

Background: Recent studies have identified thousands of associations between DNA methylation CpGs and complex diseases/traits, emphasizing the critical role of epigenetics in understanding disease aetiology and identifying biomarkers. However, association analyses based on methylation array data are susceptible to batch/slide effects, which can lead to inflated false positive rates or reduced statistical powerResults: We use multiple DNA methylation datasets based on the popular Illumina Infinium MethylationEPIC BeadChip array to describe consistent patterns and the joint distribution of slide effects across CpGs, confirming and extending previous results. The susceptible CpGs overlap with the Illumina Infinium HumanMethylation450 BeadChip array content.Conclusions: Our findings reveal systematic patterns in slide effects. The observations provide further insights into the characteristics of these effects and can improve existing adjustment approaches.

Technical Report: A Comprehensive Comparison between Different Quantification Versions of Nightingale Health's 1H-NMR Metabolomics Platform

1H-NMR metabolomics data is increasingly used to track health and disease. Nightingale Health, a major supplier of 1H-NMR metabolomics, has recently updated the quantification strategy to further align with clinical standards. Such updates, however, might influence backward replicability, particularly affecting studies with repeated measures. Using data from BBMRI-NL consortium (~28,000 samples from 28 cohorts), we compared Nightingale data, originally released in 2014 and 2016, with a re-quantified version released in 2020, of which both versions were based on the same NMR spectra. Apart from two discontinued and twenty-three new analytes, we generally observe a high concordance between quantification versions with 73 out of 222 (33%) analytes showing a mean ρ > 0.9 across all cohorts. Conversely, five analytes consistently showed lower Spearman's correlations (ρ < 0.7) between versions, namely acetoacetate, LDL-L, saturated fatty acids, S-HDL-C, and sphingomyelins. Furthermore, previously trained multi-analyte scores, such as MetaboAge or MetaboHealth, might be particularly sensitive to platform changes. Whereas MetaboHealth replicated well, the MetaboAge score had to be retrained due to use of discontinued analytes. Notably, both scores in the re-quantified data recapitulated mortality associations observed previously. Concluding, we urge caution in utilizing different platform versions to avoid mixing analytes, having different units, or simply being discontinued.

DNA Methylation Age Acceleration Is Not Associated with Age of Onset in Parkinson's Disease

BACKGROUND

Epigenetic clocks using DNA methylation (DNAm) to estimate biological age have become popular tools in the study of neurodegenerative diseases. Notably, several recent reports have shown a strikingly similar inverse relationship between accelerated biological aging, as measured by DNAm, and the age of onset of several neurodegenerative disorders, including Parkinson's disease (PD). Common to all of these studies is that they were performed without control subjects and using the exact same measure of accelerated aging: DNAm age minus chronological age.

OBJECTIVE

We aimed to assess the validity of these findings in PD, using the same dataset as in the original study, blood DNAm data from the Parkinson's Progression Markers Initiative cohort, but also including control samples in the analyses.

METHODS

We replicated the analyses and findings of the previous study and then reanalyzed the dataset incorporating control samples to account for underlying age-related biases.

RESULTS

Our reanalysis shows that there is no correlation between age of onset and DNAm age acceleration. Conversely, there is a pattern of overestimating DNAm age in younger and underestimating DNAm age in older individuals in the dataset that entirely explains the previously reported association.

CONCLUSIONS

Our findings refute the previously reported inverse relationship between DNAm age acceleration and age of onset in PD. We show that these findings are fully accounted for by an expected over/underestimation of DNAm age in younger/older individuals. Furthermore, this effect is likely to be responsible for nearly identical findings reported in other neurodegenerative diseases. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Critical evaluation of the reliability of DNA methylation probes on the Illumina MethylationEPIC BeadChip microarrays

DNA methylation (DNAm) plays a crucial role in a number of complex diseases. However, the reliability of DNAm levels measured using Illumina arrays varies across different probes. Previous research primarily assessed probe reliability by comparing duplicate samples between the 450k-450k or 450k-EPIC platforms, with limited investigations on Illumina EPIC arrays. We conducted a comprehensive assessment of the EPIC array probe reliability using 138 duplicated blood DNAm samples generated by the Alzheimer's Disease Neuroimaging Initiative study. We introduced a novel statistical measure, the modified intraclass correlation, to better account for the disagreement in duplicate measurements. We observed higher reliability in probes with average methylation beta values of 0.2 to 0.8, and lower reliability in type I probes or those within the promoter and CpG island regions. Importantly, we found that probe reliability has significant implications in the analyses of Epigenome-wide Association Studies (EWAS). Higher reliability is associated with more consistent effect sizes in different studies, the identification of differentially methylated regions (DMRs) and methylation quantitative trait locus (mQTLs), and significant correlations with downstream gene expression. Moreover, blood DNAm measurements obtained from probes with higher reliability are more likely to show concordance with brain DNAm measurements. Our findings, which provide crucial reliable information for probes on the EPIC array, will serve as a valuable resource for future DNAm studies.

Critical evaluation of the reliability of DNA methylation probes on the Illumina MethylationEPIC BeadChip microarrays

DNA methylation (DNAm) plays a crucial role in a number of complex diseases. However, the reliability of DNAm levels measured using Illumina arrays varies across different probes. Previous research primarily assessed probe reliability by comparing duplicate samples between the 450k-450k or 450k-EPIC platforms, with limited investigations on Illumina EPIC arrays. We conducted a comprehensive assessment of the EPIC array probe reliability using 138 duplicated blood DNAm samples generated by the Alzheimer's Disease Neuroimaging Initiative study. We introduced a novel statistical measure, the modified intraclass correlation, to better account for the disagreement in duplicate measurements. We observed higher reliability in probes with average methylation beta values of 0.2 to 0.8, and lower reliability in type I probes or those within the promoter and CpG island regions. Importantly, we found that probe reliability has significant implications in the analyses of Epigenome-wide Association Studies (EWAS). Higher reliability is associated with more consistent effect sizes in different studies, the identification of differentially methylated regions (DMRs) and methylation quantitative trait locus (mQTLs), and significant correlations with downstream gene expression. Moreover, blood DNAm measurements obtained from probes with higher reliability are more likely to show concordance with brain DNAm measurements. Our findings, which provide crucial reliable information for probes on the EPIC array, will serve as a valuable resource for future DNAm studies.

Genome-Wide DNA Methylation Profiles in Whole-Blood and Buccal Samples-Cross-Sectional, Longitudinal, and across Platforms

The field of DNA methylation research is rapidly evolving, focusing on disease and phenotype changes over time using methylation measurements from diverse tissue sources and multiple array platforms. Consequently, identifying the extent of longitudinal, inter-tissue, and inter-platform variation in DNA methylation is crucial for future advancement. DNA methylation was measured in 375 individuals, with 197 of those having 2 blood sample measurements ~10 years apart. Whole-blood samples were measured on Illumina Infinium 450K and EPIC methylation arrays, and buccal samples from a subset of 58 participants were measured on EPIC array. The data were analyzed with the aims to examine the correlation between methylation levels in longitudinal blood samples in 197 individuals, examine the correlation between methylation levels in the blood and buccal samples in 58 individuals, and examine the correlation between blood methylation profiles assessed on the EPIC and 450K arrays in 83 individuals. We identified 136,833, 7674, and 96,891 CpGs significantly and strongly correlated (>0.50) longitudinally, across blood and buccal samples as well as array platforms, respectively. A total of 3674 of these CpGs were shared across all three sets. Analysis of these shared CpGs identified previously found associations with aging, ancestry, and 7016 mQTLs as well.

Implication of DNA methylation during lifestyle mediated weight loss

Over the past 50 years, the number of overweight/obese people increased significantly, making obesity a global public health challenge. Apart from rare monogenic forms, obesity is a multifactorial disease, most likely resulting from a concerted interaction of genetic, epigenetic and environmental factors. Although recent studies opened new avenues in elucidating the complex genetics behind obesity, the biological mechanisms contributing to individual's risk to become obese are not yet fully understood. Non-genetic factors such as eating behaviour or physical activity are strong contributing factors for the onset of obesity. These factors may interact with genetic predispositions most likely via epigenetic mechanisms. Epigenome-wide association studies or methylome-wide association studies are measuring DNA methylation at single CpGs across thousands of genes and capture associations to obesity phenotypes such as BMI. However, they only represent a snapshot in the complex biological network and cannot distinguish between causes and consequences. Intervention studies are therefore a suitable method to control for confounding factors and to avoid possible sources of bias. In particular, intervention studies documenting changes in obesity-associated epigenetic markers during lifestyle driven weight loss, make an important contribution to a better understanding of epigenetic reprogramming in obesity. To investigate the impact of lifestyle in obesity state specific DNA methylation, especially concerning the development of new strategies for prevention and individual therapy, we reviewed 19 most recent human intervention studies. In summary, this review highlights the huge potential of targeted interventions to alter disease-associated epigenetic patterns. However, there is an urgent need for further robust and larger studies to identify the specific DNA methylation biomarkers which influence obesity.

DNA Methylation Signatures of Response to Conventional Synthetic and Biologic Disease-Modifying Antirheumatic Drugs (DMARDs) in Rheumatoid Arthritis

Rheumatoid arthritis (RA) is a complex condition that displays heterogeneity in disease severity and response to standard treatments between patients. Failure rates for conventional, target synthetic, and biologic disease-modifying rheumatic drugs (DMARDs) are significant. Although there are models for predicting patient response, they have limited accuracy, require replication/validation, or for samples to be obtained through a synovial biopsy. Thus, currently, there are no prediction methods approved for routine clinical use. Previous research has shown that genetics and environmental factors alone cannot explain the differences in response between patients. Recent studies have demonstrated that deoxyribonucleic acid (DNA) methylation plays an important role in the pathogenesis and disease progression of RA. Importantly, specific DNA methylation profiles associated with response to conventional, target synthetic, and biologic DMARDs have been found in the blood of RA patients and could potentially function as predictive biomarkers. This review will summarize and evaluate the evidence for DNA methylation signatures in treatment response mainly in blood but also learn from the progress made in the diseased tissue in cancer in comparison to RA and autoimmune diseases. We will discuss the benefits and challenges of using DNA methylation signatures as predictive markers and the potential for future progress in this area.

Identification of Cell-Specific Differential DNA Methylation Associated With Methotrexate Treatment Response in Rheumatoid Arthritis

OBJECTIVE

We undertook this study to estimate changes in cell-specific DNA methylation (DNAm) associated with methotrexate (MTX) response using whole blood samples collected from rheumatoid arthritis (RA) patients before and after initiation of MTX treatment.

METHODS

Patients included in this study were from the Rheumatoid Arthritis Medication Study (n = 66) and the University of California San Francisco Rheumatoid Arthritis study (n = 11). All patients met the American College of Rheumatology RA classification criteria. Blood samples were collected at baseline and following treatment. Disease Activity Scores in 28 joints using the C-reactive protein level were collected at baseline and after 3-6 months of treatment with MTX. Methylation profiles were generated using the Illumina Infinium HumanMethylation450 and MethylationEPIC v1.0 BeadChip arrays using DNA from whole blood. MTX response was defined using the EULAR response criteria (responders showed good/moderate response; nonresponders showed no response). Differentially methylated positions were identified using the Limma software package and Tensor Composition Analysis, which is a method for identifying cell-specific differential DNAm at the CpG level from tissue-level ("bulk") data. Differentially methylated regions were identified using Comb-p software.

RESULTS

We found evidence of differential global methylation between treatment response groups. Further, we found patterns of cell-specific differential global methylation associated with MTX response. After correction for multiple testing, 1 differentially methylated position was associated with differential DNAm between responders and nonresponders at baseline in CD4+ T cells, CD8+ T cells, and natural killer cells. Thirty-nine cell-specific differentially methylated regions associated with MTX treatment response were identified. There were no significant findings in analyses of whole blood samples.

CONCLUSION

We identified cell-specific changes in DNAm that were associated with MTX treatment response in RA patients. Future studies of DNAm and MTX treatment response should include measurements of DNAm from sorted cells.

Pooled analysis of epigenome-wide association studies of food consumption in KORA, TwinsUK and LLS

PURPOSE

Examining epigenetic patterns is a crucial step in identifying molecular changes of disease pathophysiology, with DNA methylation as the most accessible epigenetic measure. Diet is suggested to affect metabolism and health via epigenetic modifications. Thus, our aim was to explore the association between food consumption and DNA methylation.

METHODS

Epigenome-wide association studies were conducted in three cohorts: KORA FF4, TwinsUK, and Leiden Longevity Study, and 37 dietary exposures were evaluated. Food group definition was harmonized across the three cohorts. DNA methylation was measured using Infinium MethylationEPIC BeadChip in KORA and Infinium HumanMethylation450 BeadChip in the Leiden study and the TwinsUK study. Overall, data from 2293 middle-aged men and women were included. A fixed-effects meta-analysis pooled study-specific estimates. The significance threshold was set at 0.05 for false-discovery rate-adjusted p values per food group.

RESULTS

We identified significant associations between the methylation level of CpG sites and the consumption of onions and garlic (2), nuts and seeds (18), milk (1), cream (11), plant oils (4), butter (13), and alcoholic beverages (27). The signals targeted genes of metabolic health relevance, for example, GLI1, RPTOR, and DIO1, among others.

CONCLUSION

This EWAS is unique with its focus on food groups that are part of a Western diet. Significant findings were mostly related to food groups with a high-fat content.

Genetic variation in cis-regulatory domains suggests cell type-specific regulatory mechanisms in immunity

Studying the interplay between genetic variation, epigenetic changes, and regulation of gene expression is crucial to understand the modification of cellular states in various conditions, including immune diseases. In this study, we characterize the cell-specificity in three key cells of the human immune system by building cis maps of regulatory regions with coordinated activity (CRDs) from ChIP-seq peaks and methylation data. We find that only 33% of CRD-gene associations are shared between cell types, revealing how similarly located regulatory regions provide cell-specific modulation of gene activity. We emphasize important biological mechanisms, as most of our associations are enriched in cell-specific transcription factor binding sites, blood-traits, and immune disease-associated loci. Notably, we show that CRD-QTLs aid in interpreting GWAS findings and help prioritize variants for testing functional hypotheses within human complex diseases. Additionally, we map trans CRD regulatory associations, and among 207 trans-eQTLs discovered, 46 overlap with the QTLGen Consortium meta-analysis in whole blood, showing that mapping functional regulatory units using population genomics allows discovering important mechanisms in the regulation of gene expression in immune cells. Finally, we constitute a comprehensive resource describing multi-omics changes to gain a greater understanding of cell-type specific regulatory mechanisms of immunity.

A systematic evaluation of normalization methods and probe replicability using infinium EPIC methylation data

BACKGROUND

The Infinium EPIC array measures the methylation status of > 850,000 CpG sites. The EPIC BeadChip uses a two-array design: Infinium Type I and Type II probes. These probe types exhibit different technical characteristics which may confound analyses. Numerous normalization and pre-processing methods have been developed to reduce probe type bias as well as other issues such as background and dye bias.

METHODS

This study evaluates the performance of various normalization methods using 16 replicated samples and three metrics: absolute beta-value difference, overlap of non-replicated CpGs between replicate pairs, and effect on beta-value distributions. Additionally, we carried out Pearson's correlation and intraclass correlation coefficient (ICC) analyses using both raw and SeSAMe 2 normalized data.

RESULTS

The method we define as SeSAMe 2, which consists of the application of the regular SeSAMe pipeline with an additional round of QC, pOOBAH masking, was found to be the best performing normalization method, while quantile-based methods were found to be the worst performing methods. Whole-array Pearson's correlations were found to be high. However, in agreement with previous studies, a substantial proportion of the probes on the EPIC array showed poor reproducibility (ICC < 0.50). The majority of poor performing probes have beta values close to either 0 or 1, and relatively low standard deviations. These results suggest that probe reliability is largely the result of limited biological variation rather than technical measurement variation. Importantly, normalizing the data with SeSAMe 2 dramatically improved ICC estimates, with the proportion of probes with ICC values > 0.50 increasing from 45.18% (raw data) to 61.35% (SeSAMe 2).

Propensity scores as a novel method to guide sample allocation and minimize batch effects during the design of high throughput experiments

BACKGROUND

We developed a novel approach to minimize batch effects when assigning samples to batches. Our algorithm selects a batch allocation, among all possible ways of assigning samples to batches, that minimizes differences in average propensity score between batches. This strategy was compared to randomization and stratified randomization in a case-control study (30 per group) with a covariate (case vs control, represented as β1, set to be null) and two biologically relevant confounding variables (age, represented as β2, and hemoglobin A1c (HbA1c), represented as β3). Gene expression values were obtained from a publicly available dataset of expression data obtained from pancreas islet cells. Batch effects were simulated as twice the median biological variation across the gene expression dataset and were added to the publicly available dataset to simulate a batch effect condition. Bias was calculated as the absolute difference between observed betas under the batch allocation strategies and the true beta (no batch effects). Bias was also evaluated after adjustment for batch effects using ComBat as well as a linear regression model. In order to understand performance of our optimal allocation strategy under the alternative hypothesis, we also evaluated bias at a single gene associated with both age and HbA1c levels in the 'true' dataset (CAPN13 gene).

RESULTS

Pre-batch correction, under the null hypothesis (β1), maximum absolute bias and root mean square (RMS) of maximum absolute bias, were minimized using the optimal allocation strategy. Under the alternative hypothesis (β2 and β3 for the CAPN13 gene), maximum absolute bias and RMS of maximum absolute bias were also consistently lower using the optimal allocation strategy. ComBat and the regression batch adjustment methods performed well as the bias estimates moved towards the true values in all conditions under both the null and alternative hypotheses. Although the differences between methods were less pronounced following batch correction, estimates of bias (average and RMS) were consistently lower using the optimal allocation strategy under both the null and alternative hypotheses.

CONCLUSIONS

Our algorithm provides an extremely flexible and effective method for assigning samples to batches by exploiting knowledge of covariates prior to sample allocation.

Epigenetic-based age acceleration in a representative sample of older Americans: Associations with aging-related morbidity and mortality

Biomarkers developed from DNA methylation (DNAm) data are of growing interest as predictors of health outcomes and mortality in older populations. However, it is unknown how epigenetic aging fits within the context of known socioeconomic and behavioral associations with aging-related health outcomes in a large, population-based, and diverse sample. This study uses data from a representative, panel study of US older adults to examine the relationship between DNAm-based age acceleration measures in the prediction of cross-sectional and longitudinal health outcomes and mortality. We examine whether recent improvements to these scores, using principal component (PC)-based measures designed to remove some of the technical noise and unreliability in measurement, improve the predictive capability of these measures. We also examine how well DNAm-based measures perform against well-known predictors of health outcomes such as demographics, SES, and health behaviors. In our sample, age acceleration calculated using "second and third generation clocks," PhenoAge, GrimAge, and DunedinPACE, is consistently a significant predictor of health outcomes including cross-sectional cognitive dysfunction, functional limitations and chronic conditions assessed 2 y after DNAm measurement, and 4-y mortality. PC-based epigenetic age acceleration measures do not significantly change the relationship of DNAm-based age acceleration measures to health outcomes or mortality compared to earlier versions of these measures. While the usefulness of DNAm-based age acceleration as a predictor of later life health outcomes is quite clear, other factors such as demographics, SES, mental health, and health behaviors remain equally, if not more robust, predictors of later life outcomes.

Refining epigenetic prediction of chronological and biological age

BACKGROUND

Epigenetic clocks can track both chronological age (cAge) and biological age (bAge). The latter is typically defined by physiological biomarkers and risk of adverse health outcomes, including all-cause mortality. As cohort sample sizes increase, estimates of cAge and bAge become more precise. Here, we aim to develop accurate epigenetic predictors of cAge and bAge, whilst improving our understanding of their epigenomic architecture.

METHODS

First, we perform large-scale (N = 18,413) epigenome-wide association studies (EWAS) of chronological age and all-cause mortality. Next, to create a cAge predictor, we use methylation data from 24,674 participants from the Generation Scotland study, the Lothian Birth Cohorts (LBC) of 1921 and 1936, and 8 other cohorts with publicly available data. In addition, we train a predictor of time to all-cause mortality as a proxy for bAge using the Generation Scotland cohort (1214 observed deaths). For this purpose, we use epigenetic surrogates (EpiScores) for 109 plasma proteins and the 8 component parts of GrimAge, one of the current best epigenetic predictors of survival. We test this bAge predictor in four external cohorts (LBC1921, LBC1936, the Framingham Heart Study and the Women's Health Initiative study).

RESULTS

Through the inclusion of linear and non-linear age-CpG associations from the EWAS, feature pre-selection in advance of elastic net regression, and a leave-one-cohort-out (LOCO) cross-validation framework, we obtain cAge prediction with a median absolute error equal to 2.3 years. Our bAge predictor was found to slightly outperform GrimAge in terms of the strength of its association to survival (HRGrimAge = 1.47 [1.40, 1.54] with p = 1.08 × 10-52, and HRbAge = 1.52 [1.44, 1.59] with p = 2.20 × 10-60). Finally, we introduce MethylBrowsR, an online tool to visualise epigenome-wide CpG-age associations.

CONCLUSIONS

The integration of multiple large datasets, EpiScores, non-linear DNAm effects, and new approaches to feature selection has facilitated improvements to the blood-based epigenetic prediction of biological and chronological age.

Effect of mechanical unloading on genome-wide DNA methylation profile of the failing human heart

Heart failure (HF) is characterized by global alterations in myocardial DNA methylation, yet little is known about the epigenetic regulation of the noncoding genome and potential reversibility of DNA methylation with left ventricular assist device (LVAD) therapy. Genome-wide mapping of myocardial DNA methylation in 36 patients with HF at LVAD implantation, 8 patients at LVAD explantation, and 7 nonfailing (NF) donors using a high-density bead array platform identified 2,079 differentially methylated positions (DMPs) in ischemic cardiomyopathy (ICM) and 261 DMPs in nonischemic cardiomyopathy (NICM). LVAD support resulted in normalization of 3.2% of HF-associated DMPs. Methylation-expression correlation analysis yielded several protein-coding genes that are hypomethylated and upregulated (HTRA1, FBXO16, EFCAB13, and AKAP13) or hypermethylated and downregulated (TBX3) in HF. A potentially novel cardiac-specific super-enhancer long noncoding RNA (lncRNA) (LINC00881) is hypermethylated and downregulated in human HF. LINC00881 is an upstream regulator of sarcomere and calcium channel gene expression including MYH6, CACNA1C, and RYR2. LINC00881 knockdown reduces peak calcium amplitude in the beating human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs). These data suggest that HF-associated changes in myocardial DNA methylation within coding and noncoding genomes are minimally reversible with mechanical unloading. Epigenetic reprogramming strategies may be necessary to achieve sustained clinical recovery from heart failure.

Genome-wide placental DNA methylations in fetal overgrowth and associations with leptin, adiponectin and fetal growth factors

BACKGROUND

Fetal overgrowth "programs" an elevated risk of type 2 diabetes in adulthood. Epigenetic alterations may be a mechanism in programming the vulnerability. We sought to characterize genome-wide alterations in placental gene methylations in fetal overgrowth and the associations with metabolic health biomarkers including leptin, adiponectin and fetal growth factors.

RESULTS

Comparing genome-wide placental gene DNA methylations in large-for-gestational-age (LGA, an indicator of fetal overgrowth, n = 30) versus optimal-for-gestational-age (OGA, control, n = 30) infants using the Illumina Infinium Human Methylation-EPIC BeadChip, we identified 543 differential methylation positions (DMPs; 397 hypermethylated, 146 hypomethylated) at false discovery rate < 5% and absolute methylation difference > 0.05 after adjusting for placental cell-type heterogeneity, maternal age, pre-pregnancy BMI and HbA1c levels during pregnancy. Twenty-five DMPs annotated to 20 genes (QSOX1, FCHSD2, LOC101928162, ADGRB3, GCNT1, TAP1, MYO16, NAV1, ATP8A2, LBXCOR1, EN2, INCA1, CAMTA2, SORCS2, SLC4A4, RPA3, UMAD1,USP53, OR2L13 and NR3C2) could explain 80% of the birth weight variations. Pathway analyses did not detect any statistically significant pathways after correcting for multiple tests. We validated a newly discovered differentially (hyper-)methylated gene-visual system homeobox 1 (VSX1) in an independent pyrosequencing study sample (LGA 47, OGA 47). Our data confirmed a hypermethylated gene-cadherin 13 (CDH13) reported in a previous epigenome-wide association study. Adiponectin in cord blood was correlated with its gene methylation in the placenta, while leptin and fetal growth factors (insulin, IGF-1, IGF-2) were not.

CONCLUSIONS

Fetal overgrowth may be associated with a large number of altered placental gene methylations. Placental VSX1 and CDH13 genes are hypermethylated in fetal overgrowth. Placental ADIPOQ gene methylations and fetal circulating adiponectin levels were correlated, suggesting the contribution of placenta-originated adiponectin to cord blood adiponectin.

Measuring biological age using omics data

Age is the key risk factor for diseases and disabilities of the elderly. Efforts to tackle age-related diseases and increase healthspan have suggested targeting the ageing process itself to 'rejuvenate' physiological functioning. However, achieving this aim requires measures of biological age and rates of ageing at the molecular level. Spurred by recent advances in high-throughput omics technologies, a new generation of tools to measure biological ageing now enables the quantitative characterization of ageing at molecular resolution. Epigenomic, transcriptomic, proteomic and metabolomic data can be harnessed with machine learning to build 'ageing clocks' with demonstrated capacity to identify new biomarkers of biological ageing.

An evaluation of the genome-wide false positive rates of common methods for identifying differentially methylated regions using illumina methylation arrays

Differentially methylated regions (DMRs) are genomic regions with specific methylation patterns across multiple loci that are associated with a phenotype. We examined the genome-wide false positive (GFP) rates of five widely used DMR methods: comb-p, Bumphunter, DMRcate, mCSEA and coMethDMR using both Illumina HumanMethylation450 (450 K) and MethylationEPIC (EPIC) data and simulated continuous and dichotomous null phenotypes (i.e., generated independently of methylation data). coMethDMR provided well-controlled GFP rates (~5%) except when analysing skewed continuous phenotypes. DMRcate generally had well-controlled GFP rates when applied to 450 K data except for the skewed continuous phenotype and EPIC data only for the normally distributed continuous phenotype. GFP rates for mCSEA were at least 0.096 and comb-p yielded GFP rates above 0.34. Bumphunter had high GFP rates of at least 0.35 across conditions, reaching as high as 0.95. Analysis of the performance of these methods in specific regions of the genome found that regions with higher correlation across loci had higher regional false positive rates on average across methods. Based on the false positive rates, coMethDMR is the most recommended analysis method, and DMRcate had acceptable performance when analysing 450 K data. However, as both could display higher levels of FPs for skewed continuous distributions, a normalizing transformation of skewed continuous phenotypes is suggested. This study highlights the importance of genome-wide simulations when evaluating the performance of DMR-analysis methods.

Identification of influential probe types in epigenetic predictions of human traits: implications for microarray design

BACKGROUND

CpG methylation levels can help to explain inter-individual differences in phenotypic traits. Few studies have explored whether identifying probe subsets based on their biological and statistical properties can maximise predictions whilst minimising array content. Variance component analyses and penalised regression (epigenetic predictors) were used to test the influence of (i) the number of probes considered, (ii) mean probe variability and (iii) methylation QTL status on the variance captured in eighteen traits by blood DNA methylation. Training and test samples comprised ≤ 4450 and ≤ 2578 unrelated individuals from Generation Scotland, respectively.

RESULTS

As the number of probes under consideration decreased, so too did the estimates from variance components and prediction analyses. Methylation QTL status and mean probe variability did not influence variance components. However, relative effect sizes were 15% larger for epigenetic predictors based on probes with known or reported methylation QTLs compared to probes without reported methylation QTLs. Relative effect sizes were 45% larger for predictors based on probes with mean Beta-values between 10 and 90% compared to those based on hypo- or hypermethylated probes (Beta-value ≤ 10% or ≥ 90%).

CONCLUSIONS

Arrays with fewer probes could reduce costs, leading to increased sample sizes for analyses. Our results show that reducing array content can restrict prediction metrics and careful attention must be given to the biological and distribution properties of CpG probes in array content selection.

A computational solution for bolstering reliability of epigenetic clocks: Implications for clinical trials and longitudinal tracking

Epigenetic clocks are widely used aging biomarkers calculated from DNA methylation data, but this data can be surprisingly unreliable. Here we show technical noise produces deviations up to 9 years between replicates for six prominent epigenetic clocks, limiting their utility. We present a computational solution to bolster reliability, calculating principal components from CpG-level data as input for biological age prediction. Our retrained principal-component versions of six clocks show agreement between most replicates within 1.5 years, improved detection of clock associations and intervention effects, and reliable longitudinal trajectories in vivo and in vitro. This method entails only one additional step compared to traditional clocks, requires no replicates or prior knowledge of CpG reliabilities for training, and can be applied to any existing or future epigenetic biomarker. The high reliability of principal component-based clocks is critical for applications to personalized medicine, longitudinal tracking, in vitro studies, and clinical trials of aging interventions.

Low reliability of DNA methylation across Illumina Infinium platforms in cord blood: implications for replication studies and meta-analyses of prenatal exposures

Background

There is an increasing interest in the role of epigenetics in epidemiology, but the emerging research field faces several critical biological and technical challenges. In particular, recent studies have shown poor correlation of measured DNA methylation (DNAm) levels within and across Illumina Infinium platforms in various tissues. In this study, we have investigated concordance between 450 k and EPIC Infinium platforms in cord blood. We could not replicate our previous findings on the association of prenatal paracetamol exposure with cord blood DNAm, which prompted an investigation of cross-platform DNAm differences.

Results

This study is based on two DNAm data sets from cord blood samples selected from the Norwegian Mother, Father and Child Cohort Study (MoBa). DNAm of one data set was measured using the 450 k platform and the other data set was measured using the EPIC platform. Initial analyses of the EPIC data could not replicate any of our previous significant findings in the 450 k data on associations between prenatal paracetamol exposure and cord blood DNAm. A subset of the samples (n = 17) was included in both data sets, which enabled analyses of technical sources potentially contributing to the negative replication. Analyses of these 17 samples with repeated measurements revealed high per-sample correlations (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\stackrel{\mathrm{-}}{\text{R}}\hspace{0.17em}\approx$$\end{document}R-≈ 0.99), but low per-CpG correlations (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\stackrel{\mathrm{-}}{\text{R}}$$\end{document}R- ≈ 0.24) between the platforms. 1.7% of the CpGs exhibited a mean DNAm difference across platforms > 0.1. Furthermore, only 26.7% of the CpGs exhibited a moderate or better cross-platform reliability (intra-class correlation coefficient ≥ 0.5).

Conclusion

The observations of low cross-platform probe correlation and reliability corroborate previous reports in other tissues. Our study cannot determine the origin of the differences between platforms. Nevertheless, it emulates the setting in studies using data from multiple Infinium platforms, often analysed several years apart. Therefore, the findings may have important implications for future epigenome-wide association studies (EWASs), in replication, meta-analyses and longitudinal studies. Cognisance and transparency of the challenges related to cross-platform studies may enhance the interpretation, replicability and validity of EWAS results both in cord blood and other tissues, ultimately improving the clinical relevance of epigenetic epidemiology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-022-01299-3.

Longitudinal Study of DNA Methylation and Epigenetic Clocks Prior to and Following Test-Confirmed COVID-19 and mRNA Vaccination

The host epigenetic landscape rapidly changes during SARS-CoV-2 infection, and evidence suggest that severe COVID-19 is associated with durable scars to the epigenome. Specifically, aberrant DNA methylation changes in immune cells and alterations to epigenetic clocks in blood relate to severe COVID-19. However, a longitudinal assessment of DNA methylation states and epigenetic clocks in blood from healthy individuals prior to and following test-confirmed non-hospitalized COVID-19 has not been performed. Moreover, the impact of mRNA COVID-19 vaccines upon the host epigenome remains understudied. Here, we first examined DNA methylation states in the blood of 21 participants prior to and following test-confirmed COVID-19 diagnosis at a median time frame of 8.35 weeks; 756 CpGs were identified as differentially methylated following COVID-19 diagnosis in blood at an FDR adjusted p-value < 0.05. These CpGs were enriched in the gene body, and the northern and southern shelf regions of genes involved in metabolic pathways. Integrative analysis revealed overlap among genes identified in transcriptional SARS-CoV-2 infection datasets. Principal component-based epigenetic clock estimates of PhenoAge and GrimAge significantly increased in people over 50 following infection by an average of 2.1 and 0.84 years. In contrast, PCPhenoAge significantly decreased in people fewer than 50 following infection by an average of 2.06 years. This observed divergence in epigenetic clocks following COVID-19 was related to age and immune cell-type compositional changes in CD4+ T cells, B cells, granulocytes, plasmablasts, exhausted T cells, and naïve T cells. Complementary longitudinal epigenetic clock analyses of 36 participants prior to and following Pfizer and Moderna mRNA-based COVID-19 vaccination revealed that vaccination significantly reduced principal component-based Horvath epigenetic clock estimates in people over 50 by an average of 3.91 years for those who received Moderna. This reduction in epigenetic clock estimates was significantly related to chronological age and immune cell-type compositional changes in B cells and plasmablasts pre- and post-vaccination. These findings suggest the potential utility of epigenetic clocks as a biomarker of COVID-19 vaccine responses. Future research will need to unravel the significance and durability of short-term changes in epigenetic age related to COVID-19 exposure and mRNA vaccination.

Epigenetic mechanisms of lung carcinogenesis involve differentially methylated CpG sites beyond those associated with smoking

Smoking-related epigenetic changes have been linked to lung cancer, but the contribution of epigenetic alterations unrelated to smoking remains unclear. We sought for a sparse set of CpG sites predicting lung cancer and explored the role of smoking in these associations. We analysed CpGs in relation to lung cancer in participants from two nested case-control studies, using (LASSO)-penalised regression. We accounted for the effects of smoking using known smoking-related CpGs, and through conditional-independence network. We identified 29 CpGs (8 smoking-related, 21 smoking-unrelated) associated with lung cancer. Models additionally adjusted for Comprehensive Smoking Index-(CSI) selected 1 smoking-related and 49 smoking-unrelated CpGs. Selected CpGs yielded excellent discriminatory performances, outperforming information provided by CSI only. Of the 8 selected smoking-related CpGs, two captured lung cancer-relevant effects of smoking that were missed by CSI. Further, the 50 CpGs identified in the CSI-adjusted model complementarily explained lung cancer risk. These markers may provide further insight into lung cancer carcinogenesis and help improving early identification of high-risk patients.

Methylation at CpG sites related to growth differentiation factor-15 was not prospectively associated with cardiovascular death in discordant monozygotic twins

Myocardial infarction patients had decreased methylation at four growth differentiating factor-15 (GDF-15) related CpG sites (cg13033858, cg16936953, cg17150809, and cg18608055). These sites had not been studied for their association with cardiovascular disease (CVD) deaths. Thus, we aimed to assess the associations independent of genes, shared environment, and traditional CVD risk factors. Nineteen white, male, monozygotic twin pairs discordant for CVD deaths were included from the National Heart, Lung and Blood Institute Twin Study (NHLBI) initiated in 1969. Data on vital status was collected through December 31, 2014. Methylation of buffy coat DNA at exam 3 (1986–87) was measured using the Illumina HumanMethylation450 BeadChip. Principal component analysis was used to generate a score representing blood leukocyte composition and baseline CVD risk factors and predominated with natural killer cells, CD4+ T cells, and Framingham risk score. Conditional logistic regression demonstrated that methylation at the four CpG sites was not associated with CVD deaths before (all p > 0.05, bootstrapped p > 0.05) and after adjustment for the score (all p > 0.05). Joint influences of cg16936953 and the score were statistically significant (p < 0.05). In conclusion, joint influences of methylation at the site cg16936953 and the score are prospectively associated with CVD deaths independent of germline and common environment.

ClinicalTrials.gov Identifier for NHLBI Twin Study: NCT00005124.

Environmentally sensitive hotspots in the methylome of the early human embryo

In humans, DNA methylation marks inherited from gametes are largely erased following fertilisation, prior to construction of the embryonic methylome. Exploiting a natural experiment of seasonal variation including changes in diet and nutritional status in rural Gambia, we analysed three datasets covering two independent child cohorts and identified 259 CpGs showing consistent associations between season of conception (SoC) and DNA methylation. SoC effects were most apparent in early infancy, with evidence of attenuation by mid-childhood. SoC-associated CpGs were enriched for metastable epialleles, parent-of-origin-specific methylation and germline differentially methylated regions, supporting a periconceptional environmental influence. Many SoC-associated CpGs overlapped enhancers or sites of active transcription in H1 embryonic stem cells and fetal tissues. Half were influenced but not determined by measured genetic variants that were independent of SoC. Environmental 'hotspots' providing a record of environmental influence at periconception constitute a valuable resource for investigating epigenetic mechanisms linking early exposures to lifelong health and disease.

Integrating DNA Methylation Measures of Biological Aging into Social Determinants of Health Research

PURPOSE OF REVIEW

Acceleration of biological processes of aging is hypothesized to drive excess morbidity and mortality in socially disadvantaged populations. DNA methylation measures of biological aging provide tools for testing this hypothesis.

RECENT FINDINGS

Next-generation DNA methylation measures of biological aging developed to predict mortality risk and physiological decline are more predictive of morbidity and mortality than the original epigenetic clocks developed to predict chronological age. These new measures show consistent evidence of more advanced and faster biological aging in people exposed to socioeconomic disadvantage and may be able to record the emergence of socially determined health inequalities as early as childhood. Next-generation DNA methylation measures of biological aging also indicate race/ethnic disparities in biological aging. More research is needed on these measures in samples of non-Western and non-White populations. New DNA methylation measures of biological aging open opportunities for refining inference about the causes of social disparities in health and devising policies to eliminate them. Further refining measures of biological aging by including more diversity in samples used for measurement development is a critical priority for the field.

Genome-Wide Placental Gene Methylations in Gestational Diabetes Mellitus, Fetal Growth and Metabolic Health Biomarkers in Cord Blood

Gestational diabetes mellitus (GDM) "program" an elevated risk of metabolic syndrome in the offspring. Epigenetic alterations are a suspected mechanism. GDM has been associated with placental DNA methylation changes in some epigenome-wide association studies. It remains unclear which genes or pathways are affected, and whether any placental differential gene methylations are correlated to fetal growth or circulating metabolic health biomarkers. In an epigenome-wide association study using the Infinium MethylationEPIC Beadchip, we sought to identify genome-wide placental differentially methylated genes and enriched pathways in GDM, and to assess the correlations with fetal growth and metabolic health biomarkers in cord blood. The study samples were 30 pairs of term placentas in GDM vs. euglycemic pregnancies (controls) matched by infant sex and gestational age at delivery in the Shanghai Birth Cohort. Cord blood metabolic health biomarkers included insulin, C-peptide, proinsulin, IGF-I, IGF-II, leptin and adiponectin. Adjusting for maternal age, pre-pregnancy BMI, parity, mode of delivery and placental cell type heterogeneity, 256 differentially methylated positions (DMPs,130 hypermethylated and 126 hypomethylated) were detected between GDM and control groups accounting for multiple tests with false discovery rate <0.05 and beta-value difference >0.05. WSCD2 was identified as a differentially methylated gene in both site- and region-level analyses. We validated 7 hypermethylated (CYP1A2, GFRA1, HDAC4, LIMS2, NAV3, PAX6, UPK1B) and 10 hypomethylated (DPP10, CPLX1, CSMD2, GPR133, NRXN1, PCSK9, PENK, PRDM16, PTPRN2, TNXB) genes reported in previous epigenome-wide association studies. We did not find any enriched pathway accounting for multiple tests. DMPs in 11 genes (CYP2D7P1, PCDHB15, ERG, SIRPB1, DKK2, RAPGEF5, CACNA2D4, PCSK9, TSNARE1, CADM2, KCNAB2) were correlated with birth weight (z score) accounting for multiple tests. There were no significant correlations between placental gene methylations and cord blood biomarkers. In conclusions, GDM was associated with DNA methylation changes in a number of placental genes, but these placental gene methylations were uncorrelated to the observed metabolic health biomarkers (fetal growth factors, leptin and adiponectin) in cord blood. We validated 17 differentially methylated placental genes in GDM, and identified 11 differentially methylated genes relevant to fetal growth.

The SEQC2 epigenomics quality control (EpiQC) study

BACKGROUND

Cytosine modifications in DNA such as 5-methylcytosine (5mC) underlie a broad range of developmental processes, maintain cellular lineage specification, and can define or stratify types of cancer and other diseases. However, the wide variety of approaches available to interrogate these modifications has created a need for harmonized materials, methods, and rigorous benchmarking to improve genome-wide methylome sequencing applications in clinical and basic research. Here, we present a multi-platform assessment and cross-validated resource for epigenetics research from the FDA's Epigenomics Quality Control Group.

RESULTS

Each sample is processed in multiple replicates by three whole-genome bisulfite sequencing (WGBS) protocols (TruSeq DNA methylation, Accel-NGS MethylSeq, and SPLAT), oxidative bisulfite sequencing (TrueMethyl), enzymatic deamination method (EMSeq), targeted methylation sequencing (Illumina Methyl Capture EPIC), single-molecule long-read nanopore sequencing from Oxford Nanopore Technologies, and 850k Illumina methylation arrays. After rigorous quality assessment and comparison to Illumina EPIC methylation microarrays and testing on a range of algorithms (Bismark, BitmapperBS, bwa-meth, and BitMapperBS), we find overall high concordance between assays, but also differences in efficiency of read mapping, CpG capture, coverage, and platform performance, and variable performance across 26 microarray normalization algorithms.

CONCLUSIONS

The data provided herein can guide the use of these DNA reference materials in epigenomics research, as well as provide best practices for experimental design in future studies. By leveraging seven human cell lines that are designated as publicly available reference materials, these data can be used as a baseline to advance epigenomics research.

An effective processing pipeline for harmonizing DNA methylation data from Illumina's 450K and EPIC platforms for epidemiological studies

OBJECTIVE

Illumina BeadChip arrays are commonly used to generate DNA methylation data for large epidemiological studies. Updates in technology over time create challenges for data harmonization within and between studies, many of which obtained data from the older 450K and newer EPIC platforms. The pre-processing pipeline for DNA methylation is not trivial, and influences the downstream analyses. Incorporating different platforms adds a new level of technical variability that has not yet been taken into account by recommended pipelines. Our study evaluated the performance of various tools on different versions of platform data harmonization at each step of pre-processing pipeline, including quality control (QC), normalization, batch effect adjustment, and genomic inflation. We illustrate our novel approach using 450K and EPIC data from the Diabetes Autoimmunity Study in the Young (DAISY) prospective cohort.

RESULTS

We found normalization and probe filtering had the biggest effect on data harmonization. Employing a meta-analysis was an effective and easily executable method for accounting for platform variability. Correcting for genomic inflation also helped with harmonization. We present guidelines for studies seeking to harmonize data from the 450K and EPIC platforms, which includes the use of technical replicates for evaluating numerous pre-processing steps, and employing a meta-analysis.

Identifying causal models between genetically regulated methylation patterns and gene expression in healthy colon tissue

BACKGROUND

DNA methylation is involved in the regulation of gene expression and phenotypic variation, but the inter-relationship between genetic variation, DNA methylation and gene expression remains poorly understood. Here we combine the analysis of genetic variants related to methylation markers (methylation quantitative trait loci: mQTLs) and gene expression (expression quantitative trait loci: eQTLs) with methylation markers related to gene expression (expression quantitative trait methylation: eQTMs), to provide novel insights into the genetic/epigenetic architecture of colocalizing molecular markers.

RESULTS

Normal mucosa from 100 patients with colon cancer and 50 healthy donors included in the Colonomics project have been analyzed. Linear models have been used to find mQTLs and eQTMs within 1 Mb of the target gene. From 32,446 eQTLs previously detected, we found a total of 6850 SNPs, 114 CpGs and 52 genes interrelated, generating 13,987 significant combinations of co-occurring associations (meQTLs) after Bonferromi correction. Non-redundant meQTLs were 54, enriched in genes involved in metabolism of glucose and xenobiotics and immune system. SNPs in meQTLs were enriched in regulatory elements (enhancers and promoters) compared to random SNPs within 1 Mb of genes. Three colorectal cancer GWAS SNPs were related to methylation changes, and four SNPs were related to chemerin levels. Bayesian networks have been used to identify putative causal relationships among associated SNPs, CpG and gene expression triads. We identified that most of these combinations showed the canonical pathway of methylation markers causes gene expression variation (60.1%) or non-causal relationship between methylation and gene expression (33.9%); however, in up to 6% of these combinations, gene expression was causing variation in methylation markers.

CONCLUSIONS

In this study we provided a characterization of the regulation between genetic variants and inter-dependent methylation markers and gene expression in a set of 150 healthy colon tissue samples. This is an important finding for the understanding of molecular susceptibility on colon-related complex diseases.

Risky family climates presage increased cellular aging in young adulthood

A scientific consensus is emerging that children reared in risky family climates are prone to chronic diseases and premature death later in life. Few prospective data, however, are available to inform the mechanisms of these relationships. In a prospective study involving 323 Black families, we sought to determine whether, and how, childhood risky family climates are linked to a potential risk factor for later-life disease: increases in cellular aging (indexed by epigenetic aging). As hypothesized, risky family climates were associated with greater outflows of the stress hormones epinephrine and norepinephrine at ages 19 and 20 years; this, in turn, led to increases in cellular aging across ages 20-27 years. If sustained, these tendencies may place children from risky family climates on a trajectory toward the chronic diseases of aging.

CUE: CpG impUtation ensemble for DNA methylation levels across the human methylation450 (HM450) and EPIC (HM850) BeadChip platforms

DNA methylation at CpG dinucleotides is one of the most extensively studied epigenetic marks. With technological advancements, geneticists can profile DNA methylation with multiple reliable approaches. However, profiling platforms can differ substantially in the CpGs they assess, consequently hindering integrated analysis across platforms. Here, we present CpG impUtation Ensemble (CUE), which leverages multiple classical statistical and modern machine learning methods, to impute from the Illumina HumanMethylation450 (HM450) BeadChip to the Illumina HumanMethylationEPIC (HM850) BeadChip. Data were analysed from two population cohorts with methylation measured both by HM450 and HM850: the Extremely Low Gestational Age Newborns (ELGAN) study (n = 127, placenta) and the VA Boston Posttraumatic Stress Disorder (PTSD) genetics repository (n = 144, whole blood). Cross-validation results show that CUE achieves the lowest predicted root-mean-square error (RMSE) (0.026 in PTSD) and the highest accuracy (99.97% in PTSD) compared with five individual methods tested, including k-nearest-neighbours, logistic regression, penalized functional regression, random forest, and XGBoost. Finally, among all 339,033 HM850-only CpG sites shared between ELGAN and PTSD, CUE successfully imputed 289,604 (85.4%) sites, where success was defined as RMSE < 0.05 and accuracy >95% in PTSD. In summary, CUE is a valuable tool for imputing CpG methylation from the HM450 to HM850 platform.

Severity of Idiopathic Scoliosis Is Associated with Differential Methylation: An Epigenome-Wide Association Study of Monozygotic Twins with Idiopathic Scoliosis

Epigenetic mechanisms may contribute to idiopathic scoliosis (IS). We identified 8 monozygotic twin pairs with IS, 6 discordant (Cobb angle difference > 10°) and 2 concordant (Cobb angle difference ≤ 2°). Genome-wide methylation in blood was measured with the Infinium HumanMethylation EPIC Beadchip. We tested for differences in methylation and methylation variability between discordant twins and tested the association between methylation and curve severity in all twins. Differentially methylated region (DMR) analyses identified gene promoter regions. Methylation at cg12959265 (chr. 7 DPY19L1) was less variable in cases (false discovery rate (FDR) = 0.0791). We identified four probes (false discovery rate, FDR < 0.10); cg02477677 (chr. 17, RARA gene), cg12922161 (chr. 2 LOC150622 gene), cg08826461 (chr. 2), and cg16382077 (chr. 7) associated with curve severity. We identified 57 DMRs where hyper- or hypo-methylation was consistent across the region and 28 DMRs with a consistent association with curve severity. Among DMRs, 21 were correlated with bone methylation. Prioritization of regions based on methylation concordance in bone identified promoter regions for WNT10A (WNT signaling), NPY (regulator of bone and energy homeostasis), and others predicted to be relevant for bone formation/remodeling. These regions may aid in understanding the complex interplay between genetics, environment, and IS.

Epigenetic Biomarkers of Prenatal Tobacco Smoke Exposure Are Associated with Gene Deletions in Childhood Acute Lymphoblastic Leukemia

BACKGROUND

Parental smoking is implicated in the etiology of acute lymphoblastic leukemia (ALL), the most common childhood cancer. We recently reported an association between an epigenetic biomarker of early-life tobacco smoke exposure at the AHRR gene and increased frequency of somatic gene deletions among ALL cases.

METHODS

Here, we further assess this association using two epigenetic biomarkers for maternal smoking during pregnancy-DNA methylation at AHRR CpG cg05575921 and a recently established polyepigenetic smoking score-in an expanded set of 482 B-cell ALL (B-ALL) cases in the California Childhood Leukemia Study with available Illumina 450K or MethylationEPIC array data. Multivariable Poisson regression models were used to test the associations between the epigenetic biomarkers and gene deletion numbers.

RESULTS

We found an association between DNA methylation at AHRR CpG cg05575921 and deletion number among 284 childhood B-ALL cases with MethylationEPIC array data, with a ratio of means (RM) of 1.31 [95% confidence interval (CI), 1.02-1.69] for each 0.1 β value reduction in DNA methylation, an effect size similar to our previous report in an independent set of 198 B-ALL cases with 450K array data [meta-analysis summary RM (sRM) = 1.32; 95% CI, 1.10-1.57]. The polyepigenetic smoking score was positively associated with gene deletion frequency among all 482 B-ALL cases (sRM = 1.31 for each 4-unit increase in score; 95% CI, 1.09-1.57).

CONCLUSIONS

We provide further evidence that prenatal tobacco-smoke exposure may influence the generation of somatic copy-number deletions in childhood B-ALL.

IMPACT

Analyses of deletion breakpoint sequences are required to further understand the mutagenic effects of tobacco smoke in childhood ALL.

High-resolution targeted bisulfite sequencing reveals blood cell type-specific DNA methylation patterns in IL13 and ORMDL3

BACKGROUND

Methylation of DNA at CpG sites is an epigenetic modification and a potential modifier of disease risk, possibly mediating environmental effects. Currently, DNA methylation is commonly assessed using specific microarrays that sample methylation at a few % of all methylated sites.

METHODS

To understand if significant information on methylation can be added by a more comprehensive analysis of methylation, we set up a quantitative method, bisulfite oligonucleotide-selective sequencing (Bs-OS-seq), and compared the data with microarray-derived methylation data. We assessed methylation at two asthma-associated genes, IL13 and ORMDL3, in blood samples collected from children with and without asthma and fractionated white blood cell types from healthy adult controls.

RESULTS

Our results show that Bs-OS-seq can uncover vast amounts of methylation variation not detected by commonly used array methods. We found that high-density methylation information from even one gene can delineate the main white blood cell lineages.

CONCLUSIONS

We conclude that high-resolution methylation studies can yield clinically important information at selected specific loci missed by array-based methods, with potential implications for future studies of methylation-disease associations.

Epigenetic prediction of 17β-estradiol and relationship to trauma-related outcomes in women

17β-estradiol (E2) levels in women correlate with multiple neuropsychiatric symptoms, including those that are stress-related. Furthermore, prior work from our group has demonstrated that E2 status influences DNA methylation (DNAm) across the genome. We developed and validated a DNAm-based predictor of E2 (one of four naturally occurring estrogens) using a training set of 183 females and a test set of 79 females from the same traumatized cohort. We showed that predicted E2 levels were highly correlated with measured E2 concentrations in our testing set (r ​= ​0.75, p ​= ​1.8e-15). We further demonstrated that predicted E2 concentrations, in combination with measured values, negatively correlated with current post-traumatic stress disorder (PTSD) (β ​= ​−0.38, p ​= ​0.01) and major depressive disorder (MDD) diagnoses (β ​= ​−0.45, p ​= ​0.02), as well as a continuous measure of PTSD symptom severity (β ​= ​−2.3, p ​= ​0.007) in females. Finally, we tested our predictor in an independent data set (n ​= ​85) also comprised of recently traumatized female subjects to determine if the predictor would generalize to a different population than the one on which it was developed. We found that the correlation between predicted and actual E2 concentrations in the external validation data set was also high (r ​= ​0.48, p ​= ​3.0e-6). While further validation is warranted, a DNAm predictor of E2 concentrations will advance our understanding of hormone-epigenetic interactions. Furthermore, such a DNAm predictor may serve as an epigenetic proxy for E2 concentrations and thus provide an important biomarker to better evaluate the contribution of E2 to current and potentially future psychiatric symptoms in samples for which E2 is not measured.

•

Neuropsychiatric symptoms correlate with estradiol (E2) levels in females.

•

We developed a DNA methylation-based E2 predictor using machine learning approach.

•

Our predictor performed well in an external validation cohort of traumatized women.

•

Predicted E2 concentrations correlated with stress-related phenotypes.

•

Our model may serve as an epigenetic biomarker of E2 status in women.

Genetic impacts on DNA methylation: research findings and future perspectives

Multiple recent studies highlight that genetic variants can have strong impacts on a significant proportion of the human DNA methylome. Methylation quantitative trait loci, or meQTLs, allow for the exploration of biological mechanisms that underlie complex human phenotypes, with potential insights for human disease onset and progression. In this review, we summarize recent milestones in characterizing the human genetic basis of DNA methylation variation over the last decade, including heritability findings and genome-wide identification of meQTLs. We also discuss challenges in this field and future areas of research geared to generate insights into molecular processes underlying human complex traits.

Epigenetic predictors of lifestyle traits applied to the blood and brain

Modifiable lifestyle factors influence the risk of developing many neurological diseases. These factors have been extensively linked with blood-based genome-wide DNA methylation, but it is unclear if the signatures from blood translate to the target tissue of interest-the brain. To investigate this, we apply blood-derived epigenetic predictors of four lifestyle traits to genome-wide DNA methylation from five post-mortem brain regions and the last blood sample prior to death in 14 individuals in the Lothian Birth Cohort 1936. Using these matched samples, we found that correlations between blood and brain DNA methylation scores for smoking, high-density lipoprotein cholesterol, alcohol and body mass index were highly variable across brain regions. Smoking scores in the dorsolateral prefrontal cortex had the strongest correlations with smoking scores in blood (r = 0.5, n = 14, P = 0.07) and smoking behaviour (r = 0.56, n = 9, P = 0.12). This was also the brain region which exhibited the largest correlations for DNA methylation at site cg05575921 - the single strongest correlate of smoking in blood-in relation to blood (r = 0.61, n = 14, P = 0.02) and smoking behaviour (r = -0.65, n = 9, P = 0.06). This suggested a particular vulnerability to smoking-related differential methylation in this region. Our work contributes to understanding how lifestyle factors affect the brain and suggest that lifestyle-related DNA methylation is likely to be both brain region dependent and in many cases poorly proxied for by blood. Though these pilot data provide a rarely-available opportunity for the comparison of methylation patterns across multiple brain regions and the blood, due to the limited sample size available our results must be considered as preliminary and should therefore be used as a basis for further investigation.

Epigenome wide association study of response to methotrexate in early rheumatoid arthritis patients

Aim

To identify differentially methylated positions (DMPs) and regions (DMRs) that predict response to Methotrexate (MTX) in early rheumatoid arthritis (RA) patients.

Materials and methods

DNA from baseline peripheral blood mononuclear cells was extracted from 72 RA patients. DNA methylation, quantified using the Infinium MethylationEPIC, was assessed in relation to response to MTX (combination) therapy over the first 3 months.

Results

Baseline DMPs associated with response were identified; including hits previously described in RA. Additionally, 1309 DMR regions were observed. However, none of these findings were genome-wide significant. Likewise, no specific pathways were related to response, nor could we replicate associations with previously identified DMPs.

Conclusion

No baseline genome-wide significant differences were identified as biomarker for MTX (combination) therapy response; hence meta-analyses are required.

Klotho, PTSD, and advanced epigenetic age in cortical tissue

This study examined the klotho (KL) longevity gene polymorphism rs9315202 and psychopathology, including posttraumatic stress disorder (PTSD), depression, and alcohol-use disorders, in association with advanced epigenetic age in three postmortem cortical tissue regions: dorsolateral and ventromedial prefrontal cortices and motor cortex. Using data from the VA National PTSD Brain Bank (n = 117), we found that rs9315202 interacted with PTSD to predict advanced epigenetic age in motor cortex among the subset of relatively older (>=45 years), white non-Hispanic decedents (corrected p = 0.014, n = 42). An evaluation of 211 additional common KL variants revealed that only variants in linkage disequilibrium with rs9315202 showed similarly high levels of significance. Alcohol abuse was nominally associated with advanced epigenetic age in motor cortex (p = 0.039, n = 114). The rs9315202 SNP interacted with PTSD to predict decreased KL expression via DNAm age residuals in motor cortex among older white non-Hispanics decedents (indirect β = -0.198, p = 0.027). Finally, in dual-luciferase enhancer reporter system experiments, we found that inserting the minor allele of rs9315202 in a human kidney cell line HK-2 genomic DNA resulted in a change in KL transcriptional activities, likely operating via long noncoding RNA in this region. This was the first study to examine multiple forms of psychopathology in association with advanced DNA methylation age across several brain regions, to extend work concerning the association between rs9315202 and advanced epigenetic to brain tissue, and to identify the effects of rs9315202 on KL gene expression. KL augmentation holds promise as a therapeutic intervention to slow the pace of cellular aging, disease onset, and neuropathology, particularly in older, stressed populations.