Cerebrospinal fluid metabolomics identifies 19 brain-related phenotype associations

Brain development and bioenergetic changes

Bioenergetics describe the biochemical processes responsible for energy supply in organisms. When these changes become dysregulated in brain development, multiple neurodevelopmental diseases can occur, implicating bioenergetics as key regulators of neural development. Historically, the discovery of disease processes affecting individual stages of brain development has revealed critical roles that bioenergetics play in generating the nervous system. Bioenergetic-dependent neurodevelopmental disorders include neural tube closure defects, microcephaly, intellectual disability, autism spectrum disorders, epilepsy, mTORopathies, and oncogenic processes. Developmental timing and cell-type specificity of these changes determine the long-term effects of bioenergetic disease mechanisms on brain form and function. Here, we discuss key metabolic regulators of neural progenitor specification, neuronal differentiation (neurogenesis), and gliogenesis. In general, transitions between glycolysis and oxidative phosphorylation are regulated in early brain development and in oncogenesis, and reactive oxygen species (ROS) and mitochondrial maturity play key roles later in differentiation. We also discuss how bioenergetics interface with the developmental regulation of other key neural elements, including the cerebrospinal fluid brain environment. While questions remain about the interplay between bioenergetics and brain development, this review integrates the current state of known key intersections between these processes in health and disease.

A metabolome-wide Mendelian randomization study prioritizes causal circulating metabolites for reproductive disorders including primary ovarian insufficiency, polycystic ovary syndrome, and abnormal spermatozoa

BACKGROUND

Accumulating studies have highlighted the significant role of circulating metabolomics in the etiology of reproductive system disorders. However, the causal effects between genetically determined metabolites (GDMs) and reproductive diseases, including primary ovarian insufficiency (POI), polycystic ovary syndrome (PCOS), and abnormal spermatozoa (AS), still await thorough clarification.

METHODS

With the currently most comprehensive genome-wide association studies (GWAS) data of metabolomics, systematic two-sample Mendelian randomization (MR) analyses were conducted to disclose causal associations between 1,091 blood metabolites and 309 metabolite ratios with reproductive disorders. The inverse-variance weighted (IVW) method served as the primary analysis approach, and multiple effective MR methods were employed as complementary analyses including MR-Egger, weighted median, constrained maximum likelihood (cML-MA), contamination mixture method, robust adjusted profile score (MR-RAPS), and debiased inverse-variance weighted method. Heterogeneity and pleiotropy were assessed via MR-Egger intercept and Cochran's Q statistical analysis. Outliers were detected by Radial MR and MR-PRESSO methods. External replication and metabolic pathway analysis were also conducted.

RESULTS

Potential causal associations of 63 GDMs with POI were unearthed, and five metabolites with strong causal links to POI were emphasized. Two metabolic pathways related to the pathogenesis of POI were pinpointed. Suggestive causal effects of 70 GDMs on PCOS were detected, among which 7 metabolites stood out for strong causality with elevated PCOS risk. Four metabolic pathways associated with PCOS mechanisms were recognized. For AS, 64 GDMs as potential predictive biomarkers were identified, particularly highlighting two metabolites for their strong causal connections with AS. Three pathways underneath the AS mechanism were identified. Multiple assessments were conducted to further confirm the reliability and robustness of our causal inferences.

CONCLUSION

By extensively assessing the causal implications of circulating GDMs on reproductive system disorders, our study underscores the intricate and pivotal role of metabolomics in reproductive ill-health, laying a theoretical foundation for clinical strategies from metabolic insights.

Exploring the Bidirectional Effects of Gut Microbiota and Short-Chain Fatty Acids on Urticaria Subtypes Through Mendelian Randomization and Mediation Analysis

Background

Emerging evidence links gut microbiota and their by-products, notably short-chain fatty acids (SCFAs), to urticaria. This study employs multiple Mendelian Randomization (MR) analyses to unravel the complex interactions among gut microbiota, SCFAs, and different subtypes of urticaria, aiming to elucidate the underlying mechanisms and enhance future clinical research.

Methods

We analyzed published genome-wide association study (GWAS) summary statistics to identify associations between gut microbiota and three common subtypes of urticaria: spontaneous, dermatographic, and temperature-triggered. Initial two-sample and reverse MR analyses explored the causality in these relationships. Subsequent multivariate MR analyses investigated the role of SCFAs in modulating these interactions, with multiple sensitivity analyses to ensure robustness.

Findings

Specific taxa were differently associated with various urticaria subtypes. From microbiota to urticaria: one taxon was negatively associated with dermatographic urticaria; seven taxa were negatively associated and four positively associated with temperature-triggered urticaria; four taxa were negatively associated and six positively associated with spontaneous urticaria. Conversely, from urticaria to microbiota: five taxa were negatively associated with dermatographic urticaria; four were negatively and two positively associated with temperature-triggered urticaria; and two were negatively associated with spontaneous urticaria. These associations were observed at a nominal significance level (P < 0.05). After applying Bonferroni correction for multiple testing, these associations did not reach statistical significance. The observed trends, however, provide insights into potential microbiota-urticaria interactions. Multivariate MR analyses elucidated the role of SCFAs, particularly acetate, which plays a crucial role in modulating immune response. Adjusting for acetate revealed direct effects of Actinobacteria, Bifidobacteriales, and Bifidobacteriaceae on spontaneous urticaria, with corresponding mediation effects of -22%, -24.9%, and -24.9% respectively. Similarly, adjustments for Alcaligenaceae and Betaproteobacteria indicated significant negative effects of acetate on dermatographic and spontaneous urticaria, with mediation effects of -21.7% and -23.7%, respectively.

