The detailed biological functions of circular RNA (circRNA) are largely unexplored. Using circRNA sequencing, we identified 169 differentially expressed circRNA in pancreatic ductal adenocarcinoma (PDAC) cells compared with nontumor human pancreatic ductal epithelial cells. Among them, circFOXK2 was validated with significant upregulation in PDAC cells and 63% of primary tumors (53 of 84). circFOXK2 promoted cell growth, migration, and invasion and was involved in cell-cycle progression and apoptosis. circFOXK2 contained multiple miRNA binding sites, functioning as a sponge for miR-942, which in turn promoted expression of ANK1, GDNF, and PAX6. A novel and highly specific circRNA-pulldown followed by mass spectrometry analysis identified 94 circFOXK2-interacting proteins, which were involved in cell adhesion, mRNA splicing, and structural molecule activity. Of these, circFOKX2 interactions with YBX1 and hnRNPK enhanced expression of oncogenes NUF2 and PDXK. Knockdown of circFOXK2 reduced binding of YBX1 and hnRNPK to NUF2 and PDXK, in turn decreasing their expression. Collectively, our findings demonstrate that circFOXK2 in complex with YBX1 and hnRNPK promotes expression of oncogenic proteins that contribute to PDAC progression.

Significance:

This study reveals a prominent role for the circRNA circFOXK2 in PDAC progression, suggesting that circFOXK2 might be a novel diagnostic marker for PDAC.

Pancreatic ductal adenocarcinoma (PDAC) is the fourth leading cause of cancer death in the United States and Europe (1, 2). Surgery, chemotherapy, and radiotherapy are standard therapy for PDAC. Unfortunately, patients with PDAC can only extend survival up to several months (3, 4). Also, combining gemcitabine with other chemotherapy drugs such as nab-paclitaxel can only provide a limited extension of survival with developing adverse effects including leukopenia and peripheral neuropathy (5). Therefore, novel therapeutic targets are urgently needed for treating PDAC.

Circular RNA (circRNA), which was first discovered in 1976, is a type of noncoding RNA with unknown function described previously. circRNA is formed by the linking of 3′ end of exon back to its 5′ end, forming a circular structure. ciRS-7 was identified as the first functional circRNA, which acted as a miRNA sponge to inhibit miR-7 activity in brain (6, 7). Since then, evidence has been obtained regarding expression patterns and functional roles of circRNAs in various biological processes including cancer progression. F-circSR1 and F-circSR2 promoted cell migration in lung cancer (8). CircIRAK3 facilitated breast cancer metastasis through sponging miR-3607 (9). circβ-catenin promoted cell proliferation by activating Wnt pathway in liver cancer (10). However, the importance of circRNAs in PDAC development is not fully understood.

In this study, we identified 169 differentially expressed circRNAs in PDAC cells by circRNA sequencing. We validated that one of the circRNAs, circFOXK2, was significantly upregulated in both PDAC cells and primary tumors. circFOXK2 promoted cell growth, clonogenic ability, migration, invasion, and liver metastasis in PDAC. Also, circFOXK2 functioned as a miRNA sponge for miR-942, and in turn promoted that expression of ankyrin 1 (ANK1), Glial cell–derived neurotrophic factor (GDNF), and paired box 6 (PAX6). Importantly, circRNA-pulldown and mass spectrometry identified 94 circFOXK2-interacting proteins to multiple biological processes including cell adhesion, mRNA splicing, and structural molecule activity. We demonstrated that circFOXK2 complexed with Y-box binding protein 1 (YBX1) and heterogeneous nuclear ribonucleoprotein K (hnRNPK) to promote the expression of oncogenic proteins NUF2 component of NDC80 kinetochore complex (NUF2) and pyridoxal kinase (PDXK) in PDAC.

Mammalian cell lines and clinical samples

PDAC cell lines PANC1, SW1990, CAPAN-2, CFPAC1, and BxPC-3, and HEK293 cells were obtained from ATCC. The human pancreatic ductal epithelial (HPDE) cell line was a gift from Dr. Ming-Sound Tsao (University Health Network, Ontario Cancer Institute and Princess Margaret Hospital Site, Toronto, Canada; ref. 11). All cell lines were verified by short tandem repeat profiling at the GENEWIZ, Inc. within 6 months of use and were cultured under the condition as described previously (12). All cell lines undergo routine Mycoplasma testing. 84 pairs of PDAC primary tumor and adjacent nontumor tissues were obtained from patients who underwent pancreatic resection at the Prince of Wales Hospital, Hong Kong. The study was carried out according to the ethical guidelines and with the approval of the Joint CUHK-NTEC Clinical Research Ethics Committee in accordance with Declaration of Helsinki and the written informed consent was obtained from all patients recruited. All specimens were fixed and embedded into paraffin.

circRNA sequencing and identification

circRNA sequencing was performed to analyze the expression pattern of circRNAs in HPDE, PANC1, and SW1990 cells. The total RNA extracted from cells was depleted from ribosomal RNA to linear RNA by RNase R. Then RNA was fragmented and was reverse transcribed. After linking with the sequencing adaptor and PCR amplification, a library was constructed for the circRNA sequencing. The raw sequencing reads from each sample were first mapped to the reference human genome using TopHat2 (13). The unmapped reads were extracted and mapped to reference human genome using TopHat-Fusion (14). Reads were processed into two anchors from both ends of the reads. Anchors that were aligned into the same chromosome but in reversed orientation were considered as the potential back-spliced junction reads. Because the common back-spliced junctions were GT/AG, GC/AG, and AT/AC, back-spliced junction reads with these junctions were extracted. Mapped reads (from TopHat2) and back-spliced reads (from TopHat-Fusion) were used to quantify the abundance of each circRNA candidate, denoted in RPM (reads per million mapped reads). DEGseq was used to compare the expression level of each circRNA between samples. circRNA sequencing data are available in the NCBI Gene Expression Omnibus under accession number GSE135731.

Plasmid and oligonucleotide transfection

circRNA overexpression plasmid was constructed by cloning the FOXK2 exon 2 and 3 into pcDNA3.1 (+) circRNA mini vector, which was a gift from Jeremy Wilusz (Addgene plasmid, catalog no. 60648;ref. 15). pmiR-Reporter plasmid with the inclusion of circRNA sequence in the 3′untranslated region for luciferase assay was constructed by cloning circRNA sequence into region directly downstream of the firefly luciferase gene in the pmiR-Reporter (Promega). Mutation of the miRNA binding site in the pmiR-Reporter plasmid was generated using KAPA HiFi DNA Polymerase (Kapa Biosystem) and primers with mutation sites. The lentiviral vector for stable knockdown of circFOXK2 was generated by cloning the short hairpin RNA (shRNA) sequence targeting circFOXK2 into the lentiviral vector with H1 promoter as described previously (16). miRNAs and siRNAs were purchased from GenePharma. Plasmids, miRNA mimics, and siRNA transfections were performed by lipofectamine 3000 (Invitrogen), according to the manufacturer's protocol. Sequences of miRNA mimics and siRNAs were listed in Supplementary Table S1.

qRT-PCR

The cytoplasmic and nuclear fractions were extracted using NE-PER Nuclear and Cytoplasmic Extraction Reagents (Thermo Fisher Scientific). RNA from whole-cell lysate or cell fractions was isolated by TRIzol Reagent. Formalin-fixed, paraffin-embedded (FFPE) sample RNA was isolated by miRNeasy FFPE Kit (Qiagen) according to the manufacturer's protocol. Measurement of gene expression level was performed by qRT-PCR. Reverse transcription of total RNA (except miRNA) was performed by High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems). Reverse transcription of miRNA was performed by Mir-X miRNA First-Strand Synthesis Kit (Takara). qRT-PCR was performed by ABI 7900HT Real-Time PCR system using SYBR Green PCR Master Mix (Applied Biosystems). The primers used in this study were listed in Supplementary Table S2.

