Impact Factor 4.137 | CiteScore 4.28
More on impact ›

Original Research ARTICLE

Front. Oncol., 04 February 2020 | https://doi.org/10.3389/fonc.2020.00077

Prediction of Response to Preoperative Neoadjuvant Chemotherapy in Locally Advanced Cervical Cancer Using Multicenter CT-Based Radiomic Analysis

Xin Tian1, Caixia Sun2,3, Zhenyu Liu2,4, Weili Li1, Hui Duan1, Lu Wang1, Huijian Fan1, Mingwei Li1, Pengfei Li1, Lihui Wang3*, Ping Liu1*, Jie Tian2,4,5* and Chunlin Chen1*
  • 1Department of Gynaecology and Obstetrics, Nanfang Hospital, Southern Medical University, Guangzhou, China
  • 2CAS Key Laboratory of Molecular Imaging, Institute of Automation, Chinese Academy of Sciences, Beijing, China
  • 3Key Laboratory of Intelligent Medical Image Analysis and Precise Diagnosis of Guizhou Province, School of Computer Science and Technology, Guizhou University, Guiyang, China
  • 4School of Artificial Intelligence, University of Chinese Academy of Sciences, Beijing, China
  • 5Beijing Advanced Innovation Center for Big Data-Based Precision Medicine, Beihang University, Beijing, China

Objective: To investigate whether pre-treatment CT-derived radiomic features could be applied for prediction of clinical response to neoadjuvant chemotherapy (NACT) in locally advanced cervical cancer (LACC).

Patients and Methods: Two hundred and seventy-seven LACC patients treated with NACT followed by surgery/radiotherapy were included in this multi-institution retrospective study. One thousand and ninety-four radiomic features were extracted from venous contrast enhanced and non-enhanced CT imaging for each patient. Five combined methods of feature selection were used to reduce dimension of features. Radiomics signature was constructed by Random Forest (RF) method in a primary cohort of 221 patients. A combined model incorporating radiomics signature with clinical factors was developed using multivariable logistic regression. Prediction performance was then tested in a validation cohort of 56 patients.

Results: Radiomics signature containing pre- and post-contrast imaging features can adequately distinguish chemotherapeutic responders from non-responders in both primary and validation cohorts [AUCs: 0.773 (95% CI, 0.701–0.845) and 0.816 (95% CI, 0.690-0.942), respectively] and remain relatively stable across centers. The combined model has a better predictive performance with an AUC of 0.803 (95% CI, 0.734–0.872) in the primary set and an AUC of 0.821 (95% CI, 0.697–0.946) in the validation set, compared to radiomics signature alone. Both models showed good discrimination, calibration.

Conclusion: Newly developed radiomic model provided an easy-to-use predictor of chemotherapeutic response with improved predictive ability, which might facilitate optimal treatment strategies tailored for individual LACC patients.

Introduction

Locally advanced cervical cancer (LACC) suffer from high risks of treatment failure (1). Neoadjuvant chemotherapy (NACT) followed by radical hysterectomy is increasingly applied as alternative to standard chemoradiotherapy for LACC in some countries as NACT could shrink tumor volume and render unresectable tumors operable (2, 3). Recent evidence further revealed comparable oncological outcome between NACT-surgery and chemoradiotherapy for stage IB2-IIB cervical cancer, indicating optimal treatment decision relied upon quality of life (4). NACT is beneficial in controlling micro-metastasis and potentially improves patients' quality of life by preserving sexual capacity and ovarian function (5, 6). Moreover, pathological optimal responders could even exhibit survival benefit without need for post-operative adjuvant treatments (7, 8). However, chemotherapeutic non-responders suffer from disease progression and experience worse prognosis due to delay in definite therapies (9, 10). Therefore, clinical use of NACT for LACC should be more precise and individualized. Only chemotherapeutic responders can be candidates for preoperative NACT, whereas non-responders should directly undergo chemoradiotherapy avoiding ineffective therapies. Therefore, an overwhelming need to discover pretherapy predictors of NACT response to stratify LACC patients to appropriate therapies.

Immunohistochemical markers identified by measuring alterations in protein expression from tumor tissue samples (pre- and post-chemotherapy) are relatively unreliable due to lack of validation and conflicting results in various studies (11, 12). Some proposed genomic, proteomic predictive signature are limited by small sample size and high cost of examination (13, 14). Several pioneer studies demonstrated quantitative parameters derived from pre- and post-treatment functional imaging of LACC patients were powerful markers to non-invasively predict early therapeutic response to NACT (15, 16). These quantitative parameters allow for characterizing tumor biological process during NACT. However, post-treatment nature limits its extensive utility in therapy decision.

