Research Paper Volume 14, Issue 9 pp 3782—3800

Transcriptional expressions of hsa-mir-183 predicted target genes as independent indicators for prognosis in bladder urothelial carcinoma

Ming Li1, *, , Da-Ming Xu1, *, , Shu-Bin Lin1, , Zheng-Liang Yang1, , Teng-Yu Xu1, , Jin-Huan Yang1, , Ze-Xin Lin1, , Ze-Kai Huang1, , Jun Yin2,3, ,

  • 1 Division of Urological Surgery, Second Affiliated Hospital of Shantou University Medical College, Shantou, Guangdong, China
  • 2 Division of Hematology, Second Affiliated Hospital of Shantou University Medical College, Shantou, Guangdong, China
  • 3 Department of Clinical Laboratory Medicine, Second Affiliated Hospital of Shantou University Medical College, Shantou, Guangdong, China
* Equal contribution as co-first authors

Received: November 3, 2021       Accepted: April 22, 2022       Published: May 3, 2022      

https://doi.org/10.18632/aging.204040
How to Cite

Copyright: © 2022 Li et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Objective: To uncover novel prognostic and therapeutic targets for BLCA, our study is the first to investigate the role of hsa-mir-183 and its up-regulated predicted target genes in bladder urothelial carcinoma.

Methods: To address this issue, our study explored the roles of hsa-mir-183 predicted target genes in the prognosis of BLCA via UALCAN, Metascape, Kaplan-Meier plotter, Human Protein Atlas, TIMER2.0, cBioPortal and Genomics of Drug Sensitivity in Cancer databases.

Results: High transcriptional expressions of PDCD6, GNG5, PHF6 and MAL2 were markedly relevant to favorable OS in BLCA patients, whereas SLC25A15 and PTDSS1 had opposite expression significance. Additionally, high transcriptional expression of PDCD6, GNG5, PHF6, MAL2, SLC25A15 and PTDSS1 were significantly correlated with BLCA individual cancer stages and molecular subtypes. Furthermore, high mutation rate of PDCD6, MAL2, SLC25A15 and PTDSS1 were observed. Finally, TP53 mutation of PDCD6, GNG5, PHF6, MAL2, SLC25A15 and PTDSS1 has guiding significance for drug selection in BLCA.

Conclusions: PDCD6, GNG5, PHF6, MAL2, SLC25A15 and PTDSS1 could be the advanced independent indicators for prognosis of BLCA patients, and TP53-mutation might be a biomarker for drug option in BLCA patients.

Introduction

Bladder Urothelial Carcinoma (BLCA) is one of the most common urinary system tumors with the unfavorable mortality and increased incidence in the world [1, 2]. The histological features of BLCA are the transformation between urothelium and malignant cells, Involving a series process of polymorphic nuclei, hyperchromatic nuclei, nucleolar herniation, and an increased number of uncontrolled mitoses [3]. Recently, the detailed molecular subtypes of BLCA have been widely explored, namely neuronal, basal squamous, luminal, luminal infiltrated and luminal papillary, which has important prognostic and therapeutic implications. Studies on BLCA have suggested that BLCA in situ without treatment may lead to progression to muscular invasive disease (occur in 50%) and disease relapse (occur in 90%) [4, 5]. New insights on the pathogenesis, invasion and metastasis of BLCA have been ongoing in recent studies, however, the research on the gene molecular basis of bladder cancer has been limited. Providing crucial insight into the gene level characteristics of BLCA will help to better reveal the mechanism of pathogenesis and discover more significant diagnostic or prognostic biomarkers and therapeutic options.

MicroRNA (miRNA), a small noncoding RNA with the function of negative regulating gene expression, plays an important role in the occurrence and development of various tumors [68]. The functions of miRNA are not to directly block the expression of genes, but to make a slightly and specifically adjustment to the expression of target genes in post-transcriptional control [9, 10]. MiRNA plays an important regulatory role in the biological processes of normal human cells, including cell division cycle, cell growth, cell apoptosis, cell differentiation and cell response. Almost all tumors have different types of miRNA expression, which are either up-regulated or down-regulated, thus providing new opportunities for miRNA and their predicted target genes to be used as the potential diagnostic or prognostic biomarkers. Hsa-mir-183 has been shown to frequently involve methylation events in human hepatocellular carcinoma and was an unfavorable factor for survival [11]. Huang et al., found that hsa-mir-183 was an overexpressed miRNA in non-small-cell lung cancer and was positively correlated with tumor invasiveness [12]. Additionally, hsa-mir-183 has been reported to be markedly correlated with overall survival for colon adenocarcinoma patients [13]. However, the role of hsa-mir-183 and its predicted target genes in BLCA remains unclear.

To uncover novel prognostic and therapeutic targets for BLCA, our study is the first to investigate the role of hsa-mir-183 and its up-regulated predicted target genes and evaluate the significance of TP53-mutation related drug option in BLCA patients.

Materials and Methods

