Next Article in Journal
Effects of In-Hospital Exercise on Frailty in Patients with Hepatocellular Carcinoma
Next Article in Special Issue
Dedicated Axillary MRI-Based Radiomics Analysis for the Prediction of Axillary Lymph Node Metastasis in Breast Cancer
Previous Article in Journal
Stereotactic Radiotherapy for Hepatocellular Carcinoma, Radiosensitization Strategies and Radiation-Immunotherapy Combination
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

MRI-Based Radiomics Input for Prediction of 2-Year Disease Recurrence in Anal Squamous Cell Carcinoma

by
Nicolas Giraud
1,2,*,
Olivier Saut
2,
Thomas Aparicio
3,
Philippe Ronchin
4,
Louis-Arnaud Bazire
5,
Emilie Barbier
6,
Claire Lemanski
7,
Xavier Mirabel
8,
Pierre-Luc Etienne
9,
Astrid Lièvre
10,
Wulfran Cacheux
11,
Ariane Darut-Jouve
12,
Christelle De la Fouchardière
13,
Arnaud Hocquelet
14,
Hervé Trillaud
14,
Thomas Charleux
1,
Gilles Breysacher
15,
Delphine Argo-Leignel
16,
Alexandre Tessier
17,
Nicolas Magné
18,
Meher Ben Abdelghani
19,
Côme Lepage
20 and
Véronique Vendrely
1
add Show full author list remove Hide full author list
1
Radiation Oncology Department, Hôpital Haut-Lévêque, CHU Bordeaux, 33600 Pessac, France
2
Modelisation in Oncology (MOnc) Team, INRIA Bordeaux-Sud-Ouest, CNRS UMR 5251 and Université de Bordeaux, 33400 Talence, France
3
CHU St Louis, 75010 Paris, France
4
Centre Azuréen de Cancérologie, 06250 Mougins, France
5
Institut Curie, 75005 Paris, France
6
Fédération Francophone de Cancérologie Digestive, 21000 Dijon, France
7
Institut du Cancer de Montpellier, 34090 Montpellier, France
8
Centre Oscar Lambret, 59000 Lille, France
9
Centre Armoricain de Radiothérapie, 22190 Plérin, France
10
CHU Rennes, 35000 Rennes, France
11
Hôpital Privé Pays de Savoie, 74100 Annemasse, France
12
Institut de Cancérologie de Bourgogne, 21000 Dijon, France
13
Centre Léon Bérard, 69008 Lyon, France
14
Radiology Department, CHU Bordeaux, 33000 Bordeaux, France
15
Hôpitaux Civils de Colmar, 68000 Colmar, France
16
GH Bretagne Sud, 56100 Lorient, France
17
CH Annecy, 74370 Annecy, France
18
Institut de Cancérologie Lucien Neuwirth, 42270 St-Priest-en-Jarez, France
19
Institut de Cancérologie Strasbourg Europe, 67200 Strasbourg, France
20
CHU Dijon, 21000 Dijon, France
*
Author to whom correspondence should be addressed.
Cancers 2021, 13(2), 193; https://doi.org/10.3390/cancers13020193
Submission received: 15 November 2020 / Revised: 20 December 2020 / Accepted: 1 January 2021 / Published: 7 January 2021
(This article belongs to the Special Issue PET and MRI Radiomics in Cancer Predictive Modeling)

Abstract

:

Simple Summary

Exclusive chemo-radiotherapy (CRT) is the standard treatment for non-metastatic anal squamous cell carcinomas. Identifying novel prognostic factors could help to improve CRT outcomes, notably for locally advanced diseases where relapses still occur in around 35% of patients. In this study, we aim to assess the potential value of a pre-therapeutic MRI radiomic analysis added to standard clinical variables in order to build a logistic regression model predicting 2-year recurrence after CRT. In a population of 82 patients randomly divided in training (n = 54) and testing (n = 28) sets, after selection of optimal variables, a model using two radiomic (FirstOrder_Entropy and GLCM_JointEnergy) and two clinical (tumor size and CRT length) features was able to predict the 2-year recurrence with good performances in the testing set. Radiomic biomarkers provided valuable additional and independent information added to clinical data, and could help contribute to identify high risk patients amenable to treatment intensification with view of personalized medicine.

Abstract

Purpose: Chemo-radiotherapy (CRT) is the standard treatment for non-metastatic anal squamous cell carcinomas (ASCC). Despite excellent results for T1-2 stages, relapses still occur in around 35% of locally advanced tumors. Recent strategies focus on treatment intensification, but could benefit from a better patient selection. Our goal was to assess the prognostic value of pre-therapeutic MRI radiomics on 2-year disease control (DC). Methods: We retrospectively selected patients with non-metastatic ASCC treated at the CHU Bordeaux and in the French FFCD0904 multicentric trial. Radiomic features were extracted from T2-weighted pre-therapeutic MRI delineated sequences. After random division between training and testing sets on a 2:1 ratio, univariate and multivariate analysis were performed on the training cohort to select optimal features. The correlation with 2-year DC was assessed using logistic regression models, with AUC and accuracy as performance gauges, and the prediction of disease-free survival using Cox regression and Kaplan-Meier analysis. Results: A total of 82 patients were randomized in the training (n = 54) and testing sets (n = 28). At 2 years, 24 patients (29%) presented relapse. In the training set, two clinical (tumor size and CRT length) and two radiomic features (FirstOrder_Entropy and GLCM_JointEnergy) were associated with disease control in univariate analysis and included in the model. The clinical model was outperformed by the mixed (clinical and radiomic) model in both the training (AUC 0.758 versus 0.825, accuracy of 75.9% versus 87%) and testing (AUC 0.714 versus 0.898, accuracy of 78.6% versus 85.7%) sets, which led to distinctive high and low risk of disease relapse groups (HR 8.60, p = 0.005). Conclusion: A mixed model with two clinical and two radiomic features was predictive of 2-year disease control after CRT and could contribute to identify high risk patients amenable to treatment intensification with view of personalized medicine.

