Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 04 January 2023
Sec. Surgical Oncology
This article is part of the Research Topic Advances in Precision Oncology and Multimodality Approach for the Management of Central Nervous System Tumors View all 4 articles

MRI radiomic features of peritumoral edema may predict the recurrence sites of glioblastoma multiforme

Hao Long,&#x;Hao Long1,2†Ping Zhang&#x;Ping Zhang3†Yuewei Bi,&#x;Yuewei Bi2,4†Chen Yang,Chen Yang2,4Manfeng Wu,Manfeng Wu2,4Dian He,Dian He2,4Shaozhuo Huang,Shaozhuo Huang2,4Kaijun Yang,Kaijun Yang1,2Songtao Qi,Songtao Qi1,2Jun Wang,,*&#x;Jun Wang1,2,4*‡
  • 1Department of Neurosurgery, Nanfang Hospital, Southern Medical University, Guangzhou, China
  • 2The First Clinical Medicine College, Southern Medical University, Guangzhou, China
  • 3Department of oncology, Guangdong 999 Brain Hospital, Guangzhou, China
  • 4Neural Networks Surgery Team, Southern Medical University, Guangzhou, China

Background and purpose: As one of the most aggressive malignant tumor in the central nervous system, the main cause of poor outcome of glioblastoma (GBM) is recurrence, a non-invasive method which can predict the area of recurrence pre-operation is necessary.To investigate whether there is radiological heterogeneity within peritumoral edema and identify the reproducible radiomic features predictive of the sites of recurrence of glioblastoma(GBM), which may be of value to optimize patients’ management.

Materials and methods: The clinical information and MR images (contrast-enhanced T1 weighted and FLAIR sequences) of 22 patients who have been histologically proven glioblastoma, were retrospectively evaluated. Kaplan-Meier methods was used for survival analysis. Oedematous regions were manually segmented by an expert into recurrence region, non-recurrence region. A set of 94 radiomic features were obtained from each region using the function of analyzing MR image of 3D slicer. Paired t test was performed to identify the features existing significant difference. Subsequently, the data of two patients from TCGA database was used to evaluate whether these features have clinical value.

Results: Ten features with significant differences between the recurrence and non-recurrence subregions were identified and verified on two individual patients from the TCGA database with pathologically confirmed diagnosis of GBM.

Conclusions: Our results suggested that heterogeneity does exist in peritumoral edema, indicating that the radiomic features of peritumoral edema from routine MR images can be utilized to predict the sites of GBM recurrence. Our findings may further guide the surgical treatment strategy for GBM.

Highlights

Peritumoral edema has heterogeneity.

The heterogeneity of peritumoral edema is related to the recurrence of glioblastoma.

Radiomics may predict the potential site of glioblastoma recurrence.

Introduction

Glioblastoma (GBM) is the most aggressive malignant tumor in the central nervous system. In order to avoid recurrence after surgery, maximal safe resection is always performed, the effect of which is affirmative. A recent study indicated that wider surgical resection was associated with better prognosis (1). Additionally, a meta-analysis suggested that gross total resection was superior to subtotal resection in improving the patients’ prognosis (2). However, excessive resection may lead to complications. Excessive removal of functional brain areas will significantly decrease the patients’ quality of life. Accordingly, it is necessary to find a new method to optimize surgical resection.

Notably, recurrences, as the main cause of poor outcome despite the complex surgical resection and adjuvant radiochemotherapy (3), mainly occur in the area of peritumoral edema (4). Thus, in order to improve the prognosis of this aggressive tumor, it is essential to optimize surgical resection and remove the peritumoral edema area with high risk of recurrence. In reality, GBM heterogeneity, which has been the hot topic in neuro-oncology lately, is not limited to the intratumoral region but also involves the peritumoral edema (5), and this heterogeneity can be used to identify the potential recurrence region. Many studies have investigated the cellular and molecular characteristics of the peritumoral edema via genomic, transcriptomic, or proteomic analyses (68). However, considering the low sensitivity, histopathological examination is not appropriate for identifying potential recurrence regions. Furthermore, the characteristics of a sample cannot properly represent the peritumoral edema as a whole.

