Next Article in Journal
Treatment of Bone Metastasis
Previous Article in Journal
Cognitive Trajectories in Older Patients with Cancer Undergoing Radiotherapy—A Prospective Observational Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Novel Harmonization Method for Multi-Centric Radiomic Studies in Non-Small Cell Lung Cancer

by
Marco Bertolini
1,
Valeria Trojani
1,*,
Andrea Botti
1,
Noemi Cucurachi
1,
Marco Galaverni
2,
Salvatore Cozzi
3,
Paolo Borghetti
4,
Salvatore La Mattina
4,
Edoardo Pastorello
4,
Michele Avanzo
5,
Alberto Revelant
6,
Matteo Sepulcri
7,
Chiara Paronetto
7,
Stefano Ursino
8,
Giulia Malfatti
8,
Niccolò Giaj-Levra
9,
Lorenzo Falcinelli
10,
Cinzia Iotti
3,
Mauro Iori
1 and
Patrizia Ciammella
3
1
S.C. Fisica Medica, Azienda USL-IRCCS di Reggio Emilia, 42124 Reggio Emilia, Italy
2
S.C. Radioterapia, Azienda Ospedaliero-Universitaria Maggiore, 43126 Parma, Italy
3
S.C. Radioterapia, Azienda USL-IRCCS di Reggio Emilia, 42124 Reggio Emilia, Italy
4
Department of Radiation Oncology, University and Spedali Civili Hospital, 25123 Brescia, Italy
5
Medical Physics Department, Centro di Riferimento Oncologico di Aviano (CRO) IRCCS, 33081 Aviano, Italy
6
Radiation Oncology Department, Centro di Riferimento Oncologico di Aviano (CRO) IRCCS, 33081 Aviano, Italy
7
Radiotherapy, Veneto Institute of Oncology IOV—IRCCS, 35128 Padova, Italy
8
Department of Radiation Oncology, Santa Chiara University Hospital, 56100 Pisa, Italy
9
Department of Radiation Oncology, IRCCS Sacro Cuore Don Calabria Hospital Negrar, 37024 Verona, Italy
10
Radiation Oncology Section, S. Maria della Misericordia Hospital, 06129 Perugia, Italy
*
Author to whom correspondence should be addressed.
Curr. Oncol. 2022, 29(8), 5179-5194; https://doi.org/10.3390/curroncol29080410
Submission received: 2 May 2022 / Revised: 8 July 2022 / Accepted: 14 July 2022 / Published: 22 July 2022
(This article belongs to the Topic Artificial Intelligence in Cancer Diagnosis and Therapy)

Abstract

:
The purpose of this multi-centric work was to investigate the relationship between radiomic features extracted from pre-treatment computed tomography (CT), positron emission tomography (PET) imaging, and clinical outcomes for stereotactic body radiation therapy (SBRT) in early-stage non-small cell lung cancer (NSCLC). One-hundred and seventeen patients who received SBRT for early-stage NSCLC were retrospectively identified from seven Italian centers. The tumor was identified on pre-treatment free-breathing CT and PET images, from which we extracted 3004 quantitative radiomic features. The primary outcome was 24-month progression-free-survival (PFS) based on cancer recurrence (local/non-local) following SBRT. A harmonization technique was proposed for CT features considering lesion and contralateral healthy lung tissues using the LASSO algorithm as a feature selector. Models with harmonized CT features (B models) demonstrated better performances compared to the ones using only original CT features (C models). A linear support vector machine (SVM) with harmonized CT and PET features (A1 model) showed an area under the curve (AUC) of 0.77 (0.63–0.85) for predicting the primary outcome in an external validation cohort. The addition of clinical features did not enhance the model performance. This study provided the basis for validating our novel CT data harmonization strategy, involving delta radiomics. The harmonized radiomic models demonstrated the capability to properly predict patient prognosis.

1. Introduction

Non-small-cell lung cancer (NSCLC) is, overall, the second-most-common cancer and a leading cause of cancer-related death worldwide, despite recent therapeutic advances [1]. Stage I disease represents approximately 25% of the patients receiving diagnoses of NSCLC and accounts for the most curable cohort of the population [2]. Surgery is the gold standard for these patients: lobectomy with hilar and mediastinal lymph node dissection is the preferred approach, given the Lung Cancer Study Group (LCSG) trial results [3]. Instead, sublobar resection has shown inferior local control and a trend toward decreased survival. However, evaluation of sublobar resection in selected patients is currently underway. The historical standard therapy for unresectable early-stage NSCLC was conventionally fractionated radiation therapy (RT) (e.g., 2 Gy per fraction, for a total dose of 54–60 Gy). However, the reported long-term local control (LC; 30–70%) and overall survival (OS;15–30%) rates with this approach are suboptimal [4,5,6].
Advances in imaging and radiation treatment planning and delivery (e.g., with image guidance and motion management) made the delivery of “ablative doses” of radiation to small targets possible with better results in terms of local control [7,8,9,10,11].
Stereotactic Body radiation therapy (SBRT) has proved to be the first therapeutic option in inoperable stage I NSCLC patients or for those who refuse surgical treatment, with similar rates of local tumor control and overall clinical outcomes [12,13]. Recently, a meta-analysis by Li et al. reported a significant superiority in the local control rate and in 3-year and 5-year OS (54.73% and 29.30 % vs. 39.5 and 27.47) in the SBRT group compared with conventionally fractionated RT [14].
SBRT was reported to have a local control rate in excess of 85% at 3 years [14,15]. Despite consistent clinical outcomes, it is well known that dose fractionation heterogeneity and technical expertise may influence the outcome with SBRT [16,17,18]. A recent study reported that the factors affecting outcomes after SBRT for early-stage NSCLC are Biological Effective Dose (BED) and tumor size [19].
Radiomics is a recent technique introduced in medicine to describe characteristics of medical images quantitatively. Radiomics belongs to artificial intelligence (AI) applications, but it is based on the calculation of features using well-defined mathematical formulas applied directly to the image pixel values (or to a filtered version of the original images). The mathematical definitions of radiomic features are based on the distribution and the relationship between pixels and voxels in the images’ region of interest. The concept behind this method lies in the fact that the human eye cannot appreciate all the characteristics of a medical image. Haralick et al. [20] described how the textural features, highlighting the behavior of gray levels’ dependencies, can identify different areas in an image. Later, textural information was proposed as an application in medical imaging [21,22]. The improvement in hardware calculation power made these techniques able to compute a high number of medical imaging biomarkers in an acceptable span of time; those indices should help the physician during the treatment decision task, allowing a personalized care pathway for different patients. However, these biomarkers are not yet ready to be used in oncology without a robust validation or a demonstration of their reliability [23]. Among them, radiomic indices and feature signatures are increasingly present in the panorama of modern scientific literature [24,25]. The main issue and challenge up to date are to understand how to overcome the limits of this approach [26,27].
To date, in the literature, several studies have investigated the ability of radiomics features in the tumor-healthy tissue differentiation task, both for computed tomography (CT) and positron emission tomography (PET) datasets, as described by Chu et al. [28] that used feature values in a random forest classifier for diagnostic purposes. In another study [29], healthy tissues’ features were used as additional information for an automatic segmentation algorithm. More recently, feature-extracted CT images were combined with BED values to predict tumor response to SBRT [30].
Despite the great work done to date, to our knowledge, there is still no characterization of the radiomic features’ ability to give specific information about healthy tissue compared to the sick one when machine learning models for prognosis are involved. One of the main challenges in the field of radiomics, which makes its clinical application difficult, is the harmonization of the features to be analyzed.
Our present multicentric work aims to propose a novel concept of harmonizing the CT radiomic signal using a combination derived from both the tumor and the healthy contralateral tissue, to overcome the variability typical in each patient in different conditions (i.e., manufacturer/technical characteristics, acquisition, reconstruction protocol, and different anatomy).

