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

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.

The origins and functional effects of postzygotic mutations throughout the human life span

Postzygotic mutations (PZMs) begin to accrue in the human genome immediately after fertilization, but how and when PZMs affect development and lifetime health remain unclear. To study the origins and functional consequences of PZMs, we generated a multitissue atlas of PZMs spanning 54 tissue and cell types from 948 donors. Nearly half the variation in mutation burden among tissue samples can be explained by measured technical and biological effects, and 9% can be attributed to donor-specific effects. Through phylogenetic reconstruction of PZMs, we found that their type and predicted functional impact vary during prenatal development, across tissues, and through the germ cell life cycle. Thus, methods for interpreting effects across the body and the life span are needed to fully understand the consequences of genetic variants.

Multiomic signatures of body mass index identify heterogeneous health phenotypes and responses to a lifestyle intervention

Multiomic profiling can reveal population heterogeneity for both health and disease states. Obesity drives a myriad of metabolic perturbations and is a risk factor for multiple chronic diseases. Here we report an atlas of cross-sectional and longitudinal changes in 1,111 blood analytes associated with variation in body mass index (BMI), as well as multiomic associations with host polygenic risk scores and gut microbiome composition, from a cohort of 1,277 individuals enrolled in a wellness program (Arivale). Machine learning model predictions of BMI from blood multiomics captured heterogeneous phenotypic states of host metabolism and gut microbiome composition better than BMI, which was also validated in an external cohort (TwinsUK). Moreover, longitudinal analyses identified variable BMI trajectories for different omics measures in response to a healthy lifestyle intervention; metabolomics-inferred BMI decreased to a greater extent than actual BMI, whereas proteomics-inferred BMI exhibited greater resistance to change. Our analyses further identified blood analyte-analyte associations that were modified by metabolomics-inferred BMI and partially reversed in individuals with metabolic obesity during the intervention. Taken together, our findings provide a blood atlas of the molecular perturbations associated with changes in obesity status, serving as a resource to quantify metabolic health for predictive and preventive medicine.

Causal effects of circulating vitamin levels on the risk of heart failure: a Mendelian randomization study

BACKGROUND

Observational studies suggest inverse associations between serum vitamin levels and the risk of heart failure (HF). However, the causal effects of vitamins on HF have not been fully elucidated. Here, we conducted a Mendelian randomization (MR) study to investigate the causal associations between genetically determined vitamin levels and HF.

METHODS

Genetic instrumental variables for circulating vitamin levels, including vitamins A, B, C, D, and E, which were assessed as either absolute or metabolite levels were obtained from public genome-wide association studies. Summary statistics for single-nucleotide-polymorphisms and HF associations were retrieved from the HERMES Consortium (47,309 cases and 930,014 controls) and FinnGen Study (30,098 cases and 229,612 controls). Two-sample MR analyses were implemented to assess the causality between vitamin levels and HF per outcome database, and the results were subsequently combined by meta-analysis.

RESULTS

Our MR study did not find significant associations between genetically determined circulating vitamin levels and HF risk. For absolute vitamin levels, the odds ratio for HF ranged from 0.97 (95% confidence interval [CI]: 0.85-1.09, P = 0.41) for vitamin C to 1.05 (95% CI: 0.61-1.82, P = 0.85) for vitamin A. For vitamin metabolites, the odds ratio ranged between 0.94 (95% CI: 0.75-1.19, P = 0.62) for α-tocopherol and 1.11 (95% CI: 0.98-1.26, P = 0.09) for γ-tocopherol.

CONCLUSION

Evidence from our study does not support the causal effects of circulating vitamin levels on HF. Therefore, there may be no direct beneficial effects of vitamin intake on the prevention of primary HF.

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

Studies combining metabolomics and genetics, known as metabolite genome-wide association studies (mGWAS), have provided valuable insights into our understanding of the genetic control of metabolite levels. However, the biological interpretation of these associations remains challenging due to a lack of existing tools to annotate mGWAS gene-metabolite pairs beyond the use of conservative statistical significance threshold. Here, we computed the shortest reactional distance (SRD) based on the curated knowledge of the KEGG database to explore its utility in enhancing the biological interpretation of results from three independent mGWAS, including a case study on sickle cell disease patients. Results show that, in reported mGWAS pairs, there is an excess of small SRD values and that SRD values and p-values significantly correlate, even beyond the standard conservative thresholds. The added-value of SRD annotation is shown for identification of potential false negative hits, exemplified by the finding of gene-metabolite associations with SRD ≤1 that did not reach standard genome-wide significance cut-off. The wider use of this statistic as an mGWAS annotation would prevent the exclusion of biologically relevant associations and can also identify errors or gaps in current metabolic pathway databases. Our findings highlight the SRD metric as an objective, quantitative and easy-to-compute annotation for gene-metabolite pairs that can be used to integrate statistical evidence to biological networks.

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.

Uncovering the complex genetic architecture of human plasma lipidome using machine learning methods

Genetic architecture of plasma lipidome provides insights into regulation of lipid metabolism and related diseases. We applied an unsupervised machine learning method, PGMRA, to discover phenotype-genotype many-to-many relations between genotype and plasma lipidome (phenotype) in order to identify the genetic architecture of plasma lipidome profiled from 1,426 Finnish individuals aged 30-45 years. PGMRA involves biclustering genotype and lipidome data independently followed by their inter-domain integration based on hypergeometric tests of the number of shared individuals. Pathway enrichment analysis was performed on the SNP sets to identify their associated biological processes. We identified 93 statistically significant (hypergeometric p-value < 0.01) lipidome-genotype relations. Genotype biclusters in these 93 relations contained 5977 SNPs across 3164 genes. Twenty nine of the 93 relations contained genotype biclusters with more than 50% unique SNPs and participants, thus representing most distinct subgroups. We identified 30 significantly enriched biological processes among the SNPs involved in 21 of these 29 most distinct genotype-lipidome subgroups through which the identified genetic variants can influence and regulate plasma lipid related metabolism and profiles. This study identified 29 distinct genotype-lipidome subgroups in the studied Finnish population that may have distinct disease trajectories and therefore could be useful in precision medicine research.