To estimate the prognostic and therapeutic significance of hsa-mir-183 and its up-regulated predicted target genes in BLCA patients, the clinical data, transcriptional expression data of miRNA and predicted target genes were collected from UALCAN (http://ualcan.path.uab.edu). The Metascape was applied for GO and KEGG analysis of predicted target genes (https://metascape.org). The Overall Survival (OS) data of mRNA expression of predicted target genes were collected from Kaplan-Meier plotter (http://www.kmplot.com/). The protein expression data were collected from Human Protein Atlas (https://www.proteinatlas.org). The data of pan cancer analysis of predicted target genes were collected from TIMER2.0 (http://timer.cistrome.org). Gene mutations of predicted target genes in BLCA and their OS and Progression Free Survival (PFS) were analyzed by cBioPortal (http://www.cbioportal.org). Finally, the drug sensitivity and resistance of TP53 mutation in BLCA were analyzed by Genomics of Drug Sensitivity in Cancer (https://www.cancerrxgene.org).

Ethics approval

The study protocols were conducted according to the principles of the Declaration of Helsinki and were approved by the Scientific and Medical Ethical Committee of the Second Affiliated Hospital of Shantou University Medical College.

Results

Over-expression of different miRNA in BLCA patients

The result of top 10 over-expressed miRNA in BLCA patients were shown by TCGA miRNA analysis (Data from UALCAN) (Figure 1). The expression of hsa-mir-183 in patients with BLCA were significantly increased compared to normal tissue (P = 1.62E-12). The expression of hsa-mir-183 were also significantly increased based on patient’s gender or age (Figure 2). There were significant differences in the Comparison of Normal-vs-Male (P = 1.62E-12), Normal-vs-Female (P = 1.63E-8), Normal-vs-Age (41–60 Yrs) (P < 1E-12), Normal-vs-Age (61–80 Yrs) (P < 1E-12) and Normal-vs-Age (81–100 Yrs) (P = 7.55E-9).

Top 10 over-expressed miRNAs in BLCA.

Figure 1. Top 10 over-expressed miRNAs in BLCA.

Expression of hsa-mir-183 in BLCA based on patient’s gender and age.

Figure 2. Expression of hsa-mir-183 in BLCA based on patient’s gender and age.

Expression of predicted target genes of hsa-mir-183 in BLCA patients

There were 247 predicted genes of hsa-mir-183 (Data from TargetScan, microRNA.org and miRDB), including 15 up-regulated genes, 74 down-regulated genes and 158 unchanged genes in BLCA. 15 up-regulated genes were selected to study in our research (Table 1). The expression of all up-regulated predicted genes of has-mir-183 in patients with BLCA were significantly increased compared to normal tissue (Figure 3).

Table 1. Up-expressed predicted target genes of hsa-mir-183 in BLCA (TargetScan, microRNA.org and miRDB).

Target geneDescriptionExpression in BLCA
ENY2ENY2, transcription and export complex 2 subunitUp-expression
PLAG1PLAG1 zinc fingerUp-expression
PDCD6programmed cell death 6Up-expression
GNG5G protein subunit gamma 5Up-expression
SLC6A6solute carrier family 6 member 6Up-expression
CELSR3cadherin EGF LAG seven-pass G-type receptor 3Up-expression
SLC22A23solute carrier family 22 member 23Up-expression
SPATS2spermatogenesis associated serine rich 2Up-expression
CDK5R1cyclin dependent kinase 5 regulatory subunit 1Up-expression
PHF6PHD finger protein 6Up-expression
MAL2mal, T cell differentiation protein 2Up-expression
ATP2C1ATPase secretory pathway Ca2+ transporting 1Up-expression
NRASNRAS proto-oncogene, GTPaseUp-expression
SLC25A15solute carrier family 25 member 15Up-expression
PTDSS1phosphatidylserine synthase 1Up-expression
Transcriptional expression of predicted target genes of hsa-mir-183 in BLCA based on TCGA sample types. (**P *P

Figure 3. Transcriptional expression of predicted target genes of hsa-mir-183 in BLCA based on TCGA sample types. (**P < 0.001, *P < 0.01).

Predicted functions and pathways of up-expressed predicted target genes of hsa-mir-183 and their 10 frequently altered neighbor genes

Functions and pathways of up-expressed predicted target genes of hsa-mir-183 and their 10 frequently altered neighbor genes were analyzed by GO and KEGG in Metascape. As the Enrichment heatmap showed that pathways included hsa04726: Serotonergic synapse, hsa05219: Bladder cancer, GO:0007169: transmembrane receptor protein tyrosine kinase signaling pathway, GO:0048471: perinuclear region of cytoplasm, GO:0008324: cation transmembrane transporter activity, GO:1990234: transferase complex and GO:1905477: positive regulation of protein localization to membrane (Figure 4).

Predicted functions and pathways of up-expressed predicted target genes of hsa-mir-183 and their 10 frequently altered neighbor genes (Metascape).

Figure 4. Predicted functions and pathways of up-expressed predicted target genes of hsa-mir-183 and their 10 frequently altered neighbor genes (Metascape).

Prognostic value of mRNA expression of up-expressed predicted target genes of hsa-mir-183 in BLCA patients

Prognostic value of mRNA expression of up-expressed predicted target genes of hsa-mir-183 were analyzed by Kaplan-Meier plotter (Figure 5). It was showed that mRNA expressions of PDCD6, GNG5, PHF6, MAL2, SLC25A15 and PTDSS1 were significantly associated with BLCA patients’ prognosis. According to the database analysis form Kaplan-Meier plotter, PDCD6 (HR = 0.55, 95% CI: 0.4–0.75, and P = 0.00013), GNG5 (HR = 0.73, 95% CI: 0.54–0.89, and P = 0.034), PHF6 (HR = 0.68, 95% CI: 0.51–0.92, and P = 0.011) and MAL2 (HR = 0.68, 95% CI: 0.5–0.91, and P = 0.01) seem to be the favorable factor for BLCA patients’ prognosis. SLC25A15 (HR = 1.68, 95% CI: 1.24–2.27, and P = 0.00062) and PTDSS1 (HR = 1.41, 95% CI: 1.01–1.96, and P = 0.041) more likely to be the unfavorable factor for BLCA patients’ prognosis. In addition, ENY2, PLAG1, SLC6A6, CELSR3, SLC22A23, SPATS2, CDK5R1, ATP2C1 and NRAS were not prognostic in BLCA patients.

Prognostic value of transcriptional expression of up-expressed predicted target genes of hsa-mir-183 in BLCA (Kaplan-Meier Plotter).

Figure 5. Prognostic value of transcriptional expression of up-expressed predicted target genes of hsa-mir-183 in BLCA (Kaplan-Meier Plotter).

Protein expression and pan cancer analysis of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15 and PTDSS1) in BLCA patients

Protein expression of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15 and PTDSS1) in BLCA were explored by Human Protein Atlas (Figure 6). The proteins of PDCD6, GNG5 and PHF6 were not or low expressed in BLCA tissues, whereas medium expressions of them were found in normal urinary bladder tissues. It is noteworthy that the protein expression of MAL2 and SLC25A15 could not be detected in both normal urinary bladder tissues and BLCA tissues. Additionally, the protein expression of PTDSS1 in tumor tissues was stronger than that in normal tissues. Furthermore, we performed a pan-cancer analysis of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15 and PTDSS1) by using TIMER (Figure 7).

