Performance comparison of modified ComBat for harmonization of radiomic features for multicenter studies

Identifying novel data-driven subgroups in congenital heart disease using multi-modal measures of brain structure

Individuals with congenital heart disease (CHD) have an increased risk of neurodevelopmental impairments. Given the hypothesized complexity linking genomics, atypical brain structure, cardiac diagnoses and their management, and neurodevelopmental outcomes, unsupervised methods may provide unique insight into neurodevelopmental variability in CHD. Using data from the Pediatric Cardiac Genomics Consortium Brain and Genes study, we identified data-driven subgroups of individuals with CHD from measures of brain structure. Using structural magnetic resonance imaging (MRI; N = 93; cortical thickness, cortical volume, and subcortical volume), we identified subgroups that differed primarily on cardiac anatomic lesion and language ability. In contrast, using diffusion MRI (N = 88; white matter connectivity strength), we identified subgroups that were characterized by differences in associations with rare genetic variants and visual-motor function. This work provides insight into the differential impacts of cardiac lesions and genomic variation on brain growth and architecture in patients with CHD, with potentially distinct effects on neurodevelopmental outcomes.

Impact of harmonization on the reproducibility of MRI radiomic features when using different scanners, acquisition parameters, and image pre-processing techniques: a phantom study

This study investigated the impact of ComBat harmonization on the reproducibility of radiomic features extracted from magnetic resonance images (MRI) acquired on different scanners, using various data acquisition parameters and multiple image pre-processing techniques using a dedicated MRI phantom. Four scanners were used to acquire an MRI of a nonanatomic phantom as part of the TCIA RIDER database. In fast spin-echo inversion recovery (IR) sequences, several inversion durations were employed, including 50, 100, 250, 500, 750, 1000, 1500, 2000, 2500, and 3000 ms. In addition, a 3D fast spoiled gradient recalled echo (FSPGR) sequence was used to investigate several flip angles (FA): 2, 5, 10, 15, 20, 25, and 30 degrees. Nineteen phantom compartments were manually segmented. Different approaches were used to pre-process each image: Bin discretization, Wavelet filter, Laplacian of Gaussian, logarithm, square, square root, and gradient. Overall, 92 first-, second-, and higher-order statistical radiomic features were extracted. ComBat harmonization was also applied to the extracted radiomic features. Finally, the Intraclass Correlation Coefficient (ICC) and Kruskal-Wallis's (KW) tests were implemented to assess the robustness of radiomic features. The number of non-significant features in the KW test ranged between 0-5 and 29-74 for various scanners, 31-91 and 37-92 for three times tests, 0-33 to 34-90 for FAs, and 3-68 to 65-89 for IRs before and after ComBat harmonization, with different image pre-processing techniques, respectively. The number of features with ICC over 90% ranged between 0-8 and 6-60 for various scanners, 11-75 and 17-80 for three times tests, 3-83 to 9-84 for FAs, and 3-49 to 3-63 for IRs before and after ComBat harmonization, with different image pre-processing techniques, respectively. The use of various scanners, IRs, and FAs has a great impact on radiomic features. However, the majority of scanner-robust features is also robust to IR and FA. Among the effective parameters in MR images, several tests in one scanner have a negligible impact on radiomic features. Different scanners and acquisition parameters using various image pre-processing might affect radiomic features to a large extent. ComBat harmonization might significantly impact the reproducibility of MRI radiomic features.

Impact of MRI radiomic feature normalization for prognostic modelling in uterine endometrial and cervical cancers

Widespread clinical use of MRI radiomic tumor profiling for prognostication and treatment planning in cancers faces major obstacles due to limitations in standardization of radiomic features. The purpose of the current work was to assess the impact of different MRI scanning- and normalization protocols for the statistical analyses of tumor radiomic data in two patient cohorts with uterine endometrial-(EC) (n = 136) and cervical (CC) (n = 132) cancer. 1.5 T and 3 T, T1-weighted MRI 2 min post-contrast injection, T2-weighted turbo spin echo imaging, and diffusion-weighted imaging were acquired. Radiomic features were extracted from within manually segmented tumors in 3D and normalized either using z-score normalization or a linear regression model (LRM) accounting for linear dependencies with MRI acquisition parameters. Patients were clustered into two groups based on radiomic profile. Impact of MRI scanning parameters on cluster composition and prognostication were analyzed using Kruskal-Wallis tests, Kaplan-Meier plots, log-rank test, random survival forests and LASSO Cox regression with time-dependent area under curve (tdAUC) (α = 0.05). A large proportion of the radiomic features was statistically associated with MRI scanning protocol in both cohorts (EC: 162/385 [42%]; CC: 180/292 [62%]). A substantial number of EC (49/136 [36%]) and CC (50/132 [38%]) patients changed cluster when clustering was performed after z-score-versus LRM normalization. Prognostic modeling based on cluster groups yielded similar outputs for the two normalization methods in the EC/CC cohorts (log-rank test; z-score: p = 0.02/0.33; LRM: p = 0.01/0.45). Mean tdAUC for prognostic modeling of disease-specific survival (DSS) by the radiomic features in EC/CC was similar for the two normalization methods (random survival forests; z-score: mean tdAUC = 0.77/0.78; LRM: mean tdAUC = 0.80/0.75; LASSO Cox; z-score: mean tdAUC = 0.64/0.76; LRM: mean tdAUC = 0.76/0.75). Severe biases in tumor radiomics data due to MRI scanning parameters exist. Z-score normalization does not eliminate these biases, whereas LRM normalization effectively does. Still, radiomic cluster groups after z-score- and LRM normalization were similarly associated with DSS in EC and CC patients.

Automated PD-L1 status prediction in lung cancer with multi-modal PET/CT fusion

Programmed death-ligand 1 (PD-L1) expressions play a crucial role in guiding therapeutic interventions such as the use of tyrosine kinase inhibitors (TKIs) and immune checkpoint inhibitors (ICIs) in lung cancer. Conventional determination of PD-L1 status includes careful surgical or biopsied tumor specimens. These specimens are gathered through invasive procedures, representing a risk of difficulties and potential challenges in getting reliable and representative tissue samples. Using a single center cohort of 189 patients, our objective was to evaluate various fusion methods that used non-invasive computed tomography (CT) and18 F-FDG positron emission tomography (PET) images as inputs to various deep learning models to automatically predict PD-L1 in non-small cell lung cancer (NSCLC). We compared three different architectures (ResNet, DenseNet, and EfficientNet) and considered different input data (CT only, PET only, PET/CT early fusion, PET/CT late fusion without as well as with partially and fully shared weights to determine the best model performance. Models were assessed utilizing areas under the receiver operating characteristic curves (AUCs) considering their 95% confidence intervals (CI). The fusion of PET and CT images as input yielded better performance for PD-L1 classification. The different data fusion schemes systematically outperformed their individual counterparts when used as input of the various deep models. Furthermore, early fusion consistently outperformed late fusion, probably as a result of its capacity to capture more complicated patterns by merging PET and CT derived content at a lower level. When we looked more closely at the effects of weight sharing in late fusion architectures, we discovered that while it might boost model stability, it did not always result in better results. This suggests that although weight sharing could be beneficial when modality parameters are similar, the anatomical and metabolic information provided by CT and PET scans are too dissimilar to consistently lead to improved PD-L1 status predictions.

Review on radiomic analysis in 18F-fluorodeoxyglucose positron emission tomography for prediction of melanoma outcomes

Over the past decade, several strategies have revolutionized the clinical management of patients with cutaneous melanoma (CM), including immunotherapy and targeted tyrosine kinase inhibitor (TKI)-based therapies. Indeed, immune checkpoint inhibitors (ICIs), alone or in combination, represent the standard of care for patients with advanced disease without an actionable mutation. Notably BRAF combined with MEK inhibitors represent the therapeutic standard for disease disclosing BRAF mutation. At the same time, FDG PET/CT has become part of the routine staging and evaluation of patients with cutaneous melanoma. There is growing interest in using FDG PET/CT measurements to predict response to ICI therapy and/or target therapy. While semiquantitative values such as standardized uptake value (SUV) are limited for predicting outcome, new measures including tumor metabolic volume, total lesion glycolysis and radiomics seem promising as potential imaging biomarkers for nuclear medicine. The aim of this review, prepared by an interdisciplinary group of experts, is to take stock of the current literature on radiomics approaches that could improve outcomes in CM.

The effect of harmonization on the variability of PET radiomic features extracted using various segmentation methods

PURPOSE

This study aimed to examine the robustness of positron emission tomography (PET) radiomic features extracted via different segmentation methods before and after ComBat harmonization in patients with non-small cell lung cancer (NSCLC).

METHODS

We included 120 patients (positive recurrence = 46 and negative recurrence = 74) referred for PET scanning as a routine part of their care. All patients had a biopsy-proven NSCLC. Nine segmentation methods were applied to each image, including manual delineation, K-means (KM), watershed, fuzzy-C-mean, region-growing, local active contour (LAC), and iterative thresholding (IT) with 40, 45, and 50% thresholds. Diverse image discretizations, both without a filter and with different wavelet decompositions, were applied to PET images. Overall, 6741 radiomic features were extracted from each image (749 radiomic features from each segmented area). Non-parametric empirical Bayes (NPEB) ComBat harmonization was used to harmonize the features. Linear Support Vector Classifier (LinearSVC) with L1 regularization For feature selection and Support Vector Machine classifier (SVM) with fivefold nested cross-validation was performed using StratifiedKFold with 'n_splits' set to 5 to predict recurrence in NSCLC patients and assess the impact of ComBat harmonization on the outcome.

RESULTS

From 749 extracted radiomic features, 206 (27%) and 389 (51%) features showed excellent reliability (ICC ≥ 0.90) against segmentation method variation before and after NPEB ComBat harmonization, respectively. Among all, 39 features demonstrated poor reliability, which declined to 10 after ComBat harmonization. The 64 fixed bin widths (without any filter) and wavelets (LLL)-based radiomic features set achieved the best performance in terms of robustness against diverse segmentation techniques before and after ComBat harmonization. The first-order and GLRLM and also first-order and NGTDM feature families showed the largest number of robust features before and after ComBat harmonization, respectively. In terms of predicting recurrence in NSCLC, our findings indicate that using ComBat harmonization can significantly enhance machine learning outcomes, particularly improving the accuracy of watershed segmentation, which initially had fewer reliable features than manual contouring. Following the application of ComBat harmonization, the majority of cases saw substantial increase in sensitivity and specificity.

CONCLUSION

Radiomic features are vulnerable to different segmentation methods. ComBat harmonization might be considered a solution to overcome the poor reliability of radiomic features.

Identifying radiogenomic associations of breast cancer based on DCE-MRI by using Siamese Neural Network with manufacturer bias normalization

BACKGROUND AND PURPOSE

The immunohistochemical test (IHC) for Human Epidermal Growth Factor Receptor 2 (HER2) and hormone receptors (HR) provides prognostic information and guides treatment for patients with invasive breast cancer. The objective of this paper is to establish a non-invasive system for identifying HER2 and HR in breast cancer using dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI).

