Inter- and Intra-Patient Repeatability of Radiomic Features from Multiparametric Whole-Body MRI in Patients with Metastatic Prostate Cancer.

(1) Background: We assessed the test-re-test repeatability of radiomics in metastatic castration-resistant prostate cancer (mCPRC) bone disease on whole-body diffusion-weighted (DWI) and T1-weighted Dixon MRI. (2) Methods: In 10 mCRPC patients, 1.5 T MRI, including DWI and T1-weighted gradient-echo Dixon sequences, was performed twice on the same day. Apparent diffusion coefficient (ADC) and relative fat-fraction-percentage (rFF%) maps were calculated. Per study, up to 10 target bone metastases were manually delineated on DWI and Dixon images. All 106 radiomic features included in the Pyradiomics toolbox were derived for each target volume from the ADC and rFF% maps. To account for inter- and intra-patient measurement repeatability, the log-transformed individual target measurements were fitted to a hierarchical model, represented as a Bayesian network. Repeatability measurements, including the intraclass correlation coefficient (ICC), were derived. Feature ICCs were compared with mean ADC and rFF ICCs. (3) Results: A total of 65 DWI and 47 rFF% targets were analysed. There was no significant bias for any features. Pairwise correlation revealed fifteen ADC and fourteen rFF% feature sub-groups, without specific patterns between feature classes. The median intra-patient ICC was generally higher than the inter-patient ICC. Features that describe extremes in voxel values (minimum, maximum, range, skewness, and kurtosis) showed generally lower ICCs. Several mostly shape-based texture features were identified, which showed high inter- and intra-patient ICCs when compared with the mean ADC or mean rFF%, respectively. (4) Conclusions: Pyradiomics texture features of mCRPC bone metastases varied greatly in inter- and intra-patient repeatability. Several features demonstrated good repeatability, allowing for further exploration as diagnostic parameters in mCRPC bone disease.

Deep Learning for Delineation of the Spinal Canal in Whole-Body Diffusion-Weighted Imaging: Normalising Inter- and Intra-Patient Intensity Signal in Multi-Centre Datasets.

BACKGROUND: Whole-Body Diffusion-Weighted Imaging (WBDWI) is an established technique for staging and evaluating treatment response in patients with multiple myeloma (MM) and advanced prostate cancer (APC). However, WBDWI scans show inter- and intra-patient intensity signal variability. This variability poses challenges in accurately quantifying bone disease, tracking changes over follow-up scans, and developing automated tools for bone lesion delineation. Here, we propose a novel automated pipeline for inter-station, inter-scan image signal standardisation on WBDWI that utilizes robust segmentation of the spinal canal through deep learning. METHODS: We trained and validated a supervised 2D U-Net model to automatically delineate the spinal canal (both the spinal cord and surrounding cerebrospinal fluid, CSF) in an initial cohort of 40 patients who underwent WBDWI for treatment response evaluation (80 scans in total). Expert-validated contours were used as the target standard. The algorithm was further semi-quantitatively validated on four additional datasets (three internal, one external, 207 scans total) by comparing the distributions of average apparent diffusion coefficient (ADC) and volume of the spinal cord derived from a two-component Gaussian mixture model of segmented regions. Our pipeline subsequently standardises WBDWI signal intensity through two stages: (i) normalisation of signal between imaging stations within each patient through histogram equalisation of slices acquired on either side of the station gap, and (ii) inter-scan normalisation through histogram equalisation of the signal derived within segmented spinal canal regions. This approach was semi-quantitatively validated in all scans available to the study (N = 287). RESULTS: The test dice score, precision, and recall of the spinal canal segmentation model were all above 0.87 when compared to manual delineation. The average ADC for the spinal cord (1.7 × 10-3 mm2/s) showed no significant difference from the manual contours. Furthermore, no significant differences were found between the average ADC values of the spinal cord across the additional four datasets. The signal-normalised, high-b-value images were visualised using a fixed contrast window level and demonstrated qualitatively better signal homogeneity across scans than scans that were not signal-normalised. CONCLUSION: Our proposed intensity signal WBDWI normalisation pipeline successfully harmonises intensity values across multi-centre cohorts. The computational time required is less than 10 s, preserving contrast-to-noise and signal-to-noise ratios in axial diffusion-weighted images. Importantly, no changes to the clinical MRI protocol are expected, and there is no need for additional reference MRI data or follow-up scans.

