Abstract

Objective. Recently, many studies have shown that microRNAs (miRNA) exhibit altered expression in various cancers and may play an important role as prognostic biomarker of cancers. We performed a meta-analysis to evaluate the impact of miR-375 expression in solid tumors on patients’ overall survival (OS). Methods. Studies were identified by searching PubMed, Embace, and Cochrane Library (last search update was in May 2014) and were assessed by further quality evaluation. The pooled hazard ratios (HRs) with 95% confidence intervals (CIs) for total and stratified analyses were calculated to investigate the association between miR-375 expression and cancer patients OS. Results. Our analysis results indicated that downregulation of miR-375 predicted poor OS (HR = 1.91, 95% CI 1.48–2.45, ). Subgroup analyses showed that lower expression of miR-375 was significantly related with poor OS in patients with esophageal carcinoma (HR = 2.24, 95% CI 1.69–2.96, ) and non-small-cell lung cancer (NSCLC) (HR = 1.71, 95% CI 1.31–2.24, ). Conclusions. The findings from this meta-analysis suggest that miR-375 expression is associated with OS of patients with malignant tumors and could be a useful clinical prognostic biomarker.

1. Introduction

Cancer is one of the most common causes of death worldwide and has become a major public health issue [1]. Hence, newer cancer biomarkers with high sensitivity and specificity are essential for the proper detection, treatment, and prognosis of this fatal disease. MicroRNAs (miRNAs), small noncoding RNAs with a length of approximate 22 nucleotides, have been demonstrated to have high potential prognostic value in many cancers. MiRNAs biogenesis begins with transcription of primary transcripts (pri-miRNAs) in cell nucleus. Catalyzed by the ribonuclease Drosha and its essential cofactor DGCR8, pri-miRNAs release a 60–80-nucleotide precursor (pre-miRNA) which is exported to cell cytoplasm and cleaved by Dicer to generate a 22-nucleotide duplex. One strand of the duplex is incorporated into the RNA-induced silencing complex (RISC) by bounding with Argonaute, whereas the other is degraded. Base-pairing between the miRNAs and target mRNAs (usually in the 3′ untranslated region) guides RISC to complementary transcripts, leading to translation repression or the target mRNAs degradation [2]. It is estimated that miRNAs regulate 1/3 to 2/3 of human genes [3]. Thus, the dysregulation of the biogenesis and function of miRNAs is often associated with human diseases, especially malignancies [4]. Accumulating evidence has demonstrated that miRNAs act as oncogenes or tumor suppressors by targeting genes involved in cell differentiation, proliferation, survival, apoptosis, and metastasis (reviewed in [5]). The expression of numerous miRNAs is dysregulated in various cancers, which is often associated with diagnosis, staging, progression, prognosis, and response to clinical therapies (reviewed in [6]).

MicroRNA-375 (miR-375) was originally identified from murine pancreatic β-cell line MIN6 as a pancreatic islet-specific miRNA [7]. The miR-375 gene is located in an intergenic region between the beta-A2 crystallin (cryba2) and undescribed coiled-coil domain-containing protein 108 (ccdc108) genes in human chromosome 2q35 region [8], a genomic region conserving the synteny between humans and mice. Moreover, the sequences of pre-miR-375 in both species present a 100% homology, highlighting the high degree of conservation for this specific miRNA [8]. Further study revealed that miR-375 is a multifunctional miRNA participating in pancreatic islet development, glucose homeostasis, mucosal immunity, lung surfactant secretion, and more importantly tumorigenesis (reviewed in [9]). Recently, miR-375 has been found significantly dysregulated in many human cancers [1015]. In addition, miR-375 has been found to have prognostic value in a variety of tumors. Multiple studies reported a significant association between low miR-375 and poor prognosis [12, 1625]. However, insignificant or even opposite results have been found in some other studies [2630]. According to the evidence we have gained so far, it is still not enough to come to a conclusion whether miR-375 could be used as a potent biomarker for prognosis. Therefore, a systematical and comprehensive meta-analysis was carried out to investigate the relationship between miR-375 expression and the survival of patients with cancer. To our knowledge, this meta-analysis is the first to evaluate the prognostic value of miR-375 in cancer patients.

2. Materials and Methods

This meta-analysis was carried out in accordance with the guidelines of meta-analysis of observational studies in epidemiology (MOOSE) [31].

2.1. Search Strategy

Studies were performed by searching PubMed, Embace, and Cochrane Library (updated by May 2014). The search strategy was “microRNA-375 OR miR-375” AND “tumor OR neoplasm OR cancer OR carcinoma.” No language limitation was applied. The reference lists of relative articles were also screened to further identify potential studies. The comprehensive database search was carried out independently by two authors (Y. Shao and Y. Geng). The disagreements were resolved by consensus.

