ORIGINAL RESEARCH article

Front. Oncol., 14 July 2023

Sec. Gastrointestinal Cancers: Colorectal Cancer

Volume 13 - 2023 | https://doi.org/10.3389/fonc.2023.1182944

Prognostic scoring system based on eosinophil- and basophil-related markers for predicting the prognosis of patients with stage II and stage III colorectal cancer: a retrospective cohort study

  • 1. Department of Epidemiology, School of Public Health, NHC Key Laboratory of Etiology and Epidemiology (23618504), Harbin Medical University, Harbin, Heilongjiang, China

  • 2. Department of Colorectal Surgery, Harbin Medical University Cancer Hospital, Harbin, Heilongjiang, China

  • 3. Division of Chronic and Non-communicable Diseases, Harbin Center for Diseases Control and Prevention, Harbin, Heilongjiang, China

Article metrics

View details

8

Citations

2,7k

Views

781

Downloads

Abstract

Background:

Systemic inflammation is associated with the prognosis of colorectal cancer (CRC). The current study aimed to construct a comprehensively inflammatory prognostic scoring system named risk score (RS) based on eosinophil- and basophil-related markers and assess its prognostic value in patients with stage II and stage III CRC.

Patients and methods:

A total of 3,986 patients were enrolled from January 2007 to December 2013. The last follow-up time was January 2019. They were randomly assigned to the training set and testing set in a 3:2 split ratio. Least absolute shrinkage and selection operator (LASSO)–Cox regression analysis was performed to select the optimal prognostic factors in the construction of RS. The Kaplan–Meier curve, time-dependent receiver operating characteristic (ROC), and Cox analysis were used to evaluate the association between RS and overall survival (OS).

Results:

In the training set, all inflammatory markers showed certain prognostic values. Based on LASSO-Cox analysis, nine markers were integrated to construct RS. The Kaplan–Meier curve showed that a higher RS (RS > 0) had a significantly worse prognosis (log-rank p< 0.0001). RS (>0) remained an independent prognostic factor for OS (hazard ratio (HR): 1.70, 95% confidence interval (CI), 1.43–2.03, p< 0.001). The prognostic value of RS was validated in the entire cohort. Time-dependent ROC analysis showed that RS had a stable prognostic effect throughout the follow-up times and could enhance the prognostic ability of the stage by combination. Nomogram was established based on RS and clinicopathological factors for predicting OS in the training set and validated in the testing set. The area under the curve (AUC) values of the 3-year OS in the training and testing sets were 0.748 and 0.720, respectively. The nomogram had a satisfactory predictive accuracy and had better clinical application value than the tumor stage alone.

Conclusions:

RS might be an independent prognostic factor for OS in patients with stage II and III CRC, which is helpful for risk stratification of patients. Additionally, the nomogram might be used for personalized prediction and might contribute to formulating a better clinical treatment plan.

1 Introduction

Colorectal cancer (CRC), including colon cancer and rectal cancer, is one of the most common malignant tumors threatening human health. Based on the GLOBOCAN 2020 estimation (1), CRC is the second leading cause of death, only next to lung cancer, with more than 935,000 deaths. There were 51,020 CRC deaths in the United States in 2019, equivalent to 8.4% of all cancer deaths (2). In China, 191,000 patients died of CRC, which ranks fifth in cancer deaths in 2015 (3). Surgical resection is the most common therapy for patients with CRC, and tumor stage based on pathological characteristics is widely used for evaluating the prognosis of CRC patients (4, 5). However, studies have shown that the prognosis with the same stage varies greatly and is highly heterogeneous (6). In addition, the pathological characteristics are mainly obtained by biopsy or pathological reports, which are hard to represent the overall condition of the tumor. Thus, it is necessary to identify a non-invasive and more accurate marker to assess the prognosis of CRC patients.

Recently, increasing evidence has demonstrated that systemic inflammation is closely associated with the progression and prognosis of CRC (79). The inflammatory markers based on peripheral leukocytes include neutrophil, lymphocyte, and monocyte counts, and related indicators such as neutrophil-to-lymphocyte ratio and lymphocyte-to-monocyte ratio, which have been widely studied to predict the prognosis of patients with CRC as well as other types of malignant tumors (1015). Eosinophils and basophils account for a small proportion of circulating leukocytes in the bloodstream. To date, little is known about the prognostic impact of eosinophil- and basophil-related markers in CRC patients. In 2018, Wei et al. first reported that circulating hypoeosinophilia and basophilia were associated with worse prognosis in patients with CRC in 569 samples (16). In addition, studies have shown that tumor eosinophilia and basophilia infiltrations contribute to predicting the survival of cancer patients (1719). To our knowledge, the eosinophil-to-neutrophil ratio (20, 21), eosinophil-to-basophil ratio (22), monocyte-to-eosinophil ratio (23, 24), eosinophil-to-lymphocyte ratio (25), and basophil-to-lymphocyte ratio (26) have mainly been reported in inflammatory diseases. However, the relationship between these inflammatory markers and the prognosis of patients with stage II and stage III CRC is indistinct. Additionally, compared to a single type marker, a combination of them might be more valuable and could provide more accurate information for prognosis.

Thus, our study comprehensively analyzed and integrated these inflammatory markers based on eosinophil- and basophil-related markers by least absolute shrinkage and selection operator (LASSO)–Cox regression analysis and constructed a prognostic scoring system named risk score (RS) in CRC. We evaluated the prognostic value of RS for stage II and stage III CRC patients. Additionally, the development and validation of a nomogram for personalized survival prediction might contribute to formulating a better clinical treatment plan.

