Circular RNA circEMB promotes osteosarcoma progression and metastasis by sponging miR-3184-5p and regulating EGFR expression
Biomarker Research volume 11, Article number: 3 (2023)
Osteosarcoma (OSA) is the most prevalent type of bone cancer with a high rate of metastasis. Circular RNAs (CircRNAs) play an essential role in multiple aspects of tumour biology. This study aimed to elucidate the role of circEMB in OSA.
circRNAs related to OSA invasion were identified via RNA sequencing and qRT-PCR. The relationship between circEMB levels and clinicopathological features of OSA was examined using the clinical specimens and data of 53 patients with OSA. Several in vivo and in vitro experiments, including intravital imaging, whole-transcriptome sequencing, transwell assay, flow cytometry, dual-luciferase reporter assay, RIP assay, RNA pull-down assay and RNA-FISH, were performed to examine the effects of circEMB on the malignant behaviour of OSA.
A novel circRNA, named circEMB (hsa_circ_001310), was identified in this study. circEMB can promote the malignant behaviour of OSA. In vitro experiments revealed that circEMB knockdown decreased cell proliferation, inhibited tumour invasion and metastasis; increased apoptosis and resulted in G1/S phase arrest. In vivo experiments revealed that circEMB knockdown inhibited tumour growth and metastasis in xenograft-bearing mice. Mechanistically, circEMB affects the malignant behaviour of OSA by mediating EGFR as an miR-3184-5p sponge. In addition, the circEMB/miR-3184-5p/EGFR axis modulates methotrexate (MTX) resistance in OSA.
CircEMB plays a critical role in promoting cancer via the miR-3184-5p/EGFR pathway, indicating that circEMB may serve as a therapeutic target for OSA.
Osteosarcoma (OSA) is one of the most prevalent types of bone cancer, which often originates in the metaphysis of the extremities such as the distal femur and proximal tibia [1, 2]. Owing to its high degree of invasiveness, approximately 10–15% of newly diagnosed patients with OSA have metastatic disease, which reduces the 5-year survival rate from 60 to 20% and significantly affects the clinical prognosis of patients [3,4,5]. At present, therapeutic options for OSA include surgery and neoadjuvant chemotherapy. However, the survival rate of patients with metastasis and resistance to chemotherapy is relatively low . Numerous studies have reported that gene variants play a crucial role in the pathogenesis and prognosis of OSA [7, 8]. Therefore, molecular mechanisms underlying the pathogenesis and development of OSA should be elucidated, and specific therapeutic targets should be identified.
Circular RNAs (circRNAs) are a new type of endogenous RNA. They form a continuous covalent closed-loop structure via reverse splicing of precursor mRNA (pre-mRNA) transcripts and lack the 5′-cap structure and 3′-poly-A tail [9,10,11]. circRNAs are highly stable and conserved because they are resistant to exonuclease activity . Therefore, they can be used as novel biomarkers for diseases. circRNAs contribute to the occurrence and development of OSA by controlling the proliferation, invasion and metastatic abilities; apoptosis; metabolism and drug resistance of tumour cells [13,14,15]. Therefore, understanding the functions of circRNAs is important for elucidating molecular mechanisms underlying the development of OSA and for overcoming chemotherapy resistance.
In the present study, we identified a novel circRNA named circEMB (circBase ID: hsa_circ_001310 or hsa_circ_0001481) via high-throughput sequencing analysis of tissues specimens collected from patients with OSA . circEMB is formed via the splicing of exons 3–9 of the parent gene EMB (aliases: Embigin) [17, 18]. The results revealed that circEMB was upregulated in OSA tissues compared with normal tissues, and its increased expression was associated with the proliferation, invasion and metastatic abilities and apoptosis of OSA cells. Mechanistically, circEMB can affect the expression of epidermal growth factor receptor (EGFR) by regulating the sponging activity of miR-3184-5p, and the circEMB/miR-3184-5p/EGFR axis can regulate methotrexate (MTX) resistance.
Materials and methods
The human normal osteoblast cell line hFOB1.19 and human OSA cell lines 143B, U2OS, MNNG, MG63 and Saos-2 were obtained from the American Type Culture Collection (ATCC, Manassas, VA). The cells were cultured in Dulbecco’s Modified Eagle’s medium (DMEM) (Gibco, USA) supplemented with 10% fetal bovine serum (BI, Israel), 1% penicillin G and streptomycin (Gibco, USA). All OSA cell lines were cultured in an incubator with 5% CO2 at 37 °C, whereas hFOB1.19 cells were cultured at 34 °C and 5% CO2. The MTX-resistant OSA cell line U2OS was provided by Dr. M. Serra (Istituti Ortopedici Rizzoli, Bologna, Italy). U2OS cells were continuously cultured in a medium containing 300-ng/mL MTX.
RNA sequencing was performed and libraries were constructed by Gene Denovo Biotechnology (Guangzhou, China). The samples were processed and sequenced on an Illumina HiSeqTM 2500 platform. The offline data obtained were filtered and compared with the reference genome to obtain the gene expression data, which were compared in the RPM format. Ensembl_release98 served as the reference genome . RNA sequencing of circRNAs was performed after the removal of conventional ribosomal RNA and degradation of linear RNA with the RNase R enzyme. The obtained circRNAs were fragmented in a fragmentation buffer, and a library was subsequently constructed.
Real-time reverse transcription polymerase chain reaction (qRT-PCR)
Total RNA was extracted from whole cell lysates of paracancerous tissues, cancerous tissues and cells from patients with OSA using the RNAiso Plus kit (TaKaRa, Japan) according to the manufacturer’s instructions. The quality of the extracted RNA was determined on the NanoDrop 2000 spectrophotometer. Two transcription reagents were used to convert mRNA, circRNA and miRNA to cDNA. Except for miRNA, for which we used Evo M-MLV RT Kit (Accurate Biotechnology, China), we use PrimeScript RT Reagent Kit (TaKaRa, Japan). We used two reagents to perform real-time amplification for circRNA and miRNA, namely SYBR Green Premix Pro Taq HS qPCR Kit (Accurate Biotechnology, China) and SYBR Premix Ex Taq II (TaKaRa, Japan). All qRT-PCR reactions were performed utilizing the Bio-Rad CFX connect system (Bio-Rad, CA, USA). The specific program was designed as previously described . The relative gene expression was determined using the 2 –ΔΔCT method. The specific primers are listed in the Additional file 1: Table S1.
Fluorescence in situ hybridization (FISH)
FISH assay was conducted for detection of the location of circEMB and miR-3184-5p in OSA cells. Cy3-labeled circEMB and FAM-labeled miR-3184-5p probes were used to detect the localization of circEMB and miR-3184-5p in OSA cells, respectively. According to the manufacturer’s instructions, the cells were inoculated in confocal dishes and fixed using 4% paraformaldehyde. After that, the cell was incubated with oligo probe overnight at 37 °C. Finally, samples were acquired by confocal microscopy (AX confocal, Nikon, Tokyo, Japan) to obtain images.
RNA interference (RNAi)
For circEMB knockdown, two shRNAs targeting the back-splice junction of circEMB were generated by IGE Co., Ltd. (Guangzhou, China). Then, shRNA against circEMB and shRNA-NC as negative control were packaged into lentiviruses by Genechem Co., Ltd. (Shanghai, China). The miR-3184-5p mimic/inhibitor, and their corresponding negative controls (NC) were purchased from IGE Co., Ltd. (Guangzhou, China). The lentivirus-miR-3184-5p inhibitor and lentivirus-sh-EGFR were also procured from Genechem Co., Ltd. (Shanghai, China). All vectors were validated via sequencing. Lentivirus infection was performed using the HiTransG reagent (Genechem, China), as per the manufacturer’s specifications.
Cell proliferation assay
The proliferation of OSA cells was detected using the cell counting kit-8 (CCK-8, GLPBIO, USA). Cells in different groups were quantified using a cell counting method. Subsequently, 2000 OSA cells (100 uL) in the logarithmic phase were seeded in a 96-well plate, and phosphate-buffer saline (PBS) was introduced into the marginal well. The growth and proliferation of cells were detected after 6, 24, 48 and 72 h by CCK-8 assay. A total of 10 uL of CCK-8 detection solution was introduced into every well and subjected to incubation at 37 °C for 2 h.
Colony formation assay
A total of 800 OSA cells were seeded in a 6-well plate, and the medium was changed every 2 days. After 14 days of culture at 37 °C, cell colonies were stained with 4% polymethanol and 0.1% crystal violet and counted.
Migration and invasion assays
transwell chamber was used to assess the migration and invasion abilities of OSA cells. Transfected OSA cells were re-suspended in 200 uL of DMEM and added to the upper compartment, whereas a medium containing 10% FBS was added to the lower compartment. The transwell chamber was incubated at 37 °C for 24 h. Thereafter, cells in the upper compartment were wiped using cotton swabs, whereas those in the lower compartment were fixed with paraformaldehyde, stained with crystal violet and counted. For invasion assay, the upper compartment was precoated with 50 uL of Matrigel (BD Bioscience, USA) to detect intrusion.
Wound healing assays
OSA cells were seeded in 6-well plates, and a scratch was created in the cell monolayer with a 200 μL pipette tip after 90% confluency was achieved. After the scratch was created, the cells were cultured in a serum-free medium and photographed at 0 and 24 h.
Flow cytometry was used to assess apoptosis and cell cycle progression. Cells in the logarithmic growth phase were digested in trypsin without ethylenediaminetetraacetic acid (EDTA) and washed thrice with PBS. Apoptosis was determined using the Annexin V-APC Apoptosis Assay kit (BestBio, China) according to the manufacturer’s instructions. For cell cycle analysis, cells were lysed in trypsin containing EDTA, washed thrice with PBS and fixed with ethanol overnight. Cell cycle progression was determined according to the instructions on the Cell Cycle Analysis kit (BestBio, China). Each experiment was repeated at least thrice.
Dual-luciferase reporter assay
The 3′-untranslated region (UTR) sequences of circEMB and EGFR and the corresponding mutants were cloned into pmiR-RB-REPORT™ vectors, which were named circEMB-WT, circEMB-MUT, EGFR-WT and EGFR-WUT. The luciferase plasmid was co-transfected with miR-3184-5p mimics or miR-NC. The relative luciferase activity was determined according to the instructions on the Dual-Luciferase Reporter Assay System (Promega, USA).
RNA immunoprecipitation (RIP) was performed using the Magna RIP Kit (Millipore, Billerica, MA, USA) according to the manufacturer’s instructions. Briefly, OSA cells were lysed in RNA immunoprecipitation (RIP) buffer supplemented with protease inhibitors and RNase. Magnetic beads were pre-incubated with anti-IgG and anti-Ago2 antibodies and were subsequently incubated with the aforementioned cell lysates overnight at 4 °C. Finally, qRT-PCR and western blotting were performed to determine the presence of binding targets.
RNA pull-down assay using a biotin-labelled probe
C1 magnetic beads were incubated with a biotinylated circEMB probe or oligonucleotide probe (IGE Biotech Co., Ltd., Guangzhou, China) at 25 °C for 2 h. The solution was co-incubated with a cell lysis solution at 4 °C overnight. The precipitates were extracted and purified using an RNeasy Mini Kit (Qiagen, USA) and assessed via qRT-PCR according to the manufacturer’s instructions.
The PARIS kit (Invitrogen, USA) was used to isolate nuclear and cytosolic fractions from cells according to the manufacturer’s instructions.
Multiple databases, including ENCORI (http://starbase.sysu.edu.cn/), TargetScan (http://www.targetscan.org/vert_72/), miRanda (http://www.microrna.org/), RNAhybrid (http://bibiserv.techfak.uni-bielefeld.de/rnahybrid/) and RNA22 (https://cm.jefferson.edu/rna22/), were used to predict the interaction between circEMB and miRNA [21,22,23,24]. The interaction between miR-3184-5p and EGFR was predicted using TargetScan, miRanda, RNAhybrid and miRDB (http://mirdb.org/) . Pathway enrichment analysis was performed using Reactome in the KOBAS software (KOBAS, Surrey, UK).
Western blot (WB) analysis
Total protein was isolated from tissues and cells lysed in radioimmunoprecipitation assay (RIPA) buffer (Fudebio, China) and quantified via bicinchoninic acid (BCA) assay. Subsequently, SDS-PAGE (Fudebio, China) was performed to separate proteins of varying molecular weights. The separated proteins on the gel were transferred to a polyvinylidene difluoride (PVDF) membrane (Invitrogen). The membrane was blocked with 5% skim milk at room temperature for 1.5 h and incubated with primary antibodies overnight at 4 °C. The following day, the membrane was incubated with horseradish peroxidase (HRP)-conjugated secondary antibodies (1:10000) (Proteintech, China) at room temperature for 1 h, and protein bands were visualised via enhanced chemiluminescence (Millipore, USA).
Haematoxylin and eosin (HE) and immunohistochemical (IHC) staining
The extracted fresh animal tissues were washed with normal saline to remove excess blood and immediately fixed in 4% paraformaldehyde overnight. The tissues were embedded in paraffin using an automatic dehydrator and a paraffin-embedding machine and stained with haematoxylin and eosin (H&E) to visualise the histopathological features. Immunohistochemical (IHC) analysis was performed as described previously .
Construction of a xenograft-bearing nude mouse model
All animal experiments were approved by the Ethics Committee of Zhujiang Hospital, Southern Medical University. Nude mice were purchased from the Experimental Animal Center of Southern Medical University. The effects of circEMB on the proliferation of OSA cells were examined using tumour xenograft-bearing mouse models. Stable lentivirus-infected cell lines were isolated, resuspended in 100 μL of PBS (143B, 1 × 106) and injected into the right hindlimb of female nude mice (age, 4 weeks). Malignant growth was examined using an in vivo imaging system. After 4 weeks, mice were euthanised via carbon dioxide inhalation, and tumour tissues were harvested.
To determine the effects of circEMB and its signalling axis on the metastatic capability of OSA cells, 2 × 106 cells (200 uL; 143B) was injected into the tail vein of male nude mice (age, 4 weeks). After 4 weeks, tumour growth in the lung was detected using an in vivo imaging system, and mice were sacrificed via cervical dislocation. Thereafter, lung tissues were harvested for HE staining.
AutodockTools (version 1.5.6) was used for molecular docking. The 2D structure of MTX (the ligand) was downloaded from PubMed and converted into a 3D structure using Chem3D. The crystal structure of EGFR was downloaded from Protein Data Bank (PDB) (PDB code: 5UG9) and was used as the protein target for docking simulation. Docking results with the highest score were visualised and analysed using PyMOL.
The SPSS Statistics (version 22.0) software (Chicago, Illinois, USA) was used for statistical analysis. The Mann–Whitney U test or two-tailed Student’s t-test was used to determine the significance of between-group differences, whereas ANOVA followed by Tukey’s test was used to determine the significance of differences among multiple groups. Survival curves were approximated using the Kaplan–Meier approach, and survival data were compared using the log-rank test. Pearson correlation analysis was performed to determine the link between gene expression and clinicopathological features. A P-value of < 0.05 was considered significant. Data are expressed as the mean ± standard deviation (SD).
Characteristics of circEMB in OSA tissues
RNA-seq was used to determine circRNA expression in paired OSA tissues and adjacent normal tissues from 3 patients with OSA (Fig. 1a). A total of 404 differentially expressed circRNAs (LogFC > 1, P < 0.05) were identified between tumour and normal tissues; of which, 320 were included in the circBase database. Of these 320 circRNAs, 157 were upregulated and 163 were downregulated in tumour tissues. The migration and metastasis of OSA is greatly important indicators for better prognosis, which deserves further studies. To identify circRNAs associated with metastasis, qRT-PCR was performed to determine the expression patterns of the top 10 upregulated circRNAs in the metastatic and non-metastatic tissues of patients (Fig. 1b). The expression of hsa_circ_001310 was found to be significantly different between metastatic and non-metastatic tissues. Based on the results of qRT-PCR and sequencing, we hypothesised that hsa_circ_001310 is an important driver of carcinogenesis. hsa_circ_001310 (referred to as circEMB throughout the article) contains 954 nucleotides and is derived from the EMB gene. Sanger sequencing revealed that circEMB is formed via head-to-tail splicing of exons 3–9 (chromosome 5: 49694940–49,707,217) (Fig. 1c).
The expression of circEMB was found to be high in OSA cells (143B, U2OS, MNNG, MG63 and Saos-2) and hFOB1.19 cells, with the expression being highest in 143B and U2OS cells (Fig. 1d).To investigate the characteristics of the newly discovered circEMB, RNA extract was treated with RNase R, and the stability of circRNA was compared with linear RNA. The results demonstrated that the linear form of EMB was easily degraded by RNase R, whereas the circular form could resist RNase R (Fig. 1e). To rule out both trans-clipping and genome rearrangement, PCR was performed for cDNA and genomic DNA (gDNA) using divergent and convergent primers, respectively, and the findings. Conventionally, the divergent primers could amplify circEMB but not gDNA (Fig. 1f). FISH and cytoplasmic separation experiments revealed that circEMB was predominantly cytoplasmic (Fig. 1g-i). Altogether, these results indicate that circEMB was stable and upregulated in OSA cells.
circEMB expression in OSA tissues and its relationship with clinical characteristics
The relationship between circEMB expression (53 OSA tissues and 53 matched adjacent normal tissues) and the clinicopathological characteristics of patients was examined. Patients with OSA exhibited significantly upregulated circEMB expression (P = 0.00608) (Fig. 2a) (Table 1). Additionally, circEMB expression was correlated with metastasis and clinical stages of OSA (Fig. 2b-c), indicating that the degree of malignancy is correlated with circEMB expression. Receiver operating characteristic (ROC) curve analysis also revealed that circEMB could serve as a biomarker of osteosarcoma metastasis, with an AUC of 0.703 (Fig. 2e). According to the results of qRT-PCR, the samples were categorised into two groups based on the median circEMB expression as follows: low-circEMB-expression (n = 26) and high-circEMB-expression (n = 27) (Fig. 2f). Kaplan–Meier analysis revealed that survival was worse in the high-circEMB-expression group (Fig. 2d). Altogether, these results suggest that high circEMB expression is associated with the poor prognosis of OSA.
circEMB promotes the proliferation, migration and invasion capabilities of OSA cells in vitro
Several in vitro experiments were performed to examine the regulatory effects of circEMB on OSA cells. Two shRNAs targeting circEMB were constructed, and their knockout efficiency was verified in the two cell lines (143B and U2OS) with the highest expression of circEMB, without affecting the mRNA expression of EMB (Fig. 3a). The results showed that shRNA-1 not only had the highest knockout efficiency in 143B and U2OS cells but also had a significantly greater effect on OSA cells than shRNA-2. CCK-8 (Fig. 3b) and clone formation (Fig. 3c) assays revealed that the proliferation of OSA cells was inhibited after circEMB knockdown. Transwell and wound healing assays revealed that circEMB knockdown remarkably suppressed the migration and invasion abilities of OSA cells (Fig. 3d-e). Additionally, flow cytometry showed that circEMB knockdown induced apoptosis in OSA cells (Fig. 3f) and resulted in G1/S-phase arrest (Fig. 3g). shRNA-1 was selected for further analysis because its inhibitory effects on the malignant behaviour of OSA were more prominent. Altogether, these results suggest that circEMB promotes the malignant behaviour of OSA.
circEMB serves as an miR-3184-5p sponge
Given that circEMB localises in the cytoplasm, it may function as a molecular sponge of miRNAs after transcription. Therefore, RIP assay was performed using Ago2 to examine the role of circEMB in sponging miRNAs. The results revealed that circEMB expression was remarkably increased in the Ago2 group, indicating that circEMB binds to miRNAs via the Ago2 protein (Fig. 4a-b).
Two online databases, namely, RNA22 and ENCORI, and three local databases, namely, TargetScan, miRanda and RNAhybrid, were used to predict the miRNA targets of circEMB, and seven miRNAs (hsa-miR-135a-5p, hsa-miR-135b-5p, hsa-miR-2335-5p, hsa-miR-3184-5p, hsa-miR-370-5p, hsa-miR-423-5p, hsa-miR-522-5p) were eventually identified (Fig. 4c). However, the binding mechanism of the same circRNA to the presence of different miRNAs in different tumors affects cell growth [26,27,28]. Therefore, we need to explore the specific miRNA that circEMB binds in OSA. Next, we used biotinylated circEMB probes to determine the mode of interaction between circEMB and the seven predicted miRNAs to identify more stable interactions between them in OSA. Finally, we found that circEMB and miR-3184-5p may have a stronger and more stable relationship (Fig. 4d-e).
Given that miR-3184-5p plays a substantial role in tumorigenesis in OSA, we examined whether it can act as an Ago2 sponge. To ensure the validity of subsequent miR-3184-5p-related experiments, we validating the transfection efficiencies of the miR-3184-5p mimic and inhibitor (Additional file 2: Fig. S1A-B). RIP assays with anti-Ago2 and anti-IgG revealed that miR-3184-5p overexpression increased the adsorption of circEMB on anti-Ago2 instead of anti-IgG (Fig. 4f-h). The dual-luciferase reporter assay further verified the binding between circEMB and miR-3184-5p (Fig. 4i-j). Meanwhile, RNA-FISH experiment results showed that circEMB and miR-3184-5p were co-localized in the cytoplasm of OSA (Fig. 4k).
Altogether, these results indicate that circEMB forms a complex with Ago2 and affects Ago2-mediated miR-3184-5p-guided RNAi and Ago2-dependent miRNA biogenesis.
Inhibition of miR-3184-5p mediates the oncogenic function of circEMB
Owing to the lack of studies on miR-3184-5p in OSA, we validated its expression profile in clinical OSA samples (n = 53) via qPCR and examined whether miR-3184-5p expression was associated with the clinicopathological features of patients with OSA. The results revealed that the loss of miR-3184-5p was associated with advanced pathological stages, TNM stage and invasion depth (P < 0.05) (Table 2) (Additional file 3: Fig. S2A-B). We categorized the enrolled patients into the high-expression group (n = 27) and low-expression group (n = 26) according to the median value of miR-3184-5p expression in OSA tissues (Additional file 3: Fig. S2C). Kaplan-Meier survival analysis (log-rank test) exhibited that high expression levels of miR-3184-5p indicate good prognostic outcomes (Additional file 3: Fig. S2D). Additionally, a negative correlation was observed between the expression of miR-3184-5p and circEMB (Additional file 3: Fig. S2E). Consistently, miR-3184-5p expression was lower in OSA cell lines (143B, U2OS, MNNG, MG63 and Saos-2) than in the normal human osteoblast cell line hFOB1.19 (Additional file 3: Fig. S2F). These results indicate that miR-3184-5p plays an important role as a tumour suppressor.
To explore the functional role of miR-3184-5p, we transfected miR-3184-5p mimics into 143B and U2OS cells. All the results displayed that miR-3184-5p could significantly affect the ability of proliferation, migration, invasion, apoptosis and cell cycle arrest of OSA cells (Additional file 4: Fig. S3A-K). These results further demonstrate that high expression of miR-3184-5p suppresses OSA progression and metastasis.
Based on the abovementioned results, we hypothesised that circEMB affects the occurrence of OSA via miR-3184-5p. To test this hypothesis, we treated circEMB-knockdown cells with an inhibitor of miR-3184-5p and subsequently performed CCK-8 and colony formation assays (Fig. 5a-c). The miR-3184-5p inhibition was found to reverse the circEMB knockdown-induced decrease in the proliferation of 143B and U2OS cells. Transwell and wound healing assays revealed that miR-3184-5p inhibition attenuated the circEMB knockdown-induced decrease in the migration and invasion abilities of 143B and U2OS cells (Fig. 5d-e). Compared with circEMB knockdown alone, combined circEMB knockdown and miR-3184-5p inhibition significantly inhibited apoptosis and cell cycle arrest (Fig. 5f-g). Altogether, these results suggest that circEMB promoted tumour progression in OSA primarily by inhibiting the antitumour effects of miR-3184-5p.
circEMB regulates EGFR expression via miR-3184-5p
circRNAs sponge miRNAs to prevent them from exerting regulatory effects on target genes. RNA sequencing was performed after silencing circEMB to identify the downstream target genes regulated by circEMB (Fig. 6a). Because circRNAs function as miRNA sponges to attenuate the miRNA-induced suppression of downstream target genes, like cricRNAs, have co-directional expression. Additionally, TargetScan, miRanda, RNAhybrid and miRDB were used to identify the target genes of miR-3184-5p. The results revealed 394 genes as the direct targets of miR-3184-5p (Fig. 6b); of which, 50 genes were selected for subsequent experiments. These 50 genes were used for pathway enrichment analysis using KOBAS to identify significant signalling pathways involved in the regulation of the circEMB/miR-3184-5p axis (Fig. 6c). Pathway enrichment analysis and GSEA revealed that the target genes were primarily associated with the Notch signalling pathway (Fig. 6d). Thereafter, four genes involved in the Notch signalling pathway, namely, EGFR, ARRB2, PSME3 and ATP2A3, were identified as potential regulatory targets of the circEMB/miR-3184-5p axis. The expression patterns of the target genes in circEMB-knockdown OSA cells (143B and U2OS) revealed that EGFR expression was more closely associated with circEMB expression (Fig. 6e-f). Because miRNAs always target the 3′-UTR of mRNAs, the sequences of miR-3184-5p and the 3′-UTR of EGFR were compared. Dual-luciferase reporter assay revealed a specific binding site for miR-3184-5p in the 3′-UTR of EGFR (Fig. 6g-h). Pearson analysis revealed that EGFR expression was positively correlated with circEMB expression and negatively correlated with miR-3184-5p expression in OSA tissues (Fig. 6i-j). Additionally, EGFR expression was significantly higher in OSA tissues than in paracancerous tissues (Fig. 6k). Altogether, these results suggest that EGFR is a direct target of miR-3184-5p.
circEMB promotes OSA progression and chemoresistance via the circEMB/miR-3184-5p/EGFR axis
We hypothesised that circEMB can antagonise miR-3184-5p. To verify this hypothesis, we examined whether circEMB affects EGFR expression via miR-3184-5p, thus regulating the occurrence and progression of OSA. Rescue experiments were performed via circEMB knockdown, EGFR knockdown and miR-3184-5p inhibition to examine whether circEMB acts as an oncogene by modulating EGFR. CCK-8 assay, transwell migration and invasion assay and flow cytometry revealed that circEMB downregulation reversed the effects of miR-3184-5p inhibition on the proliferation, migration and invasion abilities; apoptosis and cell cycle arrest of OSA cells, whereas EGFR knockdown reversed the circEMB and miR-3184-5p knockdown-induced cell functional phenotype (Fig. 7a-f). Additionally, the results of western blotting validated that the circEMB/miR-3184-5p/EGFR axis regulated the expression of cell proliferation-related proteins (PCNA), apoptosis-related proteins (BAX and BCL2), cell cycle-related proteins (CDK4 and CCND1) and epithelial–mesenchymal transformation-related proteins (Vimentin, E-cadherin and N-cadherin) by affecting EGFR expression. Western blotting also revealed that the circEMB/miR-3184-5p/EGFR axis has the potential to influence the Notch signalling pathway (EGFR and Jagged1) (Fig. 7g-h). In vivo experiments showed that the circEMB/miR-3184-5p/EGFR axis affected OSA of capable for subcutaneous tumour formation (Fig. 7i). IHC analysis revealed that circEMB knockdown decreased the protein expression of Ki-67, PCNA, N-cadherin and Vimentin but increased the protein expression of E-cadherin. However, silencing the miR-3184-5p/EGFR axis partially reversed these effects (Fig. 7j). In vivo imaging of the tail vein (in mouse models of lung metastasis) showed that the circEMB/miR-3184-5p/EGFR axis affected OSA of capable for metastasis (Fig. 7k). HE staining further validated these results (Fig. 7l). Collectively, Altogether, these results suggest that circEMB promotes the malignant behaviour of OSA by regulating the miR-3184-5p/EGFR axis.
The Notch pathway is tightly associated with the chemosensitivity of tumour cells. Because the circEMB/miR-3184-5p/EGFR axis affects the development of OSA via the Notch signalling pathway, we examined the effects of circEMB on the chemoresistance of OSA and verified the results using MTX.
qRT-PCR revealed that circEMB expression was higher in the MTX-resistant OSA cell line U2OS (U2OS/MTX cells) (Fig. 8a). A study from Austria showed that gefitinib, an EGFR inhibitor, in combination with MTX can increase the sensitivity of OSA cells to MTX . Therefore, we examined the relationship between EGFR and MTX. Molecular docking showed that EGFR and MTX had robust interactions and strong binding affinity (− 7.9 kcal/mol). The amino acid residues ALA871, LYS913, LEU718 and GLY796 on the binding domain of EGFR were found to bind to MTX via hydrophobic interactions. The amino acid residues ALA743, LEU844, VAL726, LEU792, GLN791, PHE856, ASN842, LYS875, ARG858, PRO877, GLY721, GLY719 and GLY796LEU718 of EGFR were found to bind to MTX via van der Waals forces. Additionally, MTX formed π-π interactions with the amino acid residue PHE723 of EGFR (Fig. 8b-g). To explore the potential role of circEMB/miR-3184-5p/EGFR axis on MTX resistance in OSA, we perturbed the expression of circEMB/miR-3184-5p/EGFR axis in U2OS/MTX. We found that the IC50 of U2OS/MTX to MTX decreased after circEMB expression was reduced, and the resistance phenotype of U2OS/MTX was reversed after miR-3184-5p knockdown. Similarly, we found that EGFR as part of the biological axis also had an effect on the IC50 of MTX (Fig. 8h). Flow cytometry revealed that the circEMB/miR-3184-5p/EGFR axis regulated drug sensitivity via apoptosis and cell cycle regulation (Fig. 8i-j).
At present, the management strategies for newly diagnosed OSA include individualised high-dose neo-adjuvant chemotherapy (also called preoperative chemotherapy), surgery and adjuvant chemotherapy after surgery. However, owing to metastasis and drug resistance, the 5-year survival rate of patients with OSA is as low as 20% . Therefore, it is of clinical significance to elucidate cellular mechanisms underlying the development of OSA and develop new treatment strategies for improving the prognosis.
Epigenetic alterations, including DNA and m6A methylation, and noncoding RNAs involved in tumour development, invasion, metastasis and drug resistance play an important role in guiding clinical treatment [31,32,33,34,35]. circRNAs, noncoding RNAs, are potential biomarkers for evaluating tumour metastasis and drug resistance outcomes in several cancers, including hepatocellular carcinoma, breast cancer and glioma [36,37,38,39]. However, the mechanisms through which circRNAs affect the development, metastasis and drug resistance of OSA remain unclear.
This study revealed that the circEMB/miR-3184-5p/EGFR axis regulates the progression of OSA by regulating cell proliferation, tumour invasion and metastasis, apoptosis and cell cycle arrest. Additionally, the circEMB/miR-3184-5p/EGFR axis is involved in resistance to MTX (Fig. 9).
A fundamentally novel circRNA, named circEMB, was identified via high-throughput sequencing of tissue samples collected from patients with OSA. circEMB is upregulated in OSA tissues and promotes the proliferation, migration and invasion abilities of OSA cells. To the best of our knowledge, this study is the first to report that circEMB promotes the malignant behaviour of OSA.
circEMB localises in the cytoplasm and is significantly associated with miRNAs. Downstream miRNAs are important agents through which circRNAs affect the target genes. However, most studies investigating the relationship between circRNAs and miRNAs lack systematic experimental verification and mainly focus on identifying candidate miRNAs via bioinformatic analysis . In this study, miR-3184-5p was identified as the target miRNA of circEMB. Importantly, circEMB was verified to bind the 3’UTR region of miR-3184-5p through Ago2.
miR-3184-5p inhibits the malignant behaviour of various tumours [41,42,43,44]. However, previous studies have not reported the role of miR-3184-5p in tumorigenesis in OSA. To the best of our knowledge, this is the first study to report that miR-3184-5p is found in OSA tissues and cells and is negatively associated with circEMB expression. circEMB can sponge miR-3184-5p to rescue the functions of downstream genes that affect the malignant behaviour of OSA.
Furthermore, EGFR was identified as a potential functional target of miR-3184-5p. EGFR, an epidermal growth factor family member, is closely associated with cell proliferation, apoptosis and tumour invasion and metastasis in OSA [45,46,47,48], which is consistent with the findings of this study. Increased EGFR expression was positively correlated with circEMB expression and negatively correlated with miR-3184-5p expression in OSA tissues and cells. Functionally, after successive circEMB and miR-3184-5p knockdown, the promotion of cell proliferation and invasion triggered by EGFR knockdown was remarkably reversed. These results suggest that circEMB affects EGFR expressions by sponging miR-3184-5p, thereby accelerating OSA progression.
Clinically, MTX is a front-line chemotherapeutic agent for patients with OSA . However, many patients with OSA exhibit primary and/or acquired resistance to MTX. In this study, circEMB was found to be highly expressed in MTX-resistant cells. Inhibiting the function of EGFR can alter the resistance of OSA to MTX . In this study, molecular docking revealed an interaction between EGFR and MTX. Based on the findings of this study and those of previous studies, we speculated that suppressed EGFR expression inhibits EGFR and alters the sensitivity of tumour cells to MTX. Experimentally, we found that the circEMB/miR-3184-5p/EGFR axis can affect drug resistance in OSA.
This retrospective study has several limitations, and there is scope for future work. First, the mouse model of OSA was not used in situ; therefore, we could not completely simulate a real tumour microenvironment of OSA in vivo. Second, the molecular mechanisms through which the circEMB/miR-3184-5p/EGFR axis regulates MTX resistance warrant further investigation.
CircEMB, a previously unrecognised oncogenic driver, is upregulated in OSA tissues and is associated with cell proliferation, migration and invasion in OSA. In addition, we propose the mechanism by which circEMB can block the EGFR inhibitory activity of miR-3184-5p. The circEMB/miR-3184-5p/EGFR axis affects the sensitivity of OSA cells to MTX. Therefore, circEMB is a candidate prognostic biomarker and therapeutic target for OSA.
Availability of data and materials
The datasets used and/or analysed in the present study are available from the corresponding author upon reasonable request.
Short hairpin RNAs
Fluorescence in situ hybridization
Haematoxylin and eosin
Epidermal growth factor receptor
Proliferating cell nuclear antigen
B-cell leukaemia/lymphoma 2-associated X
B-cell leukaemia/lymphoma 2
Cyclin-dependent kinase 4
Jagged canonical Notch ligand 1
Meazza C, Scanagatta P. Metastatic osteosarcoma: a challenging multidisciplinary treatment. Expert Rev Anticanc. 2016;16:543–56.
Mirabello L, Troisi RJ, Savage SA. Osteosarcoma incidence and survival rates from 1973 to 2004. Cancer-Am Cancer Soc. 2009;115:1531–43.
Ritter J, Bielack SS. Osteosarcoma. Ann Oncol. 2010;21:i320–5.
Bielack SS, Kempf-Bielack B, Delling G, Exner GU, Flege S, Helmke K, et al. Prognostic factors in high-grade osteosarcoma of the extremities or trunk: an analysis of 1,702 patients treated on neoadjuvant cooperative osteosarcoma study group protocols. J Clin Oncol. 2002;20:776–90.
Belayneh R, Fourman MS. Bhogal S. Weiss KR. Update on Osteosarcoma. Curr Oncol Rep. 2021;23.
Isakoff MS, Bielack SS, Meltzer P, Gorlick R. Osteosarcoma: current treatment and a collaborative pathway to success. J Clin Oncol. 2015;33:3029–35.
Sayles LC, Breese MR, Koehne AL, Leung SG, Lee AG, Liu H, et al. Genome-informed targeted therapy for osteosarcoma. Cancer Discov. 2018;9:17–1152.
Kansara M, Teng MW, Smyth MJ, Thomas DM. Translational biology of osteosarcoma. Nat Rev Cancer. 2014;14:722–35.
Memczak S, Jens M, Elefsinioti A, Torti F, Krueger J, Rybak A, et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013;495:333–8.
Lasda E, Parker R. Circular RNAs: diversity of form and function. RNA. 2014;20:1829–42.
Jeck WR, Sharpless NE. Detecting and characterizing circular RNAs. Nat Biotechnol. 2014;32:453–61.
Barrett SP, Salzman J. Circular RNAs: analysis, expression and potential functions. Development. 2016;143:1838–47.
Li Z, Li X, Xu D, Chen X, Li S, Zhang L, et al. An update on the roles of circular RNAs in osteosarcoma. Cell Proliferat. 2021;54:e12936.
Li S, Liu F, Zheng K, Wang W, Qiu E, Pei Y, et al. CircDOCK1 promotes the tumorigenesis and cisplatin resistance of osteogenic sarcoma via the miR-339-3p/IGF1R axis. Mol Cancer. 2021;20:161.
Zhao Z, Shen J. Circular RNA participates in the carcinogenesis and the malignant behavior of cancer. RNA Biol. 2017;14:514–21.
Glažar P, Papavasileiou P, Rajewsky N. circBase: a database for circular RNAs. RNA. 2014;20:1666–70.
Jung DE, Kim JM, Kim C, Song SY. Embigin is overexpressed in pancreatic ductal adenocarcinoma and regulates cell motility through epithelial to mesenchymal transition via the TGF-β pathway. Mol Carcinogen. 2016;55:633–45.
Chao F, Zhang J, Zhang Y, Liu H, Yang C, Wang J, et al. Embigin, regulated by HOXC8, plays a suppressive role in breast tumorigenesis. Oncotarget. 2015;6:23496–509.
Yates AD, Achuthan P, Akanni W, Allen J, Allen J, Alvarez-Jarreta J, et al. Ensembl 2020. Nucleic Acids Res. 2020;48:682–8.
Tan J, Liang H, Yang B, Zhu S, Wu G, Li L, et al. Identification and analysis of three hub prognostic genes related to osteosarcoma metastasis. J Oncol. 2021;2021:1–16.
Rehmsmeier M, Steffen P, Hochsmann M, Giegerich R. GIEGERICH R: fast and effective prediction of microRNA/target duplexes. RNA. 2004;10:1507–17.
Li J, Liu S, Zhou H, Qu L, Yang J. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein–RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 2013;42:D92–7.
McGeary SE, Lin KS, Shi CY, Pham TM, Bisaria N. Kelley GM. Bartel DP. The biochemical basis of microRNA targeting efficacy. Science. 2019;366.
Desvignes T, Loher P, Eilbeck K, Ma J, Urgese G, Fromm B, et al. Unification of miRNA and isomiR research: the mirGFF3 format and the mirtop API. Bioinformatics. 2020;36:698–703.
Chen Y, Wang X. miRDB: an online database for prediction of functional microRNA targets. Nucleic Acids Res. 2020;48:D127–31.
Chen W, Wu M, Cui S, Zheng Y, Liu Z, Luo L. CircRNA Circ-ITCH inhibits the proliferation and invasion of glioma cells through targeting the miR-106a-5p/SASH1 Axis. Cell Transplant. 2021;30:1504089270.
Yang C, Yuan W, Yang X, Li P, Wang J, Han J, et al. Circular RNA circ-ITCH inhibits bladder cancer progression by sponging miR-17/miR-224 and regulating p21, PTEN expression. Mol Cancer. 2018;17:19.
Guo Y, Zheng H, Yin J, Wang H. Rs4911154 of circ-ITCH aggravated tumor malignancy of thyroid nodules via the circ-ITCH/miR-22-3p/CBL axis. Sci Rep-Uk. 2021;11:18491.
Sevelda F, Mayr L, Kubista B, Lötsch D, van Schoonhoven S, Windhager R, et al. EGFR is not a major driver for osteosarcoma cell growth in vitro but contributes to starvation and chemotherapy resistance. J Exp Clin Canc Res. 2015;20:34.
Meltzer PS, Helman LJ. New horizons in the treatment of osteosarcoma. New Engl J Med. 2021;385:2066–76.
Dawson MA, Kouzarides T. Cancer epigenetics: from mechanism to therapy. Cell. 2012;150:12–27.
Esteller M. Cancer epigenomics: DNA methylomes and histone-modification maps. Nat Rev Genet. 2007;8:286–98.
Kanwal R, Gupta S. Epigenetic modifications in cancer. Clin Genet. 2012;81:303–11.
Bhol CS, Panigrahi DP, Praharaj PP, Mahapatra KK, Patra S, Mishra SR, et al. Epigenetic modifications of autophagy in cancer and cancer therapeutics. Semin Cancer Biol. 2020;66:22–33.
Huang X, Huang L, Xie Q, Zhang L, Huang S, Hong M, Li J, Huang Z, Zhang H. LncRNAs serve as novel biomarkers for diagnosis and prognosis of childhood ALL. Biomark Res. 2021;9:1–11.
Sang Y, Chen B, Song X, Li Y, Liang Y, Han D, et al. circRNA_0025202 regulates tamoxifen sensitivity and tumor progression via regulating the miR-182-5p/FOXO3a Axis in breast Cancer. Mol Ther. 2019;27:1638–52.
Xu J, Ji L, Liang Y, Wan Z, Zheng W, Song X, et al. CircRNA-SORE mediates sorafenib resistance in hepatocellular carcinoma by stabilizing YBX1. Signal Transduction and Targeted Therapy. 2020;5:298.
Sun J, Li B, Shu C. Ma Q. Wang J: Functions and clinical significance of circular RNAs in glioma. Mol Cancer. 2020;19.
Huang S, Gong N, Li J, Hong M, Li L, Zhang L, Zhang H. The role of ncRNAs in neuroblastoma: mechanisms, biomarkers and therapeutic targets. Biomark Res. 2022;10:1–11.
Huang Y, Zhang Z, Wang J, Shen S, Yao T, Xu Y, et al. circSPG21 protects against intervertebral disc disease by targeting miR-1197/ATP1B3. Exp Mol Med. 2021;53:1547–58.
Liu H, Li Y, Zhang J, Wu N, Liu F, Wang L, et al. Erb-B2 receptor tyrosine kinase 2 is negatively regulated by the p53-responsive microRNA-3184-5p in cervical cancer cells. Oncol Rep. 2021;45:95–106.
Lin D, Fu Z, Yang G, Gao D, Wang T, Liu Z, et al. Exportin-5 SUMOylation promotes hepatocellular carcinoma progression. Exp Cell Res. 2020;395:112219.
Li D, Li Z, YanFei W, Wang Y, Shi J, Liu C, et al. LncRNA FOXP4-AS promotes the progression of non-small cell lung cancer by regulating the miR-3184-5p/EIF5A axis. J Tissue Eng Regen M. 2022;16:335–45.
Wu X, Xiao Y, Zhou Y, Zhou Z, Yan W. LncRNA FOXP4-AS1 is activated by PAX5 and promotes the growth of prostate cancer by sequestering miR-3184-5p to upregulate FOXP4. Cell Death Dis. 2019;10:472.
Feigin ME, Muthuswamy SK. ErbB receptors and cell polarity: new pathways and paradigms for understanding cell migration and invasion. Exp Cell Res. 2009;315:707–16.
Arteaga CL, Engelman JA. ERBB receptors: from oncogene discovery to basic science to mechanism-based cancer therapeutics. Cancer Cell. 2014;25:282–303.
Yuan Y, Wang Y, Liu Z, Sun Y, Yao Y, Yu W, et al. MAT2B promotes proliferation and inhibits apoptosis in osteosarcoma by targeting epidermal growth factor receptor and proliferating cell nuclear antigen. Int J Oncol. 2019;54:2019–29.
Lin H, Zhang C, Zhang H, Xia Y, Zhang C, Luo J, et al. Physakengose G induces apoptosis via EGFR/mTOR signaling and inhibits autophagic flux in human osteosarcoma cells. Phytomedicine. 2018;42:190–8.
Grem JL, King SA, Wittes RE. Leyland-Jones B: the role of methotrexate in osteosarcoma. J Natl Cancer Inst. 1988;80:626–55.
This study was funded by the Natural Science Foundation of Guangdong Province (Grant no.: 2022A1515010293) and the National Natural Science Foundation of China (Grant no.: 81801210).
Ethics approval and consent to participate
This study was approved by the Ethics Committee of Zhujiang Hospital, Southern Medical University (Guangzhou, China) (human OSA study: 2018-GJGBWK-002; animal experimental procedures: LAEC-2020–207). All experiments were performed in accordance with the Declaration of Helsinki. Patients with OSA who participated in this study provided informed consent.
Consent for publication
All authors have read and approved the final version of 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.
Table S1. Primer sequences for qRT-PCR in this study.
Fig. S1. Verification of transfection efficiency of miR-3184-5p in OSA cells.
(A) After the transfection of mimic and inhibitor of miR-3184-5p, qRT-PCR was used to detect the transfection efficiencies. Data are expressed as mean ± SD (n = 3) (*, p < 0.05; **, p < 0.01; ***, p < 0.001).
The relationship between the miR-3184-5p and clinicopathological characteristics. (A) Comparison of miR-3184-5p expression between non-metastatic and metastatic OSA tissues. (B) Distribution characteristics of miR-3184-5p expression in the TNM stage of OSA. (C) In the 53-patients cohort, according to the median value of miR-3184-5p expression, they were divided into high expression group and low expression group. (D) Overall survival analysis of patients with low and high expression of miR-3184-5p using log-rank test and Kaplan–Meier analysis. (E) Correlation analysis between circEMB and miR-3184-5p in OSA tissues. (F) Relative miR-3184-5p levels in OSA cell lines and normal osteoblasts were ascertained by qRT-PCR. Data are expressed as mean ± SD (n = 3) (*, p < 0.05; **, p < 0.01; ***, p < 0.001).
Overexpression of miR-3184-5p was associated with decreased OSA cell migration and invasion, nor increased apoptosis and G1/S phase arrest. (A-B) The proliferation ability of OSA cells after overexpression of miR-3184-5p was determined by CCK-8. (C-F) Transwell assays detected the changes of invasion and metastasis ability of OSA cells with or without miR-3184-5p overexpression. (C and E) Refer to 143B, while (D and F) refer to U2OS. Scale bar, 100 μm. (G-H) Flow cytometric apoptosis analysis detected the changes of apoptosis of OSA cells with or without miR-3184-5p overexpression. (I-K) Flow cytometry assay showed the regulation of cell cycle by overexpressing miR-3184-5p. Data are expressed as mean ± SD (n = 3) (*, p < 0.05; **, p < 0.01; ***, p < 0.001).
About this article
Cite this article
Tan, J., Yang, B., Zhong, H. et al. Circular RNA circEMB promotes osteosarcoma progression and metastasis by sponging miR-3184-5p and regulating EGFR expression. Biomark Res 11, 3 (2023). https://doi.org/10.1186/s40364-022-00442-9