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

Plasma proteomic characterization of motoric cognitive risk and mild cognitive impairment

INTRODUCTION

Motoric cognitive risk (MCR) is a pre-dementia syndrome characterized by mobility and cognitive dysfunction. This study conducted a proteome-wide study of MCR and compared the proteomic signatures of MCR to that of mild cognitive impairment (MCI).

METHODS

Participants were classified as MCR using a memory questionnaire and 4-meter walk. We measured 4877 plasma proteins collected during late-life and midlife. Multivariable logistic regression related each protein to late-life MCR/MCI. MCR-associated proteins were replicated internally at midlife and in an external cohort.

RESULTS

Proteome-wide analysis (n = 4076) identified 25 MCR-associated proteins. Eight of these proteins remained associated with late-life MCR when measured during midlife. Two proteins (SVEP1 and TAGLN) were externally replicated. Compared to MCI, MCR had a distinct and much stronger proteomic signature enriched for cardiometabolic and immune pathways.

DISCUSSION

Our findings highlight the divergent biology underlying two pre-dementia syndromes. Metabolic and immune dysfunction may be a primary driver of MCR.

HIGHLIGHTS

MCR is defined by concurrent cognitive and gait dysfunction. MCR protein biomarkers have key roles in cardiometabolic and vascular function. MCR biomarkers are also associated with cerebrovascular disease and dementia. MCR and MCI demonstrate overlapping but divergent proteomic signatures.

Genetic mapping of serum metabolome to chronic diseases among Han Chinese

Serum metabolites are potential regulators for chronic diseases. To explore the genetic regulation of metabolites and their roles in chronic diseases, we quantified 2,759 serum metabolites and performed genome-wide association studies (GWASs) among Han Chinese individuals. We identified 184 study-wide significant (p < 1.81 × 10-11) metabolite quantitative trait loci (metaboQTLs), 88.59% (163) of which were novel. Notably, we identified Asian-ancestry-specific metaboQTLs, including the SNP rs2296651 for taurocholic acid and taurochenodesoxycholic acid. Leveraging the GWAS for 37 clinical traits from East Asians, Mendelian randomization analyses identified 906 potential causal relationships between metabolites and clinical traits, including 27 for type 2 diabetes and 38 for coronary artery disease. Integrating genetic regulation of the transcriptome and proteome revealed putative regulators of several metabolites. In summary, we depict a landscape of the genetic architecture of the serum metabolome among Han Chinese and provide insights into the role of serum metabolites in chronic diseases.

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.

High-throughput untargeted metabolomics reveals metabolites and metabolic pathways that differentiate two divergent pig breeds

Metabolomics can describe the molecular phenome and may contribute to dissecting the biological processes linked to economically relevant traits in livestock species. Comparative analyses of metabolomic profiles in purebred pigs can provide insights into the basic biological mechanisms that may explain differences in production performances. Following this concept, this study was designed to compare, on a large scale, the plasma metabolomic profiles of two Italian heavy pig breeds (Italian Duroc and Italian Large White) to indirectly evaluate the impact of their different genetic backgrounds on the breed metabolomes. We utilised a high-throughput untargeted metabolomics approach in a total of 962 pigs that allowed us to detect and relatively quantify 722 metabolites from various biological classes. The molecular data were analysed using a bioinformatics pipeline specifically designed for identifying differentially abundant metabolites between the two breeds in a robust and statistically significant manner, including the Boruta algorithm, which is a Random Forest wrapper, and sparse Partial Least Squares Discriminant Analysis (sPLS-DA) for feature selection. After thoroughly evaluating the impact of random components on missing value imputation, 100 discriminant metabolites were selected by Boruta and 17 discriminant metabolites (all included within the previous list) were identified with sPLS-DA. About half of the 100 discriminant metabolites had a higher concentration in one or the other breed (48 in Italian Large White pigs, with a prevalence of amino acids and peptides; 52 in Italian Duroc pigs, with a prevalence of lipids). These metabolites were from seven distinct super pathways and had an absolute mean value of percentage difference between the two breeds (|Δ|%) of 39.2 ± 32.4. Six of these metabolites had |Δ|%> 100. A general correlation network analysis based on Boruta-identified metabolites consisted of 31 singletons and 69 metabolites connected by 141 edges, with two large clusters (> 15 nodes), three medium clusters (3-6 nodes) and eight additional pairs, with most metabolites belonging to the same super pathway. The major cluster representing the lipids super-pathway included 24 metabolites, primarily sphingomyelins. Overall, this study identified metabolomic differences between Italian Duroc and Italian Large White pigs explained by the specific genetic background of the two breeds. These biomarkers can explain the biological differences between these two breeds and can have potential practical applications in pig breeding and husbandry.

Coupling metabolomics and exome sequencing reveals graded effects of rare damaging heterozygous variants on gene function and human traits

Genetic studies of the metabolome can uncover enzymatic and transport processes shaping human metabolism. Using rare variant aggregation testing based on whole-exome sequencing data to detect genes associated with levels of 1,294 plasma and 1,396 urine metabolites, we discovered 235 gene-metabolite associations, many previously unreported. Complementary approaches (genetic, computational (in silico gene knockouts in whole-body models of human metabolism) and one experimental proof of principle) provided orthogonal evidence that studies of rare, damaging variants in the heterozygous state permit inferences concordant with those from inborn errors of metabolism. Allelic series of functional variants in transporters responsible for transcellular sulfate reabsorption (SLC13A1, SLC26A1) exhibited graded effects on plasma sulfate and human height and pinpointed alleles associated with increased odds of diverse musculoskeletal traits and diseases in the population. This integrative approach can identify new players in incompletely characterized human metabolic reactions and reveal metabolic readouts informative of human traits and diseases.

Genetic architecture of cerebrospinal fluid and brain metabolite levels and the genetic colocalization of metabolites with human traits

Brain metabolism perturbation can contribute to traits and diseases. We conducted a genome-wide association study for cerebrospinal fluid (CSF) and brain metabolite levels, identifying 205 independent associations (47.3% new signals, containing 11 new loci) for 139 CSF metabolites, and 32 independent associations (43.8% new signals, containing 4 new loci) for 31 brain metabolites. Of these, 96.9% (CSF) and 71.4% (brain) of the new signals belonged to previously analyzed metabolites in blood or urine. We integrated the metabolite quantitative trait loci (MQTLs) with 23 neurological, psychiatric and common human traits and diseases through colocalization to identify metabolites and biological processes implicated in these phenotypes. Combining CSF and brain, we identified 71 metabolite-trait associations, such as glycerophosphocholines with Alzheimer's disease, O-sulfo-L-tyrosine with Parkinson's disease, glycine, xanthine with waist-to-hip ratio and ergothioneine with inflammatory bowel disease. Our study expanded the knowledge of MQTLs in the central nervous system, providing insights into human traits.

Plasma proteome-wide analysis of cerebral small vessel disease identifies novel biomarkers and disease pathways

Cerebral small vessel disease (SVD), as defined by neuroimaging characteristics such as white matter hyperintensities (WMHs), cerebral microhemorrhages (CMHs), and lacunar infarcts, is highly prevalent and has been associated with dementia risk and other clinical sequelae. Although conditions such as hypertension are known to contribute to SVD, little is known about the diverse set of subclinical biological processes and molecular mediators that may also influence the development and progression of SVD. To better understand the mechanisms underlying SVD and to identify novel SVD biomarkers, we used a large-scale proteomic platform to relate 4,877 plasma proteins to MRI-defined SVD characteristics within 1,508 participants of the Atherosclerosis Risk in Communities (ARIC) Study cohort. Our proteome-wide analysis of older adults (mean age: 76) identified 13 WMH-associated plasma proteins involved in synaptic function, endothelial integrity, and angiogenesis, two of which remained associated with late-life WMH volume when measured nearly 20 years earlier, during midlife. We replicated the relationship between 9 candidate proteins and WMH volume in one or more external cohorts; we found that 11 of the 13 proteins were associated with risk for future dementia; and we leveraged publicly available proteomic data from brain tissue to demonstrate that a subset of WMH-associated proteins was differentially expressed in the context of cerebral atherosclerosis, pathologically-defined Alzheimer's disease, and cognitive decline. Bidirectional two-sample Mendelian randomization analyses examined the causal relationships between candidate proteins and WMH volume, while pathway and network analyses identified discrete biological processes (lipid/cholesterol metabolism, NF-kB signaling, hemostasis) associated with distinct forms of SVD. Finally, we synthesized these findings to identify two plasma proteins, oligodendrocyte myelin glycoprotein (OMG) and neuronal pentraxin receptor (NPTXR), as top candidate biomarkers for elevated WMH volume and its clinical manifestations.

