Comprehensive Analysis of the Expression and Prognosis Value of Chromobox Family Members in Clear Cell Renal Cell Carcinoma

Abstract

Clear cell renal cell carcinoma (ccRCC) accounts for 80% of all renal cancers and has a poor prognosis. Chromobox (CBX) family protein expression has been reported in a variety of human malignancies, but the roles of CBXs in ccRCC remain unclear. In this study, by using ONCOMINE, UALCAN, GEPIA, Kaplan-Meier Plotter, cBioPortal, and TIMER, we found the transcriptional levels of CBX3 and CBX4 in ccRCC tissues were significantly higher than those in normal kidney tissues, whereas the transcriptional levels of CBX1, CBX5, CBX6, and CBX7 were significantly reduced in ccRCC tissues. The promoters of CBX2, CBX3, CBX4, CBX5, CBX6, CBX7, and CBX8 were hypermethylated, whereas the CBX1 promoter was hypomethylated in ccRCC. The expression of CBX1, CBX3, CBX4, CBX5, CBX6, and CBX7 was significantly associated with clinicopathological parameters in ccRCC patients. ccRCC patients with high expression levels of CBX3, CBX4, and CBX8 and low expression levels of CBX1, CBX5, CBX6, and CBX7 showed a strong association with poor overall survival. Genetic alterations in CBXs were correlated with poor overall survival and disease-free survival in patients with ccRCC. Moreover, we found significant associations between the expression of CBXs and infiltration of immune cells (B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells). Our results provide novel insights into the development of CBX-based biomarkers and therapeutic targets for ccRCC.

Introduction

Renal cell carcinoma (RCC) is the most prevalent urological cancer with a steady increase in mortality and morbidity in recent years (1). Clear cell renal cell carcinoma (ccRCC), the most prevalent subtype of RCC, accounts for more than 75% of the RCC cases (2). Despite recent advances in the detection and treatment of RCC, the clinical outcomes of patients with metastatic RCC remain poor. The 5-year survival rate is less than 10% (3). In the last few decades, several studies have been carried out on the pathogenesis of ccRCC, but the detailed molecular events of ccRCC occurrence and development are not yet fully understood. Therefore, it is necessary to identify more effective therapeutic targets and prognostic markers for ccRCC.

Chromobox (CBX) family members are typical components of the polycomb repressive complex (PRC) that epigenetically modifies chromatin and participates in multiple physiological processes, such as differentiation, senescence, and DNA repair (4). CBX family proteins can be divided into two types: heterochromatin protein 1 (HP1) group (containing CBX1, CBX3, and CBX5), with an N-terminal chromodomain, and polycomb (Pc) group (containing CBX2, CBX4, CBX6, CBX7, and CBX8), with a conserved N-terminal chromodomain and a C-terminal Pc repressor box (5). Evidence has demonstrated that CBXs play complex roles in human tumors. CBXs can exert both anti-tumor and pro-tumor activities depending on the tumor type and cellular context. For example, studies have shown that CBX4 can promote angiogenesis of liver tumors through HIF-α sumoylation, which is a negative factor for the prognosis of patients with hepatocellular carcinoma (6). CBX4 also inhibits the expression of Runx2 by recruiting HDAC3, thereby inhibiting tumor metastasis in colorectal cancer (7). CBX7 shows inhibitory activity in thyroid, colon, and lung cancer (8–10), while playing an oncogenic role in gastric cancer (11). In ccRCC, CBX4 transcriptionally suppresses KLF6 by interacting with HDAC1 to exert oncogenic activities (12). However, the roles of other members of the CBX family in the development and progression of ccRCC remain unknown.

In this study, we aimed to analyze the expression and clinical relevance of CBX family members in ccRCC based on online public databases. We analyzed expression patterns, clinicopathological parameters, prognostic values, genetic alterations, and immune cell infiltration. Our findings indicate that CBXs may have complex and distinct functions in ccRCC and serve as potential therapeutic targets for the clinical intervention of ccRCC in the future.

Materials and Methods

UALCAN

