Skip to main content

MRI-based radiomic prognostic signature for locally advanced oral cavity squamous cell carcinoma: development, testing and comparison with genomic prognostic signatures



. At present, the prognostic prediction in advanced oral cavity squamous cell carcinoma (OCSCC) is based on the tumor-node-metastasis (TNM) staging system, and the most used imaging modality in these patients is magnetic resonance image (MRI). With the aim to improve the prediction, we developed an MRI-based radiomic signature as a prognostic marker for overall survival (OS) in OCSCC patients and compared it with published gene expression signatures for prognosis of OS in head and neck cancer patients, replicated herein on our OCSCC dataset.


For each patient, 1072 radiomic features were extracted from T1 and T2-weighted MRI (T1w and T2w). Features selection was performed, and an optimal set of five of them was used to fit a Cox proportional hazard regression model for OS. The radiomic signature was developed on a multi-centric locally advanced OCSCC retrospective dataset (n = 123) and validated on a prospective cohort (n = 108).


The performance of the signature was evaluated in terms of C-index (0.68 (IQR 0.66–0.70)), hazard ratio (HR 2.64 (95% CI 1.62–4.31)), and high/low risk group stratification (log-rank p < 0.001, Kaplan-Meier curves). When tested on a multi-centric prospective cohort (n = 108), the signature had a C-index of 0.62 (IQR 0.58–0.64) and outperformed the clinical and pathologic TNM stage and six out of seven gene expression prognostic signatures. In addition, the significant difference of the radiomic signature between stages III and IVa/b in patients receiving surgery suggests a potential association of MRI features with the pathologic stage.


Overall, the present study suggests that MRI signatures, containing non-invasive and cost-effective remarkable information, could be exploited as prognostic tools.


Oral cavity squamous cell carcinomas (OCSCCs), including tongue cancers, are the most common malignancies of the oral cavity, accounting for approximately half of all head and neck squamous cell carcinomas (HNSCCs) [1].

The treatment mainstay for loco-regionally advanced (i.e., clinical tumor-node-metastasis (cTNM) stage III/IVa-b according to the 8th edition of the AJCC/UICC staging system) OCSCC is surgery followed by adjuvant radiation, plus concomitant chemotherapy in case of adverse pathologic factors (i.e., residual tumor after resection and/or extracapsular spread to regional lymph nodes). In this scenario, the strongest prognostic factor currently available is pathologic tumor-node (pTNM) metastasis stage.

However, the pTNM is available only after performing surgery, and at diagnosis no robust baseline prognostic factors are available, with the only exception of the cTNM, that is obtained through clinical and radiological assessments. In this setting, magnetic resonance imaging (MRI) provides a high-quality resolution in determining soft tissue infiltration, so it is one of the most relevant tools to guide surgical planning. Nevertheless, so far images per se have not been used in the clinical practice for prognostic purposes, apart from determining the cTNM.

Thus, additional biomarkers are needed to better stratify patients, refine the stage-based clinical decision and lead to more personalized therapies. In this context, high throughput “omics” technologies have recently gained increasing interest for the identification of prognostic factors [2]. Among those, genomics and radiomics are the ones that have been more extensively investigated as far as the development of prognostic models is concerned [3,4,5,6,7,8].

Starting from early 2000, omics techniques including different comprehensive analyses on gene expression (transcriptomics), proteins and metabolites, were applied to identify biomarkers (see glossary in a very recent radiomic review [9]). In a review considering publications up to 2019 about transcriptomics and epigenomics in HNSCC, we reported 33 studies analyzing gene expression and disclosing biological and/or prognostic roles in the disease [10]. Thereafter, thanks to the availability of data from The Cancer Genome Atlas [11], gene expression signatures have been generated including all HNSCC anatomical subsites or specifically OCSCC [10, 12, 13].

Radiomics, namely the extraction and mining of quantitative features from radiological imaging, has the potential to provide non-invasive prognostic biomarkers. To date, several studies applied radiomics to characterize tumors, predict pathological features and predict clinical outcomes in HNSCC patients, as extensively reviewed in [14,15,16]. Most of the radiomic studies focusing on prognosis of overall survival (OS) in HNSCC patients were based on computer tomography or positron emission tomography [6, 7, 17,18,19,20,21,22,23,24,25,26,27,28,29,30,31], while only few prognostic models with features extracted from MRI were developed [32,33,34,35,36,37,38,39,40,41,42]. In particular, to the best of the authors’ knowledge, only three studies [32, 34, 42] developed an MRI-radiomic signature for prognosis of OS specific for OCSCC patients. Few other MRI-based radiomic studies focused on oral cavity cancer were proposed to predict tumor grading [43], pathological differentiation [44] and extracapsular nodal spread [45]. In this scenario, given the high heterogeneity of HNSCC patients [46] and the need for more personalized treatments, the development of MRI-based radiomic signature specific for OCSCC patients is fundamental.

To this aim we exploited our recently obtained database BD2Decide [47] that could enable linking together rigorously annotated patient-specific multiparameter clinical, pathologic, demographic, transcriptomics and radiomics data from the currently largest cohort of patients with locoregionally advanced HNSCC.

The purpose of this study was to develop a prognostic MRI-based radiomic signature from the retrospective OCSCC cohort, and to test it on the prospective one. Finally, the prognostic effect of the MRI-based radiomic model was compared to conventional prognostic metrics (cTNM and pTNM) and to some published gene expression prognostic signatures [48,49,50,51,52,53,54].

Materials and methods

Analyzed dataset

A subset of the BD2Decide project (NCT02832102) patients with positive pathologic diagnosis of OCSCC and loco-regionally advanced disease (cTNM III, IVa or IVb according to the 8th edition of AJCC/UICC) treated with curative intent [47] was used for this study. The protocols were approved by the Ethical Committees of the participating centers, and data acquisition followed the General Data Protection Regulation of the EU. All patients signed the informed consent. The inclusion criteria of the study were the following: (i) availability of T1-weighted (T1w) and T2-weighted (T2w) MR image sequence and (ii) images acquired with 1.5 T scanner (Fig. 1). Images from the selected patients were collected from four different clinical centers: the Fondazione IRCCS Istituto Nazionale dei Tumori of Milan, Italy; the Azienda Ospedaliero Universitaria di Parma, Italy; the Spedali Civili di Brescia, Italy; the Amsterdam VU Medisch Centrum Medical Center, the Netherlands.

Fig. 1
figure 1

