Abstract
Objective
To investigate the performance of diffusion-tensor imaging (DTI) and hydrogen proton magnetic resonance spectroscopy (1H-MRS) parameters in predicting the immunohistochemistry (IHC) biomarkers of glioma.
Methods
Patients with glioma confirmed by pathology from March 2015 to September 2019 were analyzed, the preoperative DTI and 1H-MRS images were collected, apparent diffusion coefficient (ADC) and fractional anisotropy (FA), in the lesion area were measured, the relative values relative ADC (rADC) and relative FA (rFA) were obtained by the ratio of them in the lesion area to the contralateral normal area. The peak of each metabolite in the lesion area of 1H-MRS image: N-acetylaspartate (NAA), choline (Cho), and creatine (Cr), and metabolite ratio: NAA/Cho, NAA/(Cho + Cr) were selected and calculated. The preoperative IHC data were collected including CD34, Ki-67, p53, S-100, syn, vimentin, NeuN, Nestin, and glial fibrillary acidic protein.
Results
One predicting parameter of DTI was screened, the rADC of the Ki-67 positive group was lower than that of the negative group. Two parameters of 1H-MRS were found to have significant reference values for glioma grades, the NAA and Cr decreased as the grade of glioma increased, moreover, Ki-67 Li was negatively correlated with NAA and Cr.
Conclusion
NAA and Cr have potential application value in predicting glioma grades and tumor proliferation activity. Only rADC has predictive value for Ki-67 expression among DTI parameters.
Keywords
Introduction
Glioma is the most aggressive primary brain tumor with a poor prognosis and low survival. 1 According to the World Health Organization (WHO) grading system, it is divided into 4 pathological grades (I-IV), 2 namely, grades I and II are low-grade gliomas (LGGs), and grades III and IV are high-grade gliomas (HGGs). In 2016 and 2021, WHO updated the classification of glioma by combining molecular markers with histopathology, emphasizing the importance of molecular detection.2,3 The fifth edition of the WHO Classification of central nervous system tumors (WHO CNS 5), published in 2021, has achieved substantial changes by further promoting the role of molecular diagnosis in CNS tumor classification but still rooted in other established methods of tumor characterization, including histology and immunohistochemistry (IHC). 3 Among them, the expression of Ki-67, vimentin, CD34, S-100, p53, syn, and Glial Fibrillary Acidic Protein (GFAP) are important factors in judging and predicting the biological behavior of tumor cells.4–8 Recently, molecular biomarkers have become increasingly important in providing assistance and defining diagnostic information.
Diffusion-tensor imaging (DTI) is an in vivo functional imaging technology developed on the basis of diffusion-weighted imaging (DWI). It can comprehensively evaluate the diffusion movement of water molecules, the degree of compression, infiltration, and destruction of surrounding white matter fiber bundles from the microscopic perspective. Studies have shown that the apparent diffusion coefficient (ADC) and fractional anisotropy (FA) of the tumor have high application value in the differential diagnosis, grades, and the degree of white matter infiltration of glioma.9,10
According to WHO CNS 5, in addition to major changes in the classification of glioma, it is important to provide information on the role of Hydrogen proton magnetic resonance spectroscopy (1H-MRS) in glioma diagnosis, follow-up, and research. 11 1H-MRS is a proven technique for glioma diagnosis and follow-up in clinical practice. Tumor metabolite levels can be assessed, the main metabolites are N-acetylaspartate (NAA), choline (Cho), and creatine (Cr), NAA is a marker of neuronal activity, Cr is a cue for energy metabolism in brain tissue, and Cho is a marker of cell membrane synthesis. In glioma, the general rule is that NAA and Cr decrease with increasing tumor grade, while other metabolites increase. Cho/NAA quantitatively reflects tumor cell proliferation and local neuronal destruction without complex transitions. 12
At present, only DTI and MRS parameters predicting Ki-67 Labeling Index (Ki-67 Li) tumor proliferation activity studies can be detected, 13 while vimentin, CD34, S-100, p53, Syn, and GFAP have not been discussed yet. Therefore, this study was based on the predictive efficacy of 3 T 1H-MRS and DTI parameters on glioma grades and IHC characteristics, in order to preliminarily explore the feasibility of 1H-MRS and DTI parameters in predicting molecular typing of glioma.
Materials and Method
This retrospective study was according to the World Medical Association Declaration of Helsinki statement for research involving human subjects and was approved by the institutional ethics committee of the hospital (Research Ethics Committee of our Hospital. Approval No.2021-001-01) on March 8, 2021. All patients signed informed consent for magnetic resonance imaging (MRI) safety examination before the MRI examination. All data were fully anonymized before we accessed them.
Population
WHO CNS 5 has changed all CNS WHO tumor grades to Arabic numerals: (1-4). Because CNS tumor grades correlate with overall expected clinical biological behavior, the WHO CNS 5 generally retains the grade ranges used for tumor types in previous versions, so in our study, the glioma grades I to IV promulgated by WHO in 2016 were still used.
A total of 77 patients with pathologically confirmed gliomas from March 2015 to September 2019 in the PACS system were retrospectively included. Among them, 57 cases had complete DTI parameters, including 31 males and 26 females, with an average age of (49.70 ± 14.05) years from 7 to 72 years. In the 57 patients, there were 57 Ki-67 IHC cases, 48 CD34 IHC cases, 48 S-100 IHC cases, 48 vimentin IHC cases, 36 syn IHC cases, and 29 p53 IHC cases, respectively. 28 cases had the nestin IHC test and 57 cases had the GFAP IHC test, they were all positive expressions, so nestin and GFAP were not included in the study. A total of 31 patients were tested for NeuN IHC (label 0-29 cases and label 1-2 cases), since the total weight of each case was < 5, the significance could not be calculated, so NeuN was not included. All IHC indicators were divided into label 0 group (negative expression ones) and label 1 group (positive expression ones). More information is listed in Table 1.
Demographic and Clinical Characteristics of Patients.
Abbreviations: DTI, diffusion tensor imaging; MRS, magnetic resonance spectroscopy; WHO, World Health Organization.
Among the above 57 cases, 42 cases had complete MRS parameters, including 25 males and 17 females, with an average age of (47.45 ± 13.96) from 7 to 71 years. In the 42 patients, there were 42 Ki-67 IHC cases, 35 CD34 IHC cases, 36 S-100 IHC cases, 36 vimentin IHC cases, 29 syn IHC cases, and 20 p53 IHC cases, respectively. All 42 patients underwent the GFAP IHC test, and as all of them were positive, it was not analyzed. A total of 22 cases underwent the Nestin IHC test(label 0-2 cases and label 1-20 cases), and 23 patients were tested for NeuN IHC (label 0-20 cases and label 1-3 cases), since the total weight of Nestin and NeuN was < 5, the significance could not be calculated, so no analysis was made. More information is listed in Table 1.
Inclusion criteria: Patients with glioma were diagnosed by histopathologic examination through surgery. Exclusion criteria: Intracranial decompression, chemotherapy, and radiotherapy were performed before the MRI examination. Artifacts appeared in the MRI images.
MRI Protocol
Philips 3.0 T MRI system and 8-channel orthogonal head coil were used. SE-EPI sequence was performed for DTI, and the scanning parameters were as follows: TR: 10 000 ms, TE was the smallest. 1H-MRS acquisition parameters: Chemical shift imaging (CSI) point resolved spectrum (2d_press_144 fast sense), TR = 1500 ms TE = 144 ms Voxel size = 10 mm × 10 mm, FOV = 80 mm × 110 mm.
Immunohistochemical status
All patients in this study underwent surgical resection of glioma. IHC staining was performed on formalin-fixed and paraffin-embedded sections of brain tumor resected specimens. The expressions of Ki-67, vimentin, CD34, S-100, p53, syn, and GFAP in glioma were detected by SP staining with benchmark GX automatic IHC (Figure 1). The rate of Ki-67 Li < 20% was classified as negative and ≥ 20% as positive, 14 and the cell proliferation index was measured with IHC analysis by using the rabbit anti-human monoclonal antibody against the Ki-67 protein. IHC indicators were classified as negative and positive which were represented by 0, 1. IHC staining results were evaluated by 2 pathologists.

