EWAS Atlas: a curated knowledgebase of epigenome-wide association studies

Direction-aware functional class scoring enrichment analysis of infinium DNA methylation data

Infinium Methylation BeadChip arrays remain one of the most popular platforms for epigenome-wide association studies, but tools for downstream pathway analysis have their limitations. Functional class scoring (FCS) is a group of pathway enrichment techniques that involve the ranking of genes and evaluation of their collective regulation in biological systems, but the implementations described for Infinium methylation array data do not retain direction information, which is important for mechanistic understanding of genomic regulation. Here, we evaluate several candidate FCS methods that retain directional information. According to simulation results, the best-performing method involves the mean aggregation of probe limma t-statistics by gene followed by a rank-ANOVA enrichment test using the mitch package. This method, which we call 'LAM,' outperformed an existing over-representation analysis method in simulations, and showed higher sensitivity and robustness in an analysis of real lung tumour-normal paired datasets. Using matched RNA-seq data, we examine the relationship of methylation differences at promoters and gene bodies with RNA expression at the level of pathways in lung cancer. To demonstrate the utility of our approach, we apply it to three other contexts where public data were available. First, we examine the differential pathway methylation associated with chronological age. Second, we investigate pathway methylation differences in infants conceived with in vitro fertilization. Lastly, we analyse differential pathway methylation in 19 disease states, identifying hundreds of novel associations. These results show LAM is a powerful method for the detection of differential pathway methylation complementing existing methods. A reproducible vignette is provided to illustrate how to implement this method.

DNA methylation heterogeneity attributable to a complex tumor immune microenvironment prompts prognostic risk in glioma

Gliomas are malignant tumours of the human nervous system with different World Health Organization (WHO) classifications, glioblastoma (GBM) with higher grade and are more malignant than lower-grade glioma (LGG). To dissect how the DNA methylation heterogeneity in gliomas is influenced by the complex cellular composition of the tumour immune microenvironment, we first compared the DNA methylation profiles of purified human immune cells and bulk glioma tissue, stratifying three tumour immune microenvironmental subtypes for GBM and LGG samples from The Cancer Genome Atlas (TCGA). We found that more intermediate methylation sites were enriched in glioma tumour tissues, and used the Proportion of sites with Intermediate Methylation (PIM) to compare intertumoral DNA methylation heterogeneity. A larger PIM score reflected stronger DNA methylation heterogeneity. Enhanced DNA methylation heterogeneity was associated with stronger immune cell infiltration, better survival rates, and slower tumour progression in glioma patients. We then created a Cell-type-associated DNA Methylation Heterogeneity Contribution (CMHC) score to explore the impact of different immune cell types on heterogeneous CpG site (CpGct) in glioma tissues. We identified eight prognosis-related CpGct to construct a risk score: the Cell-type-associated DNA Methylation Heterogeneity Risk (CMHR) score. CMHR was positively correlated with cytotoxic T-lymphocyte infiltration (CTL), and showed better predictive performance for IDH status (AUC = 0.96) and glioma histological phenotype (AUC = 0.81). Furthermore, DNA methylation alterations of eight CpGct might be related to drug treatments of gliomas. In conclusion, we indicated that DNA methylation heterogeneity is associated with a complex tumour immune microenvironment, glioma phenotype, and patient's prognosis.

Evaluating the association between placenta DNA methylation and cognitive functions in the offspring

The placenta plays a crucial role in protecting the fetus from environmental harm and supports the development of its brain. In fact, compromised placental function could predispose an individual to neurodevelopmental disorders. Placental epigenetic modifications, including DNA methylation, could be considered a proxy of placental function and thus plausible mediators of the association between intrauterine environmental exposures and genetics, and childhood and adult mental health. Although neurodevelopmental disorders such as autism spectrum disorder have been investigated in relation to placenta DNA methylation, no studies have addressed the association between placenta DNA methylation and child's cognitive functions. Thus, our goal here was to investigate whether the placental DNA methylation profile measured using the Illumina EPIC array is associated with three different cognitive domains (namely verbal score, perceptive performance score, and general cognitive score) assessed by the McCarthy Scales of Children's functions in childhood at age 4. To this end, we conducted epigenome-wide association analyses, including data from 255 mother-child pairs within the INMA project, and performed a follow-up functional analysis to help the interpretation of the findings. After multiple-testing correction, we found that methylation at 4 CpGs (cg1548200, cg02986379, cg00866476, and cg14113931) was significantly associated with the general cognitive score, and 2 distinct differentially methylated regions (DMRs) (including 27 CpGs) were significantly associated with each cognitive dimension. Interestingly, the genes annotated to these CpGs, such as DAB2, CEP76, PSMG2, or MECOM, are involved in placenta, fetal, and brain development. Moreover, functional enrichment analyses of suggestive CpGs (p < 1 × 10-4) revealed gene sets involved in placenta development, fetus formation, and brain growth. These findings suggest that placental DNA methylation could be a mechanism contributing to the alteration of important pathways in the placenta that have a consequence on the offspring's brain development and cognitive function.

Genetic effects on the skin methylome in healthy older twins

Whole-skin DNA methylation variation has been implicated in several diseases, including melanoma, but its genetic basis has not yet been fully characterized. Using bulk skin tissue samples from 414 healthy female UK twins, we performed twin-based heritability and methylation quantitative trait loci (meQTL) analyses for >400,000 DNA methylation sites. We find that the human skin DNA methylome is on average less heritable than previously estimated in blood and other tissues (mean heritability: 10.02%). meQTL analysis identified local genetic effects influencing DNA methylation at 18.8% (76,442) of tested CpG sites, as well as 1,775 CpG sites associated with at least one distal genetic variant. As a functional follow-up, we performed skin expression QTL (eQTL) analyses in a partially overlapping sample of 604 female twins. Colocalization analysis identified over 3,500 shared genetic effects affecting thousands of CpG sites (10,067) and genes (4,475). Mediation analysis of putative colocalized gene-CpG pairs identified 114 genes with evidence for eQTL effects being mediated by DNA methylation in skin, including in genes implicating skin disease such as ALOX12 and CSPG4. We further explored the relevance of skin meQTLs to skin disease and found that skin meQTLs and CpGs under genetic influence were enriched for multiple skin-related genome-wide and epigenome-wide association signals, including for melanoma and psoriasis. Our findings give insights into the regulatory landscape of epigenomic variation in skin.

Uncovering methylation-dependent genetic effects on regulatory element function in diverse genomes

A major goal in evolutionary biology and biomedicine is to understand the complex interactions between genetic variants, the epigenome, and gene expression. However, the causal relationships between these factors remain poorly understood. mSTARR-seq, a methylation-sensitive massively parallel reporter assay, is capable of identifying methylation-dependent regulatory activity at many thousands of genomic regions simultaneously, and allows for the testing of causal relationships between DNA methylation and gene expression on a region-by-region basis. Here, we developed a multiplexed mSTARR-seq protocol to assay naturally occurring human genetic variation from 25 individuals sampled from 10 localities in Europe and Africa. We identified 6,957 regulatory elements in either the unmethylated or methylated state, and this set was enriched for enhancer and promoter annotations, as expected. The expression of 58% of these regulatory elements was modulated by methylation, which was generally associated with decreased RNA expression. Within our set of regulatory elements, we used allele-specific expression analyses to identify 8,020 sites with genetic effects on gene regulation; further, we found that 42.3% of these genetic effects varied between methylated and unmethylated states. Sites exhibiting methylation-dependent genetic effects were enriched for GWAS and EWAS annotations, implicating them in human disease. Compared to datasets that assay DNA from a single European individual, our multiplexed assay uncovers dramatically more genetic effects and methylation-dependent genetic effects, highlighting the importance of including diverse individuals in assays which aim to understand gene regulatory processes.

Phenotype prediction using biologically interpretable neural networks on multi-cohort multi-omics data

