3.1 Effect of Genetically Predicted Immunophenotypes On CRC
To investigate the causal effect of immunophenotypes on CRC, we conducted a two-sample MR analysis using five assays. The IVW method served as the primary analytical basis, and the results, presented in Figure 3, indicate that 29 immunophenotypes were linked to colorectal cancer risk. To minimize false-positive results, we employed the FDR method and Bonferroni correction for multiple testing. After correcting the IVW-P values by FDR, we found that 26 immunophenotypes were associated with colorectal cancer risk (P-FDR<0.05). These included six immunophenotypes of the B cell panel, 5 immunophenotypes of the TBNK panel, 7 immunophenotypes of the Treg panel, 4 immunophenotypes of the maturation stages of the T cell panel, 4 immunophenotypes in the Treg panel, 4 immunophenotypes in the maturation stages of the T cell panel, one in the monocyte and cDCs panel, and two in the myeloid cell panel (Supplementary Figure 1, 2, 3, and 4).
Among the Immunophenotypes with a P-FDR<0.05, we found 16 immunophenotypes that suggest similar risk ratios (OR>1). These immunophenotypes show an increased risk of developing colorectal cancer with an increase in RC, AC, or MFI for that specific immunophenotypes. (Such as B cell panel: BAFF-R on naive-mature B cell, BAFF-R on B cell, Treg panel: CD28+ CD45RA+ CD8br AC, CD3 on CD39+ CD4+, CD3 on CD28+ CD45RA+ CD8br, CD3 on CD28- CD8br, CD25 on secreting Treg , TBNK panel: Lymphocyte %leukocyte, T cell %leukocyte, CD3 on HLA DR+ T cell , CD3 on HLA DR+ CD8br, Maturation stages of the T cell panel: CD3 on EM CD4+, CD3 on CM CD8br, CD4 on CD45RA+ CD4+, cDCs panel: CD62L- HLA DR++ monocyte AC, Myeloid cell panel: CD33br HLA DR+ CD14- AC).
There were 10 immunophenotypes with OR less than 1, indicating a relative reduction in the risk of colorectal cancer with an increase in RC, AC, or MFI for that immunophenotype. (Such as B cell panel: Naive-mature B cell AC, CD24+ CD27+ %B cell, Transitional %lymphocyte, CD20 on IgD+ CD24-, Treg panel: CD39+ CD4+ %T cell, CD25++ CD8br %T cell, TBNK panel: DN (CD4-CD8-) %leukocyte, Maturation stages of the T cell panel: HVEM on EM CD8br, Monocyte panel: CCR2 on monocyte, Myeloid cell panel: Im MDSC %CD33dim HLA DR- CD66b-). In the TBNK panel, we found that four immunophenotypes were positively causally associated with the risk of colorectal cancer (Supplementary Figure 5).
To further minimize the risk of false-positive P-values, we utilized the most rigorous Bonferroni correction method. After correction, P-Bonferroni was found to be <1.72×10-3 (0.05/29), only two immunophenotypes were significantly associated with colorectal cancer (Figure 4A): lymphocyte % leukocyte belonging to the TBNK panel, and CD3 on CM CD8br belonging to the Maturation stages of the T cell panel.
For the Lymphocyte %leukocyte in TBNK panel, the OR was estimated to be 1.0011 (95% CI=1.0005-1.0017, P=4.44×10-4) using IVW method. Additional analytical methods yielded the following results: MR-Egger (OR=1.0010, 95% CI 1.0002-10017, P=0.028), Weighted Median (OR=1.0010, 95% CI 1.0000-1.0020, P=0.042), Simple mode (OR=1.0008, 95% CI 0.9991-1.0025, P=0.36), Weighted mode (OR=1.0010, 95%CI 1.0002-1.0019, P=0.025). The MR-PRESSO global P-value and Pleiotropy P value were both >0.05 in the test of horizontal pleiotropy and the test of heterogeneity was also >0.05. In the Maturation stages of the T cell panel, CD3 on CM CD8br held the OR of 1.0013 (95% CI 1.0005-1.0021, P=9.54×10-4) by IVW method. Results from other methods are listed as follows: MR-Egger (OR=1.0006, 95% CI 0.9991-1.0021, P=0.464), Weighted Median (OR=1.0011, 95% CI 0.9999-1.0023, P=0.065), Simple mode (OR=1.0013, 95% CI 0.9995-1.0032, P=0.178), Weighted mode (OR = 1.0011, 95% CI 0.9999-1.0022, P=0.091). Similarly in the test of horizontal multiplicity, both MR-PRESSO global P-value and Pleiotropy P-value were >0.05 and the test of heterogeneity was also >0.05. Further validation was provided by the stability observed in the scatter plots, funnel plots and leave-one-out analysis, which proved the robustness of the results (Figure 5).
3.2 Effect of Genetically Predicted CRC On Immunophenotypes
Before adjustment, CRC was observed to induce increase in 2 immunophenotypes (Figure 4B). Linear relationships, funnels, and leave-one-out plots for two immune cells in reverse Mendelian randomization analysis, as shown in Figure (Supplementary Figure 6). However, after Bonferroni adjustment, none of these p values remained statistically significant.
3.3 Results of Case-Control Study
3.3.1. Baseline Characteristic
The study included 100 patients diagnosed with CRC and 100 members of the healthy control group (Figure 6). There were no significant differences between the two cohorts in terms of age, sex, and smoking versus alcohol consumption (Supplementary Table 1). However, patients with colorectal cancer had higher absolute counts of lymphocytes, CD3+T cells, CD4+T cells, CD4+ and CD8+T cells, and NK cells, as well as higher relative counts of lymphocytes and CD4+T cells compared to the healthy control group (P<0.05). However, the counts of CD8+T cell and B cells were relatively lower in colorectal cancer patients (Figure 7A and B).
3.3.2 Single-factorial Conditional Logistic Regression Analysis of Immune Cells and Risk of Colorectal Cancer Development
A single-factorial conditional logistic regression model was constructed with the occurrence of colorectal cancer as the dependent variable, stratified by absolute and relative count quartiles of peripheral blood immune cells in colorectal cancer cases and healthy controls (Supplementary Table 2 and 3).
Based on the Supplementary Table 2, the risk of colorectal cancer was 4.317 times higher for study subjects in the Lymphocyte 3rd quartile group compared to the 1st quartile group (OR=4.317, 95% CI: 1.812-10.282), The highest quartile group had a risk 9.25 times higher (OR=9.25, 95% CI: 3.744-22.852).
The study found that subjects in CD4+T cell quartile 2 had a 2.043-fold increased risk of colorectal cancer compared to the reference subgroup (OR=3.043, 95% CI: 1.328-6.977). Quartile 3 had a 1.826-fold increased (OR=2.826, 95% CI: 1.225-6.52), while the highest quartile group had a 3.211-fold increased risk (OR=4.211, 95% CI: 1.817-9.758). Study participants in the highest CD8+T cell quartile group nevertheless had a 55.6% lower risk of colorectal cancer (OR=0.444, 95% CI: 0.199-0.989). Additionally, the risk of colorectal cancer was 63.2% lower in the 3rd quartile of B-cell subgroup compared to the 1st quartile (OR=0.368, 95% CI 0.162-0.835), and even further reduced by 65.1% in the highest quartile group (OR=0.349, 95% CI: 0.156-0.781).
The Supplementary Table 3 results, stratified by absolute count group quartiles, indicate that only study subjects in CD4+T cell quartile 3 had a 1.941-fold increased risk of colorectal cancer compared to quartile 1 (OR=1.941, 95% CI: 1.081-3.485). None of the other absolute count metrics were statistically significant.
3.3.3 Multi-factorial Conditional Logistic Regression of Immune Cells and Risk of Colorectal Cancer Development
Variables that were statistically significant (P<0.05) in the single-factorial analysis were included in the multi-factorial conditional logistic regression model for further analysis (Supplementary Table 4 and Figure 7C). The relative counts of lymphocytes and CD4+T cells were positively associated with the risk of colorectal cancer, while the relative counts of B cells were negatively associated with the risk of colorectal cancer.
After adjusting for age, sex, smoking, and alcohol consumption, the risk of colorectal cancer was significantly higher for study participants in quartile 2 of the Lymphocyte relative count (OR=4.113, 95% CI: 1.415-11.95), quartile 3 (OR=13.79, 95% CI: 4.027-47.215), and the highest quartile group (OR=39.48, 95% CI: 9.384-166.097) compared to those in quartile 1. Subjects in the second quartile of CD4+T cell relative counts had a 4.949 times greater risk of colorectal cancer than those in the first quartile (OR=4.949, 95% CI: 1.639-14). Similarly, individuals in the 3rd quartile had an increased risk of 8.973 times (OR=8.973, 95% CI: 2.531-31.812), while those in the highest quartile had a 9.367-fold risk (OR=9.367, 95% CI: 2.495-35.159). The risk of colorectal cancer was significantly lower in the 2nd, 3rd, and highest quartile groups of relative B-cell count compared to the 1st quartile group (OR=0.251, 95% CI: 0.085-0.741, OR=0.141, 95% CI: 0.048-0.415, OR=0.183, 95% CI: 0.062-0.541). In the CD4+T cell absolute count subgroups, the risk of colorectal cancer was reduced by 71.4% in study subjects in the 2nd quartile group compared to the 1st quartile group (OR=0.286, 95% CI: 0.087-0.935) and by 77.8% in the highest quartile group (OR=0.222, 95% CI: 0.053-0.93). However, the remaining absolute count metrics did not show statistical significance (Figure 7D).
3.3.4 Clinical Baseline Characterization and Prognostic Analysis of Immune Cells in Colorectal Cancer
We analyzed the immune cell counts of CRC cases along with their baseline characteristics and prognosis to investigate the potential association between immune cells and CRC (Table 1). The study found that increased lymphocyte counts were linked to risk factors for metastasis, recurrence, and death in colorectal cancer. Additionally, high values of CD4+T cell relative counts were significantly associated with elevated Clinical N stage (Supplementary Table 5, 6, 7, and 8).
To investigate the clinicopathological characteristics of colorectal cancer and the effect of immune cells on DFS, and to analyze the association. In this study, we grouped the cohort of colorectal cancer patients by median, categorized the absolute and relative counts of immune cells into high and low groups, and plotted the KM curves, and the results are shown in Supplementary Table 9 and Supplementary Figure 7. In the Kaplan-Meier curves of 34 colorectal cancer patients who developed recurrence, metastasis, and death, age, Tumor differentiation, and absolute counts of CD3+T cells and CD4+T cells as well as relative counts of lymphocytes in the immune cells were statistically significantly associated with DFS (P<0.05) (Figure 8). Specifically, lower age, increased CD3+T cell, CD4+T cell absolute counts, and relative lymphocyte counts resulted in superior DFS.
Subsequently, we conducted Cox regression model analysis on cases of colorectal cancer and found that age, tumor location, absolute CD3+T cell count, absolute CD4+T cell count, and relative lymphocyte count were associated with DFS after treatment. These indexes were included in the multivariate COX regression model analysis, which showed that only left hemi-colon carcinoma was an independent factor influencing DFS after treatment (HR=0.303, 95% CI: 0.103-0.894) (Supplementary Table 10).
Table 1. Baseline clinicopathologic features of colorectal cancer cases
Characters
|
|
Age, year
|
64.48±10.48
|
Sex Male: Female, n
|
58:42
|
Smoking Yes: No, n
|
28:72
|
Alcohol Yes: No, n
|
30:70
|
Tumor location Right-side: Left-side: Rectum, n
|
46:25:29
|
Tumor differentiation Low: Middle-low: Middle, n
|
14:31:55
|
Pathological T stage 1: 2: 3: 4, n
|
10:6:76:8
|
Clinical N stage 0: 1: 2: 3, n
|
48:30:17:5
|
Microsatellite Instability MSI-L: MSI-H, n
|
96:4
|
Recurrence/metastasis/death Yes: No, n
|
34:66
|
DFS, month
|
12.44±4.28
|