Whole Genome Association Study of the Plasma Metabolome Identifies Metabolites Linked to Cardiometabolic Disease in Black Individuals

Genome-wide association study reveals shared and distinct genetic architecture underlying fatty acid and bioactive oxylipin metabolites in the Hispanic Community Health Study/Study of Latinos (HCHS/SOL)

Bioactive fatty acid-derived oxylipin molecules play key roles in mediating inflammation and oxidative stress, which underlie many chronic diseases. Circulating levels of fatty acids and oxylipins are influenced by both environmental and genetic factors; characterizing the genetic architecture of bioactive lipids could yield new insights into underlying biological pathways. Thus, we performed a genome wide association study (GWAS) of n=81 fatty acids and oxylipins in n=11,584 Hispanic Community Health Study/Study of Latinos (HCHS/SOL) participants with genetic and lipidomic data measured at study baseline (58.6% female, mean age = 46.1 years, standard deviation = 13.8 years). Additionally, given the effects of central obesity on inflammation, we examined interactions with waist circumference using two-degree-of-freedom joint tests. Heritability estimates ranged from 0% to 47.9%, and 48 of the 81oxylipins and fatty acids were significantly heritable. Moreover, 40 (49.4%) of the 81 oxylipins and fatty acids had at least one genome-wide significant (p< 6.94E-11) variant resulting in 19 independent genetic loci involved in fatty acid and oxylipin synthesis, as well as downstream pathways. Four loci (lead variant minor allele frequency [MAF] range: 0.08-0.50), including the desaturase-encoding FADS and the OATP1B1 transporter protein-encoding SLCO1B1, exhibited associations with four or more fatty acids and oxylipins. The majority of the 15 remaining loci (87.5%) (lead variant MAF range = 0.03-0.45, mean = 0.23) were only associated with one oxylipin or fatty acid, demonstrating evidence of distinct genetic effects. Finally, while most loci identified in two-degree-of-freedom tests were previously identified in our main effects analyses, we also identified an additional rare variant (MAF = 0.002) near CARS2, a locus previously implicated in inflammation. Our analyses revealed shared and distinct genetic architecture underlying fatty acids and oxylipins, providing insights into genetic factors and motivating future multi-omics work to characterize these compounds and elucidate their roles in disease pathways.

Lac-Phe mediates the effects of metformin on food intake and body weight

Metformin is a widely prescribed anti-diabetic medicine that also reduces body weight. There is ongoing debate about the mechanisms that mediate metformin's effects on energy balance. Here, we show that metformin is a powerful pharmacological inducer of the anorexigenic metabolite N-lactoyl-phenylalanine (Lac-Phe) in cells, in mice and two independent human cohorts. Metformin drives Lac-Phe biosynthesis through the inhibition of complex I, increased glycolytic flux and intracellular lactate mass action. Intestinal epithelial CNDP2+ cells, not macrophages, are the principal in vivo source of basal and metformin-inducible Lac-Phe. Genetic ablation of Lac-Phe biosynthesis in male mice renders animals resistant to the effects of metformin on food intake and body weight. Lastly, mediation analyses support a role for Lac-Phe as a downstream effector of metformin's effects on body mass index in participants of a large population-based observational cohort, the Multi-Ethnic Study of Atherosclerosis. Together, these data establish Lac-Phe as a critical mediator of the body weight-lowering effects of metformin.

Genetic factors associated with suicidal behaviors and alcohol use disorders in an American Indian population

