Abstract

Pancreatic cancer (PC) has a high degree of malignancy and poor prognosis, and countless patients have distant metastasis when diagnosed. Gemcitabine (GEM) chemotherapy is one of the main ways of treatment. However, PC cells have been displayed chemoresistance to GEM during treatment. Circular RNAs (circRNAs) have been demonstrated to be the most popular diagnostic and prognostic biomarkers in PC with GEM resistance. Here, we assessed the potential of circLMTK2 in the GEM resistance of PC cells. Functional assays were implemented to measure the impacts of circLMTK2 on the proliferation, migration/invasion, and apoptosis of GEM-resistant PC cells. Bioinformatics analysis and mechanical experiments displayed the underlying mechanism of circLMTK2 in GEM-resistant PC cells. We found that circLMTK2 was upregulated in PC and GEM-resistant PC tissues and cells. CircLMTK2 knockdown suppressed proliferation, invasion, migration, and enhanced apoptosis in GEM-resistant PC cells. Moreover, circLMTK2 silencing could decrease GEM resistance-associated tumor size in vivo. In terms of mechanism, circLMTK2 served as a sponge for miR-485-5p, and miR-485-5p bound to p21 (RAC1) activated kinase 1 (PAK1), which were clarified via the dual-luciferase assay in PC cell lines. We confirmed that circLMTK2 knockdown attenuated GEM-resistant PC cells by regulating PAK1 via miR-485-5p. Our study demonstrated that circLMTK2 may be a novel diagnostic and prognostic biomarker in GEM-resistant PC cells.

1. Introduction

Pancreatic cancer (PC), a malignant tumor with a high degree of malignancy and poor prognosis, has a 5-year survival rate less than 7% [1]. Numerous patients have been diagnosed with distant metastasis and thus lost the chance of surgery [2]. Therefore, chemotherapy is considered to be the main treatment for improving the prognosis of patients with PC. Gemcitabine (GEM) is a first-line treatment for locally advanced and metastatic PC [3]. However, GEM resistance is an important factor limiting the survival time of patients with PC [4]. Therefore, it is of great importance to investigate the molecular mechanism of GEM resistance in PC and explore the key regulatory targets that mediate GEM resistance.

Although many signaling pathways and downstream genes have been confirmed to be involved in the regulation of GEM resistance in PC, such as PI3K/AKT, Wnt, ERK, HMGB1, MDR1, and so on [58]. However, the regulatory mechanism of GEM resistance in PC is still unclear. The occurrence of chemoresistance in PC is related to the imbalance of key gene expression and epigenetic regulation network [9]. CircRNAs are a novel class of noncoding RNAs that are circularized by connecting the 3′-5′. Studies have confirmed that circRNAs have the advantages of stable structure, high conservation, and tissue-specific expression [10]. Various studies have shown that circRNAs, as an important part of gene regulatory network and epigenetic regulatory network, play important roles in chemoresistance in cancer therapy [11, 12]. For instance, circ_0058063 could regulate the cisplatin-resistance of bladder cancer cells via the miR-335-5p/B2M axis [13]. Besides, high circFARP1 levels were confirmed positively related to elevated serum LIF levels and poor survival in pancreatic ductal adenocarcinoma. Silencing circFARP1 expression could suppress pancreatic ductal adenocarcinoma growth and GEM resistance [14]. Therefore, further studies on the role and mechanism of circRNAs in chemotherapy resistance of PC are important.

In solid tumors, circRNA mainly enhances or inhibits the expression of its target genes by binding with microRNAs (miRNAs), and regulates the metabolic functions of cell proliferation, differentiation, and apoptosis [15]. Based on the above research background, we first systematically analyzed the GEM resistance-related circRNAs/miRNA/mRNA network in PC using bioinformatics analysis. Through interaction network analysis and previous study results, the circLMTK2-miR-485-5p-PAK1 axis attracted our attention.

In this study, we detected the expression of circLMTK2 in PC and evaluate its impact on prognosis in PC patients. We also explored the regulatory effects of circLMTK2 on the proliferation, apoptosis, invasion, and migration in GEM-resistant PC cells. In addition, we investigated the effect of circLMTK2 on GEM-resistant tumor growth in vivo. Finally, we clarified the potential molecular mechanism of circLMTK2 in GEM-resistant PC cells.