In vivo subcutaneous injection

Male BALB/c nude mice ages 4 to 6 weeks were acquired from Laboratory Animal Services Centre of the Chinese University of Hong Kong (Shatin, Hong Kong). Animal handling and experimental procedures were approved by the Animal Experimental Ethics Committee of the institute. For tumor growth assay, 6 × 105 cells were resuspended in 1 × PBS with 20% Matrigel (Corning) and were injected subcutaneously into the right flank of the nude mice (7 mice/group). After tumor formation, tumor growth was monitored every 3–4 days, and the tumor volume was measured by a caliper and calculated by the equation: volume = (length × width2)/2.

In vivo orthotopic injection

Male BALB/c nude mice ages 4 to 6 weeks were acquired from Laboratory Animal Services Centre of the Chinese University of Hong Kong (Shatin, Hong Kong). Animal handling and experimental procedures were approved by the Animal Experimental Ethics Committee of the institute. For tumor metastasis assay, 5 × 105 cells were resuspended in 1 × PBS with 20% Matrigel and were injected orthotopically to the head of the pancreas (17). Tumors and organs were collected and examined for metastasis.

Luciferase assay

HEK293 cells were seeded in 96-well plate. After 24 hours, cells were cotransfected with 100 ng pmiR-Reporter plasmid, 1 ng Renilla luciferase reporter plasmid, and miRNA mimics. Each miRNA mimics or miRNC was cotransfected with pmiR-Reporter with or without circRNA sequence. After 72 hours, cells were washed twice with iced PBS and luciferase activity was measured with Nano-Glo Dual-Luciferase Reporter Assay System (Promega). For comparison, Firefly luciferase activity was first normalized with Renilla luciferase activity. Then, the effect of miRNAs on pmiR-Reporter with circRNA was normalized with pmiR-Reporter without circRNA. Finally, fold change was calculated by comparing the effect of miRNA mimics with miRNC.

miRNA pulldown

miRNA pulldown assay was performed by transfecting PANC1 cells with 100 nmol/L 3′-end biotinylated miRNA mimics. After 24 hours, the cells were washed twice with iced PBS, followed by cell lysis using miRNA pulldown lysis buffer (20 mmol/L Tris-HCl at pH 7.5, 100 mmol/L KCl, 5 nmol/L MgCl2, 0.3% NP-40). In the miRNA pulldown assay, 25 μL Dynabeads MyOne Streptavidin C1 (Thermo Fisher Scientific) was washed three times with miRNA pulldown lysis buffer. Then the beads were blocked with 1 mg/mL yeast tRNA (Thermo Fisher Scientific) and 1 mg/mL BSA at 4°C for 2 hours with rotation. After that, the beads were washed twice with miRNA pulldown lysis buffer. Biotin-labeled miRNAs were isolated by incubating the beads with 100 μL cell lysate and 100 μL miRNA pulldown lysis buffer at 4°C for 4 hours with rotation. The beads were washed twice with miRNA pulldown lysis buffer. Biotin-labeled miRNAs and their interacting RNAs were isolated by TRIzol Reagent. Detection of miRNA-interacting RNAs was performed by qRT-PCR.

In vitro transcription

circRNA overexpression plasmid was digested with XhoI (NEB) to form a linearized template DNA. Then in vitro transcription of circRNAs and their parental mRNAs was performed by MEGA script T7 Transcription Kit (Thermo Fisher Scientific; ref. 16). 1 μg template DNA was mixed with 1 μL of ATP solution, 1 μL of CTP solution, 1 μL of GTP solution, 0.9 μL of UTP solution, 0.15 μL of biotin-UTP solution (Epicentre), 1 μL of enzyme mixture, and RNase-free water to 10 μL. After in vitro transcription at 37°C overnight, template DNA was digested with DNase for 15 minutes at 37°C. For circRNA transcription, parental mRNA was digested by 10 U RNase R (Epicentre) for 5 hours at 37°C. Digestion was terminated by ammonium acetate stop solution and in vitro transcribed RNAs were purified by phenol:cholorform:isoamylalcohol solution according to the manufacturer's protocol.

RNA pulldown

RNA pulldown assay for investigating RNA–protein interaction was performed by Pierce Magnetic RNA-Protein Pull-Down Kit (Thermo Fisher Scientific). 5 μg of the in vitro transcripted RNA was heated at 90°C for 2 minutes in 1× RNA capture buffer, followed by incubation on ice for 2 minutes and at room temperature for 30 minutes. Then the biotin-labeled RNA was mixed with 75-μL streptavidin magnetic beads, which were prewashed twice with Tris buffer. After incubation at room temperature for 30 minutes, the RNA-labeled beads were washed twice with Tris buffer. 200 μg PANC1 protein was added to the RNA-labeled beads in 1× Protein-RNA Binding Buffer. The mixture was incubated at 4°C with rotation overnight. Then the RNA-labeled beads with proteins were washed twice with wash buffer. RNA-interacting proteins were eluted with 50-μL elution buffer by incubation at 37°C for 30 minutes.

Mass spectrometry

Proteins eluted from the RNA pulldown assay were resolved by 12% SDS-PAGE gel. When protein bands entered separating gel, electrophoresis was stopped so the bands were “stacked.” The protein bands were stained by the staining solution (0.1% Coomassie Blue R250 in 40% ethanol and 10% acetic acid) and were destained by the destaining solution. Then the gel was washed with the washing solution (50% methanol and 10% acetic acid) to remove the potential contaminants. The gel with bands was excised and was completely destained (JTBaker). In-gel digestion was performed by reduction and alkylation by 8 mmol/L DTT (USB Chemicals) and 40 mmol/L iodoacetamide (GE Healthcare), respectively. Protein digestion was performed by overnight trypsin (1 ng/μL; Promega) incubation at 37°C. Subsequent tryptic peptides were extracted from the gel with 100% ACN/3%TFA (JT Baker) and 40% ACN/3%TFA (JT Baker). The peptide extracts were pooled together and SpeedVac dried. Peptides derived from each condition were desalted using C18 StageTips (3M Corp) for LC/MS-MS analysis. Peptides were analyzed Dionex Ultimate3000 nanoRSLC system coupled to Thermo Fisher Orbitrap Fusion Tribid Lumos. Raw mass spectrometry data were processed using Proteome discoverer 2.1 version against the Human UniProt FASTA database (July 2017). Mass spectrometry data of circRNA pulldown are available via ProteomeXchange with identifier PXD015048.

RNA immunoprecipitation

RNA immunoprecipitation was performed by Magna RIP RNA Binding Protein Immunoprecipitation Kit (Millipore) according to the manufacturer's protocol. Briefly, cells were washed twice with ice-cold PBS, followed by cell lysis using equal volume of RIP lysis buffer. Magnetic beads were washed twice with RIP wash buffer, followed by incubation with 2 μg antibody against hnRNPK (rabbit; Proteintech 11426-1-AP), YBX1 (rabbit; Proteintech 20339-1-AP), SEPT11 (rabbit; Proteintech 14672-1-AP), ILF3 (rabbit; Proteintech 19887-1-AP), ASF (rabbit; Proteintech 12929-2-AP), and RAB11FIP1 (rabbit; Proteintech 16778-1-AP) for 30 minutes at room temperature. Immunoprecipitation was performed by incubating cell lysate with the magnetic bead–antibody complex overnight at 4°C. Then the beads were washed six times with RIP wash buffer, followed by proteinase K digestion at 37°C for 30 minutes. RNA was purified by phenol:chloroform:isoamylalcohol (15593-031, Invitrogen). qRT-PCR was used to analyze the enrichment of RNAs with target proteins.

