Downregulated miR-383-5p contributes to the proliferation and migration of gastric cancer cells and is associated with poor prognosis

Abstract

Aim

The study aims to identify differentially expressed microRNAs (DEMs) in gastric cancer (GC) and explore the expression, prognosis and downstream regulation role of miR-383-5p in GC.

Methods

The GC miRNA-Seq and clinical information were downloaded from Firebrowse which stores integrated data sourced from The Cancer Genome Atlas database. The DEMs were identified with limma package in R software at the cut-off criteria of P < 0.05 and |log2 fold change| > 1.0 (|log2FC| > 1.0). The expression of miR-383-5p in GC cell lines and 54 paired GC tissues was measured by quantitative real-time polymerase chain reaction (qRT-PCR). The overall survival curve of miR-383-5p and the association between its expression and clinicopathological features were explored. Wound healing and cell counting kit-8 assays were performed to investigate the capacity of miR-383-5p in cell proliferation and migration. The downstream target genes were predicted by bioinformatics tools (miRDB, TargetScan and starBase). The consensus target genes were selected for gene functional enrichment analysis by FunRich v3.0 software. The luciferase reporter assay was performed to verify the potential targeting sites of miR-383-5p on lactate dehydrogenase A (LDHA).

Results

A total of 21 down-regulated miRNAs (including miR-383-5p) and 202 up-regulated miRNAs were identified by analyzing GC miRNA-Seq data. Survival analysis found that patients with low miR-383-5p expression had a shorter survival time (median survival time 21.1 months) than those with high expression (46.9 months). The results of qRT-PCR indicated that miR-383-5p was downregulated in GC cell lines and tissues, which was consistent with miRNA-Seq data. The expression of miR-383-5p was significantly associated with tumor size and differentiation grade. Besides, overexpression of miR-383-5p suppressed GC cells proliferation and migration. A total of 49 common target genes of miR-383-5p were obtained by bioinformatics tools and gene functional enrichment analysis showed that these predicted genes participated in PI3K, mTOR, c-MYC, TGF-beta receptor, VEGF/VEGFR and E-cadherin signaling pathways. The data showed that expression of miR-383-5p was negatively correlated with target LDHA (r = −0.203). Luciferase reporter assay suggested that LDHA was a target of miR-383-5p.

Conclusion

The present study concluded that miR-383-5p was downregulated and may act as a tumor suppressor in GC. Furthermore, its target genes were involved in important signaling pathways. It could be a prognostic biomarker and play a vital role in exploring the molecular mechanism of GC.

Introduction

Gastric cancer (GC) is a common malignancy of human digestive system which has high incidence and mortality worldwide (Bray et al., 2018). In China, both the incidence and mortality of GC rank second among malignant neoplasms (Chen et al., 2016). The number of Chinese patients with GC has increased over the years, which seriously threatens the health of people. The majority of Chinese GC patients are diagnosed at advanced stage and have an unsatisfactory 5-year overall survival rate (Correa, 2013). Most studies have found that the development of GC is associated with multiple factors, such as irregular diet, genetic and epigenetic influence (Carcas, 2014).

MicroRNAs (miRNAs, 20–24 nucleotides in length) are a series of non-coding RNAs and play important roles in the regulation of gene expression at post-transcriptional level (Zhang, Wang & Gemeinhart, 2013). Mechanistically, miRNAs negatively regulate gene expression through binding to sites in the 3′-untranslated regions of messenger RNAs (Zen & Zhang, 2012). Over the past years, mounting studies have confirmed that miRNAs can act as oncogenes or anti-oncogenes in the initiation and development of GC by regulating the downstream target genes (Chen et al., 2018; Hui et al., 2018; Wang et al., 2019b). Thus, exploring the expression and regulation role of miRNAs may be in favor of uncovering the tumorigenesis mechanism of GC.

