Silencing TAK1 alters gene expression signatures in bladder cancer cells
- Authors:
- Published online on: March 7, 2017 https://doi.org/10.3892/ol.2017.5819
- Pages: 2975-2981
-
Copyright: © Chen et al. This is an open access article distributed under the terms of Creative Commons Attribution License.
Abstract
Introduction
Bladder cancer may be classified into two stages: Non-muscle-invasive disease of low grade, and progressive muscle-invasive disease of high grade, which may further deteriorate into metastatic cancer (1). This type of cancer predominantly occurs as urothelial cell carcinoma and ranks as the fourth most common malignancy and the eighth most common cause of cancer-associated mortality in men (2). Overexpressed or mutant fibroblast growth factor receptor 3 (FGFR3) may contribute to the development of the transformed phenotype of urothelial carcinoma (3), and the targeted inhibition of FGFR3 may thus prevent the development of superficial bladder cancer (4). The transcription factor nuclear factor κB (NFκB) is an important mediator of the angiogenesis and metastasis of bladder cancer (5). Activating mutations in FGFR3 can promote NFκB transcriptional activity through the important mediator transforming growth factor-β-activated kinase 1 (TAK1), which exerts a positive regulatory effect on the activity of NFκB in urothelial cell carcinoma (6). However, the exact regulatory effect of TAK1 on the progression of bladder cancer remains unclear.
RNA interference may effectively aid in defining the roles of specific genes in the progression of cancer when specifically designed siRNAs are used to silence target genes (7,8). Using the microarray expression data of bladder cancer cells transfected with control siRNA or TAK1-specific siRNA which were deposited by Salazar et al (6), the present study aimed to identify the differentially expressed genes (DEGs) and to screen for tumor-associated DEGs, followed by identification of the biological processes or pathways implicated by DEGs and the hub nodes in the complex protein-protein interaction (PPI) network and sub-network, in an attempt to provide a deeper insight into the molecular mechanisms mediating the effect of silencing TAK1 on bladder cancer.
Materials and methods
Gene expression profiles
The gene expression profiles (accession number GSE52452) (6) from 6 samples of MGHU3 (Y375C) mutant FGFR3 bladder cancer cells that were transfected with control siRNA (GSM1267159-61) or TAK1-specific siRNA (GSM1267150-52), were downloaded from the public functional genomics data repository of Gene Expression Omnibus (9). These 6 samples (3 control samples vs. 3 TAK1 siRNA-treated samples) were previously researched using Affymetrix Human Gene 1.0 ST Array.
Data processing and screening of DEGs
The raw expression data were processed through RMA background correction (10), quantile normalization, log base 2 (log2) transformation and probeset summarization to obtain the gene expression matrix by employing Affy package in Bioconductor (11) and probe annotation files provided by Brain Array Lab (brainarray.mbni.med.umich.edu/). Subsequently, the lists of DEGs between control and TAK1-specific siRNA-treated cells were generated under the thresholds of |log2 fold change (FC)|≥1 (fold-change magnitude, >2) and P-value <0.05 using the empirical Bayes method offered by Limma package (12).
Functional enrichment analysis of DEGs
To provide insight into the biological functions or pathways involving the identified DEGs, the Gene Ontology (GO) Biological Process (13) and Kyoto Encyclopedia of Genes and Genomes (KEGG) (14) pathway enrichment analyses of DEGs were performed with the criterion of P<0.05 using the Database for Annotation, Visualization and Integrated Discovery (15).
Functional annotations of DEGs
To identify the specific functions of DEGs, transcription factors were screened from the DEGs based on the TRANSFAC database (16), accompanied by screening of tumor-suppressor genes (TSGs) using the TSGene (17) database and identification of oncogenic genes using the tumor-associated genes (TAG) database (18).
Construction of PPI network and screening of sub-network
In consideration of the gene-gene interactions in the complex biological systems, the identified DEGs were inputted into the Search Tool for the Retrieval of Interacting Genes (STRING) database (19) to identify the interacting pairs with a combined score >0.9 (data downloaded on May 9, 2014). The PPI network was constructed using Cytoscape software (20), followed by excavation of the sub-network using the BioNet tool (bionet.bioapps.biozentrum.uni-wuerzburg.de/). KEGG pathway enrichment analysis of DEGs in the sub-network was performed to identify the associated biological pathways.
Results
Screened DEGs
With the criteria of |log2 FC|≥1 and P<0.05, a total of 43 downregulated genes and 21 upregulated genes were identified in the mutant FGFR3 bladder cancer cells transfected with TAK1-specific siRNA in comparison with the cells transfected with control siRNA.
Functional enrichment results of DEGs
By performing a KEGG pathway enrichment analysis of the separate upregulated and downregulated genes, it was demonstrated that the downregulated genes were significantly enriched in five pathways, including the NOD-like receptor signaling pathway [P=0.014731; mitogen-activated protein kinase kinase kinase 7 (MAP3K7) and tumor necrosis factor α-induced protein 3 (TNFAIP3)] and the upregulated genes were enriched in the cytokine-cytokine receptor interaction pathway [P=0.018547; epidermal growth factor receptor (EGFR) and interleukin 18 (IL18)] (Table I). The GO Biological Process enrichment analysis revealed that the downregulated genes were associated with cellular responses, including regulation of catalytic activity [P=5.45×10−4; e.g. aldehyde dehydrogenase 1 family member A1 (ALDH1A1), β-arrestin (ARRB1), inositol trisphosphate receptor (ITPR1) and MAP3K7)] and cellular response to chemical stimulus [P=1.38×10−3; e.g. ALDH1A1, ARRB1, ITPR1, MAP3K7, cytochrome P450 family 3 subfamily A polypeptide 5 (CYP3A5) and glutathione S-transferase α1 (GSTA1)]; the upregulated genes were associated with cellular developmental processes [P=5.80×10−4; e.g. EGFR, early growth response 1 (EGR1) and parathyroid hormone-like hormone (PTHLH)] (Table II).
Results of functional annotation of DEGs
Using TRANSFAC, TSGene and TAG databases, the functional annotation analysis screened 1 transcription factor (upregulated EGR1) and 9 TAGs, including 1 oncogene (upregulated EGFR), 5 TSGs (downregulated TNFAIP3, ras association domain family member 5 and carcinoembryonic antigen related cell adhesion molecule 7; upregulated GLI pathogenesis related 1 and EGR1) and 3 other TAGs (downregulated breast carcinoma amplified sequence 1 and activin membrane-bound inhibitor; upregulated PTHLH).
Analysis of the PPI network
Using the STRING database and Cytoscape tool, a PPI analysis of the DEGs was performed and an integrated PPI network was obtained. Based on this network, the top 10 hub nodes representing important roles in the complex metabolic networks were EGFR (degree, 178), CYP3A5 (degree, 83), MAP3K7 (degree, 58), GSTA1 (degree, 56), PTHLH (degree, 47), ALDH1A1 (degree, 45), potassium voltage-gated channel subfamily D member 2 (degree, 43), EGR1 (degree, 42), ARRB1 (degree, 36) and ITPR1 (degree, 33) (Fig. 1). Additionally, this network revealed a correlation of MAP3K7 with TNFAIP3.
Analysis of the screened sub-network
Using the BioNet tool, a sub-network with EGFR at the core was screened from the PPI network, in which EGFR was associated with v-erb-b2 avian erythroblastic leukemia viral oncogene homolog 2 (ERBB2), growth factor receptor-bound protein 2 (GRB2) and phosphoinositide-3-kinase, regulatory subunit 1 (PIK3R1) (Fig. 2). Furthermore, the DEGs in the sub-network were found to be enriched in various pathways associated with cancer, including the ErbB signaling pathway (P=3.07×10−7; e.g. EGFR, ERBB2, GRB2 and PIK3R1), pathways in cancer (P=5.90×10−7; e.g. EGFR, ERBB2, GRB2 and PIK3R1), prostate cancer (P=1.04×10−4; e.g. EGFR, ERBB2, GRB2 and PIK3R1), endometrial cancer (P=1.64×10−4; EGFR, ERBB2, GRB2 and PIK3R1), non-small cell lung cancer (P=1.90×10−4; EGFR, ERBB2, GRB2 and PIK3R1) and glioma (P=3.90×10−4; e.g. EGFR, GRB2 and PIK3R1) (Table III).
Table III.Top 10 significantly enriched KEGG pathways of differentially expressed genes in the sub-network. |
Discussion
By re-analyzing the microarray data from specific mutant FGFR3 bladder cancer cells using bioinformatic methods, the present study identified 64 genes with significantly altered expression between TAK1-specific siRNA-treated cells and the control siRNA-treated cells, including 9 TAGs. The downregulated genes were related to the NOD-like receptor signaling pathway and cellular response functions, whereas the upregulated genes were associated with cellular developmental processes. By constructing a PPI network, the present study identified 10 hub nodes that may exert the predominant effects on the network. Furthermore, the DEGs in the sub-network with EGFR at the core were associated with various types of cancer and the ErbB signaling pathway.
TAK1 siRNA vs. control siRNA samples yielded 43 downregulated genes, of which 6 were identified to be hub nodes in the PPI network, including CYP3A5, MAP3K7, GSTA1, ALDH1A1, ARRB1 and ITPR1. CYP3A5 and GSTA1 are associated with the detoxification of chemical stimuli, electrophilic compounds or other damage stimuli (21,22). They were also enriched in the term metabolism of xenobiotics by cytochrome P450, which plays important roles in detoxification (23). MAP3K7, which correlated with TNFAIP3, was enriched in the NOD-like receptor signaling pathway which is involved in sensing intracellular microbial motifs or other damage stimuli (24). ALDH1A1 and CYP3A5 were enriched in retinol metabolism; the loss of retinol acyltransferase is inversely correlated with the invasion bladder cancer (25). ARRB1 is considered to cause a specific dampening cellular response to stimuli or sensory signals (26) and ITPR1 mediates calcium release which then amplifies apoptosis in response to specific stimuli (27,28). Accordingly, the present study further demonstrated that the 6 hub nodes in the PPI network were enriched in the GO function cellular response to chemical stimulus. It is therefore possible that silencing of TAK1 using specific siRNA transfection may alter the cellular response-associated pathways or functions in response to stimuli via regulating the expression of the 6 DEGs.
TAK1 siRNA treatment induced 21 upregulated genes in bladder cancer cells, among which EGFR, PTHLH and EGR1 were identified as hub nodes in the PPI network. These 3 genes were all enriched in cellular developmental processes. As exemplified, activated EGFR, which contributes to phenotypic characteristics in various tumor types (29), is an effective therapeutic target for the treatment of bladder cancer (30,31). Increased expression of PTHLH resulting from downregulated p38MAPK signaling is associated with metastatic lesions to the liver and lung from colon cancer cells (32). EGR1 also exerts pro-tumorigenic effects by contributing to tumor infiltration, node formation and metastasis (33). Thus, the upregulation of these genes in TAK1 siRNA-treated bladder cancer cells may suggest that silencing TAK1 has the effect of promoting cancer cell developmental processes. Furthermore, the present study also predicted EGFR, EGR1 and PTHLH as TAGs, which suggests their potential use as therapeutic targets for the diagnosis and treatment of bladder cancer.
Using BioNet software, a sub-network with EGFR at the core was screened from the PPI network. According to the PPI network, EGFR was correlated with ERBB2, GRB2 and PIK3R1. These 4 DEGs were enriched in various cancer types, including prostate cancer, endometrial cancer, non-small cell lung cancer and glioma. As previously reported, the EGFR family of four receptors including EGFR and ERBB2, is implicated in the development and progression of various human cancer types (34). The activation of ERBB2 may result in resistance to cetuximab-based therapy targeting EGFR, while the inhibition of this gene can restore cetuximab sensitivity in patients with cetuximab-resistant cancers (35). GRB2 amplification is observed in esophageal squamous cell carcinoma and is significantly involved in lymph node metastases (36). Phosphatidylinositol 3-kinase (PI3K), which may promote cancer cell survival, is an important therapeutic target in cancer (37). PIK3R1, the inhibitory subunit of PI3K, often undergoes mutations in endometrial cancer (38,39). The enrichment results also revealed the these 4 DEGs were enriched in ErbB signaling pathway, which is involved in regulating cell survival and adhesion (40). Silencing TAK1 in bladder cancer cells led to the upregulation of EGFR and downregulation of ERBB2, GRB2 and PIK3R1, suggesting an inhibitory effect of silencing TAK1 on the metastasis of bladder cancer cells via modulation of the ErbB signaling pathway. However, the explicit molecular mechanisms require further research.
In summary, the present study demonstrated that silencing TAK1 induced 43 downregulated and 21 upregulated genes. Silencing TAK1 may lead to decreased cellular response to chemical stimuli via downregulating CYP3A5, MAP3K7, GSTA1, ALDH1A1, ARRB1 and ITPR1, as well as increased cancer cell developmental processes via upregulating EGFR, EGR1 and PTHLH. In addition, silencing TAK1 may exert regulatory effects on bladder cancer metastasis and other various cancer types via regulating the expression of EGFR, ERBB2, GRB2 and PIK3R1.
Acknowledgements
The present study was supported by the Zhejiang Provincial Natural Science Foundation (grant no., LY13H050001).
References
McCarthy N: Bladder cancer: Seemingly similar. Nat Rev Cancer. 14:214–215. 2014. View Article : Google Scholar : PubMed/NCBI | |
Siegel R, Ma J, Zou Z and Jemal A: Cancer statistics, 2014. CA Cancer J Clin. 64:9–29. 2014. View Article : Google Scholar : PubMed/NCBI | |
Tomlinson D, Baldo O, Harnden P and Knowles M: FGFR3 protein expression and its relationship to mutation status and prognostic variables in bladder cancer. J Pathol. 213:91–98. 2007. View Article : Google Scholar : PubMed/NCBI | |
Tomlinson DC, Hurst CD and Knowles MA: Knockdown by shRNA identifies S249C mutant FGFR3 as a potential therapeutic target in bladder cancer. Oncogene. 26:5889–5899. 2007. View Article : Google Scholar : PubMed/NCBI | |
Karashima T, Sweeney P, Kamat A, Huang S, Kim SJ, Bar-Eli M, McConkey DJ and Dinney CP: Nuclear factor-kappaB mediates angiogenesis and metastasis of human bladder cancer through the regulation of interleukin-8. Clin Cancer Res. 9:2786–2797. 2003.PubMed/NCBI | |
Salazar L, Kashiwada T, Krejci P, Meyer AN, Casale M, Hallowell M, Wilcox WR, Donoghue DJ and Thompson LM: Fibroblast growth factor receptor 3 interacts with and activates TGFβ-activated kinase 1 tyrosine phosphorylation and NFκB signaling in multiple myeloma and bladder cancer. Plos One. 9:e864702014. View Article : Google Scholar : PubMed/NCBI | |
Reynolds A, Leake D, Boese Q, Scaringe S, Marshall WS and Khvorova A: Rational siRNA design for RNA interference. Nat Biotechnol. 22:326–330. 2004. View Article : Google Scholar : PubMed/NCBI | |
Folini M, Pennati M and Zaffaroni N: RNA interference-mediated validation of genes involved in telomere maintenance and evasion of apoptosis as cancer therapeutic targets. Method Mol Biol. 487:303–330. 2009. | |
Barrett T and Edgar R: Gene expression omnibus: Microarray data storage, submission, retrieval and analysis. Methods Enzymol. 411:352–369. 2006. View Article : Google Scholar : PubMed/NCBI | |
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U and Speed TP: Exploration, normalization and summaries of high density oligonucleotide array probe level data. Biostatistics. 4:249–264. 2003. View Article : Google Scholar : PubMed/NCBI | |
Gautier L, Cope L, Bolstad BM and Irizarry RA: Affy-analysis of affymetrix genechip data at the probe level. Bioinformatics. 20:307–315. 2004. View Article : Google Scholar : PubMed/NCBI | |
Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 3:Article32004.PubMed/NCBI | |
Cheng L, Lin H, Hu Y, Wang J and Yang Z: Gene function prediction based on the gene ontology hierarchical structure. Plos One. 9:e1071872014. View Article : Google Scholar : PubMed/NCBI | |
Altermann E and Klaenhammer TR: PathwayVoyager: Pathway mapping using the Kyoto encyclopedia of genes and genomes (KEGG) database. BMC Genomics. 6:602005. View Article : Google Scholar : PubMed/NCBI | |
Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC and Lempicki RA: DAVID: Database for annotation, visualization and integrated discovery. Genome Biol. 4:P32003. View Article : Google Scholar : PubMed/NCBI | |
Fogel GB, Weekes DG, Varga G, Dow ER, Craven AM, Harlow HB, Su EW, Onyia JE and Su C: A statistical analysis of the TRANSFAC database. Biosystems. 81:137–154. 2005. View Article : Google Scholar : PubMed/NCBI | |
Zhao M, Sun J and Zhao Z: TSGene: A web resource for tumor suppressor genes. Nucleic Acids Res. 41:D970–D976. 2013. View Article : Google Scholar : PubMed/NCBI | |
Chen JS, Hung WS, Chan HH, Tsai SJ and Sun HS: In silico identification of oncogenic potential of fyn-related kinase in hepatocellular carcinoma. Bioinformatics. 29:420–427. 2013. View Article : Google Scholar : PubMed/NCBI | |
Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, Lin J, Minguez P, Bork P, von Mering C and Jensen LJ: STRING v9.1: Protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 41:D808–D815. 2013. View Article : Google Scholar : PubMed/NCBI | |
Kohl M, Wiese S and Warscheid B: Cytoscape: Software for visualization and analysis of biological networks. Methods Mol Biol. 696:291–303. 2011. View Article : Google Scholar : PubMed/NCBI | |
Lah L, Podobnik B, Novak M, Korošec B, Berne S, Vogelsang M, Kraševec N, Zupanec N, Stojan J, Bohlmann J and Komel R: The versatility of the fungal cytochrome P450 monooxygenase system is instrumental in xenobiotic detoxification. Mol Microbiol. 81:1374–1389. 2011. View Article : Google Scholar : PubMed/NCBI | |
Karpusas M, Axarli I, Chiniadis L, Papakyriakou A, Bethanis K, Scopelitou K, Clonis YD and Labrou NE: The interaction of the chemotherapeutic drug chlorambucil with human glutathione transferase A1-1: Kinetic and structural analysis. Plos One. 8:e563372013. View Article : Google Scholar : PubMed/NCBI | |
Nebert DW and Dalton TP: The role of cytochrome P450 enzymes in endogenous signalling pathways and environmental carcinogenesis. Nat Rev Cancer. 6:947–960. 2006. View Article : Google Scholar : PubMed/NCBI | |
Fritz JH, Ferrero RL, Philpott DJ and Girardin SE: Nod-like proteins in immunity, inflammation and disease. Nat Immunol. 7:1250–1257. 2006. View Article : Google Scholar : PubMed/NCBI | |
Boorjian S, Tickoo SK, Mongan NP, Yu H, Bok D, Rando RR, Nanus DM, Scherr DS and Gudas LJ: Reduced lecithin retinol acyltransferase expression correlates with increased pathologic tumor stage in bladder cancer. Clin Cancer Res. 10:3429–3437. 2004. View Article : Google Scholar : PubMed/NCBI | |
Wincent J, Schoumans J and Anderlid BM: De novo deletion of chromosome 11q13.4-q14.3 in a boy with microcephaly, ptosis and developmental delay. Eur J Med Genet. 53:50–53. 2010. View Article : Google Scholar : PubMed/NCBI | |
Ferris CD, Huganir RL, Supattapone S and Snyder SH: Purified inositol 1,4,5-trisphosphate receptor mediates calcium flux in reconstituted lipid vesicles. Nature. 342:87–89. 1989. View Article : Google Scholar : PubMed/NCBI | |
Boehning D, Patterson RL, Sedaghat L, Glebova NO, Kurosaki T and Snyder SH: Cytochrome c binds to inositol (1,4,5) trisphosphate receptors, amplifying calcium-dependent apoptosis. Nat Cell Biol. 5:1051–1061. 2003. View Article : Google Scholar : PubMed/NCBI | |
Nicholson R, Gee J and Harper M: EGFR and cancer prognosis. Eur J Cancer. 37:(Supply 4). S9–S15. 2001. View Article : Google Scholar : PubMed/NCBI | |
Black PC, Brown GA, Dinney CP, Kassouf W, Inamoto T, Arora A, Gallagher D, Munsell MF, Bar-Eli M, McConkey DJ and Adam L: Receptor heterodimerization: A new mechanism for platelet-derived growth factor induced resistance to anti-epidermal growth factor receptor therapy for bladder cancer. J Urol. 185:693–700. 2011. View Article : Google Scholar : PubMed/NCBI | |
Chadalapaka G, Jutooru I, Burghardt R and Safe S: Drugs that target specificity proteins downregulate epidermal growth factor receptor in bladder cancer cells. Mol Cancer Res. 8:739–750. 2010. View Article : Google Scholar : PubMed/NCBI | |
Urosevic J, Garcia-Albéniz X, Planet E, Real S, Céspedes MV, Guiu M, Fernandez E, Bellmunt A, Gawrzak S, Pavlovic M, et al: Colon cancer cells colonize the lung from established liver metastases through p38 MAPK signalling and PTHLH. Nat Cell Biol. 16:685–694. 2014. View Article : Google Scholar : PubMed/NCBI | |
Zheng L, Pu J, Jiang G, Weng M, He J, Mei H, Hou X and Tong Q: Abnormal expression of early growth response 1 in gastric cancer: Association with tumor invasion, metastasis and heparanase transcription. Pathol Int. 60:268–277. 2010. View Article : Google Scholar : PubMed/NCBI | |
Yarden Y: The EGFR family and its ligands in human cancer. signalling mechanisms and therapeutic opportunities. Eur J Cancer. 37:(Supply 4). S3–S8. 2001. View Article : Google Scholar : PubMed/NCBI | |
Yonesaka K, Zejnullahu K, Okamoto I, Satoh T, Cappuzzo F, Souglakos J, Ercan D, Rogers A, Roncalli M, Takeda M, et al: Activation of ERBB2 signaling causes resistance to the EGFR-directed therapeutic antibody cetuximab. Sci Transl Med. 3:99ra862011. View Article : Google Scholar : PubMed/NCBI | |
Li LY, Li EM, Wu ZY, Cao HH, Shen JH, Xu XE, Chen B, Wu JY and Xu LY: Overexpression of GRB2 is correlated with lymph node metastasis and poor prognosis in esophageal squamous cell carcinoma. Int J Clin Exp Pathol. 7:3132–3140. 2014.PubMed/NCBI | |
Engelman JA: Targeting PI3K signalling in cancer: Opportunities, challenges and limitations. Nat Rev Cancer. 9:550–562. 2009. View Article : Google Scholar : PubMed/NCBI | |
Urick ME, Rudd ML, Godwin AK, Sgroi D, Merino M and Bell DW: PIK3R1 (p85α) is somatically mutated at high frequency in primary endometrial cancer. Cancer Res. 71:4061–4067. 2011. View Article : Google Scholar : PubMed/NCBI | |
Cheung LW, Hennessy BT, Li J, Yu S, Myers AP, Djordjevic B, Lu Y, Stemke-Hale K, Dyer MD, Zhang F, et al: High frequency of PIK3R1 and PIK3R2 mutations in endometrial cancer elucidates a novel mechanism for regulation of PTEN protein stability. Cancer Discov. 1:170–185. 2011. View Article : Google Scholar : PubMed/NCBI | |
Yarden Y and Sliwkowski MX: Untangling the ErbB signalling network. Nat Rev Mol Cell Biol. 2:127–137. 2001. View Article : Google Scholar : PubMed/NCBI |