ORIGINAL RESEARCH article

Front. Oncol., 10 May 2019

Sec. Cancer Imaging and Image-directed Interventions

Volume 9 - 2019 | https://doi.org/10.3389/fonc.2019.00328

Assessment of Overall Survival in Glioma Patients as Predicted by Metabolomic Criteria

  • 1. Department of Neurosurgery, Hospital Universitario La Paz, Madrid, Spain

  • 2. Institute of Biomedical Research “Alberto Sols” CSIC/UAM, Madrid, Spain

  • 3. Department of Statistics CSIC, Madrid, Spain

  • 4. Department of Neurosurgery, University Hospital Reina Sofía, Córdoba, Spain

Abstract

Objective: We assess the efficacy of the metabolomic profile from glioma biopsies in providing estimates of postsurgical Overall Survival in glioma patients.

Methods: Tumor biopsies from 46 patients bearing gliomas, obtained neurosurgically in the period 1992–1998, were analyzed by high resolution 1H magnetic resonance spectroscopy (HR- 1H MRS), following retrospectively individual postsurgical Overall Survival up to 720 weeks.

Results: The Overall Survival profile could be resolved in three groups; Short (shorter than 52 weeks, n = 19), Intermediate (between 53 and 364 weeks, n = 19) or Long (longer than 365 weeks, n = 8), respectively. Classical histopathological analysis assigned WHO grades II–IV to every biopsy but notably, some patients with low grade glioma depicted unexpectedly Short Overall Survival, while some patients with high grade glioma, presented unpredictably Long Overall Survival. To explore the reasons underlying these different responses, we analyzed HR-1H MRS spectra from acid extracts of the same biopsies, to characterize the metabolite patterns associated to OS predictions. Poor prognosis was found in biopsies with higher contents of alanine, acetate, glutamate, total choline, phosphorylcholine, and glycine, while more favorable prognosis was achieved in biopsies with larger contents of total creatine, glycerol-phosphorylcholine, and myo-inositol. We then implemented a multivariate analysis to identify hierarchically the influence of metabolomic biomarkers on OS predictions, using a Classification Regression Tree (CRT) approach. The CRT based in metabolomic biomarkers grew up to three branches and split into eight nodes, predicting correctly the outcome of 94.7% of the patients in the Short Overall Survival group, 78.9% of the patients in the Intermediate Overall Survival group, and 75% of the patients in the Long Overall Survival group, respectively.

Conclusion: Present results indicate that metabolic profiling by HR-1H MRS improves the Overall Survival predictions derived exclusively from classical histopathological gradings, thus favoring more precise therapeutic decisions.

Introduction

Gliomas are the most frequent primary brain tumors, currently managed through surgical resection, radiotherapy, and chemotherapy (1) approaches, but leading inevitably to large disability and mortality outcomes. The selection of the recommended therapeutic intervention in each case relies in estimates overall survival (OS) based commonly in histopathological and genetic criteria. However, current assessments of OS entail considerable uncertainties, limiting concomitantly more precise, effective, and personalized therapies. On these grounds, exploring additional criteria to improve OS predictions acquires vital relevance to improve treatment outcomes in glioma patients.

Histopathological and immunohistochemical criteria have classically provided the basis for the initial WHO classification of gliomas in grades I-IV (2, 3) determining, in general terms, the OS estimate and the recommended therapeutic intervention. More recently, the 2016 WHO classification of central nervous system tumors added an important collection of molecular signatures, restructuring the original histopathological classification of gliomas to include subgroups with specific genetic profiles (4). Although these refinements considerably improved the precision in the treatment prescribed, as well as our knowledge of glioma physiopathology and classification, the limited reproducibility of histopathological evaluations lead, not unfrequently, to imprecise histopathological classification and unreliable OS predictions at the individual level (57).

Magnetic Resonance Imaging approaches have been currently used to assess OS of gliomas. Briefly, radiomic parameters including surface area (8), shape features (9), tumor, and necrosis volumes, necrosis-tumor ratio (10) have been used to evaluate OS. However, these studies became many times limited to short OS periods as they evaluated only glioblastoma multiforme cases.

The metabolomic profiles of gliomas are able to provide an additional source of information to improve OS predictions, evaluating the down-stream metabolic alterations of aberrant cellularity and gene expression (11). Magnetic Resonance Spectroscopy (MRS) has been shown to be well-endowed to provide the metabolic profile of gliomas both in vivo and in vitro (12, 13). Briefly, in vivo1H MRS revealed non-invasively, important hallmarks of cancer, including alterations in pH homeostasis (14), energy related (15), and phospholipid metabolites (16, 17), an ensemble of valuable metabolic fingerprints to classify high grade (HGG) or low grade gliomas (LGG) (16, 18). Moreover, the metabolic profiles determined by 1H MRS in vivo reached considerable clinical prognostic relevance (1922). Alternatively, complementary in vitro HR-1H MRS approaches have proved to be able to resolve a larger number of metabolites than in vivo1H MRS, thus increasing the size of the metabolome investigated, the number of potential alterations detected, and their influence on the tumoral phenotype, at the expense of the more invasive in vitro methodology (23, 24). However, the predictive role of the metabolomic profiles obtained by HR-1H MRS in providing OS estimates, received considerably less attention.

On these grounds, we aimed here to provide a pilot study evaluating OS estimates derived from metabolomic biomarkers as detected by HR-1H MRS, using a retrospective database of human glioma biopsies.

Materials and Methods

Glioma Patients and Tumor Biopsies

