Synopsis
In this study, we investigate the
feasibility of using the recently proposed diffusion MRI method, POMACE
(Pulsed and Oscillating gradient MRI Assessment
of Cell size and Extracellular space) for characterization of tissue
changes induced by chemotherapy. Surface-to-volume ratio (S/V), cell
size (Rcell) and extracellular space volume fraction (ECS) are estimated
prior to and after 5FU and bevacizumab treatment
in GL261 murine glioma model. Preliminary results indicate an early ECS
decrease two days after 5FU/bevacizumab administration. By combining
OGSE and PGSE, POMACE provides non-invasive and quantitative metrics
regarding tumor microstructure that can be applied
to evaluate treatment response.Purpose
Diffusion weighted imaging (DWI) is a
promising tool to measure cellular structural characteristics non-invasively. However,
cellular restriction of water diffusion in cancer cells is typically expected
in 5-10 μm range which cannot be fully assessed by using either pulsed gradient
spin echo (PGSE) or oscillating gradient spin echo (OGSE) diffusion experiment
alone. Recent studies [1,2] showed that combining PGSE and OGSE DWI experiment
allowed for characterization of the tumor microstructure using quantitative parameters,
such as surface-to-volume ratio (S/V),
cell size (R
cell) or
extracellular space fraction (ECS). The
purpose of this study was to investigate the utility of such techniques to characterize
the early treatment response in a mouse glioblastoma model. The in vivo imaging
findings will be compared with histopathological assessment.
Methods
GL261 murine GBM cells (10
6
in 5μL) were injected into the subcortex of 6-8 week old C57BL/6 mice (n=9,
female) under anesthesia (air +3% isoflurane). The MR protocol consisted
of 10 OGSE and 4 PGSE measurements at 7T (Bruker Biospec Avance2, Etlingen),
probing the diffusion coefficient D(t) with diffusion times = 6/8/16/31 ms for
PGSE and frequencies in the range of 65-225 Hz for the fast ramp cos-OGSE (TR/TE=3000/70 ms,
BW=300 kHz, 1 readout segment, NA=20, NR=2, res. 250x250x1500 μm, FoV 20x20 mm,
b=[0,200,400] s/mm
2, dir. (1,1,1), TA=6 min, total time 84 min).
Mice were imaged prior and after tumor treatment: pre-treatment MRI followed by
bevacizumab (10 mg/kg) on the first day; 5FU (80 mg/kg) on the second day;
post-treatment MRI on the third day. Mice were then sacrificed for ex vivo
imaging, histology and light microscopy. S/V
and free diffusivity (D
0)
are estimated using the linear dependence of D(t) with ω
-1/2 in the
very short time regime, as in [1]. The extracellular fraction (ECS), cell size
(R
cell) and intra and extracellular diffusivities (D
intra/D
ecs) are estimated using a two-compartment geometrical
model for the tumor environment. Cells were assumed impermeable for diffusion
times t < 30 ms. Diffusion in the ECS was considered Gaussian and in the
tortuosity limit, as a first approximation [2]. This approach of using the
diffusion time range suitable for cancer cells is also shown in Figure 1 and referred to as POMACE (Pulsed and Oscillating
gradient MRI Assessment of Cell size and Extracellular space). The
fit stability was improved by setting D
ecs=2.2
μm
2/ms. Parameter estimation was performed in Matlab (The Natick,
MA) using non-linear fitting at voxel level to derive parametric maps of ECS, R
cell,
S/V, D
0, D
intra
and D
ecs in the
tumor. The POMACE parameters were assessed at tumor level for statistical
significance between pre- and post-treatment imaging sessions in treated and control groups.
Results
D(t) was highly dependent
on diffusion time / oscillation frequency inside the tumor (+80% D(t) between t=31ms
and ω=225Hz, Fig. 2). S/V, ECS and R
cell maps revealed strong microstructure heterogeneity
within the tumor (Fig. 3). Regions of low S/V
and high R
cell correlated
well with the areas of high contrast enhancement (arrows in Fig. 3). As
expected for spheres, the total S/V
correlated very well with 6*(1-ECS)/R
cell inside the GBM (R
2=0.87
after exclusion of rim voxels). A 17% decrease in S/V was found compared to POMACE-based estimations. At tumor level,
the ECS ratio post/pre-treatment was
found significantly lower than unity for the treated group (ΔECS=-5%, N=4, P<0.02, Fig. 4). No
significant difference was found between groups for S/V, R
cell, D
0, or D
intra. Changes in tumor growth (pre/post-treatment)
were not found significantly different between control and treated groups (+49%
vs. +65%, P>0.5).
Discussion and conclusion
Compared to contrast
enhanced MRI and conventional DWI, time-dependent diffusion allowed to
disentangle microstructural regional variations observed in post contrast images
into objective quantifiable geometrical metrics (S/V, R
cell, ECS) that can be assessed with different
modalities. Using POMACE, a novel way to combine OGSE and PGSE data, S/V changes can be decomposed into cell
fraction (ECS) and cell size variations using the D(t) time-dependence at short
frequencies / long times. Comparison of S/V,
ECS and R
cell measurements suggest (a) partial cellular membrane
deterioration, and/or (b) decrease of the apparent restrictive surface due to
increase cellular packing. In addition, changes in tumor structure , such as ECS,
prior to changes in tumor growth were detected in the treated group, although
additional analysis including ex vivo MRI and histopathological assessment are
required in order to validate these changes and establish the time-dependent
diffusion parameters a set of objective markers of treatment response.
Acknowledgements
This work was supported by NIH R01 CA160620.References
[1] Reynaud et al. Magn.
Reson Med (2015) doi:10.1002/mrm.25865; [2] Reynaud et al. Proc. ISMRM 2015
(program number 344)