Hallmarks of environmental insults.

Environmental insults impair human health around the world. Contaminated air, water, soil, food, and occupational and household settings expose humans of all ages to a plethora of chemicals and environmental stressors. We propose eight hallmarks of environmental insults that jointly underpin the damaging impact of environmental exposures during the lifespan. Specifically, they include oxidative stress and inflammation, genomic alterations and mutations, epigenetic alterations, mitochondrial dysfunction, endocrine disruption, altered intercellular communication, altered microbiome communities, and impaired nervous system function. They provide a framework to understand why complex mixtures of environmental exposures induce severe health effects even at relatively modest concentrations.

Loci for insulin processing and secretion provide insight into type 2 diabetes risk.

Insulin secretion is critical for glucose homeostasis, and increased levels of the precursor proinsulin relative to insulin indicate pancreatic islet beta-cell stress and insufficient insulin secretory capacity in the setting of insulin resistance. We conducted meta-analyses of genome-wide association results for fasting proinsulin from 16 European-ancestry studies in 45,861 individuals. We found 36 independent signals at 30 loci (p value < 5 × 10-8), which validated 12 previously reported loci for proinsulin and ten additional loci previously identified for another glycemic trait. Half of the alleles associated with higher proinsulin showed higher rather than lower effects on glucose levels, corresponding to different mechanisms. Proinsulin loci included genes that affect prohormone convertases, beta-cell dysfunction, vesicle trafficking, beta-cell transcriptional regulation, and lysosomes/autophagy processes. We colocalized 11 proinsulin signals with islet expression quantitative trait locus (eQTL) data, suggesting candidate genes, including ARSG, WIPI1, SLC7A14, and SIX3. The NKX6-3/ANK1 proinsulin signal colocalized with a T2D signal and an adipose ANK1 eQTL signal but not the islet NKX6-3 eQTL. Signals were enriched for islet enhancers, and we showed a plausible islet regulatory mechanism for the lead signal in the MADD locus. These results show how detailed genetic studies of an intermediate phenotype can elucidate mechanisms that may predispose one to disease.

Proteogenomic Data Integration Reveals CXCL10 as a Potentially Downstream Causal Mediator for IL-6 Signaling on Atherosclerosis.

BACKGROUND: Genetic and experimental studies support a causal involvement of IL-6 (interleukin-6) signaling in atheroprogression. Although trials targeting IL-6 signaling are underway, any benefits must be balanced against an impaired host immune response. Dissecting the mechanisms that mediate the effects of IL-6 signaling on atherosclerosis could offer insights about novel drug targets with more specific effects. METHODS: Leveraging data from 522 681 individuals, we constructed a genetic instrument of 26 variants in the gene encoding the IL-6R (IL-6 receptor) that proxied for pharmacological IL-6R inhibition. Using Mendelian randomization, we assessed its effects on 3281 plasma proteins quantified with an aptamer-based assay in the INTERVAL cohort (n=3301). Using mediation Mendelian randomization, we explored proteomic mediators of the effects of genetically proxied IL-6 signaling on coronary artery disease, large artery atherosclerotic stroke, and peripheral artery disease. For significant mediators, we tested associations of their circulating levels with incident cardiovascular events in a population-based study (n=1704) and explored the histological, transcriptomic, and cellular phenotypes correlated with their expression levels in samples from human atherosclerotic lesions. RESULTS: We found significant effects of genetically proxied IL-6 signaling on 70 circulating proteins involved in cytokine production/regulation and immune cell recruitment/differentiation, which correlated with the proteomic effects of pharmacological IL-6R inhibition in a clinical trial. Among the 70 significant proteins, genetically proxied circulating levels of CXCL10 (C-X-C motif chemokine ligand 10) were associated with risk of coronary artery disease, large artery atherosclerotic stroke, and peripheral artery disease, with up to 67% of the effects of genetically downregulated IL-6 signaling on these end points mediated by decreases in CXCL10. Higher midlife circulating CXCL10 levels were associated with a larger number of cardiovascular events over 20 years, whereas higher CXCL10 expression in human atherosclerotic lesions correlated with a larger lipid core and a transcriptomic profile reflecting immune cell infiltration, adaptive immune system activation, and cytokine signaling. CONCLUSIONS: Integrating multiomics data, we found a proteomic signature of IL-6 signaling activation and mediators of its effects on cardiovascular disease. Our analyses suggest the interferon-γ-inducible chemokine CXCL10 to be a potentially causal mediator for atherosclerosis in 3 vascular compartments and, as such, could serve as a promising drug target for atheroprotection.

Differences and commonalities in the genetic architecture of protein quantitative trait loci in European and Arab populations.

Polygenic scores (PGS) can identify individuals at risk of adverse health events and guide genetics-based personalized medicine. However, it is not clear how well PGS translate between different populations, limiting their application to well-studied ethnicities. Proteins are intermediate traits linking genetic predisposition and environmental factors to disease, with numerous blood circulating protein levels representing functional readouts of disease-related processes. We hypothesized that studying the genetic architecture of a comprehensive set of blood-circulating proteins between a European and an Arab population could shed fresh light on the translatability of PGS to understudied populations. We therefore conducted a genome-wide association study with whole-genome sequencing data using 1301 proteins measured on the SOMAscan aptamer-based affinity proteomics platform in 2935 samples of Qatar Biobank and evaluated the replication of protein quantitative traits (pQTLs) from European studies in an Arab population. Then, we investigated the colocalization of shared pQTL signals between the two populations. Finally, we compared the performance of protein PGS derived from a Caucasian population in a European and an Arab cohort. We found that the majority of shared pQTL signals (81.8%) colocalized between both populations. About one-third of the genetic protein heritability was explained by protein PGS derived from a European cohort, with protein PGS performing ~20% better in Europeans when compared to Arabs. Our results are relevant for the translation of PGS to non-Caucasian populations, as well as for future efforts to extend genetic research to understudied populations.

Nocturnal heat exposure and stroke risk.

BACKGROUND AND AIMS: In recent decades, nighttime temperatures have increased faster than daytime temperatures. The increasing prevalence of nocturnal heat exposure may pose a significant risk to cardiovascular health. This study investigated the association between nighttime heat exposure and stroke risk in the region of Augsburg, Germany, and examined its temporal variations over 15 years. METHODS: Hourly meteorological parameters, including mean temperature, relative humidity, and barometric pressure, were acquired from a local meteorological station. A data set was obtained consisting of 11 037 clinical stroke cases diagnosed during warmer months (May to October) between the years 2006 and 2020. The average age of cases was 71.3 years. Among these cases, 642 were identified as haemorrhagic strokes, 7430 were classified as ischaemic strokes, and 2947 were transient ischaemic attacks. A time-stratified case-crossover analysis with a distributed lag non-linear model was used to estimate the stroke risk associated with extreme nighttime heat, as measured by the hot night excess (HNE) index after controlling for the potential confounding effects of daily maximum temperature and other climatic variables. Subgroup analyses by age group, sex, stroke subtype, and stroke severity were performed to identify variations in susceptibility to nighttime heat. RESULTS: Results suggested a significant increase in stroke risk on days with extreme nighttime heat (97.5% percentile of HNE) (odds ratio 1.07, 95% confidence interval 1.01-1.15) during the full study period. When comparing the results for 2013-20 with the results for 2006-12, there was a significant increase (P < .05) in HNE-related risk for all strokes and specifically for ischaemic strokes during the more recent period. Furthermore, older individuals, females, and patients with mild stroke symptoms exhibited a significantly increased vulnerability to nighttime heat. CONCLUSIONS: This study found nocturnal heat exposure to be related to elevated stroke risk after controlling for maximum daytime temperature, with increasing susceptibility between 2006 and 2020. These results underscore the importance of considering nocturnal heat as a critical trigger of stroke events in a warming climate.

C-reactive protein modifies lipoprotein(a)-related risk for coronary heart disease: the BiomarCaRE project.

BACKGROUND AND AIMS: Recent investigations have suggested an interdependence of lipoprotein(a) [Lp(a)]-related risk for cardiovascular disease with background inflammatory burden. The aim the present analysis was to investigate whether high-sensitive C-reactive protein (hsCRP) modulates the association between Lp(a) and coronary heart disease (CHD) in the general population. METHODS: Data from 71 678 participants from 8 European prospective population-based cohort studies were used (65 661 without/6017 with established CHD at baseline; median follow-up 9.8/13.8 years, respectively). Fine and Gray competing risk-adjusted models were calculated according to accompanying hsCRP concentration (<2 and ≥2 mg/L). RESULTS: Among CHD-free individuals, increased Lp(a) levels were associated with incident CHD irrespective of hsCRP concentration: fully adjusted sub-distribution hazard ratios [sHRs (95% confidence interval)] for the highest vs. lowest fifth of Lp(a) distribution were 1.45 (1.23-1.72) and 1.48 (1.23-1.78) for a hsCRP group of <2 and ≥2 mg/L, respectively, with no interaction found between these two biomarkers on CHD risk (Pinteraction = 0.82). In those with established CHD, similar associations were seen only among individuals with hsCRP ≥ 2 mg/L [1.34 (1.03-1.76)], whereas among participants with a hsCRP concentration <2 mg/L, there was no clear association between Lp(a) and future CHD events [1.29 (0.98-1.71)] (highest vs. lowest fifth, fully adjusted models; Pinteraction = 0.024). CONCLUSIONS: While among CHD-free individuals Lp(a) was significantly associated with incident CHD regardless of hsCRP, in participants with CHD at baseline, Lp(a) was related to recurrent CHD events only in those with residual inflammatory risk. These findings might guide adequate selection of high-risk patients for forthcoming Lp(a)-targeting compounds.

