Main

Renal cell tumours (RCTs) account for 4% of all adult neoplasms and 90–95% of all tumours arising in the kidney, ranking 14th in incidence worldwide, with an age-standardised mortality rate of 1.6 out of 100 000 (Ferlay et al, 2010). Renal cell tumours are morphologically and genetically heterogeneous (Baldewijns et al, 2008), comprising both malignant tumours (which are subclassified mainly as clear-cell renal cell carcinoma (ccRCC, 70–80% of cases), papillary RCC (pRCC, 10–15% of cases) and chromophobe RCC (chRCC, 5–10% of cases) and benign tumours (among which oncocytoma is the most common subtype; Kovacs et al, 1997).

Because histological subtypes differ in clinical aggressiveness and prognosis (Amin et al, 2002; Ficarra et al, 2006), accurate classification is required for appropriate patient management. Moreover, most RCTs are clinically silent at their earliest stages, and 20–30% are diagnosed when metastatic spread has already occurred (Abrahams et al, 2004). Although widespread use of imaging techniques (mainly ultrasonography) has increased detection of suspicious renal masses, prompting new pre-operative diagnostic challenges as histological diagnosis using needle biopsy material meets with important limitations, hampering an accurate categorisation in many instances (Shen et al, 2012). In this setting, diagnosis relies mainly on morphologic features, which show some overlap among tumour subtypes. The discrimination between chRCC (mainly the eosinophilic variant) and oncocytoma is one of the most critical and, sometimes, difficult differential diagnosis. Although these tumour types share some morphologic characteristics, chRCC is a malignant neoplasm, capable of local invasion and metastatic spread, whereas oncocytoma is a benign tumour just requiring a more conservative management.

Over the years, several attempts have been made to assist morphology in differential diagnosis between chRCC from oncocytoma, including immunohistochemical profiles (Abrahams et al, 2004; Lin et al, 2006; Shen et al, 2012), histochemical stains (Skinnider and Jones, 1999) and gene expression analysis (Lee et al, 2011). However, sensitivity and specificity of those techniques are suboptimal and prompt the need for more accurate biomarkers. Interestingly, some recent studies have attempted to discriminate among RCC subtypes using microRNA (miRNA) expression analysis. Nevertheless, most of those studies have mainly dealt with ccRCC or, when the most relevant histological subtypes were included, only a limited number of samples of each subtype have been analysed, precluding a definitive conclusion about their accuracy (Nakada et al, 2008; Petillo et al, 2009; Jung et al, 2009; Fridman et al, 2010; Juan et al, 2010; Youssef et al, 2011; Zhao et al, 2013). MicroRNAs are small non-coding RNAs (22 nucleotides in length), which are involved in several essential biological processes such as cell differentiation, growth, apoptosis and proliferation (Esteller, 2011), and their deregulation has been implicated in tumorigenesis, including that of the kidney (Lu et al, 2005; Jung et al, 2009; Chow et al, 2010; White and Yousef, 2010). In addition to the differential expression patterns of miRNAs among RCT subtypes (Petillo et al, 2009; Jung et al, 2009; Juan et al, 2010; Valera et al, 2011; Youssef et al, 2011), altered miRNA expression might also provide relevant prognostic information (Neal et al, 2010).

In two recent reviews (Henrique et al, 2012; Jeronimo and Henrique, 2011), we found that five miRNAs (miR-21, miR-141, miR-155, miR-183 and miR-200b) had been reported as displaying diagnostic or prognostic value in RCT (Nakada et al, 2008; Petillo et al, 2009; Jung et al, 2009; Juan et al, 2010; Youssef et al, 2011). Thus, we aimed to confirm and extend those findings through expression analysis of a five miRNA panel in a single series of RCT, comprising the four major subtypes. We first assessed the expression levels of selected miRNAs in fresh-frozen tissues, focusing on the discrimination between benign and malignant tumours, as well as between chRCC and oncocytoma. In addition, the prognostic value of each miRNAs was determined. Finally, we validated our findings in a series of ex-vivo fine-needle aspiration biopsies from RCT to assess the feasibility of this approach as an ancillary tool in routine pathology.

Materials and methods

Clinical samples

A total of 130 fresh-frozen tissues were prospectively collected and included in this study, comprising 120 RCTs (30 cases of each of the four major subtypes (ccRCC, pRCC, chRCC and oncocytoma)) and 10 morphologically normal renal tissues (obtained from morphologic normal kidney tissue of patients subjected to nephrectomy due to upper urinary urothelial carcinoma; Table 1). In addition, a validation set comprising 60 ex-vivo fine-needle aspiration biopsies from RCT (15 of each subtype) was included. Samples from RCT patients were procured from patients diagnosed and treated at the Portuguese Oncology Institute – Porto (Portugal), between 2003 and 2007, who underwent partial or total nephrectomy, after obtaining informed consent. For each subtype, cases were consecutively selected until it reached 30 (for tissue samples) or 15 (for ex-vivo fine-needle aspiration biopsies) cases. This strategy was used to maximise the representation of the less common RCT types, thus ensuring that tumour heterogeneity in each subtype would be considered in the molecular analyses. Tumour tissue samples, obtained immediately after surgery, were snap-frozen, stored at −80 °C and subsequently cut in a cryostat for RNA extraction. Ex-vivo fine-needle aspiration biopsies of RCTs were obtained through 4–6 passes of a 23-gauge needle attached to a 10-ml syringe, then washed in PBS and stored at −80 °C until further use.

Table 1 Clinical and pathological features of patients included in this study, including the data for the two sets of samples (fresh-frozen tissues and ex-vivo biopsies)

Routine histopathological assessment of all surgical specimens, in formalin-fixed paraffin-embedded tissue, was performed by an expert uropathologist (RH) and included tumour classification (WHO), grading (Fuhrman) and staging (TNM; Eble et al, 2004). Relevant clinical data were collected from clinical charts. This study, as well as the use of samples and access to clinical data, was approved by the Institutional Review Board (Comissão de Ética para a Saúde) of the Portuguese Oncology Institute – Porto.

RNA extraction

Total RNA was extracted from fresh-frozen tissues and ex-vivo aspiration biopsies using Trizol reagent (Invitrogen, Carlsbad, CA, USA), according to manufacturer’s instructions. Briefly, 1500 μl of Trizol reagent was added to each 2 ml tube and samples were homogenised using a rotor shaker. Tubes were incubated for 5–10 min at room temperature and then 300 μl of chloroform (Merck, Darmstadt, Germany) were added. Regarding biopsies the protocol was similar, but the Trizol reagent and chloroform volumes were 500 and 200 μl, respectively. Tubes were vigorously hand-shaked for 15 s and incubated for 3 min at room temperature, followed by a 15-min centrifugation at 12 000 g at 4 °C. Next, the upper phase was collected. RNA was purified using the PureLink RNA Mini Kit (Invitrogen), according to manufacturer’s indications. RNA concentration and purity ratios were then evaluated using NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA). In addition, RNA quality was checked by electrophoresis in a 2% agarose gel.

