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

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

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

Developing and testing a robotic MRI/CT fusion biopsy technique using a purpose-built interventional phantom.

BACKGROUND: Magnetic resonance imaging (MRI) can be used to target tumour components in biopsy procedures, while the ability to precisely correlate histology and MRI signal is crucial for imaging biomarker validation. Robotic MRI/computed tomography (CT) fusion biopsy offers the potential for this without in-gantry biopsy, although requires development. METHODS: Test-retest T1 and T2 relaxation times, attenuation (Hounsfield units, HU), and biopsy core quality were prospectively assessed (January-December 2021) in a range of gelatin, agar, and mixed gelatin/agar solutions of differing concentrations on days 1 and 8 after manufacture. Suitable materials were chosen, and four biopsy phantoms were constructed with twelve spherical 1-3-cm diameter targets visible on MRI, but not on CT. A technical pipeline was developed, and intraoperator and interoperator reliability was tested in four operators performing a total of 96 biopsies. Statistical analysis included T1, T2, and HU repeatability using Bland-Altman analysis, Dice similarity coefficient (DSC), and intraoperator and interoperator reliability. RESULTS: T1, T2, and HU repeatability had 95% limits-of-agreement of 8.3%, 3.4%, and 17.9%, respectively. The phantom was highly reproducible, with DSC of 0.93 versus 0.92 for scanning the same or two different phantoms, respectively. Hit rate was 100% (96/96 targets), and all operators performed robotic biopsies using a single volumetric acquisition. The fastest procedure time was 32 min for all 12 targets. CONCLUSIONS: A reproducible biopsy phantom was developed, validated, and used to test robotic MRI/CT-fusion biopsy. The technique was highly accurate, reliable, and achievable in clinically acceptable timescales meaning it is suitable for clinical application

A framework for optimization of diffusion-weighted MRI protocols for large field-of-view abdominal-pelvic imaging in multicenter studies.

PURPOSE: To develop methods for optimization of diffusion-weighted MRI (DW-MRI) in the abdomen and pelvis on 1.5 T MR scanners from three manufacturers and assess repeatability of apparent diffusion coefficient (ADC) estimates in a temperature-controlled phantom and abdominal and pelvic organs in healthy volunteers. METHODS: Geometric distortion, ghosting, fat suppression, and repeatability and homogeneity of ADC estimates were assessed using phantoms and volunteers. Healthy volunteers (ten per scanner) were each scanned twice on the same scanner. One volunteer traveled to all three institutions in order to provide images for qualitative comparison. The common volunteer was excluded from quantitative analysis of the data from scanners 2 and 3 in order to ensure statistical independence, giving n = 10 on scanner 1 and n = 9 on scanners 2 and 3 for quantitative analysis. Repeatability and interscanner variation of ADC estimates in kidneys, liver, spleen, and uterus were assessed using within-patient coefficient of variation (wCV) and Kruskal-Wallis tests, respectively. RESULTS: The coefficient of variation of ADC estimates in the temperature-controlled phantom was 1%-4% for all scanners. Images of healthy volunteers from all scanners showed homogeneous fat suppression and no marked ghosting or geometric distortion. The wCV of ADC estimates was 2%-4% for kidneys, 3%-7% for liver, 6%-9% for spleen, and 7%-10% for uterus. ADC estimates in kidneys, spleen, and uterus showed no significant difference between scanners but a significant difference was observed in liver (p < 0.05). CONCLUSIONS: DW-MRI protocols can be optimized using simple phantom measurements to produce good quality images in the abdomen and pelvis at 1.5 T with repeatable quantitative measurements in a multicenter study

Development of a temperature-controlled phantom for magnetic resonance quality assurance of diffusion, dynamic, and relaxometry measurements.

