A pan-cancer analysis of USP5: its diagnostic, prognostic and immunological roles in human cancers

doi:10.21203/rs.3.rs-2604402/v1

Download PDF

Article

A pan-cancer analysis of USP5: its diagnostic, prognostic and immunological roles in human cancers

https://doi.org/10.21203/rs.3.rs-2604402/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 02 Jun, 2023

Read the published version in Scientific Reports →

You are reading this latest preprint version

Background

Posttranslational modifications (PTM) of proteins, such as acetylation, deubiquitination, and phosphorylation play important roles in various kinds of cancer progression. Ubiquitin-specific proteinase 5 (USP5), a unique member of deubiquitinating enzymes (DUBs) which recognizes unanchored polyubiquitin specifically, could regulate the stability of many tumorigenesis-associated proteins to influence cancer initiation and progression. However, the diverse biological significance of USP5 in pan-cancer has not been systematically and comprehensively studied.

Methods

Here, we explored the role of USP5 in pan-cancer using The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) database, and we also acquired and analyzed data via various software and web platforms such as R, GEPIA2.0, HPA, TISIDB, cBioPortal, UALCAN, TIMER 2.0, CancerSEA and BioGRID.

Results

USP5 expression was high in most cancers and differed significantly in different molecular and immune subtypes of cancers. In addition, USP5 had certain diagnostic value in multiple cancers, and high expression of USP5 generally predicted poor prognosis for cancer patients. We also found that the most frequent genetic alterations type of USP5 was mutation, and the DNA methylation level of USP5 decreased in various cancers. Furthermore, USP5 expression correlated with cancer-associated fibroblasts (CAFs), endothelial cells (EC) and genetic markers of immunodulators in cancers. Moreover, the result from single cell sequencing showed that USP5 could regulate several tumor biological behaviors such as apoptosis, DNA damage and metastasis. Gene enrichment analysis indicated “spliceosome” and “RNA splicing” may be the critical mechanism for USP5 to involve in cancer.

Conclusion

Taken together, our study elucidates the biological significance of USP5 in the diagnosis, prognosis and immune in human pan-cancer.

Biological sciences/Cancer/Tumour biomarkers

Biological sciences/Genetics/Cancer genetics

Biological sciences/Genetics/Epigenetics

USP5

pan-cancer

prognosis

immune

bioinformatics

Cancer is a leading cause of human death, with significant negative impacts on social health and the economy all over the world[1, 2]. Although the development of surgery, chemoradiotherapy, targeted therapy, and immunotherapy have improved the therapeutic effect on cancer, the prognosis and survival rate of cancer patients still remains poor for many reasons, including drug resistance and side effects[3]. Thus, identifying new pan-cancer biomarkers and treatment targets for cancer is critical to enhancing human health[4].

Deubiquitination is a common posttranslational modification (PTM) of proteins, which is involved in regulating various physiological functions and pathological processes, such as signal transduction and cancer progression[5, 6]. Ubiquitin-specific proteinase 5 (USP5) which is discovered and purified by Wilkinson et al, also known as Ubiquitin isopeptidase T (ISOT), belongs to the ubiquitin-specific proteinase (USP) family, the largest family of deubiquitinating enzymes (DUBs)[7]. USP5 is located near the CD4 gene on human chromosome 12p13 and encodes five separate domains with the size of 93.3-kDa[8]. It is unique in that it can specifically identify and remove ubiquitin from the proximal end of the unanchored polyubiquitin chains[9]. And many researches prove that USP5 regulates a variety of cellular activities, including the repair of DNA double-strand breaks[10], inflammatory responses[11], and stress responses[12].

As for the role of USP5 in cancer, it has attracted the attention from many researchers recently. The functional link between the aberrant USP5 expression and the development of various cancers such as breast cancer[13], lung cancer[14, 15], colorectal cancer[16] and hepatocellular carcinoma[17] has been established by many studies. Meanwhile, USP5 has been demonstrated to be closely correlated with some key molecules and pathways regulating cancer, which indicates the potential value of USP5 as a novel treatment target for cancer[18]. However, the detailed role of USP5 in pan-cancer remains elusive so far.

Nowadays, the rapid development of biological databases makes the bioinformatics analyses much more reliable and representative with large sample size and advanced algorithms. In this study, we explored USP5 expression profiles, diagnostic value, prognostic value, genetic alteration, protein methylation level, immune infiltration, functional states and functional enrichment in pan-caner by using multiple bioinformatics methods. And this thorough analysis revealed the certain value of USP5 in the diagnosis and prognosis of various cancers, the potential role of USP5 in some unexplored cancers, the underlying molecular mechanisms of USP5 in the pathogenesis of human cancers and the implications of USP5 in anti-tumor immune response.

Gene expression analysis

Human Protein Atlas (HPA) was used to get the USP5 mRNA and protein expression levels in human organs/tissues. The USP5 mRNA data across 33 cancer types and corresponding paracancer and normal samples were obtained from The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) databases. R software (v 3.6.3) was used to perform statistical analysis, and the “ggplot2” (v3.3.3) package was used for visualization. The Wilcoxon rank-sum test was used to detect the difference between groups, and p < 0.05 was considered statistically significant. Expression Profiling Interactive Analysis 2 (GEPIA2.0) was used to analyze the correlation between the expression of USP5 and patients’ pathological stage in all TCGA cancers. HPA was further used to evaluate expression differences of USP5 at the protein level.