Reverse transcription

Reverse transcription (RT) was performed using TaqMan MicroRNA Reverse Transcription Kit and Megaplex RT human pool A (Applied Biosystems, Foster City, CA, USA). The reaction mixture had a final volume of 12 μl and included the following: 3 μl of total RNA (750 ng), 1.6 μl of megaplex RT primers (10 × ), 0.4 μl of dNTPs with dTTP (100 mM), 3 μl of MultiScribe reverse transcriptase (50 U μl−1), 1.6 μl of 10 × RT buffer, 0.2 μl of RNase inhibitor (20 U μl−1) and 0.4 μl of nuclease-free water. Reactions were performed in PCR tubes according to the following conditions: 40 cycles at 16 °C for 2 min, 42 °C for 1 min and 50 °C for 1 s, with a final incubation at 85 °C for 5 min.

Quantitative real-time RT–PCR

Quantitative RT–PCR (qRT–PCR) was performed using TaqMan Small RNA Assays (Applied Biosystems) in a 7500 Real-Time PCR system (Applied Biosystems), according to the recommended protocol. Briefly, for each reaction 0.5 μl of TaqMan Small RNA Assay (20 × ), 0.75 μl of RT product, 5 μl of TaqMan Universal PCR Master Mix II no UNG (2 × ) and 3.75 μl of nuclease-free water were added. According to the manufacturer’s instructions, the protocol conditions were: 50 °C for 2 min, 95 °C for 10 min, followed by 40 cycles at 95 °C for 15 s and 60 °C for 1 min. Expression levels of the five selected miRNAs (hsa-miR-21: Tm000397; hsa-miR-141: Tm000483; hsa-miR-155: Tm002626; hsa-miR-183: Tm002269; and hsa-miR-200b: Tm002251) were assessed in triplicate for each sample and two water blanks were added to each plate as negative controls.

