Reproducible Radiomics Features from Multi-MRI-Scanner Test-Retest-Study: Influence on Performance and Generalizability of Models.

BACKGROUND: Radiomics models trained on data from one center typically show a decline of performance when applied to data from external centers, hindering their introduction into large-scale clinical practice. Current expert recommendations suggest to use only reproducible radiomics features isolated by multiscanner test-retest experiments, which might help to overcome the problem of limited generalizability to external data. PURPOSE: To evaluate the influence of using only a subset of robust radiomics features, defined in a prior in vivo multi-MRI-scanner test-retest-study, on the performance and generalizability of radiomics models. STUDY TYPE: Retrospective. POPULATION: Patients with monoclonal plasma cell disorders. Training set (117 MRIs from center 1); internal test set (42 MRIs from center 1); external test set (143 MRIs from center 2-8). FIELD STRENGTH/SEQUENCE: 1.5T and 3.0T; T1-weighted turbo spin echo. ASSESSMENT: The task for the radiomics models was to predict plasma cell infiltration, determined by bone marrow biopsy, noninvasively from MRI. Radiomics machine learning models, including linear regressor, support vector regressor (SVR), and random forest regressor (RFR), were trained on data from center 1, using either all radiomics features, or using only reproducible radiomics features. Models were tested on an internal (center 1) and a multicentric external data set (center 2-8). STATISTICAL TESTS: Pearson correlation coefficient r and mean absolute error (MAE) between predicted and actual plasma cell infiltration. Fisher's z-transformation, Wilcoxon signed-rank test, Wilcoxon rank-sum test; significance level P < 0.05. RESULTS: When using only reproducible features compared with all features, the performance of the SVR on the external test set significantly improved (r = 0.43 vs. r = 0.18 and MAE = 22.6 vs. MAE = 28.2). For the RFR, the performance on the external test set deteriorated when using only reproducible instead of all radiomics features (r = 0.33 vs. r = 0.44, P = 0.29 and MAE = 21.9 vs. MAE = 20.5, P = 0.10). CONCLUSION: Using only reproducible radiomics features improves the external performance of some, but not all machine learning models, and did not automatically lead to an improvement of the external performance of the overall best radiomics model. TECHNICAL EFFICACY: Stage 2.

Identification of focal lesion characteristics in MRI which indicate presence of corresponding osteolytic lesion in CT in patients with multiple myeloma.

PURPOSE: The presence of bone marrow focal lesions and osteolytic lesions in patients with multiple myeloma (MM) is of high prognostic significance for their individual outcome. It is not known yet why some focal lesions seen in MRI, reflecting localized bone marrow infiltration of myeloma cells, remain non-lytic, whereas others are associated with destruction of mineralized bone. In this study, we analyzed MRI characteristics of manually segmented focal lesions in MM patients to identify possible features that might discriminate lytic and non-lytic lesions. METHOD: The initial cohort included a total of 140 patients with different stages of MM who had undergone both whole-body MRI and whole-body low-dose CT within 30 days, and of which 29 satisfied the inclusion criteria for this study. Focal lesions in MRI and corresponding osteolytic areas in CT were segmented manually. Analysis of the lesions included volume, location and first order texture features analysis. RESULTS: There were significantly more lytic lesions in the axial skeleton than in the appendicular skeleton (p = 0.037). Out of 926 focal lesions in the axial skeleton seen on MRI, 544 (59.3 %) were osteolytic. Analysis of volume and first order texture features showed differences in texture and volume between focal lesions in MRI with and without local bone destruction in CT, but these findings were not statistically significant. CONCLUSIONS: Neither morphological imaging characteristics like size and location nor first order texture features could predict whether focal lesions seen in MRI would exhibit corresponding bone destruction in CT. Studies performing biopsies of such lesions are ongoing.

Same-day repeatability and Between-Sequence reproducibility of Mean ADC in PI-RADS lesions.

