Upregulation of circ_0002483 was associated with poor survival in patients with LUAC
The GEO dataset (GSE101684) was used to screen the differentially-expressed circRNAs between LUAC and non-cancerous tissue samples and top 10 upregulated or downregulated circRNAs were identified according to the FC > 2 and P < 0.01, of which hsa_circ_0002483 harbored a remarkable elevation in LUAC (Figure 1A). Then, FISH analysis validated that circ_0002483 expression levels were increased in pair-matched LUAC tissue samples as compared with the non-cancerous tissues (Figure 1B; n = 80, P = 0.01). The similar result was shown in LUAC patients with stage III-IV (n = 42) as compared with those with stage I-II (n = 38) (Figure 1C; P = 0.02).
In addition, we analyzed the association of circ_002483 with the clinicopathological characteristics in patients with LUAC and found that elevated expression of circ_002483 was related with TNM stage (Table 1; P = 0.048) rather than other parameters (P > 0.05) in LUAC. On the basis of the cutoff value of circ_002483, we divided the cases into circ_002483-high (n = 43) and circ_002483-low groups (n = 37). Kaplan-Meier analysis uncovered that the patients with circ_002483-high group possessed a poorer survival as compared with those with circ_002483-low expression (Figure 1D; P = 0.014). Univariate and Multivariate analysis unveiled that TNM stage and pathological stage rather than circ_002483 expression are independent prognostic factors of poor survival in patients with LUAC (Tables S2).
Identification of a novel circ_0002483 and its cellular localization
It is indicated that hsa_circ_0002483 (chr8:141874410-141900868) is originated from exon 2, 6 regions within protein tyrosine kinase 2 (PTK2) locus, and the spliced sequence is about 482 bp (Figure 2A). The back-spliced junction of circ_0002483 was amplified by divergent primer and circ_0002483 could be amplified from only cDNA but not gDNA in A549 cells (Figure 2B). Relative to linear PTK2, circ_0002483 produced a resistance to RNase R treatment in A549 and PC9 cell lines, and circ_0002483 harbored a loop structure in LUAC cells (Figure 2C). After A549 cells were treated by Actinomycin D, qRT-PCR displayed that circ_0002483 was more stable than PTK2 (Figure 2D). qRT-PCR and FISH analysis revealed that circ_0002483 was mainly localized in the cytoplasm of LUAC cells and tissues (Figure 2E, F).
Circ_0002483 facilitated proliferation, colony formation and invasion in vitro
Upregulation of circ_0002483 in LUAC suggested that it might be an oncogenic factor. To test this hypothesis, we evaluated the function of circ_0002483 in A549 and PC9 cells and constructed the overexpression vectors and the siRNA against circ_0002483. As indicated in Figure 3A, the overexpression vectors and the siRNA of circ_0002483 could markedly increase and decrease the expression of circ_0002483 in A549 and PC9 cell lines, respectively. MTT assay was conducted to chart the growth curve which indicated that upregulation of circ_0002483 obviously increased the proliferation viability of A549 and PC9 cells, while downregulation of circ_0002483 repressed cell growth (Figure 3B). Colony formation assay further verified that the cell colony number of A549 and PC9 was remarkably elevated by upregulation of circ_0002483 and reduced by downregulation of circ_0002483 (Figure 3C, D). In addition, the effects of circ_0002483 on LUAC cell invasion were estimated by Transwell assay, which indicated that overexpression of circ_0002483 increased invasive capabilities of A549 and PC9 cells, while knockdown of circ_0002483 harbored the opposite effects (Figure 3E, F).
Circ_0002483 was negatively associated with miR-125a-3p expression in LUAC tissue samples
To demonstrate the underlying mechanisms of circ_0002483 in LUAC cells, we identified the binding sites of 5 miRNAs with circ_0002483 3’UTR by Circular RNA profiling and Interactome (Figure 4A). TCGA cohort showed that, relative to the other 4 miRNAs, only miR-125a-3p possessed a decreased expression in paired (n = 39) and unpaired LUAC tissue samples (n = 448, Figure 4B). The downregulation and cytoplastic location of miR-125a-3p in pair-matched LUAC tissues were further confirmed by FISH analysis (n = 80, P = 0.025; Figure 4C-E). Pearson correlation analysis demonstrated that circ_0002483 harbored a negative correlation with miR-125a-3p expression in LUAC tissue samples (P = 0.036; Figure 4E).
Then, we analyzed the association of miR-125a-3p with clinicopathological characteristics in patients with LUAC (Table S3) and found that miR-125a-3p expression was related with gender (P = 0.012) and lymph node metastasis (P = 0.013) in patients with LUAC. However, the patients with miR-125a-3p-low expression showed no difference in poor survival and tumor recurrence as compared with those with miR-378a-3p-high expression (Figure S1).
Circ_0002483 acted as a sponge of miR-125a-3p
The specific binding sites of miR-125a-3p with WT or Mut circ_0002483 3’UTR were demonstrated in Figure 5A. To determine whether circ_0002483 could bind with miR-125a-3p, we co-transfected A549 and PC9 cells with WT or Mut circ_0002483 3’UTR reporter vectors and miR-125a-3p mimics or inhibitor, and found that miR-125a-3p mimics decreased the luciferase activities of WT circ_0002483 3’UTR in A549 and PC9 cells, while miR-125a-3p inhibitor enhanced such an effect (Figure 5B). But, miR-125a-3p exhibited no effects on those of Mut circ_0002483 3’UTR. Further qRT-PCR analysis indicated that upregulation of circ_0002483 reduced the expression of miR-125a-3p, and downregulation of circ_0002483 increased its expression (Figure 5C, D), whereas miR-125a-3p exhibited no effects on circ_0002483 expression in A549 and PC9 cells (Figure S2). Moreover, RIP assay was executed for Ago2 protein in A549 and PC9 cells, and the endogenous expression of circ_0002483 and miR-125a-3p pulled down from Ago2-expressed A549 and PC9 cells, indicated by qRT-PCR analysis was mainly enriched in Ago2 pellet relative to the input control (Figure 5E, F). After co-transfection with circ_0002483 and miR-125a-3p mimics or si-circ_0002483 and miR-125a-3p inhibitors into A549 and PC9 cells for 5 days, we found that miR-125a-3p suppressed the proliferation viability and attenuated circ_0002483-caused cell proliferation, while miR-125a-3p inhibitors displayed the opposite effects (Figure 5G).
MiR-125a-3p reversed circ_0002483-caused upregulation of CCL4-CCR5 axis
The targets of miR-125a-3p were identified by Targetscan7.1 and mirPathv.3 and the binding sites of miR-125a-3p with 3’UTR of CCL4 were demonstrated in Figure 6A. To determine whether miR-125a-3p could bind with 3’UTR of CCL4, we co-transfected A549 and PC9 cells with WT or Mut CCL4 3’UTR reporter vectors and miR-125a-3p mimics, and found that miR-125a-3p mimics decreased the luciferase activities of WT CCL4 3’UTR in A549 and PC9 cells, however, miR-125a-3p exhibited no effects on those of Mut CCL4 3’UTR (Figure 6B). TCGA cohort indicated that CCL4 and CCR5 were upregulated in 515 LUAC tissue samples (Figure 6C, D; P < 0.0001) and miR-125a-3p harbored a negative correlation with both of them in LUAC (Figure 6E; P < 0.05) . In addition, we found that high expression of CCL4 harbored no association with clinicopathological features in patients with LUAC with (Table S4), and the patients with CCL4-high expression possessed no difference in poor survival and tumor recurrence relative to those with CCL4-low expression (Figure S3). Moreover, qRT-PCR and Western blot analyses demonstrated that miR-125a-3p inhibited the expression of CCL4 and CCR5 and reversed circ_0002483-caused upregulation of CCL4-CCR5 axis in A549 cells (Figure 6F, G).
Knockdown of circ_0002483 repressed cell growth in vivo
To ascertain the effects of circ_0002483 on LUAC tumor growth in vivo, stably transfected A549 cells infected with si-NC or si-circ_0002484 were constructed and subcutaneously injected into 6-week old female nude mice (Figure 7A). Moreover, the tumor volume and weight in the si-circ_0002384 group were smaller and lighter than the si-NC group (Figure 7B, C). Hematoxylin and eosin (HE) staining and IHC of Ki-67 indicated that knockdown of circ_0002483 inhibited cell proliferation as compared with the si-NC group (Figure 7D).