2. Methods

2.1. Patients

The study was approved by the Ethics Committee of the Lihuili Hospital affiliated to Ningbo University and obtained the written informed consent of all subjects. The cancer tissues and adjacent tissues of 40 patients with PC were obtained from the Lihuili Hospital affiliated to Ningbo University during surgery. GEM-resistant: PC patients had no effect on GEM treatment, or the tumor progressed after previous GEM treatment (no more than 6 months).

2.2. Cell Culture and Transfection

Human PC cells (MIA-PaCa-2, PANC-1, BxPC-3, and SW 1990) were obtained from Cell Bank of Chinese Academy of Sciences (Shanghai, China). MIA-PaCa-2 and PANC-1 cells were maintained in Dulbecco’s Modified Eagle’s Medium. BxPC-3 cells were maintained in RPMI-1640 medium. SW 1990 cells were maintained in Leibovitz’s L-15 medium. All cell mediums were cultured in a 5% CO2 incubator at 37°C. PANC-1-GEM and MIA-PaCa-2-GEM were created according by previous study [16]. sh-circLMTK2, miR-485-5p mimic/inhibitor, sh-PAK1/pcDNA3.1 PAK1, or negative control (NC) were transfected into PANC-1-GEM and MIA-PaCa-2-GEM cells using Lipo 3000 (Thermo Fisher Scientific, USA) according to the manufacture’s protocol.

2.3. Real-Time Quantitative PCR Analysis

Total RNA from samples was extracted using TRizol reagent (Thermo Fisher Scientific, USA), and was further reverse transcribed into cDNA. Fast SYBR Green Master Mix (Thermo Fisher Scientific, USA) was then used to carry out qRT-PCR on Real-Time System (Bio-Rad, USA). Each sample was set to 3 repetitions; the expressions of circLMTK2, miR-485-5p, and PAK1 were analyzed using the 2−ΔΔCT method. Primer information is provided in Table 1.

2.4. Western Blot Analysis

The protein was extracted using lysis buffer, and protein concentration was measured via BCA assay kits. After electrophoresis and membrane transfer, the bands were incubated with the primary antibody including PAK1 (1/1000, ab223849, Abcam) and GAPDH (1/2500, ab9485, Abcam) at 4°C overnight, then cultured with secondary antibody at room temperature for 2 h. The chemiluminescence system (Bio-Rad, USA) imaging system was performed to measure all bands.

2.5. Subcellular Localization Assay

Nuclear, cytoplasmic, and total RNA was isolated using a PARIS™ kit (Invitrogen), followed by qRT-PCR analysis. U6 was used as endogenous controls for the nucleus, while GAPDH was used for the cytoplasm.

2.6. Cell Counting Kit-8 (CCK-8) Assay

Cells were seeded in a 96-well plate for 24 h and were transfected with the indicated plasmids. Ten microliters of CCK-8 solution (Beyotime, Shanghai, China) was then added into each well for 2 h at 37°C. The absorbency was measured at 450 nm by using a microplate reader.

2.7. Wound Healing and Transwell Assay

The adherent cells were scraped vertically using the 200-μL pipette and cultured with the medium containing 1% FBS. The wound healing was confirmed under an inverted microscope (Zeiss, Japan) at 0 and 24 hours, and photos were taken. In the transwell assay, total 1 × 104 transfected cells were added and cultured on the upper chamber with Matrigel (BD Bioscience, CA, USA) coated, and the bottom chamber was seeded 500 μL medium with 10% FBS. After 48 h incubation, the upper chamber was removed and the bottom chamber of cells was fixed using formaldehyde. The invaded cell number was counted via an inverted microscope (Zeiss, Japan).

2.8. Terminal deoxynucleotidyl transferase-mediateddUTP-nick end labeling (TUNEL)

Based on the instructions, cell apoptosis was detected by the TUNEL detection kit (Vazyme).

2.9. Luciferase Reporter Assay

First, construction of miR-485-5p and circLMTK2 wild-type (WT) or circLMTK2 mutant-type (MUT) plasmids, miR-485-5p, and PAK1 wild-type (WT) or PAK1 mutant-type (MUT) plasmids were performed. The recombinant plasmid was then transfected into PC cells with miR-485-5p mimic or mimic NC. After 48 h, luciferase activities were measured using Dual-Luciferase Reporter Assay System (Promega) according to manufacturer’s protocol.