2 Patients and methods

2.1 Patients and study design

A total of 4,144 primary stage II and III CRC patients undergoing surgical resection followed by pathological diagnosis at the Harbin Medical University Cancer Hospital between January 2007 and December 2013 were enrolled in this study. We excluded patients with neoadjuvant chemotherapy or other radiotherapy/chemotherapy before surgery (n = 2), a postoperative survival time of less than 30 days, or a follow-up time of less than 12 months (n = 156). Finally, 3,986 CRC patients were included in further analyses. Then, they were randomly assigned to the training set and testing set in a 3:2 split ratio. The training set consisting of 2,391 patients was used to train our model, while the testing set consisting of the remaining 1,595 patients was used to evaluate the performance and generalizability of the model after it had been trained. The flowchart of patient screening is shown in Figure 1. This study complied with the standards of the Declaration of Helsinki.

Figure 1

2.2 Data collection

Patients’ clinicopathological features and routine blood results were retrieved from the medical records. The clinicopathological features include sex, age, tumor location, gross appearance, differentiation degree, histological type, tumor stage, cancer nodes, perineural invasion, neoplastic thrombosis, postoperative chemotherapy, and postoperative radiotherapy. Blood routine tests were based on a single blood sample of each patient, which was measured by an autoanalyzer (Sysmex XE-2100, Kobe, Japan). Follow-up information was obtained retrospectively through electronic medical records and telephone interviews. The last time of follow-up was in January 2019. The overall survival (OS) was defined as the time from surgery to death from any cause or the last follow-up visit.

2.3 Inflammation biomarkers

Data on peripheral white blood cell (W; 109/L), neutrophil (N; 109/L), lymphocyte (L; 109/L), monocyte (M; 109/L), eosinophil (E; 109/L), basophil (B; 109/L), and platelet (P; 109/L) were extracted from the results of the first blood routine tests (limit to 30 days before surgery). When the absolute numbers of basophils and eosinophils were 0 (109/L), we added 0.0001 (109/L) to enable them to calculate ratios. Inflammatory factors based on eosinophils and basophils, including EBR, NER, MER, LER, WER, PER, NBR, MBR, LBR, WBR, and PBR were calculated as follows: EBR = E/B, NER = N/E, MER = M/E, LER = L/E, WER = W/E, PER = P/E, NBR = N/B, MBR = M/B, LBR = L/B, WBR = W/B, and PBR = P/B.

2.4 Development of an inflammatory prognostic scoring system

The optimal cutoff values of the above 11 inflammatory markers for predicting the overall survival of patients with CRC were identified by X-tile 3.6.1 software (27) (Yale University, New Haven, CT, USA), which were then classified as categorical variables according to the cutoff value. The variables below and above the cutoff values were scored 0 and 1, respectively. Pearson’s correlation method was used to calculate correlation coefficients for inflammatory biomarkers. Considering the existence of multicollinearity among them, the LASSO-Cox regression analysis with 10-fold cross-validation was performed to select the optimal prognostic factors using the R package “glmnet” and “survminer”. The inflammatory biomarkers with non-zero coefficients were incorporated to construct the novel RS, which was calculated as follows:

Here, n represents the number of inflammatory markers, Score is the score of each inflammatory marker, Coef is the coefficient of LASSO-Cox regression analysis, and risk score (RS) represents a weighted sum of the prognosis score of each marker. RS was also divided into two groups (low and high) by X-tile software in the training set. A time-dependent receiver operating characteristic (ROC) curve was performed to evaluate the predictive value of RS for OS in CRC patients by the R package “timeROC”.

2.5 Statistical analysis

Continuous data were shown as median (interquartile range [IQR]), whereas categorical variables were reported as numbers and percentages. The Wilcoxon rank sum test and Pearson’s chi-squared test were used to compare the clinicopathological characteristics of the training and testing cohorts. The Kaplan–Meier and log-rank tests were utilized to generate the survival curves and compare the survival differences among the groups. The univariate and multivariate Cox analyses were used to estimate the association between indicators and OS, and the results were presented as hazard ratio (HR) and 95% confidence interval (CI). The subgroup analysis was conducted, stratified by sex, age, tumor location, gross appearance, differentiation degree, histological type, tumor stage, perineural invasion, postoperative chemotherapy, and postoperative radiotherapy in the training set.

A nomogram was established by the independent prognostic factors according to multivariate Cox analysis in the training set through the package “rms” in R software. The performance of the nomogram was evaluated by the concordance index (C-index) and time-dependent ROC curve. The C-index >0.5 indicates that the model could discriminate the outcome. The closer the value of the C-index approached 1.0, the higher the prognostic accuracy. Finally, 1,000 bootstrap resamples were performed for internal validation. The performance of nomograms was explored also through calibration curves and decision curve analysis (DCA).

Statistical analyses were conducted with SPSS 23.0 software (SPSS, Inc., Chicago, IL, USA) and R Studio version 3.6.3. All statistical tests were two-sided, and a p-value<0.05 was considered statistically significant.

3 Results

3.1 Baseline characteristics of patients

A total of 3,986 patients with stage II and III CRC were included in this study; the population was randomly divided into the training set (2,391 patients) and the testing set (1,595 patients). The longest follow-up time was nearly 144 months. Detailed baseline characteristics of each set are described in Table 1. In the training set, the median follow-up time was 73 months (IQR, 56–95.5 months), and the median OS time was 68 months (IQR, 44–92 months). In the testing set, the median follow-up time was 71 months (IQR, 56–94 months), and the median OS time was 67 months (IQR, 45–92 months). There were no significant differences between the training set and the testing set (p > 0.05), indicating that the division of data was balanced (Table 1).

