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
(PRM
rFLAIR) 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 PRM
rFLAIR values. The calculated PRM
rFLAIR
values were categorized into three groups based on significant changes in
signal intensity: decrease (PRM
rFLAIR‑, blue voxels), increase (PRM
rFLAIR+,
red voxels) and unchanged (PRM
rFLAIR0, 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 PRM
rFLAIR 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 PRM
rFLAIR+
and PRM
rFLAIR‑ 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 PRM
rFLAIR patient samples were analyzed using the Kaplan‑Meier tumor progression curve (Figure 4). Figure 4 reveals the
metric PRM
rFLAIR+ 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 PRM
rFLAIR+,
the mean and median survival time for non‑responders
was estimated to be significantly less (p = 0.02) than for responders.
Similarly, by using PRM
rFLAIR-, 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 (PRM
rFLAIR+) were classified as non-responders with
significantly shorter OS and TTP versus responding patients with insufficient PRM
rFLAIR+
signal. Additionally, patients with a significantly decreased change in
voxel intensity (PRM
rFLAIR-), 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 PRM
rFLAIR
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.