This study was approved by the Ethics Committee of Clinical Research from the University Hospital La Paz (http://www.madrid.org/cs/Satellite?language=es&pagename=HospitalLaPaz/Page/HPAZ_home) and carried out following their recommendations. All subjects gave written informed consent in accordance with the Declaration of Helsinki. We retrospectively reviewed the database of the Neurosurgery Department of the University Hospital La Paz, selecting 66 consecutive patients with glioma (grades I-IV) who underwent neurosurgery during the period 1992–1998 (Figure 1). Briefly, solid parts of glioma tumors were extracted from the brain without the use of bipolar coagulation and divided into two adjacent and similar portions, one of them used for HR-1H MRS analysis and the other for histopathological diagnosis, following available WHO criteria (25). Tumor characteristics were evaluated by two independent radiologists and classified according to size, localization, and eloquence (26). Patients lost in follow up (n = 10) and those undergoing surgery for recurrence (n = 7) were excluded from further analysis. Grade I gliomas (n = 3) were also excluded because of their well-known physiopathological differences with the other glioma grades (27, 28). We then recovered the individual demographic, clinical, histopathological and in vitro spectroscopic 1H NMR features and gathered the OS information on postsurgical outcomes of these patients, including relevant clinical symptoms, adjuvant therapies, and OS time. No missing data for the variables of interest were found in these patients.

Figure 1

Histopathological Criteria

Histopathological grading of the biopsies was provided by the Anatomopathology Department of the Hospital, following standard WHO criteria (25), and archived until used in this study.

HR-1H MRS

The 1H MRS biopsy was immediately frozen in liquid nitrogen (−169°C) in the operating room of the Hospital and stored at −82°C until transferred in a liquid nitrogen container to the Institute of Biomedical Research CSIC/UAM for further processing and NMR analysis. Briefly, biopsies were reduced to powder in a previously chilled (methanol/dry ice) mortar, extracted with 6% perchloric acid, neutralized with KOD, lyophilized and resuspended in D2O (99.9% D) for 1H MRS analysis (12). HR-1H MRS spectra of biopsy extracts were acquired at 8.4 Tesla (360.13 MHz, pH 7.2, 22°C) in a Bruker AM-360 spectrometer equipped with a commercial 1H selective probe using 5-mm tubes and 0.5 ml of tissue extract. Acquisition conditions were: 90° pulses, 16.9 s total cycle time and 16,384 data points acquired in the time domain during of 1.901 s. The intensity of the residual water resonance was further reduced using a 5 s presaturating pulse centered on the water frequency. Prior to Fourier transformation, the free induction decays were zero-filled up to 128 K and multiplied by an exponential function resulting in 0.5 Hz artificial line broadening in the transformed spectrum. Further spectral processing, including phase and baseline corrections were performed by the same operator. Chemical shifts were referred to the methyl signal of TSP (2,2′-3,3′ tetradeutero trimethyl-sylil propionate sodium salt) at 0 ppm as an internal reference.

The following metabolites (resonances used in quantification, number of protons originating the resonance, multiplicity) could be consistently identified in the high resolution proton spectrum (29): Valine (1.09 ppm, 3H, d), Lactate (1.35 ppm, 3H, d), Alanine (1.45 ppm, 3H, d), Acetate (1.93 ppm, 3H, s), N-acetyl-aspartic acid (2.01 ppm, 3H, s), Gamma-amino butyric acid (2.31 ppm, 2H, t), Glutamate (2.43 ppm,2H, m), Glutamine (2.45 ppm, 2H, m), Aspartic acid (2.80 ppm, 2H, dd), Creatine (3.05 ppm, 3H, s) and Phosphocreatine (3.055 ppm, 3H, s), free Choline (3.20 ppm, 9H, s), Phosphorylcholine (3.22ppm, 9H, s), Glycerophophorylcholine (3.25 ppm, 9H, s), Taurine (3.45 ppm, 2H, t), Glycine (3.55 ppm, 2H, s), and Myo-inositol 4.07 ppm, 1H, dd). For every one of these metabolites, lorentzian curves were fitted to the most conveniently resolved proton resonances, and the resulting integral divided by the total number of protons of the corresponding metabolite (6). These values were further standardized by the sum of all the measured metabolites in the HR-1H NMR spectra and expressed as a molar percentage (18, 30). Assignments were performed with the aid of chemical shift values reported in the literature (29, 31) and confirmed when necessary by the addition of authentic standards.

Statistical Methods

Statistical analyses were performed using the IBM SPSS Statistics 24 package as implemented on an Intel-PC platform, operating under Windows 10 environment. Univariate statistical approaches provided means and standard errors for the molar fractions of every metabolite. To investigate statistical dependences between clinical features and groups of OS, we used asymptotic chi-square with Monte Carlo exact probability tests, and to test the differences of means within each metabolic variable through the OS groups, we used the ANOVA test and Student t-tests. Finally, to explore the hierarchical contribution of individual HR-1H MRS biomarkers to the three groups of OS, we implemented a multivariate Classification Regression Tree (CRT) (32), classifying automatically the database using hierarchical nodes and branches, selecting step-wise the optimal discriminant biomarker for each split from the collection of available HR-1H MRS variables. The dependent variable was OS, using Chi-squared Automatic Interaction Detection (CHAID) as a growing method to provide automatically the optimal splits in every branch. Finally, we used these results to generate a Classification-Confusion Matrix (CCM), summarizing the correct and incorrect classifications provided by the metabolomics CRT, yielding the global percentage of correct classifications. Statistical significance in the ANOVA and multivariate analysis was defined as p = 0.05, considering confidence intervals higher than 95%.

Results

Glioma Database

We investigated 46 patients (23 males, 23 females) with a median age of 49 years, presenting the following glioma grade distribution; Grade II (11 cases), Grade III (16 cases), and Grade IV (19 cases). The database of glioma patients (Figure 2) showed two clearly separated groups by OS, either depicting a Long Overall Survival (Long-OS) that survive more than 364 weeks (w), or less. The latter group, including patients with a wide range of survival (1–364 w), was further divided in two groups using the median of survival as a cut-point, resulting in Short Overall Survival (Short-OS) patients (1–52 w), or Intermediate Overall Survival (Intermediate-OS) patients (53–364 w), respectively.

Figure 2

Demographic, clinical, and radiological variables, therapeutic regimes and histopathological profiles of these groups are summarized in Table 1. A Chi-square Monte Carlo test was performed to compare the differences in independent clinical variables among the three survival categories. Radiotherapy, histopathologic grade, age, and localization showed statistically significant influence on OS.

Table 1