Withdrawn: Combinatorial lipidomics and proteomics underscore erythrocyte lipid membrane aberrations in the development of adverse cardio-cerebrovascular complications in maintenance hemodialysis patients

This article has been withdrawn: please see Elsevier Policy on Article Withdrawal (https://www.elsevier.com/about/policies/article-withdrawal). The authors reached out to the Publisher to alert the Publisher to incorrect text published in the article. After investigating the situation, the journal came to the conclusion that the wrong version of the file was sent by the authors to the production team during the proof stage and the misplaced text was not noticed by the authors when they approved the final version. After consulting with the Editor-in-Chief of the journal, the decision was made to withdraw the current version of the article.

The HuMet Repository: Watching human metabolism at work

Metabolism oscillates between catabolic and anabolic states depending on food intake, exercise, or stresses that change a multitude of metabolic pathways simultaneously. We present the HuMet Repository for exploring dynamic metabolic responses to oral glucose/lipid loads, mixed meals, 36-h fasting, exercise, and cold stress in healthy subjects. Metabolomics data from blood, urine, and breath of 15 young, healthy men at up to 56 time points are integrated and embedded within an interactive web application, enabling researchers with and without computational expertise to search, visualize, analyze, and contextualize the dynamic metabolite profiles of 2,656 metabolites acquired on multiple platforms. With examples, we demonstrate the utility of the resource for research into the dynamics of human metabolism, highlighting differences and similarities in systemic metabolic responses across challenges and the complementarity of metabolomics platforms. The repository, providing a reference for healthy metabolite changes to six standardized physiological challenges, is freely accessible through a web portal.

SLC17A1/3 transporters mediate renal excretion of Lac-Phe in mice and humans

N-lactoyl-phenylalanine (Lac-Phe) is a lactate-derived metabolite that suppresses food intake and body weight. Little is known about the mechanisms that mediate Lac-Phe transport across cell membranes. Here we identify SLC17A1 and SLC17A3, two kidney-restricted plasma membrane-localized solute carriers, as physiologic urine Lac-Phe transporters. In cell culture, SLC17A1/3 exhibit high Lac-Phe efflux activity. In humans, levels of Lac-Phe in urine exhibit a strong genetic association with the SLC17A1-4 locus. Urine Lac-Phe levels are increased following a Wingate sprint test. In mice, genetic ablation of either SLC17A1 or SLC17A3 reduces urine Lac-Phe levels. Despite these differences, both knockout strains have normal blood Lac-Phe and body weights, demonstrating SLC17A1/3-dependent de-coupling of urine and plasma Lac-Phe pools. Together, these data establish SLC17A1/3 family members as the physiologic urine Lac-Phe transporters and uncover a biochemical pathway for the renal excretion of this signaling metabolite.

Systematic investigation of genetically determined plasma and urinary metabolites to discover potential interventional targets for colorectal cancer

BACKGROUND

We aimed to identify plasma and urinary metabolites related to colorectal cancer (CRC) risk and elucidate their mediator role in the associations between modifiable risk factors and CRC.

METHODS

Metabolite quantitative trait loci were derived from 2 published metabolomics genome-wide association studies, and summary-level data were extracted for 651 plasma metabolites and 208 urinary metabolites. Genetic associations with CRC were obtained from a large-scale genome-wide association study meta-analysis (100 204 cases, 154 587 controls) and the FinnGen cohort (4957 cases, 304 197 controls). Mendelian randomization and colocalization analyses were performed to evaluate the causal roles of metabolites in CRC. Druggability evaluation was employed to prioritize potential therapeutic targets. Multivariable Mendelian randomization and mediation estimation were conducted to elucidate the mediating effects of metabolites on the associations between modifiable risk factors and CRC.

RESULTS

The study identified 30 plasma metabolites and 4 urinary metabolites for CRC. Plasma sphingomyelin and urinary lactose, which were positively associated with CRC risk, could be modulated by drug interventions (ie, olipudase alfa, tilactase). Thirteen modifiable risk factors were associated with 9 metabolites, and 8 of these modifiable risk factors were associated with CRC risk. These 9 metabolites mediated the effect of modifiable risk factors (Actinobacteria, body mass index, waist to hip ratio, fasting insulin, smoking initiation) on CRC.

CONCLUSION

This study identified key metabolite biomarkers associated with CRC and elucidated their mediator roles in the associations between modifiable risk factors and CRC. These findings provide new insights into the etiology and potential therapeutic targets for CRC and the etiological pathways of modifiable environmental factors with CRC.

Common genetic variants associated with urinary phthalate levels in children: A genome-wide study

INTRODUCTION

Phthalates, or dieters of phthalic acid, are a ubiquitous type of plasticizer used in a variety of common consumer and industrial products. They act as endocrine disruptors and are associated with increased risk for several diseases. Once in the body, phthalates are metabolized through partially known mechanisms, involving phase I and phase II enzymes.

OBJECTIVE

In this study we aimed to identify common single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) associated with the metabolism of phthalate compounds in children through genome-wide association studies (GWAS).

METHODS

The study used data from 1,044 children with European ancestry from the Human Early Life Exposome (HELIX) cohort. Ten phthalate metabolites were assessed in a two-void pooled urine collected at the mean age of 8 years. Six ratios between secondary and primary phthalate metabolites were calculated. Genome-wide genotyping was done with the Infinium Global Screening Array (GSA) and imputation with the Haplotype Reference Consortium (HRC) panel. PennCNV was used to estimate copy number variants (CNVs) and CNVRanger to identify consensus regions. GWAS of SNPs and CNVs were conducted using PLINK and SNPassoc, respectively. Subsequently, functional annotation of suggestive SNPs (p-value < 1E-05) was done with the FUMA web-tool.

RESULTS

We identified four genome-wide significant (p-value < 5E-08) loci at chromosome (chr) 3 (FECHP1 for oxo-MiNP_oh-MiNP ratio), chr6 (SLC17A1 for MECPP_MEHHP ratio), chr9 (RAPGEF1 for MBzP), and chr10 (CYP2C9 for MECPP_MEHHP ratio). Moreover, 115 additional loci were found at suggestive significance (p-value < 1E-05). Two CNVs located at chr11 (MRGPRX1 for oh-MiNP and SLC35F2 for MEP) were also identified. Functional annotation pointed to genes involved in phase I and phase II detoxification, molecular transfer across membranes, and renal excretion.

CONCLUSION

Through genome-wide screenings we identified known and novel loci implicated in phthalate metabolism in children. Genes annotated to these loci participate in detoxification, transmembrane transfer, and renal excretion.

European and African-specific plasma protein-QTL and metabolite-QTL analyses identify ancestry-specific T2D effector proteins and metabolites

Initially focused on the European population, multiple genome-wide association studies (GWAS) of complex diseases, such as type-2 diabetes (T2D), have now extended to other populations. However, to date, few ancestry-matched omics datasets have been generated or further integrated with the disease GWAS to nominate the key genes and/or molecular traits underlying the disease risk loci. In this study, we generated and integrated plasma proteomics and metabolomics with array-based genotype datasets of European (EUR) and African (AFR) ancestries to identify ancestry-specific muti-omics quantitative trait loci (QTLs). We further applied these QTLs to ancestry-stratified T2D risk to pinpoint key proteins and metabolites underlying the disease-associated genetic loci. We nominated five proteins and four metabolites in the European group and one protein and one metabolite in the African group to be part of the molecular pathways of T2D risk in an ancestry-stratified manner. Our study demonstrates the integration of genetic and omic studies of different ancestries can be used to identify distinct effector molecular traits underlying the same disease across diverse populations. Specifically, in the AFR proteomic findings on T2D, we prioritized the protein QSOX2; while in the AFR metabolomic findings, we pinpointed the metabolite GlcNAc sulfate conjugate of C21H34O2 steroid. Neither of these findings overlapped with the corresponding EUR results.