2. Materials and Methods

2.1. Study Design

This study was a retrospective multicentric work. It was approved by the Area Vasta Emilia Nord (AVEN) Ethics Committee (ID: 817/2018/OSS*/IRCCSRE). The study was also approved by the ethics committees of all the participating institutions; it was performed in accordance with the principles of Good Clinical Practice (GCP) in respect of the ICH GCP guidelines, the ethical principles contained in the Helsinki declaration and its subsequent updates. Each patient gave informed consent for joining the study.

2.1.1. Patient Cohort

Patients who underwent SBRT for histologically proven diagnosis of primary early-stage NSCLC were retrospectively collected from January 2010 to December 2019. A multicenter research project named “TEXture Analysis of PET/CT in lung cancer patients treated with Stereotactic body radiation therapy (TEXAS)” was designed to involve seven Italian Centers.
Inclusion criteria were: (1) histologically proven diagnosis of NSCLC; (2) early-stage T1–T3N0M0 (TNM 7th edition); (3) patients who underwent SBRT, with treatment biological effective dose BED10 ≥ 100 Gy; and (4) age > 18 years.
Exclusion criteria were: (1) lung tumor greater than 7 cm; (2) histologically proven diagnosis of small cell lung cancer or metastasis; (3) previous thoracic irradiation; (4) presence of bone, lymph node, or visceral metastatic lesions; (5) patients with secondary pulmonary nodules from non-NSCLC or NSCLC; (6) past non-NSCLC tumors with evidence of active disease at the time of SBRT and synchronous non-NSCLC tumors (arising within six months of SBRT diagnosis of NSCLC) with the exception in both cases of non-melanomatous skin tumors.
The patient cohort was divided into training (76 patients from three centers) and external validation (41 patients from the other four centers) datasets. This strategy for the distribution of centers among datasets was made to balance the two groups according to the patients’ outcomes as described in the following sections. The external validation step was a fundamental part of the study in order to confirm the performances obtained in the training phase.

2.1.2. SBRT Details

Conventional computed tomography (CT) simulation scans were obtained. The radiation oncologist contoured gross tumor volume (GTV) on the CT, as part of the therapeutic pathway. A 5–10 mm isotropic margin was added to GTV to generate the planning target volume (PTV). Intensity-modulated radiation therapy (IMRT) was delivered to all patients. The dose normalization ensured that at least 95% of PTV receives 100% of the prescribed dose with a homogeneous distribution. For all patients, ipsilateral and contralateral lung, heart, chest wall, esophagus, spinal cord, and bronchial trees were contoured as organs at risk (OARs).

2.2. Image Acquisition

All patients included in the study had PET/CT images, previously acquired as part of their care pathways, and a pre-treatment CT used for planning of SBRT. The planning CT acquisition protocols and scanning devices differed among institutions, as reported in Table 1. PET image sets, corrected for attenuation, were acquired no more than three months before the start of the treatment. Patients fasted at least 6 h before the injection 18F-FDG tracer and the serum glucose level measured at the injection time was below 160 mg/mL in all patients. PET examinations were performed 60 min after the intravenous administration of the radiotracer using a specific protocol for each institution shown in Table 1.

2.3. Image Segmentation

Computed tomography and PET image sets were exported in DICOM format into a dedicated research computer for radiomics analysis. For the present study, gross tumor volume contouring was separately performed on the CT (manually, referred to as GTVCT) and PET (automatically, hereinafter named GTVPET) images of the pre-treatment PET/CT studies.
Two radiation oncologists with experience in lung cancer contoured each lesion on every sequential slice of the planning CT using standardized window settings for parenchyma (W = 1600 and L = −600), according to EORTC guidelines [31] for all patients. Regarding GTVPET delineation, the radiation oncologists placed a region of interest (ROI) on the area of tumor FDG uptake on PET images and an automatic contour—consisting of the region encompassed by a given fixed percent intensity level relative to the maximum registered tumor activity (40% of SUV max)—was generated. We decided to use this approach as a previous study showed that GTVPET delineation using this fixed threshold was better correlated with the gross tumor (based on pathologic examination) instead of using as basis the manually delineated GTVCT [32].
In order to perform radiomic feature harmonization, we used an ROI from the healthy tissue. This was obtained by copying the GTVCT into a healthy lung region, i.e., the contralateral lung at the same level of the GTV (named Contra_Lung). The Contra Lung initial volume was also shifted by 0.6 and 0.3 cm in six directions for a total of 12 shifts, avoiding the inclusion of surrounding tissues of the healthy lung. These shifts had the aim of simulating the uncertainty in the positioning of the healthy ROI (for future reproducibility of the harmonization method). We chose the shifts in accordance with PET image resolution to account for a likely uncertainty in ROI positioning since PET imaging can be used to localize the GTV before the treatment.
An example picture showing the location of the ROIs mentioned above is shown in Figure 1.

2.4. Outcome

In this study, PFS was considered as the primary endpoint and was converted into a binary outcome, which was set to 1 for patients who were alive and without disease progression at 24 months, 0 otherwise. PFS was defined as the time from the start of the SBRT to documented relapse or death. The use of the 2-year threshold was chosen because it could properly describe the treatment effectiveness. In fact, a preliminary analysis of the Kaplan–Meier curves of PFS after SBRT for our cohort showed that the majority of the progressions occurred in a period ranging from 2 to 3 years.

2.5. Radiomics Analysis

Our analysis followed the steps defined for our radiomic study (Figure 2), which included image preprocessing. The first phase consisted of spatial resampling to an isotropic voxel size to obtain reproducible and rotationally invariant features. Then, image range re-segmentation updated the ROI voxels according to a chosen intensity range to remove all voxels for which intensity values fall outside the selected intensity range. Finally, the images were discretized by intensity, grouping the original intensity values (256) into specific ranges (bins). The aim was to reduce image noise and computational burden. The intensity discretization process fixed the width of the re-segmentation interval and the bin width, defining a new bin for each intensity interval. Selecting the bin width allowed direct control of the absolute range represented on each bin. The image preprocessing of intensity and spatial discretization is described in Supplementary Material Table S1. Intensity discretization parameters were chosen accordingly to the guidelines proposed by Orlhac et al. [33,34].

2.5.1. Feature Extraction

After the image preprocessing steps, feature calculation and their extraction was performed. Features (intensity-based, shape-based, and second-order) were extracted from original images and filtered images (using wavelets, Laplacian of Gaussian (LoG), and gamma modifier filters) [35]. Radiomic features were calculated using a homemade software employing the widely used pyRadiomics library in order to apply pre-determined filters to the original images and compute features from the edited results. The list of the extracted radiomic features can be found at https://pyradiomics.readthedocs.io/en/latest/features.html, (accessed on 15 July 2022).