PURPOSE: This study aimed to assess repeatability after repositioning (inter-scan), intra-rater, inter-rater and inter-sequence variability of mean apparent diffusion coefficient (ADC) measurements in MRI-detected prostate lesions. METHOD: Forty-three patients with suspicion for prostate cancer were included and received a clinical prostate bi-/multiparametric MRI examination with repeat scans of the T2-weighted and two DWI-weighted sequences (ssEPI and rsEPI). Two raters (R1 and R2) performed single-slice, 2D regions of interest (2D-ROIs) and 3D-segmentation-ROIs (3D-ROIs). Mean bias, corresponding limits of agreement (LoA), mean absolute difference, within-subject coefficient of variation (CoV) and repeatability/reproducibility coefficient (RC/RDC) were calculated. Bradley & Blackwood test was used for variance comparison. Linear mixed models (LMM) were used to account for multiple lesions per patient. RESULTS: Inter-scan repeatability, intra-rater and inter-sequence reproducibility analysis of ADC showed no significant bias. 3D-ROIs demonstrated significantly less variability than 2D-ROIs (p < 0.01). Inter-rater comparison demonstrated small significant systematic bias of 57 × 10-6 mm2/s for 3D-ROIs (p < 0.001). Intra-rater RC, with the lowest variation, was 145 and 189 × 10-6 mm2/s for 3D- and 2D-ROIs, respectively. For 3D-ROIs of ssEPI, RCs and RDCs were 190-198 × 10-6 mm2/s for inter-scan, inter-rater and inter-sequence variation. No significant differences were found for inter-scan, inter-rater and inter-sequence variability. CONCLUSIONS: In a single-scanner setting, single-slice ADC measurements showed considerable variation, which may be lowered using 3D-ROIs. For 3D-ROIs, we propose a cut-off of âˆ¼ 200 × 10-6 mm2/s for differences introduced by repositioning, rater or sequence effects. The results suggest that follow-up measurements should be possible by different raters or sequences.

(Smoldering) multiple myeloma: mismatch between tumor load estimated from bone marrow biopsy at iliac crest and tumor load shown by MRI.

In multiple myeloma and its precursor stages, precise quantification of tumor load is of high importance for diagnosis, risk assessment, and therapy response evaluation. Both whole-body MRI, which allows to investigate the complete bone marrow of a patient, and bone marrow biopsy, which is commonly used to assess the histologic and genetic status, are relevant methods for tumor load assessment in multiple myeloma. We report on a series of striking mismatches between the plasma cell infiltration estimating the tumor load from unguided biopsies of the bone marrow at the posterior iliac crest and the tumor load assessment from whole-body MRI.

Prediction of Bone Marrow Biopsy Results From MRI in Multiple Myeloma Patients Using Deep Learning and Radiomics.

OBJECTIVES: In multiple myeloma and its precursor stages, plasma cell infiltration (PCI) and cytogenetic aberrations are important for staging, risk stratification, and response assessment. However, invasive bone marrow (BM) biopsies cannot be performed frequently and multifocally to assess the spatially heterogenous tumor tissue. Therefore, the goal of this study was to establish an automated framework to predict local BM biopsy results from magnetic resonance imaging (MRI). MATERIALS AND METHODS: This retrospective multicentric study used data from center 1 for algorithm training and internal testing, and data from center 2 to 8 for external testing. An nnU-Net was trained for automated segmentation of pelvic BM from T1-weighted whole-body MRI. Radiomics features were extracted from these segmentations, and random forest models were trained to predict PCI and the presence or absence of cytogenetic aberrations. Pearson correlation coefficient and the area under the receiver operating characteristic were used to evaluate the prediction performance for PCI and cytogenetic aberrations, respectively. RESULTS: A total of 672 MRIs from 512 patients (median age, 61 years; interquartile range, 53-67 years; 307 men) from 8 centers and 370 corresponding BM biopsies were included. The predicted PCI from the best model was significantly correlated ( P ≤ 0.01) to the actual PCI from biopsy in all internal and external test sets (internal test set: r = 0.71 [0.51, 0.83]; center 2, high-quality test set: r = 0.45 [0.12, 0.69]; center 2, other test set: r = 0.30 [0.07, 0.49]; multicenter test set: r = 0.57 [0.30, 0.76]). The areas under the receiver operating characteristic of the prediction models for the different cytogenetic aberrations ranged from 0.57 to 0.76 for the internal test set, but no model generalized well to all 3 external test sets. CONCLUSIONS: The automated image analysis framework established in this study allows for noninvasive prediction of a surrogate parameter for PCI, which is significantly correlated to the actual PCI from BM biopsy.

Test-retest, inter- and intra-rater reproducibility of size measurements of focal bone marrow lesions in MRI in patients with multiple myeloma.