Radiomics, a recent method that aims at extracting features from medical images using data-characterization algorithms, not only has the potential to reflect the underlying characteristics of peritumoral edema, but is also easier to perform (9). Compared with traditional radiographic images, digital images can be converted into quantitative radiomic image features via high-throughput computing. These radiomic features can partially reflect the biological features and underlying pathophysiology of the relevant tissues. After decades of development, radiomics has been applied in predicting the treatment response and outcomes of cancer, tumor staging, identifying tissue, and assessment of cancer genetics (10), including GBM (11) and grade II and III gliomas (12). Radiomics can also be used for tumor regional heterogeneity analysis (13). However, the association of GBM recurrence and certain radiomic features has rarely been investigated. Moreover, it has been shown that preoperative MRI parameters contribute to prediction of the prognosis of patients with GBM (14).

Therefore, this study aimed to identify heterogeneity within the peritumoral edema region via radiomics, which may be beneficial to the optimization of surgical resection.

Materials and methods

Study population

Ethics committee approval and patient consent were obtained. We retrospectively reviewed the data of patients diagnosed with GBM from April 2008 to October 2016 in our department. Patients were included if they met all of the following criteria: (1) histopathologically confirmed diagnosis of GBM based on the World Health Organization histological grading system by at least one expert pathologist; (2) available preoperative images before the initial surgery, including contrast-enhanced T1-weighted and three-dimensional (3D) fluid-attenuated inversion recovery (FLAIR) images; (3) no residual tumor detected after the operation; (4)pathologically and image confirmed tumor recurrence,and (5) available images of tumor recurrence before reoperation (including contrast-enhanced T1-weighted and FLAIR images). The follow-up data of two patients from The Cancer Genome Atlas (TCGA) database were used to test our hypothesis.

Endpoint

We considered OS and progression-free survival (PFS) as the endpoints of our study. OS was defined as time from treatment start to death or time from treatment start to the end of our study. PFS was defined as the time from treatment start to progression or death or time from treatment start to the end of our study, whichever occurred first.

MRI acquisition and segmentation

All MR images of all eligible patients were obtained using the same MRI scanner (Discovery ST, GE Healthcare) according to the standard clinical scanning protocols at our hospital. Two types of MRI sequences, (1) contrast-enhanced T1-weighted spin-echo image (T1-CE): repetition time/echo time (TR/TE), 2283.9/24.3 msec confirming GBM recurrence; (2)FLAIR: repetition time/echo time (TR/TE), 7502.0/122.1 msec indicating peritumoral edema regions pre-operation, were selected for the study.