Escape from X-inactivation in twins exhibits intra- and inter-individual variability across tissues and is heritable

X-chromosome inactivation (XCI) silences one X in female cells to balance sex-differences in X-dosage. A subset of X-linked genes escape XCI, but the extent to which this phenomenon occurs and how it varies across tissues and in a population is as yet unclear. To characterize incidence and variability of escape across individuals and tissues, we conducted a transcriptomic study of escape in adipose, skin, lymphoblastoid cell lines and immune cells in 248 healthy individuals exhibiting skewed XCI. We quantify XCI escape from a linear model of genes' allelic fold-change and XIST-based degree of XCI skewing. We identify 62 genes, including 19 lncRNAs, with previously unknown patterns of escape. We find a range of tissue-specificity, with 11% of genes escaping XCI constitutively across tissues and 23% demonstrating tissue-restricted escape, including cell type-specific escape across immune cells of the same individual. We also detect substantial inter-individual variability in escape. Monozygotic twins share more similar escape than dizygotic twins, indicating that genetic factors may underlie inter-individual differences in escape. However, discordant escape also occurs within monozygotic co-twins, suggesting environmental factors also influence escape. Altogether, these data indicate that XCI escape is an under-appreciated source of transcriptional differences, and an intricate phenotype impacting variable trait expressivity in females.

Prognostic 7-SLC-Gene Signature Identified via Weighted Gene Co-Expression Network Analysis for Patients with Hepatocellular Carcinoma

Background

Solute carrier (SLC) proteins play an important role in tumor metabolism. But SLC-associated genes' prognostic significance in hepatocellular carcinoma (HCC) remained elusive. We identified SLC-related factors and developed an SLC-related classifier to predict and improve HCC prognosis and treatment.

Methods

From the TCGA database, corresponding clinical data and mRNA expression profiles of 371 HCC patients were acquired, and those of 231 tumor samples were derived from the ICGC database. Genes associated with clinical features were filtered using weighted gene correlation network analysis (WGCNA). Next, univariate LASSO Cox regression studies developed SLC risk profiles, with the ICGC cohort data being used in validation.

Result

Univariate Cox regression analysis revealed that 31 SLC genes (P < 0.05) were related to HCC prognosis. 7 (SLC22A25, SLC2A2, SLC41A3, SLC44A1, SLC48A1, SLC4A2, and SLC9A3R1) of these genes were applied in developing a SLC gene prognosis model. Samples were classified into the low-andhigh-risk groups by the prognostic signature, with those in the high-risk group showing a significantly worse prognosis (P < 0.001 in the TCGA cohort and P=0.0068 in the ICGC cohort). ROC analysis validated the signature's prediction power. In addition, functional analyses showed enrichment of immune-related pathways and different immune status between the two risk groups.

Conclusion

The 7-SLC-gene prognostic signature established in this study helped predict the prognosis, and was also correlated with the tumor immune status and infiltration of different immune cells in the tumor microenvironment. The current findings may provide important clinical indications for proposing a novel combination therapy consists of targeted anti-SLC therapy and immunotherapy for HCC patients.

Host metabolome and faecal microbiome shows potential interactions impacted by age and weaning times in calves

BACKGROUND

Calves undergo nutritional, metabolic, and behavioural changes from birth to the entire weaning period. An appropriate selection of weaning age is essential to reduce the negative effects caused by weaning-related dietary transitions. This study monitored the faecal microbiome and plasma metabolome of 59 female Holstein calves during different developmental stages and weaning times (early vs. late) and identified the potential associations of the measured parameters over an experimental period of 140 days.

RESULTS

A progressive development of the microbiome and metabolome was observed with significant differences according to the weaning groups (weaned at 7 or 17 weeks of age). Faecal samples of young calves were dominated by bifidobacterial and lactobacilli species, while their respective plasma samples showed high concentrations of amino acids (AAs) and biogenic amines (BAs). However, as the calves matured, the abundances of potential fiber-degrading bacteria and the plasma concentrations of sphingomyelins (SMs), few BAs and acylcarnitines (ACs) were increased. Early-weaning at 7 weeks significantly restructured the microbiome towards potential fiber-degrading bacteria and decreased plasma concentrations of most of the AAs and SMs, few BAs and ACs compared to the late-weaning event. Strong associations between faecal microbes, plasma metabolites and calf growth parameters were observed during days 42-98, where the abundances of Bacteroides, Parabacteroides, and Blautia were positively correlated with the plasma concentrations of AAs, BAs and SMs as well as the live weight gain or average daily gain in calves.

CONCLUSION

The present study reported that weaning at 17 weeks of age was beneficial due to higher growth rate of late-weaned calves during days 42-98 and a quick adaptability of microbiota to weaning-related dietary changes during day 112, suggesting an age-dependent maturation of the gastrointestinal tract. However, the respective plasma samples of late-weaned calves contained several metabolites with differential concentrations to the early-weaned group, suggesting a less abrupt but more-persistent effect of dietary changes on host metabolome compared to the microbiome.

Associations between genetically predicted levels of blood metabolites and pancreatic cancer risk