2.2. Inclusion and Exclusion Criteria

Eligible studies included in this meta-analysis met the following criteria: (1) reporting explicit methods for the detection of miR-375 expression in tumor tissue or blood; (2) investigating the association between miR-375 expression and survival outcome, and the end-points being overall survival (OS); (3) reporting sufficient data to estimate the hazard ratio (HR) and 95% confidence intervals (CI) according to miR-375 expression. If a study reporting the same patient cohort was included in several publications, only the most recent or complete study was selected. Studies of case reports, letters, reviews, and animal trails were excluded. The titles and abstracts of the identified articles were evaluated independently by two reviewers and the irrelevant articles were excluded. The full text of the extracted articles was carefully examined for comprehensive evaluation. The disagreements were resolved by a discussion with the third reviewer (W. Gu).

2.3. Data Extraction and Qualitative Assessment

The data from all eligible studies were extracted by two investigators independently, which included first authors’ surname, publication year, origin of population, sample number, tumor type, follow-up period, first-line therapies, source of miRNA, miR-375 assessment methods and the cut-off definition, and HR of miR-375 expression for OS as well as 95% corresponding interval (CI) and value. If a study reported the results by both univariate and multivariate analyses, the latter was selected since it considered the confounding factors and was therefore more precise.

The quality of each study was assessed independently by two researchers according to the Newcastle-Ottawa Quality Assessment Scale (NOS) [32]. For quality, scores ranged from 0 (lowest) to 9 (highest), and studies with scores of 6 or more were rated as high quality.

2.4. Statistical Analysis

The high or low expression of miR-375 was defined according to the cut-off values provided by the authors. HRs and their 95% CIs were combined to evaluate the association between miR-375 expression and prognosis. If the statistical variables were described in the study, we pooled them directly. Otherwise, the statistical variables were calculated from available numerical data in the articles according to the methods described by Tierney et al. [33]. To reduce reading variability, the data from Kaplan-Meier survival curves was evaluated by three independent persons as described by Engauge Digitizer version 4.1. The additional information and original data needed for the meta-analysis were acquired by contacting with the corresponding authors of eligible articles. An observed HR greater than 1 indicated a worse prognosis in patients with miR-375 downregulation. Statistical heterogeneity was assessed by visual inspection of forest plots, by performing the Chi-square test (assessing the value) and calculating the statistic [34, 35]. If the value was less than 0.05 and/or exceeded 50%, indicating the presence of heterogeneity, a random-effects model (the DerSimonian-Laird method) was used. Otherwise, the fixed-effects model (the Mante-Haenszel method) was used. Subgroup analysis was further performed to explore the source of heterogeneity. Publication bias was estimated by visually assessing the asymmetry of an inverted funnel plot. Furthermore, Begg’s test and Egger’s test were performed to provide quantitative evidence of publication bias. If a publication bias was observed, it was adjusted by the use of the Duval and Tweedie trim-and-fill method [36]. All analyses were performed using STATA vision 12.0 (Stata Corporation, College Station, TX, USA). A value less than 0.05 was considered to be statistically significant except where otherwise specified.

3. Results

3.1. Study Characteristics

According to the criteria mentioned in Materials and Methods, 504 abstracts were initially selected. However, 488 irrelevant abstracts were excluded. Sixteen full-text articles were reviewed for further evaluation. Of them, 3 were excluded because the data of HRs or OS were not available [2830]. The remaining 13 articles contained 16 studies [12, 1627], because one article included two independent cohort studies [16] and another article included three studies [12]. Thus, 16 studies were included in this meta-analysis, which were published between 2009 and 2014 (Figure 1). The total number of patients in all studies was 1,652, ranging from 37 to 249 patients. The category of cancers included esophageal carcinoma (8 studies), non-small-cell lung cancer (NSCLC, 3 studies), glioma, breast cancer, gastric cancer, head and neck squamous cell carcinoma (HNSCC), and pancreatic ductal adenocarcinoma (PDAC). Quantitative RT-PCR was used to detect miRNAs expression in all studies except one. The expression of miR-375 was detected in tumor tissues (11 studies) or blood samples (5 studies). The cut-off values of miR-375 varied in different studies. HRs were estimated in 8 studies and reported in the text of other studies. The major characteristics of the 16 eligible studies are listed in Table 1.

3.2. Qualitative Assessment

Sixteen eligible studies included in our meta-analysis were assessed for quality according to the NOS. The quality of all included studies varied from 6 to 9, with a mean of 6.6. A higher value indicated a better methodology. Therefore, all studies were included in the subsequent analysis.

3.3. Meta-Analysis Results

