Enhanced resolution profiling in twins reveals differential methylation signatures of type 2 diabetes with links to its complications

BACKGROUND

Type 2 diabetes (T2D) susceptibility is influenced by genetic and environmental factors. Previous findings suggest DNA methylation as a potential mechanism in T2D pathogenesis and progression.

METHODS

We profiled DNA methylation in 248 blood samples from participants of European ancestry from 7 twin cohorts using a methylation sequencing platform targeting regulatory genomic regions encompassing 2,048,698 CpG sites.

FINDINGS

We find and replicate 3 previously unreported T2D differentially methylated CpG positions (T2D-DMPs) at FDR 5% in RGL3, NGB and OTX2, and 20 signals at FDR 25%, of which 14 replicated. Integrating genetic variation and T2D-discordant monozygotic twin analyses, we identify both genetic-based and genetic-independent T2D-DMPs. The signals annotate to genes with established GWAS and EWAS links to T2D and its complications, including blood pressure (RGL3) and eye disease (OTX2).

INTERPRETATION

The results help to improve our understanding of T2D disease pathogenesis and progression and may provide biomarkers for its complications.

FUNDING

Funding acknowledgements for each cohort can be found in the Supplementary Note.

Detecting genetic effects on phenotype variability to capture gene-by-environment interactions: a systematic method comparison

Genetically associated phenotypic variability has been widely observed across organisms and traits, including in humans. Both gene-gene and gene-environment interactions can lead to an increase in genetically associated phenotypic variability. Therefore, detecting the underlying genetic variants, or variance Quantitative Trait Loci (vQTLs), can provide novel insights into complex traits. Established approaches to detect vQTLs apply different methodologies from variance-only approaches to mean-variance joint tests, but a comprehensive comparison of these methods is lacking. Here, we review available methods to detect vQTLs in humans, carry out a simulation study to assess their performance under different biological scenarios of gene-environment interactions, and apply the optimal approaches for vQTL identification to gene expression data. Overall, with a minor allele frequency (MAF) of less than 0.2, the squared residual value linear model (SVLM) and the deviation regression model (DRM) are optimal when the data follow normal and non-normal distributions, respectively. In addition, the Brown-Forsythe (BF) test is one of the optimal methods when the MAF is 0.2 or larger, irrespective of phenotype distribution. Additionally, a larger sample size and more balanced sample distribution in different exposure categories increase the power of BF, SVLM, and DRM. Our results highlight vQTL detection methods that perform optimally under realistic simulation settings and show that their relative performance depends on the phenotype distribution, allele frequency, sample size, and the type of exposure in the interaction model underlying the vQTL.

Cross-sectional analyses of metabolites across biological samples mediating dietary acid load and chronic kidney disease

Chronic kidney disease (CKD) is a major public health burden, with dietary acid load (DAL) and gut microbiota playing crucial roles. As DAL can affect the host metabolome, potentially via the gut microbiota, we cross-sectionally investigated the interplay between DAL, host metabolome, gut microbiota, and early-stage CKD (TwinsUK, n = 1,453). DAL was positively associated with CKD stage G1-G2 (Beta (95% confidence interval) = 0.34 (0.007; 0.7), p = 0.046). After adjusting for covariates and multiple testing, we identified 15 serum, 14 urine, 8 stool, and 7 saliva metabolites, primarily lipids and amino acids, associated with both DAL and CKD progression. Of these, 8 serum, 2 urine, and one stool metabolites were found to mediate the DAL-CKD association. Furthermore, the stool metabolite 5-methylhexanoate (i7:0) correlated with 26 gut microbial species. Our findings emphasize the gut microbiota's therapeutic potential in countering DAL's impact on CKD through the host metabolome. Interventional and longitudinal studies are needed to establish causality.

Epigenome-wide Association Study Shows Differential DNA Methylation of MDC1, KLF9, and CUTA in Autoimmune Thyroid Disease

CONTEXT

Autoimmune thyroid disease (AITD) includes Graves disease (GD) and Hashimoto disease (HD), which often run in the same family. AITD etiology is incompletely understood: Genetic factors may account for up to 75% of phenotypic variance, whereas epigenetic effects (including DNA methylation [DNAm]) may contribute to the remaining variance (eg, why some individuals develop GD and others HD).

OBJECTIVE

This work aimed to identify differentially methylated positions (DMPs) and differentially methylated regions (DMRs) comparing GD to HD.

METHODS

Whole-blood DNAm was measured across the genome using the Infinium MethylationEPIC array in 32 Australian patients with GD and 30 with HD (discovery cohort) and 32 Danish patients with GD and 32 with HD (replication cohort). Linear mixed models were used to test for differences in quantile-normalized β values of DNAm between GD and HD and data were later meta-analyzed. Comb-p software was used to identify DMRs.

RESULTS

We identified epigenome-wide significant differences (P < 9E-8) and replicated (P < .05) 2 DMPs between GD and HD (cg06315208 within MDC1 and cg00049440 within KLF9). We identified and replicated a DMR within CUTA (5 CpGs at 6p21.32). We also identified 64 DMPs and 137 DMRs in the meta-analysis.

CONCLUSION

Our study reveals differences in DNAm between GD and HD, which may help explain why some people develop GD and others HD and provide a link to environmental risk factors. Additional research is needed to advance understanding of the role of DNAm in AITD and investigate its prognostic and therapeutic potential.

Usual dietary intake and change in DNA methylation over years: EWAS in KORA FF4 and KORA fit

Introduction

Changes in DNA methylation can increase or suppress the expression of health-relevant genes. We investigated for the first time the relationship between habitual food consumption and changes in DNA methylation.

Methods

The German KORA FF4 and KORA Fit studies were used to study the change in methylation over a median follow-up of 4 years. Only subjects participating in both surveys and with available dietary and methylation data were included in the analysis (n = 465). DNA methylation was measured using the Infinium MethylationEPIC BeadChip (Illumina), resulting in 735,527 shared CpGs across both studies. Generalized estimating equation models with an interaction term of exposure and time point were used to analyze the association of 34 food groups, folic acid, and two dietary patterns with changes in DNA methylation over time.

Results

The results were corrected for genomic inflation. Significant interaction terms indicate different effects between both time points. We observed only a few significant associations between food intake and change in DNA methylation, except for cream and spirit consumption. The annotated genes include CLN3, PROM1, DLEU7, TLL2, and UGT1A10.

Discussion

We identified weak associations between food consumption and DNA methylation change. The differential results for cream and spirits, both consumed in low quantities, require replication in independent studies.

Dietary and supplemental intake of vitamins C and E is associated with altered DNA methylation in an epigenome-wide association study meta-analysis

BACKGROUND

Dietary intake of antioxidants such as vitamins C and E protect against oxidative stress, and may also be associated with altered DNA methylation patterns.

METHODS

We meta-analysed epigenome-wide association study (EWAS) results from 11,866 participants across eight population-based cohorts to evaluate the association between self-reported dietary and supplemental intake of vitamins C and E with DNA methylation. EWAS were adjusted for age, sex, BMI, caloric intake, blood cell type proportion, smoking status, alcohol consumption, and technical covariates. Significant results of the meta-analysis were subsequently evaluated in gene set enrichment analysis (GSEA) and expression quantitative trait methylation (eQTM) analysis.

RESULTS

In meta-analysis, methylation at 4,656 CpG sites was significantly associated with vitamin C intake at FDR ≤ 0.05. The most significant CpG sites associated with vitamin C (at FDR ≤ 0.01) were enriched for pathways associated with systems development and cell signalling in GSEA, and were associated with downstream expression of genes enriched in the immune response in eQTM analysis. Furthermore, methylation at 160 CpG sites was significantly associated with vitamin E intake at FDR ≤ 0.05, but GSEA and eQTM analysis of the top most significant CpG sites associated with vitamin E did not identify significant enrichment of any biological pathways investigated.

CONCLUSIONS

We identified significant associations of many CpG sites with vitamin C and E intake, and our results suggest that vitamin C intake may be associated with systems development and the immune response.

A Fecal Metabolite Signature of Impaired Fasting Glucose: Results From Two Independent Population-Based Cohorts

