Machine learning-based random forest predicts anastomotic leakage after anterior resection for rectal cancer
Introduction
Colorectal cancer (CRC) is the third leading cause of cancer-related deaths according to a latest statistical report (1). Although improvements have been successful at reducing the incidence of CRC (2), the morbidity and mortality still remain high. Anterior resection (AR) based on total mesorectal excision (TME) (3) is the major surgical treatment for rectal cancer.
Anastomotic leakage (AL) is one of the most common and serious complications after AR for rectal cancer. The incidence of AL varies from 3% to 21% according to previous reports (4-7). Patients with AL are likely to suffer from increased length of hospital stay and rates of re-operation, and even recurrence and mortality. Some of the patients may even require a de-functioning stoma that would significantly affect their quality of life. A temporary stoma is often recommended for patients with high risk of AL (8). However, the necessity of a temporary stoma is hard to decide upon since AL is associated with many factors, and is difficult to be predicted.
Though multiple studies have reported the predictors for AL (5,6,9), accurate prediction of AL remains difficult. According to the current knowledge, the healing of the anastomotic site depends on the tension and blood supply around the anastomotic site (10). The previous reports of risk factors included relatively few patients, and their results were not consistent (5,7,11-13). Moreover, the performance of the established models predicting the incidence of AL remains unsatisfactory. Random forest, a new and highly flexible machine learning algorithm, has wide application prospects, and has been demonstrated to have better performance in disease prediction (14). However, to the best of our knowledge, no article has reported the application of random forest in AL prediction so far. In this study, we aimed to analyze a large number of rectal cancer patients after AR to illustrate the risk factors of AL, and to create a random forest classifier to better predict the incidence of AL and give an advice on whether to a do temporary stoma. We present the following article in accordance with the TRIPOD reporting checklist (available at http://dx.doi.org/10.21037/jgo-20-436).
Methods
The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). The study was approved and monitored by the Ethics Committee of Changhai Hospital (No. CHEC2020-035). Because of the retrospective nature of the study, the requirement for informed consent was waived.
Patients
Patients who underwent AR for rectal cancer in Shanghai Changhai Hospital from August 2009 to June 2018, were included in the study. In addition, data of patients who underwent AR from July 2018 to June 2019, in our institution were collected as a group for external validation. Demographics, clinicopathological variables, and follow-up were extracted from the prospectively maintained CRC database.
Inclusion and exclusion criteria
Enrolled patients met the following inclusion criteria: rectal cancer patients; underwent AR with TME; complete clinical data. The exclusion criteria included: patients who underwent local excision, Miles or Hartmann surgery; patients with tumor >15 cm from the anal verge; patients with multiple primary colorectal carcinomas.
Diagnosis of AL
AL was defined as the defect of the intestinal wall at the anastomotic site (including suture and staple lines of neorectal reservoirs) leading to a communication between the intra- and extra-luminal compartments (15). A pelvic abscess close to the anastomosis was also considered as AL. Based on its impact on clinical management, AL was classified into three grades (grade A, B, and C) according to the International Study Group of Rectal Cancer: Grade A resulted in no change in patients’ management, whereas grade B leakage required active clinical intervention but was manageable without re-operation. Grade C required re-operation (15). In this study, AL was diagnosed via digital rectal examination (DRE), endoscopy, or imaging examination within 6 months. The follow-up was conducted via outpatient or telephone.
Variables
Demographic variables were defined and analyzed as follows: sex, age at operation, body mass index (BMI). Clinicopathological variables were diabetes, preoperative albumin (pALB) (<35 vs. ≥35 g/L), preoperative hemoglobin (pHGB) (<90 vs. ≥90 g/L), carcinoembryonic antigen (CEA) (>5 vs. ≤5 ng/mL), carbohydrate antigen 199 (CA199) (>37 vs. ≤37 U/mL), preoperative bowel stenosis or obstruction, surgical approach (laparoscopic or open), blood loss, blood transfusion, distal tumor distance from the anal verge, neoadjuvant chemoradiotherapy (nCRT), surgeon volume (high volume surgeon: amount of colorectal surgeries >100 in the previous year; a total of 8 well-trained surgeons were include, qualifications of surgeons can be seen in Table S1), temporary stoma, pathological T stage, pathological N stage, pathological types, and American Society of Anesthesiologists (ASA) score. The above variables were selected as candidate predictors because of previous reports and clinical experiences.
Statistical analysis
Statistical analyses were conducted with the statistical package for social sciences (SPSS version 22.0.0, IBM SPSS statistics, IBM Corporation, Armonk, NY, USA), R software (version 3.5.1; http://www.Rproject.org), and Python software (version 3.7.4; https://www.python.org). Descriptive statistics were computed for all variables. These included means and standard deviations (SD) for continuous factors, and frequencies for categorical factors. Comparisons of the distribution of clinicopathological characteristics were performed by using the two-tailed t-test (or Wilcoxon rank sum test as appropriate) for continuous variables and chi-square test (or the Fisher exact test as appropriate) for categorical variables. To evaluate whether temporary stoma was a risk factor for AL, propensity score matching (PSM) was implemented to reduce the possibility of selection bias by using a logistic regression model (16). P values of 0.05 or lower were considered statistically significant.
Candidate predictors incorporated in the prediction nomogram were based on the multivariate logistic regression analysis. In this study, the random forest classifier was also applied for which all the categorical data were transformed into numerical values in order to train the model. The pre-processed data set was then split into training set and validation set. Grid-search cross validation technique was used to tune the number of estimators in the classifier, and all trainings were conducted with 5-fold cross validation to prevent overfitting.
Receiver operating characteristic (ROC) curve was conducted to evaluate the clinical usefulness of the nomogram and the random forest classifier by comparing the area under the curve (AUC) formed by the real results and the predicted results (17).
Results
Demographic and clinicopathological characteristics of the patients
A total of 5,220 eligible patients were enrolled (Figure 1), including 3,454 (66.2%) patients with a temporary stoma. 1,824 patients (34.9%) were female and 3,396 patients (65.1%) were male. Patients with AL classified as Grade B accounted for 3.83% (200/5,220), while those classified as grade A and grade C accounted for 1.86% (97/5,220) and 0.56% (29/5,220), respectively. Additionally, 1,343 patients (25.7%) received nCRT and 896 patients (17.2%) underwent laparoscopic surgery. The overall incidence of AL was 6.25% (326/5,220) (Table 1).
Full table
In the test dataset with 836 patients, the incidence of clinical AL was 5.4% (45/836). Additionally, 26.6% of these patients received nCRT (222/836), which was higher than that in the training dataset. Further, 22.7% (190/836) of these patients underwent laparoscopic surgery, which was also higher than that in the training dataset. Similarly, the proportion of patients with temporary stoma in the test dataset was also higher than that in the training dataset (69.2% vs. 66.2%, respectively). Conversely, the distance of tumor from the anal verge was lower in the test set (6.53±3.28 vs. 6.99±3.18, respectively). The reasons for the differences regarding temporary stoma and nCRT between the two datasets were owing to the development of the concepts of nCRT, and techniques in laparoscopic surgery during the past few years. Patients who received neoadjuvant therapy or underwent low AR were more likely to receive a temporary stoma during surgery (Table S2).
Risk factors associated with AL
According to univariate analysis, of all the examined variables, sex, blood transfusion volume, temporary stoma, surgeon volume, nCRT, distance of tumor from the anal verge, bowel stenosis or obstruction, pHGB level, diabetes, and ASA score were associated with AL (Table 1). Moreover, we found that the incidence of AL was significantly higher during the surgeon’s first year of laparoscopic surgery compared with the following years (P=0.009) (Table 2). Therefore, the surgical approach was re-defined as open surgery, first year, and following years of laparoscopic surgery. The risk factors with P<0.2 in the univariate analysis were included in the multivariate analysis. The multivariate analysis showed that sex, surgeon volume, distance of tumor from the anal verge, bowel stenosis or obstruction, pHGB, nCRT, diabetes, and surgical approach were independent risk factors of AL (Table 3).
Full table
Full table
As shown in Table 1, patients with temporary stoma appeared to have a higher incidence of AL, which was contradictory to a previous report (18). To further confirm whether temporary stoma is associated with AL, we used PSM to balance the baseline data. According to the univariate analysis, BMI, sex, blood transfusion, blood loss, distance of tumor from the anal verge, nCRT, surgical approach, and pALB were considered as the matching variables (Table S3). A total of 2,218 patients were selected after PSM (ratio =1:1, caliper =0.05), and there was no significant difference between patients with and without temporary stoma regarding the incidence of AL (P=0.58) (Table 4).
Full table
Prediction model for AL
A random forest classifier that is a classic machine learning approach was constructed to predict the incidence of AL. All the risk factors with P<0.2 in univariate analysis were included in the random forest model, and their respective feature importance in the model were calculated. Since the number of positive class (AL patients) was much less than that of the negative class (that we only have 326 positive labels out of 5,220 patients), we up-sampled patients from positive class with replacement to solve the imbalanced class problem. After resampling, we obtained a new dataset containing 9,788 patients’ information with the same numbers of positive and negative classes.
The top 8 factors with the greatest impact weights were chosen into the final random forest model (Table 5). In order to prevent overfitting, we performed 5-fold cross-validation to train the model. During each training, 80% of the patients were sampled to form the training set, and the rest 20% patients were left as the validation set.
Full table
To evaluate our model, we used AUC-ROC curve, specificity, and sensitivity as our evaluation metrics. The average AUC for the training set was 0.89±0.00 (Figure 2A) while that for the validation set was 0.85±0.02 (Figure 2B). The sensitivity of training and validation set was 0.827 and 0.818, respectively, and the specificity was 0.739 and 0.67, respectively. According to the data in the test set, the AUC was 0.87, while the sensitivity and specificity were 0.844 and 0.697, respectively (Figure 2C).
In order to verify the efficacy of the random forest classifier, we also constructed a nomogram to predict AL based on the same clinical data (Figure 3). Patients from August, 2009, to June, 2018, were collected as training set while patients from July, 2018, to June, 2019, were regarded as test set. Based on the multivariate analysis, sex, surgeon volume, distance of tumor from the anal verge, bowel stenosis or obstruction, pHGB level, nCRT, diabetes, and surgical approach were included in the final nomogram. The AUC of training and test set was 0.735 and 0.724, respectively, which was inferior to that of the random forest classifier (Figure 4A,B).
The final random forest classifier has been stored and can be used for future applications. Moreover, we are constructing a website embedded with the model and user instructions so that future doctors can further validate and utilize our model to predict AL with our model as the reference (http://www.changhai-rc-al-prediction.org).
Discussion
In this study, 5,220 patients who underwent AR for rectal cancer were enrolled and analyzed to identify the predictors for AL. We collected 20 demographic and clinicopathological characteristics to identify the predictors for AL. After univariate and multivariate analysis, 8 predictors for AL were identified as follows: sex, surgeon volume, distance of tumor from the anal verge, bowel stenosis or obstruction, pHGB, nCRT, diabetes, and surgical approach. Furthermore, we created a random forest classifier, which is a part of machine learning, providing a more accurate prediction for the risk of AL (AUC =0.87) than the nomogram (AUC =0.724), which is based on logistic regression model. To the best of our knowledge, our study is the largest retrospective single-center study according to the data volume so far.
In our study, AUC and ROC curves confirmed that the random forest model, based on predictors per our findings, had greater predictive efficiency than the nomogram, which was widely used in the previous studies (5,6,11,12,19). However, the principle of the nomogram to predict AL is based on logistic regression that has limitations in the fitting of model creation (20). Nevertheless, machine learning that derives from the computer field has been widely used in different fields, and can partly overcome the limitations of the regression models (21,22). Random forest is a type of machine learning based on decision tree algorithm, and it has shown better predictive value than the traditional prediction model (23). However, to the best of our knowledge, no prediction model of random forest for AL has been reported to date. We constructed the random forest model to predict the incidence of AL for the first time and with the largest data volume (n=5,220). To prevent our model from overfitting, we performed cross-validation to generate robust results (24). Thus, the results in our study are likely to be more reliable and convincing. Moreover, compared to other studies, our model also showed a better predicting performance (11,12,19,20,25-27) (Table 6).
Full table
Of all enrolled patients who underwent AR for rectal cancer in this study (n=5,220), 326 (6.2%) patients were diagnosed with AL. In previous studies, there were significant variations in the incidence of AL after rectal resection, ranging from 3% to 21% (5-7). AL is more likely to occur in the left colon and rectal surgeries than in the right colon surgery (9). Clinical AL gains more attention not only because it is easy to be detected but also for its impact on follow-up treatment.
Whether a temporary stoma can reduce the incidence of AL is debatable (18). Some studies believed that a temporary stoma could indirectly accelerate anastomotic healing, thus reducing the incidence of AL (4). Other studies discovered that a temporary stoma could not reduce the incidence of AL, but could relieve the difficulty in managing postoperative complications (7). In this study, we used PSM to balance the differences of baseline data between the stoma and non-stoma groups. After PSM, we found that there was no correlation between stoma and AL, which indicated that a temporary stoma could not decrease the incidence of AL. However, our data showed that 37.0% (27/73) patients with AL needed re-operation in the non-stoma group, while in the stoma group, the re-operation rate was 0.8% (2/253) (Table S4). It was a hint for surgeons that a temporary stoma could reduce the incidence of re-operation necessitated by AL. Besides, the stoma rate was 66.2% in our study which was higher than several other studies (12,20,25). In contrast, in our study the average distance of tumor from the anal verge was lower than that in most previous studies (around 7.5 cm) (12,20,25,26). In our study, the overall average distance was 6.99±3.18 cm, and that in the stoma group was 5.75±2.52 cm, while in non-stoma group that was 9.41±2.90 cm, which implied that there were more low rectal cancer patients in our study.
Distance of tumor from the anal verge and sex are widely acknowledged as risk factors for AL in previous studies (28). Consistently, our study showed that male patients and those with lower distance of tumor from the anal verge were more likely to have AL. Similarly, preoperative bowel stenosis or obstruction, preoperative anemia (<90 g/L) or massive blood loss during operation, diabetes, and nCRT that were also considered as predictors in our study, have also been reported in several other studies (19,29). A proximal bowel stenosis or even obstruction may induce the proximal and distal bowel tissue edema, and thus increase the incidence of AL after AR. Patients with anemia or diabetes may have reduced blood supply, and tend to have a high risk of infection that may affect wound healing negatively (5). The nCRT may damage the tissue of gut lumen, and even impair the sphincter’s functioning, thus delaying the anastomotic healing. In addition to the above predictive factors, we also found that surgeon volume, and surgical approach was related to the incidence of AL. It was argued that low volume surgeons were more likely to incur AL (29,30), and we found the same conclusion. However, when it comes to surgical approach, except open surgery (12,29), performing laparoscopic surgery in first year could also raise the incidence of AL. This result indicates that we cannot neglect the human factor in high volume surgeons and laparoscopic learning curve may contribute to AL (30-32).
However, there were also some limitations in our study. This was a retrospective study in a single center, which might have caused selection bias for patients and lacked external validation of other centers. We are now conducting prospective multicenter study to further verify our model through website. Besides, other relevant variables, such as smoking, were not recorded in the database. Whether the status of KRAS is associated with AL is controversial (13). However, because of the missing data of KRAS in early years, we only collected 3,806 patients with the accurate status of KRAS and found that there was no correlation between the status of KARS and AL (P=0.784) (Table S5).
Conclusions
Our study suggests that eight factors are closely related to the incidence of AL. Our prediction model of random forest may be a practical tool for more accurate prediction of AL after AR for rectal cancer. Machine learning combined with big data has great application prospects in predicting AL. The results of our random forest model could also provide a rational advice on whether to do a temporary stoma, which might reduce the high rate of stoma and avoid the ensuing complications.
Acknowledgments
Funding: This study was supported by Clinical Science and Technology Innovation Project of Shanghai Shenkang Hospital Development Center (SHDC12016122), 234 Climbing the Discipline Program of first affiliated hospital of Naval Medical University (2019YXK032) and Ethicon Excellence in Surgery Grant (HZB-20181119-30).
Footnote
Reporting Checklist: The authors have completed the TRIPOD reporting checklist. Available at http://dx.doi.org/10.21037/jgo-20-436
Data Sharing Statement: Available at http://dx.doi.org/10.21037/jgo-20-436
Peer Review File: Available at http://dx.doi.org/10.21037/jgo-20-436
Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at http://dx.doi.org/10.21037/jgo-20-436). The authors have no conflicts of interest to declare.
Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). The study was approved and monitored by the Ethics Committee of Changhai Hospital (No. CHEC2020-035). Because of the retrospective nature of the study, the requirement for informed consent was waived.
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
References
- Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin 2020;70:7-30. [Crossref] [PubMed]
- Wieszczy P, Kaminski MF, Franczyk R, et al. Colorectal Cancer Incidence and Mortality After Removal of Adenomas During Screening Colonoscopies. Gastroenterology 2020;158:875-883.e5. [Crossref] [PubMed]
- Weitz J, Koch M, Debus J, et al. Colorectal cancer. Lancet 2005;365:153-65. [Crossref] [PubMed]
- Koyama M, Murata A, Sakamoto Y, et al. Risk Factors for Anastomotic Leakage After Intersphincteric Resection Without a Protective Defunctioning Stoma for Lower Rectal Cancer. Ann Surg Oncol 2016;23:S249-56. [Crossref] [PubMed]
- Rencuzogullari A, Benlice C, Valente M, et al. Predictors of Anastomotic Leak in Elderly Patients After Colectomy: Nomogram-Based Assessment from the American College of Surgeons National Surgical Quality Program Procedure-Targeted Cohort. Dis Colon Rectum 2017;60:527-36. [Crossref] [PubMed]
- Yao HH, Shao F, Huang Q, et al. Nomogram to predict anastomotic leakage after laparoscopic anterior resection with intracorporeal rectal transection and double-stapling technique anastomosis for rectal cancer. Hepatogastroenterology 2014;61:1257-61. [PubMed]
- Yun JA, Cho YB, Park YA, et al. Clinical manifestations and risk factors of anastomotic leakage after low anterior resection for rectal cancer. ANZ J Surg 2017;87:908-14. [Crossref] [PubMed]
- Tan WS, Tang CL, Shi L, et al. Meta-analysis of defunctioning stomas in low anterior resection for rectal cancer. Br J Surg 2009;96:462-72. [Crossref] [PubMed]
- Park JS, Huh JW, Park YA, et al. Risk Factors of Anastomotic Leakage and Long-Term Survival After Colorectal Surgery. Medicine (Baltimore) 2016;95:e2890 [Crossref] [PubMed]
- Kream J, Ludwig KA, Ridolfi TJ, et al. Achieving low anastomotic leak rates utilizing clinical perfusion assessment. Surgery 2016;160:960-7. [Crossref] [PubMed]
- Hoshino N, Hida K, Sakai Y, et al. Nomogram for predicting anastomotic leakage after low anterior resection for rectal cancer. Int J Colorectal Dis 2018;33:411-8. [Crossref] [PubMed]
- Zheng H, Wu Z, Wu Y, et al. Laparoscopic surgery may decrease the risk of clinical anastomotic leakage and a nomogram to predict anastomotic leakage after anterior resection for rectal cancer. Int J Colorectal Dis 2019;34:319-28. [Crossref] [PubMed]
- Zhang W, Lou Z, Liu Q, et al. Multicenter analysis of risk factors for anastomotic leakage after middle and low rectal cancer resection without diverting stoma: a retrospective study of 319 consecutive patients. Int J Colorectal Dis 2017;32:1431-7. [Crossref] [PubMed]
- Mogensen UB, Ishwaran H, Gerds TA. Evaluating Random Forests for Survival Analysis using Prediction Error Curves. J Stat Softw 2012;50:1-23. [Crossref] [PubMed]
- Rahbari NN, Weitz J, Hohenberger W, et al. Definition and grading of anastomotic leakage following anterior resection of the rectum: a proposal by the International Study Group of Rectal Cancer. Surgery 2010;147:339-51. [Crossref] [PubMed]
- Austin PC, Jembere N, Chiu M. Propensity score matching and complex surveys. Stat Methods Med Res 2018;27:1240-57. [Crossref] [PubMed]
- Carter JV, Pan J, Rai SN, et al. ROC-ing along: Evaluation and interpretation of receiver operating characteristic curves. Surgery 2016;159:1638-45. [Crossref] [PubMed]
- Shiomi A, Ito M, Maeda K, et al. Effects of a diverting stoma on symptomatic anastomotic leakage after low anterior resection for rectal cancer: a propensity score matching analysis of 1,014 consecutive patients. J Am Coll Surg 2015;220:186-94. [Crossref] [PubMed]
- Xiao C, Zhou M, Yang X, et al. Novel nomogram with microvascular density in the surgical margins can accurately predict the risk for anastomotic leakage after anterior resection for rectal cancer. J Surg Oncol 2019;120:1412-9. [Crossref] [PubMed]
- Klose J, Tarantino I, von Fournier A, et al. A Nomogram to Predict Anastomotic Leakage in Open Rectal Surgery-Hope or Hype? J Gastrointest Surg 2018;22:1619-30. [Crossref] [PubMed]
- Cai Z, Xu D, Zhang Q, et al. Classification of lung cancer using ensemble-based feature selection and machine learning methods. Mol Biosyst 2015;11:791-800. [Crossref] [PubMed]
- Kursa MB. Robustness of Random Forest-based gene selection methods. BMC Bioinformatics 2014;15:8. [Crossref] [PubMed]
- Wang Y, Xia ST, Tang Q, et al. A Novel Consistent Random Forest Framework: Bernoulli Random Forests. IEEE Trans Neural Netw Learn Syst 2018;29:3510-23. [Crossref] [PubMed]
- Poldrack RA, Huckins G, Varoquaux G. Establishment of Best Practices for Evidence for Prediction: A Review. JAMA Psychiatry 2020;77:534-40. [Crossref] [PubMed]
- Liu XH, Wu XR, Zhou C, et al. Conversion is a risk factor for postoperative anastomotic leak in rectal cancer patients - A retrospective cohort study. Int J Surg 2018;53:298-303. [Crossref] [PubMed]
- Ma T, Zhong Q, Cao W, et al. Clinical Anastomotic Leakage After Rectal Cancer Resection Can Be Predicted by Pelvic Anatomic Features on Preoperative MRI Scans: A Secondary Analysis of a Randomized Controlled Trial. Dis Colon Rectum 2019;62:1326-35. [Crossref] [PubMed]
- Yang SU, Park EJ, Baik SH, et al. Modified Colon Leakage Score to Predict Anastomotic Leakage in Patients Who Underwent Left-Sided Colorectal Surgery. J Clin Med 2019;8:1450. [Crossref] [PubMed]
- Shinji S, Ueda Y, Yamada T, et al. Male sex and history of ischemic heart disease are major risk factors for anastomotic leakage after laparoscopic anterior resection in patients with rectal cancer. BMC Gastroenterol 2018;18:117. [Crossref] [PubMed]
- Wang XT, Li L, Kong FB, et al. Surgical-related risk factors associated with anastomotic leakage after resection for rectal cancer: a meta-analysis. Jpn J Clin Oncol 2020;50:20-8. [Crossref] [PubMed]
- Curtis NJ, Foster JD, Miskovic D, et al. Association of Surgical Skill Assessment With Clinical Outcomes in Cancer Surgery. JAMA Surg 2020;155:590-8. [Crossref] [PubMed]
- Akiyoshi T, Kuroyanagi H, Ueno M, et al. Learning curve for standardized laparoscopic surgery for colorectal cancer under supervision: a single-center experience. Surg Endosc 2011;25:1409-14. [Crossref] [PubMed]
- Barrie J, Jayne DG, Wright J, et al. Attaining surgical competency and its implications in surgical clinical trial design: a systematic review of the learning curve in laparoscopic and robot-assisted laparoscopic colorectal cancer surgery. Ann Surg Oncol 2014;21:829-40. [Crossref] [PubMed]