The value of baseline 18F-sodium fluoride and 18F-choline PET activity for identifying responders to radium-223 treatment in castration-resistant prostate cancer bone metastases.

OBJECTIVES: To investigate whether baseline 18F-sodium fluoride (NaF) and 18F-choline PET activity is associated with metastatic castration-resistant prostate cancer (mCRPC) global and individual bone metastases' DWI MR imaging response to radium-223 treatment. METHODS: Thirty-six bone-only mCRPC patients were prospectively recruited from three centers. Whole-body (WB)-MRI with DWI and 18F-NaF and 18F-choline PET/CT were performed at therapy baseline and 8-week intervals. In each patient, bone disease median global (g)ADC change between baseline and follow-up was calculated. Additionally, up to five bone target lesions per patient were delineated and individual median ADC change recorded. An ADC increase > 30% defined response per-patient and per-lesion. For the same targets, baseline 18F-NaF and 18F-choline PET SUVmax were recorded. Mean SUVmax across patient targets was correlated with gADC change and lesion SUVmax with per-lesion ADC change. RESULTS: A total of 133 lesions in 36 patients (14 responders) were analyzed. 18F-NaF PET per-patient mean SUVmax was significantly higher in responders (median = 56.0 versus 38.7 in non-responders; p = 0.008), with positive correlation between SUVmax and gADC increase (rho = 0.42; p = 0.015). A 48.7 SUVmax threshold identified responders with 77% sensitivity and 75% specificity. Baseline 18F-NaF PET per-lesion SUVmax was higher in responding metastases (median = 51.6 versus 31.8 in non-responding metastases; p = 0.001), with positive correlation between baseline lesion SUVmax and ADC increase (rho = 0.39; p < 0.001). A 36.8 SUVmax threshold yielded 72% sensitivity and 63% specificity. No significant association was found between baseline 18F-choline PET SUVmax and ADC response on a per-patient (p = 0.164) or per-lesion basis (p = 0.921). CONCLUSION: 18F-NaF PET baseline SUVmax of target mCRPC bone disease showed significant association with response to radium-223 defined by ADC change. CLINICAL RELEVANCE STATEMENT: 18F-sodium fluoride PET/CT baseline maximum SUV of castration-resistant prostate cancer bone metastases could be used as a predictive biomarker for response to radium-223 therapy. KEY POINTS: • 18F-sodium fluoride PET baseline SUVmax of castration-resistant prostate cancer bone metastases showed significant association with response to radium-223. • Baseline 18F-sodium fluoride PET can improve patient selection for radium-223 therapy. • Change in whole-body DWI parameters can be used for response correlation with baseline 18F-sodium fluoride PET SUVmax in castration-resistant prostate cancer bone metastases.

Deep Learning Framework with Multi-Head Dilated Encoders for Enhanced Segmentation of Cervical Cancer on Multiparametric Magnetic Resonance Imaging.

T2-weighted magnetic resonance imaging (MRI) and diffusion-weighted imaging (DWI) are essential components of cervical cancer diagnosis. However, combining these channels for the training of deep learning models is challenging due to image misalignment. Here, we propose a novel multi-head framework that uses dilated convolutions and shared residual connections for the separate encoding of multiparametric MRI images. We employ a residual U-Net model as a baseline, and perform a series of architectural experiments to evaluate the tumor segmentation performance based on multiparametric input channels and different feature encoding configurations. All experiments were performed on a cohort of 207 patients with locally advanced cervical cancer. Our proposed multi-head model using separate dilated encoding for T2W MRI and combined b1000 DWI and apparent diffusion coefficient (ADC) maps achieved the best median Dice similarity coefficient (DSC) score, 0.823 (confidence interval (CI), 0.595-0.797), outperforming the conventional multi-channel model, DSC 0.788 (95% CI, 0.568-0.776), although the difference was not statistically significant (p > 0.05). We investigated channel sensitivity using 3D GRAD-CAM and channel dropout, and highlighted the critical importance of T2W and ADC channels for accurate tumor segmentation. However, our results showed that b1000 DWI had a minor impact on the overall segmentation performance. We demonstrated that the use of separate dilated feature extractors and independent contextual learning improved the model's ability to reduce the boundary effects and distortion of DWI, leading to improved segmentation performance. Our findings could have significant implications for the development of robust and generalizable models that can extend to other multi-modal segmentation applications.

Radium-223 in metastatic castration-resistant prostate cancer: whole-body diffusion-weighted magnetic resonance imaging scanning to assess response.

BACKGROUND: Radium-223 is a bone-seeking, ɑ-emitting radionuclide used to treat men with bone metastases from castration-resistant prostate cancer. Sclerotic bone lesions cannot be evaluated using Response Evaluation Criteria in Solid Tumors. Therefore, imaging response biomarkers are needed. METHODS: We conducted a phase 2 randomized trial to assess disease response to radium-223. Men with metastatic castration-resistant prostate cancer and bone metastases were randomly allocated to 55 or 88 kBq/kg radium-223 every 4 weeks for 6 cycles. Whole-body diffusion-weighted magnetic resonance imaging (DWI) was performed at baseline, at cycles 2 and 4, and after treatment. The primary endpoint was defined as a 30% increase in global median apparent diffusion coefficient. RESULTS: Disease response on DWI was seen in 14 of 36 evaluable patients (39%; 95% confidence interval = 23% to 56%), with marked interpatient and intrapatient heterogeneity of response. There was an association between prostate-specific antigen response and MRI response (odds ratio = 18.5, 95% confidence interval = 1.32 to 258, P = .013). Mean administered activity of radium-223 per cycle was not associated with global MRI response (P = .216) but was associated with DWI response using a 5-target-lesion evaluation (P = .007). In 26 of 36 (72%) patients, new bone metastases, not present at baseline, were seen on DWI scans during radium-223 treatment. CONCLUSIONS: DWI is useful for assessment of disease response in bone. Response to radium-223 is heterogeneous, both between patients and between different metastases in the same patient. New bone metastases appear during radium-223 treatment.The REASURE trial is registered under ISRCTN17805587.

Non-contrast CT synthesis using patch-based cycle-consistent generative adversarial network (Cycle-GAN) for radiomics and deep learning in the era of COVID-19.

Handcrafted and deep learning (DL) radiomics are popular techniques used to develop computed tomography (CT) imaging-based artificial intelligence models for COVID-19 research. However, contrast heterogeneity from real-world datasets may impair model performance. Contrast-homogenous datasets present a potential solution. We developed a 3D patch-based cycle-consistent generative adversarial network (cycle-GAN) to synthesize non-contrast images from contrast CTs, as a data homogenization tool. We used a multi-centre dataset of 2078 scans from 1,650 patients with COVID-19. Few studies have previously evaluated GAN-generated images with handcrafted radiomics, DL and human assessment tasks. We evaluated the performance of our cycle-GAN with these three approaches. In a modified Turing-test, human experts identified synthetic vs acquired images, with a false positive rate of 67% and Fleiss' Kappa 0.06, attesting to the photorealism of the synthetic images. However, on testing performance of machine learning classifiers with radiomic features, performance decreased with use of synthetic images. Marked percentage difference was noted in feature values between pre- and post-GAN non-contrast images. With DL classification, deterioration in performance was observed with synthetic images. Our results show that whilst GANs can produce images sufficient to pass human assessment, caution is advised before GAN-synthesized images are used in medical imaging applications.

Optimisation of b-values for the accurate estimation of the apparent diffusion coefficient (ADC) in whole-body diffusion-weighted MRI in patients with metastatic melanoma.

OBJECTIVE: To establish optimised diffusion weightings ('b-values') for acquisition of whole-body diffusion-weighted MRI (WB-DWI) for estimation of the apparent diffusion coefficient (ADC) in patients with metastatic melanoma (MM). Existing recommendations for WB-DWI have not been optimised for the tumour properties in MM; therefore, evaluation of acquisition parameters is essential before embarking on larger studies. METHODS: Retrospective clinical data and phantom experiments were used. Clinical data comprised 125 lesions from 14 examinations in 11 patients with multifocal MM, imaged before and/or after treatment with immunotherapy at a single institution. ADC estimates from these data were applied to a model to estimate the optimum b-value. A large non-diffusing phantom was used to assess eddy current-induced geometric distortion. RESULTS: Considering all tumour sites from pre- and post-treatment examinations together, metastases exhibited a large range of mean ADC values, [0.67-1.49] × 10-3 mm2/s, and the optimum high b-value (bhigh) for ADC estimation was 1100 (10th-90th percentile: 740-1790) s/mm2. At higher b-values, geometric distortion increased, and longer echo times were required, leading to reduced signal. CONCLUSIONS: Theoretical optimisation gave an optimum bhigh of 1100 (10th-90th percentile: 740-1790) s/mm2 for ADC estimation in MM, with the large range of optimum b-values reflecting the wide range of ADC values in these tumours. Geometric distortion and minimum echo time increase at higher b-values and are not included in the theoretical optimisation; bhigh in the range 750-1100 s/mm2 should be adopted to maintain acceptable image quality but performance should be evaluated for a specific scanner. KEY POINTS: • Theoretical optimisation gave an optimum high b-value of 1100 (10th-90th percentile: 740-1790) s/mm2 for ADC estimation in metastatic melanoma. • Considering geometric distortion and minimum echo time (TE), a b-value in the range 750-1100 s/mm2 is recommended. • Sites should evaluate the performance of specific scanners to assess the effect of geometric distortion and minimum TE.