Immunohistochemistry of positive and negative staining (100 ×).
Data Labeling and Measurement
DTI parameter measurement (Figure 2): DTI used the continuous fiber tract tracing method. In the lesion area, the ADC, FA, and fiber bundle length were measured, and the relative values rADC, rFA, and rlength were obtained by the ratio of ADC, and FA values in the lesion area to the contralateral normal brain area.

Collection of diffusion tensor imaging (DTI) patients.
The lesion area of the MRS image was selected, and the metabolite peak and the metabolite ratio were calculated (Figure 3). Regions of interest (ROIs) were delineated on the solid parts of each tumor center, ROIs on necrotic areas, cystic intertumoral areas, and areas of presumed perilesional vasogenic edema were excluded from the analysis. We selected NAA, Cho, Cr, NAA/Cho, and NAA/(Cho + Cr).

Collection of magnetic resonance spectroscopy (MRS) patients.
Among all participants, a consensus was reached between a radiologist (with 10 years of experience in radiology) and a technologist (with 5 + years of experience in MRI procedures).
Statistical Analysis
IBM SPSS Statistics.22 was used for statistical analysis. It was used to test the normality (Kolmogorov Smirnov, P > .05) and the homogeneity of variance (Levene test, P > .05) of continuous variables in the clinic, 2-independent sample t-test (meeting the normality and homogeneity of variance test) was conducted to compare the distribution differences between the 2 groups, approximate t-test was applied when equal variances not assumed, rank sum test (Mann-Whitney nonparametric test) was performed when the quantitative data were skewed. Analysis of variance (F-test) was selected to compare the mean of 3 sample groups, and the Kruskal-Wallis Test (χ2) was used for unequal variances. The correlation between the Ki-67 Li and DTI, MRS parameters were evaluated with the Pearson correlation coefficient, and the Spearman correlation coefficient was selected as the data with bivariate non-normal distribution. P < .05 was statistically significant.
Results
The Relationship Between DTI Parameters and the Negative and Positive of IHC Indicators
The relationship between DTI-specific parameters ADC, AA, rADC, rFA, rlength, and Ki-67, CD34, S-100, vimentin, p53, and syn was selected, respectively, only the rADC of the Ki-67 positive group was lower than that in negative ones (Z = −1.947, P = .052), this indicates that increased cell density tends to occur in voxels with lower ADC values. The rlength in Ki-67, CD34, S-100, and p53 positive groups were lower than those in negative ones, it shows that the damage of fiber bundles in Ki-67, CD34, S-100, and p53 positive cases is slightly more serious than that in negative ones however, there were no statistical differences. Although the ADC in the p53 positive group was higher than that in negative ones, the difference was not statistically significant (Table 2).
The Relationship Between DTI Parameters and the Negative and Positive of IHC Indicators.
Abbreviations: IHC, immunohistochemistry; DTI, diffusion tensor imaging; ADC, apparent diffusion coefficient; rADC, relative apparent diffusion coefficient; rFA, relative fractional anisotropy; rlength, relative length.
The Relationship among DTI Parameters and Glioma Grades
There was no significant difference among DTI-specific parameters ADC, FA, rADC, rFA, rlength, and glioma grades (Table 3).
The Relationship Among DTI Parameters and Glioma Grade.
Abbreviations: WHO, World Health Organization; DTI, diffusion tensor imaging; ADC, apparent diffusion coefficient; rADC, relative apparent diffusion coefficient; rFA, relative fractional anisotropy; rlength, relative length.
Correlation Between Glioma Proliferation Index Ki67 Li and DTI Parameters
There was no correlation between Ki-67 Li and DTI parameters FA, ADC, rFA, rADC, but rADC decreased with Ki-67 Li increasing (r = −0.239, P = .073), and rFA increased while Ki-67 Li increasing (r = 0.230, P = .086) (Figure 3).
MRS Parameter Kurtosis in Differentiating Negative to Positive IHC Indicators
The MRS parameters NAA, Cho, Cr, NAA/Cho, and NAA/(Cho + Cr) were not related to Ki-67, CD34, S-100, vimentin, p53, and syn (Table 4).
MRS Parameter Kurtosis in Differentiating Negative to Positive of IHC Indicators.
Abbreviations: MRS, magnetic resonance spectroscopy; IHC, immunohistochemistry; Cr, creatine; NAA, N-acetylaspartate; Cho, choline.
The Relationship Between MRS Parameters and Glioma Grades
The difference in NAA among WHO II, III, and IV grades of glioma was statistically significant (F = 3.257, P = .049), and the NAA decreased with the increased glioma grades. The difference in Cr among WHO II, III, and IV grades of glioma was statistically significant (F = 4.948, P = 0.012), with the increased grades, the value of Cr decreased. However, there was no significant difference between Cho, NAA/Cho, NAA/(Cho + Cr), and glioma grades (Table 5).
The Relationship Between MRS Parameters and Glioma Grades.
Abbreviations: MRS, magnetic resonance spectroscopy; WHO, World Health Organization; Cr, creatine; NAA, N-acetylaspartate; Cho, choline.
We selected the average value of Cr in WHO III as the reference value and calculated the specificity and sensitivity in predicting the grade of glioma, they are 75.0% and 61.1%, respectively, Cr has a relatively high specificity and low sensitivity in distinguishing the glioma grades. The selection of the cut-off point of Cr was related to it (Table 6).
The Relationship Between MRS Parameter Cr and Glioma Grades.
Abbreviations: MRS, magnetic resonance spectroscopy; Cr, creatine; WHO, World Health Organization.
Correlation Between Ki-67 Li and MRS Parameters
Ki-67 Li was negatively correlated with NAA (r = −0.322, P = .037), and the NAA decreased with Ki-67 Li increasing. Ki-67 Li was negatively correlated with Cr (r = −0.305, P = .050), and the Cr decreased with the increased Ki-67 Li. There was no correlation between Ki-67 Li and Cho, NAA/Cho, and NAA/(Cho + Cr) (Figure 4). This suggests that the proposed parameters NAA, and Cr may provide an effective tool to identify tumors with high proliferative activity (Figure 5).

