# Microstructure Characterization of Bone Metastases from Prostate Cancer with Diffusion MRI: Preliminary Findings

^{1}Centre for Medical Image Computing, University College London, London, United Kingdom^{2}Department of Radiation Oncology, Odette Cancer Centre, Toronto, Canada^{3}CR-UK and EPSRC Cancer Imaging Centre, Institute of Cancer Research, Royal Marsden NHS Foundation Trust, London, United Kingdom^{4}Radiology, Royal Marsden NHS Foundation Trust, Institute of Cancer Research, Sutton, United Kingdom^{5}Siemens Healthcare GmbH, Erlangen, Germany^{6}Clinical Imaging Research Centre, National University of Singapore, Singapore, Singapore

**Purpose:** To examine the usefulness of rich diffusion protocols with high *b*-values and varying diffusion time for probing microstructure in bone metastases. Analysis techniques including biophysical and mathematical models were compared with the clinical apparent diffusion coefficient (ADC).

**Methods:** Four patients were scanned using 13 *b*-values up to 3,000 s/mm^{2} and diffusion times ranging 18–52 ms. Data were fitted to mono-exponential ADC, intravoxel incoherent motion (IVIM), Kurtosis and Vascular, extracellular, and restricted diffusion for cytometry in tumors (VERDICT) models. Parameters from the models were compared using correlation plots.

**Results:** ADC and IVIM did not fit the data well, failing to capture the signal at high *b*-values. The Kurtosis model best explained the data in many voxels, but in voxels exhibiting a more time-dependent signal, the VERDICT model explained the data best. The ADC correlated significantly (*p* < 0.004) with the intracellular diffusion coefficient (*r* = 0.48), intracellular volume fraction (*r* = −0.21), and perfusion fraction (*r* = 0.46) parameters from VERDICT, suggesting that these factors all contribute to ADC contrast. The mean kurtosis correlated with the intracellular volume fraction parameter (*r* = 0.26) from VERDICT, consistent with the hypothesis that kurtosis relates to cellularity, but also correlated weakly with the intracellular diffusion coefficient (*r* = 0.18) and cell radius (*r* = 0.16) parameters, suggesting that it may be difficult to attribute physical meaning to kurtosis.

**Conclusion:** Both Kurtosis and VERDICT explained the diffusion signal better than ADC and IVIM, primarily due to poor fitting at high *b*-values in the latter two models. The Kurtosis and VERDICT models captured information at high *b* using parameters (Kurtosis or intracellular volume fraction and radius) that do not have a simple relationship with ADC and that may provide additional microstructural information in bone metastases.

## Introduction

Bone metastases, the most common site of metastatic disease in prostate cancer, are associated with reduced quality of life and an increased risk of complications from bone weakness. Histology is not generally feasible in bone and there are currently no imaging criteria for detecting benefit or improvement following therapy in osteoblastic bone metastases (1). There is, therefore, a need for new imaging techniques and those sensitive to the tissue microstructure, such as the apparent diffusion coefficient (ADC), have shown promise in predicting response (2).

However, ADC is a relatively broad reflection of the tissue microstructure. Studies in other tumor types demonstrate that more advanced diffusion techniques can provide additional information. For example, the use of more extensive imaging protocols, including high *b*-values, results in higher diagnostic accuracy in prostate cancer (3). Also, the inclusion of varying diffusion times reveals information about important length scales, such as the cell size and organization, as has been demonstrated in breast cancer (4).

With more extensive imaging protocols, it is possible to fit more complex models. For example, the intravoxel incoherent motion (IVIM) model uses low *b*-value data to provide information about vascular perfusion in addition to tissue diffusion (5, 6). The Kurtosis method captures some of the non-Gaussian diffusion effects at higher *b*-values. These may be caused by restricted diffusion such as that inside cells, but the Kurtosis model does not make biophysical assumptions about the source of non-Guassian diffusion (7). In contrast, Vascular, extracellular, and restricted diffusion for cytometry in tumors (VERDICT) is an adaptable framework (8) that uses a biophysical model to describe the tissue in terms of vascular, extracellular, and intracellular features. The vascular compartment is described by a pseudo-diffusion coefficient and a separate extracellular diffusion compartment, similar to the two components in IVIM, although VERDICT allows alternative compartment shapes to be considered. Diffusion in the intracellular space is modeled by a restricted compartment that captures non-Gaussian diffusion effects at high *b*-values.