FeatureOverall survivalTotal patients n (% column)Chi-square (montecarlo sig.)
Short n (% row)Intermediate n (% row)Long n (% row)
Age<25 y0 (0.0%)5 (62.5%)3(37.5%)8 (17.4%)0.01
25–54 y5 (31.3%)7 (43.8%)4(25%)16 (34.8%)
>54 y14 (63.6%)7 (31.8%)1 (4.5%)22 (47.8%)
SexFemale8 (34.8%)11(47.8%)4 (17.4%)23 (50%)0.74
Men11(47.8%)8 (34.8%)4 (17.4%)23 (50%)
ComorbilityaYes10 (34.5%)13 (44.8%)6 (20.7%)29 (63%)0.56
No9 (52.9%)6 (35.3%)2 (11.8%)17 (37%)
LocalizationbA2 (15.4%)7 (53.8%)4 (30.8%)13 (30.2%)0.04
B5 (41.7%)6 (50%)1 (8.3%)12 (27.9%)
C12 (66.7%)3 (16.7%)3 (16.7%)18 (41.9%)
Tumor volumebSmall13 (39.4%)15 (45.5%)5 (15.2%)33 (73.3%)0.49
Big6 (50%)3 (25%)3 (25%)12 (26.7%)
EloquencybYes16 (53.3%)10 (33.3%)4 (13.3%)30 (69.8%)0.19
No3 (23.1%)6 (46.2%)4(30.8 %)13 (30.2%)
ResectionPartial/Complete10 (52.6%)5 (26.3%)4 (21.1%)19 (41.3%)0.27
Biopsy9 (33.3%)14 (51.9%)4 (14.8%)27 (58.7%)
RadiotherapyYes19 (54.3%)12 (34.3%)4 (11.4%)35 (76.1%)0.003
No0 (0%)7 (63.6%)4 (36.4%)11 (23.9%)
ChemotherapyYes0 (0%)1 (33.3%)2 (66.7%)3 (6.5%)0.07
No19(44.2%)18 (41.9%)6 (14%)43 (93.5%)
Histopathology gradecII0 (0%)6 (54.5%)5 (45.5%)11 (23.9%)<0.001
III3 (18.8%)10 (62.5%)3 (18.8%)16 (34.8%)
IV16(84.2%)3 (15.8%)0 (0%)19 (41.3%)

Clinical features and overall survival of glioma patients.

a

Arterial hypertension, diabetes mellitus and/or pulmonary, renal, cardiac, oncologic or any severe disease,

b

Classification of tumors according to (26).

c

Histopathologic grade according to (3).

HR-1H MRS

An illustrative example of underdetermined histopathological OS prediction is provided in Figure 3, showing representative HR-1H MRS spectra from extracts of glioma biopsies obtained from two young male patients, assigned the same histopathological Grade II, but resulting in very different OS. Despite both patients underwent complete surgical resection without adjuvant radio- or chemotherapy, the patient represented in Figure 3A, survived <3 years (149 w), while the patient represented by Figure 3B survived more than 11 years (605 w).

Figure 3

Interestingly, HR-1H MRS spectra of these biopsies disclosed remarkable differences, particularly the relative increases in myo-inositol (MI) and glycerol-phosphorylcholine (GPC) in Figure 3B. These findings suggested that HR-1H MRS analysis of biopsy extracts could contribute additional OS criteria to those normally obtained from general histopathological classification, thus prompting further HR-1H MRS analyses of the database.

Univariate Statistics

OS was well-reflected in the metabolic profiles obtained, with evident relationships between OS and specific metabolite changes (Figure 4). Table 2 provides a more detailed analysis of the relationship between OS and specific metabolite molar fractions, highlighting the discriminant power in OS of each metabolite (F-value) and its statistical significance (p-value) as derived from ANOVA tests.

Figure 4

Table 2

MetaboliteShort (n = 19)Intermediate (n = 19)Long (n = 8)TotalFbp
OVERALL SURVIVALa
MI9.03 ± 1.2614.75 ± 2.1425.72 ± 1.5114.30 ± 1.3615.10<0.001
Ala8.47 ± 1.125.67 ± 0.951.67 ± 0.246.13 ± 0.7013.65<0.001
GPC3.17 ± 0.345.17 ± 0.727.18 ± 0.974.70 ± 0.429.44<0.001
Gly13.38 ± 1.896.95 ± 1.033.61 ± 0.689.03 ± 1.059.37<0.001
PC4.17 ± 0.482.84 ± 0.481.43 ± 0.283.14 ± 0.328.360.001
tCr13.36 ± 1.1015.03 ± 1.1920.84 ± 1.6315.35 ± 0.816.530.003
Ac4.49 ± 0.973.09 ± 0.810.94 ± 0.173.30 ± 0.555.720.006
fCho1.33 ± 0.130.86 ± 0.170.73 ± 0.241.03 ± 0.104.650.01
Glu12.20 ± 1.2712.09 ± 1.566.04 ± 1.1611.08 ± 0.913.620.03
Gln15.47 ± 1.1918.76 ± 1.3818.51 ± 1.3017.36 ± 0.812.020.14
Val2.34 ± 0.471.60± 0.431.45 ± 0.231.88 ± 0.270.8960.42
Asp0.77 ± 0.140.70 ± 0.140.50 ± 0.120.69 ± 0.090.5840.56
NAA3.35 ± 0.944.57 ± 0.954.55 ± 0.944.06 ± 0.570.5440.58
GABA1.25 ± 0.310.92 ± 0.151.17 ± 0.331.10 ± 0.150.4830.62
Tau5.69 ± 0.796.17 ± 0.714.95 ± 0.635.76 ± 0.450.4160.66
tCho8.67 ± 0.578.88 ± 0.769.33 ± 1.028.87 ± 0.420.1450.86

General linear model analysis of OS in glioma patients as revealed by the metabolic profile determined by HR-1H MRS.

a

Results are given as mean ± standard error of mean.

b

A Box Cox transform was performed on the data before running ANOVA analysis. Ac, Acetate; Ala, Alanine; Asp, aspartic acid; fCho, free choline; GABA, gamma-aminobutyric acid; Gln, glutamine; Glu, glutamate; Gly, glycine; GPC, glycerophosphocholine; MI, myo-inositol; NAA, N-acetyl-aspartic acid; PC, phosphorylcholine; Tau, taurine; tCr, total creatine; Val, valine. Bold characters indicate p < 0.05.

