Diffusion model comparison identifies distinct tumor sub-regions and tracks treatment response.

PURPOSE: MRI biomarkers of tumor response to treatment are typically obtained from parameters derived from a model applied to pre-treatment and post-treatment data. However, as tumors are spatially and temporally heterogeneous, different models may be necessary in different tumor regions, and model suitability may change over time. This work evaluates how the suitability of two diffusion-weighted (DW) MRI models varies spatially within tumors at the voxel level and in response to radiotherapy, potentially allowing inference of qualitatively different tumor microenvironments. METHODS: DW-MRI data were acquired in CT26 subcutaneous allografts before and after radiotherapy. Restricted and time-independent diffusion models were compared, with regions well-described by the former hypothesized to reflect cellular tissue, and those well-described by the latter expected to reflect necrosis or oedema. Technical and biological validation of the percentage of tissue described by the restricted diffusion microstructural model (termed %MM) was performed through simulations and histological comparison. RESULTS: Spatial and radiotherapy-related variation in model suitability was observed. %MM decreased from a mean of 64% at baseline to 44% 6 days post-radiotherapy in the treated group. %MM correlated negatively with the percentage of necrosis from histology, but overestimated it due to noise. Within MM regions, microstructural parameters were sensitive to radiotherapy-induced changes. CONCLUSIONS: There is spatial and radiotherapy-related variation in different models' suitability for describing diffusion in tumor tissue, suggesting the presence of different and changing tumor sub-regions. The biological and technical validation of the proposed %MM cancer imaging biomarker suggests it correlates with, but overestimates, the percentage of necrosis.

Towards a 'resolution limit' for DW-MRI tumor microstructural models: A simulation study investigating the feasibility of distinguishing between microstructural changes.

PURPOSE: To determine the feasibility of extracting sufficiently precise estimates of cell radius, R, and intracellular volume fraction, fi , from DW-MRI data in order to distinguish between specific microstructural changes tissue may undergo, specifically focusing on cell death in tumors. METHODS: Simulations with optimized and non-optimized clinical acquisitions were performed for a range of microstructures, using a two-compartment model. The ability to distinguish between (i) cell shrinkage with cell density constant, mimicking apoptosis, and (ii) cell size constant with cell density decreasing, mimicking loss of cells, was evaluated based on the precision of simulated parameter estimates. Relationships between parameter precision, SNR, and the magnitude of specific parameter changes, were used to infer SNR requirements for detecting changes. RESULTS: Accuracy and precision depended on microstructural properties, SNR, and the acquisition protocol. The main benefit of optimized acquisitions tended to be improved accuracy and precision of R, particularly for small cells. In most cases considered, higher SNR was required for detecting changes in R than for changes in fi . CONCLUSIONS: Given the relative changes in R and fi due to apoptosis, simulations indicate that, for a range of microstructures, detecting changes in R require higher SNR than detecting changes in fi , and that such SNR is typically not achieved in clinical data. This suggests that if apoptotic cell size decreases are to be detected in clinical settings, improved SNR is required. Comparing measurement precision with the magnitude of expected biological changes should form part of the validation process for potential biomarkers.

Mapping Hypoxia in Renal Carcinoma with Oxygen-enhanced MRI: Comparison with Intrinsic Susceptibility MRI and Pathology.

