Biochemical insights from population studies with genetics and metabolomics

Bioinformatics-driven untargeted metabolomic profiling for clinical screening of methamphetamine abuse

PURPOSE

Amphetamine-type stimulants are very common, and their usage is becoming a very big social problem all over the world. Thousands of addicts encounter several health problems including mental, metabolic, behavioral and neurological disorders. In addition to these, there are several reports about the elevated risk of tendency on committing criminal cases by addicted persons. Hence, methamphetamine addiction is not only an individual health problem but also a social problem. In our study, we aimed to investigate the pathogenesis of chronic usage of methamphetamine via untargeted metabolomics approach.

METHODS

38 plasma samples were carefully collected and extracted for untargeted metabolomics assay. A liquid-liquid extraction was performed to get as much metabolite as possible from the samples. After the extraction procedure, samples were transferred into vials and they were evaluated via time of flight mass spectrometry instrument.

RESULTS

Significantly, altered metabolites were identified by the fold analysis and Welch's test between the groups. 42 different compounds were annotated regarding to data-dependent acquisition method. Pathway analysis were also performed to understand the hazardous effect of methamphetamine on human body.

CONCLUSION

It has been reported that drug exposure may affect several metabolic pathways for amino acids, fats, energy metabolism and vitamins. An alternative bioinformatic model was also developed and validated in order to predict the chronic methamphetamine drug users in any criminal cases. This generated model passes the ROC curve analysis and permutation test and classify the controls and drug users correctly by evaluating the metabolic alterations between the groups.

Role of human plasma metabolites in prediabetes and type 2 diabetes from the IMI-DIRECT study

AIMS/HYPOTHESIS

Type 2 diabetes is a chronic condition that is caused by hyperglycaemia. Our aim was to characterise the metabolomics to find their association with the glycaemic spectrum and find a causal relationship between metabolites and type 2 diabetes.

METHODS

As part of the Innovative Medicines Initiative - Diabetes Research on Patient Stratification (IMI-DIRECT) consortium, 3000 plasma samples were measured with the Biocrates AbsoluteIDQ p150 Kit and Metabolon analytics. A total of 911 metabolites (132 targeted metabolomics, 779 untargeted metabolomics) passed the quality control. Multivariable linear and logistic regression analysis estimates were calculated from the concentration/peak areas of each metabolite as an explanatory variable and the glycaemic status as a dependent variable. This analysis was adjusted for age, sex, BMI, study centre in the basic model, and additionally for alcohol, smoking, BP, fasting HDL-cholesterol and fasting triacylglycerol in the full model. Statistical significance was Bonferroni corrected throughout. Beyond associations, we investigated the mediation effect and causal effects for which causal mediation test and two-sample Mendelian randomisation (2SMR) methods were used, respectively.

RESULTS

In the targeted metabolomics, we observed four (15), 34 (99) and 50 (108) metabolites (number of metabolites observed in untargeted metabolomics appear in parentheses) that were significantly different when comparing normal glucose regulation vs impaired glucose regulation/prediabetes, normal glucose regulation vs type 2 diabetes, and impaired glucose regulation vs type 2 diabetes, respectively. Significant metabolites were mainly branched-chain amino acids (BCAAs), with some derivatised BCAAs, lipids, xenobiotics and a few unknowns. Metabolites such as lysophosphatidylcholine a C17:0, sum of hexoses, amino acids from BCAA metabolism (including leucine, isoleucine, valine, N-lactoylvaline, N-lactoylleucine and formiminoglutamate) and lactate, as well as an unknown metabolite (X-24295), were associated with HbA1c progression rate and were significant mediators of type 2 diabetes from baseline to 18 and 48 months of follow-up. 2SMR was used to estimate the causal effect of an exposure on an outcome using summary statistics from UK Biobank genome-wide association studies. We found that type 2 diabetes had a causal effect on the levels of three metabolites (hexose, glutamate and caproate [fatty acid (FA) 6:0]), whereas lipids such as specific phosphatidylcholines (PCs) (namely PC aa C36:2, PC aa C36:5, PC ae C36:3 and PC ae C34:3) as well as the two n-3 fatty acids stearidonate (18:4n3) and docosapentaenoate (22:5n3) potentially had a causal role in the development of type 2 diabetes.

CONCLUSIONS/INTERPRETATION

Our findings identify known BCAAs and lipids, along with novel N-lactoyl-amino acid metabolites, significantly associated with prediabetes and diabetes, that mediate the effect of diabetes from baseline to follow-up (18 and 48 months). Causal inference using genetic variants shows the role of lipid metabolism and n-3 fatty acids as being causal for metabolite-to-type 2 diabetes whereas the sum of hexoses is causal for type 2 diabetes-to-metabolite. Identified metabolite markers are useful for stratifying individuals based on their risk progression and should enable targeted interventions.

Identification of a Novel Biomarker Panel for Breast Cancer Screening

Breast cancer remains a major public health concern, and early detection is crucial for improving survival rates. Metabolomics offers the potential to develop non-invasive screening and diagnostic tools based on metabolic biomarkers. However, the inherent complexity of metabolomic datasets and the high dimensionality of biomarkers complicates the identification of diagnostically relevant features, with multiple studies demonstrating limited consensus on the specific metabolites involved. Unlike previous studies that rely on singular feature selection techniques such as Partial Least Square (PLS) or LASSO regression, this research combines supervised and unsupervised machine learning methods with random sampling strategies, offering a more robust and interpretable approach to feature selection. This study aimed to identify a parsimonious and robust set of biomarkers for breast cancer diagnosis using metabolomics data. Plasma samples from 185 breast cancer patients and 53 controls (from the Cooperative Human Tissue Network, USA) were analyzed. This study also overcomes the common issue of dataset imbalance by using propensity score matching (PSM), which ensures reliable comparisons between cancer and control groups. We employed Univariate Naïve Bayes, L2-regularized Support Vector Classifier (SVC), Principal Component Analysis (PCA), and feature engineering techniques to refine and select the most informative features. Our best-performing feature set comprised 11 biomarkers, including 9 metabolites (SM(OH) C22:2, SM C18:0, C0, C3OH, C14:2OH, C16:2OH, LysoPC a C18:1, PC aa C36:0 and Asparagine), a metabolite ratio (Kynurenine-to-Tryptophan), and 1 demographic variable (Age), achieving an area under the ROC curve (AUC) of 98%. These results demonstrate the potential for a robust, cost-effective, and non-invasive breast cancer screening and diagnostic tool, offering significant clinical value for early detection and personalized patient management.

BMA-based Mendelian randomization identifies blood metabolites as causal candidates in pregnancy-induced hypertension