American Indians (AI) demonstrate the highest rates of both suicidal behaviors (SB) and alcohol use disorders (AUD) among all ethnic groups in the US. Rates of suicide and AUD vary substantially between tribal groups and across different geographical regions, underscoring a need to delineate more specific risk and resilience factors. Using data from over 740 AI living within eight contiguous reservations, we assessed genetic risk factors for SB by investigating: (1) possible genetic overlap with AUD, and (2) impacts of rare and low-frequency genomic variants. Suicidal behaviors included lifetime history of suicidal thoughts and acts, including verified suicide deaths, scored using a ranking variable for the SB phenotype (range 0-4). We identified five loci significantly associated with SB and AUD, two of which are intergenic and three intronic on genes AACSP1, ANK1, and FBXO11. Nonsynonymous rare and low-frequency mutations in four genes including SERPINF1 (PEDF), ZNF30, CD34, and SLC5A9, and non-intronic rare and low-frequency mutations in genes OPRD1, HSD17B3 and one lincRNA were significantly associated with SB. One identified pathway related to hypoxia-inducible factor (HIF) regulation, whose 83 nonsynonymous rare and low-frequency variants on 10 genes were significantly linked to SB as well. Four additional genes, and two pathways related to vasopressin-regulated water metabolism and cellular hexose transport, also were strongly associated with SB. This study represents the first investigation of genetic factors for SB in an American Indian population that has high risk for suicide. Our study suggests that bivariate association analysis between comorbid disorders can increase statistical power; and rare and low-frequency variant analysis in a high-risk population enabled by whole-genome sequencing has the potential to identify novel genetic factors. Although such findings may be population specific, rare functional mutations relating to PEDF and HIF regulation align with past reports and suggest a biological mechanism for suicide risk and a potential therapeutic target for intervention.

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.

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.

Polygenic Risk Score for Cardiovascular Diseases in Artificial Intelligence Paradigm: A Review

Cardiovascular disease (CVD) related mortality and morbidity heavily strain society. The relationship between external risk factors and our genetics have not been well established. It is widely acknowledged that environmental influence and individual behaviours play a significant role in CVD vulnerability, leading to the development of polygenic risk scores (PRS). We employed the PRISMA search method to locate pertinent research and literature to extensively review artificial intelligence (AI)-based PRS models for CVD risk prediction. Furthermore, we analyzed and compared conventional vs. AI-based solutions for PRS. We summarized the recent advances in our understanding of the use of AI-based PRS for risk prediction of CVD. Our study proposes three hypotheses: i) Multiple genetic variations and risk factors can be incorporated into AI-based PRS to improve the accuracy of CVD risk predicting. ii) AI-based PRS for CVD circumvents the drawbacks of conventional PRS calculators by incorporating a larger variety of genetic and non-genetic components, allowing for more precise and individualised risk estimations. iii) Using AI approaches, it is possible to significantly reduce the dimensionality of huge genomic datasets, resulting in more accurate and effective disease risk prediction models. Our study highlighted that the AI-PRS model outperformed traditional PRS calculators in predicting CVD risk. Furthermore, using AI-based methods to calculate PRS may increase the precision of risk predictions for CVD and have significant ramifications for individualized prevention and treatment plans.

Lac-Phe mediates the anti-obesity effect of metformin

Metformin is a widely prescribed anti-diabetic medicine that also reduces body weight. The mechanisms that mediate metformin's effects on energy balance remain incompletely defined. Here we show that metformin is a powerful pharmacological inducer of the anorexigenic metabolite Lac-Phe in mice as well as in two independent human cohorts. In cell culture, metformin drives Lac-Phe biosynthesis via inhibition of complex I, increased glycolytic flux, and intracellular lactate mass action. Other biguanides and structurally distinct inhibitors of oxidative phosphorylation also increase Lac-Phe levels in vitro. Genetic ablation of CNDP2, the principal biosynthetic enzyme for Lac-Phe, in mice renders animals resistant to metformin's anorexigenic and anti-obesity effects. Mediation analyses also support a role for Lac-Phe in metformin's effect on body mass index in humans. These data establish the CNDP2/Lac-Phe pathway as a critical mediator of the effects of metformin on energy balance.

Disparities in cardio-oncology: Implication of angiogenesis, inflammation, and chemotherapy

