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.

Image quality in whole-body MRI using the MY-RADS protocol in a prospective multi-centre multiple myeloma study.

BACKGROUND: The Myeloma Response Assessment and Diagnosis System (MY-RADS) guidelines establish a standardised acquisition and analysis pipeline for whole-body MRI (WB-MRI) in patients with myeloma. This is the first study to assess image quality in a multi-centre prospective trial using MY-RADS. METHODS: The cohort consisted of 121 examinations acquired across ten sites with a range of prior WB-MRI experience, three scanner manufacturers and two field strengths. Image quality was evaluated qualitatively by a radiologist and quantitatively using a semi-automated pipeline to quantify common artefacts and image quality issues. The intra- and inter-rater repeatability of qualitative and quantitative scoring was also assessed. RESULTS: Qualitative radiological scoring found that the image quality was generally good, with 94% of examinations rated as good or excellent and only one examination rated as non-diagnostic. There was a significant correlation between radiological and quantitative scoring for most measures, and intra- and inter-rater repeatability were generally good. When the quality of an overall examination was low, this was often due to low quality diffusion-weighted imaging (DWI), where signal to noise ratio (SNR), anterior thoracic signal loss and brain geometric distortion were found as significant predictors of examination quality. CONCLUSIONS: It is possible to successfully deliver a multi-centre WB-MRI study using the MY-RADS protocol involving scanners with a range of manufacturers, models and field strengths. Quantitative measures of image quality were developed and shown to be significantly correlated with radiological assessment. The SNR of DW images was identified as a significant factor affecting overall examination quality. TRIAL REGISTRATION: ClinicalTrials.gov, NCT03188172 , Registered on 15 June 2017. CRITICAL RELEVANCE STATEMENT: Good overall image quality, assessed both qualitatively and quantitatively, can be achieved in a multi-centre whole-body MRI study using the MY-RADS guidelines. KEY POINTS: • A prospective multi-centre WB-MRI study using MY-RADS can be successfully delivered. • Quantitative image quality metrics were developed and correlated with radiological assessment. • SNR in DWI was identified as a significant predictor of quality, allowing for rapid quality adjustment.

Repeatability of quantitative individual lesion and total disease multiparametric whole-body MRI measurements in prostate cancer bone metastases.

OBJECTIVES: To assess the repeatability of quantitative multiparametric whole-body MRI (mpWB-MRI) parameters in advanced prostate cancer (APC) bone metastases. METHODS: 1.5T MRI was performed twice on the same day in 10 APC patients. MpWB-MRI-included diffusion weighted imaging (DWI) and T1-weighted gradient-echo 2-point Dixon sequences. ADC and relative fat-fraction percentage (rFF%) maps were calculated, respectively. A radiologist delineated up to 10 target bone metastases per study. Means of ADC, b900 signal intensity(SI), normalised b900 SI, rFF% and maximum diameter (MD) for each target lesion and overall parameter averages across all targets per patient were recorded. The total disease volume (tDV in ml) was manually delineated on b900 images and mean global (g)ADC was derived. Bland-Altman analyses were performed with calculation of 95% repeatability coefficients (RC). RESULTS: Seventy-three individual targets (median MD 26 mm) were included. Lesion mean ADC RC was 12.5%, mean b900 SI RC 137%, normalised mean b900 SI RC 110%, rFF% RC 3.2 and target MD RC 5.5 mm (16.3%). Patient target lesion average mean ADC RC was 6.4%, b900 SI RC 104% and normalised mean b900 SI RC 39.6%. Target average rFF% RC was 1.8, average MD RC 1.3 mm (4.8%). tDV segmentation RC was 6.4% and mean gADC RC 5.3%. CONCLUSIONS: APC bone metastases' ADC, rFF% and maximum diameter, tDV and gADC show good repeatability. ADVANCES IN KNOWLEDGE: APC bone metastases' mean ADC and rFF% measurements of single lesions and global disease volumes are repeatable, supporting their potential role as quantitative biomarkers in metastatic bone disease.

Mechanical mapping of the prostate in vivo using Dynamic Instrumented Palpation; towards an in vivo strategy for cancer assessment.

