B7-H6 is aberrantly expressed in tumors
To explore the difference of B7-H6 expression between tumors and normal tissues, we compared the messenger RNA (mRNA) expression levels of B7-H6 between tumors and normal tissues of various cancer types from the TCGA database and GTEx. The results showed that compared to correspondingly normal tissue, B7-H6 expression was significantly higher in 9 types of tumor tissues (CHOL, COAD, ESCA, KICH, KIRC, KIRP, PAAD, READ and STAD) while significantly lower in 15 types of tumor tissues (ACC, BLCA, CESC, GBM, LAML, LIHC, LUAD, LUSC, OV, PRAD, SKCM, TGCT, THCA, UCEC and UCS) (Fig. 1a).
B7-H6 expression in ccRCC tumor samples (n = 539) and normal samples (n = 100) was also compared and the result demonstrated that B7-H6 expression was markedly increased within ccRCC cancer samples compared to normal kidney samples (Fig. 1b). Furthermore, the expression of B7-H6 in 72 paired ccRCC samples and adjacent normal samples was also analyzed, we found that the B7-H6 expression in ccRCC tissues was significantly upregulated compared with the expression in the adjacent normal tissues (Fig. 1c).
Association Between B7-h6 Expression And Clinicopathological Characteristics In Ccrcc Patients
The B7-H6 expression data and clinicopathological information from 539 patients with ccRCC were obtained from the TCGA database KIRC cohort for analysis. The ccRCC patients in TCGA KIRC cohort were divided into B7-H6 high-expression and B7-H6 low-expression groups with the median value of B7-H6 expression as a cut-off value. The relationship between the clinical parameters and expression level of B7-H6 was shown in Fig. 2. The results showed that the mRNA expression level of B7-H6 gradually decreased with the increase of tumor stage, histologic grade and pathologic stage in ccRCC (Fig. 2a, 2c and 2d). The expression of B7-H6 in patients with distant metastasis was significantly lower than patients without distant metastasis. Still, the expression level was no significant difference between ccRCC patients with or without lymphatic metastasis (Fig. 2c and 2b). The study also found that the expression of B7-H6 was significantly higher in the survive group than in the deceased group. We then analyzed the correlations between clinical parameters and B7-H6 expression and found that B7-H6 expression was significantly correlated with T (tumor) stage (p < 0.001), M (distant metastasis) stage (p < 0.001), Pathologic stage (p < 0.001), gender (p = 0.025) and as well as age (p = 0.043). (Table 1).
The up-regulated B7-H6 expression, younger age (< 60) and negative distant metastasis are independent prognostic factors of favorable prognosis.
High B7-h6 Expression Is Associated With Favorable Survival In Ccrcc
Overall survival (OS), disease specific survival (DSS), relapse free survival (RFS) and progression free survival (PFS) are four common prognostic monitoring indexes used to basically summarize the prognosis and survival of cancer patients. To explore the potential prognostic significance of B7-H6 expression in ccRCC, different datasets were used to analyze the patient survival between low and high B7-H6 expression level groups in ccRCC patients. Clinical data for ccRCC were downloaded from TCGA database and analyzed by R software. The result showed that high B7-H6 expression was significantly positively associated with the OS, DSS and progress free interval of patients with ccRCC (Fig. 3a, 3b and 3c). Using the Kaplan-Meier plotter database, we found that the high expression level of B7-H6 was significantly correlated with the favorable prognostic survival rate (HR = 0.35, p = 3.3e-10, Fig. 3d). We also generated survival curves for overall survival and disease-free survival based on GEPIA using the log-rank test and the Mantel-Cox test. High B7-H6 expression was significantly associated with favorable OS and DFS (Fig. 3e, 3f). These results indicated that low expression of B7-H6 may be one of the causes of poor prognosis, and this gene could be a prognostic biomarker for ccRCC patients.
To explore whether low expression of B7-H6 is an independent risk factor for OS in ccRCC patients, univariate and multivariate Cox analyses were performed. The univariate Cox analysis indicated that low B7-H6 expression was associated with a worse overall survival time. Other clinical parameters, such as Age, Histologic grade (G3&G4 vs. G1&G2), Pathologic stage (Stage III & Stage IV vs. Stage I & Stage II), T stage (T3 & T4 vs. T1 & T2), N stage (N1 vs. N0) and M stage (M1 vs. M0) also correlated with the overall survival time (Table 2). For verifying the prognostic value of B7-H6 in ccRCC patients, we performed the multivariate analysis. The results showed that only B7-H6 expression level and M stage (M1 vs. M0) were independently associated with the overall survival time, which means that B7-H6 expression level in evaluating patients' clinical prognosis is superior to T stage, N stage and Pathologic stage. These findings indicated that B7-H6 expression was an independent protective factor for OS in ccRCC patients (HR = 0.477, p = 0.002, Table 2).
Correlation Analysis Between B7-h6 Expression And Tumor Infiltrating Immune Cells
Tumor infiltrating immune cells affect the antitumor efficacy and survival of patients in various cancers [21, 22]. ccRCC is one type of immunogenic cancers with abundant infiltrating immune cells [1, 23]. Therefore, we analyzed the correlation of B7-H6 expression with tumor infiltrating immune cells and tumor purity with TIMER database. The results displayed that the B7-H6 expression level had obviously positive correlation with tumor purity(r = 0.184, p = 6.72e-05), infiltrating levels of B cells (r = 0.187, p = 5.44e-05), CD4+ T cells (r = 0.143, p = 2.05e-03), macrophages (r = 0.281, p = 1.38e-09), neutrophils (r = 0.224, p = 1.24e-06), and dendritic cells (r = 0.201, p = 1.49e-05) in ccRCC, but no significant association with CD8+ T cells (Fig. 4a).
To further determine the effect of B7-H6 expression on tumor immune cell infiltration, we performed ssGSEA to assess the association between the immune cell infiltration pattern and B7-H6 level based on transcriptome profiling data for 539 KIRC cohort in TCGA database. Spearman correlation analyses revealed that high B7-H6 expression was mainly associated with low infiltration of some immune cell types (Fig. 4b), especially Treg cells and cytotoxic cells. We also observed positive correlations of B7-H6 expression level with infiltration of several types of immune cells, such as Eosinophils, central Memory T (Tcm) cells, and T helper 17 (Th17) cells (Fig. 4b).
Furthermore, the infiltrating proportions and differences of immune cells between the B7-H6 high and B7-H6 low groups were estimated via the CIBERSORT tool. The results displayed that Eosinophils (p < 0.001), central Memory T (Tcm) cells (p < 0.001), T helper 17 (Th17) cells (p < 0.001), neutrophils (p < 0.001) and Th cells (p < 0.001) shared a higher proportion in B7-H6 high expression group compared with low expression group (Fig. 4c). While the proportion of Treg cells, cytotoxic cells, CD56bright NK cells, Th1 and Th2 cells were significantly lower (p < 0.001) (Fig. 4c). We further investigated if the expression of B7-H6 was associated with immune biomarker sets of different immune cells in ccRCC. As expected, a significant correlation was obtained between the expression of B7-H6 and most of the immune markers of immune cells in ccRCC after tumor purity modulation was performed (Table 3). Results showed that B7-H6 expression was strongly linked to STAT5B (Treg cell), STAT6 (Th2 cell), and STAT3 (Th17 cell) (Fig. 5n and Fig. 5o). B7-H6 was also positively linked to STAT1, TNF-α (Th1 cell), IRF5, INOS (M1 Macrophage), CD163, MS4A4A (M2 Macrophage), most markers of Neutrophils (CD66b, CD15, CD11b), some markers of dendritic cell (HLA-DPB1, HLA-DQB1, HLA-DRA, HLA-DPA1, BDCA1, BDCA4, CD11c and CD144) (Table 3). Moreover, B7-H6 levels showed a negative association with markers of natural killer cell (KIR2DL4, XCL1and CD7) and B cell (CD19 and CD79a) (Fig. 5). Since B7-H6 expression was significantly negatively correlated to infiltration of NK cells, we then evaluated the correlation between B7-H6 and its receptor NKp30 (NCR3) of NK cells. The analyses revealed that B7-H6 expression was inversely correlated with NKp30 expression (r = -0.130, p = 0.002), indicating a potential NK inhibitory role for B7-H6 in ccRCC (Fig. 4). In addition, B7-H6 expression in ccRCC showed a positive correlation with TIM-3 in T cell exhaustion but a negative association with PD-1, LAG3 and GZMB (Table 3 and Fig. 5).
Overall, the results demonstrated that B7-H6 had a complicated relationship with immune cell infiltration in ccRCC.