1. Introduction

Anal squamous cell carcinoma (ASCC) is a relatively rare yet increasing disease, reaching around 48,000 new cases worldwide in 2018 [1]. With a low rate of distant metastases, ASCC is usually amenable to loco-regional treatment. Definitive chemo-radiotherapy (CRT) is the current standard of care for non-metastatic disease based on the results of several phase III trials performed in the 1990s [2,3,4]. Since then, the improvement in radiation techniques as well as a better chemotherapy (CT) management has allowed incremental and continuous progress in disease control rates as well as reducing early and late side effects [5,6]. Nevertheless, whilst excellent for small T1-2 tumors, CRT results remain more modest for locally advanced ASCC with up to 35% local or distant recurrences, more than 80% of them occurring during the first two years [7,8,9]. Several trials with intensification strategies in terms of radiation dose or chemotherapy association have been instigated but struggled to showcase clinically significant improvements compared to the standard association, perhaps due to non-optimal patient selections [10]. Recent efforts aim to escalate the treatment by adding targeted therapies, epidermal growth factor receptor (EGFR) inhibitors, to classic CRT, showing conflicting results. Among these trials, the French phase 1–2 FFCD-0904 trial added panitumumab, an anti-EGFR drug, to standard CRT, but the expected improvement in terms of complete pathological control at 8 weeks was not met [11,12,13,14]. Novel ongoing trials involving immunotherapy seem promising, but phase III results are yet to be presented [15]. Ideally, high-risk patients could be identified at baseline in order to offer more aggressive strategies, with regards to personalized medicine.
Radiomics is a newly emerging field of study, using data-characterization algorithms and mathematical tools applied on manually or automatically segmented volumes of interest on medical images [16,17]. After several steps, each bearing its own specifications (image acquisition, data standardization, segmentation, features extraction and qualification), a great number of quantitative metrics are extracted (Figure 1). These metrics individually describe various properties of the image, ranging from first-order attributes describing voxels’ intensities or their spatial distribution to higher-order features depicting for example the relationship between two or more adjacent voxels. The subsequent step is to integrate these new quantitative items inside predictive algorithms, for instance logistic regression models, alongside more conventional data (clinical, pathological, radiological, etc). The addition of these radiomic indicators display promising results in the field of oncology for a wide range of diseases [18,19]. In ASCC specifically, few papers have been published, using mostly PET or MRI imaging, and often small cohort sizes [20].
In this study, we aimed to assess the post-CRT 2-year disease-free survival (DFS) predictive performance of a logistic regression model incorporating MRI-based pre-CRT radiomics and clinical features.

2. Results

2.1. Patients Characteristics

Between January 2010 and September 2019, 86 patients underwent exclusive CRT for non-metastatic ASCC at the Bordeaux University Hospital. Among them, six patients were treated with suboptimal radiation dose because of palliative intent and 2 received an induction CT, thus excluded. In the end, 36 patients with an available pre-therapeutic MRI and a follow-up longer than 2 years were included.
Simultaneously, patients from the FFCD-0904 trial were included in the cohort. Of the 55 patients included in the phases I and II of this trial, four did not receive an optimal dose of CRT (toxicity assessment in phase I) and were excluded. Among the remaining patients, 46 had an available pre-CRT MRI.
In the end, a total of 82 patients were selected for further examination. Patient selection is detailed in Figure 2. Clinical and histopathological did not significantly differ between the two cohorts, subsequently blended (Table 1). Patients from the joint cohort were randomly divided into training (n = 54) and testing (n = 28) sets.

2.2. Outcome

Median follow-up in the event-free population was 24 months (range, 24–101). Among the 82 patients, loco-regional or distant relapse occurred in 24 patients (29%), after a median duration of 5 months (1–14 months). Among them, 12 patients (50%) presented with in-field relapse, eight (33%) with metastatic-only relapse and four (17%) with both. Nine patients died after tumor relapse (median 11 months, range 1–40 months).

2.3. Logistic Regression