Prediabetes is a metabolic condition associated with gut microbiome composition, although mechanisms remain elusive. We searched for fecal metabolites, a readout of gut microbiome function, associated with impaired fasting glucose (IFG) in 142 individuals with IFG and 1,105 healthy individuals from the UK Adult Twin Registry (TwinsUK). We used the Cooperative Health Research in the Region of Augsburg (KORA) cohort (318 IFG individuals, 689 healthy individuals) to replicate our findings. We linearly combined eight IFG-positively associated metabolites (1-methylxantine, nicotinate, glucuronate, uridine, cholesterol, serine, caffeine, and protoporphyrin IX) into an IFG-metabolite score, which was significantly associated with higher odds ratios (ORs) for IFG (TwinsUK: OR 3.9 [95% CI 3.02-5.02], P < 0.0001, KORA: OR 1.3 [95% CI 1.16-1.52], P < 0.0001) and incident type 2 diabetes (T2D; TwinsUK: hazard ratio 4 [95% CI 1.97-8], P = 0.0002). Although these are host-produced metabolites, we found that the gut microbiome is strongly associated with their fecal levels (area under the curve >70%). Abundances of Faecalibacillus intestinalis, Dorea formicigenerans, Ruminococcus torques, and Dorea sp. AF24-7LB were positively associated with IFG, and such associations were partially mediated by 1-methylxanthine and nicotinate (variance accounted for mean 14.4% [SD 5.1], P < 0.05). Our results suggest that the gut microbiome is linked to prediabetes not only via the production of microbial metabolites but also by affecting intestinal absorption/excretion of host-produced metabolites and xenobiotics, which are correlated with the risk of IFG. Fecal metabolites enable modeling of another mechanism of gut microbiome effect on prediabetes and T2D onset.

ARTICLE HIGHLIGHTS

Prediabetes is a metabolic condition associated with gut microbiome composition, although mechanisms remain elusive. We investigated whether there is a fecal metabolite signature of impaired fasting glucose (IFG) and the possible underlying mechanisms of action. We identified a fecal metabolite signature of IFG associated with prevalent IFG in two independent cohorts and incident type 2 diabetes in a subanalysis. Although the signature consists of metabolites of nonmicrobial origin, it is strongly correlated with gut microbiome composition. Fecal metabolites enable modeling of another mechanism of gut microbiome effect on prediabetes by affecting intestinal absorption or excretion of host compounds and xenobiotics.

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 [Formula: see text]: 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.

Metabolomic biomarkers of habitual B vitamin intakes unveil novel differentially methylated positions in the human epigenome

BACKGROUND

B vitamins such as folate (B9), B6, and B12 are key in one carbon metabolism, which generates methyl donors for DNA methylation. Several studies have linked differential methylation to self-reported intakes of folate and B12, but these estimates can be imprecise, while metabolomic biomarkers can offer an objective assessment of dietary intakes. We explored blood metabolomic biomarkers of folate and vitamins B6 and B12, to carry out epigenome-wide analyses across up to three European cohorts. Associations between self-reported habitual daily B vitamin intakes and 756 metabolites (Metabolon Inc.) were assessed in serum samples from 1064 UK participants from the TwinsUK cohort. The identified B vitamin metabolomic biomarkers were then used in epigenome-wide association tests with fasting blood DNA methylation levels at 430,768 sites from the Infinium HumanMethylation450 BeadChip in blood samples from 2182 European participants from the TwinsUK and KORA cohorts. Candidate signals were explored for metabolite associations with gene expression levels in a subset of the TwinsUK sample (n = 297). Metabolomic biomarker epigenetic associations were also compared with epigenetic associations of self-reported habitual B vitamin intakes in samples from 2294 European participants.

RESULTS

Eighteen metabolites were associated with B vitamin intakes after correction for multiple testing (Bonferroni-adj. p < 0.05), of which 7 metabolites were available in both cohorts and tested for epigenome-wide association. Three metabolites - pipecolate (metabolomic biomarker of B6 and folate intakes), pyridoxate (marker of B6 and folate) and docosahexaenoate (DHA, marker of B6) - were associated with 10, 3 and 1 differentially methylated positions (DMPs), respectively. The strongest association was observed between DHA and DMP cg03440556 in the SCD gene (effect = 0.093 ± 0.016, p = 4.07E-09). Pyridoxate, a catabolic product of vitamin B6, was inversely associated with CpG methylation near the SLC1A5 gene promoter region (cg02711608 and cg22304262) and with SLC7A11 (cg06690548), but not with corresponding changes in gene expression levels. The self-reported intake of folate and vitamin B6 had consistent but non-significant associations with the epigenetic signals.

CONCLUSION

Metabolomic biomarkers are a valuable approach to investigate the effects of dietary B vitamin intake on the human epigenome.

Leveraging osteoclast genetic regulatory data to identify genes with a role in osteoarthritis

There has been a growing interest in the role of the subchondral bone and its resident osteoclasts in the progression of osteoarthritis (OA). A recent genome-wide association study (GWAS) identified 100 independent association signals for OA traits. Most of these signals are led by noncoding variants, suggesting that genetic regulatory effects may drive many of the associations. We have generated a unique human osteoclast-like cell-specific expression quantitative trait locus (eQTL) resource for studying the genetics of bone disease. Considering the potential role of osteoclasts in the pathogenesis of OA, we performed an integrative analysis of this dataset with the recently published OA GWAS results. Summary data-based Mendelian randomization (SMR) and colocalization analyses identified 38 genes with a potential role in OA, including some that have been implicated in Mendelian diseases with joint/skeletal abnormalities, such as BICRA, EIF6, CHST3, and FBN2. Several OA GWAS signals demonstrated colocalization with more than one eQTL peak, including at 19q13.32 (hip OA with BCAM, PRKD2, and BICRA eQTL). We also identified a number of eQTL signals colocalizing with more than one OA trait, including FAM53A, GCAT, HMGN1, MGAT4A, RRP7BP, and TRIOBP. An SMR analysis identified 3 loci with evidence of pleiotropic effects on OA-risk and gene expression: LINC01481, CPNE1, and EIF6. Both CPNE1 and EIF6 are located at 20q11.22, a locus harboring 2 other strong OA candidate genes, GDF5 and UQCC1, suggesting the presence of an OA-risk gene cluster. In summary, we have used our osteoclast-specific eQTL dataset to identify genes potentially involved with the pathogenesis of OA.

Hemin availability induces coordinated DNA methylation and gene expression changes in Porphyromonas gingivalis

Periodontal disease is a chronic inflammatory disease in which the oral pathogen Porphyromonas gingivalis plays an important role. Porphyromonas gingivalis expresses virulence determinants in response to higher hemin concentrations, but the underlying regulatory processes remain unclear. Bacterial DNA methylation has the potential to fulfil this mechanistic role. We characterized the methylome of P. gingivalis, and compared its variation to transcriptome changes in response to hemin availability. Porphyromonas gingivalis W50 was grown in chemostat continuous culture with excess or limited hemin, prior to whole-methylome and transcriptome profiling using Nanopore and Illumina RNA-Seq. DNA methylation was quantified for Dam/Dcm motifs and all-context N6-methyladenine (6mA) and 5-methylcytosine (5mC). Of all 1,992 genes analyzed, 161 and 268 were respectively over- and under-expressed with excess hemin. Notably, we detected differential DNA methylation signatures for the Dam "GATC" motif and both all-context 6mA and 5mC in response to hemin availability. Joint analyses identified a subset of coordinated changes in gene expression, 6mA, and 5mC methylation that target genes involved in lactate utilization and ABC transporters. The results identify altered methylation and expression responses to hemin availability in P. gingivalis, with insights into mechanisms regulating its virulence in periodontal disease. IMPORTANCE DNA methylation has important roles in bacteria, including in the regulation of transcription. Porphyromonas gingivalis, an oral pathogen in periodontitis, exhibits well-established gene expression changes in response to hemin availability. However, the regulatory processes underlying these effects remain unknown. We profiled the novel P. gingivalis epigenome, and assessed epigenetic and transcriptome variation under limited and excess hemin conditions. As expected, multiple gene expression changes were detected in response to limited and excess hemin that reflect health and disease, respectively. Notably, we also detected differential DNA methylation signatures for the Dam "GATC" motif and both all-context 6mA and 5mC in response to hemin. Joint analyses identified coordinated changes in gene expression, 6mA, and 5mC methylation that target genes involved in lactate utilization and ABC transporters. The results identify novel regulatory processes underlying the mechanism of hemin regulated gene expression in P. gingivalis, with phenotypic impacts on its virulence in periodontal disease.

Stratified genome-wide association analysis of type 2 diabetes reveals subgroups with genetic and environmental heterogeneity

