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
Original Article

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

Zhongqiang Zheng1, Xuanzi Sun2

1Department of General Surgery, The First Affiliated Hospital of Xi’an Jiaotong University, Xi’an, China; 2Department of Radiation Oncology, The First Affiliated Hospital of Xi’an Jiaotong University, Xi’an, China

Contributions: (I) Conception and design: Both authors; (II) Administrative support: X Sun; (III) Provision of study materials or patients: Both authors; (IV) Collection and assembly of data: Z Zheng; (V) Data analysis and interpretation: Z Zheng; (VI) Manuscript writing: Both authors; (VII) Final approval of manuscript: Both authors.

Correspondence to: Xuanzi Sun, MD. Department of Radiation Oncology, The First Affiliated Hospital of Xi’an Jiaotong University, No. 277 Yanta West Road, Xi’an 710061, China. Email: sunxuanzi12@163.com.

Background: Perineural invasion (PNI) in colon cancer (CC) is widely associated with poor prognosis. In this study, we aimed to develop a predictive model for PNI and to assess its prognostic value in CC patients.

Methods: Data for CC patients with or without PNI were obtained from the Surveillance, Epidemiology, and End Results (SEER) database from 2010 to 2015. Potential features were selected by stepwise logistic regression, and multivariate logistic regression was used to develop the nomogram. Nomogram performance was assessed based on its calibration curve, discrimination ability and clinical utility. The prognostic value of PNI was assessed using Kaplan-Meier analysis, a competing risk model, and a Fine-Gray multivariable regression model.

Results: A total of 51,826 subjects were included in the study. The nomogram consisted of 11 features was constructed, which provided good calibration and discrimination with area under the curve values of 0.787 vs. 0.781 (development cohort vs. validation cohort). Patients with PNI had worse CC-specific survival (P<0.001) and a higher CC-specific death rate (Gray’s test, P<0.001) than patients without PNI. Fine-Gray multivariable regression analysis showed that patients with PNI had a higher CC-specific death rate than patients without PNI [hazard ratio (HR) =1.243; 95% confidence interval (CI): 1.183–1.305; P<0.001]. Pathologic stage T4 (pT4) CC patients without PNI treated with chemotherapy (ChemT) plus radiotherapy (RT) had a lower CC-specific death rate than ChemT-treated or non-therapy patients.

Conclusions: The nomogram developed herein has certain clinical application value for predicting PNI risk in CC patients. PNI is a survival predictor for CC patients. pT4 patients without PNI might benefit from combined ChemT and RT.

Keywords: Colon cancer (CC); Surveillance, Epidemiology, and End Results (SEER); perineural invasion (PNI); nomogram; competing risk model


Submitted Jun 20, 2024. Accepted for publication Dec 04, 2024. Published online Jan 23, 2025.

doi: 10.21037/tcr-24-1030


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

Baseline clinical characteristics of patients with CC

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

Patients’ clinical characteristics in the development and validation cohort

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

Logistic regression analysis of the risk factors for patients with CC

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).

Figure 1 Nomogram for predicting the probability of PNI in CC patients. RCC, right-side colon cancer; LCC, left-side colon cancer; CEA, carcinoembryonic antigen; PNI, perineural invasion; CC, colon cancer.
Figure 2 The ROC curves of nomogram for predicting PNI in the development cohort (A) and validation cohort (B). AUC, area under the curve; ROC, receiver operating characteristic; PNI, perineural invasion.
Figure 3 The calibration curves of the nomogram for predicting PNI in the development cohort (A) and validation cohort (B). ROC, receiver operating characteristic; PNI, perineural invasion.
Figure 4 The DCA curves of the nomogram for predicting PNI in the development cohort (A) and validation cohort (B). DCA, decision curve analysis; PNI, perineural invasion.

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).

Figure 5 Kaplan-Meier survival analysis for PNI (+) and PNI (−) CC patients. (A) Overall survival curves for the PNI (+) and PNI (−) group. (B) CC-specific survival curves for the PNI (+) and PNI (−) CC patients. PNI, perineural invasion; CC, colon cancer.

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).

Figure 6 Cumulative incidence of CCSD and non-CCSD in PNI (+) and PNI (−) CC patients. PNI, perineural invasion; CCSD, colon cancer-specific death; CC, colon cancer.

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

Multivariable competing risk analysis in CC patients

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).

Figure 7 Competing risk analysis of the effect of different treatment methods after surgery on the survival in pT4 CC patients. (A) Cumulative incidence of CCSD and non-CCSD in PNI (+) CC patients of different treatments. (B) Cumulative incidence of CCSD and non-CCSD in PNI (−) CC patients of different treatments. CCSD, colon cancer-specific death; ChemT, chemotherapy; CCRP, combined chemotherapy and radiotherapy; pT4, pathologic stage T4; CC, colon cancer; PNI, perineural invasion.

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 Shaanxi Province Science Foundation for Youths (No. 2024JC-YBQN-0805).

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

  1. 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]
  2. Siegel RL, Miller KD, Fuchs HE, et al. Cancer statistics, 2022. CA Cancer J Clin 2022;72:7-33. [Crossref] [PubMed]
  3. 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]
  4. 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]
  5. 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]
  6. 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]
  7. 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]
  8. 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]
  9. 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]
  10. 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]
  11. 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]
  12. 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]
  13. 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]
  14. Wu CY, Ye K. Risk factors for the prognosis of colon cancer. World J Gastrointest Oncol 2024;16:3738-40. [Crossref] [PubMed]
  15. 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]
  16. 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]
  17. 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]
  18. 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]
  19. Liebig C, Ayala G, Wilks JA, et al. Perineural invasion in cancer: a review of the literature. Cancer 2009;115:3379-91. [Crossref] [PubMed]
  20. 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]
  21. 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]
  22. 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]
  23. Park SY. Nomogram: An analogue tool to deliver digital knowledge. J Thorac Cardiovasc Surg 2018;155:1793. [Crossref] [PubMed]
  24. 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]
  25. 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]
  26. 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]
  27. 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]
  28. 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]
  29. 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]
  30. 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]
  31. 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]
  32. Cohen AM. Surgical considerations in patients with cancer of the colon and rectum. Semin Oncol 1991;18:381-7.
  33. 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]
  34. 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]
  35. 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]
  36. 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]
  37. 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]
Cite this article as: Zheng Z, Sun X. 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. Transl Cancer Res 2025;14(1):141-158. doi: 10.21037/tcr-24-1030

Download Citation