The integrative metabolomic-transcriptomic landscape of glioblastome multiforme

Comprehensive understanding of glioblastoma molecular phenotypes: classification, characteristics, and transition

Among central nervous system-associated malignancies, glioblastoma (GBM) is the most common and has the highest mortality rate. The high heterogeneity of GBM cell types and the complex tumor microenvironment frequently lead to tumor recurrence and sudden relapse in patients treated with temozolomide. In precision medicine, research on GBM treatment is increasingly focusing on molecular subtyping to precisely characterize the cellular and molecular heterogeneity, as well as the refractory nature of GBM toward therapy. Deep understanding of the different molecular expression patterns of GBM subtypes is critical. Researchers have recently proposed tetra fractional or tripartite methods for detecting GBM molecular subtypes. The various molecular subtypes of GBM show significant differences in gene expression patterns and biological behaviors. These subtypes also exhibit high plasticity in their regulatory pathways, oncogene expression, tumor microenvironment alterations, and differential responses to standard therapy. Herein, we summarize the current molecular typing scheme of GBM and the major molecular/genetic characteristics of each subtype. Furthermore, we review the mesenchymal transition mechanisms of GBM under various regulators.

Alternative RNA splicing modulates ribosomal composition and determines the spatial phenotype of glioblastoma cells

Glioblastoma (GBM) is characterized by exceptionally high intratumoral heterogeneity. However, the molecular mechanisms underlying the origin of different GBM cell populations remain unclear. Here, we found that the compositions of ribosomes of GBM cells in the tumour core and edge differ due to alternative RNA splicing. The acidic pH in the core switches before messenger RNA splicing of the ribosomal gene RPL22L1 towards the RPL22L1b isoform. This allows cells to survive acidosis, increases stemness and correlates with worse patient outcome. Mechanistically, RPL22L1b promotes RNA splicing by interacting with lncMALAT1 in the nucleus and inducing its degradation. Contrarily, in the tumour edge region, RPL22L1a interacts with ribosomes in the cytoplasm and upregulates the translation of multiple messenger RNAs including TP53. We found that the RPL22L1 isoform switch is regulated by SRSF4 and identified a compound that inhibits this process and decreases tumour growth. These findings demonstrate how distinct GBM cell populations arise during tumour growth. Targeting this mechanism may decrease GBM heterogeneity and facilitate therapy.

Integrative Metabolomics Reveals Deep Tissue and Systemic Metabolic Remodeling in Glioblastoma

(1) Background: Glioblastoma is the most common malignant brain tumor in adults. Its etiology remains unknown in most cases. Glioblastoma pathogenesis consists of a progressive infiltration of the white matter by tumoral cells leading to progressive neurological deficit, epilepsy, and/or intracranial hypertension. The mean survival is between 15 to 17 months. Given this aggressive prognosis, there is an urgent need for a better understanding of the underlying mechanisms of glioblastoma to unveil new diagnostic strategies and therapeutic targets through a deeper understanding of its biology. (2) Methods: To systematically address this issue, we performed targeted and untargeted metabolomics-based investigations on both tissue and plasma samples from patients with glioblastoma. (3) Results: This study revealed 176 differentially expressed lipids and metabolites, 148 in plasma and 28 in tissue samples. Main biochemical classes include phospholipids, acylcarnitines, sphingomyelins, and triacylglycerols. Functional analyses revealed deep metabolic remodeling in glioblastoma lipids and energy substrates, which unveils the major role of lipids in tumor progression by modulating its own environment. (4) Conclusions: Overall, our study demonstrates in situ and systemic metabolic rewiring in glioblastoma that could shed light on its underlying biological plasticity and progression to inform diagnosis and/or therapeutic strategies.

Bioengineered Models to Study Microenvironmental Regulation of Glioblastoma Metabolism