Deep-learned estimation of uncertainty in measurements of apparent diffusion coefficient from whole-body diffusion-weighted MRI.

PURPOSE: To use deep learning to calculate the uncertainty in apparent diffusion coefficient (σADC) voxel-wise measurements to clinically impact the monitoring of treatment response and improve the quality of ADC maps. MATERIALS AND METHODS: We use a uniquely designed diffusion-weighted imaging (DWI) acquisition protocol that provides gold-standard measurements of σADC to train a deep learning model on two separate cohorts: 16 patients with prostate cancer and 28 patients with mesothelioma. Our network was trained with a novel cost function, which incorporates a perception metric and a b-value regularisation term, on ADC maps calculated by combinations of 2 or 3 b-values (e.g. 50/600/900, 50/900, 50/600, 600/900 s/mm2). We compare the accuracy of the deep-learning based approach for estimation of σADC with gold-standard measurements. RESULTS: The model accurately predicted the σADC for every b-value combination in both cohorts. Mean values of σADC within areas of active disease deviated from those measured by the gold-standard by 4.3% (range, 2.87-6.13%) for the prostate and 3.7% (range, 3.06-4.54%) for the mesothelioma cohort. We also showed that the model can easily be adapted for a different DWI protocol and field-of-view with only a few images (as little as a single patient) using transfer learning. CONCLUSION: Deep learning produces maps of σADC from standard clinical diffusion-weighted images (DWI) when 2 or more b-values are available.

Radiomic Features From Diffusion-Weighted MRI of Retroperitoneal Soft-Tissue Sarcomas Are Repeatable and Exhibit Change After Radiotherapy.

Background: Size-based assessments are inaccurate indicators of tumor response in soft-tissue sarcoma (STS), motivating the requirement for new response imaging biomarkers for this rare and heterogeneous disease. In this study, we assess the test-retest repeatability of radiomic features from MR diffusion-weighted imaging (DWI) and derived maps of apparent diffusion coefficient (ADC) in retroperitoneal STS and compare baseline repeatability with changes in radiomic features following radiotherapy (RT). Materials and Methods: Thirty patients with retroperitoneal STS received an MR examination prior to treatment, of whom 23/30 were investigated in our repeatability analysis having received repeat baseline examinations and 14/30 patients were investigated in our post-treatment analysis having received an MR examination after completing pre-operative RT. One hundred and seven radiomic features were extracted from the full manually delineated tumor region using PyRadiomics. Test-retest repeatability was assessed using an intraclass correlation coefficient (baseline ICC), and post-radiotherapy variance analysis (post-RT-IMS) was used to compare the change in radiomic feature value to baseline repeatability. Results: For the ADC maps and DWI images, 101 and 102 features demonstrated good baseline repeatability (baseline ICC > 0.85), respectively. Forty-three and 2 features demonstrated both good baseline repeatability and a high post-RT-IMS (>0.85), respectively. Pearson correlation between the baseline ICC and post-RT-IMS was weak (0.432 and 0.133, respectively). Conclusions: The ADC-based radiomic analysis shows better test-retest repeatability compared with features derived from DWI images in STS, and some of these features are sensitive to post-treatment change. However, good repeatability at baseline does not imply sensitivity to post-treatment change.

Virtual Biopsy in Soft Tissue Sarcoma. How Close Are We?