Using univariate analysis, seven variables (two clinical and five radiomic) were associated with the outcome in the training set: tumor length (p = 0.01), CRT duration (p = 0.043), FO_Entropy (p = 0.045), FO_Uniformity (p = 0.045), Glcm_JointEnergy (p = 0.047), Glcm_JointEntropy (p = 0.049) and Glrlm_GrayLevelVariance (p = 0.045). Description of these radiomic features can be found in Table 2.
Strongly intercorrelated features were investigated using a correlation matrix and Pearson tests. FO_Uniformity was found to be strongly correlated with Glrlm_GrayLevelVariance and FO_Entropy; and Glcm_JointEnergy with Glcm_JointEntropy (Figure 3). Among these, FO_Entropy for the first group and Glcm_JointEnergy for the second were most correlated with the outcome, thus retained in the model alongside the tumor length and CRT duration. In multivariate analysis, adding to these four variables the patients’ age and sex, the tumor length (p = 0.036) and Glcm_JointEntropy (p = 0.026) remained significantly associated with the 2-year DFS.
ROC curves were drawn in order to determine the optimal classification threshold for each of the four model’s variables, using the Youden method (Table 3). These thresholds allowed to dichotomize each of the four quantitative features and were applied to both the training and testing sets. A binary logistic regression model was subsequently fitted on the training set based on these four dichotomized features, and performed 2-year DFS predictions in the two sets.
In the training set, this mixed radiomic and clinical model obtained an accuracy of 87%, with an area under the curve (AUC) of the ROC curve of 0.825 (sensitivity of 58.8%, specificity of 100% to predict the recurrence). In the testing set, a prediction accuracy of 85.7% was achieved, with a ROC-AUC of 0.898 (sensitivity of 42.8%, specificity of 100% to predict the recurrence).
It allowed to distinguish two groups among the testing set according to the risk of recurrence after CRT, with a hazard ratio of 8.6 (p = 0.005), Figure 4. This model performed better in the testing set than the clinical-only (accuracy of 78.6%, ROC-AUC 0.714) and radiomics-only (accuracy of 67.9%, ROC-AUC of 0.677) models.

2.4. Inter-Reader Variability

In the relevant sample (n = 22), the two independent segmentations were relatively similar, with a mean DICE coefficient of 0.87 and an average mean Hausdorff distance of 0.68mm. In this sample, the differences in segmentation did not influence the model’s predictions.

3. Discussion

Male gender, N-positive stage and tumor length greater than 5 cm are among the few recognized ASCC clinical prognostic factors, associated with worse clinical outcomes [21,22]. Other prognostic factors have been suggested throughout the years, like an age greater than 55, increased circumferential tumor spread, skin ulceration, inguinal node development and total RT dose superior to 60 Gy for worse colostomy-free survival [23], or the baseline neutrophil to lymphocyte ratio for loco-regional recurrence [24]. Prognostic nomograms have even been elaborated, with promising results in predicting cancer specific survival and overall survival as well as risk-stratifying patients, but were only based on clinical data [25]. In that regard, radiomics features extracted from medical imaging could bring additional informative data. Adding these novel parameters to previously used clinical parameters seems to offer additive performance, paving the way to build more accurate predictive models, but have for now barely been explored in ASCC [20].
In this multicentric study, the largest to our knowledge to explore the input of MRI T2w radiomics in ASCC recurrence prediction, we showed that a simple mixed (clinical and radiomic) LR model was able to predict the 2-y DFS with good performances using optimal thresholds determined on the training set. This model incorporates two clinical and two radiomics features. Among the two studies published on this subject, Hocquelet et al. found that Skewness and Cluster Shade were correlated with recurrence or death, which we did not highlight in our study [26]. Owczarczyk et al. observed that baseline T2w Energy, a witness of voxel signal distribution (high energy being linked to homogeneous intensity distribution), was associated with 2-year DFS [27]. In our study, although FO_Energy was not statistically associated with the outcome (p = 0.38 in the training set), FO_Entropy and FO_Uniformity also depict to some extent the heterogeneity of voxels’ intensities [28]. In both trials, the homogeneity in baseline tumor voxel intensities seems to be linked to poorer outcomes.
In our LR model, thresholds were used in order to dichotomize the quantitative variables of interest using optimal values given by the Youden method on ROC curves. While this raises a possible “threshold effect”, the value found for tumor length was 49 mm, which is consistent with the literature data (50 mm usually considered as a critical prognostic threshold and the boundary between T2 and T3 in the TNM 8th edition). This value is more difficult to interpret for radiomic variables as their absolute values are harder to apprehend. Moreover, while showcasing overall good performances on prediction accuracy, specificity and ROC AUC, the results for sensitivity are poorer (42.9%), implying a possible high number of false negatives. Consequently, while reducing the risk of over-treatment (all patients predicted to present disease recurrence relapsed), the cost would be a few patients that would exhibit disease recurrence not being pre-detected (under-diagnosis). This raises a fundamental question, as one could argue that with the current CRT results in locally advanced ASCC, the priority should be given to detecting all recurring patients, even if some other patients could be over-treated. Thirdly, this model takes into account the length of CRT (removing it resulted in poorer performances), which de facto prevents from using this model in intensification trials including adjuvant chemotherapy or dose escalation. However, it highlights an important message that avoiding as much treatment interruptions as possible and reducing their duration seem to improve the patients’ outcome, pauses being increasingly questioned [29]. As a matter of fact, in a recent pooled analysis, shorter overall treatment time was associated with longer PFS [30]. Finally, a significant proportion of patients who entered in this model received an experimental drug (panitumumab) not standardly used, but this was taken into account in the model and did not seem to interfere with the outcome.
In this study, we limited our specter of imaging to T2w sequences because other modalities were unsystematically performed, but inclusion of other imaging modalities e.g., diffusion MRI sequences or PET-CT imaging could further improve the model’s predictive performances. Moreover, the benefit of adding pathological data into the model, such as the p16 status obtained on the initial biopsy or the measurement of tumor infiltrating lymphocytes (TILs), of probable prognostic value, could warrant further investigations [31,32,33].

4. Materials and Methods

4.1. Patient Selection