Although extended diffusion protocols and more complex fitting models have proven useful in revealing non-invasive microstructure indices in breast (9) and prostate (10), they have not yet been examined in bone metastases. This study explores the potential of microstructure characterization in bone metastases using a rich diffusion protocol with *b*-values up to 3,000 s/mm^{2} and a range of diffusion times, then fitting with a variety of diffusion models.

## Materials and Methods

### Patient Cohort

Eight regions of interest (ROIs) were identified by a radiologist (Nina Tunariu, 8 years experience) as bone metastases in four patients with prostate cancer. This study was carried out in accordance with the recommendations of Health Insurance Portability and Accountability Act guidelines and the local Research Ethics Committee with written informed consent from all subjects. All subjects gave written informed consent in accordance with the Declaration of Helsinki. The protocol was approved by the local Research Ethics Committee (Royal Marsden NHS Foundation Trust). One ROI contained a significant proportion of normal fatty bone marrow as result of heterogeneous response to therapy.

### Data Acquisition

Data were collected at 1.5 T (MAGNETOM Avanto, Siemens Healthcare, Erlangen, Germany) using a 24-channel spine matrix coil and 12-channel body matrix receive coil. Thirteen fat-saturated axial diffusion-weighted images (DWI) were acquired using a prototype sequence (single-shot 2D pulsed-gradient spin-echo echoplanar acquisition, 2.97 mm × 2.97 mm resolution, 38 cm × 38 cm field of view, 5 mm slice thickness, repetition time TR = 3 s, 3 signal averages). Each scan had a different *b*-value with a corresponding *b* = 0 image; the diffusion scan parameters and echo times (TEs) are listed in Table 1. Gradients were applied in a single direction aligned with the main B_{0} field (0, 0, 1).

### Data Fitting

Data were fitted voxelwise with six different models, outlined in Table 2, using a maximum likelihood procedure that accounts for Rician noise. Noise was estimated by applying the difference method (11) to repeated *b* = 0 measurements in the target ROI.

Signal equations for the models were calculated as the sum of the individual compartmental diffusion signals, *S _{j}*, weighted by their volume fraction,

*f*

_{j}:

where TE is the echo time. All models included a signal normalization constant (S_{0}) and a T2 time constant (T2). These parameters were estimated from the global fit to all data to account for differing echo times. A separate T2 fit of all *b* = 0 points confirmed that T2 decay was mono-exponential. The signal equations for the compartments, *S*_{j}, are as follows; derivations are found in Ref. (7, 12).

Kurtosis (7) is a phenomenological model that attempts to describe deviations from Gaussian diffusion using the apparent diffusion kurtosis, *K*:

with $b={(\mathrm{\gamma}g\mathrm{\delta})}^{2}(\mathrm{\Delta}-\mathrm{\delta}/3)$, where γ is the gyromagnetic ratio, g is the gradient strength, δ is the gradient duration, and Δ is the gradient separation.

The Ball compartment describes Gaussian diffusion:

The Astrosticks compartment assumes a uniform distribution of Sticks whose diffusion is restricted to a single dimension:

where $\mathrm{\varphi}(z)=2{\mathrm{\pi}}^{-1/2}{\int}_{0}^{z}\text{exp}(-{t}^{2}\text{})\text{\hspace{0.17em}}$ is the error function.

The Sphere compartment calculation uses the GPD approximation (13):

Here, $Y(x)={e}^{-{{\mathrm{\beta}}_{\text{m}}}^{2}{D}_{1}x}$, β_{m} is the mth root of ${J}_{3/2}({\mathrm{\beta}}_{\text{m}}R)-{\mathrm{\beta}}_{\text{m}}R\text{\hspace{0.17em}}{J}_{5/2}({\mathrm{\beta}}_{\text{m}}R)=0$ and J_{ν} is the Bessel function of the first kind, order ν. The summation was carried out over the first 31 roots of the equation.

