Skip to main content

Urinary expression of let-7c cluster as non-invasive tool to assess the risk of disease progression in patients with high grade non-muscle invasive bladder Cancer: a pilot study

Abstract

Background

High grade non-muscle-invasive bladder cancer (HG-NMIBC) is a heterogeneous disease with variable risk of progression. Urinary microRNAs are promising biomarkers for BC detection and surveillance. Let-7c-5p miRNA, clustered with miR-99a-5p and -125b-5p, is deregulated in cancer, including BC. The aim of this study is to evaluate urinary let-7c cluster expression in Ta/T1 HG-NMIBC patients and its impact on progression-free survival (PFS).

Methods

Quantitative Real-Time-Polymerase-Chain-Reaction (qRT-PCR) was used to analyze the let-7c cluster expression in 57 urine and 49 neoplastic paired tissue samples prospectively collected from transurethral resection (TUR) HG-NMIBC patients. Twenty urine and 10 bladder tissue samples were collected and analyzed as normal controls. QRT-PCR was also used to detect intra−/extra-cellular let-7c cluster in BC cells. Receiver Operating Characteristic (ROC) curves were used to identify urinary miRNAs cut-off values predicting T-stage and PFS. Uni/multivariable Cox regression was performed to identify predictors of PFS. A nomogram predicting progression risk and a decision curve analysis (DCA) were performed.

Results

Urinary let-7c was significantly up-regulated in patients compared with controls, while the whole cluster was down-regulated in tumor tissues. Supporting these findings, in vitro comparison of extra−/intra-cellular ratios of cluster levels between BC cells, showed a higher ratio for let-7c in HG-NMIBC versus low-grade cells. Urinary let-7c cluster expression was increased in higher T-stage and was an independent predictor of progression. Lower EORTC-score and downregulation of urinary cluster were predictors of higher PFS on univariable Cox regression, while on multivariable analysis only cluster expression was an independent progression predictor. On DCA, a benefit was evident for patients with a PFS probability > 20%.

Conclusions

Urinary let-7c cluster evaluation may improve prognosis, identifying patients at risk of progression and addressing early radical treatment.

Background

Among all bladder tumors, 70–75% are non-muscle-invasive (NMIBC: 70% Ta, 20% T1, 10% Tis) at first evaluation. NMIBC is a heterogeneous group composed of stages Ta (non-invasive papillary carcinoma), T1 (sub-epithelial connective tissue invasive carcinoma) and CIS (Tis Carcinoma in situ), and various grades including well (G1) moderately (G2) and poor (G3) differentiated tumors. Many large cohorts studies analyzed the long-term survival of NMIBC patients in general after transurethral resection (TUR) and intravesical chemotherapy [1]. Survival range varies between studies but are mostly higher than in muscle-invasive bladder cancer (MIBC), high grade (HG) Ta/T1 disease represents about 5 to 23% of NMIBC at first diagnosis, with a steady increase in incidence. It is considered as a highly malignant tumor with a variable and unpredictable biology as well as a critical aspect of management. Most of NMIBC cases are treated by TUR of the tumor, but about 50–70% undergoes recurrence within 2 years and 10–30% progresses to muscle-invasive disease [1]. The most accepted predictors of progression in these patients are tumor size, multifocality, presence of CIS and vascular invasion [2]. According to the European Organization for Research and Treatment of Cancer (EORTC) classification, any T1, HG/G3, CIS, recurrent multifocal and larger than 3 cm Ta G1-G2/Low grade (LG), are considered bladder cancer (BC) at high risk of progression to muscle-invasive disease [1, 3]. However, these features may not accurately reflect clinical outcome and only partially address a reliable risk-adjusted treatment planning. In the era of precision medicine, many efforts have been made to increase knowledge about disease biology, to find reliable biomarkers associated with progression in BC [3].

MicroRNA (miRNAs), short noncoding RNAs regulators of gene expression [4], are also present in several body fluids such as serum, plasma and urine, with a high degree of stability indicating their extensive potential as biomarkers [5, 6]. Recent literature has highlighted the dysregulation of several miRNAs, among them let-7c cluster, in BC tissue as well as in urine and blood, suggesting their potential role as diagnostic and prognostic biomarkers in BC patients [7,8,9,10,11,12].

Let-7c belongs to the let-7 miRNA family that plays an important role in several cellular processes, and it is widely viewed as a tumor suppressor miRNA [13, 14]. Let-7c coding sequence, organized in a cluster with miR-99a and miR-125b, is localized on chromosome 21 and resides within intron 6 of the poorly characterized and recently named LINC00478 gene [15]. Consistent with its tumor suppressor activity, let-7c is downregulated in many cancer, including BC, and has been related to tumor progression [16,17,18].

Here, we assessed the impact on Progression Free Survival (PFS) probabilities of urinary let-7c cluster expression to evaluate its suitability as potential non-invasive biomarker associated with NMIBC patients at higher risk of progression to muscle-invasive disease.