As the studies evaluating OS were of obvious statistical heterogeneity (), we used a random-effects model to pool the HRs. The result showed that downregulated miR-375 was significantly associated with poor OS outcome in various carcinomas, with the pooled HR of 1.91 (95% CI 1.48–2.45, ) (Table 2; Figure 2).

Considering the heterogeneity among these studies, the effect of miR-375 expression was further evaluated by subgroup analysis. The subgroups were classified according to the main characteristics such as tumor type, source of miRNA, miRNA assay method, type of method used to obtain the HR, patient origin, and analysis type. In the subgroup of tumor type, we found the downregulation of miR-375 was significantly associated with worse OS in esophageal carcinoma (HR = 2.24, 95% CI 1.69–2.96; ; fixed-effects model) and NSCLC (HR = 1.71, 95% CI 1.31–2.24; fixed-effects model), without any heterogeneity in the data (, ; ,, resp.) (Table 2 and Figure 3). There was only one study that evaluated the association between lower miR-375 expression and OS in HNSCC, gastric cancer, breast cancer, PDAC, and glioma, respectively, and therefore, these tumors were defined as “other cancers.” Combined data from these five studies showed that decreased miR-375 expression was not correlated with poor OS (HR = 1.59, 95% CI 0.71–3.58; random-effects model) and with significant statistical heterogeneity () (Table 2). The association between lower miR-375 expression and worse OS outcome was statistically significant in other subgroups, including HR reported in text (HR = 1.97, 95% CI 1.57–2.47, fixed-effects model; for heterogeneity test, = 10.8%), miR-375 assay by qRT-PCR (HR = 1.88, 95% CI 1.42–2.48, random-effects model; for heterogeneity test, = 46.2%), tissue-source of miRNA (HR = 2.16, 95% CI 1.73–2.71, fixed-effects model; for heterogeneity test, ), multivariate analysis (HR = 1.97, 95% CI 1.57–2.47, fixed-effects model; for heterogeneity test, ), Chinese patients (HR = 1.96, 95% CI 1.62–2.37, fixed-effects model; for heterogeneity test, ), and American patients (HR = 1.97, 95% CI 1.30–2.99, fixed-effects model; for heterogeneity test, ) (Table 2). The expression of miR-375 did not show prognostic impact in subgroups of HRs by data extrapolated, univariate analysis, and blood-source of miRNA. The subgroup of miR-375 detected by MISH only included one study; therefore the result related entirely to the individual study.

3.4. Heterogeneity Analysis

Sensitivity analysis was performed by sequential omission of individual studies using the fixed-effects model, demonstrating that the study by Madhavan et al. [20] apparently influenced the overall results (Figure 4). When the Galbraith plot was analyzed, one study was identified as outliers of heterogeneity [20]. By excluding this study from the analysis, similar pooled HR and significance were obtained (HR = 1.96, 95% CI 1.65–2.33, ) but heterogeneity was absent (, ). We also conducted a meta-regression to explore the potential factors responsible for the heterogeneity. As a result, all these factors including tumor type, miR-375 assay method, follow-up time, publication year, patients origin, and cut-off values did not contribute to the heterogeneity obviously.

3.5. Publication Bias

Begg’s funnel plot and Egger’s test were used to evaluate the publication bias (Figure 5). The values of Egger’s and Begg’s tests were all over 0.05 ( for Begg’s test; for Egger’s test). Hence, there was no evidence for significant publication bias in the meta-analysis.

4. Discussion

Recently, genome-wide miRNA expression profiling studies revealed that miR-375 is widely present in various tissues and organs, and its expression is significantly aberrant in malignant tumors, such as HNSCC, NSCLC, melanoma, glioma, hepatocellular, esophageal, gastric, breast, and prostate cancer [1016, 19, 27]. It is indubitable that miR-375 is an important cancer-related miRNA. Increasing evidence has demonstrated that miR-375 is frequently downregulated in multiple types of cancer and acts as a tumor suppressor by repressing many critical oncogenes (Table 3). In hepatocellular carcinoma, the restoration of miR-375 in cancer cells decreased cell proliferation, clonogenicity, migration, and invasion and induced G1 arrest and apoptosis [10]. Moreover, MTDH was directly regulated by miR-375 in both hepatocellular carcinoma and HNSCC [10, 37, 38]. In gastric cancer, miR-375 was frequently downregulated and inhibited gastric cancer cell proliferation via targeting Janus kinase 2 [11]. Tsukamoto et al. also found a tumor suppressive role of miR-375 in gastric cancer. The ectopic expression of miR-375 reduced cell viability and induced apoptosis by targeting PDK1 and YWHAZ [39]. Similarly, in esophageal squamous cell carcinoma miR-375 was confirmed to inhibit cell proliferation, colony formation, and metastasis in vitro and in vivo [23]. Another study also verified the tumor-suppressive effect of miR-375 in esophageal cancer cell lines by targeting PDK1 [40]. However, there were still some contradictory views requiring adequate attention. It was reported that miR-375 was significantly upregulated in tumor tissues or serum of prostate carcinoma patients [14, 41, 42], but the role of miR-375 in prostate cancer was unclear. miR-375 was also upregulated in estrogen receptor alpha- (ERα-) positive breast cancer cell lines, which promoted cell proliferation and induced ERα upregulation via RASD1, a negative regulator of ERα [15]. In addition, higher expression of miR-375 was reported to result in progression of invasive lobular breast carcinoma [43]. Recently, miR-375 was found to participate in the process of epithelial-to-mesenchymal transition (EMT). It could partly reverse EMT-like properties in MCF-7 cells [44], which suggested its promoting role in tumor metastasis. Therefore, whether miR-375 acts as a tumor suppressor or an oncogene is unclear. Although plenty of studies focused on the potential application of miR-375 as a prognostic biomarker, the prognostic value of miR-375 expression in cancer patients is still controversial.

