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

Causal relationship between iron status and preeclampsia-eclampsia: a Mendelian randomization analysis

BACKGROUND

Preeclampsia/eclampsia is a severe pregnancy-related disorder associated with hypertension and organ damage. While observational studies have suggested a link between maternal iron status and preeclampsia/eclampsia, the causal relationship remains unclear. The aim of this study was to investigate the genetic causality between iron status and preeclampsia/eclampsia using large-scale genome-wide association study (GWAS) summary data and Mendelian randomization (MR) analysis.

METHODS

Summary data for the GWAS on preeclampsia/eclampsia and genetic markers related to iron status were obtained from the FinnGen Consortium and the IEU genetic databases. The "TwoSampleMR" software package in R was employed to test the genetic causality between these markers and preeclampsia/eclampsia. The inverse variance weighted (IVW) method was primarily used for MR analysis. Heterogeneity, horizontal pleiotropy, and potential outliers were evaluated for the MR analysis results.

RESULTS

The random-effects IVW results showed that ferritin (OR = 1.11, 95% CI: .89-1.38, p = .341), serum iron (OR = .90, 95% CI: .75-1.09, p = .275), TIBC (OR = .98, 95% CI: .89-1.07, p = .613), and TSAT (OR = .94, 95% CI: .83-1.07, p = .354) have no genetic causal relationship with preeclampsia/eclampsia. There was no evidence of heterogeneity, horizontal pleiotropy, or possible outliers in our MR analysis (p > .05).

CONCLUSIONS

Our study did not detect a genetic causal relationship between iron status and preeclampsia/eclampsia. Nonetheless, this does not rule out a relationship between the two at other mechanistic levels.

Epigenetic Age Mediates the Association of Life's Essential 8 With Cardiovascular Disease and Mortality

BACKGROUND

Life's Essential 8 (LE8) is an enhanced metric for cardiovascular health. The interrelations among LE8, biomarkers of aging, and disease risks are unclear.

METHODS AND RESULTS

LE8 score was calculated for 5682 Framingham Heart Study participants. We implemented 4 DNA methylation-based epigenetic age biomarkers, with older epigenetic age hypothesized to represent faster biological aging, and examined whether these biomarkers mediated the associations between the LE8 score and cardiovascular disease (CVD), CVD-specific mortality, and all-cause mortality. We found that a 1 SD increase in the LE8 score was associated with a 35% (95% CI, 27-41; P=1.8E-15) lower risk of incident CVD, a 36% (95% CI, 24-47; P=7E-7) lower risk of CVD-specific mortality, and a 29% (95% CI, 22-35; P=7E-15) lower risk of all-cause mortality. These associations were partly mediated by epigenetic age biomarkers, particularly the GrimAge and the DunedinPACE scores. The potential mediation effects by epigenetic age biomarkers tended to be more profound in participants with higher genetic risk for older epigenetic age, compared with those with lower genetic risk. For example, in participants with higher GrimAge polygenic scores (greater than median), the mean proportion of mediation was 39%, 39%, and 78% for the association of the LE8 score with incident CVD, CVD-specific mortality, and all-cause mortality, respectively. No significant mediation was observed in participants with lower GrimAge polygenic score.

CONCLUSIONS

DNA methylation-based epigenetic age scores mediate the associations between the LE8 score and incident CVD, CVD-specific mortality, and all-cause mortality, particularly in individuals with higher genetic predisposition for older epigenetic age.

Causal Associations of DNA Methylation and Cardiovascular Disease: A Two-Sample Mendelian Randomization Study

Background

There is growing evidence that concentrations of DNA methylation are associated with cardiovascular disease; however, it is unclear whether this association reflects a causal relationship.

Methods

We utilized a two-sample Mendelian randomization (MR) approach to investigate whether DNA methylation can affect the risk of developing cardiovascular disease in human life. We primarily performed the inverse variance weighted (IVW) method to analyze the causal effect of DNA methylation on multiple cardiovascular diseases. Additionally, to ensure the robustness of our findings, we conducted several sensitivity analyses using alternative methodologies. These analysis methods included maximum likelihood, MR-Egger regression, weighted median method, and weighted model methods.

Results

Inverse variance weighted estimates suggested that an SD increase in DNA methylation Hannum age acceleration exposure increased the risk of cardiac arrhythmias (OR = 1.03, 95% CI 1.00-1.05, p = 0.0290) and atrial fibrillation (OR = 1.03, 95% CI 1.00-1.05, p = 0.0022). We also found that an SD increase in DNA methylation PhenoAge acceleration exposure increased the risk of heart failure (OR = 1.01, 95% CI 1.00-1.03, p = 0.0362). Exposure to DNA methylation-estimated granulocyte proportions was found to increase the risk of hypertension (OR = 1.00, 95% CI 1.00-1.0001, p = 0.0291). Exposure to DNA methylation-estimated plasminogen activator inhibitor-1 levels was found to increase the risk of heart failure (OR = 1.00, 95% CI 1.00-1.00, p = 0.0215).

Conclusion

This study reveals a causal relationship between DNA methylation and CVD. Exposed to high levels of DNA methylation Hannum age acceleration inhabitants with an increased risk of cardiac arrhythmias and atrial fibrillation. DNA methylation PhenoAge acceleration levels exposure levels were positively associated with the increased risk of developing heart failure. This has important implications for the prevention of cardiovascular diseases.

The effect of epigenetic aging on neurodegenerative diseases: a Mendelian randomization study

Background

Aging has always been considered as a risk factor for neurodegenerative diseases, but there are individual differences and its mechanism is not yet clear. Epigenetics may unveil the relationship between aging and neurodegenerative diseases.

Methods

Our study employed a bidirectional two-sample Mendelian randomization (MR) design to assess the potential causal association between epigenetic aging and neurodegenerative diseases. We utilized publicly available summary datasets from several genome-wide association studies (GWAS). Our investigation focused on multiple measures of epigenetic age as potential exposures and outcomes, while the occurrence of neurodegenerative diseases served as potential exposures and outcomes. Sensitivity analyses confirmed the accuracy of the results.

Results

The results show a significant decrease in risk of Parkinson's disease with GrimAge (OR = 0.8862, 95% CI 0.7914-0.9924, p = 0.03638). Additionally, we identified that HannumAge was linked to an increased risk of Multiple Sclerosis (OR = 1.0707, 95% CI 1.0056-1.1401, p = 0.03295). Furthermore, we also found that estimated plasminogen activator inhibitor-1(PAI-1) levels demonstrated an increased risk for Alzheimer's disease (OR = 1.0001, 95% CI 1.0000-1.0002, p = 0.04425). Beyond that, we did not observe any causal associations between epigenetic age and neurodegenerative diseases risk.

Conclusion

The findings firstly provide evidence for causal association of epigenetic aging and neurodegenerative diseases. Exploring neurodegenerative diseases from an epigenetic perspective may contribute to diagnosis, prognosis, and treatment of neurodegenerative diseases.

Editorial: Can Improving Youth Mental Health Reduce Mortality?

It is well established that mental health conditions, including substance use disorders, are associated with premature mortality. A meta-analysis1 has demonstrated that this association holds across a range of diagnoses. Although the effect is stronger for schizophrenia, depression and anxiety contribute to more deaths overall because of their high prevalence rates. Moreover, more than two-thirds of associated deaths were explained by natural causes.1 The next logical questions, then, are as follows: which mechanisms underlie this association, and can they can be mitigated? In the current issue of JAACAP, Clark et al.2 aim to tie mental health symptoms and substance use to the acceleration of biological aging.

Exploring machine learning strategies for predicting cardiovascular disease risk factors from multi-omic data

BACKGROUND

Machine learning (ML) classifiers are increasingly used for predicting cardiovascular disease (CVD) and related risk factors using omics data, although these outcomes often exhibit categorical nature and class imbalances. However, little is known about which ML classifier, omics data, or upstream dimension reduction strategy has the strongest influence on prediction quality in such settings. Our study aimed to illustrate and compare different machine learning strategies to predict CVD risk factors under different scenarios.

METHODS