Despite extensive research and aggressive therapies, glioblastoma (GBM) remains a central nervous system malignancy with poor prognosis. The varied histopathology of GBM suggests a landscape of differing microenvironments and clonal expansions, which may influence metabolism, driving tumor progression. Indeed, GBM metabolic plasticity in response to differing nutrient supply within these microenvironments has emerged as a key driver of aggressiveness. Additionally, emergent biophysical and biochemical interactions in the tumor microenvironment (TME) are offering new perspectives on GBM metabolism. Perivascular and hypoxic niches exert crucial roles in tumor maintenance and progression, facilitating metabolic relationships between stromal and tumor cells. Alterations in extracellular matrix and its biophysical characteristics, such as rigidity and topography, regulate GBM metabolism through mechanotransductive mechanisms. This review highlights insights gained from deployment of bioengineering models, including engineered cell culture and mathematical models, to study the microenvironmental regulation of GBM metabolism. Bioengineered approaches building upon histopathology measurements may uncover potential therapeutic strategies that target both TME-dependent mechanotransductive and biomolecular drivers of metabolism to tackle this challenging disease. Longer term, a concerted effort integrating in vitro and in silico models predictive of patient therapy response may offer a powerful advance toward tailoring of treatment to patient-specific GBM characteristics.

Network analysis with biological data of cancer patients: A scoping review

BACKGROUND & OBJECTIVE

Network Analysis (NA) is a mathematical method that allows exploring relations between units and representing them as a graph. Although NA was initially related to social sciences, the past two decades was introduced in Bioinformatics. The recent growth of the networks' use in biological data analysis reveals the need to further investigate this area. In this work, we attempt to identify the use of NA with biological data, and specifically: (a) what types of data are used and whether they are integrated or not, (b) what is the purpose of this analysis, predictive or descriptive, and (c) the outcome of such analyses, specifically in cancer diseases.

METHODS & MATERIALS

The literature review was conducted on two databases, PubMed & IEEE, and was restricted to journal articles of the last decade (January 2010 - December 2019). At a first level, all articles were screened by title and abstract, and at a second level the screening was conducted by reading the full text article, following the predefined inclusion & exclusion criteria leading to 131 articles of interest. A table was created with the information of interest and was used for the classification of the articles. The articles were initially classified to analysis studies and studies that propose a new algorithm or methodology. Each one of these categories was further screened by the following clustering criteria: (a) data used, (b) study purpose, (c) study outcome. Specifically for the studies proposing a new algorithm, the novelty presented in each one was detected.

RESULTS

& Conclusions: In the past five years researchers are focusing on creating new algorithms and methodologies to enhance this field. The articles' classification revealed that only 25% of the analyses are integrating multi-omics data, although 50% of the new algorithms developed follow this integrative direction. Moreover, only 20% of the analyses and 10% of the newly developed methodologies have a predictive purpose. Regarding the result of the works reviewed, 75% of the studies focus on identifying, prognostic or not, gene signatures. Concluding, this review revealed the need for deploying predictive and multi-omics integrative algorithms and methodologies that can be used to enhance cancer diagnosis, prognosis and treatment.

Targeting Immunometabolism in Glioblastoma

We have only recently begun to understand how cancer metabolism affects antitumor responses and immunotherapy outcomes. Certain immunometabolic targets have been actively pursued in other tumor types, however, glioblastoma research has been slow to exploit the therapeutic vulnerabilities of immunometabolism. In this review, we highlight the pathways that are most relevant to glioblastoma and focus on how these immunometabolic pathways influence tumor growth and immune suppression. We discuss hypoxia, glycolysis, tryptophan metabolism, arginine metabolism, 2-Hydroxyglutarate (2HG) metabolism, adenosine metabolism, and altered phospholipid metabolism, in order to provide an analysis and overview of the field of glioblastoma immunometabolism.

Radiogenomics of Gliomas

The 2016 World Health Organization brain tumor classification is based on genomic and molecular profile of tumor tissue. These characteristics have improved understanding of the brain tumor and played an important role in treatment planning and prognostication. There is an ongoing effort to develop noninvasive imaging techniques that provide insight into tissue characteristics at the cellular and molecular levels. This article focuses on the molecular characteristics of gliomas, transcriptomic subtypes, and radiogenomic studies using semantic and radiomic features. The limitations and future directions of radiogenomics as a standalone diagnostic tool also are discussed.

