Reproducibility of tumor uptake heterogeneity characterization through textural feature analysis in 18F-FDG PET

Potential and limitations of radiomics in neuro-oncology

Radiomics seeks to apply classical methods of image processing to obtain quantitative parameters from imaging. Derived features are subsequently fed into algorithmic models to aid clinical decision making. The application of radiomics and machine learning techniques to clinical medicine remains in its infancy. The great potential of radiomics lies in its objective, granular approach to investigating clinical imaging. In neuro-oncology, advanced machine learning techniques, particularly deep learning, are at the forefront of new discoveries in the field. However, despite the great promise of machine learning aided radiomic approaches, the current use remains confined to scholarly research, without real-world deployment in neuro-oncology. The paucity of data, inconsistencies in preprocessing, radiomic feature instability, and the rarity of the events of interest are critical barriers to clinical translation. In this article, we will outline the major steps in the process of radiomics, as well as review advances and challenges in the field as they pertain to neuro-oncology.

Repeatability of 18F-FDG PET Radiomic Features in Cervical Cancer

Knowledge of the intrinsic variability of radiomic features is essential to the proper interpretation of changes in these features over time. The primary aim of this study was to assess the test-retest repeatability of radiomic features extracted from 18F-FDG PET images of cervical tumors. The impact of different image preprocessing methods was also explored. Methods: Patients with cervical cancer underwent baseline and repeat 18F-FDG PET/CT imaging within 7 d. PET images were reconstructed using 2 methods: ordered-subset expectation maximization (PETOSEM) or ordered-subset expectation maximization with point-spread function (PETPSF). Tumors were segmented to produce whole-tumor volumes of interest (VOIWT) and 40% isocontours (VOI40). Voxels were either left at the default size or resampled to 3-mm isotropic voxels. SUV was discretized to a fixed number of bins (32, 64, or 128). Radiomic features were extracted from both VOIs, and repeatability was then assessed using the Lin concordance correlation coefficient (CCC). Results: Eleven patients were enrolled and completed the test-retest PET/CT imaging protocol. Shape, neighborhood gray-level difference matrix, and gray-level cooccurrence matrix features were repeatable, with a mean CCC value of 0.81. Radiomic features extracted from PETOSEM images showed significantly better repeatability than features extracted from PETPSF images (P < 0.001). Radiomic features extracted from VOI40 were more repeatable than features extracted from VOIWT (P < 0.001). For most features (78.4%), a change in bin number or voxel size resulted in less than a 10% change in feature value. All gray-level emphasis and gray-level run emphasis features showed poor repeatability (CCC values < 0.52) when extracted from VOIWT but were highly repeatable (mean CCC values > 0.96) when extracted from VOI40Conclusion: Shape, gray-level cooccurrence matrix, and neighborhood gray-level difference matrix radiomic features were consistently repeatable, whereas gray-level run length matrix and gray-level zone length matrix features were highly variable. Radiomic features extracted from VOI40 were more repeatable than features extracted from VOIWT Changes in voxel size or SUV discretization parameters typically resulted in relatively small differences in feature value, though several features were highly sensitive to these changes.

Relationship between tumor heterogeneity and volume in cervical cancer: Evidence from integrated fluorodeoxyglucose 18 PET/MR texture analysis

OBJECTIVE

The aim of this study was to evaluate the effect of cervical cancer volume on PET/magnetic resonance (MR) texture heterogeneity.

MATERIALS AND METHODS

We retrospectively analyzed the PET/MR images of 138 patients with pathologically diagnosed cervical squamous cell carcinoma, including 50 patients undergoing surgery and 88 patients receiving concurrent chemoradiotherapy. Fluorodeoxyglucose 18 (18FDG)-PET/MR examination were performed for each patient before treatment, and the PET and MR texture analysis were undertaken. The texture features of the tumor based on gray-level co-occurrence matrices were extracted, and the correlation between tumor texture features and volume parameters was analyzed using Spearman's rank correlation coefficient. Finally, the variation trend of tumor texture heterogeneity was analyzed as tumor volumes increased.

RESULTS

PET texture features were highly correlated with metabolic tumor volume (MTV), including entropy-log2, entropy-log10, energy, homogeneity, dissimilarity, contrast, correlation, and the correlation coefficients (rs) were 0.955, 0.955, -0.897, 0.883, -0.881, -0.876, and 0.847 (P < 0.001), respectively. In the range of smaller MTV, the texture heterogeneity of energy, entropy-log2, and entropy-log10 increases with an increase in tumor volume, whereas the texture heterogeneity of homogeneity, dissimilarity, contrast, and correlation decreases with an increase in tumor volume. Only homogeneity, contrast, correlation, and dissimilarity had high correlation with tumor volume on MRI. The correlation coefficients (rs) were 0.76, -0.737, 0.644, and -0.739 (P < 0.001), respectively. The texture heterogeneity of MRI features that are highly correlated with tumor volume decreases with increasing tumor volume.

CONCLUSION

In the small tumor volume range, the heterogeneity variation trend of PET texture features is inconsistent as the tumor volume increases, but the variation trend of MRI texture heterogeneity is consistent, and MRI texture heterogeneity decreases as tumor volume increases. These results suggest that MRI is a better imaging modality when compared with PET in determining tumor texture heterogeneity in the small tumor volume range.

Considerable effects of imaging sequences, feature extraction, feature selection, and classifiers on radiomics-based prediction of microvascular invasion in hepatocellular carcinoma using magnetic resonance imaging

BACKGROUND

Microvascular invasion (MVI) has a significant effect on the prognosis of hepatocellular carcinoma (HCC), but its preoperative identification is challenging. Radiomics features extracted from medical images, such as magnetic resonance (MR) images, can be used to predict MVI. In this study, we explored the effects of different imaging sequences, feature extraction and selection methods, and classifiers on the performance of HCC MVI predictive models.

METHODS

After screening against the inclusion criteria, 69 patients with HCC and preoperative gadoxetic acid-enhanced MR images were enrolled. In total, 167 features were extracted from the MR images of each sequence for each patient. Experiments were designed to investigate the effects of imaging sequence, number of gray levels (Ng), quantization algorithm, feature selection method, and classifiers on the performance of radiomics biomarkers in the prediction of HCC MVI. We trained and tested these models using leave-one-out cross-validation (LOOCV).

RESULTS

The radiomics model based on the images of the hepatobiliary phase (HBP) had better predictive performance than those based on the arterial phase (AP), portal venous phase (PVP), and pre-enhanced T1-weighted images [area under the receiver operating characteristic (ROC) curve (AUC) =0.792 vs. 0.641/0.634/0.620, P=0.041/0.021/0.010, respectively]. Compared with the equal-probability and Lloyd-Max algorithms, the radiomics features obtained using the Uniform quantization algorithm had a better performance (AUC =0.643/0.666 vs. 0.792, P=0.002/0.003, respectively). Among the values of 8, 16, 32, 64, and 128, the best predictive performance was achieved when the Ng was 64 (AUC =0.792 vs. 0.584/0.697/0.677/0.734, P<0.001/P=0.039/0.001/0.137, respectively). We used a two-stage feature selection method which combined the least absolute shrinkage and selection operator (LASSO) and recursive feature elimination (RFE) gradient boosting decision tree (GBDT), which achieved better stability than and outperformed LASSO, minimum redundancy maximum relevance (mRMR), and support vector machine (SVM)-RFE (stability =0.967 vs. 0.837/0.623/0.390, respectively; AUC =0.850 vs. 0.792/0.713/0.699, P=0.142/0.007/0.003, respectively). The model based on the radiomics features of HBP images using the GBDT classifier showed a better performance for the preoperative prediction of MVI compared with logistic regression (LR), SVM, and random forest (RF) classifiers (AUC =0.895 vs. 0.850/0.834/0.884, P=0.558/0.229/0.058, respectively). With the optimal combination of these factors, we established the best model, which had an AUC of 0.895, accuracy of 87.0%, specificity of 82.5%, and sensitivity of 93.1%.

