Contrast Media & Molecular Imaging
Volume 2018 (2018), Article ID 3574310, 8 pages
https://doi.org/10.1155/2018/3574310

Clinical Study

Radiomic Profiling of Head and Neck Cancer: ¹⁸F-FDG PET Texture Analysis as Predictor of Patient Survival

G. Feliciani,¹ F. Fioroni,² E. Grassi,² M. Bertolini,² A. Rosca,³ G. Timon,³ M. Galaverni,³ C. Iotti,³ A. Versari,⁴ M. Iori,² and P. Ciammella³

¹Medical Physics Unit, Istituto Scientifico Romagnolo per lo Studio e la Cura dei Tumori (IRST) IRCCS, Meldola, Italy
²Medical Physics Unit, AUSL-IRCCS, Reggio Emilia, Italy
³Radiation Oncology Unit, AUSL-IRCCS, Reggio Emilia, Italy
⁴Nuclear Medicine Unit, AUSL-IRCCS, Reggio Emilia, Italy

Correspondence should be addressed to F. Fioroni; fioroni.federica@ausl.re.it

Received 6 April 2018; Accepted 26 August 2018; Published 27 September 2018

Academic Editor: Claudio Losio

Copyright © 2018 G. Feliciani et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Background and Purpose. The accurate prediction of prognosis and pattern of failure is crucial for optimizing treatment strategies for patients with cancer, and early evidence suggests that image texture analysis has great potential in predicting outcome both in terms of local control and treatment toxicity. The aim of this study was to assess the value of pretreatment ¹⁸F-FDG PET texture analysis for the prediction of treatment failure in primary head and neck squamous cell carcinoma (HNSCC) treated with concurrent chemoradiation therapy. Methods. We performed a retrospective analysis of 90 patients diagnosed with primary HNSCC treated between January 2010 and June 2017 with concurrent chemo-radiotherapy. All patients underwent ¹⁸F-FDG PET/CT before treatment. ¹⁸F-FDG PET/CT texture features of the whole primary tumor were measured using an open-source texture analysis package. Least absolute shrinkage and selection operator (LASSO) was employed to select the features that are associated the most with clinical outcome, as progression-free survival and overall survival. We performed a univariate and multivariate analysis between all the relevant texture parameters and local failure, adjusting for age, sex, smoking, primary tumor site, and primary tumor stage. Harrell c-index was employed to score the predictive power of the multivariate cox regression models. Results. Twenty patients (22.2%) developed local failure, whereas the remaining 70 (77.8%) achieved durable local control. Multivariate analysis revealed that one feature, defined as low-intensity long-run emphasis (LILRE), was a significant predictor of outcome regardless of clinical variables (hazard ratio < 0.001, ).The multivariate model based on imaging biomarkers resulted superior in predicting local failure with a c-index of 0.76 against 0.65 of the model based on clinical variables alone. Conclusion. LILRE, evaluated on pretreatment ¹⁸F-FDG PET/CT, is associated with higher local failure in patients with HNSCC treated with chemoradiotherapy. Using texture analysis in addition to clinical variables may be useful in predicting local control.

1. Introduction

Concurrent chemoradiotherapy (CRT) is the mainstay of treatment for early and locally advanced head and neck squamous cell carcinoma (HNSCC) [1, 2]. The accurate prediction of prognosis and failure in these patients is crucial for optimizing treatment. Some clinical features are commonly accepted as risk factors, such as tumor size, local anatomic invasion, nodal involvement, presence of distant metastases, and HPV status [3–5].

Nowadays, imaging plays a central role in the investigation of tumor prognosis. Radiological images are acquired as routine practice for almost every patient with HNSCC and represent an immense source of potential data for decoding tumor phenotypes and tumor heterogeneity [6–10]. Image texture is defined as a complex visual pattern within an image, consisting of simpler subpatterns with characteristic features, and texture analysis allows the mathematic detection of tumor heterogeneity. In the past years, CT texture analysis has been investigated in oncologic imaging for its ability to predict treatment outcome in patients with various tumors (lung, esophageal, and renal cancer) [11–13]. New and more descriptive metrics of tumor heterogeneity based on texture analysis (TA) of ¹⁸F-FDG PET images are now used to further improve outcome prediction [14–17] in different types of cancer [18, 19], including HNSCC.