All patients treated by CRT at the CHU Bordeaux between 2010 and 2019 for histologically proven non-metastatic ASCC were retrospectively considered. All patients treated with palliative intent or who received induction CT were excluded, as well as those event-free with less than 2 years of follow-up.
Additionally, patients from the French FFCD-0904 phase I-II multicentric prospective study (NCT01581840) were included in the cohort. This trial included patients with T2 > 3 cm, T3-4 and/or N+ non-metastatic tumors. All participating patients signed an informed consent. The study protocol was approved by the Ethics Committee at site. The study was conducted according to the principles of the Declaration of Helsinki, the International Conference on Harmonization Guideline on Good Clinical Practice, the French laws and regulations.
All patients for which the pre-therapeutic MRI could not be retrieved were excluded.

4.2. Outcome

The primary endpoint was the 2-year DFS, relapse being defined as clinical or radiological loco-regional or distant progression, either directly after the CRT course (CRT resistance) or after an initial response. In doubtful cases, histology was mandatory.

4.3. MRI

The images were performed on different 1.5 or 3-T MRI machines, with protocols left to the choice of each expert center in order to adhere as much as possible to real-life conditions. Patients were scanned in supine position. Analyzed MRI sequences included axial T2-weighted (T2w) sequences, without fat suppression. We chose to focus on the texture analysis of T2w images because they are easy to acquire, less prone to artefacts, non-invasive and suitable for every patient. Importantly, T2w sequences offer a high signal-to-noise ratio, spatial resolution and soft tissue contrast images of the anal sphincter structures.

4.4. Clinical Features and Treatment-Related Data

The following clinical variables were collected from medical records: age, sex, HIV status, TNM stage, pre-CRT tumor length, pre-rectal septum infiltration. When applicable, date and type (local, distant or both) of relapse was recorded. All categorical clinical features were remapped to ordinal values and binarized.

4.5. Chemo-Radiotherapy

All patients underwent CRT as modality of treatment using at least 6 MV photons. The total duration of the CRT (total number of days between the first and final days of RT, treated as a continuous variable), the received radiation dose and the type of CT (including the prescription of panitumumab) were congregated.
In the Bordeaux cohort, radiation doses were 59.4 Gy in 33 fractions for all patients. CT consisted of Xeloda (1650 mg/m2 on RT treatment days) and mitomycin C (MMC, two injections of 10 mg/m2), 5-Fluorouracil (5FU, 800 mg/m2/day over 4 days) and Cisplatin (80 mg/m2) or Xeloda alone (1650 mg/m2 on RT treatment days).
In the FFCD cohort, all patients received a 65 Gy radiation dose in two consecutive sequences (45 Gy in 25 fractions followed by 20 Gy in 10 fractions), with no treatment pause. The CT protocol associated four injections of panitumumab 3 mg/kg and two injections of continuous 5FU 400 mg/m2/day over 4 days plus MMC 10 mg/m2, as defined by the phase I toxicity assessment [14].

4.6. Tumor Delineation

Anal canal tumors were manually delineated on all axial slices using Varian Eclipse® on pre-therapeutic T2w sequences and validated by two experts. A second segmentation was blindly performed by one expert for 27% (n = 22) patients, allowing to compute concordance indices based on DICE coefficients and average Hausdorff distances [34].

4.7. MRI Presampling

Prior to radiomic extraction, transformations were applied to the pre-CRT T2w sequences. A N4 bias correction algorithm was carried out based on the simple itk python library to correct MRI low frequency intensity non-uniformity. Then, an in-house python script was created to normalize intra-VOI voxel intensities.

4.8. Radiomic Features

A total of 100 radiomics quantitative features were extracted from each segmented volume using the PyRadiomics library, including first-order (FO) shape and intensity indicators as well as second-order textural features based on the gray level co-occurrence (GLCM), gray level run length (GLRLM), gray level size zone (GLSZM) and gray level dependence (GLDM) matrixes.

4.9. Statistical Analysis and Feature Selection

The cohort was first randomly split into two sets, 2/3 for training (n = 54) and 1/3 for testing (n = 28). Clinical and radiomic variables (listed in Supplementary Table S1) were assessed for their 2-year DFS predictive performance with univariate (Student if applicable or Wilcoxon tests for continuous variables and Fisher tests for categorical features, as well as receiver operating characteristic (ROC) curves) and multivariate (logistic regression) analyses.
These analyses in the training set allowed feature selection. After eliminating redundant variables using the Pearson correlation test (when two variables were significantly correlated, the variable most correlated with the outcome was chosen), optimal cut-off values for each elected variable were defined using the Youden method applied on each ROC curve. This threshold was then applied in the two sets. Quantitative performance evaluation was carried out using accuracy, sensitivity, and specificity regarding prediction of relapse and the log-rank test regarding DFS stratification. All statistical analyses were carried out using the RStudio software (v1.3.959).

5. Conclusions

After optimization and performance assessment, a four-variable model fitted with a LR algorithm achieved good performances to predict 2-y DFS. With a 100% positive predictive value, this model could help identifying patients at very high risk of recurrence, allowing for a better guidance of patients eligible to treatment intensification (adjuvant CT, or RT dose escalation for example). The ability to foretell relapse could guide therapeutic strategy and lead to intensification on a personalized basis. Adding radiomic features to clinical data improved the model’s results, individualizing some patients with a very high risk of relapse. Yet, another portion of patients who eventually relapsed were not identified by the model, illustrating their poor discrimination compared to relapse-free patients according to the existing clinical and radiomic features.
As Machine Learning models thrive with high-dimensional data, new biomarkers as well as other imaging modalities or repeated imaging obtained throughout the treatment for delta-radiomics (for instance with novel MRI linear accelerators) could further bring relevant information to the model. In addition, exploring the correlation between these features and the outcome with a longer follow-up is warranted, as well as including more patients. This is why there is a strong need of national or international collaboration, especially in neoplasms like ASCC with relatively low incidences. The next intended step is to challenge such models in a prospective setting, before considering using them in interventional trials.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6694/13/2/193/s1, Table S1. Complete list of variables tested in the model’s creation.