Type 2 diabetes (T2D) is a heterogeneous illness caused by genetic and environmental factors. Previous genome-wide association studies (GWAS) have identified many genetic variants associated with T2D and found evidence of differing genetic profiles by age-at-onset. This study seeks to explore further the genetic and environmental drivers of T2D by analyzing subgroups on the basis of age-at-onset of diabetes and body mass index (BMI). In the UK Biobank, 36 494 T2D cases were stratified into three subgroups, and GWAS was performed for all T2D cases and for each subgroup relative to 421 021 controls. Altogether, 18 single nucleotide polymorphisms were significantly associated with T2D genome-wide in one or more subgroups and also showed evidence of heterogeneity between the subgroups (Cochrane's Q P < 0.01), with two SNPs remaining significant after multiple testing (in CDKN2B and CYTIP). Combined risk scores, on the basis of genetic profile, BMI and age, resulted in excellent diabetes prediction [area under the ROC curve (AUC) = 0.92]. A modest improvement in prediction (AUC = 0.93) was seen when the contribution of genetic and environmental factors was evaluated separately for each subgroup. Increasing sample sizes of genetic studies enables us to stratify disease cases into subgroups, which have sufficient power to highlight areas of genetic heterogeneity. Despite some evidence that optimizing combined risk scores by subgroup improves prediction, larger sample sizes are likely needed for prediction when using a stratification approach.

Genetic impacts on DNA methylation help elucidate regulatory genomic processes

BACKGROUND

Pinpointing genetic impacts on DNA methylation can improve our understanding of pathways that underlie gene regulation and disease risk.

RESULTS

We report heritability and methylation quantitative trait locus (meQTL) analysis at 724,499 CpGs profiled with the Illumina Infinium MethylationEPIC array in 2358 blood samples from three UK cohorts. Methylation levels at 34.2% of CpGs are affected by SNPs, and 98% of effects are cis-acting or within 1 Mbp of the tested CpG. Our results are consistent with meQTL analyses based on the former Illumina Infinium HumanMethylation450 array. Both SNPs and CpGs with meQTLs are overrepresented in enhancers, which have improved coverage on this platform compared to previous approaches. Co-localisation analyses across genetic effects on DNA methylation and 56 human traits identify 1520 co-localisations across 1325 unique CpGs and 34 phenotypes, including in disease-relevant genes, such as USP1 and DOCK7 (total cholesterol levels), and ICOSLG (inflammatory bowel disease). Enrichment analysis of meQTLs and integration with expression QTLs give insights into mechanisms underlying cis-meQTLs (e.g. through disruption of transcription factor binding sites for CTCF and SMC3) and trans-meQTLs (e.g. through regulating the expression of ACD and SENP7 which can modulate DNA methylation at distal sites).

CONCLUSIONS

Our findings improve the characterisation of the mechanisms underlying DNA methylation variability and are informative for prioritisation of GWAS variants for functional follow-ups. The MeQTL EPIC Database and viewer are available online at https://epicmeqtl.kcl.ac.uk .

Integrative genomic analyses in adipocytes implicate DNA methylation in human obesity and diabetes

DNA methylation variations are prevalent in human obesity but evidence of a causative role in disease pathogenesis is limited. Here, we combine epigenome-wide association and integrative genomics to investigate the impact of adipocyte DNA methylation variations in human obesity. We discover extensive DNA methylation changes that are robustly associated with obesity (N = 190 samples, 691 loci in subcutaneous and 173 loci in visceral adipocytes, P < 1 × 10-7). We connect obesity-associated methylation variations to transcriptomic changes at >500 target genes, and identify putative methylation-transcription factor interactions. Through Mendelian Randomisation, we infer causal effects of methylation on obesity and obesity-induced metabolic disturbances at 59 independent loci. Targeted methylation sequencing, CRISPR-activation and gene silencing in adipocytes, further identifies regional methylation variations, underlying regulatory elements and novel cellular metabolic effects. Our results indicate DNA methylation is an important determinant of human obesity and its metabolic complications, and reveal mechanisms through which altered methylation may impact adipocyte functions.

DNA methylation analysis is used to identify novel genetic loci associated with circulating fibrinogen levels in blood

BACKGROUND

Fibrinogen plays an essential role in blood coagulation and inflammation. Circulating fibrinogen levels may be determined based on interindividual differences in DNA methylation at cytosine-phosphate-guanine (CpG) sites and vice versa.

OBJECTIVES

To perform an EWAS to examine an association between blood DNA methylation levels and circulating fibrinogen levels to better understand its biological and pathophysiological actions.

METHODS

We performed an epigenome-wide association study of circulating fibrinogen levels in 18 037 White, Black, American Indian, and Hispanic participants, representing 14 studies from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium. Circulating leukocyte DNA methylation was measured using the Illumina 450K array in 12 904 participants and using the EPIC array in 5133 participants. In each study, an epigenome-wide association study of fibrinogen was performed using linear mixed models adjusted for potential confounders. Study-specific results were combined using array-specific meta-analysis, followed by cross-replication of epigenome-wide significant associations. We compared models with and without CRP adjustment to examine the role of inflammation.

RESULTS

We identified 208 and 87 significant CpG sites associated with fibrinogen levels from the 450K (p < 1.03 × 10-7) and EPIC arrays (p < 5.78 × 10-8), respectively. There were 78 associations from the 450K array that replicated in the EPIC array and 26 vice versa. After accounting for overlapping sites, there were 83 replicated CpG sites located in 61 loci, of which only 4 have been previously reported for fibrinogen. The examples of genes located near these CpG sites were SOCS3 and AIM2, which are involved in inflammatory pathways. The associations of all 83 replicated CpG sites were attenuated after CRP adjustment, although many remained significant.

CONCLUSION

We identified 83 CpG sites associated with circulating fibrinogen levels. These associations are partially driven by inflammatory pathways shared by both fibrinogen and CRP.

Interplay between the (Poly)phenol Metabolome, Gut Microbiome, and Cardiovascular Health in Women: A Cross-Sectional Study from the TwinsUK Cohort

BACKGROUND

Dietary (poly)phenol consumption is inversely associated with cardiovascular disease (CVD) risk in epidemiological studies, but little is known about the role of the gut microbiome in this relationship.

METHODS

In 200 healthy females, aged 62.0 ± 10.0 years, from the TwinsUK cohort, 114 individual (poly)phenol metabolites were measured from spot urine using ultra-high-performance liquid chromatography-mass spectrometry. The associations between metabolites, the gut microbiome (alpha diversity and genera), and cardiovascular scores were investigated using linear mixed models adjusting age, BMI, fibre, energy intake, family relatedness, and multiple testing (FDR < 0.1).

RESULTS

Significant associations were found between phenolic acid metabolites, CVD risk, and the gut microbiome. A total of 35 phenolic acid metabolites were associated with the Firmicutes phylum, while 5 metabolites were associated with alpha diversity (FDR-adjusted p < 0.05). Negative associations were observed between the atherosclerotic CVD (ASCVD) risk score and five phenolic acid metabolites, two tyrosol metabolites, and daidzein with stdBeta (95% (CI)) ranging from -0.05 (-0.09, -0.01) for 3-(2,4-dihydroxyphenyl)propanoic acid to -0.04 (-0.08, -0.003) for 2-hydroxycinnamic acid (FDR-adjusted p < 0.1). The genus 5-7N15 in the Bacteroidetes phylum was positively associated with the same metabolites, including 3-(3,5-dihydroxyphenyl)propanoic acid, 3-(2,4-dihydroxyphenyl)propanoic acid, 3-(3,4-dihydroxyphenyl)propanoic acid), 3-hydroxyphenylethanol-4-sulfate, and 4-hydroxyphenylethanol-3-sulfate)(stdBeta (95% CI): 0.23 (0.09, 0.36) to 0.28 (0.15, 0.42), FDR-adjusted p < 0.05), and negatively associated with the ASCVD score (stdBeta (95% CI): -0.05 (-0.09, -0.01), FDR-adjusted p = 0.02). Mediation analysis showed that genus 5-7N15 mediated 23.8% of the total effect of 3-(3,4-dihydroxyphenyl)propanoic acid on the ASCVD score.

CONCLUSIONS

Coffee, tea, red wine, and several vegetables and fruits, especially berries, are the most abundant food sources of phenolic acids that have the strongest associations with CVD risk. We found that the gut microbiome, particularly the genus 5-7N15, partially mediates the negative association between urinary (poly)phenols and cardiovascular risk, supporting a key role of the gut microbiome in the health benefits of dietary (poly)phenols.

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.

Epigenome-Wide Association Study Reveals CpG Sites Associated with Thyroid Function and Regulatory Effects on KLF9

Background: Thyroid hormones play a key role in differentiation and metabolism and are known regulators of gene expression through both genomic and epigenetic processes including DNA methylation. The aim of this study was to examine associations between thyroid hormones and DNA methylation. Methods: We carried out a fixed-effect meta-analysis of epigenome-wide association study (EWAS) of blood DNA methylation sites from 8 cohorts from the ThyroidOmics Consortium, incorporating up to 7073 participants of both European and African ancestry, implementing a discovery and replication stage. Statistical analyses were conducted using normalized beta CpG values as dependent and log-transformed thyrotropin (TSH), free thyroxine, and free triiodothyronine levels, respectively, as independent variable in a linear model. The replicated findings were correlated with gene expression levels in whole blood and tested for causal influence of TSH and free thyroxine by two-sample Mendelian randomization (MR). Results: Epigenome-wide significant associations (p-value <1.1E-7) of three CpGs for free thyroxine, five for free triiodothyronine, and two for TSH concentrations were discovered and replicated (combined p-values = 1.5E-9 to 4.3E-28). The associations included CpG sites annotated to KLF9 (cg00049440) and DOT1L (cg04173586) that overlap with all three traits, consistent with hypothalamic-pituitary-thyroid axis physiology. Significant associations were also found for CpGs in FKBP5 for free thyroxine, and at CSNK1D/LINCO1970 and LRRC8D for free triiodothyronine. MR analyses supported a causal effect of thyroid status on DNA methylation of KLF9. DNA methylation of cg00049440 in KLF9 was inversely correlated with KLF9 gene expression in blood. The CpG at CSNK1D/LINC01970 overlapped with thyroid hormone receptor alpha binding peaks in liver cells. The total additive heritability of the methylation levels of the six significant CpG sites was between 25% and 57%. Significant methylation QTLs were identified for CpGs at KLF9, FKBP5, LRRC8D, and CSNK1D/LINC01970. Conclusions: We report novel associations between TSH, thyroid hormones, and blood-based DNA methylation. This study advances our understanding of thyroid hormone action particularly related to KLF9 and serves as a proof-of-concept that integrations of EWAS with other -omics data can provide a valuable tool for unraveling thyroid hormone signaling in humans by complementing and feeding classical in vitro and animal studies.

Accelerated epigenetic aging and DNA methylation alterations in Berardinelli-Seip congenital lipodystrophy

Berardinelli-Seip Congenital Lipodystrophy type 2 (CGL2) is a very rare human genetic disorder with potential significance to the understanding of the pathobiology of aging. CGL2 patients display characteristic progeroid features and suffer from type 2 diabetes, insulin resistance, and fatty liver. In this study, we profiled genome-wide DNA methylation levels in CGL2 patients with BSCL2 mutations to study epigenetic age acceleration and DNA methylation alterations. This analysis revealed significant age acceleration in blood DNA of CGL2 patients using both first and second generation epigenetic clocks. We also observed a shortened lifespan of C. elegans following knockdown of the BSCL2 homolog seip-1 on a daf-16/Foxo mutant background. DNA methylation analysis revealed significant differentially methylated sites enriched for lyase activity, kinase regulator activity, protein kinase regulator activity, and kinase activator activity. We could also observe significant hypomethylation in the promoter of the Dual Specificity Phosphatase 22 (DUSP22) gene when comparing CGL2 patients vs controls. We conclude that in line with the observed progeroid features, CGL2 patients exhibit significant epigenetic age acceleration and DNA methylation alterations that might affect pathways/genes of potential relevance to the disease.

Age acquired skewed X chromosome inactivation is associated with adverse health outcomes in humans

BACKGROUND

Ageing is a heterogenous process characterised by cellular and molecular hallmarks, including changes to haematopoietic stem cells and is a primary risk factor for chronic diseases. X chromosome inactivation (XCI) randomly transcriptionally silences either the maternal or paternal X in each cell of 46, XX females to balance the gene expression with 46, XY males. Age acquired XCI-skew describes the preferential selection of cells across a tissue resulting in an imbalance of XCI, which is particularly prevalent in blood tissues of ageing females, and yet its clinical consequences are unknown.

METHODS

We assayed XCI in 1575 females from the TwinsUK population cohort using DNA extracted from whole blood. We employed prospective, cross-sectional, and intra-twin study designs to characterise the relationship of XCI-skew with molecular and cellular measures of ageing, cardiovascular disease risk, and cancer diagnosis.

RESULTS

We demonstrate that XCI-skew is independent of traditional markers of biological ageing and is associated with a haematopoietic bias towards the myeloid lineage. Using an atherosclerotic cardiovascular disease risk score, which captures traditional risk factors, XCI-skew is associated with an increased cardiovascular disease risk both cross-sectionally and within XCI-skew discordant twin pairs. In a prospective 10 year follow-up study, XCI-skew is predictive of future cancer incidence.

CONCLUSIONS

Our study demonstrates that age acquired XCI-skew captures changes to the haematopoietic stem cell population and has clinical potential as a unique biomarker of chronic disease risk.

FUNDING

KSS acknowledges funding from the Medical Research Council [MR/M004422/1 and MR/R023131/1]. JTB acknowledges funding from the ESRC [ES/N000404/1]. MM acknowledges funding from the National Institute for Health Research (NIHR)-funded BioResource, Clinical Research Facility and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. TwinsUK is funded by the Wellcome Trust, Medical Research Council, European Union, Chronic Disease Research Foundation (CDRF), Zoe Global Ltd and the National Institute for Health Research (NIHR)-funded BioResource, Clinical Research Facility and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London.

Pulmonary Function and Blood DNA Methylation: A Multiancestry Epigenome-Wide Association Meta-analysis

Rationale: Methylation integrates factors present at birth and modifiable across the lifespan that can influence pulmonary function. Studies are limited in scope and replication. Objectives: To conduct large-scale epigenome-wide meta-analyses of blood DNA methylation and pulmonary function. Methods: Twelve cohorts analyzed associations of methylation at cytosine-phosphate-guanine probes (CpGs), using Illumina 450K or EPIC/850K arrays, with FEV₁, FVC, and FEV₁/FVC. We performed multiancestry epigenome-wide meta-analyses (total of 17,503 individuals; 14,761 European, 2,549 African, and 193 Hispanic/Latino ancestries) and interpreted results using integrative epigenomics. Measurements and Main Results: We identified 1,267 CpGs (1,042 genes) differentially methylated (false discovery rate, <0.025) in relation to FEV₁, FVC, or FEV₁/FVC, including 1,240 novel and 73 also related to chronic obstructive pulmonary disease (1,787 cases). We found 294 CpGs unique to European or African ancestry and 395 CpGs unique to never or ever smokers. The majority of significant CpGs correlated with nearby gene expression in blood. Findings were enriched in key regulatory elements for gene function, including accessible chromatin elements, in both blood and lung. Sixty-nine implicated genes are targets of investigational or approved drugs. One example novel gene highlighted by integrative epigenomic and druggable target analysis is TNFRSF4. Mendelian randomization and colocalization analyses suggest that epigenome-wide association study signals capture causal regulatory genomic loci. Conclusions: We identified numerous novel loci differentially methylated in relation to pulmonary function; few were detected in large genome-wide association studies. Integrative analyses highlight functional relevance and potential therapeutic targets. This comprehensive discovery of potentially modifiable, novel lung function loci expands knowledge gained from genetic studies, providing insights into lung pathogenesis.

ACE2 expression in adipose tissue is associated with cardio-metabolic risk factors and cell type composition-implications for COVID-19

BACKGROUND

COVID-19 severity varies widely. Although some demographic and cardio-metabolic factors, including age and obesity, are associated with increasing risk of severe illness, the underlying mechanism(s) are uncertain.

SUBJECTS/METHODS

In a meta-analysis of three independent studies of 1471 participants in total, we investigated phenotypic and genetic factors associated with subcutaneous adipose tissue expression of Angiotensin I Converting Enzyme 2 (ACE2), measured by RNA-Seq, which acts as a receptor for SARS-CoV-2 cellular entry.

RESULTS

Lower adipose tissue ACE2 expression was associated with multiple adverse cardio-metabolic health indices, including type 2 diabetes (T2D) (P = 9.14 × 10-6), obesity status (P = 4.81 × 10-5), higher serum fasting insulin (P = 5.32 × 10-4), BMI (P = 3.94 × 10-4), and lower serum HDL levels (P = 1.92 × 10-7). ACE2 expression was also associated with estimated proportions of cell types in adipose tissue: lower expression was associated with a lower proportion of microvascular endothelial cells (P = 4.25 × 10-4) and higher proportion of macrophages (P = 2.74 × 10-5). Despite an estimated heritability of 32%, we did not identify any proximal or distal expression quantitative trait loci (eQTLs) associated with adipose tissue ACE2 expression.

CONCLUSIONS

Our results demonstrate that individuals with cardio-metabolic features known to increase risk of severe COVID-19 have lower background ACE2 levels in this highly relevant tissue. Reduced adipose tissue ACE2 expression may contribute to the pathophysiology of cardio-metabolic diseases, as well as the associated increased risk of severe COVID-19.

Adipose methylome integrative-omic analyses reveal genetic and dietary metabolic health drivers and insulin resistance classifiers

BACKGROUND

There is considerable evidence for the importance of the DNA methylome in metabolic health, for example, a robust methylation signature has been associated with body mass index (BMI). However, visceral fat (VF) mass accumulation is a greater risk factor for metabolic disease than BMI alone. In this study, we dissect the subcutaneous adipose tissue (SAT) methylome signature relevant to metabolic health by focusing on VF as the major risk factor of metabolic disease. We integrate results with genetic, blood methylation, SAT gene expression, blood metabolomic, dietary intake and metabolic phenotype data to assess and quantify genetic and environmental drivers of the identified signals, as well as their potential functional roles.

METHODS

