Clinical characteristics in full cohort
A total of 351 surgically resected FFPE primary CRC samples were included in the study. In the whole cohort, 39.0% (137 out of 351) were infected with schistosoma (Figure 1A). The clinical and pathologic features of the cohort are summarized in Table 1. In the whole cohort, age of patients at diagnosis ranged from 33 to 91 years (median, 69 years) and were predominantly male (60.2%, 212 out of 351). By anatomic site, 27% tumors were in the rectum, 33% in left colon and 40% in right colon. Lymph node metastasis were observed in 40% of patients and 46% of patients were at late stage disease. While patients without lymph node metastasis were 60%. On the basis of the AJCC Staging Manual (seventh edition), there were very few highly differentiated cases in the follow-up data. Thus, highly differentiated and moderately differentiated cases were classified as "well differentiation", and classified poorly differentiated cases as "poor differentiation". 76% cases were well differentiated, and 24% were poorly differentiated. As shown in Table 1, lymphovascular invasion, perineura invasion, cancer node and tumor budding were prone to appear in patients with stage III-IV tumors or patients with lymph node metastasis. More poorly differentiated tumors and deeper tumor invasion depth were also mostly observed in patients with late tumor stage or patients with lymph node metastasis. The distribution trend of other clinicopathologic features, such as colonic perforation, ulceration, histological type were similar within different subgroups.
Survival analysis
The median follow-up times was 62.4 (1.25-134.4) months. During the follow up, there was 41.6% (146 out of 351) patients died. Mean and median times to overall survival were 62.54 and 62.85, respectively.
To investigate the association between schistosomiasis and clinical outcomes, we conducted Kaplan-Meier analysis according to schistosoma infection status. Result demonstrated that schistosoma infection was significantly associated with poor survival in total colorectal cancer patients (median survival time: 80.82 for CRC-S set; 119.20 for CRC-NS set. P=0.0277) (Figure 1B).
Further analysis was conducted to explore the effect of schistosoma infection on CRC patients with similar stage tumors. In stage Ⅰ-Ⅱ set (N=192), a Kaplan-Meier (K-M) curve was plotted and found that schistosoma infection (40%) was uncorrelated with survival (P=0.5018) (Figure 2A). Nevertheless, in stage Ⅲ-Ⅳ set (N=159), K-M analysis showed significant correlation between schistosoma infection and OS (P=0.0260) (Figure 2B).
In patients with lymph node metastasis (N=144), schistosoma infection was observed in 39% (56 out of 144) CRC patients and associated with poor survival (P=0.0249) (Figure 2C). In contrast, there was no statistically significant difference observed in OS between CRC-S and CRC-NS patients without lymph node metastasis (P=0.4005) (Figure 2D).
Univariate and multivariate analysis
The Cox proportional hazards model was used to determine factors that may influence OS of CRC patients. In the whole cohort, by univariate analysis and multivariate analysis (Table 2), gender (P=0.003), TNM stage (P<0.001), schistosomasis (P=0.025), lymphovascular invasion(P=0.030) and cancer node (P<0.001) were significantly independent predictors. Schistosomiasis was statistically significant associated with decreasing OS.
In late-stage (Ⅲ-Ⅳ) CRC patients(Table 2), gender (P=0.030), pathological T stage (P=0.12), tumor differentiation (P=0.016), schistosoma infection (P=0.008) and cancer node (P=0.004) were significantly independent prognostic factor for OS. While in early-stage (Ⅰ-Ⅱ), cancer node (P=0.007) was the only independent prognostic factor for OS in multivariate analysis.
In patients with lymph node metastasis (Table 2), gender (P=0.026), pathological T stage (P=0.025), schistosoma infection (P =0.023) and cancer node (P=0.003) were independently prognostic factors. In patients without lymph node metastasis (Table 2), TNM stage (P<0.001) and tumor budding (P=0.014) but not schistosoma infection were associated with OS in multivariate analysis. These results further proved that schistosoma infection may have different effects on CRC patients` clinical outcomes, especially for patients with stage Ⅲ-Ⅳ tumor and patients with lymph node metastasis.
Association of schistosomiasis with clinicopathological features
The relationship between schistosomiasis and clinicopathological features was shown in Table 3. Patients with schistosomiasis were significantly older than the patients without schistosomiasis (median age: 74.0 years vs 64.0 years, p<0.001). Clinical stage of patients with schistosomiasis and without schistosomiasis were similar (p=0.816). In the total cohort, the male/female ratio was also higher in CRC-S set (1.67 vs 1.43). Besides, in patients with lymph node metastasis, there were significant association between male sex and female sex (p<0.001). There were no significant differences of other clinicopathological characteristics between CRC-NS and CRC-S set.
In order to further investigate the effect of schistosomiasis on particular CRC population, we divided the whole cohort into different groups according to their clinical stage or the state of lymph node metastasis and further subgrouped into CRC-S and CRC-NS set based on schistosomiasis. Except age, there were no correlation between other clinicopathological features and schistosomiasis when compared between CRC-NS and CRC-S sets in different groups(Table 3).