The microbial composition of pancreatic ductal adenocarcinoma: A systematic review of 16S rRNA gene sequencing.

BACKGROUND: Pancreatic cancer, specifically pancreatic ductal adenocarcinoma (PDAC), continues to pose a significant clinical and scientific challenge. The most significant finding of recent years is that PDAC tumours harbour their specific microbiome, which differs amongst tumour entities and is distinct from healthy tissue. This review aims to evaluate and summarise all PDAC studies that have used the next-generation technique, 16S rRNA gene amplicon sequencing within each bodily compartment. As well as establishing a causal relationship between PDAC and the microbiome. MATERIALS AND METHODS: This systematic review was carried out according to the Preferred Reporting Items for Systematic Reviews and Meta-analysis (PRISMA) guidelines. A comprehensive search strategy was designed, and 1727 studies were analysed. RESULTS: In total, 38 studies were selected for qualitative analysis and summarised significant PDAC bacterial signatures. Despite the growing amount of data provided, we are not able to state a universal 16S rRNA gene microbial signature that can be used for PDAC screening. This is most certainly due to the heterogeneity of the presentation of results, lack of available datasets and the intrinsic selection bias between studies. CONCLUSION: Several key studies have begun to shed light on causality and the influence the microbiome constituents and their produced metabolites could play in tumorigenesis and influencing outcomes. The challenge in this field is to shape the available microbial data into targetable signatures. Making sequenced data readily available is critical, coupled with the coordinated standardisation of data and the need for consensus guidelines in studies investigating the microbiome in PDAC.

comorbidPGS: an R package assessing shared predisposition between Phenotypes using Polygenic Scores.

Introduction Polygenic Score (PGS) is a valuable method for assessing the estimated genetic liability to a given outcome or genetic variability contributing to a quantitative trait. While PRSs are widely used for complex traits, their application in uncovering shared genetic predisposition between phenotypes, i.e. when genetic variants influence more than one phenotype, remains limited. Methods We developed an R package, comorbidPGS, which facilitates a systematic evaluation of shared genetic effects among (cor)related phenotypes using PGSs. The comorbidPGS package takes as input a set of Single Nucleotide Polymorphisms (SNPs) along with their established effects on the original phenotype (Po), referred to as Po-PGS. It generates a comprehensive summary of effect(s) of Po-PGS on target phenotype(s) (Pt) with customisable graphical features. Results We applied comorbidPGS to investigate the shared genetic predisposition between phenotypes defining elevated blood pressure (Systolic Blood Pressure, SBP; Diastolic Blood Pressure, DBP; Pulse Pressure, PP) and several cancers (Breast Cancer, BrC; Pancreatic Cancer, PanC; Kidney Cancer, KidC; Prostate Cancer, PrC; Colorectal Cancer, CrC) using the European ancestry UK Biobank individuals and GWAS meta-analyses summary statistics from independent set of European ancestry individuals. We report a significant association between elevated DBP and the genetic risk of PrC (ß (SE)=0.066 (0.017), P-value=9.64×10^(-5)), as well as between CrC PGS and both, lower SBP (ß (SE)=-0.10 [0.029], P-value=3.83×10^(-4))) and lower DBP (ß (SE)=-0.055 [0.017], P-value=1.05×10^(-3)). Our analysis highlights two nominally significant relationships for individuals with genetic predisposition to elevated SBP leading to higher risk of KidC (OR [95%CI]=1.04 [1.0039-1.087], P-value=2.82×10^(-2)) and PrC (OR [95%CI]=1.02 [1.003-1.041], P-value=2.22×10^(-2)). Conclusion Using comorbidPGS, we underscore mechanistic relationships between blood pressure regulation and susceptibility to three comorbid malignancies. This package offers valuable means to evaluate shared genetic susceptibility between (cor)related phenotypes through polygenic scores.

Genome-wide characterization of circulating metabolic biomarkers.