Integrating multi-omics data into predictive models has the potential to enhance accuracy, which is essential for precision medicine. In this study, we developed interpretable predictive models for multi-omics data by employing neural networks informed by prior biological knowledge, referred to as visible networks. These neural networks offer insights into the decision-making process and can unveil novel perspectives on the underlying biological mechanisms associated with traits and complex diseases. We tested the performance, interpretability and generalizability for inferring smoking status, subject age and LDL levels using genome-wide RNA expression and CpG methylation data from the blood of the BIOS consortium (four population cohorts, Ntotal = 2940). In a cohort-wise cross-validation setting, the consistency of the diagnostic performance and interpretation was assessed. Performance was consistently high for predicting smoking status with an overall mean AUC of 0.95 (95% CI: 0.90-1.00) and interpretation revealed the involvement of well-replicated genes such as AHRR, GPR15 and LRRN3. LDL-level predictions were only generalized in a single cohort with an R2 of 0.07 (95% CI: 0.05-0.08). Age was inferred with a mean error of 5.16 (95% CI: 3.97-6.35) years with the genes COL11A2, AFAP1, OTUD7A, PTPRN2, ADARB2 and CD34 consistently predictive. For both regression tasks, we found that using multi-omics networks improved performance, stability and generalizability compared to interpretable single omic networks. We believe that visible neural networks have great potential for multi-omics analysis; they combine multi-omic data elegantly, are interpretable, and generalize well to data from different cohorts.

Maternal age is related to offspring DNA methylation: A meta-analysis of results from the PACE consortium

Worldwide trends to delay childbearing have increased parental ages at birth. Older parental age may harm offspring health, but mechanisms remain unclear. Alterations in offspring DNA methylation (DNAm) patterns could play a role as aging has been associated with methylation changes in gametes of older individuals. We meta-analyzed epigenome-wide associations of parental age with offspring blood DNAm of over 9500 newborns and 2000 children (5-10 years old) from the Pregnancy and Childhood Epigenetics consortium. In newborns, we identified 33 CpG sites in 13 loci with DNAm associated with maternal age (PFDR < 0.05). Eight of these CpGs were located near/in the MTNR1B gene, coding for a melatonin receptor. Regional analysis identified them together as a differentially methylated region consisting of 9 CpGs in/near MTNR1B, at which higher DNAm was associated with greater maternal age (PFDR = 6.92 × 10-8) in newborns. In childhood blood samples, these differences in blood DNAm of MTNR1B CpGs were nominally significant (p < 0.05) and retained the same positive direction, suggesting persistence of associations. Maternal age was also positively associated with higher DNA methylation at three CpGs in RTEL1-TNFRSF6B at birth (PFDR < 0.05) and nominally in childhood (p < 0.0001). Of the remaining 10 CpGs also persistent in childhood, methylation at cg26709300 in YPEL3/BOLA2B in external data was associated with expression of ITGAL, an immune regulator. While further study is needed to establish causality, particularly due to the small effect sizes observed, our results potentially support offspring DNAm as a mechanism underlying associations of maternal age with child health.

Multiomics Screening Identified CpG Sites and Genes That Mediate the Impact of Exposure to Environmental Chemicals on Cardiometabolic Traits

An understanding of the molecular mechanism whereby an environmental chemical causes a disease is important for the purposes of future applications. In this study, a multiomics workflow was designed to combine several publicly available datasets in order to identify CpG sites and genes that mediate the impact of exposure to environmental chemicals on cardiometabolic traits. Organophosphate and prenatal lead exposure were previously reported to change methylation level at the cg23627948 site. The outcome of the analyses conducted in this study revealed that, as the cg23627948 site becomes methylated, the expression of the GNA12 gene decreases, which leads to a higher body fat percentage. Prenatal perfluorooctane sulfonate exposure was reported to increase the methylation level at the cg21153102 site. Findings of this study revealed that higher methylation at this site contributes to higher diastolic blood pressure by changing the expression of CHP1 and GCHFR genes. Moreover, HKR1 mediates the impact of B12 supplementation → cg05280698 hypermethylation on higher kidney function, while CTDNEP1 mediates the impact of air pollution → cg03186999 hypomethylation on higher systolic blood pressure. This study investigates CpG sites and genes that mediate the impact of environmental chemicals on cardiometabolic traits. Furthermore, the multiomics approach described in this study provides a convenient workflow with which to investigate the impact of an environmental factor on the body's biomarkers, and, consequently, on health conditions, using publicly available data.

Immune gene regulation is associated with age and environmental adversity in a nonhuman primate

Phenotypic aging is ubiquitous across mammalian species, suggesting shared underlying mechanisms of aging. Aging is linked to molecular changes to DNA methylation and gene expression, and environmental factors, such as severe external challenges or adversities, can moderate these age-related changes. Yet, it remains unclear whether environmental adversities affect gene regulation via the same molecular pathways as chronological, or 'primary', aging. Investigating molecular aging in naturalistic animal populations can fill this gap by providing insight into shared molecular mechanisms of aging and the effects of a greater diversity of environmental adversities - particularly those that can be challenging to study in humans or laboratory organisms. Here, we characterised molecular aging - specifically, CpG methylation - in a sample of free-ranging rhesus macaques living off the coast of Puerto Rico (n samples = 571, n individuals = 499), which endured a major hurricane during our study. Age was associated with methylation at 78,661 sites (31% of all sites tested). Age-associated hypermethylation occurred more frequently in areas of active gene regulation, while hypomethylation was enriched in regions that show less activity in immune cells, suggesting these regions may become de-repressed in older individuals. Age-associated hypomethylation also co-occurred with increased chromatin accessibility while hypermethylation showed the opposite trend, hinting at a coordinated, multi-level loss of epigenetic stability during aging. We detected 32,048 CpG sites significantly associated with exposure to a hurricane, and these sites overlapped age-associated sites, most strongly in regulatory regions and most weakly in quiescent regions. Together, our results suggest that environmental adversity may contribute to aging-related molecular phenotypes in regions of active gene transcription, but that primary aging has specific signatures in non-regulatory regions.

MIAMI-AD (Methylation in Aging and Methylation in AD): an integrative knowledgebase that facilitates explorations of DNA methylation across sex, aging, and Alzheimer's disease

Alzheimer's disease (AD) is a common neurodegenerative disorder with a significant impact on aging populations. DNA methylation (DNAm) alterations have been implicated in both the aging processes and the development of AD. Given that AD affects more women than men, it is also important to explore DNAm changes that occur specifically in each sex. We created MIAMI-AD, a comprehensive knowledgebase containing manually curated summary statistics from 98 published tables in 38 studies, all of which included at least 100 participants. MIAMI-AD enables easy browsing, querying, and downloading DNAm associations at multiple levels-at individual CpG, gene, genomic regions, or genome-wide, in one or multiple studies. Moreover, it also offers tools to perform integrative analyses, such as comparing DNAm associations across different phenotypes or tissues, as well as interactive visualizations. Using several use case examples, we demonstrated that MIAMI-AD facilitates our understanding of age-associated CpGs in AD and the sex-specific roles of DNAm in AD. This open-access resource is freely available to the research community, and all the underlying data can be downloaded. MIAMI-AD facilitates integrative explorations to better understand the interplay between DNAm across aging, sex, and AD. Database URL: https://miami-ad.org/.

Comparison of Infinium MethylationEPIC v2.0 to v1.0 for human population epigenetics: considerations for addressing EPIC version differences in DNA methylation-based tools

Background

The recently launched DNA methylation profiling platform, Illumina MethylationEPIC BeadChip Infinium microarray v2.0 (EPICv2), is highly correlated with measurements obtained from its predecessor MethylationEPIC BeadChip Infinium microarray v1.0 (EPICv1). However, the concordance between the two versions in the context of DNA methylation-based tools, including cell type deconvolution algorithms, epigenetic clocks, and inflammation and lifestyle biomarkers has not yet been investigated.

Findings

We profiled DNA methylation on both EPIC versions using matched venous blood samples from individuals spanning early to late adulthood across three cohorts. On combining the DNA methylomes of the cohorts, we observed that samples primarily clustered by the EPIC version they were measured on. Within each cohort, when we calculated cell type proportions, epigenetic age acceleration (EAA), rate of aging estimates, and biomarker scores for the matched samples on each version, we noted significant differences between EPICv1 and EPICv2 in the majority of these estimates. These differences were not significant, however, when estimates were adjusted for EPIC version or when EAAs were calculated separately for each EPIC version.

Conclusions

