Medical image processing and analysis (also known as Radiomics) is a rapidly growing discipline that maps digital medical images into quantitative data, with the end goal of generating imaging biomarkers as decision support tools for clinical practice. The use of imaging data from routine clinical work-up has tremendous potential in improving cancer care by heightening understanding of tumor biology and aiding in the implementation of precision medicine. As a noninvasive method of assessing the tumor and its microenvironment in their entirety, radiomics allows the evaluation and monitoring of tumor characteristics such as temporal and spatial heterogeneity. One can observe a rapid increase in the number of computational medical imaging publications—milestones that have highlighted the utility of imaging biomarkers in oncology. Nevertheless, the use of radiomics as clinical biomarkers still necessitates amelioration and standardization in order to achieve routine clinical adoption. This Review addresses the critical issues to ensure the proper development of radiomics as a biomarker and facilitate its implementation in clinical practice.

Radiomics extracts high-throughput quantitative data from medical images to contribute to precision medicine. Radiomic shape features have been shown to correlate with patient outcomes. However, how radiomic shape features vary in function of tumor complexity and tumor volume, as well as with method used for meshing and voxel resampling, remains unknown. The aims of this study are to create tumor models with varying degrees of complexity, or spiculatedness, and evaluate their relationship with quantitatively extracted shape features. Twenty-eight tumor models were mathematically created using spherical harmonics with the spiculatedness degreedbeing increased by increments of 3 (d = 11 tod = 92). Models were 3D printed with identical bases of 5 cm, imaged with a CT scanner with two different slice thicknesses, and semi-automatically delineated. Resampling of the resulting masks on a 1 × 1 × 1 mm³grid was performed, and the voxel size of each model was then calculated to eliminate volume differences. Four MATLAB-based algorithms (isosurface (M1), isosurface filter (M2), isosurface remeshing (M3), and boundary (M4)) were used to extract nine 3D features (Volume, Surface area, Surface-to-volume, Compactness1, Compactness2, Compactness3, Spherical Disproportion, Sphericity and Fractional Concavity). To quantify the impact of 3D printing, acquisition, segmentation and meshing, features were computed directly from the stereolithography (STL) file format that was used for 3D printing, and compared to those computed. Changes in feature values between 0.6 and 2 mm slice acquisitions were also compared. Spearman’s rank-order correlation coefficients were computed to determine the relationship of each shape feature with spiculatedness for each of the four meshing algorithms. Percent changes were calculated between shape features extracted from the original and resampled contoured images to evaluate the influence of spatial resampling. Finally, the percent change in shape features when the volume was changed from 25% to 150% of their original volume was quantified for three distinct tumor models and compared to the percent change observed when modifying the spiculatedness of the model fromd = 11 tod = 92. Values extracted using isosurface remeshing method are the closest to the STL reference ones, with mean differences less than 10.8% (Compactness2) for all features. Seven of the eight features had strong significant correlations with tumor model complexity irrespective of the meshing algorithm (r > 0.98, p < 10^-4), with fractional concavity having the lowest correlation coefficient (r = 0.83, p < 10^-4, M2). Comparisons of features extracted from the 0.6 and 2 mm slice thicknesses showed that mean differences were from 2.1% (Compactness3) to 12.7% (Compactness2) for the isosurface remeshing method. Resampling on a 1 × 1 × 1 mm³grid resulted in between 1.3% (Compactness3) to 9.5% (Fractional Concavity) mean changes in feature values. Compactness2, Compactness3, Spherical Disproportion, Sphericity and Fractional Concavity were the features least affected by volume changes. Compactness1 had a 90.4% change with volume, which was greater than the change between the least and most spiculated models. This is the first methodological study that directly demonstrates the relationship of tumor spiculatedness with radiomic shape features, that also produced 3D tumor models, which may serve as reference phantoms for future radiomic studies. Surface Area, Surface-to-volume, and Spherical Disproportion had direct relationships with spiculatedness while the three formulas for Compactness, Sphericity and Fractional Concavity had inverse relationships. The features Compactness2, Compactness3, Spherical Disproportion, and Sphericity should be prioritized as these have minimal variations with volume changes, slice thickness and resampling.

