The Image Biomarker Standardization Initiative: Standardized Convolutional Filters for Reproducible Radiomics and Enhanced Clinical Insights.

Filters are commonly used to enhance specific structures and patterns in images, such as vessels or peritumoral regions, to enable clinical insights beyond the visible image using radiomics. However, their lack of standardization restricts reproducibility and clinical translation of radiomics decision support tools. In this special report, teams of researchers who developed radiomics software participated in a three-phase study (September 2020 to December 2022) to establish a standardized set of filters. The first two phases focused on finding reference filtered images and reference feature values for commonly used convolutional filters: mean, Laplacian of Gaussian, Laws and Gabor kernels, separable and nonseparable wavelets (including decomposed forms), and Riesz transformations. In the first phase, 15 teams used digital phantoms to establish 33 reference filtered images of 36 filter configurations. In phase 2, 11 teams used a chest CT image to derive reference values for 323 of 396 features computed from filtered images using 22 filter and image processing configurations. Reference filtered images and feature values for Riesz transformations were not established. Reproducibility of standardized convolutional filters was validated on a public data set of multimodal imaging (CT, fluorodeoxyglucose PET, and T1-weighted MRI) in 51 patients with soft-tissue sarcoma. At validation, reproducibility of 486 features computed from filtered images using nine configurations × three imaging modalities was assessed using the lower bounds of 95% CIs of intraclass correlation coefficients. Out of 486 features, 458 were found to be reproducible across nine teams with lower bounds of 95% CIs of intraclass correlation coefficients greater than 0.75. In conclusion, eight filter types were standardized with reference filtered images and reference feature values for verifying and calibrating radiomics software packages. A web-based tool is available for compliance checking.

Transferability of radiomic signatures from experimental to human interstitial lung disease.

Background: Interstitial lung disease (ILD) defines a group of parenchymal lung disorders, characterized by fibrosis as their common final pathophysiological stage. To improve diagnosis and treatment of ILD, there is a need for repetitive non-invasive characterization of lung tissue by quantitative parameters. In this study, we investigated whether CT image patterns found in mice with bleomycin induced lung fibrosis can be translated as prognostic factors to human patients diagnosed with ILD. Methods: Bleomycin was used to induce lung fibrosis in mice (n_control = 36, n_experimental = 55). The patient cohort consisted of 98 systemic sclerosis (SSc) patients (n_ILD = 65). Radiomic features (n_histogram = 17, n_texture = 137) were extracted from microCT (mice) and HRCT (patients) images. Predictive performance of the models was evaluated with the area under the receiver-operating characteristic curve (AUC). First, predictive performance of individual features was examined and compared between murine and patient data sets. Second, multivariate models predicting ILD were trained on murine data and tested on patient data. Additionally, the models were reoptimized on patient data to reduce the influence of the domain shift on the performance scores. Results: Predictive power of individual features in terms of AUC was highly correlated between mice and patients (r = 0.86). A model based only on mean image intensity in the lung scored AUC = 0.921 ± 0.048 in mice and AUC = 0.774 (CI95% 0.677-0.859) in patients. The best radiomic model based on three radiomic features scored AUC = 0.994 ± 0.013 in mice and validated with AUC = 0.832 (CI95% 0.745-0.907) in patients. However, reoptimization of the model weights in the patient cohort allowed to increase the model's performance to AUC = 0.912 ± 0.058. Conclusion: Radiomic signatures of experimental ILD derived from microCT scans translated to HRCT of humans with SSc-ILD. We showed that the experimental model of BLM-induced ILD is a promising system to test radiomic models for later application and validation in human cohorts.

Improved Survival Prediction by Combining Radiological Imaging and S-100B Levels Into a Multivariate Model in Metastatic Melanoma Patients Treated With Immune Checkpoint Inhibition.

Purpose: We explored imaging and blood bio-markers for survival prediction in a cohort of patients with metastatic melanoma treated with immune checkpoint inhibition. Materials and Methods: 94 consecutive metastatic melanoma patients treated with immune checkpoint inhibition were included into this study. PET/CT imaging was available at baseline (Tp0), 3 months (Tp1) and 6 months (Tp2) after start of immunotherapy. Radiological response at Tp2 was evaluated using iRECIST. Total tumor burden (TB) at each time-point was measured and relative change of TB compared to baseline was calculated. LDH, CRP and S-100B were also analyzed. Cox proportional hazards model and logistic regression were used for survival analysis. Results: iRECIST at Tp2 was significantly associated with overall survival (OS) with C-index=0.68. TB at baseline was not associated with OS, whereas TB at Tp1 and Tp2 provided similar predictive power with C-index of 0.67 and 0.71, respectively. Appearance of new metastatic lesions during follow-up was an independent prognostic factor (C-index=0.73). Elevated LDH and S-100B ratios at Tp2 were significantly associated with worse OS: C-index=0.73 for LDH and 0.73 for S-100B. Correlation of LDH with TB was weak (r=0.34). A multivariate model including TB change, S-100B, and appearance of new lesions showed the best predictive performance with C-index=0.83. Conclusion: Our analysis shows only a weak correlation between LDH and TB. Additionally, baseline TB was not a prognostic factor in our cohort. A multivariate model combining early blood and imaging biomarkers achieved the best predictive power with regard to survival, outperforming iRECIST.