Pregnancy-induced hypertension (PIH), a prominent determinant of maternal mortality and morbidity worldwide, is hindered by the absence of efficacious biomarkers for early diagnosis, contributing to suboptimal outcomes. Here, we explored potential causal relationships between blood metabolites and the risk of PIH using Mendelian randomization (MR). We employed a two-sample univariable MR approach to empirically estimate the causal relationships between 249 circulating metabolites and PIH. Inverse variance weighted, MR-egger, weight median, simple mode, and weighted mode methods were used for causal estimates. The exposure-to-outcome directionality was confirmed with the MR Steiger test. The Bayesian model averaging MR (MR-BMA) method was applied to detect the predominant causal metabolic traits with alignment for pleiotropy effects. In the primary analysis, analyzing 249 metabolites, we identified 25 causally linked to PIH, including 11 lipid-related traits and 6 associated with fatty acid (un)saturation. Importantly, MR-BMA analyses corroborated the total concentration of branched-chain amino acids(total-BCAA) to be the highest rank causal metabolite, followed by leucine (Leu), phospholipids to total lipids ratio in medium LDL (M-LDL-PL-pct), and Val (all P < 0.05). The directionality of causality predicted by univariable MR and MR-BMA for these metabolites remained consistent. This study highlights the causal connection between metabolites and PIH risk. It highlighted BCAAs as the strongest causal candidates warranting further investigation. Since PIH typically occurs in the second and third trimesters, extending these findings could inform earlier strategies to reduce its risk. Directed acyclic graph of the MR framework investigating the causal relationship between metabolites and PIH. MR: Mendelian randomization; GIVs: genetic instrument variables; SNPs: single-nucleotide polymorphism; IVW: inverse variance weighted; WM: weighted median; PIH: pregnancy-induced hypertension; SM: significant metabolite; MR-BMA: Bayesian model averaging MR.

Metabolomics and mitochondrial dysfunction in cardiometabolic disease

Circulating metabolites are indicators of systemic metabolic dysfunction and can be detected through contemporary techniques in metabolomics. These metabolites are involved in numerous mitochondrial metabolic processes including glycolysis, fatty acid β-oxidation, and amino acid catabolism, and changes in the abundance of these metabolites is implicated in the pathogenesis of cardiometabolic diseases (CMDs). Epigenetic regulation and direct metabolite-protein interactions modulate metabolism, both within cells and in the circulation. Dysfunction of multiple mitochondrial components stemming from mitochondrial DNA mutations are implicated in disease pathogenesis. This review will summarize the current state of knowledge regarding: i) the interactions between metabolites found within the mitochondrial environment during CMDs, ii) various metabolites' effects on cellular and systemic function, iii) how harnessing the power of metabolomic analyses represents the next frontier of precision medicine, and iv) how these concepts integrate to expand the clinical potential for translational cardiometabolic medicine.

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.

Metabolic Signatures Elucidate the Effect of Body Mass Index on Type 2 Diabetes

Obesity plays an important role in the development of insulin resistance and diabetes, but the molecular mechanism that links obesity and diabetes is still not completely understood. Here, we used 146 targeted metabolomic profiles from the German KORA FF4 cohort consisting of 1715 participants and associated them with obesity and type 2 diabetes. In the basic model, 83 and 51 metabolites were significantly associated with body mass index (BMI) and T2D, respectively. Those metabolites are branched-chain amino acids, acylcarnitines, lysophospholipids, or phosphatidylcholines. In the full model, 42 and 3 metabolites were significantly associated with BMI and T2D, respectively, and replicate findings in the previous studies. Sobel mediation testing suggests that the effect of BMI on T2D might be mediated via lipids such as sphingomyelin (SM) C16:1, SM C18:1 and diacylphosphatidylcholine (PC aa) C38:3. Moreover, mendelian randomization suggests a causal relationship that BMI causes the change of SM C16:1 and PC aa C38:3, and the change of SM C16:1, SM C18:1, and PC aa C38:3 contribute to T2D incident. Biological pathway analysis in combination with genetics and mice experiments indicate that downregulation of sphingolipid or upregulation of phosphatidylcholine metabolism is a causal factor in early-stage T2D pathophysiology. Our findings indicate that metabolites like SM C16:1, SM C18:1, and PC aa C38:3 mediate the effect of BMI on T2D and elucidate their role in obesity related T2D pathologies.

[The promise of omics in the precision medicine era]

The rise of omics technologies that simultaneously measure thousands of molecules in a complex biological sample represents the core of systems biology. These technologies have profoundly impacted biomarkers and therapeutic targets discovery in the precision medicine era. Systems biology aims to perform a systematic probing of complex interactions in biological systems. Powered by high-throughput omics technologies and high-performance computing, systems biology provides relevant, resolving, and multi-scale overviews from cells to populations. Precision medicine takes advantage of these conceptual and technological developments and is based on two main pillars: the generation of multimodal data and their subsequent modeling. High-throughput omics technologies enable the comprehensive and holistic extraction of biological information, while computational capabilities enable multidimensional modeling and, as a result, offer an intuitive and intelligible visualization. Despite their promise, translating these technologies into clinically actionable tools has been slow. In this contribution, we present the most recent multi-omics data generation and analysis strategies and their clinical deployment in the post-genomic era. Furthermore, medical application challenges of omics-based biomarkers are discussed.

The crosstalk of the human microbiome in breast and colon cancer: A metabolomics analysis

The human microbiome's role in colon and breast cancer is described in this review. Understanding how the human microbiome and metabolomics interact with breast and colon cancer is the chief area of this study. First, the role of the gut and distal microbiome in breast and colon cancer is investigated, and the direct relationship between microbial dysbiosis and breast and colon cancer is highlighted. This work also focuses on the many metabolomic techniques used to locate prospective biomarkers, make an accurate diagnosis, and research new therapeutic targets for cancer treatment. This review clarifies the influence of anti-tumor medications on the microbiota and the proactive measures that can be taken to treat cancer using a variety of therapies, including radiotherapy, chemotherapy, next-generation biotherapeutics, gene-based therapy, integrated omics technology, and machine learning.

Metabolomic profiling in small vessel disease identifies multiple associations with disease severity