In the present study, we carried out a meta-analysis to evaluate the prognostic value of miR-375 by combining 16 studies with 1,652 patients. Our results demonstrated that lower miR-375 expression was associated with poor survival in patients with various carcinomas (HR = 1.91 95% CI 1.48–2.54, , random-effects model), indicating that miR-375 may serve as a positive prognostic marker for solid tumors. In the subgroup analyses, the association between lower miR-375 expression and worse OS was statistically significant in most subgroups, especially in esophageal carcinoma (HR = 2.24, 95% CI 1.69–2.96; ; fixed-effects model) and NSCLC (HR = 1.71, 95% CI 1.31–2.24; fixed-effects model). Recently, there were only two studies exploring the relationship between miR-375 expression and the survival of cancer patients. Consistent with our results, they also indicated low miR-375 expression was significantly associated with poor OS of patients with esophageal carcinoma. However, the results of these two studies are not credible enough because they did not eliminate interference factors from their meta-analyses, including miR-375 expression in pericarcinoma tissues of the same patient cohort, which may possibly come to a conclusion with deviation.

In sensitivity analysis, the study of Madhavan et al. was identified as an outlier of heterogeneity. It indicated a positive correlation between upregulated miR-375 expression and poor OS in metastatic breast cancer (MBC) patients (HR = 0.31, 95% CI 0.11–0.88) [20], which may be attributed to mesenchymal-to-epithelial transition (MET), a process essential for successful colonization and establishment of metastasis induced by high miR-375 expression [44]. The significant prognostic effect of high miR-375 expression in this study with obvious heterogeneity possibly led to the negative statistical results in some subgroup analyses, including data extrapolated, univariate analysis, and blood-source of miRNA: in these subgroups, there was no association between lower expression of miR-375 and poor prognosis.

Although the predictive value of miR-375 was statistically proved by the meta-analysis in this study, it should be carefully comprehend for the following reasons. First, four studies which explored the association between miR-375 expression and prognosis were not included in the meta-analysis because they did not provide available data of HRs or OS, and three of them showed insignificant or even opposite results. These studies might influence the reliability of our results. Secondly, due to the lack of uniform cut-off value in miR-375 expression, the cut-off values applied by different researchers deviated from the actual value more or less, which may affect the validity of miR-375 as a predictive factor in cancer prognosis. Thirdly, several HRs were calculated based on the data extracted from the survival curve, bringing minor deviations. Fourthly, the statistical heterogeneity was not obvious in this meta-analysis, but the clinical and methodological heterogeneity existed in the baseline demographic characteristics, including population, tumor types and stages, the cut-off value of miR-375 expression, and duration of follow-up. Finally, although no significant publication bias was detected in this meta-analysis, the results still need to be verified by a large number of publications.

In conclusion, this meta-analysis summarized the global researches on the relationship of aberrant miR-375 expression and the prognosis of patients with cancer and clarified that downregulation of miR-375 is significantly associated with poor survival in patients with various types of carcinoma, especially in esophageal carcinoma and NSCLC. In view of the limitation of the current analysis, it should be cautious to appreciate the conclusion, and further prospective multicenter studies designed adequately with larger sample size are needed to verify the association between miR-375 and cancer prognosis as well as the efficiency of therapies.

Conflict of Interests

The authors declare that they have no conflict of interests.

Authors’ Contribution

Yingjie Shao and Yiting Geng contributed equally to this work.

Funding

This research project was supported by the National Natural Science Foundation of China (NSFC) (81171653) and the General Program of Jiangsu Province Department of Health (H201350).

Acknowledgments

The authors would like to acknowledge all the investigators, site staff, and their families involved in the study. They also thank Professor Jingting Jiang, M.D., for his editorial support.