Abstract
BACKGROUND AND PURPOSE: Accurate tumor grading is essential for treatment planning of pediatric brain tumors. We hypothesized that multiparametric analyses of a combination of permeability metrics and ADC histogram metrics would differentiate high- and low-grade tumors with high accuracy.
MATERIALS AND METHODS: DTI and dynamic contrast-enhanced MR imaging using T1-mapping with flip angles of 2°, 5°, 10°, and 15°, followed by a 0.1-mmol/kg body weight gadolinium-based bolus was performed on all patients in addition to standard MR imaging. Permeability data were processed and transfer constant from the blood plasma into the extracellular extravascular space, rate constant from the extracellular extravascular space back into blood plasma, extravascular extracellular volume fraction, and fractional blood plasma volume were calculated from 3D tumor volumes. Apparent diffusion coefficient histogram metrics were calculated for 3 separate tumor volumes derived from T2-FLAIR sequences, T1 contrast-enhanced sequences, and permeability maps, respectively.
RESULTS: Results from 41 patients (0.3–16.76 years of age; mean, 6.22 years) with newly diagnosed contrast-enhancing brain tumors (16 low-grade; 25 high-grade) were included in the institutional review board–approved retrospective analysis. Wilcoxon tests showed a higher transfer constant from blood plasma into extracellular extravascular space and rate constant from extracellular extravascular space back into blood plasma, and lower extracellular extravascular volume fraction (P < .001) in high-grade tumors. The mean ADCs of FLAIR and enhancing tumor volumes were significantly lower in high-grade tumors (P < .001). ROC analysis showed that a combination of extravascular volume fraction and mean ADC of FLAIR volume differentiated high- and low-grade tumors with high accuracy (area under receiver operating characteristic curve = 0.918).
CONCLUSIONS: ADC histogram metrics combined with permeability metrics differentiate low- and high-grade pediatric brain tumors with high accuracy.
ABBREVIATIONS:
- AROC
- area under receiver operating characteristic curve
- DCE-MRI
- dynamic contrast-enhanced MR imaging
- FL
- 3D tumor volume from T2-FLAIR images
- EES
- extracellular extravascular space
- Kep
- rate constant from EES back into blood plasma
- Ktrans
- transfer constant from blood plasma into EES
- PG
- 3D tumor volume from T1-contrast-enhanced images
- ROC
- receiver operating characteristic
- ve
- extracellular extravascular volume fraction
- vp
- fractional blood plasma volume
Primary brain tumors have an incidence rate of 5.67 per 100,0001 in children younger than 19 years of age and remain the most common cause of disease-related death in children. Accurate assessment of tumor grade is critical for diagnosis and treatment planning, especially with recent advances in genomic and molecular subtyping of tumors and targeted treatments for different subtypes. Conventional MR imaging remains the technique of choice for diagnosis and treatment; however, advanced imaging techniques such as perfusion and diffusion imaging have also proved useful for physiologic characterization of tumors.
Dynamic contrast-enhanced MR imaging (DCE-MRI) provides information on the vessel permeability and microvasculature in brain tumors and has been used extensively in adults in tumor grading,2⇓⇓–5 but less frequently in children.6,7 DTI provides information about tumor cellularity and is useful for characterization of tissue, tumor grading, tumor response to treatment, and distinction of tissue types.8 ADC histogram analysis provides a comprehensive analysis of the volumetric diffusion characteristics of the tumor and has been used to predict survival in both adults9,10 and children.11
We retrospectively analyzed DCE-MR imaging and DTI data from children presenting with a brain mass during 3 years at Boston Children's Hospital and examined the permeability and ADC histogram metrics as predictors of tumor grade. We hypothesized that a combination of DCE-MR imaging and ADC histogram metrics would improve the accuracy of tumor grading on MR imaging.
Materials and Methods
Subjects
Children who presented with a brain mass from 2013 to 2017 who had DCE-MR imaging and DTI studies were included. The study was performed with the approval of the institutional review board of the Dana Farber Cancer Institute.
MR Imaging Acquisition
All MR imaging studies were performed on a 3T scanner (Skyra; Siemens, Erlangen, Germany). Standard MR imaging in all patients consisted of sagittal T1, axial T2, precontrast axial T2-FLAIR, axial DTI (TR = 6500 ms, TE = 88 ms, slice thickness = 2 mm, b-values of 0 and 1000 s/mm2, 30 gradient directions), and postcontrast T1 images. All patients underwent a DCE-MR imaging protocol as follows: 1) variable flip angle echo-spoiled gradient-echo T1-mapping sequences using flip angles of 15°, 10°, 5°, and 2°; TR = 5 seconds; TE = minimum; FOV = 240 mm; section thickness = 5 mm; and 2) DCE-MR imaging sequence consisting of 50 phases, 7 seconds apart, flip angle = 15°, TR = 4 seconds, TE = minimum. FOV, section thickness, and scan locations were identical to those in the T1-mapping sequences. A single bolus of gadobutrol (0.1-mL/kg body weight) was injected 20 seconds after the start of scanning at an injection rate of 2 mL/s.
DCE-MR Imaging Postprocessing
DCE-MR images were transferred to a VersaVue workstation (In vivo, Gainesville, Florida) for automated processing using OmniLook software. Voxelwise T1 maps were generated from the variable flip angle images, and the 2-compartment Tofts model12 was used to generate voxelwise values of the pharmacokinetic parameters, namely the transfer constant from the blood plasma into the extracellular extravascular space (Ktrans), rate constant from extracellular extravascular space back into blood plasma (Kep), extracellular extravascular volume fraction (ve), and fractional blood plasma volume (vp).
The model of Weinmann et al13 was used for blood plasma concentration. ROIs were drawn on each section of the tumor around contrast-enhancing portions of the tumor by an imaging data analyst or by a PhD scientist and verified by a Certificate of Added Qualification–certified neuroradiologist, and the mean (over voxels) and SDs of each of the variables were recorded. Only enhancing tumors were included in the analysis.
ADC Histogram Analysis
Due to their superior anatomic contrast, axial B0 images from the DTI sequence were registered to the axial T2-FLAIR images using the mutual information algorithm in FSL (http://www.fmrib.ox.ac.uk/fsl)14 and the same transformation matrix used on the ADC maps to register them to the T2-FLAIR images. Axial postcontrast images were also separately registered to T2-FLAIR using the same methodology. 3D ROIs comprising the 3D tumor volume from T2-FLAIR images (FL) and the 3D tumor volume from T1-contrast-enhanced images (PG) were automatically generated using the thresholding feature in Fiji (http://fiji.sc/),15 an Open Source (https://opensource.org/) distribution of Java modules along with ImageJ software (National Institutes of Health, Bethesda, Maryland). The ADC values of the voxels in the 3D tumor volume from T2-FLAIR images (FL), and PG ROIs were then used to generate the FL and PG ADC volumes, respectively. These volumes were thresholded using a uniform range of 600–2600 × 10−6 mm2/s to automatically exclude cyst, necrosis, and hemorrhage. FL and PG ADC histograms were then generated using a uniform bin width of 7.8 mm2/s to facilitate comparison of shape descriptors of the histograms. Histogram metrics used for statistical analysis were the number of peaks (unimodal or bimodal), mean, SD, skewness, and kurtosis of these histograms. In the case of tumors showing bimodal ADC histograms, the lower peak was measured for mean and SD, and skewness and kurtosis were not recorded.
DCE-MR imaging and ADC histogram analyses were blinded to tumor grading results.
Statistical Methods
Low- and high-grade tumor groups were compared for permeability variables and ADC histogram metrics using the Wilcoxon test and t test using the NPAR1WAY and t test procedures of SAS (SAS Institute, Cary, North, Carolina).16 Corresponding to the t test, estimates and 95% confidence intervals for the difference between means were calculated.
Receiver operating characteristic (ROC) curves and their areas (AROC) were estimated for highly statistically significant variables in the permeability set (ve, Ktrans, Kep), ADC histogram set (FL_ADC_mean, PG_ADC_mean), and the 6 pairs of variables comprising 1 variable from each set. Focus was on the ve and FL_ADC_mean pair because these variables and their relationship were considered the most stable.
The ROC curve shows the trade-off between the true-positive rate (classifying the tumor of a patient with a high-grade tumor as high-grade) and false-positive rate (classifying the tumor of a patient with a low-grade tumor as high-grade). The AROC estimates the average true-positive rate over all possible false-positive rates and is also an estimate of the probability of correctly classifying a random pair of patients, 1 from the low-grade group and 1 from the high-grade group.
In addition to point estimates, Wald-based tests (corresponding to the null hypothesis that the true AROC is 0.5) and 95% confidence intervals for the AROC were calculated. All estimates, confidence intervals, and tests were nonparametric and were calculated using the LOGISTIC procedure of SAS.
Results
Sixty-three patients were identified. Six patients had final diagnoses that were not brain tumors, 10 had nonenhancing tumors and were excluded from DCE-MR imaging analysis, and 6 patients were excluded due to motion or susceptibility artifacts. The remaining 41 patients were included.
Twenty-three children were boys, and 18, girls (age range, 0.30–16.76 years; median age, 4.76 years; mean age, 6.22 years). There were 31 infratentorial and 10 supratentorial tumors, and all except 2 were newly diagnosed.
Forty patients had enhancing biopsy-confirmed tumors, and 1 had a diffuse midline glioma located in the pons presumed to be high-grade. World Health Organization grades 1–2 were considered low-grade, whereas grades 3–4 were considered high-grade. Of these tumors, 16 were confirmed low-grade (7 pilocytic astrocytomas, 5 low-grade gliomas, 1 mature teratoma, 1 atypical meningioma, 1 low-grade ganglioglioma, and 1 low-grade mixed germ cell tumor). Twenty-five tumors were high-grade (12 medulloblastomas, 4 glioblastomas, 4 anaplastic ependymomas, and 1 each of atypical teratoid/rhabdoid tumor, embryonal tumor not otherwise specified, choroid plexus carcinoma, embryonal tumor with rhabdoid features, and diffuse midline glioma). Not all final diagnoses were given as the World Health Organization 2016 classification because some cases did not fit well into this classification at the time of diagnosis.
Molecular subtyping was successful in 11 of the 12 patients with medulloblastoma, yielding 3 sonic hedgehog/SHH-activated, 2 wingless/WNT-activated, and 6 classified as group 3 or 4.
There was no statistically significant difference in mean ages between the low-grade (mean, 7.28 years; median, 5.31 years) and the high-grade cohort (mean, 5.54 years; median, 4.57 years) (P = .552).
Group comparisons of all permeability and ADC histogram variables are shown in Table 1. Of the pharmacokinetic variables analyzed, Ktrans and Kep were significantly higher in high-grade tumors compared with low-grade tumors (P < .001), and ve was significantly lower in the high-grade group (P < .001).
Group comparisons: low- versus high-grade tumor
Of the ADC histogram metrics analyzed, the mean ADC of the T2-FLAIR tumor volume (FL_ADC_mean) and the PG tumor volume (PG_ADC_mean) were significantly higher for low-grade tumors (Fig 1) compared with high-grade tumors (P < .001) (Fig 2). Low-grade tumors also showed a tendency for lower skewness for T2-FLAIR tumor volume (FL_ADC_skew, P = .006). None of the other variables showed statistically significant group differences between low- and high-grade tumors. Six high-grade tumors showed bimodal T2-FLAIR ADC histograms, and 4 showed bimodal PG ADC histograms, while all the ADC histograms for low-grade tumors had just a single peak.
A 6-year-old boy with a low-grade pilocytic astrocytoma. The upper row shows the T2-weighted image (A), T2-FLAIR image (B), postcontrast T1-weighted image (C), ADC map (D), and the ADC histogram of the FLAIR tumor volume (E), respectively. The lower row shows the Ktrans map (F), Kep map (G), ve map (H), vp map (I), and the ADC histogram of the enhancing tumor volume (J), respectively. The enhancing tumor volume was defined as the nodular enhancement of the solid component for pilocytic astrocytomas. Note the low Ktrans and Kep and the high ve values in this low-grade tumor. The mean ADC of the FLAIR tumor volume is 1821 × 10−6 mm2/s and 1808 × 10−6 mm2/s for the enhancing tumor volume.
A 5-year-old boy with parietal anaplastic ependymoma. The upper row shows the T2-weighted image (A), T2-FLAIR image (B), postcontrast T1-weighted image (C), ADC map (D), and the ADC histogram of the FLAIR tumor volume (E), respectively. The lower row shows the Ktrans map (F), Kep map (G), ve map (H), vp map (I), and the ADC histogram of the enhancing tumor volume (J), respectively. Note the high Ktrans and Kep and the lower ve values in this high-grade tumor. The histogram of the FLAIR tumor volume is bimodal. The mean ADC of the lower peak of the FLAIR tumor volume is 867 × 10−6 mm2/s and 1262 × 10−6 mm2/s for the enhancing tumor volume.
No group differences were found among the molecular subtypes in the medulloblastoma cohort (P ≥ .294).
ROC analysis of the 5 highly statistically significant variables (Ktrans, Kep, ve, FL_ADC_mean, and PG_ADC_mean) showed that all displayed high sensitivity and specificity in discriminating low- and high-grade tumors, with AROC values above 0.82 (P < .001) as seen in Table 2. AROC values for ve and FL_ADC_mean were 0.843 and 0.870, respectively, while the AROC estimate for the 2-variable combination of these metrics was 0.918 (Fig 3). Focus was on this pair because the individual variables and their relationship were considered the most stable.
Predictors of AROC
A, Scatterplot of ve versus FL_ADC_mean for all 41 cases analyzed. High-grade tumors are shown in red, and low-grade, in black. B, ROC curve for ve (AROC = 0.843). C, ROC curve for FL_ADC_mean (AROC = 0.870). D, ROC curve when ve and FL_ADC_mean are combined (AROC = 0.918). Note the improved accuracy shown for the multiparametric predictor, reflected in the increased AROC.
The Spearman correlation between ve and FL_ADC_mean across the 41 tumors was 0.49 (P = .001). Correlation within the low- and high-grade groups was not significant.
Discussion
Multiparametric imaging of pediatric brain tumors provides physiologic information beyond standard MR imaging. Koob et al7 recently reported increased accuracy in differentiating grades I and IV pediatric brain tumors using a combination of diffusion and DSC perfusion metrics. Similarly, Morana et al17 showed improved prognostic capability by combining DWI, arterial spin-labeling, and PET data in pediatric diffuse astrocytic tumors. Other studies have demonstrated the advantage of combinations of DWI, DSC perfusion, MR spectroscopy, and PET data in studying pediatric brain tumors.18⇓–20 To our knowledge, this is the first study to demonstrate a multiparametric approach using DCE and ADC histogram metrics in pediatric brain tumors.
While DSC and arterial spin-labeling perfusion metrics have been reported for pediatric brain tumor grading,21⇓–23 there are fewer studies using DCE-MR imaging.6 Our results show that among the DCE-MR imaging parameters, Ktrans, Kep, and ve all successfully distinguished low- and high-grade tumors, and this result is consistent with our report from a smaller cohort.6 High-grade tumors exhibited higher values of Ktrans and Kep, indicating leaky vessels. This finding is consistent with the angiogenesis seen in high-grade tumors caused by neoangiogenic proteins such as vascular endothelial growth factor. Our finding of lower ve in high-grade tumors differs from reports in adult glioblastomas,3,24,25 likely due to the multiple different high-grade tumor types in children. Jensen et al,26 however, found that adult patients with gliomas with lower ve in peritumoral edema had significantly shorter overall survival rates, supporting our hypothesis that lower ve corresponds to a higher grade and more aggressive tumor.
Of the ADC histogram metrics we tested, FL_ADC_mean and PG_ADC_mean showed significant differences between low- and high-grade tumors, with lower ADC corresponding to high-grade tumors. Lower ADC indicates high tumor cellularity and has been associated with shorter survival as well.9⇓–11 ADC histogram skewness and kurtosis did not show the same significant differences. Although the number of ADC histograms with bimodal peaks was too small for statistical analysis (6 for FL_ADC_mean and 4 for PG_ADC_mean), all tumors exhibiting bimodal peaks were high-grade. High-grade tumors often yield bimodal ADC distributions due to the presence of peritumoral edema.9⇓–11
We did not detect any group differences among the molecular subtypes of medulloblastoma, which may be due to the small sample size. Recent studies have reported other distinguishing MR imaging features of some of these subtypes,27⇓–29 and a larger future sample size may show differences in permeability and ADC histogram metrics among the subtypes. We would postulate that the group 3 tumors, which are ill-defined and have the poorest prognosis,30 would have high Ktrans and Kep and low ve with low mean ADC histogram values.
Because ve is a measure of the extracellular extravascular space and ADC is also indicative of extracellular space, we tested for a correlation between ADC and ve. ADC and ve were highly correlated across the entire cohort (P = .001). Although a previous study was unable to find voxelwise correlations between ve and ADC in adult gliomas,31 our findings suggest a link between the mean values in pediatric brain tumors, possibly due to the different tumor types seen in children.
Permeability metrics assess tumor microvasculature, whereas ADC histogram metrics assess tumor cellularity. We hypothesized that these modalities may provide complementary information to assist in tumor grading and tested combinations of permeability and ADC variables as predictors of tumor grade using ROC curves. A combination of ve and FL_ADC_mean showed a higher AROC (0.918) compared with ve (0.843) or FL_ADC_mean (0.870) alone. Thus, multiparametric analysis may yield the best predictors of tumor grade in pediatric brain tumors.
One limitation of the study is the relatively small number and heterogeneity of tumors studied. Larger numbers of individual tumor types are necessary to validate these findings across all tumor types. The addition of PET19 and MR spectroscopy7,18 has also shown promise and may benefit future analyses. Future prospective work could include other multiparametric approaches incorporating ADC histogram and T1 permeability to differentiate molecular and genomic subtypes of pediatric brain tumors.
Conclusions
The permeability metrics Ktrans, Kep, and ve and the ADC histogram metrics FL_ADC_mean and PG_ADC_mean were able to distinguish low- and high-grade pediatric brain tumors with a high degree of accuracy. A multiparametric approach combining permeability and ADC histogram metrics improves the accuracy of pediatric brain tumor grading with high sensitivity and specificity.
Footnotes
Disclosures: Tina Young Poussaint—UNRELATED: Grants/Grants Pending: Pediatric Brain Tumor Consortium Neuroimaging Center, Comments: National Institutes of Health grant with subcontract*; Royalties: Springer Science for a book. *Money paid to the institution.
Paper previously presented at: Annual Meeting of the American Society of Neuroradiology and the Foundation of the ASNR Symposium, April 22–27, 2017; Long Beach, California.
References
- Received August 18, 2017.
- Accepted after revision October 30, 2017.
- © 2018 by American Journal of Neuroradiology