A scatter plot of the relationship between Ki-67 Li and DTI parameters.

A scatter plot of the relationship between Ki-67 Li and MRS parameters.
Discussion
At present, there are sufficient studies showing that gliomas with the same or similar histological features may carry different molecular or genetic information. 3 As a biomarker of cell proliferation, Ki-67 has been incorporated into the grade and prognosis prediction of CNS tumors. 15 CD34 is a transmembrane phosphoglycoprotein, first discovered in hematopoietic stem and progenitor cells, it is known as the best marker for microvessel density studies and plays a crucial role in the regulation of glioma angiogenesis, CD34 overexpression is associated with higher-grade gliomas and can serve as a potential glioma diagnostic and prognostic marker, as well as a useful therapeutic target. 16 Pleomorphic low-grade neuroepithelial tumor of the young (PLNTY) is a glioma associated with a history of epilepsy in youth, a diffuse growth pattern, frequent oligodendroglioma-like components, calcification, CD34 immunoreactivity, and the MAPK pathway, CD34 immunostaining was often intense and diffuse in the tumor. 3 Vimentin is an intermediate filament family protein that maintains cell integrity and participates in multiple cell signaling pathways to regulate cancer cell motility and invasion, and it is an independent important prognostic factor in patients with high-grade gliomas. 17 S-100 protein 18 named for its solubility in saturated ammonium sulfate, is derived from brain tissue, is a dimer, and belongs to the thermolabile acidic calcium-binding protein. It is present in gliomas, but its amount varies with the histological type of the tumor, mostly in ependymomas, it is also present in schwannomas and neurofibromas, but not in neuronal-derived in tumors such as medulloblastoma and neuroblastoma. A recent study showed 19 that p53-IHC should be used as a complement to the morphological diagnosis of H&E and that p53-IHC diffuse or strong positivity is sufficient to diagnose astrocytoma without 1p/19q testing, and the mutational status of p53 is correlated with IHC, there appears to be an incomplete parallel between the results. Syn is synaptophysin, mainly present in neuron neuroendocrine cells in presynaptic vesicles of neurons, and this antibody is mainly used to label neuroendocrine cells and their tumors. 20 GFAP is almost exclusively expressed in astrocytes, 21 and all cases in this study were positive for GFAP, so it could not be used for dichotomous studies. Therefore, fully considering the molecular pathological characteristics of glioma can help to determine the heterogeneity and prognosis of glioma. 22
In the research of Figini et al, 23 they pointed out that diffusion MRI has the potential to provide such in vivo indicators. It provides easy access to the microstructure of the whole tumor at a rather high spatial resolution. DTI is the most often utilized approach nowadays. 24 ADC and FA are DTI indices that correspond with changes in cellular density and extracellular matrix characteristics caused by glioma invasion and proliferation. FA was substantially greater in grade II and III isocitrate dehydrogenase (IDH) wild-type gliomas than in IDH mutant gliomas. 23 The present study demonstrates that the DTI parameter rADC has a predictive value for the negative and positive of IHC marker Ki-67, although the rlength value in Ki-67, CD34, S-100, and p53 positive groups was lower than in negative ones, the ABADC value of p53-positive cases was higher than that of negative ones, but unfortunately the difference was not statistically significant. The interpretation of DTI results in clinical work was divided into 2 aspects. On the one hand, the lesions were directly observed through the ADC map and FA map calculated by DTI. On the other hand, the lesions were quantitatively evaluated by measuring parameters such as ADC and FA values of different grades of glioma, and the overall ADC and FA values were obtained by the fiber tract tracing method, which reflected the destruction of the fiber bundles in the brain tissue by the tumor. But there is disagreement about the interpretation of ADC changes in DWI-derived brain tumors, it has been reported that restricted spread was affected not only by tumor cell density and metabolic activity, but also by other factors, such as ischemia or compression, in addition, DWI-derived ADCs were sensitive to other tissue characteristics, including due to vasogenic edema or extracellular fluid arising from the tumor-induced disruption of extracellular structures. 25 Due to the poor display of glioma lesions on the ADC and FA maps, our study did not choose to outline the lesion contour on the ADC and FA maps to extract image features. Instead, DTI used a continuous fiber tract tracing method, and quantitative indicators such as ADC, FA, and other values of the target lesion were measured. The results of this study showed that the measurement of DTI quantitative indicators only rADC could be used for the prediction of glioma IHC, this indicates that increased cell density tends to occur in voxels with lower ADC values. A recent retrospective DTI study has stratified LGG according to IDH and 1p/19q status. They reported that IDH wild-type astrocytoma and oligodendroglioma had lower ADC and higher FA, but there was no significant difference between glioma with and without 1p/19q deletion.26–28
In our study, DTI parameters failed to predict the grade of glioma. No correlation was found between Ki-67 Li and DTI parameters such as ABFA, ABADC, rFA, and rADC, but rADC decreased with the increased Ki-67 Li, and rFA increased with the increased Ki-67 Li. Similar studies suggested that some features of FA may be helpful in distinguishing LGG from HGG.29,30 There were also studies that showed that FA played no role in differentiating glioma grades. 31
Preoperative grade of glioma based on conventional MRI can be difficult, especially in the absence of significant edema and/or contrast enhancement. This is a common diagnostic challenge, even in high-grade tumors. In cases of uncertain diagnosis, more information can be obtained with MRS. 32 Currently, optimized 1H-MRS can noninvasively detect 2-hydroxyglutarate, a specific metabolite of IDH gene mutation. 12 Our results did not find that the MRS parameters NAA, Cho, Cr, NAA/Cho, NAA/(Cho + Cr) could differentiate Ki-67, CD34, S-100, vimentin, p53, and syn negative to positive. However, our study showed that NAA and Cr could provide a better differentiation between grades of glioma by the estimation of tumor cellularity. Correlation of Ki-67 Li with MRS parameters using Pearson's correlations in order to assess which parameters could predict tumor proliferative activity, Ki-67 Li was negatively correlated with NAA and Cr, this suggests that NAA, Cr may provide an effective tool to identify tumors with high proliferative activity. Our study demonstrated that NAA and Cr were promising imaging modalities for the noninvasive estimation of glioma grades and proliferative activity. Moreover, we selected the average value of Cr in WHO III as the reference value, Cr has a relatively high specificity in distinguishing the glioma grades. It has been reported that the most important metabolite that differentiated tumor progression (TP) from radiation-induced pseudoprogression (PSP) was Cho, however, Cr was not considered to be affected by radiation damage. Thus, in brain tissue undergoing radiation necrosis, an increased Cho/Cr ratio was observed. 33 It would be fascinating to investigate the potential of 1H-MRS in the preoperative assessment of several molecular markers, as well as the ability of this approach to predict progression. The approach might be useful in differentiating distinct subtypes of glioma regarding the 2021 WHO Classification of CNS tumors.
There are several potential limitations in this study. The main limitation of this study is the small sample size and not doing power analysis for sample size calculation. Secondly, since DTI and MRS characterize different tumor characteristics, the 2 imaging modalities may also provide complementary information in the pretreatment evaluation of gliomas, in the current study, we focused on evaluating IHC protein markers, and further research could assess genotyping to determine whether they would provide additional clinically useful information in glioma evaluation. In addition, not all patients were tested for all of the above IHC indicators. Therefore, future studies should validate the findings of this study in a larger cohort with detailed pathology. An innovation of our study is that it combined the quantitative parameters of DTI and MRS to identify or predict IHC indicators, and maximized the value of MR research for preoperative pathological grading of gliomas.
Conclusions
Overall, the present study demonstrates that the DTI parameter rADC has predictive value for the negative and positive of Ki-67 in glioma, and the MRS parameters NAA and Cr have potential application value for predicting the grade of glioma and tumor proliferative activity. However, there was no statistically significant difference between the other DTI parameters as well as all MRS parameters and the multiple IHC features of glioma. Meanwhile, DTI parameters could not predict the glioma grades in our study.
Supplemental Material
sj-docx-1-mix-10.1177_15353508241261583 - Supplemental material for Study on the Relationship Between MRI Functional Imaging and Multiple Immunohistochemical Features of Glioma: A Noninvasive and More Precise Glioma Management
Supplemental material, sj-docx-1-mix-10.1177_15353508241261583 for Study on the Relationship Between MRI Functional Imaging and Multiple Immunohistochemical Features of Glioma: A Noninvasive and More Precise Glioma Management by Jing Li, Jingtao Sun, Ning Wang and Yan Zhang in Molecular Imaging
Footnotes
Abbreviations
Acknowledgments
We appreciate Professor Huaijun Liu's guidance and help.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
Ethical Approval
The study was approved by the institutional ethics committee of the hospital (Research Ethics Committee of our Hospital. Approval No. 2021-001-01) on March 8, 2021. All patients signed informed consent for an MRI safety examination before the MRI examination. All data were fully anonymized before we accessed them.
Supplemental material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
