Developing a Pipeline for Multiparametric MRI-Guided Radiation Therapy: Initial Results from a Phase II Clinical Trial in Newly Diagnosed Glioblastoma.

Quantitative mapping of hyperperfused and hypercellular regions of glioblastoma has been proposed to improve definition of tumor regions at risk for local recurrence following conventional radiation therapy. As the processing of the multiparametric dynamic contrast-enhanced (DCE-) and diffusion-weighted (DW-) magnetic resonance imaging (MRI) data for delineation of these subvolumes requires additional steps that go beyond the standard practices of target definition, we sought to devise a workflow to support the timely planning and treatment of patients. A phase II study implementing a multiparametric imaging biomarker for tumor hyperperfusion and hypercellularity consisting of DCE-MRI and high b-value DW-MRI to guide intensified (75 Gy/30 fractions) radiation therapy (RT) in patients with newly diagnosed glioblastoma was launched. In this report, the workflow and the initial imaging outcomes of the first 12 patients are described. Among all the first 12 patients, treatment was initiated within 6 weeks of surgery and within 2 weeks of simulation. On average, the combined hypercellular volume and high cerebral blood volume/tumor perfusion volume were 1.8 times smaller than the T1 gadolinium abnormality and 10 times smaller than the FLAIR abnormality. Hypercellular volume and high cerebral blood volume/tumor perfusion volume each identified largely distinct regions and showed 57% overlap with the enhancing abnormality, and minimal-to-no extension outside of the FLAIR. These results show the feasibility of implementing a workflow for multiparametric magnetic resonance-guided radiation therapy into clinical trials with a coordinated multidisciplinary team, and the unique and complementary tumor subregions identified by the combination of high b-value DW-MRI and DCE-MRI.

Introduction

Conventional therapies for glioblastoma (GBM) continue to rely on anatomic imaging modalities for both surgery and radiation therapy (RT), including T1 gadolinium- (T1-Gd) enhanced and T2-weighted fluid attenuated inversion recovery (T2-FLAIR) sequences that do not provide biological information about the underlying disease. Multiple studies have shown the prognostic value of physiological magnetic resonance imaging (MRI) techniques such as proton spectroscopy and perfusion and diffusion MRI and measures such as progression-free survival (PFS) and overall survival (OS) in predicting treatment response in patients with GBM (1–9). These imaging techniques may show abnormal tumor infiltration beyond the contrast-enhanced or nonenhanced areas conventionally targeted by surgery and radiation, and these may potentially be used to guide radiation treatment, reduce tumor recurrence, and improve patient outcome (10).

Dynamic contrast-enhanced (DCE)-MRI assesses relative cerebral blood volume (rCBV), cerebral blood flow, and vascular permeability, which are associated with neovascularization and tumor growth and predict PFS and OS in patients with GBM (1, 2, 5, 11). While regions of elevated rCBV often overlap with regions of contrast enhancement, the nonenhancing, infiltrating tumor beyond this region may potentially be underestimated with perfusion MRI (12). In contrast, diffusion-weighted (DW) MRI may identify tumor phenotype by estimating water mobility in the tissue microenvironment as an indicator of tumor cellularity (13). Apparent diffusion coefficient (ADC) is inversely correlated with cellularity but it may be unreliable in regions of highly cellular tumor, normal brain tissue, edema, and micronecrosis, yielding elevated ADC compared with normal tissue using standard b-values of 0–1000 s/mm2 (10). At our center, we developed a novel DW-MRI technique using high b-value (b = 3000 s/mm2) to selectively isolate solid, often nonenhancing, tumor that is predictive of PFS and often extends beyond the high-dose radiation target (14). We have shown that a combination of these imaging techniques (DCE-MRI and high b-value DW-MRI) into a multiparametric imaging signature predicts PFS with spatial correspondence with patterns of failure, representing biologically high-risk tumor subvolumes identifiable before therapy (12).