Epigenome-wide association analyses were carried out to determine visceral fat mass-associated differentially methylated positions (VF-DMPs) in SAT samples from 538 TwinsUK participants. Validation and replication were performed in 333 individuals from 3 independent cohorts. To assess functional impacts of the VF-DMPs, the association between VF and gene expression was determined at the genes annotated to the VF-DMPs and an association analysis was carried out to determine whether methylation at the VF-DMPs is associated with gene expression. Further epigenetic analyses were carried out to compare methylation levels at the VF-DMPs as the response variables and a range of different metabolic health phenotypes including android:gynoid fat ratio (AGR), lipids, blood metabolomic profiles, insulin resistance, T2D and dietary intake variables. The results from all analyses were integrated to identify signals that exhibit altered SAT function and have strong relevance to metabolic health.

RESULTS

We identified 1181 CpG positions in 788 genes to be differentially methylated with VF (VF-DMPs) with significant enrichment in the insulin signalling pathway. Follow-up cross-omic analysis of VF-DMPs integrating genetics, gene expression, metabolomics, diet, and metabolic traits highlighted VF-DMPs located in 9 genes with strong relevance to metabolic disease mechanisms, with replication of signals in FASN, SREBF1, TAGLN2, PC and CFAP410. PC methylation showed evidence for mediating effects of diet on VF. FASN DNA methylation exhibited putative causal effects on VF that were also strongly associated with insulin resistance and methylation levels in FASN better classified insulin resistance (AUC=0.91) than BMI or VF alone.

CONCLUSIONS

Our findings help characterise the adiposity-associated methylation signature of SAT, with insights for metabolic disease risk.

DNA methylation signature of chronic low-grade inflammation and its role in cardio-respiratory diseases

We performed a multi-ethnic Epigenome Wide Association study on 22,774 individuals to describe the DNA methylation signature of chronic low-grade inflammation as measured by C-Reactive protein (CRP). We find 1,511 independent differentially methylated loci associated with CRP. These CpG sites show correlation structures across chromosomes, and are primarily situated in euchromatin, depleted in CpG islands. These genomic loci are predominantly situated in transcription factor binding sites and genomic enhancer regions. Mendelian randomization analysis suggests altered CpG methylation is a consequence of increased blood CRP levels. Mediation analysis reveals obesity and smoking as important underlying driving factors for changed CpG methylation. Finally, we find that an activated CpG signature significantly increases the risk for cardiometabolic diseases and COPD.

Meta-analysis of epigenome-wide association studies of carotid intima-media thickness

Common carotid intima-media thickness (cIMT) is an index of subclinical atherosclerosis that is associated with ischemic stroke and coronary artery disease (CAD). We undertook a cross-sectional epigenome-wide association study (EWAS) of measures of cIMT in 6400 individuals. Mendelian randomization analysis was applied to investigate the potential causal role of DNA methylation in the link between atherosclerotic cardiovascular risk factors and cIMT or clinical cardiovascular disease. The CpG site cg05575921 was associated with cIMT (beta = -0.0264, p value = 3.5 × 10^-8) in the discovery panel and was replicated in replication panel (beta = -0.07, p value = 0.005). This CpG is located at chr5:81649347 in the intron 3 of the aryl hydrocarbon receptor repressor gene (AHRR). Our results indicate that DNA methylation at cg05575921 might be in the pathway between smoking, cIMT and stroke. Moreover, in a region-based analysis, 34 differentially methylated regions (DMRs) were identified of which a DMR upstream of ALOX12 showed the strongest association with cIMT (p value = 1.4 × 10^-13). In conclusion, our study suggests that DNA methylation may play a role in the link between cardiovascular risk factors, cIMT and clinical cardiovascular disease.

Validating biomarkers and models for epigenetic inference of alcohol consumption from blood

Background

Information on long-term alcohol consumption is relevant for medical and public health research, disease therapy, and other areas. Recently, DNA methylation-based inference of alcohol consumption from blood was reported with high accuracy, but these results were based on employing the same dataset for model training and testing, which can lead to accuracy overestimation. Moreover, only subsets of alcohol consumption categories were used, which makes it impossible to extrapolate such models to the general population. By using data from eight population-based European cohorts (N = 4677), we internally and externally validated the previously reported biomarkers and models for epigenetic inference of alcohol consumption from blood and developed new models comprising all data from all categories.

Results

By employing data from six European cohorts (N = 2883), we empirically tested the reproducibility of the previously suggested biomarkers and prediction models via ten-fold internal cross-validation. In contrast to previous findings, all seven models based on 144-CpGs yielded lower mean AUCs compared to the models with less CpGs. For instance, the 144-CpG heavy versus non-drinkers model gave an AUC of 0.78 ± 0.06, while the 5 and 23 CpG models achieved 0.83 ± 0.05, respectively. The transportability of the models was empirically tested via external validation in three independent European cohorts (N = 1794), revealing high AUC variance between datasets within models. For instance, the 144-CpG heavy versus non-drinkers model yielded AUCs ranging from 0.60 to 0.84 between datasets. The newly developed models that considered data from all categories showed low AUCs but gave low AUC variation in the external validation. For instance, the 144-CpG heavy and at-risk versus light and non-drinkers model achieved AUCs of 0.67 ± 0.02 in the internal cross-validation and 0.61–0.66 in the external validation datasets.

Conclusions

The outcomes of our internal and external validation demonstrate that the previously reported prediction models suffer from both overfitting and accuracy overestimation. Our results show that the previously proposed biomarkers are not yet sufficient for accurate and robust inference of alcohol consumption from blood. Overall, our findings imply that DNA methylation prediction biomarkers and models need to be improved considerably before epigenetic inference of alcohol consumption from blood can be considered for practical applications.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-021-01186-3.

DNA methylation signatures of incident coronary heart disease: findings from epigenome-wide association studies

Coronary heart disease (CHD) is a type of cardiovascular disease (CVD) that affects the coronary arteries, which provide oxygenated blood to the heart. It is a major cause of mortality worldwide. Various prediction methods have been developed to assess the likelihood of developing CHD, including those based on clinical features and genetic variation. Recent epigenome-wide studies have identified DNA methylation signatures associated with the development of CHD, indicating that DNA methylation may play a role in predicting future CHD. This narrative review summarises recent findings from DNA methylation studies of incident CHD (iCHD) events from epigenome-wide association studies (EWASs). The results suggest that DNA methylation signatures may identify new mechanisms involved in CHD progression and could prove a useful adjunct for the prediction of future CHD.

Identical twins carry a persistent epigenetic signature of early genome programming

Monozygotic (MZ) twins and higher-order multiples arise when a zygote splits during pre-implantation stages of development. The mechanisms underpinning this event have remained a mystery. Because MZ twinning rarely runs in families, the leading hypothesis is that it occurs at random. Here, we show that MZ twinning is strongly associated with a stable DNA methylation signature in adult somatic tissues. This signature spans regions near telomeres and centromeres, Polycomb-repressed regions and heterochromatin, genes involved in cell-adhesion, WNT signaling, cell fate, and putative human metastable epialleles. Our study also demonstrates a never-anticipated corollary: because identical twins keep a lifelong molecular signature, we can retrospectively diagnose if a person was conceived as monozygotic twin.

Estrogen and COVID-19 symptoms: Associations in women from the COVID Symptom Study

It has been widely observed that adult men of all ages are at higher risk of developing serious complications from COVID-19 when compared with women. This study aimed to investigate the association of COVID-19 positivity and severity with estrogen exposure in women, in a population based matched cohort study of female users of the COVID Symptom Study application in the UK. Analyses included 152,637 women for menopausal status, 295,689 women for exogenous estrogen intake in the form of the combined oral contraceptive pill (COCP), and 151,193 menopausal women for hormone replacement therapy (HRT). Data were collected using the COVID Symptom Study in May-June 2020. Analyses investigated associations between predicted or tested COVID-19 status and menopausal status, COCP use, and HRT use, adjusting for age, smoking and BMI, with follow-up age sensitivity analysis, and validation in a subset of participants from the TwinsUK cohort. Menopausal women had higher rates of predicted COVID-19 (P = 0.003). COCP-users had lower rates of predicted COVID-19 (P = 8.03E-05), with reduction in hospital attendance (P = 0.023). Menopausal women using HRT or hormonal therapies did not exhibit consistent associations, including increased rates of predicted COVID-19 (P = 2.22E-05) for HRT users alone. The findings support a protective effect of estrogen exposure on COVID-19, based on positive association between predicted COVID-19 with menopausal status, and negative association with COCP use. HRT use was positively associated with COVID-19, but the results should be considered with caution due to lack of data on HRT type, route of administration, duration of treatment, and potential unaccounted for confounders and comorbidities.

Genomic and phenotypic insights from an atlas of genetic effects on DNA methylation

Characterizing genetic influences on DNA methylation (DNAm) provides an opportunity to understand mechanisms underpinning gene regulation and disease. In the present study, we describe results of DNAm quantitative trait locus (mQTL) analyses on 32,851 participants, identifying genetic variants associated with DNAm at 420,509 DNAm sites in blood. We present a database of >270,000 independent mQTLs, of which 8.5% comprise long-range (trans) associations. Identified mQTL associations explain 15-17% of the additive genetic variance of DNAm. We show that the genetic architecture of DNAm levels is highly polygenic. Using shared genetic control between distal DNAm sites, we constructed networks, identifying 405 discrete genomic communities enriched for genomic annotations and complex traits. Shared genetic variants are associated with both DNAm levels and complex diseases, but only in a minority of cases do these associations reflect causal relationships from DNAm to trait or vice versa, indicating a more complex genotype-phenotype map than previously anticipated.

