Skip to main content

Comprehensive and integrative analysis identifies microRNA-106 as a novel non-invasive biomarker for detection of gastric cancer

Abstract

Background

Recently, accumulating evidences have revealed that microRNA-106 (miR-106) may serve as a non-invasive and cost-effective biomarker in gastric cancer (GC) detection. However, inconsistent results have prevented its application to clinical practice.

Methods

As a result of this, a comprehensive meta-analysis was conducted to evaluate the diagnostic performance of miR-106 alone and miR-106-related combination markers for GC detection. Meanwhile, an integrative bioinformatics analysis was performed to explore the function of miR-106 at the systems biology level.

Results

The results in our work showed that sensitivity of 0.71 (95% CI 0.65–0.76) and specificity of 0.82 (0.72–0.88), with the under area AUC (area under the curve) value of 0.80 (0.76–0.83) for miR-106 alone. Prospectively, miR-106-related combination markers improved the combined sensitivity, specificity and AUC, describing the discriminatory ability of 0.78 (0.65–0.87), 0.83 (0.77–0.89) and 0.88 (0.85–0.90) in the present analysis. Furthermore, targets of miR-106 were obtained and enriched by gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis, revealing their associations with the occurrence and development of GC. Hub genes and significant modules were identified from the protein–protein interaction networks constructed by miR-106 targets and found closely associated with the initiation and progression of GC again.

Conclusions

Our comprehensive and integrative analysis revealed that miR-106 may be suitable as a diagnostic biomarker for GC while microRNA combination biomarkers may provide a new alternative for clinical application. However, it is necessary to conduct large-scale population-based studies and biological experiments to further investigate the diagnostic value of miR-106.

Background

Gastric cancer (GC), a major public health challenge, is one of the leading causes of cancer death worldwide [1]. There are limited detection methods for early diagnosis and few effective screening procedures in some countries. The most reliable program for diagnosis is mainly based on endoscopy and biopsy [2]. However, this is invasive and inconvenient for patients to undergo. Consequently, most patients can only be diagnosed precisely in advanced stages when the clinical outcomes are poor [3]. Therefore, there is a great need to explore new accurate and efficient, preferentially non-invasive, markers for early detection of GC.

In recent time, accumulating evidences have suggested that microRNAs may serve as novel biomarkers for cancer detection. MicroRNAs are a class of small non-coding RNAs with intermediate posttranscriptional regulation of the target genes [4]. A large number of studies have demonstrated that microRNAs play vital roles in a wide variety of physiological processes including cancer cell growth, differentiation, invasion, and metastasis [5]. Moreover, a number of studies have indicated that circulating microRNAs have high degree of stability and tolerance even under unfavorable physiochemical conditions including extreme variations in pH, temperature and freeze–thaw cycles [6]. It is also promising that microRNAs have outstanding stability in multiple clinical samples including plasma, serum, feces and tissue, which enables them to be detectable effortlessly [7]. Given their critical involvement in the vital biological processes and perfect biomarker features mentioned above, microRNAs could be considered as good candidates for using as non-invasive markers, and the application of them as biomarkers for early detecting GC is viable [8].

As one of the most representative microRNA biomarkers, microRNA-106 (miR-106) has been extensively studied by a great number of researches in several cancers. MiR-106 belongs to the miR-17 family, one of the most common studied onco-microRNA groups, which includes miRs-17, -20a, -20b, -93, -106a and -106b. MiR-106a is a member of the miR-106a-92 cluster located on chromosome Xq26.2 while miR-106b is located at 7q21 [9, 10]. There have been several studies indicating that both of investigated miR-106 could be expressed in the same individuals of gastric tumour tissues. Several studies have previously reported that circulating miR-106 could specifically serve as a pivotal and promising biomarker for GC [11]. Nevertheless, the suitability of circulating miR-106 in early detection and diagnosis of GC remains inconsistent due to different sample sizes, disease statuses, sample sources, detection methods, and other uncontrolled factors. Moreover, the potential molecular mechanism of miR-106 is still poorly understood for the present insufficient knowledge.

In the present study, a comprehensive meta-analysis was performed to obtain a better understanding of the clinical feasibility of miR-106 as promising biomarker for early detection and diagnosis of GC. By focusing not only on a single miR-106 marker, we explored whether combination biomarkers based on miR-106 are more effective than individual miR-106. Furthermore, an integrative bioinformatics analysis was carried out to evaluate the functions of miR-106 at the systems biology level.

Methods

Literature search strategy

A comprehensive computerized literature search for articles (up to December 27, 2017) was carried out based on several electronic databases including PubMed, EMbase, Web of Science and the Cochrane Library using the following search terms: (“cancer” OR “tumor’’ OR “carcinoma” OR “neoplasm”) AND (“gastric” OR “stomach’’ OR “gastrointestinal’’ OR “digestive” OR “GC”) AND (“microRNA-106” OR “miR-106” OR “miR-106a” OR “miR-106b”). In addition, the references of identified articles were examined for all relevant studies.

Eligibility criteria

The studies qualified to be included should meet the following criteria: (1) they investigated the potential of circulating (blood, serum, and plasma) or other sources of miR-106 for detecting GC; (2) they used the gold standard to make a definitive diagnosis of GC; (3) they provided adequate data which can be used to calculate the rates of true positive (TP), false positive (FP), false negative (FN), and true negative (TN).