METHODS

In light of the absence of high-performance algorithms and external validation in previously published methods, this study utilizes 3D deep features and radiomics features to represent the information of the Region of Interest (ROI). A Siamese Neural Network was employed as the classifier, with 3D deep features and radiomics features serving as the network input. To neutralize manufacturer bias, a batch effect normalization method, ComBat, was introduced. To enhance the reliability of the study, two datasets, Predict Your Therapeutic Response with Imaging and moLecular Analysis (I-SPY 1) and I-SPY 2, were incorporated. I-SPY 2 was utilized for model training and validation, while I-SPY 1 was exclusively employed for external validation. Additionally, a breast tumor segmentation network was trained to improve radiomic feature extraction.

RESULTS

The results indicate that our approach achieved an average Area Under the Curve (AUC) of 0.632, with a Standard Error of the Mean (SEM) of 0.042 for HER2 prediction in the I-SPY 2 dataset. For HR prediction, our method attained an AUC of 0.635 (SEM 0.041), surpassing other published methods in the AUC metric. Moreover, the proposed method yielded competitive results in other metrics. In external validation using the I-SPY 1 dataset, our approach achieved an AUC of 0.567 (SEM 0.032) for HR prediction and 0.563 (SEM 0.033) for HER2 prediction.

CONCLUSION

This study proposes a non-invasive system for identifying HER2 and HR in breast cancer. Although the results do not conclusively demonstrate superiority in both tasks, they indicate that the proposed method achieved good performance and is a competitive classifier compared to other reference methods. Ablation studies demonstrate that both radiomics features and deep features for the Siamese Neural Network are beneficial for the model. The introduced manufacturer bias normalization method has been shown to enhance the method's performance. Furthermore, the external validation of the method enhances the reliability of this research. Source code, pre-trained segmentation network, Radiomics and deep features, data for statistical analysis, and Supporting Information of this article are online at: https://github.com/FORRESTHUACHEN/Siamese_Neural_Network_based_Brest_cancer_Radiogenomic.

Prediction of prognosis in glioblastoma with radiomics features extracted by synthetic MRI images using cycle-consistent GAN

To propose a style transfer model for multi-contrast magnetic resonance imaging (MRI) images with a cycle-consistent generative adversarial network (CycleGAN) and evaluate the image quality and prognosis prediction performance for glioblastoma (GBM) patients from the extracted radiomics features. Style transfer models of T1 weighted MRI image (T1w) to T2 weighted MRI image (T2w) and T2w to T1w with CycleGAN were constructed using the BraTS dataset. The style transfer model was validated with the Cancer Genome Atlas Glioblastoma Multiforme (TCGA-GBM) dataset. Moreover, imaging features were extracted from real and synthesized images. These features were transformed to rad-scores by the least absolute shrinkage and selection operator (LASSO)-Cox regression. The prognosis performance was estimated by the Kaplan-Meier method. For the accuracy of the image quality of the real and synthesized MRI images, the MI, RMSE, PSNR, and SSIM were 0.991 ± 2.10 × 10 - 4 , 2.79 ± 0.16, 40.16 ± 0.38, and 0.995 ± 2.11 × 10 - 4 , for T2w, and .992 ± 2.63 × 10 - 4 , 2.49 ± 6.89 × 10 - 2 , 40.51 ± 0.22, and 0.993 ± 3.40 × 10 - 4 for T1w, respectively. The survival time had a significant difference between good and poor prognosis groups for both real and synthesized T2w (p < 0.05). However, the survival time had no significant difference between good and poor prognosis groups for both real and synthesized T1w. On the other hand, there was no significant difference between the real and synthesized T2w in both good and poor prognoses. The results of T1w were similar in the point that there was no significant difference between the real and synthesized T1w. It was found that the synthesized image could be used for prognosis prediction. The proposed prognostic model using CycleGAN could reduce the cost and time of image scanning, leading to a promotion to build the patient's outcome prediction with multi-contrast images.

Multivariate testing and effect size measures for batch effect evaluation in radiomic features

While precision medicine applications of radiomics analysis are promising, differences in image acquisition can cause "batch effects" that reduce reproducibility and affect downstream predictive analyses. Harmonization methods such as ComBat have been developed to correct these effects, but evaluation methods for quantifying batch effects are inconsistent. In this study, we propose the use of the multivariate statistical test PERMANOVA and the Robust Effect Size Index (RESI) to better quantify and characterize batch effects in radiomics data. We evaluate these methods in both simulated and real radiomics features extracted from full-field digital mammography (FFDM) data. PERMANOVA demonstrated higher power than standard univariate statistical testing, and RESI was able to interpretably quantify the effect size of site at extremely large sample sizes. These methods show promise as more powerful and interpretable methods for the detection and quantification of batch effects in radiomics studies.

The impact of the combat method on radiomics feature compensation and analysis of scanners from different manufacturers

BACKGROUND

This study investigated whether the Combat compensation method can remove the variability of radiomic features extracted from different scanners, while also examining its impact on the subsequent predictive performance of machine learning models.

MATERIALS AND METHODS

135 CT images of Credence Cartridge Radiomic phantoms were collected and screened from three scanners manufactured by Siemens, Philips, and GE. 100 radiomic features were extracted and 20 radiomic features were screened according to the Lasso regression method. The radiomic features extracted from the rubber and resin-filled regions in the cartridges were labeled into different categories for evaluating the performance of the machine learning model. Radiomics features were divided into three groups based on the different scanner manufacturers. The radiomic features were randomly divided into training and test sets with a ratio of 8:2. Five machine learning models (lasso, logistic regression, random forest, support vector machine, neural network) were employed to evaluate the impact of Combat on radiomic features. The variability among radiomic features were assessed using analysis of variance (ANOVA) and principal component analysis (PCA). Accuracy, precision, recall, and area under the receiver curve (AUC) were used as evaluation metrics for model classification.

RESULTS

The principal component and ANOVA analysis results show that the variability of different scanner manufacturers in radiomic features was removed (P˃0.05). After harmonization with the Combat algorithm, the distributions of radiomic features were aligned in terms of location and scale. The performance of machine learning models for classification improved, with the Random Forest model showing the most significant enhancement. The AUC value increased from 0.88 to 0.92.

CONCLUSIONS

The Combat algorithm has reduced variability in radiomic features from different scanners. In the phantom CT dataset, it appears that the machine learning model's classification performance may have improved after Combat harmonization. However, further investigation and validation are required to fully comprehend Combat's impact on radiomic features in medical imaging.

Noninvasive prediction of perineural invasion in intrahepatic cholangiocarcinoma by clinicoradiological features and computed tomography radiomics based on interpretable machine learning: a multicenter cohort study

BACKGROUND

Perineural invasion (PNI) of intrahepatic cholangiocarcinoma (ICC) is a strong independent risk factor for tumour recurrence and long-term patient survival. However, there is a lack of noninvasive tools for accurately predicting the PNI status. The authors develop and validate a combined model incorporating radiomics signature and clinicoradiological features based on machine learning for predicting PNI in ICC, and used the Shapley Additive explanation (SHAP) to visualize the prediction process for clinical application.

METHODS

This retrospective and prospective study included 243 patients with pathologically diagnosed ICC (training, n =136; external validation, n =81; prospective, n =26, respectively) who underwent preoperative contrast-enhanced computed tomography between January 2012 and May 2023 at three institutions (three tertiary referral centres in Guangdong Province, China). The ElasticNet was applied to select radiomics features and construct signature derived from computed tomography images, and univariate and multivariate analyses by logistic regression were used to identify the significant clinical and radiological variables with PNI. A robust combined model incorporating radiomics signature and clinicoradiological features based on machine learning was developed and the SHAP was used to visualize the prediction process. A Kaplan-Meier survival analysis was performed to compare prognostic differences between PNI-positive and PNI-negative groups and was conducted to explore the prognostic information of the combined model.

RESULTS