OBJECTIVE: To investigate the reproducibility of size measurements of focal bone marrow lesions (FL) in MRI in patients with monoclonal plasma cell disorders under variation of patient positioning and observer. METHODS: A data set from a prospective test-retest study was used, in which 37 patients with a total of 140 FL had undergone 2 MRI scans with identical parameters after patient repositioning. Two readers measured long and short axis diameter on the initial scan in T1 weighted, T2 weighted short tau inversion recovery and diffusion-weighted imaging sequences. The first reader additionally measured FL on the retest-scan. The Bland-Altman method was used to assess limits of agreement (LoA), and the frequencies of absolute size changes were calculated. RESULTS: In the simple test-retest experiment with one identical reader, a deviation of ≥1 mm / ≥2 mm / ≥3 mm for the long axis diameter in T1 weighted images was observed in 66% / 25% / 8% of cases. When comparing measurements of one reader on the first scan to the measurement of the other reader on the retest scan, a change of ≥1 mm / ≥3 mm / ≥5 mm for the long axis diameter in T1 weighted images was observed in 78% / 21% / 5% of cases. CONCLUSION: Small deviations in FL size are common and probably due to variation in patient positioning or inter-rater variability alone, without any actual biological change of the FL. Knowledge of the uncertainty associated with size measurements of FLs is critical for radiologists and oncologists when interpreting changes in FL size in clinical practice and in clinical trials. ADVANCES IN KNOWLEDGE: According to the MY-RADs criteria, size measurements of focal lesions in MRI are now of relevance for response assessment in patients with monoclonal plasma cell disorders.Size changes of 1 or 2 mm are frequently observed due to uncertainty of the measurement only, while the actual focal lesion has not undergone any biological change.Size changes of at least 6 mm or more in T1 weighted or T2 weighted short tau inversion recovery sequences occur in only 5% or less of cases when the focal lesion has not undergone any biological change.

Deep Learning for Automatic Bone Marrow Apparent Diffusion Coefficient Measurements From Whole-Body Magnetic Resonance Imaging in Patients With Multiple Myeloma: A Retrospective Multicenter Study.

OBJECTIVES: Diffusion-weighted magnetic resonance imaging (MRI) is increasingly important in patients with multiple myeloma (MM). The objective of this study was to train and test an algorithm for automatic pelvic bone marrow analysis from whole-body apparent diffusion coefficient (ADC) maps in patients with MM, which automatically segments pelvic bones and subsequently extracts objective, representative ADC measurements from each bone. MATERIALS AND METHODS: In this retrospective multicentric study, 180 MRIs from 54 patients were annotated (semi)manually and used to train an nnU-Net for automatic, individual segmentation of the right hip bone, the left hip bone, and the sacral bone. The quality of the automatic segmentation was evaluated on 15 manually segmented whole-body MRIs from 3 centers using the dice score. In 3 independent test sets from 3 centers, which comprised a total of 312 whole-body MRIs, agreement between automatically extracted mean ADC values from the nnU-Net segmentation and manual ADC measurements from 2 independent radiologists was evaluated. Bland-Altman plots were constructed, and absolute bias, relative bias to mean, limits of agreement, and coefficients of variation were calculated. In 56 patients with newly diagnosed MM who had undergone bone marrow biopsy, ADC measurements were correlated with biopsy results using Spearman correlation. RESULTS: The ADC-nnU-Net achieved automatic segmentations with mean dice scores of 0.92, 0.93, and 0.85 for the right pelvis, the left pelvis, and the sacral bone, whereas the interrater experiment gave mean dice scores of 0.86, 0.86, and 0.77, respectively. The agreement between radiologists' manual ADC measurements and automatic ADC measurements was as follows: the bias between the first reader and the automatic approach was 49 × 10 -6 mm 2 /s, 7 × 10 -6 mm 2 /s, and -58 × 10 -6 mm 2 /s, and the bias between the second reader and the automatic approach was 12 × 10 -6 mm 2 /s, 2 × 10 -6 mm 2 /s, and -66 × 10 -6 mm 2 /s for the right pelvis, the left pelvis, and the sacral bone, respectively. The bias between reader 1 and reader 2 was 40 × 10 -6 mm 2 /s, 8 × 10 -6 mm 2 /s, and 7 × 10 -6 mm 2 /s, and the mean absolute difference between manual readers was 84 × 10 -6 mm 2 /s, 65 × 10 -6 mm 2 /s, and 75 × 10 -6 mm 2 /s. Automatically extracted ADC values significantly correlated with bone marrow plasma cell infiltration ( R = 0.36, P = 0.007). CONCLUSIONS: In this study, a nnU-Net was trained that can automatically segment pelvic bone marrow from whole-body ADC maps in multicentric data sets with a quality comparable to manual segmentations. This approach allows automatic, objective bone marrow ADC measurements, which agree well with manual ADC measurements and can help to overcome interrater variability or nonrepresentative measurements. Automatically extracted ADC values significantly correlate with bone marrow plasma cell infiltration and might be of value for automatic staging, risk stratification, or therapy response assessment.