In addition, the studies were excluded if (1) they were obviously not associated with our topic; (2) they published in forms of reviews, letters, case reports, or editorials; (3) they were non-English publications; or (4) they provided unqualified data.

Data extraction

Data were collected independently by two investigators (Peng and Shen) from the articles based on standardized forms. The following characteristics from each study were extracted: (1) first author; (2) publication year; (3) research country; (4) study population; (5) patient characteristics (age, gender, cancer type, etc.); (6) sample sources; (7) participants numbers; (8) miR-106 measuring methods, (9) diagnostic data including sensitivity, specificity, true positive (TP), false positive (FP), false negative (FN) and true negative (TN).

Quality assessment

The quality of each study enrolled in our analysis was assessed independently by two investigators on the basis of the QUADAS-2 (Quality Assessment of Diagnostic Accuracy Studies 2) [12]. The QUADAS-2 assessment tool contains 4 domains including patient selection, index test, reference standard, and flow and timing supported by signaling questions (yes, no, or unclear), risk of bias (high, low, or unclear) and concerns about applicability (high, low, or unclear). An answer of “yes” gets a score of 1, which indicates a low risk of bias, while an answer of “no” or “unclear” gets a score of 0 which means a potential high risk of bias.

Statistical analysis for meta-analysis

The overall diagnostic results were estimated using the TP, FP, FN, and TN test results extracted directly from each study or via recalculation based on sensitivity and specificity together with other data retrieved from each study. The bivariate meta-analysis model was applied to calculate the pooled sensitivity, specificity, positive likelihood ratios (PLR), negative likelihood ratios (NLR), and diagnostic odds ratio (DOR) along with their corresponding 95% confidence intervals (CIs) [13]. The summary receiver operator characteristic (SROC) curve was generated based on the sensitivity and specificity of each study [14]. In addition, we calculated the corresponding area under the SROC curve (AUC) for the quantitative assessment of diagnostic power. Furthermore, the heterogeneity across studies was examined by using Q test and I2 test [15]. A P value ≤ 0.1 from Q test and an I2 ≥ 50% from I2 test suggest the presence of significant heterogeneity among eligible study. The heterogeneity caused by threshold effect was quantified using Spearman correlation analysis. Possible sources of heterogeneity in the aspect of non-threshold effect were explored by carrying out subgroup, meta-regression, and sensitivity analyses [16]. We examined the potential publication bias through Deeks’ funnel analysis [17]. Statistical analyses were performed in STATA (version 14.0) and Meta-DiSc statistical software (version 1.4) software. Values of P < 0.05 were considered to represent statistical significance.

Integrative functional analysis of miR-106

To further explore the function of miR-106, an integrative functional analysis was performed. Targets of miR-106 were obtained from miRTarBase, a resource for experimentally validated microRNA-target interactions [18]. Updated in 2018 with plenty of integrative improvements, miRTarBase has been a more comprehensively annotated microRNA-target interactions database approved by experimental evidence compared with its previous versions and other databases in the field of microRNA-target prediction research. We conducted GO and KEGG pathway enrichment analysis by separately mapping the target mRNAs of miR-106a and miR-106b to the online tool DAVID (Database for Annotation, Visualization, and Integrated Discovery) [19,20,21]. P-value < 0.05 and count ≥ 2 were considered as the cut-off criteria.

PPI network analysis

The target mRNAs regulated by miR-106a and miR-106b were mapped to the Search Tool for the Retrieval of Interacting Genes (STRING) database to retrieve the PPI information, respectively [22]. In the present study, we only selected the PPI data with the combined score > 0.4. Network visualization was used the powerful tool Cytoscape [23]. Meanwhile, we identified several hub genes of miR-106a and miR-106b targets by using three different methods including betweenness centrality, closeness centrality and degree centrality based on the plug-in CytoNCA [24]. In addition, the significant modules in the PPI network were screened by plug-in Molecular Complex Detection (MCODE) of Cytoscape. As for a further comment, KEGG pathway analysis was performed with the hub genes and genes involved in the selected modules by DAVID.

Results

Study selection

After initial searches of the databases mentioned above, a total of 369 articles were retrieved preliminarily, and the flow diagram of the literature search was shown in Fig. 1. After careful exclusion of inappropriate ones according to the inclusion criteria step by step, ten published articles were finally included in our study, including twelve studies for miR-106 alone and five for miR-106-related combination markers [25,26,27,28,29,30,31,32,33,34].

Fig. 1
figure 1

Flow diagram of the study selection process

Baseline characteristics of included studies

Among the included publications for evaluating miR-106 alone in detecting GC, a total of 10 articles on 12 studies were analyzed involving 788 cases and 701 healthy people as the control group. There were five studies investigating serum miR-106, five studies assessing plasma miR-106, one study evaluating tissue miR-106 and one study involving miR-106 in gastric juice. Among all these studies, six focused on miR-106a and while 6 studies concentrated on miR-106b. None of the studies in Table 1 evaluated both miRNAs (miR-106a and miR-106b). With respect to value of miR-106-related combination markers, 4 publications on 5 studies involved 322 patients and 263 controls. Sample sources were all serum (n = 5). All the studies measured the expression level of miR-106 and the related combination markers by using the quantitative real-time polymerase chain reaction (qRT-PCR) assay. The detailed information was shown in Tables 1 and 2. The quality assessments of the included studies indicated that overall they were of moderate to high quality, which enhanced the reliability of our analysis.