Author Contributions

N.G.: data collection and statistical analysis; T.C.: data processing and segmentation (inter-reader variability assessment); V.V., T.A., P.R., L.-A.B., E.B., C.L. (Claire Lemanski), X.M., P.-L.E., A.L., W.C., A.D.-J., C.D.l.F., A.H., H.T., G.B., D.A.-L., A.T., N.M., M.B.A., C.L. (Côme Lepage): data collection; O.S., V.V.: supervision and manuscript correction. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Review Board of Ile-de-France (PC/AP 11/2012 emitted on the 7th of March, 2012).

Informed Consent Statement

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

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy reasons (individual patient data, partly belonging to a FFCD academic trial unpublished to date).

Acknowledgments

A special thanks to NVIDIA® for lending 2 NVIDIA GPUs.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bray, F.; Ferlay, J.; Soerjomataram, I.; Siegel, R.L.; Torre, L.A.; Jemal, A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 2018, 68, 394–424. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Bosset, J.F.; Pavy, J.J.; Roelofsen, F.; Bartelink, H. Combined radiotherapy and chemotherapy for anal cancer EORTC Radiotherapy and Gastrointestinal Cooperative Groups. Lancet 1997, 349, 205–206. [Google Scholar] [CrossRef]
  3. UKCCCR. Epidermoid anal cancer: Results from the UKCCCR randomised trial of radiotherapy alone versus radiotherapy, 5-fluorouracil, and mitomycin. UKCCCR Anal Cancer Trial Working Party. UK Coordinating Committee on Cancer Research. Lancet 1996, 348, 1049–1054. [Google Scholar] [CrossRef]
  4. Glynne-Jones, R.; Nilsson, P.J.; Aschele, C.; Goh, V.; Peiffert, D.; Cervantes, A.; Arnold, D.; European Society for Medical Oncology (ESMO), European Society of Surgical Oncology (ESSO), European Society of Radiotherapy and Oncology (ESTRO). Anal cancer: ESMO-ESSO-ESTRO clinical practice guidelines for diagnosis, treatment and follow-up. Eur. J. Surg. Oncol. 2014, 40, 1165–1176. [Google Scholar] [CrossRef] [PubMed]
  5. Ajani, J.A.; Winter, K.A.; Gunderson, L.L.; Pedersen, J.; Benson, A.B.; Thomas, C.R.; Mayer, R.J.; Haddock, M.G.; Rich, T.A.; Willett, C. Fluorouracil, mitomycin, and radiotherapy vs. fluorouracil, cisplatin, and radiotherapy for carcinoma of the anal canal: A randomized controlled trial. JAMA 2008, 299, 1914–1921. [Google Scholar] [CrossRef]
  6. James, R.D.; Glynne-Jones, R.; Meadows, H.M.; Cunningham, D.; Myint, A.S.; Saunders, M.P.; Maughan, T.; McDonald, A.; Essapen, S.; Leslie, M.; et al. Mitomycin or cisplatin chemoradiation with or without maintenance chemotherapy for treatment of squamous-cell carcinoma of the anus (ACT II): A randomised, phase 3, open-label, 2 × 2 factorial trial. Lancet Oncol. 2013, 14, 516–524. [Google Scholar] [CrossRef]
  7. Shridhar, R.; Shibata, D.; Chan, E.; Thomas, C.R. Anal cancer: Current standards in care and recent changes in practice. CA Cancer. J. Clin. 2015, 65, 139–162. [Google Scholar] [CrossRef]
  8. Renehan, A.G.; Saunders, M.P.; Schofield, P.F.; O’Dwyer, S.T. Patterns of local disease failure and outcome after salvage surgery in patients with anal cancer. Br. J. Surg. 2005, 92, 605–612. [Google Scholar] [CrossRef]
  9. Northover, J.; Glynne-Jones, R.; Sebag-Montefiore, D.; James, R.; Meadows, H.; Wan, S.; Jitlal, M.; Ledermann, J. Chemoradiation for the treatment of epidermoid anal cancer: 13- year follow-up of the first randomised UKCCCR Anal Cancer Trial (ACT I). Br. J. Cancer 2010, 102, 1123–1128. [Google Scholar] [CrossRef] [Green Version]
  10. Peiffert, D.; Tournier-Rangeard, L.; Gérard, J.P.; Lemanski, C.; Francois, E.; Giovannini, M.; Cvitkovic, F.; Mirabel, X.; Bouche, O.; Luporsi, E.; et al. Induction chemotherapy and dose intensification of the radiation boost in locally advanced anal canal carcinoma: Final analysis of the randomized UNICANCER ACCORD 03 trial. J. Clin. Oncol. 2012, 30, 1941–1948. [Google Scholar] [CrossRef]
  11. Garg, M.K.; Zhao, F.; Sparano, J.A.; Palefsky, J.; Whittington, R.; Mitchell, E.P.; Mulcahy, M.F.; Armstrong, K.I.; Nabbout, N.H.; Kalnicki, S.; et al. Cetuximab Plus Chemoradiotherapy in Immunocompetent Patients with Anal Carcinoma: A Phase II Eastern Cooperative Oncology Group–American College of Radiology Imaging Network Cancer Research Group Trial (E3205). J. Clin. Oncol. 2017, 35, 718. [Google Scholar] [CrossRef] [PubMed]
  12. Sparano, J.A.; Lee, J.Y.; Palefsky, J.; Henry, D.H.; Wachsman, W.; Rajdev, L.; Aboulafia, D.; Ratner, L.; Fitzgerald, T.J.; Kachnic, L.; et al. Cetuximab Plus Chemoradiotherapy for HIV-Associated Anal Carcinoma: A Phase II AIDS Malignancy Consortium Trial. J. Clin. Oncol. 2017, 35, 727. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Sclafani, F.; Rao, S. Systemic Therapies for Advanced Squamous Cell Anal Cancer. Curr. Oncol. Rep. 2018, 20, 53. [Google Scholar] [CrossRef] [PubMed]
  14. Vendrely, V.; Lemanski, C.; Gnep, K.; Barbier, E.; Hajbi, F.E.; Lledo, G.; Dahan, L.; Terrebonne, E.; Manfredi, S.; Mirabel, X.; et al. Anti-epidermal growth factor receptor therapy in combination with chemoradiotherapy for the treatment of locally advanced anal canal carcinoma: Results of a phase I dose-escalation study with panitumumab (FFCD 0904). Radiother. Oncol. 2019, 140, 84–89. [Google Scholar] [CrossRef] [PubMed]
  15. Bian, J.J.; Almhanna, K. Anal cancer and immunotherapy-are we there yet? Transl. Gastroenterol. Hepatol. 2019, 4, 57. [Google Scholar] [CrossRef]
  16. Bibault, J.E.; Xing, L.; Giraud, P.; El Ayachy, R.; Giraud, N.; Decazes, P.; Burgun, A.; Giraud, P. Radiomics: A primer for the radiation oncologist. Cancer Radiother. 2020, 24, 403–410. [Google Scholar] [CrossRef]
  17. Lambin, P.; Leijenaar, R.T.H.; Deist, T.M.; Peerlings, J.; de Jong, E.E.C.; van Timmeren, J.; Sanduleanu, S.; Larue, R.T.H.M.; Even, A.J.G.; Jochems, A.; et al. Radiomics: The bridge between medical imaging and personalized medicine. Nat. Rev. Clin. Oncol. 2017, 14, 749–762. [Google Scholar] [CrossRef]
  18. Kourou, K. Machine learning applications in cancer prognosis and prediction. Comput. Struct. Biotechnol. J. 2015, 13, 8–17. [Google Scholar] [CrossRef] [Green Version]
  19. Schick, U.; Lucia, F.; Dissaux, G.; Visvikis, D.; Badic, B.; Masson, I.; Pradier, O.; Bourbonne, V.; Hatt, M. MRI-derived radiomics: Methodology and clinical applications in the field of pelvic oncology. BJR Br. J. Radiol. 2019, 92, 20190105. [Google Scholar] [CrossRef]
  20. Giraud, N.; Sargos, P.; Leduc, N.; Saut, O.; Vuong, T.; Vendrely, V. Radiomics and Machine Learning in Anal Squamous Cell Carcinoma: A New Step for Personalized Medicine? Appl. Sci. 2020, 10, 1988. [Google Scholar] [CrossRef] [Green Version]
  21. Bartelink, H.; Roelofsen, F.; Eschwege, F.; Rougier, P.; Bosset, J.F.; Gonzalez, D.G.; Peiffert, D.; Van Glabbeke, M.; Pierart, M. Concomitant radiotherapy and chemotherapy is superior to radiotherapy alone in the treatment of locally advanced anal cancer: Results of a phase III randomized trial of the European Organization for Research and Treatment of Cancer Radiotherapy and Gastrointestinal Cooperative Groups. J. Clin. Oncol. 1997, 15, 2040–2049. [Google Scholar] [PubMed]
  22. Myerson, R.J.; Kong, F.; Birnbaum, E.H.; Fleshman, J.W.; Kodner, I.J.; Picus, J.; Ratkin, G.A.; Read, T.E.; Walz, B.J. Radiation therapy for epidermoid carcinoma of the anal canal, clinical and treatment factors associated with outcome. Radiother. Oncol. J. Eur. Soc. Radiol. Oncol. 2001, 61, 15–22. [Google Scholar] [CrossRef]
  23. Faivre, J.C.; Peiffert, D.; Vendrely, V.; Lemanski, C.; Hannoun-Levi, J.M.; Mirabel, X.; Stanbury, T.; Salleron, J.; Guillemin, F. Prognostic factors of colostomy free survival in patients presenting with locally advanced anal canal carcinoma: A pooled analysis of two prospective trials (KANAL 2 and ACCORD 03). Radiother. Oncol. 2018, 129, 463–470. [Google Scholar] [CrossRef] [PubMed]
  24. Toh, E.; Wilson, J.; Sebag-Montefiore, D.; Botterill, I. Neutrophil:lymphocyte ratio as a simple and novel biomarker for prediction of locoregional recurrence after chemoradiotherapy for squamous cell carcinoma of the anus. Colorectal Dis. 2014, 16, O90–O97. [Google Scholar] [CrossRef] [PubMed]
  25. Wu, Y.; Han, X.; Li, Y.; Zhu, K.; Yu, J. Survival prediction models for patients with anal carcinoma receiving definitive chemoradiation: A population based study. Oncol. Lett. 2019, 19, 1443–1451. [Google Scholar] [CrossRef] [Green Version]
  26. Hocquelet, A.; Auriac, T.; Perier, C.; Dromain, C.; Meyer, M.; Pinaquy, J.B.; Denys, A.; Trillaud, H.; De Senneville, B.D.; Vendrely, V. Pre-treatment magnetic resonance-based texture features as potential imaging biomarkers for predicting event free survival in anal cancer treated by chemoradiotherapy. Eur. Radiol. 2018, 28, 2801–2811. [Google Scholar] [CrossRef]
  27. Owczarczyk, K.; Prezzi, D.; Cascino, M.; Kozarski, R.; Gaya, A.; Siddique, M.; Cook, G.J.; Glynne-Jones, R.; Goh, V. MRI heterogeneity analysis for prediction of recurrence and disease free survival in anal cancer. Radiother. Oncol. 2019, 134, 119–126. [Google Scholar] [CrossRef] [Green Version]
  28. Haralick, R.M.; Shanmugam, K.; Dinstein, I. Textural Features for Image Classification. IEEE Trans. Syst. Man. Cybern. 1973, SMC-3, 610–621. [Google Scholar] [CrossRef] [Green Version]
  29. Mehta, S.; Ramey, S.J.; Kwon, D.; Rich, B.J.; Ahmed, A.A.; Wolfson, A.; Yechieli, R.; Portelance, L.; Mellon, E.A. Impact of radiotherapy duration on overall survival in squamous cell carcinoma of the anus. J. Gastrointest. Oncol. 2020, 11, 277–290. [Google Scholar] [CrossRef]
  30. del Campo, E.R.; Matzinger, O.; Haustermans, K.; Peiffert, D.; Glynne-Jones, R.; Winter, K.A.; Konski, A.A.; Ajani, J.A.; Bosset, J.F.; Hannoun-Levi, J.M.; et al. Pooled Analysis of external-beam RADiotherapy parameters in phase II and phase III trials in radiochemotherapy in Anal Cancer (PARADAC). Eur. J. Cancer 2019, 121, 130–143. [Google Scholar] [CrossRef]
  31. Serup-Hansen, E.; Linnemann, D.; Skovrider-Ruminski, W.; Høgdall, E.; Geertsen, P.F.; Havsteen, H. Human Papillomavirus Genotyping and p16 Expression as Prognostic Factors for Patients with American Joint Committee on Cancer Stages I to III Carcinoma of the Anal Canal. JCO J. Clin. Oncol. 2014, 32, 1812–1817. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Gilbert, D.C.; Williams, A.; Allan, K.; Stokoe, J.; Jackson, T.; Linsdall, S.; Bailey, C.M.; Summers, J. p16INK4A, p53, EGFR expression and KRAS mutation status in squamous cell cancers of the anus: Correlation with outcomes following chemo-radiotherapy. Radiother. Oncol. 2013, 109, 146–151. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Jones, C.M.; Goh, V.; Sebag-Montefiore, D.; Gilbert, D.C. Biomarkers in anal cancer: From biological understanding to stratified treatment. Br. J. Cancer 2017, 116, 156–162. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Zou, K.H.; Warfield, S.K.; Bharatha, A.; Tempany, C.M.; Kaus, M.R.; Haker, S.J.; Wells, W.M., 3rd; Jolesz, F.A.; Kikinis, R. Statistical validation of image segmentation quality based on a spatial overlap index. Acad. Radiol. 2004, 11, 178–189. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Radiomics workflow and integration into Machine Learning models.