Statistical analysis

Statistical analysis was performed by GraphPad Prism 7. Student t test, χ2t test and Pearson correlation were used as appropriate. Data were shown in mean ± SD. Statistically significant was considered when P value (two-sided) was less than 0.05.

Data Availability

Supplementary methods are described in Supplementary information.

Identification of circRNAs in PDAC cells

To identify the differentially expressed circRNAs in PDAC, we performed circRNA sequencing of nontumor HPDE and PDAC cells PANC1 and SW1990 (accession number GSE135731). In total, 17,158 circRNAs were identified in which 84% were exonic (Supplementary Fig. S1A). Among them, 83 upregulated and 86 downregulated circRNAs were revealed in PANC1 and SW1990 cells in comparison with nontumor HPDE cells (Supplementary Fig. S1B). To validate the results from circRNA sequencing and the clinical significances of the upregulated circRNAs, we measured the expressions of selected circRNAs in a panel of PDAC cells and primary tumors (Supplementary Fig. S1C and S1D). circFOXK2 was found with significant upregulation in PDAC cells (Fig. 1A), as well as 63% of PDAC primary tumors (53 of 84; Fig. 1B).

Figure 1.

Characterization of circFOXK2 in PDAC. A, The expression level of circFOXK2 in PDAC cells. The expression level was compared with nontumor HPDE cells. B, The expression level of circFOXK2 in PDAC primary tumor samples. The expression level of circFOXK2 in the tumor sample was compared with respective adjacent nontumor tissue. C, Genomic location of circFOXK2. circFOXK2 was formed by the back-splicing of exon 2 and 3 of FOXK2. D, PCR of gDNA and cDNA using divergent and convergent primers of circFOXK2. E, circFOXK2 was amplified using divergent primer, followed by Sanger sequencing of back-splicing junction. F, Expression level of circFOXK2 and GAPDH after RNase R treatment of total RNAs from PANC1 cells. G, circFOXK2 expression level after reverse transcription with random or oligo-dT primers using total RNAs from PANC1 cells. H, Expression level of circFOXK2 and FOXK2 after actinomycin D treatment for 0, 4, 12, and 24 hours in PANC1 cells. I, Cellular localization of circFOXK2 in PANC1 cells. Nuclear and cytoplasmic fraction was separated, followed by RNA extraction. circFOXK2, U6, and GAPDH levels were analyzed by qRT-PCR. Data are from at least three independent experiments. Mean ± SD. *, P < 0.05; **, P < 0.01; ***, P < 0.001.

Figure 1.

Characterization of circFOXK2 in PDAC. A, The expression level of circFOXK2 in PDAC cells. The expression level was compared with nontumor HPDE cells. B, The expression level of circFOXK2 in PDAC primary tumor samples. The expression level of circFOXK2 in the tumor sample was compared with respective adjacent nontumor tissue. C, Genomic location of circFOXK2. circFOXK2 was formed by the back-splicing of exon 2 and 3 of FOXK2. D, PCR of gDNA and cDNA using divergent and convergent primers of circFOXK2. E, circFOXK2 was amplified using divergent primer, followed by Sanger sequencing of back-splicing junction. F, Expression level of circFOXK2 and GAPDH after RNase R treatment of total RNAs from PANC1 cells. G, circFOXK2 expression level after reverse transcription with random or oligo-dT primers using total RNAs from PANC1 cells. H, Expression level of circFOXK2 and FOXK2 after actinomycin D treatment for 0, 4, 12, and 24 hours in PANC1 cells. I, Cellular localization of circFOXK2 in PANC1 cells. Nuclear and cytoplasmic fraction was separated, followed by RNA extraction. circFOXK2, U6, and GAPDH levels were analyzed by qRT-PCR. Data are from at least three independent experiments. Mean ± SD. *, P < 0.05; **, P < 0.01; ***, P < 0.001.

Close modal

Characterization of circFOXK2

Next, we examined the physical circular structure of circFOXK2. circFOXK2 is formed by the back-splicing of exon 2 and 3 of Forkhead box protein K2 (FOXK2; Fig. 1C). Outward-facing divergent primer and inward-facing convergent primer were designed to validate the formation of circFOXK2. Both primers amplified a product of expected size from cDNA, while only convergent primers resulted in an amplification from genomic DNA (gDNA; Fig. 1D). The presence of back-splicing junction was confirmed by Sanger sequencing (Fig. 1E). Also, circFOXK2 was observed to resist the digestion by RNase R, which specifically degraded linear RNAs but not circRNAs (Fig. 1F). Reduction in reverse-transcription efficiency by oligo-dT primers due to the lack of polyA tail also demonstrated the circularity of circFOXK2 (Fig. 1G). These results suggested the formation of circFOXK2 was not due to genomic rearrangement. Owing to its circular structure, we found that circFOXK2 was more stable than FOXK2 (Fig. 1H). Localization of circFOXK2 was examined by measuring its level in cytoplasm and nucleus, which demonstrated the enrichment of circFOXK2 in the cytoplasm of PDAC cells (Fig. 1I). Coding potential analysis suggested the lack of protein coding ability of circFOXK2 (Supplementary Table S3). Our results validated the circularity of circFOXK2.

circFOXK2 promotes cell growth and invasion in vitro

To elucidate the functions of circFOXK2 in PDAC, we employed siRNA, which specifically targeted the back-splicing junction of circFOXK2, without altering the expression of its parental gene (Supplementary Fig. S2A). Knockdown of circFOXK2 inhibited cell growth and clonogenic ability of PDAC cells (Fig. 2AC). The repressed cell growth by knockdown of circFOXK2 was due to cell-cycle arrest in G0–G1 phase and induction of apoptosis (Fig. 2D). Annexin V staining, cleavage of PARP and caspase-3, and decrease in Bcl2 level also demonstrated the induction of apoptosis after knockdown of circFOXK2 in PDAC cells (Fig. 2E and F). Also, knockdown of circFOXK2 inhibited cell migration and invasion (Fig. 2G and H). Furthermore, we constructed stable knockdown CFPAC1 cells using sh-circFOXK2 lentiviral system. Stable knockdown of circFOXK2 inhibited cell growth, clonogenic ability, and invasiveness (Supplementary Fig. S2B–S2E). On the other hand, overexpression of circFOXK2 in HPDE cells significantly promoted cell growth, migration, and invasion. (Supplementary Fig. S3A–S3E). Collectively, our results demonstrated the importance of circFOXK2 in promoting PDAC cell growth and invasion.

Figure 2.

circFOXK2 promotes cell growth, migration, and invasion in PDAC cells. A, Cell growth analyzed by MTT assay after knockdown of circFOXK2 by siRNA for 24, 72, and 120 hours in PDAC cells. B, Anchorage-dependent colony formation after knockdown of circFOXK2 in PDAC cells. Cells were stained with crystal violet. C, Anchorage-independent colony formation after knockdown of circFOXK2 in PANC1 cells. Cells were stained with crystal violet. D, Analysis of cell-cycle distribution by flow cytometry after knockdown of circFOXK2 in CFPAC1 (left) and PANC1 (right) cells for 120 hours. E, Western blot analysis of apoptosis marker PARP, Bcl2, and caspase-3 after knockdown of circFOXK2 in CFPAC1 cells. F, Cell apoptosis assessed by Annexin V staining (red) after knockdown of circFOXK2 in PADC cells for 72 hours. G and H, Wound healing cell migration assay (G) and Transwell cell invasion assay (H) after knockdown of circFOXK2 in PDAC cells. Cells in invasion assay were stained by crystal violet. Data are from at least three independent experiments. Mean ± SD. *, P < 0.05; **, P < 0.01; ***, P < 0.001.