Genome-wide association analyses using high-throughput metabolomics platforms have led to novel insights into the biology of human metabolism1-7. This detailed knowledge of the genetic determinants of systemic metabolism has been pivotal for uncovering how genetic pathways influence biological mechanisms and complex diseases8-11. Here we present a genome-wide association study for 233 circulating metabolic traits quantified by nuclear magnetic resonance spectroscopy in up to 136,016 participants from 33 cohorts. We identify more than 400 independent loci and assign probable causal genes at two-thirds of these using manual curation of plausible biological candidates. We highlight the importance of sample and participant characteristics that can have significant effects on genetic associations. We use detailed metabolic profiling of lipoprotein- and lipid-associated variants to better characterize how known lipid loci and novel loci affect lipoprotein metabolism at a granular level. We demonstrate the translational utility of comprehensively phenotyped molecular data, characterizing the metabolic associations of intrahepatic cholestasis of pregnancy. Finally, we observe substantial genetic pleiotropy for multiple metabolic pathways and illustrate the importance of careful instrument selection in Mendelian randomization analysis, revealing a putative causal relationship between acetone and hypertension. Our publicly available results provide a foundational resource for the community to examine the role of metabolism across diverse diseases.

Bile Microbiome Signatures Associated with Pancreatic Ductal Adenocarcinoma Compared to Benign Disease: A UK Pilot Study.

Pancreatic ductal adenocarcinoma (PDAC) has a very poor survival. The intra-tumoural microbiome can influence pancreatic tumourigenesis and chemoresistance and, therefore, patient survival. The role played by bile microbiota in PDAC is unknown. We aimed to define bile microbiome signatures that can effectively distinguish malignant from benign tumours in patients presenting with obstructive jaundice caused by benign and malignant pancreaticobiliary disease. Prospective bile samples were obtained from 31 patients who underwent either Endoscopic Retrograde Cholangiopancreatography (ERCP) or Percutaneous Transhepatic Cholangiogram (PTC). Variable regions (V3-V4) of the 16S rRNA genes of microorganisms present in the samples were amplified by Polymerase Chain Reaction (PCR) and sequenced. The cohort consisted of 12 PDAC, 10 choledocholithiasis, seven gallstone pancreatitis and two primary sclerosing cholangitis patients. Using the 16S rRNA method, we identified a total of 135 genera from 29 individuals (12 PDAC and 17 benign). The bile microbial beta diversity significantly differed between patients with PDAC vs. benign disease (Permanova p = 0.0173). The separation of PDAC from benign samples is clearly seen through unsupervised clustering of Aitchison distance. We found three genera to be of significantly lower abundance among PDAC samples vs. benign, adjusting for false discovery rate (FDR). These were Escherichia (FDR = 0.002) and two unclassified genera, one from Proteobacteria (FDR = 0.002) and one from Enterobacteriaceae (FDR = 0.011). In the same samples, the genus Streptococcus (FDR = 0.033) was found to be of increased abundance in the PDAC group. We show that patients with obstructive jaundice caused by PDAC have an altered microbiome composition in the bile compared to those with benign disease. These bile-based microbes could be developed into potential diagnostic and prognostic biomarkers for PDAC and warrant further investigation.

51st European Mathematical Genetics Meeting (EMGM) 2023.

NA.

Machine Learning to Advance Human Genome-Wide Association Studies.

Machine learning, including deep learning, reinforcement learning, and generative artificial intelligence are revolutionising every area of our lives when data are made available. With the help of these methods, we can decipher information from larger datasets while addressing the complex nature of biological systems in a more efficient way. Although machine learning methods have been introduced to human genetic epidemiological research as early as 2004, those were never used to their full capacity. In this review, we outline some of the main applications of machine learning to assigning human genetic loci to health outcomes. We summarise widely used methods and discuss their advantages and challenges. We also identify several tools, such as Combi, GenNet, and GMSTool, specifically designed to integrate these methods for hypothesis-free analysis of genetic variation data. We elaborate on the additional value and limitations of these tools from a geneticist's perspective. Finally, we discuss the fast-moving field of foundation models and large multi-modal omics biobank initiatives.

Metagenomics of Parkinson's disease implicates the gut microbiome in multiple disease mechanisms.

Parkinson's disease (PD) may start in the gut and spread to the brain. To investigate the role of gut microbiome, we conducted a large-scale study, at high taxonomic resolution, using uniform standardized methods from start to end. We enrolled 490 PD and 234 control individuals, conducted deep shotgun sequencing of fecal DNA, followed by metagenome-wide association studies requiring significance by two methods (ANCOM-BC and MaAsLin2) to declare disease association, network analysis to identify polymicrobial clusters, and functional profiling. Here we show that over 30% of species, genes and pathways tested have altered abundances in PD, depicting a widespread dysbiosis. PD-associated species form polymicrobial clusters that grow or shrink together, and some compete. PD microbiome is disease permissive, evidenced by overabundance of pathogens and immunogenic components, dysregulated neuroactive signaling, preponderance of molecules that induce alpha-synuclein pathology, and over-production of toxicants; with the reduction in anti-inflammatory and neuroprotective factors limiting the capacity to recover. We validate, in human PD, findings that were observed in experimental models; reconcile and resolve human PD microbiome literature; and provide a broad foundation with a wealth of concrete testable hypotheses to discern the role of the gut microbiome in PD.