Briefly, OS increased with increasing levels of MI, GPC, and total Creatine (tCr), and decreased with increasing levels of alanine (Ala), free choline (fCho), glutamate (Glu), phosphorylcholine (PC), and glycine (Gly). The remaining metabolites detected by HR-1H NMR, including acetate (Ac), glutamine (Gln), valine (Val), aspartate (Asp), N-acetyl-aspartic acid (NAA), Gamma-amino butyric (GABA), taurine (Tau), and total choline (tCho) were not found to influence significantly OS. The most powerful discriminators of OS were (Metabolite/ F value/ p-value) in decreasing order; MI/15.1/0.000, Ala/13.6/0.000, GPC/9.5/0.001, Gly/9.4/0.000, tCr/6.5/0.003, Ac/5.7/0.006, fCho/4.6/0.015, Glu/3.6/0.035.

We also investigated the relationship between histopathological grade and the molar fractions of metabolites detectable in extracts glioma biopsies using ANOVA tests (Table 3). Notably, the priority of metabolites providing optimal glioma grade discrimination, was different from the one yielding optimal OS discriminative power (Table 2). The following molar ratios of metabolites were found to provide optimal discriminant power between histopathological grades (metabolite/F-value/p-value); MI/21.7/0.000, PC/11.7/0.000, GPC/9.2/0.000, Ala/8.8/0.001, Gly/6.5/0.003, tCr/6.3/0.004, Glu 5.9/0.005 and Gln/3.5/0.040.

Table 3

IVIIIIITotalFp
n = 11n = 16n = 19
HISTOPATHOLOGICAL GRADEa
MI8.41 ± 1.09b13.94 ± 1.9524.99 ± 2.314.3 ± 1.3621.71c<0.000
PC4.36 ± 0.562.90 ± 0.361.40 ± 0.253.14 ± 0.3211.74<0.000
GPC3.01 ± 0.195.32 ± 0.846.70 ± 0.874.70 ± 0.429.16<0.000
Ala8.36 ± 1.155.77 ± 1.092.82 ± 0.696.13 ± 0.78.850.001
Gly12.74 ± 1.97.98 ± 1.424.14 ± 0.539.03 ± 1.056.540.003
tCr13.54 ± 1.1714.35 ± 1.3119.94 ± 1.2815.35 ± 0.816.300.004
Glu13.85 ± 1.5510.89 ± 1.326.59 ± 0.9511.08 ± 0.915.900.005
Succ1.36 ± 0.491.11 ± 0.20.58 ± 0.081.09 ± 0.223.990.03
Gln14.95 ± 1.1419.05 ± 1.3119.05 ± 1.6717.36 ± 0.813.470.04
fCho1.23 ± 0.131.05 ± 0.210.67 ± 0.161.03 ± 0.13.100.05
Val2.32 ± 0.481.97 ± 0.490.98 ± 0.111.88 ± 0.272.510.09
Asp0.88 ± 0.160.63 ± 0.120.47 ± 0.120.69 ± 0.091.930.16
Ac3.75 ± 0.963.68 ± 0.981.95 ± 0.673.30 ± 0.551.490.24
Tau6.09 ± 0.856.05 ± 0.744.77 ± 0.465.76 ± 0.450.420.66
tCho8.60 ± 0.579.27 ± 0.868.76 ± 0.838.87 ± 0.420.240.79
GABA1.03 ± 0.191.25 ± 0.361.02 ± 0.231.10 ± 0.150.230.80
NAA4.14 ± 1.174.06 ± 0.863.94 ± 0.494.06 ± 0.570.010.99

General linear model analysis of the HR- 1H NMR metabolic profiles associated to different glioma histopathological grades.

a

According to Louis et al. (3)

b

Results are given as mean ± standard error of mean.

c

A Box Cox transform was performed on the data before running ANOVA analysis. Bold characters indicate p < 0.05.

Together, these results show that relevant metabolites contribute with different strengths either to histopathological grading or to OS predictions.

Classification Regression Trees (CRT)

To investigate the hierarchical contribution of these metabolites to the OS observed, we implemented a multivariate CRT (Figure 5) (32). Starting with the complete patient database (Node 0), high MI levels (Branch 1) provided the most powerful biomarker to predict Long-OS survival within the three OS groups. MI levels ≤ 23.35 were found in 37 biopsies (Node 1). Of these, only one patient survived more than 7 years (Long-OS), while the rest of the patients depicted either Intermediate- (17 patients) or Short-OS (19 patients). In contrast, MI levels >23.35 were detected in nine biopsies (Node 2), of which seven (77.8%) depicted Long-OS, two showed Intermediate-OS, and none had a Short-OS, suggesting that high MI levels dismiss a Short-OS prediction.

Figure 5

Patients with MI levels ≤ 23.35 could be further split in two groups using GPC (Brach 2). Sixty-nine percent of the patients with GPC levels ≤ 4.36 (Node 3) were classified as Short-OS, while 91% of the patients with GPC > 4.36 (Node 4) were classified as Intermediate-OS. GPC provided thus a convenient biomarker to distinguish between Short- and Intermediate-OS estimates.

The metabolomic CRT grew beyond Branch 2, improving the classification using either Gly or Ala splits (Branch 3). All patients of Node 3 with Gly levels higher 13.87 (Node 6), depicted Short-OS, suggesting that high Gly levels are predictive of a negative outcome. Patients from Node 4, could be further stratified by their Ala levels. Those having Ala levels lower than 8.48 (Node 7) indicated Intermediate-OS (100%). Summarizing, the metabolomic CRT indicated dominant roles of MI, GPC, Ala, and Gly OS prediction of glioma patients.

Finally, we compared the number of correct predictions derived from the metabolomic CRT with those observed clinically, in the Classification-Confusion Matrix of Table 4. Out of 19 patients in the Short-OS group, the metabolomic approach classified correctly 18 patients. In the Intermediate-OS group, the metabolomic approach classified correctly 15 of 19 patients. Finally, out of the 8 patients identified with Long-OS, the metabolomic approach correctly classified 6 patients.

