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

Advancing the evaluation of radiation-induced vaginal toxicity using ultrasound radiomics: Phantom validation and pilot clinical study

BACKGROUND

Radiation-induced long-term toxicities, such as vaginal stenosis, severely impact the quality of life for patients undergoing pelvic radiotherapy (RT) for gynecologic (GYN) malignancies. However, current methods for assessing these toxicities rely on subjective physical examinations and patient-reported symptoms, leading to inconsistencies in grading and suboptimal management.

PURPOSE

This pilot study investigates the potential of ultrasound-based radiomics, specifically gray level co-occurrence matrix (GLCM) texture metrics, as objective and quantitative biomarkers for evaluating long-term radiation-induced vaginal toxicity.

METHODS

A two-phase study was conducted. First, a phantom study was performed to identify robust GLCM texture features with low variability [coefficient of variance (COV) < 10%] across ultrasound brightness settings. In a subsequent clinical pilot study, 22 female participants were recruited: 10 had received pelvic radiotherapy (RT) with follow-up times ranging from 8 to 23 months, while 12 served as non-RT controls. All participants underwent transvaginal ultrasound imaging, and GLCM texture features were extracted for analysis. A Mann-Whitney U test was used to assess between-group differences of distribution, with a p value < 0.05 identified as statistically significance. Cohen's d values were calculated to quantify effect sizes, with a value of greater than 0.8 indicating large effects.

RESULTS

Seventeen GLCM features demonstrated robustness (COVs < 10%) across brightness settings in the phantom study, including two with COVs < 1%, 10 with COVs between 1% and 5%, and five with COVs between 5% and 10%. In the clinical study, four texture features showed significant differences between the treated group and controls (p < 0.05). Specifically, the treated group exhibited a 15.5% increase in correlation (p = 0.03), a 35.8% decrease in contrast (p = 0.03), a 10.1% decrease in difference entropy (p = 0.04), and a 17.9% decrease in dissimilarity (p = 0.07).

CONCLUSION

This phantom and pilot study demonstrated that ultrasound GLCM features can serve as reliable quantitative biomarkers for assessing radiation-induced vaginal toxicity in female patients receiving pelvic RT for GYN cancers. Implementing these biomarkers in clinical practice could enhance the objectivity of toxicity evaluations, leading to more consistent grading and better-informed follow-up care for patients.

Practical use of radiomic features as a metric for image quality discrimination in [18F] FDG-PET: a pilot study

PURPOSE

Radiomics features have been utilised as group metrics of image quality in many areas of diagnostic radiology. In this pilot study, the relationship between technical metrics used in image quality assurance and visual grading scores provided by a radiologist were evaluated. Image dataset harmonisation allowed comparison between the two and allowed trends to be extracted. We propose a reproducible technique to identify the metrics.

METHODS

A retrospective chart review of 30 [18F] FDG-PET/CT performed in a nuclear medicine referral centre was performed. Image datasets were reprocessed to correspond to a bed duration of 180, 120, 60, 30 s per bed position and were analysed according to a pre-set bank of semi-quantitative features by a radiology resident. The extraction of radiomic features in PET images was performed using SLICER-RADIOMICS Module version 5.2.2. To facilitate the comparison of radiomic features and radiologist scoring data, normalisation was performed on both data sets. Fréchet distance analysis, Mean Square Error and Mean Absolute Error display the level of agreement between features and radiologist following the rescale of the data.

RESULTS

Of the 120 reprocessed image datasets, 115 were included in the study. We focused on overall image quality score rather than individual radiomic metrics as this identified the most robust trend. A significant difference in the 30 s image dataset with respect to each group individually and combined for the radiologist overall score was observed.

CONCLUSION

Our results show that a large percentage change in certain features can indicate a significant change in quality in clinically processed images.

Impact of harmonization and oversampling methods on radiomics analysis of multi-center imbalanced datasets: application to PET-based prediction of lung cancer subtypes

BACKGROUND

Medical imaging data frequently encounter image-generation heterogeneity and class imbalance properties, challenging strong generalized predictive performances with data-driven machine-learning methods. The purpose of this study was to investigate the impact of harmonization and oversampling methods on multi-center imbalanced datasets, with specific application to PET-based radiomics modeling for histologic subtype prediction in non-small cell lung cancer (NSCLC).

METHODS

The retrospective study included 245 patients with adenocarcinoma (ADC) and 78 patients with squamous cell carcinoma (SCC) from 4 centers. Utilizing 1502 radiomics features per patient, we trained, validated, and tested 4 machine-learning classifiers, to investigate the effect of no harmonization (NoH) or 4 feature harmonization methods, paired with no oversampling (NoO) or 5 oversampling methods on subtype prediction. Model performance was evaluated using the average area under the ROC curve (AUROC) and G-mean via 5 times 5-fold cross-validations. Statistical comparisons of the combined models against baseline (NoH + NoO) were performed for each fold of cross-validation using the DeLong test.

RESULTS

The number of cross-combinations with both AUROC and G-mean outperforming baseline in validation and testing was 15, 4, 2, and 7 (out of 29) for random forest (RF), linear discriminant analysis (LDA), logistic regression (LR), and support vector machine (SVM), respectively. ComBat harmonization combined with oversampling (SMOTE) via RF yielded better performance than baseline (AUROC and G-mean of validation: 0.725 vs. 0.608 and 0.625 vs. 0.398; testing: 0.637 vs. 0.567 and 0.506 vs. 0.287), though statistical significances were not observed.

CONCLUSIONS

Applying harmonization and oversampling methods in multi-center imbalanced datasets can improve NSCLC-subtype prediction, but the effect varies widely across classifiers. We have created open-source comparisons of harmonization and oversampling on different classifiers for comprehensive evaluations in different studies.

A deep learning method for the recovery of standard-dose imaging quality from ultra-low-dose PET on wavelet domain

PURPOSE

Recent development in positron emission tomography (PET) dramatically increased the effective sensitivity by increasing the geometric coverage leading to total-body PET imaging. This encouraging breakthrough brings the hope of ultra-low dose PET imaging equivalent to transatlantic flight with the assistance of deep learning (DL)-based methods. However, conventional DL approaches face limitations in addressing the heterogeneous domain of PET imaging. This study aims to develop a wavelet-based DL method capable of restoring high-quality imaging from ultra-low-dose PET scans.

MATERIALS AND METHODS

In contrast to conventional DL techniques that denoise images in the spatial domain, we introduce WaveNet, a novel approach that inputs wavelet-decomposed frequency components of PET imaging to perform denoising in the frequency domain. A dataset comprising total-body 18F -FDG PET images of 1447, acquired using total-body PET scanners including Biograph Vision Quadra (Siemens Healthineers) and uEXPLORER (United Imaging) in Bern and Shanghai, was utilized for developing and testing the proposed method. The quality of enhanced images was assessed using a customized scoring system, which incorporated weighted global physical metrics and local indices.

RESULTS

Our proposed WaveNet consistently outperforms the baseline UNet model across all levels of dose reduction factors (DRF), with greater improvements observed as image quality decreases. Statistical analysis (p < 0.05) and visual inspection validated the superiority of WaveNet. Moreover, WaveNet demonstrated superior generalizability when applied to two cross-scanner datasets (p < 0.05).

CONCLUSION

WaveNet developed with total-body PET scanners may offer a computational-friendly and robust approach to recover image quality from ultra-low-dose PET imaging. Its adoption may enhance the reliability and clinical acceptance of DL-based dose reduction techniques.

Robustness of 18F-FDG PET Radiomic Features in Lung Cancer: Impact of Advanced Reconstruction Algorithm

18F-FDG PET radiomics is emerging as a promising tool to identify imaging biomarkers for quantifying intratumor heterogeneity in lung cancer. However, the robustness of PET radiomic features (RFs) is influenced by factors such as image reconstruction algorithms, tumor segmentation, and discretization. Although the impact of these factors on RFs has been explored, the specific influence of the advanced block sequential regularized expectation maximization (BSREM) reconstruction algorithm remains unclear. This study investigated the potential variations in PET RFs associated with different factors when using BSREM. Methods: Retrospective data of 18F-FDG PET from 120 lung cancer patients were reconstructed twice using advanced BSREM and conventional ordered-subset expectation maximization methods. For each reconstruction set, 3 tumor segmentations were performed, including manual, 40% threshold, and Nestle methods. Two discretization methods using absolute and relative settings were applied for each dataset before RF extraction. Stable and robust RFs were assessed by the coefficient of variance and intraclass correlation coefficient, respectively. Results: High instability was exhibited by 19%, 33%, and 36% of RFs, with a coefficient of variation of more than 20% for reconstruction, segmentation, and discretization, respectively. Conversely, 60%, 19%, and 35% of RFs demonstrated robustness against these factors, with an intraclass correlation coefficient of more than 0.90. The comparative evaluation revealed significantly greater robustness for most RF subtypes in BSREM than in ordered-subset expectation maximization under varying segmentation and discretization conditions (P < 0.05). Conclusion: The stability and robustness of PET RFs are enhanced if BSREM is applied rather than the conventional method. Study results suggest that the advanced reconstruction method could offer potential benefits in providing consistent PET-based radiomic analysis for improving diagnostic and prognostic value.

Robust and generalizable artificial intelligence for multi-organ segmentation in ultra-low-dose total-body PET imaging: a multi-center and cross-tracer study

PURPOSE

Positron Emission Tomography (PET) is a powerful molecular imaging tool that visualizes radiotracer distribution to reveal physiological processes. Recent advances in total-body PET have enabled low-dose, CT-free imaging; however, accurate organ segmentation using PET-only data remains challenging. This study develops and validates a deep learning model for multi-organ PET segmentation across varied imaging conditions and tracers, addressing critical needs for fully PET-based quantitative analysis.

MATERIALS AND METHODS

This retrospective study employed a 3D deep learning-based model for automated multi-organ segmentation on PET images acquired under diverse conditions, including low-dose and non-attenuation-corrected scans. Using a dataset of 798 patients from multiple centers with varied tracers, model robustness and generalizability were evaluated via multi-center and cross-tracer tests. Ground-truth labels for 23 organs were generated from CT images, and segmentation accuracy was assessed using the Dice similarity coefficient (DSC).

RESULTS

In the multi-center dataset from four different institutions, our model achieved average DSC values of 0.834, 0.825, 0.819, and 0.816 across varying dose reduction factors and correction conditions for FDG PET images. In the cross-tracer dataset, the model reached average DSC values of 0.737, 0.573, 0.830, 0.661, and 0.708 for DOTATATE, FAPI, FDG, Grazytracer, and PSMA, respectively.

CONCLUSION

The proposed model demonstrated effective, fully PET-based multi-organ segmentation across a range of imaging conditions, centers, and tracers, achieving high robustness and generalizability. These findings underscore the model's potential to enhance clinical diagnostic workflows by supporting ultra-low dose PET imaging.

