Abstract MicroRNA-101 has been reported as an important factor in carcinogenesis of several malignant tumors. However, its actual role in prognosis among solid malignancies remains unclear. Accordingly, we performed this meta-analysis aiming to identify prognostic significance of miR-101 in solid tumor. Pooled hazard ratios (HRs) with 95% confidence intervals (CIs) for overall survival (OS) or disease-free survival (DFS)/metastasis-free survival (MFS)/progression-free survival (PFS)/relapse-free survival (RFS)/time-to progression (TTP) were estimated with random effects or fixed effects models on the basis of heterogeneity. Subgroup analysis, sensitive analysis and meta-regression analysis were also conducted to clarify the possible confounding factors and investigate the source of heterogeneity. Publication bias was evaluated by using Begg’s and Egger’s tests. A total of 21 studies containing 3753 cases were selected into our quantitative analysis via electronic database search. A lower expression of miR-101 was significantly associated with worse OS (HR = 0.66, 95%CI [0.52–0.85], P = 0.001) and PFS (HR = 0.70, 95%CI [0.51–0.95], P = 0.023) in patients with solid tumor. The under-expression of miRNA-101 is a credible indicator of poorer prognosis in several of solid malignancies.

Citation: Ma X, Bai J, Xie G, Liu Y, Shuai X, Tao K (2017) Prognostic significance of microRNA-101 in solid tumor: A meta-analysis. PLoS ONE 12(7): e0180173. https://doi.org/10.1371/journal.pone.0180173 Editor: Burton B. Yang, University of Toronto, CANADA Received: May 3, 2017; Accepted: June 12, 2017; Published: July 25, 2017 Copyright: © 2017 Ma et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Data Availability: All relevant data are within the paper and its Supporting Information files. Funding: This work was supported by National Natural Science Foundation of China (http://www.nsfc.gov.cn/publish/portal1/)(No. 81572413) and Scientific and Technological Application Foundation Project of Wuhan (No.2015060101010044). No. 81572413, KT, role in conceptualization, investigation, resources, writing review & editing, visualization, supervision, project administration and funding acquisition; No.2015060101010044, KT, role in conceptualization, investigation, resources, writing review & editing, visualization, supervision, project administration and funding acquisition. Competing interests: The authors have declared that no competing interests exist.

Introduction MicroRNAs (miRNA, miRs) are a subset of small non-coding RNA molecules that are approximately 18–22 nucleotides in length. MiRNAs play crucial regulatory roles in gene expression at the post-transcriptional level [1, 2]. The major mechanism of miRNA action is the interaction with the 3’-UTR of the targeted gene mRNA, followed by degradation of the mRNA or inhibition of mRNA protein translation. In human cancers, numerous studies have shown that the expression of miRNAs is deregulated and these miRNAs act as regulatory molecules in many biological processes, including differentiation, proliferation, and apoptosis of tumor cells [3–5]. Several miRNAs are downregulated in many tumors and appear to function as tumor suppressor genes [6]. Among these downregulated miRNAs, miR-101 is one of the most downregulated miRNAs in human cancers, multiple research studies have been exploring the prognostic function of miR-101 in cancer patients in order to find a reliable biomarker to guide for cancer treatment[7, 8]. MiRNA-101, located on chromosome 1[65058434–65058508], is widely known as a tumor suppressive miRNA that is strongly downregulated in several cancers including neuroblastoma, gastric cancer, prostate cancer, nasopharyngeal carcinoma and hepatocellular carcinoma [9–11]. The aberrant expression of miR-101 not only has diagnostic implications but also can predict cancer patient survival [12]. Although an overwhelming majority of evidence has explored a negative prognostic value of miR-101 under-expression across miscellaneous neoplasms, the prognostic impact of miR-101 in malignancies remains controversial. Failed to draw a similar conclusion, Slattery et al[13] and Lv et al [14] presented a worse survival status of patients under stronger miR-101 expression, suggesting an astonishing positive prognostic significance in circumstance of miR-101 under-expression. Therefore, in the present study through gathering available evidence, we carried out an integrated meta-analysis as well as subgroups analysis to identify the relationship between miR-101 expression level and survival of cancer patients by pooling the hazard ratio (HR) from studies addressing the correlation between miR-101 and OS/PFS of patients with malignancies, aiming to provide more theoretical supports for targeted treatment.

Materials and methods Search strategy We performed a thorough search for available literatures in electronic databases of Pubmed, Embase and Web of science until April 2017, using the following words “(microRNA-101 OR miR-101 OR miR101 OR miRNA-101) AND (tumor OR neoplasm OR cancer OR carcinoma OR malignancy)”. In order to avoid missing the potentially related articles, reference lists were also screened. Two authors independently carried out this procedure and any discrepancy was resolved by mutual discussion. Selection criteria Inclusion criteria were as follows:1) studies exploring any of the solid tumor; 2) studies dealing with miR-101 expression and OS/DFS/PFS/RFS/MFS/TTP; 3) studies that categorized patients into low- and high-expression groups based on the miR-101 expression; 4) studies providing HR directly or key information to calculate HR indirectly, such as Kaplan-Meier curves and original survival data; 5) studies assessing miR-101 expression in tissue or blood. The following were the exclusion criteria: 1) studies on myelomas, lymphomas, or leukemia; 2) duplicated or overlapped studies; non-original articles, such as reviews, articles or letters; 3) laboratory studies on cell lines or animals level; 4) studies on a set of microRNAs but not miR-101 alone; 5) studies with a sample-size less than 20 participants. Qualitative assessment Newcastle-Ottawa Scale(NOS)[15] was adopted to evaluate the quality of each eligible article. The scale was revised with certain adaptive modifications to match the practical needs of the pooled analysis. There are three aspects contained in the scale: selection, comparability, and outcome. Stars awarded for each quality item serve as a quick visual assessment. Stars are awarded such that the highest quality studies are awarded up to nine stars. Studies with more than 6 stars were considered as of high quality. Otherwise, studies were excluded from the final meta-analysis. Data extraction All eligible publications were reviewed by Shuai and Bai, The following details of each article were recorded: first author’s name, publication year, cancer type, treatment, sample size, stage of disease, miR-101 test method, the cutoff value to discriminate high or low expression of miR-101, sample sources, follow-up time, extracting method of HR, outcome, NOS and et al. The HR value was extracted directly if it was calculated by a multivariate analysis. Otherwise, the results from univariate analysis were also allowed in the meta-analysis. If both multivariate analysis and univariate analysis were not available, Kaplan-Meier curves were used to extract HR value by using the described method[16]. Statistical analysis The heterogeneity of the studies included in this meta-analysis was assessed by the Q statistic test and the I2 statistic test, where I2 more than 50% indicated evidence of heterogeneity[17]. The random-effects model was selected when I2 was significant (>50%); otherwise, the fixed-effects model was selected. Publication bias was examined using Begg’s funnel plot and Egger’s linear regression test. P<0.05 was considered significant[18]. All analyses were performed with STATA version 12.0 software (Stata Corporation, College Station, TX).

