Genome-Wide Association Study with Targeted and Non-targeted NMR Metabolomics Identifies 15 Novel Loci of Urinary Human Metabolic Individuality

Genetic Architecture and Analysis Practices of Circulating Metabolites in the NHLBI Trans-Omics for Precision Medicine (TOPMed) Program

Circulating metabolite levels partly reflect the state of human health and diseases and can be impacted by genetic determinants. Hundreds of loci associated with circulating metabolites have been identified; however, most findings focus on predominantly European ancestry or single-study analyses. Leveraging the rich metabolomics resources generated by the NHLBI Trans-Omics for Precision Medicine (TOPMed) Program, we harmonized and accessibly cataloged 1,729 circulating metabolites among 25,058 ancestrally diverse samples. We provided a set of reasonable strategies for outlier and imputation handling to process metabolite data. Following the practical analysis framework, we further performed a genome-wide association analysis on 1,135 selected metabolites using whole genome sequencing data from 16,359 individuals passing the quality control filters, and discovered 1,778 independent loci associated with 667 metabolites. Among 108 novel locus-metabolite pairs, we detected not only novel loci within previously implicated metabolite associated genes but also novel genes (such as GAB3 and VSIG4 located in the X chromosome) that have putative roles in metabolic regulation. In the sex-stratified analysis, we revealed 85 independent locus-metabolite pairs with evidence of sexual dimorphism, including well-known metabolic genes such as FADS2 , D2HGDH , SUGP1 , UTG2B17 , strongly supporting the importance of exploring sex difference in the human metabolome. Taken together, our study depicted the genetic contribution to circulating metabolite levels, providing additional insight into the understanding of human health.

Genome-wide characterization of 54 urinary metabolites reveals molecular impact of kidney function

Dissecting the genetic mechanisms underlying urinary metabolite concentrations can provide molecular insights into kidney function and open possibilities for causal assessment of urinary metabolites with risk factors and disease outcomes. Proton nuclear magnetic resonance metabolomics provides a high-throughput means for urinary metabolite profiling, as widely applied for blood biomarker studies. Here we report a genome-wide association study meta-analysed for 3 European cohorts comprising 8,011 individuals, covering both people with type 1 diabetes and general population settings. We identify 54 associations (p < 9.3 × 10-10) for 19 of 54 studied metabolite concentrations. Out of these, 33 were not reported previously for relevant urinary or blood metabolite traits. Subsequent two-sample Mendelian randomization analysis suggests that estimated glomerular filtration rate causally affects 13 urinary metabolite concentrations whereas urinary ethanolamine, an initial precursor for phosphatidylcholine and phosphatidylethanolamine, was associated with higher eGFR lending support for a potential protective role. Our study provides a catalogue of genetic associations for 53 metabolites, enabling further investigation on how urinary metabolites are linked to human health.

Restricting lysine normalizes toxic catabolites associated with ALDH7A1 deficiency in cells and mice

Lysine metabolism converges at α-aminoadipic semialdehyde dehydrogenase (ALDH7A1). Rare loss-of-function mutations in ALDH7A1 cause a toxic accumulation of lysine catabolites, including piperideine-6-carboxylate (P6C), that are thought to cause fatal seizures in children unless strictly managed with dietary lysine reduction. In this study, we perform metabolomics and expression analysis of tissues from Aldh7a1-deficient mice, which reveal tissue-specific differences in lysine metabolism and other metabolic pathways. We also develop a fluorescent biosensor to characterize lysine transporter activity and identify competitive substrates that reduce the accumulation of lysine catabolites in ALDH7A1-deficient HEK293 cells. Lastly, we show that intravenous administration of lysine α-oxidase from Trichoderma viride reduces lysine and P6C levels by >80% in mice. Our results improve our understanding of lysine metabolism and make inroads toward improving therapeutic strategies for lysine catabolic disorders.

Biological sex affects functional variation across the human genome

Humans display sexual dimorphism across many traits, but little is known about underlying genetic mechanisms and impacts on disease. We utilized single-cell RNA-seq of 480 lymphoblastoid cell lines to reveal that the vast majority (79%) of sex-biased genes are targets of transcription factors that display sex-biased expression. Further, we developed a two-step regression method that identified sex-biased expression quantitative trait loci (sb-eQTL) across the genome. In contrast to previous work, these sb-eQTL are abundant (n=10,754; FDR 5%) and reproducible (replication up to π1=0.56). These sb-eQTL are enriched in over 600 GWAS phenotypes, including 120 sb-eQTL associated with the female-biased autoimmune disease multiple sclerosis. Our results demonstrate widespread genetic impacts on sexual dimorphism and identify possible mechanisms and clinical targets for sex differences in diverse diseases.

Genetic variants affecting mitochondrial function provide further insights for kidney disease

BACKGROUND

Chronic kidney disease (CKD) is a complex disorder that has become a high prevalence global health problem, with diabetes being its predominant pathophysiologic driver. Autosomal genetic variation only explains some of the predisposition to kidney disease. Variations in the mitochondrial genome (mtDNA) and nuclear-encoded mitochondrial genes (NEMG) are implicated in susceptibility to kidney disease and CKD progression, but they have not been thoroughly explored. Our aim was to investigate the association of variation in both mtDNA and NEMG with CKD (and related traits), with a particular focus on diabetes.

METHODS

We used the UK Biobank (UKB) and UK-ROI, an independent collection of individuals with type 1 diabetes mellitus (T1DM) patients.

RESULTS

Fourteen mitochondrial variants were associated with estimated glomerular filtration rate (eGFR) in UKB. Mitochondrial variants and haplogroups U, H and J were associated with eGFR and serum variables. Mitochondrial haplogroup H was associated with all the serum variables regardless of the presence of diabetes. Mitochondrial haplogroup X was associated with end-stage kidney disease (ESKD) in UKB. We confirmed the influence of several known NEMG on kidney disease and function and found novel associations for SLC39A13, CFL1, ACP2 or ATP5G1 with serum variables and kidney damage, and for SLC4A1, NUP210 and MYH14 with ESKD. The G allele of TBC1D32-rs113987180 was associated with higher risk of ESKD in patients with diabetes (OR:9.879; CI95%:4.440-21.980; P = 2.0E-08). In UK-ROI, AGXT2-rs71615838 and SURF1-rs183853102 were associated with diabetic nephropathies, and TFB1M-rs869120 with eGFR.

CONCLUSIONS

We identified novel variants both in mtDNA and NEMG which may explain some of the missing heritability for CKD and kidney phenotypes. We confirmed the role of MT-ND5 and mitochondrial haplogroup H on renal disease (serum variables), and identified the MT-ND5-rs41535848G variant, along with mitochondrial haplogroup X, associated with higher risk of ESKD. Despite most of the associations were independent of diabetes, we also showed potential roles for NEMG in T1DM.

Causal relationships of metabolites with allergic diseases: a trans-ethnic Mendelian randomization study

BACKGROUND

Allergic diseases exert a considerable impact on global health, thus necessitating investigations into their etiology and pathophysiology for devising effective prevention and treatment strategies. This study employs a Mendelian randomization (MR) analysis and meta-analysis to identify metabolite targets potentially associated with allergic diseases.

METHODS

A two-sample MR analysis was conducted to explore potential causal relationships between circulating and urinary metabolites and allergic diseases. Exposures were derived from a genome-wide association study (GWAS) of 486 circulating metabolites and a GWAS of 55 targeted urinary metabolites. Outcome data for allergic diseases, including atopic dermatitis (AD), allergic rhinitis (AR), and asthma, were obtained from the FinnGen biobank in Europe (cohort 1) and the Biobank Japan in Asia (cohort 2). MR results from both cohorts were combined using a meta-analysis.

RESULTS

MR analysis identified 50 circulating metabolites and 6 urinary metabolites in cohort 1 and 54 circulating metabolites and 2 urinary metabolites in cohort 2 as potentially causally related to allergic diseases. A meta-analysis of the MR results revealed stearoylcarnitine (OR 8.654; 95% CI 4.399-17.025; P = 4.06E-10) and 1-arachidonoylglycerophosphoinositol (OR 2.178; 95% CI 1.388-3.419; P = 7.15E-04) as the most reliable causal circulating metabolites for asthma and AR, respectively. Further, histidine (OR 0.734; 95% CI: 0.594-0.907; P = 0.004), tyrosine (OR 0.601; 95% CI: 0.380-0.952; P = 0.030), and alanine (OR 0.280; 95% CI: 0.125-0.628; P = 0.002) emerged as urinary metabolites with the greatest protective effects against asthma, AD, and AR, respectively.

CONCLUSIONS

Imbalances in numerous circulating and urinary metabolites may be implicated in the development and progression of allergic diseases. These findings have significant implications for the development of targeted strategies for the prevention and treatment of allergic diseases.

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

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

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Multi-omics research strategies in ischemic stroke: A multidimensional perspective

Ischemic stroke (IS) is a multifactorial and heterogeneous neurological disorder with high rate of death and long-term impairment. Despite years of studies, there are still no stroke biomarkers for clinical practice, and the molecular mechanisms of stroke remain largely unclear. The high-throughput omics approach provides new avenues for discovering biomarkers of IS and explaining its pathological mechanisms. However, single-omics approaches only provide a limited understanding of the biological pathways of diseases. The integration of multiple omics data means the simultaneous analysis of thousands of genes, RNAs, proteins and metabolites, revealing networks of interactions between multiple molecular levels. Integrated analysis of multi-omics approaches will provide helpful insights into stroke pathogenesis, therapeutic target identification and biomarker discovery. Here, we consider advances in genomics, transcriptomics, proteomics and metabolomics and outline their use in discovering the biomarkers and pathological mechanisms of IS. We then delineate strategies for achieving integration at the multi-omics level and discuss how integrative omics and systems biology can contribute to our understanding and management of IS.

GABA signaling enforces intestinal germinal center B cell differentiation

Recent compelling results indicate possible links between neurotransmitters, intestinal mucosal IgA + B cell responses, and immunoglobulin A nephropathy (IgAN) pathogenesis. Here, we demonstrated that γ-amino butyric acid (GABA) transporter-2 (GAT-2) deficiency induces intestinal germinal center (GC) B cell differentiation and worsens the symptoms of IgAN in a mouse model. Mechanistically, GAT-2 deficiency enhances GC B cell differentiation through activation of GABA–mammalian target of rapamycin complex 1 (mTORC1) signaling. In addition, IgAN patients have lower GAT-2 expression but higher activation of mTORC1 in blood B cells, and both are correlated with kidney function in IgAN patients. Collectively, this study describes GABA signaling–mediated intestinal mucosal immunity as a previously unstudied pathogenesis mechanism of IgAN and challenges the current paradigms of IgAN.

mGWAS-Explorer: Linking SNPs, Genes, Metabolites, and Diseases for Functional Insights

Tens of thousands of single-nucleotide polymorphisms (SNPs) have been identified to be significantly associated with metabolite abundance in over 65 genome-wide association studies with metabolomics (mGWAS) to date. Obtaining mechanistic or functional insights from these associations for translational applications has become a key research area in the mGWAS community. Here, we introduce mGWAS-Explorer, a user-friendly web-based platform to help connect SNPs, metabolites, genes, and their known disease associations via powerful network visual analytics. The application of the mGWAS-Explorer was demonstrated using a COVID-19 and a type 2 diabetes case studies.

Network Approaches to Integrate Analyses of Genetics and Metabolomics Data with Applications to Fetal Programming Studies

The integration of genetics and metabolomics data demands careful accounting of complex dependencies, particularly when modelling familial omics data, e.g., to study fetal programming of related maternal–offspring phenotypes. Efforts to identify genetically determined metabotypes using classic genome wide association approaches have proven useful for characterizing complex disease, but conclusions are often limited to a series of variant–metabolite associations. We adapt Bayesian network models to integrate metabotypes with maternal–offspring genetic dependencies and metabolic profile correlations in order to investigate mechanisms underlying maternal–offspring phenotypic associations. Using data from the multiethnic Hyperglycemia and Adverse Pregnancy Outcome (HAPO) study, we demonstrate that the strategic specification of ordered dependencies, pre-filtering of candidate metabotypes, incorporation of metabolite dependencies, and penalized network estimation methods clarify potential mechanisms for fetal programming of newborn adiposity and metabolic outcomes. The exploration of Bayesian network growth over a range of penalty parameters, coupled with interactive plotting, facilitate the interpretation of network edges. These methods are broadly applicable to integration of diverse omics data for related individuals.

Metabolome Genome-Wide Association Study Identifies 74 Novel Genomic Regions Influencing Plasma Metabolites Levels

Metabolites are small products of metabolism that provide a snapshot of the wellbeing of an organism and the mechanisms that control key physiological processes involved in health and disease. Here we report the results of a genome-wide association study of 722 circulating metabolite levels in 8809 subjects of European origin, providing both breadth and depth. These analyses identified 202 unique genomic regions whose variations are associated with the circulating levels of 478 different metabolites. Replication with a subset of 208 metabolites that were available in an independent dataset for a cohort of 1768 European subjects confirmed the robust associations, including 74 novel genomic regions not associated with any metabolites in previous works. This study enhances our knowledge of genetic mechanisms controlling human metabolism. Our findings have major potential for identifying novel targets and developing new therapeutic strategies.

Serum Glycine Levels Are Associated With Cortical Bone Properties and Fracture Risk in Men

CONTEXT

In a recent study a pattern of 27 metabolites, including serum glycine, associated with bone mineral density (BMD).

OBJECTIVE

To investigate associations for serum and urinary glycine levels with BMD, bone microstructure, and fracture risk in men.

METHODS

In the population-based Osteoporotic Fractures in Men (MrOS) Sweden study (men, 69-81 years) serum glycine and BMD were measured at baseline (n = 965) and 5-year follow-up (n = 546). Cortical and trabecular bone parameters of the distal tibia were measured at follow-up using high-resolution peripheral quantitative computed tomography. Urinary (n = 2682) glycine was analyzed at baseline. X-ray-validated fractures (n = 594) were ascertained during a median follow-up of 9.6 years. Associations were evaluated using linear regression (bone parameters) or Cox regression (fractures).

RESULTS

Circulating glycine levels were inversely associated with femoral neck (FN)-BMD. A meta-analysis (n = 7543) combining MrOS Sweden data with data from 3 other cohorts confirmed a robust inverse association between serum glycine levels and FN-BMD (P = 7.7 × 10-9). Serum glycine was inversely associated with the bone strength parameter failure load in the distal tibia (P = 0.002), mainly as a consequence of an inverse association with cortical cross-sectional area and a direct association with cortical porosity. Both serum and urinary glycine levels predicted major osteoporotic fractures (serum: hazard ratio [HR] per SD increase = 1.22, 95% CI, 1.05-1.43; urine: HR = 1.13, 95% CI, 1.02-1.24). These fracture associations were only marginally reduced in models adjusted by FRAX with BMD.

