MiR-96-5p inhibition induces cell apoptosis in gastric adenocarcinoma

He-Ying Zhou, Chun-Qi Wu, En-Xu Bi

Abstract BACKGROUND Gastric adenocarcinoma (GAC) mortality rates have remained relatively changed over the past 30 years, and it continues to be one of the leading causes of cancerrelated death.AIM To search for novel miRNAs related to GAC prognosis and further investigate the effect of miR-96-5p on MGC-803 cells.METHODS The miRNA expression profile data of GAC based on The Cancer Genome Atlas were obtained and used to screen differently expressed miRNAs (DEMs) and DEMs related to GAC prognosis. Then, the expression of DEMs related to GAC prognosis was identified in GAC tumor samples and adjacent normal samples by qRT-PCR. The target gene, ZDHHC5, of miR-96-5p was predicted using TargetScan, miRTarBase, and miRDB databases and confirmed by luciferase reporter assay. Furthermore, MGC-803 cells were transfected with inhibitor NC,miR-96-5p inhibitor, si-ZDHHC5, or miR-96-5p inhibitor + si-ZDHHC5, and then cell apoptosis was detected by flow cytometry. The expression of ZDHHC5, Bcl-2, and COX-2 was detected using western blotting.RESULTS A total of 299 DEMs and 35 DEMs related to GAC prognosis were screened based on The Cancer Genome Atlas. Then compared with adjacent normal samples, the levels of miR-96-5p, miR-222-5p, and miR-652-5p were remarkably increased,while miR-125-5p, miR-145-3p, and miR-379-3p levels were reduced in GAC tumor samples (P ＜ 0.01), which were consistent with bioinformatics analysis.Furthermore, ZDHHC5 was defined as a direct target gene of miR-96-5p. miR-96-5p inhibition increased the number of apoptotic cells as well as promoted the expression of ZDHHC5, Bcl-2, and COX-2 in MGC-803 cells (P ＜ 0.01). After ZDHHC5 inhibition, the number of apoptotic cells and the expression of ZDHHC5, Bcl-2, and COX-2 were reduced. The addition of an miR-96-5p inhibitor partly reversed these effects (P ＜ 0.01).CONCLUSION Our findings identified six miRNAs related to GAC prognosis and suggested that downregulated miR-96-5p might induce cell apoptosis via upregulating ZDHHC5 expression in MGC-803 cells.

Key words: Gastric adenocarcinoma; Differently expressed miRNAs; Prognosis;MicroRNA-96-5p; Cell apoptosis; The Cancer Genome Atlas

INTRODUCTION

Gastric adenocarcinoma (GAC) is the most common malignant tumor originating in the stomach and is counted as one of the top ten common cancers worldwide, with approximately 951000 diagnosed cases and 723000 deaths in 2012[1,2]. Currently, the common and effective therapeutic method is the combination of surgery and adjuvant radiation therapy or chemotherapy, which have improved the 5-year survival rate of GAC[3]. However, delayed diagnosis occurs in most patients with proximal or distant metastasis due to the nontypical symptoms of early GAC, which results in poor treatment and prognosis[4]. Therefore, it is important to further reveal novel diagnostic and therapeutic methods as well as the underlying molecular mechanism of GAC.

It is widely known that a major challenge for the treatment of GAC is poor prognosis, and environmental exposure and gene mutation have been identified to be associated with this outcome[5]. Plenty of evidence indicates that the poor prognosis of GAC is significantly related to many molecular biomarkers, such as microRNAs(miRNA)[6,7]. miRNAs, as endogenous noncoding small-molecule RNAs, widely exist in severe conditions[8]. It is known to function in post-transcriptional regulation of gene expression through binding 3'-untranslated region of their target mRNA,accordingly modulating various key cell biological processes, such as embryonic development, tumor cell proliferation, differentiation, and apoptosis[8,9].