CONCLUSIONS

Imaging sequences, feature extraction and selection methods, and classifiers can have a considerable effect on the predictive performance of radiomics models for HCC MVI.

Impact of segmentation and discretization on radiomic features in 68Ga-DOTA-TOC PET/CT images of neuroendocrine tumor

OBJECTIVE

To identify the impact of segmentation methods and intensity discretization on radiomic features (RFs) extraction from 68Ga-DOTA-TOC PET images in patients with neuroendocrine tumors.

METHODS

Forty-nine patients were retrospectively analyzed. Tumor contouring was performed manually by four different operators and with a semi-automatic edge-based segmentation (SAEB) algorithm. Three SUVmax fixed thresholds (20, 30, 40%) were applied. Fifty-one RFs were extracted applying two different intensity rescale factors for gray-level discretization: one absolute (AR60 = SUV from 0 to 60) and one relative (RR = min-max of the VOI SUV). Dice similarity coefficient (DSC) was calculated to quantify segmentation agreement between different segmentation methods. The impact of segmentation and discretization on RFs was assessed by intra-class correlation coefficients (ICC) and the coefficient of variance (COVL). The RFs' correlation with volume and SUVmax was analyzed by calculating Pearson's correlation coefficients.

RESULTS

DSC mean value was 0.75 ± 0.11 (0.45-0.92) between SAEB and operators and 0.78 ± 0.09 (0.36-0.97), among the four manual segmentations. The study showed high robustness (ICC > 0.9): (a) in 64.7% of RFs for segmentation methods using AR60, improved by applying SUVmax threshold of 40% (86.5%); (b) in 50.9% of RFs for different SUVmax thresholds using AR60; and (c) in 37% of RFs for discretization settings using different segmentation methods. Several RFs were not correlated with volume and SUVmax.

CONCLUSIONS

RFs robustness to manual segmentation resulted higher in NET 68Ga-DOTA-TOC images compared to 18F-FDG PET/CT images. Forty percent SUVmax thresholds yield superior RFs stability among operators, however leading to a possible loss of biological information. SAEB segmentation appears to be an optimal alternative to manual segmentation, but further validations are needed. Finally, discretization settings highly impacted on RFs robustness and should always be stated.

Magnetic resonance imaging histogram analysis of corpus callosum in a functional neurological disorder

Background/aim

The aim of the present study was to examine and compare the corpus callosum (CC) via histogram analysis (HA) on T1-weighted MR images of patients diagnosed with Functional Neurological Disorder (FND) and healthy controls.

Materials and methods

The study group included 19 female patients diagnosed with FND, and the control group included 20 healthy subjects. All participants were scanned with a 1.5 T MR scanner. A high-resolution structural image of the entire brain was obtained with sagittal 3D spiral fast spin echo T1-weighted images. Gray level intensity, standard deviation of the histogram, entropy, uniformity, skewness, and kurtosis values were determined with texture analysis. A student’s t-test was used to compare the group data. P < 0.05 was accepted as statistically significant.

Results

It was determined that the mean gray level intensity, standard deviation of the histogram, entropy calculated by the maximum, median and variance and size M percentage values were higher in patients with FND. Kurtosis and size U percentages values were lower in patients with FND.

Conclusion

In the present study, analysis of CC with T1-weighted MR image HA demonstrated significant differences between FND patients and healthy controls. The study findings indicated that HA is a beneficial technique for demonstrating textural variations between the CCs of patients with FND and healthy controls using MR images.

Role of Multimodality Imaging in Gastroesophageal Reflux Disease and Its Complications, with Clinical and Pathologic Correlation

Gastroesophageal reflux disease (GERD) is a common condition and impairs the quality of life for millions of patients, accounts for considerable health care spending, and is a primary risk factor for esophageal adenocarcinoma. There have been substantial advances in understanding the pathogenesis of GERD and its complications and much progress in diagnosis and management of GERD; however, these have not been comprehensively discussed in the recent radiology literature. Understanding the role of imaging in GERD and its complications is important to aid in multidisciplinary treatment of GERD. GERD results from prolonged or recurrent reflux of gastric contents into the esophagus. Common symptoms include heartburn or regurgitation. Prolonged reflux of gastric contents into the esophagus can cause erosive esophagitis. Over time, the inflammatory response related to esophagitis can lead to deposition of fibrous tissue and development of strictures. Alternatively, the esophageal mucosa can undergo metaplasia (Barrett esophagus), a precursor to dysplasia (which can lead to adenocarcinoma). Conventional barium esophagography has long been considered the primary imaging modality for the esophagus, and the fluoroscopic findings for diagnosis of GERD have been well established. Multimodality imaging has a clear role in detection and assessment of the complications of GERD, specifically reflux esophagitis and Barrett esophagus; differentiation of benign and malignant strictures; and detection, staging, and posttreatment surveillance of esophageal adenocarcinoma. Given the dramatic changes in utilization of abdominal imaging during the past 2 decades, with significantly declining volume of fluoroscopic procedures and concomitant increase in CT and MRI studies, it is crucial that modern radiologists appreciate the value of barium esophagography in the workup of GERD and recognize the key imaging features of GERD and its complications at CT and MRI.

Prediction of local recurrence and distant metastasis using radiomics analysis of pretreatment nasopharyngeal [18F]FDG PET/CT images

OBJECTIVES

To develop a radiomics signature to predict locoregional recurrence (LR) and distant metastasis (DM), as extracted from pretreatment 2-deoxy-2-[¹⁸F]fluoro-D-glucose ([18F]FDG) positron emission tomography/X-ray computed tomography (PET/CT) images in locally advanced nasopharyngeal carcinoma (NPC).

MATERIALS AND METHODS

Eighty-five patients with Stage III-IVB NPC underwent pretreatment [¹⁸F]FDG PET/CT scans and received radiotherapy or chemoradiotherapy. 53 of them achieved disease control, and 32 of them failed after treatment (15: LR, 17: DM). A total of 114 radiomic features were extracted from PET/CT images. For univariate analysis, Wilcoxon test and Chi-square test were used to compare median values of features between different treatment outcomes and predict the risk of treatment failure, respectively. For multivariate analysis, all features were grouped into clusters based on Pearson correlation using hierarchical clustering, and the representative feature of each cluster was chosen by the Relief algorithm. Then sequential floating forward selection (SFFS) coupled with a support vector machine (SVM) classifier were used to derive the optimized feature set in terms of the area under receiver operating characteristic (ROC) curve (AUC). The performance of the model was evaluated by leave-one-out-cross-validation, fivefold cross-validation, tenfold cross-validation.

RESULTS

Twenty features had significant differences between disease control and treatment failure. NPC patients with values of Compactness1, Compactness2, Coarseness_NGTDM or SGE_GLGLM above the median as well as patients with values of Irregularity, RLN_GLRLM or GLV_GLSZM below the median, showed a significant (p < 0.05) higher risk of treatment failure (about 50% vs. 25%). The derived radiomics signature consisted of 5 features with the highest AUC value of 0.8290 (sensitivity: 0.8438, specificity: 0.7736) using leave-one-out-cross-validation.

CONCLUSION

Locoregional recurrence (LR) and DM of locally advanced NPC can be predicted using radiomics analysis of pretreatment [¹⁸F]FDG PET/CT. The SFFS feature selection coupled with SVM classifier can derive the optimized feature set with correspondingly highest AUC value for pretreatment prediction of LR and/or DM of NPC.

Tumor heterogeneity for differentiation between liver tumors and normal liver tissue in 18F-FDG PET/CT

AIM

 Malignancies show higher spatial heterogeneity than normal tissue. We investigated, if textural parameters from FDG PET describing the heterogeneity function as tool to differentiate between tumor and normal liver tissue.

METHODS

 FDG PET/CT scans of 80 patients with liver metastases and 80 patients with results negative upper abdominal organs were analyzed. Metastases and normal liver tissue were analyzed drawing up to three VOIs with a diameter of 25 mm in healthy liver tissue of the tumoral affected and results negative liver, whilst up to 3 metastases per patient were delineated. Within these VOIs 30 different textural parameters were calculated as well as SUV. The parameters were compared in terms of intra-patient and inter-patient variability (2-sided t test). ROC analysis was performed to analyze predictive power and cut-off values.

RESULTS

 28 textural parameters differentiated healthy and pathological tissue (p < 0.05) with high sensitivity and specificity. SUV showed ability to differentiate but with a lower significance. 15 textural parameters as well as SUV showed a significant variation between healthy tissues out of tumour infested and negative livers. Mean intra- and inter-patient variability of metastases were found comparable or lower for 6 of the textural features than the ones of SUV. They also showed good values of mean intra- and inter-patient variability of VOIs drawn in liver tissue of patients with metastases and of results negative ones.

CONCLUSION

 Heterogeneity parameters assessed in FDG PET are promising to classify tissue and differentiate malignant lesions usable for more personalized treatment planning, therapy response evaluation and precise delineation of tumors for target volume determination as part of radiation therapy planning.

MRI texture feature repeatability and image acquisition factor robustness, a phantom study and in silico study

Purpose

To evaluate robustness and repeatability of magnetic resonance imaging (MRI) texture features in water and tissue phantom test-retest study.

Materials and methods

Separate water and tissue phantoms were imaged twice with the same protocol in a test-retest experiment using a 1.5-T scanner. Protocols were acquired to favour signal-to-noise ratio and resolution. Forty-six features including first order statistics and second-order texture features were extracted, and repeatability was assessed by calculating the concordance correlation coefficient. Separately, base image noise and resolution were manipulated in an in silico experiment, and robustness of features was calculated by assessing percentage coefficient of variation and linear correlation of features with noise and resolution. These simulation data were compared with the acquired data. Features were classified by their degree (high, intermediate, or low) of robustness and repeatability.

Results

Eighty percent of the MRI features were repeatable (concordance correlation coefficient > 0.9) in the phantom test-retest experiment. The majority (approximately 90%) demonstrated a strong or intermediate correlation with image acquisition parameter, and 19/46 (41%) and 13/46 (28%) of features were highly robust to noise and resolution, respectively (coefficient of variation < 5%). Agreement between the acquired and simulation data varied, with the range of agreement within feature classes between 11 and 92%.

Conclusion

Most MRI features were repeatable in a phantom test-retest study. This phantom data may serve as a lower limit of feature MRI repeatability. Robustness of features varies with acquisition parameter, and appropriate features can be selected for clinical validation studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s41747-020-00199-6.

A Phantom Study to Investigate Robustness and Reproducibility of Grey Level Co-Occurrence Matrix (GLCM)-Based Radiomics Features for PET

Quantification and classification of heterogeneous radiotracer uptake in Positron Emission Tomography (PET) using textural features (termed as radiomics) and artificial intelligence (AI) has the potential to be used as a biomarker of diagnosis and prognosis. However, textural features have been predicted to be strongly correlated with volume, segmentation and quantization, while the impact of image contrast and noise has not been assessed systematically. Further continuous investigations are required to update the existing standardization initiatives. This study aimed to investigate the relationships between textural features and these factors with 18F filled torso NEMA phantom to yield different contrasts and reconstructed with different durations to represent varying levels of noise. The phantom was also scanned with heterogeneous spherical inserts fabricated with 3D printing technology. All spheres were delineated using: (1) the exact boundaries based on their known diameters; (2) 40% fixed; and (3) adaptive threshold. Six textural features were derived from the gray level co-occurrence matrix (GLCM) using different quantization levels. The results indicate that homogeneity and dissimilarity are the most suitable for measuring PET tumor heterogeneity with quantization 64 provided that the segmentation method is robust to noise and contrast variations. To use these textural features as prognostic biomarkers, changes in textural features between baseline and treatment scans should always be reported along with the changes in volumes.

Radiomics and Digital Image Texture Analysis in Oncology (Review)

One of the most promising areas of diagnosis and prognosis of diseases is radiomics, a science combining radiology, mathematical modeling, and deep machine learning. The main concept of radiomics is image biomarkers (IBMs), the parameters characterizing various pathological changes and calculated based on the analysis of digital image texture. IBMs are used for quantitative assessment of digital imaging results (CT, MRI, ultrasound, PET). The use of IBMs in the form of "virtual biopsy" is of particular relevance in oncology. The article provides the basic concepts of radiomics identifying the main stages of obtaining IBMs: data collection and preprocessing, tumor segmentation, data detection and extraction, modeling, statistical processing, and data validation. The authors have analyzed the possibilities of using IBMs in oncology, describing the currently known features and advantages of using radiomics and image texture analysis in the diagnosis and prognosis of cancer. The limitations and problems associated with the use of radiomics data are considered. Although the novel effective tool for performing virtual biopsy of human tissue is at the development stage, quite a few projects have already been implemented, and medical software packages for radiomics analysis of digital images have been created.

Computed tomography-based deep-learning prediction of neoadjuvant chemoradiotherapy treatment response in esophageal squamous cell carcinoma

BACKGROUND

Deep learning is promising to predict treatment response. We aimed to evaluate and validate the predictive performance of the CT-based model using deep learning features for predicting pathologic complete response to neoadjuvant chemoradiotherapy (nCRT) in esophageal squamous cell carcinoma (ESCC).

MATERIALS AND METHODS

Patients were retrospectively enrolled between April 2007 and December 2018 from two institutions. We extracted deep learning features of six pre-trained convolutional neural networks, respectively, from pretreatment CT images in the training cohort (n = 161). Support vector machine was adopted as the classifier. Validation was performed in an external testing cohort (n = 70). We assessed the performance using the area under the receiver operating characteristics curve (AUC) and selected an optimal model, which was compared with a radiomics model developed from the training cohort. A clinical model consisting of clinical factors only was also built for baseline comparison. We further conducted a radiogenomics analysis using gene expression profiles to reveal underlying biology associated with radiological prediction.

RESULTS

The optimal model with features extracted from ResNet50 achieved an AUC and accuracy of 0.805 (95% CI, 0.696-0.913) and 77.1% (65.6%-86.3%) in the testing cohort, compared with 0.725 (0.605-0.846)) and 67.1% (54.9%-77.9%) for the radiomics model. All the radiological models showed better predictive performance than the clinical model. Radiogenomics analysis suggested a potential association mainly with WNT signaling pathway and tumor microenvironment.

CONCLUSIONS

The novel and noninvasive deep learning approach could provide efficient and accurate prediction of treatment response to nCRT in ESCC, and benefit clinical decision making of therapeutic strategy.

Radiomics Analysis of [18F]-Fluorodeoxyglucose-Avid Thyroid Incidentalomas Improves Risk Stratification and Selection for Clinical Assessment

Background: [¹⁸F]-Fluorodeoxyglucose (FDG)-avid thyroid lesions incidentally detected on positron emission tomography/computed tomography (PET/CT) scans represent a tumor lesion in about 30% of cases. The present study evaluated the ability of PET metrics and radiomics features to predict final diagnosis of [¹⁸F]FDG thyroid incidentalomas (TIs). Methods: A total of 104 patients with 107 TIs were retrospectively studied; 30 nodules (28%) were diagnosed as malignant. After volumetric segmentation of each thyroid lesion, metabolic tumor volume (MTV), total lesion glycolysis (TLG), standardized uptake values (SUVs), and metabolic heterogeneity were estimated, and 107 radiomics features were extracted following a standard protocol. Results: MTV, TLG, SUVmax, SUVmean, and SUVpeak among functional PET parameters, and gray-level co-occurrence matrix (GLCM)_InverseDifferenceMoment, shape_Sphericity, GLCM_SumSquares, firstorder_Maximum2DDiameterSlice, firstorder_Energy, and GLCM_Contrast among nonredundant radiomics features, showed significantly different values between malignant and benign TIs (Mann-Whitney U-test, p < 0.01 for all). Univariate logistic regression revealed that these parameters demonstrated good ability to predict final diagnosis of TIs (p < 0.02 for all). Shape_Sphericity was the best predictor classifying 82% of TIs correctly (p < 0.0001). Only TLG, SUVmax, and shape_Sphericity retained significance (p < 0.0001) by multivariate analysis. Malignant lesion prevalence increased from 7% to 100% in accordance with the number (score, 0-3) of the three positive parameters present (χ² trend, p < 0.0001). A score of 0 excludes malignant TIs with a negative predictive value of 93%, while a score of 3 predicted malignancy with a positive predictive value of 100%. Conclusions: PET metrics and radiomics analysis can improve identification of [¹⁸F]FDG-avid TIs at high risk of malignancy. A model based on TLG, SUVmax, and shape_Sphericity may allow prediction of a final diagnosis, providing useful information for the management of TIs.

Harmonization strategies for multicenter radiomics investigations

Carrying out large multicenter studies is one of the key goals to be achieved towards a faster transfer of the radiomics approach in the clinical setting. This requires large-scale radiomics data analysis, hence the need for integrating radiomic features extracted from images acquired in different centers. This is challenging as radiomic features exhibit variable sensitivity to differences in scanner model, acquisition protocols and reconstruction settings, which is similar to the so-called 'batch-effects' in genomics studies. In this review we discuss existing methods to perform data integration with the aid of reducing the unwanted variation associated with batch effects. We also discuss the future potential role of deep learning methods in providing solutions for addressing radiomic multicentre studies.

Is FDG-PET texture analysis related to intratumor biological heterogeneity in lung cancer?

OBJECTIVES

We aimed at investigating the origin of the correlations between tumor volume and ¹⁸F-FDG-PET texture indices in lung cancer.

METHODS

Eighty-five consecutive patients with newly diagnosed non-small cell lung cancer (NSCLC) underwent a ¹⁸F-FDG-PET/CT scan before treatment. Seven phantom spheres uniformly filled with ¹⁸F-FDG, and covering a range of activities and volumes similar to that found in lung tumors, were also scanned. Established texture indices were computed for lung tumors and homogeneous spheres. The dependence between textural indices and volume in homogeneous spheres was modeled and then used to predict texture indices in lung tumors. Correlation analyses were carried out between predicted and texture features measured in lung tumors. Cox proportional hazards regression was used to investigate the associations between overall survival and volume-adjusted textural features.

RESULTS

All textural features showed strong, non-linear correlations with volume, both in tumors and homogeneous spheres. Correlations between predicted versus measured texture features were very high for contrast (r² = 0.91), dissimilarity (r² = 0.90), ZP (r² = 0.90), GLNN (r² = 0.86), and homogeneity (r² = 0.82); high for entropy (r² = 0.50) and HILAE (r² = 0.53); and low for energy (r² = 0.30). Cox regressions showed that among volume-adjusted features, only HILAE was associated with overall survival (b = - 0.35, p = 0.008).

CONCLUSION

We have shown that texture indices previously found to be correlated with a number of clinically relevant outcomes might not provide independent information apart from that driven by their correlation with tumor volume, suggesting that these metrics might not be suitable as intratumor heterogeneity markers.

KEY POINTS

• Associations between texture FDG-PET indices and overall survival have been widely reported in lung cancer, with tumor volume also being associated with overall survival, and therefore, it is still unclear whether the predictive power of textural indices is simply driven by this correlation. • Our results demonstrated strong non-linear correlations between textural indices and volume, showing an analogous behavior for lung tumors from patients and homogeneous spheres inserted in phantoms. • Our findings showed that texture FDG-PET indices might not provide independent information apart from that driven by their correlation with tumor volume.

Radiomics in Breast Imaging from Techniques to Clinical Applications: A Review

Recent advances in computer technology have generated a new area of research known as radiomics. Radiomics is defined as the high throughput extraction and analysis of quantitative features from imaging data. Radiomic features provide information on the gray-scale patterns, inter-pixel relationships, as well as shape and spectral properties of radiological images. Moreover, these features can be used to develop computational models that may serve as a tool for personalized diagnosis and treatment guidance. Although radiomics is becoming popular and widely used in oncology, many problems such as overfitting and reproducibility issues remain unresolved. In this review, we will outline the steps of radiomics used for oncology, specifically addressing applications for breast cancer patients and focusing on technical issues.

Radiomics in PET/CT: Current Status and Future AI-Based Evolutions

This short review aims at providing the readers with an update on the current status, as well as future perspectives in the quickly evolving field of radiomics applied to the field of PET/CT imaging. Numerous pitfalls have been identified in study design, data acquisition, segmentation, features calculation and modeling by the radiomics community, and these are often the same issues across all image modalities and clinical applications, however some of these are specific to PET/CT (and SPECT/CT) imaging and therefore the present paper focuses on those. In most cases, recommendations and potential methodological solutions do exist and should therefore be followed to improve the overall quality and reproducibility of published studies. In terms of future evolutions, the techniques from the larger field of artificial intelligence (AI), including those relying on deep neural networks (also known as deep learning) have already shown impressive potential to provide solutions, especially in terms of automation, but also to maybe fully replace the tools the radiomics community has been using until now in order to build the usual radiomics workflow. Some important challenges remain to be addressed before the full impact of AI may be realized but overall the field has made striking advances over the last few years and it is expected advances will continue at a rapid pace.