Radiomics, an emerging approach for image interpretation, can non-invasively predict treatment response solely by radiographic examination (17). Radiomics can profile tumor heterogeneity of different responders by extracting extensively quantitative features that reflect underlying pathophysiology from medical images (18). In fact, radiomics has proved to perform well in prognosis prediction of cervical cancer treated with chemoradiotherapy using PET/CT and MRI images (19). Moreover, pretherapy computed tomography (CT) radiomics also managed to predict response to neoadjuvant chemoradiotherapy for locally advanced rectal cancer (20). Furthermore, Li et al. pilot study successfully applied pre-treatment CT radiomics in predicting pathological response after neoadjuvant chemotherapy in locally advanced gastric cancers (21). In light of this, radiomics is promising for pre-therapy prediction of response to NACT for LACC.

CT is widely used for pre-treatment staging of cervical cancer (22, 23). Investigation regarding the feasibility of pretreatment CT-based radiomics for prediction of NACT response in LACC has not yet reported. Meantime, limited reproducibility in multi-center context remain a crucial bottleneck in widespread use of radiomics in clinic. Thus, the aim of the study is to establish and validate an effective pretherapy CT-based radiomic model for predicting response to NACT while testify stability and reproducibility of radiomic features and models on a multi-center basis.

Materials and Methods

Patients

This retrospective study was approved by institutional review board of Nanfang Hospital, Southern Medical University with informed consent obtained from all participants. The study was conducted according to the Declaration of Helsinki and was based on data collected in 10 centers from January 2009 to June 2018. In total, 277 consecutive patients with locally advanced cervical cancer underwent neoadjuvant chemotherapy with pre-treatment CT data were included (Figure 1). Eligibility criteria include: (1) Histologically diagnosis of squamous cell carcinoma of the uterine cervix. (2) Patients with clinical stage IB2–III cervical cancer (FIGO 2009) (24). (3) Initial treatment with NACT followed by radical surgery or radiotherapy without any previous therapy. (4) Pre-treatment pelvic CT images. Exclusion criteria are as follows: (1) Underwent preoperative neo-adjuvant concurrent radiotherapy together with chemotherapy. (2) Patients had other malignancy besides cervical cancer. (3) Patients lacking all necessary CT sequences, either venous phase or non-contrast CT images. (4) Patients received only one cycle of preoperative chemotherapy but had insufficient tumor shrinkage that recognized as “stable disease” by RECIST1.1 Criteria. The enrolled 277 patients were randomly divided into a primary cohort of 221 patients (mean age, 48.75 ± 7.94 years; range, 27–69 years) and an independent validation cohort of 56 patients (mean age, 46.36 ± 9.29 years; range, 27–67 years). Both cohorts have comparable ratios of non-responders versus responders. Baseline clinic-pathologic data including age, FIGO stage, HPV infection status, tumor size, and level of SCC antigen before/after chemotherapy were derived from medical records.

FIGURE 1
www.frontiersin.org

Figure 1. Patient population and exclusions.

Neoadjuvant Chemotherapy

All the patients received preoperative neoadjuvant chemotherapy intra-arterially or intravenously. Pre-treatment CT examinations were performed within 1 week before the first cycle of NACT. Tumor volume was carefully assessed by gynecological examination together with either transvaginal ultrasound, pelvic MRI, or CT scan at the start of each chemotherapy cycle and 3 weeks after treatment. Final tumor size was confirmed in postoperative tumor pathology. Compared to pre-treatment evaluation, patients with sufficient tumor shrinkage to qualify for operation underwent radical hysterectomy after one cycle of chemotherapy. Patients with insufficient tumor shrinkage for operation had additional 1–2 cycles of chemotherapy at 21-days intervals. Surgery was performed within 21–28 days after the completion of the last cycle. Radiotherapy was performed for patients who had progressive diseases or adverse side effects due to the toxicity of anti-tumor agents during or after NACT.

At the present study, multicentric patients had distinct courses of NACT. Thirty-six patients received only one cycle of intra-arterial chemotherapy, which cisplatin (60–75 mg/m2)/carboplatin (300 mg/2)+bleomycin (45 mg) was the most adopted program(61.1%). Two hundred and forty-one patients underwent platinum-based combination chemotherapy intravenously, consisting of 1 (103 cases), 2 (124 cases), or 3 (14 cases) courses of treatment. Two hundred and eleven of Two hundred and forty-one (87.5%) patients were treated with paclitaxel (135–175 mg/m2) or docetaxel (60–70 mg/m2) plus any kinds of platinum regimens. Details are presented in Table A6.

Assessment of Response

The short-term response was evaluated by change in tumor size according to Response Evaluation Criteria In Solid Tumors (RECIST v. 1.1) as follows: complete response (CR) was defined as eradication of cervical lesion; partial response (PR) (a ≥30% decrease in the longest tumor diameter); progressive disease (PD) (a ≥20% increase in the longest diameter (LD) of tumor) and stable disease (SD) as the decrease or increase of LD of the cervical lesion was less than PR or PD (25). Many studies proposed two cycles of NACT were adequate for LACC patients to obtain optimal efficacy with relatively low adverse reactions (9, 26). Based on various treatment courses in the study, we set up new response evaluation criteria especially for upcoming radiomic analysis as following: For patients with CR or PR by RECIST (v.1.1) within two cycles of treatment are recognized as “responder” (chemo-sensitive). Whereas, those still with PD or SD after two cycles of chemotherapy were deemed as “non-responder” (chemo-insensitive) regardless of final treatment cycles. Thus, there were 201 responders and 76 non-responders among the qualified 277 consecutive patients (Table A1).

CT Feature Extraction and Selection

Pre-treatment non-contrast and venous contrasted enhanced CT images at 1.5–3 mm thickness for each case were retrieved from the image-archiving workstation. The whole tumor volume of each patient was manually segmented as region of interest (ROI) by a radiologist via the ITK-SNAP software (www.itksnap.org). As contour of large tumor mass were better visible on contrast-enhanced images. ROI were first drawn on each transverse slice of pelvic venous phase images then applied to non-contrast images. The reproducibility of extracted features from ROIs were ensured through inter- and intra-observer reproducibility evaluation (detailed in the Appendix 1).

In order to get a standard normal distribution of image intensities, each slice of CT images was normalized with z-score. Features were extracted from ROIs of non-contrast and venous enhanced CT using MATLAB R2016b software, respectively. Four groups of features were extracted for each modality: 17 first order statistical features, 8 shape features, 54 texture features and 568 wavelet features (Appendix 2; Table A3). The extracted features are reproducible and match the benchmarks of IBSI (27, 28). Ultimately, a total of 1,294 features were extracted from each patient. Each feature of patients was also normalized with z-score to remove the effect of different CT scanners.

The strategy of feature selection was as follow: Firstly, we selected features with an inter-class correlation coefficients (ICC) greater than 0.85, which usually have good reproducibility. Then, five methods were used to further select predictive features of response, including recursive feature elimination based on a support vector machine (SVM-RFE), least absolute shrinkage and selection operator (LASSO), extremely randomized trees (ET), random forest (RF), and ridge regression. The parameters of each method were determined by grid search with five-fold cross-validation. Ultimately, we selected ones that were repeatedly appeared as the most predictive features among all five methods to be the optimal radiomic features. Feature selection was first performed in the primary cohort and stored the name of the selected features. Then, we selected the features with stored names in validation cohort.

Building and Validation of Radiomic Signature

For constructing radiomics signature, a RF model was used to predict NACT response. The model was trained on the primary cohort, and parameters of the model were set by grid search with five-fold cross-validation. In the process of training the model, the synthetic minority over-sampling technique (SMOTE) was applied to expand the ratio between responders and non-responders to 1, so that the model was not affected by the sample rate and does not bias to either side during the learning experience. The association of radiomic features with the response of NACT was first assessed in primary cohort, and then validated in validation cohort.

Construction and Validation of Combined Model

To verify whether prediction performance will be improved when adding clinical information, a multivariable logistic regression model was constructed combining radiomics signature with clinical factors. Clinical candidates include age, FIGO stage and maximum diameter of tumor. Akaikes information criterion (AIC) was adopted to select optimal clinical factors.

After completing model construction, a nomogram was employed to provide a quantitative tool to predict individualized probability of response. For assessing the performance of the model, calibration curves were applied to assess capabilities in primary and validation cohorts, where the curves were used to describe the difference in predictive and actual probability of response. Hosmer-Lemeshow test was adopted to calculate p-value to represent the fit degree of predictive and actual curves. Besides, clinical usefulness of the combined model is determined by decision curve analysis. It is created by plotting net benefit at various threshold probability.

For comparing the performance between the combined model and radiomics signature, Integrated Discrimination Improvement (IDI) was calculated to quantificationally describe the increment of accuracy between the two models. Study flowchart is displayed in Figure 2A.

FIGURE 2
www.frontiersin.org

Figure 2. A schema for radiomics pipeline. (A) Flowchart of the study. One thousand and twenty-four features were extracted from pre-therapy CT scans. Both extracted radiomic features of non-contrast and venous enhanced images are pooled as part of following feature selection analysis. After chosen by five different feature selection methods, the final top predictive features were then constructed as radiomics signature by Random Forest. Finally, radiomics signature and clinical factors were combined into a nomogram. (B) Difference of quantitative radiomic features and predictive probability between chemotherapy responder vs. non-responder. Left images are their corresponding pre-therapeutic CT non-enhanced and venous enhanced images. ROIs are drawn in a purple circle. Right graph indicated the value of radiomic features and predictive probability derived from radiomics signature.

Statistical Analysis