Conclusion

This study confirms the interconnected roles of gut microbiota, SCFAs, and urticaria. It highlights SCFAs' potential mediating role in influencing urticaria through microbiota, providing insights for future therapeutic strategies.

Pharmacological targets of SGLT2 inhibition on prostate cancer mediated by circulating metabolites: a drug-target Mendelian randomization study

Background

The relationship between sodium-glucose cotransporter 2 (SGLT2) inhibitors and prostate cancer is still unknown. Although these inhibitors can influence tumor glycolysis, the underlying mechanism requires further exploration.

Methods

A two-sample two-step MR was used to determine 1) causal effects of SGLT2 inhibition on prostate cancer; 2) causal effects of 1,400 circulating metabolites or metabolite ratios on prostate cancer; and 3) mediation effects of these circulating metabolites. Genetic proxies for SGLT2 inhibition were identified as variants in the SLC5A2 gene and glycated hemoglobin level (HbA1c). Additionally, positive control analysis on type 2 diabetes mellitus (T2DM) was conducted to test the selection of genetic proxies. Phenome Wide Association Study (PheWAS) and MR-PheWAS analysis were used to explore potential treatable diseases and adverse outcomes of SGLT2 inhibitors.

Results

Genetically predicted SGLT2 inhibition (per 1 SD decrement in HbA1c) was associated with reduced risk of T2DM [odds ratio (OR) = 0.66 (95% CI 0.53, 0.82), P = 1.57 × 10-4]; prostate cancer [0.34 (0.23, 0.49), P = 2.21 × 10-8] and prostate-specific antigen [0.26 (0.08, 0.81), P = 2.07 × 10-2]. The effect of SGLT2 inhibition on prostate cancer was mediated by uridine level, with a mediated proportion of 9.34% of the total effect. In MR-PheWAS, 65 traits were found to be associated with SLGT2 inhibitors (P < 1.78 × 10-5), and among them, 13 were related to diabetes.

Conclusion

Our study suggested that SGLT2 inhibition could lower prostate cancer risk through uridine mediation. More mechanistic and clinical research is necessary to explore how uridine mediates the link between SGLT2 inhibition and prostate cancer.

Genome-wide association study of exercise-induced skeletal muscle hypertrophy and the construction of predictive model

The aim of the current study was to investigate interindividual differences in muscle thickness of the rectus femoris (MTRF) following 12 wk of resistance training (RT) or high-intensity interval training (HIIT) to explore the genetic architecture underlying skeletal muscle hypertrophy and to construct predictive models. We conducted musculoskeletal ultrasound assessments of the MTRF response in 440 physically inactive adults after the 12-wk exercise period. A genome-wide association study was used to identify variants associated with the MTRF response, separately for RT and HIIT. Using the polygenic predictor score (PPS), we estimated the genetic contribution to exercise-induced hypertrophy. Predictive models for the MTRF response were constructed using random forest (RF), support vector mac (SVM), and generalized linear model (GLM) in 10 cross-validated approaches. MTRF increased significantly after both RT (8.8%, P < 0.05) and HIIT (5.3%, P < 0.05), but with considerable interindividual differences (RT: -13.5 to 38.4%, HIIT: -14.2 to 30.7%). Eleven lead single-nucleotide polymorphisms in RT and eight lead single-nucleotide polymorphisms in HIIT were identified at a significance level of P < 1 × 10-5. The PPS was associated with the MTRF response, explaining 47.2% of the variation in response to RT and 38.3% of the variation in response to HIIT. Notably, the GLM and SVM predictive models exhibited superior performance compared with RF models (P < 0.05), and the GLM demonstrated optimal performance with an area under curve of 0.809 (95% confidence interval: 0.669-0.949). Factors such as PPS, baseline MTRF, and exercise protocol exerted influence on the MTRF response to exercise, with PPS being the primary contributor. The GLM and SVM predictive model, incorporating both genetic and phenotypic factors, emerged as promising tools for predicting exercise-induced skeletal muscle hypertrophy.NEW & NOTEWORTHY The interindividual variability induced muscle hypertrophy by resistance training (RT) or high-intensity interval training (HIIT) and the associated genetic architecture remain uncertain. We identified genetic variants that underlie RT- or HIIT-induced muscle hypertrophy and established them as pivotal factors influencing the response regardless of the training type. The genetic-phenotype predictive model developed has the potential to identify nonresponders or individuals with low responsiveness before engaging in exercise training.