A Role for FDG PET Radiomics in Personalized Medicine?

Radiomics describes the extraction of multiple features from medical images, including molecular imaging modalities, that with bioinformatic approaches, provide additional clinically relevant information that may be invisible to the human eye. This information may complement standard radiological interpretation with data that may better characterize a disease or that may provide predictive or prognostic information. Progressing from predefined image features, often describing heterogeneity of voxel intensities within a volume of interest, there is increasing use of machine learning to classify disease characteristics and deep learning methods based on artificial neural networks that can learn features without a priori definition and without the need for preprocessing of images. There have been advances in standardization and harmonization of methods to a level that should support multicenter studies. However, in this relatively early phase of research in the field, there are limited aspects that have been adopted into routine practice. Most of the reports in the molecular imaging field describe radiomic approaches in cancer using 18F-fluorodeoxyglucose positron emission tomography (18F-FDG-PET). In this review, we will describe radiomics in molecular imaging and summarize the pertinent literature in lung cancer where reports are most prevalent and mature.

Texture Feature Comparison Between Step-and-Shoot and Continuous-Bed-Motion 18F-FDG PET

Our objective was to investigate the differences in texture features between step-and-shoot (SS) and continuous-bed-motion (CBM) imaging in phantom and clinical studies. Methods: A National Electrical Manufacturers Association body phantom was filled with ¹⁸F-FDG solution at a sphere-to-background ratio of 4:1. SS and CBM were performed using the same acquisition duration, and the data were reconstructed using 3-dimensional ordered-subset expectation maximization with time-of-flight algorithms. Texture features were extracted using the software LIFEx. A volume of interest was delineated on the 22-, 28-, and 37-mm spheres with a threshold of 42% of the maximum SUV. The voxel intensities were discretized using 2 resampling methods, namely a fixed bin size and a fixed bin number discretization. The discrete resampling values were set to 64 and 128. In total, 31 texture features were calculated with gray-level cooccurrence matrix (GLCM), gray-level run length matrix, neighborhood gray-level different matrix, and gray-level zone length matrix. The texture features of the SS and CBM images were compared for all settings using the paired t test and the coefficient of variation. In a clinical study, 27 lesions from 20 patients were examined using the same acquisition and image processing as were used during the phantom study. The percentage difference (%Diff) and correlation between the texture features from SS and CBM images were calculated to evaluate agreement between the 2 scanning techniques. Results: In the phantom study, the 11 features exhibited no significant difference between SS and CBM images, and the coefficient of variation was no more than 10%, depending on resampling conditions, whereas entropy and dissimilarity from GLCM fulfilled the criteria for all settings. In the clinical study, the entropy and dissimilarity from GLCM exhibited a low %Diff and excellent correlation in all resampling conditions. The %Diff of entropy was lower than that of dissimilarity. Conclusion: Differences between the texture features of SS and CBM images varied depending on the type of feature. Because entropy for GLCM exhibits minimal differences between SS and CBM images irrespective of resampling conditions, entropy may be the optimal feature to reduce the differences between the 2 scanning techniques.

Training and validation of a robust PET radiomic-based index to predict distant-relapse-free-survival after radio-chemotherapy for locally advanced pancreatic cancer

PURPOSE

To assess the value of 18F-Fluorodeoxyglucose (18F-FDG) PET Radiomic Features (RF) in predicting Distant Relapse Free Survival (DRFS) in patients with Locally AdvancedPancreaticCancer (LAPC) treated with radio-chemotherapy.

MATERIALS & METHODS

One-hundred-ninety-eight RFs were extracted using IBSI (Image Biomarker Standardization Initiative) consistent software from pre-radiotherapy images of 176 LAPC patients treated with moderate hypo-fractionation (44.25 Gy, 2.95 Gy/fr). Tumors were segmented by applying a previously validated semi-automatic method. One-hundred-twenty-six RFs were excluded due to poor reproducibility and/or repeatability and/or inter-scanner variability. The original cohort was randomly split into a training (n = 116) and a validation (n = 60) group. Multi-variable Cox regression was applied to the training group, including only independent RFs in the model. The resulting radiomic index was tested in the validation cohort. The impact of selected clinical variables was also investigated.

RESULTS

The resulting Cox model included two first order RFs: Center of Mass Shift (COMshift) and 10th Intensity percentile (P₁₀) (p = 0.0005, HR = 2.72, 95%CI = 1.54-4.80), showing worse outcomes for patients with lower COMshift and higher P₁₀. Once stratified by quartile values (<lowest quartile vs >highest quartile vs the remaining), the index properly stratified patients according to their DRFS (p = 0.0024, log-rank test). Performances were confirmed in the validation cohort (p = 0.03, HR = 2.53, 95%CI = 0.96-6.65). The addition of clinical factors did not significantly improve the models' performance.

CONCLUSIONS

A radiomic-based index including only two robust PET-RFs predicted DRFS of LAPC patients after radio-chemotherapy. The current results could find relevant applications in the treatment personalization of LAPC. A multi-institution independent validation has been planned.

Radiomics for precision medicine: Current challenges, future prospects, and the proposal of a new framework

The advancement of artificial intelligence concurrent with the development of medical imaging techniques provided a unique opportunity to turn medical imaging from mostly qualitative, to further quantitative and mineable data that can be explored for the development of clinical decision support systems (cDSS). Radiomics, a method for the high throughput extraction of hand-crafted features from medical images, and deep learning -the data driven modeling techniques based on the principles of simplified brain neuron interactions, are the most researched quantitative imaging techniques. Many studies reported on the potential of such techniques in the context of cDSS. Such techniques could be highly appealing due to the reuse of existing data, automation of clinical workflows, minimal invasiveness, three-dimensional volumetric characterization, and the promise of high accuracy and reproducibility of results and cost-effectiveness. Nevertheless, there are several challenges that quantitative imaging techniques face, and need to be addressed before the translation to clinical use. These challenges include, but are not limited to, the explainability of the models, the reproducibility of the quantitative imaging features, and their sensitivity to variations in image acquisition and reconstruction parameters. In this narrative review, we report on the status of quantitative medical image analysis using radiomics and deep learning, the challenges the field is facing, propose a framework for robust radiomics analysis, and discuss future prospects.

Towards guidelines to harmonize textural features in PET: Haralick textural features vary with image noise, but exposure-invariant domains enable comparable PET radiomics

Purpose

Image texture is increasingly used to discriminate tissues and lesions in PET/CT. For quantification or in computer-aided diagnosis, textural feature analysis must produce robust and comparable values. Because statistical feature values depend on image count statistics, we investigated in depth the stability of Haralick features values as functions of acquisition duration, and for common image resolutions and reconstructions.

Methods

A homogeneous cylindrical phantom containing 9.6 kBq/ml Ge-68 was repeatedly imaged on a Siemens Biograph mCT, with acquisition durations ranging from three seconds to three hours. Images with 1.5, 2, and 4 mm isometrically spaced voxels were reconstructed with filtered back-projection (FBP), ordered subset expectation maximization (OSEM), and the Siemens TrueX algorithm. We analysed Haralick features derived from differently quantized (3 to 8-bit) grey level co-occurrence matrices (GLCMs) as functions of exposure E, which we defined as the product of activity concentration in a volume of interest (VOI) and acquisition duration. The VOI was a 50 mm wide cube at the centre of the phantom. Feature stability was defined for df/dE → 0.

