Impact of Image Reconstruction Settings on Texture Features in 18F-FDG PET

Introduction to Radiomics

Radiomics is a rapidly evolving field of research concerned with the extraction of quantitative metrics-the so-called radiomic features-within medical images. Radiomic features capture tissue and lesion characteristics such as heterogeneity and shape and may, alone or in combination with demographic, histologic, genomic, or proteomic data, be used for clinical problem solving. The goal of this continuing education article is to provide an introduction to the field, covering the basic radiomics workflow: feature calculation and selection, dimensionality reduction, and data processing. Potential clinical applications in nuclear medicine that include PET radiomics-based prediction of treatment response and survival will be discussed. Current limitations of radiomics, such as sensitivity to acquisition parameter variations, and common pitfalls will also be covered.

FLT PET Radiomics for Response Prediction to Chemoradiation Therapy in Head and Neck Squamous Cell Cancer

Radiomics is an image analysis approach for extracting large amounts of quantitative information from medical images using a variety of computational methods. Our goal was to evaluate the utility of radiomic feature analysis from ¹⁸F-fluorothymidine positron emission tomography (FLT PET) obtained at baseline in prediction of treatment response in patients with head and neck cancer. Thirty patients with advanced-stage oropharyngeal or laryngeal cancer, treated with definitive chemoradiation therapy, underwent FLT PET imaging before treatment. In total, 377 radiomic features of FLT uptake and feature variants were extracted from volumes of interest; these features variants were defined by either the primary tumor or the total lesion burden, which consisted of the primary tumor and all FLT-avid nodes. Feature variants included normalized measurements of uptake, which were calculated by dividing lesion uptake values by the mean uptake value in the bone marrow. Feature reduction was performed using clustering to remove redundancy, leaving 172 representative features. Effects of these features on progression-free survival were modeled with Cox regression and P-values corrected for multiple comparisons. In total, 9 features were considered significant. Our results suggest that smaller, more homogenous lesions at baseline were associated with better prognosis. In addition, features extracted from total lesion burden had a higher concordance index than primary tumor features for 8 of the 9 significant features. Furthermore, total lesion burden features showed lower interobserver variability.

Application of radiomics for the prediction of HPV status for patients with head and neck cancers

PURPOSE

To perform radiomic analysis of primary tumors extracted from pretreatment contrast-enhanced computed tomography (CE-CT) images for patients with oropharyngeal cancers to identify discriminant features and construct an optimal classifier for the characterization and prediction of human papilloma virus (HPV) status.

MATERIALS AND METHODS

One hundred and eighty seven patients with oropharyngeal cancers with known HPV status (confirmed by immunohistochemistry-p16 protein testing) were retrospectively studied as follows: Group A: 95 patients (19HPV- and 76HPV+) from the MICAII grand challenge. Group B: 92 patients (52HPV- and 40HPV+) from our institution. Radiomic features (172) were extracted from pretreatment diagnostic CE-CT images of the gross tumor volume (GTV). Levene and Kolmogorov-Smirnov's tests with absolute biserial correlation (>0.48) were used to identify the discriminant features between the HPV+ and HPV- groups. The discriminant features were used to train and test eight different classifiers. Area under receiver operating characteristic (AUC), positive predictive and negative predictive values (PPV and NPV, respectively) were used to evaluate the performance of the classifiers. Principal component analysis (PCA) was applied on the discriminant feature set and seven PCs were used to train and test a generalized linear model (GLM) classifier.

RESULTS

Among 172 radiomic features only 12 radiomic features (from 3 categories) were significantly different (P < 0.05, |BSC| > 0.48) between the HPV+ and HPV- groups. Among the eight classifiers trained and applied for prediction of HPV status, the GLM showed the highest performance for each discriminant feature and the combined 12 features: AUC/PPV/NPV = 0.878/0.834/0.811. The GLM high prediction power was AUC/PPV/NPV = 0.849/0.731/0.788 and AUC/PPV/NPV = 0.869/0.807/0.870 for unseen test datasets for groups A and B, respectively. After eliminating the correlation among discriminant features by applying PCA analysis, the performance of the GLM was improved by 3.3%, 2.2%, and 1.8% for AUC, PPV, and NPV, respectively.

CONCLUSIONS

Results imply that GTV's for HPV+ patients exhibit higher intensities, smaller lesion size, greater sphericity/roundness, and higher spatial intensity-variation/heterogeneity. Results are suggestive that radiomic features primarily associated with the spatial arrangement and morphological appearance of the tumor on contrast-enhanced diagnostic CT datasets may be potentially used for classification of HPV status.

The impact of iterative reconstruction protocol, signal-to-background ratio and background activity on measurement of PET spatial resolution

OBJECTIVES

The present study aims to assess the impact of acquisition time, different iterative reconstruction protocols as well as image context (including contrast levels and background activities) on the measured spatial resolution in PET images.

METHODS

Discovery 690 PET/CT scanner was used to quantify spatial resolutions in terms of full width half maximum (FWHM) as derived (i) directly from capillary tubes embedded in air and (ii) indirectly from 10 mm-diameter sphere of the NEMA phantom. Different signal-to-background ratios (SBRs), background activity levels and acquisition times were applied. The emission data were reconstructed using iterative reconstruction protocols. Various combinations of iterations and subsets (it × sub) were evaluated.

RESULTS

For capillary tubes, improved FWHM values were obtained for higher it × sub, with improved performance for PSF algorithms relative to non-PSF algorithms. For the NEMA phantom, by increasing acquisition times from 1 to 5 min, intrinsic FWHM for reconstructions with it × sub 32 (54) was improved by 15.3% (13.2%), 15.1% (13.8%), 14.5% (12.8%) and 13.7% (12.7%) for OSEM, OSEM + PSF, OSEM + TOF and OSEM + PSF + TOF, respectively. Furthermore, for all reconstruction protocols, the FWHM improved with more impact for higher it × sub.

CONCLUSION

Our results indicate that PET spatial resolution is greatly affected by SBR, background activity and the choice of the reconstruction protocols.

Optimization of treatment strategy by using a machine learning model to predict survival time of patients with malignant glioma after radiotherapy