CONCLUSIONS

Serum and urinary glycine are indirectly associated with FN-BMD and cortical bone strength, and directly associated with fracture risk in men.

Untargeted Metabolome- and Transcriptome-Wide Association Study Suggests Causal Genes Modulating Metabolite Concentrations in Urine

Gene products can affect the concentrations of small molecules (aka “metabolites”), and conversely, some metabolites can modulate the concentrations of gene transcripts. While many specific instances of this interplay have been revealed, a global approach to systematically uncover human gene-metabolite interactions is still lacking. We performed a metabolome- and transcriptome-wide association study to identify genes influencing the human metabolome using untargeted metabolome features, extracted from ¹H nuclear magnetic resonance spectroscopy (NMR) of urine samples, and gene expression levels, quantified from RNA-Seq of lymphoblastoid cell lines (LCL) from 555 healthy individuals. We identified 20 study-wide significant associations corresponding to 15 genes, of which 5 associations (with 2 genes) were confirmed with follow-up NMR data. Using metabomatching, we identified the metabolites corresponding to metabolome features associated with the genes, namely, N-acetylated compounds with ALMS1 and trimethylamine (TMA) with HPS1. Finally, Mendelian randomization analysis supported a potential causal link between the expression of genes in both the ALMS1- and HPS1-loci and their associated metabolite concentrations. In the case of HPS1, we additionally observed that TMA concentration likely exhibits a reverse causal effect on HPS1 expression levels, indicating a negative feedback loop. Our study highlights how the integration of metabolomics, gene expression, and genetic data can pinpoint causal genes modulating metabolite concentrations.

Metabolomics in Parkinson's disease

Parkinson's disease (PD) is a multifactorial neurodegenerative disorder in which environmental (lifestyle, dietary, infectious disease) factors as well as genetic make-up play a role. Metabolomics, an evolving research field combining biomarker discovery and pathogenetics, is particularly useful in studying complex pathophysiology in general and Parkinson's disease (PD) specifically. PD, the second most frequent neurodegenerative disorder, is characterized by the loss of dopaminergic neurons in the substantia nigra and the presence of intraneural inclusions of α-synuclein aggregates. Although considered a predominantly movement disorder, PD is also associated with number of non-motor features. Metabolomics has provided useful information regarding this neurodegenerative process with the aim of identifying a disease-specific fingerprint. Unfortunately, many disease variables such as clinical presentation, motor system involvement, disease stage and duration substantially affect biomarker relevance. As such, metabolomics provides a unique approach to studying this multifactorial neurodegenerative disorder.

An atlas connecting shared genetic architecture of human diseases and molecular phenotypes provides insight into COVID-19 susceptibility

BACKGROUND

While genome-wide associations studies (GWAS) have successfully elucidated the genetic architecture of complex human traits and diseases, understanding mechanisms that lead from genetic variation to pathophysiology remains an important challenge. Methods are needed to systematically bridge this crucial gap to facilitate experimental testing of hypotheses and translation to clinical utility.

RESULTS

Here, we leveraged cross-phenotype associations to identify traits with shared genetic architecture, using linkage disequilibrium (LD) information to accurately capture shared SNPs by proxy, and calculate significance of enrichment. This shared genetic architecture was examined across differing biological scales through incorporating data from catalogs of clinical, cellular, and molecular GWAS. We have created an interactive web database (interactive Cross-Phenotype Analysis of GWAS database (iCPAGdb)) to facilitate exploration and allow rapid analysis of user-uploaded GWAS summary statistics. This database revealed well-known relationships among phenotypes, as well as the generation of novel hypotheses to explain the pathophysiology of common diseases. Application of iCPAGdb to a recent GWAS of severe COVID-19 demonstrated unexpected overlap of GWAS signals between COVID-19 and human diseases, including with idiopathic pulmonary fibrosis driven by the DPP9 locus. Transcriptomics from peripheral blood of COVID-19 patients demonstrated that DPP9 was induced in SARS-CoV-2 compared to healthy controls or those with bacterial infection. Further investigation of cross-phenotype SNPs associated with both severe COVID-19 and other human traits demonstrated colocalization of the GWAS signal at the ABO locus with plasma protein levels of a reported receptor of SARS-CoV-2, CD209 (DC-SIGN). This finding points to a possible mechanism whereby glycosylation of CD209 by ABO may regulate COVID-19 disease severity.

CONCLUSIONS

Thus, connecting genetically related traits across phenotypic scales links human diseases to molecular and cellular measurements that can reveal mechanisms and lead to novel biomarkers and therapeutic approaches. The iCPAGdb web portal is accessible at http://cpag.oit.duke.edu and the software code at https://github.com/tbalmat/iCPAGdb .

Causal Inference Methods to Integrate Omics and Complex Traits

Major biotechnological advances have facilitated a tremendous boost to the collection of (gen-/transcript-/prote-/methyl-/metabol-)omics data in very large sample sizes worldwide. Coordinated efforts have yielded a deluge of studies associating diseases with genetic markers (genome-wide association studies) or with molecular phenotypes. Whereas omics-disease associations have led to biologically meaningful and coherent mechanisms, the identified (non-germline) disease biomarkers may simply be correlates or consequences of the explored diseases. To move beyond this realm, Mendelian randomization provides a principled framework to integrate information on omics- and disease-associated genetic variants to pinpoint molecular traits causally driving disease development. In this review, we show the latest advances in this field, flag up key challenges for the future, and propose potential solutions.

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.

Cerebrospinal fluid metabolomics identifies 19 brain-related phenotype associations

The study of metabolomics and disease has enabled the discovery of new risk factors, diagnostic markers, and drug targets. For neurological and psychiatric phenotypes, the cerebrospinal fluid (CSF) is of particular importance. However, the CSF metabolome is difficult to study on a large scale due to the relative complexity of the procedure needed to collect the fluid. Here, we present a metabolome-wide association study (MWAS), which uses genetic and metabolomic data to impute metabolites into large samples with genome-wide association summary statistics. We conduct a metabolome-wide, genome-wide association analysis with 338 CSF metabolites, identifying 16 genotype-metabolite associations (metabolite quantitative trait loci, or mQTLs). We then build prediction models for all available CSF metabolites and test for associations with 27 neurological and psychiatric phenotypes, identifying 19 significant CSF metabolite-phenotype associations. Our results demonstrate the feasibility of MWAS to study omic data in scarce sample types.

Multi-omics integration in biomedical research - A metabolomics-centric review

Recent advances in high-throughput technologies have enabled the profiling of multiple layers of a biological system, including DNA sequence data (genomics), RNA expression levels (transcriptomics), and metabolite levels (metabolomics). This has led to the generation of vast amounts of biological data that can be integrated in so-called multi-omics studies to examine the complex molecular underpinnings of health and disease. Integrative analysis of such datasets is not straightforward and is particularly complicated by the high dimensionality and heterogeneity of the data and by the lack of universal analysis protocols. Previous reviews have discussed various strategies to address the challenges of data integration, elaborating on specific aspects, such as network inference or feature selection techniques. Thereby, the main focus has been on the integration of two omics layers in their relation to a phenotype of interest. In this review we provide an overview over a typical multi-omics workflow, focusing on integration methods that have the potential to combine metabolomics data with two or more omics. We discuss multiple integration concepts including data-driven, knowledge-based, simultaneous and step-wise approaches. We highlight the application of these methods in recent multi-omics studies, including large-scale integration efforts aiming at a global depiction of the complex relationships within and between different biological layers without focusing on a particular phenotype.

