Nomograms for predicting pathological response to neoadjuvant treatments in patients with rectal cancer

Dong-Lin Ren, Juan Li, Hui-Chuan Yu, Shao-Yong Peng, Wei-Da Lin, Xiao-Lin Wang, Roshan Ara Ghoorun,Yan-Xin Luo

Abstract BACKGROUND In recent decades, neoadjuvant therapy (NT) has been the standardized treatment for locally advanced rectal cancer (LARC). Approximately 8%-35% of patients with LARC who received NT were reported to have achieved a complete pathological response (pCR). If the pathological response (PR) can be accurately predicted, these patients may not need surgery. In addition, no response after NT implies that the tumor is destructive, resistant to both chemotherapy and radiotherapy, and prone to having a high metastatic potential. Therefore,developing accurate models to predict PR has great clinical significance and can help achieve individualized treatment in LARC patients.AIM To establish nomograms for predicting PR to different NT regimens based on pretreatment parameters for patients with LARC.METHODS Rectal cancer patients were identified from the database of The Sixth Affiliated Hospital, Sun Yat-sen University from January 2012 to December 2016. Logistic regression and nomograms were developed to predict the probability of pCR and good downstaging to ypT0-2N0M0 (ypTNM 0-I), respectively, based on pretreatment parameters for all LARC patients. Nomograms were also developed for three NT regimens (capecitabine/deGramont-RT, mFOLFOX6, and mFOLFOX6-RT) to predict pCR probability.RESULTS Four hundred and three patients were included in this study; 72 (17.9%) had pCR at the final pathology report, and 177 (43.9%) achieved good downstaging to ypT0-2N0M0 (ypTNM 0-I). The nomogram for predicting pCR probability showed that NT regimens, tumor differentiation, mesorectal fascia (MRF) status,and tumor length significantly influenced pCR probability. When predicting the probability of good downstaging, tumor differentiation, MRF status, and clinical T stage were the significant factors. Nomograms were developed based on NT regimens. For the capecitabine/de Gramont-RT group, the multivariate analysis showed that the neutrophil-lymphocyte ratio (NLR) was the only significant factor, thus we could not develop a nomogram for this regimen. For the mFOLFOX6-RT group, the analysis showed that the significant factors were tumor length and MRF status; and for the mFOLFOX6 group, the significant factors were tumor length and tumor differentiation.CONCLUSION We established accurate nomograms for predicting the PR to preoperative NT regimens based on pretreatment parameters for LARC patients.

Key words: Neoadjuvant therapy; Locally advanced rectal cancer; Nomogram; Prediction of pathological response; Complete pathological response; Good downstaging

INTRODUCTION

In recent decades, neoadjuvant therapy (NT) has been the standardized treatment for locally advanced rectal cancer (LARC)[1]. NT was reported to decrease the risk of local recurrence and have reduced toxicity[2,3]. Pathological complete response (pCR) is characterized as complete elimination of malignant cells in a resected specimen[4,5].Approximately 8%-35% of patients with LARC who received NT was reported to have achieved pCR[6-9]. Researchers have also found that good pathological response are associated with a longer disease-free survival (DFS) and lower local and distant recurrence rates[10-15].

Individualized treatment for LARC patients can be achieved by developing an accurate model to predict the probability of pCR or good downstaging. Some authors suggest that if pCR can be accurately predicted, these patients can be strictly followed without requiring surgery[16-18]. Radical surgery can drastically reduce the quality of life by impairing normal intestinal and genitourinary functions[19]. However, other authors argue that follow-up alone is unsafe and that the pathology cannot be accurately assessed without surgery after NT[20]. In addition, no tumor response or progression after NT implies that the tumor is destructive, resistant to both chemotherapy and radiotherapy, and prone to having a high metastatic potential.Thus, identifying potential responders and non-responders may aid in predicting treatment outcomes and choices.

Previous studies have reported that low carcinoembryonic antigen (CEA)levels[21,22], high pretreatment hemoglobin (HB) levels, early clinical T stage (cT)[23],early clinical N stage (cN), small tumor size, and long radiation surgery interval[6,23-25]are related to pCR probability. However, few models or nomograms have been established and even fewer are used clinically to predict a good pathological response after NT for LARC. Additionally, few models are available to predict neoadjuvant treatments. Therefore, developing accurate models to predict pathological responses has great clinical significance and remains a great challenge.