Previous studies have demonstrated that miRNA dysregulation significantly influences the prognosis of gastric cancer patients (e.g., miRNA-203[10], miR-21[11], and miR-25[12]. Imaokaet al[10]reported that a low serum miR-203 expression is associated with poor prognosis and may be a noninvasive biomarker for prognosis of gastric cancer patients. Simonianet al[11]observed that circulating miR-21 may be considered as a diagnostic and prognostic biomarker in gastric cancer. In addition, Liet al[12]revealed that miR-25 is associated with the prognosis of gastric cancer and can induce cell migration and proliferation by targeting transducer of ERBB2.1. Thus, it is essential to search for more novel miRNAs related to GAC prognosis, which may contribute to the development of GAC diagnosis.

In the current study, the miRNA expression profile data of GAC based on The Cancer Genome Atlas (TCGA) were analyzed to screen differently expressed miRNAs(DEMs) and DEMs related to GAC prognosis. Furthermore, DEMs were identified in clinical samples, and the mechanism of DEM was investigatedin vitro. According to this, we aimed to search for new therapeutic targets for GAC and provide some useful insights in improving the prognosis of GAC patients.

MATERIALS AND METHODS

Data extraction and DEM screening

The miRNA expression profile data (level 3, processed and standardized data) and the corresponding clinical information of GAC were downloaded from TCGA(https://portal.gdc.cancer.gov/) on February 11, 2019 based on the platform of Illumina HiSeq 2000 RNA Sequencing platform. A total of 452 samples were obtained from this dataset, including 410 GAC tumor samples and 42 normal control samples.The edgeR package in R was utilized to screen DEMs between GAC samples and normal samples. The thresholds were defined as false discovery rate ＜ 0.05 and |log fold change| ＞ 1. Meanwhile, volcano plots and heat maps were generated based on the obtained DEMs.

DEMs screening related to prognosis

The overall survival time was individually extracted from clinical information. Then,combined with the overall survival times and the expression levels of DEMs, DEMs related to prognosis were screened using KMsurv package of R, with the threshold of log-rankP＜ 0.05.

Clinical validation sample collection

This study obtained ethical approval from the ethics committee of Jinan Seventh People's Hospital, and the study was performed according to the Helsinki Declaration. A total of 20 paired tumor tissues and adjacent normal tissues (distance of 3-4 cm from the tumor tissue) were collected from GAC patients who underwent surgery in Jinan Seventh People's Hospital between September 2018 to September 2019. The specimens were confirmed by hematoxylin eosin staining and stored in RNA later. In addition, 5 mL peripheral blood was obtained from these 20 GAC patients. Meanwhile, the same amount of peripheral blood was extracted from 20 paired healthy subjects. Written informed consent from all participants was obtained,and the clinical information, including age, weight, gender, distant metastasis, lymph node metastasis, depth of invasion, and TNM stage are shown in Table 1.

Table 1 Distribution of characteristics in gastric adenocarcinoma patients and healthy subjects

Predicting target genes of DEMs

Target genes of DEMs related to prognosis were predicted using the three online analysis databases, including miRDB, miRTarBase, and TargetScan. Overlapping target genes among the three tools were selected to make the bioinformatic analysis more reliable. Then, the intersection of the predicting target genes among the three databases was obtained using a Venn diagram online tool, and the target genes that overlapped in the three databases were considered as a potential target gene of DEM.

Cell culture and transfection

Human gastric carcinoma cell line MGC-803 was obtained from Shanghai Obio Technology Co., Ltd. The cells were maintained in Dulbecco's Modified Eagle Media(DMEM, Gibco, Carlsbad, CA, United States) supplemented with 10% fetal bovine serum (FBS, Gibco, Carlsbad, CA, United States). The construction of ZDHHC5 silence vector (si-ZDHHC5) was performed by GenePharma (Shanghai, China). The miR-96-5p inhibitor and inhibitor NC were purchased from Thermo (Waltham, MA,United States). MGC-803 cells were inoculated in six-well plates for 24 h with approximately 5 × 105cells in each well, and then inhibitor NC, miR-96-5p inhibitor,si-ZDHHC5, or miR-96-5p inhibitor + si-ZDHHC5 was transfected into MGC-803 cells by Lipofectamine 2000 (Invitrogen, Carlsbad, CA, United States) per the manufacturer's instructions. Meanwhile, MGC-803 cells without transfection served as the control group. Cells were harvested after 48 h of transfection to perform follow-up experiments.

Real-time quantitative polymerase chain reaction

Total RNA from tissues and peripheral blood was obtained by Trizol (Invitrogen) and RNeasy Plus Mini Kit (Qiagen, Valencia, CA, United States) per the manufacturer's instructions. The concentration of RNA was detected using NanoDrop ND-2000(Invitrogen). To generate cDNA, Mir-X? miRNA FirstStrand Synthesis Kit (Takara,Dalian, China) was used. The quantitative polymerase chain reaction (qPCR) was carried out using the SYBR Premix ExTaqTM II (Takara) by ABI 7900 qRT-PCR System (Applied Biosystems, Foster City, CA, United States). The primer sequences are listed in Table 2. U6 and glyceraldehyde-phosphate dehydrogenase were used as the internal control of measuring miR-19a andADIPOR2expression. Data were analyzed by the 2-ΔΔCtmethod.

Table 2 Primers used for the quantitative real-time polymerase chain reaction

Luciferase reporter assay

The target gene of miR-96-5p was verified using the luciferase assay. The 3'-untranslated region of ZDHHC5 was cloned into a pGL3-basic vector, named as Luc-ZDHHC5. Luc-ZDHHC5 and phRL-TK plasmid were co-transfected with miR-96-5p mimic, miR-96-5p NC (negative control), or siZDHHC5 (positive control) (synthesized by Biosyntech, Suzhou, China) into 293T cells. After 48 h of transfection, the relative luciferase activity was measured by the Dual-Glo Luciferase Assay System (Promega)in accordance to the manufacturer's introductions. Renilla luciferase activity was used to normalize luciferase activity.

Western blotting

Total proteins were isolated by RIPA Lysis Buffer (Beyotime, Shanghai, China).Proteins concentrations were tested by bicinchoninic acid kit (Beyotime). The protein sample was separated on SDS-PAGE gel, transferred to polyvinylidene fluoride membranes, and followed by the blockage with 5% nonfat milk for 1 h. Next, the membranes were probed with primary antibodies of ZDHHC5 (1:1000, Proteintech),Bcl-2 (1:1000, Abcam), COX-2 (1:1000, Abcam), and GAPDH (1:1000, Beyotime)overnight at 4 °C. Then, membranes were incubated with secondary antibody (1:1000,Beyotime) for 2 h in a dark room at room temperature. GAPDH was used as the control protein. Enhanced chemiluminescence Plus reagent (Beyotime) was used to image blots. The band quantification was performed using Image J software.

Flow cytometry assay

Annexin V-FITC Apoptosis Detection kit was used to evaluate cell apoptosis. MGC-803 cells were grown in 6-well plates for 24 h and then transfected with inhibitor NC,miR-96-5p inhibitor, si-ZDHHC5, or miR-96-5p inhibitor + si-ZDHHC5 for 48 h. Next,cells were digested with trypsin and washed with PBS, followed by resuspending in 1× Binding Buffer, and stained with propidium iodide and FITC-Annexin V for 15 min at 25 °C in the dark. Cells were finally detected using a flow cytometer (Beckman Coulter, Fullerton, CA, United States).

Statistical analysis

Statistical analysis was conducted using SPSS Statistics software 22.0 (Chicago, IL,United States). Continuous variables were expressed as mean ± standard deviation and analyzed by independent-samplesttest. Categorical variables were expressed as percentages and assessed by two-sided chi-square test. The differences of multiple groups were performed by one-way ANOVA following with post-hoc of Dunnettttest.P＜ 0.05 was considered to be statistically significant.

RESULTS

DEMs between GAC sample and normal sample based on TCGA

Based on the selective criteria, a total of 299 DEMs were identified between GAC and normal control samples, including 225 upregulated and 74 downregulated miRNAs.As shown in Figure 1A and 1B, volcano plots and heat maps were conducted for these 299 DEMs.

DEMs related to prognosis based on TCGA

Based on these 299 DEMs, the relationships between patient overall survival and miRNA expression were evaluated, and the results showed that 35 DEMs were significantly related to the prognosis of GAC patients (P＜ 0.05). Among these DEMs,seven miRNAs had a higher association with GAC prognosis (P＜ 0.01), including miR-96-5p (P= 8.049 × 10-3), miR-125-5p (P= 9.638 × 10-4), miR-145-3p (P= 6.002 ×10-3), miR-222-5p (P= 1.812 × 10-3), miR-379-3p (P= 5.032 × 10-3), miR-652-5p (P= 3.145× 10-3), and miR-708-3p (P= 7.984 × 10-3) (Figure 2).

Figure 2 Association of differentially expressed miRNAs with overall survival of gastric adenocarcinoma. Dark lines represent high expression, and lighter lines represent low expression.

DEMs identification in clinical samples

A total of 20 GAC patients and 20 healthy subjects were included in this study. No significant difference was found in age, weight, and gender between GAC patients and healthy subjects (Table 1). Based on the above survival analysis, six miRNAs were selected for identification in GAC tumor samples and adjacent normal samples. qRTPCR revealed that compared with adjacent normal samples, the levels of miR-96-5p,miR-222-5p, and miR-652-5p were remarkably increased, while miR-125-5p, miR-145-3p, and miR-379-3p levels were obviously reduced in GAC tumor samples (P＜ 0.01,Figure 3A), which was consistent with bioinformatics analysis results by TCGA.Moreover, miR-96-5p level was detected in the blood of GAC patients and healthy subjects, but no significant difference was found (Figure 3B).

Figure 3 The identification of differentially expressed miRNA levels in gastric adenocarcinoma patients and healthy subjects. A: The levels of miR-96-5p,miR-222-5p, miR-652-5p, miR-125-5p, miR-145-3p, and miR-379-3p in tumor samples and adjacent normal samples by qRT-PCR; B: The miR-96-5p level in the blood of gastric adenocarcinoma patients and healthy subjects by qRT-PCR (bP ＜ 0.01). Values are mean ± SD.

Target gene prediction and identification of miR-96-5p

Considering miR-96-5p had the highest association with GAC prognosis, the function of miR-96-5p was investigated in the following experiments. The results found that a total of 39 overlapped target genes existed in TargetScan, miRTarBase, and miRDB databases (Figure 4A). Based on this bioinformatics analysis,ZDHHC5was considered as a potential target gene of miR-96-5p (Figure 4B). Luciferase receptor assay showed that the relative luciferase activity was reduced after co-transfection with miR-96-5p mimic or siZDHHC5 compared with co-transfection with miR-96-5p NC (Figure 4C), which suggestedZDHHC5was a direct target gene of miR-96-5p.

Figure 4 Target gene prediction and identification of miR-96-5p. A: The intersection of the predicting target genes of miR-96-5p among the TargetScan,miRTarBase and miRDB databases using a Venn diagram; B: Prediction of the binding site of miR-96-5p to zinc finger DHHC-type palmitoyltransferase 5 (ZDHHC5)by bioinformatics; C: Target regulation of miR-96-5p to ZDHHC5 observed by the luciferase reporter system. ZDHHC5: Zinc finger DHHC-type palmitoyltransferase 5.

Effect of miR-96-5p on apoptosis in MGC-803 cells

To further investigate the effects of miR-96-5p on GAC, the miR-96-5p inhibitor was used to inhibit miR-96-5p in MGC-803 cells. Flow cytometry assay showed that the number of apoptotic cells increased in MGC-803 cells transfected with the miR-96-5p inhibitor, while inhibition ofZDHHC5decreased cell apoptosis compared with cells with NC. Notably, co-transfection of the miR-96-5p inhibitor and si-ZDHHC5 partly reversed the effect of inhibitingZDHHC5on cell apoptosis (P＜ 0.01, Figure 5A). In addition, western blotting revealed that compared with MGC-803 cells without treatment, miR-96-5p inhibition promoted the expression of ZDHHC5, Bcl-2, and COX-2 (apoptosis proteins) in MGC-803 cells (P＜ 0.01, Figure 5B). However,inhibitingZDHHC5decreased the expression of ZDHHC5, Bcl-2, and COX-2. The addition of the miR-96-5p inhibitor increased the expression of these three proteins (P＜ 0.01, Figure 5B).

Figure 5 miR-96-5p inhibition might induce cell apoptosis in MGC-803 cells. A: The number of apoptotic cells after transfection with inhibitor NC, miR-96-5pinhibitor, si-zinc finger DHHC-type palmitoyltransferase 5 (ZDHHC5), or miR-96-5p inhibitor + si-ZDHHC5 in MGC-803 cells by flow cytometry analysis; B: Theexpression of ZDHHC5, Bcl-2, and COX2 after cells were treated with inhibitor NC, miR-96-5p inhibitor, si-ZDHHC5, or miR-96-5p inhibitor + si-ZDHHC5 in MGC-803 cells by western blotting (aP ＜ 0.01). Values are mean ± SD. ZDHHC5: Zinc finger DHHC-type palmitoyltransferase 5.

DISCUSSION

In the present study, a total of 299 DEMs and 35 DEMs related to GAC prognosis were screened based on the miRNA expression profile data from TCGA. Then, six miRNAs were selected for identification in GAC tumor samples and adjacent normal samples.The results were consistent with bioinformatics analysis. Furthermore, miR-96-5p was considered as an important biomarker and investigated in thein vitroexperiments.Our results revealed thatZDHHC5was a direct target gene of miR-96-5p, and miR-96-5p inhibition increased the expression of Bcl-2 and COX-2.

Six miRNAs were identified in this study, and the results showed that the levels of miR-96-5p, miR-222-5p, and miR-652-5p were overexpressed, while miR-125-5p, miR-145-3p, and miR-379-3p levels were downregulated in GAC sample. Several studies have demonstrated that miR-96-5p is overexpressed in various cancers, including colorectal cancer[13], pancreatic carcinoma[14], prostate cancer[15], hepatocellular carcinoma[16], and breast cancer[17], and it is an oncogene by promoting cell proliferation. miR-652-5p was reported to be associated with non-small cell lung cancer[18], esophageal adenocarcinoma[19], and breast cancer[20], while the mechanism of miR-652-5p was unknown. Current studies of miR-222-5p are focused on the role of angiogenesis in endothelium[21,22], and few studies investigated the effect of miR-222-5p in cancers. miR-125-5p was identified as a tumor suppressor in glioblastoma[23],cervical cancer[24], and renal cell carcinoma[25], and it was involved in proliferation,migration, and apoptosis. Many have investigated the role of miR-145-3p in cancers,such as bladder cancer[26], lung squamous cell carcinoma[27], gallbladder cancer[28], and head and neck squamous cell carcinoma[29]. It is also considered a tumor suppressor.Few studies have investigated the role of miRNA-379-3p; only one recent study reported that miRNA-379-5p exerted an antitumor effect by regulating tumor invasion and metastasis in hepatocellular carcinoma[30]. Unfortunately, the effects of these miRNAs on GAC have not been reported. Therefore, it is essential to further reveal the mechanism and prognostic significance of these miRNAs in GAC.

Due to the highest association of miR-96-5p with GAC prognosis, the effects of miR-96-5p on MGC-803 cells were investigated in this study. Notably, a previous study has shown that miR-96-5p exerts an inhibiting role in cell proliferation and migration by downregulation of FoxQ1 in gastric cancer cells[31]. Contradictorily, a recent study has demonstrated that miR-96-5p exerts a promoting effect on cell progression by directly targeting FOXO3 in gastric cancer. Consistent with this recent study[32], this study found that the miR-96-5p inhibitor induced cell apoptosis in MGC-803 cells.

It is generally acknowledged that miRNAs develop biological functions by impeding translation of target mRNAs. In agreement with the bioinformatics prediction, our study revealed thatZDHHC5was identified as a target gene of miR-96-5p.ZDHHC5, encoding zinc finger DHHC-type containing 5, is one member of the family of ZDHHC proteins and was identified as a putative palmitoyl S-acyltransferases[33]. It has been suggested that S-palmitoylation is closely associated with cancer development, and ZDHHC enzymes are the key enzymes responsible for palmitoylation[34].

Individual ZDHHC enzymes exert different effects on various cancers, either as tumor suppressors or oncoproteins[34]. A previous study documented that high expression of ZDHHC5 is associated with a poor prognosis in glioma[35]. In addition,the report of Tianet al[36]has suggested that DHHC5 knockdown can dramatically inhibit cell proliferation and invasion in non-small cell lung cancer. The present study revealed that miR-96-5p inhibition increased the number of apoptotic cells as well as promoted the expression of ZDHHC5, Bcl-2, and COX-2 in MGC-803 cells, while inhibitingZDHHC5decreased cell apoptosis. Co-transfection of the miR-96-5p inhibitor and si-ZDHHC5 partly reversed the effect ofZDHHC5inhibition on cell apoptosis. These results indicated that miR-96-5p inhibition induced cell apoptosis by upregulating ZDHHC5 expression. This result was inconsistent with previous studies,which may be due to ZDHHC5 having different functions in different cancer types[37].

Figure 4

ARTICLE HIGHLIGHTS

molecular biomarkers, such as microRNA (miRNA). Previous studies have demonstrated that miRNA dysregulation significantly influences the prognosis of gastric cancer patients (e.g.,miRNA-203, miR-21, and miR-25). Thus, it is essential to search for novel miRNAs related to GAC prognosis, which may contribute to the development of GAC diagnosis.

Research objectives

This study aimed to search for new miRNA therapeutic targets for GAC and investigate the mechanism of differentially expressed miRNA (DEM)in vitro, which might provide some useful insights in improving the prognosis of GAC patients.

Research methods

First, the miRNA expression profile data of GAC based on The Cancer Genome Atlas were obtained and used to screen DEMs and DEMs related to GAC prognosis by bioinformatics methods. Then, the expression of DEMs related to GAC prognosis was identified in GAC tumor samples and adjacent normal samples by qRT-PCR.ZDHHC5, a target gene of miR-96-5p, was predicted and confirmed by the luciferase reporter assay. Furthermore, MGC-803 cells were transfected with inhibitor NC, miR-96-5p inhibitor, si-ZDHHC5, or miR-96-5p inhibitor + si-ZDHHC5. Cell apoptosis was detected by flow cytometry. The expression of ZDHHC5, Bcl-2,and COX-2 was detected using western blotting.

Research results

A total of 299 DEMs and 35 DEMs related to GAC prognosis were screened based on the miRNA expression profile data from The Cancer Genome Atlas. Six miRNAs, including miR-96-5p, miR-222-5p, miR-652-5p, miR-125-5p, miR-145-3p, and miR-379-3p, were selected for identification in GAC tumor samples and adjacent normal samples. The results were consistent with bioinformatics analysis. Furthermore, miR-96-5p was considered as an important biomarker and investigated inin vitroexperiments. Our results revealed thatZDHHC5was a direct target gene of miR-96-5p, and miR-96-5p inhibition induced cell apoptosis and increased the expression of Bcl-2 and COX-2.

Research conclusions

In conclusion, this work identified six miRNAs related to GAC prognosis, including miR-96-5p,miR-125-5p, miR-145-3p, miR-222-5p, miR-379-3p, and miR-652-5p. Furthermore, downregulated miR-96-5p markedly induced cell apoptosis through targetingZDHHC5.

Research perspectives

Current findings provide a potential molecular mechanism of miR-96-5p in GAC. However,further studies are needed to investigate the mechanism and prognostic significance of these miRNAs in GAC.