Based on these findings, we wished to develop a phase II study to evaluate the feasibility and efficacy of using a multiparametric hypervascular/hypercellular MRI signature to identify areas at highest risk of failure before radiation treatment in patients with newly diagnosed GBM (NCT02805179). Building on a prior phase I/II study showing the safety and efficacy of radiation dose-escalation with concurrent temozolomide (15), this multiparametric advanced imaging technique was used to prospectively guide the boost volume for dose-intensified radiation. To conduct this trial, the development of a workflow was required to permit the integration of an advanced, multiparametric imaging technique into the radiation treatment planning process. Here, we report the workflow and imaging characteristics of the initial patients treated on this prospective clinical trial.

Methodology

Patient Population

Adult patients of ≥18 years of age with newly diagnosed, pathologically confirmed suptratentorial GBM following any extent of resection were enrolled on this University of Michigan IRB-approved clinical trial following study-specific informed consent. Research was conducted in compliance with the World Medical Association Declaration of Helsinki-Ethical Principles for Medical Research Involving Human Subjects. Eligibility included Karnofsky performance status ≥70, minimal life expectancy of 12 weeks, adequate organ function, and maximal contiguous volume of tumor based on advanced imaging-defined boost volume of <1/3 volume of brain. Patients unable to undergo MRI scans or with prior overlapping radiation therapy were excluded. All patients were treated with standard concurrent daily (75 mg/m2) and adjuvant monthly (150–200 mg/m2) temozolomide.

MRI and Computed Tomography Simulation

All patients underwent an MRI simulation and computed tomography (CT) simulation after surgery for radiation planning, within 14 days of commencing chemoradiation. Rigid alignment of the T1-weighted contrast-enhanced and T2 FLAIR MRI to the CT image volumes in the Eclipse image registration workspace was performed by the medical physicist and verified by the radiation oncologist.

Commissioning of Hardware and Software QA

All MRI scans were performed on a 3 T scanner (Skyra, Siemens Healthineers, Erlangen, Germany) in the Radiation Oncology Department. Routine quality assurance of this scanner consists of daily checks of intensity uniformity as well as weekly checks following the ACR phantom accreditation scanning protocol (16). T1 mapping is a critical element of DCE-MRI analysis. To assess the accuracy, repeatability, and interplatform reproducibility of T1 quantification from variable flip angles, we scanned a National Institute of Standards and Technology (NIST) T1 water phantom on our system, provided by our participation in an NCI Quantitative Imaging Network (QIN) multicenter collaborative project (17). We used the extended Tofts model to quantify DCE-MRI, which was implemented in an in-house functional image analysis tool (imFIAT) (18). The performance of our implementation of the extended Tofts model was evaluated using digital reference objects, that is, synthesized DCE phantoms with and without noise, which was fully reported previously (19). In addition, we participated in an NCI QIN multicenter arterial input function (AIF) challenge to validate and compare our AIF delineation procedure to others' (20). On the basis of these evaluations and validations, imFIAT has been granted a level-2 benchmark by the NCI QIN (21).

MRI Acquisition

All images were acquired on a 3 T scanner (Skyra) using a 20-channel head coil. Conventional images, such as 2-dimensional (2D) T2-FLAIR images, and 3-dimensional pre- and postcontrast T1-weighted images, were acquired. In addition, physiological image acquisitions are described in the following subsections.

Diffusion-Weighted Imaging.

DW images were acquired using a 2D RESOLVE pulse sequence with diffusion weighting in 3 orthogonal directions and b-values of 0 and 3000 s/mm2 (1 and 4, respectively) to reduce geometric distortion required for radiation treatment planning. RESOLVE is a multishot technique that uses 2D navigator correction with readout-segmented echo planar imaging (22). Thirty slices were acquired to cover the whole brain with echo time (TE)/repetition time (TR) = 81/7650 milliseconds, matrix size = 160 × 160, and slice thickness/gap = 4.0/1.2 mm for ∼4.23 minutes. DW images acquired with b = 3000 s/mm2 were used for target definition. In addition, DW images were acquired by a 2D echo planar spin echo pulse sequence with diffusion weighting in 3 orthogonal directions and 11 b-values from 0 to 2500 s/mm2 as a backup scan. Thirty slices were acquired with TE/TR = 93/8200 milliseconds, matrix size = 192 × 192, slice thickness/gap = 4.0/1.2 mm, parallel imaging factor of 4 and a single average for 5 minutes. Parallel imaging factor of 4 was used to reduce the echo training time and thereby reduce geometric distortion. A full characterization of geometric accuracy of DW images with these acquisition parameters was previously reported (14).