Table 1

CharacteristicOverall, n = 3,986Training set, n = 2,391Testing set, n = 1,595p-Value1
Time, median (IQR)68.00 (45.00, 92.00)68.00 (44.00, 92.00)67.00 (45.00, 92.00)>0.90
Sex, n (%)0.09
 Female2,361 (59%)1,442 (60%)919 (58%)
 Male1,625 (41%)949 (40%)676 (42%)
Age, n (%)0.03
 <602,173 (55%)1,270 (53%)903 (57%)
 ≥601,813 (45%)1,121 (47%)692 (43%)
Location, n (%)>0.90
 Colon cancer1,895 (48%)1,136 (48%)759 (48%)
 Rectal cancer2,091 (52%)1,255 (52%)836 (52%)
Gross appearance, n (%)0.20
 Bulge2,654 (67%)1,575 (66%)1,079 (68%)
 Infiltration or ulcer1,332 (33%)816 (34%)516 (32%)
Differentiation degree, n (%)0.13
 Poor532 (13%)335 (14%)197 (12%)
 Moderate or well3,454 (87%)2,056 (86%)1,398 (88%)
Histological type, n (%)0.30
 Adenocarcinoma3,022 (76%)1,798 (75%)1,224 (77%)
 Mucinous adenocarcinoma or signet ring cell cancer964 (24%)593 (25%)371 (23%)
Tumor stage, n (%)0.40
 Stage II2,348 (59%)1,396 (58%)952 (60%)
 Stage III1,638 (41%)995 (42%)643 (40%)
Cancer nodes, n (%)0.90
 No3,723 (93%)2,232 (93%)1,491 (93%)
 Yes263 (6.6%)159 (6.6%)104 (6.5%)
Perineural invasion, n (%)0.08
 No3,688 (93%)2,198 (92%)1,490 (93%)
 Yes298 (7.5%)193 (8.1%)105 (6.6%)
Neoplastic thrombosis, n (%)0.40
 No3,857 (97%)2,309 (97%)1,548 (97%)
 Yes129 (3.2%)82 (3.4%)47 (2.9%)
Postoperative chemotherapy, n (%)0.70
 No2,279 (57%)1,372 (57%)907 (57%)
 Yes1,707 (43%)1,019 (43%)688 (43%)
Postoperative radiotherapy, n (%)0.80
 No3,807 (96%)2,282 (95%)1,525 (96%)
 Yes179 (4.5%)109 (4.6%)70 (4.4%)
WER, median (IQR)53.71 (31.83, 96.97)52.60 (31.98, 96.93)54.31 (31.53, 98.85)>0.90
PER, median (IQR)2,129.46 (1,189.93, 4,087.11)2,101.56 (1,192.09, 4,023.50)2,160.08 (1,180.57, 4,242.35)0.70
LER, median (IQR)16.15 (9.53, 29.17)16.18 (9.56, 28.93)16.10 (9.52, 29.77)0.70
NER, median (IQR)30.99 (17.91, 59.47)30.65 (18.06, 59.78)31.60 (17.85, 58.69)>0.90
MER, median (IQR)3.57 (2.09, 6.58)3.55 (2.07, 6.58)3.59 (2.13, 6.59)0.50
EBR, median (IQR)3.30 (1.83, 6.60)3.36 (1.89, 6.61)3.21 (1.76, 6.51)0.40
WBR, median (IQR)164.11 (108.51, 301.65)166.17 (108.58, 300.83)161.35 (108.12, 302.22)0.80
PBR, median (IQR)6,361.04 (4,192.57, 12,006.08)6,425.99 (4,235.56, 11,969.20)6,283.52 (4,174.27, 12,160.95)0.80
LBR, median (IQR)46.96 (31.85, 84.18)46.78 (31.77, 83.96)47.27 (31.98, 84.56)0.80
NBR, median (IQR)99.00 (61.83, 192.43)100.47 (62.30, 187.65)96.50 (61.41, 200.11)0.60
MBR, median (IQR)11.22 (7.44, 20.70)11.37 (7.37, 20.19)11.13 (7.56, 21.25)0.60

Baseline characteristics of CRC patients in the training set and testing set.

Wilcoxon rank sum test; Pearson’s chi-squared test.

CRC, colorectal cancer; IQR, interquartile range; WER, white blood cell-to-eosinophil ratio; PER, platelet-to-eosinophil ratio; LER, lymphocyte-to-eosinophil ratio; NER, neutrophil-to-eosinophil ratio; MER, monocyte-to-eosinophil ratio; EBR, eosinophil-to-basophil ratio; WBR, white blood cell-to-basophil ratio; PBR, platelet-to-basophil ratio; LBR, lymphocyte-to-basophil ratio; NBR, neutrophil-to-basophil ratio; MBR, monocyte-to-basophil ratio.

1Bold indicates significance (p-value<0.05).

3.2 Optimal cutoff values of inflammatory biomarkers for predicting the overall survival of CRC

We determined the optimal cutoff values of inflammatory biomarkers in the training set by X-tile software (Supplementary Figures 1–3), in which patients were divided into low and high groups. Univariate Cox analyses for OS of inflammatory factors in the training set showed that higher pretreatment EBR (p = 0.01), PBR (p = 0.003), and LBR (p< 0.001) had significantly favorable OS probability than patients in the low groups, whereas others were accompanied by inferior OS (Supplementary Table 1).

