Research | Open | Published:
Upregulation of microRNA-106b is associated with poor prognosis in hepatocellular carcinoma
Diagnostic Pathologyvolume 9, Article number: 226 (2014)
MicroRNA-106b (miR-106b) is a member of the miR-106b ~ 25 cluster. It has been reported that miR-106b acts as an oncogene and is upregulated in many human cancers. However, the prognostic value of miR-106b in hepatocellular carcinoma (HCC) remains unclear. The aim of this study was to investigate the clinical significance of miR-106b expression in HCC.
We determined the expression level of miR-106b in 104 cases of paired HCC and adjacent non-tumor tissues by quantitative real-time PCR (qRT-PCR). The correlation between miR-106b expression and prognosis of HCC was studied by univariate and multivariate analysis. Multivariate analysis of the prognostic factors was performed with Cox proportional hazards model.
MiR-106b expression was significantly upregulated in as high as 76.0% of HCC tissues, compared with their non-tumor counterparts (P < 0.001). High miR-106b expression was significantly associated with large tumor size (P = 0.019) and vascular invasion (P = 0.016). Kaplan-Meier analysis showed that patients with high miR-106b expression had a worse overall survival than patients with low miR-106b expression (log-rank P = 0.004). The multivariate Cox regression analysis indicated that miR-106b expression was an independent prognostic factor for overall survival (HR, 2.002; 95% CI, 1.130-6.977; P = 0.027).
Our data indicated that miR-106b expression was significantly upregulated in HCC and could serve as a potential unfavorable prognostic biomarker.
The virtual slide(s) for this article can be found here: http://www.diagnosticpathology.diagnomx.eu/vs/13000_2014_226
Hepatocellular carcinoma (HCC) is one of the most common cancers worldwide with poor prognosis . Accurate prediction of prognosis and patient stratification are crucial for guiding patients’ personalized clinical treatment. These are currently performed by clinical and/or pathological staging systems -. Recently important progress has been made with using comprehensive approaches to identify the molecular diversity in HCC. A large number of genetic and epigenetic abnormalities were found during the process of HCC development -. Consequently, many new prognostic biomarkers of HCC have been described ,. Addition of new biomarkers in the current staging systems is likely to improve the prognostic prediction of HCC patients and identify specific subgroups of tumor.
MicroRNAs (miRNAs) are highly conserved, small non-coding RNA molecules of approximately 22 nucleotides in length that function as posttranscriptional gene regulators ,. Many studies have shown that miRNAs are directly involved in cancer initiation and progression. Indeed, miRNAs could act as oncogenes or tumor suppressors ,. In addition, some miRNA deregulations were found to be associated with specific tumor phenotype . These observations suggest that miRNA could be used as putative biomarkers to classify tumors.
MiR-106b is a member of the miR-106b ~ 25 cluster, which consists of miR-106b, miR-93 and miR-25, that is encoded within intron 13 of the minichromosome maintenance complex component 7 (MCM7) gene on chromosome 7q22.1 . MiR-106b has been reported as an oncogene and upregulated in many human cancers, including stomach, prostate and kidney -. Previous studies suggest that miR-106b also plays an important role in hepatocarcinogenesis -. In our previous study, using miRNA array, we also found miR-106b is one of the upregulated miRNAs in HCC . However, the prognostic value of miR-106b in HCC has not been fully clarified yet.
In the present study, we investigated the expression levels of miR-106b in resected HCC specimens and adjacent non-tumor tissues. We also analyzed the association of miR-106b expression with clinicopathological characteristics and overall survival of the patients and determined whether miR-106b is potentially predictive of prognosis in HCC.
Patients and tissue samples
A total of 104 pairs of tumor and matched adjacent non-tumor liver tissues were collected from patients with HCC who underwent hepatic resections between July 2003 and July 2007 at Sun Yat-sen University Cancer Center. These patients were selected for hepatic resection in accordance with our defined criteria . For all cases, tissue samples were obtained immediately after resection and were snap-frozen in liquid nitrogen and stored at −80°C until use. All tissues were archived in the institution’s liver tumor bank. The diagnosis was confirmed histologically in all cases, based on detailed examination of sections stained with H&E. None of the patients had received any other therapy, including percutaneous ablation and chemo-embolization, before surgery. The research protocol was approved by the Institutional Review Board of Sun Yat-sen University Cancer Center. Written informed consent was obtained from every patient to participate in this study.
Clinical follow-up was performed for all the patients. All patients were monitored by physical examination, serum alpha-fetoprotein (AFP), ultrasonography, and x-ray every 1–3 months in the first year, and every 3–6 months thereafter for surveillance of recurrence or metastases. Computed tomography and/or magnetic resonance imaging and/or positron emission tomography were used whenever needed to confirm the clinical findings. The clinical end point in the present study was overall survival. Overall survival was defined as the interval from curative surgery to the date of death or the date of last contact if the patient was still alive. Survival information of all patients was updated by telephone visit and questionnaire letters.
Quantitative real-time RT-PCR (qRT-PCR)
To detect miRNA expression by qRT-PCR, total RNA was extracted from tissue samples using the Trizol reagent (Invitrogen, Carlsbad, CA, USA). qRT-PCR of miR-106b was performed using Bulge-Loop™ miRNA qRT-PCR kits (Ribobio, Guangzhou, China) according to the manufacturer’s instructions. Briefly, reverse transcription was performed using a specific miR-106b stem-loop primer and the reverse primer for U6 small nuclear RNA, followed by real-time PCR with a miR-106b-specific forward primer and a universal reverse primer using an ABI-Prism 7900HT system (Applied Biosystems, Foster City, CA). All samples were carried out in triplicate. U6 was used as an endogenous control. The expression level of miR-106b for each sample was calculated, represented by the ΔCt value (Ct of miR-106b - Ct of U6). The relative miR-106b expression levels in paired tissues that were collected from the same patients were analyzed by the 2-ΔΔCt method, represented by the -ΔΔCt value [−(ΔCt of tumor tissues - ΔCt of non-tumor tissues)]. The -ΔΔCt value represents the value of log2 T/NT for each patient.
Statistical comparisons were performed using either the Student’s t-test, chi-square test or Fisher’s exact test, as appropriate. In a univariate analysis, survival curves were assessed by the Kaplan-Meier method, and differences between curves were analyzed using the log-rank test. Covariates with P values < 0.05 in the univariate analysis were subjected to multivariate analysis. Multivariate Cox regression models were constructed to estimate the hazard ratios (HRs) of independent factors for survival after controlling for potential confounding factors. All statistical analyses were carried out using the SPSS (Statistical Package for the Social Sciences) 13.0 software (SPSS, Chicago, IL) and GraphPad Prism 5.0 software (GraphPad Software, San Diego, CA). Two-tailed P values < 0.05 were considered statistically significant.
MiR-106b expression in HCC
We performed quantitative real-time RT-PCR to examine the miR-106b expression level in 104 pairs of HCC and adjacent non-tumor tissues. As shown in Figure 1, the expression level of miR-106b in HCC tissues was significantly higher than that in adjacent non-tumor tissues after normalization (median fold change of T/NT = 2.27, P < 0.001). MiR-106b was upregulated in 79 of 104 patients (76.0%) totally (Figure 2). The median fold change of miR-106b was used as a cutoff value to divide all 104 patients into two groups: the low expression group (n = 52) and the high expression group (n = 52).
Correlation of miR-106b expression with clinicopathological characteristics of HCC patients
We next analyzed the correlation between miR-106b expression and the clinicopathological characteristics of HCC, including patients’ age, gender, HBsAg, Child-Pugh classification, serum AFP level, tumor size, tumor number, vascular invasion, histological grade (Edmondson-Steiner) and TNM stage. As summarized in Table 1, miR-106b expression was significantly higher in HCC patients with large tumor than those with small tumor (P = 0.019). Additionally, a Pearson correlation analysis also showed that the miR-106b level and tumor size were positively correlated (r = 0.2894, P = 0.0029; Figure 3). Moreover, miR-106b was expressed at significantly higher levels in patients with vascular invasion than in patients without vascular invasion (P = 0.016). However, no significant correlation was observed between miR-106b expression and other clinicopathological characteristics.
Prognostic analysis of miR-106b expression and clinicopathological factors
The association between miR-106b expression and prognosis of HCC patients was investigated by Kaplan-Meier analysis and log-rank test. As shown in Figure 4, HCC patients with high miR-106b expression had shorter overall survival than those with low miR-106b expression. The 1, 3, and 5-year overall survival rate in the high expression group was 84.0%, 51.6%, and 36.5%, respectively, compared with 84.4%, 60.2%, and 56.2%, respectively, in the low expression group (log-rank test, P = 0.004).
Univariate analysis demonstrated that Serum AFP level (P = 0.041), tumor size (P < 0.001), vascular invasion (P = 0.001), histological grade (P = 0.048) and TNM stage (P = 0.004), and miR-106b expression (P = 0.004) were significantly associated with overall survival of HCC patients (Table 2). No significant associations were found for age at diagnosis, gender, HBsAg status, Child-Pugh classification, and tumor number. Multivariate analysis using the Cox proportional hazards model for variables that were significant in the univariate analysis showed that tumor size (P = 0.022), vascular invasion (P = 0.028) and miR-106b expression (P = 0.027) were independent prognostic factors for patients with HCC (Table 2).
Discussion and conclusions
MiRNAs are considered to be ideal biomarkers for the diagnosis and prognosis of cancer because they are easy to detect and strongly associated with clinical prognoses -. So far, increasing evidence shows that more and more deregulated miRNAs are involved in HCC carcinogenesis and patient prognosis ,,-,. However, due to the complicated and diverse roles of miRNAs in HCC, clarifying the clinical significance and function of an individual miRNA may enable deeper insights into the complexity of HCC progression and disease management.
MiR-106b is a member of the oncogenic miR-106b-25 cluster . The biological function of miR-106b in HCC has been investigated in previous studies. It has been demonstrated that ectopic expression of miR-106b can enhance the proliferation and anchorage-independent growth of HCC cells, whereas inhibition of miR-106b had the opposite effect ,. It has also been shown that miR-106b contributed to metastasis by activating the EMT process and promoting cell migration in vitro and metastasis in vivo in HCC, indicating that upregulation of miR-106b may correlate with HCC progression . However, the clinical significance and prognostic value of miR-106b in HCC remain unclear.
In the present study, we analyzed the clinical relevance of miR-106b expression in HCC patients. Our qRT-PCR data confirmed that miR-106b expression was upregulated in tumor tissues compared with the adjacent non-tumor tissues, which validated our previous microarray high-throughput profiling results . These results are also consistent with previous findings in other studies ,. In addition, we correlated clinicopathological characteristics of the patients with miR-106b expression and observed that overexpression of miR-106b correlated with large tumor size and vascular invasion, which strongly indicated that this miRNA plays oncogenic roles in HCC, including promoting cell growth, cell invasion and tumor metastasis.These observations could be attributed to the miR-106b target genes, which include p21/CDKN1A, adenomatous polyposis coli (APC), transforming growth factor-β type II receptor (TGF-βRII), and RhoA and RhoC ,,,. Hence, by regulating its target genes, miR-106b could promote HCC cell cycle progression, cell proliferation, and cell migration and invasion.
More importantly, we demonstrated that miR-106b expression was significantly associated with overall survival of patients with HCC. Kaplan-Meier analysis showed that patients with high miR-106b expression level had a significantly shorter overall survival than those with low miR-106b expression level. Furthermore, multivariate Cox analysis proved that miR-106b was a prognostic factor independent of adjusted well-known prognostic variables for HCC including serum AFP level, tumor size, vascular invasion, histological grade, and TNM stage. Thus, miR-106b could be used as a potential prognostic biomarker in addition to other known prognostic indicator, in order to identify a subgroup of patients who have higher risk of death, thus, should receive monitoring more frequently and effective adjuvant treatment. Many studies have shown that miRNA profiles are maintained and are not significantly affected in formalin-fixed paraffin-embedded tissues ,. Therefore, although our study is based on frozen tissue samples, it offers a significant opportunity for using archived samples to detect changes in miR-106b expression levels by using in situ hybridization methodology with miRNA localization in routine clinical setting . In summary, our data indicated that miR-106b expression was significantly upregulated and associated with poor prognosis in HCC. MiR-106b was identified as an independent biomarker for predicting the clinical prognosis of HCC patients. Further studies are needed to validate the prognostic value of miR-106b expression in other cohorts.
Quantitative real-time PCR
Minichromosome maintenance complex component 7
Adenomatous polyposis coli
Transforming growth factor-β type II receptor
Siegel R, Ma J, Zou Z, Jemal A: Cancer statistics, 2014. CA Cancer J Clin. 2014, 64: 9-29. 10.3322/caac.21208.
Sobin LH, Compton CC: TNM seventh edition: what’s new, what’s changed: communication from the International Union against cancer and the American joint committee on cancer. Cancer. 2010, 116: 5336-5339. 10.1002/cncr.25537.
Llovet JM, Bruix J: Prospective validation of the Cancer of the Liver Italian Program (CLIP) score: a new prognostic system for patients with cirrhosis and hepatocellular carcinoma. Hepatology. 2000, 32: 679-680. 10.1053/jhep.2000.16475.
Forner A, Reig ME, de Lope CR, Bruix J: Current strategy for staging and treatment: the BCLC update and future prospects. Semin Liver Dis. 2010, 30: 61-74. 10.1055/s-0030-1247133.
Schmilovitz-Weiss H, Tobar A, Halpern M, Levy I, Shabtai E, Ben-Ari Z: Tissue expression of squamous cellular carcinoma antigen and Ki67 in hepatocellular carcinoma-correlation with prognosis: a historical prospective study. Diagn Pathol. 2011, 6: 121-10.1186/1746-1596-6-121.
Guo X, Xiong L, Sun T, Peng R, Zou L, Zhu H, Zhang J, Li H, Zhao J: Expression features of SOX9 associate with tumor progression and poor prognosis of hepatocellular carcinoma. Diagn Pathol. 2012, 7: 44-10.1186/1746-1596-7-44.
Su ZX, Zhao J, Rong ZH, Geng WM, Wu YG, Qin CK: Upregulation of microRNA-25 associates with prognosis in hepatocellular carcinoma. Diagn Pathol. 2014, 9: 47-10.1186/1746-1596-9-47.
Wang WY, Zhang HF, Wang L, Ma YP, Gao F, Zhang SJ, Wang LC: High expression of microRNA-130b correlates with poor prognosis of patients with hepatocellular carcinoma. Diagn Pathol. 2014, 9: 160-10.1186/s13000-014-0160-5.
Chaiteerakij R, Addissie BD, Roberts LR: Update on biomarkers of hepatocellular carcinoma.Clin Gastroenterol Hepatol 2013, S1542-3565(13)01789-8.,
Zhao YJ, Ju Q, Li GC: Tumor markers for hepatocellular carcinoma. Mol Clin Oncol. 2013, 1: 593-598.
Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116: 281-297. 10.1016/S0092-8674(04)00045-5.
He L, Hannon GJ: MicroRNAs: small RNAs with a big role in gene regulation. Nat Rev Genet. 2004, 5: 522-531. 10.1038/nrg1379.
Calin GA, Croce CM: MicroRNA signatures in human cancers. Nat Rev Cancer. 2006, 6: 857-866. 10.1038/nrc1997.
Ling H, Zhang W, Calin GA: Principles of microRNA involvement in human cancers. Chin J Cancer. 2011, 30: 739-748. 10.5732/cjc.011.10243.
Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, Sweet-Cordero A, Ebert BL, Mak RH, Ferrando AA, Downing JR, Jacks T, Horvitz HR, Golub TR: MicroRNA expression profiles classify human cancers. Nature. 2005, 435: 834-838. 10.1038/nature03702.
Petrocca F, Vecchione A, Croce CM: Emerging role of miR-106b-25/miR-17-92 clusters in the control of transforming growth factor beta signaling. Cancer Res. 2008, 68: 8191-8194. 10.1158/0008-5472.CAN-08-1768.
Li F, Liu J, Li S: MicorRNA 106b approximately 25 cluster and gastric cancer. Surg Oncol. 2013, 22: e7-e10. 10.1016/j.suronc.2013.01.003.
Hudson RS, Yi M, Esposito D, Glynn SA, Starks AM, Yang Y, Schetter AJ, Watkins SK, Hurwitz AA, Dorsey TH, Stephens RM, Croce CM, Ambs S: MicroRNA-106b-25 cluster expression is associated with early disease recurrence and targets caspase-7 and focal adhesion in human prostate cancer. Oncogene. 2013, 32: 4139-4147. 10.1038/onc.2012.424.
Slaby O, Jancovicova J, Lakomy R, Svoboda M, Poprach A, Fabian P, Kren L, Michalek J, Vyzula R: Expression of miRNA-106b in conventional renal cell carcinoma is a potential marker for prediction of early metastasis after nephrectomy. J Exp Clin Cancer Res. 2010, 29: 90-10.1186/1756-9966-29-90.
Li Y, Tan W, Neo TW, Aung MO, Wasser S, Lim SG, Tan TM: Role of the miR-106b-25 microRNA cluster in hepatocellular carcinoma. Cancer Sci. 2009, 100: 1234-1242. 10.1111/j.1349-7006.2009.01164.x.
Yau WL, Lam CS, Ng L, Chow AK, Chan ST, Chan JY, Wo JY, Ng KT, Man K, Poon RT, Pang RW: Over-expression of miR-106b promotes cell migration and metastasis in hepatocellular carcinoma by activating epithelial-mesenchymal transition process. PLoS One. 2013, 8: e57882-10.1371/journal.pone.0057882.
Shen G, Jia H, Tai Q, Li Y, Chen D: miR-106b downregulates adenomatous polyposis coli and promotes cell proliferation in human hepatocellular carcinoma. Carcinogenesis. 2013, 34: 211-219. 10.1093/carcin/bgs320.
Wei R, Huang GL, Zhang MY, Li BK, Zhang HZ, Shi M, Chen XQ, Huang L, Zhou QM, Jia WH, Zheng XF, Yuan YF, Wang HY: Clinical significance and prognostic value of microRNA expression signatures in hepatocellular carcinoma. Clin Cancer Res. 2013, 19: 4780-4791. 10.1158/1078-0432.CCR-12-2728.
Li B, Liu W, Wang L, Li M, Wang J, Huang L, Huang P, Yuan Y: CpG island methylator phenotype associated with tumor recurrence in tumor-node-metastasis stage I hepatocellular carcinoma. Ann Surg Oncol. 2010, 17: 1917-1926. 10.1245/s10434-010-0921-7.
Calin GA, Ferracin M, Cimmino A, Di Leva G, Shimizu M, Wojcik SE, Iorio MV, Visone R, Sever NI, Fabbri M, Iuliano R, Palumbo T, Pichiorri F, Roldo C, Garzon R, Sevignani C, Rassenti L, Alder H, Volinia S, Liu CG, Kipps TJ, Negrini M, Croce CM: A MicroRNA signature associated with prognosis and progression in chronic lymphocytic leukemia. N Engl J Med. 2005, 353: 1793-1801. 10.1056/NEJMoa050995.
Yanaihara N, Caplen N, Bowman E, Seike M, Kumamoto K, Yi M, Stephens RM, Okamoto A, Yokota J, Tanaka T, Calin GA, Liu CG, Croce CM, Harris CC: Unique microRNA molecular profiles in lung cancer diagnosis and prognosis. Cancer Cell. 2006, 9: 189-198. 10.1016/j.ccr.2006.01.025.
Kong YW, Ferland-McCollough D, Jackson TJ, Bushell M: microRNAs in cancer management. Lancet Oncol. 2012, 13: e249-e258. 10.1016/S1470-2045(12)70073-6.
Giordano S, Columbano A: MicroRNAs: new tools for diagnosis, prognosis, and therapy in hepatocellular carcinoma?. Hepatology. 2013, 57: 840-847. 10.1002/hep.26095.
Kan T, Sato F, Ito T, Matsumura N, David S, Cheng Y, Agarwal R, Paun BC, Jin Z, Olaru AV, Selaru FM, Hamilton JP, Yang J, Abraham JM, Mori Y, Meltzer SJ: The miR-106b-25 polycistron, activated by genomic amplification, functions as an oncogene by suppressing p21 and Bim. Gastroenterology. 2009, 136: 1689-1700. 10.1053/j.gastro.2009.02.002.
Li J, Smyth P, Flavin R, Cahill S, Denning K, Aherne S, Guenther SM, O’Leary JJ, Sheils O: Comparison of miRNA expression patterns using total RNA extracted from matched samples of formalin-fixed paraffin-embedded (FFPE) cells and snap frozen cells. BMC Biotechnol. 2007, 7: 36-10.1186/1472-6750-7-36.
Klopfleisch R, Weiss AT, Gruber AD: Excavation of a buried treasure–DNA, mRNA, miRNA and protein analysis in formalin fixed, paraffin embedded tissues. Histol Histopathol. 2011, 26: 797-810.
Renwick N, Cekan P, Bognanni C, Tuschl T: Multiplexed miRNA fluorescence in situ hybridization for formalin-fixed paraffin-embedded tissues. Methods Mol Biol. 2014, 1211: 171-187. 10.1007/978-1-4939-1459-3_14.
This work was supported in part by the National Natural Science Foundation of China (81101863 and 81172344), Guangdong Department of Science & Technology Translational Medicine Center grant (2011A080300002), Guangdong Natural Science Foundation (S2011040004570), and the Fundamental Research Funds for the Central Universities (11ykpy60).
The authors declare that they have no competing interests.
B-KL carried out the quantitative real-time RT-PCR assay and drafted the manuscript; P-ZH, G-LH, Y-DL and JH participated in the real-time quantitative RT-PCR assay and data analysis; Y-FY designed the study and revised the manuscript. All authors read and approved the final manuscript.