A calibrated palpation sensor has been developed for making instrumented Digital Rectal Examinations (iDREs) with a view to assessing patients for prostate cancer. The instrument measures the dynamic stiffness of the palpable surface of the prostate, and has been trialled on 12 patients in vivo. The patients had been diagnosed with prostate cancer and were scheduled for radical prostatectomy. As far as possible, patients with asymmetric disease were chosen so as to give a variation in gland condition over the palpable surface. The device works by applying an oscillating pressure (force) to a flexible probe whose displacement into the tissue is also measured in order to yield a dynamic stiffness, the static stiffness being incidentally measured at the mean oscillatory force. The device was deployed mounted on the index finger of a urologist and measurements taken at 12-16 positions on each patient using light and firm pressure and palpation frequencies of 1 or 5 Hz. In parallel, conventional DRE assessments were made by a consultant urologist for cancer. After in vivo measurement, the glands were removed and examined histologically with each palpation point being classified as cancerous (C) or not (NC). The work has established the first measurements of static modulus of living prostate tissue to be: 26.8 (13.3) kPa for tissue affected by prostate cancer (C classification), and 24.8 kPa (11.9) for tissue unaffected by cancer (NC classification), values quoted as median (interquartile range). The dynamic properties were characterised by: dynamic modulus, 5.15 kPa (4.86) for the C classification and 4.61 kPa (3.08) for the NC classification and the time lag between force and displacement at 5 Hz palpation frequency, 0.0175 s (0.0078) for the C classification and 0.0186 s (0.0397) for the NC classification, values again quoted as median (interquartile range). With the limited set of features that could be generated, an Artificial Neural Network (ANN) classification yielded a sensitivity of 97%, negative predictive value of 86%, positive predictive value of 67% and accuracy of 70% but with relatively poor specificity (30%). Besides extending the feature set, there are a number of changes in probe design, probing strategy and in mechanics analysis, which are expected to improve the diagnostic capabilities of the method.

Locating and sizing tumor nodules in human prostate using instrumented probing - computational framework and experimental validation.

Detection of tumor nodules is key to early cancer diagnosis. This study investigates the potential of using the mechanical data, acquired from probing the prostate for detecting the existence, and, more importantly, characterizing the size and depth, from the posterior surface, of the prostate cancer (PCa) nodules. A computational approach is developed to quantify the uncertainty of nodule detectability and is based on identifying stiffness anomalies in the profiles of point force measurements across transverse sections of the prostate. The capability of the proposed method was assessed firstly using a 'training' dataset of in silico models including PCa nodules with random size, depth and location, followed by a clinical feasibility study, involving experimental data from 13 ex vivo prostates from patients who had undergone radical prostatatectomy. Promising levels of sensitivity and specificity were obtained for detecting the PCa nodules in a total of 44 prostate sections. This study has shown that the proposed methods could be a useful complementary tool to exisiting diagnostic methods of PCa. The future study will involve implementing the proposed measurement and detection strategies in vivo, with the help of a miniturized medical device.

Biomechanical response of lumbar intervertebral disc in daily sitting postures: a poroelastic finite element analysis.

This study aims to establish and validate a poroelastic L4-L5 finite element model to evaluate the effect of different sitting postures and their durations on the mechanical responses of the disc. During the sustained loading conditions, the height loss, fluid loss and von-Mises stress gradually increased, but the intradiscal pressure decreased. The varying rates of aforementioned parameters were more significant at the initial loading stage and less so at the end. The predicted values in the flexed sitting posture were significantly greater than other postures. The extended sitting posture caused an obvious von-Mises stress concentration in the posterior region of the inter-lamellar matrix. From the biomechanical perspective, prolonged sitting may pose a high risk of lumbar disc degeneration, and therefore adjusting the posture properly in the early stage of sitting time may be useful to mitigate that. Additionally, upright sitting is a safer posture, while flexed sitting posture is more harmful.

Biomechanical effect of age-related structural changes on cervical intervertebral disc: A finite element study.

