- Open Access
circNFIB1 inhibits lymphangiogenesis and lymphatic metastasis via the miR-486-5p/PIK3R1/VEGF-C axis in pancreatic cancer
Molecular Cancer volume 19, Article number: 82 (2020)
Patients with lymph node (LN)-positive pancreatic ductal adenocarcinoma (PDAC) have extremely poor survival rates. Circular RNAs (circRNAs), a newly discovered type of endogenous noncoding RNAs, have been proposed to mediate the progression of diverse types of tumors. However, the role and underlying regulatory mechanisms of circRNAs in the LN metastasis of PDAC remain unknown.
Next-generation sequencing was used to identify differentially expressed circRNAs between PDAC and normal adjacent tissues. In vitro and in vivo experiments were conducted to evaluate the functional role of circNFIB1. RNA pulldown and luciferase assays were performed to examine the binding of circNFIB1 and miR-486-5p.
In the present study, we identified that a novel circRNA (circNFIB1, hsa_circ_0086375) was downregulated in PDAC and negatively associated with LN metastasis in PDAC patients. Functionally, circNFIB1 knockdown promoted lymphangiogenesis and LN metastasis of PDAC both in vitro and in vivo. Mechanistically, circNFIB1 functioned as a sponge of miR-486-5p, and partially reversed the effect of miR-486-5p. Moreover, circNFIB1 attenuated the oncogenic effect of miR-486-5p and consequently upregulated PIK3R1 expression, which further downregulated VEGF-C expression through inhibition of the PI3K/Akt pathway, and ultimately suppressed lymphangiogenesis and LN metastasis in PDAC.
Our findings provide novel insight into the underlying mechanism of circRNA-mediated LN metastasis of PDAC and suggest that circNFIB1 may serve as a potential therapeutic target for LN metastasis in PDAC.
Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal malignant cancers and currently ranks fourth in cancer-related death worldwide [1,2,3]. Despite the improved clinical technology, the prognosis of PDAC remains poor and its average five-year survival rate is less than 5% [3, 4]. Previous studies have demonstrated that the poor prognosis of PDAC is predominantly caused by its early invasion of the lymph nodes, which contributes to devasting dissemination into the peripheral tissues and metastasis to distant places [5, 6]. Despite the obvious clinical importance of lymphatic metastasis in PDAC, the mechanisms that lead to PDAC spread via the lymphatic vessels remain largely unknown.
Lymphatic vessels serve as the dissemination routes for cancer cells from the local tumor to the lymph nodes, and subsequently to distant locations via the blood . The expansion of lymphatic vessels, termed lymphangiogenesis, is considered to be the most critical and the rate-limiting step among the processes by which cancer cells metastasize via the LNs [8,9,10]. Previous studies have demonstrated that lymphangiogenesis is positively regulated by vascular endothelial growth factor-C (VEGF-C) and blocking the secretion of VEGF-C via its inhibitor suppresses the formation of new lymphatic vessels, which subsequently limits the LN metastasis of tumors . Nevertheless, although VEGF-C plays a crucial role in lymphangiogenesis and lymphatic metastasis, little is known of the regulatory mechanism of VEGF-C expression and its induction of lymphangiogenesis in PDAC.
Circular RNAs (circRNAs) represent a type of novel noncoding RNAs with a covalently circular structure arising from the non-canonical splicing of pre-mRNAs. Moreover, circRNAs lack 5′ caps and 3 poly-A tails and are much more stable compared with their linear counterparts since they are resistant to exonuclease . CircRNAs are widely expressed in several cancers and play essential roles in regulating cancer progression . Yang et al. found that circ-ITCH functioned as a tumor suppressor in bladder cancer through the circ-ITCH/miR-17, miR-224/p21, PTEN axis . In addition, a study conducted by Wei et al. demonstrated that circCDYL promoted the malignant transformation of hepatoma by controlling the tumor-initiating properties of hepatocellular carcinoma cells via activating the PI3K-AKT-mTORC1/β-catenin and NOTCH2 pathways . Nevertheless, the specific regulatory mechanisms of circRNAs in the LN metastasis of PDAC have not been fully explored.
MicroRNAs (miRNAs) are a kind of small endogenous non-coding RNAs comprised of 18 ~ 22 nucleotides, which mainly affect the expression of their target genes post-transcriptionally to regulate the progression of human cancers . As a member of miRNA family, miR-486-5p has been reported to function as an oncogene in multiple cancers. Lopez-Bertoni et al. found that the inhibition of miR-486-5p reduced tumor size via a PTEN-dependent mechanism in vivo . Hanna et al. reported that miR-486-5p promoted the proliferation and invasion of alveolar rhabdomyosarcoma . Although miR-486-5p plays a crucial role in human cancers, the precise regulatory mechanism of miR-486-5p in LN metastasis of PDAC is still not fully investigated.
In this study, we identified a novel circRNA (circNFIB1, has_circ_0086375), which was downregulated in PDAC and negatively associated with LN metastasis in PDAC patients. Functionally, circNFIB1 knockdown was found to promote the lymphangiogenesis and LN metastasis of PDAC. Mechanistically, circNFIB1 functioned as an miR-486-5p sponge, and partially reversed the effect of miR-486-5p. Moreover, circNFIB1 attenuated the oncogenic effect of miR-486-5p, consequently upregulating PIK3R1 expression, which further downregulated VEGF-C expression through inhibition of the PI3K/Akt pathway and ultimately resulted in the suppression of lymphangiogenesis and LN metastasis in PDAC.
With the approval of the Protection of Human Subjects Committees at Sun Yat-sen Memorial Hospital, Sun Yat-sen University (Guangzhou, China), paired PDAC tissues and non-cancerous tissues (NATs) were collected from 160 PDAC patients who underwent surgery at Sun Yat-Sen Memorial Hospital, Sun Yat-Sen University from February 2014 to February 2019. Each sample was confirmed by two independent professional pathologists. For RNA extraction, the tissues were immediately snap-frozen in liquid nitrogen and stored at − 80 °C until further investigation. For immunohistochemistry (IHC), the tissues were fixed in 10% (v/v) neutral-buffered formalin followed by 70% ethanol dehydration and then embedded in paraffin. Tumor staging was determined following the instruction of the latest edition of the tumor node metastasis (TNM) system of the American Joint Committee. Written informed consent was obtained from all the patients before performing the study. The detailed clinicopathological characteristics of the patients are summarized in Additional file 1.
Cell lines and cell culture
The human PDAC cell lines (PANC1, Capan-2, SW1990) and human pancreatic ductal endothelial cells (HPDE) were purchased from the American Type Culture Collection (ATCC, Manassas, VA). The PANC-1, Capan-2 and SW1990 cell lines were maintained in high glucose Dulbecco’s Modified Eagle’s Medium (DMEM, Gibco, USA), supplemented with 10% FBS (BI, Israel). The HPDE cell line was cultured in RPMI 1640 (Gibco, USA), supplemented with 10% FBS. Human lymphatic endothelial cells (HLECs) were obtained from the ATCC and maintained in the Endothelial Cell Medium (ECM, Gibco, USA). All cell lines were cultured under a humidified atmosphere containing 5% CO2.
Total RNA was isolated using TRIzol reagent. DNase I was applied to eliminate the genomic DNA and Ribo-Zero rRNA Removal Kits (MRZMB126, Epicentre) was applied to deplete ribosomal RNA. We then purified Poly(A) + RNA with oligo  magnetic beads and fragmented it into short sequences. The quality of the purified RNA was measured with a BioAnalyzer 2100 (AgilenN Technology, Santa Clara, CA). The library was prepared in accordance with the instructions of the TruSeq Stranded mRNA LT 11 Sample Prep Kit (Cat. No.15032612, Illumina). Each library was sequenced on an Illumina HiSeq2500 in 125PE mode 13 (Illumina, San Diego, CA, USA).
RNA pull-down assay
For the preparation of streptavidin magnetic beads, the beads were washed three times with an RNase-free lysis buffer, and blocked with 10 μg/μL BSA and yeast tRNA on a rotator at 4 °C for 3 h. The blocked beads were incubated with a biotinylated circNFIB1 probe or oligo probe at 25 °C for 2 h. Later, approximately 2 × 107 PANC1 and Capan-2 cells were harvested and fixed with 1% formaldehyde before being lysed in lysis buffer and incubated with prepared pull-down probes, 50 μL of the cell lysate was collected as the input. The pull-down complexes were harvested and tested by quantitative real-time PCR (qRT-PCR). The probes used in RNA pull-down assay were list in Additional file 2.
Popliteal LN metastasis study
All animal experiments were conducted following the guidelines approved by the Institutional Animal Care and Use Committee of Sun Yat-sen Memorial Hospital, Sun Yat-sen University. four to five weeks of age, 18 g - 20 g female BALB/c nude mice were purchased from the Experimental Animal Center. Six mice were included in each group, and 100 μL PBS suspension of 5 × 106 PANC-1 cells transduced with sh-circNFIB1-luc or sh-NC-luc were incubated into the footpads of the mice. The mice were imaged with an IVIS Spectrum Imaging System (PerkinElmer) to monitor lymphatic metastasis. When the tumors reached 200 mm3, the primary tumor and popliteal LNs were enucleated, embedded in paraffin, and analyzed with an anti-LYVE-1 antibody by IHC. Images were captured with a Nikon Eclipse 80i system with NIS-Elements software (Nikon, Japan).
All quantitative data are presented as the mean ± standard deviations from at least three independent experiments. The cumulative survival time was calculated using the Kaplan-Meier method and analyzed with a log-rank test. A multivariate Cox proportional hazards model was used to estimate the adjusted hazard ratios and 95% confidence intervals, as well as identify independent prognostic factors. Chi-square tests (χ2 tests) were used to assess the relationships between non-parametric variables, and two-tailed Student’s t-test or one-way ANOVA was used to evaluate the relationship between parametric variables. The threshold for statistical significance was set at p < 0.05. For multiple comparisons, the p-value was corrected by using Bonferroni correction. All statistical tests were performed with IBM. (Version 19.0.; IBM Corp., Armonk, NY).
Further applied methods
Additional lentivirus infection and cell transfection, qRT-PCR, HLECs tube formation assay and Transwell assay, subcellular fractionation assay, fluorescence in situ hybridization (FISH), Colocalization of circNFIB1 and miR-486a-5p, luciferase activity assay, Western Blotting, IHC, ELISA-based quantification of secreted VEGF-C, gel electrophoresis analysis are further described in the Additional file 3.
Identification and the characteristics of circNFIB1 in PDAC
To identify the essential circRNAs that contribute to the progression of PDAC, next-generation sequencing was performed to analyze five pairs of PDAC tissues and the corresponding normal adjacent tissues (NATs)(GSE136569). 13 circRNAs were differentially downregulated (fold-change < 0.4 and P < 0.05) in PDAC tissues and further validated in a larger cohort of 160-case of PDAC tissues and paired NATs by qRT-PCR analysis (Fig. 1a, b). We then focused on circNFIB1 (hsa_circ_0086375), which was expressed at low level in PDAC tissues, for further study (Fig. 1c).
Since circRNAs differ from their counterpart genes regarding their secondary structure, we next evaluated the circular structure of circNFIB1, which was aroused from exons 16 to 18 of the NFIB gene (chr9: 14146687–14,155,892) (Fig. 1d). Sanger sequencing validated the back-spliced junction of circNFIB1 (Fig. 1e), which was consistent with the sequences obtained from circBase, and the distinct product of circNFIB1 could only be amplified from the cDNA but not from the gDNA (Fig. 1f, g). Furthermore, we found that circNFIB1 expression was significantly decreased when oligo-dT primers were used for reverse transcription compared with the use of random primers, indicating the absence of poly-A tail in circRNA, which in construct, exists in linear RNAs (Fig. 1h). Using treatment with RNase R, an exoribonuclease which could only degrade linear mRNA, NFIB mRNA expression dramatically decreased, whereas circNFIB1 expression remained unaffected, further confirming that circNFIB1 consisted of a closed circular structure (Fig. 1i). Since RNAs that exist in a circular form are more stable than the linear form, we further analyzed the stability of circNFIB1. We observed a longer half-life of circNFIB1 compared with that of NFIB mRNA in PDAC cells treated with Actinomycin D, a transcription inhibitor (Fig. 1j, k). Taken together, these results indicate that circNFIB1 is a highly stable circRNA derives from exons 16 to 18 of the NFIB gene locus.
circNFIB1 is negatively associated with LN metastasis in PDAC patients
To study the role of circNFIB1 in PDAC, we analyzed circNFIB1 expression in a cohort of 160 PDAC patients. We found that patients with lymphatic metastasis and a highly pathological TMN stage had lower circNFIB1 expression (Fig. 2a, b and Additional file 1: Table S1). In addition, qRT-PCR analysis showed that there were lower levels of circNFIB1 in metastatic tumor cells in the LNs compared with the paired primary tumors (Fig. 2c). Consistent with the results obtained from clinical samples, a lower expression of circNFIB1 was detected in PDAC cells compared with HPDE (Fig. 2d). Moreover, IHC analysis showed a negative correlation between circNFIB1 and lymphatic vessel density (Fig. 2e, f). Together, our findings suggest that circNFIB1 is essential for the inhibition of the LN metastasis of PDAC.
circNFIB1 inhibits lymphangiogenesis in vitro
Given that lymphangiogenesis is a key factor for LN metastasis, we further explored whether circNFIB1 inhibits lymphangiogenesis in vitro. The expression of circNFIB1 was downregulated by transfection with siRNAs that specifically targeted the back-spliced region of circNFIB1 without affecting NFIB expression (Fig. 3a, b). Similarly, circNFIB1 was successfully overexpressed by transfection with circNFIB1 plasmids, and no obvious changes were observed with NFIB (Fig. 3c, d). In vitro assays showed that the conditioned media from circNFIB1-silencing PDAC dramatically promoted the tube formation and migration ability of HLECs (Fig. 3e-j). In contrast, conditioned media from circNFIB1-transduced PDAC significantly inhibited the HLEC tube formation and migration compared with the control group (Fig. 3k-p). Moreover, given that enhanced migration ability of cancer cells contributes to LN metastasis, we further examined the effect of circNFIB1 on the migration of PDAC cells. Transwell assays showed that circNFIB1 knockdown promoted the migration ability of PANC-1, Capan-2 and SW1990 cells, suggesting that circNFIB1 inhibited the migration of PDAC cells (Additional file 4: Fig. S1a, b). Together, these results indicate that circNFIB1 suppresses the lymphangiogenesis of PDAC in vitro.
circNFIB1 suppresses LN metastasis of PDAC in vivo
Next, we employed a popliteal LN metastasis model to examine the role of circNFIB1 in the LN metastasis of PDAC in vivo. Briefly, sh-circNFIB1#1/luc- and sh-NC/luc-transfected PDAC cells were inoculated into the footpads of nude mice. The effects of circNFIB1 on LN metastasis were analyzed when the footpad tumor size reached 200 mm3. Strikingly, the in vivo imaging showed that circNFIB1 knockdown significantly promoted the popliteal LN metastasis of PDAC cells (Fig. 4a, b). The volume of the popliteal LN was larger in the sh-circNFIB1#1 group compared with the control (Fig. 4c-e). Moreover, sh-circNFIB1#1 group exhibited higher LN metastatic rate than the control (Fig. 4f). Additionally, we found that the primary tumor size in sh-NC group was lower than the sh-circNFIB1#1 group at the same time point, suggesting that circNFIB1 suppressed the tumorigenesis of PDAC (Additional file 4: Fig. S1c). Collectively, our results indicate that circNFIB1 inhibits the LN metastasis of PDAC in vivo.
circNFIB1 functions as a miR-486-5p sponge in PDAC cells
Since the localization of RNAs within cells are crucial for their functional display, we further determined the subcellular localization of circNFIB1 in PDAC cells. FISH and subcellular fractionation assays revealed that circNFIB1 was predominantly located in the cytoplasm (Fig. 5a, b). Cytoplasm-localized circRNAs primarily function as competitive endogenous RNAs (ceRNAs) and contribute to post-transcriptional regulations . To determine whether circNFIB1 functions as an miRNA sponge in PDAC, nine miRNAs that potentially bound to circNFIB1 were predicted by CircInteractome  (Fig. 5c). Pull-down assays revealed that only miR-486-5p was enriched by circNFIB1 both in Capan-2 and PANC-1 cells (Fig. 5d, e). Moreover, luciferase reporter assays showed that miR-486-5p overexpression significantly decreased the luciferase activity of circNFIB1 compared with the NC group (Fig. 5f-h). After site-directed mutagenesis of the predicted complementary binding sites on circNFIB1, miR-486-5p failed to affect the luciferase activity of circNFIB1, which supported the sponge effect of circNFIB1 by binding to miR-486-5p on specific sequences (Fig. 5h). Furthermore, pull-down assays with biotin-labeled miR-486-5p and control showed that circNFIB1 was captured by miR-486-5p in PDAC cells, validating the interaction between circNFIB1 and miR-486-5p (Fig. 5i). Consistently, FISH assays confirmed the co-localization of circNFIB1 and miR-486-5p in the cytoplasm of PDAC cells (Fig. 5j).
Since an interaction between circNFIB1 and miR-486-5p was determined, we further investigated whether miR-486-5p mediated lymphangiogenesis in PDAC. miR-486-5p knockdown suppressed PDAC cell-induced tube formation and the migration ability of HELCs (Fig. 6a, b), whereas miR-486-5p overexpression enhanced the ability of PDAC to induce HLEC tube formation and migration (Fig. 6c, d), suggesting that miR-486-5p promoted lymphangiogenesis in PDAC cells. Taken together, these results indicate that circNFIB1 directly binds to miR-486-5p and acts as an miR-486-5p sponge in PDAC cells.
circNFIB1 antagonizes the miR-486-5p-mediated suppression of PIK3R1 expression
To identify the downstream targets regulated by miR-486-5p, the potential target genes of miR-486-5p were predicted by the intersection of TargetScan  and miRanda  (Additional file 4: Fig. S1d). Subsequently, the 95 predicted genes were subjected to pathway analysis using the PANTHER database  and 10 of the most enriched signaling pathways were recognized (Additional file 4: Fig. S1e). We further analyzed whether the predicted genes involved in these pathways were affected by miR-486-5p and found that PIK3R1 was negatively co-expressed with miR-486-5p in PDAC cells (Additional file 4: Fig. S1f-i). Moreover, miRNAs suppress target gene expression through binding to the 3’UTR of their mRNAs . We found that the 3’UTR of PIK3R1 harbored sequences complementary to miR-486-5p sequences (Fig. 6e). The luciferase assays demonstrated that miR-486-5p reduced the luciferase activity of the PIK3R1 3’UTR luciferase-construct rather than the 3’UTR luciferase-construct with mutant sequences in the miR-486-5p binding site, indicating that miR-486-5p degraded PIK3R1 by directly targeting the 3’UTR region of PIK3R1 (Fig. 6f).
Next, we investigated whether circNFIB1 could induce PIK3R1 expression as a sponge of miR-486-5p. We found that circNFIB1 silencing significantly decreased the level of PI3KR1 expression in PDAC (Fig. 6g), whereas circNFIB1 overexpression increased the level of PIK3R1 (Fig. 6h), supporting the regulatory role of circNFIB1 on PIK3R1 expression. Importantly, miR-486-5p silencing upregulated PIK3R1 in PDAC, and downregulating circNFIB1 dramatically reversed this effect (Fig. 6i-l). Collectively, these findings suggest that circNFIB1 antagonizes the inhibitory effect of miR-486-5p on PIK3R1 expression.
circNFIB1 inhibits LN metastasis via the miR-486-5p/PI3KR1/VEGF-C axis in PDAC
PIK3R1 (also known as PI3K p85α), a regulatory subunit of PI3K, has been reported to play an inhibitory role in the activation of the PI3K/Akt signaling pathway and suppress tumor progression [24, 25]. Thus, we explored whether circNFIB1 extensively regulates the activation of the PI3K/Akt signaling pathway in PDAC. The overexpression of circNFIB1 inhibited the phosphorylation of Akt and increased the level of PIK3R1 and phosphorylated GSK3β (Fig. 7a, b). Conversely, the downregulation of circNFIB1 expression facilitated the phosphorylation of Akt and reduced the level of PIK3R1 and phosphorylated GSK3β, suggesting that circNFIB1 upregulated PIK3R1 and inhibited activation of the PI3K/Akt signaling pathway in PDAC cells (Fig. 7c, d).
Since the role of circNFIB1 on the activation of the PI3K/Akt signaling pathway was clarified, the downstream targets of the PI3K/Akt signaling pathway associated with lymphangiogenesis in PDAC should be further determined. VEGF-C is required for lymphangiogenesis and considered to be a downstream regulator of the PI3K/Akt signaling pathway . qRT-PCR analysis revealed that VEGF-C expression was increased by silencing circNFIB1 and decreased by ectopic circNFIB1 in PDAC cells (Fig. 7e, f). Consistently, Western blot and ELISA assay verified that the secretion of VEGF-C by PDAC cells was inhibited by circNFIB1 (Fig. 7g, h and Additional file 4: Fig. S1j-m). Moreover, we found that circNFIB1 depletion promoted VEGF-C expression and secretion; however, following treatment with LY294002, a PI3K/Akt signaling pathway inhibitor, the expression and secretion of VEGF-C were dramatically depleted (Fig. 7i-l), further confirming that circNFIB1 induced VEGF-C suppression via inactivation of the PI3K/Akt signaling pathway.
We next evaluated whether circNFIB1-mediated VEGF-C suppression contributed to the inhibition of lymphangiogenesis and LN metastasis in PDAC. Blocking VEGF-C with neutralizing antibody, pV1006R-r, dramatically reversed circNFIB1 depletion-induced HLEC tube formation and migration in PDAC (Fig. 8a-f). Consistently, in vivo assays showed that circNFIB1 knockdown increased tumor burden in the LNs and treatment with pV1006R-r inhibited circNFIB1-depletion-mediated LN metastasis (Fig. 8g). Moreover, administration of pV1006R-r decreased the LN metastatic rate of sh-circNFIB1#1 tumor bearing mice, which significantly prolonged the survival times (Fig. 8h, i). Together, these results reveal that circNFIB1 suppresses the LN metastasis of PDAC via the miR-486-5p/PI3K/VEGF-C axis.
Clinical relevance of circNFIB1 as miR-486-5p sponge in PDAC patients
Given that the critical role of circNFIB1 in the inhibition of LN metastasis via sponging miR-486-5p was determined, we next assessed the clinical significance of circNFIB1 and miR-486-5p in a cohort of 160 PDAC patients. We found that miR-486-5p was overexpressed in PDAC tissues compared with NATs (Fig. 9a). Moreover, miR-486-5p expression in PDAC patients with LN metastasis was higher than those without LN metastasis (Fig. 9b). Importantly, Kaplan-Meier analysis revealed that high miR-486-5p expression or low circNFIB1 expression was associated with poor overall survival (OS) and disease-free survival (DFS) in 160 PDAC patients (Fig. 9c-f). Univariate and multivariate analyses demonstrated that circNFIB1 functioned as an independent factor for a good prognosis of PDAC patients (Additional file 5: Table S3 and Additional file 6: Table S4). Furthermore, correlation analysis demonstrated that VEGF-C levels were negatively correlated with circNFIB1 expression and positively associated with miR-486-5p expression in PDAC (Fig. 9g-h). In conclusion, our results indicate that circNFIB1 acts as a miR-486-5p sponge and mediates the inhibition of LN metastasis via VEGF-C suppression in PDAC patients.
CircRNAs represent a class of noncoding RNAs with a high degree of stability and participate in multiple biological processes and the regulation of gene expression in human cells . However, little is known about the role and underlying mechanisms of circRNAs in the LN metastasis of PDAC. Herein, circNFIB1 was identified to be downregulated in PDAC and negatively correlated with LN metastasis in PDAC. Functionally, circNFIB1 suppressed lymphangiogenesis and LN metastasis in PDAC. Mechanistically, circNFIB1 functioned as a ceRNA and modulated the miR-486-5p/PI3KR1/VEGF-C axis, which inhibited the lymphangiogenesis and LN metastasis of PDAC. To our knowledge, this is the first report to provide insight into the regulatory mechanism of circRNA-mediated LN metastasis inhibition of PDAC and highlight that circNFIB1 may serve as a potential therapeutic target for LN metastasis in PDAC.
Lymphatic metastasis confers a poor prognosis for patients with cancer, including those with breast cancer and bladder cancer [28, 29]. Among the multiple factors associated with LN metastasis, lymphangiogenesis is considered to be essential for LN metastasis . Yang et al. found that soluble VEGFR3 inhibited lymphangiogenesis, leading to further inhibition of lymphatic metastasis of bladder cancer . Ji et al. showed that the TNFR1-mediated VEGF-C/VEGFR3 signaling pathway promoted lymphangiogenesis, resulting in LN metastasis . However, few studies have focused on the regulation of lymphangiogenesis in PDAC. In this study, we demonstrated that circNFIB1 overexpression was negatively correlated with the lymphatic vessel density in a 160-case cohort of PDAC patients. Furthermore, circNFIB1 inhibited lymphangiogenesis in PDAC through reducing the expression of VEGF-C and ultimately limited LN metastasis in PDAC. Thus, our findings clarified the molecular mechanism underlying circNFIB1 overexpression and the inhibition LN metastasis, which provides novel insight into the regulatory mechanism of LN metastasis in PDAC.
VEGF-C is a VEGFR3 receptor ligand, which is expressed on lymphatic vascular endothelial cells and is critical for lymphangiogenesis . The disfunction of VEGF-C blocks the development of lymphatic vessels and inhibits lymphatic metastasis of cancer cells . Chen et al. found that LNMAT1 promoted the expression of CCL2 and recruited macrophages to increase VEGF-C and facilitate the lymphangiogenesis of bladder cancer . Similarly, Liu et al. demonstrated that TBL1XR1 bound to the VEGF-C promotor to induce VEGF-C expression and promoted lymphangiogenesis and LN metastasis in esophageal squamous cell carcinoma . Although the importance of VEGF-C in lymphangiogenesis and LN metastasis has been determined in multiple cancers, the mechanism underlying the regulation of VEGF-C in PDAC remains unclear. In the present study, we reported that circNFIB1 downregulated VEGF-C expression through directly sponging miR-486-5p and inactivating the PI3K/AKT signaling pathway, which further inhibited lymphangiogenesis and LN metastasis of PDAC. Therefore, our study revealed a molecular mechanism for circRNA-induced VEGF-C suppression and indicated that circNFIB1 may be a potential target by which to block VEGF-C-induced LN metastasis in PDAC patients.
The PI3K/Akt signaling pathway has been thoroughly demonstrated to be a key regulator in cancer progression. Overactivation of this pathway has been observed in a variety of cancers, including PDAC . Moreover, it has been reported that activation of the PI3K/Akt signaling pathway contributed to LN metastasis in several cancers. Yoo et al. revealed that Shh signaling promoted the lymphatic metastasis of gastric cancer through activation of the PI3K/Akt pathway . In addition, Saal et al. reported that a PIK3CA mutation induced the activation of the PI3K/Akt signaling pathway and facilitated LN metastasis in breast cancer . Nevertheless, little has been described about the PI3K/Akt pathway and LN metastasis of PDAC. In the present study, we found that circNFIB1 induced the upregulation of PIK3R1, and the regulatory subunit of PI3K could inhibit activation of the PI3K/Akt signaling pathway and reduce VEGFC expression, thereby limiting lymphangiogenesis and LN metastasis of PDAC. Moreover, blocking the PI3K/Akt signaling pathway with its inhibitor, LY294002, significantly decreased circNFIB1 depletion-induced VEGF-C, which serves as the key regulator for lymphangiogenesis. Thus, inhibition of PI3K/Akt signaling pathway may provide a novel intervention strategy for LN metastasis in PDAC.
In summary, we uncovered a novel mechanism by which circRNA-induced VEGF-C inhibition suppressed LN metastasis of PDAC. We clinically demonstrated that the downregulation of circNFIB1 is associated with the poor prognosis of PDAC patients. We further indicated that circNFIB1 overexpression could effectively inhibit LN metastasis via the miR-486-5p/PIK3R1/VEGF-C axis. Identification of circNFIB1, which inhibits LN metastasis, both expands our knowledge of LN metastasis regulation and develops a potential therapeutic target for LN metastasis in PDAC.
Availability of data and materials
Our circRNAs next-generation sequencing data used in this study have been deposited in the NCBI’s Gene Expression Omnibus and are accessible through GEO accession number: GSE136569 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE136569).
Pancreatic ductal adenocarcinoma
Vascular endothelial growth factor-C
Normal adjacent tissues
Quantitative real-time PCR
Human pancreatic ductal endothelial cells
Fluorescence in situ hybridization
Human lymphatic endothelial cells
Ellis C, Ramzy A, Kieffer TJ. Regenerative medicine and cell-based approaches to restore pancreatic function. Nat Rev Gastroenterol Hepatol. 2017;14:612–28.
Makohon-Moore A, Iacobuzio-Donahue CA. Pancreatic cancer biology and genetics from an evolutionary perspective. Nat Rev Cancer. 2016;16:553–65.
Vincent A, Herman J, Schulick R, Hruban RH, Goggins M. Pancreatic cancer. Lancet. 2011;378:607–20.
Philip PA, Mooney M, Jaffe D, Eckhardt G, Moore M, Meropol N, Emens L, O'Reilly E, Korc M, Ellis L, et al. Consensus report of the national cancer institute clinical trials planning meeting on pancreas cancer treatment. J Clin Oncol. 2009;27:5660–9.
Groot VP, Gemenetzis G, Blair AB, Rivero-Soto RJ, Yu J, Javed AA, Burkhart RA, Rinkes I, Molenaar IQ, Cameron JL, et al. Defining and predicting early recurrence in 957 patients with resected pancreatic ductal adenocarcinoma. Ann Surg. 2019;269:1154–62.
Paniccia A, Hosokawa P, Henderson W, Schulick RD, Edil BH, McCarter MD, Gajdos C. Characteristics of 10-year survivors of pancreatic ductal adenocarcinoma. JAMA Surg. 2015;150:701–10.
Sundar SS, Ganesan TS. Role of lymphangiogenesis in cancer. J Clin Oncol. 2007;25:4298–307.
Achen MG, McColl BK, Stacker SA. Focus on lymphangiogenesis in tumor metastasis. Cancer Cell. 2005;7:121–7.
Cao Y. Opinion: emerging mechanisms of tumour lymphangiogenesis and lymphatic metastasis. Nat Rev Cancer. 2005;5:735–43.
Stacker SA, Williams SP, Karnezis T, Shayan R, Fox SB, Achen MG. Lymphangiogenesis and lymphatic vessel remodelling in cancer. Nat Rev Cancer. 2014;14:159–72.
Jussila L, Alitalo K. Vascular growth factors and lymphangiogenesis. Physiol Rev. 2002;82:673–700.
Kristensen LS, Andersen MS, Stagsted LVW, Ebbesen KK, Hansen TB, Kjems J. The biogenesis, biology and characterization of circular RNAs. Nat Rev Genet. 2019;20:675–91.
Beermann J, Piccoli MT, Viereck J, Thum T. Non-coding RNAs in development and disease: background, mechanisms, and therapeutic approaches. Physiol Rev. 2016;96:1297–325.
Yang C, Yuan W, Yang X, Li P, Wang J, Han J, Tao J, Li P, Yang H, Lv Q, Zhang W. Circular RNA circ-ITCH inhibits bladder cancer progression by sponging miR-17/miR-224 and regulating p21, PTEN expression. Mol Cancer. 2018;17:19.
Wei Y, Chen X, Liang C, Ling Y, Yang X, Ye X, Zhang H, Yang P, Cui X, Ren Y, et al. A noncoding regulatory RNAs network driven by Circ-CDYL acts specifically in the early stages hepatocellular carcinoma. Hepatology. 2019;71:130–47.
Rupaimoole R, Slack FJ. MicroRNA therapeutics: towards a new era for the management of cancer and other diseases. Nat Rev Drug Discov. 2017;16:203–22.
Lopez-Bertoni H, Kotchetkov IS, Mihelson N, Lal B, Rui Y, Ames H, Lugo-Fagundo M, Guerrero-Cazares H, Quinones-Hinojosa A, Green JJ, Laterra J. A Sox2/miR-486-5p axis regulates survival of GBM cells by inhibiting tumor suppressor networks. Cancer Res. 2020;80:1644–55.
Hanna JA, Garcia MR, Lardennois A, Leavey PJ, Maglic D, Fagnan A, Go JC, Roach J, Wang YD, Finkelstein D, Hatley ME. PAX3-FOXO1 drives miR-486-5p and represses miR-221 contributing to pathogenesis of alveolar rhabdomyosarcoma. Oncogene. 2018;37:1991–2007.
Dudekula DB, Panda AC, Grammatikakis I, De S, Abdelmohsen K, Gorospe M. CircInteractome: A web tool for exploring circular RNAs and their interacting proteins and microRNAs. RNA Biol. 2016;13:34–42.
Agarwal V, Bell GW, Nam JW, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. Elife. 2015;4:e05005.
Betel D, Wilson M, Gabow A, Marks DS, Sander C. The microRNA.org resource: targets and expression. Nucleic Acids Res. 2008;36:D149–53.
Mi H, Huang X, Muruganujan A, Tang H, Mills C, Kang D, Thomas PD. PANTHER version 11: expanded annotation data from Gene Ontology and Reactome pathways, and data analysis tool enhancements. Nucleic Acids Res. 2017;45:D183–9.
Mori MA, Ludwig RG, Garcia-Martin R, Brandao BB, Kahn CR. Extracellular miRNAs: from biomarkers to mediators of physiology and disease. Cell Metab. 2019;30:656–73.
Taniguchi CM, Winnay J, Kondo T, Bronson RT, Guimaraes AR, Aleman JO, Luo J, Stephanopoulos G, Weissleder R, Cantley LC, Kahn CR. The phosphoinositide 3-kinase regulatory subunit p85alpha can exert tumor suppressor properties through negative regulation of growth factor signaling. Cancer Res. 2010;70:5305–15.
Luo J, Sobkiw CL, Logsdon NM, Watt JM, Signoretti S, O'Connell F, Shin E, Shim Y, Pao L, Neel BG, et al. Modulation of epithelial neoplasia and lymphoid hyperplasia in PTEN+/− mice by the p85 regulatory subunits of phosphoinositide 3-kinase. Proc Natl Acad Sci U S A. 2005;102:10238–43.
Brouillard P, Boon L, Vikkula M. Genetics of lymphatic anomalies. J Clin Invest. 2014;124:898–904.
Patop IL, Wust S, Kadener S. Past, present, and future of circRNAs. EMBO J. 2019;38:e100836.
McAllaster JD, Cohen MS. Role of the lymphatics in cancer metastasis and chemotherapy applications. Adv Drug Deliv Rev. 2011;63:867–75.
Chen C, Luo Y, He W, Zhao Y, Kong Y, Liu H, Zhong G, Li Y, Li J, Huang J, et al. Exosomal long noncoding RNA LNMAT2 promotes lymphatic metastasis in bladder cancer. J Clin Invest. 2020;130:404–21.
Zheng W, Aspelund A, Alitalo K. Lymphangiogenic factors, mechanisms, and applications. J Clin Invest. 2014;124:878–87.
Yang H, Kim C, Kim MJ, Schwendener RA, Alitalo K, Heston W, Kim I, Kim WJ, Koh GY. Soluble vascular endothelial growth factor receptor-3 suppresses lymphangiogenesis and lymphatic metastasis in bladder cancer. Mol Cancer. 2011;10:36.
Ji H, Cao R, Yang Y, Zhang Y, Iwamoto H, Lim S, Nakamura M, Andersson P, Wang J, Sun Y, et al. TNFR1 mediates TNF-alpha-induced tumour lymphangiogenesis and metastasis by modulating VEGF-C-VEGFR3 signalling. Nat Commun. 2014;5:4944.
Karaman S, Detmar M. Mechanisms of lymphatic metastasis. J Clin Invest. 2014;124:922–8.
Chen C, He W, Huang J, Wang B, Li H, Cai Q, Su F, Bi J, Liu H, Zhang B, et al. LNMAT1 promotes lymphatic metastasis of bladder cancer via CCL2 dependent macrophage recruitment. Nat Commun. 2018;9:3826.
Liu L, Lin C, Liang W, Wu S, Liu A, Wu J, Zhang X, Ren P, Li M, Song L. TBL1XR1 promotes lymphangiogenesis and lymphatic metastasis in esophageal squamous cell carcinoma. Gut. 2015;64:26–36.
LoRusso PM. Inhibition of the PI3K/AKT/mTOR pathway in solid tumors. J Clin Oncol. 2016;34:3803–15.
Yoo YA, Kang MH, Lee HJ, Kim BH, Park JK, Kim HK, Kim JS, Oh SC. Sonic hedgehog pathway promotes metastasis and lymphangiogenesis via activation of Akt, EMT, and MMP-9 pathway in gastric cancer. Cancer Res. 2011;71:7061–70.
Saal LH, Holm K, Maurer M, Memeo L, Su T, Wang X, Yu JS, Malmstrom PO, Mansukhani M, Enoksson J, et al. PIK3CA mutations correlate with hormone receptors, node metastasis, and ERBB2, and are mutually exclusive with PTEN loss in human breast carcinoma. Cancer Res. 2005;65:2554–9.
The authors thank Prof. J.X. Zhang, Department of Medical Statistics and Epidemiology, School of Public Health, Sun Yat-Sen University, Guangzhou, China, for statistical advice and research comments.
This study was supported by grants from the National Natural Science Foundation of China (Grant No. 81802530, 81672395, 81702951, 81672395, 81672807, 81702417, 81000917, 81402213, 81370059 and 81701715), the Guangdong Science and Technology Department (Grant No. 2018A030313564, S2012010008934, 2017A020215072, 2014A030313044, 2014A030311047, 2016A030313340, 2016A030313296, 2017A030313880 and 2017A030310200), Young Teacher Training Funding of Sun Yat-sen University (Grant No. 19ykpy121).
Ethics approval and consent to participate
Ethical consent was approved by the Committees for Ethical Review of Research involving Human Subjects at Sun Yat-sen University. Written informed consent was obtained from each patient prior to sample collection. The animal experiments were approved by the Use Committee for Animal Care at Sun Yat-sen University.
Consent for publication
The authors declare no potential conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1 Table S1.
Correlation between circNFIB1 expression and clinicopathologic characteristics of PDAC patients
Additional file 2 Table S2.
Primers and probes used in the experiments.
Additional file 3.
Additional file 4 Figure S1.
The identification of the downstream targets of miR-486-5p.
Additional file 5 Table S3.
Univariate and multivariate analyses of OS for circNFIB1 expression in PDAC patients.
Additional file 6 Table S4.
Univariate and multivariate analyses of DFS for circNFIB1 expression in PDAC patients.
Additional file 7 Table S5.
Antibodies used in the experiments.
Additional file 8 Figure S2.
Full uncut original pictures.
Additional file 9 Figure S3.
Full uncut original gels.
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, visit http://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.
About this article
Cite this article
Kong, Y., Li, Y., Luo, Y. et al. circNFIB1 inhibits lymphangiogenesis and lymphatic metastasis via the miR-486-5p/PIK3R1/VEGF-C axis in pancreatic cancer. Mol Cancer 19, 82 (2020). https://doi.org/10.1186/s12943-020-01205-6
- PI3K/Akt signaling pathway
- Lymphatic metastasis
- Pancreatic cancer