Hierarchical joint analysis of marginal summary statistics-Part II: High-dimensional instrumental analysis of omics data

Instrumental variable (IV) analysis has been widely applied in epidemiology to infer causal relationships using observational data. Genetic variants can also be viewed as valid IVs in Mendelian randomization and transcriptome-wide association studies. However, most multivariate IV approaches cannot scale to high-throughput experimental data. Here, we leverage the flexibility of our previous work, a hierarchical model that jointly analyzes marginal summary statistics (hJAM), to a scalable framework (SHA-JAM) that can be applied to a large number of intermediates and a large number of correlated genetic variants-situations often encountered in modern experiments leveraging omic technologies. SHA-JAM aims to estimate the conditional effect for high-dimensional risk factors on an outcome by incorporating estimates from association analyses of single-nucleotide polymorphism (SNP)-intermediate or SNP-gene expression as prior information in a hierarchical model. Results from extensive simulation studies demonstrate that SHA-JAM yields a higher area under the receiver operating characteristics curve (AUC), a lower mean-squared error of the estimates, and a much faster computation speed, compared to an existing approach for similar analyses. In two applied examples for prostate cancer, we investigated metabolite and transcriptome associations, respectively, using summary statistics from a GWAS for prostate cancer with more than 140,000 men and high dimensional publicly available summary data for metabolites and transcriptomes.

Metabolome-Wide Mendelian Randomization Assessing the Causal Role of Serum and Cerebrospinal Metabolites in Traumatic Brain Injury

Previous studies have identified metabolites as biomarkers or potential therapeutic targets for traumatic brain injury (TBI). However, the causal association between them remains unknown. Therefore, we investigated the causal effect of serum metabolites and cerebrospinal fluid (CSF) metabolites on TBI susceptibility through Mendelian randomization (MR). Genetic variants related to metabolites and TBI were extracted from a corresponding genome-wide association study (GWAS). Causal effects were estimated through the inverse variance weighted approach, supplemented by a weighted median, weight mode, and the MR-Egger test. In addition, sensitivity analyses were further performed to evaluate the stability of the MR results, including the MR-Egger intercept, leave-one-out analysis, Cochrane's Q-test, and the MR-PRESSO global test. Metabolic pathway analysis was applied to uncover the underlying pathways of the significant metabolites in TBI. In blood metabolites, substances such as 4-acetaminophen sulfate and kynurenine showed positive links, whereas beta-hydroxyisovalerate and creatinine exhibited negative correlations. CSF metabolites such as N-formylanthranilic acid were positively related, while kynurenate showed negative associations. The metabolic pathway analysis highlighted the potential biological pathways involved in TBI. Of these 16 serum metabolites, 11 CSF metabolites and metabolic pathways may serve as useful circulating biomarkers in clinical screening and prevention, and may be candidate molecules for the exploration of mechanisms and drug targets.

Deciphering the causal relationship between plasma and cerebrospinal fluid metabolites and glioblastoma multiforme: a Mendelian Randomization study

BACKGROUND

Glioblastoma Multiforme (GBM) is one of the most aggressive and fatal brain cancers. The study of metabolites could be crucial for understanding GBM's biology and reveal new treatment strategies.

METHODS

The GWAS data for GBM were sourced from the FinnGen database. A total of 1400 plasma metabolites were collected from the GWAS Catalog dataset. The cerebrospinal fluid (CSF) metabolites data were collected from subsets of participants in the WADRC and WRAP studies. We utilized the inverse variance weighting (IVW) method as the primary tool to explore the causal relationship between metabolites in plasma and CSF and glioblastoma, ensuring the exclusion of instances with horizontal pleiotropy. Additionally, four supplementary analytical methods were applied to reinforce our findings. Aberrant results were identified and omitted based on the outcomes of the leave-one-out sensitivity analysis. Conclusively, a reverse Mendelian Randomization analysis was also conducted to further substantiate our results.

RESULTS

The study identified 69 plasma metabolites associated with GBM. Of these, 40 metabolites demonstrated a significant positive causal relationship with GBM, while 29 exhibited a significant negative causal association. Notably, Trimethylamine N-oxide (TMAO) levels in plasma, not CSF, were found to be a significant exposure factor for GBM (OR = 3.1627, 95% CI = (1.6347, 6.1189), P = 0.0006). The study did not find a reverse causal relationship between GBM and plasma TMAO levels.

CONCLUSIONS

This research has identified 69 plasma metabolites potentially associated with the incidence of GBM, among which TMAO stands out as a promising candidate for an early detectable biomarker for GBM.

Circulating Metabolite Abundances Associated With Risks of Bipolar Disorder, Schizophrenia, and Depression: A Mendelian Randomization Study

BACKGROUND

Preventive measures and treatments for psychiatric disorders are limited. Circulating metabolites are potential candidates for biomarker and therapeutic target identification, given their measurability and essential roles in biological processes.

METHODS