Results from the qRT–PCR were analysed using the 7500 Software version 2.0.5 (Applied Biosystems). Levels of miRNA expression were determined using the relative standard curve method (Biosystems, 2004). In each sample, the mean quantity of each miRNA was normalised to the mean quantity for the endogenous controls RNU48 and RNU6B, according to the following formula: miRNA expression=candidate miRNA expression mean quantity/((RNU48 mean quantity+RNU6B mean quantity)/2). Results were then multiplied by 10 000 for easier tabulation.

Statistical analysis

Differences in expression levels of the candidate miRNAs among the different histological subtypes were first analysed using a non-parametric Kruskal–Wallis test, followed by pairwise comparisons using non-parametric Mann–Whitney U-test, when appropriate. The relationship between miRNA expression and clinicopathological variables (gender, Fuhrman grade (recoded into two groups: grades 1–2 vs 3–4) and pathological stage (recoded into two groups: pT1–pT2 vs pT3–pT4)) was evaluated using Mann–Whitney U-test. Spearman’s non-parametric correlation tests were additionally carried out to ascertain correlations between age and miRNA expression levels. Receiver operator characteristic (ROC) curves were constructed by plotting the true-positive rate (sensitivity) against the false-positive rate (1−specificity) for each miRNA and for the best combination of miRNAs. The selection of the best miRNA panel was achieved using logistic regression, and the areas under the curve (AUCs) were calculated to assess the panel’s diagnostic performance. Disease-specific survival (DSS) curves (Kaplan–Meier with log-rank test) were computed for clinical variables (age, gender, histological subtype, Fuhrman grade and pathological stage) and miRNA expression levels. A Cox regression model comprising all significant variables (multivariate test) was computed to assess the relative contribution of each variable to the follow-up status. For the purpose of survival analyses, all cases were coded based on each miRNA expression levels, using the median value as the cut-off value. Statistical analysis was performed using SPSS for Windows, version 20.0 (SPSS, Chicago, IL, USA) and differences were considered statistically significant when P<0.05. For multiple comparisons, the P-value was adjusted according to Bonferroni’s method (i.e., the level of significance was adjusted to P<0.05/n, in which n represents the number of groups under comparison).

Results

MicroRNA expression levels and clinicopathological correlates

The relative expression levels of miR-21, miR-141, miR-155, miR-183 and miR-200b were determined in fresh-frozen tissues of 120 RCTs and 10 normal renal tissue samples. Relevant clinical and histopathological data are displayed in Table 1. No significant differences in age or gender between RCTs patients and normal tissue donors were apparent. No statistically significant associations were disclosed between miR expression levels and any of the clinicopathological features (age, gender, Fuhrman grade categories or pathological stage). Renal cell tumours showed significantly lower expression levels of miR-21, miR-141 and miR-200b compared with that of normal tissues (P<0.001 for all; Figure 1A and Supplementary Table 1). Moreover, expression levels of all candidate miRNAs differed significantly between benign and malignant RCTs. Oncocytomas displayed lower expression levels for all tested miRs, except miR-183 (Figure 1B and Supplementary Table 2).

Figure 1
figure 1

Distribution of miRNA expression levels in kidney tissues. (A) Normal vs tumour tissues. (B) Benign vs malignant tumour tissues. Statistically significant differences are represented as ***P<0.001, **P<0.01 and *P<0.02.