In this nomenclature, the ADC is a parameter in the single-compartment Ball fit. The IVIM model is the sum of two Ball compartments, one representing vasculature with pseudodiffusion coefficient *D*_{p} and fraction *f*_{p} and a second representing tissue with diffusion coefficient *D*_{t}. VERDICT models are constructed as the sum of signals from vascular pseudodiffusion (p), extracellular (E), and restricted intracellular sphere (I) compartments with fractions *f*_{p}, *f*_{E}, and *f*_{I} and diffusion coefficients *D*_{p}, *D*_{E}, and *D*_{I}, respectively, with *D*_{E} = *D*_{I} assumed. The average radius of the spherical intracellular compartment was modeled by the parameter *R*. Different compartment shapes were tested for the vascular and extracellular spaces in VERDICT, and these models are named in vascular–extracellular–intracellular order using the terminology in Ref. (12): Astrosticks–Ball–Sphere (ABS), Ball–Ball–Sphere (BBS), and Ball–Astrosticks–Sphere (BAS).

The pseudodiffusion coefficient was found to have little influence on the fit and was fixed to *D*_{p} = 1 × 10^{−2} mm^{2}/s and the intra- and extracellular diffusion coefficients were assumed to be equal and less than 3 × 10^{−3} mm^{2}/s. Volume fractions summed to one. Sphere radius *R* was constrained between 0.1 and 20 μm.

### Statistical Analysis

Models were compared using the Akaike Information Criterion (AIC) (14), where lower AICs indicate models that explain the data better after accounting for varying numbers of model parameters. Graphs present the relative difference in AIC: ΔAIC=AIC–AIC_{min}, where AIC_{min} is the lowest AIC of all six models in each voxel.

The relationship between model parameters was examined using correlation plots. Linear regression was performed on all voxels from the seven non-fatty ROIs pooled together and the Pearson correlation coefficient was calculated. Correlations were considered significant if *p* < 0.004 (*p* < 0.05/12, to account for the 12 pairs of variables compared).

## Results

Figure 1 shows a representative DWI for *b* = 3,000 s/mm^{2} with its corresponding *b* = 0 image. Four ROIs defining bone metastases are outlined in cyan. Heterogeneity is visible in the *b* = 0 and *b* = 3,000 s/mm^{2} images, but with different patterns. The lowest average SNR in an ROI at *b* = 3,000 s/mm^{2} was 10.

**Figure 1**. **(A)** *b* = 0 and **(B)** *b* = 3,000 s/mm^{2} diffusion-weighted images for one patient with four regions of interest outlined in cyan.

Figure 2 summarizes the model selection results. Figure 2A presents a map of the model that best explained the data in each voxel in four ROIs of one patient. There were very few voxels where ADC and IVIM models best explained the data. Many voxels were sufficiently explained by the Kurtosis model, but a small number of voxels were better explained by the more complex VERDICT models. This was consistent across all ROIs, as shown by the boxplots in Figure 2B, which demonstrate that ADC and IVIM models have much higher median ΔAIC across all voxels. Kurtosis and two of the VERDICT models (ABS and BBS) both tend to have low median ΔAIC, although the interquartile range of the Kurtosis ΔAIC is higher because some voxels were very poorly explained by Kurtosis. The BAS VERDICT model performed consistently worse than the other VERDICT models. Since the ABS and BBS VERDICT models yielded similar fits and parameters, results for BBS are presented in subsequent figures, but ABS results were similar.

**Figure 2**. Summary of model selection results. **(A)** Map showing the model that best explained the data [lowest Akaike Information Criterion (AIC)] in each voxel for one patient. **(B)** Boxplots of ΔAIC, where the red median lines indicate that Kurtosis, Astrosticks–Ball–Sphere (ABS), and Ball–Ball–Sphere (BBS) models tend to have similarly low AICs in most voxels, but the larger upper range of the Kurtosis box indicates that some voxels are better explained using VERDICT.

