Abstract
BACKGROUND AND PURPOSE: Current imaging assessment of high-grade brain tumors relies on the Response Assessment in Neuro-Oncology criteria, which measure gross volume of enhancing and nonenhancing lesions from conventional MRI sequences. These assessments may fail to reliably distinguish tumor and nontumor. This study aimed to classify enhancing and nonenhancing lesion areas into tumor-versus-nontumor components.
MATERIALS AND METHODS: A total of 140 MRI scans obtained from 32 patients with high-grade gliomas and 6 patients with brain metastases were included. Classification of lesion areas was performed using a support vector machine classifier trained on 4 components: enhancing and nonenhancing, tumor and nontumor, based on T1-weighted, FLAIR, and dynamic-contrast-enhancing MRI parameters. Classification results were evaluated by 2-fold cross-validation analysis of the training set and MR spectroscopy. Longitudinal changes of the component volumes were compared with Response Assessment in Neuro-Oncology criteria.
RESULTS: Normalized T1-weighted values, FLAIR, plasma volume, volume transfer constant, and bolus-arrival-time parameters differentiated components. High sensitivity and specificity (100%) were obtained within the enhancing and nonenhancing areas. Longitudinal changes in component volumes correlated with the Response Assessment in Neuro-Oncology criteria in 27 patients; 5 patients (16%) demonstrated an increase in tumor component volumes indicating tumor progression. These changes preceded Response Assessment in Neuro-Oncology assessments by several months. Seven patients treated with bevacizumab showed a shift to an infiltrative pattern of progression.
CONCLUSIONS: This study proposes an automatic classification method: segmented Response Assessment in Neuro-Oncology criteria based on advanced imaging that reliably differentiates tumor and nontumor components in high-grade gliomas. The segmented Response Assessment in Neuro-Oncology criteria may improve therapy-response assessment and provide earlier indication of progression.
ABBREVIATIONS:
- BAT
- bolus arrival time
- DCE
- dynamic contrast-enhanced
- GB
- glioblastoma
- HGG
- high-grade gliomas
- kep
- interstitium-to-plasma rate constant
- Ktrans
- volume transfer constant
- NAWM
- normal-appearing white matter
- nFLAIR
- normalized FLAIR images
- nT1WI+Gd
- normalized T1-weighted images post-contrast agent
- ve
- interstitial volume
- vp
- plasma volume
- RANO
- Response Assessment in Neuro-Oncology
- SPGR
- spoiled gradient-recalled
- sRANO
- segmented RANO
- SVM
- support vector machine
High-grade gliomas (HGG), specifically glioblastoma (GB), remain the most common and aggressive brain tumors in adults. Despite recent advances in treatment, long-term survival remains low.1,2 Treatment of newly diagnosed and recurrent HGG consists of a combination of surgery, radiation therapy, and chemotherapy and, more recently, biologic and immunotherapies, some of which are under clinical investigation. Unique imaging phenomena characterize therapeutic responses to these therapies, challenging conventional radiologic interpretation.
MRI is the method of choice for the initial diagnosis, follow-up, and therapy response assessment of brain tumors. HGG are typically characterized by contrast-enhanced areas on T1-weighted imaging and hyperintensity on FLAIR and T2WI. The current standard for the radiologist's assessment in patients with HGG relies on the Response Assessment in Neuro-Oncology (RANO) criteria,3 which expand upon the earlier Macdonald4 criteria, to incorporate the nonenhancing component of the tumor, as this component may indicate infiltrative or diffuse tumor growth. However, the current RANO criteria, which rely on conventional imaging, may fail to distinguish different tissue components that have similar imaging patterns (ie, contrast-enhanced T1WI) yet different prognostic implications (tumor-versus-treatment effects), the identification of which is crucial for appropriate clinical management.3,5
Glioblastoma shows high heterogeneity within the lesion area, including areas of active tumor, necrosis, infiltrative tumor cells, and vasogenic edema. Enhancing lesion area can be measured on conventional T1WI; however, it may represent a variety of both tumor and nontumor processes, including tumor progression, pseudoprogression, postsurgical changes, and radiation effects.3 Differentiation of the nonenhancing area between infiltrative tumor and nontumor area (edema and/or gliosis) is challenging, as both are characterized by hyperintense FLAIR/T2WI signal. Hyperintense FLAIR/T2WI signal can be expected to represent pure vasogenic edema in metastatic tumors (characteristically with well-defined borders between malignant cells and normal brain),6 while in HGG, the hyperintensity may reflect areas of infiltrating tumor cells and/or vasogenic edema.
Glioblastoma is characterized by a highly complex neovascularization process. The increased vascularity of GB allows analysis and characterization of the tumor using advanced vascular MRI methods, including perfusion and permeability. Dynamic contrast-enhanced (DCE) imaging has been widely used for characterization of tumor biology and therapy-response assessment. DCE relies on dynamic T1WI sequences during bolus injection of contrast agent. Hemodynamic parameters, including plasma volume (vp) and tissue transfer constants (specifically, volume transfer constant [Ktrans] and interstitium-to-plasma rate constant [kep]), can be extracted from DCE data.
Several studies have suggested the utility of DCE parameters as important imaging markers for brain tumor diagnosis, grading, and therapy-response monitoring7⇓–9; to differentiate between active and nonactive tumor components10⇓⇓⇓–14 or infiltrative and noninfiltrative tumor10,15; and for early-stage assessment of antiangiogenic therapy response in HGG.16,17 However, none of these studies have aimed to classify the entire lesion area, both enhancing and nonenhancing, into tumor and nontumor areas at the voxel level.
Various techniques have been proposed for segmentation and classification of brain lesions.18 Support vector machine (SVM) is a supervised binary classifier, originally aimed at classifying data into 2 classes. Based on a training dataset, the algorithm finds the hyperplane that maximally separates between points. New data are then mapped onto the same space and assigned to a particular class. SVM has been widely used as a classifier of medical images and has proved advantageous over other algorithms of its kind.18,19
The aim of this study was to classify lesion areas in patients with HGG using SVM into 4 distinct components: 1) enhancing tumor, 2) enhancing nontumor, 3) nonenhancing tumor, and 4) nonenhancing nontumor. Segmentation of lesion areas into tumor and nontumor components refines the RANO criteria and may improve therapy response MRI assessment.
Materials and Methods
Patients
A total of 140 MRI scans were obtained from 38 patients, who were scanned every 2–3 months. Ninety-five scans were from 26 patients with GB; 31 scans (5 patients) were of anaplastic astrocytomas; 3 scans (1 patient) were of anaplastic oligodendroglias; and 11 scans (6 patients) were of brain metastases (4 breast and 2 lung cancers) (Table 1). Inclusion criteria required a normal glomerular filtration rate and no contraindication to MRI. All patients with HGG who were scanned longitudinally (n = 26) had recurrent HGG and had been treated initially with standard chemoradiation.1 Additional therapies are detailed in the On-line Appendix. Nine patients received bevacizumab at recurrence and had pre- and post bevacizumab MRIs.
Patient characteristics
The study was approved by the Tel Aviv Sourasky Medical Center review board, and written informed consent was obtained from all patients.
MR Imaging Protocol
Scans were performed on 3T MRI scanners: 91 scans were performed on a Signa Excite scanner (GE Healthcare) using an 8-channel head coil, and 49 scans were performed on a Magnetom Prisma scanner (Siemens) using a 20-channel head coil. The protocol included conventional imaging: spoiled gradient-recalled (SPGR)/FLASH T1WI performed before and after contrast agent injection and FLAIR. DCE data were acquired using multiphase 3D T1WI SPGR/FLASH before, during, and after contrast agent injection (FOV 250 mm; matrix 256X256/256X184, section thickness of 5 mm, TR/TE = 5/2.2 ms, and flip angle = 20°. For the T1 maps, variable flip angle SPGR/FLASH data were acquired with nominal flip angles of 5°, 10°, 15°, 20°, and 30°. Dynamic data were acquired with a temporal resolution of 6 seconds and scan duration of 6 minutes. A power injector was used to infuse a single dose (0.2 mL/kg) of contrast agent (gadoterate meglumine, 0.5 mol/L, Dotarem; Guerbet, Aulnay-sous-Bois, France) followed by 20 mL saline flush, both at a constant rate of 5 mL/s second. Fourteen-to-twenty sections were centered on the tumor area as identified in the conventional images, providing brain coverage of 70–100 mm. MR spectroscopy data were acquired with a single-voxel point-resolved sequence by using a cubic volume of ∼8 mL located at the lesion area and normal-appearing white matter (NAWM) in the contralateral hemisphere (TR/TE = 1500/35 ms).
Lesion Classification
Lesion classification involved several steps as detailed in the On-line Appendix and shown in On-line Fig 1. Following preprocessing, the enhancing and nonenhancing lesion areas, designated as targets for classification, were automatically segmented from the T1WI and FLAIR (respectively) images by using an adaptive threshold as previously described.20 Manual correction was performed when necessary. Normalized images of FLAIR (nFLAIR) and postcontrast T1WI (nT1WI+Gd) were calculated relative to NAWM. The DCE pharmacokinetic parameters, vp, Ktrans, interstitial volume (ve), kep, and bolus arrival time (BAT), were estimated by using DCE-Up-Sampled-Temporal-Resolution (DUSTER).21
Classification of the lesion area into 4 components was performed automatically by using the SVM. For the training data, volumes of interest of 5–10 voxels each (30–60 cc) were manually defined with the FMRIB Software Library (FSL; http://www.fmrib.ox.ac.uk/fsl) in areas with known pathology for the 4 tissue types: 1) enhancing tumor, defined retrospectively in patients with progressive HGG based on the RANO assessment on follow-up scans; 2) enhancing nontumor, defined retrospectively in patients with nonprogressive HGG, with treatment-related changes based on stable longitudinal (>6-month) RANO assessments; 3) nonenhancing tumor (infiltrative), defined as a peritumoral area (<2 cm from the enhanced tumor margins) in a nonenhancing FLAIR hyperintense region in patients with HGG with progressive disease; and 4) nonenhancing, nontumor (edema), defined in patients with brain metastasis, in whom the hyperintense FLAIR signal represented pure vasogenic edema, without tumor cell infiltration.6
Seven MRI parameters were measured within and compared between tissue types: nT1WI+Gd, nFLAIR, vp, Ktrans, ve, kep, and BAT. SVM classifier with a linear kernel was trained only on the MRI parameters that were found to be significantly different (P < .05) among the tissue types. Voxelwise classification using SVM was performed separately for the enhancing and nonenhancing lesion areas in all patients.
Evaluation
Two-Fold Cross-Validation Analysis.
To test the sensitivity and specificity of the method with validating confidence intervals, we performed 300 iterations of the training and testing data based on different random selections of the even and odd partitions. The group was divided into 2 subgroups (randomized, even-versus-odd) with 1 subgroup as training data, the other as the test data, and vice versa.
MR Spectroscopy.
Classification results were validated based on MR spectroscopy results, performed off-line using LCModel (http://www.lcmodel.com/).22 Metabolite ratios of glycerophosphocholine + phosphocholine (GPC+PCho) and of N-acetyl aspartate + N-acetyl aspartylglutamate (NAA+NAAG) relative to creatine + phosphocreatine (Cr+PCr) were calculated. While many MR spectra were acquired for clinical purposes, MR spectroscopy data were used in this study when voxel location included only 1 tissue component or were defined within the NAWM.
Longitudinal Assessment.
Longitudinal assessment of the classification results was performed in patients with primary brain tumors (32 patients). Longitudinal changes in the volume of each component were evaluated relative to RANO criteria assessed by a senior neuroradiologist defining 4 categories: complete response, partial response, stable disease, or progression.3
Statistical Analysis
To identify differences between tumor and nontumor components for the various MRI parameters (separately for the enhancing and nonenhancing area), we used the Wilcoxon test (because parameters were not normally distributed).
Results
Classification Results
Training Data.
We defined VOIs for the 4 tissue types: the enhancing tumor component in 28 scans from 5 patients with progressive HGG (4 GB and 1 anaplastic astrocytoma); the enhancing nontumor component in 10 scans from 3 patients with high-grade tumors (2 GB and 1 oligodendroglioma); the nonenhancing tumor component in 24 scans from 5 patients with GB; and the nonenhancing, nontumor component in 11 scans from 6 patients with brain metastasis.
Mean values of the 7 MRI parameters, nT1WI+Gd, nFLAIR, vp, Ktrans, ve, kep, and BAT, measured in the 4 components (tissue types) defined for the training data and NAWM, are provided in Table 2. For the enhancing components, significant differences were detected for the nT1WI+Gd, vp, ktrans, kep, and ve parameters, demonstrating higher vascularity (vp) and impaired permeability (nT1WI+Gd, Ktrans) in the tumor component relative to the nontumor component. For the nonenhancing components, significant differences between components were detected for all parameters, demonstrating higher vascularity (vp) and slightly impaired permeability (nT1WI+Gd, Ktrans) in the tumor component relative to the nontumor component. The ve and kep parameters showed significant group differences yet a substantially high SD (>150%) and were therefore not used in the training set.
Mean and SD of the extracted MR imaging parameters between lesion components
Classification of Lesion Area
The SVM classifier was trained on the nT1WI+Gd, vp, and Ktrans parameters for enhancing components and on the nT1WI+Gd, nFLAIR, vp, Ktrans, and BAT parameters for nonenhancing components. Voxelwise classification was performed on all scans of all patients (n = 140), separately for enhancing and nonenhancing lesion areas. Figure 1 shows classification results obtained in 3 patients. In the first patient (Fig 1, case 1) with breast cancer metastasis, the nonenhancing lesion area was classified as nonenhancing, nontumor (ie, edema) and the enhancing lesion area was classified as enhancing tumor. The second patient (Fig 1, case 2) with GB was an example in whom most of the nonenhancing lesion was classified as tumor component. This patient was diagnosed with progressive disease on a subsequent follow-up scan (2 months later). In the third patient (Fig 1, case 3), also with GB, the enhancing lesion and most of the nonenhancing area were classified as tumor component, consistent with a diagnosis of progressive disease as diagnosed on subsequent follow-up.
Classification results obtained in a patient with active breast cancer metastasis (case 1) and in 2 patients with glioblastoma (cases 2 and 3).
Longitudinal Assessment
Longitudinal assessment of the volume of each component and the correlating radiologic assessment based on RANO, for 26 patients with high-grade lesions scanned longitudinally >3 times (a total of 118 scans) are shown in On-line Fig 2. Consistent results were obtained between changes in the volume of the lesion components and the radiologist's assessment in 27 patients (of 32, 84%). However, in 5 (16%; patients 1, 17, 18, 20, and 26), increased volume of the nonenhancing tumor component indicating tumor progression preceded the radiologic diagnosis based on RANO criteria by several months.
Figure 2 shows longitudinal data obtained from a 54-year-old patient with anaplastic astrocytoma (patient 18). While the conventional imaging showed a pattern of stability at scans 3 and 5, with reduction/no substantial changes in the enhancing and nonenhancing lesion areas, classification results revealed a different clinical scenario–progressive tumor growth. A continual increase in the nonenhancing tumor component indicated a pattern of infiltrative disease progression through scans 2–5 (particularly from scan 2 to 3). This patient was diagnosed with progressive disease at scan 4 and died 71 days following scan 5.
Longitudinal classification and radiologic results obtained in a 54-year-old patient with anaplastic astrocytoma (patient 18) scanned longitudinally every 2 months.
In 7 of 9 patients who received bevacizumab, substantial reductions (mean, 56%) were detected mainly in the volume of the nonenhancing, nontumor component (interpreted as edema). Three patients of 7 (patients 12,13, and 23), though showing a major reduction in the volume of the nonenhancing, nontumor component, had an increase in the nonenhancing tumor component, suggesting a shift to an infiltrative pattern of tumor progression during bevacizumab treatment.
Evaluation
Sensitivity and Specificity of the Training Data.
Analysis demonstrated 100% sensitivity (85%/100% for the 5th/90th percentiles) and 100% specificity (100%/100%) for the identification of enhancing tumor and enhancing nontumor components, and 100% sensitivity (100%/100%) and 100% specificity (100%/100%) for identification of the nonenhancing (infiltrative) tumor and nonenhancing nontumor components.
MR Spectroscopy Results.
Forty-one spectra were obtained with good quality, with voxel location consisting of 1 of the 4 components or within NAWM: 14 from NAWM; 6 from nonenhancing, nontumor component; 1 from nonenhancing tumor; and 19 from enhancing tumor components. MR spectroscopy results are provided in Table 2. Spectra from the enhancing, nontumor components showed primarily lipid peaks, without a detectable metabolite ratio. MR spectroscopy supported the classification results, demonstrating a >1.5-fold increase in glycerolphosphocholine + phosphocholine/creatine + phosphocholine (GPC+PCh/Cr+PCr) ratio in the tumor components relative to the nontumor component and NAWM.
Discussion
In this study, longitudinal classification of lesion area into tumor and nontumor components was performed in patients with HGG based on multi-MRI parameters extracted from conventional imaging and DCE MRI. RANO criteria improved the earlier Macdonald4 criteria by incorporating the nonenhancing component of the lesion. Our proposed method, referred to as segmented RANO (sRANO), provides a logical next step in the evolution of MR clinical imaging, to better assess tumor growth and therapy response. Both enhancing and nonenhancing lesion areas may include tumor and nontumor components. The proposed sRANO classifies the lesion areas and defines each component separately. In 16% of our cases followed longitudinally, we were able to identify tumor progression several months in advance of RANO criteria.
Conventional MRI lacks the precision to reliably differentiate tumor and nontumor areas with similar imaging patterns. However, the implications of such differentiation regarding treatment are substantial. Within the enhancing lesion area, recognition of dynamic, active tumor typically requires a change in therapy. Likewise, the identification of radiation necrosis or treatment-related changes can prevent unnecessary interventions. While salvage re-irradiation of brain tumors can typically be performed with a minimum of resultant treatment necrosis,23 re-irradiation to an area of preexisting treatment-related necrosis could exacerbate existing tissue damage and lead to potentially life-threatening neurologic toxicity.
Classification between tumor and nontumor areas within the nonenhancing lesion may provide new guidelines for radiation. Standards of treatment-planning techniques differ geographically.24 Accurate delineation of tumor area may either limit the size of the radiation field and decrease radiation exposure to nontumor brain tissue or define more accurately the optimal radiation field to include the entire tumor area and enable more efficient treatment.
In addition, as the use of immune-mediating therapies in glioma increases, the complexity of interpreting inflammatory responses to treatment (“pseudoprogression”) in enhancing and nonenhancing lesion areas will also increase. Misinterpretation of the response to a potentially helpful therapy could lead to premature cessation of a useful therapy and compromise patient outcome.25
Within the enhancing lesion area, all measured vascular parameters (except the BAT) clearly differentiated tumor and nontumor components. While both components demonstrated increased tissue permeability, the tumor component was characterized by much higher permeability and vp compared with the nontumor component and relative to NAWM. GB is characterized by a complex neovascularization process that results in formation of new, abnormal blood vessels and an impaired blood-brain barrier, and thus manifests as a hyperperfused area with increased blood flow, volume, and permeability, as was seen in our cohort study.26,27 The nontumor component demonstrated reduced vp compared with NAWM, consistent with data that have shown areas of treatment-related changes characterized by impaired blood brain barrier yet with reduced blood flow and volume, due to treatment-induced vascular endothelial damage and coagulative necrosis.5,28
Within the nonenhancing lesion area, despite a similar appearance between tumor and nontumor components on conventional imaging (FLAIR and T2WI), significant differences were found in all measured MRI parameters. The nonenhancing, nontumor component was characterized by reduced vp, in comparison with the NAWM and the nonenhancing tumor component. This pattern is consistent with previous reports of tissues proximal to brain tumors in patients10,29,30 and animal models31 and can be explained by compression of regional capillaries caused by vasogenic edema.32 The tumor component in our cohort was characterized by increased tissue permeability and increased perfusion as demonstrated by Ktrans and vp parameters, consistent with findings in previous studies.10,29,33
In this work, classification of lesion area was performed with a linear SVM classifier, a simple form of SVM, highly suitable for a limited training set size, for which the separating hyperplane is simply a plane in the feature space and the relative contribution of the features can be easily obtained. Various machine learning– and computational intelligence–based methods have been proposed for segmentation of brain tumors.34 In this study, SVM was found to be a robust, rapid method, tailored to the clinical data, easily implemented, and most important, it differentiated lesion components with high sensitivity and specificity.
The success of a supervised algorithm is directly determined by the selected training set. Thus, only clear-cut cases were used for this study while training the classifier, based on retrospective radiologic assessment for several months. Classification into tumor and nontumor components was performed based on conventional parameters (nT1WI+Gd, and nFLAIR) and parameters extracted from DCE (vp, Ktrans, and BAT). While perfusion parameters can be obtained using dynamic susceptibility contrast imaging and have been shown to differentiate between tumor and nontumor components,11⇓–13,35,36 DCE imaging is preferable due to its higher spatial resolution, less sensitivity to susceptibility artifacts, and provision of quantitative parameters, including permeability (Ktrans).
Several issues are important when interpreting the clinical relevance of classification results. The relatively small sample size and absence of histology-proved diagnosis in our cohort call for additional studies. Tumor components were found to be highly vascular while the nontumor components were found to have reduced vascularity. However, the enhancing, nontumor component may represent a more complex situation. It may reflect treatment-related changes (such as radiation necrosis) but may also represent tumor-associated hypoxic necrosis, which is one of the histologic hallmarks of GB.37 Thus, in such cases, it is important to consider the nontumor component in the context of active tumor surrounding the necrotic area.
Conclusions
The current study proposes a model of segmented RANO criteria, sRANO, that classifies tumor and nontumor components within a lesion area, with high sensitivity and specificity. Longitudinal assessment in patients demonstrated consistency between the classification results and radiologist's assessment in most cases; in 16% of patients, the segmentation results identified growth of highly vascular components and preceded the conventional radiologic diagnosis of tumor progression by several months. The proposed sRANO method and results presented in this study demonstrate the importance and contribution of segmentation of the enhancing and nonenhancing lesion areas into tumor and nontumor components, to improve therapy-response assessment of patients with malignant brain tumors.
Acknowledgments
We are grateful to Vicki Myers for editorial assistance, to Faina Vitenstein and Tuvia Ganot for their assistance in patient recruitment and performing MRI scans, and to Dr Ricardo Tarrasch for statistical consulting.
Footnotes
D.T. Blumenthal and M. Artzi contributed equally to this work.
Disclosures: Deborah T. Blumenthal—UNRELATED: Board Membership: Medical Advisory Board for Vascular Biogenics; Travel/Accommodations/Meeting Expenses Unrelated to Activities Listed: Roche, Bristol-Myers Squibb, Comments: travel grants to international meetings.* Money paid to the institution.
Paper previously presented in part at: Annual Meeting and Exhibition of the International Society for Magnetic Resonance in Medicine, May 7–13, 2016; Singapore.
Indicates open access to non-subscribers at www.ajnr.org
REFERENCES
- Received October 6, 2016.
- Accepted after revision January 3, 2017.
- © 2017 by American Journal of Neuroradiology
Jump to section
Related Articles
Cited By...
- Multi-Metric Quantitative MRI Identifies Spatially Distinct Age-Related Brain Differences in Female Bonnet Macaques
- Efficacy of intra-arterial carboplatin and bevacizumab in the C6 rat glioma model of glioblastoma multiforme
- Support Vector Machine for Lung Adenocarcinoma Staging Through Variant Pathways