3.3 Risk score construction for overall survival

Pearson’s correlation method was used to calculate correlation coefficients for the above 11 inflammatory markers, which showed a high correlation among the inflammatory markers (Figure 2). According to the results of univariate Cox analyses in the training set, inflammatory factors with p< 0.05 were included in the LASSO-Cox regression model; among the 11 candidate inflammatory biomarkers, PER, LER, NER, MER, WBR, PBR, LBR, NBR, and MBR were non-zero coefficients, and the optimal λ value = 0.0028, log (λ) = −5.8721 (Figure 3). The risk score was calculated based on the corresponding coefficient from LASSO, calculated as follows: RS = 0.0969 × PER + 0.0577 × LER + 0.0833 × NER + 0.2716 × MER + 0.0494 × WBR + (−0.2983) × PBR + (−0.3039) × LBR +0.0896 × NBR + 0.0577 × MBR. X-tile 3.6.1 software was also used to determine the optimal cutoff values for RS, which was 0 (Supplementary Figure 4). Patients were separated into the low-risk group (RS ≤ 0) and high-risk group (RS > 0) for further study. The Kaplan–Meier survival curve showed that the OS probability in the low-risk group was significantly higher than that in the high-risk group (log-rank p< 0.0001, Figure 4A). The prognostic accuracy of RS was evaluated by area under the curve (AUC) in the time-dependent ROC, yielding AUC values with 1-, 3-, 5-, and 10-year OS rates of 0.628, 0.587, 0.571, and 0.491, respectively, in the training set (Figure 4B).

Figure 2

Figure 3

Figure 4

3.4 Independent prognostic factors for CRC patients

In the training set, univariate Cox analysis showed that sex, age, tumor location, gross appearance, degree of differentiation, histological type, tumor stage, cancer nodes, perineural invasion, neoplastic thrombosis, postoperative chemotherapy, postoperative radiotherapy, eosinophils, basophils, and RS were associated with the prognosis of OS (all p< 0.1). All these statistically significant factors were then subjected to the multivariate Cox analysis. After adjustment of clinicopathological characteristics, eosinophils (HR: 0.75, 95% CI, 0.62–0.91) and basophils (HR: 0.81, 95% CI, 0.67–0.98) were still significantly associated with the OS of CRC (Supplementary Table 2). RS (high-risk vs. low-risk, HR: 1.70, 95% CI, 1.43–2.03) remained as an independent prognostic factor for poor OS (Table 2), which was further verified in the entire set (n = 3,986, HR: 1.44, 95% CI, 1.24–1.66). However, the result in the testing set was not significant (n = 1,595, Supplementary Table 3).

Table 2

CharacteristicUnivariate analysesMultivariate analyses
HR95% CIp-ValueHR95% CIp-Value
Sex
 Male1.001.00
 Female0.870.75–1.020.080.830.71–0.970.021
Age
 <601.001.00
 ≥601.641.41–1.91<0.0011.671.43–1.96<0.001
Location
 Colon cancer1.001.00
 Rectal cancer1.311.13–1.52<0.0011.231.05–1.440.010
Gross appearance
 Bulge1.001.00
 Infiltration or ulcer1.471.26–1.71<0.0011.341.14–1.56<0.001
Differentiation degree
 Poor1.001.00
 Moderate or well0.600.50–0.73<0.0010.700.58–0.85<0.001
Histological type
 Adenocarcinoma1.001.00
 Mucinous adenocarcinoma or signet ring cell cancer1.231.04–1.460.0151.351.14–1.60<0.001
Tumor stage
 Stage II1.001.00
 Stage III2.271.95–2.64<0.0012.201.86–2.60<0.001
Cancer nodes
 No1.001.00
 Yes2.091.61–2.70<0.0011.220.93–1.610.150
Perineural invasion
 No1.001.00
 Yes1.761.38–2.25<0.0011.441.11–1.860.006
Neoplastic thrombosis
 No1.001.00
 Yes2.151.53–3.00<0.0011.421.00–2.010.052
Postoperative chemotherapy
 No1.001.00
 Yes0.760.65–0.89<0.0010.660.56–0.79<0.001
Postoperative radiotherapy
 No1.001.00
 Yes1.961.47–2.60<0.0011.791.33–2.41<0.001
Risk score
 Low1.001.00
 High1.681.41–2.00<0.0011.701.43–2.03<0.001

Univariate and multivariate Cox analyses of baseline characteristics and risk score on overall survival in stage II and stage III colorectal cancer patients.

All analyses were adjusted for sex, age, tumor location, gross appearance, differentiation degree, histologic type, tumor stage, cancer nodes, perineural invasion, neoplastic thrombosis, postoperative chemotherapy, and postoperative radiotherapy.

3.5 Subgroup analysis

We investigated the prognostic effect of eosinophils, basophils, and RS in different subgroups stratified by sex, age, tumor location, gross appearance, differentiation degree, histological type, tumor stage, perineural invasion, postoperative chemotherapy, and postoperative radiotherapy in the training set. The results of subgroup analysis showed that both eosinophils and basophils were significantly associated with OS in the subgroups of<60 years, bulge, adenocarcinoma, no perineural invasion, no postoperative chemotherapy, and no postoperative radiotherapy (Supplementary Figures 5, 6). RS was still an independent prognostic factor for OS in all subgroups, except the poorly differentiated group, perineural invasion group, and postoperative radiotherapy group in the training set (Supplementary Figure 7).