CONSORT flow diagram

Clinical endpoint

The clinical endpoint analyzed in this study was OS, defined as the time between the primary tumor diagnosis and the day of death or last follow-up. In details, OS times were calculated in months from the date of diagnosis to the date of death of any cause (event), and censored at the date of last follow-up for patients that were still alive.

MRI image acquisition

For each patient T1w and T2w MRI images were acquired using scanners with a field strength of 1.5 T. The T1w and T2w images were acquired using a turbo spin-echo pulse sequence. Other image acquisition parameters, such as time of repetition, time of echo, pixel spacing, slice thickness, were not standardized. Table 1 gives an overview of the image acquisition parameters used to acquire the images in the study.

Table 1 Image acquisition parameters for the patients of the BD2Decide dataset

Image segmentation

The gross tumor volume was segmented at the clinical centers using a semi-automatic segmentation software based on coupled shape modeling [55]. The segmentation of the region of interest (ROI), corresponding to the primary tumor, was performed manually slice by slice by expert radiologists (one for each center) dedicated to head and neck cancers. The tumor boundary was delineated by considering T2w sequence as reference, and checked and corrected based on the T1w sequence. An example of segmented images is displayed in Fig. 2.

Fig. 2
figure 2

Type of magnetic resonance images acquired for the study: (A) T1-weighted image; (B) T2-weighted image. The segmented region of interest is also displayed

Image preprocessing

Image preprocessing was applied to the MRI images to reduce all the imaging-related sources of variability. Specifically, four steps of image preprocessing were applied [56]: (i) a 3D Gaussian filter with a 3 × 3 × 3 voxel kernel and σ = 0.5 was used to denoise the images; (ii) the N4ITK algorithm [57] was used for the correction of intensity non-uniformities due to local variations of the magnetic field; (iii) intensity standardization was performed using Z-score to ensure that each MRI image had similar ranges of signals; (iv) voxel size resampling to an isotropic resolution of 2 mm was performed with B-spline interpolation [58].

Radiomic features extraction

The extraction of radiomic features was performed using Pyradiomics 2.2.0 (open-source, available at and run on Python) [59]. A total of 1072 radiomic features, 536 per image type (T1w, T2w) were extracted. The features belonged to different categories: shape and size (14 features), first order statistics (18 features), textural (40 features), wavelet (464 features). Textural features were computed using the grey level co-occurrence matrix (GLCM) and the grey level run length matrix (GLRLM). The full list of radiomic features is available in Pyradiomics documentation [60]. A fixed-bin histogram discretization (32 bins) was used prior to features extraction.

Radiomic features postprocessing and radiomic model development

Features selection and survival model training was performed on the 123 retrospective patients (training dataset). First, features were Z-score normalized to ensure comparable ranges for the feature values. The mean and standard deviation used for the feature-wise normalization of the training set were then applied to normalize the features of the prospective set. The normalization was performed to improve the convergence of the optimization algorithms used during model fitting [61]. Then, as the number of extracted features was much larger than the number of available patients, a process of features selection was required to reduce the dimensionality of the radiomic dataset and avoid overfitting.

The features selection process was performed on the training set and comprised (i) stability analysis and (ii) supervised feature selection based on both univariate and multivariate Cox regression, similarly to a previous MRI-based radiomic study [62]. As regards the first step of the features selection, the stability of the features to small translations of the ROI (i.e., ± 10% of the length of the bounding box surrounding the ROI along both the x and y directions), used to mimic the effect of inter-reader variability in the segmentation, was evaluated as described elsewhere [56, 62, 63]. Features with intra-class correlation coefficient (ICC) above 0.75 were considered stable and selected. The second step of features selection consisted in the evaluation of features prognostic significance and Harrel’s concordance index (C-index) [64] in both univariate and multivariate Cox regression. Specifically, 100 bootstrap iterations were performed in which the training dataset was divided in independent subsets of training and validation (70% and 30%, respectively) and, at each iteration, potentially different feature sets can be found. For each iteration, univariate Cox models were used to select the features that were significantly prognostic for OS (p < 0.05). The 20 features associated with lower p in the univariate analysis were combined in a multivariate Cox model using a wrapper forward feature selection algorithm [62, 65]. The combination of features associated with the highest C-index was the selected feature set of the specific bootstrap iteration. At the end of the 100 bootstrap iterations the 5 most selected features were considered. The overall process (i.e., second step of features selection) was repeated 20 times with different seeds of the random number generator to increase the robustness, and the final features set contained the 5 most selected features over the 20 repetitions. Figure 3 schematically describes the 2-step process of features selection.

Fig. 3
figure 3

Features selection process

The final feature set was used to train a multivariate Cox model for OS on the retrospective set. The radiomic signature was thus obtained for each patient as the linear combination of the features and the corresponding regression coefficients. The median value of the signature in the training set was used as a threshold to classify the patients into high and low risk groups. Patients with a signature higher than a threshold are classified as high-risk patients, while patients with signature below the threshold are classified as low-risk patients.

All the steps of features postprocessing, feature selection and model training were performed in Matlab 2022a (Mathworks, Natick, MA, USA).

Model testing and comparison

The independent prospective set comprised 108 patients. First the Z-score normalization of the features of the prospective set was performed, based on the mean and standard deviation of Z-score normalization estimated from the retrospective set. The radiomic signature of the prospective set was computed as linear combination of the features and the corresponding regression coefficients of the trained Cox model, and high/low risk classification was based on the threshold estimated from the training set (i.e., median value of the signature on the training set).

The radiomic signature prognostic performance was evaluated on both the retrospective (used for model training) and prospective (used for model testing) sets. First, C-index between the signature and the OS [64] was computed through 100 bootstrap iterations to obtain a distribution of C-indexes. Second, the hazard ratio (HR) was computed. Third, the p-value of the log-rank test [66] comparing the Kaplan-Meier curves [67] for high and low risk groups was evaluated.

The prognostic performance of the radiomic signature was compared with that of the cTNM and the pTNM on the prospective subset 1 (n = 75 patients as shown in Fig. 1), and with that of 7 published prognostic gene expression signatures (detailed below) on the prospective subset 2 (n = 90 patients as shown in Fig. 1). Finally, the relationship of the radiomic signature with the pTNM was also evaluated, by analyzing the distribution of the radiomic signature at different pTNM stages in the dataset of 177 patients (i.e., retrospective and prospective patients for whom the pTNM was available, Fig. 1).

