Introduction

Molecular Medicine Reports

1791-2997 1791-3004

D.A. Spandidos

10.3892/mmr.2017.7726

mmr-16-06-8709

Articles

Identification and bioinformatics analysis of miRNAs associated with human muscle invasive bladder cancer

Chen

Liang

1 * Yuan

Lushun

1 * Wang

Gang

1 Cao

Rui

1 Peng

Jianping

1 Shu

1 2 Qian

Guofeng

3 Wang

Xinghuan

1 Xiao

1 4

1Department of Urology, Zhongnan Hospital of Wuhan University, Wuhan, Hubei 430071, P.R. China 2Department of Urology, The Central Hospital of Wuhan, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei 430071, P.R. China 3Department of Endocrinology, The First Affiliated Hospital of Zhejiang University, Hangzhou, Zhejiang 310000, P.R. China 4Laboratory of Precision Medicine, Zhongnan Hospital of Wuhan University, Wuhan, Hubei 430071, P.R. China

Correspondence to: Dr Xinghuan Wang or Dr Yu Xiao, Department of Urology, Zhongnan Hospital of Wuhan University, 169 Donghu Road, Wuhan, Hubei 430071, P.R. China, E-mail: wangxinghuan@whu.edu.cn, E-mail: yu.xiao@whu.edu.cn *

Contributed equally

122017

04102017

16 6 8709 8720 14022017 11082017

2017

This is an open access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations are made.

Accumulated evidence has indicated that micro (mi)RNAs play vital roles in the occurrence and development of human muscle invasive bladder cancer (MIBC), however, little is known about the miRNAs' regulatory networks. In the present study, the authors aimed to use bioinformatics analysis to identify the key miRNAs and potential target genes, as well as studying the underlying mechanisms for MIBC. They collected several human MIBC tissues to generate a miRNA expression analysis by microarray analysis comparing with normal bladder tissues, identifying 104 differentially expressed miRNAs (102 were downregulated and 2 were upregulated) and predicted 11,884 putative target genes of the dysregulated miRNAs. To understand the function of dysregulated miRNAs in the development of MIBC, networks among miRNAs and genes, gene ontologies and pathways were built. The subsequent bioinformatics analysis indicated that the mitogen-associated protein kinase (MAPK) signaling pathway, apoptosis and pathways in cancer and the cell cycle, were significantly enriched Overall, these results provided comprehensive information on the biological function of dysregulated miRNAs in the development of MIBC. The identification of miRNAs and their putative targets may offer new diagnostic and therapeutic strategies for human muscle invasive bladder cancer.

bioinformatics analysis muscle invasive bladder cancer microRNA microarray miRNAs-gene-networks miRNAs-GO-networks

Introduction

Human bladder cancer (BCa), one of the most common genitourinary malignancies arising from mucous membrane, accounts for >3% of all malignant tumors around the world and more and more people are diagnosed each year (1). Early stage of bladder cancers do not usually cause symptoms and, with regard to those of late stage, it still frequently recurs and gradually progresses into muscle invasive BCa (2). The gold standard of BCa diagnosis is based on cystoscopy, which is invasive and relatively expensive. Currently, non-invasive and specific markers are used for urinary cytology, however, this method is not sensitive to detection of low-grade BCa. Thus, new, highly sensitive and specific urine-based diagnostic tools are particularly attractive because urine is a promising and readily available source for molecular markers, including RNA.

As small nucleotides of RNA, by binding to complementary sequences in the 3′-untranslated regions of specific mRNAs, micro (mi)RNAs inhibit the translation of specific target genes (3). In previous years, many cardinal cellular and physiological processes have been demonstrated to be associated with altered change in miRNAs expression levels (4–6). miRNAs are involved in the regulation of various cellular processes including cellular differentiation, cell cycle progression and apoptosis. Moreover, miRNAs, as important factors in tumorigenesis and metastasis and their expression signatures are associated with the prognosis and progression in a variety of cancers (7–9). Therefore, the authors could suppose that those miRNAs of which the expressions are changed significantly in tumors relative to normal tissues may have influence on tumor progression (10). Moreover, identifying target genes united with differentially expressed miRNAs might illustrate the roles of miRNAs in cancer biology (11). Previous studies have reported that we could distinguish malignant or normal tissue, as well as various tumor entities by miRNA expression profiles (12–16). As highly stable molecules, the authors could quantify miRNAs in tissues and body fluids, which makes them considered as promising cancer biomarkers. It has been reported that diagnosis, stage and sub-classification of cancer can be explored by differentially expressed miRNAs, which also can predict treatment efficacy and prognosis (4). However, the molecular mechanisms are not yet clearly elucidated, which makes it necessary to identify novel miRNAs.

In the present study, the group performed an MIBC-related miRNA expression profile, so as to address the functional roles of miRNAs in MIBC. Several human MIBC tissues and normal bladder tissues were conducted to conduct a microarray analysis (GEO accession No. GSE76211), revealing several differentially expressed miRNAs. The putative target genes of dysregulated miRNAs were predicted using miRNA databases. Then the networks among miRNAs and genes, gene ontologies (GOs) and pathways were built. The purpose of the current study was to identify candidate predictive tumor-associated miRNAs in MIBC patients.

Materials and methods <sec> <title>Preparation for human bladder samples

Three MIBC tissues samples were collected from patients after surgery at Zhongnan Hospital of Wuhan University (Wuhan, China). Three normal bladder tissue samples were collected from donors by accidental death. Information of the MIBC patients and donors was listed in Table I. Those samples used in the study have been described in previous publication (17–19). Briefly, two pathologists confirmed the histology diagnosis independently and all the tissues were snap-frozen for total RNA isolation at liquid nitrogen after excision from operation room. Informed consent was obtained for surgery patients and normal donors from the patients and their relatives, respectively. The Ethics Committee at Zhongnan Hospital of Wuhan University (Wuhan, China) approved the experiments using human bladder tissue samples for RNA isolation analysis (approval number: 2015029). All methods used for human bladder tissue samples were performed in accordance with the approved guidelines and regulations.