Preprocessing of the images was conducted before analysis to reduce bias. The regions of peritumoral edema were delineated on the FLAIR images. The lesions of recurrence were delineated on contrast-enhanced T1-weighted images. All these were conducted by an experienced neurosurgeon using 3D Slicer 5.2.0 (http://www.slicer.org) operating on Mac OS X. ‘Registration’ function of 3D Slicer was conducted firstly to eliminate different sequences shift induced by the movement of the patient’s head during MR scanning. That way, the two MRI images, contrast-enhanced T1-weighted and FLAIR, of each patient were guaranteed to locate on the same slicer. Subsequently, with reference to the location of the recurrent tumor on contrast-enhanced T1-weighted images, the peritumoral edema regions were divided into two subregions on the FLAIR images: recurrence and non-recurrence regions. Figure 1 shows representative contrast-enhanced T1-weighted and FLAIR images with annotations for the tumor and the recurrence and non-recurrence regions. The area in the green coil represents the primary tumor in T1-CE image, the area in the red and blue coil represents the whole peritumoral edema regions before first operation, the area in the red coil represents the location tumor recurrence in peritumoral edema while the area in the blue coil represents other peritumoral edema regions in FLAIR sequences.

FIGURE 1
www.frontiersin.org

Figure 1 Segmentation of peritumoral edema. The annotations for tumors as delineated on contrast-enhanced T1 weighted images are outlined in green, while the annotations for recurrence and non-recurrence regions within the peritumoral edema as delineated on fluid-attenuated inversion recovery images are outlined in blue and red, respectively.

Radiomic feature extraction and selection

Quantitative radiomic features were extracted from the recurrence and non-recurrence regions on the FLAIR images by 3D Slicer using the “Slicer Radiomics” plugin. 94 texture features that could be divided into five groups were extracted from MR images. Group 1 (shape) consisted of 16 descriptors that reflect the morphological characteristics of tumor and recurrence and non-recurrence regions. Group 2 (first order) consisted of 19 descriptors that describe the distribution of voxel intensities by frequently-used basic metrics within the MR image. Group 3 (GLCM) consisted of 26 descriptors. Group 4 (GLRLM) consisted of 16 descriptors. Group 5 (GLSZM) consisted of 16 descriptors. The descriptors in group 3,group 4 and group 5 were calculated by specific matrix algorithm.

To standardize the feature value, we converted it into standardized value according to the formula |a−b|/a×100%, where “a” stands for the value of the entire peritumoral edema region, while in analysis radiomics features from normal brain region, ‘‘a’’ represents the feature values extracted from the entire normal brain region, and “b” stands for the value of the subregions (recurrence or non-recurrence). Subsequently, the standardized values were used to identify the features with a significant difference between the recurrence and non-recurrence regions using the paired t test.

Statistical analysis

Survival analysis was performed using Kaplan-Meier’s method. Radiomics features were examined whether followed a gaussian distribution by Shapiro-Wilk test. Then, the radiomic characteristics of the two groups were compared using paired t test. The differences were considered significant at p < 0.05. All analyses were performed using IBM SPSS Statistics v23.0 (IBM Corp.).

Results

Clinicopathological characteristics

A total of 22 patients with recurrent GBM were included in this study,among which 19 (86%) was male while 3 (14%) were female. The median age at diagnosis was 49 (range, 23–67).

Regarding the tumor location, 3 (14%) patients had tumors in the parietal lobe, 6 (27%) in the frontal lobe, 7 (32%) in the temporal lobe, 2 (9%) in the occipital lobe, and 3 (14%) in other regions. The Ki-67 proliferation index was lower than 50% in 13 (59%) patients, whereas 2 (9%) patients had no available data. The O6-methylguanine-DNA methyltransferase and IDH1 status was positive in 11 (50%) and 12 (55%) patients, and negative in 10 (46%) and 9 (41%) patients, respectively (Table 1).

TABLE 1
www.frontiersin.org

Table 1 The clinical characteristics of the 22 patients with recurrent glioblastoma multiforme.

PFS and OS

It is well known that patients with a recurrence GBM have a much poor prognosis, PFS and OS were calculated to indicate its prognosis. The 1-year progression-free survival was 22.7% (95% confidence interval, 13.8–31.6), and the 1-year overall survival was 60% (95% confidence interval, 48.2–71.8) for the entire study population, with a median follow-up of 529 (range 0–1,319) days (Figure 2).

FIGURE 2
www.frontiersin.org

Figure 2 Kaplan-Meier curves for overall survival (A) and progression-free of the included patients (B).

Radiomic feature selection

A set of 94 radiomic features was extracted from the recurrence and non-recurrence regions. The following 10 features demonstrated a significant difference between the recurrence and non-recurrence regions: difference entropy, short run low gray-level emphasis (15), long run low gray-level emphasis (15), run variance (16), long-run emphasis (17), low gray-level run emphasis (18), size-zone non-uniformity normalized (16), small-area emphasis (17), low gray-level zone emphasis (18),and major axis. These features belong to four different texture indexes, namely, gray-level co-occurence matrix(GLCM), gray-level run-length matrix (GLRLM), gray-level size-zone matrix (GLSZM), and shape (Table 2; Figure 3). A schematic diagram presenting the findings of our study and how the radiomic features with significant difference were calculated is shown in Figure 4.

TABLE 2
www.frontiersin.org

Table 2 Radiomic features that may be relevant for predicting the site of tumor recurrence.

FIGURE 3
www.frontiersin.org

Figure 3 Ten features with a significant difference between the recurrence and non-recurrence regions within the peritumoral edema.

FIGURE 4
www.frontiersin.org

Figure 4 A schematic diagram presenting the findings of our study.

Verification

In order to verify the effectiveness of the 10 selected features to distinguish between recurrence and non-recurrence regions, we analyzed the FLAIR images of two patients from TCGA-GBM database with the same method as described above. Based on the hypothesis that the radiomic features of non-recurrence regions are identical to those of normal tissue outside tumor and edema, a two-sample t-test was performed to determine whether there was a significant difference in the 10 selected features between the radiomic data extracted from normal brain regions of patients in TCGA database and the non-recurrence regions of our patients. We found that seven of the 10 features showed no significant difference (P > 0.05),the 7 characteristics that were significantly different are difference entropy, short run low gray-level emphasis, run variance, long-run emphasis, small-area emphasis, low gray-level zone emphasis and major axis.

The acquisition of radiomic features based on GLRLM and GLSZM requires partition of the MR images (20). Some data of the recurrence regions were inevitably involved in the calculation of limbic non-recurrence regions after partition, due to the small sample size. Furthermore, the limbic voxels of the field of interest may contain normal tissue, and parts of the tumor region. For all the above reason, we believe some bias may be introduced into the calculating of radiomic features, which can be amplified by small sample patients sizes. Despite the fact that three features showed a significant difference, we still suggest that the result of our verification supports an association between the 10 identified radiomic features and the recurrence of GBM.

Discussion

In terms of preoperative prediction for tumor recurrence, although various techniques and methods have been developed to identify different subregions (with heterogeneity) within the peritumoral edema, none of them has been routinely used with high efficacy in clinical practice (21, 22). For better prognosis of GBM, it is necessary to develop techniques of noninvasive preoperative examination for tumor recurrence. Among these, radiomics has superiority in pathological grading of gliomas and genotype definition. The molecular marker status could be reflected by imaging phenotypes, which contributed to the promotion of treatment planning before surgery and medical management (23). Similarly, radiomic features based on multi-contrast MRI can provide information correlating with the glioma grade and Ki-67 labeling index (24). It was also shown that subtypes of GBM can be identified by anatomic MRI with the biomarker volume ratio (25). Since maximum safe resection is the normal surgical treatment, a more accurate definition of the tumor volume and regions in high risk of infiltration can be developed on the base of spectroscopic MRI biomarkers, which associated radiomics with GBM treatment (26). In addition, quantitative analysis of features from MR images could reveal predictive imaging biomarkers via intratumoral segmentation and comparison of multimodality MR images (27). Thus, radiomics has been proven to be a bridge between imaging features and clinical treatment.

Accordingly (28), the radiomics features mainly belong to four subgroups, shape, GLCM, GLRLM and GLSZM. The GLCM depicts the spatial relationship between pixels or voxel pairs with predefined intensities of gray-level in different directions (2D analysis is in horizontal, vertical or diagonal direction, 3D analysis is in 13 directions), and distance between pixels or voxels is predefined. The GLRLM captures information of the 2- or 3-dimensional spatial distribution of continuous pixels with the same gray-level in one or more directions. Based on a calculation principle similar to GLRLM, the GLSZM captures information about spatial distribution of interconnected neighboring pixels or voxels with the same gray-level in one or more directions, in 2 or 3 dimensions. In our study, a set of 94 radiomic features was extracted from the recurrence and non-recurrence regions. Radiomics has also been applied in many studies as a preoperative noninvasive method to confirm the relationship between peritumoral edema and the prognosis of glioma (29, 30). As conventional MRI has some limitations in detection, Lemercier et al. verified the effectiveness of the apparent diffusion coefficient value gradient in differentiating glioma from metastatic lesions by means of diffusion-weighted MRI (30).

In this context, we used radiomics to investigate the heterogeneity of the peritumoral edema, considering that gliomas potentially recur in this area. Prasanna et al. extracted features from multiparametric MR images of GBM and found 10 radioactive subsets of “peritumoral” MR features that suggested intensity heterogeneity and texture patterns in the peritumoral brain region (31). In the present study, we linked radiomics with predicting the potential sites of tumor recurrence within the peritumoral edema and successfully verified that 7 radiomic features of non-recurrence regions and normal brain tissue share a relatively high similarity. Our study exhibited the potential of applying radiomics to the optimization of surgical resection for GBM.

In this study, we investigated whether there was heterogeneity within the peritumoral edema region and whether using radiomic features that display heterogeneity could roughly predict the site of GBM recurrence. Our results suggested that heterogeneity does exist in peritumoral edema, indicating that radiomic features of the peritumoral edema from routine MRI can be utilized to predict sites of GBM recurrence.

We first extracted 94 features from the peritumoral edema regions on radiological images, and through comparison between the recurrence and non-recurrence regions in the peritumoral edema, we selected 10 features exhibiting significant differences, suggesting that the peritumoral edema has certain heterogeneity. This finding is consistent with those of previous studies. At histological and cellular level, through computer-guided stereotaxic biopsies, Clavreul et al. isolated a new type of stromal cells from the peritumoral zone, named glioblastoma-associated stromal cells (1, 2). Another cell type presenting specifically in the periphery of GBM are reactive astrocytes, which are located around the tumor and may exert the function of promoting tumor growth and division (32, 33). A significant difference between the peritumoral tissue and normal brain tissue also exists at molecular level. Piwecka et al. used microarrays and sequencing to investigate the micro-ribonucleic acid (miRNA) profile in glioma tissue and peritumoral brain regions, and compare it with that of normal brain tissue. The results showed differences in the expression of 97 miRNAs between the glioma tissue and normal brain tissue (34). Similarly, Mangiola et al. used an array to detect genomic alterations and found that, compared with normal brain tissue, 15 genes were over-expressed, while 42 genes were down-regulated (35). These results suggest that although they are not as obvious as those between tumors and normal brain tissue, the differences between the peritumoral tissue and normal brain tissue cannot be ignored as well.

We also verified our results using images of two patients from TCGA database, the result of our verification supports an association between the 10 identified radiomic features and the recurrence of GBM. The purpose of our study was to find radiomics characteristics with significant difference for further prediction model establishment. Previous studies have also suggested that peritumoral edema is associated with tumor recurrence (1, 2). The downregulation of tumor suppressor genes in the peritumoral edema also suggested a correlation between peritumoral edema and tumor recurrence (35). Basen on these radiomic features acquired in this study, in another research conducted by our team, we successfully identified that satellite lesions within the peritumoral edema region of dysembryoplastic neuroepithelial tumor patients may elicit epilepsy recurrence, and radiomics could be used for detecting and evaluating these epilepsy-associated lesion (36). The differences between the peritumoral edema and normal brain tissue and those within the peritumoral edema suggest that the identification of specific regions in the peritumoral edema with relapse-related biomarkers is the basis for the prediction of tumor recurrence.

Limitations

The major limitations of this study were the limited sample size and the single-institution design, which may reduce the efficacy of the statistical results. As a malignant central nervous system tumor, patients with GBM have a poor prognosis and low PFS and OS. Thus, the total number of patients meet our main selection criteria: pathologically confirmed diagnosis of GBM and pathologically and image confirmed tumor recurrence meanwhile accepted secondary surgey is limited. With regard to the verification, in TCGA database, the patients meet our main selection criteria mentioned above is limited and we only found two subjects available for further validation. In addition, all ROIs were conducted by one experienced neurosurgeon in this study, however, interobserver variability should be taken into account. Therefore, multicenter studies with a larger, double-blind cohort meet our selection criteria standards are required to further verify and complement our results. In addition, this was a retrospective study; thus, further prospective studies are needed to validate the clinical value of our prediction model for GBM management. In the new classification, all glioblastomas are wild-type, so there is no such entity as IDH mutant positive for glioblastoma. These are reclassified as Astrocytomas, WHO III or IV, based on histology, the patients with GBM who were diagnosis as IDH mutant type may need to be separately analyzed further.

Conclusions

In summary, radiomics is a valuable method to determine the recurrence-related heterogeneity of peritumoral edema in order to better predict the site of tumor recurrence and help to determine the surgical resection range more accurately. However, further prospective multicenter cohort studies are required to evaluate the sensitivity and specificity of the radiomic features we have identified, which may be helpful to further guide surgical treatment to improve the prognosis of GBM.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving human participants were reviewed and approved by Ethics approval were obtained from the Ethics Committee of Nanfang Hospital Southern Mecial University. All the experiment protocol for involving human data was in accordance with the guidelines of Declaration of Helsinki. And the written consent was obtained from study participants. The patients/participants provided their written informed consent to participate in this study.

Author contributions

HL: Research project execution, statistical analysis execution. PZ: Research project organization, statistical analysis design, manuscript review and critique. YB Statistical analysis design & execution, manuscript writing of the first draft & review and critique. CY: Statistical analysis design & review and critique. MW: Research project execution, Statistical analysis execution. DH: Research project execution, Statistical analysis review and critique. SH: Statistical analysis review and critique. KY: Research project organization, manuscript review and critique. SQ: Research project organization, manuscript review and critique. JW: Research project conception & organization, statistical analysis review and critique, manuscript review and critique. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by grants from the Natural Science Fund of Guangdong Province (No. 2017A030313597), “Climbing Program” Special Fund of Guangdong Province (No. pdjh2019b0100, No. pdjh2020b0112) and Southern Medical University (No. LX2016N006, No. KJ20161102, No.201912121004S, No.201912121013, No. S202012121088, No. X202012121354).

Acknowledgments

The authors thank the patients and their families.

Conflict of interest

The authors of this manuscript declare no relationships with any companies, whose products or services may be related to the subject matter of the article.

Publisher’s note

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.

Abbreviations

3D, three-dimensional; FLAIR, fluid-attenuated inversion recovery; GBM, glioblastoma multiforme; GLRLM, gray-level run-length matrix; GLSZM, gray-level size zone matrix; miRNAs, micro-ribonucleic acids; MRI, magnetic resonance imaging; TCGA, The Cancer Genome Atlas.

References

1. Clavreul A, Guette C, Faguer R, Tétaud C, Boissard A, Lemaire L, et al. Glioblastoma-associated stromal cells (GASCs) from histologically normal surgical margins have a myofibroblast phenotype and angiogenic properties. J Pathol (2014) 233:74–88. doi: 10.1002/path.4332

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Clavreul A, Etcheverry A, Chassevent A, Quillien V, Avril T, Jourdan ML, et al. Isolation of a new cell population in the glioblastoma microenvironment. J Neurooncol (2012) 106:493–504. doi: 10.1007/s11060-011-0701-7

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Stupp R, Mason WP, Van Den Bent MJ, Weller M, Fisher B, Taphoorn MJ, et al. Radiotherapy plus concomitant and adjuvant temozolomide for glioblastoma. New Engl J Med (2005) 352:987–96. doi: 10.1056/NEJMoa043330

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Petrecca K, Guiot M, Panet-Raymond V, Souhami L. Failure pattern following complete resection plus radiotherapy and temozolomide is at the resection margin in patients with glioblastoma. J Neuro Oncol (2013) 111:19–23. doi: 10.1007/s11060-012-0983-4

CrossRef Full Text | Google Scholar

5. Lemée J, Clavreul A, Menei P. Intratumoral heterogeneity in glioblastoma: don’t forget the peritumoral brain zone. Neuro-oncology (2015) 17:1322–32. doi: 10.1093/neuonc/nov119

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Engelhorn T, Savaskan NE, Schwarz MA, Kreutzer J, Meyer EP, Hahnen E, et al. Cellular characterization of the peritumoral edema zone in malignant brain tumors. Cancer Sci (2009) 100:1856–62. doi: 10.1111/j.1349-7006.2009.01259.x

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Dehnhardt M, Zoriy MV, Khan Z, Reifenberger G, Ekström TJ, Sabine Becker J, et al. Element distribution is altered in a zone surrounding human glioblastoma multiforme. J Trace Elements Med Biol (2008) 22:17–23. doi: 10.1016/j.jtemb.2007.08.002

CrossRef Full Text | Google Scholar

8. Lemée JM, Clavreul A, Aubry M, Com E, de Tayrac M, Eliat PA, et al. Characterizing the peritumoral brain zone in glioblastoma: a multidisciplinary analysis. J Neurooncol (2015) 122:53–61. doi: 10.1007/s11060-014-1695-8

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Abrol S, Kotrotsou A, Salem A, Zinn PO, Colen RR. Radiomic phenotyping in brain cancer to unravel hidden information in medical images. Top Magn Reson Imaging (2017) 26:43–53. doi: 10.1097/RMR.0000000000000117

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Yip SS, Aerts HJ. Applications and limitations of radiomics. Phys Med Biol (2016) 61:R150–66. doi: 10.1088/0031-9155/61/13/R150

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Wang K, Wang Y, Fan X, Wang J, Li G, Ma J, et al. Radiological features combined with IDH1 status for predicting the survival outcome of glioblastoma patients. Neuro Oncol (2016) 18:589–97. doi: 10.1093/neuonc/nov239

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Yu J, Shi Z, Lian Y, Li Z, Liu T, Gao Y, et al. Noninvasive IDH1 mutation estimation based on a quantitative radiomics approach for grade II glioma. Eur Radiol (2017) 27:3509–22. doi: 10.1007/s00330-016-4653-3

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Zhou M, Chaudhury B, Hall LO, Goldgof DB, Gillies RJ, Gatenby RA. Identifying spatial imaging biomarkers of glioblastoma multiforme for survival group prediction. J Magn Reson Imaging (2017) 46:115–23. doi: 10.1002/jmri.25497

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Wangaryattawanich P, Hatami M, Wang J, Thomas G, Flanders A, Kirby J, et al. Multicenter imaging outcomes study of the cancer genome atlas glioblastoma patient cohort: imaging predictors of overall and progression-free survival. Neuro Oncol (2015) 17:1525–37. doi: 10.1093/neuonc/nov117

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Dasarathy BV, Holder EB. Image characterizations based on joint gray level–run length distributions. Pattern Recognition Lett (1991) 12:497–502. doi: 10.1016/0167-8655(91)80014-2

CrossRef Full Text | Google Scholar

16. Thibault G, Fertil B, Navarro C, Pereira S, Mari JL. Texture indexes and Gray level size zone matrix application to cell nuclei classification, in: 10th International Conference on Pattern Recognition and Information Processing, PRIP, 2009 (Minsk, Belarus). (2009) 140–5.

Google Scholar

17. Galloway MM. (1975). Texture analysis using gray level run lengths. Computer Graphics and Image Processing 4:172–9. doi: 10.1016/S0146-664X(75)80008-6

CrossRef Full Text | Google Scholar

18. Chu A, Sehgal CM, Greenleaf JF. Use of gray value distribution of run lengths for texture analysis. Pattern Recognition Lett (1990) 11:415–9. doi: 10.1016/0167-8655(90)90112-F

CrossRef Full Text | Google Scholar

19. Haralick RM, Shanmugam K, Dinstein I. Textural features for image classification. IEEE Transactions on Systems, Man, and Cybernetics (1973) SMC-3:610–21. doi: 10.1109/TSMC.1973.4309314

CrossRef Full Text | Google Scholar

20. Krishnan KR, Sudhakar R. Automatic classification of liver diseases from ultrasound images using GLRLM texture features. (Berlin Heidelberg: Springer) (2013).

Google Scholar

21. Martirosyan NL, Cavalcanti DD, Eschbacher JM, Delaney PM, Scheck AC, Abdelwahab MG, et al. Use of in vivo near-infrared laser confocal endomicroscopy with indocyanine green to detect the boundary of infiltrative tumor. J Neurosurg (2011) 115:1131–8. doi: 10.3171/2011.8.JNS11559

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Böhringer HJ, Lankenau E, Stellmacher F, Reusche E, Hüttmann G, Giese A. Imaging of human brain tumor tissue by near-infrared laser coherence tomography. Acta Neurochir (Wien) (2009) 151:507–17. doi: 10.1007/s00701-009-0248-y

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Han Y, Xie Z, Zang Y, Zhang S, Gu D, Zhou M, et al. Non-invasive genotype prediction of chromosome 1p/19q co-deletion by development and validation of an MRI-based radiomics signature in lower-grade gliomas. J Neurooncol (2018) 140:297–306. doi: 10.1007/s11060-018-2953-y

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Li Y, Liu X, Qian Z, Sun Z, Xu K, Wang K, et al. Genotype prediction of ATRX mutation in lower-grade gliomas using an MRI radiomics signature. Eur Radiol (2018) 28:2960–8. doi: 10.1007/s00330-017-5267-0

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Naeini KM, Pope WB, Cloughesy TF, Harris RJ, Lai A, Eskin A, et al. Identifying the mesenchymal molecular subtype of glioblastoma using quantitative volumetric analysis of anatomic magnetic resonance images. Neuro Oncol (2013) 15:626–34. doi: 10.1093/neuonc/not008

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Cordova JS, Shu HK, Liang Z, Gurbani SS, Cooper LA, Holder CA, et al. Whole-brain spectroscopic MRI biomarkers identify infiltrating margins in glioblastoma patients. Neuro Oncol (2016) 18:1180–9. doi: 10.1093/neuonc/now036

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Cui Y, Tha KK, Terasaka S, Yamaguchi S, Wang J, Kudo K, et al. Prognostic imaging biomarkers in glioblastoma: Development and independent validation on the basis of multiregion and quantitative analysis of MR images. Radiology (2016) 278:546–53. doi: 10.1148/radiol.2015150358

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Mayerhoefer ME, Materka A, Langs G, Häggström I, Szczypiński P, Gibbs P, et al. Introduction to radiomics. J Nucl Med (2020) 61:488–95. doi: 10.2967/jnumed.118.222893

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Schoenegger K, Oberndorfer S, Wuschitz B, Struhal W, Hainfellner J, Prayer D, et al. Peritumoral edema on MRI at initial diagnosis: an independent prognostic factor for glioblastoma? Eur J Neurol (2009) 16:874–8. doi: 10.1111/j.1468-1331.2009.02613.x

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Wu CX, Lin GS, Lin ZX, Zhang JD, Liu SY, Zhou CF, et al. Peritumoral edema shown by MRI predicts poor clinical outcome in glioblastoma. World J Surg Oncol (2015) 13:97. doi: 10.1186/s12957-015-0496-7

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Prasanna P, Patel J, Partovi S, Madabhushi A, Tiwari P. Radiomic features from the peritumoral brain parenchyma on treatment-naïve multi-parametric MR imaging predict long versus short-term survival in glioblastoma multiforme: Preliminary findings. Eur Radiol (2017) 27:4188–97. doi: 10.1007/s00330-016-4637-3

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Lin Q, Balasubramanian K, Fan D, Kim SJ, Guo L, Wang H, et al. Reactive astrocytes protect melanoma cells from chemotherapy by sequestering intracellular calcium through gap junction communication channels. Neoplasia (2010) 12:748–54. doi: 10.1593/neo.10602

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Fitzgerald DP, Palmieri D, Hua E, Hargrave E, Herring JM, Qian Y, et al. Reactive glia are recruited by highly proliferative brain metastases of breast cancer and promote tumor cell colonization. Clin Exp Metastasis (2008) 25:799–810. doi: 10.1007/s10585-008-9193-z

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Piwecka M, Rolle K, Belter A, Barciszewska AM, Żywicki M, Michalak M, et al. Comprehensive analysis of microRNA expression profile in malignant glioma tissues. Mol Oncol (2015) 9:1324–40. doi: 10.1016/j.molonc.2015.03.007

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Mangiola A, Saulnier N, De Bonis P, Orteschi D, Sica G, Lama G, et al. Gene expression profile of glioblastoma peritumoral tissue: an ex vivo study. PloS One (2013) 8:e57145.

PubMed Abstract | Google Scholar

36. Wang J, Luo X, Chen C, Deng J, Long H, Yang K, et al. Preoperative MRI for postoperative seizure prediction: a radiomics study of dysembryoplastic neuroepithelial tumor and a systematic review. Neurosurg Focus (2022) 53:E7.

PubMed Abstract | Google Scholar

Keywords: glioblastoma, recurrence, radiomics, peritumoral edema, heterogeneity

Citation: Long H, Zhang P, Bi Y, Yang C, Wu M, He D, Huang S, Yang K, Qi S and Wang J (2023) MRI radiomic features of peritumoral edema may predict the recurrence sites of glioblastoma multiforme. Front. Oncol. 12:1042498. doi: 10.3389/fonc.2022.1042498

Received: 12 September 2022; Accepted: 02 December 2022;
Published: 04 January 2023.

Edited by:

Yu Zhu, Zhejiang University, China

Reviewed by:

YI YANG, Central South University, China
Xiaorong Yan, First Affiliated Hospital of Fujian Medical University, China
Chen Jiang, University of Science and Technology of China, China

Copyright © 2023 Long, Zhang, Bi, Yang, Wu, He, Huang, Yang, Qi and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jun Wang, smuwangjun@163.com

These authors have contributed equally to this work

ORCID: Jun Wang, orcid.org/0000-0003-3509-0790

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.