A three-phase trans-ethnic study reveals B7-H3 expression is a significant and independent biomarker associated with colon cancer overall survival
Introduction
Colon cancer is a common malignant tumor of the digestive system, which is the third and second most commonly diagnosed cancer in men and women worldwide, respectively, with 1.15 million new cases and 576,858 deaths in 2020 (1). Even though the therapeutic effect of colon cancer has greatly improved with the development of medical technology, the prognosis of colon cancer is still poor (2). Immunotherapy, which boosts the immune system works to find and attack cancer cells, envisions a novel approach to cancer therapy (3). The immune system consists of a complex set of molecular interactions that are regulated by immune checkpoints and are named differentiated clusters based on their order of discovery (4,5). Immune checkpoint pathways are often activated to inhibit nascent anti-tumor immune responses to tumor cells, and immune checkpoint therapy enhances the immune activity against tumors by blocking or stimulating these pathways (6,7), including the most widely studied and well-recognized pathways of cytotoxic T lymphocyte-associated molecule-4 (CTLA-4), programmed cell death receptor-1 (PD-1), and programmed cell death ligand 1 (PD-L1). Inhibitory immune checkpoints (e.g., PD-L1, PD-1, CTLA-4, and B7-H3) bind with their associated ligands to induce an inhibitory response and prevent the cascade of stimulatory immune checkpoint signals that activate T cell activity (8,9). Therefore, those immune checkpoints are consequentially associated with tumor prognosis and response to treatment.
Nevertheless, the association between B7-H3 (B7 Homolog 3) and colon cancer survival still remains unclear. Also known as CD276, B7-H3 is a costimulatory molecule belonging to the B7-CD28 family. It is a membrane protein encoded on chromosome 15 (10,11), and is widely expressed in heart, thymus, prostate, testis, uterus, placenta, spleen, liver, pancreas, small intestine, and colon in both normal and tumor tissues (12-14). An array of studies have indicated that highly expressed B7-H3 promotes tumor progression and metastasis, as well as being associated with poor prognosis in various of cancers, including glioma (15), hepatocellular carcinoma (16), lung cancer (17), breast cancer (18), osteosarcoma (19), cutaneous melanoma (20), and pancreatic cancer (21). However, the results of existing prognostic studies of B7-H3 expression and colon cancer (22-37) have been inconsistent and have presented conflicting conclusions; thus, we believe there is an urgent need to elaborate the effect of B7-H3 on colon cancer survival. Besides, various technical bottlenecks of previous studies need to be addressed: (I) there have been very few European studies (22,23,29), and majority of these studies were conducted among Asian populations (24-28,30-37); (II) no study so far has focused on the non-linear effects of B7-H3 on colon cancer survival, with previous studies simply having tested its linear effect; (III) there has been no independent validation of B7-H3, the association has consistently been evaluated in a single population; and (IV) there have been no trans-ethnic population studies of B7-H3 thus far.
Hence, we performed a 3-phase designed and trans-ethnic study to test and confirm the prognostic effect of B7-H3 expression on colon cancer survival using 3 independent cohorts comprising European and Asian populations, followed by a series of Kaplan-Meier (K-M) analysis, bioinformatics analysis of tumor immune microenvironment (TIME), and immune checkpoint prognostic prediction analysis, along with sensitivity analysis. We present the following article in accordance with the REMARK reporting checklist (available at https://dx.doi.org/10.21037/jgo-21-821) (38).
Methods
Study populations
We utilized European colon cancer patients from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) cohorts, as well as Asian colon cancer patients recruited from Fujian Medical University Union Hospital (UNION) to evaluate the association between gene expression level of B7-H3 and overall survival (OS) of colon cancer (Table 1).
Table 1
Variable | TCGA (n=433) | GEO (n=811) | UNION (n=179) |
---|---|---|---|
Age (years), mean ± SD | 66.33±12.83 | 66.28±13.31 | 55.87±15.22 |
Age group, n (%) | |||
<65 years | 169 (39.03) | 325 (40.07) | 126 (70.39) |
≥65 years | 264 (60.97) | 485 (59.80) | 53 (29.61) |
Unknown | 18 | 1 | 0 |
Gender, n (%) | |||
Female | 200 (46.19) | 370 (45.62) | 89 (49.72) |
Male | 233 (53.81) | 441 (54.38) | 90 (50.28) |
Height (cm), mean ± SD | 168.52±12.32 | – | 163.09±7.82 |
Weight (kg), mean ± SD | 81.34±20.72 | – | 59.63±10.57 |
T stage, n (%) | |||
T0/Tis | 1 (0.23) | 4 (0.69) | 0 (0.00) |
T1 | 11 (2.54) | 12 (2.07) | 0 (0.00) |
T2 | 75 (17.32) | 48 (8.29) | 0 (0.00) |
T3 | 296 (68.36) | 376 (64.94) | 34 (18.99) |
T4 | 50 (11.55) | 119 (20.55) | 145 (81.01) |
Unknown | 0 | 252 | 0 |
N stage, n (%) | |||
N0 | 254 (58.66) | 311 (53.71) | 32 (17.88) |
N1 | 102 (23.56) | 136 (23.49) | 50 (27.93) |
N2 | 77 (17.78) | 100 (17.27) | 97 (54.19) |
N3 | 0 (0.00) | 6 (1.04) | 0 (0.00) |
NX | 0 (0.00) | 6 (1.04) | 0 (0.00) |
Unknown | 0 | 252 | 0 |
M stage, n (%) | |||
M0 | 320 (73.90) | 496 (85.66) | 0 (0.00) |
M1 | 61 (14.09) | 61 (10.54) | 179 (100.00) |
MX | 45 (10.39) | 2 (0.35) | 0 (0.00) |
Unknown | 7 | 252 | 0 |
Clinical stage, n (%) | |||
I | 73 (16.86) | 65 (8.01) | 0 (0.00) |
II | 165 (38.11) | 341 (42.05) | 0 (0.00) |
III | 123 (28.41) | 285 (35.14) | 0 (0.00) |
IV | 61 (14.09) | 116 (14.30) | 179 (100.00) |
Unknown | 11 | 4 | 0 |
Adenocarcinoma, n (%) | |||
Yes | 367 (84.76) | 579 (100.00) | 135 (75.42) |
No | 66 (15.24) | 0 (0.00) | 44 (24.58) |
Unknown | 0 | 232 | 0 |
Tumor location, n (%) | |||
Left | 135 (31.18) | 228 (39.38) | 111 (62.01) |
Right | 204 (47.11) | 351 (60.62) | 68 (37.99) |
Unknown | 94 | 232 | 0 |
Complete tumor excision, n (%) | |||
Yes | 53 (12.24) | – | 57 (31.84) |
No | 13 (3.00) | – | 122 (68.16) |
Unknown | 367 | – | 0 |
Carcinoma cell embolus, n (%) | |||
Negative | 226 (52.19) | – | 84 (46.93) |
Positive | 164 (37.88) | – | 95 (53.07) |
Unknown | 43 | 811 | 0 |
Tumor size, n (%) | |||
<5 cm | – | – | 66 (36.87) |
≥5 cm | – | – | 113 (63.13) |
Unknown | 433 | 811 | 0 |
Postoperative therapy, n (%) | |||
Yes | 145 (33.48) | 240 (41.45) | 130 (72.63) |
No | 235 (54.27) | 323 (55.79) | 48 (26.82) |
Unknown | 53 | 248 | 1 |
CEA (median, Q1−Q3) | 3.00 (1.80–7.52) | – | 8.00 (2.70–33.4) |
CEA group, n (%) | |||
Normal | 178 (41.11) | – | 73 (40.78) |
Elevated | 97 (22.40) | – | 106 (59.22) |
Unknown | 158 | 811 | 0 |
CA199, n (%) | 136.37±277.48 | ||
Normal | – | – | 106 (59.22) |
Elevated | – | – | 73 (40.78) |
Unknown | 433 | 811 | 0 |
Median survival months | |||
Median (95% CI) | 93 (66–NA) | 135 (106–NA) | 16 (11–20) |
Deaths, n (%) | 95 (21.94) | 287 (35.39) | 137 (76.54) |
NA, not available; T, tumor; N, node; M, metastasis; CEA, carcinoembryonic antigen; TCGA, The Cancer Genome Atlas; GEO, Gene Expression Omnibus; UNION, Fujian Medical University Union Hospital.
TCGA
Gene expression profiles (platform: Illumina HiSeq 2000 RNA Sequencing; San Diego, CA, USA) and clinical data of colon cancer patients were obtained from TCGA (https://portal.gdc.cancer.gov/) database in July 2021, including 471 tumor tissues and 41 adjacent-normal tissues. A total of 433 participants (95 deceased and 338 alive) with complete clinical and OS time were reserved for subsequent association analysis. The gene expression level was measured by fragments per kilobase of transcript per million fragments (FPKM) value and log2 transformed before analysis. Unqualified probes were excluded if they meet any of the quality control (QC) criteria: (I) high missing rates (>30%); or (II) coefficient of variance (CV) <5% (Figure S1). The data used in this study comply with the requirements of TCGA official published data and are publicly available.
GEO
A total of 811 colon cancer patients with complete clinical and OS information, as well as gene expression data, were acquired from the GEO (https://www.ncbi.nlm.nih.gov/geo/) database, including GSE39582, GSE17536, and GSE17537. Gene expression data was profiled by GPL570 Affymetrix Human Genome U133 Plus 2.0 Array platform (Affymetrix, Santa Clara, CA, USA). The raw intensity values were background corrected, log2 transformed, and then quantile normalized. Next, the expression values were derived by robust multi-array average (RMA) method, an algorithm exclusively designed to create an expression matrix from Affymetrix platform. Redundant probes were collapsed, and then annotated to human gene symbols prior to analysis. We applied the same QC criteria for the GEO data.
UNION
Colon cancer patients in the UNION cohort were recruited from January 2010 to January 2018. Patients were admitted for surgical treatment of colorectal cancer. Routine preoperative examinations were performed to exclude patients with organic lesions of the heart, lung, liver, and others. Additionally, patients (I) who did not receive surgery but only a simple ileostomy; (II) who had intestinal obstruction, intestinal perforation, intestinal bleeding, or required emergency surgical resection; (III) who had simultaneous polygenic carcinoma of the colon, heterogeneous polygenic carcinoma of the colon, and familial adenomatous polyposis; (IV) who presented other malignant tumors or died due to postoperative complications, were further excluded. Finally, 179 patients were included in the study. The expression level of B7-H3 in tumor cells was quantified by immunohistochemistry (IHC). Briefly, Slides (4-µm thick consecutive paraffin sections) from the blocks with the highest tumor content for each sample were used for IHC staining, immersed in xylene, rehydrated through graded concentrations of ethanol followed by phosphate-buffered saline (PBS) buffer, and deionized water for 5 min each. Slides were then heated to 100 °C for 20 min in a pH 9 Tris-based solution. All slides were incubated with the primary antibodies for 60 min at 37 °C for 1 h (dilutions: Human B7-H3 Antibody 1:300, R&D Systems, Minneapolis, MN, USA) and were then washed. A secondary antibody followed by incubation with rabbit anti-goat IgG and an avidin-biotin complex (Boster, Wuhan, China) was added for 30 min and the slides were again washed. The sections were processed with the universal SP Elivision-plus kit (Maixin Bio, Fuzhou, China). Finally, the sections were counterstained with hematoxylin (39). The percentages of B7-H3 in 3 representative high-power fields of individual samples were analyzed for intensity of B7-H3 membranous expression and were scored as 0 (<5% of B7-H3), 1 (5–25%), 2 (>25–50%), or 3 (>50%). Individual samples were evaluated by at least 2 pathologists in a blinded manner, and those samples with inconsistent scores were further discussed and decided. A score of 0 or 1 with B7-H3 on IHC was regarded as low B7-H3 expression group, and 2 or 3 staining as high B7-H3 expression group. This study was approved by the institutional review board of Fujian Medical University Union Hospital (No. 2018YF024-01). All participants or their surrogates provided their written informed consent. All procedures performed in this study involving human participants were in accordance with the Declaration of Helsinki (as revised in 2013).
Statistical analysis
Three-phase study of B7-H3 expression and colon cancer OS
The analysis workflow is presented in Figure 1. We adopted a 3-phase trans-ethnic study to evaluate the association between expression of B7-H3 and colon cancer OS. In the discovery phase, we used a restricted cubic spline (RCS) regression model by setting the median expression of B7-H3 as a reference point to investigate the potential non-linear effect of B7-H3 on OS using a European population from the TCGA cohort. The RCS are a smoothly joined sum of polynomial functions that do not assume linearity of the relationship between variable and the response (40). This technique provides great flexibility for fitting data and modeling complex relationships between survival outcome and the variable of interest, while adjusting for other covariates. Furthermore, an RCS regression model can avoid arbitrary categorization of the continuous measures and permit identification of the risk function inflexion point. Then, the B7-H3 expression was stratified by knot information into different groups (41). Finally, the association between B7-H3 and survival was tested using a Cox proportional hazards model adjusted for covariates. These covariates were significant variables screened out by step forward regression model with P value of entry ≤0.05 and P value of remove >0.05, among a series of demographic and clinical variables. As a result, age, gender, and clinical stage were adjusted in the Cox proportional hazards model (Tables S1,S2). Proportional hazard assumption was confirmed based on Schoen’s method (42). The association strength was measured by hazard ratio (HR) and 95% confidence interval (CI).
In the validation phase I, we used the RCS regression model again to evaluate the non-linear effect of B7-H3 on colon cancer survival using a European population from the GEO cohort. The Cox proportional hazards model was utilized to test the association between B7-H3 and survival, adjusted for the same covariates aforementioned in the discovery phase.
In the validation phase II, we directly tested the grade level of B7-H3 using Cox proportional hazards regression model using an Asian population from the UNION cohort. Since all participants were in M1 stage, the model was adjusted for age, gender, tumor (T) stage, and node (N) stage. Finally, meta-analyses using both fixed-effects and random-effects models were applied in our study to combine results from 3 phases by R package meta (https://cran.r-project.org/web/packages/meta/index.html).
Supplementary analyses of B7-H3
Furthermore, we performed a series of supplementary analyses. First, K-M survival curves were drawn to represent the survival difference between patients with different expression levels (43). The cutoff values of B7-H3 expression were defined based on using R package survMisc (https://cran.r-project.org/web/packages/survMisc/index.html), which used log-rank test to find the optimal cut-off value, in TCGA and GEO cohorts, respectively. According to these cutoff values, these participants were divided into a high and low expression group, their survival differences were further compared using Cox proportional hazards regression model adjusted for covariates.
Second, bioinformatics analysis of B7-H3 including immune infiltration analysis, somatic copy number alteration (SCNA) analysis, and correlation analysis were performed using TIMER2.0, which is a freely available web server to the research community (44). The TIMER2.0 resource uses RNA-seq expression profiles to detect and quantify the situation of immune cell infiltration in tumor tissues, so as to determine the relationship between tumor and immune cells in TCGA colon cancer database (45). (I) In immune infiltration analysis, the scatterplots were displayed, showing the association between B7-H3 expression and immune infiltration level in colon cancer tumors using Spearman correlation analysis. (II) In SCNA analysis, SCNAs were categorized into different groups, including deep deletion (−2), arm-level deletion (−1), diploid/normal (0), arm-level gain (+1), and high amplification (+2). Box plots were presented to show the distributions of each immune subset (B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells) at each copy number status in colon cancer. The infiltration level for each category was compared with the normal using a 2-sided Wilcoxon rank-sum test. (III) In correlation analysis, we further evaluated the association between B7-H3 and 27 well-recognized immune checkpoint genes (46).
Finally, we proposed an immune checkpoint prognostic score of colon cancer using overall participants from TCGA and GEO cohorts to distinguish patients at high risk of mortality. These well-recognized immune checkpoint genes associated with colon cancer survival were screened out by forward stepwise regression model adjusted for age, gender, clinical stage, study site and B7-H3, with P value of entry ≤0.05 and P value of remove >0.05. Patients were then categorized into high and low risk groups by median value, and K-M survival curves between 2 groups were tested. Sensitivity analyses were also performed stratified by age, gender, and clinical stage. The time-dependent receiver operating characteristic (ROC) curve and area under the ROC curve (AUC) were used to measure the prediction ability (47). Statistical analyses were performed using R version 3.4.4 (The R Foundation of Statistical Computing, Vienna, Austria).
Results
After quality control, 3 datasets were composed of with 433 participants in TCGA cohort, 811 participants in the GEO cohort, and 179 participants in the UNION cohort for further association analysis. Demographic and clinical characteristics of these participants were described as mean ± SD for continuous variables and number (%) for categorical variables in Table 1.
Clinical experimental studies of B7-H3 and colon cancer
First, the expression level of B7-H3 in colon cancer and adjacent normal tissues were evaluated by differential analysis of gene expression data in TCGA and GEO cohort. Compared with the adjacent normal tissue, tumor tissue had high-expressed level of B7-H3 (P=1.18×10−21). Besides, there was significant association between B7-H3 expression and clinicopathological factors (age: P=0.029; clinical stage: P=0.002), which were shown in boxplot (Figure S2).
B7-H3 was significantly and robustly associated colon cancer OS
In the discovery phase, we observed the apparent non-linear association pattern between B7-H3 and OS time of colon cancer patients from cubic spline regression model using European population in TCGA (Figure 2). This suggested that the inflection points were 3.7, 4.1, and 4.6. The covariates-adjusted Cox proportional hazards regression model indicated B7-H3 was a significant risk factor of colon cancer survival (HRTCGA =4.60, 95% CI: 2.15 to 9.83, P=8.37×10−05). In the validation phase I, there was an approximate linear association pattern between B7-H3 and colon cancer survival derived from cubic spline regression model using European population in GEO. By Cox model, B7-H3 was still significantly correlated with a worse prognosis (HRGEO =1.47, 95% CI: 1.12 to 1.94, P=0.0056). In validation phase II, we again confirmed that B7-H3 was a significant prognostic factor of colon cancer survival using an Asian population in UNION (HRUNION =1.63, 95% CI: 1.36 to 1.95, P=7.91×10−08). Finally, meta-analysis of 3 phases ensured the significant association regardless of a fixed-effects model (HRfixed =1.65, 95% CI: 1.42 to 1.91, P=2.85×10−11) or random-effects model (HRrandom =1.85, 95% CI: 1.28 to 2.66, P=0.0010) (Table 2).
Table 2
Variable | Unadjusted model | Adjusted model | |||||
---|---|---|---|---|---|---|---|
HR | 95% CI | P value | HR | 95% CI | P value | ||
TCGA | 3.240 | 1.514–6.933 | 0.0024 | 4.598 | 2.150–9.833 | 8.37E–05 | |
GEO | 1.332 | 1.016–1.746 | 0.0380 | 1.472 | 1.119–1.937 | 0.0056 | |
UNION | 1.553 | 1.317–1.831 | 1.7E–07 | 1.631 | 1.364–1.950 | 7.91E–08 | |
Meta | |||||||
Fixed-effects model | 1.529 | 1.331–1.756 | 1.91E–09 | 1.646 | 1.421–1.907 | 2.85E–11 | |
Random-effects model | 1.585 | 1.203–2.087 | 0.0011 | 1.846 | 1.281–2.661 | 0.0010 | |
Test of heterogeneity | 0.0919 | 0.0217 |
HR, hazard ratio; CI, confidence interval; TCGA, The Cancer Genome Atlas; GEO, Gene Expression Omnibus; UNION, Fujian Medical University Union Hospital.
K-M survival analysis
The optimal cutoff value of B7-H3 which corresponded to the highest log-rank test score was 4.317 and 5.864 in TCGA and GEO, respectively. Therefore, participants in TCGA and GEO were divided into high and low expression groups. In UNION, participants with the grade levels (0 and 1) of B7-H3 expression were categorized into low expression group and the others were in high expression group. The survival curves were significantly separated between the 2 groups across 3 datasets (Figure 3), indicating that high expression of B7-H3 had harm effect on colon cancer survival (HRTCGA =1.63, 95% CI: 1.07 to 2.49, P=2.22×10−02; HRGEO =1.75, 95% CI: 1.26 to 2.42, P=7.62×10−04; HRUNION =2.36, 95% CI: 1.64 to 3.40, P=4.21×10−06).
Bioinformatics analysis
Cancer cells within the TIME and neighboring tumor-associated noncancerous cells play an important role in tumor biology. The immune infiltration analysis indicated that B7-H3 expression had a significant negative correlation with immune infiltration level for tumor purity and B cells, but a significant positive correlation with that for CD4+ T cells, macrophages, neutrophils, and dendritic cells (Figure 4). Besides, the infiltration levels were significantly distributed among different groups of SCNAs of B7-H3 (Figure 5), including B cell (P<0.05), CD8+ T cell (P<0.001), neutrophil (P<0.01), and dendritic cell (P<0.001). In addition, B7-H3 was significantly and weakly correlated with 6 out of 27 well-recognized immune checkpoint genes (Figure S3), including PVRL2 (ρ=0.47, P=1.90×10−05), PVR (ρ =0.43, P=8.31×10−5), IGSF11 (ρ=0.29, P=9.74×10−03), CEACAM1 (ρ=0.35, P=1.90×10−03), PTDSS1 (ρ=0.43, P=1.04×10−04), and HMGB1 (ρ=0.49, P=6.61×10−05).
Immune checkpoint prognostic score of colon cancer
We performed stepwise regression model to screen other immune checkpoint genes using Cox model adjusted for covariates and B7-H3 expression and identified 7 more genes associated with colon cancer survival (Table S3). In this multivariable Cox model, B7-H3 was still significant (HR =1.47, 95% CI: 1.07 to 2.02, P=0.0184), indicating it was an independent risk factor of colon cancer prognosis. Then, we proposed an immune checkpoint prognostic score of colon cancer, which was a weighted linear combination of 8 significant genes. The weight of each gene was defined as the corresponding estimated log hazard ratio [ln(HR)] derived from the multivariable Cox model. The immune checkpoint prognostic score could distinguish patients at high risk of mortality in all subpopulations stratified by age (Figure S4), gender (Figure S5), and clinical stages (Figure S6). These 8 immune checkpoint genes slightly improved the prognostic prediction ability for 3-year survival (2.38%) and 5-year survival (1.63%), respectively. Hence, the prognostic model of colon cancer achieved an acceptable accuracy (AUC3-year =0.74 and AUC3-year =0.72) (Figure S7).
Discussion
In recent years, previous studies of B7-H3 have mainly focused on the its expression and role in different types of cancers (48). Elevated expression of B7H3 has been found in serum of patients with malignant tumors, which is related to clinical stage and progression of cancer (49,50), since B7-H3 appears to be correlated with different proteins that contribute to tumor migration, invasion and angiogenesis (51,52). Besides, B7-H3 blocking in combination with chemotherapy is a promising treatment option based on preclinical animal studies (53).
Currently, some previous studies have described a significant correlation between tumor expression of B7-H3 and prognosis (14,54), whereas, majority of them reported that high tumor B7-H3 expression was associated with more advanced disease (31), increased risk of recurrence (35) or shorter survival time (55), which was consistent with our results. Effective therapeutic strategies against colon cancer are scarce. Therefore, it is imperative to seek a promising therapeutic target for anti-tumor drugs and also important to develop a prognostic tool to identify patients at high risk of mortality who require more attention and aggressive treatment. By leveraging publicly available transcriptional data for European colon cancers in TCGA and GEO cohorts and our in-door data for Asian colon cancers in the UNION cohort, we performed a three-phase study, which perhaps was the first trans-ethnic prognostic study of B7-H3 expression on this population level. We first tested the prognostic effect of B7-H3 on survival in the European population, and then continuously validated its effect using another 2 independent European and Asian populations. Highly expressed B7-H3 was a significant and independent harmful factor to colon cancer survival, even after adjusting for effects of other immune checkpoint genes, correlating with these genes, and joint involvement in TIME. We also proposed an immune checkpoint prognostic risk score which significantly stratified patients into high and low risk groups, regardless of their characteristics.
The mechanism of B7-H3 expression affecting colon cancer survival still remains to be fully investigated. B7-H3 expression is often induced on immune cells (56,57), which can mediate chemotherapeutic resistance and sensitivity to survival of patients (58,59). It was also shown that B7-H3 promotes VEGFA expression and angiogenesis, which dependents on the NF-κB pathway in colorectal cancer (CRC), and CRC cell recruits regulatory T cells to promote chemoresistance via NF-κB signaling pathway (60). Besides, B7-H3 is also regulated by Th1, IL-4, IFN-gamma and TNF-alpha, which is an important cancer-promoting inflammatory molecule and significantly increases the release of soluble B7-H3 in colon cancer cell lines (61,62). Meanwhile, the B7-H3 pathway has a dual role in contributing to the regulation of innate immune responses. On the one hand, B7-H3 co-stimulates innate immunity by augmenting proinflammatory cytokines release from lipopolysaccharide-stimulated monocytes or macrophages in both a Toll like receptor 4- and 2-dependent manner (49,50,63). On the other hand, B7-H3 provides an additional mechanism allowing neuroblastoma cells, which protect them from natural killer cell-mediated lysis, to escape the control of immune response which may be associated with the expression of B7-H3 and multiple possibly related molecules, such as programmed cell death protein 1 (PD-1),cytotoxic T-lymphocyte-associated protein 4 (CTAL-4) or the presence of different types of cells associated with immune responses (11,64). Nevertheless, the role of B7-H3 in controlling the innate immunity is quite complex and requires more elucidation (65). Emerging evidence indicates that immune checkpoint blockade combined with radiotherapy acquires considerable success in multiple tumors (66,67). On account of its important role in immune evasion and cancer progression, B7-H3 also possesses the capability to be a potential novel therapeutic target for colon cancer patients (68). Besides, B7-H3 also plays a non-immunological role in promoting tumor invasion and progression (69,70). It was observed that the 3' UTR of B7-H3 transcript contains a target site for miR29 microRNA, and there is an inverse correlation between the expression of B7-H3 protein and miR29 levels, suggesting regulation of expression of this gene product by miR29. Experimental studies in mice showed that induction of colitis and water avoidance stress affect levels of Mir29a and Mir29b and intestinal permeability in wild-type mice (71). Despite inconsistent results and conflicting conclusions of existing prognostic studies of B7-H3 (22,24-30,32-37), our study, with the largest sample size, confirmed that it was a significant and independent biomarker to colon cancer survival by several trans-ethnical validations.
Our study has several strengths. First, this was perhaps the largest prognostic study of B7-H3 in colon cancer patients on a population level. A 3-phase design was applied to control false positives and ensure the robustness of our results. The significant association observed in the discovery phase was further replicated in the validation phases I and II. Therefore, the probability of reporting the false positive result was strictly controlled under 0.0125% (5%×5%×5%). Second, though population heterogeneity exists even among those of the same ethnicity, taking advantage of the RCS model, we successfully elucidated the nonlinear effect of B7-H3 expression on colon cancer survival in the TCGA European population and its approximate linear effect in GEO European population. Third, by trans-ethnic validation using an Asian population, we confirmed the robust association between B7-H3 expression and colon cancer survival again, meaning that this potential therapeutic target can be uniformly applied to all colon cancer patients worldwide. Fourth, by adjusting for other well-recognized immune checkpoint genes, B7-H3 was still a significant biomarker to colon cancer survival, indicating it was an independent prognostic factor. Finally, we proposed an immune checkpoint prognostic risk score of colon cancer using 8 genes, which can identify patients at high risk of mortality regardless of age, gender, and clinical stages, possible providing an effective tool aiding in clinical interventions.
We also acknowledge some limitations. First, the censoring rate is high in the TCGA and GEO cohorts, which resulted in a low statistical power in time-to-event data analysis. Regardless, we still observed significant results, indicating the association between B7-H3 expression and colon cancer survival was quite conservative. Second, several important clinical variables (e.g., status of complete tumor excision and tumor size) were not available in the TCGA and GEO cohorts, therefore a limited number of covariates were adjusted in the Cox proportional hazards regression model. We envision a better model adjusted for more clinical variables to achieve more precise estimates of B7-H3 effect and more accurate prognostic prediction ability in future studies. Finally, functional studies of B7-H3 are warranted to confirm its biological effect.
Conclusions
In summary, in the largest prognostic study of B7-H3 in colon cancer patients on population levels, we utilized a 3-phase design to confirm that B7-H3 was a significant and independent biomarker of colon cancer prognosis. The proposed immune checkpoint prognostic risk score has the capability to identify colon cancers with high risk of mortality.
Acknowledgements
The authors express great thanks to patients and investigators who participated in TCGA and GEO projects for sharing the raw data.
Funding: This work was supported by the National Natural Science Foundation of China (81902378 to P Chi); Construction Project of Fujian Province Minimally Invasive Medical Center (2017-171 to P Chi); Startup Fund for Scientific Research, Fujian Medical University (2018XQ2025 to Y Gao); Joint Funds for the Innovation of Science and Technology of Fujian Province (2017Y9038, 2017Y9104) and Young and Middle-aged Backbone Training Project in the Health System of Fujian Province (2019-ZQN-45).
Footnote
Reporting Checklist: The authors have completed the REMARK reporting checklist. Available at https://dx.doi.org/10.21037/jgo-21-821
Data Sharing Statement: Available at https://dx.doi.org/10.21037/jgo-21-821
Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://dx.doi.org/10.21037/jgo-21-821). 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. All procedures performed in this study involving human participants were in accordance with the Declaration of Helsinki (as revised in 2013). This study was approved by the institutional review board of Fujian Medical University Union Hospital (No. 2018YF024-01). All participants or their surrogates provided their written informed consent.
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
- 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]
- Xie YH, Chen YX, Fang JY. Comprehensive review of targeted therapy for colorectal cancer. Signal Transduct Target Ther 2020;5:22. [Crossref] [PubMed]
- Waldman AD, Fritz JM, Lenardo MJ. A guide to cancer immunotherapy: from T cell basic science to clinical practice. Nat Rev Immunol 2020;20:651-68. [Crossref] [PubMed]
- Momtaz P, Postow MA. Immunologic checkpoints in cancer therapy: focus on the programmed death-1 (PD-1) receptor pathway. Pharmgenomics Pers Med 2014;7:357-65. [PubMed]
- Pardoll DM. The blockade of immune checkpoints in cancer immunotherapy. Nat Rev Cancer 2012;12:252-64. [Crossref] [PubMed]
- Marin-Acevedo JA, Dholaria B, Soyano AE, et al. Next generation of immune checkpoint therapy in cancer: new developments and challenges. J Hematol Oncol 2018;11:39. [Crossref] [PubMed]
- Dougan M, Dranoff G, Dougan SK. Cancer immunotherapy: beyond checkpoint blockade. Annu Rev Cancer Biol 2019;3:55-75. [Crossref]
- Rotte A. Combination of CTLA-4 and PD-1 blockers for treatment of cancer. J Exp Clin Cancer Res 2019;38:255. [Crossref] [PubMed]
- He X, Xu C. Immune checkpoint signaling and cancer immunotherapy. Cell Res 2020;30:660-9. [Crossref] [PubMed]
- Janakiram M, Shah UA, Liu W, et al. The third group of the B7-CD28 immune checkpoint family: HHLA2, TMIGD2, B7x, and B7-H3. Immunol Rev 2017;276:26-39. [Crossref] [PubMed]
- Chapoval AI, Ni J, Lau JS, et al. B7-H3: a costimulatory molecule for T cell activation and IFN-γ production. Nat Immunol 2001;2:269-74. [Crossref] [PubMed]
- Liu F, Zhang T, Zou S, et al. B7-H3 promotes cell migration and invasion through the Jak2/Stat3/MMP9 signaling pathway in colorectal cancer. Mol Med Rep 2015;12:5455-60. [Crossref] [PubMed]
- Du H, Hirabayashi K, Ahn S, et al. Antitumor Responses in the Absence of Toxicity in Solid Tumors by Targeting B7-H3 via Chimeric Antigen Receptor T Cells. Cancer Cell 2019;35:221-237.e8. [Crossref] [PubMed]
- Yang M, Tang X, Zhang Z, et al. Tandem CAR-T cells targeting CD70 and B7-H3 exhibit potent preclinical activity against multiple solid tumors. Theranostics 2020;10:7622-34. [Crossref] [PubMed]
- Souweidane MM, Kramer K, Pandit-Taskar N, et al. Convection-enhanced delivery for diffuse intrinsic pontine glioma: a single-centre, dose-escalation, phase 1 trial. Lancet Oncol 2018;19:1040-50. [Crossref] [PubMed]
- Kang FB, Wang L, Jia HC, et al. B7-H3 promotes aggression and invasion of hepatocellular carcinoma by targeting epithelial-to-mesenchymal transition via JAK2/STAT3/Slug signaling pathway. Cancer Cell Int 2015;15:45. [Crossref] [PubMed]
- Yang S, Wei W, Zhao Q. B7-H3, a checkpoint molecule, as a target for cancer immunotherapy. Int J Biol Sci 2020;16:1767-73. [Crossref] [PubMed]
- Arigami T, Narita N, Mizuno R, et al. B7-h3 ligand expression by primary breast cancer and associated with regional nodal metastasis. Ann Surg 2010;252:1044-51. [Crossref] [PubMed]
- Majzner RG, Theruvath JL, Nellan A, et al. CAR T Cells Targeting B7-H3, a Pan-Cancer Antigen, Demonstrate Potent Preclinical Activity Against Pediatric Solid Tumors and Brain Tumors. Clin Cancer Res 2019;25:2560-74. [Crossref] [PubMed]
- Wang J, Chong KK, Nakamura Y, et al. B7-H3 associated with tumor progression and epigenetic regulatory activity in cutaneous melanoma. J Invest Dermatol 2013;133:2050-8. [Crossref] [PubMed]
- Yamato I, Sho M, Nomi T, et al. Clinical importance of B7-H3 expression in human pancreatic cancer. Br J Cancer 2009;101:1709-16. [Crossref] [PubMed]
- Tang J, Jiang W, Liu D, et al. The comprehensive molecular landscape of the immunologic co-stimulator B7 and TNFR ligand receptor families in colorectal cancer: immunotherapeutic implications with microsatellite instability. Oncoimmunology 2018;7:e1488566. [Crossref] [PubMed]
- Vogl UM, Öhler L, Rasic M, et al. Evaluation of Prognostic Immune Signatures in Patients with Breast, Colorectal and Pancreatic Cancer Receiving Chemotherapy. Anticancer Res 2017;37:1947-55. [Crossref] [PubMed]
- Bin Z, Guangbo Z, Yan G, et al. Overexpression of B7-H3 in CD133+ colorectal cancer cells is associated with cancer progression and survival in human patients. J Surg Res 2014;188:396-403. [Crossref] [PubMed]
- Hu X, Xu M, Hu Y, et al. B7-H3, Negatively Regulated by miR-128, Promotes Colorectal Cancer Cell Proliferation and Migration. Cell Biochem Biophys 2021;79:397-405. [Crossref] [PubMed]
- Wu J, Wang F, Liu X, et al. Correlation of IDH1 and B7H3 expression with prognosis of CRC patients. Eur J Surg Oncol 2018;44:1254-60. [Crossref] [PubMed]
- Mao Y, Chen L, Wang F, et al. Cancer cell-expressed B7-H3 regulates the differentiation of tumor-associated macrophages in human colorectal carcinoma. Oncol Lett 2017;14:6177-83. [Crossref] [PubMed]
- Fan H, Zhu JH, Yao XQ. Prognostic significance of B7-H3 expression in patients with colorectal cancer: A meta-analysis. Pak J Med Sci 2016;32:1568-73. [Crossref] [PubMed]
- Zhang W, Acuna-Villaorduna A, Kuan K, et al. B7-H3 and PD-L1 Expression Are Prognostic Biomarkers in a Multi-racial Cohort of Patients with Colorectal Cancer. Clin Colorectal Cancer 2021;20:161-9. [Crossref] [PubMed]
- Zhang T, Jin Y, Jiang X, et al. Clinical and Prognostic Relevance of B7-H3 and Indicators of Glucose Metabolism in Colorectal Cancer. Front Oncol 2020;10:546110. [Crossref] [PubMed]
- Ingebrigtsen VA, Boye K, Tekle C, et al. B7-H3 expression in colorectal cancer: nuclear localization strongly predicts poor outcome in colon cancer. Int J Cancer 2012;131:2528-36. [Crossref] [PubMed]
- Tang YC, Zhang Y, Zhou J, et al. Ginsenoside Rg3 targets cancer stem cells and tumor angiogenesis to inhibit colorectal cancer progression in vivo. Int J Oncol 2018;52:127-38. [PubMed]
- Zhang C, Chen Y, Li F, et al. B7-H3 is spliced by SRSF3 in colorectal cancer. Cancer Immunol Immunother 2021;70:311-21. [Crossref] [PubMed]
- Lu Z, Zhao ZX, Cheng P, et al. B7-H3 immune checkpoint expression is a poor prognostic factor in colorectal carcinoma. Mod Pathol 2020;33:2330-40. [Crossref] [PubMed]
- Ingebrigtsen VA, Boye K, Nesland JM, et al. B7-H3 expression in colorectal cancer: associations with clinicopathological parameters and patient outcome. BMC Cancer 2014;14:602. [Crossref] [PubMed]
- Mao Y, Sun J, Wang WP, et al. Clinical significance of costimulatory molecule B7-H3 expression on CD3(+) T cells in colorectal carcinoma. Chin Med J (Engl) 2013;126:3035-8. [PubMed]
- Liu X, Wang F, Wu J, et al. Expression of CYP1B1 and B7-H3 significantly correlates with poor prognosis in colorectal cancer patients. Int J Clin Exp Pathol 2018;11:2654-64. [PubMed]
- McShane LM, Altman DG, Sauerbrei W, et al. REporting recommendations for tumour MARKer prognostic studies (REMARK). Br J Cancer 2005;93:387-91. [Crossref] [PubMed]
- Liao Q, Zeng Z, Guo X, et al. LPLUNC1 suppresses IL-6-induced nasopharyngeal carcinoma cell proliferation via inhibiting the Stat3 activation. Oncogene 2014;33:2098-109. [Crossref] [PubMed]
- Gauthier J, Wu QV, Gooley TA. Cubic splines to model relationships between continuous variables and outcomes: a guide for clinicians. Bone Marrow Transplant 2020;55:675-80. [Crossref] [PubMed]
- Durrleman S, Simon R. Flexible regression models with cubic splines. Stat Med 1989;8:551-61. [Crossref] [PubMed]
- Schoenfeld D. Partial residuals for the proportional hazards regression model. Biometrika 1982;69:239-41. [Crossref]
- Rich JT, Neely JG, Paniello RC, et al. A practical guide to understanding Kaplan-Meier curves. Otolaryngol Head Neck Surg 2010;143:331-6. [Crossref] [PubMed]
- Li T, Fu J, Zeng Z, et al. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res 2020;48:W509-14. [Crossref] [PubMed]
- Li T, Fan J, Wang B, et al. TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res 2017;77:e108-10. [Crossref] [PubMed]
- Qin S, Xu L, Yi M, et al. Novel immune checkpoint targets: moving beyond PD-1 and CTLA-4. Mol Cancer 2019;18:155. [Crossref] [PubMed]
- Hajian-Tilaki K. Receiver Operating Characteristic (ROC) Curve Analysis for Medical Diagnostic Test Evaluation. Caspian J Intern Med 2013;4:627-35. [PubMed]
- Hofmeyer KA, Ray A, Zang X. The contrasting role of B7-H3. Proc Natl Acad Sci U S A 2008;105:10277-8. [Crossref] [PubMed]
- Zhang G, Wang J, Kelly J, et al. B7-H3 augments the inflammatory response and is associated with human sepsis. J Immunol 2010;185:3677-84. [Crossref] [PubMed]
- Chen X, Quinn EM, Ni H, et al. B7-H3 participates in the development of experimental pneumococcal meningitis by augmentation of the inflammatory response via a TLR2-dependent mechanism. J Immunol 2012;189:347-55. [Crossref] [PubMed]
- Kontos F, Michelakos T, Kurokawa T, et al. B7-H3: An Attractive Target for Antibody-based Immunotherapy. Clin Cancer Res 2021;27:1227-35. [Crossref] [PubMed]
- Wang L, Zhao Y, Xu M, et al. Serum miR-1301-3p, miR-335-5p, miR-28-5p, and their target B7-H3 may serve as novel biomarkers for colorectal cancer. J BUON 2019;24:1120-7. [PubMed]
- Yasui K, Kondou R, Iizuka A, et al. Effect of preoperative chemoradiotherapy on the immunological status of rectal cancer patients. J Radiat Res 2020;61:766-75. [Crossref] [PubMed]
- Wang J, Chen X, Xie C, et al. MicroRNA miR-29a Inhibits Colon Cancer Progression by Downregulating B7-H3 Expression: Potential Molecular Targets for Colon Cancer Therapy. Mol Biotechnol 2021;63:849-61. [Crossref] [PubMed]
- Shi T, Ma Y, Cao L, et al. B7-H3 promotes aerobic glycolysis and chemoresistance in colorectal cancer cells by regulating HK2. Cell Death Dis 2019;10:308. [Crossref] [PubMed]
- Picarda E, Ohaegbulam KC, Zang X. Molecular Pathways: Targeting B7-H3 (CD276) for Human Cancer Immunotherapy. Clin Cancer Res 2016;22:3425-31. [Crossref] [PubMed]
- Wang Y, Zhang B, Huang H, et al. An improved method to build lung cancer PDX models by surgical resection samples and its association with B7-H3 expression. Transl Cancer Res 2019;8:2848-57. [Crossref]
- Yin Y, Yao S, Hu Y, et al. The Immune-microenvironment Confers Chemoresistance of Colorectal Cancer through Macrophage-Derived IL6. Clin Cancer Res 2017;23:7375-87. [Crossref] [PubMed]
- Ye L, Zhang T, Kang Z, et al. Tumor-Infiltrating Immune Cells Act as a Marker for Prognosis in Colorectal Cancer. Front Immunol 2019;10:2368. [Crossref] [PubMed]
- Wang R, Ma Y, Zhan S, et al. B7-H3 promotes colorectal cancer angiogenesis through activating the NF-kappaB pathway to induce VEGFA expression. Cell Death Dis 2020;11:55. [Crossref] [PubMed]
- Sun J, Chen LJ, Zhang GB, et al. Clinical significance and regulation of the costimulatory molecule B7-H3 in human colorectal carcinoma. Cancer Immunol Immunother 2010;59:1163-71. [Crossref] [PubMed]
- Meng F, Yang M, Chen Y, et al. miR-34a induces immunosuppression in colorectal carcinoma through modulating a SIRT1/NF-kappaB/B7-H3/TNF-alpha axis. Cancer Immunol Immunother 2021;70:2247-59. [Crossref] [PubMed]
- Li Y, Huang J, Foley NM, et al. B7H3 ameliorates LPS-induced acute lung injury via attenuation of neutrophil migration and infiltration. Sci Rep 2016;6:31284. [Crossref] [PubMed]
- Castriconi R, Dondero A, Augugliaro R, et al. Identification of 4Ig-B7-H3 as a neuroblastoma-associated molecule that exerts a protective role from an NK cell-mediated lysis. Proc Natl Acad Sci U S A 2004;101:12640-5. [Crossref] [PubMed]
- Yi KH, Chen L. Fine tuning the immune response through B7-H3 and B7-H4. Immunol Rev 2009;229:145-51. [Crossref] [PubMed]
- Patel SA, Minn AJ. Combination Cancer Therapy with Immune Checkpoint Blockade: Mechanisms and Strategies. Immunity 2018;48:417-33. [Crossref] [PubMed]
- Qi Y, Liu B, Sun Q, et al. Immune Checkpoint Targeted Therapy in Glioma: Status and Hopes. Front Immunol 2020;11:578877. [Crossref] [PubMed]
- Ye Q, Liu J, Xie K. B7 Family Proteins in Cancer Progression: Immunological and Non-Immunological Functions. J Cancer Treatment Diagn 2019;3:1-6. [Crossref]
- Zhou WT, Jin WL. B7-H3/CD276: An Emerging Cancer Immunotherapy. Front Immunol 2021;12:701006. [Crossref] [PubMed]
- McCartney DL, Walker RM, Morris SW, et al. Identification of polymorphic and off-target probe binding sites on the Illumina Infinium MethylationEPIC BeadChip. Genom Data 2016;9:22-4. [Crossref] [PubMed]
- Zhou Q, Costinean S, Croce CM, et al. MicroRNA 29 targets nuclear factor-kappaB-repressing factor and Claudin 1 to increase intestinal permeability. Gastroenterology 2015;148:158-69.e8. [Crossref] [PubMed]
(English Language Editor: J. Jones)