Leveraging large-scale genome-wide association studies, we conducted Mendelian randomization analyses to assess the associations between circulating metabolite abundances and the risks of bipolar disorder, schizophrenia, and depression. Genetic instruments were selected for 94 metabolites measured in the Canadian Longitudinal Study on Aging cohort (N = 8299). We repeated Mendelian randomization analyses based on the UK Biobank, INTERVAL, and EPIC (European Prospective Investigation into Cancer)-Norfolk studies.

RESULTS

After validating Mendelian randomization assumptions and colocalization evidence, we found that a 1 SD increase in genetically predicted circulating abundances of eicosapentaenoate and docosapentaenoate was associated with odds ratios of 0.72 (95% CI, 0.65-0.79) and 0.63 (95% CI, 0.55-0.72), respectively, for bipolar disorder. Genetically increased Ω-3 unsaturated fatty acids abundance and Ω-3-to-total fatty acids ratio, as well as genetically decreased Ω-6-to-Ω-3 ratio, were negatively associated with the risk of bipolar disorder in the UK Biobank. Genetically increased circulating abundances of 3 N-acetyl-amino acids were associated with an increased risk of schizophrenia with a maximum odds ratio of 1.31 (95% CI, 1.18-1.44) per 1 SD increase. Furthermore, a 1 SD increase in genetically predicted circulating abundance of hypotaurine was associated with an odds ratio of 0.85 (95% CI, 0.78-0.93) for depression.

CONCLUSIONS

The biological mechanisms that underlie Ω-3 unsaturated fatty acids, NAT8-catalyzed N-acetyl-amino acids, and hypotaurine warrant exploration to identify new biomarkers and potential therapeutic targets.

Hippocampal metabolic profile during epileptogenesis in the pilocarpine model of epilepsy

Temporal lobe epilepsy (TLE) is a common form of refractory epilepsy in adulthood. The metabolic profile of epileptogenesis is still poorly investigated. Elucidation of such a metabolic profile using animal models of epilepsy could help identify new metabolites and pathways involved in the mechanisms of epileptogenesis process. In this study, we evaluated the metabolic profile during the epileptogenesis periods. Using a pilocarpine model of epilepsy, we analyzed the global metabolic profile of hippocampal extracts by untargeted metabolomics based on ultra-performance liquid chromatography-high-resolution mass spectrometry, at three time points (3 h, 1 week, and 2 weeks) after status epilepticus (SE) induction. We demonstrated that epileptogenesis periods presented different hippocampal metabolic profiles, including alterations of metabolic pathways of amino acids and lipid metabolism. Six putative metabolites (tryptophan, N-acetylornithine, N-acetyl-L-aspartate, glutamine, adenosine, and cholesterol) showed significant different levels during epileptogenesis compared to their respective controls. These putative metabolites could be associated with the imbalance of neurotransmitters, mitochondrial dysfunction, and cell loss observed during both epileptogenesis and epilepsy. With these findings, we provided an overview of hippocampal metabolic profiles during different stages of epileptogenesis that could help investigate pathways and respective metabolites as predictive tools in epilepsy.

The genomic alterations in glioblastoma influence the levels of CSF metabolites

Cerebrospinal fluid (CSF) analysis is underutilized in patients with glioblastoma (GBM), partly due to a lack of studies demonstrating the clinical utility of CSF biomarkers. While some studies show the utility of CSF cell-free DNA analysis, studies analyzing CSF metabolites in patients with glioblastoma are limited. Diffuse gliomas have altered cellular metabolism. For example, mutations in isocitrate dehydrogenase enzymes (e.g., IDH1 and IDH2) are common in diffuse gliomas and lead to increased levels of D-2-hydroxyglutarate in CSF. However, there is a poor understanding of changes CSF metabolites in GBM patients. In this study, we performed targeted metabolomic analysis of CSF from n = 31 patients with GBM and n = 13 individuals with non-neoplastic conditions (controls), by mass spectrometry. Hierarchical clustering and sparse partial least square-discriminant analysis (sPLS-DA) revealed differences in CSF metabolites between GBM and control CSF, including metabolites associated with fatty acid oxidation and the gut microbiome (i.e., carnitine, 2-methylbutyrylcarnitine, shikimate, aminobutanal, uridine, N-acetylputrescine, and farnesyl diphosphate). In addition, we identified differences in CSF metabolites in GBM patients based on the presence/absence of TP53 or PTEN mutations, consistent with the idea that different mutations have different effects on tumor metabolism. In summary, our results increase the understanding of CSF metabolites in patients with diffuse gliomas and highlight several metabolites that could be informative biomarkers in patients with GBM.

Critical Factors in Sample Collection and Preparation for Clinical Metabolomics of Underexplored Biological Specimens

This review article compiles critical pre-analytical factors for sample collection and extraction of eight uncommon or underexplored biological specimens (human breast milk, ocular fluids, sebum, seminal plasma, sweat, hair, saliva, and cerebrospinal fluid) under the perspective of clinical metabolomics. These samples are interesting for metabolomics studies as they reflect the status of living organisms and can be applied for diagnostic purposes and biomarker discovery. Pre-collection and collection procedures are critical, requiring protocols to be standardized to avoid contamination and bias. Such procedures must consider cleaning the collection area, sample stimulation, diet, and food and drug intake, among other factors that impact the lack of homogeneity of the sample group. Precipitation of proteins and removal of salts and cell debris are the most used sample preparation procedures. This review intends to provide a global view of the practical aspects that most impact results, serving as a starting point for the designing of metabolomic experiments.

