Bioinformatics analysis of the gene expression profile in Bladder carcinoma

Bioinformatics analysis of the gene expression profile in Bladder carcinoma

Jing Xiao; Qiu Yiqing

Department of Urology, The Second Affiliated Hospital of Zhejiang University, Hangzhou, Zhejiang Province, China

^{Send correspondence to}

ABSTRACT

Bladder carcinoma, which has the ninth highest incidence among malignant tumors in the world, is a complex, multifactorial disease. The malignant transformation of bladder cells results from DNA mutations and alterations in gene expression levels. In this work, we used a bioinformatics approach to investigate the molecular mechanisms of bladder carcinoma. Biochips downloaded from the Gene Expression Omnibus (GEO) were used to analyze the gene expression profile in urinary bladder cells from individuals with carcinoma. The gene expression profile of normal genomes was used as a control. The analysis of gene expression revealed important alterations in genes involved in biological processes and metabolic pathways. We also identified some small molecules capable of reversing the altered gene expression in bladder carcinoma; these molecules could provide a basis for future therapies for the treatment of this disease.

Keywords: Bladder carcinoma, differential expression, expression profile, small molecular mimic, susceptibility forecast.

Introduction

Bladder carcinoma is one of the most common malignant tumors and has a high death rate in China (Dai et al., 2011). World-wide, bladder carcinoma has the ninth highest incidence among malignant tumors (Parkin, 2008). Bladder carcinoma results from cumulative, long-term interactions between genetic and environmental elements and is a complex, multifactorial process involving gene mutations and progressive cellular damage. Bladder carcinoma is much more common in men than in women and predisposing factors include cigarette smoking, long-term contact with chemical products and race (Parkin, 2008).

The malignant transformation of bladder cells originates in DNA mutations and alterations in gene expression. Most studies of the etiology of bladder carcinoma have focused on genetic mutations, particularly in oncogenes such as H-Ras, C-Myc, HER-2, BcL-2, erbB-2 and fibroblast growth factor receptor 3 (FGFR3) (Sanchez-Carbayo and Cordon-Cardo, 2007; Cordon-Cardo, 2008), as well as the inactive of anti-oncogenes and the deletion or loss of heterozygosity (LOH) in associated chromosomal regions for genes such as p53, pRb, p21 and p16 (Shariat et al., 2004). The invasiveness and metastasis of bladder carcinoma involve the abnormal expression of growth factors and related receptors such as vascular endothelial growth factor (VEGF) (Wang et al., 2000) and fibroblast growth factor receptor (FGFR) (Knowles, 2008).

A detailed understanding of multi-gene hereditary diseases, such as cancer, is no longer feasible with traditional single gene studies but requires more powerful tools such as large-scale gene expression analysis provided by recent technological advances. The application of such technologies to the study of bladder carcinoma should be useful in identifying new tumor marker genes that will allow the early diagnosis of bladder cancer.

In this study, we used the gene expression profile of GSE27448 to investigate differential gene expression in bladder carcinoma and to search for small molecules that may reverse the altered gene expression. The availability of high-throughput gene expression data and their analysis by computational bioinformatics may be useful in identifying molecules that could provide a basis for developing new therapies for treating bladder carcinoma.

Materials and Methods

Gene expression in bladder carcinoma cells and normal cells

To investigate the differential gene expression in bladder carcinoma compared to normal cells and the possible mechanisms involved in the initiation of bladder cancer, we analyzed the gene expression profile of bladder carcinoma and normal cells available from the Gene Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo/); accession number: GSE27448). This dataset was deposited by Zaravinos et al. (2011) and obtained by cell culture, RNA extraction and individual on-chip analysis. The dataset contains five normal cell chips and nine bladder carcinoma cell chips. The platform used was GPL2895 (GE Healthcare/Amersham Biosciences CodeLink Human Whole Genome Bio array). The original files and the platform probe annotation files were downloaded.

Identification of differentially expressed genes (DEGs)

The original data were classified as bladder carcinoma and control (i.e., normal) groups and were analyzed using R software (v.2.13.0) (R Development Core Team, 2008). Initially, the RMA (robust multichip averaging) method (Irizarry et al., 2003) was applied to normalize the data on different chips, and then Limma (Smyth, 2004), a linear regression model software, was used to compare the differential expression on different classes of chips. A value of p < 0.05 was used as the cut-off criterion to select genes that were differentially expressed in bladder carcinoma relative to the control group.