Association of Integrated Proteomic and Metabolomic Modules with Risk of Kidney Disease Progression

Key Points

Integrated analysis of proteome and metabolome identifies modules associated with CKD progression and kidney failure. Ephrin transmembrane proteins and podocyte-expressed CRIM1 and NPNT emerged as central components and warrant experimental and clinical investigation.

Background

Proteins and metabolites play crucial roles in various biological functions and are frequently interconnected through enzymatic or transport processes.

Methods

We present an integrated analysis of 4091 proteins and 630 metabolites in the Chronic Renal Insufficiency Cohort study (N =1708; average follow-up for kidney failure, 9.5 years, with 537 events). Proteins and metabolites were integrated using an unsupervised clustering method, and we assessed associations between clusters and CKD progression and kidney failure using Cox proportional hazards models. Analyses were adjusted for demographics and risk factors, including the eGFR and urine protein–creatinine ratio. Associations were identified in a discovery sample (random two thirds, n =1139) and then evaluated in a replication sample (one third, n =569).

Results

We identified 139 modules of correlated proteins and metabolites, which were represented by their principal components. Modules and principal component loadings were projected onto the replication sample, which demonstrated a consistent network structure. Two modules, representing a total of 236 proteins and 82 metabolites, were robustly associated with both CKD progression and kidney failure in both discovery and validation samples. Using gene set enrichment, several transmembrane-related terms were identified as overrepresented in these modules. Transmembrane–ephrin receptor activity displayed the largest odds (odds ratio=13.2, P value = 5.5×10−5). A module containing CRIM1 and NPNT expressed in podocytes demonstrated particularly strong associations with kidney failure (P value = 2.6×10−5).

Conclusions

This study demonstrates that integration of the proteome and metabolome can identify functions of pathophysiologic importance in kidney disease.

The discovery of GGT1 as a novel gene for ischemic stroke conferring protection against disease risk in non-smokers and non-abusers of alcohol

OBJECTIVES

Increased plasma gamma-glutamyl transferase (GGT1) has been identified as a robust and independent risk factor for ischemic stroke (IS), but the molecular mechanisms of the enzyme-disease association are unclear. The present study investigated whether polymorphisms in the GGT1 gene contribute to IS susceptibility.

MATERIALS AND METHODS

DNA samples obtained from 1288 unrelated individuals (600 IS patients and 688 controls) were genotyped for common single nucleotide polymorphisms of GGT1 using the MassArray-4 platform.

RESULTS

The rs5751909 polymorphism was significantly associated with decreased risk of ischemic stroke regardless sex and age (Pperm ≤ 0.01, dominant genetic model). The haplotype rs4820599A-rs5760489A-rs5751909A showed strong protection against ischemic stroke (OR 0.53, 95 %CI 0.36 - 0.77, Pperm ≤ 0.0001). The protective effect of SNP rs5751909 in the stroke phenotype was successfully replicated in the UK Biobank, SiGN, and ISGC cohorts (P ≤ 0.01). GGT1 polymorphisms showed joint (epistatic) effects on the risk of ischemic stroke, with some known IS-associated GWAS loci (e.g., rs4322086 and rs12646447) investigated in our population. In addition, SNP rs5751909 was found to be strongly associated with a decreased risk of ischemic stroke in non-smokers (OR 0.54 95 %CI 0.39-0.75, Pperm = 0.0002) and non-alcohol abusers (OR 0.43 95 %CI 0.30-0.61, Pperm = 2.0 × 10-6), whereas no protective effects of this SNP against disease risk were observed in smokers and alcohol abusers (Pperm < 0.05).

CONCLUSIONS

We propose mechanisms underlying the observed associations between GGT1 polymorphisms and ischemic stroke risk. This pilot study is the first to demonstrate that GGT1 is a novel susceptibility gene for ischemic stroke and provides additional evidence of the genetic contribution to impaired redox homeostasis underlying disease pathogenesis.

SLC17 transporters mediate renal excretion of Lac-Phe in mice and humans

N-lactoyl-phenylalanine (Lac-Phe) is a lactate-derived metabolite that suppresses food intake and body weight. Little is known about the mechanisms that mediate Lac-Phe transport across cell membranes. Here we identify SLC17A1 and SLC17A3, two kidney-restricted plasma membrane-localized solute carriers, as physiologic urine Lac-Phe transporters. In cell culture, SLC17A1/3 exhibit high Lac-Phe efflux activity. In humans, levels of Lac-Phe in urine exhibit a strong genetic association with the SLC17A1-4 locus. Urine Lac-Phe levels are also increased following a Wingate sprint test. In mice, genetic ablation of either SLC17A1 or SLC17A3 reduces urine Lac-Phe levels. Despite these differences, both knockout strains have normal blood Lac-Phe and body weights, demonstrating that urine and plasma Lac-Phe pools are functionally and biochemically de-coupled. Together, these data establish SLC17 family members as the physiologic urine transporters for Lac-Phe and uncover a biochemical pathway for the renal excretion of this signaling metabolite.

A systematic review of metabolomic findings in adult and pediatric renal disease

Chronic kidney disease (CKD) affects over 0.5 billion people worldwide across their lifetimes. Despite a growingly ageing world population, an increase in all-age prevalence of kidney disease persists. Adult-onset forms of kidney disease often result from lifestyle-modifiable metabolic illnesses such as type 2 diabetes. Pediatric and adolescent forms of renal disease are primarily caused by morphological abnormalities of the kidney, as well as immunological, infectious and inherited metabolic disorders. Alterations in energy metabolism are observed in CKD of varying causes, albeit the molecular mechanisms underlying pathology are unclear. A systematic indexing of metabolites identified in plasma and urine of patients with kidney disease alongside disease enrichment analysis uncovered inborn errors of metabolism as a framework that links features of adult and pediatric kidney disease. The relationship of genetics and metabolism in kidney disease could be classified into three distinct landscapes: (i) Normal genotypes that develop renal damage because of lifestyle and / or comorbidities; (ii) Heterozygous genetic variants and polymorphisms that result in unique metabotypes that may predispose to the development of kidney disease via synergistic heterozygosity, and (iii) Homozygous genetic variants that cause renal impairment by perturbing metabolism, as found in children with monogenic inborn errors of metabolism. Interest in the identification of early biomarkers of onset and progression of CKD has grown steadily in the last years, though it has not translated into clinical routine yet. This systematic review indexes findings of differential concentration of metabolites and energy pathway dysregulation in kidney disease and appraises their potential use as biomarkers.

The Physiological and Pathological Role of Acyl-CoA Oxidation

Fatty acid metabolism, including β-oxidation (βOX), plays an important role in human physiology and pathology. βOX is an essential process in the energy metabolism of most human cells. Moreover, βOX is also the source of acetyl-CoA, the substrate for (a) ketone bodies synthesis, (b) cholesterol synthesis, (c) phase II detoxication, (d) protein acetylation, and (d) the synthesis of many other compounds, including N-acetylglutamate-an important regulator of urea synthesis. This review describes the current knowledge on the importance of the mitochondrial and peroxisomal βOX in various organs, including the liver, heart, kidney, lung, gastrointestinal tract, peripheral white blood cells, and other cells. In addition, the diseases associated with a disturbance of fatty acid oxidation (FAO) in the liver, heart, kidney, lung, alimentary tract, and other organs or cells are presented. Special attention was paid to abnormalities of FAO in cancer cells and the diseases caused by mutations in gene-encoding enzymes involved in FAO. Finally, issues related to α- and ω- fatty acid oxidation are discussed.

Quantitative trait loci mapping of circulating metabolites in cerebrospinal fluid to uncover biological mechanisms involved in brain-related phenotypes