Figure 3 shows the fit in more detail for two voxels (columns) that are representative of (a, c, e) voxels where Kurtosis best explained the data and (b, d, f) voxels where VERDICT best explained the data. The simple ADC model (a, b) failed at high *b*-values in both voxels. The Kurtosis model (c, d) is better able to capture non-Gaussian behavior at high *b*. Where Kurtosis best explained the data, the VERDICT fit was similar (c versus e), but Kurtosis was selected as the simpler model. The gray insets separate out three data points with *b* = 800 s/mm^{2}, but with different gradient separations and diffusion times. Where VERDICT best explained the data (d versus f), there was greater diffusion time-dependence relative to the noise (indicated by the red bars). The diffusion time-dependence at constant *b*-value cannot be explained using the Kurtosis model (flat green line in c, d insets), but is captured by the restricted sphere component in VERDICT (f inset).

**Figure 3**. Normalized signal and selected fits are shown for **(A,C,E)** a voxel where Kurtosis best described the data and **(B,D,F)** a voxel where Vascular, extracellular, and restricted diffusion for cytometry in tumors (VERDICT) [Ball–Ball–Sphere (BBS)] best described the data. The apparent diffusion coefficient (ADC) model **(A,B)** failed at high *b*-values in both voxels. In the first voxel **(C,E)**, the Kurtosis and VERDICT models provided similar fits, so, Kurtosis was selected as the simpler model. In the second voxel **(D,F)**, the gray insets showed the signal for three points with *b* = 800 s/mm^{2}, where Kurtosis is unable to capture signal changes with different gradient separations, Δ, and VERDICT better explained this time-dependence. The dashed red lines and the red bars in the insets both indicate the noise level, demonstrating that it is difficult to detect a time-dependent trend in the first voxel beyond the larger relative noise level.

Figure 4 shows parameter maps in one patient for the ADC, Kurtosis, and VERDICT models (maps for other patients are presented in the Supplementary Material). The diffusion coefficients for all models (ADC, *D*_{kurt}, and *D*_{I}) appeared heterogeneous. However, the potential heterogeneity in the intracellular and pseudodiffusion maps, *f*_{I} and *f*_{p}, respectively, had a different pattern. The radius parameter varied between ~6 and 12 μm, although several voxels (not plotted) reached the lower bound allowed by the fitting; these had low *f*_{I} so that *R* was poorly determined.

**Figure 4**. Parameter maps for the apparent diffusion coefficient (ADC), kurtosis, and vascular, extracellular, and restricted diffusion for cytometry in tumors (VERDICT) [Ball–Ball–Sphere (BBS)] models. There was some apparent heterogeneity in the ADC map and the *D*_{kurt} and *D*_{I} parameters from the kurtosis and VERDICT models showed similar patterns. For VERDICT, the intracellular volume fraction (*f*_{I}) and perfusion fraction (*f*_{p}) maps appeared heterogeneous, but with a different pattern. Several voxels in the radius (*R*) maps hit the fitting lower bound and were not plotted. These correspond to regions where the Kurtosis model provided better fit in Figure 2A and had low f_{I}, meaning that the radius was poorly determined. Parameter maps for the remaining patients can be seen in the Supplementary Material.

Figure 5 plots the correlation between the VERDICT (BBS) model parameters and those from the ADC and Kurtosis models. The data from the fatty ROI (ROI8) do not align exactly with the trends from other ROIs. The fatty ROI exhibited the highest ADC, *D*_{kurt} and *D*_{I}, but had low *f*_{I} and *f*_{p}. Low *f*_{I} also resulted in a poorly determined *R* value that defaulted to the lower fitting limit for almost all voxels in this ROI. Due to the deviation from the trends followed by the less fatty voxels and the influence of points at the extreme ends of the parameter space, voxels from the fatty ROI were not included in the linear regression and Pearson correlation calculation.

**Figure 5**. Correlation plots combining voxel data for all patients, with regions of interest (ROIs) separated by color. Linear regression on combined data from the seven non-fatty ROIs demonstrated significant correlation between apparent diffusion coefficient (ADC) and vascular, extracellular, and restricted diffusion for cytometry in tumors (VERDICT) [Ball–Ball–Sphere (BBS)] parameters *D*_{I}, *f*_{I}, and *f*_{p}. These same three VERDICT parameters correlate with *D*_{kurt}, whereas *K* correlates with *f*_{I} and, more weakly, with *R* and *D*_{I}.