Pancreatic ductal adenocarcinoma (PDAC) is one of the most aggressive solid malignancies, which is featured by systematic metabolism. Thus, a better understanding of metabolic dysregulation in PDAC is important to better characterize its etiology. Here, we performed a large metabolome-wide association study (MWAS) to systematically explore associations between genetically predicted metabolite levels in blood and PDAC risk. Using data from 881 subjects of European descent in the TwinsUK Project, comprehensive genetic models were built to predict serum metabolite levels. These prediction models were applied to the genetic data of 8275 cases and 6723 controls included in the PanScan (I, II and III) and PanC4 consortia. After assessing the metabolite-PDAC risk associations by a slightly modified TWAS/FUSION framework, we identified five metabolites (including two dipeptides) showing significant associations with PDAC risk at false discovery rate (FDR) <0.05. Integrated with gut microbial information, two-sample Mendelian randomization (MR) analyses were further performed to investigate the relationship among serum metabolites, gut microbiome features and PDAC. The flavonoid-degrading bacteria Flavonifractor sp90199495 was found to be associated with metabolite X-21849 and it was also shown to be associated with PDAC risk. Collectively, our study identified novel candidate metabolites for PDAC risk, which could lead to new insights into the etiology of PDAC and improved treatment options.

Combined plasma metabolomic and transcriptomic analysis identify histidine as a biomarker and potential contributor in SLE pathogenesis

OBJECTIVES

To investigate metabolite alterations in the plasma of SLE patients to identify novel biomarkers and provide insight into SLE pathogenesis.

METHODS

Patients with SLE (n = 41, discovery cohort and n = 37, replication cohort), healthy controls (n = 30 and n = 29) and patients with RA (n = 19, disease control) were recruited. Metabolic profiles of the plasma samples were analysed using liquid chromatography-time-of-flight mass spectrometry and capillary electrophoresis-time-of-flight mass spectrometry. Transcriptome data was analysed using RNA-sequencing for 18 immune cell subsets. The importance of histidine (His) in plasmablast differentiation was investigated by using mouse splenic B cells.

RESULTS

We demonstrate that a specific amino acid combination including His can effectively distinguish between SLE patients and healthy controls. Random forest and partial least squares-discriminant analysis identified His as an effective classifier for SLE patients. A decrease in His plasma levels correlated with damage accrual independent of prednisolone dosage and type I IFN signature. The oxidative phosphorylation signature in plasmablasts negatively correlated with His levels. We also showed that plasmablast differentiation induced by innate immune signals was dependent on His.

CONCLUSIONS

Plasma His levels are a potential biomarker for SLE patients and are associated with damage accrual. Our data suggest the importance of His as a pathogenic metabolite in SLE pathogenesis.

Genetic influences on human blood metabolites in the Japanese population

An increase in ethnic diversity in genetic studies has the potential to provide unprecedented insights into how genetic variations influence human phenotypes. In this study, we conducted a quantitative trait locus (QTL) analysis of 121 metabolites measured using gas chromatography-mass spectrometry with plasma samples from 4,888 Japanese individuals. We found 60 metabolite-gene associations, of which 13 have not been previously reported. Meta-analyses with another Japanese and a European study identified six and two additional unreported loci, respectively. Genetic variants influencing metabolite levels were more enriched in protein-coding regions than in the regulatory regions while being associated with the risk of various diseases. Finally, we identified a signature of strong negative selection for uric acid ( S ˆ  = -1.53, p = 6.2 × 10^-18). Our study expanded the knowledge of genetic influences on human blood metabolites, providing valuable insights into their physiological, pathological, and selective properties.

Dissecting Causal Associations of Diet-Derived Circulating Antioxidants with Six Major Mental Disorders: A Mendelian Randomization Study

Although observational studies have suggested associations between circulating antioxidants and many mental disorders, causal inferences have not been confirmed. Mendelian randomization (MR) analyses were conducted using summary-level statistics from genome-wide association studies (GWASs) to explore whether genetically determined absolute circulating antioxidants (i.e., ascorbate, retinol, β-carotene, and lycopene) and metabolites (i.e., α- and γ-tocopherol, ascorbate, and retinol) were causally associated with the risk of six major mental disorders, including anxiety disorders (AD), major depressive disorder (MDD), bipolar disorder (BIP), schizophrenia (SCZ), post-traumatic stress disorder (PTSD), and obsessive-compulsive disorder (OCD). MR analyses were performed per specific-outcome databases, including the largest GWAS published to date (from 9725 for OCD to 413,466 for BIP participants), UK Biobank (over 370,000 participants), and FinnGen (over 270,000 participants), followed by meta-analyses. We found no significant evidence that genetically determined diet-derived circulating antioxidants were significantly causally associated with the risk of the six above-mentioned major mental disorders. For absolute antioxidant levels, the odds ratios (ORs) ranged from 0.91 (95% CI, 0.67-1.23) for the effect of β-carotene on OCD to 1.18 (95% CI, 0.90-1.54) for the effect of ascorbate on OCD. Similarly, for antioxidant metabolites, ORs ranged from 0.87 (95% CI, 0.55-1.38) for the effect of ascorbate on MDD to 1.08 (95% CI, 0.88-1.33) for the effect of ascorbate on OCD. Our study does not support significant causal associations of genetically determined diet-derived circulating antioxidants with the risk of major mental disorders.

Genomic atlas of the plasma metabolome prioritizes metabolites implicated in human diseases

Metabolic processes can influence disease risk and provide therapeutic targets. By conducting genome-wide association studies of 1,091 blood metabolites and 309 metabolite ratios, we identified associations with 690 metabolites at 248 loci and associations with 143 metabolite ratios at 69 loci. Integrating metabolite-gene and gene expression information identified 94 effector genes for 109 metabolites and 48 metabolite ratios. Using Mendelian randomization (MR), we identified 22 metabolites and 20 metabolite ratios having estimated causal effect on 12 traits and diseases, including orotate for estimated bone mineral density, α-hydroxyisovalerate for body mass index and ergothioneine for inflammatory bowel disease and asthma. We further measured the orotate level in a separate cohort and demonstrated that, consistent with MR, orotate levels were positively associated with incident hip fractures. This study provides a valuable resource describing the genetic architecture of metabolites and delivers insights into their roles in common diseases, thereby offering opportunities for therapeutic targets.

Comprehensive characterization of putative genetic influences on plasma metabolome in a pediatric cohort

BACKGROUND

The human exposome is composed of diverse metabolites and small chemical compounds originated from endogenous and exogenous sources, respectively. Genetic and environmental factors influence metabolite levels, while the extent of genetic contributions across metabolic pathways is not yet known. Untargeted profiling of human metabolome using high-resolution mass spectrometry (HRMS) combined with genome-wide genotyping allows comprehensive identification of genetically influenced metabolites. As such previous studies of adults discovered and replicated genotype-metabotype associations. However, these associations have not been characterized in children.

RESULTS

We conducted the largest genome by metabolome-wide association study to date of children (N = 441) using 619,688 common genetic variants and 14,342 features measured by HRMS. Narrow-sense heritability (h²) estimates of plasma metabolite concentrations using genomic relatedness matrix restricted maximum likelihood (GREML) method showed a bimodal distribution with high h² (> 0.8) for 15.9% of features and low h² (< 0.2) for most of features (62.0%). The features with high h² were enriched for amino acid and nucleic acid metabolism, while carbohydrate and lipid concentrations showed low h². For each feature, a metabolite quantitative trait loci (mQTL) analysis was performed to identify genetic variants that were potentially associated with plasma levels. Fifty-four associations among 29 features and 43 genetic variants were identified at a genome-wide significance threshold p < 3.5 × 10^-12 (= 5 × 10^-8/14,342 features). Previously reported associations such as UGT1A1 and bilirubin; PYROXD2 and methyl lysine; and ACADS and butyrylcarnitine were successfully replicated in our pediatric cohort. We found potential candidates for novel associations including CSMD1 and a monostearyl alcohol triglyceride (m/z 781.7483, retention time (RT) 89.3 s); CALN1 and Tridecanol (m/z 283.2741, RT 27.6). A gene-level enrichment analysis using MAGMA revealed highly interconnected modules for dADP biosynthesis, sterol synthesis, and long-chain fatty acid transport in the gene-feature network.

CONCLUSION

Comprehensive profiling of plasma metabolome across age groups combined with genome-wide genotyping revealed a wide range of genetic influence on diverse chemical species and metabolic pathways. The developmental trajectory of a biological system is shaped by gene-environment interaction especially in early life. Therefore, continuous efforts on generating metabolomics data in diverse human tissue types across age groups are required to understand gene-environment interaction toward healthy aging trajectories.

Rare and common genetic determinants of metabolic individuality and their effects on human health

Garrod's concept of 'chemical individuality' has contributed to comprehension of the molecular origins of human diseases. Untargeted high-throughput metabolomic technologies provide an in-depth snapshot of human metabolism at scale. We studied the genetic architecture of the human plasma metabolome using 913 metabolites assayed in 19,994 individuals and identified 2,599 variant-metabolite associations (P < 1.25 × 10-11) within 330 genomic regions, with rare variants (minor allele frequency ≤ 1%) explaining 9.4% of associations. Jointly modeling metabolites in each region, we identified 423 regional, co-regulated, variant-metabolite clusters called genetically influenced metabotypes. We assigned causal genes for 62.4% of these genetically influenced metabotypes, providing new insights into fundamental metabolite physiology and clinical relevance, including metabolite-guided discovery of potential adverse drug effects (DPYD and SRD5A2). We show strong enrichment of inborn errors of metabolism-causing genes, with examples of metabolite associations and clinical phenotypes of non-pathogenic variant carriers matching characteristics of the inborn errors of metabolism. Systematic, phenotypic follow-up of metabolite-specific genetic scores revealed multiple potential etiological relationships.

APOE alleles modulate associations of plasma metabolites with variants from multiple genes on chromosome 19q13.3

The APOE ε2, ε3, and ε4 alleles differentially impact various complex diseases and traits. We examined whether these alleles modulated associations of 94 single-nucleotide polymorphisms (SNPs) harbored by 26 genes in 19q13.3 region with 217 plasma metabolites using Framingham Heart Study data. The analyses were performed in the E2 (ε2ε2 or ε2ε3 genotype), E3 (ε3ε3 genotype), and E4 (ε3ε4 or ε4ε4 genotype) groups separately. We identified 31, 17, and 22 polymorphism-metabolite associations in the E2, E3, and E4 groups, respectively, at a false discovery rate P FDR < 0.05. These entailed 51 and 19 associations with 20 lipid and 12 polar analytes. Contrasting the effect sizes between the analyzed groups showed 20 associations with group-specific effects at Bonferroni-adjusted P < 7.14E-04. Three associations with glutamic acid or dimethylglycine had significantly larger effects in the E2 than E3 group and 12 associations with triacylglycerol 56:5, lysophosphatidylethanolamines 16:0, 18:0, 20:4, or phosphatidylcholine 38:6 had significantly larger effects in the E2 than E4 group. Two associations with isocitrate or propionate and three associations with phosphatidylcholines 32:0, 32:1, or 34:0 had significantly larger effects in the E4 than E3 group. Nine of 70 SNP-metabolite associations identified in either E2, E3, or E4 groups attained P FDR < 0.05 in the pooled sample of these groups. However, none of them were among the 20 group-specific associations. Consistent with the evolutionary history of the APOE alleles, plasma metabolites showed higher APOE-cluster-related variations in the E4 than E2 and E3 groups. Pathway enrichment mainly highlighted lipids and amino acids metabolism and citrate cycle, which can be differentially impacted by the APOE alleles. These novel findings expand insights into the genetic heterogeneity of plasma metabolites and highlight the importance of the APOE-allele-stratified genetic analyses of the APOE-related diseases and traits.

The metabolomics of human aging: Advances, challenges, and opportunities

