Circular RNA EPB41 expression predicts unfavorable prognoses in NSCLC by regulating miR-486-3p/eIF5A axis-mediated stemness
Cancer Cell International volume 22, Article number: 219 (2022)
Dysregulation of circular RNAs (circRNAs) has recently been found to play an important role in the progression and development of cancers such as non-small cell lung cancer (NSCLC). Yet the functions of many circRNAs in NSCLC remain unclear. In this study, the circRNA expression profiles in NSCLC tumor tissues and adjacent non-tumorous tissues were detected by high-throughput sequencing. Bioinformatics analyses, the dual-luciferase reporter system, fluorescence in situ hybridization (FISH) and miRNA/mRNA high-throughput sequencing were used to identify circ-EPB41 and its downstream target. The subcutaneous tumor/caudal vein transfer mouse model was used for tumor growth and invasion analysis. The results show that the circ-EPB41 was upregulated in NSCLC tissues and cell lines. Increased circ-EPB41 expression in NSCLC was significantly correlated with malignant characteristics, and positive to post-surgical overall survival of NSCLC patients. Reduced circ-EPB41 expression in NSCLC decreased cell proliferation and invasion in both in vitro and in vivo experiments. The miRNA/mRNA high-throughput sequencing suggested that downregulation of circ-EPB41 promoted microRNA (miR)-486-3p and suppressed eukaryotic translation initiation factor 5A (eIF5A) expression. Luciferase reporter experiments confirmed that miR-486-3p/eIF5A were downstream targets of circ-EPB41. In addition, we also found that downregulation of circ-EPB41 suppressed self-renewal and decreased expression of stemness markers SOX2, OCT-4, Nanog and CD133 by sponging miR-486-3p to enhance eIF5A expression. Taken togeter, these data revealed the important role of circ-EPB41 in regulating NSCLC cell invasion and proliferation by modifying miR-486-3p/eIF5A axis-mediated stemness. We believe our study provides a novel perspective regarding the role of circRNAs in NSCLC progression.
Lung cancer (LC) is among the most prevalent cancers, and its prevalence is a major factor affecting cancer mortality worldwide. There were about 0.73 million newly diagnosed cases and 0.61 million deaths because of LC in 2015 in China [1, 2]. Based on the annual predictions of the American Cancer Society, we can expect about 1.762 million new cancer patients and 606,880 cancer deaths in the US in 2019, of which LC will be among the most prevalent causes of both mortality and incidence . Approximately 80–85% all human LCs are non-small cell lung cancer (NSCLC), consisting of two main subtypes, lung squamous cell carcinoma and lung adenocarcinoma [4, 5]. The 5-year relative survival rate is 57% for patients with stage I disease and declines to 4% for patients with stage IV disease . Because of the lack of effective diagnostic methods and markers, the majority of NSCLC patients are diagnosed in the advanced stage, which is characterized by poor prognosis and distant metastasis . Therefore, a better understanding of the molecular mechanisms that underlie NSCLC progression is indispensable for developing more effective therapies.
Circular RNAs (circRNAs) belong to a class of noncoding and protein-coding RNAs localized to the cytoplasm which originate from both single and multiple exons [7, 8]. Many studies have shown that circRNAs are involved in various pathological and physiological processes including cancer progression and tumorigenesis [9, 10]. For instance, circFGFR1 may interact directly with miR-381-3p and function as a miRNA sponge to upregulate expression of miR-381-3p target gene C-X-C motif chemokine receptor 4, which enhances NSCLC progression as well as resistance to anti-programmed cell death 1-based therapy . Circ_0000003 promotes NSCLC cell metastasis and proliferation through miR-338-3p/insulin receptor substrate 2 . Circ_0002483 has been verified to inhibit NSCLC progression in vivo and in vitro, and increases NSCLC cell sensitivity to Taxol through sponging of miR-182-5p to inhibit expression of FOXO1, FOXO3 and GRB2 mRNAs . Our previous studies have found that circARHGAP10 suppresses NSCLC progression by functioning like a miR-150-5p sponge to enhance GLUT1 expression . In NSCLC, previous investigations have discovered the existence of various circRNAs such as circ-ZKSCAN1, circ-CMPK1, circ-HIPK3 and circ-Foxo3 [15,16,17,18]. While the circRNA function and mechanism in NSCLC are not fully explored.
Current research has examined circRNA profiles expressed in human NSCLC tissues and shown that hsa_circ_0000042 (circ-EPB41) is a conserved circRNA, which is significantly upregulated in NSCLC tissues. In addition, circ-EPB41 expression is closely associated with NSCLC patient prognosis. We report here that circ-EPB41 promoted NSCLC progression by regulating miR-486-3p and eukaryotic translation initiation factor 5A (eIF5A)-mediated stemness. Thus, circ-EPB41 may serve as a novel marker for NSCLC prognosis and may be a promising therapeutic target.
Materials and methods
In total, we collected 90 fresh NSCLC tissues with paired adjacent non-cancerous lung tissues after obtaining informed patient consent at Renji Hospital affiliated with Shanghai Jiaotong University, China. We evaluated pathological and histological diagnostics of NSCLC based upon the Revised International System for Staging Lung Cancer. Patients received no radiotherapy or chemotherapy prior to tissue sampling. We snap-froze samples in liquid nitrogen and stored them at − 80 °C before RNA extraction. The Ethics Committee of Shanghai University of Medicine and Health Sciences approved the study.
Strand-specific RNA-Seq library construction and high-throughput RNA-Seq
We extracted total RNA from three paired NSCLC tissues as well as adjacent non-cancerous lung tissues using TRIzol Reagent (Invitrogen, Carlsbad, CA, USA). We subjected 3 μg of total RNA from each sample to the VAHTS Total RNA-Seq (H/M/R) Library Prep Kit from Illumina (Vazyme Biotech Co., Ltd, Nanjing, China) to remove ribosomal RNA and different classes of RNA such as non-coding RNA and mRNA. We treated purified RNA with RNase R (Epicentre Technologies, Madison, WI, USA; 40 U, 37 °C, 3 h), followed by TRIzol purification. We prepared RNA-Seq libraries using the KAPA Stranded RNA-Seq Library Prep Kit (Roche, Basel, Switzerland) and subjected them to deep sequencing with the Illumina HiSeq 4000 at Aksomics, Inc. (Shanghai, China; accession code: H1712024).
For miRNA and mRNA analysis, A549 cells transfected with siRNA against circ-EPB41 or a negative control (NC) vector were used for high-throughput RNA-Seq as previously described (accession code: H1712024).
Cell lines and culture
We obtained the human normal lung epithelial cell line BEAS-2B and the NSCLC cell lines PC9, A549, H1975 and H1650 from the Cell Bank of the Chinese Academy of Sciences (Shanghai, China) and cultivated them in Dulbecco’s Modified Eagle’s Medium (DMEM; Life Technologies, Carlsbad, CA, USA) supplemented with 100 IU/mL penicillin, 100 μg/mL streptomycin and 10% fetal bovine serum (FBS; Invitrogen) at 37 °C in a humidified atmosphere with 5% CO2.
Fluorescence in situ hybridization (FISH)
We obtained specific probes for circ-EPB41 (Dig-5ʹ-CATACTGCTCTTGCTCCATCGAGGCTCC-3ʹ-Dig) and miR-486-3p (Dig-5ʹ-CGGGGCAGCUCAGUACAGGAUA-3ʹ-Dig) from Geneseed Biotech (Guangzhou, China). Cells were immunostained with Cy3-conjugated anti-digoxin and FITC-conjugated anti-biotin antibodies (Jackson ImmunoResearch Inc., West Grove, PA, USA). We counterstained nuclei with 4,6-diamidino-2-phenylindole (DAPI). Staining was visualized with a Zeiss LSM 700 confocal microscope (Carl Zeiss, Oberkochen, Germany).
We imputed circRNA/miRNA target genes through the Interactome (https://circinteractome.nia.nih.gov/), miRanda (http://www.microrna.org/microrna/home.domiRanda) and circBank (http://www.circbank.cn/) websites. We calculated interactive correlations between miR-486-3p and eIF5A via the Targetscan website.
Cell cycle detection
We dissociated cells in the logarithmic growth phase with 0.25% trypsin (Invitrogen, Carlsbad, CA, USA), resuspended them in PBS and fixed them in 70% ice-cold ethyl alcohol overnight at 4 °C. Afterwards, we centrifuged cells at 1200g for 5 min, resuspended them in 50 μL RNase A (Invitrogen, Carlsbad, CA, USA) and incubated them at 37 °C for 30 min. Then, we added 400 μL propidium iodide to the suspension for 30 min, followed by flow cytometric detection (BD Biosciences, San Jose, CA, USA).
Plasmid construction and stable transfection
We synthesized human eIF5A and circ-EPB41 cDNA and then cloned it into the luciferase-labelled pcDNA3.1 vector (Invitrogen). We synthesized wild-type (WT) and mutant (MUT) circ-EPB41 and eIF5A cDNAs and cloned them into pZW1 vectors (Shanghai Institutes for Biological Sciences, Shanghai, China). We also transfected H1650 and A549 cells with the aforementioned plasmids.
To further validate the effects of circ-EPB41 in the in vivo experiments, a lentiviral-based small hairpin RNA (shRNA) (lentiviral vector: CMV-MCS-EF1α-GFP-PURO; Novagen) targeting circ-EPB41 was constructed. GFP detection was performed 72 h after infection into A549, and at a green fluorescence > 95%, the transfection was considered successful. For in vivo metastatic detection, pLVX-Luc2-P2A-AcGFP1-puro lentiviral vector (Novagen) was used for A549 infection and it generated the luc-A549 cell line.
Total RNA isolation and the quantitative reverse transcription-polymerase chain reaction (RT-qPCR)
We isolated total RNA from tumor tissues and cells with TRIzol reagent following standard procedures. We examined the concentration and purity of RNA samples spectrophotometrically by capturing the absorbance at 230 nm, 260 nm and 280 nm with a NanoDrop ND-1000 (Thermo Fisher Scientific, Wilmington, DE, USA). We deemed OD260/OD280 ratios between 1.8 and 2.1 as acceptable, and also deemed OD260/OD230 ratios > 1.8 as acceptable.
We reverse transcribed total RNA before RT-qPCR detection. We obtained primers specific for circ-EPB41, eIF5A and miR-486-3p from GenePharma (Shanghai, China). We performed RT-qPCR using an AB7300 thermo recycler (Applied Biosystems, Carlsbad, CA, USA) with the primers listed below and TaqMan Universal PCR Master Mix. We utilized GAPDH as a reference gene for circRNAs and mRNAs and U6 RNA as an internal control for miRNA expression level. We quantified gene expression via the 2−ΔΔCt method. The primers utilized to assay circ-EPB41 expression were forward, 5′-AGGATCCAAATTTCGATACAGTGGC-3′ and reverse, 5′-ATTTCTTAGCTGCTCGGTAACTGGG-3′. The primers for miR-486-3p were forward, 5′-GCGGGGCAGCTCAGTA-3′ and reverse, 5′-GTTGGCTCTGGTGCAGGGTCCGAGGTATTCGCACCAGAGCCAACATCCTG-3′. The eIF5A primers were forward, 5′-GAGATGCAGGGGCCTCAGCCACC-3′ and reverse, 5′-GTGATAGGTACCCATCCTGGATG-3′. U6 primers were forward, 5′-CTCGCTTCGGCAGCACA-3′ and reverse, 5′-AACGCTTCACGAATTTGCGT-3′. The GAPDH primers were forward, 5′-GCACCGTCAAGGCTGAGAAC-3′ and reverse, 5′-GGATCTCGCTCCTGGAAGATG-3′.
RNA interference and overexpression
The miR-486-3p inhibitors and siRNA against circ-EPB41 (sicirc-EPB41) were purchased from GenePharma. Transfections were performed following the supplier's protocol. In brief, we transferred cells to 6-well culture plates and transfected them using Lipofectamine 2000 (Invitrogen). To achieve eIF5A overexpression, we transfected the pCDNA3.0 vector described above. For xenograft experiments, a lentiviral-mediated circ-EPB41-silencing vector (sh-circ-EPB41) was transfected into A549 cells.
Tumor sphere formation assay
We harvested A549 and H1650 cells and resuspended them as single cells in serum-free medium. After precise cell counting, we added 200 cells/well in 200 μL of serum-free medium to a 96-well plate, 10 wells per group. We changed medium every 2 days. We took images of five randomly selected regions in each group of wells with a camera-equipped microplate reader (Leica, Wetzlar, Germany). We calculated sphere percentage by the number of spheres/200.
Cell proliferation assay
We employed Cell Counting Kit-8 assays (CCK-8; Gibco, Logan, Utah, USA) to quantify cellular proliferation. We seeded transfected cells into 96-well plates at a density of 5,000 cells/well in triplicate. We detected cell viability with the CCK-8 system at 0, 24, 48 and 72 h after seeding, following standard procedure.
For colony formation assays, we seeded transfected cells into 6-well plates at a density of 2000 cells/well and maintained them in DMEM containing 10% FBS for 10 days. We imaged colonies and counted them after fixing and staining them.
For invasion assays, we placed Transwell assay inserts (Millipore, Billerica, MA, USA) into 24-well plates. The membrane in the upper Transwell chamber was coated with Matrigel (BD Biosciences). We placed 500 μL of DMEM containing 10% FBS into the bottom chamber and seeded 10,000 cells in 200 μL serum-free DMEM in the upper chamber. After 1–2 days, we utilized methanol to fix cells on the membrane and stained them with Crystal Violet. Finally, we observed the cells via a microscope (Leica).
Western blot assay
We extracted protein from tissues or cells with RIPA lysis buffer and performed western blot assays as previously described . We obtained primary antibodies against SOX2, OCT-4, Nanog, CD133 and GAPDH from Cell Signaling Technology (Beverly, MA, USA) and stained protein blots following standard procedures. We visualized immunoreactivity with a chemiluminescence detection kit (Western Blotting Substrate, Donghuan Biotech, China).
We fixed cells or tumor tissues in 4% paraformaldehyde for 20 min at room temperature (RT). We then incubated the fixed cells or tissues with 0.5% Triton-X-100 in PBS for five min at RT and blocked them with 5% bovine serum albumin in PBS for 1 h at RT. The cells were incubated with primary antibodies against SOX2 (1:200; Abcam, Cambridge, MA, USA) at 4 °C overnight, followed by incubation with fluorochrome-conjugated secondary antibody (1:200; Abcam) and Nestin primary antibody conjugated to FITC (1:200; Millipore) at RT for 1 h. We stained nuclei in DAPI for 15 min at RT and visualized staining via fluorescence microscope.
The 5-ethynyl-2′-deoxyuridine (EdU) assay
We purchased an EdU assay kit (RiboBio, China) to measure cell proliferation and DNA synthesis. We seeded 10,000 A549 or H1650 cells in 96-well plates overnight, then added Edu solution (25 μM) to the wells for 1 day. We applied 4% formalin to fix the cells at RT for 2 h and permeabilized the cells with 0.5% TritonX-100 for 10 min. We then added 200 μL of Apollo reaction solution for 30 min to stain EdU and 200 μL of DAPI to stain the nuclei. Lastly, we utilized a Nikon microscope (Nikon, Japan) to detect cell proliferation and DNA synthesis, reflected by blue and red signals, respectively.
Dual-luciferase reporter assay
We constructed WT/mut-circ-EPB41 or WT/mut-eIF5A 3'UTR fragments and inserted them downstream into a pMIR-REPORT plasmid luciferase reporter gene (H306; Obio Technology, Shanghai, China). We employed Lipofectamine 2000 to transfect the reporter plasmid into A549 cells and co-transfected the miR-486-3p mimic with the reporter gene into HEK293T cells. Afterwards, we employed the Dual Luciferase Reporter System Kit (E1910; Promega, Madison, WI, USA) to measure firefly and Renilla luciferase activity.
To detect the role of circ-EPB41 in a LC metastasis model, we intravenously injected 1 × 106 stable lentiviral-mediated circ-EPB41-silenced A549 cells (sh-circ-EPB41) or A549 NC cells into male nude mice (Chinese Science Academy, Shanghai, China) via the tail vein. After 1 month, we analyzed A549 cell metastasis by bioluminescence imaging following an intravenous luciferin injection (150 mg luciferin/kg body weight) into the tails. The numbers of metastatic foci in lung tissues were caculation according to the HE staining.
For xenograft assays, we injected 1 × 106 modified (circ-EPB41) or control (WT) A549 cells subcutaneously into the right flank of male nude mice. We calculated tumor volumes (length × width2 × 0.5) at the indicated time points and excised the tumors 4 weeks after injection.
The mice were killed by cervical dislocation and the tumors were collected. The animal protocols complied with the rule of the ethics committee of Shanghai University of Medicine and Health Sciences.
We fixed tumor tissue samples in 10% formalin and embedded them in paraffin. We stained sections (5-μm thick) with Ki67 to evaluate proliferation. We examined sections with an Axiophot light microscope and imaged them by digital camera.
We assessed differences between any two groups by paired/unpaired two-tailed t-tests. We used Pearson’s correlation test to detect correlations between groups. We expressed data as means ± SD. A P-value < 0.05 was regarded as significant. We performed statistical analysis using the GraphPad Prism package (GraphPad Inc., San Diego, CA, USA).
High circ-EPB41 expression in NSCLC is correlated with poor prognosis.
To verify the correlation between abnormal circRNA expression and NSCLC progression, we performed RNA-seq analyses on ribosomal RNA-depleted total RNA from three NSCLC and three adjacent normal tissue samples. The data indicated that 22 circRNAs were upregulated and 23 circRNAs were downregulated (Additional file 1; Fig. 1A). To validate the RNA-seq results, we selected six upregulated circRNAs (hsa_circ_0001721 [circ-CDK14], hsa_circ_0000130 [SNX27], hsa_circ_0001747 [circ-MKLN1], hsa_circ_0001869 [circ-ZCCHC6], hsa_circ_0008992 [circ-PAK1] and hsa_circ_0000042 [circ-EPB41]) for RT-qPCR analysis using ten paired NSCLC tissues. Results were consistent with RNA-seq data and we also found that circ-EPB41 expression was significantly enhanced in NSCLC tissues comparing with non-tumorous tissues (increased 15 times, the median value was 3.5; Fig. 1B).
Among these specific candidates, circ-EPB41 with higher expression levels was examined in more detail. We found that circ-EPB41 originated by circularization of two exons from the EPB41 gene, which is located at chr1:29362337-29365938. EPB41 consists of 3601 bp and the spliced mature circRNA is 271 bp (Fig. 1C). In order to further identify the subcellular localization of circ-EPB41, specific probes for circ-EPB41 were constructed across the back-splicing site. FISH assays detected that circ-EPB41 was localized mainly to the cytoplasm, and circ-EPB41 expression in NSCLC tissues was significantly increased compared to matched non-tumorous tissues (Fig. 1D). RT-qPCR detection also showed that circ-EPB41 expression was increased in NSCLC tissues using 90 paired NSCLC tissues (Fig. 1E). Furthermore, RT-qPCR detection showed that the expression of circ-EPB4 in A549, PC9, H1650 and H1975 cells was increased compared to normal lung epithelial cells (BEAS-2B; Fig. 1F). A549 and H1650 cells had the highest expression of circ-EPB41 and were selected for further study. Both RT-qPCR and FISH assays confirmed that circ-EPB41 was a circular RNA, because both promer and FISH probes were constructed across the back-splicing site.
We divided 90 paired NSCLC tissues into relatively high-expressing (n = 46) and relatively low-expressing (n = 44) tissues compared to normal circ-EPB41 expression. There was no correlation between circ-EPB41 expression and demographic factors such as sex or patient age (≤ 60 years vs > 60 years). However, high expression of circ-EPB41 was positively related to lymph node metastasis (negative vs positive), TNM stage (I/II or III/IV vs high) and tumor size (≤ 3 cm vs > 3 cm; Table 1). Also, Gehan-Breslow-Wilcoxon test survival curves demonstrated that NSCLC patients with high circ-EPB41 expression exhibited poor overall survival (Fig. 1G). These data suggested that the circ-EPB41 expression has an important function in NSCLC progression.
The circ-EPB41 downregulation suppresses NSCLC invasion and proliferation in vivo and in vitro
To identify circ-EPB41’s role in NSCLC proliferation, we used H1650 and A549 cells. RT-qPCR analyses suggested that circ-EPB41 expression significantly decreased after transfection with sicirc-EPB41 in H1650 and A549 cells compared with the NC group (Fig. 2A). Cell cycle analysis by flow cytometry conducted on A549 cells illustrated that circ-EPB41 knockdown resulted in mild cell cycle arrest at the G0/G1 phase in nearly 14% of cells (Fig. 2B). CCK8 detection (Fig. 2C, D) and colony formation assays (Fig. 2E, F) demonstrated that circ-EPB41 silencing suppressed cell proliferation of H1650 and A549 cells. We then examined the effect of circ-EPB41 knockdown on tumor formation in a stable lentiviral strain (sh-circEPB41) or in sh-NC A549 cells. Xenograft results confirmed that circ-EPB4 knockdown suppressed tumor growth, as determined by tumor volume and immunohistochemistry (Fig. 2G) using Ki67 staining (Fig. 2H), compared with the NC group.
Transwell assays demonstrated that circ-EPB41 downregulation suppressed NSCLC invasion in H1650 and A549 cells (Fig. 3A, B). We also employed a small animal live imaging system to obtain fluorescence images from nude mice 30 days after tail-vein inoculation of circ-EPB41-silenced and control cells (Fig. 3C). The fluorescence intensity in the circ-EPB41-silenced group was consistently weaker than the NC group. HE staining analysis show that circ-EPB41 silence decreased the pulmonary metastasis ability by decreased the numbers of metastatic foci in lung tissues after HE staining analysis (Fig. 3D–F).
The miR-486-3p and eIF5A are downstream targets of circ-EPB41
Next, we investigated how circ-EPB41 regulation affects NSCLC proliferation and invasion. Previous studies have found that circRNA can act as a miRNA sponge. Given that circ-EPB41 is mainly localized to the cytoplasm, we hypothesized that circ-EPB41 might regulate tumor biological behavior by sponging miRNAs. We employed Interactome, miRanda and circBank to predict miRNAs targeted by circ-EPB41. The results showed that three miRNAs could be sponged by circ-EPB41 (Fig. 4A). Bioinformatics analysis further found that miR-587, miR-486-3p and miR-595 were circ-EPB41 targets (Fig. 4B; Additional file 2). High-throughput sequencing for miRNA expression using A549 cells found that downregulation of circ-EPB41 promoted miR-486-3p expression (Additional file 3). To verify this finding, we constructed a dual-luciferase reporter system by inserting the circ-EPB41 sequence into the psiCHECK2 plasmid (WT). The data showed that when co-transfected with WT or NC and miRNAs, only mimic miR-486-3p significantly decreased luciferase activity (Fig. 4C), suggesting that miR-486-3p was the direct target of circ-EPB41. We then cloned circ-EPB41 sequences containing mutated miR-486-3p binding sites into the psiCHECK2 plasmid and observed no obvious alterations in luciferase activity after co-transfection with MUT and corresponding miR-486-3p mimic (Fig. 4D). FISH detection showed that circ-EPB41 and miR-486-3p interacted and were co-localized in the cytoplasm of both A549 and H1650 cells (Fig. 4E). RT-qPCR detection also found that miR-486-3p expression was increased NSCLC cancer tissues when compared with tumor-adjacent tissue (Fig. 4F). High-throughput sequencing of A549 mRNA expression found that downregulation of circ-EPB41 decreased eIF5A expression (Fig. 4G). To verify that eIF5A was a miR-486-3p target, WT and MUT eIF5A 3'UTR sequences were cloned into the psiCHECK2 plasmid. The data revealed that while miR-486-3p mimic and reporter plasmid co-transfection visibly suppressed luciferase activity, miR-486-3p mimic and mutated eIF5A 3'UTR vector co-transfection did not have a significant effect on luciferase activity (Fig. 4H). Therefore, the data demonstrated that miR-486-3p directly targeted eIF5A.
The circ-EPB41 promotes NSCLC progression by regulating miR-486-3p/eIF5A axis-mediated stemness
Previous studies have found that eIF5A can regulate cancer stem cells (CSCs) [20,21,22,23]. In the present study we examined if circ-EPB41 can regulate CSCs in NSCLC. Immunofluorescence detection of OCT-4 staining found that circ-EPB41 silencing suppressed stemness marker OCT-4 expression in tumor tissues (Fig. 5A, B). Western blot detection also found that expression of stemness markers SOX2, OCT-4, Nanog and CD133 was decreased after circ-EPB41 silencing in tumor tissues (Fig. 5C, D), suggesting that circ-EPB41 promoted the progression of NSCLC by regulating stemness.
In order to identify the interaction among circ-EPB41, miR-486-3p and eIF5A that regulates stemness, we transfected H1650 and A549 cells with negative vector (NC), sicirc-EPB41, miR-486-3p inhibitor and eIF5A overexpression vector singly or in combination. RT-qPCR detection showed that downregulation of circ-EPB41 suppressed circ-EPB41 expression, while downregulation of miR-486-3p or overexpression of eIF5A could not restore circ-EPB41 expression after circ-EPB41 silencing in A549 and H1650 cells (Fig. 6A, B). The results also found that downregulation of circ-EPB41 promoted miR-486-3p expression, but after treatment with miR-486-3p inhibitor, miR-486-3p expression was suppressed in A549 and H1650 cells. Upregulation of eIF5A could not suppress miR-486-3p expression after circ-EPB41 was silenced (Fig. 6C, D). RT-qPCR detection of eIF5A indicated that circ-EPB41 silencing suppressed eIF5A expression. Meanwhile, downregulation of miR-486-3p restored eIF5A expression, and eIF5A overexpression increased eIF5A expression in H1650 and A549 cells (Fig. 6E, F).
Immunofluorescence detection with OCT-4 staining found that silencing of circ-EPB41 suppressed OCT-4 expression in both A549 and H1650 cells, while downregulation of miR-486-3p or overexpression eIF5A restored stemness relative protein OCT-4 expression after circ-EPB41 silencing (Fig. 6G, H). We then conducted tumor sphere formation assays. These results indicated that downregulation of miR-486-3p or overexpression of eIF5A restored the sphere percentage in H1650 and A549 cells (Fig. 7A–C), suggesting that circ-EPB41 expression enhanced NSCLC progression by regulating miR-486-3p/eIF5A axis-mediated stemness.
Downregulating miR-486-3p or overexpressing eIF5A restores cell proliferation and invasion after circ-EPB41 silencing
Our study also found using colony formation assays and Edu assays in H1650 and A549 cells that downregulation of miR-486-3p or overexpression of eIF5A renewed proliferative ability after eIF5A was silenced (Fig. 8A–F). Transwell invasion analysis also demonstrated that downregulating miR-486-3p or overexpressing eIF5A restored invasive ability after eIF5A was silenced in A549 and H1650 cells (Fig. 8G–I).
During last few decades, various studies have discovered abnormal non-coding RNA expression profiles in a number of cancers. Investigations have determined that many miRNAs and circRNAs have vital functions in modulating cancer metastasis and tumor growth [24, 25], although the mechanisms by which circRNAs participate in cancer progression and development are still unknown [26, 27]. To date, a few circRNAs have been shown to act critically in NSCLC [11, 13, 28,29,30]. Here, we report that circ-EPB41 is an essential circRNA that is frequently upregulated in NSCLC tissues, and that high expression predicts poor prognosis for lymph node metastasis, high TNM stage and large tumor size. Downregulation of circ-EPB41 induced G2/M arrest in NSCLC cells, while circ-EPB41 silencing suppressed NSCLC proliferation and metastasis in in vitro and in vivo experiments, suggesting that circ-EPB41 functions in NSCLC progression.
Circ-EPB41 originates from two EPB41 gene exons. FISH detection indicated that circ-EPB41 is localized to cytoplasm. Additional published evidence implies that circRNAs function as miRNA sponges or as endogenous RNA (ceRNA) competitors to decrease functional dysregulation of miRNAs and their target genes, causing tumor proliferation and invasion, including in NSCLC [31, 32]. In the present study, A549 cells with or without circ-EPB41 silencing were used for high-throughput sequencing analysis. The results showed that downregulating circ-EPB41 resulted in a series of abnormally expressed mRNAs and miRNAs. Bioinformatics analysis found that miR-486-3p, miR-587 and miR-595 were the circ-EPB41 targets. Luciferase reporter assays and FISH analysis further confirmed that miR-486-3p was the direct target of circ-EPB41. Previous studies have confirmed that miR-486-3p expression was decreased in tumor tissues of LC patients . Increased miR-486-3p represses cells metastasis and proliferation in many tumors including cervical cancer , oral cancer  and squamous cell carcinoma , suggesting that circ-EPB41 overexpression also suppressed NSCLC cell metastasis and proliferation via adsorption of miR-486-3p.
Further study with high-throughput sequencing and bioinformatics analysis found that miR-486-3p could target the 3'UTR of eIF5A. Previous studies have also found that eIF5A overexpression was associated with poor survival in many tumors [37, 38]. In the current study, we discovered that downregulation of miR-486-3p or overexpression of eIF5A restored cell proliferation and invasion after circ-EPB41 silencing, suggesting that circ-EPB41 expression predicts unfavorable prognoses in NSCLC by regulating the miR-486-3p/eIF5A axis.
Previous studies have also found that eIF5A expression can influence tumor progression by regulating tumor stem cell differentiation . CSCs are tumor cells with principal characteristics of clonal tumor initiation capacity, self-renewal and clonal long-term repopulation potential. In general, stem cell marker expression, such as SOX2, OCT-4, Nanog and CD133, with transcriptional signatures particular to CSCs, correlate functionally with aggressive behavior and are closely predictive of patient overall survival. Clinical data indicate that CSCs may be indispensable therapy targets [39, 40]. In the current investigation, we discovered that downregulation of circ-EPB41 suppressed stemness of NSCLC in in vivo and in vitro experiments. Downregulation of miR-486-3p or overexpression of eIF5A restored CSC differentiation after circ-EPB41 silencing, which suggested that circ-EPB41 expression predicts unfavorable prognoses in NSCLC by regulating miR-486-3p/eIF5A axis-mediated stemness. However, the specific regulatory mechanism controlling stemness in NSCLC needs further study.
In summary, the present study demonstrated that circ-EPB41 was upregulated in NSCLC tissues compared with normal lung tissue, which was related to high proliferation and NSCLC cell invasive activity through sponging of miR-486-3p. Moreover, circ-EPB41 sponged miR-486–3 to promote stemness in NSCLC by enhancing eIF5A expression. Therefore, circ-EPB41 represents as a prognostic marker and a target for NSCLC.
Availability of data and materials
The data generated or analyzed during this study are included in this article, or if absent are available from the corresponding author upon reasonable request.
Short hairpin RNA
Non-small cell lung cancer
Quantitative reverse transcription-polymerase chain reaction
Fluorescence in situ hybridization
Eukaryotic translation initiation factor 5A
Cancer stem cell
Cell Counting Kit-8
Miller KD, Nogueira L, Mariotto AB, et al. Cancer treatment and survivorship statistics, 2019. CA Cancer J Clin. 2019;69:363–85.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2016. CA Cancer J Clin. 2016;66:7–30.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2019. CA Cancer J Clin. 2019;69:7–34.
Qin H, Wang F, Liu H, et al. New advances in immunotherapy for non-small cell lung cancer. Am J Transl Res. 2018;10:2234–45.
Lovly CM, Carbone DP. Lung cancer in 2010: one size does not fit all. Nat Rev Clin Oncol. 2011;8:68–70.
Goldstraw P, Chansky K, Crowley J, et al. The IASLC Lung Cancer Staging Project: proposals for revision of the TNM Stage Groupings in the Forthcoming (Eighth) Edition of the TNM Classification for Lung Cancer. J Thorac Oncol. 2016;11:39–51.
An integrated encyclopedia of DNA elements in the human genome. Nature. 2012; 489: 57–74.
Zhao J, Lee EE, Kim J, et al. Transforming activity of an oncoprotein-encoding circular RNA from human papillomavirus. Nat Commun. 2019;10:2300.
Zhang PF, Wei CY, Huang XY, et al. Circular RNA circTRIM33-12 acts as the sponge of MicroRNA-191 to suppress hepatocellular carcinoma progression. Mol Cancer. 2019;18:105.
Wang L, Long H, Zheng Q, Bo X, Xiao X, Li B. Circular RNA circRHOT1 promotes hepatocellular carcinoma progression by initiation of NR2F6 expression. Mol Cancer. 2019;18:119.
Zhang PF, Pei X, Li KS, et al. Circular RNA circFGFR1 promotes progression and anti-PD-1 resistance by sponging miR-381-3p in non-small cell lung cancer cells. Mol Cancer. 2019;18:179.
Li S, Niu X, Li H, Liang Y, Sun Z, Yan Y. Circ_0000003 promotes the proliferation and metastasis of non-small cell lung cancer cells via miR-338-3p/insulin receptor substrate 2. Cell Cycle. 2019;18:3525–39.
Li X, Yang B, Ren H, et al. Hsa_circ_0002483 inhibited the progression and enhanced the Taxol sensitivity of non-small cell lung cancer by targeting miR-182-5p. Cell Death Dis. 2019;10:953.
Jin M, Shi C, Yang C, Liu J, Huang G. Upregulated circRNA ARHGAP10 predicts an unfavorable prognosis in NSCLC through regulation of the miR-150-5p/GLUT-1 axis. Mol Ther Nucleic Acids. 2019;18:219–31.
Wang Y, Xu R, Zhang D, et al. Circ-ZKSCAN1 regulates FAM83A expression and inactivates MAPK signaling by targeting miR-330-5p to promote non-small cell lung cancer progression. Transl Lung Cancer Res. 2019;8:862–75.
Cui D, Qian R, Li Y. Circular RNA circ-CMPK1 contributes to cell proliferation of non-small cell lung cancer by elevating cyclin D1 via sponging miR-302e. Mol Genet Genomic Med. 2020;8: e999.
Lu H, Han X, Ren J, Ren K, Li Z, Sun Z. Circular RNA HIPK3 induces cell proliferation and inhibits apoptosis in non-small cell lung cancer through sponging miR-149. Cancer Biol Ther. 2020;21:113–21.
Shen Z, Zhou L, Zhang C, Xu J. Reduction of circular RNA Foxo3 promotes prostate cancer progression and chemoresistance to docetaxel. Cancer Lett. 2020;468:88–101.
Zhang X, Wang S, Wang H, et al. Circular RNA circNRIP1 acts as a microRNA-149-5p sponge to promote gastric cancer progression via the AKT1/mTOR pathway. Mol Cancer. 2019;18:20.
Strnadel J, Choi S, Fujimura K, et al. eIF5A-PEAK1 signaling regulates YAP1/TAZ protein expression and pancreatic cancer cell growth. Cancer Res. 2017;77:1997–2007.
Xu J, Su Q, Gao M, Liang Q, Li J, Chen X. Differential expression and effects of peroxiredoxin-6 on drug resistance and cancer stem cell-like properties in non-small cell lung cancer. Onco Targets Ther. 2019;12:10477–86.
Heng WS, Cheah SC. Chelerythrine chloride downregulates beta-catenin and inhibits stem cell properties of non-small cell lung carcinoma. Molecules. 2020;25:224.
Masciale V, Grisendi G, Banchelli F, et al. Correlating tumor-infiltrating lymphocytes and lung cancer stem cells: a cross-sectional study. Ann Transl Med. 2019;7:619.
Liang WC, Wong CW, Liang PP, et al. Translation of the circular RNA circbeta-catenin promotes liver cancer cell growth through activation of the Wnt pathway. Genome Biol. 2019;20:84.
Tan S, Gou Q, Pu W, et al. Circular RNA F-circEA produced from EML4-ALK fusion gene as a novel liquid biopsy biomarker for non-small cell lung cancer. Cell Res. 2018;28:693–5.
Wang Y, Liu J, Ma J, et al. Exosomal circRNAs: biogenesis, effect and application in human diseases. Mol Cancer. 2019;18:116.
Su M, Xiao Y, Ma J, et al. Circular RNAs in Cancer: emerging functions in hallmarks, stemness, resistance and roles as potential biomarkers. Mol Cancer. 2019;18:90.
Wu K, Liao X, Gong Y, et al. Circular RNA F-circSR derived from SLC34A2-ROS1 fusion gene promotes cell migration in non-small cell lung cancer. Mol Cancer. 2019;18:98.
Chen L, Nan A, Zhang N, et al. Circular RNA 100146 functions as an oncogene through direct binding to miR-361-3p and miR-615-5p in non-small cell lung cancer. Mol Cancer. 2019;18:13.
Tan S, Sun D, Pu W, et al. Circular RNA F-circEA-2a derived from EML4-ALK fusion gene promotes cell migration and invasion in non-small cell lung cancer. Mol Cancer. 2018;17:138.
Dong Y, Xu T, Zhong S, et al. Circ_0076305 regulates cisplatin resistance of non-small cell lung cancer via positively modulating STAT3 by sponging miR-296-5p. Life Sci. 2019;239: 116984.
Liu YT, Han XH, Xing PY, et al. Circular RNA profiling identified as a biomarker for predicting the efficacy of Gefitinib therapy for non-small cell lung cancer. J Thorac Dis. 2019;11:1779–87.
Xu Y, Wang Y, Yao A, et al. Low frequency magnetic fields induce autophagy-associated cell death in lung cancer through miR-486-mediated Inhibition of Akt/mTOR signaling pathway. Sci Rep. 2017;7:11776.
Ye H, Yu X, Xia J, Tang X, Tang L, Chen F. MiR-486-3p targeting ECM1 represses cell proliferation and metastasis in cervical cancer. Biomed Pharmacother. 2016;80:109–14.
Chou ST, Peng HY, Mo KC, et al. MicroRNA-486-3p functions as a tumor suppressor in oral cancer by targeting DDR1. J Exp Clin Cancer Res. 2019;38:281.
Wang JX, Liu Y, Jia XJ, et al. Upregulation of circFLNA contributes to laryngeal squamous cell carcinoma migration by circFLNA-miR-486-3p-FLNA axis. Cancer Cell Int. 2019;19:196.
Wang Z, Jiang J, Qin T, Xiao Y, Han L. EIF5A regulates proliferation and chemoresistance in pancreatic cancer through the sHH signalling pathway. J Cell Mol Med. 2019;23:2678–88.
Liu Z, Teng L, Gao L, Wang H, Su Y, Li J. The role of eukaryotic translation initiation factor 5A–1 (eIF5A-1) gene in HPV 16 E6 induces cell growth in human cervical squamous carcinoma cells. Biochem Biophys Res Commun. 2018;504:6–12.
Plaks V, Kong N, Werb Z. The cancer stem cell niche: how essential is the niche in regulating stemness of tumor cells? Cell Stem Cell. 2015;16:225–38.
Fendler A, Bauer D, Busch J, et al. Inhibiting WNT and NOTCH in renal cancer stem cells and the implications for human patients. Nat Commun. 2020;11:929.
We thank the supports of Prof. Jianjun Liu (Renji Hospital, Shanghai Jiaotong University School of Medicine).
We thank the National Natural Science Foundation of China (Grant No. 81830052, 82003142 and 81530053) and the Shanghai Key Laboratory of Molecular Imaging (No. 18DZ2260400) for supporting this research.
Ethical approval and consent to participate
The study has been examined and certified by the Ethics Committee of Shanghai University of Medicine and Health Sciences, and informed consent was obtained from all participants included in the study, in agreement with institutional guidelines.
Consent for publication
All authors have agreed to publish this manuscript.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
: High-throughput sequencing detection show the different expression of circRNA from three NSCLC and three adjacent normal tissue samples.
: Bioinformatics analysis prediction the relationship among miR-587, miR-486-3p and miR-595 to circ-EPB41.
: High-throughput sequencing for miRNA expression analysis in A549 cells between NC and sh-circEPB41 group.
: High-throughput sequencing for mRNA expression analysis in A549 cells between NC and sh-circEPB41 group.
About this article
Cite this article
Jin, M., Liu, X., Wu, Y. et al. Circular RNA EPB41 expression predicts unfavorable prognoses in NSCLC by regulating miR-486-3p/eIF5A axis-mediated stemness. Cancer Cell Int 22, 219 (2022). https://doi.org/10.1186/s12935-022-02618-7