In Vivo Repeatability and Multiscanner Reproducibility of MRI Radiomics Features in Patients With Monoclonal Plasma Cell Disorders: A Prospective Bi-institutional Study.

OBJECTIVES: Despite the extensive number of publications in the field of radiomics, radiomics algorithms barely enter large-scale clinical application. Supposedly, the low external generalizability of radiomics models is one of the main reasons, which hinders the translation from research to clinical application. The objectives of this study were to investigate reproducibility of radiomics features (RFs) in vivo under variation of patient positioning, magnetic resonance imaging (MRI) sequence, and MRI scanners, and to identify a subgroup of RFs that shows acceptable reproducibility across all different acquisition scenarios. MATERIALS AND METHODS: Between November 30, 2020 and February 16, 2021, 55 patients with monoclonal plasma cell disorders were included in this prospective, bi-institutional, single-vendor study. Participants underwent one reference scan at a 1.5 T MRI scanner and several retest scans: once after simple repositioning, once with a second MRI protocol, once at another 1.5 T scanner, and once at a 3 T scanner. Radiomics feature from the bone marrow of the left hip bone were extracted, both from original scans and after different image normalizations. Intraclass correlation coefficient (ICC) was used to assess RF repeatability and reproducibility. RESULTS: Fifty-five participants (mean age, 59 ± 7 years; 36 men) were enrolled. For T1-weighted images after muscle normalization, in the simple test-retest experiment, 110 (37%) of 295 RFs showed an ICC ≥0.8: 54 (61%) of 89 first-order features (FOFs), 35 (95%) of 37 volume and shape features, and 21 (12%) of 169 texture features (TFs). When the retest was performed with different technical settings, even after muscle normalization, the number of FOF/TF with an ICC ≥0.8 declined to 58/13 for the second protocol, 29/7 for the second 1.5 T scanner, and 49/7 for the 3 T scanner, respectively. Twenty-five (28%) of the 89 FOFs and 6 (4%) of the 169 TFs from muscle-normalized T1-weighted images showed an ICC ≥0.8 throughout all repeatability and reproducibility experiments. CONCLUSIONS: In vivo, only few RFs are reproducible with different MRI sequences or different MRI scanners, even after application of a simple image normalization. Radiomics features selected by a repeatability experiment only are not necessarily suited to build radiomics models for multicenter clinical application. This study isolated a subset of RFs, which are robust to variations in MRI acquisition observed in scanners from 1 vendor, and therefore are candidates to build reproducible radiomics models for monoclonal plasma cell disorders for multicentric applications, at least when centers are equipped with scanners from this vendor.

Combining Deep Learning and Radiomics for Automated, Objective, Comprehensive Bone Marrow Characterization From Whole-Body MRI: A Multicentric Feasibility Study.

OBJECTIVES: Disseminated bone marrow (BM) involvement is frequent in multiple myeloma (MM). Whole-body magnetic resonance imaging (wb-MRI) enables to evaluate the whole BM. Reading of such whole-body scans is time-consuming, and yet radiologists can transfer only a small fraction of the information of the imaging data set to the report. This limits the influence that imaging can have on clinical decision-making and in research toward precision oncology. The objective of this feasibility study was to implement a concept for automatic, comprehensive characterization of the BM from wb-MRI, by automatic BM segmentation and subsequent radiomics analysis of 30 different BM spaces (BMS). MATERIALS AND METHODS: This retrospective multicentric pilot study used a total of 106 wb-MRI from 102 patients with (smoldering) MM from 8 centers. Fifty wb-MRI from center 1 were used for training of segmentation algorithms (nnU-Nets) and radiomics algorithms. Fifty-six wb-MRI from 8 centers, acquired with a variety of different MRI scanners and protocols, were used for independent testing. Manual segmentations of 2700 BMS from 90 wb-MRI were performed for training and testing of the segmentation algorithms. For each BMS, 296 radiomics features were calculated individually. Dice score was used to assess similarity between automatic segmentations and manual reference segmentations. RESULTS: The "multilabel nnU-Net" segmentation algorithm, which performs segmentation of 30 BMS and labels them individually, reached mean dice scores of 0.88 ± 0.06/0.87 ± 0.06/0.83 ± 0.11 in independent test sets from center 1/center 2/center 3-8 (interrater variability between radiologists, 0.88 ± 0.01). The subset from the multicenter, multivendor test set (center 3-8) that was of high imaging quality was segmented with high precision (mean dice score, 0.87), comparable to the internal test data from center 1. The radiomic BM phenotype consisting of 8880 descriptive parameters per patient, which result from calculation of 296 radiomics features for each of the 30 BMS, was calculated for all patients. Exemplary cases demonstrated connections between typical BM patterns in MM and radiomic signatures of the respective BMS. In plausibility tests, predicted size and weight based on radiomics models of the radiomic BM phenotype significantly correlated with patients' actual size and weight ( P = 0.002 and P = 0.003, respectively). CONCLUSIONS: This pilot study demonstrates the feasibility of automatic, objective, comprehensive BM characterization from wb-MRI in multicentric data sets. This concept allows the extraction of high-dimensional phenotypes to capture the complexity of disseminated BM disorders from imaging. Further studies need to assess the clinical potential of this method for automatic staging, therapy response assessment, or prediction of biopsy results.