Among 243 patients (mean age, 61.2 years ± 11.0 (SD); 152 men and 91 women), 108 (44.4%) were diagnosed as PNI-positive. The radiomics signature was constructed by seven radiomics features, with areas under the curves of 0.792, 0.748, and 0.729 in the training, external validation, and prospective cohorts, respectively. Three significant clinicoradiological features were selected and combined with radiomics signature to construct a combined model using machine learning. The eXtreme Gradient Boosting exhibited improved accuracy and robustness (areas under the curves of 0.884, 0.831, and 0.831, respectively). Survival analysis showed the construction combined model could be used to stratify relapse-free survival (hazard ratio, 1.933; 95% CI: 1.093-3.418; P =0.021).

CONCLUSIONS

We developed and validated a robust combined model incorporating radiomics signature and clinicoradiological features based on machine learning to accurately identify the PNI statuses of ICC, and visualize the prediction process through SHAP for clinical application.

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.

Radiological artificial intelligence - predicting personalized immunotherapy outcomes in lung cancer

Personalized medicine has revolutionized approaches to treatment in the field of lung cancer by enabling therapies to be specific to each patient. However, physicians encounter an immense number of challenges in providing the optimal treatment regimen for the individual given the sheer complexity of clinical aspects such as tumor molecular profile, tumor microenvironment, expected adverse events, acquired or inherent resistance mechanisms, the development of brain metastases, the limited availability of biomarkers and the choice of combination therapy. The integration of innovative next-generation technologies such as deep learning-a subset of machine learning-and radiomics has the potential to transform the field by supporting clinical decision making in cancer treatment and the delivery of precision therapies while integrating numerous clinical considerations. In this review, we present a brief explanation of the available technologies, the benefits of using these technologies in predicting immunotherapy response in lung cancer, and the expected future challenges in the context of precision medicine.

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.

Enhancing histopathological image classification of invasive ductal carcinoma using hybrid harmonization techniques

This study aims to develop a robust pipeline for classifying invasive ductal carcinomas and benign tumors in histopathological images, addressing variability within and between centers. We specifically tackle the challenge of detecting atypical data and variability between common clusters within the same database. Our feature engineering-based pipeline comprises a feature extraction step, followed by multiple harmonization techniques to rectify intra- and inter-center batch effects resulting from image acquisition variability and diverse patient clinical characteristics. These harmonization steps facilitate the construction of more robust and efficient models. We assess the proposed pipeline's performance on two public breast cancer databases, BreaKHIS and IDCDB, utilizing recall, precision, and accuracy metrics. Our pipeline outperforms recent models, achieving 90-95% accuracy in classifying benign and malignant tumors. We demonstrate the advantage of harmonization for classifying patches from different databases. Our top model scored 94.7% for IDCDB and 95.2% for BreaKHis, surpassing existing feature engineering-based models (92.1% for IDCDB and 87.7% for BreaKHIS) and attaining comparable performance to deep learning models. The proposed feature-engineering-based pipeline effectively classifies malignant and benign tumors while addressing variability within and between centers through the incorporation of various harmonization techniques. Our findings reveal that harmonizing variabilities between patches from different batches directly impacts the learning and testing performance of classification models. This pipeline has the potential to enhance breast cancer diagnosis and treatment and may be applicable to other diseases.

A common [18F]-FDG PET radiomic signature to predict survival in patients with HPV-induced cancers

Locally advanced cervical cancer (LACC) and anal and oropharyngeal squamous cell carcinoma (ASCC and OPSCC) are mostly caused by oncogenic human papillomaviruses (HPV). In this paper, we developed machine learning (ML) models based on clinical, biological, and radiomic features extracted from pre-treatment fluorine-18-fluorodeoxyglucose positron emission tomography ([18F]-FDG PET) images to predict the survival of patients with HPV-induced cancers. For this purpose, cohorts from five institutions were used: two cohorts of patients treated for LACC including 104 patients from Gustave Roussy Campus Cancer (Center 1) and 90 patients from Leeds Teaching Hospitals NHS Trust (Center 2), two datasets of patients treated for ASCC composed of 66 patients from Institut du Cancer de Montpellier (Center 3) and 67 patients from Oslo University Hospital (Center 4), and one dataset of 45 OPSCC patients from the University Hospital of Zurich (Center 5). Radiomic features were extracted from baseline [18F]-FDG PET images. The ComBat technique was applied to mitigate intra-scanner variability. A modified consensus nested cross-validation for feature selection and hyperparameter tuning was applied on four ML models to predict progression-free survival (PFS) and overall survival (OS) using harmonized imaging features and/or clinical and biological variables as inputs. Each model was trained and optimized on Center 1 and Center 3 cohorts and tested on Center 2, Center 4, and Center 5 cohorts. The radiomic-based CoxNet model achieved C-index values of 0.75 and 0.78 for PFS and 0.76, 0.74, and 0.75 for OS on the test sets. Radiomic feature-based models had superior performance compared to the bioclinical ones, and combining radiomic and bioclinical variables did not improve the performances. Metabolic tumor volume (MTV)-based models obtained lower C-index values for a majority of the tested configurations but quite equivalent performance in terms of time-dependent AUCs (td-AUC). The results demonstrate the possibility of identifying common PET-based image signatures for predicting the response of patients with induced HPV pathology, validated on multi-center multiconstructor data.

Dual-Centre Harmonised Multimodal Positron Emission Tomography/Computed Tomography Image Radiomic Features and Machine Learning Algorithms for Non-small Cell Lung Cancer Histopathological Subtype Phenotype Decoding

AIMS

We aimed to build radiomic models for classifying non-small cell lung cancer (NSCLC) histopathological subtypes through a dual-centre dataset and comprehensively evaluate the effect of ComBat harmonisation on the performance of single- and multimodality radiomic models.

MATERIALS AND METHODS

A public dataset of NSCLC patients from two independent centres was used. Two image fusion methods, namely guided filtering-based fusion and image fusion based on visual saliency map and weighted least square optimisation, were used. Radiomic features were extracted from each scan, including first-order, texture and moment-invariant features. Subsequently, ComBat harmonisation was applied to the extracted features from computed tomography (CT), positron emission tomography (PET) and fused images to correct the centre effect. For feature selection, least absolute shrinkage and selection operator (Lasso) and recursive feature elimination (RFE) were investigated. For machine learning, logistic regression (LR), support vector machine (SVM) and AdaBoost were evaluated for classifying NSCLC subtypes. Training and evaluation of the models were carried out in a robust framework to offset plausible errors and performance was reported using area under the curve, balanced accuracy, sensitivity and specificity before and after harmonisation. N-way ANOVA was used to assess the effect of different factors on the performance of the models.

RESULTS

Support vector machine fed with selected features by recursive feature elimination from a harmonised PET feature set achieved the highest performance (area under the curve = 0.82) in classifying NSCLC histopathological subtypes. Although the performance of the models did not significantly improve for CT images after harmonisation, the performance of PET and guided filtering-based fusion feature signatures significantly improved for almost all models. Although the selection of the image modality and feature selection methods was effective on the performance of the model (ANOVA P-values <0.001), machine learning and harmonisation did not change the performance significantly (ANOVA P-values = 0.839 and 0.292, respectively).

CONCLUSION

This study confirmed the potential of radiomic analysis on PET, CT and hybrid images for histopathological classification of NSCLC subtypes.

Computed Tomography-Based Radiomics Using Tumor and Vessel Features to Assess Resectability in Cancer of the Pancreatic Head

The preoperative prediction of resectability pancreatic ductal adenocarcinoma (PDAC) is challenging. This retrospective single-center study examined tumor and vessel radiomics to predict the resectability of PDAC in chemo-naïve patients. The tumor and adjacent arteries and veins were segmented in the portal-venous phase of contrast-enhanced CT scans, and radiomic features were extracted. Features were selected via stability and collinearity testing, and least absolute shrinkage and selection operator application (LASSO). Three models, using tumor features, vessel features, and a combination of both, were trained with the training set (N = 86) to predict resectability. The results were validated with the test set (N = 15) and compared to the multidisciplinary team's (MDT) performance. The vessel-features-only model performed best, with an AUC of 0.92 and sensitivity and specificity of 97% and 73%, respectively. Test set validation showed a sensitivity and specificity of 100% and 88%, respectively. The combined model was as good as the vessel model (AUC = 0.91), whereas the tumor model showed poor performance (AUC = 0.76). The MDT's prediction reached a sensitivity and specificity of 97% and 84% for the training set and 88% and 100% for the test set, respectively. Our clinician-independent vessel-based radiomics model can aid in predicting resectability and shows performance comparable to that of the MDT. With these encouraging results, improved, automated, and generalizable models can be developed that reduce workload and can be applied in non-expert hospitals.

Optimization and validation of voxel size-related radiomics variability by combatting batch effect harmonization in pulmonary nodules: a phantom and clinical study

Background

Broad generalization of radiomics-assisted models may be impeded by concerns about variability. This study aimed to evaluate the merit of combatting batch effect (ComBat) harmonization in reducing the variability of voxel size-related radiomics in both phantom and clinical study in comparison with image resampling correction method.

Methods

A pulmonary phantom with 22 different types of nodules was scanned by computed tomography (CT) with different voxel sizes. The variability of voxel size-related radiomics features was evaluated using concordance correlation coefficient (CCC), dynamic range (DR), and intraclass correlation coefficient (ICC). ComBat and image resampling compensation methods were used to reduce variability of voxel size-related radiomics. The percentage of robust radiomics features was compared before and after optimization. Pathologically differential diagnosis of invasive adenocarcinoma (IAC) from adenocarcinoma in situ (AIS) and minimally invasive adenocarcinoma (MIA) (AIS-MIA group) was used for clinical validation in 134 patients.

Results