Our findings indicate that EPIC version differences predominantly explain DNA methylation variation and influence estimates of DNA methylation-based tools, and therefore we recommend caution when combining cohorts run on different versions. We demonstrate the importance of calculating DNA methylation-based estimates separately for each EPIC version or accounting for EPIC version either as a covariate in statistical models or by using version correction algorithms.

mLiftOver: harmonizing data across Infinium DNA methylation platforms

MOTIVATION

Infinium DNA methylation BeadChips are widely used for genome-wide DNA methylation profiling at the population scale. Recent updates to probe content and naming conventions in the EPIC version 2 (EPICv2) arrays have complicated integrating new data with previous Infinium array platforms, such as the MethylationEPIC (EPIC) and the HumanMethylation450 (HM450) BeadChip.

RESULTS

We present mLiftOver, a user-friendly tool that harmonizes probe ID, methylation level, and signal intensity data across different Infinium platforms. It manages probe replicates, missing data imputation, and platform-specific bias for accurate data conversion. We validated the tool by applying HM450-based cancer classifiers to EPICv2 cancer data, achieving high accuracy. Additionally, we successfully integrated EPICv2 healthy tissue data with legacy HM450 data for tissue identity analysis and produced consistent copy number profiles in cancer cells.

AVAILABILITY AND IMPLEMENTATION

mLiftOver is implemented R and available in the Bioconductor package SeSAMe (version 1.21.13+): https://bioconductor.org/packages/release/bioc/html/sesame.html. Analysis of EPIC and EPICv2 platform-specific bias and high-confidence mapping is available at https://github.com/zhou-lab/InfiniumAnnotationV1/raw/main/Anno/EPICv2/EPICv2ToEPIC_conversion.tsv.gz. The source code is available at https://github.com/zwdzwd/sesame/blob/devel/R/mLiftOver.R under the MIT license.

Development and Validation of a Novel Placental DNA Methylation Biomarker of Maternal Smoking during Pregnancy in the ECHO Program

BACKGROUND

Maternal cigarette smoking during pregnancy (MSDP) is associated with numerous adverse health outcomes in infants and children with potential lifelong consequences. Negative effects of MSDP on placental DNA methylation (DNAm), placental structure, and function are well established.

OBJECTIVE

Our aim was to develop biomarkers of MSDP using DNAm measured in placentas (N = 96 ), collected as part of the Vitamin C to Decrease the Effects of Smoking in Pregnancy on Infant Lung Function double-blind, placebo-controlled randomized clinical trial conducted between 2012 and 2016. We also aimed to develop a digital polymerase chain reaction (PCR) assay for the top ranking cytosine-guanine dinucleotide (CpG) so that large numbers of samples can be screened for exposure at low cost.

METHODS

We compared the ability of four machine learning methods [logistic least absolute shrinkage and selection operator (LASSO) regression, logistic elastic net regression, random forest, and gradient boosting machine] to classify MSDP based on placental DNAm signatures. We developed separate models using the complete EPIC array dataset and on the subset of probes also found on the 450K array so that models exist for both platforms. For comparison, we developed a model using CpGs previously associated with MSDP in placenta. For each final model, we used model coefficients and normalized beta values to calculate placental smoking index (PSI) scores for each sample. Final models were validated in two external datasets: the Extremely Low Gestational Age Newborn observational study, N = 426 ; and the Rhode Island Children's Health Study, N = 237 .

RESULTS

Logistic LASSO regression demonstrated the highest performance in cross-validation testing with the lowest number of input CpGs. Accuracy was greatest in external datasets when using models developed for the same platform. PSI scores in smokers only (n = 72 ) were moderately correlated with maternal plasma cotinine levels. One CpG (cg27402634), with the largest coefficient in two models, was measured accurately by digital PCR compared with measurement by EPIC array (R 2 = 0.98 ).

DISCUSSION

To our knowledge, we have developed the first placental DNAm-based biomarkers of MSDP with broad utility to studies of prenatal disease origins. https://doi.org/10.1289/EHP13838.

Longitudinal analysis of epigenome-wide DNA methylation reveals novel loci associated with BMI change in East Asians

BACKGROUND

Obesity is a global public health concern linked to chronic diseases such as cardiovascular disease and type 2 diabetes (T2D). Emerging evidence suggests that epigenetic modifications, particularly DNA methylation, may contribute to obesity. However, the molecular mechanism underlying the longitudinal change of BMI has not been well-explored, especially in East Asian populations.

METHODS

This study performed a longitudinal epigenome-wide association analysis of DNA methylation to uncover novel loci associated with BMI change in 533 individuals across two Chinese cohorts with repeated DNA methylation and BMI measurements over four years.

RESULTS

We identified three novel CpG sites (cg14671384, cg25540824, and cg10848724) significantly associated with BMI change. Two of the identified CpG sites were located in regions previously associated with body shape and basal metabolic rate. Annotation of the top 20 BMI change-associated CpGs revealed strong connections to obesity and T2D. Notably, these CpGs exhibited active regulatory roles and located in genes with high expression in the liver and digestive tract, suggesting a potential regulatory pathway from genome to phenotypes of energy metabolism and absorption via DNA methylation. Cross-sectional and longitudinal EWAS comparisons indicated different mechanisms between CpGs related to BMI and BMI change.

CONCLUSION

This study enhances our understanding of the epigenetic dynamics underlying BMI change and emphasizes the value of longitudinal analyses in deciphering the complex interplay between epigenetics and obesity.

MSA: scalable DNA methylation screening BeadChip for high-throughput trait association studies

The Infinium DNA Methylation BeadChips have significantly contributed to population-scale epigenetics research by enabling epigenome-wide trait association discoveries. Here, we design, describe, and experimentally verify a new iteration of this technology, the Methylation Screening Array (MSA), to focus on human trait screening and discovery. This array utilizes extensive data from previous Infinium platform-based epigenome-wide association studies (EWAS). It incorporates knowledge from the latest single-cell and cell type-resolution whole genome methylome profiles. The MSA is engineered to achieve scalable screening of epigenetics-trait association in an ultra-high sample throughput. Our design encompassed diverse human trait associations, including those with genetic, cellular, environmental, and demographical variables and human diseases such as genetic, neurodegenerative, cardiovascular, infectious, and immune diseases. We comprehensively evaluated this array's reproducibility, accuracy, and capacity for cell-type deconvolution and supporting 5-hydroxymethylation profiling in diverse human tissues. Our first atlas data using this platform uncovered the complex chromatin and tissue contexts of DNA modification variations and genetic variants linked to human phenotypes.

DNA methylation signatures of youth-onset type 2 diabetes and exposure to maternal diabetes

OBJECTIVE

Youth-onset type 2 diabetes (T2D) is physiologically distinct from adult-onset, but it is not clear how the two diseases differ at a molecular level. In utero exposure to maternal type 2 diabetes (T2D) is known to be a specific risk factor for youth-onset T2D. DNA methylation (DNAm) changes associated with T2D but which differ between youth- and adult-onset might delineate the impacts of T2D development at different ages and could also determine the contribution of exposure to in utero diabetes.

METHODS

We performed an epigenome-wide analysis of DNAm on whole blood from 218 youth with T2D and 77 normoglycemic controls from the iCARE (improving renal Complications in Adolescents with type 2 diabetes through REsearch) cohort. Associations were tested using multiple linear regression models while adjusting for maternal diabetes, sex, age, BMI, smoking status, second-hand smoking exposure, cell-type proportions and genetic ancestry.

RESULTS

We identified 3830 differentially methylated sites associated with youth T2D onset, of which 3794 were moderately (adjusted p-value < 0.05 and effect size estimate > 0.01) associated and 36 were strongly (adjusted p-value < 0.05 and effect size estimate > 0.05) associated. A total of 3725 of these sites were not previously reported in the EWAS Atlas as associated with T2D, adult obesity or youth obesity. Moreover, three CpGs associated with youth-onset T2D in the PFKFB3 gene were also associated with maternal T2D exposure (FDR < 0.05 and effect size > 0.01). This is the first study to link PFKFB3 and T2D in youth.

CONCLUSION

Our findings support that T2D in youth has different impacts on DNAm than adult-onset, and suggests that changes in DNAm could provide an important link between in utero exposure to maternal diabetes and the onset of T2D.

