- 1Department of Radiology, Xiangya Hospital, Central South University, Changsha, China
- 2National Clinical Research Center for Geriatric Disorders, Changsha, China
- 3Hunan Engineering Research Center of Skin Health and Disease, Changsha, China
- 4Hunan Key Laboratory of Skin Cancer and Psoriasis, Changsha, China
- 5National Clinical Research Center for Geriatric Disorders (Xiangya Hospital), Central South University, Changsha, China
- 6Department of Pharmaceuticals Diagnosis, GE Healthcare, Hangzhou, China
- 7Department of Pathology, Xiangya Hospital, Central South University, Changsha, China
- 8Department of Neurosurgery, Xiangya Hospital, Central South University, Changsha, China
- 9Department of Diagnostic Radiology, City of Hope National Medical Center, Duarte, CA, United States
Purpose: Early recurrence of glioblastoma after standard treatment makes patient care challenging. This study aimed to assess preoperative magnetic resonance imaging (MRI) radiomics for predicting early recurrence of glioblastoma.
Patients and Methods: A total of 122 patients (training cohort: n = 86; validation cohort: n = 36) with pathologically confirmed glioblastoma were included in this retrospective study. Preoperative brain MRI images were analyzed for both radiomics and the Visually Accessible Rembrandt Image (VASARI) features of glioblastoma. Models incorporating MRI radiomics, the VASARI parameters, and clinical variables were developed and presented in a nomogram. Performance was assessed based on calibration, discrimination, and clinical usefulness.
Results: The nomogram consisting of the radiomic signatures, the VASARI parameters, and blood urea nitrogen (BUN) values showed good discrimination between the patients with early recurrence and those with later recurrence, with an area under the curve of 0.85 (95% CI, 0.77-0.94) in the training cohort and 0.84 [95% CI, 0.71-0.97] in the validation cohort. Decision curve analysis demonstrated favorable clinical application of the nomogram.
Conclusion: This study showed the potential usefulness of preoperative brain MRI radiomics in predicting the early recurrence of glioblastoma, which should be helpful in personalized management of glioblastoma.
Introduction
Glioblastoma multiforme (GBM) is the most malignant primary brain tumor (1) and represents one third of primary brain tumors with 79,000 new cases worldwide per year (2). The standard treatment for newly diagnosed GBM is maximal surgical resection followed by radiotherapy plus temozolomide (3), which takes about six months (4–6). Nevertheless, early recurrence may occur due to the aggressiveness and diffuse infiltrative growth of GBM (7, 8).
Pretreatment identification of patients at risk for early GBM recurrence has several benefits (9). First, a more aggressive treatment strategy, such as a more extensive resection (10), with concurrent extended individualized radiotherapy or new radiotherapy-based methods, may be warranted. Second, patients may need to choose new treatment early in the course of the disease to avoid delays in treatment (11, 12). Finally, patients with early occurrence may need additional in-depth testing, such as gene sequencing, to assist in clinical decision-making. Previous studies have identified various indicators associated with risk of GBM recurrence such as age, gender, Karnofsky performance status (KPS), and pathology findings (13–17), but no reliable tool is available to accurately predict early recurrence. Molecular pathology and genotyping may be helpful in assessing recurrence (18, 19). However, tissue collection is invasive and may lead to misdiagnosis due to errors in the sampling of heterogenous tumors. Clinicians typically use brain MRI to evaluate the visual radiological features of GBM such as the size, location, edema and enhancement characteristics. However, conventional MRI evaluation is not adequate for predicting early recurrence in GBM (19, 20). Therefore, there is a need to assess additional imaging biomarkers analyzed by computational methods for predicting GBM recurrence (21, 22).
Radiomics quantifies high-dimensional imaging features of tumors (23, 24), which may reflect tumor heterogeneity and molecular pathology. Radiomics has been used to evaluate GBM recurrence because a large number of image features can be extracted from brain MRI data, including the contrast enhanced sequence, fluid attenuation inversion recovery (FLAIR), T1-weighted imaging, and T2-weighted imaging (23–26). In addition, GBM is a highly vascularized tumor, and the extent of vascularization is directly associated with the prognosis and recurrence for patients with GBM (24, 25, 27, 28). Hence, it is reasonable to speculate that MRI radiomics, taking into consideration of the various imaging features including enhancing characteristics reflecting tumor vascularity, could be a potentially useful tool to predict early recurrence of GBM.
In this study, we evaluated preoperative brain MRI scans for both radiomic features analyzed by computational methods and the conventional radiological characteristics included in the Visually Accessible Rembrandt Images (VASARI) feature set assessed by neuroradiologists through visual inspection (6). We built a nomogram incorporating radiomic features, clinical variables, and the VASARI parameters to predict early recurrence of GBM. We hypothesized that integration of conventional radiological parameters and clinical variables into the computational radiomic model could improve the model performance for differentiating patients with and without early recurrence of GBM
Patients and Methods
Patients
Ethical approval was obtained from the Ethics Committee and Institutional Review Board of Xiangya Hospital of Central South University, P. R. China (IRB number: 201607831). Informed consent was waived due to the retrospective nature of the study.
The study cohort consisted of all patients between July 2010 and April 2018 who had GBM that was pathologically confirmed from surgical specimens. The patients were identified by searching our institutional data base and medical chart review. We included only the patients with preoperative brain MRI scans obtained less than 14 days prior to surgery. Patients with missing clinical data and those who received treatment before surgery, such as radiotherapy, chemotherapy, or chemoradiotherapy, were excluded. Patients with a progression-free survival (PFS) less than 7 months were assigned to the early recurrence group (7, 8), and the remaining patients were assigned to the later recurrence group. Figure 1 presents a study flow diagram.
Re-Assessment of Pathology Results
All pathology slides for the study cohort were re-assessed by two pathologists specializing in brain tumors to confirm GBM diagnosis (G.G. and Z.H., with 10 and 25 years of experience, respectively). Both pathologists were blinded to the clinicopathological data. If any discrepancies arose, agreements were reached by consensus.
Tumor Segmentation and Radiomic Feature Extraction
Brain MRI imaging was performed on a 3T MRI scanner (Discovery MR750w, GE Healthcare) or a 1.5T MRI scanner (MAGNETOM Avanto, Siemens Healthineers, Erlangen, Germany). All images were retrieved from the Picture Archiving and Communication System (PACS) at our hospital. All brain MRI scans were reviewed independently by two neuroradiologists [reader 1 (H.D.) and reader 2 (Z.H.), with 5 and 25 years of experience, respectively] who were blinded to the patient information. Any disagreements were resolved in a panel format with two additional researchers (X.Y. and H.T.).
All conventional radiological findings from the brain MRI images, as assessed by the neuroradiologists visually, were evaluated according to the VASARI parameters (https://wiki.nci.nih.gov/display/CIP/VASARI). There were 29 VASARI parameters recorded as F1-F29 (6).
Image pre-processing, tumor segmentation, and radiomic feature extraction were performed as described previously (29). Wavelet transform and Laplacian of Gaussian (LoG) filtering were used for image denoising. For each GBM tumor, manual contouring and segmentation were performed on the axial, sagittal, coronal images of the tumor at the level with the largest tumor dimension, as determined by both neuroradiologists (readers 1 and 2) who defined the margins of the tumors. The segmented tumor with regions of interest (ROI) was saved for subsequent radiomic feature extraction. We extracted reliable features from the original image and its corresponding filtered image. A total of 1204 quantitative radiomic features were extracted from each MRI image including first order statistics, shape, gray-level run-length matrix (GLRLM), gray-level co-occurrence matrix (GCLM), neighboring gray tone difference matrix (NGTDM), gray level dependence matrix (GLDM), and gray-level size zone matrix (GLSZM), using an open-source python package PyRadiomics (2.2.0) (http://www.radiomics.io/pyradiomics.html) (30). To remove the potential differences between MRI images acquired from the two different MRI scanners, normalization with the final 256 bins was performed on all original MRI images using the gray-scale discretization method before extracting the radiomic features (Analysis Kit software, version V3.0.0.R, GE Healthcare). Figure 2 presents the workflow for tumor segmentation, radiomic feature extraction and selection, and predictive modeling.
Radiomic Feature Selection, Predictive Modeling and Statistical Analysis
Max-Relevance and Min-Redundancy (mRMR) was performed to eliminate the redundant and irrelevant features (31), and 20 features were retained. Then, the least absolute shrinkage and selection operator (LASSO) method was used to select the optimized subset of features to construct the radiomic signature (Rad-score) and to build the models. If the Pearson correlation coefficient of the feature pairs was larger than 0.8, one of them was deleted. A 10-fold cross-validation was conducted to avoid model overfitting. The optimal parameter λ was obtained in terms of the largest value of lambda so that the error was within 1 standard error of the minimum. According to λ, the features corresponding to the non-zero LASSO coefficient were selected (see Supplementary Figure 1S). Multivariate logistic regression analysis was applied to develop a model for predicting early recurrence of GBM and the nomogram was generated (Figure 2).
Statistical analysis was conducted with R software (version 3.0.1; http://www.Rproject.org). The reported statistical significance levels were all two-sided, with statistical significance set at 0.05.
Model Performance and the Nomogram
Calibration curves (Hosmer-Lemeshow H test) were used to evaluate the calibration of the models, and receiver-operating characteristic (ROC) curves were performed to access the differentiation efficiency. The performance of the internally validated nomogram was tested in the validation cohort. The logistic regression formula formed in the primary cohort was applied to all patients of the validation cohort. Then, the calibration curve was derived on the basis of the regression analysis, and the clinical usefulness of the models was evaluated by using decision curve analysis in the validation cohort. Decision curve analysis was conducted to determine the clinical usefulness of the radiomic nomogram by quantifying the net benefits at different threshold probabilities in the validation dataset. Additionally, a correlation matrix analysis was performed to evaluate the correlations among all selected features that were integrated into the final predictive model.
Results
Patient Information
A total of 122 patients with GBM were included in this study, including 65 with early recurrence. Median PFS for all patients was 189 ± 19 (range, 26 -1639) days after surgery, with the PFS for early recurrence group and later recurrence group being 97 ± 8 (range, 26 - 203) and 326 ± 31 (range, 217 - 1639) days, respectively. Each patient was randomly assigned to the training cohort (n=86) or the validation cohort (n=36) at a ratio of 5:2. Patient information and comparison between the training and validation cohorts are summarized in Table 1. No variables were statistically different between the training and the validation cohorts (P > 0.05), ensuring a reasonable classification.
Table 1 Demographic, clinical, laboratory, VASARI parameters and radiomic score (Rad-score) of the 122 patients with glioblastoma (GBM).
Inter-Observer and Intra-Observer Reproducibility of Radiomic Features
Satisfactory inter- and intra-observer reproducibility of the texture feature extraction was achieved. The inter-observer intraclass correlation coefficients (ICCs) calculated based on features extracted by reader 1 (first extraction) and reader 2 ranged from 0.791 to 0.897. The intra-observer ICCs, calculated based on reader 1’s feature extraction results from two attempts, ranged from 0.804 to 0.901. Therefore, all subsequent analyses were based on the radiomic features extracted by reader 1.
Predictive Model Building and Model Performance
Of all the radiomic features, 20 potential predictors were retained based on data from the 86 patients in the training cohort. Features with nonzero coefficients in the LASSO logistic regression model were used (Supplementary Figure 1S). The most predictive subset of features was selected and the corresponding coefficients were evaluated (Supplementary Figure 2S).
A significant increase in the Rad-scores for the early recurrence group compared to the later recurrence group was found in the training cohort (P < 0.001), which was further confirmed in the validation cohort (P < 0.005). The radiomic signature yielded an area under the curve (AUC) of 0.81 (95%CI 0.71–0.90) in the training cohort and 0.79 (95%CI 0.64–0.93) in the validation cohort. The ROC curves are presented in Figure 3.
Figure 3 Rad-scores and receiver-operating characteristic (ROC) curves for the early recurrence group and the later recurrence group. (A) Box plots showing the Rad-scores for the early recurrence group and the later recurrence group. The label 0 indicates the later recurrence group and 1 indicates the early recurrence group. The left panel shows the training cohort and the right panel shows the validation cohort. (B) ROC curve for the training cohort. (C) ROC curve for the validation cohort.
Regarding the conventional radiological findings according to the VASARI parameters, we found that F5 for tumor enhancement (Odds Ratio=1.61) and F21 for deep white matter invasion (Odds Ratio=3.45) were risk factors for predicting early recurrence, and F26 indicating the extent of tumor resection (Odds Ratio=0.51) was a protective factor against early recurrence. In our study, F5 and F21 were predictors of early recurrence, and F26 predicted the lack of early recurrence.
Based on the radiomic signature and clinical variables, a model was constructed to predict early recurrence. A logistic regression analysis identified the radiomic signature (F5, F21, F26) and blood urea nitrogen (BUN) concentration as independent predictors (Figures 4A, B), which was presented as a nomogram (Figure 4C). The model with the nomogram provided an AUC of 0.85 (0.77-0.94) in the training cohort and 0.84 (0.71-0.97) in the validation cohort (Figures 4A, B). The calibration curve for the probability of resistance in the training and validation cohorts showed good agreement between prediction and observation.
Figure 4 The receiver operating characteristic (ROC) curves of the predictive models and the corresponding nomogram. ROC curve for the combined model (Nomogram, red), radiomic model (Radiomics, blue), and clinical model (Clinics, green) for (A) the training cohort and (B) the validation cohort. (C) Nomogram with significant indicators. F5, F21, and F26 are part of the Visually Accessible Rembrandt Image (VASARI) feature set. F5, enhancing tumor; F21, deep white matter invasion; F26, extent of tumor resection.
Clinical Applications
The decision curve analysis for both the radiomic model, clinical model, and the nomogram is presented in Figure 5. The decision curve showed that if the threshold probability is greater than 20%, using the nomogram in the current study to predict the recurrence time added more benefit than other forecast schemes or the forecast-none scheme. Within this range, net benefits were comparable on the basis of the radiomic nomogram and clinical model.
Figure 5 Decision curve analysis for the models built with the nomogram (red), radiomic model (Radiomics, blue), clinical model (Clinics, green), and a combination of all models (Nomogram, Red). The y-axis measures the net benefit.
Discussion
In this study, we found that radiomic features derived from preoperative brain MR images were associated with early recurrence in patients with GBM. Our predictive model combining radiomic features, VASARI parameters, and clinical variables could efficiently differentiate the patients with early recurrence from those with later recurrence. Thus, our study results demonstrate the potential usefulness of a non-invasive radiomic approach for predicting early recurrence of GBM.
Radiomics has been used to predict treatment response through assessment of tumor heterogeneity (32). Tumor heterogeneity, which includes variations in histological markers and in the presence of genetic alterations, has been associated with poor clinical outcomes in various diseases, including GBM (33). Heterogeneous tumors are more likely to contain cancer cells that may proliferate faster, may more readily metastasize, and may be more resistant to treatment (33, 34). Therefore, standard treatment with radiotherapy plus temozolomide may not be effective. At the same time, tumor heterogeneity often leads to vascular proliferation (33). A prior study showed that tumor vascularity is a prognostic factor for newly diagnosed glioblastoma (25). Higher tumor heterogeneity may reflect a more aggressive tumor and higher probability to recur (35). Therefore, it was not surprising that our radiomic model, which reflected tumor heterogeneity, achieved robust performance.
Our study found a renal function indicator, BUN, being relevant for predicting early recurrence in patients with GBM. BUN is a commonly used marker of renal function and its serum concentration varies according to glomerular filtration rate (32). High BUN indicates accumulation of ammonia in tissues and blood (33). In general, chemotherapy drugs such as cisplatin have strong renal toxicity and the dosage is determined and adjusted according to the BUN values. Therefore, abnormally elevated BUN may lead to a decrease in the dosage of chemotherapy drugs, which may result in inadequate treatment of the tumor and thus facilitate recurrence (34). In addition, elevated BUN is an indicator of poor physical condition of the patient, who may not be able to fight off tumor recurrence.
Our data showed that several VASARI parameters were independent factors for predicting early recurrence. The VASARI lexicon contains visually analyzable features extracted from routine medically indicated brain MRI and provides standardized visual grading of brain MRI findings for GBM. In our study, F5 and F21 were predictors of early recurrence, and F26 predicted the lack of early recurrence. The F5 parameter in the VASARI feature set indicates tumor enhancement, and a higher score means a higher enhancement ratio with more tumor blood supply and more enhancing tumor. Mathivet et al. (25) observed a progressive increase in vessel diameter during GBM development in an orthotopic mouse model. Therefore, it is reasonable to speculate that increased GBM tumor vascularity may support early recurrence. The VASARI parameter F21 indicates deep white matter invasion, which may pose challenges for complete tumor resection during surgery and render the tumor prone to recur. The VARASI parameter F26 indicates the extent of tumor resection when comparing the preoperative and postoperative MRI images. It is understandable that a more complete tumor resection should improve prognosis and made it less likely to recur.
This study had several limitations. First, this was a retrospective study at a single institution, limiting the generalizability of our study results. In addition, although there were 122 patients with GBM included in our study, our sample size was still modest for a radiomic study, given the heterogeneous nature of GBM. Second, we used 2D texture features of the brain MRI images. A 3D approach for textural features may offer more information about the entire tumor, which may improve predictive model performance (36). Lastly, we did not assess the overall survival rate for this cohort. Clinical data for some patients after recurrence were incomplete because they went back to their local hospitals to continue treatment. In addition, some patients were lost to long-term follow up. It has also been challenging to assess overall survival as patients with recurrent GBM underwent various therapies such as additional surgery, with or without concurrent chemoradiation, or radiation only or chemotherapy with temozolamide only or with addition of bevacizumab, immunotherapy, herbal holistic remedies, etc. With our modest sample size, we did not have the statistical power to tease out the survival rate for patients undergoing different treatments. A future multi-center prospective study of glioblastoma recurrence will be necessary to properly assess survival rate.
In summary, our MRI radiomic analysis and nomogram showed potential value for predicting early recurrence of GBM, which may assist in personalized treatment planning. Future prospective multicenter study with a larger sample size will be needed to validate our study result and to optimize the prediction models for clinical practice.
Data Availability Statement
The original data will be made available to qualified researchers upon request.
Ethics Statement
The studies involving human participants were reviewed and approved by Ethics Committee and Institutional Review Board of Xiangya Hospital of Central South University, P. R. China (IRB number: 201607831). Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author Contributions
Conceived and designed the study: XY, JN, and ZT. Collected and analyzed the data: XY, YF, PP, HD, HT, ZDH, HL, GG, and ZLH. Prepared the manuscript: XY, JW, and BC. All authors contributed to the article and approved the submitted version.
Funding
This study was funded in part by Natural Science Foundation of Hunan Province, P. R. China (2018JJ2641), and China Post-doctoral Science Foundation (2018M632997).
Conflict of Interest
Author PP was employed by company GE Healthcare.
The remaining 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.
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.
Acknowledgments
We thank Dr. Jiezhong Wang and Xiaomei Gao (Xiangya Hospital, Central South University, P.R. China) for helpful discussion and assistance in data analysis. We thank staff members in the Departments of Radiology, Neurosurgery, and Pathology at Xiangya Hospital for their efforts in collecting the information used in this study. We thank our statisticians Dr. Minxue Shen and Dr. Tingfan Wu for their expert statistical review of our manuscript. Editing assistance was provided by our scientific writer Nancy Linford, PhD.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.769188/full#supplementary-material
References
1. Tan SK, Pastori C, Penas C, Komotar RJ, Ivan ME, Wahlestedt C, et al. Serum Long Noncoding RNA HOTAIR as a Novel Diagnostic and Prognostic Biomarker in Glioblastoma Multiforme. Mol Cancer (2018) 17(1):74. doi: 10.1186/s12943-018-0822-0
2. Aldape K, Zadeh G, Mansouri S, Reifenberger G, von Deimling A. Glioblastoma: Pathology, Molecular Mechanisms and Markers. Acta Neuropathol (2015) 129(6):829–48. doi: 10.1007/s00401-015-1432-1
3. Sulman EP, Ismaila N, Armstrong TS, Tsien C, Batchelor TT, Cloughesy T, et al. Radiation Therapy for Glioblastoma: American Society of Clinical Oncology Clinical Practice Guideline Endorsement of the American Society for Radiation Oncology Guideline. J Clin Oncol (2017) 35(3):361–9. doi: 10.1200/JCO.2016.70.7562
4. Shergalis A, Bankhead A 3rd, Luesakul U, Muangsin N, Neamati N. Current Challenges and Opportunities in Treating Glioblastoma. Pharmacol Rev (2018) 70(3):412–45. doi: 10.1124/pr.117.014944
5. Allen BG, Bodeker KL, Smith MC, Monga V, Sandhu S, Hohl R, et al. First-In-Human Phase I Clinical Trial of Pharmacologic Ascorbate Combined With Radiation and Temozolomide for Newly Diagnosed Glioblastoma. Clin Cancer Res (2019) 25(22):6590–7. doi: 10.1158/1078-
6. 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(11):1525–37. doi: 10.1093/neuonc/nov117
7. Brown TJ, Brennan MC, Li M, Church EW, Brandmeir NJ, Rakszawski KL, et al. Association of the Extent of Resection With Survival in Glioblastoma: A Systematic Review and Meta-Analysis. JAMA Oncol (2016) 2(11):1460–9. doi: 10.1001/jamaoncol.2016.1373
8. Gravendeel LA, Kouwenhoven MC, Gevaert O, de Rooi JJ, Stubbs AP, Duijm JE, et al. Intrinsic Gene Expression Profiles of Gliomas are a Better Predictor of Survival Than Histology. Cancer Res (2009) 69(23):9065–72. doi: 10.1158/0008-5472.CAN-09-2307
9. Mansouri A, Hachem LD, Mansouri S, Nassiri F, Laperriere NJ, Xia D, et al. MGMT Promoter Methylation Status Testing to Guide Therapy for Glioblastoma: Refining the Approach Based on Emerging Evidence and Current Challenges. Neuro Oncol (2019) 21(2):167–78. doi: 10.1093/neuonc/noy132
10. Blumenthal DT, Gorlia T, Gilbert MR, Kim MM, Burt Nabors L, Mason WP, et al. Is More Better? The Impact of Extended Adjuvant Temozolomide in Newly Diagnosed Glioblastoma: A Secondary Analysis of EORTC and NRG Oncology/RTOG. Neuro Oncol (2017) 19(8):1119–26. doi: 10.1093/neuonc/nox025
11. Reichel D, Sagong B, Teh J, Zhang Y, Wagner S, Wang H, et al. Near Infrared Fluorescent Nanoplatform for Targeted Intraoperative Resection and Chemotherapeutic Treatment of Glioblastoma. ACS Nano (2020) 14(7):8392–408. doi: 10.1021/acsnano.0c02509
12. Bastiancich C, Bianco J, Vanvarenberg K, Ucakar B, Joudiou N, Gallez B, et al. Injectable Nanomedicine Hydrogel for Local Chemotherapy of Glioblastoma After Surgical Resection. J Control Release (2017) 264:45–54. doi: 10.1016/j.jconrel.2017.08.019
13. Alexander BM, Ba S, Berger MS, Berry DA, Cavenee WK, Chang SM, et al. Adaptive Global Innovative Learning Environment for Glioblastoma: GBM AGILE. Clin Cancer Res (2018) 24(4):737–43. doi: 10.1158/1078-0432.CCR-17-0764
14. Tanguturi SK, Trippa L, Ramkissoon SH, Pelton K, Knoff D, Sandak D, et al. Leveraging Molecular Datasets for Biomarker-Based Clinical Trial Design in Glioblastoma. Neuro Oncol (2017) 19(7):908–17. doi: 10.1093/neuonc/now312
15. Stupp R, Mason WP, van den Bent MJ, Weller M, Fisher B, Taphoorn MJ, et al. Radiotherapy Plus Concomitant and Adjuvant Temozolomide for Glioblastoma. N Engl J Med (2005) 352(10):987–96. doi: 10.1056/NEJMoa043330
16. Gately L, McLachlan SA, Philip J, Ruben J, Dowling A. Long-Term Survivors of Glioblastoma: A Closer Look. J Neurooncol (2018) 136(1):155–62. doi: 10.1007/s11060-017-2635-1
17. Kickingereder P, Burth S, Wick A, Götz M, Eidel O, Schlemmer HP, et al. Radiomic Profiling of Glioblastoma: Identifying an Imaging Predictor of Patient Survival With Improved Performance Over Established Clinical and Radiologic Risk Models. Radiology (2016) 280(3):880–9. doi: 10.1148/radiol.2016160845
18. Phi JH, Park AK, Lee S, Choi SA, Baek IP, Kim P, et al. Genomic Analysis Reveals Secondary Glioblastoma After Radiotherapy in a Subset of Recurrent Medulloblastomas. Acta Neuropathol (2018) 135(6):939–53. doi: 10.1007/s00401-018-1845-8
19. Barboriak DP, Zhang Z, Desai P, Snyder BS, Safriel Y, McKinstry RC, et al. Interreader Variability of Dynamic Contrast-Enhanced MRI of Recurrent Glioblastoma: The Multicenter ACRIN 6677/RTOG 0625 Study. Radiology (2019) 290(2):467–76. doi: 10.1148/radiol.2019181296
20. Lundemann M, Munck Af Rosenschöld P, Muhic A, Larsen VA, Poulsen HS, Engelholm SA, et al. Feasibility of Multi-Parametric PET and MRI for Prediction of Tumour Recurrence in Patients With Glioblastoma. Eur J Nucl Med Mol Imaging (2019) 46(3):603–13. doi: 10.1007/s00259-018-4180-3
21. Tien AC, Li J, Bao X, Derogatis A, Kim S, Mehta S, et al. A Phase 0 Trial of Ribociclib in Recurrent Glioblastoma Patients Incorporating a Tumor Pharmacodynamic- and Pharmacokinetic-Guided Expansion Cohort. Clin Cancer Res (2019) 25(19):5777–86. doi: 10.1158/1078-0432.CCR-19-0133
22. Kim JY, Park JE, Jo Y, Shim WH, Nam SJ, Kim JH, et al. Incorporating Diffusion- and Perfusion-Weighted MRI Into a Radiomics Model Improves Diagnostic Performance for Pseudoprogression in Glioblastoma Patients. Neuro Oncol (2019) 21(3):404–14. doi: 10.1093/neuonc/noy133
23. Choi YS, Ahn SS, Kim DW, Chang JH, Kang SG, Kim EH, et al. Incremental Prognostic Value of ADC Histogram Analysis Over MGMT Promoter Methylation Status in Patients With Glioblastoma. Radiology (2016) 281(1):175–84. doi: 10.1148/radiol.2016151913
24. Yang F, Xie Y, Tang J, Liu B, Luo Y, He Q, et al. Uncovering a Distinct Gene Signature in Endothelial Cells Associated With Contrast Enhancement in Glioblastoma. Front Oncol (2021) 11:683367. doi: 10.3389/fonc.2021.683367
25. Park JE, Kim HS, Jo Y, Yoo RE, Choi SH, Nam SJ, et al. Radiomics Prognostication Model in Glioblastoma Using Diffusion- and Perfusion-Weighted MRI. Sci Rep (2020) 10(1):4250. doi: 10.1038/s41598-020-61178-w
26. Kang D, Park JE, Kim YH, Kim JH, Oh JY, Kim J, et al. Diffusion Radiomics as a Diagnostic Model for Atypical Manifestation of Primary Central Nervous System Lymphoma: Development and Multicenter External Validation. Neuro Oncol (2018) 20(9):1251–61. doi: 10.1093/neuonc/noy021
27. Stadlbauer A, Eyüpoglu I, Buchfelder M, Dörfler A, Zimmermann M, Heinz G, et al. Vascular Architecture Mapping for Early Detection of Glioblastoma Recurrence. Neurosurg Focus (2019) 47(6):E14. doi: 10.3171/2019.9.FOCUS19613
28. Mathivet T, Bouleti C, Van Woensel M, Stanchi F, Verschuere T, Phng LK, et al. Dynamic Stroma Reorganization Drives Blood Vessel Dysmorphia During Glioma Growth. EMBO Mol Med (2017) 9(12):1629–45. doi: 10.15252/emmm.201607445
29. Ma L, Zheng X, Kenli L, Shengli L, Jianlin L, Xiaoping Y. Game Theoretic Interpretability for Learning Based Preoperative Gliomas Grading Future Gener. Comput Syst (2020) 112:1–10. doi: 10.1016/j.future.2020.04.038
30. van Griethuysen JJM, Fedorov A, Parmar C, Hosny A, Aucoin N, Narayan V, et al. Computational Radiomics System to Decode the Radiographic Phenotype. Cancer Res (2017) 77(21):e104–7. doi: 10.1158/0008-5472.CAN-17-0339
31. Pei Q, Yi X, Chen C, Pang P, Fu Y, Lei G, et al. Pre-Treatment CT-Based Radiomics Nomogram for Predicting Microsatellite Instability Status in Colorectal Cancer. Eur Radiol (2021) 10.1007/s00330-021-08167-3. doi: 10.1007/s00330-021-08167-3
32. Sala E, Mema E, Himoto Y, Veeraraghavan H, Brenton JD, Snyder A, et al. Unravelling Tumour Heterogeneity Using Next-Generation Imaging: Radiomics, Radiogenomics, and Habitat Imaging. Clin Radiol (2017) 72(1):3–10. doi: 10.1016/j.crad.2016.09.013
33. Smith SJ, Diksin M, Chhaya S, Sairam S, Estevez-Cebrero MA, Rahman R. The Invasive Region of Glioblastoma Defined by 5ALA Guided Surgery Has an Altered Cancer Stem Cell Marker Profile Compared to Central Tumour. Int J Mol Sci (2017) 18(11):2452. doi: 10.3390/ijms18112452
34. Mroz EA, Tward AD, Pickering CR, Myers JN, Ferris RL, Rocco JW. High Intratumor Genetic Heterogeneity Is Related to Worse Outcome in Patients With Head and Neck Squamous Cell Carcinoma. Cancer (2013) 119(16):3034–42. doi: 10.1002/cncr.28150
35. Andor N, Graham TA, Jansen M, Xia LC, Aktipis CA, Petritsch C, et al. Pan-Cancer Analysis of the Extent and Consequences of Intratumor Heterogeneity. Nat Med (2016) 22(1):105–13. doi: 10.1038/nm.3984
Keywords: blood urea nitrogen, glioblastoma, magnetic resonance imaging, nomogram, preoperative, radiomics, recurrence, Visually Accessible Rembrandt Images (VASARI)
Citation: Wang J, Yi X, Fu Y, Pang P, Deng H, Tang H, Han Z, Li H, Nie J, Gong G, Hu Z, Tan Z and Chen BT (2021) Preoperative Magnetic Resonance Imaging Radiomics for Predicting Early Recurrence of Glioblastoma. Front. Oncol. 11:769188. doi: 10.3389/fonc.2021.769188
Received: 01 September 2021; Accepted: 11 October 2021;
Published: 27 October 2021.
Edited by:
Kevin Camphausen, National Cancer Institute (NCI), United StatesReviewed by:
Weiwei Zong, Henry Ford Health System, United StatesKhaled Elsayad, University of Münster, Germany
Copyright © 2021 Wang, Yi, Fu, Pang, Deng, Tang, Han, Li, Nie, Gong, Hu, Tan and Chen. 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: Xiaoping Yi, yixiaoping@csu.edu.cn; Jilin Nie, doctornjl@163.com; Zeming Tan, doctortzm@163.com
†These authors have contributed equally to this work and share first authorship