Pan-Cancer Analyses of the Tumor Microenvironment Reveal That Ubiquitin-Conjugating Enzyme E2C Might Be a Potential Immunotherapy Target

Huang, Guang-zhao; Chen, Ze-qun; Wu, Juan; Shao, Ting-ru; Zou, Chen; Ai, Yi-long; Lv, Xiao-zhi

doi:https://doi.org/10.1155/2021/9250207

Journal of Immunology Research

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusions Abbreviations Data Availability Additional Points Ethical Approval Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2021 | Article ID 9250207 | https://doi.org/10.1155/2021/9250207

Pan-Cancer Analyses of the Tumor Microenvironment Reveal That Ubiquitin-Conjugating Enzyme E2C Might Be a Potential Immunotherapy Target

Guang-zhao Huang,¹Ze-qun Chen,²Juan Wu,³Ting-ru Shao,¹Chen Zou,⁴Yi-long Ai,⁴and Xiao-zhi Lv¹

Academic Editor: Xiaofeng Yang

Received21 Sept 2021

Revised28 Oct 2021

Accepted15 Nov 2021

Published13 Dec 2021

Abstract

Increasing evidence indicated that the tumor microenvironment (TME) played a crucial role in cancer initiation and progression. Ubiquitin-conjugating enzyme E2C (UBE2C) was differentially expressed in many cancer types. However, the immunological and prognostic roles of UBE2C were unclear. Differentially expressed genes (DEGs) of 29 cancer types were downloaded from GEPIA2 and 4 cancer types failed to download owing to no DEGs. Furthermore, the gene expression profiles, mutation data, and survival data of 33 cancer types were obtained from UCSC Xena. Clinical stage relevance, tumor mutational burden (TMB), TME relevance analysis, and gene set enrichment analysis (GSEA) of DEGs in 33 cancer types were performed. And DEGs were identified in oral squamous cell carcinoma (OSCC) by biological experiments. Previous studies indicated that UBE2C was related to the prognosis of many cancers. In our study, the higher UBE2C expression level meant a terminal clinical stage in 8 cancer types and the expression level of UBE2C was related to TMB in 20 cancer types. In addition, both immune relevance analysis and GSEA showed that UBE2C might participate in immune response in many cancers. Furthermore, the UBE2C mRNA level and protein level were all identified as upregulated in OSCC cell lines and tissues. UBE2C was differentially expressed in many cancer types and related to the pathogenesis and TME of many cancers, which might be a potential diagnostic and therapeutic biomarker.

1. Introduction

Cancers are severe public health problems all over the world and one of the most common leading causes of death. Annual cases are only expected to increase owing to the aging population [1]. A study showed that an estimated 1,762,450 new cancer cases and 606,880 cancer deaths occurred in the United States in 2019 [2]. In China, there were approximately 4.3 million new cancer cases and 2.9 million new cancer deaths in 2018 [3]. With the increase of various cancer morbidity, it is imperious to understand the potential mechanism of tumorigenesis and tumor progression and identify effective biomarkers for cancer diagnosis and therapy. For the past 20 years, studies have indicated that inflammatory immune cells played a crucial role in cancer-related inflammation [4]. For example, Coussens and Werb identified that during tumor formation, the tissue structure evolved into a highly specialized microenvironment characterized by chronic inflammation [5]. Hanahan and Weinberg revealed that tumor-associated inflammation exists in different stages of tumorigenesis and contributes to genomic instability, epigenetic modification, induction of cancer cell proliferation, enhancement of tumor antiapoptotic pathways, stimulation of angiogenesis, and ultimately tumor dissemination [6]. In addition, patients with Helicobacter pylori were more likely to develop gastric cancer [7]. Human papillomavirus (HPV) infection was associated with head and neck cancer tumorigeneses [8]. Immune dysregulation also increased colorectal cancer incidence [9]. Therefore, the immune system may be relevant to tumorigenesis and tumor progression significantly. Surgery, radiotherapy, and chemotherapy were conventional strategies in the treatment of cancers. In recent years, accumulating studies have shown that immunotherapy is also a targeted systemic option [10–12]. For example, nivolumab may restore antitumor immunity via disrupting PD-1-mediated signaling. Studies also indicated that nivolumab played a crucial role in the treatment of nonsquamous non-small-cell lung cancer and melanoma [12, 13]. In addition, immunotherapies on the basis of PD-1/PDL-1 showed promising clinical responses in hepatocellular carcinoma [14] and lung cancer [15]. Therefore, looking for effective immunotherapeutic targets may contribute to cancer therapy.

Ubiquitin-conjugating enzyme E2C (UBE2C), a crucial member of the E2 family, is involved in the cell cycle process through interacting with the late-promoting complex/annular port (APC/C) and identified as a novel potential tumor biomarker [16]. UBE2C is an integral component of the ubiquitin proteasome system. Studies showed that the expression levels of UBE2C were dysregulated in many cancers with unsatisfied clinical outcomes. For instance, activation of the oncogene UBE2C was concurrently underlying the initiation, progression, and metastasis of lung cancer [17]. Mechanism analysis showed that UBE2C played a significant role in the development and progression of pancreatic ductal adenocarcinoma (PDAC) by regulating cell proliferation and EMT. UBE2C was identified as a novel potential therapeutic target for pancreatic cancer [18]. Recently, increasing evidences demonstrated that immune cell infiltration and the tumor microenvironment (TME) were associated with cancer initiation and progression [6]. A study have showed that UBE2C was related to the prognosis of many cancers [19]. However, the relevance between UBE2C and immune cell infiltration and TME was still unclear.