A shift in radiology to a data-driven specialty has been unlocked by synergistic developments in imaging biomarkers (IB) and computational science. This is advancing the capability to deliver "virtual biopsies" within oncology. The ability to non-invasively probe tumour biology both spatially and temporally would fulfil the potential of imaging to inform management of complex tumours; improving diagnostic accuracy, providing new insights into inter- and intra-tumoral heterogeneity and individualised treatment planning and monitoring. Soft tissue sarcomas (STS) are rare tumours of mesenchymal origin with over 150 histological subtypes and notorious heterogeneity. The combination of inter- and intra-tumoural heterogeneity and the rarity of the disease remain major barriers to effective treatments. We provide an overview of the process of successful IB development, the key imaging and computational advancements in STS including quantitative magnetic resonance imaging, radiomics and artificial intelligence, and the studies to date that have explored the potential biological surrogates to imaging metrics. We discuss the promising future directions of IBs in STS and illustrate how the routine clinical implementation of a virtual biopsy has the potential to revolutionise the management of this group of complex cancers and improve clinical outcomes.

Automatic Segmentation of Pelvic Cancers Using Deep Learning: State-of-the-Art Approaches and Challenges.

The recent rise of deep learning (DL) and its promising capabilities in capturing non-explicit detail from large datasets have attracted substantial research attention in the field of medical image processing. DL provides grounds for technological development of computer-aided diagnosis and segmentation in radiology and radiation oncology. Amongst the anatomical locations where recent auto-segmentation algorithms have been employed, the pelvis remains one of the most challenging due to large intra- and inter-patient soft-tissue variabilities. This review provides a comprehensive, non-systematic and clinically-oriented overview of 74 DL-based segmentation studies, published between January 2016 and December 2020, for bladder, prostate, cervical and rectal cancers on computed tomography (CT) and magnetic resonance imaging (MRI), highlighting the key findings, challenges and limitations.

Accelerating Whole-Body Diffusion-weighted MRI with Deep Learning-based Denoising Image Filters.

PURPOSE: To use deep learning to improve the image quality of subsampled images (number of acquisitions = 1 [NOA1]) to reduce whole-body diffusion-weighted MRI (WBDWI) acquisition times. MATERIALS AND METHODS: Both retrospective and prospective patient groups were used to develop a deep learning-based denoising image filter (DNIF) model. For initial model training and validation, 17 patients with metastatic prostate cancer with acquired WBDWI NOA1 and NOA9 images (acquisition period, 2015-2017) were retrospectively included. An additional 22 prospective patients with advanced prostate cancer, myeloma, and advanced breast cancer were used for model testing (2019), and the radiologic quality of DNIF-processed NOA1 (NOA1-DNIF) images were compared with NOA1 images and clinical NOA16 images by using a three-point Likert scale (good, average, or poor; statistical significance was calculated by using a Wilcoxon signed ranked test). The model was also retrained and tested in 28 patients with malignant pleural mesothelioma (MPM) who underwent lung MRI (2015-2017) to demonstrate feasibility in other body regions. RESULTS: The model visually improved the quality of NOA1 images in all test patients, with the majority of NOA1-DNIF and NOA16 images being graded as either "average" or "good" across all image-quality criteria. From validation data, the mean apparent diffusion coefficient (ADC) values within NOA1-DNIF images of bone disease deviated from those within NOA9 images by an average of 1.9% (range, 1.1%-2.6%). The model was also successfully applied in the context of MPM; the mean ADCs from NOA1-DNIF images of MPM deviated from those measured by using clinical-standard images (NOA12) by 3.7% (range, 0.2%-10.6%). CONCLUSION: Clinical-standard images were generated from subsampled images by using a DNIF.Keywords: Image Postprocessing, MR-Diffusion-weighted Imaging, Neural Networks, Oncology, Whole-Body Imaging, Supervised Learning, MR-Functional Imaging, Metastases, Prostate, Lung Supplemental material is available for this article. Published under a CC BY 4.0 license.

CT-Based Pelvic T1-Weighted MR Image Synthesis Using UNet, UNet++ and Cycle-Consistent Generative Adversarial Network (Cycle-GAN).