Table 4

Observed overall survival (n = 46)Predicted overall survivalPercent correct classifications
ShortIntermediateLong
Short (n = 19)18a1094.7
Intermediate (n = 19)415078.9
Long (n = 8)11675.0
All patients (n = 46)84.7

Classification confusion matrix of correct/incorrect classifications of overall survival in patients bearing gliomas using metabolomic criteria.

a

Numbers in bold indicate number and percentages of correct classifications. Growing Methods: CRT, dependent variable OS.

In summary, the metabolomic classification reached defined OS predictions in all three groups, separating well the longer OS groups (Intermediate-OS and Long-OS). This entails considerable relevance, since the prediction of Long-OS in glioma patients remains currently a vital challenge for neurosurgeons, with important implications in the definition of the recommended therapeutic strategy.

Discussion

Previous OS Studies

The present study complements and extends earlier OS predictors based on the WHO classification (24) or in vivo MRI/MRS studies (1921, 3340), contributing a novel array of metabolomic biomarkers organized hierarchically by a multivariate CRT. Earlier studies implementing in vivo MRS/MRI approaches investigated mainly glioblastoma patients (Table 5), associating Short-OS to increases in PC (4) and tCho/NAA (20, 3640, 42), and longer OS to higher contents in GPC (35). Additionally, Cho/Cr ratio has been proposed as a biomarker of cellular proliferation and prognosis (22). More recently, a correlation between 2-hydroxiglutarate and the IDH1 mutation (43), suggested that in vivo detection of 2-hydroxiglutarate could become a useful prognostic biomarker. However, routinely and regular detection in vivo of 2-hydroxyglutarate still remains an important technological challenge in most imaging centers, limiting wider applications. In summary, the present study contributes a larger cohort of patients examined by routinely available in vitro MRS, followed during a longer period of time, including also a collection of both Low Grade Gliomas (LGG) and High Grade Gliomas (HGG). Notably, some of the metabolomic biomarkers found valuable here in OS prediction, like MI or GPC (Table 2), are difficult to resolve, or not even detectable in vivo. Consequently, our results suggest that postsurgical HR-1H NMR analysis of extracted tumor biopsies may provide a useful complement to available in vivo MRS explorations when addressing OS predictions.

Table 5

Authors & yearPatientsGlioma gradeFollow-up perioda, b
Li et al. (39)72HGGc17.2 ma
Reijneveld et al. (40)14LGGd30 m (9–40)b
Hattingen et al. (41)45LGG37 m (52.1–260.5)
Chang et al. (38)143LGG & HGGn.s.
Yamasaki et al. (33)HGG26.1 m (6.5–83.8)
Steffen-Smith et al. (37)39HGG7.1 m (1.6–61.6)
Quon et al. (36)26HGG22.9 m (5–37)
Hattingen et al. (35)32HGG (recidives)8.1 m
Tolia et al. (21)12HGGn.s.
Steidl et al. (34)37HGG (recidives)n.s.
Roldán et al. (20)28HGG3–98 m
Present serie46LGG & HGG14.9 m (0.24–170.4)

Overview of literature correlating OS and MRS biomarkers.

a

Median of the duration of the study in months,

b

Parenthesis includes the range of the study in months (m),

c

HGG: High grade glioma,

d

LGG: Low grade glioma, n.s.: not specified.

Metabolomic CRT

We implemented a CRT methodology to find, hierarchically, the best classification of OS estimates, using metabolomic criteria. These results complement the earlier decision tree of Li et al. (39), who considered age, MRI features (T2, T1, Volume Contrast Enhancement) and in vivo spectroscopic biomarkers (Lip+Lac, Cho/Cr, and Cr/NAA), as the main determinants predicting glioblastoma OS.

Our metabolomic CRT grew up to three branches. In the first branch, MI became the most robust biomarker of survival, with larger MI contents revealing longer survivals. Increased MI levels have been reported in inflammatory diseases as Alzheimer (41), renal failure (44), diabetes mellitus (44), and traumatic brain injury (45, 46), suggesting a universal role of this osmolite in pathophysiolological volume regulation. However, MI levels were reported previously to increase (4749) or decrease (44) with increasing glioma grade in vivo. Present results reveal that higher MI levels are associated to longer OS, and lower tumoral grades (Table 3). Since MI occurs primarily in normal astrocytes (44), we hypothesize that the MI resonances detected in the tumor biopsies reveal the healthy astrocyte content within the tumor mass. Larger relative MI contents reveal indirectly relatively larger normal astrocyte populations and smaller tumoral cell burdens, supporting consequently longer OS (41, 50, 51).

The second hierarchical branch classifying the lower MI content group, is GPC, with higher GPC levels predicting longer OS. We, and others, have previously reported that relatively higher contributions of GPC and PC are associated to low and high grade gliomas, respectively (16, 52). Interestingly, GPC and PC levels are thought to reflect the balance between phospholipid degradation and phospholipid synthesis, respectively, with increased GPC levels revealing relatively a negative balance between synthesis and degradation, lower tumoral proliferation, and more prolonged survivals.

Gly and Ala provided the third branch of OS discrimination for patient groups with low or high levels of GPC, respectively. High Gly and Ala levels revealed poor prognosis, associated to shorter OS. Indeed, Ala and Gly levels previously reported as hypoxia and redox stress biomarkers (53), revealing tumoral progression to hypoxia, redox stress, and fatal energy failure. They can now be associated to shorter OS predictions.

Limitations

The time span of biopsy collection in this study preceded some of the advances in the characterization of glioma genetics and their influence in malignancy and OS. This circumstance, and the long survival period investigated, precluded the use of genetic biomarkers validated later in the coverage of the present retrospective study. Thus, the correlations between metabolomic and genomic biomarkers of OS in glioma deserve further investigation. Finally, the number of patients involved in the present pilot study is admittedly small but sufficiently robust to support the use of metabolomic biomarkers detected by in vitro HR-1HMRS in OS predictions of postsurgical survival from glioma patients. A multicenter study to extend the number of patients and hospitals involved, is currently being implemented.

