04 March 2021: Database Analysis
Identification of Biomarkers Based on Bioinformatics Analysis: The Expression of Ubiquitin-Conjugating Enzyme E2T (UBE2T) in the Carcinogenesis and Progression of Hepatocellular CarcinomaXuemiao Zhao 1ABEF , Wei Weng 2ACDE , Miao Jin 3CD , Sunjian Li 4AB , Qingwei Chen 5EF , Bingrong Li 2BC , Zhou Zhou 6BC , Chuanqiang Lan 2BE , Yunjun Yang 1AEF*
Med Sci Monit 2021; 27:e929023
BACKGROUND: The purpose of this study was to screen and identify key genes in the occurrence and development of hepatocellular carcinoma (HCC) based on bioinformatics analysis.
MATERIAL AND METHODS: Three Gene Expression Omnibus (GEO) series (GSE) – GSE121248, GSE87630, and GSE84598 – were downloaded from the GEO database. GEO2R was used to screen different genes and a Venn diagram was drawn to screen coexpressed differentially expressed genes (DEGs). Coexpressed DEGs were obtained by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis, a protein-protein interaction network diagram was produced by Cytoscape, and module genes were calculated by the Molecular Complex Detection Cytoscape plug-in. Finally, overall survival, progression-free survival, and relapse-free survival analysis of the key genes selected were performed using the online Kaplan-Meier plotter. For the target genes, the online network UCSC Cancer Genome Browser was used to analyze the gene expression profiles of the grade and vascular invasion of HCC.
RESULTS: A total of 296 coexpressed DEGs were obtained from the 3 GSEs and 12 key genes were obtained from the modular analysis. Survival analysis showed that the upregulated genes UBE2T and FBLN5 were involved in the poor prognosis of HCC. Furthermore, the expression of UBE2T was significantly related to the grade and vascular invasion of HCC.
CONCLUSIONS: The expression of the UBE2T gene was significantly upregulated in HCC tissue compared to in normal liver tissue. UBE2T may be a new marker for the diagnosis and subsequent therapy of HCC.
Keywords: Gene Expression Profiling, Protein Interaction Mapping
Hepatocellular carcinoma (HCC) accounts for 85% to 90% of primary liver cancer and is the third leading cause of cancer death in the world . With clinical characteristics such as insidious onset and rapid progression, most patients with HCC are already in the middle and advanced stages when they are diagnosed, thereby missing the optimum opportunity for surgical resection . Making the situation for patients even worse, resistance to sorafenib and lenvatinib, the first-line drugs for the treatment of advanced HCC, develops within a few months of use . Therefore, there is an urgent need to develop more effective clinical treatments for HCC. The drug resistance to sorafenib and lenvatinib forced researchers to pay more attention to gene therapy. Gene therapy is a treatment method that intervenes in and regulates gene expression at the molecular level to achieve therapeutic goals [4,5]. Brown et al reported that blocking carnitine palmitoyltransferase decreased the apoptosis of intrahepatic CD4+ T cells and inhibited HCC tumor formation . Shi et al reported that the FOXP3 gene can inhibit HCC progression through the TGF-β/Smad2/3 signaling pathway . The above studies revealed encouraging potential new therapies for HCC.
Nevertheless, owing to the complexity of the carcinogenesis and progression of HCC, its genetic development mechanism is still not understood. To explore the gene expression characteristics of HCC and find new diagnostic markers and therapeutic targets, we downloaded gene expression profile data from the Gene Expression Omnibus (GEO) database, and analyzed differentially expressed genes (DEGs) between HCC tissue and normal liver tissue based on biometric analysis. The selection of DEGs was performed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis, a protein-protein interaction (PPI) network was drawn, and key genes were selected out. The survival analysis of key genes was calculated through the Kaplan-Meier online tool, and the most characteristic genes in different conditions were analyzed using the University of California Santa Cruz (UCSC) Cancer Genome Browser. The results showed that the ubiquitin-conjugating enzyme E2T (UBE2T) gene had high clinical value and might contribute to the diagnosis and treatment of HCC.
UBE2T is well known as a conspicuous gene that plays a crucial role in many kinds of malignancies, including renal cell carcinoma, lung cancer, cervical cancer, multiple myeloma, osteosarcoma, and, especially, HCC [8,9]. The present study is one of many to illustrate that UBE2T is an independent biomarker for HCC [10,11].
Material and Methods
IDENTIFICATION OF DEGS:
The online tool GEO2R (
GO AND KEGG PATHWAY ANALYSIS:
The Database for Annotation, Visualization and Integrated Discovery (DAVID) database provides a comprehensive set of functional annotation tools for researchers to better understand the biological meaning underlying the innumerable existing genes. GO enrichment analysis is a common method to discriminate the enrichment degree of GO terms of different genes. GO enrichment analysis is divided into 3 aspects: biological process (BP), cell component (CC), and molecular function (MF). KEGG is a genome deciphering database and an enrichment analysis tool for the research of biological regulatory networks. In this study, the enrichment analysis of DEGs was completed using the online tool DAVID (
PPI NETWORK CONSTRUCTION AND MODULE ANALYSIS:
The PPI network information of DEGs was analyzed using the online tool STRING (
HUB GENES ANALYSIS:
Key genes consisted of seed genes selected from each module. The Biological Networks Gene Oncology tool (BINGO) plug-in in Cytoscape was used to analyze and visualize the BP of key genes. Hierarchical cluster analysis of key genes was done by the UCSC Cancer Genome Browser (
IDENTIFICATION OF DEGS IN HCC: The differential genes in each data series were confirmed after screening, with 953 DEGs in GSE121248, 1163 in GSE87630, and 1957 in GSE84598 (Figure 1A–1C). The overlap of the 3 data sets is shown in the Venn diagram (Figure 1D), in which 296 DEGs were coexpressed, consisting of 54 upregulated genes and 242 downregulated genes.
KEGG AND GO ENRICHMENT ANALYSES OF DEGS: The results of GO enrichment showed the downregulated genes of DEGs were mainly enriched in the epoxygenase P450 pathway, oxidation-reduction process complement activation, xenobiotic metabolic process, cellular response to tumor necrosis factor, and proteolysis on BP; extracellular region, organelle membrane, extracellular exosome, extracellular space, and blood microparticle on CC; heme binding, oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen, monooxygenase activity, iron ion binding, oxygen binding, serine-type endopeptidase activity, and catalytic activity on MF; and metabolic pathways, complement and coagulation cascades, chemical carcinogenesis, carbon metabolism, and biosynthesis of antibiotics on KEGG. The upregulated genes of DEGs were mainly enriched in mitotic nuclear division and cell division on BP, and nucleoplasm, nucleus, and cytosol on CC (Table 1).
PPI NETWORK CONSTRUCTION AND MODULE ANALYSIS: The PPI network of DEGs was constructed using Cytoscape (Figure 2A), and the most significant module was obtained (Figure 2B). Enrichment analysis of key modules using DAVID showed that the key module genes were mainly enriched in cell division, mitotic nuclear division, nucleoplasm, nucleus, protein binding, and ATP binding (Table 2).
HUB GENE SELECTION AND ANALYSIS: A total of 12 key genes with degree ≥10 were screened. Their abbreviations, full names, and functions are shown in Table 3. The visualized BP process diagram of key genes is shown in Figure 3A. Hierarchical cluster analysis demonstrated that the UBE2T gene was upregulated in tumor samples and downregulated in normal tissue, in contrast to other genes (Figure 3B). Subsequently, OS survival analysis of these key genes was performed using Kaplan-Meier curve analysis. The results showed that patients with high expression of UBE2T and FBLN5 genes tended to have poorer OS (Figure 4A). The upregulation of UBE2T, CFP, and FBLN5 genes was significantly correlated with poor RFS (Figure 4B), while the high expression of UBE2T, CFP, and FBLN5 was associated with poor progression-free survival (PFS) (Figure 4C). Based on the Kaplan-Meier online database, we found that the high expression of UBE2T and FBLN5 in patients with HCC was involved in poor prognosis in multiple survival analyses, suggesting that they may have an important role in the development of HCC. Oncomine analysis showed that the TOP2A gene was clearly highly expressed in HCC among different data sets, while FBLN5 did not show a significant correlation (Figure 5). In the Wurmbach liver data set, a high level of UBE2T was related to tumor grade, sample type, satellite lesions, and vascular invasion (Figure 6).
For patients with advanced HCC, drug therapy based on sorafenib and lenvatinib is the main palliative treatment, but resistance to the therapy will develop within a few months [12,13]. Therefore, the discovery of new gene therapy sites is crucial. Based on GSE121248, GSE87630, and GSE84598 from the GEO public database, we screened out 296 DEGs in HCC and normal liver tissues, of which 54 were upregulated genes and 242 were downregulated genes. These genes are mainly enriched in the epoxygenase P450 pathway, oxidation-reduction process, and complement activation by GO enrichment analysis, and in mitotic nuclear division, cell division, nucleoplasm, nucleus, and cytosol by KEGG enrichment analysis. It has been reported that cytochrome CYP450 cyclooxygenase can catalyze the epoxidation of unsaturated fatty acids such as arachidonic acid and that arachidonic acid-derived products can effectively promote angiogenesis of lipids and tumor development . Zhang et al reported that the growth and prognosis of cancer can be regulated by affecting the balance of redox reactions . These findings support our enrichment analysis results. A PPI network diagram of DEGs was drawn to explore the interrelationships among different genes, and the key modules were identified by the MCODE plug-in. Enrichment analysis showed that the genes of the key modules were mainly enriched in cell division, mitotic nuclear division, nucleoplasm, nucleus, protein binding, and ATP binding. Neil et al found that the division error of the formation of chromosome bridges of a single cell could rapidly increase the complexity of genes and promote the continuous evolution and subclonal heterogeneity of many cancers . Fletcher et al found that ATP binding plays a considerable role in tumor progression and metastases . All of these reports further verified the enrichment analysis results of the key genes in our research. In our present study, the key genes were extracted from all modules, and the hierarchical cluster analysis of key genes showed that UBE2T gene was significantly highly expressed in HCC. The results of the OS survival curve revealed that the high expression of UBE2T and FBLN5 genes was markedly involved in the poor prognosis of patients with HCC, and the low expression of the CFP gene was apparently involved in the poor prognosis of patients with HCC. The RFS and PFS of these 3 genes showed that the high expression of UBE2T, FBLN5, and CFP were significantly related to the prognosis of patients, yet only UBE2T and FBLN5 were statistically related to the prognosis of HCC in various survival analyses. Meta-analysis of UBE2T and FBLN5 using the Oncomine database showed that UBE2T was evidently higher in HCC tissue than in normal liver tissue. The results produced by the UCSC Cancer Genomics Browser showed that UBE2T was notably positively correlated with sample type, HCC grade, HCC satellite lesions, and vascular invasion, while FBLN5 showed no obvious correlations. Consequently, the high expression of UBE2T may be a poor prognostic factor that promotes the occurrence and development of HCC.
The protein encoded by the UBE2T (ubiquitin-conjugating enzyme E2T) gene catalyzes the covalent attachment of ubiquitin to protein substrates [18,19]. Many studies have confirmed that UBE2T plays an essential role in tumor cell proliferation, invasion, and metastasis [20,21]. The high expression of UBE2T is associated with poor prognosis of a variety of cancers. Hao et al found that in kidney cancer, UBE2T promotes the proliferation of renal cell carcinoma cells by regulating Pi3K/akt signaling . The authors found that UBE2T was highly expressed in lung cancer tissues and cell lines, indicating that UBE2T is involved in the formation of cancer cells . Alagpulinsa et al found that the amplification and overexpression of UBE2T promoted homologous recombination in multiple myeloma . Liu et al found that the overexpression of UBE2T was related to tumor proliferation and invasion and was an independent prognostic factor of HCC. The high expression of UBE2T is associated with higher pathological grades of HCC, advanced TMN staging, poor overall survival rate, and disease-free survival rate, confirming that UBE2T is an independent prognostic factor for the overall survival of patients with HCC . Wei et al found that miR-1305 targeted UBE2T to inhibit the Akt signaling pathway, thereby inhibiting the self-renewal and tumorigenicity of HCC stem cells . In the present study, we found that UBE2T interacted with genes including aurora kinase A (AURKA), TPX2, and CDC20 through the PPI network. AURKA is a mitotic serine/threonine kinase that acts as an oncogene and plays a key role in the development of HCC . Overexpression of TPX2 regulates the cell cycle and thus promotes the development of HCC, ensuring the principal role of the UBE2T gene in the occurrence and development of HCC . In our study, the expression level of UBE2T in HCC tissue was significantly higher than that in normal liver tissue, and the high expression of UBE2T was involved in the poor prognosis of patients with HCC. Ubiquitin signaling is a fundamental eukaryotic regulatory system, controlling diverse cellular functions, and once there are mutations or impairment of the E2 genes, severe disease states can occur . Various studies reported that UBE2T influences tumors through the Fanconi anemia pathway [25,26].
This study has limitations. We have not conducted further experimental verification of the UBE2T gene at the cellular and molecular level; therefore, further research should be done to identify the entire mechanism of UBE2T’s biological effects.
We obtained 296 DEGs in HCC and normal liver tissues through bioinformatic analysis, from which the independent predictor UBE2T was screened out. We determined that UBE2T may be a potential treatment target and prognostic marker for HCC.
FiguresFigure 1. Volcano map and Venn diagram of GSE121248, GSE87630, and GSE84598. (A–C) In the volcano map, the adjusted P value <0.01 and |log FC| ≥1 was considered statistically significant. Upregulated genes are marked in red and downregulated genes are marked in blue. (D) An overlap of 296 differentially expressed genes (DEGs) were finally screened on the Venn diagram. Figure 2. Protein-protein interaction (PPI) network and the most significant modules. (A) PPI network was drawn by Cytoscape. Upregulated genes are marked in red, downregulated genes are marked in blue. (B) The most significant module, including 29 nodes and 388 edges, was obtained from the PPI network. Figure 3. BP analysis and hierarchical cluster analysis of key gene. (A) In the BP graph drawn by the BINGO plug-in in Cytoscape, the size of the node was decided by the number of corresponding genes. (B) Hierarchical cluster analysis graph shows the expression of 12 hub genes in primary hepatocellular carcinoma (HCC) tissue and normal tissue. Upregulated genes are marked in red, downregulated genes are marked in blue. Figure 4. The survival analysis chart was drawn by the Kaplan-Meier plotter online tool; P<0.01 was considered statistically significant. (A) Overall survival curve analysis of key genes. (B) Progression-free survival curve analysis of UBE2T, CFP, and FBLN5. (C) Relapse-free survival curve analysis of UBE2T, CFP, and FBLN5. Figure 5. Analysis of (A) FBLN5 and (B) UBE2T in hepatocellular carcinoma (HCC) and normal liver tissues based on the Oncomine database. (1) Hepatocellular carcinoma vs normal Guichard liver, Nat Genet, 2012. (2) Hepatocellular carcinoma vs normal Guichard liver 2, Nat Genet, 2012. (3) Hepatocellular carcinoma vs normal TCGA liver, no associated paper, 2012. (4) Hepatocellular carcinoma vs normal Wurmbach liver, Hepatology, 2007. Figure 6. The Box plot of gene expression. The relationships between UBE2T and (A) sample type, (B) HCC grade, (C) satellite lesions, and (D) vascular invasion are shown. The relationships between FBLN5 and (E) sample type, (F) HCC grade, (G) satellite lesions, and (H) vascular invasion are shown.
1. Ozakyol A: Global epidemiology of hepatocellular carcinoma (HCC Epidemiology): J Gastrointest Cancer, 2017; 48(3); 238-40
2. El-Serag HB, Marrero JA, Rudolph L, Reddy KR: Diagnosis and treatment of hepatocellular carcinoma: Gastroenterology, 2008; 134(6); 1752-63
3. Xu Y, Xu H, Li M: KIAA1199 promotes sorafenib tolerance and the metastasis of hepatocellular carcinoma by activating the EGF/EGFR-dependent epithelial-mesenchymal transition program: Cancer Lett, 2019; 454(December 2018); 78-89
4. Dunbar CE, High KA, Joung JK: Gene therapy comes of age: Science, 2018; 359(6372); aan4672
5. Dunbar CE, High KA, Joung JK: Gene therapy clinical trials worldwide to 2017: An update (Journal of Gene Medicine, 2018, 20(5):e3015): J Gene Med, 2019; 21(9); 3124
6. Brown ZJ, Fu Q, Ma C: Carnitine palmitoyltransferase gene upregulation by linoleic acid induces CD4+ T cell apoptosis promoting HCC development: Cell Death Dis, 2018; 9(6); 620
7. Shi JY, Ma LJ, Zhang JW: FOXP3 is a HCC suppressor gene and acts through regulating the TGF-β/Smad2/3 signaling pathway: BMC Cancer, 2017; 17(1); 1-10
8. Hao P, Kang B, Li Y: UBE2T promotes proliferation and regulates PI3K/Akt signaling in renal cell carcinoma: Mol Med Rep, 2019; 20(2); 1212-20
9. Wang Y, Leng H, Chen H: Knockdown of UBE2T inhibits osteosarcoma cell proliferation, migration, and invasion by suppressing the PI3K/Akt signaling pathway: Oncol Res, 2016; 24(5); 361-69
10. Liu LL, Zhu JM, Yu XN: UBE2T promotes proliferation via G2/M checkpoint in hepatocellular carcinoma: Cancer Manag Res, 2019; 11; 8359-70
11. Guo J, Wang M, Wang JP, Wu CX: Ubiquitin-conjugating enzyme E2T knockdown suppresses hepatocellular tumorigenesis via inducing cell cycle arrest and apoptosis: World J Gastroenterol, 2019; 9327(43); 6386-403
12. Chen KF, Chen HL, Tai WT: Activation of phosphatidylinositol 3-kinase/Akt signaling pathway mediates acquired resistance to sorafenib in hepatocellular carcinoma cells: J Pharmacol Exp Ther, 2011; 337(1); 155-61
13. Fu R, Jiang S, Li J: Activation of the HGF/c-MET axis promotes lenvatinib resistance in hepatocellular carcinoma cells with high c-MET expression: Med Oncol, 2020; 37(4); 1-7
14. Sausville LN, Williams SM, Pozzi A: Cytochrome P450 epoxygenases and cancer: A genetic and a molecular perspective: Pharmacol Ther, 2019; 196; 183-94
15. Zhang Y, Li J, Dong X: Psat1 regulated oxidation-reduction balance affects the growth and prognosis of epithelial ovarian cancer: Onco Targets Ther, 2020; 13; 5443-53
16. Wei X, You X, Zhang J, Zhou C: MicroRNA-1305 Inhibits the stemness of LCSCs and tumorigenesis by repressing the UBE2T-dependent Akt-signaling pathway: Mol Ther Nucleic Acid, 2019; 16(519); 721-32
17. Fletcher JI, Haber M, Henderson MJ, Norris MD: ABC transporters in cancer: More than just drug efflux pumps: Nat Rev Cancer, 2010; 10(2); 147-56
18. Alpi AF, Pace PE, Babu MM, Patel KJ: Mechanistic insight into site-restricted monoubiquitination of FANCD2 by Ube2t, FANCL, and FANCI: Mol Cell, 2008; 32(6); 767-77
19. Morreale FE, Bortoluzzi A, Chaugule VK: Allosteric targeting of the fanconi anemia ubiquitin-conjugating enzyme Ube2T by fragment screening: J Med Chem, 2017; 60(9); 4093-98
20. Hu W, Xiao L, Cao C: UBE2T promotes nasopharyngeal carcinoma cell proliferation, invasion, and metastasis by activating the AKT/GSKβ/β-catenin pathway: Oncotarget, 2016; 7(12); 15161-72
21. Alagpulinsa DA, Kumar S, Talluri S: Amplification and overexpression of E2 ubiquitin conjugase UBE2T promotes homologous recombination in multiple myeloma: Blood Adv, 2019; 3(23); 3968-72
22. Hao J, Xu A, Xie X: Elevated expression of UBE2T in lung cancer tumors and cell lines: Tumor Biol, 2008; 29(3); 195-203
23. Wang B, Hsu CJ, Chou CH: Variations in the AURKA gene: Biomarkers for the development and progression of hepatocellular carcinoma: Int J Med Sci, 2018; 15(2); 170-75
24. Bai L, Ren Y, Cui T: Overexpression of CDCA5, KIF4A, TPX2, and FOXM1 coregulated cell cycle and promoted hepatocellular carcinoma development: J Comput Biol, 2020; 27(6); 965-74
25. Longerich S, Kwon Y, Tsai M: Regulation of FANCD2 and FANCI monoubiquitination by their interaction and by DNA: Nucleic Acids Res, 2014; 42(9); 5657-70
26. Lewis TW, Barthelemy JR, Virts EL: Deficiency of the Fanconi anemia E2 ubiqitin conjugase UBE2T only partially abrogates Alu-mediated recombination in a new model of homology dependent recombination: Nucleic Acids Res, 2019; 47(7); 3503-20
08 September 2021 : Clinical ResearchAcceptance of COVID-19 Vaccination and Its Associated Factors Among Cancer Patients Attending the Oncology ...
Med Sci Monit In Press; DOI: 10.12659/MSM.932788
16 July 2021 : Review articleSilent Hypoxemia in Patients with COVID-19 Pneumonia: A Review
Med Sci Monit In Press; DOI: 10.12659/MSM.930776
09 July 2021 : Database AnalysisA Simple Clinical Prediction Tool for COVID-19 in Primary Care with Epidemiology: Temperature-Leukocytes-CT...
Med Sci Monit In Press; DOI: 10.12659/MSM.931467
20 Sep 2021 : Clinical ResearchComparison of the Use of Magnetic Resonance Imaging of Partial Anterior Cruciate Ligament Tears Using Maxim...
Med Sci Monit In Press; DOI: 10.12659/MSM.932228
20 Sep 2021 : Clinical ResearchAutomated Boluses and Delayed-Start Timers Prolong Perineural Local Anesthetic Infusions and Analgesia Foll...
Med Sci Monit In Press; DOI: 10.12659/MSM.933190
16 Sep 2021 : Clinical ResearchA Retrospective Evaluation of Operative and Postoperative Outcomes in Patients with Spinal Metastases from ...
Med Sci Monit In Press; DOI: 10.12659/MSM.932995
20 Mar 2020 : Clinical ResearchSocial Capital and Sleep Quality in Individuals Who Self-Isolated for 14 Days During the Coronavirus Diseas...
Med Sci Monit 2020; 26:e923921
15 Apr 2020 : Clinical ResearchPsychological Impact and Coping Strategies of Frontline Medical Staff in Hunan Between January and March 20...
Med Sci Monit 2020; 26:e924171
05 May 2020 : Review articleAn Evidence Based Perspective on mRNA-SARS-CoV-2 Vaccine Development
Med Sci Monit 2020; 26:e924700