Figure 1. Radiomics workflow and integration into Machine Learning models.
Cancers 13 00193 g001
Figure 2. Flowchart of patient selection. (Abbreviations: CRT = chemo-radiotherapy, NM = non-metastatic, ASCC = anal squamous cell cancer, CT = chemotherapy, MRI = magnetic resonance imaging).
Figure 2. Flowchart of patient selection. (Abbreviations: CRT = chemo-radiotherapy, NM = non-metastatic, ASCC = anal squamous cell cancer, CT = chemotherapy, MRI = magnetic resonance imaging).
Cancers 13 00193 g002
Figure 3. Correlation matrix between features significantly associated with the outcome in the training cohort. (Abbreviations: JEnergy = Glcm_JointEnergy, JEntropy = Glcm_JointEntropy, GLVariance = Glrlm_GrayLevelVariance, CRTDuration = Chemo-radiotherapy duration).
Figure 3. Correlation matrix between features significantly associated with the outcome in the training cohort. (Abbreviations: JEnergy = Glcm_JointEnergy, JEntropy = Glcm_JointEntropy, GLVariance = Glrlm_GrayLevelVariance, CRTDuration = Chemo-radiotherapy duration).
Cancers 13 00193 g003
Figure 4. Kaplan–Meier estimates of disease-free survival using the combined (clinical and radiomics) model for (A) training and (B) testing set, and ROC curve in the testing set (C).
Figure 4. Kaplan–Meier estimates of disease-free survival using the combined (clinical and radiomics) model for (A) training and (B) testing set, and ROC curve in the testing set (C).
Cancers 13 00193 g004
Table 1. Patient and tumor characteristics in the two subsequently blended cohorts. (Abbreviations: F = female, M = male, HIV = human immunodeficiency virus, DFS = disease-free survival).
Table 1. Patient and tumor characteristics in the two subsequently blended cohorts. (Abbreviations: F = female, M = male, HIV = human immunodeficiency virus, DFS = disease-free survival).
Bordeaux Cohort (n = 36)FFCD Cohort (n = 46)p-Value
Age at diagnosis (mean, years)62.360.10.39
Sex: 0.96
-  F30/36 (83.3%)37/46 (80.4%)
-  M6/36 (16.7%)9/46 (19.6%)
HIV status: 0.79
-  Positive1/36 (2.8%)3/46 (6.5%)
-  Negative35/36 (97.2%)43/46 (93.5%)
Tumor length (mean, mm)51.247.50.40
Pre-rectal tumor infiltration16/36 (44.4%)13/46 (28.2%)0.20
T-stage: 0.44
-  cT1-211/36 (30.6%)19/46 (41.2%)
-  cT3-425/36 (69.4%)27/46 (58.7%)
N-stage: 1
-  cN010/36 (27.8%)12/46 (26.1%)
-  cN+26/36 (72.2%)34/46 (73.9%)
AJCC groupings: 0.20
-  I1/36 (2.8%)0/46 (0%)
-  IIa4/36 (11.1%)5/46 (10.9%)
-  IIb1/36 (2.8%)4/46 (8.7%)
-  IIIa6/36 (16.7%)14/46 (30.4%)
-  IIIb4/36 (11.1%)3/46 (6.5%)
-  IIIc20/36 (55.5%)20/46 (43.5%)
2-year DFS23/36 (63.9%)35/46 (76.1%)0.34
Table 2. Signification of the radiomics features significantly associated with the outcome in the training cohort.
Table 2. Signification of the radiomics features significantly associated with the outcome in the training cohort.
FO_Uniformity i = 1 N g p ( i ) 2 Sum of the squares of each intensity value. A greater uniformity implies a greater homogeneity or a smaller range of discrete intensity values.
FO_Entropy i = 1 N g p ( i ) l o g 2 ( p ( i ) + ε ) Randomness in the image values, with ϵ an arbitrarily small positive number (≈2.2 × 10−16).
GLCM_JointEnergy i = 1 N g j = 1 N g ( p ( i , j ) 2 ) Measure of homogeneous patterns in the image. A greater energy implies more intensity value pairs that neighbor each other at higher frequencies.
GLCM_JointEntropy i = 1 N g j = 1 N g p ( i , j ) l o g 2 ( p ( i , j ) + ε ) Measure of the variability in neighborhood intensity values.
GLRLM_GrayLevelVariance i = 1 N g j = 1 N s p ( i , j ) ( i μ ) 2 Variance in gray level intensities for the zones.
Table 3. Correlation between variables and 2-year disease recurrence. (Abbreviations: Se = sensitivity, Sp = specificity).
Table 3. Correlation between variables and 2-year disease recurrence. (Abbreviations: Se = sensitivity, Sp = specificity).
VariableUnivariate AnalysisBest Cut-Offp-ValueMultivariate Analysis
AUCSe (%)Sp (%)Odds-Ratiop-Value
Tumor length (mm)0.71 (0.56–0.86)7165490.015.750.036
CRT duration (days)0.67 (0.52–0.83)825150.50.043
FO_Entropy0.67 (0.51–0.84)59870.9920.0450.120.026
GLCM_JointEnergy0.67 (0.50–0.84)41950.3170.047
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Giraud, N.; Saut, O.; Aparicio, T.; Ronchin, P.; Bazire, L.-A.; Barbier, E.; Lemanski, C.; Mirabel, X.; Etienne, P.-L.; Lièvre, A.; et al. MRI-Based Radiomics Input for Prediction of 2-Year Disease Recurrence in Anal Squamous Cell Carcinoma. Cancers 2021, 13, 193. https://doi.org/10.3390/cancers13020193