Metabolomics of Mouse Embryonic CSF Following Maternal Immune Activation

The cerebrospinal fluid (CSF) serves various roles in the developing central nervous system (CNS), from neurogenesis to lifelong cognitive functions. Changes in CSF composition due to inflammation can impact brain function. We recently identified an abnormal cytokine signature in embryonic CSF (eCSF) following maternal immune activation (MIA), a mouse model of autism spectrum disorder (ASD). We hypothesized that MIA leads to other alterations in eCSF composition and employed untargeted metabolomics to profile changes in the eCSF metabolome in mice after inducing MIA with polyI:C. We report these data here as a resource, include a comprehensive MS1 and MS2 reference dataset, and present additional datasets comparing two mouse strains (CD-1 and C57Bl/6) and two developmental time points (E12.5 and E14.5). Targeted metabolomics further validated changes upon MIA. We show a significant elevation of glucocorticoids and kynurenine pathway related metabolites. Both pathways are relevant for suppressing inflammation or could be informative as disease biomarkers. Our resource should inform future mechanistic studies regarding the etiology of MIA neuropathology and roles and contributions of eCSF metabolites to brain development.

Rapid Metabolic Profiling of 1 μL Crude Cerebrospinal Fluid by Matrix-Assisted Laser Desorption/Ionization Mass Spectrometry Imaging Can Differentiate De Novo Parkinson's Disease

Parkinson's disease (PD) is a highly prevalent neurodegenerative disorder affecting the motor system. However, the correct diagnosis of PD and atypical parkinsonism may be difficult with high clinical uncertainty. There is an urgent need to identify reliable biomarkers using high-throughput, molecular-specific methods to improve current diagnostics. Here, we present a matrix-assisted laser desorption/ionization mass spectrometry imaging method that requires minimal sample preparation and only 1 μL of crude cerebrospinal fluid (CSF). The method enables analysis of hundreds of samples in a single experiment while simultaneously detecting numerous metabolites with subppm mass accuracy. To test the method, we analyzed CSF samples from 12 de novo PD patients (that is, newly diagnosed and previously untreated) and 12 age-matched controls. Within the identified molecules, we found neurotransmitters and their metabolites such as γ-aminobutyric acid, 3-methoxytyramine, homovanillic acid, serotonin, histamine, amino acids, and metabolic intermediates. Limits of detection were estimated for multiple neurotransmitters with high linearity (R2 > 0.99) and sensitivity (as low as 16 pg/μL). Application of multivariate classification led to a highly significant (P < 0.001) model of PD prediction with a 100% classification rate, which was further thoroughly validated with a permutation test and univariate analysis. Molecules related to the neuromelanin pathway were found to be significantly increased in the PD group, indicated by their elevated relative intensities compared to the control group. Our method enables rapid detection of PD-related biomarkers in low sample volumes and could serve as a valuable tool in the development of robust PD diagnostics.

Advances in mass spectrometry-based multi-scale metabolomic methodologies and their applications in biological and clinical investigations

Metabolomics is a nascent field of inquiry that emerged in the late 20th century. It encompasses the comprehensive profiling of metabolites across a spectrum of organisms, ranging from bacteria and cells to tissues. The rapid evolution of analytical methods and data analysis has greatly accelerated progress in this dynamic discipline over recent decades. Sophisticated techniques such as liquid chromatograph mass spectrometry (MS), gas chromatograph MS, capillary electrophoresis MS, and nuclear magnetic resonance serve as the cornerstone of metabolomic analysis. Building upon these methods, a plethora of modifications and combinations have emerged to propel the advancement of metabolomics. Despite this progress, scrutinizing metabolism at the single-cell or single-organelle level remains an arduous task over the decades. Some of the most thrilling advancements, such as single-cell and single-organelle metabolic profiling techniques, offer profound insights into the intricate mechanisms within cells and organelles. This allows for a comprehensive study of metabolic heterogeneity and its pivotal role in multiple biological processes. The progress made in MS imaging has enabled high-resolution in situ metabolic profiling of tissue sections and even individual cells. Spatial reconstruction techniques enable the direct representation of metabolic distribution and alteration in three-dimensional space. The application of novel metabolomic techniques has led to significant breakthroughs in biological and clinical studies, including the discovery of novel metabolic pathways, determination of cell fate in differentiation, anti-aging intervention through modulating metabolism, metabolomics-based clinicopathologic analysis, and surgical decision-making based on on-site intraoperative metabolic analysis. This review presents a comprehensive overview of both conventional and innovative metabolomic techniques, highlighting their applications in groundbreaking biological and clinical studies.

CSF metabolites associated with biomarkers of Alzheimer's disease pathology

Introduction