In this study, by analyzing pretreatment parameters in LARC patients before NT at our institution, we established accurate models and nomograms to predict the probability of pCR and good downstaging, respectively, with currently available pretreatment parameters that can be easily used in clinical decision-making.

MATERIALS AND METHODS

Patients

Rectal cancer patients were identified from the database of The Sixth Affiliated Hospital, Sun Yat-sen University from January 2012 to December 2016. Four hundred and three patients who met the following criteria were included: histopathologically confirmed rectal adenocarcinoma, ＞18 years old, tumor located no more than 12 cm above the anal verge, clinical stage of cT3/4 or lymph node (+), and non-metastatic.All patients were assessed via abdominal-pelvic computed tomography (CT) and pelvic magnetic resonance imaging (MRI), and 44 (11%) patients received transrectal ultrasound testing. All received NT followed by total mesorectal excision (TME)radical surgery.

We collected all available clinical information before treatment: gender, age, body mass index (BMI), cT, cN, mesorectal fascia (MRF) status, tumor differentiation, tumor length (TL), distance of tumor from the anal verge (DTAV), tumor circumferential extent (TCE), serum tumor marker CEA, HB, neutrophil-lymphocyte ratio (NLR),platelets (PLT), apolipoprotein A-1 (ApoA1), apolipoprotein B (ApoB), and NT regimen type. All tumor-related parameters such as cT, cN, MRF status, DTAV, and TCE were assessed by MRI. Tumor length was also assessed by using MRI to measure the maximum diameter of tumor. CT, transrectal ultrasound, and endoscopy provided additional verification. Tumor differentiation was identified by enteroscopic pathology.

This retrospective study was approved by the Institutional Review Board of The Sixth Affiliated Hospital, Sun Yat-sen University.

Therapy

During the period we identified patients for the current study, a clinical trial(FOWARC) was conducted at our institution comparing the effectiveness and safety of administering only chemotherapy with mFOLFOX6 or mFOLFOX6 plus radiotherapy to LARC patients with the effectiveness and safety in patients undergoing a standard NT regimen with fluorouracil plus radiotherapy.Consequently, 273 (67.7%) patients in our study were included in the FOWARC trial.The NT regimens included in our study were capecitabine /fluorouracil plus radiotherapy (standard group, capecitabine/deGramont-RT), mFOLFOX6 without radiotherapy (mFOLFOX6), and mFOLFOX6 plus radiotherapy (mFOLFOX6-RT).Details of all these treatments have been reported in previous studies[26,27]. The radiation dose for the radiotherapy was 46.0-50.4 Gy, delivered as 1.8-2.0 Gy/d, and the dose was the same in the capecitabine/deGramont-RT and mFOLFOX6-RT groups. Patients in the capecitabine/deGramont-RT and mFOLFOX6-RT groups underwent standard TME radical surgery after NT. The interval between radiation and surgery was 6-12 wk in mFOLFOX6-RT and de Gramont-RT groups. The interval between chemotherapy and TME radical surgery was about 2-4 wk in the mFOLFOX6 group.

Pathological assessment

All resected specimens were examined to determine the post-TN staging according to the American Joint Committee on Cancer-International Union Against Cancer(seventh edition), which is currently considered the most accurate and standard staging system in this period[28]. pCR was defined as no malignant cells found in the resected specimens, including the primary tumor and lymph nodes, and ypT0-2N0M0(ypTNM 0-I) was classified as good downstaging.

Statistical analysis

Chi-square analysis was selected for the univariate logistic regression analysis of the countable data. Normal distribution tests were performed for the metrological data,nonparametric test was used for the indicators that were not normally distributed,and the expression form of the median (upper quartile to lower quartile) is used.

Parameters such as age (≤ 60 years vs ＞ 60 years), BMI (＜ 25 kg/cm2vs ≥ 25 kg/cm2),CEA (＞ 5 ng/mL vs ≤ 5 ng/mL), HB (≤ 125 g/L vs ＞ 125 g/L), NLR (＞ 3 vs ≤ 3), DTAV(＜ 5 cm vs ≥ 5 cm) and TL (＞ 3 cm vs ≤ 3 cm) were dichotomized according to previous studies[24,29,30]. PLT, ApoA1, ApoB, and the interval were used as continuous variables,however, all these variables were not normally distributed, so a nonparametric test was used.

Univariate logistic regression analysis was used to analyze variables related to the probability of pCR or good downstaging. Variables that achieved significance at P ≤0.05 in the univariate logistic regression analysis were further analyzed into the forward stepwise multivariable logistic regression. Multivariate logistic regression analysis was used to construct nomograms. Because the NT regimen was a statistically significant factor for predicting pCR probability, all patients were divided into three subgroups (the capecitabine/deGramont-RT, mFOLFOX6-RT, or mFOLFOX6 groups) based on the NT regimens. We then attempted to develop three nomograms based on the different NT regimens to predict pCR probability. The C-index was acquired for the nomogram, and internal validation using the bootstrap method was performed to determine the adjusted C-index. Calibration curves of the nomograms were generated to show the relationship between the predicted and observed outcomes. All statistical analyses were performed using SPSS 24.0 and R 3.5.1.

RESULTS

Of the 403 patients in our study, 281 (69%) were men. As assessed pathologically, 72(17.86%) individuals achieved pCR; 177 (43.9%) patients achieved ypTNM 0-I and were classified as having good downstaging.

The interval between radiation and surgery was 52 (47-59) d in the mFOLFOX6-RT group and 54 (49-58.25) d in the deGramont-RT group, and there was no significance difference between the two groups. The interval between chemotherapy and surgery was 22 (18-25.75) d in the mFOLFOX6 group, which was much shorter than those in the other two groups.

All patients received TME surgery (28 underwent APR and 375 underwent sphincter-saving surgery).

pCR patients and non-pCR patients did not differ significantly in terms of gender,BMI, CEA, NLR, HB, PLT, ApoA1, ApoB, cT, cN, or TCE in the univariate analysis (P＞ 0.05); however, significant differences were found for age, tumor differentiation, TL,DTAV, MRF status, interval, and NT regimen (Table 1). Statistically significant factors in the univariate logistic regression analysis (P ≤ 0.05) to predict pCR were entered into a multivariate analysis, in which NT regimen types (aP ＜ 0.05), tumor differentiation (bP ＜ 0.05), TL (cP ＜ 0.05), and MRF status (dP ＜ 0.05) were significantly associated with pCR probability (Table 2). For the NT regimens, the odds ratio (OR)was 5.339 [95% confidence interval (CI): 2.394-11.903] for the mFOLFOX6-RT regimen compared with the capecitabine/deGramont-RT regimen. The mFOLFOX6 regimen and capecitabine/deGramont-RT regimen did not differ significantly. For tumor differentiation, the OR was 2.966 (95%CI: 1.449-6.069) for well tumor differentiation compared with moderate-poor differentiation. For TL (＞ 3 cm) compared with TL (≤ 3 cm), the OR was 2.608 (95%CI: 1.347-5.052), and for MRF(-) compared with MRF(+),the OR was 2.729 (95%CI: 1.199-6.211).

Patients with good downstaging and bad downstaging did not significantly differ in terms of age, gender, BMI, NLR, HB, PLT, ApoA1, ApoB, cN, TCE, the interval, or NT regimen in the univariate analysis (P ＞ 0.05); however, significant differences were found for CEA, tumor differentiation, DTAV, TL, cT, and MRF status in the univariate logistic regression analysis for good downstaging (Table 3). In the multivariate analysis, tumor differentiation (eP ＜ 0.05), MRF statuses (fP ＜ 0.05), and cT (gP ＜ 0.05)were significantly associated with the probability of good downstaging (Table 4). The OR was 4.814 (95%CI: 2.343-9.892) for well differentiation compared with moderatepoor differentiation, 4.226 (95%CI: 1.894-9.426) for MRF(-) compared with MRF(+),and 0.248 (95%CI: 0.063-0.974) for cT3 compared with cT2.

Table 1 Predictive factors for complete pathological response in univariate logistic regression analysis for all patients