Analysis of blood methylation quantitative trait loci in East Asians reveals ancestry-specific impacts on complex traits

Methylation quantitative trait loci (mQTLs) are essential for understanding the role of DNA methylation changes in genetic predisposition, yet they have not been fully characterized in East Asians (EAs). Here we identified mQTLs in whole blood from 3,523 Chinese individuals and replicated them in additional 1,858 Chinese individuals from two cohorts. Over 9% of mQTLs displayed specificity to EAs, facilitating the fine-mapping of EA-specific genetic associations, as shown for variants associated with height. Trans-mQTL hotspots revealed biological pathways contributing to EA-specific genetic associations, including an ERG-mediated 233 trans-mCpG network, implicated in hematopoietic cell differentiation, which likely reflects binding efficiency modulation of the ERG protein complex. More than 90% of mQTLs were shared between different blood cell lineages, with a smaller fraction of lineage-specific mQTLs displaying preferential hypomethylation in the respective lineages. Our study provides new insights into the mQTL landscape across genetic ancestries and their downstream effects on cellular processes and diseases/traits.

Low-input and single-cell methods for Infinium DNA methylation BeadChips

The Infinium BeadChip is the most widely used DNA methylome assay technology for population-scale epigenome profiling. However, the standard workflow requires over 200 ng of input DNA, hindering its application to small cell-number samples, such as primordial germ cells. We developed experimental and analysis workflows to extend this technology to suboptimal input DNA conditions, including ultra-low input down to single cells. DNA preamplification significantly enhanced detection rates to over 50% in five-cell samples and ∼25% in single cells. Enzymatic conversion also substantially improved data quality. Computationally, we developed a method to model the background signal's influence on the DNA methylation level readings. The modified detection P-value calculation achieved higher sensitivities for low-input datasets and was validated in over 100 000 public diverse methylome profiles. We employed the optimized workflow to query the demethylation dynamics in mouse primordial germ cells available at low cell numbers. Our data revealed nuanced chromatin states, sex disparities, and the role of DNA methylation in transposable element regulation during germ cell development. Collectively, we present comprehensive experimental and computational solutions to extend this widely used methylation assay technology to applications with limited DNA.

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.

Epigenomic insights into common human disease pathology

The epigenome-the chemical modifications and chromatin-related packaging of the genome-enables the same genetic template to be activated or repressed in different cellular settings. This multi-layered mechanism facilitates cell-type specific function by setting the local sequence and 3D interactive activity level. Gene transcription is further modulated through the interplay with transcription factors and co-regulators. The human body requires this epigenomic apparatus to be precisely installed throughout development and then adequately maintained during the lifespan. The causal role of the epigenome in human pathology, beyond imprinting disorders and specific tumour suppressor genes, was further brought into the spotlight by large-scale sequencing projects identifying that mutations in epigenomic machinery genes could be critical drivers in both cancer and developmental disorders. Abrogation of this cellular mechanism is providing new molecular insights into pathogenesis. However, deciphering the full breadth and implications of these epigenomic changes remains challenging. Knowledge is accruing regarding disease mechanisms and clinical biomarkers, through pathogenically relevant and surrogate tissue analyses, respectively. Advances include consortia generated cell-type specific reference epigenomes, high-throughput DNA methylome association studies, as well as insights into ageing-related diseases from biological 'clocks' constructed by machine learning algorithms. Also, 3rd-generation sequencing is beginning to disentangle the complexity of genetic and DNA modification haplotypes. Cell-free DNA methylation as a cancer biomarker has clear clinical utility and further potential to assess organ damage across many disorders. Finally, molecular understanding of disease aetiology brings with it the opportunity for exact therapeutic alteration of the epigenome through CRISPR-activation or inhibition.

Analytical challenges in omics research on asthma and allergy: A National Institute of Allergy and Infectious Diseases workshop

Studies of asthma and allergy are generating increasing volumes of omics data for analysis and interpretation. The National Institute of Allergy and Infectious Diseases (NIAID) assembled a workshop comprising investigators studying asthma and allergic diseases using omics approaches, omics investigators from outside the field, and NIAID medical and scientific officers to discuss the following areas in asthma and allergy research: genomics, epigenomics, transcriptomics, microbiomics, metabolomics, proteomics, lipidomics, integrative omics, systems biology, and causal inference. Current states of the art, present challenges, novel and emerging strategies, and priorities for progress were presented and discussed for each area. This workshop report summarizes the major points and conclusions from this NIAID workshop. As a group, the investigators underscored the imperatives for rigorous analytic frameworks, integration of different omics data types, cross-disciplinary interaction, strategies for overcoming current limitations, and the overarching goal to improve scientific understanding and care of asthma and allergic diseases.

DNA methylation marker identification and poly-methylation risk score in prediction of healthspan termination

Aim: To elucidate the epigenetic consequences of DNA methylation in healthspan termination (HST), considering the current limited understanding. Materials & methods: Genetically predicted DNA methylation models were established (n = 2478). These models were applied to genome-wide association study data on HST. Then, a poly-methylation risk score (PMRS) was established in 241,008 individuals from the UK Biobank. Results: Of the 63,046 CpGs from the prediction models, 13 novel CpGs were associated with HST. Furthermore, people with high PMRSs showed higher HST risk (hazard ratio: 1.18; 95% CI: 1.13-1.25). Conclusion: The study indicates that DNA methylation may influence HST by regulating the expression of genes (e.g., PRMT6, CTSK). PMRSs have a promising application in discriminating subpopulations to facilitate early prevention.

Associations of measured and genetically predicted leukocyte telomere length with vascular phenotypes: a population-based study

Shorter leukocyte telomere length (LTL) is associated with cardiovascular dysfunction. Whether this association differs between measured and genetically predicted LTL is still unclear. Moreover, the molecular processes underlying the association remain largely unknown. We used baseline data of the Rhineland Study, an ongoing population-based cohort study in Bonn, Germany [56.2% women, age: 55.5 ± 14.0 years (range 30 - 95 years)]. We calculated genetically predicted LTL in 4180 participants and measured LTL in a subset of 1828 participants with qPCR. Using multivariable regression, we examined the association of measured and genetically predicted LTL, and the difference between measured and genetically predicted LTL (ΔLTL), with four vascular functional domains and the overall vascular health. Moreover, we performed epigenome-wide association studies of three LTL measures. Longer measured LTL was associated with better microvascular and cardiac function. Longer predicted LTL was associated with better cardiac function. Larger ΔLTL was associated with better microvascular and cardiac function and overall vascular health, independent of genetically predicted LTL. Several CpGs were associated (p < 1e-05) with measured LTL (n = 5), genetically predicted LTL (n = 8), and ΔLTL (n = 27). Genes whose methylation status was associated with ΔLTL were enriched in vascular endothelial signaling pathways and have been linked to environmental exposures, cardiovascular diseases, and mortality. Our findings suggest that non-genetic causes of LTL contribute to microvascular and cardiac function and overall vascular health, through an effect on the vascular endothelial signaling pathway. Interventions that counteract LTL may thus improve vascular function.

Identification of genetically predicted DNA methylation markers associated with non-small cell lung cancer risk among 34,964 cases and 448,579 controls

BACKGROUND

Although the associations between genetic variations and lung cancer risk have been explored, the epigenetic consequences of DNA methylation in lung cancer development are largely unknown. Here, the genetically predicted DNA methylation markers associated with non-small cell lung cancer (NSCLC) risk by a two-stage case-control design were investigated.

METHODS

The genetic prediction models for methylation levels based on genetic and methylation data of 1595 subjects from the Framingham Heart Study were established. The prediction models were applied to a fixed-effect meta-analysis of screening data sets with 27,120 NSCLC cases and 27,355 controls to identify the methylation markers, which were then replicated in independent data sets with 7844 lung cancer cases and 421,224 controls. Also performed was a multi-omics functional annotation for the identified CpGs by integrating genomics, epigenomics, and transcriptomics and investigation of the potential regulation pathways.

RESULTS