Genomic studies of molecular traits have provided mechanistic insights into complex disease, though these lag behind for brain-related traits due to the inaccessibility of brain tissue. We leveraged cerebrospinal fluid (CSF) to study neurobiological mechanisms in vivo , measuring 5,543 CSF metabolites, the largest panel in CSF to date, in 977 individuals of European ancestry. Individuals originated from two separate cohorts including cognitively healthy subjects (n=490) and a well-characterized memory clinic sample, the Amsterdam Dementia Cohort (ADC, n=487). We performed metabolite quantitative trait loci (mQTL) mapping on CSF metabolomics and found 126 significant mQTLs, representing 65 unique CSF metabolites across 51 independent loci. To better understand the role of CSF mQTLs in brain-related disorders, we performed a metabolome-wide association study (MWAS), identifying 40 associations between CSF metabolites and brain traits. Similarly, over 90% of significant mQTLs demonstrated colocalized associations with brain-specific gene expression, unveiling potential neurobiological pathways.

Linkage analysis using whole exome sequencing data implicates SLC17A1, SLC17A3, TATDN2 and TMEM131L in type 1 diabetes in Kuwaiti families

Type 1 diabetes (T1D) is characterized by the progressive destruction of pancreatic β-cells, leading to insulin deficiency and lifelong dependency on exogenous insulin. Higher estimates of heritability rates in monozygotic twins, followed by dizygotic twins and sib-pairs, indicate the role of genetics in the pathogenesis of T1D. The incidence and prevalence of T1D are alarmingly high in Kuwait. Consanguineous marriages account for 50-70% of all marriages in Kuwait, leading to an excessive burden of recessive allele enrichment and clustering of familial disorders. Thus, genetic studies from this Arab region are expected to lead to the identification of novel gene loci for T1D. In this study, we performed linkage analyses to identify the recurrent genetic variants segregating in high-risk Kuwaiti families with T1D. We studied 18 unrelated Kuwaiti native T1D families using whole exome sequencing data from 86 individuals, of whom 37 were diagnosed with T1D. The study identified three potential loci with a LOD score of ≥ 3, spanning across four candidate genes, namely SLC17A1 (rs1165196:pT269I), SLC17A3 (rs942379: p.S370S), TATDN2 (rs394558:p.V256I), and TMEM131L (rs6848033:p.R190R). Upon examination of missense variants from these genes in the familial T1D dataset, we observed a significantly increased enrichment of the genotype homozygous for the minor allele at SLC17A3 rs56027330_p.G279R accounting for 16.2% in affected children from 6 unrelated Kuwaiti T1D families compared to 1000 genomes Phase 3 data (0.9%). Data from the NephQTL database revealed that the rs1165196, rs942379, rs394558, and rs56027330 SNPs exhibited genotype-based differential expression in either glomerular or tubular tissues. Data from the GTEx database revealed rs942379 and rs394558 as QTL variants altering the expression of TRIM38 and IRAK2 respectively. Global genome-wide association studies indicated that SLC17A1 rs1165196 and other variants from SLC17A3 are associated with uric acid concentrations and gout. Further evidence from the T1D Knowledge portal supported the role of shortlisted variants in T1D pathogenesis and urate metabolism. Our study suggests the involvement of SLC17A1, SLC17A3, TATDN2, and TMEM131L genes in familial T1D in Kuwait. An enrichment selection of genotype homozygous for the minor allele is observed at SLC17A3 rs56027330_p.G279R variant in affected members of Kuwaiti T1D families. Future studies may focus on replicating the findings in a larger T1D cohort and delineate the mechanistic details of the impact of these novel candidate genes on the pathophysiology of T1D.

metGWAS 1.0: an R workflow for network-driven over-representation analysis between independent metabolomic and meta-genome-wide association studies

MOTIVATION

The method of genome-wide association studies (GWAS) and metabolomics combined provide an quantitative approach to pinpoint metabolic pathways and genes linked to specific diseases; however, such analyses require both genomics and metabolomics datasets from the same individuals/samples. In most cases, this approach is not feasible due to high costs, lack of technical infrastructure, unavailability of samples, and other factors. Therefore, an unmet need exists for a bioinformatics tool that can identify gene loci-associated polymorphic variants for metabolite alterations seen in disease states using standalone metabolomics.

RESULTS

Here, we developed a bioinformatics tool, metGWAS 1.0, that integrates independent GWAS data from the GWAS database and standalone metabolomics data using a network-based systems biology approach to identify novel disease/trait-specific metabolite-gene associations. The tool was evaluated using standalone metabolomics datasets extracted from two metabolomics-GWAS case studies. It discovered both the observed and novel gene loci with known single nucleotide polymorphisms when compared to the original studies.

AVAILABILITY AND IMPLEMENTATION

The developed metGWAS 1.0 framework is implemented in an R pipeline and available at: https://github.com/saifurbd28/metGWAS-1.0.

Genome-Wide Association Study of CKD Progression

SIGNIFICANCE STATEMENT

Rapid progression of CKD is associated with poor clinical outcomes. Most previous studies looking for genetic factors associated with low eGFR have used cross-sectional data. The authors conducted a meta-analysis of genome-wide association studies of eGFR decline among 116,870 participants with CKD, focusing on longitudinal data. They identified three loci (two of them novel) associated with longitudinal eGFR decline. In addition to the known UMOD/PDILT locus, variants within BICC1 were associated with significant differences in longitudinal eGFR slope. Variants within HEATR4 also were associated with differences in eGFR decline, but only among Black/African American individuals without diabetes. These findings help characterize molecular mechanisms of eGFR decline in CKD and may inform new therapeutic approaches for progressive kidney disease.

BACKGROUND

Rapid progression of CKD is associated with poor clinical outcomes. Despite extensive study of the genetics of cross-sectional eGFR, only a few loci associated with eGFR decline over time have been identified.

METHODS

We performed a meta-analysis of genome-wide association studies of eGFR decline among 116,870 participants with CKD-defined by two outpatient eGFR measurements of <60 ml/min per 1.73 m 2 , obtained 90-365 days apart-from the Million Veteran Program and Vanderbilt University Medical Center's DNA biobank. The primary outcome was the annualized relative slope in outpatient eGFR. Analyses were stratified by ethnicity and diabetes status and meta-analyzed thereafter.

RESULTS

In cross-ancestry meta-analysis, the strongest association was rs77924615, near UMOD / PDILT ; each copy of the G allele was associated with a 0.30%/yr faster eGFR decline ( P = 4.9×10 -27 ). We also observed an association within BICC1 (rs11592748), where every additional minor allele was associated with a 0.13%/yr slower eGFR decline ( P = 5.6×10 -9 ). Among participants without diabetes, the strongest association was the UMOD/PDILT variant rs36060036, associated with a 0.27%/yr faster eGFR decline per copy of the C allele ( P = 1.9×10 -17 ). Among Black participants, a significantly faster eGFR decline was associated with variant rs16996674 near APOL1 (R 2 =0.29 with the G1 high-risk genotype); among Black participants with diabetes, lead variant rs11624911 near HEATR4 also was associated with a significantly faster eGFR decline. We also nominally replicated loci with known associations with eGFR decline, near PRKAG2, FGF5, and C15ORF54.

CONCLUSIONS

Three loci were significantly associated with longitudinal eGFR change at genome-wide significance. These findings help characterize molecular mechanisms of eGFR decline and may contribute to the development of new therapeutic approaches for progressive CKD.

Targeted Metabolomics Analysis of Individuals Carrying the ANGPTL8 R59W Variant

ANGPTL8 is recognized as a regulator of lipid metabolism through its role in inhibiting lipoprotein lipase activity. ANGPTL8 gene variants, particularly rs2278426 leading to the R59W variant in the protein, have been associated with lipid traits in various ethnicities. We aimed to use metabolomics to understand the impact of the ANGPTL8 R59W variant on metabolites in humans. We used the Biocrates-p400 kit to quantify 408 plasma metabolites in 60 adult male Arab individuals from Kuwait and identify differences in metabolite levels between individuals carrying reference genotypes and those with carrier genotypes at ANGPTL8 rs2278426. Individuals with carrier genotypes (CT+TT) compared to those carrying the reference genotype (CC) showed statistically significant differences in the following metabolites: acylcarnitine (perturbs metabolic pathways), phosphatidylcholine (supports liver function and cholesterol levels), cholesteryl ester (brings chronic inflammatory response to lipoprotein depositions in arteries), α-aminoadipic acid (modulates glucose homeostasis), histamine (regulates glucose/lipid metabolism), sarcosine (links amino acid and lipid metabolism), diacylglycerol 42:1 (regulates homeostasis of cellular lipid stores), and lysophosphatidylcholine (regulates oxidative stress and inflammatory response). Functional aspects attributed to these metabolites indicate that the ANGPTL8 R59W variant influences the concentrations of lipid- and inflammation-related metabolites. This observation further highlights the role of ANGPTL8 in lipid metabolism.

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.