Although there was a wide expression range within the four RCT subtypes, with a significant degree of overlap, expression levels of all miRs differed significantly among them (P<0.001 for all, Kruskal–Wallis test; Table 2). Pairwise comparisons are shown in Table 3 and graphically illustrated in Figure 2. In general, oncocytomas displayed the lower miR expression levels, significantly differing from pRCC or ccRCC regarding four miRs (miR-21, miR-155, miR-183 and miR-200b), and from chRCC in two miRs (miR-141 and miR-200b). Interestingly, ccRCC and pRCC only differed for miR-155 expression levels, whereas chRCC differed from ccRCC and pRCC for miR-21, miR-141 and miR-155 expression levels. In addition, miR-183 expression levels were also different between chRCC and ccRCC (Table 3). Thus, reduced expression of miR-200b surfaced as the most discriminative feature between oncocytomas and RCCs.

Table 2 Distribution of microRNA expression levels among different histological subtypes in fresh-frozen tissues
Table 3 Comparison of microRNA expression among renal cell tumour subtypes in fresh-frozen tissues
Figure 2
figure 2

Distribution expression levels of ( A ) miR-21, ( B ) miR-141, ( C ) miR-155, ( D ) miR-183 and ( E ) miR-200b according with the histological subtype of RCTs. Statistically significant differences are represented as ***P<0.001, **P<0.003 and *P<0.0125.

Diagnostic performance of miRNA expression levels in tissue samples

Performance of the five studied miRs was assessed in three different settings: identification of RCTs (vs normal renal tissue), discrimination of malignant from benign tumours and distinction of chRCC from oncocytoma. For that purpose, the cut-off value corresponded to the best performance of each miRNA according to the respective ROC curve analysis. Validity and information estimates for each marker and for the best combination of markers are displayed in Table 4. Receiver operator characteristic curve analysis showed that a panel comprising expressions of miR-141 and miR-200b allowed for the discrimination between RCT and normal renal tissue with 99.2% sensitivity and 100% specificity, corresponding to an AUC of 0.991. In addition, the same panel allowed for the differentiation between benign and malignant tumours with 85.6% sensitivity and 100% specificity, displaying an AUC of 0.912. Furthermore, expression levels of miR-141 and miR-200b also distinguished chRCC from oncocytoma with 90% sensitivity and 100% specificity, corresponding to an AUC of 0.90 (Figure 3A and B).

Table 4 Validity estimates for each tested miR and for the best combination of miRs in each diagnostic setting, in fresh-frozen tissues
Figure 3
figure 3

ROC curves. ROC curves evaluating the performance of the gene panel (miR-141 and miR-200b) as a biomarker for malignant renal tumours (A and C) and as a biomarker of chRCC (B and D). (A and B) Performed in kidney tissue samples; (C and D) performed in ex-vivo aspiration renal biopsies.

Survival analysis

The median follow-up of this series of RCT patients was 65 months (range: 1–120 months). A total of 12 patients (13.3%) have died from RCC during this period. Disease-specific survival analysis showed that tumour subtype ccRCC or pRCC and higher pathological tumour stage (pT3–T4) were significantly associated with worse outcome (P=0.011 and P<0.001, respectively; Figure 4A and B). Although age at diagnosis over 62 years was associated with worse DSS (P=0.035), gender and Fuhrman grade did not disclose any prognostic value within the available follow-up time. Concerning miRNA expression levels, miR-200b and miR-183 did not exhibit any prognostic value. However, higher expression levels of miR-21 and miR-155, and lower expression levels of miR-141 were associated with worse DSS (P=0.006, P=0.037 and P=0.024, respectively; Figure 4C–E). However, in multivariate analysis only pathological stage independently predicted prognosis, whereas miRNA expression levels did not retain an independent prognostic value (Supplementary Table 3).

Figure 4
figure 4

Disease-specific survival according to pathological and molecular parameters. (A) Histopathological classification; (B) pathological stage; (CE) miR expression levels.

Validation of the miRNA panel in ex-vivo aspiration biopsies

The two best-performing miRNA in tissue samples, miR-141 and miR-200b, were then selected for analysis in ex-vivo samples. This set comprised 60 ex-vivo fine-needle aspiration biopsies. Relevant clinical and histopathological data are summarised in Table 1 and the relative expression levels for each miR are depicted in Supplementary Table 4 and Supplementary Figures S1 and S2.

