Development of a nomogram predicting perineural invasion risk and assessment of the prognostic value of perineural invasion in colon cancer: a population study based on the Surveillance, Epidemiology, and End Results database
Highlight box
Key findings
• Perineural invasion (PNI) is a survival predictor for colon cancer (CC) patients. Our nomogram has clinical application value for predicting PNI risk in CC patients. Pathologic stage T4 (pT4) CC patients without PNI might benefit from combined chemotherapy and radiotherapy (CCRP).
What is known and what is new?
• PNI in CC is widely associated with poor prognosis.
• Our study develops a nomogram for predicting the risk of PNI in patients with CC before operation. Our study also demonstrates that pT4 patients without PNI might benefit from CCRP.
What is the implication, and what should change now?
• Our predictive model can provide a certain clinical application value for predicting PNI risk in patients with CC before operation.
Introduction
Colon cancer (CC) is one of the most common malignant cancers, with 1.15 million new cases reported every year and 576,858 deaths related to it (1,2). Currently, the American Joint Committee on Cancer (AJCC) tumor-node-metastasis (TNM) staging system is considered the clinical treatment guide, and TNM stage is the primary prognostic factor for CC. However, previous studies have reported that various histopathological factors, such as perineural invasion (PNI) and vascular invasion, are associated with a poor prognosis in CC patients, indicating that TNM staging is not accurate enough to predict the prognosis of all CC patients (3-5).
PNI, also called neurotropic carcinomatous spread, is defined as invasion into the neural tissue and/or nerve sheaths, which may prove to be another metastatic route (6). Growing through nerves has become a common mode of metastasis in CC, in addition to direct infiltration and dissemination through blood and lymph channels (7). In colorectal cancer (CRC), PNI is widely associated with poor outcome (8,9). According to the National Comprehensive Cancer Network (NCCN) clinical practice guidelines, PNI is also defined as a high-risk feature (10). However, despite its strong clinical importance, PNI can only be diagnosed by histopathology after surgery. Preoperative prediction of PNI in CC patients is very useful for the determination of treatment plans. Therefore, further studies are urgently needed to construct a prediction model for PNI.
At present, radical surgery combined with adjuvant chemotherapy (ChemT) is the cornerstone of treatment for locally advanced CC. However, the prognosis remains dismal (11,12). Despite radical surgery and adjuvant ChemT, nearly 30% of CC patients will have recurrence, and the prognosis remains poor (13,14). While adjuvant radiotherapy (RT) can only be considered in several clinical scenarios, such as positive margins or advanced local disease [pathologic stage T4 (pT4)] (15,16), combined chemotherapy and radiotherapy (CCRP) has been less explored in CC patients. And whether there is an influence of PNI on the outcomes of CCRP has not been reported. Thus, additional studies are needed to confirm the efficacy of CCRP in pT4 CC patients with or without PNI.
In the present study, we developed a nomogram for predicting PNI in CC patients using the population-based Surveillance, Epidemiology, and End Results (SEER) database. The effect of PNI on CC prognosis was further investigated using Kaplan-Meier analysis, a competing risk model, and the Fine-Gray multivariable regression model. Finally, we investigated the efficiency of CCRP in pT4 CC patients with or without PNI. We present this article in accordance with the TRIPOD reporting checklist (available at https://tcr.amegroups.com/article/view/10.21037/tcr-24-1030/rc).
Methods
Patient cohort
The SEER Program of the National Cancer Institute (NCI) is a cancer database based on the United States (US) population that consists of clinical data of cancer patients from 18 registries and covers approximately 26% of the US population (17). The data of CC patients in the SEER database from January 2010 to December 2015 were queried using SEER*Stat software version 8.3.9 (www.seer.cancer.gov). All procedures were performed according to approved guidelines. As SEER is an anonymized public database, there was no ethical approval needed. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013).
The inclusion criteria included the following: (I) primary CC (histology codes: 8140–8147, 8210, 8211, 8220, 8221, 8260–8263, 8480, 8481, 8490; site codes: C18.0 and C18.2–18.9); (II) aged ≥18 years; (III) complete follow-up; and (IV) CC as the only primary malignancy. In our study, the following variables were included: race, age at diagnosis, sex, tumor size, tumor grade, tumor site, histology, T status according to AJCC 7th edition, N status according to AJCC 7th edition, carcinoembryonic antigen (CEA), PNI, metastasis sites (liver, lung, brain, bone), and information on therapy (surgery status, RT status, ChemT status), vital status, survival months, cause of death, and marital status.
Following preliminary selection, patients were excluded if (I) they had a second primary tumor; (II) they had an overall survival time <1 month; or (III) they had incomplete patient records.
A total of 51,826 CC patients were included. To ensure the accuracy of the nomogram, all included CC patients were randomly divided into a development cohort (70%) and an internal validation cohort (30%). To evaluate the effect of PNI on prognosis, all patients were also divided into two groups by PNI: the PNI (−) group and PNI (+) group. The tumor sites were grouped as follows: right-side colon cancer (RCC; cecum, ascending colon, hepatic flexure and transverse colon) and left-side colon cancer (LCC; splenic flexure, descending colon and sigmoid colon). “No/unknown” RT records were considered as no RT. “No/unknown” ChemT records were considered as no ChemT.
End points
Final follow-up was conducted in November 2020, with a median follow-up of 41 months (ranging from 1 to 95 months). Overall survival was defined as the time from diagnosis to death. CC-specific death and CC-specific survival were defined as the time from the initiation of therapy to death or survival from CC.
Statistical analysis
The statistical analysis in this study was performed in SPSS 24.0 (IBM Corp, Armonk, NY, USA) or R software (Version 4.1.3, https://www.r-project.org/). The data were summarized using descriptive statistics and frequency tables. Categorical data were evaluated with the Chi-squared test. Independent risk factors for PNI in CC patients were examined by univariate and multivariate logistic regression analyses. Stepwise regression analysis was used to introduce variables into the model. Nomogram construction was performed based on the results of stepwise regression analysis. The receiver operating characteristic (ROC) curves and the corresponding area under the curve (AUC) values were used to evaluate the discrimination of the nomogram. The ROC curves, with AUC values quantified with the survival ROC package. Additionally, the calibration of the nomogram was evaluated by the corresponding calibration curves. The clinical benefits of the nomogram were assessed using decision curve analysis (DCA). Furthermore, Kaplan-Meier survival curve analysis and the log-rank test were used to determine whether there was a difference in overall survival and CC-specific survival between the PNI (−) group and the PNI (+) group. We used the R package cmprsk to identify statistical differences between CC-specific death and non-CC-specific death due to PNI by using competing risk model analysis and Gray’s test. We examined associations between factors and CC-specific death, including non-CC-specific death as a competing risk, using the Fine-Gray multivariable regression model. Statistical significance was established for a two-sided P value of 0.05.
Results
Patient characteristics
Among the 51,826 patients, 6,610 (12.75%) exhibited PNI, while 45,216 (87.25%) did not. By comparing the PNI (+) group and PNI (−) group, there were significant differences (P<0.05) in age, race, tumor site, tumor size, tumor grade, histology, T status, N status, CEA level, surgery status, RT status, ChemT status, bone metastasis, brain metastasis, liver metastasis, lung metastasis and marital status. The key information of the variables is shown in Table 1. The CC patients included in this study were randomly divided into a development cohort (n=36,280) and an internal validation cohort (n=15,546). There were no significant differences between the two groups according to the Chi-squared test (Table 2).
Table 1
Characteristics | Total (n=51,826), n (%) | PNI (−) (n=45,216), n (%) | PNI (+) (n=6,610), n (%) | P |
---|---|---|---|---|
Age (years) | <0.001 | |||
<40 | 1,328 (2.6) | 1,077 (2.4) | 251 (3.8) | |
40–59 | 13,118 (25.3) | 11,057 (24.4) | 2,061 (31.2) | |
60–79 | 26,060 (50.3) | 22,920 (50.7) | 3,140 (47.5) | |
≥80 | 11,320 (21.8) | 10,162 (22.5) | 1,158 (17.5) | |
Race | <0.001 | |||
White | 40,599 (78.3) | 35,541 (78.6) | 5,058 (76.5) | |
Black | 6,508 (12.6) | 5,562 (12.3) | 946 (14.3) | |
Other races | 4,719 (9.1) | 4,113 (9.1) | 606 (9.2) | |
Sex | 0.82 | |||
Male | 25,849 (49.9) | 22,543 (49.9) | 3,306 (50.0) | |
Female | 25,977 (50.1) | 22,673 (50.1) | 3,304 (50.0) | |
Tumor site | <0.001 | |||
Right colon | 32,764 (63.2) | 28,941 (64.0) | 3,823 (57.8) | |
Left colon | 19,062 (36.8) | 16,275 (36.0) | 2,787 (42.2) | |
Tumor size (cm) | <0.001 | |||
≤3 | 13,862 (26.7) | 12,607 (27.9) | 1,255 (19.0) | |
>3 to ≤5 | 18,938 (36.5) | 16,267 (36.0) | 2,671 (40.4) | |
>5 | 19,026 (36.7) | 16,342 (36.1) | 2,684 (40.6) | |
Grade | <0.001 | |||
Grade I | 3,636 (7.0) | 3,414 (7.6) | 222 (3.4) | |
Grade II | 37,560 (72.5) | 33,442 (74.0) | 4,118 (62.3) | |
Grade III | 8,780 (16.9) | 6,932 (15.3) | 1,848 (28.0) | |
Grade IV | 1,850 (3.6) | 1,428 (3.2) | 422 (6.4) | |
Histology | <0.001 | |||
Adenocarcinoma | 46,429 (89.6) | 40,471 (89.5) | 5,958 (90.1) | |
Mucinous adenocarcinoma | 4,826 (9.3) | 4,333 (9.6) | 493 (7.5) | |
Signet ring cell carcinoma | 571 (1.1) | 412 (0.9) | 159 (2.4) | |
T | <0.001 | |||
T1 | 4,649 (9.0) | 4,581 (10.1) | 68 (1.0) | |
T2 | 7,649 (14.8) | 7,442 (16.5) | 207 (3.1) | |
T3 | 29,670 (57.2) | 26,072 (57.7) | 3,598 (54.4) | |
T4 | 9,858 (19.0) | 7,121 (15.7) | 2,737 (41.4) | |
N | <0.001 | |||
N0 | 28,497 (55.0) | 27,039 (59.8) | 1,458 (22.1) | |
N1 | 14,110 (27.2) | 11,761 (26.0) | 2,349 (35.5) | |
N2 | 9,219 (17.8) | 6,416 (14.2) | 2,803 (42.4) | |
CEA | <0.001 | |||
Negative | 29,843 (57.6) | 27,089 (59.9) | 2,754 (41.7) | |
Positive | 21,983 (42.4) | 18,127 (40.1) | 3,856 (58.3) | |
Surgery | <0.001 | |||
Yes | 51,436 (99.2) | 44,835 (99.2) | 6,601 (99.9) | |
No | 390 (0.8) | 381 (0.8) | 9 (0.1) | |
RT | <0.001 | |||
Yes | 892 (1.7) | 714 (1.6) | 178 (2.7) | |
No | 50,934 (98.3) | 44,502 (98.4) | 6,432 (97.3) | |
ChemT | <0.001 | |||
Yes | 20,380 (39.3) | 16,217 (35.9) | 4,163 (63.0) | |
No | 31,446 (60.7) | 28,999 (64.1) | 2,447 (37.0) | |
Bone | <0.001 | |||
Yes | 199 (0.4) | 143 (0.3) | 56 (0.8) | |
No | 51,627 (99.6) | 45,073 (99.7) | 6,554 (99.2) | |
Brain | <0.001 | |||
Yes | 65 (0.1) | 45 (0.1) | 20 (0.3) | |
No | 51,761 (99.9) | 45,171 (99.9) | 6,590 (99.7) | |
Liver | <0.001 | |||
Yes | 5,611 (10.8) | 3,931 (8.7) | 1,680 (25.4) | |
No | 46,215 (89.2) | 41,285 (91.3) | 4,930 (74.6) | |
Lung | <0.001 | |||
Yes | 1,236 (2.4) | 873 (1.9) | 363 (5.5) | |
No | 50,590 (97.6) | 44,343 (98.1) | 6,247 (94.5) | |
Marital status | 0.04 | |||
Married | 28,904 (55.8) | 25,296 (55.9) | 3,608 (54.6) | |
Unmarried/SDW | 22,922 (44.2) | 19,920 (44.1) | 3,002 (45.4) |
CC, colon cancer; PNI, perineural invasion; CEA, carcinoembryonic antigen; RT, radiotherapy; ChemT, chemotherapy; SDW, separated and divorced and widowed.
Table 2
Characteristics | Development cohort (n=36,280), n (%) | Validation cohort (n=15,546), n (%) | P |
---|---|---|---|
Age (years) | 0.66 | ||
<40 | 917 (2.5) | 411 (2.6) | |
40–59 | 9,142 (25.2) | 3,976 (25.6) | |
60–79 | 18,272 (50.4) | 7,788 (50.1) | |
≥80 | 7,949 (21.9) | 3,371 (21.7) | |
Race | 0.36 | ||
White | 28,435 (78.4) | 12,164 (78.2) | |
Black | 4,514 (12.4) | 1,994 (12.8) | |
Other races | 3,331 (9.2) | 1,388 (8.9) | |
Sex | 0.058 | ||
Male | 17,996 (49.6) | 7,853 (50.5) | |
Female | 18,284 (50.4) | 7,693 (49.5) | |
Tumor site | 0.39 | ||
Right colon | 22,892 (63.1) | 9,872 (63.5) | |
Left colon | 13,388 (36.9) | 5,674 (36.5) | |
Tumor size (cm) | 0.11 | ||
≤3 | 9,654 (26.6) | 4,208 (27.1) | |
>3 to ≤5 | 13,362 (36.8) | 5,576 (35.9) | |
>5 | 13,264 (36.6) | 5,762 (37.1) | |
Grade | 0.41 | ||
Grade I | 2,521 (6.9) | 1,115 (7.2) | |
Grade II | 26,264 (72.4) | 11,296 (72.7) | |
Grade III | 6,206 (17.1) | 2,574 (16.6) | |
Grade IV | 1,289 (3.6) | 561 (3.6) | |
Histology | 0.14 | ||
Adenocarcinoma | 32,518 (89.6) | 13,911 (89.5) | |
Mucinous adenocarcinoma | 3,384 (9.3) | 1,442 (9.3) | |
Signet ring cell carcinoma | 378 (1.0) | 193 (1.2) | |
T | 0.25 | ||
T1 | 3,225 (8.9) | 1,424 (9.2) | |
T2 | 5,353 (14.8) | 2,296 (14.8) | |
T3 | 20,724 (57.1) | 8,946 (57.5) | |
T4 | 6,978 (19.2) | 2,880 (18.5) | |
N | 0.96 | ||
N0 | 19,938 (55.0) | 8,559 (55.1) | |
N1 | 9,878 (27.2) | 4,232 (27.2) | |
N2 | 6,464 (17.8) | 2,755 (17.7) | |
CEA | 0.48 | ||
Negative | 20,854 (57.5) | 8,989 (57.8) | |
Positive | 15,426 (42.5) | 6,557 (42.2) | |
PNI | 0.07 | ||
No | 31,588 (87.1) | 13,628 (87.7) | |
Yes | 4,692 (12.9) | 1,918 (12.3) | |
Surgery | 0.87 | ||
Yes | 36,009 (99.3) | 15,427 (99.2) | |
No | 271 (0.7) | 119 (0.8) | |
RT | 0.27 | ||
Yes | 609 (1.7) | 283 (1.8) | |
No | 35,671 (98.3) | 15,263 (98.2) | |
ChemT | 0.83 | ||
Yes | 14,278 (39.4) | 6,102 (39.3) | |
No | 22,002 (60.6) | 9,444 (60.7) | |
Bone | 0.90 | ||
Yes | 138 (0.4) | 61 (0.4) | |
No | 36,142 (99.6) | 15,485 (99.6) | |
Brain | 0.79 | ||
Yes | 44 (0.1) | 21 (0.1) | |
No | 36,236 (99.9) | 15,525 (99.9) | |
Liver | 0.29 | ||
Yes | 3,963 (10.9) | 1,648 (10.6) | |
No | 32,317 (89.1) | 13,898 (89.4) | |
Lung | 0.69 | ||
Yes | 872 (2.4) | 364 (2.3) | |
No | 35,408 (97.6) | 15,182 (97.7) | |
Marital status | 0.63 | ||
Married | 20,208 (55.7) | 8,696 (55.9) | |
Unmarried/SDW | 16,072 (44.3) | 6,850 (44.1) |
CEA, carcinoembryonic antigen; PNI, perineural invasion; RT, radiotherapy; ChemT, chemotherapy; SDW, separated and divorced and widowed.
Identification of risk factors of the patients with PNI in the training cohort
A univariate logistic regression analysis was performed to identify the risk factors for PNI in CC patients. The results revealed that age, race, tumor site, tumor size, grade, histology, T status, N status, CEA level, bone metastasis, brain metastasis, liver metastasis, and lung metastasis were independent risk factors for PNI (Table 3). Then, multivariate logistic regression analysis was performed, and the results showed that younger age, black race, tumor in the left colon, larger tumor size, higher tumor grade, higher T status, higher N status, higher CEA levels and liver metastasis were independently associated with PNI. Conversely, age ≥80 years and a histological type of mucinous adenocarcinoma were identified as protective factors against PNI (Table 3).
Table 3
Characteristics | Univariate analysis | Multivariate analysis | |||
---|---|---|---|---|---|
OR [95% CI] | P | OR [95% CI] | P | ||
Age (years) | |||||
<40 | Ref | Ref | |||
40–59 | 0.810 [0.683–0.966] | <0.001 | 0.975 [0.810–1.180] | 0.79 | |
60–79 | 0.582 [0.492–0.691] | 0.02 | 0.842 [0.701–1.016] | 0.07 | |
≥80 | 0.467 [0.391–0.560] | <0.001 | 0.761 [0.626–0.927] | 0.006 | |
Race | |||||
White | Ref | Ref | |||
Black | 1.235 [1.129–1.349] | <0.001 | 1.159 [1.050–1.277] | 0.003 | |
Other races | 1.062 [0.954–1.180] | 0.26 | 0.950 [0.846–1.064] | 0.38 | |
Sex | |||||
Male | Ref | NA | NA | ||
Female | 0.973 [0.915–1.035] | 0.38 | NA | NA | |
Tumor site | |||||
Right colon | Ref | Ref | |||
Left colon | 1.319 [1.239–1.404] | <0.001 | 1.274 [1.187–1.366] | <0.001 | |
Tumor size (cm) | |||||
≤3 | Ref | Ref | |||
>3 to ≤5 | 1.673 [1.538–1.822] | <0.001 | 0.869 [0.791–0.955] | 0.004 | |
>5 | 1.685 [1.548–1.834] | <0.001 | 0.700 [0.636–0.772] | <0.001 | |
Grade | |||||
Grade I | Ref | Ref | |||
Grade II | 1.902 [1.616–2.255] | <0.001 | 1.294 [1.089–1.547] | 0.004 | |
Grade III | 4.064 [3.429–4.849] | <0.001 | 1.882 [1.569–2.271] | <0.001 | |
Grade IV | 4.519 [3.678–5.572] | <0.001 | 1.886 [1.510–2.362] | <0.001 | |
Histology | |||||
Adenocarcinoma | Ref | Ref | |||
Mucinous adenocarcinoma | 0.817 [0.729–0.914] | <0.001 | 0.716 [0.633–0.807] | <0.001 | |
Signet ring cell carcinoma | 2.792 [2.225–3.482] | <0.001 | 1.220 [0.940–1.558] | 0.12 | |
T | |||||
T1 | Ref | Ref | |||
T2 | 2.156 [1.553–3.056] | <0.001 | 2.084 [1.496–2.963] | <0.001 | |
T3 | 9.949 [7.467–13.636] | <0.001 | 6.464 [4.807–8.927] | <0.001 | |
T4 | 28.762 [21.547–39.478] | <0.001 | 13.560 [10.034–18.807] | <0.001 | |
N | |||||
N0 | Ref | Ref | |||
N1 | 3.754 [3.459–4.076] | <0.001 | 2.383 [2.187–2.598] | <0.001 | |
N2 | 8.296 [7.644–9.008] | <0.001 | 4.094 [3.740–4.484] | <0.001 | |
CEA | |||||
Negative | Ref | Ref | |||
Positive | 2.085 [1.959–2.219] | <0.001 | 1.194 [1.112–1.283] | <0.001 | |
Bone | |||||
Yes | Ref | Ref | |||
No | 2.478 [1.678–3.581] | <0.001 | 0.708 [0.466–1.056] | 0.10 | |
Brain | |||||
Yes | Ref | Ref | |||
No | 3.148 [1.619–5.829] | <0.001 | 1.390 [0.684–2.708] | 0.35 | |
Liver | |||||
Yes | Ref | Ref | |||
No | 3.706 [3.434–3.998] | <0.001 | 1.696 [1.551–1.855] | <0.001 | |
Lung | |||||
Yes | Ref | Ref | |||
No | 2.868 [2.466–3.326] | <0.001 | 1.082 [0.914–1.279] | 0.36 | |
Marital status | |||||
Married | Ref | NA | NA | ||
Unmarried/SDW | 1.026 [0.964–1.091] | 0.42 | NA | NA |
CC, colon cancer; OR, odds ratio; CI, confidence interval; Ref, reference; NA, not applicable; CEA, carcinoembryonic antigen; SDW, separated and divorced and widowed.
Development and validation of the nomogram to predict PNI probability
To further predict the probability of PNI in CC patients, a nomogram was developed with the 11 significant factors associated with PNI in the stepwise multivariate logistic regression (Figure 1). The AUCs for the development and validation cohorts were 0.787 [95% confidence interval (CI): 0.780–0.793] and 0.781 (95% CI: 0.771–0.791), respectively (Figure 2). In addition, the calibration curves illustrated that our nomogram has strong calibration (Figure 3). Furthermore, the DCA curves also indicated that the nomogram has good clinical practicability in the development cohort and internal validation cohort (Figure 4).
![Click on image to zoom](http://cdn.amegroups.cn/journals/pbpc/files/journals/3/articles/95697/public/95697-PB12-9710-R1.jpg/w300)
![Click on image to zoom](http://cdn.amegroups.cn/journals/pbpc/files/journals/3/articles/95697/public/95697-PB7-1680-R1.jpg/w300)
![Click on image to zoom](http://cdn.amegroups.cn/journals/pbpc/files/journals/3/articles/95697/public/95697-PB8-6327-R1.jpg/w300)
![Click on image to zoom](http://cdn.amegroups.cn/journals/pbpc/files/journals/3/articles/95697/public/95697-PB9-5426-R1.jpg/w300)
Kaplan-Meier survival analysis
After analyzing the risk factors for PNI in patients with CC, overall survival and CC-specific survival in CC patients were explored using Kaplan-Meier analysis. Among all 51,826 patients included in this study, 19,165 (36.98%) died. The cumulative incidence rate of CC-specific death was 24.69% (12,794/51,826), and the cumulative non-CC-specific death incidence rate was 12.29% (6,371/51,826). Compared with patients in the PNI (−) group, patients in the PNI (+) group had worse overall survival and CC-specific survival (Figure 5).
![Click on image to zoom](http://cdn.amegroups.cn/journals/pbpc/files/journals/3/articles/95697/public/95697-PB10-7792-R1.jpg/w300)
Competing risk model of CC-specific death and non-CC-specific death
A total of 19,165 patients who died were included in our cohort, of which 66.76% (12,794/19,165) were CC-specific death and 33.24% (6,371/19,165) were non-CC-specific death. Compared with patients in the PNI (−) group, patients in the PNI (+) group had a higher cumulative CC-specific death rate (Gray’s test, P<0.001), and there was no significant difference in the cumulative non-CC-specific death rate between the two groups (Gray’s test, P=0.63) (Figure 6).
![Click on image to zoom](http://cdn.amegroups.cn/journals/pbpc/files/journals/3/articles/95697/public/95697-PB11-5806-R1.jpg/w300)
Multivariable competing risk analysis of survival
The Fine-Gray multivariable regression model was constructed to investigate independent prognostic factors of CC-specific death (Table 4). The results revealed that patients with PNI had a higher CC-specific death rate [hazard ratio (HR) =1.243; 95% CI: 1.183–1.305; P<0.001] than patients without PNI. Compared with the corresponding subgroups, patients with female sex, tumor in the left colon, and no RT tended to have significantly lower CC-specific death (P<0.05). Conversely, patients with older age, black race, highly differentiated tumors, mucinous adenocarcinoma and signet ring cell carcinoma, higher T status, higher N status, elevated CEA levels, bone metastasis, liver metastasis, lung metastasis, no surgery, no ChemT tended to have significantly higher CC-specific death (P<0.05). Additionally, age, race, sex, tumor site, histology, AJCC T status, CEA levels, liver metastasis, lung metastasis, marital status, ChemT status and RT status were also related to non-CC-specific death (P<0.05).
Table 4
Characteristics | CCSD (n=12,794) | Non-CCSD (n=6,371) | |||
---|---|---|---|---|---|
HR [95% CI] | P | HR [95% CI] | P | ||
Age (years) | |||||
<40 | Ref | Ref | |||
40–59 | 1.083 [0.971–1.206] | 0.15 | 1.921 [1.264–2.920] | 0.002 | |
60–79 | 1.375 [1.236–1.529] | <0.001 | 4.946 [3.278–7.464] | <0.001 | |
≥80 | 1.978 [1.766–2.215] | <0.001 | 10.972 [7.263–16.575] | <0.001 | |
Race | |||||
White | Ref | Ref | |||
Black | 1.139 [1.078–1.203] | <0.001 | 1.101 [1.017–1.192] | 0.02 | |
Other races | 0.909 [0.849–0.973] | 0.006 | 0.827 [0.750–0.911] | <0.001 | |
Sex | |||||
Male | Ref | Ref | |||
Female | 0.927 [0.892–0.962] | <0.001 | 0.706 [0.671–0.744] | <0.001 | |
Tumor site | |||||
Right colon | Ref | Ref | |||
Left colon | 0.919 [0.884–0.956] | <0.001 | 0.915 [0.865–0.967] | 0.002 | |
Tumor size (cm) | |||||
≤3 | Ref | Ref | |||
>3 to ≤5 | 0.972 [0.921–1.026] | 0.31 | 1.016 [0.951–1.086] | 0.63 | |
>5 | 1.050 [0.994–1.109] | 0.08 | 0.933 [0.868–1.003] | 0.06 | |
Grade | |||||
Grade I | Ref | Ref | |||
Grade II | 1.084 [0.997–1.179] | 0.059 | 0.975 [0.889–1.069] | 0.59 | |
Grade III | 1.384 [1.263–1.516] | <0.001 | 0.965 [0.864–1.078] | 0.53 | |
Grade IV | 1.534 [1.366–1.722] | <0.001 | 1.023 [0.869–1.204] | 0.79 | |
Histology | |||||
Adenocarcinoma | Ref | Ref | |||
Mucinous adenocarcinoma | 1.149 [1.081–1.221] | <0.001 | 1.096 [1.010–1.188] | 0.03 | |
Signet ring cell carcinoma | 1.640[1.428–1.884] | <0.001 | 1.001 [0.772–1.298] | >0.99 | |
PNI | |||||
No | Ref | Ref | |||
Yes | 1.243 [1.183–1.305] | <0.001 | 0.920 [0.838–1.009] | 0.08 | |
T | |||||
T1 | Ref | Ref | |||
T2 | 0.913 [0.800–1.042] | 0.18 | 1.073 [0.971–1.186] | 0.17 | |
T3 | 1.749 [1.559–1.962] | <0.001 | 1.038 [0.942–1.145] | 0.45 | |
T4 | 3.162 [2.800–3.572] | <0.001 | 0.869 [0.770–0.982] | 0.02 | |
N | |||||
N0 | Ref | Ref | |||
N1 | 2.136 [2.024–2.254] | <0.001 | 0.996 [0.929–1.068] | 0.91 | |
N2 | 3.352 [3.163–3.552] | <0.001 | 0.910 [0.829–1.000] | 0.050 | |
CEA | |||||
Negative | Ref | Ref | |||
Positive | 1.531 [1.471–1.594] | <0.001 | 1.274 [1.208–1.343] | <0.001 | |
Bone | |||||
Yes | Ref | Ref | |||
No | 1.336 [1.061.682] | 0.01 | 1.421 [0.856–2.359] | 0.17 | |
Brain | |||||
Yes | Ref | Ref | |||
No | 1.179 [0.589–2.361] | 0.64 | 1.241 [0.579–2.660] | 0.58 | |
Liver | |||||
Yes | Ref | Ref | |||
No | 3.340 [3.179–3.510] | <0.001 | 0.548 [0.481–0.625] | <0.001 | |
Lung | |||||
Yes | Ref | Ref | |||
No | 1.576 [1.444–1.720] | <0.001 | 0.728 [0.561–0.945] | 0.02 | |
Marital status | |||||
Married | Ref | Ref | |||
Unmarried/SDW | 1.194 [1.149–1.241] | <0.001 | 1.290 [1.225–1.359] | <0.001 | |
Surgery | |||||
Yes | Ref | Ref | |||
No | 4.299 [3.410–5.421] | <0.001 | 1.195 [0.859–1.663] | 0.29 | |
ChemT | |||||
Yes | Ref | Ref | |||
No | 1.321 [1.256–1.390] | <0.001 | 1.955 [1.810–2.112] | <0.001 | |
RT | |||||
Yes | Ref | Ref | |||
No | 0.835 [0.744–0.936] | 0.002 | 0.738 [0.582–0.937] | 0.01 |
CC, colon cancer; CCSD, colon cancer-specific death; HR, hazard ratio; CI, confidence interval; non-CCSD, non-colon cancer-specific death; Ref, reference; PNI, perineural invasion; CEA, carcinoembryonic antigen; SDW, separated and divorced and widowed; ChemT, chemotherapy; RT, radiotherapy.
Competing risk model of the effect of different treatment methods after surgery on the survival of pT4 CC patients
Since adjuvant RT is not traditionally used for early CC, we then explored the efficiency of adjuvant therapy on the prognosis of pT4 patients with or without PNI. A total of 5,832 patients with pT4 stage disease who died were included in this analysis, of which 32.61% (1,901/5,830) were in the PNI (+) subgroup, and 67.39% (3,929/5,830) were in the PNI (−) subgroup. As shown in Figure 7, PNI (+) pT4 CC patients treated with CCRP or ChemT had a better cumulative CC-specific death incidence (Gray’s test, P<0.001) than those in the non-therapy group (Figure 7A). In contrast, patients in the CCRP group had a better cumulative CC-specific death incidence (Gray’s test, P<0.001) than did the ChemT group and non-therapy group in PNI (−) pT4 CC patients (Figure 7B).
![Click on image to zoom](http://cdn.amegroups.cn/journals/pbpc/files/journals/3/articles/95697/public/95697-PB12-8038-R1.jpg/w300)
Discussion
In the present study, based on the analysis of a cohort of 51,826 patients with CC in the SEER database from January 2010 to December 2015, we developed a nomogram to predict PNI based on the clinical features of patients, and we also validated that PNI could significantly reduce CC-specific survival for CC patients. pT4 CC patients without PNI treated with CCRP had a lower CC-specific death rate than patients with ChemT or non-therapy. The nomogram in this study demonstrated a high degree of discrimination and calibration, and this study is the first real-world study to examine the impact of PNI on CC patients by using the Fine-Gray multivariable regression model.
PNI was reported for the first time in 1800s, and its prognostic value has gradually attracted attention since the 1970s (18). PNI is described as tumor invasion in, around and through the nerves (6). According to a previous study, 33% of patients with CRC exhibit PNI (19). In this study, the prevalence of PNI in CC patients was approximately up to 12.75%. Previous studies have reported that PNI is associated with more severe clinicopathologic characteristics (7,10,20). Our study found that PNI was significantly correlated with tumor site, tumor size, tumor grade, tumor histology, T status, N status, CEA levels, bone metastasis, brain metastasis, liver metastasis and lung metastasis.
PNI is related to poor prognosis for CRC patients, and early identification of PNI in patients with CRC is necessary for the development of individualized clinical therapeutic strategies (9,21). Currently, PNI in CC can only be assessed by histopathological examination after surgery (22). Therefore, it is essential to construct a prediction model for PNI. As a clinical tool, nomograms have become widely used as prediction models in recent years (23,24). Huang et al. selected clinical features by least absolute shrinkage and selection operator (LASSO) logistic regression analysis and constructed a nomogram based on a population in a single center to predict PNI status in CRC patients with AUCs of 0.704 and 0.692 in the training and validation cohorts, respectively (25). In our study, to construct a nomogram for predicting PNI in CC patients from the SEER database, significant clinical features were selected using stepwise logistic regression. The nomogram was constructed using eleven features: age, race, tumor grade, tumor site, histology, T status, N status, CEA levels, liver metastasis, bone metastasis and tumor size. The nomogram demonstrated good discrimination and calibration with AUCs of 0.787 and 0.781 in the training cohort and validation cohort, respectively. The difference between the two nomograms could be mainly attributed to the different features used for model construction. In our study, more clinical features were included, and tumor stage was the main feature of the model. Thus, we think that the nomogram has certain application value in the clinic for developing therapeutic strategies.
Researchers have found that PNI is related to a poor prognosis in CRC patients (7,21,26). PNI has been included in the AJCC cancer staging system as a negative prognostic factor, and PNI has been presented in pathology reports for years (27). Our study validated the prognostic value of PNI based on the data of CC patients from the SEER database. Kaplan-Meier curve analysis indicated that patients in the PNI (+) group had poorer overall survival and CC-specific survival than patients in the PNI (−) group, which is consistent with findings from past studies (7,28). Tu et al. found that PNI was a negative prognostic factor for CC in stage I (HR =1.590; 95% CI: 0.951–2.658) and stage II (HR =1.607; 95% CI: 1.426–1.812) (7). In this study, we demonstrated that the PNI (+) group had a higher CC-specific death rate (HR =1.243; 95% CI: 1.183–1.305) than the PNI (−) group, and there was no statistically difference in the non-CC-specific death rate between the two groups. The difference in HR between the two studies could be attributed to the fact that the occurrence of competitive events might affect CC-specific survival and preclude the primary event of CC-specific death. Compared to Fine-Gray regression, Cox regression overestimated the cumulative event rate and yielded a higher hazard ratio. Additionally, the differences might be occurred because of the different stage of tumors. Tu et al. performed their study mainly in CC patients in stage I and stage II, while in our study, CC patients in all stages were included (7).
Currently, CC can be treated with a combination of surgery, RT, ChemT and targeted therapy (29,30). The role of RT for the treatment of CC remains controversial (12,31). In our study, after multivariable competing risk analysis, CC patients who underwent surgery and ChemT had a lower cumulative CC-specific death rate than others, which is consistent with previous studies (22,32). Some previous studies reported that additional RT after surgery conferred a survival benefit in T4 CC patients (15,33), while other studies reported that RT had little effect on prognosis (31,34). In this study, the results indicated that CC patients who underwent RT had a higher cumulative CC-specific death rate than patients without RT. This may be due to selection bias and the preference of patients with RT who may have a higher tumor stage and a poorer prognosis. Therefore, we further explored the impact of adjuvant therapy on the prognosis of pT4 stage patients with or without PNI. The results showed that patients with CCRP had lower cumulative CC-specific death incidence than patients treated with ChemT or non-therapy in the PNI (−) pT4 CC patients, and there was no statistically difference between the CCRP group and the ChemT group in PNI (+) pT4 CC patients. Collectively, these findings indicate that CC patients can benefit from surgery and ChemT, and pT4 patients without PNI might benefit from CCRP.
Inevitably, this study has also some notable limitations that need to be addressed. First, as a retrospective study design, unregistered confounding factors and selection bias cannot be avoided, and study external validity may be limited. Second, it is possible that risk estimates could be influenced by potential confounders, such as family history, comorbidities, surgery methods, or the exact scope of lymph node dissection. Third, new biomarkers have been widely reported to have prognostic values in CC recently (35-37), but they could not be included in this study because of missing related data. Finally, our model has internal validation, but does not have external validation to confirm its accuracy. An external multicenter dataset needs to be used to validate our results. Although the data from SEER database indicates that RT was conducted after radical resection, yet we are unable to determine the sequence of RT and ChemT. New clinical studies should be conducted to see if CCRP after surgery can improve prognostic of pT4 CC patients without PNI.
Conclusions
Our study develops a nomogram for predicting the risk of PNI in patients with CC before operation. Our study also demonstrates that the presence of PNI is independently associated with an increase in the CC-specific death rate and could be considered a reliable prognostic factor in CC patients. Finally, our study indicates that pT4 patients without PNI might benefit from CCRP after surgery.
Acknowledgments
None.
Footnote
Reporting Checklist: The authors have completed the TRIPOD reporting checklist. Available at https://tcr.amegroups.com/article/view/10.21037/tcr-24-1030/rc
Peer Review File: Available at https://tcr.amegroups.com/article/view/10.21037/tcr-24-1030/prf
Funding: This work was supported by
Conflicts of Interest: Both authors have completed the ICMJE uniform disclosure form (available at https://tcr.amegroups.com/article/view/10.21037/tcr-24-1030/coif). 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).
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
- Sung H, Ferlay J, Siegel RL, et al. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin 2021;71:209-49. [Crossref] [PubMed]
- Siegel RL, Miller KD, Fuchs HE, et al. Cancer statistics, 2022. CA Cancer J Clin 2022;72:7-33. [Crossref] [PubMed]
- Cho SS, Park JW, Kang GH, et al. Prognostic Impact of Extramural Lymphatic, Vascular, and Perineural Invasion in Stage II Colon Cancer: A Comparison With Intramural Invasion. Dis Colon Rectum 2023;66:366-73. [Crossref] [PubMed]
- Nikberg M, Chabok A, Letocha H, et al. Lymphovascular and perineural invasion in stage II rectal cancer: a report from the Swedish colorectal cancer registry. Acta Oncol 2016;55:1418-24. [Crossref] [PubMed]
- Parnaby CN, Scott NW, Ramsay G, et al. Prognostic value of lymph node ratio and extramural vascular invasion on survival for patients undergoing curative colon cancer resection. Br J Cancer 2015;113:212-9. [Crossref] [PubMed]
- Liu Z, Luo C, Chen X, et al. Noninvasive prediction of perineural invasion in intrahepatic cholangiocarcinoma by clinicoradiological features and computed tomography radiomics based on interpretable machine learning: a multicenter cohort study. Int J Surg 2024;110:1039-51. [Crossref] [PubMed]
- Tu J, Yao Z, Wu W, et al. Perineural Invasion Is a Strong Prognostic Factor but Not a Predictive Factor of Response to Adjuvant Chemotherapy in Node-Negative Colon Cancer. Front Oncol 2021;11:663154. [Crossref] [PubMed]
- Cao Y, Deng S, Yan L, et al. Perineural invasion is associated with poor prognosis of colorectal cancer: a retrospective cohort study. Int J Colorectal Dis 2020;35:1067-75. [Crossref] [PubMed]
- Wang H, Huo R, He K, et al. Perineural invasion in colorectal cancer: mechanisms of action and clinical relevance. Cell Oncol (Dordr) 2024;47:1-17. [Crossref] [PubMed]
- Benson AB, Venook AP, Adam M, et al. Colon Cancer, Version 3.2024, NCCN Clinical Practice Guidelines in Oncology. J Natl Compr Canc Netw 2024;22:e240029. [Crossref] [PubMed]
- Agas RAF, Co LBA, Sogono PG, et al. Assessing the Effect of Radiotherapy in Addition to Surgery in Colon Adenocarcinomas: a Systematic Review and Meta-analysis of Contemporary Evidence. J Gastrointest Cancer 2020;51:445-60. [Crossref] [PubMed]
- Lin Y, Wang L, Shao L, et al. Prognostic analysis and beneficiary identification of adjuvant external beam radiotherapy for stage pT4b sigmoid colon cancer. Sci Rep 2021;11:11782. [Crossref] [PubMed]
- Oki E, Ando K, Taniguchi H, et al. Sustainable Clinical Development of Adjuvant Chemotherapy for Colon Cancer. Ann Gastroenterol Surg 2022;6:37-45. [Crossref] [PubMed]
- Wu CY, Ye K. Risk factors for the prognosis of colon cancer. World J Gastrointest Oncol 2024;16:3738-40. [Crossref] [PubMed]
- Wegner RE, Abel S, Monga D, et al. Utilization of Adjuvant Radiotherapy for Resected Colon Cancer and Its Effect on Outcome. Ann Surg Oncol 2020;27:825-32. [Crossref] [PubMed]
- Zhang H, Li Y, Xia F, et al. Study protocol of short-course radiotherapy combined with CAPOX and PD-1 inhibitor for locally advanced colon cancer: a randomised, prospective, multicentre, phase II trial (TORCH-C). BMJ Open 2024;14:e079442. [Crossref] [PubMed]
- Wingo PA, Jamison PM, Hiatt RA, et al. Building the infrastructure for nationwide cancer surveillance and control--a comparison between the National Program of Cancer Registries (NPCR) and the Surveillance, Epidemiology, and End Results (SEER) Program (United States). Cancer Causes Control 2003;14:175-93. [Crossref] [PubMed]
- Liebig C, Ayala G, Wilks J, et al. Perineural invasion is an independent predictor of outcome in colorectal cancer. J Clin Oncol 2009;27:5131-7. [Crossref] [PubMed]
- Liebig C, Ayala G, Wilks JA, et al. Perineural invasion in cancer: a review of the literature. Cancer 2009;115:3379-91. [Crossref] [PubMed]
- van Wyk HC, Going J, Horgan P, et al. The role of perineural invasion in predicting survival in patients with primary operable colorectal cancer: A systematic review. Crit Rev Oncol Hematol 2017;112:11-20. [Crossref] [PubMed]
- Knijn N, Mogk SC, Teerenstra S, et al. Perineural Invasion is a Strong Prognostic Factor in Colorectal Cancer: A Systematic Review. Am J Surg Pathol 2016;40:103-12. [Crossref] [PubMed]
- Leijssen LGJ, Dinaux AM, Taylor MS, et al. Perineural Invasion Is a Prognostic but not a Predictive Factor in Nonmetastatic Colon Cancer. Dis Colon Rectum 2019;62:1212-21. [Crossref] [PubMed]
- Park SY. Nomogram: An analogue tool to deliver digital knowledge. J Thorac Cardiovasc Surg 2018;155:1793. [Crossref] [PubMed]
- Iasonos A, Schrag D, Raj GV, et al. How to build and interpret a nomogram for cancer prognosis. J Clin Oncol 2008;26:1364-70. [Crossref] [PubMed]
- Huang X, Liu J, Wu G, et al. Development and Validation of a Nomogram for Preoperative Prediction of Perineural Invasion in Colorectal Cancer. Med Sci Monit 2019;25:1709-17. [Crossref] [PubMed]
- Yang Y, Huang X, Sun J, et al. Prognostic value of perineural invasion in colorectal cancer: a meta-analysis. J Gastrointest Surg 2015;19:1113-22. [Crossref] [PubMed]
- Compton C, Fenoglio-Preiser CM, Pettigrew N, et al. American Joint Committee on Cancer Prognostic Factors Consensus Conference: Colorectal Working Group. Cancer 2000;88:1739-57. [Crossref] [PubMed]
- Alotaibi AM, Lee JL, Kim J, et al. Prognostic and Oncologic Significance of Perineural Invasion in Sporadic Colorectal Cancer. Ann Surg Oncol 2017;24:1626-34. [Crossref] [PubMed]
- Long X, Deng Z, Li G, et al. Identification of critical genes to predict recurrence and death in colon cancer: integrating gene expression and bioinformatics analysis. Cancer Cell Int 2018;18:139. [Crossref] [PubMed]
- Yuan Y, Chen J, Wang J, et al. Development and Clinical Validation of a Novel 4-Gene Prognostic Signature Predicting Survival in Colorectal Cancer. Front Oncol 2020;10:595. [Crossref] [PubMed]
- Sebastian NT, Tan Y, Miller ED, et al. Surgery with and without adjuvant radiotherapy is associated with similar survival in T4 colon cancer. Colorectal Dis 2020;22:779-89. [Crossref] [PubMed]
- Cohen AM. Surgical considerations in patients with cancer of the colon and rectum. Semin Oncol 1991;18:381-7.
- Huang Y, Gu X, Ge K, et al. The survival benefit of adjuvant radiotherapy for pathological T4N2M0 colon cancer in the Modern Chemotherapy Era: evidence from the SEER database 2004-2015. Artif Cells Nanomed Biotechnol 2020;48:834-40. [Crossref] [PubMed]
- Xiao J, Caan BJ, Cespedes Feliciano EM, et al. Association of Low Muscle Mass and Low Muscle Radiodensity With Morbidity and Mortality for Colon Cancer Surgery. JAMA Surg 2020;155:942-9. [Crossref] [PubMed]
- Barras D, Missiaglia E, Wirapati P, et al. BRAF V600E Mutant Colorectal Cancer Subtypes Based on Gene Expression. Clin Cancer Res 2017;23:104-15. [Crossref] [PubMed]
- Xu J, Zhang R, Peng Q, et al. The profile and prognostic value of circulating lymphocyte subsets in metastatic colon cancer. Int Immunopharmacol 2023;117:109848. [Crossref] [PubMed]
- Lv J, Li W, Wang X, et al. Identification of MKI67, TPR, and TCHH Mutations as Prognostic Biomarkers for Patients With Defective Mismatch Repair Colon Cancer Stage II/III. Dis Colon Rectum 2023;66:1481-91. [Crossref] [PubMed]