Skip to main content

Prognostic value of ubiquitin E2 UBE2W and its correlation with tumor-infiltrating immune cells in breast cancer

Abstract

Background

Ubiquitin-conjugating enzyme E2W (UBE2W) is a protein-coding gene that has an important role in ubiquitination and may be vital in the repair of DNA damage. However, studies on the prognostic value of UBE2W and its correlation with tumor-infiltrating immune cells in multiple cancers have not been addressed.

Methods

We investigated UBE2W expression in the Oncomine database, the Tumor Immune Estimation Resource (TIMER), TNMplot database. Then, the clinical prognostic value of UBE2W was analyzed via online public databases. Meanwhile, we explored the correlation between UBE2W and DNA repair associate genes expression and DNA methyltransferase expression by TIMER and Gene Expression Profiling Interactive Analysis (GEPIA). By using the same method, the correlation between UBE2W and tumor-infiltrating immune cells was explored. Genomic Profiles of UBE2W in breast cancer (BRCA) were accessed in cBioPortal (v3.5.0). Functional proteins associated network was analyzed by STRING database (v11.0).

Results

UBE2W was abnormally expressed and significantly correlated with mismatch repair (MMR) gene mutation levels, DNA methyltransferase, and BRCA1/2 expression in breast cancer. High expression of UBE2W may promote the tumor immunosuppression and metastasis, induce endocrine therapy resistance and deteriorate outcomes of patients with breast cancer. These findings suggest that UBE2W could be a potential biomarker of prognosis and tumor-infiltrating immune cells. Besides, RBX1 may be a new E3 that was regulated by UBE2W.

Conclusions

Ubiquitin E2 UBE2W is a potential prognostic biomarker and is correlated with immune infiltration in BRCA.

Peer Review reports

Introduction

Cancer is a major public health problem with increasing rapid incidence and mortality in the worldwide [1]. Over the last few decades, the various molecular mechanism of carcinogenesis has been identified, including gene mutation and epigenetic modification [2]. As an important post-translational modification, ubiquitin could participate in the oncogenesis and tumor resistance through protein degradation pathways [3].

Ubiquitin–proteasome system (UPS) mediated 80% ~ 85% protein degradation in eukaryotes [4]. The system regulates a variety of cellular activities including cell cycle, apoptosis, transcriptional regulation, DNA repair, and immune response [5, 6]. Recent findings indicate that UPS plays a key role in the carcinogenesis and progression of different types of tumors [7, 8]. The natural corollary of the above is the prediction that UPS components may be major contributory factors to human cancer. Ubiquitin-conjugating enzyme E2W (UBE2W) is an E2 ubiquitin-conjugating enzyme with direct protein N-terminal monoubiquitylation activity [9]. UBE2W could participate in the regulation of DNA repair via different pathways [10, 11]. A few recent studies disclosed that loss of ubiquitin E2 UBE2W rescues hypersensitivity of Rnf4 mutant cells to DNA damage [12]. Moreover, UBE2W was confirmed directly to participate in the ubiquitination of BRCA1/BRAD1 [13]. However, the biology and potential functions of UBE2W in tumorigenesis and mechanism are poorly understood. We speculated that UBE2W might be involved in tumorigenesis by ubiquitinating substrate N-terminal.

Additionally, UPS exerts the ability to regulate protein stability and function which is critical for the regulation of immunity [14]. For instance, ubiquitin-protein ligase E3 component N-recognin 5 (UBR5) could promote tumor proliferation by inhibiting the immune response by decreasing the CD8+ T cell [15]. Ubiquitin-conjugating enzyme UBE2D3 was also reported to regulate immunity via the NF-kB signaling pathway [16]. Thus, we speculated UBE2W may be involved in the tumor immune microenvironment.

In the present study, we primarily focus on the association between UBE2W expression and patients’ prognosis in different cancer types. The correlation between UBE2W expression and the immune environment was also explored.

Methods

Gene expression analysis