2.5.2. Harmonization Process

The harmonization process consisted of, for our two available image modalities, calculating features for 14 different ROIs. One of them coincided with GTVCT, the other 13 with the duplicated GTVCT positioned in the healthy region and its shifts, as described in Section 2.3. This allows us to consider operators’ variability in the positioning of the healthy ROI. The general idea of employing this harmonization formula was inspired by another work [36], and it is shown in Equation (1):
f H A R M ( i ) = f G T V ( i ) f H E A L T H Y ( i ) σ ( i )  
where: f H A R M ( i ) is the ith harmonized feature, f H E A L T H Y ( i ) is the median ith feature value calculated on the 13 healthy tissue samples for each patient and modality, and σ ( i ) is the difference between the 75th and the 25th percentile of the ith feature distribution. Shape features were not harmonized. We applied this harmonization only to CT data because of its intrinsic dependence on the protocol acquisition parameters. Furthermore, CT is used for a morphologic and anatomical characterization and pixel values are related to a physical characteristic of the tissue (the linear attenuation coefficient). On the other hand, PET, being a functional imaging modality, is less sensitive to low signal changes in spatial coordinates. Especially in this case, for the healthy lung, in PET pixel values, there is no useful physical information regarding a region where we do not register a signal from the radiotracer absorption.

2.5.3. Feature Selection

LASSO feature selection was applied, in which the following function is minimized (Equation (2))
i = 1 n ( y i β 0 j = 1 p β j x i j ) 2 + λ j = 1 p | β j |  
where: y i is the observed value, β j is the LASSO coefficients, and λ j = 1 p | β j | is the shrinkage penalty [37].
The parameter λ was chosen using 10-fold cross-validation (CV) computing its error. LASSO penalty brings to zero the weight coefficients ( β j ) of irrelevant features not predictive of the chosen outcome. In addition, LASSO handles sets of collinear features by increasing the weight of one of them while setting the other weights to zero. Because the considered outcome was binary, we used a binomial function for LASSO regression. In Table S2 we show the shrinkage penalties for our trained models.

2.5.4. Model Building

The original and harmonized features were used to develop a supervised machine learning binary classifier. A linear support vector machine (SVM, Model 1) [38] and an Ensemble Subspace Discriminant (ESD, Model 2) [39] were trained by optimizing their performance in 10-fold cross-validation in the training dataset.
Linear SVM classifiers provide low generalization error, even with small learning sample datasets. ESD classifiers are used to decide an explicit discriminant subspace of low dimension.
The two described model types were applied to five different combinations of input features: (A) harmonized CT + PET, (B) harmonized CT, (C) original CT, (D) only original PET, and (E) harmonized CT + PET + selected clinical variables in order to assess the effect of harmonization on the performance of the predictive models. The interested reader can find more information in Text S1 in Supplementary Materials.
The clinical variables in method (E) were chosen among the available ones by using Kaplan–Meier survival curves as described in Section 2.5.5.

2.5.5. Statistical Analysis

For each model, a 95% confidence interval (CI) of the AUCs was calculated for the training and external validation sets. Furthermore, accuracy (95% CI are reported), precision, and recall were calculated.
Subsequently, Kaplan–Meier survival curves were computed using the PFS to select the clinical features. A clinical feature exhibiting a p-value from a log-rank test less than 0.05 was considered significant and included in model E. Matlab R2021b (Mathworks, Natick, MA) and R (Vienna, Austria), available at https://www.R-project.org (accessed on 15 July 2022), were used to perform the statistical analysis.
The p-values related to statistical differences among the AUC values of each model were calculated using two-sided DeLong test.

3. Results

3.1. Clinical Results

One-hundred and seventeen early-stage NSCLC patients met the inclusion criteria. The baseline characteristics of the patients are summarized in Table 2. The median age was 78 years and there were more male (72.6%) than female patients. With a median follow-up of 29.8 months, the median PFS was 24.2 months, and 2-year PFS percentage was 51.2%. Median OS and 2-year OS percentage were 28.5 months and 64%, respectively. The clinical characteristics, including age, gender, Charlson comorbidity index (CCI), diffusing capacity of carbon monoxide (DLCO), tumor size, Eastern Cooperative Oncology Group (ECOG) performance status, and biological equivalent dose to PTV, showed no significant differences between the training and external validation cohorts (Table 2).
No clinical or treatment-related features were shown to be significantly related to PFS in the univariate analysis of the whole population, except for gender (p = 0.04 in favor of female) and lung site (right vs. left in favor of the right one, p = 0.006).

3.2. PFS Models

The PFS predictive performance of the linear SVM and ESD models using radiomic features and clinical features are reported in Table 3. In Figure 3, all the models are graphically compared considering their confidential intervals. Models using harmonized features and PET (A,E) achieved AUCs greater than 0.70, both in training and validation. The performances of models using CT-only harmonized features (B models) are not confirmed on the validation dataset (AUC training > 0.75, AUC validation <0.60), while adding PET features leads to better stability between the training and validation sets. Only C-type models (original CT-only features) showed a low-mean AUC (< 0.62), both in training and validation.
It is worth noting that both A1 and A2 models significantly outperformed C1 and C2 models, both in the training and external validation datasets (p = 0.0001, p = 0.01 and p = 0.02, p = 0.046, respectively, for linear SVM and subspace discriminant models), likewise for E and C models (E1: p < 0.0001, and p = 0.02, and E2: p = 0.01, and p = 0.02, respectively, for training and external validation datasets). C models outperformed B models, but only in the training dataset (p < 0.0001 and p = 0.01, respectively, for linear SVM and subspace discriminant models) and D and C models had the same performance metrics. In summary, models using clinical information (E models) do not add a significant improvement to A models. This effect is also appreciable in Figure 4, where all the p-values among the models are represented.

4. Discussion