As the global population becomes older, understanding the impact of aging on health and disease becomes paramount. Recent advancements in multiomic technology have allowed for the high-throughput molecular characterization of aging at the population level. Metabolomics studies that analyze the small molecules in the body can provide biological information across a diversity of aging processes. Here, we review the growing body of population-scale metabolomics research on aging in humans, identifying the major trends in the field, implicated biological pathways, and how these pathways relate to health and aging. We conclude by assessing the main challenges in the research to date, opportunities for advancing the field, and the outlook for precision health applications.

Differences in Metabolomic Profiles Between Black and White Women and Risk of Coronary Heart Disease: an Observational Study of Women From Four US Cohorts

BACKGROUND

Racial differences in metabolomic profiles may reflect underlying differences in social determinants of health by self-reported race and may be related to racial disparities in coronary heart disease (CHD) among women in the United States. However, the magnitude of differences in metabolomic profiles between Black and White women in the United States has not been well-described. It also remains unknown whether such differences are related to differences in CHD risk.

METHODS

Plasma metabolomic profiles were analyzed using liquid chromatography-tandem mass spectrometry in the WHI-OS (Women's Health Initiative-Observational Study; 138 Black and 696 White women), WHI-HT trials (WHI-Hormone Therapy; 156 Black and 1138 White women), MESA (Multi-Ethnic Study of Atherosclerosis; 114 Black and 219 White women), JHS (Jackson Heart Study; 1465 Black women with 107 incident CHD cases), and NHS (Nurses' Health Study; 2506 White women with 136 incident CHD cases). First, linear regression models were used to estimate associations between self-reported race and 472 metabolites in WHI-OS (discovery); findings were replicated in WHI-HT and validated in MESA. Second, we used elastic net regression to construct a racial difference metabolomic pattern (RDMP) representing differences in the metabolomic patterns between Black and White women in the WHI-OS; the RDMP was validated in the WHI-HT and MESA. Third, using conditional logistic regressions in the WHI (717 CHD cases and 719 matched controls), we examined associations of metabolites with large differences in levels by race and the RDMP with risk of CHD, and the results were replicated in Black women from the JHS and White women from the NHS.

RESULTS

Of the 472 tested metabolites, levels of 259 (54.9%) metabolites, mostly lipid metabolites and amino acids, significantly differed between Black and White women in both WHI-OS and WHI-HT after adjusting for baseline characteristics, socioeconomic status, lifestyle factors, baseline health conditions, and medication use (false discovery rate <0.05); similar trends were observed in MESA. The RDMP, composed of 152 metabolites, was identified in the WHI-OS and showed significantly different distributions between Black and White women in the WHI-HT and MESA. Higher RDMP quartiles were associated with an increased risk of incident CHD (odds ratio=1.51 [0.97-2.37] for the highest quartile comparing to the lowest; Ptrend=0.02), independent of self-reported race and known CHD risk factors. In race-stratified analyses, the RDMP-CHD associations were more pronounced in White women. Similar patterns were observed in Black women from the JHS and White women from the NHS.

CONCLUSIONS

Metabolomic profiles significantly and substantially differ between Black and White women and may be associated with CHD risk and racial disparities in US women.

Gut microbiome of multiple sclerosis patients and paired household healthy controls reveal associations with disease risk and course

Changes in gut microbiota have been associated with several diseases. Here, the International Multiple Sclerosis Microbiome Study (iMSMS) studied the gut microbiome of 576 MS patients (36% untreated) and genetically unrelated household healthy controls (1,152 total subjects). We observed a significantly increased proportion of Akkermansia muciniphila, Ruthenibacterium lactatiformans, Hungatella hathewayi, and Eisenbergiella tayi and decreased Faecalibacterium prausnitzii and Blautia species. The phytate degradation pathway was over-represented in untreated MS, while pyruvate-producing carbohydrate metabolism pathways were significantly reduced. Microbiome composition, function, and derived metabolites also differed in response to disease-modifying treatments. The therapeutic activity of interferon-β may in part be associated with upregulation of short-chain fatty acid transporters. Distinct microbial networks were observed in untreated MS and healthy controls. These results strongly support specific gut microbiome associations with MS risk, course and progression, and functional changes in response to treatment.

The UK Biobank: A Shining Example of Genome-Wide Association Study Science with the Power to Detect the Murky Complications of Real-World Epidemiology

Genome-wide association studies (GWASs) have successfully identified thousands of genetic variants that are reliably associated with human traits. Although GWASs are restricted to certain variant frequencies, they have improved our understanding of the genetic architecture of complex traits and diseases. The UK Biobank (UKBB) has brought substantial analytical opportunity and performance to association studies. The dramatic expansion of many GWAS sample sizes afforded by the inclusion of UKBB data has improved the power of estimation of effect sizes but, critically, has done so in a context where phenotypic depth and precision enable outcome dissection and the application of epidemiological approaches. However, at the same time, the availability of such a large, well-curated, and deeply measured population-based collection has the capacity to increase our exposure to the many complications and inferential complexities associated with GWASs and other analyses. In this review, we discuss the impact that UKBB has had in the GWAS era, some of the opportunities that it brings, and exemplar challenges that illustrate the reality of using data from this world-leading resource.

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

Integrating genetic information with metabolomics has provided new insights into genes affecting human metabolism. However, gene-metabolite integration has been primarily studied in individuals of European Ancestry, limiting the opportunity to leverage genomic diversity for discovery. In addition, these analyses have principally involved known metabolites, with the majority of the profiled peaks left unannotated. Here, we perform a whole genome association study of 2,291 metabolite peaks (known and unknown features) in 2,466 Black individuals from the Jackson Heart Study. We identify 519 locus-metabolite associations for 427 metabolite peaks and validate our findings in two multi-ethnic cohorts. A significant proportion of these associations are in ancestry specific alleles including findings in APOE, TTR and CD36. We leverage tandem mass spectrometry to annotate unknown metabolites, providing new insight into hereditary diseases including transthyretin amyloidosis and sickle cell disease. Our integrative omics approach leverages genomic diversity to provide novel insights into diverse cardiometabolic diseases.

