Whole-genome sequencing identifies common-to-rare variants associated with human blood metabolites

Multi-scale variational autoencoder for imputation of missing values in untargeted metabolomics using whole-genome sequencing data

BACKGROUND

Missing data is a common challenge in mass spectrometry-based metabolomics, which can lead to biased and incomplete analyses. The integration of whole-genome sequencing (WGS) data with metabolomics data has emerged as a promising approach to enhance the accuracy of data imputation in metabolomics studies.

METHOD

In this study, we propose a novel method that leverages the information from WGS data and reference metabolites to impute unknown metabolites. Our approach utilizes a multi-scale variational autoencoder to jointly model the burden score, polygenetic risk score (PGS), and linkage disequilibrium (LD) pruned single nucleotide polymorphisms (SNPs) for feature extraction and missing metabolomics data imputation. By learning the latent representations of both omics data, our method can effectively impute missing metabolomics values based on genomic information.

RESULTS

We evaluate the performance of our method on empirical metabolomics datasets with missing values and demonstrate its superiority compared to conventional imputation techniques. Using 35 template metabolites derived burden scores, PGS and LD-pruned SNPs, the proposed methods achieved R2-scores > 0.01 for 71.55 % of metabolites.

CONCLUSION

The integration of WGS data in metabolomics imputation not only improves data completeness but also enhances downstream analyses, paving the way for more comprehensive and accurate investigations of metabolic pathways and disease associations. Our findings offer valuable insights into the potential benefits of utilizing WGS data for metabolomics data imputation and underscore the importance of leveraging multi-modal data integration in precision medicine research.

Causal association between plasma metabolites and neurodegenerative diseases

BACKGROUND

Establishing causal relationships between metabolic biomarkers and neurodegenerative diseases such as Alzheimer's disease (AD) and Parkinson's disease (PD) is a challenge faced by observational studies. In this study, our aim was to investigate the causal associations between plasma metabolites and neurodegenerative diseases using Mendelian Randomization (MR) methods.

METHODS

We utilized genetic associations with 1400 plasma metabolic traits as exposures. We used large-scale genome-wide association study (GWAS) summary statistics for AD and PD as our discovery datasets. For validation, we performed repeated analyses using different GWAS datasets. The main statistical method employed was inverse variance-weighted (IVW). We also conducted enrichment pathway analysis for IVW-identified metabolites.

RESULTS

In the discovered dataset, there are a total of 69 metabolites (36 negatively, 33 positively) potentially associated with AD, and 47 metabolites (24 negatively， 23 positively) potentially associated with PD. Among these, 4 significant metabolites overlap with significant metabolites (PIVW < 0.05)in the validation dataset for AD, and 1 metabolite overlaps with significant metabolites in the validation dataset for PD. Three metabolites serve as common potential metabolic markers for both AD and PD, including Tryptophan betaine, Palmitoleoylcarnitine (C16:1), and X-23655 levels. Further pathway enrichment analysis suggests that the SLC-mediated transmembrane transport pathway, involving tryptophan betaine and carnitine metabolites, may represent potential intervention targets for treating AD and PD.

CONCLUSION

This study offers novel insights into the causal effects of plasma metabolites on degenerative diseases through the integration of genomics and metabolomics. The identification of metabolites and metabolic pathways linked to AD and PD enhances our comprehension of the underlying biological mechanisms and presents promising targets for future therapeutic interventions in AD and PD.

Genetics and precision genomics approaches to pulmonary hypertension

Considerable progress has been made in the genomics of pulmonary arterial hypertension (PAH) since the 6th World Symposium on Pulmonary Hypertension, with the identification of rare variants in several novel genes, as well as common variants that confer a modest increase in PAH risk. Gene and variant curation by an expert panel now provides a robust framework for knowing which genes to test and how to interpret variants in clinical practice. We recommend that genetic testing be offered to specific subgroups of symptomatic patients with PAH, and to children with certain types of group 3 pulmonary hypertension (PH). Testing of asymptomatic family members and the use of genetics in reproductive decision-making require the involvement of genetics experts. Large cohorts of PAH patients with biospecimens now exist and extension to non-group 1 PH has begun. However, these cohorts are largely of European origin; greater diversity will be essential to characterise the full extent of genomic variation contributing to PH risk and treatment responses. Other types of omics data are also being incorporated. Furthermore, to advance gene- and pathway-specific care and targeted therapies, gene-specific registries will be essential to support patients and their families and to lay the foundation for genetically informed clinical trials. This will require international outreach and collaboration between patients/families, clinicians and researchers. Ultimately, harmonisation of patient-derived biospecimens, clinical and omic information, and analytic approaches will advance the field.

The mediating role of metabolites between gut microbiome and Hirschsprung disease: a bidirectional two-step Mendelian randomization study

Background

Gut microbiome (GM) was observed to be associated with the incidence of Hirschsprung disease (HD). However, the effect and mechanism of GM in HD is still unclear. To investigate the relationship between GM and HD and the effect of metabolites as mediators, a bidirectional two-step Mendelian randomization (MR) study was conducted.

Methods

The study selected instrument variables (IVs) from summary-level genome-wide association studies (GWAS). The MiBioGen consortium provided the GWAS data for GM, while the GWAS data for metabolites and HD were obtained from the GWAS Catalog consortium. Two-sample MR analyses were performed to estimate bidirectional correlations between IVs associated with GM and HD. Then, genetic variants related to 1,400 metabolite traits were selected for further mediation analyses using the Product method.

Results

This study found that seven genus bacteria had a significant causal relationship with the incidence of HD but not vice versa. 27 metabolite traits were significantly correlated with HD. After combining the significant results, three significant GM-metabolites-HD lines have been identified. In the Peptococcus-Stearoyl sphingomyelin (d18:1/18:0)-HD line, the Stearoyl sphingomyelin (d18:1/18:0) levels showed a mediation proportion of 14.5%, while in the Peptococcus-lysine-HD line, the lysine levels had a mediation proportion of 12.9%. Additionally, in the Roseburia-X-21733-HD line, the X-21733 levels played a mediation proportion of 23.5%.

Conclusion

Our MR study indicates a protective effect of Peptococcus on HD risk that is partially mediated through serum levels of stearoyl sphingomyelin (d18:1/18:0) and lysine, and a risk effect of Roseburia on HD that is partially mediated by X-21733 levels. These findings could serve as novel biomarkers and therapeutic targets for HD.

Genetic Architecture and Analysis Practices of Circulating Metabolites in the NHLBI Trans-Omics for Precision Medicine (TOPMed) Program

Circulating metabolite levels partly reflect the state of human health and diseases, and can be impacted by genetic determinants. Hundreds of loci associated with circulating metabolites have been identified; however, most findings focus on predominantly European ancestry or single study analyses. Leveraging the rich metabolomics resources generated by the NHLBI Trans-Omics for Precision Medicine (TOPMed) Program, we harmonized and accessibly cataloged 1,729 circulating metabolites among 25,058 ancestrally-diverse samples. We provided recommendations for outlier and imputation handling to process metabolite data, as well as a general analytical framework. We further performed a pooled analysis following our practical recommendations and discovered 1,778 independent loci associated with 667 metabolites. Among 108 novel locus - metabolite pairs, we detected not only novel loci within previously implicated metabolite associated genes, but also novel genes (such as GAB3 and VSIG4 located in the X chromosome) that have putative roles in metabolic regulation. In the sex-stratified analysis, we revealed 85 independent locus-metabolite pairs with evidence of sexual dimorphism, including well-known metabolic genes such as FADS2 , D2HGDH , SUGP1 , UTG2B17 , strongly supporting the importance of exploring sex difference in the human metabolome. Taken together, our study depicted the genetic contribution to circulating metabolite levels, providing additional insight into the understanding of human health.

Epigenome-wide association study on the plasma metabolome suggests self-regulation of the glycine and serine pathway through DNA methylation

BACKGROUND

The plasma metabolome reflects the physiological state of various biological processes and can serve as a proxy for disease risk. Plasma metabolite variation, influenced by genetic and epigenetic mechanisms, can also affect the cellular microenvironment and blood cell epigenetics. The interplay between the plasma metabolome and the blood cell epigenome remains elusive. In this study, we performed an epigenome-wide association study (EWAS) of 1183 plasma metabolites in 693 participants from the LifeLines-DEEP cohort and investigated the causal relationships in DNA methylation-metabolite associations using bidirectional Mendelian randomization and mediation analysis.

RESULTS

After rigorously adjusting for potential confounders, including genetics, we identified five robust associations between two plasma metabolites (L-serine and glycine) and three CpG sites located in two independent genomic regions (cg14476101 and cg16246545 in PHGDH and cg02711608 in SLC1A5) at a false discovery rate of less than 0.05. Further analysis revealed a complex bidirectional relationship between plasma glycine/serine levels and DNA methylation. Moreover, we observed a strong mediating role of DNA methylation in the effect of glycine/serine on the expression of their metabolism/transport genes, with the proportion of the mediated effect ranging from 11.8 to 54.3%. This result was also replicated in an independent population-based cohort, the Rotterdam Study. To validate our findings, we conducted in vitro cell studies which confirmed the mediating role of DNA methylation in the regulation of PHGDH gene expression.