CLINICAL TRIAL NUMBER

Not applicable. This is a retrospective study based on collected data, which has been approved by the Research Ethics Committee of Ruijin Hospital affiliated to Shanghai Jiao Tong University School of Medicine.

Myocardial perfusion SPECT radiomic features reproducibility assessment: Impact of image reconstruction and harmonization

BACKGROUND

Coronary artery disease (CAD) has one of the highest mortality rates in humans worldwide. Single-photon emission computed tomography (SPECT) myocardial perfusion imaging (MPI) provides clinicians with myocardial metabolic information non-invasively. However, there are some limitations to interpreting SPECT images performed by physicians or automatic quantitative approaches. Radiomics analyzes images objectively by extracting quantitative features and can potentially reveal biological characteristics that the human eye cannot detect. However, the reproducibility and repeatability of some radiomic features can be highly susceptible to segmentation and imaging conditions.

PURPOSE

We aimed to assess the reproducibility of radiomic features extracted from uncorrected MPI-SPECT images reconstructed with 15 different settings before and after ComBat harmonization, along with evaluating the effectiveness of ComBat in realigning feature distributions.

MATERIALS AND METHODS

A total of 200 patients (50% normal and 50% abnormal) including rest and stress (without attenuation and scatter corrections) MPI-SPECT images were included. Images were reconstructed using 15 combinations of filter cut-off frequencies, filter orders, filter types, reconstruction algorithms, number of iterations and subsets resulting in 6000 images. Image segmentation was performed on the left ventricle in the first reconstruction for each patient and applied to 14 others. A total of 93 radiomic features were extracted from the segmented area, and ComBat was used to harmonize them. The intraclass correlation coefficient (ICC) and overall concordance correlation coefficient (OCCC) tests were performed before and after ComBat to examine the impact of each parameter on feature robustness and to assess harmonization efficiency. The ANOVA and the Kruskal-Wallis tests were performed to evaluate the effectiveness of ComBat in correcting feature distributions. In addition, the Student's t-test, Wilcoxon rank-sum, and signed-rank tests were implemented to assess the significance level of the impacts made by each parameter of different batches and patient groups (normal vs. abnormal) on radiomic features.

RESULTS

Before applying ComBat, the majority of features (ICC: 82, OCCC: 61) achieved high reproducibility (ICC/OCCC ≥ 0.900) under every batch except Reconstruction. The largest and smallest number of poor features (ICC/OCCC < 0.500) were obtained by IterationSubset and Order batches, respectively. The most reliable features were from the first-order (FO) and gray-level co-occurrence matrix (GLCM) families. Following harmonization, the minimum number of robust features increased (ICC: 84, OCCC: 78). Applying ComBat showed that Order and Reconstruction were the least and the most responsive batches, respectively. The most robust families, in a descending order, were found to be FO, neighborhood gray-tone difference matrix (NGTDM), GLCM, gray-level run length matrix (GLRLM), gray-level size zone matrix (GLSZM), and gray-level dependence matrix (GLDM) under Cut-off, Filter, and Order batches. The Wilcoxon rank-sum test showed that the number of robust features significantly differed under most batches in the Normal and Abnormal groups.

CONCLUSION

The majority of radiomic features show high levels of robustness across different OSEM reconstruction parameters in uncorrected MPI-SPECT. ComBat is effective in realigning feature distributions and enhancing radiomic features reproducibility.

[Investigation of the Influence of Image Reconstruction Parameters to Improve the Ability to Depict Internal Tumor Necrosis]

PURPOSE

There are attempts to assess tumor heterogeneity by texture analysis. However, the ordered subsets-expectation maximization (OSEM) reconstruction method has problems depicting heterogeneities. The aim of this study was to identify image reconstruction parameters that improve the ability to depict internal tumor necrosis using a self-made phantom that simulates internal necrosis.

METHODS

Self-made phantoms were prepared using polypropylene cylinders with inner diameters of 18.0 mm and 6.0 mm. The concentration ratios of the simulated tumor : tumor interior were 4 : 0 and 4 : 1. For each reconstruction method, the iteration for OSEM and OSEM+point spread function (PSF) were 1 to 25 and the subset was 12. The β values for block sequential regularized expectation maximization (BSREM) were set between 10 and 400. We evaluated the features of the profile curve, contrast-to-noise ratio, and grey-level co-occurrence matrix (GLCM).

RESULTS

In the phantom study, OSEM and OSEM+PSF showed a better delineation of the differences between the inside and outside of the cylinder as iteration was increased and BSREM showed a better delineation as β was decreased. The highest value for each feature, both 4 : 0 and 4 : 1, was BSREM β 10 for angular second moment (ASM) and inverse differential moment (IDM), OSEM iteration 25 for contrast and entropy.

CONCLUSION

We have identified image reconstruction parameters that improve the ability to visualize internal tumor necrosis. The parameter was BRSEM β 10.

Mitigating Interobserver Variability in Radiomics with ComBat: A Feasibility Study

This study investigates Intraobserver Features Variability (IFV) in radiomics studies and assesses the effectiveness of the ComBat harmonization method in mitigating these effects. Methods: This study utilizes data from the NSCLC-Radiomics-Interobserver1 dataset, comprising CT scans of 22 Non-Small Cell Lung Cancer (NSCLC) patients, with multiple Gross Tumor Volume (GTV) delineations performed by five radiation oncologists. Segmentation was completed manually ("vis") or by autosegmentation with manual editing ("auto"). A total of 1229 radiomic features were extracted from each GTV, segmentation method, and oncologist. Features extracted included first order, shape, GLCM, GLRLM, GLSZM, and GLDM from original, wavelet-filtered, and LoG-filtered images. Results: Before implementing ComBat harmonization, 83% of features exhibited p-values below 0.05 in the "vis" approach; this percentage decreased to 34% post-harmonization. Similarly, for the "auto" approach, 75% of features demonstrated statistical significance prior to ComBat, but this figure declined to 33% after its application. Among a subset of three expert radiation oncologists, percentages changed from 77% to 25% for "vis" contouring and from 64% to 23% for "auto" contouring. This study demonstrates that ComBat harmonization could effectively reduce IFV, enhancing the feasibility of multicenter radiomics studies. It also highlights the significant impact of physician experience on radiomics analysis outcomes.

Multicenter PET image harmonization using generative adversarial networks

PURPOSE

To improve reproducibility and predictive performance of PET radiomic features in multicentric studies by cycle-consistent generative adversarial network (GAN) harmonization approaches.

METHODS

GAN-harmonization was developed to harmonize whole-body PET scans to perform image style and texture translation between different centers and scanners. GAN-harmonization was evaluated by application to two retrospectively collected open datasets and different tasks. First, GAN-harmonization was performed on a dual-center lung cancer cohort (127 female, 138 male) where the reproducibility of radiomic features in healthy liver tissue was evaluated. Second, GAN-harmonization was applied to a head and neck cancer cohort (43 female, 154 male) acquired from three centers. Here, the clinical impact of GAN-harmonization was analyzed by predicting the development of distant metastases using a logistic regression model incorporating first-order statistics and texture features from baseline 18F-FDG PET before and after harmonization.

RESULTS

Image quality remained high (structural similarity: left kidney ≥ 0.800, right kidney ≥ 0.806, liver ≥ 0.780, lung ≥ 0.838, spleen ≥ 0.793, whole-body ≥ 0.832) after image harmonization across all utilized datasets. Using GAN-harmonization, inter-site reproducibility of radiomic features in healthy liver tissue increased at least by ≥ 5 ± 14% (first-order), ≥ 16 ± 7% (GLCM), ≥ 19 ± 5% (GLRLM), ≥ 16 ± 8% (GLSZM), ≥ 17 ± 6% (GLDM), and ≥ 23 ± 14% (NGTDM). In the head and neck cancer cohort, the outcome prediction improved from AUC 0.68 (95% CI 0.66-0.71) to AUC 0.73 (0.71-0.75) by application of GAN-harmonization.

CONCLUSIONS

GANs are capable of performing image harmonization and increase reproducibility and predictive performance of radiomic features derived from different centers and scanners.

An 18F-FDG PET/CT and Mean Lung Dose Model to Predict Early Radiation Pneumonitis in Stage III Non-Small Cell Lung Cancer Patients Treated with Chemoradiation and Immunotherapy

Radiation pneumonitis (RP) that develops early (i.e., within 3 mo) (RPEarly) after completion of concurrent chemoradiation (cCRT) leads to treatment discontinuation and poorer survival for patients with stage III non-small cell lung cancer. Since no RPEarly risk model exists, we explored whether published RP models and pretreatment 18F-FDG PET/CT-derived features predict RPEarly Methods: One hundred sixty patients with stage III non-small cell lung cancer treated with cCRT and consolidative immunotherapy were analyzed for RPEarly Three published RP models that included the mean lung dose (MLD) and patient characteristics were examined. Pretreatment 18F-FDG PET/CT normal-lung SUV featured included the following: 10th percentile of SUV (SUVP10), 90th percentile of SUV (SUVP90), SUVmax, SUVmean, minimum SUV, and SD. Associations between models/features and RPEarly were assessed using area under the receiver-operating characteristic curve (AUC), P values, and the Hosmer-Lemeshow test (pHL). The cohort was randomly split, with similar RPEarly rates, into a 70%/30% derivation/internal validation subset. Results: Twenty (13%) patients developed RPEarly Predictors for RPEarly were MLD alone (AUC, 0.72; P = 0.02; pHL, 0.87), SUVP10, SUVP90, and SUVmean (AUC, 0.70-0.74; P = 0.003-0.006; pHL, 0.67-0.70). The combined MLD and SUVP90 model generalized in the validation subset and was deemed the final RPEarly model (RPEarly risk = 1/[1+e(- x )]; x = -6.08 + [0.17 × MLD] + [1.63 × SUVP90]). The final model refitted in the 160 patients indicated improvement over the published MLD-alone model (AUC, 0.77 vs. 0.72; P = 0.0001 vs. 0.02; pHL, 0.65 vs. 0.87). Conclusion: Patients at risk for RPEarly can be detected with high certainty by combining the normal lung's MLD and pretreatment 18F-FDG PET/CT SUVP90 This refined model can be used to identify patients at an elevated risk for premature immunotherapy discontinuation due to RPEarly and could allow for interventions to improve treatment outcomes.

Impact of different reconstruction algorithms and setting parameters on radiomics features of PSMA PET images: A preliminary study

PURPOSE

Radiomics analysis of oncologic positron emission tomography (PET) images is an area of significant activity and potential. The reproducibility of radiomics features is an important consideration for routine clinical use. This preliminary study investigates the robustness of radiomics features in PSMA-PET images across penalized-likelihood (Q.Clear) and standard ordered subset expectation maximization (OSEM) reconstruction algorithms and their setting parameters in phantom and prostate cancer (PCa) patients.

METHOD

A NEMA image quality (IQ) phantom and 8 PCa patients were selected for phantom and patient analyses, respectively. PET images were reconstructed using Q.Clear (reconstruction β-value: 100-700, at intervals of 100 for both NEMA IQ phantom and patients) and OSEM (duration: 15sec, 30sec, 1 min, 2 min, 3 min, 4 min and 5 min for NEMA phantom and duration: 30 s, 1 min and 2 min for patients) reconstruction methods. Subsequently, 129 radiomic features were extracted from the reconstructed images. The coefficient of variation (COV) of each feature across reconstruction methods and their parameters was calculated to determine feature robustness.

RESULTS

The extracted radiomics features showed a different range of variability, depending on the reconstruction algorithms and setting parameters. Specifically, 23.0 % and 53.5 % of features were found as robust against β-value variations in Q.Clear and different durations in OSEM reconstruction algorithms, respectively. Taking into account the two algorithms and their parameters, eleven features (8.5 %) showed COV ≤ 5 % and eighteen (14 %) showed 5 % 20 %. The mean COVs of the extracted radiomics features were significantly different between the two reconstruction methods (p < 0.05) except for the phantom morphological features.

CONCLUSIONS

All radiomics features were affected by reconstruction methods and parameters, but features with small or very small variations are considered better candidates for reproducible quantification of either tumor or metastatic tissues in clinical trials. There is a need for standardization before the implementation of PET radiomics in clinical practice.

Robustness of [18F]FDG PET/CT radiomic analysis in the setting of drug-induced cardiotoxicity

BACKGROUND AND OBJECTIVES

Standardization of radiomic data acquisition protocols is still at a very early stage, revealing a strong need to work towards the definition of uniform image processing methodologies The aim of this study is to identify sources of variability in radiomic data derived from image discretization and resampling methodologies prior to image feature extraction. Furthermore, to identify robust potential image-based biomarkers for the early detection of cardiotoxicity.

METHODS

Image post-acquisition processing, interpolation, and volume of interest (VOI) segmentation were performed. Four experiments were conducted to assess the reliability in terms of the intraclass correlation coefficient (ICC) of the radiomic features and the effects of the variation of voxel size and gray level discretization. Statistical analysis was performed separating the patients according to cardiotoxicity diagnosis. Differences of texture features were studied with Mann-Whitney U test. P-values <0.05 after multiple testing correction were considered statistically significant. Additionally, a non-supervised k-Means clustering algorithm was evaluated.

RESULTS

The effect of the variation in the voxel size demonstrated a non-dependency relationship with the values of the radiomic features, regardless of the chosen discretization method. The median ICC values were 0.306 and 0.872 for absolute agreement and consistency, respectively, when varying the discretization bin number. The median ICC values were 0.678 and 0.878 for absolute agreement and consistency, respectively, when varying the discretization bin size. A total of 16 first order, 6 Gray Level Co-occurrence Matrix (GLCM), 4 Gray Level Dependence Matrix (GLDM) and 4 Gray Level Run Length Matrix (GLRLM) features demonstrated statistically significant differences between the diagnosis groups for interim scans (P<0.05) for the fixed bin size (FBS) discretization methodology. However, no statistically significant differences between diagnostic groups were found for the fixed bin number (FBN) discretization methodology. Two clusters based on the radiomic features were identified.

CONCLUSIONS

Gray level discretization has a major impact on the repeatability of the radiomic features. The selection of the optimal processing methodology has led to the identification of texture-based patterns for the differentiation of early cardiac damage profiles.

Exploring published and novel pre-treatment CT and PET radiomics to stratify risk of progression among early-stage non-small cell lung cancer patients treated with stereotactic radiation

PURPOSE

Disease progression after definitive stereotactic body radiation therapy (SBRT) for early-stage non-small cell lung cancer (NSCLC) occurs in 20-40% of patients. Here, we explored published and novel pre-treatment CT and PET radiomics features to identify patients at risk of progression.

MATERIALS/METHODS

Published CT and PET features were identified and explored along with 15 other CT and PET features in 408 consecutively treated early-stage NSCLC patients having CT and PET < 3 months pre-SBRT (training/set-aside validation subsets: n = 286/122). Features were associated with progression-free survival (PFS) using bootstrapped Cox regression (Bonferroni-corrected univariate predictor: p ≤ 0.002) and only non-strongly correlated predictors were retained (|Rs|<0.70) in forward-stepwise multivariate analysis.

RESULTS

Tumor diameter and SUVmax were the two most frequently reported features associated with progression/survival (in 6/20 and 10/20 identified studies). These two features and 12 of the 15 additional features (CT: 6; PET: 6) were candidate PFS predictors. A re-fitted model including diameter and SUVmax presented with the best performance (c-index: 0.78; log-rank p-value < 0.0001). A model built with the two best additional features (CTspiculation1 and SUVentropy) had a c-index of 0.75 (log-rank p-value < 0.0001).

CONCLUSIONS

A re-fitted pre-treatment model using the two most frequently published features - tumor diameter and SUVmax - successfully stratified early-stage NSCLC patients by PFS after receiving SBRT.

Machine Learning Radiomics Signature for Differentiating Lymphoma versus Benign Splenomegaly on CT

BACKGROUND

We aimed to develop and validate a preoperative CT-based radiomics signature for differentiating lymphoma versus benign splenomegaly.

METHODS

We retrospectively analyzed CT studies from 139 patients (age range 26-93 years, 43% female) between 2011 and 2019 with histopathological diagnosis of the spleen (19 lymphoma, 120 benign) and divided them into developing (n = 79) and testing (n = 60) datasets. The volumetric radiomic features were extracted from manual segmentation of the whole spleen on venous-phase CT imaging using PyRadiomics package. LASSO regression was applied for feature selection and development of the radiomic signature, which was interrogated with the complete blood cell count and differential count. All p values < 0.05 were considered to be significant.

RESULTS

Seven features were selected for constructing the radiomic signature after feature selection, including first-order statistics (10th percentile and Robust Mean Absolute Deviation), shape-based (Surface Area), and texture features (Correlation, MCC, Small Area Low Gray-level Emphasis and Low Gray-level Zone Emphasis). The radiomic signature achieved an excellent diagnostic accuracy of 97%, sensitivity of 89%, and specificity of 98%, distinguishing lymphoma versus benign splenomegaly in the testing dataset. The radiomic signature significantly correlated with the platelet and segmented neutrophil percentage.

CONCLUSIONS

CT-based radiomics signature can be useful in distinguishing lymphoma versus benign splenomegaly and can reflect the changes in underlying blood profiles.

Enhancing Interoperability and Harmonisation of Nuclear Medicine Image Data and Associated Clinical Data

Nuclear imaging techniques such as positron emission tomography (PET) and single photon emission computed tomography (SPECT) in combination with computed tomography (CT) are established imaging modalities in clinical practice, particularly for oncological problems. Due to a multitude of manufacturers, different measurement protocols, local demographic or clinical workflow variations as well as various available reconstruction and analysis software, very heterogeneous datasets are generated. This review article examines the current state of interoperability and harmonisation of image data and related clinical data in the field of nuclear medicine. Various approaches and standards to improve data compatibility and integration are discussed. These include, for example, structured clinical history, standardisation of image acquisition and reconstruction as well as standardised preparation of image data for evaluation. Approaches to improve data acquisition, storage and analysis will be presented. Furthermore, approaches are presented to prepare the datasets in such a way that they become usable for projects applying artificial intelligence (AI) (machine learning, deep learning, etc.). This review article concludes with an outlook on future developments and trends related to AI in nuclear medicine, including a brief research of commercial solutions.

Optimal batch determination for improved harmonization and prognostication of multi-center PET/CT radiomics feature in head and neck cancer

Objective. To determine the optimal approach for identifying and mitigating batch effects in PET/CT radiomics features, and further improve the prognosis of patients with head and neck cancer (HNC), this study investigated the performance of three batch harmonization methods.Approach. Unsupervised harmonization identified the batch labels by K-means clustering. Supervised harmonization regarding the image acquisition factors (center, manufacturer, scanner, filter kernel) as known/given batch labels, and Combat harmonization was then implemented separately and sequentially based on the batch labels, i.e. harmonizing features among batches determined by each factor individually or harmonizing features among batches determined by multiple factors successively. Extensive experiments were conducted to predict overall survival (OS) on public PET/CT datasets that contain 800 patients from 9 centers.Main results. In the external validation cohort, results show that compared to original models without harmonization, Combat harmonization would be beneficial in OS prediction with C-index of 0.687-0.740 versus 0.684-0.767. Supervised harmonization slightly outperformed unsupervised harmonization in all models (C-index: 0.692-0.767 versus 0.684-0.750). Separate harmonization outperformed sequential harmonization in CT_m+clinic and CT_cm+clinic models with C-index of 0.752 and 0.722, respectively, while sequential harmonization involved clinical features in PET_rs+clinic model further improving the performance and achieving the highest C-index of 0.767.Significance. Optimal batch determination especially sequential harmonization for Combat holds the potential to improve the prognostic power of radiomics model in multi-center HNC dataset with PET/CT imaging.

Comparative analysis of batch correction methods for FDG PET/CT using metabolic radiogenomic data of lung cancer patients

In radiomics research, the issue of different instruments being used is significant. In this study, we compared three correction methods to reduce the batch effects in radiogenomic data from fluorodeoxyglucose (FDG) PET/CT images of lung cancer patients. Texture features of the FDG PET/CT images and genomic data were retrospectively obtained. The features were corrected with different methods: phantom correction, ComBat method, and Limma method. Batch effects were estimated using three analytic tools: principal component analysis (PCA), the k-nearest neighbor batch effect test (kBET), and the silhouette score. Finally, the associations of features and gene mutations were compared between each correction method. Although the kBET rejection rate and silhouette score were lower in the phantom-corrected data than in the uncorrected data, a PCA plot showed a similar variance. ComBat and Limma methods provided correction with low batch effects, and there was no significant difference in the results of the two methods. In ComBat- and Limma-corrected data, more texture features exhibited a significant association with the TP53 mutation than in those in the phantom-corrected data. This study suggests that correction with ComBat or Limma methods can be more effective or equally as effective as the phantom method in reducing batch effects.

Fusion-based tensor radiomics using reproducible features: Application to survival prediction in head and neck cancer

BACKGROUND

Numerous features are commonly generated in radiomics applications as applied to medical imaging, and identification of robust radiomics features (RFs) can be an important step to derivation of reliable, reproducible solutions. In this work, we utilize a tensor radiomics (TR) framework, where numerous fusions are explored, to generate different flavours of RFs, and we aimed to identify RFs that are robust to fusion techniques in head and neck cancer. Overall, we aimed to predict progression-free survival (PFS) using Hybrid Machine Learning Systems (HMLS) and reproducible RFs.

METHODS

The study was performed on 408 patients with head and neck cancer from The Cancer Imaging Archive. After image preprocessing, 15 fusion techniques were employed to combine Positron Emission Tomography (PET) and Computed Tomography (CT) images. Subsequently, 215 RFs were extracted through a standardized radiomics software, with 17 'flavours' generated using PET-only, CT-only, and 15 fused PET&CT images. The variability of RFs across flavours was studied using the Intraclass Correlation Coefficient (ICC). Furthermore, the features were categorized into seven reliability groups, 106 reproducible RFs with ICC>0.75 were selected, highly correlated flavours were removed, Principal Component Analysis was used to convert 17 flavours to 1 attribute, the polynomial function was utilized to increase RFs, and Analysis of variance (ANOVA) was used to select the relevant attributes. Finally, 3 classifiers including Random Forest (RFC), Logistic regression (LR), and Multi-layer perceptron were applied to the preselected relevant attributes to predict binary PFS. In 5-fold cross-validation, 80% of 4 divisions were utilized to train the model, and the remaining 20% was utilized to evaluate the model. Further, the remaining fold was used for external nested testing.

RESULTS

Reliability analysis indicated that most morphological features belong to the high-reliability category. By contrast, local intensity and statistical features extracted from images belong to the low-reliability category. In the tensor framework, the highest 5-fold cross-validation accuracy of 76.7%±3.3% with an external nested testing of 70.6%±6.7% resulted from the reproducible TR+polynomial function+ANOVA+LR algorithm while the accuracy of 70.0%±4.2% with the external nested testing of 67.7%±4.9% was achieved through the PCA fusion+RFC (non-tensor paradigm).

CONCLUSIONS

This study demonstrated that using reproducible RFs as utilized within a tensor fusion radiomics framework, linked with ANOVA and LR, added value to prediction of progression-free survival outcome in head and neck cancer patients.

Development of a radiomic-clinical nomogram for prediction of survival in patients with diffuse large B-cell lymphoma treated with chimeric antigen receptor T cells

BACKGROUND

In our current work, an 18F-FDG PET/CT radiomics-based model was developed to assess the progression-free survival (PFS) and overall survival (OS) of patients with relapsed or refractory (R/R) diffuse large B-cell lymphoma (DLBCL) who received chimeric antigen receptor (CAR)-T cell therapy.

METHODS

A total of 61 DLBCL cases receiving 18F-FDG PET/CT before CAR-T cell infusion were included in the current analysis, and these patients were randomly assigned to a training cohort (n = 42) and a validation cohort (n = 19). Radiomic features from PET and CT images were obtained using LIFEx software, and radiomics signatures (R-signatures) were then constructed by choosing the optimal parameters according to their PFS and OS. Subsequently, the radiomics model and clinical model were constructed and validated.

RESULTS

The radiomics model that integrated R-signatures and clinical risk factors showed superior prognostic performance compared with the clinical models in terms of both PFS (C-index: 0.710 vs. 0.716; AUC: 0.776 vs. 0.712) and OS (C-index: 0.780 vs. 0.762; AUC: 0.828 vs. 0.728). For validation, the C-index of the two approaches was 0.640 vs. 0.619 and 0.676 vs. 0.699 for predicting PFS and OS, respectively. Moreover, the AUC was 0.886 vs. 0.635 and 0.778 vs. 0.705, respectively. The calibration curves indicated good agreement, and the decision curve analysis suggested that the net benefit of radiomics models was higher than that of clinical models.

CONCLUSIONS

PET/CT-derived R-signature could be a potential prognostic biomarker for R/R DLBCL patients undergoing CAR-T cell therapy. Moreover, the risk stratification could be further enhanced when the PET/CT-derived R-signature was combined with clinical factors.

Impact of Tracer Dose Reduction in [18 F]-Labelled Fluorodeoxyglucose-Positron Emission Tomography ([18 F]-FDG)-PET) on Texture Features and Histogram Indices: A Study in Homogeneous Tissues of Phantom and Patient

BACKGROUND

Histogram indices (HIs) and texture features (TFs) are considered to play an important role in future oncologic PET-imaging and it is unknown how these indices are affected by changes of tracer doses. A randomized undersampling of PET list mode data enables a simulation of tracer dose reduction. We performed a phantom study to compare HIs/TFs of simulated and measured tracer dose reductions and evaluated changes of HIs/TFs in the liver of patients with PETs from simulated reduced tracer doses. Overall, 42 HIs/TFs were evaluated in a NEMA phantom at measured and simulated doses (stepwise reduction of [18 F] from 100% to 25% of the measured dose). [18 F]-FDG-PET datasets of 15 patients were simulated from 3.0 down to 0.5 MBq/kgBW in intervals of 0.25 MBq/kgBW. HIs/TFs were calculated from two VOIs placed in physiological tissue of the right and left liver lobe and linear correlations and coefficients of variation analysis were performed.

RESULTS

All 42 TFs did not differ significantly in measured and simulated doses (p > 0.05). Also, 40 TFs showed the same behaviour over dose reduction regarding differences in the same group (measured or simulated), and for 26 TFs a linear behaviour over dose reduction for measured and simulated doses could be validated. Out of these, 13 TFs could be identified, which showed a linear change in TF value in both the NEMA phantom and patient data and therefore should maintain the same informative value when transferred in a dose reduction setting. Out of this Homogeneity 2, Entropy and Zone size non-uniformity are of special interest because they have been described as preferentially considerable for tumour heterogeneity characterization.

CONCLUSIONS

We could show that there was no significant difference of measured and simulated HIs/TFs in the phantom study and most TFs reveal a linear behaviour over dose reduction, when tested in homogeneous tissue. This indicates that texture analysis in PET might be robust to dose modulations.

Development and External Validation of 18F-FDG PET-Based Radiomic Model for Predicting Pathologic Complete Response after Neoadjuvant Chemotherapy in Breast Cancer

The aim of our retrospective study is to develop and externally validate an 18F-FDG PET-derived radiomics model for predicting pathologic complete response (pCR) after neoadjuvant chemotherapy (NAC) in breast cancer patients. A total of 87 breast cancer patients underwent curative surgery after NAC at Soonchunhyang University Seoul Hospital and were randomly assigned to a training cohort and an internal validation cohort. Radiomic features were extracted from pretreatment PET images. A radiomic-score model was generated using the LASSO method. A combination model incorporating significant clinical variables was constructed. These models were externally validated in a separate cohort of 28 patients from Soonchunhyang University Buscheon Hospital. The model performances were assessed using area under the receiver operating characteristic (AUC). Seven radiomic features were selected to calculate the radiomic-score. Among clinical variables, human epidermal growth factor receptor 2 status was an independent predictor of pCR. The radiomic-score model achieved good discriminability, with AUCs of 0.963, 0.731, and 0.729 for the training, internal validation, and external validation cohorts, respectively. The combination model showed improved predictive performance compared to the radiomic-score model alone, with AUCs of 0.993, 0.772, and 0.906 in three cohorts, respectively. The 18F-FDG PET-derived radiomic-based model is useful for predicting pCR after NAC in breast cancer.

Radiomic signatures based on pretreatment 18F-FDG PET/CT, combined with clinicopathological characteristics, as early prognostic biomarkers among patients with invasive breast cancer

Purpose

The aim of this study was to investigate the predictive role of fluorine-18 fluorodeoxyglucose positron emission tomography/computed tomography (18F-FDG PET/CT) in the prognostic risk stratification of patients with invasive breast cancer (IBC). To achieve this, we developed a clinicopathologic-radiomic-based model (C-R model) and established a nomogram that could be utilized in clinical practice.

Methods

We retrospectively enrolled a total of 91 patients who underwent preoperative 18F-FDG PET/CT and randomly divided them into training (n=63) and testing cohorts (n=28). Radiomic signatures (RSs) were identified using the least absolute shrinkage and selection operator (LASSO) regression algorithm and used to compute the radiomic score (Rad-score). Patients were assigned to high- and low-risk groups based on the optimal cut-off value of the receiver operating characteristic (ROC) curve analysis for both Rad-score and clinicopathological risk factors. Univariate and multivariate Cox regression analyses were performed to determine the association between these variables and progression-free survival (PFS) or overall survival (OS). We then plotted a nomogram integrating all these factors to validate the predictive performance of survival status.

Results

The Rad-score, age, clinical M stage, and minimum standardized uptake value (SUVmin) were identified as independent prognostic factors for predicting PFS, while only Rad-score, age, and clinical M stage were found to be prognostic factors for OS in the training cohort. In the testing cohort, the C-R model showed superior performance compared to single clinical or radiomic models. The concordance index (C-index) values for the C-R model, clinical model, and radiomic model were 0.816, 0.772, and 0.647 for predicting PFS, and 0.882, 0.824, and 0.754 for OS, respectively. Furthermore, decision curve analysis (DCA) and calibration curves demonstrated that the C-R model had a good ability for both clinical net benefit and application.

Conclusion

The combination of clinicopathological risks and baseline PET/CT-derived Rad-score could be used to evaluate the prognosis in patients with IBC. The predictive nomogram based on the C-R model further enhanced individualized estimation and allowed for more accurate prediction of patient outcomes.

Pretherapy 18F-FDG PET/CT in Predicting Disease Relapse in Patients With Immunoglobulin G4-Related Disease: A Prospective Study

OBJECTIVE

To investigate whether the levels of inflammation detected by ¹⁸F-fluorodeoxyglucose (FDG) positron emission tomography (PET)/computed tomography (CT) can predict disease relapse in immunoglobulin G4-related disease (IgG4-RD) patients receiving standard induction steroid therapy.

MATERIALS AND METHODS

This prospective study analyzed pretherapy FDG PET/CT images from 48 patients (mean age, 63 ± 12.9 years; 45 males and 3 females) diagnosed with IgG4-RD between September 2008 and February 2018, who subsequently received standard induction steroid therapy as the first-line treatment. Multivariable Cox proportional hazards models were used to identify the potential prognostic factors associated with relapse-free survival (RFS).

RESULTS

The median follow-up time for the entire cohort was 1913 days (interquartile range [IQR], 803-2929 days). Relapse occurred in 81.3% (39/48) patients during the follow-up period. The median time to relapse was 210 days (IQR, 140-308 days) after completion of standardized induction steroid therapy. Among the 17 parameters analyzed, Cox proportional hazard analysis identified whole-body total lesion glycolysis (WTLG) > 600 on FDG-PET as an independent risk factor for disease relapse (median RFS, 175 vs. 308 days; adjusted hazard ratio, 2.196 [95% confidence interval: 1.080-4.374]; P = 0.030).

CONCLUSION

WTLG on pretherapy FDG PET/CT was the only significant factor associated with RFS in IgG-RD patients receiving standard steroid induction therapy.