DCE Imaging.

DCE images were acquired by a 3D gradient echo pulse sequence, called TWIST, in the sagittal orientation to avoid the in-flow effect and ensure artery coverage for an input function delineation. To cover the whole brain, a field of view of 250 × 256 × 187 mm3 was used with a matrix of 128 × 128 × 104 to obtain an isotropic voxel size of ∼1.9, which allows reformatting of the images in an axial plane or other planes as desired. Other acquisition parameters included flip angle = ∼10°, TE/TR = ∼0.95/2.65 milliseconds, temporal resolution = ∼3 s, dynamic phase volumes = 60, and total acquisition time = 3 minutes. Contrast was injected after 5 dynamic image volumes to achieve sufficient baseline data points.

Acquisition for T1 Quantification.

3D gradient echo images with 4 flip angles (3°, 7°, 12°, and 16°), TE/TR = 2.27/5.34 milliseconds, a voxel size = ∼2 mm, and total acquisition time = 1:45 minutes before contrast injection were acquired to quantify native T1. Low spatial resolution B1 maps were acquired to correct flip angle errors in T1 quantification, with an acquisition time of ∼12 seconds.

Target Volume Definition and Data Transfer

Physician-defined volumes were delineated in the Eclipse treatment planning system (Varian Medical Systems) directly on the MRI scans. T2/FLAIR abnormality was defined on the FLAIR MRI (FLAIR∧GTV). The surgical cavity (Cavity∧GTV), residual contrast enhancement (Gd∧GTV), and combination of cavity and contrast enhancement (GTV_Low) were delineated on the T1-Gd MRI. Volumes were then exported from the treatment planning system to the image analysis software (functional image analysis tool or imFIAT), for creation of DCE-MRI tumor volumes (high CBV [hCBV]), and high b-value DW-MRI tumor volumes (hypercellular volume [HCV]).

Image Analysis

DCE Analysis.

Three-parameter Tofts model was used to quantify the fractional plasma volume (Vp), transfer constant of contrast (Ktrans), and the fractional volume of extravascular extracellular space (ve) (23). The model was programed using C++ with a GUI in a functional image analysis tool (imFIAT). A full characterization of performance of software using digital reference objects with a large range of physiological parameters, acquisition parameters, and added Gaussian noise has been previously published (19).

In brief, we used the general assumption that where C is a contrast concentration in a voxel, and ΔR is a change in longitudinal relaxation rates after and before the contrast injection. If TR × R1 ≪1, which is generally satisfactory for brain normal tissue and tumors, where ΔS is a change in gradient echo intensities after and before the contrast injection, S is the averaged baseline gradient echo intensity before contrast injection, and R is the longitudinal relaxation rate before contrast injection. To obtain an AIF, 20 voxels with maximum intensity differences in a dynamic frame that was 1–2 time frames (∼4–7 seconds) before the enhancement peak were delineated. We participated in an NCI QIN challenge project of AIF delineation using our approach and our software (24). The parameters quantified from our AIF were well correlated with others (24).

T1 Calculation.

T1 maps were derived by fitting where α is a flip angle, TR is repetition time, and S is margination amplitude, to the 4-flip angle T1-weighted images using Simplex algorithm.

Using equations (1) and (2), AIF, T10 and the 3-parameter Tofts model, Vp maps were calculated. Then, hematocrit of 0.45 was used to convert Vp to CBV as CBV = [Vp/(1 − 0.45) × 100 (ml/100 g)], where blood density (1 mL/g) was used.

HCV Delineation

