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.

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 Caenorhabditis elegans following knockdown of the BSCL2 homolog seip-1 on a daf-16/forkhead box, class O 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 gene when comparing CGL2 patients versus 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.

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.

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 .

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.

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.

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.

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.

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.

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.

Author Correction: Bayesian reassessment of the epigenetic architecture of complex traits.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Bayesian reassessment of the epigenetic architecture of complex traits.

Linking epigenetic marks to clinical outcomes improves insight into molecular processes, disease prediction, and therapeutic target identification. Here, a statistical approach is presented to infer the epigenetic architecture of complex disease, determine the variation captured by epigenetic effects, and estimate phenotype-epigenetic probe associations jointly. Implicitly adjusting for probe correlations, data structure (cell-count or relatedness), and single-nucleotide polymorphism (SNP) marker effects, improves association estimates and in 9,448 individuals, 75.7% (95% CI 71.70-79.3) of body mass index (BMI) variation and 45.6% (95% CI 37.3-51.9) of cigarette consumption variation was captured by whole blood methylation array data. Pathway-linked probes of blood cholesterol, lipid transport and sterol metabolism for BMI, and xenobiotic stimuli response for smoking, showed >1.5 times larger associations with >95% posterior inclusion probability. Prediction accuracy improved by 28.7% for BMI and 10.2% for smoking over a LASSO model, with age-, and tissue-specificity, implying associations are a phenotypic consequence rather than causal.