Multi-ancestry genome-wide association study of gestational diabetes mellitus highlights genetic links with type 2 diabetes.

Gestational diabetes mellitus (GDM) is associated with increased risk of pregnancy complications and adverse perinatal outcomes. GDM often reoccurs and is associated with increased risk of subsequent diagnosis of type 2 diabetes (T2D). To improve our understanding of the aetiological factors and molecular processes driving the occurrence of GDM, including the extent to which these overlap with T2D pathophysiology, the GENetics of Diabetes In Pregnancy Consortium assembled genome-wide association studies of diverse ancestry in a total of 5485 women with GDM and 347 856 without GDM. Through multi-ancestry meta-analysis, we identified five loci with genome-wide significant association (P < 5 × 10-8) with GDM, mapping to/near MTNR1B (P = 4.3 × 10-54), TCF7L2 (P = 4.0 × 10-16), CDKAL1 (P = 1.6 × 10-14), CDKN2A-CDKN2B (P = 4.1 × 10-9) and HKDC1 (P = 2.9 × 10-8). Multiple lines of evidence pointed to the shared pathophysiology of GDM and T2D: (i) four of the five GDM loci (not HKDC1) have been previously reported at genome-wide significance for T2D; (ii) significant enrichment for associations with GDM at previously reported T2D loci; (iii) strong genetic correlation between GDM and T2D and (iv) enrichment of GDM associations mapping to genomic annotations in diabetes-relevant tissues and transcription factor binding sites. Mendelian randomization analyses demonstrated significant causal association (5% false discovery rate) of higher body mass index on increased GDM risk. Our results provide support for the hypothesis that GDM and T2D are part of the same underlying pathology but that, as exemplified by the HKDC1 locus, there are genetic determinants of GDM that are specific to glucose regulation in pregnancy.

Bifidobacterium Is Enriched in Gut Microbiome of Kashmiri Women with Polycystic Ovary Syndrome.

Polycystic ovary syndrome (PCOS) is a very common endocrine condition in women in India. Gut microbiome alterations were shown to be involved in PCOS, yet it is remarkably understudied in Indian women who have a higher incidence of PCOS as compared to other ethnic populations. During the regional PCOS screening program among young women, we recruited 19 drug naive women with PCOS and 20 control women at the Sher-i-Kashmir Institute of Medical Sciences, Kashmir, North India. We profiled the gut microbiome in faecal samples by 16S rRNA sequencing and included 40/58 operational taxonomic units (OTUs) detected in at least 1/3 of the subjects with relative abundance (RA) ≥ 0.1%. We compared the RAs at a family/genus level in PCOS/non-PCOS groups and their correlation with 33 metabolic and hormonal factors, and corrected for multiple testing, while taking the variation in day of menstrual cycle at sample collection, age and BMI into account. Five genera were significantly enriched in PCOS cases: Sarcina, Megasphaera, and previously reported for PCOS Bifidobacterium, Collinsella and Paraprevotella confirmed by different statistical models. At the family level, the relative abundance of Bifidobacteriaceae was enriched, whereas Peptococcaceae was decreased among cases. We observed increased relative abundance of Collinsella and Paraprevotella with higher fasting blood glucose levels, and Paraprevotella and Alkalibacterium with larger hip, waist circumference, weight, and Peptococcaceae with lower prolactin levels. We also detected a novel association between Eubacterium and follicle-stimulating hormone levels and between Bifidobacterium and alkaline phosphatase, independently of the BMI of the participants. Our report supports that there is a relationship between gut microbiome composition and PCOS with links to specific reproductive health metabolic and hormonal predictors in Indian women.

A multi-omics study of circulating phospholipid markers of blood pressure.