Metabolomics technology facilitates studying associations between small molecules and disease processes. Correlating metabolites in cerebrospinal fluid (CSF) with Alzheimer's disease (AD) CSF biomarkers may elucidate additional changes that are associated with early AD pathology and enhance our knowledge of the disease.

Methods

The relative abundance of untargeted metabolites was assessed in 161 individuals from the Wisconsin Registry for Alzheimer's Prevention. A metabolome-wide association study (MWAS) was conducted between 269 CSF metabolites and protein biomarkers reflecting brain amyloidosis, tau pathology, neuronal and synaptic degeneration, and astrocyte or microglial activation and neuroinflammation. Linear mixed-effects regression analyses were performed with random intercepts for sample relatedness and repeated measurements and fixed effects for age, sex, and years of education. The metabolome-wide significance was determined by a false discovery rate threshold of 0.05. The significant metabolites were replicated in 154 independent individuals from then Wisconsin Alzheimer's Disease Research Center. Mendelian randomization was performed using genome-wide significant single nucleotide polymorphisms from a CSF metabolites genome-wide association study.

Results

Metabolome-wide association study results showed several significantly associated metabolites for all the biomarkers except Aβ42/40 and IL-6. Genetic variants associated with metabolites and Mendelian randomization analysis provided evidence for a causal association of metabolites for soluble triggering receptor expressed on myeloid cells 2 (sTREM2), amyloid β (Aβ40), α-synuclein, total tau, phosphorylated tau, and neurogranin, for example, palmitoyl sphingomyelin (d18:1/16:0) for sTREM2, and erythritol for Aβ40 and α-synuclein.

Discussion

This study provides evidence that CSF metabolites are associated with AD-related pathology, and many of these associations may be causal.

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.

Protective Alzheimer's disease-associated APP A673T variant predominantly decreases sAPPβ levels in cerebrospinal fluid and 2D/3D cell culture models

The rare A673T variant was the first variant found within the amyloid precursor protein (APP) gene conferring protection against Alzheimer's disease (AD). Thereafter, different studies have discovered that the carriers of the APP A673T variant show reduced levels of amyloid beta (Aβ) in the plasma and better cognitive performance at high age. Here, we analyzed cerebrospinal fluid (CSF) and plasma of APP A673T carriers and control individuals using a mass spectrometry-based proteomics approach to identify differentially regulated targets in an unbiased manner. Furthermore, the APP A673T variant was introduced into 2D and 3D neuronal cell culture models together with the pathogenic APP Swedish and London mutations. Consequently, we now report for the first time the protective effects of the APP A673T variant against AD-related alterations in the CSF, plasma, and brain biopsy samples from the frontal cortex. The CSF levels of soluble APPβ (sAPPβ) and Aβ42 were significantly decreased on average 9-26% among three APP A673T carriers as compared to three well-matched controls not carrying the protective variant. Consistent with these CSF findings, immunohistochemical assessment of cortical biopsy samples from the same APP A673T carriers did not reveal Aβ, phospho-tau, or p62 pathologies. We identified differentially regulated targets involved in protein phosphorylation, inflammation, and mitochondrial function in the CSF and plasma samples of APP A673T carriers. Some of the identified targets showed inverse levels in AD brain tissue with respect to increased AD-associated neurofibrillary pathology. In 2D and 3D neuronal cell culture models expressing APP with the Swedish and London mutations, the introduction of the APP A673T variant resulted in lower sAPPβ levels. Concomitantly, the levels of sAPPα were increased, while decreased levels of CTFβ and Aβ42 were detected in some of these models. Our findings emphasize the important role of APP-derived peptides in the pathogenesis of AD and demonstrate the effectiveness of the protective APP A673T variant to shift APP processing towards the non-amyloidogenic pathway in vitro even in the presence of two pathogenic mutations.

Molecular Landscape of Tourette's Disorder

Tourette's disorder (TD) is a highly heritable childhood-onset neurodevelopmental disorder and is caused by a complex interplay of multiple genetic and environmental factors. Yet, the molecular mechanisms underlying the disorder remain largely elusive. In this study, we used the available omics data to compile a list of TD candidate genes, and we subsequently conducted tissue/cell type specificity and functional enrichment analyses of this list. Using genomic data, we also investigated genetic sharing between TD and blood and cerebrospinal fluid (CSF) metabolite levels. Lastly, we built a molecular landscape of TD through integrating the results from these analyses with an extensive literature search to identify the interactions between the TD candidate genes/proteins and metabolites. We found evidence for an enriched expression of the TD candidate genes in four brain regions and the pituitary. The functional enrichment analyses implicated two pathways ('cAMP-mediated signaling' and 'Endocannabinoid Neuronal Synapse Pathway') and multiple biological functions related to brain development and synaptic transmission in TD etiology. Furthermore, we found genetic sharing between TD and the blood and CSF levels of 39 metabolites. The landscape of TD not only provides insights into the (altered) molecular processes that underlie the disease but, through the identification of potential drug targets (such as FLT3, NAALAD2, CX3CL1-CX3CR1, OPRM1, and HRH2), it also yields clues for developing novel TD treatments.