HNSCC/OCSCC gene expression signatures

Formalin-fixed, paraffin-embedded samples of prospective subset 2 (n = 90) were collected, and after histopathological revision, selected tumor areas were manually macro-dissected. RNA extraction was performed using the Qiagen RNeasy Mini Kit with the QIAcube robotic station (Qiagen, Düsseldorf, Germany), according to the manufacturer’s recommendations. Quantification and quality check were performed using the Qubit 3.0 fluorometer (Life Technologies, Carlsbad, CA, USA) and the TapeStation 4200 system (Agilent Technologies, Santa Clara, CA, USA).

Affymetrix human Clariom D arrays (Affymetrix, Santa Clara, CA, US) were used for the gene expression profiling in BD2decide cohort, targeting 540,000 transcripts. Probe synthesis was performed from total RNA using the GeneChip WT Pico Reagent Kit and WT Labeling Kit (Affymetrix) [68]. A total of 6 cycles pre-in vitro transcription amplification was performed according to the manufacturer’s protocol. Biotinylated and fragmented single-stranded cDNAs were hybridized to the arrays that were washed and stained using an FS-450 fluidics station (Affymetrix, fluidics protocol FS450_0001). Signal intensities were detected by a 30007G gene array scanner. The scanned images were processed using the Affymetrix GeneChip Command Console software. CEL data files were processed with Transcriptome Analysis Console Software v4.0.1 and further filtering procedures were performed in R software (version 4.0.3).