High-throughput techniques allow us to measure a wide-range of phospholipids which can provide insight into the mechanisms of hypertension. We aimed to conduct an in-depth multi-omics study of various phospholipids with systolic blood pressure (SBP) and diastolic blood pressure (DBP). The associations of blood pressure and 151 plasma phospholipids measured by electrospray ionization tandem mass spectrometry were performed by linear regression in five European cohorts (n = 2786 in discovery and n = 1185 in replication). We further explored the blood pressure-related phospholipids in Erasmus Rucphen Family (ERF) study by associating them with multiple cardiometabolic traits (linear regression) and predicting incident hypertension (Cox regression). Mendelian Randomization (MR) and phenome-wide association study (Phewas) were also explored to further investigate these association results. We identified six phosphatidylethanolamines (PE 38:3, PE 38:4, PE 38:6, PE 40:4, PE 40:5 and PE 40:6) and two phosphatidylcholines (PC 32:1 and PC 40:5) which together predicted incident hypertension with an area under the ROC curve (AUC) of 0.61. The identified eight phospholipids are strongly associated with triglycerides, obesity related traits (e.g. waist, waist-hip ratio, total fat percentage, body mass index, lipid-lowering medication, and leptin), diabetes related traits (e.g. glucose, insulin resistance and insulin) and prevalent type 2 diabetes. The genetic determinants of these phospholipids also associated with many lipoproteins, heart rate, pulse rate and blood cell counts. No significant association was identified by bi-directional MR approach. We identified eight blood pressure-related circulating phospholipids that have a predictive value for incident hypertension. Our cross-omics analyses show that phospholipid metabolites in the circulation may yield insight into blood pressure regulation and raise a number of testable hypothesis for future research.

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

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

Genetic Studies of Leptin Concentrations Implicate Leptin in the Regulation of Early Adiposity.

Leptin influences food intake by informing the brain about the status of body fat stores. Rare LEP mutations associated with congenital leptin deficiency cause severe early-onset obesity that can be mitigated by administering leptin. However, the role of genetic regulation of leptin in polygenic obesity remains poorly understood. We performed an exome-based analysis in up to 57,232 individuals of diverse ancestries to identify genetic variants that influence adiposity-adjusted leptin concentrations. We identify five novel variants, including four missense variants, in LEP, ZNF800, KLHL31, and ACTL9, and one intergenic variant near KLF14. The missense variant Val94Met (rs17151919) in LEP was common in individuals of African ancestry only, and its association with lower leptin concentrations was specific to this ancestry (P = 2 × 10-16, n = 3,901). Using in vitro analyses, we show that the Met94 allele decreases leptin secretion. We also show that the Met94 allele is associated with higher BMI in young African-ancestry children but not in adults, suggesting that leptin regulates early adiposity.

CDH6 and HAGH protein levels in plasma associate with Alzheimer's disease in APOE ε4 carriers.

Many Alzheimer's disease (AD) genes including Apolipoprotein E (APOE) are found to be expressed in blood-derived macrophages and thus may alter blood protein levels. We measured 91 neuro-proteins in plasma from 316 participants of the Rotterdam Study (incident AD = 161) using Proximity Extension Ligation assay. We studied the association of plasma proteins with AD in the overall sample and stratified by APOE. Findings from the Rotterdam study were replicated in 186 AD patients of the BioFINDER study. We further evaluated the correlation of these protein biomarkers with total tau (t-tau), phosphorylated tau (p-tau) and amyloid-beta (Aß) 42 levels in cerebrospinal fluid (CSF) in the Amsterdam Dementia Cohort (N = 441). Finally, we conducted a genome-wide association study (GWAS) to identify the genetic variants determining the blood levels of AD-associated proteins. Plasma levels of the proteins, CDH6 (ß = 0.638, P = 3.33 × 10-4) and HAGH (ß = 0.481, P = 7.20 × 10-4), were significantly elevated in APOE ε4 carrier AD patients. The findings in the Rotterdam Study were replicated in the BioFINDER study for both CDH6 (ß = 1.365, P = 3.97 × 10-3) and HAGH proteins (ß = 0.506, P = 9.31 × 10-7) when comparing cases and controls in APOE ε4 carriers. In the CSF, CDH6 levels were positively correlated with t-tau and p-tau in the total sample as well as in APOE ε4 stratum (P < 1 × 10-3). The HAGH protein was not detected in CSF. GWAS of plasma CDH6 protein levels showed significant association with a cis-regulatory locus (rs111283466, P = 1.92 × 10-9). CDH6 protein is implicated in cell adhesion and synaptogenesis while HAGH protein is related to the oxidative stress pathway. Our findings suggest that these pathways may be altered during presymptomatic AD and that CDH6 and HAGH may be new blood-based biomarkers.