Conclusion

We used a multivariate CRT to assess postsurgical OS predictions based in the 1H HR-MRS analysis of the metabolomic profile from neurosurgical biopsies of glioma patients. Present results show that the metabolic profiles of glioma biopsies constitute accurate and independent biomarkers of OS in glioma patients.

Statements

Data availability statement

All datasets generated for this study are included in the manuscript.

Ethics statement

This study was carried out in accordance with the recommendations of the Ethics Committee of the Hospital La Paz with the approval number PI-2097 with written informed consent from all subjects. All subjects gave written informed consent in accordance with the Declaration of Helsinki.

Author contributions

MG-G collected and integrated retrospective patient data, investigated survival patterns, and wrote the first draft. SC analyzed HR-1H NMR spectra. LB provided the univariate and multivariate statistical analyses. PL-L acquired HR-1H NMR spectra. PF and AP validated demographic, radiological, and histopathological assessments. JR performed many of the neurosurgical procedures providing integrated clinical information, and JS conceived the study and wrote the final draft with all authors commenting.

Funding

This work was supported in part by grants PI2017/00361 from Instituto de Investigación Carlos III to JR, grant B2017/BMD-3688 from the Community of Madrid to JR and SC, and grant PI-0143-2016 from the Regional Ministry of Health of the Regional Government of Andalucía to JS.

Acknowledgments

Authors are grateful to Mr. Javier Pérez CSIC for the professional drafting of the illustrations, to Mrs. Maria José Guillén CSIC for skillful processing of the biopsies and to Mrs. Teresa Navarro CSIC for granting access to the HR-1HMRS facility of the Institute of Biomedical Research Alberto Sols.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

    Abbreviations

  • Ac

    Acetate

  • Ala

    Alanine

  • Asp

    aspartic acid

  • Cr

    creatine

  • CCM

    classification confusion matrix

  • CRT

    Classification Regression Tree

  • d

    doublet multiplicity

  • dd

    doublet of doublets multiplicity

  • fCho

    free choline

  • GABA

    gamma-aminobutyric acid

  • GLM

    General Lineal Model

  • Gln

    glutamine

  • Glu

    glutamate

  • Gly

    glycine

  • GPC

    glycerophosphocholine

  • GroPEtn

    glycerophosphoethanolamine

  • HGG, high grade glioma 1H MRS

    proton magnetic resonance spectroscopy

  • Lac

    lactate

  • LGG

    low grade glioma

  • Lip

    lipid

  • m

    multiplet

  • MRS

    magnetic resonance spectroscopy

  • MI

    myo-inositol

  • NAA

    N-acetyl-aspartic acid

  • OS

    overall survival

  • PC

    phosphorylcholine

  • PCr

    phosphocreatine

  • PDE

    phosphodiester

  • PE

    phosphoethanolamine

  • PtdCho

    phosphatidylcholine

  • PtdEtn

    phosphatidylethanolamine

  • PtdIns

    phosphatidylinositol

  • s, singlet multiplicity SEM

    Standard error of mean

  • t

    triplet

  • Tau

    taurine

  • tCho

    total choline

  • tCr

    total creatine

  • Val

    valine

  • w

    weeks.