Childhood growth and development and DNA methylation age in mid-life

BACKGROUND

In the first study of its kind, we examine the association between growth and development in early life and DNAm age biomarkers in mid-life.

METHODS

Participants were from the Medical Research Council National Survey of Health and Development (n = 1376). Four DNAm age acceleration (AgeAccel) biomarkers were measured when participants were aged 53 years: AgeAccelHannum; AgeAccelHorvath; AgeAccelLevine; and AgeAccelGrim. Exposure variables included: relative weight gain (standardised residuals from models of current weight z-score on current height, and previous weight and height z-scores); and linear growth (standardised residuals from models of current height z-score on previous height and weight z-scores) during infancy (0-2 years, weight gain only), early childhood (2-4 years), middle childhood (4-7 years) and late childhood to adolescence (7-15 years); age at menarche; and pubertal stage for men at 14-15 years. The relationship between relative weight gain and linear growth and AgeAccel was investigated using conditional growth models. We replicated analyses from the late childhood to adolescence period and pubertal timing among 240 participants from The National Child and Development Study (NCDS).

RESULTS

A 1SD increase in relative weight gain in late childhood to adolescence was associated with 0.50 years (95% CI 0.20, 0.79) higher AgeAccelGrim. Although the CI includes the null, the estimate was similar in NCDS [0.57 years (95% CI - 0.01, 1.16)] There was no strong evidence that relative weight gain and linear growth in childhood was associated with any other AgeAccel biomarker. There was no relationship between pubertal timing in men and AgeAccel biomarkers. Women who reached menarche ≥ 12 years had 1.20 years (95% CI 0.15, 2.24) higher AgeAccelGrim on average than women who reached menarche < 12 years; however, this was not replicated in NCDS and was not statistically significant after Bonferroni correction.

CONCLUSIONS

Our findings generally do not support an association between growth and AgeAccel biomarkers in mid-life. However, we found rapid weight gain during pubertal development, previously related to higher cardiovascular disease risk, to be associated with older AgeAccelGrim. Given this is an exploratory study, this finding requires replication.

Epigenome-wide association study of diet quality in the Women's Health Initiative and TwinsUK cohort

BACKGROUND

Diet quality is a risk factor for chronic disease and mortality. Differential DNA methylation across the epigenome has been associated with chronic disease risk. Whether diet quality is associated with differential methylation is unknown. This study assessed whether diet quality was associated with differential DNA methylation measured across 445 548 loci in the Women's Health Initiative (WHI) and the TwinsUK cohort.

DESIGN

The discovery cohort consisted of 4355 women from the WHI. The replication cohort consisted of 571 mono- and dizygotic twins from the TwinsUK cohort. DNA methylation was measured in whole blood using the Illumina Infinium HumanMethylation450 Beadchip. Diet quality was assessed using the Alternative Healthy Eating Index 2010 (AHEI-2010). A meta-analysis, stratified by study cohort, was performed using generalized linear models that regressed methylation on AHEI-2010, adjusting for cell composition, chip number and location, study characteristics, principal components of genetic relatedness, age, smoking status, race/ethnicity and body mass index (BMI). Statistical significance was defined as a false discovery rate < 0.05. Significant sites were tested for replication in the TwinsUK cohort, with significant replication defined by P < 0.05 and a consistent direction.

RESULTS

Diet quality was significantly associated with differential DNA methylation at 428 cytosine-phosphate-guanine (CpG) sites in the discovery cohort. A total of 24 CpG sites were consistent with replication in the TwinsUK cohort, more than would be expected by chance (P = 2.7x10-4), with one site replicated in both the blood and adipose tissue (cg16379999 located in the body of SEL1L).

CONCLUSIONS

Diet quality was associated with methylation at 24 CpG sites, several of which have been associated with adiposity, inflammation and dysglycaemia. These findings may provide insight into pathways through which diet influences chronic disease.

Genome-wide association studies identify 137 genetic loci for DNA methylation biomarkers of aging

BACKGROUND

Biological aging estimators derived from DNA methylation data are heritable and correlate with morbidity and mortality. Consequently, identification of genetic and environmental contributors to the variation in these measures in populations has become a major goal in the field.

RESULTS

Leveraging DNA methylation and SNP data from more than 40,000 individuals, we identify 137 genome-wide significant loci, of which 113 are novel, from genome-wide association study (GWAS) meta-analyses of four epigenetic clocks and epigenetic surrogate markers for granulocyte proportions and plasminogen activator inhibitor 1 levels, respectively. We find evidence for shared genetic loci associated with the Horvath clock and expression of transcripts encoding genes linked to lipid metabolism and immune function. Notably, these loci are independent of those reported to regulate DNA methylation levels at constituent clock CpGs. A polygenic score for GrimAge acceleration showed strong associations with adiposity-related traits, educational attainment, parental longevity, and C-reactive protein levels.

CONCLUSION

This study illuminates the genetic architecture underlying epigenetic aging and its shared genetic contributions with lifestyle factors and longevity.

A multi-ethnic epigenome-wide association study of leukocyte DNA methylation and blood lipids

Here we examine the association between DNA methylation in circulating leukocytes and blood lipids in a multi-ethnic sample of 16,265 subjects. We identify 148, 35, and 4 novel associations among Europeans, African Americans, and Hispanics, respectively, and an additional 186 novel associations through a trans-ethnic meta-analysis. We observe a high concordance in the direction of effects across racial/ethnic groups, a high correlation of effect sizes between high-density lipoprotein and triglycerides, a modest overlap of associations with epigenome-wide association studies of other cardio-metabolic traits, and a largely non-overlap with lipid loci identified to date through genome-wide association studies. Thirty CpGs reached significance in at least 2 racial/ethnic groups including 7 that showed association with the expression of an annotated gene. CpGs annotated to CPT1A showed evidence of being influenced by triglycerides levels. DNA methylation levels of circulating leukocytes show robust and consistent association with blood lipid levels across multiple racial/ethnic groups.

Epigenome-wide association meta-analysis of DNA methylation with coffee and tea consumption

Coffee and tea are extensively consumed beverages worldwide which have received considerable attention regarding health. Intake of these beverages is consistently linked to, among others, reduced risk of diabetes and liver diseases; however, the mechanisms of action remain elusive. Epigenetics is suggested as a mechanism mediating the effects of dietary and lifestyle factors on disease onset. Here we report the results from epigenome-wide association studies (EWAS) on coffee and tea consumption in 15,789 participants of European and African-American ancestries from 15 cohorts. EWAS meta-analysis of coffee consumption reveals 11 CpGs surpassing the epigenome-wide significance threshold (P-value <1.1×10-7), which annotated to the AHRR, F2RL3, FLJ43663, HDAC4, GFI1 and PHGDH genes. Among them, cg14476101 is significantly associated with expression of the PHGDH and risk of fatty liver disease. Knockdown of PHGDH expression in liver cells shows a correlation with expression levels of genes associated with circulating lipids, suggesting a role of PHGDH in hepatic-lipid metabolism. EWAS meta-analysis on tea consumption reveals no significant association, only two CpGs annotated to CACNA1A and PRDM16 genes show suggestive association (P-value <5.0×10-6). These findings indicate that coffee-associated changes in DNA methylation levels may explain the mechanism of action of coffee consumption in conferring risk of diseases.

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.

Epigenome-Wide Association Study of Thyroid Function Traits Identifies Novel Associations of fT3 With KLF9 and DOT1L

CONTEXT

Circulating concentrations of free triiodothyronine (fT3), free thyroxine (fT4), and thyrotropin (TSH) are partly heritable traits. Recent studies have advanced knowledge of their genetic architecture. Epigenetic modifications, such as DNA methylation (DNAm), may be important in pituitary-thyroid axis regulation and action, but data are limited.

OBJECTIVE

To identify novel associations between fT3, fT4, and TSH and differentially methylated positions (DMPs) in the genome in subjects from 2 Australian cohorts.

METHOD

We performed an epigenome-wide association study (EWAS) of thyroid function parameters and DNAm using participants from: Brisbane Systems Genetics Study (median age 14.2 years, n = 563) and the Raine Study (median age 17.0 years, n = 863). Plasma fT3, fT4, and TSH were measured by immunoassay. DNAm levels in blood were assessed using Illumina HumanMethylation450 BeadChip arrays. Analyses employed generalized linear mixed models to test association between DNAm and thyroid function parameters. Data from the 2 cohorts were meta-analyzed.

RESULTS