Gating has a negligible impact on dose delivered in MRI-guided online adaptive radiotherapy of prostate cancer.

BACKGROUND AND PURPOSE: MR-guided radiotherapy (MRgRT) allows real-time beam-gating to compensate for intra-fractional target position variations. This study investigates the dosimetric impact of beam-gating and the impact of PTV margin on prostate coverage for prostate cancer patients treated with online-adaptive MRgRT. MATERIALS AND METHODS: 20 consecutive prostate cancer patients were treated with online-adaptive MRgRT SBRT with 36.25 Gy in 5 fractions (PTV D95% ≥ 95% (N = 5) and PTV D95% ≥ 100% (N = 15)). Sagittal 2D cine MRIs were used for gating on the prostate with a 3 mm expansion as the gating window. We computed motion-compensated dose distributions for (i) all prostate positions during treatment (simulating non-gated treatments) and (ii) for prostate positions within the gating window (gated treatments). To evaluate the impact of PTV margin on prostate coverage, we simulated coverage with smaller margins than clinically applied both for gated and non-gated treatments. Motion-compensated fraction doses were accumulated and dose metrics were compared. RESULTS: We found a negligible dosimetric impact of beam-gating on prostate coverage (median of 0.00 Gy for both D95% and Dmean). For 18/20 patients, prostate coverage (D95% ≥ 100%) would have been ensured with a prostate-to-PTV margin of 3 mm, even without gating. The same was true for all but one fraction. CONCLUSION: Beam-gating has negligible dosimetric impact in online-adaptive MRgRT of prostate cancer. Accounting for motion, the clinically used prostate-to-PTV margin could potentially be reduced from 5 mm to 3 mm for 18/20 patients.

Computed tomography-based radiomics decodes prognostic and molecular differences in interstitial lung disease related to systemic sclerosis.

BACKGROUND: Radiomic features calculated from routine medical images show great potential for personalised medicine in cancer. Patients with systemic sclerosis (SSc), a rare, multiorgan autoimmune disorder, have a similarly poor prognosis due to interstitial lung disease (ILD). Here, our objectives were to explore computed tomography (CT)-based high-dimensional image analysis ("radiomics") for disease characterisation, risk stratification and relaying information on lung pathophysiology in SSc-ILD. METHODS: We investigated two independent, prospectively followed SSc-ILD cohorts (Zurich, derivation cohort, n=90; Oslo, validation cohort, n=66). For every subject, we defined 1355 robust radiomic features from standard-of-care CT images. We performed unsupervised clustering to identify and characterise imaging-based patient clusters. A clinically applicable prognostic quantitative radiomic risk score (qRISSc) for progression-free survival (PFS) was derived from radiomic profiles using supervised analysis. The biological basis of qRISSc was assessed in a cross-species approach by correlation with lung proteomic, histological and gene expression data derived from mice with bleomycin-induced lung fibrosis. RESULTS: Radiomic profiling identified two clinically and prognostically distinct SSc-ILD patient clusters. To evaluate the clinical applicability, we derived and externally validated a binary, quantitative radiomic risk score (qRISSc) composed of 26 features that accurately predicted PFS and significantly improved upon clinical risk stratification parameters in multivariable Cox regression analyses in the pooled cohorts. A high qRISSc score, which identifies patients at risk for progression, was reverse translatable from human to experimental ILD and correlated with fibrotic pathway activation. CONCLUSIONS: Radiomics-based risk stratification using routine CT images provides complementary phenotypic, clinical and prognostic information significantly impacting clinical decision making in SSc-ILD.

Comparison of robust to standardized CT radiomics models to predict overall survival for non-small cell lung cancer patients.

BACKGROUND: Radiomics is a promising tool for the identification of new prognostic biomarkers. Radiomic features can be affected by different scanning protocols, often present in retrospective and prospective clinical data. We compared a computed tomography (CT) radiomics model based on a large but highly heterogeneous multicentric image dataset with robust feature pre-selection to a model based on a smaller but standardized image dataset without pre-selection. MATERIALS AND METHODS: Primary tumor radiomics was extracted from pre-treatment CTs of IIIA/N2/IIIB NSCLC patients from a prospective Swiss multicentric randomized trial (npatient  = 124, ninstitution  = 14, SAKK 16/00) and a validation dataset (npatient  = 31, ninstitution  = 1). Four robustness studies investigating inter-observer delineation variation, motion, convolution kernel, and contrast were conducted to identify robust features using an intraclass correlation coefficient threshold >0.9. Two 12-months overall survival (OS) logistic regression models were trained: (a) on the entire multicentric heterogeneous dataset but with robust feature pre-selection (MCR) and (b) on a smaller standardized subset using all features (STD). Both models were validated on the validation dataset acquired with similar reconstruction parameters as the STD dataset. The model performances were compared using the DeLong test. RESULTS: In total, 113 stable features were identified (nshape  = 8, nintensity  = 0, ntexture  = 7, nwavelet  = 98). The convolution kernel had the strongest influence on the feature robustness (<20% stable features). The final models of MCR and STD consisted of one and two features respectively. Both features of the STD model were identified as non-robust. MCR did not show performance significantly different from STD on the validation cohort (AUC [95%CI] = 0.72 [0.48-0.95] and 0.79 [0.63-0.95], p = 0.59). CONCLUSION: Prognostic OS CT radiomics model for NSCLC based on a heterogeneous multicentric imaging dataset with robust feature pre-selection performed equally well as a model on a standardized dataset.

Radiomics, Tumor Volume, and Blood Biomarkers for Early Prediction of Pseudoprogression in Patients with Metastatic Melanoma Treated with Immune Checkpoint Inhibition.

PURPOSE: We assessed the predictive potential of positron emission tomography (PET)/CT-based radiomics, lesion volume, and routine blood markers for early differentiation of pseudoprogression from true progression at 3 months. EXPERIMENTAL DESIGN: 112 patients with metastatic melanoma treated with immune checkpoint inhibition were included in our study. Median follow-up duration was 22 months. 716 metastases were segmented individually on CT and 2[18F]fluoro-2-deoxy-D-glucose (FDG)-PET imaging at three timepoints: baseline (TP0), 3 months (TP1), and 6 months (TP2). Response was defined on a lesion-individual level (RECIST 1.1) and retrospectively correlated with FDG-PET/CT radiomic features and the blood markers LDH/S100. Seven multivariate prediction model classes were generated. RESULTS: Two-year (median) overall survival, progression-free survival, and immune progression-free survival were 69% (not reached), 24% (6 months), and 42% (16 months), respectively. At 3 months, 106 (16%) lesions had progressed, of which 30 (5%) were identified as pseudoprogression at 6 months. Patients with pseudoprogressive lesions and without true progressive lesions had a similar outcome to responding patients and a significantly better 2-year overall survival of 100% (30 months), compared with 15% (10 months) in patients with true progressions/without pseudoprogression (P = 0.002). Patients with mixed progressive/pseudoprogressive lesions were in between at 53% (25 months). The blood prediction model (LDH+S100) achieved an AUC = 0.71. Higher LDH/S100 values indicated a low chance of pseudoprogression. Volume-based models: AUC = 0.72 (TP1) and AUC = 0.80 (delta-volume between TP0/TP1). Radiomics models (including/excluding volume-related features): AUC = 0.79/0.78. Combined blood/volume model: AUC = 0.79. Combined blood/radiomics model (including volume-related features): AUC = 0.78. The combined blood/radiomics model (excluding volume-related features) performed best: AUC = 0.82. CONCLUSIONS: Noninvasive PET/CT-based radiomics, especially in combination with blood parameters, are promising biomarkers for early differentiation of pseudoprogression, potentially avoiding added toxicity or delayed treatment switch.

CT radiomics and PET radiomics: ready for clinical implementation?

INTRODUCTION: Today, rapid technical and clinical developments result in an increasing number of treatment options for oncological diseases. Thus, decision support systems are needed to offer the right treatment to the right patient. Imaging biomarkers hold great promise in patient-individual treatment guidance. Routinely performed for diagnosis and staging, imaging datasets are expected to hold more information than used in the clinical practice. Radiomics describes the extraction of a large number of meaningful quantitative features from medical images, such as computed tomography (CT) and positron emission tomography (PET). Due to the non-invasive nature and ability to capture 3D image-based heterogeneity, radiomic features are potential surrogate markers of the cancer phenotype. Several radiomic studies are published per day, owing to encouraging results of many radiomics-based patient outcome models. Despite this comparably large number of studies, radiomics is mainly studied in proof of principle concept. Hence, a translation of radiomics from a hot topic research field into an essential clinical decision-making tool is lacking, but of high clinical interest. EVIDENCE ACQUISITION: Herein, we present a literature review addressing the clinical evidence of CT and PET radiomics. An extensive literature review was conducted in PubMed, including papers on robustness and clinical applications. EVIDENCE SYNTHESIS: We summarize image-modality related influences on the robustness of radiomic features and provide an overview of clinical evidence reported in the literature. Today, more evidence has been provided for CT imaging, however, PET imaging offers the promise of direct imaging of biological processes and functions. We provide a summary of future research directions, which needs to be addressed in order to successfully introduce radiomics into clinical medicine. In comparison to CT, more focus should be directed towards harmonization of PET acquisition and reconstruction protocols, which is important for transferable modelling. CONCLUSIONS: Both CT and PET radiomics are promising pre-treatment and intra-treatment biomarkers for outcome prediction. Most studies are performed in retrospective setting, however their validation in prospective data collections is ongoing.