The Pearson correlation coefficients are shown above each plot and the asterisk indicates a statistically significant correlation (*p* < 0.004). Both ADC and *D*_{kurt} correlated significantly with *D*_{I} (*p* = 7 × 10^{−25}, 4 × 10^{−35} for ADC and *D*_{kurt}, respectively), *f*_{I} (*p* = 10^{−5}, 2 × 10^{−4}) and *f*_{p} (*p* = 3 × 10^{−23}, 10^{−26}). The kurtosis parameter, K, correlated most strongly with the intracellular water fraction, *f*_{I} (*p* = 10^{−7}), and more weakly with the diffusion coefficient, *D*_{I} (*p* = 3 × 10^{−4}), and radius parameter, *R* (*p* = 10^{−3}). *R* did not significantly correlate with the ADC (*p* = 0.09) or *D*_{kurt} (*p* = 0.01).

## Discussion

Although high *b*-values are beginning to be adopted in primary prostate tumors (10, 15) and brain glioma (16), along with models to characterize cellularity and microstructure, these methods had not been studied in bone metastases, where microstructure differs. Microstructure methods in bone would be particularly valuable due to the challenge of acquiring histological information from these areas. This study acquired diffusion data in bone metastases from protocols including high *b*-values and varying diffusion times, demonstrating that these data were not adequately described by conventional ADC and IVIM models. This is consistent with findings from other tumor types (8, 15, 17) and expected in regions with complex and heterogeneous microstructure, including cellular regions.

The kurtosis and VERDICT models captured the signal at high *b*-values and non-Gaussian effects, likely due to intracellular water in these voxels. Kurtosis was often sufficient to characterize the data. However, there were some voxels where the more complex VERDICT model was required to fully describe the signal, particularly variations with diffusion time. These voxels tended to have the largest intracellular volume fraction estimates in parameter maps, which suggest that VERDICT explains the data better than Kurtosis in regions with a larger intracellular signal contribution, resulting in both higher SNR and stronger time-dependence. Kurtosis is not typically performed with varying diffusion time and may be sufficient to describe the data from such protocols.

Of the three VERDICT models tested, BAS did not fit as well as the ABS or BBS models, suggesting that Astrosticks is not a good description of the extracellular space. BBS and ABS demonstrated similar fits and parameter values. These models differ only in the shape of the vascular compartment, and it is likely that low vascular fraction (*f*_{p} < 0.20) makes it difficult to describe the shape accurately.

The median ADC was calculated in each ROI and the mean across ROIs was 0.85 × 10^{−3} mm^{2}/s (range 0.43–1.78 × 10^{−3} mm^{2}/s; the mean was 0.71 × 10^{−3} mm^{2}/s when the fatty ROI was excluded as an outlier at the upper end of the range). This average is within the range reported in previous bone metastasis studies that used a more limited diffusion protocol (2), although the range in this study skews lower, as expected with the inclusion of higher *b*-values. In cellular regions, images at high *b*-values will be less attenuated, resulting in lower ADCs when a simple exponential model is used to fit all data. Previous studies have noted that the relationship between ADC and cell density is complex in bone marrow due to displacement of fat cells as tumor cells infiltrate and the return of fat cells during treatment response (18). Fat cells may have low water content and thus low signal, complicating the ADC fit. The presence of a wider range of *b*-values and careful accounting for noise in this study may explain the high ADC in the fatty ROI. The VERDICT model parameters also support this: the ADC was negatively correlated with the intracellular fraction *f*_{I} and the estimate of *f*_{I} in the fatty ROI was very low, suggesting little intracellular water influenced the diffusion signal in fat-containing voxels.

