Urothelial carcinoma of the bladder (UCC) is a common disease that arises by at least two different molecular pathways. The biology of UCC is incompletely understood, making the management of this disease difficult. Recent evidence implicates a regulatory role for microRNA in cancer. We hypothesized that altered microRNA expression contributes to UCC carcinogenesis. To test this hypothesis, we examined the expression of 322 microRNAs and their processing machinery in 78 normal and malignant urothelial samples using real-time rtPCR. Genes targeted by differentially expressed microRNA were investigated using real-time quantification and microRNA knockdown. We also examined the role of aberrant DNA hypermethylation in microRNA downregulation. We found that altered microRNA expression is common in UCC and occurs early in tumorogenesis. In normal urothelium from patients with UCC, 11% of microRNAs had altered expression when compared with disease-free controls. This was associated with upregulation of Dicer, Drosha, and Exportin 5. In UCC, microRNA alterations occur in a tumor phenotype–specific manner and can predict disease progression. High-grade UCC were characterized by microRNA upregulation, including microRNA-21 that suppresses p53 function. In low-grade UCC, there was downregulation of many microRNA molecules. In particular, loss of microRNAs-99a/100 leads to upregulation of FGFR3 before its mutation. Promoter hypermethylation is partly responsible for microRNA downregulation. In conclusion, distinct microRNA alterations characterize UCC and target genes in a pathway-specific manner. These data reveal new insights into the disease biology and have implications regarding tumor diagnosis, prognosis and therapy. [Cancer Res 2009;69(21):8472–81]

Bladder cancer is the fifth commonest malignancy in the United States with an incidence of 67,160 new cases and 13,750 deaths in 2007 (1). The majority of tumors are urothelial cell carcinoma (UCC). Clinical and molecular evidence suggests there are at least two distinct varieties of this tumor. Most UCC belong to a low-grade pathway characterized by FGFR3 mutation, chromosome 9 loss, and an indolent clinical phenotype. Around one-third of UCC are high grade in differentiation and arise as lesions initially confined to the bladder mucosa (nonmuscle invasive, NMI). Progression to muscle invasion occurs in around 50% of high-grade lesions and is associated with an ominous prognosis despite radical treatment (2, 3). Although numerous genetic and epigenetic alterations are observed in high-grade UCC, loss of p53 function seems most critical (4, 5).

MicroRNAs (miR) are short noncoding RNA molecules that posttranscriptionally modulate protein expression. They are transcribed as hairpin pri-miRs and processed into pre-miRs by Drosha, an RNase III endonuclease complexed with DGCR8. Pre-MiRs are exported into the cytoplasm by Exportin 5 before cleavage by Dicer into mature miRs. MiRs are directed to mRNAs with a complementary seed sequence to their first 1 to 8 nucleotides (6). Alterations in mRNA expression seem important for carcinogenesis (7, 8). Their profile may be used to identify key tumorogenic pathways (9) or clinical outcome (10, 11).

To date, few data detail miR in UCC (9, 12, 13). We hypothesized that altered miR expression occurs in UCC and contributes to carcinogenesis. To test this hypothesis, we investigated the expression of numerous miRs and their processing molecules in urothelial tissues. Our studies reveal that UCC are characterized by widespread alterations in miR expression. Changes in expression occur early in the tumorogenic pathway, in a phenotype-specific manner, and are associated with altered expression of their processing machinery. Of particular note, low-grade tumors are characterized by miR-mediated FGFR3 upregulation before its mutation.

Patients, samples, and cell lines

We studied 72 urothelial samples and six cell lines (Table 1). Tumors were classified using the 2004 WHO/ISUP criteria, and selected from three cancer groups: (a) low-grade NMI, (b) high-grade NMI, and (c) invasive UCC. Twenty normal urothelial samples were obtained from 10 patients with UCC [NU(UCC) taken distant to any tumor] and 10 disease-free controls [NU(Controls)]. Tissues were immediately frozen in liquid nitrogen and histologic confirmation was obtained before use. Patients were treated according to tumor stage and grade following endoscopic resection (14). Adjuvant intravesical chemotherapy was used for low-grade and maintenance intravesical Bacillus Calmette-Guerin (BCG) for high-grade NMI tumors. Patients underwent surveillance stratified by their disease. Radical cystectomy with pelvic lymphadenectomy was used for invasive or BCG-refractory high-grade NMI tumors (4, 14). We analyzed UCC cell lines representing the three tumor groups (RT4, RT112, and EJ/T24, respectively, purchased from American Type Culture Collection) grown in Dulbecco's medium with 10% FCS, and three normal human urothelial (NHU) cell lines maintained in keratinocyte serum–free medium containing bovine pituitary extract, epidermal growth factor (Invitrogen), and cholera toxin. Nonimmortalized NHU cells were derived from histologically confirmed normal urothelium obtained from patients without a history of UCC, using standard methods (15). DNA Methyltransferase inhibition experiments were performed in quadruplet by adding 2 μmol/L 5-azacytidine to the media for 5 to 7 d before harvesting.

Table 1.

Details of patients and tissue samples analyzed in this study

UCCNormal urothelium
Low-grade NMIHigh-grade NMIInvasiveNon-UCCUCC patients
Total 22 12 18 10 10 
Gender 
    Male 14 12 14 10 
    Female 
Age 
    Mean 75 71.9 71.4 67 66.4 
    Range 56–83 55–81 46–90 54–81 46–90 
Source of urothelium* 
    TURP     
    RP    
    TURBT     
    RC     
Grade 
    1 11  1 
    2 11  4 
    3 12 17  5 
Stage 
    pTa 22  2 
    pTis   
    pT1   
    pT2-4 18  8 
Smoker 
    Current   
    Ex 13 10 12   
    Never   
Recurrence 
    Yes 14   
    No 13   
Progression 
    Yes 11   
    No 18   
Follow up (mo) 
    Mean 43.4 57.0 23.5   
    Range 4–93 15–93 2–84   
UCCNormal urothelium
Low-grade NMIHigh-grade NMIInvasiveNon-UCCUCC patients
Total 22 12 18 10 10 
Gender 
    Male 14 12 14 10 
    Female 
