TRIB3 Promotes the Malignant Progression of Bladder Cancer: An Integrated Analysis of Bioinformatics and in vitro Experiments

Background: Bladder cancer is a common malignant tumor characterized by high mortality and high management costs; however, it lacks useful molecular prognostic markers. Tribbles pseudokinase 3 (TRIB3) is a pseudokinase that participates in cell tumor progression and metabolism and whose function in bladder cancer is not precisely known.

Main Methods: We downloaded transcriptome data and clinical data of bladder cancer from associated databases and extracted the expression matrix of TRIB3 for multiple bioinformatics analysis. RT-PCR detected the expression of TRIB3 in bladder cancer cells. After knockdown of TRIB3 with siRNA, we investigated TRIB3 function using CCK8, Cell Cycle and Transwell assays.

Key Findings: Kaplan–Meier analysis of TRIB3 in the four cohorts showed that high expression of TRIB3 correlated with poor outcome. Expression of TRIB3 positively correlated with stage and grade and down-regulation of TRIB3 expression significantly inhibited proliferation, migration and cell cycle of bladder cancer cells.

Significance: TRIB3 is a potential prognostic marker and therapeutic target. It can be used to individualize the treatment of bladder cancer.

Introduction

Bladder cancer is common cancer worldwide. It is more common in men than in women, and its prevalence increases sharply with age (Yee et al., 2011). Treatment depends on the stage (Botteman et al., 2003; Steven and Poulsen, 2007), however, treatment responses and overall outcomes in patients at similar stages vary (Ye et al., 2014). The tumor tends to recur and metastasize (Holmäng et al., 1995; Steven and Poulsen, 2007). Therefore, there is an urgent need to identify molecular markers to predict treatment responses and to design individualized treatment (Ye et al., 2014).

Tribbles pseudokinase 3 (TRIB3) is a member of the pseudokinase gene family (Grosshans and Wieschaus, 2000). It lacks kinase activity due to a lack of ATP-specific binding sites and catalytic cores (Manning et al., 2002); however, it can compete with kinases for peptide substrates to regulate biological function (Qi et al., 2006). TRIB3 regulates a series of physiological processes, including cell proliferation, differentiation, apoptosis, glucose metabolism, lipid metabolism, and cell stress (Yokoyama and Nakamura, 2011). TRIB3 protein is also involved in signal pathways such as the MAPK pathway (Kiss-Toth et al., 2004), the PI3K pathway (Ding et al., 2008), the NF-κB pathway (Wu et al., 2003), and TGF-beta (Hua et al., 2011). These findings suggest that TRIB3 participates in various physiological and pathological processes. The expression of TRIB3 is higher in breast cancer (Wennemers et al., 2011), ovarian cancer (Wang et al., 2020), lung cancer (Zhang X. et al., 2019), and colorectal cancer (Hua et al., 2019), and outcomes are poor. However, few studies reported the role of TRIB3 expression in bladder cancer.

In the present study, we used several bladder cancer cohorts and online databases to explore the prognostic potential of TRIB3 and to explore the mechanisms of TRIB3 in bladder cancer cells through functional experiments in vitro.

Materials and Methods

Downloading and Filtering Raw Data

The Cancer Genome Atlas (TCGA¹) is the most authoritative oncology database. We downloaded the transcriptome and sample information, Gene Ontology Consortium (GEO²) and ArrayExpress from TCGA.³ also contain substantial amounts of sequencing data. We obtained transcriptome and sample information from these two databases for GSE32548 (n = 146) (Lindgren et al., 2012), GSE32894 (n = 308) (Sjödahl et al., 2012) and E-MTAB-1803 (n = 85) (Rebouissou et al., 2014). In the R language environment, we standardized all RNA expression matrices and calculated log values using the R package “edge.” Finally, we selected patients with complete transcriptional sequencing and follow-up information. The final numbers of cases included in the analyses were 403 (TCGA-BLCA), 130 (GSE32548), 224 (GSE32894), and 73 (E-MTAB-1803). Specific clinical information is displayed in Table 1.

TABLE 1

Table 1. Summary of basic clinical information of four bladder cancer cohorts.

Kaplan–Meier Analysis

We used the R packages “survival” and “survminer” for survival analysis. We calculated expression values using the function “res. cat.” We used the best cut-off value to divide patients into high-risk and low-risk groups. We performed Kaplan–Meier analysis on the two groups when p < 0.05, using the online tool OSblca⁴ (Zhang G. et al., 2019). The database integrates multiple bladder cancer data sets and is a powerful tool for evaluating bladder cancer prognostic markers. We downloaded transcriptome and survival data containing 195 bladder cancer samples from the attachment to the article for survival validation (Mariathasan et al., 2018).

Univariate and Multivariate Cox Regression Analysis

We used R-package “survival” for univariate and multivariate Cox regression analysis. Multivariate Cox analysis included age, sex, primary tumor extent (T stage), grade and prognostic factor gene.

Receiver Operating Characteristic (ROC) Curve

We plotted the ROC curve of overall survival using the R package “survivalROC,” and calculated the area under the curve (AUC). AUC > 0.5 generally indicates specific prediction ability. Higher AUC suggests more accurate prediction results of the model.

Concordance Index (C-Index)

C-index is used to calculate the discrimination between the predicted and real values of the model in survival analysis. We use R package “survcomp” and “survival” to fit multiple factors for survival analysis and calculate the c-index.

UALCAN Database

UALCAN is a comprehensive and interactive internet resource for analyzing cancer OMICS data⁵ (Chandrashekar et al., 2017). Using this tool, one can query the mRNA distribution of genes in normal and tumor tissues.

Metabase Enrichment Analysis

The internet-based tool cBioPortal queries the co-expression genes of target genes. Spearman coefficients greater than 0.3 indicate that there is co-expression. We analyzed these genes using pathway and process enrichment analysis of the Metabase database⁶ (Zhou et al., 2019). The enrichment analysis included KEGG Pathway, GO Biological Processes, Reactome Gene Sets, Canonical Pathways, and CORUM to evaluate these genes’ potential biological functions and pathways. We selected term enrichment significance according to p-value < 0.01, minimum count of 3, and enrichment factor >1.5.

BioGRID and STRING Database

BioGRID⁷ is a biomedical interaction repository. The database can be used for searching publications for protein and genetic interactions, chemical interactions, and post-translational modifications from significant model organism species. String⁸ database is used for searches of known protein-protein interactions, and it predicts protein-protein interactions.

Human Samples and Ethics Statements

A total of 40 pairs of bladder cancer tissues and para-cancer tissues were collected in the Department of Urology, the First Hospital of China Medical University (Shenyang). The age range of these patients is 38–80 years. The protein extracted from 40 pairs of samples was analyzed by western blot to detect the difference of TRIB3 expression between cancer tissues and para-cancer tissues. The Ethics Committee of the First Hospital of China Medical University approved this study according to the Declaration of Helsinki, and written informed consent was obtained from all the patients.

Cell Culture and Transfection

We cultured T24 and UMUC3 cells in DMEM (Gibco, United States) containing 10% FBS (Gibco) at 37°C and saturated humidity with 5% CO₂. The day before transfection, we seeded exponentially growing cells onto 6-well plates. To increase the transfection efficiency, we added serum-free medium to each well before transfection. We diluted 6 μl Lipofectamine RNAiMAX (Invitrogen, Carlsbad, CA, United States) with 125 μl Opti-MEM (Gibco) and added diluted 3 μl TRIB3-siRNA or TRIB3-nc (RiboBio Co., Ltd., CHINA) with 125 μl Opti-MEM (Gibco). We mixed for 5 min and added them to 6-well plates. We incubated the cells at 37°C in 5% CO₂ was incubated for 12 h. We then replaced the medium with DMEM containing 10% FBS. We performed CCK8 and Transwell assays 24 h after transfection and extracted RNA 48 h after transfection. The TRIB3-siRNA sequences were genOFFTM st-h-TRIB3_001: GGAGTTGGATGACAACTTA and genOFFTM st-h-TRIB3_002: CTACGTGGGACCTGAGATA.

Western Blot Experiment

RIPA buffer (Sigma, St. Louis, MO, United States) was added to the paracancerous-tissue and cancerous-tissue, which was lysed by sonication for 30 min, collecting the tissue lysates at 13,000 rpm for 30 min, at 4°C. And each protein concentration was calculated by the BCA method (Solarbio, Beijing, China). Protein (40 μg) were separated by 10% SDS-PAGE and transferred onto PVDF membranes (Merck Millipore, Billerica, MA, United States), and then blocked by 5% skim milk the membranes were incubated with the dilution (1:1,000) of specific primary antibodies at 4°C overnight. TRIB3 (ImmunoWay Biotechnology, United States, YN1887), β-Actin (Bioss Antibodies, United States) Then, the membranes were incubated with the secondary antibodies for 2 h at room temperature. Goat anti-rabbit IgG (Santa Cruz Biotechnology, United States). the immune-reactive proteins were visualized using Pierce^TM ECL Plus western blotting substrate (Thermo Fisher Scientific).

RT-PCR

We extracted total RNA using TRIzol^® reagent (Thermo Fisher Scientific) according to the manufacturer’s protocols and determined the purity and concentration of RNA solution using a NanoDrop 2000 spectrophotometer (NanoDrop Technologies; Thermo Fisher Scientific). We used the RT-PCR kit (Takara, China) of SuperScript III First-Strand Synthesis System to reverse transcribe RNA into cDNA. We performed RT–PCR using the SYBR premix ExTaq^TM kit (Takara, China). The amplification conditions were as follows: denaturation at 95°C for 30 s, 95°C for 5 s, and 60°C for 30 s, for a total of 40 cycles. β-actin primer: Forward: 5′- CAGCATTGGCAATGAGGGGTTC-3′, Reverse: 5′- AGGTCTTTGCGGATGTCCACGT-3′. TRIB3 primer: Forward: 5′- TGCCCTACAGGCACTGAGTA-3′, Reverse: 5′-GTCCGAGTGAAAAAGGCGTA-3′. We performed quantitative real-time PCR amplification using the LightCycler^TM480 II system (Roche Diagnostics, Switzerland) and calculated results using the 2^–ΔΔCT method.

Cell Proliferation Assay

We prepared exponentially growing cells as single-cell suspensions. We seeded cells in 96-well plates (2,000 cells/well) and assessed them at 0, 24, 48, 72, and 96 h. Time 0 indicated 6 h post-plating. Following the addition of 10 μl of CCK8 solution (Vazyme, Nanjing, China) to each well, we incubated cells for 3 h at 37°C in a 5% CO₂ atmosphere. We incubated the plate with shaking for 5 min. We measured optical density (OD) at 490 nm, using a Microplate reader (Bio-Rad, CA, United States). We performed these experiments in triplicate.

Transwell Assay

We used the Transwell assay to determine the effect of TRIB3 on the migration abilities of T24 and UMUC3 cells. Twenty-four hours after transfection, we collected cells and assayed them. We digested cells and prepared cell suspension. We added 600 μl of medium containing 10% FBS to the lower chamber and added 200 μl of the cell suspension (1 × 10⁵ cells/ml) in serum-free medium to the upper chamber. We incubated the Transwell at 37°C for 24 h. We removed the medium from the upper chamber and fixed the cells with 4% paraformaldehyde for 30 min. Then, we washed the cells three times with PBS and stained them with 0.2% crystal violet for 30 min. Finally, we photographed and counted the cells in five independent microscopic fields. We performed these experiments in triplicate.

Cell Cycle Assay

The bladder cancer cells transfected by TRIB3-NC and TRIB3-siRNA for 24 h were collected, and digested with trypsin, and washed three times with PBS, and centrifuged at 2000 g for 5 min, then removing the supernatant. Adjust the cell concentration to 2 × 106 cells/ml. A 700 μl of 100% chilled ethanol and 300 μl pre-cool PBS was added to the cells, which was fixed for 2 h at −20°C. And centrifuged at 2000 g for 5 min and removed the ethanol. 500 μl of the PI/RNase (Solarbio, Beijing, China) staining working solution (PI:RNase A was prepared at 4:1) was added. Following staining at 37°C for 30 min, flow cytometry was performed.

Results

TRIB3 Predicts Overall Survival in Bladder Cancer Cohorts

In four cohorts, we divided samples into two groups according to the expression of TRIB3 for Kaplan–Meier analysis. There were significant differences in E-MTAB-1803 (p = 0.048) (Figure 1A), TCGA_BLCA (p < 0.001) (Figure 1B), GSE32548 (p = 0.004) (Figure 1C) and GSE32894 (p = 0.001) (Figure 1D). The outcomes of the high-expression groups in all four cohorts were worse than those of the low-expression group. We verified the ability of TRIB3 to predict overall survival using OSblca. The tool integrated the data of TCGA-BLCA, GSE13507, GSE31684, GSE48075, and GSE48276 cohorts (n = 934) to draw Kaplan-Meier overall survival curves for total bladder cancer (Figure 1E) and muscle-invasive bladder cancer (Figure 1F). Outcomes in the high expression TRIB3 group were poor (p < 0.05).

FIGURE 1

Figure 1. Kaplan–Meier analysis of TRIB3 expression. (A) E-MTAB-1803. (B) TCGA_BLCA. (C) GSE32548. (D) GSE32894. (E) Kaplan–Meier curve of TRIB3 expression in a meta-cohort of bladder cancer. (F) Kaplan–Meier curve of TRIB3 expression in a meta-cohort of muscle-invasive bladder cancer. Red represents the TRIB3 high expression group, green represents the TRIB3 low expression group, and CI represents the confidence interval.

ROC Curve and AUC Value of TRIB3 on Overall Survival Rate

Receiver operating characteristic curve and AUC values return diagnostic values of markers. We drew 1-, 3-, and 5-year ROC curves of TRIB3 expression on the overall survival rate in four cohorts and calculated the AUC corresponding to each curve. The AUCs of the 1-, 3-, and 5-year TCGA-BLCA cohorts were 0.664, 0.592, and 0.550, respectively (Figure 2A). The AUCs of the 1-, 3-, and 5-year GSE32894 cohorts were 0.578, 0.712, and 0.720, respectively (Figure 2B). The AUC values of the 1-, 3-, and 5-year GSE32548 cohorts were 0.752, 0.704, and 0.665, respectively (Figure 2C). The AUC values of the 1-, 3-, and 5-year E-MTAB-1803 cohorts were 0.620, 0.447, and 0.523 (Figure 2D).

FIGURE 2

Figure 2. ROC curve. (A) TCGA-BLCA. (B) GSE32894. (C) GSE32548. (D) E-MTAB-1803. The blue, green, and red lines represent the ROC curve of 1-year, 3-year, and 5-year overall survival rates.

Univariate and Multivariate COX Analysis

We performed univariate and multivariate analysis of clinical factors (Age, sex, T stage, grade) and TRIB3 (Table 2). Univariate COX analysis showed that TRIB3 and T-stage had prognostic ability in all four cohorts. Multivariate COX analysis showed that T stage still had independent prognostic ability in four cohorts, while TRIB3 had independent prognostic ability in cohort TCGA-BLCA, GSE32548, and GSE32894. To sum up, the prognostic ability of TRIB3 is stronger than that of age, sex, and grade. Than we calculated the C-index of TRIB3, T stage, and integrative model (Combination of TRIB3 and T stage). The C-index of them are all more 0.6 in the four cohort. However, the C-index of T stage is more than the C-index of TRIB3. In the TCGA-BLCA and GSE32894 cohorts, the C-index of integrative model is more than T stage, and it’s the opposite in the GSE32548 and E-MTAB-1803 cohorts. Therefore, whether we can combine TRIB3 and T stage for survival prediction needs a larger cohort to explore (Supplementary Table 1).

TABLE 2

Table 2. Univariate and multivariate Cox regression analysis of clinical-factors/TRIB3 with overall survival rate in patients.

TRIB3 Is Positively Correlated With the Grade of Bladder Cancer

We compared the difference of TRIB3 expression in various grades and T stages. We found that expression levels of TRIB3 differed according to grades and T stage of TCGA-BLCA, GSE32894, GSE32548, and E-MTAB-1803 (p < 0.05, Figures 3A,B). High grade correlated with a trend toward increased expression levels of TRIB3. In cohort GSE32548 and GSE32894, the expression of TRIB3 was different only in different grades and T stages (p < 0.05, Figures 3A,B).

FIGURE 3

Figure 3. The expression of TRIB3 in various clinical types. (A) The expression scatters the plot of TRIB3 under different grades in four cohorts. (B) The expression scatters the plot of TRIB3 under different T stages in four cohorts.

There Is Differential Expression of TRIB3 in Various Types of Bladder Cancer

We obtained data on bladder cancer typing from UCSC Xena⁹. The expression of TRIB3 in various immune subtypes also differed (p < 0.05) (Figure 4A). Among them, the expression level of TRIB3 of basal squamous subtype is the highest. The expression of TRIB3 differed depending on the mRNA cluster (p < 0.001, Figure 4B). Among them, the expression level of TRIB3 of C2 (IFN-gamma dominant) is the highest. Using web tool “UALCAN” to query the data of TRIB3 expression in bladder cancer tissues and normal tissues, we found that TRIB3 was highly expressed in tumor tissues (p < 0.001, Figure 4C).

FIGURE 4

Figure 4. The differential expression of TRIB3 in other typing schemes. (A) Expression levels of the five immunophenotypes: C1 (wound healing), C2 (IFN-gamma dominant), C3 (inflammatory), C4 (lymphocyte depleted), and C6 (TGF-b dominant). (B) Expression levels in five mRNA types: basal squamous, luminal, luminal infiltrated, luminal papillary, and neuronal subtype. (C) Comparison of expression levels of TRIB3 in normal tissue and tumor tissue. p < 0.001 indicates statistical significance.

TRIB3 Expression Level in Bladder Cancer Tissues and Its Clinical Significance

To confirm the above results, we performed western blot experiment on tissues samples obtained from 40 bladder cancer patients operated at the First Hospital of China Medical University. The results show that the expression of TRIB3 in cancer tissues was significantly higher than the paracancer tissues (Figure 5A), and confirmed by Image J (NIH, United States) which was used to densitometer quantitative analysis of TRIB3 Western blot results of 40 pairs of tissues (Figure 5B). The analysis results showed that TRIB3 was relatively highly expressed in 35 of the 40 pairs of bladder cancer tissue samples (Figures 5A,B). Moreover, to assess the clinical significance of TRIB3 expressions in bladder cancer, the clinical characteristics between TRIB3 and bladder cancer were evaluated in tissues. It was revealed that high expression levels of TRIB3 were significantly associated with higher tumor TNM stage (Table 3). However, there are no significant statistical significance between the expression of trib3 and whether distant metastasis occurs. We believe it is due to the insufficient number of bladder cancer patients with metastasis. Taken together, our observations strongly suggest that TRIB3 might have a pro-oncogenic role in bladder cancer. Finally, we used Sanjeev’s bladder cancer cohort for verification, and it was found that TRIB3 still had the ability of prognosis (p < 0.05, Figure 5C).

FIGURE 5

Figure 5. High TRIB3 expression in bladder cancer tissues of patients. (A) The expression of TRIB3 protein in bladder cancer tissues and paracancer tissues by western blot (B) The relative gray values of TRIB3 protein in bladder cancer tissues. (C) Kaplan–Meier analysis of TRIB3 expression in validation cohort.

TABLE 3

Table 3. The clinical significance of TRIB3 expressions in bladder cancer.

Down-Regulation of TRIB3 Induces Cell Cycle Arrest in G0/G1 Phase

We searched for genes that positively correlate to TRIB3 expression in the BLCA (TCGA 2017) cohort of the cBioPortal website (Spearman’s Correlation >2017). We introduced these genes into the Metascape database for enrichment analysis. The three items with significant enrichment were Cell Cycle, mitotic nuclear division, and Cytosolic tRNA aminoacylation (Figure 6A), all of which relate to the cell cycle. Then flow cytometry verified that the proportion of the bladder cancer cells transfected by TRIB3-siRNA in the G0/G1 phase was significantly increased (p < 0.05), as compared to the cells transfected by TRIB3-NC (Figure 7D).

FIGURE 6

Figure 6. Metascape enrichment analysis and protein interactions. (A) The Metascape enrichment analysis results show that larger log10 (P) correlates with the significance of the enrichment. (B) Protein-protein interactions in the STRING database. (C) Protein-protein interactions in the BioGRID database.

FIGURE 7

Figure 7. TRIB3-siRNA inhibits the proliferation and migration and induces cell arrest in bladder cancer cells. (A) Relative mRNA levels of control and si-TRIB3 in T24 and UMUC3 cells. (B) The CCK-8 assay results showed that T24 and UMUC3 cells’ proliferation ability decreased when treated with si-TRIB3. (C) The invasion assay shows that the invasive ability of T24 and UMUC3 cells treated with si-TRIB3 decreased. (D) G0/G1 phase arrest of bladder cancer after TRIB3 siRNA transfection, confirmed by statistical analysis. *P < 0.05 and **P < 0.01.

Proteins Interacting With TRIB3

We found that TRIB3 and AKT1/AKT2 are the most closely related (Figure 6B). The BioGRID database showed Physical Edges between TRIB3 and a large number of proteins (Figure 6C). The most significant correlations were AES, SQSTM1, MDF1, SIRT1, RELA, SMAD3, PML, KAT5, and AKT1.

Further Study of TRIB3 in a Bladder Cancer Cell Knockdown Test

Because of the high expression of TRIB3 in bladder cancer and poor outcome, we knocked down TRIB3 in bladder cancer cells and conducted functional experiments further to explore the function of TRIB3 (Figure 7A). CCK8 analysis showed that the proliferation ability of si-TRIB3 cells was significantly lower (p < 0.05) (Figure 7B). The Transwell analysis showed that si-TRIB3 bladder cancer cells’ migration ability was significantly lower (p < 0.05) (Figure 7C).

Discussion

Bladder cancer is one of the most common malignant tumors worldwide. Because of the lack of specific characteristics in the early stages, most patients present with advanced bladder cancer. Bladder cancer treatment includes surgery, chemotherapy, radiotherapy, and poor outcomes (Marta et al., 2012; Racioppi et al., 2012). Mining robust molecular markers and therapeutic targets stratify patients at risk to facilitate treatment decision-making and individualized treatment. In the present study, we combined computational biology methods and experimental verification methods to identify and characterize the role of pseudokinase-TRIB3 in poor outcomes and malignant progression of bladder cancer. We demonstrated the utility of integrating bioinformatics and in vitro experiments to explore potential diagnostic and therapeutic targets for bladder cancer.

One of the main findings of this study was that TRIB3 is a robust prognostic biomarker of bladder cancer. We obtained datasets from TCGA-BLCA, GSE32548, GSE32894, and E-MTAB-18034. Kaplan–Meier analysis of all four cohorts showed that high expression levels of TRIB3 correlated with shorter survival rates. The AUC values of 1BI 3 and 5 years in the cohort were almost all greater than 0.5, suggesting that TRIB3 expression levels can predict risk at these time points. TRIB3 successfully predicted outcome in patients with bladder cancer in two meta-cohorts. These results suggest that expression levels of TRIB3 may serve as an index to predict outcomes for patients with bladder cancer. The four queues used in this study derive from GEO, TCGA, and ArrayExpress databases. These three databases are currently the largest three transcriptome databases. Data upload has strict requirements and high reliability. The four queues’ samples derive from several regions, which suggests that TRIB3 expression a wide range of generalizability. Our study’s scale was polycentric, cross-regional, and large-scale, which increases the reliability of our results.

Another important finding was that TRIB3 is a potential oncogene in bladder cancer. We identified TRIB3 expression in all four cohorts. Expression levels of TRIB3 increased with increased stage and grade, suggesting that TRIB3 promotes the progression of bladder cancer and increases the risk of invasion and deterioration of bladder cancer. We confirmed high expression levels of TRIB3 mRNA in bladder cancer cells using qRT-PCR. After inhibiting the expression of TRIB3 in bladder cancer cells, we found that the proliferation and migration of bladder cancer cells decreased. This experiment suggests that expression of TRIB3 plays a tumor-promoting role in bladder cancer. Studies reported that TRIB3 promotes acute promyelocytic leukemia by stabilizing oncoprotein PML-RAR α and inhibiting p53-mediated aging (Li K. et al., 2017). Knockout of TRIB3 inhibits tumorigenesis and cancer progression. It disrupts the TRIB3-SQSTM1 or TRIB3-PML/RAR α interaction through specific peptides, and this is a potential strategy for treating certain solid cancers and acute promyelocytic leukemia (Hua et al., 2011; Li K. et al., 2017). Studies confirmed that inhibition of the expression of TRIB3 in ovarian cancer cells causes cell cycle arrest in the G0/G1 phase in ovarian cancer cells. Proliferation, migration, and invasion of ovarian cancer cells reduce, and apoptosis increases (Wang et al., 2020). Enrichment analysis of co-expressed genes of TRIB3 showed that TRIB3 was related to the cell cycle (Figure 6A). Flow cytometry analysis revealed that arrested G0/G1 cell cycle could be induced by TRIB3-siRNA, confirmed by statistical analysis (Figure 7D), thereby inhibiting the proliferation and migration of bladder cancer cells. Other studies confirmed that the expression of TRIB3 is greater in solid tumors such as colorectal cancer, breast cancer, lung cancer, ovarian cancer, melanoma, liver cancer, and leukemia (Wennemers et al., 2011; Salazar et al., 2015; Hua et al., 2019; Zhang X. et al., 2019; Wang et al., 2020; Zhu et al., 2020), suggesting that TRIB3 is a potential target for cancer treatment. Some studies have shown that biological scaffold coordination complexes can be used to design new target-specific therapeutic agents, and TRIB3 is also a scaffold protein that acts as a platform response. In the future, we can try to combine TRIB3 with new compounds for combined therapy (Mirgayazova et al., 2019). We conclude from these results that TRIB3 is a potential oncogene and a potential therapeutic target for bladder cancer.

Conclusion

Overall, TRIB3 successfully predicted outcomes in patients with bladder cancer in four cohorts. It may serve as a prognostic marker for individualized treatment of bladder cancer in the future. TRIB3 expression promotes the proliferation and invasion of bladder cancer cells. Although more clinical samples and more in-depth experiments are needed to explore the mechanism of TRIB3 in the future, our results suggest that TRIB3 is a potential prognostic marker and therapeutic target in bladder cancer.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

Author Contributions

JY and JL contributed to conception and design. JY, JL, and MY carried out analysis and wrote the manuscript. JA, YKZ, and SJ collected and processed the data. YYZ and YY edited the manuscript and provided constructive comments. All authors read and approved the final manuscript.

Funding

This work was supported in part by the National Natural Science Foundation of China (grants 81672523, 81472404, 81472403, 81272834, and 31000572), the 2018 Support Plan for innovative talents in Colleges and Universities of Liaoning Province (China), the 2018 “million talents Project” funded Project of Liaoning Province (China), the 2019 Key Research and Development Projects of Shenyang (China), and Start-up Fund for Returnee of Chinese Scholars by Shenyang Medical College (grant 20204039).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We appreciate the free use of TCGA, GEO, ArrayExpress databases.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.649208/full#supplementary-material

Supplementary Table 1 | C-index of the TRIB3, T stage and Integrative model.

Abbreviations

AUC, area under the curve; GEO, Gene Expression Omnibus; ROC, receiver operating characteristic; TCGA, The Cancer Genome Atlas.

Footnotes

1. ^ https://cancergenome.nih.gov/
2. ^ http://www.ncbi.nlm.nih.gov/geo/
3. ^ https://www.ebi.ac.uk/arrayexpress/
4. ^ http://bioinfo.henu.edu.cn/BLCA/BLCAList.jsp
5. ^ http://ualcan.path.uab.edu/index.html
6. ^ http://metascape.org
7. ^ https://thebiogrid.org/
8. ^ https://string-db.org/
9. ^ http://xena.ucsc.edu/

References

Botteman, M. F., Pashos, C. L., Redaelli, A., Laskin, B., and Hauser, R. (2003). The health economics of bladder cancer: a comprehensive review of the published literature. PharmacoEconomics 18, 1315–1330. doi: 10.1007/bf03262330

PubMed Abstract | CrossRef Full Text | Google Scholar

Chandrashekar, D. S., Bashel, B., Balasubramanya, S. A. H., Creighton, C. J., Ponce-Rodriguez, I., Chakravarthi, B. V. S. K., et al. (2017). Ualcan: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia 8, 649–658. doi: 10.1016/j.neo.2017.05.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, J., Kato, S., and Du, K. (2008). Pi3k activates negative and positive signals to regulate trb3 expression in hepatic cells. Exp. Cell Res. 7, 1566–1574. doi: 10.1016/j.yexcr.2008.01.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Grosshans, J., and Wieschaus, E. (2000). A genetic link between morphogenesis and cell division during formation of the ventral furrow in drosophila. Cell 5, 523–531. doi: 10.1016/s0092-8674(00)80862-4

CrossRef Full Text | Google Scholar

Holmäng, S., Hedelin, H., Anderstr, C., and Johansson, S. L. (1995). The relationship among multiple recurrences, progression and prognosis of patients with stages ta and t1 transitional cell cancer of the bladder followed for at least 20 years. J. Urol. 6, 1823–1826; discussion1826–7.

Google Scholar

Hua, F., Mu, R., Liu, J., Xue, J., Wang, Z., Lin, H., et al. (2011). Trb3 interacts with smad3 promoting tumor cell migration and invasion. J. Cell Sci. 124, 3235–3246. doi: 10.1242/jcs.082875

PubMed Abstract | CrossRef Full Text | Google Scholar

Hua, F., Shang, S., Yang, Y. W., Zhang, H. Z., Xu, T. L., Yu, J. J., et al. (2019). Trib3 interacts with β-catenin and tcf4 to increase stem cell features of colorectal cancer stem cells and tumorigenesis. Gastroenterology 3, 708–721.e15.

Google Scholar

Kiss-Toth, E., Bagstaff, S. M., Sung, H. Y., Jozsa, V., Dempsey, C., Caunt, J. C., et al. (2004). Human tribbles, a protein family controlling mitogen-activated protein kinase cascades. J. Biol. Chem. 41, 42703–42708. doi: 10.1074/jbc.m407732200

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, K., Wang, F., Cao, W. B., Lv, X. X., Hua, F., Cui, B., et al. (2017). Trib3 promotes apl progression through stabilization of the oncoprotein pml-rarα and inhibition of p53-mediated senescence. Cancer Cell 31, 697–710.e7.

Google Scholar

Lindgren, D., Sjödahl, G., Lauss, M., Staaf, J., Chebil, G., Lövgren, K., et al. (2012). Integrated genomic and gene expression profiling identifies two major genomic circuits in urothelial carcinoma. PLoS One 6:e38863. doi: 10.1371/journal.pone.0038863

PubMed Abstract | CrossRef Full Text | Google Scholar

Manning, G., Whyte, D. B., Martinez, R., Hunter, T., and Sudarsanam, S. (2002). The protein kinase complement of the human genome. Science 5600, 1912–1934. doi: 10.1126/science.1075762

PubMed Abstract | CrossRef Full Text | Google Scholar

Mariathasan, S., Turley, S. J., Nickles, D., Castiglioni, A., Yuen, K., Wang, Y., et al. (2018). Tgfβ attenuates tumour response to pd-l1 blockade by contributing to exclusion of t cells. Nature 7693, 544–548.

Google Scholar

Marta, G. N., Hanna, S. A., Gadia, R., Correa, S. F., Silva, J. L., and Carvalho Hde, A. (2012). The role of radiotherapy in urinary bladder cancer: current status. Int. Braz. J. Urol. 2, 144–153; discussion153–4. doi: 10.1590/s1677-55382012000200002

PubMed Abstract | CrossRef Full Text | Google Scholar

Mirgayazova, R., Khadiullina, R., Mingaleeva, R., Chasov, V., Gomzikova, M., Garanina, E., et al. (2019). Novel isatin-based activator of p53 transcriptional functions in tumor cells. Mol. Biol. Res. Commun. 8, 119–128.

Google Scholar

Qi, L., Heredia, J. E., Altarejos, J. Y., Screaton, R., Goebel, N., Niessen, S., et al. (2006). Trb3 links the e3 ubiquitin ligase cop1 to lipid metabolism. Science 5781, 1763–1766. doi: 10.1126/science.1123374

PubMed Abstract | CrossRef Full Text | Google Scholar

Racioppi, M., D’Agostino, D., Totaro, A., Pinto, F., Sacco, E., D’Addessi, A., et al. (2012). Value of current chemotherapy and surgery in advanced and metastatic bladder cancer. Urol. Int. 3, 249–258. doi: 10.1159/000335556

PubMed Abstract | CrossRef Full Text | Google Scholar

Rebouissou, S., Bernard-Pierrot, I., de Reyniès, A., Lepage, M. L., Krucker, C., Chapeaublanc, E., et al. (2014). Egfr as a potential therapeutic target for a subset of muscle-invasive bladder cancers presenting a basal-like phenotype. Sci. Transl. Med. 6:244ra91. doi: 10.1126/scitranslmed.3008970

PubMed Abstract | CrossRef Full Text | Google Scholar

Salazar, M., Lorente, M., García-Taboada, E., Pérez Gómez, E., Dávila, D., Zúñiga-García, P., et al. (2015). Loss of tribbles pseudokinase-3 promotes akt-driven tumorigenesis via foxo inactivation. Cell Death Differ. 22, 131–44. doi: 10.1038/cdd.2014.133

PubMed Abstract | CrossRef Full Text | Google Scholar

Sjödahl, G., Lauss, M., Lövgren, K., Chebil, G., Gudjonsson, S., Veerla, S., et al. (2012). A molecular taxonomy for urothelial carcinoma. Clin. Cancer Res. 12, 3377–3386. doi: 10.1158/1078-0432.ccr-12-0077-t

PubMed Abstract | CrossRef Full Text | Google Scholar

Steven, K., and Poulsen, A. L. (2007). Radical cystectomy and extended pelvic lymphadenectomy: survival of patients with lymph node metastasis above the bifurcation of the common iliac vessels treated with surgery only. J. Urol. 178, 1218–1223. doi: 10.1016/j.juro.2007.05.160

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., Wang, C., Li, X., Hu, Y., Gou, R., Guo, Q., et al. (2020). Down-regulation of trib3 inhibits the progression of ovarian cancer via mek/erk signaling pathway. Cancer Cell Int. 20:418.

Google Scholar

Wennemers, M., Bussink, J., Scheijen, B., Nagtegaal, I. D., van Laarhoven, H. W., Raleigh, J. A., et al. (2011). Tribbles homolog 3 denotes a poor prognosis in breast cancer and is involved in hypoxia response. Breast Cancer Res. 13:R82.

Google Scholar

Wu, M., Xu, L. G., Zhai, Z., and Shu, H. B. (2003). Sink is a p65-interacting negative regulator of nf-kappab-dependent transcription. J. Biol. Chem. 29, 27072–27079. doi: 10.1074/jbc.m209814200

PubMed Abstract | CrossRef Full Text | Google Scholar

Ye, F., Wang, L., Castillo-Martin, M., McBride, R., Galsky, M. D., Zhu, J., et al. (2014). Biomarkers for bladder cancer management: present and future. Am. J. Clin. Exp. Urol. 1, 1–14. doi: 10.3109/1354750x.2012.717105

PubMed Abstract | CrossRef Full Text | Google Scholar

Yee, D. S., Ishill, N. M., Lowrance, W. T., Herr, H. W., and Elkin, E. B. (2011). Ethnic differences in bladder cancer survival. Urology 78, 544–549. doi: 10.1016/j.urology.2011.02.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Yokoyama, T., and Nakamura, T. (2011). Tribbles in disease: signaling pathways important for cellular function and neoplastic transformation. Cancer Sci. 102, 1115–1122. doi: 10.1111/j.1349-7006.2011.01914.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, G., Wang, Q., Yang, M., Yuan, Q., Dang, Y., Sun, X., et al. (2019). Osblca: a web server for investigating prognostic biomarkers of bladder cancer patients. Front. Oncol. 9:466. doi: 10.3389/fonc.2019.00466

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X., Zhong, N., Li, X., and Chen, M. B. (2019). Trib3 promotes lung cancer progression by activating β-catenin signaling. Eur. J. Pharmacol. 863:172697. doi: 10.1016/j.ejphar.2019.172697

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Y., Zhou, B., Pache, L., Chang, M., Khodabakhshi, A. H., Tanaseichuk, O., et al. (2019). Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat. commun. 10:1523. doi: 10.1038/s41467-019-09234-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, Y., Cui, J., Liu, J., Hua, W., Wei, W., and Sun, G. (2020). Sp2 promotes invasion and metastasis of hepatocellular carcinoma by targeting trib3 protein. Cancer Med. 9, 3592–3603. doi: 10.1002/cam4.2977

PubMed Abstract | CrossRef Full Text | Google Scholar