Synopsis
This work investigates the use of optimised diffusion-weighted acquisitions for distinguishing between different microstructural changes relevant to characterising tumour tissue. Optimised protocols are found for a 'baseline' microstructure, and for two distinct changes which would lead to an ADC increase: (1)
volume fraction decrease
with cell size constant (therefore a decrease in cell density), (2) cell size decrease and
coupled volume fraction decrease (therefore a constant cell density). Model fitting simulations are performed with optimised and non-optimised protocols, demonstrating that the improved precision achieved with optimised protocols may be beneficial in terms of distinguishing between these microstructural changes.Purpose
Applying biophysical models to diffusion MRI data has the
potential to provide more specific information about tissue microstructure,
beyond indices such as the apparent diffusion coefficient (ADC)
1,2.
In particular, where different microstructural changes result in similar ADC
changes, the use of biophysical models may allow these situations (e.g. change in
cell size vs. change in volume fraction) to be distinguished. However, fitting
such models to noisy data can yield imprecise and inaccurate
parameter estimates with artificial correlations between parameters
3,
potentially hampering the utility of such approaches. Here,
optimum experimental design concepts
4,5 are used to investigate the effects
of acquisition parameters on estimates from a simple model of tumour tissue
undergoing distinct microstructural changes. Optimum acquisitions
are found, and used in model-fitting simulations,
comparing results with those obtained using non-optimum acquisitions.
Methods
Model: The normalised
PGSE signal,
S/
S0,
was modelled analytically by combining restricted
diffusion inside a sphere with hindered extracellular diffusion, with four
model parameters: cell radius,
R, intracellular volume fraction,
fi, and
intra- and extra-cellular diffusivities,
Di,
and
De;
S0=1.
Starting from a ‘baseline’
with
R=10 µm,
fi=0.6,
Di=
De=1.5 µm
2/ms,
representing a plausible model of tumour tissue, two
possible microstructural changes were considered:
(1) a decrease in
fi to 0.24 (a cell
density decrease, mimicking complete cell death), (2) a
decrease in
R to 7.37 µm with
an associated decrease in
fi to
0.24 (cell density remains constant, mimicking apoptotic cell shrinkage, with
a 60% cell volume decrease
6). Changes (1) and (2) both give
an ADC increase for a typical multi-
b-value
clinical protocol; see Figure 1, top.
Optimum design: Optimum PGSE
parameters (gradient strength,
G,
separation, ∆, and duration, δ) are those that maximise or minimise some summary statistic of the
signal model’s information matrix,
M4,5.
Sensitivities for
M (∂
S/∂
R,
∂
S/∂
fi, ∂
S/∂
Di, and ∂
S/∂
De) were
calculated numerically for 57016 combinations of {
G, ∆, δ} satisfying
typical clinical constraints:
Gmax=60
mT/m, TE
max=100 ms,
bmin=150 s/mm
2 (avoiding
perfusion effects). Sensitivities were used in a Fedorov exchange algorithm
7 to
find four optimum {
G(mT/m), ∆(ms), δ(ms)} combinations for estimating
R,
fi,
Di, and
De for each microstructure. D-optimum
designs
4 were calculated, corresponding to maximising the
determinant of
M.
Simulations:
Synthetic signals were generated for all three microstructures,
for three acquisition strategies: (a)
D-optimum for each specific microstructure, (b)
D-optimum for the ‘baseline’ microstructure, and
(c) non-optimum, taken as {30,20,15},
{30,72,15}, {60,20,15}, {60,72,15}; i.e. fixed δ, variable
G and ∆. Gaussian noise (SD=
S0/SNR,
5000 instances, SNR=10
6, 50) was added, and the model
was fit using least-squares.
Results and discussion
Figure 1 (bottom) lists D-optimum designs for each
microstructure. Figure 2 plots the correlations between the model parameters for the ‘baseline’ microstructure at SNR=10
6
for non-optimum and D-optimum designs. The D-optimum
acquisition tends to reduce the correlation between parameters (most evident on
the top row) and improves precision. Behaviour at a more realistic, but
still reasonably high, SNR of 50 (Figure
3), shows that
while similar correlation patterns are
observed for the two designs, there
is benefit in using the D-optimum acquisition:
the second peak of the
R distribution
is suppressed,
Di
is estimated with better precision, and
De
and
fi estimates
are more accurate. Figure 4 shows
R
and
fi histograms for the
three acquisition strategies, for each microstructure. The main
benefit of D-optimum designs appears to be in
R estimates: distributions for the two changes are similar and
broad when using non-optimum designs, while D-optimum designs considerably
improve the precision, offering a greater chance of detecting the change in
cell size. This is true for the
best-case scenario where each microstructure is assessed with its own specific
D-optimum design (Figure 4, middle
column), as well as the more realistic scenario where a
single optimised protocol is used for all microstructures (Figure 4, right
column). Note that
the change in
R considered here is relatively large
6 and
sensitivity to smaller changes will be lower, especially as SNR decreases. Precision
of
fi estimates is
generally improved with D-optimum designs, while
Di estimates are generally poor for each strategy
and
De distributions are
similar (
Di and
De results not shown).
Conclusion
D-optimum
designs of diffusion MRI
acquisitions may prove beneficial for distinguishing between
different changes tumour tissue may undergo. Note
that some knowledge of likely underlying microstructure is
required to enable optimisation,
and further work will assess the extent to which a single
optimised protocol can be applied to a
wider range of physiologically relevant microstructures.
Joint consideration of optimised acquisition parameters and the magnitude of
specific biological changes is likely to be important when evaluating the
utility of microstructural models.
Acknowledgements
This is a contribution from the Cancer
Imaging Centre in Cambridge & Manchester, which is funded by the EPSRC and Cancer Research UK (C8742/A18097).
The authors would like to acknowledge the assistance given by IT
Services and the use of the Computational Shared Facility at The
University of Manchester.References
1Assaf et al. Magn Res Med. 2004;52:965-978. 2Panagiotaki et al. Cancer
Res. 2014;74:1-11. 3Jelescu et al. ISMRM. 2015, p.
1024. 4Atkinson et al. Optimum Experimental Designs, with SAS.
Oxford University Press, 2007. 5Alexander. Magn Res Med. 2008;60:439-448
. 6Poulsen
et al. Am J Physiol Cell Physiol. 2010;298: C14–C25. 7Wheeler 2004. optFederov. AlgDesign. The R project for statistical
computing http://www.r-project.org/.