Genome-Wide Association Study of Response to Selenium Supplementation and Circulating Selenium Concentrations in Adults of European Descent

BACKGROUND

Selenium (Se) is a trace element that has been linked to many health conditions. Genome-wide association studies (GWAS) have identified variants for blood and toenail Se levels, but no GWAS has been conducted to date on responses to Se supplementation.

OBJECTIVES

A GWAS was performed to identify the single nucleotide polymorphisms (SNPs) associated with changes in Se concentrations after 1 year of supplementation. A GWAS of basal plasma Se concentrations at study entry was conducted to evaluate whether SNPs for Se responses overlap with SNPs for basal Se levels.

METHODS

A total of 428 participants aged 40-80 years of European descent from the Selenium and Celecoxib Trial (Sel/Cel Trial) who received daily supplementation with 200 µg of selenized yeast were included for the GWAS of responses to supplementation. Plasma Se concentrations were measured from blood samples collected at the time of recruitment and after 1 year of supplementation. Linear regression analyses were performed to assess the relationship between each SNP and changes in Se concentrations. We further examined whether the identified SNPs overlapped with those related to basal Se concentrations.

RESULTS

No SNP was significantly associated with changes in Se concentration at a genome-wide significance level. However, rs56856693, located upstream of the NEK6, was nominally associated with changes in Se concentrations after supplementation (P = 4.41 × 10-7), as were 2 additional SNPs, rs11960388 and rs6887869, located in the dimethylglycine dehydrogenase (DMGDH)/betaine-homocysteine S-methyltransferase (BHMT) region (P = 0.01). Alleles of 2 SNPs in the DMGDH/BHMT region associated with greater increases in Se concentrations after supplementation were also strongly associated with higher basal Se concentrations (P = 8.67 × 10-8).

CONCLUSIONS

This first GWAS of responses to Se supplementation in participants of European descent from the Sel/Cel Trial suggests that SNPs in the NEK6 and DMGDH/BHMT regions influence responses to supplementation.

An atlas connecting shared genetic architecture of human diseases and molecular phenotypes provides insight into COVID-19 susceptibility