Gene ontology (GO) analysis of DEGs

To investigate the DEGs at a functional level, DAVID (Database for Annotation, Visualization and Integrated Discovery) was used to cluster the genes according to the Gene Ontology (GO) (Ashburner et al., 2000) categories of cellular component, biological process and molecular function (Huang da et al., 2009a,b).

Pathway enrichment analysis

To gain more insights of the changes at a functional level, we investigated the dysregulated biological pathways in bladder carcinoma. All of the metabolic and non-metabolic pathways available from the canonical KEGG PATHWAY DATABASE were used as DAVID inputs for KEGG PATHWAY cluster analysis (Huang da et al., 2009a,b). A value of p < 0.1 and at least two DEGs contained in a pathway were chosen as the cut-off criteria.

Identification of small molecules

The connectivity map (CMap) database, which contains whole genomic expression profiles for small active molecular inferences, consists 6,100 classes of small molecular interference experiments and 7,056 expression profiles (Lamb et al., 2006). We analyzed genes that were differentially expressed between normal and bladder carcinoma cells and contrasted them with genes involved in small molecular interference in the CMap database in an attempt to identify small molecules associated with these DEGs. The DEGs were classified into up-regulated and down-regulated groups and the top 500 significant probes in each group were chosen for GSEA (Gene Set Enrichment Analysis)analysis; these were then compared to the

Xiao and Yiqing

DEGs after small molecule treatment. Finally, the enrichment values were calculated. These values varied between -1 and 1 and reflected the similarity among the genes: values closer to -1 indicated greater similarity between the genes, i.e., the small molecule could imitate the status of normal cells, whereas values closer to 1 indicated that the small molecule could imitate the status of bladder carcinoma cells.

Results

Identification of DEGs

The classic t-test in the Limma package was used to analyze the gene expression profiles of bladder carcinoma and normal cells and identify the DEGs in bladder carcinoma cells. A value of p < 0.05 was used as the significant threshold for DEGs. Based on these criteria, 6562 gene probes (corresponding to 2779 genes) were found to have an altered expression.

GO clustering of DEGs

Functional classification of these 2779 DEGs with the online biological classification tool DAVID and a statistical cut-off criterion of p < 0.001 indicated significant enrichment of these genes in various GO categories. Table 1 shows that in bladder carcinoma DEGs were enriched in cellular components related to the cytoskeleton and locomotion (Z discs, I band, myofibrils, contractile fibers, sarcolemma and sarcomere) and signal transduction. There were also changes in cellular components related to protein expression, e.g., the ribonucleoprotein complex, ribosomes, nucleoli, heterochromatin, nuclear lumen and cytosolic large ribosomal subunit.

Table 2 shows the clusters obtained when these DEGs were classified according to biological process (p < 0.01). This analysis revealed inter-and intracellular changes in components related to transportation, such as the regulation of intracellular transport, protein localization, the regulation of nucleocytoplasmic transport, the regulation of intracellular protein transport, negative regulation of intracellular transport and negative regulation of nucleocytoplasmic transport. These changes meant that signal molecules could not be transported to their target site, thus altering intercellular signal transduction. Similar changes were observed in other biological processes, such as intracellular signaling cascades, small GTPase-mediated signal transduction and the protein kinase cascade. Other changes affected protein translation and post-translational modifications, such as amino acid auto-phosphorylation, translation, post-translational elongation, protein complex assembly and protein complex biogenesis. Biological processes related to cell development, e.g., peripheral nervous system development, auditory receptor cell differentiation, mechanoreceptor differentiation and auditory receptor cell development, were also affected.

Table 3 shows the DEGs clustered according to molecular function (p < 0.01). The changes in gene expression affected mainly gene transcription and translation, e.g., nucleoside-triphosphatase regulator activity, purine nucleotide binding and double-stranded RNA binding, and signal transduction, e.g., insulin receptor binding, GTPase regulator activity and GTPase activator activity.

Altered biological pathways in bladder carcinoma cells

The gene expression profile changed significantly during the progression of bladder carcinoma. Some genes changed significantly under pathological conditions. These DEGs were selected for KEGG pathway enrichment analysis and subsequent identification of the altered pathways in bladder carcinoma cells. Based on a threshold value of p < 0.1 and a gene count>2we identified eight dysregulated pathways (Table 4). The most significant enrichment involved the complement and coagulation cascades (p = 0.002326). Some of the significant pathways were related to signaling transduction, such as PPAR signaling pathway (p = 0.052499) and SNARE interactions in vesicular transport (p = 0.075573).

Identification of related active small molecules

The DEGs were classified into up-regulated and down-regulated gene groups followed by GSEA analysis and matching after small molecule treatment in the CMap database. This analysis identified some small molecules that could reverse the status of bladder carcinoma cells. Table 5 shows the 20 most significant small molecules identified by this procedure.

Thioguanosine (enrichment = 0.865) and tyloxapol (enrichment = 0.857) can partially imitate the carcinoma status of bladder cells, i.e., these small molecules may be strong inducers of bladder cancer. In contrast, adiphenine (enrichment = -0.962) and viomycin (enrichment = -0.947) can imitate the status of normal cells, i.e., they can reverse the abnormal status and are therefore potentially lead molecules for developing new therapeutic drugs for treating bladder carcinoma. Further investigation of these small molecules may clarify their role in the pathogenesis of bladder cancer.

Discussion

Bladder cancer involves the proliferation of malignant intrabladder cells and is the fourth most common solid tumor in men and the seventh most common in women worldwide, with 350,000 new cases each year. The American Cancer Society (AGS) statistics indicate that in 2006 there were 61,420 new cases of bladder cancer and 13,060 deaths. In China, bladder cancer is also commonly associated with urinary tract malignant tumors. The morbidity in 2005 was 4/100 in males and 1.5/100 in females. The steady increase in the morbidity of bladder cancer in some Chinese cities in recent years indicates a need for additional research on this disease.

The main biological processes that showed DEGs based on GO cluster analysis were signal transduction, physical immunity and disease initiation. Nearly 2800 genes were differentially expressed in bladder carcinoma compared with normal cells. These DEGs may be important for investigating the mechanism of disease development and may be useful targets for treating bladder carcinoma, hence the potential interest in studying these genes. The results of GO cluster analysis and pathway enrichment analysis of DEGs suggested that the development of bladder carcinoma may involve important changes in cellular signal transduction.

As a signal molecular receptor, SNARE plays a key role in signal transduction (Skalski et al., 2011). Changes in SNARE interactions in vesicular transport directly affect cellular transduction (Fasshauer, 2003; Han et al., 2009). Peroxisome proliferator-activated receptors (PPAR), which are nuclear hormone receptors (Umemoto and Fujiki, 2012), regulate gene expression and changes in the PPAR signaling pathway result in abnormal gene expression. These signaling pathways regulate cell growth (Sertznig et al., 2008), differentiation and proliferation (Peters et al., 2012) by modulating the expression of downstream genes. The changes mediated by these signaling pathways ultimately affect protein synthesis, such as the ribosomal pathway and development (progesteronemediated oocyte maturation) (Fauconnet et al., 2002). Such changes are closely related to the onset of bladder carcinoma (Yoshimura et al., 2003).

Based on the DEGs and small molecule interference data we identified a set of small molecules that can imitate the normal cell status, i.e., they can reverse the abnormal gene expression of bladder carcinoma. These small molecules may provide lead compounds for developing new therapies for treating bladder carcinoma. For example, adiphenine inhibits nicotinic receptors, thereby decreasing the frequency of acetylcholine-induced single-channel currents and cluster duration (the latter by 36-fold at 100 µM) while increasing the decay rate (IC₅₀ = 15 µM) of macroscopic currents and accelerating desensitization from the open state; however, prior application of the drug to resting receptors is required for these effects to be seen.(Spitzmaul et al., 2009). The nicotinic receptor is reported to have an important role in various types of cancer (Russo et al., 2012). Hence, adiphenine is a potential therapeutic antagonist for bladder carcinoma. The small molecules that imitate the normal cell status identified in this study may therefore offer new fields of research in bladder cancer therapy.

Received: July 26, 2012; Accepted: December 23, 2012.

Associate Editor: Jeremy A. Squire

License information: This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

• Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al. (2000) Gene ontology: Tool for the unification of biology. The Gene Ontology Consortium. Nat Genet 25:25-29.
• Cordon-Cardo C (2008) Molecular alterations associated with bladder cancer initiation and progression. Scand J Urol Nephrol 218:154-165.
• Dai QS, He HC, Cai C, Chen JH, Han ZD, Qin GQ, Liang YX and Zhong WD (2011) Multicenter case-control study of the relationship between smoking and bladder cancer in China. Zhonghua Yi Xue Za Zhi 91:2407-2410 [In Chinese]
• Fasshauer D (2003) Structural insights into the SNARE mechanism. Biochim Biophys Acta 1641:87-97.
• Fauconnet S, Lascombe I, Chabannes E, Adessi GL, Desvergne B, Wahli W and Bittard H (2002) Differential regulation of vascular endothelial growth factor expression by peroxisome proliferator-activated receptors in bladder cancer cells. J Biol Chem 277:23534-23543.
• Han C, Chen T, Yang M, Li N, Liu H and Cao X (2009) Human SCAMP5, a novel secretory carrier membrane protein, facilitates calcium-triggered cytokine secretion by interaction with SNARE machinery. J Immunol 182:2986-2996.
• Huang da W, Sherman BT and Lempicki RA (2009a) Bioinformatics enrichment tools: Paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res 37:1-13.
• Huang da W, Sherman BT and Lempicki RA (2009b) Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc 4:44-57.
• Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U and Speed TP (2003) Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics 4:249-264.
• Knowles MA (2008) Novel therapeutic targets in bladder cancer: Mutation and expression of FGF receptors. Future Oncol 4:71-83.
• Lamb J, Crawford ED, Peck D, Modell JW, Blat IC, Wrobel MJ, Lerner J, Brunet JP, Subramanian A, Ross KN, et al. (2006) The Connectivity Map: Using gene-expression signatures to connect small molecules, genes, and disease. Science 313:1929-1935.
• Parkin DM (2008) The global burden of urinary bladder cancer. Scand J Urol Nephrol 218:12-20.
• Peters JM, Shah YM and Gonzalez FJ (2012) The role of peroxisome proliferator-activated receptors in carcinogenesis and chemoprevention. Nat Rev Cancer 12:181-195.
• R Development Core Team (2008) R: A language and environment for statistical computing. Royal Foundation for Statistical Computing, Austria, 1697 pp.
• Russo P, Cardinale A, Margaritora S and Cesario A (2012) Nicotinic receptor and tobacco-related cancer. Life Sci 91:1087-1092.
• Sanchez-Carbayo M and Cordon-Cardo C (2007) Molecular alterations associated with bladder cancer progression. Semin Oncol 34:75-84.
• Sertznig P, Seifert M, Tilgen W and Reichrath J (2008) Peroxisome proliferator-activated receptors (PPARs) and the human skin: Importance of PPARs in skin physiology and dermatologic diseases. Am J Clin Dermatol 9:15-31.
• Shariat SF, Tokunaga H, Zhou J, Kim J, Ayala GE, Benedict WF and Lerner SP (2004) p53, p21, pRB, and p16 expression predict clinical outcome in cystectomy with bladder cancer. J Clin Oncol 22:1014-1024.
• Skalski M, Sharma N, Williams K, Kruspe A and Coppolino MG (2011) SNARE-mediated membrane traffic is required for focal adhesion kinase signaling and Src-regulated focal adhesion turnover. Biochim Biophys Acta 1813:148-158.
• Smyth GK (2004) Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol 3:e3.
• Spitzmaul G, Gumilar F, Dilger JP and Bouzat C (2009) The local anaesthetics proadifen and adiphenine inhibit nicotinic receptors by different molecular mechanisms. Br J Pharmacol 157:804-817.
• Umemoto T and Fujiki Y (2012) Ligand-dependent nucleo-cytoplasmic shuttling of peroxisome proliferator-activated receptors, PPARa and PPARy. Genes Cells 17:576-596.
• Wang S, Xia T, Zhang Z, Kong X, Zeng L, Mi P and Xue Z (2000) Expression of VEGF and tumor angiogenesis in bladder cancer. Zhonghua Wai Ke Za Zhi 38:34-36 [In Chinese]
• Yoshimura R, Matsuyama M, Segawa Y, Hase T, Mitsuhashi M, Tsuchida K, Wada S, Kawahito Y, Sano H and Nakatani T (2003) Expression of peroxisome proliferator-activated receptors (PPARs) in human urinary bladder carcinoma and growth inhibition by its agonists. Int J Cancer 104:597-602.
• Zaravinos A, Lambrou GI, Boulalas I, Delakas D and Spandidos DA (2011) Identification of common DEGs in urinary bladder cancer. PLoS One 6:e18135.

Publication Dates

• Publication in this collection
  12 Apr 2013
• Date of issue
  2013

History

• Received
  26 July 2012
• Accepted
  23 Dec 2012