Table 1 The main features of the included studies on individual miR-106
Table 2 The main features of the included studies on miR-106-related combination markers

Diagnostic value of miR-106 in GC

As indicated in the forest plot (Fig. 2a, b), the existence of significant heterogeneity among individual studies was observed as the Q value was 31.46 (P < 0.001) and I2 value was 65.04% (95% CI 43.54–86.53) for sensitivity, and the Q value was 94.70 (P < 0.001) and I2 value was 88.38% (95% CI 83.04–93.72) for specificity. Thus, a random-effects model was applied to evaluate the pooled results. Overall, the pooled assessment outcomes were as follows: sensitivity, 0.71 (95% CI 0.65–0.76); specificity, 0.82 (0.72–0.88); PLR, 3.86 (2.48–6.01); NLR, 0.36 (0.28–0.45); and DOR, 10.84 (5.88–20.00), respectively. The pooled PLR meant that patients with GC had a nearly four-fold greater chance of being miR-106 positive than that for patients without GC while the combined NLR of 0.36 indicated that expected proportion of patients having GC is 36% if the miR-106 is negative. The DOR value suggested that someone who was screened to be positive for GC with a high expression of miR-106 had a 10.84-fold higher possibility of actually suffering from GC than someone with a negative GC result. Finally, the SROC curve (Fig. 3a) was plotted and the corresponding AUC was calculated of 0.80 (0.76–0.83), indicating moderate diagnostic accuracy overall.

Fig. 2
figure 2

Forest plots of sensitivities and specificities from test accuracy studies in the diagnosis of GC. a, b Forest plots of sensitivities and specificities for miR-106 alone; c, d forest plots of sensitivities and specificities for miR-106-related combination markers

Fig. 3
figure 3

The SROC curves in the diagnosis of GC. a SROC curve of overall including the outliers for miR-106 alone; b SROC curve of outliers excluded for miR-106 alone; c SROC curve for miR-106 alone in serum samples; d SROC curve for miR-106-related combination markers in serum samples. SROC summary receiver operator characteristic, GC gastric cancer

The threshold effect may result from the differences between sensitivity and specificity. The representative way of assessing the threshold effect is estimating Spearman’s correlation coefficient of sensitivity and specificity. According to our results, Spearman’s correlation coefficient was calculated to be − 0.049, with a P value of 0.88 (P > 0.05), suggesting that no obvious heterogeneity generating from the threshold effect.

In the present study, subgroup analysis was applied in order to identify the potential sources of heterogeneity (see Table 3). The results implied that plasma miR-106 had relatively high diagnostic accuracy for GC detection compared with serum miR-106, with sensitivity of 0.77 (95% CI 0.69–0.83) vs. 0.64 (0.56–0.72), specificity of 0.82 (0.67–0.91) vs. 0.82 (0.66–0.91), and AUC of 0.83 (0.80–0.86) vs. 0.73 (0.69–0.77). Among the twelve studies, ten studies detected the miR-106 in circulating samples. Hence, subgroup analysis was also performed by circulating samples. The pooled sensitivity, specificity and AUC for circulating miR-106 were 0.71 (0.64–0.78), 0.82 (0.72–0.89), and 0.81 (0.77–0.84). We found no significant difference in the diagnostic accuracy between studies with miR-106a and miR-106b, with AUC of 0.78 (95% CI 0.74–0.82) vs. 0.81 (0.77–0.84) although miR-106b exhibited higher diagnostic sensitivity of 0.74 (0.65–0.81) compared with studies with miR-106a, for which the value was 0.68 (0.60–0.75) while miR-106a was associated with specificity of 0.85 (0.74–0.92) compared with the value for miR-106b of 0.78 (0.63–0.88). Subgroup analysis by the sample size suggested that a large sample size had higher specificity and AUC of 0.85 (0.75–0.91) and 0.82 (0.78–0.85) compared to a small sample size with specificity of 0.73 (0.56–0.85) and AUC of 0.75 (0.71–0.78), indicating large sample size had relatively overall high diagnostic accuracy than small sample size, although a small sample size exhibited higher sensitivity of 0.74 (0.66–0.80) compared with a large sample size of 0.70 (0.62–0.77).

Table 3 Pooled results of diagnostic accuracy of miR-106 and combination biomarkers in gastric cancer

Meta-regression analysis was conducted to reveal the potential sources of the heterogeneity. We considered 5 covariates (publication year, country, sample size, sample source and miR-106 classification) may contribute to the heterogeneity. The results showed that neither publication year, nor sample size or sample source or miR-106 classification was the source of heterogeneity, but the country has influence in specificity (P < 0.05).

The goodness of fit and bivariate normality analyses suggested that the bivariate random-effects model was robust for the meta-analysis (Fig. 4). Besides that, one deviated study that may affect the robustness of the meta-analysis was identified with the method of influence analysis and outlier detection. After excluding the deviated study, no significant changes in sensitivity (0.71 vs. 0.69), specificity (0.82 vs. 0.80), PLR (3.86 vs. 3.49), NLR (0.36 vs. 0.38), DOR (10.84 vs. 9.07), and AUC (0.80 vs. 0.77) were observed between the overall analysis with and without outliers, which indicated that there was high robustness in our meta-analysis.