RNA extraction

Based on the manufacturer's protocol, total RNA was extracted from the frozen tissue block using RNeasy Mini kit (cat. no. 74101, Qiagen GmbH, Hilden, Germany), combined with QIAshredder (cat. no. 79654, Qiagen GmbH) using a centrifuge (cat. no. 5424, Eppendorf, Hamburg, Germany). In order to remove genomic DNA, DNase I digestion (cat. no. 79253, Qiagen GmbH), DNase I digestion (cat. no. 79254, Qiagen GmbH) was used in each RNA preparation.

miRNA microarray

After assessing RNA quality and quantity, the miRNAs microarray analysis (Affymetrix microRNA 4.0 Array, Affymetrix, Inc., Santa Clara, CA, USA) was performed according to the manufacturer's instructions. Briefly, 1 µg of total RNA was labeled with Biotin using the FlashTag Biotin HSR RNA Labeling kit (Genisphere LLC, Hatfield, PA, USA) and then hybridized overnight with the array, which was washed, stained, and read by an GeneChip Scanner 3000 7G (Affymetrix, Inc.). MiRNA microarray data GSE40355 used for validation were obtained from Gene Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo/). This dataset included eight normal bladder tissues samples, eight low grade BCa tissues samples and eight high grade BCa tissue samples. Among them, significant expressed miRNAs were screened out from high grade BCa tissues compared with normal bladder tissues.

Data analysis of miRNAs microarray

CEL-files of the raw data were first exported by Affymetrix GeneChip Command Console Software Version 4.0 (Affymetrix, Inc.) and then uploaded to the website of Gminix-Cloud Biotechnology Information (GCBI) by Genminix Informatics Co., Ltd. (Shanghai, China; http://www.gcbi.com.cn/gclib/html/index) for further analysis, including difference analysis of miRNAs profiles, prediction of miRNAs target genes, GO/pathway enrichment analysis, miRNAs-gene-network and miRNAs-GO-network analysis. The miRNAs array data used in the present paper has been uploaded to the NCBI Gene Expression Omnibus and the GEO accession number is GSE76211.

According to the GCBI online method description for difference analysis, the procedure for candidate miRNAs selection is as follows: When the number of samples in each group is no less than 3, SAM method is used for difference analysis. The authors implemented a series of steps to obtain the estimation of significance of difference and false discovery rate for every filtered gene:

i) Calculate the exchange factor s₀: Firstly, calculate the standard deviation for all genes, denote s^α as the α percentile for s_i. For the percentile value q₁ < q₂ … q₁₀₀ of the s_i, calculate the statistic:

vj=mad{diα=ri/(si+sα)|si=[qj,qj+1)}

Where mad denotes the mean absolute deviation. At last, α (denote as ᾶ) was selected to make the CV (coefficient of variation) of the v_j achieve minimum. Then, the exchange factor s₀ is s^ᾶ used.

ii) Calculate the statistic value (d Score) for every gene:

di=ri/(si+s0)

Where r_i reflects the difference in average level among different groups, s_i reflects the variation of sample population. See details in references (20,21)

iii) Calculate the order statistic:

d(1)≤d(2)≤···≤d(i)≤···≤d(p)

iv) In order to get the above statistic's estimate, the authors made a permutation method (a loop strategy through every sample, the total number no less than 1,000; the detail are omitted) the expected distribution of d score. The estimated statistic values are denoted as follows:

d(1)*≤d(2)*≤···≤d(i)*≤···≤d(p)*

v) The authors obtain the order statistic value under the permutation:

d¯(i)=∑i=11000d(i)*1000

vi) By calculating the maximum distance between the order statistic d_(i) and the expected order statistic d¯(i), the authors constructed a series of rejection regions for q-value. In fact, a grid of delta values was obtained by dividing 50 equivalent delta value for the above distance.

vii) For a fixed delta value, by computing the difference ∆_(i) ₌ d_(i) -d¯(i) the authors identified the nearest ∆_(i) for gene i. The cut-up is marked as min {∆_(i) _{≥ delta}} for positive gene and the cut-down max {∆_(i) _{≤ delta}} as for negative gene. The genes with differences above the cut-up value (we denote the number of these genes as R_(p)) were considered as significantly positive genes. While the genes with differences lower than the cut-down value were considered as significantly negative genes.

viii) Under the above cut-up and cut-down thresholds, the simulation of step VII was performed respectively on the statistics obtained from step V, such that the number of positive genes could be obtained under random state (≥1,000 permutations). The median of the 1,000 positive genes was estimated as the number of false positive genes, to allow the false discovery rate (FDR) to be estimated (the positive FDR=V(p)R(p) and the negative FDR was similar), and thus the proportion of false positive genes in the full set of positive genes.

ix) Finally, according to the definition of q-value (22), the authors obtained the q-value for the gene, i, by selecting the minimum of the FDR for the 50 delta values determined in step VII (every delta as a rejection region).

Results <sec> <title>Identification of differentially expressed miRNAs (DE-miRNAs) in BCa tissues

The obtained miRNA expression profiles (GSE76211) of BCa and normal bladder tissues were analyzed by the Affymetrix microRNA 4.0 Array, which contains 2,578 probes and can interrogate all mature miRNAs sequences in miRBase Release 20. The results revealed that 104 miRNAs were dysregulated in BCa group under the condition of ‘P<0.05 and fold change >1.5’, compared with normal bladder group. Among them, 102 miRNAs were downregulated and 2 were upregulated (Fig. 1). All of the dysregulated miRNA were listed in Table II.

Identification of putative target genes