Before optimization, the number of excellent features in the phantom and clinical data was 26.12% and 32.31%, respectively. The excellent features were increased after image resampling and ComBat correction. For clinical optimization, the effect of the ComBat compensation method was significantly better than that of image resampling, with excellent features reaching 90.96% and poor features only amounting to 4.96%. In addition, the hierarchical clustering analysis showed that the first-order and shape features had better robustness than did texture features. In clinical validation, the area under the curve (AUC) of the testing set was 0.865 after ComBat correction.

Conclusions

The ComBat harmonization can optimize voxel size-related CT radiomics variability in pulmonary nodules more efficiently than image resampling harmonization.

Combining ultrasound radiomics, complete blood count, and serum biochemical biomarkers for diagnosing intestinal disorders in cats using machine learning

This retrospective analytical observational cohort study aimed to model and predict the classification of feline intestinal diseases from segmentations of a transverse section from small intestine ultrasound (US) image, complete blood count (CBC), and serum biochemical profile data using a variety of machine-learning approaches. In 149 cats from three institutions, images were obtained from cats with biopsy-confirmed small cell epitheliotropic lymphoma (lymphoma), inflammatory bowel disease (IBD), no pathology ("healthy"), and other conditions (warrant a biopsy for further diagnosis). CBC, blood serum chemistry, small intestinal ultrasound, and small intestinal biopsy were obtained within a 2-week interval. CBC and serum biomarkers and radiomic features were combined for modeling. Four classification schemes were investigated: (1) normal versus abnormal; (2) warranting or not warranting a biopsy; (3) lymphoma, IBD, healthy, or other conditions; and (4) lymphoma, IBD, or other conditions. Two feature selection methods were used to identify the top 3, 5, 10, and 20 features, and six machine learning models were trained. The average (95% CI) performance of models for all combinations of features, numbers of features, and types of classifiers was 0.886 (0.871-0.912) for Model 1 (normal vs. abnormal), 0.751 (0.735-0.818) for Model 2 (biopsy vs. no biopsy), 0.504 (0.450-0.556) for Model 3 (lymphoma, IBD, healthy, or other), and 0.531 (0.426-0.589), for Model 4 (lymphoma, IBD, or other). Our findings suggest model accuracies above 0.85 can be achieved in Model 1 and 2, and that including CBC and biochemistry data with US radiomics data did not significantly improve accuracy in our models.

ComBat Harmonization of Myocardial Radiomic Features Sensitive to Cardiac MRI Acquisition Parameters

Purpose

To investigate the effect of ComBat harmonization methods on the robustness of cardiac MRI-derived radiomic features to variations in imaging parameters.

Materials and Methods

This Health Insurance Portability and Accountability Act-compliant retrospective study used a publicly available data set of 11 healthy controls (mean age, 33 years ± 16 [SD]; six men) and five patients (mean age, 52 years ± 16; four men). A single midventricular short-axis section was acquired with 3-T MRI using cine balanced steady-state free precision, T1-weighted, T2-weighted, T1 mapping, and T2 mapping imaging sequences. Each sequence was acquired using baseline parameters and after variations in flip angle, spatial resolution, section thickness, and parallel imaging. Image registration was performed for all sequences at a per-individual level. Manual myocardial contouring was performed, and 1652 radiomic features per sequence were extracted using baseline and variations in imaging parameters. Radiomic feature stability to change in imaging parameters was assessed using Cohen d sensitivity. The stability of radiomic features was assessed both without and after ComBat harmonization of radiomic features. Three ComBat methods were studied: parametric, nonparametric, and Gaussian mixture model (GMM).

Results

For all sequences combined, 51.4% of features were robust to changes in imaging parameters when no ComBat method was applied. ComBat harmonization substantially increased the number of stable features to 95.1% (95% CI: 94.9, 95.3) when parametric ComBat was used and 90.9% (95% CI: 90.6, 91.2) when nonparametric ComBat was used. GMM combat resulted in only 52.6% stable features.

Conclusion

ComBat harmonization improved the stability of radiomic features to changes in imaging parameters across all cardiac MRI sequences.Keywords: Cardiac MRI, Radiomics, ComBat, Harmonization Supplemental material is available for this article. © RSNA, 2023.

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.

Harmonization of multi-site functional MRI data with dual-projection based ICA model

Modern neuroimaging studies frequently merge magnetic resonance imaging (MRI) data from multiple sites. A larger and more diverse group of participants can increase the statistical power, enhance the reliability and reproducibility of neuroimaging research, and obtain findings more representative of the general population. However, measurement biases caused by site differences in scanners represent a barrier when pooling data collected from different sites. The existence of site effects can mask biological effects and lead to spurious findings. We recently proposed a powerful denoising strategy that implements dual-projection (DP) theory based on ICA to remove site-related effects from pooled data, demonstrating the method for simulated and in vivo structural MRI data. This study investigates the use of our DP-based ICA denoising method for harmonizing functional MRI (fMRI) data collected from the Autism Brain Imaging Data Exchange II. After frequency-domain and regional homogeneity analyses, two modalities, including amplitude of low frequency fluctuation (ALFF) and regional homogeneity (ReHo), were used to validate our method. The results indicate that DP-based ICA denoising method removes unwanted site effects for both two fMRI modalities, with increases in the significance of the associations between non-imaging variables (age, sex, etc.) and fMRI measures. In conclusion, our DP method can be applied to fMRI data in multi-site studies, enabling more accurate and reliable neuroimaging research findings.

Generalization optimizing machine learning to improve CT scan radiomics and assess immune checkpoint inhibitors' response in non-small cell lung cancer: a multicenter cohort study

Background

Recent developments in artificial intelligence suggest that radiomics may represent a promising non-invasive biomarker to predict response to immune checkpoint inhibitors (ICIs). Nevertheless, validation of radiomics algorithms in independent cohorts remains a challenge due to variations in image acquisition and reconstruction. Using radiomics, we investigated the importance of scan normalization as part of a broader machine learning framework to enable model external generalizability to predict ICI response in non-small cell lung cancer (NSCLC) patients across different centers.

Methods

Radiomics features were extracted and compared from 642 advanced NSCLC patients on pre-ICI scans using established open-source PyRadiomics and a proprietary DeepRadiomics deep learning technology. The population was separated into two groups: a discovery cohort of 512 NSCLC patients from three academic centers and a validation cohort that included 130 NSCLC patients from a fourth center. We harmonized images to account for variations in reconstruction kernel, slice thicknesses, and device manufacturers. Multivariable models, evaluated using cross-validation, were used to estimate the predictive value of clinical variables, PD-L1 expression, and PyRadiomics or DeepRadiomics for progression-free survival at 6 months (PFS-6).

Results

The best prognostic factor for PFS-6, excluding radiomics features, was obtained with the combination of Clinical + PD-L1 expression (AUC = 0.66 in the discovery and 0.62 in the validation cohort). Without image harmonization, combining Clinical + PyRadiomics or DeepRadiomics delivered an AUC = 0.69 and 0.69, respectively, in the discovery cohort, but dropped to 0.57 and 0.52, in the validation cohort. This lack of generalizability was consistent with observations in principal component analysis clustered by CT scan parameters. Subsequently, image harmonization eliminated these clusters. The combination of Clinical + DeepRadiomics reached an AUC = 0.67 and 0.63 in the discovery and validation cohort, respectively. Conversely, the combination of Clinical + PyRadiomics failed generalizability validations, with AUC = 0.66 and 0.59.

Conclusion

We demonstrated that a risk prediction model combining Clinical + DeepRadiomics was generalizable following CT scan harmonization and machine learning generalization methods. These results had similar performances to routine oncology practice using Clinical + PD-L1. This study supports the strong potential of radiomics as a future non-invasive strategy to predict ICI response in advanced NSCLC.

Image harmonization: A review of statistical and deep learning methods for removing batch effects and evaluation metrics for effective harmonization

Magnetic resonance imaging and computed tomography from multiple batches (e.g. sites, scanners, datasets, etc.) are increasingly used alongside complex downstream analyses to obtain new insights into the human brain. However, significant confounding due to batch-related technical variation, called batch effects, is present in this data; direct application of downstream analyses to the data may lead to biased results. Image harmonization methods seek to remove these batch effects and enable increased generalizability and reproducibility of downstream results. In this review, we describe and categorize current approaches in statistical and deep learning harmonization methods. We also describe current evaluation metrics used to assess harmonization methods and provide a standardized framework to evaluate newly-proposed methods for effective harmonization and preservation of biological information. Finally, we provide recommendations to end-users to advocate for more effective use of current methods and to methodologists to direct future efforts and accelerate development of the field.

Machine Learning-Based Detection of Parkinson's Disease From Resting-State EEG: A Multi-Center Study

Resting-state EEG (rs-EEG) has been demonstrated to aid in Parkinson's disease (PD) diagnosis. In particular, the power spectral density (PSD) of low-frequency bands (δ and θ) and high-frequency bands (α and β) has been shown to be significantly different in patients with PD as compared to subjects without PD (non-PD). However, rs-EEG feature extraction and the interpretation thereof can be time-intensive and prone to examiner variability. Machine learning (ML) has the potential to automatize the analysis of rs-EEG recordings and provides a supportive tool for clinicians to ease their workload. In this work, we use rs-EEG recordings of 84 PD and 85 non-PD subjects pooled from four datasets obtained at different centers. We propose an end-to-end pipeline consisting of preprocessing, extraction of PSD features from clinically-validated frequency bands, and feature selection. Following, we assess the classification ability of the features via ML algorithms to stratify between PD and non-PD subjects. Further, we evaluate the effect of feature harmonization, given the multi-center nature of the datasets. Our validation results show, on average, an improvement in PD detection ability (69.6% vs. 75.5% accuracy) by logistic regression when harmonizing the features and performing univariate feature selection (k = 202 features). Our final results show an average global accuracy of 72.2% with balanced accuracy results for all the centers included in the study: 60.6%, 68.7%, 77.7%, and 82.2%, respectively.Clinical relevance- We present an end-to-end pipeline to extract clinically relevant features from rs-EEG recordings that can facilitate the analysis and detection of PD. Further, we provide an ML system that shows a good performance in detecting PD, even in the presence of centers with different acquisition protocols. Our results show the relevance of harmonizing features and provide a good starting point for future studies focusing on rs-EEG analysis and multi-center data.