Cerebral small vessel disease is a major cause of vascular cognitive impairment and dementia. There are few treatments, largely reflecting limited understanding of the underlying pathophysiology. Metabolomics can be used to identify novel risk factors to better understand pathogenesis and to predict disease progression and severity. We analysed data from 624 patients with symptomatic cerebral small vessel disease from two prospective cohort studies. Serum samples were collected at baseline and patients underwent MRI scans and cognitive testing at regular intervals with up to 14 years of follow-up. Using ultra-performance liquid chromatography-mass spectrometry and nuclear magnetic resonance spectroscopy, we obtained metabolic and lipidomic profiles from 369 annotated metabolites and 54 764 unannotated features and examined their association with respect to disease severity, assessed using MRI small vessel disease markers, cognition and future risk of all-cause dementia. Our analysis identified 28 metabolites that were significantly associated with small vessel disease imaging markers and cognition. Decreased levels of multiple glycerophospholipids and sphingolipids were associated with increased small vessel disease load as evidenced by higher white matter hyperintensity volume, lower mean diffusivity normalized peak height, greater brain atrophy and impaired cognition. Higher levels of creatine, FA(18:2(OH)) and SM(d18:2/24:1) were associated with increased lacune count, higher white matter hyperintensity volume and impaired cognition. Lower baseline levels of carnitines and creatinine were associated with higher annualized change in peak width of skeletonized mean diffusivity, and 25 metabolites, including lipoprotein subclasses, amino acids and xenobiotics, were associated with future dementia incidence. Our results show multiple distinct metabolic signatures that are associated with imaging markers of small vessel disease, cognition and conversion to dementia. Further research should assess causality and the use of metabolomic screening to improve the ability to predict future disease severity and dementia risk in small vessel disease. The metabolomic profiles may also provide novel insights into disease pathogenesis and help identify novel treatment approaches.

Genetic Architecture of Untargeted Lipidomics in Cardiometabolic-Disease Patients Combines Strong Polygenic Control and Pleiotropy

Analysis of the genetic control of small metabolites provides powerful information on the regulation of the endpoints of genome expression. We carried out untargeted liquid chromatography−high-resolution mass spectrometry in 273 individuals characterized for pathophysiological elements of the cardiometabolic syndrome. We quantified 3013 serum lipidomic features, which we used in both genome-wide association studies (GWAS), using a panel of over 2.5 M imputed single-nucleotide polymorphisms (SNPs), and metabolome-wide association studies (MWAS) with phenotypes. Genetic analyses showed that 926 SNPs at 551 genetic loci significantly (q-value < 10−8) regulate the abundance of 74 lipidomic features in the group, with evidence of monogenic control for only 22 of these. In addition to this strong polygenic control of serum lipids, our results underscore instances of pleiotropy, when a single genetic locus controls the abundance of several distinct lipid features. Using the LIPID MAPS database, we assigned putative lipids, predominantly fatty acyls and sterol lipids, to 77% of the lipidome signals mapped to the genome. We identified significant correlations between lipids and clinical and biochemical phenotypes. These results demonstrate the power of untargeted lipidomic profiling for high-density quantitative molecular phenotyping in human-genetic studies and illustrate the complex genetic control of lipid metabolism.

Metabolomic Fingerprinting for the Detection of Early-Stage Lung Cancer: From the Genome to the Metabolome

The five-year survival rate of lung cancer patients is very low, mainly because most newly diagnosed patients present with locally advanced or metastatic disease. Therefore, early diagnosis is key to the successful treatment and management of lung cancer. Unfortunately, early detection methods of lung cancer are not ideal. In this brief review, we described early detection methods such as chest X-rays followed by bronchoscopy, sputum analysis followed by cytological analysis, and low-dose computed tomography (LDCT). In addition, we discussed the potential of metabolomic fingerprinting, compared to that of other biomarkers, including molecular targets, as a low-cost, high-throughput blood-based test that is both feasible and affordable for early-stage lung cancer screening of at-risk populations. Accordingly, we proposed a paradigm shift to metabolomics as an alternative to molecular and proteomic-based markers in lung cancer screening, which will enable blood-based routine testing and be accessible to those patients at the highest risk for lung cancer.

FLEXiGUT: Rationale for exposomics associations with chronic low-grade gut inflammation

FLEXiGUT is the first large-scale exposomics study focused on chronic low-grade inflammation. It aims to characterize human life course environmental exposure to assess and validate its impact on gut inflammation and related biological processes and diseases. The cumulative influences of environmental and food contaminants throughout the lifespan on certain biological responses related to chronic gut inflammation will be investigated in two Flemish prospective cohorts, namely the "ENVIRONAGE birth cohort", which provides follow-up from gestation to early childhood, and the "Flemish Gut Flora Project longitudinal cohort", a cohort of adults. The exposome will be characterised through biomonitoring of legacy and emerging contaminants, mycotoxins and markers of air pollution, by analysing the available metadata on nutrition, location and activity, and by applying state-of-the-art -omics techniques, including metagenomics, metabolomics and DNA adductomics, as well as the assessment of telomere length and measurement of inflammatory markers, to encompass both exposure and effect. Associations between exposures and health outcomes will be uncovered using an integrated -omics data analysis framework comprising data exploration, pre-processing, dimensionality reduction and data mining, combined with machine learning-based pathway analysis approaches. This is expected to lead to a more profound insight in mechanisms underlying disease progression (e.g. metabolic disorders, food allergies, gastrointestinal cancers) and/or accelerated biological ageing.

Leukocyte proteomics coupled with serum metabolomics identifies novel biomarkers and abnormal amino acid metabolism in Kawasaki disease

Kawasaki disease (KD) is a systemic vasculitis that can lead to severe cardiovascular complications, whereas the development and clinical usage of specific biomarkers might help diagnose KD and avoid certain complications. To this end, the molecular profiles of acute KD patients with coronary artery lesions (CAL) were first investigated through leukocyte proteomics and serum metabolomics assays. A total of 269 differentially abundant proteins and 35 differentially abundant metabolites with the top fold-changed levels were identified in acute KD patients compared to those in the healthy controls. Among them, several highly promising candidate marker proteins and metabolites indicative of KD progression were further analysed, such as the increased proteins ALPL, NAMPT, and S100P, as well as the decreased proteins C1QB and apolipoprotein family members. Moreover, metabolites, including succinic acid, dGMP, hyaluronic acid, L-tryptophan, propionylcarnitine, inosine, and phosphorylcholine, were found to be highly accurate at distinguishing between KD patients and healthy controls. Interestingly, the abnormal expression levels of a distinct set of proteins and metabolites in acute KD patients can be restored to normal levels upon intravenous immunoglobulin (IVIG) treatment. Overall, this work has revealed novel biomarkers and abnormal amino-acid metabolism as a prominent feature involved in KD patients with CAL. SIGNIFICANCE: KD is frequently concomitant with the development of life-threatening coronary vasculitis. Here, the profiles of leukocyte proteomics and serum metabolomics in acute KD patients with CALs were first investigated, and several hub molecules identified here could be used as supplemental biomarkers for KD diagnosis. Moreover, the metabolomic abnormalities especially the amino acids are particularly prominent in KD patients. Interestingly, the abnormal expression levels of a distinct set of proteins and metabolites in acute KD patients can be restored to normal levels upon IVIG treatment. Therefore, these findings might help understand the IVIG activities and also the underlying mechanisms of IVIG-resistant patients, thereby providing a new perspective for the exploration of mechanisms related to KD.

