Exome Sequencing Analysis Identifies Rare Variants in ATM and RPL8 That Are Associated With Shorter Telomere Length

Whole exome sequencing analyses reveal novel genes in telomere length and their biomedical implications

Telomere length is a putative biomarker of aging and is associated with multiple age-related diseases. There are limited data on the landscape of rare genetic variations in telomere length. Here, we systematically characterize the rare variant associations with leukocyte telomere length (LTL) through exome-wide association study (ExWAS) among 390,231 individuals in the UK Biobank. We identified 18 robust rare-variant genes for LTL, most of which estimated effects on LTL were significant (> 0.2 standard deviation per allele). The biological functions of the rare-variant genes were associated with telomere maintenance and capping and several genes were specifically expressed in the testis. Three novel genes (ASXL1, CFAP58, and TET2) associated with LTL were identified. Phenotypic association analyses indicated significant associations of ASXL1 and TET2 with cancers, age-related diseases, blood assays, and cardiovascular traits. Survival analyses suggested that carriers of ASXL1 or TET2 variants were at increased risk for cancers; diseases of the circulatory, respiratory, and genitourinary systems; and all-cause and cause-specific deaths. The CFAP58 carriers were at elevated risk of deaths due to cancers. Collectively, the present whole exome sequencing study provides novel insights into the genetic landscape of LTL, identifying novel genes associated with LTL and their implications on human health and facilitating a better understanding of aging, thus pinpointing the genetic relevance of LTL with clonal hematopoiesis, biomedical traits, and health-related outcomes.

Identification of prospective aging drug targets via Mendelian randomization analysis

Aging represents a multifaceted process culminating in the deterioration of biological functions. Despite the introduction of numerous anti-aging strategies, their therapeutic outcomes have often been less than optimal. Consequently, discovering new targets to mitigate aging effects is of critical importance. We applied Mendelian randomization (MR) to identify potential pharmacological targets against aging, drawing upon summary statistics from both the Decode and FinnGen cohorts, with further validation in an additional cohort. To address potential reverse causality, bidirectional MR analysis with Steiger filtering was utilized. Additionally, Bayesian co-localization and phenotype scanning were implemented to investigate previous associations between genetic variants and traits. Summary-data-based Mendelian randomization (SMR) analysis was conducted to assess the impact of genetic variants on aging via their effects on protein expression. Additionally, mediation analysis was orchestrated to uncover potential intermediaries in these associations. Finally, we probed the systemic implications of drug-target protein expression across diverse indications by MR-PheWas analysis. Utilizing a Bonferroni-corrected threshold, our MR examination identified 10 protein-aging associations. Within this cohort of proteins, MST1, LCT, GMPR2, PSMB4, ECM1, EFEMP1, and ISLR2 appear to exacerbate aging risks, while MAX, B3GNT8, and USP8 may exert protective influences. None of these proteins displayed reverse causality except EFEMP1. Bayesian co-localization inferred shared variants between aging and proteins such as B3GNT8 (rs11670143), ECM1 (rs61819393), and others listed. Mediator analysis pinpointed 1,5-anhydroglucitol as a partial intermediary in the influence LCT exhibits on telomere length. Circulating proteins play a pivotal role in influencing the aging process, making them promising candidates for therapeutic intervention. The implications of these proteins in aging warrant further investigation in future clinical research.

Differential Methylation of Telomere-Related Genes Is Associated with Kidney Disease in Individuals with Type 1 Diabetes

Diabetic kidney disease (DKD) represents a major global health problem. Accelerated ageing is a key feature of DKD and, therefore, characteristics of accelerated ageing may provide useful biomarkers or therapeutic targets. Harnessing multi-omics, features affecting telomere biology and any associated methylome dysregulation in DKD were explored. Genotype data for nuclear genome polymorphisms in telomere-related genes were extracted from genome-wide case-control association data (n = 823 DKD/903 controls; n = 247 end-stage kidney disease (ESKD)/1479 controls). Telomere length was established using quantitative polymerase chain reaction. Quantitative methylation values for 1091 CpG sites in telomere-related genes were extracted from epigenome-wide case-control association data (n = 150 DKD/100 controls). Telomere length was significantly shorter in older age groups (p = 7.6 × 10^-6). Telomere length was also significantly reduced (p = 6.6 × 10^-5) in DKD versus control individuals, with significance remaining after covariate adjustment (p = 0.028). DKD and ESKD were nominally associated with telomere-related genetic variation, with Mendelian randomisation highlighting no significant association between genetically predicted telomere length and kidney disease. A total of 496 CpG sites in 212 genes reached epigenome-wide significance (p ≤ 10^-8) for DKD association, and 412 CpG sites in 193 genes for ESKD. Functional prediction revealed differentially methylated genes were enriched for Wnt signalling involvement. Harnessing previously published RNA-sequencing datasets, potential targets where epigenetic dysregulation may result in altered gene expression were revealed, useful as potential diagnostic and therapeutic targets for intervention.