Purpose To cross-validate T1-weighted oxygen-enhanced (OE) MRI measurements of tumor hypoxia with intrinsic susceptibility MRI measurements and to demonstrate the feasibility of translation of the technique for patients. Materials and Methods Preclinical studies in nine 786-0-R renal cell carcinoma (RCC) xenografts and prospective clinical studies in eight patients with RCC were performed. Longitudinal relaxation rate changes (∆R1) after 100% oxygen inhalation were quantified, reflecting the paramagnetic effect on tissue protons because of the presence of molecular oxygen. Native transverse relaxation rate (R2*) and oxygen-induced R2* change (∆R2*) were measured, reflecting presence of deoxygenated hemoglobin molecules. Median and voxel-wise values of ∆R1 were compared with values of R2* and ∆R2*. Tumor regions with dynamic contrast agent-enhanced MRI perfusion, refractory to signal change at OE MRI (referred to as perfused Oxy-R), were distinguished from perfused oxygen-enhancing (perfused Oxy-E) and nonperfused regions. R2* and ∆R2* values in each tumor subregion were compared by using one-way analysis of variance. Results Tumor-wise and voxel-wise ∆R1 and ∆R2* comparisons did not show correlative relationships. In xenografts, parcellation analysis revealed that perfused Oxy-R regions had faster native R2* (102.4 sec-1 vs 81.7 sec-1) and greater negative ∆R2* (-22.9 sec-1 vs -5.4 sec-1), compared with perfused Oxy-E and nonperfused subregions (all P < .001), respectively. Similar findings were present in human tumors (P < .001). Further, perfused Oxy-R helped identify tumor hypoxia, measured at pathologic analysis, in both xenografts (P = .002) and human tumors (P = .003). Conclusion Intrinsic susceptibility biomarkers provide cross validation of the OE MRI biomarker perfused Oxy-R. Consistent relationship to pathologic analyses was found in xenografts and human tumors, demonstrating biomarker translation. Published under a CC BY 4.0 license. Online supplemental material is available for this article.

A biomimetic tumor tissue phantom for validating diffusion-weighted MRI measurements.

PURPOSE: To develop a biomimetic tumor tissue phantom which more closely reflects water diffusion in biological tissue than previously used phantoms, and to evaluate the stability of the phantom and its potential as a tool for validating diffusion-weighted (DW) MRI measurements. METHODS: Coaxial-electrospraying was used to generate micron-sized hollow polymer spheres, which mimic cells. The bulk structure was immersed in water, providing a DW-MRI phantom whose apparent diffusion coefficient (ADC) and microstructural properties were evaluated over a period of 10 months. Independent characterization of the phantom's microstructure was performed using scanning electron microscopy (SEM). The repeatability of the construction process was investigated by generating a second phantom, which underwent high resolution synchrotron-CT as well as SEM and MR scans. RESULTS: ADC values were stable (coefficients of variation (CoVs) < 5%), and varied with diffusion time, with average values of 1.44 ± 0.03 µm2 /ms (Δ = 12 ms) and 1.20 ± 0.05 µm2 /ms (Δ = 45 ms). Microstructural parameters showed greater variability (CoVs up to 13%), with evidence of bias in sphere size estimates. Similar trends were observed in the second phantom. CONCLUSION: A novel biomimetic phantom has been developed and shown to be stable over 10 months. It is envisaged that such phantoms will be used for further investigation of microstructural models relevant to characterizing tumor tissue, and may also find application in evaluating acquisition protocols and comparing DW-MRI-derived biomarkers obtained from different scanners at different sites. Magn Reson Med 80:147-158, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

First-in-human technique translation of oxygen-enhanced MRI to an MR Linac system in patients with head and neck cancer.

BACKGROUND AND PURPOSE: Tumour hypoxia is prognostic in head and neck cancer (HNC), associated with poor loco-regional control, poor survival and treatment resistance. The advent of hybrid MRI - radiotherapy linear accelerator or 'MR Linac' systems - could permit imaging for treatment adaptation based on hypoxic status. We sought to develop oxygen-enhanced MRI (OE-MRI) in HNC and translate the technique onto an MR Linac system. MATERIALS AND METHODS: MRI sequences were developed in phantoms and 15 healthy participants. Next, 14 HNC patients (with 21 primary or local nodal tumours) were evaluated. Baseline tissue longitudinal relaxation time (T1) was measured alongside the change in 1/T1 (termed ΔR1) between air and oxygen gas breathing phases. We compared results from 1.5 T diagnostic MR and MR Linac systems. RESULTS: Baseline T1 had excellent repeatability in phantoms, healthy participants and patients on both systems. Cohort nasal concha oxygen-induced ΔR1 significantly increased (p < 0.0001) in healthy participants demonstrating OE-MRI feasibility. ΔR1 repeatability coefficients (RC) were 0.023-0.040 s-1 across both MR systems. The tumour ΔR1 RC was 0.013 s-1 and the within-subject coefficient of variation (wCV) was 25% on the diagnostic MR. Tumour ΔR1 RC was 0.020 s-1 and wCV was 33% on the MR Linac. ΔR1 magnitude and time-course trends were similar on both systems. CONCLUSION: We demonstrate first-in-human translation of volumetric, dynamic OE-MRI onto an MR Linac system, yielding repeatable hypoxia biomarkers. Data were equivalent on the diagnostic MR and MR Linac systems. OE-MRI has potential to guide future clinical trials of biology guided adaptive radiotherapy.

Coaxial electrospun biomimetic copolymer fibres for application in diffusion magnetic resonance imaging.

Objective. The use of diffusion magnetic resonance imaging (dMRI) opens the door to characterizing brain microstructure because water diffusion is anisotropic in axonal fibres in brain white matter and is sensitive to tissue microstructural changes. As dMRI becomes more sophisticated and microstructurally informative, it has become increasingly important to use a reference object (usually called an imaging phantom) for validation of dMRI. This study aims to develop axon-mimicking physical phantoms from biocopolymers and assess their feasibility for validating dMRI measurements.Approach. We employed a simple and one-step method-coaxial electrospinning-to prepare axon-mimicking hollow microfibres from polycaprolactone-b-polyethylene glycol (PCL-b-PEG) and poly(D, L-lactide-co-glycolic) acid (PLGA), and used them as building elements to create axon-mimicking phantoms. Electrospinning was firstly conducted using two types of PCL-b-PEG and two types of PLGA with different molecular weights in various solvents, with different polymer concentrations, for determining their spinnability. Polymer/solvent concentration combinations with good fibre spinnability were used as the shell material in the following co-electrospinning process in which the polyethylene oxide polymer was used as the core material. Following the microstructural characterization of both electrospun and co-electrospun fibres using optical and electron microscopy, two prototype phantoms were constructed from co-electrospun anisotropic hollow microfibres after inserting them into water-filled test tubes.Main results. Hollow microfibres that mimic the axon microstructure were successfully prepared from the appropriate core and shell material combinations. dMRI measurements of two phantoms on a 7 tesla (T) pre-clinical scanner revealed that diffusivity and anisotropy measurements are in the range of brain white matter.Significance. This feasibility study showed that co-electrospun PCL-b-PEG and PLGA microfibre-based axon-mimicking phantoms could be used in the validation of dMRI methods which seek to characterize white matter microstructure.

Image Contrast, Image Pre-Processing, and T1 Mapping Affect MRI Radiomic Feature Repeatability in Patients with Colorectal Cancer Liver Metastases.

Imaging biomarkers require technical, biological, and clinical validation to be translated into robust tools in research or clinical settings. This study contributes to the technical validation of radiomic features from magnetic resonance imaging (MRI) by evaluating the repeatability of features from four MR sequences: pre-contrast T1- and T2-weighted images, pre-contrast quantitative T1 maps (qT1), and contrast-enhanced T1-weighted images. Fifty-one patients with colorectal cancer liver metastases were scanned twice, up to 7 days apart. Repeatability was quantified using the intraclass correlation coefficient (ICC) and repeatability coefficient (RC), and the impact of non-Gaussian feature distributions and image normalisation was evaluated. Most radiomic features had non-Gaussian distributions, but Box-Cox transformations enabled ICCs and RCs to be calculated appropriately for an average of 97% of features across sequences. ICCs ranged from 0.30 to 0.99, with volume and other shape features tending to be most repeatable; volume ICC > 0.98 for all sequences. 19% of features from non-normalised images exhibited significantly different ICCs in pair-wise sequence comparisons. Normalisation tended to increase ICCs for pre-contrast T1- and T2-weighted images, and decrease ICCs for qT1 maps. RCs tended to vary more between sequences than ICCs, showing that evaluations of feature performance depend on the chosen metric. This work suggests that feature-specific repeatability, from specific combinations of MR sequence and pre-processing steps, should be evaluated to select robust radiomic features as biomarkers in specific studies. In addition, as different repeatability metrics can provide different insights into a specific feature, consideration of the appropriate metric should be taken in a study-specific context.

Detection of early changes in the post-radiosurgery vestibular schwannoma microenvironment using multinuclear MRI.