Figure 2.

circFOXK2 promotes cell growth, migration, and invasion in PDAC cells. A, Cell growth analyzed by MTT assay after knockdown of circFOXK2 by siRNA for 24, 72, and 120 hours in PDAC cells. B, Anchorage-dependent colony formation after knockdown of circFOXK2 in PDAC cells. Cells were stained with crystal violet. C, Anchorage-independent colony formation after knockdown of circFOXK2 in PANC1 cells. Cells were stained with crystal violet. D, Analysis of cell-cycle distribution by flow cytometry after knockdown of circFOXK2 in CFPAC1 (left) and PANC1 (right) cells for 120 hours. E, Western blot analysis of apoptosis marker PARP, Bcl2, and caspase-3 after knockdown of circFOXK2 in CFPAC1 cells. F, Cell apoptosis assessed by Annexin V staining (red) after knockdown of circFOXK2 in PADC cells for 72 hours. G and H, Wound healing cell migration assay (G) and Transwell cell invasion assay (H) after knockdown of circFOXK2 in PDAC cells. Cells in invasion assay were stained by crystal violet. Data are from at least three independent experiments. Mean ± SD. *, P < 0.05; **, P < 0.01; ***, P < 0.001.

Close modal

circFOXK2 promotes tumor growth and metastasis in vivo

Knockdown of circFOXK2 was also employed for in vivo study. Xenograft mice model was generated by subcutaneous injection of sh-circFOXK2 CFPAC1 cells. Knockdown of circFOXK2 by shRNA significantly inhibited tumor growth (Fig. 3AE). Because results from in vitro experiments revealed that circFOXK2 promoted PDAC cell migration and invasion, we constructed PDAC metastatic mice model through orthotopic injection of sh-circFOXK2 CFPAC1 cells to the pancreas. Knockdown of circFOXK2 inhibited the PDAC metastasis to liver (Fig. 3F). Our results found that that circFOXK2 played critical roles in promoting PDAC growth and liver metastasis.

Figure 3.

circFOXK2 promotes tumor growth and liver metastasis in PDAC cells. A, Photographs of mice xenograft at day 34 after knockdown of circFOXK2 in CFPAC1 cells. B, Tumor volume of xenograft after knockdown of circFOXK2. C, Tumor mass of xenograft after knockdown of circFOXK2. D, Hematoxylin and eosin staining of mice subcutaneous tumors after knockdown of circFOXK2. E, IHC staining of Ki67 and Bcl2 in mice subcutaneous tumors after knockdown of circFOXK2. F, Analysis of tumor metastasis to liver after knockdown of circFOXK2 in CFPAC1 cells. Data are from at least three independent experiments. Mean ± SD. *, P < 0.05; **, P < 0.01; ***, P < 0.001.

Figure 3.

circFOXK2 promotes tumor growth and liver metastasis in PDAC cells. A, Photographs of mice xenograft at day 34 after knockdown of circFOXK2 in CFPAC1 cells. B, Tumor volume of xenograft after knockdown of circFOXK2. C, Tumor mass of xenograft after knockdown of circFOXK2. D, Hematoxylin and eosin staining of mice subcutaneous tumors after knockdown of circFOXK2. E, IHC staining of Ki67 and Bcl2 in mice subcutaneous tumors after knockdown of circFOXK2. F, Analysis of tumor metastasis to liver after knockdown of circFOXK2 in CFPAC1 cells. Data are from at least three independent experiments. Mean ± SD. *, P < 0.05; **, P < 0.01; ***, P < 0.001.

Close modal

Identification of circFOXK2-interacting miRNAs

We next investigated the detailed mechanism of circFOXK2 in PDAC. Because many cytoplasmic circRNAs serve as “microRNA sponges” in regulating cancer progression, including liver, lung, and breast cancer (8–10), we examined whether circFOXK2 could function as a miRNA sponge. First, we performed the bioinformatic analysis by TargetScan (18) to reveal potential miRNA binding sites on circFOXK2. To validate the circFOXK2–miRNA interaction, luciferase reporter assay using circFOXK2-luciferase reporter and miRNA mimics was performed (Fig. 4A), which demonstrated a reduction in luciferase activity with miR-942 mimics (Fig. 4B). Mutating the miR-942 binding sites restored the luciferase activity (Fig. 4C). circFOXK2–miRNA interaction was further validated by miRNA pulldown. Biotin-labeled miR-942 mimics could significantly enriched circFOXK2 in PANC1 cells (Fig. 4D). These results suggested the interaction between circFOXK2 and miR-942 in PDAC cells.

Figure 4.

circFOXK2 functions as miR-942 sponge. A, Bioinformatic analysis of potential miRNA binding sites on circFOXK2. Potential binding sites were mutated for luciferase assay. B, Luciferase reporter assay for luciferase activity of circFOXK2-reporter in HEK293 cells cotransfected with miRNA mimics. C, Luciferase reporter assay for luciferase activity of circFOXK2-reporter with mutated miR-942 binding sites in HEK293 cells cotransfected with miR-942 mimics. D, miRNAs pulldown assay in PANC1 cells cotransfected with biotin-labeled miR-942 mimics, which was analyzed by qRT-PCR. E, miR-942 expression level after knockdown of circFOXK2 in PANC1 cells. F, circFOXK2 expression level after transfecting miR-942 mimics in PANC1 cells. Data are from at least three independent experiments. Mean ± SD. **, P < 0.01; ***, P < 0.001.

Figure 4.

circFOXK2 functions as miR-942 sponge. A, Bioinformatic analysis of potential miRNA binding sites on circFOXK2. Potential binding sites were mutated for luciferase assay. B, Luciferase reporter assay for luciferase activity of circFOXK2-reporter in HEK293 cells cotransfected with miRNA mimics. C, Luciferase reporter assay for luciferase activity of circFOXK2-reporter with mutated miR-942 binding sites in HEK293 cells cotransfected with miR-942 mimics. D, miRNAs pulldown assay in PANC1 cells cotransfected with biotin-labeled miR-942 mimics, which was analyzed by qRT-PCR. E, miR-942 expression level after knockdown of circFOXK2 in PANC1 cells. F, circFOXK2 expression level after transfecting miR-942 mimics in PANC1 cells. Data are from at least three independent experiments. Mean ± SD. **, P < 0.01; ***, P < 0.001.

Close modal

To study the functions of circFOXK2–miR-942 interaction in PDAC, we first knocked down circFOXK2, which leads to a subsequent increase in miR-942 expression level (Fig. 4E). Conversely, transfecting miR-942 mimics did not alter the expression of circFOXK2, but downregulated the targets of miR-942: ANK1, GDNF, and PAX6 (Fig. 4F; Supplementary Fig. S4A and S4B). These results suggested that circFOXK2 inhibited the activity of miR-942 in PDAC. As the role of miR-942 in PDAC is unclear, we attempted to study the importance of miR-942 in PDAC. We found that miR-942 was frequently downregulated in PDAC cells and primary tumors (Supplementary Fig. S5A–S5C). Also, we found that transfecting miR-942 mimics inhibited cell growth and clonogenic ability in PDAC cells (Supplementary Fig. S5D and S5E). Our results illustrated the role of miR-942 in promoting PDAC cell growth.