Harnessing Genomic Analysis to Explore the Role of Telomeres in the Pathogenesis and Progression of Diabetic Kidney Disease

The prevalence of diabetes is increasing globally, and this trend is predicted to continue for future decades. Research is needed to uncover new ways to manage diabetes and its co-morbidities. A significant secondary complication of diabetes is kidney disease, which can ultimately result in the need for renal replacement therapy, via dialysis or transplantation. Diabetic kidney disease presents a substantial burden to patients, their families and global healthcare services. This review highlights studies that have harnessed genomic, epigenomic and functional prediction tools to uncover novel genes and pathways associated with DKD that are useful for the identification of therapeutic targets or novel biomarkers for risk stratification. Telomere length regulation is a specific pathway gaining attention recently because of its association with DKD. Researchers are employing both observational and genetics-based studies to identify telomere-related genes associated with kidney function decline in diabetes. Studies have also uncovered novel functions for telomere-related genes beyond the immediate regulation of telomere length, such as transcriptional regulation and inflammation. This review summarises studies that have revealed the potential to harness therapeutics that modulate telomere length, or the associated epigenetic modifications, for the treatment of DKD, to potentially slow renal function decline and reduce the global burden of this disease.

Fat metabolism is associated with telomere length in six population-based studies

Telomeres are repetitive DNA sequences located at the end of chromosomes, which are associated to biological aging, cardiovascular disease, cancer, and mortality. Lipid and fatty acid metabolism have been associated with telomere shortening. We have conducted an in-depth study investigating the association of metabolic biomarkers with telomere length (LTL). We performed an association analysis of 226 metabolic biomarkers with LTL using data from 11 775 individuals from six independent population-based cohorts (BBMRI-NL consortium). Metabolic biomarkers include lipoprotein lipids and subclasses, fatty acids, amino acids, glycolysis measures and ketone bodies. LTL was measured by quantitative polymerase chain reaction or FlowFISH. Linear regression analysis was performed adjusting for age, sex, lipid-lowering medication and cohort-specific covariates (model 1) and additionally for body mass index (BMI) and smoking (model 2), followed by inverse variance-weighted meta-analyses (significance threshold pmeta = 6.5x10-4). We identified four metabolic biomarkers positively associated with LTL, including two cholesterol to lipid ratios in small VLDL (S-VLDL-C % and S-VLDL-ce %) and two omega-6 fatty acid ratios (FAw6/FA and LA/FA). After additionally adjusting for BMI and smoking, these metabolic biomarkers remained associated with LTL with similar effect estimates. In addition, cholesterol esters in very small VLDL (XS-VLDL-ce) became significantly associated with LTL (p = 3.6x10-4). We replicated the association of FAw6/FA with LTL in an independent dataset of 7845 individuals (p = 1.9x10-4). To conclude, we identified multiple metabolic biomarkers involved in lipid and fatty acid metabolism that may be involved in LTL biology. Longitudinal studies are needed to exclude reversed causation.

Low frequency variants associated with leukocyte telomere length in the Singapore Chinese population

The role of low frequency variants associated with telomere length homeostasis in chronic diseases and mortalities is relatively understudied in the East-Asian population. Here we evaluated low frequency variants, including 1,915,154 Asian specific variants, for leukocyte telomere length (LTL) associations among 25,533 Singapore Chinese samples. Three East Asian specific variants in/near POT1, TERF1 and STN1 genes are associated with LTL (Meta-analysis P 2.49×10-14-6.94×10-10). Rs79314063, a missense variant (p.Asp410His) at POT1, shows effect 5.3 fold higher and independent of a previous common index SNP. TERF1 (rs79617270) and STN1 (rs139620151) are linked to LTL-associated common index SNPs at these loci. Rs79617270 is associated with cancer mortality [HR95%CI = 1.544 (1.173, 2.032), PAdj = 0.018] and 4.76% of the association between the rs79617270 and colon cancer is mediated through LTL. Overall, genetically determined LTL is particularly associated with lung adenocarcinoma [HR95%CI = 1.123 (1.051, 1.201), Padj = 0.007]. Ethnicity-specific low frequency variants may affect LTL homeostasis and associate with certain cancers.