Because the type of NT regimen was a statistically significant factor for predicting pCR probability, patients were divided into three subgroups (the capecitabine/deGramont-RT, mFOLFOX6, and mFOLFOX6-RT groups) based on the NT regimen. Table 5 shows the distribution of pretreatment clinical parameters in the NT regimen groups. No differences were found in any factors between the three groups except age (hP ＜ 0.05) and DTAV (iP ＜ 0.05). In the univariate analysis of the capecitabine/deGramont-RT group, NLR was the only significant factor for predicting pCR probability (Table 6). NLR (＞ 3) (jP ＜ 0.05) was the only significant factor, with an OR of 4.278 (95%CI: 1.051-17.413) compared with NLR ≤ 3 in the further multivariate analysis (Table 7). We could not develop a nomogram to predict pCR probability in this case.

Table 8 shows that TL and MRF status were significant factors predicting pCRprobability in the univariate analysis of the mFOLFOX6-RT regimen. TL (kP ＜ 0.05)and MRF(+) (lP ＜ 0.05) were significant factors, with an OR of 2.452 (95%CI: 1.015-5.926) for TL(≤ 3 cm) compared with TL(＞ 3 cm) and an OR of 3.829 (95%CI: 1.42-10.325) for MRF(-) compared with MRF(+) in the further multivariate analysis (Table 9).

Table 2 Predictive factors for complete pathological response in multivariate logistic regression analysis for all patients

In the univariate analysis of the mFOLFOX6 regimen, tumor differentiation and TL were significant factors for predicting pCR probability (Table 10). Further multivariate analysis showed that differentiation (mP ＜ 0.05) and TL (nP ＜ 0.05) were significant factors, with an OR of 8.881 (95%CI: 2.263-34.85) for well tumor differentiation compared with moderate-poor differentiation and an OR of 4.805 (95%CI: 1.25-18.466)for TL (≤ 3 cm) compared with TL (＞ 3 cm) (Table 11).

Predictive nomograms established for pCR and good downstaging

Nomograms were developed based on the significant factors in the multivariate logistic regression analysis. The nomogram for predicting pCR probability showed that NT regimen and tumor differentiation influenced the probability of pCR,followed by TL and MRF status (Figure 1). When developing the nomogram to predict the probability of good downstaging, tumor differentiation and MRF status were the most important, followed by cT (Figure 2). We attempted to develop three nomograms to predict pCR probability based on NT regimens, because only one significant factor was found for the capecitabine/deGramont-RT regimen, we could not develop a nomogram. MRF status and TL were the significant factors for the mFOLFOX6-RT group (Figure 3). For the mFOLFOX6 group, tumor differentiation and TL were the significant factors in the nomogram for predicting pCR probability(Figure 4). Using the nomograms, we could easily calculate the probability of pCR and ypTNM (0-I), and we calculated pCR probabilities based on the NT regimens.

We used 1000 bootstrap resamples to compute an adjusted C-index, which was 79.34% for predicting pCR (95%CI: 73.48%-85.21%) for all patients, with a C-index of 69.85% (95%CI: 60.94%-78.76%) for the mFOLFOX6-RT group and 83.39% (95%CI:67.26%-93.52%) for the mFOLFOX6 group. For predicting good downstaging, the adjusted C-index was 68.08% (95%CI: 63.08%-73.07%) for all patients. Calibration curves between predicted and actual observations by internal validation demonstrated that these nomograms showed good statistical performance for predicting the probability of pCR and good downstaging. Figures 5-8 show the calibration curve between the predicted and actual observations by internal validation and demonstrates that these nomograms showed good statistical performance for predicting the probability of pCR and good downstaging.

Table 3 Predictive factors for good downstaging in univariate logistic regression analysis for all patients

DISCUSSION

At present, preoperative NT is the standard treatment for patients with LARC.Patients who respond well to preoperative treatment have shown to have an excellent long-term prognosis. Knowledge of these factors ultimately leads to individualized treatment strategies; for example, patients who do not respond to the usual management can choose an aggressive preoperative regimen before NT. Conversely,to accurately determine an excellent pathological response after NT, surgeons may choose to perform local excision or a “watch and wait” strategy. In some cases, radical surgical resection may not benefit for some patients who achieve a good response because radical surgical resection may be associated with high rates of temporary or permanent stomas, defecatory disorders, urinary and sexual dysfunction, and unnecessary mortality[31,32]. pCR after NT is reported to have an excellent long-term prognosis irrespective of the treatment strategy, so noninvasive treatment strategies,such as a “watch and wait” strategy, have become more popular for patients whoachieve a good response[33,34]. Thus, understanding the factors that predict the pathological response to NT is becoming crucial.

Table 4 Predictive factors for good downstaging in multivariate logistic regression analysis for all patients

Our study identified clinical variables related to the pCR and good downstaging of LARC patients after NT. In the nomogram, we demonstrated that type of NT regimen,tumor differentiation, MRF status, and TL predicted pCR, whereas tumor differentiation, MRF status, and cT predicted good downstaging.

In our model, the mFOLFOX6-RT group had a higher probability of pCR compared with the capecitabine/de Gramont-RT group. We acknowledge that a potential selection bias may contribute to this high pCR rate. The data missing were more frequently in patients not reaching pCR than those with pCR, possibly resulting from more attentions pCR-patients got in clinical practice, follow-up, or research work. The pCR rate was 35.7% for mFOLFOX6-RT, which is higher than that of FOWARC[26,27]. It is expected since this is a single-center statistic result, while the FORWARC trial is a muli-center study. Although the benefits of oxaliplatin have not been demonstrated and it is not part of standard NT regimens, oxaliplatin is a standard component of chemotherapy for treating colon cancer[35]. Importantly, it has been reported in more and more studies[36,37]that the regimen combining mFOLFOX 6 with RT is getting a higher pCR rate of 38% in a clinical trial published in Lancet Oncology[38]. However,the role of oxaliplatin adding to fluorouracil-based neoadjuvant chemoradiotherapy(CRT) is unclear for LARC patients, and more studies are needed in the future.

pCR probabilities did not significantly differ between the capecitabine/deGramont-RT and themFOLFOX6 groups. Additionally, NT regimen was not a significant factor for predicting the probability of good downstaging. To avoid radiotherapeutic harm to LARC patients, the use of neoadjuvant chemotherapy alone has been proposed.Our model showed that patients treated with the mFOLFOX6 regimen alone had an acceptable probability of pCR and good downstaging. Thus, some chemosensitive patients can avoid radiation therapy.

Tumor differentiation was associated with both pCR and good downstaging. Well differentiation was associated with a higher pCR probability, which is consistent with previous studies[23,39], and was related to good downstaging compared with moderatepoor differentiation. Patients with a well differentiated tumor have a higher pCR probability indicating that a mild NT regimen, local resection, or “watch and wait”strategy can be considered after NT. Patients with moderate-poor differentiation may have a poor likelihood of pCR and good downstaging indicating that a “watch and wait” strategy requires careful selection.

Factors associated with pCR and good downstaging both included MRF status.MRF(+) implies that the tumor is aggressive, and even after NT, patients with MRF(+)may have a poor likelihood of pCR and good downstaging, indicating that an enhanced NT regimen and radical surgery are needed and a “watch and wait”strategy requires careful selection. While patients with MRF(-) may have a higher pCR probability indicating that a mild NT regimen, local resection, or “watch and wait”strategy can be considered after NT.

Table 5 Distribution of pretreatment clinical parameters in different neoadjuvant therapy regimen groups

TL was also a significant factor in the multivariate logistic regression analysis for predicting pCR probability in all patients. Van Stiphout et al[40]reported that TL was related to the probability of pCR after NT, although this study was based on data from positron emission tomography (PET)-CT results. TL (＞ 3 cm) implies an aggressive tumor, and even after NT, patients with TL (＞ 3 cm) may have a lower pCR probability indicating that an enhanced NT regimen and radical surgery are needed.While for patients with TL (≤ 3 cm) may have a higher pCR probability indicating that a mild NT regimen, local resection, or “watch and wait” strategy can be considered after NT.