Artificial intelligence-driven radiomics study in cancer: the role of feature engineering and modeling

Modern medicine is reliant on various medical imaging technologies for non-invasively observing patients' anatomy. However, the interpretation of medical images can be highly subjective and dependent on the expertise of clinicians. Moreover, some potentially useful quantitative information in medical images, especially that which is not visible to the naked eye, is often ignored during clinical practice. In contrast, radiomics performs high-throughput feature extraction from medical images, which enables quantitative analysis of medical images and prediction of various clinical endpoints. Studies have reported that radiomics exhibits promising performance in diagnosis and predicting treatment responses and prognosis, demonstrating its potential to be a non-invasive auxiliary tool for personalized medicine. However, radiomics remains in a developmental phase as numerous technical challenges have yet to be solved, especially in feature engineering and statistical modeling. In this review, we introduce the current utility of radiomics by summarizing research on its application in the diagnosis, prognosis, and prediction of treatment responses in patients with cancer. We focus on machine learning approaches, for feature extraction and selection during feature engineering and for imbalanced datasets and multi-modality fusion during statistical modeling. Furthermore, we introduce the stability, reproducibility, and interpretability of features, and the generalizability and interpretability of models. Finally, we offer possible solutions to current challenges in radiomics research.

The added values of 18F-FDG PET/CT in differentiating cancer recurrence and osteoradionecrosis of mandible in patients with treated oral squamous cell carcinoma

BACKGROUND

Osteoradionecrosis (ORN) of the jaw requires a differential diagnosis to exclude cancer recurrence. Here, we sought to develop a scoring system comprising ¹⁸F-FDG PET/CT parameters for distinguishing between the two conditions in patients with oral squamous cell carcinoma (OSCC).

METHODS

The study consisted of 103 OSCC patients with suspected ORN of the jaw. All participants underwent ¹⁸F-FDG PET/CT imaging within 6 months of diagnostic histopathology. Following extraction of PET parameters, we identified clinical and imaging predictors of mandibular recurrence-free survival (MRFS) using receiver operating characteristic curve analysis and multivariate Cox regression models.

RESULTS

The results of histopathology revealed mandibular cancer recurrence in 24 patients (23.3%). Multivariate Cox regression analyses identified an age at diagnosis ≤ 52 years (P = 0.013), a location of the SUVmax voxel with soft tissue predominance (P = 0.019), and mandibular total lesion glycolysis (TLG) > 62.68 g (P < 0.001) as independent risk factors for MRFS. A scoring system was devised with scores from 0 (no risk factor) to 3 (presence of all three risk factors). High-risk patients with a score of 2-3 compared with score of 0-1 had a significantly higher likelihood of mandibular cancer recurrence (hazard ratio: 32.50, 95% confidence interval: 8.51-124.18, P < 0.001). The scoring system had a sensitivity of 87.50%, a specificity of 82.28%, and an accuracy of 83.50% for identifying mandibular cancer recurrence.

CONCLUSIONS

The scoring system of our study is clinically useful for identifying mandibular cancer recurrence in patients with suspected ORN of the jaw.

Influence of Image Processing on Radiomic Features From Magnetic Resonance Imaging

OBJECTIVE

Before implementing radiomics in routine clinical practice, comprehensive knowledge about the repeatability and reproducibility of radiomic features is required. The aim of this study was to systematically investigate the influence of image processing parameters on radiomic features from magnetic resonance imaging (MRI) in terms of feature values as well as test-retest repeatability.

MATERIALS AND METHODS

Utilizing a phantom consisting of 4 onions, 4 limes, 4 kiwifruits, and 4 apples, we acquired a test-retest dataset featuring 3 of the most commonly used MRI sequences on a 3 T scanner, namely, a T1-weighted, a T2-weighted, and a fluid-attenuated inversion recovery sequence, each at high and low resolution. After semiautomatic image segmentation, image processing with systematic variation of image processing parameters was performed, including spatial resampling, intensity discretization, and intensity rescaling. For each respective image processing setting, a total of 45 radiomic features were extracted, corresponding to the following 7 matrices/feature classes: conventional indices, histogram matrix, shape matrix, gray-level zone length matrix, gray-level run length matrix, neighboring gray-level dependence matrix, and gray-level cooccurrence matrix. Systematic differences of individual features between different resampling steps were assessed using 1-way analysis of variance with Tukey-type post hoc comparisons to adjust for multiple testing. Test-retest repeatability of radiomic features was measured using the concordance correlation coefficient, dynamic range, and intraclass correlation coefficient.

RESULTS

Image processing influenced radiological feature values. Regardless of the acquired sequence and feature class, significant differences ( P < 0.05) in feature values were found when the size of the resampled voxels was too large, that is, bigger than 3 mm. Almost all higher-order features depended strongly on intensity discretization. The effects of intensity rescaling were negligible except for some features derived from T1-weighted sequences. For all sequences, the percentage of repeatable features (concordance correlation coefficient and dynamic range ≥ 0.9) varied considerably depending on the image processing settings. The optimal image processing setting to achieve the highest percentage of stable features varied per sequence. Irrespective of image processing, the fluid-attenuated inversion recovery sequence in high-resolution overall yielded the highest number of stable features in comparison with the other sequences (89% vs 64%-78% for the respective optimal image processing settings). Across all sequences, the most repeatable features were generally obtained for a spatial resampling close to the originally acquired voxel size and an intensity discretization to at least 32 bins.

CONCLUSION

Variation of image processing parameters has a significant impact on the values of radiomic features as well as their repeatability. Furthermore, the optimal image processing parameters differ for each MRI sequence. Therefore, it is recommended that these processing parameters be determined in corresponding test-retest scans before clinical application. Extensive repeatability, reproducibility, and validation studies as well as standardization are required before quantitative image analysis and radiomics can be reliably translated into routine clinical care.

A heterogeneous phantom study for investigating the stability of PET images radiomic features with varying reconstruction settings

The purpose of this work was to assess the capability of radiomic features in distinguishing PET image regions with different uptake patterns. Furthermore, we assessed the stability of PET radiomic features with varying image reconstruction settings. An in-house phantom was designed and constructed, consisting of homogenous and heterogenous artificial phantom inserts. Four artificially constructed inserts were placed into a water filled phantom and filled with varying levels of radioactivity to simulate homogeneous and heterogeneous uptake patterns. The phantom was imaged for 80 min. PET images were reconstructed whilst varying reconstruction parameters. The parameters adjusted included, number of ordered subsets, number of iterations, use of time-of-flight and filter cut off. Regions of interest (ROI) were established by segmentation of the phantom inserts from the reconstructed images. In total seventy eight 3D radiomic features for each ROI with unique reconstructed parameters were extracted. The Friedman test was used to determine the statistical power of each radiomic feature in differentiating phantom inserts with different hetero/homogeneous configurations. The Coefficient of Variation (COV) of each feature, with respect to the reconstruction setting was used to determine feature stability. Forty three out of seventy eight radiomic features were found to be stable (COV ≤ 5%) against all reconstruction settings. To provide any utility, stable features are required to differentiate between regions with different hetro/homogeneity. Of the forty three stable features, fifteen (35%) features showed a statistically significant difference between the artificially constructed inserts. Such features included GLCM (Difference average, Difference entropy, Dissimilarity and Inverse difference), GLRL (Long run emphasis, Grey level non uniformity and Run percentage) and NGTDM (Complexity and Strength). The finding of this work suggests that radiomic features are capable of distinguishing between radioactive distribution patterns that demonstrate different levels of heterogeneity. Therefore, radiomic features could serve as an adjuvant diagnostic tool along with traditional imaging. However, the choice of the radiomic features needs to account for variability introduced when different reconstruction settings are used. Standardization of PET image reconstruction settings across sites performing radiomic analysis in multi-centre trials should be considered.

Impact of Aggregation Methods for Texture Features on Their Robustness Performance: Application to Nasopharyngeal 18F-FDG PET/CT

PURPOSE

This study aims to investigate the impact of aggregation methods used for the generation of texture features on their robustness of nasopharyngeal carcinoma (NPC) based on ¹⁸F-FDG PET/CT images.

METHODS

128 NPC patients were enrolled and 95 texture features were extracted for each patient including six feature families under different aggregation methods. For GLCM and GLRLM features, six aggregation methods were considered. For GLSZM, GLDZM, NGTDM and NGLDM features, three aggregation methods were considered. The robustness of the features affected by aggregation methods was assessed by the pair-wise intra-class correlation coefficient (ICC). Furthermore, the effects of discretization and partial volume correction (PVC) on the percent of ICC categories of all texture features were evaluated by overall ICC instead of the pair-wise ICC.

RESULTS

There were 12 features with excellent pair-wise ICCs varying aggregation methods, namely joint average, sum average, autocorrelation, long run emphasis, high grey level run emphasis, short run high grey level emphasis, long run high grey level emphasis, run length variance, SZM high grey level emphasis, DZM high grey level emphasis, high grey level count emphasis and dependence count percentage. For GLCM and GLRLM features, 19/25 and 14/16 features showed excellent pair-wise ICCs varying aggregation methods (averaged and merged) on the same dimensional features (2D, 2.5D or 3D). Different discretization levels and partial volume corrections lead to consistent robustness of textural features affected by aggregation methods.

CONCLUSION

Different dimensional features with the same aggregation methods showed worse robustness compared with the same dimensional features with different aggregation methods. Different discretization levels and PVC algorithms had a negligible effect on the percent of ICC categories of all texture features.

Feasibility of radiomic feature harmonization for pooling of [18F]FET or [18F]GE-180 PET images of gliomas

INTRODUCTION

Large datasets are required to ensure reliable non-invasive glioma assessment with radiomics-based machine learning methods. This can often only be achieved by pooling images from different centers. Moreover, trained models should perform with high accuracy when applied to data from different centers. In this study, the impact of reconstruction settings and segmentation methods on radiomic features derived from amino acid and TSPO PET images of glioma patients was examined. Additionally, the ability to model and thus reduce feature differences was investigated.

METHODS

[¹⁸F]FET and [¹⁸F]GE-180 PET data were acquired from 19 glioma patients. For each acquisition, 10 reconstruction settings and 9 segmentation methods were included to emulate multicentric data. Statistical robustness measures were calculated before and after ComBat harmonization. Differences between features due to setting variations were assessed using Friedman test, coefficient of variation (CV) and inter-rater reliability measures, including intraclass and Spearman's rank correlation coefficients and Fleiss' Kappa.

RESULTS

According to Friedman analyses, most features (>60%) showed significant differences. Yet, CV and inter-rater reliability measures indicated higher robustness. ComBat resulted in almost complete harmonization (>87%) according to Friedman test and little to no improvement according to CV and inter-rater reliability measures. [¹⁸F]GE-180 features were more sensitive to reconstruction settings than [¹⁸F]FET features.

CONCLUSIONS