Cancers and cardiovascular diseases are the top two causes of death in the United States. Over the past decades, novel therapies have slowed the cancer mortality rate, yet cardiac failures have risen due to the toxicity of cancer treatments. The mechanisms behind this relationship are poorly understood and it is crucial that we properly treat patients at risk of developing cardiac failure in response to cancer treatments. Currently, we rely on early-stage biomarkers of inflammation and angiogenesis to detect cardiotoxicity before it becomes irreversible. Identification of such biomarkers allows healthcare professionals to decrease the adverse effects of cancer therapies. Angiogenesis and inflammation have a systemic influence on the heart and vasculature following cancer therapy. In the field of cardio-oncology, there has been a recent emphasis on gender and racial disparities in cardiotoxicity and the impact of these disparities on disease outcomes, but there is a scarcity of data on how cardiotoxicity varies across diverse populations. Here, we will discuss how current markers of angiogenesis and inflammation induced by cancer therapy are related to disparities in cardiovascular health.

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.

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.

Genetic Factors Associated with Suicidal Behaviors and Alcohol Use Disorders in an American Indian Population

American Indians (AI) demonstrate the highest rates of both suicidal behaviors (SB) and alcohol use disorders (AUD) among all ethnic groups in the US. Rates of suicide and AUD vary substantially between tribal groups and across different geographical regions, underscoring a need to delineate more specific risk and resilience factors. Using data from over 740 AI living within eight contiguous reservations, we assessed genetic risk factors for SB by investigating: (1) possible genetic overlap with AUD, and (2) impacts of rare and low frequency genomic variants. Suicidal behaviors included lifetime history of suicidal thoughts and acts, including verified suicide deaths, scored using a ranking variable for the SB phenotype (range 0-4). We identified five loci significantly associated with SB and AUD, two of which are intergenic and three intronic on genes AACSP1, ANK1, and FBXO11. Nonsynonymous rare mutations in four genes including SERPINF1 (PEDF), ZNF30, CD34, and SLC5A9, and non-intronic rare mutations in genes OPRD1, HSD17B3 and one lincRNA were significantly associated with SB. One identified pathway related to hypoxia-inducible factor (HIF) regulation, whose 83 nonsynonymous rare variants on 10 genes were significantly linked to SB as well. Four additional genes, and two pathways related to vasopressin-regulated water metabolism and cellular hexose transport, also were strongly associated with SB. This study represents the first investigation of genetic factors for SB in an American Indian population that has high risk for suicide. Our study suggests that bivariate association analysis between comorbid disorders can increase statistical power; and rare variant analysis in a high-risk population enabled by whole-genome sequencing has the potential to identify novel genetic factors. Although such findings may be population specific, rare functional mutations relating to PEDF and HIF regulation align with past reports and suggest a biological mechanism for suicide risk and a potential therapeutic target for intervention.

CYP4F2 is a human-specific determinant of circulating N-acyl amino acid levels

N-acyl amino acids are a large family of circulating lipid metabolites that modulate energy expenditure and fat mass in rodents. However, little is known about the regulation and potential cardiometabolic functions of N-acyl amino acids in humans. Here, we analyze the cardiometabolic phenotype associations and genomic associations of four plasma N-acyl amino acids (N-oleoyl-leucine, N-oleoyl-phenylalanine, N-oleoyl-serine, and N-oleoyl-glycine) in 2,351 individuals from the Jackson Heart Study. We find that plasma levels of specific N-acyl amino acids are associated with cardiometabolic disease endpoints independent of free amino acid plasma levels and in patterns according to the amino acid head group. By integrating whole genome sequencing data with N-acyl amino acid levels, we identify that the genetic determinants of N-acyl amino acid levels also cluster according to amino acid head group. Furthermore, we identify the CYP4F2 locus as a genetic determinant of plasma N-oleoyl-leucine and N-oleoyl-phenylalanine levels in human plasma. In experimental studies, we demonstrate that CYP4F2-mediated hydroxylation of N-oleoyl-leucine and N-oleoyl-phenylalanine results in metabolic diversification and production of many previously unknown lipid metabolites with varying characteristics of the fatty acid tail group, including several that structurally resemble fatty acid hydroxy fatty acids (FAHFAs). These studies provide a structural framework for understanding the regulation and disease-associations of N-acyl amino acids in humans and identify that the diversity of this lipid signaling family can be significantly expanded through CYP4F-mediated ω-hydroxylation.