Diet-Derived Circulating Antioxidants and Risk of Digestive System Tumors: A Mendelian Randomization Study

Previous observational case-control studies have shown significant controversy over the impact of dietary intake-related circulating antioxidants on the risk of digestive system tumors. We conducted a two-sample Mendelian randomized (MR) analysis to determine whether there was a significant causal relationship between increased levels of circulating antioxidants and digestive system tumors. Our circulating antioxidants (vitamin C, carotenoids, vitamin A, and vitamin E) were derived from absolute circulating antioxidants and circulating antioxidant metabolites, and their corresponding instrumental variables were screened from published studies. The digestive system tumors we studied included colorectal, gastric, pancreatic, liver, and esophageal cancer, and the corresponding summary GAWS (genome-wide association study) data were obtained from the UK Biobank database. We first evaluated the causal relationship between each tumor and circulating antioxidants and then used meta-analysis to summarize the results of MR analysis of different tumors. No significant associations were noted for genetically predicted circulating antioxidants and higher risk of digestive system tumors in our study. The pooled ORs (odds ratio) are 0.72 (95% CI: 0.46-1.11; β-carotene), 0.93 (95% CI: 0.81-1.08; lycopene), 2.12 (95% CI: 0.31-14.66; retinol), and 0.99 (95% CI: 0.96-1.02; ascorbate) for absolute circulating antioxidants; for circulating antioxidant metabolites, the pooled ORs for digestive system tumors risk per unit increase of antioxidants were 1.29 (95% CI: 0.39-4.28; α-tocopherol), 1.72 (95% CI: 0.85-3.49; γ-tocopherol), 1.05 (95% CI: 0.96-1.14; retinol), and 1.21 (95% CI: 0.97-1.51; ascorbate), respectively. Our study suggested that increased levels of dietary-derived circulating antioxidants did not reduce the risk of digestive system tumors.

Investigating Causal Associations of Diet-Derived Circulating Antioxidants with the Risk of Digestive System Cancers: A Mendelian Randomization Study

Molecular mechanisms and observational studies have found that diet-derived antioxidants are associated with digestive system cancers, whereas there is a lack of causal evidence from randomized clinical trials. In this study, we aimed to assess the causality of these associations through a Mendelian randomization (MR) study. Single nucleotide polymorphisms of diet-derived circulating antioxidants (i.e., α- and γ-tocopherol, ascorbate, retinol, β-carotene, lycopene, and urate), accessed by absolute levels and relative metabolite concentrations, were used as genetic instruments. Summary statistics for digestive system cancers were obtained from the UK Biobank and FinnGen studies. Two-sample MR analyses were performed in each of the two outcome databases, followed by a meta-analysis. The inverse-variance weighted MR was adopted as the primary analysis. Five additional MR methods (likelihood-based MR, MR-Egger, weighted median, penalized weighted median, and MR-PRESSO) and replicate MR analyses for outcomes from different sources were used as sensitivity analyses. Genetically determined antioxidants were not significantly associated with five digestive system cancers, after correcting for multiple tests. However, we found suggestive evidence that absolute ascorbate levels were negatively associated with colon cancer in UK Biobank-the odds ratio (OR) per unit increase in ascorbate was 0.774 (95% confidence interval [CI] 0.608-0.985, p = 0.037), which was consistent with the results in FinnGen, and the combined OR was 0.764 (95% CI 0.623-0.936, p = 0.010). Likewise, higher absolute retinol levels suggestively reduced the pancreatic cancer risk in FinnGen-the OR per 10% unit increase in ln-transformed retinol was 0.705 (95% CI 0.529-0.940, p = 0.017), which was consistent with the results in UK Biobank and the combined OR was 0.747 (95% CI, 0.584-0.955, p = 0.020). Sensitivity analyses verified the above suggestive evidence. Our findings suggest that higher levels of antioxidants are unlikely to be a causal protective factor for most digestive system cancers, except for the suggestive protective effects of ascorbate on colon cancer and of retinol on pancreatic cancer.

Association of Physical Activity With Bioactive Lipids and Cardiovascular Events

BACKGROUND

To clarify the mechanisms underlying physical activity (PA)-related cardioprotection, we examined the association of PA with plasma bioactive lipids (BALs) and cardiovascular disease (CVD) events. We additionally performed genome-wide associations.

METHODS

PA-bioactive lipid associations were examined in VITAL (VITamin D and OmegA-3 TriaL)-clinical translational science center (REGISTRATION: URL: https://www.

CLINICALTRIALS

gov; Unique identifier: NCT01169259; N=1032) and validated in JUPITER (Justification for the Use of statins in Prevention: an Intervention Trial Evaluating Rosuvastatin)-NC (REGISTRATION: URL: https://www.

CLINICALTRIALS

gov; Unique identifier: NCT00239681; N=589), using linear models adjusted for age, sex, race, low-density lipoprotein-cholesterol, total-C, and smoking. Significant BALs were carried over to examine associations with incident CVD in 2 nested CVD case-control studies: VITAL-CVD (741 case-control pairs) and JUPITER-CVD (415 case-control pairs; validation).

RESULTS

We detected 145 PA-bioactive lipid validated associations (false discovery rate <0.1). Annotations were found for 6 of these BALs: 12,13-diHOME, 9,10-diHOME, lysoPC(15:0), oxymorphone-3b-D-glucuronide, cortisone, and oleoyl-glycerol. Genetic analysis within JUPITER-NC showed associations of 32 PA-related BALs with 22 single-nucleotide polymorphisms. From PA-related BALs, 12 are associated with CVD.

CONCLUSIONS

We identified a PA-related bioactive lipidome profile out of which 12 BALs also had opposite associations with incident CVD events.

Crosstalk between Host Genome and Metabolome among People with HIV in South Africa