According to Friedman test, feature distributions could be successfully aligned using ComBat. However, depending on settings, changes in patient ranks were observed for some features and could not be eliminated by harmonization. Thus, for clinical utilization it is recommended to exclude affected features.

Radiomics prognostic analysis of PET/CT images in a multicenter head and neck cancer cohort: investigating ComBat strategies, sub-volume characterization, and automatic segmentation

PURPOSE

This study aimed to investigate the impact of several ComBat harmonization strategies, intra-tumoral sub-volume characterization, and automatic segmentations for progression-free survival (PFS) prediction through radiomics modeling for patients with head and neck cancer (HNC) in PET/CT images.

METHODS

The HECKTOR MICCAI 2021 challenge set containing PET/CT images and clinical data of 325 oropharynx HNC patients was exploited. A total of 346 IBSI-compliant radiomic features were extracted for each patient's primary tumor volume defined by the reference manual contours. Modeling relied on least absolute shrinkage Cox regression (Lasso-Cox) for feature selection (FS) and Cox proportional-hazards (CoxPH) models were built to predict PFS. Within this methodological framework, 8 different strategies for ComBat harmonization were compared, including before or after FS, in feature groups separately or all features directly, and with center or clustering-determined labels. Features extracted from tumor sub-volume clustering were also investigated for their prognostic additional value. Finally, 3 automatic segmentations (2 threshold-based and a 3D U-Net) were also compared. All results were evaluated with the concordance index (C-index).

RESULTS

Radiomics features without harmonization, combined with clinical factors, led to models with C-index values of 0.69 in the testing set. The best version of ComBat harmonization, i.e., after FS, for feature groups separately and relying on clustering-determined labels, achieved a C-index of 0.71. The use of features extracted from tumor sub-volumes further improved the C-index to 0.72. Models that relied on the automatic segmentations yielded close but slightly lower prognostic performance (0.67-0.70) compared to reference contours.

CONCLUSION

A standard radiomics pipeline allowed for prediction of PFS in a multicenter HNC cohort. Applying a specific strategy of ComBat harmonization improved the performance. The extraction of intra-tumoral sub-volume features and automatic segmentation could contribute to the improvement and automation of prognosis modeling, respectively.

Preoperative prediction of mediastinal lymph node metastasis in non-small cell lung cancer based on 18F-FDG PET/CT radiomics

AIM

To establish and verify a 2-[¹⁸F]-fluoro-2-deoxy-d-glucose (FDG) positron-emission tomography (PET)/computed tomography (CT)-based radiomics nomogram to predict mediastinal lymph node metastasis (LNM) in non-small cell lung cancer (NSCLC) patients preoperatively.

MATERIALS AND METHODS

This retrospective study enrolled 155 NSCLC patients (primary cohort, n=93; validation cohort, n=62). For each patient, 2,704 radiomic features were extracted from the primary lung cancer regions. Four procedures including the Mann-Whitney U-test, Spearman's correlation analysis, minimum redundancy-maximum relevance (mRMR), and least absolute shrinkage and selection operator (LASSO) binary logistic regression were utilised for determining essential features and establishing a radiomics signature. After that, a nomogram was established. The nomogram's potential was assessed based on its discrimination, calibration, and clinical usefulness. The radiomics signature and nomogram predictive performances were evaluated with respect to the area under the receiver operating characteristic curve (AUC), specificity, accuracy, and sensitivity.

RESULTS

The radiomics signature composed of eight selected features had good discriminatory performance of LNM versus non-LNM groups an AUC of 0.851 and 0.826 in primary and validation cohorts, respectively. The nomogram also indicated good discrimination with an AUC of 0.869 and 0.847 in the primary and validation cohorts, respectively. Furthermore, good calibration was demonstrated utilising the nomogram.

CONCLUSIONS

An ¹⁸F-FDG PET/CT-based radiomics nomogram that integrates the radiomics signature and age was promoted to predict mediastinal LNM within NSCLC patients, which could potentially facilitate individualised therapy for mediastinal LNM before treatment. The nomogram was beneficial in clinical practice, as illustrated by decision curve analysis.

Harmonization of robust radiomic features in the submandibular gland using multi-ultrasound systems: a preliminary study

OBJECTIVE

This study aimed to identify robust radiomic features in multiultrasonography of the submandibular gland and normalize the interdevice discrepancies by applying a machine-learning-based harmonization method.

METHODS

Ultrasonographic images of normal submandibular gland of young healthy adults, aged between 20 and 40 years, were selected from two different devices. In a total of 30 images, the region of interest was determined along the border of gland parenchyma, and 103 radiomic features were extracted using A-VIEW. The coefficient of variation (CV) was obtained for individual features, and the features showing CV less than 10% were selected. For the selected features, the interdevice discrepancy was normalized using machine-learning method, called the ComBat harmonization. Median differences of the features between the two scanners, before and after harmonization, were compared using Mann-Whitney U-test; confidence interval of 95%.

RESULTS

Among total 103 radiomic features, 17 features were selected as robust, showing CV less than 10% in both scanners. All values of selected features, except two, showed a statistical difference between the two devices. After applying the ComBat harmonization method, the median and distribution of the 16 features were harmonized to show no significant difference between the two scanners (p > 0.05). One feature remained different (p ≤ 0.05).

CONCLUSION

On ultrasonographic examination, robust radiomic features for normal submandibular gland were obtained and interdevice normalization was efficiently conducted using ComBat harmonization. Our findings would be useful for multidevices or multicenter studies based on clinical ultrasonographic imaging data to improve the accuracy of the overall diagnostic model.

Role of Machine Learning in Precision Oncology: Applications in Gastrointestinal Cancers

Gastrointestinal (GI) cancers, consisting of a wide spectrum of pathologies, have become a prominent health issue globally. Despite medical imaging playing a crucial role in the clinical workflow of cancers, standard evaluation of different imaging modalities may provide limited information. Accurate tumor detection, characterization, and monitoring remain a challenge. Progress in quantitative imaging analysis techniques resulted in "radiomics", a promising methodical tool that helps to personalize diagnosis and treatment optimization. Radiomics, a sub-field of computer vision analysis, is a bourgeoning area of interest, especially in this era of precision medicine. In the field of oncology, radiomics has been described as a tool to aid in the diagnosis, classification, and categorization of malignancies and to predict outcomes using various endpoints. In addition, machine learning is a technique for analyzing and predicting by learning from sample data, finding patterns in it, and applying it to new data. Machine learning has been increasingly applied in this field, where it is being studied in image diagnosis. This review assesses the current landscape of radiomics and methodological processes in GI cancers (including gastric, colorectal, liver, pancreatic, neuroendocrine, GI stromal, and rectal cancers). We explain in a stepwise fashion the process from data acquisition and curation to segmentation and feature extraction. Furthermore, the applications of radiomics for diagnosis, staging, assessment of tumor prognosis and treatment response according to different GI cancer types are explored. Finally, we discussed the existing challenges and limitations of radiomics in abdominal cancers and investigate future opportunities.

Impact of image quality on radiomics applications

Radiomics features extracted from medical images have been widely reported to be useful in the patient specific outcome modeling for variety of assessment and prediction purposes. Successful application of radiomics features as imaging biomarkers, however, is dependent on the robustness of the approach to the variation in each step of the modeling workflow. Variation in the input image quality is one of the main sources that impacts the reproducibility of radiomics analysis when a model is applied to broader range of medical imaging data. The quality of medical image is generally affected by both the scanner related factors such as image acquisition/reconstruction settings and the patient related factors such as patient motion. This article aimed to review the published literatures in this field that reported the impact of various imaging factors on the radiomics features through the change in image quality. The literatures were categorized by different imaging modalities and also tabulated based on the imaging parameters and the class of radiomics features included in the study. Strategies for image quality standardization were discussed based on the relevant literatures and recommendations for reducing the impact of image quality variation on the radiomics in multi-institutional clinical trial were summarized at the end of this article.

The latest application progress of radiomics in prediction and diagnosis of liver diseases

INTRODUCTION

Early detection and individualized treatment of patients with liver disease is the key to survival. Radiomics can extract high-throughput quantitative features by multimode imaging, which has good application prospects for the diagnosis, staging and prognosis of benign and malignant liver diseases. Therefore, this paper summarizes the current research status in the field of liver disease, in order to help these patients achieve personalized and precision medical care.

AREAS COVERED

This paper uses several keywords on the PubMed database to search the references, and reviews the workflow of traditional radiomics, as well as the characteristics and influencing factors of different imaging modes. At the same time, the references on the application of imaging in different benign and malignant liver diseases were also summarized.

EXPERT OPINION

For patients with liver disease, the traditional imaging evaluation can only provide limited information. Radiomics exploits the characteristics of high-throughput and high-dimensional extraction, enabling liver imaging capabilities far beyond the scope of traditional visual image analysis. Recent studies have demonstrated the prospect of this technology in personalized diagnosis and treatment decision in various fields of the liver. However, further clinical validation is needed in its application and practice.

Combatting the effect of image reconstruction settings on lymphoma [18F]FDG PET metabolic tumor volume assessment using various segmentation methods

Background

[¹⁸F]FDG PET-based metabolic tumor volume (MTV) is a promising prognostic marker for lymphoma patients. The aim of this study is to assess the sensitivity of several MTV segmentation methods to variations in image reconstruction methods and the ability of ComBat to improve MTV reproducibility.

Methods

Fifty-six lesions were segmented from baseline [¹⁸F]FDG PET scans of 19 lymphoma patients. For each scan, EARL1 and EARL2 standards and locally clinically preferred reconstruction protocols were applied. Lesions were delineated using 9 semiautomatic segmentation methods: fixed threshold based on standardized uptake value (SUV), (SUV = 4, SUV = 2.5), relative threshold (41% of SUVmax [41M], 50% of SUVpeak [A50P]), majority vote-based methods that select voxels detected by at least 2 (MV2) and 3 (MV3) out of the latter 4 methods, Nestle thresholding, and methods that identify the optimal method based on SUVmax (L2A, L2B). MTVs from EARL2 and locally clinically preferred reconstructions were compared to those from EARL1. Finally, different versions of ComBat were explored to harmonize the data.

Results

MTVs from the SUV4.0 method were least sensitive to the use of different reconstructions (MTV ratio: median = 1.01, interquartile range = [0.96–1.10]). After ComBat harmonization, an improved agreement of MTVs among different reconstructions was found for most segmentation methods. The regular implementation of ComBat (‘Regular ComBat’) using non-transformed distributions resulted in less accurate and precise MTV alignments than a version using log-transformed datasets (‘Log-transformed ComBat’).

Conclusion

MTV depends on both segmentation method and reconstruction methods. ComBat reduces reconstruction dependent MTV variability, especially when log-transformation is used to account for the non-normal distribution of MTVs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13550-022-00916-9.

Repeatability and reproducibility of magnetic resonance imaging-based radiomic features in rectal cancer