In present study, we first identified that miR-383-5p was down-regulated in GC tissues by analyzing GC miRNA-Seq data. In order to confirm this finding, quantitative real-time polymerase chain reaction (qRT-PCR) experiment was performed to measure the expression of miR-383-5p in GC cell lines and tissues. Kaplan–Meier survival analysis also found that patients in high miR-383-5p expression group have longer overall survival time than those in low miR-383-5p expression group. Overexpression of miR-383-5p suppressed GC cells proliferation and migration. All the results showed that miR-383-5p was downregulated and it may play an anti-oncogene role in GC. The potential target genes of miR-383-5p were predicted through online bioinformatics tools. The functional enrichment analysis of target genes indicated that miR-383-5p may take part in PI3K, mTOR, c-MYC, TGF-beta receptor, VEGF/VEGFR and E-cadherin signaling pathways. LDHA, one of the 49 common target genes of miR-383-5p, was selected for target validation. The data showed that the expression of miR-383-5p was negatively correlated with LDHA (r = −0.203). Furthermore, a luciferase reporter assay suggested that the luciferase activity in wild-type (WT) LDHA-3′UTR group was significantly decreased by the miR-383-5p mimics, and there were no differences in the mutant LDHA-3′UTR group. Above all, miR-383-5p can be a meaningful target in understanding the potential molecular mechanism of GC tumorigenesis and progression.

Materials and Methods MiRNA-Seq data and clinical information

The GC clinical information and miRNA-Seq data, which contain 389 cancer tissue samples and 41 gastric normal tissue samples, were downloaded from Firebrowse website (The Cancer Genome Atlas (TCGA) data version 2016_01_28). At the cut-off criterion of |log2 fold change| > 1.0 (|log2FC| > 1.0) and P < 0.05, the differentially expressed microRNAs (DEMs) were identified using the limma package in R software. The follow-up days and vital status of patients were extracted from clinical information data. Patients meeting the following criteria were included for overall survival: (1) patients have integrated follow-up days and vital status; (2) patients have both follow-up days and expression value of miR-383-5p. In total, 382 patients were respectively divided into low and high expression group according to the median value of miR-383-5p expression. The overall survival curve of low and high miR-383-5p expression groups were analyzed with the method of Kaplan–Meier and log-rank test.

GC cell lines culture

Human GC cell lines (SGC-7901, BGC-823, MGC-803, and MKN-45) and a normal gastric mucous membrane cell line (GES-1) were purchased from the Institute of Biochemistry and Cell Biology of the Chinese Academy of Sciences (Shanghai, China), and all cell lines were cultivated in RPMI 1640 medium (GIBCO-BRL) supplemented with 10% fetal bovine serum (FBS; Gibco, Grand Island, NY, USA), 100 U/mL penicillin and 100 mg/mL streptomycin. All cells were cultured in humidified air at 37 °C and 5% CO₂.

Patient tissues collection

All the GC tissues and the corresponding adjacent normal tissues were collected from 54 patients who received surgical resection at the No. 967 Hospital of PLA Joint Logistics Support Force and the Northern Theater Command General Hospital. Tissues were histologically confirmed and immediately stored at −80 °C after resection. The clinicopathological features of 54 GC patients were recorded and preserved. All patients signed the informed consent and this study was approved by the Research Ethics Committee of Jinzhou Medical University (IRB No: EC-2018-JZ-016).

RNA extraction and qRT-PCR assays

The TRIzol reagent (Invitrogen, Carlsbad, CA, USA) was utilized for extracting total RNA of tissues and cells. Reverse Transcription Kit (GenePharma, Shanghai, China) was used for obtaining cDNA reverse transcribed from RNA. qRT-PCR assay was performed with SYBR-Green Hairpin-it™ MicroRNAs Kit (GenePharma, Shanghai, China), which was conducted on ABI 7500 FAST Real-Time PCR System. The expression level was determined using 2^−ΔΔCt method and normalized to U6. All the sequences of primers used in present study were summarized in Table 1.

Table 1: Sequence of primers used for PCR.