A survey of literature was performed to retrieve gene expression signatures generated for the prognosis in HNSCC/OCSCC tumor cohorts and subsequently reproduce the signatures’ score on the OCSCC prospective subset 2 (n = 90 patients). First, we created a database of signatures and the following inclusion criteria were applied: (i) paper analyzing whole gene expression data; (ii) paper reporting prognostic gene expression signatures; (iii) complete description of the bioinformatics methods; (iv) availability of the gene list, weights, and algorithm to compute signature’s score. Inclusion criteria were applied for a survey in literature in public database such as “Pubmed” ( and “EMBASE” ( imposing a selection using keywords. Filters includes MeSH terms such as “HNSCC”, “signatures”, and “gene expression”. Consequently, each paper was carefully evaluated to exclude those based on proteins or immunohistochemistry. The genes were re-annotated based on EntrezID [69] and the bioinformatics methods were retrieved from the original manuscripts to reproduce the signature’s algorithms allowing applying them on external datasets along with thresholds for patient stratification.

Statistical analysis

The characteristics of retrospective and prospective cohorts were compared using Mann-Whitney U tests for continuous variables and chi-squared test for categorical variables. Results of the Cox proportional hazard regression model are presented as Harrel’s C-index (median and IQR), HR (median and 95% confidence intervals (CIs)) and p-value of the log-rank test. To compare the distributions of C-index the Mann-Whitney U test was applied in case of comparison between two groups and the Kruskal-Wallis test with Tukey-Kramer p-value correction was applied in case of multiple comparison. All the statistical analyses were performed in Matlab 2022a.


Patient characteristics

A dataset of 231 patients affected by stage III/IVa-b OCSCC, included in the BD2Decide database [47] for whom an MRI evaluation was available was considered for this study. The patient dataset comprised 123 retrospective patients and 108 prospective ones. Table 2 shows the clinical data of the selected patients included in the study.

The two sets, beside the expected difference in median follow-up, significantly differed in the higher percent of stage IVa-b and smoking status in prospective patients.

Table 2 Clinical data of the patients used for the study

Features selection and survival model training and testing

Features selection and survival model training was performed on the retrospective dataset composed of 123 patients. From the initial 1072 radiomic features, 222 stable features (67 T1w and 155 T2w) were considered. Of these, the final 5 selected features are listed in Table 3, together with the respective means and standard deviations used for the Z-score normalization and the corresponding regression coefficients of the trained Cox model. All the 5 selected features were extracted from the waveletLLL transform of the T2w images, with 4 textural-related features and 1 first order statistics feature. The trained Cox model demonstrated a significant prognostic value of the radiomic signature for OS on the retrospective cohort: C-index 0.68 (IQR 0.66–0.70), HR 2.64 (95% CI 1.62–4.31), log-rank p < 0.001. The Cox model applied on the 108 prospective patients maintained for the OS a significantly prognostic value, even if less relevant, due to the shorter follow-up: C-index 0.62 (IQR 0.58–0.64), HR 2.12 (95% CI 1.04–4.47), log-rank p = 0.05. The Kaplan-Meier curves for the retrospective and prospective datasets (Fig. 4A and B respectively) in both cases significantly separated the high- from the low-risk groups. In both the retrospective and prospective cohorts, the radiomic signature outperformed the cTNM staging, with higher stratification capability between high- and low-risk groups and higher C-index (Supplementary figures S1 and S2).

Table 3 Selected features and their mean (standard deviation) and regression coefficients
Fig. 4
figure 4

Radiomic model. (A) Kaplan-Meier curves for the retrospective dataset (n = 123 patients, used for model training). (B) Kaplan-Meier curves for the prospective dataset (n = 108 patients, used for model testing). For comparison, a follow-up time of 36 months was displayed. Shadows represent 95% confidence interval

Model comparison and analysis

The radiomic prognostic model outperformed the cTNM and pTNM staging, both in terms of C-index and high/low risk patient stratification. Table 4 lists the results of the radiomic signature and the cTNM and pTNM staging in the entire prospective set of 108 patients and the prospective subset 1 (75 patients for whom the pTNM was available).

As shown in Fig. 5, on the prospective subset 1 of 75 patients, the radiomic signature successfully separated the low/high risk group patients (log-rank p = 0.03), while the differences between cTNM and pTNM classes were not statistically significant (Fig. 5A and B C). Moreover, the C-index of the radiomic prognostic model was 0.62 (IQR 0.58–0.65) compared to 0.54 (IQR 0.52–0.55) of the cTNM staging and of 0.56 (IQR 0.53–0.61) of the pTNM staging (Kruskal-Wallis p < 0.001, with Tukey-Kramer correction) (Table 4; Fig. 5D).

Fig. 5
figure 5

Comparison of the prognostic performance of the radiomic model with the clinical tumor-node-metastasis (cTNM) stage and pathological tumor-node-metastasis (pTNM) stage on the prospective subset 1 (n = 75 prospective patients). (A/B) Kaplan-Meier curves for the cTNM/pTNM, with low-risk corresponding to cTNM/pTNM stage = III and high-risk corresponding to cTNM/pTNM stage = IV. Shadows represent 95% confidence interval. (C) Kaplan-Meier curves for the radiomic signature. Shadows represent 95% confidence interval. (D) Concordance indexes (C-index) for the cTNM, pTNM and radiomic signatures. *p < 0.05 (Kruskal-Wallis with Tukey-Kramer correction)

From the literature survey, after filtering for the eligibility criteria, 7 gene expression prognostic signatures [48,49,50,51,52,53,54] (Supplementary Table 1) were identified and applied to the prospective subset 2 (n = 90 prospective patients). The prognostic performance of the radiomic model on the prospective subset 2 (patients with available gene expression) was compared with that of the 7 gene expression signatures. Figure 6 shows the C-index of the 7 genomic signatures, the cTNM staging and the radiomic signature. In particular, the radiomic model performed significantly better than 6 over 7 genomic signatures and the cTNM staging (Kruskal-Wallis p < 0.001, with Tukey-Kramer correction). For G-3 genomic signature the C-index (0.64, IQR 0.60–0.68) did not present statistical difference from the radiomic model (0.66, IQR 0.60–0.69).

Fig. 6
figure 6

Comparison of the prognostic performance of the radiomic model with 7 prognostic genomic signatures and the clinical tumor-node-metastasis stage (cTNM) in terms of concordance index (C-index), on the prospective subset 2 (n = 90 prospective patients). Grey boxplots: C-index of the 7 genomic signatures; blue boxplot: C-index of the cTNM staging; red boxplot: C-index of the radiomic signature (Rad). The C-index of the radiomic model was significantly different (Kruskal-Wallis p < 0.001, with Tukey-Kramer correction) from that of the cTNM staging and of 6 over 7 genomic signatures (G-1, G-2, G-4, G-5, G-6 and G-7). No statistical difference was found between the C-index of the radiomic signature and G-3

Finally, as shown in Fig. 7, the radiomic signature was significantly different between pTNM stages III and IVa/b (Kruskal-Wallis p < 0.001, with Tukey-Kramer correction). No significant differences in the radiomic signature distribution were observed between the pTNM stage IVa and IVb. By combining the radiomic signature and the pTNM staging in a multivariate Cox model applied to the dataset of 177 patients, we found that both are independent prognostic factors of the OS (p < 0.001, HR 1.85 and 1.73, respectively).

Fig. 7
figure 7

Radiomic signature distributions in patients classified according to the pathologic tumor-node-metastasis stage (pTNM = III – IVa – IVb). *p < 0.05 (Kruskal-Wallis, with Tukey-Kramer correction). n = 37 patients presented pTNM = III; n = 76 patients presented pTNM = IVa; n = 64 patients presented pTNM = IVb.

Table 4 Performances of radiomic, clinical and pathologic tumor-node-metastasis signatures


An MRI radiomic prognostic signature, based on features extracted from T1w and T2w images of primary tumor and composed of 5 features coming from T2w images, was developed using locally advanced OCSCC from a multi-centric retrospective cohort and tested on a multi-centric prospective cohort. The radiomic signature demonstrated a significant prognostic power for OS and successfully stratified patients in low/high risk groups by Kaplan-Meier curves in both cohorts, even if the follow-up of prospective cohort was significantly shorter than that of retrospective cohort.

To date, only three studies developed MRI-radiomic prognostic signatures (based on T1w and/or T2w images) in OCSCC [32, 34, 42]. Mes et al. [32] developed a T1w-based radiomic signature on a mono-centric retrospective cohort (102 patients) and tested it on an external mono-centric retrospective cohort (76 patients). Noteworthy, in both cohorts 30–40% of stage I and II were included and the median follow-up was similar to our retrospective cohort. Despite these differences, the prognostic power and the low/high risk patient stratification performance obtained from their MRI-radiomic signature was comparable to the ones achieved herein. Wang et al. [34] developed a T2w-based radiomic signature to predict lymph node metastasis on a mono-centric retrospective cohort of 236 patients with tongue cancer (training set of 157 patients and test set of 79 patients). The combined clinical-radiomic signature was found to be an independent prognostic factor for poor OS in a multivariate Cox regression analysis (HR of 17.46). Moreover, Mossinelli et al. [42] developed several MRI-based radiomic signature for the prognosis of OS in 79 tongue cancer patients, by considering different MRI sequences. Clinical, radiomic and clinical-radiomic models were developed. When considering the radiomic models for OS, C-indexes of 0.79, 0.73, 0.72 and 0.75 were obtained with the contrast-enhanced T1, T2w, apparent diffusion coefficient map and diffusion-weighted images, respectively. From their study, the contrast-enhanced T1 sequence provided the best performance in predicting OS, while comparable results were obtained with T2w images, apparent diffusion coefficient map and diffusion-weighted images. Although the promising results, the study is limited by the low number of patients and the lack of internal and external validation. Similar to [32, 34, 42], our study demonstrated the contribution of MRI-based radiomics in OS prediction for OCSCC patients, compared to the traditionally metrics used in clinics, namely cTNM and pTNM. Besides T1w and T2w, also diffusion-weighted-based radiomics demonstrated to be useful for predicting OS, loco-regional recurrence, cause-specific mortality [42] and the histological tumor grade [70] in OCSCC patients. However, the diffusion-weighted sequence is not performed in the clinical routine, thus limiting its potential application, and was not available for the cohorts entered in the present study. Although from one side the inclusion of radiomic features from other sequences might have increased the prognostic performance of the developed radiomic model by providing complementary information, from the other side the use of standard-of-care T1w and T2w images widens its applicability, thus representing a strength of the study.

Parallel to the radiomic approach proposed herein, end-to-end radiomic-based deep learning models have been proposed (e.g., [71, 72]) and used to predict OS in HNSCC patients [73]. However, differently from the deep features, radiomic features can be linked to tissue properties, as the shape, size and texture, thus allowing for an easier interpretability of the model.

When compared to cTNM and pTNM, our radiomic model showed a higher performance in forecasting survival but the two prognostic factors (radiomics and pTNM) resulted independent predictors of outcome at multivariable analysis. Interestingly, even if evaluated in a limited number of cases (75 patients), our radiomic signature was significantly different between pathologic III and IVa/b stages suggesting that it may be useful to predict pathologic stage in patients receiving surgery. This may have an impact on treatment personalization.

The establishment of BD2Decide database provides numerous opportunities for conducting different analyses of clinical/radiomics/transcriptomics data to explore their single or combined prognostic role (see further data on [47]). As a preliminary analysis, we selected 7 gene expression prognostic signatures from the available literature and compared their performance with that of radiomic signature using the prospective dataset 2 (90 patients). Only one gene expression signature (G-3, [50]) presented, in terms of C-index, a prognostic performance comparable to that of the developed radiomic model. Also in this case, the difference in the cohorts where the signatures were developed, the anatomical subsite origin of the analyzed tumors (OCSCC versus all HNSCC) and the length of follow-up should be considered as confounding aspects. Deeper analyses and possibly a biological oriented characterization of the radiomic features will hopefully enable to better interpret these data.

The present study is not exempt from limitations. First, the small sample size and the different follow-up between the retrospective and prospective cohorts may have limited the prognostic power. Second, the prognostic performance was higher in the retrospective cohort (used for training) than in the retrospective cohort (used for validation). This can be related to the differences between the two cohorts. Cross-validation is a potential solution to address overfitting. However, model training on a retrospective cohort combined with validation on a prospective one is recommended whenever possible. Moreover, compared to cross-validation, train/validation split allows to obtained a unique prognostic model which can be applied to new external datasets. Third, only a feature selection pipeline, based on both univariate and multivariate Cox regression, was considered, based on previous studies [62]. However, as reported by Parmar et al. [31], in which 13 feature selection methods and 11 machine learning classifiers were compared in terms of prognostic performance in HNSCC patients, the feature selection method accounted for about the 14% of the total variance in the area under receiver operator characteristic curve. In future, other feature selection methods, as the least absolute shrinkage and selection operator [37, 39] or the minimum redundancy maximum relevance [31] methods will be explored as well as regularization methods for feature selection for Cox regression model [74]. Fourth, only T1w and T2w image sequences were considered because the contrast-enhanced T1 sequence was available only for a subset.

In future, the potential of radiomics in the prediction of OS in HPV+ OCSCC patients should be evaluated. Moreover, radiomic features extracted from metastatic regional lymph nodes, which were not included herein, can be considered, as done in previous works [62]. Finally, combined radiomic, clinical and genomic models can be developed to explore the additive prognostic information compared to the single models.


All taken together, our results further demonstrate that MRI, which is the most used imaging modality in HNSCC patients, contains remarkable prognostic information, and can provide a non-invasive and cost-effective prognostic factor, especially given that MRI scans (and especially T1w and T2w sequences) are performed routinely in clinical practice.

Furthermore, if confirmed in other patient cohorts, the developed radiomic signature, by potentially predicting the pathologic stage in patients receiving surgery, could be applied to support the clinical decision process.

Data Availability

All relevant data and materials have been included in the article and its supplementary data files. Further inquiries can be directed to the corresponding authors.



Concordance index


Confidence interval


Clinical tumor-node-metastasis


Grey level co-occurrence matrix


Grey level run length matrix


Head and neck squamous cell carcinomas


Hazard ratio


Intra-class correlation coefficient


Interquartile range


Magnetic resonance imaging


Oral cavity squamous cell carcinoma


Overall survival


Pathological tumor-node-metastasis


Region of interest






  1. Ferlay J, Colombet M, Soerjomataram I, Dyba T, Randi G, Bettio M, et al. Cancer incidence and mortality patterns in Europe: estimates for 40 countries and 25 major cancers in 2018. Eur J Cancer Elsevier Ltd. 2018;103:356–87.

    Article  CAS  Google Scholar 

  2. Ribeiro IP, Barroso L, Marques F, Melo JB, Carreira IM. Early detection and personalized treatment in oral cancer: the impact of omics approaches. Mol Cytogenet Molecular Cytogenetics. 2016;9:1–7.

    Google Scholar 

  3. Méndez E, Houck JR, Doody DR, Fan W, Lohavanichbutr P, Rue TC, et al. A genetic expression profile associated with oral cancer identifies a group of patients at high risk of poor survival. Clin Cancer Res. 2009;15:1353–61.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Lohavanichbutr P, Méndez E, Holsinger FC, Tessa C. A 13-gene signature prognostic of HPV-negative OSCC: discovery and external validation. Clin Cancer Res. 2009;377:364–77.

    Google Scholar 

  5. De Cecco L, Bossi P, Locati L, Canevari S, Licitra L. Comprehensive gene expression meta-analysis of head and neck squamous cell carcinoma microarray data defines a robust survival predictor. Ann Oncol. 2014;25:1628–35.

    Article  PubMed  Google Scholar 

  6. Aerts HJWL, Velazquez ER, Leijenaar RTH, Parmar C, Grossmann P, Cavalho S, et al. Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach. Nat Commun. 2014;5:4006.

    Article  CAS  PubMed  Google Scholar 

  7. Leijenaar RTH, Carvalho S, Hoebers FJP, Aerts HJWL, Van Elmpt WJC, Huang SH, et al. External validation of a prognostic CT-based radiomic signature in oropharyngeal squamous cell carcinoma. Acta Oncol (Madr). 2015;54:1423–9.

    Article  CAS  Google Scholar 

  8. Zhang B, Tian J, Dong D, Gu D, Dong Y, Zhang L, et al. Radiomics features of multiparametric MRI as novel prognostic factors in advanced nasopharyngeal carcinoma. Clin Cancer Res. 2017;23:4259–69.

    Article  PubMed  Google Scholar 

  9. Huang EP, O’Connor JPB, McShane LM, Giger ML, Lambin P, Kinahan PE, et al. Criteria for the translation of radiomics into clinically useful tests. Nat Rev Clin Oncol England. 2023;20:69–82.

    Article  Google Scholar 

  10. Serafini MS, Lopez-Perez L, Fico G, Licitra L, De Cecco L, Resteghini C. Transcriptomics and epigenomics in head and neck cancer: available repositories and molecular signatures. Cancers head neck. England. 2020;5:2.

    Google Scholar 

  11. Lawrence MS, Sougnez C, Lichtenstein L, Cibulskis K, Lander E, Gabriel SB, et al. Comprehensive genomic characterization of head and neck squamous cell carcinomas. Nature [Internet]. 2015;517:576–82.

  12. Tonella L, Giannoccaro M, Alfieri S, Canevari S, De Cecco L. Gene expression signatures for Head and Neck Cancer patient stratification: are results ready for clinical application? Curr treat options Oncol. United States. 2017;18:32.

    Google Scholar 

  13. Chai AWY, Lim KP, Cheong SC. Translational genomics and recent advances in oral squamous cell carcinoma. Semin Cancer Biol England. 2020;61:71–83.

    Article  CAS  Google Scholar 

  14. Bruixola G, Remacha E, Jiménez-Pastor A, Dualde D, Viala A, Montón JV, et al. Radiomics and radiogenomics in head and neck squamous cell carcinoma: potential contribution to patient management and challenges. Cancer Treat Rev Netherlands. 2021;99:102263.

    Article  Google Scholar 

  15. Tanadini-Lang S, Balermpas P, Guckenberger M, Pavic M, Riesterer O, Vuong D, et al. Radiomic biomarkers for head and neck squamous cell carcinoma. Strahlentherapie und Onkol [Internet]. 2020;196:868–78.

  16. Peng Z, Wang Y, Wang Y, Jiang S, Fan R, Zhang H, et al. Application of radiomics and machine learning in head and neck cancers. Int J Biol Sci Australia. 2021;17:475–86.

    Article  Google Scholar 

  17. Gonçalves M, Gsaxner C, Ferreira A, Li J, Puladi B, Kleesiek J, et al. Radiomics in Head and Neck Cancer Outcome predictions. Basel, Switzerland: Diagnostics; Switzerland; 2022. p. 12.

  18. Zhai T, Dijk LV, Van, Huang B, Lin Z, Ribeiro CO, Brouwer CL, et al. Improving the prediction of overall survival for head and neck cancer patients using image biomarkers in combination with clinical parameters. Radiother Oncol The Authors. 2017;124:256–62.

    Article  Google Scholar 

  19. Liu Z, Cao Y, Diao W, Cheng Y, Jia Z, Peng X. Radiomics-based prediction of survival in patients with head and neck squamous cell carcinoma based on pre- and post-treatment (18)F-PET/CT. Aging (Albany NY). United States. 2020;12:14593–619.

    CAS  Google Scholar 

  20. Liao KY-K, Chiu C-C, Chiang W-C, Chiou Y-R, Zhang G, Yang S-N, et al. Radiomics features analysis of PET images in oropharyngeal and hypopharyngeal cancer. Med (Baltimore) United States. 2019;98:e15446.

    Article  Google Scholar 

  21. Lv W, Ashrafinia S, Ma J, Lu L, Rahmim A. Multi-level Multi-Modality Fusion Radiomics: application to PET and CT imaging for prognostication of Head and Neck Cancer. IEEE J Biomed Heal informatics United States. 2020;24:2268–77.

    Article  Google Scholar 

  22. Xie C, Du R, Ho JW, Pang HH, Chiu KW, Lee EY, et al. Effect of machine learning re-sampling techniques for imbalanced datasets in (18)F-FDG PET-based radiomics model on prognostication performance in cohorts of head and neck cancer patients. Eur J Nucl Med Mol Imaging Germany. 2020;47:2826–35.

    Article  Google Scholar 

  23. Haider SP, Zeevi T, Baumeister P, Reichel C, Sharaf K, Forghani R et al. Potential Added Value of PET/CT Radiomics for Survival Prognostication beyond AJCC 8th Edition Staging in Oropharyngeal Squamous Cell Carcinoma. Cancers (Basel). Switzerland; 2020;12.

  24. Leger S, Zwanenburg A, Pilz K, Lohaus F, Linge A, Zöphel K, et al. A comparative study of machine learning methods for time-To-event survival data for radiomics risk modelling. Sci Rep. 2017;7:1–11.

    Article  CAS  Google Scholar 

  25. Bogowicz M, Jochems A, Deist TM, Tanadini-Lang S, Huang SH, Chan B, et al. Privacy-preserving distributed learning of radiomics to predict overall survival and HPV status in head and neck cancer. Sci Rep England. 2020;10:4542.

    Article  Google Scholar 

  26. Ou D, Blanchard P, Rosellini S, Levy A, Nguyen F, Leijenaar RTH, et al. Predictive and prognostic value of CT based radiomics signature in locally advanced head and neck cancers patients treated with concurrent chemoradiotherapy or bioradiotherapy and its added value to human papillomavirus status. Oral Oncol England. 2017;71:150–5.

    Article  Google Scholar 

  27. Cozzi L, Franzese C, Fogliata A, Franceschini D, Navarria P, Tomatis S, et al. Predicting survival and local control after radiochemotherapy in locally advanced head and neck cancer by means of computed tomography based radiomics. Strahlentherapie und Onkol Organ der Dtsch Rontgengesellschaft. [et al]. Germany. 2019;195:805–18.

    Google Scholar 

  28. Chen L, Wang H, Zeng H, Zhang Y, Ma X. Evaluation of CT-based radiomics signature and nomogram as prognostic markers in patients with laryngeal squamous cell carcinoma. Cancer imaging Off Publ Int Cancer Imaging Soc England. 2020;20:28.

    CAS  Google Scholar 

  29. Ger RB, Zhou S, Elgohari B, Elhalawani H, Mackin DM, Meier JG, et al. Radiomics features of the primary tumor fail to improve prediction of overall survival in large cohorts of CT- and PET-imaged head and neck cancer patients. PLoS One United States. 2019;14:e0222509.

    Article  CAS  Google Scholar 

  30. Foley KG, Shi Z, Whybra P, Kalendralis P, Larue R, Berbee M, et al. External validation of a prognostic model incorporating quantitative PET image features in oesophageal cancer. Radiother Oncol J Eur Soc Ther Radiol Oncol Ireland. 2019;133:205–12.

    Article  Google Scholar 

  31. Parmar C, Grossmann P, Rietveld D, Rietbergen MM, Lambin P, Aerts HJWL. Radiomic Machine-Learning Classifiers for prognostic biomarkers of Head and Neck Cancer. Front Oncol Switzerland. 2015;5:272.

    Google Scholar 

  32. Mes SW, van Velden FHP, Peltenburg B, Peeters CFW, te Beest DE, van de Wiel MA, et al. Outcome prediction of head and neck squamous cell carcinoma by MRI radiomic signatures. Eur Radiol European Radiology. 2020;30:6311–21.

    Article  PubMed  Google Scholar 

  33. Bos P, van den Brekel MWM, Gouw ZAR, Al-Mamgani A, Taghavi M, Waktola S, et al. Improved outcome prediction of oropharyngeal cancer by combining clinical and MRI features in machine learning models. Eur J Radiol Ireland. 2021;139:109701.

    Article  Google Scholar 

  34. Wang F, Tan R, Feng K, Hu J, Zhuang Z, Wang C, et al. Magnetic resonance imaging-based Radiomics features Associated with depth of Invasion Predicted Lymph Node Metastasis and Prognosis in Tongue Cancer. J Magn Reson Imaging United States. 2022;56:196–209.

    Article  Google Scholar 

  35. Jethanandani A, Lin TA, Volpe S, Elhalawani H, Mohamed ASR, Yang P, et al. Exploring applications of Radiomics in magnetic resonance imaging of Head and Neck Cancer: a systematic review. Front Oncol. 2018;8:131.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Yuan Y, Ren J, Shi Y, Tao X. MRI-based radiomic signature as predictive marker for patients with head and neck squamous cell carcinoma. Eur J Radiol Ireland. 2019;117:193–8.

    Article  Google Scholar 

  37. Alfieri S, Romanò R, Bologna M, Calareso G, Corino V, Mirabile A, et al. Prognostic role of pre-treatment magnetic resonance imaging (MRI)-based radiomic analysis in effectively cured head and neck squamous cell carcinoma (HNSCC) patients. Acta Oncol England. 2021;60:1192–200.

    Article  CAS  Google Scholar 

  38. Siow TY, Yeh C-H, Lin G, Lin C-Y, Wang H-M, Liao C-T et al. MRI Radiomics for Predicting Survival in Patients with Locally Advanced Hypopharyngeal Cancer Treated with Concurrent Chemoradiotherapy. Cancers (Basel). Switzerland; 2022;14.

  39. Chen J, Lu S, Mao Y, Tan L, Li G, Gao Y, et al. An MRI-based radiomics-clinical nomogram for the overall survival prediction in patients with hypopharyngeal squamous cell carcinoma: a multi-cohort study. Eur Radiol Germany. 2022;32:1548–57.

    Article  Google Scholar 

  40. Boot PA, Mes SW, de Bloeme CM, Martens RM, Leemans CR, Boellaard R, et al. Magnetic resonance imaging based radiomics prediction of human papillomavirus infection status and overall survival in oropharyngeal squamous cell carcinoma. Oral Oncol England. 2023;137:106307.

    Article  CAS  Google Scholar 

  41. Bologna M, Corino V, Cavalieri S, Calareso G, Gazzani SE, Poli T, et al. Prognostic radiomic signature for head and neck cancer: development and validation on a multi-centric MRI dataset. Radiother Oncol J Eur Soc Ther Radiol Oncol Ireland. 2023;183:109638.

    Article  Google Scholar 

  42. Mossinelli C, Tagliabue M, Ruju F, Cammarata G, Volpe S, Raimondi S, et al. The role of radiomics in tongue cancer: a new tool for prognosis prediction. Head Neck. United States; 2023.

  43. Ren J, Tian J, Yuan Y, Dong D, Li X, Shi Y, et al. Magnetic resonance imaging based radiomics signature for the preoperative discrimination of stage I-II and III-IV head and neck squamous cell carcinoma. Eur J Radiol Ireland. 2018;106:1–6.

    Article  Google Scholar 

  44. Yu B, Huang C, Xu J, Liu S, Guan Y, Li T, et al. Prediction of the degree of pathological differentiation in tongue squamous cell carcinoma based on radiomics analysis of magnetic resonance images. BMC Oral Health England. 2021;21:585.

    Article  CAS  Google Scholar 

  45. Frood R, Palkhi E, Barnfield M, Prestwich R, Vaidyanathan S, Scarsbrook A. Can MR textural analysis improve the prediction of extracapsular nodal spread in patients with oral cavity cancer? Eur Radiol. Germany. 2018;28:5010–8.

    Google Scholar 

  46. Machiels J-P, René Leemans C, Golusinski W, Grau C, Licitra L, Gregoire V. Squamous cell carcinoma of the oral cavity, larynx, oropharynx and hypopharynx: EHNS-ESMO-ESTRO Clinical Practice Guidelines for diagnosis, treatment and follow-up. Ann. Oncol. Off. J. Eur. Soc. Med. Oncol. England; 2020. p. 1462–75.

  47. Cavalieri S, De Cecco L, Brakenhoff RH, Serafini MS, Canevari S, Rossi S et al. Development of a multiomics database for personalized prognostic forecasting in head and neck cancer: The Big Data to Decide EU Project. Head Neck. 2020;1–12.

  48. Schomberg J, Ziogas A, Anton-Culver H, Norden-Krichmar T. Identification of a gene expression signature predicting survival in oral cavity squamous cell carcinoma using Monte Carlo cross validation. Oral Oncol England. 2018;78:72–9.

    Article  CAS  Google Scholar 

  49. Guo W, Chen X, Zhu L, Wang Q. A six-mRNA signature model for the prognosis of head and neck squamous cell carcinoma. Oncotarget United States. 2017;8:94528–38.

    Article  Google Scholar 

  50. Wintergerst L, Selmansberger M, Maihoefer C, Schüttrumpf L, Walch A, Wilke C, et al. A prognostic mRNA expression signature of four 16q24.3 genes in radio(chemo)therapy-treated head and neck squamous cell carcinoma (HNSCC). Mol Oncol United States. 2018;12:2085–101.

    Article  CAS  Google Scholar 

  51. Lohavanichbutr P, Méndez E, Holsinger FC, Rue TC, Zhang Y, Houck J, et al. A 13-gene signature prognostic of HPV-negative OSCC: discovery and external validation. Clin cancer Res an Off J Am Assoc Cancer Res United States. 2013;19:1197–203.

    Article  CAS  Google Scholar 

  52. Liu J, Han M, Yue Z, Dong C, Wen P, Zhao G, et al. Prediction of Radiosensitivity in Head and Neck squamous cell Carcinoma based on multiple Omics Data. Front Genet Switzerland. 2020;11:960.

    Article  CAS  Google Scholar 

  53. Wu F, Du Y, Hou X, Cheng W. A prognostic model for oral squamous cell carcinoma using 7 genes related to tumor mutational burden. BMC Oral Health England. 2022;22:152.

    Article  Google Scholar 

  54. Chen Y, Feng Y, Yan F, Zhao Y, Zhao H, Guo Y. A Novel Immune-Related gene signature to identify the Tumor Microenvironment and Prognose Disease among patients with oral squamous cell carcinoma patients using ssGSEA: a Bioinformatics and Biological Validation Study. Front Immunol Switzerland. 2022;13:922195.

    Article  CAS  Google Scholar 

  55. Jung F, Steger S, Knapp O, Noll M, Wesarg S. COSMO-coupled shape model for radiation therapy planning of head and neck cancer. Work Clin Image-Based Proced Transl Res Med Imaging CLIP 2014. 2014. p. 25–32.

  56. Bologna M, Corino V, Mainardi L. Technical note: virtual phantom analyses for preprocessing evaluation and detection of a robust feature set for MRI-radiomics of the brain. Med Phys. 2019;46:5116–23.

    Article  PubMed  Google Scholar 

  57. Tustison NJ, Cook PA, Gee JC. N4ITK: improved N3 bias correction. IEEE Trans Med Imaging. 2010;29:1310–20.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Leijenaar RT, Bogowicz M, Jochems A, Hoebers FJ, Wesseling FW, Huang SH, et al. Development and validation of a radiomic signature to predict HPV (p16) status from standard CT imaging: a multicenter study. Br J Radiol. 2018;91:20170498.

    Article  PubMed  PubMed Central  Google Scholar 

  59. van Griethuysen JJM, Fedorov A, Parmar C, Hosny A, Aucoin N, Narayan V, et al. Computational Radiomics System to Decode the Radiographic phenotype. Cancer Res. 2017;77:e104–7.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Pyradiomics. features description. 2018.

  61. Ioffe S, Szegedy C. Batch normalization: Accelerating deep network training by reducing internal covariate shift. 32nd Int Conf Mach Learn ICML 2015. 2015;1:448–56.

  62. Bologna M, Corino V, Calareso G, Tenconi C, Alfieri S, Iacovelli NA, et al. Baseline mri-radiomics can predict overall survival in non-endemic ebv-related nasopharyngeal carcinoma patients. Cancers (Basel). 2020;12:1–20.

    Article  Google Scholar 

  63. Bologna M, Corino VDA, Montin E, Messina A, Calareso G, Greco G, et al. Assessment of Stability and discrimination capacity of Radiomic features on apparent diffusion coefficient images. J Digit Imaging. 2018;31:879–94.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Harrell FE, Kerry LL, Mark DB. Tutorial in biostatistics multivariable prognostic models: issues in developing models, evaluating assumptions and adequacy, and measuring and reducing errors. Stat Med. 1996;15:361–87.

    Article  PubMed  Google Scholar 

  65. James G, Witten D, Hastie T, Tibshirani R. An introduction to statistical learning - with applications in R. Springer; 2013.

  66. Peto R, Peto J. Asymptotically efficient rank invariant test procedures. J R Stat Soc. 1972;135:185–207.

    Google Scholar 

  67. Kaplan EL, Meier P. Nonparametric estimation from incomplete samples. J Am Stat Assoc. 1958;73:457–81.

    Article  Google Scholar 

  68. Cavalieri S, Serafini MS, Carenzo A, Canevari S, Brakenhoff RH, Leemans CR, et al. Clinical validity of a prognostic gene expression cluster-based model in human papillomavirus-positive Oropharyngeal Carcinoma. Volume 5. JCO Precis Oncol. United States; 2021.

  69. Maglott D, Ostell J, Pruitt KD, Tatusova T. Entrez Gene: gene-centered information at NCBI. Nucleic Acids Res England. 2011;39:D52–7.

    Article  CAS  Google Scholar 

  70. Ren J, Qi M, Yuan Y, Tao X. Radiomics of apparent diffusion coefficient maps to predict histologic grade in squamous cell carcinoma of the oral tongue and floor of mouth: a preliminary study. Acta Radiol England. 2021;62:453–61.

    Article  Google Scholar 

  71. Thammasorn P, Chaovalitwongse WA, Hippe DS, Wootton LS, Ford EC, Spraker MB, et al. Nearest neighbor-based strategy to optimize Multi-View Triplet Network for classification of Small-Sample Medical Imaging Data. IEEE Trans Neural Networks Learn Syst. 2023;34:586–600.

    Article  Google Scholar 

  72. Thammasorn P, Schaub SK, Hippe DS, Spraker MB, Peeken JC, Wootton LS, et al. Regularizing the Deepsurv Network using projection loss for Medical Risk Assessment. IEEE Access. 2022;10:8005–20.

    Article  Google Scholar 

  73. Tortora M, Gemini L, Scaravilli A, Ugga L, Ponsiglione A, Stanzione A et al. Radiomics Applications in Head and Neck Tumor Imaging: A Narrative Review. Cancers (Basel). Switzerland; 2023;15.

  74. Simon N, Friedman J, Hastie T, Tibshirani R. Regularization Paths for Cox’s proportional hazards Model via Coordinate Descent. J Stat Softw United States. 2011;39:1–13.

    Google Scholar 

Download references


The authors acknowledge Silvana Canevari for the helpful discussion.


Not applicable.

Author information

Authors and Affiliations



AC: methodology, formal analysis, results interpretation, manuscript writing and revision; LDC: clinical data collection, gene expression analysis, result interpretation, manuscript revision; SC: clinical data collection, results interpretation, manuscript revision; DL, FP: clinical data collection, gene expression analysis; GC, DM, PG, CRL, RHB, MR: clinical data collection, image segmentation; TP: clinical data collection, manuscript revision; LL: clinical data collection, results interpretation, manuscript revision; VC and LM: conceptualization, supervision, results interpretation, manuscript revision. All authors have read and agreed to the published version of the manuscript.

Corresponding author

Correspondence to Anna Corti.

Ethics declarations

Ethics approval and consent to participate

The protocols were approved by the Ethical Committees of the participating centers, and data acquisition followed the General Data Protection Regulation of the EU. All patients signed the informed consent.

Consent for publication

All authors have read and approved the manuscript and agree with submission to Biomarker Research.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Corti, A., De Cecco, L., Cavalieri, S. et al. MRI-based radiomic prognostic signature for locally advanced oral cavity squamous cell carcinoma: development, testing and comparison with genomic prognostic signatures. Biomark Res 11, 69 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Radiomics
  • Radiogenomics
  • Magnetic resonance imaging
  • Head and neck cancer
  • Survival models
  • Overall survival
  • Oral cavity squamous cell carcinoma