A genome-wide association study identifies 41 loci associated with eicosanoid levels

Eicosanoids are biologically active derivatives of polyunsaturated fatty acids with broad relevance to health and disease. We report a genome-wide association study in 8406 participants of the Atherosclerosis Risk in Communities Study, identifying 41 loci associated with 92 eicosanoids and related metabolites. These findings highlight loci required for eicosanoid biosynthesis, including FADS1-3, ELOVL2, and numerous CYP450 loci. In addition, significant associations implicate a range of non-oxidative lipid metabolic processes in eicosanoid regulation, including at PKD2L1/SCD and several loci involved in fatty acyl-CoA metabolism. Further, our findings highlight select clearance mechanisms, for example, through the hepatic transporter encoded by SLCO1B1. Finally, we identify eicosanoids associated with aspirin and non-steroidal anti-inflammatory drug use and demonstrate the substantial impact of genetic variants even for medication-associated eicosanoids. These findings shed light on both known and unknown aspects of eicosanoid metabolism and motivate interest in several gene-eicosanoid associations as potential functional participants in human disease.

Proteomics analysis of plasma from middle-aged adults identifies protein markers of dementia risk in later life

A diverse set of biological processes have been implicated in the pathophysiology of Alzheimer's disease (AD) and related dementias. However, there is limited understanding of the peripheral biological mechanisms relevant in the earliest phases of the disease. Here, we used a large-scale proteomics platform to examine the association of 4877 plasma proteins with 25-year dementia risk in 10,981 middle-aged adults. We found 32 dementia-associated plasma proteins that were involved in proteostasis, immunity, synaptic function, and extracellular matrix organization. We then replicated the association between 15 of these proteins and clinically relevant neurocognitive outcomes in two independent cohorts. We demonstrated that 12 of these 32 dementia-associated proteins were associated with cerebrospinal fluid (CSF) biomarkers of AD, neurodegeneration, or neuroinflammation. We found that eight of these candidate protein markers were abnormally expressed in human postmortem brain tissue from patients with AD, although some of the proteins that were most strongly associated with dementia risk, such as GDF15, were not detected in these brain tissue samples. Using network analyses, we found a protein signature for dementia risk that was characterized by dysregulation of specific immune and proteostasis/autophagy pathways in adults in midlife ~20 years before dementia onset, as well as abnormal coagulation and complement signaling ~10 years before dementia onset. Bidirectional two-sample Mendelian randomization genetically validated nine of our candidate proteins as markers of AD in midlife and inferred causality of SERPINA3 in AD pathogenesis. Last, we prioritized a set of candidate markers for AD and dementia risk prediction in midlife.

Integrating genetics and metabolomics from multi-ethnic and multi-fluid data reveals putative mechanisms for age-related macular degeneration

Age-related macular degeneration (AMD) is a leading cause of blindness in older adults. Investigating shared genetic components between metabolites and AMD can enhance our understanding of its pathogenesis. We conduct metabolite genome-wide association studies (mGWASs) using multi-ethnic genetic and metabolomic data from up to 28,000 participants. With bidirectional Mendelian randomization analysis involving 16,144 advanced AMD cases and 17,832 controls, we identify 108 putatively causal relationships between plasma metabolites and advanced AMD. These metabolites are enriched in glycerophospholipid metabolism, lysophospholipid, triradylcglycerol, and long chain polyunsaturated fatty acid pathways. Bayesian genetic colocalization analysis and a customized metabolome-wide association approach prioritize putative causal AMD-associated metabolites. We find limited evidence linking urine metabolites to AMD risk. Our study emphasizes the contribution of plasma metabolites, particularly lipid-related pathways and genes, to AMD risk and uncovers numerous putative causal associations between metabolites and AMD risk.

Genome-Wide Association Studies for Albuminuria of Nondiabetic Taiwanese Population

INTRODUCTION

Chronic kidney disease, which is defined by a reduced estimated glomerular filtration rate and albuminuria, imposes a large health burden worldwide. Ethnicity-specific associations are frequently observed in genome-wide association studies (GWAS). This study conducts a GWAS of albuminuria in the nondiabetic population of Taiwan.

METHODS

Nondiabetic individuals aged 30-70 years without a history of cancer were enrolled from the Taiwan Biobank. A total of 6,768 subjects were subjected to a spot urine examination. After quality control using PLINK and imputation using SHAPEIT and IMPUTE2, a total of 3,638,350 single-nucleotide polymorphisms (SNPs) remained for testing. SNPs with a minor allele frequency of less than 0.1% were excluded. Linear regression was used to determine the relationship between SNPs and log urine albumin-to-creatinine ratio.

RESULTS

Six suggestive loci are identified in or near the FCRL3 (p = 2.56 × 10-6), TMEM161 (p = 4.43 × 10-6), EFCAB1 (p = 2.03 × 10-6), ELMOD1 (p = 2.97 × 10-6), RYR3 (p = 1.34 × 10-6), and PIEZO2 (p = 2.19 × 10-7). Genetic variants in the FCRL3 gene that encode a secretory IgA receptor are found to be associated with IgA nephropathy, which can manifest as proteinuria. The PIEZO2 gene encodes a sensor for mechanical forces in mesangial cells and renin-producing cells. Five SNPs with a p-value between 5 × 10-6 and 5 × 10-5 are also identified in five genes that may have a biological role in the development of albuminuria.

CONCLUSION

Five new loci and one known suggestive locus for albuminuria are identified in the nondiabetic Taiwanese population.

Transcriptome- and proteome-wide association studies nominate determinants of kidney function and damage

BACKGROUND

The pathophysiological causes of kidney disease are not fully understood. Here we show that the integration of genome-wide genetic, transcriptomic, and proteomic association studies can nominate causal determinants of kidney function and damage.

RESULTS

Through transcriptome-wide association studies (TWAS) in kidney cortex, kidney tubule, liver, and whole blood and proteome-wide association studies (PWAS) in plasma, we assess for effects of 12,893 genes and 1342 proteins on kidney filtration (glomerular filtration rate (GFR) estimated by creatinine; GFR estimated by cystatin C; and blood urea nitrogen) and kidney damage (albuminuria). We find 1561 associations distributed among 260 genomic regions that are supported as putatively causal. We then prioritize 153 of these genomic regions using additional colocalization analyses. Our genome-wide findings are supported by existing knowledge (animal models for MANBA, DACH1, SH3YL1, INHBB), exceed the underlying GWAS signals (28 region-trait combinations without significant GWAS hit), identify independent gene/protein-trait associations within the same genomic region (INHBC, SPRYD4), nominate tissues underlying the associations (tubule expression of NRBP1), and distinguish markers of kidney filtration from those with a role in creatinine and cystatin C metabolism. Furthermore, we follow up on members of the TGF-beta superfamily of proteins and find a prognostic value of INHBC for kidney disease progression even after adjustment for measured glomerular filtration rate (GFR).

CONCLUSION

In summary, this study combines multimodal, genome-wide association studies to generate a catalog of putatively causal target genes and proteins relevant to kidney function and damage which can guide follow-up studies in physiology, basic science, and clinical medicine.

Unique genetic architecture of CSF and brain metabolites pinpoints the novel targets for the traits of human wellness

Brain metabolism perturbation can contribute to traits and diseases. We conducted the first large-scale CSF and brain genome-wide association studies, which identified 219 independent associations (59.8% novel) for 144 CSF metabolites and 36 independent associations (55.6% novel) for 34 brain metabolites. Most of the novel signals (97.7% and 70.0% in CSF and brain) were tissue specific. We also integrated MWAS-FUSION approaches with Mendelian Randomization and colocalization to identify causal metabolites for 27 brain and human wellness phenotypes and identified eight metabolites to be causal for eight traits (11 relationships). Low mannose level was causal to bipolar disorder and as dietary supplement it may provide therapeutic benefits. Low galactosylglycerol level was found causal to Parkinson's Disease (PD). Our study expanded the knowledge of MQTL in central nervous system, provided insights into human wellness, and successfully demonstrates the utility of combined statistical approaches to inform interventions.