Whole-body magnetic resonance imaging plus serological follow-up for early identification of progression in smouldering myeloma patients to prevent development of end-organ damage.

The definition of multiple myeloma (MM) was updated in 2014, with the intent to enable earlier treatment and thereby avoid appearance of end-organ damage at progression from smouldering multiple myeloma (SMM) to MM. The purpose of this study was to investigate to which extent the development of end-organ damage at progression to MM was reduced under the updated guidelines. In this prospective observational cohort study (ClinicalTrials.gov Identifier: NCT01374412), between 2014 and 2020, 96 SMM patients prospectively underwent whole-body magnetic resonance imaging (wb-MRI) and serological follow-up at baseline and every 6 months thereafter. A total of 22 patients progressed into MM during follow-up, of which seven (32%) showed SLiM-criteria only but no end-organ damage. Four (57%) of the seven patients who progressed by SLiM-criteria only progressed with >1 focal lesion (FL) or a growing FL, and three (43%) due to serum free light-chain-ratio ≥100. Fifteen (68%) out of 22 patients who progressed still suffered from end-organ damage at progression. The updated disease definition reduced the proportion of SMM patients suffering from end-organ damage at progression to MM by one third. wb-MRI is an important tool for detection of SMM patients who progress to MM without end-organ damage.

Pseudoprospective Paraclinical Interaction of Radiology Residents With a Deep Learning System for Prostate Cancer Detection: Experience, Performance, and Identification of the Need for Intermittent Recalibration.

OBJECTIVES: The aim of this study was to estimate the prospective utility of a previously retrospectively validated convolutional neural network (CNN) for prostate cancer (PC) detection on prostate magnetic resonance imaging (MRI). MATERIALS AND METHODS: The biparametric (T2-weighted and diffusion-weighted) portion of clinical multiparametric prostate MRI from consecutive men included between November 2019 and September 2020 was fully automatically and individually analyzed by a CNN briefly after image acquisition (pseudoprospective design). Radiology residents performed 2 research Prostate Imaging Reporting and Data System (PI-RADS) assessments of the multiparametric dataset independent from clinical reporting (paraclinical design) before and after review of the CNN results and completed a survey. Presence of clinically significant PC was determined by the presence of an International Society of Urological Pathology grade 2 or higher PC on combined targeted and extended systematic transperineal MRI/transrectal ultrasound fusion biopsy. Sensitivities and specificities on a patient and prostate sextant basis were compared using the McNemar test and compared with the receiver operating characteristic (ROC) curve of CNN. Survey results were summarized as absolute counts and percentages. RESULTS: A total of 201 men were included. The CNN achieved an ROC area under the curve of 0.77 on a patient basis. Using PI-RADS ≥3-emulating probability threshold (c3), CNN had a patient-based sensitivity of 81.8% and specificity of 54.8%, not statistically different from the current clinical routine PI-RADS ≥4 assessment at 90.9% and 54.8%, respectively ( P = 0.30/ P = 1.0). In general, residents achieved similar sensitivity and specificity before and after CNN review. On a prostate sextant basis, clinical assessment possessed the highest ROC area under the curve of 0.82, higher than CNN (AUC = 0.76, P = 0.21) and significantly higher than resident performance before and after CNN review (AUC = 0.76 / 0.76, P ≤ 0.03). The resident survey indicated CNN to be helpful and clinically useful. CONCLUSIONS: Pseudoprospective paraclinical integration of fully automated CNN-based detection of suspicious lesions on prostate multiparametric MRI was demonstrated and showed good acceptance among residents, whereas no significant improvement in resident performance was found. General CNN performance was preserved despite an observed shift in CNN calibration, identifying the requirement for continuous quality control and recalibration.