Optimal 18F-FDG PET/CT radiomics model development for predicting EGFR mutation status and prognosis in lung adenocarcinoma: a multicentric study

Purpose

To develop and interpret optimal predictive models to identify epidermal growth factor receptor (EGFR) mutation status and subtypes in patients with lung adenocarcinoma based on multicentric ¹⁸F-FDG PET/CT data, and further construct a prognostic model to predict their clinical outcome.

Methods

The ¹⁸F-FDG PET/CT imaging and clinical characters of 767 patients with lung adenocarcinoma from 4 cohorts were collected. Seventy-six radiomics candidates using cross-combination method to identity EGFR mutation status and subtypes were built. Further, Shapley additive explanations and local interpretable model-agnostic explanations were used for optimal models' interpretation. Moreover, in order to predict the overall survival, a multivariate Cox proportional hazard model based on handcrafted radiomics features and clinical characteristics was constructed. The predictive performance and clinical net benefit of the models were evaluated via area under receiver operating characteristic (AUC), C-index and decision curve analysis.

Results

Among the 76 radiomics candidates, light gradient boosting machine classifier (LGBM) combined with recursive feature elimination wrapped LGBM feature selection method achieved best performance in predicting EGFR mutation status (AUC reached 0.80, 0.61, 0.71 in the internal test cohort and two external test cohorts, respectively). And extreme gradient boosting classifier combined with support vector machine feature selection method achieved best performance in predicting EGFR subtypes (AUC reached 0.76, 0.63, 0.61 in the internal test cohort and two external test cohorts, respectively). The C-index of the Cox proportional hazard model achieved 0.863.

Conclusions

The integration of cross-combination method and the external validation from multi-center data achieved a good prediction and generalization performance in predicting EGFR mutation status and its subtypes. The combination of handcrafted radiomics features and clinical factors achieved good performance in predicting prognosis. With the urgent needs of multicentric ¹⁸F-FDG PET/CT trails, robust and explainable radiomics models have great potential in decision making and prognosis prediction of lung adenocarcinoma.

Better Together: Data Harmonization and Cross-Study Analysis of Abdominal MRI Data From UK Biobank and the German National Cohort

OBJECTIVES

The UK Biobank (UKBB) and German National Cohort (NAKO) are among the largest cohort studies, capturing a wide range of health-related data from the general population, including comprehensive magnetic resonance imaging (MRI) examinations. The purpose of this study was to demonstrate how MRI data from these large-scale studies can be jointly analyzed and to derive comprehensive quantitative image-based phenotypes across the general adult population.

MATERIALS AND METHODS

Image-derived features of abdominal organs (volumes of liver, spleen, kidneys, and pancreas; volumes of kidney hilum adipose tissue; and fat fractions of liver and pancreas) were extracted from T1-weighted Dixon MRI data of 17,996 participants of UKBB and NAKO based on quality-controlled deep learning generated organ segmentations. To enable valid cross-study analysis, we first analyzed the data generating process using methods of causal discovery. We subsequently harmonized data from UKBB and NAKO using the ComBat approach for batch effect correction. We finally performed quantile regression on harmonized data across studies providing quantitative models for the variation of image-derived features stratified for sex and dependent on age, height, and weight.

RESULTS

Data from 8791 UKBB participants (49.9% female; age, 63 ± 7.5 years) and 9205 NAKO participants (49.1% female, age: 51.8 ± 11.4 years) were analyzed. Analysis of the data generating process revealed direct effects of age, sex, height, weight, and the data source (UKBB vs NAKO) on image-derived features. Correction of data source-related effects resulted in markedly improved alignment of image-derived features between UKBB and NAKO. Cross-study analysis on harmonized data revealed comprehensive quantitative models for the phenotypic variation of abdominal organs across the general adult population.

CONCLUSIONS

Cross-study analysis of MRI data from UKBB and NAKO as proposed in this work can be helpful for future joint data analyses across cohorts linking genetic, environmental, and behavioral risk factors to MRI-derived phenotypes and provide reference values for clinical diagnostics.

Normalization Strategies in Multi-Center Radiomics Abdominal MRI: Systematic Review and Meta-Analyses

Goal: Artificial intelligence applied to medical image analysis has been extensively used to develop non-invasive diagnostic and prognostic signatures. However, these imaging biomarkers should be largely validated on multi-center datasets to prove their robustness before they can be introduced into clinical practice. The main challenge is represented by the great and unavoidable image variability which is usually addressed using different pre-processing techniques including spatial, intensity and feature normalization. The purpose of this study is to systematically summarize normalization methods and to evaluate their correlation with the radiomics model performances through meta-analyses. This review is carried out according to the PRISMA statement: 4777 papers were collected, but only 74 were included. Two meta-analyses were carried out according to two clinical aims: characterization and prediction of response. Findings of this review demonstrated that there are some commonly used normalization approaches, but not a commonly agreed pipeline that can allow to improve performance and to bridge the gap between bench and bedside.

Cross-scanner harmonization methods for structural MRI may need further work: A comparison study

The clinical usefulness MRI biomarkers for aging and dementia studies relies on precise brain morphological measurements; however, scanner and/or protocol variations may introduce noise or bias. One approach to address this is post-acquisition scan harmonization. In this work, we evaluate deep learning (neural style transfer, CycleGAN and CGAN), histogram matching, and statistical (ComBat and LongComBat) methods. Participants who had been scanned on both GE and Siemens scanners (cross-sectional participants, known as Crossover (n = 113), and longitudinally scanned participants on both scanners (n = 454)) were used. The goal was to match GE MPRAGE (T1-weighted) scans to Siemens improved resolution MPRAGE scans. Harmonization was performed on raw native and preprocessed (resampled, affine transformed to template space) scans. Cortical thicknesses were measured using FreeSurfer (v.7.1.1). Distributions were checked using Kolmogorov-Smirnov tests. Intra-class correlation (ICC) was used to assess the degree of agreement in the Crossover datasets and annualized percent change in cortical thickness was calculated to evaluate the Longitudinal datasets. Prior to harmonization, the least agreement was found at the frontal pole (ICC = 0.72) for the raw native scans, and at caudal anterior cingulate (0.76) and frontal pole (0.54) for the preprocessed scans. Harmonization with NST, CycleGAN, and HM improved the ICCs of the preprocessed scans at the caudal anterior cingulate (>0.81) and frontal poles (>0.67). In the Longitudinal raw native scans, over- and under-estimations of cortical thickness were observed due to the changing of the scanners. ComBat matched the cortical thickness distributions throughout but was not able to increase the ICCs or remove the effects of scanner changeover in the Longitudinal datasets. CycleGAN and NST performed slightly better to address the cortical thickness variations between scanner change. However, none of the methods succeeded in harmonizing the Longitudinal dataset. CGAN was the worst performer for both datasets. In conclusion, the performance of the methods was overall similar and region dependent. Future research is needed to improve the existing approaches since none of them outperformed each other in terms of harmonizing the datasets at all ROIs. The findings of this study establish framework for future research into the scan harmonization problem.

Prognostic analysis of curatively resected pancreatic cancer using harmonized positron emission tomography radiomic features

BACKGROUND

Texture features reflecting tumour heterogeneity enable us to investigate prognostic factors. The R package ComBat can harmonize the quantitative texture features among several positron emission tomography (PET) scanners. We aimed to identify prognostic factors among harmonized PET radiomic features and clinical information from pancreatic cancer patients who underwent curative surgery.

METHODS

Fifty-eight patients underwent preoperative enhanced dynamic computed tomography (CT) scanning and fluorodeoxyglucose PET/CT using four PET scanners. Using LIFEx software, we measured PET radiomic parameters including texture features with higher order and harmonized these PET parameters. For progression-free survival (PFS) and overall survival (OS), we evaluated clinical information, including age, TNM stage, and neural invasion, and the harmonized PET radiomic features based on univariate Cox proportional hazard regression. Next, we analysed the prognostic indices by multivariate Cox proportional hazard regression (1) by using either significant (p < 0.05) or borderline significant (p = 0.05-0.10) indices in the univariate analysis (first multivariate analysis) or (2) by using the selected features with random forest algorithms (second multivariate analysis). Finally, we checked these multivariate results by log-rank test.

RESULTS

Regarding the first multivariate analysis for PFS after univariate analysis, age was the significant prognostic factor (p = 0.020), and MTV and GLCM contrast were borderline significant (p = 0.051 and 0.075, respectively). Regarding the first multivariate analysis of OS, neural invasion, Shape sphericity and GLZLM LZLGE were significant (p = 0.019, 0.042 and 0.0076). In the second multivariate analysis, only MTV was significant (p = 0.046) for PFS, whereas GLZLM LZLGE was significant (p = 0.047), and Shape sphericity was borderline significant (p = 0.088) for OS. In the log-rank test, age, MTV and GLCM contrast were borderline significant for PFS (p = 0.08, 0.06 and 0.07, respectively), whereas neural invasion and Shape sphericity were significant (p = 0.03 and 0.04, respectively), and GLZLM LZLGE was borderline significant for OS (p = 0.08).

CONCLUSIONS