3.6 The prognostic accuracy of risk score, TNM, their combination, and previously reported markers

We evaluated the prognostic accuracy of RS, TNM staging, their combination, and previously reported markers (neutrophil-to-lymphocyte ratio (NLR), lymphocyte-to-monocyte ratio (LMR), and platelet-to-lymphocyte ratio (PLR)) by AUC in the time-dependent ROC in the training set; the details of the AUC values are listed in Supplementary Table 4. RS had a stable prognostic effect, which tended to be higher than the NLR, LMR, and PLR throughout the follow-up times and could enhance the prognostic effect of the stage by the combination (Supplementary Figure 8).

3.7 Development and validation of the nomogram model for predicting the overall survival

A nomogram based on multivariate regression of the training set was further built; the nomogram was established by the above clinicopathological characteristics and RS (Figure 5A) by assigning points to each variable at the top line and then calculating the total points to predict 3-year OS probability or 5-year OS probability. For the training set, the nomogram’s C-index was 0.692 (95% CI, 0.672–0.712), which was the same as the bootstrapping method used in internal validation (0.692). For the testing set, the C-index was 0.691 (95% CI, 0.667–0.716) and also 0.691 in the internal validation. Moreover, the calibration curves of the nomogram indicated good agreement between the nomogram-predicted probability of 3 years and the actual 3-year OS proportion in the training set and testing set (Figures 5B, C).

Figure 5

Furthermore, the AUC values of the nomogram were higher than tumor stage in both training (3-year AUC: 0.748 vs. 0.680, Supplementary Figure 9A) and testing (3-year AUC: 0.720 vs. 0.648, Supplementary Figure 9B) cohorts for 3-year OS. Finally, DCA for 3-year OS prediction in the training (Supplementary Figure 9C) and testing (Supplementary Figure 9D) sets also showed favorable effects and had better clinical application value than the tumor stage alone.

4 Discussion

In this large retrospective cohort study, we constructed a novel RS by integrating the inflammatory markers selected by the LASSO analysis (PER, LER, NER, MER, WBR, PBR, LBR, NBR, and MBR). The multivariate Cox regression analysis revealed that RS was an independent prognostic factor for OS, and the high-risk group showed a significantly worse outcome in stage II and stage III CRC patients. RS had a stable prognostic ability at different follow-up times and could enhance the prognostic effect of tumor stage by combination. Furthermore, the nomogram constructed by RS and clinicopathological characteristics might be used for personalized prediction and help clinicians identify high-risk patients.

For CRC patients, the tumor stage is widely used for prognostication (4, 5). However, the system ignores other clinical features, which makes it difficult to represent the overall condition of the tumor and is highly heterogeneous. In recent years, studies have reported that molecular genetic markers, such as microsatellite instability and K-ras/BRAF mutation, are also related to the prognosis of CRC (2830). These molecular genetic markers usually require complex and expensive laboratory techniques. Inflammation also plays a critical role in all stages of tumor progression (3133). Multiple researchers have indicated the prognostic value of inflammation-related factors in CRC patients with different stages (3436). The inflammatory process frequently causes changes in numerous hematological parameters, such as peripheral blood cell counts and levels of C-reactive protein and albumin. In comparison, peripheral blood cell counts are easy to measure, inexpensive, and widely available in routine clinical practice.

Several studies have examined the effects of markers, such as NLR, LMR, and PLR, with the results showing that high NLR, low LMR, and high PLR exhibited the worst OS in CRC (1013, 3739). However, these single markers were only based on the ratios of two types of blood cell counts, which might be influenced by various systemic factors and not accurately provide information on the process of inflammation. Therefore, a comprehensive blood biomarker is urgently needed in clinical practice.

Eosinophils and basophils, as rare sets of peripheral blood leukocytes, play important roles in tumors. Eosinophils are becoming recognized as a powerful immune effector and immunomodulator in the tumor microenvironment and have a potential role in tumor treatment (40). Basophils play a key role in various IgE-mediated and IgE-independent allergic inflammation (41). Studies showed that basophils released several angiogenic factors that play a pivotal role in inflammatory and tumor angiogenesis; histamine is released by basophils and has been implicated in CRC (42, 43). The cancer-changed immune cells in the tumor microenvironment have been reported to be closely related to the markers in peripheral blood (44). Eosinophils and basophils can be found not only in the tumor microenvironment but also in the blood. Previous studies have shown that tumor eosinophilia and basophilia infiltrations contribute to improving the survival of cancer patients (17, 45). Our study found that higher levels of circulating eosinophils and basophils in CRC tumors might be associated with better prognosis and survival, which is consistent with previous results. Therefore, we speculated that the combination of them in the blood might have a great predictive significance.

Our study, for the first time, comprehensively analyzed the 11 common inflammatory biomarkers based on eosinophil- and basophil-related markers, such as WER, PER, LER, NER, MER, EBR, WBR, PBR, LBR, NBR, and MBR; all provide certain prognostic values. To avoid the influence of multicollinearity, we performed the LASSO-Cox regression analysis and identified the nine valuable inflammatory markers to construct the RS. We found that RS was a significant independent prognostic factor in the training set. The prognostic value was validated in the entire cohort.