Repeatability and Reproducibility of ADC Measurements and MRI Signal Intensity Measurements of Bone Marrow in Monoclonal Plasma Cell Disorders: A Prospective Bi-institutional Multiscanner, Multiprotocol Study.

BACKGROUND/OBJECTIVES: Apparent diffusion coefficient (ADC) and signal intensity (SI) measurements play an increasing role in magnetic resonance imaging (MRI) of monoclonal plasma cell disorders. The purpose of this study was to assess interrater variability, repeatability, and reproducibility of ADC and SI measurements from bone marrow (BM) under variation of MRI protocols and scanners. PATIENTS AND METHODS: Fifty-five patients with suspected or confirmed monoclonal plasma cell disorder were prospectively included in this institutional review board-approved study and underwent several measurements after the standard clinical whole-body MR scan, including repeated scan after repositioning, scan with a second MRI protocol, scan at a second 1.5 T scanner with a harmonized MRI protocol, and scan at a 3 T scanner. For T1-weighted, T2-weighted STIR, B800 images, and ADC maps, regions of interest were placed in the BM of the iliac crest and sacral bone, and in muscle tissue for image normalization. Bland-Altman plots were constructed, and absolute bias, relative bias to mean, limits of agreement, and coefficients of variation were calculated. RESULTS: Interrater variability and repeatability experiments showed a maximal relative bias of -0.077 and a maximal coefficient of variation of 16.2% for all sequences. Although the deviations at the second 1.5 T scanner with harmonized MRI protocol to the first 1.5 T scanner showed a maximal relative bias of 0.124 for all sequences, the variation of the MRI protocol and scan at the 3 T scanner led to large relative biases of up to -0.357 and -0.526, respectively. When comparing the 3 T scanner to the 1.5 T scanner, normalization to muscle reduced the bias of T1-weighted and T2-weighted sequences, but not of ADC maps. CONCLUSIONS: The MRI scanners with identical field strength and harmonized MRI protocols can provide relatively stable quantitative measurements of BM ADC and SI. Deviations in MRI field strength and MRI protocol should be avoided when applying ADC cutoff values, which were established at other scanners or when performing multicentric imaging trials.

[Potential of radiomics and artificial intelligence in myeloma imaging : Development of automatic, comprehensive, objective skeletal analyses from whole-body imaging data]. / Potenzial von Radiomics und künstlicher Intelligenz in der Myelombildgebung : Entwicklung automatischer, umfassender, objektiver Skelettanalysen aus Ganzkörperbildgebungsdaten.

CLINICAL/METHODICAL ISSUE: Multiple myeloma can affect the complete skeleton, which makes whole-body imaging necessary. With the current assessment of these complex datasets by radiologists, only a small part of the accessible information is assessed and reported. STANDARD RADIOLOGICAL METHODS: Depending on the question and availability, computed tomography (CT), magnetic resonance imaging (MRI), or positron emission tomography (PET) is performed and the results are then visually examined by radiologists. METHODOLOGICAL INNOVATIONS: A combination of automatic skeletal segmentation using artificial intelligence and subsequent radiomics analysis of each individual bone have the potential to provide automatic, comprehensive, and objective skeletal analyses. PERFORMANCE: A few automatic skeletal segmentation algorithms for CT already show promising results. In addition, first studies indicate correlations between radiomics features of bone and bone marrow with established disease markers and therapy response. ACHIEVEMENTS: Artificial intelligence (AI) and radiomics algorithms for automatic skeletal analysis from whole-body imaging are currently in an early phase of development.

Local Radiation Therapy Before and During Induction Delays Stem Cell Mobilization and Collection in Multiple Myeloma Patients.