Genome-wide association studies (GWAS) of circulating metabolites have revealed the role of genetic regulation on the human metabolome. Most previous investigations focused on European ancestry, and few studies have been conducted among populations of African descent living in Africa, where the infectious disease burden is high (e.g., human immunodeficiency virus (HIV)). It is important to understand the genetic associations of the metabolome in diverse at-risk populations including people with HIV (PWH) living in Africa. After a thorough literature review, the reported significant gene−metabolite associations were tested among 490 PWH in South Africa. Linear regression was used to test associations between the candidate metabolites and genetic variants. GWAS of 154 plasma metabolites were performed to identify novel genetic associations. Among the 29 gene−metabolite associations identified in the literature, we replicated 10 in South Africans with HIV. The UGT1A cluster was associated with plasma levels of biliverdin and bilirubin; SLC16A9 and CPS1 were associated with carnitine and creatine, respectively. We also identified 22 genetic associations with metabolites using a genome-wide significance threshold (p-value < 5 × 10−8). In a GWAS of plasma metabolites in South African PWH, we replicated reported genetic associations across ancestries, and identified novel genetic associations using a metabolomics approach.

Successful use of empagliflozin to treat neutropenia in two G6PC3-deficient children: Impact of a mutation in SGLT5

Neutropenia and neutrophil dysfunction found in deficiencies in G6PC3 and in the glucose-6-phosphate transporter (G6PT/SLC37A4) are due to accumulation of 1,5-anhydroglucitol-6-phosphate (1,5-AG6P), an inhibitor of hexokinase made from 1,5-anhydroglucitol (1,5-AG), an abundant polyol present in blood. Lowering blood 1,5-AG with an SGLT2 inhibitor greatly improved neutrophil counts and function in G6PC3-deficient mice and in patients with G6PT-deficiency. We evaluate this treatment in two G6PC3-deficient children. While neutropenia was severe in one child (PT1), which was dependent on granulocyte cololony-stimulating factor (GCSF), it was significantly milder in the other one (PT2), which had low blood 1,5-AG levels and only required GCSF during severe infections. Treatment with the SGLT2-inhibitor empagliflozin decreased 1,5-AG in blood and 1,5-AG6P in neutrophils and improved (PT1) or normalized (PT2) neutrophil counts, allowing to stop GCSF. On empagliflozin, both children remained infection-free (>1 year - PT2; >2 years - PT1) and no side effects were reported. Remarkably, sequencing of SGLT5, the gene encoding the putative renal transporter for 1,5-AG, disclosed a rare heterozygous missense mutation in PT2, replacing the extremely conserved Arg401 by a histidine. The higher urinary clearance of 1,5-AG explains the more benign neutropenia and the outstanding response to empagliflozin treatment found in this child. Our data shows that SGLT2 inhibitors are an excellent alternative to treat the neutropenia present in G6PC3-deficiency.

Genomics-based identification of a potential causal role for acylcarnitine metabolism in depression

BACKGROUND

Altered metabolism of acylcarnitines - transporting fatty acids to mitochondria - may link cellular energy dysfunction to depression. We examined the potential causal role of acylcarnitine metabolism in depression by leveraging genomics and Mendelian randomization.

METHODS

Summary statistics were obtained from large GWAS: the Fenland Study (N = 9363), and the Psychiatric Genomics Consortium (246,363 depression cases and 561,190 controls). Two-sample Mendelian randomization analyses tested the potential causal link of 15 endogenous acylcarnitines with depression.

RESULTS

In univariable analyses, genetically-predicted lower levels of short-chain acylcarnitines C2 (odds ratio [OR] 0.97, 95% confidence intervals [CIs] 0.95-1.00) and C3 (OR 0.97, 95%CIs 0.96-0.99) and higher levels of medium-chain acylcarnitines C8 (OR 1.04, 95%CIs 1.01-1.06) and C10 (OR 1.04, 95%CIs 1.02-1.06) were associated with increased depression risk. No reverse potential causal role of depression genetic liability on acylcarnitines levels was found. Multivariable analyses showed that the association with depression was driven by the medium-chain acylcarnitines C8 (OR 1.04, 95%CIs 1.02-1.06) and C10 (OR 1.04, 95%CIs 1.02-1.06), suggesting a potential causal role in the risk of depression. Causal estimates for C8 (OR = 1.05, 95%CIs = 1.02-1.07) and C10 (OR = 1.05, 95%CIs = 1.02-1.08) were confirmed in follow-up analyses using genetic instruments derived from a GWAS meta-analysis including up to 16,841 samples.

DISCUSSION

Accumulation of medium-chain acylcarnitines is a signature of inborn errors of fatty acid metabolism and age-related metabolic conditions. Our findings point to a link between altered mitochondrial energy production and depression pathogenesis. Acylcarnitine metabolism represents a promising access point for the development of novel therapeutic approaches for depression.

Genetic Architecture of Plasma Alpha-Aminoadipic Acid Reveals a Relationship With High-Density Lipoprotein Cholesterol

Background Elevated plasma levels of alpha-aminoadipic acid (2-AAA) have been associated with the development of type 2 diabetes and atherosclerosis. However, the nature of the association remains unknown. Methods and Results We identified genetic determinants of plasma 2-AAA through meta-analysis of genome-wide association study data in 5456 individuals of European, African, and Asian ancestry from the Framingham Heart Study, Diabetes Prevention Program, Jackson Heart Study, and Shanghai Women's and Men's Health Studies. No single nucleotide polymorphisms reached genome-wide significance across all samples. However, the top associations from the meta-analysis included single-nucleotide polymorphisms in the known 2-AAA pathway gene DHTKD1, and single-nucleotide polymorphisms in genes involved in mitochondrial respiration (NDUFS4) and macrophage function (MSR1). We used a Mendelian randomization instrumental variable approach to evaluate relationships between 2-AAA and cardiometabolic phenotypes in large disease genome-wide association studies. Mendelian randomization identified a suggestive inverse association between increased 2-AAA and lower high-density lipoprotein cholesterol (P=0.005). We further characterized the genetically predicted relationship through measurement of plasma 2-AAA and high-density lipoprotein cholesterol in 2 separate samples of individuals with and without cardiometabolic disease (N=98), and confirmed a significant negative correlation between 2-AAA and high-density lipoprotein (rs=-0.53, P<0.0001). Conclusions 2-AAA levels in plasma may be regulated, in part, by common variants in genes involved in mitochondrial and macrophage function. Elevated plasma 2-AAA associates with reduced levels of high-density lipoprotein cholesterol. Further mechanistic studies are required to probe this as a possible mechanism linking 2-AAA to future cardiometabolic risk.