Fig. 4
figure 4

Sensitivity analysis results. a Goodness of fit; b bivariate normality; c influence analysis; d outlier detection

Diagnostic value of miR-106-related combination marker in GC

As shown in Fig. 2c, d, I2 values for sensitivity and specificity were 76.97% (95% CI 56.57–97.38%, P < 0.001) and 49.95% (0.00–100.00%, P < 0.001), respectively, indicating significant heterogeneity for sensitivity and moderate heterogeneity for specificity. The pooled sensitivity, specificity, PLR, NLR, and DOR are 0.78 (0.65–0.87), 0.83 (0.77–0.89), 4.69 (2.98–7.37), 0.27 (0.15–0.46), and 17.68 (6.78–46.09), respectively. The SROC curve is shown in Fig. 3d with AUC = 0.88 (0.85–0.90), indicating a relatively higher accuracy in GC detection. Since all the combination markers were detected in serum, we compared them with serum miR-106 alone (Fig. 3c). Serum miR-106-related combination marker had a higher level of predictive power than serum miR-106 alone, with sensitivity of 0.78 (0.65–0.87) vs. 0.64 (0.56–0.72), specificity of 0.83 (0.77–0.89) vs. 0.82 (0.66–0.91), and AUC of 0.88 (0.85–0.90) vs. 0.73 (0.69–0.77).

In the present study for combination biomarkers, Spearman’s rank correlation was also evaluated to explore the potential heterogeneity from threshold effect. It was revealed from the results that heterogeneity may generate from the threshold effect, from a Spearman’s correlation coefficient of − 0.90 with P = 0.037.

Publication bias

Begg’s funnel plot and Egger’s test were applied to evaluate the presence of publication bias (Fig. 5). The funnel plots indicated no symmetry for all enrolled studies and Deeks’ test returned P values of 0.56 and 0.45 for miR-106 and miR-106-related combination markers, respectively, revealing no obvious publication bias in the present study. However, it is difficult for judging publication bias whether or not exists due to the limited number of studies enrolled in the analysis for miR-106-related combination markers.

Fig. 5
figure 5

Deeks’ funnel plots for the assessment of potential bias in the meta-analysis for diagnosis. a Funnel plot of the studies on miR-106 alone; b funnel plot of the studies on miR-106-related combination markers

Integrative functional analysis results of miR-106

Based on the above results, we wondered why miR-106 could serve as a promising biomarker for detecting GC. We supposed that if the identified miR-106 could be a biomarker of GC, the genes regulated by miR-106 should also be involved in the initiation and progression of GC. Therefore, we performed GO and pathway enrichment analyses on targets of miR-106a and miR-106b to explore the function and pathogenesis of them, respectively.

We performed the GO enrichment analysis by separately mapping the target genes of miR-106a and miR-106b to the online software DAVID at three different levels: molecular function (MF), cell component (CC) and biological processes (BP). The top 10 items of each GO level that were significantly enriched by the target genes were illustrated at Fig. 6. The enriched GO terms in BP for miR-106a mainly included the regulation of transcription, apoptotic process while the enriched GO terms in BP for miR-106b were also associated with regulation of transcription and apoptotic process. For the CC items, the target genes of miR-106a were enriched in the hallmarks of a cell such as nucleoplasm, cytosol, cytoplasm, nucleus and the target genes of miR-106b were also enriched in the same components including nucleoplasm, cytosol, cytoplasm and nucleus. Most GO MF items for miR-106a converged on the binding functions such as protein binding, DNA binding, transcription factor binding while most GO MF items for miR-106b were also related to protein binding, DNA binding and transcription factor binding.

Fig. 6
figure 6

GO annotation of miR-106 target genes. a Top 10 GO items for target genes of miR-106a; b top 10 GO items for target genes of miR-106b. GO gene ontology, BP biological processes, CC cell component, MF molecular function

The top 20 significantly enriched pathways of miR-106a and miR-106b were outlined in Fig. 7a, b, respectively. To our surprise, the enriched KEGG terms of miR-106a and miR-106b were approximately the same. What’s more, we identified several pathways from the top 20 enriched KEGG terms, namely pathways in cancer, p53 signaling pathway, cell cycle, TGF-beta signaling pathway and Proteoglycans in cancer, which were related to the occurrence and development of GC.

Fig. 7
figure 7

Pathway enrichment results for miR-106 target genes. a Top 20 pathways enriched by target genes of miR-106a; b top 20 pathways enriched by target genes of miR-106

The consistent enrichment results of miR-106a and miR-106b including GO and KEGG pathway analysis indicated that although the different location on chromosome, miR-106a and miR-106b may have the similar functions and could be both used in detecting GC. Meanwhile, the GO and KEGG pathway enrichment results explained why miR-106 why miR-106 could serve as a promising biomarker for detecting GC to some extent.

PPI network construction and identification of key target nodes