We identified 2 DMPs with epigenome-wide significant (P < 2.4E-7) associations with TSH and 6 with fT3, including cg00049440 in KLF9 (P = 2.88E-10) and cg04173586 in DOT1L (P = 2.09E-16), both genes known to be induced by fT3. All DMPs had a positive association between DNAm and TSH and a negative association between DNAm and fT3. There were no DMPs significantly associated with fT4. We identified 23 differentially methylated regions associated with fT3, fT4, or TSH.

CONCLUSIONS

This study has demonstrated associations between blood-based DNAm and both fT3 and TSH. This may provide insight into mechanisms underlying thyroid hormone action and/or pituitary-thyroid axis function.

Novel DNA methylation signatures of tobacco smoking with trans-ethnic effects

BACKGROUND

Smoking remains one of the leading preventable causes of death. Smoking leaves a strong signature on the blood methylome as shown in multiple studies using the Infinium HumanMethylation450 BeadChip. Here, we explore novel blood methylation smoking signals on the Illumina MethylationEPIC BeadChip (EPIC) array, which also targets novel CpG-sites in enhancers.

METHOD

A smoking-methylation meta-analysis was carried out using EPIC DNA methylation profiles in 1407 blood samples from four UK population-based cohorts, including the MRC National Survey for Health and Development (NSHD) or 1946 British birth cohort, the National Child Development Study (NCDS) or 1958 birth cohort, the 1970 British Cohort Study (BCS70), and the TwinsUK cohort (TwinsUK). The overall discovery sample included 269 current, 497 former, and 643 never smokers. Replication was pursued in 3425 trans-ethnic samples, including 2325 American Indian individuals participating in the Strong Heart Study (SHS) in 1989-1991 and 1100 African-American participants in the Genetic Epidemiology Network of Arteriopathy Study (GENOA).

RESULTS

Altogether 952 CpG-sites in 500 genes were differentially methylated between smokers and never smokers after Bonferroni correction. There were 526 novel smoking-associated CpG-sites only profiled by the EPIC array, of which 486 (92%) replicated in a meta-analysis of the American Indian and African-American samples. Novel CpG sites mapped both to genes containing previously identified smoking-methylation signals and to 80 novel genes not previously linked to smoking, with the strongest novel signal in SLAMF7. Comparison of former versus never smokers identified that 37 of these sites were persistently differentially methylated after cessation, where 16 represented novel signals only profiled by the EPIC array. We observed a depletion of smoking-associated signals in CpG islands and an enrichment in enhancer regions, consistent with previous results.

CONCLUSION

This study identified novel smoking-associated signals as possible biomarkers of exposure to smoking and may help improve our understanding of smoking-related disease risk.

Large-scale association analyses identify host factors influencing human gut microbiome composition

To study the effect of host genetics on gut microbiome composition, the MiBioGen consortium curated and analyzed genome-wide genotypes and 16S fecal microbiome data from 18,340 individuals (24 cohorts). Microbial composition showed high variability across cohorts: only 9 of 410 genera were detected in more than 95% of samples. A genome-wide association study of host genetic variation regarding microbial taxa identified 31 loci affecting the microbiome at a genome-wide significant (P < 5 × 10-8) threshold. One locus, the lactase (LCT) gene locus, reached study-wide significance (genome-wide association study signal: P = 1.28 × 10-20), and it showed an age-dependent association with Bifidobacterium abundance. Other associations were suggestive (1.95 × 10-10 < P < 5 × 10-8) but enriched for taxa showing high heritability and for genes expressed in the intestine and brain. A phenome-wide association study and Mendelian randomization identified enrichment of microbiome trait loci in the metabolic, nutrition and environment domains and suggested the microbiome might have causal effects in ulcerative colitis and rheumatoid arthritis.

Equivalent DNA methylation variation between monozygotic co-twins and unrelated individuals reveals universal epigenetic inter-individual dissimilarity

BACKGROUND

Although the genomes of monozygotic twins are practically identical, their methylomes may evolve divergently throughout their lifetime as a consequence of factors such as the environment or aging. Particularly for young and healthy monozygotic twins, DNA methylation divergence, if any, may be restricted to stochastic processes occurring post-twinning during embryonic development and early life. However, to what extent such stochastic mechanisms can systematically provide a stable source of inter-individual epigenetic variation remains uncertain until now.

RESULTS

We enriched for inter-individual stochastic variation by using an equivalence testing-based statistical approach on whole blood methylation microarray data from healthy adolescent monozygotic twins. As a result, we identified 333 CpGs displaying similarly large methylation variation between monozygotic co-twins and unrelated individuals. Although their methylation variation surpasses measurement error and is stable in a short timescale, susceptibility to aging is apparent in the long term. Additionally, 46% of these CpGs were replicated in adipose tissue. The identified sites are significantly enriched at the clustered protocadherin loci, known for stochastic methylation in developing neurons. We also confirmed an enrichment in monozygotic twin DNA methylation discordance at these loci in whole genome bisulfite sequencing data from blood and adipose tissue.

CONCLUSIONS

We have isolated a component of stochastic methylation variation, distinct from genetic influence, measurement error, and epigenetic drift. Biomarkers enriched in this component may serve in the future as the basis for universal epigenetic fingerprinting, relevant for instance in the discrimination of monozygotic twin individuals in forensic applications, currently impossible with standard DNA profiling.

Serum metabolites reflecting gut microbiome alpha diversity predict type 2 diabetes

Type 2 diabetes (T2D) is associated with reduced gut microbiome diversity, although the cause is unclear. Metabolites generated by gut microbes also appear to be causative factors in T2D. We therefore searched for serum metabolites predictive of gut microbiome diversity in 1018 females from TwinsUK with concurrent metabolomic profiling and microbiome composition. We generated a Microbial Metabolites Diversity (MMD) score of six circulating metabolites that explained over 18% of the variance in microbiome alpha diversity. Moreover, the MMD score was associated with a significantly lower odds of prevalent (OR[95%CI] = 0.22[0.07;0.70], P = .01) and incident T2D (HR[95%CI] = 0.31[0.11,0.90], P = .03). We replicated our results in 1522 individuals from the ARIC study (prevalent T2D: OR[95%CI] = 0.79[0.64,0.96], P = .02, incident T2D: HR[95%CI] = 0.87[0.79,0.95], P = .003). The MMD score mediated 28%[15%,94%] of the total effect of gut microbiome on T2D after adjusting for confounders. Metabolites predicting higher microbiome diversity included 3-phenylpropionate(hydrocinnamate), indolepropionate, cinnamoylglycine and 5-alpha-pregnan-3beta,20 alpha-diol monosulfate(2) of which indolepropionate and phenylpropionate have already been linked to lower incidence of T2D. Metabolites correlating with lower microbial diversity included glutarate and imidazole propionate, of which the latter has been implicated in insulin resistance. Our results suggest that the effect of gut microbiome diversity on T2D is largely mediated by microbial metabolites, which might be modifiable by diet.

ACE2 expression in adipose tissue is associated with COVID-19 cardio-metabolic risk factors and cell type composition

COVID-19 severity has varied widely, with demographic and cardio-metabolic factors increasing risk of severe reactions to SARS-CoV-2 infection, but the underlying mechanisms for this remain uncertain. We investigated phenotypic and genetic factors associated with subcutaneous adipose tissue expression of Angiotensin I Converting Enzyme 2 ( ACE2 ), which has been shown to act as a receptor for SARS-CoV-2 cellular entry. In a meta-analysis of three independent studies including up to 1,471 participants, lower adipose tissue ACE2 expression was associated with adverse cardio-metabolic health indices including type 2 diabetes (T2D) and obesity status, higher serum fasting insulin and BMI, and lower serum HDL levels (P<5.32x10 -4 ). ACE2 expression levels were also associated with estimated proportions of cell types in adipose tissue; lower ACE2 expression was associated with a lower proportion of microvascular endothelial cells (P=4.25x10 -4 ) and higher macrophage proportion (P=2.74x10 -5 ), suggesting a link to inflammation. Despite an estimated heritability of 32%, we did not identify any proximal or distal genetic variants (eQTLs) associated with adipose tissue ACE2 expression. Our results demonstrate that at-risk individuals have lower background ACE2 levels in this highly relevant tissue. Further studies will be required to establish how this may contribute to increased COVID-19 severity.

Blood DNA methylation sites predict death risk in a longitudinal study of 12, 300 individuals