Rare genetic variants affecting urine metabolite levels link population variation to inborn errors of metabolism

Metabolite levels in urine may provide insights into genetic mechanisms shaping their related pathways. We therefore investigate the cumulative contribution of rare, exonic genetic variants on urine levels of 1487 metabolites and 53,714 metabolite ratios among 4864 GCKD study participants. Here we report the detection of 128 significant associations involving 30 unique genes, 16 of which are known to underlie inborn errors of metabolism. The 30 genes are strongly enriched for shared expression in liver and kidney (odds ratio = 65, p-FDR = 3e-7), with hepatocytes and proximal tubule cells as driving cell types. Use of UK Biobank whole-exome sequencing data links genes to diseases connected to the identified metabolites. In silico constraint-based modeling of gene knockouts in a virtual whole-body, organ-resolved metabolic human correctly predicts the observed direction of metabolite changes, highlighting the potential of linking population genetics to modeling. Our study implicates candidate variants and genes for inborn errors of metabolism.

Urinary metabolite quantitative trait loci in children and their interaction with dietary factors

Human metabolism is influenced by genetic and environmental factors. Previous studies have identified over 23 loci associated with more than 26 urine metabolites levels in adults, which are known as urinary metabolite quantitative trait loci (metabQTLs). The aim of the present study is the identification for the first time of urinary metabQTLs in children and their interaction with dietary patterns. Association between genome-wide genotyping data and 44 urine metabolite levels measured by proton nuclear magnetic resonance spectroscopy was tested in 996 children from the Human Early Life Exposome project. Twelve statistically significant urine metabQTLs were identified, involving 11 unique loci and 10 different metabolites. Comparison with previous findings in adults revealed that six metabQTLs were already known, and one had been described in serum and three were involved the same locus as other reported metabQTLs but had different urinary metabolites. The remaining two metabQTLs represent novel urine metabolite-locus associations, which are reported for the first time in this study [single nucleotide polymorphism (SNP) rs12575496 for taurine, and the missense SNP rs2274870 for 3-hydroxyisobutyrate]. Moreover, it was found that urinary taurine levels were affected by the combined action of genetic variation and dietary patterns of meat intake as well as by the interaction of this SNP with beverage intake dietary patterns. Overall, we identified 12 urinary metabQTLs in children, including two novel associations. While a substantial part of the identified loci affected urinary metabolite levels both in children and in adults, the metabQTL for taurine seemed to be specific to children and interacted with dietary patterns.

Genetics and Not Shared Environment Explains Familial Resemblance in Adult Metabolomics Data

Metabolites are small molecules involved in cellular metabolism where they act as reaction substrates or products. The term 'metabolomics' refers to the comprehensive study of these molecules. The concentrations of metabolites in biological tissues are under genetic control, but this is limited by environmental factors such as diet. In adult mono- and dizygotic twin pairs, we estimated the contribution of genetic and shared environmental influences on metabolite levels by structural equation modeling and tested whether the familial resemblance for metabolite levels is mainly explained by genetic or by environmental factors that are shared by family members. Metabolites were measured across three platforms: two based on proton nuclear magnetic resonance techniques and one employing mass spectrometry. These three platforms comprised 237 single metabolic traits of several chemical classes. For the three platforms, metabolites were assessed in 1407, 1037 and 1116 twin pairs, respectively. We carried out power calculations to establish what percentage of shared environmental variance could be detected given these sample sizes. Our study did not find evidence for a systematic contribution of shared environment, defined as the influence of growing up together in the same household, on metabolites assessed in adulthood. Significant heritability was observed for nearly all 237 metabolites; significant contribution of the shared environment was limited to 6 metabolites. The top quartile of the heritability distribution was populated by 5 of the 11 investigated chemical classes. In this quartile, metabolites of the class lipoprotein were significantly overrepresented, whereas metabolites of classes glycerophospholipids and glycerolipids were significantly underrepresented.

Updates in Rhea: SPARQLing biochemical reaction data