Age 
    Mean 75 71.9 71.4 67 66.4 
    Range 56–83 55–81 46–90 54–81 46–90 
Source of urothelium* 
    TURP     
    RP    
    TURBT     
    RC     
Grade 
    1 11  1 
    2 11  4 
    3 12 17  5 
Stage 
    pTa 22  2 
    pTis   
    pT1   
    pT2-4 18  8 
Smoker 
    Current   
    Ex 13 10 12   
    Never   
Recurrence 
    Yes 14   
    No 13   
Progression 
    Yes 11   
    No 18   
Follow up (mo) 
    Mean 43.4 57.0 23.5   
    Range 4–93 15–93 2–84   

*RC, radical cystectomy; RP, radical prostatectomy; TURP/BT, transurethral resection of prostate/bladder tumor.

Details of associated tumour.

Ex-smoking status defined as stopping >1 y before tumor.

MiR expression profiling

For both normal and malignant samples, ten 10 μmol/L sections were microdissected to obtain >90% pure cell populations. Enriched small and total RNA were extracted using the mirVana kit (Ambion) according to manufacturer's protocol. RNA concentrations were measured using a 2100 Bioanalyzer (Agilent, Cheshire). The expression of 354 mature miRs and 3 small nucleolar RNA molecules was determined using real-time PCR with preprinted microfluidic cards (Human miR v1.0, Applied Biosystems; ref. 16). Initially, reverse transcription using stem loop primers was performed with 50 ng RNA, MultiScribe Reverse Transcriptase (Applied Biosystems), RNase inhibitor, 100 nm deoxynucleotide triphosphates (dNTP) and nuclease-free water. To generate 357 miR and reference rtPCR products, 8 multiplex pools were used per sample [each pool containing ingredients and primers for 48 mature miRs: total volume 10 μL (Human Multiplex RT Pool 1-8 v1.0, Applied Biosystems)]. The multiplexed RT assays were cooled to 16°C for 30 min, heated to 42°C for 30 min, and 86°C for 5 min before diluting 62.5-fold in nuclease-free water. From each diluted assay, 50 μL were mixed equally with 2× Taqman Universal PCR MasterMix and loaded separately into the eight reservoirs of a microfluidic card (Human miR v1.0, Applied Biosystems). Following careful centrifugation (2 min at 331 grams) and preparation (well sealing and reservoir removal), the loaded card was analyzed on an ABI 7900HT real-time PCR system using the recommended PCR conditions. Each microfluidic card contained MGB-labeled probes specific to 354 mature miRs, 11 duplicate assays, 2 empty wells, and 17 replicates of 3 endogenous small nucleolar RNA molecules for relative quantification (RNU6B, RNU44, and RNU48). MiR concentrations were calculated from the cycle number that each reaction crossed an arbitrarily threshold (Ct = 0.2). The amplification plots were checked manually (SDS 2.2.1, Applied Biosystems) to confirm the Ct value corresponded with the midpoint of logarithmic amplification. To test the reproducibility of this real-time quantification, we analyzed each cell line and 10 tissue samples in duplicate.

Identification of genes targeted by aberrantly expressed miR

We identified potential protein targets of miRs using TargetScan7

(version 4.2; refs. 17, 18). Target genes were ranked according to context score and the average was used for those with multiple targeting miRs. mRNA transcripts of target and miR-processing genes were measured in those samples with sufficient material [n = 71 (93%)]. cDNA synthesis from 100 ng whole RNA was performed using random primers, RT buffer, dNTP (100 mmol/L), RNase inhibitor, and MultiScribe Reverse Transcriptase (cDNA Reverse Transcription kit, Applied Biosystems). Real-time quantified PCR using 2 μL of synthesized cDNA, gene-specific primers with FAM-TAMRA–labeled probes (Supplementary Table S1), water, and 2× Taqman Universal PCR MasterMix (Applied Biosystems) in a total volume of 10 μL was performed on the ABI 7900HT system according to manufacturers guidelines. Relative mRNA quantification was measured with respect to the mean of glyceraldehyde-3-phosphate dehydrogenase and β-actin (and PCNA as a proliferation reference gene).

miR knockdown in NHU cells

MiRs-99a/100 targeting of FGFR3 was investigated using protein expression following miR knockdown and a Luciferase reporter assay. All experiments were replicated between three and six times using nonimmortalized NHU cell lines at 70% confluence reverse transfected with anti–miR molecules specific to miR-99a, miR-100, and a negative control (scrambled sequence). To determine target protein expression, NHU cells were transfected with 50 nmol/L anti-miR or control sequences using siPORT NeoFX transfection reagent (Ambion). The mixture was dispensed into a six-well dish, seeded at 2 × 105 cells per well, and incubated for 48 h. Transfection was determined by real-time-rtPCR (Applied Biosystems) and FGFR3 expression measured by Western Blotting using anti-FGFR3 primary antibody (1:1,000; Cell Signaling Technologies, Inc.). Cells were lysed in radioimmunoprecipitation assay buffer (20 mmol/L Tris·HCl, 135 mmol/L NaCl, 10% glycerol, 1% Igepal, 0.1% SDS, 0.5% deoxycholic acid, 2 mmol/L EDTA) containing protease and phosphatase inhibitors (Complete EDTA-free protease inhibitor cocktail and PhosSTOP phosphatase inhibitor cocktail; Roche), and protein content was quantified using the DC-protein assay reagent (Bio-Rad). Protein lysates (50 μg) were loaded onto 8% gels, fractionated, and electroblotted onto nitrocellulose membranes. After blocking with 5% nonfat milk powder and 0.1% Tween, the membranes were incubated overnight with anti-FGFR3 primary antibody (1:1,000; Cell Signaling Technologies, Inc.) at 4°C, washed, and incubated at room temperature for 1 h with goat anti-rabbit peroxidase-labeled secondary antibody (1:1,000; Cell Signaling Technologies, Inc.). β-Actin expression was also measured as a loading control (Sigma, 1:2,000). Immune complexes were quantified using ImageJ for band densitometry. To validate direct targeting of FGFR3, NHU cells were cotransfected with either anti–miR-99a, anti–miR-100, or negative control, and the pSSG_3 untranslated region (UTR) plasmid containing the luc2P Luciferase reporter and either the partial 3′UTR for FGFR3 (from 27–1,762 bp downstream of the gene) or a control (scrambled sequence; SwitchGear Genomics). Cotransfection was performed in opaque 96-well plates using 40% to 80% confluent cells and FuGene transfection reagent (Roche) and Opti-MEM serum–free media. After incubation for 48 h, luciferase activity was determined following the addition of Promega Steady-Glo Luciferase Assay Reagent (Promega) in a luminometer. Negative controls using scrambled anti-miR sequences or a scrambled 3′ UTR sequence were included.

Generation of FGFR3 mutant and wild-type paired cell lines

Telomerase-immortalized NHU cells were used to investigate miR expression following FGFR3 mutation (19). Site-directed mutagenesis on FGFR3 IIIb cDNA created the S249C mutation. The presence of this mutation was verified by sequencing. Wild-type and mutant FGFR3 were cloned into a retroviral expression vector (pFB; Stratagene) containing a hygromycin resistance cassette (20). The expression vectors were transfected into Phoenix A cells using siPORT XP-1 (Ambion). Telomerase-immortalized NHU were incubated for 8 h with retroviral supernatant containing 8 μg/mL polybrene and selected with hygromycin 48 h after transduction. FGF1 cell stimulation was also performed using 20 ng/mL recombinant human FGF1 and 10 μg/mL heparin (R&D Systems Europe Ltd.) after 3 h of incubation in depleted medium (KGM2 without added growth factors).

Statistical methods

Relative miR concentrations were calculated with respect to the median of three reference RNA molecules (ΔCt = Ct miR − Ct median control). The median was chosen after analysis revealed variation in concentrations between samples (data not shown). Expression fold changes were computed using −ΔΔCt calculations (21). Median data centering was performed before analysis with Significance of Analysis of Microarray software (22) and BRB-ArrayTools8

(version 3.7) developed by Dr. Richard Simon and BRB-ArrayTools Development Team. Class-specific miR signatures were established using Prediction of Microarray software (23). MiR expression and tumor outcome were investigated using the log-rank test and plotted by the Kaplan-Meier method within SPSS (version 14.0, SPSS, Inc.). Tumor progression was defined as the presence of pathologic, radiological, or clinical evidence of an increase in tumor stage and measured from the time of surgery to the time of proven event. Changes in miR expression and statistical significance were illustrated using Volcano plots. Significant expression difference was defined as both a P value of < 0.05 and a false discovery rate of <0.05. Hierarchical clustering was performed after removal of those miRs with a low detection frequencies in all groups (<20%: defined as “noise”) using Cluster 3.0 and visualized in Tree view (Eisen Lab9). Correlation between variables was assessed using Pearson's correlation coefficient within SPSS.

MiR expression in normal urothelium

Following assay optimization (Supplementary Results), we compared the expression of 322 miRs in normal urothelium from patients with UCC and disease-free controls. Thirty-six miRs were differentially expressed between the two groups (P < 0.05 and FDR 0.0; Supplementary Table S2) and were always upregulated in the NU(UCC). There was clustering of differential expression according to RNA family, e.g., let-7a/b/c/d/e/g, miRs-382/487b, miRs-17/93, miRs-181b/d, miRs-34a/c, and chromosomal location, e.g., 9q22.32 (let-7a/d), 14q32.31 (miRs-376a/487b/382/412), and 19q13.41 (miRs-520a/b/c/d/e/518a/519e). Unsupervised hierarchical clustering revealed two branches corresponding to the two groups (Supplementary Fig. S1). Three miRs were detected only in NU(UCC) [miR-569 (detected in 50% of samples), miR-633 (40%), miR-507(30%)] and two by only the NU(controls) [miR-412 (50%), miR551a (30%); data not shown].

Differential miR expression in malignant urothelium

We compared malignant miR expression with that from disease-free normal urothelium (Fig. 1). The annotated raw data can be downloaded from our website.10

Twelve miRs were expressed only in the malignant tissues and five to seven of these could identify cancer with a sensitivity of 90% to 100% and a specificity of 100% to 80% (Table 2A), suggesting roles as diagnostic biomarkers. Quantitative analyses revealed that 16 miRs were differentially expressed in the malignant and normal urothelia (P < 0.05 and FDR<0.05; Table 2B). Thirteen of these were upregulated in cancer and seven of these were also aberrantly expressed in the normal urothelium from UCC cases, suggesting early alteration in the disease pathway. Hierarchical clustering revealed miR expression stratified tissues mostly according to phenotype.

Figure 1.

MiR expression in normal and malignant urothelium. Centroid linkage hierarchical clustering was performed after selecting miRs with >90% expression frequency and mean centering data using city block similarity in Cluster 3.0. First dendrogram row, subsequent progression (black box, progression; white box, no progression (also in brackets after sample ID (0, no; 1, yes)]. Boxes outline miRs with similar expression profiles.

Figure 1.

MiR expression in normal and malignant urothelium. Centroid linkage hierarchical clustering was performed after selecting miRs with >90% expression frequency and mean centering data using city block similarity in Cluster 3.0. First dendrogram row, subsequent progression (black box, progression; white box, no progression (also in brackets after sample ID (0, no; 1, yes)]. Boxes outline miRs with similar expression profiles.

Close modal
Table 2.

Differences in miR expression between normal and malignant urothelium

A. Absolute expression of a few miRs allows the distinction of malignant and normal urothelium
graphic
 
B. Sixteen MiRs have significant differential expression (P < 0.05 and FDR<0.05) in normal and malignant urothelium when compared* 
miR Chr Family Normal Urothelium ΔCt UCC ΔCt Fold change FDR q % P 
Present Mean St Dev Present Mean St Dev 
miR-649 22 miR-649 90% 18.59 4.11 90% 12.03 2.93 94.00 0.0 0.0019 
miR-135b miR-135 100% 7.45 1.53 98% 3.18 2.42 19.41 0.0 0.0037 
miR-133b miR-133 100% 0.59 1.93 98% 5.41 3.05 0.04 0.0 0.0041 
miR-520b 19 miR-93 20% 15.00 4.94 77% 11.64 3.94 10.28 0.0 0.0056 
miR-204 miR-204 100% 4.42 1.85 92% 9.20 3.40 0.04 0.0 0.0063 
miR-601 miR-601 90% 14.57 3.51 94% 9.92 2.79 24.98 0.0 0.0091 
miR-646 20 70% 14.58 3.39 94% 10.04 2.89 23.24 0.0 0.0160 
miR-639 19 miR-639 80% 13.29 2.63 90% 8.82 3.16 22.18 0.0 0.0225 
miR-380-5p 14 70% 4.41 4.69 75% 10.50 4.34 0.01 2.6 0.0285 
miR-21 17 miR-21 100% 2.45 1.56 100% −0.97 2.80 10.71 0.0 0.0296 
miR-644 20 miR-644 70% 17.57 3.63 92% 12.66 4.05 30.18 0.0 0.0305 
miR-93 miR-93 100% 2.38 2.05 100% −0.45 1.90 7.15 0.0 0.0308 
miR-15a 13 miR-15 100% 5.98 1.80 96% 2.89 2.39 8.50 0.0 0.0354 
miR-520d 19 miR-93 30% 21.04 1.02 75% 11.47 4.58 758.94 0.0 0.0394 
miR-549 15 miR-549 0%   56% 14.43 3.49 0.00 4.3 0.0427 
miR-449b miR-34 100% 12.17 2.24 90% 9.11 2.69 8.36 0.0 0.0439 
A. Absolute expression of a few miRs allows the distinction of malignant and normal urothelium
graphic
 
B. Sixteen MiRs have significant differential expression (P < 0.05 and FDR<0.05) in normal and malignant urothelium when compared* 
miR Chr Family Normal Urothelium ΔCt UCC ΔCt Fold change FDR q % P 
Present Mean St Dev Present Mean St Dev 
miR-649 22 miR-649 90% 18.59 4.11 90% 12.03 2.93 94.00 0.0 0.0019 
miR-135b miR-135 100% 7.45 1.53 98% 3.18 2.42 19.41 0.0 0.0037 
miR-133b miR-133 100% 0.59 1.93 98% 5.41 3.05 0.04 0.0 0.0041 
miR-520b 19 miR-93 20% 15.00 4.94 77% 11.64 3.94 10.28 0.0 0.0056 
miR-204 miR-204 100% 4.42 1.85 92% 9.20 3.40 0.04 0.0 0.0063 
miR-601 miR-601 90% 14.57 3.51 94% 9.92 2.79 24.98 0.0 0.0091 
miR-646 20 70% 14.58 3.39 94% 10.04 2.89 23.24 0.0 0.0160 
miR-639 19 miR-639 80% 13.29 2.63 90% 8.82 3.16 22.18 0.0 0.0225 
miR-380-5p 14 70% 4.41 4.69 75% 10.50 4.34 0.01 2.6 0.0285 
miR-21 17 miR-21 100% 2.45 1.56 100% −0.97 2.80 10.71 0.0 0.0296 
miR-644 20 miR-644 70% 17.57 3.63 92% 12.66 4.05 30.18 0.0 0.0305 
miR-93 miR-93 100% 2.38 2.05 100% −0.45 1.90 7.15 0.0 0.0308 
miR-15a 13 miR-15 100% 5.98 1.80 96% 2.89 2.39 8.50 0.0 0.0354 
miR-520d 19 miR-93 30% 21.04 1.02 75% 11.47 4.58 758.94 0.0 0.0394 
miR-549 15 miR-549 0%   56% 14.43 3.49 0.00 4.3 0.0427 
miR-449b miR-34 100% 12.17 2.24 90% 9.11 2.69 8.36 0.0 0.0439 

Abbreviations: FN, false negative; FP, false positive; Sens, sensitivity; Spec, specificity.

*The ΔCt mean and SD values are shown, together with fold change expression associated with carcinogenesis (calculated using the 2-ΔΔCt method).

miR-93.hd/291-3p/294/295/302/372/373/520.

Genes targeted by aberrantly expressed miR

With TargetScan, we identified potential target genes for the 16 differentially expressed miRs and selected the top one-third, as this fraction is most associated with altered protein expression (18). Of 1,095 predicted genes, 276 were targeted by more than one miR, including PLAG1 by 8 miRs, and E2F7 and DMTF1 by 5. Target ranking by context score revealed many carcinogenic genes with high binding affinity (Supplementary Fig. S2). We selected nine genes, from those with highest binding affinity and potentially carcinogenic roles, and measured their expression in our urothelial samples. For six of nine genes, there was differential expression between normal and malignant urothelial (Supplementary Fig. S2C, P < 0.02 t test). For five of these genes, this differential expression was inverse to that of their targeting miR, suggesting they play a regulatory role in expression (17). From UCC gene expression array data (24), we identified 143 putative target genes. For 76 (54%), an inverse relationship between miR and mRNA expression was present, and in 23 (11%), this difference reached significance (P < 0.05, t test, Supplementary Fig. S3).

Expression of miR with respect to tumor phenotype

Comparisons of the three tumor groups with disease-free normal urothelium revealed global and specific differential miR expression (Figs. 2 and 3). In general, there was a downregulation of miR in the low-grade samples (18 of 25 differentially expressed miRs) and an upregulation in high-grade and invasive tumors (28 of 30 and 20 of 20 miRs, respectively; χ2 P < 0.0001). Few aberrantly expressed miRs were common to the low-grade and high-grade/invasive cohorts [4 of 59 (6%)], whereas many were shared between the high-grade and invasive tumors [12 of 38 (32%), χ2 P = 0.003]. The magnitude of differential miR expression varied among tumor groups [average fold change, 4.32 (low grade), 28.24 (high-grade NMI), and 50.18 (invasive); ANOVA P = 0.006]. When the three groups were compared with each other, the low-grade tumors were distinct from the high-grade NMI (n = 9 miRs differentially expressed) and invasive cohorts (n = 45 miRs differentially expressed), whereas no significant differences in expression were present between the latter two (Fig. 3A).

Figure 2.

MiR expression in low- and high-grade bladder cancer. A, changes in expression (log scale) and statistical significance (SAM P values plotted as −log 10) of low-grade NMI tumors when compared with normal urothelium from non-UCC cases reveals both upregulation and downregulation of numerous miRs. B, for high-grade NMI and invasive tumors, there is a general upregulation in miR expression. Although changes in some miRs are shared between the two tumor phenotypes (e.g., miR-133b/204), others seem specific (miR-21 or miRs-100/99a). Expression is relative to the median of three RNA molecules and linear median centered.

Figure 2.

MiR expression in low- and high-grade bladder cancer. A, changes in expression (log scale) and statistical significance (SAM P values plotted as −log 10) of low-grade NMI tumors when compared with normal urothelium from non-UCC cases reveals both upregulation and downregulation of numerous miRs. B, for high-grade NMI and invasive tumors, there is a general upregulation in miR expression. Although changes in some miRs are shared between the two tumor phenotypes (e.g., miR-133b/204), others seem specific (miR-21 or miRs-100/99a). Expression is relative to the median of three RNA molecules and linear median centered.

Close modal
Figure 3.

MiR expression and tumor phenotype. A, differentially expressed miRs specific to each phenotype. B, signature panels of miR for each tissue. C, progression analysis in all 52 UCC revealed that dichotomous aberrant expression of miRs-21/100/99a was associated with progression to more advanced UCC (Kaplan-Meier method and tested using log-rank analysis).

Figure 3.

MiR expression and tumor phenotype. A, differentially expressed miRs specific to each phenotype. B, signature panels of miR for each tissue. C, progression analysis in all 52 UCC revealed that dichotomous aberrant expression of miRs-21/100/99a was associated with progression to more advanced UCC (Kaplan-Meier method and tested using log-rank analysis).

Close modal

We defined a specific signature for each tumor group by selecting those miRs that were differentially expressed when compared with both normal urothelium and one or more of the other two tumor groups (Fig. 3B). As many miRs were shared between the high-grade NMI and invasive tumor groups, we merged these to create a single defining signature. We also identified four miRs that were altered in UCC regardless of tumor phenotype.

Phenotype-specific miR targets

We were interested by the difference in miR expression between the low-grade and high-grade tumors. MiR upregulation was a feature of the high-grade pathway, and miRs-21/373 were among those with most prominence. Recent data reveal that miR-21 negatively regulates the p53 tumor suppressor pathway, leading to a loss of cellular control (25). MiR-373 is known to regulate prometastatic pathways (26), possibly through LATS2 suppression (of note, LATS2 was downregulated in our tumors; Supplementary Fig. S2; ref. 8).

Low-grade tumors were characterized by downregulation of many miRs. Of the 16 miRs with significant downregulation, 7 are predicted to target FGFR3 (miRs-99a/100/214/145/30a/125b/507) and 1 is the only highly conserved miR to target H-Ras (miR-218). We focused upon miRs-99a/100 as they target a highly conserved 8mer region within the FGFR3 UTR and their expression was inversely correlated with FGFR3 mRNA (r = −0.48 and −0.52, respectively; P < 0.0004; Fig. 4A). To investigate whether FGFR3 is a target of miRs-99a/100, we manipulated their expression in NHU cells. These cells represent normal urothelium, have a 200-fold higher miR99a/100 expression than RT4 cells and a low FGFR3 expression (data not shown). Following transfection with anti-miRs, a 70% to 80% reduction in miR expression (data not shown) resulted in an average 3-fold (miR-99a) and 6-fold (miR-100) upregulation of FGFR3 protein (Fig. 4). Direct targeting of the 3′ UTR region of the FGFR3 gene was confirmed by a Luciferase reporter assay (Fig. 4). Transfection of NHU cells with anti-miRs to 99a/100 and a plasmid containing Luciferase and a partial FGFR3 3′UTR resulted in a 2-fold (miR-99a) and 3.4-fold (miR-100) increase in Luciferase expression, when compared with cells transfected with a scrambled UTR sequence. To investigate the reverse relationship, we produced matching NHU-TERT cell lines with and without the S249C FGFR3 mutation. In each cell line regardless of FGF exposure, confluence or passage number the expression of both miRs was unchanged by the presence of the mutation (Fig. 4D), suggesting the epigenetic event is upstream of the gene mutation.

Figure 4.

Analyses of miR 99a/100 in UCC. A, expression of miRs 100/99a and FGFR3 are inversely correlated (DCt values shown). B, we found that a 70% to 80% miR-99a/100 knockdown induced by anti-miRs in NHU cells transcribed with Luciferase and the 3′UTR of FGFR3 was associated with increased Luciferase expression when compared with NHU cells transfected with a scrambled UTR sequence. C, upregulation of FGFR3 protein expression was found in NHU cells transfected with anti-miRs to 99a/100 when compared with negative controls (scrambled anti-miR sequence; Western blot images [FGFR3 is top band (splice variant seen below)] and photometric relative quantification standardized for β-actin and measured using ImageJ). D, the presence of the activating S249C FGFR3 mutation does not alter miR-99a/100 expression (relative miR and mRNA quantification). E, reversal of DNA hypermethylation leads to re-expression of 5 downregulated miRs in the low-grade cell line (RT4). Relatively few events were seen in the cell lines representative of high-grade and invasive bladder cancer (RT112 and EJ/T24). Four miRs (mir-127/124a/373/517c) reported to be silenced by DNA hypermethylation are also shown (28, 29).

Figure 4.

Analyses of miR 99a/100 in UCC. A, expression of miRs 100/99a and FGFR3 are inversely correlated (DCt values shown). B, we found that a 70% to 80% miR-99a/100 knockdown induced by anti-miRs in NHU cells transcribed with Luciferase and the 3′UTR of FGFR3 was associated with increased Luciferase expression when compared with NHU cells transfected with a scrambled UTR sequence. C, upregulation of FGFR3 protein expression was found in NHU cells transfected with anti-miRs to 99a/100 when compared with negative controls (scrambled anti-miR sequence; Western blot images [FGFR3 is top band (splice variant seen below)] and photometric relative quantification standardized for β-actin and measured using ImageJ). D, the presence of the activating S249C FGFR3 mutation does not alter miR-99a/100 expression (relative miR and mRNA quantification). E, reversal of DNA hypermethylation leads to re-expression of 5 downregulated miRs in the low-grade cell line (RT4). Relatively few events were seen in the cell lines representative of high-grade and invasive bladder cancer (RT112 and EJ/T24). Four miRs (mir-127/124a/373/517c) reported to be silenced by DNA hypermethylation are also shown (28, 29).

Close modal

MiR expression and tumor progression

If miRs-99a/100 and 21/373 characterize low-grade and high-grade UCC, one would expect their expression to be associated with tumor behavior. For each miR, there was a significant change in tumor progression rate according to their expression, in a manner similar to their associated genetic traits (log-rank P < 0.05; Fig. 3C). MiRs-99a/100 downregulation was associated with a better outcome (analogous to FGFR3 mutation), whereas miRs-21/373 upregulation was associated with a worse outcome (as with p53 mutation), when compared with contrasting tumors. Multivariate analysis revealed that miR-99a/100 expression was associated with progression when analyzed together with tumor stage, grade, and miR-21 expression (Cox MVA P = 0.03 and P = 0.006, respectively).

Epigenetic regulation of miR

The role of aberrant promoter hypermethylation for miR downregulation in low-grade UCC was analyzed using DNA methyltransferase inhibition. Following 5-azactytidine treatment, eight miRs (Fig. 4E) had >2-fold increase in expression, including six in the low-grade cell line (RT4). The lower rates of miR re-expression in the high-grade tumor lines (RT112 and EJ/T24) suggest that regulation by hypermethylation of these miRs is relatively low-grade specific. Both RT112 and EJ/T24 have higher rates of DNA hypermethylation than RT4 at known tumor suppressor gene promoters (27). For comparison, we analyzed four hypermethylated miRs reported in UCC (miRs-127/124a/373/517c; refs. 28, 29). For these, re-expression was usually shared between two or more cell lines rather than just isolated to RT4. MiRs-99a/100 are located at 21q21.1 and 11q24.1, respectively, together with two other miRs and one protein coding genes. Neither region contains a canonical CpG island. To investigate long-range epigenetic silencing of miRs-99a/100, we looked at the expression of neighboring genes following 5-azactytidine treatment (Supplementary Fig. S4). There was an increase in expression of mir-125b (located in duplicate genes, at each locus) but not of either protein coding gene nor the other miR in RT4.

Processing molecules and miR expression

We analyzed the expression of molecules known to be important for miR processing. The expression of each varied significantly between tissue type (ANOVA P < 0.05; Supplementary Fig. S5) but not between tumor phenotype. For Dicer, Drosha (RNASEN), and Exportin 5, there was an upregulation in the normal urothelium from UCC cases that was reversed once malignancy appeared. The largest changes were for Dicer and Drosha (7.4- and 6.1-fold upregulation, P = 0.00006 and 0.00001, respectively), whose expression was closely correlated (Pearson's correlation = 0.79, P = 7 × 10−19). For RAN and DGCR8 the opposite was seen, with an initial downregulation followed by upregulation in the tumors. Concentrations of these two mRNA were almost identical in each tissue (Pearson's correlation = 0.48, P = 3 × 10−6). When expression of these molecules was compared with that of the 322 miRs, several associations were seen. Dicer is targeted by several miRs and expression of miR-130b (r = −0.28, P = 0.02) and let-7g (r = −0.25, P = 0.05) were inversely correlated with Dicer expression.

Here, we have shown that altered miR expression occurs commonly in UCC, in a phenotype-specific manner and targets key pathways considered “hallmarks” of the disease. Our first finding was that protumorogenic miR alterations occur before the histologic onset of malignancy, supporting observations of promoter hypermethylation and genetic instability (27, 30). Almost half of our tumor-associated miRs were upregulated in the normal urothelium from patients with the disease, and this was associated with increased Dicer, Drosha, and Exportin 5 expression. Dicer overexpression in premalignant lesions of the prostate and lung, and reduced Dicer and Drosha expression in ovarian cancer are reported (3133), supporting our observation. Expression of the miR machinery was not associated with specific miRs, suggesting a global miR upregulation rather than the selection of tumor-specific species.

When analyzing all UCC cases together, we found the aberrant expression of 16 miRs, including some generic to carcinogenesis [e.g., miR-21 (ref. 25), mir-133b (ref. 34)]. Numerous interesting observations were apparent. For example, reduced miR-204 expression was frequent. MiR-204 is located on chromosome 9q21.11 within an intron of TRPM3, a gene downregulated in UCC (35, 36). MiR-204 potentially targets FRS2 an important member of the FGFR3 pathway (37). Observed changes in miR expression occurred in isolation or were clustered. Isolated differences arose from solitary located miRs such as miRs-649/135b or miR-601. The latter is located within intron 1 of DENND1A, a gene upregulated in UCC (38). Isolated differences also arose when a single member of a miR cluster was differentially expressed, including downregulation of miR-133b (adjacent to mir-206) or upregulation of miR-449b (located with miR-449a) and miR-93 (clustered with miRs-25/106b within the MCM7 gene). Clustered alterations included downregulation of miRs-133a/1 (located in duplicated clusters on Chr18q11.1 and Chr20q13.33), miRs-99a/125b, and miRs-143/145. Familial clustering was also present with three members of the miR-93 family (miR-93/520b/520d), miR-133a, and 133b, and three of the miR-34 family (miR-449b/34c/34a) having aberrant expression. Bioinformatic predictions identified around 1,000 genes as putative targets of these miRs. Of these genes, 20% were targeted by more than one miR. We confirmed consistent expression changes for some of these targets including new observations suggesting a carcinogenic role for LATS2, YOD1, and RAB22a in UCC. We identified miR as a potential mechanism for SOX4 upregulation in UCC (39).

Although the comparison of UCC with normal urothelium revealed interesting data, it did not detail events within the disease well. To obtain a detailed picture, it was necessary to analyze miR expression between tumor phenotypes. High-grade tumors were characterized by miR upregulation and for most, altered expression occurred before the onset of muscle invasion. Low-grade UCCs were characterized by miR downregulation and this affected many molecules targeting FGFR3 or H-Ras. Confirmation of FGFR3 miR targeting was obtained by the inverse correlation between miR and mRNA expression, and functional miR knockdown experiments. As aberrant expression of miRs-99a/100 was almost ubiquitous in low-grade UCC, occurs more frequently than FGFR3 mutation, and the presence of the commonest FGFR3 mutation (S249C) did not alter miR expression, we concluded the miR alteration leads to FGFR3 upregulation before the acquisition of gene mutation. This is an intriguing finding and suggests that epigenetic upregulation of FGFR3 facilitates the acquisition of mutation by either increasing cell turnover, reducing regulation to make cells more susceptible to carcinogens, or providing a selection advantage for mutant cells. To investigate potential causes of reduced miR expression in low-grade UCC, we examined DNA methylation. For six miRs, there was significant upregulation in RT4 following DNA methyltransferase inhibition, including miR-1 for which aberrant hypermethylation has recently been reported in hepatocellular carcinoma (40). Our findings with respect to miR-s99a/100 were less clear. There are no canonical CpG islands close to either gene, but two of the four genes clustered around miR-99a or miR-100 loci had increased expression following 5-azacytidine. This finding has been reported by others (40, 41) and suggests long-range epigenetic regulation (42). However, the expression of neighboring protein-coding genes did not alter with DNA Methyltransferase inhibition.

In summary, we have found new mechanisms related to biology and progression of UCC. These epigenetic events precede histologic changes of malignancy or disease-progression and occur in a phenotype-specific manner. These data have implications regarding out understanding of UCC tumor biology, and suggest miR as a novel diagnostic or prognostic biomarker, and therapeutic target.

No potential conflicts of interest were disclosed.

Grant support: GSK Clinician Scientist fellowship (J.W.F. Catto) and project grants (J.W.F. Catto) from Yorkshire Cancer Research, Sheffield Hospitals Charitable trust, and the European Union (Framework 7).

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

We thank Messrs. Anderson, Chapple, Hastie, Hall, Inman, Oakley, and Smith for allowing us to study their patients, Sister Louise Goodwin and Leila Ayandi for their help in tissue collection, Richard Bryant and Paul Heath for array advice, and the staff and patients of the Royal Hallamshire hospital.

1
Jemal
A
,
Siegel
R
,
Ward
E
, et al
. 
Cancer statistics, 2008
.
CA Cancer J Clin
2008
;
58
:
71
96
.
2
von der Maase
H
,
Sengelov
L
,
Roberts
JT
, et al
. 
Long-term survival results of a randomized trial comparing gemcitabine plus cisplatin, with methotrexate, vinblastine, doxorubicin, plus cisplatin in patients with bladder cancer
.
J Clin Oncol
2005
;
23
:
4602
8
.
3
Herr
HW
,
Dotan
Z
,
Donat
SM
,
Bajorin
DF
. 
Defining optimal therapy for muscle invasive bladder cancer
.
J Urol
2007
;
177
:
437
43
.
4
Catto
JW
,
Azzouzi
AR
,
Rehman
I
, et al
. 
Promoter hypermethylation is associated with tumor location, stage, and subsequent progression in transitional cell carcinoma
.
J Clin Oncol
2005
;
23
:
2903
10
.
5
Knowles
MA
. 
Molecular subtypes of bladder cancer: Jekyll and Hyde or chalk and cheese?
Carcinogenesis
2006
;
27
:
361
73
.
6
Stefani
G
,
Slack
FJ
. 
Small non-coding RNAs in animal development
.
Nat Rev Mol Cell Biol
2008
;
9
:
219
30
.
7
He
L
,
Thomson
JM
,
Hemann
MT
, et al
. 
A microRNA polycistron as a potential human oncogene
.
Nature
2005
;
435
:
828
33
.
8
Voorhoeve
PM
,
le Sage
C
,
Schrier
M
, et al
. 
A genetic screen implicates miRNA-372 and miRNA-373 as oncogenes in testicular germ cell tumors
.
Cell
2006
;
124
:
1169
81
.
9
Volinia
S
,
Calin
GA
,
Liu
CG
, et al
. 
A microRNA expression signature of human solid tumors defines cancer gene targets
.
Proc Natl Acad Sci U S A
2006
;
103
:
2257
61
.
10
Schetter
AJ
,
Leung
SY
,
Sohn
JJ
, et al
. 
MicroRNA expression profiles associated with prognosis and therapeutic outcome in colon adenocarcinoma
.
JAMA
2008
;
299
:
425
36
.
11
Yu
SL
,
Chen
HY
,
Chang
GC
, et al
. 
MicroRNA signature predicts survival and relapse in lung cancer
.
Cancer Cell
2008
;
13
:
48
57
.
12
Gottardo
F
,
Liu
CG
,
Ferracin
M
, et al
. 
Micro-RNA profiling in kidney and bladder cancers
.
Urol Oncol
2007
;
25
:
387
92
.
13
Neely
LA
,
Rieger-Christ
KM
,
Neto
BS
, et al
. 
A microRNA expression ratio defining the invasive phenotype in bladder tumors
.
Urol Oncol
2008
,
Epub
.
14
Catto
JW
,
Rosario
DJ
. 
The road to cystectomy: who, when and why?
Eur Urol, EAU Update series
2005
;
3
:
171
9
.
15
Southgate
J
,
Hutton
KA
,
Thomas
DF
,
Trejdosiewicz
LK
. 
Normal human urothelial cells in vitro: proliferation and induction of stratification
.
Lab Invest
1994
;
71
:
583
94
.
16
Chen
C
,
Ridzon
DA
,
Broomer
AJ
, et al
. 
Real-time quantification of microRNAs by stem-loop RT-PCR
.
Nucleic Acids Res
2005
;
33
:
e179
.
17
Selbach
M
,
Schwanhausser
B
,
Thierfelder
N
,
Fang
Z
,
Khanin
R
,
Rajewsky
N
. 
Widespread changes in protein synthesis induced by microRNAs
.
Nature
2008
;
455
:
58
63
.
18
Baek
D
,
Villen
J
,
Shin
C
,
Camargo
FD
,
Gygi
SP
,
Bartel
DP
. 
The impact of microRNAs on protein output
.
Nature
2008
;
455
:
64
71
.
19
Chapman
EJ
,
Hurst
CD
,
Pitt
E
,
Chambers
P
,
Aveyard
JS
,
Knowles
MA
. 
Expression of hTERT immortalises normal human urothelial cells without inactivation of the p16/Rb pathway
.
Oncogene
2006
;
25
:
5037
45
.
20
Tomlinson
DC
,
L'Hote
CG
,
Kennedy
W
,
Pitt
E
,
Knowles
MA
. 
Alternative splicing of fibroblast growth factor receptor 3 produces a secreted isoform that inhibits fibroblast growth factor-induced proliferation and is repressed in urothelial carcinoma cell lines
.
Cancer Res
2005
;
65
:
10441
9
.
21
Schmittgen
TD
,
Lee
EJ
,
Jiang
J
, et al
. 
Real-time PCR quantification of precursor and mature microRNA
.
Methods
2008
;
44
:
31
8
.
22
Tusher
VG
,
Tibshirani
R
,
Chu
G
. 
Significance analysis of microarrays applied to the ionizing radiation response
.
Proc Natl Acad Sci U S A
2001
;
98
:
5116
21
.
23
Tibshirani
R
,
Hastie
T
,
Narasimhan
B
,
Chu
G
. 
Diagnosis of multiple cancer types by shrunken centroids of gene expression
.
Proc Natl Acad Sci U S A
2002
;
99
:
6567
72
.
24
Wild
PJ
,
Herr
A
,
Wissmann
C
, et al
. 
Gene expression profiling of progressive papillary noninvasive carcinomas of the urinary bladder
.
Clin Cancer Res
2005
;
11
:
4415
29
.
25
Papagiannakopoulos
T
,
Shapiro
A
,
Kosik
KS
. 
MicroRNA-21 targets a network of key tumor-suppressive pathways in glioblastoma cells
.
Cancer Res
2008
;
68
:
8164
72
.
26
Huang
Q
,
Gumireddy
K
,
Schrier
M
, et al
. 
The microRNAs miR-373 and miR-520c promote tumour invasion and metastasis
.
Nat Cell Biol
2008
;
10
:
202
10
.
27
Dhawan
D
,
Hamdy
FC
,
Rehman
I
, et al
. 
Evidence for the early onset of aberrant promoter methylation in urothelial carcinoma
.
J Pathol
2006
;
209
:
336
43
.
28
Lujambio
A
,
Ropero
S
,
Ballestar
E
, et al
. 
Genetic unmasking of an epigenetically silenced microRNA in human cancer cells
.
Cancer Res
2007
;
67
:
1424
9
.
29
Saito
Y
,
Liang
G
,
Egger
G
, et al
. 
Specific activation of microRNA-127 with downregulation of the proto-oncogene BCL6 by chromatin-modifying drugs in human cancer cells
.
Cancer Cell
2006
;
9
:
435
43
.
30
Catto
JW
,
Hartmann
A
,
Stoehr
R
, et al
. 
Multifocal urothelial cancers with the mutator phenotype are of monoclonal origin and require panurothelial treatment for tumor clearance
.
J Urol
2006
;
175
:
2323
30
.
31
Chiosea
S
,
Jelezcova
E
,
Chandran
U
, et al
. 
Up-regulation of dicer, a component of the MicroRNA machinery, in prostate adenocarcinoma
.
Am J Pathol
2006
;
169
:
1812
20
.
32
Chiosea
S
,
Jelezcova
E
,
Chandran
U
, et al
. 
Overexpression of Dicer in precursor lesions of lung adenocarcinoma
.
Cancer Res
2007
;
67
:
2345
50
.
33
Merritt
WM
,
Lin
YG
,
Han
LY
, et al
. 
Dicer, Drosha, and outcomes in patients with ovarian cancer
.
N Engl J Med
2008
;
359
:
2641
50
.
34
Bandres
E
,
Cubedo
E
,
Agirre
X
, et al
. 
Identification by Real-time PCR of 13 mature microRNAs differentially expressed in colorectal cancer and non-tumoral tissues
.
Mol Cancer
2006
;
5
:
29
.
35
Blaveri
E
,
Brewer
JL
,
Roydasgupta
R
, et al
. 
Bladder cancer stage and outcome by array-based comparative genomic hybridization
.
Clin Cancer Res
2005
;
11
:
7012
22
.
36
Dyrskjot
L
,
Kruhoffer
M
,
Thykjaer
T
, et al
. 
Gene expression in the urinary bladder: a common carcinoma in situ gene expression signature exists disregarding histopathological classification
.
Cancer Res
2004
;
64
:
4040
8
.
37
L'Hote
CG
,
Knowles
MA
. 
Cell responses to FGFR3 signalling: growth, differentiation and apoptosis
.
Exp Cell Res
2005
;
304
:
417
31
.
38
Sanchez-Carbayo
M
,
Socci
ND
,
Lozano
J
,
Saint
F
,
Cordon-Cardo
C
. 
Defining molecular profiles of poor outcome in patients with invasive bladder cancer using oligonucleotide microarrays
.
J Clin Oncol
2006
;
24
:
778
89
.
39
Aaboe
M
,
Birkenkamp-Demtroder
K
,
Wiuf
C
, et al
. 
SOX4 expression in bladder carcinoma: clinical aspects and in vitro functional characterization
.
Cancer Res
2006
;
66
:
3434
42
.
40
Datta
J
,
Kutay
H
,
Nasser
MW
, et al
. 
Methylation mediated silencing of MicroRNA-1 gene and its role in hepatocellular carcinogenesis
.
Cancer Res
2008
;
68
:
5049
58
.
41
Lujambio
A
,
Calin
GA
,
Villanueva
A
, et al
. 
A microRNA DNA methylation signature for human cancer metastasis
.
Proc Natl Acad Sci U S A
2008
;
105
:
13556
61
.
42
Frigola
J
,
Song
J
,
Stirzaker
C
,
Hinshelwood
RA
,
Peinado
MA
,
Clark
SJ
. 
Epigenetic remodeling in colorectal cancer results in coordinate gene suppression across an entire chromosome band
.
Nat Genet
2006
;
38
:
540
9
.