Comprehensive genetic analysis of the human lipidome identifies loci associated with lipid homeostasis with links to coronary artery disease

We integrated lipidomics and genomics to unravel the genetic architecture of lipid metabolism and identify genetic variants associated with lipid species putatively in the mechanistic pathway for coronary artery disease (CAD). We quantified 596 lipid species in serum from 4,492 individuals from the Busselton Health Study. The discovery GWAS identified 3,361 independent lipid-loci associations, involving 667 genomic regions (479 previously unreported), with validation in two independent cohorts. A meta-analysis revealed an additional 70 independent genomic regions associated with lipid species. We identified 134 lipid endophenotypes for CAD associated with 186 genomic loci. Associations between independent lipid-loci with coronary atherosclerosis were assessed in ∼456,000 individuals from the UK Biobank. Of the 53 lipid-loci that showed evidence of association (P < 1 × 10-3), 43 loci were associated with at least one lipid endophenotype. These findings illustrate the value of integrative biology to investigate the aetiology of atherosclerosis and CAD, with implications for other complex diseases.

Whole-exome sequencing identifies rare genetic variants associated with human plasma metabolites

Metabolite levels measured in the human population are endophenotypes for biological processes. We combined sequencing data for 3,924 (whole-exome sequencing, WES, discovery) and 2,805 (whole-genome sequencing, WGS, replication) donors from a prospective cohort of blood donors in England. We used multiple approaches to select and aggregate rare genetic variants (minor allele frequency [MAF] < 0.1%) in protein-coding regions and tested their associations with 995 metabolites measured in plasma by using ultra-high-performance liquid chromatography-tandem mass spectrometry. We identified 40 novel associations implicating rare coding variants (27 genes and 38 metabolites), of which 28 (15 genes and 28 metabolites) were replicated. We developed algorithms to prioritize putative driver variants at each locus and used mediation and Mendelian randomization analyses to test directionality at associations of metabolite and protein levels at the ACY1 locus. Overall, 66% of reported associations implicate gene targets of approved drugs or bioactive drug-like compounds, contributing to drug targets' validating efforts.

Host and gut microbial tryptophan metabolism and type 2 diabetes: an integrative analysis of host genetics, diet, gut microbiome and circulating metabolites in cohort studies

OBJECTIVE

Tryptophan can be catabolised to various metabolites through host kynurenine and microbial indole pathways. We aimed to examine relationships of host and microbial tryptophan metabolites with incident type 2 diabetes (T2D), host genetics, diet and gut microbiota.

METHOD

We analysed associations between circulating levels of 11 tryptophan metabolites and incident T2D in 9180 participants of diverse racial/ethnic backgrounds from five cohorts. We examined host genome-wide variants, dietary intake and gut microbiome associated with these metabolites.

RESULTS

Tryptophan, four kynurenine-pathway metabolites (kynurenine, kynurenate, xanthurenate and quinolinate) and indolelactate were positively associated with T2D risk, while indolepropionate was inversely associated with T2D risk. We identified multiple host genetic variants, dietary factors, gut bacteria and their potential interplay associated with these T2D-relaetd metabolites. Intakes of fibre-rich foods, but not protein/tryptophan-rich foods, were the dietary factors most strongly associated with tryptophan metabolites. The fibre-indolepropionate association was partially explained by indolepropionate-associated gut bacteria, mostly fibre-using Firmicutes. We identified a novel association between a host functional LCT variant (determining lactase persistence) and serum indolepropionate, which might be related to a host gene-diet interaction on gut Bifidobacterium, a probiotic bacterium significantly associated with indolepropionate independent of other fibre-related bacteria. Higher milk intake was associated with higher levels of gut Bifidobacterium and serum indolepropionate only among genetically lactase non-persistent individuals.

CONCLUSION

Higher milk intake among lactase non-persistent individuals, and higher fibre intake were associated with a favourable profile of circulating tryptophan metabolites for T2D, potentially through the host-microbial cross-talk shifting tryptophan metabolism toward gut microbial indolepropionate production.

Ratios of Acetaminophen Metabolites Identify New Loci of Pharmacogenetic Relevance in a Genome-Wide Association Study

Genome-wide association studies (GWAS) with non-targeted metabolomics have identified many genetic loci of biomedical interest. However, metabolites with a high degree of missingness, such as drug metabolites and xenobiotics, are often excluded from such studies due to a lack of statistical power and higher uncertainty in their quantification. Here we propose ratios between related drug metabolites as GWAS phenotypes that can drastically increase power to detect genetic associations between pairs of biochemically related molecules. As a proof-of-concept we conducted a GWAS with 520 individuals from the Qatar Biobank for who at least five of the nine available acetaminophen metabolites have been detected. We identified compelling evidence for genetic variance in acetaminophen glucuronidation and methylation by UGT2A15 and COMT, respectively. Based on the metabolite ratio association profiles of these two loci we hypothesized the chemical structure of one of their products or substrates as being 3-methoxyacetaminophen, which we then confirmed experimentally. Taken together, our study suggests a novel approach to analyze metabolites with a high degree of missingness in a GWAS setting with ratios, and it also demonstrates how pharmacological pathways can be mapped out using non-targeted metabolomics measurements in large population-based studies.