Representative immunohistochemistry images of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15 and PTDSS1) in normal urinary bladder tissues and BLCA tissues (Human Protein Atlas).

Figure 6. Representative immunohistochemistry images of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15 and PTDSS1) in normal urinary bladder tissues and BLCA tissues (Human Protein Atlas).

Pan cancer analysis of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15 and PTDSS1) (TIMER). (***P **P *P

Figure 7. Pan cancer analysis of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15 and PTDSS1) (TIMER). (***P < 0.001, **P < 0.01, *P < 0.05).

Relationship between mRNA expression of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15, PTDSS1) and clinicopathological parameters in BLCA patients

The association of mRNA expression of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15, PTDSS1) with clinicopathological parameters of BLCA patients were analyzed by UALCAN, containing BLCA individual cancer stages and molecular subtypes (Figures 8 and 9). Figure 8 showed that the mRNA expression of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15, PTDSS1) were significantly correlated with BLCA individual cancer stages. PDCD6, GNG5, PHF6 and MAL2 as favorable factor for BLCA patients, the high mRNA expressions of them tended to be in stage 1 or 2, whereas the high mRNA expressions of unfavorable factor of SLC25A15 and PTDSS1 tended to be in stage 3 or 4. Similarly, Figure 9 showed that the mRNA expression of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15, PTDSS1) were significantly correlated with BLCA molecular subtypes. The high mRNA expressions of PDCD6, GNG5, PHF6 and MAL2 tended to be the molecular subtypes of luminal and luminal papillary, whereas the high mRNA expressions of SLC25A15 and PTDSS1 tended to be the molecular subtypes of neuronal and basal squamous.

Relationship between mRNA expression of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15, PTDSS1) and BLCA individual cancer stages. (***P **P *P

Figure 8. Relationship between mRNA expression of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15, PTDSS1) and BLCA individual cancer stages. (***P < 0.001, **P < 0.01, *P < 0.05).

Relationship between mRNA expression of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15, PTDSS1) and BLCA molecular subtypes. (***P **P *P

Figure 9. Relationship between mRNA expression of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15, PTDSS1) and BLCA molecular subtypes. (***P < 0.001, **P < 0.01, *P < 0.05).

Gene mutations of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15, PTDSS1) in BLCA

Gene mutations of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15, PTDSS1) in BLCA and their OS and PFS were analyzed by cBioPortal (Figure 10). Figure 10A showed that PDCD6, MAL2, SLC25A15, PTDSS1 were prone to mutation in BLCA patients. Figure 10B and 10C revealed the association between the gene mutation and BLCA patients’ OS and DFS. In OS analyze, median months overall (95% CI) of unaltered group was 36.43, while it was significantly shorter in PHF6 (1.87) and SLC25A15 (17.62). Similarly, in PFS analyze, Median Months Progression Free (95% CI) of unaltered group was 32.58, while it was significantly shorter in PHF6 (1.87), MAL2 (5.52) and SLC25A15 (17.52).

Gene mutations of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15, PTDSS1) and their association with OS and PFS in BLCA patients (A: Genetic alterations; B: Overall survival; C: Progress free survival).

Figure 10. Gene mutations of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15, PTDSS1) and their association with OS and PFS in BLCA patients (A: Genetic alterations; B: Overall survival; C: Progress free survival).

TP53 mutation of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15, PTDSS1) and its association with drug selection in BLCA

TP53 mutation of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15, PTDSS1) were analyzed by UALCAN. The expression of TP53 mutation of GNG5, MAL2, SLC25A15, PTDSS1 in patients with BLCA were significantly increased compared to non-mutation group, and all predicted target genes were markedly increased compared to normal group (Figure 11). In addition, types of TP53 mutation were analyzed by cBioPortal (Figure 12) from Bladder Cancer (MSK/TCGA, 2020). Mutation frequency of TP53 was 51.7% from NM_000546 | ENST00000269305 CCDS11118 | P53_HUMAN.

TP53 mutation of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15, PTDSS1) in BLCA.

Figure 11. TP53 mutation of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15, PTDSS1) in BLCA.