In this work, a multi-centric cohort of early-stage NSCLC patients treated with SBRT was used to build and validate predictive models of PFS greater than 24 months using radiomic features from CT and PET exams and clinical information. Several existing works in the literature [40] describe acquisition protocol variabilities in multi-centric studies, which could affect the performance of radiomic models. Since radiomics computes features from the pixel values in the images, differences in acquisition protocol can lead to biased results. The rationale behind our harmonization method lies in the fact that retrospective multicenter radiomic studies are challenging but necessary, as gathering data from several centers for a centralized analysis is complex for legal, ethical, administrative, and technical reasons. Most of the time, the different centers involved do not follow standardized acquisition and reconstruction protocols; therefore, the collected data suffer from intra-, and inter-variability, making radiomic features sensitive to multicenter variability. Our novel harmonization technique aims to reduce the bias caused by the absence of standardized protocols. Generally, feature analysis is performed by calculating them inside an ROI that coincides with the lesion target. Our study aims to tackle this issue by attempting to reduce this effect using the healthy region of the patient as the baseline from which to harmonize the radiomic data computed from the lesion. From our results, the harmonization improves models’ performance when it is used on CT image sets. On the other hand, we expected that for PET-only images, the harmonization method is not easily applicable due to the functional aim of this imaging modality. In such a modality, several healthy tissues (i.e., lungs) are not 18-FDG-avid, while a harmonization based on healthy radiotracer accumulation, to our knowledge, has yet to be studied. Our work investigated and evaluated the feasibility of this technique, which could be employed and better analyzed in future studies.
When using original feature values, PET features were preferred over CT during feature selection, resulting in an only PET-based model. Furthermore, in A models, two features from CT were included in the final prediction score (Log_Sigma30mm_GLDM_Small-Dependence-High-Gray-Level-Emphasis (SDHGLE) and Wavelet_LHH_NGTDM_Busyness), which were also selected in B models. In the same way, a subset of selected features in the A method (Log_sigma10mm_GLSZM_Size-Zone-Low-Gray-Level-Zone-Emphasis (LGLZE), Exponential_FIRST ORDER_Median, Square_GLSZM_ZoneEntropy (ZE)) is also present among the selected features in the PET-only-based method. This could mean that our approach, also given the higher performance metric of the A method, was able to merge the information hidden in the CT and PET image sets in a multi-centric cohort of patients, highlighting the importance of properly handling hybrid imaging in radiomic models.
Other harmonization techniques were previously described in the literature [41]. Among these, the ComBat harmonization method, which removes batch effects, mostly based on an empirical Bayes framework, is one of the most used. Conversely, the ComBat method has some limitations: for instance, the dimension of the homogeneous group cannot be too few (in our study, it would have not been applicable as most of the centers provided less than 15 patients). Indeed, recently, some methods have been proposed to overcome these limitations, e.g., using bootstrap and Monte Carlo technique to improve robustness in the estimation [42].
Even if Monte Carlo and bootstrap strategies aimed at overcoming the cohort size limitations, the objective of this method still focuses on removing differences in radiomic feature distributions among different labels (corresponding to different centers). ComBat, thus, relies on the individual distributions and changes made to feature values are dependent on a group of patients. While we know that there are data supporting the effectiveness of this method (especially in making the feature distributions uniform), we wanted to tackle the multicenter studies issue from a different angle, which is to account for the individual patients’ differences (caused both by the scanner/institution protocols and their anatomy) taken directly from the lesion imaging. This renders the method easier, both computationally and for cohort eligibility reasons (which, in the ComBat method, is needed for representing the single center in terms of homogeneity being an assumption of the method). In fact, if validated further, our method can be applied even in heterogenous cohorts since it uses only the single image set of the patient.
Our approach aimed to use all the information present in the CT data, both from cancer lesions and the contralateral healthy tissue, simulating the radiologists’ skill in subjectively evaluating a lesion and adding this information in quantifiable and statistical terms (through the radiomic features).
In our work, the well-known and studied concept of delta radiomics was implemented not in a temporal sense but spatially (cancer vs. healthy tissue), which is an approach that, to our knowledge, was not applied in other prognostic oncological works. Traditional radiomics uses absolute values extracted from regions of interest to predict a clinical outcome. On the other hand, delta radiomics predicts a clinical outcome through the combination of radiomic values computed from image sets acquired at different time points (i.e., radiographs to monitor follow-ups or differences between basal PET and interim PET), which is a rationale also used in clinical practice to assess lesion progression (i.e., PERCIST). In our manuscript, we decided to apply delta radiomics not between different time points but between different anatomic locations (healthy vs. tumor tissues). The assumption behind this use of delta radiomics is that each patient can have an intrinsic “baseline” value for a certain radiomic feature (caused by individual anatomy and institution protocols) that needs to be accounted for when building predictive models. Comparison between normal and tumor tissue behavior (even in terms of pixel values) is also common in clinical practice (i.e., SUV values typical of physiological metabolism or HU/density values of healthy tissue). Some authors [43,44,45] explain that delta radiomics—which is the use of textural indices associated with different time points or anatomical regions—is more successful than traditional radiomics. Our work aims to provide the basis for a framework where the study of simple absolute feature values can make room for the analysis of their relationship to a reference, used as a threshold or as a comparison.
There are several limitations to the current work. Our study suffered a restricted number of patients selected retrospectively. Nonetheless, the patients’ number seemed reasonable at the current phase of our study. It assures the homogeneity in terms of pathologies, as including only NSCLC lesions prevented possible biases created by evaluating different diseases, even in the lung anatomical district. Our study highlighted the necessity to monitor and carefully use features related to pixel values and their relationships. In this case, we can assume that the importance of the radiomic features is not only held in their numerical value but also in the intrinsic relationship among those values. The textural indices’ ability to perform more complex clinical tasks (i.e., predicting toxicity and its grade) could be further examined in the next phase of our work or in a prospective study design, which could also assess the robustness of our method. We believe that a prospective study will be able to validate these models within a cohort gathered with a better strategy.
Interestingly enough, we found that the improved performance in models employing harmonized features in the training phase was also confirmed in the validation dataset; the use of an external dataset is becoming more and more crucial to radiomics studies to assure and facilitate their introduction in clinical practice.
In our study, no clinical or treatment-related features were shown to be significantly related to PFS, except for gender and lung site. It is well known in the literature that gender is a prognostic factor for PFS [46,47,48]. Due to the size of our cohort, we did not find significant correlations between PFS and other studied clinical prognostic variables, such as age or histology (also due to the inclusion criteria). To our knowledge, we did not find any other study reporting a significant correlation between PFS and tumor laterality; thus, we will investigate this finding together with our model generalization power in a future prospective study. Indeed, in the literature, many studies showed a significant correlation between some clinical or treatment-related characteristics and outcomes (PFS and OS) and some created predictive models. Among the various statistical prediction models, nomograms can be accurate and feasible prognostic instruments with high utility in estimating individual patient risk and may, thus, help guide treatment decisions in clinical practice. At present, there are some nomograms, based on clinical features, developed for early-stage NSCLC treated with SBRT [49,50,51], but there is still need for validation of the clinical variables found in those studies and their experimental results in more robust cohorts, such as prospective ones. Therefore, a need exists for a robust recurrence-related prediction model to help select high-risk candidates who may benefit from additional systemic therapies.
In this scenario, a predictive model based on radiomics and clinical and treatment-related characteristics can improve the prediction of clinical outcomes, as already demonstrated by other works. We also found out that clinical variables did not improve the radiomics models, but only the proposed harmonization process statistically significantly improves the model’s performance.
As previously stated, our future aim is to apply our method to a prospective multicentric cohort to further validate the framework’s stability. In addition, other anatomical regions should be explored to generalize the harmonization, even when the healthy area is not so easily defined as in the lung case. Regarding the employment of this method also in PET datasets, it could be interesting to explore the feasibility of applying our harmonization to 18F-FDG-avid anatomical regions, such as the liver or the brain, which are, however, not related to a pathologic response. Such a method could be especially useful where CT-PET is the only exam included in the care pathway of the patient.

5. Conclusions

A novel strategy of CT data harmonization involving delta radiomics, considering both cancer and healthy tissue in the contralateral lung, was tested and externally validated in a multi-centric study for NSCLC patients, to initially assess its feasibility.
The radiomics models with harmonized features can predict better the selected patient outcome in our cohort, providing valuable additional information to the clinician.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/curroncol29080410/s1, Table S1: The image preprocessing of intensity and spatial discretization, Table S2: Lasso feature selection, Text S1: Models’ description.

Author Contributions