Population-based reference values for kidney function and kidney function decline in 25- to 95-year-old Germans without and with diabetes.

Understanding normal aging of kidney function is pivotal to help distinguish individuals at particular risk for chronic kidney disease. Glomerular filtration rate (GFR) is typically estimated via serum creatinine (eGFRcrea) or cystatin C (eGFRcys). Since population-based age-group-specific reference values for eGFR and eGFR-decline are scarce, we aimed to provide such reference values from population-based data of a wide age range. In four German population-based cohorts (KORA-3, KORA-4, AugUR, DIACORE), participants underwent medical exams, interview, and blood draw up to five times within up to 25 years. We analyzed eGFRcrea and eGFRcys cross-sectionally and longitudinally (12,000 individuals, age 25-95 years). Cross-sectionally, we found age-group-specific eGFRcrea to decrease approximately linearly across the full age range, for eGFRcys up to the age of 60 years. Within age-groups, there was little difference by sex or diabetes status. Longitudinally, linear mixed models estimated an annual eGFRcrea decline of -0.80 [95% confidence interval -0.82, -0.77], -0.79 [-0.83, -0.76], and -1.20 mL/min/1.73m2 [-1.33, -1.08] for the general population, "healthy" individuals, or individuals with diabetes, respectively. Reference values for eGFR using cross-sectional data were shown as percentile curves for "healthy" individuals and for individuals with diabetes. Reference values for eGFR-decline using longitudinal data were presented as 95% prediction intervals for "healthy" individuals and for individuals with diabetes, obesity, and/or albuminuria. Thus, our results can help clinicians to judge eGFR values in individuals seen in clinical practice according to their age and to understand the expected range of annual eGFR-decline based on their risk profile.

Associations of carotid intima media thickness with gene expression in whole blood and genetically predicted gene expression across 48 tissues.

Carotid intima media thickness (cIMT) is a biomarker of subclinical atherosclerosis and a predictor of future cardiovascular events. Identifying associations between gene expression levels and cIMT may provide insight to atherosclerosis etiology. Here, we use two approaches to identify associations between mRNA levels and cIMT: differential gene expression analysis in whole blood and S-PrediXcan. We used microarrays to measure genome-wide whole blood mRNA levels of 5647 European individuals from four studies. We examined the association of mRNA levels with cIMT adjusted for various potential confounders. Significant associations were tested for replication in three studies totaling 3943 participants. Next, we applied S-PrediXcan to summary statistics from a cIMT genome-wide association study (GWAS) of 71 128 individuals to estimate the association between genetically determined mRNA levels and cIMT and replicated these analyses using S-PrediXcan on an independent GWAS on cIMT that included 22 179 individuals from the UK Biobank. mRNA levels of TNFAIP3, CEBPD and METRNL were inversely associated with cIMT, but these associations were not significant in the replication analysis. S-PrediXcan identified associations between cIMT and genetically determined mRNA levels for 36 genes, of which six were significant in the replication analysis, including TLN2, which had not been previously reported for cIMT. There was weak correlation between our results using differential gene expression analysis and S-PrediXcan. Differential expression analysis and S-PrediXcan represent complementary approaches for the discovery of associations between phenotypes and gene expression. Using these approaches, we prioritize TNFAIP3, CEBPD, METRNL and TLN2 as new candidate genes whose differential expression might modulate cIMT.

Genetic and clinical determinants of abdominal aortic diameter: genome-wide association studies, exome array data and Mendelian randomization study.