Other than the clinical factors, MTV and GLCM contrast for PFS and Shape sphericity and GLZLM LZLGE for OS may be prognostic PET parameters. A prospective multicentre study with a larger sample size may be warranted.

Machine learning for normal tissue complication probability prediction: Predictive power with versatility and easy implementation

Background and purpose

A popular Normal tissue Complication (NTCP) model deployed to predict radiotherapy (RT) toxicity is the Lyman-Burman Kutcher (LKB) model of tissue complication. Despite the LKB model's popularity, it can suffer from numerical instability and considers only the generalized mean dose (GMD) to an organ. Machine learning (ML) algorithms can potentially offer superior predictive power of the LKB model, and with fewer drawbacks. Here we examine the numerical characteristics and predictive power of the LKB model and compare these with those of ML.

Materials and methods

Both an LKB model and ML models were used to predict G2 Xerostomia on patients following RT for head and neck cancer, using the dose volume histogram of parotid glands as the input feature. Model speed, convergence characteristics and predictive power was evaluated on an independent training set.

Results

We found that only global optimization algorithms could guarantee a convergent and predictive LKB model. At the same time our results showed that ML models remained unconditionally convergent and predictive, while staying robust to gradient descent optimization. ML models outperform LKB in Brier score and accuracy but compare to LKB in ROC-AUC.

Conclusion

We have demonstrated that ML models can quantify NTCP better than or as well as LKB models, even for a toxicity that the LKB model is particularly well suited to predict. ML models can offer this performance while offering fundamental advantages in model convergence, speed, and flexibility, and so could offer an alternative to the LKB model that could potentially be used in clinical RT planning decisions.

MR Intensity Normalization Methods Impact Sequence Specific Radiomics Prognostic Model Performance in Primary and Recurrent High-Grade Glioma

PURPOSE

This study investigates the impact of different intensity normalization (IN) methods on the overall survival (OS) radiomics models' performance of MR sequences in primary (pHGG) and recurrent high-grade glioma (rHGG).

METHODS

MR scans acquired before radiotherapy were retrieved from two independent cohorts (rHGG C1: 197, pHGG C2: 141) from multiple scanners (15, 14). The sequences are T1 weighted (w), contrast-enhanced T1w (T1wce), T2w, and T2w-FLAIR. Sequence-specific significant features (SF) associated with OS, extracted from the tumour volume, were derived after applying 15 different IN methods. Survival analyses were conducted using Cox proportional hazard (CPH) and Poisson regression (POI) models. A ranking score was assigned based on the 10-fold cross-validated (CV) concordance index (C-I), mean square error (MSE), and the Akaike information criterion (AICs), to evaluate the methods' performance.

RESULTS

Scatter plots of the 10-CV C-I and MSE against the AIC showed an impact on the survival predictions between the IN methods and MR sequences (C1/C2 C-I range: 0.62-0.71/0.61-0.72, MSE range: 0.20-0.42/0.13-0.22). White stripe showed stable results for T1wce (C1/C2 C-I: 0.71/0.65, MSE: 0.21/0.14). Combat (0.68/0.62, 0.22/0.15) and histogram matching (HM, 0.67/0.64, 0.22/0.15) showed consistent prediction results for T2w models. They were also the top-performing methods for T1w in C2 (Combat: 0.67, 0.13; HM: 0.67, 0.13); however, only HM achieved high predictions in C1 (0.66, 0.22). After eliminating IN impacted SF using Spearman's rank-order correlation coefficient, a mean decrease in the C-I and MSE of 0.05 and 0.03 was observed in all four sequences.

CONCLUSION

The IN method impacted the predictive power of survival models; thus, performance is sequence-dependent.

Multisequence magnetic resonance imaging-based radiomics models for the prediction of microsatellite instability in endometrial cancer

PURPOSE

To evaluate the performance of multisequence magnetic resonance imaging (MRI)-based radiomics models in the assessment of microsatellite instability (MSI) status in endometrial cancer (EC).

MATERIALS AND METHODS

This retrospective multicentre study included 338 EC patients with available MSI status and preoperative MRI scans, divided into training (37 MSI, 123 microsatellite stability [MSS]), internal validation (15 MSI, 52 MSS), and external validation cohorts (30 MSI, 81 MSS). Radiomics features were extracted from T2-weighted images, diffusion-weighted images, and contrast-enhanced T1-weighted images. The ComBat harmonisation method was applied to remove intrascanner variability. The Boruta wrapper algorithm was used for key feature selection. Three classification algorithms, logistic regression (LR), random forest (RF), and support vector machine (SVM), were applied to build the radiomics models. The area under the receiver operating characteristic curve (AUC) was calculated to compare the diagnostic performance of the models. Decision curve analysis (DCA) was conducted to determine the clinical usefulness of the models.

RESULTS

Among the 1980 features, Boruta finally selected nine radiomics features. A higher MSI prediction performance was achieved after running the ComBat harmonisation method. The SVM algorithm had the best performance, with AUCs of 0.921, 0.903, and 0.937 in the training, internal validation, and external validation cohorts, respectively. The DCA results showed that the SVM algorithm achieved higher net benefits than the other classifiers over a threshold range of 0.581-0.783.

CONCLUSION

The multisequence MRI-based radiomics models showed promise in preoperatively predicting the MSI status in EC in this multicentre setting.

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.

Joint EANM/SNMMI guideline on radiomics in nuclear medicine : Jointly supported by the EANM Physics Committee and the SNMMI Physics, Instrumentation and Data Sciences Council

PURPOSE

The purpose of this guideline is to provide comprehensive information on best practices for robust radiomics analyses for both hand-crafted and deep learning-based approaches.

METHODS

In a cooperative effort between the EANM and SNMMI, we agreed upon current best practices and recommendations for relevant aspects of radiomics analyses, including study design, quality assurance, data collection, impact of acquisition and reconstruction, detection and segmentation, feature standardization and implementation, as well as appropriate modelling schemes, model evaluation, and interpretation. We also offer an outlook for future perspectives.

CONCLUSION

Radiomics is a very quickly evolving field of research. The present guideline focused on established findings as well as recommendations based on the state of the art. Though this guideline recognizes both hand-crafted and deep learning-based radiomics approaches, it primarily focuses on the former as this field is more mature. This guideline will be updated once more studies and results have contributed to improved consensus regarding the application of deep learning methods for radiomics. Although methodological recommendations in the present document are valid for most medical image modalities, we focus here on nuclear medicine, and specific recommendations when necessary are made for PET/CT, PET/MR, and quantitative SPECT.

Combination of pre-treatment dynamic [18F]FET PET radiomics and conventional clinical parameters for the survival stratification in patients with IDH-wildtype glioblastoma

PURPOSE

The aim of this study was to build and evaluate a prediction model which incorporates clinical parameters and radiomic features extracted from static as well as dynamic [18F]FET PET for the survival stratification in patients with newly diagnosed IDH-wildtype glioblastoma.

METHODS

A total of 141 patients with newly diagnosed IDH-wildtype glioblastoma and dynamic [18F]FET PET prior to surgical intervention were included. Patients with a survival time ≤ 12 months were classified as short-term survivors. First order, shape, and texture radiomic features were extracted from pre-treatment static (tumor-to-background ratio; TBR) and dynamic (time-to-peak; TTP) images, respectively, and randomly divided into a training (n = 99) and a testing cohort (n = 42). After feature normalization, recursive feature elimination was applied for feature selection using 5-fold cross-validation on the training cohort, and a machine learning model was constructed to compare radiomic models and combined clinical-radiomic models with selected radiomic features and clinical parameters. The area under the ROC curve (AUC), accuracy, sensitivity, specificity, and positive and negative predictive values were calculated to assess the predictive performance for identifying short-term survivors in both the training and testing cohort.

RESULTS

A combined clinical-radiomic model comprising six clinical parameters and six selected dynamic radiomic features achieved highest predictability of short-term survival with an AUC of 0.74 (95% confidence interval, 0.60-0.88) in the independent testing cohort.

CONCLUSIONS

This study successfully built and evaluated prediction models using [18F]FET PET-based radiomic features and clinical parameters for the individualized assessment of short-term survival in patients with a newly diagnosed IDH-wildtype glioblastoma. The combination of both clinical parameters and dynamic [18F]FET PET-based radiomic features reached highest accuracy in identifying patients at risk. Although the achieved accuracy level remained moderate, our data shows that the integration of dynamic [18F]FET PET radiomic data into clinical prediction models may improve patient stratification beyond established prognostic markers.

A Framework of Analysis to Facilitate the Harmonization of Multicenter Radiomic Features in Prostate Cancer

Pooling radiomic features coming from different centers in a statistical framework is challenging due to the variability in scanner models, acquisition protocols, and reconstruction settings. To remove technical variability, commonly called batch effects, different statistical harmonization strategies have been widely used in genomics but less considered in radiomics. The aim of this work was to develop a framework of analysis to facilitate the harmonization of multicenter radiomic features extracted from prostate T2-weighted magnetic resonance imaging (MRI) and to improve the power of radiomics for prostate cancer (PCa) management in order to develop robust non-invasive biomarkers translating into clinical practice. To remove technical variability and correct for batch effects, we investigated four different statistical methods (ComBat, SVA, Arsynseq, and mixed effect). The proposed approaches were evaluated using a dataset of 210 prostate cancer (PCa) patients from two centers. The impacts of the different statistical approaches were evaluated by principal component analysis and classification methods (LogitBoost, random forest, K-nearest neighbors, and decision tree). The ComBat method outperformed all other methods by achieving 70% accuracy and 78% AUC with the random forest method to automatically classify patients affected by PCa. The proposed statistical framework enabled us to define and develop a standardized pipeline of analysis to harmonize multicenter T2W radiomic features, yielding great promise to support PCa clinical practice.

