Skip to main content

Identification of invasion-metastasis associated MiRNAs in gallbladder cancer by bioinformatics and experimental validation

Abstract

Background

Recent studies exploring the roles of invasion-metastasis associated miRNAs in gallbladder cancer (GBC) are limited. In the study, we aimed to identify the invasion-metastasis associated miRNAs in GBC by bioinformatics and experimental validation.

Methods

MiRNAs of different expression were identified by comparing GBC tumor samples with different survival from Gene Expression Omnibus database. MiRTarBase was used for identifying the potential target genes of miRNAs. Then, we performed Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. And miRNA-gene and protein–protein interaction (PPI) network were constructed for hub genes evaluation. We further explored and compared miR-642a-3p and miR-145-5p expression in both The Cancer Genome Atlas database and our hospital data. Finally, quantitative real-time PCR, wound healing assay, and Transwell assay were conducted to validate the invasion-metastasis associated miRNAs in GBC.

Results

In GSE104165 database, 25 up-regulated and 97 down-regulated miRNAs were detected with significantly different expression in GBC tumor samples. Then, 477 potential target genes were identified from the 2 most up-regulated miRNAs (miR-4430 and miR-642a-3p) and 268 genes from the 2 most down-regulated miRNAs (miR-451a and miR-145-5p). After GO and KEGG analysis, mTOR and PI3K-Akt signaling pathways were found associated with the potential target genes. Based on PPI network, the top 10 highest degree hub nodes were selected for hub genes. Furthermore, the miRNA-hub gene network showed significant miR-642a-3p up-regulation and miR-145-5p down-regulation in both GBC tissues and cell lines. In the experimental validation, miR-145-5p up-regulation and miR-642a-3p down-regulation were confirmed to suppress GBC invasion and metastasis.

Conclusions

MiR-642a-3p and miR-145-5p were identified as invasion-metastasis associated miRNAs via bioinformatics and experimental validation, and both up-regulation of miR-642a-3p and down-regulation of miR-145-5p would be served as novel treatment options for GBC in the future.

Background

Gallbladder cancer (GBC) is the most common type of biliary tract cancer (BTC), which has a high mortality and a poor prognosis [1,2,3,4]. According to the 8th American Joint Committee on Cancer (AJCC) guideline [5], surgical resection is the best potential treatment for GBC at an early stage, while chemotherapy/radiotherapy/immunotherapy/targeted therapy is recommended for GBC at an advanced stage. Patients with advanced GBC have limited response to treatment and poor prognosis due to the highly invasive and metastatic characteristics of GBC, including local tumor growth, hepatic invasion, and lymph nodes metastasis [6, 7]. Several studies reported the crucial role of multiple biological processes in cancer invasion and metastasis, such as circulating cancer cells, immune evasion [8], epithelial-mesenchymal transition (EMT) [9, 10], and cancer stem cells [11], while the multi-step process and the molecular mechanism of GBC remained unclear. Therefore, it is urgent to explore the novel biomarkers associated with invasion and metastasis in GBC patients to improve the prognosis.

MiRNAs, the small single-stranded non-coding RNAs (ranging from 18 to 25 bp in length), are critical in tumor development, such as invasion, migration, and metastasis [12, 13]. Different expression of miRNAs has been detected in various cancers. For example, miR-15a/16-1 could attenuate immunosuppression and represent a potential immunotherapy against hepatocellular carcinoma [14], meanwhile, miR-124-3p could promote cell proliferation in intrahepatic cholangiocarcinoma [15]. Additionally, by modulating STAT3 signaling, miR-124-3p was identified to suppress programmed cell death-ligand 1 (PD-L1) expression and inhibit tumorigenesis of colorectal cancer [16]. Nevertheless, as a highly aggressive cancer, studies exploring the potential roles of invasion-metastasis associated miRNAs in GBC are limited.

In the study, we aimed to conduct bioinformatic analysis to identify the invasion-metastasis associated miRNAs in GBC, and experimental validation were further performed to confirm the role of miRNAs. In turn, the study would provide potential effective treatment for improving the prognosis of GBC patients.

Methods

Data collection