The purpose of this study was to predict the survival time of patients with malignant glioma after radiotherapy with high accuracy by considering additional clinical factors and optimize the prescription dose and treatment duration for individual patient by using a machine learning model. A total of 35 patients with malignant glioma were included in this study. The candidate features included 12 clinical features and 192 dose-volume histogram (DVH) features. The appropriate input features and parameters of the support vector machine (SVM) were selected using the genetic algorithm based on Akaike's information criterion, i.e. clinical, DVH, and both clinical and DVH features. The prediction accuracy of the SVM models was evaluated through a leave-one-out cross-validation test with residual error, which was defined as the absolute difference between the actual and predicted survival times after radiotherapy. Moreover, the influences of various values of prescription dose and treatment duration on the predicted survival time were evaluated. The prediction accuracy was significantly improved with the combined use of clinical and DVH features compared with the separate use of both features (P < 0.01, Wilcoxon signed rank test). Mean ± standard deviation of the leave-one-out cross-validation using the combined clinical and DVH features, only clinical features and only DVH features were 104.7 ± 96.5, 144.2 ± 126.1 and 204.5 ± 186.0 days, respectively. The prediction accuracy could be improved with the combination of clinical and DVH features, and our results show the potential to optimize the treatment strategy for individual patients based on a machine learning model.

11C-methionine-PET for differentiating recurrent brain tumor from radiation necrosis: radiomics approach with random forest classifier

Differentiating recurrent brain tumor from radiation necrosis is often difficult. This study aims to investigate the efficacy of 11C-methionine (MET)-PET radiomics for distinguishing recurrent brain tumor from radiation necrosis, as compared with conventional tumor-to-normal cortex (T/N) ratio evaluation. We enrolled 41 patients with metastatic brain tumor or glioma treated using radiation therapy who underwent MET-PET. The area with a standardized uptake value > 1.3 times that of the normal brain cortex was contoured. Forty-two PET features were extracted and used in a random forest classifier and the diagnostic performance was evaluated using a 10-fold cross-validation scheme. Gini index was measured to identify relevant PET parameters for classification. The reference standard was surgical histopathological analysis or more than 6 months of follow-up with MRI. Forty-four lesions were used for the analysis. Thirty-three and 11 lesions were confirmed as recurrent brain tumor and radiation necrosis, respectively. Radiomics and T/N ratio evaluation showed sensitivities of 90.1% and 60.6%, and specificities of 93.9% and 72.7% with areas under the curve of 0.98 and 0.73, respectively. Gray level co-occurrence matrix dissimilarity was the most pertinent feature for diagnosis. MET-PET radiomics yielded excellent outcome for differentiating recurrent brain tumor from radiation necrosis, which outperformed T/N ratio evaluation.

Usefulness of gradient tree boosting for predicting histological subtype and EGFR mutation status of non-small cell lung cancer on 18F FDG-PET/CT

OBJECTIVE

To develop and evaluate a radiomics approach for classifying histological subtypes and epidermal growth factor receptor (EGFR) mutation status in lung cancer on PET/CT images.

METHODS

PET/CT images of lung cancer patients were obtained from public databases and used to establish two datasets, respectively to classify histological subtypes (156 adenocarcinomas and 32 squamous cell carcinomas) and EGFR mutation status (38 mutant and 100 wild-type samples). Seven types of imaging features were obtained from PET/CT images of lung cancer. Two types of machine learning algorithms were used to predict histological subtypes and EGFR mutation status: random forest (RF) and gradient tree boosting (XGB). The classifiers used either a single type or multiple types of imaging features. In the latter case, the optimal combination of the seven types of imaging features was selected by Bayesian optimization. Receiver operating characteristic analysis, area under the curve (AUC), and tenfold cross validation were used to assess the performance of the approach.

RESULTS

In the classification of histological subtypes, the AUC values of the various classifiers were as follows: RF, single type: 0.759; XGB, single type: 0.760; RF, multiple types: 0.720; XGB, multiple types: 0.843. In the classification of EGFR mutation status, the AUC values were: RF, single type: 0.625; XGB, single type: 0.617; RF, multiple types: 0.577; XGB, multiple types: 0.659.

CONCLUSIONS

The radiomics approach to PET/CT images, together with XGB and Bayesian optimization, is useful for classifying histological subtypes and EGFR mutation status in lung cancer.

Radiomics features of the primary tumor fail to improve prediction of overall survival in large cohorts of CT- and PET-imaged head and neck cancer patients

Radiomics studies require many patients in order to power them, thus patients are often combined from different institutions and using different imaging protocols. Various studies have shown that imaging protocols affect radiomics feature values. We examined whether using data from cohorts with controlled imaging protocols improved patient outcome models. We retrospectively reviewed 726 CT and 686 PET images from head and neck cancer patients, who were divided into training or independent testing cohorts. For each patient, radiomics features with different preprocessing were calculated and two clinical variables-HPV status and tumor volume-were also included. A Cox proportional hazards model was built on the training data by using bootstrapped Lasso regression to predict overall survival. The effect of controlled imaging protocols on model performance was evaluated by subsetting the original training and independent testing cohorts to include only patients whose images were obtained using the same imaging protocol and vendor. Tumor volume, HPV status, and two radiomics covariates were selected for the CT model, resulting in an AUC of 0.72. However, volume alone produced a higher AUC, whereas adding radiomics features reduced the AUC. HPV status and one radiomics feature were selected as covariates for the PET model, resulting in an AUC of 0.59, but neither covariate was significantly associated with survival. Limiting the training and independent testing to patients with the same imaging protocol reduced the AUC for CT patients to 0.55, and no covariates were selected for PET patients. Radiomics features were not consistently associated with survival in CT or PET images of head and neck patients, even within patients with the same imaging protocol.

Effects of alterations in positron emission tomography imaging parameters on radiomics features

Radiomics studies require large patient cohorts, which often include patients imaged using different imaging protocols. We aimed to determine the impact of variability in imaging protocol parameters and interscanner variability using a phantom that produced feature values similar to those of patients. Positron emission tomography (PET) scans of a Hoffman brain phantom were acquired on GE Discovery 710, Siemens mCT, and Philips Vereos scanners. A standard-protocol scan was acquired on each machine, and then each parameter that could be changed was altered individually. The phantom was contoured with 10 regions of interest (ROIs). Values for 45 features with 2 different preprocessing techniques were extracted for each image. To determine the impact of each parameter on the reliability of each radiomics feature, the intraclass correlation coefficient (ICC) was calculated with the ROIs as the subjects and the parameter values as the raters. For interscanner comparisons, we compared the standard deviation of each radiomics feature value from the standard-protocol images to the standard deviation of the same radiomics feature from PET scans of 224 patients with non-small cell lung cancer. When the pixel size was resampled prior to feature extraction, all features had good reliability (ICC > 0.75) for the field of view and matrix size. The time per bed position had excellent reliability (ICC > 0.9) on all features. When the filter cutoff was restricted to values below 6 mm, all features had good reliability. Similarly, when subsets and iterations were restricted to reasonable values used in clinics, almost all features had good reliability. The average ratio of the standard deviation of features on the phantom scans to that of the NSCLC patient scans was 0.73 using fixed-bin-width preprocessing and 0.92 using 64-level preprocessing. Most radiomics feature values had at least good reliability when imaging protocol parameters were within clinically used ranges. However, interscanner variability was about equal to interpatient variability; therefore, caution must be used when combining patients scanned on equipment from different vendors in radiomics data sets.

Metabolic Biomarker-Based BRAFV600 Mutation Association and Prediction in Melanoma

The aim of this study was to associate and predict B-rapidly accelerated fibrosarcoma valine 600 (BRAFV600) mutation status with both conventional and radiomics ¹⁸F-FDG PET/CT features, while exploring several methods of feature selection in melanoma radiomics. Methods: Seventy unresectable stage III-IV melanoma patients who underwent a baseline ¹⁸F-FDG PET/CT scan were identified. Patients were assigned to the BRAFV600 group or BRAF wild-type group according to mutational status. ¹⁸F-FDG uptake quantification was performed by semiautomatic lesion delineation. Four hundred eighty radiomics features and 4 conventional PET features (SUV_max, SUV_mean, SUV_peak, and total lesion glycolysis) were extracted per lesion. Six different methods of feature selection were implemented, and 10-fold cross-validated predictive models were built for each. Model performances were evaluated with areas under the curve (AUCs) for the receiver operating characteristic curves. Results: Thirty-five BRAFV600 mutated patients (100 lesions) and 35 BRAF wild-type patients (79 lesions) were analyzed. AUCs predicting the BRAFV600 mutation varied from 0.54 to 0.62 and were susceptible to feature selection method. The best AUCs were achieved by feature selection based on literature, a penalized binary logistic regression model, and random forest model. No significant difference was found between the BRAFV600 and BRAF wild-type group in conventional PET features or predictive value. Conclusion: BRAFV600 mutation status is not associated with, nor can it be predicted with, conventional PET features, whereas radiomics features were of low predictive value (AUC = 0.62). We showed feature selection methods to influence predictive model performance, describing and evaluating 6 unique methods. Detecting BRAFV600 status in melanoma based on ¹⁸F-FDG PET/CT alone does not yet provide clinically relevant knowledge.

The stability of imaging biomarkers in radiomics: a framework for evaluation

This paper studies the sensitivity of a range of image texture parameters used in radiomics to: (i) the number of intensity levels, (ii) the method of quantisation to select the intensity levels and (iii) the use of an intensity threshold. 43 commonly used texture features were studied for the gross target volume outlined on the CT component of PET/CT scans of 50 patients with non-small cell lung carcinoma (NSCLC). All cases were quantised for all values between 4 and 128 intensity levels using four commonly used quantisation methods. All results were analysed with and without a threshold range of  -200 HU to 300 HU. Cases were ranked for each texture feature and for all quantisation methods with the Spearman's rank correlation coefficient determined to evaluate stability. Results showed large fluctuations in ranking, particularly for low numbers of levels, differences between quantisation methods and with the use of a threshold, with values Spearman's Rank Correlation for many parameters below 0.2. Our results demonstrated the sensitivity of radiomics features to the parameters used during analysis and highlight the risk of low reproducibility comparing studies with slightly different parameters. In terms of the lung cancer CT datasets, this study supports the use of 128 intensity levels, the same uniform quantiser applied to all scans and thresholding of the data. It also supports several of the features recommended in the literature for such studies such as skewness and kurtosis. A recommended framework is presented for curation of the data analysis process to ensure stability of results.

Experimental Multicenter and Multivendor Evaluation of the Performance of PET Radiomic Features Using 3-Dimensionally Printed Phantom Inserts

The sensitivity of radiomic features to several confounding factors, such as reconstruction settings, makes clinical use challenging. To investigate the impact of harmonized image reconstructions on feature consistency, a multicenter phantom study was performed using 3-dimensionally printed phantom inserts reflecting realistic tumor shapes and heterogeneity uptakes. Methods: Tumors extracted from real PET/CT scans of patients with non-small cell lung cancer served as model for three 3-dimensionally printed inserts. Different heterogeneity pattern were realized by printing separate compartments that could be filled with different activity solutions. The inserts were placed in the National Electrical Manufacturers Association image-quality phantom and scanned various times. First, a list-mode scan was acquired and 5 statistically equal replicates were reconstructed. Second, the phantom was scanned 4 times on the same scanner. Third, the phantom was scanned on 6 PET/CT systems. All images were reconstructed using EANM Research Ltd. (EARL)-compliant and locally clinically preferred reconstructions. EARL-compliant reconstructions were performed without (EARL1) or with (EARL2) point-spread function. Images were analyzed with and without resampling to 2-mm cubic voxels. Images were discretized with a fixed bin width (FBW) of 0.25 and a fixed bin number (FBN) of 64. The intraclass correlation coefficient (ICC) of each scan setup was calculated and compared across reconstruction settings. An ICC above 0.75 was regarded as high. Results: The percentage of features yielding a high ICC was largest for the statistically equal replicates (70%-91% for FBN; 90%-96% for FBW discretization). For scans acquired on the same system, the percentage decreased, but most features still resulted in a high ICC (FBN, 52%-63%; FBW, 75%-85%). The percentage of features yielding a high ICC decreased more in the multicenter setting. In this case, the percentage of features yielding a high ICC was larger for images reconstructed with EARL-compliant reconstructions: for example, 40% for EARL1 and 60% for EARL2 versus 21% for the clinically preferred setting for FBW discretization. When discretized with FBW and resampled to isotropic voxels, this benefit was more pronounced. Conclusion: EARL-compliant reconstructions harmonize a wide range of radiomic features. FBW discretization and a sampling to isotropic voxels enhances the benefits of EARL-compliant reconstructions.

Radiomics in nuclear medicine: robustness, reproducibility, standardization, and how to avoid data analysis traps and replication crisis

Radiomics in nuclear medicine is rapidly expanding. Reproducibility of radiomics studies in multicentre settings is an important criterion for clinical translation. We therefore performed a meta-analysis to investigate reproducibility of radiomics biomarkers in PET imaging and to obtain quantitative information regarding their sensitivity to variations in various imaging and radiomics-related factors as well as their inherent sensitivity. Additionally, we identify and describe data analysis pitfalls that affect the reproducibility and generalizability of radiomics studies. After a systematic literature search, 42 studies were included in the qualitative synthesis, and data from 21 were used for the quantitative meta-analysis. Data concerning measurement agreement and reliability were collected for 21 of 38 different factors associated with image acquisition, reconstruction, segmentation and radiomics-specific processing steps. Variations in voxel size, segmentation and several reconstruction parameters strongly affected reproducibility, but the level of evidence remained weak. Based on the meta-analysis, we also assessed inherent sensitivity to variations of 110 PET image biomarkers. SUV_mean and SUV_max were found to be reliable, whereas image biomarkers based on the neighbourhood grey tone difference matrix and most biomarkers based on the size zone matrix were found to be highly sensitive to variations, and should be used with care in multicentre settings. Lastly, we identify 11 data analysis pitfalls. These pitfalls concern model validation and information leakage during model development, but also relate to reporting and the software used for data analysis. Avoiding such pitfalls is essential for minimizing bias in the results and to enable reproduction and validation of radiomics studies.