Predicting subacute hematological toxicity of 177Lu-DOTATATE therapy using healthy organs' uptake on post-treatment quantitative SPECT: A pilot study

The aim is to investigate the usefulness of 177Lu-DOTA-0-Tyr3-Octreotate (DOTATATE) healthy organs' (spleen, kidneys, bone marrow) standard uptake value for the prediction of subacute hematological toxicity in patients undergoing 177Lu-DOTATATE treatment. All patients referred from January 2021 to May 2022 for 177Lu-DOTATATE treatment were retrospectively screened. For each treatment session, baseline clinical data including age, sex, weight, delay between 177Lu-DOTATATE treatment and last cold somatostatin analogue intake were collected. Mean standardized uptake value (SUVmean) of healthy organs was measured and analyzed by generalized linear mixed effect models. Outcomes (significant decrease of platelets, hemoglobin levels and neutrophils) were assessed 1 month later, considering their within-subject biological coefficient of variation, published by the European Federation of Clinical Chemistry and Laboratory Medicine. A total of 9 patients (33 treatment sessions) were included. No predictive factors were identified for platelet and neutrophil decrease. Splenic SUVmean was found to be a significant predictor of hemoglobin levels decrease. Using an optimal threshold of ≥6.22, derived sensitivity and specificity to predict hemoglobin decrease were 85.7% [46.4; 99.0] and 76.9% [57.5; 89.2] respectively with an accuracy of 82.4%. Although not significantly predictive of hematological toxicity, bone marrow SUVmean and renal SUVmean were correlated with splenic SUVmean. Quantitative single photon emission computed tomography and healthy organs analysis might help to foresee hematological subacute toxicity in patients undergoing 177Lu-DOTATATE treatment and improve patient management.

Performance sensitivity analysis of brain metastasis stereotactic radiosurgery outcome prediction using MRI radiomics

Recent studies have used T1w contrast-enhanced (T1w-CE) magnetic resonance imaging (MRI) radiomic features and machine learning to predict post-stereotactic radiosurgery (SRS) brain metastasis (BM) progression, but have not examined the effects of combining clinical and radiomic features, BM primary cancer, BM volume effects, and using multiple scanner models. To investigate these effects, a dataset of n = 123 BMs from 99 SRS patients with 12 clinical features, 107 pre-treatment T1w-CE radiomic features, and BM progression determined by follow-up MRI was used with a random decision forest model and 250 bootstrapped repetitions. Repeat experiments assessed the relative accuracy across primary cancer sites, BM volume groups, and scanner model pairings. Correction for accuracy imbalances across volume groups was investigated by removing volume-correlated features. We found that using clinical and radiomic features together produced the most accurate model with a bootstrap-corrected area under the receiver operating characteristic curve of 0.77. Accuracy also varied by primary cancer site, BM volume, and scanner model pairings. The effect of BM volume was eliminated by removing features at a volume-correlation coefficient threshold of 0.25. These results show that feature type, primary cancer, volume, and scanner model are all critical factors in the accuracy of radiomics-based prognostic models for BM SRS that must be characterised and controlled for before clinical translation.

A methodological framework for AI-assisted diagnosis of active aortitis using radiomic analysis of FDG PET-CT images: Initial analysis

BACKGROUND

The aim of this study was to explore the feasibility of assisted diagnosis of active (peri-)aortitis using radiomic imaging biomarkers derived from [18F]-Fluorodeoxyglucose Positron Emission Tomography-Computed Tomography (FDG PET-CT) images.

METHODS

The aorta was manually segmented on FDG PET-CT in 50 patients with aortitis and 25 controls. Radiomic features (RF) (n = 107), including SUV (Standardized Uptake Value) metrics, were extracted from the segmented data and harmonized using the ComBat technique. Individual RFs and groups of RFs (i.e., signatures) were used as input in Machine Learning classifiers. The diagnostic utility of these classifiers was evaluated with area under the receiver operating characteristic curve (AUC) and accuracy using the clinical diagnosis as the ground truth.

RESULTS

Several RFs had high accuracy, 84% to 86%, and AUC scores 0.83 to 0.97 when used individually. Radiomic signatures performed similarly, AUC 0.80 to 1.00.

CONCLUSION

A methodological framework for a radiomic-based approach to support diagnosis of aortitis was outlined. Selected RFs, individually or in combination, showed similar performance to the current standard of qualitative assessment in terms of AUC for identifying active aortitis. This framework could support development of a clinical decision-making tool for a more objective and standardized assessment of aortitis.

Harmonization Strategies in Multicenter MRI-Based Radiomics

Radiomics analysis is a powerful tool aiming to provide diagnostic and prognostic patient information directly from images that are decoded into handcrafted features, comprising descriptors of shape, size and textural patterns. Although radiomics is gaining momentum since it holds great promise for accelerating digital diagnostics, it is susceptible to bias and variation due to numerous inter-patient factors (e.g., patient age and gender) as well as inter-scanner ones (different protocol acquisition depending on the scanner center). A variety of image and feature based harmonization methods has been developed to compensate for these effects; however, to the best of our knowledge, none of these techniques has been established as the most effective in the analysis pipeline so far. To this end, this review provides an overview of the challenges in optimizing radiomics analysis, and a concise summary of the most relevant harmonization techniques, aiming to provide a thorough guide to the radiomics harmonization process.

A Framework for Harmonization of Radiomics Data for Multicenter Studies and Clinical Trials

PURPOSE

Variability in computed tomography images intrinsic to individual scanners limits the application of radiomics in clinical and research settings. The development of reproducible and generalizable radiomics-based models to assess lesions requires harmonization of data. The purpose of this study was to develop, test, and analyze the efficacy of a radiomics data harmonization model.

MATERIALS AND METHODS

Radiomic features from biopsy-proven untreated hepatic metastasis (N = 380) acquired from 167 unique patients with pancreatic, colon, and breast cancers were analyzed. Radiomic features from volume-match 551 samples of normal liver tissue and 188 hepatic cysts were included as references. A novel linear mixed effect model was used to identify effects associated with lesion size, tissue type, and scanner model. Six separate machine learning models were then used to test the effectiveness of radiomic feature harmonization using multivariate analysis.

RESULTS

Proposed model identifies and removes scanner-associated effects while preserving cancer-specific functional dependence of radiomic features on the tumor size. Data harmonization improves the performance of classification models by reducing the scanner-associated variability. For example, the multiclass logistic regression model, LogitBoost, demonstrated the improvement in sensitivity in the range from 15% to 40% for each type of liver metastasis, whereas the overall model accuracy and the kappa coefficient increased by 5% and 8% accordingly.

CONCLUSION

The model removed scanner-associated effects while preserving cancer-specific functional dependence of radiomic features.

Site effects how-to and when: An overview of retrospective techniques to accommodate site effects in multi-site neuroimaging analyses

Site differences, or systematic differences in feature distributions across multiple data-acquisition sites, are a known source of heterogeneity that may adversely affect large-scale meta- and mega-analyses of independently collected neuroimaging data. They influence nearly all multi-site imaging modalities and biomarkers, and methods to compensate for them can improve reliability and generalizability in the analysis of genetics, omics, and clinical data. The origins of statistical site effects are complex and involve both technical differences (scanner vendor, head coil, acquisition parameters, imaging processing) and differences in sample characteristics (inclusion/exclusion criteria, sample size, ancestry) between sites. In an age of expanding international consortium research, there is a growing need to disentangle technical site effects from sample characteristics of interest. Numerous statistical and machine learning methods have been developed to control for, model, or attenuate site effects - yet to date, no comprehensive review has discussed the benefits and drawbacks of each for different use cases. Here, we provide an overview of the different existing statistical and machine learning methods developed to remove unwanted site effects from independently collected neuroimaging samples. We focus on linear mixed effect models, the ComBat technique and its variants, adjustments based on image quality metrics, normative modeling, and deep learning approaches such as generative adversarial networks. For each method, we outline the statistical foundation and summarize strengths and weaknesses, including their assumptions and conditions of use. We provide information on software availability and comment on the ease of use and the applicability of these methods to different types of data. We discuss validation and comparative reports, mention caveats and provide guidance on when to use each method, depending on context and specific research questions.

The impact of harmonization on radiomic features in Parkinson’s disease and healthy controls: A multicenter study

Radiomics is a challenging development area in imaging field that is greatly capturing interest of radiologists and neuroscientists. However, radiomics features show a strong non-biological variability determined by different facilities and imaging protocols, limiting the reproducibility and generalizability of analysis frameworks. Our study aimed to investigate the usefulness of harmonization to reduce site-effects on radiomics features over specific brain regions. We selected T1-weighted magnetic resonance imaging (MRI) by using the MRI dataset Parkinson’s Progression Markers Initiative (PPMI) from different sites with healthy controls (HC) and Parkinson’s disease (PD) patients. First, the investigation of radiomics measure discrepancies were assessed on healthy brain regions-of-interest (ROIs) via a classification pipeline based on LASSO feature selection and support vector machine (SVM) model. Then, a ComBat-based harmonization approach was applied to correct site-effects. Finally, a validation step on PD subjects evaluated diagnostic accuracy before and after harmonization of radiomics data. Results on healthy subjects demonstrated a dependence from site-effects that could be corrected with ComBat harmonization. LASSO regressor after harmonization was unable to select any feature to distinguish controls by site. Moreover, harmonized radiomics features achieved an area under the receiving operating characteristic curve (AUC) of 0.77 (compared to AUC of 0.71 for raw radiomics measures) in distinguish Parkinson’s patients from HC. We found a not-negligible site-effect studying radiomics of HC pre- and post-harmonization of features. Our validation study on PD patients demonstrated a significant influence of non-biological noise source in diagnostic performances. Finally, harmonization of multicenter radiomic data represent a necessary step to make analysis pipelines reliable and replicable for multisite neuroimaging studies.