Because circFOXK2 functioned as a miRNA sponge for miR-942, the expression level of ANK1, GDNF, and PAX6 were measured after knockdown or overexpression of circFOXK2. Knockdown of circFOXK2 decreased the expression levels of ANK1, GDNF, and PAX6 in PDAC cells, while overexpression of circFOXK2 promoted the expression of ANK1, GDNF, and PAX6 in HPDE cells (Fig. 5A and B; Supplementary Fig. S6A–S6D). The effect of circFOXK2 overexpression on GDNF and PAX6 could be rescued by the introduction of miR-942 mimics (Fig. 5C). Also, in the mice tumors with circFOXK2 knocked down, ANK1, GDNF, and PAX6 were downregulated (Fig. 5D; Supplementary Fig. S6E). Analysis of PDAC primary tumors also revealed the positive correlation between circFOXK2, ANK1, and PAX6 levels (Fig. 5E). Collectively, our results suggested that circFOXK2 inhibited miR-942, and in turn promoted the expression of ANK1, GDNF, and PAX6 in PDAC.

Figure 5.

circFOXK2 promotes the expression of miR-942 targets in PDAC. A and B, qRT-PCR analysis of mRNA and Western blot analysis of protein levels of ANK1, GDNF, and PAX6 after knockdown of circFOXK2 in PDAC cells. C, GDNF and PAX6 expression levels in circFOXK2-overexpressing HPDE cells after transfecting miR-942 mimics. D, IHC staining of ANK1, GDNF, and PAX6 in mice subcutaneous tumors after knockdown of circFOXK2. E, Correlation between circFOXK2 and ANK1, and PAX6 levels in PDAC primary tumors. Data are from at least three independent experiments. Mean ± SD. *, P < 0.05; **, P < 0.01; ***, P < 0.001.

Figure 5.

circFOXK2 promotes the expression of miR-942 targets in PDAC. A and B, qRT-PCR analysis of mRNA and Western blot analysis of protein levels of ANK1, GDNF, and PAX6 after knockdown of circFOXK2 in PDAC cells. C, GDNF and PAX6 expression levels in circFOXK2-overexpressing HPDE cells after transfecting miR-942 mimics. D, IHC staining of ANK1, GDNF, and PAX6 in mice subcutaneous tumors after knockdown of circFOXK2. E, Correlation between circFOXK2 and ANK1, and PAX6 levels in PDAC primary tumors. Data are from at least three independent experiments. Mean ± SD. *, P < 0.05; **, P < 0.01; ***, P < 0.001.

Close modal

Identification of circFOXK2-interacting proteins

Many studies have demonstrated the major role of circRNA through functioning as miRNA sponges, whereas the binding capability of circRNA on protein in regulating gene expression remains uncertain. Therefore, we next studied the importance of circRNA–protein interaction in PDAC. We in vitro transcribed circFOXK2, which functioned as a probe to pull-down circFOXK2-interacting proteins. To eliminate the “false positive” of the pull-down proteins, we constructed additional two negative controls: (i) RNA probe with an identical sequence as circFOXK2, without the formation of circRNAs; (ii) RNA probe with similar secondary structure as circRNA, which was circGFP (Fig. 6A). Therefore, the proteins pulled down by circFOXK2 were due to the unique sequence and secondary structure of circFOXK2. Mass spectrometry analysis identified 94 circFOXK2-interacting proteins in PANC1 cells (accession number PXD015048; Fig. 6B). Also, protein–protein interactions (PPI) of the circRNA-interacting proteins were analyzed in the STRING database (19). 101 PPIs were observed for circFOXK2-interacting proteins (Supplementary Fig. S7A). Gene ontology analysis also revealed that the circFOXK2-interacting proteins play important roles in many biological processes, including cell adhesion, mRNA splicing, and structural molecule activity (Supplementary Fig. S7B). Our results suggested that circFOXK2 interacted with multiple proteins with critical roles in PDAC.

Figure 6.

Identification of circFOXK2-interacting proteins in PDAC. A, Schematic diagram showing the process of circFOXK2-pulldown. circFOXK2, FOXK2, and circGFP probes were generated by in vitro transcription with biotin-UTP. The biotin-labeled probes were used to pull down interacting proteins in PANC1 cells. The pulldown proteins were identified by mass spectrometry. B, Mass spectrometry analysis of circFOXK2-interacting proteins in PANC1 cells. C, RNA immunoprecipitation for YBX1, hnRNPK, SEPT11, ILF3, RAB11FIP1, and ASF in PDAC cells. The interactions between circFOXK2 and YBX1, hnRNPK, SEPT11, ILF3, RABIIFIP1, and ASF were analyzed by qRT-PCR. D, Expression levels of YBX1 and hnRNPK targets CORO1C, NUF2, PDXK, and PPP2R2C after knockdown of circFOXK2 in PDAC cells. E, circFOXK2-pulldown assay in PANC1 cells using biotin-labeled circFOXK2 probe. The interaction between circFOXK2, NUF2, and PDXK was analyzed by qRT-PCR. F, The expressions of NUF2 and PDXK in circFOXK2-overexpressed HPDE cells after knockdown of YBX1 and hnRNPK. G and H, RNA immunoprecipitation for YBX1 (G) and hnRNPK (H) after knockdown of circFOXK2 in PANC1 cells. The NUF2 and PDXK enrichments were analyzed by qRT-PCR. I, Expression levels of PDXK and NUF2 and the correlations between NUF2 and circFOXK2 level in PDAC primary tumors. Data are from at least three independent experiments. Mean ± SD. *, P < 0.05; **, P < 0.01; ***, P < 0.001.

Figure 6.

Identification of circFOXK2-interacting proteins in PDAC. A, Schematic diagram showing the process of circFOXK2-pulldown. circFOXK2, FOXK2, and circGFP probes were generated by in vitro transcription with biotin-UTP. The biotin-labeled probes were used to pull down interacting proteins in PANC1 cells. The pulldown proteins were identified by mass spectrometry. B, Mass spectrometry analysis of circFOXK2-interacting proteins in PANC1 cells. C, RNA immunoprecipitation for YBX1, hnRNPK, SEPT11, ILF3, RAB11FIP1, and ASF in PDAC cells. The interactions between circFOXK2 and YBX1, hnRNPK, SEPT11, ILF3, RABIIFIP1, and ASF were analyzed by qRT-PCR. D, Expression levels of YBX1 and hnRNPK targets CORO1C, NUF2, PDXK, and PPP2R2C after knockdown of circFOXK2 in PDAC cells. E, circFOXK2-pulldown assay in PANC1 cells using biotin-labeled circFOXK2 probe. The interaction between circFOXK2, NUF2, and PDXK was analyzed by qRT-PCR. F, The expressions of NUF2 and PDXK in circFOXK2-overexpressed HPDE cells after knockdown of YBX1 and hnRNPK. G and H, RNA immunoprecipitation for YBX1 (G) and hnRNPK (H) after knockdown of circFOXK2 in PANC1 cells. The NUF2 and PDXK enrichments were analyzed by qRT-PCR. I, Expression levels of PDXK and NUF2 and the correlations between NUF2 and circFOXK2 level in PDAC primary tumors. Data are from at least three independent experiments. Mean ± SD. *, P < 0.05; **, P < 0.01; ***, P < 0.001.

Close modal

circFOXK2 promotes the expression of oncogenic proteins via complexing YBX1 and hnRNPK