Smoking-by-genotype interaction in type 2 diabetes risk and fasting glucose.

Smoking is a potentially causal behavioral risk factor for type 2 diabetes (T2D), but not all smokers develop T2D. It is unknown whether genetic factors partially explain this variation. We performed genome-environment-wide interaction studies to identify loci exhibiting potential interaction with baseline smoking status (ever vs. never) on incident T2D and fasting glucose (FG). Analyses were performed in participants of European (EA) and African ancestry (AA) separately. Discovery analyses were conducted using genotype data from the 50,000-single-nucleotide polymorphism (SNP) ITMAT-Broad-CARe (IBC) array in 5 cohorts from from the Candidate Gene Association Resource Consortium (n = 23,189). Replication was performed in up to 16 studies from the Cohorts for Heart Aging Research in Genomic Epidemiology Consortium (n = 74,584). In meta-analysis of discovery and replication estimates, 5 SNPs met at least one criterion for potential interaction with smoking on incident T2D at p<1x10-7 (adjusted for multiple hypothesis-testing with the IBC array). Two SNPs had significant joint effects in the overall model and significant main effects only in one smoking stratum: rs140637 (FBN1) in AA individuals had a significant main effect only among smokers, and rs1444261 (closest gene C2orf63) in EA individuals had a significant main effect only among nonsmokers. Three additional SNPs were identified as having potential interaction by exhibiting a significant main effects only in smokers: rs1801232 (CUBN) in AA individuals, rs12243326 (TCF7L2) in EA individuals, and rs4132670 (TCF7L2) in EA individuals. No SNP met significance for potential interaction with smoking on baseline FG. The identification of these loci provides evidence for genetic interactions with smoking exposure that may explain some of the heterogeneity in the association between smoking and T2D.

Integration of epidemiologic, pharmacologic, genetic and gut microbiome data in a drug-metabolite atlas.