In our study, differentially expressed genes were screened from GEPIA 2 database (http://gepia2.cancer-pku.cn/#degenes) and UBE2C was identified to be differentially expressed across cancers. Gene expression, mutation, and survival data of 33 cancer types were downloaded from UCSC Xena (https://xena.ucsc.edu/). The purpose of this study was to systematically analyze the clinical relevance and tumor immune microenvironment of UBE2C across 33 cancers. Then, clinical stage and TME relevance analyses were performed. Finally, GSEA was used to annotate the potential function of UBE2C across 33 cancers.

2. Materials and Methods

2.1. Differentially Expressed Gene Acquisition

29 cancer-type differentially expressed genes (DEGs) were downloaded from GEPIA 2 database (http://gepia2.cancer-pku.cn/#degenes) with the cutoff criteria and value < 0.01, and LIMMA was selected as a differential method. PCPG, SARC, MESO, and UVM failed to download owing to no DEGs. And then, 33 cancer-type gene expression data, mutation data, and survival data were further obtained from UCSC Xena (https://xena.ucsc.edu/).

2.2. Univariate Cox Regression Analysis

Overall survival (OS), disease-specific survival (DSS), disease-free interval (DFI), and progression-free interval (PFI) data were gained from UCSC Xena (https://xena.ucsc.edu/) TCGA pan-cancer. Univariate cox regression analyses according to OS, DSS, DFI, and PFI were used to explore the prognostic value of UBE2C in 33 cancer types.

2.3. PPI Network Construction and GO and KEGG Enrichment Analysis

The protein-protein interaction network was constructed in STRING database (https://www.string-db.org/). UBE2C was imported to STRING database to explore the mutual regulation relationship between other genes, and the cutoff of confidence was 0.9. And the top 20 proteins were listed. Furthermore, the proteins enrolled in the PPI network were performed to Gene Ontology (GO), biological processes (BP), cellular component (CC), molecular function (MF), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis in DAVID database (https://david.ncifcrf.gov/).

2.4. Clinical Stage Relevance Analysis

Clinical stage data were also obtained from UCSC Xena (https://xena.ucsc.edu/) TCGA pan-cancer (PANCAN). According to clinical stage, the expression level of UBE2C was compared between stage and stage with the aim at exploring UBE2C expression in the early stage and terminal stage in 33 cancer types.

2.5. Tumor Mutational Burden TMB) Relevance Analysis

TMB is defined as the number of mutations that exist within a tumor and related to the emergence of neoantigens that trigger antitumor immunity. In addition, TMB is also regarded as a new biomarker for prediction of response to immunotherapy [20]. Therefore, TMB data were extracted from mutation files which were downloaded from UCSC Xena (https://xena.ucsc.edu/) and the prognostic value of TMB was explored across 33 cancer types. In addition, relevance analyses between TMB and UBE2C expression levels were also performed.

2.6. TME Relevance Analysis

TME cells consist of a crucial element of tumor tissues. Accumulating evidences demonstrated their clinicopathologic significance in predicting outcomes and therapeutic efficacy [21, 22]. Therefore, the relevance between UBE2C and TME was explored in our study. The “ESTIMATE” package in R software was used to calculate the “immune” and “stromal” scores in 33 cancer types [23]. And the “CIBERSORT” algorithm was performed to assess the 22 immune cell infiltration [23] including naive B cells, memory B cells, plasma cells, CD8 T cells, naive CD4 T cells, resting memory CD4 T cells, memory-activated CD4 T cells, follicular helper T cells, T regulatory cells (Tregs), gamma delta T cells, resting NK cells, activated NK cells, monocytes, M0 macrophages, M1 macrophages, M2 macrophages, resting dendritic cells, activated dendritic cells, resting mast cells, activated mast cells, eosinophils, and neutrophils in 33 cancer types. Finally, the correlation analyses of “immune” and “stromal” scores and immune cell infiltration and UBE2C expression levels were explored in R software.

2.7. Coexpression Analysis

According to gene expression profiles in 33 cancer types, the gene expression levels of immune markers were extracted based on previous research [24–26]. Spearman’s correlation analysis was performed to explore the relevance between UBE2C and these immune-related genes to reveal the potential function of UBE2C in the immune process. And the value and Spearman’s correlation coefficient were calculated. And a heat map was plotted.

2.8. Gene Set Enrichment Analysis (GSEA)

Gene set enrichment analysis (GSEA) derives its power by focusing on gene sets, that is, groups of genes that share common biological function, chromosomal location, and regulation. GSEA was used to investigate the effect of UBE2C on gene enrichment across 33 cancer types in GSEA (version 4.1.0) software.

2.9. Cell Culture

The normal oral epithelial cell line human oral keratinocytes (HOK) and OSCC cell lines SCC9, SCC25, and CAL27 were obtained from the Institute of Antibody Engineering, Southern Medical University (Guangzhou, China). Cell lines HOK, SCC25, and CAL27 were seeded in DMEM and SCC9 in DMEM/F12 containing 10% foetal bovine serum (FBS) (ExCell Bio Inc.) and incubated at 37°C with 5% CO₂.

2.10. Collection of OSCC Tissue Samples

40 cases of OSCC tissue samples and adjacent normal tissues (ANTs) were collected from Nanfang Hospital, Southern Medical University. Respectively, all tumor specimens and ANTs were confirmed as squamous cell carcinoma and normal tissues pathologically. The Nanfang Hospital ethics committee (AF/SC-09/03.2) approved this study.

2.11. RNA Extraction and Quantitative Real-Time PCR (qRT-PCR)

Total RNAs of cell lines and sample tissues were extracted with TRIzol reagent (TaKaRa, Cat# 9109), and the same amount of total RNAs was reversed to cDNA on the basis of the Reverse Transcription Kit according to manufacturer’s protocol (Vazyme). QRT-PCR assay was used to detect the UBE2C mRNA level with ChamQ Universal SYBR qPCR Master Mix (Vazyme Biotech Co. Ltd.). The sequences of the PCR primers were as follows: UBE2C, forward 5-GACCTGAGGTATAAGCTCTCGC-3 and reverse 5-TTACCCTGGGTGTCCACGTT-3; GAPDH, forward 5-CGCTGA GTACGTCGTGGAGTC-3 and reverse 5-GCTGATGATCTT GAGGCTGTTGTC-3.

2.12. Western Blot Assay

The protein of OSCC cell lines and 3 pairs of tissues were extracted by RIPA lysis buffer. Polyacrylamide gel electrophoresis (PAGE) was used to separate protein samples. Furthermore, proteins were transferred to PVDF membranes and then sealed with 5% skim milk. Subsequently, primary antibodies were incubated at 4°C overnight and second antibodies were incubated at room temperature for 1 hour. Finally, the protein level was quantified by ECL. The antibody information was as follows: UBE2C (ABclonal, 1 : 1000), GAPDH (ProteinTech, 1 : 5000), and goat anti-rabbit (ProteinTech, 1 : 10 000).

2.13. Immunohistochemistry

Sample tissues collected from Nanfang Hospital were fixed with 4% formaldehyde, dehydrated, immersed in wax, embedded in paraffin, and cut into 4 μm sections. All tissue sections were dewaxed with xylene and rehydrated in graded ethanol including 100%, 95%, 90%, 80%, 70%, and 50%. Then, endogenous peroxidase was blocked by 3% hydrogen peroxide for 10 mins. Antigen retrieval assay was performed in a pressure cooker for 10–12 minutes with a 0.01 M citrate buffer (pH 6.0). All tissue sections were sealed with 5% BSA. Subsequently, UBE2C antibodies (ABclonal,1 : 100) were incubated at 4°C overnight and secondary antibody was incubated at room temperature for 1 hour. Finally, the sections were visualized with 3,3-diaminobenzidine (DAB). The staining extent was scored from 4 to 0 (76–100%, 26%–5%, 6%–25%, 1%–5%, and 0%). The staining intensity was scored as strong (), weak (), or negative (). A score ranging from 0 to 8 was calculated by multiplying the staining extent and intensity [27].

3. Results

3.1. UBE2C Is Differentially Expressed in Many Cancer Types

By comparison, UBE2C was the only differentially expressed gene in many cancer types (Table 1) according to the GEPIA 2 database. 4 cancer types including PCPG, SARC, MESO, and UVM failed to download owing to no DEGs in GEPIA 2. Survival analysis indicated that patients with a low UBE2C expression level had better prognosis (Figure 1(a)). In addition, Dastsooz et al. explored the prognostic value of UBE2C in various cancers [28]. In our study, univariate Cox regression analysis in 33 cancer types also showed that UBE2C might be a prognostic biomarker in many cancers according to OS, DSS, DFI, and PFI (Figure S1). Furthermore, gene expression profiles were downloaded from UCSC Xena. Differentially expressed analysis indicated that UBE2C was differentially expressed in many types (Figure 1(b)). The value of differential analysis was presented in Supplementary Table 1. There were 9 cancer types including ACC, DLBC, LAML, LGG, MESO, OV, TGCT, UCS, and UVM which were excluded owing to no normal control samples in TCGA database. In addition, SKCM was not performed differentially expressed analysis owing to no normal sample. And THYM had no significance between normal and tumor samples because of high standard deviation. However, UBE2C was upregulated in THYM according to GEPIA 2 database. To sum up, UBE2C was identified to be differentially expressed in many cancer types.

(a)

(b)

3.2. Functional Analysis of Proteins Enrolled in the PPI Network

The top 20 molecules interacting with UBE2C were listed (Figure 2(a)). Furthermore, GO enrichment analysis indicated that these genes were mainly relevant to the proteasome-mediated ubiquitin-dependent protein catabolic process in BP, ubiquitin ligase complex in CC, and histone kinase activity in MF (Figure 2(b)). In addition, KEGG pathway analysis showed that most of proteins enrolled in the PPI network are enriched in the cell cycle and p53 signaling pathway which were associated with cancer development significantly (Figure 2(c)).

(a)

(b)

(c)

3.3. High UBE2C Level Was Related to the Terminal Clinical Stage in Several Cancer Types

Compared early clinical stage to , the UBE2C expression level was significantly increased in the terminal stage in 8 cancer types including ACC (Figure 3(a)), HNSC (Figure 3(b)), KICH (Figure 3(c)), KIRC (Figure 3(d)), KIRP (Figure 3(e)), LIHC (Figure 3(f)), LUSC (Figure 3(g)), and TGCT (Figure 3(h)) indicating that UBE2C might be a terminal biomarker in these cancer types.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

3.4. Prognostic Value Analyses of TMB

According to survival data and the TMB level of each patient, TMB was associated with the survival across 13 cancer types including ACC, BLCA, BRCA, KICH, KIRC, LGG, OV, PCPG, SKCM, STAD, TGCT, THCA, and THYM. Among them, the low TMB level means poor prognosis in BLCA, SKCM, STAD, TGCT, and OV. And in the other 8 cancer types, the lower the TMB level, the better the prognosis in patients (Figures 4(a)–4(m)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

(m)

3.5. UBE2C Was Associated with TMB and Immune Markers

According to TMB data and the UBE2C expression level of each patient, relevance analysis in each cancer type was performed (Table 2). Then, the radar map of TMB was plotted (Figure 5(a)) according to relevance analysis. The results showed that UBE2C was related to the TMB of 20 cancer types including BRCA, PRAD, LUAD, LGG, THYM, STAD, SARC, BLCA, PAAD, LUSC, ACC, SKCM, KICH, COAD, CESC, KIRC, HNSC, MESO, UCEC, and OV. Each gray circle represents the value of the correlation coefficient from −0.7–0.7, and each red spot represents the correlation coefficient between UBE2C and TMB. As shown in Figure 5(a), among the 20 cancer types, UBE2C has only negative associations with THYM and COAD. UBE2C had the strongest correlation with THYM (). And in positive correlation, UBE2C had the strongest correlation with ACC (). Furthermore, the coexpression analysis between UBE2C and the immune marker was performed. According to gene expression profiles in 33 cancer types, UBE2C were related to many immune markers, such as CD200 which is expressed by various cell types, including B cells, a subset of T cells, and thymocytes. And the encoded protein played a significant role in immunosuppression and regulation of antitumor activity. In addition, UBE2C was correlated with CD276 in many cancers. And CD276 was thought to participate in the regulation of T cell-mediated immune response indicating that UBE2C also played a crucial role in immune response. All of the results were shown in Figure 5(b). And the Spearman’s correlation coefficient and value were presented in supplementary Table 2 and supplementary Table 3.

(a)

(b)

Figure 5

Tumor mutational burden (TMB) relevance analysis. (a) UBE2C was associated with the TMB of many cancers including BRCA, PRAD, LUAD, LGG, THYM, STAD, SARC, BLCA, PAAD, LUSC, ACC, SKCM, KICH, COAD, CESC, KIRC, HNSC, MESO, UCEC, and OV. Each gray circle represents the value of the correlation coefficient from −0.7–0.7, and each red spot represents the correlation coefficient between UBE2C and TMB. (b) The heat map of correlation analysis between UBE2C and many immune markers. Blue represents the value, while red stands for Spearman’s correlation coefficient.

3.6. TME Relevance Analysis

The StromalScore, ImmuneScore, and 22 immune cell levels of each patient were calculated in R software. According to the UBE2C expression level, the relevance between StromalScore, ImmuneScore, and UBE2C were explored (Table 3). The results showed that the higher UBE2C expression level, the lower the StromalScore and ImmuneScore in CESC (Figure 6(a)), COAD (Figure 6(b)), GBM (Figure 6(c)), LUSC (Figure 6(d)), PAAD (Figure 6(e)), READ (Figure 6(f)), STAD (Figure 6(g)), and UCEC (Figure 6(h)). In addition, UBE2C was also related to the StromalScore in BRCA, HNSC, LIHC, LUAD, TGCT, THYM, and ImmuneScore in KIRC (Table 3). BRCA (Figure 6(i)) and HNSC (Figure 6(j)) were presented as examples. And the other results of StromalScore and ImmuneScore relevance analysis were displayed in supplementary file1. In these cancers, the high UBE2C expression level meant high ImmuneScore only in KIRC. Furthermore, the relationship between UBE2C expression and 22 immune cells was also explored in R software on the basis of 22 immune cell levels across 33 cancer types and the values were presented in Supplementary Table 4. Majority of cancer types were associated with dysregulation of immune cell levels. For instance, in STAD, PAAD, BRCA, SARC, and LIHC, the higher UBE2C expression level meant lower B cell naïve level. In addition, about 20 cancer types were relevant to resting memory CD4 T cells, such as KIRC which UBE2C, and resting memory CD4 T cells were negatively correlated. At last, all relevance analysis results were visualized in R software (supplementary file2).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

3.7. Gene Set Enrichment Analysis

GSEA was performed with the purpose of exploring the potential function of UBE2C in 33 cancer types. The results showed that UBE2C was related to cell differentiation in a variety of cancers, such as BRCA, UCEC, KIRP, and UVM. In addition, UBE2C participated in immune response in BRCA, CHOL, KIRC, LAML, PAAD, SKCM, THCA, and UCS. UBE2C was associated with the immunoglobulin complex in CHOL and SKCM. In UCS, UBE2C was mainly responsible for immunoglobulin production and regulation of B cell activation. And UBE2C also participated in humoral immune response mediated by circulating immunoglobulin. UBE2C also played a significant role in the metabolic process across BLCA, LIHC, and READ. In MESO, UBE2C was relevant to the cell cycle and MIRNA catabolic process. Besides, UBE2C could promote the cell migration and proliferation in BLCA, CESC, and LUAD. Several cancer types were listed as examples (Figure 7). And the top, 5 GSEA results of other cancer types were presented in supplementary file3.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

3.8. UBE2C Was Upregulated in OSCC Cell Lines and Tissues

First, OSCC data were downloaded from TCGA. Bioinformatics analysis indicated that UBE2C was overexpressed in TCGA (Figure 8(a)), and the higher expression level meant high grade level in OSCC (Figure 8(b)). Furthermore, the expression level of UBE2C was validated in OSCC cell lines and tissues. Our results indicated that UBE2C mRNA was upregulation in OSCC cell lines (Figure 8(c)) and in 40 OSCC patients than ANTs (Figure 8(d)), which were similar with the results in TCGA database. In addition, the protein levels of UBE2C were also upregulated in OSCC cell lines and tissues (Figures 8(e)–8(g)). Moreover, the immunohistochemistry assay also showed that UBE2C was overexpressed in OSCC tissues (Figures 8(h) and 8(i)). In addition, however, the relevance analysis had no significance between the UBE2C expression level and clinical parameters in our study (Table 4).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

Figure 8

UBE2C was upregulated in OSCC. (a) The mRNA level of UBE2C in TCGA OSCC samples (32 normal samples and 319 tumor samples). (b) The high UBE2C expression level meant a high grade level in OSCC. (c, d) The mRNA and protein levels of UBE2C in SCC9 (), SCC25 (), CAL27 (), and OSCC tissues (, ). (e, f) The UBE2C protein level in cell lines and OSCC tissues and adjacent normal tissues. (g) Statistical analysis of the UBE2C protein level in OSCC cell lines. (h) IHC staining based on OSCC tissues. (i) Statistic analysis of IHC assay (, ). T: tumor tissues; ANT: adjacent normal tissues; n: normal adjacent control tissues; IHC: immunohistochemistry.

4. Discussion

The tumorigenesis and development of cancers involve multiple dysregulated processes leading to uncontrolled cell growth. According to the differentially expressed genes obtained from GEPIA2 and gene expression profiles gained from UCSC Xena, UBE2C was identified as a differentially expressed gene in 33 cancer types, which indicated that UCE2C might play an important role in carcinogenesis. The functions of UBE2C are broadly distributed across various cellular processes, including immune response, cell proliferation, migration, cell cycle, and cell differentiation. The constructed PPI network showed that some nodes were identified to play an important role in cancer development. For example, CDK1 was associated with tumor initiation in human melanoma via interacting with SOX2 [29]. In addition, CDK1 also was regarded as a potential prognostic biomarker in lung cancer [30]. Furthermore, GO enrichment analysis demonstrated that these molecules were related to the proteasome-mediated ubiquitin-dependent protein catabolic process and histone kinase activity. And according to KEGG pathway analysis, UBE2C and these genes were enrolled in the PPI network and enriched in the cell cycle and p53 signaling pathway which played a significant role in cancer initiation, progression, and therapy [31, 32]. However, a systematic analysis of the molecular characterization and functional effects of UBE2C on 33 cancer types is still lacking. Presently, accumulating evidences suggested that the immune system played an important role in the tumorigenic process [4]. Therefore, it might be helpful for understanding the potential mechanism of UBE2C in tumorigenesis and screening the potential therapeutic targets for pan-cancer to explore the role of UBE2C in TME.

In this study, we investigated the molecular characterization and TME relevance analyses of UBE2C in more than 10,000 tumor samples of 33 cancer types. The results showed that UBE2C was a differentially expressed genes in 33 cancer types. Recently, studies also showed that UBE2C was upregulated and promoted epithelial-mesenchymal transition via p53 in endometrial cancer [33]. In hepatocellular carcinoma, high UBE2C expression induced by overexpression of DNA primase subunit 1 (PRIM1) also could cause the ubiquitination and degradation of p53 [34]. Recently, a study has shown that not only did p53 become dysfunctional in most cancers but mutant p53 also acquired dominant negative activity and carcinogenic properties. And p53 remains an attractive therapeutic target for cancer therapy [35] indicating that UBE2C might be a promising therapeutic biomarker. Zhang et al. also demonstrated that the high UBE2C expression level resulted in chromosomal instability that disturbed the cell cycle and led to poor prognosis of intestinal-type gastric cancer [36]. UBE2C was regarded as an important marker of chromosomal instability and has been associated with malignant growth [37]. In addition, growing researches suggested that UBE2C was dysregulated in prostate cancer [38], cervical cancer [39], breast cancer [40], and so on. In addition, Jin et al. demonstrated that UBE2C could promote the progression of head and neck squamous cell carcinoma (HNSCC) [19]. High UBE2C expression tended to lymph node metastasis in HNSCC. Our study also indicated that UBE2C was upregulated in OSCC cell lines and tissues by RT-qPCR and IHC assay. However, due to the small sample size of this experiment, our results showed that the UBE2C expression level in OSCC had no significant relationship with the clinical parameters. In total, UBE2C could be remarkably associated with cancer development and a potential biomarker for pan-cancer. Recently, a study has demonstrated that UBE2C upregulation was relevant to the advanced histologic grade, FIGO stage, and recurrence and might be a new biomarker for the diagnosis and therapy of endometrial cancer [33]. Studies also showed that UBE2C played a crucial role in hepatocellular carcinoma [41, 42]. For instance, the lower UBE2C expression level indicated higher sensitivity for the therapy of chemotherapeutic drug, including adriamycin and 5-fluorouracil, and knockdown of UBE2C also increased the sensitivity of hepatocellular carcinoma cell to sorafenib [41]. UBE2C, as an oncogene in these cancers, played a vital role in cancer development and therapy which indicated that it might be a potential effective therapeutic target for pan-cancer. Furthermore, our study showed that TMB, a new biomarker for the prediction of immunotherapeutic response, was related to the prognosis in 13 cancer types. And the differences in prognosis between cancers may be due to tumor heterogeneity. In addition, the UBE2C expression level was remarkably related to TMB. A study demonstrated that TMB was also a useful biomarker for immune checkpoint blockade (ICB) selection across some cancer types [43]. Therefore, UBE2C might be associated with immunotherapy and might be an immunotherapeutic target in these 20 cancer types. Furthermore, TME relevance analysis also indicated that the higher UBE2C expression level meant lower StromalScore and ImmuneScore levels in many cancer types like PAAD, UCEC, and BRCA. Microenvironment-mediated drug resistance can be induced by soluble factors secreted by tumor or stromal cells [44], and the TME was significantly related to the therapeutic response and clinical outcome. Recently, a study suggested that IL-17 promoted the progression of invasive prostate adenocarcinomas under castrate conditions via creating an immunotolerant and proangiogenic TME [45]. Castrate mice without IL-17 receptor C (IL-17RC) prostate had lower rates of cellular proliferation as well as lower UBE2C protein level, indicating that UBE2C might play a role in TME of prostate cancer. In addition, our study showed that UBE2C was significantly relevant to TME in many cancer types and the TME played a crucial role in tumor initiation, progression, and metastasis process and also threw light on therapeutic efficacy [44] suggesting that UBE2C might a potential immunotherapeutic target. Unfortunately, there are no other reports about the relationship between UBE2C and TME.

At last, GSEA results showed that UBE2C was associated with the immunoglobulin complex in CHOL and SKCM and participated in the humoral immune response mediated by circulating immunoglobulin in LAML and immunoglobulin production and regulation of B cell activation in UCS. Moreover, coexpression analysis also showed that the UBE2C expression level was associated with many immune markers, such as CD200, CD276, and TMIGD2. These immune markers were related to immunosuppression and regulation of antitumor activity and participated in the regulation of T cell-mediated immune response indicating that UBE2C played a vital role in the immune response in a majority of cancer types. For example, Luo et al. demonstrated that UBE2C was involved in receptor ligand activity and the hormone activity high UBE2C expression profile meant worse biochemical recurrence-free survival in prostate cancer [46]. In addition, UBE2C might be involved in B cell receptor signaling pathways and other functions associated with immune system activation in nasopharyngeal carcinoma [47]. These results threw light on the dysregulation of UBE2C in cancers and uncover that UBE2C might be a potential diagnostic and therapeutic biomarker for various cancers.

5. Conclusions

In our study, UBE2C was identified to be differentially expressed in many cancer types and related to the initiation, progression, and prognosis of many cancers. TME relevance analysis indicated that UBE2C also played a crucial role in the immune response of various cancers and the results are consistent with GSEA results. In summary, the systematic analysis of UBE2C revealed that it might be a diagnostic and therapeutic biomarker across many cancers.

Though our study might have vital clinical significance, there are a few limitations in this study. First, UBE2C has not been identified as an effective biomarker in a variety of cancer samples. Second, the potential mechanism of UBE2C promoting some cancer initiation and progression should be investigated in further study. Third, the relevance between UBE2C and immune response needs further investigation.

Abbreviations

DEGs:	Differentially expressed genes
OSCC:	Oral squamous cell carcinoma
UBE2C:	Ubiquitin-conjugating enzyme E2C
TCGA:	The Cancer Genome Atlas
DEGs:	Differentially expressed genes
GSEA:	Gene set enrichment analysis
OS:	Overall survival
DSS:	Disease-specific survival
DFI:	Disease-free interval
PFI:	Progression-free interval
PPI:	Protein-protein interaction
TMB:	Tumor mutational burden
TME:	Tumor microenvironment
GO:	Gene Ontology
KEGG:	Kyoto Encyclopedia of Genes and Genomes
BP:	Biological processes
CC:	Cellular component
MF:	Molecular function
ANTs:	Adjacent normal tissues.

Data Availability

All data for this study are available from the corresponding authors if required.

Additional Points

Statement. A preprint has previously been published according to this manuscript [48].

Ethical Approval

The Nanfang Hospital ethics committee (AF/SC-09/03.2) approved this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Data curation was performed by Guang-zhao Huang, Ze-qun Chen, and Juan Wu; formal analysis was performed by Guang-zhao Huang and Juan Wu; the methodology was done by Guang-zhao Huang, Juan Wu, and Ting-ru Shao; the investigation was done by Ting-ru Shao and Chen Zou; project administration was performed by Yi-long Ai and Xiao-zhi Lv; supervision was done by Yi-long Ai and Xiao-zhi Lv; funding acquisition was performed by Xiao-zhi Lv; writing of the original draft and the review & editing of the writing were done by Guang-zhao Huang. Guang-zhao Huang, Ze-qun Chen, and Juan Wu contributed equally to this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (81472536), the Science and Technology Planning Project of Guangdong Province (2017A020215181 and 2014A020212440), the Project of Educational Commission of Guangdong Province of China (2018KTSCX026), the Scientific Research Planning Project of Southern Medical University (CX2018N016), and the Presidential Foundation of the Nanfang Hospital (2014027 and 2019Z030).

Supplementary Materials

Supplementary 1. Figure S1: univariate Cox regression analysis of UBE2C according to OS (A), DSS (B), DFI (C), and PFI (D).

Supplementary 2. Supplementary Table1: the value of differentially expressed analysis in Figure 1(b). Supplementary Table2: the Spearman’s correlation coefficient between UBE2C and immune markers. Supplementary Table3: the value of correlation analysis between UBE2C and immune markers. Supplementary Table4: the value between UBE2C expression and 22 immune cells correlation analysis.

Supplementary 3. Supplementary file1: tumor microenvironment (TME) relevance analysis was listed. Supplementary file2: visualization of relevance analysis between UBE2C expression and 22 immune cell levels. Supplementary file3: the gene set enrichment analysis (GSEA) results of the other 23 cancer types.

References

J. Ferlay, I. Soerjomataram, R. Dikshit et al., “Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012,” International Journal of Cancer, vol. 136, no. 5, pp. E359–E386, 2015.
View at: Publisher Site | Google Scholar
R. L. Siegel, K. D. Miller, and A. Jemal, “Cancer statistics, 2019,” CA: A Cancer Journal for Clinicians, vol. 69, no. 1, pp. 7–34, 2019.
View at: Publisher Site | Google Scholar
R. M. Feng, Y. N. Zong, S. M. Cao, and R. H. Xu, “Current cancer situation in China: good or bad news from the 2018 Global Cancer Statistics?” Cancer Commun (Lond)., vol. 39, no. 1, p. 22, 2019.
View at: Publisher Site | Google Scholar
H. Gonzalez, C. Hagerling, and Z. Werb, “Roles of the immune system in cancer: from tumor initiation to metastatic progression,” Genes & Development, vol. 32, no. 19-20, pp. 1267–1284, 2018.
View at: Publisher Site | Google Scholar
L. M. Coussens and Z. Werb, “Inflammation and cancer,” Nature, vol. 420, no. 6917, pp. 860–867, 2002.
View at: Publisher Site | Google Scholar
D. Hanahan and R. A. Weinberg, “Hallmarks of cancer: the next generation,” Cell, vol. 144, no. 5, pp. 646–674, 2011.
View at: Publisher Site | Google Scholar
S. P. Hussain, P. Amstad, K. Raja et al., “Increased p53 mutation load in noncancerous colon tissue from ulcerative colitis: a cancer-prone chronic inflammatory disease,” Cancer Research, vol. 60, no. 13, pp. 3333–3337, 2000.
View at: Google Scholar
C. R. Leemans, B. J. Braakhuis, and R. H. Brakenhoff, “The molecular biology of head and neck cancer,” Nature Reviews. Cancer, vol. 11, no. 1, pp. 9–22, 2011.
View at: Publisher Site | Google Scholar
P. L. Lakatos and L. Lakatos, “Risk for colorectal cancer in ulcerative colitis: changes, causes and management strategies,” World Journal of Gastroenterology, vol. 14, no. 25, pp. 3937–3947, 2008.
View at: Publisher Site | Google Scholar
J. R. Brahmer, S. S. Tykodi, L. Q. Chow et al., “Safety and activity of anti-PD-L1 antibody in patients with advanced cancer,” The New England Journal of Medicine, vol. 366, no. 26, pp. 2455–2465, 2012.
View at: Publisher Site | Google Scholar
P. Sharma, M. K. Callahan, P. Bono et al., “Nivolumab monotherapy in recurrent metastatic urothelial carcinoma (CheckMate 032): a multicentre, open-label, two-stage, multi-arm, phase 1/2 trial,” The Lancet Oncology, vol. 17, no. 11, pp. 1590–1598, 2016.
View at: Publisher Site | Google Scholar
H. Borghaei, L. Paz-Ares, L. Horn et al., “Nivolumab versus docetaxel in advanced nonsquamous non-small-cell lung cancer,” The New England Journal of Medicine, vol. 373, no. 17, pp. 1627–1639, 2015.
View at: Publisher Site | Google Scholar
C. Robert, G. V. Long, B. Brady et al., “Nivolumab in previously untreated melanoma withoutBRAFmutation,” The New England Journal of Medicine, vol. 372, no. 4, pp. 320–330, 2015.
View at: Publisher Site | Google Scholar
J. J. Harding, D. I. El, and G. K. Abou-Alfa, “Immunotherapy in hepatocellular carcinoma: primed to make a difference?” Cancer, vol. 122, no. 3, pp. 367–377, 2016.
View at: Publisher Site | Google Scholar
A. Bocanegra, G. Fernandez-Hinojal, M. Zuazo-Ibarra et al., “PD-L1 expression in systemic immune cell populations as a potential predictive biomarker of responses to PD-L1/PD-1 blockade therapy in lung cancer,” International Journal of Molecular Sciences, vol. 20, no. 7, p. 1631, 2019.
View at: Publisher Site | Google Scholar
C. Xie, C. Powell, M. Yao, J. Wu, and Q. Dong, “Ubiquitin-conjugating enzyme E2C: a potential cancer biomarker,” The International Journal of Biochemistry & Cell Biology, vol. 47, pp. 113–117, 2014.
View at: Publisher Site | Google Scholar
J. Guo, Y. Wu, J. du et al., “Deregulation of UBE2C-mediated autophagy repression aggravates NSCLC progression,” Oncogene, vol. 7, no. 6, p. 49, 2018.
View at: Publisher Site | Google Scholar
X. Wang, L. Yin, L. Yang et al., “Silencing ubiquitin-conjugating enzyme 2C inhibits proliferation and epithelial-mesenchymal transition in pancreatic ductal adenocarcinoma,” The FEBS Journal, vol. 286, no. 24, pp. 4889–4909, 2019.
View at: Publisher Site | Google Scholar
Z. Jin, X. Zhao, L. Cui et al., “UBE2C promotes the progression of head and neck squamous cell carcinoma,” Biochemical and Biophysical Research Communications, vol. 523, no. 2, pp. 389–397, 2020.
View at: Publisher Site | Google Scholar
M. Allgäuer, J. Budczies, P. Christopoulos et al., “Implementing tumor mutational burden (TMB) analysis in routine diagnostics-a primer for molecular pathologists and clinicians,” Transl Lung Cancer Res., vol. 7, no. 5, pp. 703–715, 2018.
View at: Publisher Site | Google Scholar
D. Zeng, R. Zhou, Y. Yu et al., “Gene expression profiles for a prognostic immunoscore in gastric cancer,” The British Journal of Surgery, vol. 105, no. 10, pp. 1338–1348, 2018.
View at: Publisher Site | Google Scholar
Y. Jiang, Q. Zhang, Y. Hu et al., “ImmunoScore signature: a prognostic and predictive tool in gastric cancer,” Annals of Surgery, vol. 267, no. 3, pp. 504–513, 2018.
View at: Publisher Site | Google Scholar
N. Zhu and J. Hou, “Assessing immune infiltration and the tumor microenvironment for the diagnosis and prognosis of sarcoma,” Cancer Cell International, vol. 20, no. 1, p. 577, 2020.
View at: Publisher Site | Google Scholar
P. Danaher, S. Warren, L. Dennis et al., “Gene expression markers of tumor infiltrating leukocytes,” Journal for Immunotherapy of Cancer, vol. 5, no. 1, p. 18, 2017.
View at: Publisher Site | Google Scholar
J. H. Pan, H. Zhou, L. Cooper et al., “LAYN is a prognostic biomarker and correlated with immune infiltrates in gastric and colon cancers,” Frontiers in Immunology, vol. 10, p. 6, 2019.
View at: Publisher Site | Google Scholar
N. O. Siemers, J. L. Holloway, H. Chang et al., “Genome-wide association analysis identifies genetic correlates of immune infiltrates in solid tumors,” PLoS One, vol. 12, no. 7, article e179726, 2017.
View at: Publisher Site | Google Scholar
G. Z. Huang, Z. Y. Lu, Y. Rao, H. Gao, and X. Z. Lv, “Screening and identification of autophagy-related biomarkers for oral squamous cell carcinoma (OSCC) via integrated bioinformatics analysis,” Journal of Cellular and Molecular Medicine, vol. 25, no. 9, pp. 4444–4454, 2021.
View at: Publisher Site | Google Scholar
H. Dastsooz, M. Cereda, D. Donna, and S. Oliviero, “A comprehensive bioinformatics analysis of UBE2C in cancers,” International Journal of Molecular Sciences, vol. 20, no. 9, p. 2228, 2019.
View at: Publisher Site | Google Scholar
D. Ravindran Menon, Y. Luo, J. J. Arcaroli et al., “CDK1 interacts with Sox2 and promotes tumor initiation in human melanoma,” Cancer Research, vol. 78, no. 23, pp. 6561–6574, 2018.
View at: Publisher Site | Google Scholar
G. I. Evan and K. H. Vousden, “Proliferation, cell cycle and apoptosis in cancer,” Nature, vol. 411, no. 6835, pp. 342–348, 2001.
View at: Publisher Site | Google Scholar
M. Li, F. He, Z. Zhang, Z. Xiang, and D. Hu, “CDK1 serves as a potential prognostic biomarker and target for lung cancer,” The Journal of International Medical Research, vol. 48, no. 2, p. 030006051989750, 2020.
View at: Publisher Site | Google Scholar
A. H. Stegh, “Targeting the p53 signaling pathway in cancer therapy - the promises, challenges and perils,” Expert Opinion on Therapeutic Targets, vol. 16, no. 1, pp. 67–83, 2012.
View at: Publisher Site | Google Scholar
Y. Liu, R. Zhao, S. Chi et al., “UBE2CIs upregulated by estrogen and promotes epithelial-mesenchymal transition via p53 in endometrial cancer,” Molecular Cancer Research, vol. 18, no. 2, pp. 204–215, 2020.
View at: Publisher Site | Google Scholar
M. Zhu, M. Wu, S. Bian et al., “DNA primase subunit 1 deteriorated progression of hepatocellular carcinoma by activating AKT/mTOR signaling and UBE2C-mediated P53 ubiquitination,” Cell & Bioscience, vol. 11, no. 1, p. 42, 2021.
View at: Publisher Site | Google Scholar
B. Hong, A. P. van den Heuvel, V. V. Prabhu, S. Zhang, and W. el-Deiry, “Targeting tumor suppressor p53 for cancer therapy: strategies, challenges and opportunities,” Current Drug Targets, vol. 15, no. 1, pp. 80–89, 2014.
View at: Publisher Site | Google Scholar
J. Zhang, X. Liu, G. Yu et al., “UBE2C is a potential biomarker of intestinal-type gastric cancer with chromosomal instability,” Frontiers in Pharmacology, vol. 9, p. 847, 2018.
View at: Publisher Site | Google Scholar
S. Bajaj, S. K. Alam, K. S. Roy, A. Datta, S. Nath, and S. Roychoudhury, “E2 ubiquitin-conjugating enzyme, UBE2C gene, is reciprocally regulated by wild-type and gain-of-function mutant p53,” The Journal of Biological Chemistry, vol. 291, no. 27, pp. 14231–14247, 2016.
View at: Publisher Site | Google Scholar
J. Hu, X. Wu, C. Yang et al., “Anticancer effect of icaritin on prostate cancer via regulating miR-381-3p and its target gene UBE2C,” Cancer Medicine, vol. 8, no. 18, pp. 7833–7845, 2019.
View at: Publisher Site | Google Scholar
A. J. Chiang, C. J. Li, K. H. Tsui et al., “UBE2C drives human cervical cancer progression and is positively modulated by mTOR,” Biomolecules, vol. 11, no. 1, p. 37, 2021.
View at: Publisher Site | Google Scholar
T. Qin, G. Huang, L. Chi et al., “Exceptionally high UBE2C expression is a unique phenomenon in basal-like type breast cancer and is regulated by BRCA1,” Biomedicine & Pharmacotherapy, vol. 95, pp. 649–655, 2017.
View at: Publisher Site | Google Scholar
Y. Xiong, J. Lu, Q. Fang et al., “UBE2C functions as a potential oncogene by enhancing cell proliferation, migration, invasion, and drug resistance in hepatocellular carcinoma cells,” Bioscience Reports, vol. 39, no. 4, 2019.
View at: Publisher Site | Google Scholar
G. Zhang, Z. Kang, H. Mei, Z. Huang, and H. Li, “Promising diagnostic and prognostic value of six genes in human hepatocellular carcinoma,” American Journal of Translational Research, vol. 12, no. 4, pp. 1239–1254, 2020.
View at: Google Scholar
T. A. Chan, M. Yarchoan, E. Jaffee et al., “Development of tumor mutation burden as an immunotherapy biomarker: utility for the oncology clinic,” Annals of Oncology, vol. 30, no. 1, pp. 44–56, 2019.
View at: Publisher Site | Google Scholar
T. Wu and Y. Dai, “Tumor microenvironment and therapeutic response,” Cancer Letters., vol. 387, pp. 61–68, 2017.
View at: Publisher Site | Google Scholar
Q. Zhang, S. Liu, Q. Zhang et al., “Interleukin-17 promotes development of castration-resistant prostate cancer potentially through creating an immunotolerant and pro-angiogenic tumor microenvironment,” The Prostate, vol. 74, no. 8, pp. 869–879, 2014.
View at: Publisher Site | Google Scholar
C. Luo, J. Chen, and L. Chen, “Exploration of gene expression profiles and immune microenvironment between high and low tumor mutation burden groups in prostate cancer,” International Immunopharmacology, vol. 86, 2020.
View at: Publisher Site | Google Scholar
W. Fang, X. Li, Q. Jiang et al., “Transcriptional patterns, biomarkers and pathways characterizing nasopharyngeal carcinoma of Southern China,” Journal of Translational Medicine, vol. 6, no. 1, p. 32, 2008.
View at: Publisher Site | Google Scholar
G. Huang, Z. Q. Chen, T. R. Shao et al., Ubiquitin conjugating enzyme E2C correlates with tumor microenvironment and impacts with patient prognosis: a pan-cancer analysis, https://papers.ssrn.com/sol3/papers.cfm?abstract_id=3823129.

Copyright

Copyright © 2021 Guang-zhao Huang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

488

Downloads

701

Citations