In multiple myeloma, local radiation therapy (RT) of osseous lesions before peripheral blood stem cell (PBSC) mobilization is assumed to impair the PBSC mobilization and collection. However, the results of previously published studies are inconsistent and do not evaluate detailed metrics of RT and PBSC outcome parameters. In total, 352 patients undergoing PBSC mobilizations and RT in first-line treatment were evaluated. Patients were grouped into RT (n = 283) and no RT (n = 69) before PBSC mobilization. Except for the International Staging System score, both groups were homogeneous regarding the first diagnosis characteristics, first-line treatments, and response parameters. RT metrics (RT yes versus no, volume of irradiated hematopoietic bone marrow [BM], biologically equivalent doses in 2 Gy fractions [EQD2]) were analyzed for the following PBSC outcome parameters: achievement of the PBSC collection goal, CD34+ cell collection yield, duration of the mobilization phase, and number of leukapheresis (LP) sessions to reach the collection goal. No statistically significant differences in the percentage of collection failures to reach at least 3 sufficient PBSC transplants were identified comparing patients with (n = 32 [11%]) and without RT (n = 4 [6%]) before PBSC mobilization (P = .265). However, patients with RT before PBSC mobilization showed a significant prolongation of the PBSC mobilization (median 1 day, P =.026) and required a higher number of LP sessions to reach the collection goal (median 1 LP, P < .001) compared with patients who received RT after PBSC mobilization. Moreover, patients with RT before PBSC mobilization reached a significantly lower CD34+ cell collection result (mean 8.94 versus 9.81 × 106/kg body weight [bw], P = .002). No correlation was identified between the overall CD34+ cell yield and the volume of irradiated hematopoietic BM or EQD2, respectively. In the RT before PBSC mobilization group, patients who required more than 1 LP session to reach the PBSC collection goal after RT had a significantly higher percentage of radiated hematopoietic BM compared to those who required only 1 LP session (mean 9.7% versus 7.2%, P = .002). Overall, our study indicates a negative impact of RT on PBSC mobilization and collection. Apart from emergency settings, it might be beneficial to postpone RT to a post-PBSC collection time point. © 2021 American Society for Blood and Marrow Transplantation. Published by Elsevier Inc. All rights reserved.

Analyzing Longitudinal wb-MRI Data and Clinical Course in a Cohort of Former Smoldering Multiple Myeloma Patients: Connections between MRI Findings and Clinical Progression Patterns.

The purpose of this study was to analyze size and growth dynamics of focal lesions (FL) as well as to quantify diffuse infiltration (DI) in untreated smoldering multiple myeloma (SMM) patients and correlate those MRI features with timepoint and cause of progression. We investigated 199 whole-body magnetic resonance imaging (wb-MRI) scans originating from longitudinal imaging of 60 SMM patients and 39 computed tomography (CT) scans for corresponding osteolytic lesions (OL) in 17 patients. All FLs >5 mm were manually segmented to quantify volume and growth dynamics, and DI was scored, rating four compartments separately in T1- and fat-saturated T2-weighted images. The majority of patients with at least two FLs showed substantial spatial heterogeneity in growth dynamics. The volume of the largest FL (p = 0.001, c-index 0.72), the speed of growth of the fastest growing FL (p = 0.003, c-index 0.75), the DI score (DIS, p = 0.014, c-index 0.67), and its dynamic over time (DIS dynamic, p < 0.001, c-index 0.67) all significantly correlated with the time to progression. Size and growth dynamics of FLs correlated significantly with presence/appearance of OL in CT within 2 years after the respective MRI assessment (p = 0.016 and p = 0.022). DIS correlated with decrease of hemoglobin (p < 0.001). In conclusion, size and growth dynamics of FLs correlate with prognosis and local bone destruction. Connections between MRI findings and progression patterns (fast growing FL-OL; DIS-hemoglobin decrease) might enable more precise diagnostic and therapeutic approaches for SMM patients in the future.

Spatial Distribution of Focal Lesions in Whole-Body MRI and Influence of MRI Protocol on Staging in Patients with Smoldering Multiple Myeloma According to the New SLiM-CRAB-Criteria.