Previous literature has investigated the biomechanical response of healthy and degenerative discs, but the biomechanical response of suboptimal healthy intervertebral discs received less attention. The purpose was to compare the biomechanical responses and risk of herniation of young healthy, suboptimal healthy, and degenerative intervertebral discs. A cervical spine model was established and validated using the finite element method. Suboptimal healthy, mildly, moderately, and severely degenerative disc models were developed. Disc height deformation, range of motion, intradiscal pressure, and von Mises stress in annulus fibrosus were analyzed by applying a moment of 4 Nm in flexion, extension, lateral bending, and axial rotation with 100 N compressive loads. Disc height deformation in young healthy, suboptimal healthy, mildly, moderately, and severely degenerative discs was 40%, 37%, 21%, 12%, and 8%, respectively. The decreasing order of the range of motion was young healthy spine > suboptimal healthy spine > mildly degenerative spine > moderately degenerative spine > severely degenerative spine. The mean stress of annulus ground substance in the suboptimal healthy disc was higher than in the young healthy disc. The mean stress of inter-lamellar matrix and annulus ground substance in moderately and severely degenerative discs was higher than in other discs. Age-related structural changes and degenerative changes increased the stiffness and reduced the elastic deformation of intervertebral discs. Decreased range of motion due to the effects of aging or degeneration on the intervertebral disc, may cause compensation of adjacent segments and lead to progressive degeneration of multiple segments. The effect of aging on the intervertebral disc increased the risk of annulus fibrosus damage from the biomechanical point of view. Moderately and severely degenerative discs may have a higher risk of herniation due to the higher risk of damage and layers separation of annulus fibrosus caused by increased stress in the annulus ground substance and inter-lamellar matrix.

Implementation of Whole-Body MRI (MY-RADS) within the OPTIMUM/MUKnine multi-centre clinical trial for patients with myeloma.

BACKGROUND: Whole-body (WB) MRI, which includes diffusion-weighted imaging (DWI) and T1-w Dixon, permits sensitive detection of marrow disease in addition to qualitative and quantitative measurements of disease and response to treatment of bone marrow. We report on the first study to embed standardised WB-MRI within a prospective, multi-centre myeloma clinical trial (IMAGIMM trial, sub-study of OPTIMUM/MUKnine) to explore the use of WB-MRI to detect minimal residual disease after treatment. METHODS: The standardised MY-RADS WB-MRI protocol was set up on a local 1.5 T scanner. An imaging manual describing the MR protocol, quality assurance/control procedures and data transfer was produced and provided to sites. For non-identical scanners (different vendor or magnet strength), site visits from our physics team were organised to support protocol optimisation. The site qualification process included review of phantom and volunteer data acquired at each site and a teleconference to brief the multidisciplinary team. Image quality of initial patients at each site was assessed. RESULTS: WB-MRI was successfully set up at 12 UK sites involving 3 vendor systems and two field strengths. Four main protocols (1.5 T Siemens, 3 T Siemens, 1.5 T Philips and 3 T GE scanners) were generated. Scanner limitations (hardware and software) and scanning time constraint required protocol modifications for 4 sites. Nevertheless, shared methodology and imaging protocols enabled other centres to obtain images suitable for qualitative and quantitative analysis. CONCLUSIONS: Standardised WB-MRI protocols can be implemented and supported in prospective multi-centre clinical trials. Trial registration NCT03188172 clinicaltrials.gov; registration date 15th June 2017 https://clinicaltrials.gov/ct2/show/study/NCT03188172.

Identification of tumor nodule in soft tissue: An inverse finite-element framework based on mechanical characterization.

Identification and characterization of nodules in soft tissue, including their size, shape, and location, provide a basis for tumor identification. This study proposes an inverse finite-element (FE) based computational framework, for characterizing the size of examined tissue sample and detecting the presence of embedded tumor nodules using instrumented palpation, without a priori anatomical knowledge. The inverse analysis was applied to a model system, the human prostate, and was based on the reaction forces which can be obtained by trans-rectal mechanical probing and those from an equivalent FE model, which was optimized iteratively, by minimizing an error function between the two cases, toward the target solution. The tumor nodule can be identified through its influence on the stress state of the prostate. The effectiveness of the proposed method was further verified using a realistic prostate model reconstructed from magnetic resonance (MR) images. The results show the proposed framework to be capable of characterizing the key geometrical indices of the prostate and identifying the presence of cancerous nodules. Therefore, it has potential, when combined with instrumented palpation, for primary diagnosis of prostate cancer, and, potentially, solid tumors in other types of soft tissue.