UBE2W mRNA expression levels were analyzed through the Oncomine database (https://www.oncomine.org) [17]. The threshold settings were as follows: gene ranking of the top 10%, fold change of 2.0, and P-value of 0.01. The Tumor Immune Estimation Resource (TIMER) database (https://timer.comp-genomics) and TNMplot database (http://www.tnmplot.com) were used to verify the expression of UBE2W in various cancers [18,19,20,21]. TNMplot database (http://www.tnmplot.com) with RNA-sequence data and gene-chip data were used to explore the UBE2W mRNA expression in tumor, normal, and metastatic tissues. Immunohistochemical staining pictures of UBE2W protein expression were obtained from the Human Protein Atlas (HPA) database (https://www.proteinatlas.org/) [22].

Survival and prognosis analysis

The association of UBE2W with prognosis in different tumors via PrognoScan database (http://dna.bio.kyutech.ac.jp/PrognoScan) and Kaplan-Meier-plotter (KM plotter) database (http://kmplot.com/analysis/) [23, 24]. Subgroups survival analysis was provided in Kaplan-Meier survival plots. P-value < 0.05 was considered statistically significant. The ROC plotter (http://www.rocplot.org) was used to predict the UBE2W expression and response to therapy in 3104 BRCA patients [25].

Correlation analysis between UBE2W expression and immune characteristics

The correlation between the expression of UBE2W with the abundances of six tumor-infiltrating immune cells (TIICs) (B cells, CD4+ T cells, CD8+ T cells, macrophages, neutrophils, and dendritic cells) was explored via the TIMER database. Sixty-six related gene markers of TIICs derived from the CellMarker database (http://biocc.hrbmu.edu.cn/CellMarker/) were used for the analysis [26]. Gene Expression Profiling Interactive Analysis (GEPIA) (http://gepia.cancer-pku.cn/index.html) was used to validate the gene correlation analysis in TIMER [27].

Correlation analysis between UBE2W and BRCA1/2 genes, MMR genes, DNA methyltransferase analysis

Both the TIMER database and GEPIA database was used to evaluate the relationship between UBE2W expression and DNA repair gene (BRCA1/2 genes, five MMR genes), DNA methyltransferase genes. MMR genes contain MLH1, MSH2, MSH6, PMS2, and EPCAM. Meanwhile, DNA methyltransferase genes include DNMT1, DNMT2, DNMT3A, and DNMT3B.

Genomic analysis

cBioPortal for Cancer Genomics (v3.5.0) (http://www.cbioportal.org/) is a powerful platform that contains DNA copy number data, mRNA and microRNA expression data, non-synonymous mutations, protein level, and clinical data of various tumors [28, 29]. The genomics features of UBE2W on a specific cancer type were analyzed via the cBioPortal database. Besides, the association between the CNV alteration of UBE2W with survival was also displayed.

Functional protein association network analysis

STRING database (https://string-db.org/) is a powerful tool for the analysis of protein interactions involving 2031 species, 9,643,763 proteins, and 1,380,838,440 interactions [30]. The co-expression protein of UBE2W was obtained from the cBioPortal database analysis. The proteins network was established via the STRING database.

Statistical analysis

Analysis of UBE2W expression differences based on different databases, so the test methods are different. In the Oncomine database, Student’s t-test was used to analyze the two-class differential expression. The statistical significance was computed by the Wilcoxon test in the TIMER database. While the significant differences in the TNM plot database were analyzed by the Mann-Whitney U test. Survival analysis with p values was based on a log-rank test. Spearman’s correlation coefficient and the P-value were used to evaluate the gene correlation. P-value < 0.05 was considered statistically significant. The flow diagram is shown in Supplementary Figure 1.

Results

The high mRNA expression level of UBE2W in breast Cancer

First, we analyze the expression level of UBE2W via the Oncomine database (Fig. 1a, p-value = 0.05, fold change = 2). Several datasets showed high expression levels of UBE2W was significantly in cancer samples verse in normal samples, including breast cancer (BRCA), head and neck cancers (HNSC), lymphoma, colorectal cancer, melanoma (Supplementary Table 1). Next, we further investigate the mRNA expression levels of UBE2W in The Cancer Genome Atlas (TCGA) based on TIMER database (Fig. 1b). The expression of UBE2W was absolutely higher in breast cancer than normal tissues similarly. The higher mRNA expression level of UBE2W was also detected in the cholangiocarcinoma (CHOL), colon adenocarcinoma (COAD), esophageal carcinoma (ESCA), head and neck tumor (HNSC), liver hepatocellular carcinoma (LIHC), lung adenocarcinoma (LUAD), lung squamous cell carcinoma (LUSC), stomach adenocarcinoma (STAD) than in the normal tissues. In addition, TNMplot database with RNA-seq data from TCGA, Genotype-Tissue Expression GTEX, Therapeutically Applicable Research to Generate Effective Treatments (TARGET) were used to verify the expression of UBE2W (Fig. 1c). Taken together, these results strongly demonstrated the UBE2W gene was abnormally regulated in multiple cancers than matched normal tissues, especially in breast cancer.

Fig. 1
figure1

Expression of UBE2W in various human tumors. a Increased expression of UBE2W in different tumors compared to normal tissues in the Oncomine database. b UBE2W expression of different tumor types from the TCGA database was investigated by TIMER (*P < 0.05, **P < 0.01, ***P < 0.001). c UBE2W was abnormally expressed in pan cancers by TNMplot. Significant differences by Mann-Whitney U test are marked with red color. Abbreviations: ACC: adrenocortical carcinoma; BLCA: bladder urothelial carcinoma; BRCA: breast invasive carcinoma; CESC: cervical squamous cell carcinoma; CHOL: cholangiocarcinoma; COAD: colon adenocarcinoma; DLBC: lymphoid neoplasm diffuse large B cell lymphoma; ESCA: esophageal carcinoma; GBM: glioblastoma multiforme; LGG: brain lower grade glioma; HNSC: head and neck squamous cell carcinoma; KICH: kidney chromophobe; KIRC: kidney renal clear cell carcinoma; KIRP: kidney renal papillary cell carcinoma; LAML: acute myeloid leukemia; LIHC: liver hepatocellular carcinoma; LUAD: lung adenocarcinoma; LUSC: lung squamous cell carcinoma; MESO: mesothelioma; OV: ovarian serous cystadenocarcinoma; PAAD: pancreatic adenocarcinoma; PCPG: pheochromocytoma and paraganglioma; PRAD: prostate adenocarcinoma; READ: rectum adenocarcinoma; SARC: sarcoma; SKCM: skin cutaneous melanoma; STAD: stomach adenocarcinoma; TGCT: testicular germ cell tumors; THCA: thyroid carcinoma; THYM: thymoma; UCEC: uterine corpus endometrial carcinoma; UCS: uterine carcinosarcoma; and UVM: uveal melanoma; AML: ; Lung_AC: lung adenocarcinoma; Lung_SC: lung squamous cell carcinoma; Renal_CC: renal clear cell carcinoma; Renal_CH: renal chromophobe cell carcinoma; Renal_PA: renal papillary cell carcinoma; Uterus_CS: uterine carcinosarcoma; Uterus_EC: uterine corpus endometrial carcinoma

High UBE2W expression correlated with poor prognosis in breast Cancer

PrognoScan database was used to analyze the correlation of UBE2W expression with prognosis. Six out of eleven cancers showed a potential relevance between the UBE2W expression with prognosis (Supplementary Table 2). A higher expression level of UBE2W indicated a bad overall survival (OS) (Fig. 2d) and relapse free survival (RFS) (Fig. 2e) in breast cancer. The same trend was also seen in melanoma, meningioma, colorectal cancer (Supplementary Figure 2A-D). While, higher UBE2W expression was correlated to the better survival prognosis among the other three solid cancers (astrocytoma, glioma, lung cancer) (Supplementary Figure 2E-I).

Fig. 2
figure2

Survival analysis and predictive value analysis of UBE2W in public databases. a-c Overall survival (OS), Relapse free survival (RFS), and Distant metastasis-free survival curves of BRCA by different expression levels of UBE2W in the K-M plotter database. d-e OS and RFS curves of BRCA by different expression levels of UBE2W in the PrognoScan database. f Gene chip data of UBE2W expression in normal, tumor, metastatic tissues. g RNA-seq data of UBE2W expression in normal, tumor, metastatic tissues. h UBE2W expression between responders and non-responders in endocrine therapy of BRCA. i The ROC plot of UBE2W expression in predicting the effect of endocrine therapy in BRCA. Abbreviations: OS: overall survival; RFS: relapse-free survival; DMFS: distant metastasis free survival

Next, we further analyze the prognostic value of UBE2W in the KM plotter database. Notably, higher UBE2W expression was still correlated to the bad prognosis in breast cancer not only in OS (Fig. 2a) and RFS (Fig. 2b) but also in distant-metastasis free survival (DMFS) (Fig. 2c, HR = 1.24, p = 0.029) (p-value < 0.05). Besides, it is also significantly related to the shorter survival in ovarian cancer (Supplementary Figure 2 J), LIHC (Supplementary Figure 2 K), rectum adenocarcinoma (Supplementary Figure 2 L), esophageal adenocarcinoma, UCEC (data has not shown). Consistent with the results of PrognoScan analysis, a higher expression level of UBE2W was correlated with better prognostic in lung cancer (Supplementary Figure 2D), gastric cancer, KIRC, and esophageal squamous cell carcinoma. In general, UBE2W has an impotent prognostic value in various tumors, especially in breast cancer. And UBE2W may play a crucial role in the progress of breast cancer.

High expression UBE2W deteriorates the outcomes of BRCA by promoting metastasis and endocrine resistance

Besides, we further explored the mechanism of why the UBE2W expression could affect the prognosis in breast cancer. First, the difference of UBE2W expression levels in the normal, tumor, metastatic tissues was analyzed via the TNMplot database. Both gene-chip data and RNA-seq data reveal a higher expression level of UBE2W in metastatic samples than in tumors (Fig. 2f, g). Survival analysis further verified the high expression UBE2W was correlated with poor distant metastasis-free survival (Fig. 2c). Second, the correlation between the expression level of UBE2W and clinical variables was analyzed by KM plotter analysis (Table 1). Specifically, compared with low UBE2W mRNA expression level, high UBE2W mRNA expression level was related to poor prognosis only in luminal A subtype disease (OS HR = 1.43, p = 0.046, RFS HR = 1.34, p < 0.001), which was defined as human epidermal growth factor receptor (HER2) negative, progesterone receptor (PR) positive, estrogen receptor (ER) positive, Ki-67 < 14%. Then, we explored the predictive value of UBE2W expression level in response to endocrine therapy of breast cancer. By the ROC plotter analysis, the expression level of UBE2W was higher in non-responder patients than in responder patients (Fig. 2h). The expression level of UBE2W could predict the response to endocrine therapy with an AUC value of 0.663 and p < 0.05 (Fig. 2i). The above results indicate that high UBE2W expression could promote tumor metastasis and cause resistance to endocrine therapy in breast cancer.

Table 1 Correlation of UBE2W mRNA expression and clinicopathological factors in Breast cancer by Kaplan-Meier plotter database

Association between UBE2W expression with DNA repair and DNA methyltransferase in breast cancer

DNA mismatch repair (MMR) plays a critical role in DNA repair pathways [31]. In similar, BRCA genes contribute to DNA repair and regulate the response of the organism to DNA damage [32]. DNA methyltransferase could affect the expression of tumor-suppressing genes by the means of epigenetic modification [33]. All those genes were certified to relevant with the tamoxifen resistance and the prognosis of breast cancer [34,35,36]. Therefore, a correlation analysis of UBE2W expression with DNA repair genes and DNA methyltransferase genes in BRCA was pursued using the TIMER and verified via the GEPIA database. Results showed that UBE2W expression exhibited a positive relationship to those genes (DNMT1, DNMT2, DNMT3A, DNMT3B, MLH1, MSH2, MSH6, PMS2, EPCAM, BRCA1, BRCA2) in breast cancers (Fig. 3, Table 2). These results indicate that UBE2W may mediate therapy resistance by regulating DNA damage or methylation in breast cancer.

Fig. 3
figure3

Correlation of UBE2W expression with DNA repair and DNA methyltransferase genes in BRCA

Table 2 Correlation of UBE2W expression and various genes in GEPIA database

Correlation of UBE2W expression with tumor immune microenvironment in breast cancer

Antiestrogens could induce the immunosuppression in the tumor microenvironment (TME) of BRCA which caused the drug resistance of endocrine therapy [37]. Then, we explored the correlation of the UBE2W expression with the tumor-infiltrating immune cells in BRCA via TIMER analysis (Fig. 4). We observed that UBE2W expression levels were significantly positively associated with tumor purity (r = 0.176, P = 2.23e-08) in breast cancer. In addition, UBE2W expression was significantly positive correlated with the level of CD8+ T cells (r = 0.205, P = 6.46e-11), neutrophil (r = 0.222, P = 1.56e-12) and macrophage cells (r = 0.342, P = 1.27e-28). However, there were significantly negative correlation with level of CD4+ T cells (r = − 0.198, P = 2.95e-10) and Myeloid dendritic cells (r = − 0.106, P = 8.41e-04). And no significant correlation between the expression level of UBE2W with the level of B cells (r = − 0.064, P = 4.35e-02). These results strongly exhibited that UBE2W could recruit the immune cells in the tumor microenvironment in BRCA.

Fig. 4
figure4

Correlation of UBE2W expression with immune infiltration level in BRCA

Furthermore, we analyzed the relationships between UBE2W expression with various markers of immune cells in BRCA via a public database. In the TIMER database, it was conspicuous that UBE2W expression was significantly correlated with the gene markers of monocytes/macrophages (monocytes, and M2 macrophages) than gene markers of other immune cells (Supplementary Table 3). Elevated UBE2W expression was positively related to the genes of CD206, CD163, IL10 which were the specific markers of M2 macrophages. However, in the GEPIA database, UBE2W expression was a negative correlation with the gene markers of CD8+ T cells and general T cells in BRCA tumors (Supplementary Table 4). Significant relationships were also detected in the gene markers of T cell exhaustion, monocyte with the expression level of UBE2W. In brief, those results hinted that UBE2W may play a promoting role in regulating the polarization of macrophages and decreasing the infiltration of CD8+ T cells in the TME.

Genomic alterations and the co-expression genes of UBE2W in the BRCA

A cohort of 1904 breast cancer patients with complete cancer genomics data was selected from 1 of 18 studies in the cBioPortal database [38,39,40]. The OncoPrint shows that the UBE2W gene was altered in 20% BRCA patients and the mRNA expression level of UBE2W for all patients (Fig. 5a). We further accessed that UBE2W mRNA expression was elevated in amplified cases (Fig. 5b) and that cases with UBE2W alteration have worse overall survival (Fig. 5c) and relapse-free survival (Fig. 5d) than cases without UBE2W alteration (OS: p = 4.790e-3, log-rank test; RFS: p = 3.790e-4, log-rank test).

Fig. 5
figure5

Genomic alterations of UBE2W of BRCA in cBioPortal. a Oncoprint of UBE2W alteration in BRCA. b Mutation details for UBE2W. c, d Survival analysis. Abbreviations: OS: overall survival

To find the interaction E3 ubiquitin ligase of UBE2W, we further accessed the co-expression genes of UBE2W in BRCA and performed a network analysis in the public database. Co-expression genes were obtained from the cBioPortal database (data has not shown). Then a network analysis was constructed via STRING analysis (Supplementary Figure 3B) which implied the RBX1 may be the E3 ubiquitin ligase of UBE2W. To confirm this result, we analyzed the expression level and prognosis value of the RBX1 gene in the BRCA (Supplementary Figure 4). Taken these results together, it revealed that UBE2W could affect the prognosis of BRCA through the interaction of E3 ubiquitin-protein ligase RBX1.

Discussion

UBE2W is a protein-coding gene localized on chromosome 8q21.11 and was first reported in 2006 [41]. At present, researches on UBE2W are mainly about neurodegenerative diseases, Fanconi anemia, Huntington’s disease, and male reproductive systems [12, 42, 43]. UBE2W downregulation could promote spermatogenic cell apoptosis through activating the P53/Bcl-2/caspase 6/caspase 9 signal pathways [44]. The loss of UBE2W could sensitized cells to the DNA damage agents with the absence of Rnf4 [12]. However, the roles of UBE2W in human pan-cancer has not been identified. Whether UBE2W could be used as a biomarker is still unknown. An integrated analysis that we performed reveals that UBE2W was a potential oncogenic gene with prognostic value in BRCA, which is the first study of the UBE2W gene with cancer. According to the results from the Human Protein Atlas (HPA) database (https://www.proteinatlas.org/), UBE2W exhibits a significantly higher expression profile in PTC tissues than in normal tissues (NT) (Supplementary Figure 3A). Notably, we found higher expression of UBE2W in metastatic tissues than primary tumors were observed in BRCA. Meanwhile, the high expression level of UBE2W is significantly correlated with the non-response of endocrine therapy. We further detected that UBE2W expression was correlated with DNA repair and methylation. Immune cell infiltration levels and various immune-related markers were also employed to analyze the relationship with UBE2W in specific cancer. Also, an exploration of the genomic alteration of UBE2W was performed. Therefore, our study provides clues to shed light on the potential effects of UBE2W in cancer research.

In the present research, we verified different expression levels of UBE2W in normal and tumor from multi-omics data integration and analysis in specific cancer. Firstly, the difference mRNA expression level of UBE2W between tumor and normal was performed in various public datasets. Oncomine datasets shown UBE2W expression was upregulated in breast cancer, colorectal cancer, head and neck cancer, lung cancer. In the TIMER dataset, higher UBE2W expression was detected in BRCA, CHOL, COAD, ESCA, GBM, HNSC, LUAD, LUSC, STAD. While in the TNM plot database which contained includes 56,938 unique multilevel quality control samples, results of UBE2W expression were similar to the TIMER analysis. Three public analysis platforms confirmed the higher expression level of UBE2W in BRCA than the normal tissues. Compare the mRNA expression level of UBE2W in the metastatic tumors to the primary tumors, a rising trend was detected. Secondly, we found UBE2W exhibited a significantly higher expression profile in BRCA tumor tissues than in normal tissues. Finally, the genomic analysis revealed that higher UBE2W mRNA expression was related to the copy number amplification in BRCA. In summary, it confirms that UBE2W expression was significantly higher in tumors rather than normal tissues in BRCA. A high expression level of UBE2W may promote tumor metastasis in BRCA. Then, a consistent result regarding the prognosis of UBE2W was found in BRCA via PrognoScan and KM plotter analysis. Meanwhile, UBE2W expression could be a predictive biomarker in the response of endocrine therapy in BRCA which was confirmed by the Cox regression analysis in the KM plotter database. Moreover, UBE2W genomic alteration was also correlated with the poor prognosis survival in BRCA. In brief, those results prompt that UBE2W is a prognostic biomarker for breast cancer. Combined co-expression analysis and survival analysis, we found RBX1 was a gene co-expression with UBE2W and high RBX1 expression was correlated with the poor survival in BRCA. Previous studies pointed out that RBX1 expression could regulate the tumor suppressor degradation in liver cancer and be a prognostic factor in lung cancer [45, 46]. Therefore, we guess that RBX1 may be a new E3-ubiquitin ligase of UBE2W in BRCA.

Dysregulation of DNA damage repair (DDR) genes was reported to be a factor that promoting endocrine therapy resistance [34, 47]. Both BRCA1 and BRCA2 are well-known genes concerning homologous recombination in DDR [48]. A novel resistance mechanism with re-express BRCA1 causes the acquisition of therapy resistance [36]. And BRCA2 mutation was associated with disease progression in breast cancer [49]. While MMRs consist of multiple genes (MLH1, MSH2, MSH6, PML2, EPCAM) were demonstrated to induce endocrine therapy resistance via suppressing the CDK4 activity [34]. And the function loss of those genes brings about a mutator phenotype which causes an increasing tumor risk [50]. In this study, we found that UBE2W expression was closely related to the expression level of BRCA1/2 and 5 MMR genes in breast cancer. Previous studies have shown that DNMTs (DNMT1, DNMT3A, DNMT3B) were overexpressed in tamoxifen-resistance tumors and DNMT3B remained as an independent prognostic survival factor for BRCA [35]. In this study, UBE2W expression was positively related to the expression of DNMTs in breast cancer. Furthermore, we analyzed the correlation of UBE2W with various genes (EGFR, ERBB2, ESR1, FOXM1, GATA3, GREB1, IGF1R, NCOR1, PIK3CA, PTEN) which were reported to be a promoter of endocrine therapy resistance [51, 52]. Results showed that UBE2W expression was positively correlated to all those genes in breast cancer, especially in the luminal A subtype (Supplementary Figure 3C). All those results demonstrated that a high expression level of UBE2W may induce poor prognosis survival through promoting the endocrine therapy resistance in breast cancer, especially in the subtype of luminal A.

It was reported that antiestrogen treatment could induce immunosuppression in TME which may contribute to the endocrine therapy in BRCA [37]. In this study, we found that UBE2W expression was associated with tumor infiltration immune cells in BRCA. Interestingly, UBE2W expression was positive correlated with the expression of CD206, CD163, and IL10 in the TIMER database. While, in the GEPIA database, UBE2W was negative correlated with the immune-specific genes of CD8+ T cells. These results exhibited that UBE2W expression could promote immunosuppression in TME.

Nevertheless, our research has some limitations. Considering that our study was based on public databases, systematic bias cannot be avoided. The mechanism of UBE2W in promoting cancer progression and therapy resistance needs a more rigorous exploration. Although our results discover the correlation of UBE2W expression with TME, further efforts are required to verify whether UBE2W expression could affect the prognosis of BRCA across the tumor infiltration immune cells.

Conclusions

In conclusion, UBE2W was an unfavorable prognostic factor and could be a predictive biomarker of endocrine therapy in BRCA. UBE2W might promote tumor metastasis and the resistance of endocrine therapy though influencing the process of DNA repair and DNA methyltransferase, inducing the tumor immunosuppression. Besides, RBX1 may be a new ubiquitinated ligase of UBE2W in BRCA.

Availability of data and materials

All data are available via the corresponding author.

References

  1. 1.

    Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424. https://doi.org/10.3322/caac.21492.

    Article  Google Scholar 

  2. 2.

    Hartman M, Loy EY, Ku CS, Chia KS. Molecular epidemiology and its current clinical use in cancer. Lancet Oncol. 2010;11(4):383–90. https://doi.org/10.1016/S1470-2045(10)70005-X.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Hoeller D, Hecker CM, Dikic I. Ubiquitin and ubiquitin-like proteins in cancer pathogenesis. Nat Rev Cancer. 2006;6(10):776–88. https://doi.org/10.1038/nrc1994.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Kleiger G, Mayor T. Perilous journey: a tour of the ubiquitin–proteasome system. Trends Cell Biol. 2014;24(6):352–9. https://doi.org/10.1016/j.tcb.2013.12.003.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  5. 5.

    Ebstein F, Keller M, Paschen A, Walden P, Seeger M, Burger E, et al. Exposure to Melan-a/MART-126-35 tumor epitope specific CD8(+)T cells reveals immune escape by affecting the ubiquitin-proteasome system (UPS). Sci Rep. 2016;6(1):25208. https://doi.org/10.1038/srep25208.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Bader M, Steller H. Regulation of cell death by the ubiquitin–proteasome system. Curr Opin Cell Biol. 2009;21(6):878–84. https://doi.org/10.1016/j.ceb.2009.09.005.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Voutsadakis IA. The ubiquitin-proteasome system in colorectal cancer. Biochim Biophys Acta (BBA) - Mol Basis Dis. 2008;1782(12):800–8. https://doi.org/10.1016/j.bbadis.2008.06.007.

    CAS  Article  Google Scholar 

  8. 8.

    Dawson SP. Hepatocellular carcinoma and the ubiquitin–proteasome system. Biochim Biophys Acta (BBA) - Mol Basis Dis. 2008;1782(12):775–84. https://doi.org/10.1016/j.bbadis.2008.08.003.

    CAS  Article  Google Scholar 

  9. 9.

    Tatham MH, Plechanovova A, Jaffray EG, Salmen H, Hay RT. Ube2W conjugates ubiquitin to alpha-amino groups of protein N-termini. Biochem J. 2013;453(1):137–45. https://doi.org/10.1042/BJ20130244.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Alpi A, Pace P, Babu M, Patel K. Mechanistic insight into site-restricted monoubiquitination of FANCD2 by Ube2t, FANCL, and FANCI. Mol Cell. 2008;32(6):767–77. https://doi.org/10.1016/j.molcel.2008.12.003.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Zhang Y, Zhou X, Zhao L, Li C, Zhu H, Xu L, et al. UBE2W interacts with FANCL and regulates the monoubiquitination of fanconi anemia protein FANCD2. Molecules and Cells. 2010;31(2):113–22. https://doi.org/10.1007/s10059-011-0015-9.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Maure J-F, Moser SC, Jaffray EG, Alpi AF, Hay RT. Loss of ubiquitin E2 Ube2w rescues hypersensitivity of Rnf4 mutant cells to DNA damage. Sci Rep. 2016;6(1):26178.

  13. 13.

    Christensen DE, Brzovic PS, Klevit RE. E2–BRCA1 RING interactions dictate synthesis of mono- or specific polyubiquitin chain linkages. Nat Struct Mol Biol. 2007;14(10):941–8. https://doi.org/10.1038/nsmb1295.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Fujita Y, Tinoco R, Li Y, Senft D, Ronai ZA. Ubiquitin ligases in Cancer immunotherapy - balancing antitumor and autoimmunity. Trends Mol Med. 2019;25(5):428–43. https://doi.org/10.1016/j.molmed.2019.02.002.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Song M, Wang C, Wang H, Zhang T, Li J, Benezra R, et al. Targeting ubiquitin protein ligase E3 component N-recognin 5 in cancer cells induces a CD8+ T cell mediated immune response. Oncoimmunology. 2020;9(1):1746148. https://doi.org/10.1080/2162402X.2020.1746148.

    Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Zhou Y, Chen R, Luo X, Zhang WD, Qin JJ. The E2 ubiquitin-conjugating enzyme UbcH5c: an emerging target in cancer and immune disorders. Drug Discov Today. 2020;25(11):1988–97. https://doi.org/10.1016/j.drudis.2020.09.015.

    CAS  Article  Google Scholar 

  17. 17.

    Rhodes DR, Kalyana-Sundaram S, Mahavisno V, Varambally R, Yu J, Briggs BB, et al. Oncomine 3.0: genes, pathways, and networks in a collection of 18,000 cancer gene expression profiles. Neoplasia. 2007;9(2):166–80. https://doi.org/10.1593/neo.07112.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Li T, Fu J, Zeng Z, Cohen D, Li J, Chen Q, et al. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 2020;48(W1):W509–14. https://doi.org/10.1093/nar/gkaa407.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017;77(21):e108–10. https://doi.org/10.1158/0008-5472.CAN-17-0307.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Li B, Severson E, Pignon JC, Zhao H, Li T, Novak J, et al. Comprehensive analyses of tumor immunity: implications for cancer immunotherapy. Genome Biol. 2016;17(1):174. https://doi.org/10.1186/s13059-016-1028-7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  21. 21.

    Bartha Á, Győrffy B. TNMplot.com: a web tool for the comparison of gene expression in normal, tumor and metastatic tissues. Int J Mol Sci. 2021;22(5):2622.

  22. 22.

    Uhlén M, Fagerberg L, Hallström B, Lindskog C, Oksvold P, Mardinoglu A, et al. Proteomics. Tissue-based map of the human proteome. Science (New York, NY). 2015;347(6220):1260419.

    Article  Google Scholar 

  23. 23.

    Gyorffy B, Lanczky A, Eklund AC, Denkert C, Budczies J, Li Q, et al. An online survival analysis tool to rapidly assess the effect of 22,277 genes on breast cancer prognosis using microarray data of 1,809 patients. Breast Cancer Res Treat. 2010;123(3):725–31. https://doi.org/10.1007/s10549-009-0674-9.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Mizuno H, Kitada K, Nakai K, Sarai A. PrognoScan: a new database for meta-analysis of the prognostic value of genes. BMC Med Genet. 2009;2(1):18. https://doi.org/10.1186/1755-8794-2-18.

    CAS  Article  Google Scholar 

  25. 25.

    Fekete J, Gyorffy B. ROCplot.org: validating predictive biomarkers of chemotherapy/hormonal therapy/anti-HER2 therapy using transcriptomic data of 3,104 breast cancer patients. Int J Cancer. 2019;145(11):3140–51. https://doi.org/10.1002/ijc.32369.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Zhang X, Lan Y, Xu J, Quan F, Zhao E, Deng C, et al. CellMarker: a manually curated resource of cell markers in human and mouse. Nucleic Acids Res. 2019;47(D1):D721–8. https://doi.org/10.1093/nar/gky900.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–W102. https://doi.org/10.1093/nar/gkx247.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6(269):pl1. https://doi.org/10.1126/scisignal.2004088.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, et al. The cBio Cancer genomics portal: an open platform for exploring multidimensional Cancer genomics data: figure 1. Cancer Discovery. 2012;2(5):401–4. https://doi.org/10.1158/2159-8290.CD-12-0095.

    Article  Google Scholar 

  30. 30.

    Szklarczyk D, Gable AL, Nastou KC, et al. The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 2021;49:D605–D612.

  31. 31.

    Schofield MJ, Hsieh P. DNA mismatch repair: molecular mechanisms and biological function. Annu Rev Microbiol. 2003;57(1):579–608. https://doi.org/10.1146/annurev.micro.57.030502.090847.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Yoshida K, Miki Y. Role of BRCA1 and BRCA2 as regulators of DNA repair, transcription, and cell cycle in response to DNA damage. Cancer Sci. 2004;95(11):866–71. https://doi.org/10.1111/j.1349-7006.2004.tb02195.x.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Pechalrieu D, Etievant C, Arimondo PB. DNA methyltransferase inhibitors in cancer: from pharmacology to translational studies. Biochem Pharmacol. 2017;129:1–13. https://doi.org/10.1016/j.bcp.2016.12.004.

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    Haricharan S, Punturi N, Singh P, Holloway KR, Anurag M, Schmelz J, et al. Loss of MutL disrupts CHK2-dependent cell-cycle control through CDK4/6 to promote intrinsic endocrine therapy resistance in primary breast Cancer. Cancer Discov. 2017;7(10):1168–83. https://doi.org/10.1158/2159-8290.CD-16-1179.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Jahangiri R, Jamialahmadi K, Gharib M, Emami Razavi A, Mosaffa F. Expression and clinicopathological significance of DNA methyltransferase 1, 3A and 3B in tamoxifen-treated breast cancer patients. Gene. 2019;685:24–31. https://doi.org/10.1016/j.gene.2018.10.060.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    Ter Brugge P, Kristel P, van der Burg E, Boon U, de Maaker M, Lips E, et al. Mechanisms of Therapy Resistance in Patient-Derived Xenograft Models of BRCA1-Deficient Breast Cancer. J Natl Cancer Inst. 2016;108(11).

  37. 37.

    Joffroy CM, Buck MB, Stope MB, Popp SL, Pfizenmaier K, Knabbe C. Antiestrogens induce transforming growth factor beta-mediated immunosuppression in breast cancer. Cancer Res. 2010;70(4):1314–22. https://doi.org/10.1158/0008-5472.CAN-09-3292.

    CAS  Article  PubMed  Google Scholar 

  38. 38.

    Pereira B, Chin SF, Rueda OM, Vollan HK, Provenzano E, Bardwell HA, et al. The somatic mutation profiles of 2,433 breast cancers refines their genomic and transcriptomic landscapes. Nat Commun. 2016;7(1):11479. https://doi.org/10.1038/ncomms11479.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Curtis C, Shah SP, Chin SF, Turashvili G, Rueda OM, Dunning MJ, et al. The genomic and transcriptomic architecture of 2,000 breast tumours reveals novel subgroups. Nature. 2012;486(7403):346–52. https://doi.org/10.1038/nature10983.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Rueda OM, Sammut SJ, Seoane JA, Chin SF, Caswell-Jin JL, Callari M, et al. Dynamics of breast-cancer relapse reveal late-recurring ER-positive genomic subgroups. Nature. 2019;567(7748):399–404. https://doi.org/10.1038/s41586-019-1007-8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Yin G, Ji C, Wu T, Shen Z, Xu X, Xie Y, et al. Cloning, characterization and subcellular localization of a gene encoding a human ubiquitin-conjugating enzyme (E2) homologous to the Arabidopsis thaliana UBC-16 gene product. Front Biosci. 2006;11(1):1500–7. https://doi.org/10.2741/1899.

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Scaglione KM, Zavodszky E, Todi SV, Patury S, Xu P, Rodriguez-Lebron E, et al. Ube2w and ataxin-3 coordinately regulate the ubiquitin ligase CHIP. Mol Cell. 2011;43(4):599–612. https://doi.org/10.1016/j.molcel.2011.05.036.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  43. 43.

    Wang B, Zeng L, Merillat SA, Fischer S, Ochaba J, Thompson LM, et al. The ubiquitin conjugating enzyme Ube2W regulates solubility of the Huntington's disease protein, huntingtin. Neurobiol Dis. 2018;109(Pt A):127–36.

    CAS  Article  Google Scholar 

  44. 44.

    Lei B, Xie L, Zhang S, Lv D, Shu F, Deng Y. UBE2W down-regulation promotes cell apoptosis and correlates with hypospermatogenesis. Andrologia. 2020;52(1):e13474. https://doi.org/10.1111/and.13474.

    Article  PubMed  Google Scholar 

  45. 45.

    Xu J, Li L, Yu G, Ying W, Gao Q, Zhang W, et al. The neddylation-cullin 2-RBX1 E3 ligase axis targets tumor suppressor RhoB for degradation in liver cancer. Mol Cell Proteomics. 2015;14(3):499–509. https://doi.org/10.1074/mcp.M114.045211.

    CAS  Article  PubMed  Google Scholar 

  46. 46.

    Xing R, Chen KB, Xuan Y, Feng C, Xue M, Zeng YC. RBX1 expression is an unfavorable prognostic factor in patients with non-small cell lung cancer. Surg Oncol. 2016;25(3):147–51. https://doi.org/10.1016/j.suronc.2016.05.006.

    Article  PubMed  Google Scholar 

  47. 47.

    Haricharan S, Bainbridge MN, Scheet P, Brown PH. Somatic mutation load of estrogen receptor-positive breast tumors predicts overall survival: an analysis of genome sequence data. Breast Cancer Res Treat. 2014;146(1):211–20. https://doi.org/10.1007/s10549-014-2991-x.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  48. 48.

    Rehman FL, Lord CJ, Ashworth A. Synthetic lethal approaches to breast cancer therapy. Nat Rev Clin Oncol. 2010;7(12):718–24. https://doi.org/10.1038/nrclinonc.2010.172.

    CAS  Article  PubMed  Google Scholar 

  49. 49.

    Gornstein EL, Sandefur S, Chung JH, Gay LM, Holmes O, Erlich RL, et al. BRCA2 reversion mutation associated with acquired resistance to Olaparib in estrogen receptor-positive breast Cancer detected by genomic profiling of tissue and liquid biopsy. Clin Breast Cancer. 2018;18(2):184–8. https://doi.org/10.1016/j.clbc.2017.12.010.

    CAS  Article  PubMed  Google Scholar 

  50. 50.

    Jiricny J. The multifaceted mismatch-repair system. Nat Rev Mol Cell Biol. 2006;7(5):335–46. https://doi.org/10.1038/nrm1907.

    CAS  Article  PubMed  Google Scholar 

  51. 51.

    De Marchi T, Foekens JA, Umar A, Martens JW. Endocrine therapy resistance in estrogen receptor (ER)-positive breast cancer. Drug Discov Today. 2016;21(7):1181–8. https://doi.org/10.1016/j.drudis.2016.05.012.

    CAS  Article  PubMed  Google Scholar 

  52. 52.

    Tryfonidis K, Zardavas D, Katzenellenbogen BS, Piccart M. Endocrine treatment in breast cancer: cure, resistance and beyond. Cancer Treat Rev. 2016;50:68–81. https://doi.org/10.1016/j.ctrv.2016.08.008.

    Article  PubMed  Google Scholar 

Download references

Acknowledgments

We thank all the clinicians and patients who participated in the study.

Funding

This work was supported by the National Natural Science Foundation of China (81672987).

Author information

Affiliations

Authors

Contributions

Y.H.G. had the original idea for the study. Y.H.G. had access to all data in the study and was responsible for the integrity of the data and the accuracy of the data analyses. Y.Y., W.W.X., and W.H.X. participated in analyzing the data and writing of the manuscript. R.Z.L prepared figures and tables. All authors have read, edited, and approved the final version of the manuscript.

Corresponding author

Correspondence to Yuan-Hong Gao.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors report no conflicts of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Figure 1

. Flow diagram. Supplementary Figure 2. Impact of gene expression of UBE2W on various cancer survival in PrognoScan database. Overall survival (OS) of skin cancer(A), eye cancer (B), brain cancer (C, F, G), colorectal cancer (D), ovarian cancer (J), liver hepatocellular carcinoma (K), rectum adenocarcinoma (L), lung cancer (H). Relapse-free survival (RFS) of lung cancer (I). Disease-free survival (DFS) of colorectal cancer (E). Supplementary Figure 3. (A) UBE2W protein expression profile in the HPA database. (B) Network view of theUBE2W co-expression genes in breast cancer. (C) Correlation of UBE2W expression with various genes in BRCA in TIMER. Significant positive correlations by the Spearman test are marked with red color (p <0.05). No significant correlations are marked with gray color (p>0.05). Supplementary Figure 4. Expression and prognosis value of RBX1. (A) RBX1 expression in different cancers in TIMER (***P < 0.001). (B) RBX1 expression level in BRCA in Oncomine. (C) UBE2W expression in normal, tumor, metastatic tissues in TNMplot. (E, F) Overall survival and relapse-free survival analysis of RBX1in Kaplan-Meier Plotter. (G, H) Overall survival and relapse-free survival analysis of RBX1in PrognoScan. Supplementary Table 1. UBE2W expression in cancers vs normal tissue in Oncomine database (positive). Supplementary Table 2. Positive results associated with UBE2W expression in different cancers from the Prognoscan database. Supplementary Table 3. Correlation analysis between UBE2W and markers of infiltrating immune cells in TIMER. Supplementary Table 4. Correlation analysis between UBE2W and markers of infiltrating immune cells in GEPIA.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Yuan, Y., Xiao, WW., Xie, WH. et al. Prognostic value of ubiquitin E2 UBE2W and its correlation with tumor-infiltrating immune cells in breast cancer. BMC Cancer 21, 479 (2021). https://doi.org/10.1186/s12885-021-08234-4

Download citation

Keywords

  • UBE2W
  • Breast cancer
  • Prognosis
  • Tumor immune environment
  • Endocrine therapy resistance