Impact of intensity discretization on textural indices of [18F]FDG-PET tumour heterogeneity in lung cancer patients

Quantifying tumour heterogeneity from [¹⁸F]FDG-PET images promises benefits for treatment selection of cancer patients. Here, the calculation of texture parameters mandates an initial discretization step (binning) to reduce the number of intensity levels. Typically, three types of discrimination methods are used: lesion relative resampling (LRR) with fixed bin number, lesion absolute resampling (LAR) and absolute resampling (AR) with fixed bin widths. We investigated the effects of varying bin widths or bin number using 27 commonly cited local and regional texture indices (TIs) applied on lung tumour volumes. The data set were extracted from 58 lung cancer patients, with three different and robust tumour segmentation methods. In our cohort, the variations of the mean value as the function of the bin widths were similar for TIs calculated with LAR and AR quantification. The TI histograms calculated by LRR method showed distinct behaviour and its numerical values substantially effected by the selected bin number. The correlations of the AR and LAR based TIs demonstrated no principal differences between these methods. However, no correlation was found for the interrelationship between the TIs calculated by LRR and LAR (or AR) discretization method. Visual classification of the texture was also performed for each lesion. This classification analysis revealed that the parameters show statistically significant correlation with the visual score, if LAR or AR discretization method is considered, in contrast to LRR. Moreover, all the resulted tendencies were similar regardless the segmentation methods and the type of textural features involved in this work.

Neural network training for cross-protocol radiomic feature standardization in computed tomography

Radiomics has shown promising results in several medical studies, yet it suffers from a limited discrimination and informative capability as well as a high variation and correlation with the tomographic scanner types, pixel spacing, acquisition protocol, and reconstruction parameters. We propose and compare two methods to transform quantitative image features in order to improve their stability across varying image acquisition parameters while preserving the texture discrimination abilities. In this way, variations in extracted features are representative of true physiopathological tissue changes in the scanned patients. A first approach is based on a two-layer neural network that can learn a nonlinear standardization transformation of various types of features including handcrafted and deep features. Second, domain adversarial training is explored to increase the invariance of the transformed features to the scanner of origin. The generalization of the proposed approach to unseen textures and unseen scanners is demonstrated by a set of experiments using a publicly available computed tomography texture phantom dataset scanned with various imaging devices and parameters.

Challenges and caveats of a multi-center retrospective radiomics study: an example of early treatment response assessment for NSCLC patients using FDG-PET/CT radiomics

Background

Prognostic models based on individual patient characteristics can improve treatment decisions and outcome in the future. In many (radiomic) studies, small size and heterogeneity of datasets is a challenge that often limits performance and potential clinical applicability of these models. The current study is example of a retrospective multi-centric study with challenges and caveats. To highlight common issues and emphasize potential pitfalls, we aimed for an extensive analysis of these multi-center pre-treatment datasets, with an additional ¹⁸F-fluorodeoxyglucose (FDG) positron emission tomography/computed tomography (PET/CT) scan acquired during treatment.

Methods

The dataset consisted of 138 stage II-IV non-small cell lung cancer (NSCLC) patients from four different cohorts acquired from three different institutes. The differences between the cohorts were compared in terms of clinical characteristics and using the so-called ‘cohort differences model’ approach. Moreover, the potential prognostic performances for overall survival of radiomic features extracted from CT or FDG-PET, or relative or absolute differences between the scans at the two time points, were assessed using the LASSO regression method. Furthermore, the performances of five different classifiers were evaluated for all image sets.

Results

The individual cohorts substantially differed in terms of patient characteristics. Moreover, the cohort differences model indicated statistically significant differences between the cohorts. Neither LASSO nor any of the tested classifiers resulted in a clinical relevant prognostic model that could be validated on the available datasets.

Conclusion

The results imply that the study might have been influenced by a limited sample size, heterogeneous patient characteristics, and inconsistent imaging parameters. No prognostic performance of FDG-PET or CT based radiomics models can be reported. This study highlights the necessity of extensive evaluations of cohorts and of validation datasets, especially in retrospective multi-centric datasets.

Neo-adjuvant chemoradiotherapy response prediction using MRI based ensemble learning method in rectal cancer patients

OBJECTIVES

The aim of this study was to investigate and validate the performance of individual and ensemble machine learning models (EMLMs) based on magnetic resonance imaging (MRI) to predict neo-adjuvant chemoradiation therapy (nCRT) response in rectal cancer patients. We also aimed to study the effect of Laplacian of Gaussian (LOG) filter on EMLMs predictive performance.

METHODS

98 rectal cancer patients were divided into a training (n = 53) and a validation set (n = 45). All patients underwent MRI a week before nCRT. Several features from intensity, shape and texture feature sets were extracted from MR images. SVM, Bayesian network, neural network and KNN classifiers were used individually and together for response prediction. Predictive performance was evaluated using the area under the receiver operator characteristic (ROC) curve (AUC).

RESULTS

Patients' nCRT responses included 17 patients with Grade 0, 28 with Grade 1, 34 with Grade 2, and 19 with Grade 3 according to AJCC/CAP pathologic grading. In without preprocessing MR Image the best result was for Bayesian network classifier with AUC and accuracy of 75.2% and 80.9% respectively, which was confirmed in the validation set with an AUC and accuracy of 74% and 79% respectively. In EMLMs the best result was for 4 (SVM.NN.BN.KNN) classifier EMLM with AUC and accuracy of 97.8% and 92.8% in testing and 95% and 90% in validation set respectively.

CONCLUSIONS

In conclusion, we observed that machine learning methods can used to predict nCRT response in patients with rectal cancer. Preprocessing LOG filters and EL models can improve the prediction process.

Radiomic feature stability across 4D respiratory phases and its impact on lung tumor prognosis prediction