RNA immunoprecipitation assay was performed in PDAC cells to validate the results from mass spectrometry analysis. circFOXK2 was significantly enriched by YBX1 and hnRNPK (Fig. 6C). Because YBX1 complexed with hnRNPK in regulating gene expression to promote PDAC progression (20; Supplementary Fig. S8), we hypothesized that circFOXK2 interacted with YBX1 and hnRNPK complex to promote PDAC development. Publicly available datasets of knockdown of YBX1 were first employed to identify genes that were consistently regulated by YBX1 in multiple cancers (Supplementary Fig. S9A). RNA immunoprecipitation assay with YBX1 and hnRNPK demonstrated the enrichment of Coronin 1C (CORO1C), NUF2, PDXK, and Protein Phosphatase 2 Regulatory Subunit B gamma (PPP2R2C), indicating they were the direct targets of YBX1 and hnRNPK complex (Supplementary Fig. S9B). Consistently, knockdown of YBX1 significantly downregulated these targets in CFPAC1 cells (Supplementary Fig. S9C). Analysis of The Cancer Genome Atlas dataset also revealed the upregulation of PPP2R2C, PDXK, and NUF2, and their positive correlations to YBX1 expression in PDAC (Supplementary Fig. S9D–S9F). Importantly, NUF2 and PDXK were downregulated after the knockdown of circFOXK2 while they were significantly enriched by circFOXK2 pulldown (Fig. 6D and E). The effect of circFOXK2 overexpression on NUF2 and PDXK expressions could be rescued by the knockdown of YBX1 and hnRNPK (Fig. 6F). Knockdown of circFOXK2 decreased the interaction of YBX1 and hnRNPK to NUF2 and PDXK (Fig. 6G and H). Analysis of PDAC cells and primary tumors revealed the upregulation of PDXK and NUF2, and the positive correlation between NUF2 and circFOXK2 expression (Fig. 6I; Supplementary Fig. S10A and S10B). Also, knockdown of NUF and PDXK significantly inhibited PDAC cell growth and invasion (Supplementary Fig. S10C–S10E). Collectively, our results suggested that circFOXK2, YBX1, and hnRNPK complex interacted and promoted the expressions of oncogenic proteins NUF2 and PDXK in PDAC.

In this study, we profiled circRNAs expressions in nontumor HPDE and PDAC cells by circRNA sequencing. We identified 169 differentially expressed circRNAs in PDAC. Further validation experiments demonstrated the upregulation of circFOXK2 in PDAC cells and primary tumors. Loss-of-function by both siRNA and shRNA and gain-of-function studies suggested the critical roles of circFOXK2 in PDAC progression by enhancing cell growth, migration, invasion, and metastasis in vitro and in vivo. Mechanistically, circFOXK2 functioned as a miRNA sponge and interacted with multiple proteins to regulate gene expressions.

Many studies have illustrated the importance of circRNAs in cancer progression (8–10). Previous circRNA profiling in 20 PDAC tissues identified 289 dysregulated circRNAs (21). Another study using 6 PDAC tissues discovered 351 differentially expressed circRNAs in PDAC (22). Recently, analysis using publicly available circRNA and miRNA datasets identified 256 differentially expressed circRNAs with 51 circRNA–miRNA interactions in PDAC (23). In general, these differentially expressed circRNAs were involved in protein localization, RNA binding, chromatin modification, and protein binding (21). circRNA_100782 functioned as a miR-124 sponge to promote PDAC cell proliferation via IL6-STAT3 pathway (24). hsa_circ_0006215 promoted migration via inhibiting miR-378a-3p in BxPC-3 cells (25). circ-PDE8A stimulated MACC1 expression by sponging miR-338 to promote cell invasion (26). However, the detailed roles of circRNAs in PDAC are still very unclear. Therefore, our study is the first to perform circRNA sequencing to identify differentially expressed circRNAs in PDAC. We identified 83 upregulated and 86 downregulated circRNAs in PDAC cells. Also, we demonstrated upregulation of circFOXK2 in PDAC cells and primary tumors played important roles in promoting cancer cell growth, migration, invasion, and liver metastasis.

The functional role of circRNAs was first described as miRNA sponges (6). CDR1as harbored 63 binding sites for miR-77. The interaction between CDR1as and miR-7 blocked the miR-7 activity. Since that, the functions of circRNAs as a “miRNA sponges” were comprehensively studied in many biological processes. Also, the importance of circRNA–miRNA interactions was also documented in many cancers, including PDAC. circLARP4 inhibited miR-424-5p and regulated LASTS1 expression in gastric cancer (27). circITCH absorbed miR-17 and miR-224 to promote p21 and PTEN expression in bladder cancer (28). ciRS-7 targeted miR-7 to promote EGFR/STAT3 signaling, while circZMYM2 sponged miR335-5p and in turn favored the expression of JMJD2C in PDAC (29, 30). Here, we performed the bioinformatic analysis, luciferase assay, and miRNA pulldown assay to identify circFOXK2-interacting miRNAs in PDAC cells. We observed that circFOXK2 absorbed miR-942, and in turn promoted the expression of ANK1, GDNF, and PAX6. Studies indicated the potential of miR-942 as a diagnostic marker for PDAC (31), nevertheless, the roles of miR-942 in PDAC is poorly understood. We found that miR-942 inhibited PDAC cell growth and was the negative regulator of ANK1, GDNF, and PAX6, which were involved in PDAC cell growth, invasion, and metastasis (32–34).

Apart from functioning as miRNA sponge, there remained limited evidence on circRNA–protein interactions. circ-Amotl1 interacted with c-myc and promoted c-myc translocation to the nucleus in breast cancer (35). Circ-CTNNB1 interacted with DDX3 to transactivate YY1 and YY1 downstream targets in gastric cancer (36). Importantly, circRNA–protein interaction has not been documented in PDAC. Here, we developed a novel and highly specific circRNA-pulldown assay and identified the interacting proteins of circFOXK2 by mass spectrometry. We identified 94 circFOXK2-interacting proteins in PANC1 cells. Also, analysis using STRING database revealed that the circFOXK2-interacting proteins formed complexes with 101 protein–protein interactions. Gene ontology analysis indicated these proteins were associated with multiple biological processes, including cell adhesion, structural molecule activity, and mRNA splicing. These may suggest that circFOXK2 directly or indirectly interacts with multiple proteins with important roles in biological processes.

We then validated the results from mass spectrometry by RNA immunoprecipitation assay. Notably, we found that YBX1 and hnRNPK interacted with circFOXK2 in both PANC1 and CFPAC1 cells. YBX1 functions as RNA binding protein that regulates transcription, mRNA processing, and translation. YBX1 was frequently upregulated and involved in cell invasion and drug-resistance development in many cancers (37–39). In PDAC, studies reported that YBX1 promoted cell invasion (40). hnRNPK, as a partner of YBX1 in gene regulation, was also aberrantly expressed in various cancers (41). hnRNPK regulated cell migration via promoting the expression of GSN in lung adenocarcinoma (42). The upregulation of hnRNPK also favored cancer metastasis by promoting the expression of MMP-3 (43). Also, YBX1/hnRNPK complex was involved in cancer progression by regulating the splicing process (44). In this study, we found that YBX1 and hnRNPK promoted PDAC cell invasion. Also, we found that CORO1C, NUF2, PDXK, and PPP2R2C were directly regulated by YBX1 and hnRNPK complex in PDAC cells. Reported studies have demonstrated the growth-promoting role of NUF2 in cancer, including PDAC (45). Also, we observed the upregulation of NUF2 and PDXK in promoting PDAC cell growth and invasion. Importantly, among these YBX1 and hnRNPK targets, NUF2 and PDXK were also regulated by circFOXK2. Notably, knockdown of circFOXK2 reduced the interaction of YBX1 and hnRNPK to NUF2 and PDXK, suggesting circFOXK2 played an important role in YBX1- and hnRNPK-mediated gene regulation. Collectively, circFOXK2 complexed with YBX1 and hnRNPK to promote the expression of oncogenic proteins NUF2 and PDXK.