As was showed in Figure 12, TP53 driver had the highest probability of mutation, including Missense (197), Truncating (80) and Splice (17). Furthermore, the drug sensitivity and resistance of TP53 mutation in BLCA were analyzed by Genomics of Drug Sensitivity in Cancer (Figure 13). Nutlin-3a (−) drug target was MDM2 and Nutlin-3a (−) drug target pathway was p53 pathway in BLCA. As was showed in Figure 13 and Table 2, Nutlin-3a (−) had a significant increased resistance in TP53 mutation (P = 2.12e-06) of BLCA patients and the IC50 of Nutlin-3a (−) was significantly increased in TP53 mutation group (P < 0.001). TWS119 drug target was GSK3 and TWS119 drug target pathway was WNT signaling in BLCA. As was showed in Figure 13 and Table 2, TWS119 had a significant increased sensitivity in TP53 mutation (P = 5.95e-06) of BLCA patients and the IC50 of TWS119 was significant decreased in TP53 mutation group (P < 0.001). Mitomycin-C drug target was DNA crosslinker and Mitomycin-C drug target pathway was DNA replication in BLCA. As was showed in Figure 13 and Table 2, Mitomycin-C had a significant increased sensitivity in TP53 mutation (P = 0.000322) of BLCA patients and the IC50 of Mitomycin-C was significant decreased in TP53 mutation group (P = 0.001).

Types of TP53 mutation from bladder cancer (MSK/TCGA, 2020).

Figure 12. Types of TP53 mutation from bladder cancer (MSK/TCGA, 2020).

The IC50 of TWS119, Mitomycin-C and Nutlin-3a (−) in TP53 mutation of BLCA patients.

Figure 13. The IC50 of TWS119, Mitomycin-C and Nutlin-3a (−) in TP53 mutation of BLCA patients.

Table 2. Association between TP53 mutation and drug option of BLCA.

DrugDrug TargetIC50 effectEffect sizeP-valueFDR%Tissue
TWS119GSK3Increased sensitivity−3.635.95e-060.00714BLCA
Mitomycin-CDNA crosslinkerIncreased sensitivity−2.530.0003220BLCA
Nutlin-3a (−)MDM2Increased resistance4.082.12e-060.00255BLCA

Discussion

BLCA was the most common histopathological type of bladder cancer, and its occurrence and development were affected by multiple complex factors such as gene transcription, gene regulation, gene mutation and epigenetic inheritance [14, 15]. Each miRNA had a series of related target genes, and abnormal expression of miRNA often leaded to abnormal transcription of target genes, which was related to the occurrence of various cancers [16]. Although hsa-mir-183 have been reported to play an essential role in liver cancer, lung cancer and other tumors, the specific role of hsa-mir-183 and its predicted target genes in BLCA remain to be explored. In our study, the transcriptional expressions, mutations and Prognosis of hsa-mir-183 and its predicted target genes in BLCA have been elucidated.

Our study found that up-expressed hsa-mir-183 and over transcriptional expressions of hsa-mir-183 predicted target genes (ENY2, PLAG1, PDCD6, GNG5, SLC6A6, CELSR3, SLC22A23, SPATS2, CDK5R1, PHF6, MAL2, ATP2C1, NRAS, SLC25A15 and PTDSS1) were observed in BLCA patients. Moreover, predicted functions and pathways of up-expressed predicted target genes of hsa-mir-183 and their 10 frequently altered neighbor genes indicated a link with bladder cancer. Besides, high transcriptional expressions of PDCD6, GNG5, PHF6 and MAL2 were markedly relevant to favorable OS in BLCA patients, whereas high transcriptional expression of SLC25A15 and PTDSS1 were markedly relevant to unfavorable OS in BLCA patients. Additionally, ENY2, PLAG1, SLC6A6, CELSR3, SLC22A23, SPATS2, CDK5R1, ATP2C1 and NRAS were not prognostic in BLCA patients. However, protein expression of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15 and PTDSS1) were inconsistent with the transcriptional expression. The mRNA levels were significantly increased than that observed in protein levels. There are two possible explanations for this phenomenon. miRNA do not limit gene translation but lead to rapid destruction of newborn peptides [17]. Alternatively, miRNA may negatively affect the initiation, elongation and termination of transcription, while the number of ribosomes remains unchanged, resulting in a decreased rate of protein completion [17]. Multivariate analysis indicated that transcriptional expression of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15 and PTDSS1) were significantly correlated with BLCA individual cancer stages and molecular subtypes. Furthermore, high mutation rates of PDCD6, MAL2, SLC25A15 and PTDSS1 were observed in BLCA patients and some genetic alteration were significantly relevant to OS and PFS. Finally, TP53 mutation of predicted target genes (PDCD6, GNG5, PHF6, MAL2, SLC25A15 and PTDSS1) has guiding significance for drug selection in BLCA.

PDCD6

Programmed cell death 6, can encode the calcium binding protein, which is belonging to the penta-EF-hand protein family. High expression of PDCD6 was relevant to unfavorable prognosis in colorectal cancer patients via activating c-Raf, MEK and ERK pathway [18]. Inhibition of PDCD6 expression can reduce the rate of tumor metastasis [19]. PDCD6 can induce migration and invasion of tumor cell in ovarian cancer causing unfavorable clinical outcome to patients [20]. However, the expression of GG genotype of PDCD6 has been reported to reduce the risk of lung cancer [21]. Similarly, up-expressed PDCD6 played the role of inhibiting cell growth in cervical carcinoma cells [22]. These results suggested that PDCD6 may play different roles in tumorigenesis and prognosis. In our study, the transcriptional expression of PDCD6 was significantly increased in BLCA and suggested a favorable prognosis. The function of MAPK pathway is to induce cell hyperplasia and invasion, but it can also increase the occurrence of cell apoptosis or inhibit cell invasion. Different effects of PDCD6 predominantly depend on the signal strength and the microenvironment of abnormal signal activation [23]. In addition, the downstream MAPK signal is mainly regulated by the upstream c-Raf signal. It has been proved that PDCD6-c-Raf-MAPK pathway exists in tumors in recent research [18]. Therefore, we speculated that tumor microenvironment or signal specificity lead to different roles of PDCD6 in tumors.