The VERDICT model parameters fell within plausible biological ranges. The median *D*_{I} of (1.1 ± 0.3) × 10^{−3} mm^{2}/s was higher than the ADC because restriction is accounted for and in agreement with previous studies (8). There are challenges to obtaining histology in bone metastases, but the prostate cell size range observed in cell culture of 13 ± 3 μm (19) is in agreement with the BBS VERDICT size estimate (2 × *R* = 15 ± 3 μm). The intracellular fraction parameter *f*_{I} was 0.3 ± 0.1, somewhat lower than observed in primary prostate tumors (10), which may be due to the heterogeneous way in which the metastasis infiltrates the bone and partial volume effects. The low perfusion fraction, *f*_{p} = 0.04 ± 0.03, was low, as expected in bone metastases. The maps also suggested patterns of heterogeneity: ADC, *D*_{kurt}, and *D*_{I} appeared similar, with lower values toward the center of the bone, but *f*_{p} and *f*_{I} demonstrate different patterns. The meaning and biological importance of these patterns is of interest for future study.

The correlations were in agreement with previous findings from prostate tumors (10). Both ADC and *D*_{kurt} correlated with *D*_{I}, *f*_{I}, and *f*_{p}, indicating that the parameters from the biophysical VERDICT models all contribute to the values of diffusion coefficients from simpler models and thus these diffusion coefficients are less likely to be specific to microstructural properties. The apparent kurtosis *K* correlated most strongly with *f*_{I}, but showed weaker correlation with *R* and *D*_{I}. This supports the hypothesis that *K* increases with increasing tissue complexity (20), but also indicates that it may be difficult to separate out different microscopic effects using kurtosis. The radius, *R*, had no correlation with ADC, which is expected given that the effects of restriction are strongest at high *b*-values and the ADC model fits these *b*-values poorly.

There are several limitations to this study. It involves a relatively small number of patients. However, even this small number demonstrates that high *b*-values and varying diffusion times provide information not captured by current clinical methods that may be valuable in measuring tumor response. Future work will explore whether VERDICT or Kurtosis parameters improve prediction of patient outcome in a larger cohort.

The models are all necessarily simplifications of complex biological systems to make data fitting tractable. The T2 was modeled as a mono-exponential decay, an assumption confirmed by fitting all *b* = 0 data. This confirms the efficacy of the fat saturation, which eliminates most signal from aliphatic protons, although the use of longer echo times may also attenuate signal from fat in this study. The number of VERDICT parameters was reduced by assuming the intra- and extracellular diffusion coefficients were equal. It has been demonstrated that the use of short diffusion times or specialized oscillating gradient sequences and high SNR are needed to accurately determine the intracellular diffusion coefficient (21). This is not generally feasible with a clinical protocol, necessitating an assumption about the intracellular diffusion and the equality assumption has been successfully applied in other work (10). Alternative assumptions, such as a tortuosity approximation for the extracellular diffusion coefficient, are also possible. Incorrect assumptions about intracellular diffusion would cause systematic bias in the other model parameters, particularly the radius.

Finally, this study was conducted at 1.5 T with gradient separations corresponding to diffusion times ranging from 17.7 to 51.7 ms. Higher-field scanners often have stronger gradients, which would allow shorter diffusion times to be examined. A wider range of scan parameters, along with protocol optimization techniques (22), would allow design of protocols that are more sensitive to restriction effects and may also allow relaxation of the constraints used to make fitting tractable.

In conclusion, this study showed that microstructure methods capture additional information at high *b*-values at 1.5 T in bone metastases. The kurtosis, intracellular volume fraction, and radius parameters from kurtosis and VERDICT fitting are not simply related to the ADC and may yield new information. Obtaining histological samples from bone is very challenging and, therefore, non-invasive methods of inferring tissue characteristics are needed for treatment monitoring. The usefulness of high *b*-value, multi-diffusion time protocols, and more complex models should be further explored for characterizing response to therapy in metastatic bone disease.

## Ethics Statement

This study was carried out in accordance with the recommendations of HIPAA (Health Insurance Portability and Accountability Act) guidelines and the local Research Ethics Committee with written informed consent from all subjects. All subjects gave written informed consent in accordance with the Declaration of Helsinki. The protocol was approved by the local Research Ethics Committee (Royal Marsden NHS Foundation Trust).

## Author Contributions

Study design: CB, EP, and DC. Pulse sequence design: TF. Data acquisition: VM, NT, MO, and DC. Data analysis and interpretation: CB, NT, MO, DC, and EP. Project management: EP, DA, DH, and ML. Manuscript preparation: All.