AMA Style

Giraud N, Saut O, Aparicio T, Ronchin P, Bazire L-A, Barbier E, Lemanski C, Mirabel X, Etienne P-L, Lièvre A, et al. MRI-Based Radiomics Input for Prediction of 2-Year Disease Recurrence in Anal Squamous Cell Carcinoma. Cancers. 2021; 13(2):193. https://doi.org/10.3390/cancers13020193

Chicago/Turabian Style

Giraud, Nicolas, Olivier Saut, Thomas Aparicio, Philippe Ronchin, Louis-Arnaud Bazire, Emilie Barbier, Claire Lemanski, Xavier Mirabel, Pierre-Luc Etienne, Astrid Lièvre, and et al. 2021. "MRI-Based Radiomics Input for Prediction of 2-Year Disease Recurrence in Anal Squamous Cell Carcinoma" Cancers 13, no. 2: 193. https://doi.org/10.3390/cancers13020193

APA Style

Giraud, N., Saut, O., Aparicio, T., Ronchin, P., Bazire, L. -A., Barbier, E., Lemanski, C., Mirabel, X., Etienne, P. -L., Lièvre, A., Cacheux, W., Darut-Jouve, A., De la Fouchardière, C., Hocquelet, A., Trillaud, H., Charleux, T., Breysacher, G., Argo-Leignel, D., Tessier, A., ... Vendrely, V. (2021). MRI-Based Radiomics Input for Prediction of 2-Year Disease Recurrence in Anal Squamous Cell Carcinoma. Cancers, 13(2), 193. https://doi.org/10.3390/cancers13020193

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop