Distinguishing between different microstructural changes using optimised diffusion-weighted acquisitions
Damien J. McHugh1,2 and Geoff J.M. Parker1,2,3

1Centre for Imaging Sciences, The University of Manchester, Manchester, United Kingdom, 2CRUK & EPSRC Cancer Imaging Centre in Cambridge & Manchester, United Kingdom, 3Bioxydyn Ltd., Manchester, United Kingdom

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 parameters3, potentially hampering the utility of such approaches. Here, optimum experimental design concepts4,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, i, 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 µm2/ms, representing a plausible model of tumour tissue, two possible microstructural changes were considered: (1) a decrease in i 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 f­i to 0.24 (cell density remains constant, mimicking apoptotic cell shrinkage, with a 60% cell volume decrease6). 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, TEmax=100 ms, bmin=150 s/mm2 (avoiding perfusion effects). Sensitivities were used in a Fedorov exchange algorithm7 to find four optimum {G(mT/m), ∆(ms), δ(ms)} combinations for estimating R, fi, Di, and De for each microstructure. D-optimum designs4 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=106, 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=106 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 large6 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/.

Figures

Fig. 1. Top: Starting from a ‘baseline’ microstructure, two distinct changes are considered: (1) volume fraction decrease with cell size constant (decrease in cell density), (2) cell size decrease and coupled volume fraction decrease (constant cell density). Diffusion coefficients are assumed constant. Bottom: D-optimum acquisition parameters for each microstructure.

Fig. 2. Correlations between model parameters for non-optimum (left) and D-optimum (right) acquisitions, for the ‘baseline’ microstructure with SNR = 106. Red crosses indicate ground truth.

Fig. 3. Top: Correlations between model parameters for non-optimum (left) and D-optimum (right) acquisitions, for the ‘baseline’ microstructure with SNR = 50. Red crosses indicate ground truth. Bottom: Model parameter histograms for the data shown above. Red lines indicate ground truth.

Fig. 4. Histograms for R and fi for the ‘baseline’ microstructure and the two changes (top, middle, and bottom rows) with SNR = 50. Left column: Non-optimum acquisition. Middle column: Each microstructure assessed with its own specific D-optimum acquisition. Right column: Each microstructure assessed with D-optimum acquisition for ‘baseline’.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
2030