The current study has identified 104 miRNAs that were significantly dysregulated in BCa tissues compared with normal bladder tissues. As miRNAs play their functional roles by regulating target genes expression at the posttranscriptional level, the authors predicted the target genes of dysregulated miRNAs using GCBI online tools, which were mainly based on the algorithms of miRanda and TargetScan. A total of 11,884 genes were predicted as putative target genes of dysregulated miRNAs.

GO/pathway enrichment analysis of putative target genes of dysregulated miRNAs

To understand the role of miRNAs in cancer development, GO and pathway enrichment analysis were performed. The data in Fig. 2A indicated that the top 10 dysregulated GOs were ‘transcription, DNA-dependent’, ‘signal transduction, positive regulation of transcription from RNA polymerase II promoter’, ‘regulation of transcription, DNA-dependent’, ‘axon guidance’, ‘positive regulation of transcription, DNA-dependent’, ‘apoptotic process’, ‘negative regulation of transcription from RNA polymerase II promoter’, ‘synaptic transmission’ and ‘negative regulation of transcription, DNA-dependent’. GO analysis obviously suggests that many dysregulated miRNAs may contribute to tumorigenesis of bladder through many important functions such as transcription regulation, signal transduction as well as apoptotic process. Combined with the KEGG database, the authors analyzed the pathways involving the putative target genes. As illustrated in Fig. 2B, the top ten dysregulated pathways were the mitogen-associated protein kinase (MAPK) signaling pathway, apoptosis, pathways in cancer, cell cycle, p53 signaling pathway, calcium signaling pathway, Wnt signaling pathway, adherens junction, focal adhesion and ErbB signaling pathway.

Pathway network analysis

Then, pathway network (Path-net) analysis was performed to draw an interaction network covering 36 significantly changed pathways (Fig. 3). Among them, the MAPK signaling pathway (degree=44), apoptosis (degree=30), pathways in cancer (degree=29) and cell cycle (degree=24) showed the highest degree, suggesting that these four pathways might play a core role in regulation of bladder cancer development.

miRNAs-gene-networks and miRNAs-GO-networks

Based on the significantly regulated GOs and pathways, the authors selected intersected genes and further constructed miRNAs-gene-networks and miRNAs-GO-networks to screen the key regulatory functions of the identified miRNAs and their target genes, respectively. As shown in Figs. 4 and 5, and Table III, the top rated nine miRNAs were hsa-miR-497-5p, hsa-miR-29a-3p, hsa-miR-124-3p, hsa-miR-4269, hsa-miR-145-5p, hsa-miR-204-5p, hsa-miR-4510, hsa-miR-6840-3p and hsa-miR-6722-3p. All of these miRNAs were downregulated in BCa tissues compared with normal bladder tissues. The microarray analysis demonstrated that the dysregulated miRNAs primarily play vital roles in various biological processes, including transcription regulation, apoptotic process, gene expression and signal transduction. Taken together, deregulation of certain miRNAs and several important pathways may be closely associated with human bladder cancer development.

Validation of candidate miRNAs

To confirm that the top rated nine miRNAs identified were indeed dysregulated in BCa tissues, the authors used another miRNA expression profiling GSE40355 including 8 normal bladder tissues samples, 8 low grade BCa tissues samples and 8 high grade BCa tissue samples for validation. In the dataset of GSE40355, the expression values of candidate miRNAs between normal bladder tissue samples with high grade BCa tissue samples were extracted for t test, and P<0.05 were considered statistically significant. In addition, 157 dysregulated miRNAs including 69 upregulated and 88 downregulated miRNAs were obtained by analysis of GSE40355. As presented in Fig. 6 and Table IV, among the top rated nine miRNAs screened from our miRNA microarray, four miRNAs involving hsa-miR-497-5p, hsa-miR-29a-3p, hsa-miR-145-5p and hsa-miR-204-5p were also significantly altered in GSE40355.

Discussion

In the current study, by comparing array-based miRNA expression profiling performed on three MIBC tissues and three normal bladder tissues, the integrated bioinformatic analysis from miRNA expression profiling identified that 104 miRNAs were differentially expressed (P<0.05, fold change >1.5), of which, 102 were downregulated and 2 were upregulated (Fig. 1 and Table II). Among the top 10 fold change miRNAs, many of them were closely linked to occurrence and development of bladder cancer and might play vital role as oncogenes and tumor suppressor genes, which has been published in several studies (23–27). Several studies have showed that miR-490-5p was significantly downregulated in BCa tissue samples compared with adjacent normal tissues (28). Low expression of miR-139-3p, miR-133, miR-145 and miR-1 in MIBC tissue samples compared with normal tissues was also reported by studies and played a functional role in bladder cancer cell lines (29–32).