HCV was determined on DW images with b = 3000 s/mm2. A threshold was used from the normal tissue volume of interest (VOI) that was most contralateral to GBM. To obtain the normal brain VOI, an automated process was used to first extract the brain surface and find the middle line near the central fissure of the brain on T2-weighted images (b = 0). Then, the FLAIR abnormality volume was mirrored to the opposite side of the brain surface through the middle line. To remove CSF influence on the signals from the normal VOI, we remove all voxels with strong CSF signals by classifying CSF on T2-weighted images (b = 0) using fuzzy c-means. The VOI was eroded at least 5 mm from GTV∧FLAIR, and had ∼600 pixels per slice. Then, voxels within the GTV∧FLAIR on each slice were thresholded using mean + 2SD of the intensities in the VOI on the slice to account for DW intensity variations across slices. All these processes are fully automated. If a visual inspection of the normal brain VOI indicated the VOI inadequate, HCV could be recreated after adjusting the normal brain VOI by physician coauthors.

hCBV Delineation

hCBV was delineated on the CBV images with a threshold that was established in a previous study (12). Because normal white matter (WM) and gray matter (GM) have intrinsically different CBV values, the threshold value obtained from an uninvolved contralateral volume would vary depending on the ratio of GM to WM in that volume. We therefore segmented uninvolved contralateral GM in the frontal lobe (which has a higher CBV than uninvolved WM) and defined the hCBV tumor volume as the volume of tumor with CBV >1 SD above GM. This definition resulted in hCBV tumor volumes that predicted PFS and OS (12). Therefore, we used this definition in this pilot study of the clinical trial. This threshold was applied to GTV_Gd with 0–3 mm extension on CBV maps.

Volume Review

HCV and hCBV volumes were reviewed by the physician, neuroradiologist, and MRI physicist. Volumes were slightly edited during central review to remove components outside of the brain parenchyma for both HCV and hCBV, or overlap with blood vessels rather than parenchymal tumor volume for hCBV. Finalized tumor volumes were then imported from imFIAT back into the treatment planning system as binary image volumes associated with the HCV and hCBV image set. Images were automatically registered in the Eclipse Image Registration workspace to the original T1-post-Gd scan and checked by the clinical physicist. The clinical physicist then copied the HCV and hCBV volumes to the CT data set in the treatment planning system.

Image and Volume Registration and Delineation

The physician reviewed the HCV and hCBV volumes in the treatment planning system. CTV and PTV structures were then created as follows: CTV_Low was defined as a 1.7 cm expansion from GTV_Low delimited by normal anatomic boundaries. PTV_Low was defined 0.3 cm (0.2 cm positioning uncertainty with daily CBCT plus 0.1 cm MRI to CT registration uncertainty) as per institutional standard. For the advanced MRI (HCV and hCBV) boost target volumes, no CTV margin was used. PTV_High was defined as a 0.5 cm expansion from the HCV/hCBV volumes (0.2 cm positioning uncertainty + 0.2 cm RESOLVE DWI uncertainty + 0.1 cm MR to CT registration uncertainty).

PTV_Low was prescribed 60 Gy in 30 fractions and PTV_High was prescribed 75 Gy in 30 fractions using a simultaneous integrated boost technique. Volumetric modulated arc therapy using the Eclipse treatment planning system was used in all cases. The goal was to cover 100% of the target volumes with 95% of the prescribed dose, while maintaining conventional dose limits as utilized on cooperative group trials for high-grade glioma. This included maintaining optic chiasm and optic nerves ≤54 Gy, brainstem surface (ventral 3 mm of brainstem) ≤64 Gy, and brainstem core ≤55 Gy.

Results

Patient Characteristics and Imaging Subvolumes

The initial 12 patients enrolled between September 2016 and June 2017 were included in this analysis. Baseline characteristics of patients are described in Table 1. All patients had IDH1 wild-type tumors by immunohistochemistry. Fifty percent of patients underwent gross total resection, 33% underwent subtotal resection, and the remainder underwent biopsy alone. The workflow for image acquisition, volume delineation and data transfer, and treatment planning is depicted in Figure 1. All patients initiated radiation within 6 weeks of surgery and within 2 weeks of simulation. Advanced volume processing was generally done within 24–36 hours of simulation.

Table 1.

Baseline Patient Characteristics