2.10. Mouse Xenograft Assay

Total 100 μl of sh-circLMTK2 or sh-NCGEM-resistant PANC-1 cells (2 × 107/mL) were injected into the abdominal cavity of BALB/c nude mice. The tumor formation could be touched after 1 week, and then the tumor size was measured every week, and the tumor body was cut and weighed at the fourth week. The animal experiments were approved by the Institutional Animal Care and Use Committee of Lihuili Hospital affiliated to Ningbo University and were implemented following the institutional guidelines.

2.11. Statistical Analysis

The SPSS 22.0 software was conducted for data analysis in this study. Mean ± SD was used to represent the measurement data. Student’s t-test was performed to analyze statistical analysis between two groups. ANOVAs with Bonferroni post hoc tests were used for more than three groups. The survival rate was analyzed by Kaplan–Meier curve. was considered statistical significance.

3. Results

3.1. circLMTK2 is Increased in PC and GEM-Resistant PC Tissues and Cells

From the GEO database, we analyzed the GSE110580 dataset, the GSE79234 dataset, and the GSE80617 dataset. The screening threshold for significant differences in gene expression was and |log2fc| > 0.585, and the heat map and volcano map or scatter plot of genes with significant differences in this analysis were as follows (Figure 1(a)1(f)).

After obtaining significantly different circRNA, miRNA, and mRNA, we tried to construct a circRNA-miRNA-mRNA interaction network which is displayed in Figures 1(g)-1(h). For the significantly different genes located in the two ceRNA networks, the protein-protein interaction (PPI) interaction network was analyzed by string database and visualized by Cytoscape software (Figure 1(i)). The enriched gene ontology (GO) terms and enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of the differentially expressed genes are displayed in Figures 1(j) and 1(k). Finally, the differential genes were mapped in significantly enriched pathways using Pathview package (Figure 1(l)).

Considering that PAK1 was abnormally expressed in pancreatic cancer and could inhibit DNA damage and inhibit cell apoptosis via the NF-κB pathway, thereby making pancreatic cancer cells resistant to GEM in the previous study, we hypothesized that circLMTK2 might regulate the expression of PAK1 by binding to miR-485-5p through the circRNA-miRNA-mRNA interaction network. Thus, we used qRT-PCR to investigate the circLMTK2 expression in human PC and GEM-resistant PC tissues. The results indicated that circLMTK2 was significantly increased in PC and GEM-resistant PC tissues (Figure 2(a) and 2(b)). Besides, the Kaplan–Meier survival analysis displayed that higher circLMTK2 expression led to poorer survival rate (Figure 2(c)). Interestingly, we also found that circLMTK2 expression was increased in PC cell lines (Figure 2(d)). PANC-1 and MIA-PaCa-2 cell lines were selected for the following study due that both cell lines harbored higher circLMTK2 expression than other PC cell lines. After creating successful PANC-1 and MIA-PaCa-2-GEM-resistant cells, we detected circLMTK2 expression was further increased in GEM-resistant PC cell lines (Figure 2(e)). In addition, the amplified product was sequenced using Sanger sequencing to validate the circularized junction of circLMTK2 (Figure 2(f)). To study the location of circLMTK2, we conducted the subcellular fractionation assay. The result indicated that circLMTK2 was mainly expressed in cytoplasm (Figure 2(g)).

3.2. circLMTK2 Knockdown Suppresses Proliferation, Invasion, and Migration, and Enhances Apoptosis in GEM-Resistant PC Cells

To study the function of circLMTK2 in GEM-resistant PC cells, we first successfully inhibited circLMTK2 expression via transfecting sh-circLMTK2-1 and sh-circLMTK2-2 into GEM-resistant PC cells. The CCK-8 assay demonstrated that the knockdown of circLMTK2 significantly attenuated cell proliferation in GEM-resistant PC cells (Figure 3(a)). Besides, wound healing and transwell assays indicated circLMTK2 knockdown could significantly suppress cell migration and invasion in GEM-resistant PC cells (Figure 3(b)3(g)). In addition, the knockdown of circLMTK2 promoted GEM-resistant PC cells’ apoptosis as exhibited by the Tunel assay (Figure 3(h)3(j)).