Background Combining radiotherapy (RT) with immuno-oncology (IO) therapy (IORT) may enhance IO-induced antitumor response. Quantitative imaging biomarkers can be used to provide prognosis, predict tumor response in a non-invasive fashion and improve patient selection for IORT. A biologically inspired CD8 T-cells-associated radiomics signature has been developed on previous cohorts. We evaluated here whether this CD8 radiomic signature is associated with lesion response, whether it may help to assess disease spatial heterogeneity for predicting outcomes of patients treated with IORT. We also evaluated differences between irradiated and non-irradiated lesions. Methods Clinical data from patients with advanced solid tumors in six independent clinical studies of IORT were investigated. Immunotherapy consisted of 4 different drugs (antiprogrammed death-ligand 1 or anticytotoxic T-lymphocyte-associated protein 4 in monotherapy). Most patients received stereotactic RT to one lesion. Irradiated and non-irradiated lesions were delineated from baseline and the first evaluation CT scans. Radiomic features were extracted from contrast-enhanced CT images and the CD8 radiomics signature was applied. A responding lesion was defined by a decrease in lesion size of at least 30%. Dispersion metrices of the radiomics signature were estimated to evaluate the impact of tumor heterogeneity in patient’s response. Results A total of 94 patients involving multiple lesions (100 irradiated and 189 non-irradiated lesions) were considered for a statistical interpretation. Lesions with high CD8 radiomics score at baseline were associated with significantly higher tumor response (area under the receiving operating characteristic curve (AUC)=0.63, p=0.0020). Entropy of the radiomics scores distribution on all lesions was shown to be associated with progression-free survival (HR=1.67, p=0.040), out-of-field abscopal response (AUC=0.70, p=0.014) and overall survival (HR=2.08, p=0.023), which remained significant in a multivariate analysis including clinical and biological variables. Conclusions These results enhance the predictive value of the biologically inspired CD8 radiomics score and suggests that tumor heterogeneity should be systematically considered in patients treated with IORT. This CD8 radiomics signature may help select patients who are most likely to benefit from IORT.

Background: The predictive power of novel biological markers for treatment response to immune checkpoint inhibitors (ICI) is still not satisfactory for the majority of patients with cancer. One should identify valid predictive markers in the peripheral blood, as this is easily available before and during treatment. The current interim analysis of patients of the ST-ICI cohort therefore focuses on the development and validation of a liquid immune profile-based signature (LIPS) to predict response of patients with metastatic cancer to ICI targeting the programmed cell death protein 1 (PD-1)/programmed cell death-ligand 1 (PD-L1) axis. Methods: A total of 104 patients were prospectively enrolled. 54 immune cell subsets were prospectively analyzed in patients’ peripheral blood by multicolor flow cytometry before treatment with ICI (pre-ICI; n=89), and after the first application of ICI (n=65). Pre-ICI, patients were randomly allocated to a training (n=56) and a validation cohort (n=33). Univariate Cox proportional hazards regression analysis and least absolute shrinkage and selection operator Cox model were used to create a predictive immune signature, which was also checked after the first ICI, to consider the dynamics of changes in the immune status. Results: Whole blood samples were provided by 89 patients pre-ICI and by 65 patients after the first ICI. We identified a LIPS which is based on five immune cell subtypes: CD14^high monocytes, CD8+/PD-1⁺ T cells, plasmacytoid dendritic cells, neutrophils, and CD3⁺/CD56⁺/CD16⁺ natural killer (NK)T cells. The signature achieved a high accuracy (C-index 0.74 vs 0.71) for predicting overall survival (OS) benefit in both the training and the validation cohort. In both cohorts, the low-risk group had significantly longer OS than the high-risk group (HR 0.26, 95% CI 0.12 to 0.56, p=0.00025; HR 0.30, 95% CI 0.10 to 0.91, p=0.024, respectively). Regarding the whole cohort, LIPS also predicted progression-free survival (PFS). The identified LIPS was not affected by clinicopathological features with the exception of brain metastases. NKT cells and neutrophils of the LIPS can be used as dynamic predictive biomarkers for OS and PFS after first administration of the ICI. Conclusion: Our study identified a predictive LIPS for survival of patients with cancer treated with PD-1/PD-L1 ICI, which is based on immune cell subsets in the peripheral whole blood. Trial registration number: NCT03453892.