The purpose of this study was to assess how different MRI protocols (spinal vs. spinal plus pelvic vs. whole-body (wb)-MRI) affect staging in patients with smoldering multiple myeloma (SMM), according to the SLiM-CRAB-criterion '>1 focal lesion (FL) in MRI'. In this retrospective study, a baseline cohort of 147 SMM patients with wb-MRI at initial diagnosis was investigated, including prognostic data regarding development of CRAB-criteria. Fifty-two patients formed a follow-up cohort with a median of three wb-MRIs. The locations of all FLs were determined and it was calculated how staging decisions regarding the criterion '>1 FL in MRI' would have been made if only a limited anatomic area (spine vs. spine plus pelvis) would have been covered by the MRI protocol. Furthermore, subgroups of patients selected by different cutoff-protocol-combinations were compared regarding their prognosis for development of CRAB-criteria. With an MRI protocol limited to spine/spine plus pelvis, only 28%/64% of patients who actually had >1 FL in wb-MRI would have been rated correctly as having '>1 FL in MRI'. Fifty-four percent/36% of patients with exactly 1 FL in spine/spine plus pelvis revealed >1 FL when the entire wb-MRI was analyzed. During follow-up, four more patients developed >1 FL in wb-MRI; both limited MRI protocols would have detected only one of these four patients as having >1 FL at the correct timepoint. Having >1 FL in spine/in spine plus pelvis/in the whole body was associated with a 43%/57%/49% probability of developing CRAB-criteria within 2 years. Patients with >3 FL in spine plus pelvis and patients with >4 FL in the whole body had an 80% probability to develop CRAB-criteria within 2 years. MRI protocols limited to the spine or to spine plus pelvis lead to substantial underdiagnoses of patients who actually have >1 FL in wb-MRI at baseline and during follow-up, which influences staging and treatment decisions according to the current SLiM-CRAB criteria. However, given the spatial distribution of FLs and the analysis on clinical course of patients indicates that the cutoff for the number of FLs should be adopted according to the MRI protocol when using MRI for staging in SMM.

Towards quantitative imaging biomarkers of tumor dissemination: A multi-scale parametric modeling of multiple myeloma.

The advent of medical imaging and automatic image analysis is bringing the full quantitative assessment of lesions and tumor burden at every clinical examination within reach. This opens avenues for the development and testing of functional disease models, as well as their use in the clinical practice for personalized medicine. In this paper, we introduce a Bayesian statistical framework, based on mixed-effects models, to quantitatively test and learn functional disease models at different scales, on population longitudinal data. We also derive an effective mathematical model for the crossover between initially detected lesions and tumor dissemination, based on the Iwata-Kawasaki-Shigesada model. We finally propose to leverage this descriptive disease progression model into model-aware biomarkers for personalized risk-assessment, taking all available examinations and relevant covariates into account. As a use case, we study Multiple Myeloma, a disseminated plasma cell cancer, in which proper diagnostics is essential, to differentiate frequent precursor state without end-organ damage from the rapidly developing disease requiring therapy. After learning the best biological models for local lesion growth and global tumor burden evolution on clinical data, and computing corresponding population priors, we use individual model parameters as biomarkers, and can study them systematically for correlation with external covariates, such as sex or location of the lesion. On our cohort of 63 patients with smoldering Multiple Myeloma, we show that they perform substantially better than other radiological criteria, to predict progression into symptomatic Multiple Myeloma. Our study paves the way for modeling disease progression patterns for Multiple Myeloma, but also for other metastatic and disseminated tumor growth processes, and for analyzing large longitudinal image data sets acquired in oncological imaging. It shows the unprecedented potential of model-based biomarkers for better and more personalized treatment decisions and deserves being validated on larger cohorts to establish its role in clinical decision making.

Volumetry based biomarker speed of growth: Quantifying the change of total tumor volume in whole-body magnetic resonance imaging over time improves risk stratification of smoldering multiple myeloma patients.

The purpose of this study was to improve risk stratification of smoldering multiple myeloma patients, introducing new 3D-volumetry based imaging biomarkers derived from whole-body MRI. Two-hundred twenty whole-body MRIs from 63 patients with smoldering multiple myeloma were retrospectively analyzed and all focal lesions >5mm were manually segmented for volume quantification. The imaging biomarkers total tumor volume, speed of growth (development of the total tumor volume over time), number of focal lesions, development of the number of focal lesions over time and the recent imaging biomarker '>1 focal lesion' of the International Myeloma Working Group were compared, taking 2-year progression rate, sensitivity and false positive rate into account. Speed of growth, using a cutoff of 114mm3/month, was able to isolate a high-risk group with a 2-year progression rate of 82.5%. Additionally, it showed by far the highest sensitivity in this study and in comparison to other biomarkers in the literature, detecting 63.2% of patients who progress within 2 years. Furthermore, its false positive rate (8.7%) was much lower compared to the recent imaging biomarker '>1 focal lesion' of the International Myeloma Working Group. Therefore, speed of growth is the preferable imaging biomarker for risk stratification of smoldering multiple myeloma patients.