M.B.: study, design, literature search, data analysis/interpretation, writing original draft, supervision, V.T.: study, design, literature search, data analysis/interpretation, writing original draft, A.B.: Data analysis/interpretation, data validation, reviewing draft, N.C.: literature search, data analysis/interpretation, data validation, figures, writing, M.G.: study design, data collection, S.C.: literature search, reviewing draft, P.B.: resources, data collection, S.L.M.: resources, data collection, E.P.: resources, data collection, M.A.: resources, data collection, reviewing draft, A.R.: resources, data collection, M.S.: resources, data collection, C.P.: resources, data collection, S.U.: resources, data collection, G.M.: resources, data collection, N.G.-L.: resources, data collection, L.F.: resources, data collection, C.I.: resources, data collection, M.I.: resources, data collection, P.C.: resources, data collection, supervision. All authors have read and agreed to the published version of the manuscript.

Funding

This study was partially supported by the Italian Ministry of Health—Ricerca Corrente.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Ethics Committee of by the Area Vasta Emilia Nord (AVEN) (protocol code 817/2018/OSS*/IRCCSRE approved on 16 July 2019).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The training weights of the models proposed in this work (namely, models E and A described in Section 2.5.4) are available in a GitHub repository at the link: https://github.com/ausl-re/TEXAS, accessed on 15 July 2022. Further instructions on how to perform the harmonization and to use the model will be added.

Acknowledgments