Rhea (http://www.rhea-db.org) is a comprehensive and non-redundant resource of over 11 000 expert-curated biochemical reactions that uses chemical entities from the ChEBI ontology to represent reaction participants. Originally designed as an annotation vocabulary for the UniProt Knowledgebase (UniProtKB), Rhea also provides reaction data for a range of other core knowledgebases and data repositories including ChEBI and MetaboLights. Here we describe recent developments in Rhea, focusing on a new resource description framework representation of Rhea reaction data and an SPARQL endpoint (https://sparql.rhea-db.org/sparql) that provides access to it. We demonstrate how federated queries that combine the Rhea SPARQL endpoint and other SPARQL endpoints such as that of UniProt can provide improved metabolite annotation and support integrative analyses that link the metabolome through the proteome to the transcriptome and genome. These developments will significantly boost the utility of Rhea as a means to link chemistry and biology for a more holistic understanding of biological systems and their function in health and disease.

mGWAS: next generation genetic prediction in kidney disease

A recent metabolite genome-wide association study (mGWAS) investigated the relationship between genetic factors and the urine metabolome in kidney disease. The findings demonstrate that mGWAS hold promise for identifying novel genetic factors involved in adsorption, distribution, metabolism and excretion of metabolites and pharmaceuticals, as well as biomarkers for disease progression.

Genetic studies of urinary metabolites illuminate mechanisms of detoxification and excretion in humans

The kidneys integrate information from continuous systemic processes related to the absorption, distribution, metabolism and excretion (ADME) of metabolites. To identify underlying molecular mechanisms, we performed genome-wide association studies of the urinary concentrations of 1,172 metabolites among 1,627 patients with reduced kidney function. The 240 unique metabolite-locus associations (metabolite quantitative trait loci, mQTLs) that were identified and replicated highlight novel candidate substrates for transport proteins. The identified genes are enriched in ADME-relevant tissues and cell types, and they reveal novel candidates for biotransformation and detoxification reactions. Fine mapping of mQTLs and integration with single-cell gene expression permitted the prioritization of causal genes, functional variants and target cell types. The combination of mQTLs with genetic and health information from 450,000 UK Biobank participants illuminated metabolic mediators, and hence, novel urinary biomarkers of disease risk. This comprehensive resource of genetic targets and their substrates is informative for ADME processes in humans and is relevant to basic science, clinical medicine and pharmaceutical research.

Metabolomics for Investigating Physiological and Pathophysiological Processes

Metabolomics uses advanced analytical chemistry techniques to enable the high-throughput characterization of metabolites from cells, organs, tissues, or biofluids. The rapid growth in metabolomics is leading to a renewed interest in metabolism and the role that small molecule metabolites play in many biological processes. As a result, traditional views of metabolites as being simply the "bricks and mortar" of cells or just the fuel for cellular energetics are being upended. Indeed, metabolites appear to have much more varied and far more important roles as signaling molecules, immune modulators, endogenous toxins, and environmental sensors. This review explores how metabolomics is yielding important new insights into a number of important biological and physiological processes. In particular, a major focus is on illustrating how metabolomics and discoveries made through metabolomics are improving our understanding of both normal physiology and the pathophysiology of many diseases. These discoveries are yielding new insights into how metabolites influence organ function, immune function, nutrient sensing, and gut physiology. Collectively, this work is leading to a much more unified and system-wide perspective of biology wherein metabolites, proteins, and genes are understood to interact synergistically to modify the actions and functions of organelles, organs, and organisms.

Metabolomics in genetic testing

Metabolomics is an intriguing field of study providing a new readout of the biochemical activities taking place at the moment of sampling within a subject's biofluid or tissue. Metabolite concentrations are influenced by several factors including disease, environment, drugs, diet and, importantly, genetics. Metabolomics signatures, which describe a subject's phenotype, are useful for disease diagnosis and prognosis, as well as for predicting and monitoring the effectiveness of treatments. Metabolomics is conventionally divided into targeted (i.e., the quantitative analysis of a predetermined group of metabolites) and untargeted studies (i.e., analysis of the complete set of small-molecule metabolites contained in a biofluid without a pre-imposed metabolites-selection). Both approaches have demonstrated high value in the investigation and understanding of several monogenic and multigenic conditions. Due to low costs per sample and relatively short analysis times, metabolomics can be a useful and robust complement to genetic sequencing.

Genome-Wide Association Studies of Metabolite Concentrations (mGWAS): Relevance for Nephrology

Metabolites are small molecules that are intermediates or products of metabolism, many of which are freely filtered by the kidneys. In addition, the kidneys have a central role in metabolite anabolism and catabolism, as well as in active metabolite reabsorption and/or secretion during tubular passage. This review article illustrates how the coupling of genomics and metabolomics in genome-wide association analyses of metabolites can be used to illuminate mechanisms underlying human metabolism, with a special focus on insights relevant to nephrology. First, genetic susceptibility loci for reduced kidney function and chronic kidney disease (CKD) were reviewed systematically for their associations with metabolite concentrations in metabolomics studies of blood and urine. Second, kidney function and CKD-associated metabolites reported from observational studies were interrogated for metabolite-associated genetic variants to generate and discuss complementary insights. Finally, insights originating from the simultaneous study of both blood and urine or by modeling intermetabolite relationships are summarized. We also discuss methodologic questions related to the study of metabolite concentrations in urine as well as among CKD patients. In summary, genome-wide association analyses of metabolites using metabolite concentrations quantified from blood and/or urine are a promising avenue of research to illuminate physiological and pathophysiological functions of the kidney.

Association between serum haptoglobin and carotid arterial functions: usefulness of a targeted metabolomics approach

Background

Serum haptoglobin (Hp) has been closely associated with cardio-cerebrovascular diseases. We investigated a metabolic profile associated with circulating Hp and carotid arterial functions via a targeted metabolomics approach to provide insight into potential mechanisms.

Methods

A total of 240 participants, including 120 patients with type 2 diabetes mellitus (T2DM) and 120 non-diabetes mellitus (non-DM) subjects were recruited in this study. Targeted metabolic profiles of serum metabolites were determined using an AbsoluteIDQ™ p180 Kit (BIOCRATES Life Sciences AG, Innsbruck, Austria). Ultrasound of the bilateral common carotid artery was used to measure intima-media thickness and inter-adventitial diameter. Serum Hp levels were tested by enzyme-linked immunosorbent assay.

Results

Serum Hp levels in T2DM patients and non-DM subjects were 103.40 (72.46, 131.99) mg/dL and 100.20 (53.99, 140.66) mg/dL, respectively. Significant differences of 19 metabolites and 17 metabolites were found among serum Hp tertiles in T2DM patients and non-DM subjects, respectively (P < 0.05). Of these, phosphatidylcholine acyl-alkyl C32:2 (PC ae C32:2) was the common metabolite observed in two populations, which was associated with the serum Hp groups and lipid traits (P < 0.05). Furthermore, the metabolite ratios of two acidic amino acids, including aspartate to PC ae C32:2 (Asp/PC ae C32:2) and glutamate to PC ae C32:2 (Glu/PC ae C32:2) were correlated with serum Hp, carotid arterial functions and other biochemical index in both populations significantly (P < 0.05).

Conclusions

Targeted metabolomics analyses might provide a new insight into the potential mechanisms underlying the association between serum Hp and carotid arterial functions.

Electronic supplementary material

The online version of this article (10.1186/s12933-019-0808-2) contains supplementary material, which is available to authorized users.

Advances in metabolome information retrieval: turning chemistry into biology. Part II: biological information recovery

This work reports the second part of a review intending to give the state of the art of major metabolic phenotyping strategies. It particularly deals with inherent advantages and limits regarding data analysis issues and biological information retrieval tools along with translational challenges. This Part starts with introducing the main data preprocessing strategies of the different metabolomics data. Then, it describes the main data analysis techniques including univariate and multivariate aspects. It also addresses the challenges related to metabolite annotation and characterization. Finally, functional analysis including pathway and network strategies are discussed. The last section of this review is devoted to practical considerations and current challenges and pathways to bring metabolomics into clinical environments.

Whole-genome sequencing study of serum peptide levels: the Atherosclerosis Risk in Communities study

Oligopeptides are important markers of protein metabolism, as they are cleaved from larger polypeptides and proteins. Genetic association studies may help elucidate their origin and function. In 1,552 European Americans and 1,872 African Americans of the Atherosclerosis Risk in Communities study, we performed whole-genome and whole-exome sequencing and measured serum levels of 25 peptides. Common variants (minor allele frequency > 5%) were analysed individually. We grouped low-frequency variants (minor allele frequency ≤ 5%) by a genome-wide sliding window using region-based aggregate tests. Furthermore, low-frequency regulatory variants were grouped by gene, as were functional coding variants. All analyses were performed separately in each ancestry group and then meta-analysed. We identified 22 common variant associations with peptide levels (P-value < 4.2 × 10-10), including 16 novel gene-peptide pairs. Notably, variants in kinin-kallikrein genes KNG1, F12, KLKB1, and ACE were associated with several different peptides. Variants in KLKB1 and ACE were associated with a fragment of complement component 3f. Both common variants and low-frequency coding variants in CPN1 were associated with a fibrinogen cleavage peptide. Four sliding windows were significantly associated with peptide levels (P-value < 4.2 × 10-10). Our results highlight the importance of the kinin-kallikrein system in the regulation of serum peptide levels, strengthen the evidence for a broad link between the kinin-kallikrein and complement systems, and suggest a role of CPN1 in the conversion of fibrinogen to fibrin.

Metabolite ratios as potential biomarkers for type 2 diabetes: a DIRECT study

AIMS/HYPOTHESIS

Circulating metabolites have been shown to reflect metabolic changes during the development of type 2 diabetes. In this study we examined the association of metabolite levels and pairwise metabolite ratios with insulin responses after glucose, glucagon-like peptide-1 (GLP-1) and arginine stimulation. We then investigated if the identified metabolite ratios were associated with measures of OGTT-derived beta cell function and with prevalent and incident type 2 diabetes.

METHODS

We measured the levels of 188 metabolites in plasma samples from 130 healthy members of twin families (from the Netherlands Twin Register) at five time points during a modified 3 h hyperglycaemic clamp with glucose, GLP-1 and arginine stimulation. We validated our results in cohorts with OGTT data (n = 340) and epidemiological case-control studies of prevalent (n = 4925) and incident (n = 4277) diabetes. The data were analysed using regression models with adjustment for potential confounders.

RESULTS

There were dynamic changes in metabolite levels in response to the different secretagogues. Furthermore, several fasting pairwise metabolite ratios were associated with one or multiple clamp-derived measures of insulin secretion (all p < 9.2 × 10^-7). These associations were significantly stronger compared with the individual metabolite components. One of the ratios, valine to phosphatidylcholine acyl-alkyl C32:2 (PC ae C32:2), in addition showed a directionally consistent positive association with OGTT-derived measures of insulin secretion and resistance (p ≤ 5.4 × 10^-3) and prevalent type 2 diabetes (OR_{Val_PC ae C32:2} 2.64 [β 0.97 ± 0.09], p = 1.0 × 10^-27). Furthermore, Val_PC ae C32:2 predicted incident diabetes independent of established risk factors in two epidemiological cohort studies (HR_{Val_PC ae C32:2} 1.57 [β 0.45 ± 0.06]; p = 1.3 × 10^-15), leading to modest improvements in the receiver operating characteristics when added to a model containing a set of established risk factors in both cohorts (increases from 0.780 to 0.801 and from 0.862 to 0.865 respectively, when added to the model containing traditional risk factors + glucose).

CONCLUSIONS/INTERPRETATION

In this study we have shown that the Val_PC ae C32:2 metabolite ratio is associated with an increased risk of type 2 diabetes and measures of insulin secretion and resistance. The observed effects were stronger than that of the individual metabolites and independent of known risk factors.

Metabomatching: Using genetic association to identify metabolites in proton NMR spectroscopy

A metabolome-wide genome-wide association study (mGWAS) aims to discover the effects of genetic variants on metabolome phenotypes. Most mGWASes use as phenotypes concentrations of limited sets of metabolites that can be identified and quantified from spectral information. In contrast, in an untargeted mGWAS both identification and quantification are forgone and, instead, all measured metabolome features are tested for association with genetic variants. While the untargeted approach does not discard data that may have eluded identification, the interpretation of associated features remains a challenge. To address this issue, we developed metabomatching to identify the metabolites underlying significant associations observed in untargeted mGWASes on proton NMR metabolome data. Metabomatching capitalizes on genetic spiking, the concept that because metabolome features associated with a genetic variant tend to correspond to the peaks of the NMR spectrum of the underlying metabolite, genetic association can allow for identification. Applied to the untargeted mGWASes in the SHIP and CoLaus cohorts and using 180 reference NMR spectra of the urine metabolome database, metabomatching successfully identified the underlying metabolite in 14 of 19, and 8 of 9 associations, respectively. The accuracy and efficiency of our method make it a strong contender for facilitating or complementing metabolomics analyses in large cohorts, where the availability of genetic, or other data, enables our approach, but targeted quantification is limited.

Quantitative Serum Nuclear Magnetic Resonance Metabolomics in Large-Scale Epidemiology: A Primer on -Omic Technologies

Detailed metabolic profiling in large-scale epidemiologic studies has uncovered novel biomarkers for cardiometabolic diseases and clarified the molecular associations of established risk factors. A quantitative metabolomics platform based on nuclear magnetic resonance spectroscopy has found widespread use, already profiling over 400,000 blood samples. Over 200 metabolic measures are quantified per sample; in addition to many biomarkers routinely used in epidemiology, the method simultaneously provides fine-grained lipoprotein subclass profiling and quantification of circulating fatty acids, amino acids, gluconeogenesis-related metabolites, and many other molecules from multiple metabolic pathways. Here we focus on applications of magnetic resonance metabolomics for quantifying circulating biomarkers in large-scale epidemiology. We highlight the molecular characterization of risk factors, use of Mendelian randomization, and the key issues of study design and analyses of metabolic profiling for epidemiology. We also detail how integration of metabolic profiling data with genetics can enhance drug development. We discuss why quantitative metabolic profiling is becoming widespread in epidemiology and biobanking. Although large-scale applications of metabolic profiling are still novel, it seems likely that comprehensive biomarker data will contribute to etiologic understanding of various diseases and abilities to predict disease risks, with the potential to translate into multiple clinical settings.

Endocrinology Meets Metabolomics: Achievements, Pitfalls, and Challenges

The metabolome, although very dynamic, is sufficiently stable to provide specific quantitative traits related to health and disease. Metabolomics requires balanced use of state-of-the-art study design, chemical analytics, biostatistics, and bioinformatics to deliver meaningful answers to contemporary questions in human disease research. The technology is now frequently employed for biomarker discovery and for elucidating the mechanisms underlying endocrine-related diseases. Metabolomics has also enriched genome-wide association studies (GWAS) in this area by providing functional data. The contributions of rare genetic variants to metabolome variance and to the human phenotype have been underestimated until now.

Extending the MR-Egger method for multivariable Mendelian randomization to correct for both measured and unmeasured pleiotropy

Methods have been developed for Mendelian randomization that can obtain consistent causal estimates while relaxing the instrumental variable assumptions. These include multivariable Mendelian randomization, in which a genetic variant may be associated with multiple risk factors so long as any association with the outcome is via the measured risk factors (measured pleiotropy), and the MR‐Egger (Mendelian randomization‐Egger) method, in which a genetic variant may be directly associated with the outcome not via the risk factor of interest, so long as the direct effects of the variants on the outcome are uncorrelated with their associations with the risk factor (unmeasured pleiotropy). In this paper, we extend the MR‐Egger method to a multivariable setting to correct for both measured and unmeasured pleiotropy. We show, through theoretical arguments and a simulation study, that the multivariable MR‐Egger method has advantages over its univariable counterpart in terms of plausibility of the assumption needed for consistent causal estimation and power to detect a causal effect when this assumption is satisfied. The methods are compared in an applied analysis to investigate the causal effect of high‐density lipoprotein cholesterol on coronary heart disease risk. The multivariable MR‐Egger method will be useful to analyse high‐dimensional data in situations where the risk factors are highly related and it is difficult to find genetic variants specifically associated with the risk factor of interest (multivariable by design), and as a sensitivity analysis when the genetic variants are known to have pleiotropic effects on measured risk factors.

Metabolomic Biomarkers in the Progression to Type 1 Diabetes

Metabolomics is the snapshot of all detectable metabolites and lipids in biological materials and has potential in reflecting genetic and environmental factors contributing to the development of complex diseases, such as type 1 diabetes. The progression to seroconversion to development of type 1 diabetes has been studied using this technique, although in relatively small cohorts and at limited time points. Overall, three observations have been consistently reported; phospholipids at birth are lower in children developing type 1 diabetes early in childhood, methionine levels are lower in children at seroconversion, and triglycerides are increased at seroconversion and associated to microbiome diversity, indicating an association between the metabolome and microbiome in type 1 diabetes progression.

Non-targeted metabolomics combined with genetic analyses identifies bile acid synthesis and phospholipid metabolism as being associated with incident type 2 diabetes

AIMS/HYPOTHESIS

Identification of novel biomarkers for type 2 diabetes and their genetic determinants could lead to improved understanding of causal pathways and improve risk prediction.

METHODS

In this study, we used data from non-targeted metabolomics performed using liquid chromatography coupled with tandem mass spectrometry in three Swedish cohorts (Uppsala Longitudinal Study of Adult Men [ULSAM], n = 1138; Prospective Investigation of the Vasculature in Uppsala Seniors [PIVUS], n = 970; TwinGene, n = 1630). Metabolites associated with impaired fasting glucose (IFG) and/or prevalent type 2 diabetes were assessed for associations with incident type 2 diabetes in the three cohorts followed by replication attempts in the Cooperative Health Research in the Region of Augsburg (KORA) S4 cohort (n = 855). Assessment of the association of metabolite-regulating genetic variants with type 2 diabetes was done using data from a meta-analysis of genome-wide association studies.

RESULTS

Out of 5961 investigated metabolic features, 1120 were associated with prevalent type 2 diabetes and IFG and 70 were annotated to metabolites and replicated in the three cohorts. Fifteen metabolites were associated with incident type 2 diabetes in the four cohorts combined (358 events) following adjustment for age, sex, BMI, waist circumference and fasting glucose. Novel findings included associations of higher values of the bile acid deoxycholic acid and monoacylglyceride 18:2 and lower concentrations of cortisol with type 2 diabetes risk. However, adding metabolites to an existing risk score improved model fit only marginally. A genetic variant within the CYP7A1 locus, encoding the rate-limiting enzyme in bile acid synthesis, was found to be associated with lower concentrations of deoxycholic acid, higher concentrations of LDL-cholesterol and lower type 2 diabetes risk. Variants in or near SGPP1, GCKR and FADS1/2 were associated with diabetes-associated phospholipids and type 2 diabetes.

CONCLUSIONS/INTERPRETATION

We found evidence that the metabolism of bile acids and phospholipids shares some common genetic origin with type 2 diabetes.

ACCESS TO RESEARCH MATERIALS

Metabolomics data have been deposited in the Metabolights database, with accession numbers MTBLS93 (TwinGene), MTBLS124 (ULSAM) and MTBLS90 (PIVUS).

Omics-Based Strategies in Precision Medicine: Toward a Paradigm Shift in Inborn Errors of Metabolism Investigations

The rise of technologies that simultaneously measure thousands of data points represents the heart of systems biology. These technologies have had a huge impact on the discovery of next-generation diagnostics, biomarkers, and drugs in the precision medicine era. Systems biology aims to achieve systemic exploration of complex interactions in biological systems. Driven by high-throughput omics technologies and the computational surge, it enables multi-scale and insightful overviews of cells, organisms, and populations. Precision medicine capitalizes on these conceptual and technological advancements and stands on two main pillars: data generation and data modeling. High-throughput omics technologies allow the retrieval of comprehensive and holistic biological information, whereas computational capabilities enable high-dimensional data modeling and, therefore, accessible and user-friendly visualization. Furthermore, bioinformatics has enabled comprehensive multi-omics and clinical data integration for insightful interpretation. Despite their promise, the translation of these technologies into clinically actionable tools has been slow. In this review, we present state-of-the-art multi-omics data analysis strategies in a clinical context. The challenges of omics-based biomarker translation are discussed. Perspectives regarding the use of multi-omics approaches for inborn errors of metabolism (IEM) are presented by introducing a new paradigm shift in addressing IEM investigations in the post-genomic era.

Metabolomics enables precision medicine: "A White Paper, Community Perspective"

INTRODUCTION BACKGROUND TO METABOLOMICS

Metabolomics is the comprehensive study of the metabolome, the repertoire of biochemicals (or small molecules) present in cells, tissues, and body fluids. The study of metabolism at the global or "-omics" level is a rapidly growing field that has the potential to have a profound impact upon medical practice. At the center of metabolomics, is the concept that a person's metabolic state provides a close representation of that individual's overall health status. This metabolic state reflects what has been encoded by the genome, and modified by diet, environmental factors, and the gut microbiome. The metabolic profile provides a quantifiable readout of biochemical state from normal physiology to diverse pathophysiologies in a manner that is often not obvious from gene expression analyses. Today, clinicians capture only a very small part of the information contained in the metabolome, as they routinely measure only a narrow set of blood chemistry analytes to assess health and disease states. Examples include measuring glucose to monitor diabetes, measuring cholesterol and high density lipoprotein/low density lipoprotein ratio to assess cardiovascular health, BUN and creatinine for renal disorders, and measuring a panel of metabolites to diagnose potential inborn errors of metabolism in neonates.

OBJECTIVES OF WHITE PAPER—EXPECTED TREATMENT OUTCOMES AND METABOLOMICS ENABLING TOOL FOR PRECISION MEDICINE

We anticipate that the narrow range of chemical analyses in current use by the medical community today will be replaced in the future by analyses that reveal a far more comprehensive metabolic signature. This signature is expected to describe global biochemical aberrations that reflect patterns of variance in states of wellness, more accurately describe specific diseases and their progression, and greatly aid in differential diagnosis. Such future metabolic signatures will: (1) provide predictive, prognostic, diagnostic, and surrogate markers of diverse disease states; (2) inform on underlying molecular mechanisms of diseases; (3) allow for sub-classification of diseases, and stratification of patients based on metabolic pathways impacted; (4) reveal biomarkers for drug response phenotypes, providing an effective means to predict variation in a subject's response to treatment (pharmacometabolomics); (5) define a metabotype for each specific genotype, offering a functional read-out for genetic variants: (6) provide a means to monitor response and recurrence of diseases, such as cancers: (7) describe the molecular landscape in human performance applications and extreme environments. Importantly, sophisticated metabolomic analytical platforms and informatics tools have recently been developed that make it possible to measure thousands of metabolites in blood, other body fluids, and tissues. Such tools also enable more robust analysis of response to treatment. New insights have been gained about mechanisms of diseases, including neuropsychiatric disorders, cardiovascular disease, cancers, diabetes and a range of pathologies. A series of ground breaking studies supported by National Institute of Health (NIH) through the Pharmacometabolomics Research Network and its partnership with the Pharmacogenomics Research Network illustrate how a patient's metabotype at baseline, prior to treatment, during treatment, and post-treatment, can inform about treatment outcomes and variations in responsiveness to drugs (e.g., statins, antidepressants, antihypertensives and antiplatelet therapies). These studies along with several others also exemplify how metabolomics data can complement and inform genetic data in defining ethnic, sex, and gender basis for variation in responses to treatment, which illustrates how pharmacometabolomics and pharmacogenomics are complementary and powerful tools for precision medicine.

CONCLUSIONS KEY SCIENTIFIC CONCEPTS AND RECOMMENDATIONS FOR PRECISION MEDICINE

Our metabolomics community believes that inclusion of metabolomics data in precision medicine initiatives is timely and will provide an extremely valuable layer of data that compliments and informs other data obtained by these important initiatives. Our Metabolomics Society, through its "Precision Medicine and Pharmacometabolomics Task Group", with input from our metabolomics community at large, has developed this White Paper where we discuss the value and approaches for including metabolomics data in large precision medicine initiatives. This White Paper offers recommendations for the selection of state of-the-art metabolomics platforms and approaches that offer the widest biochemical coverage, considers critical sample collection and preservation, as well as standardization of measurements, among other important topics. We anticipate that our metabolomics community will have representation in large precision medicine initiatives to provide input with regard to sample acquisition/preservation, selection of optimal omics technologies, and key issues regarding data collection, interpretation, and dissemination. We strongly recommend the collection and biobanking of samples for precision medicine initiatives that will take into consideration needs for large-scale metabolic phenotyping studies.

Genome-Wide Association of the Laboratory-Based Nicotine Metabolite Ratio in Three Ancestries

Introduction:

Metabolic enzyme variation and other patient and environmental characteristics influence smoking behaviors, treatment success, and risk of related disease. Population-specific variation in metabolic genes contributes to challenges in developing and optimizing pharmacogenetic interventions. We applied a custom genome-wide genotyping array for addiction research (Smokescreen), to three laboratory-based studies of nicotine metabolism with oral or venous administration of labeled nicotine and cotinine, to model nicotine metabolism in multiple populations. The trans-3′-hydroxycotinine/cotinine ratio, the nicotine metabolite ratio (NMR), was the nicotine metabolism measure analyzed.

Methods:

Three hundred twelve individuals of self-identified European, African, and Asian American ancestry were genotyped and included in ancestry-specific genome-wide association scans (GWAS) and a meta-GWAS analysis of the NMR. We modeled natural-log transformed NMR with covariates: principal components of genetic ancestry, age, sex, body mass index, and smoking status.

Results:

African and Asian American NMRs were statistically significantly (P values ≤ 5E-5) lower than European American NMRs. Meta-GWAS analysis identified 36 genome-wide significant variants over a 43 kilobase pair region at CYP2A6 with minimum P = 2.46E-18 at rs12459249, proximal to CYP2A6. Additional minima were located in intron 4 (rs56113850, P = 6.61E-18) and in the CYP2A6-CYP2A7 intergenic region (rs34226463, P = 1.45E-12). Most (34/36) genome-wide significant variants suggested reduced CYP2A6 activity; functional mechanisms were identified and tested in knowledge-bases. Conditional analysis resulted in intergenic variants of possible interest (P values < 5E-5).

Conclusions:

This meta-GWAS of the NMR identifies CYP2A6 variants, replicates the top-ranked single nucleotide polymorphism from a recent Finnish meta-GWAS of the NMR, identifies functional mechanisms, and provides pan-continental population biomarkers for nicotine metabolism.

Implications:

This multiple ancestry meta-GWAS of the laboratory study-based NMR provides novel evidence and replication for genome-wide association of CYP2A6 single nucleotide and insertion–deletion polymorphisms. We identify three regions of genome-wide significance: proximal, intronic, and distal to CYP2A6. We replicate the top-ranking single nucleotide polymorphism from a recent GWAS of the NMR in Finnish smokers, identify a functional mechanism for this intronic variant from in silico analyses of RNA-seq data that is consistent with CYP2A6 expression measured in postmortem lung and liver, and provide additional support for the intergenic region between CYP2A6 and CYP2A7.

Mendelian randomization to assess causal effects of blood lipids on coronary heart disease: lessons from the past and applications to the future

PURPOSE OF REVIEW

Mendelian randomization is a technique for judging the causal impact of a risk factor on an outcome from observational data using genetic variants. Although evidence from Mendelian randomization for the effects of major lipids and lipoproteins on coronary heart disease (CHD) risk has been around for a relatively long time, new data resources and new methodological approaches have given fresh insight into these relationships. The lessons from these analyses are likely to be highly relevant when it comes to lipidomics and the analyses of lipid subspecies whose biology is less well understood.

RECENT FINDINGS

Although analyses of low-density lipoprotein cholesterol and lipoprotein(a) are unambiguous as there are genetic variants that associate exclusively with these risk factors and have well understood biology, analyses for triglycerides, and high-density lipoprotein cholesterol (HDL-c) are less clear. For example, a subset of genetic variants having specific associations with HDL-c are not associated with CHD risk, but an allele score including all variants associated with HDL-c does associate with CHD risk. Recently developed methods, such as multivariable Mendelian randomization, Mendelian randomization-Egger, and a weighted median method, suggest that the relationship between HDL-c and CHD risk is null, thus confirming experimental evidence.

SUMMARY

Robust methods for Mendelian randomization have important utility for understanding the causal relationships between major lipids and CHD risk, and are likely to play an important role in judging the causal relevance of lipid subspecies and other metabolites measured on high-dimensional phenotyping platforms.

Translational Bioinformatics: Past, Present, and Future

Though a relatively young discipline, translational bioinformatics (TBI) has become a key component of biomedical research in the era of precision medicine. Development of high-throughput technologies and electronic health records has caused a paradigm shift in both healthcare and biomedical research. Novel tools and methods are required to convert increasingly voluminous datasets into information and actionable knowledge. This review provides a definition and contextualization of the term TBI, describes the discipline’s brief history and past accomplishments, as well as current foci, and concludes with predictions of future directions in the field.