The process of statistical analysis was performed on R software (version 3.5.1; https://www.r-project.org). P < 0.05 was considered to indicate a statistically significant difference and the p-value was bilateral.

Results

Clinical Characteristics

Patient characteristics in primary and validation cohorts are given in Table 1. Of the 277 enrolled patients, there were 201 responders (72.56%) and 76 non-responders (27.43%). Median age of the patients was 48 [40–56] years. There are no significant differences between the two cohorts in age and maximum tumor diameter. Primary cohort consisted of 160 responders (72.40%) and 61 non-responders (27.60%) while validation cohort had 41 responders (73.21%) and 15 non-responders (26.79%). CT based LN metastasis positivity was 28.75 and 48.78% in primary and validation cohorts, respectively. Each cohort has a relative equal ratio of different FIGO stages.

TABLE 1
www.frontiersin.org

Table 1. Characteristics of patients in the primary and validation cohorts.

Feature Selection

A total of 1294 features were extracted from non-contrast and venous enhanced CT images. Eight hundred and sixty-three features of ICC > 0.85 was reserved for further feature selection. The details of selected features by five methods were presented in Table A4. Finally, six features that repeatedly appeared in five series of the subset of features were selected (Table A5). The comparison of these features between different centers was shown in Figure 3. Because the number of patients is scattered among 10 hospitals, we only selected two hospitals with the largest number of patients for comparison. Results revealed that there are no significant differences in the distribution of five characteristics extracted from non-contrast CT images between hospitals (P > 0.05) while the distribution of “HHL_glrlm_SRHGLE” feature derived from venous-enhanced images showed slight difference (P = 0.02).

FIGURE 3
www.frontiersin.org

Figure 3. Probability density function of features between centers. Comparison of the six derived features between two hospitals with largest numbers of patients revealed no significant difference exist in the five features from non-contrast CT images (glcm_correlation, LLH_glcm_entropy, HLL_glrlm_GLN, LHH_glcm_IDMN, glcm_homogeneity1) (P > 0.05). “HHL_glrlm_SRHGLE” features from venous-enhanced images showed slight difference. P values are for Mann-Whitney U-test.

Construction and Validation of Radiomic Signature

An RF model with Gini criterion was constructed using the top features. The performance of model was first assessed in primary cohort, and then validated in the validation cohort. As shown in Table 2, radiomics signature yielded good performance with an AUC of 0.773 (95% CI, 0.701–0.845) and an accuracy of 78.28% (95% CI, 73.30–83.26%) in primary cohort while achieved an AUC of 0.816 (95% CI, 0.690–0.942) and an accuracy of 80.36% (95% CI, 69.64–89.29%) in validation cohort. There were difference in corresponding quantitative values of radiomic features between non-responders and responders extracted from CT images (Figure 2B).

TABLE 2
www.frontiersin.org

Table 2. Performance of radiomics signature and the combined model.

Construction and Validation of Combined Model

The combined model was constructed incorporating radiomics signature, age and FIGO stage while presented as a nomogram (Figure 4A). The calibration curves of the combined model in both primary and validation cohorts showed a good fit between predictive probability of response and actual response rate (Figure 4B). Non-significant statistics of the accompanied Hosmer-Lemeshow test (primary cohort: p = 0.396; validation cohort: p = 0.604) implied the model was adequately calibrated without departure from the ideal model. Decision curve demonstrated that the application of the model in clinical decision-making could achieve greater benefits (Figure 4C).

FIGURE 4
www.frontiersin.org

Figure 4. Nomogram developed with the combined model and calibration curves, decision curve analysis for the combined model. (A) The developed nomogram. (B) Calibration curve of the combined model in the primary and validation cohorts. The middle gray line represents a perfect prediction. The blue line represents the performance of the combined model. Better prediction is demonstrated by a closer fit of the blue line to the grayline. (C) Decision curve analysis for the combined model. The y-axis depicts the net benefit. The blue line represents the combined model. The grayline represents the assumption that all patients have response of NACT. The black line is the opposite.

Comparison of Radiomics Signature and Combined Model

The combined model achieved adequate discrimination performance with AUC of 0.803 (95% CI, 0.734–0.872) and AUC of 0.821 (95% CI, 0.697–0.946) in the primary and validation cohorts respectively. In terms of predictive performance, the combined model yielded significant improvement compared to radiomics signature alone (Figure 5; Table 2). Besides, the classification accuracy of the combined model estimating response probability was also superior to sole radiomics signature in both the primary [80.54% (95% CI, 76.02–85.07%)] and validation cohorts [82.14% (95% CI, 73.21–89.29%)], respectively.

FIGURE 5
www.frontiersin.org

Figure 5. Receiver Operating Curve (ROC) of the models in each cohort. (A) ROC in primary cohort. (B) ROC in the validation cohort. The middle gray curve represents the dividing line with AUC of 0.5.

Discussion

In the present study, we successfully established and validated a radiomic signature by pooling together non-enhanced and venous enhanced CT imaging features for predicting response to NACT in LACC using multicenter data set. This six-features radiomics signature can discriminate responders in both the primary (AUC: 0.773) and validation cohort (AUC: 0.816). Moreover, the combined model incorporating radiomics signature and two clinical factors (age, FIGO stage) achieved better performance than radiomic signature alone.

Genomic heterogeneity affect treatment response in LACC. Different responders of neoadjuvant chemotherapy in LACC exhibit diverse genomic profile (29). Twenty-seven-predictive gene signature has been identified to be qualified for distinguishing responders from non-responders in LACC patients treated with concurrent chemoradiotherapy (30). Previous research argued that radiomics could quantify intra-tumor heterogeneity, reflecting underlying tumor gene-expression patterns through the spatial arrangement of imaging voxels (31). Indeed, in the presented study, the established radiomics signature was comprised of six key features describing intratumor voxel patterns (textures) that non-invasively captured intratumor heterogeneity of different responders. “Gray-level co-occurrence matrix” (GLCM) -correlation is a measure of linear dependency of grayscale while GLCM-homogeneity measures local gray level uniformity of the image. As shown in the Figure 2B, local gray level of responder's images is relatively uniform with a larger value of GLCM-homogeneity while signal-intensity variations seem more apparent in non-responder images. “glrlm_GLN” feature is informative for gray level non-uniformity, which shows significant higher value in chemotherapeutic non-responder (32). Note that radiomic features (LHH_glcm_IDMN, LLH_glcm_entropy, HHL_glrlm_SRHGLE) were generated from decomposing original images through using a coiflet wavelet transformation, which at multiple scales could further reflect tumor spatial heterogeneity. This quantitative radiomic features may reflect tumor hetergenous histopathological patterns that are associated with different treatment response such as intratumoral necrosis, vascularization, cellular density (33).

In addition, in order to evaluate the stability and consistency of the established radiomic model, we combined patients from different hospitals to construct two extra radiomic signature as independent validation (Table A2). Only <4% difference was observed between two groups in both training and validation sets (no significant power). Meantime both groups showed good predictive performance, consistent to the model established in the article (Figure A1). This result revealed that our model performed stable between different centers, and it was not affected by patient composition of the training set. Moreover, the reproducibility of derived radiomic features was also robust on a multi-center basis, backed by evidence that after comparison between two hospitals with largest numbers of patients, five of the six features remained stable. Noticeably, the features extracted from non-enhanced CT images were still reliable despite scan-related parameters variation.

Evaluated by AIC methods, adding FIGO stage to radiomics signature enhanced predictive power whereas primary tumor maximal diameter failed. Strikingly, age, a characteristic that rarely proposed to influence treatment response, showing value in model improvement. Consistent with Zhou et al. argument that age might affect NACT efficacy to some extent in cervical cancer (34). The combined model uses a nomogram as a noninvasive and easy-to-use tool for individualized prediction of response to NACT. It simply requires conventional pretreatment CT scans and basic clinical information that readily accessible before the initiation of definitive therapy. Moreover, the features extracted from images in this study were strictly implemented according to the IBSI reference standard to improve the reproducibility.

With its convenience and high predictive accuracy of NACT response, the model could have several clinical implications. The model could facilitate the optimization of individualized LACC initial treatment planning. In fact, according to the identification of chemo-sensitivity, physicians can select chemotherapy-sensitive patients to undergo NACT followed by surgery, whereas, chemotherapy-insensitive patients are subject to radiotherapy instead. This not only benefits young patients in preserving ovarian function and sexual capacity but also improves survival outcome. Recently, NACT has been explored as an approach to fasten the initiation of treatment before access to radiation in some countries with insufficient radiation equipment (35). Radiomic model may assist in identifying qualified patients to undergo NACT before radiation. Moreover, this pre-treatment radiomic model might be of value in predicting response to chemotherapeutic response in metastatic cervical cancer, helping to determine whether alternative treatment approaches might be more appropriate. Further clinical implementation of this convenient model could potentially accelerate the progress of personalized medicine in hope of improving the prognosis of LACC patients.

The major shortcoming in our model is the lack of important clinicopathologic predictive biomarkers such as Ki-67, P53, P protein, especially pretreatment serum SCC antigen. It is largely due to heterogeneous clinical data that was retrospectively collected from different centers. Further researches incorporating pretreatment biomarkers in radiomic model development is needed. Another limitation is that the model can only be effective in squamous cell carcinoma. As various histological type of cervical cancer contributes to distinct chemosensitivity, training cohort only consisted of cases of squamous cell carcinoma for assuring data homogeneity for radiomic analysis. Other histologic subtypes of cervical cancer such as adenocarcinoma, neuroendocrine are unsuitable for model implementation. External validation from a larger dataset of multicentric LACC patients undergoing NACT ought to be applied in justifying the reproducibility and robustness of our proposed radiomic model.

Conclusion

Pretreatment CT-based radiomic combined model has the best predictive performance for NACT response with radiomic signature remain relatively stable across centers. It may serve as effective and convenient imaging-based predictors that help in patient risk stratification and improved candidate selection for neoadjuvant chemotherapy, facilitating the individualized treatment strategies tailored for LACC patients.

Data Availability Statement

All datasets for this study are included in the article/Supplementary Material.

Ethics Statement

This retrospective study was approved by institutional review board of Nanfang Hospital, Southern Medical University with written informed consent obtained from all participants. The study was conducted according to the Declaration of Helsinki.

Author Contributions

CC and JT: conception and design. XT, WL, HD, LuW, HF, ML, and PLi: collection and assembly of data. CC, PLiu, WL, ZL, CS, XT, HD, and LiW: development of methodology. CS, XT, ZL, WL, HD, LuW, LiW, PLiu, CC, and JT: data analysis and interpretation. All authors: manuscript writing and final approval of manuscript.

Funding

This paper was financially supported by the National Natural Science Foundation of China under Grant No. 81772012, the Beijing Natural Science Foundation under Grant No. 7182109, the National Key Research and Development Plan of China under Grant Nos. 2017YFA0205200, 2016YFA0100900, and 2016YFA0100902, Chinese Academy of Sciences under Grant No. GJJSTD20170004 and QYZDJ-SSW-JSC00. The study was supported by the National Science and Technology Pillar Program during the Twelfth Five-year Plan Period (Grant No. 2014BAI05B03), the Major Program of Natural Science Foundation of Guangdong Province (Grant No. 2015A030311024), the Medical Scientific Research Foundation of Guangdong Province (Grant No. A2015063).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors would like to thank department of gynaecology and obstetrics in Jiangmen Central Hospital, the First Affiliated Hospital of Zhengzhou University, Fourth Hospital of Hebei Medical University, Henan Provincial people's Hospital for their great support in the study. We also express gratitude to all the patients for generous contribution.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2020.00077/full#supplementary-material

References

1. Elit L, Fyles AW, Devries MC, Oliver TK, Fung-Kee-Fung M, Gynecology Cancer Disease Site Group. Follow-up for women after treatment for cervical cancer: a systematic review. Gynecol Oncol. (2009) 1143:528–35. doi: 10.1016/j.ygyno.2009.06.001

CrossRef Full Text | Google Scholar

2. deSouza NM, Soutter WP, Rustin G, Mahon MM, Jones B, Dina R, et al. Use of neoadjuvant chemotherapy prior to radical hysterectomy in cervical cancer: monitoring tumour shrinkage and molecular profile on magnetic resonance and assessment of 3-year outcome. Br J Cancer. (2004) 9012:2326–31. doi: 10.1038/sj.bjc.6601870

CrossRef Full Text | Google Scholar

3. Zhao H, He Y, Yang SL, Zhao Q, Wu YM. Neoadjuvant chemotherapy with radical surgery vs radical surgery alone for cervical cancer: a systematic review and meta-analysis. Onco Targets Ther. (2019) 12:1881–91. doi: 10.2147/OTT.S186451

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Kenter GG, Greggi S, Vergote I, Katsaros D, Kobierski J, Massuger L, et al., Results from neoadjuvant chemotherapy followed by surgery compared to chemoradiation for stage Ib2-IIb cervical cancer. J Clin Oncol. (2019) 37:5503. doi: 10.1200/JCO.2019.37.15_suppl.5503

CrossRef Full Text | Google Scholar

5. Benedetti-Panici P, Greggi S, Colombo A, Amoroso M, Smaniotto D, Giannarelli D, et al. Neoadjuvant chemotherapy and radical surgery versus exclusive radiotherapy in locally advanced squamous cell cervical cancer: results from the italian multicenter randomized study. J Clin Oncol. (2002) 201:179–88. doi: 10.1200/JCO.2002.20.1.179

CrossRef Full Text | Google Scholar

6. Yin M, Zhao F, Lou G, Zhang H, Sun M, Li C, et al. The long-term efficacy of neoadjuvant chemotherapy followed by radical hysterectomy compared with radical surgery alone or concurrent chemoradiotherapy on locally advanced-stage cervical cancer. Int J Gynecol Cancer. (2011) 211:92–9. doi: 10.1111/IGC.0b013e3181fe8b6e

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Gadducci A, Sartori E, Maggino T, Zola P, Cosio S, Zizioli V, et al. Pathological response on surgical samples is an independent prognostic variable for patients with Stage Ib2-IIb cervical cancer treated with neoadjuvant chemotherapy and radical hysterectomy: an Italian multicenter retrospective study (CTF Study). Gynecol Oncol. (2013) 1313:640–4. doi: 10.1016/j.ygyno.2013.09.029

CrossRef Full Text | Google Scholar

8. Huang K, Sun H, Chen Z, Li X, Wang S, Zhao X, et al. Optimal pathological response indicated better long-term outcome among patients with stage IB2 to IIB cervical cancer submitted to neoadjuvant chemotherapy. Sci Rep. (2016) 6:28278. doi: 10.1038/srep28278

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Hu T1, Li S, Chen Y, Shen J, Li X, Huang K, et al. Matched-case comparison of neoadjuvant chemotherapy in patients with FIGO stage IB1-IIB cervical cancer to establish selection criteria. Eur J Cancer. (2012) 4815:2353–60. doi: 10.1016/j.ejca.2012.03.015

CrossRef Full Text | Google Scholar

10. Ye Q, Yuan HX, Chen HL. Responsiveness of neoadjuvant chemotherapy before surgery predicts favorable prognosis for cervical cancer patients: a meta-analysis. J Cancer Res Clin Oncol. (2013) 13911:1887–98. doi: 10.1007/s00432-013-1509-y

CrossRef Full Text | Google Scholar

11. Jin L, Shen Q, Ding S, Jiang W, Jiang L, Zhu X, et al. Immunohistochemical expression of Annexin A2 and S100A proteins in patients with bulky stage IB-IIA cervical cancer treated with neoadjuvant chemotherapy. Gynecol Oncol. (2012) 1261:140–6. doi: 10.1016/j.ygyno.2012.04.005

CrossRef Full Text | Google Scholar

12. Zhu H, Chen A, Li S, Tao X, Sheng B, Chetry M, et al. Predictive role of galectin-1 and integrin alpha5beta1 in cisplatin-based neoadjuvant chemotherapy of bulky squamous cervical cancer. Biosci Rep. (2017) 37:BSR20170958. doi: 10.1042/BSR20170958

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Guo L, Zhang C, Zhu J, Yang Y, Lan J, Su G, et al. Proteomic identification of predictive tissue biomarkers of sensitive to neoadjuvant chemotherapy in squamous cervical cancer. Life Sci. (2016) 151:102–8. doi: 10.1016/j.lfs.2016.03.006

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Chung HH, Kim MK, Kim JW, Park NH, Song YS, Kang SB, et al. XRCC1 R399Q polymorphism is associated with response to platinum-based neoadjuvant chemotherapy in bulky cervical cancer. Gynecol Oncol. (2006) 1033:1031–7. doi: 10.1016/j.ygyno.2006.06.016

CrossRef Full Text | Google Scholar

15. Himoto Y, Fujimoto K, Kido A, Matsumura N, Baba T, Daido S, et al. Assessment of the early predictive power of quantitative magnetic resonance imaging parameters during neoadjuvant chemotherapy for uterine cervical cancer. Int J Gynecol Cancer. (2014) 244:751–7. doi: 10.1097/IGC.0000000000000124

CrossRef Full Text | Google Scholar

16. Prescott JW, Zhang D, Wang JZ, Mayr NA, Yuh WT, Saltz J, et al. Temporal analysis of tumor heterogeneity and volume for cervical cancer treatment outcome prediction: preliminary evaluation. J Digit Imaging. (2010) 233:342–57. doi: 10.1007/s10278-009-9179-7

CrossRef Full Text | Google Scholar

17. Liu Z, Wang S, Dong D, Wei J, Fang C, Zhou X, et al. The applications of radiomics in precision diagnosis and treatment of oncology: opportunities and challenges. Theranostics. (2019) 95:1303–22. doi: 10.7150/thno.30309

CrossRef Full Text | Google Scholar

18. Lambin P, Rios-Velazquez E, Leijenaar R, Carvalho S, van Stiphout RG, Granton P, et al. Radiomics: extracting more information from medical images using advanced feature analysis. Eur J Cancer. (2012) 484: 441–6. doi: 10.1016/j.ejca.2011.11.036

CrossRef Full Text | Google Scholar

19. Lucia F, Visvikis D, Desseroit MC, Miranda O, Malhaire JP, Robin P, et al. Prediction of outcome using pretreatment (18)F-FDG PET/CT and MRI radiomics in locally advanced cervical cancer treated with chemoradiotherapy. Eur J Nucl Med Mol Imaging. (2018) 455:768–86. doi: 10.1007/s00259-017-3898-7

CrossRef Full Text | Google Scholar

20. Chee CG, Kim YH, Lee KH, Lee YJ, Park JH, Lee HS, et al. CT texture analysis in patients with locally advanced rectal cancer treated with neoadjuvant chemoradiotherapy: a potential imaging biomarker for treatment response and prognosis. PLoS ONE. (2017) 128:e0182883. doi: 10.1371/journal.pone.0182883

CrossRef Full Text | Google Scholar

21. Li Z, Zhang D, Dai Y, Dong J, Wu L, Li Y, et al. Computed tomography-based radiomics for prediction of neoadjuvant chemotherapy outcomes in locally advanced gastric cancer: a pilot study. Chin J Cancer Res. (2018) 304:406–14. doi: 10.21147/j.issn.1000-9604.2018.04.03

CrossRef Full Text | Google Scholar

22. Bipat S, Glas AS, van der Velden J, Zwinderman AH, Bossuyt PM, Stoker J. Computed tomography and magnetic resonance imaging in staging of uterine cervical carcinoma: a systematic review. Gynecol Oncol. (2003) 911:59–66. doi: 10.1016/s0090-8258(03)00409-8

CrossRef Full Text | Google Scholar

23. Olpin J, Chuang L, Berek J, Gaffney D. Imaging and cancer of the cervix in low- and middle-income countries. Gynecol Oncol Rep. (2018) 25:115–21. doi: 10.1016/j.gore.2018.07.001

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Pecorelli S. Revised FIGO staging for carcinoma of the vulva, cervix, and endometrium. Int J Gynaecol Obstet. (2009) 1052:103–4. doi: 10.1016/j.ijgo.2009.02.012

CrossRef Full Text | Google Scholar

25. Eisenhauer EA, Therasse P, Bogaerts J, Schwartz LH, Sargent D, Ford R, et al. New response evaluation criteria in solid tumours: revised RECIST guideline (version 1.1). Eur J Cancer. (2009) 452:228–47. doi: 10.1016/j.ejca.2008.10.026

CrossRef Full Text | Google Scholar

26. Yang Z, Chen D, Zhang J, Yao D, Gao K, Wang H, et al. The efficacy and safety of neoadjuvant chemotherapy in the treatment of locally advanced cervical cancer: a randomized multicenter study. Gynecol Oncol. (2016) 1412:231–9. doi: 10.1016/j.ygyno.2015.06.027

CrossRef Full Text | Google Scholar

27. Zwanenburg A, Leger S, Vallières M, Löck S. Image biomarker standardisation initiative. ArXiv Prepr arXiv:161207003. (2019)

Google Scholar

28. Vallières M, Zwanenburg A, Badic B, Cheze Le Rest C, Visvikis D, Hatt M. Responsible radiomics research for faster clinical translation. J Nucl Med. (2018) 592:189–93. doi: 10.2967/jnumed.117.200501

CrossRef Full Text | Google Scholar

29. Fernandez-Retana J, Lasa-Gonsebatt F, Lopez-Urrutia E, Coronel-Martínez J, Cantu De Leon D, Jacobo-Herrera N, et al. Transcript profiling distinguishes complete treatment responders with locally advanced cervical cancer. Transl Oncol. (2015) 82:77–84. doi: 10.1016/j.tranon.2015.01.003

CrossRef Full Text | Google Scholar

30. Horikawa N, Baba T, Matsumura N, Murakami R, Abiko K, Hamanishi J, et al. Genomic profile predicts the efficacy of neoadjuvant chemotherapy for cervical cancer patients. BMC Cancer. (2015) 15:739. doi: 10.1186/s12885-015-1703-1

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Aerts HJ, Velazquez ER, Leijenaar RT, Parmar C, Grossmann P, Carvalho S, et al. Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach. Nat Commun. (2014) 5:4006. doi: 10.1038/ncomms5006

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Filipczuk P, Fevens T, Krzyzak A, Obuchowicz A. Glcm and glrlm based texture features for computer-aided breast cancer diagnosis. J Med Info Technol. (2012) 19:109–115.

Google Scholar

33. Karlo CA, Di Paolo PL, Chaim J, Hakimi AA, Ostrovnaya I, Russo P, et al. Radiogenomics of clear cell renal cell carcinoma: associations between CT imaging features and mutations. Radiology. (2014) 270:464–71. doi: 10.1148/radiol.13130663

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Zhou J, Li X, Huang K, Jia Y, Tang F, Sun H, et al. Young cervical cancer patients may be more responsive than older patients to neoadjuvant chemotherapy followed by radical surgery. PLoS ONE. (2016) 112:e0149534. doi: 10.1371/journal.pone.0149534

CrossRef Full Text | Google Scholar

35. Ferrandina G, Palluzzi E, Gallotta V, Gambacorta MA, Autorino R, Turco LC, et al. Neo-adjuvant platinum-based chemotherapy followed by chemoradiation and radical surgery in locally advanced cervical cancer (Lacc) patients: a phase II study. Eur J Surg Oncol. (2018) 447:1062–68. doi: 10.1016/j.ejso.2018.04.019

CrossRef Full Text | Google Scholar

Keywords: locally advanced cervical cancer (LACC), radiomics, neoadjuvant chemotherapy, response prediction, CT

Citation: Tian X, Sun C, Liu Z, Li W, Duan H, Wang L, Fan H, Li M, Li P, Wang L, Liu P, Tian J and Chen C (2020) Prediction of Response to Preoperative Neoadjuvant Chemotherapy in Locally Advanced Cervical Cancer Using Multicenter CT-Based Radiomic Analysis. Front. Oncol. 10:77. doi: 10.3389/fonc.2020.00077

Received: 19 September 2019; Accepted: 15 January 2020;
Published: 04 February 2020.

Edited by:

San-Gang Wu, First Affiliated Hospital of Xiamen University, China

Reviewed by:

Zujun Hou, Suzhou Institute of Biomedical Engineering and Technology (CAS), China
Umberto Malapelle, University of Naples Federico II, Italy

Copyright © 2020 Tian, Sun, Liu, Li, Duan, Wang, Fan, Li, Li, Wang, Liu, Tian and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Lihui Wang, wlh1984@gmail.com; Ping Liu, lpivy@126.com; Jie Tian, jie.tian@ia.ac.cn; Chunlin Chen, ccl1@smu.edu.cn

These authors share first authorship