For predicting the probability of good downstaging, cT was also a significant factor in the multivariate logistic regression analysis. In a study performed by Joye et al[41], a low cT stage was linked with ypT0-1N0, and it, together with other factors, could beused as a selection tool for organ-preserving strategies. Our study also showed that low cT stage was more likely to achieve good downstaging with NT, and indicated that less invasive surgery can be selected.

Table 6 Predictive factors for complete pathological response in univariate logistic regression analysis for the capecitabine/deGramont-RT regimen

For the capecitabine/de Gramont-RT regimen, the only significant factor was the NLR. Kim et al[42]showed that an elevated NLR before CRT can be used to predict poor tumor response and adverse prognostic factors. As lymphocytes decrease and neutrophils increase, NLR affects the adverse tumor reaction and adverse prognosis.Our study showed that the NLR before NT was related to better pathological responses to the capecitabine/deGramont-RT regimen; thus, further studies are needed to validate the relationship between NLR and pathological response to NT.

Do you know why? In fact, she felt very sad when she read the letter for the first time. But she read it for a few more times and , at last, she found the key – only every other line should be read, that is the first line, the third, the fifth … and so on to the end.

For the mFOLFOX6-RT regimen, the significant factors for predicting pCR probability were MRF status and TL. MRF(+) and long TL indicated that the tumor was aggressive and patients have heavy tumor load, and were related to poor neoadjuvant pathological responses. Patients with moderate-poor differentiation and TL (＞ 3 cm) have a lower pCR probability indicating that the efficacy of CRT is poor for these patients, and radical surgery can be directly selected without NT to avoid complications caused by CRT.

For the mFOLFOX6 regimen, the nomogram for predicting pCR probability showed that differentiation and TL were significant factors. Poor differentiation andlong TL indicated an aggressive tumor, and they were related to a poor neoadjuvant pathological response. Patients with moderate-poor differentiation and TL (＞ 3 cm)will have a lower probability of pCR indicating that radical surgery after NT is needed, or mFOLFOX6-RT regimen is chosen to increase pCR probability. However,good differentiation and short TL were related to a good neoadjuvant pathological response and high probability of pCR indicating that local resection or a “watch and wait” strategy can be chosen.

Table 7 Predictive factors for complete pathological response in multivariate logistic regression analysis for the capecitabine/deGramont-RT regimen

To the best of our knowledge, our study is the first to use different NT regimen types to predict a pathological response. We established an accurate model with easily obtained variables to predict the probability of pCR and good downstaging. Our analysis was also strengthened through cross-validation. These models can be used to assist with individualized therapy as follows. For LARC patients expected to have a poor pathological response, NT and NT-related harm can be avoided. For patients expected to have good pathological responses to chemotherapy alone, radiotherapy can be avoided. For patients who are not expected to have good pathological response from a standard NT regimen, an enhanced mFOFOLX6-RT regimen can be considered. For patients who are not expected to have good pathological response from an enhanced regimen, radical surgery can be directly chosen without NT to avoid complications caused by CRT. For patients with a high probability of pCR after NT, local resection or a “watch and wait” strategy can be used to avoid complications.

Our analysis had several limitations. First, this was a retrospective study, in which some factors associated with pCR were unavailable, such as smoking status,molecular subtype and so on. Second, mFOLFOX6 and mFOLFOX6-RT are not the standard regimens for LARC, and both regimens remain in the clinical trial phase.Finally, our nomograms are based on the experience of our single institution. These results must be validated in a group of independent external institutions.

The nomograms established in our study can be used to evaluate the probability of a pathological responses before NT and after NT. However, additional studies are required to answer clinical questions, regarding which patients can be treated only with neoadjuvant chemotherapy, which patients need oxaliplatin added to the neoadjuvant CRT, which patients need radical surgery, which patients can undergo local excision, and which patients can be managed with a “watch and wait” strategy after achieving a good response.

We established accurate nomograms to predicting the pathological responses to different preoperative NT regimens based on pretreatment parameters for LARC patients. These nomograms can be used to distinguish patient types and facilitate developing individualized treatments.

Table 8 Predictive factors for complete pathological response in univariate logistic regression analysis for the mFOLFOX6-RT regimen

Table 9 Predictive factors for complete pathological response in multivariate logistic regression analysis for the mFOLFOX6-RT regimen