A heartfelt thanks to all the people involved in the various centers in collecting and managing the data, especially Simona Marani, Maria Paola Ruggieri, Giulia Mascari, and Cinthia Aristei. A special thanks to the AIRO Lung Group, especially Vieri Scotti, Stefano Vagge, and Alessio Bruni.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Siegel, R.; Naishadham, D.; Jemal, A. Cancer statistics, 2013. CA Cancer J. Clin. 2013, 63, 11–30. [Google Scholar] [CrossRef] [Green Version]
  2. Edge, S.B.; Byrd, D.R.; Compton, C.C.; Fritz, A.G.; Greene, F.L.; Trotti, A. AJCC Cancer Staging Manual, 7th ed.; Springer: New York, NY, USA, 2010; pp. 253–264. [Google Scholar]
  3. Ginsberg, R.J.; Rubinstein, L.V. Randomized trial of lobectomy versus limited resection for T1 N0 non-small cell lung cancer. Lung Cancer Study Group. Ann. Thorac. Surg. 1995, 60, 615–622. [Google Scholar] [CrossRef]
  4. Bogart, J.A. Fractionated radiotherapy for high-risk patients with early-stage non-small cell lung cancer. Semin. Thorac. Cardiovasc. Surg. 2010, 22, 44–52. [Google Scholar] [CrossRef]
  5. Hayman, J.A.; Martel, M.K.; Ten Haken, R.K.; Normolle, D.P.; Todd, R.F., III; Littles, J.F.; Sullivan, M.A.; Possert, P.W.; Turrisi, A.T.; Lichter, A.S. Dose escalation in non-small-cell lung cancer using three-dimensional conformal radiation therapy: Update of a phase I trial. J. Clin. Oncol. 2001, 19, 127–136. [Google Scholar] [CrossRef]
  6. Qiao, X.; Tullgren, O.; Lax, I.; Sirzén, F.; Lewensohn, R. The role of radiotherapy in the treatment of stage I non-small cell lung cancer. Lung Cancer 2003, 41, 1–11. [Google Scholar] [CrossRef]
  7. Papiez, L.; Timmerman, R.; DesRosiers, C.; Randall, M. Extracranial stereotactic radioablation: Physical principles. Acta Oncol. 2003, 42, 882–894. [Google Scholar] [CrossRef] [Green Version]
  8. Timmerman, R.; Papiez, L.; McGarry, R.; Likes, L.; DesRosiers, C.; Frost, S.; Williams, M. Extracranial stereotactic radioablation: Results of a phase I study in medically inoperable stage I non-small cell lung cancer. Chest 2003, 124, 1946–1955. [Google Scholar] [CrossRef]
  9. Timmerman, R.; Paulus, R.; Galvin, J.; Michalski, J.; Straube, W.; Bradley, J.; Fakiris, A.; Bezjak, A.; Videtic, G.; Johnstone, D.; et al. Stereotactic body radiation therapy for inoperable early-stage lung cancer. JAMA 2010, 303, 1070–1076. [Google Scholar] [CrossRef] [Green Version]
  10. Timmerman, R.D.; Herman, J.; Cho, L.C. Emergence of stereotactic body radiation therapy and its impact on current and future clinical practice. J. Clin. Oncol. 2014, 32, 2847–2854. [Google Scholar] [CrossRef]
  11. Videtic, G.M.; Stephans, K.L.; Woody, N.M.; Reddy, C.A.; Zhuang, T.; Magnelli, A.; Djemil, T. 30 Gy or 34 Gy? Comparing 2 single-fraction SBRT dose schedules for stage I medically inoperable non-small cell lung cancer. Int. J. Radiat. Oncol. Biol. Phys. 2014, 90, 203–208. [Google Scholar] [CrossRef]
  12. Simone, C.B., 2nd; Wildt, B.; Haas, A.R.; Pope, G.; Rengan, R.; Hahn, S.M. Stereotactic body radiation therapy for lung cancer. Chest 2013, 143, 1784–1790. [Google Scholar] [CrossRef] [Green Version]
  13. Ackerson, B.G.; Tong, B.C.; Hong, J.C.; Gu, L.; Chino, J.; Trotter, J.W.; D’Amico, T.A.; Torok, J.A.; Lafata, K.; Chang, C.; et al. Stereotactic body radiation therapy versus sublobar resection for stage I NSCLC. Lung Cancer 2018, 125, 185–191. [Google Scholar] [CrossRef]
  14. Li, C.; Wang, L.; Wu, Q.; Zhao, J.; Yi, F.; Xu, J.; Wei, Y.; Zhang, W. A meta-analysis comparing stereotactic body radiotherapy vs. conventional radiotherapy in inoperable stage I non-small cell lung cancer. Medicine 2020, 99, e21715. [Google Scholar] [CrossRef]
  15. Fakiris, A.J.; McGarry, R.C.; Yiannoutsos, C.T.; Papiez, L.; Williams, M.; Henderson, M.A.; Timmerman, R. Stereotactic body radiation therapy for early-stage non-small-cell lung carcinoma: Four-year results of a prospective phase II study. Int. J. Radiat. Oncol. Biol. Phys. 2009, 75, 677. [Google Scholar] [CrossRef]
  16. Koshy, M.; Malik, R.; Mahmood, U.; Husain, Z.; Sher, D.J. Stereotactic body radiotherapy and treatment at a high volume facility is associated with improved survival in patients with inoperable stage I non-small cell lung cancer. Radiother. Oncol. 2015, 114, 148–154. [Google Scholar] [CrossRef]
  17. Olsen, J.R.; Robinson, C.G.; El Naqa, I.; Creach, K.M.; Drzymala, R.E.; Bloch, C.; Parikh, P.J.; Bradley, J.D. Dose Response for stereotactic body radiotherapy in early-stage nonsmall-cell lung cancer. Int. J. Radiat. Oncol. Biol. Phys. 2013, 81, e299–e303. [Google Scholar] [CrossRef]
  18. Onishi, H.; Shirato, H.; Nagata, Y.; Hiraoka, M.; Fujino, M.; Gomi, K.; Niibe, Y.; Karasawa, K.; Hayakawa, K.; Takai, Y.; et al. Hypofractionated Stereotactic Radiotherapy (HypoFXSRT) for stage I non-small cell lung cancer: Updated results of 257 patients in a Japanese multi-institutional study. J. Thorac. Oncol. 2011, 2, S94–S100. [Google Scholar] [CrossRef] [Green Version]
  19. Lee, P.; Loo, B.W., Jr.; Biswas, T.; Ding, G.X.; El Naqa, I.M.; Jackson, A.; Kong, F.M.; LaCouture, T.; Miften, M.; Solberg, T.; et al. Local Control After Stereotactic Body Radiation Therapy for Stage I Non-Small Cell Lung Cancer. Int. J. Radiat. Oncol. Biol. Phys. 2021, 110, 160–171. [Google Scholar] [CrossRef]
  20. Haralick, R.M.; Shanmugam, K.; Dinstein, I. Textural Features for Image Classification. IEEE Trans. Syst. Man. Cybern. 1973, SMC-3, 610–621. [Google Scholar] [CrossRef] [Green Version]
  21. Castellano, G.; Bonilha, L.; Li, L.M.; Cendes, F. Texture analysis of medical images. Clin. Radiol. 2004, 59, 1061–1069. [Google Scholar] [CrossRef]
  22. Tourassi, G.D. Journey toward computer-aided diagnosis: Role of image texture analysis. Radiology 1999, 213, 317–320. [Google Scholar] [CrossRef] [PubMed]
  23. Yip, S.S.; Aerts, H.J. Applications and limitations of radiomics. Phys. Med. Biol. 2016, 61, R150–R166. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Lambin, P.; Rios-Velazquez, E.; Leijenaar, R.; Carvalho, S.; van Stiphout, R.G.; Granton, P.; Zegers, C.M.; Gillies, R.; Boellard, R.; Dekker, A.; et al. Radiomics: Extracting more information from medical images using advanced feature analysis. Eur. J. Cancer 2012, 48, 441–446. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Lambin, P.; Leijenaar, R.T.-H.; Deist, T.M.; Peerlings, J.; de Jong, E.E.C.; van Timmeren, J.; Sanduleanu, S.; Larue, R.T.H.M.; Even, A.J.G.; Jochems, A.; et al. Radiomics: The bridge between medical imaging and personalized medicine. Nat. Rev. Clin. Oncol. 2017, 14, 749–762. [Google Scholar] [CrossRef]
  26. Drucker, E.; Krapfenbauer, K. Pitfalls and limitations in translation from biomarker discovery to clinical utility in predictive and personalised medicine. EPMA J. 2013, 4, 7. [Google Scholar] [CrossRef] [Green Version]
  27. Limkin, E.J.; Sun, R.; Dercle, L.; Zacharaki, E.I.; Robert, C.; Reuzé, S.; Schernberg, A.; Paragios, N.; Deutsch, E.; Ferté, C. Promises and challenges for the implementation of computational medical imaging (radiomics) in oncology. Ann. Oncol. 2017, 28, 1191–1206. [Google Scholar] [CrossRef]
  28. Chu, L.C.; Park, S.; Kawamoto, S.; Fouladi, D.F.; Shayesteh, S.; Zinreich, E.S.; Graves, J.S.; Horton, K.M.; Hruban, R.H.; Yuille, A.L.; et al. Utility of CT Radiomics Features in Differentiation of Pancreatic Ductal Adenocarcinoma from Normal Pancreatic Tissue. Am. J. Roentgenol. 2019, 213, 349–357. [Google Scholar] [CrossRef]
  29. Markel, D.; Caldwell, C.; Alasti, H.; Soliman, H.; Ung, Y.; Lee, J.; Sun, A. Automatic Segmentation of Lung Carcinoma Using 3D Texture Features in 18-FDG PET/CT. Int. J. Mol. Imaging. 2013, 2013, 980769. [Google Scholar] [CrossRef] [Green Version]
  30. Avanzo, M.; Gagliardi, V.; Stancanello, J.; Blanck, O.; Pirrone, G.; El Naqa, I.; Revelant, A.; Sartor, G. Combining computed tomography and biologically effective dose in radiomics and deep learning improves prediction of tumor response to robotic lung stereotactic body radiation therapy. Med. Phys. 2021, 48, 6257–6269. [Google Scholar] [CrossRef]
  31. De Ruysscher, D.; Faivre-Finn, C.; Moeller, D.; Nestle, U.; Hurkmans, C.W.; le Péchoux, C.; Belderbos, J.; Guckenberger, M.; Senan, S. Lung Group and the Radiation Oncology Group of the European Organization for Research and Treatment of Cancer (EORTC). European Organization for Research and Treatment of Cancer (EORTC) recommendations for planning and delivery of high-dose, high precision radiotherapy for lung cancer. Radiother. Oncol. 2017, 124, 1–10. [Google Scholar]
  32. Van Loon, J.; Siedschlag, C.; Stroom, J.; Blauwgeers, H.; van Suylen, R.J.; Knegjens, J.; Rossi, M.; van Baardwijk, A.; Boersma, L.; Klomp, H.; et al. Microscopic disease extension in three dimensions for non–small-cell lung cancer: Development of a prediction model using pathology-validated positron emission tomography and computed tomography features. Int. J. Radiat. Oncol. Biol. Phys. 2012, 82, 448–456. [Google Scholar] [CrossRef] [PubMed]
  33. Reuzé, S.; Orlhac, F.; Chargari, C.; Nioche, C.; Limkin, E.; Riet, F.; Escande, A.; Haie-Meder, C.; Dercle, L.; Gouy, S.; et al. Prediction of cervical cancer recurrence using textural features extracted from 18F-FDG PET images acquired with different scanners. Oncotarget 2017, 8, 43169–43179. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Orlhac, F.; Soussan, M.; Chouahnia, K.; Martinod, E.; Buvat, I. 18F-FDG PET-Derived Textural Indices Reflect Tissue-Specific Uptake Pattern in Non-Small Cell Lung Cancer. PLoS ONE 2015, 10, e0145063. [Google Scholar] [CrossRef]
  35. Van Griethuysen, J.J.M.; Fedorov, A.; Parmar, C.; Hosny, A.; Aucoin, N.; Narayan, V.; Beets-Tan, R.G.H.; Fillon-Robin, J.C.; Pieper, S.; Aerts, H.J.W.L. Computational Radiomics System to Decode the Radiographic Phenotype. Cancer Res. 2017, 77, e104–e107. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Haga, A.; Takahashi, W.; Aoki, S.; Nawa, K.; Yamashita, H.; Abe, O.; Nakagawa, K. Standardization of imaging features for radiomics analysis. J. Med. Investig. 2019, 66, 35–37. [Google Scholar] [CrossRef]
  37. James, G.; Witten, D.; Hastie, T.; Tibshirani, R. Support vector machines. In An Introduction to Statistical Learning: Applications in R; Springer: New York, NY, USA, 2013; p. 359. [Google Scholar]
  38. Cortes, C.; Vapnik, V. Support-vector networks. Mach. Learn. 1995, 20, 273–297. [Google Scholar] [CrossRef]
  39. Polikar, R. Ensemble based systems in decision making. IEEE Circuits Syst. Mag. 2006, 6, 21–45. [Google Scholar] [CrossRef]
  40. Kothari, G.; Korte, J.; Lehrer, E.J.; Zaorsky, N.G.; Lazarakis, S.; Kron, T.; Hardcastle, N.; Siva, S. A systematic review and meta-analysis of the prognostic value of radiomics based models in non-small cell lung cancer treated with curative radiotherapy. Radiother. Oncol. 2021, 155, 188–203. [Google Scholar] [CrossRef]
  41. Mali, S.A.; Ibrahim, A.; Woodruff, H.C.; Andrearczyk, V.; Müller, H.; Primakov, S.; Salahuddin, Z.; Chatterjee, A.; Lambin, P. Making Radiomics More Reproducible across Scanner and Imaging Protocol Variations: A Review of Harmonization Methods. J. Pers. Med. 2021, 11, 842. [Google Scholar] [CrossRef]
  42. Da-ano, R.; Masson, I.; Lucia, F.; Doré, M.; Robin, P.; Alfieri, J.; Rousseau, C.; Mervoyer, A.; Reinhold, C.; Castelli, J.; et al. Performance comparison of modified ComBat for harmonization of radiomic features for multicenter studies. Sci. Rep. 2020, 10, 10248. [Google Scholar] [CrossRef]
  43. Fave, X.; Zhang, L.; Yang, J.; Mackin, D.; Balter, P.; Gomez, D.; Followill, D.; Jones, A.K.; Stingo, F.; Liao, Z.; et al. Delta-radiomics features for the prediction of patient outcomes in non-small cell lung cancer. Sci. Rep. 2017, 7, 588. [Google Scholar] [CrossRef] [PubMed]
  44. Alahmari, S.S.; Cherezov, D.; Goldgof, D.B.; Hall, L.O.; Gillies, R.J.; Schabath, M.B. Delta Radiomics Improves Pulmonary Nodule Malignancy Prediction in Lung Cancer Screening. IEEE Access 2018, 6, 77796–77806. [Google Scholar] [CrossRef] [PubMed]
  45. Nardone, V.; Reginelli, A.; Guida, C.; Belfiore, M.P.; Biondi, M.; Mormile, M.; Banci Buonamici, F.; di Giorgio, E.; Spadafora, M.; Tini, P.; et al. Delta-radiomics increases multicentre reproducibility: A phantom study. Med. Oncol. 2020, 37, 38. [Google Scholar] [CrossRef] [PubMed]
  46. Pinto, J.A.; Vallejos, C.S.; Raez, L.E.; Mas, L.A.; Ruiz, R.; Torres-Roman, J.S.; Morante, Z.; Araujo, J.M.; Gómez, H.L.; Aguilar, A.; et al. Gender and outcomes in non-small cell lung cancer: An old prognostic variable comes back for targeted therapy and immunotherapy? ESMO Open 2018, 3, e000344. [Google Scholar] [CrossRef] [Green Version]
  47. De Perrot, M.; Licker, M.; Bouchardy, C.; Usel, M.; Robert, J.; Spiliopoulos, A. Sex differences in presentation, management, and prognosis of patients with non-small cell lung carcinoma. J. Thorac. Cardiovasc. Surg. 2000, 119, 21–26. [Google Scholar] [CrossRef] [Green Version]
  48. Hsu, L.H.; Chu, N.M.; Liu, C.C.; Tsai, S.Y.; You, D.L.; Ko, J.S.; Lu, M.C.; Feng, A.C. Sex-associated differences in non-small cell lung cancer in the new era: Is gender an independent prognostic factor? Lung Cancer 2009, 66, 262–267. [Google Scholar] [CrossRef]
  49. Louie, A.V.; Haasbeek, C.J.; Mokhles, S.; Rodrigues, G.B.; Stephans, K.L.; Lagerwaard, F.J.; Palma, D.A.; Videtic, G.M.; Warner, A.; Takkenberg, J.J.; et al. Predicting Overall Survival After Stereotactic Ablative Radiation Therapy in Early-Stage Lung Cancer: Development and External Validation of the Amsterdam Prognostic Model. Int. J. Radiat. Oncol. Biol. Phys. 2015, 93, 82–90. [Google Scholar] [CrossRef]
  50. Ye, L.; Shi, S.; Zeng, Z.; Huang, Y.; Hu, Y.; He, J. Nomograms for predicting disease progression in patients of Stage I non-small cell lung cancer treated with stereotactic body radiotherapy. Jpn. J. Clin. Oncol. 2018, 48, 160–166. [Google Scholar] [CrossRef]
  51. Kang, J.; Ning, M.S.; Feng, H.; Li, H.; Bahig, H.; Brooks, E.D.; Welsh, J.W.; Ye, R.; Miao, H.; Chang, J.Y. Predicting 5-Year Progression and Survival Outcomes for Early Stage Non-small Cell Lung Cancer Treated with Stereotactic Ablative Radiation Therapy: Development and Validation of Robust Prognostic Nomograms. Int. J. Radiat. Oncol. Biol. Phys. 2020, 106, 90–99. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Visualization of the CT ROIs in a patient. The contralateral ROI was shifted in 12 different positions (shown in red).