Purpose Diffusion-weighted (DW) and dynamic contrast-enhanced magnetic resonance imaging (MRI) are increasingly applied for the assessment of functional tissue biomarkers for diagnosis, lesion characterization, or for monitoring of treatment response. However, these techniques are vulnerable to the influence of various factors, so there is a necessity for a standardized MR quality assurance procedure utilizing a phantom to facilitate the reliable estimation of repeatability of these quantitative biomarkers arising from technical factors (e.g., B1 variation) affecting acquisition on scanners of different vendors and field strengths. The purpose of this study is to present a novel phantom designed for use in quality assurance for multicenter trials, and the associated repeatability measurements of functional and quantitative imaging protocols across different MR vendors and field strengths.Methods A cylindrical acrylic phantom was manufactured containing 7 vials of polyvinylpyrrolidone (PVP) solutions of different concentrations, ranging from 0% (distilled water) to 25% w/w, to create a range of different MR contrast parameters. Temperature control was achieved by equilibration with ice-water. Repeated MR imaging measurements of the phantom were performed on four clinical scanners (two at 1.5 T, two at 3.0 T; two vendors) using the same scanning protocol to assess the long-term and short-term repeatability. The scanning protocol consisted of DW measurements, inversion recovery (IR) T1 measurements, multiecho T2 measurement, and dynamic T1-weighted sequence allowing multiple variable flip angle (VFA) estimation of T1 values over time. For each measurement, the corresponding calculated parameter maps were produced. On each calculated map, regions of interest (ROIs) were drawn within each vial and the median value of these voxels was assessed. For the dynamic data, the autocorrelation function and their variance were calculated; for the assessment of the repeatability, the coefficients of variation (CoV) were calculated.Results For both field strengths across the available vendors, the apparent diffusion coefficient (ADC) at 0 °C ranged from (1.12 ± 0.01) × 10(-3) mm(2)/s for pure water to (0.48 ± 0.02) × 10(-3) mm(2)/s for the 25% w/w PVP concentration, presenting a minor variability between the vendors and the field strengths. T2 and IR-T1 relaxation time results demonstrated variability between the field strengths and the vendors across the different acquisitions. Moreover, the T1 values derived from the VFA method exhibited a large variation compared with the IR-T1 values across all the scanners for all repeated measurements, although the calculation of the standard deviation of the VFA-T1 estimate across each ROI and the autocorrelation showed a stability of the signal for three scanners, with autocorrelation of the signal over the dynamic series revealing a periodic variation in one scanner. Finally, the ADC, the T2, and the IR-T1 values exhibited an excellent repeatability across the scanners, whereas for the dynamic data, the CoVs were higher.Conclusions The combination of a novel PVP phantom, with multiple compartments to give a physiologically relevant range of ADC and T1 values, together with ice-water as a temperature-controlled medium, allows reliable quality assurance measurements that can be used to measure agreement between MRI scanners, critical in multicenter functional and quantitative imaging studies

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

Implementing diffusion-weighted MRI for body imaging in prospective multicentre trials: current considerations and future perspectives.

For body imaging, diffusion-weighted MRI may be used for tumour detection, staging, prognostic information, assessing response and follow-up. Disease detection and staging involve qualitative, subjective assessment of images, whereas for prognosis, progression or response, quantitative evaluation of the apparent diffusion coefficient (ADC) is required. Validation and qualification of ADC in multicentre trials involves examination of i) technical performance to determine biomarker bias and reproducibility and ii) biological performance to interrogate a specific aspect of biology or to forecast outcome. Unfortunately, the variety of acquisition and analysis methodologies employed at different centres make ADC values non-comparable between them. This invalidates implementation in multicentre trials and limits utility of ADC as a biomarker. This article reviews the factors contributing to ADC variability in terms of data acquisition and analysis. Hardware and software considerations are discussed when implementing standardised protocols across multi-vendor platforms together with methods for quality assurance and quality control. Processes of data collection, archiving, curation, analysis, central reading and handling incidental findings are considered in the conduct of multicentre trials. Data protection and good clinical practice are essential prerequisites. Developing international consensus of procedures is critical to successful validation if ADC is to become a useful biomarker in oncology. KEY POINTS:• Standardised acquisition/analysis allows quantification of imaging biomarkers in multicentre trials. • Establishing "precision" of the measurement in the multicentre context is essential. • A repository with traceable data of known provenance promotes further research

Utility of Multi-Parametric Quantitative Magnetic Resonance Imaging for Characterization and Radiotherapy Response Assessment in Soft-Tissue Sarcomas and Correlation With Histopathology.

Purpose: To evaluate repeatability of quantitative multi-parametric MRI in retroperitoneal sarcomas, assess parameter changes with radiotherapy, and correlate pre-operative values with histopathological findings in the surgical specimens. Materials and Methods: Thirty patients with retroperitoneal sarcoma were imaged at baseline, of whom 27 also underwent a second baseline examination for repeatability assessment. 14/30 patients were treated with pre-operative radiotherapy and were imaged again after completing radiotherapy (50.4 Gy in 28 daily fractions, over 5.5 weeks). The following parameter estimates were assessed in the whole tumor volume at baseline and following radiotherapy: apparent diffusion coefficient (ADC), parameters of the intra-voxel incoherent motion model of diffusion-weighted MRI (D, f, D*), transverse relaxation rate, fat fraction, and enhancing fraction after gadolinium-based contrast injection. Correlation was evaluated between pre-operative quantitative parameters and histopathological assessments of cellularity and fat fraction in post-surgical specimens (ClinicalTrials.gov, registration number NCT01902667). Results: Upper and lower 95% limits of agreement were 7.1 and -6.6%, respectively for median ADC at baseline. Median ADC increased significantly post-radiotherapy. Pre-operative ADC and D were negatively correlated with cellularity (r = -0.42, p = 0.01, 95% confidence interval (CI) -0.22 to -0.59 for ADC; r = -0.45, p = 0.005, 95% CI -0.25 to -0.62 for D), and fat fraction from Dixon MRI showed strong correlation with histopathological assessment of fat fraction (r = 0.79, p = 10-7, 95% CI 0.69-0.86). Conclusion: Fat fraction on MRI corresponded to fat content on histology and therefore contributes to lesion characterization. Measurement repeatability was excellent for ADC; this parameter increased significantly post-radiotherapy even in disease categorized as stable by size criteria, and corresponded to cellularity on histology. ADC can be utilized for characterizing and assessing response in heterogeneous retroperitoneal sarcomas