## Conflict of Interest Statement

TF is an employee of Siemens Healthcare. The other authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Funding

We acknowledge CRUK and EPSRC support to the Cancer Imaging Centre at KCL/UCL and the Cancer Imaging Centre at ICR and RMH in association with MRC and Department of Health C1060/A10334, C1060/A16464 and NHS funding to the NIHR Biomedical Research Centre, Clinical Research Facility in Imaging. This work was also supported by funding from the EPSRC MIMIC project (EP/K020439/1) and the VPH-PRISM EU project (FP7-ICT-2011-9, 601040). MOL is an NIHR Emeritus Senior Investigator. EP is supported by EPSRC fellowship EP/N021967/1. EPSRC grant EP/M020533/1 further supports DA’s and DH’s work on this topic.

## Supplementary Material

The Supplementary Material for this article can be found online at http://www.frontiersin.org/articles/10.3389/fonc.2018.00026/full#supplementary-material.

## References

1. Padhani AR, Lecouvet FE, Tunariu N, Koh D-M, De Keyzer F, Collins DJ, et al. METastasis reporting and data system for prostate cancer: practical guidelines for acquisition, interpretation, and reporting of whole-body magnetic resonance imaging-based evaluations of multiorgan involvement in advanced prostate cancer. *Eur Urol* (2017) 71:81–92. doi:10.1016/j.eururo.2016.05.033

2. Perez-Lopez R, Mateo J, Mossop H, Blackledge MD, Collins DJ, Rata M, et al. Diffusion-weighted imaging as a treatment response biomarker for evaluating bone metastases in prostate cancer: a pilot study. *Radiology* (2017) 283:168–77. doi:10.1148/radiol.2016160646

3. Katahira K, Takahara T, Kwee TC, Oda S, Suzuki Y, Morishita S, et al. Ultra-high-b-value diffusion-weighted MR imaging for the detection of prostate cancer: evaluation in 201 cases with histopathological correlation. *Eur Radiol* (2011) 21:188–96. doi:10.1007/s00330-010-1883-7

4. Teruel JR, Cho GY, Moccaldi RTM, Goa PE, Bathen TF, Feiweier T, et al. Stimulated echo diffusion tensor imaging (STEAM-DTI) with varying diffusion times as a probe of breast tissue. *J Magn Reson Imaging* (2017) 45:84–93. doi:10.1002/jmri.25376

5. 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. doi:10.1148/radiology.168.2.3393671

6. Koh DM, Collins DJ, Orton MR. Intravoxel incoherent motion in body diffusion-weighted MRI: reality and challenges. *AJR Am J Roentgenol* (2011) 196:1351–61. doi:10.2214/AJR.10.5515

7. Jensen JH, Helpern JA, Ramani A, Lu H, Kaczynski K. Diffusional kurtosis imaging: the quantification of non-Gaussian water diffusion by means of magnetic resonance imaging. *Magn Reson Med* (2005) 53:1432–40. doi:10.1002/mrm.20508

8. Panagiotaki E, Walker-Samuel S, Siow B, Johnson SP, Rajkumar V, Pedley RB, et al. Noninvasive quantification of solid tumor microstructure using VERDICT MRI. *Cancer Res* (2014) 74:1902–12. doi:10.1158/0008-5472.CAN-13-2511

9. Bailey C, Siow B, Panagiotaki E, Hipwell JH, Mertzanidou T, Owen J, et al. Microstructural models for diffusion MRI in breast cancer and surrounding stroma: an ex vivo study. *NMR Biomed* (2017) 30:e3679. doi:10.1002/nbm.3679

10. Panagiotaki E, Chan RW, Dikaios N, Ahmed HU, O’Callaghan J, Freeman A, et al. Microstructural characterization of normal and malignant human prostate tissue with vascular, extracellular, and restricted diffusion for cytometry in tumours magnetic resonance imaging. *Invest Radiol* (2015) 50:218–27. doi:10.1097/RLI.0000000000000115

