The HuMet Repository: Watching human metabolism at work.

The human metabolism constantly responds to stimuli such as food intake, fasting, exercise, and stress, triggering adaptive biochemical processes across multiple metabolic pathways. To understand the role of these processes and disruptions thereof in health and disease, detailed documentation of healthy metabolic responses is needed but still scarce on a time-resolved metabolome-wide level. Here, we present the HuMet Repository, a web-based resource for exploring dynamic metabolic responses to six physiological challenges (exercise, 36 h fasting, oral glucose and lipid loads, mixed meal, cold stress) in healthy subjects. For building this resource, we integrated existing and newly derived metabolomics data measured in blood, urine, and breath samples of 15 young healthy men at up to 56 time points during the six highly standardized challenge tests conducted over four days. The data comprise 1.1 million data points acquired on multiple platforms with temporal profiles of 2,656 metabolites from a broad range of biochemical pathways. By embedding the dataset into an interactive web application, we enable users to easily access, search, filter, analyze, and visualize the time-resolved metabolomic readouts and derived results. Users can put metabolites into their larger context by identifying metabolites with similar trajectories or by visualizing metabolites within holistic metabolic networks to pinpoint pathways of interest. In three showcases, we outline the value of the repository for gaining biological insights and generating hypotheses by analyzing the wash-out of dietary markers, the complementarity of metabolomics platforms in dynamic versus cross-sectional data, and similarities and differences in systemic metabolic responses across challenges. With its comprehensive collection of time-resolved metabolomics data, the HuMet Repository, freely accessible at https://humet.org/, is a reference for normal, healthy responses to metabolic challenges in young males. It will enable researchers with and without computational expertise, to flexibly query the data for their own research into the dynamics of human metabolism.

Artificial intelligence and digital biomarker in precision pathology guiding immune therapy selection and precision oncology.

BACKGROUND: The currently available immunotherapies already changed the strategy how many cancers are treated from first to last line. Understanding even the most complex heterogeneity in tumor tissue and mapping the spatial cartography of the tumor immunity allows the best and optimized selection of immune modulating agents to (re-)activate the patient's immune system and direct it against the individual cancer in the most effective way. RECENT FINDINGS: Primary cancer and metastases maintain a high degree of plasticity to escape any immune surveillance and continue to evolve depending on many intrinsic and extrinsic factors In the field of immune-oncology (IO) immune modulating agents are recognized as practice changing therapeutic modalities. Recent studies have shown that an optimal and lasting efficacy of IO therapeutics depends on the understanding of the spatial communication network and functional context of immune and cancer cells within the tumor microenvironment. Artificial intelligence (AI) provides an insight into the immune-cancer-network through the visualization of very complex tumor and immune interactions in cancer tissue specimens and allows the computer-assisted development and clinical validation of such digital biomarker. CONCLUSIONS: The successful implementation of AI-supported digital biomarker solutions guides the clinical selection of effective immune therapeutics based on the retrieval and visualization of spatial and contextual information from cancer tissue images and standardized data. As such, computational pathology (CP) turns into "precision pathology" delivering individual therapy response prediction. Precision Pathology does not only include digital and computational solutions but also high levels of standardized processes in the routine histopathology workflow and the use of mathematical tools to support clinical and diagnostic decisions as the basic principle of a "precision oncology".

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.

Dynamic patterns of postprandial metabolic responses to three dietary challenges.

Food intake triggers extensive changes in the blood metabolome. The kinetics of these changes depend on meal composition and on intrinsic, health-related characteristics of each individual, making the assessment of changes in the postprandial metabolome an opportunity to assess someone's metabolic status. To enable the usage of dietary challenges as diagnostic tools, profound knowledge about changes that occur in the postprandial period in healthy individuals is needed. In this study, we characterize the time-resolved changes in plasma levels of 634 metabolites in response to an oral glucose tolerance test (OGTT), an oral lipid tolerance test (OLTT), and a mixed meal (SLD) in healthy young males (n = 15). Metabolite levels for samples taken at different time points (20 per individual) during the challenges were available from targeted (132 metabolites) and non-targeted (502 metabolites) metabolomics. Almost half of the profiled metabolites (n = 308) showed a significant change in at least one challenge, thereof 111 metabolites responded exclusively to one particular challenge. Examples include azelate, which is linked to ω-oxidation and increased only in OLTT, and a fibrinogen cleavage peptide that has been linked to a higher risk of cardiovascular events in diabetes patients and increased only in OGTT, making its postprandial dynamics a potential target for risk management. A pool of 89 metabolites changed their plasma levels during all three challenges and represents the core postprandial response to food intake regardless of macronutrient composition. We used fuzzy c-means clustering to group these metabolites into eight clusters based on commonalities of their dynamic response patterns, with each cluster following one of four primary response patterns: (i) "decrease-increase" (valley-like) with fatty acids and acylcarnitines indicating the suppression of lipolysis, (ii) "increase-decrease" (mountain-like) including a cluster of conjugated bile acids and the glucose/insulin cluster, (iii) "steady decrease" with metabolites reflecting a carryover from meals prior to the study, and (iv) "mixed" decreasing after the glucose challenge and increasing otherwise. Despite the small number of subjects, the diversity of the challenges and the wealth of metabolomic data make this study an important step toward the characterization of postprandial responses and the identification of markers of metabolic processes regulated by food intake.