CONCLUSIONS

Our findings reveal a potential feedback mechanism in which glycine and serine regulate gene expression through DNA methylation.

The consumption of animal products is associated with plasma levels of alpha-aminoadipic acid (2-AAA)

BACKGROUND AND AIMS

The cardiometabolic disease-associated metabolite, alpha-aminoadipic acid (2-AAA) is formed from the breakdown of the essential dietary amino acid lysine. However, it was not known whether elevated plasma levels of 2-AAA are related to dietary nutrient intake. We aimed to determine whether diet is a determinant of circulating 2-AAA in healthy individuals, and whether 2-AAA is altered in response to dietary modification.

METHODS AND RESULTS

We investigated the association between 2-AAA and dietary nutrient intake in a cross-sectional study of healthy individuals (N = 254). We then performed a randomized cross-over dietary intervention trial to investigate the effect of lysine supplementation (1 week) on 2-AAA in healthy individuals (N = 40). We further assessed the effect of a vegetarian diet on 2-AAA in a short-term (4-day) dietary intervention trial in healthy omnivorous women (N = 35). We found that self-reported dietary intake of animal products, including meat, poultry, and seafood, was associated with higher plasma 2-AAA cross-sectionally (P < 0.0001). Supplementary dietary lysine (5g/day) caused no significant increase in plasma 2-AAA; however, plasma 2-AAA was altered by general dietary modification. Further, plasma 2-AAA was significantly reduced by a short-term vegetarian diet (P = 0.003).

CONCLUSION

We identified associations between plasma 2-AAA and consumption of animal products, which were validated in a vegetarian dietary intervention trial, but not in a trial designed to specifically increase the 2-AAA amino acid precursor lysine. Further studies are warranted to investigate whether implementation of a vegetarian diet improves cardiometabolic risk in individuals with elevated 2-AAA.

Rare variants at KCNJ2 are associated with LDL-cholesterol levels in a cross-population study

Leveraging whole genome sequencing data of 1751 individuals from the UK and 2587 Qatari subjects, we suggest here an association of rare variants mapping to the sour taste-associated gene KCNJ2 with reduced low-density lipoprotein cholesterol (LDL-C, P = 2.10 × 10-12) and with a 22% decreased dietary trans-fat intake. This study identifies a novel candidate rare locus for LDL-C, adding insights into the genetic architecture of a complex trait implicated in cardiovascular disease.

Elucidating causal relationships of diet-derived circulating antioxidants and the risk of non-scarring alopecia: A Mendelian randomization study

Previous observational studies revealed controversy about the effect of circulating antioxidants on risk of alopecia. In the present study, we investigated the causal relationships between diet-derived circulating antioxidants and 2 non-scarring alopecia using Mendelian randomization (MR). Instrumental variables for antioxidants (lycopene, retinol, ascorbate, β-carotene, α-tocopherol, and γ-tocopherol) were selected from published studies. Data for alopecia areata (AA) and androgenetic alopecia (AGA) was obtained from the FinnGen study project (R9 released in 2023), including 195 cases and 201,019 controls for AGA and 682 cases and 361,140 controls for AA. We used the inverse variance weighted method as the primary MR method. Three additional methods were used as sensitivity analysis to validate the robustness of the results. We found a causal relationship between absolute β-carotene levels and AGA risk (P = .039), but not with AA (P = .283). The results of Wald ratio showed a protective effect of absolute β-carotene levels against AGA, with per 0.1 ln-transformed β-carotene being associated with a 76% lower risk of AGA (OR: 0.24, 95% CI: 0.06-0.93). Based on the fixed effects inverse variance weighting results, we found that α-tocopherol was protective against both AGA (P = .026) and AA (P = .018). For each unit increase in α-tocopherol, the effects of change in AGA and AA were 0.02 (95% CI: 0.00-0.61) and 0.10 (95% CI: 0.01-0.67), respectively. The results did not reveal any other causal relationships. Our study identified 3 causal associations of antioxidants with the risk of non-scarring alopecia. These results provide new insights into the prevention of non-scarring alopecia through diet.

Elucidating causal relationships of diet-derived circulating antioxidants and the risk of osteoporosis: A Mendelian randomization study

Background

Osteoporosis (OP) is typically diagnosed by evaluating bone mineral density (BMD), and it frequently results in fractures. Here, we investigated the causal relationships between diet-derived circulating antioxidants and the risk of OP using Mendelian randomization (MR).

Methods

Published studies were used to identify instrumental variables related to absolute levels of circulating antioxidants like lycopene, retinol, ascorbate, and β-carotene, as well as antioxidant metabolites such as ascorbate, retinol, α-tocopherol, and γ-tocopherol. Outcome variables included BMD (in femoral neck, lumbar spine, forearm, heel, total body, total body (age over 60), total body (age 45-60), total body (age 30-45), total body (age 15-30), and total body (age 0-15)), fractures (in arm, spine, leg, heel, and osteoporotic fractures), and OP. Inverse variance weighted or Wald ratio was chosen as the main method for MR analysis based on the number of single nucleotide polymorphisms (SNPs). Furthermore, we performed sensitivity analyses to confirm the reliability of the findings.

Results

We found a causal relationship between absolute retinol levels and heel BMD (p = 7.6E-05). The results of fixed effects IVW showed a protective effect of absolute retinol levels against heel BMD, with per 0.1 ln-transformed retinol being associated with a 28% increase in heel BMD (OR: 1.28, 95% CI: 1.13-1.44). In addition, a sex-specific effect of the absolute circulating retinol levels on the heel BMD has been observed in men. No other significant causal relationship was found.

Conclusion

There is a positive causal relationship between absolute retinol levels and heel BMD. The implications of our results should be taken into account in future studies and in the creation of public health policies and OP prevention tactics.

Exogenous and endogenous antioxidants in osteoporosis risk: causal associations unveiled by Mendelian Randomization analysis

Background

Recent epidemiological studies and animal experiments have highlighted the significant role of oxidative stress in the development of osteoporosis (OP). The provision of antioxidants is widely considered a fundamental strategy to combat free radical-induced stress, inhibit oxidative damage, and potentially reverse the adverse effects of oxidative stress on bone health. However, there is no consensus in the scientific literature regarding the practical effectiveness of antioxidants in OP prevention and treatment. Some studies have not shown a clear connection between antioxidant supplementation and decreased OP risk. Therefore, it is essential to clarify the potential causal relationship between antioxidants and the development of OP.

Methods

The study utilized the inverse variance weighted (IVW) approach as the primary analytical method in the Mendelian Randomization (MR) framework to investigate the causal effects of five exogenous and six endogenous antioxidants on the risk of OP. To thoroughly assess potential pleiotropic effects and heterogeneity among the data analyzed, the MR-Egger intercept test was employed, and Cochran's Q statistic was calculated.

Results

In the evaluation of exogenous antioxidants, single-directional two-sample MR analyses did not reveal any statistically significant relationship between these agents and the risk of OP. Regarding endogenous antioxidants, bidirectional two-sample MR analyses were conducted, which generally indicated that most genetically regulated endogenous antioxidants had no significant association with the onset risk of OP. A significant causal relationship was found between OP and serum albumin levels (β: -0.0552, 95%CI: -0.0879 to -0.0225, p < 0.0011 after Bonferroni adjustment, power = 100%).

Conclusion

The research uncovers OP as a possible determinant contributing to a decrement in serum albumin levels, and further suggests a potentially intimate relationship between the downward trajectory of serum albumin concentrations and the advancement of the OP disease process.

Metabolome-wide Mendelian randomization for age at menarche and age at natural menopause

BACKGROUND

The role of metabolism in the variation of age at menarche (AAM) and age at natural menopause (ANM) in the female population is not entirely known. We aimed to investigate the causal role of circulating metabolites in AAM and ANM using Mendelian randomization (MR).

METHODS

We combined MR with genetic colocalization to investigate potential causal associations between 658 metabolites and AAM and between 684 metabolites and ANM. We extracted genetic instruments for our exposures from four genome-wide association studies (GWAS) on circulating metabolites and queried the effects of these variants on the outcomes in two large GWAS from the ReproGen consortium. Additionally, we assessed the mediating role of the body mass index (BMI) in these associations, identified metabolic pathways implicated in AAM and ANM, and sought validation for selected metabolites in the Avon Longitudinal Study of Parents and Children (ALSPAC).

RESULTS

Our analysis identified 10 candidate metabolites for AAM, but none of them colocalized with AAM. For ANM, 76 metabolites were prioritized (FDR-adjusted MR P-value ≤ 0.05), with 17 colocalizing, primarily in the glycerophosphocholines class, including the omega-3 fatty acid and phosphatidylcholine (PC) categories. Pathway analyses and validation in ALSPAC mothers also highlighted the role of omega and polyunsaturated fatty acids levels in delaying age at menopause.

CONCLUSIONS

Our study suggests that metabolites from the glycerophosphocholine and fatty acid families play a causal role in the timing of both menarche and menopause. This underscores the significance of specific metabolic pathways in the biology of female reproductive longevity.