Results

The most stable feature values occurred in low resolution FBPs, whereas some feature values from 1.5 mm TrueX reconstructions ranged over two orders of magnitude. Within the same reconstructions, most feature value-exposure curves reached stable plateaus at similar exposures, regardless of GLCM quantization. With 8-bit GLCM, median time to stability was 16 s and 22 s for FBPs, 18 s and 125 s for OSEM, and 23 s, 45 s, and 76 s for PSF reconstructions, with longer durations for higher resolutions. Stable exposures coincided in OSEM and TrueX reconstructions with image noise distributions converging to a Gaussian. In FBP, the occurrence of stable values coincided the disappearance of negatives image values in the VOI.

Conclusions

Haralick feature values depend strongly on exposure, but invariance exists within defined domains of exposure. Here, we present an easily replicable procedure to identify said stable exposure domains, where image noise does not substantially add to textural feature values. Only by imaging at predetermined feature-invariant exposure levels and by adjusting exposure to expected activity concentrations, can textural features have a quantitative use in PET/CT. The necessary exposure levels are attainable by modern PET/CT systems in clinical routine.

Reliability of CT radiomic features reflecting tumour heterogeneity according to image quality and image processing parameters

The reliability of radiomics features (RFs) is crucial for quantifying tumour heterogeneity. We assessed the influence of imaging, segmentation, and processing conditions (quantization range, bin number, signal-to-noise ratio [SNR], and unintended outliers) on RF measurement. Low SNR and unintended outliers increased the standard deviation and mean values of histograms to calculate the first-order RFs. Variations in imaging processing conditions significantly altered the shape of the probability distribution (centre of distribution, extent of dispersion, and segmentation of probability clusters) in second-order RF matrices (i.e. grey-level co-occurrence and grey-level run length), thereby eventually causing fluctuations in RF estimation. Inconsistent imaging and processing conditions decreased the number of reliably measured RFs in terms of individual RF values (intraclass correlation coefficient ≥0.75) and inter-lesion RF ratios (coefficient of variation <15%). No RF could be reliably estimated under inconsistent SNR and inclusion of outlier conditions. By contrast, with high SNR and no outliers, all first-order RFs, 11 (42%) grey-level co-occurrence RFs and five (42%) grey-level run length RFs showed acceptable reliability. Our study suggests that optimization of SNR, exclusion of outliers, and application of relevant quantization range and bin number should be performed to ensure the robustness of radiomics studies assessing tumor heterogeneity.

Heterogeneity in tumours: Validating the use of radiomic features on 18F-FDG PET/CT scans of lung cancer patients as a prognostic tool

AIM

The aim was to validate promising radiomic features (RFs)1 on 18F-flourodeoxyglucose positron emission tomography/computed tomography-scans (18F-FDG PET/CT) of non-small cell lung cancer (NSCLC) patients undergoing definitive chemo-radiotherapy.

METHODS

18F-FDG PET/CT scans performed for radiotherapy (RT) planning were retrieved. Auto-segmentation with visual adaption was used to define the primary tumour on PET images. Six pre-selected prognostic and reproducible PET texture -and shape-features were calculated using texture respectively shape analysis. The correlation between these RFs and metabolic active tumour volume (MTV)3, gross tumour volume (GTV)4 and maximum and mean of standardized uptake value (SUV)5 was tested with a Spearman's Rank test. The prognostic value of RFs was tested in a univariate cox regression analysis and a multivariate cox regression analysis with GTV, clinical stage and histology. P-value ≤ 0.05 were considered significant.

RESULTS

Image analysis was performed for 233 patients: 145 males and 88 females, mean age of 65.7 and clinical stage II-IV. Mean GTV was 129.87 cm3 (SD 130.30 cm3). Texture and shape-features correlated more strongly to MTV and GTV compared to SUV-measurements. Four RFs predicted PFS in the univariate analysis. No RFs predicted PFS in the multivariate analysis, whereas GTV and clinical stage predicted PFS (p = 0.001 and p = 0.008 respectively).

CONCLUSION

The pre-selected RFs were insignificant in predicting PFS in combination with GTV, clinical stage and histology. These results might be due to variations in technical parameters. However, it is relevant to question whether RFs are stable enough to provide clinically useful information.

Differential Diagnosis of Nonhypervascular Pancreatic Neuroendocrine Neoplasms From Pancreatic Ductal Adenocarcinomas, Based on Computed Tomography Radiological Features and Texture Analysis

RATIONALE AND OBJECTIVES

To determine computed tomography (CT) radiological features and texture features that are rewarding in differentiating nonhypervascular pancreatic neuroendocrine neoplasms (PNENs) from pancreatic ductal adenocarcinomas (PDACs).

MATERIALS AND METHODS

We compared patients to pathologically proven nonhypervascular PNENs and age-matched controls with pathologically proven PDACs in a 1:2 ratio. Preoperative CT images in the arterial phase (AP) and portal vein phase (PVP) were obtained. Two radiologists independently reviewed the morphological characteristics of each tumor. Three-dimensional regions of interest (ROIs), drawn using ITK-SNAP software, were input into AK software (Artificial Intelligent Kit, GE) to extract texture features from AP and PVP images. Differences between PNENs and PDACs were analyzed with the chi-squared test, least absolute shrinkage and selection operator, kappa statistics, and uni- and multivariate logistic regression analyses.

RESULTS

In total, 40 nonhypervascular PNENs and 80 PDACs were evaluated. Maximum diameter on axial section, margin, calcification, vascularity in the tumor, and tumor heterogeneity were significantly different between PDACs and nonhypervascular PNENs. Multivariate analysis showed well-defined tumor margin (odds ratio: 21.0) and presence of calcification (odds ratio: 4.4) were significant predictors of nonhypervascular PNENs. The area under the receiver operating characteristic curve of the radiological feature model, AP texture model, and PVP texture model were 0.780, 0.855, and 0.929, respectively, based on logistic regression.

CONCLUSION

A well-defined margin and calcification in the tumor were helpful in discriminating nonhypervascular PNENs from PDACs. Texture analysis of contrast-enhanced CT images could be beneficial in differentially diagnosing nonhypervascular PNENs and PDACs.

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.

The Emerging Role of Radiomics in COPD and Lung Cancer

Medical imaging plays a key role in evaluating and monitoring lung diseases such as chronic obstructive pulmonary disease (COPD) and lung cancer. The application of artificial intelligence in medical imaging has transformed medical images into mineable data, by extracting and correlating quantitative imaging features with patients' outcomes and tumor phenotype - a process termed radiomics. While this process has already been widely researched in lung oncology, the evaluation of COPD in this fashion remains in its infancy. Here we outline the main applications of radiomics in lung cancer and briefly review the workflow from image acquisition to the evaluation of model performance. Finally, we discuss the current assessments of COPD and the potential application of radiomics in COPD.

Assessment of Stability and Discrimination Capacity of Radiomic Features on Apparent Diffusion Coefficient Images

