Differential and survival analyses for miR-1301-3p
We extracted miR-1301-3p expression data from three PC-miRNA expression profiles (GSE31568, GSE41372, and GSE32688) of the Gene Expression Omnibus (GEO) database. Next, we compared the difference of miR-1301-3p between PC sample tissues and healthy tissues. To verify the prognostic importance, the clinical data and miR-1301-3p expression value were also obtained based on the PC cohort of The Cancer Genome Atlas (TCGA) database. According to the median value of miR-1301-3p, we divided the PC patients into high and low-level groups and performed Kaplan-Meier survival analyses. A P-value was calculated by the log-rank test.
Patients and samples in the validation cohort
Between February 2018 and August 2020, a total of 72 PC patients were enrolled in the validation cohort, who were admitted to the Department of Hepatobiliary Surgery in Beijing Chao-Yang Hospital. These patients had not accepted radiotherapy or chemotherapy preoperatively, and the final diagnosis was determined by pathological results. Paired PC sample tissues and adjacent normal tissues were immediately stored in liquid nitrogen for two hours and then transferred into − 80 °C refrigerator for storage. Postoperatively, these patients have followed an average of 12 months, ranging from two to 29 months. The ethics committee of Beijing Chao-Yang Hospital approved this study, and all the patients signed the informed consent form.
Function annotation and signaling pathway enrichment for miR-1301-3p
We first applied the miRWalk2.0 database to predict the binding genes of miR-1301-3p and then performed correlation analyses between mRNA expression of these genes and miR-1301-3p value based on the PC cohort of TCGA. Finally, the combinative genes negatively correlated with miR-1301-3p were regarded as the target genes of miR-1301-3p.
To understand the functions of these target genes, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses using the clusterProfiler R package[10]. A protein-protein interaction (PPI) network was constructed using the Search Tool for the Retrieval of Interacting Genes database and visualized with Cytoscape software. The correlation of each PPI relationship pair was represented by a combined score ranging from 0 (low) to 1 (high). In our study, an interaction > 0.4 (moderate) was applied as the cut-off value. The Molecular Complex Detection plug-in in Cytoscape software was used to identify the hub genes among the PPI network. The screening conditions were set as degree cut-off = 2, K-Core = 2, and Node Score Cutoff = 0.2. Besides, we verified the hub gene expression using the PC-mRNA data of GSE16515.
Quantitative real-time polymerase chain reaction
TRIzol (Invitrogen, USA) was used for total RNA extraction and qualified using the NanoDrop ND-1000 (Thermo Fisher, USA). Total RNA was converted to first-strand cDNA according to the manufacturer’s protocol (rtStarTM First Strand cDNA Synthesis Kit, Arraystar lnc.). Specific primers for miR-1301-3p was designed by RiboBio (Guangzhou, China).The sequence of forward and reverse primers for miR-1301-3p was 5’- ACACTCCAGCTGGGTTGCAGCTGCCTGGGAGT-3’ and 5’-CTCAACTGGTGTCGTGGAGTCGGCAATTCAGTTGAGGAAGTCAC-3’. qRT-PCR was performed using Arraystar SYBR® Green Real-time qPCR Master Mix (Arraystar lnc) according to the manufacturer’s instruction. The relative expression of miR-1301-3p was calculated using the 2−ΔΔCt method and normalized to β-actin expression levels.
Cell culture and transfection
Five PC cell lines were selected to test the expression of miR1301-3p, including SW1990, AsPC-1, CFPAC-1, PANC-1, and Patu-8988. These cells were purchased from the American Type Culture Collection (Manassas, VA, USA). We cultured the PC cells with DMEM medium (Biological Industries) containing penicillin/streptomycin and 10% fetal bovine serum (FBS) at 37 °C with 5% CO2. Since miR-1301-3p high expressed in SW1990 and PANC-1 cells (Fig. 2A), we selected these two cell lines to conduct further experiments.
MiR1301-3p mimics, inhibitor, and negative control (NC) were designed and synthesized by RiboBio (Guangzhou, China). Small interference RNA against RhoA (si-RhoA)-535, 575, and 628 were designed and synthesized by GenePharma (Shanghai, China). In brief, miR-1301-3p mimics and inhibitor (100 nM) were transfected into PANC-1 and SW1990 cells, using Lipofectamine 3000 reagent (Invitrogen, USA), following the manufacturer's instruction.
Cell Counting Kit-8 (CCK-8) assays
According to the manufacturer’s instruction, we performed CCK-8 assays (Sigma, Aldrich) to examine the proliferation ability of PC cells. Approximately 2 × 103 cells were added to each well of the 96-well plate, and then, the plate was cultured for 24 h at 37 °C. Next, we added 50 µl of the miR-1301-3p mimics, inhibitor, and NC to wells of the 96-well plate. Then the plate was placed in the 37 °C incubator again for 24 h. At 0, 24, 48, 72, and 96 h, 10 µl CCK-8 solution was added into each well. After 2 h, optical density (OD) 450 nm values were measured using the enzyme-labeled instrument (Bio-Rad, United States). Cells were tested three times for each group.
Transwell migration and invasion assays
We conducted transwell migration assays with a 6.5 chamber with 8 um pores (Corning, NY, USA). A suspension containing 1 × 104 PANC-1 and SW1990 cells were prepared and suspended separately in serum-free DMEM and added into the upper chamber. After that, 500 ul of 10% serum-containing DMEM was added into the lower chamber of the well and incubated 24 h at 37 °C. After 24 h, PC cells in the upper chamber were removed. For transwell invasion assays, 1 × 104 cells were seeded on the Matrigel-coated membrane inserts. We used 0.1% crystal violet to stain the cells that migrated into the underside of the membrane. Four random fields were selected at 4 × magnification for counting cell numbers. Each experiment was performed three times.
Wound healing assays
SW1990 and PANC-1 cells (1 × 105) were incubated in six-well culture plates for 48 h until the cells were 80–90% confluent. Cells were maintained in 10% FBS containing DMEM media for 24 hours. PBS was used to wash away the non-adherent cells. A sterile 200-µl pipet was used to make a scratch in the center of the cell monolayer. The monolayer was washed three times with PBS, and fresh media was added. After 0 h, 24 h, and 48 h, the wound width was measured at 2.5 × magnification. Each assay was performed three times.
Protein extraction and western blotting assays
Total protein was extracted from PC cells after 72 h transfection, and the BCA protein assay kit (Beyotime, China) was used to measure protein concentration followed by the manufacturer's instruction. Briefly, 12% SDS-PAGE was used for electrophoresis, and then the proteins were transferred to PVDF membranes. Glyceraldehyde-3-Phosphate Dehydrogenase (GAPDH), RhoA, E-cadherin, N-cadherin, and Vimentin antibodies were used to analyze total protein. Specific primary rabbit-anti-human antibodies (CST, 1:1000) were used to incubate the membranes at 4 °C overnight. On the second day, the membranes were incubated with HRP-conjugated anti-rabbit IgG antibodies (1:2000) at room temperature for 1 h. An enhanced chemiluminescence detection system was used to visualize the bands. GAPDH was used as an internal control. Rabbit anti-GAPDH, RHOA, N-cadherin, E-cadherin, and vimentin antibodies (Cell Signaling Technology, Danvers, MA, USA) were used to analyze cell lysates.
Luciferase reporter assays
According to starBase network tool, RhoA is a predicted target for miR1301-3p. The binding site between miR1301-3p and 3’UTR of RhoA was evaluated by using the pmirGLO dual-luciferase miRNA expression vector containing wild type (WT) or mutant (MUT) 3’UTR of RhoA. The WT or MUT 3’UTR of RhoA and miR-1301-3p mimics were co-transfected into PANC-1 cells. After 48 h, the luciferase reporter assay system was used to examine the luciferase activity. Each experiment was performed in triplicate.
Statistical analysis
R software version 3.6.0 was used to perform statistical analyses. Continuous variables between the two groups were compared by a paired sample t-test. The data were presented as the mean ± standard deviation. Qualitative data were analyzed by the chi-square test. Fisher's exact test was employed to compare the categorical variables among groups. GraphPad Prism 8.0 (GraphPad Software, lnc. La Jolla, CA, USA) was applied to produce figures. P-value < 0.05 was considered statistically significant.