Stereotactic radiosurgery (SRS) is an established, effective therapy against vestibular schwannoma (VS). The mechanisms of tumour response are, however, unknown and in this study we sought to evaluate changes in the irradiated VS tumour microenvironment through a multinuclear MRI approach. Five patients with growing sporadic VS underwent a multi-timepoint comprehensive MRI protocol, which included diffusion tensor imaging (DTI), dynamic contrast-enhanced (DCE) MRI and a spiral 23Na-MRI acquisition for total sodium concentration (TSC) quantification. Post-treatment voxelwise changes in TSC, DTI metrics and DCE-MRI derived microvascular biomarkers (Ktrans, ve and vp) were evaluated and compared against pre-treatment values. Changes in tumour TSC and microvascular parameters were observable as early as 2 weeks post-treatment, preceding changes in structural imaging. At 6 months post-treatment there were significant voxelwise increases in tumour TSC (p < 0.001) and mean diffusivity (p < 0.001, repeated-measures ANOVA) with marked decreases in tumour microvascular parameters (p < 0.001, repeated-measures ANOVA). This study presents the first in vivo evaluation of alterations in the VS tumour microenvironment following SRS, demonstrating that changes in tumour sodium homeostasis and microvascular parameters can be imaged as early as 2 weeks following treatment. Future studies should seek to investigate these clinically relevant MRI metrics as early biomarkers of SRS response.

Co-electrospraying of tumour cell mimicking hollow polymeric microspheres for diffusion magnetic resonance imaging.

Diffusion magnetic resonance imaging (dMRI) is considered as a useful tool to study solid tumours. However, the interpretation of dMRI signal and validation of quantitative measurements of is challenging. One way to address these challenges is by using a standard reference material that can mimic tumour cell microstructure. There is a growing interest in using hollow polymeric microspheres, mainly prepared by multiple steps, as mimics of cells in healthy and diseased tissue. The present work reports on tumour cell-mimicking materials composed of hollow microspheres for application as a standard material in dMRI. These microspheres were prepared via one-step co-electrospraying process. The shell material was poly(d,l-lactic-co-glycolic acid) (PLGA) polymers with different molecule weights and/or ratios of glycolic acid-to-lactic, while the core was polyethylene glycol (PEG) or ethylene glycol. The resultant co-electrosprayed products were characterised by optical microscopy, scanning electron microscopy (SEM) and synchrotron X-ray micro-CT. These products were found to have variable structures and morphologies, e.g. from spherical particles with/without surface hole, through beaded fibres to smooth fibres, which mainly depend on PLGA composition and core materials. Only the shell material of PLGA polymer with ester terminated, Mw 50,000-75,000 g mol-1, and lactide:glycolide 85:15 formed hollow microspheres via the co-electrospraying process using the core material of 8 wt% PEG/chloroform as the core. A water-filled test object (or phantom) was designed and constructed from samples of the material generated from co-electrosprayed PLGA microspheres and tested on a 7 T MRI scanner. The preliminary MRI results provide evidence that hollow PLGA microspheres can restrict/hinder water diffusion as cells do in tumour tissue, implying that the phantom may be suitable for use as a quantitative validation and calibration tool for dMRI.

Ground Truth for Diffusion MRI in Cancer: A Model-Based Investigation of a Novel Tissue-Mimetic Material.

This work presents preliminary results on the development, characterisation, and use of a novel physical phantom designed as a simple mimic of tumour cellular structure, for diffusion-weighted magnetic resonance imaging (DW-MRI) applications. The phantom consists of a collection of roughly spherical, micron-sized core-shell polymer 'cells', providing a system whose ground truth microstructural properties can be determined and compared with those obtained from modelling the DW-MRI signal. A two-compartment analytic model combining restricted diffusion inside a sphere with hindered extracellular diffusion was initially investigated through Monte Carlo diffusion simulations, allowing a comparison between analytic and simulated signals. The model was then fitted to DW-MRI data acquired from the phantom over a range of gradient strengths and diffusion times, yielding estimates of 'cell' size, intracellular volume fraction and the free diffusion coefficient. An initial assessment of the accuracy and precision of these estimates is provided, using independent scanning electron microscope measurements and bootstrap-style simulations. Such phantoms may be useful for testing microstructural models relevant to the characterisation of tumour tissue.