11. Murphy BW, Carson PL, Ellis JH, Zhang YT, Hyde RJ, Chenevert TL. Signal-to-noise measures for magnetic resonance imagers. *Magn Reson Imaging* (1993) 11:425–8. doi:10.1016/0730-725X(93)90076-P

12. Panagiotaki E, Schneider T, Siow B, Hall MG, Lythgoe MF, Alexander DC. Compartment models of the diffusion MR signal in brain white matter: a taxonomy and comparison. *Neuroimage* (2012) 59:2241–54. doi:10.1016/j.neuroimage.2011.09.081

13. Stepišnik J. Time-dependent self-diffusion by NMR spin-echo. *Phys B Condens Matter* (1993) 183:343–50. doi:10.1016/0921-4526(93)90124-O

14. Akaike H. A new look at the statistical model identification. *IEEE Trans Automat Contr* (1974) 19:716–23. doi:10.1109/TAC.1974.1100705

15. Jambor I, Merisaari H, Taimen P, Boström P, Minn H, Pesola M, et al. Evaluation of different mathematical models for diffusion-weighted imaging of normal prostate and prostate cancer using high b-values: a repeatability study. *Magn Reson Med* (2015) 73:1988–98. doi:10.1002/mrm.25323

16. Pramanik PP, Parmar HA, Mammoser AG, Junck LR, Kim MM, Tsien CI, et al. Hypercellularity components of glioblastoma identified by high b-value diffusion-weighted imaging. *Int J Radiat Oncol Biol Phys* (2015) 92:811–9. doi:10.1016/j.ijrobp.2015.02.058

17. Rosenkrantz AB, Sigmund EE, Johnson G, Babb JS, Mussi TC, Melamed J, et al. Prostate cancer: feasibility and preliminary experience of a diffusional kurtosis model for detection and assessment of aggressiveness of peripheral zone cancer. *Radiology* (2012) 264:126–35. doi:10.1148/radiol.12112290

18. Padhani AR, Gogbashian A. Bony metastases: assessing response to therapy with whole-body diffusion MRI. *Cancer Imaging* (2011) 11:S129–45. doi:10.1102/1470-7330.2011.9034

19. Park S, Ang RR, Duffy SP, Bazov J, Chi KN, Black PC, et al. Morphological differences between circulating tumor cells from prostate cancer patients and cultured prostate cancer cells. *PLoS One* (2014) 9:e85264. doi:10.1371/journal.pone.0085264

20. Rosenkrantz AB, Triolo MJ, Melamed J, Rusinek H, Taneja SS, Deng F-M. Whole-lesion apparent diffusion coefficient metrics as a marker of percentage Gleason 4 component within Gleason 7 prostate cancer at radical prostatectomy. *J Magn Reson Imaging* (2015) 41:708–14. doi:10.1002/jmri.24598

21. Li H, Jiang X, Xie J, Gore JC, Xu J. Impact of transcytolemmal water exchange on estimates of tissue microstructural properties derived from diffusion MRI. *Magn Reson Med* (2017) 77:2239–49. doi:10.1002/mrm.26309

Keywords: diffusion MRI, bone metastases, kurtosis, intravoxel incoherent motion, vascular, extracellular, and restricted diffusion for cytometry in tumors

Citation: Bailey C, Collins DJ, Tunariu N, Orton MR, Morgan VA, Feiweier T, Hawkes DJ, Leach MO, Alexander DC and Panagiotaki E (2018) Microstructure Characterization of Bone Metastases from Prostate Cancer with Diffusion MRI: Preliminary Findings. *Front. Oncol.* 8:26. doi: 10.3389/fonc.2018.00026

Received: 30 October 2017; Accepted: 29 January 2018;

Published: 16 February 2018

Edited by:

Andreia Vasconcellos Faria, Johns Hopkins University, United StatesReviewed by:

Piotr Kozlowski, University of British Columbia, CanadaPilar López-Larrubia, Consejo Superior de Investigaciones Científicas (CSIC), Spain

Copyright: © 2018 Bailey, Collins, Tunariu, Orton, Morgan, Feiweier, Hawkes, Leach, Alexander and Panagiotaki. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Colleen Bailey, colleen.bailey@sunnybrook.ca