While genome-wide associations studies (GWAS) have successfully elucidated the genetic architecture of complex human traits and diseases, understanding mechanisms that lead from genetic variation to pathophysiology remains an important challenge. Methods are needed to systematically bridge this crucial gap to facilitate experimental testing of hypotheses and translation to clinical utility. Here, we leveraged cross-phenotype associations to identify traits with shared genetic architecture, using linkage disequilibrium (LD) information to accurately capture shared SNPs by proxy, and calculate significance of enrichment. This shared genetic architecture was examined across differing biological scales through incorporating data from catalogs of clinical, cellular, and molecular GWAS. We have created an interactive web database (interactive Cross-Phenotype Analysis of GWAS database (iCPAGdb); http://cpag.oit.duke.edu) to facilitate exploration and allow rapid analysis of user-uploaded GWAS summary statistics. This database revealed well-known relationships among phenotypes, as well as the generation of novel hypotheses to explain the pathophysiology of common diseases. Application of iCPAGdb to a recent GWAS of severe COVID-19 demonstrated unexpected overlap of GWAS signals between COVID-19 and human diseases, including with idiopathic pulmonary fibrosis driven by the DPP9 locus. Transcriptomics from peripheral blood of COVID-19 patients demonstrated that DPP9 was induced in SARS-CoV-2 compared to healthy controls or those with bacterial infection. Further investigation of cross-phenotype SNPs with severe COVID-19 demonstrated colocalization of the GWAS signal of the ABO locus with plasma protein levels of a reported receptor of SARS-CoV-2, CD209 (DC-SIGN), pointing to a possible mechanism whereby glycosylation of CD209 by ABO may regulate COVID-19 disease severity. Thus, connecting genetically related traits across phenotypic scales links human diseases to molecular and cellular measurements that can reveal mechanisms and lead to novel biomarkers and therapeutic approaches.

Biochemical phenotyping of multiple myeloma patients at diagnosis reveals a disorder of mitochondrial complexes I and II and a Hartnup-like disturbance as underlying conditions, also influencing different stages of the disease

The aim of this study was to identify novel plasma metabolic signatures with possible relevance during multiple myeloma (MM) development and progression. A biochemical quantitative phenotyping platform based on targeted electrospray ionization tandem mass spectrometry technology was used to aid in the identification of any eventual perturbed biochemical pathway in peripheral blood plasma from 36 MM patients and 73 healthy controls. Our results showed that MM cases present an increase in short and medium/long-chain species of acylcarnitines resembling Multiple AcylCoA Dehydrogenase Deficiency (MADD), particularly, associated with MM advanced International Staging System (ISS). Lipids profile showed lower concentrations of phosphatidylcholine (PC), lysophosphatidylcholine (LPC) and sphingomyelins (SM) in the MM patients and its respective ISS groups. MM cases were accompanied by a drop in the concentration of essential amino acids, especially tryptophan, with a significant inverse correlation between the progressive drop in tryptophan with the elevation of β2-microglobulin, with the increase in systemic methylation levels (Symmetric Arginine Dimethylation, SDMA) and with the accumulation of esterified carnitines in relation to free carnitine (AcylC/C0). Serotonin was significantly elevated in cases of MM, without a clear association with ISS. Kynurenine/tryptophan ratio demonstrates that the activity of dioxigenases is even higher in the cases classified as ISS 3. In conclusion, our study showed that MM patients at diagnosis showed metabolic disorders resembling both mitochondrial complexes I and II and Hartnup-like disturbances as underlying conditions, also influencing different stages of the disease.

Urinary metabolic phenotyping for Alzheimer's disease

Finding early disease markers using non-invasive and widely available methods is essential to develop a successful therapy for Alzheimer’s Disease. Few studies to date have examined urine, the most readily available biofluid. Here we report the largest study to date using comprehensive metabolic phenotyping platforms (NMR spectroscopy and UHPLC-MS) to probe the urinary metabolome in-depth in people with Alzheimer’s Disease and Mild Cognitive Impairment. Feature reduction was performed using metabolomic Quantitative Trait Loci, resulting in the list of metabolites associated with the genetic variants. This approach helps accuracy in identification of disease states and provides a route to a plausible mechanistic link to pathological processes. Using these mQTLs we built a Random Forests model, which not only correctly discriminates between people with Alzheimer’s Disease and age-matched controls, but also between individuals with Mild Cognitive Impairment who were later diagnosed with Alzheimer’s Disease and those who were not. Further annotation of top-ranking metabolic features nominated by the trained model revealed the involvement of cholesterol-derived metabolites and small-molecules that were linked to Alzheimer’s pathology in previous studies.

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.

A Genome-wide Association Study Discovers 46 Loci of the Human Metabolome in the Hispanic Community Health Study/Study of Latinos

Variation in levels of the human metabolome reflect changes in homeostasis, providing a window into health and disease. The genetic impact on circulating metabolites in Hispanics, a population with high cardiometabolic disease burden, is largely unknown. We conducted genome-wide association analyses on 640 circulating metabolites in 3,926 Hispanic Community Health Study/Study of Latinos participants. The estimated heritability for 640 metabolites ranged between 0%-54% with a median at 2.5%. We discovered 46 variant-metabolite pairs (p value < 1.2 × 10-10, minor allele frequency ≥ 1%, proportion of variance explained [PEV] mean = 3.4%, PEVrange = 1%-22%) with generalized effects in two population-based studies and confirmed 301 known locus-metabolite associations. Half of the identified variants with generalized effect were located in genes, including five nonsynonymous variants. We identified co-localization with the expression quantitative trait loci at 105 discovered and 151 known loci-metabolites sets. rs5855544, upstream of SLC51A, was associated with higher levels of three steroid sulfates and co-localized with expression levels of SLC51A in several tissues. Mendelian randomization (MR) analysis identified several metabolites associated with coronary heart disease (CHD) and type 2 diabetes. For example, two variants located in or near CYP4F2 (rs2108622 and rs79400241, respectively), involved in vitamin E metabolism, were associated with the levels of octadecanedioate and vitamin E metabolites (gamma-CEHC and gamma-CEHC glucuronide); MR analysis showed that genetically high levels of these metabolites were associated with lower odds of CHD. Our findings document the genetic architecture of circulating metabolites in an underrepresented Hispanic/Latino community, shedding light on disease etiology.

Non-targeted urine metabolomics and associations with prevalent and incident type 2 diabetes

Better risk prediction and new molecular targets are key priorities in type 2 diabetes (T2D) research. Little is known about the role of the urine metabolome in predicting the risk of T2D. We aimed to use non-targeted urine metabolomics to discover biomarkers and improve risk prediction for T2D. Urine samples from two community cohorts of 1,424 adults were analyzed by ultra-performance liquid chromatography/mass spectrometry (UPLC-MS). In a discovery/replication design, three out of 62 annotated metabolites were associated with prevalent T2D, notably lower urine levels of 3-hydroxyundecanoyl-carnitine. In participants without diabetes at baseline, LASSO regression in the training set selected six metabolites that improved prediction of T2D beyond established risk factors risk over up to 12 years' follow-up in the test sample, from C-statistic 0.866 to 0.892. Our results in one of the largest non-targeted urinary metabolomics study to date demonstrate the role of the urine metabolome in identifying at-risk persons for T2D and suggest urine 3-hydroxyundecanoyl-carnitine as a biomarker candidate.

Genome-wide associations of human gut microbiome variation and implications for causal inference analyses

Recent population-based1-4 and clinical studies5 have identified a range of factors associated with human gut microbiome variation. Murine quantitative trait loci6, human twin studies7 and microbiome genome-wide association studies1,3,8-12 have provided evidence for genetic contributions to microbiome composition. Despite this, there is still poor overlap in genetic association across human studies. Using appropriate taxon-specific models, along with support from independent cohorts, we show an association between human host genotype and gut microbiome variation. We also suggest that interpretation of applied analyses using genetic associations is complicated by the probable overlap between genetic contributions and heritable components of host environment. Using faecal 16S ribosomal RNA gene sequences and host genotype data from the Flemish Gut Flora Project (n = 2,223) and two German cohorts (FoCus, n = 950; PopGen, n = 717), we identify genetic associations involving multiple microbial traits. Two of these associations achieved a study-level threshold of P = 1.57 × 10-10; an association between Ruminococcus and rs150018970 near RAPGEF1 on chromosome 9, and between Coprococcus and rs561177583 within LINC01787 on chromosome 1. Exploratory analyses were undertaken using 11 other genome-wide associations with strong evidence for association (P < 2.5 × 10-8) and a previously reported signal of association between rs4988235 (MCM6/LCT) and Bifidobacterium. Across these 14 single-nucleotide polymorphisms there was evidence of signal overlap with other genome-wide association studies, including those for age at menarche and cardiometabolic traits. Mendelian randomization analysis was able to estimate associations between microbial traits and disease (including Bifidobacterium and body composition); however, in the absence of clear microbiome-driven effects, caution is needed in interpretation. Overall, this work marks a growing catalogue of genetic associations that will provide insight into the contribution of host genotype to gut microbiome. Despite this, the uncertain origin of association signals will likely complicate future work looking to dissect function or use associations for causal inference analysis.

Urine 6-Bromotryptophan: Associations with Genetic Variants and Incident End-Stage Kidney Disease

Higher serum 6-bromotryptophan has been associated with lower risk of chronic kidney disease (CKD) progression, implicating mechanisms beyond renal clearance. We studied genetic determinants of urine 6-bromotryptophan and its association with CKD risk factors and incident end-stage kidney disease (ESKD) in 4,843 participants of the German Chronic Kidney Disease (GCKD) study. 6-bromotryptophan was measured from urine samples using mass spectrometry. Patients with higher levels of urine 6-bromotryptophan had higher baseline estimated glomerular filtration rate (eGFR, p < 0.001). A genome-wide association study of urine 6-bromotryptophan identified two significant loci possibly related to its tubular reabsorption, SLC6A19, and its production, ERO1A, which was also associated with serum 6-bromotryptophan in an independent study. The association between urine 6-bromotryptophan and time to ESKD was assessed using Cox regression. There were 216 ESKD events after four years of follow-up. Compared with patients with undetectable levels, higher 6-bromotryptophan levels were associated with lower risk of ESKD in models unadjusted and adjusted for ESKD risk factors other than eGFR ( Collapse

Key Words

• prognostic markers
• kidney diseases
• biomarkers
• epidemiology
• genetics research
• outcomes research
• risk factors

Collapse

MESH Headings

• Amino Acid Transport Systems, Neutral/genetics
• Female
• Genetic Predisposition to Disease/genetics
• Genetic Variation/genetics
• Glomerular Filtration Rate/genetics
• Humans
• Kidney Failure, Chronic/epidemiology
• Kidney Failure, Chronic/genetics
• Kidney Failure, Chronic/urine
• Male
• Membrane Glycoproteins/genetics
• Middle Aged
• Oxidoreductases/genetics
• Proportional Hazards Models
• Risk Factors
• Tryptophan/analogs & derivatives
• Tryptophan/urine

Collapse

Grants

• R21 DK112087 NIDDK NIH HHS
• UL1 RR025005 NCRR NIH HHS
• R01 HL086694 NHLBI NIH HHS
• U01 HG004402 NHGRI NIH HHS
• HHSN268201700002C NHLBI NIH HHS
• HHSN268201700001I NHLBI NIH HHS
• HHSN268201700004I NHLBI NIH HHS
• HHSN268201700004C NHLBI NIH HHS
• R01 HL087641 NHLBI NIH HHS
• HHSN268201700003I NHLBI NIH HHS
• HHSN268201700005C NHLBI NIH HHS
• HHSN268201700001C NHLBI NIH HHS
• HHSN268201700003C NHLBI NIH HHS
• HHSN268201700002I NHLBI NIH HHS
• HHSN268201700005I NHLBI NIH HHS

Collapse

Affiliation(s)

Peggy Sekula Institute of Genetic Epidemiology, Faculty of Medicine and Medical Center - University of Freiburg, Freiburg, Germany
Adrienne Tin Division of Nephrology, Department of Medicine, University of Mississippi Medical Center, Jackson, MS, USA The Memory Impairment and Neurodegenerative Dementia Center, University of Mississippi Medical Center, Jackson, MS, USA
Ulla T Schultheiss Institute of Genetic Epidemiology, Faculty of Medicine and Medical Center - University of Freiburg, Freiburg, Germany Division of Nephrology, Department of Medicine, Faculty of Medicine and Medical Center - University of Freiburg, Freiburg, Germany
Seema Baid-Agrawal Department of Nephrology and Transplant Center, Sahlgrenska University Hospital, University of Gothenburg, Gothenburg, Sweden
Robert P Mohney Metabolon Inc., Durham, NC, USA
Inga Steinbrenner Institute of Genetic Epidemiology, Faculty of Medicine and Medical Center - University of Freiburg, Freiburg, Germany
Bing Yu School of Public Health, The University of Texas Health Science Center at Houston, Houston, USA
Shengyuan Luo Department of Epidemiology, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, USA Welch Center for Prevention, Epidemiology, and Clinical Research, Johns Hopkins University, Baltimore, MD, USA
Eric Boerwinkle School of Public Health, The University of Texas Health Science Center at Houston, Houston, USA
Kai-Uwe Eckardt Department of Nephrology and Medical Intensive Care, Charité - Universitätsmedizin Berlin, Berlin, Germany Department of Nephrology and Hypertension, Friedrich-Alexander-Universität Erlangen-Nürnberg, 91054, Erlangen, Germany
Josef Coresh Department of Epidemiology, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, USA Welch Center for Prevention, Epidemiology, and Clinical Research, Johns Hopkins University, Baltimore, MD, USA
Morgan E Grams Department of Epidemiology, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, USA Welch Center for Prevention, Epidemiology, and Clinical Research, Johns Hopkins University, Baltimore, MD, USA Division of Nephrology, Johns Hopkins School of Medicine, Baltimore, MD, USA
Anna Kӧttgen Institute of Genetic Epidemiology, Faculty of Medicine and Medical Center - University of Freiburg, Freiburg, Germany. Department of Epidemiology, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, USA.

Collapse

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.

Linking glycemic dysregulation in diabetes to symptoms, comorbidities, and genetics through EHR data mining

Diabetes is a diverse and complex disease, with considerable variation in phenotypic manifestation and severity. This variation hampers the study of etiological differences and reduces the statistical power of analyses of associations to genetics, treatment outcomes, and complications. We address these issues through deep, fine-grained phenotypic stratification of a diabetes cohort. Text mining the electronic health records of 14,017 patients, we matched two controlled vocabularies (ICD-10 and a custom vocabulary developed at the clinical center Steno Diabetes Center Copenhagen) to clinical narratives spanning a 19 year period. The two matched vocabularies comprise over 20,000 medical terms describing symptoms, other diagnoses, and lifestyle factors. The cohort is genetically homogeneous (Caucasian diabetes patients from Denmark) so the resulting stratification is not driven by ethnic differences, but rather by inherently dissimilar progression patterns and lifestyle related risk factors. Using unsupervised Markov clustering, we defined 71 clusters of at least 50 individuals within the diabetes spectrum. The clusters display both distinct and shared longitudinal glycemic dysregulation patterns, temporal co-occurrences of comorbidities, and associations to single nucleotide polymorphisms in or near genes relevant for diabetes comorbidities.

Omics studies for comprehensive understanding of immunoglobulin A nephropathy: state-of-the-art and future directions

Immunoglobulin A nephropathy (IgAN) is the most common worldwide primary glomerulonephritis with a strong autoimmune component. The disease shows variability in both clinical phenotypes and endpoints and can be potentially subdivided into more homogeneous subtypes through the identification of specific molecular biomarkers. This review focuses on the role of omics in driving the identification of potential molecular subtypes of the disease through the integration of multilevel data from genomics, transcriptomics, epigenomics, proteomics and metabolomics. First, the identification of molecular biomarkers, including mapping of the full spectrum of common and rare IgAN risk alleles, could permit a more precise stratification of IgAN patients. Second, the analysis of transcriptomic patterns and their modulation by epigenetic factors like microRNAs has the potential to increase our understanding in the pathogenic mechanisms of the disease. Third, the specificity of urinary proteomic and metabolomic signatures and the understanding of their functional relevance may contribute to the development of new non-invasive biomarkers for a better molecular characterization of the renal damage and its follow-up. All these approaches can give information for targeted therapeutic decisions and will support novel clinical decision making. In conclusion, we offer a framework of omic studies and outline barriers and potential solutions that should be used for improving the diagnosis and treatment of the disease. The ongoing decade is exploiting novel high-throughput molecular technologies and computational analyses for improving the diagnosis (precision nephrology) and treatment (personalized therapy) of the IgAN subtypes.

Bivariate Genome-Wide Association Scan Identifies 6 Novel Loci Associated With Lipid Levels and Coronary Artery Disease

BACKGROUND

Plasma lipid levels are heritable and genetically associated with risk of coronary artery disease (CAD). However, genome-wide association studies (GWAS) routinely analyze these traits independent of one another. Joint GWAS for two related phenotypes can lead to a higher powered analysis to detect variants contributing to both traits.

METHODS

We performed a bivariate GWAS to discover novel loci associated with both heart disease, using a CAD meta-analysis (122 733 cases and 424 528 controls), and lipid traits, using results from the Global Lipid Genetics Consortium (188 577 subjects).

RESULTS

We identified six previously unreported loci at genome-wide significance ( P<5×10^-8), three which were associated with triglycerides and CAD, two which were associated with LDL (low-density lipoprotein) cholesterol and CAD, and one associated with total cholesterol and CAD. At several of our loci, the GWAS signals jointly localize with genetic variants associated with expression level changes for more than one neighboring genes, indicating that these loci may be affecting disease risk through regulatory activity.

CONCLUSIONS

We discovered six novel variants individually associated with both lipids and CAD.

Target genes, variants, tissues and transcriptional pathways influencing human serum urate levels

Elevated serum urate levels cause gout and correlate with cardiometabolic diseases via poorly understood mechanisms. We performed a trans-ancestry genome-wide association study of serum urate in 457,690 individuals, identifying 183 loci (147 previously unknown) that improve the prediction of gout in an independent cohort of 334,880 individuals. Serum urate showed significant genetic correlations with many cardiometabolic traits, with genetic causality analyses supporting a substantial role for pleiotropy. Enrichment analysis, fine-mapping of urate-associated loci and colocalization with gene expression in 47 tissues implicated the kidney and liver as the main target organs and prioritized potentially causal genes and variants, including the transcriptional master regulators in the liver and kidney, HNF1A and HNF4A. Experimental validation showed that HNF4A transactivated the promoter of ABCG2, encoding a major urate transporter, in kidney cells, and that HNF4A p.Thr139Ile is a functional variant. Transcriptional coregulation within and across organs may be a general mechanism underlying the observed pleiotropy between urate and cardiometabolic traits.

Automated Analysis of Large-Scale NMR Data Generates Metabolomic Signatures and Links Them to Candidate Metabolites

Identification of metabolites in large-scale ¹H NMR data from human biofluids remains challenging due to the complexity of the spectra and their sensitivity to pH and ionic concentrations. In this work, we tested the capacity of three analysis tools to extract metabolite signatures from 968 NMR profiles of human urine samples. Specifically, we studied sets of covarying features derived from principal component analysis (PCA), the iterative signature algorithm (ISA), and averaged correlation profiles (ACP), a new method we devised inspired by the STOCSY approach. We used our previously developed metabomatching method to match the sets generated by these algorithms to NMR spectra of individual metabolites available in public databases. On the basis of the number and quality of the matches, we concluded that ISA and ACP can robustly identify ten and nine metabolites, respectively, half of which were shared, while PCA did not produce any signatures with robust matches.

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.

Proof of concept for quantitative urine NMR metabolomics pipeline for large-scale epidemiology and genetics

BACKGROUND

Quantitative molecular data from urine are rare in epidemiology and genetics. NMR spectroscopy could provide these data in high throughput, and it has already been applied in epidemiological settings to analyse urine samples. However, quantitative protocols for large-scale applications are not available.

METHODS

We describe in detail how to prepare urine samples and perform NMR experiments to obtain quantitative metabolic information. Semi-automated quantitative line shape fitting analyses were set up for 43 metabolites and applied to data from various analytical test samples and from 1004 individuals from a population-based epidemiological cohort. Novel analyses on how urine metabolites associate with quantitative serum NMR metabolomics data (61 metabolic measures; n = 995) were performed. In addition, confirmatory genome-wide analyses of urine metabolites were conducted (n = 578). The fully automated quantitative regression-based spectral analysis is demonstrated for creatinine and glucose (n = 4548).

RESULTS

Intra-assay metabolite variations were mostly <5%, indicating high robustness and accuracy of urine NMR spectroscopy methodology per se. Intra-individual metabolite variations were large, ranging from 6% to 194%. However, population-based inter-individual metabolite variations were even larger (from 14% to 1655%), providing a sound base for epidemiological applications. Metabolic associations between urine and serum were found to be clearly weaker than those within serum and within urine, indicating that urinary metabolomics data provide independent metabolic information. Two previous genome-wide hits for formate and 2-hydroxyisobutyrate were replicated at genome-wide significance.

CONCLUSION

Quantitative urine metabolomics data suggest broad novelty for systems epidemiology. A roadmap for an open access methodology is provided.

A catalog of genetic loci associated with kidney function from analyses of a million individuals

Chronic kidney disease (CKD) is responsible for a public health burden with multi-systemic complications. Through trans-ancestry meta-analysis of genome-wide association studies of estimated glomerular filtration rate (eGFR) and independent replication (n = 1,046,070), we identified 264 associated loci (166 new). Of these, 147 were likely to be relevant for kidney function on the basis of associations with the alternative kidney function marker blood urea nitrogen (n = 416,178). Pathway and enrichment analyses, including mouse models with renal phenotypes, support the kidney as the main target organ. A genetic risk score for lower eGFR was associated with clinically diagnosed CKD in 452,264 independent individuals. Colocalization analyses of associations with eGFR among 783,978 European-ancestry individuals and gene expression across 46 human tissues, including tubulo-interstitial and glomerular kidney compartments, identified 17 genes differentially expressed in kidney. Fine-mapping highlighted missense driver variants in 11 genes and kidney-specific regulatory variants. These results provide a comprehensive priority list of molecular targets for translational research.

Integrated analysis of genomics, longitudinal metabolomics, and Alzheimer's risk factors among 1,111 cohort participants

Although Alzheimer's disease (AD) is highly heritable, genetic variants are known to be associated with AD only explain a small proportion of its heritability. Genetic factors may only convey disease risk in individuals with certain environmental exposures, suggesting that a multiomics approach could reveal underlying mechanisms contributing to complex traits, such as AD. We developed an integrated network to investigate relationships between metabolomics, genomics, and AD risk factors using Wisconsin Registry for Alzheimer's Prevention participants. Analyses included 1,111 non-Hispanic Caucasian participants with whole blood expression for 11,376 genes (imputed from dense genome-wide genotyping), 1,097 fasting plasma metabolites, and 17 AD risk factors. A subset of 155 individuals also had 364 fastings cerebral spinal fluid (CSF) metabolites. After adjusting each of these 12,854 variables for potential confounders, we developed an undirected graphical network, representing all significant pairwise correlations upon adjusting for multiple testing. There were many instances of genes being indirectly linked to AD risk factors through metabolites, suggesting that genes may influence AD risk through particular metabolites. Follow-up analyses suggested that glycine mediates the relationship between carbamoyl-phosphate synthase 1 and measures of cardiovascular and diabetes risk, including body mass index, waist-hip ratio, inflammation, and insulin resistance. Further, 38 CSF metabolites explained more than 60% of the variance of CSF levels of tau, a detrimental protein that accumulates in the brain of AD patients and is necessary for its diagnosis. These results further our understanding of underlying mechanisms contributing to AD risk while demonstrating the utility of generating and integrating multiple omics data types.

Genome-wide association analyses identify 139 loci associated with macular thickness in the UK Biobank cohort

The macula, located near the center of the retina in the human eye, is responsible for providing critical functions, such as central, sharp vision. Structural changes in the macula are associated with many ocular diseases, including age-related macular degeneration (AMD) and glaucoma. Although macular thickness is a highly heritable trait, there are no prior reported genome-wide association studies (GWASs) of it. Here we describe the first GWAS of macular thickness, which was measured by spectral-domain optical coherence tomography using 68 423 participants from the UK Biobank cohort. We identified 139 genetic loci associated with macular thickness at genome-wide significance (P < 5 × 10-8). The most significant loci were LINC00461 (P = 5.1 × 10-120), TSPAN10 (P = 1.2 × 10-118), RDH5 (P = 9.2 × 10-105) and SLC6A20 (P = 1.4 × 10-71). Results from gene expression demonstrated that these genes are highly expressed in the retina. Other hits included many previously reported AMD genes, such as NPLOC4 (P = 1.7 × 10-103), RAD51B (P = 9.1 × 10-14) and SLC16A8 (P = 1.7 × 10-8), further providing functional significance of the identified loci. Through cross-phenotype analysis, these genetic loci also exhibited pleiotropic effects with myopia, neurodegenerative diseases (e.g. Parkinson's disease, schizophrenia and Alzheimer's disease), cancer (e.g. breast, ovarian and lung cancers) and metabolic traits (e.g. body mass index, waist circumference and type 2 diabetes). Our findings provide the first insight into the genetic architecture of macular thickness and may further elucidate the pathogenesis of related ocular diseases, such as AMD.

Identification of seven novel loci associated with amino acid levels using single-variant and gene-based tests in 8545 Finnish men from the METSIM study

Comprehensive metabolite profiling captures many highly heritable traits, including amino acid levels, which are potentially sensitive biomarkers for disease pathogenesis. To better understand the contribution of genetic variation to amino acid levels, we performed single variant and gene-based tests of association between nine serum amino acids (alanine, glutamine, glycine, histidine, isoleucine, leucine, phenylalanine, tyrosine, and valine) and 16.6 million genotyped and imputed variants in 8545 non-diabetic Finnish men from the METabolic Syndrome In Men (METSIM) study with replication in Northern Finland Birth Cohort (NFBC1966). We identified five novel loci associated with amino acid levels (P = < 5×10-8): LOC157273/PPP1R3B with glycine (rs9987289, P = 2.3×10-26); ZFHX3 (chr16:73326579, minor allele frequency (MAF) = 0.42%, P = 3.6×10-9), LIPC (rs10468017, P = 1.5×10-8), and WWOX (rs9937914, P = 3.8×10-8) with alanine; and TRIB1 with tyrosine (rs28601761, P = 8×10-9). Gene-based tests identified two novel genes harboring missense variants of MAF <1% that show aggregate association with amino acid levels: PYCR1 with glycine (Pgene = 1.5×10-6) and BCAT2 with valine (Pgene = 7.4×10-7); neither gene was implicated by single variant association tests. These findings are among the first applications of gene-based tests to identify new loci for amino acid levels. In addition to the seven novel gene associations, we identified five independent signals at established amino acid loci, including two rare variant signals at GLDC (rs138640017, MAF=0.95%, Pconditional = 5.8×10-40) with glycine levels and HAL (rs141635447, MAF = 0.46%, Pconditional = 9.4×10-11) with histidine levels. Examination of all single variant association results in our data revealed a strong inverse relationship between effect size and MAF (Ptrend<0.001). These novel signals provide further insight into the molecular mechanisms of amino acid metabolism and potentially, their perturbations in disease.

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.

The metabolomic quest for a biomarker in chronic kidney disease

Chronic kidney disease (CKD) is a growing burden on people and on healthcare for which the diagnostics are niether disease-specific nor indicative of progression. Biomarkers are sought to enable clinicians to offer more appropriate patient-centred treatments, which could come to fruition by using a metabolomics approach. This mini-review highlights the current literature of metabolomics and CKD, and suggests additional factors that need to be considered in this quest for a biomarker, namely the diet and the gut microbiome, for more meaningful advances to be made.

Amino acid transporters revisited: New views in health and disease

Amino acid transporters (AATs) are membrane-bound transport proteins that mediate transfer of amino acids into and out of cells or cellular organelles. AATs have diverse functional roles ranging from neurotransmission to acid-base balance, intracellular energy metabolism, and anabolic and catabolic reactions. In cancer cells and diabetes, dysregulation of AATs leads to metabolic reprogramming, which changes intracellular amino acid levels, contributing to the pathogenesis of cancer, obesity and diabetes. Indeed, the neutral amino acid transporters (NATs) SLC7A5/LAT1 and SLC1A5/ASCT2 are likely involved in several human malignancies. However, a clinical therapy that directly targets AATs has not yet been developed. The purpose of this review is to highlight the structural and functional diversity of AATs, their diverse physiological roles in different tissues and organs, their wide-ranging implications in human diseases and the emerging strategies and tools that will be necessary to target AATs therapeutically.

Inborn-like errors of metabolism are determinants of breast cancer risk, clinical response and survival: a study of human biochemical individuality

Breast cancer remains a leading cause of morbidity and mortality worldwide yet methods for early detection remain elusive. We describe the discovery and validation of biochemical signatures measured by mass spectrometry, performed upon blood samples from patients and controls that accurately identify (>95%) the presence of clinical breast cancer. Targeted quantitative MS/MS conducted upon 1225 individuals, including patients with breast and other cancers, normal controls as well as individuals with a variety of metabolic disorders provide a biochemical phenotype that accurately identifies the presence of breast cancer and predicts response and survival following the administration of neoadjuvant chemotherapy. The metabolic changes identified are consistent with inborn-like errors of metabolism and define a continuum from normal controls to elevated risk to invasive breast cancer. Similar results were observed in other adenocarcinomas but were not found in squamous cell cancers or hematologic neoplasms. The findings describe a new early detection platform for breast cancer and support a role for pre-existing, inborn-like errors of metabolism in the process of breast carcinogenesis that may also extend to other glandular malignancies.

Statement of Significance: Findings provide a powerful tool for early detection and the assessment of prognosis in breast cancer and define a novel concept of breast carcinogenesis that characterizes malignant transformation as the clinical manifestation of underlying metabolic insufficiencies.

From lipid locus to drug target through human genomics

In the last decade, over 175 genetic loci have robustly been associated to levels of major circulating blood lipids. Most loci are specific to one or two lipids, whereas some (SUGP1, ZPR1, TRIB1, HERPUD1, and FADS1) are associated to all. While exposing the polygenic architecture of circulating lipids and the underpinnings of dyslipidaemia, these genome-wide association studies (GWAS) have provided further evidence of the critical role that lipids play in coronary heart disease (CHD) risk, as indicated by the 2.7-fold enrichment for macrophage gene expression in atherosclerotic plaques and the association of 25 loci (such as PCSK9, APOB, ABCG5-G8, KCNK5, LPL, HMGCR, NPC1L1, CETP, TRIB1, ABO, PMAIP1-MC4R, and LDLR) with CHD. These GWAS also confirmed known and commonly used therapeutic targets, including HMGCR (statins), PCSK9 (antibodies), and NPC1L1 (ezetimibe). As we head into the post-GWAS era, we offer suggestions for how to move forward beyond genetic risk loci, towards refining the biology behind the associations and identifying causal genes and therapeutic targets. Deep phenotyping through lipidomics and metabolomics will refine and increase the resolution to find causal and druggable targets, and studies aimed at demonstrating gene transcriptional and regulatory effects of lipid associated loci will further aid in identifying these targets. Thus, we argue the need for deeply phenotyped, large genetic association studies to reduce costs and failures and increase the efficiency of the drug discovery pipeline. We conjecture that in the next decade a paradigm shift will tip the balance towards a data-driven approach to therapeutic target development and the application of precision medicine where human genomics takes centre stage.

The role of genetic variation of human metabolism for BMI, mental traits and mental disorders

OBJECTIVE

The aim was to assess whether loci associated with metabolic traits also have a significant role in BMI and mental traits/disorders METHODS: We first assessed the number of single nucleotide polymorphisms (SNPs) with genome-wide significance for human metabolism (NHGRI-EBI Catalog). These 516 SNPs (216 independent loci) were looked-up in genome-wide association studies for association with body mass index (BMI) and the mental traits/disorders educational attainment, neuroticism, schizophrenia, well-being, anxiety, depressive symptoms, major depressive disorder, autism-spectrum disorder, attention-deficit/hyperactivity disorder, Alzheimer's disease, bipolar disorder, aggressive behavior, and internalizing problems. A strict significance threshold of p < 6.92 × 10-6 was based on the correction for 516 SNPs and all 14 phenotypes, a second less conservative threshold (p < 9.69 × 10-5) on the correction for the 516 SNPs only.

RESULTS

19 SNPs located in nine independent loci revealed p-values < 6.92 × 10-6; the less strict criterion was met by 41 SNPs in 24 independent loci. BMI and schizophrenia showed the most pronounced genetic overlap with human metabolism with three loci each meeting the strict significance threshold. Overall, genetic variation associated with estimated glomerular filtration rate showed up frequently; single metabolite SNPs were associated with more than one phenotype. Replications in independent samples were obtained for BMI and educational attainment.

CONCLUSIONS

Approximately 5-10% of the regions involved in the regulation of blood/urine metabolite levels seem to also play a role in BMI and mental traits/disorders and related phenotypes. If validated in metabolomic studies of the respective phenotypes, the associated blood/urine metabolites may enable novel preventive and therapeutic strategies.

The fecal metabolome as a functional readout of the gut microbiome

The human gut microbiome plays a key role in human health 1 , but 16S characterization lacks quantitative functional annotation 2 . The fecal metabolome provides a functional readout of microbial activity and can be used as an intermediate phenotype mediating host-microbiome interactions 3 . In this comprehensive description of the fecal metabolome, examining 1,116 metabolites from 786 individuals from a population-based twin study (TwinsUK), the fecal metabolome was found to be only modestly influenced by host genetics (heritability (H2) = 17.9%). One replicated locus at the NAT2 gene was associated with fecal metabolic traits. The fecal metabolome largely reflects gut microbial composition, explaining on average 67.7% (±18.8%) of its variance. It is strongly associated with visceral-fat mass, thereby illustrating potential mechanisms underlying the well-established microbial influence on abdominal obesity. Fecal metabolic profiling thus is a novel tool to explore links among microbiome composition, host phenotypes, and heritable complex traits.

New perspectives: systems medicine in cardiovascular disease

Background

Cardiovascular diseases (CVD) represent one of the most important causes of morbidity and mortality worldwide. Innovative approaches to increase the understanding of the underpinnings of CVD promise to enhance CVD risk assessment and might pave the way to tailored therapies. Within the last years, systems medicine has emerged as a novel tool to study the genetic, molecular and physiological interactions.

Conclusion

In this review, we provide an overview of the current molecular-experimental, epidemiological and bioinformatical tools applied in systems medicine in the cardiovascular field. We will discuss the status and challenges in implementing interdisciplinary systems medicine approaches in CVD.