The full spectrum of SLC22 OCT1 mutations illuminates the bridge between drug transporter biophysics and pharmacogenomics

Mutations in transporters can impact an individual's response to drugs and cause many diseases. Few variants in transporters have been evaluated for their functional impact. Here, we combine saturation mutagenesis and multi-phenotypic screening to dissect the impact of 11,213 missense single-amino-acid deletions, and synonymous variants across the 554 residues of OCT1, a key liver xenobiotic transporter. By quantifying in parallel expression and substrate uptake, we find that most variants exert their primary effect on protein abundance, a phenotype not commonly measured alongside function. Using our mutagenesis results combined with structure prediction and molecular dynamic simulations, we develop accurate structure-function models of the entire transport cycle, providing biophysical characterization of all known and possible human OCT1 polymorphisms. This work provides a complete functional map of OCT1 variants along with a framework for integrating functional genomics, biophysical modeling, and human genetics to predict variant effects on disease and drug efficacy.

Intake of different types of seafood and meat and risk of type 2 diabetes in women: a prospective study supported by a dietary intervention in mice

Detailed knowledge regarding the associations between intake of different types of seafood and meat and the risk of type 2 diabetes (T2D), and insight into possible mechanisms are warranted. In this study we aimed to evaluate the associations between intake of different types of seafood and meat and the subsequent risk of T2D using the Norwegian Mother, Father, and Child Cohort Study (MoBa), and furthermore, by using a mouse model to gain further insight into possible molecular mechanisms contributing to the associated metabolic changes. Women in MoBa who were free of pharmacologically treated diabetes at baseline (n = 60,777) were prospectively evaluated for incident T2D, identified on the basis of medication usages > 90 days after delivery, ascertained by the Norwegian Prescription Database. Dietary intake was obtained with a validated 255-item food frequency questionnaire which assessed habitual diet during the first 4-5 months of pregnancy. Metabolic phenotypes and plasma metabolome were investigated in female mice fed isocaloric diets with different types of seafood and meat mimicking the dietary intake in the human cohort. During maximum 10-year and mean (SD) 7.2 (1.6) years follow-up time, 681 (1.1%) women developed pharmacologically treated T2D. All statistical models identified a higher risk of T2D with increased shellfish intake, whereas no associations were observed for total seafood, fatty fish, total meat and red meat in the adjusted models. In mice, the shellfish-based western diet induced reduced glucose tolerance and insulin secretion compared to the diet based on lean fish, and we identified a number of metabolites elevated in plasma from shellfish-fed mice that correlated with glucose intolerance. Mice fed a western diet based on meat also exhibited reduced glucose tolerance in comparison to lean fish fed mice, whereas mice fed fatty fish, total seafood or red meat did not differ from lean fish fed mice. We observed a diet-specific metabolic signature in plasma demonstrating five distinct metabolite profiles in mice fed shellfish, fatty fish, total seafood/lean fish, a mixed diet and meat. In conclusion, these findings demonstrate that different types of seafood have different outcome on T2D risk. In women, intake of shellfish was associated with higher risk of T2D. In female mice, a shellfish enriched diet reduced glucose tolerance and altered the abundance of several distinct plasma metabolites correlating with glucose tolerance.

Essential Nutrients and White Matter Hyperintensities: A Two-Sample Mendelian Randomization Study

Stroke and dementia have been linked to the appearance of white matter hyperintensities (WMHs). Meanwhile, diffusion tensor imaging (DTI) might capture the microstructural change in white matter early. Specific dietary interventions may help to reduce the risk of WMHs. However, research on the relationship between specific nutrients and white matter changes is still lacking. We aimed to investigate the causal effects of essential nutrients (amino acids, fatty acids, mineral elements, and vitamins) on WMHs and DTI measures, including fraction anisotropy (FA) and mean diffusivity (MD), by a Mendelian randomization analysis. We selected single nucleotide polymorphisms (SNPs) associated with each nutrient as instrumental variables to assess the causal effects of nutrient-related exposures on WMHs, FA, and MD. The outcome was from a recently published large-scale European Genome Wide Association Studies pooled dataset, including WMHs (N = 18,381), FA (N = 17,663), and MD (N = 17,467) data. We used the inverse variance weighting (IVW) method as the primary method, and sensitivity analyses were conducted using the simple median, weighted median, and MR-Egger methods. Genetically predicted serum calcium level was positively associated with WMHs risk, with an 8.1% increase in WMHs risk per standard deviation unit increase in calcium concentration (OR = 1.081, 95% CI = 1.006-1.161, p = 0.035). The plasma linoleic acid level was negatively associated with FA (OR = 0.776, 95% CI = 0.616-0.978, p = 0.032). Our study demonstrated that genetically predicted calcium was a potential risk factor for WMHs, and linoleic acid may be negatively associated with FA, providing evidence for interventions from the perspective of gene-environment interactions.

Diet-derived circulating antioxidants, periodontitis and dental caries: A Mendelian randomization study

BACKGROUND AND OBJECTIVE

Given the potential association between oxidative stress, periodontitis and dental caries, whether dietary supplementation with antioxidants is beneficial for periodontitis and dental caries has been widely reported, but remains controversial. This study aims to clarify these relationships through two-sample Mendelian randomization (MR) analysis.

METHODS

Circulating antioxidants (copper, selenium, zinc, ascorbate, β-carotene, lycopene, retinol and vitamin E) were derived from absolute circulating antioxidants and circulating antioxidant metabolites. Summary data of periodontitis and dental caries were obtained from two separate databases, respectively. We performed inverse-variance weighted (IVW) analysis separately in different databases, followed by meta-analysis. The robustness of results was examined by sensitivity analyses, including three complementary MR methods, heterogeneity and pleiotropy tests, and PhenoScanner query.

RESULTS

IVW analysis showed that elevated levels of absolute circulating retinol reduced the risk of periodontitis (GLIDE: OR = 0.41, 95% CI = 0.18-0.95, p = .038, power = 100%; FinnGen: OR = 0.15, 95% CI = 0.04-0.54, p = .004, power = 100%). The pooled OR for periodontitis risk per unit increase of retinol is 0.30 (95% CI = 0.15-0.61, p = .001, I2 = 40.3%, power = 100%). No significant associations were noted for genetically predicted circulating antioxidants and dental caries risk. The sensitivity analyses yielded similar estimates.

CONCLUSION

This study demonstrates that a negative causality between circulating retinol and periodontitis risk, and null linkage between circulating antioxidants and dental caries risk, suggesting potential strategies for the prevention and control of periodontitis.

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.

Identification of blood metabolites associated with risk of Alzheimer's disease by integrating genomics and metabolomics data

Specific metabolites have been reported to be potentially associated with Alzheimer's disease (AD) risk. However, the comprehensive understanding of roles of metabolite biomarkers in AD etiology remains elusive. We performed a large AD metabolome-wide association study (MWAS) by developing blood metabolite genetic prediction models. We evaluated associations between genetically predicted levels of metabolites and AD risk in 39,106 clinically diagnosed AD cases, 46,828 proxy AD and related dementia (proxy-ADD) cases, and 401,577 controls. We further conducted analyses to determine microbiome features associated with the detected metabolites and characterize associations between predicted microbiome feature levels and AD risk. We identified fourteen metabolites showing an association with AD risk. Five microbiome features were further identified to be potentially related to associations of five of the metabolites. Our study provides new insights into the etiology of AD that involves blood metabolites and gut microbiome, which warrants further investigation.

Diet-derived circulating antioxidants and risk of epilepsy: A study combining metabolomics and mendelian randomization

Background

Previous studies offer inconclusive results on the association between diet-derived circulating antioxidants and epilepsy.

Objective

This study aims to assess oxidative stress presence in epilepsy patients' circulation and investigate the causal link between diet-derived circulating antioxidants and epilepsy.

Methods

Untargeted metabolomics analysis was conducted on plasma samples from 62 epileptic patients and 20 healthy individuals to evaluate oxidative stress based on metabolite alterations in epilepsy patients' circulation. Two-sample Mendelian Randomization (MR) analysis examined the causation between diet-derived circulating antioxidants (measured by absolute levels and relative metabolite concentrations) and epilepsy, utilizing the inverse-variance weighted (IVW) method as the primary outcome, with complementary MR analysis methods (MR Egger, weighted median, weighted mode, and simple mode).

Results

Untargeted metabolomics analysis revealed elevated circulating oxidizing metabolites (palmitic acid, oleic acid, linoleic acid, and myristic acid) and reduced reducing metabolites (glutamine) in epilepsy patients, providing robust evidence of oxidative stress. The IVW analysis indicated significantly reduced epilepsy risk (odds ratio: 0.552; 95% confidence interval: 0.335-0.905, P = 0.018) with genetically determined higher absolute circulating β-carotene. However, other diet-derived circulating antioxidants (lycopene, retinol, ascorbic acid, and selenium) and antioxidant metabolites (α-tocopherol, γ-tocopherol, ascorbic acid, and retinol) did not significantly associate with epilepsy risk. Additional MR analysis methods and heterogeneity assessments confirmed the results' robustness.

Conclusion

This study provides compelling evidence of oxidative stress in epilepsy patients' circulation. However, the majority of diet-derived circulating antioxidants (lycopene, retinol, ascorbic acid, vitamin E, and selenium) are unlikely to causally associate with reduced epilepsy risk, except for β-carotene.

Prognostic Performance of Hematological and Serum Iron and Metabolite Indices for Detection of Early Iron Deficiency Induced Metabolic Brain Dysfunction in Infant Rhesus Monkeys

BACKGROUND

The current pediatric practice of monitoring for infantile iron deficiency (ID) via hemoglobin (Hgb) screening at one y of age does not identify preanemic ID nor protect against later neurocognitive deficits.

OBJECTIVES

To identify biomarkers of iron-related metabolic alterations in the serum and brain and determine the sensitivity of conventional iron and heme indices for predicting risk of brain metabolic dysfunction using a nonhuman primate model of infantile ID.

METHODS

Simultaneous serum iron and RBC indices, and serum and cerebrospinal fluid (CSF) metabolomic profiles were determined in 20 rhesus infants, comparing iron sufficient (IS; N = 10) and ID (N = 10) infants at 2 and 4 mo of age.

RESULTS

Reticulocyte hemoglobin (RET-He) was lower at 2 wk in the ID group. Significant IS compared with ID differences in serum iron indices were present at 2 mo, but Hgb and RBC indices differed only at 4 mo (P < 0.05). Serum and CSF metabolomic profiles of the ID and IS groups differed at 2 and 4 mo (P < 0.05). Key metabolites, including homostachydrine and stachydrine (4-5-fold lower at 4 mo in ID group, P < 0.05), were altered in both serum and CSF. Iron indices and RET-He at 2 mo, but not Hgb or other RBC indices, were correlated with altered CSF metabolic profile at 4 mo and had comparable predictive accuracy (area under the receiver operating characteristic curve scores, 0.75-0.80).

CONCLUSIONS

Preanemic ID at 2 mo was associated with metabolic alterations in serum and CSF in infant monkeys. Among the RBC indices, only RET-He predicted the future risk of abnormal CSF metabolic profile with a predictive accuracy comparable to serum iron indices. The concordance of homostachydrine and stachydrine changes in serum and CSF indicates their potential use as early biomarkers of brain metabolic dysfunction in infantile ID.

Diet-derived circulating antioxidants and risk of inflammatory bowel disease: a Mendelian randomization study and meta-analysis

Background

Previous studies have shown conflicting results regarding the impact of circulating antioxidants on the risk of inflammatory bowel disease (IBD). In this study, our intent was to investigate the causal relationship between circulating antioxidants and IBD using Mendelian randomization (MR).

Methods