The information retrieved by STRING was integrated and set up the PPI network of miR-106a and miR-106b targets, respectively. As a result, a PPI network with statistical significance made up of 639 nodes was identified with the set of 896 target genes of miR-106a. Meanwhile, with the set of 1156 target genes of miR-106b, a statistically significant network consisting of 960 nodes were screened. We estimated three network parameters including degree centrality, betweenness centrality, and closeness centrality of the constructed network to explore the relationships between microRNA targets. Each method screened the hub genes in the network. Top 50 genes evaluated by the three methods in the PPI network were screened and intersected. Finally, 30 and 29 genes which could be screened by all the three methods were identified as hub genes for miR-106a and miR-106b, respectively. The regulatory networks were reconstructed with miR-106a and miR-106b and their identified hub nodes, plotted in Fig. 8. Furthermore, we evaluated the biological function of the selected key miR-106 targets, we found that these key target nodes regulated by miR-106a and miR-106b both played a role in FoxO signaling pathway, pathways in cancer, PI3K-Akt signaling pathway, cell cycle, and p53 signaling pathway (Fig. 8).

Fig. 8
figure 8

PPI network construction results. a Degree distributions of nodes for network constructed with miR-106a targets; b degree distributions of nodes for network set up with miR-106b targets; c hub genes of network for miR-106a targets; d hub genes of network for miR-106b targets; e pathway enrichment results for the selected hub genes of miR-106a targets network; f pathway enrichment results for the selected hub genes of miR-106b targets network. PPI protein–protein interaction

Next, active modules were identified (Fig. 9). According to KEGG pathway enrichment analysis, the genes involved in the significant module of miR-106a targets network were mainly associated with p53 signaling pathway, pathways in cancer, FoxO signaling pathway, microRNAs in cancer, PI3K-Akt signaling pathway and cell cycle while the genes involved in the significant module of miR-106b targets network were related to pathways in cancer, p53 signaling pathway, FoxO signaling pathway, microRNAs in cancer and cell cycle.

Fig. 9
figure 9

The significant modules from the PPI network. a The significant module in the PPI network for miR-106a targets; b the significant module in the PPI network for miR-106b targets; c pathways enriched by all the nodes involved in the identified module for miR-106a; d pathways enriched by all the nodes involved in the identified module for miR-106a. PPI protein–protein interaction

Discussion

Over the decades, microRNAs have gained great attention in scientific researches for cancer detection, diagnosis, and treatment as they possess perfect biomarker characteristics and are highly involved in the occurrence and development of a variety of cancers [35]. As one of the most researched microRNAs, miR-106 has been suggested by emerging evidences that it could be a novel potential biomarker for GC detection. However, the detection accuracy was inconsistent among a series of quantitative analyses. These conflicting conclusions prompted us to employ this comprehensive and up-to-date research so that a conclusion on the diagnostic power of miR-106 for monitoring GC can be drawn. Meanwhile, we conducted an integrative functional analysis of miR-106 to understand the question why it could help distinguish GC patients from normal controls.

In the present study, we found that miR-106 achieved the overall pooled sensitivity of 0.71, specificity of 0.82, and AUC of 0.80, indicating a moderate overall accuracy. Circulating miR-106 as a more researched diagnostic marker in GC detection compared with other sample sources yielded a pooled sensitivity of 0.71, specificity of 0.82, and AUC of 0.81. Subgroup analysis indicated that sample type and sample size may influence the diagnostic accuracy. Specifically, it was revealed that plasma-based assays and miR-106 assays with a large sample size had significantly better overall diagnostic accuracy than serum-based ones and miR-106 assays with a small sample size, respectively. Meanwhile, although located in different chromosomes, miR-106a yielded a similar diagnostic accuracy compared with miR-106b.

Up to now, most attention on biomarker prediction in GC has been absorbed in single biomarkers. Actually, single biomarker is hard to reveal GC evolutionary process at the systems biology level as GC is a highly heterogeneous disease. On the contrary, combination markers may be more reliable with greater power for explaining the internal mechanisms of GC [36]. Therefore, we performed an analysis for miR-106-related combination markers in GC to investigate whether they were more powerful than miR-106 alone in detecting GC. It is worth noting that serum miR-106-related combination markers had a higher level of predictive power than serum miR-106 alone. However, it is difficult for us to conduct further investigations due to the limited number of studies enrolled in the analysis for miR-106-related combination markers.

We also performed integrative and comprehensive bioinformatics analysis to explore the function of miR-106 at the systems biology level. Most GO terms enriched by miR-106a and miR-106b target genes were both significantly associated with the regulation of transcription, apoptotic process at the BP level, basic cell structures at CC level along with the binding functions such as protein binding, DNA binding, transcription factor binding at MF level. Furthermore, the enriched KEGG pathways of miR-106a and miR-106b target genes were approximately the same, including pathways in cancer, p53 signaling pathway, cell cycle, TGF-beta signaling pathway and proteoglycans in cancer, which were highly associated with the occurrence and development of GC. Interestingly, although the different locations on chromosome, the targets of miR-106a and miR-106b fall into the similar functional modules, pathways or networks and then become more consistent when enriched to systems biology levels. In general, functionally concerned genes often emerge a coordinated expression to exert their roles in the same functional modules, indicating that miR-106a and miR-106b may have a synergistic effect in the initiation and progression of GC. The above results not only demonstrated the robustness of our study but explained why miR-106 could serve as a promising biomarker for detecting GC to some extent.