The objectives of the study are to develop a new way to assess stability and discrimination capacity of radiomic features without the need of test-retest or multiple delineations and to use information obtained to perform a preliminary feature selection. Apparent diffusion coefficient (ADC) maps were computed from diffusion-weighted magnetic resonance images (DW-MRI) of two groups of patients: 18 with soft tissue sarcomas (STS) and 18 with oropharyngeal cancers (OPC). Sixty-nine radiomic features were computed, using three different histogram discretizations (16, 32, and 64 bins). Geometrical transformations (translations) of increasing entity were applied to the regions of interest (ROIs), and the intra-class correlation coefficient (ICC) was used to compare the features computed on the original and modified ROIs. The distribution of ICC values for minimal and maximal entity translations (ICC₁₀ and ICC₁₀₀, respectively) was used to adjust thresholds of ICC (ICC_min and ICC_max) used to discriminate between good, unstable (ICC₁₀ < ICC_min), and non-discriminative features (ICC₁₀₀ > ICC_max). Fifty-four and 59 radiomic features passed the stability-based selection for all the three histogram discretizations for the OPC and STS datasets, respectively. The excluded features were similar across the different histogram discretizations (Jaccard’s index 0.77 ± 0.13 and 0.9 ± 0.1 for OPC and STS, respectively) but different between datasets (Jaccard’s index 0.19 ± 0.02). The results suggest that the observed radiomic features are mainly stable and discriminative, but the stability depends on the region of the body under observation. The method provides a way to assess stability without the need of test-retest or multiple delineations.

Can Laws Be a Potential PET Image Texture Analysis Approach for Evaluation of Tumor Heterogeneity and Histopathological Characteristics in NSCLC?

We investigated the association between the textural features obtained from ¹⁸F-FDG images, metabolic parameters (SUVmax_, SUVmean, MTV, TLG), and tumor histopathological characteristics (stage and Ki-67 proliferation index) in non-small cell lung cancer (NSCLC). The FDG-PET images of 67 patients with NSCLC were evaluated. MATLAB technical computing language was employed in the extraction of 137 features by using first order statistics (FOS), gray-level co-occurrence matrix (GLCM), gray-level run length matrix (GLRLM), and Laws' texture filters. Textural features and metabolic parameters were statistically analyzed in terms of good discrimination power between tumor stages, and selected features/parameters were used in the automatic classification by k-nearest neighbors (k-NN) and support vector machines (SVM). We showed that one textural feature (gray-level nonuniformity, GLN) obtained using GLRLM approach and nine textural features using Laws' approach were successful in discriminating all tumor stages, unlike metabolic parameters. There were significant correlations between Ki-67 index and some of the textural features computed using Laws' method (r = 0.6, p = 0.013). In terms of automatic classification of tumor stage, the accuracy was approximately 84% with k-NN classifier (k = 3) and SVM, using selected five features. Texture analysis of FDG-PET images has a potential to be an objective tool to assess tumor histopathological characteristics. The textural features obtained using Laws' approach could be useful in the discrimination of tumor stage.

Prognostic Value of Functional Parameters of 18F-FDG-PET Images in Patients with Primary Renal/Adrenal Lymphoma

Objectives

The aim of this study is to explore the textural features that may identify the morphological changes in the lymphoma region and predict the prognosis of patients with primary renal lymphoma (PRL) and primary adrenal lymphoma (PAL).

Methods

This retrospective study comprised nineteen non-Hodgkin's lymphoma (NHL) patients undergoing ¹⁸F-FDG-PET/CT at West China Hospital from December 2013 to May 2017. ¹⁸F-FDG-PET images were reviewed independently by two board certificated radiologists of nuclear medicine, and the texture features were extracted from LifeX packages. The prognostic value of PET FDG-uptake parameters, patients' baseline characteristics, and textural parameters were analyzed using Kaplan–Meier analysis. Cox regression analysis was used to identify the independent prognostic factors among the imaging and clinical features.

Results

The overall survival of included patients was 18.84 ± 13.40 (mean ± SD) months. Univariate Cox analyses found that the tumor stage, GLCM (gray-level co-occurrence matrix) entropy, GLZLM_GLNU (gray-level nonuniformity), and GLZLM_ZLNU (zone length nonuniformity), values were significant predictors for OS. Among them, GLRLM_RLNU ≥216.6 demonstrated association with worse OS at multivariate analysis (HR 9.016, 95% CI 1.041–78.112, p=0.046).

Conclusions

The texture analysis of ¹⁸F-FDG-PET images could potentially serve as a noninvasive strategy to predict the overall survival of patients with PRL and PAL.

AI-based applications in hybrid imaging: how to build smart and truly multi-parametric decision models for radiomics

INTRODUCTION

The quantitative imaging features (radiomics) that can be obtained from the different modalities of current-generation hybrid imaging can give complementary information with regard to the tumour environment, as they measure different morphologic and functional imaging properties. These multi-parametric image descriptors can be combined with artificial intelligence applications into predictive models. It is now the time for hybrid PET/CT and PET/MRI to take the advantage offered by radiomics to assess the added clinical benefit of using multi-parametric models for the personalized diagnosis and prognosis of different disease phenotypes.

OBJECTIVE

The aim of the paper is to provide an overview of current challenges and available solutions to translate radiomics into hybrid PET-CT and PET-MRI imaging for a smart and truly multi-parametric decision model.

Radiogenomics-based cancer prognosis in colorectal cancer

Radiogenomics aims at investigating the relationship between imaging radiomic features and gene expression alterations. This study addressed the potential prognostic complementary value of contrast enhanced computed tomography (CE-CT) radiomic features and gene expression data in primary colorectal cancers (CRC). Sixty-four patients underwent CT scans and radiomic features were extracted from the delineated tumor volume. Gene expression analysis of a small set of genes, previously identified as relevant for CRC, was conducted on surgical samples from the same tumors. The relationships between radiomic and gene expression data was assessed using the Kruskal–Wallis test. Multiple testing was not performed, as this was a pilot study. Cox regression was used to identify variables related to overall survival (OS) and progression free survival (PFS). ABCC2 gene expression was correlated with N (p = 0.016) and M stages (p = 0.022). Expression changes of ABCC2, CD166, CDKNV1 and INHBB genes exhibited significant correlations with some radiomic features. OS was associated with Ratio 3D Surface/volume (p = 0.022) and ALDH1A1 expression (p = 0.042), whereas clinical stage (p = 0.004), ABCC2 expression (p = 0.035), and Entropy_{GLCM_E} (p = 0.0031), were prognostic factors for PFS. Combining CE-CT radiomics with gene expression analysis and histopathological examination of primary CRC could provide higher prognostic stratification power, leading to improved patient management.

Reproducibility and Generalizability in Radiomics Modeling: Possible Strategies in Radiologic and Statistical Perspectives

Radiomics, which involves the use of high-dimensional quantitative imaging features for predictive purposes, is a powerful tool for developing and testing medical hypotheses. Radiologic and statistical challenges in radiomics include those related to the reproducibility of imaging data, control of overfitting due to high dimensionality, and the generalizability of modeling. The aims of this review article are to clarify the distinctions between radiomics features and other omics and imaging data, to describe the challenges and potential strategies in reproducibility and feature selection, and to reveal the epidemiological background of modeling, thereby facilitating and promoting more reproducible and generalizable radiomics research.