DNA methylation has fundamental roles in gene programming and aging that may help predict mortality. However, no large-scale study has investigated whether site-specific DNA methylation predicts all-cause mortality. We used the Illumina-HumanMethylation450-BeadChip to identify blood DNA methylation sites associated with all-cause mortality for 12, 300 participants in 12 Cohorts of the Heart and Aging Research in Genetic Epidemiology (CHARGE) Consortium. Over an average 10-year follow-up, there were 2,561 deaths across the cohorts. Nine sites mapping to three intergenic and six gene-specific regions were associated with mortality (P < 9.3x10^-7) independently of age and other mortality predictors. Six sites (cg14866069, cg23666362, cg20045320, cg07839457, cg07677157, cg09615688)—mapping respectively to BMPR1B, MIR1973, IFITM3, NLRC5, and two intergenic regions—were associated with reduced mortality risk. The remaining three sites (cg17086398, cg12619262, cg18424841)—mapping respectively to SERINC2, CHST12, and an intergenic region—were associated with increased mortality risk. DNA methylation at each site predicted 5%–15% of all deaths. We also assessed the causal association of those sites to age-related chronic diseases by using Mendelian randomization, identifying weak causal relationship between cg18424841 and cg09615688 with coronary heart disease. Of the nine sites, three (cg20045320, cg07839457, cg07677157) were associated with lower incidence of heart disease risk and two (cg20045320, cg07839457) with smoking and inflammation in prior CHARGE analyses. Methylation of cg20045320, cg07839457, and cg17086398 was associated with decreased expression of nearby genes (IFITM3, IRF, NLRC5, MT1, MT2, MARCKSL1) linked to immune responses and cardiometabolic diseases. These sites may serve as useful clinical tools for mortality risk assessment and preventative care.

Dysregulated Antibody, Natural Killer Cell and Immune Mediator Profiles in Autoimmune Thyroid Diseases

The pathogenesis of autoimmune thyroid diseases (AITD) is poorly understood and the association between different immune features and the germline variants involved in AITD are yet unclear. We previously observed systemic depletion of IgG core fucosylation and antennary α1,2 fucosylation in peripheral blood mononuclear cells in AITD, correlated with anti-thyroid peroxidase antibody (TPOAb) levels. Fucose depletion is known to potentiate strong antibody-mediated NK cell activation and enhanced target antigen-expressing cell killing. In autoimmunity, this may translate to autoantibody-mediated immune cell recruitment and attack of self-antigen expressing normal tissues. Hence, we investigated the crosstalk between immune cell traits, secreted proteins, genetic variants and the glycosylation patterns of serum IgG, in a multi-omic and cross-sectional study of 622 individuals from the TwinsUK cohort, 172 of whom were diagnosed with AITD. We observed associations between two genetic variants (rs505922 and rs687621), AITD status, the secretion of Desmoglein-2 protein, and the profile of two IgG N-glycan traits in AITD, but further studies need to be performed to better understand their crosstalk in AITD. On the other side, enhanced afucosylated IgG was positively associated with activatory CD335- CD314+ CD158b+ NK cell subsets. Increased levels of the apoptosis and inflammation markers Caspase-2 and Interleukin-1α positively associated with AITD. Two genetic variants associated with AITD, rs1521 and rs3094228, were also associated with altered expression of the thyrocyte-expressed ligands known to recognize the NK cell immunoreceptors CD314 and CD158b. Our analyses reveal a combination of heightened Fc-active IgG antibodies, effector cells, cytokines and apoptotic signals in AITD, and AITD genetic variants associated with altered expression of thyrocyte-expressed ligands to NK cell immunoreceptors. Together, TPOAb responses, dysregulated immune features, germline variants associated with immunoactivity profiles, are consistent with a positive autoreactive antibody-dependent NK cell-mediated immune response likely drawn to the thyroid gland in AITD.

DNA Methylation Age and Physical and Cognitive Aging

BACKGROUND

DNA methylation (DNAm) age acceleration (AgeAccel) has been shown to be predictive of all-cause mortality but it is unclear what functional aspect(s) of aging it captures. We examine associations between four measures of AgeAccel in adults aged 45-87 years and physical and cognitive performance and their age-related decline.

METHODS

AgeAccelHannum, AgeAccelHorvath, AgeAccelPheno, and AgeAccelGrim were calculated in the Medical Research Council National Survey of Health and Development (NSHD), National Child Development Study (NCDS) and TwinsUK. Three measures of physical (grip strength, chair rise speed, and forced expiratory volume in one second [FEV1]) and two measures of cognitive (episodic memory and mental speed) performance were assessed.

RESULTS

AgeAccelPheno and AgeAccelGrim, but not AgeAccelHannum and AgeAccelHorvath were related to performance in random effects meta-analyses (n = 1,388-1,685). For example, a 1-year increase in AgeAccelPheno or AgeAccelGrim was associated with a 0.01 mL (95% confidence interval [CI]: 0.01, 0.02) or 0.03 mL (95% CI: 0.01, 0.05) lower mean FEV1 respectively. In NSHD, AgeAccelPheno and AgeAccelGrim at 53 years were associated with age-related decline in performance between 53 and 69 years as tested by linear mixed models (p < .05). In a subset of NSHD participants (n = 482), there was little evidence that change in any AgeAccel measure was associated with change in performance conditional on baseline performance.

CONCLUSIONS

We found little evidence to support associations between the first generation of DNAm-based biomarkers of aging and age-related physical or cognitive performance in midlife to early old age. However, there was evidence that the second generation biomarkers, AgeAccelPheno and AgeAccelGrim, could act as makers of an individual's healthspan as proposed.

DNA methylation aging clocks: challenges and recommendations

Epigenetic clocks comprise a set of CpG sites whose DNA methylation levels measure subject age. These clocks are acknowledged as a highly accurate molecular correlate of chronological age in humans and other vertebrates. Also, extensive research is aimed at their potential to quantify biological aging rates and test longevity or rejuvenating interventions. Here, we discuss key challenges to understand clock mechanisms and biomarker utility. This requires dissecting the drivers and regulators of age-related changes in single-cell, tissue- and disease-specific models, as well as exploring other epigenomic marks, longitudinal and diverse population studies, and non-human models. We also highlight important ethical issues in forensic age determination and predicting the trajectory of biological aging in an individual.

DNA methylation-based estimator of telomere length

Telomere length (TL) is associated with several aging-related diseases. Here, we present a DNA methylation estimator of TL (DNAmTL) based on 140 CpGs. Leukocyte DNAmTL is applicable across the entire age spectrum and is more strongly associated with age than measured leukocyte TL (LTL) (r ~-0.75 for DNAmTL versus r ~ -0.35 for LTL). Leukocyte DNAmTL outperforms LTL in predicting: i) time-to-death (p=2.5E-20), ii) time-to-coronary heart disease (p=6.6E-5), iii) time-to-congestive heart failure (p=3.5E-6), and iv) association with smoking history (p=1.21E-17). These associations are further validated in large scale methylation data (n=10k samples) from the Framingham Heart Study, Women's Health Initiative, Jackson Heart Study, InChianti, Lothian Birth Cohorts, Twins UK, and Bogalusa Heart Study. Leukocyte DNAmTL is also associated with measures of physical fitness/functioning (p=0.029), age-at-menopause (p=0.039), dietary variables (omega 3, fish, vegetable intake), educational attainment (p=3.3E-8) and income (p=3.1E-5). Experiments in cultured somatic cells show that DNAmTL dynamics reflect in part cell replication rather than TL per se. DNAmTL is not only an epigenetic biomarker of replicative history of cells, but a useful marker of age-related pathologies that are associated with it.

Association of dietary folate and vitamin B-12 intake with genome-wide DNA methylation in blood: a large-scale epigenome-wide association analysis in 5841 individuals

BACKGROUND

Folate and vitamin B-12 are essential micronutrients involved in the donation of methyl groups in cellular metabolism. However, associations between intake of these nutrients and genome-wide DNA methylation levels have not been studied comprehensively in humans.

OBJECTIVE

The aim of this study was to assess whether folate and/or vitamin B-12 intake are asssociated with genome-wide changes in DNA methylation in leukocytes.

METHODS

A large-scale epigenome-wide association study of folate and vitamin B-12 intake was performed on DNA from 5841 participants from 10 cohorts using Illumina 450k arrays. Folate and vitamin B-12 intakes were calculated from food-frequency questionnaires (FFQs). Continuous and categorical (low compared with high intake) linear regression mixed models were applied per cohort, controlling for confounders. A meta-analysis was performed to identify significant differentially methylated positions (DMPs) and regions (DMRs), and a pathway analysis was performed on the DMR annotated genes.

RESULTS

The categorical model resulted in 6 DMPs, which are all negatively associated with folate intake, annotated to FAM64A, WRAP73, FRMD8, CUX1, and LCN8 genes, which have a role in cellular processes including centrosome localization, cell proliferation, and tumorigenesis. Regional analysis showed 74 folate-associated DMRs, of which 73 were negatively associated with folate intake. The most significant folate-associated DMR was a 400-base pair (bp) spanning region annotated to the LGALS3BP gene. In the categorical model, vitamin B-12 intake was associated with 29 DMRs annotated to 48 genes, of which the most significant was a 1100-bp spanning region annotated to the calcium-binding tyrosine phosphorylation-regulated gene (CABYR). Vitamin B-12 intake was not associated with DMPs.

CONCLUSIONS

We identified novel epigenetic loci that are associated with folate and vitamin B-12 intake. Interestingly, we found a negative association between folate and DNA methylation. Replication of these methylation loci is necessary in future studies.