Table 10 Predictive factors for complete pathological response in univariate logistic regression analysis for the mFOLFOX6 regimen

Platelet, apolipoprotein A-1, apolipoprotein B, and the interval were calculated as metrological data, and others are counting data. pCR: Complete pathological response; NLR: Neutrophil-lymphocyte ratio; DTAV: Distance of tumor from the anal verge; TL: Tumor length; TCE: Tumor circumferential extent; PLT: Platelet; ApoA1: Apolipoprotein A-1; ApoB: Apolipoprotein B; MRF: Mesorectal fascia; CEA: Carcinoembryonic antigen.

Table 11 Predictive factors for complete pathological response in multivariate logistic regression analysis for the mFOLFOX6 regimen

Figure 1 Nomogram for predicting the probability of pathological complete response for all patients. MRF: Mesorectal fascia.

Figure 2 Nomogram for predicting the probability of good downstaging (ypTNM stage 0-I) for all patients. MRF: Mesorectal fascia.

Figure 3 Nomogram for predicting the probability of pathological complete response for the mFOLFOX6-RT regimen. MRF: Mesorectal fascia.

Figure 4 Nomogram for predicting the probability of pathological complete response for the mFOLFOX6 regimen.

Figure 5 Calibration curve of the predicted and observed probabilities of pathological complete response for all patients.

Figure 6 Calibration curve of the predicted and observed probabilities of good downstaging for all patients.

Figure 7 Calibration curve of the predicted and observed probabilities of pathological complete response for the mFOLFOX6-RT regimen.

Figure 8 Calibration curve of the predicted and observed probabilities of pathological complete response for the mFOLFOX6 regimen.

ARTICLE HIGHLIGHTS

Research background

In recent decades, neoadjuvant therapy (NT) has been the standardized treatment for locally advanced rectal cancer (LARC). Approximately 8-35% of patients with LARC who received NT were reported to have achieved a complete pathological response (pCR). If the pathological response can be accurately predicted, these patients may not need surgery. In addition, no response after NT implies that the tumor is destructive, resistant to both chemotherapy and radiotherapy, and prone to having a high metastatic potential.

Few models or nomograms have been established and even fewer are used clinically to predict a good pathological response after NT for LARC. Therefore, developing accurate models to predict pathological response (PR) has great clinical significance and can help achieve individualized treatment in LARC patients.

Research motivation

Our goal was to establish nomograms that can be used to assist with individualized therapy as follows: for which patients NT and NT-related harm can be avoided; which patients will have good pathological responses to chemotherapy alone and radiotherapy can be avoided; which patients will have a good pathological response from a standard NT regimen, which patients need an enhanced mFOFOLX6-RT regimen; and which patients can use local resection or a“watch and wait” strategy to avoid complications. Solving these problems may aid in clinical treatment choices.

Research objectives

Our main objective was to establish nomograms for predicting a pathological response to different NT regimens based on pretreatment parameters for patients with LARC. We established accurate nomograms for predicting the pathological response to preoperative NT regimens based on pretreatment parameters for LARC patients. These nomograms can be used to distinguish patient types and facilitate developing individualized treatments.

Research methods

Rectal cancer patients were identified from the database of The Sixth Affiliated Hospital, Sun Yat-sen University from January 2012 to December 2016. Four hundred and three patients who met the criteria were included. We collected all available clinical information before treatment.

The NT regimens included in our study were capecitabine/fluorouracil plus radiotherapy(standard group, capecitabine/deGramont-RT), mFOLFOX6 without radiotherapy(mFOLFOX6), and mFOLFOX6 plus radiotherapy (mFOLFOX6-RT). The radiation dose for the radiotherapy was 46.0-50.4 Gy, delivered as 1.8-2.0 Gy/d.

pCR was defined as no malignant cells found in the resected specimens, including the primary tumor and lymph nodes, and ypT0-2N0M0 (ypTNM 0-I) was classified as good downstaging.

