About the Authors:
Josep Puig
* E-mail: [email protected]
Affiliation: Imaging Research Unit, Institut de Diagnostic per la Imatge (IDI), Girona Biomedical Research Institute (IDIBGI), Hospital Universitari Dr Josep Trueta, Girona, Spain
Javier Sánchez-González
Affiliation: Philips Healthcare Ibérica, Madrid, Spain
Gerard Blasco
Affiliation: Imaging Research Unit, Institut de Diagnostic per la Imatge (IDI), Girona Biomedical Research Institute (IDIBGI), Hospital Universitari Dr Josep Trueta, Girona, Spain
Pepus Daunis-i-Estadella
Affiliation: Department of Computer Science, Applied Mathematics, and Statistics, University of Girona, Girona, Spain
Christian Federau
Affiliations Department of Diagnostic and Interventional Radiology, Centre Hospitalier Universitaire Vaudois (CHUV) and University of Lausanne, Lausanne, Switzerland, Department of Radiology, Neuroradiology Division, Stanford University, Palo Alto, United States of America
Ángel Alberich-Bayarri
Affiliation: Biomedical Imaging Research Group (GIBI230), La Fe Polytechnics and University Hospital, La Fe Health Research Institute, Valencia, Spain
Carles Biarnes
Affiliation: Imaging Research Unit, Institut de Diagnostic per la Imatge (IDI), Girona Biomedical Research Institute (IDIBGI), Hospital Universitari Dr Josep Trueta, Girona, Spain
Kambiz Nael
Affiliation: Department of Radiology, Icahn School of Medicine at Mount Sinai, New York, United States of America
Marco Essig
Affiliation: Department of Radiology, University of Manitoba, Winnipeg, Canada
Rajan Jain
Affiliation: Division of Neuroradiology, Department of Radiology, NYU Langone Medical Center, New York, United States of America
Max Wintermark
Affiliation: Department of Radiology, Neuroradiology Division, Stanford University, Palo Alto, United States of America
Salvador Pedraza
Affiliation: Imaging Research Unit, Institut de Diagnostic per la Imatge (IDI), Girona Biomedical Research Institute (IDIBGI), Hospital Universitari Dr Josep Trueta, Girona, Spain
Introduction
High-grade neoplasms produce a complex microvascular network to satisfy a growing need for nutriments and oxygen [1], and glioblastoma is among the most angiogenic tumors [2]. Cerebral blood volume (CBV) correlates with the degree of angiogenesis and increased local perfusion correlates with tumor grading and survival [3–5]. Therefore, hemodynamic parameters influenced by vessel density and perfusion of the microvasculature, such as CBV and cerebral blood flow (CBF), can be used as surrogate biomarkers [6–8]. As these perfusion parameters can be measured by dynamic susceptibility contrast (DSC) MRI, this technique plays an important role in the baseline evaluation and follow-up of brain tumors. Recently, intravoxel incoherent motion (IVIM) has been proposed as an alternative perfusion MRI technique [9–18]. IVIM uses endogenous tracers to obtain perfusion-related indexes from diffusion-weighted imaging (DWI) datasets without contrast agents.
Considering the vascular bed as a random network of vessels where blood flows freely, Le Bihan et al. [19] demonstrated that IVIM could distinguish between water diffusion and the microcirculation of blood in the capillary network. In biological tissues, diffusion and perfusion are physically different phenomena, and the incoherent motion of spins, which can be understood as the spatial “mixing” of spins during the time of measurement in each image voxel, bi-exponentially reduces the signal amplitude observed when different diffusion b-values are applied [20]. In other words, DWI is also sensitive to perfusion because the flow of blood in randomly oriented capillaries mimics a diffusion process. IVIM modeling allows the extraction of two diffusion coefficients, one related to molecular diffusion restriction, called the diffusion coefficient (D), and another related to movements of blood in the microvasculature, called the pseudodiffusion coefficient (D*). A third parameter, the perfusion fraction (f), describes the fraction of incoherent signal arising from the vascular compartment in each voxel. In recent years, advances in MR hardware have allowed short-time acquisitions with multiple b-values and sufficient signal-to-noise ratio, reviving interest in IVIM for imaging tumors in the brain [8–12] and in body tissues where vascularity is important [18,21–25] for characterizing tumors and predicting or monitoring the response to treatment [9,12,26,27]. Recent evidence suggests that f values can help differentiate between low- and high-grade gliomas [10–12], and f correlates moderately with DSC rCBV [11]. In addition, in rodent models of glioma, f correlates positively with vessel density at histology [28].
To our knowledge, no data about the usefulness of IVIM-metrics to predict survival in glioblastoma patients are available. Therefore, we determined whether IVIM-metrics D, D*, and f are useful in predicting tumor response to treatment and survival in newly diagnosed glioblastoma, by analyzing them in contrast-enhancing regions (CER) and non-enhancing regions (NCER) surrounding the CER.
Materials and Methods
Patient characteristics
The ethics committee at Hospital Dr Josep Trueta approved this prospective study, and all patients provided written informed consent to participate in it. From November 2009 to March 2011, we enrolled 23 consecutive patients with newly diagnosed glioblastoma suspected on CT or MRI. After imaging, all lesions were biopsied. Eight (34.78%) patients were excluded from the study: five because histology ruled out glioblastoma and three because motion artifacts hindered image evaluation. Therefore, 15 patients (7 women; mean age, 66±11 years) were included. Patients were managed according to published guidelines [29]. The combination of surgery, radiotherapy, and chemotherapy with concomitant and adjuvant temozolomide was considered standard treatment. Patients did not receive corticosteroids before MRI. Survival was measured from the pretreatment MRI study to death.
Conventional MRI
MRI was performed on a 1.5-T MR scanner (Gyroscan Intera 1.5T Master; Philips Healthcare, Best, the Netherlands) using an eight-channel head coil. Before contrast administration, we acquired axial T1WI SE (TR536ms, TE15ms), axial T2WI fast SE (TR4400ms, TE110ms), and axial FLAIR (TR8000ms, TE115ms, TI2200ms) sequences. We used a 230-mm field of view, 5-mm section thickness, and 256x192 matrix for these sequences. Five minutes after gadobutrol injection, we obtained axial T1WI SE (TR600ms, TE10ms) images parallel to the bicommissural line.
IVIM MRI
24 axial DW images (TR3000/TE76ms) were obtained using single-shot spin-echo echo-planar imaging (EPI) before contrast-enhanced MRI. The EPI factor was 41, and the sensitivity-encoding factor was two. We used a 200-mm field of view, 5-mm section thickness, and 96 x 77 matrix. The measured pixel size was 2.4x2.9x5mm. We used 13 b-values: 0, 10, 20, 30, 50, 100, 150, 200, 350, 500, 650, 800, and 1000 sec/mm2. The total acquisition time was 3 minutes 48 seconds per patient.
DSC-MRI
Using the same section orientations used for DW images, we acquired dynamic T2*-weighted gradient-echo echo-planar images (TR1800ms; TE25ms) during the first pass of a standard dose (0.1 mmol/kg) bolus of gadobutrol injected via an antecubital vein at 5 ml/s followed by 30 ml saline solution. To assure that steady-state magnetization was reached, a five series of dummy scans (i.e., the pulse sequence is run, but data are not acquired) were inserted immediately before the start of each perfusion series consisted of 50 dynamic acquisitions. Between 8 and 10 per-Gd baseline images were acquired. Based on T2WI and FLAIR images, we selected seven to ten sections through the tumor for PWI in a single TR with an in-plane resolution of 1.95x1.95 mm2 and slice thickness of 7 mm. The methods used for acquiring data and the algorithm for calculating CBV-corrected maps for contrast agent extravasation are described elsewhere [7].
IVIM Image Processing
The IVIM model considers that two compartments exist in biologic tissue: a slow-moving compartment, where particles diffuse in a Brownian fashion as a consequence of thermal energy, and a fast-moving compartment (the vascular compartment), where water molecules moves as a consequence of forced blood circulation [19]. In the vascular compartment, D* describes the displacement of blood on a macroscopic level in an assumed randomly laid vascular network. For the perfusion to be physiologically meaningful, D* must be greater than D. Therefore:(1)where f is the fraction of the diffusion linked to microcirculation (perfusion fraction), D is the diffusion parameter representing pure molecular diffusion (diffusion coefficient), and D* is the diffusion-related incoherent microcirculation; S0 is the signal intensity at a b-value of zero (i.e., without diffusion weighting), and Sb is the signal intensity for each b-value (i.e., at each diffusion gradient).
Regions of Interest
A neuroradiologist with 20 years’ experience (S.P.) reviewed the anatomical images, using Olea Sphere V.2.0 software (Olea Medical, La Ciotat, France). NCER was defined as the hyperintense area surrounding the CER on FLAIR [30]. A fully automated deconvolution analysis was performed to create parametric images of CBV and CBF [31] in the MR Extended Workspace (Philips Healthcare, Best, the Netherlands). Due to the technique’s limitations in obtaining absolute CBF values, an extra ROI was placed in healthy gray matter as a reference [32]. To scale all CBF values, gray matter flow was established at the same level (65ml/100g/min) for all patients [33]. DW data were registered to the image with b-value = 0 s/mm2 using an affine transformation and a mutual information algorithm to avoid image distortion due to eddy currents. Images were analyzed with a computer program developed within the research group on the Philips Research Imaging Development Environment research platform using Interactive Data Language 6.3 (Research Systems Inc.; Boulder, CO, USA); this program fits every pixel to the three parameters in the model described in eq 1 using a Levenberg-Marquardt least-squares minimization algorithm [34]. The software generates three IVIM maps (one for D, one for D*, and one for f) in about 3 minutes. A reader (G.B. with 11 years’ experience) manually placed ROIs in the CER and NCER for each tumor, with maximal f on three contiguous axial sections. Mean ROI size was 32±12 mm2. Large vessels and cystic or necrotic tumor areas were excluded. The corresponding ROIs were then drawn on the DSC images (mean size, 30±14mm2) and the results were averaged for CER and NCER. To enable intraobserver reliabilities to be calculate, all measurements were repeated 1 month after the first determination. The observer was blinded to the clinical and outcome data of the patients.
Statistical Analysis
Means and standard deviations were calculated for all parameters. Data were evaluated through the significance of the Pearson product-moment correlation coefficient. Linear regression analysis was performed using rCBV and rCBF values from DSC-MRI and the f, D, and D* values from IVIM. Receiver operating characteristic analysis was used to determine the optimal perfusion MRI parameter cutoffs for predicting 6-month survival. Prognostic factors included age, sex, Karnofsky Performance Score, treatment, volume of CER, volume of NCER, IVIM-metrics, and DSC-MRI parameters. Survival curves were calculated using the Kaplan-Meier method. We used the multivariate Cox proportional hazards model to adjust for the influence of prognostic factors. We used intraclass correlation coefficients (ICC) to compare measurements of rCBV, rCBF, f, D, and D* in CER and NCER, classifying intraobserver reliability as fair (ICC = 0.5–0.7), good (0.7–0.9), or almost perfect (>0.90). We also analyzed the variability of the measurements by Bland-Altman plots showing the mean difference between two methods of measurement, and 95% limits of agreement as the mean difference [35]. Minitab version 16.2.1 was used for statistical analyses (Minitab Inc.;State College,PA,USA). Significance was set at P<0.05.
Results
Patient Data
Table 1 summarizes patients’ clinical and imaging characteristics. All 15 patients died during the observation period. Survival was 10.6±6.23 months (range, 5–21 m) in patients receiving standard treatment and 4.7±2.81 months (range,1–8.5 m) in patients not receiving standard treatment.
[Figure omitted. See PDF.]
Table 1. Patient characteristics.
https://doi.org/10.1371/journal.pone.0158887.t001
Associations and correlations between DSC-MRI parameters and IVIM-metrics
Table 2 shows the associations for the values of DWI parameters, DSC-MRI parameters, and IVIM-metrics and the correlations between them for CER and NCER. The IVIM-metric fCER correlated moderately with rCBFCER and rCBVCER (R = 0.65;P = 0.01 and R = 0.49;P = 0.04, respectively) (Fig 1) and with DCER and D*CER (R = -0.68;P = 0.01 and R = 0.80;P<0.001, respectively). D*CER also correlated with rCBFCER (R = 0.71;P<0.001). rCBFCER significantly correlated with rCBVCER (R = 0.82;P<0.001). Moreover, fCER correlated with D*NCER (R = 0.85;P = 0.04) and with rCBFNCER (R = 0.56;P = 0.03) and rCBVNCER (R = 0.56;P = 0.03).
[Figure omitted. See PDF.]
Table 2. Correlations between IVIM-metrics and DSC-MRI parameters for CER and NCERa.
https://doi.org/10.1371/journal.pone.0158887.t002
[Figure omitted. See PDF.]
Fig 1. Glioblastoma in a 62-year-old woman.
(A) Axial FLAIR image. (B) Contrast-enhanced T1-weighted image shows a rim-enhancing tumor. (C) Cerebral blood flow map obtained with DSC-MRI shows hyperperfusion signals predominantly in the left margin of the tumor (arrows). (D) DSC cerebral blood volume map. The hyperperfusion signal intensities correlate with those seen on cerebral blood flow map. (E) f map clearly highlights the area with high perfusion in the margins of the tumor (arrows), which is more evident than in C and D. (F) D map shows small restricted diffusion area (average D value = 0.895 x 10−3 mm2/s) predominantly in the anterior tumor margin (arrows). (G) D* map shows increased fast-diffusion values in the tumor tissue (arrows).
https://doi.org/10.1371/journal.pone.0158887.g001
Survival analysis
Table 3 compares the IVIM-metrics and DSC-MRI parameters for patients who survived more than 6 months versus those who survived less than 6 months. Values for fCER, D*CER, rCBFCER, and rCBVCER were significantly higher in patients who survived less than 6 months (Fig 2). The proportion of patients that received standard treatment was higher in the group that survived more than 6 months. Table 4 shows the cutoff values for fCER, D*CER, rCBFCER, and rCBVCER for predicting survival. The cutoff fCER>9.86% had the highest AUC for predicting 6-month survival (100% sensitivity, 71.4% specificity, 100% positive predictive value (PPV), and 80% negative predictive value (NPV); AUC 0.893). The cutoff D*CER>21.712x10−3mm2/s yielded 100% sensitivity, 71.4% specificity, 80% PPV, and 100% NPV with an AUC of 0.857. The treatment was a significant predictor for 6-month survival (85.7% sensitivity, 57.1% specificity, 70% PPV and 80% NPV; AUC 0.723). In the Cox regression analysis, treatment was the most important factor (hazard ratio 5.484, 95% confidence interval 1.162–25.88, P = 0.031). In the multivariate analysis, only fCER combined with treatment predicted survival 100% (Table 4). Survival rate was significantly shorter in patients with high values of fCER (P = 0.008) and D*CER (P = 0.007) independently of the treatment received (Fig 3).
[Figure omitted. See PDF.]
Table 3. Clinical data, and diffusion and perfusion parameters in contrast-enhancing and non-enhancing regions according to survivala.
https://doi.org/10.1371/journal.pone.0158887.t003
[Figure omitted. See PDF.]
Table 4. Survival prediction: summary of class performance and hazard ratios for associations between imaging features and overall survivala.
https://doi.org/10.1371/journal.pone.0158887.t004
[Figure omitted. See PDF.]
Fig 2. Boxplots of DSC-MRI parameters and IVIM-metrics for CER according 6-month survival.
https://doi.org/10.1371/journal.pone.0158887.g002
[Figure omitted. See PDF.]
Fig 3. Kaplan-Meier survival curves comparing survival rates for treatment and for pre-specified cutoff values of fCER and D*CER (upper row) and for these cutoffs according to treatment received (lower row) Surgery, radiotherapy, and chemotherapy with concomitant and adjuvant temozolomide was considered standard treatment (tx).
https://doi.org/10.1371/journal.pone.0158887.g003
Intraobserver Reliability
Intraobserver agreements were almost perfect for DSC-MRI CER indexes (ICC = 0.916), DSC-MRI NCER indexes (ICC = 0.949), fCER (ICC = 0.979), fNCER (ICC = 0.983), DCER (ICC = 0.930), and DNCER (ICC = 0.98); intraobserver agreement was good for DSC gray matter indexes (ICC = 0.731). The Bland-Altman plots confirmed the high intraobserver reliability (S1 Fig).
Discussion
Our preliminary results are an important step in collecting evidence about the feasibility and usefulness of IVIM imaging as a quantitative method to measure perfusion in glioblastoma. To our knowledge, this is the first study to demonstrate the usefulness of IVIM-metrics in predicting survival in patients with newly diagnosed glioblastoma. We found that patients with increased fCER and D*CER had significantly shorter survival independently of the treatment they received. Although a detailed analysis will require more data, the following can be deduced from this small cohort: the cutoffs fCER = 9.86% and D*CER = 21.712 x10−3mm2/s on pretreatment MRI yielded the highest predictive power for 6-month survival (AUC 0.893 and 0.857, respectively). However, in the Cox regression models, treatment was the only significant variable (P = 0.031), although D*CER was nearly significant (P = 0.068).
The standard of care for newly diagnosed glioblastoma is now maximum safe surgical resection followed by radiotherapy plus concomitant and adjuvant chemotherapy with temozolomide [36]. This approach is based on a landmark phase III trial that reported median survival after surgery of 14.6 months in patients randomized to receive radiotherapy plus temozolomide compared to 12.1 months in those that receive radiotherapy alone [37], and other studies have corroborated this survival benefit [36,38]. Our results are consistent with these reports; importantly, however, we also found that adding fCER data to treatment data enable survival to be predicted with an accuracy of 100%. Our preliminary results suggest that patients treated with the standard of care who had fCER or D*CER values over a pre-specified cutoff had worse survival than those who had fCER or D*CER values below the cutoff. Therefore, IVIM-metrics may help tailor the therapeutic approach in upcoming studies.
We found a negative correlation between fCER and DCER, probably because regions with highest tumor cellularity almost certainly correspond to regions with highest vascularity. Bisdas et al. [12] revealed that IVIM fitting of the diffusion data allowed the contribution of perfusion to be separated from the contribution of true diffusion, thus providing better information than the apparent diffusion coefficient (ADC) for discriminating between low- and high-grade gliomas. Although one study found lower ADC values in high-grade gliomas than in lower-grade gliomas [39], another study reported considerable overlap in ADC values between low- and high-grade gliomas [40].
We found a strong positive correlation between fCER and D*CER and moderate positive correlations between fCER and rCBFCER and rCBVCER. The correlation with rCBFCER was stronger than the correlation with rCBVCER because rCBV is sensitive to vessel wall permeability, whereas f reflects only blood flowing in small vessels. Our results are in line with those recently reported by Federau et al. [11] in 21 gliomas (16 high-grade and 5 low-grade). They found that f correlated moderately with rCBV (r = 0.59) and, in the regions of maximum f, was significantly higher in the high-grade group.
In a recent study, Iima et al. [28] used a 17.2-T MR scanner to investigate the IVIM perfusion model and 2 non-Gaussian diffusion models for evaluating tissue characteristics in rodent gliomas [41,42]. IVIM maps highlighted tumor areas as generally heterogeneous, as confirmed by histology, and f was significantly higher in tumors than in contralateral tissue (P<0.001), as would be expected given neovascularization. Indeed, there was a significant positive correlation between f and microvessel density (R = 0.56, P<0.05), and a negative correlation was found between cellularity and D (R = -0.70, P<0.01).
Infiltrating tumor cells are present in the perivascular spaces in areas of vasogenic edema around the CER [43]. Higher rCBV ratios have been found in NCER surrounding gliomas than in the NCER surrounding metastatic lesions [44]. As would be expected given the vascularity of the CER and NCER, we found positive correlations between fCER and rCBFNCER and rCBVNCER. Furthermore, we found that fCER strongly correlated with D*NCER. The D*value is considered proportional to mean capillary segment length and average blood velocity [18]; like the f value, the D*value may depend on the attenuation of the tumor microvessels and may correlate with the degree of angiogenesis with intact vessels, probably more frequent in the NCER, in terms of basement membrane thickness and pericyte coverage [45]. Further research could shed more light on potential IVIM-metrics to characterize the NCER of high-grade gliomas.
Several authors have used IVIM to separate the diffusion and perfusion components of DWI data, highlighting its potential value in clinical practice[12,18,46]. In healthy volunteers, Wirestam et al. [46] demonstrated modest but significant correlations between f and CBV (R = 0.56;P<0.001) and between CBF obtained from the median value of D* in IVIM and CBF obtained from DSC-MRI (R = 0.35;P<0.001). Federau et al. [18] demonstrated that f, D*, and fD* parameters change gradually under a hypercapnia and hyperoxygenation challenge in the brain. Bisdas et al. [14] recently reported that f was significantly correlated with DSC-derived vascular plasma volume and vascular transit time in healthy brain tissue, whereas in tumor regions, DSC-derived plasma flow was positively correlated with D* and inversely correlated with f.
The IVIM method has many theoretical advantages over DSC-MRI. Unlike DSC-MRI, IVIM perfusion-related parameters can be obtained using DWI datasets without the need for intravenous contrast agents, an important advantage considering that some agents are contraindicated in some patients due to the risk of nephrogenic systemic fibrosis [47]. Moreover, whereas DSC-MRI requires knowledge of the arterial input function, IVIM is intrinsically quantitative and the intravoxel excitation and readout obviates the need for this measurement [48,49]. Likewise, DSC-MRI requires a contralateral normalization measurement that can be difficult to obtain because of anatomical distortion, whereas IVIM-metrics are obtained through direct assessment of the tumor. IVIM diffusion and perfusion parameters might also be useful for guiding biopsy within gliomas [50]. Federau et al. [13] recently reported that T2-prepared IVIM inversion recovery acquisition seems to increase the quantitative blood volume contrast and contrast-to-noise ratio compared to standard IVIM acquisition and DSC-MRI, improving subjective lesion detection, contrast quality, and diagnostic confidence.
Some important limitations of this study merit comment. This pilot study was done at a single center, and the patient sample was too small to draw any definite conclusions about the usefulness of IVIM-metrics for patient management and survival prediction. The reported cutoffs most probably are not optimal and need to be validated. Although the large vessels try to be excluded, even small vessels could potentially affect the signal intensity or calculated IVIM map. IVIM is sensitive to motion if the curve is fitted on a voxel-by-voxel basis, so unavoidable patient movements may be problematic. Susceptibility inhomogeneities, as might occur around the petrous apex or the paranasal sinuses or due to the presence of blood postoperatively, for example, can harm the IVIM signal; however, they are also problematic in DSC-MRI. Analyzing the non-Gaussian diffusion behavior of water (kurtosis or biexponential model) can potentially provide information on microcirculation and tissue microstructure [28], but our diffusion images were acquired with maximum b-value of 1000 s/mm2, limiting our ability to go deeper in the analysis of non-Gaussian water movement due to cellular membrane boundaries [51,52]. The manual placement of the ROIs was subjective. Differences in slice thickness and spatial resolution in DSC, IVIM and T1WI would raise issues of systematic bias to match across different modalities, In our study, f and D* maps was particularly noisy and did not show clear anatomical findings (Fig 1), which might be expected since these parameters should be tissue specific (e.g. GM has much higher perfusion than WM). Better signal to noise ratio from higher field [14] or more averages must be sought for future studies. Partial-volume contamination from cerebrospinal fluid or necrotic areas may have varied during the study. IVIM and DSC have similar spatial resolution; however, DSC data probably has higher SNR for deriving perfusion parameters, whereas the SNR in IVIM can make it difficult to extract the perfusion component reliably. Moreover, IVIM suffers just as much as DSC from the effects of large-vessel partial volumes, but the IVIM equation is still an approximation of the signal model to the data and is not really quantitative in the sense that the ADC measurement can be (when data from low b-values are excluded). Another limitation of IVIM is that the level of perfusion must be high enough before IVIM can reliably detect and measure a perfusion-related parameter. Grech-Sollars et al [53] found that the IVIM parameter f had a poorer inter-scanner coefficient of variation when scanners of different field strengths were combined for normal, and the parameter was also affected by the scan acquisition resolution, for which perfusion is lower than the GBM angiogenic core.
Conclusions
In summary, IVIM seems feasible for evaluating the diffusion and perfusion characteristics of glioblastoma, and fCER and D*CER correlate well with response to therapy and survival. Knowing which patients will respond better to treatment is important for individualizing care, so these parameters might help improve outcomes. Further studies are warranted to test the generalizability of our findings to other cohorts to determine whether IVIM-metrics can be used as perfusion biomarkers in gliomas.
Supporting Information
[Figure omitted. See PDF.]
S1 Fig. Bland-Altman plots.
The plots show that practically all the values are within the confidence limits.
https://doi.org/10.1371/journal.pone.0158887.s001
(TIF)
Author Contributions
Conceived and designed the experiments: JP JS-G GB SP. Performed the experiments: JP JS-G GB. Analyzed the data: JP JS-G GB PD-i-E. Contributed reagents/materials/analysis tools: JP JS-G GB PD-i-E CB. Wrote the paper: JP JS-G GB PD-i-E AA-B CF KN ME RJ MW SP. Statistical analysis: PD-i-E.
Citation: Puig J, Sánchez-González J, Blasco G, Daunis-i-Estadella P, Federau C, Alberich-Bayarri Á, et al. (2016) Intravoxel Incoherent Motion Metrics as Potential Biomarkers for Survival in Glioblastoma. PLoS ONE 11(7): e0158887. https://doi.org/10.1371/journal.pone.0158887
1. Goldbrunner RH, Wagner S, Roosen K, Tonn JC. Models for assessment of angiogenesis in gliomas. J Neurooncol. 2000;50:53–62 pmid:11245281
2. Linkous AG, Yazlovitskaya EM. Angiogenesis in glioblastoma multiforme: navigating the maze. Anticancer Agents Med Chem. 2011;11:712–718 pmid:21707499
3. Aronen HJ, Gazit IE, Louis DN, Buchbinder BR, Pardo FS, Weisskoff RM, et al. Cerebral blood volume maps of gliomas: comparison with tumor grade and histologic findings. Radiology. 1994;191:41–51 pmid:8134596
4. Law M, Yang S, Wang H, Babb JS, Johnson G, Cha S, et al. Glioma grading: sensitivity, specificity, and predictive values of perfusion MR imaging and proton MR spectroscopic imaging compared with conventional MR imaging. AJNR Am J Neuroradiol. 2003;24:1989–1998 pmid:14625221
5. Law M, Young RJ, Babb JS, Peccerelli N, Chheang S, Gruber ML, et al. Gliomas: predicting time to progression or survival with cerebral blood volume measurements at dynamic susceptibility-weighted contrast-enhanced perfusion MR imaging. Radiology. 2008;247:490–498. pmid:18349315
6. Cha S, Knopp EA, Johnson G, Wetzel SG, Litt AW, Zagzag D. Intracranial mass lesions: dynamic contrast-enhanced susceptibility-weighted echo-planar perfusion MR imaging. Radiology. 2002; 223:11–29 pmid:11930044
7. Paulson ES, Schmainda KM. Comparison of dynamic susceptibility-weighted contrast-enhanced MR methods: recommendations for measuring relative cerebral blood volume in brain tumors. Radiology. 2008; 249:601–613 pmid:18780827
8. Shin JH, Lee HK, Kwun BD, Kim JS, Kang W, Choi CG, et al. Using relative cerebral blood flow and volume to evaluate the histopathologic grade of cerebral gliomas: preliminary results. AJR Am J Roentgenol. 2002;179:783–789 pmid:12185064
9. Federau C, O'Brien K, Meuli R, Hagmann P, Maeder P. Measuring brain perfusion with intravoxel incoherent motion (IVIM): initial clinical experience. J Magn Reson Imaging. 2014; 39:624–632 pmid:24068649
10. Hu YC, Yan LF, Wu L, Du P, Chen BY, Wang L, et al. Intravoxel incoherent motion diffusion-weighted MR imaging of gliomas: efficacy in preoperative grading. Sci Rep. 2014 1;4:7208 pmid:25434593
11. Federau C, Meuli R, O'Brien K, Maeder P, Hagmann P. Perfusion measurement in brain gliomas with intravoxel incoherent motion MRI. AJNR Am J Neuroradiol. 2014; 35:256–262 pmid:23928134
12. Bisdas S, Koh TS, Roder C, Braun C, Schittenhelm J, Ernemann U. Intravoxel incoherent motion diffusion-weighted MR imaging of gliomas: feasibility of the method and initial results. Neuroradiology. 2013; 55:1189–1196 pmid:23852430
13. Federau C, O'Brien K. Increased brain perfusion contrast with T₂-prepared intravoxel incoherent motion (T2prep IVIM) MRI. NMR Biomed. 2015; 28:9–16 pmid:25303668
14. Bisdas S, Braun C, Skardelly M, Schittenhelm J, Teo TH, Thng CH, et al. Correlative assessment of tumor microcirculation using contrast-enhanced perfusion MRI and intravoxel incoherent motion diffusion-weighted MRI: is there a link between them? NMR Biomed. 2014;27:1184–1191 pmid:25088433
15. Sasaki M, Sumi M, Eida S, Katayama I, Hotokezaka Y, Nakamura T. Simple and reliable determination of intravoxel incoherent motion parameters for the differential diagnosis of head and neck tumors. PLoS One. 2014; 9:e112866 pmid:25402436
16. Sumi M, Nakamura T. Head and neck tumors: assessment of perfusion-related parameters and diffusion coefficients based on the intravoxel incoherent motion model. AJNR Am J Neuroradiol. 2013; 34:410–416 pmid:22859281
17. Sumi M, Van Cauteren M, Sumi T, Obara M, Ichikawa Y, Nakamura T. Salivary gland tumors: use of intravoxel incoherent motion MR imaging for assessment of diffusion and perfusion for the differentiation of benign from malignant tumors. Radiology. 2012; 263:770–777 pmid:22447854
18. Federau C, Maeder P, O'Brien K, Browaeys P, Meuli R, Hagmann P. Quantitative measurement of brain perfusion with intravoxel incoherent motion MR imaging. Radiology. 2012; 265:874–881 pmid:23074258
19. Le Bihan D, Breton E, Lallemand D, Aubin ML, Vignaud J, Laval-Jeantet M. Separation of diffusion and perfusion in intravoxel incoherent motion MR imaging. Radiology. 1988; 168:497–505 pmid:3393671
20. Latour LL, Svoboda K, Mitra PP, Sotak CH. Time-dependent diffusion of water in a biological model system. Proc Natl Acad Sci. 1994; 91:1229–1233 pmid:8108392
21. Chandarana H, Kang SK, Wong S, Rusinek H, Zhang JL, Arizono S, et al. Diffusion-weighted intravoxel incoherent motion imaging of renal tumors with histopathologic correlation. Invest Radiol. 2012; 47: 688–696 pmid:22996315
22. Wang LL, Lin J, Liu K, Chen CZ, Liu H, Lv P, et al. Intravoxel incoherent motion diffusion-weighted MR imaging in differentiation of lung cancer from obstructive lung consolidation: comparison and correlation with pharmacokinetic analysis from dynamic contrast-enhanced MR imaging. Eur Radiol. 2014; 24:1914–1922 pmid:24788038
23. Liu C, Liang C, Liu Z, Zhang S, Huang B. Intravoxel incoherent motion (IVIM) in evaluation of breast lesions: comparison with conventional DWI. Eur J Radiol. 2013; 82:e782–789 pmid:24034833
24. Pang Y, Turkbey B, Bernardo M, Kruecker J, Kadoury S, Merino MJ, et al. Intravoxel incoherent motion MR imaging for prostate cancer: an evaluation of perfusion fraction and diffusion coefficient derived from different b-value combinations. Magn Reson Med. 2013; 69:553–562 pmid:22488794
25. Andreou A, Koh DM, Collins DJ, Blackledge M, Wallace T, Leach MO, et al. Measurement reproducibility of perfusion fraction and pseudodiffusion coefficient derived by intravoxel incoherent motion diffusion-weighted MR imaging in normal liver and metastases. Eur Radiol. 2013; 23:428–434 pmid:23052642
26. Koh DM, Collins DJ, Orton MR. Intravoxel incoherent motion in body diffusion-weighted MRI: reality and challenges. AJR Am J Roentgenol. 2011; 196:1351–1561 pmid:21606299
27. Kim HS, Suh CH, Kim N, Choi CG, Kim SJ. Histogram analysis of intravoxel incoherent motion for differentiating recurrent tumor from treatment effect in patients with glioblastoma: initial clinical experience. AJNR Am J Neuroradiol. 2014; 35:490–497 pmid:23969343
28. Iima M, Reynaud O, Tsurugizawa T, Ciobanu L, Li JR, Geffroy F, et al. Characterization of glioma microcirculation and tissue features using intravoxel incoherent motion magnetic resonance imaging in a rat brain model. Invest Radiol. 2014; 49:485–490 pmid:24619211
29. Olson JJ, Fadul CE, Brat DJ, Mukundan S, Ryken TC. Management of newly diagnosed glioblastoma: guidelines development, value and application. J Neurooncol. 2009; 93:1–23 pmid:19430879
30. Pope WB, Sayre J, Perlina A, Villablanca JP, Mischel PS, Cloughesy TF. MR imaging correlates of survival in patients with high-grade gliomas. AJNR Am J Neuroradiol. 2005; 26:2466–2474 pmid:16286386
31. Carroll TJ, Rowley HA, Haughton VM. Automatic calculation of the arterial input function for cerebral perfusion imaging with MR imaging. Radiology. 2003; 227:593–600 pmid:12663823
32. Ostergaard L, Johannsen P, Host-Poulsen P, Vestergaard-Poulsen P, Asboe H, Gee AD, et al. Cerebral blood flow measurements by magnetic resonance imaging bolus tracking: comparison with [(15)O]H2O positron emission tomography in humans. J Cereb Blood Flow Metab. 1998; 18:935–940 pmid:9740096
33. Boxerman JL, Schmainda KM, Weisskoff RM. Relative cerebral blood volume maps corrected for contrast agent extravasation significantly correlate with glioma tumor grade, whereas uncorrected maps do not. AJNR Am J Neuroradiol. 2006; 27:859–867 pmid:16611779
34. More J, Garbow B, Hillstrom K. User Guide for MINPACK-1’ Technical Report ANL-80-74, Argonne National Laboratory. 1980
35. Myles PS, Cui J. Using the Bland-Altman method to measure agreement with repeated measures. Br J Anaesth. 2007; 99:309–311 pmid:17702826
36. Rock K, McArdle O, Forde P, Dunne M, Fitzpatrick D, O'Neill B, et al. A clinical review of treatment outcomes in glioblastoma multiforme—the validation in a non-trial population of the results of a randomised phase III clinical trial: has a more radical approach improved survival? Br J Radiol. 2012; 85:e729–e733 pmid:22215883
37. Stupp R, Mason WP, van den Bent MJ, Weller M, Fisher B, Taphoorn MJ, et al. Radiotherapy plus concomitant and adjuvant temozolomide for glioblastoma. N Engl J Med. 2005; 352:987–996 pmid:15758009
38. Ronning PA, Helseth E, Meling TR, Johannesen TB. A population-based study on the effect of temozolomide in the treatment of glioblastoma multiforme. Neuro Oncol. 2012;14:1178–1184 pmid:22869622
39. Kono K, Inoue Y, Nakayama K, Shakudo M, Morino M, Ohata K, et al. The role of diffusion-weighted imaging in patients with brain tumors. AJNR Am J Neuroradiol. 2001; 22:1081–1088 pmid:11415902
40. Sugahara T, Korogi Y, Kochi M, Ikushima I, Shigematu Y, Hirai T, et al. Usefulness of diffusion-weighted MRI with echo-planar technique in the evaluation of cellularity in gliomas. J Magn Reson Imaging. 1999; 9:53–60 pmid:10030650
41. Jensen JH, Helpern JA. MRI quantification of non-Gaussian water diffusion by kurtosis analysis. NMR Biomed. 2010; 23:698–710 pmid:20632416
42. Niendorf T, Dijkhuizen RM, Norris DG, van Lookeren Campagne M, Nicolay K. Biexponential diffusion attenuation in various states of brain tissue: implications for diffusion-weighted imaging. Magn Reson Med. 1996; 36:847–857 pmid:8946350
43. Giese A, Bjerkvig R, Berens ME, Westphal M. Cost of migration: invasion of malignant gliomas and implications for treatment. J Clin Oncol. 2003; 21:1624–1636 pmid:12697889
44. Hakyemez B, Erdogan C, Gokalp G, Dusak A, Parlak M. Solitary metastases and high-grade gliomas: radiological differentiation by morphometric analysis and perfusion-weighted MRI. Clin Radiol. 2010; 65:15–20 pmid:20103416
45. Jain R, Gutierrez J, Narang J, Scarpace L, Schultz LR, Lemke N, et al. In vivo correlation of tumor blood volume and permeability with histological and molecular angiogenic markers in gliomas. AJNR Am J Neuroradiol. 2011; 32:388–394 pmid:21071537
46. Wirestam R, Borg M, Brockstedt S, Lindgren A, Holtås S, Ståhlberg F. Perfusion-related parameters in intravoxel incoherent motion MR imaging compared with CBV and CBF measured by dynamic susceptibility-contrast MR technique. Acta Radiol. 2001; 42:123–128 pmid:11281143
47. Prince MR, Zhang HL, Roditi GH, Leiner T, Kucharczyk W. Risk factors for NSF: a literature review. J Magn Reson Imaging. 2009; 30:1298–1308 pmid:19937930
48. Calamante F, Gadian DG, Connelly A. Quantification of perfusion using bolus tracking magnetic resonance imaging in stroke: assumptions, limitations, and potential implications for clinical use. Stroke. 2002; 33:1146–1151 pmid:11935075
49. Knutsson L, Stahlberg F, Wirestam R. Absolute quantification of perfusion using dynamic susceptibility contrast MRI: pitfalls and possibilities. MAGMA. 2010; 23:1–21 pmid:19960361
50. Weber MA, Henze M, Tüttenberg J, Stieltjes B, Meissner M, Zimmer F, et al. Biopsy targeting gliomas: do functional imaging techniques identify similar target areas? Invest Radiol. 2010; 45:755–768 pmid:20829706
51. Le Bihan D.Apparent diffusion coefficient and beyond: what diffusion MR imaging can tell us about tissue structure. Radiology. 2013; 268:318–322 pmid:23882093
52. Yuan J, Yeung DK, Mok GS, Bhatia KS, Wang YX, Ahuja AT, et al. Non-Gaussian analysis of diffusion weighted imaging in head and neck at 3T: a pilot study in patients with nasopharyngeal carcinoma. PLoS One. 2014 23;9:e87024 pmid:24466318
53. Grech-Sollars M, Hales PW, Miyazaki K, Raschke F, Rodriguez D, Wilson M, et al. Multi-centre reproducibility of diffusion MRI parameters for clinical sequences in the brain. NMR Biomed. 2015;28:468–485. pmid:25802212
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2016 Puig et al. This is an open access article distributed under the terms of the Creative Commons Attribution License: http://creativecommons.org/licenses/by/4.0/ (the “License”), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
Objective
Intravoxel incoherent motion (IVIM) is an MRI technique with potential applications in measuring brain tumor perfusion, but its clinical impact remains to be determined. We assessed the usefulness of IVIM-metrics in predicting survival in newly diagnosed glioblastoma.
Methods
Fifteen patients with glioblastoma underwent MRI including spin-echo echo-planar DWI using 13 b-values ranging from 0 to 1000 s/mm2. Parametric maps for diffusion coefficient (D), pseudodiffusion coefficient (D*), and perfusion fraction (f) were generated for contrast-enhancing regions (CER) and non-enhancing regions (NCER). Regions of interest were manually drawn in regions of maximum f and on the corresponding dynamic susceptibility contrast images. Prognostic factors were evaluated by Kaplan-Meier survival and Cox proportional hazards analyses.
Results
We found that fCER and D*CER correlated with rCBFCER. The best cutoffs for 6-month survival were fCER>9.86% and D*CER>21.712 x10−3mm2/s (100% sensitivity, 71.4% specificity, 100% and 80% positive predictive values, and 80% and 100% negative predictive values; AUC:0.893 and 0.857, respectively). Treatment yielded the highest hazard ratio (5.484; 95% CI: 1.162–25.88; AUC: 0.723; P = 0.031); fCER combined with treatment predicted survival with 100% accuracy.
Conclusions
The IVIM-metrics fCER and D*CER are promising biomarkers of 6-month survival in newly diagnosed glioblastoma.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer