Research Paper Volume 13, Issue 14 pp 18978—18992

CircRNA circFOXK2 facilitates oncogenesis in breast cancer via IGF2BP3/miR-370 axis

Wei Zhang1, , Hui Liu1, , Junjie Jiang1, , Yunyun Yang2, , Wenjie Wang3, , Zhengyan Jia4, ,

  • 1 Department of Thyroid and Breast I, Cangzhou Central Hospital, Cangzhou, Hebei Province, China
  • 2 Outpatient Comprehensive Treatment, Cangzhou Central Hospital, Cangzhou, Hebei Province, China
  • 3 Department of General Surgery, Botou Hospital, Botou, Hebei Province, China
  • 4 Department of General Surgery, Qingxian People’s Hospital, Qingxian, Hebei Province, China

Received: April 3, 2021       Accepted: June 22, 2021       Published: July 30, 2021      

https://doi.org/10.18632/aging.203347
How to Cite

Copyright: © 2021 Zhang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Metastasis is the leading cause of breast cancer (BC)-related deaths. Circular RNAs (circRNAs) have emerged as essential regulators for cancer progression and metastasis. Therefore, the objective of this study was to investigate the role of circRNAs in BC metastasis and related mechanism. In this study, we established the BC cell line with high or low potential of metastasis. RNA sequencing, migration and invasion assay, Fluorescence in situ hybridization, luciferase report assay, circRNA pulldown, and transmission electron microscopy were performed to elucidate the molecular mechanism. The results showed that circRNA circFOXK2 was significantly increased in BC cells with high metastatic ability, and the upregulation of circFOXK2 was correlated with poor clinicopathological characteristics. Functional experiments demonstrated that overexpression of circFOXK2 promoted migration and invasion of BC cells. Also. circFOXK2 could act with IGF2BP3, an RNA-binding protein, and miR-370 to synergistically promote BC metastasis. Moreover, miR-370 could be transferred through exosomes to enhance the metastatic ability of recipient cells. In conclusion, circFOXK2 functions as a key regulator in BC metastasis, and the role of circFOXK2 on BC metastasis is tightly associated with the involvement of IGF2BP3 and miR-370. CircFOXK2 might serve as a potential biomarker for the diagnosis and treatment of BC.

Introduction

The incidence of cancer and cancer-associated death has been on the rise at 3.5% in China since 2000, making cancer a significant burden on public health care [1, 2]. Among diverse cancer types, breast cancer (BC) is the most common malignancy among females in China [3]. It is estimated that 70–80% of patients with early-stage and non-metastatic BC are curable, whereas BC with advanced-stage is not curable due to current diagnostic and therapeutic strategies [4]. Over the past decades, the treatment of BC has dramatically evolved in many aspects, including surgery, radiotherapy, chemotherapy, hormonal manipulation, or a combinational treatment [5]. The survival rate of BC, however, is still low and displays heterogeneous patterns in different regions thanks to the lack of early diagnosis and cost-effectiveness of treatments [6].

It has been demonstrated that the majority of BC-related deaths result from metastasis to other organs rather than the primary tumor itself [7]. To date, identification of BC metastasis is to detect the clinical manifestations of the metastatic organs, biopsies of metastatic organs, radiological assessments, medical imaging, as well as molecular markers [8]. Although these screenings have dramatically lowered the BC metastasis-associated mortality, these strategies mentioned above are still not sufficient and accurate for diagnosing BC at the earliest stage [7]. Therefore, it is an urgent need to explore the molecular mechanism underlying BC metastasis. Growing evidence demonstrated that metastasis is a complicated multi-step process, including invasion, intravasation, survival in blood, extravasation, as well as colonization at the distant metastatic organs [9]. As the first and most critical step of metastasis, the invasion is regarded as a migratory process of cancer cells from the primary tumor to distant other organs [10]. However, the initiation of invasion during BC progression is not fully understood.

Over the past decades, RNAs, particularly the non-coding RNAs (ncRNAs), have drawn increasing attention from both clinical and academic researchers due to the rapid development of high-throughput RNA sequencing techniques and bioinformatics [11]. Most RNAs in eukaryotic cells are ncRNAs, not messenger RNAs (mRNAs) [12]. Accumulating evidence revealed that ncRNAs, including long non-coding RNAs (lncRNAs), microRNAs (miRNAs), and circular RNAs (circRNAs), play an essential role in physiological and pathological processes, including cancers [13]. Among these ncRNAs, circRNAs have been reported to interact with miRNAs or RNA-binding proteins (RBPs) to regulate the development and progression of cancers [14, 15]. Also, circRNAs are considered as promising biomarkers for early diagnosis of cancers, such as serum cricRNAs [16]. In BC, several circRNAs are found to be aberrantly expressed in BC and participate in the carcinogenesis of BC [1719]. However, the function of circRNAs in BC metastasis has not been comprehensively investigated.

In this study, by establishing BC cells with high or low potential of metastasis, we aimed to investigate the functional cricRNA associated with BC metastasis and related molecular mechanism. This study provides insight into the molecular basis of BC metastasis and the role of circRNA in the metastatic process.

Results

CircFOXK2 is upregulated in highly metastatic BC cells

To investigate the mechanism by which BC cells are highly metastatic, we established two BT-549 cell populations with high and low potential of metastasis (BT-549-H and BT-549-L), as previously described [20, 21]. Then, we determined the characteristics of BT-549-H and BT-549-L cells. The results showed that BT-549-H cells had a higher level of migratory and invasive abilities compared with BT-549-L and BT-549 cells (Figure 1A and 1B). Also, injecting through the tail vein of mice, we found more lung metastasis and higher expression of Ki67 in mice injected with BT-549-H cells relative to BT-549-L (Figure 1C and 1D). As showed in Figure 1E, both cell populations showed similar cell viability. Next, RNA sequencing analysis was carried out to determine the differentially expressed circRNAs associated with metastasis of BC cells, as showed in a volcano plot (Figure 1F). By performing qRT-PCR, we measured the expressions of the top 10 upregulated and downregulated circRNAs, respectively. The results showed that circFOXK2 displayed the highest upregulation in BT-549-H, compared with BT-549-L (Figure 1G). Though analysis in database circBase [22], CircAtlas [23], and CircFunBase [24], circFOXK2, also named hsa_circ_0000816, was 343 nt in length and located in FOXK2 2–3 exons. Moreover, circFOXK2 was primarily expressed in the nuclei of BT-549 cells, as determined by RNA fluorescence in situ hybridization (FISH) assay (Figure 1H).

CircFOXK2 is upregulated in highly metastatic BC cells. (A) Migration ability of BT-549 cells with high and low potential of metastasis (BT-549-H and BT-549-L). Scale bar: 20 μm. (B) Invasion ability of BT-549 cells with high and low potential of metastasis (BT-549-H and BT-549-L). Scale bar: 20 μm. (C) Lung metastasis of mice injected with BT-549 or BT-549-H cells through the tail vein, as determined by H&E stain. Scale bar: 100 μm. (D) The expression of Ki67 in lung tissues of mice injected with BT-549 or BT-549-H cells through the tail vein, as determined by IHC assay. Scale bar: 40 μm. (E) Cell viability of BT-549-H and BT-549-L cells. (F) Volcano plot of differentially expressed circRNAs between BT-549-H and BT-549-L cells. (G) Heatmap for expressions of top 10 upregulated and downregulated circRNAs between BT-549-H and BT-549-L cells. (H) Cellular distribution of circFOXK2 in BT-549 cells, as detected by RNA FISH assay. Data were represented as mean ± SD. Each experimental group had at least three replicates. *p **p ***p

Figure 1. CircFOXK2 is upregulated in highly metastatic BC cells. (A) Migration ability of BT-549 cells with high and low potential of metastasis (BT-549-H and BT-549-L). Scale bar: 20 μm. (B) Invasion ability of BT-549 cells with high and low potential of metastasis (BT-549-H and BT-549-L). Scale bar: 20 μm. (C) Lung metastasis of mice injected with BT-549 or BT-549-H cells through the tail vein, as determined by H&E stain. Scale bar: 100 μm. (D) The expression of Ki67 in lung tissues of mice injected with BT-549 or BT-549-H cells through the tail vein, as determined by IHC assay. Scale bar: 40 μm. (E) Cell viability of BT-549-H and BT-549-L cells. (F) Volcano plot of differentially expressed circRNAs between BT-549-H and BT-549-L cells. (G) Heatmap for expressions of top 10 upregulated and downregulated circRNAs between BT-549-H and BT-549-L cells. (H) Cellular distribution of circFOXK2 in BT-549 cells, as detected by RNA FISH assay. Data were represented as mean ± SD. Each experimental group had at least three replicates. *p < 0.05, **p < 0.01, ***p < 0.001.

CircFOXK2 is involved in metastasis of BC

To explore the role of circFOXK2 in metastasis of BC, we measured the expression of circFOXK2 in BC tissues and found that the level of circFOXK2 was significantly higher in advanced-stage tissues compared with early-stage (Figure 2A). Correlation analysis suggested that the upregulation of circFOXK2 was correlated with invasive histological type, lymph node metastasis, and advanced stage. For liver metastasis, the expression of circFOXK2 was higher in metastatic liver tissues compared with paired BC tissues (Figure 2B), which was verified through FISH assay (Figure 2C). These results together indicated that circFOXK2 is associated with metastasis of BC.

CircFOXK2 is involved in the metastasis of BC. (A) Expression of circFOXK2 in BC tissues in the early or advanced stage. (B) Expression of circFOXK2 in BC tissues or paired liver metastasis. (C) Expression of circFOXK2 in liver metastasis, as detected by RNA FISH assay. Scale bar: 30 μm. Data were represented as mean ± SD. Each experimental group had at least three replicates. *p **p ***p

Figure 2. CircFOXK2 is involved in the metastasis of BC. (A) Expression of circFOXK2 in BC tissues in the early or advanced stage. (B) Expression of circFOXK2 in BC tissues or paired liver metastasis. (C) Expression of circFOXK2 in liver metastasis, as detected by RNA FISH assay. Scale bar: 30 μm. Data were represented as mean ± SD. Each experimental group had at least three replicates. *p < 0.05, **p < 0.01, ***p < 0.001.

CircFOXK2 is required for BC metastasis

To further explore the effect of circFOXK2 on BC metastasis, we applied the plasmid of circLONP2 and circFOXK2-specific anti-sense oligonucleotide (ASO) to overexpress and knockdown the expression of circFOXK2, respectively (Figure 3A and 3F). Meanwhile, the overexpress and knockdown circFOXK2 did not influence the expression of FOXK2 mRNA (Figure 3A and 3F). In migration and invasion assays, BT-549 cells with overexpression of circFOXK2 displayed significantly increased migratory and invasive abilities compared with control cells (Figure 3B and 3C). Also, the upregulation of circFOXK2 led to a higher level of lung metastasis and Ki67 in the mouse model (Figure 3D and 3E). On the other hand, the downregulation of circFOXK2 displayed the opposite roles (Figure 3G3J). Therefore, the results collectively demonstrated that circFOXK2 plays an essential role in BC metastasis.

CircFOXK2 is required for BC metastasis. (A) Expressions of circFOXK2 and FOXK2 mRNA in BT-549 cells transfected with circFOXK2-expressing plasmids and the negative control. (B) Migration ability of circFOXK2-overexpressing BT-549 cells. Scale bar: 20 μm. (C) Invasion ability of circFOXK2-overexpressing BT-549 cells. Scale bar: 20 μm. (D) Lung metastasis of mice injected with circFOXK2-overexpressing BT-549 cells through the tail vein, as determined by H&E stain. Scale bar: 50 μm. (E) The expression of Ki67 in lung tissues of mice injected with circFOXK2-overexpressing BT-549 cells through the tail vein, as determined by IHC assay. Scale bar: 40 μm. (F) Expressions of circFOXK2 and FOXK2 mRNA in BT-549 cells transfected with circLONP2-specific anti-sense oligonucleotide (ASO). (G) Migration ability of circFOXK2-Knockdown BT-549 cells. Scale bar: 20 μm. (H) Invasion ability of circFOXK2-Knockdown BT-549 cells. Scale bar: 20 μm. (I) Lung metastasis of mice injected with circFOXK2-Knockdown BT-549 cells through the tail vein, as determined by H&E stain. Scale bar: 50 μm. (J) The expression of Ki67 in lung tissues of mice injected with circFOXK2-Knockdown BT-549 cells through the tail vein, as determined by IHC assay. Scale bar: 40 μm. Data were represented as mean ± SD. Each experimental group had at least three replicates. *p **p ***p

Figure 3. CircFOXK2 is required for BC metastasis. (A) Expressions of circFOXK2 and FOXK2 mRNA in BT-549 cells transfected with circFOXK2-expressing plasmids and the negative control. (B) Migration ability of circFOXK2-overexpressing BT-549 cells. Scale bar: 20 μm. (C) Invasion ability of circFOXK2-overexpressing BT-549 cells. Scale bar: 20 μm. (D) Lung metastasis of mice injected with circFOXK2-overexpressing BT-549 cells through the tail vein, as determined by H&E stain. Scale bar: 50 μm. (E) The expression of Ki67 in lung tissues of mice injected with circFOXK2-overexpressing BT-549 cells through the tail vein, as determined by IHC assay. Scale bar: 40 μm. (F) Expressions of circFOXK2 and FOXK2 mRNA in BT-549 cells transfected with circLONP2-specific anti-sense oligonucleotide (ASO). (G) Migration ability of circFOXK2-Knockdown BT-549 cells. Scale bar: 20 μm. (H) Invasion ability of circFOXK2-Knockdown BT-549 cells. Scale bar: 20 μm. (I) Lung metastasis of mice injected with circFOXK2-Knockdown BT-549 cells through the tail vein, as determined by H&E stain. Scale bar: 50 μm. (J) The expression of Ki67 in lung tissues of mice injected with circFOXK2-Knockdown BT-549 cells through the tail vein, as determined by IHC assay. Scale bar: 40 μm. Data were represented as mean ± SD. Each experimental group had at least three replicates. *p < 0.05, **p < 0.01, ***p < 0.001.

IGF2BP3 is critical for the effect of circFOXK2 on BC metastasis

It has been demonstrated that RNA-binding proteins (RBPs) are critical for the function of circRNAs [25]. According to predictions of CircFunBase [24], several potential RBPs might be involved in the role of circFOXK2, including EIF4A3, FMRP, HuR, AGO2, and IGF2BP1-3 (Figure 4A). Of which, IGF2BP family has been reported to be essential to cancer development and progression [2628]. By using RNA pulldown and Western blotting assays, the results showed that circFOXK2 physically interacted with three IGF2BP family members, including IGF2BP1, IGF2BP2, and IGF2BP3, but not with EIF4A3, FMRP, HuR, and AGO2 (Figure 4B). Among three, IGF2BP3 acts as a critical factor in the regulation of cancers, such as tumor cell proliferation, invasion, and chemoresistance [29]. Thus, we speculated that IGF2BP3 might be an essential RBP mediating the effect of circFOXK2 in BC metastasis. In our rescue experiments, IGF2BP3 was overexpressed and inhibited successfully in BT-549 cells (Figure 4C). The expression of circFOXK2 positively correlated with the protein level of IGF2BP3 (Figure 4D). Also, the overexpression of IGF2BP3 promoted cell migration and invasion and reversed the effect of knockdown of circFOXK2 in BT-549 cells (Figure 4E and 4F). Meanwhile, the downregulation of IGF2BP3 attenuated the effect of overexpression of circFOXK2 on migration and invasion of BC cells (Figure 4G and 4H). These results showed that IGF2BP3 is essential for the effect of circFOXK2 on BC metastasis.

IGF2BP3 is critical for the effect of circFOXK2 on BC metastasis. (A) Prediction of RNA-binding proteins of circFOXK2. (B) Interaction between circFOXK2 and IGF2BP1, IGF2BP2, IGF2BP3, EIF4A3, FMRP, HuR, and AGO2, as determined by RNA pulldown and Western blotting assay. (C) Protein expression of IGF2BP3 in BT-549 cells transfected with IGF2BP3-expressing plasmids or IGF2BP3-specific small interfering RNA (siRNA). (D) Pearson correlation between the expressions of circFOXK2 and IGF2BP3. (E and G) Rescue experiments for the migration ability of BT-549 cells treated as indicated. Scale bar: 20 μm. (F and H) Rescue experiments for the invasion ability of BT-549 cells treated as indicated. Scale bar: 20 μm. Data were represented as mean ± SD. Each experimental group had at least three replicates. *p **p ***p

Figure 4. IGF2BP3 is critical for the effect of circFOXK2 on BC metastasis. (A) Prediction of RNA-binding proteins of circFOXK2. (B) Interaction between circFOXK2 and IGF2BP1, IGF2BP2, IGF2BP3, EIF4A3, FMRP, HuR, and AGO2, as determined by RNA pulldown and Western blotting assay. (C) Protein expression of IGF2BP3 in BT-549 cells transfected with IGF2BP3-expressing plasmids or IGF2BP3-specific small interfering RNA (siRNA). (D) Pearson correlation between the expressions of circFOXK2 and IGF2BP3. (E and G) Rescue experiments for the migration ability of BT-549 cells treated as indicated. Scale bar: 20 μm. (F and H) Rescue experiments for the invasion ability of BT-549 cells treated as indicated. Scale bar: 20 μm. Data were represented as mean ± SD. Each experimental group had at least three replicates. *p < 0.05, **p < 0.01, ***p < 0.001.

CircFOXK2-miR-370 interaction is essential for BC metastasis

The role of circRNA-miRNAs interaction has been well studied in various cancers [30, 31]. In this study, we used Circular RNA Interactome Database to predict the potential miRNAs interacting with circFOXK2 [32]. The results showed that circFOXK2 had a binding sequence of miR-370 (Figure 5A). Then, luciferase assay revealed that cells transfected with miR-370 mimic plus plasmids carrying wildtype binding sequence showed significantly decreased luciferase activity than those treated with plasmids carrying mutant binding sequence (Figure 5B). Also, circRIP assay demonstrated that miR-370 was enriched in circFOXK2-specific probes compared with those in the control group (Figure 5C). Moreover, the overexpression and knockdown of circFOXK2 increased and decreased the expression of miR-370, respectively (Figure 5D). Also, the expression of circFOXK2 displayed a positive correlation with miR-370 (Figure 5E). Collectively, these results further verified the interaction between circFOXK2 and miR-370. To further determine the role of miR-370 in BC metastasis, we applied miRNA mimic and inhibitor to overexpress and knockdown miR-370, respectively (Figure 5F and 5G). As showed in rescue experiments, the overexpression of miR-370 restored the inhibited migratory and invasive abilities in BT-549 cells with downregulation of circFOXK2, whereas the knockdown of miR-370 exerted opposite roles (Figure 5H5K). As such, these results indicated that circFOXK2-miR-370 interaction is essential for BC metastasis.

CircFOXK2-miR-370 interaction is essential for BC metastasis. (A) Putative binding site of miR-370 in the circFOXK2 sequence. (B) Luciferase activity in BT-549 cells treated as indicated. (C) circRIP assay. (D) Expression of miR-370 in circFOXK2-overexpressing/knockdown BT-549 cells. (E) Pearson correlation between the expressions of circFOXK2 and miR-370. (F) Expression of miR-370 in BT-549 cells treated with miR-370 mimic. (G) Expression of miR-370 in BT-549 cells treated with miR-370 inhibitor. (H and J) Rescue experiments for the migration ability of BT-549 cells treated as indicated. Scale bar: 20 μm. (I and K) Rescue experiments for the invasion ability of BT-549 cells treated as indicated. Scale bar: 20 μm. Data were represented as mean ± SD. Each experimental group had at least three replicates. *p **p ***p

Figure 5. CircFOXK2-miR-370 interaction is essential for BC metastasis. (A) Putative binding site of miR-370 in the circFOXK2 sequence. (B) Luciferase activity in BT-549 cells treated as indicated. (C) circRIP assay. (D) Expression of miR-370 in circFOXK2-overexpressing/knockdown BT-549 cells. (E) Pearson correlation between the expressions of circFOXK2 and miR-370. (F) Expression of miR-370 in BT-549 cells treated with miR-370 mimic. (G) Expression of miR-370 in BT-549 cells treated with miR-370 inhibitor. (H and J) Rescue experiments for the migration ability of BT-549 cells treated as indicated. Scale bar: 20 μm. (I and K) Rescue experiments for the invasion ability of BT-549 cells treated as indicated. Scale bar: 20 μm. Data were represented as mean ± SD. Each experimental group had at least three replicates. *p < 0.05, **p < 0.01, ***p < 0.001.

Exosomal miR-370 derived from BT-549-H promote BC metastasis

Exosomes play a critical role in intercellular communication in cancers; in particular, cancer cells with high metastatic potential can impact neighboring cells through exosomes [33, 34]. As shown in Figure 6A, BT-549 cells cocultured with BT-549-H displayed significantly increased migratory and invasive abilities, which was reversed the blockade of exosome generation by GW4869 [35] (Figure 6A and 6B). Then, we isolated exosomes from BT-549 (BT-549-Exo) and BT-549-H (BT-549-H-Exo) cells, respectively. The morphology and size of exosomes were identified through transmission electron microscopy, nanoparticle tracking analysis (Figure 6C). Exosomal markers CD63 and Tsg101 were measured by Western blotting assay (Figure 6D). Moreover, we found that circFOXK2 rarely expressed in exosomes, whereas miR-370 was highly expressed in BT-549-H-Exo, compared with BT-549-Exo (Figure 6E). Meanwhile, the level of miR-370 was higher in BT-549-H-Exo than BT-549-Exo (Figure 6F). Fluorescence staining assay demonstrated that exosomal miR-370 derived from BT-549-H-Exo was taken up by recipient BT-549 cells (Figure 6G). Furthermore, the overexpression and knockdown of circFOXK2 could increase and decrease both exosome generation and expression of exosomal miR-370, respectively (Figure 6H and 6I). Intriguingly, exosomes with overexpression of miR-370 promoted migration and invasion abilities in BT-549 cells, while the opposite effect was found in BT-549 cells treated with exosomes with knockdown of miR-370 (Figure 6J and 6K). These results together suggested that circFOXK2 induces metastasis of recipient cells by promoting the transfer of exosomal miR-370.

Exosomal miR-370 derived from BT-549-H promote BC metastasis. (A) Migration ability of BT-549 cells cocultured with BT-549-H cells or BT-549-H plus GW4869 treatment. Scale bar: 20 μm. (B) Invasion ability of BT-549 cells cocultured with BT-549-H cells or BT-549-H plus GW4869 treatment. Scale bar: 20 μm. (C) Morphology and size distribution of exosomes derived from BT-549 (BT-549-Exo) or BT-549-H (BT-549-H-Exo) cells, as determined by transmission electron microscopy, nanoparticle tracking analysis. (D) Expressions of exosomal markers CD63 and Tsg101, as detected by Western blotting. (E) Expressions of circFOXK2 and miR-370 in BT-549-Exo and BT-549-H-Exo. (F) Expression of miR-370 in BT-549 cells treated with BT-549-Exo and BT-549-H-Exo. (G) Exosomal miR-370 was taken up by recipient BT-549 cells, as determined by fluorescence staining assay. Scale bar: 20 μm. (H) Number of exosomes in BT-549 cells with overexpression or knockdown of circFOXK2. (I) Expression of exosomal miR-370 in BT-549 cells with overexpression or knockdown of circFOXK2. (J) Migration ability of BT-549 cells treated with exosomes containing miR-370 mimic or inhibitor. Scale bar: 20 μm. (K) Invasion ability of BT-549 cells treated with exosomes containing miR-370 mimic or inhibitor. Scale bar: 20 μm. Data were represented as mean ± SD. Each experimental group had at least three replicates. *p **p ***p

Figure 6. Exosomal miR-370 derived from BT-549-H promote BC metastasis. (A) Migration ability of BT-549 cells cocultured with BT-549-H cells or BT-549-H plus GW4869 treatment. Scale bar: 20 μm. (B) Invasion ability of BT-549 cells cocultured with BT-549-H cells or BT-549-H plus GW4869 treatment. Scale bar: 20 μm. (C) Morphology and size distribution of exosomes derived from BT-549 (BT-549-Exo) or BT-549-H (BT-549-H-Exo) cells, as determined by transmission electron microscopy, nanoparticle tracking analysis. (D) Expressions of exosomal markers CD63 and Tsg101, as detected by Western blotting. (E) Expressions of circFOXK2 and miR-370 in BT-549-Exo and BT-549-H-Exo. (F) Expression of miR-370 in BT-549 cells treated with BT-549-Exo and BT-549-H-Exo. (G) Exosomal miR-370 was taken up by recipient BT-549 cells, as determined by fluorescence staining assay. Scale bar: 20 μm. (H) Number of exosomes in BT-549 cells with overexpression or knockdown of circFOXK2. (I) Expression of exosomal miR-370 in BT-549 cells with overexpression or knockdown of circFOXK2. (J) Migration ability of BT-549 cells treated with exosomes containing miR-370 mimic or inhibitor. Scale bar: 20 μm. (K) Invasion ability of BT-549 cells treated with exosomes containing miR-370 mimic or inhibitor. Scale bar: 20 μm. Data were represented as mean ± SD. Each experimental group had at least three replicates. *p < 0.05, **p < 0.01, ***p < 0.001.

Discussion

BC initiates as a local disease and then metastasize to the lymph nodes and other distant organs [7]. Previous studies demonstrated that metastasis is the leading cause of BC-associated death [36]. Based on the widely recognized model of metastasis, a population of cancer cells in the primary tumor acquire genetic modifications over time, allowing these cells are capable of metastasizing and forming a new tumor in other distant organs [37]. Given this model, exploring the detailed mechanism of initiating cancer cell metastasizing would significantly extend our understanding of metastasis. In this study, we established a BC BT-549 cell line with a high potential of metastasis (BT-549-H) and then screened the potential circRNAs involved in the acquired metastatic ability of BT-549 cells. Also, we found that circFOXK2 was increased in BT-549-H cells and correlated with clinicopathological features of patients with BC. Moreover, the effect of circFOXK2 on BC metastasis was mediated through RNA-binding protein IGF2BP3 and miR-370. Meanwhile, miR-370 could be transferred from BT-549-H cells to recipient cells through exosomes, promoting the invasive ability of recipient cells. These results together elucidate the role of circFOXK2 in BC metastasis and related functional pathways.

In this study, by establishing a BT-549 cell population with a high potential of metastasis and RNA sequencing, we demonstrated that circFOXK2, mainly expressed in the nuclei, was associated with enhanced metastatic ability of BT-549 cell. Also, functional experiments showed that the overexpression of circFOXK2 significantly improved the migratory and invasive abilities of BT-549 cells, whereas the knockdown of circFOXK2 exerted the opposite role. According to analysis in the databases [2224], circFOXK2, also named hsa_circ_0000816, was located in 2–3 exons of the FOXK2 gene. Qiao et al. reported that the upregulation of circFOXK2 is associated with periventricular white matter damage (PWMD) of premature infants [38]. Also, circFOXK2 may participate in the regulation of myotonic dystrophy [39]. Moreover, circFOXK2, a sponge of miR-206, is predicted to be involved in asthenospermia, a common cause of human male infertility [40]. In cancers, the role of circRNAs in metastasis has been widely recognized. For example, circ_0067934 is upregulated in both tissues and cells of hepatocellular carcinoma and promotes metastasis through miR-1324/FZD5/Wnt/β-catenin pathway [41]. Also, the downregulation of circ_100395 was associated with enhanced metastasis and poor prognosis in lung cancer [42]. Wong et al. demonstrated that circFOXK2 promotes progression and metastasis of pancreatic cancer through binding with miR-942 and RBPs complex [43]. In this study, we first that circFOXK2 is a crucial regulator for the regulation of metastasis, which provides an avenue to explore the mechanism mediating the effect of circFOXK2 on BC metastasis.

Given the multiple mechanisms underlying the function of circRNAs, the interaction between circRNAs and RBPs play a critical role in transcriptional modulation, translation, and extracellular transportation [44]. In general, RBPs are a group of proteins associated with the regulation of gene expression at either the transcriptional or translational level [25]. In this study, we demonstrated that circFOXK2 could interact with IGF2BP family members, including IGF2BP1, IGF2BP2, and IGF2BP3, in BC cells. of which IGF2BP3 and circFOXK2 synergistically regulate the metastatic ability of BC cells. Accumulating studies have been shown that IGF2BP3 functions as an essential regulator in BC progression. For example, the aberrant level of IGF2BP3 is detected in the majority of invasive triple-negative breast carcinomas [45, 46], while the expression of IGF2BP3 is only elevated in adenoid cystic carcinomas in basal-like BC [47, 48]. Furthermore, a tight correlation between IGF2BP3 expression and lymph node metastasis is found in colorectal adenocarcinoma [49], oral squamous cell carcinoma [50], and gastric cancer [51]. Collectively, IGF2BP3 functions as a fine-tuner regulating the expression of genes related to cancer progression and metastasis.

In the present study, we applied online tools to predict potential miRNA interacting with circFOXK2 and identified that circFOXK2 could directly target miR-370. As an oncogenic factor, miR-370 functions as a promoter for cancer progression through targeting TGFβ-RII [52] or FOXO1 [53, 54]. Unlike the well-studied relation that circRNAs act as miRNA sponges, we observed that the expression of circFOXK2 and miR-370 displayed a similar pattern. Functionally, our rescue experiments showed that the miR-370 and circFOXK2 synergistically regulated BC cell invasion and migration. These results indicate that circFOXK2 might interact with miR-370, thereby promoting the biogenesis of miR-370. It has been reported that circRNA ciRS-7 participates in the biogenesis of mature miR-7 [14, 55]. On the other hand, the expression of miR-370 was found to be increased in exosomes derived from BT-549-H cells, and exosomal miR-370 could be internalized by recipient cells, eventually promote invasion and migration. Collectively, our observations suggested that circFOXK2 works with miR-370 to promote BC metastasis cooperatively. However, the detailed mechanism by which how circFOXK2 interacts with miR-370 should be investigated in future studies.

In conclusion, the results suggested that circFOXK2 was upregulated in metastatic BC cells and is correlated to poor clinicopathological features of BC patients. Functionally, circFOXK2 promotes invasion and migration of BC cells, and the effect of circFOXK2 on BC metastasis is associated with the involvement of IGF2BP3 and miR-370. This study indicates that circFOXK2 might serve as a biomarker for the diagnosis and treatment of BC.

Materials and Methods

Ethics statement

All patients were informed before inclusion, and the written consents were given. All experiments were approved by the ethics committee of our hospital. All animal experiments complied with the guidelines of the Animal Ethics Committee for the care and use of our hospital.

Patients

BC primary tissues were collected from patients with BC who underwent operation between March 2015 and March 2017 at our hospital. The exclusion criteria were as follows: 1) suffering from other malignancies; 2) patients had previous treatment; 3) histologic diagnosis was not BC, and 4) patients had not complete data of analysis. All samples were stored at –80°C until use.

Cell culture

BC cell line BT-549 and HEK293T cells were purchased from the cell bank of the Chinese Academy of Sciences (Shanghai, China). BT-549 cells were cultured in RPMI 1640 media supplemented with 10% fetal bovine serum (FBS) (Gibco, USA). HEK293T cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM) with 10% FBS (Gibco, USA). Cells were maintained in an incubator (Thermo Fisher Scientific, USA) at 37°C with 5% CO2.

Migration and invasion assay

To establish BT-549 cell population with high and low potential of metastasis (BT-549-H and BT-549-L), BT-549 cells were subjected to repetitive invasion assay, as previously described [20, 21]. After incubation for 36 hours at 37°C invaded cells underneath the membrane and uninvaded cells were collected and expanded for the next round of screen. After 10 and 30 screen rounds, cell populations were classified as BT-549-L and BT-549-H cells, respectively. Migratory and invasive abilities were determined using QCM Chemotaxis Cell Migration Assay (24-well, 8 μm) and QCM ECMatrix Cell Invasion Assay (24-well, 8 μm) (Sigma-Aldrich, USA) according to the manufacturer’s instruction. The data were quantified in 5 random places under a microscope.

CCK-8 assay

Cell viability was determined using Cell Counting Kit 8 (WST-8/CCK8) (Abcam, Japan) according to the manufacturer’s instructions. The BT-549-H and BT-549-L cells were seeded (1 × 105 cells/well) in a 96-well dish. OD values were measured by absorbance at 460 nm at 1, 2, 3, 4, and 5 day.

RNA sequencing

RNA sequencing assay was carried out between BT-549-H and BT-549-L cells, as previously described [56]. In brief, total RNAs were isolated from cells, exosomes, tissues using Trizol (Invitrogen, USA). RNA quality and concentration were determined by NanoDrop™ 2000 (Thermo Scientific, USA). RNA sequencing libraries were established and sequenced by Beyotime Biotechnology (Shanghai, China). Ribosomal RNAs (rRNAs) were removed from total RNAs (5 μg) using Ribo-Zero Plus rRNA Depletion Kit (Illumina, USA). Linear RNAs were digested using RNase R (New England Biolabs Inc, USA). Sequencing libraries were established using NEBNext Ultra RNA Library Prep Kit for Illumina (New England Biolabs Inc, USA) according to the manufacturer’s instructions. RNA samples were fragmented into pieces of ~ 300 bp in length, and the first-strand cDNAs were synthesized by reverse transcription and random hexamer primers. Afterward, the second-strand cDNAs were synthesized using Second Strand Synthesis Reaction Buffer. The final cDNA fragments were applied to the end modification processes, such as the addition of a single “A” base and the ligation of the adapters. Then, the chain specific libraries were constructed using USER Enzyme (New England Biolabs Inc, USA) and amplified by PCR. The libraries were qualified by NEBNext® Library Quant Kit for Illumina (New England Biolabs Inc, USA). Lastly, the libraries were subjected to sequencing assay on an Illumina HiSeq sequencer system (Illumina, USA). The sequencing quality of raw data was evaluated by FastQC software [56]. High-quality reads were aligned to the human reference genome (GRCh38/hg38) using Tophat2 software [57] with default parameters. Unaligned reads were used for subsequent circRNA analysis using CIRCexplorer2 [58] and Find_circ [59]. Differential expression analysis between BT-549-H and BT-549-L was carried out using Limma (v3.32.10) R package [60]. CricRNA with Q value > 0.01 was defined as significantly differential expression.

Quantitative real-time PCR

Total RNAs were isolated from cells, exosomes, tissues using Trizol (Invitrogen, USA). RNA quality and concentration were determined by NanoDrop™ 2000 (Thermo Scientific, USA). cDNA was synthesized using PrimeScript™ RT reagent Kit (TaKaRa, China). Quantitative real-time PCR was carried out using TB Green™ Premix Ex Taq™ II (TaKaRa, China) on ABI 7500 real-time PCR system (Applied Biosystems, USA). The relative expression was calculated using 2−ΔΔCt method [61]. β-actin and U6 were used as reference genes.

Western blotting

Total protein of cells or exosomes was isolated by using the cell lysis buffer (Beyotime Institute of Biotechnology, China). The western blotting assay was performed as previously reported [62]. The primary antibodies for IGF2BP3, CD63, and Tsg101 were obtained from Santa Cruz Biotechnology (Santa Cruz Biotechnology, USA). Optical density was quantified by the Uvitec Alliance software (Eppendorf, Germany).

Cell transfection

Anti-sense oligonucleotide (ASO) for circFOXK2 knockdown and pcDNA3.1vectors for circFOXK2 overexpression were obtained from RiboBio (Guangzhou, China). The sequence information were as following: ASO-circFOXK2: 5′-GAAGGUGCACAUUCAGGUUTT-3′; ASO-negative control: 5′-TTCTCCGAACGTGTCACGT-3′; pcDNA3.1 vector-circFOXK2 (forward): 5′-GCGATATCGTGCACATTCAGGTTCCCGAG-3′, pcDNA3.1 vector-circFOXK2 (reverse): 5′-GCCCGCGGCTTCGGGCTGTCTCCA-3′. ASO negative control and pcDNA3.1vectors were used as the negative control, respectively. Lentiviral miR-370 mimic, inhibitor, respective negative control were obtained from Qiagen (Hilden, Germany). Plasmids for IGF2BP3 overexpression and siRNAs for IGF2BP3 knockdown were purchased from Qiagen (Hilden, Germany). Cell transfection was performed using Lipofectamine™ 3000 Transfection Reagent (Invitrogen, USA) according to the manufacturer’s instructions.

Fluorescence in situ hybridization (FISH)

FISH assay was carried out as previously described [63]. In brief, circFOXK2-specific probes marked with Digoxigenin (DIG) -11-uridine triphosphate (UTP) (Roche, USA) was used. Cells were fixed with 4% paraformaldehyde for 10 min and permeabilized in 0.5% Triton X-100 in PBS solution for 5 min. Cells were then hybridized with the probe at 37°C for 16 hours. Then, the cells were washed with sodium citrate containing 0.1% Tween-20 for 5 min and then saline-sodium citrate (SSC) buffer for 5 min. Cells were stained with 4′,6-diamidino-2-phenylindole (DAPI) (Invitrogen, USA) for 10 min. Images were taken using SP8 laser confocal microscopy (Leica, Germany).

Luciferase report assay

The conserved sequences containing the putative binding site of miR-370 were synthesized from the circFOXK2 sequence and then were cloned into pGL3-enhancer vector (Promega Corporation, USA). The mutated binding site of miR-370 was also cloned into the same luciferase reporter. Luciferase reporter plasmids and miR-370 mimic/negative control were transfected into HEK293T cell using Lipofectamine™ 3000 Transfection Reagent (Invitrogen, USA) according to the manufacturer’s instructions. After 24 hours of transfection, relative luciferase activity was determined using Dual-Luciferase® Reporter Assay System (Promega Corporation, USA) according to the manufacturer’s instructions.

RNA and circRNA pulldown

Biotin-labeled circFOXK2-specific probe and negative control probe were (Sangon Biotech, China) were used for RNA pulldown assay. The assay was performed as previously described [64, 65]. RNAs attached to the beads were isolated using Trizol (Invitrogen, USA) and measured by qRT-PCR. Proteins attached to the beads were measured by Western blotting.

Exosome isolation and identification

Exosomes were isolated from BT-549 and BT-549-H cultured medium with exosome-free PBS using Total Exosome Isolation Reagent (Invitrogen, USA) according to the manufacturer’s instructions. Then, transmission electron microscopy (TEM) and nanoparticle tracking analysis (NTS) were performed to determine exosome morphology, and size, and concentration, as previously described [66]. The concentration of exosomal proteins was determined using Pierce BCA Protein Assay Kit (Thermo Fisher Scientific, USA).

Exosome electroporation

MiR-370 mimic or inhibitor was loaded into exosomes through electroporation assay using Gene Pulser Xcell Electroporation Systems (BioRad, USA) as previously described [67].

Cell-exosome coculture

BT-549 and BT-549-H cells (1 × 105) were placed in the inner chamber of 24 transwell plates. Recipient cells were placed in the outer chamber. Twenty-four hours later, cells in the outer chamber were used to determine the migration and invasion abilities.

Animal study

Male BALB/c athymic nude mice (6–8 weeks) were purchased from Shanghai Laboratory Animal Center (Shanghai, China) and maintained in standard conditions in the animal facility at our hospital. CircFOXK2-overexpressing/knockdown BT-549 or BT-549 cells (1 × 105) were injected into mice through the lateral tail vein (n = 8 per group). Eight weeks post-injection, mice were sacrificed, and lung and liver tissues were collected. The number of tumor nodules was quantified through hematoxylin and eosin (HE) staining and immunohistochemistry (IHC) assay, as previously described [68, 69].

Statistical analysis

Data were represented as mean ± SD. Data were analyzed by SPSS 18.0 software (SPSS Inc, USA). Pearson correlation was used to analyze the correlation between circFOXK2 and IGF2BP3, and circFOXK2 and miR-370. Mean differences between groups were analyzed using the Tukey’s test. In this study, differences were regarded to be significant at p < 0.05.

Author Contributions

All authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by Wei Zhang, Hui Liu, Junjie Jiang, Yunyun Yang. The first draft of the manuscript was written by Wenjie Wang, Zhengyan Jia and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Funding

This research was supported by our hospital: Cangzhou Central Hospital.

References

  • 1. Chen W, Zheng R, Baade PD, Zhang S, Zeng H, Bray F, Jemal A, Yu XQ, He J. Cancer statistics in China, 2015. CA Cancer J Clin. 2016; 66:115–32. https://doi.org/10.3322/caac.21338 [PubMed]
  • 2. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018; 68:394–424. https://doi.org/10.3322/caac.21492 [PubMed]
  • 3. Fan L, Strasser-Weippl K, Li JJ, St Louis J, Finkelstein DM, Yu KD, Chen WQ, Shao ZM, Goss PE. Breast cancer in China. Lancet Oncol. 2014; 15:e279–89. https://doi.org/10.1016/S1470-2045(13)70567-9 [PubMed]
  • 4. Harbeck N, Penault-Llorca F, Cortes J, Gnant M, Houssami N, Poortmans P, Ruddy K, Tsang J, Cardoso F. Breast cancer. Nat Rev Dis Primers. 2019; 5:66. https://doi.org/10.1038/s41572-019-0111-2 [PubMed]
  • 5. Kutanzi KR, Yurchenko OV, Beland FA, Checkhun VF, Pogribny IP. MicroRNA-mediated drug resistance in breast cancer. Clin Epigenetics. 2011; 2:171–85. https://doi.org/10.1007/s13148-011-0040-8 [PubMed]
  • 6. Sankaranarayanan R, Swaminathan R, Brenner H, Chen K, Chia KS, Chen JG, Law SC, Ahn YO, Xiang YB, Yeole BB, Shin HR, Shanta V, Woo ZH, et al. Cancer survival in Africa, Asia, and Central America: a population-based study. Lancet Oncol. 2010; 11:165–73. https://doi.org/10.1016/S1470-2045(09)70335-3 [PubMed]
  • 7. Weigelt B, Peterse JL, van't Veer LJ. Breast cancer metastasis: markers and models. Nat Rev Cancer. 2005; 5:591–602. https://doi.org/10.1038/nrc1670 [PubMed]
  • 8. Sun YF, Yang XR, Zhou J, Qiu SJ, Fan J, Xu Y. Circulating tumor cells: advances in detection methods, biological issues, and clinical relevance. J Cancer Res Clin Oncol. 2011; 137:1151–73. https://doi.org/10.1007/s00432-011-0988-y [PubMed]
  • 9. Tang DJ, Dong SS, Ma NF, Xie D, Chen L, Fu L, Lau SH, Li Y, Li Y, Guan XY. Overexpression of eukaryotic initiation factor 5A2 enhances cell motility and promotes tumor metastasis in hepatocellular carcinoma. Hepatology. 2010; 51:1255–63. https://doi.org/10.1002/hep.23451 [PubMed]
  • 10. Clark AG, Vignjevic DM. Modes of cancer cell invasion and the role of the microenvironment. Curr Opin Cell Biol. 2015; 36:13–22. https://doi.org/10.1016/j.ceb.2015.06.004 [PubMed]
  • 11. Salzman J, Gawad C, Wang PL, Lacayo N, Brown PO. Circular RNAs are the predominant transcript isoform from hundreds of human genes in diverse cell types. PLoS One. 2012; 7:e30733. https://doi.org/10.1371/journal.pone.0030733 [PubMed]
  • 12. Patop IL, Kadener S. circRNAs in Cancer. Curr Opin Genet Dev. 2018; 48:121–27. https://doi.org/10.1016/j.gde.2017.11.007 [PubMed]
  • 13. Shang Q, Yang Z, Jia R, Ge S. The novel roles of circRNAs in human cancer. Mol Cancer. 2019; 18:6. https://doi.org/10.1186/s12943-018-0934-6 [PubMed]
  • 14. Memczak S, Jens M, Elefsinioti A, Torti F, Krueger J, Rybak A, Maier L, Mackowiak SD, Gregersen LH, Munschauer M, Loewer A, Ziebold U, Landthaler M, et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013; 495:333–38. https://doi.org/10.1038/nature11928 [PubMed]
  • 15. Arnaiz E, Sole C, Manterola L, Iparraguirre L, Otaegui D, Lawrie CH. CircRNAs and cancer: Biomarkers and master regulators. Semin Cancer Biol. 2019; 58:90–99. https://doi.org/10.1016/j.semcancer.2018.12.002 [PubMed]
  • 16. Zhang HD, Jiang LH, Sun DW, Hou JC, Ji ZL. CircRNA: a novel type of biomarker for cancer. Breast Cancer. 2018; 25:1–7. https://doi.org/10.1007/s12282-017-0793-9 [PubMed]
  • 17. Yang R, Xing L, Zheng X, Sun Y, Wang X, Chen J. The circRNA circAGFG1 acts as a sponge of miR-195-5p to promote triple-negative breast cancer progression through regulating CCNE1 expression. Mol Cancer. 2019; 18:4. https://doi.org/10.1186/s12943-018-0933-7 [PubMed]
  • 18. Wang H, Xiao Y, Wu L, Ma D. Comprehensive circular RNA profiling reveals the regulatory role of the circRNA-000911/miR-449a pathway in breast carcinogenesis. Int J Oncol. 2018; 52:743–54. https://doi.org/10.3892/ijo.2018.4265 [PubMed]
  • 19. Tang YY, Zhao P, Zou TN, Duan JJ, Zhi R, Yang SY, Yang DC, Wang XL. Circular RNA hsa_circ_0001982 Promotes Breast Cancer Cell Carcinogenesis Through Decreasing miR-143. DNA Cell Biol. 2017; 36:901–08. https://doi.org/10.1089/dna.2017.3862 [PubMed]
  • 20. Feng Y, Zhu J, Ou C, Deng Z, Chen M, Huang W, Li L. MicroRNA-145 inhibits tumour growth and metastasis in colorectal cancer by targeting fascin-1. Br J Cancer. 2014; 110:2300–09. https://doi.org/10.1038/bjc.2014.122 [PubMed]
  • 21. Wang CC, Tsai MF, Dai TH, Hong TM, Chan WK, Chen JJ, Yang PC. Synergistic activation of the tumor suppressor, HLJ1, by the transcription factors YY1 and activator protein 1. Cancer Res. 2007; 67:4816–26. https://doi.org/10.1158/0008-5472.CAN-07-0504 [PubMed]
  • 22. Maass PG, Glažar P, Memczak S, Dittmar G, Hollfinger I, Schreyer L, Sauer AV, Toka O, Aiuti A, Luft FC, Rajewsky N. A map of human circular RNAs in clinically relevant tissues. J Mol Med (Berl). 2017; 95:1179–89. https://doi.org/10.1007/s00109-017-1582-9 [PubMed]
  • 23. Wu W, Ji P, Zhao F. CircAtlas: an integrated resource of one million highly accurate circular RNAs from 1070 vertebrate transcriptomes. Genome Biol. 2020; 21:101. https://doi.org/10.1186/s13059-020-02018-y [PubMed]
  • 24. Meng X, Hu D, Zhang P, Chen Q, Chen M. CircFunBase: a database for functional circular RNAs. Database (Oxford). 2019; 2019:baz003. https://doi.org/10.1093/database/baz003 [PubMed]
  • 25. Zang J, Lu D, Xu A. The interaction of circRNAs and RNA binding proteins: An important part of circRNA maintenance and function. J Neurosci Res. 2020; 98:87–97. https://doi.org/10.1002/jnr.24356 [PubMed]
  • 26. Bell JL, Wächter K, Mühleck B, Pazaitis N, Köhn M, Lederer M, Hüttelmaier S. Insulin-like growth factor 2 mRNA-binding proteins (IGF2BPs): post-transcriptional drivers of cancer progression? Cell Mol Life Sci. 2013; 70:2657–75. https://doi.org/10.1007/s00018-012-1186-z [PubMed]
  • 27. Huang X, Huang M, Kong L, Li Y. miR-372 suppresses tumour proliferation and invasion by targeting IGF2BP1 in renal cell carcinoma. Cell Prolif. 2015; 48:593–99. https://doi.org/10.1111/cpr.12207 [PubMed]
  • 28. Fawzy IO, Hamza MT, Hosny KA, Esmat G, Abdelaziz AI. Abrogating the interplay between IGF2BP1, 2 and 3 and IGF1R by let-7i arrests hepatocellular carcinoma growth. Growth Factors. 2016; 34:42–50. https://doi.org/10.3109/08977194.2016.1169532 [PubMed]
  • 29. Lederer M, Bley N, Schleifer C, Hüttelmaier S. The role of the oncofetal IGF2 mRNA-binding protein 3 (IGF2BP3) in cancer. Semin Cancer Biol. 2014; 29:3–12. https://doi.org/10.1016/j.semcancer.2014.07.006 [PubMed]
  • 30. Cheng J, Zhuo H, Xu M, Wang L, Xu H, Peng J, Hou J, Lin L, Cai J. Regulatory network of circRNA-miRNA-mRNA contributes to the histological classification and disease progression in gastric cancer. J Transl Med. 2018; 16:216. https://doi.org/10.1186/s12967-018-1582-8 [PubMed]
  • 31. Verduci L, Strano S, Yarden Y, Blandino G. The circRNA-microRNA code: emerging implications for cancer diagnosis and treatment. Mol Oncol. 2019; 13:669–80. https://doi.org/10.1002/1878-0261.12468 [PubMed]
  • 32. 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. https://doi.org/10.1080/15476286.2015.1128065 [PubMed]
  • 33. Melo SA, Sugimoto H, O'Connell JT, Kato N, Villanueva A, Vidal A, Qiu L, Vitkin E, Perelman LT, Melo CA, Lucci A, Ivan C, Calin GA, Kalluri R. Cancer exosomes perform cell-independent microRNA biogenesis and promote tumorigenesis. Cancer Cell. 2014; 26:707–21. https://doi.org/10.1016/j.ccell.2014.09.005 [PubMed]
  • 34. Steinbichler TB, Dudás J, Riechelmann H, Skvortsova II. The role of exosomes in cancer metastasis. Semin Cancer Biol. 2017; 44:170–81. https://doi.org/10.1016/j.semcancer.2017.02.006 [PubMed]
  • 35. Essandoh K, Yang L, Wang X, Huang W, Qin D, Hao J, Wang Y, Zingarelli B, Peng T, Fan GC. Blockade of exosome generation with GW4869 dampens the sepsis-induced inflammation and cardiac dysfunction. Biochim Biophys Acta. 2015; 1852:2362–71. https://doi.org/10.1016/j.bbadis.2015.08.010 [PubMed]
  • 36. Giavazzi R, Alessandri G, Spreafico F, Garattini S, Mantovani A. Metastasizing capacity of tumour cells from spontaneous metastases of transplanted murine tumours. Br J Cancer. 1980; 42:462–72. https://doi.org/10.1038/bjc.1980.259 [PubMed]
  • 37. Fidler IJ, Kripke ML. Metastasis results from preexisting variant cells within a malignant tumor. Science. 1977; 197:893–95. https://doi.org/10.1126/science.887927 [PubMed]
  • 38. Qiao L, Mo S, Zhou Y, Zhang Y, Li B, Wu S, Lin L, Zhu L, Zhao R. Circular RNA expression alteration in whole blood of premature infants with periventricular white matter damage. Genomics. 2020; 112:2875–85. https://doi.org/10.1016/j.ygeno.2020.03.027 [PubMed]
  • 39. Czubak K, Taylor K, Piasecka A, Sobczak K, Kozlowska K, Philips A, Sedehizadeh S, Brook JD, Wojciechowska M, Kozlowski P. Global Increase in Circular RNA Levels in Myotonic Dystrophy. Front Genet. 2019; 10:649. https://doi.org/10.3389/fgene.2019.00649 [PubMed]
  • 40. Liu L, Li T, Li F, Zhao X, Zhang R, Liu J, Zhang W, Lu J, Zhang X, Ma X. The influence of l-carnitine on the expression of miRNAs in asthenospermia spermatozoa and the network regulation of the associated molecules. Andrologia. 2020; 52:e13478. https://doi.org/10.1111/and.13478 [PubMed]
  • 41. Zhu Q, Lu G, Luo Z, Gui F, Wu J, Zhang D, Ni Y. CircRNA circ_0067934 promotes tumor growth and metastasis in hepatocellular carcinoma through regulation of miR-1324/FZD5/Wnt/β-catenin axis. Biochem Biophys Res Commun. 2018; 497:626–32. https://doi.org/10.1016/j.bbrc.2018.02.119 [PubMed]
  • 42. Chen D, Ma W, Ke Z, Xie F. CircRNA hsa_circ_100395 regulates miR-1228/TCF21 pathway to inhibit lung cancer progression. Cell Cycle. 2018; 17:2080–90. https://doi.org/10.1080/15384101.2018.1515553 [PubMed]
  • 43. Wong CH, Lou UK, Li Y, Chan SL, Tong JH, To KF, Chen Y. CircFOXK2 Promotes Growth and Metastasis of Pancreatic Ductal Adenocarcinoma by Complexing with RNA-Binding Proteins and Sponging MiR-942. Cancer Res. 2020; 80:2138–49. https://doi.org/10.1158/0008-5472.CAN-19-3268 [PubMed]
  • 44. Abdelmohsen K, Panda AC, Munk R, Grammatikakis I, Dudekula DB, De S, Kim J, Noh JH, Kim KM, Martindale JL, Gorospe M. Identification of HuR target circular RNAs uncovers suppression of PABPN1 translation by CircPABPN1. RNA Biol. 2017; 14:361–69. https://doi.org/10.1080/15476286.2017.1279788 [PubMed]
  • 45. Walter O, Prasad M, Lu S, Quinlan RM, Edmiston KL, Khan A. IMP3 is a novel biomarker for triple negative invasive mammary carcinoma associated with a more aggressive phenotype. Hum Pathol. 2009; 40:1528–33. https://doi.org/10.1016/j.humpath.2009.05.005 [PubMed]
  • 46. Sidoni A, Cartaginese F. IMP3 expression in triple-negative breast carcinoma. Hum Pathol. 2010; 41:1355–56. https://doi.org/10.1016/j.humpath.2010.05.004 [PubMed]
  • 47. Won JR, Gao D, Chow C, Cheng J, Lau SY, Ellis MJ, Perou CM, Bernard PS, Nielsen TO. A survey of immunohistochemical biomarkers for basal-like breast cancer against a gene expression profile gold standard. Mod Pathol. 2013; 26:1438–50. https://doi.org/10.1038/modpathol.2013.97 [PubMed]
  • 48. Vranic S, Gurjeva O, Frkovic-Grazio S, Palazzo J, Tawfik O, Gatalica Z. IMP3, a proposed novel basal phenotype marker, is commonly overexpressed in adenoid cystic carcinomas but not in apocrine carcinomas of the breast. Appl Immunohistochem Mol Morphol. 2011; 19:413–16. https://doi.org/10.1097/PAI.0b013e3182143399 [PubMed]
  • 49. Yuan RH, Wang CC, Chou CC, Chang KJ, Lee PH, Jeng YM. Diffuse expression of RNA-binding protein IMP3 predicts high-stage lymph node metastasis and poor prognosis in colorectal adenocarcinoma. Ann Surg Oncol. 2009; 16:1711–19. https://doi.org/10.1245/s10434-009-0446-0 [PubMed]
  • 50. Li S, Cha J, Kim J, Kim KY, Kim HJ, Nam W, Cha IH. Insulin-like growth factor II mRNA-binding protein 3: a novel prognostic biomarker for oral squamous cell carcinoma. Head Neck. 2011; 33:368–74. https://doi.org/10.1002/hed.21457 [PubMed]
  • 51. Kim HJ, Kim GE, Lee JS, Lee JH, Nam JH, Choi C. Insulin-like growth factor-II mRNA-binding protein 3 expression in effusion cytology: a marker for metastatic adenocarcinoma cells and a potential prognostic indicator in gastric adenocarcinoma. Acta Cytol. 2014; 58:167–73. https://doi.org/10.1159/000357199 [PubMed]
  • 52. Lo SS, Hung PS, Chen JH, Tu HF, Fang WL, Chen CY, Chen WT, Gong NR, Wu CW. Overexpression of miR-370 and downregulation of its novel target TGFβ-RII contribute to the progression of gastric carcinoma. Oncogene. 2012; 31:226–37. https://doi.org/10.1038/onc.2011.226 [PubMed]
  • 53. Fan C, Liu S, Zhao Y, Han Y, Yang L, Tao G, Li Q, Zhang L. Upregulation of miR-370 contributes to the progression of gastric carcinoma via suppression of FOXO1. Biomed Pharmacother. 2013; 67:521–26. https://doi.org/10.1016/j.biopha.2013.04.014 [PubMed]
  • 54. Wu Z, Sun H, Zeng W, He J, Mao X. Upregulation of MircoRNA-370 induces proliferation in human prostate cancer cells by downregulating the transcription factor FOXO1. PLoS One. 2012; 7:e45825. https://doi.org/10.1371/journal.pone.0045825 [PubMed]
  • 55. Hansen TB, Jensen TI, Clausen BH, Bramsen JB, Finsen B, Damgaard CK, Kjems J. Natural RNA circles function as efficient microRNA sponges. Nature. 2013; 495:384–88. https://doi.org/10.1038/nature11993 [PubMed]
  • 56. Li XN, Wang ZJ, Ye CX, Zhao BC, Li ZL, Yang Y. RNA sequencing reveals the expression profiles of circRNA and indicates that circDDX17 acts as a tumor suppressor in colorectal cancer. J Exp Clin Cancer Res. 2018; 37:325. https://doi.org/10.1186/s13046-018-1006-x [PubMed]
  • 57. 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. https://doi.org/10.1186/gb-2013-14-4-r36 [PubMed]
  • 58. Zhang XO, Wang HB, Zhang Y, Lu X, Chen LL, Yang L. Complementary sequence-mediated exon circularization. Cell. 2014; 159:134–47. https://doi.org/10.1016/j.cell.2014.09.001 [PubMed]
  • 59. Jeck WR, Sorrentino JA, Wang K, Slevin MK, Burd CE, Liu J, Marzluff WF, Sharpless NE. Circular RNAs are abundant, conserved, and associated with ALU repeats. RNA. 2013; 19:141–57. https://doi.org/10.1261/rna.035667.112 [PubMed]
  • 60. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015; 43:e47. https://doi.org/10.1093/nar/gkv007 [PubMed]
  • 61. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001; 25:402–08. https://doi.org/10.1006/meth.2001.1262 [PubMed]
  • 62. Pulito C, Mori F, Sacconi A, Goeman F, Ferraiuolo M, Pasanisi P, Campagnoli C, Berrino F, Fanciulli M, Ford RJ, Levrero M, Pediconi N, Ciuffreda L, et al. Metformin-induced ablation of microRNA 21-5p releases Sestrin-1 and CAB39L antitumoral activities. Cell Discov. 2017; 3:17022. https://doi.org/10.1038/celldisc.2017.22 [PubMed]
  • 63. Zhang J, Liu H, Hou L, Wang G, Zhang R, Huang Y, Chen X, Zhu J. 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. https://doi.org/10.1186/s12943-017-0719-3 [PubMed]
  • 64. Zeng K, Chen X, Xu M, Liu X, Hu X, Xu T, Sun H, Pan Y, He B, Wang S. CircHIPK3 promotes colorectal cancer growth and metastasis by sponging miR-7. Cell Death Dis. 2018; 9:417. https://doi.org/10.1038/s41419-018-0454-8 [PubMed]
  • 65. Han D, Li J, Wang H, Su X, Hou J, Gu Y, Qian C, Lin Y, Liu X, Huang M, Li N, Zhou W, Yu Y, Cao X. Circular RNA circMTO1 acts as the sponge of microRNA-9 to suppress hepatocellular carcinoma progression. Hepatology. 2017; 66:1151–64. https://doi.org/10.1002/hep.29270 [PubMed]
  • 66. Zhang X, Sai B, Wang F, Wang L, Wang Y, Zheng L, Li G, Tang J, Xiang J. Hypoxic BMSC-derived exosomal miRNAs promote metastasis of lung cancer cells via STAT3-induced EMT. Mol Cancer. 2019; 18:40. https://doi.org/10.1186/s12943-019-0959-5 [PubMed]
  • 67. Yang J, Zhang X, Chen X, Wang L, Yang G. Exosome Mediated Delivery of miR-124 Promotes Neurogenesis after Ischemia. Mol Ther Nucleic Acids. 2017; 7:278–87. https://doi.org/10.1016/j.omtn.2017.04.010 [PubMed]
  • 68. Wang XP, Zhou J, Han M, Chen CB, Zheng YT, He XS, Yuan XP. MicroRNA-34a regulates liver regeneration and the development of liver cancer in rats by targeting Notch signaling pathway. Oncotarget. 2017; 8:13264–76. https://doi.org/10.18632/oncotarget.14807 [PubMed]
  • 69. Hui F, Qin X, Zhang Q, Li R, Liu M, Ren T, Zhao M, Zhao Q. Alpinia oxyphylla oil induces apoptosis of hepatocellular carcinoma cells via PI3K/Akt pathway in vitro and in vivo. Biomed Pharmacother. 2019; 109:2365–74. https://doi.org/10.1016/j.biopha.2018.11.124 [PubMed]