Discussion Numerous profiling studies have demonstrated that miRNA expression levels failed to agree in various types of cancers, miRNAs can be potential biomarkers for cancer prognosis. Increasing data favor the potential use of miR-101 as a cancer prognostic predictor. Recently, genome-wide miRNA expression profiling studies revealed that miR-101 is widely present in various tissues and organs, and its aberrant expression was reported in various cancers including HCC[14, 19, 20, 25], CRC[13, 30], breast cancer[23, 26], NSCLC[22], gliomas[28], and et al. It is indubitable that miR-101 is an important cancer-related miRNA. More and more evidence has demonstrated that miR-101 is frequently downregulated in multiple types of cancer and acts as a tumor suppressor by repressing many critical oncogenes. In hepatocellular carcinoma, Wang et al found c-Myc collaborates with EZH2-containing PRC2 complex in silencing miRNA-101 during hepatocarcinogenesis and lower expression of miR-101 is positively correlated with poorer prognosis [38]. In CRC, it is reported that loss of miR-101 expression promotes Wnt/β-catenin signaling pathway activation and malignancy in colon cancer cells[39]. Similarly, in glioblastoma, Liu et al demonstrated that miRNA-101 inhibits proliferation, migration and invasion of glioblastoma by targeting SOX9 [40], Michiel et al also found that miRNA-101 is down-regulated in glioblastoma resulting in EZH2-induced proliferation, migration, and angiogenesis [41]. Moreover, miRNA-101 reverses temozolomide resistance by inhibition of GSK3β in glioblastoma[42]. In nasopharyngeal carcinoma, MicroRNA-101 inhibits invasion and angiogenesis through targeting ITGA3 and its systemic delivery inhibits lung metastasis[43]. On the other hand, 175 targeted genes validated by experiment and 5206 targeted genes predicted by miRanda software can be found in GCBI website based on classical miRNA-3’-UTR pathway (https://www.gcbi.com.cn/gclib/html/dictSearchAct/MI0000103/miRNA), indicating miRNA-101 may play a complicated role in many gene ontology functions and pathways networks. Further experiments need to be conducted to elucidate the role of miR-101 in carcinogenesis. However, among all the studies referring to the relationship between miRNA-101 and OS/PFS, there were still some contradictory views requiring adequate attention, a comprehensive study is therefore in urgent. To the best of our knowledge, few studies have systemically explored the possible prognostic role of miR-101 down-regulation in solid malignancies before. In order to get a more convincing outcome, HRs for both OS and PFS were calculated independently. On the whole, our quantitative results strongly supported the current mainstream viewpoint that an undesirable impact of miR-101 low expression was related with poor OS and PFS, taking no account of confounding factors. Among all the included studies, on the contrary, two studies [13, 14] highlighted that obvious advantage on survival duration was obtained in miR-101 under-expression cases and no significant HR was found in other three studies [26–28, 35], whose HR value and 95%CI were extracted from K-M survival curve, indicating that this indirect method may impose slightly bias on the HR we calculated. It’s worth mentioning that 4 studies providing K-M curve was excluded due to the lack of clear categorization [9, 44–46]. The appropriate HRs can't be obtained until the K-M survival curve and the exact number of each group are available simultaneously. Based on the outcome of subgroup analysis, significant HR was not found in all subgroup even significant HR was found in each studies, which may partly attribute to a relative small number of studies and high heterogeneity and need further elucidate. Apart from the inspiring outcomes, limitations still exist in this quantitative meta-analysis. First of all, despite the usage of random-effects model and subgroup analysis, the heterogeneity across studies failed to be eliminated completely, which could result in bias of the outcome in certain extent. Secondly, due to lack of direct HR and 95%CI data, we merely extracted the data by using Engauge software indirectly, which may bring about slight error in HR and 95%CI. Thirdly, lack of abundant miR-101 expression data in the global population makes it difficult to set a standard cut-off value for measurement of miR-101 expression levels, categorization between studies did not get in consensus. Additionally. Other parameters that may partially contribute to the heterogeneity were not explored, such as pathological grade and body mass index. In spite of the limitations mentioned above, there are still numerous valuable implications in this comprehensive meta-analysis, which reveals that low expression of miR-101 is associated with unfavorable survival outcomes in patients with various types of carcinomas, particularly with regard to OS. Further large-scale, well-designed and multi-center prospective studies should be conducted to confirm these findings before the application of miR-101 for the prognosis of cancers.

Acknowledgments This study was generously supported by the National Natural Science Foundation of China (No. 81572413) and Scientific and Technological Application Foundation Project of Wuhan (No. 2015060101010044)