Univariate logistic regression analysis was used to analyze variables related to the probability of pCR or good downstaging. Variables that achieved significance at P ≤ 0.05 in the univariate logistic regression analysis were further analyzed into the forward stepwise multivariable logistic regression. Multivariate logistic regression analysis was used to construct the nomograms. Because the NT regimen was a statistically significant factor for predicting pCR probability, we then attempted to develop three nomograms based on the different NT regimens to predict pCR probability. The C-index was acquired for the nomogram, and internally validated using the bootstrap method to determine the adjusted C-index. Calibration curves of the nomograms were generated to show the relationship between the predicted and observed outcomes.

All statistical analyses were performed using SPSS 24.0 and R 3.5.1.

Research results

Of the 403 patients in our study, 281 (69%) were men. As assessed pathologically, 72 (17.86%)individuals achieved a pCR; 177 (43.9%) patients achieved ypTNM 0-I and were classified as having good downstaging.

Significant differences were found for age, tumor differentiation, TL, DTAV, mesorectal fascia(MRF) status, interval, and NT regimen in the univariate analysis. In the multivariate analysis,NT regimen types, tumor differentiation, TL, and MRF status were significantly associated with pCR probability.

Significant differences were found for carcinoembryonic antigen (CEA), tumor differentiation,distance of tumor from the anal verge (DTAV), tumor length (TL), cT, and MRF status in the univariate logistic regression analysis for good downstaging. In the multivariate analysis, tumor differentiation, MRF statuses, and cT were significantly associated with the probability of good downstaging.

Table 5 shows the distribution of pretreatment clinical parameters in the NT regimen groups.No differences were found in any factors between the three groups except age and DTAV.

In the univariate analysis of the capecitabine/deGramont-RT group, NLR was the only significant factor for predicting pCR probability. NLR (＞ 3) was the only significant factor compared with NLR ≤ 3 in the further multivariate analysis. We could not develop a nomogram to predict pCR probability in this case.

In the univariate analysis of the mFOLFOX6-RT regimen, TL and MRF status were significant factors predicting pCR probability. TL and MRF(+) were significant factors in multivariate analysis.

In the univariate analysis of the mFOLFOX6 regimen, tumor differentiation and TL were significant factors for predicting pCR probability. Further multivariate analysis showed that differentiation and TL were significant factors.

Nomograms were developed based on the significant factors in the multivariate logistic regression analysis. We used 1000 bootstrap resamples to compute an adjusted C-index.Calibration curves between predicted and actual observations by internal validation demonstrated that these nomograms showed good statistical performance for predicting the probability of pCR and good downstaging.

Research conclusions

We established accurate nomograms to predicting the pathological responses to different preoperative NT regimens based on pretreatment parameters for LARC patients. These nomograms can be used to distinguish patient types and facilitate developing individualized treatments.

To the best of our knowledge, our study is the first to use different NT regimen types to predict a pathological response. We established an accurate model with easily obtained variables to predict the probability of pCR and good downstaging. Our analysis was also strengthened through cross-validation. These models can be used to assist with individualized therapy as follows. For LARC patients expected to have a poor pathological response, NT and NT-related harm can be avoided. For patients expected to have good pathological responses to chemotherapy alone, radiotherapy can be avoided. For patients who are not expected to have good pathological response from a standard NT regimen, an enhanced mFOFOLX6-RT regimen can be considered. For patients with a high probability of pCR after NT, local resection or a“watch and wait” strategy can be used to avoid complications.

Our analysis had several limitations. First, this was a retrospective study, in which some factors associated with pCR were unavailable, such as smoking status, molecular subtypes and so on. Second, mFOLFOX6 and mFOLFOX6-RT are not the standard regimens for LARC, and both regimens remain in the clinical trial phase. Finally, our nomograms are based on the experience of our single institution. These results must be validated in a group of independent external institutions.

The nomograms established in our study can be used to evaluate the probability of a pathological responses before NT and after NT. However, additional studies are required to answer clinical questions regarding which patients can be treated only with neoadjuvant chemotherapy, which patients need oxaliplatin added to the neoadjuvant chemoradiotherapy,which patients need radical surgery, which patients can undergo local excision, and which patients can be managed with a “watch and wait” strategy after achieving a good response.

Research perspectives

In the future, we plan to include a larger number of patients to enhance the accuracy of the prediction. On the other hand, we plan to add a second external cohort for validation to strengthen the reliability of the nomogram.