3.3. circLMTK2 Directly Binds to miR-485-5p

After then, the binding sites of miR-485-5p in circLMTK2 were predicted and mutated accordingly (Figure 4(a)), and we found miR-485-5p was downregulated significant in PC and GEM-resistant PC tissues, as well as PC cell lines (Figure 4(b)4(d)). Then, the luciferase reporter assay was used to determine the physical connection between circLMTK2 and miR-485-5p in PANC-1 and MIA-PaCa-2 GEM-resistant cells. We found that miR-485-5p mimics reduced the luciferase activity of circLMTK2 wild-type but had no impacts on that of circLMTK2 mutant sites (Figure 4(e)).

To investigate the role of miR-485-5p in PANC-1 and MIA-PaCa-2 GEM-resistant cells, we tried to knock down miR-485-5p expression via miR-485-5p inhibitor (Figure 4(f)). Interestingly, the CCK-8 assay indicated miR-485-5p inhibitor could reverse the antiproliferation effect of sh-circLMTK2 (Figure 4(g)). Besides, the suppression of migration and invasion of PANC-1 and MIA-PaCa-2 GEM-resistant cells induced by sh-circLMTK2 could be revised by miR-485-5p inhibitor (Figure 4(h)4(j)). In addition, the Tunel assay showed that miR-485-5p inhibition could restore the enhanced apoptosis caused by silencing circLMTK2 (Figure 4(k)).

3.4. PAK1 is a Direct Target of miR-485-5p

Subsequently, the binding site between miR-485-5p and the 3′-UTR of PAK1 is shown in Figure 5(a). Interestingly, PAK1 expression was elevated in PC cell lines (Figure 5(b)). Then, the luciferase reporter assay was used to determine the physical connection between PAK1 and miR-485-5p in PANC-1 and MIA-PaCa-2 GEM-resistant cells. We found that miR-485-5p mimics reduced the luciferase activity of PAK1 wild-type but had no impacts on that of PAK1 mutant sites (Figure 5(c)). We next treated PANC-1 and MIA-PaCa-2 cells with miR-485-5p mimic and found that the expression of PAK1 was significantly downregulated in miR-485-5p mimic groups than mimic NC groups (Figure 5(d)). To investigate the role of PAK1 in PANC-1 and MIA-PaCa-2 GEM-resistant cells, we tried to inhibit PAK1 via transfecting with sh-PAK1 (Figure 5(e)). The CCK-8 assay indicated that PAK1 inhibition significantly weakened cell proliferation in PANC-1 and MIA-PaCa-2 GEM-resistant cells (Figure 5(f)). In addition, silencing PAK1 could suppress the migration and invasion of PANC-1 and MIA-PaCa-2 GEM-resistant cells (Figure 5(g)5(i)). Moreover, the TUNEL assay demonstrated that knockdown of PAK1 could enhance apoptosis in PANC-1 and MIA-PaCa-2 GEM-resistant cells (Figure 5(j)).

3.5. PAK1 Overexpression Eliminates the Effect of circLMTK2 Knockdown

Next, we verified whether circLMTK2 exerted its role in PANC-1 and MIA-PaCa-2 GEM-resistant cells via the miR-485-5p/PAK1 pathway through functional rescue experiment. The qRT-PCR assay was used to confirm the efficiency of pcDNA3.1/PAK1, and the inhibitory effect of silencing circLMTK2 on PAK1 expression could be revised by pcDNA3.1/PAK1 transfection (Figure 6(a) and 6(b)). The CCK-8 assay showed that cell proliferation attenuated via sh-circLMTK2 could be restored by PAK1 overexpression (Figure 6(c)). Besides, the suppression of migration and invasion of PANC-1 and MIA-PaCa-2 GEM-resistant cells induced by sh-circLMTK2 could be revised by PAK1 overexpression (Figure 6(d)6(f)). In addition, the TUNEL assay showed that PAK1 overexpression could restore the enhanced apoptosis caused by silencing circLMTK2 (Figure 6(g)).

3.6. circLMTK2 Knockdown Inhibits the Growth of PC Tumor In Vivo