Name	Sequence (5′-3′)
miR-383-5p (RT)	GTCGTATCCAGTGCGTGTCGTGGAGTCGGCAATTGCACTGGATACGACAGCCAC
miR-383-5p (forward)	GGGAGATCAGAAGGTGATTGTGGCT
miR-383-5p (reverse)	CAGTGCGTGTCGTGGAGT
U6 (forward)	CTCGCTTCGGCAGCACA
U6 (reverse)	AACGCTTCACGAATTTGCGT

DOI: 10.7717/peerj.7882/table-1

RNA oligonucleotide and cell transfection

The miR-383-5p mimics and mimics NC were designed and synthesized by GenePharma Co., Ltd (Shanghai, China). According to the manufacturer’s protocols, the GC cells were transfected using Lipofectamine™ 3000 reagent (Invitrogen, Carlsbad, CA, USA; Thermo Fisher Scientific, Waltham, MA, USA). All the cells were cultivated for 48 h after transfection.

Cell proliferation assay

After transfection, the MGC-803 or MKN-45 cells were seeded into a 96-well plate with 4 × 10³ cells per well in triplicate. At 0, 24, 48, and 72 h, each well was added with 10 μL of cell counting kit-8 (CCK-8) reagent (Dojindo, Kumamoto, Japan) and then incubated at 37 °C for 3 h. The absorbance at 450 nm was measured using a spectrophotometer.

Wound healing assay

The MGC-803 or MKN-45 cells were seeded into six-well plate. When the cells were cultured to a density of 90%, a 100 μL pipette tip was used to draw a straight wound. Then, the cells were cultured with serum-free medium in the humidified incubator. At 0 and 48 h, an inverted microscope was utilized to visualize the wound healing and photograph. Before observing the healing status, mitomycin C (10 ug/mL) was added into each well for 2 h to exclude the influence of proliferation on cell migration.

Target genes prediction of miR-383-5p and functional enrichment analysis

