- Open Access
Hsa_circ_0001485 promoted osteogenic differentiation by targeting BMPR2 to activate the TGFβ-BMP pathway
Stem Cell Research & Therapy volume 13, Article number: 453 (2022)
Circular RNAs (circRNAs) are a new type of stable noncoding RNA and have been proven to play a crucial role in osteoporosis. This study explored the role and mechanism of hsa_circ_0001485 in osteogenic differentiation.
Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis and Gene Ontology (GO) enrichment analysis were performed according to the previous sequencing data in human bone marrow mesenchymal stem cells (BMSC) before and after the induction of osteogenic differentiation on the differentially expressed circRNAs, to screen out signaling pathways associated with osteogenic differentiation. The hFOB 1.19 cells were used to verify the function and mechanism of specific circRNAs in osteogenic differentiation. Additionally, small interfering fragments and overexpression plasmids were used to determine the role of specific circRNAs during osteogenic differentiation. Furthermore, pull-down experiments and mass spectrometry were performed to determine the proteins that bind to specific circRNAs.
The KEGG and GO enrichment analyses showed that the TGFβ-BMP signaling pathway was related to the osteogenic differentiation process, and four circRNAs were associated with the pathway. The quantitative polymerase chain reaction analysis revealed that hsa_circ_0001485 expression was increased during the osteogenic differentiation process of BMSCs. Knockdown of hsa_circ_0001485 suppressed the activity of the alkaline phosphatase enzyme and the expression of RUNX2, osteopontin, and osteocalcin in the osteogenic hFOB 1.19 cells, whereas overexpression of hsa_circ_0001485 promoted their expression. Additionally, we found that hsa_circ_0001485 and BMPR2 targeted binding to activate the TGFβ-BMP signaling pathway and promoted osteogenic differentiation through mass spectrometry analysis.
This study demonstrates that hsa_circ_0001485 is highly expressed in the osteogenic hFOB 1.19 cells, which activate the TGFβ-BMP pathway through targeted binding of BMPR2, and plays a positive role in regulating osteogenic differentiation.
Osteoporosis (OP) is a systemic bone disease in which bone density and mass decrease, leading to fractures [1, 2]. With the rapid improvement in the social economy, the aging of the population is developing rapidly, increasing the incidence of chronic diseases related to aging, such as osteoporosis, causing a heavy economic burden . More seriously, osteoporotic fractures are closely related to the high risk of subsequent fractures and mortality . However, the etiology and pathogenesis of OP are not fully understood.
Studies have manifested that the main pathogenesis is the disruption of the coupling balance between osteoblast-mediated bone formation and osteoclast-mediated bone resorption, resulting in more significant bone resorption than bone formation [5, 6]. Therefore, promoting osteoblast proliferation and differentiation can effectively prevent and treat OP. The BMSCs, as a class of multipotential stem cells, are the precursor cells of osteoblasts and can differentiate into osteoblasts under certain induction conditions [7, 8]. Several studies have also revealed that BMSCs play a key role in OP. For example, Let-7a-5p attenuates osteogenesis of BMSCs in postmenopausal OP mice ; downregulation of TP53INP2 prevents osteogenic differentiation of BMSCs during OP ; the miR-291a-3p accelerates osteogenic differentiation of BMSCs in dexamethasone-induced OP ; the SIRT3 affects BMSCs senescence and senescence OP . Therefore, further investigation of the molecular mechanisms involving the differentiation of BMSCs into osteoblasts will provide new insights into preventing OP and its therapy.
Circular RNA (circRNA) is a noncoding RNA molecule with a closed-loop structure that is stable in expression, not easily degraded, and widely exists in various eukaryotic organisms [12, 13]. The circRNAs were shown to play an influential role in biological growth and development, disease occurrence and development through miRNA sponging, protein binding, and gene transcription and coding functions [14,15,16]. Recently, researchers have paid attention to the role of circRNAs in osteogenic differentiation . The decreased bone formation ability of the bone marrow mesenchymal stem cells (BMSCs) is one of the major causes of OP. However, the specific mechanism of BMSCs in osteoblast differentiation is unclear. Studies have revealed that hsa_circ_0006215 competitively binds to miR-942-5p and promotes the differentiation of BMSCs and osteogenesis–angiogenesis coupling by regulating RUNX2 . The RUNX2 is an essential molecule in osteodifferentiation and bone development and is responsible for activating osteoblast differentiation marker genes .
Additionally, osteopontin (OPN), alkaline phosphatase (ALP), and osteocalcin (OCN) are the most investigated osteogenic-related genes . Liu et al. found that circRNA AFF4 activated the SMAD1/5 signaling pathway by binding to miR-135a-5p and played a role in the bone differentiation of BMSCs . Han et al. revealed that hsa_circ_0076690 was significantly correlated with bone mineral density (BMD) in patients with OP, with a sensitivity of 79% and a specificity of 85%, and can be used as a potential diagnostic biomarker of OP . Moreover, hsa_circ_0024097 acts as a ceRNA and sponges with miR-376b-3p targeting Yap1 and passes through the Wnt/β-catenin pathway to promote osteogenic differentiation . Although the effect of circRNA on osteogenic differentiation has become a hot topic for researchers, it has only revealed the tip of the iceberg. The sponge mechanism has been used in many previous investigations to explain the role of circRNA in osteogenic differentiation since circRNA functions like another noncoding LncRNA, which can affect different physiological and pathological conditions by directly binding to proteins [24, 25]. However, further elucidation of whether circRNA can directly bind to essential osteogenic proteins to affect osteogenic differentiation still lacks research.
In this study, the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analysis of differentially expressed circRNAs showed that the TGFβ-BMP signaling pathway was involved in the osteogenic differentiation process, in which hsa_circ_0001485 played an indispensable role. The TGFβ-BMP signaling pathway was proven to play a vital function in bone development and postnatal bone homeostasis. Furthermore, it has clinical significance for treating osteoporosis, fracture healing, and osteoarthritis . We also applied mass spectrometry (MS) to analyze the target proteins associated with hsa_circ_0001485 and TGFβ-BMP pathway. We also discovered that hsa_circ_0001485 might promote osteogenic differentiation by activating the TGFβ-BMP pathway through targeted binding with the BMPR2 protein. This study revealed a new mechanism by which circRNAs regulate osteogenic differentiation, laying a foundation for targeted therapy of OP.
Methods and materials
The BMSCs (SALIAI, Guangzhou, China; G02001), MG63 (CRL-1427, ATCC), and hFOB 1.19 (Cellcook, Guangzhou, China; CC4005) were incubated in DMEM/F12 (GIBCO; A4192001) supplemented 10% fetal bovine serum (FBS; GIBCO, 10099141C) in a 5% CO2 incubator at 37 °C. When the cell density reached approximately 80%, they were dissociated with trypsin and inoculated into 6-well plates for other experiments.
The BMSCs were inoculated into 6-well plates and cultured to adhere to the wall. Osteogenic differentiation of BMSCs was induced using an osteoblast differentiation kit (#HuxMA-90021, Cyagen, Guangzhou, China) according to the instructions , and the cells were collected at 0, 7, 14, and 21 days of induction. The medium was changed once in 3 days. The kit included OriCell®Basal Medium For Cell Culture (177 mL; BLDM-03011), OriCell®Fetal Bovine Serum (Superior-Quality) (20 mL; FBSSR-01021), and OriCell®Supplement For Human BMSC Osteogenic Differentiation (3 mL; HUXMX-04021).
Cell transfection and infection
We designed hsa_circ_0001485 interference fragments (AACAGTTTCAGAGTGGCAGAA) to suppress the expression of hsa_circ_0001485 and synthesized them using GenePharma Co. Ltd. (Shanghai, China). First, serum-free medium (DMEM/F12) and interference fragments were successively added to the transfection tube and thoroughly mixed. Subsequently, the transfection agent, Lipofectamine 3000 (Invitrogen, Carlsbad, CA, USA), was added to the same tube. After thorough mixing, the mixture was left to stand for 10–15 min. Next, the mixture was added to the culture plate and replaced after 48 h with a medium containing 10% FBS. Additionally, for the hsa_circ_0001485 overexpression, hsa_circ_0001485 overexpression plasmids were synthesized by RiboBio (Guangzhou, China) and transfected into hFOB 1.19 cells for overexpressing hsa_circ_0001485 (OE-circ1485), and the vector LV003 kindly provided by Forevergen (Guangzhou, China) was used as a negative control. Finally, the transfection efficiency was measured by quantitative polymerase chain reaction (qPCR), and cells were collected for subsequent experiments.
Moreover, Guangzhou Forevergen Bioscience Company further constructed hsa_circ_0001485 knockdown and overexpression lentiviruses and their vector controls (NC1, Sh-hsa_circ_0001485, NC2, and OE-hsa_circ_0001485). The sequence for knockdown was transformed from the siRNA sequence, which was 5′-AACAGTTTCAGAGTGGCAGAA-3′, and the corresponding NC1 sequence was 5′-TTCTCCGAACGTGTCACGTTTC-3′. Lentiviruses (109 TU/mL) and polybrene (5 μg/mL, Sigma) were added to the medium and incubated with BMSCs together for 24 h at a multiplicity of infection of 30. The qPCR measured the transfection efficiency.
Quantitative real-time polymerase chain reaction
Total RNA was extracted using TRIzol reagent (Invitrogen, TR118-500) according to the manufacturer’s instructions. Reverse transcriptase reactions were performed using the Prime Script RT kit (Takara, Japan, RR047A). The reverse transcription polymerase chain reaction (RT–PCR) was performed using qPCR and SYBR Premix Ex Taq TM (Takara, Japan, RR82WR) and an ABI 7500 sequence detection polymerase chain reaction (PCR) system (Applied Biosystems, USA). Different primers were designed to identify circRNA expression (hsa_circ_0006952, hsa_circ_0008532, hsa_circ_0005564, and hsa_circ_0001485). Convergence and divergence primers for hsa_circ_0001485 and GAPDH were used to verify the circRNA circular formed using agarose gel electrophoresis, and then, the divergence primers were used to perform Sanger sequencing for hsa_circ_000148 products. The GAPDH was used as an internal control to normalize the levels of circRNAs. Table 1 presents the primer sequences.
Alizarin red (AR) staining
The method of AR staining has been described previously . First, the cell culture medium was removed and washed with PBS twice. Subsequently, 2 ml of 4% formaldehyde was added and fixed for 30 min. Next, the cells were washed twice with PBS, and 1 ml of 0.2% AR solution (Solarbio, China, G1450) was added to each well for staining for 5 min. After washing with PBS, the images of osteogenesis staining were captured under a microscope. Subsequently, the samples were treated with 1 ml of 10% cetylpyridinium chloride (CPC; Sigma) for 30 min to elute the stain. Finally, 100 μL of the eluted stain was added to a 96-well plate and tested using a spectrophotometer at 550 nm. The calcium content of each sample was analyzed based on the standard curve prepared by ALZ and CPC.
The ALP staining was performed to test the osteogenic ability of the cells. The osteogenic hFOB 1.19 cells were washed with PBS, fixed with 10% formaldehyde for 20 min, and washed with double steaming water twice. The ALP staining (Sigma–Aldrich, C3206) was added and placed in a horizontal shaker for incubation for 1 h in the dark. The cells were rinsed with double steaming water four times for 5 min. The excess double steam water was removed, and the cells were observed under an inverted microscope and photographed for analysis.
ALP enzyme activity determination
The osteogenic hFOB 1.19 cells were rinsed with PBS three times and lysed with RIPA lysis buffer (Beyotime), and the supernatant was collected and added to 96-well plates. Reagent from the ALP Assay kit (Beyotime, China, P0321S) was added to a 96-well plate. After incubation at 37 °C for 30 min, absorbance was measured at 405 nm, and the ALP activity was assessed according to the instructions.
Western blot analysis
In the presence of a protease inhibitor (Beyotime, China, P1005), the osteogenic hFOB 1.19 cells were lysed using RIPA lysis buffer (Beyotime, China, P0013B), and centrifugation at 4 °C (12,000 rpm) for 15 min collected the supernatant. Protein samples were separated by SDS–PAGE electrophoresis and then transferred to polyvinylidene fluoride (PVDF; Millipore, IPVH20200) membranes. The TBST containing 5% skim milk (Beyotime, China, P0216) was used to block the nonspecific binding site for 1 h, and the antibodies were incubated overnight at 4 °C. After incubation with the antimouse IgG (Jackson, 115-035-003) or antirabbit IgG (Jackson, 111-035-003) at room temperature for 1 h, the signal was detected with an enhanced chemiluminescence reagent (Beyotime, China, P0018FM). Table 2 presents the list of the primary antibodies.
RNA pull-down and silver staining
The CircRNA pull-down was performed using the Magnetic RNA–protein pull-down kit (ThermoFisher, 20164, USA), as previously reported . The osteogenic hFOB 1.19 cells transfected with hsa_circ_0001485 overexpression plasmids were lysed and hybridized with the probe and incubated at 70 °C for 5 min. Streptavidin magnetic beads were then incubated with the beads at room temperature for 30 min. Next, the protein binding buffer and total protein were added to a tube containing streptavidin magnetic beads. After rotary incubation at 4 °C for 1.5 h, streptavidin magnetic beads were washed thrice with washing buffer and then incubated with the elution buffer at 37 °C for 15 min. After the RNA was pulled down, the same amount of protein was added to the 10% polypropylene gel. After electrophoresis for 1.5 h, silver staining (Beyotime, China, P0017S) was performed according to the instructions.
The differential bands displayed by silver staining were cut off for MS analysis. The protein bands were cut out and washed using Milli-Q water for 1 min. Next, the protein bands were added to the decolorization solution in a 37 °C constant temperature box for 30 min. After centrifugation, the supernatant was added with 100% ACN and shaken for 30 s until the colloidal particles turned white. Next, the collected liquid was added to 25 mM DTT/50 mM NH4HCO3 and reacted at 56 °C for 30 min. Subsequently, DTT was collected and treated with 55 mM IAA/50 mM NH4HCO3 in the darkness for 30 min. After three washes, the sample was dehydrated with 100% ACN until colloidal particles turned white. After enzyme digestion, the peptide was purified, and the sample was dissolved in solution (0.1% formic acid, 2% acetonitrile). After centrifugation (13,000 rpm at 4 °C for 20 min), the supernatant was collected and identified through MS. The MS analysis was conducted using a nano-ultra‐performance liquid chromatography (nanoUPLC) system (Waters, Milford, MA, USA) and Synapt High Definition Mass Spectrometry (Waters, Milford, MA, USA), which was similar to the previous study .
Bone formation of BMSCs in vivo
The BMSCs infected with lentiviruses (NC1, Sh-hsa_circ_0001485, NC2, and OE-hsa_circ_0001485) were cultured in osteoblast differentiation solution for 7 days before the in vivo study. After being trypsinized and resuspended directly in DMEM (Thermo Fisher Scientific, C11965500BT), BMSCs (5 × 105) were loaded onto 40 mg hydroxyapatite/tricalcium phosphate (HA/TCP; Zimmer) and then implanted into the dorsal subcutaneous space on the two symmetrical sites of 8-week-old BALB/c-nu/nu female mice (Guangzhou Forevergen Bioscience Company; n = 5 per group). Specimens were harvested 8 weeks after implantation, and the animals were euthanized by CO2 asphyxiation. The bone constructs were fixed in 4% paraformaldehyde and then decalcified for 10 days in 10% EDTA (pH 7.4). After decalcification, the specimens were dehydrated and subsequently embedded in paraffin. Sections (5 mm thickness) were stained with hematoxylin, eosin (H&E), and Masson’s trichrome staining according to the standard protocol. The ethical committee approved our experiment with Guangzhou Forevergen Biosciences (Guangzhou, China; Approval No. IACUC-AEWC-F2010008).
Data are presented as the mean ± standard deviation. Statistical analysis was performed using GraphPad Prism 8.0 (GraphPad Software, La Jolla, CA, USA) and SPSS software Version 22.0 (IBM Corporation, New York, USA). All procedures were repeated thrice. A student’s t test was performed to test for statistical significance between two groups; three or more groups were compared using one-way analysis of variance. A value of P < 0.05 and P < 0.01 was considered statistically significant.
Hsa_circ_0001485 was upregulated in the osteogenic differentiation of BMSCs
We induced osteogenic differentiation of BMSCs to investigate the function of circRNAs in the occurrence and development of OP. Next, we examined the differentially expressed circRNAs between the induced and control groups through RNA sequencing. In previous sequencing data, there were 3440 circRNAs in the control group and 3893 circRNAs in the BMSC-14d group, among which 2191 circRNAs were specific to the BMSC-14d group . The KEGG and GO analyses were performed to conduct a bioenrichment analysis of the differential circRNAs unique to the BMSC-14d group. The signaling pathways related to osteogenic differentiation were screened, from which the TGFβ-BMP signaling pathway was selected as the object of exploration (Fig. 1A, B). The GSEA of the TGFβ and BMP signaling pathways showed that both were activated at BMSC-14d compared with BMSC-0d (Fig. 1C). Subsequently, we selected four relevant circRNAs from the TGFβ-BMP signaling pathway, hsa_circ_0006952, hsa_circ_0008532, hsa_circ_0005564, and hsa_circ_0001485 (Table 3), and detected their expressions in BMSCs at 0, 3, 7, and 14 d by qPCR. The data showed that relative to the control (BMSC-0d) group, the level of hsa_circ_0001485 was significantly upregulated in BMSC osteogenic differentiation groups, especially during the induction of 14 days (Fig. 1D), so hsa_circ_0001485 was used as the target circRNA for research. Figure 1E shows the structure of hsa_circ_0001485, located on chromosome 5, MAP3K1, and exons 4–6. The PCR production of cDNA and gDNA demonstrated that the divergent primer circ_0001485 could amplify the loop in cDNA from hFOB1.19 cells (Fig. 1F). Then, the amplified product of the divergent primer in cDNA was sequenced by Sanger sequencing, and the cyclic site was identified to confirm the cyclic structure of hsa_circ_0001485 (Fig. 1G). We investigated the expression of hsa_circ_0001485 in two human osteoblasts (MG63 cells and hFOB 1.19 cells), and the results revealed that hsa_circ_0001485 was significantly increased in hFOB cells compared with undifferentiated BMSC and MG63 cells (Fig. 1H), so hFOB 1.19 cells were used for the function and mechanism of hsa_circ_0001485 in osteogenic differentiation.
Hsa_circ_0001485 facilitated osteogenic differentiation in vitro
The hFOB 1.19 cells were transfected with hsa_circ_0001485 small interfering RNA (si-circ1485) for 48 h. Subsequently, hFOB 1.19 cells were cultured in a differentiation medium to induce osteogenic differentiation for 14 days. The data revealed that the expression of hsa_circ_0001485 was potently suppressed by hsa_circ_0001485 silencing in the osteogenic hFOB 1.19 cells (Fig. 2A). Through AR and ALR staining, we found that the number of calcified nodules and ALP activity in the osteogenic hFOB 1.19 cells was reduced in the si-circ1485 group compared with that in the si-NC group (Fig. 2B, C). Meanwhile, we also disclosed that relative to the si-NC transfection group, transfection of si-circ1485 could dramatically attenuate the expressions of RUNX2, OPN, and OCN proteins in the osteogenic hFOB 1.19 cells (Fig. 2D, E). Then, we constructed the hsa_circ_0001485 overexpression vector (OE-circ1485) and transfected the osteogenic hFOB 1.19 cells. We discovered that the expression level of hsa_circ_0001485 in the overexpression group was higher than that in the control (NC) group (Fig. 2F). Besides, we proved that hsa_circ_0001485 overexpression could potently increase the number of calcified nodules and the activity of ALP in the osteogenic hFOB 1.19 cells (Fig. 2G, H). Consistent with expectations, hsa_circ_0001485 overexpression also markedly enhanced the expression levels of RUNX2, OPN, and OCN proteins in the osteogenic hFOB 1.19 cells (Fig. 2I, J). Based on the above results, it was shown that hsa_circ_0001485 could induce osteogenic differentiation in hFOB 1.19 cells in vitro.
Hsa_circ_0001485 activated the TGFβ-BMP signaling pathway by combining with BMPR2 in the osteogenic hFOB 1.19 cells
To confirm the function of hsa_circ_0001485 in the TGFβ-BMP signaling pathway, we designed a probe to conduct a pull-down experiment on osteogenic hFOB 1.19 cells transfected with OE-1485 (Fig. 3A). A silver staining test was performed on the pull-down products. Differences were found between the experimental probe group (circ1485) and the control probe group (NC; Fig. 3B). We adopted the gelatinization method. We cut the parts with apparent differences in bands for MS analysis. The MS identified 523 proteins, and we noticed a target protein, BMPR2 (the green and orange peaks are the detected BMPR2 peaks), which played an important role in BMP signaling pathway activation (Fig. 3C). As a major type of BMP receptor, BMPR2 combines with multiple BMPs for the downstream BMP signaling pathway activation . We measured the expression of downstream proteins in the BMP signaling pathway after transfection of OE-circ1485 and si-circ1485. The results confirmed that the expression of BMPR2, p-Smad1/5/9 and Smad4, was markedly increased in the hsa_circ_0001485 overexpression group relative to that in the NC group, and these protein levels were dramatically decreased in hsa_circ_0001485 silencing group comparable to that in the si-NC group (Fig. 3D–F).
Hsa_circ_0001485 promoted osteogenic differentiation through the TGFβ-BMP signaling pathway in the osteogenic hFOB1.19 cells
To further query the potential function of BMPR2, we designed three interference fragments of BMPR2. After transfection for 48 h in the osteogenic hFOB 1.19 cells, we found that si-BMPR transfection had strong interference effects, especially si-BMPR2-2 (Fig. 4A). Western blot data revealed that the protein expression level of BMPR2 was notably decreased in the BMPR2 silencing group compared with the si-NC group, further indicating the effectiveness of BMPR2-2 interference in the osteogenic hFOB1.19 cells (Fig. 4B, C). Through ALP and AR staining, it was observed that after interference with BMPR2, the number of calcified nodules and the expression of ALP in the osteogenic hFOB 1.19 cells decreased (Fig. 4D). Additionally, the protein expression of RUNX2, OPN, and OCN was observably attenuated in the si-BMPR2 group relative to that in the si-NC group (Fig. 4E, F), indicating that interfering with BMPR2 suppressed the osteogenic differentiation ability of the osteogenic hFOB 1.19 cells. We further transfected OE-circ1485 into the hFOB 1.19 cells, which were then processed as si-BMPR2 or BMPR2 inhibitor (DMH1) simultaneously to inhibit the expression of BMPR2. The results denoted that after osteogenic induction, hsa_circ_0001485 overexpression significantly upregulated hsa_circ_0001485, and the upregulation of hsa_circ_0001485 expression mediated by hsa_circ_0001485 overexpression was not modulated by BMPR2 silencing, indicating that BMPR2 is a downstream regulatory gene of hsa_circ_0001485. We also discovered that hsa_circ_0001485 overexpression did not affect the mRNA level of BMPR2, suggesting that hsa_circ_0001485 only affected the protein level of BMPR2. (Fig. 4G). Besides, overexpression of hsa_circ_0001485 prominently promoted the number of calcified nodules and the expression of ALP, while the promotion mediated by hsa_circ_0001485 overexpression could also be memorably inhibited by si-BMPR2 or DMH1 (Fig. 4H). Moreover, we discovered that inhibition of BMPR2 notably weakened the increase in BMPR2, p-Smad1/5/9, Smad4, RUNX2, OPN, and OCN expressions induced by hsa_circ_0001485 overexpression (Fig. 4I, J), which indicated that hsa_circ_0001485 facilitated osteogenic differentiation through the TGFβ-BMP pathway.
Hsa_circ_0001485 facilitated osteogenic differentiation in vivo
Finally, we constructed a classic bone formation model in vivo to determine the role of hsa_circ_0001485 in osteoblast differentiation. The BMSCs infected with lentiviruses (NC1, sh-hsa_circ_0001485, NC2, and OE-hsa_circ_0001485) were cultured in an osteogenic induction medium for 7 days before the in vivo study. The BMSCs were then loaded onto scaffolds and implanted into the subcutaneous space of nude mice. In the H&E and Masson staining results, the new bone formation ability was significantly attenuated in the sh-hsa_circ_0001485 group relative to that in the NC1 group. Simultaneously, the new bone formation ability was enhanced considerably in the hsa_circ_0001485 overexpression group comparable to that in the NC2 group (Fig. 5).
As the aging of the population increases, the size of the OP population increases, and it gradually becomes a clinical disease with a large probability of occurrence and poor treatment effects [2, 37]. Multiple signaling pathways modulate osteoblast and osteoclast differentiation and formation. Developing high-throughput sequencing technology provides strong technical support for studying complex diseases, such as osteoporosis. The CircRNA is widely expressed in human tissues and mediates various physiological events, such as organogenesis, tumorigenesis, and organ development . Recent studies have confirmed that the expression profile of circRNAs is changed during the formation of osteoclasts . Besides, circRNA plays a key role in the formation of osteoblasts and osteoclasts and the differentiation of BMSCs through various pathways, thus affecting the bone formation and bone resorption . Also, multiple studies have proved that circRNA plays a crucial role in OP [41,42,43]. In our previous sequencing results, 219 unique circRNAs were identified after osteogenic differentiation of BMSCs . This study performed KEGG and GO analyses to identify the pathways related to bone differentiation. The TGFβ and BMP were considered most closely related to osteogenic differentiation. Focusing on the TGFβ-BMP pathway, hsa_circ_0001485 was the most significantly increased differentially expressed circRNA.
In our study, we first used human primary BMSC cells for sequencing because the cells are closer to the actual human condition. However, due to the great difficulty of primary cell culture, the limited number of passages, and the high cost, we chose human osteogenic-related cells (MG63 and hFOB 1.19 cells) for our functional study. The BMSCs were used as the controls for the study. According to these results, we finally selected hFOB 1.19 cells to validate the function and mechanism of hsa_circ_0001485 during osteogenic differentiation. The hFOB 1.19 is a human immortalized osteoblast cell line that expresses various osteoblast-specific markers, such as ALP, OCN, and type I collagen . In vivo experiments have confirmed that hFOB 1.19 has osteogenic activity and can be used to study the proliferation and differentiation of normal human osteoblasts and the related cytokines [45, 46].
High expression of ALP activity is an early sign of osteoblast differentiation and maturation. Increased ALP activity enhances bone formation and promotes the formation of bone matrix mineralization . As an essential transcription factor for osteoblast differentiation, RUNX2 encourages the expression of osteoblast secretion proteins, OCN and OPN . After inhibiting hsa_circ_0001485 with small interfering RNA, we found that the ALP enzyme activity in human osteoblast hFOB 1.19 cells was significantly reduced, and the number of calcified nodules was also reduced. Simultaneously, RUNX2, OCN, and OPN levels were also reduced. The overexpression of hsa_circ_0001485 in the osteogenic hFOB 1.19 cells enhanced osteogenic differentiation. However, circRNAs, such as circ-ITCH, circ_0006766, and circ_0000020, among others, have been proven to promote the osteogenic differentiation of BMSCs [49,50,51]. Our research highlighted that in osteoblast hFOB 1.19 cells, overexpression of hsa_circ_0001485 improved osteoblast differentiation, which was undoubtedly a good therapeutic direction for treating OP.
The TGFβ-BMP signaling is vital in embryonic bone development and postnatal bone homeostasis . The TGF-βs and BMPs act on the receptor complex and transduce signals to classical Smad-dependent signaling pathways (TGF-β/BMP ligands, receptors, and Smads), and nonclassical Smad-independent signaling pathways (p38 mitogen-activated protein kinase/p38 MAPK) regulate the differentiation of mesenchymal stem cells during bone development, formation, and homeostasis [53, 54]. The TGFβ binds to the TGFβ receptor complex to phosphorylate Smad2/3 and regulates cell functions, while the BMPs bind to the BMP receptor complex to phosphorylate Smad1/5/8 and regulate cell functions, including cell differentiation and growth . Our MS data showed that 523 proteins, such as BMPR2, RRAS, FKBP1A, MAPK3, and MAPK1, are bound to hsa_circ_0001485 and were related to the TGFβ signaling pathway. The BMPR2 was selected for follow-up study as an important receptor in the BMP family, which could promote the binding and enucleation of Smad1/5/9 and Smad4 by phosphorylation, thus enabling osteogenic differentiation and angiogenesis and inhibiting lipid synthesis . In nontraumatic osteonecrosis of the femoral head, miR-100-5p inhibited the BMPR2-Smad1/5/9 signaling pathway by reducing the expression of BMPR2 and destroying the bone formation of hBMSCs . In human dental pulp stem cells, overexpression of circRFWD2 activated the BMP-Smad pathway through the BMPR2 receptor and targeted miR-6817-5p to enhance osteogenic differentiation . The literature review shows no reports on the regulatory relationship between hsa_circ_0001485 and BMPR2. In our research, inhibition of BMPR2 by small interfering RNA and inhibitor significantly inhibited osteogenic differentiation induced by overexpression of hsa_circ_0001485. This occurrence suggested that overexpression of hsa_circ_0001485 activated the BMP-Smad signaling pathway by targeting the binding of BMPR2, thereby promoting the protein levels of RUNX2, OPN, and OCN and enhancing the osteogenic differentiation of the osteogenic hFOB 1.19 cells. Besides, we discovered that hsa_circ_0001485 can only upregulate BMPR2 protein level but has no effect on mRNA level. Recently, it is seen that circRNAs can function by interacting with proteins . Thus, the circRNA-protein interaction modes mainly include altering protein interactions, tethering or isolating proteins, recruiting proteins to chromatin, forming circular RNA–protein–mRNA complexes, and protein translocation or redistribution . In comparison, the specific mechanisms of hsa_circ_0001485 and BMPR2 protein have not been explored, which is also one of this study’s limitations. In future studies, we will investigate further the regulatory relationship between hsa_circ_0001485 and BMPR2 protein.
Currently, the study of OP-related circRNAs is still preliminary due to limitations in bioinformatics, transcriptomics, and proteomics, the small sample size of population studies, inconsistent findings, and lack of corresponding in vivo experimental validation. Additionally, our study was conducted to investigate the possible molecular mechanisms of OP treatment with BMSCs. Studies proved that transplantation of BMSCs, as a new approach, has great potential and application that can not only bypass the side effects of drug therapy but also treat OP at its root. However, circRNA is stable and not easily degraded; therefore, it can be applied as a target molecule to treat OP. Therefore, we speculate that BMSCs transplantation modified by the hsa_circ_0001485 may have potential therapeutic implications for OP in the elderly, which is the focal path of our future research.
Using RNA sequencing in this study, we identified has_circ_0001485 as closely related to osteogenic differentiation. This finding confirms that has_circ_0001485 binds to BMPR2 to activate the TGF-BMP signaling pathway, thereby promoting osteogenic differentiation of osteogenic hFOB 1.19 cells. Targeting the expression of has_circ_0001485 by small molecule drugs may play a crucial role in bone formation, providing a new strategy for improving OP bone loss.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Bone mineral density
Fetal bovine serum
Quantitative real-time polymerase chain reaction
Harding AT, Beck BR. Exercise, osteoporosis, and bone geometry. Sports (Basel). 2017;5(2):29.
Lee CY, Chen PE, Tung TH. Clinical epidemiology of osteoporosis among elderly fishing and agricultural population in Taipei, Taiwan. J Men’s Health. 2020;16(1):e53–62.
Chen J, Hong Z, Zhao C, Bi Q, Qiu B. Associations between polymorphisms of the PDLIM4 gene and susceptibility to osteoporotic fracture in an elderly population of Han Chinese. Biosci Rep. 2019;39(1):BSR20181505.
Liu R, Chao A, Wang K, Wu J. Incidence and risk factors of medical complications and direct medical costs after osteoporotic fracture among patients in China. Arch Osteoporos. 2018;13(1):12.
Blakely KK, Johnson C. New osteoporosis treatment means new bone formation. Nurs Womens Health. 2020;24(1):52–7.
Dempster DW, Zhou H, Rao SD, Recknor C, Miller PD, Leder BZ, et al. Early effects of abaloparatide on bone formation and resorption indices in postmenopausal women with osteoporosis. J Bone Miner Res. 2021;36(4):644–53.
Ma W, Dou Q, Ha X. Let-7a-5p inhibits BMSCs osteogenesis in postmenopausal osteoporosis mice. Biochem Biophys Res Commun. 2019;510(1):53–8.
Qadir A, Liang S, Wu Z, Chen Z, Hu L, Qian A. Senile osteoporosis: the involvement of differentiation and senescence of bone marrow stromal cells. Int J Mol Sci. 2020;21(1):349.
Yang Y, Sun Y, Mao WW, Zhang H, Ni B, Jiang L. Oxidative stress induces downregulation of TP53INP2 and suppresses osteogenic differentiation of BMSCs during osteoporosis through the autophagy degradation pathway. Free Radic Biol Med. 2021;166:226–37.
Li ZH, Hu H, Zhang XY, Liu GD, Ran B, Zhang PG, et al. MiR-291a-3p regulates the BMSCs differentiation via targeting DKK1 in dexamethasone-induced osteoporosis. Kaohsiung J Med Sci. 2020;36(1):35–42.
Guo Y, Jia X, Cui Y, Song Y, Wang S, Geng Y, et al. Sirt3-mediated mitophagy regulates AGEs-induced BMSCs senescence and senile osteoporosis. Redox Biol. 2021;41: 101915.
Li Y, Zheng Q, Bao C, Li S, Guo W, Zhao J, et al. Circular RNA is enriched and stable in exosomes: a promising biomarker for cancer diagnosis. Cell Res. 2015;25(8):981–4.
Salzman J, Circular RNA. Expression: its potential regulation and function. Trends Genet. 2016;32(5):309–16.
Rong D, Sun H, Li Z, Liu S, Dong C, Fu K, et al. An emerging function of circRNA-miRNAs-mRNA axis in human diseases. Oncotarget. 2017;8(42):73271–81.
Zheng Q, Bao C, Guo W, Li S, Chen J, Chen B, et al. Circular RNA profiling reveals an abundant circHIPK3 that regulates cell growth by sponging multiple miRNAs. Nat Commun. 2016;7:11215.
Han B, Chao J, Yao H. Circular RNA and its mechanisms in disease: from the bench to the clinic. Pharmacol Ther. 2018;187:31–44.
Dou C, Cao Z, Yang B, Ding N, Hou T, Luo F, et al. Changing expression profiles of lncRNAs, mRNAs, circRNAs and miRNAs during osteoclastogenesis. Sci Rep. 2016;6:21499.
Ji H, Cui X, Yang Y, Zhou X. CircRNA hsa_circ_0006215 promotes osteogenic differentiation of BMSCs and enhances osteogenesis-angiogenesis coupling by competitively binding to miR-942-5p and regulating RUNX2 and VEGF. Aging (Albany NY). 2021;13(7):10275–88.
Vimalraj S, Arumugam B, Miranda PJ, Selvamurugan N. Runx2: Structure, function, and phosphorylation in osteoblast differentiation. Int J Biol Macromol. 2015;78:202–8.
Bernardi S, Re F, Bosio K, Dey K, Almici C, Malagola M, et al. Chitosan-hydrogel polymeric scaffold acts as an independent primary inducer of osteogenic differentiation in human mesenchymal stromal cells. Materials (Basel). 2020;13(16):3546.
Liu C, Liu AS, Zhong D, Wang CG, Yu M, Zhang HW, et al. Circular RNA AFF4 modulates osteogenic differentiation in BM-MSCs by activating SMAD1/5 pathway through miR-135a-5p/FNDC5/Irisin axis. Cell Death Dis. 2021;12(7):631.
Han S, Kuang M, Sun C, Wang H, Wang D, Liu Q. Circular RNA hsa_circ_0076690 acts as a prognostic biomarker in osteoporosis and regulates osteogenic differentiation of hBMSCs via sponging miR-152. Aging (Albany NY). 2020;12(14):15011–20.
Huang Y, Xiao D, Huang S, Zhuang J, Zheng X, Chang Y, et al. Circular RNA YAP1 attenuates osteoporosis through up-regulation of YAP1 and activation of Wnt/β-catenin pathway. Biomed Pharmacother. 2020;129: 110365.
Guo C, Gao YY, Ju QQ, Zhang CX, Gong M, Li ZL. LINC00649 underexpression is an adverse prognostic marker in acute myeloid leukemia. BMC Cancer. 2020;20(1):841.
Yuan P, Lei L, Dong S, Liu D. Circular RNA hsa_circ_0068033 acts as a diagnostic biomarker and suppresses the progression of breast cancer through sponging miR-659. Onco Targets Ther. 2020;13:1921–9.
Wu M, Chen G, Li YP. TGF-β and BMP signaling in osteoblast, skeletal development, and bone formation, homeostasis and disease. Bone Res. 2016;4:16009.
Zhang S, Liu Y, Zheng Z, Zeng X, Liu D, Wang C, et al. MicroRNA-223 suppresses osteoblast differentiation by inhibiting DHRS3. Cell Physiol Biochem. 2018;47(2):667–79.
Wang X, Mi Y, He W, Hu X, Yang S, Zhao L, et al. Down-regulation of miR-340-5p promoted osteogenic differentiation through regulation of runt-related transcription factor-2 (RUNX2) in MC3T3-E1 cells. Bioengineered. 2021;12(1):1126–37.
Chambonnier G, Roux L, Redelberger D, Fadel F, Filloux A, Sivaneson M, et al. The hybrid histidine kinase LadS forms a multicomponent signal transduction system with the GacS/GacA two-component system in Pseudomonas aeruginosa. PLoS Genet. 2016;12(5): e1006032.
Wu YJ, Huang CH, Hsieh TJ, Tseng WL, Lu CY. Identification of significant protein markers by mass spectrometry after co-treatment of cells with different drugs: an in vitro survey platform. Rapid Commun Mass Spectrom. 2020;34(Suppl 1): e8582.
Liu Z, Liu Q, Chen S, Su H, Jiang T. Circular RNA Circ_0005564 promotes osteogenic differentiation of bone marrow mesenchymal cells in osteoporosis. Bioengineered. 2021;12(1):4911–23.
Qin X-B, Wen K, Wu X-X, et al. MiR-183 regulates the differentiation of osteoblasts in the development of osteoporosis by targeting Smad4. Acta Histochem. 2021;123(7): 151786.
Huang H, Dou L, Song J, et al. CBFA2T2 is required for BMP-2-induced osteogenic differentiation of mesenchymal stem cells. Biochem Biophys Res Commun. 2018;496(4):1095–101.
Liu Z, Liu Q, Chen S, et al. Circular RNA Circ_0005564 promotes osteogenic differentiation of bone marrow mesenchymal cells in osteoporosis. Bioengineered. 2021;12(1):4911–23.
Cao Z, Yang Z, Wei S, et al. Comprehensive circRNA expression profile and function network in osteoblast-like cells under simulated microgravity. Gene. 2020;764: 145106.
Luo L, Zheng W, Lian G, Chen H, Li L, Xu C, et al. Combination treatment of adipose-derived stem cells and adiponectin attenuates pulmonary arterial hypertension in rats by inhibiting pulmonary arterial smooth muscle cell proliferation and regulating the AMPK/BMP/Smad pathway. Int J Mol Med. 2018;41(1):51–60.
Xie Y, Zhang L, Xiong Q, Gao Y, Ge W, Tang P. Bench-to-bedside strategies for osteoporotic fracture: from osteoimmunology to mechanosensation. Bone Res. 2019;7:25.
Cheung CL, Ang SB, Chadha M, Chow ES, Chung YS, Hew FL, et al. An updated hip fracture projection in Asia: the Asian Federation of Osteoporosis Societies study. Osteoporos Sarcopenia. 2018;4(1):16–21.
Choksi P, Jepsen KJ, Clines GA. The challenges of diagnosing osteoporosis and the limitations of currently available tools. Clin Diabetes Endocrinol. 2018;4:12.
Yang Y, Yujiao W, Fang W, Linhui Y, Ziqi G, Zhichen W, et al. The roles of miRNA, lncRNA and circRNA in the development of osteoporosis. Biol Res. 2020;53(1):40.
He T, Liu W, Cao L, Liu Y, Zou Z, Zhong Y, et al. CircRNAs and LncRNAs in osteoporosis. Differentiation. 2020;116:16–25.
Wang H, Zhou K, Xiao F, Huang Z, Xu J, Chen G, et al. Identification of circRNA-associated ceRNA network in BMSCs of OVX models for postmenopausal osteoporosis. Sci Rep. 2020;10(1):10896.
Yu L, Liu Y. circRNA_0016624 could sponge miR-98 to regulate BMP2 expression in postmenopausal osteoporosis. Biochem Biophys Res Commun. 2019;516(2):546–50.
Lu M, Zhao XH. The growth proliferation, apoptotic prevention, and differentiation induction of the gelatin hydrolysates from three sources to human fetal osteoblasts (hFOB 11.9 Cells). Molecules. 2018;23(6):1287.
Cai T, Wu Y, Ke R, Yang J, Ghanem A, Mu X. Overexpression of glypican 3 promotes proliferation, regulates cell cycle progression, and inhibits apoptosis of human fetal osteoblastic cell line 119. J Craniofac Surg. 2017;28(6):1481–5.
Zhang Z, Zhang L, Hua Y, Jia X, Li J, Hu S, et al. Comparative proteomic analysis of plasma membrane proteins between human osteosarcoma and normal osteoblastic cell lines. BMC Cancer. 2010;10:206.
Zhang W, Dong Z, Li D, Li B, Liu Y, Zheng X, et al. Cathepsin K deficiency promotes alveolar bone regeneration by promoting jaw bone marrow mesenchymal stem cells proliferation and differentiation via glycolysis pathway. Cell Prolif. 2021;54(7): e13058.
Park KR, Kim S, Cho M, Yun HM. Limonoid triterpene, obacunone increases runt-related transcription factor 2 to promote osteoblast differentiation and function. Int J Mol Sci. 2021;22(5):2483.
Zhong D, Xu GZ, Wu JZ, Liu H, Tang JY, Wang CG. Circ-ITCH sponges miR-214 to promote the osteogenic differentiation in osteoporosis via upregulating YAP1. Cell Death Dis. 2021;12(4):340.
Guo Z, Xie M, Zou Y, Liang Q, Liu F, Su J, et al. Circular RNA Hsa_circ_0006766 targets microRNA miR-4739 to regulate osteogenic differentiation of human bone marrow mesenchymal stem cells. Bioengineered. 2021;12(1):5679–87.
Zhou R, Miao S, Xu J, Sun L, Chen Y. Circular RNA circ_0000020 promotes osteogenic differentiation to reduce osteoporosis via sponging microRNA miR-142-5p to up-regulate Bone Morphogenetic Protein BMP2. Bioengineered. 2021;12(1):3824–36.
Chen G, Deng C, Li YP. TGF-β and BMP signaling in osteoblast differentiation and bone formation. Int J Biol Sci. 2012;8(2):272–88.
Lei SS, Su J, Zhang Y, Huang XW, Wang XP, Huang MC, et al. Benefits and mechanisms of polysaccharides from Chinese medicinal herbs for anti-osteoporosis therapy: a review. Int J Biol Macromol. 2021;193(Pt B):1996–2005.
Han H, Li M, Liu H, Li H. Electroacupuncture regulates inflammation, collagen deposition and macrophage function in skeletal muscle through the TGF-β1/Smad3/p38/ERK1/2 pathway. Exp Ther Med. 2021;22(6):1457.
Hayrapetyan A, Jansen JA, van den Beucken JJ. Signaling pathways involved in osteogenesis and their application for bone regenerative medicine. Tissue Eng Part B Rev. 2015;21(1):75–87.
Yang W, Zhu W, Yang Y, Guo M, Qian H, Jiang W, et al. Exosomal miR-100-5p inhibits osteogenesis of hBMSCs and angiogenesis of HUVECs by suppressing the BMPR2/Smad1/5/9 signalling pathway. Stem Cell Res Ther. 2021;12(1):390.
Huang X, Pan X, Zhang B, Huang W, Cen X, Liu J, et al. CircRFWD2 promotes osteogenic differentiation of human dental pulp stem cells by targeting miR-6817-5p through BMP-Smad and p38 MAPK pathway. Cell Transplant. 2021;30:9636897211052960.
Du WW, Zhang C, Yang W, Yong T, Awan FM, Yang BB. Identifying and characterizing circRNA–protein interaction. Theranostics. 2017;7(17):4183–91.
Huang A, Zheng H, Wu Z, Chen M, Huang Y. Circular RNA-protein interactions: functions, mechanisms, and identification. Theranostics. 2020;10(8):3503–17.
The Natural Science Foundation of Guangdong Province of China (Grant number 2018A030313694) supported our study.
Ethics approval and consent to participate
The ethical committee approved our experiment with Guangzhou Forevergen Biosciences (Guangzhou, China) (Approval No. IACUC-AEWC-F2010008).
Consent for publication
The authors reported that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Chen, SC., Jiang, T., Liu, QY. et al. Hsa_circ_0001485 promoted osteogenic differentiation by targeting BMPR2 to activate the TGFβ-BMP pathway. Stem Cell Res Ther 13, 453 (2022). https://doi.org/10.1186/s13287-022-03150-1
- Osteogenic differentiation
- TGFβ-BMP pathway