Datasets of miRNA expression of GBC tissues with survival data from the National Center for Biotechnology Information (NCBI) (https://www.ncbi.nlm.nih.gov/geo) online Gene Expression Omnibus (GEO) database were included in the study. The search strategy was “gallbladder cancer” [MeSH Terms] AND “miRNA” [MeSH Terms]. After screening the datasets, the GSE104165 dataset, which contained 40 human GBC tumor samples and 8 normal gallbladder samples with survival data, was selected for further analysis. The sample platform was GPL18402 Agilent-046064 Unrestricted_Human_miRNA_V19.0_Microarray.

Screening miRNAs of different expression

The miRNA expression data profile in txt format were preprocessed. Limma package was applied to calculate the miRNA of different expression from the Bioconductor package (http://www.bioconductor.org/). MiRNAs of different expression were found by comparing GBC tumor samples with long survival and short survival in R Studio (Version 4.0.4). Different expression was determined as significant when P-value < 0.05 and |fold change (FC)|> 1.

Target genes prediction and gene ontology (GO) analysis

MiRTarBase (http://mirtarbase.mbc.nctu.edu.tw/php/index.php), an online miRNA-target interactions database, was used for identifying the potential target genes of the up-regulated or down-regulated miRNAs. The data were further validated by microarray, western blot, reporter assay, and next-generation sequencing experiments. GO analysis was extensively used for the determination of unique biological attributes of genes, gene products, and sequences, such as molecular function, cell components, and biological process [17]. Enrichr (http://amp.pharm.mssm.edu/Enrichr/), a 2016 updated webserver with a comprehensive enrichment analysis of the gene set, was introduced to manipulate enrichment analysis of the predicted target genes pathway and functional annotation, which was further analyzed for GO and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways [18].

MiRNA-gene and protein–protein interaction (PPI) network construction

Tools were used to determine the interactions between genes or proteins to establish the miRNA-gene and PPI network. The target genes were first entered to assess functional associations in the STRING database (http://string-db.org) [19]. Next, the connectivity degree of PPI network was calculated in Cytoscape software (Version 3.6.0) for further hub genes evaluation. Finally, a miRNA-hub gene regulatory network was generated.

The cancer genome atlas (TCGA) database and survival data of GBC patients

The BTC miRNA-sequencing data and clinical information were retrieved from TCGA database (https://portal.gdc.cancer.gov) to explore miR-642a-3p and miR-145-5p expression levels in BTC tissues. The data of overall survival (OS) were matched with miR-642a-3p and miR-145-5p expression values. A total of 41 GBC specimens were randomly selected from curative resection GBC patients retrospectively in Sir Run-Run Shaw Hospital (SRRSH), Hangzhou, Zhejiang Province, China. The interval between surgery and death was defined as OS. The follow-up period was defined as the duration between surgery and recurrence or death. The final follow-up was in December 2019.

Validation of potential target genes using TCGA database

Based on short and long survival of GBC patients in TCGA database, the related mRNA expression (SYK, SH3GL1, CDKN1A, MYC, VEGFA, and EGFR) were compared between two groups. Moreover, the expression of some potential mRNA (CDKN1A, MYC, VEGFA) was further compared based on different T stages, which represented tumor invasion.

Cell culture and tissue collection

Human GBC cell line GBC-SD and human intrahepatic biliary epithelial cells (HiBEC) were purchased from a cell bank (Chinese Academy of Sciences). The SGC-996 cell line was provided by Prof. Anonymize (Blinded Per Author Guidelines). The GBC-SD and SGC-996 cells were given 25 units/ml of penicillin, 25 g/ml of streptomycin, and 10% fetal bovine serum (FBS) (v/v) in a 5% (v/v) CO2 humidified incubator at 37 °C, and they were cultivated in DMEM and RPMI-1640, respectively. Our hospital provided the GBC tumor tissues and matched noncancerous tissues.

RNA extraction and quantitative real-time PCR (qRT-PCR)

TRIzol reagent (Ambion, USA) was used to extract total RNA from cells, and miRNeasy FFPE kits (Qiagen, USA) were used to isolate the FFPE tissues based on the manufacturer’s protocols. All-in-One™ miRNA qRT-PCR Detection Kit (GeneCopoeia, USA) or Hifair® II 1st Strand cDNA Synthesis SuperMix (Yeasen, China) were used to make complementary DNA from 1 μg of RNA. qRT-PCR was performed using a LightCycler® 480 (Roche Molecular Systems, Inc., USA). The primers used for miRNA qRT-PCR were shown in Additional file 1: Table S1, and 5s small nuclear RNA was the miRNA assays internal control. Based on the cycle thresholds (CT), qRT-PCR results were examined and calculated as the relative RNA levels. The FC in RNA levels between each sample was examined by the 2−ΔΔCT method.

Cell transfection

Hsa-miRNA mimics and cognate negative control RNAs were obtained from RiBo (Guangzhou, China) and transfected into GBC cell lines at the concentration of 50 nM using Lipofectamine™ 3000 (Invitrogen, USA) based on the user instructions.

Wound healing assay and transwell assay

After transfection, 5 × 104 GBC cells were seeded in ibidi Culture-Insert (ibidi GmbH, Martinsried, Germany) on a 12-well plate. The Culture-Insert was gently removed after appropriate cell attachment. Then the cells could migrate into the wound area. Photomicrographs were taken with a microscope at 0 h and 48 h after wounding, respectively. The cell invasion assay was performed by 24-well Transwell chambers (Corning, USA). Next, 1 × 105 cells were evenly suspended in a 200 µl serum-free medium and added to the up-layer inserts after transfection. 600 µl of medium with 20% FBS was then added to the low-layer compartment for chemoattractant. The cells on the membrane top surface were excised using a cotton-wool bud, whereas the cells on the bottom surface were fixed with 4% paraformaldehyde and stained with 0.1% crystal violet after 24-h incubation. Five random visual fields in each chamber were captured and counted under a microscope (Zeiss, German).

Immunohistochemistry (IHC)

IHC was further conducted to indirectly validate the potential target genes expression of miR-642a-3p and miR-145-5p. Briefly, 4-μm sections from tissue blocks of the included GBC patients in our hospital were taken on coated slides. Then the slides were deparaffinized, hydrated, and blocked with hydrogen peroxide. The appropriate buffer for each antibody was utilized for heat induced antigen retrieval. The slides were incubated with primary antibody, following with secondary antibody (anti-SYK antibody, anti-SH3GL1 antibody, anti-CDKN1A antibody, anti-MYC antibody, anti-VEGFA antibody, anti-EGFR antibody). Diaminobenzidine (DAB) chromogen and haematoxylin were used for staining and counterstaining, respectively.

Statistical analysis

Unpaired tailed student’s t tests, Kaplan–Meier survival analysis, and log-rank tests were used for statistical analyses with GraphPad Prism 8 (GraphPad Software, Inc., La Jolla, CA, USA). The results were shown as mean ± standard deviation (SD). P < 0.05 was considered statistically significant.

Results

Identification of miRNAs of different expression and their target genes

MiRNA array GSE104165 was retrieved from the GEO dataset (Additional file 1: Fig. S1). The volcano plot and heatmap with miRNAs clustering of different expression were portrayed in Fig. 1. A total of 122 miRNAs in GBC tumor samples were identified to be differentially expressed, among which the expression of 25 miRNAs was found to be significantly up-regulated and 97 miRNAs were down-regulated. The top 10 most up-regulated and down-regulated miRNAs were summarized in Tables 1 and 2, respectively. Then, MiRTarBase was applied to identify the possible target genes of the up- and down-regulated miRNAs. Notably, the 3rd up-regulated miRNA, miR-3676-5p, did not successfully predict target genes in the database, thus the top 2 miRNAs, including miR-4430 and miR-642-3p, were selected for further analysis. Similarly, miR-451a and miR-145-5p were the top 2 down-regulated miRNAs. MiRTarBase was used to identify the 477 potential target genes from the 2 most up-regulated miRNAs and 268 genes from the 2 most down-regulated miRNAs.

Fig. 1
figure 1

MiRNAs clustering of different expression. A Volcano plots of differentially expressed miRNAs. Green data points indicated down-regulated miRNAs and red data points indicated up-regulated miRNAs. The differences was set as |log FC|> 1. Black dots indicated miRNAs that were not differentially expressed. B The heatmap of differentially expressed miRNAs. Log FC of the genes was represented by the top-right legend. FC fold change

Table 1 The most up-regulated expression of miRNAs in GBC
Table 2. The most down-regulated expression of miRNAs in GBC

GO functional annotation analysis

Based on the potential target genes, 3 GO functional annotation analysis categories involving biological process (BP), cellular component (CC), and molecular function (MF) were conducted. The followings were enriched GO functions for the target genes of up-regulated miRNA: leukocyte degranulation, regulation of homotypic cell–cell adhesion, and centromeric sister chromatid cohesion in BP category; nucleolus, 90S pre-ribosome and nuclear origin of replication recognition complex in CC category; and phosphatidylinositol-5-phosphate binding, DNA replication origin binding and poly in MF category (Fig. 2A–C). The enriched GO functions for the 2 down-regulated miRNAs target genes were shown in Fig. 2D–F, involving the transcription from RNA polymerase II (RNAP II) promoter positive regulation, transcription positive reregulation, DNA-templated, transforming growth factor beta2 production in BP category regulation; nuclear chromatin, chromatin and protein kinase complex in CC category; and RNAP II transcription regulatory region sequence-specific binding, transcriptional activator activity, RNAP II core promoter proximal region sequence-specific binding and interleukin-6 receptor binding in MF category.

Fig. 2
figure 2

GO functions for the target genes of the up-regulated and down-regulated miRNAs. A Up-regulated miRNA enriched biological process. B Up-regulated miRNAs enriched cellular component. C Up-regulated miRNAs enriched molecular function. D Down-regulated miRNAs enriched biological process. E Down-regulated miRNAs enriched cellular component. F Down-regulated miRNAs enriched molecular function. GO gene ontology

Enrichment analysis in KEGG pathway

We subsequently performed a KEGG pathway enrichment analysis to investigate the target genes enriched pathways. As shown in Fig. 3A, the enriched KEGG pathways for up-regulation were in HTLV-1 infection, MAPK signaling pathway, and viral carcinogenesis. The enriched KEGG for down-regulated miRNAs included proteoglycans and microRNA in cancer, and the PI3K-Akt signaling pathway (Fig. 3B). mTOR pathway plays the central role in tumor progression of GBC [20, 21]. Therefore, 5 genes (AKT1S1, PTEN, MAPK1, PIK3R2, PRKCA) of up-regulated miRNAs and 13 genes (CAB39, IRS1, BRAF, TSC1, HIF1A, VEGFA, IKBKB, RPS6KA3, RRAGC, RPS6KB1, AKT1, MAPK1, EIF4E) of down-regulated miRNAs, which were involved in mTOR pathway and interacted with the selected miRNAs, may be associated with GBC progression.

Fig. 3
figure 3

The distribution of target genes of four selected miRNAs of different expression in GO-enriched functions. A Up-regulated miRNAs and (B) down-regulated miRNAs. Symbols of target genes were displayed on the left side of the graph. Gene involvement under GO terms was established by the colored connecting lines to the right. GO gene ontology. *P < 0.05; **P < 0.01

MiRNA-hub and PPI network construction

STRING online database was applied to generate a PPI network, based on the 4 selected miRNAs (top 2 most up-regulated and top 2 most down-regulated) with different expression as well as potential target genes. The top 10 hub genes were screened out according to node degrees (Table 3). The hub genes for the up-regulated miRNAs were identified as MARPK1, PTEN, SYK, TAF1, CDKN1A, POLR2E, PRKCA, CDK4, SH3GL1, and DNAJC10. Meanwhile, the hub genes for the down-regulated miRNAs were AKT1, VEGFA, MAPK1, EGFR, MYC, IL6, BCL2, ACTB, SMAD4, ESR1.

Table 3. Identification of hub genes by PPI network

Cytoscape software was utilized to construct the miRNA-hub gene network. We found that 6 hub genes (CDK4, DNAJC10, PRKCA, TAF1, POLR2E, and PTEN) could be regulated by up-regulated miR-4430, and 3 hub genes (SYK, SH3GL1, and CDKN1A) could be regulated by miR-642a-3p (Fig. 4A). In addition, miR-451a could potentially target 7 genes (SMAD4, ESR1, BCL2, MARPK1, IL6, MYC, and ACTB), and 3 hub genes (MYC, VEGFA, and EGFR) could be regulated by miR-145-5p (Fig. 4B).

Fig. 4
figure 4

The regulatory network between dysregulated miRNAs and hub genes. A Up-regulated miRNAs and (B) down-regulated miRNAs. The mRNA expression of predicted targets for miRNAs were from the TCGA database. C MiR-642a-3p: SYK expression, SH3GL1 expression, CDKN1A expression. D MiR-145-5p: MYC expression, VEGFA expression, EGFR expression. *P < 0.05

Validation of potential target genes using TCGA database

Next, we evaluated the potential targets of hub genes with the TCGA database. Based on short and long survival, SYK and SH3GL1, but CDKN1A, had a similar trend of overall survival like miR-642a-3p (Fig. 4C). The further analysis showed CDKN1A had a similar result with miR-642a-3p based on different T stages, which represented tumor invasion (Additional file 1: Fig. S2A). In the meantime, MYC and VEGFA, like miR-145-5p, were related to T stages (invasion-related) (Additional file 1: Fig. S2B-C; Fig. 4D).

Prognostic roles of miR-642a-3p and miR-145-5p and validation of predicted hub genes using TCGA database and SRRSH database

Further investigation was applied to analyze the expression and clinical significance of both miR-642-3p and miR-145-5p in GBC. GBC-SD and SGC-996 cell lines, primary GBC tissues samples were detected for the expression of miR-642-3p and miR-145-5p. It was confirmed that miR-642-3p were up-regulated and miR-145-5p were down-regulated in GBC cell lines (Fig. 5A, B). Using SRRSH database, we found that miR-642-3p and miR-145-5p were respectively up-regulated and down-regulated in GBC patients with poor prognosis via RNA extraction and qRT-PCR based on their tissue slides (Fig. 5C, D). IHC was further conducted to indirectly validate the potential target genes expression of miR-642a-3p and miR-145-5p. Similar results were observed in related target genes, such as SYK, SH3GL1, CDKN1A, MYC, VEGFA and so on (Additional file 1: Fig. S3).

Fig. 5
figure 5

The expression and prognostic roles of miR-642a-3p and miR-145-5p in GBC. A MiR-642a-3p expression and (B) miR-145-5p expression in 2 GBC cell lines (SGC-996, GBC-SD) were compared with HiBEC. C The expression of miR-642a-3p and (D) miR-145-5p in tumor tissues of GBC patients with different prognosis. E Kaplan–Meier curve of OS in CHOL patients with high vs. low expression of miR-642a-3p from TCGA database. F Kaplan–Meier curve of OS in CHOL patients with high vs. low expression of miR-145-5p for patients from TCGA database. G Kaplan–Meier curve of OS with high vs. low expression of miR-642a-3p in GBC patients from clinical data. H Kaplan–Meier curve of OS with high vs. low expression of miR-145-5p in GBC patients from clinical data. GBC gallbladder cancer, OS overall survival, CHOL cholangiocarcinoma, TCGA the Cancer Genome Atlas. *P < 0.05; **P < 0.01; ***P < 0.001; ****P < 0.0001

Survival analysis was performed to assess the prognostic value of miR-642-3p and miR-145-5p in GBC. Although miR-642-3p was not associated with worse OS in TCGA (Fig. 5E), high miR-642-3p presented decreased survival in our hospital data (Fig. 5G). It was implied that miR-642-3p may be a potential oncogene. By contrast, miR-145-5p was indicated as a favorable prognostic factor both in TCGA (Fig. 5F) and our data (Fig. 5H). Collectively, these data indicated that miR-642-3p and miR-145-5p were unusually expressed, which could be prognostic biomarkers and associated with invasion and metastasis in GBC.

Validation of potential target genes of miR-642a-3p and miR-145-5p in vitro

We used miR-642a-3p inhibitor and miR-145-5p mimics for further validation of potential target genes of miR-642a-3p and miR-145-5p, respectively. When using miR-642a-3p inhibitor, we found mRNA level of SH3GL1 and CDKN1A were decreased (Additional file 1: Fig. S4A). Moreover, when using miR-145-5p mimics, we found mRNA level of MYC, VEGFA, and EGFR were significantly regulated (Additional file 1: Fig. S4B). Besides, the relationship of the above-mentioned genes with tumor invasion have been reported by previous studies [22,23,24,25].

MiR-642a-3p and miR-145-5p on GBC migration and invasion in vitro

To examine the role of miR-642a-3p and miR-145-5p in GBC migration, a wound-healing assay was employed. RT-qPCR revealed significantly decreased miR-642a-3p and increased miR-145-5p levels of GBC cell lines after transfection, respectively (Additional file 1: Fig. S5). Our results demonstrated that decreased miR-642a-3p expression and increased miR-145-5p expression markedly suppress the migration ability of GBCs compared to negative control cells (Fig. 6A). Transwell invasion assay was performed to evaluate the influence of both miRNAs on GBC invasion. The results demonstrated the decrease of miR-642a-3p and over-expressed miR-145-5p suppressed cell invasion, respectively (P < 0.05, Fig. 6B–D). Hence, all data confirmed that both 2 miRNAs (miR-642a-3p and miR-145-5p) were key regulators relevant to invasion and metastasis of GBC, and approaches of targeting them may represent novel directions to improve the GBC prognosis.

Fig. 6
figure 6

MiR-642a-3p under-expression and miR-145-5p over-expression regulated GBC invasion and migration. AC After being transfected with miR-642a-3p inhibitor and miR-145-5p mimic, SGC-996 and GBC-SD showed weaker invasion than the control cells, respectively. D SGC-996 and GBC-SD were transfected with inhibitor NC, miR-642a-3p inhibitor, mimic NC, miR-145-5p mimic, respectively. Wound healing assay was performed in GBC cell with 48 h of recovery. GBC: Gallbladder cancer. *P < 0.05; **P < 0.01; ***P < 0.001; ****P < 0.0001. Error bars represented SD for n = 3

Discussion

In the study, miR-642a-3p and miR-145-5p were identified as invasion-metastasis associated miRNAs according to the bioinformatic analysis and experimental validation. As both two miRNAs were crucial in GBC invasion and metastasis, they would be utilized as promising targets for effective treatment to improve the prognosis of GBC patients.

Currently, miRNAs were gradually considered as essential regulators for tumor initiation, promotion, and progression transcriptional dynamics in GBC [26, 27], and thus researchers should pay more attention to systemically analyzing the invasion-metastasis associated miRNA in GBC. MiR-7-2-3p and miR-29c-3p were identified as metastasis suppressors for GBC via high-throughput screening, and they were also related to the pathogenesis of GBC [28]. Meanwhile, Ma et al. [29] demonstrated that miR-663a could regulate epithelial membrane protein-3 to suppresses GBC progression via interfering the MAPK/ERK pathway, indicating that the miR-663a/EMP3/MAPK/ERK axis would be potential treatment for GBC. Besides, after being activated by TGF-β1, miR-20a could play a crucial role in the pathogenesis and trigger metastasis of GBC through the miR-20a/Smad7/β-catenin axis [26]. Despite the identified miRNAs, further studies would be performed to offer more insights into improving the prognosis of GBC.

Indeed, a total of 4 miRNAs, including the top 2 up-regulated miRNAs (miR-642a-3p and miR-4430) and the top 2 down-regulated miRNAs (miR-145-5p and miR-451a), were selected for experimental validation. It was shown that miR-145-5p up-regulation and miR-642a-3p down-regulation regulated cell invasion and migration in GBC. Interestingly, the decrease in miR-4430 expression suppressed cell invasion and wound healing migration to enhance GBC loss-of-function, in the meantime, up-regulation of miR-451a significantly suppressed the invasion (Additional file 1: Fig. S6). Similar to our results, Ueta et al. [30] found that serum EVs miR-451a were significantly down-regulated and miR-451a inhibited GBC cell proliferation and induced apoptosis. Taken together, miR-451a would also be a novel therapeutic target for GBC patients.

In the miRNA-gene hub network, the hub genes were potentially regulated by miR-145-5p and miR-642a-3p. Various studies demonstrated that miR-145-5p was linked with various types of cancer, such as hepatocellular carcinoma [31], gastric cancer [32, 33], and upper tract urothelial carcinoma [34]. MiR-642a-3p was also linked to cell migration and invasion of hepatocellular carcinoma [35]. Besides, miR-29c-3p was found to be down-regulated (Table 2), which was in accordance with the results that up-regulation of miR-29c-3p could reverse EMT and decrease the metastasis ability in vitro and in vivo [28]. Similar to our results in Table 2, Jin et al. [36] also found that miR-143-3p suppressed tumor angiogenesis and growth of GBC through the ITGA6/PI3K/AKT/PLGF pathways. The aberrant expression of miRNAs have been relevant to GBC tumorigenesis and progression.

qRT-PCR results indicated that miR-145-5p and miR-642a-3p were significantly up- and down-regulated among GBC cell lines and clinical samples, respectively. MiR-145-5p up-regulation and miR-642a-3p down-regulation could significantly suppress in vitro activation, migration, and invasion of GBC using wound-healing and Transwell invasion assay. To explore the possible pathways, we performed GO annotation and KEGG pathway analysis for the predicted target genes of the top 2 most up- and down-regulated miRNAs using the Enrichr tool. Several target genes were revealed to be enriched in cell–cell adhesion, which was closely correlated to cell migration and invasion. Therefore, it would be a promising mechanism that can relieve the impacts of miR-145-5p and miR-642a-3p for GBC by targeting these genes.

The study has several limitations that need to be addressed. First, the public dataset with both miRNA and mRNA profiling in GBC was limited. Therefore, we tried to validate our results using different GBC databases (TCGA and SRRSH databases). Unfortunately, some verification results would not be obvious because of the small sample size. And further validation experiments should be performed with a large sample size. Second, merely 2 most up- and down-regulated miRNAs and the relevant target genes were included in the enrichment analysis. Further studies would be focused on more miRNAs of different expression to explore more potential treatment for GBC. Third, we selected 2 up- and down-regulated miRNAs for in vitro experiment validation, which would lead to the omittance of some functional miRNAs. And in vitro and in vivo experiments should be performed to reveal other functional miRNAs with a high validation accuracy. Finally, the deeper molecular mechanisms of GBC invasion and metastasis should also be explored in the future.

Conclusions

In the present study, miR-642a-3p and miR-145-5p were identified as invasion-metastasis associated miRNAs via bioinformatic analysis and experimental validation. Both miR-642a-3p (down-regulation) and miR-145-5p (up-regulation) would be served as potential therapeutic targets for GBC in the future.

Availability of data and materials

The original data of the study are available from the corresponding authors upon reasonable request.

Abbreviations

GBC:

Gallbladder cancer

BTC:

Biliary tract cancer

AJCC:

American Joint Committee on Cancer

EMT:

Epithelial–mesenchymal transition

PD-L1:

Programmed cell death-ligand 1

NCBI:

National Center for Biotechnology Information

GEO:

Gene Expression Omnibus

FC:

Fold change

GO:

Gene Ontology

KEGG:

Kyoto Encyclopedia of Genes and Genomes

PPI:

Protein–protein interaction

TCGA:

The Cancer Genome Atlas

OS:

Overall survival

HiBEC:

Human intrahepatic biliary epithelial cells

FBS:

Fetal bovine serum

qRT-PCR:

Quantitative real-time PCR

CT:

Cycle thresholds

IHC:

Immunohistochemistry

DAB:

Diaminobenzidine

SEM:

Standard error of mean

BP:

Biological process

CC:

Cellular component

MF:

Molecular function

RNAP II:

RNA polymerase II

References:

  1. Hundal R, Shaffer EA. Gallbladder cancer: epidemiology and outcome. Clin Epidemiol. 2014;6:99–109.

    PubMed  PubMed Central  Google Scholar 

  2. Rakić M, Patrlj L, Kopljar M, Kliček R, Kolovrat M, Loncar B, et al. Gallbladder cancer hepatobiliary. Surg Nutr. 2014;3:221–6.

    Google Scholar 

  3. Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics, 2021. CA Cancer J Clin. 2021;71:7–33.

    Article  Google Scholar 

  4. Valle JW, Kelley RK, Nervi B, Oh DY, Zhu AX. Biliary tract cancer. Lancet. 2021;397:428–44.

    CAS  Article  Google Scholar 

  5. Chun YS, Pawlik TM, Vauthey JN. 8th Edition of the AJCC cancer staging manual: pancreas and hepatobiliary cancers. Ann Surg Oncol. 2018;25:845–7.

    Article  Google Scholar 

  6. Mantripragada KC, Hamid F, Shafqat H, Olszewski AJ. Adjuvant therapy for resected gallbladder cancer: analysis of the national cancer data base. J Natl Cancer Inst. 2017;109.

  7. Chen M, Cao J, Bai Y, Tong C, Lin J, Jindal V, et al. Development and validation of a nomogram for early detection of malignant gallbladder lesions. Clin Transl Gastroenterol. 2019;10: e00098.

    Article  Google Scholar 

  8. Tauriello DVF, Palomo-Ponce S, Stork D, Berenguer-Llergo A, Badia-Ramentol J, Iglesias M, et al. TGFbeta drives immune evasion in genetically reconstituted colon cancer metastasis. Nature. 2018;554:538–43.

    CAS  Article  Google Scholar 

  9. Wu MJ, Chen YS, Kim MR, Chang CC, Gampala S, Zhang Y, et al. Epithelial–mesenchymal transition directs stem cell polarity via regulation of mitofusin. Cell Metab. 2019;29(993–1002): e6.

    Google Scholar 

  10. Zheng L, Xu M, Xu J, Wu K, Fang Q, Liang Y, et al. ELF3 promotes epithelial-mesenchymal transition by protecting ZEB1 from miR-141-3p-mediated silencing in hepatocellular carcinoma. Cell Death Dis. 2018;9:387.

    Article  Google Scholar 

  11. Civenni G, Bosotti R, Timpanaro A, Vazquez R, Merulla J, Pandit S, et al. Epigenetic control of mitochondrial fission enables self-renewal of stem-like tumor cells in human prostate cancer. Cell Metab. 2019;30(303–18): e6.

    Google Scholar 

  12. Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136:215–33.

    CAS  Article  Google Scholar 

  13. Shukla GC, Singh J, Barik S. MicroRNAs: processing, maturation, target recognition and regulatory functions. Mol Cell Pharmacol. 2011;3:83–92.

    CAS  PubMed  PubMed Central  Google Scholar 

  14. Liu N, Chang CW, Steer CJ, Wang XW, Song G. MicroRNA-15a/16-1 prevents hepatocellular carcinoma by disrupting the communication between Kupffer cells and regulatory T cells. Gastroenterology. 2021;162:575.

    Article  Google Scholar 

  15. Zhu M, Wei C, Lin J, Dong S, Gao D, Chen J, et al. UHRF1 is regulated by miR-124-3p and promotes cell proliferation in intrahepatic cholangiocarcinoma. J Cell Physiol. 2019;234:19875–85.

    CAS  Article  Google Scholar 

  16. Roshani Asl E, Rasmi Y, Baradaran B. MicroRNA-124-3p suppresses PD-L1 expression and inhibits tumorigenesis of colorectal cancer cells via modulating STAT3 signaling. J Cell Physiol. 2021;236:7071–87.

    CAS  Article  Google Scholar 

  17. Gene OC. The Gene Ontology (GO) project in 2006. Nucleic Acids Res. 2006;34:D322–6.

    Article  Google Scholar 

  18. Chen EY, Tan CM, Kou Y, Duan Q, Wang Z, Meirelles GV, et al. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinformatics. 2013;14:128.

    Article  Google Scholar 

  19. Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47:D607–13.

    CAS  Article  Google Scholar 

  20. Mohri D, Ijichi H, Miyabayashi K, Takahashi R, Kudo Y, Sasaki T, et al. A potent therapeutics for gallbladder cancer by combinatorial inhibition of the MAPK and mTOR signaling networks. J Gastroenterol. 2016;51:711–21.

    CAS  Article  Google Scholar 

  21. Zhang Y, Wang H, Chen T, Wang H, Liang X, Zhang Y, et al. C24-ceramide drives gallbladder cancer progression through directly targeting PIP4K2C to facilitate mTOR signaling activation. Hepatology. 2020;73:692.

    Article  Google Scholar 

  22. Takahashi T, Shivapurkar N, Riquelme E, Shigematsu H, Reddy J, Suzuki M, et al. Aberrant promoter hypermethylation of multiple genes in gallbladder carcinoma and chronic cholecystitis. Clin Cancer Res. 2004;10:6126–33.

    CAS  Article  Google Scholar 

  23. Wisnieski F, Calcagno DQ, Leal MF, Chen ES, Gigek CO, Santos LC, et al. Differential expression of histone deacetylase and acetyltransferase genes in gastric cancer and their modulation by trichostatin A. Tumour Biol. 2014;35:6373–81.

    CAS  Article  Google Scholar 

  24. Singh P, Jain SL, Sakhuja P, Agarwal A. Expression of VEGF-A, HER2/neu, and KRAS in gall bladder carcinoma and their correlation with clinico-pathological parameters. Indian J Pathol Microbiol. 2021;64:687–92.

    Article  Google Scholar 

  25. Shen H, He M, Lin R, Zhan M, Xu S, Huang X, et al. PLEK2 promotes gallbladder cancer invasion and metastasis through EGFR/CCL2 pathway. J Exp Clin Cancer Res. 2019;38:247.

    Article  Google Scholar 

  26. Chang Y, Liu C, Yang J, Liu G, Feng F, Tang J, et al. MiR-20a triggers metastasis of gallbladder carcinoma. J Hepatol. 2013;59:518–27.

    CAS  Article  Google Scholar 

  27. Srivastava K, Srivastava A, Mittal B. Common genetic variants in pre-microRNAs and risk of gallbladder cancer in North Indian population. J Hum Genet. 2010;55:495–9.

    CAS  Article  Google Scholar 

  28. Lu K, Feng F, Yang Y, Liu K, Duan J, Liu H, et al. High-throughput screening identified miR-7-2-3p and miR-29c-3p as metastasis suppressors in gallbladder carcinoma. J Gastroenterol. 2020;55:51–66.

    Article  Google Scholar 

  29. Ma Q, Zhang Y, Liang H, Zhang F, Liu F, Chen S, et al. EMP3, which is regulated by miR-663a, suppresses gallbladder cancer progression via interference with the MAPK/ERK pathway. Cancer Lett. 2018;430:97–108.

    CAS  Article  Google Scholar 

  30. Ueta E, Tsutsumi K, Kato H, Matsushita H, Shiraha H, Fujii M, et al. Extracellular vesicle-shuttled miRNAs as a diagnostic and prognostic biomarker and their potential roles in gallbladder cancer patients. Sci Rep. 2021;11:12298.

    CAS  Article  Google Scholar 

  31. Wang W, Ding B, Lou W, Lin S. Promoter hypomethylation and miR-145-5p downregulation-mediated HDAC11 overexpression promotes sorafenib resistance and metastasis of hepatocellular carcinoma cells. Front Cell Dev Biol. 2020;8:724.

    Article  Google Scholar 

  32. Zhou K, Song B, Wei M, Fang J, Xu Y. ANGPT2/MiR-145-5p suppresses the proliferation, migration and invasion of gastric cancer epithelial cells via the NOD_LIKE_RECEPTOR axis. Cancer Cell Int. 2020;20:416.

    CAS  Article  Google Scholar 

  33. Zhou T, Chen S, Mao X. miR-145-5p affects the differentiation of gastric cancer by targeting KLF5 directly. J Cell Physiol. 2019;234:7634–44.

    CAS  Article  Google Scholar 

  34. Hsu W, Li W, Lee Y, Huang A, Chang L, Lin H, et al. MicroRNA-145 suppresses cell migration and invasion in upper tract urothelial carcinoma by targeting ARF6. FASEB J Off Publ Fed Am Soc Exp Biol. 2020;34:5975–92.

    CAS  Google Scholar 

  35. Yu Z, Du Y, Li H, Huang J, Jiang D, Fan J, et al. miR-642 serves as a tumor suppressor in hepatocellular carcinoma by regulating SEMA4C and p38 MAPK signaling pathway. Oncol Lett. 2020;20:74.

    CAS  PubMed  PubMed Central  Google Scholar 

  36. Jin YP, Hu YP, Wu XS, Wu YS, Ye YY, Li HF, et al. miR-143-3p targeting of ITGA6 suppresses tumour growth and angiogenesis by downregulating PLGF expression via the PI3K/AKT pathway in gallbladder carcinoma. Cell Death Dis. 2018;9:182.

    Article  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by Scientific Research Fund of Zhejiang Provincial Education Department (Y202148325), and National Natural Science Foundation of China (81827804 and 81802330).

Author information

Affiliations

Authors

Contributions

JSC, HJS, JHH, XJC, and MYC designed the study and collected the data; JSC, AYF, BZ, RAJ, SJL, and TEC analyzed and interpreted the data; JSC, SJ, WT, YTT, and ZYL wrote the manuscript; XJC and MYC revised the manuscript. All authors read and approved the final version of the manuscript.

Corresponding authors

Correspondence to Xiujun Cai or Mingyu Chen.

Ethics declarations

Ethics approval and consent to participate

The current study was approved by the institutional ethics review board of Sir Run-Run Shaw Hospital, Zhejiang University, Hangzhou, Zhejiang Province, China (No. 20210625-32) (Additional file 2: Ethics file). The patients provided their written informed consent to participate in the study.

Consent for publication

Not applicable.

Competing interests

The authors deny any conflict of interest.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1: Table S1

. List of the primers used for miRNA quantitative real-time PCR. Figure S1. Data GSE104165 was shown. Figure S2. The comparison of mRNA expression of hub genes based on different T stages. (A) MiR-642a-3p: CDKN1A expression. (B) MiR-145-5p: MYC expression. (C) MiR-145-5p: VEGFA expression. *P. Figure S3. The comparison of immunohistochemistry of GBC tissue samples based on different survival. (A) SYK. (B) SH3GL1. (C) CDKN1A. (D) MYC. (E) VEGFA. (F) EGFR. GBC: Gallbladder cancer. *P. Figure S4. The mRNA expression of predicted hub genes in GBC-SD cell. (A) The expression of SH3GL1, CDKN1A, and SYK in negative control vs. miR-642a-3p inhibitor. (B) The expression of MYC, EGFR, and VEGFA in negative control vs. miR-145-5p mimics. GBC: Gallbladder cancer. *P. Figure S5. (A-B) After being transfected with 50 nM mimics, the expression of miR-642a-3p and miR-145-5p were significantly elevated in SGC-996 and GBC-SD cell lines. qRT-PCR was applied to detect the miR-642a-3p and miR-145-5p expression levels after 48h transfection. GBC: Gallbladder cancer; qRT-PCR: Quantitative real-time PCR. *P . Figure S6. Down-regulation of miR-4430 and up-regulation of miR-451a regulated GBC cell invasion and metastasis. (A) SGC-996 and GBC-SD, which were transfected with miR-4430 inhibitor and miR-451a mimic, invaded less versus control cancer cells, respectively. (B) Quantification of SGC-996 and GBC-SD after miR-4430 inhibitor. (C) Quantification of SGC-996 and GBC-SD after miR-451a. (D) SGC-996 and GBC-SD were transfected with inhibitor NC, miR-4430 inhibitor, mimic NC, miR-451a mimic, respectively. Wound healing assay was performed in GBC cell with 48h of recovery. GBC: Gallbladder cancer. NC: Negative control. ***P.

Additional file 2

. Ethics file.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visithttp://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Cao, J., Shao, H., Hu, J. et al. Identification of invasion-metastasis associated MiRNAs in gallbladder cancer by bioinformatics and experimental validation. J Transl Med 20, 188 (2022). https://doi.org/10.1186/s12967-022-03394-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-022-03394-8

Keywords

  • Gallbladder cancer
  • MiRNA
  • Invasion
  • Metastasis
  • Bioinformatics