BACKGROUND: Computed tomography (CT) and magnetic resonance imaging (MRI) are the mainstay imaging modalities in radiotherapy planning. In MR-Linac treatment, manual annotation of organs-at-risk (OARs) and clinical volumes requires a significant clinician interaction and is a major challenge. Currently, there is a lack of available pre-annotated MRI data for training supervised segmentation algorithms. This study aimed to develop a deep learning (DL)-based framework to synthesize pelvic T1-weighted MRI from a pre-existing repository of clinical planning CTs. METHODS: MRI synthesis was performed using UNet++ and cycle-consistent generative adversarial network (Cycle-GAN), and the predictions were compared qualitatively and quantitatively against a baseline UNet model using pixel-wise and perceptual loss functions. Additionally, the Cycle-GAN predictions were evaluated through qualitative expert testing (4 radiologists), and a pelvic bone segmentation routine based on a UNet architecture was trained on synthetic MRI using CT-propagated contours and subsequently tested on real pelvic T1 weighted MRI scans. RESULTS: In our experiments, Cycle-GAN generated sharp images for all pelvic slices whilst UNet and UNet++ predictions suffered from poorer spatial resolution within deformable soft-tissues (e.g. bladder, bowel). Qualitative radiologist assessment showed inter-expert variabilities in the test scores; each of the four radiologists correctly identified images as acquired/synthetic with 67%, 100%, 86% and 94% accuracy. Unsupervised segmentation of pelvic bone on T1-weighted images was successful in a number of test cases. CONCLUSION: Pelvic MRI synthesis is a challenging task due to the absence of soft-tissue contrast on CT. Our study showed the potential of deep learning models for synthesizing realistic MR images from CT, and transferring cross-domain knowledge which may help to expand training datasets for 21 development of MR-only segmentation models.

Whole-body diffusion-weighted MRI in lymphoma-comparison of global apparent diffusion coefficient histogram parameters for differentiation of diseased nodes of lymphoma patients from normal lymph nodes of healthy individuals.

BACKGROUND: Morphologic features yield low diagnostic accuracy to distinguish between diseased and normal lymph nodes. The purpose of this study was to compare diseased lymphomatous and normal lymph nodes using global apparent diffusion coefficient (gADC) histogram parameters derived from whole-body diffusion-weighted MRI (WB-DWI). METHODS: 1.5 Tesla WB-DWI of 23 lymphoma patients and 20 healthy volunteers performed between 09/2010 and 07/2015 were retrospectively reviewed. All diseased nodal groups in the lymphoma cohort and all nodes visible on b900 images in healthy volunteers were segmented from neck to groin to generate a total diffusion volume (tDV). A connected component-labelling algorithm separated spatially distinct nodes. Mean, median, skewness, kurtosis, minimum, maximum, interquartile range (IQR), standard deviation (SD), 10th and 90th centile of the gADC distribution were derived from the tDV of each patient/volunteer and from spatially distinct nodes. gADC and regional nodal ADC parameters were compared between malignant and normal nodes using t-tests and ROC curve analyses. A P value ≤0.05 was deemed statistically significant. RESULTS: Mean, median, IQR, 10th and 90th centiles of gADC and regional nodal ADC values were significantly lower in diseased compared with normal lymph nodes. Skewness, kurtosis and tDV were significantly higher in lymphoma. The SD, min and max gADC showed no significant difference between the two groups (P>0.128). The diagnostic accuracies of gADC parameters by AUC from highest to lowest were: 10th centile, mean, median, 90th centile, skewness, kurtosis and IQR. A 10th centile gADC threshold of 0.68×10-3 mm2/s identified diseased lymphomatous nodes with 91% sensitivity and 95% specificity. CONCLUSIONS: WB-DWI derived gADC histogram parameters can distinguish between malignant lymph nodes of lymphoma patients and normal lymph nodes of healthy individuals.

Noninvasive MRI Native T1 Mapping Detects Response to MYCN-targeted Therapies in the Th-MYCN Model of Neuroblastoma.