In conclusion, we profiled the circRNAs expressions in HPDE, PANC1, and SW1990 cells by circRNA sequencing. circFOXK2 was significantly upregulated in both PDAC cells and primary tumors. Gain-of-function and loss-of-function studies demonstrated important roles of circFOXK2 in promoting cancer cell growth, migration, invasion, and liver metastasis. We found that circFOXK2 functioned as a sponge of miR-942, and in turn promoted the expressions of ANK1, GDNF, and PAX6 (Fig. 7). More importantly, circFOXK2 interacted with YBX1 and hnRNPK to promote the expressions of oncogenic proteins NUF2 and PDXK (Fig. 7). These revealed a novel mechanism of circRNA in regulating gene expressions by interacting with YBX1/hnRNPK complex.

Figure 7.

Schematic diagram describing the role of circFOXK2 in the development and progression of PDAC. circFOXK2 functioned as a sponge for miR-942 and in turn promoted the expression of ANK1, GDNF, and PAX6. Also, circFOXK2 interacted with YBX1 and hnRNPK to promote the expression of NUF2 and PDXK. These pathways favored cell growth, migration, invasion, and liver metastasis in PDAC.

Figure 7.

Schematic diagram describing the role of circFOXK2 in the development and progression of PDAC. circFOXK2 functioned as a sponge for miR-942 and in turn promoted the expression of ANK1, GDNF, and PAX6. Also, circFOXK2 interacted with YBX1 and hnRNPK to promote the expression of NUF2 and PDXK. These pathways favored cell growth, migration, invasion, and liver metastasis in PDAC.

Close modal

No potential conflicts of interest were disclosed.

Conception and design: C.H. Wong, Y. Chen

Development of methodology: C.H. Wong, Y. Li

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): C.H. Wong, U.K. Lou, Y. Li, S.L. Chan, J.H.M. Tong, K.-F. To

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): C.H. Wong, Y. Li, Y. Chen

Writing, review, and/or revision of the manuscript: C.H. Wong, S.L. Chan, K.-F. To, Y. Chen

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): C.H. Wong, Y. Chen

Study supervision: Y. Chen

This work was supported by General Research Fund, Research Grants Council of Hong Kong (4171217 and 14120618); National Natural Science Foundation of China (81672323); and Direct Grant from CUHK to Y. Chen.

Mass spectrometry was performed by The Proteomics and Metabolomics Core Facility and Li Ka Shing Faculty of Medicine, The University of Hong Kong.

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Siegel
RL
,
Miller
KD
,
Jemal
A
. 
Cancer statistics, 2018
.
CA Cancer J Clin
2018
;
68
:
7
30
.
2.
Malvezzi
M
,
Bertuccio
P
,
Rosso
T
,
Rota
M
,
Levi
F
,
La Vecchia
C
, et al
European cancer mortality predictions for the year 2015: does lung cancer have the highest death rate in EU women?
Ann Oncol
2015
;
26
:
779
86
.
3.
Cress
RD
,
Yin
D
,
Clarke
L
,
Bold
R
,
Holly
EA
. 
Survival among patients with adenocarcinoma of the pancreas: a population-based study (United States)
.
Cancer Causes Control
2006
;
17
:
403
9
.
4.
Springfeld
C
,
Jäger
D
,
Büchler
MW
,
Strobel
O
,
Hackert
T
,
Palmer
DH
, et al
Chemotherapy for pancreatic cancer
.
Presse Med
2019
;
48
:
e159
74
.
5.
Von Hoff
DD
,
Ervin
T
,
Arena
FP
,
Chiorean
EG
,
Infante
J
,
Moore
M
, et al
Increased survival in pancreatic cancer with nab-paclitaxel plus gemcitabine
.
N Engl J Med
2013
;
369
:
1691
703
.
6.
Memczak
S
,
Jens
M
,
Elefsinioti
A
,
Torti
F
,
Krueger
J
,
Rybak
A
, et al
Circular RNAs are a large class of animal RNAs with regulatory potency
.
Nature
2013
;
495
:
333
8
.
7.
Hansen
TB
,
Jensen
TI
,
Clausen
BH
,
Bramsen
JB
,
Finsen
B
,
Damgaard
CK
, et al
Natural RNA circles function as efficient microRNA sponges
.
Nature
2013
;
495
:
384
8
.
8.
Wu
K
,
Liao
X
,
Gong
Y
,
He
J
,
Zhou
JK
,
Tan
S
, et al
Circular RNA F-circSR derived from SLC34A2-ROS1 fusion gene promotes cell migration in non-small cell lung cancer
.
Mol Cancer
2019
;
18
:
98
.
9.
Wu
J
,
Jiang
Z
,
Chen
C
,
Hu
Q
,
Fu
Z
,
Chen
J
, et al
CircIRAK3 sponges miR-3607 to facilitate breast cancer metastasis
.
Cancer Lett
2018
;
430
:
179
92
.
10.
Liang
WC
,
Wong
CW
,
Liang
PP
,
Shi
M
,
Cao
Y
,
Rao
ST
, et al
Translation of the circular RNA circβ-catenin promotes liver cancer cell growth through activation of the Wnt pathway
.
Genome Biol
2019
;
20
:
84
.
11.
Ouyang
H
,
Mou
Lj
,
Luk
C
,
Liu
N
,
Karaskova
J
,
Squire
J
, et al
Immortal human pancreatic duct epithelial cell lines with near normal genotype and phenotype
.
Am J Pathol
2000
;
157
:
1623
31
.
12.
Li
CH
,
Xiao
Z
,
Tong
JH
,
To
KF
,
Fang
X
,
Cheng
AS
, et al
EZH2 coupled with HOTAIR to silence MicroRNA-34a by the induction of heterochromatin formation in human pancreatic ductal adenocarcinoma
.
Int J Cancer
2017
;
140
:
120
9
.
13.
Kim
D
,
Pertea
G
,
Trapnell
C
,
Pimentel
H
,
Kelley
R
,
Salzberg
SL
. 
TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions
.
Genome Biol
2013
;
14
:
R36
.
14.
Kim
D
,
Salzberg
SL
. 
TopHat-Fusion: an algorithm for discovery of novel fusion transcripts
.
Genome Biol
2011
;
12
:
R72
.
15.
Liang
D
,
Wilusz
JE
. 
Short intronic repeat sequences facilitate circular RNA production
.
Genes Dev
2014
;
28
:
2233
47
.
16.
Xu
F
,
Li
CH
,
Wong
CH
,
Chen
GG
,
Lai
PBS
,
Shao
S
, et al
Genome-wide screening and functional analysis identifies tumor suppressor long noncoding RNAs epigenetically silenced in hepatocellular carcinoma
.
Cancer Res
2019
;
79
:
1305
17
.
17.
Li
CH
,
To
KF
,
Tong
JH
,
Xiao
Z
,
Xia
T
,
Lai
PB
, et al
Enhancer of zeste homolog 2 silences microRNA-218 in human pancreatic ductal adenocarcinoma cells by inducing formation of heterochromatin
.
Gastroenterology
2013
;
144
:
1086
97
.
18.
Grimson
A
,
Farh
KK
,
Johnston
WK
,
Garrett-Engele
P
,
Lim
LP
,
Bartel
DP
. 
MicroRNA targeting specificity in mammals: determinants beyond seed pairing
.
Mol Cell
2007
;
27
:
91
105
.
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
.
20.
Shnyreva
M
,
Schullery
DS
,
Suzuki
H
,
Higaki
Y
,
Bomsztyk
K
. 
Interaction of two multifunctional proteins. Heterogeneous nuclear ribonucleoprotein K and Y-box-binding protein
.
J Biol Chem
2000
;
275
:
15498
503
.
21.
Guo
S
,
Xu
X
,
Ouyang
Y
,
Wang
Y
,
Yang
J
,
Yin
L
, et al
Microarray expression profile analysis of circular RNAs in pancreatic cancer
.
Mol Med Rep
2018
;
17
:
7661
71
.
22.
Li
H
,
Hao
X
,
Wang
H
,
Liu
Z
,
He
Y
,
Pu
M
, et al
Circular RNA expression profile of pancreatic ductal adenocarcinoma revealed by microarray
.
Cell Physiol Biochem
2016
;
40
:
1334
44
.
23.
Zhang
Q
,
Wang
JY
,
Zhou
SY
,
Yang
SJ
,
Zhong
SL
. 
Circular RNA expression in pancreatic ductal adenocarcinoma
.
Oncol Lett
2019
;
18
:
2923
30
.
24.
Chen
G
,
Shi
Y
,
Zhang
Y
,
Sun
J
. 
CircRNA_100782 regulates pancreatic carcinoma proliferation through the IL6-STAT3 pathway
.
Onco Targets Ther
2017
;
10
:
5783
94
.
25.
Zhu
P
,
Ge
N
,
Liu
D
,
Yang
F
,
Zhang
K
,
Guo
J
, et al
Preliminary investigation of the function of hsa_circ_0006215 in pancreatic cancer
.
Oncol Lett
2018
;
16
:
603
11
.
26.
Li
Z
,
Yanfang
W
,
Li
J
,
Jiang
P
,
Peng
T
,
Chen
K
, et al
Tumor-released exosomal circular RNA PDE8A promotes invasive growth via the miR-338/MACC1/MET pathway in pancreatic cancer
.
Cancer Lett
2018
;
432
:
237
50
.
27.
Zhang
J
,
Liu
H
,
Hou
L
,
Wang
G
,
Zhang
R
,
Huang
Y
, et al
Circular RNA_LARP4 inhibits cell proliferation and invasion of gastric cancer by sponging miR-424-5p and regulating LATS1 expression
.
Mol Cancer
2017
;
16
:
151
.
28.
Yang
C
,
Yuan
W
,
Yang
X
,
Li
P
,
Wang
J
,
Han
J
, et al
Circular RNA circ-ITCH inhibits bladder cancer progression by sponging miR-17/miR-224 and regulating p21, PTEN expression
.
Mol Cancer
2018
;
17
:
19
.
29.
Liu
L
,
Liu
FB
,
Huang
M
,
Xie
K
,
Xie
QS
,
Liu
CH
, et al
Circular RNA ciRS-7 promotes the proliferation and metastasis of pancreatic cancer by regulating miR-7-mediated EGFR/STAT3 signaling pathway
.
Hepatobiliary Pancreat Dis Int
2019
;
18
:
580
6
.
30.
An
Y
,
Cai
H
,
Zhang
Y
,
Liu
S
,
Duan
Y
,
Sun
D
, et al
circZMYM2 competed endogenously with miR-335-5p to regulate JMJD2C in pancreatic cancer
.
Cell Physiol Biochem
2018
;
51
:
2224
36
.
31.
Cao
Z
,
Liu
C
,
Xu
J
,
You
L
,
Wang
C
,
Lou
W
, et al
Plasma microRNA panels to diagnose pancreatic cancer: results from a multicenter study
.
Oncotarget
2016
;
7
:
41575
83
.
32.
Omura
N
,
Mizuma
M
,
MacGregor
A
,
Hong
SM
,
Ayars
M
,
Almario
JA
, et al
Overexpression of ankyrin1 promotes pancreatic cancer cell growth
.
Oncotarget
2016
;
7
:
34977
87
.
33.
Okada
Y
,
Eibl
G
,
Duffy
JP
,
Reber
HA
,
Hines
OJ
. 
Glial cell-derived neurotrophic factor upregulates the expression and activation of matrix metalloproteinase-9 in human pancreatic cancer
.
Surgery
2003
;
134
:
293
9
.
34.
Mascarenhas
JB
,
Young
KP
,
Littlejohn
EL
,
Yoo
BK
,
Salgia
R
,
Lang
D
. 
PAX6 is expressed in pancreatic cancer and actively participates in cancer progression through activation of the MET tyrosine kinase receptor gene
.
J Biol Chem
2006
;
284
:
27524
32
.
35.
Yang
Q
,
Du
WW
,
Wu
N
,
Yang
W
,
Awan
FM
,
Fang
L
, et al
A circular RNA promotes tumorigenesis by inducing c-myc nuclear translocation
.
Cell Death Differ
2017
;
24
:
1609
20
.
36.
Yang
F
,
Fang
E
,
Mei
H
,
Chen
Y
,
Li
H
,
Li
D
, et al
Cis-Acting circ-CTNNB1 promotes β-catenin signaling and cancer progression via DDX3-mediated transactivation of YY1
.
Cancer Res
2019
;
79
:
557
71
.
37.
Heumann
A
,
Kaya
Ö
,
Burdelski
C
,
Hube-Magg
C
,
Kluth
M
,
Lang
DS
, et al
Up regulation and nuclear translocation of Y-box binding protein 1 (YB-1) is linked to poor prognosis in ERG-negative prostate cancer
.
Sci Rep
2017
;
7
:
2056
.
38.
Lim
JP
,
Nair
S
,
Shyamasundar
S
,
Chua
PJ
,
Muniasamy
U
,
Matsumoto
K
, et al
Silencing Y-box binding protein-1 inhibits triple-negative breast cancer cell invasiveness via regulation of MMP1 and beta-catenin expression
.
Cancer Lett
2019l
;
452
:
119
31
.
39.
Chua
PJ
,
Lim
JP
,
Guo
TT
,
Khanna
P
,
Hu
Q
,
Bay
BH
, et al
Y-box binding protein-1 and STAT3 independently regulate ATP-binding cassette transporters in the chemoresistance of gastric cancer cells
.
Int J Oncol
2018
;
53
:
2579
89
.
40.
Lu
J
,
Li
X
,
Wang
F
,
Guo
Y
,
Huang
Y
,
Zhu
H
, et al
YB-1 expression promotes pancreatic cancer metastasis that is inhibited by microRNA-216a
.
Exp Cell Res
2017
;
359
:
319
26
.
41.
Gallardo
M
,
Hornbaker
MJ
,
Zhang
X
,
Hu
P
,
Bueso-Ramos
C
,
Post
SM
. 
Aberrant hnRNP K expression: all roads lead to cancer
.
Cell Cycle
2016
;
15
:
1552
7
.
42.
Liu
XH
,
Ma
J
,
Feng
JX
,
Feng
Y
,
Zhang
YF
,
Liu
LX
. 
Regulation and related mechanism of GSN mRNA level by hnRNPK in lung adenocarcinoma cells
.
Biol Chem
2019
;
400
:
951
63
.
43.
Gao
R
,
Yu
Y
,
Inoue
A
,
Widodo
N
,
Kaul
SC
,
Wadhwa
R
. 
Heterogeneous nuclear ribonucleoprotein K (hnRNP-K) promotes tumor metastasis by induction of genes involved in extracellular matrix, cell movement, and angiogenesis
.
J Biol Chem
2013
;
288
:
15046
56
.
44.
Lu
J
,
Gao
FH
. 
Role and molecular mechanism of heterogeneous nuclear ribonucleoprotein K in tumor development and progression
.
Biomed Rep
2016
;
4
:
657
63
.
45.
Hu
P
,
Chen
X
,
Sun
J
,
Bie
P
,
Zhang
LD
. 
siRNA-mediated knockdown against NUF2 suppresses pancreatic cancer proliferation in vitro and in vivo
.
Biosci Rep
2015
;
35
.
DOI: 10.1042/BSR20140124
.