CYP4F2 is a human-specific determinant of circulating N-acyl amino acid levels

N-acyl amino acids are a large family of circulating lipid metabolites that modulate energy expenditure and fat mass in rodents. However, little is known about the regulation and potential cardiometabolic functions of N-acyl amino acids in humans. Here, we analyze the cardiometabolic phenotype associations and genetic regulation of four plasma N-fatty acyl amino acids (N-oleoyl-leucine, N-oleoyl-phenylalanine, N-oleoyl-serine, and N-oleoyl-glycine) in 2,351 individuals from the Jackson Heart Study. N-oleoyl-leucine and N-oleoyl-phenylalanine were positively associated with traits related to energy balance, including body mass index, waist circumference, and subcutaneous adipose tissue. In addition, we identify the CYP4F2 locus as a human-specific genetic determinant of plasma N-oleoyl-leucine and N-oleoyl-phenylalanine levels. In vitro, CYP4F2-mediated hydroxylation of N-oleoyl-leucine and N-oleoyl-phenylalanine results in metabolic diversification and production of many previously unknown lipid metabolites with varying characteristics of the fatty acid tail group, including several that structurally resemble fatty acid hydroxy fatty acids (FAHFAs). By contrast, FAAH-regulated N-oleoyl-glycine and N-oleoyl-serine were inversely associated with traits related to glucose and lipid homeostasis. These data uncover a human-specific enzymatic node for the metabolism of a subset of N-fatty acyl amino acids and establish a framework for understanding the cardiometabolic roles of individual N-fatty acyl amino acids in humans.

Multivariate GWAS analysis reveals loci associated with liver functions in continental African populations

BACKGROUND

Liver disease is any condition that causes liver damage and inflammation and may likely affect the function of the liver. Vital biochemical screening tools that can be used to evaluate the health of the liver and help diagnose, prevent, monitor, and control the development of liver disease are known as liver function tests (LFT). LFTs are performed to estimate the level of liver biomarkers in the blood. Several factors are associated with differences in concentration levels of LFTs in individuals, such as genetic and environmental factors. The aim of our study was to identify genetic loci associated with liver biomarker levels with a shared genetic basis in continental Africans, using a multivariate genome-wide association study (GWAS) approach.

METHODS

We used two distinct African populations, the Ugandan Genome Resource (UGR = 6,407) and South African Zulu cohort (SZC = 2,598). The six LFTs used in our analysis were: aspartate transaminase (AST), alanine transaminase (ALT), alkaline phosphatase (ALP), gamma-glutamyl transferase (GGT), total bilirubin, and albumin. A multivariate GWAS of LFTs was conducted using the exact linear mixed model (mvLMM) approach implemented in GEMMA and the resulting P-values were presented in Manhattan and quantile-quantile (QQ) plots. First, we attempted to replicate the findings of the UGR cohort in SZC. Secondly, given that the genetic architecture of UGR is different from that of SZC, we further undertook similar analysis in the SZC and discussed the results separately.

RESULTS

A total of 59 SNPs reached genome-wide significance (P = 5x10-8) in the UGR cohort and with 13 SNPs successfully replicated in SZC. These included a novel lead SNP near the RHPN1 locus (lead SNP rs374279268, P-value = 4.79x10-9, Effect Allele Frequency (EAF) = 0.989) and a lead SNP at the RGS11 locus (lead SNP rs148110594, P-value = 2.34x10-8, EAF = 0.928). 17 SNPs were significant in the SZC, while all the SNPs fall within a signal on chromosome 2, rs1976391 mapped to UGT1A was identified as the lead SNP within this region.

CONCLUSIONS

Using multivariate GWAS method improves the power to detect novel genotype-phenotype associations for liver functions not found with the standard univariate GWAS in the same dataset.