The bioinformatics websites of TargetScan (Agarwal et al., 2015) (http://www.targetscan.org/vert_71/), miRDB (Wang, 2016) (http://www.mirdb.org/) and starBase (Li et al., 2014) (http://starbase.sysu.edu.cn/) were applied for predicting potential target genes of miR-383-5p. The consensus results of the three tools were selected for further analysis. Genes functional enrichment analysis was performed by FunRich v3.0 software which was a widely used tool for the gene functional enrichment and interaction network analysis (Pathan et al., 2015). All the procedures were conducted according to official protocols and default parameters.

Luciferase reporter assay

The luciferase reporter vectors containing WT or mutant (MUT) 3′UTR of LDHA and miR-383-5p mimics or negative control were transfected into MKN-45 cells. Luciferase activity was tested using the Dual-Luciferase Reporter Assay System (Promega, Madison, WI, USA). Activities were normalized to Renilla luciferase. All experiments were performed three times.

Statistical analysis

MiRNA-Seq data was processed by the limma (Ritchie et al., 2015) package in R software. The survival curve was described by Kaplan–Meier survival plot and analyzed with log-rank test. The differences between the two groups were analyzed by paired or unpaired Student’s t-test. The Chi-square test was used for exploring the association between miR-383-5p expression and clinical features (such as: age, gender, tumor size, lymph node metastasis, TNM stage, and differentiation grade). P < 0.05 was recognized as statistically significant and all statistical analysis were conducted by IBM SPSS software 19.0.

Results Identification of DEMs in GC

At the cut-off criterion of |log2FC| > 1.0 and P < 0.05, 223 DEMs were identified by screening GC miRNA-Seq data. A volcano plot was drawn to visualize the 21 down-regulated and 202 up-regulated miRNAs (Fig. 1). The top 20 of down- and up-regulated miRNAs ranked by FC was listed in Table 2. We found that miR-383-5p was significantly down-regulated with a log2FC of −1.12, which indicated that it may act as a tumor-suppressor in GC.

Figure 1: The volcano plot. The red spots represent 202 up-regulated miRNAs and the green spots represent 21 down-regulated miRNAs (blue imaginary lines represent Fold change: ±2 and P-value: 0.05). DOI: 10.7717/peerj.7882/fig-1

Table 2: A total of 42 DEMs identified between GC and adjacent normal tissues.

Down-regulated	logFC	P-value	Up-regulated	logFC	P-value
miR-1-3p	−3.13	1.56E-19	miR-196a-5p	5.55	3.01E-39
miR-133a-3p	−3.13	9.66E-20	miR-196b-5p	4.94	3.33E-28
miR-133b	−3.03	9.52E-20	miR-767-5p	4.05	3.28E-06
miR-802	−2.46	0.000768	miR-552-3p	3.55	0.000579
miR-490-3p	−2.46	2.63E-08	miR-105-5p	3.49	7.65E-06
miR-1265	−2.43	0.014022	miR-135b-5p	3.46	2.25E-21
miR-204-5p	−2.10	3.67E-11	miR-767-3p	3.44	0.003757
miR-145-5p	−1.87	3.04E-17	miR-194-5p	2.97	2.20E-22
miR-139-3p	−1.76	2.84E-21	miR-200a-5p	2.94	2.40E-23
miR-145-3p	−1.75	2.37E-19	miR-192-5p	2.87	1.71E-18
miR-129-5p	−1.70	4.41E-08	miR-200a-3p	2.81	2.26E-19
miR-139-5p	−1.59	2.49E-21	miR-200b-3p	2.80	1.08E-22
miR-30a-3p	−1.41	4.25E-13	miR-141-5p	2.79	5.05E-23
miR-490-5p	−1.40	6.09E-07	miR-1269a	2.78	0.003274
miR-551b-3p	−1.32	1.73E-06	miR-183-5p	2.71	4.27E-21
miR-143-3p	−1.27	1.05E-11	miR-194-3p	2.64	1.03E-15
miR-486-5p	−1.22	2.70E-07	miR-429	2.58	2.15E-14
miR-29c-3p	−1.14	7.86E-12	miR-141-3p	2.57	5.25E-17
miR-383-5p	−1.12	1.13E-05	miR-675-5p	2.55	0.011142
miR-195-3p	−1.06	1.83E-09	miR-146b-5p	2.50	4.41E-38

DOI: 10.7717/peerj.7882/table-2

miR-383-5p was confirmed to be downregulated in GC

The expression value of miR-383-5p was extracted from miRNA-Seq data. Totally, there were 387 GC samples and 41 normal samples. The expression of miR-383-5p in GC was significantly lower than normal tissues (Fig. 2A). The results of qRT-PCR showed that miR-383-5p was significantly down-regulated in GC at the level of cell and tissue (Figs. 2B and 2C). Combining the miRNA-Seq data and qRT-PCR assay, we confirmed that miR-383-5p was down-regulated and might be a novel tumor suppressor gene in GC.

Figure 2: The expression of miR-383-5p in GC tissues and cells. (A) miRNA-Seq data indicated that miR-383-5p was down-regulated in GC tissues; (B) and (C) qRT-PCR showed that miR-383-5p was significantly down-regulated in GC cells and tissues. *P [less than] 0.05 compared with the control group (GES-1 or normal tissues). DOI: 10.7717/peerj.7882/fig-2

Association between miR-383-5p and prognosis, clinicopathological features

A total of 382 GC TCGA samples with necessary data were selected to investigate the prognostic role of miR-383-5p. According to the median expression value of miR-383-5p, patients were equally divided to the low and high expression groups. The Kaplan–Meier survival analysis indicated that patients with low miR-383-5p expression had a shorter survival time (median survival time 21.1 months) than those with high expression (46.9 months) (Fig. 3). Furthermore, we explored the association between miR-383-5p expression and clinicopathological features. The results manifested that low miR-383-5p expression was significantly associated with large tumor size and poor differentiation grade (Table 3). Nevertheless, the features of age, gender, lymph node metastasis, and TNM stage were found to be of no significant difference.

Figure 3: The Kaplan–Meier overall survival curve. The group with low miR-383-5p expression had a significantly less survival time (median survival time 21.1 months) than that with high expression (46.9 months). DOI: 10.7717/peerj.7882/fig-3

[ Table omitted. See PDF. ]

DOI: 10.7717/peerj.7882/table-3

Note:

*P < 0.05, statistically significant.

Overexpression of miR-383-5p inhibits GC cell proliferation and migration

The CCK-8 and wound healing assays were performed to assess the effect of miR-383-5p on the proliferation and migration of GC cells. Compared with NC group, transfection with miR-383-5p mimics weakened the migration capacity of MGC-803 and MKN-45 cells (Figs. 4A–4H). Besides, the CCK-8 assay showed that miR-383-5p mimics inhibited GC cells proliferation (Figs. 4I and 4J).

Figure 4: The effect of miR-383-5p on the proliferation and migration of GC cells. (A–H) miR-383-5p mimics inhibited GC cells migration; (I) and (J) miR-383-5p mimics inhibited GC cells proliferation. *P [less than] 0.05 compared with NC group. DOI: 10.7717/peerj.7882/fig-4

Target prediction and genes functional enrichment analyses

Three bioinformatics websites (TargetScan, miRDB, and starBase) were selected for predicting target genes of miR-383-5p. In view of that each website had diverse bioinformatics algorithm, we took the consensus results of different predictions. As described in the venn plot, 49 consensus target genes were obtained (Fig. 5A). To comprehend the function of miR-383-5p target genes, the 49 genes were used for functional enrichment analysis by FunRich. At the aspect of biological pathway analysis, we found that these genes participated in PI3K, mTOR, c-MYC, TGF-beat receptor, VEGF/VEGFR and E-cadherin signaling pathways (Fig. 5B).

Figure 5: The venn plot of target genes and functional enrichment analysis. (A) A total of 49 consensus genes were obtained from TargetScan, miRDB, and starBase websites; (B) Biological pathway analysis revealed that these genes participated in PI3K, mTOR, c-MYC, TGF-beta receptor, VEGF/VEGFR and E-cadherin signaling pathways. DOI: 10.7717/peerj.7882/fig-5

miR-383-5p targets LDHA

LDHA is one of the 49 common target genes and it has a quite higher target prediction score. The miRNA-target prediction from starBase demonstrated that miR-383-5p had potential binding site on the 3′UTR of LDHA (Fig. 6A). The association regression analysis based on starBase project indicated that miR-383-5p expression was negatively correlated with LDHA (Fig. 6B). Furthermore, luciferase reporter assay suggested that the luciferase activity in WT LDHA-3′UTR group was significantly decreased by the miR-383-5p mimics, and there were no differences in the mutant LDHA-3′UTR group (Fig. 6C). All the results indicated that LDHA was a target of miR-383-5p.

Figure 6: miR-383-5p directly targets LDHA. (A) The predicted miR-383-5p binding sites on the 3′-UTR of LDHA; (B) The data based on starBase project showed that expression of miR-383-5p was negatively correlated with target LDHA (r = −0.203); (C) Luciferase reporter assay indicated that the luciferase activity in wild-type LDHA-3′UTR group was significantly decreased by the miR-383-5p mimics. *P [less than] 0.05. DOI: 10.7717/peerj.7882/fig-6

Discussion

Gastric cancer has the malignant features of terrible proliferation, invasion, metastasis, and multiple drug resistance, which lead to high mortality and poor prognosis. Increasing studies have proved that miRNAs are aberrantly expressed and involved in the initiation and development of GC (Chen et al., 2019; Kang et al., 2018; Maruyama et al., 2018; Wang et al., 2019a). Thus, identifying DEMs and exploring the biological function of miRNAs can be useful for finding novel biomarkers and understanding the mechanism of GC progression.

In this study, we first downloaded and analyzed the GC miRNA-Seq and clinical data from Firebrowse website (Deng et al., 2017), which conserved integrated gene expression profiles and clinical information data from TCGA. Through screening for DEMs, we found that miR-383-5p was down-regulated in GC tissues. Besides, we also investigated the prognostic role of miR-383-5p, and Kaplan–Meier survival analysis indicated that patients with low miR-383-5p expression had a shorter survival time than those with high expression. All these results inspired us that miR-383-5p may play an important role in GC.

Furthermore, we retrieved literatures published worldwide to comprehend the studies about miR-383-5p. Zhao et al. (2017) found that miR-383-5p was significantly decreased in lung adenocarcinoma and overexpression of miR-383-5p inhibited cell proliferation by G1 cell cycle phase arrest and induced apoptosis in vitro. In hepatocellular carcinoma, miR-383-5p was proved to be a tumor suppressor and to modulate hepatocellular carcinoma tumorigenesis and progress by targeting AKR1B10 (Wang et al., 2018) and LDHA (Fang et al., 2017). Besides, Jiang et al.’s (2019) study reported that overexpression of miR-383-5p could inhibit ovarian cancer cell proliferation and enhance chemosensitivity of cells by regulating TRIM27. MiR-383-5p could also suppress ovarian cancer cell proliferation, invasion and aerobic glycolysis through regulating LDHA (Han et al., 2017). Azarbarzin et al. (2017) found that miR-383-5p was downregulated in intestinal-type GC and could be used a diagnostic biomarker. However, the molecular function and clinical significance of miR-383-5p in GC have not been studied. The present study investigated that miR-383-5p was decreased in GC and its expression was associated with tumor size and differentiation grade. Furthermore, the CCK-8 and wound healing assays demonstrated that overexpression of miR-383-5p could inhibit GC cells proliferation and migration. All the results indicated that miR-383-5p could act as a tumor suppressor in GC and it also had vital clinical value.

We further explored the downstream regulation role by predicting the potential target genes of miR-383-5p. The consensus target genes were obtained by integrating the results from three bioinformatics tools, which improved the accuracy of prediction. The functional enrichment analysis demonstrated that miR-383-5p may be involved in PI3K, mTOR, c-MYC, TGF-beta receptor, VEGF/VEGFR, and E-cadherin signaling pathways through regulating the target genes. It is well known that the mTOR pathway regulates tumor growth and metastasis by mediating tumor metabolic homeostasis (Xia & Xu, 2015). Multiple miRNAs were reported to participate in the regulation of PI3K/AKT/mTOR signaling pathway (Riquelme et al., 2016). Yu’s study showed that miR-106b was overexpressed in CD44(+) GC stem-like cells and could retain cancer stem cell characteristics through modulating TGF-β/Smad signaling pathway (Yu et al., 2014). MiR-372 negatively targets KIF26B to suppresses GC cells proliferation and metastasis by regulating VEGF pathway (Zhang et al., 2017). Above all, miR-383-5p may act as a novel tumor suppressor in taking part in the biological function of GC.

LDHA acts as a glycolytic enzyme in the process of catalyzing the formation of lactate from pyruvate. LDHA is not only involved in normal cells metabolism, but closely related to tumor malignancy. Studies have found that LDHA was upregulated in kinds of malignant tumors and can play an oncogene role in GC (Wang et al., 2017; Zhu et al., 2018). The association regression analysis based on starBase project indicated that miR-383-5p expression was negatively correlated with LDHA. Besides, luciferase reporter assay showed that the luciferase activity in WT LDHA-3′UTR group was significantly decreased by the miR-383-5p mimics, and there were no differences in the mutant LDHA-3′UTR group. All the results suggested that LDHA was a downstream target of miR-383-5p.

Conclusions

In summary, we found that miR-383-5p may act as a tumor suppressor in GC. It is of important clinical significance and prognostic value, which could contribute to revealing the molecular mechanism of GC tumorigenesis and progress.

Additional Information and Declarations

Competing Interests

The authors declare that they have no competing interests.

Author Contributions

Chao Wei conceived and designed the experiments, performed the experiments, analyzed the data, contributed reagents/materials/analysis tools, prepared figures and/or tables, authored or reviewed drafts of the paper, approved the final draft.

Jian-Jun Gao authored or reviewed drafts of the paper, approved the final draft, design the study.

Human Ethics

The following information was supplied relating to ethical approvals (i.e., approving body and any reference numbers):

The Clinical Research Ethics Committee of Jinzhou Medical University approved this research (No: EC-2018-JZ-016).

Data Availability

The following information was supplied regarding data availability:

Data is available at Firehose, http://gdac.broadinstitute.org/: Stomach adenocarcinoma; Cohort: STAD.

Funding

The authors received no funding for this work.

References

Agarwal V, Bell GW, Nam J-W, Bartel DP. 2015. Predicting effective microRNA target sites in mammalian mRNAs. eLife 4:e05005

Azarbarzin S, Feizi MAH, Safaralizadeh R, Kazemzadeh M, Fateh A. 2017. The value of MiR-383, an intronic MiRNA, as a diagnostic and prognostic biomarker in intestinal-type gastric cancer. Biochemical Genetics 55(3):244-252

Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. 2018. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA: A Cancer Journal for Clinicians 68(6):394-424

Carcas LP. 2014. Gastric cancer review. Journal of Carcinogenesis 13(1):14

Chen Z, Li Z, Soutto M, Wang W, Piazuelo MB, Zhu S, Guo Y, Maturana MJ, Corvalan AH, Chen X, Xu Z, El-Rifai WM+2 more. 2019. Integrated analysis of mouse and human gastric neoplasms identifies conserved microRNA networks in gastric carcinogenesis. Gastroenterology 156(4):1127-1139.e8

Chen D, Wang H, Chen J, Li Z, Li S, Hu Z, Huang S, Zhao Y, He X. 2018. MicroRNA-129-5p regulates glycolysis and cell proliferation by targeting the glucose transporter SLC2A3 in gastric cancer cells. Frontiers in Pharmacology 9:502

Chen W, Zheng R, Baade PD, Zhang S, Zeng H, Bray F, Jemal A, Yu XQ, He J. 2016. Cancer statistics in China, 2015. CA: A Cancer Journal for Clinicians 66(2):115-132

Correa P. 2013. Gastric cancer: overview. Gastroenterology Clinics of North America 42(2):211-217

Deng M, Bragelmann J, Kryukov I, Saraiva-Agostinho N, Perner S. 2017. FirebrowseR: an R client to the Broad Institute’s Firehose pipeline. Database 2017:baw160

Fang Z, He L, Jia H, Huang Q, Chen D, Zhang Z. 2017. The miR-383-LDHA axis regulates cell proliferation, invasion and glycolysis in hepatocellular cancer. Iranian Journal of Basic Medical Sciences 20:187-192

Han RL, Wang FP, Zhang PA, Zhou XY, Li Y. 2017. miR-383 inhibits ovarian cancer cell proliferation, invasion and aerobic glycolysis by targeting LDHA. Neoplasma 64(2):244-252

Hui W, Ma X, Zan Y, Song L, Zhang S, Dong L. 2018. MicroRNA-1292-5p inhibits cell growth, migration and invasion of gastric carcinoma by targeting DEK. American Journal of Cancer Research 8:1228-1238

Jiang J, Xie C, Liu Y, Shi Q, Chen Y. 2019. Up-regulation of miR-383-5p suppresses proliferation and enhances chemosensitivity in ovarian cancer cells by targeting TRIM27. Biomedicine & Pharmacotherapy 109:595-601

Kang W, Huang T, Zhou Y, Zhang J, Lung RWM, Tong JHM, Chan AWH, Zhang B, Wong CC, Wu F, Dong Y, Wang S, Yang W, Pan Y, Chak WP, Cheung AHK, Pang JCS, Yu J, Cheng ASL, To KF+10 more. 2018. miR-375 is involved in Hippo pathway by targeting YAP1/TEAD4-CTGF axis in gastric carcinogenesis. Cell Death & Disease 9(2):92

Li J-H, Liu S, Zhou H, Qu L-H, Yang J-H. 2014. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein–RNA interaction networks from large-scale CLIP-seq data. Nucleic Acids Research 42(D1):D92-D97

Maruyama S, Furuya S, Shiraishi K, Shimizu H, Akaike H, Hosomura N, Kawaguchi Y, Amemiya H, Kawaida H, Sudo M, Inoue S, Kono H, Ichikawa D+3 more. 2018. miR-122-5p as a novel biomarker for alpha-fetoprotein-producing gastric cancer. World Journal of Gastrointestinal Oncology 10(10):344-350

Pathan M, Keerthikumar S, Ang CS, Gangoda L, Quek CY, Williamson NA, Mouradov D, Sieber OM, Simpson RJ, Salim A, Bacic A, Hill AF, Stroud DA, Ryan MT, Agbinya JI, Mariadason JM, Burgess AW, Mathivanan S+8 more. 2015. FunRich: an open access standalone functional enrichment and interaction network analysis tool. PROTEOMICS 15(15):2597-2601

Riquelme I, Tapia O, Leal P, Sandoval A, Varga MG, Letelier P, Buchegger K, Bizama C, Espinoza JA, Peek RM, Araya JC, Roa JC+2 more. 2016. miR-101-2, miR-125b-2 and miR-451a act as potential tumor suppressors in gastric cancer through regulation of the PI3K/AKT/mTOR pathway. Cell Oncology 39:23-33

Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. 2015. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Research 43(7):e47

Wang X. 2016. Improving microRNA target prediction by modeling with unambiguously identified microRNA-target pairs from CLIP-ligation studies. Bioinformatics 32(9):1316-1322

Wang R, Sun Y, Yu W, Yan Y, Qiao M, Jiang R, Guan W, Wang L. 2019a. Downregulation of miRNA-214 in cancer-associated fibroblasts contributes to migration and invasion of gastric cancer cells through targeting FGF9 and inducing EMT. Journal of Experimental & Clinical Cancer Research 38(1):20

Wang Y-N, Xu F, Zhang P, Wang P, Wei Y-N, Wu C, Cheng S-J. 2019b. MicroRNA-575 regulates development of gastric cancer by targeting PTEN. Biomedicine & Pharmacotherapy 113:108716

Wang J, Zhou Y, Fei X, Chen X, Chen Y. 2018. Biostatistics mining associated method identifies AKR1B10 enhancing hepatocellular carcinoma cell growth and degenerated by miR-383-5p. Scientific Reports 8(1):11094

Wang H, Zhou R, Sun L, Xia J, Yang X, Pan C, Huang N, Shi M, Bin J, Liao Y, Liao W+1 more. 2017. TOP1MT deficiency promotes GC invasion and migration via the enhancements of LDHA expression and aerobic glycolysis. Endocrine-Related Cancer 24(11):565-578

Xia P, Xu XY. 2015. PI3K/Akt/mTOR signaling pathway in cancer stem cells: from basic research to clinical application. American Journal of Cancer Research 5:1602-1609

Yu D, Shin H-S, Lee YS, Lee YC. 2014. miR-106b modulates cancer stem cell characteristics through TGF-β/smad signaling in CD44-positive gastric cancer cells. Laboratory Investigation 94(12):1370-1381

Zen K, Zhang C-Y. 2012. Circulating microRNAs: a novel class of biomarkers to diagnose and monitor human cancers. Medicinal Research Reviews 32(2):326-348

Zhang H, Ma R-R, Wang X-J, Su Z-X, Chen X, Shi D-B, Guo X-Y, Liu H-T, Gao P. 2017. KIF26B, a novel oncogene, promotes proliferation and metastasis by activating the VEGF pathway in gastric cancer. Oncogene 36(40):5609-5619

Zhang Y, Wang Z, Gemeinhart RA. 2013. Progress in microRNA delivery. Journal of Controlled Release 172(3):962-974

Zhao S, Gao X, Zang S, Li Y, Feng X, Yuan X. 2017. MicroRNA-383-5p acts as a prognostic marker and inhibitor of cell proliferation in lung adenocarcinoma by cancerous inhibitor of protein phosphatase 2A. Oncology Letters 14(3):3573-3579

Zhu W, Ma L, Qian J, Xu J, Xu T, Pang L, Zhou H, Shu Y, Zhou J. 2018. The molecular mechanism and clinical significance of LDHA in HER2-mediated progression of gastric cancer. American Journal of Translational Research 10:2055-2067