Clinical characterization of cases
These original cases were selected from diagnosed CRC patients, and these cases with missing the key parameters or clinical markers were excluded for next step general and specific level analysis. The process and characteristic analyses were demonstrated by flow chart (Figure.1). Characterization of these cases was summarized, including corresponding clinical features, 5-year survival rate and tumor distribution of primary site (Table.1). First of all, the qualified the patients were divided into four age-groups (Table.1), including below 50-years-old group (124 cases, 12.3%), 50-60-years-old group (273 cases, 27.1%), 61~70-years-old group (326 cases, 32.4%) and over 70-years-old group (282 cases, 28.0%). The majority of cases (882 cases, 87.7%, p<0.001) were significantly distributed in the range of over 50 -years-old (Figure.S1A). The 5-year survival rate of over 60-year-old patients (median 55.5 months, 44%) was significantly lower than below 60-year-old patients [median 62 months, 52%, log rank p=0.0004] (Figure.S1B). It suggests that population of over 60-year older patients have higher risk to develop into CRC and lower survival rate. There is not significantly differences in the five-year survival rates between gender groups.
Characterization of tumor primary sites
Totally, 737 cases (73.3% incident rate) were distributed in the terminal tumor site (median survival period = 57 months, 5-year survival rate = 48.4%), while the other cases (n=269, 26.7%) were allocated in others tumor sites (median survival period = 59 months, and 5-year survival rate =48.3%) (Table.1). Altogether, 11 primary tumor sites were identified as rectum (59.3%), sigmoid (13.9%), right hemicolon (7.6%), colon (6.3%), ileocecal (4.5%), descending colon (1.1%), hepatic flexure colon (0.8%), ascending colon (2.4%), transverse colon (1.7%), junction of the rectum and sigmoid (0.6%), and left hemicolon (1.9%). Rectum and sigmoid were classified as the terminal tumor sites 17, and determined with higher incidence rate compared to other tumor sites (Figure 3 A&B).
Combined gender bias with primary tumor sites analysis, male patients (58%) had significantly higher incident rate than female patients (42%), especially in the terminal tumor sites. The distribution ratio of male patients (61%, 86 cases; 57%, 341 cases) is higher than female patients (39%, 54 cases; 43%, 256 cases) in sigmoid cancer and rectal cancer, respectively (Figure 3C&D). Combined age with primary tumor sites analysis, 74.8% of 50-70-years-old cases were distributed in terminal tumor sites, 25.2% of patients allocated in others tumor sites. The diagnosis rate of early stage terminal colorectal cancer (41.1%) was significantly higher than that of other primary tumor sites (31.2%) (Table 2), while the advance stage of terminal colorectal cancer (58.9%) was significantly lower than that other primary tumor sites (68.8%). The OS and five-year survival rates of terminal tumor sites (sigmoid and rectum, 48.4%) were similar with other primary tumor sites (48.3%).
In this study, the percentage of tumor non-infiltrated cases was significantly higher than tumor infiltrated cases. The OS of tumor non-infiltrated cases (median 54 months, 40%, log rank p=2e-15, Figure.S4A) was significantly lower than tumor-infiltrated cases (median 67 months, 57%, p<0.01). The majority of tumor non-infiltrated cases were similarly distributed in the terminal tumor sites (441/656, 67%) and other tumor sites (163/229, 71%), but 98% of tubular adenocarcinoma cases were identified as non-infiltrated tumor (Figure.S4B).
Tumor stage, lymph node metastasis, adenocarcinoma subtypes and survival outcomes
Patients with advanced stage (median 55 months, 43%), high degree of differentiation (median 55 months, 43%), lymph node metastasis (median 54 months, 42%) and severe lymph node metastasis (over number 5, median 44 months, 27%) have significantly shorter OS and lower 5-years survival rate than those with early stage (median 64 months, 55%, log rank p=7e-10, Figure 2A), lower degree of differentiation grade (median 67 months, 56% , log rank p=0.024, Figure 2B), without lymph node metastasis (median 61 months, 51% , Figure 2C, log rank p=0.007) and mild lymph node metastasis (below number 5, median 60 months, 49%, log rank p =1e-05, Figure 2D). More interesting, the tumor stage (r= -0.18, p = 5.77e-09), lymph node metastasis (r= -0.14, p = 1.47e-05) and the number of lymph node metastasis (r= -0.17, p = 3.79e-08) were significantly negative correlated with patient’s OS (Figure S2A). The lymph node metastasis (r= 0.67, p < 2.2e-16) and the number of lymph node metastasis (r= 0.55, p < 2.2e-16) showed significantly positive correlation with tumor stage. The lymph node metastasis (r= 0.82, p < 2.2e-16) was significantly positive correlated with the number of lymph node metastasis (Figure S2A). The patient’s 5-year survival rate of different lymph node metastasis number dropped from 52% to 26% (Figure S2B), and the patients with lymph node metastasis number over number 5 have significantly worse survival outcome than those below 5 (Figure S2C, p = 0.016).
Through tumor subtype analysis, 988 patients were identified with pathohistological feature and classified as subtype of tubular adenocarcinoma, mucinous adenocarcinoma, adenocarcinoma, adenocarcinoma with necrosis and mixed type (without tubular adenocarcinoma, contained more than two other subtypes), respectively (Figure.S1C). Among four subtypes of adenocarcinoma, the survival rate of tubular adenocarcinoma patients (median 38 months, 5%) was significantly lower than others histological tumors (median 60.5 months, 50%, log rank p =2e-15, Figure.4A). Furthermore, the majority of tubular adenocarcinoma cases were mainly distributed in the terminal tumor site (78%), and especially occurred in the sites of sigmoid (14/58, 24%) and rectum (31/58, 53%) (Figure.4B, p = 2e-15).
Correlation of serum markers with the CRC primary tumor
Multiple serum indicators of CRC patients were tested and summarized (Table. S2). The correlation analysis revealed that different patterns of these physiological indicators were correlated with primary tumor sites (Figure 3C, Figure.S3A), tumor stage (Figure.S3B) and overall survival time. First of all, the serum level of CA19-9 was positive correlated with CRC tumor sites (r=0.88, p=2.2e-16), Overall Survival time (r=0.11, p value=0.05454) and negative correlated with tumor stage (r=-0.13, p=0.0214) (Figure 3C). The serum level of CEA was significantly positive correlated with CRC tumor stage (r=0.19, p =3.53-e06) and very weak negative correlated with CRC tumor site (r=-0.06, p=0.16) and overall survival time (r=-0.06, p value = 0.1546) (Figure 3C). In the meanwhile, the other two markers CRP (r= -0.21, p = 2.65e-7) and ALB (r= -0.19, p=3.54e-6) were significantly negative correlated with CRC tumor sites (Figure S3C), and showed weak correlation with CRC tumor stage (CRP, r=0.06, p =0.1544; ALB, r= -0.05, p =0.23) and Overall survival time (CRP, r=-0.08, p value = 0.03882; ALB, r=0.00159, p value=0.969) (Figure. S3C). The others physiological indicators show little correlation with CRC stage and tumor sites. The advanced stage tumor patients, with higher abnormal level of CEA (>5 mg/L, 41%), CRP (>8 mg/L, 33%) and ALB (< 40g/L, 57%), were mainly distributed in terminal tumor sites (Figure.4D). The others seven physiological indicators (Table.S1), including leukocyte, blood platelet, neutrophil, lymphocyte, monocyte, did not show significantly difference between normal and tumor patients. In addition, the patients (99%) with higher level of CA19-9 (CA19-9-2) were mainly distributed in the terminal tumor sites of colorectal cancer, and only small portion allocated in other primary tumor sites (Table.1, Figure.4E). However, there was no significant difference in overall survival time of the two different CA19-9 levels patient groups. Combined with the tumor stage analysis, lower CA19-9 level patients (73%) were more likely to be associated with advanced stage than higher CA19-9 level patients (58%). The serum level of CA19-9 was weak positive correlated with tumor non-infiltrating state (Figure.S4C, r=0.07, p=1.2e-6). Furthermore, the majority of patients diagnosed as tubular adenocarcinoma (81%) were associated with the higher level of CA19-9 (Figure.4F). Combined several markers’ analysis with primary tumor sites (Figure.4D), CA19-9 may work as a promising biomarker for terminal bowel cancer and tubular adenocarcinoma of CRC.