UALCAN (http://ualcan.path.uab.edu) is a comprehensive and interactive website that provides cancer data based on TCGA database (13). In this study, data on the mRNA expression and promoter methylation levels as well as the association of the transcriptional expression of CBXs with clinicopathologic parameters in ccRCC patients were obtained using UALCAN. The Student’s t-test was used to compare differences, and statistical significance was set at P < 0.05.

ONCOMINE

ONCOMINE (www.oncomine.org) is an online cancer database used for powerful genome-wide expression analysis (14). Data on significant changes in the transcription levels of CBXs between ccRCC and normal kidney tissues were obtained from the ONCOMINE database. The cutoffs of the P value and fold change were 0.01 and 1.5, respectively.

Human Protein Atlas

The Human Protein Atlas (https://www.proteinatlas.org) is a website that provides the distribution and expression of protein in normal human tissues, tumor tissues, cell lines, and blood cells by immunohistochemistry (15). In this study, immunohistochemistry images were used to directly compare the protein expression of different CBX family members in human normal and ccRCC tissues.

GEPIA

GEPIA (http://gepia.cancer-pku.cn/index.html) is a newly developed website for interactive gene expression profiling analysis (16). In this study, disease-free survival (DFS) analysis of CBXs in ccRCC patients was performed using a Kaplan-Meier curve, with hazard ratios (HRs) and 95% confidence intervals (CIs). In addition, the relative expression levels of CBXs in ccRCC tissues were analyzed using GEPIA.

Kaplan-Meier Plotter

The prognostic value of CBXs in patients with ccRCC was analyzed using the Kaplan-Meier plotter (http://kmplot.com/analysis/) (17). Based on CBX expression levels, ccRCC patients were divided into two groups and the data were assessed using a Kaplan-Meier survival plot. Data were considered statistically significant at P < 0.05.

cBioPortal

cBioPortal (http://www.cbioportal.org/) is a comprehensive web resource that can be used to explore multidimensional cancer genomics data (18). We obtained the genetic alterations of CBXs in ccRCC from cBioPortal (TCGA, PanCancer Atlas). mRNA expression z-scores (RNA Seq V2 RSEM) were obtained using a z-score threshold of ±1.8. Moreover, genetic alterations in CBXs and their relationships with overall survival (OS) and the disease-free survival (DFS) in patients with ccRCC were analyzed using Kaplan-Meier curves. Statistical significance was set at P < 0.05.

TIMER

Timer (https://cistrome.shinyapps.io/timer/) is a website for the comprehensive analysis of tumor-infiltrating immune cells and their impact on clinical outcomes (19). In this study, the “Gene module” was used to explore the correlation between CBX expression and the abundance of immune infiltrates, and the “Survival module” was used to analyze the association between the infiltration of immune cells and CBXs with the clinical outcome of ccRCC patients.

GSEA Analysis

Gene Set Enrichment Analysis (GSEA) is a computational method that determines whether a set of priori defined genes shows statistically significant, concordant differences between two biological states. GSEA was obtained from the Broad Institute (http://software.broadinstitute.org/gsea/index.jsp) (20). The cut-off values were predefined as P < 0.05 and FDR < 0.25.

Cell Culture

The ccRCC cell lines (786-O and 769-P) were purchased from Chinese Academy of Sciences Cell Bank (Shanghai, China) and HK2 cell was provided by NHC Key Laboratory of Cancer Proteomics, Xiangya Hospital, Central South University (Changsha, China). 786-O and 769-P cells were cultured in RPMI-1640 medium (Gibco, USA), and HK2 cell was cultured in Dulbecco’s modified Eagle’s medium (DMEM, Gibco, USA), both were supplemented with 10% fetal bovine serum (Gibco, USA) and 1% penicillin-streptomycin (Gibco, USA) in a humidified incubator at 37°C containing 5% CO₂.

RNA Extraction and qRT-PCR

Total RNA was extracted from cell using Trizol reagent (Takara, Dalian, China) and the RNA quality and concentration were assessed using a NanoDrop™ 1000 Spectrophotometer (Thermo Fisher Scientific, USA). The OD260/OD280 ratios were 1.8-2.0, and the OD260/OD230 ratios were 2.0-2.2. Then the complementary DNA was synthesized by reverse transcription reaction using the TaqMan MicroRNA Reverse Transcription Kit (Applied Biosystems, USA). The qRT-PCR was performed with miScript SYBR Green PCR kit (Qiagen) using an ABI7500 instrument (Applied Biosystems). The PCR conditions were 95°C for 30 sec, 45 cycles at 94°C for 15 sec and one step at 60°C for 31 sec, and 55°C for 30 sec. ACTB was used as a normalizer for RNA quantification. The data were analyzed using the 2^–ΔΔCt method.

Western Blotting

Cells were lysed by lysis buffer (150 mM NaCl, 50 mM Tris-HCl, pH7.5) supplemented with inhibitors (1 mM phenylmethylsulphonyl fluoride, 1 mg/ml aprotinin, 1 mg/ml leupeptin, 1 mg/ml pepstatin, 1 mM Na₃VO₄, 1 mM NaF). Then cell lysates were centrifuged at 4°C, 14000 g for 20 minutes to remove cell debris. Bicinchonininc acid (BCA) assay kit (BCA Assay Reagent, Pierce Chemical, America) was used to determine the protein content. For western blotting, protein samples were separated by SDS-PAGE and transferred to the polyvinylidene fluoride (PVDF) membrane. Using primary antibodies incubated the target proteins overnight at 4°C, then the membranes were incubated with corresponding secondary antibodies, including anti-rabbit and anti-mouse horseradish peroxidase (HRP)-linked IgG (Proteintech,1:5000) for 2 hours. The working dilution of primary antibodies were: CBX3: 1:1000 (Immunoway, YM3741); CBX6: 1:1000 (Immunoway, YT0699); CBX7: 1:1000 (Sangon Biotech, D262675); α-Tubulin: 1:5000 (Proteintech, 66031-1-Ig).

Immunohistochemistry

Tissue microarrays containing 9 human ccRCC tissues and paired adjacent normal kidney tissues were purchased from Avilabio (DC-Kid11008). The experimental operation was conducted according to the protocol of the immunohistochemical staining kit (ZSGB-BIO, PV-9000). The working dilution of primary antibodies were: CBX3: 1:200 (Immunoway, YM3741); CBX6: 1:200 (Immunoway, YT0699); CBX7: 1:50 (Sangon Biotech, D262675). The research was approved by the local ethics committee. To quantify the level of protein expression, five random fields in each tissue sample were photographed with a microscope and analyzed by two experienced pathologists. The expression levels were scored as proportion of immunopositive staining area (0 for ≤5% positive cells, 1 for 6–25%, 2 for 26–50%, 3 for 51–75%, and 4 for ≥76%) multiplied by intensity of staining (0 for negative staining, 1 for weak staining, 2 for moderate staining, and 3 for strong staining).

Statistical Analysis

Data were presented as the means ± standard error of the mean (SEM). Statistical analyses were done using SPSS. Student’s t-test were used for comparison between groups. P < 0.05 were considered statistically significant.

Results

Expression Levels and Promoter Methylation Status of CBXs in Patients With ccRCC

We firstly analyzed the mRNA expression of different CBX family members in ccRCC patients using the UALCAN and ONCOMINE databases. As shown in Figure 1, the transcriptional levels of CBX3 and CBX4 in ccRCC tissues were significantly higher than those in normal kidney tissues, whereas the transcriptional levels of CBX1, CBX5, CBX6, and CBX7 were significantly reduced in ccRCC tissues. CBX2 and CBX8 were expressed at similar levels in ccRCC and normal kidney tissues. Next, the transcription levels of the CBX family members were measured using the ONCOMINE database. The expression levels of CBX1, CBX3, CBX4, CBX5, and CBX7 also showed significant changes in ccRCC (Table 1). And the mRNA expressions of CBXs in normal kidney cell HK2 and ccRCC cells 786-O and 769-P were analyzed using qRT-PCR. The transcriptional levels of CBX3 and CBX4 in ccRCC cells were significantly higher than those in normal kidney cell, whereas the transcriptional levels of CBX1, CBX5, CBX6, and CBX7 were significantly reduced in ccRCC cells (Figure 2), which was consistent with the database. The abundance of CBX2 and CBX8 in both normal kidney cell and ccRCC cells were very low, and the qRT-PCR cycle number were more than 45.

Figure 1

Figure 2

We also explored the protein expression of CBXs in ccRCC. Firstly, we examined the protein expression levels of CBX3, CBX5 and CBX6 in HK2, 786-O and 769-P cell lines. The results showed that CBX3 protein was much higher in ccRCC cell lines than that in HK2 cells, while CBX6 and CBX7 were significantly reduced in ccRCC cells (Figure 3). And Human Protein Atlas data showed that low protein expressions of CBX3 and CBX4 were found in normal kidney tissues, while medium and high protein expressions of them were observed in ccRCC tissues (Figures 4C, D). However, higher protein expressions of CBX1, CBX5, CBX6 and CBX7 were observed in normal kidney tissues than in ccRCC tissues (Figures 4A, E–G). CBX2 protein was not expressed both in normal kidney tissues and ccRCC tissues (Figure 4B). CBX8 protein was expressed at similar level in ccRCC and normal kidney tissues (Figure 4H). Our own immunohistochemical results suggestd that CBX3 was overexpressed in ccRCC tissues, while expression levels of CBX6 and CBX7 were reduced in ccRCC tissues (Figure 5), which were consistent with the database. We also compared the relative expression levels of CBXs in ccRCC tissues and found that the relative expression of CBX3 was the highest among all CBXs (Figure 6). Moreover, given the important role of promoter methylation in the process of renal tumors, we used UALCAN to investigate the promoter methylation levels of CBXs in ccRCC and found that the promoters of CBX2, CBX3, CBX4, CBX5, CBX6, CBX7, and CBX8 were hypermethylated in ccRCC tissues, while the CBX1 promoter was hypomethylated in ccRCC tissues, compared to those in the normal kidney tissues (Figure 7).

Figure 3

Figure 4

Figure 5

Figure 6

Figure 7

The Relationship of mRNA Levels of CBXs With Clinicopathological Parameters in Patients With ccRCC

After mRNA expression and promoter methylation were analyzed, we next explored the relationship between mRNA levels of CBXs and clinicopathological parameters in ccRCC patients using UALCAN; this included the patients’ cancer stages and tumor grades. The mRNA expression of CBXs was significantly correlated with cancer stage and tumor grade. As the tumor progressed, patients with more advanced tumor grades tended to express higher mRNA expression of CBX3 and CBX4 but lower mRNA expression of CBX1, CBX5, CBX6, and CBX7 (Figure 8). In addition, as shown in Figure 9, low expression of CBX1, CBX5, CBX6, and CBX7 and high expression of CBX3 and CBX4 were significantly correlated with advanced cancer stages. CBX2 and CBX8 expression did not differ significantly with cancer stage or tumor grade. In conclusion, mRNA expression of CBX1, CBX3, CBX4, CBX5, CBX6, and CBX7 was significantly associated with clinicopathological parameters in ccRCC patients.

Figure 8

Figure 9

The Prognostic Value of CBXs in Patients With ccRCC

Furthermore, we analyzed the prognostic value of CBXs in patients with ccRCC using the Kaplan-Meier plotter and GEPIA database. ccRCC patients with high expression levels of CBX3, CBX4, and CBX8, and low expression levels of CBX1, CBX5, CBX6, and CBX7 showed a strong association with poor OS (Figure 10). The value of differentially expressed CBXs in DFS of ccRCC patients was also evaluated. The results showed that low expression levels of CBX1 and CBX7 were significantly correlated with poor DFS (Figure 11). These results indicated that the expression levels of CBX1, CBX3, CBX4, CBX5, CBX6, CBX7, and CBX8 were significantly associated with ccRCC patient prognosis, and they might be exploited as useful biomarkers for predicting the survival of patients with ccRCC.

Figure 10

Figure 11

Genetic Alteration of CBXs and Their Relationship With OS and DFS in Patients With ccRCC

Next, we applied the cBioPortal online tool (TCGA, PanCancer Atlas) to investigate genetic alterations in CBXs and their relationship with OS and DFS in ccRCC. In the 510 sequenced ccRCC patients, genetic alterations were found in 192 patients, and the alteration rate was 38%. The OncoPrints included missense mutation, truncating mutation, fusion, deletion, amplification, as well as high and low levels of mRNA. The rates of genetic alterations of CBX1, CBX2, CBX3, CBX4, CBX5, CBX6, CBX7, and CBX8 were 6%, 5%, 11%, 8%, 9%, 7%, 9%, and 4%, respectively (Figures 12A, B). Moreover, Kaplan-Meier plots and log-rank tests indicated that genetic alterations in CBXs were correlated with poor OS and DFS in patients with ccRCC (Figures 12C, D)

Figure 12

Immune Cell Infiltration of CBXs in Patients With ccRCC

Inflammatory responses and immune cell infiltration can affect the clinical outcomes of cancer patients. Therefore, we comprehensively analyzed the association between CBX expression levels and immune cell infiltration using the TIMER database. CBX1, CBX2, CBX4, CBX5, CBX6, and CBX7 expression was positively correlated with the infiltration of B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells (Figures 13A, B, D–G). There was a positive correlation between CBX3 expression and the infiltration of B cells, CD8+ T cells, macrophages, neutrophils, and dendritic cells (Figure 13C). Similarly, the expression of CBX8 was positively associated with the infiltration of CD4+ T cells and negatively associated with the infiltration of CD8+ T cells (Figure 13H). Moreover, we used the Cox proportional hazard model and found that CD8+ T cells, CBX4, CBX7, and CBX8 were significantly correlated with the clinical outcome of ccRCC patients after adjustment for confounding factors (Table 2).

Figure 13

GSEA Analysis of the Gene Expression Files in CBXs mRNA High Subgroup Versus Low Subgroup

We also used GSEA analysis to determine the gene expression files in CBX family mRNA high versus low. GSEA results showed that several gene sets that were significantly enriched in CBXs mRNA high subgroup compared to low subgroup, including Wnt signaling pathway, oxidative phosphorylation, tight junction, endocytosis, spliceosome. While proximal tubule bicarbonate reclamation, cytosolic DNA-sensing pathway and pathways in cancer were significantly enriched in CBXs mRNA low subgroup compared to high subgroup (Supplementary Figure 1).

Discussion

Epigenetic modification has been proven to be related to multiple physiological processes, including differentiation, senescence, and DNA repair (27–30). As important components of epigenetic regulation complexes, CBX family proteins can mediate the recruitment of PRC1 to chromatin and participate in the occurrence and progression of various tumors (31, 32). CBX family proteins function as both oncogenes and tumor suppressors, depending on the tumor type and cellular context. However, to date, the prognostic value and biological function of CBX family members in the progression of ccRCC remain unknown. In this study, we found that six out of eight CBX family members were differentially expressed in ccRCC, and seven CBX family members were significantly associated with the survival of ccRCC patients. In addition, we showed that the expression of six CBX family members was closely related to cancer stage and tumor grade in ccRCC. The expression of CBX family members was associated with the infiltration of six immune cells and influenced the outcome of ccRCC patients. Moreover, we found that promoter methylation status and genetic alteration might cause dysregulation of CBX family members in ccRCC. For most members of the CBX family, our study was the first to systematically analyze the their expression and prognostic significance in ccRCC.

In previous studies, overexpression of CBX1 has been found in many cancers such as hepatocellular carcinoma (33, 34), gastric cancer (35, 36), colorectal cancer (37), and pituitary cancer (38). In HCC, CBX1 functions as an oncogene by interacting with the transcription factor HMGA2 to activate the Wnt/β-catenin signaling pathway (33). High expression of CBX1 is significantly associated with a large tumor size, tumor vascular invasion, poor tumor differentiation, and prognosis in HCC patients (33, 34). High CBX1 expression was found to be significantly associated with poor prognosis in gastric cancer patients receiving adjuvant 5-fluorouracil-based chemotherapy (35). However, little is known about the role of CBX1 in ccRCC. In this study, we found that the mRNA expression of CBX1 was downregulated in ccRCC tissues compared to that in normal tissues. Low CBX1 levels are associated with advanced cancer stage and tumor grade. Moreover, the Kaplan-Meier curve showed that low expression of CBX1 was significantly associated with shorter OS and DFS. Our results are the first to indicate the anti-tumor effect of CBX1 in ccRCC, which is contrary to its pro-tumor effect in some cancers, such as liver cancer.

In 2014, Clermont et al. conducted genotranscriptomic meta-analysis and reported that the expression of CBX2 was higher in many tumors compared to that in normal tissues (39). The most represented cancer types originated from the colon, breast, stomach, and lungs (39). CBX2 overexpression and amplification are closely associated with metastatic progression and shorter OS, especially in breast cancer (40, 41). However, it has also been reported that CBX2 is expressed at similar levels in ccRCC tissues and normal tissues, and the expression of CBX2 is not associated with the prognosis of ccRCC (42). In our study, we found that CBX2 was not related to the expression, tumor characteristics, or prognosis of ccRCC, which is consistent with the results of previous research. Our results showed that CBX2 may not play a role in ccRCC.

Overexpression of CBX3 has also been observed in many malignancies such as lung cancer (43, 44), gastric cancer (36, 45), and tongue squamous cell carcinoma (46, 47). High CBX3 targets CDKN1A to promote the proliferation of U87 cells and predicts poor recurrence-free survival and OS in glioma patients (48). Moreover, CBX3 facilitates the progression of lung adenocarcinoma by directly repressing NCOR2 and ZBTB7A expression (43). However, a study has also shown that the expression of CBX3 can mediate the tumor suppression effect of lncRNA LINC00998 in malignant glioma (49). These studies suggest that CBX3 functions as both an oncogene and a tumor suppressor. In our study, CBX3 expression was increased in ccRCC, and the expression of CBX3 was correlated with advanced cancer stage and tumor grade. In addition, high CBX3 expression was significantly associated with poor OS in ccRCC patients, indicating that CBX3 is an oncogene in ccRCC.

CBX4 is the most studied member of the CBX family. CBX4 has been reported to up-regulate in multiple human tumors including lung adenocarcinoma (50), osteosarcoma (51), breast cancer (52), and cervical cancer (53). And the overexpression of CBX4 has been reported to predict shorter OS in several cancers (54, 55). In ccRCC, CBX4 interacts with HDAC1 to transcriptionally inhibit KLF6, acting as an oncogene with prognostic potential (12). In addition, CBX4 expression can mediate the pro-tumor effect of circTLK1 in RCC (56). However, it has also been reported that CBX4 suppresses cell migration, invasion, and cancer metastasis in colorectal cancer, and high expression of CBX4 is associated with better OS in CRC patients (7). In our study, the expression of CBX4 was upregulated in ccRCC tissues compared to that in normal tissues, and the expression of CBX4 was associated with advanced cancer stage and tumor grade. Moreover, high CBX4 expression was significantly associated with poor OS in ccRCC patients, suggesting that CBX4 promotes ccRCC progression.

CBX5 expression has been found in multiple cancers including breast cancer (57), gastric cancer (58), and lung cancer (59). Upregulation of CBX5 is associated with increased cell proliferation and poor clinical prognosis (59). However, in metastatic colon cancer, thyroid carcinoma, and breast cancer cells, CBX5 expression was found to be downregulated compared to that in poorly invasive or non-metastatic cells (57). Therefore, CBX5 has been characterized as a suppressor of metastasis. In our study, the expression of CBX5 was downregulated in ccRCC tissues compared to that in normal tissues. Low CBX5 levels were associated with advanced cancer stage and tumor grade. In addition, the survival analysis curve showed that low expression of CBX5 was significantly associated with shorter OS in ccRCC patients. Our results revealed the anti-tumor effect of CBX5 in ccRCC.

CBX6 and CBX7 have been found to play contradictory roles in human cancers. The present study suggested that CBX6 and CBX7 could function as both oncogenes and tumor suppressors, depending on the tumor type and cellular context. For example, high CBX6 expression promotes cell growth by regulating the S100A9/NF-κB/MAPK pathway in hepatocellular carcinoma (60). Meanwhile, CBX6 suppressed the progression of breast cancer by significantly downregulating bone marrow stromal cell antigen-2 (61). Similarly, CBX7 acts as a tumor suppressor in several cancers including thyroid cancer (8), colorectal cancer (62). In contrast, in prostate cancer and ovarian cancer, CBX7 is upregulated and acts as a tumor promoter (63, 64). In our study, CBX6 and CBX7 were favorable prognostic factors for ccRCC, and low CBX6 and CBX7 expression was positively associated with advanced cancer stage and tumor grade in ccRCC patients. Our results are consistent with the tumor suppression effects of CBX6 and CBX7.

Cumulative studies have shown that CBX8 can serve as a therapeutic target and a valuable prognostic marker for multiple cancers. In hepatocellular carcinoma, CBX8 upregulates EGR1 and miR-365-3p to stimulate the AKT/β-catenin pathway and shows oncogenic activity (65). Overexpression of CBX8 induces cell proliferation but inhibits cell migration, invasion, and metastasis in esophageal squamous cell carcinoma by repressing Snail (66). Moreover, CBX8 promotes cell proliferation by repressing the p53 pathway and serves as a predictor for muscle-invasive bladder cancer (67). In our study, CBX8 was expressed at similar levels in ccRCC tissues and normal kidney tissues, and CBX8 expression was not associated with advanced tumor characteristics. The survival analysis curve showed that CBX8 expression was an unfavorable prognostic factor for patients with ccRCC.

Several studies have suggested that tumor-infiltrating immune cells can influence tumor progression and recurrence and are closely related to the response to immunotherapy and clinical outcomes of cancer patients (68, 69). In our study, we found a significant relationship between the expression of CBXs and the infiltration of immune cells (B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells), indicating that CBXs might influence the immune status of ccRCC patients.

Conclusions

In conclusion, we firstly systematically investigated the expression and prognostic value of CBXs in ccRCC using an online database. Our results indicated that increased expression of CBX3 and CBX4 promoted ccRCC progression, while CBX1, CBX5, CBX6, and CBX7 were favorable factors in ccRCC. Although more studies should be performed to validate our results, our work provides new insights for selecting therapeutic targets and appropriate drugs for ccRCC patients, and our findings will help determine more accurate predictors for the survival of patients with ccRCC.

Funding

This work was supported by National Natural Science Foundation of China (81602572, 81873574), Natural Science Foundation of Hunan Province, China (2020JJ4905) and Scientific Research Project of Hunan Provincial Health Commission, China (20200489).

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

• 1

  SungHFerlayJSiegelRLLaversanneMSoerjomataramIJemalAet al. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin (2021) 71(3):209–49. doi: 10.3322/caac.21660

  • CrossRef
  • Google Scholar

• 2

  MochH. An Overview of Renal Cell Cancer: Pathology and Genetics. Semin Cancer Biol (2013) 23(1):3–9. doi: 10.1016/j.semcancer

  • CrossRef
  • Google Scholar

• 3

  BarataPCRiniBI. Treatment of Renal Cell Carcinoma: Current Status and Future Directions. CA Cancer J Clin (2017) 67(6):507–24. doi: 10.3322/caac.21411

  • CrossRef
  • Google Scholar

• 4

  GilJO’LoghlenA. PRC1 Complex Diversity: Where Is it Taking Us? Trends Cell Biol (2014) 24(11):632–41. doi: 10.1016/j.tcb.2014.06.005

  • CrossRef
  • Google Scholar

• 5

  WottonDMerrillJC. Pc2 and Sumoylation. Biochem Soc Trans (2007) 35(Pt 6):1401–4. doi: 10.1042/BST0351401

  • CrossRef
  • Google Scholar

• 6

  LiJXuYLongXDWangWJiaoHKMeiZet al. Cbx4 Governs HIF-1α to Potentiate Angiogenesis of Hepatocellular Carcinoma by its SUMO E3 Ligase Activity. Cancer Cell (2014) 25(1):118–31. doi: 10.1016/j.ccr.2013.12.008

  • CrossRef
  • Google Scholar

• 7

  WangXLiLWuYZhangRZhangMLiaoDet al. CBX4 Suppresses Metastasis via Recruitment of HDAC3 to the Runx2 Promoter in Colorectal Carcinoma. Cancer Res (2016) 76(24):7277–89. doi: 10.1158/0008-5472.CAN-16-2100

  • CrossRef
  • Google Scholar

• 8

  PallantePFedericoABerlingieriMTBiancoMFerraroAForzatiFet al. Loss of the CBX7 Gene Expression Correlates With a Highly Malignant Phenotype in Thyroid Cancer. Cancer Res (2008) 68(16):6770–8. doi: 10.1158/0008-5472.CAN-08-0695

  • CrossRef
  • Google Scholar

• 9

  PallantePTerraccianoLCarafaVSchneiderSZlobecILugliAet al. The Loss of the CBX7 Gene Expression Represents an Adverse Prognostic Marker for Survival of Colon Carcinoma Patients. Eur J Cancer (2010) 46(12):2304–13. doi: 10.1016/j.ejca.2010.05.011

  • CrossRef
  • Google Scholar

• 10

  CacciolaNASepeRForzatiFFedericoAPellecchiaSMalapelleUet al. Restoration of CBX7 Expression Increases the Susceptibility of Human Lung Carcinoma Cells to Irinotecan Treatment. Naunyn Schmiedebergs Arch Pharmacol (2015) 388(11):1179–86. doi: 10.1007/s00210-015-1153-y

  • CrossRef
  • Google Scholar

• 11

  ZhangXWZhangLQinWYaoXHZhengLZLiuXet al. Oncogenic Role of the Chromobox Protein CBX7 in Gastric Cancer. J Exp Clin Cancer Res (2010) 29(1):114. doi: 10.1186/1756-9966-29-114

  • CrossRef
  • Google Scholar

• 12

  JiangNNiuGPanYHPanWZhangMFZhangCZet al. CBX4 Transcriptionally Suppresses KLF6 via Interaction With HDAC1 to Exert Oncogenic Activities in Clear Cell Renal Cell Carcinoma. EBioMedicine (2020) 53:102692. doi: 10.1016/j.ebiom.2020.102692

  • CrossRef
  • Google Scholar

• 13

  ChandrashekarDSBashelBBalasubramanyaSAHCreightonCJPonce-RodriguezIChakravarthiBVSKet al. UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses. Neoplasia (2017) 19(8):649–58. doi: 10.1016/j.neo.2017.05.002

  • CrossRef
  • Google Scholar

• 14

  RhodesDRYuJShankerKDeshpandeNVaramballyRGhoshDet al. ONCOMINE: A Cancer Microarray Database and Integrated Data-Mining Platform. Neoplasia (2004) 6(1):1–6. doi: 10.1016/s1476-5586(04)80047-2

  • CrossRef
  • Google Scholar

• 15

  AsplundAEdqvistPHSchwenkJMPonténF. Antibodies for Profiling the Human Proteome-the Human Protein Atlas as a Resource for Cancer Research. Proteomics (2012) 12(13):2067–77. doi: 10.1002/pmic.201100504

  • CrossRef
  • Google Scholar

• 16

  TangZLiCKangBGaoGLiCZhangZ. GEPIA: A Web Server for Cancer and Normal Gene Expression Profiling and Interactive Analyses. Nucleic Acids Res (2017) 45(W1):W98–W102. doi: 10.1093/nar/gkx247

  • CrossRef
  • Google Scholar

• 17

  NagyÁLánczkyAMenyhártOGyőrffyB. Validation of MiRNA Prognostic Power in Hepatocellular Carcinoma Using Expression Data of Independent Datasets. Sci Rep (2018) 8(1):9227. doi: 10.1038/s41598-018-27521-y

  • CrossRef
  • Google Scholar

• 18

  GaoJAksoyBADogrusozUDresdnerGGrossBSumerSOet al. Integrative Analysis of Complex Cancer Genomics and Clinical Profiles Using the Cbioportal. Sci Signal (2013) 6(269):pl1. doi: 10.1126/scisignal.2004088

  • CrossRef
  • Google Scholar

• 19

  LiTFanJWangBTraughNChenQLiuJSet al. TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res (2017) 77(21):e108–10. doi: 10.1158/0008-5472.CAN-17-0307

  • CrossRef
  • Google Scholar

• 20

  SubramanianATamayoPMoothaVKMukherjeeSEbertBLGilletteMAet al. Gene Set Enrichment Analysis: A Knowledge-Based Approach for Interpreting Genome-Wide Expression Profiles. Proc Natl Acad Sci USA (2005) 102(43):15545–50. doi: 10.1073/pnas.0506580102

  • CrossRef
  • Google Scholar

• 21

  BeroukhimRBrunetJPDi NapoliAMertzKDSeeleyAPiresMMet al. Patterns of Gene Expression and Copy-Number Alterations in Von-Hippel Lindau Disease-Associated and Sporadic Clear Cell Carcinoma of the Kidney. Cancer Res (2009) 69(11):4674–81. doi: 10.1158/0008-5472.CAN-09-0146

  • CrossRef
  • Google Scholar

• 22

  JonesJOtuHSpentzosDKoliaSInanMBeeckenWDet al. Gene Signatures of Progression and Metastasis in Renal Cell Cancer. Clin Cancer Res (2005) 11(16):5730–9. doi: 10.1158/1078-0432.CCR-04-2225

  • CrossRef
  • Google Scholar

• 23

  GumzMLZouHKreinestPAChildsACBelmonteLSLeGrandSNet al. Secreted Frizzled-Related Protein 1 Loss Contributes to Tumor Phenotype of Clear Cell Renal Cell Carcinoma. Clin Cancer Res (2007) 13(16):4740–9. doi: 10.1158/1078-0432.CCR-07-0143

  • CrossRef
  • Google Scholar

• 24

  LenburgMELiouLSGerryNPFramptonGMCohenHTChristmanMF. Previously Unidentified Changes in Renal Cell Carcinoma Gene Expression Identified by Parametric Analysis of Microarray Data. BMC Cancer (2003) 3:31. doi: 10.1186/1471-2407-3-31

  • CrossRef
  • Google Scholar

• 25

  VasselliJRShihJHIyengarSRMaranchieJRissJWorrellRet al. Predicting Survival in Patients With Metastatic Kidney Cancer by Gene-Expression Profiling in the Primary Tumor. Proc Natl Acad Sci USA (2003) 100(12):6958–63. doi: 10.1073/pnas.1131754100

  • CrossRef
  • Google Scholar

• 26

  YusenkoMVKuiperRPBoetheTLjungbergBvan KesselAGKovacsG. High-Resolution DNA Copy Number and Gene Expression Analyses Distinguish Chromophobe Renal Cell Carcinomas and Renal Oncocytomas. BMC Cancer (2009) 9:152. doi: 10.1186/1471-2407-9-152

  • CrossRef
  • Google Scholar

• 27

  Álvarez-ErricoDVento-TormoRSiewekeMBallestarE. Epigenetic Control of Myeloid Cell Differentiation, Identity and Function. Nat Rev Immunol (2015) 15(1):7–17. doi: 10.1038/nri3777

  • CrossRef
  • Google Scholar

• 28

  PalSTylerJK. Epigenetics and Aging. Sci Adv (2016) 2(7):e1600584. doi: 10.1126/sciadv.1600584

  • CrossRef
  • Google Scholar

• 29

  KimJH. Chromatin Remodeling and Epigenetic Regulation in Plant DNA Damage Repair. Int J Mol Sci (2019) 20(17):4093. doi: 10.3390/ijms20174093

  • CrossRef
  • Google Scholar

• 30

  CavalliGHeardE. Advances in Epigenetics Link Genetics to the Environment and Disease. Nature (2019) 571(7766):489–99. doi: 10.1038/s41586-019-1411-0

  • CrossRef
  • Google Scholar

• 31

  MaRGZhangYSunTTChengB. Epigenetic Regulation by Polycomb Group Complexes: Focus on Roles of CBX Proteins. J Zhejiang Univ Sci B (2014) 15(5):412–28. doi: 10.1631/jzus.B1400077

  • CrossRef
  • Google Scholar

• 32

  van WijnenAJBagheriLBadreldinAALarsonANDudakovicAThalerRet al. Biological Functions of Chromobox (CBX) Proteins in Stem Cell Self-Renewal, Lineage-Commitment, Cancer and Development. Bone (2021) 143:115659. doi: 10.1016/j.bone.2020.115659

  • CrossRef
  • Google Scholar

• 33

  YangYFPanYHTianQHWuDCSuSG. CBX1 Indicates Poor Outcomes and Exerts Oncogenic Activity in Hepatocellular Carcinoma. Transl Oncol (2018) 11(5):1110–8. doi: 10.1016/j.tranon.2018.07.002

  • CrossRef
  • Google Scholar

• 34

  NingGHuangYLZhenLMXuWXJiaoQYangFJet al. Transcriptional Expressions of Chromobox 1/2/3/6/8 as Independent Indicators for Survivals in Hepatocellular Carcinoma Patients. Aging (Albany NY) (2018) 10(11):3450–73. doi: 10.18632/aging.101658

  • CrossRef
  • Google Scholar

• 35

  MaTMaNChenJLTangFXZongZYuZMet al. Expression and Prognostic Value of Chromobox Family Members in Gastric Cancer. J Gastrointest Oncol (2020) 11(5):983–98. doi: 10.21037/jgo-20-223

  • CrossRef
  • Google Scholar

• 36

  LinKZhuJHuCBuFLuoCZhuXet al. Comprehensive Analysis of the Prognosis for Chromobox Family in Gastric Cancer. J Gastrointest Oncol (2020) 11(5):932–51. doi: 10.21037/jgo-20-208

  • CrossRef
  • Google Scholar

• 37

  LiQPanYCaoZZhaoS. Comprehensive Analysis of Prognostic Value and Immune Infiltration of Chromobox Family Members in Colorectal Cancer. Front Oncol (2020) 10:582667. doi: 10.3389/fonc.2020.582667

  • CrossRef
  • Google Scholar

• 38

  HuAZhangYZhaoXLiJYingY. CBX1 is a Direct Target of miR-205-5p and Contributes to the Progression of Pituitary Tumor. Pharmazie (2019) 74(3):154–6. doi: 10.1691/ph.2019.8908

  • CrossRef
  • Google Scholar

• 39

  ClermontPLSunLCreaFThuKLZhangAParoliaAet al. Genotranscriptomic Meta-Analysis of the Polycomb Gene CBX2 in Human Cancers: Initial Evidence of an Oncogenic Role. Br J Cancer (2014) 111(8):1663–72. doi: 10.1038/bjc.2014.474

  • CrossRef
  • Google Scholar

• 40

  ZhengSLvPSuJMiaoKXuHLiM. Overexpression of CBX2 in Breast Cancer Promotes Tumor Progression Through the PI3K/AKT Signaling Pathway. Am J Transl Res (2019) 11(3):1668–82.

  • Google Scholar

• 41

  IqbalMASiddiquiSUr RehmanAAhmad SiddiquiFSinghPKumarBet al. Multiomics Integrative Analysis Reveals Antagonistic Roles of CBX2 and CBX7 in Metabolic Reprogramming of Breast Cancer. Mol Oncol (2021) 15(5):1450–65. doi: 10.1002/1878-0261.12894

  • CrossRef
  • Google Scholar

• 42

  MengJGaoLZhangMGaoSFanSLiangC. Systematic Investigation of the Prognostic Value of Cell Division Cycle-Associated Proteins for Clear Cell Renal Cell Carcinoma Patients. Biomark Med (2020) 14(3):223–38. doi: 10.2217/bmm-2019-0498

  • CrossRef
  • Google Scholar

• 43

  AlamHLiNDharSSWuSJLvJChenKet al. Hp1γ Promotes Lung Adenocarcinoma by Downregulating the Transcription-Repressive Regulators NCOR2 and ZBTB7A. Cancer Res (2018) 78(14):3834–48. doi: 10.1158/0008-5472.CAN-17-3571

  • CrossRef
  • Google Scholar

• 44

  HanSSKimWJHongYHongSHLeeSJRyuDRet al. RNA Sequencing Identifies Novel Markers of Non-Small Cell Lung Cancer. Lung Cancer (2014) 84(3):229–35. doi: 10.1016/j.lungcan.2014.03.018

  • CrossRef
  • Google Scholar

• 45

  LinHLianJXiaLGuanGYouJ. CBX3 Promotes Gastric Cancer Progression and Affects Factors Related to Immunotherapeutic Responses. Cancer Manag Res (2020) 12:10113–25. doi: 10.2147/CMAR.S271807

  • CrossRef
  • Google Scholar

• 46

  ZhangHChenWFuXSuXYangA. CBX3 Promotes Tumor Proliferation by Regulating G1/s Phase via P21 Downregulation and Associates With Poor Prognosis in Tongue Squamous Cell Carcinoma. Gene (2018) 654:49–56. doi: 10.1016/j.gene.2018.02.043

  • CrossRef
  • Google Scholar

• 47

  LinHZhaoXXiaLLianJYouJ. Clinicopathological and Prognostic Significance of CBX3 Expression in Human Cancer: A Systematic Review and Meta-Analysis. Dis Markers (2020) 2020:2412741. doi: 10.1155/2020/2412741

  • CrossRef
  • Google Scholar

• 48

  ZhaoSPWangFYangMWangXYJinCLJiQKet al. CBX3 Promotes Glioma U87 Cell Proliferation and Predicts an Unfavorable Prognosis. J Neurooncol (2019) 145(1):35–48. doi: 10.1007/s11060-019-03286-w

  • CrossRef
  • Google Scholar

• 49

  CaiHYuYNiXLiCHuYWangJet al. Lncrna LINC00998 Inhibits the Malignant Glioma Phenotype via the CBX3-Mediated C-Met/Akt/mTOR Axis. Cell Death Dis (2020) 11(12):1032. doi: 10.1038/s41419-020-03247-6

  • CrossRef
  • Google Scholar

• 50

  WangZFangZChenGLiuBXuJLiFet al. Chromobox 4 Facilitates Tumorigenesis of Lung Adenocarcinoma Through the Wnt/β-Catenin Pathway. Neoplasia (2021) 23(2):222–33. doi: 10.1016/j.neo.2020.12.005

  • CrossRef
  • Google Scholar

• 51

  WangXQinGLiangXWangWWangZLiaoD. Targeting the CK1α/CBX4 Axis for Metastasis in Osteosarcoma. Nat Commun (2020) 11(1):1141. doi: 10.1038/s41467-020-14870-4

  • CrossRef
  • Google Scholar

• 52

  HuangFJDangJQZhangSChengZY. Circular RNA Hsa_Circ_0008039 Promotes Proliferation, Migration and Invasion of Breast Cancer Cells Through Upregulating CBX4 via Sponging miR-515-5p. Eur Rev Med Pharmacol Sci (2020) 24(4):1887–98. doi: 10.26355/eurrev_202002_20367

  • CrossRef
  • Google Scholar

• 53

  ChenYDuJWangYShiHJiangQWangYet al. MicroRNA-497-5p Induces Cell Cycle Arrest of Cervical Cancer Cells in s Phase by Targeting CBX4. Onco Targets Ther (2019) 12:10535–45. doi: 10.2147/OTT.S210059

  • CrossRef
  • Google Scholar

• 54

  ZengJSZhangZDPeiLBaiZZYangYYangHet al. CBX4 Exhibits Oncogenic Activities in Breast Cancer via Notch1 Signaling. Int J Biochem Cell Biol (2018) 95:1–8. doi: 10.1016/j.biocel.2017.12.006

  • CrossRef
  • Google Scholar

• 55

  JiaoHKXuYLiJWangWMeiZLongXDet al. Prognostic Significance of Cbx4 Expression and Its Beneficial Effect for Transarterial Chemoembolization in Hepatocellular Carcinoma. Cell Death Dis (2015) 6(3):e1689. doi: 10.1038/cddis.2015.57

  • CrossRef
  • Google Scholar

• 56

  LiJHuangCZouYYeJYuJGuiY. CircTLK1 Promotes the Proliferation and Metastasis of Renal Cell Carcinoma by Sponging miR-136-5p. Mol Cancer (2020) 19(1):103. doi: 10.1186/s12943-020-01225-2

  • CrossRef
  • Google Scholar

• 57

  Vad-NielsenJNielsenAL. Beyond the Histone Tale: HP1α Deregulation in Breast Cancer Epigenetics. Cancer Biol Ther (2015) 16(2):189–200. doi: 10.1080/15384047.2014.1001277

  • CrossRef
  • Google Scholar

• 58

  GuoJZhangZPanLZhouY. Identification of miR-758-3p as Potential Modulator of CBX5 Expression in Gastric Cancer. Technol Cancer Res Treat (2018) 17:1533033818816061. doi: 10.1177/1533033818816061

  • CrossRef
  • Google Scholar

• 59

  YuYHChiouGYHuangPILoWLWangCYLuKHet al. Network Biology of Tumor Stem-Like Cells Identified a Regulatory Role of CBX5 in Lung Cancer. Sci Rep (2012) 2:584. doi: 10.1038/srep00584

  • CrossRef
  • Google Scholar

• 60

  ZhengHJiangWHTianTTanHSChenYQiaoGLet al. CBX6 Overexpression Contributes to Tumor Progression and Is Predictive of a Poor Prognosis in Hepatocellular Carcinoma. Oncotarget (2017) 8(12):18872–84. doi: 10.18632/oncotarget.14770

  • CrossRef
  • Google Scholar

• 61

  DengHGuanXGongLZengJZhangHChenMYet al. CBX6 is Negatively Regulated by EZH2 and Plays a Potential Tumor Suppressor Role in Breast Cancer. Sci Rep (2019) 9(1):197. doi: 10.1038/s41598-018-36560-4

  • CrossRef
  • Google Scholar

• 62

  ZhengXZhouJZhangBZhangJWilsonJGuLet al. Critical Evaluation of Cbx7 Downregulation in Primary Colon Carcinomas and Its Clinical Significance in Chinese Patients. BMC Cancer (2015) 15:145. doi: 10.1186/s12885-015-1172-6

  • CrossRef
  • Google Scholar

• 63

  GongLTangYJiangLTangWLuoS. Regulation of Circgolph3 and its Binding Protein CBX7 on the Proliferation and Apoptosis of Prostate Cancer Cells. Biosci Rep (2020) 40(12):BSR20200936. doi: 10.1042/BSR20200936

  • CrossRef
  • Google Scholar

• 64

  ShinjoKYamashitaYYamamotoEAkatsukaSUnoNKamiyaAet al. Expression of Chromobox Homolog 7 (CBX7) Is Associated With Poor Prognosis in Ovarian Clear Cell Adenocarcinoma via TRAIL-Induced Apoptotic Pathway Regulation. Int J Cancer (2014) 135(2):308–18. doi: 10.1002/ijc.28692

  • CrossRef
  • Google Scholar

• 65

  ZhangCZChenSLWangCHHeYFYangXXieDet al. CBX8 Exhibits Oncogenic Activity via AKT/β-Catenin Activation in Hepatocellular Carcinoma. Cancer Res (2018) 78(1):51–63. doi: 10.1158/0008-5472.CAN-17-0700

  • CrossRef
  • Google Scholar

• 66

  WangGTangJZhanWZhangRZhangMLiaoDet al. CBX8 Suppresses Tumor Metastasis via Repressing Snail in Esophageal Squamous Cell Carcinoma. Theranostics (2017) 7(14):3478–88. doi: 10.7150/thno.20717

  • CrossRef
  • Google Scholar

• 67

  YuanGJChenXLuJFengZHChenSLChenRXet al. Chromobox Homolog 8 Is a Predictor of Muscle Invasive Bladder Cancer and Promotes Cell Proliferation by Repressing the P53 Pathway. Cancer Sci (2017) 108(11):2166–75. doi: 10.1111/cas.13383

  • CrossRef
  • Google Scholar

• 68

  WangSSLiuWLyDXuHQuLZhangL. Tumor-Infiltrating B Cells: Their Role and Application in Anti-Tumor Immunity in Lung Cancer. Cell Mol Immunol (2019) 16(1):6–18. doi: 10.1038/s41423-018-0027-x

  • CrossRef
  • Google Scholar

• 69

  LiBSeversonEPignonJCZhaoHLiTNovakJet al. Comprehensive Analyses of Tumor Immunity: Implications for Cancer Immunotherapy. Genome Biol (2016) 17(1):174. doi: 10.1186/s13059-016-1028-7

  • CrossRef
  • Google Scholar