Remarkably, expression levels of this panel of miRNAs were able not only to distinguish benign from malignant RCT with 73.3% sensitivity and a 100% specificity (AUC of 90.4%), but also oncocytoma from chRCC with 100% sensitivity and 100% specificity (AUC of 100%; Figure 3C and D; Table 5).

Table 5 Validity estimates for each tested miR and for the best combination of miRs in each diagnostic setting, in ex-vivo aspiration biopsies

Discussion

In this study, we aimed to define a small set of miRs that might allow for accurate identification of RCTs, as well as for discrimination between oncocytoma and RCCs, especially chRCC. This would be of clinical relevance, as diagnostic workup of suspicious renal masses incidentally found by abdominal ultrasonography is increasingly more frequent and demanding. Indeed, each RCT subtype displays quite dissimilar clinical behaviour, ranging from totally benign to overtly malignant, and successful pretherapeutic cytological or histological assessment is limited (Amin et al, 2002; Ficarra et al, 2006). Only a few studies addressed the use of miRNA expression as biomarkers for RCTs detection, and these have been mainly restricted to the ccRCC subtype, or have only analysed a very limited number of samples (Jung et al, 2009; Petillo et al, 2009; Juan et al, 2010; Youssef et al, 2011; Redova et al, 2012; Zhao et al, 2013). After an extensive review of published literature, we selected five miRNAs (miR-21, miR-141, miR-155, miR-183 and miR-200b) with putative diagnostic and prognostic value (Jeronimo and Henrique, 2011; Henrique et al, 2012), and tested them in a relatively large set of tissue samples that comprised the major histological subtypes. To ascertain their clinical and pathological relevance, a validation study was subsequently performed in a set of ex-vivo fine-needle aspiration biopsies.

Of the five miRs tested, three (miR-21, miR-141 and miR-200b) were significantly downregulated in RCTs compared with normal renal tissue. In previous reports, miR-21 was found to be upregulated in RCT (Juan et al, 2010; Faragalla et al, 2012; Zaman et al, 2012), which apparently contradicts our results. However, in those studies, normal renal tissue was obtained from nephrectomy specimens harbouring RCT, which did not occur in our study. This is an important issue, as we have previously shown that morphologically normal renal tissue from kidneys harbouring RCT display epigenetic alterations in line with the respective tumours (Costa et al, 2007). Remarkably, variations in miR-21 expression among RCT subtypes observed in our study matches that reported by Faragalla et al (2012), with ccRCC depicting the highest median levels, followed by pRCC, chRCC and oncocytoma. Indeed, only miR-21 expression levels of ‘normal renal tissue’ are notably different between our results and their study (Faragalla et al, 2012). These findings prompt the need for an adequate definition of ‘normal tissue’, as the interpretation of results in tumours might be considerably biased.

Concerning miR-141, our results corroborate those of two previous reports (Nakada et al, 2008; Fridman et al, 2010). Thus, higher miR-141 expression levels seem to be a hallmark of chRCC and might constitute a valuable biomarker for discrimination from oncocytoma. Strikingly, a miRNA profiling of ccRCC also identified miR-141 (and 200b) as being downregulated in ccRCC, although with concurrent upregulation of miR-155 (Jung et al, 2009). These results are in line with ours, as we found that the highest miR-155 expression levels in ccRCC and pRCC significantly differed from those of oncocytoma and chRCC. Our findings concerning miR-200b mirror those of Youssef et al (2011), although in a smaller data set. Interestingly, in our larger data set we were able to demonstrate that miR-200b expression levels were significantly lower in oncocytomas compared with all RCC subtypes. Overall, the comparisons of miR expression levels among RCT subtypes also denote the common origin (segment of the nephron) of ccRCC and pRCC on one hand, and of chRCC and oncocytoma on the other hand, emphasising the importance of searching for discriminative biomarkers, which might enable accurate identification of each RCT subtype.

Interestingly, a panel comprising miR-141 and miR-200b demonstrated the best performance in frozen-tissue samples, displaying AUC values ranging from 90.0 to 99.1. Although these results are interesting per se, its clinical usefulness depends on the possibility of using it in diagnostic samples. For that purpose, we further validated this biomarker panel in a set of fine-needle aspiration biopsies performed ex vivo. Although this procedure is not completely equivalent to an imaging-guided diagnostic fine-needle aspiration biopsy performed in a patient (which may yield lower amounts of tumour cells), it is, nonetheless, the best approximation without jeopardising patients’ diagnosis. On the other hand, because the nephrectomy specimen is already available its histopathological characterisation is guaranteed, whereas diagnostic biopsies may not be followed by surgical excision, thus precluding accurate tumour classification for comparison purposes. Remarkably, the biomarker panel performance in ex-vivo biopsies was comparable to that of fresh-frozen tissues. To the best of our knowledge, this is the first attempt to demonstrate the feasibility of using miRs as tumour biomarkers in renal tumour biopsies, and may thus constitute a significant step forward in the development of epigenetic-based biomarkers for management of RCC suspects.

The clinical significance of our findings could be extended if miRNA expression levels might convey prognostic information. Thus, we performed DSS analysis using expression levels determined in fresh-frozen-tissue samples. As expected, tumour subtype and pathological stage were of prognostic value in univariate analysis, although only the later showed independent prognostic value in multivariate analysis. Remarkably, miR-21, miR-141 and miR-155 expression levels also displayed prognostic significance in RCC, although only in univariate analysis. A possible explanation for these findings may lie in the association between specific miR expression levels and tumour subtypes. Indeed, whereas for miR-21 and miR-155 the association with poorer DSS was observed for higher (> median) expression levels, the opposite was verified for miR-141. Interestingly, higher miR-21 and miR-155 expression levels and lower miR-141 expression levels were associated with pRCC and ccRCC subtypes, which displayed the worse prognosis compared with that of chRCC. The fact that tumour subtype did not surfaced as independent prognostic parameter for DSS in multivariate analysis is most likely due to the association between tumour subtype and pathological stage, as pT3–4 tumours were mostly of pRCC or ccRCC subtype. Our findings concerning miR-21 and miR-141 are corroborated by previous reports, although with generally smaller patient cohorts (Jung et al, 2009; Faragalla et al, 2012; Zaman et al, 2012). In addition, the prognostic value of miR-155 expression levels has been reported for breast cancer (Song et al, 2012) and non-small cell lung cancer (Yanaihara et al, 2006; Yang et al, 2013), whereas miR-21 and miR-141 expression seem to be of prognostic significance in non-small cell lung cancer (Yanaihara et al, 2006; Yang et al, 2013) and colon cancer (Cheng et al, 2011), respectively.

The aforementioned association of specific miRs altered expression and RCT subtype might also provide clues concerning the cause of miR dysregulation. Renal cell tumour subtypes display characteristic chromosomal aberrations, including whole or partial deletions and duplications (Baldewijns et al, 2008). Strikingly, some of those alterations might explain the altered pattern of miR expression. For instance, miR-200b is mapped at 1p36.33 and loss of 1p or of the whole chromosome 1 is frequently observed in oncocytoma and chRCC. On the other hand, miR-21 and miR-155 are mapped at 17q23.1 and 21q21.2–21.3, which are frequently lost chromosomal regions in chRCC. Conversely, pRCC, which commonly show gain of chromosome 17, are among the RCT subtypes with higher miR-21 expression levels. However, other variations in miR expression might not be explained by chromosomal-level alterations and the respective cause(s) remain to be investigated.

Conclusions

Herein we demonstrate that expression levels of a panel of two miRNAs (miR-141 or miR-200b) allows for accurate distinction of normal kidney from RCT tissue samples, as well as for accurate discrimination among RCT subtypes, including the separation of benign from malignant RCT. Furthermore, the selected miR panel is able to convey prognostic information, although not independent of tumour subtype or pathological stage. Importantly, the same panel displays an impressive performance for accurate detection of RCC in clinical samples obtained from ex-vivo fine-needle aspiration biopsies, demonstrating the feasibility of this approach in routine diagnostic practice.