Clinical Characteristic	No (%)
Median age (range)	65–(51-77)
Male	8 (67%)
Extent of resection
Gross total resection	6 (50%)
Subtotal resection	4 (33%)
Biopsy	2 (17%)
MGMT methylation status
Positive	3 (25%)
Negative	9 (75%)
Tumor location
Frontal lobe	4 (33%)
Temporal lobe	5 (42%)
Parietal lobe	2 (17%)
Occipital lobe	1 (8%)

Figure 1.

Integrated workflow diagram for the implementation of an advanced dynamic contrast-enhanced (DCE)- and high b-value magnetic resonance (MR) imaging signature to guide dose-intensified radiotherapy. TPS = treatment planning system; Gd = T1-Gd-enhanced MRI; GTV = gross tumor volume; CTV = clinical target volume; PTV = planning target volume.

Characteristics and distributions of CBV in normal frontal GM, normal WM, and the hCBV tumor volumes are shown in Figure 2. Note that the CBV value + 2SD in the frontal WM was below the mean value in the normal frontal GM, and thus, it cannot be used to define the elevated CBV in the tumor volume.

Figure 2.

The averaged CBV values in the hCBV tumor volumes, normal frontal white matter (WM) and normal frontal gray matter (GM) in the 12 study patients. The error bars depict the averaged standard deviations of cerebral blood volume (CBV) in the 3 volumes of interest across the 12 patients. Note that the mean CBV value + 2SDs in the frontal WM is smaller than the mean value in the frontal GM, and thus, it cannot be used as a threshold value to define the elevated CBV in the tumor volume.

Characteristics of conventional and advanced MRI target volumes are listed in Table 2. On average, the Gd-enhanced volume was >3 times larger than those of either HCV or hCBV, and 1.8 times larger than the union (combination) of HCV and hCBV. HCV and hCBV identified largely distinct volumes with only ∼1 cc of overlap between the 2 (range, 0.002–6.8 cc). The enhancing component of the union of HCV and hCBV was only 57% (range, 0.3–0.9) of the volume, with an average of 4.1 cc (range, 1.4–7.9) extending outside of the enhancing region. Only 2 patients showed minimal extension of the union HCV and hCBV volume beyond FLAIR (0.06 cc and 0.04 cc, respectively), and FLAIR volumes were ∼10 times larger than the union of HCV and hCBV. An example of HCV and hCBV volumes that are largely nonoverlapping overlaid on the corresponding T1-Gd MRI is shown in Figure 3. Two examples of representative radiation plans for 2 different patients are shown in Figure 4. For comparison, example plans without the advanced MRI boost are also depicted. As showed, the advanced MRI boost volume was generally a smaller tumor subregion contained within the conventional target volume.

Table 2.

Volume and Overlap of Conventional and Advanced Imaging Subvolumes

Target	Mean Volume (cc)	Range
GTV^Gd	23.9	3.9–49.9
GTV^FLAIR	128.9	39.2–248.5
GTV^HCV	7.5	1.7–20.4
GTV^hCBV	6.6	0.5–18.2
Union of HCV and hCBV	13.1	2.3–31.8
Overlap of HCV and hCBV	0.9	0.002–6.8
Overlap Gd and HCV	5.3	0.4–17.4
Overlap of Gd and hCBV	4.5	0.3–15.0
Overlap Gd and Union	8.9	0.9–26.0
Overlap FLAIR and HCV	7.5	1.7–20.4
Overlap FLAIR and hCBV	6.5	0.5–18.2
Overlap FLAIR and Union	13.1	2.3–31.7
HCV outside of Gd	2.2	0.8–3.6
hCBV outside of Gd	2.0	0.0–6.1
Union outside of Gd	4.1	1.4–7.9

GTV˄Gd = Gadolinium enhanced target volume; GTV^FLAIR = FLAIR target volume; GTV˄HCV = Hypercellular high b-value DW-MRI target volume; GTV^hCBV = Hyperperfused DCE-MRI target volume.

Figure 3.

An example of a patient with largely nonoverlapping hypercellular tumor regions (TVHCV) identified by high b-value diffusion-weighted (DW)-magnetic resonance imaging (MRI) (cyan, left panel) and hyperperfused tumor regions (TVHCBV) identified by DCE-MRI (red, middle panel). Significant extension of TVHCV is showed beyond the T1 Gd-enhanced region (overlay on T1 Gd-enhanced image, right panel).