Integrated Metabolomics and Transcriptomics Using an Optimised Dual Extraction Process to Study Human Brain Cancer Cells and Tissues

The integration of untargeted metabolomics and transcriptomics from the same population of cells or tissue enhances the confidence in the identified metabolic pathways and understanding of the enzyme–metabolite relationship. Here, we optimised a simultaneous extraction method of metabolites/lipids and RNA from ependymoma cells (BXD-1425). Relative to established RNA (mirVana kit) or metabolite (sequential solvent addition and shaking) single extraction methods, four dual-extraction techniques were evaluated and compared (methanol:water:chloroform ratios): cryomill/mirVana (1:1:2); cryomill-wash/Econospin (5:1:2); rotation/phenol-chloroform (9:10:1); Sequential/mirVana (1:1:3). All methods extracted the same metabolites, yet rotation/phenol-chloroform did not extract lipids. Cryomill/mirVana and sequential/mirVana recovered the highest amounts of RNA, at 70 and 68% of that recovered with mirVana kit alone. sequential/mirVana, involving RNA extraction from the interphase of our established sequential solvent addition and shaking metabolomics-lipidomics extraction method, was the most efficient approach overall. Sequential/mirVana was applied to study a) the biological effect caused by acute serum starvation in BXD-1425 cells and b) primary ependymoma tumour tissue. We found (a) 64 differentially abundant metabolites and 28 differentially expressed metabolic genes, discovering four gene-metabolite interactions, and (b) all metabolites and 62% lipids were above the limit of detection, and RNA yield was sufficient for transcriptomics, in just 10 mg of tissue.

Perspective of mesenchymal transformation in glioblastoma

Despite aggressive multimodal treatment, glioblastoma (GBM), a grade IV primary brain tumor, still portends a poor prognosis with a median overall survival of 12–16 months. The complexity of GBM treatment mainly lies in the inter- and intra-tumoral heterogeneity, which largely contributes to the treatment-refractory and recurrent nature of GBM. By paving the road towards the development of personalized medicine for GBM patients, the cancer genome atlas classification scheme of GBM into distinct transcriptional subtypes has been considered an invaluable approach to overcoming this heterogeneity. Among the identified transcriptional subtypes, the mesenchymal subtype has been found associated with more aggressive, invasive, angiogenic, hypoxic, necrotic, inflammatory, and multitherapy-resistant features than other transcriptional subtypes. Accordingly, mesenchymal GBM patients were found to exhibit worse prognosis than other subtypes when patients with high transcriptional heterogeneity were excluded. Furthermore, identification of the master mesenchymal regulators and their downstream signaling pathways has not only increased our understanding of the complex regulatory transcriptional networks of mesenchymal GBM, but also has generated a list of potent inhibitors for clinical trials. Importantly, the mesenchymal transition of GBM has been found to be tightly associated with treatment-induced phenotypic changes in recurrence. Together, these findings indicate that elucidating the governing and plastic transcriptomic natures of mesenchymal GBM is critical in order to develop novel and selective therapeutic strategies that can improve both patient care and clinical outcomes. Thus, the focus of our review will be on the recent advances in the understanding of the transcriptome of mesenchymal GBM and discuss microenvironmental, metabolic, and treatment-related factors as critical components through which the mesenchymal signature may be acquired. We also take into consideration the transcriptomic plasticity of GBM to discuss the future perspectives in employing selective therapeutic strategies against mesenchymal GBM.

Metabolic alterations in meningioma reflect the clinical course

Background

Meningiomas are common brain tumours that are usually defined by benign clinical course. However, some meningiomas undergo a malignant transformation and recur within a short time period regardless of their World Health Organization (WHO) grade. The current study aimed to identify potential markers that can discriminate between benign and malignant meningioma courses.

Methods

We profiled the metabolites from 43 patients with low- and high-grade meningiomas. Tumour specimens were analyzed by nuclear magnetic resonance analysis; 270 metabolites were identified and clustered with the AutoPipe algorithm.

Results

We observed two distinct clusters marked by alterations in glycine/serine and choline/tryptophan metabolism. Glycine/serine cluster showed significantly lower WHO grades and proliferation rates. Also progression-free survival was significantly longer in the glycine/serine cluster.