Incorporation of Dosimetric Gradients and Parotid Gland Migration Into Xerostomia Prediction.

Purpose: Due to the sharp gradients of intensity-modulated radiotherapy (IMRT) dose distributions, treatment uncertainties may induce substantial deviations from the planned dose during irradiation. Here, we investigate if the planned mean dose to parotid glands in combination with the dose gradient and information about anatomical changes during the treatment improves xerostomia prediction in head and neck cancer patients. Materials and methods: Eighty eight patients were retrospectively analyzed. Three features of the contralateral parotid gland were studied in terms of their association with the outcome, i.e., grade ≥ 2 (G2) xerostomia between 6 months and 2 years after radiotherapy (RT): planned mean dose (MD), average lateral dose gradient (GRADX), and parotid gland migration toward medial (PGM). PGM was estimated using daily megavoltage computed tomography (MVCT) images. Three logistic regression models where analyzed: based on (1) MD only, (2) MD and GRADX, and (3) MD, GRADX, and PGM. Additionally, the cohort was stratified based on the median value of GRADX, and a univariate analysis was performed to study the association of the MD with the outcome for patients in low- and high-GRADX domains. Results: The planned MD failed to recognize G2 xerostomia patients (AUC = 0.57). By adding the information of GRADX (second model), the model performance increased to AUC = 0.72. The addition of PGM (third model) led to further improvement in the recognition of the outcome (AUC = 0.79). Remarkably, xerostomia patients in the low-GRADX domain were successfully identified (AUC = 0.88) by the MD alone. Conclusions: Our results indicate that GRADX and PGM, which together serve as a proxy of dosimetric changes, provide valuable information for xerostomia prediction.

Design and Selection of Machine Learning Methods Using Radiomics and Dosiomics for Normal Tissue Complication Probability Modeling of Xerostomia.

PURPOSE: The purpose of this study is to investigate whether machine learning with dosiomic, radiomic, and demographic features allows for xerostomia risk assessment more precise than normal tissue complication probability (NTCP) models based on the mean radiation dose to parotid glands. MATERIAL AND METHODS: A cohort of 153 head-and-neck cancer patients was used to model xerostomia at 0-6 months (early), 6-15 months (late), 15-24 months (long-term), and at any time (a longitudinal model) after radiotherapy. Predictive power of the features was evaluated by the area under the receiver operating characteristic curve (AUC) of univariate logistic regression models. The multivariate NTCP models were tuned and tested with single and nested cross-validation, respectively. We compared predictive performance of seven classification algorithms, six feature selection methods, and ten data cleaning/class balancing techniques using the Friedman test and the Nemenyi post hoc analysis. RESULTS: NTCP models based on the parotid mean dose failed to predict xerostomia (AUCs < 0.60). The most informative predictors were found for late and long-term xerostomia. Late xerostomia correlated with the contralateral dose gradient in the anterior-posterior (AUC = 0.72) and the right-left (AUC = 0.68) direction, whereas long-term xerostomia was associated with parotid volumes (AUCs > 0.85), dose gradients in the right-left (AUCs > 0.78), and the anterior-posterior (AUCs > 0.72) direction. Multivariate models of long-term xerostomia were typically based on the parotid volume, the parotid eccentricity, and the dose-volume histogram (DVH) spread with the generalization AUCs ranging from 0.74 to 0.88. On average, support vector machines and extra-trees were the top performing classifiers, whereas the algorithms based on logistic regression were the best choice for feature selection. We found no advantage in using data cleaning or class balancing methods. CONCLUSION: We demonstrated that incorporation of organ- and dose-shape descriptors is beneficial for xerostomia prediction in highly conformal radiotherapy treatments. Due to strong reliance on patient-specific, dose-independent factors, our results underscore the need for development of personalized data-driven risk profiles for NTCP models of xerostomia. The facilitated machine learning pipeline is described in detail and can serve as a valuable reference for future work in radiomic and dosiomic NTCP modeling.