Early prediction of radiotherapy-induced parotid shrinkage and toxicity based on CT radiomics and fuzzy classification

MOTIVATION Patients under radiotherapy for head-and-neck cancer often suffer of long-term xerostomia, and/or consistent shrinkage of parotid glands. In order to avoid these drawbacks, adaptive therapy can be planned for patients at risk, if the prediction is obtained timely, before or during the early phase of treatment. Artificial intelligence can address the problem, by learning from examples and building classification models. In particular, fuzzy logic has shown its suitability for medical applications, in order to manage uncertain data, and to build transparent rule-based classifiers. In previous works, clinical, dosimetric and image-based features were considered separately, to find different possible predictors of parotid shrinkage. On the other hand, a few works reported possible image-based predictors of xerostomia, while the combination of different types of features has been little addressed. OBJECTIVE This paper proposes the application of a novel machine learning approach, based on both statistics and fuzzy logic, aimed at the classification of patients at risk of i) parotid gland shrinkage and ii) 12-months xerostomia. Both problems are addressed with the aim of individuating predictors and models to classify respective outcomes. METHODS Knowledge is extracted from a real dataset of radiotherapy patients, by means of a recently developed method named Likelihood-Fuzzy Analysis, based on the representation of statistical information by fuzzy rule-based models. This method enables to manage heterogeneous variables and missing data, and to obtain interpretable fuzzy models presenting good generalization power (thus high performance), and to measure classification confidence. Numerous features are extracted to characterize patients, coming from different sources, i.e. clinical features, dosimetric parameters, and radiomics-based measures obtained by texture analysis of Computed Tomography images. A learning approach based on the composition of simple models in a more complicated one allows to consider the features separately, in order to identify predictors and models to use when only some data source is available, and obtaining more accurate results when more information can be combined. RESULTS Regarding parotid shrinkage, a number of good predictors is detected, some already known and confirmed here, and some others found here, in particular among radiomics-based features. A number of models are also designed, some using single features and others involving models composition to improve classification accuracy. In particular, the best model to be used at the initial treatment stage, and another one applicable at the half treatment stage are identified. Regarding 12-months toxicity, some possible predictors are detected, in particular among radiomics-based features. Moreover, the relation between final parotid shrinkage rate and 12-months xerostomia is evaluated. The method is compared to the naïve Bayes classifier, which reveals similar results in terms of classification accuracy and best predictors. The interpretable fuzzy rule-based models are explicitly presented, and the dependence between predictors and outcome is explained, thus furnishing in some cases helpful insights about the considered problems. CONCLUSION Thanks to the performance and interpretability of the fuzzy classification method employed, predictors of both parotid shrinkage and xerostomia are detected, and their influence on each outcome is revealed. Moreover, models for predicting parotid shrinkage at initial and half radiotherapy stages are found.

[1]  Motoyasu Nakamura,et al.  Changes in parotid gland morphology and function in patients treated with intensity-modulated radiotherapy for nasopharyngeal and oropharyngeal tumors , 2013, Oral Radiology.

[2]  Giovanna Rizzo,et al.  Texture analysis for the assessment of structural changes in parotid glands induced by radiotherapy. , 2013, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[3]  M. Hoogeman,et al.  Local anatomic changes in parotid and submandibular glands during radiotherapy for oropharynx cancer and correlation with dose, studied in detail with nonrigid registration. , 2008, International journal of radiation oncology, biology, physics.

[4]  Toshinori Hirai,et al.  Radiation-induced parotid gland changes in oral cancer patients: correlation between parotid volume and saliva production. , 2010, Japanese journal of clinical oncology.

[5]  R Calandrino,et al.  An automatic contour propagation method to follow parotid gland deformation during head-and-neck cancer tomotherapy , 2011, Physics in medicine and biology.

[6]  Toshinori Hirai,et al.  Longitudinal changes over 2 years in parotid glands of patients treated with preoperative 30-Gy irradiation for oral cancer. , 2011, Japanese journal of clinical oncology.

[7]  Giovanna Rizzo,et al.  Early classification of parotid glands shrinkage in radiotherapy patients: A comparative study , 2015 .

[8]  Paul Kinahan,et al.  Radiomics: Images Are More than Pictures, They Are Data , 2015, Radiology.

[9]  Giovanna Rizzo,et al.  Introducing the Jacobian-volume-histogram of deforming organs: application to parotid shrinkage evaluation , 2011, Physics in medicine and biology.

[10]  Massimo Esposito,et al.  A Composite Model for Classifying Parotid Shrinkage in Radiotherapy Patients Using Heterogeneous Data , 2015, AIME.

[11]  Lotfi A. Zadeh,et al.  Fuzzy Sets , 1996, Inf. Control..

[12]  L. Dawson,et al.  Xerostomia and its predictors following parotid-sparing irradiation of head-and-neck cancer. , 2001, International journal of radiation oncology, biology, physics.

[13]  Nir Friedman,et al.  Bayesian Network Classifiers , 1997, Machine Learning.

[14]  A. Fiorentino,et al.  Parotid gland volumetric changes during intensity-modulated radiotherapy in head and neck cancer. , 2012, The British journal of radiology.

[15]  Kevin J Harrington,et al.  Weekly volume and dosimetric changes during chemoradiotherapy with intensity-modulated radiation therapy for head and neck cancer: a prospective observational study. , 2010, International journal of radiation oncology, biology, physics.

[16]  G Rizzo,et al.  The Shape of Parotid DVH Predicts the Entity of Gland Deformation During IMRT for Head and Neck Cancers , 2015, Technology in cancer research & treatment.

[17]  John H. Holmes,et al.  Artificial Intelligence in Medicine , 2015, Lecture Notes in Computer Science.

[18]  Radhe Mohan,et al.  Quantification of volumetric and geometric changes occurring during fractionated radiotherapy for head-and-neck cancer using an integrated CT/linear accelerator system. , 2004, International journal of radiation oncology, biology, physics.

[19]  Joseph O Deasy,et al.  Radiotherapy dose-volume effects on salivary gland function. , 2010, International journal of radiation oncology, biology, physics.

[20]  Giovanna Rizzo,et al.  Early changes of parotid density and volume predict modifications at the end of therapy and intensity of acute xerostomia , 2014, Strahlentherapie und Onkologie.

[21]  Liangxiao Jiang,et al.  Learning Tree Augmented Naive Bayes for Ranking , 2005, DASFAA.

[22]  Massimo Esposito,et al.  Likelihood-Fuzzy Analysis of Parotid Gland Shrinkage in Radiotherapy Patients , 2014, InMed.

[23]  Ian H. Witten,et al.  The WEKA data mining software: an update , 2009, SKDD.

[24]  Giovanna Rizzo,et al.  Density variation of parotid glands during IMRT for head-neck cancer: correlation with treatment and anatomical parameters. , 2012, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[25]  Toshinori Hirai,et al.  Histopathological changes in parotid and submandibular glands of patients treated with preoperative chemoradiation therapy for oral cancer. , 2012, Journal of radiation research.

[26]  Marta Paiusco,et al.  A two-variable linear model of parotid shrinkage during IMRT for head and neck cancer. , 2009, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[27]  Binbin Wu,et al.  Parotid gland shrinkage during IMRT predicts the time to Xerostomia resolution , 2015, Radiation Oncology.

[28]  Giuseppe De Pietro,et al.  Best Fuzzy Partitions to Build Interpretable DSSs for Classification in Medicine , 2013, HAIS.