Despite the promising results achieved in preliminary studies, TA still suffers from standardization issues, and the scientific community agreement is missing [20, 21]. The absolute values of the TA features used for stratifying patients depend on a number of variables, starting from image acquisition parameters to the contouring algorithm employed for tumor segmentation. Previous studies on phantoms and patients exploited these dependencies trying to find a common way to apply TA and standardize results [20, 22–25]. In the present study, we retrospectively apply TA to ¹⁸F-FDG PET/CT images of 90 patients affected by HNSCC and treated with CRT in order to identify imaging biomarkers able to predict patients’ outcomes such as overall survival (OS) and progression-free survival (PFS). We compared TA features prognostic values with standard clinical variables such as age, sex, size and site of primary tumors, clinical stage, and nodal involvement.

2. Materials and Methods

2.1. Patients

We enrolled 129 patients (median age 60 years, range 22–87 years) with diagnosis of HNSCC consecutively treated between January 2010 and June 2017 with concurrent CRT with curative intent at the Radiation Oncology Unit of our institution. The ethics committee of our hospital approved this retrospective study, and informed consent was obtained from all participants. All patients underwent ¹⁸F-FDG PET/CT before treatment for initial staging and were followed up for at least 5 years after treatment or until death. Patients’ medical records were retrospectively analyzed to extract information for outcome assessment such as details about therapy, last follow up date, disease status, pattern of recurrence, death and cause of death. Patients’ demographics and clinical characteristics, including age, sex, site of primary tumor, clinical stage, and type of treatment and are summarized in Table 1. Tumor sites were oral cavity, larynx, oropharynx, and hypopharynx in 5%, 20%, 60%, and 15% of patients, respectively.

Table 1: Patients’ clinical and demographic data.

Number of patients (%)

N
90

Mean age (range)
60 (22–87)

≥60 years
52 (58%)

<60 years
38 (42%)

Gender

Male
68 (75%)

Female
22 (15%)

Primary site

Oral cavity
4 (4%)

Oropharynx
49 (55%)

Hypopharynx
10 (10%)

Nasopharynx
13 (15%)

Larynx
14 (16%)

Clinical T stage

T1
15 (17%)

T2
30 (35%)

T3
29 (33%)

T4
13 (15%)

Clinical N stage

N0
5 (5%)

N1
21 (24%)

N2
59 (66%)

N3
5 (5%)

Clinical stage

III
32 (35%)

IV
58 (65%)

Patient cohort after follow-up (PFS)

NED
62 (69%)

RD
28 (31%)

Patient cohort after follow-up (OS)

Alive
65 (72%)

Dead
25 (28%)

2.2. Treatment Regimens and Follow-Up

All patients were uniformly treated with intensity-modulated radiation therapy (IMRT) and concurrent platinum-based chemotherapy. Weekly cisplatin was the primary choice of chemotherapeutic agent. Forty-three patients received neoadjuvant chemotherapy before CRT. The decision regarding the use of adjuvant chemotherapy was individualized, based on the extent of disease, medical conditions of the patient, and the radiation oncologist opinion. Most of the patients were treated with a simultaneous integrated boost (SIB) technique, with a total dose on gross tumor volume up to 66 Gy or equivalent in 30–33 fractions (mainly 66 Gy/2.2 Gy per fraction); the high-intermediate risk cervical lymphatic areas received 54–60 Gy or equivalent. Patients were followed after the conclusion of treatment to evaluate local control. All patients were followed clinically for at least 5 years after completion of CRT, every 1–3 months for the first two years, every 4–6 months for the next three years, and annually thereafter.

The follow-up evaluation included physical and endoscopic examinations jointly performed by radiation oncologist and otorhinolaryngologist. In addition, CT, MRI, and ¹⁸F-FDG PET/CT imaging were used to assess the clinical response and were usually performed within 6 months from the end of the treatment. Recurrence or distant metastasis were diagnosed based on either a positive biopsy or clinical/radiographic evidence of progression. The follow-up period was designated as the total time of follow-up, starting at treatment initiation and ending either at histologically confirmed local failure, or at radiologically systemic recurrence, or at last failure-free patient contact.

2.3. ¹⁸F-FDG PET/CT Imaging Protocol

Patients underwent imaging on a GE Discovery STE16 PET/CT scanner before treatment using a standard PET/CT clinical protocol. Patients fasted for at least 6 h before the intravenous administration of ¹⁸F-FDG (3.7 MBq/kg). Serum glucose concentrations were measured before FDG injection and were less than 150 mg/dl, or if between 150 and 200 mg/dl, the patient was hydrated and glucose concentration measured after 40–50 min. The scans were acquired with the patients immobilized in radiation treatment position using a dedicated flat table and personalized thermoplastic masks. PET images were corrected for random and scatter noise components and then reconstructed on a 256 × 256 image matrix using 3-D VUE Point HD algorithm (two iterations, 28 subsets, postfilter 6 mm) corrected for attenuation. Pixel spacing was of 2.73 × 2.73 mm with 3.27 mm slice thickness.

The qualitative and quantitative PET imaging evaluations were performed for each patient by two-blinded expert nuclear medicine physicians using the PET/CT fused images in transverse, coronal, and sagittal planes. A consensus was then reached by comparison of the two evaluations. The presence of abnormal FDG uptake, excluding the areas of physiologically increased uptake, was considered suspicious for malignancy. The readers had knowledge of all available clinical and imaging information related to the patients.

An expert nuclear medicine physician together with a medical physicist and a radiation oncologist, using a semiautomatic segmentation technique, outlined the primary tumor. The image segmentation of the whole primary lesion for each section was performed with a dedicated Advantage Workstation v 4.4 (GE Medical System, USA), using the 40% SUVmax isocontour algorithm to avoid operator dependence in contouring and any other manual adjustments. We assumed 40% threshold to acquire standardization with other works that perform prediction analysis [15]. Contours were then checked by visual inspection.

2.4. Texture Analysis

We extracted 75 features from each segmented tissue volume using an Open-Source Package “CGITA” version 1.3 [26]. Features were derived from data contained in the voxels of the segmented structures and can be grouped into different categories. First-order features were derived from the histogram of voxel intensities (SUVmean, SUVmax, skewness, kurtosis, etc.). Second-order textural features were based on matrices that contained information about the regional spatial arrangement of the voxels such as their homogeneity, contrast, and coarseness simulating the human perception of the image. Higher-order features such as Grey-level run length features focused on local collinear voxels with the same grey level.

The TA features analyzed are reported in Table 1 of reference [26]. A detailed description of the features is also reported in the supplemental material of the paper of Aerts et al. [27]. Images were digitized in 64 digitization bins according to the minimum and maximum values in the segmented volume. Tumor volumes smaller than 2.6 ml were excluded from the statistics due to insufficient number of voxels to perform texture analysis. This choice was based on previous evaluation about feature stability on NEMA IEC quality assurance phantom (2.6 ml corresponds to the volume of the 3rd sphere of the NEMA IEC phantom in the increasing order) [24, 28].

2.5. Statistical Analysis

First, associations of demographic and clinical characteristics of patients, such as sex, age, clinical stage (III vs IV), and tumor site (oropharynx vs other sites), with local control, progression-free survival (PFS), and overall survival (OS) were tested with univariate and multivariate analysis (Cox proportional hazard model). In particular, given the high heterogeneity of the population, we decided to create two tumor groups that are oropharynx and other sites (hypopharynx, larynx, oral cavity, and nasopharynx) in order to test its effect on survival. Furthermore, correlation between clinical parameters and predictive imaging biomarkers was investigated through spearman ρ.

Texture parameters were extracted and then compared in patients with local control against patients with local failure. To select both clinical and imaging variables that are more related to clinical outcomes, we used a “least absolute shrinkage and selection operator” (LASSO) method. LASSO can “shrink” the effect of unimportant features and can set their effects to zero together with removing redundancy among the features.

Given the iterative nature of LASSO algorithm, it was run one thousand times in order to have a statistics of the most descriptive clinical and imaging features. The most occurring features (more than 500 times over 1000 runs) were selected to build the final image biomarker-based multivariate Cox regression analysis. An internal 10-fold cross-validation algorithm was applied for validation. Finally, hazard ratios (HR) and confidence intervals (CI) were calculated. Discrimination, reflecting a correct ordering of the relative predictions for individuals, and the model’s ability to distinguish patients with local control against patients with local failure, were determined by the Harrell’s concordance-index (c-index) with R statistical package-based code (https://www.r-project.org/). c-index is an extension of ROC curves for multivariable models with 0.5 value indicating random discrimination and 1 perfect discrimination capabilities. The difference in prognostic value between the Cox model developed with or without imaging biomarkers was evaluated based on the comparison between Harrell’s concordance indexes. Finally, Kaplan–Meier curves have been calculated in order to show selected predictor stratification capabilities.

3. Results

At a median follow-up of 38 months (range 24–848 months), one hundred patients (79%) were alive, 13 died due to disease progression, and 12 due to other causes. Four patients were lost at follow-up.

Eighteen early-stage patients were excluded due to the small tumor volume (<2.6 ml) as textural features becomes unstable or even undefined due to small number of voxels involved; 13 patients due to technical problems regarding PET/CT acquisition (different acquisition, reconstruction protocols, or missing uptake data). The final number of evaluated patients for statistical analysis was 90. Of these 90 patients, 62 showed no evidence of disease (NED) at last follow-up; 28 patients (30%) had recurrent disease (RD), of which 21 patients developed locoregional recurrences and 7 isolated distant metastases.

The pattern of recurrence was analyzed in 21 patients who had locoregional relapse: “in-field” (if >80% of the tumor recurrence resided within the prescription 95% isodose surface) in 20 patients (95%), at RT field margin (if 20–80% of the lesion was inside the 95% isodose surface) in no patients (0%), and “out-field” in one patient (5%).

The estimated 2- and 5-year OS rates were 86% (95% CI 79–92%) and 68 % (95% CI 58–78%), respectively. The corresponding 2- and 5-year PFS rates were 70% (95% CI 58–77%) and 66% (95% CI 55–75%), respectively.

Results of univariate analysis are showed in Table 2. No standard clinical and demographics parameters showed a statistically significant correlation with PFS. Gender and age showed significant correlation for OS. In particular, gender seems to have a strong impact on survival in our dataset (HR = 8.316). Clinical stage was correlated with both PFS and OS with an HR of 1.76 and 1.92, respectively, without any statistical significance.

Table 2: Univariate model of patients’ clinical data with progression-free survival (PFS) and overall survival (OS). The number of patients is 90 for both PFS and OS, whereas the number of events is 28 and 25, respectively.

value
HR
95% CI for HR

Lower
Upper

Univariate Cox regression for PFS

Gender (F vs M)
0.698
1.196
0.485
2.951

Age
0.099
1.027
0.995
1.060

Stage (III vs IV)
0.197
1.758
0.747
4.137

Tumor site (oroph vs other)
0.914
1.043
0.488
2.227

CHT (No vs Yes)
0.101
0.470
0.190
1.159

Univariate Cox regression for OS

Gender (F vs M)
0.038
8.495
1.149
62.820

Age
0.004
1.052
1.016
1.089

Stage (III vs IV)
0.169
2.010
0.797
5.070

Tumor site (oropharynx vs others)
0.772
0.886
0.391
2.009

CHT (No vs Yes)
0.838
0.895
0.307
2.608

At multivariate analysis, shown in Table 3(A), no standard clinical parameters were correlated with PFS, and the c-index resulted in 0.65 in the ability of prediction of the model, without statistical significance. Regarding OS, the multivariable Cox model showed a statistically significant correlation with both gender and age with a c-index of 0.73 as can be seen in Table 4(A). It is noticeable that in any Cox regression, tumor site seems to have no effects on survival prediction for our population. We can then state that this parameter will not influence our further models with imaging biomarkers.

Table 3: Multivariate analysis of the patients’ dataset for PFS without (A) and including (B) imaging biomarkers. Harrel c-indexes of the models, which score their prognostic power, are 0.65 and 0.76, respectively. The two c-indexes are significantly different with a value of 0.01.

value
HR
95% CI for HR

Lower
Upper

(A)

Gender (F vs M)
0.713
1.191
0.469
3.026

Age
0.426
1.014
0.979
1.050

Stage (III vs IV)
0.160
1.894
0.777
4.617

Tumor site (oropharynx vs others)
0.902
0.953
0.445
2.042

CHT (No vs Yes)
0.176
0.489
0.173
1.378

(B)

Age
0.254
1.02E + 00
9.86E − 01
1.06E + 00

Stage (III vs IV)
0.440
1.43E + 00
5.77E − 01
3.55E + 00

CHT (No vs Yes)
0.012
2.03E − 01
5.84E − 02
7.06E − 01

Run percentage
0.176
1.97E − 01
1.87E − 02
2.07E + 00

LILRE
<0.001
1.41E − 83
4.05E − 128
4.90E − 39

Coarseness
0.970
4.12E − 01
2.32E − 21
7.33E + 19

Code similarity
0.129
4.27E − 23
5.77E − 52
3.16E + 06

Table 4: Multivariate analysis of the patients dataset for OS without (A) and including (B) imaging biomarkers.

value
HR
95% CI for HR

Lower
Upper

(A)

Gender (F vs M)
0.042
8.036
1.082
59.770

Age
0.017
1.057
1.012
1.109

Stage (III vs IV)
0.262
1.734
0.663
4.535

OvsL
0.895
0.937
0.395
2.029

Novscht
0.317
1.916
0.521
7.051

(B)

Gender (F vs M)
0.090
5.75E + 00
7.90E − 01
4.81E + 01

Age
0.011
1.05E + 00
1.01E + 00
1.09E + 00

Stage (III vs IV)
0.341
1.59E + 00
5.82E − 01
3.83E + 00

LILRE
0.635
7.08E − 05
9.72E − 25
8.35E + 9

SUL peak
0.068
1.10E + 00
1.009E + 00
1.192E + 00

The LASSO algorithm revealed that four parameters, including run percentage, low-intensity long-run emphasis (LILRE), coarseness, and code similarity, showed significant differences between local failure and local control groups, and these parameters were selected, together with two clinical variables (age and stage), to build a multivariate Cox regression for PFS. No significant correlation exists between the analyzed parameters as spearman rho coefficient is always below 0.6.

In Table 3(B) the resulting model for PFS is shown; the Harrel c-index is of 0.76 with a statistical significance of at cross-validation. Also the comparison of the two c-indexes underlined a difference between the two models with statistical significance and . In particular LILRE, which is the strongest predictor in the model, is associated with the presence of long string of low grey level pixel in the tumor contours. In Figure 1 the values of the feature in RD and NED patients are shown. An evident difference in median values appears between the populations. The Wilcoxon Test scored a value of 0.006 between RD and NED patients. Furthermore, in Figure 2, the difference in a 2D render of the morphology of the 2 populations (representative patients) is revealed. In RD patients, it is evident how the borders of the tumors are sharper in comparison to NED patients. Finally, in Figure 3(a), Kaplan–Meier survival curves for LILRE, where survival curves are split by the median value of LILRE, and in Figure 3(b), the effect of the use of chemotherapy (CHT) are shown. Pearson correlation coefficient ρ calculated between LILRE and the other clinical variables was always below 0.2, confirming the independence of this imaging biomarker.

Figure 1: Low-intensity long-run emphasis values calculated on different populations in arbitrary unit of measure: recurrent disease (RD) and no evidence of disease (NED) patients.

(a)

(b)

Figure 2: Detail of the uptake of 2D slices centered in the middle of 2 tumor VOIs. On the left, the slice is taken from a RD patient and on the right from a NED patient. It can be appreciated that a lower value in LILRE produces a more speculated uptake in the tumor, whereas in NED patient, the uptake is more Gaussian like. (a) RD patient. (b) NED patient.

(a)

(b)

Figure 3: Kaplan–Meier survival curves for the two statistically significant predictors emerged from the multivariate Cox regression. In (a), the median value of low-intensity long-run emphasis feature is employed to split the curves of the patients, whereas in (b), the split is performed by using chemotherapy usage (Yes or No).

By adding imaging biomarkers and processing variables in the same way of PFS dataset, we end up with the model shown in Table 4(B). This model, which includes five variables (gender, age, staging, low-intensity large-zone emphasis, and SUL peak), achieves a c-index of 0.76, but it is not significantly different from the model with only the clinical variables.

4. Discussion

To date, there is increased evidence suggesting that genomic heterogeneity of aggressive tumors could translate into intratumoral spatial heterogeneity, which can be represented on anatomical and functional scales [18, 27, 29, 30].

This is the central idea that underlies radiomics, in which large amounts of information through advanced quantitative analysis of radiological images are used as noninvasive means to characterize intratumoral heterogeneity and to identify important prognostic features of cancer [27, 31–34]. Nowadays, medical imaging plays a central role in the investigation of intratumoural heterogeneity, as radiological images are acquired as routine practice for almost every patient with cancer. Medical images such as ¹⁸F-FDG PET/CT and CT are minimally invasive and include an immense source of potential data for decoding tumour phenotypes.

Texture feature-based analysis of clinical outcomes is actively being investigated as a prognostic tool in radiation oncology in order to identify potentially predictive radiomic biomarkers for either clinical outcomes (local control and/or treatment toxicity) or distinguishing local recurrence from radiation-induced injury. Some studies have also been performed to identify radiomic signatures for HNSCC. Aerts et al. analyzed radiomic values of CTs from 1,019 patients suffering from non-small-cell lung cancer or head and neck cancer in relation to prognostic features and showed that combining radiomic signature with TNM staging improved the predictive power in all groups of patients [27]. Parmar et al. investigated 440 radiomic features extracted from the CTs of 878 lung cancers and HNSCC patients and showed that, at multivariate analysis, the radiomic feature clusters highly correlate with tumor stage and moderately correlate with HPV status [34]. In an additional study, the same group established a reliable machine-learning method for the prediction of OS in HNSCC patients based on CT scan-derived radiomic features. Here, training was performed on 440 radiomic features among a cohort of 101 HNSCC patients, while another 95 patients served as the validation cohort. In this study, authors showed that machine-learning methods had a high predictive power with a good stability and believed that this technique could improve the application of radiomics in cancer [35]. These 3 publications all suggest a relevance of radiomics in HNSCC and a potential future role in HNSCC classification and treatment to improve clinical decision-making. Recently, it was reported that a 4-feature-based radiomic signature showed strong correlation with survival and distant metastasis in oropharyngeal squamous cell carcinoma [36]. The concordance index (CI) in training cohorts and validation cohorts were ranging from 0.55 to 0.69.

Vallières et al. recently investigated the possibility of constructing a prediction model integrating clinical information with radiomics analysis, using advanced machine learning, in patients with HNSCC treated with chemoradiation. 1,615 different radiomic features were extracted from PET and CT pretreatment images of 300 patients, and the results demonstrated the potential of radiomics for assessing the risk of specific tumor outcomes using multiple stratification groups [10].

In this study, we analyzed a cohort of 90 patients with diagnosis of HNSCC using a multivariate Cox model comparing features based on clinical variables and imaging biomarkers. Univariate analysis showed how gender and age are important predictors for OS. The effect of gender on overall survival is a known phenomenon, as shown also recently in the large European multicenter study [37]. This trial reported that a significant effect on OPC survival was apparent for female sex (aHR: 0.50; 95% CI: 0.29–0.85).

Models including imaging biomarkers were always superior to those with only clinical variables, even if only in the case of PFS, we had a statistically significant difference. The multivariable model developed with clinical variables alone scored a c-index of 0.65, whereas including imaging biomarkers, a c-index of 0.76 was reached in predicting patient outcome as PFS. Furthermore, the comparison of the c-index underlined the significance of the difference between the two models. This improvement in patient stratification can be addressed to the introduction of the textural feature LILRE which is strongly associated with patient disease recurrence, as shown in the box plot in Figure 1 and in Kaplan–Meier curves in Figure 3(a) employing the median value of LILRE over the entire population as split value. It is also interesting to notice that this feature can be associated to morphological properties of the image itself. These findings underline the importance of further studies to validate the use of imaging biomarkers to improve patients’ stratification and prognostic accuracy.

Our study has several limitations, and should be considered as a preliminary experience with need for possible methodological and technical refinements, beyond independent external validation. It is a retrospective study with a relatively small population. The characteristics of enrolled patients were heterogeneous, as well as disease sites, and the selection of treatment might be biased, so that the significance of features can be underestimated.

Despite this, it seems that heterogeneity is not influencing the results as no standard Cox regression models gives any statistically significant association with survival in this cohort. Furthermore, no correlation was found between our main predictor (LILRE) and tumor site. Additionally, HPV status was not determined in all patients, and it has not been considered as a predictive factor. Although especially in oropharyngeal cancer infection with HPV has a predictive value for treatment response [4, 38], no different radiomic features were identified between HPV-positive and HPV-negative patients.

In the other two analyses, statistically significant differences in some texture features between HPV+ and HPV− HNSCC were found, even if a single radiomic feature may not have enough predictive power [39, 40].

Despite the limitations of this study, our results appear to be promising. In order to confirm these findings, external validation will be performed, and for the case of oropharynx, tumor HPV infection data will be included.

5. Conclusion

We compared a clinical variables alone statistical model with another one which includes imaging biomarkers. Superiority of the imaging biomarkers model was proven in the case of PFS. Further studies will be aimed at confirming our findings with an external dataset and to demonstrate the superiority of models comprising imaging biomarkers also for overall survival.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

J. Setton, N. Caria, J. Romanyshyn et al., “Intensity-modulated radiotherapy in the treatment of oropharyngeal cancer: an update of the memorial sloan-kettering cancer center experience,” International Journal of Radiation Oncology∗Biology∗Physics, vol. 82, no. 1, pp. 291–298, 2012.

J. T. Parsons, W. M. Mendenhall, S. P. Stringer et al., “Squamous cell carcinoma of the oropharynx,” Cancer, vol. 94, no. 11, pp. 2967–2980, 2002.

K. K. Ang, J. Harris, R. Wheeler et al., “Human papillomavirus and survival of patients with oropharyngeal cancer,” New England Journal of Medicine, vol. 363, no. 1, pp. 24–35, 2010.

C. Fakhry, W. H. Westra, S. Li et al., “Improved survival of patients with human papillomavirus-positive head and neck squamous cell carcinoma in a prospective clinical trial,” JNCI Journal of the National Cancer Institute, vol. 100, no. 4, pp. 261–269, 2008.

I. Martinez, J. Wang, K. F. Hobson, R. L. Ferris, and S. A. Khan, “Identification of differentially expressed genes in HPV-positive and HPV-negative oropharyngeal squamous cell carcinomas,” European Journal of Cancer, vol. 43, no. 2, pp. 415–432, 2007.

R. Lim, A. Eaton, N. Y. Lee et al., “¹⁸F-FDG PET/CT metabolic tumor volume and total lesion glycolysis predict outcome in oropharyngeal squamous cell carcinoma,” Journal of Nuclear Medicine, vol. 53, no. 10, pp. 1506–1513, 2012.

N.-M. Cheng, J. T.-C. Chang, C.-G. Huang et al., “Prognostic value of pretreatment ¹⁸F-FDG PET/CT and human papillomavirus type 16 testing in locally advanced oropharyngeal squamous cell carcinoma,” European Journal of Nuclear Medicine and Molecular Imaging, vol. 39, no. 11, pp. 1673–1684, 2012.

M. L. Belli, M. Mori, S. Broggi et al., “Quantifying the robustness of [¹⁸F] FDG-PET/CT radiomic features with respect to tumor delineation in head and neck and pancreatic cancer patients,” Physica Medica, vol. 49, pp. 105–111, 2018.

M. Picchio, M. Kirienko, P. Mapelli et al., “Predictive value of pre-therapy ¹⁸F-FDG PET/CT for the outcome of ¹⁸F-FDG PET-guided radiotherapy in patients with head and neck cancer,” European Journal of Nuclear Medicine and Molecular Imaging, vol. 41, no. 1, pp. 21–31, 2014.

M. Vallières, E. Kay-rivest, L. J. Perrin et al., “Radiomics strategies for risk assessment of tumour failure in head-and-neck cancer,” Scientific Reports, vol. 7, no. 1, 2017.

R. J. Gillies, P. E. Kinahan, and H. Hricak, “Radiomics: images are more than pictures, they are data,” Radiology, vol. 278, no. 2, pp. 563–577, 2016.

B. Ganeshan, S. Abaleke, R. C. D. Young, C. R. Chatwin, and K. A. Miles, “Texture analysis of non-small cell lung cancer on unenhanced computed tomography: initial evidence for a relationship with tumour glucose metabolism and stage,” Cancer Imaging, vol. 10, pp. 137–143, 2010.

B. Ganeshan and K. A. Miles, “Quantifying tumour heterogeneity with CT,” Cancer Imaging, vol. 13, no. 1, pp. 140–149, 2013.

S. Chicklore, V. Goh, M. Siddique, A. Roy, P. K. Marsden, and G. J. R. Cook, “Quantifying tumour heterogeneity in ¹⁸F-FDG PET/CT imaging by texture analysis,” European Journal of Nuclear Medicine and Molecular Imaging, vol. 40, no. 1, pp. 133–140, 2013.

N. Cheng, Y. D. Fang, L. Lee et al., “Zone-size nonuniformity of ¹⁸F-FDG PET regional textural features predicts survival in patients with oropharyngeal cancer,” European Journal of Nuclear Medicine and Molecular Imaging, vol. 42, no. 3, pp. 419–428, 2015.

N.-M. Cheng, Y.-H. D. Fang, J. T.-C. Chang et al., “Textural features of pretreatment ¹⁸F-FDG PET/CT images: prognostic significance in patients with advanced T-stage oropharyngeal squamous cell carcinoma,” Journal of Nuclear Medicine, vol. 54, no. 10, pp. 1703–1709, 2013.

G. J. R. Cook, C. Yip, M. Siddique et al., “Are pretreatment ¹⁸F-FDG PET tumor textural features in non-small cell lung cancer associated with response and survival after chemoradiotherapy?” Journal of Nuclear Medicine, vol. 54, no. 1, pp. 19–26, 2013.

I. El Naqa, P. Grigsby, A. Apte et al., “Exploring feature-based approaches in PET images for predicting cancer treatment outcomes,” Pattern Recognition, vol. 42, no. 6, pp. 1162–1171, 2009.

M. Hatt, M. Majdoub, M. Vallières et al., “¹⁸F-FDG PET uptake characterization through texture analysis: investigating the complementary nature of heterogeneity and functional tumor volume in a multi-cancer site patient cohort,” Journal of Nuclear Medicine, vol. 56, no. 1, pp. 38–44, 2015.

R. T. H. Leijenaar, G. Nalbantov, S. Carvalho et al., “The effect of SUV discretization in quantitative FDG-PET radiomics: the need for standardized methodology in tumor texture analysis,” Scientific Reports, vol. 5, no. 1, 2015.

A. Chalkidou, M. J. O’Doherty, and P. K. Marsden, “False discovery rates in PET and CT studies with texture features: a systematic review,” PLoS One, vol. 10, no. 5, Article ID e0124165, 2015.

P. E. Galavis, C. Hollensen, N. Jallow, B. Paliwal, and R. Jeraj, “Variability of textural features in FDG PET images due to different acquisition modes and reconstruction parameters,” Acta Oncologica, vol. 49, no. 7, pp. 1012–1016, 2010.

R. T. H. Leijenaar, S. Carvalho, E. R. Velazquez et al., “Stability of FDG-PET radiomics features: an integrated analysis of test-retest and inter-observer variability,” Acta Oncologica, vol. 52, no. 7, pp. 1391–1397, 2013.

M. J. Nyflot, F. Yang, D. Byrd, S. R. Bowen, G. A. Sandison, and P. E. Kinahan, “Quantitative radiomics: impact of stochastic effects on textural feature analysis implies the need for standards,” Journal of Medical Imaging, vol. 2, no. 4, Article ID 41002, 2015.

G. Doumou, M. Siddique, C. Tsoumpas, V. Goh, and G. J. Cook, “The precision of textural analysis in ¹⁸F-FDG-PET scans of oesophageal cancer,” European Radiology, vol. 25, no. 9, pp. 2805–2812, 2015.

Y. H. D. Fang, C. Y. Lin, M. J. Shih et al., “Development and evaluation of an open-source software package “cGITA” for quantifying tumor heterogeneity with molecular images,” BioMed Research International, vol. 2014, Article ID 248505, 9 pages, 2014.

H. J. W. L. Aerts, E. R. Velazquez, R. T. H. Leijenaar et al., “Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach,” Nature Communications, vol. 5, no. 1, 2014.

I. Shiri, A. Rahmim, P. Ghaffarian, P. Geramifar, H. Abdollahi, and B.-R. Ahmad, “The impact of image reconstruction settings on ¹⁸F-FDG PET radiomic features : multi-scanner phantom and patient studies,” European Radiology, vol. 27, no. 11, pp. 4498–4509, 2017.

E. Segal, C. B. Sirlin, C. Ooi et al., “Decoding global gene expression programs in liver cancer by noninvasive imaging,” Nature Biotechnology, vol. 25, no. 6, pp. 675–680, 2007.

M. Diehn, C. Nardini, D. S. Wang et al., “Identification of noninvasive imaging surrogates for brain tumor gene-expression modules,” Proceedings of the National Academy of Sciences, vol. 105, no. 13, pp. 5213–5218, 2008.

R. J. Gillies, A. R. Anderson, R. A. Gatenby, and D. L. Morse, “The biology underlying molecular imaging in oncology: from genome to anatome and back again,” Clinical Radiology, vol. 65, no. 7, pp. 517–521, 2010.

P. Lambin, E. Rios-Velazquez, R. Leijenaar et al., “Radiomics: extracting more information from medical images using advanced feature analysis,” European Journal of Cancer, vol. 48, no. 4, pp. 441–446, 2012.

V. Kumar, Y. Gu, S. Basu et al., “Radiomics: the process and the challenges,” Magnetic Resonance Imaging, vol. 30, no. 9, pp. 1234–1248, 2012.

C. Parmar, R. T. H. Leijenaar, P. Grossmann et al., “Radiomic feature clusters and prognostic signatures specific for lung and head & neck cancer,” Scientific Reports, vol. 5, no. 1, 2015.

C. Parmar, P. Grossmann, D. Rietveld, M. M. Rietbergen, P. Lambin, and H. W. L. Aerts, “Radiomic machine-learning classifiers for prognostic biomarkers of head and neck cancer,” Frontiers in Oncology, vol. 5, 2015.

R. T. H. Leijenaar, S. Carvalho, F. J. P. Hoebers et al., “External validation of a prognostic CT-based radiomic signature in oropharyngeal squamous cell carcinoma,” Acta Oncologica, vol. 54, no. 9, pp. 1423–1429, 2015.

D. Anantharaman, A. Billot, T. Waterboer et al., “Predictors of oropharyngeal cancer survival in Europe,” Oral Oncology, vol. 81, pp. 89–94, 2018.

C. Fakhry, Q. Zhang, P. F. Nguyen-Tan et al., “Human papillomavirus and overall survival after progression of oropharyngeal squamous cell carcinoma,” Journal of Clinical Oncology, vol. 32, no. 30, pp. 3365–3373, 2014.

A. Fujita, K. Buch, B. Li, Y. Kawashima, M. M. Qureshi, and O. Sakai, “Difference between HPV-positive and HPV-negative non-oropharyngeal head and neck cancer: texture analysis features on CT,” Journal of Computer Assisted Tomography, vol. 40, no. 1, pp. 43–47, 2016.

K. Buch, A. Fujita, B. Li et al., “Using texture analysis to determine human papillomavirus status of oropharyngeal squamous cell carcinomas on CT,” American Journal of Neuroradiology, vol. 36, no. 7, pp. 1343–1348, 2015.