Of the 29,894 CpG sites passing the quality control, 39 CpGs associated with NSCLC risk (Bonferroni-corrected p ≤ 1.67 × 10-6 ) were originally identified. Of these, 16 CpGs remained significant in the validation stage (Bonferroni-corrected p ≤ 1.28 × 10-3 ), including four novel CpGs. Multi-omics functional annotation showed nine of 16 CpGs were potentially functional biomarkers for NSCLC risk. Thirty-five genes within a 1-Mb window of 12 CpGs that might be involved in regulatory pathways of NSCLC risk were identified.

CONCLUSIONS

Sixteen promising DNA methylation markers associated with NSCLC were identified. Changes of the methylation level at these CpGs might influence the development of NSCLC by regulating the expression of genes nearby.

PLAIN LANGUAGE SUMMARY

The epigenetic consequences of DNA methylation in lung cancer development are still largely unknown. This study used summary data of large-scale genome-wide association studies to investigate the associations between genetically predicted levels of methylation biomarkers and non-small cell lung cancer risk at the first time. This study looked at how well larotrectinib worked in adult patients with sarcomas caused by TRK fusion proteins. These findings will provide a unique insight into the epigenetic susceptibility mechanisms of lung cancer.

Epigenome-wide association study of total nicotine equivalents in multiethnic current smokers from three prospective cohorts

The impact of tobacco exposure on health varies by race and ethnicity and is closely tied to internal nicotine dose, a marker of carcinogen uptake. DNA methylation is strongly responsive to smoking status and may mediate health effects, but study of associations with internal dose is limited. We performed a blood leukocyte epigenome-wide association study (EWAS) of urinary total nicotine equivalents (TNEs; a measure of nicotine uptake) and DNA methylation measured using the MethylationEPIC v1.0 BeadChip (EPIC) in six racial and ethnic groups across three cohort studies. In the Multiethnic Cohort Study (discovery, n = 1994), TNEs were associated with differential methylation at 408 CpG sites across >250 genomic regions (p < 9 × 10-8). The top significant sites were annotated to AHRR, F2RL3, RARA, GPR15, PRSS23, and 2q37.1, all of which had decreasing methylation with increasing TNEs. We identified 45 novel CpG sites, of which 42 were unique to the EPIC array and eight annotated to genes not previously linked with smoking-related DNA methylation. The most significant signal in a novel gene was cg03748458 in MIR383;SGCZ. Fifty-one of the 408 discovery sites were validated in the Singapore Chinese Health Study (n = 340) and the Southern Community Cohort Study (n = 394) (Bonferroni corrected p < 1.23 × 10-4). Significant heterogeneity by race and ethnicity was detected for CpG sites in MYO1G and CYTH1. Furthermore, TNEs significantly mediated the association between cigarettes per day and DNA methylation at 15 sites (average 22.5%-44.3% proportion mediated). Our multiethnic study highlights the transethnic and ethnic-specific methylation associations with internal nicotine dose, a strong predictor of smoking-related morbidities.

DNA methylation variation after a parenting program for child conduct problems: Findings from a randomized controlled trial

This study investigated associations of the Incredible Years (IY) parenting program with children's DNA methylation. Participants were 289 Dutch children aged 3-9 years (75% European ancestry, 48% female) with above-average conduct problems. Saliva was collected 2.5 years after families were randomized to IY or care as usual (CAU). Using an intention-to-treat approach, confirmatory multiple-regression analyses revealed no significant differences between the IY and CAU groups in children's methylation levels at the NR3C1 and FKBP5 genes. However, exploratory epigenome-wide analyses revealed nine differentially methylated regions between groups, coinciding with SLAMF1, MITF, FAM200B, PSD3, SNX31, and CELSR1. The study provides preliminary evidence for associations of IY with children's salivary methylation levels and highlights the need for further research into biological outcomes of parenting programs.

Epigenome-wide association study of dietary fatty acid intake

BACKGROUND

Dietary intake of n-3 polyunsaturated fatty acids (PUFA) may have a protective effect on the development of cardiovascular diseases, diabetes, depression and cancer, while a high intake of n-6 PUFA was often reported to be associated with inflammation-related traits. The effect of PUFAs on health outcomes might be mediated by DNA methylation (DNAm). The aim of our study is to identify the impact of PUFA intake on DNAm in the Cooperative Health Research in the Region of Augsburg (KORA) FF4 cohort and the Leiden Longevity Study (LLS).

RESULTS

DNA methylation levels were measured in whole blood from the population-based KORA FF4 study (N = 1354) and LLS (N = 448), using the Illumina MethylationEPIC BeadChip and Illumina HumanMethylation450 array, respectively. We assessed associations between DNAm and intake of eight and four PUFAs in KORA and LLS, respectively. Where possible, results were meta-analyzed. Below the Bonferroni correction threshold (p < 7.17 × 10-8), we identified two differentially methylated positions (DMPs) associated with PUFA intake in the KORA study. The DMP cg19937480, annotated to gene PRDX1, was positively associated with docosahexaenoic acid (DHA) in model 1 (beta: 2.00 × 10-5, 95%CI: 1.28 × 10-5-2.73 × 10-5, P value: 6.98 × 10-8), while cg05041783, annotated to gene MARK2, was positively associated with docosapentaenoic acid (DPA) in our fully adjusted model (beta: 9.80 × 10-5, 95%CI: 6.25 × 10-5-1.33 × 10-4, P value: 6.75 × 10-8). In the meta-analysis, we identified the CpG site (cg15951061), annotated to gene CDCA7L below Bonferroni correction (1.23 × 10-7) associated with eicosapentaenoic acid (EPA) intake in model 1 (beta: 2.00 × 10-5, 95% CI: 1.27 × 10-5-2.73 × 10-5, P value = 5.99 × 10-8) and we confirmed the association of cg19937480 with DHA in both models 1 and 2 (beta: 2.07 × 10-5, 95% CI: 1.31 × 10-5-2.83 × 10-5, P value = 1.00 × 10-7 and beta: 2.19 × 10-5, 95% CI: 1.41 × 10-5-2.97 × 10-5, P value = 5.91 × 10-8 respectively).

CONCLUSIONS

Our study identified three CpG sites associated with PUFA intake. The mechanisms of these sites remain largely unexplored, highlighting the novelty of our findings. Further research is essential to understand the links between CpG site methylation and PUFA outcomes.

14-Week exercise training modifies the DNA methylation levels at gene sites in non-Alzheimer's disease women aged 50 to 70 years

Exercise training emerges as a key strategy in lifestyle modification, capable of reducing the risk of developing Alzheimer's disease (AD) due to risk factors such as age, family history, genetics and low level of education associated with AD. We aim to analyze the effect of a 14-week combined exercise training (CT) on the methylation of genes associated with AD in non-alzheimer's disease women. CT sessions lasted 60 min, occurring three times a week for 14 weeks. Forty non-Alzheimer's disease women aged 50 to 70 years (60.7 ± 4.1 years) with a mean height of 1.6 ± 0.1 m, mean weight of 73.12 ± 9.0 kg and a mean body mass index of 29.69 ± 3.5 kg/m2, underwent two physical assessments: pre and post the 14 weeks. DNA methylation assays utilized the EPIC Infinium Methylation BeadChip from Illumina. We observed that 14 weeks of CT led to reductions in systolic (p = 0.001) and diastolic (p = 0.017) blood pressure and improved motor skills post-intervention. Among 25 genes linked to AD, CT induced differentially methylated sites in 12 genes, predominantly showing hypomethylated sites (negative β values). Interestingly, despite hypomethylated sites, some genes exhibited hypermethylated sites (positive β values), such as ABCA7, BDNF, and WWOX. A 14-week CT regimen was adequate to induce differential methylation in 12 CE-related genes in healthy older women, alongside improvements in motor skills and blood pressure. In conclusion, this study suggest that combined training can be a strategy to improve physical fitness in older individuals, especially able to induce methylation alterations in genes sites related to development of AD. It is important to highlight that training should act as protective factor in older adults.

Novel markers and networks related to restored skeletal muscle transcriptome after bariatric surgery

OBJECTIVE

The aim of this study was to discover novel markers underlying the improvement of skeletal muscle metabolism after bariatric surgery.

METHODS

Skeletal muscle transcriptome data of lean people and people with obesity, before and 1 year after bariatric surgery, were subjected to weighted gene co-expression network analysis (WGCNA) and least absolute shrinkage and selection operator (LASSO) regression. Results of LASSO were confirmed in a replication cohort.