Noninvasive early indicators of treatment response are crucial to the successful delivery of precision medicine in children with cancer. Neuroblastoma is a common solid tumor of young children that arises from anomalies in neural crest development. Therapeutic approaches aiming to destabilize MYCN protein, such as small-molecule inhibitors of Aurora A and mTOR, are currently being evaluated in early phase clinical trials in children with high-risk MYCN-driven disease, with limited ability to evaluate conventional pharmacodynamic biomarkers of response. T1 mapping is an MRI scan that measures the proton spin-lattice relaxation time T1. Using a multiparametric MRI-pathologic cross-correlative approach and computational pathology methodologies including a machine learning-based algorithm for the automatic detection and classification of neuroblasts, we show here that T1 mapping is sensitive to the rich histopathologic heterogeneity of neuroblastoma in the Th-MYCN transgenic model. Regions with high native T1 corresponded to regions dense in proliferative undifferentiated neuroblasts, whereas regions characterized by low T1 were rich in apoptotic or differentiating neuroblasts. Reductions in tumor-native T1 represented a sensitive biomarker of response to treatment-induced apoptosis with two MYCN-targeted small-molecule inhibitors, Aurora A kinase inhibitor alisertib (MLN8237) and mTOR inhibitor vistusertib (AZD2014). Overall, we demonstrate the potential of T1 mapping, a scan readily available on most clinical MRI scanners, to assess response to therapy and guide clinical trials for children with neuroblastoma. The study reinforces the potential role of MRI-based functional imaging in delivering precision medicine to children with neuroblastoma. SIGNIFICANCE: This study shows that MRI-based functional imaging can detect apoptotic responses to MYCN-targeted small-molecule inhibitors in a genetically engineered murine model of MYCN-driven neuroblastoma.

Noise-Corrected, Exponentially Weighted, Diffusion-Weighted MRI (niceDWI) Improves Image Signal Uniformity in Whole-Body Imaging of Metastatic Prostate Cancer.

Purpose: To characterize the voxel-wise uncertainties of Apparent Diffusion Coefficient (ADC) estimation from whole-body diffusion-weighted imaging (WBDWI). This enables the calculation of a new parametric map based on estimates of ADC and ADC uncertainty to improve WBDWI imaging standardization and interpretation: NoIse-Corrected Exponentially-weighted diffusion-weighted MRI (niceDWI). Methods: Three approaches to the joint modeling of voxel-wise ADC and ADC uncertainty (σADC) are evaluated: (i) direct weighted least squares (DWLS), (ii) iterative linear-weighted least-squares (IWLS), and (iii) smoothed IWLS (SIWLS). The statistical properties of these approaches in terms of ADC/σADC accuracy and precision is compared using Monte Carlo simulations. Our proposed post-processing methodology (niceDWI) is evaluated using an ice-water phantom, by comparing the contrast-to-noise ratio (CNR) with conventional exponentially-weighted DWI. We present the clinical feasibility of niceDWI in a pilot cohort of 16 patients with metastatic prostate cancer. Results: The statistical properties of ADC and σADC conformed closely to the theoretical predictions for DWLS, IWLS, and SIWLS fitting routines (a minor bias in parameter estimation is observed with DWLS). Ice-water phantom experiments demonstrated that a range of CNR could be generated using the niceDWI approach, and could improve CNR compared to conventional methods. We successfully implemented the niceDWI technique in our patient cohort, which visually improved the in-plane bias field compared with conventional WBDWI. Conclusions: Measurement of the statistical uncertainty in ADC estimation provides a practical way to standardize WBDWI across different scanners, by providing quantitative image signals that improve its reliability. Our proposed method can overcome inter-scanner and intra-scanner WBDWI signal variations that can confound image interpretation.

Supervised Machine-Learning Enables Segmentation and Evaluation of Heterogeneous Post-treatment Changes in Multi-Parametric MRI of Soft-Tissue Sarcoma.

Background: Multi-parametric MRI provides non-invasive methods for response assessment of soft-tissue sarcoma (STS) from non-surgical treatments. However, evaluation of MRI parameters over the whole tumor volume may not reveal the full extent of post-treatment changes as STS tumors are often highly heterogeneous, including cellular tumor, fat, necrosis, and cystic tissue compartments. In this pilot study, we investigate the use of machine-learning approaches to automatically delineate tissue compartments in STS, and use this approach to monitor post-radiotherapy changes. Methods: Eighteen patients with retroperitoneal sarcoma were imaged using multi-parametric MRI; 8/18 received a follow-up imaging study 2-4 weeks after pre-operative radiotherapy. Eight commonly-used supervised machine-learning techniques were optimized for classifying pixels into one of five tissue sub-types using an exhaustive cross-validation approach and expert-defined regions of interest as a gold standard. Final pixel classification was smoothed using a Markov Random Field (MRF) prior distribution on the final machine-learning models. Findings: 5/8 machine-learning techniques demonstrated high median cross-validation accuracies (82.2%, range 80.5-82.5%) with no significant difference between these five methods. One technique was selected (Naïve-Bayes) due to its relatively short training and class-prediction times (median 0.73 and 0.69 ms, respectively on a 3.5 GHz personal machine). When combined with the MRF-prior, this approach was successfully applied in all eight post-radiotherapy imaging studies and provided visualization and quantification of changes to independent STS sub-regions following radiotherapy for heterogeneous response assessment. Interpretation: Supervised machine-learning approaches to tissue classification in multi-parametric MRI of soft-tissue sarcomas provide quantitative evaluation of heterogeneous tissue changes following radiotherapy.