Genetic studies of paired metabolomes reveal enzymatic and transport processes at the interface of plasma and urine

The kidneys operate at the interface of plasma and urine by clearing molecular waste products while retaining valuable solutes. Genetic studies of paired plasma and urine metabolomes may identify underlying processes. We conducted genome-wide studies of 1,916 plasma and urine metabolites and detected 1,299 significant associations. Associations with 40% of implicated metabolites would have been missed by studying plasma alone. We detected urine-specific findings that provide information about metabolite reabsorption in the kidney, such as aquaporin (AQP)-7-mediated glycerol transport, and different metabolomic footprints of kidney-expressed proteins in plasma and urine that are consistent with their localization and function, including the transporters NaDC3 (SLC13A3) and ASBT (SLC10A2). Shared genetic determinants of 7,073 metabolite-disease combinations represent a resource to better understand metabolic diseases and revealed connections of dipeptidase 1 with circulating digestive enzymes and with hypertension. Extending genetic studies of the metabolome beyond plasma yields unique insights into processes at the interface of body compartments.

Non-coding and intergenic genetic variants of human arylamine N-acetyltransferase 2 (NAT2) gene are associated with differential plasma lipid and cholesterol levels and cardiometabolic disorders

Arylamine N-acetyltransferase 2 (NAT2) is a phase II metabolic enzyme, best known for metabolism of aromatic amines and hydrazines. Genetic variants occurring in the NAT2 coding region have been well-defined and are known to affect the enzyme activity or protein stability. Individuals can be categorized into rapid, intermediate, and slow acetylator phenotypes that significantly alter their ability to metabolize arylamines, including drugs (e.g., isoniazid) and carcinogens (e.g., 4-aminobiphenyl). However, functional studies on non-coding or intergenic variants of NAT2 are lacking. Multiple, independent genome wide association studies (GWAS) have reported that non-coding or intergenic variants of NAT2 are associated with elevated plasma lipid and cholesterol levels, as well as cardiometabolic disorders, suggesting a novel cellular role of NAT2 in lipid and cholesterol homeostasis. The current review highlights and summarizes GWAS reports that are relevant to this association. We also present a new finding that seven, non-coding, intergenic NAT2 variants (i.e., rs4921913, rs4921914, rs4921915, rs146812806, rs35246381, rs35570672, and rs1495741), which have been associated with plasma lipid and cholesterol levels, are in linkage disequilibrium with one another, and thus form a novel haplotype. The dyslipidemia risk alleles of non-coding NAT2 variants are associated with rapid NAT2 acetylator phenotype, suggesting that differential systemic NAT2 activity might be a risk factor for developing dyslipidemia. The current review also discusses the findings of recent reports that are supportive of the role of NAT2 in lipid or cholesterol synthesis and transport. In summary, we review data suggesting that human NAT2 is a novel genetic factor that influences plasma lipid and cholesterol levels and alters the risk of cardiometabolic disorders. The proposed novel role of NAT2 merits further investigations.

Precision sirolimus dosing in children: The potential for model-informed dosing and novel drug monitoring

The mTOR inhibitor sirolimus is prescribed to treat children with varying diseases, ranging from vascular anomalies to sporadic lymphangioleiomyomatosis to transplantation (solid organ or hematopoietic cell). Precision dosing of sirolimus using therapeutic drug monitoring (TDM) of sirolimus concentrations in whole blood drawn at the trough (before the next dose) time-point is the current standard of care. For sirolimus, trough concentrations are only modestly correlated with the area under the curve, with R 2 values ranging from 0.52 to 0.84. Thus, it should not be surprising, even with the use of sirolimus TDM, that patients treated with sirolimus have variable pharmacokinetics, toxicity, and effectiveness. Model-informed precision dosing (MIPD) will be beneficial and should be implemented. The data do not suggest dried blood spots point-of-care sampling of sirolimus concentrations for precision dosing of sirolimus. Future research on precision dosing of sirolimus should focus on pharmacogenomic and pharmacometabolomic tools to predict sirolimus pharmacokinetics and wearables for point-of-care quantitation and MIPD.

Cross-Ancestry Genome-Wide Association Study Defines the Extended CYP2D6 Locus as the Principal Genetic Determinant of Endoxifen Plasma Concentrations

The therapeutic efficacy of tamoxifen is predominantly mediated by its active metabolites 4-hydroxy-tamoxifen and endoxifen, whose formation is catalyzed by the polymorphic cytochrome P450 2D6 (CYP2D6). Yet, known CYP2D6 polymorphisms only partially determine metabolite concentrations in vivo. We performed the first cross-ancestry genome-wide association study with well-characterized patients of European, Middle-Eastern, and Asian descent (n = 497) to identify genetic factors impacting active and parent metabolite formation. Genome-wide significant variants were functionally evaluated in an independent liver cohort (n = 149) and in silico. Metabolite prediction models were validated in two independent European breast cancer cohorts (n = 287, n = 189). Within a single 1-megabase (Mb) region of chromosome 22q13 encompassing the CYP2D6 gene, 589 variants were significantly associated with tamoxifen metabolite concentrations, particularly endoxifen and metabolic ratio (MR) endoxifen/N-desmethyltamoxifen (minimal P = 5.4E-35 and 2.5E-65, respectively). Previously suggested other loci were not confirmed. Functional analyses revealed 66% of associated, mostly intergenic variants to be significantly correlated with hepatic CYP2D6 activity or expression (ρ = 0.35 to -0.52), and six hotspot regions in the extended 22q13 locus impacting gene regulatory function. Machine learning models based on hotspot variants (n = 12) plus CYP2D6 activity score (AS) increased the explained variability (~ 9%) compared with AS alone, explaining up to 49% (median R² ) and 72% of the variability in endoxifen and MR endoxifen/N-desmethyltamoxifen, respectively. Our findings suggest that the extended CYP2D6 locus at 22q13 is the principal genetic determinant of endoxifen plasma concentration. Long-distance haplotypes connecting CYP2D6 with adjacent regulatory sites and nongenetic factors may account for the unexplained portion of variability.

Metabolic profiling of the fluorinated liquid-crystal monomer 1-ethoxy-2,3-difluoro-4-(trans-4-propylcyclohexyl)benzene

1-Ethoxy-2,3-difluoro-4-(trans-4-propylcyclohexyl)benzene (EDPrB) is a typical fluorinated liquid-crystal monomer (LCM). LCMs contaminants are becoming increasingly concerning due to their potential persistence, bioaccumulation, toxicity, and broad prevalence in environmental and human samples. However, LCM metabolism is poorly understood. Herein, by introducing selected EDPrB into the appropriate liver microsomes in vitro, we examined the metabolic pathways of LCM in humans, rats, pigs, Cyprinus carpio, crucian carp, and Channa argus. A total of 20 species-dependent metabolites were identified and structurally elucidated by gas and liquid chromatography-high resolution mass spectrometry for the first time. Dealkylation, H-abstraction, and hydroxylation reactions are the primary metabolic pathways. Half of these in vitro metabolites were found in the urine, serum, and fecal samples of Sprague-Dawley rats exposed to EDPrB. Toxicity predictions indicate that 17 metabolites can be classified as toxic. According to the Ecological Structure Activity Relationships (ECOSAR), a number of metabolites exhibit equivalent or greater aquatic toxicity to that of EDPrB. Toxicity Estimation Software Tool (T.E.S.T.) predicts that some metabolites exhibit developmental toxicity and mutagenicity in rats. These findings suggest that biotransformation should be particularly emphasized, and more toxicological and monitoring studies should be performed to assess the ecological and human safety of LCMs.

Molecular pathways that drive diabetic kidney disease

Kidney disease is a major driver of mortality among patients with diabetes and diabetic kidney disease (DKD) is responsible for close to half of all chronic kidney disease cases. DKD usually develops in a genetically susceptible individual as a result of poor metabolic (glycemic) control. Molecular and genetic studies indicate the key role of podocytes and endothelial cells in driving albuminuria and early kidney disease in diabetes. Proximal tubule changes show a strong association with the glomerular filtration rate. Hyperglycemia represents a key cellular stress in the kidney by altering cellular metabolism in endothelial cells and podocytes and by imposing an excess workload requiring energy and oxygen for proximal tubule cells. Changes in metabolism induce early adaptive cellular hypertrophy and reorganization of the actin cytoskeleton. Later, mitochondrial defects contribute to increased oxidative stress and activation of inflammatory pathways, causing progressive kidney function decline and fibrosis. Blockade of the renin-angiotensin system or the sodium-glucose cotransporter is associated with cellular protection and slowing kidney function decline. Newly identified molecular pathways could provide the basis for the development of much-needed novel therapeutics.

SLC26A1 is a major determinant of sulfate homeostasis in humans

Sulfate plays a pivotal role in numerous physiological processes in the human body, including bone and cartilage health. A role of the anion transporter SLC26A1 (Sat1) for sulfate reabsorption in the kidney is supported by the observation of hyposulfatemia and hypersulfaturia in Slc26a1-knockout mice. The impact of SLC26A1 on sulfate homeostasis in humans remains to be defined. By combining clinical genetics, functional expression assays, and population exome analysis, we identify SLC26A1 as a sulfate transporter in humans and experimentally validate several loss-of-function alleles. Whole-exome sequencing from a patient presenting with painful perichondritis, hyposulfatemia, and renal sulfate wasting revealed a homozygous mutation in SLC26A1, which has not been previously described to the best of our knowledge. Whole-exome data analysis of more than 5,000 individuals confirmed that rare, putatively damaging SCL26A1 variants were significantly associated with lower plasma sulfate at the population level. Functional expression assays confirmed a substantial reduction in sulfate transport for the SLC26A1 mutation of our patient, which we consider to be novel, as well as for the additional variants detected in the population study. In conclusion, combined evidence from 3 complementary approaches supports SLC26A1 activity as a major determinant of sulfate homeostasis in humans. In view of recent evidence linking sulfate homeostasis with back pain and intervertebral disc disorder, our study identifies SLC26A1 as a potential target for modulation of musculoskeletal health.

Calcium-channel blockers: Clinical outcome associations with reported pharmacogenetics variants in 32 000 patients

AIMS

Pharmacogenetic variants impact dihydropyridine calcium-channel blockers (dCCBs; e.g., amlodipine) treatment efficacy, yet evidence on clinical outcomes in routine primary care is limited. Reported associations in pharmacogenomics knowledge base PharmGKB have weak supporting evidence. We aimed to estimate associations between reported pharmacogenetic variants and incident adverse events in a community-based cohort prescribed dCCB.

METHODS

We analysed up to 32 360 UK Biobank participants prescribed dCCB in primary care (from UK general practices, 1990-2017). We investigated 23 genetic variants. Outcomes were incident diagnosis of coronary heart disease, heart failure (HF), chronic kidney disease, oedema and switching antihypertensive medication.

RESULTS

Participants were aged 40-79 years at first dCCB prescription. Carriers of rs877087 T allele in RYR3 had increased risk of hazard ratio (HF 1.13: 95% confidence interval 1.02 to 1.25, P = .02). Although nonsignificant after multiple testing correction, the association is consistent with prior evidence. We estimated that if rs877087 T allele could experience the same treatment effect as noncarriers, the incidence of HF in patients prescribed dCCB would reduce by 9.2% (95% confidence interval 3.1 to 15.4). In patients with a history of heart disease prior to dCCB (n = 2296), rs877087 homozygotes had increased risk of new coronary heart disease or HF compared to CC variant. rs10898815 in NUMA1 and rs776746 in CYP3A5 increased likelihood of switching to an alternative antihypertensive. The remaining variants were not strongly or consistently associated with studied outcomes.

CONCLUSION

Patients with common genetic variants in NUMA1, CYP3A5 and RYR3 had increased adverse clinical outcomes. Work is needed to establish whether outcomes of dCCB prescribing could be improved by prior knowledge of pharmacogenetics variants supported by clinical evidence of association with adverse events.

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.

Identification of Novel Genetic Risk Factors for Focal Segmental Glomerulosclerosis in Children: Results From the Chronic Kidney Disease in Children (CKiD) Cohort

RATIONALE & OBJECTIVE

Focal segmental glomerulosclerosis (FSGS) is a major cause of pediatric nephrotic syndrome, and African Americans exhibit an increased risk for developing FSGS compared with other populations. Predisposing genetic factors have previously been described in adults. Here we performed genomic screening of primary FSGS in a pediatric African American population.

STUDY DESIGN

Prospective cohort with case-control genetic association study design.

SETTING & PARTICIPANTS

140 African American children with chronic kidney disease from the Chronic Kidney Disease in Children (CKiD) cohort, including 32 cases with FSGS.

PREDICTORS

Over 680,000 common single-nucleotide polymorphisms (SNPs) were tested for association. We also ran a pathway enrichment analysis and a human leucocyte antigen (HLA)-focused association study.

OUTCOME

Primary biopsy-proven pediatric FSGS.

ANALYTICAL APPROACH

Multivariate logistic regression models.

RESULTS

The genome-wide association study revealed 169 SNPs from 14 independent loci significantly associated with FSGS (false discovery rate [FDR]<5%). We observed notable signals for genetic variants within the APOL1 (P=8.6×10^-7; OR, 25.8 [95% CI, 7.1-94.0]), ALMS1 (P=1.3×10^-7; 13.0% in FSGS cases vs 0% in controls), and FGFR4 (P=4.3×10^-6; OR, 24.8 [95% CI, 6.3-97.7]) genes, all of which had previously been associated with adult FSGS, kidney function, or chronic kidney disease. We also highlighted novel, functionally relevant genes, including GRB2 (which encodes a slit diaphragm protein promoting podocyte structure through actin polymerization) and ITGB1 (which is linked to renal injuries). Our results suggest a major role for immune responses and antigen presentation in pediatric FSGS through (1) associations with SNPs in PTPRJ (or CD148, P=3.5×10^-7), which plays a role in T-cell receptor signaling, (2) HLA-DRB1∗11:01 association (P=6.1×10^-3; OR, 4.5 [95% CI, 1.5-13.0]), and (3) signaling pathway enrichment (P=1.3×10^-6).

LIMITATIONS

Sample size and no independent replication cohort with genomic data readily available.

CONCLUSIONS

Our genetic study has identified functionally relevant risk factors and the importance of immune regulation for pediatric primary FSGS, which contributes to a better description of its molecular pathophysiological mechanisms.

NetSHy: network summarization via a hybrid approach leveraging topological properties

MOTIVATION

Biological networks can provide a system-level understanding of underlying processes. In many contexts, networks have a high degree of modularity, i.e. they consist of subsets of nodes, often known as subnetworks or modules, which are highly interconnected and may perform separate functions. In order to perform subsequent analyses to investigate the association between the identified module and a variable of interest, a module summarization, that best explains the module's information and reduces dimensionality is often needed. Conventional approaches for obtaining network representation typically rely only on the profiles of the nodes within the network while disregarding the inherent network topological information.

RESULTS

In this article, we propose NetSHy, a hybrid approach which is capable of reducing the dimension of a network while incorporating topological properties to aid the interpretation of the downstream analyses. In particular, NetSHy applies principal component analysis (PCA) on a combination of the node profiles and the well-known Laplacian matrix derived directly from the network similarity matrix to extract a summarization at a subject level. Simulation scenarios based on random and empirical networks at varying network sizes and sparsity levels show that NetSHy outperforms the conventional PCA approach applied directly on node profiles, in terms of recovering the true correlation with a phenotype of interest and maintaining a higher amount of explained variation in the data when networks are relatively sparse. The robustness of NetSHy is also demonstrated by a more consistent correlation with the observed phenotype as the sample size decreases. Lastly, a genome-wide association study is performed as an application of a downstream analysis, where NetSHy summarization scores on the biological networks identify more significant single nucleotide polymorphisms than the conventional network representation.

AVAILABILITY AND IMPLEMENTATION