Figure 1. Visualization of the CT ROIs in a patient. The contralateral ROI was shifted in 12 different positions (shown in red).
Curroncol 29 00410 g001
Figure 2. Radiomic pipeline description of the implemented steps in our evaluation process.
Figure 2. Radiomic pipeline description of the implemented steps in our evaluation process.
Curroncol 29 00410 g002
Figure 3. Performances (AUC) of the studied models. The boxplot shows the minimum, maximum, and average values of the bootstrapped 95% CIs.
Figure 3. Performances (AUC) of the studied models. The boxplot shows the minimum, maximum, and average values of the bootstrapped 95% CIs.
Curroncol 29 00410 g003
Figure 4. p-values calculated using the two-sided DeLong test. Numbers in bold mean the statistical significance.
Figure 4. p-values calculated using the two-sided DeLong test. Numbers in bold mean the statistical significance.
Curroncol 29 00410 g004
Table 1. Protocol acquisition parameters for simulation CT and PET examinations stratified for centers. Whenever two scanners were used, a “|” indicated the different configurations.
Table 1. Protocol acquisition parameters for simulation CT and PET examinations stratified for centers. Whenever two scanners were used, a “|” indicated the different configurations.
CT
CenterkVmAs (Min–Max)Slice Thickness (mm)Manufacturer (s)Convolution KernelRecon Diameter
TRAINBS120191–4013.0PHILIPSB500
RE12083–3553.0GESTD +500
PD12070–3632.5GEBODY FILTER500
EXT VALAV120108–1382.5PHILIPSB500
NE12040–733.0SIEMENSB30f500
PI12027–2362.0SIEMENSB30f–B31s500
PG12080–2002.5–3GESTD +500
PET
CenterSlice thickness (mm)Manufacturer (s)Recon diameterRecon method
TRAINBS3.27GE700–8153D IR/VPFXS
RE3.27GE700–7003D IR/VPFXS
PD2–4PHILIPS|SIEMENS576–8153D-RAMLA/BLOB-OS-TF(PHILIPS)|PSF 3i21s/(SIEMENS)
EXT VALAV4PHILIPS|GE500–700BLOB-OS-TF/VPFXS
NE2–5SIEMENS576–700PSF+TOF 3i21s
PI3.27GE|PHILIPS576–7003D IR (GE)|BLOB-OS-TF(PHILIPS)
PG3.27GE|SIEMENS600–700OSEM|OSEM 2i8s
Table 2. Statical analysis of clinical variables. Abbreviations: PS: Performance status according to ECOG scale, BPCO: chronic obstructive pulmonary disease; ADK: Adenocarcinoma, SCC: squamous cell carcinoma, Fr: fraction; RT: radiotherapy VMAT: volumetric arc-therapy; IMRT: intensity modulated radiotherapy, TOMO: Tomotherapy, PTV: planning target volume. p-values in bold mean the statistical significance.
Table 2. Statical analysis of clinical variables. Abbreviations: PS: Performance status according to ECOG scale, BPCO: chronic obstructive pulmonary disease; ADK: Adenocarcinoma, SCC: squamous cell carcinoma, Fr: fraction; RT: radiotherapy VMAT: volumetric arc-therapy; IMRT: intensity modulated radiotherapy, TOMO: Tomotherapy, PTV: planning target volume. p-values in bold mean the statistical significance.
CharacteristicsTraining Cohort
(N = 76)
External Validation Co#Hort
(N = 41)
p
Gender
Male61240.04
Female1517
Age (years)78 [51–87]79 [57–88]0.72
Smoking Status
Yes50270.22
No2614
Performance Status
037180.75
13515
247
BMI25.2 [16.4–37.1]24.8 [18.3–44.7]0.17
Diabetes mellitus
Yes16120.58
No6029
BPCO
Yes43170.54
No1924
Charlson Comorbidity Index (CCI)
Median6.560.55
Range[3–13][4–10]
T diameter
Median2.352.30.58
Range[0.6–5.5][0.72–27]
Lesion type
Subsolid540.42
Solid7137
Lung Side
Lung right42220.006
Lung left3419
Lobe Site
Upper Lobe44230.89
Lower Lobe3015
Middle Lobe21
Lesion Site
Peripheral55340.92
Central217
BED10
Median115.51000.64
Range[100–180][100–132]
Table 3. Models’ results in terms of AUC, accuracy, precision, and recall.
Table 3. Models’ results in terms of AUC, accuracy, precision, and recall.
Harmo CT + Original PET Features (A)
Linear SVM (A1)
AUC *AccuracyPrecision **Recall **p ***
Training dataset0.77 [0.66–0.87]0.72 ± 0.020.670.831.0 × 10–4
External validation dataset0.75 [0.55–0.88]0.66 ± 0.010.680.650.01
Subspace Discriminant (A2)
AUC *AccuracyPrecision **Recall **p ***
Training dataset0.79 [0.67–0.87]0.71 ± 0.010.690.830.02
External validation dataset0.71 [0.52–0.86]0.63 ± 0.020.680.650.046
Harmo CT features (B)
Linear SVM (B1)
AUCAccuracyPrecision **Recall **p ***
Training dataset0.77 [0.63–0.85]0.67 ± 0.020.740.581.0 × 10−4
External validation dataset0.56 [0.39–0.74]0.58 ± 0.010.670.520.5
Subspace Discriminant (B2)
AUCAccuracyPrecision **Recall **p ***
Training dataset0.76 [0.66–0.87]0.71 ± 0.020.730.60.01
External validation dataset0.57 [0.4–0.75]0.58 ± 0.010.670.520.50
Original CT features (C)
Linear SVM (C1)
AUCAccuracyPrecision **Recall **
Training dataset0.56 [0.42–0.68]0.52 ± 0.030.490.45
External validation dataset0.50 [0.34–0.68]0.43 ± 0.020.540.65
Subspace Discriminant (C2)
AUCAccuracyPrecision **Recall **
Training dataset0.63 [0.48–0.72]0.56 ± 0.030.580.56
External validation dataset0.51 [0.39–0.74]0.54 ± 0.010.580.65
PET features only (D)
Linear SVM (D1)
AUCAccuracyPrecision **Recall **p ***
Training dataset0.68 [0.53-0.78]0.64 ± 0.030.640.800.09
External validation dataset0.65 [0.43-0.82]0.64 ± 0.010.670.780.18
Subspace Discriminant (D2)
AUCAccuracyPrecision **Recall **p ***
Training dataset0.71 [0.59–0.82]0.69 ± 0.010.670.80.10
External validation dataset0.68 [0.51–0.84]0.60 ± 0.010.670.610.08
Harmo CT + Original PET + Clinical features (E)
Linear SVM (E1)
AUC *AccuracyPrecision **Recall **p ***
Training dataset0.79 [0.67–0.87]0.73 ± 0.020.720.836.0 × 10−5
External validation dataset0.73 [0.54–0.87]0.73 ± 0.010.770.740.02
Subspace Discriminant (E2)
AUC *AccuracyPrecision **Recall **p ***
Training dataset0.76 [0.65–0.86]0.74 ± 0.010.720.830.01
External validation dataset0.75 [0.54–0.88]0.68 ± 0.020.730.700.02
* AUCs in square brackets are their bootstrapped 95% CIs. ** Precision and recall are presented for class 1. *** p-values are calculated with respect to the conditions C1 and C2 for linear SVM and ESD models, respectively. Values in bold mean the statistical significance.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bertolini, M.; Trojani, V.; Botti, A.; Cucurachi, N.; Galaverni, M.; Cozzi, S.; Borghetti, P.; La Mattina, S.; Pastorello, E.; Avanzo, M.; et al. Novel Harmonization Method for Multi-Centric Radiomic Studies in Non-Small Cell Lung Cancer. Curr. Oncol. 2022, 29, 5179-5194. https://doi.org/10.3390/curroncol29080410