Progressive dilation of the infrarenal aortic diameter is a consequence of the ageing process and is considered the main determinant of abdominal aortic aneurysm (AAA). We aimed to investigate the genetic and clinical determinants of abdominal aortic diameter (AAD). We conducted a meta-analysis of genome-wide association studies in 10 cohorts (n = 13 542) imputed to the 1000 Genome Project reference panel including 12 815 subjects in the discovery phase and 727 subjects [Partners Biobank cohort 1 (PBIO)] as replication. Maximum anterior-posterior diameter of the infrarenal aorta was used as AAD. We also included exome array data (n = 14 480) from seven epidemiologic studies. Single-variant and gene-based associations were done using SeqMeta package. A Mendelian randomization analysis was applied to investigate the causal effect of a number of clinical risk factors on AAD. In genome-wide association study (GWAS) on AAD, rs74448815 in the intronic region of LDLRAD4 reached genome-wide significance (beta = -0.02, SE = 0.004, P-value = 2.10 × 10-8). The association replicated in the PBIO1 cohort (P-value = 8.19 × 10-4). In exome-array single-variant analysis (P-value threshold = 9 × 10-7), the lowest P-value was found for rs239259 located in SLC22A20 (beta = 0.007, P-value = 1.2 × 10-5). In the gene-based analysis (P-value threshold = 1.85 × 10-6), PCSK5 showed an association with AAD (P-value = 8.03 × 10-7). Furthermore, in Mendelian randomization analyses, we found evidence for genetic association of pulse pressure (beta = -0.003, P-value = 0.02), triglycerides (beta = -0.16, P-value = 0.008) and height (beta = 0.03, P-value < 0.0001), known risk factors for AAA, consistent with a causal association with AAD. Our findings point to new biology as well as highlighting gene regions in mechanisms that have previously been implicated in the genetics of other vascular diseases.

Association of plasma proteomics with mortality in individuals with and without type 2 diabetes: Results from two population-based KORA cohort studies.

BACKGROUND: Protein biomarkers may contribute to the identification of vulnerable subgroups for premature mortality. This study aimed to investigate the association of plasma proteins with all-cause and cause-specific mortality among individuals with and without baseline type 2 diabetes (T2D) and evaluate their impact on the prediction of all-cause mortality in two prospective Cooperative Health Research in the Region of Augsburg (KORA) studies. METHODS: The discovery cohort comprised 1545 participants (median follow-up 15.6 years; 244 with T2D: 116 total, 62 cardiovascular, 31 cancer-related and 23 other-cause deaths; 1301 without T2D: 321 total, 114 cardiovascular, 120 cancer-related and 87 other-cause deaths). The validation cohort comprised 1031 participants (median follow-up 6.9 years; 203 with T2D: 76 total, 45 cardiovascular, 19 cancer-related and 12 other-cause deaths; 828 without T2D: 169 total, 74 cardiovascular, 39 cancer-related and 56 other-cause deaths). We used Cox regression to examine associations of 233 plasma proteins with all-cause and cause-specific mortality and Lasso regression to construct prediction models for all-cause mortality stratifying by baseline T2D. C-index, category-free net reclassification index (cfNRI), and integrated discrimination improvement (IDI) were conducted to evaluate the predictive performance of built prediction models. RESULTS: Thirty-five and 62 proteins, with 29 overlapping, were positively associated with all-cause mortality in the group with and without T2D, respectively. Out of these, in the group with T2D, 35, eight, and 26 were positively associated with cardiovascular, cancer-related, and other-cause mortality, while in the group without T2D, 55, 41, and 47 were positively associated with respective cause-specific outcomes in the pooled analysis of both cohorts. Regulation of insulin-like growth factor (IGF) transport and uptake by IGF-binding proteins emerged as a unique pathway enriched for all-cause and cardiovascular mortality in individuals with T2D. The combined model containing the selected proteins (five and 12 proteins, with four overlapping, in the group with and without T2D, respectively) and clinical risk factors improved the prediction of all-cause mortality by C-index, cfNRI, and IDI. CONCLUSIONS: This study uncovered shared and unique mortality-related proteins in persons with and without T2D and emphasized the role of proteins in improving the prediction of mortality in different T2D subgroups.

TIGER: technical variation elimination for metabolomics data using ensemble learning architecture.