GNG5

G protein subunit gamma 5, can regulate the information from a variety of receptors. High expression of GNG5 was relevant to unfavorable prognosis in gliomas [24]. At present, there are few studies on GNG5 in other tumors. It was reported that GNG5 may regulate the occurrence and development of tumors by stimulating the PI3K/Akt signal pathway [25]. Akt is activated by upstream PI3K, leading to phosphorylation, which finally promotes tumor cell growth. To our knowledge, this is the first study to explore the role of GNG5 in BLCA. Transcriptional expression of GNG5 was found to be significantly associated with BLCA individual cancer stages, BLCA molecular subtypes and OS. Of note, higher expression of GNG5 indicated a longer possibility of survival. Therefore, GNG5 may act as an advanced indicator of the stages, types and prognosis for BLCA which has not been reported in previous researches.

PHF6

PHD finger protein 6, has a function of transcriptional regulation by encoding domain with 2 PHD-type zinc fingers in protein. The high expression of PHF6 can play a role of tumor suppressor in patients with acute lymphoblastic leukemia and as a favorable prognostic factor [2628]. However, other study has come up with an inconsistent view that PHF6 was involved in proliferation, migration, apoptosis and metabolism in hepatocellular carcinoma and its deficiency can inhibit above mentioned processes in tumor cells [29]. Wang et al., found that PHF6-null can inhibit cell proliferation and arrest cells in the G(2)/M phase, suggesting that loss of PHF6 may result in accumulation of DNA damage [30]. Hsu et al., reported that the deficiency of PHF6 can enhance cell self-renewal and carcinogenic potential [31]. Together, these findings indicated that whether PHF6 was a tumor suppressor or a tumor protein may depend on the specific context in which it acted. In our study, high transcriptional expression of PHF6 was significantly related to high probability of survival in BLCA which were consistent with the research results in acute lymphoblastic leukemia.

MAL2

Mal, T cell differentiation protein 2, can encode multitransmembrane protein, which is a member of MAL proteolipid family. MAL2 was over-expressed in many human cancers, including breast, stomach, liver, pancreatic, colon, kidney, prostate and ovarian cancers [3236]. Zheng et al., found that MAL2 was markedly relevant to grade, stage and the Gleason score of prostate cancer and it can promote tumor progression via the Notch pathway [33]. Fang et al., showed that Inhibition of MAL2 expression can significantly increase the cytotoxicity of CD8+ T lymphocytes, thereby inhibiting the growth and metastasis of breast cancer [35]. Zhang et al., reported that overexpressed MAL2 can phosphorylate ERK1/2, which may promote pancreatic cancer progression [37]. In BLCA patients, MAL2 was also highly expressed, but it was associated with a good prognosis. At present, the specific signaling pathway that MAL2 participates in bladder cancer is not clear, so we speculate that up-expressed MAL2 may not be the initiating factor of tumorigenesis, but may lead to changes in cell morphology, signal transduction, migration and metabolism by altering protein distribution, thus promoting tumor transformation [38], which is bidirectional.

SLC25A15

Solute carrier family 25 member 15, can encode the protein transporter ornithine, which is a member of mitochondrial carrier family. Ji et al., found that SLC25A15 was up-expressed in melanoma patients and was negatively associated with OS and DFS [39]. SLC25A15 was involved in the regulation of cell proliferation and apoptosis in cells and over-expression can reverse the role of tumor suppressor genes in prostate cancer [40]. At present, most of the studies on SLC25A15 are focused on hyperornithinemia-hyperammonemia-homocitrullinuria [41], and there are few reports on tumor. SLC25A15 is primarily responsible for transporting ornithine from the cytoplasm to the mitochondrial matrix, which is an important step in regulating the urea cycle [42]. Whether SLC25A15 has other potential signaling pathways involved in tumor progression. However, by comparing recent researches with our present study, it can be preliminarily speculated that SLC25A15 is an unfavorable prognostic signal for BLCA, which may act as a new biomarker to assist the diagnosis and prognosis of BLCA.

PTDSS1

Phosphatidylserine synthase 1, has a function to promote the synthesis of phosphatidylserine. Li et al., found that over-expression of PTDSS1 in astroglioma involved in lipid metabolism [43]. Wang et al., shown that PTDSS1 was an oncogene of lung adenocarcinoma, and overexpression was significantly associated with low survival [44]. PTDSS1 can participate in the biosynthesis of phosphatidylserine by directing the synthesis of a key enzyme, phosphatidylserine synthase, which was an important part of phosphatidylserine metabolism [44]. PTDSS1-mediated phosphatidylserine signaling has been shown to be severely disrupted in metabolism of tumor cell and pathogenesis of autoimmune disorders [45]. Tumor progression and metastasis are greatly influenced by the interaction between tumor cells and their living microenvironment. TAM receptors exist on the surface of tumor cells, and an environment rich in phosphatidylserine can provide a platform for activation of TAM receptors [46]. In our study, high transcriptional expression of PTDSS1 was significantly related to low probability of survival in BLCA which indicated that PTDSS1-mediated phosphatidylserine signaling was one of the pathogenesis of BLCA.

Conclusion

To sum up, we concluded that PDCD6, GNG5, PHF6, MAL2, SLC25A15 and PTDSS1 could be the advanced independent indicators for prognosis of BLCA patients, and TP53-mutation might be a biomarker for drug option in BLCA patients.

Author Contributions

All authors made a significant contribution to the work reported, whether that is in the conception, study design, execution, acquisition of data, analysis and interpretation, or in all these areas; took part in drafting, revising or critically reviewing the article; gave final approval of the version to be published; have agreed on the journal to which the article has been submitted; and agree to be accountable for all aspects of the work.

Conflicts of Interest

The authors declare no conflicts of interest related to this study.

References

  • 1. Wullweber A, Strick R, Lange F, Sikic D, Taubert H, Wach S, Wullich B, Bertz S, Weyerer V, Stoehr R, Breyer J, Burger M, Hartmann A, et al. Bladder Tumor Subtype Commitment Occurs in Carcinoma In Situ Driven by Key Signaling Pathways Including ECM Remodeling. Cancer Res. 2021; 81:1552–66. https://doi.org/10.1158/0008-5472.CAN-20-2336 [PubMed]
  • 2. Feng RM, Zong YN, Cao SM, Xu RH. Current cancer situation in China: good or bad news from the 2018 Global Cancer Statistics? Cancer Commun (Lond). 2019; 39:22. https://doi.org/10.1186/s40880-019-0368-6 [PubMed]
  • 3. Akhtar M, Al-Bozom IA, Ben Gashir M, Taha NM, Rashid S, Al-Nabet ADM. Urothelial Carcinoma In Situ (CIS): New Insights. Adv Anat Pathol. 2019; 26:313–9. https://doi.org/10.1097/PAP.0000000000000239 [PubMed]
  • 4. Williamson SR, Montironi R, Lopez-Beltran A, MacLennan GT, Davidson DD, Cheng L. Diagnosis, evaluation and treatment of carcinoma in situ of the urinary bladder: the state of the art. Crit Rev Oncol Hematol. 2010; 76:112–26. https://doi.org/10.1016/j.critrevonc.2010.01.005 [PubMed]
  • 5. Casey RG, Catto JW, Cheng L, Cookson MS, Herr H, Shariat S, Witjes JA, Black PC. Diagnosis and management of urothelial carcinoma in situ of the lower urinary tract: a systematic review. Eur Urol. 2015; 67:876–88. https://doi.org/10.1016/j.eururo.2014.10.040 [PubMed]
  • 6. Calin GA, Sevignani C, Dumitru CD, Hyslop T, Noch E, Yendamuri S, Shimizu M, Rattan S, Bullrich F, Negrini M, Croce CM. Human microRNA genes are frequently located at fragile sites and genomic regions involved in cancers. Proc Natl Acad Sci U S A. 2004; 101:2999–3004. https://doi.org/10.1073/pnas.0307323101 [PubMed]
  • 7. Melo SA, Esteller M. Dysregulation of microRNAs in cancer: playing with fire. FEBS Lett. 2011; 585:2087–99. https://doi.org/10.1016/j.febslet.2010.08.009 [PubMed]
  • 8. Chan JJ, Tay Y. Noncoding RNA:RNA Regulatory Networks in Cancer. Int J Mol Sci. 2018; 19:E1310. https://doi.org/10.3390/ijms19051310 [PubMed]
  • 9. Fabbri M, Calin GA. Epigenetics and miRNAs in human cancer. Adv Genet. 2010; 70:87–99. https://doi.org/10.1016/B978-0-12-380866-0.60004-6 [PubMed]
  • 10. Kunej T, Godnic I, Ferdin J, Horvat S, Dovc P, Calin GA. Epigenetic regulation of microRNAs in cancer: an integrated review of literature. Mutat Res. 2011; 717:77–84. https://doi.org/10.1016/j.mrfmmm.2011.03.008 [PubMed]
  • 11. Anwar SL, Krech T, Hasemeier B, Schipper E, Schweitzer N, Vogel A, Kreipe H, Buurman R, Skawran B, Lehmann U. hsa-mir-183 is frequently methylated and related to poor survival in human hepatocellular carcinoma. World J Gastroenterol. 2017; 23:1568–75. https://doi.org/10.3748/wjg.v23.i9.1568 [PubMed]
  • 12. Huang R, Cho WC, Sun Y, Katie Chan KH. The Lung Cancer Associated MicroRNAs and Single Nucleotides Polymorphisms: a Mendelian Randomization Analysis. Annu Int Conf IEEE Eng Med Biol Soc. 2020; 2020:2346–52. https://doi.org/10.1109/EMBC44109.2020.9176344 [PubMed]
  • 13. Wang WJ, Li HT, Yu JP, Han XP, Xu ZP, Li YM, Jiao ZY, Liu HB. A Competing Endogenous RNA Network Reveals Novel Potential lncRNA, miRNA, and mRNA Biomarkers in the Prognosis of Human Colon Adenocarcinoma. J Surg Res. 2019; 235:22–33. https://doi.org/10.1016/j.jss.2018.09.053 [PubMed]
  • 14. Kim WJ, Bae SC. Molecular biomarkers in urothelial bladder cancer. Cancer Sci. 2008; 99:646–52. https://doi.org/10.1111/j.1349-7006.2008.00735.x [PubMed]
  • 15. Han Y, Chen J, Zhao X, Liang C, Wang Y, Sun L, Jiang Z, Zhang Z, Yang R, Chen J, Li Z, Tang A, Li X, et al. MicroRNA expression signatures of bladder cancer revealed by deep sequencing. PLoS One. 2011; 6:e18286. https://doi.org/10.1371/journal.pone.0018286 [PubMed]
  • 16. Mishra S, Yadav T, Rani V. Exploring miRNA based approaches in cancer diagnostics and therapeutics. Crit Rev Oncol Hematol. 2016; 98:12–23. https://doi.org/10.1016/j.critrevonc.2015.10.003 [PubMed]
  • 17. Valencia-Sanchez MA, Liu J, Hannon GJ, Parker R. Control of translation and mRNA degradation by miRNAs and siRNAs. Genes Dev. 2006; 20:515–24. https://doi.org/10.1101/gad.1399806 [PubMed]
  • 18. Wang X, Wu F, Wang H, Duan X, Huang R, Tuersuntuoheti A, Su L, Yan S, Zhao Y, Lu Y, Li K, Yao J, Luo Z, et al. PDCD6 cooperates with C-Raf to facilitate colorectal cancer progression via Raf/MEK/ERK activation. J Exp Clin Cancer Res. 2020; 39:147. https://doi.org/10.1186/s13046-020-01632-9 [PubMed]
  • 19. Zhang L, Chen X, Liu B, Han J. MicroRNA-124-3p directly targets PDCD6 to inhibit metastasis in breast cancer. Oncol Lett. 2018; 15:984–90. https://doi.org/10.3892/ol.2017.7358 [PubMed]
  • 20. Zhou B, Zhang P, Tang T, Zhang K, Wang Y, Song Y, Liao H, Zhang L. Prognostic value of PDCD6 polymorphisms and the susceptibility to bladder cancer. Tumour Biol. 2014; 35:7547–54. https://doi.org/10.1007/s13277-014-2010-1 [PubMed]
  • 21. Shen TC, Chang WS, Hsia TC, Li HT, Chen WC, Tsai CW, Bau DT. Contribution of programmed cell death 6 genetic variations, gender, and smoking status to lung cancer. Onco Targets Ther. 2019; 12:6237–44. https://doi.org/10.2147/OTT.S205544 [PubMed]
  • 22. Liu X. Up-regulation of miR-20a by HPV16 E6 exerts growth-promoting effects by targeting PDCD6 in cervical carcinoma cells. Biomed Pharmacother. 2018; 102:996–1002. https://doi.org/10.1016/j.biopha.2018.03.154 [PubMed]
  • 23. Burotto M, Chiou VL, Lee JM, Kohn EC. The MAPK pathway across different malignancies: a new perspective. Cancer. 2014; 120:3446–56. https://doi.org/10.1002/cncr.28864 [PubMed]
  • 24. Yang B, Han ZY, Wang WJ, Ma YB, Chu SH. GNG5 is an unfavourable independent prognostic indicator of gliomas. J Cell Mol Med. 2020; 24:12873–8. https://doi.org/10.1111/jcmm.15923 [PubMed]
  • 25. Jin R, Jin YY, Tang YL, Yang HJ, Zhou XQ, Lei Z. GPNMB silencing suppresses the proliferation and metastasis of osteosarcoma cells by blocking the PI3K/Akt/mTOR signaling pathway. Oncol Rep. 2018; 39:3034–40. https://doi.org/10.3892/or.2018.6346 [PubMed]
  • 26. Alexander TB, Gu Z, Iacobucci I, Dickerson K, Choi JK, Xu B, Payne-Turner D, Yoshihara H, Loh ML, Horan J, Buldini B, Basso G, Elitzur S, et al. The genetic basis and cell of origin of mixed phenotype acute leukaemia. Nature. 2018; 562:373–9. https://doi.org/10.1038/s41586-018-0436-0 [PubMed]
  • 27. Xiao W, Bharadwaj M, Levine M, Farnhoud N, Pastore F, Getta BM, Hultquist A, Famulare C, Medina JS, Patel MA, Gao Q, Lewis N, Pichardo J, et al. PHF6 and DNMT3A mutations are enriched in distinct subgroups of mixed phenotype acute leukemia with T-lineage differentiation. Blood Adv. 2018;2(23):3526–3539. Blood Adv. 2019; 3:956. https://doi.org/10.1182/bloodadvances.2018023531 [PubMed]
  • 28. Kurzer JH, Weinberg OK. PHF6 Mutations in Hematologic Malignancies. Front Oncol. 2021; 11:704471. https://doi.org/10.3389/fonc.2021.704471 [PubMed]
  • 29. Yu Q, Zhou J, Jian Y, Xiu Z, Xiang L, Yang D, Zeng W. MicroRNA-214 suppresses cell proliferation and migration and cell metabolism by targeting PDK2 and PHF6 in hepatocellular carcinoma. Cell Biol Int. 2019. [Epub ahead of print]. https://doi.org/10.1002/cbin.11207 [PubMed]
  • 30. Wang J, Leung JW, Gong Z, Feng L, Shi X, Chen J. PHF6 regulates cell cycle progression by suppressing ribosomal RNA synthesis. J Biol Chem. 2013; 288:3174–83. https://doi.org/10.1074/jbc.M112.414839 [PubMed]
  • 31. Hsu YC, Chen TC, Lin CC, Yuan CT, Hsu CL, Hou HA, Kao CJ, Chuang PH, Chen YR, Chou WC, Tien HF. Phf6-null hematopoietic stem cells have enhanced self-renewal capacity and oncogenic potentials. Blood Adv. 2019; 3:2355–67. https://doi.org/10.1182/bloodadvances.2019000391 [PubMed]
  • 32. López-Coral A, Del Vecchio GJ, Chahine JJ, Kallakury BV, Tuma PL. MAL2-Induced Actin-Based Protrusion Formation is Anti-Oncogenic in Hepatocellular Carcinoma. Cancers (Basel). 2020; 12:E422. https://doi.org/10.3390/cancers12020422 [PubMed]
  • 33. Zheng C, Wang J, Zhang J, Hou S, Zheng Y, Wang Q. Myelin and lymphocyte protein 2 regulates cell proliferation and metastasis through the Notch pathway in prostate adenocarcinoma. Transl Androl Urol. 2021; 10:2067–77. https://doi.org/10.21037/tau-21-244 [PubMed]
  • 34. Eguchi D, Ohuchida K, Kozono S, Ikenaga N, Shindo K, Cui L, Fujiwara K, Akagawa S, Ohtsuka T, Takahata S, Tokunaga S, Mizumoto K, Tanaka M. MAL2 expression predicts distant metastasis and short survival in pancreatic cancer. Surgery. 2013; 154:573–82. https://doi.org/10.1016/j.surg.2013.03.010 [PubMed]
  • 35. Fang Y, Wang L, Wan C, Sun Y, Van der Jeught K, Zhou Z, Dong T, So KM, Yu T, Li Y, Eyvani H, Colter AB, Dong E, et al. MAL2 drives immune evasion in breast cancer by suppressing tumor antigen presentation. J Clin Invest. 2021; 131:140837. https://doi.org/10.1172/JCI140837 [PubMed]
  • 36. Li J, Li Y, Liu H, Liu Y, Cui B. The four-transmembrane protein MAL2 and tumor protein D52 (TPD52) are highly expressed in colorectal cancer and correlated with poor prognosis. PLoS One. 2017; 12:e0178515. https://doi.org/10.1371/journal.pone.0178515 [PubMed]
  • 37. Zhang B, Xiao J, Cheng X, Liu T. MAL2 interacts with IQGAP1 to promote pancreatic cancer progression by increasing ERK1/2 phosphorylation. Biochem Biophys Res Commun. 2021; 554:63–70. https://doi.org/10.1016/j.bbrc.2021.02.146 [PubMed]
  • 38. Goldenring JR. A central role for vesicle trafficking in epithelial neoplasia: intracellular highways to carcinogenesis. Nat Rev Cancer. 2013; 13:813–20. https://doi.org/10.1038/nrc3601 [PubMed]
  • 39. Ji SM. Overexpression of SLC25A15 is involved in the proliferation of cutaneous melanoma and leads to poor prognosis. Med Sci (Paris). 2018; 34:74–80. https://doi.org/10.1051/medsci/201834f113 [PubMed]
  • 40. Kong Z, Wan X, Lu Y, Zhang Y, Huang Y, Xu Y, Liu Y, Zhao P, Xiang X, Li L, Li Y. Circular RNA circFOXO3 promotes prostate cancer progression through sponging miR-29a-3p. J Cell Mol Med. 2020; 24:799–813. https://doi.org/10.1111/jcmm.14791 [PubMed]
  • 41. Martinelli D, Diodato D, Ponzi E, Monné M, Boenzi S, Bertini E, Fiermonte G, Dionisi-Vici C. The hyperornithinemia-hyperammonemia-homocitrullinuria syndrome. Orphanet J Rare Dis. 2015; 10:29. https://doi.org/10.1186/s13023-015-0242-9 [PubMed]
  • 42. Palmieri F. The mitochondrial transporter family SLC25: identification, properties and physiopathology. Mol Aspects Med. 2013; 34:465–84. https://doi.org/10.1016/j.mam.2012.05.005 [PubMed]
  • 43. Li R, Tang X, Xu C, Guo Y, Qi L, Li S, Ren Q, Wu J, Chen D. Circular RNA NF1-419 Inhibits Proliferation and Induces Apoptosis by Regulating Lipid Metabolism in Astroglioma Cells. Recent Pat Anticancer Drug Discov. 2021. [Epub ahead of print]. https://doi.org/10.2174/1574892816666210729125802 [PubMed]
  • 44. Wang YT, Lin MR, Chen WC, Wu WH, Wang FS. Optimization of a modeling platform to predict oncogenes from genome-scale metabolic networks of non-small-cell lung cancers. FEBS Open Bio. 2021; 11:2078–94. https://doi.org/10.1002/2211-5463.13231 [PubMed]
  • 45. Birge RB, Boeltz S, Kumar S, Carlson J, Wanderley J, Calianese D, Barcinski M, Brekken RA, Huang X, Hutchins JT, Freimark B, Empig C, Mercer J, et al. Phosphatidylserine is a global immunosuppressive signal in efferocytosis, infectious disease, and cancer. Cell Death Differ. 2016; 23:962–78. https://doi.org/10.1038/cdd.2016.11 [PubMed]
  • 46. Burstyn-Cohen T, Maimon A. TAM receptors, Phosphatidylserine, inflammation, and Cancer. Cell Commun Signal. 2019; 17:156. https://doi.org/10.1186/s12964-019-0461-0 [PubMed]