Int J Med Sci 2022; 19(9):1377-1387. doi:10.7150/ijms.69501 This issue
1. Department of Orthopedics, Affiliated Hospital of Guangdong Medical University, Zhanjiang, Guangdong, 524001, China.
2. Oncology Hospital, Affiliated Hospital of Guangdong Medical University, Zhanjiang, Guangdong, 524001, China.
3. Stem Cell Research and Cellular Therapy Center, Affifiliated Hospital of Guangdong Medical University, Zhanjiang, Guangdong, 524001, China.
*These authors contributed equally to this work.
Osteosarcoma is the most common malignant tumor in adolescent bone malignancies. It has the characteristics of a high metastasis rate, high mortality and poor prognosis. As a subclass of endogenous noncoding RNAs, circRNAs have been identified to be related to the occurrence, development and prognosis of different kinds of cancers, but the mechanism of their effect on osteosarcoma is not clear. In the present study, we identified a novel circRNA, hsa_circ_0087302, by RNA-seq, and we found that it was expressed at low levels in osteosarcoma. Using RT-PCR, we confirmed that the expression of hsa_circ_0087302 in osteosarcoma cells was lower than that in osteoblasts. Functional validation experiments revealed that hsa_circ_0087302 overexpression inhibited proliferation, cell cycle, migration, and invasion in osteosarcoma cells. Furthermore, Western blotting experiments demonstrated that hsa_circ_0087302 affected the expression of cell cycle- and Wnt/β-catenin signaling pathway-related proteins. For the first time, we identified that hsa_circ_0087302 may affect the malignant biological behavior of osteosarcoma cells through the Wnt/β-catenin signaling pathway. In summary, hsa_circ_0087302 may provide a new direction for the diagnosis and treatment of osteosarcoma.
Keywords: Hsa_circ_0087302, Osteosarcoma, Oncogenic circRNA, Wnt/β-catenin
Osteosarcoma is the most common malignant bone tumor in adolescents . According to statistics, the incidence rate of osteosarcoma is three parts per million in China, accounting for 5%  of young people's malignant tumors. The incidence rate of osteosarcoma is increasing year by year and is the main cause of cancer death among adolescents. Because of the lack of specific diagnostic markers and detection indicators, approximately 20% of patients have imperceptible micrometastasis at the first diagnosis [3,4]. In recent decades, a comprehensive treatment strategy based on surgery combined with chemotherapy has become the standard treatment for osteosarcoma, and the 5-year survival rate has increased from 70% to 82% . However, surgical resection often affects the motor function of patients and has a high disability rate. Mainstream chemotherapy drugs have large side effects and lack of specificity, resulting in poor treatment effect and prognosis. The high degree of malignancy, metastasis in the early stage and recurrence after treatment are important factors affecting the overall survival of osteosarcoma. Experiments have demonstrated that the development of osteosarcoma is closely related to risk factors, such as age, sex, height and family inheritance [6-8]; however, the exact etiology of osteosarcoma is not clear, and the pathogenesis of osteosarcoma needs to be further clarified. Therefore, for the diagnosis and therapy of osteosarcoma, it is essential to study the molecules related to the pathogenesis of osteosarcoma, identify specific biomarkers and screen clinical early warning and treatment targets.
As an emerging, widespread and evolutionarily conserved non-coding RNA, circular RNA (circRNA) has become a hotspot in life science research in recent years. In 1976, the first circular RNA was found in plant-infected viroid and parainfluenza virus particles , and it was later detected in humans and animals . According to the difference of origin, circular RNAs are divided into the following three categories: exon circRNAs (ecRNAs) [11,12], exon-intron circRNAs (elciRNAs)  and circular intron RNAs (circRNAs) . Existing studies have shown that circRNA has various biological functions, including acting as a miRNA sponge, acting as a protein sponge, interacting with RNA-binding proteins, selective splicing, transcriptional regulation, encoding proteins and encoding tumor-related functional peptides .Among them, the function of circRNA-encoded proteins on the occurrence and development of diseases is controversial, and studies have found emerging functions of circRNA-encoded proteins in human cancers . Circular RNA may be related to the occurrence and progression of many diseases, and it is expected to be a new biomarker of diseases. Recently, a growing amount of evidence has shown that circRNAs play an essential role in the occurrence and development of tumors. Circular RNAs are related to a variety of cancers and diseases, such as laryngeal cancer, colon cancer, lung cancer, liver cancer, renal cancer, cervical cancer, gastric cancer, esophageal cancer, bladder cancer, breast cancer, glioma, leukemia and lymphoma [17-29]. Although an increasing number of circular RNAs have been identified, little is known about the specific mechanism of circular RNA in osteosarcoma. In the present study, we identified a novel circular RNA, namely hsa_circ_0087302, from RNA sequencing data, and we found that it was downregulated in osteosarcoma cells by RT-PCR. Overexpression of hsa_circ_0087302 inhibited cell proliferation, migration and invasion. Through the enrichment analysis of KEGG pathways, we found that hsa_circ_0087302 was closely related to lysine degradation, bacterial invasion of epithelial cells, the ErbB signaling pathway, protein processing in the endoplasmic reticulum, transcriptional misregulation in the cancer pathway and central carbon metabolism in cancer. Cell cycle experiments and Western blot analysis showed that most osteosarcoma cells overexpressing hsa_circ_0087302 were arrested in G0/G1 phase, and hsa_circ_0087302 overexpression affected the expression of cyclin D1, a cycle-associated protein. Because cyclin D1 is the downstream target gene of the Wnt/β-catenin signaling pathway and the Wnt/β-catenin signaling pathway has been confirmed to be involved in the malignant biological behavior of osteosarcoma, we detected the key proteins in the Wnt/β-catenin signaling pathway by Western blot analysis, and we found that overexpression of hsa_circ_0087302 affected their expression level. Thus, hsa_circ_0087302 may affect the biological behavior of osteosarcoma through the Wnt/β-catenin signaling pathway.
Human osteoblasts (hFOB1.19) and human osteosarcoma cells (HOS, MG63, U2OS and 143B) were purchased from the Chinese Academy of Sciences Cell Bank. hFOB1.19 cells were cultured in D-MEM/F-12 medium (Gibco, Carlsbad, CA, USA) containing 10% fetal bovine serum (FBS) (Gibco, Sydney, Australia) at 33.5 °C and 5% CO2. HOS and MG63 cells were cultured in MEM (Gibco, Carlsbad, CA, USA), and U2OS cells were cultured in RPMI-1640 medium (Gibco, Carlsbad, CA, USA). In addition, 143B cells were grown in D-MEM/F-12 medium. All human osteosarcoma cells were cultured in medium containing 10% FBS at 5% CO2 and 37 °C. Since human osteoblasts (hFOB1.19) and human osteosarcoma cells (HOS, MG63, U2OS and 143B) were grew adherently, the cells were digested with trypsin containing 0.25% Trypsin-EDTA.
Total RNA was isolated from cells according to the instructions of the HiPure Total RNA Mini Kit (Magen). The integrity of total RNA was determined by an Agilent 2100 Bioanalyzer (Applied Biosystems, Carlsbad, CA, USA), and the concentration of total RNA was detected by a Qubit 3.0 Fluorometer (Invitrogen, Carlsbad, CA, USA).
The KAPA RNA HyperPrep Kit with RiboErase (HMR) for Illumina® (Kapa Biosystems, Inc., Woburn, MA) was used to remove rRNA, and fragmentation buffer was added to the obtained mRNA to make short fragments. The fragmented mRNA was then used as a template, and six-base random primers (random hexamers) were used to synthesize the first strand of cDNA. To synthesize the second strand of cDNA, buffer, dNTPs, RNase H and DNA polymerase I was added, and it was purified by QiaQuick PCR kit and eluted with EB buffer. End repair was performed, and an adenine base was added. Sequencing adapters were added, and the target size fragments were recovered by agarose gel electrophoresis. A qPCR-based KAPA Biosystems Library Quantification kit (Kapa Biosystems, Inc., Woburn, MA, USA) was used to determine the accurate quantification of the sequencing application. The library preparation work was then completed, and the constructed library was sequenced.
The KEGG pathway enrichment site (http://www.genome.jp/KEGG/) was used to identify the pathway and functional enrichment of the target gene as well as to draw a pathway map.
Total RNA was extracted from cells by using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's instructions, and NanoDrop equipment (Thermo Fisher Scientific, Waltham, MA, USA) was used to quantitatively analyze the extracted total RNA. First-strand cDNA was generated using a cDNA synthesis kit (Takara, Otsu, Japan). TB Green PCR Master Mix (Takara) and ABI StepOne Real-Time PCR System (Applied Biosystems 7500, Foster City, CA, USA) were used for real-time PCR analysis, and the results were analyzed by comparing Ct (2 -ΔΔCt) values. The primer sequences were as follows: hsa_circ_0087302 (forward, 5′-GCTGCTGCTGCTGCCTATGG-3′ and reverse, 5′-CTGGATGCCTGTCAGGTTTGGAG-3′); and GAPDH (forward, 5′-CACCCACTCCTCCACCTTTG-3′ and reverse, 5′-CCACCACCCTGTTGCTGTAG-3′).
Human lentivirus-circ-0087302 and Human lentivirus-NC was purchased from GeneChem (Shanghai, China). Lentiviral vector containing hsa_circ_0087302 was used to construct hsa_circ_0087302 overexpression model (OE-circ-0087302), Lentiviral vector carrying nonsense sequence as negative control (NC). Two lentivirus suspension was added to the six-well plate respectively, which was at a 50% cell density, at a MOI of 50 and cultured for three days. Puromycin (Gibco, Grand Island, NY, USA) was used to screen the successfully transfected osteosarcoma cells for subculture.
Approximately 2×103 stably transfected osteosarcoma cells were seeded into 96-well plates and allowed to adhere. At five time points (0, 24, 48, 72 and 96 h), 10 µl of CCK-8 solution was added to each well and incubated at 37 °C for two hours. The optical density (OD) value of the sample was measured at 450 nm using a spectrophotometer.
Stably transfected OS cells (500/well) were inoculated into each well of a 6-well plate and incubated in a 37 °C incubator. After MG63 cells were incubated for 5 days and U2OS cells were incubated for 7 days, cells were fixed with 20% methanol for 20 minutes and stained with 0.1% crystal violet for 30 minutes. Colonies with more than 50 cells were counted.
Stably transfected osteosarcoma cells in logarithmic growth were collected. For the migration assays, 100 µl of cell suspension (2×104 cells) in culture medium without FBS was inoculated into the upper chamber (Costar Inc., USA). For the invasion assays, 100 μl of Matrigel (Corning) was added at a concentration of 250 μg/ml to the upper chamber followed by incubation at 37 °C for 1 hour, and the cell suspension (4.0×104 cells) in culture medium without FBS was then added to the upper chamber. To induce cell migration, 500 μl of medium containing 20% FBS was added to the lower chamber. After MG63 and U2OS cells were incubated for 24 hours in the Migration assays and 36 hours in the invasion assays, cells in the upper chamber were wiped with cotton swabs, fixed with methanol for 20 minutes and stained with crystal violet for 30 minutes. Cells were observed and imaged using an inverted microscope (Olympus, Tokyo, Japan).
Stably transfected osteosarcoma cells were cultured in six-well plates in a 37 °C incubator for 48 h. Cells were washed three times with PBS after collection and then incubated in 70% ethanol at 4 °C overnight. After washing the cells with PBS, the cells were centrifuged to obtain pellets. 25 μL of propidium iodide (C1052-2; Beyotime, Shanghai, China) and 10 μL of RNase (C1052-3; Beyotime, China) were added to each sample for staining. After incubation at 37 °C for 30 minutes in the dark, the cell cycle was detected by BD FACSCalibur flow cytometer (FACSCalibur, BD, San Jose, CA, USA).
Stably transfected OS cells were incubated in six-well plates for 48 hours and lysed using RIPA buffer (Solarbio, Beijing, China). Determination of protein concentration using a bicinchoninic acid (BCA) protein assay kit (Beyotime). Then samples were subjected to SDS-PAGE (Solarbio, Beijing, China) at 80 V for 180 minutes and subsequently transferred to a PVDF membrane at 250 mA in transfer solution for 150 minutes. The membrane was blocked with 5% skim milk at room temperature for 1 h and then washed for 30 minutes. The membrane was then incubated with primary antibody at 4 °C overnight and then washed for 30 minutes. The membrane was then incubated with the secondary antibody at room temperature for 1 h. After the membrane was washed, a chemiluminescence and fluorescence imaging system (Tanon 5200, Shanghai, China) was used to detect immune response bands. Rabbit monoclonal antibodies against GAPDH [Cell Signaling Technology (CST), Inc., Danvers, MA, USA] were diluted at 1:1000. Mouse monoclonal antibodies against β-catenin, cyclin D1, Dvl 2, Dvl 3, CD44, LRP 6, LEF 1, Axin, c-Myc and c-Jun [Cell Signaling Technology (CST), Inc., Danvers, MA, USA] were diluted at 1:1000. The secondary antibodies (anti-mouse IgG/HRP) were diluted at 1:5000 [Cell Signaling Technology (CST), Inc., Danvers, MA, USA]. The protein levels were normalized by GAPDH for each sample.
GraphPad Prism 7 was used for all statistical analyses. The measurement data were evaluated by unpaired Student's t-test and are expressed as the mean ± standard deviation (SD). A P value < 0.05 was considered statistically significant.
In order to study the circRNA expression in osteosarcoma cells, we compared the circRNA expression of osteosarcoma cell lines (HOS, MG63 and U2OS) and normal human osteoblast cell lines (hFOB1.19) by RNA-sequencing analysis. After data screening (fold change ≥ 1.5, P value < 0.05, and FDR < 0.05), there were 29 upregulated circRNAs and 68 downregulated circRNAs in osteosarcoma cells (Figure 1B). A hierarchical cluster diagram was generated for the differentially expressed circRNAs (Figure 1A). KEGG analysis showed that the differentially expressed circRNAs were related to six different signaling pathways, including Lysine degradation, Bacterial inwasion of epithelial cells, ErbB signaling pathway, Protein processing in endoplasmic reticulum, Transcriptional misregulation in cancer and Central carbon metabolism in cancer. (Figure 1C-D). By combining the list of downregulated expression circRNAs in the high-throughput sequencing results with the analysis results of circRNA microarray probe data in circbase database (http://www.circbase.org/), we found that hsa_circ_0087302 is formed from the head and tail of exon 13-16 of its parent gene TLE4, which is located on chromosome 9 (chr9: 82319697-82324614), with a length of 731 BP. The presence and structure of hsa_circ_0087302 was confirmed by primer-specific amplification followed by Sanger sequencing of the PCR product. Therefore, we select hsa_circ_0087302 as the research object of this experiment (Figure 1E).
To study the effect of hsa_circ_0087302 on osteosarcoma, we first confirmed the reliability of sequencing results by verified the expression of hsa_circ_0087302 in four osteosarcoma cell lines and human osteoblasts by qRT-PCR. As expected, hsa_circ_0087302 was expressed at low levels in osteosarcoma cells (Figure 2A). To clarify the function of hsa_circ_0087302 in vitro, we constructed the hsa_circ_0087302 overexpression model by lentivirus transduction of osteosarcoma cells (HOS, MG63, U2OS and 143B). Subsequently, using puromycin to screen the successful transfected osteosarcoma cells. The bright green fluorescence of transfected osteosarcoma cells confirmed the successful transduction (Figure 2C). Additionally, qRT-PCR experiments was also employed to further quantify the transduction efficiency of osteosarcoma cells (Figure 2B). qRT-PCR results showed that the expression level of hsa_circ_0087302 in osteosarcoma cells after lentivirus transduction was significantly higher than that in control cells (Figure 2B). We selected MG63 and U2OS cells for subsequent analyses.
To explore the effect of hsa_circ_0087302 on the proliferation and clonogenicity of osteosarcoma cells, the CCK-8 assay and colony formation assays were performed when the osteosarcoma cells were in the logarithmic growth phase. The CCK-8 assay results showed that the increased expression of hsa_circ_0087302 decreased the viability of osteosarcoma cells compared with that of control cells after four days of culture (Figure 3A-B) (P value < 0.05). The colony formation assay results were consistent with those of the CCK-8 assay. The number of colony formation decreased with the increase of hsa_circ_0087302 expression level compared with that of control cells (Figure 3C) (P value < 0.05). These results suggested that hsa_circ_0087302 may inhibit tumor growth.
Differentially expressed circRNAs in osteosarcoma cells and osteoblasts as well as enrichment analysis of related pathways. (A) Hierarchical cluster analysis of differentially expressed circRNAs. In the heatmap, red represents high expression, and green represents low expression. The columns represent cells, and the rows represent different circRNAs. (B) Volcano map showing differentially expressed circRNAs. The vertical line represents a multiple of 1.5, and the horizontal line indicates a p value equal to 0.05. (C-D) The KEGG database was used for pathway enrichment analysis of circRNAs. (E) Schematic illustration of hsa_circ_0087302 formation and the results of Sanger sequencing.
The expression of hsa_circ_0087302 is decreased in osteosarcoma cells, and construction of the hsa_circ_0087302 overexpression model. (A) The expression of hsa_circ_0087302 in OS cells was lower than that in osteoblasts. (B-C) The overexpression model was constructed by lentiviral transduction, and the transduction efficiency was verified by fluorescence microscopy and RT-PCR. NC is the negative control. Data represent the mean ± SD (n = 3) (*P < 0.05, **P < 0.01 and ***P < 0.001).
Because decreased proliferation may be caused by an aberrant cell cycle, flow cytometry was used to analyze the influence of hsa_circ_0087302 on the OS cell cycle. Flow cytometry analysis indicated that hsa_circ_0087302 overexpression led to G0/G1 phase accumulation (Figure 4). These findings suggested that hsa_circ_0087302 may participate in the regulation of the cell cycle in osteosarcoma cells.
The high metastasis rate of osteosarcoma patients may be closely related to cell migration and invasion. To verify the effect of hsa_circ_0087302 on the migration and invasion of osteosarcoma cells, we performed Transwell assays. Overexpression of hsa_circ_0087302 decreased the number of migratory and invasive cells compared with the controls (Figure 5A-B) (P value < 0.05), indicating that hsa_circ_0087302 may be involved in the mechanism of tumor metastasis.
Effect of hsa_circ_0087302 overexpression on the proliferation of osteosarcoma cells. (A-B) A CCK-8 assay was used to detect cell viability. (C-D) A colony formation assay was used to detect cell colony formation. The bar chart shows the number of colonies with more than 50 cells. NC is the negative control. Data represent the mean ± SD (n = 3) (*P < 0.05, **P < 0.01 and ***P < 0.001).
Effect of hsa_circ_0087302 on the osteosarcoma cell cycle. Changes in the cell cycle of osteosarcoma cells transduced with lentivirus were detected by flow cytometry. The bar chart shows the change in the G0/G1 phase. NC is the negative control. Data represent the mean ± SD (n = 3) (*P < 0.05, **P < 0.01 and ***P < 0.001).
Overexpression of hsa_circ_0087302 reduced cell migration and invasion. Lentiviral transduction in MG63 and U2OS cells reduced cell migration (A) and invasion (B). All images were acquired using 10x magnification. The bar chart shows the number of cells passing through the surface of the Transwell. NC is the negative control. Data represent the mean ± SD (n = 3) (*P < 0.05, **P < 0.01, ***P < 0.001).
To explore the molecular mechanism of hsa_circ_0087302 on the proliferation, migration and invasion of osteosarcoma cells, according to the results of KEGG database, we detected the key node proteins in Wnt/β-catenin signaling pathway which was closely related to lysine degradation, the results were as expected. We found that the protein levels of β-catenin, cyclin D1, Dvl 2, Dvl 3, CD44, LRP 6, LEF 1, c-Myc and c-Jun decreased but that the protein level of Axin increased with the increase of hsa_circ_0087302 expression (Figure 6A-B). These findings suggested that hsa_circ_0087302 may participate in the biological behavior of osteosarcoma via the Wnt/β-catenin signaling pathway.
In the present study, we identified differentially expressed circRNAs in osteosarcoma cells and normal human osteoblasts through RNA sequencing. Because the relationship between most differentially expressed circRNAs and osteosarcoma is unclear, we studied the specific mechanism between circRNA and osteosarcoma by focusing on hsa_circ_0087302. Compared with human osteoblasts, we verified the low expression of hsa_circ_0087302 in osteosarcoma cells. Hsa_circ_0087302 overexpression inhibited the rapid proliferation of osteosarcoma cells by affecting cell proliferation, colony formation, and at the same time causes changes in the cell cycle. Studies have shown that cell cycle dysregulation is an important factor in the abnormal proliferation of tumor cell . TP53 is a transcription factor that stabilizes and induces cell cycle arrest after genotoxic stress to function as a tumor suppressor, however, mutated TP53 not only causes cell cycle dysregulation, but also enhances the proliferative capacity of osteosarcoma cells, the study found that in a pig model, the mutated TP53 and its circular counterpart play an important role in the occurrence and progression of osteosarcoma and various cancers . In addition, overexpression of hsa_circ_0087302 inhibited cell migration and invasion. Studies have shown that the high metastasis rate of osteosarcoma is closely related to the enhancement of cell migration and invasion ability , which indicates that the abnormal expression of hsa_circ_0087302 is involved in the biological behavior of osteosarcoma metastasis.
Expression level of Wnt/β-catenin signaling pathway-related proteins after overexpression of hsa_circ_0087302. (A) The protein levels of β-catenin, cyclin D1, Dvl 2, Dvl 3, CD44, LRP 6, LEF 1, c-Myc and c-Jun were downregulated by lentiviral transduction in MG63 and U2OS cells, while the protein level of Axin was increased. GAPDH was used as a loading control. The bar chart shows the expression level of each protein in MG63 (B) and U2OS (C) cells compared to control cells. NC: negative control. Data represent the mean ± SD (n = 3) (*P < 0.05, **P < 0.01 and ***P < 0.001).
The occurrence and development of osteosarcoma involves a variety of complex molecular mechanisms and different signal transduction pathways [33,34]. KEGG pathway enrichment analysis showed that the differentially expressed circRNAs were closely related to lysine degradation. Lysine degradation mainly through four pathways, including pipecolic acid pathway, saccharopine pathway, acetyllysine pathway and the lysine-urea cycle . The relationship between lysine acetylation and Wnt/β-catenin signaling pathway is of great significance to the occurrence and progression of tumors. Histone lysine acetylation can regulate the transcription of important biomolecules in the Wnt/β-catenin signaling pathway, while non-histones Lysine acetylation directly alters the function of core molecules in Wnt/β-catenin signaling pathway . Previous studies have shown that the Wnt/β-catenin signaling pathway plays an important role in the occurrence and progression of osteosarcoma [37,39]. Therefore, we detected the expression of several key proteins in the Wnt/β-catenin signaling pathway, including the upstream regulators (LRP 6, Dvl 2 and Dvl 3) and downstream targets (CD44, cyclin D1, LEF 1, c-Myc and c-Jun) of β-catenin. In addition, we also detected Axin which associated with β-catenin degradation. Hsa_circ_0087302 overexpression decreased the expression levels of key proteins, except Axin, in osteosarcoma cells. These findings suggested that tumor cells overexpressing hsa_circ_0087302 inhibit tumor progression by inhibiting the activation of the Wnt/β-catenin signaling pathway. As a single-pass transmembrane receptor, LRP6 interacts with the frizzled (FZD) family and activates the Wnt/β-catenin signaling pathway, which has been demonstrated in many cancers, including breast cancer, prostate cancer, hepatocellular carcinoma, and retinoblastoma [40,42]. Dvl2 and Dvl3 are important mediators of the Wnt/β-catenin signaling pathway signaling transduction pathway. Studies have shown that DVL overexpression is important for the activation of the Wnt/β-catenin signaling pathway and cell growth in tumor cells [43,44]. The decrease of LRP6, Dvl2 and Dvl3 expression inhibited Wnt signal transduction. Axin is involved in the degradation of β-catenin and negatively regulates the Wnt/β-catenin signaling pathway signal transduction pathway and inhibits the proliferation, invasion and tumorigenicity of osteosarcoma . The increase of axin expression leads to the increased degradation of β-catenin and then inhibits the transduction of Wnt/β-catenin signaling pathway. LEF1 is a member of the high mobility group (HMG) DNA-binding protein family of transcription factors, which mainly regulates the transcription of downstream target genes by β-catenin . As a type I transmembrane glycoprotein, CD44 participates in the migration and invasion of osteosarcoma . CyclinD1, as a regulator of the cell cycle, affects the proliferation of tumor cells through cell cycle arrest . As an oncogene, c-Myc has been found to be involved in many biological behaviors of tumor cells, such as cell proliferation, apoptosis and the cell cycle . Moreover, as the downstream target gene of the ErbB signaling pathway and central carbon metabolism in cancer, c-Myc participates in the malignant biological behavior of tumors, which is also in accordance with the KEGG pathway enrichment analysis. C-Jun has been reported to be involved in the biological behavior of migration and invasion of osteosarcoma . Here, we found that hsa_circ_0087302 overexpression may reduce Wnt signal transduction by down regulating the expression of LRP6, Dvl2 and Dvl3, increase the degradation of β-catenin and inhibit the transmission of Wnt signal by causing the upregulation of Axin expression. The downregulation of β-catenin content and abnormal Wnt signal transduction leads to the decrease of β-catenin entry into the nucleus, resulting in the decrease of the expression of transcription factor LEF1 and downstream target gene CD44, cyclin D1, c-Myc and c-Jun by β-catenin. Thus hsa_circ_0087302 mainly affects the progression of osteosarcoma via the Wnt/β-catenin signaling pathway. However, the downstream target and specific binding protein of hsa_circ_0087302 remain unclear, and the specific mechanism of the Wnt/β-catenin signaling pathway in osteosarcoma needs to be further clarified.
In summary, the present study showed that hsa_circ_0087302 may be directly involved in the malignant biological behavior of osteosarcoma, such as proliferation, cell cycle, migration and invasion, through the Wnt/β-catenin signaling pathway. In the near future, hsa_circ_0087302 may become a specific therapeutic target for osteosarcoma, effectively improving the comprehensive diagnosis and treatment of osteosarcoma patients.
This study was supported by National Natural Science Foundation of China (81600698), Science Fund of Guangdong Medical University (GDMUM201809), and Doctoral Research Fund Project of Affiliated Hospital of Guangdong Medical University (BJ201802).
The authors have declared that no competing interest exists.
1. Nagarajan R, Kamruzzaman A, Ness KK. et al. Twenty years of follow-up of survivors of childhood osteosarcoma: A report from the Childhood Cancer Survivor Study. Cancer. 2011;117:625-34
2. Bielack SS, Kempf-Bielack B, Delling G. 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
3. Goorin AM, Schwartzentruber DJ, Devidas M. et al. Presurgical chemotherapy compared with immediate surgery and adjuvant chemotherapy for nonmetastatic osteosarcoma: Pediatric Oncology Group Study POG-8651. J Clin Oncol. 2003;21(8):1574-1580
4. Duchman KR, Gao Y, Miller BJ. Prognostic factors for survival in patients with high-grade osteosarcoma using the Surveillance, Epidemiology, and End Results (SEER) program database. Cancer Epidemiol. 2015;39:593-599
5. Chiappetta C, Mancini M, Lessi F. et al. Whole-exome analysis in osteosarcoma to identify a personalized therapy [J]. Oncotarget. 2017;8:80416-80428
6. Ottaviani G, Jaffe N. The epidemiology of osteosarcoma. Cancer Treat Res. 2009;152:3-13
7. Zhang C, Morimoto LM, de Smith AJ. et al. Genetic determinants of childhood and adult height associated with osteosarcoma risk. Cancer. 2018;124(18):3742-3752
8. Miller CW, Ikezoe T, Krug U. et al. Mutations of the CHK2 gene are found in some osteosarcomas, but are rare in breast, lung, and ovarian tumors. Genes Chromosomes Cancer. 2002;33(1):17-21
9. Sanger HL, et al. Viroids are Single-Stranded Covalently Closed Circular RNA Molecules Existing as Highly Base-Paired Rod-Like Structures.Proceedings of the National Academy of Sciences of the United States of America. 1976: 73(11):3852-3856
10. Memczaks JM, et al. Circular RNAs are a large class of animal RNAs with regulatory potenc.Nature.2013.495(7441):333-338
11. Zhang XO, Wang HB, Zhang Y. et al. Complementary sequence-mediated exon circularization. Cell. 2014;159(1):134-147
12. Chen I, Chen CY, Chuang TJ. Biogenesis, identification, and function of exonic circular RNAs. Wiley Interdiscip Rev RNA. 2015;6(5):563-579
13. Li Z, Huang C, Bao C. et al. Exon-intron circular RNAs regulate transcription in the nucleus. Nat Struct Mol Biol. 2015;22(3):256-264
14. Zhang Y, Zhang XO, Chen T. et al. Circular intronic long noncoding RNAs. Mol Cell. 2013;51(6):792-806 doi: 10.1016/j.molcel.2013.08.017
15. Barrett S P, Salzman J. Circular RNAs: analysis, expression and potential function [J]. Development. 2016;143:1838-1847
16. Lei M, Zheng G, Ning Q. et al. Translation and functional roles of circular RNAs in human cancer. Mol Cancer. 2020;19(1):30. Published 2020 Feb 15
17. Su H, Lin F, Deng X. et al. Profiling and bioinformatics analyses reveal differential circular RNA expression in radioresistant esophageal cancer cells. J Transl Med. 2016;14:225
18. Guo JN, Li J, Zhu CL, et al.Comprehensive profile of differentially expressed circular RNAs reveals that hsa_circ_0000069 is upregulated, promotes cell proliferation, migration, invasion in colorectal cancer. Onco Targets Ther. 2016;9:7451-8.
19. Yao JT, Zhao SH, Liu QP. et al. Over-expression of circRNA_100876 in non-small cell lung cancer and its prognostic value. Pathol Res Pract. 2017;213:453-Wnt /β-catenin signaling pathway
20. Shang X, Li G, Liu H. et al. Comprehensive Circular RNA Profiling Reveals That hsa_circ_0005075, a New Circular RNA Biomarker, Is Involved in Hepatocellular Crcinoma Development. Medicine (Baltimore). 2016;95:e3811
21. Wang K, Sun Y, Tao W. et al. Androgen receptor (AR) promotes clear cell renal cell carcinoma (ccRCC) migration and invasion via altering the circHIAT1/miR-195-5p/29a-3p/29c-3p/CDC42 signals. Cancer Lett. 2017;394:1-12
22. Hsiao KY, Lin YC, Gupta SK. et al. Non-coding effects of circular RNA CCDC66 promote colon cancer growth and metastasis. Cancer Res. 2017;77:2339-2350
23. Chen J, Li Y, Zheng Q. et al. Circular RNA profile identifies circPVT1 as a proliferative factor and prognostic marker in gastric cancer. Cancer Lett. 2017;388:208-219
24. Su H, Lin F, Deng X. et al. Profiling and bioinformatics analyses reveal differential circular RNA expression in radioresistant esophageal cancer cells. J Transl Med. 2016;14:225
25. Zhong Z, Lv M, Chen J. Screening differential circular RNA expression profiles reveals the regulatory role of circTCF25-miR-103a-3p/miR-107-CDK6 pathway in bladder carcinoma. Sci Rep. 2016;6:30919
26. Du WW, Fang L, Yang W. et al. Induction of tumor apoptosis through a circular RNA enhancing Foxo3 activity. Cell Death Differ. 2017;24(2):357-370
27. Yang P, Qiu Z, Jiang Y. et al. Silencing of cZNF292 circular RNA suppresses human glioma tube formation via the Wnt/beta-catenin signaling pathway. Oncotarget. 2016;7:63449-63455
28. Li W, Zhong C, Jiao J. et al. Characterization of hsa_circ_0004277 as a New Biomarker for Acute Myeloid Leukemia via Circular RNA Profile and Bioinformatics Analysis. Int J Mol Sci. 2017:18.
29. Hu Y, Zhao Y, Shi C, et al.A circular RNA from APC inhibits the proliferation of diffuse large B-cell lymphoma by inactivating Wnt/β-catenin signaling via interacting with TET1 and miR-888.Aging (Albany NY). 2019;11(19):8068-8084
30. Williams GH, Stoeber K. The cell cycle and cancer. J Pathol. 2012;226(2):352-364
31. Niu G, Hellmuth I, Flisikowska T. et al. Porcine model elucidates function of p53 isoform in carcinogenesis and reveals novel circTP53 RNA. Oncogene. 2021;40(10):1896-1908
32. Wang X, Qin G, Liang X. et al. Targeting the CK1α/CBX4 axis for metastasis in osteosarcoma. Nat Commun. 2020;11(1):1141. Published 2020 Feb 28
33. Adamopoulos C, Gargalionis AN, Basdra EK. et al. Deciphering signaling networks in osteosarcoma pathobiology. Exp Biol Med (Maywood). 2016;241:1296-305
34. Shaikh AB, Li F, Li M. et al. Present Advances and Future Perspectives of Molecular Targeted Therapy for Osteosarcoma. Int J Mol Sci. 2016;17:506
35. Vianey-Liaud C, Divry P, Poinas C. et al. Lysine metabolism in man. Ann Biol Clin (Paris). 1991;49(1):18-26
36. You H, Li Q, Kong D. et al. The interaction of canonical Wnt/β-catenin signaling with protein lysine acetylation. Cell Mol Biol Lett. 2022;27(1):7. Published 2022 Jan 15
37. Fang F, VanCleave A, Helmuth R. et al. Targeting the Wnt/β-catenin pathway in human osteosarcoma cells. Oncotarget. 2018;9(95):36780-36792
38. Lin CH, Ji T, Chen CF. et al. Wnt signaling in osteosarcoma. Adv Exp Med Biol. 2014;804:33-45
39. Li C, Shi X, Zhou G, et al.The canonical Wnt-beta-catenin pathway in development, chemotherapy of osteosarcoma. Front Biosci (Landmark Ed). 2013;18:1384-1391.
40. Roslan Z, Muhamad M, Selvaratnam L. et al. The Roles of Low-Density Lipoprotein Receptor-Related Proteins 5, 6, and 8 in Cancer: A Review. J Oncol. 2019;2019:4536302
41. Lu W, Li Y. Salinomycin suppresses LRP6 expression and inhibits both Wnt/β-catenin and mTORC1 signaling in breast and prostate cancer cells. Journal of Cellular Biochemistry. 2014;115(10):1799-1807
42. Ma J, Lu W, Chen D, et al.Role of Wnt Co-Receptor LRP6 in Triple Negative Breast Cancer Cell Migration, Invasion. J Cell Biochem. 2017;118(9):2968-2976.
43. Uematsu K, He B, You L. et al. Activation of the Wnt pathway in non-small cell lung cancer: evidence of dishevelled overexpression. Oncogene. 2003;22(46):7218-7221
44. Zhang H, Yang X, Zhu L. et al. ASPM promotes hepatocellular carcinoma progression by activating Wnt/β-catenin signaling through antagonizing autophagy-mediated Dvl2 degradation. FEBS Open Bio. 2021;11(10):2784-2799
45. Lv YF, Dai H, Yan GN. et al. Downregulation of tumor suppressing STF cDNA 3 promotes epithelial-mesenchymal transition and tumor metastasis of osteosarcoma by the Wnt/GSK-3β/β-catenin/Snail signaling pathway. Cancer Lett. 2016;373(2):164-173
46. Chen J, Liu G, Wu Y. et al. CircMYO10 promotes osteosarcoma progression by regulating miR-370-3p/RUVBL1 axis to enhance the transcriptional activity of β-catenin/LEF1 complex via effects on chromatin remodeling. Mol Cancer. 2019;18(1):150
47. Chen X, Zhong L, Li X. et al. Down-regulation of microRNA-31-5p inhibits proliferation and invasion of osteosarcoma cells through Wnt/β-catenin signaling pathway by enhancing AXIN1. Exp Mol Pathol. 2019;108:32-41
48. Dai G, Zheng D, Wang Q. et al. Baicalein inhibits progression of osteosarcoma cells through inactivation of the Wnt/β-catenin signaling pathway. Oncotarget. 2017;8(49):86098-86116
49. Hua Q, Jin M, Mi B. et al. LINC01123, a c-Myc-activated long non-coding RNA, promotes proliferation and aerobic glycolysis of non-small cell lung cancer through miR-199a-5p/c-Myc axis. J Hematol Oncol. 2019;12:91
50. Su Z, Yang B, Zeng Z. et al. Metastasis-associated gene MAPK15 promotes the migration and invasion of osteosarcoma cells via the c-Jun/MMPs pathway. Oncol Lett. 2020;20(1):99-112
Corresponding authors: Jing Huang, Ph.D., Oncology Center, Affiliated Hospital of Guangdong Medical University, Zhanjiang, Guangdong 524001, China, E-mail: hj841023com; and Hao Lin, Ph.D., Department of Orthopedics, Affiliated Hospital of Guangdong Medical University, Zhanjiang, Guangdong 524001, China, E-mail: puningworkercom.