RESULTS

The expression levels of 440 genes differing between individuals with and without obesity were no longer different 1 year after surgery, indicating restoration. WGCNA clustered 116 genes with normalized expression in one major module, particularly correlating to weight loss and decreased plasma free fatty acids (FFA), 44 of which showed an obesity-related phenotype upon deletion in mice. Among the genes of the major module, 105 represented prominent markers for reduced FFA concentration, including 55 marker genes for decreased BMI in both the discovery and replication cohorts.

CONCLUSIONS

Previously unknown gene networks and marker genes underlined the important role of FFA in restoring muscle gene expression after bariatric surgery and further suggest novel therapeutic targets for obesity.

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.

A Pregnancy and Childhood Epigenetics Consortium (PACE) meta-analysis highlights potential relationships between birth order and neonatal blood DNA methylation

Higher birth order is associated with altered risk of many disease states. Changes in placentation and exposures to in utero growth factors with successive pregnancies may impact later life disease risk via persistent DNA methylation alterations. We investigated birth order with Illumina DNA methylation array data in each of 16 birth cohorts (8164 newborns) with European, African, and Latino ancestries from the Pregnancy and Childhood Epigenetics Consortium. Meta-analyzed data demonstrated systematic DNA methylation variation in 341 CpGs (FDR adjusted P < 0.05) and 1107 regions. Forty CpGs were located within known quantitative trait loci for gene expression traits in blood, and trait enrichment analysis suggested a strong association with immune-related, transcriptional control, and blood pressure regulation phenotypes. Decreasing fertility rates worldwide with the concomitant increased proportion of first-born children highlights a potential reflection of birth order-related epigenomic states on changing disease incidence trends.

Database Resources of the National Genomics Data Center, China National Center for Bioinformation in 2024

The National Genomics Data Center (NGDC), which is a part of the China National Center for Bioinformation (CNCB), provides a family of database resources to support the global academic and industrial communities. With the rapid accumulation of multi-omics data at an unprecedented pace, CNCB-NGDC continuously expands and updates core database resources through big data archiving, integrative analysis and value-added curation. Importantly, NGDC collaborates closely with major international databases and initiatives to ensure seamless data exchange and interoperability. Over the past year, significant efforts have been dedicated to integrating diverse omics data, synthesizing expanding knowledge, developing new resources, and upgrading major existing resources. Particularly, several database resources are newly developed for the biodiversity of protists (P10K), bacteria (NTM-DB, MPA) as well as plant (PPGR, SoyOmics, PlantPan) and disease/trait association (CROST, HervD Atlas, HALL, MACdb, BioKA, BioKA, RePoS, PGG.SV, NAFLDkb). All the resources and services are publicly accessible at https://ngdc.cncb.ac.cn.

Trans-ancestry epigenome-wide association meta-analysis of DNA methylation with lifetime cannabis use

Cannabis is widely used worldwide, yet its links to health outcomes are not fully understood. DNA methylation can serve as a mediator to link environmental exposures to health outcomes. We conducted an epigenome-wide association study (EWAS) of peripheral blood-based DNA methylation and lifetime cannabis use (ever vs. never) in a meta-analysis including 9436 participants (7795 European and 1641 African ancestry) from seven cohorts. Accounting for effects of cigarette smoking, our trans-ancestry EWAS meta-analysis revealed four CpG sites significantly associated with lifetime cannabis use at a false discovery rate of 0.05 ( p < 5.85 × 10 - 7 ) : cg22572071 near gene ADGRF1, cg15280358 in ADAM12, cg00813162 in ACTN1, and cg01101459 near LINC01132. Additionally, our EWAS analysis in participants who never smoked cigarettes identified another epigenome-wide significant CpG site, cg14237301 annotated to APOBR. We used a leave-one-out approach to evaluate methylation scores constructed as a weighted sum of the significant CpGs. The best model can explain 3.79% of the variance in lifetime cannabis use. These findings unravel the DNA methylation changes associated with lifetime cannabis use that are independent of cigarette smoking and may serve as a starting point for further research on the mechanisms through which cannabis exposure impacts health outcomes.

Longitudinal changes and variation in human DNA methylation analysed with the Illumina MethylationEPIC BeadChip assay and their implications on forensic age prediction

DNA methylation, a pivotal epigenetic modification, plays a crucial role in regulating gene expression and is known to undergo dynamic changes with age. The present study investigated epigenome-wide methylation profiles in 64 individuals over two time points, 15 years apart, using the Illumina EPIC850k arrays. A mixed-effects model identified 2821 age-associated differentially methylated CpG positions (aDMPs) with a median rate of change of 0.18% per year, consistent with a 10-15% change during a human lifespan. Significant variation in the baseline DNA methylation levels between individuals of similar ages as well as inconsistent direction of change with time across individuals were observed for all the aDMPs. Twenty-three of the 2821 aDMPs were previously incorporated into forensic age prediction models. These markers displayed larger changes in DNA methylation with age compared to all the aDMPs and less variation among individuals. Nevertheless, the forensic aDMPs also showed inter-individual variations in the direction of DNA methylation changes. Only cg16867657 in ELOVL2 exhibited a uniform direction of the age-related change among the investigated individuals, which supports the current knowledge that CpG sites in ELOVL2 are the best markers for age prediction.

MIAMI-AD (Methylation in Aging and Methylation in AD): an integrative knowledgebase that facilitates explorations of DNA methylation across sex, aging, and Alzheimer's disease

Alzheimer's disease (AD) is a common neurodegenerative disorder with a significant impact on aging populations. DNA methylation (DNAm) alterations have been implicated in both the aging processes and the development of AD. Given that AD affects more women than men, it is also important to explore DNAm changes that occur specifically in each sex. We created MIAMI-AD, a comprehensive knowledge base containing manually curated summary statistics from 97 published tables in 37 studies, all of which included at least 100 participants. MIAMI-AD enables easy browsing, querying, and downloading DNAm associations at multiple levels - at individual CpG, gene, genomic regions, or genome-wide, in one or multiple studies. Moreover, it also offers tools to perform integrative analyses, such as comparing DNAm associations across different phenotypes or tissues, as well as interactive visualizations. Using several use case examples, we demonstrated that MIAMI-AD facilitates our understanding of age-associated CpGs in AD and the sex-specific roles of DNAm in AD. This open-access resource is freely available to the research community, and all the underlying data can be downloaded. MIAMI-AD (https://miami-ad.org/) facilitates integrative explorations to better understand the interplay between DNAm across aging, sex, and AD.

Waterpipe and cigarette epigenome analysis reveals markers implicated in addiction and smoking type inference

Waterpipe smoking is frequent in the Middle East and Africa with emerging prevalence worldwide. The epigenome acts as a molecular sensor to exposures and a crucial driver in several diseases. With the widespread use of waterpipe smoking, it is timely to investigate its epigenomic markers and their role in addiction, as a central player in disease prevention and therapeutic strategies. DNA methylome-wide profiling was performed on an exposure-rich population from the Middle East, constituting of 216 blood samples split equally between never, cigarette-only and waterpipe-only smokers. Waterpipe smokers showed predominantly distinct epigenetic markers from cigarette smokers, even though both smoking forms are tobacco-based. Moreover, each smoking form could be accurately (∼90 %) inferred from the DNA methylome using machine learning. Top markers showed dose-response relationship with extent of smoking and were validated using independent technologies and additional samples (total N = 284). Smoking markers were enriched in regulatory regions and several biological pathways, primarily addiction. The epigenetically altered genes were not associated with genetic etiology of tobacco use, and the methylation levels of addiction genes, in particular, were more likely to reverse after smoking cessation. In contrast, other epigenetic markers continued to feature smoking exposure after cessation, which may explain long-term health effects observed in former smokers. This study reports, for the first time, blood epigenome-wide markers of waterpipe smokers and reveals new markers of cigarette smoking, with implications in mechanisms of addiction and the capacity to discriminate between different smoking types. These markers may offer actionable targets to reverse the epigenetic memory of addiction and can guide future prevention strategies for tobacco smoking as the most preventable cause of illnesses worldwide.

Pan-precancer and cancer DNA methylation profiles revealed significant tissue specificity of interrupted biological processes in tumorigenesis