Figure 4.

Representative images from radiation plans from 2 different patients. The top row depicts images of radiation plans using advanced MRI to boost tumor subregions to 75 Gy. High-risk tumor targets are identified by advanced MRI (cyan) beyond the abnormal regions seen on T1 Gd-enhanced conventional MRI (green). The conformal 75-Gy isodose line targeting the advanced imaging tumor volume is depicted in red, and the larger 60-Gy isodose line targeting the anatomic T1-Gd-enhanced region with standard clinical margins is depicted in gray-white. The bottom row depicts images of comparison standard radiation plans for the same patients based on anatomic T1-Gd-enhanced MRI with standard clinical margins prescribed to 60 Gy. As showed, advanced MRI-identified boost regions prescribed to 75 Gy were often contained within standard anatomic MRI regions prescribed to 60 Gy.

Discussion

While the limitations of anatomic MRI for radiation therapy have been reinforced by multiple studies showing that tumor identified by advanced MRI techniques extending outside of conventionally defined volumes predicts patient prognosis independent of T1-Gd, T2-FLAIR, and other clinical factors, advanced imaging techniques have not been incorporated into routine radiation planning (10, 25). In this initial report of a prospective, single-arm phase II trial for patients with newly diagnosed GBM from a single institution, we describe the end-to-end process of delivery of dose-intensified RT to predicted, high-risk tumor subregions identified by multiparametric MRI. The advanced hypercellular and hyperperfusion tumor subvolumes were significantly smaller than the conventionally defined T1-Gd and T2/FLAIR abnormalities standardly targeted for radiation treatment planning, and identified distinct regions that were frequently nonenhancing and therefore excluded from standard radiation boost volume definition. Using physician-defined volumes on conventional T1-Gd and T2-FLAIR images, the semiautomated creation of advanced MR boost volumes was accomplished for real-time planning, yielding successful delivery of advanced imaging-defined dose-intensified RT in all patients beginning within 2 weeks of simulation.

Given the known limitations of conventional MRI for defining tumor extent and predicting outcome in patients with GBM, the use of advanced imaging including perfusion and diffusion-weighted MRI has been studied for more than a decade to assess physiologic phenotypes of prognostic significance in this disease. Dynamic susceptibility contrast (DSC) and DCE-MRI permit quantitative estimation of parameters reflective of tumor neovascularization that has been associated with tumor growth in GBM, including CBV, cerebral blood flow, and Ktrans (23, 26). Maximum CBV and pathologically verified tumor vascularity are correlated, and elevated mean relative CBV (rCBV) >1.75 in patients with low- and high-grade gliomas is associated with shorter time to progression (2). Additionally, these perfusion parameters are predictive of overall survival in patients with malignant gliomas (10, 11).

Given the well-known limitation in geometric accuracy of DSC images, we selected T1-weighted DCE images to estimate CBV for the purpose of radiation boost target definition. Several previous studies as well as ours have directly compared CBV values estimated by DSC and DCE images (27–31). An early study reported good correlation of median values in tumors between 2 CBV estimates (r = 0.67) and excellent pixel-by-pixel correlation between the 2 estimates in normal brain tissues (r = 0.96) in 9 patients with intraaxial cerebral tumors (27). Another study of 32 patients with high-grade glioma showed a weak-but-significant correlation between the 2 estimates in the enhancing tumor volumes in a pixel-by-pixel comparison (28). Another study including 17 healthy subjects and 9 patients with glioblastoma reported excellent correlations of the 2 CBV estimates in normal GM and WM (r = 0.9 and r = 0.89, respectively) and a good and significant correlation in the tumor (r = 0.67) (29). A recent study examined the diagnostic accuracy of glioma grades in 26 patients using the median tumor values of the 2 CBV estimates, which achieved a similar diagnostic accuracy (30). We have performed similar analysis in 20 patients with brain metastases, in which both DCE and DSC images were acquired in a single session with 2 contrast injections. We found good correlation between the 2 CBV estimates in both normal brain tissue and brain tumor volumes (r = 0.66–0.71) (unpublished data). These similarities and discrepancies depend upon several factors. The 2 imaging methods rely on considerably different contrast mechanisms and model theories, which can be affected by different physical and physiological conditions. Also, different acquisition parameters and modeling implementations, for example, correcting T1 and vascular leakage effects in DSC analysis, can affect the results.