Conclusion

Our findings suggest that alterations in glycine/serine metabolism are associated with lower proliferation and more recurrent tumours. Altered choline/tryptophan metabolism was associated with increases proliferation, and recurrence. Our results suggest that tumour malignancy can be reflected by metabolic alterations, which may support histological classifications to predict the clinical outcome of patients with meningiomas.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-07887-5.

Metabolomics and Multi-Omics Integration: A Survey of Computational Methods and Resources

As researchers are increasingly able to collect data on a large scale from multiple clinical and omics modalities, multi-omics integration is becoming a critical component of metabolomics research. This introduces a need for increased understanding by the metabolomics researcher of computational and statistical analysis methods relevant to multi-omics studies. In this review, we discuss common types of analyses performed in multi-omics studies and the computational and statistical methods that can be used for each type of analysis. We pinpoint the caveats and considerations for analysis methods, including required parameters, sample size and data distribution requirements, sources of a priori knowledge, and techniques for the evaluation of model accuracy. Finally, for the types of analyses discussed, we provide examples of the applications of corresponding methods to clinical and basic research. We intend that our review may be used as a guide for metabolomics researchers to choose effective techniques for multi-omics analyses relevant to their field of study.

Metabolic Reprogramming in Cancer Is Induced to Increase Proton Production

Considerable metabolic reprogramming has been observed in a conserved manner across multiple cancer types, but their true causes remain elusive. We present an analysis of around 50 such reprogrammed metabolisms (RM) including the Warburg effect, nucleotide de novo synthesis, and sialic acid biosynthesis in cancer. Analyses of the biochemical reactions conducted by these RMs, coupled with gene expression data of their catalyzing enzymes, in 7,011 tissues of 14 cancer types, revealed that all RMs produce more H⁺ than their original metabolisms. These data strongly support a model that these RMs are induced or selected to neutralize a persistent intracellular alkaline stress due to chronic inflammation and local iron overload. To sustain these RMs for survival, cells must find metabolic exits for the nonproton products of these RMs in a continuous manner, some of which pose major challenges, such as nucleotides and sialic acids, because they are electrically charged. This analysis strongly suggests that continuous cell division and other cancerous behaviors are ways for the affected cells to remove such products in a timely and sustained manner. As supporting evidence, this model can offer simple and natural explanations to a range of long-standing open questions in cancer research including the cause of the Warburg effect. SIGNIFICANCE: Inhibiting acidifying metabolic reprogramming could be a novel strategy for treating cancer.

Current and Future Trends on Diagnosis and Prognosis of Glioblastoma: From Molecular Biology to Proteomics

Glioblastoma multiforme is the most aggressive malignant tumor of the central nervous system. Due to the absence of effective pharmacological and surgical treatments, the identification of early diagnostic and prognostic biomarkers is of key importance to improve the survival rate of patients and to develop new personalized treatments. On these bases, the aim of this review article is to summarize the current knowledge regarding the application of molecular biology and proteomics techniques for the identification of novel biomarkers through the analysis of different biological samples obtained from glioblastoma patients, including DNA, microRNAs, proteins, small molecules, circulating tumor cells, extracellular vesicles, etc. Both benefits and pitfalls of molecular biology and proteomics analyses are discussed, including the different mass spectrometry-based analytical techniques, highlighting how these investigation strategies are powerful tools to study the biology of glioblastoma, as well as to develop advanced methods for the management of this pathology.

Integration of metabolomic and transcriptomic profiles to identify biomarkers in serum of lung cancer