Radiomic analysis has recently demonstrated versatile uses in improving diagnostic and prognostic prediction accuracy for lung cancer. However, since lung tumors are subject to substantial motion due to respiration, the stability of radiomic features over the respiratory cycle of the patient needs to be investigated to better evaluate the robustness of the inter-patient feature variability for clinical applications, and its impact in such applications needs to be assessed. A full panel of 841 radiomic features, including tumor intensity, shape, texture, and wavelet features, were extracted from individual phases of a four-dimensional (4D) computed tomography on 20 early-stage non-small-cell lung cancer (NSCLC) patients. The stability of each radiomic feature was assessed across different phase images of the same patient using the coefficient of variation (COV). The relationship between individual COVs and tumor motion magnitude was inspected. Population COVs, the mean COVs of all 20 patients, were used to evaluate feature motion stability and categorize the radiomic features into 4 different groups. The two extremes, the Very Small group (COV≤5%) and the Large group (COV>20%), each accounted for about a quarter of the features. Shape features were the most stable, with COV≤10% for all features. A clinical study was subsequently conducted using 140 early-stage NSCLC patients. Radiomic features were employed to predict the overall survival with a 500-round bootstrapping. Identical multiple regression model development process was applied, and the model performance was compared between models with and without a feature pre-selection step based on 4D COV to pre-exclude unstable features. Among the systematically tested cutoff values, feature pre-selection with 4D COV≤5% achieved the optimal model performance. The resulting 3-feature radiomic model significantly outperformed its counterpart with no 4D COV pre-selection, with P = 2.16x10^-27 in the one-tailed t-test comparing the prediction performances of the two models.

Radiomic Feature Robustness and Reproducibility in Quantitative Bone Radiography: A Study on Radiologic Parameter Changes

The purpose of this study was to investigate the robustness of different radiography radiomic features over different radiologic parameters including kV, mAs, filtration, tube angles, and source skin distance (SSD). A tibia bone phantom was prepared and all imaging studies was conducted on this phantom. Different radiologic parameters including kV, mAs, filtration, tube angles, and SSD were studied. A region of interest was drawn on the images and many features from different feature sets including histogram, gradient, run-length matrix, co-occurrence matrix, autoregressive model, and wavelet derived parameters were extracted. All radiomic features were categorized based on coefficient of variation (COV). Bland-Altman analysis also was used to evaluate the mean, standard deviation, and upper/lower reproducibility limits for radiomic features in response to variation in each testing parameters. Results on COV in all features showed that 22%, 34%, and 45% of features were most robust (COV ≤ 5%) against kV, mAs, and SSD respectively and there was no robust features against filtration and tube angle. Also, all features (100%) and 76% of which showed large variations (COV > 20%) against filtrations and tube angle respectively. Autoregressive model feature set has no robust features against all radiologic parameters. Features including sum-average, sum-entropy, correlation, mean, and percentile (50, 90, and 99) belong to co-occurrence matrix and histogram feature sets were found as most robust features. Bland-Altman analysis showed the high reproducibity of some feature sets against radiologic parameter changes. The results presented here indicated that radiologic parameters have great impacts on radiomic feature values and caution should be taken into account when work with these features. In quantitative bone studies, robust features with low COV can be selected for clinical or research applications. Reproducible features also can be obtained using Bland-Altman analysis.

Radiomics: Principles and radiotherapy applications

Radiomics is defined as the extraction of a large quantity of quantitative image features. The different radiomic indexes that have been proposed in the literature are described as well as the various factors that have an impact on the robustness of these indexes. We will see that several hundred quantitative features can be extracted per lesion and imaging modality. The ever-growing number of features studied raises the question of the statistical method of analysis used. This review addresses the research supporting the clinical use of radiomics in oncology in the staging of disease, discrimination between healthy and pathological tissues, the identification of genetic features, the prediction of patient survival, the response to treatment, the recurrence after radiotherapy and chemoradiotherapy and the side effects. Based on the existing literature, it remains difficult to identify features that should be used for current clinical practice.

Impact of acquisition count statistics reduction and SUV discretization on PET radiomic features in pediatric 18F-FDG-PET/MRI examinations

PURPOSE

The evaluation of features robustness with respect to acquisition and post-processing parameter changes is fundamental for the reliability of radiomics studies. The aim of this study was to investigate the sensitivity of PET radiomic features to acquisition statistics reduction and standardized-uptake-volume (SUV) discretization in PET/MRI pediatric examinations.

METHODS

Twenty-seven lesions were detected from the analysis of twenty-one 18F-FDG-PET/MRI pediatric examinations. By decreasing the count-statistics of the original list-mode data (3 MBq/kg), injected activity reduction was simulated. Two SUV discretization approaches were applied: 1) resampling lesion SUV range into fixed bins numbers (FBN); 2) rounding lesion SUV into fixed bin size (FBS). One hundred and six radiomic features were extracted. Intraclass Correlation Coefficient (ICC), Spearman correlation coefficient and coefficient-of-variation (COV) were calculated to assess feature reproducibility between low tracer activities and full tracer activity feature values.

RESULTS

More than 70% of Shape and first order features, and around 70% and 40% of textural features, when using FBS and FBN methods respectively, resulted robust till 1.2 MBk/kg. Differences in median features reproducibility (ICC) between FBS and FBN datasets were statistically significant for every activity level independently from bin number/size, with higher values for FBS. Differences in median Spearman coefficient (i.e. patient ranking according to feature values) were not statistically significant, varying the intensity resolution (i.e. bin number/size) for either FBS and FBN methods.

CONCLUSIONS

For each simulated count-statistic level, robust PET radiomic features were determined for pediatric PET/MRI examinations. A larger number of robust features were detected when using FBS methods.

The Applications of Radiomics in Precision Diagnosis and Treatment of Oncology: Opportunities and Challenges

Medical imaging can assess the tumor and its environment in their entirety, which makes it suitable for monitoring the temporal and spatial characteristics of the tumor. Progress in computational methods, especially in artificial intelligence for medical image process and analysis, has converted these images into quantitative and minable data associated with clinical events in oncology management. This concept was first described as radiomics in 2012. Since then, computer scientists, radiologists, and oncologists have gravitated towards this new tool and exploited advanced methodologies to mine the information behind medical images. On the basis of a great quantity of radiographic images and novel computational technologies, researchers developed and validated radiomic models that may improve the accuracy of diagnoses and therapy response assessments. Here, we review the recent methodological developments in radiomics, including data acquisition, tumor segmentation, feature extraction, and modelling, as well as the rapidly developing deep learning technology. Moreover, we outline the main applications of radiomics in diagnosis, treatment planning and evaluations in the field of oncology with the aim of developing quantitative and personalized medicine. Finally, we discuss the challenges in the field of radiomics and the scope and clinical applicability of these methods.

Repeatability of 18 F-FDG PET radiomic features: A phantom study to explore sensitivity to image reconstruction settings, noise, and delineation method

BACKGROUND

18 F-fluoro-2-deoxy-D-Glucose positron emission tomography (18 F-FDG PET) radiomics has the potential to guide the clinical decision making in cancer patients, but validation is required before radiomics can be implemented in the clinical setting. The aim of this study was to explore how feature space reduction and repeatability of 18 F-FDG PET radiomic features are affected by various sources of variation such as underlying data (e.g., object size and uptake), image reconstruction methods and settings, noise, discretization method, and delineation method.

METHODS

The NEMA image quality phantom was scanned with various sphere-to-background ratios (SBR), simulating different activity uptakes, including spheres with low uptake, that is, SBR smaller than 1. Furthermore, images of a phantom containing 3D printed inserts reflecting realistic heterogeneity uptake patterns were acquired. Data were reconstructed using various matrix sizes, reconstruction algorithms, and scan durations (noise). For every specific reconstruction and noise level, ten statistically equal replicates were generated. The phantom inserts were delineated using CT and PET-based segmentation methods. A total of 246 radiomic features was extracted from each image dataset. Images were discretized with a fixed number of 64 bins (FBN) and a fixed bin width (FBW) of 0.25 for the high and a FBW of 0.05 for the low uptake data. In terms of feature reduction, we determined the impact of these factors on the composition of feature clusters, which were defined on the basis of Spearman's correlation matrices. To assess feature repeatability, the intraclass correlation coefficient was calculated over the ten replicates.

RESULTS

In general, larger spheres with high uptake resulted in better repeatability compared to smaller low uptake spheres. In terms of repeatability, features extracted from heterogeneous phantom inserts were comparable to features extracted from bigger high uptake spheres. For example, for an EARL-compliant reconstruction, larger and smaller high uptake spheres yielded good repeatability for 32% and 30% of the features, while the heterogeneous inserts resulted in 34% repeatable features. For the low uptake spheres, this was the case for 22% and 20% of the features for bigger and smaller spheres, respectively. Images reconstructed with point-spread-function (PSF) resulted in the highest repeatability when compared with OSEM or time-of-flight, for example, 53%, 30%, and 32% of repeatable features, respectively (for unsmoothed data, discretized with FBN, 300 s scan duration). Reducing image noise (increasing scan duration and smoothing) and using CT-based segmentation for the low uptake spheres yielded improved repeatability. FBW discretization resulted in higher repeatability than FBN discretization, for example, 89% and 35% of the features, respectively (for the EARL-compliant reconstruction and larger high uptake spheres).

CONCLUSION

Feature space reduction and repeatability of 18 F-FDG PET radiomic features depended on all studied factors. The high sensitivity of PET radiomic features to image quality suggests that a high level of image acquisition and preprocessing standardization is required to be used as clinical imaging biomarker.

Radiomics in Oncological PET/CT: a Methodological Overview

Radiomics is a medical imaging analysis approach based on computer-vision. Metabolic radiomics in particular analyses the spatial distribution patterns of molecular metabolism on PET images. Measuring intratumoral heterogeneity via image is one of the main targets of radiomics research, and it aims to build a image-based model for better patient management. The workflow of radiomics using texture analysis follows these steps: 1) imaging (image acquisition and reconstruction); 2) preprocessing (segmentation & quantization); 3) quantification (texture matrix design & texture feature extraction); and 4) analysis (statistics and/or machine learning). The parameters or conditions at each of these steps are effect on the results. In statistical testing or modeling, problems such as multiple comparisons, dependence on other variables, and high dimensionality of small sample size data should be considered. Standardization of methodology and harmonization of image quality are one of the most important challenges with radiomics methodology. Even though there are current issues in radiomics methodology, it is expected that radiomics will be clinically useful in personalized medicine for oncology.

Texture Analysis of Imaging: What Radiologists Need to Know

OBJECTIVE

Radiologic texture is the variation in image intensities within an image and is an important part of radiomics. The objective of this article is to discuss some parameters that affect the performance of texture metrics and propose recommendations that can guide both the design and evaluation of future radiomics studies.

CONCLUSION

A variety of texture-extraction techniques are used to assess clinical imaging data. Currently, no consensus exists regarding workflow, including acquisition, extraction, or reporting of variable settings leading to poor reproducibility.

Challenges and Promises of PET Radiomics

PURPOSE

Radiomics describes the extraction of multiple, otherwise invisible, features from medical images that, with bioinformatic approaches, can be used to provide additional information that can predict underlying tumor biology and behavior.

METHODS AND MATERIALS

Radiomic signatures can be used alone or with other patient-specific data to improve tumor phenotyping, treatment response prediction, and prognosis, noninvasively. The data describing 18F-fluorodeoxyglucose positron emission tomography radiomics, often using texture or heterogeneity parameters, are increasing rapidly.

RESULTS

In relation to radiation therapy practice, early data have reported the use of radiomic approaches to better define tumor volumes and predict radiation toxicity and treatment response.

CONCLUSIONS

Although at an early stage of development, with many technical challenges remaining and a need for standardization, promise nevertheless exists that PET radiomics will contribute to personalized medicine, especially with the availability of increased computing power and the development of machine-learning approaches for imaging.

Variability and Repeatability of Quantitative Uptake Metrics in 18F-FDG PET/CT of Non-Small Cell Lung Cancer: Impact of Segmentation Method, Uptake Interval, and Reconstruction Protocol

There is increased interest in various new quantitative uptake metrics beyond SUV in oncologic PET/CT studies. The purpose of this study was to investigate the variability and test-retest ratio (TRT) of metabolically active tumor volume (MATV) measurements and several other new quantitative metrics in non-small cell lung cancer using ¹⁸F-FDG PET/CT with different segmentation methods, user interactions, uptake intervals, and reconstruction protocols. Methods: Ten patients with advanced non-small cell lung cancer received 2 series of 2 whole-body ¹⁸F-FDG PET/CT scans at 60 min after injection and at 90 min after injection. PET data were reconstructed with 4 different protocols. Eight segmentation methods were applied to delineate lesions with and without a tumor mask. MATV, SUV_max, SUV_mean, total lesion glycolysis, and intralesional heterogeneity features were derived. Variability and repeatability were evaluated using a generalized-estimating-equation statistical model with Bonferroni adjustment for multiple comparisons. The statistical model, including interaction between uptake interval and reconstruction protocol, was applied individually to the data obtained from each segmentation method. Results: Without masking, none of the segmentation methods could delineate all lesions correctly. MATV was affected by both uptake interval and reconstruction settings for most segmentation methods. Similar observations were obtained for the uptake metrics SUV_max, SUV_mean, total lesion glycolysis, homogeneity, entropy, and zone percentage. No effect of uptake interval was observed on TRT metrics, whereas the reconstruction protocol affected the TRT of SUV_max Overall, segmentation methods showing poor quantitative performance in one condition showed better performance in other (combined) conditions. For some metrics, a clear statistical interaction was found between the segmentation method and both uptake interval and reconstruction protocol. Conclusion: All segmentation results need to be reviewed critically. MATV and other quantitative uptake metrics, as well as their TRT, depend on segmentation method, uptake interval, and reconstruction protocol. To obtain quantitative reliable metrics, with good TRT performance, the optimal segmentation method depends on local imaging procedure, the PET/CT system, or reconstruction protocol. Rigid harmonization of imaging procedure and PET/CT performance will be helpful in mitigating this variability.