References

  • 1.

    BushNAChangSMBergerMS. Current and future strategies for treatment of glioma. Neurosurg Rev. (2017) 40:114. 10.1007/s10143-016-0709-8

  • 2.

    BiernatW. 2000 World Health Organization classification of tumors of the nervous system. Pol J Pathol. (2000) 51:107114.

  • 3.

    LouisDNOhgakiHWiestlerODCaveneeWKBurgerPCJouvetAet al. The 2007 WHO classification of tumours of the central nervous system. Acta Neuropathol. (2007) 114:97109. 10.1007/s00401-007-0243-4

  • 4.

    LouisDNPerryAReifenbergerGvon DeimlingAFigarella-BrangerDCaveneeWKet al. The 2016 World Health Organization classification of tumors of the central nervous system: a summary. Acta Neuropathol. (2016) 131:80320. 10.1007/s00401-016-1545-1

  • 5.

    MittlerMAWaltersBCStopaEG. Observer reliability in histological grading of astrocytoma stereotactic biopsies. J Neurosurg. (1996) 85:10914. 10.3171/jns.1996.85.6.1091

  • 6.

    PraysonRAAgamanolisDPCohenMLEstesMLKleinschmidt-DeMastersBKAbdul-KarimFet al. Interobserver reproducibility among neuropathologists and surgical pathologists in fibrillary astrocytoma grading. J Neurol Sci. (2000) 175:339. 10.1016/S0022-510X(00)00274-4

  • 7.

    CastilloMSDavisFGSurawiczTBrunerJMBignerSCoonsSet al. Consistency of primary brain tumor diagnoses and codes in cancer surveillance systems. Neuroepidemiology. (2004) 23:8593. 10.1159/000073980

  • 8.

    CuiYThaKKTerasakaSYamaguchiSWangJKudoKet al. Prognostic imaging biomarkers in glioblastoma: development and independent validation on the basis of multiregion and quantitative analysis of MR images. Radiology. (2016) 278:54653. 10.1148/radiol.2015150358

  • 9.

    SanghaniPAngBTKingNKKRenH. Overall survival prediction in glioblastoma multiforme patients from volumetric, shape and texture features using machine learning. Surg Oncol. (2018) 27:70914. 10.1016/j.suronc.2018.09.002

  • 10.

    HenkerCKriesenTGlassÄSchneiderBPiekJ. Volumetric quantification of glioblastoma: experiences with different measurement techniques and impact on survival. J Neurooncol. (2017) 135:391402. 10.1007/s11060-017-2587-5

  • 11.

    Cuperlovic-CulfMFergusonDCulfAMorinPTouaibiaM. 1H NMR metabolomics analysis of glioblastoma subtypes: correlation between metabolomics and gene expression characteristics. J Biol Chem. (2012) 287:2016475. 10.1074/jbc.M111.337196

  • 12.

    RodaJMPascualJMCarcellerFGonzález-LlanosFPérez-HiguerasASoliveraJet al. Nonhistological diagnosis of human cerebral tumors by 1H magnetic resonance spectroscopy and amino acid analysis. Clin Cancer Res. (2000) 6:398393.

  • 13.

    OpstadKSWrightAJBellBAGriffithsJRHoweFA. Correlations between in vivo (1)H MRS and ex vivo (1)H HRMAS metabolite measurements in adult human gliomas. J Magn Reson Imaging. (2010) 31:28997. 10.1002/jmri.22039

  • 14.

    García-MartínMLHérigaultGRémyCFarionRBallesterosPColesJAet al. Mapping extracellular pH in rat brain gliomas in vivo by 1H magnetic resonance spectroscopic imaging: comparison with maps of metabolites. Cancer Res. (2001) 61:652431.

  • 15.

    Vander HeidenMGCantleyLCThompsonCB. Understanding the Warburg effect: the metabolic requirements of cell proliferation. Science. (2009) 324:102933. 10.1126/science.1160809

  • 16.

    RighiVRodaJMPazJMucciATugnoliVRodriguez-TarduchyGet al. 1H HR-MAS and genomic analysis of human tumor biopsies discriminate between high and low grade astrocytomas. NMR Biomed. (2009) 22:62937. 10.1002/nbm.1377

  • 17.

    GlundeKBhujwallaZMRonenSM. Choline metabolism in malignant transformation. Nat Rev Cancer. (2011) 11:83548. 10.1038/nrc3162

  • 18.

    SoliveraJCerdánSPascualJMBarriosLRodaJM. Assessment of 31P-NMR analysis of phospholipid profiles for potential differential diagnosis of human cerebral tumors. NMR Biomed. (2009) 22:66374. 10.1002/nbm.1387

  • 19.

    HattingenERaabPFranzKLanfermannHSetzerMGerlachRet al. Prognostic value of choline and creatine in WHO grade II gliomas. Neuroradiology. (2008) 50:75967. 10.1007/s00234-008-0409-3

  • 20.

    Roldan-ValadezERiosCMotola-KubaDMatus-SantosJVillaARMoreno-JimenezS. Choline-to-N-acetyl aspartate and lipids-lactate-to-creatine ratios together with age assemble a significant Cox's proportional-hazards regression model for prediction of survival in high-grade gliomas. Br J Radiol. (2016) 89:20150502. 10.1259/bjr.20150502

  • 21.

    ToliaMVerganelakisDTsoukalasNKyrgiasGPapathanasiouMMosaEet al. Prognostic value of MRS metabolites in postoperative irradiated high grade gliomas. Biomed Res Int. (2015) 2015:341042. 10.1155/2015/341042

  • 22.

    GaoWWangXLiFShiWLiHZengQ. Cho/Cr ratio at MR spectroscopy as a biomarker for cellular proliferation activity and prognosis in glioma: correlation with the expression of minichromosome maintenance protein 2. Acta Radiol. (2019) 60:10612. 10.1177/0284185118770899

  • 23.

    ShaoWGuJHuangCLiuDHuangHHuangZet al. Malignancy-associated metabolic profiling of human glioma cell lines using 1H NMR spectroscopy. Mol Cancer. (2014) 13:197. 10.1186/1476-4598-13-197

  • 24.

    GuidoniLRicci-VitianiLRosiAPalmaAGrandeSLucianiAMet al. 1H NMR detects different metabolic profiles in glioblastoma stem-like cells. NMR Biomed. (2014) 27:12945. 10.1002/nbm.3044

  • 25.

    KleihuesPBurgerPCScheithauerBW. The new WHO classification of brain tumours. Brain Pathol. (1993) 3:25568. 10.1111/j.1750-3639.1993.tb00752.x

  • 26.

    ShinodaJSakaiNMuraseSYanoHMatsuhisaTFunakoshiT. Selection of eligible patients with supratentorial glioblastoma multiforme for gross total resection. J Neurooncol. (2001) 52:16171. 10.1023/A:1010624504311

  • 27.

    MatyjaEGrajkowskaWStȩpieńKNaganskaE. Heterogeneity of histopathological presentation of pilocytic astrocytoma - diagnostic pitfalls. A review. Folia Neuropathol. (2016) 54:197211. 10.5114/fn.2016.62530

  • 28.

    CollinsVPJonesDTGianniniC. Pilocytic astrocytoma: pathology, molecular mechanisms and markers. Acta Neuropathol. (2015) 129:77588. 10.1007/s00401-015-1410-7

  • 29.

    CerdànSParrillaRSantoroJRicoM. 1H NMR detection of cerebral myo-inositol. FEBS Lett. (1985) 187:16772. 10.1016/0014-5793(85)81235-7

  • 30.

    KlunkWEXuCJPanchalingamKMcClureRJPettegrewJW. Analysis of magnetic resonance spectra by mole percent: comparison to absolute units. Neurobiol Aging. (1994) 15:13340. 10.1016/0197-4580(94)90153-8

  • 31.

    GovindarajuVYoungKMaudsleyAA. Proton NMR chemical shifts and coupling constants for brain metabolites. NMR Biomed. (2000) 13:12953. 10.1002/1099-1492(200005)13:3<129::AID-NBM619>3.0.CO;2-V

  • 32.

    BreimanL. Classification and Regression Trees. Belmont, CA: Wadsworth International Group (1984).

  • 33.

    YamasakiFKurisuKKajiwaraYWatanabeYTakayasuTAkiyamaYet al. Magnetic resonance spectroscopic detection of lactate is predictive of a poor prognosis in patients with diffuse intrinsic pontine glioma. Neuro Oncol. (2011) 13:791801. 10.1093/neuonc/nor038

  • 34.

    SteidlEPilatusUHattingenESteinbachJPZanellaFRonellenfitschMWet al. Myoinositol as a biomarker in recurrent glioblastoma treated with bevacizumab: A 1H-magnetic resonance spectroscopy study. PLoS ONE. (2016) 11:e0168113. 10.1371/journal.pone.0168113

  • 35.

    HattingenEBährORiegerJBlaselSSteinbachJPilatusU. Phospholipid metabolites in recurrent glioblastoma: in vivo markers detect different tumor phenotypes before and under antiangiogenic therapy. PLoS ONE. (2013) 8:e56439. 10.1371/journal.pone.0056439

  • 36.

    QuonHBrunetBAlexanderAMurthaAAbdulkarimBFultonDet al. Changes in serial magnetic resonance spectroscopy predict outcome in high-grade glioma during and after postoperative radiotherapy. Anticancer Res. (2011) 31:355965.

  • 37.

    Steffen-SmithEAShihJHHippSJBentRWarrenKE. Proton magnetic resonance spectroscopy predicts survival in children with diffuse intrinsic pontine glioma. J Neurooncol. (2011) 105:36573. 10.1007/s11060-011-0601-x

  • 38.

    ChangSMNelsonSVandenbergSChaSPradosMButowskiNet al. Integration of preoperative anatomic and metabolic physiologic imaging of newly diagnosed glioma. J Neurooncol. (2009) 92:40115. 10.1007/s11060-009-9845-0

  • 39.

    LiXJinHLuYOhJChangSNelsonSJ. Identification of MRI and 1H MRSI parameters that may predict survival for patients with malignant gliomas. NMR Biomed. (2004) 17:1020. 10.1002/nbm.858

  • 40.

    ReijneveldJCvan der GrondJRamosLMBrombergJETaphoornMJ. Proton MRS imaging in the follow-up of patients with suspected low-grade gliomas. Neuroradiology. (2005) 47:88791. 10.1007/s00234-005-1435-z

  • 41.

    HattingenERaabPFranzKZanellaFELanfermannHPilatusU. Myo-inositol: a marker of reactive astrogliosis in glial tumors?NMR Biomed. (2008) 21:23341. 10.1002/nbm.1186

  • 42.

    JaskólskiDJFortuniakJMajosAGajewiczWPapierzWLiberskiPPet al. Magnetic resonance spectroscopy in intracranial tumours of glial origin. Neurol Neurochir Pol. (2013) 47:43849. 10.5114/ninp.2013.32999

  • 43.

    LeatherTJenkinsonMDDasKPoptaniH. Magnetic resonance spectroscopy for detection of 2-hydroxyglutarate as a biomarker for IDH mutation in gliomas. Metabolites. (2017) 7:E29. 10.3390/metabo7020029

  • 44.

    CastilloMSmithJKKwockL. Correlation of myo-inositol levels and grading of cerebral astrocytomas. AJNR Am J Neuroradiol. (2000) 21:16459.

  • 45.

    PascualJMSoliveraJPrietoRBarriosLLópez-LarrubiaPCerdánSet al. Time course of early metabolic changes following diffuse traumatic brain injury in rats as detected by (1)H NMR spectroscopy. J Neurotrauma. (2007) 24:94459. 10.1089/neu.2006.0190

  • 46.

    CroallISmithFEBlamireAM. Magnetic resonance spectroscopy for traumatic brain injury. Top Magn Reson Imaging. (2015) 24:26774. 10.1097/RMR.0000000000000063

  • 47.

    KimJHChangKHNaDGSongICKwonBJHanMHet al. 3T 1H-MR spectroscopy in grading of cerebral gliomas: comparison of short and intermediate echo time sequences. AJNR Am J Neuroradiol. (2006) 27:14128.

  • 48.

    NatsumedaMIgarashiHNomuraTOguraRTsukamotoYKobayashiTet al. Accumulation of 2-hydroxyglutarate in gliomas correlates with survival: a study by 3.0-tesla magnetic resonance spectroscopy. Acta Neuropathol Commun. (2014) 2:158. 10.1186/s40478-014-0158-y

  • 49.

    FanGSunBWuZGuoQGuoY. In vivo single-voxel proton MR spectroscopy in the differentiation of high-grade gliomas and solitary metastases. Clin Radiol. (2004) 59:7785. 10.1016/j.crad.2003.08.006

  • 50.

    GalanaudDNicoliFChinotOConfort-GounySFigarella-BrangerDRochePet al. Noninvasive diagnostic assessment of brain tumors using combined in vivo MR imaging and spectroscopy. Magn Reson Med. (2006) 55:123645. 10.1002/mrm.20886

  • 51.

    OpstadKSLadroueCBellBAGriffithsJRHoweFA. Linear discriminant analysis of brain tumour (1)H MR spectra: a comparison of classification using whole spectra versus metabolite quantification. NMR Biomed. (2007) 20:76370. 10.1002/nbm.1147

  • 52.

    SabatierJGilardVMalet-MartinoMRanjevaJPTerralCBreilSet al. Characterization of choline compounds with in vitro 1H magnetic resonance spectroscopy for the discrimination of primary brain tumors. Invest Radiol. (1999) 34:2305. 10.1097/00004424-199903000-00013

  • 53.

    TsunZYPossematoR. Amino acid management in cancer. Semin Cell Dev Biol. (2015) 43:2232. 10.1016/j.semcdb.2015.08.002

Summary

Keywords

classification decision tree, glioma, metabolomic profile, high resolution proton magnetic resonance spectroscopy, overall survival

Citation

Gandía-González ML, Cerdán S, Barrios L, López-Larrubia P, Feijoó PG, Palpan Jr. A, Roda JM and Solivera J (2019) Assessment of Overall Survival in Glioma Patients as Predicted by Metabolomic Criteria. Front. Oncol. 9:328. doi: 10.3389/fonc.2019.00328

Received

13 February 2019

Accepted

11 April 2019

Published

10 May 2019

Volume

9 - 2019

Edited by

Bo Gao, Affiliated Hospital of Guizhou Medical University, China

Reviewed by

Bihong T. Chen, City of Hope National Medical Center, United States; Ru Jia Wang, Tangshan Gongren Hospital, China

Updates

Copyright

*Correspondence: María L. Gandía-González

This article was submitted to Cancer Imaging and Image-directed Interventions, a section of the journal Frontiers in Oncology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics