Integrative analysis of metabolite GWAS illuminates the molecular basis of pleiotropy and genetic correlation

Haplotype Analysis Reveals Pleiotropic Disease Associations in the HLA Region

The human leukocyte antigen (HLA) region plays an important role in human health through involvement in immune cell recognition and maturation. While genetic variation in the HLA region is associated with many diseases, the pleiotropic patterns of these associations have not been systematically investigated. Here, we developed a haplotype approach to investigate disease associations phenome-wide for 412,181 Finnish individuals and 2,459 traits. Across the 1,035 diseases with a GWAS association, we found a 17-fold average per-SNP enrichment of hits in the HLA region. Altogether, we identified 7,649 HLA associations across 647 traits, including 1,750 associations uncovered by haplotype analysis. We find some haplotypes show trade-offs between diseases, while others consistently increase risk across traits, indicating a complex pleiotropic landscape involving a range of diseases. This study highlights the extensive impact of HLA variation on disease risk, and underscores the importance of classical and non-classical genes, as well as non-coding variation.

Adults prenatally exposed to the Dutch Famine exhibit a metabolic signature associated with a broad spectrum of common diseases

BACKGROUND

Exposure to famine in the prenatal period is associated with an increased risk of metabolic disease, including obesity and type 2 diabetes. We employed nuclear magnetic resonance (NMR) metabolomic profiling to identify the metabolic changes that are associated with survival of prenatal famine exposure during the Dutch Famine at the end of World War II and subsequently assess their link to disease.

METHODS

NMR metabolomics data were generated from serum in 480 individuals prenatally exposed to famine (mean 58.8 years, 0.5 SD) and 464 controls (mean 57.9 years, 5.4 SD). We tested associations of prenatal famine exposure with levels of 168 individual metabolic biomarkers and compared the metabolic biomarker signature of famine exposure with those of 154 common diseases.

RESULTS

Prenatal famine exposure was associated with higher concentrations of branched-chain amino acids ((iso)-leucine), aromatic amino acid (tyrosine), and glucose in later life (0.2-0.3 SD, p < 3 × 10-3). The metabolic biomarker signature of prenatal famine exposure was positively correlated to that of incident type 2 diabetes from the UK Biobank (r = 0.77, p = 3 × 10-27), also when re-estimating the signature of prenatal famine exposure among individuals without diabetes (r = 0.67, p = 1 × 10-18). Remarkably, this association extended to 115 common diseases for which signatures were available (0.3 ≤ r ≤ 0.9, p < 3.2 × 10-4). Correlations among metabolic signatures of famine exposure and disease outcomes were attenuated when the famine signature was adjusted for body mass index.

CONCLUSIONS

Prenatal famine exposure is associated with a metabolic biomarker signature that strongly resembles signatures of a diverse set of diseases, an observation that can in part be attributed to a shared involvement of obesity.

MRBEE: A bias-corrected multivariable Mendelian randomization method

Mendelian randomization (MR) is an instrumental variable approach used to infer causal relationships between exposures and outcomes, which is becoming increasingly popular because of its ability to handle summary statistics from genome-wide association studies. However, existing MR approaches often suffer the bias from weak instrumental variables, horizontal pleiotropy and sample overlap. We introduce MRBEE (MR using bias-corrected estimating equation), a multivariable MR method capable of simultaneously removing weak instrument and sample overlap bias and identifying horizontal pleiotropy. Our extensive simulations and real data analyses reveal that MRBEE provides nearly unbiased estimates of causal effects, well-controlled type I error rates and higher power than comparably robust methods and is computationally efficient. Our real data analyses result in consistent causal effect estimates and offer valuable guidance for conducting multivariable MR studies, elucidating the roles of pleiotropy, and identifying total 42 horizontal pleiotropic loci missed previously that are associated with myopia, schizophrenia, and coronary artery disease.

Associations of Circulating Biomarkers with Disease Risks: A Two-Sample Mendelian Randomization Study

Circulating biomarkers play a pivotal role in personalized medicine, offering potential for disease screening, prevention, and treatment. Despite established associations between numerous biomarkers and diseases, elucidating their causal relationships is challenging. Mendelian Randomization (MR) can address this issue by employing genetic instruments to discern causal links. Additionally, using multiple MR methods with overlapping results enhances the reliability of discovered relationships. Here, we report an MR study using multiple methods, including inverse variance weighted, simple mode, weighted mode, weighted median, and MR-Egger. We use the MR-base resource (v0.5.6) from Hemani et al. 2018 to evaluate causal relationships between 212 circulating biomarkers (curated from UK Biobank analyses by Neale lab and from Shin et al. 2014, Roederer et al. 2015, and Kettunen et al. 2016 and 99 complex diseases (curated from several consortia by MRC IEU and Biobank Japan). We report novel causal relationships found by four or more MR methods between glucose and bipolar disorder (Mean Effect Size estimate across methods: 0.39) and between cystatin C and bipolar disorder (Mean Effect Size: -0.31). Based on agreement in four or more methods, we also identify previously known links between urate with gout and creatine with chronic kidney disease, as well as biomarkers that may be causal of cardiovascular conditions: apolipoprotein B, cholesterol, LDL, lipoprotein A, and triglycerides in coronary heart disease, as well as lipoprotein A, LDL, cholesterol, and apolipoprotein B in myocardial infarction. This Mendelian Randomization study not only corroborates known causal relationships between circulating biomarkers and diseases but also uncovers two novel biomarkers associated with bipolar disorder that warrant further investigation. Our findings provide insight into understanding how biological processes reflecting circulating biomarkers and their associated effects may contribute to disease etiology, which can eventually help improve precision diagnostics and intervention.

Associations of Circulating Biomarkers with Disease Risks: a Two-Sample Mendelian Randomization Study

Background

Circulating biomarkers play a pivotal role in personalized medicine, offering potential for disease screening, prevention, and treatment. Despite established associations between numerous biomarkers and diseases, elucidating their causal relationships is challenging. Mendelian Randomization (MR) can address this issue by employing genetic instruments to discern causal links. Additionally, using multiple MR methods with overlapping results enhances the reliability of discovered relationships.

Methods

Here we report an MR study using multiple methods, including inverse variance weighted, simple mode, weighted mode, weighted median, and MR Egger. We use the MR-base resource (v0.5.6)1 to evaluate causal relationships between 212 circulating biomarkers (curated from UK Biobank analyses by Neale lab and from Shin et al. 2014, Roederer et al. 2015, and Kettunen et al. 2016)2-4 and 99 complex diseases (curated from several consortia by MRC IEU and Biobank Japan).

Results

We report novel causal relationships found by 4 or more MR methods between glucose and bipolar disorder (Mean Effect Size estimate across methods: 0.39) and between cystatin C and bipolar disorder (Mean Effect Size: -0.31). Based on agreement in 4 or more methods, we also identify previously known links between urate with gout and creatine with chronic kidney disease, as well as biomarkers that may be causal of cardiovascular conditions: apolipoprotein B, cholesterol, LDL, lipoprotein A, and triglycerides in coronary heart disease, as well as lipoprotein A, LDL, cholesterol, and apolipoprotein B in myocardial infarction.

Conclusions

This Mendelian Randomization study not only corroborates known causal relationships between circulating biomarkers and diseases but also uncovers two novel biomarkers associated with bipolar disorder that warrant further investigation. Our findings provide insight into understanding how biological processes reflecting circulating biomarkers and their associated effects may contribute to disease etiology, which can eventually help improve precision diagnostics and intervention.

Epistasis and pleiotropy-induced variation for plant breeding

Epistasis refers to nonallelic interaction between genes that cause bias in estimates of genetic parameters for a phenotype with interactions of two or more genes affecting the same trait. Partitioning of epistatic effects allows true estimation of the genetic parameters affecting phenotypes. Multigenic variation plays a central role in the evolution of complex characteristics, among which pleiotropy, where a single gene affects several phenotypic characters, has a large influence. While pleiotropic interactions provide functional specificity, they increase the challenge of gene discovery and functional analysis. Overcoming pleiotropy-based phenotypic trade-offs offers potential for assisting breeding for complex traits. Modelling higher order nonallelic epistatic interaction, pleiotropy and non-pleiotropy-induced variation, and genotype × environment interaction in genomic selection may provide new paths to increase the productivity and stress tolerance for next generation of crop cultivars. Advances in statistical models, software and algorithm developments, and genomic research have facilitated dissecting the nature and extent of pleiotropy and epistasis. We overview emerging approaches to exploit positive (and avoid negative) epistatic and pleiotropic interactions in a plant breeding context, including developing avenues of artificial intelligence, novel exploitation of large-scale genomics and phenomics data, and involvement of genes with minor effects to analyse epistatic interactions and pleiotropic quantitative trait loci, including missing heritability.

Interplay Between Plasma Glycine and Branched-Chain Amino Acids Contributes to the Development of Hypertension and Coronary Heart Disease

BACKGROUND

Higher levels of plasma glycine are linked to a reduced risk, while increased levels of total branched-chain amino acids (tBCAAs) are associated with a higher risk of essential hypertension and coronary heart disease (CHD). As these metabolic components are interconnected, analyzing the tBCAAs/glycine ratio may help to understand their interplay in the pathogenesis of cardiovascular disease.

METHODS

The Cox regression approach was combined with the development of novel genetic tools for assessments of associations between plasma metabolomic data (glycine, tBCAAs, and tBCAAs/glycine ratio) from the UK Biobank and the development of hypertension and CHD. Genome-wide association study was performed on 186 523 White UK Biobank participants to identify new independent genetic instruments for the 2-sample Mendelian randomization analyses. P-gain statistic >10 identified instruments associated with tBCAAs/glycine ratio significantly stronger compared with individual amino acids. Outcomes of genome-wide association study on hypertension and CHD were derived from the UK Biobank (nonoverlapping sample), FinnGen, and CARDIoGRAMplusC4D.

RESULTS

The tBCAAs/glycine ratio was prospectively associated with a higher risk of developing hypertension and CHD (hazard ratio quintile Q5 versus Q1, 1.196 [95% CI, 1.109-1.289] and 1.226 [95% CI, 1.160-1.296], respectively). Mendelian randomization analysis demonstrated that tBCAAs/glycine ratio (P-gain >10) was a risk factor for hypertension (meta-analyzed inverse-variance weighted causal estimate 0.45 log odds ratio/SD (95% CI, 0.26-0.64) and CHD (0.48 [95% CI, 0.29-0.67]) with an absolute effect significantly larger compared with the effect of glycine (-0.06 [95% CI, -0.1 to -0.03] and -0.08 [95% CI, -0.11 to -0.05], respectively) or tBCAAs (0.22 [95% CI, 0.09-0.34] and 0.12 [95% CI, 0.01-0.24], respectively).

CONCLUSIONS

The total BCAAs/glycine ratio is a key element of the metabolic signature contributing to hypertension and CHD, which may reflect biological pathways shared by glycine and tBCAAs.

MANOCCA: a robust and computationally efficient test of covariance in high-dimension multivariate omics data

Multivariate analysis is becoming central in studies investigating high-throughput molecular data, yet, some important features of these data are seldom explored. Here, we present MANOCCA (Multivariate Analysis of Conditional CovAriance), a powerful method to test for the effect of a predictor on the covariance matrix of a multivariate outcome. The proposed test is by construction orthogonal to tests based on the mean and variance and is able to capture effects that are missed by both approaches. We first compare the performances of MANOCCA with existing correlation-based methods and show that MANOCCA is the only test correctly calibrated in simulation mimicking omics data. We then investigate the impact of reducing the dimensionality of the data using principal component analysis when the sample size is smaller than the number of pairwise covariance terms analysed. We show that, in many realistic scenarios, the maximum power can be achieved with a limited number of components. Finally, we apply MANOCCA to 1000 healthy individuals from the Milieu Interieur cohort, to assess the effect of health, lifestyle and genetic factors on the covariance of two sets of phenotypes, blood biomarkers and flow cytometry-based immune phenotypes. Our analyses identify significant associations between multiple factors and the covariance of both omics data.

Adults prenatally exposed to the Dutch Famine exhibit a metabolic signature associated with a broad spectrum of common diseases

Background

Exposure to famine in the prenatal period is associated with an increased risk of metabolic disease, including obesity and type-2 diabetes. We employed nuclear magnetic resonance (NMR) metabolomic profiling to provide a deeper insight into the metabolic changes associated with survival of prenatal famine exposure during the Dutch Famine at the end of World War II and explore their link to disease.

Methods

NMR metabolomics data were generated from serum in 480 individuals prenatally exposed to famine (mean 58.8 years, 0.5 SD) and 464 controls (mean 57.9 years, 5.4 SD). We tested associations of prenatal famine exposure with levels of 168 individual metabolic biomarkers and compared the metabolic biomarker signature of famine exposure with those of 154 common diseases.

Results

Prenatal famine exposure was associated with higher concentrations of branched-chain amino acids ((iso)-leucine), aromatic amino acid (tyrosine), and glucose in later life (0.2-0.3 SD, p < 3x10-3). The metabolic biomarker signature of prenatal famine exposure was positively correlated to that of incident type-2 diabetes (r = 0.77, p = 3x10-27), also when re-estimating the signature of prenatal famine exposure among individuals without diabetes (r = 0.67, p = 1x10-18). Remarkably, this association extended to 115 common diseases for which signatures were available (0.3 ≤ r ≤ 0.9, p < 3.2x10-4). Correlations among metabolic signatures of famine exposure and disease outcomes were attenuated when the famine signature was adjusted for body mass index.

Conclusions

Prenatal famine exposure is associated with a metabolic biomarker signature that strongly resembles signatures of a diverse set of diseases, an observation that can in part be attributed to a shared involvement of obesity.

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.

Uncovering genetic associations in the human diseasome using an endophenotype-augmented disease network

MOTIVATION

Many diseases, particularly cardiometabolic disorders, exhibit complex multimorbidities with one another. An intuitive way to model the connections between phenotypes is with a disease-disease network (DDN), where nodes represent diseases and edges represent associations, such as shared single-nucleotide polymorphisms (SNPs), between pairs of diseases. To gain further genetic understanding of molecular contributors to disease associations, we propose a novel version of the shared-SNP DDN (ssDDN), denoted as ssDDN+, which includes connections between diseases derived from genetic correlations with intermediate endophenotypes. We hypothesize that a ssDDN+ can provide complementary information to the disease connections in a ssDDN, yielding insight into the role of clinical laboratory measurements in disease interactions.

RESULTS

Using PheWAS summary statistics from the UK Biobank, we constructed a ssDDN+ revealing hundreds of genetic correlations between diseases and quantitative traits. Our augmented network uncovers genetic associations across different disease categories, connects relevant cardiometabolic diseases, and highlights specific biomarkers that are associated with cross-phenotype associations. Out of the 31 clinical measurements under consideration, HDL-C connects the greatest number of diseases and is strongly associated with both type 2 diabetes and heart failure. Triglycerides, another blood lipid with known genetic causes in non-mendelian diseases, also adds a substantial number of edges to the ssDDN. This work demonstrates how association with clinical biomarkers can better explain the shared genetics between cardiometabolic disorders. Our study can facilitate future network-based investigations of cross-phenotype associations involving pleiotropy and genetic heterogeneity, potentially uncovering sources of missing heritability in multimorbidities.

AVAILABILITY AND IMPLEMENTATION

The generated ssDDN+ can be explored at https://hdpm.biomedinfolab.com/ddn/biomarkerDDN.

Metabolomic associations of impaired awareness of hypoglycaemia in type 1 diabetes

This study investigates impaired awareness of hypoglycaemia (IAH), a complication of insulin therapy affecting 20-40% of individuals with type 1 diabetes. The exact pathophysiology is unclear, therefore we sought to identify metabolic signatures in IAH to elucidate potential pathophysiological pathways. Plasma samples from 578 individuals of the Dutch type 1 diabetes biomarker cohort, 67 with IAH and 108 without IAH (NAH) were analysed using the targeted metabolomics Biocrates AbsoluteIDQ p180 assay. Eleven metabolites were significantly associated with IAH. Genome-wide association studies of these 11 metabolites identified significant single nucleotide polymorphisms (SNPs) in C22:1-OH and phosphatidylcholine diacyl C36:6. After adjusting for the SNPs, 11 sphingomyelins and phosphatidylcholines were significantly higher in the IAH group in comparison to NAH. These metabolites are important components of the cell membrane and have been implicated to play a role in cell signalling in diabetes. These findings demonstrate the potential role of phosphatidylcholine and sphingomyelins in IAH.

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.

Assessing the causal effect of genetically predicted metabolites and metabolic pathways on stroke

BACKGROUND

Stroke is a common neurological disorder that disproportionately affects middle-aged and elderly individuals, leading to significant disability and mortality. Recently, human blood metabolites have been discovered to be useful in unraveling the underlying biological mechanisms of neurological disorders. Therefore, we aimed to evaluate the causal relationship between human blood metabolites and susceptibility to stroke.

METHODS

Summary data from genome-wide association studies (GWASs) of serum metabolites and stroke and its subtypes were obtained separately. A total of 486 serum metabolites were used as the exposure. Simultaneously, 11 different stroke phenotypes were set as the outcomes, including any stroke (AS), any ischemic stroke (AIS), large artery stroke (LAS), cardioembolic stroke (CES), small vessel stroke (SVS), lacunar stroke (LS), white matter hyperintensities (WMH), intracerebral hemorrhage (ICH), subarachnoid hemorrhage (SAH), transient ischemic attack (TIA), and brain microbleeds (BMB). A two-sample Mendelian randomization (MR) study was conducted to investigate the causal effects of serum metabolites on stroke and its subtypes. The inverse variance-weighted MR analyses were conducted as causal estimates, accompanied by a series of sensitivity analyses to evaluate the robustness of the results. Furthermore, a reverse MR analysis was conducted to assess the potential for reverse causation. Additionally, metabolic pathway analysis was performed using the web-based MetOrigin.

RESULTS

After correcting for the false discovery rate (FDR), MR analysis results revealed remarkable causative associations with 25 metabolites. Further sensitivity analyses confirmed that only four causative associations involving three specific metabolites passed all sensitivity tests, namely ADpSGEGDFXAEGGGVR* for AS (OR: 1.599, 95% CI 1.283-1.993, p = 2.92 × 10-5) and AIS (OR: 1.776, 95% CI 1.380-2.285, p = 8.05 × 10-6), 1-linoleoylglycerophosph-oethanolamine* for LAS (OR: 0.198, 95% CI 0.091-0.428, p = 3.92 × 10-5), and gamma-glutamylmethionine* for SAH (OR: 3.251, 95% CI 1.876-5.635, p = 2.66 × 10-5), thereby demonstrating a high degree of stability. Moreover, eight causative associations involving seven other metabolites passed both sensitivity tests and were considered robust. The association result of one metabolite (glutamate for LAS) was considered non-robust. As for the remaining metabolites, we speculate that they may potentially possess underlying causal relationships. Notably, no common metabolites emerged from the reverse MR analysis. Moreover, after FDR correction, metabolic pathway analysis identified 40 significant pathways across 11 stroke phenotypes.

CONCLUSIONS

The identified metabolites and their associated metabolic pathways are promising circulating metabolic biomarkers, holding potential for their application in stroke screening and preventive strategies within clinical settings.

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.

Uncovering genetic associations in the human diseasome using an endophenotype-augmented disease network

Many diseases exhibit complex multimorbidities with one another. An intuitive way to model the connections between phenotypes is with a disease-disease network (DDN), where nodes represent diseases and edges represent associations, such as shared single-nucleotide polymorphisms (SNPs), between pairs of diseases. To gain further genetic understanding of molecular contributors to disease associations, we propose a novel version of the shared-SNP DDN (ssDDN), denoted as ssDDN+, which includes connections between diseases derived from genetic correlations with endophenotypes. We hypothesize that a ssDDN+ can provide complementary information to the disease connections in a ssDDN, yielding insight into the role of clinical laboratory measurements in disease interactions. Using PheWAS summary statistics from the UK Biobank, we constructed a ssDDN+ revealing hundreds of genetic correlations between disease phenotypes and quantitative traits. Our augmented network uncovers genetic associations across different disease categories, connects relevant cardiometabolic diseases, and highlights specific biomarkers that are associated with cross-phenotype associations. Out of the 31 clinical measurements under consideration, HDL-C connects the greatest number of diseases and is strongly associated with both type 2 diabetes and diabetic retinopathy. Triglycerides, another blood lipid with known genetics causes in non-mendelian diseases, also adds a substantial number of edges to the ssDDN. Our study can facilitate future network-based investigations of cross-phenotype associations involving pleiotropy and genetic heterogeneity, potentially uncovering sources of missing heritability in multimorbidities.

GWAS and genetic and phenotypic correlations of plasma metabolites with complete blood count traits in healthy young pigs reveal implications for pig immune response

Introduction: In this study estimated genetic and phenotypic correlations between fifteen complete blood count (CBC) traits and thirty-three heritable plasma metabolites in young healthy nursery pigs. In addition, it provided an opportunity to identify candidate genes associated with variation in metabolite concentration and their potential association with immune response, disease resilience, and production traits.

Methods: The blood samples were collected from healthy young pigs and Nuclear Magnetic Resonance (NMR) was used to quantify plasma metabolites. CBC was determined using the ADVIA® 2120i Hematology System. Genetic correlations of metabolite with CBC traits and single step genome-wide association study (ssGWAS) were estimated using the BLUPF90 programs.

Results: Results showed low phenotypic correlation estimates between plasma metabolites and CBC traits. The highest phenotypic correlation was observed between lactic acid and plasma basophil concentration (0.36 ± 0.04; p &lt; 0.05). Several significant genetic correlations were found between metabolites and CBC traits. The plasma concentration of proline was genetically positively correlated with hemoglobin concentration (0.94 ± 0.03; p &lt; 0.05) and L-tyrosine was negatively correlated with mean corpuscular hemoglobin (MCH; −0.92 ± 0.74; p &lt; 0.05). The genomic regions identified in this study only explained a small percentage of the genetic variance of metabolites levels that were genetically correlated with CBC, resilience, and production traits.

Discussion: The results of this systems approach suggest that several plasma metabolite phenotypes are phenotypically and genetically correlated with CBC traits, suggesting that they may be potential genetic indicators of immune response following disease challenge. Genomic analysis revealed genes and pathways that might interact to modulate CBC, resilience, and production traits.

Exploiting the mediating role of the metabolome to unravel transcript-to-phenotype associations

Despite the success of genome-wide association studies (GWASs) in identifying genetic variants associated with complex traits, understanding the mechanisms behind these statistical associations remains challenging. Several methods that integrate methylation, gene expression, and protein quantitative trait loci (QTLs) with GWAS data to determine their causal role in the path from genotype to phenotype have been proposed. Here, we developed and applied a multi-omics Mendelian randomization (MR) framework to study how metabolites mediate the effect of gene expression on complex traits. We identified 216 transcript-metabolite-trait causal triplets involving 26 medically relevant phenotypes. Among these associations, 58% were missed by classical transcriptome-wide MR, which only uses gene expression and GWAS data. This allowed the identification of biologically relevant pathways, such as between ANKH and calcium levels mediated by citrate levels and SLC6A12 and serum creatinine through modulation of the levels of the renal osmolyte betaine. We show that the signals missed by transcriptome-wide MR are found, thanks to the increase in power conferred by integrating multiple omics layer. Simulation analyses show that with larger molecular QTL studies and in case of mediated effects, our multi-omics MR framework outperforms classical MR approaches designed to detect causal relationships between single molecular traits and complex phenotypes.

Effects of protein-coding variants on blood metabolite measurements and clinical biomarkers in the UK Biobank

Genome-wide association studies (GWASs) have established the contribution of common and low-frequency variants to metabolic blood measurements in the UK Biobank (UKB). To complement existing GWAS findings, we assessed the contribution of rare protein-coding variants in relation to 355 metabolic blood measurements-including 325 predominantly lipid-related nuclear magnetic resonance (NMR)-derived blood metabolite measurements (Nightingale Health Plc) and 30 clinical blood biomarkers-using 412,393 exome sequences from four genetically diverse ancestries in the UKB. Gene-level collapsing analyses were conducted to evaluate a diverse range of rare-variant architectures for the metabolic blood measurements. Altogether, we identified significant associations (p < 1 × 10^-8) for 205 distinct genes that involved 1,968 significant relationships for the Nightingale blood metabolite measurements and 331 for the clinical blood biomarkers. These include associations for rare non-synonymous variants in PLIN1 and CREB3L3 with lipid metabolite measurements and SYT7 with creatinine, among others, which may not only provide insights into novel biology but also deepen our understanding of established disease mechanisms. Of the study-wide significant clinical biomarker associations, 40% were not previously detected on analyzing coding variants in a GWAS in the same cohort, reinforcing the importance of studying rare variation to fully understand the genetic architecture of metabolic blood measurements.

Atlas of plasma NMR biomarkers for health and disease in 118,461 individuals from the UK Biobank

Blood lipids and metabolites are markers of current health and future disease risk. Here, we describe plasma nuclear magnetic resonance (NMR) biomarker data for 118,461 participants in the UK Biobank. The biomarkers cover 249 measures of lipoprotein lipids, fatty acids, and small molecules such as amino acids, ketones, and glycolysis metabolites. We provide an atlas of associations of these biomarkers to prevalence, incidence, and mortality of over 700 common diseases ( nightingalehealth.com/atlas ). The results reveal a plethora of biomarker associations, including susceptibility to infectious diseases and risk of various cancers, joint disorders, and mental health outcomes, indicating that abundant circulating lipids and metabolites are risk markers beyond cardiometabolic diseases. Clustering analyses indicate similar biomarker association patterns across different disease types, suggesting latent systemic connectivity in the susceptibility to a diverse set of diseases. This work highlights the value of NMR based metabolic biomarker profiling in large biobanks for public health research and translation.

Regulation of Human Endogenous Metabolites by Drug Transporters and Drug Metabolizing Enzymes: An Analysis of Targeted SNP-Metabolite Associations

Drug transporters and drug-metabolizing enzymes are primarily known for their role in the absorption, distribution, metabolism, and excretion (ADME) of small molecule drugs, but they also play a key role in handling endogenous metabolites. Recent cross-tissue co-expression network analyses have revealed a "Remote Sensing and Signaling Network" of multispecific, oligo-specific, and monospecific transporters and enzymes involved in endogenous metabolism. This includes many proteins from families involved in ADME (e.g., SLC22, SLCO, ABCC, CYP, UGT). Focusing on the gut-liver-kidney axis, we identified the endogenous metabolites potentially regulated by this network of ~1000 proteins by associating SNPs in these genes with the circulating levels of thousands of small, polar, bioactive metabolites, including free fatty acids, eicosanoids, bile acids, and other signaling metabolites that act in part via G-protein coupled receptors (GPCRs), nuclear receptors, and kinases. We identified 77 genomic loci associated with 7236 unique metabolites. This included metabolites that were associated with multiple, distinct loci, indicating coordinated regulation between multiple genes (including drug transporters and drug-metabolizing enzymes) of specific metabolites. We analyzed existing pharmacogenomic data and noted SNPs implicated in endogenous metabolite handling (e.g., rs4149056 in SLCO1B1) also affecting drug ADME. The overall results support the existence of close relationships, via interactions with signaling metabolites, between drug transporters and drug-metabolizing enzymes that are part of the Remote Sensing and Signaling Network, and with GPCRs and nuclear receptors. These analyses highlight the potential for drug-metabolite interactions at the interfaces of the Remote Sensing and Signaling Network and the ADME protein network.