Purpose: Radiomics of magnetic resonance images (MRIs) in rectal cancer can non-invasively characterize tumor heterogeneity with potential to discover new imaging biomarkers. However, for radiomics to be reliable, the imaging features measured must be stable and reproducible. The aim of this study is to quantify the repeatability and reproducibility of MRI-based radiomic features in rectal cancer. Approach: An MRI radiomics phantom was used to measure the longitudinal repeatability of radiomic features and the impact of post-processing changes related to image resolution and noise. Repeatability measurements in rectal cancers were also quantified in a cohort of 10 patients with test-retest imaging among two observers. Results: We found that many radiomic features, particularly from texture classes, were highly sensitive to changes in image resolution and noise. About 49% of features had coefficient of variations ≤ 10 % in longitudinal phantom measurements. About 75% of radiomic features in in vivo test-retest measurements had an intraclass correlation coefficient of ≥ 0.8 . We saw excellent interobserver agreement with mean Dice similarity coefficient of 0.95 ± 0.04 for test and retest scans. Conclusions: The results of this study show that even when using a consistent imaging protocol many radiomic features were unstable. Therefore, caution must be taken when selecting features for potential imaging biomarkers.

Impact of Interfractional Error on Dosiomic Features

Objectives

The purpose of this study was to investigate the stability of dosiomic features under random interfractional error. We investigated the differences in the values of features with different fractions and the error in the values of dosiomic features under interfractional error.

Material and Methods

The isocenters of the treatment plans of 15 lung cancer patients were translated by a maximum of ±3 mm in each axis with a mean of (0, 0, 0) and a standard deviation of (1.2, 1.2, 1.2) mm in the x, y, and z directions for each fraction. A total of 81 dose distributions for each patient were then calculated considering four fraction number groups (2, 10, 20, and 30). A total of 93 dosiomic features were extracted from each dose distribution in four different regions of interest (ROIs): gross tumor volume (GTV), planning target volume (PTV), heart, and both lungs. The stability of dosiomic features was analyzed for each fraction number group by the coefficient of variation (CV) and intraclass correlation coefficient (ICC). The agreements in the means of dosiomic features among the four fraction number groups were tested by ICC. The percent differences (PD) between the dosiomic features extracted from the original dose distribution and the dosiomic features extracted from the dose distribution with interfractional error were calculated.

Results

Eleven out of 93 dosiomic features demonstrated a large CV (CV ≥ 20%). Overall CV values were highest in GTV ROIs and lowest in lung ROIs. The stability of dosiomic features decreased as the total number of fractions decreased. The ICC results showed that five out of 93 dosiomic features had an ICC lower than 0.75, which indicates intermediate or poor stability under interfractional error. The mean dosiomic feature values were shown to be consistent with different numbers of fractions (ICC ≥ 0.9). Some of the dosiomic features had PD greater than 50% and showed different PD values with different numbers of fractions.

Conclusion

Some dosiomic features have low stability under interfractional error. The stability and values of the dosiomic features were affected by the total number of fractions. The effect of interfractional error on dosiomic features should be considered in further studies regarding dosiomics for reproducible results.

Development and validation of a radiomic model for the diagnosis of dopaminergic denervation on [18F]FDOPA PET/CT

Purpose

FDOPA PET shows good performance for the diagnosis of striatal dopaminergic denervation, making it a valuable tool for the differential diagnosis of Parkinsonism. Textural features are image biomarkers that could potentially improve the early diagnosis and monitoring of neurodegenerative parkinsonian syndromes. We explored the performances of textural features for binary classification of FDOPA scans.

Methods

We used two FDOPA PET datasets: 443 scans for feature selection, and 100 scans from a different PET/CT system for model testing. Scans were labelled according to expert interpretation (dopaminergic denervation versus no dopaminergic denervation). We built LASSO logistic regression models using 43 biomarkers including 32 textural features. Clinical data were also collected using a shortened UPDRS scale.

Results

The model built from the clinical data alone had a mean area under the receiver operating characteristics (AUROC) of 63.91. Conventional imaging features reached a maximum score of 93.47 but the addition of textural features significantly improved the AUROC to 95.73 (p < 0.001), and 96.10 (p < 0.001) when limiting the model to the top three features: GLCM_Correlation, Skewness and Compacity. Testing the model on the external dataset yielded an AUROC of 96.00, with 95% sensitivity and 97% specificity. GLCM_Correlation was one of the most independent features on correlation analysis, and systematically had the heaviest weight in the classification model.

Conclusion

A simple model with three radiomic features can identify pathologic FDOPA PET scans with excellent sensitivity and specificity. Textural features show promise for the diagnosis of parkinsonian syndromes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00259-022-05816-7.

A cross-scanner and cross-tracer deep learning method for the recovery of standard-dose imaging quality from low-dose PET

PURPOSE

A critical bottleneck for the credibility of artificial intelligence (AI) is replicating the results in the diversity of clinical practice. We aimed to develop an AI that can be independently applied to recover high-quality imaging from low-dose scans on different scanners and tracers.

METHODS

Brain [18F]FDG PET imaging of 237 patients scanned with one scanner was used for the development of AI technology. The developed algorithm was then tested on [18F]FDG PET images of 45 patients scanned with three different scanners, [18F]FET PET images of 18 patients scanned with two different scanners, as well as [18F]Florbetapir images of 10 patients. A conditional generative adversarial network (GAN) was customized for cross-scanner and cross-tracer optimization. Three nuclear medicine physicians independently assessed the utility of the results in a clinical setting.

RESULTS

The improvement achieved by AI recovery significantly correlated with the baseline image quality indicated by structural similarity index measurement (SSIM) (r = -0.71, p < 0.05) and normalized dose acquisition (r = -0.60, p < 0.05). Our cross-scanner and cross-tracer AI methodology showed utility based on both physical and clinical image assessment (p < 0.05).

CONCLUSION

The deep learning development for extensible application on unknown scanners and tracers may improve the trustworthiness and clinical acceptability of AI-based dose reduction.

A tri-modal tissue-equivalent anthropomorphic phantom for PET, CT and multi-parametric MRI radiomics

PURPOSE

Radiomics has emerged as an advanced image processing methodology to define quantitative imaging biomarkers for prognosis and prediction of treatment response and outcome. The development of quantitative imaging biomarkers requires careful analysis to define their accuracy, stability and reproducibility through phantom measurements. Few efforts were devoted to develop realistic anthropomorphic phantoms. In this work, we developed a multimodality image phantom suitable for PET, CT and multiparametric MRI imaging.

METHODS

A tissue-equivalent gel-based mixture was designed and tested for compatibility with different imaging modalities. Calibration measurements allowed to assess gel composition to simulate PET, CT and MRI contrasts of oncological lesions. The characterized gel mixture was used to create realistic synthetic lesions (e.g. lesions with irregular shape and non-uniform image contrast), to be inserted in a standard anthropomorphic phantom. In order to show phantom usefulness, issues related to accuracy, stability and reproducibility of radiomic biomarkers were addressed as proofs-of-concept.

RESULTS

The procedure for gel preparation was straightforward and the characterized gel mixture allowed to mime simultaneously oncological lesion contrast in CT, PET and MRI imaging. Proofs-of-concept studies suggested that phantom measurements can be customized for specific clinical situations and radiomic protocols.

CONCLUSIONS

We developed a strategy to manufacture an anthropomorphic, tissue-equivalent, multimodal phantom to be customized on specific radiomics protocols, for addressing specific methodological issues both in mono and multicentric studies.

Prediction of Non-Response to Neoadjuvant Chemoradiotherapy in Esophageal Cancer Patients with 18F-FDG PET Radiomics Based Machine Learning Classification

Background: Approximately 26% of esophageal cancer (EC) patients do not respond to neoadjuvant chemoradiotherapy (nCRT), emphasizing the need for pre-treatment selection. The aim of this study was to predict non-response using a radiomic model on baseline ¹⁸F-FDG PET. Methods: Retrospectively, 143 ¹⁸F-FDG PET radiomic features were extracted from 199 EC patients (T1N1-3M0/T2–4aN0-3M0) treated between 2009 and 2019. Non-response (n = 57; 29%) was defined as Mandard Tumor Regression Grade 4–5 (n = 44; 22%) or interval progression (n = 13; 7%). Randomly, 139 patients (70%) were allocated to explore all combinations of 24 feature selection strategies and 6 classification methods towards the cross-validated average precision (AP). The predictive value of the best-performing model, i.e AP and area under the ROC curve analysis (AUC), was evaluated on an independent test subset of 60 patients (30%). Results: The best performing model had an AP (mean ± SD) of 0.47 ± 0.06 on the training subset, achieved by a support vector machine classifier trained on five principal components of relevant clinical and radiomic features. The model was externally validated with an AP of 0.66 and an AUC of 0.67. Conclusion: In the present study, the best-performing model on pre-treatment ¹⁸F-FDG PET radiomics and clinical features had a small clinical benefit to identify non-responders to nCRT in EC.

Beads phantom for evaluating heterogeneity of SUV on 18F-FDG PET images

PURPOSE

This study aimed to develop a dedicated phantom using acrylic beads for texture analysis and to represent heterogeneous ¹⁸F-fluorodeoxyglucose (FDG) distributions in various acquisition periods.

METHODS

Images of acrylic spherical beads with or without diameters of 5- and 10-mm representing heterogeneous and homogeneous ¹⁸F-FDG distribution in phantoms, respectively, were collected for 20 min in list mode. Phantom data were reconstructed using three-dimensional ordered subset expectation maximization with attenuation and scatter corrections, and the time-of-flight algorithm. The beads phantom images were acquired twice to evaluate the robustness of texture features. Thirty-one texture features were extracted, and the robustness of texture feature values was evaluated by calculating the percentage of coefficient of variation (%COV) and intraclass coefficient of correlation (ICC). Cross-correlation coefficients among texture feature values were clustered to classify the characteristics of these features.

RESULTS

Heterogeneous ¹⁸F-FDG distribution was represented by the beads phantom images. The agreements of %COV between two measurements were acceptable (ICC ≥ 0.71). All texture features were classified into four groups. Among 31 texture features, 24 exhibited significant different values between phantoms with and without beads in 1-, 2-, 3-, 4-, 5-, 20-min image acquisitions. Whereas, the homogeneous and heterogeneous ¹⁸F-FDG distribution could not be discriminated by seven texture features: low gray-level run emphasis, high gray-level run emphasis, short-run low gray-level emphasis, low gray-level zone emphasis, high gray-level zone emphasis, short-zone low gray-level emphasis, and coarseness.

CONCLUSIONS

We have developed the acrylic beads phantom for texture analysis that could represent heterogeneous ¹⁸F-FDG distributions in various acquisition periods. Most texture features could discriminate homogeneous and heterogeneous ¹⁸F-FDG distributions in the beads phantom images.

Synergistic impact of motion and acquisition/reconstruction parameters on 18 F-FDG PET radiomic features in non-small cell lung cancer: phantom and clinical studies

Objectives

This study is aimed at examining the synergistic impact of motion and acquisition/reconstruction parameters on ¹⁸F‐FDG PET image radiomic features in non‐small cell lung cancer (NSCLC) patients, and investigating the robustness of features performance in differentiating NSCLC histopathology subtypes.

Methods

An in‐house developed thoracic phantom incorporating lesions with different sizes was used with different reconstruction settings, including various reconstruction algorithms, number of subsets and iterations, full‐width at half‐maximum of post‐reconstruction smoothing filter and acquisition parameters, including injected activity and test–retest with and without motion simulation. To simulate motion, a special motor was manufactured to simulate respiratory motion based on a normal patient in two directions. The lesions were delineated semi‐automatically to extract 174 radiomic features. All radiomic features were categorized according to the coefficient of variation (COV) to select robust features. A cohort consisting of 40 NSCLC patients with adenocarcinoma (n = 20) and squamous cell carcinoma (n = 20) was retrospectively analyzed. Statistical analysis was performed to discriminate robust features in differentiating histopathology subtypes of NSCLC lesions.

Results

Overall, 29% of radiomic features showed a COV ≤5% against motion. Forty‐five percent and 76% of the features showed a COV ≤ 5% against the test–retest with and without motion in large lesions, respectively. Thirty‐three percent and 45% of the features showed a COV ≤ 5% against different reconstruction parameters with and without motion, respectively. For NSCLC histopathological subtype differentiation, statistical analysis showed that 31 features were significant (p‐value < 0.05). Two out of the 31 significant features, namely, the joint entropy of GLCM (AUC = 0.71, COV = 0.019) and median absolute deviation of intensity histogram (AUC = 0.7, COV = 0.046), were robust against the motion (same reconstruction setting).

Conclusions

Motion, acquisition, and reconstruction parameters significantly impact radiomic features, just as their synergies. Radiomic features with high predictive performance (statistically significant) in differentiating histopathological subtype of NSCLC may be eliminated due to non‐reproducibility.

Multitask Learning Radiomics on Longitudinal Imaging to Predict Survival Outcomes following Risk-Adaptive Chemoradiation for Non-Small Cell Lung Cancer

Medical imaging provides quantitative and spatial information to evaluate treatment response in the management of patients with non-small cell lung cancer (NSCLC). High throughput extraction of radiomic features on these images can potentially phenotype tumors non-invasively and support risk stratification based on survival outcome prediction. The prognostic value of radiomics from different imaging modalities and time points prior to and during chemoradiation therapy of NSCLC, relative to conventional imaging biomarker or delta radiomics models, remains uncharacterized. We investigated the utility of multitask learning of multi-time point radiomic features, as opposed to single-task learning, for improving survival outcome prediction relative to conventional clinical imaging feature model benchmarks. Survival outcomes were prospectively collected for 45 patients with unresectable NSCLC enrolled on the FLARE-RT phase II trial of risk-adaptive chemoradiation and optional consolidation PD-L1 checkpoint blockade (NCT02773238). FDG-PET, CT, and perfusion SPECT imaging pretreatment and week 3 mid-treatment was performed and 110 IBSI-compliant pyradiomics shape-/intensity-/texture-based features from the metabolic tumor volume were extracted. Outcome modeling consisted of a fused Laplacian sparse group LASSO with component-wise gradient boosting survival regression in a multitask learning framework. Testing performance under stratified 10-fold cross-validation was evaluated for multitask learning radiomics of different imaging modalities and time points. Multitask learning models were benchmarked against conventional clinical imaging and delta radiomics models and evaluated with the concordance index (c-index) and index of prediction accuracy (IPA). FDG-PET radiomics had higher prognostic value for overall survival in test folds (c-index 0.71 [0.67, 0.75]) than CT radiomics (c-index 0.64 [0.60, 0.71]) or perfusion SPECT radiomics (c-index 0.60 [0.57, 0.63]). Multitask learning of pre-/mid-treatment FDG-PET radiomics (c-index 0.71 [0.67, 0.75]) outperformed benchmark clinical imaging (c-index 0.65 [0.59, 0.71]) and FDG-PET delta radiomics (c-index 0.52 [0.48, 0.58]) models. Similarly, the IPA for multitask learning FDG-PET radiomics (30%) was higher than clinical imaging (26%) and delta radiomics (15%) models. Radiomics models performed consistently under different voxel resampling conditions. Multitask learning radiomics for outcome modeling provides a clinical decision support platform that leverages longitudinal imaging information. This framework can reveal the relative importance of different imaging modalities and time points when designing risk-adaptive cancer treatment strategies.

Robustness of CT radiomics features: consistency within and between single-energy CT and dual-energy CT

Objectives

To evaluate inter- and intra- scan mode and scanner repeatability and reproducibility of radiomics features within and between single-energy CT (SECT) and dual-energy CT (DECT).

Methods

A standardized phantom with sixteen rods of clinical-relevant densities was scanned on seven DECT-capable scanners and three SECT-only scanners. The acquisition parameters were selected to present typical abdomen-pelvic examinations with the same voxel size. Images of SECT at 120 kVp and corresponding 120 kVp-like virtual monochromatic images (VMIs) in DECT which were generated according to scanners were analyzed. Regions of interest were drawn with rigid registrations to avoid variations due to segmentation. Radiomics features were extracted via Pyradiomics platform. Test-retest repeatability was evaluated by Bland-Altman analysis for repeated scans. Intra-scanner reproducibility for different scan modes was tested by intraclass correlation coefficient (ICC) and concordance correlation coefficient (CCC). Inter-scanner reproducibility among different scanners for same scan mode was assessed by coefficient of variation (CV) and quartile coefficient of dispersion (QCD).

Results

The test-retest analysis presented that 92.91% and 87.02% of the 94 assessed features were repeatable for SECT 120kVp and DECT 120 kVp-like VMIs, respectively. The intra-scanner analysis for SECT 120kVp vs DECT 120 kVp-like VMIs demonstrated that 10.76% and 10.28% of features were with ICC > 0.90 and CCC > 0.90, respectively. The inter-scanner analysis showed that 17.09% and 27.73% of features for SECT 120kVp were with CV < 10% and QCD < 10%, and 15.16% and 32.78% for DECT 120 kVp-like VMIs, respectively.

Conclusions

The majority of radiomics features were non-reproducible within and between SECT and DECT.

Key Points

• Although the test-retest analysis showed high repeatability for radiomics features, the overall reproducibility of radiomics features within and between SECT and DECT was low.

• Only about one-tenth of radiomics features extracted from SECT images and corresponding DECT images did match each other, even their average photon energy levels were considered alike, indicating that the scan mode potentially altered the radiomics features.

• Less than one-fifth of radiomics features were reproducible among multiple SECT and DECT scanners, regardless of their fixed acquisition and reconstruction parameters, suggesting the necessity of scanning protocol adjustment and post-scan harmonization process.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00330-022-08628-3.

Design and evaluation of a modular multimodality imaging phantom to simulate heterogeneous uptake and enhancement patterns for radiomic quantification in hybrid imaging; a feasibility study

Background

Accuracy and precision assessment in radiomic features is important for the determination of their potential to characterize cancer lesions. In this regard, simulation of different imaging conditions using specialized phantoms is increasingly being investigated. In this study, the design and evaluation of a modular multimodality imaging phantom to simulate heterogeneous uptake and enhancement patterns for radiomics quantification in hybrid imaging is presented.

Methods

A modular multimodality imaging phantom was constructed that could simulate different patterns of heterogeneous uptake and enhancement patterns in positron emission tomography (PET), single‐photon emission computed tomography (SPECT), computed tomography (CT), and magnetic resonance (MR) imaging. The phantom was designed to be used as an insert in the standard NEMA‐NU2 IEC body phantom casing. The entire phantom insert is composed of three segments, each containing three separately fillable compartments. The fillable compartments between segments had different sizes in order to simulate heterogeneous patterns at different spatial scales. The compartments were separately filled with different ratios of ^99mTc‐pertechnetate, ¹⁸F‐fluorodeoxyglucose ([¹⁸F]FDG), iodine‐ and gadolinium‐based contrast agents for SPECT, PET, CT, and T₁‐weighted MR imaging respectively. Image acquisition was performed using standard oncological protocols on all modalities and repeated five times for repeatability assessment. A total of 93 radiomic features were calculated. Variability was assessed by determining the coefficient of quartile variation (CQV) of the features. Comparison of feature repeatability at different modalities and spatial scales was performed using Kruskal‐Wallis‐, Mann‐Whitney U‐, one‐way ANOVA‐ and independent t‐tests.

Results

Heterogeneous uptake and enhancement could be simulated on all four imaging modalities. Radiomic features in SPECT were significantly less stable than in all other modalities. Features in PET were significantly less stable than in MR and CT. A total of 20 features, particularly in the gray‐level co‐occurrence matrix (GLCM) and gray‐level run‐length matrix (GLRLM) class, were found to be relatively stable in all four modalities for all three spatial scales of heterogeneous patterns (with CQV < 10%).

Conclusion

The phantom was suitable for simulating heterogeneous uptake and enhancement patterns in [¹⁸F]FDG‐PET, ^99mTc‐SPECT, CT, and T₁‐weighted MR images. The results of this work indicate that the phantom might be useful for the further development and optimization of imaging protocols for radiomic quantification in hybrid imaging modalities.

A review on advances in 18F-FDG PET/CT radiomics standardisation and application in lung disease management

Radiomics analysis quantifies the interpolation of multiple and invisible molecular features present in diagnostic and therapeutic images. Implementation of 18-fluorine-fluorodeoxyglucose positron emission tomography/computed tomography (¹⁸F-FDG PET/CT) radiomics captures various disorders in non-invasive and high-throughput manner. ¹⁸F-FDG PET/CT accurately identifies the metabolic and anatomical changes during cancer progression. Therefore, the application of ¹⁸F-FDG PET/CT in the field of oncology is well established. Clinical application of ¹⁸F-FDG PET/CT radiomics in lung infection and inflammation is also an emerging field. Combination of bioinformatics approaches or textual analysis allows radiomics to extract additional information to predict cell biology at the micro-level. However, radiomics texture analysis is affected by several factors associated with image acquisition and processing. At present, researchers are working on mitigating these interrupters and developing standardised workflow for texture biomarker establishment. This review article focuses on the application of ¹⁸F-FDG PET/CT in detecting lung diseases specifically on cancer, infection and inflammation. An overview of different approaches and challenges encountered on standardisation of ¹⁸F-FDG PET/CT technique has also been highlighted. The review article provides insights about radiomics standardisation and application of ¹⁸F-FDG PET/CT in lung disease management.