Liver-Specific Polygenic Risk Score Is Associated with Alzheimer's Disease Diagnosis

BACKGROUND

Our understanding of the pathophysiology underlying Alzheimer's disease (AD) has benefited from genomic analyses, including those that leverage polygenic risk score (PRS) models of disease. The use of functional annotation has been able to improve the power of genomic models.

OBJECTIVE

We sought to leverage genomic functional annotations to build tissue-specific AD PRS models and study their relationship with AD and its biomarkers.

METHODS

We built 13 tissue-specific AD PRS and studied the scores' relationships with AD diagnosis, cerebrospinal fluid (CSF) amyloid, CSF tau, and other CSF biomarkers in two longitudinal cohort studies of AD.

RESULTS

The AD PRS model that was most predictive of AD diagnosis (even without APOE) was the liver AD PRS: n = 1,115; odds ratio = 2.15 (1.67-2.78), p = 3.62×10-9. The liver AD PRS was also statistically significantly associated with cerebrospinal fluid biomarker evidence of amyloid-β (Aβ42:Aβ40 ratio, p = 3.53×10-6) and the phosphorylated tau:amyloid-β ratio (p = 1.45×10-5).

CONCLUSION

These findings provide further evidence of the role of the liver-functional genome in AD and the benefits of incorporating functional annotation into genomic research.

Effect of Pathway-Specific Polygenic Risk Scores for Alzheimer's Disease (AD) on Rate of Change in Cognitive Function and AD-Related Biomarkers Among Asymptomatic Individuals

BACKGROUND

Genetic scores for late-onset Alzheimer's disease (LOAD) have been associated with preclinical cognitive decline and biomarker variations. Compared with an overall polygenic risk score (PRS), a pathway-specific PRS (p-PRS) may be more appropriate in predicting a specific biomarker or cognitive component underlying LOAD pathology earlier in the lifespan.

OBJECTIVE

In this study, we leveraged longitudinal data from the Wisconsin Registry for Alzheimer's Prevention and explored changing patterns in cognition and biomarkers at various age points along six biological pathways.

METHODS

PRS and p-PRSs with and without APOE were constructed separately based on the significant SNPs associated with LOAD in a recent genome-wide association study meta-analysis and compared to APOE alone. We used a linear mixed-effects model to assess the association between PRS/p-PRSs and cognitive trajectories among 1,175 individuals. We also applied the model to the outcomes of cerebrospinal fluid biomarkers in a subset. Replication analyses were performed in an independent sample.

RESULTS

We found p-PRSs and the overall PRS can predict preclinical changes in cognition and biomarkers. The effects of PRS/p-PRSs on rate of change in cognition, amyloid-β, and tau outcomes are dependent on age and appear earlier in the lifespan when APOE is included in these risk scores compared to when APOE is excluded.

CONCLUSION

In addition to APOE, the p-PRSs can predict age-dependent changes in amyloid-β, tau, and cognition. Once validated, they could be used to identify individuals with an elevated genetic risk of accumulating amyloid-β and tau, long before the onset of clinical symptoms.

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.

Genome-wide metabolite quantitative trait loci analysis (mQTL) in red blood cells from volunteer blood donors

The red blood cell (RBC)-Omics study, part of the larger NHLBI-funded Recipient Epidemiology and Donor Evaluation Study (REDS-III), aims to understand the genetic contribution to blood donor RBC characteristics. Previous work identified donor demographic, behavioral, genetic, and metabolic underpinnings to blood donation, storage, and (to a lesser extent) transfusion outcomes, but none have yet linked the genetic and metabolic bodies of work. We performed a genome-wide association (GWA) analysis using RBC-Omics study participants with generated untargeted metabolomics data to identify metabolite quantitative trait loci in RBCs. We performed GWA analyses of 382 metabolites in 243 individuals imputed using the 1000 Genomes Project phase 3 all-ancestry reference panel. Analyses were conducted using ProbABEL and adjusted for sex, age, donation center, number of whole blood donations in the past 2 years, and first 10 principal components of ancestry. Our results identified 423 independent genetic loci associated with 132 metabolites (p < 5×10^-8). Potentially novel locus-metabolite associations were identified for the region encoding heme transporter FLVCR1 and choline and for lysophosphatidylcholine acetyltransferase LPCAT3 and lysophosphatidylserine 16.0, 18.0, 18.1, and 18.2; these associations are supported by published rare disease and mouse studies. We also confirmed previous metabolite GWA results for associations, including N(6)-methyl-L-lysine and protein PYROXD2 and various carnitines and transporter SLC22A16. Association between pyruvate levels and G6PD polymorphisms was validated in an independent cohort and novel murine models of G6PD deficiency (African and Mediterranean variants). We demonstrate that it is possible to perform metabolomics-scale GWA analyses with a modest, trans-ancestry sample size.

GWAS Central: an expanding resource for finding and visualising genotype and phenotype data from genome-wide association studies