We used blood serum samples collected from 31 lung cancer (LC) patients and 29 healthy volunteers in this study. Levels of serum metabolites were qualitative quantified with gas chromatography-mass spectrometry (GC-MS), and the data were analyzed by partial least-squares discrimination analysis (PLS-DA). Based on the Kyoto Encyclopedia of Genes and Genomes database, we performed pathway-based analysis utilizing metabolites presented at differential abundance between the LC serum samples and the normal healthy serum samples for systematical investigation on the metabolic alterations associated with LC pathogenesis. Finally, we analyzed the significantly enriched pathways as well as their relevant differentially expressed messenger RNAs, and drawn a correlation network plot to identify the serum metabolic biomarkers and the significantly altered metabolic pathways for LC. GC-MS analysis showed that 23 of the 169 metabolites identified were significantly different. PLS-DA model revealed that 13 of these metabolites were with variable importance > 1, and particularly five were with area under curve > 0.9. Pathway-based analysis demonstrated that five of eight enriched metabolic pathways were statistically significant with false discovery rate < 0.05. Lastly, the correlation networks between these pathways and their related genes suggested that 29 genes had correlation degree > 10, which were mainly engaged in the purine metabolism. In conclusion, we identified indole-3-lactate, erythritol, adenosine-5-phosphate, paracetamol and threitol as serum metabolic biomarkers for LC through metabolomics analysis. Besides, we identified the purine metabolism as the significantly altered metabolic pathway in LC with the help of transcriptomics analysis.

Diffusion Tensor Imaging Reveals Microstructural Heterogeneity of Normal-Appearing White Matter and Related Cognitive Dysfunction in Glioma Patients

Immunohistochemical data based on isocitrate–dehydrogenase (IDH) mutation status have redefined glioma as a whole-brain disease, while occult tumor cell invasion along white matter fibers is inapparent in conventional magnetic resonance imaging (MRI). The functional and prognostic impact of focal glioma may however relate to the extent of white matter involvement. We used diffusion tensor imaging (DTI) to investigate microstructural characteristics of whole-brain normal-appearing white matter (NAWM) in relation to cognitive functions as potential surrogates for occult white matter involvement in glioma. Twenty patients (12 IDH-mutated) and 20 individually matched controls were preoperatively examined using DTI combined with a standardized neuropsychological examination. Tumor lesions including perifocal edema were masked, and fractional anisotropy (FA) as well as mean, radial, and axial diffusivity (MD, RD, and AD, respectively) of the remaining whole-brain NAWM were determined by using Tract-Based Spatial Statistics and histogram analyses. The relationship between extratumoral white matter integrity and cognitive performance was examined using partial correlation analyses controlling for age, education, and lesion volumes. In patients, mean FA and AD were decreased as compared to controls, which agrees with the notion of microstructural impairment of NAWM in glioma patients. Patients performed worse in all cognitive domains tested, and higher anisotropy and lower MD and RD values of NAWM were associated with better cognitive performance. In additional analyses, IDH-mutated and IDH-wildtype patients were compared. Patients with IDH-mutation showed higher FA, but lower MD, AD, and RD values as compared to IDH-wildtype patients, suggesting a better preserved microstructural integrity of NAWM, which may relate to a less infiltrative nature of IDH-mutated gliomas. Diffusion-based phenotyping and monitoring microstructural integrity of extratumoral whole-brain NAWM may aid in estimating occult white matter involvement and should be considered as a complementary biomarker in glioma.

Comprehensive and Reproducible Untargeted Lipidomic Workflow Using LC-QTOF Validated for Human Plasma Analysis

The goal of this work was to develop a label-free, comprehensive, and reproducible high-resolution liquid chromatography-mass spectrometry (LC-MS)-based untargeted lipidomic workflow using a single instrument, which could be applied to biomarker discovery in both basic and clinical studies. For this, we have (i) optimized lipid extraction and elution to enhance coverage of polar and nonpolar lipids as well as resolution of their isomers, (ii) ensured MS signal reproducibility and linearity, and (iii) developed a bioinformatic pipeline to correct remaining biases. Workflow validation is reported for 48 replicates of a single human plasma sample: 1124 reproducible LC-MS signals were extracted (median signal intensity RSD = 10%), 50% of which are redundant due to adducts, dimers, in-source fragmentation, contaminations, or positive and negative ion duplicates. From the resulting 578 unique compounds, 428 lipids were identified by MS/MS, including acyl chain composition, of which 394 had RSD < 30% inside their linear intensity range, thereby enabling robust semiquantitation. MS signal intensity spanned 4 orders of magnitude, covering 16 lipid subclasses. Finally, the power of our workflow is illustrated by a proof-of-concept study in which 100 samples from healthy human subjects were analyzed and the data set was investigated using three different statistical testing strategies in order to compare their capacity in identifying the impact of sex and age on circulating lipids.