Instrumental variables for absolute circulating antioxidants (ascorbate, retinol, lycopene, and β-carotene) and circulating antioxidant metabolites (α-tocopherol, γ-tocopherol, ascorbate, and retinol) were screened from published studies. We obtained outcome data from two genome-wide association study (GWAS) databases, including the international inflammatory bowel disease genetics consortium (IIBDGC, 14,927 controls and 5,956 cases for Crohn's disease (CD), 20,464 controls and 6,968 cases for ulcerative colitis (UC), and 21,770 controls and 12,882 cases for IBD) and the FinnGen study (375,445 controls and 1,665 cases for CD, 371,530 controls and 5,034 cases for UC, and 369,652 controls and 7,625 cases for IBD). MR analysis was performed in each of the two databases and those results were pooled using meta-analysis to assess the overall effect of exposure on each phenotype. In order to confirm the strength of the findings, we additionally conducted a replication analysis using the UK Biobank.

Results

In the meta-analysis of the IIBDGC and FinnGen, we found that each unit increase in absolute circulating level of retinol was associated with a 72% reduction in the risk of UC (OR: 0.28, 95% CI: 0.10 to 0.78, P=0.015). The UC GWAS data from the UK Biobank also confirmed this causal relationship (OR: 0.99, 95% CI: 0.97 to 1.00, P=0.016). In addition, there was suggestive evidence that absolute retinol level was negatively associated with IBD (OR: 0.41, 95% CI: 0.18 to 0.92, P=0.031). No other causal relationship was found.

Conclusion

Our results provide strong evidence that the absolute circulating level of retinol is associated with a reduction in the risk of UC. Further MR studies with more instrumental variables on circulating antioxidants, especially absolute circulating antioxidants, are needed to confirm our results.

Genetic influences on circulating retinol and its relationship to human health

Retinol is a fat-soluble vitamin that plays an essential role in many biological processes throughout the human lifespan. Here, we perform the largest genome-wide association study (GWAS) of retinol to date in up to 22,274 participants. We identify eight common variant loci associated with retinol, as well as a rare-variant signal. An integrative gene prioritisation pipeline supports novel retinol-associated genes outside of the main retinol transport complex (RBP4:TTR) related to lipid biology, energy homoeostasis, and endocrine signalling. Genetic proxies of circulating retinol were then used to estimate causal relationships with almost 20,000 clinical phenotypes via a phenome-wide Mendelian randomisation study (MR-pheWAS). The MR-pheWAS suggests that retinol may exert causal effects on inflammation, adiposity, ocular measures, the microbiome, and MRI-derived brain phenotypes, amongst several others. Conversely, circulating retinol may be causally influenced by factors including lipids and serum creatinine. Finally, we demonstrate how a retinol polygenic score could identify individuals more likely to fall outside of the normative range of circulating retinol for a given age. In summary, this study provides a comprehensive evaluation of the genetics of circulating retinol, as well as revealing traits which should be prioritised for further investigation with respect to retinol related therapies or nutritional intervention.

Role of Transporters and Enzymes in Metabolism and Distribution of 4-Chlorokynurenine (AV-101)

4-Chlorokynurenine (4-Cl-KYN, AV-101) is a prodrug of a NMDA receptor antagonist and is in clinical development for potential CNS indications. We sought to further understand the distribution and metabolism of 4-Cl-KYN, as this information might provide a strategy to enhance the clinical development of this drug. We used excretion studies in rats, in vitro transporter assays, and pharmacogenetic analysis of clinical trial data to determine how 4-Cl-KYN and metabolites are distributed. Our data indicated that a novel acetylated metabolite (N-acetyl-4-Cl-KYN) did not affect the uptake of 4-Cl-KYN across the blood-brain barrier via LAT1. 4-Cl-KYN and its metabolites were found to be renally excreted in rodents. In addition, we found that N-acetyl-4-Cl-KYN inhibited renal and hepatic transporters involved in excretion. Thus, this metabolite has the potential to limit the excretion of a range of compounds. Our pharmacogenetic analysis found that a SNP in N-acetyltransferase 8 (NAT8, rs13538) was linked to levels of N-acetyl-4-Cl-KYN relative to 4-Cl-KYN found in the plasma and that a SNP in SLC7A5 (rs28582913) was associated with the plasma levels of the active metabolite, 7-Cl-KYNA. Thus, we have a pharmacogenetics-based association for plasma drug level that could aid in the drug development of 4-Cl-KYN and have investigated the interaction of a novel metabolite with drug transporters.

The Relationship of Circulating Choline and Choline-Related Metabolite Levels with Health Outcomes: A Scoping Review of Genome-Wide Association Studies and Mendelian Randomization Studies

Choline is essential for proper liver, muscle, brain, lipid metabolism, cellular membrane composition, and repair. Understanding genetic determinants of circulating choline metabolites can help identify new determinants of choline metabolism, requirements, and their link to disease endpoints. We conducted a scoping review to identify studies assessing the association of genetic polymorphisms on circulating choline and choline-related metabolite concentrations and subsequent associations with health outcomes. This study follows the Preferred Reporting Items for Systematic Reviews and Meta-Analyses statement scoping review extension. Literature was searched to September 28, 2022, in 4 databases: Embase, MEDLINE, Web of Science, and the Biological Science Index. Studies of any duration in humans were considered. Any genome-wide association study (GWAS) investigating genetic variant associations with circulating choline and/or choline-related metabolites and any Mendelian randomization (MR) study investigating the association of genetically predicted circulating choline and/or choline-related metabolites with any health outcome were considered. Qualitative evidence is presented in summary tables. From 1248 total reviewed articles, 53 were included (GWAS = 27; MR = 26). Forty-two circulating choline-related metabolites were tested in association with genetic variants in GWAS studies, primarily trimethylamine N-oxide, betaine, sphingomyelins, lysophosphatidylcholines, and phosphatidylcholines. MR studies investigated associations between 52 total unique choline metabolites and 66 unique health outcomes. Of these, 47 significant associations were reported between 16 metabolites (primarily choline, lysophosphatidylcholines, phosphatidylcholines, betaine, and sphingomyelins) and 27 health outcomes including cancer, cardiovascular, metabolic, bone, and brain-related outcomes. Some articles reported significant associations between multiple choline types and the same health outcome. Genetically predicted circulating choline and choline-related metabolite concentrations are associated with a wide variety of health outcomes. Further research is needed to assess how genetic variability influences choline metabolism and whether individuals with lower genetically predicted circulating choline and choline-related metabolite concentrations would benefit from a dietary intervention or supplementation.

Individual bioenergetic capacity as a potential source of resilience to Alzheimer's disease

Impaired glucose uptake in the brain is one of the earliest presymptomatic manifestations of Alzheimer's disease (AD). The absence of symptoms for extended periods of time suggests that compensatory metabolic mechanisms can provide resilience. Here, we introduce the concept of a systemic 'bioenergetic capacity' as the innate ability to maintain energy homeostasis under pathological conditions, potentially serving as such a compensatory mechanism. We argue that fasting blood acylcarnitine profiles provide an approximate peripheral measure for this capacity that mirrors bioenergetic dysregulation in the brain. Using unsupervised subgroup identification, we show that fasting serum acylcarnitine profiles of participants from the AD Neuroimaging Initiative yields bioenergetically distinct subgroups with significant differences in AD biomarker profiles and cognitive function. To assess the potential clinical relevance of this finding, we examined factors that may offer diagnostic and therapeutic opportunities. First, we identified a genotype affecting the bioenergetic capacity which was linked to succinylcarnitine metabolism and significantly modulated the rate of future cognitive decline. Second, a potentially modifiable influence of beta-oxidation efficiency seemed to decelerate bioenergetic aging and disease progression. Our findings, which are supported by data from more than 9,000 individuals, suggest that interventions tailored to enhance energetic health and to slow bioenergetic aging could mitigate the risk of symptomatic AD, especially in individuals with specific mitochondrial genotypes.

Metabolic changes before and after weaning in Dezhou donkey foals in relation to gut microbiota

Weaning is undoubtedly one of the most crucial stages in the growth and development of all mammalian animals, including donkey foals. Weaning is a dynamic and coordinated process of the body, which is closely associated with the health, nutrition, and metabolism of the host. Many studies have shown that the intestinal microbiota and serum metabolites of mammals exhibit different changes during lactation, weaning, and postweaning. However, the alterations in serum metabolites in donkey foals before and postweaning and the correlation between serum metabolites and intestinal microbiota are largely unknown. This study is based on the fecal 16S rRNA and serum metabolomes of Dezhou donkey foals. In total, 10 samples (fecal and serum) were collected during the following three stages: before weaning (F.M.1), during weaning (F.M.3), and postweaning (F.M.6). To study the alterations in intestinal microflora, serum metabolites, and their correlation before and postweaning. We found that with the growth and weaning progress of donkey foals, the intestinal microbiota of donkey foals underwent obvious changes, and the diversity of fecal bacteria increased (Chao1 and Shannon indexes). The main intestinal microbial flora of donkey foals include Bacteroides and Firmicutes. We found many microbiota that are associated with immunity and digestion in the postweaning group, such as Verrucomicrobiales, Clostridia, Oscillospiraceae, Akkermansia, and Rikenellaceae, which can be considered microbial markers for the transition from liquid milk to solid pellet feed. Clostridia and Oscillospiraceae can produce organic acids, including butyric acid and acetic acid, which are crucial for regulating the intestinal microecological balance of donkeys. Furthermore, the metabolome showed that the serum metabolites enriched before and postweaning were mainly related to arachidonic acid metabolism and riboflavin metabolism. Riboflavin was associated with the development of the small intestine and affected the absorption of the small intestine. We also found that the changes in the gut microbiome of the foals were significantly correlated with changes in serum metabolites, including lysophosphatidylcholine (LPC; 12,0) and positively correlated with Lachnoclostridium and Roseburia. To summarize, this study provides theoretical data for the changes in the intestinal microbiome and serum metabolism during the entire weaning period of donkey foals.

A Genome-Wide Association Study of Serum Metabolite Profiles in Septic Shock Patients

OBJECTIVES

We sought to assess whether genetic associations with metabolite concentrations in septic shock patients could be used to identify pathways of potential importance for understanding sepsis pathophysiology.

DESIGN

Retrospective multicenter cohort studies of septic shock patients.

SETTING

All participants who were admitted to 27 participating hospital sites in three countries (Australia, New Zealand, and the United Kingdom) were eligible for inclusion.

PATIENTS

Adult, critically ill, mechanically ventilated patients with septic shock (n = 230) who were a subset of the Adjunctive Corticosteroid Treatment in Critically Ill Patients with Septic Shock trial (ClinicalTrials.gov number: NCT01448109).

INTERVENTIONS

None.

MEASUREMENTS AND MAIN RESULTS

A genome-wide association study was conducted for a range of serum metabolite levels for participants. Genome-wide significant associations (p ≤ 5 × 10-8) were found for the two major ketone bodies (3-hydroxybutyrate [rs2456680] and acetoacetate [rs2213037] and creatinine (rs6851961). One of these single-nucleotide polymorphisms (SNPs) (rs2213037) was located in the alcohol dehydrogenase cluster of genes, which code for enzymes related to the metabolism of acetoacetate and, therefore, presents a plausible association for this metabolite. None of the three SNPs showed strong associations with risk of sepsis, 28- or 90-day mortality, or Acute Physiology and Chronic Health Evaluation score (a measure of sepsis severity).

CONCLUSIONS

We suggest that the genetic associations with metabolites may reflect a starvation response rather than processes involved in sepsis pathophysiology. However, our results require further investigation and replication in both healthy and diseased cohorts including those of different ancestry.

Genetic liability for diet-derived circulating antioxidants, oxidative stress, and risk of osteoarthritis: a Mendelian randomization study

Background

Although well-documented, the causal relationships between diet-derived circulating antioxidants, oxidative stress, and osteoarthritis (OA) are equivocal. The objective of this study is to employ two-sample Mendelian randomization (MR) to investigate possible causal relationships among dietary-derived circulating antioxidants, oxidative stress damage indicators, and OA risk.

Methods

Single-nucleotide polymorphisms for diet-derived circulating antioxidants (ascorbate, β-carotene, lycopene, retinol, and α-and γ-tocopherol), assessed as absolute levels and metabolites, as well as oxidative stress injury biomarkers (GSH, GPX, CAT, SOD, albumin, and total bilirubin), were retrieved from the published data and were used as genetic instrumental variables. Summary statistics for gene-OA associations were obtained from publicly available and two relatively large-scale GWAS meta-analyses to date. The inverse-variance weighting method was utilized as the primary MR analysis. Moreover, multivariable MR was used to determine if mediators (BMI and smoking) causally mediated any connection. Furthermore, for each exposure, MR analyses were conducted per outcome database and then meta-analyzed.

Results

Genetically predicted absolute retinol level was causally associated with hip OA risk [odds ratios (ORs) = 0.40, 95% confidence interval (CI) = 0.24-0.68, FDR-corrected p = 0.042]. Moreover, genetically predicted albumin level was causally associated with total OA risk (OR = 0.80, 95% CI = 0.75-0.86, FDR-corrected p = 2.20E-11), as well as the risk of hip OA (OR = 0.75, 95% CI = 0.68-0.84, FDR-corrected p = 1.38E-06) and knee OA (OR = 0.82, 95% CI = 0.76-0.89, FDR-corrected p = 4.49E-06). In addition, MVMR confirmed that the effect of albumin on hip OA is independent of smoking initiation, alcoholic drinks per week, and moderate-to-vigorous physical activity levels but may be influenced by BMI.

Conclusion

Evidence from our study supports a potentially protective effect of high levels of retinol and albumin on OA risk.

Gene-metabolite annotation with shortest reactional distance enhances metabolite genome-wide association studies results

Metabolite genome-wide association studies (mGWAS) have advanced our understanding of the genetic control of metabolite levels. However, interpreting these associations remains challenging due to a lack of tools to annotate gene-metabolite pairs beyond the use of conservative statistical significance threshold. Here, we introduce the shortest reactional distance (SRD) metric, drawing from the comprehensive KEGG database, to enhance the biological interpretation of mGWAS results. We applied this approach to three independent mGWAS, including a case study on sickle cell disease patients. Our analysis reveals an enrichment of small SRD values in reported mGWAS pairs, with SRD values significantly correlating with mGWAS p values, even beyond the standard conservative thresholds. We demonstrate the utility of SRD annotation in identifying potential false negatives and inaccuracies within current metabolic pathway databases. Our findings highlight the SRD metric as an objective, quantitative and easy-to-compute annotation for gene-metabolite pairs, suitable to integrate statistical evidence to biological networks.

Systems biology and machine learning approaches identify metabolites that influence dietary lifespan and healthspan responses across flies and humans

Dietary restriction (DR) is a potent method to enhance lifespan and healthspan, but individual responses are influenced by genetic variations. Understanding how metabolism-related genetic differences impact longevity and healthspan are unclear. To investigate this, we used metabolites as markers to reveal how different genotypes respond to diet to influence longevity and healthspan traits. We analyzed data from Drosophila Genetic Reference Panel strains raised under AL and DR conditions, combining metabolomic, phenotypic, and genome-wide information. Employing two computational methods across species-random forest modeling within the DGRP and Mendelian randomization in the UK Biobank-we pinpointed key traits with cross-species relevance that influence lifespan and healthspan. Notably, orotate was linked to parental age at death in humans and counteracted DR effects in flies, while threonine extended lifespan, in a strain- and sex-specific manner. Thus, utilizing natural genetic variation data from flies and humans, we employed a systems biology approach to elucidate potential therapeutic pathways and metabolomic targets for diet-dependent changes in lifespan and healthspan.

DBS are suitable for 1,5-anhydroglucitol monitoring in GSD1b and G6PC3-deficient patients taking SGLT2 inhibitors to treat neutropenia

Glycogen storage disease type Ib (GSD1b) and G6PC3-deficiency are rare autosomal recessive diseases caused by inactivating mutations in SLC37A4 (coding for G6PT) and G6PC3, respectively. Both diseases are characterized by neutropenia and neutrophil dysfunction due to the intracellular accumulation of 1,5-anhydroglucitol-6-phosphate (1,5-AG6P), a potent inhibitor of hexokinases. We recently showed that the use of SGLT2 inhibitor therapy to reduce tubular reabsorption of its precursor, 1,5-anhydroglucitol (1,5-AG), a glucose analog present in blood, successfully restored the neutropenia and neutrophil function in G6PC3-deficient and GSD1b patients. The intra-individual variability of response to the treatment and the need to adjust the dose during treatment, especially in pediatric populations, can only be efficiently optimized if the concentration of 1,5-AG in blood is monitored during treatment, together with the patients' clinical signs and symptoms. Monitoring the 1,5-AG levels would be greatly simplified if it could be performed on dry blood spots (DBS) which are easy to collect, store and transport. The challenge is to know if a suitable method can be developed to perform accurate and reproducible assays for 1,5-AG using DBS. Here, we describe and validate an assay that quantifies 1,5-AG in DBS using isotopic dilution quantitation by LC-MS/MS that should greatly facilitate patients' follow-up. 1,5-AG levels measured in plasma and DBS give comparable values. This assay was used to monitor the levels of 1,5-AG in DBS from 3 G6PC3-deficient and 6 GSD1b patients during treatment with SGLT2 inhibitors. We recommend this approach to verify the adequate therapeutical response and compliance to the treatment in G6PC3-deficient and GSD1b patients treated with SGLT2 inhibitors.

Plasma metabolomics profiles in Black and White participants of the Adventist Health Study-2 cohort

BACKGROUND

Black Americans suffer disparities in risk for cardiometabolic and other chronic diseases. Findings from the Adventist Health Study-2 (AHS-2) cohort have shown associations of plant-based dietary patterns and healthy lifestyle factors with prevention of such diseases. Hence, it is likely that racial differences in metabolic profiles correlating with disparities in chronic diseases are explained largely by diet and lifestyle, besides social determinants of health.

METHODS

Untargeted plasma metabolomics screening was performed on plasma samples from 350 participants of the AHS-2, including 171 Black and 179 White participants, using ultrahigh-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) and a global platform of 892 metabolites. Differences in metabolites or biochemical subclasses by race were analyzed using linear regression, considering various models adjusted for known confounders, dietary and/or other lifestyle behaviors, social vulnerability, and psychosocial stress. The Storey permutation approach was used to adjust for false discovery at FDR < 0.05.

RESULTS

Linear regression revealed differential abundance of over 40% of individual metabolites or biochemical subclasses when comparing Black with White participants after adjustment for false discovery (FDR < 0.05), with the vast majority showing lower abundance in Blacks. Associations were not appreciably altered with adjustment for dietary patterns and socioeconomic or psychosocial stress. Metabolite subclasses showing consistently lower abundance in Black participants included various lipids, such as lysophospholipids, phosphatidylethanolamines, monoacylglycerols, diacylglycerols, and long-chain monounsaturated fatty acids, among other subclasses or lipid categories. Among all biochemical subclasses, creatine metabolism exclusively showed higher abundance in Black participants, although among metabolites within this subclass, only creatine showed differential abundance after adjustment for glomerular filtration rate. Notable metabolites in higher abundance in Black participants included methyl and propyl paraben sulfates, piperine metabolites, and a considerable proportion of acetylated amino acids, including many previously found associated with glomerular filtration rate.

CONCLUSIONS

Differences in metabolic profiles were evident when comparing Black and White participants of the AHS-2 cohort. These differences are likely attributed in part to dietary behaviors not adequately explained by dietary pattern covariates, besides other environmental or genetic factors. Alterations in these metabolites and associated subclasses may have implications for the prevention of chronic diseases in Black Americans.

Genome-wide association analysis of plasma lipidome identifies 495 genetic associations

The human plasma lipidome captures risk for cardiometabolic diseases. To discover new lipid-associated variants and understand the link between lipid species and cardiometabolic disorders, we perform univariate and multivariate genome-wide analyses of 179 lipid species in 7174 Finnish individuals. We fine-map the associated loci, prioritize genes, and examine their disease links in 377,277 FinnGen participants. We identify 495 genome-trait associations in 56 genetic loci including 8 novel loci, with a considerable boost provided by the multivariate analysis. For 26 loci, fine-mapping identifies variants with a high causal probability, including 14 coding variants indicating likely causal genes. A phenome-wide analysis across 953 disease endpoints reveals disease associations for 40 lipid loci. For 11 coronary artery disease risk variants, we detect strong associations with lipid species. Our study demonstrates the power of multivariate genetic analysis in correlated lipidomics data and reveals genetic links between diseases and lipid species beyond the standard lipids.

Metabolomic epidemiology offers insights into disease aetiology

Metabolomic epidemiology is the high-throughput study of the relationship between metabolites and health-related traits. This emerging and rapidly growing field has improved our understanding of disease aetiology and contributed to advances in precision medicine. As the field continues to develop, metabolomic epidemiology could lead to the discovery of diagnostic biomarkers predictive of disease risk, aiding in earlier disease detection and better prognosis. In this Review, we discuss key advances facilitated by the field of metabolomic epidemiology for a range of conditions, including cardiometabolic diseases, cancer, Alzheimer's disease and COVID-19, with a focus on potential clinical utility. Core principles in metabolomic epidemiology, including study design, causal inference methods and multi-omic integration, are briefly discussed. Future directions required for clinical translation of metabolomic epidemiology findings are summarized, emphasizing public health implications. Further work is needed to establish which metabolites reproducibly improve clinical risk prediction in diverse populations and are causally related to disease progression.

Quantitative trait loci mapping of circulating metabolites in cerebrospinal fluid to uncover biological mechanisms involved in brain-related phenotypes

Genomic studies of molecular traits have provided mechanistic insights into complex disease, though these lag behind for brain-related traits due to the inaccessibility of brain tissue. We leveraged cerebrospinal fluid (CSF) to study neurobiological mechanisms in vivo , measuring 5,543 CSF metabolites, the largest panel in CSF to date, in 977 individuals of European ancestry. Individuals originated from two separate cohorts including cognitively healthy subjects (n=490) and a well-characterized memory clinic sample, the Amsterdam Dementia Cohort (ADC, n=487). We performed metabolite quantitative trait loci (mQTL) mapping on CSF metabolomics and found 126 significant mQTLs, representing 65 unique CSF metabolites across 51 independent loci. To better understand the role of CSF mQTLs in brain-related disorders, we performed a metabolome-wide association study (MWAS), identifying 40 associations between CSF metabolites and brain traits. Similarly, over 90% of significant mQTLs demonstrated colocalized associations with brain-specific gene expression, unveiling potential neurobiological pathways.

Protein-metabolite association studies identify novel proteomic determinants of metabolite levels in human plasma

Although many novel gene-metabolite and gene-protein associations have been identified using high-throughput biochemical profiling, systematic studies that leverage human genetics to illuminate causal relationships between circulating proteins and metabolites are lacking. Here, we performed protein-metabolite association studies in 3,626 plasma samples from three human cohorts. We detected 171,800 significant protein-metabolite pairwise correlations between 1,265 proteins and 365 metabolites, including established relationships in metabolic and signaling pathways such as the protein thyroxine-binding globulin and the metabolite thyroxine, as well as thousands of new findings. In Mendelian randomization (MR) analyses, we identified putative causal protein-to-metabolite associations. We experimentally validated top MR associations in proof-of-concept plasma metabolomics studies in three murine knockout strains of key protein regulators. These analyses identified previously unrecognized associations between bioactive proteins and metabolites in human plasma. We provide publicly available data to be leveraged for studies in human metabolism and disease.

Relationship between food-derived antioxidant vitamin intake and breast cancer risk: a mendelian randomized study

BACKGROUND

In previous observational studies, food-derived antioxidant vitamins have been suggested to be associated with breast cancer. However, the findings were inconsistent and the causal relationship could not be clearly elucidated. To confirm the potential causal relationship between food-derived antioxidants (retinol, carotene, vitamin C and vitamin E) and the risk of breast cancer, we conducted a two-sample Mendelian randomization (MR) study.

METHODS

The instrumental variables (IVs) as proxies of genetic liability to food-derived antioxidant vitamins were obtained from the UK Biobank Database. We extracted breast cancer data (122,977 cases and 105,974 controls) from the Breast Cancer Consortium (BCAC). In addition, we studied estrogen expression status categorically, including estrogen receptor positive (ER+) breast cancer (69,501 cases and 105,974 controls) and versus estrogen receptor (ER-) negative breast cancer (21,468 cases and 105,974 controls). We performed two-sample Mendelian randomization study, and inverse variance-weighted (IVW) test was regarded as main analysis. Sensitivity analyses were further conducted to assess heterogeneity and horizontal pleiotropy.

RESULTS

The results of IVW showed that among the four food-derived antioxidants, only vitamin E had protective effect on the risk of overall breast cancer (OR = 0.837, 95% CI 0.757-0.926, P = 0.001) and ER+ breast cancer (OR = 0.823, 95% CI 0.693-0.977, P = 0.026). However, we found no association between food-derived vitamin E and ER- breast cancer.

CONCLUSIONS

Our study suggested food-derived vitamin E can decrease the risk of breast cancer overall and ER+ breast cancer, and the robustness of our results was confirmed by sensitivity analyses.

Genetic analysis of blood molecular phenotypes reveals common properties in the regulatory networks affecting complex traits

We evaluate the shared genetic regulation of mRNA molecules, proteins and metabolites derived from whole blood from 3029 human donors. We find abundant allelic heterogeneity, where multiple variants regulate a particular molecular phenotype, and pleiotropy, where a single variant associates with multiple molecular phenotypes over multiple genomic regions. The highest proportion of share genetic regulation is detected between gene expression and proteins (66.6%), with a further median shared genetic associations across 49 different tissues of 78.3% and 62.4% between plasma proteins and gene expression. We represent the genetic and molecular associations in networks including 2828 known GWAS variants, showing that GWAS variants are more often connected to gene expression in trans than other molecular phenotypes in the network. Our work provides a roadmap to understanding molecular networks and deriving the underlying mechanism of action of GWAS variants using different molecular phenotypes in an accessible tissue.

SGLT5 is the renal transporter for 1,5-anhydroglucitol, a major player in two rare forms of neutropenia

Neutropenia and neutrophil dysfunction in glycogen storage disease type 1b (GSD1b) and severe congenital neutropenia type 4 (SCN4), associated with deficiencies of the glucose-6-phosphate transporter (G6PT/SLC37A4) and the phosphatase G6PC3, respectively, are the result of the accumulation of 1,5-anhydroglucitol-6-phosphate in neutrophils. This is an inhibitor of hexokinase made from 1,5-anhydroglucitol (1,5-AG), an abundant polyol in blood. 1,5-AG is presumed to be reabsorbed in the kidney by a sodium-dependent-transporter of uncertain identity, possibly SGLT4/SLC5A9 or SGLT5/SLC5A10. Lowering blood 1,5-AG with an SGLT2-inhibitor greatly improved neutrophil counts and function in G6PC3-deficient and GSD1b patients. Yet, this effect is most likely mediated indirectly, through the inhibition of the renal 1,5-AG transporter by glucose, when its concentration rises in the renal tubule following inhibition of SGLT2. To identify the 1,5-AG transporter, both human and mouse SGLT4 and SGLT5 were expressed in HEK293T cells and transport measurements were performed with radiolabelled compounds. We found that SGLT5 is a better carrier for 1,5-AG than for mannose, while the opposite is true for human SGLT4. Heterozygous variants in SGLT5, associated with a low level of blood 1,5-AG in humans cause a 50-100% reduction in 1,5-AG transport activity tested in model cell lines, indicating that SGLT5 is the predominant kidney 1,5-AG transporter. These and other findings led to the conclusion that (1) SGLT5 is the main renal transporter of 1,5-AG; (2) frequent heterozygous mutations (allelic frequency > 1%) in SGLT5 lower blood 1,5-AG, favourably influencing neutropenia in G6PC3 or G6PT deficiency; (3) the effect of SGLT2-inhibitors on blood 1,5-AG level is largely indirect; (4) specific SGLT5-inhibitors would be more efficient to treat these neutropenias than SGLT2-inhibitors.

Integrative metabolomics science in Alzheimer's disease: Relevance and future perspectives

Alzheimer's disease (AD) is determined by various pathophysiological mechanisms starting 10-25 years before the onset of clinical symptoms. As multiple functionally interconnected molecular/cellular pathways appear disrupted in AD, the exploitation of high-throughput unbiased omics sciences is critical to elucidating the precise pathogenesis of AD. Among different omics, metabolomics is a fast-growing discipline allowing for the simultaneous detection and quantification of hundreds/thousands of perturbed metabolites in tissues or biofluids, reproducing the fluctuations of multiple networks affected by a disease. Here, we seek to critically depict the main metabolomics methodologies with the aim of identifying new potential AD biomarkers and further elucidating AD pathophysiological mechanisms. From a systems biology perspective, as metabolic alterations can occur before the development of clinical signs, metabolomics - coupled with existing accessible biomarkers used for AD screening and diagnosis - can support early disease diagnosis and help develop individualized treatment plans. Presently, the majority of metabolomic analyses emphasized that lipid metabolism is the most consistently altered pathway in AD pathogenesis. The possibility that metabolomics may reveal crucial steps in AD pathogenesis is undermined by the difficulty in discriminating between the causal or epiphenomenal or compensatory nature of metabolic findings.

A genome-wide association study identifies 41 loci associated with eicosanoid levels

Eicosanoids are biologically active derivatives of polyunsaturated fatty acids with broad relevance to health and disease. We report a genome-wide association study in 8406 participants of the Atherosclerosis Risk in Communities Study, identifying 41 loci associated with 92 eicosanoids and related metabolites. These findings highlight loci required for eicosanoid biosynthesis, including FADS1-3, ELOVL2, and numerous CYP450 loci. In addition, significant associations implicate a range of non-oxidative lipid metabolic processes in eicosanoid regulation, including at PKD2L1/SCD and several loci involved in fatty acyl-CoA metabolism. Further, our findings highlight select clearance mechanisms, for example, through the hepatic transporter encoded by SLCO1B1. Finally, we identify eicosanoids associated with aspirin and non-steroidal anti-inflammatory drug use and demonstrate the substantial impact of genetic variants even for medication-associated eicosanoids. These findings shed light on both known and unknown aspects of eicosanoid metabolism and motivate interest in several gene-eicosanoid associations as potential functional participants in human disease.

Ancestry-driven metabolite variation provides insights into disease states in admixed populations

BACKGROUND

Metabolic pathways are related to physiological functions and disease states and are influenced by genetic variation and environmental factors. Hispanics/Latino individuals have ancestry-derived genomic regions (local ancestry) from their recent admixture that have been less characterized for associations with metabolite abundance and disease risk.

METHODS

We performed admixture mapping of 640 circulating metabolites in 3887 Hispanic/Latino individuals from the Hispanic Community Health Study/Study of Latinos (HCHS/SOL). Metabolites were quantified in fasting serum through non-targeted mass spectrometry (MS) analysis using ultra-performance liquid chromatography-MS/MS. Replication was performed in 1856 nonoverlapping HCHS/SOL participants with metabolomic data.

RESULTS

By leveraging local ancestry, this study identified significant ancestry-enriched associations for 78 circulating metabolites at 484 independent regions, including 116 novel metabolite-genomic region associations that replicated in an independent sample. Among the main findings, we identified Native American enriched genomic regions at chromosomes 11 and 15, mapping to FADS1/FADS2 and LIPC, respectively, associated with reduced long-chain polyunsaturated fatty acid metabolites implicated in metabolic and inflammatory pathways. An African-derived genomic region at chromosome 2 was associated with N-acetylated amino acid metabolites. This region, mapped to ALMS1, is associated with chronic kidney disease, a disease that disproportionately burdens individuals of African descent.

CONCLUSIONS

Our findings provide important insights into differences in metabolite quantities related to ancestry in admixed populations including metabolites related to regulation of lipid polyunsaturated fatty acids and N-acetylated amino acids, which may have implications for common diseases in populations.

mGWAS-Explorer 2.0: Causal Analysis and Interpretation of Metabolite-Phenotype Associations

Metabolomics-based genome-wide association studies (mGWAS) are key to understanding the genetic regulations of metabolites in complex phenotypes. We previously developed mGWAS-Explorer 1.0 to link single-nucleotide polymorphisms (SNPs), metabolites, genes and phenotypes for hypothesis generation. It has become clear that identifying potential causal relationships between metabolites and phenotypes, as well as providing deep functional insights, are crucial for further downstream applications. Here, we introduce mGWAS-Explorer 2.0 to support the causal analysis between >4000 metabolites and various phenotypes. The results can be interpreted within the context of semantic triples and molecular quantitative trait loci (QTL) data. The underlying R package is released for reproducible analysis. Using two case studies, we demonstrate that mGWAS-Explorer 2.0 is able to detect potential causal relationships between arachidonic acid and Crohn's disease, as well as between glycine and coronary heart disease.

Diet-derived circulating antioxidants and risk of knee osteoarthritis, hip osteoarthritis and rheumatoid arthritis: a two-sample Mendelian randomization study

Objective

To examine the causal associations of diet-derived circulating antioxidants with knee osteoarthritis (OA), hip OA, and rheumatoid arthritis (RA) within the two-sample Mendelian randomization (MR) framework.

Method

Independent single-nucleotide polymorphisms (SNPs) significantly associated with circulating levels of diet-derived antioxidants (retinol, β-carotene, lycopene, vitamin C and vitamin E) were extracted as genetic instruments. Summary statistics of genetic instruments associated with knee OA, hip OA, and RA were obtained from corresponding genome-wide association studies (GWASs). The inverse-variance weighted (IVW) was applied as the primary analysis method, with four sensitivity analysis approaches employed to evaluate the robustness of the primary results.

Results

Genetically determined per unit increment of absolute circulating levels of retinol was significantly associated with a reduced risk of hip OA [odds ratio (OR) = 0.45, 95% confidence interval (CI) 0.26-0.78, p = 4.43 × 10^-3], while genetically determined per unit increase in absolute circulating levels of β-carotene was suggestively associated with increased risk of RA (OR = 1.32, 95% CI 1.07-1.62, p = 9.10 × 10^-3). No other causal association was found. Significant evidence for heterogeneity and pleiotropic outlier was only identified when absolute circulating vitamin C was evaluated as the exposure, whereas all sensitive analysis provided consistently non-significant results.

Conclusion

Our results demonstrated that genetically determined lifelong higher exposure to absolute circulating levels of retinol is associated with a decreased risk of hip OA. Further MR study with more genetic instruments for absolute circulating levels of antioxidants are needed to confirm our results.

Multi-omic longitudinal study reveals immune correlates of clinical course among hospitalized COVID-19 patients

The IMPACC cohort, composed of >1,000 hospitalized COVID-19 participants, contains five illness trajectory groups (TGs) during acute infection (first 28 days), ranging from milder (TG1-3) to more severe disease course (TG4) and death (TG5). Here, we report deep immunophenotyping, profiling of >15,000 longitudinal blood and nasal samples from 540 participants of the IMPACC cohort, using 14 distinct assays. These unbiased analyses identify cellular and molecular signatures present within 72 h of hospital admission that distinguish moderate from severe and fatal COVID-19 disease. Importantly, cellular and molecular states also distinguish participants with more severe disease that recover or stabilize within 28 days from those that progress to fatal outcomes (TG4 vs. TG5). Furthermore, our longitudinal design reveals that these biologic states display distinct temporal patterns associated with clinical outcomes. Characterizing host immune responses in relation to heterogeneity in disease course may inform clinical prognosis and opportunities for intervention.

Unique genetic architecture of CSF and brain metabolites pinpoints the novel targets for the traits of human wellness

Brain metabolism perturbation can contribute to traits and diseases. We conducted the first large-scale CSF and brain genome-wide association studies, which identified 219 independent associations (59.8% novel) for 144 CSF metabolites and 36 independent associations (55.6% novel) for 34 brain metabolites. Most of the novel signals (97.7% and 70.0% in CSF and brain) were tissue specific. We also integrated MWAS-FUSION approaches with Mendelian Randomization and colocalization to identify causal metabolites for 27 brain and human wellness phenotypes and identified eight metabolites to be causal for eight traits (11 relationships). Low mannose level was causal to bipolar disorder and as dietary supplement it may provide therapeutic benefits. Low galactosylglycerol level was found causal to Parkinson's Disease (PD). Our study expanded the knowledge of MQTL in central nervous system, provided insights into human wellness, and successfully demonstrates the utility of combined statistical approaches to inform interventions.

The full spectrum of OCT1 (SLC22A1) mutations bridges transporter biophysics to drug pharmacogenomics

Membrane transporters play a fundamental role in the tissue distribution of endogenous compounds and xenobiotics and are major determinants of efficacy and side effects profiles. Polymorphisms within these drug transporters result in inter-individual variation in drug response, with some patients not responding to the recommended dosage of drug whereas others experience catastrophic side effects. For example, variants within the major hepatic Human organic cation transporter OCT1 (SLC22A1) can change endogenous organic cations and many prescription drug levels. To understand how variants mechanistically impact drug uptake, we systematically study how all known and possible single missense and single amino acid deletion variants impact expression and substrate uptake of OCT1. We find that human variants primarily disrupt function via folding rather than substrate uptake. Our study revealed that the major determinants of folding reside in the first 300 amino acids, including the first 6 transmembrane domains and the extracellular domain (ECD) with a stabilizing and highly conserved stabilizing helical motif making key interactions between the ECD and transmembrane domains. Using the functional data combined with computational approaches, we determine and validate a structure-function model of OCT1s conformational ensemble without experimental structures. Using this model and molecular dynamic simulations of key mutants, we determine biophysical mechanisms for how specific human variants alter transport phenotypes. We identify differences in frequencies of reduced function alleles across populations with East Asians vs European populations having the lowest and highest frequency of reduced function variants, respectively. Mining human population databases reveals that reduced function alleles of OCT1 identified in this study associate significantly with high LDL cholesterol levels. Our general approach broadly applied could transform the landscape of precision medicine by producing a mechanistic basis for understanding the effects of human mutations on disease and drug response.

Genetic studies of paired metabolomes reveal enzymatic and transport processes at the interface of plasma and urine

The kidneys operate at the interface of plasma and urine by clearing molecular waste products while retaining valuable solutes. Genetic studies of paired plasma and urine metabolomes may identify underlying processes. We conducted genome-wide studies of 1,916 plasma and urine metabolites and detected 1,299 significant associations. Associations with 40% of implicated metabolites would have been missed by studying plasma alone. We detected urine-specific findings that provide information about metabolite reabsorption in the kidney, such as aquaporin (AQP)-7-mediated glycerol transport, and different metabolomic footprints of kidney-expressed proteins in plasma and urine that are consistent with their localization and function, including the transporters NaDC3 (SLC13A3) and ASBT (SLC10A2). Shared genetic determinants of 7,073 metabolite-disease combinations represent a resource to better understand metabolic diseases and revealed connections of dipeptidase 1 with circulating digestive enzymes and with hypertension. Extending genetic studies of the metabolome beyond plasma yields unique insights into processes at the interface of body compartments.

Whole-Genome Sequencing Analysis of Human Metabolome in Multi-Ethnic Populations

Circulating metabolite levels may reflect the state of the human organism in health and disease, however, the genetic architecture of metabolites is not fully understood. We have performed a whole-genome sequencing association analysis of both common and rare variants in up to 11,840 multi-ethnic participants from five studies with up to 1666 circulating metabolites. We have discovered 1985 novel variant-metabolite associations, and validated 761 locus-metabolite associations reported previously. Seventy-nine novel variant-metabolite associations have been replicated, including three genetic loci located on the X chromosome that have demonstrated its involvement in metabolic regulation. Gene-based analysis have provided further support for seven metabolite-replicated loci pairs and their biologically plausible genes. Among those novel replicated variant-metabolite pairs, follow-up analyses have revealed that 26 metabolites have colocalized with 21 tissues, seven metabolite-disease outcome associations have been putatively causal, and 7 metabolites might be regulated by plasma protein levels. Our results have depicted the genetic contribution to circulating metabolite levels, providing additional insights into understanding human disease.

Treatment of the Neutropenia Associated with GSD1b and G6PC3 Deficiency with SGLT2 Inhibitors

Glycogen storage disease type Ib (GSD1b) is due to a defect in the glucose-6-phosphate transporter (G6PT) of the endoplasmic reticulum, which is encoded by the SLC37A4 gene. This transporter allows the glucose-6-phosphate that is made in the cytosol to cross the endoplasmic reticulum (ER) membrane and be hydrolyzed by glucose-6-phosphatase (G6PC1), a membrane enzyme whose catalytic site faces the lumen of the ER. Logically, G6PT deficiency causes the same metabolic symptoms (hepatorenal glycogenosis, lactic acidosis, hypoglycemia) as deficiency in G6PC1 (GSD1a). Unlike GSD1a, GSD1b is accompanied by low neutrophil counts and impaired neutrophil function, which is also observed, independently of any metabolic problem, in G6PC3 deficiency. Neutrophil dysfunction is, in both diseases, due to the accumulation of 1,5-anhydroglucitol-6-phosphate (1,5-AG6P), a potent inhibitor of hexokinases, which is slowly formed in the cells from 1,5-anhydroglucitol (1,5-AG), a glucose analog that is normally present in blood. Healthy neutrophils prevent the accumulation of 1,5-AG6P due to its hydrolysis by G6PC3 following transport into the ER by G6PT. An understanding of this mechanism has led to a treatment aimed at lowering the concentration of 1,5-AG in blood by treating patients with inhibitors of SGLT2, which inhibits renal glucose reabsorption. The enhanced urinary excretion of glucose inhibits the 1,5-AG transporter, SGLT5, causing a substantial decrease in the concentration of this polyol in blood, an increase in neutrophil counts and function and a remarkable improvement in neutropenia-associated clinical signs and symptoms.