Evaluation of FLAIR maps by PRM provides for glioma response assessment
Deborah Sharon Honrado Guest1, Craig Galbán1, Gary Luker1, Thomas Chenevert1, Benjamin Lemasson2, Robin Johannes Marius Navest3, Klaas Nicolaij3, and Brian Ross1

1Radiology, University of Michigan, Ann Arbor, MI, United States, 2Institut des Neurosciences, Université Grenoble Alpes, Grenoble, France, 3Department of Biomedical Engineering, Technische Universiteit Eindhoven, EINDHOVEN, Netherlands

Synopsis

This study investigates the possibility of adapting the PRM method for use with normalized FLAIR images to predict OS and TTP for indication of tumor recurrence. Glioma patients were separated into non-responders and responders to treatment. Voxels present in the union of the VOIs for the rFLAIR images were used to evaluate the PRMrFLAIR values and categorize patients into groups based on changes in signal intensity. This study shows that predicting TTP and OS is achievable using PRM with rFLAIR maps for patients treated with TMZ/IR and provides the first demonstration of quantifying FLAIR signals in patients over time.

Purpose

The Parametric Response Map (PRM) compares spatially paired voxels of co‑registered images at different time‑points to quantitatively analyze significant changes in the histology of the tumor over time, which provided an opportunity to quantify changes in tissue properties in response to therapy. Studies show that the PRM metric applied to apparent diffusion coefficient (ADC) maps are predictive of clinical progression for head and neck squamous cell carcinoma 1 and one-year survival of glioma patients 2–5, whereas percentage‑change in mean ADC failed to predict outcome. Thus, temporal dynamics within heterogeneous tumors are often lost in scalar metrics 6. The PRM method has been extended across modalities and serves as a generalizable approach for quantifying and spatially visualizing tissue changes longitudinally 7,8. The purpose of this study was to investigate the possibility of adapting the PRM method for use with normalized axial fluid‑attenuated inversion recovery (FLAIR) images (PRMrFLAIR) to predict overall survival (OS) and time to progression (TTP), thus functioning as an early indicator of tumor recurrence.

Methods

Patients (n=52) with grade 3/4 glioma were recruited for a prospective single-center clinical trial. Patients were separated into groups; those surviving less than one year (non-responders) and those surviving at least one year (responders) after treatment. Each patient was treated with a total radiation dose of 60‑70 Gy over 30 fractions. Radiotherapy is typically combined with concurrent TMZ, followed by 6 months of adjuvant chemotherapy. FLAIR images were acquired post-surgical resection, prior to treatment (0‑week image) and 10 weeks after commencement of treatment. MRI parameters can be found in Figure 5. The 10‑week FLAIR images were then co‑registered to the 0-week images using Elastix 9 (Figure 1). For both time‑points, volumes of interest (VOI) were used to outline the tumor volume. Both FLAIR images were normalized using the cerebellum to obtain a relative FLAIR (rFLAIR) image, ensuring signal intensities were comparable between time‑points. Voxels present in the union of the two VOIs for the rFLAIR images were used to evaluate the PRMrFLAIR values. The calculated PRMrFLAIR values were categorized into three groups based on significant changes in signal intensity: decrease (PRMrFLAIR‑, blue voxels), increase (PRMrFLAIR+, red voxels) and unchanged (PRMrFLAIR0, green voxels) (Figure 1). Significant changes were determined using a 95% CI threshold calculated from a linear least squares regression over contralateral healthy brain tissue for ten randomly selected patients. The mean voxel intensity within the tumor VOIs were calculated for both 0 and 10‑week rFLAIR images. Patients were categorized into the PRMrFLAIR groups (red, blue and green voxels) using receiver operating characteristic curve analysis to calculate the optimal signal cut‑off (Figure 3). Within each cut‑off, a Kaplan‑Meier estimate was used to evaluate the survival probability, after treatment, of both responders and non‑responders and significance was tested using a Log Rank (Mantel‑Cox) test.

Results

The optimal cut‑off for voxel signal intensities for both the PRMrFLAIR+ and PRMrFLAIR‑ to determine one-year survival is shown in Figure 3. From this cut-off, patients were characterized as either responders or non‑responders. The categorized PRMrFLAIR patient samples were analyzed using the Kaplan‑Meier tumor progression curve (Figure 4). Figure 4 reveals the metric PRMrFLAIR+ allows for predicting mean and median TTP for non‑survivors which was significantly less (p < 0.001) than the estimated mean and median TTP for survivors. OS was also analyzed using Kaplan‑Meier estimation. Using the metric PRMrFLAIR+, the mean and median survival time for non‑responders was estimated to be significantly less (p = 0.02) than for responders. Similarly, by using PRMrFLAIR-, the predicted mean and median survival time for non-responders was also significantly less (p = 0.008) than the mean and median survival time for responders.

Discussion

Results show patients who had a significantly increased change in voxel signal intensity (PRMrFLAIR+) were classified as non-responders with significantly shorter OS and TTP versus responding patients with insufficient PRMrFLAIR+ signal. Additionally, patients with a significantly decreased change in voxel intensity (PRMrFLAIR-), had a significantly longer OS. This study shows that predicting TTP and OS is achievable using PRM with rFLAIR maps for glioma patients treated with TMZ/IR.

Conclusion

FLAIR images are routinely acquired for glioma patient clinical diagnosis and management. We demonstrated that these FLAIR images could be normalized to generate maps that could be used to provide quantitative PRM maps. This study provides the first demonstration of spatially following and quantifying FLAIR signals in patients over time and opens up new possible diagnostic uses for PRMrFLAIR in the clinic.

Acknowledgements

This work was supported by the US National Institutes of Health grant P01CA085878.

References

1. Galbán CJ, Mukherji SK, Chenevert TL, et al. A feasibility study of parametric response map analysis of diffusion-weighted magnetic resonance imaging scans of head and neck cancer patients for providing early detection of therapeutic efficacy. Transl Oncol. 2009;2(3):184-190.

2. Ellingson BM, Malkin MG, Rand SD, et al. Validation of functional diffusion maps (fDMs) as a biomarker for human glioma cellularity. J Magn Reson Imaging. 2010;31(3):538-548.

3. Moffat B a, Chenevert TL, Lawrence TS, et al. Functional diffusion map: a noninvasive MRI biomarker for early stratification of clinical brain tumor response. Proc Natl Acad Sci U S A. 2005;102(15):5524-5529.

4. Lemasson B, Chenevert TL, Lawrence TS, et al. Impact of perfusion map analysis on early survival prediction accuracy in glioma patients. Transl Oncol. 2013;6(6):766-774.

5. Moffat B a, Chenevert TL, Meyer CR, et al. The functional diffusion map: an imaging biomarker for the early prediction of cancer treatment outcome. Neoplasia. 2006;8(4):259-267.

6. Boes JL, Hoff B a, Hylton N, et al. Image registration for quantitative parametric response mapping of cancer treatment response. Transl Oncol. 2014;7(1):101-110.

7. Galbán CJ, Chenevert TL, Meyer CR, et al. The parametric response map is an imaging biomarker for early cancer treatment outcome. Nat Med. 2009;15(5):572-576.

8. Galbán CJ, Han MK, Boes JL, et al. Computed tomography–based biomarker provides unique signature for diagnosis of COPD phenotypes and disease progression. Nat Med. 2012;18(11):1711-1716.

9. Klein S, Staring M, Murphy K, Viergever MA, Pluim JPW. Elastix: A toolbox for intensity-based medical image registration. IEEE Trans Med Imaging. 2010;29(1):196-205.

Figures

Schematic of the PRM process. 10-week FLAIR images were warped to the 0-week images, then both images were normalized to obtain a rFLAIR image. Voxels in the union of the VOIs were used to evaluate the PRMrFLAIR values to categorized the patients into PRMrFLAIR‑, PRMrFLAIR+ or PRMrFLAIR0.

ROC curve for PRMrFLAIR. Voxel signal threshold to determine one‑year survival is: PRMrFLAIR+ = 20%, AUC = 0.74±0.08, p = 0.02. PRMrFLAIR- = 26%, AUC = 0.69±0.08, p = 0.03. No significance for PRMrFLAIR0 or D% in voxel intensity between the 0 and 10‑week rFLAIR images.

Mean and median OS and TTP for survivors and non‑survivors using PRMrFLAIR+ and PRMrFLAIR‑ metrics

MRI parameters. Images were acquired using a 1.5T Signa system (General Electric, USA) or an Achieva 3T system (Philips, the Netherlands)



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
0129