Finally, to investigate the role of circLMTK2 in PC in vivo, we established the xenograft model by injecting PANC-1-GEM-resistant cells transfected with sh-circLMTK2 or sh-NC into nude mice. The results indicated that the tumor size and tumor volume in the sh-circLMTK2 group were decreased significantly than that in sh-NC groups (Figure 7(a) and 7(b)). In addition, the circLMTK2 and PAK1 expression were decreased while the miR-485-5p expression was increased in tumor tissues in the sh-circLMTK2 group compared with those in the sh-NC group, respectively (Figure 7(c)7(f)).

4. Discussion

Although GEM has been considered as the first-line therapeutic agent for pancreatic cancer, the drug resistance of GEM has become an important factor affecting the chemotherapy effect of pancreatic cancer [17, 18]. Using bioinformatics analysis, we analyzed previous high throughput datasets of GEM resistance, revealed their network relationships, and obtained the circLMTK2/miR-485-5p/PAK1 axis that might play a regulatory role in pancreatic cancer.

Studies have shown that the dysregulated expression of circRNAs in tumors can not only be used as biological markers for diagnosis and prognosis but also as potential intervention targets for tumor treatment and chemoresistance. For example, circBFAR was upregulated in pancreatic ductal adenocarcinoma, and the ectopic expression of circBFAR was positively associated with poorer prognosis of pancreatic ductal adenocarcinoma patients [19]. Similarly, circRNAs are also playing important roles in the field of chemotherapy resistance, Liu et al. revealed that circHIPK3 could enhance GEM resistance in PC cells through targeting RASSF1 via sponging miR-330-5p [20]. In this study, we confirmed that circLMTK2 was significantly increased in PC and GEM-resistant PC tissues. Kaplan–Meier survival analysis showed that circLMTK2 expression was significantly negatively correlated with survival. Furthermore, circLMTK2 knockdown could inhibit cell proliferation, invasion/migration, and enhances apoptosis in GEM-resistant PC cells. Our study suggests that elevated expression of circLMTK2 is an independent risk factor for poor prognosis of PC and may be a prognostic marker for PC, and a potential therapeutic target for GEM chemotherapy in PC.

In terms of the regulatory mechanism of circRNAs on tumor cell proliferation, differentiation, and chemoresistance, a large number of studies have confirmed that circRNAs can affect their transcription and translation by sponging miRNAs, and then regulate the expression of related tumor genes, leading the occurrence and development of cancers through various related signal transduction pathways [21]. In this study, biological information analysis and experimental verification reveal that miR-485-5p bound to circLMTK2. miR-485-5p has been confirmed to play an important role in the development of cancer in previous studies [22, 23]. However, the function role of miR-485-5p in PC is still unknown. Through validation, we found that circLMTK2 could negatively regulate the expression of miR-485-5p, thereby mediating the proliferation, invasion/migration, and apoptosis of PC cells. Next, we further revealed the downstream target PAK1 of the circLMTK2/miR-485-5p axis, which has been proved to be an important therapeutic target for PC [24]. PAK1 inhibitors have also been clarified that can inhibit the progression of PC [25]. More importantly, PAK1 has been reported to promote gefitinib resistance in lung cancer cells [26]. Ghosh et al. have pointed that PAK1 contributes to tamoxifen resistance in breast cancer [27]. Consistently, our research also unveiled that silencing of PAK1 could improve GEM resistance in PC cells.

In conclusion, our study confirmed that silencing of circLMTK2 could markedly attenuate cell proliferation, invasion/migration, and apoptosis in GEM-resistant PC cells. We also demonstrated that circLMTK2 knockdown inhibited GEM-resistant associated cell proliferation, migration, invasion, and promoted cell apoptosis through the miR-485-5p/PAK1 axis. Furthermore, compared with the control group, the GEM-resistant tumor xenografts of circLMTK2 knockdown group were significantly smaller. Our results suggest that the circLMTK2/miR-485-5p/PAK1 axis may be a potential therapeutic target to improve the sensitivity response of PC to GEM.

Data Availability

The original data of this manuscript can be obtained from the corresponding author under reasonable request.

Conflicts of Interest

The authors declare that they do not have any conflicts of interest.

Authors’ Contributions

Yeting Lu and Shuping Zhou contributed equally to this work.

Acknowledgments

This research was supported by the Zhejiang Provincial Natural Science Foundation of China under Grant no. LQ21H160012.