Methods

Patients and specimens

From January 2015 to March 2018, 57 urine and 49 neoplastic paired tissue samples were prospectively collected from TUR of HG-NMIBC patients at the “Regina Elena” National Cancer Institute. For 8 patients the tissue sampling was not available for experimental purposes. Twenty urine and ten bladder tissue samples were collected as healthy controls (HC), matched for age, sex, smoking habit and comorbidities. All patients received Bacillus Calmette-Guerin (BCG) induction and maintenance for a minimum period of 1 year according to Kamat’s protocol [19].

Urine cytology and cystoscopy were performed at 3 months’ intervals for the first two years and at 6 months’ intervals thereafter. A whole-body CT scan was performed yearly. TUR was performed in case of BC recurrence. Disease progression was defined as evidence of muscle-invasion at any subsequent TUR or evidence of metastatic spread at conventional imaging.

Urine and tissue samples processing and RNA extraction

Urine was analyzed from each HC recruited at the same Institute from individuals seeking a routine health check-up, with no evidence of disease and with age-, sex- and ethnicity-matched to the patients. Urine samples obtained immediately before TUR procedure were analyzed from BC patients. Patients with positive urine culture and/or leukocyturia were excluded from the study. Normal tissue samples were obtained from bio-bank specimens of patients who underwent radical cystectomy (RC) [20].

Urine samples were collected in sterile recipient, stored at 4 °C, then centrifuged at 4.500 rpm for 30 min at 4 °C, followed by centrifugation at 8.900 rpm for 5 min to separate the supernatant from the urinary cell sediment. Urine supernatant was aliquoted in 2 mL-cryovials and stored at − 80 °C. After surgery, BC tissue samples were collected and stored in liquid nitrogen.

Total RNA was extracted from urine supernatants (starting from 2 mL) using the Urine cell free circulating RNA purification midi kit (#57000, Norgen Biotek, Canada). During extraction we added spike-in mix of non-human synthetic miRNAs (RNA spike-in mix: UniSp2, UniSp4 and UniSp5; #203203, Exiqon) into the lysis buffer before incorporating the urine sample. Frozen tissue specimens were homogenized and total RNA was obtained using Animal Tissues RNA Purification Kit (#25700, Norgen Biotek, Canada). RNA quantification from tissues was done by NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Wilmington, DE U.S.A.).

Cell culture and RNA extraction

RT-112 and T24 BC cell lines were kindly provided by Prof. Claudio Sette (“Università Cattolica del Sacro Cuore”, Rome, Italy). The RT-4 cell line was purchased from ATCC (Manassas, VA, USA). Human RT-4, RT-112 and T24 BC cell lines corresponding to LG-NMIBC, HG-NMIBC and HG-MIBC cell phenotype respectively, were grown in L-Glutamine-RPMI-1640 with 10% fetal bovine serum (Gibco®-Thermo Fisher), penicillin/streptomycin, gentamicin and MEM non-essential aminoacids (Sigma-Aldrich) at 37 °C in 5% CO2. Total RNA extraction from the culture medium (extracellular fraction), starting from 400 μL,was performed using the miRNeasy Serum/Plasma kit Adavanced (#217204, Qiagen). As technical control, we added spike-in mix of non-human synthetic miRNAs (RNA spike-in mix: UniSp2, UniSp4 and UniSp5; #203203, Exiqon) into the lysis buffer, before adding it to the supernatant sample. For culture medium collection, cells were seeded and after 24 h the medium was replaced with Serum-Free Medium (SFM). Culture medium was centrifuged at 2000 x g for 10 min, filtered using 0.22 μm polyether-sulfonate low-protein binding filters and aliquoted. Total RNA was isolated from cells (intracellular fraction) after 24 h from SFM (along with the cell supernatant collection). Cell pellet from BC cells was lysed in TRIsure™ (Bioline, Meridian Bioscience) for total RNA extraction. RNA quantification from cells was done by NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Wilmington, DE U.S.A.).

Complementary DNA (cDNA) synthesis and qRT-PCR

RNA reverse transcription (starting from 4 μL of RNA from urine and culture medium or 20 ng of tissue and cell RNA) from urine, tissues, cells and culture medium was performed using the Universal cDNA synthesis kit II (#203301, Exiqon, Denmark), including in the reaction spike-in Unisp6, to control cDNA synthesis reaction quality.

Mature cell-free (urine and culture medium) and tissue miRNAs quantification was performed by miRNA-specific LNA™-based system using SYBR® Green (miRCURY LNA™ Universal RT micro-RNA PCR; Exiqon # 203301, Vedbaek- Denmark). For qRT-PCR, cDNA was diluted 1:40 in nuclease-free water. ROX passive reference dye was added in diluted cDNA samples to obtain a robust read over the entire array of wells (ROX, Thermo Scientific, Waltham, Massachusetts, USA). The mix of qRT-PCR was composed of 5 μL ExiLENT SYBR Green Mastermix, 1 μL miRNA probe specific for hsa-let-7c-5p (LNA UniRT primer mix #204767, Exiqon), hsa-miR-99a-5p (LNA UniRT primer mix #205713, Exiqon) and hsa-miR-125b-5p (LNA UniRT primer mix #204521, Exiqon) miRNAs and 4 μL cDNA. QRT-PCR was performed using the ABI 7900 Real Time PCR System and SDS 2.2.2 software (Applied Biosystems, Foster City, CA). All the samples were tested in triplicate. We normalized urinary and culture medium miRNA levels with the expression of UniSP2 synthetic spike-in (UniSP2 LNA control primer set UniRT, #203950, Exiqon), according to previous report [21]. MiRNAs expression levels from tissue and cells were normalized using U6 snRNA (U6 snRNA LNA primer set UniRT, #203907, Exiqon). Relative quantification of miRNAs expression was determined using the comparative threshold cycle (ΔCt) method. MiRNA expression levels from culture medium was additionally normalized with the cell count.

Statistical analysis

Descriptive statistics were used to summarize study data. The association between clinical variables and miRNA levels was tested by the Pearson Chi-Square test or Fisher’s Exact test, when appropriate. Differences of urinary and tissue miRNA levels between neoplastic and HC groups were compared using the non-parametric Mann-Whitney U test. MiRNAs expression levels were plotted into box-plots. Receiver operative characteristics (ROC) curves were used to identify the most informative miRNA expression thresholds predicting PFS. The prognostic role of miRNA cluster was tested with the Kaplan-Meier method. Survival rates were computed at 6, 12, 18 and 24 months after surgery, and the log-rank test was applied to assess statistical significance between the two groups. Univariable and multivariable Cox regression analyses were performed to identify predictors of PFS probabilities. A nomogram predicting disease progression risk was built. The discrimination accuracy was measured by concordance index (C-index). Calibration plot was generated with 200 bootstrap resampling. A decision curve analysis (DCA) was performed to assess the net benefit of the predictive model. Statistical analysis was performed with the Statistical Package for Social Science (SPSS Inc., v 21.0, Chicago, IL, USA), MedCalc (v. 14.10.2) and with STATA (Stata Statistical Software: Release 15. College Station, TX: StataCorp LLC).

Results

Expression levels of urinary let-7c/miR-99a/miR-125b cluster in HG-NMIBC patients

To investigate the diagnostic role of cell-free urinary let-7c/miR-99a/miR-125b cluster in BC we assessed, by qRT-PCR, their expression levels in a selected cohort of HG Ta/T1 NMIBC which includes newly diagnosed and first recurrence patient samples (n = 57), compared to urine let-7c cluster expression levels of non-tumor HC (n = 20) matched for age and sex. Demographic, clinical and pathologic features of the cohort and general characteristics of HC are reported in Table 1. As shown in Fig. 1 panel A, HG-NMIBC patients displayed a significant up-regulation of urinary let-7c (+ 2.57 fold; p = 0.047), comparable levels of miR-99a and a trend toward significantly higher levels of miR-125b when compared to HC. We also evaluated the cluster expression levels in paired tissue (n = 49) of HG-NMIBC patients, compared to age- and sex- matched non-tumor bladder tissues of HC (n = 10). The results showed that the expression levels of the cluster were significantly decreased in tumor tissues compared to HC (let-7c, − 27.75 fold, p < 0.001; miR-99a, − 50.40 fold, p < 0.001; miR-125b, − 27.91 fold, p < 0.001; Fig. 1 panel B).

Table 1 Demographic, clinical and pathologic features of patients and general characteristics of healthy controls
Fig. 1
figure 1

Expression levels of urinary and tissue let-7c miRNA cluster in a cohort of HG-NMIBC patients. Box-Plot diagrams of relative miRNA expression levels, performed by qRT-PCR in urine supernatant (a) and in tissue (b). Boxes define the 25th and 75th percentiles; the horizontal lines into within the boxes indicates the median, and bars define the minimum and maximum values. Relative miRNAs expression was normalized to Unisp2 spike-in for urine and U6 snRNA for tissue and calculated using the comparative ∆Ct method. P-value was determined using the Mann-Whitney rank sum test

Let-7c cluster expression levels and clinic-pathological features of HG-NMIBC patients

We evaluated the relationship between urinary and tissue expression levels of the let-7c cluster and the clinic-pathological parameters of the HG-NMIBC patients (Table 2). The results revealed a significant association between urinary let-7c cluster and tumor stage. In particular, higher expression of urinary let-7c and, to a lesser extent, of miR-99a/−125b were found in patients with T1 stage, compared to Ta stage (let-7c, + 2.46 fold, p = 0.0006; miR-99a, + 1.87 fold, p = 0.022; miR-125b, + 2.73 fold, p = 0.012; Fig. 2 panel A). In contrast, comparable expression of let-7c miRNA cluster in tumor tissues was found in Ta and T1 stage disease (p > 0.05; Figure S1 panel A). Next, we evaluated, by ROC curves analysis, the diagnostic performance and accuracy of urinary let-7c cluster to discriminate the tumor stage in HG-NMIBC patients. We found that urinary let-7c expression levels with an area under the curve (AUC) value of 0.80 (95% C-index = 0.670–0.892) were robust in discriminating T1 from Ta stage. In contrast, our analysis indicates a lower diagnostic performance for the other members of the cluster (miR-99a: AUC = 0.635, miR-125b: AUC = 0.691, Fig. 2 panel B). As expected, no substantial diagnostic performance was observed, by ROC curves, in paired tissues of the entire let-7c cluster (Figure S1 panel B; let-7c and miR-99a: AUC 0.60; miR-125b: 0.68 respectively). These data suggest that urinary let-7c cluster could be a reliable biomarker to distinguish T1 from the Ta stage in HG-NMIBC. The levels of let-7c cluster were not significantly associated with age, sex, tumor history, tumor size, smoking status, diabetes and hypertension (Table 2).

Table 2 Relationship between let-7c cluster and clinic-pathological features of NMIBC patients
Fig. 2
figure 2

Association analysis between urinary let-7c cluster and tumor stage. a Box-plot diagrams of urinary let-7c cluster stratified into pathological Ta and T1 stages. Relative expression levels were normalized to Unisp2 spike-in and calculated using the comparative ΔCt method. P-value were determined using the Mann-Whitney rank sum test. Boxes define the 25th and 75th percentiles; the horizontal lines into in the boxes indicates the median, and bars define the minimum and maximum values. b ROC curves plotted for diagnostic potential and discriminatory accuracy of urinary let-7c cluster to distinguish T1 from Ta stage. The corresponding Area Under the Curve (AUC) values are reported

Evaluation of let-7c cluster expression levels in the extracellular fraction of BC cell lines

Cultured BC cells were used to investigate let-7c cluster expression levels in intracellular (cells) and extracellular (culture medium) fractions of BC cell line representative of the HG-NMIBC phenotype. We detected let-7c and miR-125b, but not miR-99a, both in intracellular and extracellular fractions although, as expected, with a relatively lower baseline amount in the extracellular fraction (Fig. 3 panel A). Next, we evaluated the extracellular to intracellular ratio of let-7c and miR-125b levels, comparing BC cells representative of different grade and invasiveness capability. Interestingly, we observed a higher ratio for let-7c in HG-NMIBC (RT-112) cells compared to the LG-NMIBC (RT-4) cell line, while no differences were detected compared to HG-MIBC (T24; + 2.61 fold, p = 0.015; Fig. 3 panel B).

Fig. 3
figure 3

Analysis of let-7c miRNA cluster levels in BC cell lines. a Levels of let-7c cluster analyzed by qRT-PCR in intracellular (cells) and extracellular (culture medium) fractions of high grade (HG) NMIBC cells (RT-112). Relative expression was normalized with U6 snRNA and Unisp2 spike-in for intracellular and extracellular fraction respectively, and calculated using the comparative ∆Ct method. The extracellular values were normalized also to the cell count. N.D. indicates not determined miR-99a levels. b qRT-PCR of let-7c cluster of extracellular versus intracellular fraction (EC/IC ratio) in low grade (LG) NMIBC (RT-4), HG-NMIBC (RT-112) and HG-MIBC (T24) cells. Data are the average of three independent experiments; error bars indicate SEM. P-value was determined using the Mann-Whitney rank sum test

Urinary let-7c cluster as predictor of PFS probability

On ROC analysis, the most informative urinary levels cut-off values to discriminate PFS probabilities were 0.010 (AUC = 0.76), 0.035 (AUC = 0.77) and 0.020 (AUC 0.67) for let-7c, miR-125b and miR-99a, respectively (Figure S2). At a median follow up of 17 months (Interquartile range- IQR 12.5–22), 12 patients (21%) experienced disease progression. Upon univariable Cox regression analysis, lower EORTC scores and downregulation of the entire cluster were significant predictors of PFS probability (p=0.049 and p< 0.001 respectively; Table 3). Kaplan Meier analysis indicated that downregulation of the entire cluster was associated with significantly higher PFS even after stratifying for EORTC Risk Class 2 or 3 (both log rank p < 0.001, Fig. 4). Multivariable Cox regression analysis revealed that, only the cluster expression was an independent predictor of disease progression (HR 0.07; [0.02–0.26] p < 0.001; Table 3). The developed nomogram had a high predictive accuracy (C-index = 0.87) and was perfectly calibrated (Figure S3 panel A and B). On DCA, the net benefit of using the model was evident for all patients with a PFS probability > 20% (Fig. 5).

Table 3 Univariable and Multivariable Cox Regression Analysis for PFS
Fig. 4
figure 4

PFS probability in NMIBC according to urinary let-7c cluster expression. Log-Rank Curves showing PFS probability of patients with NMIBC according to the expression of let-7c cluster in urine and after stratifying for EORTC Risk Class

Fig. 5
figure 5

Decision Curve Analysis. DCA (Decision Curve Analysis) providing the benefit of using the nomogram including let-7c cluster expression in urine and the EORTC Score. The net benefit was evident for all the patients with a Progression Free Survival PFS probability Pr (pfs) > 20%

Discussion

Pathological and clinical features of BC, such as the tumor stage/grade, metastasis, recurrence and chemo-sensitivity are well known predictors of cancer specific survival [3]. Since HG-NMIBC patients display a non-negligible risk of recurrence and progression, BCG bladder instillation (induction and maintenance) is recommended due to its proven efficacy on reducing both recurrence and progression rates [3]. A recent EORTC study confirmed the role and the optimal schedule of BCG treatment following TUR [1]. Up to now, given the high risk of disease progression in BCG-failure patients, RC has been the gold standard treatment for these patients [3]. Early detection of patients with NMIBC at higher risk of progression by a non-invasive molecular marker, avoiding unsuccessful BCG treatment, and early identification of patients who may really benefit from an early RC, are of critical importance especially with regards to precision medicine.

Analysis of miRNAs seems to be the most promising strategy for the development of non-invasive early diagnostic/prognostic biomarkers because they are highly stable within urine, are more resistant against nuclease degradation and require little handling care [5,6,7].

Let-7c miRNA is widely accepted as a tumor suppressor miRNA involved in the regulation of oncogenic pathways and downregulated in many types of tumor tissue, including BC [13,14,15,16,17,18]. In line with let-7c tumor-suppressor role, our analysis of the cluster expression levels in HG-NMIBC patient’s tissue showed that let-7c, as well as miR-99a and miR-125b, were highly downregulated compared to normal tissue samples. In urine, let-7c, as well as the other let-7c cluster members, miR-99a and miR-125b has been found down-regulated in BC versus HC [11, 22, 23]. In the current study, we investigated the potential value of urinary let-7c cluster as non-invasive diagnostic/prognostic biomarker of HG-NMIBC patients. Our data showed a significant upregulation of urinary let-7c cluster in BC patients’ samples compared to HC. This observation is in contrast with previous findings showing decreased levels of urinary let-7c cluster in BC patients’ samples [11, 22, 23].

Potential reasons for divergences with previous studies may resides in the several challenges related to studying circulating miRNAs, especially the urinary ones. Among these a high degree of inter-individual variability of extracellular miRNAs levels [6], and technical sources of variation between studies (e.g., the RNA extraction method, the normalization, the technological platforms used etc.). For these reasons, to avoid technical biases, we maintained constant sample processing through the whole study. Moreover, it has to be considered that the debate about the relative contribution of organs and/or tissue to miRNA in biological fluids is still open since detailed knowledge of molecular mechanisms governing miRNA release from normal and tumor tissues is still lacking [6, 24, 25]. However, several authors in previous reports have shown opposite miRNA expression profile in bio-fluids versus matched tumour tissues [26, 27]. Our results showing that urinary let-7c expression is significantly upregulated and, coherently with its role of tumour suppressor, down-regulated in matched tissue samples compared to normal controls, fit with these cases. Moreover, we investigated the let-7c cluster levels in the extracellular fraction of BC cells characterized by the phenotype of HG-NMIBC.

It has been hypothesized that extracellular miRNAs may be released through both active and passive mechanisms [24]. In vitro studies also proved that miRNAs are exported into the cell supernatant [27]. Considering this, we wondered whether BC cells, with different grades and invasive power, could secrete let-7c cluster miRNAs in the cell medium, and examined its levels originating directly from the tumor. Interestingly, we found that let-7c cluster is detected in the culture medium (extracellular fraction) of BC cells with HG-NMIBC phenotype. Of note, by comparing the extra- to intra-cellular ratio between let-7c cluster levels in different grades and invasive power BC cells, we found let-7c ratio significantly higher in HG-NMIBC compared the other cell phenotypes. Accordingly, a previous study reported an enrichment of let-7 miRNA family in the culture medium of gastric cells in relation with different tumorigenic and invasiveness propensities [27]. These data suggest that the increase of extracellular let-7c in HG-NMIBC cells could be a marker of disease at higher risk of progression to muscle-invasive BC. Moreover, the increased release of the tumor-suppressor let-7c in the culture medium found in the HG- compared to the LG-NMIBC cells could suggest a reduction of anti-cancer effect of let-7c in the cells, to maintain their oncogenic potential, as already observed in other tumor systems [27].

Several reports associate miRNA expression to cancer stage, also in BC [8, 11, 28]. In the present study we found a relevant association of increased urinary let-7c cluster levels in T1 stage disease. These data, even if provided by histological examination from TUR, suggests that the increased urinary let-7c cluster levels could early discriminate patients at higher risk of tumor progression as T1 HG-NMIBC. In matched tumour tissues, the T1 versus Ta tumour stages comparison did not lead to substantial statistically differences, thus, only the urinary miRNAs expression shows the capability of discriminating the higher and more advanced T1 stage with a greater risk of tumor progression in HG-NMIBC patients from the Ta stage. These data, in liquid biopsy vs tissue specimens, suggest the importance of let-7c cluster in view of the long-term goal to potentially reduce or avoid all forms of invasive procedures. Regarding the prognostic stratification of NMIBC patients, EORTC tables, based on the six most significant clinical and pathological features, are well-known tools to calculate the risk of disease progression in NMIBC patients [29]. Risk stratification by EORTC progression score was confirmed in many reports. May et al. reported, in a series of 521 RC performed for MIBC, that baseline high risk NMIBC had significantly worse cancer survival after RC when compared with baseline low-to-intermediate EORTC progression risk score patients [30]. In our study, upregulation of urinary let-7c cluster has an impact on PFS probabilities thus suggesting that its expression may improve prognosis estimation, identifying patients at risk of progression. EORTC score and let-7c cluster expression under the threshold value, determined in ROC curve analysis, were significant predictors of PFS on univariable Cox regression analysis. On multivariable analysis, miRNA profiling was the only independent predictor of PFS in HG-NMIBC patients. In fact, in both intermediate and high risk EORTC score cohorts, patients with urinary expression for at least one member of the cluster displayed worse PFS (both log rank p < 0.001). On DCA, the benefit of using this prognostic model was evident in patients displaying a PFS > 20%. This cohort essentially includes all cases where the clinical decision making process between conservative options (intravesical BCG) and RC is controversial, being those patients with a PFS probability < 20% with a significant risk of disease progression, where the effectiveness of conservative treatments should be carefully balanced with the significant risk of progression.

The limitations of our study are the small sample size, and the tumor heterogeneity and multifocality that may hide variability in miRNA expression. Therefore, external validation of these findings with larger prospective studies is mandatory to confirm findings and to consider the clinical application of miRNA assays for diagnostic and prognostic purposes in HG-NMIBC.

Conclusions

Our study provides a proof for the potential use of cell-free miRNA biomarkers in clinical practice. We found that urinary expression of let-7c miRNA cluster could be a reliable non-invasive clinical tool to improve prognosis prediction of HG-NMIBC and to counsel patients according to an individual assessment of disease progression risk. Thus, from a practical point of view the importance of this potential assay is to intercept patients that undergo TUR of the tumor and are at risk of progression to address a more radical treatment.

Availability of data and materials

Data and materials are available for sharing if needed

Abbreviations

AUC:

Area Under Curve

BC:

Bladder Cancer

BCG:

Bacillus Calmette-Guerin

cDNA:

Complementary DNA

CI:

Interval of Confidence

C-index:

Concordance Index

CIS:

Carcinoma in Situ

Ct:

Cycle threshold

DCA:

Decision Curve Analysis

DCA:

Decision Curve Analysis

EC:

Extracellular

EORTC:

European Organization for Research and Treatment of Cancer

HC:

Healthy Controls

HG:

High Grade

HR:

Hazard Ratio

IC:

Intracellular

IQR:

Interquartile Range

LG:

Low Grade

miRNA:

microRNA

NMIBC:

Non-Muscle-Invasive Bladder Cancer

PFS:

Progression Free Survival

qRT-PCR:

quantitative Real time-Polymerase Chain Reaction

RC:

Radical Cystectomy

ROC:

Receiver Operating Characteristic

SFM:

Serum Free Medium

TUR:

Transurethral Resection

References

  1. Cambier S, Sylvester RJ, Collette L, Gontero P, Brausi MA, van Andel G, et al. EORTC Nomograms and risk groups for predicting recurrence, progression, and disease-specific and overall survival in non-muscle-invasive stage ta-T1 Urothelial bladder Cancer patients treated with 1-3 years of maintenance Bacillus Calmette-Guérin. Eur Urol. 2016;69(1):60–9.

    Article  Google Scholar 

  2. Siegel RL, Miller KD, Jemal A. Cancer statistics 2018. CA Cancer J Clin. 2018;68(1):7–30.

    Article  Google Scholar 

  3. Babjuk M, Böhle A, Burger M, Capoun O, Cohen D, Compérat EM, et al. EAU guidelines on non-muscle-invasive Urothelial carcinoma of the bladder: update 2016. Eur Urol. 2017;71(3):447–61.

    Article  Google Scholar 

  4. Peng Y, Croce CM. The role of MicroRNAs in human cancer. Signal Transduct Target Ther. 2016;1:15004.

    Article  Google Scholar 

  5. Lodewijk I, Dueñas M, Rubio C, Munera-Maravilla E, Segovia C, Bernardini A, et al. Liquid biopsy biomarkers in bladder Cancer: a current need for patient diagnosis and monitoring. Int J Mol Sci. 2018;19(9):2514.

    Article  Google Scholar 

  6. Larrea E, Sole C, Manterola L, Goicoechea I, Armesto M, Arestin M, et al. New Concepts in Cancer Biomarkers: Circulating miRNAs in Liquid Biopsies. Int J Mol Sci. 2016;17(5).

  7. Fendler A, Stephan C, Yousef GM, Kristiansen G, Jung K. The translational potential of microRNAs as biofluid markers of urological tumours. Nat Rev Urol. 2016;13(12):734–52.

    Article  CAS  Google Scholar 

  8. Yamada Y, Enokida H, Kojima S, Kawakami K, Chiyomaru T, Tatarano S, et al. MiR-96 and miR-183 detection in urine serve as potential tumour markers of urothelial carcinoma: correlation with stage and grade, and comparison with urinary cytology. Cancer Sci. 2011;102(3):522–9.

    Article  CAS  Google Scholar 

  9. Braicu C, Buiga R, Cojocneanu R, Buse M, Raduly L, Pop LA, et al. Connecting the dots between different networks: miRNAs associated with bladder cancer risk and progression. J Exp Clin Cancer Res. 2019;38:433.

    Article  Google Scholar 

  10. Yin XH, Jin YH, Cao Y, Wong Y, Weng H, Sun C, et al. Development of a 21-miRNA signature associated with the prognosis of patients with bladder cancer. Front Oncol. 2019;9:729.

    Article  Google Scholar 

  11. Hofbauer SL, de Martino M, Lucca I, Haitel A, Susani M, Shariat SF, et al. A urinary microRNA (miR) signature for diagnosis of bladder cancer. Urol Oncol. 2018;36(12):531.e1–8.

    Article  CAS  Google Scholar 

  12. Pardini B, Cordero F, Naccarati A, Viberti C, Birolo G, Oderda M, et al. MicroRNA profiles in urine by next-generation sequencing can stratify bladder cancer subtypes. Oncotarget. 2018;9(29):20658–69.

    Article  Google Scholar 

  13. Boyerinas B, Park SM, Hau A, Murmann AE, Peter ME. The role of let-7 in cell differentiation and cancer. Endocr Relat Cancer. 2010;17(1):F19–36.

    Article  CAS  Google Scholar 

  14. Tang H, Ma M, Dai J, Cui C, Si L, Sheng X, et al. miR-let-7b and miR-let-7c suppress tumourigenesis of human mucosal melanoma and enhance the sensitivity to chemotherapy. J Exp Clin Cancer Res. 2019;38:212.

    Article  Google Scholar 

  15. Pelosi A, Careccia S, Sagrestani G, Nanni S, Manni I, Schinzari V, et al. Dual promoter usage as regulatory mechanism of let-7c expression in leukemic and solid tumours. Mol Cancer Res. 2014;12(6):878–89.

    Article  CAS  Google Scholar 

  16. Nadiminty N, Tummala R, Lou W, Zhu Y, Shi XB, Zou JX, et al. MicroRNA let-7c is downregulated in prostate cancer and suppresses prostate cancer growth. PLoS One. 2012;7(3):e32832.

    Article  CAS  Google Scholar 

  17. Zhou H, Tang K, Xiao H, Zeng J, Guan W, Guo X, et al. A panel of eight-miRNA signature as a potential biomarker for predicting survival in bladder cancer. J Exp Clin Cancer Res. 2015;34:53.

    Article  Google Scholar 

  18. Vinall RL, Tepper CG, Ripoll AA, Gandour-Edwards RF, Durbin-Johnson BP, Yap SA, et al. Decreased expression of let-7c is associated with non-response of muscle-invasive bladder cancer patients to neoadjuvant chemotherapy. Genes Cancer. 2016;7(3–4):86–97.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. Kamat AM, Flaig TW, Grossman HB, Konety B, Lamm D, O'Donnell MA, et al. Expert consensus document: consensus statement on best practice management regarding the use of intravesical immunotherapy with BCG for bladder cancer. Nat Rev Urol. 2015;12(4):225–35.

    Article  Google Scholar 

  20. Bravo E, Calzolari A, De Castro P, Mabile L, Napolitani F, Rossi AM, et al. Developing a guideline to standardize the citation of bioresources in journal articles (CoBRA). BMC Med. 2015;13:33.

    Article  Google Scholar 

  21. Marchesi F, Regazzo G, Palombi F, Terrenato I, Sacconi A, Spagnuolo M, et al. Serum miR-22 as potential non-invasive predictor of poor clinical outcome in newly diagnosed, uniformly treated patients with diffuse large B-cell lymphoma: an explorative pilot study. J Exp Clin Cancer Res. 2018;37:95.

    Article  Google Scholar 

  22. Pospisilova S, Pazourkova E, Horinek A, Brisuda A, Svobodova I, Soukup V, et al. MicroRNAs in urine supernatant as potential non-invasive markers for bladder cancer detection. Neoplasma. 2016;63(5):799–808.

    Article  CAS  Google Scholar 

  23. Snowdon J, Boag S, Feilotter H, Izard J, Siemens DR. A pilot study of urinary microRNA as a biomarker for urothelial cancer. Can Urol Assoc J. 2013;7(1–2):28–32.

    Article  Google Scholar 

  24. Turchinovich A, Weiz L, Burwinkel B. Extracellular miRNAs: the mystery of their origin and function. Trends Biochem Sci. 2012;37(11):460–5.

    Article  CAS  Google Scholar 

  25. Healy NA, Heneghan HM, Miller N, Osborne CK, Schiff R, Kerin MJ. Systemic mirnas as potential biomarkers for malignancy. Int J Cancer. 2012;131(10):2215–22.

    Article  CAS  Google Scholar 

  26. Tanaka M, Oikawa K, Takanashi M, Kudo M, Ohyashiki J, Ohyashiki K, et al. Down-regulation of miR-92 in human plasma is a novel marker for acute leukemia patients. PLoS One. 2009;4:e5532.

    Article  Google Scholar 

  27. Ohshima K, Inoue K, Fujiwara A, Hatakeyama K, Kanto K, Watanabe Y, et al. Let-7 microRNA family is selectively secreted into the extracellular environment via exosomes in a metastatic gastric cancer cell line. PLoS One. 2010;5(10):e13247.

    Article  Google Scholar 

  28. Gao JM, Huang LZ, Huang ZG, He RQ. Clinical value and potential pathways of miR-183-5p in bladder cancer: a study based on miRNA-seq data and bioinformatics analysis. Oncol Lett. 2018;15(4):5056–70.

    PubMed  PubMed Central  Google Scholar 

  29. Sylvester RJ, van der Meijden AP, Oosterlinck W, Witjes JA, Bouffioux C, Denis L, et al. Predicting recurrence and progression in individual patients with stage ta T1 bladder cancer using EORTC risk tables: a combined analysis of 2596 patients from seven EORTC trials. Eur Urol. 2006;49(3):466–5.

    Article  Google Scholar 

  30. May M, Burger M, Brookman-May S, Stief CG, Fritsche HM, Roigas J, et al. EORTC progression score identifies patients at high risk of cancer-specific mortality after radical cystectomy for secondary muscle-invasive bladder cancer. Clin Genitourin Cancer. 2014;12(4):278–86.

    Article  Google Scholar 

Download references

Acknowledgments

Biobank IRCCS- Regina Elena National Cancer Institute (BBIRE), Rome, Italy. European Community for supporting A.B. Díaz Méndez

Funding

Supported by grants: Intramural research funds from Italian Ministry of Health (MoH), Ricerca Corrente (MGR, VP); MS is an Intramural funds MoH Fellow. GR is an “Associazione Mia Neri Foundation” fellow, supported by “Cassa Sovvenzioni e Risparmio Dipendenti Banca d’Italia”. ABD is a PhD student at University of Catania and an AiPBAND (H2020-MSCA-ITN-2017) fellow.

Author information

Authors and Affiliations

Authors

Contributions

MS performed experiments, with contribution of MC, GR, ABD, LC and VL; MCF, and IS performed and analyzed results of statistical and bioinformatics analyses; MC, MV, MCF, VP, RSF, UA, AB, AB, GT, RCF and GC collected urine and tissues and management of patient’s samples data base; MC, MS, MCF, GP, AR contributed to the interpretation of the results. MG, GS and MGR designed the study, critically interpreted results and wrote the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Giuseppe Simone or Maria Giulia Rizzo.

Ethics declarations

Ethics approval and consent to participate

The research protocol was reviewed and approved by the Ethics Committee (R.S 1179/19) of the “Regina Elena” National Cancer Institute, and informed consent was obtained from all participants to the study, in agreement with institutional guidelines. Demographic, clinical and pathological data were collected in a prospectively maintained database and were treated according to ethical and legal standards adopted by the Declaration of Helsinki.

Consent for publication

Informed consent for publication was obtained from all participants.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Manuela Spagnuolo, Present address: Curie Institute, PSL Research University, Sorbonne University, CNRS UMR 3244, F-75005 Paris, France.

Supplementary information

Additional file 1

: Figure S1. Association analysis between tissue let-7c cluster and tumor stage.

Additional file 2

: Figure S2. Let-7c cluster cut-off values for PFS probability.

Additional file 3

: Figure S3. Predictive nomogram of 12-months and 24 months of PFS Probability.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Spagnuolo, M., Costantini, M., Ferriero, M. et al. Urinary expression of let-7c cluster as non-invasive tool to assess the risk of disease progression in patients with high grade non-muscle invasive bladder Cancer: a pilot study. J Exp Clin Cancer Res 39, 68 (2020). https://doi.org/10.1186/s13046-020-01550-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13046-020-01550-w

Keywords