Progress in high-throughput metabolic profiling provides unprecedented opportunities to obtain insights into the effects of drugs on human metabolism. The Biobanking BioMolecular Research Infrastructure of the Netherlands has constructed an atlas of drug-metabolite associations for 87 commonly prescribed drugs and 150 clinically relevant plasma-based metabolites assessed by proton nuclear magnetic resonance. The atlas includes a meta-analysis of ten cohorts (18,873 persons) and uncovers 1,071 drug-metabolite associations after evaluation of confounders including co-treatment. We show that the effect estimates of statins on metabolites from the cross-sectional study are comparable to those from intervention and genetic observational studies. Further data integration links proton pump inhibitors to circulating metabolites, liver function, hepatic steatosis and the gut microbiome. Our atlas provides a tool for targeted experimental pharmaceutical research and clinical trials to improve drug efficacy, safety and repurposing. We provide a web-based resource for visualization of the atlas (http://bbmri.researchlumc.nl/atlas/).

Metabolomics Profile in Depression: A Pooled Analysis of 230 Metabolic Markers in 5283 Cases With Depression and 10,145 Controls.

BACKGROUND: Depression has been associated with metabolic alterations, which adversely impact cardiometabolic health. Here, a comprehensive set of metabolic markers, predominantly lipids, was compared between depressed and nondepressed persons. METHODS: Nine Dutch cohorts were included, comprising 10,145 control subjects and 5283 persons with depression, established with diagnostic interviews or questionnaires. A proton nuclear magnetic resonance metabolomics platform provided 230 metabolite measures: 51 lipids, fatty acids, and low-molecular-weight metabolites; 98 lipid composition and particle concentration measures of lipoprotein subclasses; and 81 lipid and fatty acids ratios. For each metabolite measure, logistic regression analyses adjusted for gender, age, smoking, fasting status, and lipid-modifying medication were performed within cohort, followed by random-effects meta-analyses. RESULTS: Of the 51 lipids, fatty acids, and low-molecular-weight metabolites, 21 were significantly related to depression (false discovery rate q < .05). Higher levels of apolipoprotein B, very-low-density lipoprotein cholesterol, triglycerides, diglycerides, total and monounsaturated fatty acids, fatty acid chain length, glycoprotein acetyls, tyrosine, and isoleucine and lower levels of high-density lipoprotein cholesterol, acetate, and apolipoprotein A1 were associated with increased odds of depression. Analyses of lipid composition indicators confirmed a shift toward less high-density lipoprotein and more very-low-density lipoprotein and triglyceride particles in depression. Associations appeared generally consistent across gender, age, and body mass index strata and across cohorts with depressive diagnoses versus symptoms. CONCLUSIONS: This large-scale meta-analysis indicates a clear distinctive profile of circulating lipid metabolites associated with depression, potentially opening new prevention or treatment avenues for depression and its associated cardiometabolic comorbidity.

Lipidomic profiling identifies signatures of metabolic risk.

BACKGROUND: Metabolic syndrome (MetS), the clustering of metabolic risk factors, is associated with cardiovascular disease risk. We sought to determine if dysregulation of the lipidome may contribute to metabolic risk factors. METHODS: We measured 154 circulating lipid species in 658 participants from the Framingham Heart Study (FHS) using liquid chromatography-tandem mass spectrometry and tested for associations with obesity, dysglycemia, and dyslipidemia. Independent external validation was sought in three independent cohorts. Follow-up data from the FHS were used to test for lipid metabolites associated with longitudinal changes in metabolic risk factors. RESULTS: Thirty-nine lipids were associated with obesity and eight with dysglycemia in the FHS. Of 32 lipids that were available for replication for obesity and six for dyslipidemia, 28 (88%) replicated for obesity and five (83%) for dysglycemia. Four lipids were associated with longitudinal changes in body mass index and four were associated with changes in fasting blood glucose in the FHS. CONCLUSIONS: We identified and replicated several novel lipid biomarkers of key metabolic traits. The lipid moieties identified in this study are involved in biological pathways of metabolic risk and can be explored for prognostic and therapeutic utility.

A metabolic profile of all-cause mortality risk identified in an observational study of 44,168 individuals.

Predicting longer-term mortality risk requires collection of clinical data, which is often cumbersome. Therefore, we use a well-standardized metabolomics platform to identify metabolic predictors of long-term mortality in the circulation of 44,168 individuals (age at baseline 18-109), of whom 5512 died during follow-up. We apply a stepwise (forward-backward) procedure based on meta-analysis results and identify 14 circulating biomarkers independently associating with all-cause mortality. Overall, these associations are similar in men and women and across different age strata. We subsequently show that the prediction accuracy of 5- and 10-year mortality based on a model containing the identified biomarkers and sex (C-statistic = 0.837 and 0.830, respectively) is better than that of a model containing conventional risk factors for mortality (C-statistic = 0.772 and 0.790, respectively). The use of the identified metabolic profile as a predictor of mortality or surrogate endpoint in clinical studies needs further investigation.

Metabolomics reveals a link between homocysteine and lipid metabolism and leukocyte telomere length: the ENGAGE consortium.

Telomere shortening has been associated with multiple age-related diseases such as cardiovascular disease, diabetes, and dementia. However, the biological mechanisms responsible for these associations remain largely unknown. In order to gain insight into the metabolic processes driving the association of leukocyte telomere length (LTL) with age-related diseases, we investigated the association between LTL and serum metabolite levels in 7,853 individuals from seven independent cohorts. LTL was determined by quantitative polymerase chain reaction and the levels of 131 serum metabolites were measured with mass spectrometry in biological samples from the same blood draw. With partial correlation analysis, we identified six metabolites that were significantly associated with LTL after adjustment for multiple testing: lysophosphatidylcholine acyl C17:0 (lysoPC a C17:0, p-value = 7.1 × 10-6), methionine (p-value = 9.2 × 10-5), tyrosine (p-value = 2.1 × 10-4), phosphatidylcholine diacyl C32:1 (PC aa C32:1, p-value = 2.4 × 10-4), hydroxypropionylcarnitine (C3-OH, p-value = 2.6 × 10-4), and phosphatidylcholine acyl-alkyl C38:4 (PC ae C38:4, p-value = 9.0 × 10-4). Pathway analysis showed that the three phosphatidylcholines and methionine are involved in homocysteine metabolism and we found supporting evidence for an association of lipid metabolism with LTL. In conclusion, we found longer LTL associated with higher levels of lysoPC a C17:0 and PC ae C38:4, and with lower levels of methionine, tyrosine, PC aa C32:1, and C3-OH. These metabolites have been implicated in inflammation, oxidative stress, homocysteine metabolism, and in cardiovascular disease and diabetes, two major drivers of morbidity and mortality.