To further reveal the correlations among the target genes of miR-106a and miR-106b, we performed the PPI network analysis. Through PPI network construction, a series of hub genes were screened by three different network analysis methods. In our study, it was revealed that these key target nodes regulated by miR-106a and miR-106b both participated in FoxO signaling pathway, pathways in cancer, PI3K-Akt signaling pathway, cell cycle, and p53 signaling pathway. What’s more, module analysis of the PPI network revealed that the most significant modules of miR-106a and miR-106b targets network were both associated with p53 signaling pathway, pathways in cancer, FoxO signaling pathway, microRNAs in cancer, PI3K-Akt signaling pathway and cell cycle. PPI network analysis including hub genes identification and module analysis revealed the function of miR-106 again.

Based on above results, we found that several pathways were repeatedly mentioned in KEGG pathway analysis enriched by all the miR-106a and miR-106b targets, key hub targets and network modules, including p53 signaling pathway, pathways in cancer, FoxO signaling pathway, PI3 K-Akt signaling pathway and cell cycle. Pathways in cancer consist of several well-known signaling pathways including TGF-β, MAPK, Wnt and p53, which play important roles in cell apoptosis, proliferation, differentiation, invasion and metastasis. The well-studied p53 pathway, perhaps the most vital determinant of carcinogenesis, has been inextricably linked to establishment and progression of almost all types of cancer including GC [37, 38]. Cell cycle, another very important signaling pathway, contributes to the malignant progression of various human cancers including GC due to the involvement in cell growth, differentiation and apoptosis, as well as cancer development and metastasis [39]. Recent studies have proposed that the activation of the PI3K/Akt pathway may be responsible for the tumorigenesis by playing a pivotal role in control of cell cycle and survival of cell [40]. The information gathered so far indicates that FoxO signaling pathway could play vital roles in mediating apoptosis and thus determines cell death and survival [41]. In short, all the above pathways have been verified by the published literatures involved in the tumorigenesis and progression of GC, which may provide new ideas for the molecular mechanisms of miR-106 in GC.

Although mounting evidence from diagnostic tests indicated miR-106 as a promising GC marker, difficulty still remains for its application to clinical practice. There are several points we can do to optimize the miR-106 assay. Firstly, an appropriate standard cut-off value, consistent detection and normalization methods for miR-106 expression are required. Secondly, it was revealed from our results that plasma miR-106 may be a more powerful marker for detecting GC compared with serum miR-106. So plasma could be selected as the suitable sample source for further detection. Thirdly, as indicated in our study, sample size influenced the sensitivity and specificity. Larger sample size exhibited higher diagnostic accuracy. Thus, further large-scale prospective studies are warranted to develop integrative diagnostic models with more appropriate and better prediction capacity. Fourthly, although miR-106a and miR-106b are members of different paralogous clusters and located on different chromosomes, there has been some evidence in the literature, that both these two microRNAs can co-expressed in gastric tumor tissues. Based on our results, both miR-106a and miR-106b could be evaluated in diagnostic samples for diagnostic purposes of gastric cancer. In addition, combination biomarkers, which are combinations of several markers, have been shown to improve the prediction accuracy compared with single biomarker. According to our findings, single miR-106 was significant but not strong enough to undertake early diagnosis, while miR-106-related combination markers improved the diagnostic accuracy. The combination of miR-106 and other microRNAs may be the right way to solve the limited accuracy. Moreover, it has been reported that combination of protein-biomarkers and microRNAs may be an effective way to improve the diagnostic accuracy [42]. So more attempts are required for evaluating the combination biomarkers in the further study.

Our study had several important strengths. First, we carried out a relatively thorough systematic search and applied a comprehensive analytic approach to investigate the diagnostic power of miR-106 in patients with GC. Next, we evaluated the diagnostic value of miR-106-related combination markers in GC for the first time. It was suggested that the combination of miR-106 with other microRNAs improved the diagnostic accuracy, which may provide a novel potential tool for progress in a clinical context. Moreover, we performed integrative and comprehensive bioinformatics analysis to explore the function of miR-106 at the systems biology level, explaining the reason why miR-106 could be used in the diagnosis of GC. However, the power of our study was limited by a few factors. Firstly, most studies in the diagnostic tests enrolled healthy participants as controls and were not blind in design, which limits the diagnostic performance. Secondly, some key information including stage of cancer, sex proportions and age distributions was not known, so further analysis could not be carried out. Thirdly, there were no studies investigating the non-Asian population, which may cause potential heterogeneity from ethnicity. Fourthly, the sources of sample were inconsistent including plasma (n = 5), serum (n = 5), gastric juice (n = 1) and tissue (n = 1). Accordingly, subgroup analysis by specimen could not be performed for the limited individual sample size. In addition, the number of studies and sample sizes enrolled in the analysis for miR-106-related combination markers are limited, which make it difficult for us to conduct further investigations. As the miR-106 combination markers are all different in all six studies, it still remains an open question which should be combined with miR-106 for improving the diagnostic power.

Conclusion

Taken together, in this study, it is concluded that miR-106 is a useful biomarker for GC detection. Prospectively, combining miR-106 and other microRNAs may be considered as more powerful diagnostic tools for clinical application than individual miR-106. Integrative and comprehensive bioinformatics analysis was performed to explore the function of miR-106 at the systems biology level. Nonetheless, further large-scale prospective studies are needed to create integrative diagnostic models with more pronounced accuracy.

Abbreviations

GC:

gastric cancer

AUC:

area under the curve

GO:

gene ontology

KEGG:

Kyoto Encyclopedia of Genes and Genomes

PPI:

protein–protein interaction

QUADAS-2:

Quality Assessment of Diagnostic Accuracy Studies 2

TP:

true positive

FP:

false positive

FN:

false negative

TN:

true negative

PLR:

positive likelihood ratios

NLR:

negative likelihood ratios

DOR:

diagnostic odds ratio

CI:

confidence interval

SROC:

summary receiver operator characteristic

DAVID:

Database for Annotation, Visualization, and Integrated Discovery

STRING:

Search Tool for the Retrieval of Interacting Genes

MCODE:

Molecular Complex Detection

BP:

biological process

CC:

cellular component

MF:

molecular function

NA:

not available

qPCR:

quantitative real-time PCR

References

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

    Article  PubMed  Google Scholar 

  2. Leung WK, Wu MS, Kakugawa Y, et al. Screening for gastric cancer in Asia: current evidence and practice. Lancet Oncol. 2008;9(3):279–87.

    Article  PubMed  Google Scholar 

  3. Cancer Genome Atlas Research N. Comprehensive molecular characterization of gastric adenocarcinoma. Nature. 2014;513(7517):202–9.

    Article  CAS  Google Scholar 

  4. Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116(2):281–97.

    Article  PubMed  CAS  Google Scholar 

  5. Lin S, Gregory RI. MicroRNA biogenesis pathways in cancer. Nat Rev Cancer. 2015;15(6):321–33.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  6. Redova M, Sana J, Slaby O. Circulating miRNAs as new blood-based biomarkers for solid cancers. Future Oncol. 2013;9(3):387–402.

    Article  PubMed  CAS  Google Scholar 

  7. Koberle V, Pleli T, Schmithals C, et al. Differential stability of cell-free circulating microRNAs: implications for their utilization as biomarkers. PLoS ONE. 2013;8(9):e75184.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  8. Zhang Z, Li Z, Li Y, Zang A. MicroRNA and signaling pathways in gastric cancer. Cancer Gene Ther. 2014;21(8):305–16.

    Article  PubMed  CAS  Google Scholar 

  9. Tong AW, Nemunaitis J. Modulation of miRNA activity in human cancer: a new paradigm for cancer gene therapy? Cancer Gene Ther. 2008;15(6):341–55.

    Article  PubMed  CAS  Google Scholar 

  10. Li F, Liu J, Li S. MicorRNA 106b approximately 25 cluster and gastric cancer. Surg Oncol. 2013;22(2):e7–10.

    Article  PubMed  Google Scholar 

  11. Pan YJ, Zhuang Y, Zheng JN, Pei DS. MiR-106a: promising biomarker for cancer. Bioorg Med Chem Lett. 2016;26(22):5373–7.

    Article  PubMed  CAS  Google Scholar 

  12. Whiting PF, Rutjes AW, Westwood ME, et al. QUADAS-2: a revised tool for the quality assessment of diagnostic accuracy studies. Ann Intern Med. 2011;155(8):529–36.

    Article  PubMed  Google Scholar 

  13. Reitsma JB, Glas AS, Rutjes AW, Scholten RJ, Bossuyt PM, Zwinderman AH. Bivariate analysis of sensitivity and specificity produces informative summary measures in diagnostic reviews. J Clin Epidemiol. 2005;58(10):982–90.

    Article  PubMed  Google Scholar 

  14. Jones CM, Athanasiou T. Summary receiver operating characteristic curve analysis techniques in the evaluation of diagnostic tests. Ann Thorac Surg. 2005;79(1):16–20.

    Article  PubMed  Google Scholar 

  15. Higgins JP, Thompson SG, Deeks JJ, Altman DG. Measuring inconsistency in meta-analyses. BMJ. 2003;327(7414):557–60.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Petitti DB. Approaches to heterogeneity in meta-analysis. Stat Med. 2001;20(23):3625–33.

    Article  PubMed  CAS  Google Scholar 

  17. Deeks JJ, Macaskill P, Irwig L. The performance of tests of publication bias and other sample size effects in systematic reviews of diagnostic test accuracy was assessed. J Clin Epidemiol. 2005;58(9):882–93.

    Article  PubMed  Google Scholar 

  18. Chou CH, Shrestha S, Yang CD, et al. miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions. Nucleic Acids Res. 2018;46(D1):D296–302.

    Article  PubMed  Google Scholar 

  19. Gene Ontology C. The gene ontology (GO) project in 2006. Nucleic Acids Res. 2006;34(Database issue):D322–6.

    Google Scholar 

  20. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  21. Dennis G Jr, Sherman BT, Hosack DA, et al. DAVID: database for annotation, visualization, and integrated discovery. Genome Biol. 2003;4(5):P3.

    Article  PubMed  Google Scholar 

  22. Szklarczyk D, Morris JH, Cook H, et al. The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. 2017;45(D1):D362–8.

    Article  PubMed  CAS  Google Scholar 

  23. Shannon P, Markiel A, Ozier O, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. Tang Y, Li M, Wang J, Pan Y, Wu FX. CytoNCA: a cytoscape plugin for centrality analysis and evaluation of protein interaction networks. BioSystems. 2015;127:67–72.

    Article  PubMed  CAS  Google Scholar 

  25. Tsujiura M, Ichikawa D, Komatsu S, et al. Circulating microRNAs in plasma of patients with gastric cancers. Br J Cancer. 2010;102(7):1174–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Zhou H, Guo JM, Lou YR, et al. Detection of circulating tumor cells in peripheral blood from patients with gastric cancer using microRNA as a marker. J Mol Med. 2010;88(7):709–17.

    Article  PubMed  CAS  Google Scholar 

  27. Cai H, Yuan Y, Hao YF, Guo TK, Wei X, Zhang YM. Plasma microRNAs serve as novel potential biomarkers for early detection of gastric cancer. Med Oncol. 2013;30(1):452.

    Article  PubMed  CAS  Google Scholar 

  28. Cui L, Zhang X, Ye G, et al. Gastric juice MicroRNAs as potential biomarkers for the screening of gastric cancer. Cancer. 2013;119(9):1618–26.

    Article  PubMed  CAS  Google Scholar 

  29. Shiotani A, Murao T, Kimura Y, et al. Identification of serum miRNAs as novel non-invasive biomarkers for detection of high risk for early gastric cancer. Br J Cancer. 2013;109(9):2323–30.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  30. Zeng Q, Jin C, Chen W, et al. Downregulation of serum miR-17 and miR-106b levels in gastric cancer and benign gastric diseases. Chin J Cancer Res = Chung-kuo yen cheng yen chiu. 2014;26(6):711–6.

    PubMed  Google Scholar 

  31. Hou X, Zhang M, Qiao H. Diagnostic significance of miR-106a in gastric cancer. Int J Clin Exp Pathol. 2015;8(10):13096–101.

    PubMed  PubMed Central  CAS  Google Scholar 

  32. Li F, Guo Y, Liu J, Zhang R. The significance of elevated plasma expression of microRNA 106b ~ 25 clusters in gastric cancer. PLoS ONE. 2017;12(5):e0178427.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Wang N, Wang L, Yang Y, Gong L, Xiao B, Liu X. A serum exosomal microRNA panel as a potential biomarker test for gastric cancer. Biochem Biophys Res Commun. 2017;493(3):1322–8.

    Article  PubMed  CAS  Google Scholar 

  34. Yuan R, Wang G, Xu Z, et al. Up-regulated circulating miR-106a by DNA methylation promised a potential diagnostic and prognostic marker for gastric cancer. Anticancer Agents Med Chem. 2016;16(9):1093–100.

    Article  PubMed  CAS  Google Scholar 

  35. Di Leva G, Garofalo M, Croce CM. MicroRNAs in cancer. Ann Rev Pathol. 2014;9:287–314.

    Article  CAS  Google Scholar 

  36. Peng Q, Zhang X, Min M, Zou L, Shen P, Zhu Y. The clinical role of microRNA-21 as a promising biomarker in the diagnosis and prognosis of colorectal cancer: a systematic review and meta-analysis. Oncotarget. 2017;8(27):44893–909.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Bykov VJN, Eriksson SE, Bianchi J, Wiman KG. Targeting mutant p53 for efficient cancer therapy. Nat Rev Cancer. 2018;18(2):89–102.

    Article  PubMed  CAS  Google Scholar 

  38. Kunizaki M, Fukuda A, Wakata K, et al. Clinical significance of serum p53 antibody in the early detection and poor prognosis of gastric cancer. Anticancer Res. 2017;37(4):1979–84.

    Article  PubMed  Google Scholar 

  39. Otto T, Sicinski P. Cell cycle proteins as promising targets in cancer therapy. Nat Rev Cancer. 2017;17(2):93–115.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. Spangle JM, Roberts TM, Zhao JJ. The emerging role of PI3K/AKT-mediated epigenetic regulation in cancer. Biochem Biophys Acta. 2017;1868(1):123–31.

    PubMed  CAS  Google Scholar 

  41. Yadav RK, Chauhan AS, Zhuang L, Gan B. FoxO transcription factors in cancer metabolism. Semin Cancer Biol. 2018. https://doi.org/10.1016/j.semcancer.2018.01.004.

    Article  PubMed  Google Scholar 

  42. Cheng H, Zhang L, Cogdell DE, et al. Circulating plasma MiR-141 is a novel biomarker for metastatic colon cancer and predicts poor prognosis. PLoS ONE. 2011;6(3):e17745.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Authors’ contributions

QP and YS performed the statistical analysis, carried out the functional enrichment analysis, and drafted the manuscript. KL participated in the statistical analysis of meta-analysis. LZ and YS participated in functional enrichment analysis, and draft the manuscript. YZ conceived of the study, and participated in its design and coordination. All authors read and approved the final manuscript.

Acknowledgements

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The data supporting the conclusions of this article is within the article.

Consent for publication

All authors approve the publication of this manuscript.

Ethics approval and consent to participate

Not applicable.

Funding

This work was supported by the Second Affiliated Hospital of Soochow University Preponderant Clinic Discipline Group Project funding (XKQ2015005), Suzhou Science and Technology Development Program (SZS201509), Suzhou Clinical Medical Center Construction Project (Szzxj201503) and Jiangsu Medical Innovation Team (CXDT-37).

Publisher’s Note

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

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Yaqun Zhu.

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

Peng, Q., Shen, Y., Lin, K. et al. Comprehensive and integrative analysis identifies microRNA-106 as a novel non-invasive biomarker for detection of gastric cancer. J Transl Med 16, 127 (2018). https://doi.org/10.1186/s12967-018-1510-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-018-1510-y

Keywords