Large metabolomics datasets inevitably contain unwanted technical variations which can obscure meaningful biological signals and affect how this information is applied to personalized healthcare. Many methods have been developed to handle unwanted variations. However, the underlying assumptions of many existing methods only hold for a few specific scenarios. Some tools remove technical variations with models trained on quality control (QC) samples which may not generalize well on subject samples. Additionally, almost none of the existing methods supports datasets with multiple types of QC samples, which greatly limits their performance and flexibility. To address these issues, a non-parametric method TIGER (Technical variation elImination with ensemble learninG architEctuRe) is developed in this study and released as an R package (https://CRAN.R-project.org/package=TIGERr). TIGER integrates the random forest algorithm into an adaptable ensemble learning architecture. Evaluation results show that TIGER outperforms four popular methods with respect to robustness and reliability on three human cohort datasets constructed with targeted or untargeted metabolomics data. Additionally, a case study aiming to identify age-associated metabolites is performed to illustrate how TIGER can be used for cross-kit adjustment in a longitudinal analysis with experimental data of three time-points generated by different analytical kits. A dynamic website is developed to help evaluate the performance of TIGER and examine the patterns revealed in our longitudinal analysis (https://han-siyu.github.io/TIGER_web/). Overall, TIGER is expected to be a powerful tool for metabolomics data analysis.

Sex-specific trends in incidence of first myocardial infarction among people with and without diabetes between 1985 and 2016 in a German region.

BACKGROUND: The reduction of myocardial infarction (MI) and narrowing the gap between the populations with and without diabetes are important goals of diabetes care. We analyzed time trends for sex-specific incidence rates (IR) of first MI (both non-fatal MI and fatal MI) as well as separately for first non-fatal MI and fatal MI in the population with and without diabetes. METHODS: Using data from the KORA myocardial infarction registry (Augsburg, Germany), we estimated age-adjusted IR in people with and without diabetes, corresponding relative risks (RR), and time trends from 1985 to 2016 using Poisson regression. RESULTS: There were 19,683 people with first MI (34% fatal MI, 71% men, 30% with diabetes) between 1985 and 2016. In the entire study population, the IR of first MI decreased from 359 (95% CI: 345-374) to 236 (226-245) per 100,000 person years. In men with diabetes, IR decreased only in 2013-2016. This was due to first non-fatal MI, where IR in men with diabetes increased until 2009-2012, and slightly decreased in 2013-2016. Overall, fatal MI declined stronger than first non-fatal MI corresponding to IRs. The RR of first MI substantially increased among men from 1.40 (1.22-1.61) in 1985-1988 to 2.60 (2.26-2.99) in 1997-2000 and moderately decreased in 2013-2016: RR: 1.75 (1.47-2.09). Among women no consistent time trend for RR was observed. Time trends for RR were similar regarding first non-fatal MI and fatal MI. CONCLUSIONS: Over the study period, we found a decreased incidence of first MI and fatal MI in the entire study population. The initial increase of first non-fatal MI in men with diabetes needs further research. The gap between populations with and without diabetes remained.

Association of plasma proteomics with incident coronary heart disease in individuals with and without type 2 diabetes: results from the population-based KORA study.

BACKGROUND: Coronary heart disease (CHD) is a major global health concern, especially among individuals with type 2 diabetes (T2D). Given the crucial role of proteins in various biological processes, this study aimed to elucidate the aetiological role and predictive performance of protein biomarkers on incident CHD in individuals with and without T2D. METHODS: The discovery cohort included 1492 participants from the Cooperative Health Research in the Region of Augsburg (KORA) S4 study with 147 incident CHD cases (45 vs. 102 cases in the group with T2D and without T2D, respectively) during 15.6 years of follow-up. The validation cohort included 888 participants from the KORA-Age1 study with 70 incident CHD cases (19 vs. 51 cases in the group with T2D and without T2D, respectively) during 6.9 years of follow-up. We measured 233 plasma proteins related to cardiovascular disease and inflammation using proximity extension assay technology. Associations of proteins with incident CHD were assessed using Cox regression and Mendelian randomization (MR) analysis. Predictive models were developed using priority-Lasso and were evaluated on top of Framingham risk score variables using the C-index, category-free net reclassification index (cfNRI), and relative integrated discrimination improvement (IDI). RESULTS: We identified two proteins associated with incident CHD in individuals with and 29 in those without baseline T2D, respectively. Six of these proteins are novel candidates for incident CHD. MR suggested a potential causal role for hepatocyte growth factor in CHD development. The developed four-protein-enriched model for individuals with baseline T2D (ΔC-index: 0.017; cfNRI: 0.253; IDI: 0.051) and the 12-protein-enriched model for individuals without baseline T2D (ΔC-index: 0.054; cfNRI: 0.462; IDI: 0.024) consistently improved CHD prediction in the discovery cohort, while in the validation cohort, significant improvements were only observed for selected performance measures (with T2D: cfNRI: 0.633; without T2D: ΔC-index: 0.038; cfNRI: 0.465). CONCLUSIONS: This study identified novel protein biomarkers associated with incident CHD in individuals with and without T2D and reaffirmed previously reported protein candidates. These findings enhance our understanding of CHD pathophysiology and provide potential targets for prevention and treatment.

Plasma proteome association with coronary heart disease and carotid intima media thickness: results from the KORA F4 study.

BACKGROUND AND AIMS: Atherosclerosis is the main cause of stroke and coronary heart disease (CHD), both leading mortality causes worldwide. Proteomics, as a high-throughput method, could provide helpful insights into the pathological mechanisms underlying atherosclerosis. In this study, we characterized the associations of plasma protein levels with CHD and with carotid intima-media thickness (CIMT), as a surrogate measure of atherosclerosis. METHODS: The discovery phase included 1000 participants from the KORA F4 study, whose plasma protein levels were quantified using the aptamer-based SOMAscan proteomics platform. We evaluated the associations of plasma protein levels with CHD using logistic regression, and with CIMT using linear regression. For both outcomes we applied two models: an age-sex adjusted model, and a model additionally adjusted for body mass index, smoking status, physical activity, diabetes status, hypertension status, low density lipoprotein, high density lipoprotein, and triglyceride levels (fully-adjusted model). The replication phase included a matched case-control sample from the independent KORA F3 study, using ELISA-based measurements of galectin-4. Pathway analysis was performed with nominally associated proteins (p-value < 0.05) from the fully-adjusted model. RESULTS: In the KORA F4 sample, after Bonferroni correction, we found CHD to be associated with five proteins using the age-sex adjusted model: galectin-4 (LGALS4), renin (REN), cathepsin H (CTSH), and coagulation factors X and Xa (F10). The fully-adjusted model yielded only the positive association of galectin-4 (OR = 1.58, 95% CI = 1.30-1.93), which was successfully replicated in the KORA F3 sample (OR = 1.40, 95% CI = 1.09-1.88). For CIMT, we found four proteins to be associated using the age-sex adjusted model namely: cytoplasmic protein NCK1 (NCK1), insulin-like growth factor-binding protein 2 (IGFBP2), growth hormone receptor (GHR), and GDNF family receptor alpha-1 (GFRA1). After assessing the fully-adjusted model, only NCK1 remained significant (ß = 0.017, p-value = 1.39e-06). Upstream regulators of galectin-4 and NCK1 identified from pathway analysis were predicted to be involved in inflammation pathways. CONCLUSIONS: Our proteome-wide association study identified galectin-4 to be associated with CHD and NCK1 to be associated with CIMT. Inflammatory pathways underlying the identified associations highlight the importance of inflammation in the development and progression of CHD.

Bidirectional modulation of TCA cycle metabolites and anaplerosis by metformin and its combination with SGLT2i.

BACKGROUND: Metformin and sodium-glucose-cotransporter-2 inhibitors (SGLT2i) are cornerstone therapies for managing hyperglycemia in diabetes. However, their detailed impacts on metabolic processes, particularly within the citric acid (TCA) cycle and its anaplerotic pathways, remain unclear. This study investigates the tissue-specific metabolic effects of metformin, both as a monotherapy and in combination with SGLT2i, on the TCA cycle and associated anaplerotic reactions in both mice and humans. METHODS: Metformin-specific metabolic changes were initially identified by comparing metformin-treated diabetic mice (MET) with vehicle-treated db/db mice (VG). These findings were then assessed in two human cohorts (KORA and QBB) and a longitudinal KORA study of metformin-naïve patients with Type 2 Diabetes (T2D). We also compared MET with db/db mice on combination therapy (SGLT2i + MET). Metabolic profiling analyzed 716 metabolites from plasma, liver, and kidney tissues post-treatment, using linear regression and Bonferroni correction for statistical analysis, complemented by pathway analyses to explore the pathophysiological implications. RESULTS: Metformin monotherapy significantly upregulated TCA cycle intermediates such as malate, fumarate, and α-ketoglutarate (α-KG) in plasma, and anaplerotic substrates including hepatic glutamate and renal 2-hydroxyglutarate (2-HG) in diabetic mice. Downregulated hepatic taurine was also observed. The addition of SGLT2i, however, reversed these effects, such as downregulating circulating malate and α-KG, and hepatic glutamate and renal 2-HG, but upregulated hepatic taurine. In human T2D patients on metformin therapy, significant systemic alterations in metabolites were observed, including increased malate but decreased citrulline. The bidirectional modulation of TCA cycle intermediates in mice influenced key anaplerotic pathways linked to glutaminolysis, tumorigenesis, immune regulation, and antioxidative responses. CONCLUSION: This study elucidates the specific metabolic consequences of metformin and SGLT2i on the TCA cycle, reflecting potential impacts on the immune system. Metformin shows promise for its anti-inflammatory properties, while the addition of SGLT2i may provide liver protection in conditions like metabolic dysfunction-associated steatotic liver disease (MASLD). These observations underscore the importance of personalized treatment strategies.

Stool and blood metabolomics in the metabolic syndrome: a cross-sectional study.

INTRODUCTION/OBJECTIVES: Changes in the stool metabolome have been poorly studied in the metabolic syndrome (MetS). Moreover, few studies have explored the relationship of stool metabolites with circulating metabolites. Here, we investigated the associations between stool and blood metabolites, the MetS and systemic inflammation. METHODS: We analyzed data from 1,370 participants of the KORA FF4 study (Germany). Metabolites were measured by Metabolon, Inc. (untargeted) in stool, and using the AbsoluteIDQ® p180 kit (targeted) in blood. Multiple linear regression models, adjusted for dietary pattern, age, sex, physical activity, smoking status and alcohol intake, were used to estimate the associations of metabolites with the MetS, its components and high-sensitivity C-reactive protein (hsCRP) levels. Partial correlation and Multi-Omics Factor Analysis (MOFA) were used to investigate the relationship between stool and blood metabolites. RESULTS: The MetS was significantly associated with 170 stool and 82 blood metabolites. The MetS components with the highest number of associations were triglyceride levels (stool) and HDL levels (blood). Additionally, 107 and 27 MetS-associated metabolites (in stool and blood, respectively) showed significant associations with hsCRP levels. We found low partial correlation coefficients between stool and blood metabolites. MOFA did not detect shared variation across the two datasets. CONCLUSIONS: The MetS, particularly dyslipidemia, is associated with multiple stool and blood metabolites that are also associated with systemic inflammation. Further studies are necessary to validate our findings and to characterize metabolic alterations in the MetS. Although our analyses point to weak correlations between stool and blood metabolites, additional studies using integrative approaches are warranted.

Associations between multiple neurological biomarkers and distal sensorimotor polyneuropathy: KORA F4/FF4 study.

AIMS: The aim of this study was to assess associations between neurological biomarkers and distal sensorimotor polyneuropathy (DSPN). MATERIALS AND METHODS: Cross-sectional analyses were based on 1032 participants aged 61-82 years from the population-based KORA F4 survey, 177 of whom had DSPN at baseline. The prevalence of type 2 diabetes was 20%. Prospective analyses used data from 505 participants without DSPN at baseline, of whom 125 had developed DSPN until the KORA FF4 survey. DSPN was defined based on the examination part of the Michigan Neuropathy Screening Instrument. Serum levels of neurological biomarkers were measured using proximity extension assay technology. Associations between 88 biomarkers and prevalent or incident DSPN were estimated using Poisson regression with robust error variance and are expressed as risk ratios (RR) and 95% CI per 1-SD increase. Results were adjusted for multiple confounders and multiple testing using the Benjamini-Hochberg procedure. RESULTS: Higher serum levels of CTSC (cathepsin C; RR [95% CI] 1.23 (1.08; 1.39), pB-H = 0.044) and PDGFRα (platelet-derived growth factor receptor A; RR [95% CI] 1.21 (1.08; 1.35), pB-H = 0.044) were associated with prevalent DSPN in the total study sample. CDH3, JAM-B, LAYN, RGMA and SCARA5 were positively associated with DSPN in the diabetes subgroup, whereas GCP5 was positively associated with DSPN in people without diabetes (all pB-H for interaction <0.05). None of the biomarkers showed an association with incident DSPN (all pB-H>0.05). CONCLUSIONS: This study identified multiple novel associations between neurological biomarkers and prevalent DSPN, which may be attributable to functions of these proteins in neuroinflammation, neural development and myelination.

Association of transketolase polymorphisms with diabetic polyneuropathy in the general population: The KORA F4 study.

AIMS: We recently reported that genetic variability in the TKT gene encoding transketolase, a key enzyme in the pentose phosphate pathway, is associated with measures of diabetic sensorimotor polyneuropathy (DSPN) in recent-onset diabetes. Here, we aimed to substantiate these findings in a population-based KORA F4 study. MATERIALS AND METHODS: In this cross-sectional study, we assessed seven single nucleotide polymorphisms (SNPs) in the transketolase gene in 952 participants from the KORA F4 study with normal glucose tolerance (NGT; n = 394), prediabetes (n = 411), and type 2 diabetes (n = 147). DSPN was defined by the examination part of the Michigan Neuropathy Screening Instrument (MNSI) using the original MNSI > 2 cut-off and two alternative versions extended by touch/pressure perception (TPP) (MNSI > 3) and by TPP plus cold perception (MNSI > 4). RESULTS: After adjustment for sex, age, BMI, and HbA1c, in type 2 diabetes participants, four out of seven transketolase SNPs were associated with DSPN for all three MNSI versions (all p ≤ 0.004). The odds ratios of these associations increased with extending the MNSI score, for example, OR (95% CI) for SNP rs62255988 with MNSI > 2: 1.99 (1.16-3.41), MNSI > 3: 2.27 (1.26-4.09), and MNSI > 4: 4.78 (2.22-10.26); SNP rs9284890 with MNSI > 2: 2.43 (1.42-4.16), MNSI > 3: 3.46 (1.82-6.59), and MNSI > 4: 4.75 (2.15-10.51). In contrast, no associations were found between transketolase SNPs and the three MNSI versions in the NGT and prediabetes groups. CONCLUSIONS: The link of genetic variation in transketolase enzyme to diabetic polyneuropathy corroborated at the population level strengthens the concept suggesting an important role of pathways metabolising glycolytic intermediates in the evolution of diabetic polyneuropathy.

Circulating metabolites modulated by diet are associated with depression.

Metabolome reflects the interplay of genome and exposome at molecular level and thus can provide deep insights into the pathogenesis of a complex disease like major depression. To identify metabolites associated with depression we performed a metabolome-wide association analysis in 13,596 participants from five European population-based cohorts characterized for depression, and circulating metabolites using ultra high-performance liquid chromatography/tandem accurate mass spectrometry (UHPLC/MS/MS) based Metabolon platform. We tested 806 metabolites covering a wide range of biochemical processes including those involved in lipid, amino-acid, energy, carbohydrate, xenobiotic and vitamin metabolism for their association with depression. In a conservative model adjusting for life style factors and cardiovascular and antidepressant medication use we identified 8 metabolites, including 6 novel, significantly associated with depression. In individuals with depression, increased levels of retinol (vitamin A), 1-palmitoyl-2-palmitoleoyl-GPC (16:0/16:1) (lecithin) and mannitol/sorbitol and lower levels of hippurate, 4-hydroxycoumarin, 2-aminooctanoate (alpha-aminocaprylic acid), 10-undecenoate (11:1n1) (undecylenic acid), 1-linoleoyl-GPA (18:2) (lysophosphatidic acid; LPA 18:2) are observed. These metabolites are either directly food derived or are products of host and gut microbial metabolism of food-derived products. Our Mendelian randomization analysis suggests that low hippurate levels may be in the causal pathway leading towards depression. Our findings highlight putative actionable targets for depression prevention that are easily modifiable through diet interventions.

Association of serum magnesium with metabolic syndrome and the role of chronic kidney disease: A population-based cohort study with Mendelian randomization.

OBJECTIVES: To assess the association of serum magnesium with prevalent and incident metabolic syndrome (MetS) and its individual components in the general population and to examine any effect modification by chronic kidney disease (CKD) status. METHODS: We analysed longitudinal data from the population-based KORA F4/FF4 study, including 2996 participants (387 with CKD) for cross-sectional analysis and 1446 participants (88 with CKD) for longitudinal analysis. Associations with MetS, as well as single components of MetS, were assessed by adjusted regression models. Nonlinearity was tested by restricted cubic splines and analyses were stratified by CKD. Causality was evaluated by two-sample Mendelian randomization (MR). RESULTS: Serum magnesium (1 SD) was inversely associated with prevalent MetS (odds ratio [OR] 0.90, 95% confidence interval [CI] 0.83, 0.98). The association was more pronounced in individuals with CKD (OR 0.75, 95% CI 0.59, 0.94). Among MetS components, serum magnesium was negatively associated with elevated fasting glucose (OR 0.78, 95% CI 0.71, 0.88) and, again, this association was more pronounced in individuals with CKD (OR 0.67, 95% CI 0.53, 0.84). Serum magnesium was not associated with incident MetS or its components. Restricted cubic spline analysis revealed a significant nonlinear inverse relationship of serum magnesium with MetS and elevated fasting glucose. MR analysis suggested an inverse causal effect of serum magnesium on MetS (OR 0.91, 95% CI 0.85, 0.97). CONCLUSION: Serum magnesium is associated with prevalent, but not incident MetS, and this effect is stronger in individuals with CKD. MR analysis implies a potential, albeit weak, causal role of magnesium in MetS.