MRI Imaging of the Hemodynamic Vasculature of Neuroblastoma Predicts Response to Antiangiogenic Treatment.

Childhood neuroblastoma is a hypervascular tumor of neural origin, for which antiangiogenic drugs are currently being evaluated; however, predictive biomarkers of treatment response, crucial for successful delivery of precision therapeutics, are lacking. We describe an MRI-pathologic cross-correlative approach using intrinsic susceptibility (IS) and susceptibility contrast (SC) MRI to noninvasively map the vascular phenotype in neuroblastoma Th-MYCN transgenic mice treated with the vascular endothelial growth factor receptor inhibitor cediranib. We showed that the transverse MRI relaxation rate R 2* (second-1) and fractional blood volume (fBV, %) were sensitive imaging biomarkers of hemorrhage and vascular density, respectively, and were also predictive biomarkers of response to cediranib. Comparison with MRI and pathology from patients with MYCN-amplified neuroblastoma confirmed the high degree to which the Th-MYCN model vascular phenotype recapitulated that of the clinical phenotype, thereby supporting further evaluation of IS- and SC-MRI in the clinic. This study reinforces the potential role of functional MRI in delivering precision medicine to children with neuroblastoma. SIGNIFICANCE: This study shows that functional MRI predicts response to vascular-targeted therapy in a genetically engineered murine model of neuroblastoma.

Guidelines for Acquisition, Interpretation, and Reporting of Whole-Body MRI in Myeloma: Myeloma Response Assessment and Diagnosis System (MY-RADS).

Acknowledging the increasingly important role of whole-body MRI for directing patient care in myeloma, a multidisciplinary, international, and expert panel of radiologists, medical physicists, and hematologists with specific expertise in whole-body MRI in myeloma convened to discuss the technical performance standards, merits, and limitations of currently available imaging methods. Following guidance from the International Myeloma Working Group and the National Institute for Clinical Excellence in the United Kingdom, the Myeloma Response Assessment and Diagnosis System (or MY-RADS) imaging recommendations are designed to promote standardization and diminish variations in the acquisition, interpretation, and reporting of whole-body MRI in myeloma and allow response assessment. This consensus proposes a core clinical protocol for whole-body MRI and an extended protocol for advanced assessments. Published under a CC BY 4.0 license. Online supplemental material is available for this article.

Apparent diffusion coefficient of vertebral haemangiomas allows differentiation from malignant focal deposits in whole-body diffusion-weighted MRI.

OBJECTIVES: The aim of this study was to identify apparent diffusion coefficient (ADC) values for typical haemangiomas in the spine and to compare them with active malignant focal deposits. METHODS: This was a retrospective single-institution study. Whole-body magnetic resonance imaging (MRI) scans of 106 successive patients with active multiple myeloma, metastatic prostate or breast cancer were analysed. ADC values of typical vertebral haemangiomas and malignant focal deposits were recorded. RESULTS: The ADC of haemangiomas (72 ROIs, median ADC 1,085×10-6mm2s-1, interquartile range 927-1,295×10-6mm2s-1) was significantly higher than the ADC of malignant focal deposits (97 ROIs, median ADC 682×10-6mm2s-1, interquartile range 583-781×10-6mm2s-1) with a p-value < 10-6. Receiver operating characteristic (ROC) analysis produced an area under the curve of 0.93. An ADC threshold of 872×10-6mm2s-1 separated haemangiomas from malignant focal deposits with a sensitivity of 84.7 % and specificity of 91.8 %. CONCLUSIONS: ADC values of classical vertebral haemangiomas are significantly higher than malignant focal deposits. The high ADC of vertebral haemangiomas allows them to be distinguished visually and quantitatively from active sites of disease, which show restricted diffusion. KEY POINTS: • Whole-body diffusion-weighted MRI is becoming widely used in myeloma and bone metastases. • ADC values of vertebral haemangiomas are significantly higher than malignant focal deposits. • High ADCs of haemangiomas allows them to be distinguished from active disease.