Then, 11,884 genes were predicted as putative target genes of dysregulated miRNAs using target prediction method in GCBI online tools. According to the GO analysis, the predicted target genes were mainly involved in ‘transcription, DNA-dependent’, signal transduction, positive regulation of transcription from RNA polymerase II promoter, ‘regulation of transcription, DNA-dependent’ and axon guidance (Fig. 2A). Interestingly, the authors noticed two pairs of opposite GOs (Negative regulation of transcription, DNA-dependent vs. Positive regulation of transcription, DNA-dependent and ‘Negative regulation of transcription from RNA polymerase II promoter vs. Positive regulation of transcription from RNA polymerase II promoter’). miRNAs were involved in post-transcriptional regulation of gene expression (33), which played a key role in various cellular processes including cellular differentiation, cell cycle progression and apoptosis. Among the predictive target genes of miRNAs, there may be lots of oncogenes and tumor suppressor genes. Oncogenes could promote BCa cell proliferation by positive regulation of transcription, in contrast, tumor suppressor genes could suppress BCa cell viability by negatively regulation of transcription. Therefore, it can be inferred that miRNA regulation disorders might account for these biological behaviors. As for biological pathways, the MAPK signaling pathway, apoptosis, pathways in cancer, cell cycle, p53 signaling pathway, calcium signaling pathway, Wnt signaling pathway, adherens junction, focal adhesion and Erbb signaling pathway were the top 10 enriched pathways of the predicted target genes (Fig. 2B). Many studies including our group (17–19) have ever reported that these pathways such as the MAPK signaling pathway, Wnt signaling pathway, p53 signaling pathway and Erbb signaling pathway play a functional role in human bladder cancer cells (34,35). Besides, the Path-net analysis covering 36 significant pathways also showed that the MAPK signaling pathway, cell cycle, p53 signaling pathway, Wnt signaling pathway and calcium signaling pathway have a close correlation with bladder cancer, indicating that these pathways might play a key role in the development of human bladder cancer. In order to find out the key miRNAs, the authors conducted regulatory network analysis by overlapping significant miRNAs, pathways and GO analysis, revealing that the top nine miRNAs were hsa-miR-497-5p, hsa-miR-29a-3p, hsa-miR-124-3p, hsa-miR-4269, hsa-miR-145-5p, hsa-miR-204-5p, hsa-miR-4510, hsa-miR-6840-3p and hsa-miR-6722-3p. In order to validate that the nine candidate miRNAs identified were indeed dysregulated in BCa tissues, another miRNA microarray was performed by overlapping different expressed miRNAs and the top nine candidate miRNAs suggesting that four miRNAs involving hsa-miR-497-5p, hsa-miR-29a-3p, hsa-miR-145-5p and hsa-miR-204-5p were significantly altered. Moreover, as illustrated in Figs. 4 and 5, the top 10 target genes were TPM3, GXYLT1, SSH2, UBE2G1, CREB5, PARVA, TRAF3, IQSEC2, VANGL1 and IGF1R. Above all, these data showed that the regulatory network consisted of key miRNAs and genes may regulate biological process such as cell cycle, apoptosis and proliferation of human bladder cancer. Among the top ten target genes, tropomyosin 3 is a member of the tropomyosin family of actin-binding protein, which has been reported to relate to malignant transformation in BCa (36). VANGL planar cell polarity 1 (VANGL1), as an oncogene, is associated with many cancers. Park et al (37) revealed that miR-124 targeting VANGL could suppress colorectal cancer and Oh et al (38) reported that VANGL1 has correlation with tumor progression in human colorectal cancer. Insulin-like growth factor 1 receptor (IGF1R), with tyrosine kinase activity, binds insulin-like growth factor with a high affinity, playing a critical role in tumorigenesis and chemosensitivity (39,40). In the present study, IGF1R has high correlation with has-miR-378 family. Some studies indicated that miR-378 deficiency played a key role in the development of cardiac hypertrophy by targeting IGF1R through negatively regulated Ras signaling pathway (41). Parvin alpha, a member of the parvin family of actin-binding proteins, playing a role in cell adhesion, motility and survival, is connected with different cancers, such as colorectal cancer and lung cancer (42).

The four miRNAs screened out from the miRNA microarray and verified by another miRNA microarray indeed had a functional role in MIBC, which has been demonstrated by previous studies. Zhang et al (43) reported that miR-497 was downregulated in BCa tissues compared with normal bladder tissues and may represent a novel prognostic biomarker for the early detection of metastasis of bladder cancer. Du et al (44) also reported that miR-497 was decreased in plasma of bladder cancer patients compared with healthy patients and could be a promising novel circulating biomarkers in clinical detection of bladder cancer. Chiyomaru et al (45) found that miR-145 was a tumor suppressor and inhibited cell viability by targeting FSCN1 in BCa cells. In addition, Avgeris et al (46) suggested that miR-145 could act as a novel marker helpful for prediction of oncologic outcome for bladder cancer patients. However, there are no related reports about miR-29a-3p and miR-204-5p in BCa, therefore, the authors would like to investigate and confirm the two miRNAs using human BCa cell lines and mouse model in our next research article. However, predicting the miRNAs targets merely by bioinformatics analysis is not sufficient. Since the size of the MIBC samples used in the present study is small, these results may have many limitations. Thus, functional experiments should be performed strictly to verify the miRNAs and its targets in the further studies. As a result of that, the author group will select some of the significantly dysregulated miRNAs and perform verification experiments to confirm their targets and then figure out the functional role of miRNAs and the underlying mechanisms in MIBC.

Acknowledgements

The excellent technical assistance of Yuan Zhu, Shanshan Zhang and Danni Shan is gratefully acknowledged. The present study was supported in part by grants from the Zhongnan Hospital of Wuhan University Science, Technology and Innovation Seed Fund (grant no. cxpy20160031 and cxpy20160077) and the National Natural Science Foundation of China (grant no. 81300578). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References 1

Griffiths

Action on Bladder Cancer

Current perspectives in bladder cancer management

Int J Clin Pract674354482013

10.1111/ijcp.12075

23137019

Pastorelli

Saletta

Carpi

Campagna

del l'Osta

Schiarea

Vineis

Airoldi

Matullo

Proteome characterization of a human urothelial cell line resistant to the bladder carcinogen 4-aminobiphenyl

Proteome Sci562007

10.1186/1477-5956-5-6

17477866

1871571

Jin

Zhang

Sarnow

Kay

Biological basis for restriction of microRNA targets to the 3′ untranslated region in mammalian mRNAs

Nat Struct Mol Biol161441502009

10.1038/nsmb.1552

19182800

2713750

Bartel

Chen

Micromanagers of gene expression: The potentially widespread influence of metazoan microRNAs

Nat Rev Genet53964002004

10.1038/nrg1328

15143321

Hwang

Mendell

MicroRNAs in cell proliferation, cell death, and tumorigenesis

Br J Cancer947767802006

10.1038/sj.bjc.6603023

16495913

2361377

Nicoloso

Spizzo

Shimizu

Rossi

Calin

MicroRNAs-the micro steering wheel of tumour metastases

Nat Rev Cancer92933022009

10.1038/nrc2619

19262572

Adlakha

Saini

Brain microRNAs and insights into biological functions and therapeutic potential of brain enriched miRNA-128

Mol Cancer13332014

10.1186/1476-4598-13-33

24555688

3936914

Lin

Kong

Huang

Kim

Etheridge

Luo

Ding

Wang

Current State of circulating MicroRNAs as cancer biomarkers

Clin Chem61113811552015

10.1373/clinchem.2015.241190

26319452

Chawla

Iyer

Soodan

Sharma

Khurana

Priyadarshni

Role of miRNA in cancer diagnosis, prognosis, therapy and regulation of its expression by Epstein-Barr virus and human papillomaviruses: With special reference to oral cancer

Oral Oncol517317372015

10.1016/j.oraloncology.2015.05.008

26093389

Andrew

Marsit

Schned

Seigne

Kelsey

Moore

Perreard

Karagas

Sempere

Expression of tumor suppressive microRNA-34a is associated with a reduced risk of bladder cancer recurrence

Int J Cancer137115811662015

10.1002/ijc.29413

25556547

4485975

Bandres

Agirre

Ramirez

Zarate

Garcia-Foncillas

MicroRNAs as cancer players: Potential clinical and biological effects

DNA Cell Biol262732822007

10.1089/dna.2006.0544

17504023

Catto

Miah

Owen

Bryant

Myers

Dudziec

Larré

Milo

Rehman

Rosario

Distinct microRNA alterations characterize high- and low-grade bladder cancer

Cancer Res69847284812009

10.1158/0008-5472.CAN-09-0744

19843843

2871298

Pignot

Cizeron-Clairac

Vacher

Susini

Tozlu

Vieillefond

Zerbib

Lidereau

Debre

Amsellem-Ouazana

Bieche

microRNA expression profile in a large series of bladder tumors: Identification of a 3-miRNA signature associated with aggressiveness of muscle-invasive bladder cancer

Int J Cancer132247924912013

10.1002/ijc.27949

23169479

Gottardo

Liu

Ferracin

Calin

Fassan

Bassi

Sevignani

Byrne

Negrini

Pagano

Micro-RNA profiling in kidney and bladder cancers

Urol Oncol253873922007

10.1016/j.urolonc.2007.01.019

17826655

Hidaka

Seki

Yoshino

Yamasaki

Yamada

Nohata

Fuse

Nakagawa

Enokida

Tumor suppressive microRNA-1285 regulates novel molecular targets: Aberrant expression and functional significance in renal cell carcinoma

Oncotarget344572012

10.18632/oncotarget.417

22294552

3292891

Yoshino

Seki

Itesako

Chiyomaru

Nakagawa

Enokida

Aberrant expression of microRNAs in bladder cancer

Nat Rev Urol103964042013

10.1038/nrurol.2013.113

23712207

Cao

Meng

Liu

Wang

Qian

Cao

Guan

Dan

Xiao

Wang

Decreased TRPM7 inhibits activities and induces apoptosis of bladder cancer cells via ERK1/2 pathway

Oncotarget772941729602016

27662662

5341955

Chen

Wang

Luo

Wang

Xie

Jiang

Xiao

Qian

Wang

Downregulation of LAPTM5 suppresses cell proliferation and viability inducing cell cycle arrest at G0/G1 phase of bladder cancer cells

Int J Oncol502632712017

10.3892/ijo.2016.3788

27922670

Wang

Cao

Wang

Qian

Dan

Jiang

Xiao

Wang

Simvastatin induces cell cycle arrest and inhibits proliferation of bladder cancer cells via PPARgamma signalling pathway

Sci Rep6357832016

10.1038/srep35783

27779188

5078845

Tusher

Tibshirani

Chu

Significance analysis of microarrays applied to the ionizing radiation response

Proc Natl Acad Sci USA9851165121

2001

10.1073/pnas.091062498

11309499

33173

Chu

Seo

Narasimhan

Tibshirani

Tusher

SAM: Significance Analysis of Microarrays

Users Guide and technical documenthttp://statweb.stanford.edu/~tibs/SAM/sam.pdf

Storey

A direct approach to false discovery rates

J Roy Stat Soc Ser644794982002

10.1111/1467-9868.00346

Pan

The dual regulatory role of miR-204 in cancer

Tumour Biol3711667116772016

10.1007/s13277-016-5144-5

27438705

5080331

Yang

Xiong

Cao

Zheng

You

Zhang

Zhao

miR-497 expression, function and clinical application in cancer

Oncotarget755900559112016

10.18632/oncotarget.10152

27344185

5342460

Wang

Zhang

Ren

The role of miRNA-29 family in cancer

Eur J Cell Biol921231282013

10.1016/j.ejcb.2012.11.004

23357522

Wang

Fan

Cai

Meng

MiR-124 exerts tumor suppressive functions on the cell proliferation, motility and angiogenesis of bladder cancer by fine-tuning UHRF1

FEBS J282437643882015

10.1111/febs.13502

26310391

Shaham

Binder

Gefen

Borkhardt

Izraeli

miR-125 in normal and malignant hematopoiesis

Leukemia26201120182012

10.1038/leu.2012.90

22456625

Zhu

Chen

Mao

Lin

Luo

MicroRNA-490-5p inhibits proliferation of bladder cancer by targeting c-Fos

Biochem Biophys Res Commun4419769812013

10.1016/j.bbrc.2013.11.006

24220339

Yonemori

Seki

Yoshino

Matsushita

Miyamoto

Nakagawa

Enokida

Dual tumor-suppressors miR-139-5p and miR-139-3p targeting matrix metalloprotease 11 in bladder cancer

Cancer Sci107123312422016

10.1111/cas.13002

27355528

5021030

Tao

Qian

Yin

Zhang

microRNA-133 inhibits cell proliferation, migration and invasion in prostate cancer cells by targeting the epidermal growth factor receptor

Oncol Rep27196719752012

22407299

Kou

Gao

Shi

Wang

Guo

miR-145 inhibits invasion of bladder cancer cells by targeting PAK1

Urol Oncol328468542014

10.1016/j.urolonc.2014.01.003

24954107

Yoshino

Chiyomaru

Enokida

Kawakami

Tatarano

Nishiyama

Nohata

Seki

Nakagawa

The tumour-suppressive function of miR-1 and miR-133a targeting TAGLN2 in bladder cancer

Br J Cancer1048088182011

10.1038/bjc.2011.23

21304530

3048214

Ambros

The functions of animal microRNAs

Nature4313503552004

10.1038/nature02871

15372042

Pierzynski

Hildebrandt

Kamat

Lin

Dinney

Genetic Variants in the Wnt/β-Catenin signaling pathway as indicators of bladder cancer risk

J Urol194177117762015

10.1016/j.juro.2015.07.032

26173102

5087323

Lee

Lim

Choi

Kim

Moon

Interleukin-28A triggers wound healing migration of bladder cancer cells via NF-κB-mediated MMP-9 expression inducing the MAPK pathway

Cell Signal24173417422012

10.1016/j.cellsig.2012.04.013

22560878

Pawlak

McGarvey

Nguyen

Tomaszewski

Puthiyaveettil

Malkowicz

Helfman

Alterations in tropomyosin isoform expression in human transitional cell carcinoma of the urinary bladder

Int J Cancer1103683732004

10.1002/ijc.20151

15095301

Park

Kim

Yoon

Bae

Choi

Jung

Kim

KITENIN-targeting microRNA-124 suppresses colorectal cancer cell motility and tumorigenesis

Mol Ther22165316642014

10.1038/mt.2014.105

24909917

4435482

Park

Kim

Park

Oak

Myung

Cho

Lee

Kim

Joo

Impact of KITENIN on tumor angiogenesis and lymphangiogenesis in colorectal cancer

Oncol Rep352532602016

10.3892/or.2015.4337

26496979

XMZhu

Bai

Wang

Chen

Yang

Jia

Epigenetic silencing of miR-375 induces trastuzumab resistance in HER2-positive breast cancer by targeting IGF1R

BMC Cancer141342014

10.1186/1471-2407-14-134

24571711

3974046

Taliaferro-Smith

Oberlick

Liu

McGlothen

Alcaide

Tobin

Donnelly

Commander

Kline

Nagaraju

FAK activation is required for IGF1R-mediated regulation of EMT, migration, and invasion in mesenchymal triple negative breast cancer cells

Oncotarget6475747722015

10.18632/oncotarget.3023

25749031

4467113

Nagalingam

Sundaresan

Gupta

Geenen

Solaro

Gupta

A cardiac-enriched microRNA, miR-378, blocks cardiac hypertrophy by targeting Ras signaling

J Biol Chem28811216112322013

10.1074/jbc.M112.442384

23447532

3630862

Huang

Pan

Chang

Hong

Chen

PARVA promotes metastasis by modulating ILK signalling pathway in lung adenocarcinoma

PLoS One10e01185302015

10.1371/journal.pone.0118530

25738875

4349696

Zhang

Gong

Zhan

Man

Kong

MicroRNA-497 inhibits the proliferation, migration and invasion of human bladder transitional cell carcinoma cells by targeting E2F3

Oncol Rep36129313002016

10.3892/or.2016.4923

27430325

Shi

Yuan

Chu

Qin

Yin

Zhang

Wang

Circulating miR-497 and miR-663b in plasma are potential novel biomarkers for bladder cancer

Sci Rep5104372015

10.1038/srep10437

26014226

4444850

Chiyomaru

Enokida

Tatarano

Kawahara

Uchida

Nishiyama

Fujimura

Kikkawa

Seki

Nakagawa

miR-145 and miR-133a function as tumour suppressors and directly regulate FSCN1 expression in bladder cancer

Br J Cancer1028838912010

10.1038/sj.bjc.6605570

20160723

2833258

Avgeris

Mavridis

Tokas

Stravodimos

Fragoulis

Scorilas

Uncovering the clinical utility of miR-143, miR-145 and miR-224 for predicting the survival of bladder cancer patients following treatment

Carcinogenesis365285372015

10.1093/carcin/bgv024

25804644

Figure 1.

Differential expression of miRNAs between normal bladder and MIBC tissues. Heatmap of 104 microRNAs differentially expressed between three normal bladder and three MIBC tissues. MIBC, muscle invasive bladder cancer.

Figure 2.

Significantly changed GO/pathways of predicted target genes of dysregulated miRNAs. (A) Significantly changed GOs of predicted target genes. The y-axis shows GO category and the x-axis shows -log10 (P-value). (B) Significantly changed pathways of predicted target genes. The γ-axis shows significantly changed pathways. GO, gene ontology.

Figure 3.

Pathway network (Path-net). Significantly changed pathways were connected in a Path-net to show the interaction network among these pathways. Each pathway in the network was measured by counting the upstream and downstream pathways. The blue circle represents pathways involving upregulated miRNAs, while the yellow circle represents pathways involving both upregulated and downregulated miRNAs. The size of the circle represents the degree value and the lines show the interaction between pathways. A higher degree of pathway indicates that it plays a more important role in the signaling network. miRNAs, microRNAs.

Figure 4.

miRNAs-gene-network. According to the interactions between miRNAs and the intersected target genes, miRNAs-gene-network was constructed to illustrate the key regulatory functions of the identified miRNAs and their target genes. The blue circles represent genes, while blue square nodes represent downregulated miRNAs. The size of the circle or square node represents the degree value. A higher degree of gene/miRNAs indicates that it plays a more important role in the signaling network. miRNA/miR, microRNA.

Figure 5.

miRNAs-GO-network. The miRNAs-GO-network was generated according to the relationship of significant biological functions and miRNAs. The yellow and blue circles represent GOs, red square nodes represent upregulated miRNAs, and blue square nodes represent downregulated miRNAs. The size of the circle or square node represents the degree value. A higher degree of GO/miRNAs indicates that it plays a more important role in the signaling network. GO, gene ontology; miRNA/miR, microRNA.

Figure 6.

Validation of candidate miRNAs. Among the top rated nine miRNAs screened from our miRNA microarray, four miRNAs involving (A) hsa-miR-145-5p, (B) hsa-miR-497-5p, (C) hsa-miR-29a-3p and (D) hsa-miR-204-5p were also significantly altered in GSE40355. miRNA, microRNA.

Table I.

Information of the patients and donors.

Characteristics	MIBC patients	Donors
Number	3	3
Age, years (mean ± SD)	62±1.581	37±2.327
Gender	Male	Male
BCa stage	Stage II	–
Surgical method	Radical resection	–

Stage II denotes that BCa enters the muscle layer of the bladder. SD, standard deviation; MIBC, muscle invasive bladder cancer; BCa, bladder cancer.

Table II.

Differently expressed miRNAs in human MIBC tissues.

MIBC tissues vs. normal bladder tissues. (P<0.05, Fold-change >1.5)

miRNA	Change	P-value	Feature^a	Rank^b
hsa-miR-4786-5p	−10.205618	0.00124	Down	1
hsa-miR-490-3p	−35.773813	0.001358	Down	2
hsa-miR-3617-5p	−42.647285	0.001476	Down	3
hsa-miR-490-5p	−75.331847	0.001594	Down	4
hsa-miR-139-3p	−32.552806	0.001712	Down	5
hsa-miR-133b	−121.234878	0.00183	Down	6
hsa-miR-145-3p	−9.840205	0.001948	Down	7
hsa-miR-29b-1-5p	−10.581222	0.002066	Down	8
hsa-miR-155-5p	8.720112	0.002184	UP	9
hsa-miR-1	−40.712728	0.002302	Down	10
hsa-miR-548q	−5.779384	0.00242	Down	11
hsa-miR-133a-5p	−17.075477	0.002774	Down	12
hsa-miR-146b-5p	4.140864	0.002893	UP	13
hsa-miR-28-3p	−3.834537	0.003129	Down	14
hsa-miR-30c-2-3p	−11.293849	0.003247	Down	15
hsa-miR-143-5p	−15.820855	0.003365	Down	16
hsa-miR-6511b-3p	−6.432579	0.003601	Down	17
hsa-miR-320e	−3.630028	0.003719	Down	18
hsa-miR-99a-3p	−2.325859	0.003837	Down	19
hsa-miR-30a-3p	−18.900354	0.003955	Down	20
hsa-miR-133a-3p	−63.82892	0.004073	Down	21
hsa-miR-4510	−2.805811	0.004545	Down	22
hsa-miR-139-5p	−12.808018	0.004782	Down	23
hsa-miR-4324	−19.035251	0.005136	Down	24
hsa-miR-145-5p	−3.499059	0.005372	Down	25
hsa-miR-125b-1-3p	−13.117689	0.005608	Down	26
hsa-miR-143-3p	−3.894089	0.005962	Down	27
hsa-miR-193a-5p	−4.050292	0.006198	Down	28
hsa-miR-5684	−4.681531	0.007261	Down	29
hsa-miR-328-3p	−8.31845	0.007379	Down	30
hsa-miR-4429	−2.76178	0.007497	Down	31
hsa-miR-1287-5p	−7.357196	0.007615	Down	32
hsa-miR-3605-5p	−2.739437	0.007969	Down	33
hsa-miR-338-5p	−4.852428	0.008323	Down	34
hsa-miR-193a-3p	−9.317892	0.00856	Down	35
hsa-miR-4257	−3.696435	0.008796	Down	36
hsa-miR-6507-5p	−6.813232	0.010331	Down	37
hsa-miR-6722-3p	−2.904131	0.011275	Down	38
hsa-miR-497-5p	−3.273151	0.011865	Down	39
hsa-miR-125b-2-3p	−7.613095	0.012574	Down	40
hsa-miR-29b-2-5p	−2.995467	0.013046	Down	41
hsa-miR-6511a-3p	−8.396761	0.013282	Down	42
hsa-miR-628-3p	−4.613832	0.013636	Down	43
hsa-miR-378b	−6.024561	0.014581	Down	44
hsa-miR-664a-5p	−6.119249	0.014699	Down	45
hsa-miR-204-5p	−6.006885	0.015525	Down	46
hsa-miR-7641	−4.116459	0.015762	Down	47
hsa-miR-320d	−2.149376	0.016706	Down	48
hsa-miR-3656	−2.093619	0.017414	Down	49
hsa-miR-1273g-3p	−2.156204	0.017769	Down	50
hsa-miR-378g	−7.277976	0.018359	Down	51
hsa-miR-1225-5p	−2.552195	0.019067	Down	52
hsa-miR-3156-5p	−6.520017	0.019421	Down	53
hsa-miR-383-5p	−3.70084	0.019658	Down	54
hsa-miR-3064-3p	−2.07381	0.020012	Down	55
hsa-miR-378d	−5.5391	0.02013	Down	56
hsa-miR-23b-5p	−9.113128	0.020366	Down	57
hsa-miR-3195	−10.384316	0.020484	Down	58
hsa-miR-100-5p	−3.57042	0.020956	Down	59
hsa-miR-378e	−14.003579	0.021547	Down	60
hsa-miR-4649-5p	−2.337674	0.022019	Down	61
hsa-miR-3622b-5p	−3.04368	0.022373	Down	62
hsa-miR-6840-3p	−6.817266	0.024616	Down	63
hsa-miR-4322	−6.233632	0.024852	Down	64
hsa-miR-574-5p	−5.490031	0.025207	Down	65
hsa-miR-4770	−4.247774	0.025561	Down	66
hsa-miR-513a-5p	−3.907964	0.025679	Down	67
hsa-miR-124-3p	−1.707013	0.025797	Down	68
hsa-miR-6819-5p	−2.385821	0.026151	Down	69
hsa-miR-339-3p	−5.134485	0.027332	Down	70
hsa-miR-26b-3p	−3.129563	0.02745	Down	71
hsa-miR-29c-5p	−7.931159	0.027922	Down	72
hsa-miR-29a-3p	−2.906434	0.02804	Down	73
hsa-miR-1227-5p	−2.287321	0.028276	Down	74
hsa-miR-4327	−2.710984	0.028512	Down	75
hsa-miR-99a-5p	−4.947917	0.02863	Down	76
hsa-miR-1296-5p	−9.401257	0.029103	Down	77
hsa-miR-378i	−4.649838	0.030283	Down	78
hsa-miR-6889-5p	−4.810854	0.030401	Down	79
hsa-miR-186-5p	−3.276383	0.031464	Down	80
hsa-miR-6824-5p	−8.197661	0.031936	Down	81
hsa-miR-195-3p	−4.796311	0.034416	Down	82
hsa-miR-3622a-5p	−4.145204	0.034888	Down	83
hsa-miR-5572	−6.053895	0.03536	Down	84
hsa-miR-6776-5p	−3.050977	0.035832	Down	85
hsa-miR-6790-5p	−2.868339	0.036895	Down	86
hsa-miR-6740-5p	−2.087867	0.037721	Down	87
hsa-miR-4269	−8.226931	0.038312	Down	88
hsa-miR-125b-5p	−2.694152	0.039138	Down	89
hsa-miR-3188	−8.008159	0.039256	Down	90
hsa-miR-6787-5p	−2.714687	0.03961	Down	91
hsa-miR-99b-5p	−1.811603	0.039847	Down	92
hsa-miR-30a-5p	−3.787576	0.040673	Down	93
hsa-miR-6765-5p	−2.365181	0.041027	Down	94
hsa-miR-1909-3p	−7.648533	0.041972	Down	95
hsa-miR-4507	−2.786848	0.042326	Down	96
hsa-miR-28-5p	−2.529154	0.042916	Down	97
hsa-miR-6892-5p	−2.081233	0.046458	Down	98
hsa-miR-4433b-3p	−1.923165	0.046694	Down	99
hsa-miR-6737-5p	−3.166103	0.04693	Down	100
hsa-miR-324-3p	−3.939764	0.048465	Down	101
hsa-miR-211-3p	−3.140159	0.048819	Down	102
hsa-miR-6798-5p	−2.080154	0.049174	Down	103
hsa-miR-4299	−4.93377	0.049882	Down	104

Changes of miRNA expression in MIBC compared with normal bladder.

Ranking method is in descending order according to statistical significance (P-value) between the two groups. miRNA/miR, microRNA; MIBC, muscle invasive bladder cancer; Up, upregulated in MIBC; Down, downregulated in MIBC.

Table III.

The top 10 miRNAs with high degrees of miRNAs-gene-networks and miRNAs-GO-networks.

Rank	miRNAs	miRNA-gene-networks degree^a	Feature^b	miRNAs	miRNA-GO-networks degree^a	Feature^b
1	hsa-miR-497-5p	123	Down	hsa-miR-497-5p	637	Down
2	hsa-miR-29a-3p	88	Down	hsa-miR-29a-3p	518	Down
3	hsa-miR-124-3p	84	Down	hsa-miR-124-3p	516	Down
4	hsa-miR-4269	64	Down	hsa-miR-204-5p	494	Down
5	hsa-miR-145-5p	63	Down	hsa-miR-4269	458	Down
6	hsa-miR-204-5p	62	Down	hsa-miR-145-5p	430	Down
7	hsa-miR-4510	61	Down	hsa-miR-6840-3p	406	Down
8	hsa-miR-125b-5p	57	Down	hsa-miR-4510	404	Down
9	hsa-miR-6840-3p	50	Down	hsa-miR-6722-3p	387	Down
10	hsa-miR-6722-3p	47	Down	hsa-miR-1	383	Down

Changes of miRNA expression in MIBC compared with normal bladder.

Total number of connections around the miRNA. miRNA/miR, microRNA; GO, gene ontology; Down, downregulated in muscle invasive bladder cancer.

Table IV.

The alteration of top rated nine miRNAs in GSE76211 and GSE40355.

	GSE76211			GSE40355

miRNAs	Fold-change	P-value	Feature^a	Fold-change	P-value	Feature^a
hsa-miR-497-5p	−3.27	0.011865	Down	−7.94	5.94×10⁻⁵	Down
hsa-miR-29a-3p	−2.90	0.02804	Down	−3.61	0.0003	Down
hsa-miR-124-3p	−1.70	0.025797	Down	–	–	Down
hsa-miR-4269	−8.22	0.038312	Down	–	–	Down
hsa-miR-145-5p	−3.49	0.005372	Down	−52.69	4.39×10⁻⁸	Down
hsa-miR-204-5p	−6.00	0.015525	Down	−463.33	1.89×10⁻¹⁰	Down
hsa-miR-4510	−2.80	0.004545	Down	–	–	Down
hsa-miR-6840-3p	−6.81	0.024616	Down	–	–	Down
hsa-miR-6722-3p	−2.90	0.011275	Down	–	–	Down

Changes of miRNA expression in MIBC compared with normal bladder. miRNA/miR, microRNA; MIBC, muscle invasive bladder cancer; Down, downregulated in MIBC.