The GWAS Central resource gathers and curates extensive summary-level genome-wide association study (GWAS) data and puts a range of user-friendly but powerful website tools for the comparison and visualisation of GWAS data at the fingertips of researchers. Through our continued efforts to harmonise and import data received from GWAS authors and consortia, and data sets actively collected from public sources, the database now contains over 72.5 million P-values for over 5000 studies testing over 7.4 million unique genetic markers investigating over 1700 unique phenotypes. Here, we describe an update to integrate this extensive data collection with mouse disease model data to support insights into the functional impact of human genetic variation. GWAS Central has expanded to include mouse gene-phenotype associations observed during mouse gene knockout screens. To allow similar cross-species phenotypes to be compared, terms from mammalian and human phenotype ontologies have been mapped. New interactive interfaces to find, correlate and view human and mouse genotype-phenotype associations are included in the website toolkit. Additionally, the integrated browser for interrogating multiple association data sets has been updated and a GA4GH Beacon API endpoint has been added for discovering variants tested in GWAS. The GWAS Central resource is accessible at https://www.gwascentral.org/.

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

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

Advances in brain barriers and brain fluids research in 2021: great progress in a time of adversity

This editorial highlights advances in brain barrier and brain fluid research in 2021. It covers research on components of the blood–brain barrier, neurovascular unit and brain fluid systems; how brain barriers and brain fluid systems are impacted by neurological disorders and their role in disease progression; and advances in strategies for treating such disorders.

Vertical Transfer of Metabolites Detectable from Newborn's Dried Blood Spot Samples Using UPLC-MS: A Chemometric Study

The pregnancy period and first days of a newborn's life is an important time window to ensure a healthy development of the baby. This is also the time when the mother and her baby are exposed to the same environmental conditions and intake of nutrients, which can be determined by assessing the blood metabolome. For this purpose, dried blood spots (DBS) of newborns are a valuable sampling technique to characterize what happens during this important mother-child time window. We used metabolomics profiles from DBS of newborns (age 2-3 days) and maternal plasma samples at gestation week 24 and postpartum week 1 from n=664 mother-child pairs of the Copenhagen Prospective Studies on Asthma in Childhood 2010 (COPSAC2010) cohort, to study the vertical mother-child transfer of metabolites. Further, we investigated how persistent the metabolites are from the newborn and up to 6 months, 18 months, and 6 years of age. Two hundred seventy two metabolites from UPLC-MS (Ultra Performance Liquid Chromatography-Mass Spectrometry) analysis of DBS and maternal plasma were analyzed using correlation analysis. A total of 11 metabolites exhibited evidence of transfer (R>0.3), including tryptophan betaine, ergothioneine, cotinine, theobromine, paraxanthine, and N6-methyllysine. Of these, 7 were also found to show persistence in their levels in the child from birth to age 6 years. In conclusion, this study documents vertical transfer of environmental and food-derived metabolites from mother to child and tracking of those metabolites through childhood, which may be of importance for the child's later health and disease.

Metabolic Sex Dimorphism of the Brain at the Gene, Cell, and Tissue Level

The palpable observation in the sex bias of disease prevalence in the CNS has fascinated scientists for several generations. Brain sex dimorphism has been visualized by imaging and analytical tools at the tissue, cellular, and molecular levels. Recent work highlighted the specificity of such sex bias in the brain and its subregions, offering a unique lens through which disease pathogenesis can be investigated. The brain is the largest consumer of energy in the body and provides a unique metabolic environment for diverse lineages of cells. Immune cells are increasingly recognized as an integral part of brain physiology, and their function depends on metabolic homeostasis. This review focuses on metabolic sex dimorphism in brain tissue, resident, and infiltrating immune cells. In this context, we highlight the relevance of recent advances in metabolomics and RNA sequencing technologies at the single cell resolution and the development of novel computational approaches.

MetaboAnalyst 5.0: narrowing the gap between raw spectra and functional insights

Since its first release over a decade ago, the MetaboAnalyst web-based platform has become widely used for comprehensive metabolomics data analysis and interpretation. Here we introduce MetaboAnalyst version 5.0, aiming to narrow the gap from raw data to functional insights for global metabolomics based on high-resolution mass spectrometry (HRMS). Three modules have been developed to help achieve this goal, including: (i) a LC-MS Spectra Processing module which offers an easy-to-use pipeline that can perform automated parameter optimization and resumable analysis to significantly lower the barriers to LC-MS1 spectra processing; (ii) a Functional Analysis module which expands the previous MS Peaks to Pathways module to allow users to intuitively select any peak groups of interest and evaluate their enrichment of potential functions as defined by metabolic pathways and metabolite sets; (iii) a Functional Meta-Analysis module to combine multiple global metabolomics datasets obtained under complementary conditions or from similar studies to arrive at comprehensive functional insights. There are many other new functions including weighted joint-pathway analysis, data-driven network analysis, batch effect correction, merging technical replicates, improved compound name matching, etc. The web interface, graphics and underlying codebase have also been refactored to improve performance and user experience. At the end of an analysis session, users can now easily switch to other compatible modules for a more streamlined data analysis. MetaboAnalyst 5.0 is freely available at https://www.metaboanalyst.ca.