We compared the use of six ML classifiers in predicting CVD risk factors using blood-derived metabolomics, epigenetics and transcriptomics data. Upstream omic dimension reduction was performed using either unsupervised or semi-supervised autoencoders, whose downstream ML classifier performance we compared. CVD risk factors included systolic and diastolic blood pressure measurements and ultrasound-based biomarkers of left ventricular diastolic dysfunction (LVDD; E/e' ratio, E/A ratio, LAVI) collected from 1,249 Finnish participants, of which 80% were used for model fitting. We predicted individuals with low, high or average levels of CVD risk factors, the latter class being the most common. We constructed multi-omic predictions using a meta-learner that weighted single-omic predictions. Model performance comparisons were based on the F1 score. Finally, we investigated whether learned omic representations from pre-trained semi-supervised autoencoders could improve outcome prediction in an external cohort using transfer learning.

RESULTS

Depending on the ML classifier or omic used, the quality of single-omic predictions varied. Multi-omics predictions outperformed single-omics predictions in most cases, particularly in the prediction of individuals with high or low CVD risk factor levels. Semi-supervised autoencoders improved downstream predictions compared to the use of unsupervised autoencoders. In addition, median gains in Area Under the Curve by transfer learning compared to modelling from scratch ranged from 0.09 to 0.14 and 0.07 to 0.11 units for transcriptomic and metabolomic data, respectively.

CONCLUSIONS

By illustrating the use of different machine learning strategies in different scenarios, our study provides a platform for researchers to evaluate how the choice of omics, ML classifiers, and dimension reduction can influence the quality of CVD risk factor predictions.

Heterogeneity of aging and mortality risk among individuals with hypertension: Insights from phenotypic age and phenotypic age acceleration

OBJECTIVES

Hypertension, a key contributor to mortality, is impacted by biological aging. We investigated the relationship between novel biological aging metrics - Phenotypic Age (PA) and Phenotypic Age Acceleration (PAA) - and mortality in individuals with hypertension, exploring the mediating effects of arterial stiffness (estimated Pulse Wave Velocity, ePWV), and Heart/Vascular Age (HVA).

METHODS

Using data from 62,160 National Health and Nutrition Examination Survey (NHANES) participants (1999-2010), we selected 4,228 individuals with hypertension and computed PA, PAA, HVA, and ePWV. Weighted, multivariable Cox regression analysis yielded Hazard Ratios (HRs) relating PA, PAA to mortality, and mediation roles of ePWV, PAA, HVA were evaluated. Mendelian randomization (MR) analysis was employed to investigate causality between genetically inferred PAA and hypertension.

RESULTS

Over a 12-year median follow-up, PA and PAA were tied to increased mortality risks in individuals with hypertension. All-cause mortality hazard ratios per 10-year PA and PAA increments were 1.96 (95% CI, 1.81-2.11) and 1.67 (95% CI, 1.52-1.85), respectively. Cardiovascular mortality HRs were 2.32 (95% CI, 1.97-2.73) and 1.93 (95% CI, 1.65-2.26) for PA and PAA, respectively. ePWV, PAA, and HVA mediated 42%, 30.3%, and 6.9% of PA's impact on mortality, respectively. Mendelian randomization highlighted a causal link between PAA genetics and hypertension (OR = 1.002; 95% CI, 1.000-1.003).

CONCLUSION

PA and PAA, enhancing cardiovascular risk scores by integrating diverse biomarkers, offer vital insights for aging and mortality evaluation in individuals with hypertension, suggesting avenues for intensified aging mitigation and cardiovascular issue prevention. Validations in varied populations and explorations of underlying mechanisms are warranted.

The causal impact of gut microbiota on circulating adipokine concentrations: a two-sample Mendelian randomization study

PURPOSE

Evidence from previous experimental and observational research demonstrates that the gut microbiota is related to circulating adipokine concentrations. Nevertheless, the debate as to whether gut microbiome composition causally influences circulating adipokine concentrations remains unresolved. This study aimed to take an essential step in elucidating this issue.

METHODS

We used two-sample Mendelian randomization (MR) to causally analyze genetic variation statistics for gut microbiota and four adipokines (including adiponectin, leptin, soluble leptin receptor [sOB-R], and plasminogen activator inhibitor-1 [PAI-1]) from large-scale genome-wide association studies (GWAS) datasets. A range of sensitivity analyses was also conducted to assess the stability and reliability of the results.

RESULTS

The composite results of the MR and sensitivity analyses revealed 22 significant causal associations. In particular, there is a suggestive causality between the family Clostridiaceae1 (IVW: β = 0.063, P = 0.034), the genus Butyrivibrio (IVW: β = 0.029, P = 0.031), and the family Alcaligenaceae (IVW: β=-0.070, P = 0.014) and adiponectin. Stronger causal effects with leptin were found for the genus Enterorhabdus (IVW: β=-0.073, P = 0.038) and the genus Lachnospiraceae (NK4A136 group) (IVW: β=-0.076, P = 0.01). Eight candidate bacterial groups were found to be associated with sOB-R, with the phylum Firmicutes (IVW: β = 0.235, P = 0.03) and the order Clostridiales (IVW: β = 0.267, P = 0.028) being of more interest. In addition, the genus Roseburia (IVW: β = 0.953, P = 0.022) and the order Lactobacillales (IVW: β=-0.806, P = 0.042) were suggestive of an association with PAI-1.

CONCLUSION

This study reveals a causal relationship between the gut microbiota and circulating adipokines and may help to offer novel insights into the prevention of abnormal concentrations of circulating adipokines and obesity-related diseases.

Analysis of epigenetic clocks links yoga, sleep, education, reduced meat intake, coffee, and a SOCS2 gene variant to slower epigenetic aging

DNA methylation (DNAm) clocks hold promise for measuring biological age, useful for guiding clinical interventions and forensic identification. This study compared the commonly used DNAm clocks, using DNA methylation and SNP data generated from nearly 1000 human blood or buccal swab samples. We evaluated different preprocessing methods for age estimation, investigated the association of epigenetic age acceleration (EAA) with various lifestyle and sociodemographic factors, and undertook a series of novel genome-wide association analyses for different EAA measures to find associated genetic variants. Our results highlighted the Skin&Blood clock with ssNoob normalization as the most accurate predictor of chronological age. We provided novel evidence for an association between the practice of yoga and a reduction in the pace of aging (DunedinPACE). Increased sleep and physical activity were associated with lower mortality risk score (MRS) in our dataset. University degree, vegetable consumption, and coffee intake were associated with reduced levels of epigenetic aging, whereas smoking, higher BMI, meat consumption, and manual occupation correlated well with faster epigenetic aging, with FitAge, GrimAge, and DunedinPACE clocks showing the most robust associations. In addition, we found a novel association signal for SOCS2 rs73218878 (p = 2.87 × 10-8) and accelerated GrimAge. Our study emphasizes the importance of an optimized DNAm analysis workflow for accurate estimation of epigenetic age, which may influence downstream analyses. The results support the influence of genetic background on EAA. The associated SOCS2 is a member of the suppressor of cytokine signaling family known for its role in human longevity. The reported association between various risk factors and EAA has practical implications for the development of health programs to improve quality of life and reduce premature mortality associated with age-related diseases.

Causal relationship between gestational diabetes and preeclampsia: A bidirectional mendelian randomization analysis

AIMS

The study aimed to explore the potential causal link between gestational diabetes mellitus (GDM) and preeclampsia (PE) using a bidirectional mendelian randomization (MR) analysis.

MATERIALS

We conducted a bidirectional MR analysis to investigate the causal relationship between GDM and PE. Data from public genome-wide association studies (GWAS) for GDM and PE were obtained from the FinnGen consortium. Various MR methods were employed, including inverse-variance weighted (IVW), MR-Egger, and sensitivity analyses. Additionally, a knowledge-based approach identified genes underlying this potential connection.

RESULTS

The IVW method revealed a lack of significant association between GDM and PE (OR: 1.04, 95 % CI: 0.96-1.14; p = 0.275). Conversely, IVW analysis indicated a causal connection from PE to GDM (OR: 1.14, 95 % CI: 1.06-1.23; p < 0.001). Molecular pathway analysis identified 20 key genes, including ASAP2, central to the PE-GDM relationship. Tissue enrichment analysis showed pertinent gene expression in significant tissues. Moreover, lower ASAP2 expression was detected in PE patients' placentas.

CONCLUSIONS

Our bidirectional MR analysis offers evidence supporting a causal link between PE and GDM, elucidating their interconnected pathogenesis. Genetic and knowledge-based insights facilitate a deeper comprehension of these complex pregnancy complications.

A genome-wide association study of neutrophil count in individuals associated to an African continental ancestry group facilitates studies of malaria pathogenesis

BACKGROUND

'Benign ethnic neutropenia' (BEN) is a heritable condition characterized by lower neutrophil counts, predominantly observed in individuals of African ancestry, and the genetic basis of BEN remains a subject of extensive research. In this study, we aimed to dissect the genetic architecture underlying neutrophil count variation through a linear-mixed model genome-wide association study (GWAS) in a population of African ancestry (N = 5976). Malaria caused by P. falciparum imposes a tremendous public health burden on people living in sub-Saharan Africa. Individuals living in malaria endemic regions often have a reduced circulating neutrophil count due to BEN, raising the possibility that reduced neutrophil counts modulate severity of malaria in susceptible populations. As a follow-up, we tested this hypothesis by conducting a Mendelian randomization (MR) analysis of neutrophil counts on severe malaria (MalariaGEN, N = 17,056).

RESULTS

We carried out a GWAS of neutrophil count in individuals associated to an African continental ancestry group within UK Biobank, identifying 73 loci (r2 = 0.1) and 10 index SNPs (GCTA-COJO loci) associated with neutrophil count, including previously unknown rare loci regulating neutrophil count in a non-European population. BOLT-LMM was reliable when conducted in a non-European population, and additional covariates added to the model did not largely alter the results of the top loci or index SNPs. The two-sample bi-directional MR analysis between neutrophil count and severe malaria showed the greatest evidence for an effect between neutrophil count and severe anaemia, although the confidence intervals crossed the null.

CONCLUSION

Our GWAS of neutrophil count revealed unique loci present in individuals of African ancestry. We note that a small sample-size reduced our power to identify variants with low allele frequencies and/or low effect sizes in our GWAS. Our work highlights the need for conducting large-scale biobank studies in Africa and for further exploring the link between neutrophils and severe malaria.

Epigenetic age acceleration and risk of aortic valve stenosis: a bidirectional Mendelian randomization study

BACKGROUND

Aortic valve stenosis (AVS) is the most prevalent cardiac valve lesion in developed countries, and pathogenesis is closely related to aging. DNA methylation-based epigenetic clock is now recognized as highly accurate predictor of the aging process and associated health outcomes. This study aimed to explore the causal relationship between epigenetic clock and AVS by conducting a bidirectional Mendelian randomization (MR) analysis.

METHODS

Summary genome-wide association study statistics of epigenetic clocks (HannumAge, HorvathAge, PhenoAge, and GrimAge) and AVS were obtained and assessed for significant instrumental variables from Edinburgh DataShare (n = 34,710) and FinnGen biobank (cases = 9870 and controls = 402,311). The causal association between epigenetic clock and AVS was evaluated using inverse variance weighted (IVW), weighted median (WM), and MR-Egger methods. Multiple analyses (heterogeneity analysis, pleiotropy analysis, and sensitivity analysis) were performed for quality control assessment.

RESULTS

The MR analysis showed that the epigenetic age acceleration of HorvathAge and PhenoAge was associated with an increased risk of AVS (HorvathAge: OR = 1.043, P = 0.016 by IVW, OR = 1.058, P = 0.018 by WM; PhenoAge: OR = 1.058, P = 0.005 by IVW, OR = 1.053, P = 0.039 by WM). Quality control assessment proved our findings were reliable and robust. However, there was a lack of evidence supporting a causal link from AVS to epigenetic aging.

CONCLUSION

The present MR analysis unveiled a causal association between epigenetic clocks, especially HorvathAge and PhenoAge, with AVS. Further research is required to elucidate the underlying mechanisms and develop strategies for potential interventions.

The relationship between 11 different polygenic longevity scores, parental lifespan, and disease diagnosis in the UK Biobank

Large-scale genome-wide association studies (GWAS) strongly suggest that most traits and diseases have a polygenic component. This observation has motivated the development of disease-specific "polygenic scores (PGS)" that are weighted sums of the effects of disease-associated variants identified from GWAS that correlate with an individual's likelihood of expressing a specific phenotype. Although most GWAS have been pursued on disease traits, leading to the creation of refined "Polygenic Risk Scores" (PRS) that quantify risk to diseases, many GWAS have also been pursued on extreme human longevity, general fitness, health span, and other health-positive traits. These GWAS have discovered many genetic variants seemingly protective from disease and are often different from disease-associated variants (i.e., they are not just alternative alleles at disease-associated loci) and suggest that many health-positive traits also have a polygenic basis. This observation has led to an interest in "polygenic longevity scores (PLS)" that quantify the "risk" or genetic predisposition of an individual towards health. We derived 11 different PLS from 4 different available GWAS on lifespan and then investigated the properties of these PLS using data from the UK Biobank (UKB). Tests of association between the PLS and population structure, parental lifespan, and several cancerous and non-cancerous diseases, including death from COVID-19, were performed. Based on the results of our analyses, we argue that PLS are made up of variants not only robustly associated with parental lifespan, but that also contribute to the genetic architecture of disease susceptibility, morbidity, and mortality.

Smoking and BMI mediate the causal effect of education on lower back pain: observational and Mendelian randomization analyses

Objective

Low back pain (LBP) has been associated with education in previous observational studies, but the causality remains unclear. This study aims to assess the impact of education on LBP and to explore mediation by multiple lifestyle factors.

Design

Univariable Mendelian randomization (MR) was performed to examine the overall effect of education on LBP. Subsequently, multivariable MR was conducted to assess both the direct effect of education on LBP and the influence of potential mediators. Indirect effects were estimated using either the coefficient product method or the difference method, and the proportion of mediation was calculated by dividing the indirect effect by the total effect. The observational study utilized data from the NHANES database collected between 1999 and 2004, and included 15,580 participants aged 20 years and above.

Results

Increasing education by 4.2 years leads to a 48% reduction in the risk of LBP (OR=0.52; 95% CI: 0.46 to 0.59). Compared to individuals with less than a high school education, those with education beyond high school have a 28% lower risk of LBP (OR=0.72; 95% CI: 0.63 to 0.83). In the MR study, smoking accounts for 12.8% (95% CI: 1.04% to 20.8%) of the total effect, while BMI accounts for 5.9% (95% CI: 2.99% to 8.55%). The combined mediation effect of smoking and BMI is 27.6% (95% CI: 23.99% to 32.7%). In the NHANES study, only smoking exhibits a mediating effect, accounting for 34.3% (95% CI: 21.07% to 41.65%) of the effect, while BMI does not demonstrate a mediating role.

Conclusions

Higher levels of education provide a protective effect against the risk of LBP. Additionally, implementing interventions to reduce smoking and promote weight loss among individuals with lower levels of education can also decrease this risk.

Genetic relationship between ageing and coronary heart disease: a Mendelian randomization study

PURPOSE

Genetic relationship between ageing and coronary heart disease has not been well investigated. The aim of the study was to explore the association of several ageing biomarkers with the risk of several types of coronary heart disease using the Mendelian randomization approach.

METHODS

Summary data for telomere length, four epigenetic clocks (such as intrinsic epigenetic age acceleration), four types of coronary heart disease (such as myocardial infarction) were collected from the most updated and available genome-wide association studies. Instrumental variables were extracted from the exposure-related summary data according to correlation, independence and exclusivity assumptions. Three Mendelian randomization methods (such as inverse variance weighted) were used for causal inference. Four sensitivity analyses (such as MR-Egger intercept) were performed to prevent horizontal pleiotropy.

RESULTS

Inverse variance weighted reported that longer telomere length was related to the lower risk of myocardial infarction, angina pectoris, unstable angina pectoris and coronary atherosclerosis (P = 8.840e-11, P = 9.830e-04, P = 1.539e-05, P = 2.607e-09). Inverse variance weighted also reported that four epigenetic clocks might be not implicated in the risk of these coronary heart diseases. Furthermore, there was not enough evidence to confirm the effect of coronary heart disease on these ageing biomarkers.

CONCLUSION

Longer telomere length, but not the epigenetic clock changes, genetically decreased the risk of coronary heart disease. Considering that telomere length and epigenetic clocks were two independent ageing biomarkers, the correlation between ageing and coronary heart disease might be redefined at the genetic level.

Age-associated cortical similarity networks correlate with cell type-specific transcriptional signatures

Human brain structure shows heterogeneous patterns of change across adults aging and is associated with cognition. However, the relationship between cortical structural changes during aging and gene transcription signatures remains unclear. Here, using structural magnetic resonance imaging data of two separate cohorts of healthy participants from the Cambridge Centre for Aging and Neuroscience (n = 454, 18-87 years) and Dallas Lifespan Brain Study (n = 304, 20-89 years) and a transcriptome dataset, we investigated the link between cortical morphometric similarity network and brain-wide gene transcription. In two cohorts, we found reproducible morphometric similarity network change patterns of decreased morphological similarity with age in cognitive related areas (mainly located in superior frontal and temporal cortices), and increased morphological similarity in sensorimotor related areas (postcentral and lateral occipital cortices). Changes in morphometric similarity network showed significant spatial correlation with the expression of age-related genes that enriched to synaptic-related biological processes, synaptic abnormalities likely accounting for cognitive decline. Transcription changes in astrocytes, microglia, and neuronal cells interpreted most of the age-related morphometric similarity network changes, which suggest potential intervention and therapeutic targets for cognitive decline. Taken together, by linking gene transcription signatures to cortical morphometric similarity network, our findings might provide molecular and cellular substrates for cortical structural changes related to cognitive decline across adults aging.

Assessing causal associations between neurodegenerative diseases and neurological tumors with biological aging: a bidirectional Mendelian randomization study

Background

Aging is a significant risk factor for many neurodegenerative diseases and neurological tumors. Previous studies indicate that the frailty index, facial aging, telomere length (TL), and epigenetic aging clock acceleration are commonly used biological aging proxy indicators. This study aims to comprehensively explore potential relationships between biological aging and neurodegenerative diseases and neurological tumors by integrating various biological aging proxy indicators, employing Mendelian randomization (MR) analysis.

Methods

Two-sample bidirectional MR analyses were conducted using genome-wide association study (GWAS) data. Summary statistics for various neurodegenerative diseases and neurological tumors, along with biological aging proxy indicators, were obtained from extensive meta-analyses of GWAS. Genetic single-nucleotide polymorphisms (SNPs) associated with the exposures were used as instrumental variables, assessing causal relationships between three neurodegenerative diseases (Alzheimer's disease, Parkinson's disease, amyotrophic lateral sclerosis), two benign neurological tumors (vestibular schwannoma and meningioma), one malignant neurological tumor (glioma), and four biological aging indicators (frailty index, facial aging, TL, and epigenetic aging clock acceleration). Sensitivity analyses were also performed.

Results

Our analysis revealed that genetically predicted longer TL reduces the risk of Alzheimer's disease but increases the risk of vestibular schwannoma and glioma (All Glioma, GBM, non-GBM). In addition, there is a suggestive causal relationship between some diseases (PD and GBM) and DNA methylation GrimAge acceleration. Causal relationships between biological aging proxy indicators and other neurodegenerative diseases and neurological tumors were not observed.

Conclusion

Building upon prior investigations into the causal relationships between telomeres and neurodegenerative diseases and neurological tumors, our study validates these findings using larger GWAS data and demonstrates, for the first time, that Parkinson's disease and GBM may promote epigenetic age acceleration. Our research provides new insights and evidence into the causal relationships between biological aging and the risk of neurodegenerative diseases and neurological tumors.

Somatic mutation as an explanation for epigenetic aging

DNA methylation marks have recently been used to build models known as "epigenetic clocks" which predict calendar age. As methylation of cytosine promotes C-to-T mutations, we hypothesized that the methylation changes observed with age should reflect the accrual of somatic mutations, and the two should yield analogous aging estimates. In analysis of multimodal data from 9,331 human individuals, we find that CpG mutations indeed coincide with changes in methylation, not only at the mutated site but also with pervasive remodeling of the methylome out to ±10 kilobases. This one-to-many mapping enables mutation-based predictions of age that agree with epigenetic clocks, including which individuals are aging faster or slower than expected. Moreover, genomic loci where mutations accumulate with age also tend to have methylation patterns that are especially predictive of age. These results suggest a close coupling between the accumulation of sporadic somatic mutations and the widespread changes in methylation observed over the course of life.

ExplaiNAble BioLogical Age (ENABL Age): an artificial intelligence framework for interpretable biological age

BACKGROUND

Biological age is a measure of health that offers insights into ageing. The existing age clocks, although valuable, often trade off accuracy and interpretability. We introduce ExplaiNAble BioLogical Age (ENABL Age), a computational framework that combines machine-learning models with explainable artificial intelligence (XAI) methods to accurately estimate biological age with individualised explanations.

METHODS

To construct the ENABL Age clock, we first predicted an age-related outcome (eg, all-cause or cause-specific mortality), and then rescaled these predictions to estimate biological age, using UK Biobank and National Health and Nutrition Examination Survey (NHANES) datasets. We adapted existing XAI methods to decompose individual ENABL Ages into contributing risk factors. For broad accessibility, we developed two versions: ENABL Age-L, based on blood tests, and ENABL Age-Q, based on questionnaire characteristics. Finally, we validated diverse ageing mechanisms captured by each ENABL Age clock through genome-wide association studies (GWAS) association analyses.

FINDINGS

Our ENABL Age clock was significantly correlated with chronological age (r=0·7867, p<0·0001 for UK Biobank; r=0·7126, p<0·0001 for NHANES). These clocks distinguish individuals who are healthy (ie, their ENABL Age is lower than their chronological age) from those who are unhealthy (ie, their ENABL Age is higher than their chronological age), predicting mortality more effectively than existing clocks. Groups of individuals who were unhealthy showed approximately three to 12 times higher log hazard ratio than healthy groups, as per ENABL Age. The clocks achieved high mortality prediction power with an area under the receiver operating characteristic curve of 0·8179 for 5-year mortality and 0·8115 for 10-year mortality on the UK Biobank dataset, and 0·8935 for 5-year mortality and 0·9107 for 10-year mortality on the NHANES dataset. The individualised explanations that revealed the contribution of specific characteristics to ENABL Age provided insights into the important characteristics for ageing. An association analysis with risk factors and ageing-related morbidities and GWAS results on ENABL Age clocks trained on different mortality causes showed that each clock captures distinct ageing mechanisms.

INTERPRETATION

ENABL Age brings an important leap forward in the application of XAI for interpreting biological age clocks. ENABL Age also carries substantial potential in practical settings, assisting medical professionals in untangling the complexity of ageing mechanisms, and potentially becoming a valuable tool in informed clinical decision-making processes.

FUNDING

National Science Foundation and National Institutes of Health.

Causal relationship between depression and aging: a bidirectional two-sample Mendelian randomization study

BACKGROUND

The causal relationship and the direction of the effect between depression and aging remain controversial.

METHODS

We used a bidirectional two-sample Mendelian randomization analysis to examine the relationship between depression and age proxy indicators. We obtained pooled statistics from genome-wide association studies (GWAS) on depression and the age proxy indicators. We employed five MR analysis methods to address potential biases and ensure robustness of our results, with the inverse variance weighted (IVW) method being the primary outcome. We also conducted outlier exclusion using Radial MR, MRPRESSO, and MR Steiger filters. Additionally, sensitivity analyses were performed to assess heterogeneity and pleiotropy.

RESULTS

Our MR analysis revealed that depression causally leads to shortened telomere length (β = - 0.014; P = 0.038), increased frailty index (β = 0.076; P = 0.000), and accelerated GrimAge (β = 0.249; P = 0.024). Furthermore, our findings showed that the frailty index (OR = 1.679; P = 0.001) was causally associated with an increased risk of depression. Additionally, we found that appendicular lean mass (OR = 0.929; P = 0.000) and left-hand grip strength (OR = 0.836; P = 0.014) were causally associated with a reduced risk of depression. Sensitivity analyses demonstrated the robustness of our findings.

CONCLUSIONS

Our study provides evidence that depression contributes to the accelerated aging process, resulting in decreased telomere length, increased frailty index, and accelerated GrimAge. Additionally, we found that the frailty index increases the risk of depression, while appendicular lean mass and left-handed grip strength reduce the risk of depression.

Causal association of epigenetic aging and osteoporosis: a bidirectional Mendelian randomization study

BACKGROUND

The relationship between aging and osteoporosis is well established. However, the relationship between the body's physiological age, i.e. epigenetic age, and osteoporosis is not known. Our goal is to analyze the bidirectional causal relationship between epigenetic clocks and osteoporosis using a bidirectional Mendelian randomization study.

METHODS

We used SNPs closely associated with GrimAge, Hannum, PhenoAge, and HorvathAge in epigenetic age and SNPs closely associated with femoral neck bone mineral density, lumbar spine bone mineral density, and forearm bone mineral density as instrumental variables, respectively, using the inverse variance weighting method and several other MR methods to assess the bidirectional causal relationship between epigenetic age and osteoporosis.

RESULT

There was no evidence of a clear causal relationship of epigenetic age (GrimAge, Hannum, PhenoAge, and HorvathAge) on femoral neck bone mineral density, lumbar spine bone mineral density, and forearm bone mineral density. In reverse Mendelian randomization analysis showed a significant causal effect of lumbar spine bone mineral density on GrimAge: odds ratio (OR) = 0.692, 95% confidence interval (CI) = (0.538-0.890), p = 0.004. The results suggest that a decrease in lumbar spine bone mineral density promotes an acceleration of GrimAge.

CONCLUSION

There was no significant bidirectional causal relationship between epigenetic age and osteoporosis A decrease in lumbar spine bone density may lead to an acceleration of the epigenetic clock "GrimAge". Our study provides partial evidence for a bidirectional causal effect between epigenetic age and Osteoporosis.

Causal effect of gut microbiota on DNA methylation phenotypic age acceleration: a two-sample Mendelian randomization study

The causal relationship between gut microbiota and DNA methylation phenotypic age acceleration remains unclear. This study aims to examine the causal effect of gut microbiota on the acceleration of DNA methylation phenotypic age using Mendelian randomization. A total of 212 gut microbiota were included in this study, and their 16S rRNA sequencing data were obtained from the Genome-wide Association Study (GWAS) database. The GWAS data corresponding to DNA methylation phenotypic age acceleration were selected as the outcome variable. Two-sample Mendelian randomization (TSMR) was conducted using R software. During the analysis process, careful consideration was given to address potential biases arising from linkage disequilibrium and weak instrumental variables. The results from inverse-variance weighting (IVW) analysis revealed significant associations (P < 0.05) between single nucleotide polymorphisms (SNPs) corresponding to 16 gut microbiota species and DNA methylation phenotypic age acceleration. Out of the total, 12 gut microbiota species exhibited consistent and robust causal effects. Among them, 7 displayed a significant positive correlation with the outcome while 5 species showed a significant negative correlation with the outcome. This study utilized Mendelian randomization to unravel the intricate causal effects of various gut microbiota species on DNA methylation phenotypic age acceleration.

Epigenetic aging differences between Wichí and Criollos from Argentina: Insights from genomic history and ecology

Background and objectives

Epigenetic estimators based on DNA methylation levels have emerged as promising biomarkers of human aging. These estimators exhibit natural variations across human groups, but data about indigenous populations remain underrepresented in research. This study aims to investigate differences in epigenetic estimators between two distinct human populations, both residing in the Gran Chaco region of Argentina, the Native-American Wichí, and admixed Criollos who are descendants of intermarriages between Native Americans and the first European colonizers, using a population genetic approach.

Methodology

We analyzed 24 Wichí (mean age: 39.2 ± 12.9 yo) and 24 Criollos (mean age: 41.1 ± 14.0 yo) for DNA methylation levels using the Infinium MethylationEPIC (Illumina) to calculate 16 epigenetic estimators. Additionally, we examined genome-wide genetic variation using the HumanOmniExpress BeadChip (Illumina) to gain insights into the genetic history of these populations.

Results

Our results indicate that Native-American Wichí are epigenetically older compared to Criollos according to five epigenetic estimators. Analyses within the Criollos population reveal that global ancestry does not influence the differences observed, while local (chromosomal) ancestry shows positive associations between specific SNPs located in genomic regions over-represented by Native-American ancestry and measures of epigenetic age acceleration (AgeAccelHannum). Furthermore, we demonstrate that differences in population ecologies also contribute to observed epigenetic differences.

Conclusions and implications

Overall, our study suggests that while the genomic history may partially account for the observed epigenetic differences, non-genetic factors, such as lifestyle and ecological factors, play a substantial role in the variability of epigenetic estimators, thereby contributing to variations in human epigenetic aging.

Effects of iron homeostasis on epigenetic age acceleration: a two-sample Mendelian randomization study

BACKGROUND

Epigenetic clocks constructed from DNA methylation patterns have emerged as excellent predictors of aging and aging-related health outcomes. Iron, a crucial element, is meticulously regulated within organisms, a phenomenon referred as iron homeostasis. Previous researches have demonstrated the sophisticated connection between aging and iron homeostasis. However, their causal relationship remains relatively unexplored.

RESULTS

Through two-sample Mendelian randomization (MR) utilizing the random effect inverse variance weighted (IVW) method, each standard deviation (SD) increase in serum iron was associated with increased GrimAge acceleration (GrimAA, BetaIVW = 0.27, P = 8.54E-03 in 2014 datasets; BetaIVW = 0.31, P = 1.25E-02 in 2021 datasets), HannumAge acceleration (HannumAA, BetaIVW = 0.32, P = 4.50E-03 in 2014 datasets; BetaIVW = 0.32, P = 8.03E-03 in 2021 datasets) and Intrinsic epigenetic age acceleration (IEAA, BetaIVW = 0.34, P = 5.33E-04 in 2014 datasets; BetaIVW = 0.49, P = 9.94E-04 in 2021 datasets). Similar results were also observed in transferrin saturation. While transferrin manifested a negative association with epigenetic age accelerations (EAAs) sensitivity analyses. Besides, lack of solid evidence to support a causal relationship from EAAs to iron-related biomarkers.

CONCLUSIONS

The results of present investigation unveiled the causality of iron overload on acceleration of epigenetic clocks. Researches are warranted to illuminate the underlying mechanisms and formulate strategies for potential interventions.

Causality between heart failure and epigenetic age: a bidirectional Mendelian randomization study

AIMS

Heart failure (HF) is a prevalent age-related cardiovascular disease with poor prognosis in the elderly population. This study aimed to establish the causal relationship between ageing and HF by conducting a bidirectional Mendelian randomization (MR) analysis on epigenetic age (a marker of ageing) and HF.

METHODS AND RESULTS

Genome-wide association study data for epigenetic age (GrimAge, HorvathAge, HannumAge, and PhenoAge) and HF were collected and assessed for significant genetic variables. A bidirectional MR analysis was carried out using the random-effects inverse-variance weighted (IVW) method as the primary approach, while other methods (MR-Egger, weighted median, simple mode, and weighted mode) and multiple sensitivity analyses (heterogeneity analysis, leave-one-out sensitivity analysis, and horizontal pleiotropy analysis) were employed to evaluate the impact of epigenetic age on HF and vice versa. Bidirectional MR analysis of two samples revealed that the epigenetic PhenoAge clock increased the risk of HF [IVW odds ratio (OR) 1.015, 95% confidence interval (CI) 1.002-1.028, P = 0.028 and weighted median OR 1.020, 95% CI 1.001-1.038, P = 0.039]. Other results were not statistically significant.

CONCLUSIONS

The bidirectional MR analysis demonstrated a causal link between genetically predicted epigenetic age and HF in individuals of European descent. Further research into epigenetic age in other populations and additional genetic information related to HF is warranted.

Effect of women's fertility and sexual development on epigenetic clock: Mendelian randomization study

BACKGROUND AND OBJECTIVES

In observational studies, women's fertility and sexual development traits may have implications for DNA methylation patterns, and pregnancy-related risk factors can also affect maternal DNA methylation patterns. The aim of our study is to disentangle any potential causal associations between women's fertility and sexual development traits and epigenetic clocks, as well as to search for probable mediators by using the Mendelian randomization (MR) method.

METHODS

Instrumental variables for exposures, mediators, and outcomes were adopted from genome-wide association studies data of European ancestry individuals. The potential causal relationship between women's fertility and sexual development traits and four epigenetic clocks were evaluated by inverse variance weighted method and verified by other two methods. Furthermore, we employed multivariable MR (MVMR) adjusting for hypertension, hyperglycemia, BMI changes, and insomnia. Then, combining the MVMR results and previous research, we performed two-step MR to explore the mediating effects of BMI, AFS, and AFB. Multiple sensitivity analyses were further performed to verify the robustness of our findings.

RESULTS

Leveraging two-sample MR analysis, we observed statistically significant associations between earlier age at first birth (AFB) with a higher HannumAge, PhenoAge and GrimAge acceleration(β = - 0.429, 95% CI [- 0.781 to - 0.077], p = 0.017 for HannumAge; β = - 0.571, 95% CI [- 1.006 to - 0.136], p = 0.010 for PhenoAge, and β = - 1.136, 95% CI [- 1.508 to - 0.765], p = 2.03E-09 for GrimAge respectively) and age at first sexual intercourse (AFS) with a higher HannumAge and GrimAge acceleration(β = - 0.175, 95% CI [- 0.336 to - 0.014], p = 0.033 for HannumAge; β = - 0.210, 95% CI [- 0.350 to - 0.070], p = 0.003 for GrimAge, respectively). Further analyses indicated that BMI, AFB and AFS played mediator roles in the path from women's fertility and sexual development traits to epigenetic aging.

CONCLUSIONS

Our study suggested that AFS and AFB are associated with epigenetic aging. These findings may prove valuable in informing the development of prevention strategies and interventions targeted towards women's fertility and sexual development experiences and their relationship with epigenetic aging-related diseases.

Using proteomics and metabolomics to identify therapeutic targets for senescence mediated cancer: genetic complementarity method

Background

Senescence have emerged as potential factors of lung cancer risk based on findings from many studies. However, the underlying pathogenesis of lung cancer caused by senescence is not clear. In this study, we try to explain the potential pathogenesis between senescence and lung cancer through proteomics and metabonomics. And try to find new potential therapeutic targets in lung cancer patients through network mendelian randomization (MR).

Methods

The genome-wide association data of this study was mainly obtained from a meta-analysis and the Transdisciplinary Research in Cancer of the Lung Consortium (TRICL), respectively.And in this study, we mainly used genetic complementarity methods to explore the susceptibility of aging to lung cancer. Additionally, a mediation analysis was performed to explore the potential mediating role of proteomics and metabonomics, using a network MR design.

Results

GNOVA analysis revealed a shared genetic structure between HannumAge and lung cancer with a significant genetic correlation estimated at 0.141 and 0.135, respectively. MR analysis showed a relationship between HannumAge and lung cancer, regardless of smoking status. Furthermore, genetically predicted HannumAge was consistently associated with the proteins C-type lectin domain family 4 member D (CLEC4D) and Retinoic acid receptor responder protein 1 (RARR-1), indicating their potential role as mediators in the causal pathway.

Conclusion

HannumAge acceleration may increase the risk of lung cancer, some of which may be mediated by CLEC4D and RARR-1, suggestion that CLEC4D and RARR-1 may serve as potential drug targets for the treatment of lung cancer.

Multivariate genome-wide analysis of aging-related traits identifies novel loci and new drug targets for healthy aging

The concept of aging is complex, including many related phenotypes such as healthspan, lifespan, extreme longevity, frailty and epigenetic aging, suggesting shared biological underpinnings; however, aging-related endpoints have been primarily assessed individually. Using data from these traits and multivariate genome-wide association study methods, we modeled their underlying genetic factor ('mvAge'). mvAge (effective n = ~1.9 million participants of European ancestry) identified 52 independent variants in 38 genomic loci. Twenty variants were novel (not reported in input genome-wide association studies). Transcriptomic imputation identified age-relevant genes, including VEGFA and PHB1. Drug-target Mendelian randomization with metformin target genes showed a beneficial impact on mvAge (P value = 8.41 × 10-5). Similarly, genetically proxied thiazolidinediones (P value = 3.50 × 10-10), proprotein convertase subtilisin/kexin 9 inhibition (P value = 1.62 × 10-6), angiopoietin-like protein 4, beta blockers and calcium channel blockers also had beneficial Mendelian randomization estimates. Extending the drug-target Mendelian randomization framework to 3,947 protein-coding genes prioritized 122 targets. Together, these findings will inform future studies aimed at improving healthy aging.

Improving comparative analyses of Hi-C data via contrastive self-supervised learning

Hi-C is a widely applied chromosome conformation capture (3C)-based technique, which has produced a large number of genomic contact maps with high sequencing depths for a wide range of cell types, enabling comprehensive analyses of the relationships between biological functionalities (e.g. gene regulation and expression) and the three-dimensional genome structure. Comparative analyses play significant roles in Hi-C data studies, which are designed to make comparisons between Hi-C contact maps, thus evaluating the consistency of replicate Hi-C experiments (i.e. reproducibility measurement) and detecting statistically differential interacting regions with biological significance (i.e. differential chromatin interaction detection). However, due to the complex and hierarchical nature of Hi-C contact maps, it remains challenging to conduct systematic and reliable comparative analyses of Hi-C data. Here, we proposed sslHiC, a contrastive self-supervised representation learning framework, for precisely modeling the multi-level features of chromosome conformation and automatically producing informative feature embeddings for genomic loci and their interactions to facilitate comparative analyses of Hi-C contact maps. Comprehensive computational experiments on both simulated and real datasets demonstrated that our method consistently outperformed the state-of-the-art baseline methods in providing reliable measurements of reproducibility and detecting differential interactions with biological meanings.

Genetic Evidence for Causal Effects of Socioeconomic, Lifestyle, and Cardiometabolic Factors on Epigenetic-Age Acceleration

GrimAge acceleration (GrimAgeAccel) and PhenoAge acceleration (PhenoAgeAccel) are DNA methylation-based markers of accelerated biological aging, standing out in predicting mortality and age-related cardiometabolic morbidities. Causal risk factors for GrimAgeAccel and PhenoAgeAccel are unclear. In this study, we performed 2-sample univariable and multivariable Mendelian randomization (MR) to investigate causal associations of 19 modifiable socioeconomic, lifestyle, and cardiometabolic factors with GrimAgeAccel and PhenoAgeAccel. Instrument variants representing 19 modifiable factors were extracted from genome-wide association studies (GWASs) with up to 1 million Europeans. Summary statistics for GrimAgeAccel and PhenoAgeAccel were derived from a GWAS of 34 710 Europeans. We identified 12 and 8 factors causally associated with GrimAgeAccel and PhenoAgeAccel, respectively. Smoking was the strongest risk factor (β [standard error {SE}]: 1.299 [0.107] year) for GrimAgeAccel, followed by higher alcohol intake, higher waist circumference, daytime napping, higher body fat percentage, higher body mass index, higher C-reactive protein, higher triglycerides, childhood obesity, and type 2 diabetes; whereas education was the strongest protective factor (β [SE]: -1.143 [0.121] year), followed by household income. Furthermore, higher waist circumference (β [SE]: 0.850 [0.269] year) and education (β [SE]: -0.718 [0.151] year) were the leading causal risk and protective factors for PhenoAgeAccel, respectively. Sensitivity analyses strengthened the robustness of these causal associations. Multivariable MR analyses further demonstrated independent effects of the strongest risk and protective factors on GrimAgeAccel and PhenoAgeAccel, respectively. In conclusion, our findings provide novel quantitative evidence on modifiable causal risk factors for accelerated epigenetic aging, suggesting promising intervention targets against age-related morbidity and improving healthy longevity.

Multi-ancestry meta-analysis and fine-mapping in Alzheimer's disease

Genome-wide association studies (GWAS) of Alzheimer's disease are predominantly carried out in European ancestry individuals despite the known variation in genetic architecture and disease prevalence across global populations. We leveraged published GWAS summary statistics from European, East Asian, and African American populations, and an additional GWAS from a Caribbean Hispanic population using previously reported genotype data to perform the largest multi-ancestry GWAS meta-analysis of Alzheimer's disease and related dementias to date. This method allowed us to identify two independent novel disease-associated loci on chromosome 3. We also leveraged diverse haplotype structures to fine-map nine loci with a posterior probability >0.8 and globally assessed the heterogeneity of known risk factors across populations. Additionally, we compared the generalizability of multi-ancestry- and single-ancestry-derived polygenic risk scores in a three-way admixed Colombian population. Our findings highlight the importance of multi-ancestry representation in uncovering and understanding putative factors that contribute to risk of Alzheimer's disease and related dementias.

Biomarkers selection and mathematical modeling in biological age estimation

Biological age (BA) is important for clinical monitoring and preventing aging-related disorders and disabilities. Clinical and/or cellular biomarkers are measured and integrated in years using mathematical models to display an individual's BA. To date, there is not yet a single or set of biomarker(s) and technique(s) that is validated as providing the BA that reflects the best real aging status of individuals. Herein, a comprehensive overview of aging biomarkers is provided and the potential of genetic variations as proxy indicators of the aging state is highlighted. A comprehensive overview of BA estimation methods is also provided as well as a discussion of their performances, advantages, limitations, and potential approaches to overcome these limitations.

Emerging evidence on the role of clonal hematopoiesis of indeterminate potential in chronic kidney disease

Chronic kidney disease (CKD) was responsible for 1.2 million deaths globally in 2016. Despite the large and growing burden of CKD, treatment options are limited and generally only preserve kidney function. Characterizing molecular precursors to incident and progressive CKD could point to critically needed prevention and treatment strategies. Clonal hematopoiesis of indeterminate potential (CHIP) is typically characterized by the clonal expansion of blood cells carrying somatic mutations in specific driver genes. An age-related disorder, CHIP is rare in the young but common in older adults. Recent studies have identified causal associations between CHIP and atherosclerotic cardiovascular disease which are most likely mediated by inflammation, a hallmark of CKD. Animal evidence has supported causal effects of CHIP on kidney injury, inflammation, and fibrosis, providing impetus for human research. Although prospective epidemiologic studies investigating associations of CHIP with development and progression of CKD are few, intriguing findings have been reported. CHIP was significantly associated with kidney function decline and end stage kidney disease in the general population, although effect sizes were modest. Recent work suggests larger associations of CHIP with kidney disease progression in CKD patients, but further investigations in this area are needed. In addition, the accumulating literature has identified some heterogeneity in associations between CHIP and kidney endpoints across study populations, but reasons for these differences remain unclear. The current review provides an in-depth exploration into this nascent area of research, develops a conceptual framework linking CHIP to CKD, and discusses the clinical and public health implications of this work.

Phenome-wide analyses identify an association between the parent-of-origin effects dependent methylome and the rate of aging in humans

BACKGROUND

The variation in the rate at which humans age may be rooted in early events acting through the genomic regions that are influenced by such events and subsequently are related to health phenotypes in later life. The parent-of-origin-effect (POE)-regulated methylome includes regions enriched for genetically controlled imprinting effects (the typical type of POE) and regions influenced by environmental effects associated with parents (the atypical POE). This part of the methylome is heavily influenced by early events, making it a potential route connecting early exposures, the epigenome, and aging. We aim to test the association of POE-CpGs with early and later exposures and subsequently with health-related phenotypes and adult aging.

RESULTS

We perform a phenome-wide association analysis for the POE-influenced methylome using GS:SFHS (Ndiscovery = 5087, Nreplication = 4450). We identify and replicate 92 POE-CpG-phenotype associations. Most of the associations are contributed by the POE-CpGs belonging to the atypical class where the most strongly enriched associations are with aging (DNAmTL acceleration), intelligence, and parental (maternal) smoking exposure phenotypes. A proportion of the atypical POE-CpGs form co-methylation networks (modules) which are associated with these phenotypes, with one of the aging-associated modules displaying increased within-module methylation connectivity with age. The atypical POE-CpGs also display high levels of methylation heterogeneity, fast information loss with age, and a strong correlation with CpGs contained within epigenetic clocks.

CONCLUSIONS

These results identify the association between the atypical POE-influenced methylome and aging and provide new evidence for the "early development of origin" hypothesis for aging in humans.

Appraising the associations between systemic iron status and epigenetic clocks: A genetic correlation and bidirectional Mendelian Randomization study

BACKGROUND

Genetic correlations and bidirectional causal effects between systemic iron status and epigenetic clocks have not been fully investigated, although observational studies have suggested systemic iron status is associated with human aging.

OBJECTIVES

We explored the genetic correlations and bidirectional causal effects between systemic iron status and epigenetic clocks.

METHODS

Leveraging large-scale genome-wide association study summary-level statistics for four systemic iron status biomarkers (ferritin, serum iron, transferrin, transferrin saturation) (N = 48,972) and four measures for epigenetic age (GrimAge, PhenoAge, IEAA, HannumAge) (N = 34,710), genetic correlations and bidirectional causal effects were estimated between them mainly by applying linkage disequilibrium score (LDSC) regression, Mendelian randomization (MR), and MR based on Bayesian model averaging (MR-BMA). The main analyses were conducted employing multiplicative random effects inverse variance weighted MR. MR-Egger, weighted median, weighted mode, and MR-PRESSO were performed as sensitivity analyses to support the robustness of causal effects.

RESULTS

The LDSC results illustrated genetic correlations (Rg) between serum iron and PhenoAge (Rg = 0.1971, p = 0.048) and between transferrin saturation and PhenoAge (Rg = 0.196, p = 0.0469). We found that increased ferritin and transferrin saturation significantly increased all four measures of epigenetic age acceleration (all p < 0.0125, beta > 0). Each standard deviation genetically increases in serum iron only significantly associated with increased IEAA acceleration (beta = 0.36, 95% CI 0.16-0.57, p = 6.01E-04) and increased HannumAge acceleration (beta = 0.32, 95% CI 0.11-0.52, p = 2.69E-03). Evidence showed a suggestively significant causal effect of transferrin on epigenetic age acceleration (all 0.0125 < p <0.05). Additionally, reverse MR study indicated no significant causal effect of epigenetic clocks on systemic iron status.

CONCLUSIONS

All four iron status biomarkers had a significant or suggestively significant causal effect on epigenetic clocks, whereas reverse MR studies did not.

Assessing the Causal Association between Biological Aging Biomarkers and the Development of Cerebral Small Vessel Disease: A Mendelian Randomization Study

Biological aging biomarkers, such as leukocyte telomere length (LTL) and epigenetic clocks, have been associated with the risk of cerebral small vessel disease (CSVD) in several observational studies. However, it is unclear whether LTL or epigenetic clocks play causal roles as prognostic biomarkers in the development of CSVD. We performed a Mendelian randomization (MR) study of LTL and four epigenetic clocks on ten subclinical and clinical CSVD measures. We obtained genome-wide association (GWAS) data for LTL from the UK Biobank (N = 472,174). Data on epigenetic clocks were derived from a meta-analysis (N = 34,710), and CSVD data (N cases =1293-18,381; N controls = 25,806-105,974) were extracted from the Cerebrovascular Disease Knowledge Portal. We found that genetically determined LTL and epigenetic clocks were not individually associated with ten measures of CSVD (IVW p > 0.05), and this result was consistent across sensitivity analyses. Our findings imply that LTL and epigenetic clocks may not help in predicting CSVD development as causal prognostic biomarkers. Further studies are needed to illustrate the potential of reverse biological aging in serving as an effective form of preventive therapy for CSVD.

Multi-omic underpinnings of epigenetic aging and human longevity

Biological aging is accompanied by increasing morbidity, mortality, and healthcare costs; however, its molecular mechanisms are poorly understood. Here, we use multi-omic methods to integrate genomic, transcriptomic, and metabolomic data and identify biological associations with four measures of epigenetic age acceleration and a human longevity phenotype comprising healthspan, lifespan, and exceptional longevity (multivariate longevity). Using transcriptomic imputation, fine-mapping, and conditional analysis, we identify 22 high confidence associations with epigenetic age acceleration and seven with multivariate longevity. FLOT1, KPNA4, and TMX2 are novel, high confidence genes associated with epigenetic age acceleration. In parallel, cis-instrument Mendelian randomization of the druggable genome associates TPMT and NHLRC1 with epigenetic aging, supporting transcriptomic imputation findings. Metabolomics Mendelian randomization identifies a negative effect of non-high-density lipoprotein cholesterol and associated lipoproteins on multivariate longevity, but not epigenetic age acceleration. Finally, cell-type enrichment analysis implicates immune cells and precursors in epigenetic age acceleration and, more modestly, multivariate longevity. Follow-up Mendelian randomization of immune cell traits suggests lymphocyte subpopulations and lymphocytic surface molecules affect multivariate longevity and epigenetic age acceleration. Our results highlight druggable targets and biological pathways involved in aging and facilitate multi-omic comparisons of epigenetic clocks and human longevity.

Effects of epigenetic age acceleration on kidney function: a Mendelian randomization study

BACKGROUND

Previous studies have reported cross-sectional associations between measures of epigenetic age acceleration (EAA) and kidney function phenotypes. However, the temporal and potentially causal relationships between these variables remain unclear. We conducted a bidirectional two-sample Mendelian randomization study of EAA and kidney function. Genetic instruments for EAA and estimate glomerular filtration rate (eGFR) were identified from previous genome-wide association study (GWAS) meta-analyses of European-ancestry participants. Causal effects of EAA on kidney function and kidney function on EAA were assessed through summary-based Mendelian randomization utilizing data from the CKDGen GWAS meta-analysis of log-transformed estimated glomerular filtration rate (log-eGFR; n = 5,67,460) and GWAS meta-analyses of EAA (n = 34,710). An allele score-based Mendelian randomization leveraging individual-level data from UK Biobank participants (n = 4,33,462) further examined the effects of EAA on kidney function.

RESULTS

Using summary-based Mendelian randomization, we found that each 5 year increase in intrinsic EAA (IEAA) and GrimAge acceleration (GrimAA) was associated with - 0.01 and - 0.02 unit decreases in log-eGFR, respectively (P = 0.02 and P = 0.09, respectively), findings which were strongly supported by allele-based Mendelian randomization study (both P < 0.001). Summary-based Mendelian randomization identified 24% increased odds of CKD with each 5-unit increase in IEAA (P = 0.05), with consistent findings observed in allele score-based analysis (P = 0.07). Reverse-direction Mendelian randomization identified potentially causal effects of decreased kidney function on HannumAge acceleration (HannumAA), GrimAA, and PhenoAge acceleration (PhenoAA), conferring 3.14, 1.99, and 2.88 year decreases in HanumAA, GrimAA, and PhenoAA, respectively (P = 0.003, 0.05, and 0.002, respectively) with each 1-unit increase in log-eGFR.

CONCLUSION

This study supports bidirectional causal relationships between EAA and kidney function, pointing to potential prevention and therapeutic strategies.

Circulating adipokine concentrations and the risk of venous thromboembolism: A Mendelian randomization and mediation analysis

Background: Previous observational studies have suggested that circulating adipokine concentrations are related to a greater risk of venous thromboembolism (VTE). However, it remained unclear whether these observations reflect causality.

Objective: This study aimed to investigate the causal relationship between circulating adipokine concentrations (including adiponectin, leptin, PAI-1, MCP-1, leptin receptor, and RETN) and the risk of VTE and its subtypes (DVT and PE) and to determine whether circulating adipokine concentrations are a mediator of venous thromboembolic events in obese patients.

Methods: We used Mendelian randomization (MR) analyses to determine the effects of the body mass index (BMI), adiponectin, leptin, PAI-1, MCP-1, leptin receptor, and RETN levels on VTE, DVT, and PE in a cohort of 11,288 VTE cases, 5,632 DVT cases, 5,130 PE cases, and 254,771 controls. We then assessed the proportion of the effect of obesity on VTE, DVT, and PE explained by circulating leptin levels.

Result: Genetically predicted higher BMI was related to increased VTE (OR = 1.45, p &lt; 0.001), DVT (OR = 1.63, p &lt; 0.001), and PE (OR = 1.37, p &lt; 0.001) risk, and higher circulating leptin levels increase odds of VTE (OR = 1.96, q &lt; 0.001), DVT (OR = 2.52, q &lt; 0.001), and PE (OR = 2.26, q = 0.005). In addition, we found that the causal effect between elevated serum adiponectin and the decreased risk of VTE (OR = 0.85, p = 0.013, q = 0.053) and PE (OR = 0.81, p = 0.032, q = 0.083) and between MCP-1 and the reduced risk of VTE (OR = 0.88, p = 0.048, q = 0.143) is no longer significant after FDR adjustment. In MR mediation analysis, the mediation effect of circulating leptin levels in the causal pathway from BMI to PE was estimated to be 1.28 (0.95–1.71, p = 0.10), accounting for 39.14% of the total effect.

Conclusion: The circulating leptin level is a risk factor for VTE, DVT, and PE, but it might be a potential mediator of BMI on the risk of PE, and thus, interventions on the circulating leptin level in obesity might reduce the risk of PE. Adiponectin is a potential protective factor for both VTE and PE.

Identifying COPD subtypes using multi-trait genetics

Chronic Obstructive Pulmonary Disease (COPD) has a simple physiological diagnostic criterion but a wide range of clinical characteristics. The mechanisms underlying this variability in COPD phenotypes are unclear. To investigate the potential contribution of genetic variants to phenotypic heterogeneity, we examined the association of genome-wide associated lung function, COPD, and asthma variants with other phenotypes using phenome-wide association results derived in the UK Biobank. Our clustering analysis of the variants-phenotypes association matrix identified three clusters of genetic variants with different effects on white blood cell counts, height, and body mass index (BMI). To assess the potential clinical and molecular effects of these groups of variants, we investigated the association between cluster-specific genetic risk scores and phenotypes in the COPDGene cohort. We observed differences in steroid use, BMI, lymphocyte counts, chronic bronchitis, and differential gene and protein expression across the three genetic risk scores. Our results suggest that multi-phenotype analysis of obstructive lung disease-related risk variants may identify genetically driven phenotypic patterns in COPD.

Investigation of genetic variants and causal biomarkers associated with brain aging

Delta age is a biomarker of brain aging that captures differences between the chronological age and the predicted biological brain age. Using multimodal data of brain MRI, genomics, and blood-based biomarkers and metabolomics in UK Biobank, this study investigates an explainable and causal basis of high delta age. A visual saliency map of brain regions showed that lower volumes in the fornix and the lower part of the thalamus are key predictors of high delta age. Genome-wide association analysis of the delta age using the SNP array data identified associated variants in gene regions such as KLF3-AS1 and STX1. GWAS was also performed on the volumes in the fornix and the lower part of the thalamus, showing a high genetic correlation with delta age, indicating that they share a genetic basis. Mendelian randomization (MR) for all metabolomic biomarkers and blood-related phenotypes showed that immune-related phenotypes have a causal impact on increasing delta age. Our analysis revealed regions in the brain that are susceptible to the aging process and provided evidence of the causal and genetic connections between immune responses and brain aging.

Phenome-wide analysis identifies parent-of-origin effects on the human methylome associated with changes in the rate of aging

Variation in the rate at which humans age may be rooted in early life events acting through genomic regions that are influenced by such events and subsequently are related to health phenotypes in later life. The parent-of-origin-effect (POE)-regulated methylome includes regions either enriched for genetically controlled imprinting effects (the typical type of POE) or atypical POE introduced by environmental effects associated with parents. This part of the methylome is heavily influenced by early life events, making it a potential route connecting early environmental exposures, the epigenome and the rate of aging. Here, we aim to test the association of POE-influenced methylation of CpG dinucleotides (POE-CpG sites) with early and later environmental exposures and subsequently with health-related phenotypes and adult aging phenotypes. We do this by performing phenome-wide association analyses of the POE-influenced methylome using a large family-based population cohort (GS:SFHS, Ndiscovery=5,087, Nreplication=4,450). At the single CpG level, 92 associations of POE-CpGs with phenotypic variation were identified and replicated. Most of the associations were contributed by POE-CpGs belonging to the atypical class and the most strongly enriched associations were with aging (DNAmTL acceleration), intelligence and parental (maternal) smoking exposure phenotypes. We further found that a proportion of the atypical-POE-CpGs formed co-methylation networks (modules) which are associated with these phenotypes, with one of the aging-associated modules displaying increased internal module connectivity (strength of methylation correlation across constituent CpGs) with age. Atypical POE-CpGs also displayed high levels of methylation heterogeneity and epigenetic drift (i.e. information loss with age) and a strong correlation with CpGs contained within epigenetic clocks. These results identified associations between the atypical-POE-influenced methylome and aging and provided new evidence for the "early development of origin" hypothesis for aging in humans.

Association between Dried Fruit Intake and DNA Methylation: A Multivariable Mendelian Randomization Analysis

OBJECTIVES

Observational studies have reported associations between dried fruit intake and DNA methylation(DNAm). However, inherent flaws in observational study designs make them susceptible to confounding and reverse causality bias. Consequently, it is unclear whether a causal association exists. In the present study, we aimed to investigate the causal associations between dried fruit intake and DNAm.

METHODS

We performed two-sample Mendelian randomization (MR) using the IEU Open GWAS database aggregated data. Forty-three single nucleotide polymorphisms (SNPs) associated with dried fruit intake as instrumental variables (IVs) were selected as exposure. DNAm outcomes include Gran (estimated granulocyte proportions); AgeAccelGrim(GrimAge acceleration); Hannum (Hannum age acceleration); IEAA(Intrinsic epigenetic age acceleration), AgeAccelPheno( PhenoAge acceleration), and DNAmPAIadjAge (DNAm-estimated plasminogen activator inhibitor-1 levels). We used the MR pleiotropy residual sum and outlier test (MRPRESSO) and Radial-MR test to identify any level of multi-effect outliers and assessed the causal effect estimates(after removing outliers). The primary causal effects were estimated using inverse-variance weighted (IVW) method and undertook sensitivity analyses using MR methods robust to horizontal pleiotropy.The direct effects of dried fruit intake on DNAm were estimated using multivariable mendelian randomization (MVMR).

RESULTS

Leveraging two-sample MR analysis, we observed statistically significant associations between dried fruit intake with a lower AgeAccelGrim(β=-1.365, 95% confidence intervals [CI] -2.266 to -0.464, PIVW=2.985×10-3) and AgeAccelPheno (β= -1.933, 95% CI -3.068 to -0.798, PIVW=8.371×10-4). By contrast, the effects level on Gran (β=0.008, PIVW=0.430), Hannum(β=-0.430, PIVW=0.357), IEAA(β=-0.184, PIVW=0.700), and DNAmPAIadjAge (β=-1.861, PIVW=0.093) were not statistically significant. MVMR results adjusting for the potential effects of confounders showed that the causal relationship between dried fruit intake and AgeAccelGrim(β= -1.315, 95% CI -2.373 to -0.258, PIVW=1.480×10-2) and AgeAccelPheno(β= -1.595, 95% CI -2.987 to -0.202, PIVW=2.483×10-2) persisted. No significant horizontal polymorphism was found in the sensitivity analysis.

CONCLUSION

Our MR study suggested that increased dried fruit intake is associated with slower AgeAccelGrim and AgeAccelPheno. It can providing a promising avenue for exploring the beneficial effects of dried fruit intake on lifespan extension.

Genetic Variants in Telomerase Reverse Transcriptase Contribute to Solar Lentigines

Solar lentigines (SLs) are a hallmark of human skin aging. They result from chronic exposure to sunlight and other environmental stressors. Recent studies also imply genetic factors, but findings are partially conflicting and lack of replication. Through a multi-trait based analysis strategy, we discovered that genetic variants in telomerase reverse transcriptase were significantly associated with non-facial SL in two East Asian (Taizhou longitudinal cohort, n = 2,964 and National Survey of Physical Traits, n = 2,954) and one Caucasian population (SALIA, n = 462), top SNP rs2853672 (P-value for Taizhou longitudinal cohort = 1.32 × 10^‒28 and P-value for National Survey of Physical Traits = 3.66 × 10^‒17 and P-value for SALIA = 0.0007 and P_meta = 4.93 × 10^‒44). The same variants were nominally associated with facial SL but not with other skin aging or skin pigmentation traits. The SL-enhanced allele/haplotype upregulated the transcription of the telomerase reverse transcriptase gene. Of note, well-known telomerase reverse transcriptase‒related aging markers such as leukocyte telomere length and intrinsic epigenetic age acceleration were not associated with SL. Our results indicate a previously unrecognized role of telomerase reverse transcriptase in skin aging‒related lentigines formation.

DNA methylation GrimAge version 2

We previously described a DNA methylation (DNAm) based biomarker of human mortality risk DNAm GrimAge. Here we describe version 2 of GrimAge (trained on individuals aged between 40 and 92) which leverages two new DNAm based estimators of (log transformed) plasma proteins: high sensitivity C-reactive protein (logCRP) and hemoglobin A1C (logA1C). We evaluate GrimAge2 in 13,399 blood samples across nine study cohorts. After adjustment for age and sex, GrimAge2 outperforms GrimAge in predicting mortality across multiple racial/ethnic groups (meta P=3.6x10^-167 versus P=2.6x10^-144) and in terms of associations with age related conditions such as coronary heart disease, lung function measurement FEV1 (correlation= -0.31, P=1.1x10^-136), computed tomography based measurements of fatty liver disease. We present evidence that GrimAge version 2 also applies to younger individuals and to saliva samples where it tracks markers of metabolic syndrome. DNAm logCRP is positively correlated with morbidity count (P=1.3x10^-54). DNAm logA1C is highly associated with type 2 diabetes (P=5.8x10^-155). DNAm PAI-1 outperforms the other age-adjusted DNAm biomarkers including GrimAge2 in correlating with triglyceride (cor=0.34, P=9.6x10^-267) and visceral fat (cor=0.41, P=4.7x10^-41). Overall, we demonstrate that GrimAge version 2 is an attractive epigenetic biomarker of human mortality and morbidity risk.

The association of accelerated epigenetic age with all-cause mortality in cardiac catheterization patients as mediated by vascular and cardiometabolic outcomes

BACKGROUND

Epigenetic age is a DNA methylation-based biomarker of aging that is accurate across the lifespan and a range of cell types. The difference between epigenetic age and chronological age, termed age acceleration (AA), is a strong predictor of lifespan and healthspan. The predictive capabilities of AA for all-cause mortality have been evaluated in the general population; however, its utility is less well evaluated in those with chronic conditions. Additionally, the pathophysiologic pathways whereby AA predicts mortality are unclear. We hypothesized that AA predicts mortality in individuals with underlying cardiovascular disease; and the association between AA and mortality is mediated, in part, by vascular and cardiometabolic measures.

METHODS

We evaluated 562 participants in an urban, three-county area of central North Carolina from the CATHGEN cohort, all of whom received a cardiac catheterization procedure. We analyzed three AA biomarkers, Horvath epigenetic age acceleration (HAA), phenotypic age acceleration (PhenoAA), and Grim age acceleration (GrimAA), by Cox regression models, to assess whether AAs were associated with all-cause mortality. We also evaluated if these associations were mediated by vascular and cardiometabolic outcomes, including left ventricular ejection fraction (LVEF), blood cholesterol concentrations, angiopoietin-2 (ANG2) protein concentration, peripheral artery disease, coronary artery disease, diabetes, and hypertension. The total effect, direct effect, indirect effect, and percentage mediated were estimated using pathway mediation tests with a regression adjustment approach.

RESULTS

PhenoAA (HR = 1.05, P < 0.0001), GrimAA (HR = 1.10, P < 0.0001) and HAA (HR = 1.03, P = 0.01) were all associated with all-cause mortality. The association of mortality and PhenoAA was partially mediated by ANG2, a marker of vascular function (19.8%, P = 0.016), and by diabetes (8.2%, P = 0.043). The GrimAA-mortality association was mediated by ANG2 (12.3%, P = 0.014), and showed weaker evidence for mediation by LVEF (5.3%, P = 0.065).

CONCLUSIONS

Epigenetic age acceleration remains strongly predictive of mortality even in individuals already burdened with cardiovascular disease. Mortality associations were mediated by ANG2, which regulates endothelial permeability and angiogenic functions, suggesting that specific vascular pathophysiology may link accelerated epigenetic aging with increased mortality risks.

Measuring biological age using omics data

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