Parameters Influencing PET Imaging Features: A Phantom Study with Irregular and Heterogeneous Synthetic Lesions

Aim

To evaluate reproducibility and stability of radiomic features as effects of the use of different volume segmentation methods and reconstruction settings. The potential of radiomics in really capturing the presence of heterogeneous tumor uptake and irregular shape was also investigated.

Materials and Methods

An anthropomorphic phantom miming real clinical situations including synthetic lesions with irregular shape and nonuniform radiotracer uptake was used. ¹⁸F-FDG PET/CT measurements of the phantom were performed including 38 lesions of different shape, size, lesion-to-background ratio, and radiotracer uptake distribution. Different reconstruction parameters and segmentation methods were considered. COVs were calculated to quantify feature variations over the different reconstruction settings. Friedman test was applied to the values of the radiomic features obtained for the considered segmentation approaches. Two sets of test-retest measurement were acquired and the pairwise intraclass correlation coefficient was calculated. Fifty-eight morphological and statistical features were extracted from the segmented lesion volumes. A Mann–Whitney test was used to evaluate significant differences among each feature when calculated from heterogeneous versus homogeneous uptake. The significance of each radiomic feature in terms of capturing heterogeneity was evaluated also by testing correlation with gold standard indexes of heterogeneity and sphericity.

Results

The choice of the segmentation method has a strong impact on the stability of radiomic features (less than 20% can be considered stable features). Reconstruction affects the estimate of radiomic features (only 26% are stable). Thirty-one radiomic features (53%) resulted to be reproducible, 11 of them are able to discriminate heterogeneity. Among these, we found a subset of 3 radiomic features strongly correlated with GS heterogeneity index that can be suggested as good features for retrospective evaluations.

The developing role of FDG PET imaging for prognostication and radiotherapy target volume delineation in non-small cell lung cancer

Advancements in functional imaging technology have allowed new possibilities in contouring of target volumes, monitoring therapy, and predicting treatment outcome in non-small cell lung cancer (NSCLC). Consequently, the role of ¹⁸F-fluorodeoxyglucose positron emission tomography (FDG PET) has expanded in the last decades from a stand-alone diagnostic tool to a versatile instrument integrated with computed tomography (CT), with a prominent role in lung cancer radiotherapy. This review outlines the most recent literature on developments in FDG PET imaging for prognostication and radiotherapy target volume delineation (TVD) in NSCLC. We also describe the challenges facing the clinical implementation of these developments and present new ideas for future research.

Gold nanoparticle-induced sonosensitization enhances the antitumor activity of ultrasound in colon tumor-bearing mice

PURPOSE

As a noninvasive and nonionizing radiation, ultrasound can be focused remotely, transferring acoustic energy deep in the body, thereby addressing the penetration depth barrier of the light-based therapies. In cancer therapy, the effectiveness of ultrasound can be enhanced by utilizing nanomaterials that exhibit sonosensitizing properties called as nanosonosensitizers. The gold nanoparticle (AuNP) has been recently presented as a potent nanosonosensitizer with the potential to simultaneously enhance both the thermal and mechanical interactions of ultrasound with the tissue of the human body. Accordingly, this paper attempts to evaluate the in vivo antitumor efficiency of ultrasound in combination with AuNP.

METHODS

BALB/c mice-bearing CT26 colorectal tumor model was intraperitoneally injected with AuNPs and then subjected to ultrasound irradiation (1 MHz; 2 W/cm2 ; 10 min) for three sessions. Furthermore, [18 F]FDG (2-deoxy-2-[18 F]fluoro-d-glucose) positron-emission tomography (PET) imaging was performed and the radiomic features from different feature categorizes were extracted to quantify the tumors' phenotype.

RESULTS

The tumors were dramatically shrunk and the mice appeared healthy over 21 days of study span without the evidence of relapse. The animals treated with AuNP + ultrasound exhibited an obvious decline in tumor metabolic parameters such as standard uptake value (SUV), total lesion glycolysis (TLG), and metabolic tumor volume (MTV) compared to other treatment groups.

CONCLUSION

These findings support the use of AuNP as a potent sonosensitizing agent with the potential to use the thermal and mechanical effects of ultrasound so as to cause damage to the focused tumor site, resulting in an improved antitumor efficacy.

LIFEx: A Freeware for Radiomic Feature Calculation in Multimodality Imaging to Accelerate Advances in the Characterization of Tumor Heterogeneity

Textural and shape analysis is gaining considerable interest in medical imaging, particularly to identify parameters characterizing tumor heterogeneity and to feed radiomic models. Here, we present a free, multiplatform, and easy-to-use freeware called LIFEx, which enables the calculation of conventional, histogram-based, textural, and shape features from PET, SPECT, MR, CT, and US images, or from any combination of imaging modalities. The application does not require any programming skills and was developed for medical imaging professionals. The goal is that independent and multicenter evidence of the usefulness and limitations of radiomic features for characterization of tumor heterogeneity and subsequent patient management can be gathered. Many options are offered for interactive textural index calculation and for increasing the reproducibility among centers. The software already benefits from a large user community (more than 800 registered users), and interactions within that community are part of the development strategy.Significance: This study presents a user-friendly, multi-platform freeware to extract radiomic features from PET, SPECT, MR, CT, and US images, or any combination of imaging modalities. Cancer Res; 78(16); 4786-9. ©2018 AACR.

Phantom, clinical, and texture indices evaluation and optimization of a penalized-likelihood image reconstruction method (Q.Clear) on a BGO PET/CT scanner

INTRODUCTION

The aim of this study was to evaluate the behavior of a penalized-likelihood image reconstruction method (Q.Clear) under different count statistics and lesion-to-background ratios (LBR) on a BGO scanner, in order to obtain an optimum penalization factor (β value) to study and optimize for different acquisition protocols and clinical goals.

METHODS

Both phantom and patient images were evaluated. Data from an image quality phantom were acquired using different Lesion-to-Background ratios and acquisition times. Then, each series of the phantom was reconstructed using β values between 50 and 500, at intervals of 50. Hot and cold contrasts were obtained, as well as background variability and contrast-to-noise ratio (CNR). Fifteen ¹⁸ F-FDG patients (five brain scans and 10 torso acquisitions) were acquired and reconstructed using the same β values as in the phantom reconstructions. From each lesion in the torso acquisition, noise, contrast, and signal-to-noise ratio (SNR) were computed. Image quality was assessed by two different nuclear medicine physicians. Additionally, the behaviors of 12 different textural indices were studied over 20 different lesions.