AMA Style

Bertolini M, Trojani V, Botti A, Cucurachi N, Galaverni M, Cozzi S, Borghetti P, La Mattina S, Pastorello E, Avanzo M, et al. Novel Harmonization Method for Multi-Centric Radiomic Studies in Non-Small Cell Lung Cancer. Current Oncology. 2022; 29(8):5179-5194. https://doi.org/10.3390/curroncol29080410

Chicago/Turabian Style

Bertolini, Marco, Valeria Trojani, Andrea Botti, Noemi Cucurachi, Marco Galaverni, Salvatore Cozzi, Paolo Borghetti, Salvatore La Mattina, Edoardo Pastorello, Michele Avanzo, and et al. 2022. "Novel Harmonization Method for Multi-Centric Radiomic Studies in Non-Small Cell Lung Cancer" Current Oncology 29, no. 8: 5179-5194. https://doi.org/10.3390/curroncol29080410

APA Style

Bertolini, M., Trojani, V., Botti, A., Cucurachi, N., Galaverni, M., Cozzi, S., Borghetti, P., La Mattina, S., Pastorello, E., Avanzo, M., Revelant, A., Sepulcri, M., Paronetto, C., Ursino, S., Malfatti, G., Giaj-Levra, N., Falcinelli, L., Iotti, C., Iori, M., & Ciammella, P. (2022). Novel Harmonization Method for Multi-Centric Radiomic Studies in Non-Small Cell Lung Cancer. Current Oncology, 29(8), 5179-5194. https://doi.org/10.3390/curroncol29080410

Article Metrics

Back to TopTop