Gene Expression Analysis Reveals Novel Gene Signatures Between Young and Old Adults in Human Prefrontal Cortex

Human neurons function over an entire lifetime, yet the molecular mechanisms which perform their functions and protecting against neurodegenerative disease during aging are still elusive. Here, we conducted a systematic study on the human brain aging by using the weighted gene correlation network analysis (WGCNA) method to identify meaningful modules or representative biomarkers for human brain aging. Significantly, 19 distinct gene modules were detected based on the dataset GSE53890; among them, six modules related to the feature of brain aging were highly preserved in diverse independent datasets. Interestingly, network feature analysis confirmed that the blue modules demonstrated a remarkably correlation with human brain aging progress. Besides, the top hub genes including PPP3CB, CAMSAP1, ACTR3B, and GNG3 were identified and characterized by high connectivity, module membership, or gene significance in the blue module. Furthermore, these genes were validated in mice of different ages. Mechanically, the potential regulators of blue module were investigated. These findings highlight an important role of the blue module and its affiliated genes in the control of normal brain aging, which may lead to potential therapeutic interventions for brain aging by targeting the hub genes.

MR-spectroscopic imaging of glial tumors in the spotlight of the 2016 WHO classification

BACKGROUND

The purpose of this study is to map spatial metabolite differences across three molecular subgroups of glial tumors, defined by the IDH1/2 mutation and 1p19q-co-deletion, using magnetic resonance spectroscopy. This work reports a new MR spectroscopy based classification algorithm by applying a radiomics analytics pipeline.

MATERIALS

65 patients received anatomical and chemical shift imaging (5 × 5 × 20 mm voxel size). Tumor regions were segmented and registered to corresponding spectroscopic voxels. Spectroscopic features were computed (n = 860) in a radiomic approach and selected by a classification algorithm. Finally, a random forest machine-learning model was trained to predict the molecular subtypes.

RESULTS

A cluster analysis identified three robust spectroscopic clusters based on the mean silhouette widths. Molecular subgroups were significantly associated with the computed spectroscopic clusters (Fisher's Exact test p < 0.01). A machine-learning model was trained and validated by public available MRS data (n = 19). The analysis showed an accuracy rate in the Random Forest model by 93.8%.

CONCLUSIONS

MR spectroscopy is a robust tool for predicting the molecular subtype in gliomas and adds important diagnostic information to the preoperative diagnostic work-up of glial tumor patients. MR-spectroscopy could improve radiological diagnostics in the future and potentially influence clinical and surgical decisions to improve individual tumor treatment.

Microenvironment-Derived Regulation of HIF Signaling Drives Transcriptional Heterogeneity in Glioblastoma Multiforme

The evolving and highly heterogeneous nature of malignant brain tumors underlies their limited response to therapy and poor prognosis. In addition to genetic alterations, highly dynamic processes, such as transcriptional and metabolic reprogramming, play an important role in the development of tumor heterogeneity. The current study reports an adaptive mechanism in which the metabolic environment of malignant glioma drives transcriptional reprogramming. Multiregional analysis of a glioblastoma patient biopsy revealed a metabolic landscape marked by varying stages of hypoxia and creatine enrichment. Creatine treatment and metabolism was further shown to promote a synergistic effect through upregulation of the glycine cleavage system and chemical regulation of prolyl-hydroxylase domain. Consequently, creatine maintained a reduction of reactive oxygen species and change of the α-ketoglutarate/succinate ratio, leading to an inhibition of HIF signaling in primary tumor cell lines. These effects shifted the transcriptional pattern toward a proneural subtype and reduced the rate of cell migration and invasion in vitroImplications: Transcriptional subclasses of glioblastoma multiforme are heterogeneously distributed within the same tumor. This study uncovered a regulatory function of the tumor microenvironment by metabolism-driven transcriptional reprogramming in infiltrating glioma cells. Mol Cancer Res; 16(4); 655-68. ©2018 AACR.