Usp5 Expression In Molecular And Immune Subtypes Of Cancers

TISIDB database (http://cis.hku.hk/TISIDB/) which composed of many data types to evaluate the interaction between cancer and immune system was used to analyze the relationship between USP5 expression and molecular or immune subtypes in pan-cancer.

Diagnostic Value Analysis

The receiver operating characteristic (ROC) curve was used to estimate the diagnostic value of USP5 in pan-cancer, via using the data of the mRNA expression of USP5 in cancer and normal tissues in TCGA and GTEx. Package “pROC” (v1.17.0.1) was used to calculate the ROC curves, and the “ggplot2” (v3.3.3) package was used for plotting. The closer the area under the curve (AUC) is to 1, the better the diagnostic accuracy is. AUC in 0.5–0.7 means low accuracy, AUC in 0.7–0.9 means good accuracy, and AUC in 0.9-1 means high accuracy.

Survival Prognosis Analysis

Kaplan–Meier plots conducting by the “survival” (v3.2-10) package was used to explore the relationship between USP5 expression and prognosis (overall survival (OS), disease-specific survival (DSS) and progress free interval (PFI)) of cancers. The “survminer” (v0.4.9) package was used for visualization. And the forest plot was conducted for summarizing and presenting the results of univariate Cox regression.

Genetic Alteration And Dna Methylation Analysis

The cBioPortal (https://www.cbioportal.org/) was searched for the gene alternations of USP5 in TCGA PanCancer Atlas Studies. The genetic alterations and mutation site information were explored with the “Oncoprint”, “Cancer Type Summary” and “Mutations” modules. And the effect of the gene alterations of USP5 on clinical prognosis, including progression-free survival (PFS), DSS, disease-free survival (DFS), and OS, for all TCGA cancer types was analyzed in the “Comparison” module. Methylation level of USP5 in cancers and corresponding normal tissues was investigated in the UALCAN database (http://ualcan.path.uab.edu/analysis.html).

Immunogenomic Analyses

Various algorithms, such as EPIC, MCPCOUNTER, QUANTISEQ, TIDE, TIMER and XCELL were applied to analyze the relationship between USP5 expression and immune infiltration levels across all TCGA cancers, using TIMER2.0 tool. And we also investigated the correlations between USP5 expression and immunodulators, MHC molecules, chemokines, and chemokine receptors in pan-cancer from the TISIDB database.

Single Cell Sequencing

Using CancerSEA, we explored the correlation between USP5 expression and different functional status of cancer cells at the single cell levels. The threshold for correlation between USP5 and cancer functional status was set at correlation ≥ 0.3 and p-value < 0.05. Expression profiles of USP5 at single cells was showed by the T-SNE diagrams.

Gene Enrichment Analysis

BioGRID was used to explore potential protein interactions with USP5. GEPIA2.0 was applied to collect the top 100 USP5-correlated genes from all TCGA cancer and normal tissues. Then pairwise gene-gene Pearson correlation analysis was performed between USP5 and the selected genes. Heatmap was used to represent the expression status of selected genes containing the partial correlation (cor) and p-value. Gene ontology (GO) and Kyoto encyclopedia of genes and genome (KEGG) enrichment analyses for USP5-correlated genes were conducted via the “clusterProfiler” (v3.14.3) and “org.Hs.eg.db” (v3.10.0) package. “ggplot2” (3.3.3) was used to present the results. P-value < 0.05 was considered to be statistically significant.

USP5 expression in human organs/tissues and pan-cancer

Initially, we examined the mRNA and protein expression levels of USP5 in various organs or tissues. The results obtained from consensus dataset created by combing the HPA and GTEx transcriptomics datasets showed that mRNA of USP5 mainly expressed in skeletal muscle, skin, tongue, parathyroid gland, testis, pons, esophagus, cerebral cortex, adrenal gland and thymus (Supplementary Fig. 1A, B). In addition, the protein of USP5 primarily expressed in cerebellum, testis, cerebral cortex, hippocampus, caudate, thyroid gland, parathyroid gland, adrenal gland, bronchus and lung demonstrated by the HPA dataset (Supplementary Fig. 1A, C).

Then, we analyzed the mRNA expression level of USP5 in pan-cancer. Unpaired analysis of USP5 mRNA expression between paracancerous tissues and cancers revealed that USP5 expressed significantly higher in BRCA (Breast invasive carcinoma), CHOL (Cholangiocarcinoma), COAD (Colon adenocarcinoma), ESCA (Esophageal carcinoma), HNSC (Head and Neck squamous cell carcinoma), KIRP (Kidney renal papillary cell carcinoma), LIHC (Liver hepatocellular carcinoma), LUAD (Lung adenocarcinoma), LUSC (Lung squamous cell carcinoma), PCPG (Pheochromocytoma and Paraganglioma), STAD (Stomach adenocarcinoma) and UCEC (Uterine Corpus Endometrial Carcinoma), and significantly lower in GBM (Glioblastoma multiforme), KICH (Kidney Chromophobe) and PRAD (Prostate adenocarcinoma). There was no obvious difference shown in BLCA (Bladder Urothelial Carcinoma), CESC (Cervical squamous cell carcinoma and endocervical adenocarcinoma), KIRC (Kidney renal clear cell carcinoma), PAAD (Pancreatic adenocarcinoma), READ (Rectum adenocarcinoma) and THCA (Thyroid carcinoma). ACC (Adrenocortical carcinoma), DLBC (Lymphoid Neoplasm Diffuse Large B-cell Lymphoma), LAML (Acute Myeloid Leukemia), LGG (Brain Lower Grade Glioma), MESO (Mesothelioma), OV (Ovarian serous cystadenocarcinoma), SARC (Sarcoma), SKCM (Skin Cutaneous Melanoma), THYM (Thymoma), TGCT (Testicular Germ Cell Tumors), UCS (Uterine Carcinosarcoma) and UVM (Uveal Melanoma) were unable to be analyzed due to the lack of sufficient paracancerous samples (Fig. 1A). Moreover, the paired sample analysis showed that compared with paracancerous tissues, USP5 was overexpressed in BLCA, BRCA, CHOL, ESCA, HNSC, KIRC, KIRP, LIHC, LUAD, LUSC and STAD. On the contrary, USP5 was decreased in KICH and PRAD. And no significantly differential expression of USP5 observed in COAD, PAAD, READ, THCA and UCEC (Fig. 1B). Given the lack of paracancerous tissues in some analysis, we further detected the expression differences of USP5 using the combination of TCGA and GTEx. And the results showed that low expression of USP5 was only observed in LAML, and high expression of USP5 was observed in ACC, BRCA, CESC, CHOL, COAD, DLBC, ESCA, GBM, HNSC, KIRP, LGG, LIHC, LUAD, LUSC, OV, PAAD, PCPG, PRAD, READ, SKCM, STAD, TGCT, THCA, THYM, UCEC and UCS. There was no difference shown in BLCA, KICH and KIRC. Analysis of MESO, SARC and UVM was not possible due to lack of sufficient normal samples (Fig. 1C). Furthermore, we used GEPIA2.0 to explore the effect of USP5 mRNA expression on patient’s pathological stages. And we found that the expression of USP5 was significantly correlated with the pathological stages of CESC, KIRC, LIHC, LUAD, OV and PAAD (Fig. 1D).

At last, we further explored the protein expression level of USP5 in pan-cancer using the National Cancer Institute’s CPTAC dataset and the IHC results provided by the HPA dataset. The result of CPTAC analysis indicated that the protein expression of USP5 was up-regulated and correlated with pathological stages in clear cell RCC (Renal cell carcinoma) and OV (Supplementary Fig. 2A, B). And the IHC results of the HPA showed that staining intensity of USP5 was greater in many cancers, mainly including BRCA, LIHC, OV, PRAD, READ and UCEC, which was consistent with the analysis result of the mRNA expression level of USP5 form TCGA + GTEx (Fig. 2A-F). Overall, USP5 was overexpressed in most cancers.

Usp5 Expression In Different Molecular Subtypes And Immune Subtypes Of Pan-cancer

We analyzed the correlation between USP5 expression and molecular or immune subtypes in pan-cancer from the TISIDB database. The results indicated that USP5 was expressed differently in 10 of 17 cancers for molecular subtypes, including BRCA (5 subtypes), COAD (4 subtypes), ESCA (5 subtypes), HNSC (4 subtypes), KIRP (4 subtypes), LGG (6 subtypes), LUSC (4 subtypes), OV (4 subtypes), PCPG (4 subtypes) and UCEC (4 subtypes) (Fig. 3A-J). Meanwhile, for immune subtypes (C1: wound healing, C2: IFN-gamma dominant, C3: inflammatory, C4: lymphocyte depleted, C5: immunologically quiet, C6: TGF-b dominant), we found that USP5 expression was significantly different in 14 of 30 cancers, including BLCA, BRCA, HNSC, KICH, KIRC, KIRP, LIHC, LUAD, LUSC, MESO, OV, PCPG, SKCM and STAD (Fig. 4A-N).

Diagnostic Value Of Usp5 In Pan-cancer

We plotted the receiver operating curve (ROC) to investigate the diagnostic value of USP5 in pan-cancer. And the ROC curve analysis demonstrated that USP5 had certain diagnostic accuracy (AUC > 0.7) in 20 cancer types, including ACC, BRCA, CESC, CHOL, COAD, DLBC, GBM, HNSC, KIRP, LGG, LIHC, LUAD, LUSC, OV, PAAD, READ, STAD, TGCT, THYM and UCS (Fig. 5A-T). Among them, USP5 had great diagnostic performance (AUC > 0.9) in BRCA, CHOL, DLBC, LGG, LUSC, PAAD and THYM.

Prognostic Value Of Usp5 In Pan-cancer

To evaluate the prognostic assessment value of USP5 in pan-cancer, we carried out the Cox proportional hazards model and Kaplan-Meier analysis. And the result showed that the high level of USP5 predicted poor OS of LAML, LIHC, LUAD, MESO, SKCM and UVM (Fig. 6A, B). For DSS, USP5 played a risk role for BLCA, COAD, LUAD, MESO, SKCM and UVM (Fig. 6A, C). Furthermore, patients with high expression of USP5 had shortened PFI in ACC, COAD, MESO and UVM (Fig. 6A, D).

Genetic Alteration Of Usp5 In Pan-cancer

To investigate the genetic mutations of USP5 in pan-cancer, we used cBioPortal online platform based on TCGA data. The highest frequency of USP5 alteration appeared in AML (Acute Myeloid Leukemia), ACC, BLCA, LGG, BRCA, CESC, CHOL and COADREAD. And mutation, amplification and deep deletion were the most common genetic alterations types of USP5 (Fig. 7A). In addition, we found 149 mutation sites with missense mutation as the main alteration type in USP5. For instance, a missense mutation within the UCH domain, P650L/S alteration was detected in one case of UCEC and two cases of SKCM (Fig. 7B). Then the correlation between the putative CNA of USP5 and its gene expression in pan-cancer was shown in Fig. 7C, D. Moreover, compared with the unaltered group, the gene alteration of AGAP10P, CHD4, VWF, NCAPD2, GPR162, LRRC23, PTPN6, ATN1, LAG3 and CD4 was more predominant in group with USP5 alteration (Fig. 7E). Last, we studied the effect of USP5 genetic alteration on the prognosis of patients in pan-cancer, and the result indicated that patients with USP5 alteration had poor PFS in pan-cancer (Fig. 7F), but not OS, DFS and DSS (Supplementary Fig. 3A-C).

Analysis Of The Methylation Level Of Usp5 In Pan-cancer

DNA methylation has been proved to play an essential role in the occurrence and progression of cancers. Using UALCAN database, we compared the methylation level of USP5 between normal and cancer tissues. We found that the promoter methylation level of USP5 decreased in most cancers, including BLCA, BRCA, CHOL, COAD, ESCA, HNSC, LIHC, LUAD, LUSC, PAAD, PRAD, READ, TGCT and UCEC (Fig. 8A-N). An obvious increase in the methylation level of USP5 was showed in KIRP, KIRC and THCA (Supplementary Fig. 4A). And the difference of the USP5 methylation level was not significant in CESC, GBM, PCPG, SARC, STAD and THYM (Supplementary Fig. 4B).

Immunogenomic Analyses Of Usp5 In Pan-cancer

Considering the critical role of immune infiltration and immune regulation in the oncology progress, we first applied CIBERSORT, CIBERSORT-ABS, EPIC, MCPCOUNTER, QUANTISEQ, TIDE, TIMER and XCELL algorithms to explore the correlation between USP5 expression and the infiltration level of different immune and endothelial cells in pan-cancer of TCGA. The result showed that USP5 expression was positively correlated with the infiltration of cancer-associated fibroblasts in CESC, HNSC and HNSC-HPV- (Fig. 9A). In addition, we discovered a positive correlation between USP5 expression and endothelial cell infiltration in COAD, HNSC-HPV+, SKCM-Metastasis and THCA, while negatively correlated with endothelial cell infiltration in BRCA, BrCA-basal and THYM (Fig. 9B).

Also, we observed that USP5 was correlated with most immune inhibitors and immune stimulators except for KIR2DL1, KIR2DL3 and TNFSF18 in pan-cancer (Supplementary Fig. 5A, B). In terms of MHCs, USP5 was positively associated with most MHCs in KIRC, KIRP and UVM, and negatively associated with most MHCs in ESCA, KICH, LUSC and TGCT (Supplementary Fig. 5C). Moreover, we found that USP5 showed certain correlation with majority of chemokines with the exception of CCL1, CCL16, CCL27, CCL24 and CCL25 in pan-cancer (Supplementary Fig. 5D). Meanwhile, a negative correlation between USP5 and most chemokine receptors could be found in the majority of malignant tumors especially in ESCA, KICH, LUSC and TGCT (Supplementary Fig. 5E).

Functional States Analysis Of Usp5 At Single Cell Levels

Using the CancerSEA, we investigated the functional states of USP5 at single cell levels in various cancers. The results indicated that USP5 was positively correlated with angiogenesis, differentiation, hypoxia, inflammation and metastasis, and negatively correlated with apoptosis, cell cycle, DNA damage, DNA repair, invasion, metastasis, quiescence and stemness (Fig. 10A). Then, the association between USP5 and specific cancer types was further examined. And we observed that USP5 showed positive correlation with hypoxia in LUAD; with metastasis in RCC; with differentiation, angiogenesis and inflammation in RB (Retinoblestoma). In contrast, USP5 negative correlated with DNA repair in ALL (Acute Lymphoblastic Leukemia); with stemness in PC (Prostate cancer); with cell cycle in CRC (Colorectal cancer); with DNA repair, cell cycle and DNA damage in RB (Retinoblestoma); with DNA repair, DNA damage, apoptosis, invasion, metastasis and quiescence in UM (Uveal Melanoma) (Fig. 10B-H). Additionally, T-SNE diagrams were used to display USP5 expression profiles at single cell levels from ALL, LUAD, RCC, PC, CRC, RB and UM (Fig. 10I-O).

Functional Enrichment Analysis Of Usp5 In Pan-cancer

Finally, we screened out the USP5 co-expressed genes for a series of pathway enrichment analyses to understand the molecular mechanism of the USP5 gene in carcinogenesis and development. We first collected 179 molecules that interacted with USP5 via the BioGRID web service (Fig. 11A). Then we used GEPIA to acquire the top 100 USP5 co-expressed genes in pan-cancer. Among these, MLF2, COPS7A, PEX5, DDX47, STRAP and MRPL51 displayed strong correlations with USP5 in most cancer types (Fig. 11B, C). Furthermore, GO and KEGG enrichment analyses were used to reveal that USP5 co-expressed genes played a critical role in the regulation of spliceosome, RNA splicing, catalytic activity acting on RNA and histone binding in tumor pathogenesis (Fig. 11D).

Emerging studies indicated that USP5, a unique member of DUBs that can specifically recognize unanchored polyubiquitin, played an essential role in regulating the repair of DNA double-strand breaks[10], inflammatory responses[11], and stress responses[12]. Meanwhile, USP5 could regulate the stability of many tumorigenesis-associated proteins to influence the progression of various cancers such as, hepatocellular carcinoma[17], pancreatic ductal adenocarcinoma[19, 20], and non-small cell lung cancer (NSCLC)[14, 15]. However, the significance of USP5 in pan-cancer has not been explored until now. In the present study, using multiple bioinformatics approach, we first revealed the abnormal expression of USP5 in human cancers and its expression level in different molecular and immune subtypes of cancers, then we explored the diagnostic and prognostic values of USP5 in various cancers. In addition, we analyzed the gene mutations and methylation levels of USP5 in pan-cancer. Furthermore, the correlation between USP5 expression and the infiltration levels of immune cells and regulators was investigated, and the underlying functions of USP5 at single-cell levels was also identified. Finally, we implemented the functional enrichment analysis to recognize the potential mechanisms for USP5 to influence the pathogenesis of cancers.

Many studies revealed that USP5 was overexpressed and closely correlated with occurrence and progression of various cancers[18]. In line with their researches, our findings from TCGA and GTEx also demonstrated that the expression level of USP5 was significantly higher in most cancers compared with their regular counterparts. In addition, we found that there were meaningful correlations between USP5 expression level and the different molecular or immune subtypes of cancers, which suggested us to get a deeper understanding of USP5’s function in cancer by targeting specific molecular or immune subtypes.

Based on our results, USP5 had a certain diagnostic accuracy (AUC > 0.7) in 20 cancer types, especially in predicting BRCA, CHOL, DLBC, LGG, LUSC, PAAD and THYM (AUC > 0.9), indicating the potential clinical application value of USP5 as a reliable diagnostic biomarker. Meanwhile, as a unique member of DUBs, USP5 could regulate the stability of some key tumor-regulating molecular such as, HIF2α in BRCA[13], CyclinD1 in GBM[21], SLUG in LIHC[17], HDAC2 in OV[22], STAT3[23], WT1[20] and FOXM1[24] in PAAD and β-catenin[25], PD-L1[26] and CCND1[15] in NSCLC to promote cancer progression, and our ROC curve analysis showed the AUC of 0.905 in BRCA, 0.852 in GBM, 0.891 in LIHC, 0.820 in OV, 0.980 in PAAD and 0.946 in LUSC. Thus, the combination of USP5 with these tumor-related molecular separately may significantly improve the diagnostic accuracy for the above cancers. Using cox proportional hazards model and Kaplan-Meier analysis, we found that USP5 was negatively correlated with patients’ prognosis generally. Specifically, USP5 expression indicated poor OS in patients with LAML, LIHC, LUAD, MESO, SKCM and UVM. In addition, we further analyzed the association between the expression level of USP5 and DSS or PFI of cancer patients, and we proved that on the whole, USP5 exhibited risk role in MESO and UVM for OS, DSS and PFI, in LUAD and SKCM for OS and DSS and in COAD for DSS and PFI. In addition to the previous reported negative association between USP5 and the progression of BRCA[13], BLCA[27], CRC[28], GBM[21], LIHC[17], melanoma[29], NSCLC[15, 26], OV[22] and PAAD[23], our result first showed that USP5 may emerge as a novel biomarker for predicting the prognosis of ACC, LAML and MESO, especially MESO. These results suggested that USP5 had important diagnostic and prognostic implications in various cancers, and may serve as a therapeutic target for precision oncology.

USP5 gene was located in 12p13.31. The mutation of USP5 had been reported to cause several tissue disorders in drosophila, including severe defects in the eye development[30, 31]. And the active site (C335A) mutation was proved to prevent the deubiquitination activity of USP5[32]. However, the studies on the USP5 gene alteration in human cancers were still rare. Here, we observed that USP5 genetic alterations, including mutation and amplification could be found in various cancer types. And the frequency of CHD4, VWF, NCAPD2, GPR162, LRRC23, PTPN6, ATN1, LAG3 and CD4 alterations was obviously higher in USP5 alteration group. Moreover, the patients with USP5 alteration had shorten PFS in pan-caner. DNA methylation, a common type of epigenetic modifications, generally inhibiting gene expression via altering chromatin structure, DNA stability, and DNA conformation, played an essential role in multiple types of tumorigenesis[33, 34]. In this study, evidence indicated that DNA methylation level of USP5 was down-regulated in the majority of common malignancies, which was consistent with the elevation of USP5 expression. Further studies on the gene alteration of USP5 and the relationship between DNA methylation and USP5 expression in caner are needed.

Tumor immune microenvironment (TIME), an essential part of tumor microenvironment (TME), mainly composed of immune cells, played a critical role in cancer progression[35–37]. Identifying new targets for immunotherapy was important for improving clinical outcomes, and the impact of USP5 on the TIME was rarely explored so far. The infiltrating immune cells were closely correlated with tumor growth, metastasis and invasion[38, 39]. For example, cancer-associated fibroblasts, tumor-activated fibroblasts, could promote tumor development by secreting various cytokines or metabolites and forming barrier by shaping external-cellular matrix to inhibit the function of drugs and immune cells[40, 41]. In addition, the proliferation of tumor endothelial cells had a protective function of tumor cells by preventing the blood lymphocytes from leaking out of the blood vessels and transporting them to the tumor[42]. And various immune cells recruited by progressed tumors can affect tumor growth, invasion and pathological angiogenesis by promoting the secretion of cytokines and chemokines[43, 44]. In this study, we conducted analysis to assess the impact of USP5 on immune infiltration. And the result obtained through a variety of immune deconvolution methods showed that USP5 was significantly associated with the infiltration of immune cells, including CAFs and EC in certain tumors. Meanwhile, certain correlations between USP5 and various immunoregulation-related genes were found in many cancer types. Generally, our study suggested the potential value of USP5 as an effective target for immunotherapy to enhance the health of cancer patients. More preclinical and clinical trials are needed to explore the relationship between USP5 expression and immune checkpoints.

There was no doubt that USP5 played an essential role in tumorigenesis, but the underlying mechanism still remains elusive. Single-cell transcriptome sequencing was the key technique to analyze the potential functions of molecules at single-cell levels[45]. Using CancerSEA, we found that USP5 significantly correlated with many biological behaviors of cancers such as apoptosis, cell cycle, DNA damage, metastasis and invasion in several cancer types at single cell levels. Additionally, via the functional enrichment analyses of USP5 co-expressed genes, we showed that “spliceosome” and “RNA splicing” may be the critical mechanism for USP5 to involve in pan-cancer. Previous studies had demonstrated that USP5 could regulate cancers by mediating EMT, such as, in LIHC by stabilizing SLUG[17], in NSCLC by stabilizing β-catenin[25] and in BRCA by stabilizing HIF2α[13]. In addition, it had been found that silencing of USP5 may increase apoptosis and DNA damage to suppress the progression of PAAD[46]. Moreover, USP5 could promote cell cycle progression by preventing HERC4-mediated polyubiquitination of c-Maf in multiple myeloma[47], and the downregulation of USP5 lead to the cell cycle arrest in UCEC[48]. Aberrant RNA splicing was thought to be critical in tumorigenesis[49]. In GBM, USP5 isoform2 was closely associated with the aberrant expression of polypyrimidine tract-binding protein 1 (PTBP1), an RNA splicing factor in GBM, and the forced expression of USP5 isoform1 inhibited cell growth and migration in two GBM cell lines, implying an essential role of individual USP5 isoforms generated by alternative splicing in gliomagenesis[50]. The potential molecular mechanism of USP5 related to tumorigenesis and whether USP5 could be a target for cancer therapy still need more experimental exploration.

In summary, using comprehensive bioinformatics analysis methods, our study explored the expression levels, potential diagnostic and prognostic value, genetic mutation, protein methylation, immunomodulatory effects and relevant signaling pathways of USP5 in pan-cancer. The results indicate that USP5 is overexpressed and has certain diagnostic value in various cancer types. In addition, USP5 may be a potential prognostic and immune-related biomarker for cancer patients. This study clarifies the role of USP5 in tumorigenesis from multiple perspectives, providing some bases for further research on the specific mechanisms of USP5 in the progression and treatment of cancers.

Acknowledgements

We would like to thank all members for their contributions to this work. In addition, we sincerely thank the public databases, including TCGA and GTEx for providing the comprehensive data, and all bioinformatics tool for data analysis.

Fundings

This work was supported by Science and Technology Project of Zhuzhou (2019-005[D.C.]) Special socialized investment project of the Innovative city construction of Zhuzhou (2021-047[M.H.]).

Ethical Approval

Not applicable.

Author contributions

Conception and Design: MH; Acquisition and analysis of Data: BY; Writing the Manuscript: BY and JG; Supervision and edit: SD and DC. This manuscript was approved by all authors.

Competing interests

The authors declare no conflict of interests.

Data availability statement

The names of original contributions in this study are included in the article/supplementary material. Further inquiries can be directed to the corresponding authors.

Sung, H., et al., Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin, 2021. 71(3): p. 209-249.
Tao, S., et al., Construction and Clinical Translation of Causal Pan-Cancer Gene Score Across Cancer Types. Front Genet, 2021. 12: p. 784775.
Siegel, R.L., K.D. Miller, and A. Jemal, Cancer statistics, 2019. CA Cancer J Clin, 2019. 69(1): p. 7-34.
Shi, X., et al., Comprehensive Analyses of the Expression, Genetic Alteration, Prognosis Significance, and Interaction Networks of m(6)A Regulators Across Human Cancers. Front Genet, 2021. 12: p. 771853.
Lin, J., X. Bao, and X.D. Li, A tri-functional amino acid enables mapping of binding sites for posttranslational-modification-mediated protein-protein interactions. Mol Cell, 2021. 81(12): p. 2669-2681.e9.
Kaushal, K., et al., Deubiquitinating enzymes in cancer stem cells: functions and targeted inhibition for cancer therapy. Drug Discov Today, 2018. 23(12): p. 1974-1982.
Wilkinson, K.D., et al., Metabolism of the polyubiquitin degradation signal: structure, mechanism, and role of isopeptidase T. Biochemistry, 1995. 34(44): p. 14535-46.
Ansari-Lari, M.A., et al., A gene-rich cluster between the CD4 and triosephosphate isomerase genes at human chromosome 12p13. Genome Res, 1996. 6(4): p. 314-26.
Clague, M.J., S. Urbé, and D. Komander, Breaking the chains: deubiquitylating enzyme specificity begets function. Nat Rev Mol Cell Biol, 2019. 20(6): p. 338-352.
Nakajima, S., et al., Ubiquitin-specific protease 5 is required for the efficient repair of DNA double-strand breaks. PLoS One, 2014. 9(1): p. e84899.
Liu, Q., et al., Broad and diverse mechanisms used by deubiquitinase family members in regulating the type I interferon signaling pathway during antiviral responses. Sci Adv, 2018. 4(5): p. eaar2824.
Nostramo, R., et al., The Catalytic Activity of the Ubp3 Deubiquitinating Protease Is Required for Efficient Stress Granule Assembly in Saccharomyces cerevisiae. Mol Cell Biol, 2016. 36(1): p. 173-83.
Huang, W., et al., USP5 promotes breast cancer cell proliferation and metastasis by stabilizing HIF2α. 2022. 237(4): p. 2211-2219.
Zheng, Y., et al., EOAI, a ubiquitin-specific peptidase 5 inhibitor, prevents non-small cell lung cancer progression by inducing DNA damage. BMC Cancer, 2023. 23(1): p. 28.
Zhang, Z., et al., Deubiquitinase USP5 promotes non-small cell lung cancer cell proliferation by stabilizing cyclin D1. 2021. 10(10): p. 3995-4011.
Xu, X., et al., Ubiquitin specific peptidase 5 regulates colorectal cancer cell growth by stabilizing Tu translation elongation factor. 2019. 9(14): p. 4208-4220.
Meng, J., et al., USP5 promotes epithelial-mesenchymal transition by stabilizing SLUG in hepatocellular carcinoma. 2019. 9(2): p. 573-587.
Ning, F., et al., Structure and function of USP5: Insight into physiological and pathophysiological roles. Pharmacol Res, 2020. 157: p. 104557.
Lian, J., et al., Ubiquitin specific peptidase 5 enhances STAT3 signaling and promotes migration and invasion in Pancreatic Cancer. 2020. 11(23): p. 6802-6811.
Li, J., et al., Deubiquitinase inhibitor degrasyn suppresses metastasis by targeting USP5-WT1-E-cadherin signalling pathway in pancreatic ductal adenocarcinoma. 2020. 24(2): p. 1370-1382.
Li, G., et al., USP5 Sustains the Proliferation of Glioblastoma Through Stabilization of CyclinD1. 2021. 12: p. 720307.
Du, Y., et al., Ubiquitin specific peptidase 5 promotes ovarian cancer cell proliferation through deubiquitinating HDAC2. Aging (Albany NY), 2019. 11(21): p. 9778-9793.
Lian, J., et al., Ubiquitin specific peptidase 5 enhances STAT3 signaling and promotes migration and invasion in Pancreatic Cancer. J Cancer, 2020. 11(23): p. 6802-6811.
Li, X.Y., et al., USP5 promotes tumorigenesis and progression of pancreatic cancer by stabilizing FoxM1 protein. Biochem Biophys Res Commun, 2017. 492(1): p. 48-54.
Xue, S., et al., viaUSP5 Promotes Metastasis in Non-Small Cell Lung Cancer by Inducing Epithelial-Mesenchymal Transition Wnt/β-Catenin Pathway. 2020. 11: p. 668.
Pan, J., et al., USP5 facilitates non-small cell lung cancer progression through stabilization of PD-L1. 2021. 12(11): p. 1051.
Xiao, M., et al., MAFG-AS1 promotes tumor progression via regulation of the HuR/PTBP1 axis in bladder urothelial carcinoma. Clin Transl Med, 2020. 10(8): p. e241.
Xu, X., et al., Ubiquitin specific peptidase 5 regulates colorectal cancer cell growth by stabilizing Tu translation elongation factor. Theranostics, 2019. 9(14): p. 4208-4220.
Chen, X.H., et al., Mebendazole elicits potent antimyeloma activity by inhibiting the USP5/c-Maf axis. Acta Pharmacol Sin, 2019. 40(12): p. 1568-1577.
Fan, X., et al., Drosophila USP5 controls the activation of apoptosis and the Jun N-terminal kinase pathway during eye development. PLoS One, 2014. 9(3): p. e92250.
Kovács, L., et al., Role of the deubiquitylating enzyme DmUsp5 in coupling ubiquitin equilibrium to development and apoptosis in Drosophila melanogaster. PLoS One, 2015. 10(3): p. e0120875.
Scott, D., R. Layfield, and N.J. Oldham, Ion mobility-mass spectrometry reveals conformational flexibility in the deubiquitinating enzyme USP5. Proteomics, 2015. 15(16): p. 2835-41.
Mehdi, A. and S.A. Rabbani, Role of Methylation in Pro- and Anti-Cancer Immunity. Cancers (Basel), 2021. 13(3).
Wang, M., V. Ngo, and W. Wang, Deciphering the genetic code of DNA methylation. Brief Bioinform, 2021. 22(5).
Lei, X., et al., Immune cells within the tumor microenvironment: Biological functions and roles in cancer immunotherapy. Cancer Lett, 2020. 470: p. 126-133.
Hinshaw, D.C. and L.A. Shevde, The Tumor Microenvironment Innately Modulates Cancer Progression. Cancer Res, 2019. 79(18): p. 4557-4566.
Talty, R. and K. Olino, Metabolism of Innate Immune Cells in Cancer. Cancers (Basel), 2021. 13(4).
Steven, A. and B. Seliger, The Role of Immune Escape and Immune Cell Infiltration in Breast Cancer. Breast Care (Basel), 2018. 13(1): p. 16-21.
Cui, X., et al., A pan-cancer analysis of the oncogenic role of staphylococcal nuclease domain-containing protein 1 (SND1) in human tumors. Genomics, 2020. 112(6): p. 3958-3967.
Miao, L., et al., Targeting Tumor-Associated Fibroblasts for Therapeutic Delivery in Desmoplastic Tumors. Cancer Res, 2017. 77(3): p. 719-731.
Curtis, M., et al., Fibroblasts Mobilize Tumor Cell Glycogen to Promote Proliferation and Metastasis. Cell Metab, 2019. 29(1): p. 141-155.e9.
Sobierajska, K., et al., Endothelial Cells in the Tumor Microenvironment. Adv Exp Med Biol, 2020. 1234: p. 71-86.
Hamarsheh, S., et al., Immune modulatory effects of oncogenic KRAS in cancer. Nat Commun, 2020. 11(1): p. 5439.
Benner, B., et al., Generation of monocyte-derived tumor-associated macrophages using tumor-conditioned media provides a novel method to study tumor-associated macrophages in vitro. J Immunother Cancer, 2019. 7(1): p. 140.
He, J., et al., Intra-Tumoral Expression of SLC7A11 Is Associated with Immune Microenvironment, Drug Resistance, and Prognosis in Cancers: A Pan-Cancer Analysis. Front Genet, 2021. 12: p. 770857.
Kaistha, B.P., et al., The deubiquitinating enzyme USP5 promotes pancreatic cancer via modulating cell cycle regulators. Oncotarget, 2017. 8(39): p. 66215-66225.
Zhang, Z., et al., The ubiquitin ligase HERC4 mediates c-Maf ubiquitination and delays the growth of multiple myeloma xenografts in nude mice. Blood, 2016. 127(13): p. 1676-86.
Li, Y., J.J.C.m. Zhou, and research, USP5 Promotes Uterine Corpus Endometrial Carcinoma Cell Growth and Migration via mTOR/4EBP1 Activation. 2021. 13: p. 3913-3924.
Wang, E. and I.J.T.i.c. Aifantis, RNA Splicing and Cancer. 2020. 6(8): p. 631-644.
Izaguirre, D.I., et al., PTBP1-dependent regulation of USP5 alternative RNA splicing plays a role in glioblastoma tumorigenesis. Mol Carcinog, 2012. 51(11): p. 895-906.

No competing interests reported.

SupplementaryMaterialsUSP5pancancer.docx

Download PDF

Journal Publication

published 02 Jun, 2023

Read the published version in Scientific Reports →

Editorial decision: Major revision
09 May, 2023
Reviews received at journal
25 Apr, 2023
Reviewers agreed at journal
24 Apr, 2023
Reviewers agreed at journal
17 Apr, 2023
Reviewers invited by journal
17 Apr, 2023
Editor assigned by journal
17 Apr, 2023
Editor invited by journal
21 Feb, 2023
Submission checks completed at journal
21 Feb, 2023
First submitted to journal
19 Feb, 2023

You are reading this latest preprint version

A pan-cancer analysis of USP5: its diagnostic, prognostic and immunological roles in human cancers

Status:

Journal Publication

Version 1

Abstract

Background

Methods

Results

Conclusion

Figures

Introduction

Materials And Methods

Gene expression analysis

Usp5 Expression In Molecular And Immune Subtypes Of Cancers

Diagnostic Value Analysis

Survival Prognosis Analysis

Genetic Alteration And Dna Methylation Analysis

Immunogenomic Analyses

Single Cell Sequencing

Gene Enrichment Analysis

Results

USP5 expression in human organs/tissues and pan-cancer

Usp5 Expression In Different Molecular Subtypes And Immune Subtypes Of Pan-cancer

Diagnostic Value Of Usp5 In Pan-cancer

Prognostic Value Of Usp5 In Pan-cancer

Genetic Alteration Of Usp5 In Pan-cancer

Analysis Of The Methylation Level Of Usp5 In Pan-cancer

Immunogenomic Analyses Of Usp5 In Pan-cancer

Functional States Analysis Of Usp5 At Single Cell Levels

Functional Enrichment Analysis Of Usp5 In Pan-cancer

Discussion

Declarations

References

Additional Declarations

Supplementary Files

Status:

Journal Publication

Version 1