Reliability of multi-site UK Biobank MRI brain phenotypes for the assessment of neuropsychiatric complications of SARS-CoV-2 infection: The COVID-CNS travelling heads study

INTRODUCTION

Magnetic resonance imaging (MRI) of the brain could be a key diagnostic and research tool for understanding the neuropsychiatric complications of COVID-19. For maximum impact, multi-modal MRI protocols will be needed to measure the effects of SARS-CoV-2 infection on the brain by diverse potentially pathogenic mechanisms, and with high reliability across multiple sites and scanner manufacturers. Here we describe the development of such a protocol, based upon the UK Biobank, and its validation with a travelling heads study. A multi-modal brain MRI protocol comprising sequences for T1-weighted MRI, T2-FLAIR, diffusion MRI (dMRI), resting-state functional MRI (fMRI), susceptibility-weighted imaging (swMRI), and arterial spin labelling (ASL), was defined in close approximation to prior UK Biobank (UKB) and C-MORE protocols for Siemens 3T systems. We iteratively defined a comparable set of sequences for General Electric (GE) 3T systems. To assess multi-site feasibility and between-site variability of this protocol, N = 8 healthy participants were each scanned at 4 UK sites: 3 using Siemens PRISMA scanners (Cambridge, Liverpool, Oxford) and 1 using a GE scanner (King's College London). Over 2,000 Imaging Derived Phenotypes (IDPs), measuring both data quality and regional image properties of interest, were automatically estimated by customised UKB image processing pipelines (S2 File). Components of variance and intra-class correlations (ICCs) were estimated for each IDP by linear mixed effects models and benchmarked by comparison to repeated measurements of the same IDPs from UKB participants. Intra-class correlations for many IDPs indicated good-to-excellent between-site reliability. Considering only data from the Siemens sites, between-site reliability generally matched the high levels of test-retest reliability of the same IDPs estimated in repeated, within-site, within-subject scans from UK Biobank. Inclusion of the GE site resulted in good-to-excellent reliability for many IDPs, although there were significant between-site differences in mean and scaling, and reduced ICCs, for some classes of IDP, especially T1 contrast and some dMRI-derived measures. We also identified high reliability of quantitative susceptibility mapping (QSM) IDPs derived from swMRI images, multi-network ICA-based IDPs from resting-state fMRI, and olfactory bulb structure IDPs from T1, T2-FLAIR and dMRI data.

CONCLUSION

These results give confidence that large, multi-site MRI datasets can be collected reliably at different sites across the diverse range of MRI modalities and IDPs that could be mechanistically informative in COVID brain research. We discuss limitations of the study and strategies for further harmonisation of data collected from sites using scanners supplied by different manufacturers. These acquisition and analysis protocols are now in use for MRI assessments of post-COVID patients (N = 700) as part of the ongoing COVID-CNS study.

Striatal D1 and D2 receptor availability are selectively associated with eye-blink rates after methylphenidate treatment

Eye-blink rate has been proposed as a biomarker of the brain dopamine system, however, findings have not been consistent. This study assessed the relationship between blink rates, measured after oral placebo) (PL) and after a challenge with oral methylphenidate (MP; 60 mg) and striatal D1 receptor (D1R) (measured at baseline) and D2 receptor (D2R) availability (measured after PL and after MP) in healthy participants. PET measures of baseline D1R ([¹¹C]NNC112) (BL-D1R) and D2R availability ([¹¹C]raclopride) after PL (PL-D2R) and after MP (MP-D2R) were quantified in the striatum as non-displaceable binding potential. MP reduced the number of blinks and increased the time participants kept their eyes open. Correlations with dopamine receptors were only significant for the eye blink measures obtained after MP; being positive for BL-D1R in putamen and MP-D2R in caudate (PL-D2R were not significant). MP-induced changes in blink rates (PL minus MP) were negatively correlated with BL-D1R in caudate and putamen. Our findings suggest that eye blink measures obtained while stressing the dopamine system might provide a more sensitive behavioral biomarker of striatal D1R or D2R in healthy volunteers than that obtained at baseline or after placebo.

PET imaging in human participants revealed that D1 and D2 dopamine receptor availability was associated with eye-blink rates following treatment with oral methylphenidate, but not a placebo.

Transfer learning approach based on computed tomography images for predicting late xerostomia after radiotherapy in patients with oropharyngeal cancer

Background and purpose

Although the latest breakthroughs in radiotherapy (RT) techniques have led to a decrease in adverse event rates, these techniques are still associated with substantial toxicity, including xerostomia. Imaging biomarkers could be useful to predict the toxicity risk related to each individual patient. Our preliminary work aims to develop a radiomic-based support tool exploiting pre-treatment CT images to predict late xerostomia risk in 3 months after RT in patients with oropharyngeal cancer (OPC).

Materials and methods

We performed a multicenter data collection. We enrolled 61 patients referred to three care centers in Apulia, Italy, out of which 22 patients experienced at least mild xerostomia 3 months after the end of the RT cycle. Pre-treatment CT images, clinical and dose features, and alcohol-smoking habits were collected. We proposed a transfer learning approach to extract quantitative imaging features from CT images by means of a pre-trained convolutional neural network (CNN) architecture. An optimal feature subset was then identified to train an SVM classifier. To evaluate the robustness of the proposed model with respect to different manual contouring practices on CTs, we repeated the same image analysis pipeline on “fake” parotid contours.

Results

The best performances were achieved by the model exploiting the radiomic features alone. On the independent test, the model reached median AUC, accuracy, sensitivity, and specificity values of 81.17, 83.33, 71.43, and 90.91%, respectively. The model was robust with respect to diverse manual parotid contouring procedures.

Conclusion

Radiomic analysis could help to develop a valid support tool for clinicians in planning radiotherapy treatment, by providing a risk score of the toxicity development for each individual patient, thus improving the quality of life of the same patient, without compromising patient care.

Capecitabine—A “Permanent Mission” in Head and Neck Cancers “War Council”?

Capecitabine, an oral pro-drug that is metabolized to 5-FU, has been used in clinical practice for more than 20 years, being part of the therapeutic standard for digestive and breast cancers. The use of capecitabine has been evaluated in many trials including cases diagnosed in recurrent or metastatic settings. Induction regimens or a combination with radiation therapy were evaluated in head and neck cancers, but 5-FU still remained the fluoropyrimidine used as a part of the current therapeutic standard. Quantifications of levels or ratios for enzymes are involved in the capecitabine metabolism to 5-FU but are also involved in its conversion and elimination that may lead to discontinuation, dose reduction or escalation of treatment in order to obtain the best therapeutic ratio. These strategies based on biomarkers may be relevant in the context of the implementation of precision oncology. In particular for head and neck cancers, the identification of biomarkers to select possible cases of severe toxicity requiring discontinuation of treatment, including “multi-omics” approaches, evaluate not only serological biomarkers, but also miRNAs, imaging and radiomics which will ensure capecitabine a role in both induction and concomitant or even adjuvant and palliative settings. An approach including routine testing of dihydropyrimidine dehydrogenase (DPD) or even the thymidine phosphorylase (TP)/DPD ratio and the inclusion of miRNAs, imaging and radiomics parameters in multi-omics models will help implement “precision chemotherapy” in HNC, a concept supported by the importance of avoiding interruptions or treatment delays in this type of cancer. The chemosensitivity and prognostic features of HPV-OPC cancers open new horizons for the use of capecitabine in heavily pretreated metastatic cases. Vorinostat and lapatinib are agents that can be associated with capecitabine in future clinical trials to increase the therapeutic ratio.

Feasibility of using CT radiomic signatures for predicting CD8-T cell infiltration and PD-L1 expression in renal cell carcinoma

Objectives

To identify computed tomography (CT)-based radiomic signatures of cluster of differentiation 8 (CD8)-T cell infiltration and programmed cell death ligand 1 (PD-L1) expression levels in patients with clear-cell renal cell carcinoma (ccRCC).

Methods

Seventy-eight patients with pathologically confirmed localized ccRCC, preoperative multiphase CT and tumor resection specimens were enrolled in this retrospective study. Regions of interest (ROI) of the ccRCC volume were manually segmented from the CT images and processed using a radiomics panel comprising of 1708 metrics. The extracted metrics were used as inputs to three machine learning classifiers: Random Forest, AdaBoost, and ElasticNet to create radiomic signatures for CD8-T cell infiltration and PD-L1 expression, respectively.

Results

Using a cut-off of 80 lymphocytes per high power field, 59 % were classified to CD8 highly infiltrated tumors and 41 % were CD8 non highly infiltrated tumors, respectively. An ElasticNet classifier discriminated between these two groups of CD8-T cells with an AUC of 0.68 (95 % CI, 0.55-0.80). In addition, based on tumor proportion score with a cut-off of > 1 % tumor cells expressing PD-L1, 76 % were PD-L1 positive and 24 % were PD-L1 negative. An Adaboost classifier discriminated between PD-L1 positive and PD-L1 negative tumors with an AUC of 0.8 95 % CI: (0.66, 0.95). 3D radiomics metrics of graylevel co-occurrence matrix (GLCM) and graylevel run-length matrix (GLRLM) metrics drove the performance for CD8-Tcell and PD-L1 classification, respectively.

Conclusions

CT-radiomic signatures can differentiate tumors with high CD8-T cell infiltration with moderate accuracy and positive PD-L1 expression with good accuracy in ccRCC.