RESULTS

Q.Clear quantification and optimization in patient studies depends on the activity concentration as well as on the lesion size. In the studied range, an increase on β is translated in a decrease in lesion contrast and noise. The net product is an overall increase in the SNR, presenting a tendency to a steady value similar to the CNR in phantom data. As the activity concentration or the sphere size increase the optimal β increases, similar results are obtained from clinical data. From the subjective quality assessment, the optimal β value for torso scans is in a range between 300 and 400, and from 100 to 200 for brain scans. For the recommended torso β values, texture indices present coefficients of variation below 10%.

CONCLUSIONS

Our phantom and patients demonstrate that improvement of CNR and SNR of Q.Clear algorithm which depends on the studied conditions and the penalization factor. Using the Q.Clear reconstruction algorithm in a BGO scanner, a β value of 350 and 200 appears to be the optimal value for 18F-FDG oncology and brain PET/CT, respectively.

Repeatability and Reproducibility of Radiomic Features: A Systematic Review

Purpose:

An ever-growing number of predictive models used to inform clinical decision making have included quantitative, computer-extracted imaging biomarkers, or “radiomic features.” Broadly generalizable validity of radiomics-assisted models may be impeded by concerns about reproducibility. We offer a qualitative synthesis of 41 studies that specifically investigated the repeatability and reproducibility of radiomic features, derived from a systematic review of published peer-reviewed literature.

Methods and Materials:

The PubMed electronic database was searched using combinations of the broad Haynes and Ingui filters along with a set of text words specific to cancer, radiomics (including texture analyses), reproducibility, and repeatability. This review has been reported in compliance with Preferred Reporting Items for Systematic Reviews and Meta-Analyses guidelines. From each full-text article, information was extracted regarding cancer type, class of radiomic feature examined, reporting quality of key processing steps, and statistical metric used to segregate stable features.

Results:

Among 624 unique records, 41 full-text articles were subjected to review. The studies primarily addressed non-small cell lung cancer and oropharyngeal cancer. Only 7 studies addressed in detail every methodologic aspect related to image acquisition, preprocessing, and feature extraction. The repeatability and reproducibility of radiomic features are sensitive at various degrees to processing details such as image acquisition settings, image reconstruction algorithm, digital image preprocessing, and software used to extract radiomic features. First-order features were overall more reproducible than shape metrics and textural features. Entropy was consistently reported as one of the most stable first-order features. There was no emergent consensus regarding either shape metrics or textural features; however, coarseness and contrast appeared among the least reproducible.

Conclusions:

Investigations of feature repeatability and reproducibility are currently limited to a small number of cancer types. Reporting quality could be improved regarding details of feature extraction software, digital image manipulation (preprocessing), and the cutoff value used to distinguish stable features.

We offer a qualitative synthesis of 41 studies that specifically investigated the repeatability and reproducibility of radiomic features. The repeatability and reproducibility of radiomic features are sensitive at various degrees to image quality and to software used to extract radiomic features. Investigations of feature repeatability and reproducibility are currently limited to a small number of cancer types. No consensus was found regarding the most repeatable and reproducible features with respect to different settings.

Radiomics in Oncological PET/CT: Clinical Applications

18F-fluorodeoxyglucose (FDG) positron emission tomography/computed tomography (PET/CT) is widely used for staging, evaluating treatment response, and predicting prognosis in malignant diseases. FDG uptake and volumetric PET parameters such as metabolic tumor volume have been used and are still used as conventional PET parameters to assess biological characteristics of tumors. However, in recent years, additional features derived from PET images by computational processing have been found to reflect intratumoral heterogeneity, which is related to biological tumor features, and to provide additional predictive and prognostic information, which leads to the concept of radiomics. In this review, we focus on recent clinical studies of malignant diseases that investigated intratumoral heterogeneity on PET/CT, and we discuss its clinical role in various cancers.

Positron emission tomography-based assessment of metabolic gradient and other prognostic features in sarcoma

Intratumoral heterogeneity biomarkers derived from positron emission tomography (PET) imaging with fluorodeoxyglucose (FDG) are of interest for a number of cancers, including sarcoma. A range of radiomic texture variables, adapted from general methodologies for image analysis, has shown promise in the setting. In the context of sarcoma, our group introduced an alternative model-based approach to the measurement of heterogeneity. In this approach, the heterogeneity of a tumor is characterized by the extent to which the 3-D FDG uptake pattern deviates from a simple elliptically contoured structure. By using a nonparametric analysis of the uptake profile obtained from this spatial model, a variable assessing the metabolic gradient of the tumor is developed. The work explores the prognostic potential of this new variable in the context of FDG-PET imaging of sarcoma. A mature clinical series involving 197 patients, 88 of whom have complete time-to-death information, is used. Texture variables based on the imaging data are also evaluated in this series and a range of appropriate machine learning methodologies are then used to explore the complementary prognostic roles for structure and texture variables. We conclude that both texture-based and model-based variables can be combined to achieve enhanced prognostic assessments of outcome for patients with sarcoma based on FDG-PET imaging information.

Radiomics in Nuclear Medicine Applied to Radiation Therapy: Methods, Pitfalls, and Challenges

Radiomics is a recent area of research in precision medicine and is based on the extraction of a large variety of features from medical images. In the field of radiation oncology, comprehensive image analysis is crucial to personalization of treatments. A better characterization of local heterogeneity and the shape of the tumor, depicting individual cancer aggressiveness, could guide dose planning and suggest volumes in which a higher dose is needed for better tumor control. In addition, noninvasive imaging features that could predict treatment outcome from baseline scans could help the radiation oncologist to determine the best treatment strategies and to stratify patients as at low risk or high risk of recurrence. Nuclear medicine molecular imaging reflects information regarding biological processes in the tumor thanks to a wide range of radiotracers. Many studies involving ¹⁸F-fluorodeoxyglucose positron emission tomography suggest an added value of radiomics compared with the use of conventional PET metrics such as standardized uptake value for both tumor diagnosis and prediction of recurrence or treatment outcome. However, these promising results should not hide technical difficulties that still currently prevent the approach from being widely studied or clinically used. These difficulties mostly pertain to the variability of the imaging features as a function of the acquisition device and protocol, the robustness of the models with respect to that variability, and the interpretation of the radiomic models. Addressing the impact of the variability in acquisition and reconstruction protocols is needed, as is harmonizing the radiomic feature calculation methods, to ensure the reproducibility of studies in a multicenter context and their implementation in a clinical workflow. In this review, we explain the potential impact of positron emission tomography radiomics for radiation therapy and underline the various aspects that need to be carefully addressed to make the most of this promising approach.