A Predictive Model for Progression of CKD to Kidney Failure Based on Routine Laboratory Tests.

RATIONALE & OBJECTIVE: Stratification of chronic kidney disease (CKD) patients at risk for progressing to kidney failure requiring kidney replacement therapy (KFRT) is important for clinical decision-making and trial enrollment. STUDY DESIGN: Four independent prospective observational cohort studies. SETTING & PARTICIPANTS: The development cohort comprised 4,915 CKD patients, and 3 independent validation cohorts comprised a total of 3,063. Patients were observed for approximately 5 years. EXPOSURE: 22 demographic, anthropometric, and laboratory variables commonly assessed in CKD patients. OUTCOME: Progression to KFRT. ANALYTICAL APPROACH: A least absolute shrinkage and selection operator (LASSO) Cox proportional hazards model was fit to select laboratory variables that best identified patients at high risk for KFRT. Model discrimination and calibration were assessed and compared against the 4-variable Tangri (T4) risk equation both in a resampling approach within the development cohort and in the validation cohorts using cause-specific concordance (C) statistics, net reclassification improvement, and calibration graphs. RESULTS: The newly derived 6-variable risk score (Z6) included serum creatinine, albumin, cystatin C, and urea, as well as hemoglobin and the urinary albumin-creatinine ratio. In the the resampling approach, Z6 achieved a median C statistic of 0.909 (95% CI, 0.868-0.937) at 2 years after the baseline visit, whereas the T4 achieved a median C statistic of 0.855 (95% CI, 0.799-0.915). In the 3 independent validation cohorts, the Z6C statistics were 0.894, 0.921, and 0.891, whereas the T4C statistics were 0.882, 0.913, and 0.862. LIMITATIONS: The Z6 was both derived and tested only in White European cohorts. CONCLUSIONS: A new risk equation based on 6 routinely available laboratory tests facilitates identification of patients with CKD who are at high risk of progressing to KFRT.

Mapping the proteo-genomic convergence of human diseases.

Characterization of the genetic regulation of proteins is essential for understanding disease etiology and developing therapies. We identified 10,674 genetic associations for 3892 plasma proteins to create a cis-anchored gene-protein-disease map of 1859 connections that highlights strong cross-disease biological convergence. This proteo-genomic map provides a framework to connect etiologically related diseases, to provide biological context for new or emerging disorders, and to integrate different biological domains to establish mechanisms for known gene-disease links. Our results identify proteo-genomic connections within and between diseases and establish the value of cis-protein variants for annotation of likely causal disease genes at loci identified in genome-wide association studies, thereby addressing a major barrier to experimental validation and clinical translation of genetic discoveries.

Plasma metabolites to profile pathways in noncommunicable disease multimorbidity.

Multimorbidity, the simultaneous presence of multiple chronic conditions, is an increasing global health problem and research into its determinants is of high priority. We used baseline untargeted plasma metabolomics profiling covering >1,000 metabolites as a comprehensive readout of human physiology to characterize pathways associated with and across 27 incident noncommunicable diseases (NCDs) assessed using electronic health record hospitalization and cancer registry data from over 11,000 participants (219,415 person years). We identified 420 metabolites shared between at least 2 NCDs, representing 65.5% of all 640 significant metabolite-disease associations. We integrated baseline data on over 50 diverse clinical risk factors and characteristics to identify actionable shared pathways represented by those metabolites. Our study highlights liver and kidney function, lipid and glucose metabolism, low-grade inflammation, surrogates of gut microbial diversity and specific health-related behaviors as antecedents of common NCD multimorbidity with potential for early prevention. We integrated results into an open-access webserver ( https://omicscience.org/apps/mwasdisease/ ) to facilitate future research and meta-analyses.

A cross-platform approach identifies genetic regulators of human metabolism and health.

In cross-platform analyses of 174 metabolites, we identify 499 associations (P < 4.9 × 10-10) characterized by pleiotropy, allelic heterogeneity, large and nonlinear effects and enrichment for nonsynonymous variation. We identify a signal at GLP2R (p.Asp470Asn) shared among higher citrulline levels, body mass index, fasting glucose-dependent insulinotropic peptide and type 2 diabetes, with ß-arrestin signaling as the underlying mechanism. Genetically higher serine levels are shown to reduce the likelihood (by 95%) and predict development of macular telangiectasia type 2, a rare degenerative retinal disease. Integration of genomic and small molecule data across platforms enables the discovery of regulators of human metabolism and translation into clinical insights.

Genetic architecture of host proteins involved in SARS-CoV-2 infection.

Understanding the genetic architecture of host proteins interacting with SARS-CoV-2 or mediating the maladaptive host response to COVID-19 can help to identify new or repurpose existing drugs targeting those proteins. We present a genetic discovery study of 179 such host proteins among 10,708 individuals using an aptamer-based technique. We identify 220 host DNA sequence variants acting in cis (MAF 0.01-49.9%) and explaining 0.3-70.9% of the variance of 97 of these proteins, including 45 with no previously known protein quantitative trait loci (pQTL) and 38 encoding current drug targets. Systematic characterization of pQTLs across the phenome identified protein-drug-disease links and evidence that putative viral interaction partners such as MARK3 affect immune response. Our results accelerate the evaluation and prioritization of new drug development programmes and repurposing of trials to prevent, treat or reduce adverse outcomes. Rapid sharing and detailed interrogation of results is facilitated through an interactive webserver ( https://omicscience.org/apps/covidpgwas/ ).

Genetic architecture of host proteins interacting with SARS-CoV-2.

Strategies to develop therapeutics for SARS-CoV-2 infection may be informed by experimental identification of viral-host protein interactions in cellular assays and measurement of host response proteins in COVID-19 patients. Identification of genetic variants that influence the level or activity of these proteins in the host could enable rapid 'in silico' assessment in human genetic studies of their causal relevance as molecular targets for new or repurposed drugs to treat COVID-19. We integrated large-scale genomic and aptamer-based plasma proteomic data from 10,708 individuals to characterize the genetic architecture of 179 host proteins reported to interact with SARS-CoV-2 proteins or to participate in the host response to COVID-19. We identified 220 host DNA sequence variants acting in cis (MAF 0.01-49.9%) and explaining 0.3-70.9% of the variance of 97 of these proteins, including 45 with no previously known protein quantitative trait loci (pQTL) and 38 encoding current drug targets. Systematic characterization of pQTLs across the phenome identified protein-drug-disease links, evidence that putative viral interaction partners such as MARK3 affect immune response, and establish the first link between a recently reported variant for respiratory failure of COVID-19 patients at the ABO locus and hypercoagulation, i.e. maladaptive host response. Our results accelerate the evaluation and prioritization of new drug development programmes and repurposing of trials to prevent, treat or reduce adverse outcomes. Rapid sharing and dynamic and detailed interrogation of results is facilitated through an interactive webserver ( https://omicscience.org/apps/covidpgwas/ ).

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.

Genome-wide scan identifies novel genetic loci regulating salivary metabolite levels.

Saliva, as a biofluid, is inexpensive and non-invasive to obtain, and provides a vital tool to investigate oral health and its interaction with systemic health conditions. There is growing interest in salivary biomarkers for systemic diseases, notably cardiovascular disease. Whereas hundreds of genetic loci have been shown to be involved in the regulation of blood metabolites, leading to significant insights into the pathogenesis of complex human diseases, little is known about the impact of host genetics on salivary metabolites. Here we report the first genome-wide association study exploring 476 salivary metabolites in 1419 subjects from the TwinsUK cohort (discovery phase), followed by replication in the Study of Health in Pomerania (SHIP-2) cohort. A total of 14 distinct locus-metabolite associations were identified in the discovery phase, most of which were replicated in SHIP-2. While only a limited number of the loci that are known to regulate blood metabolites were also associated with salivary metabolites in our study, we identified several novel saliva-specific locus-metabolite associations, including associations for the AGMAT (with the metabolites 4-guanidinobutanoate and beta-guanidinopropanoate), ATP13A5 (with the metabolite creatinine) and DPYS (with the metabolites 3-ureidopropionate and 3-ureidoisobutyrate) loci. Our study suggests that there may be regulatory pathways of particular relevance to the salivary metabolome. In addition, some of our findings may have clinical significance, such as the utility of the pyrimidine (uracil) degradation metabolites in predicting 5-fluorouracil toxicity and the role of the agmatine pathway metabolites as biomarkers of oral health.

HitPickV2: a web server to predict targets of chemical compounds.

MOTIVATION: The identification of protein targets of novel compounds is essential to understand compounds' mechanisms of action leading to biological effects. Experimental methods to determine these protein targets are usually slow, costly and time consuming. Computational tools have recently emerged as cheaper and faster alternatives that allow the prediction of targets for a large number of compounds. RESULTS: Here, we present HitPickV2, a novel ligand-based approach for the prediction of human druggable protein targets of multiple compounds. For each query compound, HitPickV2 predicts up to 10 targets out of 2739 human druggable proteins. To that aim, HitPickV2 identifies the closest, structurally similar compounds in a restricted space within a vast chemical-protein interaction area, until 10 distinct protein targets are found. Then, HitPickV2 scores these 10 targets based on three parameters of the targets in such space: the Tanimoto coefficient (Tc) between the query and the most similar compound interacting with the target, a target rank that considers Tc and Laplacian-modified naïve Bayesian target models scores and a novel parameter introduced in HitPickV2, the number of compounds interacting with each target (occur). We present the performance results of HitPickV2 in cross-validation as well as in an external dataset. AVAILABILITY AND IMPLEMENTATION: HitPickV2 is available in www.hitpickv2.com. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

Characterization of missing values in untargeted MS-based metabolomics data and evaluation of missing data handling strategies.

BACKGROUND: Untargeted mass spectrometry (MS)-based metabolomics data often contain missing values that reduce statistical power and can introduce bias in biomedical studies. However, a systematic assessment of the various sources of missing values and strategies to handle these data has received little attention. Missing data can occur systematically, e.g. from run day-dependent effects due to limits of detection (LOD); or it can be random as, for instance, a consequence of sample preparation. METHODS: We investigated patterns of missing data in an MS-based metabolomics experiment of serum samples from the German KORA F4 cohort (n = 1750). We then evaluated 31 imputation methods in a simulation framework and biologically validated the results by applying all imputation approaches to real metabolomics data. We examined the ability of each method to reconstruct biochemical pathways from data-driven correlation networks, and the ability of the method to increase statistical power while preserving the strength of established metabolic quantitative trait loci. RESULTS: Run day-dependent LOD-based missing data accounts for most missing values in the metabolomics dataset. Although multiple imputation by chained equations performed well in many scenarios, it is computationally and statistically challenging. K-nearest neighbors (KNN) imputation on observations with variable pre-selection showed robust performance across all evaluation schemes and is computationally more tractable. CONCLUSION: Missing data in untargeted MS-based metabolomics data occur for various reasons. Based on our results, we recommend that KNN-based imputation is performed on observations with variable pre-selection since it showed robust results in all evaluation schemes.

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.

Fine Mapping and Functional Analysis Reveal a Role of SLC22A1 in Acylcarnitine Transport.

Genome-wide association studies have identified a signal at the SLC22A1 locus for serum acylcarnitines, intermediate metabolites of mitochondrial oxidation whose plasma levels associate with metabolic diseases. Here, we refined the association signal, performed conditional analyses, and examined the linkage structure to find coding variants of SLC22A1 that mediate independent association signals at the locus. We also employed allele-specific expression analysis to find potential regulatory variants of SLC22A1 and demonstrated the effect of one variant on the splicing of SLC22A1. SLC22A1 encodes a hepatic plasma membrane transporter whose role in acylcarnitine physiology has not been described. By targeted metabolomics and isotope tracing experiments in loss- and gain-of-function cell and mouse models of Slc22a1, we uncovered a role of SLC22A1 in the efflux of acylcarnitines from the liver to the circulation. We further validated the impacts of human variants on SLC22A1-mediated acylcarnitine efflux in vitro, explaining their association with serum acylcarnitine levels. Our findings provide the detailed molecular mechanisms of the GWAS association for serum acylcarnitines at the SLC22A1 locus by functionally validating the impact of SLC22A1 and its variants on acylcarnitine transport.

Connecting genetic risk to disease end points through the human blood plasma proteome.

Genome-wide association studies (GWAS) with intermediate phenotypes, like changes in metabolite and protein levels, provide functional evidence to map disease associations and translate them into clinical applications. However, although hundreds of genetic variants have been associated with complex disorders, the underlying molecular pathways often remain elusive. Associations with intermediate traits are key in establishing functional links between GWAS-identified risk-variants and disease end points. Here we describe a GWAS using a highly multiplexed aptamer-based affinity proteomics platform. We quantify 539 associations between protein levels and gene variants (pQTLs) in a German cohort and replicate over half of them in an Arab and Asian cohort. Fifty-five of the replicated pQTLs are located in trans. Our associations overlap with 57 genetic risk loci for 42 unique disease end points. We integrate this information into a genome-proteome network and provide an interactive web-tool for interrogations. Our results provide a basis for novel approaches to pharmaceutical and diagnostic applications.