Repeatability of Multiparametric Prostate MRI Radiomics Features

In this study we assessed the repeatability of radiomics features on small prostate tumors using test-retest Multiparametric Magnetic Resonance Imaging (mpMRI). The premise of radiomics is that quantitative image-based features can serve as biomarkers for detecting and characterizing disease. For such biomarkers to be useful, repeatability is a basic requirement, meaning its value must remain stable between two scans, if the conditions remain stable. We investigated repeatability of radiomics features under various preprocessing and extraction configurations including various image normalization schemes, different image pre-filtering, and different bin widths for image discretization. Although we found many radiomics features and preprocessing combinations with high repeatability (Intraclass Correlation Coefficient > 0.85), our results indicate that overall the repeatability is highly sensitive to the processing parameters. Neither image normalization, using a variety of approaches, nor the use of pre-filtering options resulted in consistent improvements in repeatability. We urge caution when interpreting radiomics features and advise paying close attention to the processing configuration details of reported results. Furthermore, we advocate reporting all processing details in radiomics studies and strongly recommend the use of open source implementations.

Use of Precision Imaging in the Evaluation of Pancreas Cancer

Pancreas cancer is an aggressive and fatal disease that will become one of the leading causes of cancer mortality by 2030. An all-out effort is underway to better understand the basic biologic mechanisms of this disease ranging from early development to metastatic disease. In order to change the course of this disease, diagnostic radiology imaging may play a vital role in providing a precise, noninvasive method for early diagnosis and assessment of treatment response. Recent progress in combining medical imaging, advanced image analysis and artificial intelligence, termed radiomics, can offer an innovate approach in detecting the earliest changes of tumor development as well as a rapid method for the detection of response. In this chapter, we introduce the principles of radiomics and demonstrate how it can provide additional information into tumor biology, early detection, and response assessments advancing the goals of precision imaging to deliver the right treatment to the right person at the right time.

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.

Application of Gray Level Co-Occurrence Matrix Analysis as a New Method for Enzyme Histochemistry Quantification

Enzyme histochemistry is a valuable histological method which provides a connection between morphology, activity, and spatial localization of investigated enzymes. Even though the method relies purely on arbitrary evaluations performed by the human eye, it is still wildly accepted and used in histo(patho)logy. Texture analysis emerged as an excellent tool for image quantification of subtle differences reflected in both spatial discrepancies and gray level values of pixels. The current study of texture analysis utilizes the gray-level co-occurrence matrix as a method for quantification of differences between ecto-5'-nucleotidase activities in healthy hippocampal tissue and tissue with marked neurodegeneration. We used the angular second moment, contrast (CON), correlation, inverse difference moment (INV), and entropy for texture analysis and receiver operating characteristic analysis with immunoblot and qualitative assessment of enzyme histochemistry as a validation. Our results strongly argue that co-occurrence matrix analysis could be used for the determination of fine differences in the enzyme activities with the possibility to ascribe those differences to regions or specific cell types. In addition, it emerged that INV and CON are especially useful parameters for this type of enzyme histochemistry analysis. We concluded that texture analysis is a reliable method for quantification of this descriptive technique, thus removing biases and adding it a quantitative dimension.

Radiomic features from PSMA PET for non-invasive intraprostatic tumor discrimination and characterization in patients with intermediate- and high-risk prostate cancer - a comparison study with histology reference

Purpose: To evaluate the performance of radiomic features (RF) derived from PSMA PET for intraprostatic tumor discrimination and non-invasive characterization of Gleason score (GS) and pelvic lymph node status. Patients and methods: Patients with prostate cancer (PCa) who underwent [⁶⁸Ga]-PSMA-11 PET/CT followed by radical prostatectomy and pelvic lymph node dissection were prospectively enrolled (n=20). Coregistered histopathological gross tumor volume (GTV-Histo) in the prostate served as reference. 133 RF were derived from GTV-Histo and from manually created segmentations of the intraprostatic tumor volume (GTV-Exp). Spearman´s correlation coefficients (ρ) were assessed between RF derived from the different GTVs. We additionally analyzed the differences in RF values for PCa and non-PCa tissues. Furthermore, areas under receiver-operating characteristics curves (AUC) were calculated and uni- and multivariate analyses were performed to evaluate the RF based discrimination of GS 7 and ≥8 disease and of patients with nodal spread (pN1) and non-nodal spread (pN0) in surgical specimen. The results found in the latter analyses were validated by a retrospective cohort of 40 patients. Results: Most RF from GTV-Exp showed strong correlations with RF from GTV-Histo (86% with ρ>0.7). 81% and 76% of RF from GTV-Exp and GTV-Histo significantly discriminated between PCa and non-PCa tissue. The texture feature QSZHGE discriminated between GS 7 and ≥8 considering GTV-Histo (AUC=0.93) and GTV-Exp (prospective cohort: AUC=0.91 / validation cohort: AUC=0.84). QSZHGE also discriminated between pN1 and pN0 disease considering GTV-Histo (AUC=0.85) and GTV-Exp (prospective cohort: AUC=0.87 / validation cohort: AUC=0.85). In uni- and multivariate analyses including patients of both cohorts QSZHGE was a statistically significant (p<0.01) predictor for PCa patients with GS ≥8 tumors and pN1 status. Conclusion: RF derived from PSMA PET discriminated between PCa and non-PCa tissue within the prostate. Additionally, the texture feature QSZHGE discriminated between GS 7 and GS ≥8 tumors and between patients with pN1 and pN0 disease. Our results support the role of RF in PSMA PET as a new tool for non-invasive PCa discrimination and characterization of its biological properties.

Potential Complementary Value of Noncontrast and Contrast Enhanced CT Radiomics in Colorectal Cancers

RATIONALE AND OBJECTIVES

The aim of our study was to assess the relationships between textural features extracted from contrast enhanced (CE) and noncontrast enhanced (NCE) computed tomography (CT) images of primary colorectal cancer, in order to identify radiomics features more likely to provide potential complementary information regarding outcome.

MATERIALS AND METHODS

Sixty-one patients with primary colorectal cancer underwent both CE-CT and NCE-CT scans within the same acquisition. First-order and textural features (with three different methods for grey-level discretization) were extracted from the tumor volume in both modalities and their correlation was assessed with Spearman's rank correlation (rs). Significance was assessed at p < 0.05 with correction for multiple comparisons. Kaplan-Meier estimation and log-rank tests were used to identify features associated with long term patient survival.

RESULTS

Moderate positive correlations were observed between CE-CT and NCE-CT histogram-derived entropy (Entropy_Hist) and area under the curve (CH_AUC) (rs = 0.49, p < 0.001 and rs= 0.45, p < 0.001, respectively). Some second and third order textural features were found highly correlated between CE-CT and NCE-CT, such as small zone-size emphasis SZSE (rs = 0.729, p < 0.001) and zone-size percentage (rs = 0.770, p < 0.001). Grey-levels discretization methods influenced these correlations. A few of the third order NCE-CT and CE-CT features were significantly associated with survival.

CONCLUSION

Some radiomics features with moderate correlations between nonenhanced and enhanced CT images were found to be associated with survival, thus suggesting that complementary prognostic value may be extracted from both modalities when available.

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.