Previous research had demonstrated that 5-year relative survival rates for CRC patients range from more than 90% in stage I to slightly more than 10% in stage IV (4, 46). We could conclude that the prognosis of patients with stage I CRC was excellent, while patients with stage IV CRC have extremely poor prognoses; moreover, there were fewer patients with stage I and stage IV CRC. Therefore, the current study with a large scale and long follow-up time focused on stage II and stage III CRC patients. However, only relying on the tumor stage could not exactly predict the outcome for individual patients. It is necessary to construct novel prognostic markers with good performance. These markers could help quantify the risk of stage II and III CRC patients accurately. We constructed the RS based on eosinophil- and basophil-related markers in a large sample, and patients were classified into high-risk and low-risk groups, which is important for individualized risk stratification and timely intervention to improve prognosis. After the stratification of patients by subgroup, RS was still an independent prognostic factor for OS. According to the results of time-dependent ROC analysis, RS is superior to NLR, LMR, and PLR in prediction, and the combination of RS and tumor staging can improve the staging effect.

Nomogram is a practical graphical tool that is relatively easy to use and can assess the prognosis of individual patients (47). We tried to establish a prognostic nomogram to make it more intuitive and convenient to evaluate the prognosis of OS during clinical practice, and it can help CRC patients with poor prognosis to obtain better suitable treatment in advance. We developed a nomogram by incorporating RS and significant clinicopathological characteristics. The C-index and AUC of the nomogram were both higher than those for the staging system, which indicates a better prediction effect in the training and testing sets. Furthermore, the calibration curves of the 3-year and the 5-year probability of survival also demonstrated that the nomogram had good discrimination and calibration. Compared with staging, the nomogram also had a higher net benefit, which implies better clinical applicability of the nomogram.

Our study has several limitations that should be acknowledged. First, all the patients were selected from a single Chinese hospital; it will be better to validate the predictive accuracy of the model with external multi-center validation. Additionally, we mainly focused on the correlation between markers and OS in CRC; however, the dynamic changes and specificity of inflammatory markers were ignored. Furthermore, our study was a retrospective cohort, which comes with a limitation that some data on clinicopathological characteristics are lacking, such as other immune cells, tumor-infiltration eosinophils and basophils, microsatellite instability, and mismatch repair gene detection.

We integrated the accessible inflammatory markers based on eosinophils and basophils and constructed a novel RS in patients with stage II and III CRC. RS was shown to remain an independent factor for predicting the prognosis of CRC patients. Additionally, the nomogram developed by RS and clinicopathological characteristics might be used for the personalized prediction of CRC patients.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author contributions

LG, YSZ, and YLL conceived and designed the study. LG, CYu, JF, TT, LZ, YL, SM, YPL, YYZ, JX, CJ, DZ, TZ, QF, ST, LL, and CYa performed the data collection. LG, JF, HH, DL, and LZ performed the statistical analysis. LG performed original draft writing, and JF, YZ, and YLL helped to review the manuscript. All authors read and critically revised the manuscript for intellectual content and approved the final manuscript.

Funding

This work was supported by the Innovative Research Group Project of the National Natural Science Foundation of China (no. 81473055 to YSZ).

Conflict of interest

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

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

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