DNA methylation (DNAme) alterations are known to initiate from the precancerous stage of tumorigenesis. Herein, we investigated the global and local patterns of DNAme perturbations in tumorigenesis by analysing the genome-wide DNAme profiles of the cervix, colorectum, stomach, prostate, and liver at precancerous and cancer stages. We observed global hypomethylation in tissues of both two stages, except for the cervix, whose global DNAme level in normal tissue was lower than that of the other four tumour types. For alterations shared by both stages, there were common hyper-methylation (sHyperMethyl) and hypo-methylation (sHypoMethyl) changes, of which the latter type was more frequently identified in all tissues. Biological pathways interrupted by sHyperMethyl and sHypoMethyl alterations demonstrated significant tissue specificity. DNAme bidirectional chaos indicated by the enrichment of both sHyperMethyl and sHypoMethyl changes in the same pathway was observed in most tissues and was a common phenomenon, particularly in liver lesions. Moreover, for the same enriched pathways, different tissues may be affected by distinct DNAme types. For the PI3K-Akt signalling pathway, sHyperMethyl enrichment was observed in the prostate dataset, but sHypoMethyl enrichment was observed in the colorectum and liver datasets. Nevertheless, they did not show an increased possibility in survival prediction of patients in comparison with other DNAme types. Additionally, our study demonstrated that gene-body DNAme changes of tumour suppressor genes and oncogenes may persist from precancerous lesions to the tumour. Overall, we demonstrate the tissue specificity and commonality of cross-stage alterations in DNA methylation profiles in multi-tissue tumorigenesis.

Novel insights into the whole-blood DNA methylome of asthma in ethnically diverse children and youth

BACKGROUND

The epigenetic mechanisms of asthma remain largely understudied in African Americans and Hispanics/Latinos, two populations disproportionately affected by asthma. We aimed to identify markers, regions and processes with differential patterns of DNA methylation (DNAm) in whole blood by asthma status in ethnically diverse children and youth, and to assess their functional consequences.

METHODS

DNAm levels were profiled with the Infinium MethylationEPIC or HumanMethylation450 BeadChip arrays among 1226 African Americans or Hispanics/Latinos and assessed for differential methylation per asthma status at the CpG and region (differentially methylated region (DMR)) level. Novel associations were validated in blood and/or nasal epithelium from ethnically diverse children and youth. The functional and biological implications of the markers identified were investigated by combining epigenomics with transcriptomics from study participants.

RESULTS

128 CpGs and 196 DMRs were differentially methylated after multiple testing corrections, including 92.3% and 92.8% novel associations, respectively. 41 CpGs were replicated in other Hispanics/Latinos, prioritising cg17647904 (NCOR2) and cg16412914 (AXIN1) as asthma DNAm markers. Significant DNAm markers were enriched in previous associations for asthma, fractional exhaled nitric oxide, bacterial infections, immune regulation or eosinophilia. Functional annotation highlighted epigenetically regulated gene networks involved in corticosteroid response, host defence and immune regulation. Several implicated genes are targets for approved or experimental drugs, including TNNC1 and NDUFA12. Many differentially methylated loci previously associated with asthma were validated in our study.

CONCLUSIONS

We report novel whole-blood DNAm markers for asthma underlying key processes of the disease pathophysiology and confirm the transferability of previous asthma DNAm associations to ethnically diverse populations.

The methylome of buccal epithelial cells is influenced by age, sex, and physiological properties

Epigenetic modifications, particularly DNA methylation, have emerged as regulators of gene expression and are implicated in various biological processes and disease states. Understanding the factors influencing the epigenome is essential for unraveling its complexity. In this study, we aimed to identify how the methylome of buccal epithelial cells, a noninvasive and easily accessible tissue, is associated with demographic and health-related variables commonly used in clinical settings, such as age, sex, blood immune composition, hemoglobin levels, and others. We developed a model to assess the association of multiple factors with the human methylome and identify the genomic loci significantly impacted by each trait. We demonstrated that DNA methylation variation is accurately modeled by several factors. We confirmed the well-known impact of age and sex and unveiled novel clinical factors associated with DNA methylation, such as blood neutrophils, hemoglobin, red blood cell distribution width, high-density lipoprotein cholesterol, and urea. Genomic regions significantly associated with these traits were enriched in relevant transcription factors, drugs, and diseases. Among our findings, we showed that neutrophil-impacted loci were involved in neutrophil functionality and maturation. Similarly, hemoglobin-influenced sites were associated with several diseases, including aplastic anemia, and the genomic loci affected by urea were related to congenital anomalies of the kidney and urinary tract. Our findings contribute to a better understanding of the human methylome plasticity and provide insights into novel factors shaping DNA methylation patterns, highlighting their potential clinical implications as biomarkers and the importance of considering these physiological traits in future medical epigenomic investigations.NEW & NOTEWORTHY We have developed a quantitative model to assess how the human methylome is associated with several factors and to identify the genomic loci significantly impacted by each trait. We reported novel health-related factors driving DNA methylation patterns and new site-specific regulations that further elucidate methylome dynamics. Our study contributes to a better understanding of the plasticity of the human methylome and unveils novel physiological traits with a potential role in future medical epigenomic investigations.

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.

Systems genetics approaches for understanding complex traits with relevance for human disease

Quantitative traits are often complex because of the contribution of many loci, with further complexity added by environmental factors. In medical research, systems genetics is a powerful approach for the study of complex traits, as it integrates intermediate phenotypes, such as RNA, protein, and metabolite levels, to understand molecular and physiological phenotypes linking discrete DNA sequence variation to complex clinical and physiological traits. The primary purpose of this review is to describe some of the resources and tools of systems genetics in humans and rodent models, so that researchers in many areas of biology and medicine can make use of the data.

Epigenome-wide methylation and progression to high-grade cervical intraepithelial neoplasia (CIN2+): a prospective cohort study in the United States

BACKGROUND

Methylation levels may be associated with and serve as markers to predict risk of progression of precancerous cervical lesions. We conducted an epigenome-wide association study (EWAS) of CpG methylation and progression to high-grade cervical intraepithelial neoplasia (CIN2 +) following an abnormal screening test.

METHODS

A prospective US cohort of 289 colposcopy patients with normal or CIN1 enrollment histology was assessed. Baseline cervical sample DNA was analyzed using Illumina HumanMethylation 450K (n = 76) or EPIC 850K (n = 213) arrays. Participants returned at provider-recommended intervals and were followed up to 5 years via medical records. We assessed continuous CpG M values for 9 cervical cancer-associated genes and time-to-progression to CIN2+. We estimated CpG-specific time-to-event ratios (TTER) and hazard ratios using adjusted, interval-censored Weibull accelerated failure time models. We also conducted an exploratory EWAS to identify novel CpGs with false discovery rate (FDR) < 0.05.

RESULTS

At enrollment, median age was 29.2 years; 64.0% were high-risk HPV-positive, and 54.3% were non-white. During follow-up (median 24.4 months), 15 participants progressed to CIN2+. Greater methylation levels were associated with a shorter time-to-CIN2+ for CADM1 cg03505501 (TTER = 0.28; 95%CI 0.12, 0.63; FDR = 0.03) and RARB Cluster 1 (TTER = 0.46; 95% CI 0.29, 0.71; FDR = 0.01). There was evidence of similar trends for DAPK1 cg14286732, PAX1 cg07213060, and PAX1 Cluster 1. The EWAS detected 336 novel progression-associated CpGs, including those located in CpG islands associated with genes FGF22, TOX, COL18A1, GPM6A, XAB2, TIMP2, GSPT1, NR4A2, and APBB1IP.

CONCLUSIONS

Using prospective time-to-event data, we detected associations between CADM1-, DAPK1-, PAX1-, and RARB-related CpGs and cervical disease progression, and we identified novel progression-associated CpGs.

IMPACT

Methylation levels at novel CpG sites may help identify individuals with ≤CIN1 histology at higher risk of progression to CIN2+ and inform risk-based cervical cancer screening guidelines.

Associations between epigenome-wide DNA methylation and height-related traits among Sub-Saharan Africans: the RODAM study

Human height and related traits are highly complex, and extensively research has shown that these traits are determined by both genetic and environmental factors. Such factors may partially affect these traits through epigenetic programing. Epigenetic programing is dynamic and plays an important role in controlling gene expression and cell differentiation during (early) development. DNA methylation (DNAm) is the most commonly studied epigenetic feature. In this study we conducted an epigenome-wide DNAm association analysis on height-related traits in a Sub-Saharan African population, in order to detect DNAm biomarkers across four height-related traits. DNAm profiles were acquired in whole blood samples of 704 Ghanaians, sourced from the Research on Obesity and Diabetes among African Migrants study, using the Illumina Infinium HumanMethylation450 BeadChip. Linear models were fitted to detect differentially methylated positions (DMPs) and regions (DMRs) associated with height, leg-to-height ratio (LHR), leg length, and sitting height. No epigenome-wide significant DMPs were recorded. However we did observe among our top DMPs five informative probes associated with the height-related traits: cg26905768 (leg length), cg13268132 (leg length), cg19776793 (height), cg23072383 (LHR), and cg24625894 (sitting height). All five DMPs are annotated to genes whose functions were linked to bone cell regulation and development. DMR analysis identified overlapping DMRs within the gene body of HLA-DPB1 gene, and the HOXA gene cluster. In this first epigenome-wide association studies of these traits, our findings suggest DNAm associations with height-related heights, and might influence development and maintenance of these traits. Further studies are needed to replicate our findings, and to elucidate the molecular mechanism underlying human height-related traits.

Comprehensive Evaluation of The Infinium Human MethylationEPIC v2 BeadChip

Infinium Methylation BeadChips are widely used to profile DNA cytosine modifications in large cohort studies for reasons of cost-effectiveness, accurate quantification, and user-friendly data analysis in characterizing these canonical epigenetic marks. In this work, we conducted a comprehensive evaluation of the updated Infinium MethylationEPIC v2 BeadChip (EPICv2). Our evaluation revealed that EPICv2 offers significant improvements over its predecessors, including expanded enhancer coverage, applicability to diverse ancestry groups, support for low-input DNA down to one nanogram, coverage of existing epigenetic clocks, cell type deconvolution panels, and human trait associations, while maintaining accuracy and reproducibility. Using EPICv2, we were able to identify epigenome and sequence signatures in cell line models of DNMT and SETD2 loss and/or hypomorphism. Furthermore, we provided probe-wise evaluation and annotation to facilitate the use of new features on this array for studying the interplay between somatic mutations and epigenetic landscape in cancer genomics. In conclusion, EPICv2 provides researchers with a valuable tool for studying epigenetic modifications and their role in development and disease.

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.

Low-input and single-cell methods for Infinium DNA methylation BeadChips

The Infinium BeadChip is the most widely used DNA methylome assay technology for population-scale epigenome profiling. However, the standard workflow requires over 200 ng of input DNA, hindering its application to small cell-number samples, such as primordial germ cells. We developed experimental and analysis workflows to extend this technology to suboptimal input DNA conditions, including ultra-low input down to single cells. DNA preamplification significantly enhanced detection rates to over 50% in five-cell samples and ∼25% in single cells. Enzymatic conversion also substantially improved data quality. Computationally, we developed a method to model the background signal's influence on the DNA methylation level readings. The modified detection p -values calculation achieved higher sensitivities for low-input datasets and was validated in over 100,000 public datasets with diverse methylation profiles. We employed the optimized workflow to query the demethylation dynamics in mouse primordial germ cells available at low cell numbers. Our data revealed nuanced chromatin states, sex disparities, and the role of DNA methylation in transposable element regulation during germ cell development. Collectively, we present comprehensive experimental and computational solutions to extend this widely used methylation assay technology to applications with limited DNA.

Prioritization of therapeutic targets for dyslipidemia using integrative multi-omics and multi-trait analysis

Drug targets with genetic support are several-fold more likely to succeed in clinical trials. We introduce a genetic-driven approach based on causal inferences that can inform drug target prioritization, repurposing, and adverse effects of using lipid-lowering agents. Given that a multi-trait approach increases the power to detect meaningful variants/genes, we conduct multi-omics and multi-trait analyses, followed by network connectivity investigations, and prioritize 30 potential therapeutic targets for dyslipidemia, including SORT1, PSRC1, CELSR2, PCSK9, HMGCR, APOB, GRN, HFE2, FJX1, C1QTNF1, and SLC5A8. 20% (6/30) of prioritized targets from our hypothesis-free drug target search are either approved or under investigation for dyslipidemia. The prioritized targets are 22-fold higher in likelihood of being approved or under investigation in clinical trials than genome-wide association study (GWAS)-curated targets. Our results demonstrate that the genetic-driven approach used in this study is a promising strategy for prioritizing targets while informing about the potential adverse effects and repurposing opportunities.

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.

Fathers' preconception smoking and offspring DNA methylation

BACKGROUND

Experimental studies suggest that exposures may impact respiratory health across generations via epigenetic changes transmitted specifically through male germ cells. Studies in humans are, however, limited. We aim to identify epigenetic marks in offspring associated with father's preconception smoking.

METHODS

We conducted epigenome-wide association studies (EWAS) in the RHINESSA cohort (7-50 years) on father's any preconception smoking (n = 875 offspring) and father's pubertal onset smoking < 15 years (n = 304), using Infinium MethylationEPIC Beadchip arrays, adjusting for offspring age, own smoking and maternal smoking. EWAS of maternal and offspring personal smoking were performed for comparison. Father's smoking-associated dmCpGs were checked in subpopulations of offspring who reported no personal smoking and no maternal smoking exposure.

RESULTS

Father's smoking commencing preconception was associated with methylation of blood DNA in offspring at two cytosine-phosphate-guanine sites (CpGs) (false discovery rate (FDR) < 0.05) in PRR5 and CENPP. Father's pubertal onset smoking was associated with 19 CpGs (FDR < 0.05) mapped to 14 genes (TLR9, DNTT, FAM53B, NCAPG2, PSTPIP2, MBIP, C2orf39, NTRK2, DNAJC14, CDO1, PRAP1, TPCN1, IRS1 and CSF1R). These differentially methylated sites were hypermethylated and associated with promoter regions capable of gene silencing. Some of these sites were associated with offspring outcomes in this cohort including ever-asthma (NTRK2), ever-wheezing (DNAJC14, TPCN1), weight (FAM53B, NTRK2) and BMI (FAM53B, NTRK2) (p < 0.05). Pathway analysis showed enrichment for gene ontology pathways including regulation of gene expression, inflammation and innate immune responses. Father's smoking-associated sites did not overlap with dmCpGs identified in EWAS of personal and maternal smoking (FDR < 0.05), and all sites remained significant (p < 0.05) in analyses of offspring with no personal smoking and no maternal smoking exposure.

CONCLUSION

Father's preconception smoking, particularly in puberty, is associated with offspring DNA methylation, providing evidence that epigenetic mechanisms may underlie epidemiological observations that pubertal paternal smoking increases risk of offspring asthma, low lung function and obesity.

Multiomics Data Analysis Identified CpG Sites That Mediate the Impact of Smoking on Cardiometabolic Traits

Understanding the epigenome paths through which smoking contributes to cardiometabolic traits is important for downstream applications. In this study, an SNP-based analytical pipeline was used to integrate several publicly available datasets in order to identify CpG sites that mediate the impact of smoking on cardiometabolic traits and to investigate the underlying molecular mechanisms. After applying stringent statistical criteria, 11 CpG sites were detected that showed significant association (p < 5 × 10-8) with cardiometabolic traits at both the discovery and replication stages. By integrating eQTL data, I found genes behind a number of these associations. cg05228408 was hypomethylated in smokers and contributed to higher blood pressure by lowering the expression of the CLCN6 gene. cg08639339 was hypermethylated in smokers and lowered the metabolic rate by increasing the expression of RAB29; furthermore, I noted TMEM120A mediated the impact of smoking-cg17325771 on LDL, and LTBP3 mediated the smoking-cg07029024 effect on heart rate. The pathway analysis identified processes through which the identified genes impact their traits. This study provides a list of CpG sites that mediates the impact of smoking on cardiometabolic traits and a framework to investigate the underlying molecular paths using publicly available data.