Additional physiologic properties of malignant gliomas may be assessed with DW-MRI, which has been used to assess the mobility of water molecules in the tissue microenvironment as a surrogate for tumor cellularity, and is a method for therapeutic response assessment as first shown in patients with glioma (6, 7, 13). An inverse correlation is observed between ADC and brain tumor cellularity in preclinical studies (13). However, known limitations of this approach for isolating tumor cellularity from normal brain tissue, edema, and micronecrosis in the heterogeneous GBM microenvironment may lead to unpredictably elevated ADC compared with normal tissue. This limitation may be mitigated through the use of high b-value DW-MRI (3000–4000 s/mm2) versus 0 and 800–1000 s/mm2 to attenuate signals due to edema (14, 32–34). We investigated the prognostic value of this approach, showing that the hypercellular tumor region (HCV) identified before RT using high b-value DW-MRI correlates with worse PFS in patients with newly diagnosed GBM treated with standard chemoradiation (14). We determined that in contrast to DCE-MRI alone, the combined use of high b-value MRI with DCE-MRI identifies largely spatially distinct regions with mean overall of only 21% (12). Moreover, the combination of these modalities correlated with patterns of failure and progression, and therefore, these are rationally targeted for intensified radiation treatment (12).

A limited number of studies are prospectively evaluating the incorporation of advanced imaging for the radiation treatment of patients with GBM. These include ongoing studies using proton MR spectroscopic imaging and amino acid positron emission tomography to guide radiation treatment in patients with GBM. Proton MR spectroscopic imaging detects chemical compounds reflective of cellular turnover and proliferation and correlates with histologic tumor cell density and survival in patients with GBM (35–38), although its use for radiation treatment has been limited to select centers with imaging expertise. Amino acid positron emission tomography including 11C-Methionine and 18F-radiolabeled 3,4-dihydroxy-6-[18F]fluoro-L-phenylalanine [18F]F-DOPA tracers is also under evaluation for targeting of potentially aggressive tumor regions beyond conventional MRI in ongoing trials in the United States and Europe. Studies have shown significant correlation in the standard uptake values of 11C-MET and 18F-FDOPA with nearly identical patterns of spatial uptake (39); both have been shown to be prognostic for survival and recurrence and potentially complementary to MRI, although not widely adopted and limited to research centers with expertise in complex radiotracer synthesis or on-site cyclotrons (40–42).

Our study represents the first report of the prospective implementation of a multiparametric imaging signature that is integrated in the RT workflow to guide intensified RT against distinct, poor prognosis phenotypes in patients with GBM. Initial implementation of the real-time use of a multiparametric MR signature for radiation planning involved QA and commissioning of DW- and DCE-MRI for clinical usage before clinical implementation. Implementation of an advanced MR biomarker for radiation treatment requires close coordination between the radiation oncologist, imaging physics, and clinical physics teams to delineate tumor volumes, process and transfer data between treatment planning and advanced imaging software, and ensure timely initiation of treatment. Limitations of this approach include the phenotypic and biologic diversity of GBM, and whether a multiparametric signature is sufficient to characterize this heterogeneity and guide treatment in this disease. To address this, we are acquiring and correlating other physiological imaging modalities with advanced MRI, as well as acquiring longitudinal imaging to evaluate whether temporal changes in advanced imaging features may be used to predict outcome and further tailor therapy.

In this report, we show the feasibility of the real-time use of a multiparametric MR signature to guide radiation treatment against prognostic, unique tumor subregions that substantially differ from the T1-Gd-enhancing high-risk boost volumes standardly defined by conventional MRI. Survival outcomes are awaited from this study, and future directions include translation of this workflow to a second site to validate the generalizability of this novel radiotherapeutic approach for patients with GBM.