R code implementation of NetSHy is available at https://github.com/thaovu1/NetSHy.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Combing Genome-Wide Association Studies and Single-Cell Analysis to Elucidate the Mechanisms of Kidney Disease: Proceedings of the Henry Shavelle Professorship

Background

Kidney diseases pose a significant global health burden; there is an urgent need to deepen our understanding of their underlying mechanisms.

Summary

This review focuses on new innovative approaches that merge genome-wide association studies (GWAS) and single-cell omics (including transcriptomics) in kidney disease research. We begin by detailing how GWAS has identified numerous genetic risk factors, offering valuable insight into disease susceptibility. Then, we explore the application of scRNA-seq, highlighting its ability to unravel how genetic variants influence cellular phenotypes. Through a synthesis of recent studies, we illuminate the synergy between these two powerful methodologies, demonstrating their potential in elucidating the complex etiology of kidney diseases. Moreover, we discuss how this integrative approach could pave the way for precise diagnostics and personalized treatments.

Key Message

This review underscores the transformative potential of combining GWAS and scRNA-seq in the journey toward a deeper understanding of kidney diseases.

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.

Inherited genetic effects on arsenic metabolism: A comparison of effects on arsenic species measured in urine and in blood

Inorganic arsenic (iAs) is a carcinogen, and chronic exposure is associated with adverse health outcomes, including cancer and cardiovascular disease. Consumed iAs can undergo two methylation reactions catalyzed by arsenic methyltransferase (AS3MT), producing monomethylated and dimethylated forms of arsenic (MMA and DMA). Methylation of iAs helps facilitate excretion of arsenic in urine, with DMA composing the majority of arsenic species excreted. Past studies have identified genetic variation in the AS3MT (10q24.32) and FTCD (21q22.3) regions associated with arsenic metabolism efficiency (AME), measured as the proportion of each species present in urine (iAs%, MMA%, and DMA%), but their association with arsenic species present in blood has not been examined. We use data from three studies nested within the Health Effects and Longitudinal Study (HEALS)-the Nutritional Influences on Arsenic Toxicity Study, the Folate and Oxidative Stress study, and the Folic Acid and Creatine Trial-to examine the association of previously identified genetic variants with arsenic species in both urine and blood of 334 individuals. We confirm that the genetic variants in AS3MT and FTCD known to effect arsenic species composition in urine (an excreted byproduct of metabolism) have similar effects on arsenic species in blood (a tissue type that directly interacts with many organs, including those prone to arsenic toxicity). This consistency we observe provides further support for the hypothesis the AME SNPs identified to date impact the efficiency of arsenic metabolism and elimination, thereby influencing internal dose of arsenic and the dose delivered to toxicity-prone organs and tissues.

Sex Hormone Candidate Gene Polymorphisms Are Associated with Endometriosis

The present study was designed to examine whether sex hormone polymorphisms proven by GWAS are associated with endometriosis risk. Unrelated female participants totaling 1376 in number (395 endometriosis patients and 981 controls) were recruited into the study. Nine single-nucleotide polymorphisms (SNPs) which GWAS correlated with circulating levels of sex hormones were genotyped using a TaqMan allelic discrimination assay. FSH-lowering, and LH- and testosterone-heightening polymorphisms of the FSHB promoter (allelic variants A rs11031002 and C rs11031005) exhibit a protective effect for endometriosis (OR = 0.60-0.68). By contrast, the TT haplotype loci that were GWAS correlated with higher FSH levels and lower LH and testosterone concentrations determined an increased risk for endometriosis (OR = 2.03). Endometriosis-involved epistatic interactions were found between eight loci of sex hormone genes (without rs148982377 ZNF789) within twelve genetic simulation models. In silico examination established that 8 disorder-related loci and 80 proxy SNPs are genome variants affecting the expression, splicing, epigenetic and amino acid conformation of the 34 genes which enrich the organic anion transport and secondary carrier transporter pathways. In conclusion, the present study showed that sex hormone polymorphisms proven by GWAS are associated with endometriosis risk and involved in the molecular pathophysiology of the disease due to their functionality.

Integrated proteomic and metabolomic modules identified as biomarkers of mortality in the Atherosclerosis Risk in Communities study and the African American Study of Kidney Disease and Hypertension

BACKGROUND

Proteins and metabolites are essential for many biological functions and often linked through enzymatic or transport reactions. Individual molecules have been associated with all-cause mortality. Many of these are correlated and might jointly represent pathways or endophenotypes involved in diseases.

RESULTS

We present an integrated analysis of proteomics and metabolomics via a local dimensionality reduction clustering method. We identified 224 modules of correlated proteins and metabolites in the Atherosclerosis Risk in Communities (ARIC) study, a general population cohort of older adults (N = 4046, mean age 75.7, mean eGFR 65). Many of the modules displayed strong cross-sectional associations with demographic and clinical characteristics. In comprehensively adjusted analyses, including fasting plasma glucose, history of cardiovascular disease, systolic blood pressure and kidney function among others, 60 modules were associated with mortality. We transferred the network structure to the African American Study of Kidney Disease and Hypertension (AASK) (N = 694, mean age 54.5, mean mGFR 46) and identified mortality associated modules relevant in this disease specific cohort. The four mortality modules relevant in both the general population and CKD were all a combination of proteins and metabolites and were related to diabetes / insulin secretion, cardiovascular disease and kidney function. Key components of these modules included N-terminal (NT)-pro hormone BNP (NT-proBNP), Sushi, Von Willebrand Factor Type A, EGF And Pentraxin (SVEP1), and several kallikrein proteases.

CONCLUSION

Through integrated biomarkers of the proteome and metabolome we identified functions of (patho-) physiologic importance related to diabetes, cardiovascular disease and kidney function.

Influence of the microbiome, diet and genetics on inter-individual variation in the human plasma metabolome

The levels of the thousands of metabolites in the human plasma metabolome are strongly influenced by an individual's genetics and the composition of their diet and gut microbiome. Here, by assessing 1,183 plasma metabolites in 1,368 extensively phenotyped individuals from the Lifelines DEEP and Genome of the Netherlands cohorts, we quantified the proportion of inter-individual variation in the plasma metabolome explained by different factors, characterizing 610, 85 and 38 metabolites as dominantly associated with diet, the gut microbiome and genetics, respectively. Moreover, a diet quality score derived from metabolite levels was significantly associated with diet quality, as assessed by a detailed food frequency questionnaire. Through Mendelian randomization and mediation analyses, we revealed putative causal relationships between diet, the gut microbiome and metabolites. For example, Mendelian randomization analyses support a potential causal effect of Eubacterium rectale in decreasing plasma levels of hydrogen sulfite-a toxin that affects cardiovascular function. Lastly, based on analysis of the plasma metabolome of 311 individuals at two time points separated by 4 years, we observed a positive correlation between the stability of metabolite levels and the amount of variance in the levels of that metabolite that could be explained in our analysis. Altogether, characterization of factors that explain inter-individual variation in the plasma metabolome can help design approaches for modulating diet or the gut microbiome to shape a healthy metabolome.

Identification of 969 protein quantitative trait loci in an African American population with kidney disease attributed to hypertension

Investigations into the causal underpinnings of disease processes can be aided by the incorporation of genetic information. Genetic studies require populations varied in both ancestry and prevalent disease in order to optimize discovery and ensure generalizability of findings to the global population. Here, we report the genetic determinants of the serum proteome in 466 African Americans with chronic kidney disease attributed to hypertension from the richly phenotyped African American Study of Kidney Disease and Hypertension (AASK) study. Using the largest aptamer-based protein profiling platform to date (6,790 proteins or protein complexes), we identified 969 genetic associations with 900 unique proteins; including 52 novel cis (local) associations and 379 novel trans (distant) associations. The genetic effects of previously published cis-protein quantitative trait loci (pQTLs) were found to be highly reproducible, and we found evidence that our novel genetic signals colocalize with gene expression and disease processes. Many trans- pQTLs were found to reflect associations mediated by the circulating cis protein, and the common trans-pQTLs are enriched for processes involving extracellular vesicles, highlighting a plausible mechanism for distal regulation of the levels of secreted proteins. Thus, our study generates a valuable resource of genetic associations linking variants to protein levels and disease in an understudied patient population to inform future studies of drug targets and physiology.