- Research
- Open access
- Published:
The genomic alterations in glioblastoma influence the levels of CSF metabolites
Acta Neuropathologica Communications volume 12, Article number: 13 (2024)
Abstract
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.
Introduction
Glioblastoma, isocitrate dehydrogenase (IDH)-wildtype (GBM) is a highly aggressive brain tumor with poor prognosis and limited treatment options. A major challenge in the treatment of patients with GBM is the lack of reliable, non-invasive diagnostic and prognostic biomarkers for early detection, monitoring, and therapy selection. Cerebrospinal fluid (CSF) is in direct contact with the central nervous system (CNS) and CSF metabolites have the potential to reflect the metabolic alterations associated with tumor presence and progression [42, 46]. Analysis of CSF metabolites is less invasive than tissue biopsy and allows for serial collection at different timepoints over the course of disease. Yet, CSF biomarkers remain underutilized in patients with GBM and have not yet been incorporated into routine clinical practice due to technical challenges and the need for more studies to understand CSF metabolic changes.
Metabolomic analysis has revealed differences in metabolites between viable tumor, necrotic regions, and non-neoplastic regions in GBM tissue samples [12, 18]. Several studies have explored differences in plasma metabolites between GBM patients and controls [1, 14, 33, 38]. Machine learning analysis has also demonstrated changes in blood metabolites before surgery, after surgery, and following chemo-radiation in GBM patients [1]. Another study identified seven plasma molecules, associated with energy metabolism and signaling pathways catalyzing tumor proliferation and invasion, to be biomarkers that are elevated in GBM patients [14].
Although studies of CSF metabolites from GBM patients are limited, metabolomic analysis of CSF has shown promise in the characterization of CNS tumors. For example, differences in CSF metabolites between high-grade gliomas and healthy controls have been identified [34, 46]. Moreover, we have previously identified differences in CSF biomarkers in patients with diffuse gliomas based on IDH1 mutation status, and in the CSF of patients with different types of CNS germ cell tumors [2, 17, 52]. Also, analysis of metabolites in CSF has been shown to predict malignant transformation and leptomeningeal metastasis in glioma patients [22].
In addition to the studies above, associating changes in CSF metabolites to brain tumors, studies suggest that gut microbiome-derived metabolites can influence the blood brain barrier [43, 50, 59]. Imbalances in the relative abundance of gut microbes (dysbiosis) can lead to disease in the CNS [3]. Studies have found differences in the gut microbiome of patients with multiple sclerosis (MS) versus healthy controls [10, 41]. In addition, gut dysbiosis have been described in the setting of brain tumors [13, 32]. In particular, alterations in gut levels of Akkermansia sp. have been associated with the presence of gliomas in humans and mice models, and metabolites of bacterial origin have been associated with neurotoxicity [41, 44].
In this study, we explored the metabolic profile of CSF in patients with GBM and investigated the correlation between CSF metabolites and genomic alterations. Our objective was to identify metabolic signatures in CSF that are associated with tumor presence, and to understand how these metabolites could serve as biomarkers for the diagnosis and monitoring of GBM patients.
Methods
Patients
CSF samples from 44 patients were included in the study: n = 31 from GBM patients and n = 13 from patients without a history of cancer (controls). Samples were collected via lumbar puncture (n = 18), sulcus sampling during resection (n = 17), ventricular catheter (n = 7), reservoir/shunt (n = 1), or cisterna magna (n = 1) (Fig. 1A).
CSF processing
CSF samples were processed within 3 h of collection and centrifuged twice at 1,000 × g for 10 min at 4 °C. The cell pellet was discarded, and the supernatant was immediately stored at -80 °C until the time of analysis. 100 µL of each CSF sample were used for metabolomics analysis. CSF via lumbar puncture was collected utilizing the adult LP tray (Cat. 4306 C, CareFusion, IL, USA).
Targeted metabolomics
Liquid chromatography-mass spectrometry (LC-MS) Single Reaction Monitoring (SRM) was used to characterize metabolites in CSF. Metabolites were measured using three different chromatographic methods. For each method, metabolites were normalized against the spiked internal standards and the data were log2-transformed (see Sup. Material 6 for details).
Statistical analysis
The relative abundance of 125 metabolites was measured as part of a targeted metabolite profiling analysis. The resulting counts were normalized to an internal standard and log2-transformed. The data were additionally z-transformed for each metabolite across all samples for heatmap generation. Welch’s t-test was used for comparisons of metabolite levels between two groups. Kruskal-Wallis test with post hoc Dunn’s test for multiple comparisons was used for comparisons of metabolite levels between 3 groups. All analyses were performed using Python 3.8 and R Statistical Software (v4.2.2; R Core Team 2021). Significance calculations were performed using the Python library scikit-learn [45]. The University of Texas MD Anderson Cancer Center Department of Bioinformatics next-generation clustered heat map (NGCHM) software in R was used to generate heatmaps [6]. Sparse partial least squares-discriminant analysis (sPLS-DA) plots were generated using the mixOmics R toolkit [47]. Boxplots were generated using the R package ggpubr [28] and volcano plots using Matplotlib [21] in Python. Kaplan-Meier plots were generated using the kaplanmeier Python library [53], and cutoff values for determining high and low groups were found using the function get.cutoff() from the CutoffFinder R file [7].
Mutation analysis
A subset of matched tumor samples was analyzed for genomic alterations by a targeted next-generation sequencing (NGS) panel interrogating 205 cancer-related genes for mutations and 26 genes for rearrangements (FoundationOne; Foundation Medicine, Inc.). This information was obtained from the patient’s electronical medical record.
Results
Patients and CSF samples
This study included CSF from 31 patients with IDH-wildtype glioblastoma (GBM), 23 men and 8 women, ranging in age from 18 to 76 years. Thirteen (13) patients had CSF samples collected pre-treatment, while 18 patients had CSF samples collected post-treatment (Fig. 1A, 1B). Somatic mutation analysis of tumor tissue was only available for the 18 patients with post-treatment CSF samples (Fig. 1B). In addition, CSF from 13 patients, (4 male and 9 female), ranging in age from 25 to 77 years, with non-neoplastic CNS diseases (i.e., hydrocephalus, stroke, aqueduct stenosis, migraine, arachnoid cyst, colloid cyst, and trauma) was used as control (Fig. 1A). The survival and event timeline of GBM patients in the study are illustrated in Sup. Material 1.
Differences in metabolites between pre-treatment GBM and post-treatment GBM samples
125 metabolites were identified in CSF as part of the targeted metabolic analysis. Statistical comparisons of metabolite levels were performed between GBM pre-treatment (GBM pre-tx) and control patients, as well as between GBM pre-treatment and post-treatment samples. sPLS-DA demonstrated differences between pre-treatment and control, as well as between pre-treatment and post-treatment GBM CSF samples (Fig. 2A), (Fig. 2B). Therefore, pre- and post-treatment CSF samples were analyzed separately.
Differences in CSF metabolites based on sex
There is a difference in the proportion of male and female patients between the pre- and post-treatment GBM cohorts, with a near-even split (46% female) in the pre-treatment cohort and disproportionately more male patients (11% female) in the post-treatment cohort. To confirm that the metabolic differences between pre- and post-treatment CSF samples are not due to differences in the proportion of female patients, we performed significance and fold change calculations on the GBM cohort based on sex. We identified a small number of metabolites (2-methylbutyrylcarnitine, methyl methoxyacetate, 2-hydroxyglutarate, cytosine, cholesterol sulfate, methylparaben) that are present at different levels between male and female patients in the GBM cohort (pre- and post-treatment samples) (Sup. Material 2A). Similarly, we identified a small number of metabolites (2-methylbutyrylcarnitine, 2-hydroxyglutarate, allantoin, glycine) that are different between male and female patients in the pre-treatment GBM cohort (Sup. Material 2B). To exclude variability due to patient’s sex, we compared pre-treatment and post-treatment CSF from male patients only, and we identified several metabolites that are significantly different between male pre- and post-treatment GBM CSF (3-Hydroxy-3-methylglutaric acid, stearic acid, uridine, spaglumic acid, hydroxyisocaproic acid, myristoleic acid, methyl hippurate, aminobutanal, alanine, kynurenine, s-ribosyl-L-homocysteine, alpha-ketoglutarate) (Sup. Material 2C). Therefore, the differences in CSF metabolites between pre-treatment and post-treatment GBM samples are not due only to sex.
Differences in metabolites between pre-treatment GBM and control CSF samples
sPLS-DA demonstrated differences in CSF metabolites between GBM (pre-treatment) and control samples (Fig. 2C). Supervised hierarchical clustering based on the squared Euclidean distance (Ward’s method) revealed differences in overall metabolite abundance between GBM pre-tx and controls (Fig. 2D). Volcano plots illustrated fold-change and probability values for each identified metabolite and highlighted differentially abundant metabolites. Several metabolites showed statistically significantly different levels between GBM pre-tx and control CSF (Fig. 2E). In particular, carnitine, 2-methylbutyrylcarnitine, shikimate, aminobutanal, uridine, N-acetylputrescine, and farnesyl diphosphate were present at significantly higher levels in the CSF of GBM pre-tx patients compared to controls (Fig. 3).
Correlation between CSF metabolite levels and somatic mutations in GBM
Somatic mutation analysis of tumor tissue was available for the 18 GBM patients with CSF collected post-treatment (Fig. 1B). The most common mutations identified were TERT promoter (TERTp), CDKN2A/B, NF1, PTEN, TP53, and EGFR. Global analysis of metabolites shows a statistically significant difference in CSF metabolite levels between patients with TP53-wildtype and TP53-mutant GBM (Fig. 4), as well as PTEN-mutant and PTEN-wildtype GBM (Fig. 5). Only minor significant differences in metabolites were identified according to TERTp, CDKN2A/B, NF1 or EGFR mutation status.
Analysis of the levels of carnitine compounds revealed lower levels of these metabolites (i.e., carnitine, propionylcarnitine, 2-methylbutyrylcarnitine, isobutyryl-L-carnitine, and deoxycarnitine) in CSF from patients with TP53-mutant versus TP53-wildtype GBM (Fig. 6). When compared to the control CSF, TP53-mutant samples did not show significantly different levels for four of the five carnitine compounds (exception: 2-methylbutyrylcarnitine). In contrast, the level of all carnitine compounds in the CSF from patients with TP53-wildtype GBM was significantly higher across the board, compared to that in CSF from control patients or TP53-mutant GBM patients (Fig. 6).
We next analyzed the levels of clinically-relevant metabolites that are routinely measured by magnetic resonance spectroscopy (MRS) in patients with GBM, including lactate, γ-aminobutyric acid (GABA), and choline. The data showed significantly lower levels of lactate, GABA and choline in CSF from patients with TP53-mutant GBM compared to TP53-wildtype. The data also showed higher levels of lactate, GABA, and choline in the CSF of patients with PTEN-mutant GBM compared to PTEN-wildtype GBM (Fig. 7).
Correlation between CSF metabolites and survival in patients with GBM
We analyzed the correlation between metabolite levels in CSF and overall survival (OS) for the 9 pre-treatment GBM samples (Sup. Material 4). To evaluate if metabolites levels correlate with overall survival, we separated patients into “high” or “low” groups for each metabolite based on the cutoff value calculated from CutoffFinder’s fit of mixture model [7]. Kaplan-Meier plots were generated for each metabolite. OS was compared between patients with “high” or “low” levels for each of the 125 metabolites measured in the study. Three metabolites showed a statistically significant association with OS in patients with GBM (2-methylbutyrylcarnitine, aminobutanal, and acetylcholine) (Fig. 8).
Metabolites and age
We performed a linear regression analysis between age and metabolite levels, for each metabolite, for the entire cohort of glioblastoma samples (pre- and post-treatment) and no significant correlation was observed (R^2 value ranged from 0.15 to 0.000055). We performed the same analysis with the pre-treatment samples and did not observe a significant correlation between age and metabolite levels for any of the metabolites analyzed (R^2 value ranged from 0.52 to 8.12E-07). We also performed the same analysis with the post-treatment cohort and did not observe a significant correlation between age and metabolite levels for any of the metabolites analyzed (R^2 value ranged from 0.30 to 1.76E-06).
Discussion
Our results in this study and previous studies demonstrate differences in CSF metabolites between patients with diffuse gliomas and individuals without brain cancer [2, 22, 34]. Differences in levels of CSF metabolites may be due to the metabolic activity of cancer cells and/or the effect of tumor cells on overall brain activity [40, 42]. In this study, we found several metabolites that are present at significantly different levels in CSF of patients with GBM compared to controls. The data also show that GBM mutations may influence the levels of metabolites in CSF, including metabolites that are routinely analyzed by MRS with advanced brain MRI scans. Moreover, our current and previous data show that a metabolite considered to be of bacterial origin is significantly altered in the CSF of patients with GBM [2].
Carnitine association between GBM and lipid metabolism
Carnitine and 2-methylbutyrylcarnitine were found to be present at significantly higher levels in the CSF of GBM patients compared to controls (Fig. 3A, 3B). GBM cells have been shown to exhibit elevated levels of carnitine [15], and although this was found in our study, it is unclear what role the GBM size, GBM location, and other factors may play. Carnitine and acyl-carnitines are crucial components of the carnitine shuttle system (CSS), which facilitates the transport of fatty acids into the mitochondria for fatty acid oxidation (FAO) (Sup. Material 3A) [8, 20, 35]. Recent studies have implicated FAO as an important energy production pathway in GBM as well as a key mediator of GBM plasticity and adaptability [26, 31, 49]. An increase in the CSF levels of carnitine found in this study may indicate increased reliance on lipid metabolism in GBM. Organic cation transporter OCTN2 (SLC22A5) is the primary transporter of dietary carnitine throughout the body, and OCTN2 is highly expressed in the brain of GBM patients compared to healthy brains [15, 24, 29]. Consequently, the higher carnitine levels observed in the CSF of the GBM patients evaluated here may be facilitated by upregulation of OCTN2, which in turn could potentially be a response to greater demand for carnitine by enhanced fatty acid oxidation in GBM.
Shikimate levels and the gut microbiome
In our previous study of CSF metabolites, we found an association between higher CSF levels of shikimate and IDH-wildtype gliomas compared to controls [2], a finding recapitulated in Fig. 3C. CSF shikimate levels have also been found to be elevated in patients with autism spectrum disorder [5]. Shikimate is an intermediate in the shikimate pathway. This pathway is not present in animals but it is used by plants, bacteria, algae, and other microorganisms in the biosynthesis of aromatic compounds (Sup. Material 3B) [19]. While the shikimate pathway is active in some bacteria present in the human gut microbiome, transcriptome analysis shows that most gut bacteria do not possess a complete shikimate pathway [36]. However, computational modeling reveals that the gut bacteria Akkermansia muciniphila expresses all the genes encoding the necessary enzymes in the shikimate pathway and accounts for a significant portion of the gene expression associated with the shikimate pathway in the human gut microbiome. These data suggest that A. muciniphila could be responsible for a significant amount of total shikimate produced by human gut microbiome [36]. We have previously identified higher levels of the Akkermansia genus in the gut microbiome of glioma patients and a glioma mouse model, compared to controls [44]. Therefore, the elevated CSF levels of shikimate in GBM patients identified in the current study may be influenced by the gut microbiome, in particular Akkermansia sp., suggesting a possible relationship between the gut-brain axis and CSF metabolites. Additionally, recent studies demonstrate unique microbial signatures within cancer types, including GBM [39]. To determine the source of the CSF shikimate identified in this study, future studies could determine whether shikimate-producing bacteria are present in human GBM tissue [39]. In order to detect remnants of shikimate-producing bacteria in the CSF, metagenomic sequencing could be employed, in a manner similar to but more specific than that used clinically in patients with meningitis or encephalitis [58].
Uridine levels correlate with tumor samples
We found uridine CSF levels to be greater in pre-treatment GBM samples compared to control samples. Uridine is a pyrimidine molecule that is an intermediate in nucleoside synthesis, particularly in the catabolism of L-glutamine to B-alanine (Sup. Material 3C) [62]. Consistent with our results, uridine concentration was found to be higher in tissue samples of GBM, IDH-wildtype tumors compared to controls [25].
Tumor mutations influence CSF metabolite levels
Previous studies have demonstrated differences in the levels of metabolites based on genetic alterations in tumors. For example, diffuse gliomas have been shown to exhibit distinct tissue metabolite profiles depending on IDH1 mutation status [55]. Moreover, our previous study revealed higher levels of D-2-hydroxyglutarate in the CSF of patients with IDH-mutant gliomas compared to patients with IDH-wildtype gliomas [17]. In the current study, we have identified novel differences in CSF metabolites between patients with TP53-mutant and TP53-wildtype GBM. CSF from patients with TP53-wildtype GBM contained higher levels of several metabolites including lactate, GABA, choline, carnitine, and carnitine derivates. In contrast, the levels of lactate, GABA, and choline were increased in the CSF of patients with GBM-PTEN-mutant, compared to GBM-PTEN-wildtype (Fig. 5).
Carnitine and derived compounds are more abundant in CSF from patients with TP53-wildtype GBM
We identified carnitine and carnitine derivatives as key metabolic biomarkers in the CSF of patients with GBM. Carnitine and its acyl derivatives emerged as significant factors in our dataset, with increased abundance in the CSF of GBM patients, particularly those with TP53-wildtype status. Carnitine, acylcarnitines 2-methylbutyrylcarnitine, propionylcarnitine, isobutyryl-L-carnitine, and deoxycarnitine were significantly more abundant in TP53-wildtype, compared to TP53-mutant, GBM patients (Fig. 6). In GBM patients, TP53 potentially plays an indirect yet critical role in regulating FAO and may subsequently affect carnitine and acylcarnitine levels in CSF. More specifically, p53-responsive elements have been identified in the first intron of the carnitine palmitoyltransferase IC (CPT1C) gene; studies in cell lines have shown that p53 can regulate CPT1C expression [48]. The effects of increased CPT1C expression on tumor progression and FAO activity is well-documented in the literature [11, 57, 60, 63]. While the exact mechanism for reduced levels of carnitine compounds in the CSF of patients with TP53-mutant GBM remains to be elucidated, impaired TP53 regulation of CPT1C expression in TP53-mutant GBM may be a contributing factor. When the TP53 mutations occur, FAO is reduced [48], and the tumor cells’ ability to use FAO as an alternative energy production pathway might be affected. This may result in a reduced demand for carnitine in TP53-mutant GBM and lower concentrations of carnitine and carnitine derivatives in the CSF. The present findings support this hypothesis, as the relative levels of carnitine, propionylcarnitine, and isobutyryl-L-carnitine were significantly increased in TP53-wildtype GBM but were not significantly different from control patients or patients with TP53-mutant GBM (Fig. 6).
MRS biomarkers for TP53 and PTEN mutation status
The metabolites lactate, GABA, and choline are clinically relevant because they are measured in imaging studies through magnetic resonance spectroscopy (MRS) [56]. Some studies show that tumor mutation status can potentially influence MRS results [30]. Our results reveal that the levels of these metabolites can be influenced by the mutations present in the tumor (Fig. 7). Lactate, GABA, and choline were present at significantly lower levels in the CSF of TP53-mutant GBM patients when compared to that of TP53-wildtype GBM patients (Fig. 7A and 7C). On the contrary, the relative values of each of the three MRS-measurable metabolites were significantly higher in the CSF of PTEN-mutant GBM patients compared to that of PTEN-wildtype GBM patients (with the exception of lactate, which follows a similar trend) (Fig. 7D and 7F). Absence of p53 has been associated with downregulation of choline levels in human embryonic stem cells [65], suggesting a potential link between choline metabolism and p53.
TP53 and the warburg effect
In addition, TP53 is known to regulate glycolysis and influence the Warburg effect by negatively regulating lactate dehydrogenase [64]. Similarly, it has been proposed that mutant TP53 can stimulate the Warburg effect in cultured cells [61]. Interestingly, we observed lower lactate levels in the CSF of patients with mutant TP53. The relationship between lower lactate levels in CSF and the predicted effects of mutant TP53 in regulating the Warburg effect is unclear and requires further investigation. Given that the levels of lactate, GABA, and choline levels appear to be influenced by GBM mutations, it may be important to consider tumor mutations when interpreting MRS results. To evaluate this further, it would be necessary to correlate the CSF metabolite findings with MRS imaging and perform receiver-operator curve analysis to detect the sensitivity and specificity of a given cut-off on MRS that reliably predicts TP53 and PTEN mutation status when GBM tissue is not available.
Metabolites associated with sex
Differences between male and female patients with glioblastoma have been described. For example, the incidence of glioblastoma is 1.6 times higher in men than women [51]. Also, men have worse OS than women. Although differences in the proportion of MGMT promoter methylation between male and female and female with glioblastoma have been proposed [9], it appears that the difference is not statistically significant [16]. In addition, studies have shown differences in the expression of glycolytic genes between males and females, with potential association to overall survival [23]. In our analysis, we identified some CSF metabolites that were significantly different between male and patients, consistent with the idea that physiological and hormonal differences associated with sex can influence the levels of metabolites in CSF in patients with GBM. This is important when trying to identify metabolites that can serve as GBM biomarkers, because it might be required to use different cutoffs for males and females, when using a CSF metabolite as biomarker to identify patients with GBM.
Metabolites associated with overall survival
We identified four metabolites that showed a statistically significant association with OS (Fig. 8). Increased levels of 2-methylbutyrylcarnitine and aminobutanal were found to be associated with the presence of GBM, and higher levels were associated with worse OS (Figs. 3B and 3D and 8A and 8B). Aminobutanal is involved in the synthesis of the neurotransmitter GABA [27] and GABA metabolism has been linked to GBM survival [4]. Higher levels of acetylcholine in CSF were found to be associated with worse OS (Fig. 8C). Acetylcholine is a well-defined neurotransmitter with a wide variety of roles in the CNS [37]. GBM invasion was found to be significantly enhanced in brain regions with activation of acetylcholine receptors, and it has been suggested that acetylcholine autocrine signaling facilitates GBM invasion through brain tissue [54]. The same study also found that increased acetylcholine receptor expression in GBM correlated with worse OS.
Limitations
One limitation of this study is that only CSF samples acquired post-treatment have tumor mutation data available for analysis. Another limitation is the small number of female patients in the post-treatment group. Although, the lower proportion of female patients in the post-treatment GBM group is a limitation, comparison of only the male GBM patients before and after treatment showed significant differences in CSF metabolites. This supports the result that there are differences between the pre- and post-treatment GBM samples that are not due to differences in the number of female patients between the groups. Additional studies with a larger number of patients, including CSF samples acquired before therapy, will provide further confirmation for some of our results. Also, how the time between treatment and CSF collection influences the levels of CSF metabolites is unknown. Despite these limitations, the results of this study increase our understanding of CSF metabolites in GBM patients.
Conclusion
This study demonstrates differences in CSF metabolites between GBM patients and controls. In particular, the elevated levels of carnitine and 2-methylbutyrylcarnitine point to the role of lipid metabolism in GBM biology. Moreover, the elevated CSF shikimate levels suggest a potentially novel link between GBM and the gut microbiome. In addition, the findings highlight the influence of tumor mutations on CSF metabolite levels, including those that can be independently assessed with non-invasive magnetic resonance spectroscopy.
Data availability
All data generated or analyzed during this study are included in this published article and its supplementary information files.
References
Aboud O, Liu YA, Fiehn O, Brydges C, Fragoso R, Lee HS, Riess J, Hodeify R, Bloch O (2023) Application of machine learning to Metabolomic Profile characterization in Glioblastoma patients undergoing concurrent chemoradiation. Metabolites 13. https://doi.org/10.3390/metabo13020299
Ballester LY, Lu G, Zorofchian S, Vantaku V, Putluri V, Yan Y, Arevalo O, Zhu P, Riascos RF, Sreekumar A et al (2018) Analysis of cerebrospinal fluid metabolites in patients with primary or metastatic central nervous system tumors. Acta Neuropathol Commun 6: 85 https://doi.org/10.1186/s40478-018-0588-z
Belizario JE, Faintuch J (2018) Microbiome and Gut Dysbiosis. Exp Suppl 109:459–476. https://doi.org/10.1007/978-3-319-74932-7_13
Blanchart A, Fernando R, Haring M, Assaife-Lopes N, Romanov RA, Andang M, Harkany T, Ernfors P (2017) Endogenous GAB(AA) receptor activity suppresses glioma growth. Oncogene 36:777–786. https://doi.org/10.1038/onc.2016.245
Brister D, Werner BA, Gideon G, McCarty PJ, Lane A, Burrows BT, McLees S, Adelson PD, Arango JI, Marsh Wet al et al (2022) Central Nervous System Metabolism in Autism, Epilepsy and Developmental Delays: A Cerebrospinal Fluid Analysis. Metabolites 12: https://doi.org/10.3390/metabo12050371
Broom BM, Rohrdanz MA (2022) NGCHM: Utilities for Creating Next Generation Clustered Heat Maps. R package version 0.13.1 edn, City
Budczies J, Klauschen F, Sinn BV, Győrffy B, Schmitt WD, Darb-Esfahani S, Denkert C (2012) Cutoff finder: a comprehensive and straightforward web application enabling rapid biomarker cutoff optimization. PLoS ONE 7:e51862. https://doi.org/10.1371/journal.pone.0051862
Carracedo A, Cantley LC, Pandolfi PP (2013) Cancer metabolism: fatty acid oxidation in the limelight. Nat Rev Cancer 13:227–232. https://doi.org/10.1038/nrc3483
Carrano A, Juarez JJ, Incontri D, Ibarra A, Guerrero Cazares H (2021) Sex-specific differences in Glioblastoma. Cells 10. https://doi.org/10.3390/cells10071783
Chen J, Chia N, Kalari KR, Yao JZ, Novotna M, Paz Soldan MM, Luckey DH, Marietta EV, Jeraldo PR, Chen Xet al et al (2016) Multiple sclerosis patients have a distinct gut microbiota compared to healthy controls. Sci Rep 6:28484. https://doi.org/10.1038/srep28484
Chen T, Wu G, Hu H, Wu C (2020) Enhanced fatty acid oxidation mediated by CPT1C promotes gastric cancer progression. J Gastrointest Oncol 11:695–707. https://doi.org/10.21037/jgo-20-157
Clish CB (2015) Metabolomics: an emerging but powerful tool for precision medicine. Cold Spring Harb Mol Case Stud 1:a000588. https://doi.org/10.1101/mcs.a000588
Dono A, Nickles J, Rodriguez-Armendariz AG, McFarland BC, Ajami NJ, Ballester LY, Wargo JA, Esquenazi Y (2022) Glioma and the gut-brain axis: opportunities and future perspectives. Neurooncol Adv 4:vdac054. https://doi.org/10.1093/noajnl/vdac054
Ferrasi AC, Puttini R, Galvani AF, Hamamoto Filho PT, Delafiori J, Argente VD, de Oliveira AN, Dias-Audibert FL, Catharino RR, Silva OCet al et al (2023) Metabolomics Approach reveals important Glioblastoma plasma biomarkers for Tumor Biology. Int J Mol Sci 24. https://doi.org/10.3390/ijms24108813
Fink MA, Paland H, Herzog S, Grube M, Vogelgesang S, Weitmann K, Bialke A, Hoffmann W, Rauch BH, Schroeder HWS et al (2019) -Carnitine-Mediated Tumor Cell Protection and Poor Patient Survival Associated with OCTN2 Overexpression in Glioblastoma Multiforme. Clin Cancer Res 25: 2874–2886 Doi 10.1158/1078– 0432.CCR-18-2380
Franceschi E, Tosoni A, Minichillo S, Depenni R, Paccapelo A, Bartolini S, Michiara M, Pavesi G, Urbini B, Crisi Get al et al (2018) The prognostic roles of gender and O6-Methylguanine-DNA methyltransferase methylation status in Glioblastoma patients: the female power. World Neurosurg 112:e342–e347. https://doi.org/10.1016/j.wneu.2018.01.045
Fujita Y, Nunez-Rubiano L, Dono A, Bellman A, Shah M, Rodriguez JC, Putluri V, Kamal AHM, Putluri N, Riascos RF al (2022) IDH1 p.R132H ctDNA and D-2-hydroxyglutarate as CSF biomarkers in patients with IDH-mutant gliomas. J Neurooncol 159:261–270. https://doi.org/10.1007/s11060-022-04060-1
He W, Edney MK, Paine SML, Griffiths RL, Scurr DJ, Rahman R, Kim DH (2023) Untargeted metabolomic characterization of Glioblastoma Intra-tumor Heterogeneity using OrbiSIMS. Anal Chem 95:5994–6001. https://doi.org/10.1021/acs.analchem.2c05807
Herrmann KM, Weaver LM (1999) THE SHIKIMATE PATHWAY. Annu Rev Plant Physiol Plant Mol Biol 50:473–503. https://doi.org/10.1146/annurev.arplant.50.1.473
Houten SM, Wanders RJ (2010) A general introduction to the biochemistry of mitochondrial fatty acid β-oxidation. J Inherit Metab Dis 33:469–477. https://doi.org/10.1007/s10545-010-9061-2
Hunter JD (2007) Matplotlib: a 2D graphics environment. IEEE COMPUTER SOC, City, pp 90–95
Im JH, Yoo BC, Lee JH, Kim KH, Kim TH, Lee KY, Kim JH, Park JB, Kwon JW, Shin SH et al (2020) Comparative cerebrospinal fluid metabolites profiling in glioma patients to predict malignant transformation and leptomeningeal metastasis with a potential for preventive personalized medicine. EPMA J 11: 469–484 https://doi.org/10.1007/s13167-020-00211-4
Ippolito JE, Yim AK, Luo J, Chinnaiyan P, Rubin JB (2017) Sexual dimorphism in glioma glycolysis underlies sex differences in survival. JCI Insight 2. https://doi.org/10.1172/jci.insight.92142
Juraszek B, Czarnecka-Herok J, Nałęcz KA (2021) Glioma cells survival depends both on fatty acid oxidation and on functional carnitine transport by SLC22A5. J Neurochem 156:642–657. https://doi.org/10.1111/jnc.15124
Kampa JM, Kellner U, Marsching C, Ramallo Guevara C, Knappe UJ, Sahin M, Giampa M, Niehaus K, Bednarz H (2020) Glioblastoma multiforme: metabolic differences to peritumoral tissue and IDH-mutated gliomas revealed by mass spectrometry imaging. Neuropathology 40:546–558. https://doi.org/10.1111/neup.12671
Kant S, Kesarwani P, Prabhu A, Graham SF, Buelow KL, Nakano I, Chinnaiyan P (2020) Enhanced fatty acid oxidation provides glioblastoma cells metabolic plasticity to accommodate to its dynamic nutrient microenvironment. Cell Death Dis 11:253. https://doi.org/10.1038/s41419-020-2449-5
Kanwal S, Incharoensakdi A (2020) GABA synthesis mediated by γ-aminobutanal dehydrogenase in Synechocystis sp. PCC6803 with disrupted glutamate and α-ketoglutarate decarboxylase genes. Plant Sci 290:110287. https://doi.org/10.1016/j.plantsci.2019.110287
Kassambara A (2023) ggpubr: ‘ggplot2’ Based Publication Ready Plots. City
Kido Y, Tamai I, Ohnari A, Sai Y, Kagami T, Nezu J, Nikaido H, Hashimoto N, Asano M, Tsuji A (2001) Functional relevance of carnitine transporter OCTN2 to brain distribution of L-carnitine and acetyl-L-carnitine across the blood-brain barrier. J Neurochem 79:959–969. https://doi.org/10.1046/j.1471-4159.2001.00621.x
Li Y, Ji F, Jiang Y, Zhao T, Xu C (2018) Correlation analysis of expressions of PTEN and p53 with the value obtained by magnetic resonance spectroscopy and apparent diffusion coefficient in the tumor and the tumor-adjacent area in magnetic resonance imaging for glioblastoma. J BUON 23:391–397
Lin H, Patel S, Affleck VS, Wilson I, Turnbull DM, Joshi AR, Maxwell R, Stoll EA (2017) Fatty acid oxidation is required for the respiration and proliferation of malignant glioma cells. Neuro Oncol 19:43–54. https://doi.org/10.1093/neuonc/now128
Lin B, Ye Z, Ye Z, Wang M, Cao Z, Gao R, Zhang Y (2023) Gut microbiota in brain tumors: an emerging crucial player. CNS Neurosci Ther 29(Suppl 1):84–97. https://doi.org/10.1111/cns.14081
Liu A, Aboud O, Dahabiyeh LA, Bloch O, Fiehn O (2023) A pilot study on metabolomic characterization of human glioblastomas and patient plasma. Res Sq: Doi. https://doi.org/10.21203/rs.3.rs-2662020/v1
Locasale JW, Melman T, Song S, Yang X, Swanson KD, Cantley LC, Wong ET, Asara JM (2012) Metabolomics of human cerebrospinal fluid identifies signatures of malignant glioma. Mol Cell Proteomics 11(014688):M111. https://doi.org/10.1074/mcp.M111.014688
Longo N, Frigeni M, Pasquali M (2016) Carnitine transport and fatty acid oxidation. Biochim Biophys Acta 1863:2422–2435. https://doi.org/10.1016/j.bbamcr.2016.01.023
Mesnage R, Antoniou MN (2020) Computational modelling provides insight into the effects of glyphosate on the shikimate pathway in the human gut microbiome. Curr Res Toxicol 1:25–33. https://doi.org/10.1016/j.crtox.2020.04.001
Mesulam MM (2009) Acetylcholine neurotransmission in CNS. In: Larry RS (ed) Encyclopedia of Neuroscience. Academic Press, City, pp 1–4
Muller Bark J, Karpe AV, Doecke JD, Leo P, Jeffree RL, Chua B, Day BW, Beale DJ, Punyadeera C (2023) A pilot study: metabolic profiling of plasma and saliva samples from newly diagnosed glioblastoma patients. Cancer Med 12:11427–11437. https://doi.org/10.1002/cam4.5857
Nejman D, Livyatan I, Fuks G, Gavert N, Zwang Y, Geller LT, Rotter-Maskowitz A, Weiser R, Mallel G, Gigi Eet al et al (2020) The human tumor microbiome is composed of tumor type-specific intracellular bacteria. Science 368:973–980. https://doi.org/10.1126/science.aay9189
Niu D, Sun P, Zhang F, Song F (2022) Metabonomic analysis of cerebrospinal fluid in epilepsy. Ann Transl Med 10:449. https://doi.org/10.21037/atm-22-1219
Ntranos A, Park HJ, Wentling M, Tolstikov V, Amatruda M, Inbar B, Kim-Schulze S, Frazier C, Button J, Kiebish MA et al (2022) Bacterial neurotoxic metabolites in multiple sclerosis cerebrospinal fluid and plasma. Brain 145: 569–583 https://doi.org/10.1093/brain/awab320
Panyard DJ, Kim KM, Darst BF, Deming YK, Zhong X, Wu Y, Kang H, Carlsson CM, Johnson SC, Asthana Set al et al (2021) Cerebrospinal fluid metabolomics identifies 19 brain-related phenotype associations. Commun Biol 4:63. https://doi.org/10.1038/s42003-020-01583-z
Parker A, Fonseca S, Carding SR (2020) Gut microbes and metabolites as modulators of blood-brain barrier integrity and brain health. Gut Microbes 11:135–157. https://doi.org/10.1080/19490976.2019.1638722
Patrizz A, Dono A, Zorofchian S, Hines G, Takayasu T, Husein N, Otani Y, Arevalo O, Choi HA, Savarraj Jet al et al (2020) Glioma and temozolomide induced alterations in gut microbiome. Sci Rep 10:21002. https://doi.org/10.1038/s41598-020-77919-w
Pedregosa F, Varoquaux G, Gramfort A, Michel V, Thirion B, Grisel O, Blondel M, Prettenhofer P, Weiss R Dubourg V (2011) Scikit-learn: machine learning in Python. J Mach Learn Res 12: 2825–2830
Rogachev AD, Alemasov NA, Ivanisenko VA, Ivanisenko NV, Gaisler EV, Oleshko OS, Cheresiz SV, Mishinov SV, Stupak VV, Pokrovsky AG (2021) Correlation of metabolic profiles of plasma and Cerebrospinal Fluid of High-Grade Glioma patients. Metabolites 11. https://doi.org/10.3390/metabo11030133
Rohart F, Gautier B, Singh A, Lê Cao KA (2017) mixOmics: an R package for ‘omics feature selection and multiple data integration. PLoS Comput Biol 13:e1005752. https://doi.org/10.1371/journal.pcbi.1005752
Sanchez-Macedo N, Feng J, Faubert B, Chang N, Elia A, Rushing EJ, Tsuchihara K, Bungard D, Berger SL, Jones RG al (2013) Depletion of the novel p53-target gene carnitine palmitoyltransferase 1 C delays tumor growth in the neurofibromatosis type I tumor model. Cell Death Differ 20:659–668. https://doi.org/10.1038/cdd.2012.168
Shim JK, Choi S, Yoon SJ, Choi RJ, Park J, Lee EH, Cho HJ, Lee S, Teo WY, Moon JH al (2022) Etomoxir, a carnitine palmitoyltransferase 1 inhibitor, combined with temozolomide reduces stemness and invasiveness in patient-derived glioblastoma tumorspheres. Cancer Cell Int 22:309. https://doi.org/10.1186/s12935-022-02731-7
Shreiner AB, Kao JY, Young VB (2015) The gut microbiome in health and in disease. Curr Opin Gastroenterol 31:69–75. https://doi.org/10.1097/MOG.0000000000000139
Sun T, Plutynski A, Ward S, Rubin JB (2015) An integrative view on sex differences in brain tumors. Cell Mol Life Sci 72:3323–3342. https://doi.org/10.1007/s00018-015-1930-2
Takayasu T, Shah M, Dono A, Yan Y, Borkar R, Putluri N, Zhu JJ, Hama S, Yamasaki F Tahara H (2020) cerebrospinal fluid ctDNA and metabolites are informative biomarkers for the evaluation of CNS germ cell tumors. Sci Rep 10: 14326 https://doi.org/10.1038/s41598-020-71161-0
Taskesen E (2020) Kaplanmeier is an python library to create survival curves using kaplan-meier, and compute the log-rank test. 0.1.0 edn, City
Thompson EG, Sontheimer H (2019) Acetylcholine receptor activation as a modulator of Glioblastoma Invasion. Cells 8. https://doi.org/10.3390/cells8101203
Trautwein C, Zizmare L, Mäurer I, Bender B, Bayer B, Ernemann U, Tatagiba M, Grau SJ, Pichler BJ Skardelly M (2022) tissue metabolites in diffuse glioma and their modulations by IDH1 mutation, histology, and treatment. JCI Insight 7: https://doi.org/10.1172/jci.insight.153526
Verma A, Kumar I, Verma N, Aggarwal P, Ojha R (2016) Magnetic resonance spectroscopy - revisiting the biochemical and molecular milieu of brain tumors. BBA Clin 5:170–178. https://doi.org/10.1016/j.bbacli.2016.04.002
Wang CY, Wang CH, Mai RT, Chen TW, Li CW, Chao CH (2022) Mutant p53-microRNA-200c-ZEB2-Axis-Induced CPT1C Elevation contributes to metabolic reprogramming and tumor progression in basal-like breast cancers. Front Oncol 12:940402. https://doi.org/10.3389/fonc.2022.940402
Wilson MR, Sample HA, Zorn KC, Arevalo S, Yu G, Neuhaus J, Federman S, Stryke D, Briggs B, Langelier Cet al et al (2019) Clinical metagenomic sequencing for diagnosis of Meningitis and Encephalitis. N Engl J Med 380:2327–2340. https://doi.org/10.1056/NEJMoa1803396
Xie J, Bruggeman A, De Nolf C, Vandendriessche C, Van Imschoot G, Van Wonterghem E, Vereecke L, Vandenbroucke RE (2023) Gut microbiota regulates blood-cerebrospinal fluid barrier function and abeta pathology. EMBO J 42:e111515. https://doi.org/10.15252/embj.2022111515
Zaugg K, Yao Y, Reilly PT, Kannan K, Kiarash R, Mason J, Huang P, Sawyer SK, Fuerth B, Faubert Bet al et al (2011) Carnitine palmitoyltransferase 1 C promotes cell survival and tumor growth under conditions of metabolic stress. Genes Dev 25:1041–1051. https://doi.org/10.1101/gad.1987211
Zhang C, Liu J, Liang Y, Wu R, Zhao Y, Hong X, Lin M, Yu H, Liu L, Levine AJ al (2013) Tumour-associated mutant p53 drives the Warburg effect. Nat Commun 4:2935. https://doi.org/10.1038/ncomms3935
Zhang Y, Guo S, Xie C, Fang J (2020) Uridine Metabolism and Its Role in Glucose, Lipid, and Amino Acid Homeostasis. Biomed Res Int 2020: 7091718 https://doi.org/10.1155/2020/7091718
Zhang T, Zhang Y, Liu J, Ma Y, Ye Q, Yan X, Ding L (2022) MicroRNA-377-3p inhibits hepatocellular carcinoma growth and metastasis through negative regulation of CPT1C-mediated fatty acid oxidation. Cancer Metab 10:2. https://doi.org/10.1186/s40170-021-00276-3
Zhou Y, Niu W, Luo Y, Li H, Xie Y, Wang H, Liu Y, Fan S, Li Z, Xiong W et al (2019) p53/Lactate dehydrogenase A axis negatively regulates aerobic glycolysis and tumor progression in breast cancer expressing wild-type p53. Cancer Sci 110: 939–949 https://doi.org/10.1111/cas.13928
Zhu G, Ying Y, Ji K, Duan X, Mai T, Kim J, Li Q, Yu L, Xu Y (2020) p53 coordinates glucose and choline metabolism during the mesendoderm differentiation of human embryonic stem cells. Stem Cell Res 49:102067. https://doi.org/10.1016/j.scr.2020.102067
Acknowledgements
Not applicable.
Funding
Research reported in this publication was supported by the National Cancer Institute of the National Institutes of Health under award number K08CA241651 (LYB). CBP is a McNair Scholar supported by the McNair Medical Institute at The Robert and Janice McNair Foundation. This work was partly supported by the generous philanthropic contributions to The University of Texas MD Anderson Cancer Center Glioblastoma Moon Shots Program™. The metabolomics core was supported by CPRIT Core Facility Support Award RP210227 (“Proteomic and Metabolomic Core Facility”), NCI Cancer Center Support Grant P30CA125123, NIH/NCI grant R01CA282282, NIH/NCI grant R01CA216426, and intramural funds from the Dan L. Duncan Cancer Center.
Author information
Authors and Affiliations
Contributions
DW and YF performed data analysis and interpretation and manuscript preparation. ARA contributed to data analysis and interpretation. BP, JTH and JJZ contributed to patient consent, sample annotation, collection, storage and distribution. AD, MS were responsible for clinical data collection and sample processing. NP performed metabolomic analysis. YE, FL, CBP and LYB are responsible for study conception and design. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Sample collection and processing and data collection and use were performed in accordance with local ethics regulations and approvals as well as the 1964 Helsinki declaration and its later amendments. All procedures were performed in accordance with local institutional review board guidelines. Patients consented to the use of their samples in this research project.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Supplementary Material 1:
Swimmer plot illustrating a simplified survival and event timeline of GBM patients. Events include surgery, chemotherapy, radiation treatment (RT), progression, last follow-up, death, and CSF collection. Some events were omitted or combined for clarity. All available dates for known events are included in Supplementary Material 4. Exact dates for administration of chemotherapy and radiation are not available for patients 29 and 37. (A) Swimmer plot of GBM patients who had CSF collected prior to treatment. (B) Swimmer plot of GBM patients who had CSF collected after initial treatment
Supplementary Material 2:
Volcano plots for analysis of metabolic differences based on sex in the GBM cohort. Colored points represent metabolites that are present at significantly different levels (-log10(p-value) > 1.3, log2(Group 1/Group 2) > + 1 or < -1). P-values and fold change calculations for all three figures are available in Supplementary Material 5. (A) Comparison of metabolites in the CSF of GBM patients between male and female cohorts. Six CSF metabolites are significantly different between male and female GBM patients. (B) Comparison of metabolites in the CSF of pre-treatment GBM patients between male and female cohorts. Four CSF metabolites are significantly different between male and female pre-treatment GBM patients. (C) Comparison of metabolites in the CSF of male GBM patients between the pre-treatment and post-treatment cohorts. 12 CSF metabolites are significantly different between pre-treatment and post-treatment male GBM patients
Supplementary Material 3:
Metabolic pathways of differentially abundant metabolites between pre-treatment GBM and control samples. (A) Carnitine shuttle pathway in fatty acid oxidation (FAO). (B) Shikimate pathway used by plants, bacteria, algae, and other microorganisms in the biosynthesis of aromatic compounds. (C) Catabolism of L-glutamine to B-alanine using uridine as an intermediate. (D) Butanoate metabolism with N-acetylputrescine and aminobutanal as intermediate compounds. (E) Cholesterol synthesis pathway with farnesyl diphosphate as an intermediate compound
Supplementary Material 4:
Data used in the study, including clinical data and metabolite levels for all patient samples; treatment and disease progression dates for the GBM cohort; survival data; mutation data for the most common mutated genes; and mutation frequencies for all mutated genes
Supplementary Material 5:
Significance values and fold change values for: metabolites that are significantly different and meet the fold change cutoff (+/-1) between pre-treatment and post-treatment GBM samples; all measured CSF metabolites between male and female patients in the GBM cohort; all measured CSF metabolites between male and female patients in the pre-treatment GBM cohort; and all measured CSF metabolites between pre-treatment and post-treatment samples in the GBM-male-only cohort
Supplementary Material 6:
Methods for metabolomic analysis of CSF samples
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Wang, D.H., Fujita, Y., Dono, A. et al. The genomic alterations in glioblastoma influence the levels of CSF metabolites. acta neuropathol commun 12, 13 (2024). https://doi.org/10.1186/s40478-024-01722-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40478-024-01722-1