References

  • 1

    SungHFerlayJSiegelRLaversanneMSoerjomataramIJemalAet al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin (2021) 71(3):209–49. doi: 10.3322/caac.21660

  • 2

    SiegelRLMillerKDJemalA. Cancer statistics, 2019. CA Cancer J Clin (2019) 69:734. doi: 10.3322/caac.21551

  • 3

    ChenWZhengRBaadePDZhangSZengHBrayFet al. Cancer statistics in China, 2015. CA Cancer J Clin (2016) 66:115–32. doi: 10.3322/caac.21338

  • 4

    MillerKDNogueiraLMariottoABRowlandJHYabroffKRAlfanoCMet al. Cancer treatment and survivorship statistics, 2019. CA Cancer J Clin (2019) 69:363–85. doi: 10.3322/caac.21565

  • 5

    DeSantisCELinCCMariottoABSiegelRLSteinKDKramerJLet al. Cancer treatment and survivorship statistics, 2014. CA Cancer J Clin (2014) 64:252–71. doi: 10.3322/caac.21235

  • 6

    HariDMLeungAMLeeJHSimMSVuongBChiuCGet al. AJCC cancer staging manual 7th edition criteria for colon cancer: do the complex modifications improve prognostic assessment? J Am Coll Surg (2013) 217:181–90. doi: 10.1016/j.jamcollsurg.2013.04.018

  • 7

    MeiZLiuYLiuCCuiALiangZWangGet al. Tumour-infiltrating inflammation and prognosis in colorectal cancer: systematic review and meta-analysis. Br J Cancer (2014) 110:1595–605. doi: 10.1038/bjc.2014.46

  • 8

    WuJGeXXZhuWZhiQXuMDDuanWet al. Values of applying white blood cell counts in the prognostic evaluation of resectable colorectal cancer. Mol Med Rep (2019) 19:2330–40. doi: 10.3892/mmr.2019.9844

  • 9

    FuJZhuJDuFZhangLLiDHuangHet al. Prognostic inflammatory index based on preoperative peripheral blood for predicting the prognosis of colorectal cancer patients. Cancers (Basel) (2020) 13:3. doi: 10.3390/cancers13010003

  • 10

    ShaulMEFridlenderZG. Tumour-associated neutrophils in patients with cancer. Nat Rev Clin Oncol (2019) 16:601–20. doi: 10.1038/s41571-019-0222-4

  • 11

    LeeKHKimEYYunJSParkYLDoSIChaeSWet al. The prognostic and predictive value of tumor-infiltrating lymphocytes and hematologic parameters in patients with breast cancer. BMC Cancer (2018) 18:938. doi: 10.1186/s12885-018-4832-5

  • 12

    ZhangYYLiWQLiZFGuoXHZhouSKLinAet al. Higher levels of pre-operative peripheral lymphocyte count is a favorable prognostic factor for patients with stage I and II rectal cancer. Front Oncol (2019) 9:960. doi: 10.3389/fonc.2019.00960

  • 13

    ShimuraTShibataMGondaKHayaseSSakamotoWOkayamaHet al. Prognostic impact of preoperative lymphocyte-to-monocyte ratio in patients with colorectal cancer with special reference to myeloid-derived suppressor cells. Fukushima J Med Sci (2018) 64:6472. doi: 10.5387/fms.2018-10

  • 14

    TaoYDingLYangGGQiuJMWangDWangHet al. Predictive impact of the inflammation-based indices in colorectal cancer patients with adjuvant chemotherapy. Cancer Med (2018) 7:2876–86. doi: 10.1002/cam4.1542

  • 15

    LiYJiaHYuWXuYLiXLiQet al. Nomograms for predicting prognostic value of inflammatory biomarkers in colorectal cancer patients after radical resection. Int J Cancer (2016) 139:220–31. doi: 10.1002/ijc.30071

  • 16

    WeiYZhangXWangGZhouYLuoMWangSet al. The impacts of pretreatment circulating eosinophils and basophils on prognosis of stage I–III colorectal cancer. Asia Pac J Clin Oncol (2018) 14:e243–e51. doi: 10.1111/ajco.12871

  • 17

    PrizmentAEVierkantRASmyrkTCTillmansLSLeeJJSriramaraoPet al. Tumor eosinophil infiltration and improved survival of colorectal cancer patients: Iowa women's health study. Mod Pathol (2016) 29:516–27. doi: 10.1038/modpathol.2016.42

  • 18

    De MonteLWormannSBrunettoEHeltaiSMagliacaneGReniMet al. Basophil recruitment into tumor-draining lymph nodes correlates with Th2 inflammation and reduced survival in pancreatic cancer patients. Cancer Res (2016) 76:1792–803. doi: 10.1158/0008-5472.CAN-15-1801-T

  • 19

    RigoniAColomboMPPucilloC. Mast cells, basophils and eosinophils: from allergy to cancer. Semin Immunol (2018) 35:2934. doi: 10.1016/j.smim.2018.02.001

  • 20

    ZhangXYSimpsonJLPowellHYangIAUphamJWReynoldsPNet al. Full blood count parameters for the detection of asthma inflammatory phenotypes. Clin Exp Allergy (2014) 44:1137–45. doi: 10.1111/cea.12345

  • 21

    GaoJChenBWuSWuF. Blood cell for the differentiation of airway inflammatory phenotypes in COPD exacerbations. BMC Pulm Med (2020) 20:50. doi: 10.1186/s12890-020-1086-1

  • 22

    BresciaGBarionUZanottiCCinettoFGiacomelliLMartiniAet al. Blood eosinophil-to-basophil ratio in patients with sinonasal polyps: does it have a clinical role? Ann Allergy Asthma Immunol (2017) 119:223–6. doi: 10.1016/j.anai.2017.06.008

  • 23

    YucelCKeskinMZCakmakOErganiBKoseCCelikOet al. Predictive value of pre-operative inflammation-based prognostic scores (neutrophil-to-lymphocyte ratio, platelet-to-lymphocyte ratio, and monocyte-to-eosinophil ratio) in testicular sperm extraction: a pilot study. Andrology (2017) 5:1100–4. doi: 10.1111/andr.12417

  • 24

    OzbirSDegirmentepeRBAtalayHAAlkanICakirSSOtunctemurAet al. The role of inflammatory parameters (neutrophil-to-lymphocyte ratio, platelet-to-lymphocyte ratio, and monocyte-to-eosinophil ratio) in patients with peyronie's disease. Andrology (2020) 8:348–52. doi: 10.1111/andr.12702

  • 25

    BranickaORogalaBGluckJ. Eosinophil/Neutrophil/Platelet-to-Lymphocyte ratios in various types of immediate hypersensitivity to NSAIDs: a preliminary study. Int Arch Allergy Immunol (2020) 181:774–82. doi: 10.1159/000509116

  • 26

    BresciaGPedruzziBBarionUCinettoFGiacomelliLMartiniAet al. Are neutrophil-, eosinophil-, and basophil-to-lymphocyte ratios useful markers for pinpointing patients at higher risk of recurrent sinonasal polyps? Am J Otolaryngol (2016) 37:339–45. doi: 10.1016/j.amjoto.2016.02.002

  • 27

    CampRLDolled-FilhartMRimmDL. X-Tile: a new bio-informatics tool for biomarker assessment and outcome-based cut-point optimization. Clin Cancer Res (2004) 10:7252–9. doi: 10.1158/1078-0432.CCR-04-0713

  • 28

    PopatSHubnerRHoulstonRS. Systematic review of microsatellite instability and colorectal cancer prognosis. J Clin Oncol (2005) 23:609–18. doi: 10.1200/JCO.2005.01.086

  • 29

    KarapetisCS K-FSJonkerDJO'CallaghanCJTuDTebbuttNCSimesRJet al. K-Ras mutations and benefit from cetuximab in advanced colorectal cancer. N Engl J Med (2008) 359:1757–65. doi: 10.1056/NEJMoa0804385

  • 30

    Tol JNIPuntCJ. BRAF mutation in metastatic colorectal cancer. N Engl J Med (2009) 361:98–9. doi: 10.1056/NEJMc0904160

  • 31

    CoussensLMWerbZ. Inflammation and cancer. Nature (2002) 420:860–7. doi: 10.1038/nature01322

  • 32

    BalkwillFMantovaniA. Inflammation and cancer: back to virchow? Lancet (2001) 357:539–45. doi: 10.1016/S0140-6736(00)04046-0

  • 33

    MantovaniAAllavenaPSicaABalkwillF. Cancer-related inflammation. Nature (2008) 454:436–44. doi: 10.1038/nature07205

  • 34

    ChenJHZhaiETYuanYJWuKMXuJBPengJJet al. Systemic immune-inflammation index for predicting prognosis of colorectal cancer. World J Gastroenterol (2017) 23:6261–72. doi: 10.3748/wjg.v23.i34.6261

  • 35

    SuzukiYOkabayashiKHasegawaHTsurutaMShigetaKKondoTet al. Comparison of preoperative inflammation-based prognostic scores in patients with colorectal cancer. Ann Surg (2018) 267:527–31. doi: 10.1097/SLA.0000000000002115

  • 36

    YamamotoTKawadaKObamaK. Inflammation-related biomarkers for the prediction of prognosis in colorectal cancer patients. Int J Mol Sci (2021) 22(15):s8002. doi: 10.3390/ijms22158002

  • 37

    KimJHLeeJYKimHKLeeJWJungSGJungKet al. Prognostic significance of the neutrophil-to-lymphocyte ratio and platelet-to-lymphocyte ratio in patients with stage III and IV colorectal cancer. World J Gastroenterol (2017) 23:505–15. doi: 10.3748/wjg.v23.i3.505

  • 38

    ChanJCYChanDLDiakosCIEngelAPavlakisNGillAet al. The lymphocyte-to-Monocyte ratio is a superior predictor of overall survival in comparison to established biomarkers of resectable colorectal cancer. Ann Surg (2017) 265:539–46. doi: 10.1097/sla.0000000000001743

  • 39

    JiaWYuanLNiHXuBZhaoP. Prognostic value of platelet-to-Lymphocyte ratio, neutrophil-to-Lymphocyte ratio, and lymphocyte-to-White blood cell ratio in colorectal cancer patients who received neoadjuvant chemotherapy. Technol Cancer Res Treat (2021) 20:15330338211034291. doi: 10.1177/15330338211034291

  • 40

    ReichmanHKaro-AtarDMunitzA. Emerging roles for eosinophils in the tumor microenvironment. Trends Cancer (2016) 2:664–75. doi: 10.1016/j.trecan.2016.10.002

  • 41

    MiyakeKKarasuyamaH. Emerging roles of basophils in allergic inflammation. Allergol Int (2017) 66:382–91. doi: 10.1016/j.alit.2017.04.007

  • 42

    LiuQLuoDCaiSLiQLiX. Circulating basophil count as a prognostic marker of tumor aggressiveness and survival outcomes in colorectal cancer. Clin Transl Med (2020) 9:6. doi: 10.1186/s40169-019-0255-4

  • 43

    PotoRGambardellaARMaroneGSchroederJTMatteiFSchiavoniGet al. Basophils from allergy to cancer. Front Immunol (2022) 13:1056838. doi: 10.3389/fimmu.2022.1056838

  • 44

    GuoGWangYZhouYQuanQZhangYWangHet al. Immune cell concentrations among the primary tumor microenvironment in colorectal cancer patients predicted by clinicopathologic characteristics and blood indexes. J Immunother Cancer (2019) 7:179. doi: 10.1186/s40425-019-0656-3

  • 45

    SektiogluIMCarreteroRBulbucNBaldTTutingTRudenskyAYet al. Basophils promote tumor rejection via chemotaxis and infiltration of CD8+ T cells. Cancer Res (2017) 77:291302. doi: 10.1158/0008-5472.CAN-16-0993

  • 46

    BrennerHKloorMPoxCJL. Colorectal cancer. Lancet (2014) 383:1490–502. doi: 10.1016/s0140-6736(13)61649-9

  • 47

    KattanMWScardinoPT. Evidence for the usefulness of nomograms. Nat Clin Pract Urol (2007) 4:638–9. doi: 10.1038/ncpuro0968

Summary

Keywords

colorectal cancer, prognosis, biomarker, inflammation, eosinophils, basophils, risk score

Citation

Gao L, Yuan C, Fu J, Tian T, Huang H, Zhang L, Li D, Liu Y, Meng S, Liu Y, Zhang Y, Xu J, Jia C, Zhang D, Zheng T, Fu Q, Tan S, Lan L, Yang C, Zhao Y and Liu Y (2023) Prognostic scoring system based on eosinophil- and basophil-related markers for predicting the prognosis of patients with stage II and stage III colorectal cancer: a retrospective cohort study. Front. Oncol. 13:1182944. doi: 10.3389/fonc.2023.1182944

Received

31 March 2023

Accepted

26 June 2023

Published

14 July 2023

Volume

13 - 2023

Edited by

Maria Julia Lamberti, National University of Río Cuarto, Argentina

Reviewed by

Eleonora Lai, University Hospital and University of Cagliari, Italy; Matias Caverzan, National University of Río Cuarto, Argentina

Updates

Copyright

*Correspondence: Yashuang Zhao, ; Yanlong Liu,

†These authors have contributed equally to this work and share first authorship

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics