Integrated miRNA-risk gene-pathway pair network analysis provides prognostic biomarkers for gastric cancer

Introduction

Gastric cancer ranks as the fifth most common cancer and the third leading cause for cancer death, with diverse histological and molecular subtypes.^1,2 Over the past decades, due to the remarkable advances in diagnosis and treatment, the incidence and mortality of gastric cancer have decreased significantly in most countries globally.³ However, gastric cancer is still a challenge to be tackled, demanding intensive studies.

One of the prerequisites for the effective therapy of gastric cancer is to understand its prognosis status. Recently, several studies have been performed to investigate prognostic biomarkers for gastric cancer and their prognostic effects.^4,5 For example, overexpression of stomatin-like protein 2 is reported to be associated with invasion depth, tumor node metastasis stage, as well as lymph node and distant metastases, leading to poor prognosis in gastric cancer.⁶ Notch-1 has also been proven to be correlated with gastric cancer progression and to exert a predictive role for the poor clinical outcome of gastric cancer.⁷ Li et al⁸ have found by analyzing microarray data that a seven-miRNA signature (miR-10b, miR-21, miR-223, miR-338, let-7a, miR-30a-5p, and miR-126) is related to the survival and relapse of gastric cancer. Furthermore, analysis of oncogenic signaling pathways through the genomic and proteomic expression profiles has revealed the effect of their abnormalities on the prognosis of gastric cancer.⁹ For instance, the vascular endothelial growth factor (VEGF) pathway was reported as a possible prognostic factor for gastric cancer.¹⁰ The expression of VEGF is increased in gastric cancer cell lines, facilitating increases in their proliferation, motility, and adhesion.¹¹ Targeted inhibition of the VEGF pathway has been an important strategy to improve prognosis in gastric cancer therapy.¹² By constructing an interconnected network, Cheng et al¹³ showed that deregulation of WNT, NF-κB, and RTK-Ras-PIK3-AKT pathways may be involved in the development of gastric cancer by causing excessive cell proliferation and sustained angiogenesis. However, the prognostic biomarkers that can be applied in the clinic remain rare and thus further investigation of the etiology of gastric cancer is still necessary.

In this study, we aim to screen for highly potential prognostic molecules for gastric cancer by comprehensively analyzing genes, miRNAs, and pathways and thereafter integrating them, to construct an miRNA–target–pathway pairs network, the results of which may be believed to be more credible compared with the separate analysis of each. To our knowledge, this has not been previously reported.

Methods

mRNA and miRNA expression profiles

The mRNA microarray profiles of human gastric cancer samples were retrieved from Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/geo). After the first screening, 154 related studies underwent a second screening using the following criteria: 1) the studies should have involved no less than ten samples; 2) tissue samples instead of cell samples should have been studied; 3) gastric tissue samples should be the only research topics; 4) the cancer should not have been caused by chemotherapy and other artificial stimulations; 5) the samples should have only undergone the treatment for microarray analysis; 6) normal control samples or normal adjacent tissues should be available; 7) the platforms should include data on >10,000 genes; 8) the annotation of the samples should be available; 9) the data should not be two-panel data; and 10) the profile should have many overlapping genes with the other profiles. Next, different subtypes of gastric cancer samples were merged as the disease samples, while the normal tissue and adjacent tissue samples were recognized as the control samples. Finally, a total of 13 data sets, comprising 497 disease samples and 425 control samples, were selected for further analysis (Table 1).

Table 1 Summary of mRNA expression profiles for gastric cancer

In addition, the miRNA expression profile for gastric cancer (accession number GSE26595) was also downloaded from GEO. A total of 60 gastric samples and eight normal samples were available, based on the GPL8197 platform.

Screening for differentially expressed genes and miRNAs

As the mRNA microarray data were merged from 13 individual data sets, MetaDE¹⁴ in R language was used to identify the genes whose expressions differed between cancer samples and control samples. Genes with adjusted P-value (false discovery rate or FDR) <0.01 were considered differentially expressed genes (DEGs). For the screening of the differentially expressed miRNAs (DEmiRs), limma¹⁵ in R language was utilized with the threshold of FDR <0.05.

Potential gastric cancer-related genes

The targets of the DEmiRs were first identified from the following four databases: PicTar (http://pictar.mdc-berlin.de/), TargetScan (http://www.targetscan.org/), PITA (http://genie.weizmann.ac.il/pubs/mir07/mir07_data.html), and miRanda (http://www.microrna.org/microrna/getDownloads.do). The targets of a certain miRNA were selected as those predicted by at least three of these four databases. Additionally, miRecords,¹⁶ miRTarBase,¹⁷ and TarBase¹⁸ were also searched to identify the verified targets of the DEmiRs. Targets from these screenings were merged as the final targets of the DEmiRs. Finally, the targets that overlapped with the DEGs were recognized as the potential gastric cancer-related genes.

Functional enrichment

The potential gastric cancer-related genes were subjected to functional interpretation using the online tools of the bioinformatics resource Database for Annotation, Visualization and Integrated Discovery for function and pathway enrichment. The terms with FDR <0.05 were considered the significant terms of these genes.

Pathway expression profile

The Kyoto Encyclopedia of Genes and Genomes (KEGG; http://www.genome.jp/kegg/pathway.htm) is a database for the functional analysis of genes. All the human pathways and genes in these pathways were downloaded from KEGG for further analysis. Next, the expression level of each pathway in the samples was calculated by computing the expression levels of its genes in each sample, to obtain the pathway expression profile. The expression level of the pathway was calculated using the following formula:

Pathik = Mediank (g1, g2, g3,..., gn)

(1)

Here, Path_ik is the expression level of pathway i in sample k; g₁, g₂, g₃,..., g_n are the expression values in sample k of all genes in pathway i. Path_ik remains the same when there are observation errors, abnormal values, and heteroskedastic variances; therefore, it is chosen as the repression of expression levels of pathways.

Screening of co-pathway pairs

To obtain a comprehensive understanding about the roles of abnormal pathways in gastric cancer, the co-pathway pairs that are closely connected were identified by calculating their expression correlations in all samples. The correlation was measured by the value of the Pearson correlation coefficient, which was computed using the following formula:

(2)

Here, P_X,Y is the Pearson correlation coefficient between pathways X and Y in all samples. and are the average expression values of pathways X and Y, respectively.

Subsequently, the co-pathway pairs highly related to gastric cancer were screened using Equation 3 based on the hypergeometric distribution algorithm:

(3)

Here, n represents the common gastric cancer-related genes in the two pathways. N represents all the genes of the co-pathway pairs, M represents cancer-related genes of the co-pathway pairs and C represents the combination of genes of the co-pathway pairs. The co-pathway pairs with P<0.05 were considered to be the significant gastric cancer-related co-pathway pairs.

Construction of miRNA–risk gene–pathway pair integrated network

The common cancer-related genes in the co-pathway pairs were abstracted to merge with the target genes of DEmiRs from the miRNA–risk gene pairs. Subsequently, an integrated network of miRNA–risk gene–pathway pair was constructed, for which the network analysis was carried out based on the Cytoscape software¹⁹ to calculate the network’s topological features. Additionally, classification of the most important miRNAs in the network was achieved by cross-validation by calculating the receiver operating characteristic (ROC) value.

Results

Screening for DEGs and DEmiRs

The expression profiles of the normal and gastric cancer samples from 13 studies were analyzed for the DEGs using the MetaDE method. Under the threshold of FDR <0.01, a total of 1,260 DEGs were identified, including 1,139 upregulated ones and 121 downregulated ones (Table 2). The DEmiRs were identified using limma package; totally, 144 DEmiRs were obtained, comprising 69 upregulated and 75 downregulated ones (Table 3).

Table 2 Top ten differentially expressed genes Abbreviation: FDR, false discovery rate.

Table 3 Top ten differentially expressed miRNAs Abbreviation: FDR, false discovery rate.

miRNA–target pairs

In all, a total of 8,642 predicted target genes for DEmiRs were identified from PicTar, PITA, miRanda, and TargetScan (Figure 1A), while a total of 912 verified target genes were identified from miRecords, miRTarBase, and Tarbase (Figure 1B). Finally, the common genes of potential targets and verified targets were selected for the eventual miRNA–target pairs (Figure 1).

Figure 1 The numbers of target genes for DEmiRs. Notes: In all, (A) a total of 8,642 predicted target genes for DEmiRs were identified from PicTar, PITA, miRanda, and TargetScan, while (B) a total of 912 verified target genes were identified from miRecords, miRTarBase, and Tarbase. Finally, (C) the common genes of potential targets and verified targets were selected for the eventual miRNA–target pairs. Abbreviation: DEmiR, differentially expressed miRNA.

Potential gastric cancer-related genes

The eventual targets for DEmiRs were merged with DEGs to select the potential gastric cancer-related genes (risk genes), and the miRNA–risk gene pairs were believed to play vital roles in the development of gastric cancer. Totally, 570 miRNA–risk gene pairs were obtained with 35 miRNAs and 336 risk genes (Figure 2).

Figure 2 The number of risk genes. Abbreviations: DEG, differentially expressed gene; DEmiR, differentially expressed microRNA.

Functional enrichment

The enrichment analysis showed that DNA metabolic process, trabecula formation, and epidermis development were the most dysfunctional biological processes; nuclear part, nuclear lumen, and nucleoplasm were the most dysfunctional cellular components; adenyl ribonucleotide binding, calcium ion binding, and pyrophosphatase activity were the most dysfunctional molecular functions; extracellular matrix (ECM)–receptor interaction, focal adhesion, and cell cycle were the most dysfunctional pathways (Figure 3).

Figure 3 Enriched functions and pathways of risk genes. Notes: (A) Biological process; (B) cellular components; (C) molecular functions; (D) pathways. The terms are listed in the Y-axis, while the significances are listed on the X-axis. The higher the log(P-value), the more significant is the term. Abbreviation: ECM, extracellular matrix.

Pathway expression profile and co-pathway pairs

Totally, 239 human pathways were extracted from KEGG to construct the pathway expression profile (Figure 4). According to the Pearson correlation coefficients, 45 co-pathway pairs were selected as the most closely related pairs. Next, 15 of these co-pathway pairs were chosen as the most significant ones connected with risk genes. The risk gene–pathway pairs were next combined with the miRNA–risk gene pairs for the construction of miRNA–risk gene–pathway pair network. The integrated network consisted of 1,398 interaction pairs and 480 nodes (miRNA, risk genes, co-pathway pairs) (Figure 5).

Figure 4 Heat map of the pathway expression profile. Notes: The X-axis refers to the pathways, and the Y-axis refers to the samples. The colors from green to red indicate the pathway expression level from low to high. Abbreviations: GABA, gamma-aminobutyric acid; TGF-β, transforming growth factor-β; VEGF, vascular endothelial growth factor; cAMP, 3′,5′-cyclic adenosine monophosphate; TNF, tumor necrosis factors; GnRH, gonadotropin-releasing hormone; TRP channels, transient receptor potential channels; AMPK, AMP-activated protein kinase; NOD, nucleotide-binding, oligomerization domain; RIG-I, retinoic acid-inducible gene 1; NF-κB, nuclear factor kappa-light-chain-enhancer of activated B cells; ECM, extracellular matrix; PI3K, phosphatidylinositol 3-kinase; ABC transporters, ATP-binding cassette transporters; HTLV-I, human T-cell lymphotropic virus type 1; cGMP, 3′,5′-cyclic guanosine monophosphate; PPAR, peroxisome proliferator-activated receptors; Jak-STAT, janus kinase-signal transducer and activator of transcription; MAPK: mitogen-activated protein kinase; mTOR, mammalian target of rapamycin.

Figure 5 miRNA–risk gene–pathway pair complex network. Notes: The rhombus refers to the miRNA, the rotundity refers to the mRNA, and yellow refers to the co-pathway pairs. Red indicates upregulation and green indicates downregulation. Abbreviation: ECM, extracellular matrix.

Using Cytoscape, the topological features of the network were calculated (Table 4). The co-pathway pairs “Focal adhesion/ECM–receptor interaction”, “Proteoglycans in cancer/MicroRNAs in cancer”, and “3′,5′-cyclic guanosine monophosphate (cGMP)–protein kinase G (PKG) signaling pathway/3′,5′-cyclic adenosine monophosphate (cAMP) signaling pathway” were closely associated with risk genes. cGMP–PKG signaling pathway/cAMP signaling pathway pair was the co-pathway pair with the highest degree. CALM2 was the most important gene in the integrated network, with the highest connecting degree with miRNAs and pathway pairs (Figure 6; Table 5), such as cGMP–PKG signaling pathway/cAMP signaling pathway, Estrogen signaling pathway/Melanogenesis, and Ras signaling pathway/Rap1 signaling pathway. Moreover, CALM2 was regulated by hsa-miR-193a-3p, hsa-miR-1246, hsa-miR-193b, hsa-miR-19b, hsa-miR-196b, and hsa-miR-181b. The topological features of the top 15 miRNAs with highest degrees are listed in Table 6, with hsa-miR-19b and hsa-miR-181b as the top two miRNAs. As these two miRNAs had differential expressions in gastric cancer samples compared with the normal samples (Figure 7A), their ability to classify the disease samples was detected by cross-validation. The results showed that hsa-miR-19b had an ROC value of 0.7281, while hsa-miR-181b had an ROC value of 0.9125 (Figure 7B).

Table 4 Network topological features of 15 pathways Abbreviations: cAMP, 3′,5′-cyclic adenosine monophosphate; cGMP, 3′,5′-cyclic guanosine monophosphate; ECM, extracellular matrix; PKG, protein kinase G; MAPK, mitogen-activated protein kinase; HIF-1, hypoxia-inducible factor-1.

Figure 6 Interactive nodes of risk gene CALM2 in the integrated network of miRNA–risk gene–pathway pair. Notes: The rhombus refers to the miRNA, the rotundity refers to the mRNA, and yellow refers to the co-pathway pairs. Red indicates upregulation and green indicates downregulation. Abbreviation: ECM, extracellular matrix.

Table 5 Network topological features of the top 15 risk genes with highest degrees

Table 6 Network topological features of the top 15 differentially expressed miRNAs with highest degrees Abbreviation: DEmiR, differentially expressed miRNA.

Figure 7 Expressions and classification efficacy of two miRNAs with highest degrees in the integrated network. Notes: (A) Expression levels between gastric cancer samples and normal samples; (B) ROC curves of the miRNAs for the gastric cancer samples. Abbreviation: ROC, receiver operating characteristic.

Discussion

The pathways with higher degrees in the integrated network were considered to be the pathways significantly associated with the occurrence and development of gastric cancer, and the co-pathway pair “Focal adhesion/ECM–receptor interaction” was found to interact with the DEGs of gastric cancer samples. Focal adhesion kinase regulates the intercellular signaling network, which plays a key role in the differentiation and metastasis of cells. Interruption of this pathway inhibits the spread of gastric cancer cells.^20,21 It is reported that the connection of Focal adhesion/ECM–receptor interaction induces cascade reactions, and they are both involved in gastric carcinoma.^22–24

Calmodulin is a highly conserved Ca²⁺-binding protein in the calcium signaling process, participating in signaling pathways including motility, proliferation, and cell cycle progression and thus plays important roles in cancer development.²⁵ CALM2, a risk gene with the highest degree in our integrated network, encodes calmodulin 2, which has also been found to be overexpressed in mammary cancer cells.^26–28 High expression of CALM2 may promote the activation of a calcium/calmodulin-dependent kinase 2 (CaMKII), which then enhances the metastasis of gastric cancer by upregulating NF-κB and Akt-mediated matrix metalloproteinase-9 production.^29,30

miRNAs are evolutionarily conserved small noncoding RNAs, which have shown their roles in the onset and progression of cancers.³¹ Many miRNAs have been proven to be associated with the prognosis of various solid tumors.³² As a member in the miR-17-92 cluster, miR-19b is recognized as an oncomiR influencing multiple aspects of the phenotypes of gastric cancer.³³ The expression of miR-19b is high in gastric cancer tissues and it is significantly associated with its metastasis.³⁴ In gastric cancer, the expression of miR-181 is upregulated.³⁵ Guo et al³⁶ have found that miR-181b is a potential therapeutic target for gastric neoplasms. Moreover, miR-181b also holds great potential to be used as a prognostic biomarker for gastric cancer in the late stage.³⁷ In this study, these two miRNAs showed satisfactory classification of the gastric cancers, with ROC values of 0.7281 and 0.9125, which further confirmed their prognostic roles in gastric cancer. Additionally, this consistency also proved the accuracy of the miRNA–targets–pathway pair network.

It is reported that miRNAs may function in cancer by repressing the expression of target genes via binding to the 3′-untranslated region. A previous study³⁸ has observed that CALM2 is a target anti-correlated with the DEmiR hsa-miR-130b. However, no researchers have identified the correlations between CALM2 and hsa-miR-19b/hsa-miR-181b. Considering the roles of these two genes and the association between these and CALM2, we assumed that CALM2 might be an important prognostic molecule for gastric cancer.

Conclusion

By constructing the miRNA–target–pathway pair network, our findings indicate that Focal adhesion/ECM–receptor interaction pathway, CALM2, hsa-miR-19b, and hsa-miR-181b might be used as potential prognostic targets for gastric cancer. However, further experimental studies are needed to confirm our conclusion.

Acknowledgment

This study was supported by the China Postdoctoral Science Foundation (grant number: 2014M552576) and the National Natural Science Foundation of China (grant number: 81402433).

Disclosure

The authors report no conflicts of interest in this work.

---

References

1.		Wang K, Yuen ST, Xu J, et al. Whole-genome sequencing and comprehensive molecular profiling identify new driver mutations in gastric cancer. Nat Genet. 2014;46(6):573–582.
2.		Fock K. Review article: the epidemiology and prevention of gastric cancer. Aliment Pharmacol Ther. 2014;40(3):250–260.
3.		Ferro A, Peleteiro B, Malvezzi M, et al. Worldwide trends in gastric cancer mortality (1980–2011), with predictions to 2015, and incidence by subtype. Eur J Cancer. 2014;50(7):1330–1344.
4.		Zhao W, Zhu H, Zhang S, et al. Trop2 is overexpressed in gastric cancer and predicts poor prognosis. Oncotarget. 2016;7(5):6136–6145.
5.		Ying J, Xu Q, Liu B, Zhang G, Chen L, Pan H. The expression of the PI3K/AKT/mTOR pathway in gastric cancer and its role in gastric cancer prognosis. Onco Targets Ther. 2015;8:2427–2433.
6.		Liu D, Zhang L, Shen Z, et al. Increased levels of SLP-2 correlate with poor prognosis in gastric cancer. Gastric Cancer. 2013;16(4):498–504.
7.		Zhang H, Wang X, Xu J, Sun Y. Notch1 activation is a poor prognostic factor in patients with gastric cancer. Br J Cancer. 2014;110(9):2283–2290.
8.		Li X, Zhang Y, Zhang Y, Ding J, Wu K, Fan D. Survival prediction of gastric cancer by a seven-microRNA signature. Gut. 2010;59(5):579–585.
9.		Lim SM, Lim JY, Cho JY. Targeted therapy in gastric cancer: personalizing cancer treatment based on patient genome. World J Gastroenterol. 2014;20(8):2042–2050.
10.		Penon D, Cito L, Giordano A. Novel findings about management of gastric cancer: a summary from 10(th) IGCC. World J Gastroenterol. 2014;20(27):8986–8992.
11.		Liu K, Chen H, You Q, Shi H, Wang Z. The siRNA cocktail targeting VEGF and HER2 inhibition on the proliferation and induced apoptosis of gastric cancer cell. Mol Cell Biochem. 2014;386(1–2):117–124.
12.		Abdel-Rahman O. Targeting vascular endothelial growth factor (VEGF) pathway in gastric cancer: preclinical and clinical aspects. Crit Rev Oncol Hematol. 2015;93(1):18–27.
13.		Cheng L, Yang S, Yang Y, et al. Global gene expression and functional network analysis of gastric cancer identify extended pathway maps and GPRC5A as a potential biomarker. Cancer Lett. 2012;326(1):105–113.
14.		Li J, Tseng GC. An adaptively weighted statistic for detecting differential gene expression when combining multiple transcriptomic studies. Ann Appl Stat. 2011;5:994–1019.
15.		Patnaik SK, Kannisto E, Knudsen S, et al. Evaluation of microRNA expression profiles that may predict recurrence of localized stage I non-small cell lung cancer after surgical resection. Cancer Res. 2010;70(1):36–45.
16.		Xiao F, Zuo Z, Cai G, Kang S, Gao X, Li T. miRecords: an integrated resource for microRNA-target interactions. Nucleic Acids Res. 2009;37(Database issue):D105–D110.
17.		Hsu SD, Lin FM, Wu WY, et al. miRTarBase: a database curates experimentally validated microRNA-target interactions. Nucleic Acids Res. 2011;39(Database issue):D163–D169.
18.		Sethupathy P, Corda B, Hatzigeorgiou AG. TarBase: a comprehensive database of experimentally supported animal microRNA targets. RNA. 2006;12(2):192–197.
19.		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–2504.
20.		Su GQ, Zhang FX, Mao HH, et al. Research of shRNAmir inhibitory effects towards focal adhesion kinase expression in the treatment of gastric cancer. Oncol Lett. 2015;9(2):595–603.
21.		Cheng YJ, Zhu ZX, Zhou JS, et al. Silencing profilin-1 inhibits gastric cancer progression via integrin beta1/focal adhesion kinase pathway modulation. World J Gastroenterol. 2015;21(8):2323–2335.
22.		Hu K, Chen F. Identification of significant pathways in gastric cancer based on protein-protein interaction networks and cluster analysis. Genet Mol Biol. 2012;35(3):701–708.
23.		Yang S, Shin J, Park KH, et al. Molecular basis of the differences between normal and tumor tissues of gastric cancer. Biochim Biophys Acta. 2007;1772(9):1033–1040.
24.		Liu P, Wang X, Hu CH, Hu TH. Bioinformatics analysis with graph-based clustering to detect gastric cancer-related pathways. Genet Mol Res. 2012;11(3):3497–3504.
25.		Coticchia CM, Revankar CM, Deb TB, Dickson RB, Johnson MD. Calmodulin modulates Akt activity in human breast cancer cell lines. Breast Cancer Res Treat. 2009;115(3):545–560.
26.		Król M, Skierski J, Rao N, et al. Transcriptomic profile of two canine mammary cancer cell lines with different proliferative and anti-apoptotic potential. J Physiol Pharmacol. 2009;60(Suppl 1):95–106.
27.		Calaluce R, Gubin MM, Davis JW, et al. The RNA binding protein HuR differentially regulates unique subsets of mRNAs in estrogen receptor negative and estrogen receptor positive breast cancer. BMC Cancer. 2010;10:126.
28.		Rust R, Visser L, Leij J, et al. High expression of calcium-binding proteins, S100A10, S100A11 and CALM2 in anaplastic large cell lymphoma. Br J Haematol. 2005;131(5):596–608.
29.		Liu Z, Han G, Cao Y, Wang Y, Gong H. Calcium/calmodulin-dependent protein kinase II enhances metastasis of human gastric cancer by upregulating nuclear factor-κb and Akt-mediated matrix metalloproteinase-9 production. Mol Med Rep. 2014;10(5):2459–2464.
30.		Subbannayya Y, Syed N, Barbhuiya MA, et al. Calcium calmodulin dependent kinase kinase 2 – a novel therapeutic target for gastric adenocarcinoma. Cancer Biol Ther. 2015;16(2):336–345.
31.		Garofalo M, Di Leva G, Croce CM. MicroRNAs as anti-cancer therapy. Curr Pharm Des. 2014;20(33):5328–5335.
32.		Ueda T, Volinia S, Okumura H, et al. Relation between microRNA expression and progression and prognosis of gastric cancer: a microRNA expression analysis. Lancet Oncol. 2010;11:136–146.
33.		Wang F, Li T, Zhang B, et al. MicroRNA-19a/b regulates multidrug resistance in human gastric cancer cells by targeting PTEN. Biochem Biophys Res Commun. 2013;434(3):688–694.
34.		Wu Q, Yang Z, An Y, et al. miR-19a/b modulate the metastasis of gastric cancer cells by targeting the tumour suppressor MXD1. Cell Death Dis. 2014;5:e1144.
35.		Wang J, Wang Q, Liu H, Hu B, Zhou W, Cheng Y. MicroRNA expression and its implication for the diagnosis and therapeutic strategies of gastric cancer. Cancer Lett. 2010;297(2):137–143.
36.		Guo JX, Tao QS, Lou PR, Chen XC, Chen J, Yuan GB. miR-181b as a potential molecular target for anticancer therapy of gastric neoplasms. Asian Pac J Cancer Prev. 2012;13(5):2263–2267.
37.		Jiang J, Zheng X, Xu X, et al. Prognostic significance of miR-181b and miR-21 in gastric cancer patients treated with S-1/oxaliplatin or doxifluridine/oxaliplatin. PLoS One. 2011;6(8):e23271.
38.		Yang J, Zeng Y. Identification of miRNA-mRNA crosstalk in pancreatic cancer by integrating transcriptome analysis. Eur Rev Med Pharmacol Sci. 2015;19(5):825–834.