Skip to main content

Circular RNA CDR1as regulates osteoblastic differentiation of periodontal ligament stem cells via the miR-7/GDF5/SMAD and p38 MAPK signaling pathway

Abstract

Background

Periodontal ligament stem cells (PDLSCs) are considered as candidate cells for the regeneration of periodontal and alveolar bone tissues. Antisense to the cerebellar degeneration-related protein 1 transcript (CDR1as), which is a newly discovered circular RNA (circRNA), has been reported to act as an miR-7 sponge and to be involved in many biological processes. Here, we investigated the potential roles of CDR1as and miR-7 in the osteogenic differentiation of PDLSCs.

Methods

The expression pattern of CDR1as and miR-7 in PDLSCs during osteogenesis was detected by quantitative reverse-transcription polymerase chain reaction (qRT-PCR). Then we overexpressed or knocked down CDR1as or miR-7 to confirm whether they were involved in the regulation of osteoblast differentiation in PDLSCs. Alkaline phosphatase (ALP) and alizarin red S (ARS) staining were used to detect the activity of osteoblasts and mineral deposition. Furthermore, a dual luciferase reporter assay was conducted to analyze the binding of miR-7 to growth differentiation factor (GDF)5. To further verify the role of CDR1as in osteoblast differentiation, we conducted animal experiments in vivo. New bone formation in specimens was analyzed by microcomputed tomography (micro-CT), hematoxylin and eosin staining, and immunofluorescence staining.

Results

We observed that CDR1as was significantly upregulated during the osteogenic differentiation, whereas miR-7 was significantly downregulated. Moreover, knockdown of CDR1as and overexpression of miR-7 inhibited the ALP activity, ARS staining, and expression of osteogenic genes. Overexpression of miR-7 significantly reduced the activity of luciferase reporter vectors containing the wild-type, but not the mutant, 3’ untranslated region (UTR) sequence of GDF5. Furthermore, knockdown of GDF5 partially reversed the effects of miR-7 inhibitor on osteoblast differentiation. Downregulation of CDR1as or GDF5 subsequently inhibited phosphorylation of Smad1/5/8 and p38 mitogen-activated protein kinases (MAPK), while upregulation of miR-7 decreased the level of phosphorylated Smad1/5/8 and p38 MAPK. In vivo, CDR1as knockdown lead to less bone formation compared with the control group as revealed by micro-CT and the histological analysis.

Conclusions

Our results demonstrated that CDR1as acts as a miR-7 inhibitor, triggering the upregulation of GDF5 and subsequent Smad1/5/8 and p38 MAPK phosphorylation to promote osteogenic differentiation of PDLSCs. This study provides a novel understanding of the mechanisms of osteogenic differentiation, and suggests a potential method for promoting bone formation.

Background

The periodontal ligament is a soft tissue located in the alveolar socket that supports the teeth, buffers masticatory force, and contributes to tooth nutrition, alveolar bone remodeling, and orthodontic tooth movement [1, 2]. Periodontal ligament stem cells (PDLSCs), which are derived from periodontal ligament tissue, possess the properties of mesenchymal stem cells as indicated by their capacity for self-renewal and multilineage differentiation [3]. PDLSCs have demonstrated the ability to form cementum/PDL-like structures, bone tissues, periodontal ligaments, peripheral nerves, and blood vessels in vivo [4, 5]. Owing to the advantages of their availability, culture expansion, and strong osteogenesis capacity, PDLSCs are considered good candidates for periodontal and alveolar bone regeneration [6, 7]. Thus, understanding the molecular mechanisms underlying the osteogenic differentiation of PDLSCs could facilitate the development of regenerative therapies for periodontal and bone diseases.

As novel noncoding RNAs, circular RNAs (circRNAs), unlike linear RNAs, can form a covalently closed continuous loop with their 3′ heads and 5′ tails bound together [8, 9]. CircRNAs have been detected in multiple organisms; their unique closed loop structure makes them more resistant to digestion by RNAse R, and they are ideal biomarkers for diagnosis [9]. CircRNAs are reportedly expressed in a spatial- and temporal-specific manner [10, 11]. A growing body of research suggests that circRNAs participate in embryonic development, cellular activities, and many human diseases [12, 13]. However, few studies have focused on the relationship between circRNAs and osteoblastic differentiation of stem cells. In our previous study, using RNA sequencing, we found that the expression patterns of many circRNAs were significantly altered during osteoblastic differentiation of PDLSCs [14], indicating that circRNAs play a critical role in the osteogenesis of PDLSCs.

Functionally, some circRNAs harbor microRNA (miRNA) binding sites and usually act as miRNA sponges [15]. Recently, antisense to the cerebellar degeneration-related protein 1 transcript (CDR1as), also called circular RNA sponge for miR-7 (ciRS-7), was confirmed to have approximately 70 conserved miR-7 binding sites and to act as an miR-7 sponge to perform biological functions [11, 15, 16]. Overexpression of CDR1as caused developmental defects in the midbrain of embryonic zebrafish; miR-7 inhibition also induced a similar phenotype, whereas overexpression of miR-7 partially rescued these defects [11]. CDR1as was also found to regulate the transcription and secretion of insulin by blocking miR-7 [17]. Moreover, CDR1as promotes the proliferation and invasion of several types of carcinoma cells and is associated with many human diseases [18,19,20,21]. However, there have been no studies on the potential role of CDR1as in the osteoblastic differentiation of stem cells.

Our previous RNA-sequencing results showed that CDR1as was significantly upregulated in PDLSCs undergoing osteoblast differentiation, demonstrating that CDR1as played an important role in the osteoblast differentiation process. Here, we investigated the role of CDR1as in the osteogenic differentiation of PDLSCs. CDR1as knockdown significantly inhibited osteoblast differentiation of PDLSCs in vitro and in vivo, while overexpression of miR-7 also induced a similar change. Mechanistically, CDR1as acts as an miR-7 sponge, leading to the upregulation of growth differentiation factor (GDF)5 and the phosphorylation of Smad1/5/8 and p38 mitogen-activated protein kinases (p38 MAPK). Our results reveal novel mechanisms underlying osteogenic differentiation, and suggest potential therapeutic strategies for periodontal tissue and bone regeneration.

Methods

Cell culture and induction

Healthy premolars were collected from three donors who underwent extractions for orthodontic reasons. This study was approved by the Ethics Committee of Peking University School of Stomatology (PKUSSIRB-201837096). Extracted teeth were rinsed with phosphate-buffered saline (PBS) supplemented with 10% penicillin/streptomycin (Gibco, Grand Island, NY, USA). The periodontal ligament was carefully scraped from the middle third of tooth roots and digested with collagenase and trypsin (Gibco) for 60 min, then incubated in growth medium (GM) containing α-modified Eagle’s medium (α-MEM) supplemented with 10% fetal bovine serum (Gibco) and 1% penicillin/streptomycin in a humidified 5% CO2 atmosphere at 37 °C. The primary PDLSCs migrated outwards from the tissues approximately 1 week later. The cells were passaged using 0.25% trypsin and further expanded until passage 4 for subsequent experiments. Osteogenic differentiation of PDLSCs was induced after they reached 70–80% confluence. The osteogenic induction medium (OM) was standard GM supplemented with 100 nM dexamethasone, 200 μM l-ascorbic acid, and 10 mM β-glycerophosphate (Sigma-Aldrich, St. Louis, MO, USA). The medium was changed every 2 days.

RNA oligoribonucleotides

The RNA oligoribonucleotides used in this study, including miR-7 mimic, miR-7 inhibitor, the small interfering RNAs (siRNAs) targeting CDR1as (si-CDR1as) or GDF5 (si-GDF5), and the corresponding miRNA control (miR-NC) and siRNA control (si-NC), were purchased from GenePharma Co. (Shanghai, China). The sequences of these RNA oligoribonucleotides are listed in Additional file 1 (Table S1).

Transfection assay

Transfection was conducted when cells reached 70–80% confluence using Lipofectamine 3000 (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s protocol. The miR-7 mimic, miR-7 inhibitor, si-CDR1as, si-GDF5, and corresponding negative controls were transfected separately at 100 nM. The cells were collected 48 and 72 h after transfection for mRNA and protein analysis, respectively.

Alkaline phosphatase (ALP) staining and activity

ALP staining was performed using a NBT/BCIP staining kit (CoWin Biotech, Beijing, China) as previously described [22]. After osteogenic induction for 7 days, cells were fixed in 4% paraformaldehyde for 10 min, and then ALP staining was performed following the manufacturer’s instructions.

A commercialized ALP activity colorimetric assay kit (BioVision, Milpitas, CA, USA) was used to analyze ALP activity. The cultured cells were washed with cold PBS, then lysed with 1% Triton X-100 (Sigma-Aldrich) and scraped into distilled water. The absorbance was detected at 405 nm. Total protein concentrations were determined by the bicinchoninic acid (BCA) method using the Pierce protein assay kit (Thermo Fisher Scientific, Waltham, MA, USA). ALP activity was calculated from the absorbance levels relative to the protein concentration.

Alizarin red S (ARS) staining and quantification

Mineralized nodule formation was determined by ARS staining, as described previously [14]. After osteogenic incubation for 14 days, the cells were fixed in 95% ethanol for 30 min at room temperature. Subsequently, the cells were washed with distilled water and then stained with 0.1% ARS (pH = 4.2; Sigma-Aldrich) for 20 min. To quantitatively evaluate the mineralized nodules, the stain was dissolved in 1 mL 10% cetylpyridinium chloride (Sigma-Aldrich) for 1 h and the absorbance at 570 nm was detected by spectrophotometric methods. The intensity of ARS was normalized to the total protein concentration.

RNA extraction and quantitative reverse-transcription polymerase chain reaction (qRT-PCR)

Total RNA was isolated using TRIzol reagent (Invitrogen) according to the manufacturer’s protocol. Total RNA (1 μg) was reverse-transcribed into cDNA using a cDNA Reverse Transcription Kit (Applied Biosystems, Foster City, CA, USA). Random primers were used to analyze circRNAs. qRT-PCR was performed using SYBR Green Master Mix on the ABI Prism 7500 Real-Time PCR System (Applied Biosystems), as described previously [23]. The primers used for CDR1as, miR-7, GDF5, bone morphogenetic protein 2 (BMP2), runt-related transcription factor 2 (RUNX2), osteocalcin (OCN), U6 (internal control for miRNAs), and glyceraldehyde 3-phosphate dehydrogenase (GAPDH; internal control for mRNAs and circRNAs) are listed in Additional file 1 (Table S2). The 2–ΔΔCT method was used to analyze gene expression.

Western blot analysis

Radio immunoprecipitation assay (RIPA) lysis buffer was used to extract total protein from cultured cells. Equal quantities of proteins were separated by 12% sodium dodecyl sulfate–polyacrylamide gel electrophoresis and transferred onto polyvinylidene difluoride (PVDF) membranes (Millipore, Billerica, MA, USA). Primary antibodies against RUNX2 (Abcam, Cambridge, UK), GDF5 (Abcam), phosphorylated Smad1/5/8 (pSmad1/5/8; Cell Signaling Technology, Beverly, MA, USA), p38 MAPK (Cell Signaling Technology), phosphorylated p38 MAPK (p-p38 MAPK; Cell Signaling Technology), and β-ACTIN (Abcam) were diluted 1:1000. The PVDF membranes were incubated with these primary antibodies overnight at 4 °C. After three washes with TBST, the PVDF membranes were incubated with corresponding secondary antibodies (1:10,000, Cell Signaling Technology). ImageJ software (http://rsb.info.nih.gov/ij/) was used to quantify the band intensity obtained by Western blot analysis. The signal of each target band was normalized to that of the β-ACTIN band.

Dual luciferase reporter assay

Luciferase assays were performed as described previously [24]. Briefly, 40 ng GDF5 luciferase reporter plasmids were transfected into HEK 293 T cells together with 100 nM miR-NC or miR-7 mimic using Lipofectamine 3000. After transfection for 24 h, the Renilla and firefly luciferase activities were measured separately using the Dual Luciferase Reporter Assay System (Promega, Beijing, China) following the manufacturer’s instructions. The light intensity from Renilla luciferase was normalized to that of firefly luciferase and expressed as fold-induction relative to the basal activity.

In-vivo bone formation assay

We used polylactic-co-glycolic acid (PLGA; Melone, Dalian, China) as a scaffold material. The scaffolds were prepared as thin circular slices approximately 4 mm in diameter. PDLSCs were seeded in the six-well plate at a density of 2 × 105 cells per well. We transiently transfected si-CDR1as and si-NC into PDLSCs. The cells were seeded in the scaffold material at a density of 5 × 105 cells each group and then transplanted into the mouse skull defect area.

The critical-sized mouse calvarial defect model was constructed as previously described [25, 26]. The operation was performed on 60-day-old adult male nude mice (BALB/c) under general anesthesia. After removing the pericranium, nonhealing, critical-sized calvarial defects (4 mm) were created using a dental drill. The seeded scaffolds were gently placed in the defect area, and the skin incision was closed with 5–0 Vicryl sutures. After 8 weeks, the skull tissues of animals were harvested and fixed in 4% polyoxymethylene at 4 °C.

Microcomputed tomography (CT) analysis

New bone formation in specimens was measured by micro-CT analysis using a high-resolution Inveon Micro-CT (Siemens, Munich, Germany). Images were acquired at an effective pixel size of 8.99 μm, voltage of 80 kV, current of 500 μA, and exposure time of 1500 ms. All samples were placed in the same container and scanned with uniform parameters. The specimens were scanned through a 360° rotation in 360 equiangular steps. Files were reconstructed using Inveon Research Workplace 3.0 software (Siemens). Relevant parameters, including bone mineral density (BMD, mg/cm3) and the ratio of new bone volume to existing tissue volume (BV/TV), were calculated.

Hematoxylin and eosin (H&E) staining and immunofluorescence staining analysis

The skull specimens were decalcified for 1 month using 10% ethylenediaminetetraacetic acid (EDTA; pH = 7.4), as previously described [27]. The EDTA solution was changed every 2 days. The specimens were then washed, dehydrated, and embedded in paraffin. Sections were cut at 7 μm and stained with standard H&E. The images were captured using an Olympus BX51 light microscope equipped with an Olympus DP70 camera (Olympus Co., Tokyo, Japan). Sections were also assessed by immunofluorescence analysis, as described previously [22]. The specimens were blocked with 3% goat serum albumin (ZSGB-BIO, Beijing, China) for 30 min and then incubated with primary OCN antibody (Abcam) at 4 °C overnight. Then, sections were incubated with the corresponding secondary antibodies for 1 h. DAPI was used to stain nuclei, and coverslips were mounted. Images were captured with a LSM 5 EXCITER confocal imaging system (Carl Zeiss, Jena, Germany).

Statistical analysis

SPSS software (version 16.0; SPSS, Inc., Chicago, IL, USA) was used for statistical analyses. All data are expressed as the mean ± SD of at least three independent experiments. Student’s t test was used to analyze differences between groups. One-way analysis of variance (ANOVA) was employed for multiple group testing. A two-tailed P < 0.05 was considered statistically significant.

Results

CircRNA CDR1as is upregulated and miR-7 is downregulated during osteoblast differentiation of PDLSCs

PDLSCs were cultured in osteogenic medium. qRT-PCR results showed that the expression levels of the osteogenic markers ALP, BMP2, RUNX2, and OCN were markedly increased during osteogenesis (Fig. 1b), indicating successful osteogenic induction of PDLSCs. We detected dynamic expression profiles of CDR1as and miR-7 in PDLSCs during osteogenic induction. The expression of CDR1as was gradually upregulated during osteogenic induction. After induction for 14 days, its expression was more than 10-fold higher than the initial level (Fig. 1a). The expression of miR-7 was markedly downregulated in PDLSCs during osteogenic differentiation, and showed an inverse trend to that of CDR1as (Fig. 1a).

Fig. 1
figure 1

Expression patterns of CDR1as and miR-7 during osteogenic differentiation of periodontal ligament stem cells (PDLSCs). a Relative expression of CDR1as as determined by qRT-PCR analysis during osteoblast differentiation of PDLSCs at days 0, 3, 7, and 14 (left panel). RNA expression at indicated time points was normalized to day 0. GAPDH was used as an internal control. Relative expression of miR-7 at the indicated time points (right panel). U6 was used as an internal control. b Relative mRNA expression of ALP, BMP2, RUNX2, and OCN at the indicated time points. GAPDH was used as an internal control, relative to day 0 groups. *p < 0.05, **p < 0.001. GAPDH glyceraldehyde 3-phosphate dehydrogenase, ALP alkaline phosphatase, BMP2 bone morphogenetic protein 2, RUNX2 Runt-related transcription factor 2, OCN osteocalcin

Knockdown of CDR1as inhibits osteogenic differentiation of PDLSCs

Transfection was conducted to knock down the expression of CDR1as. qRT-PCR analysis confirmed that the expression of CDR1as was decreased by approximately 90 and 70% in the two CDR1as knockdown groups (Fig. 2b). The staining and activity of ALP were decreased in the CDR1as knockdown groups after induction for 7 days (Fig. 2a). Following osteogenic induction for 14 days, the intensity of ARS staining was significantly decreased in the CDR1as knockdown groups, indicating that matrix mineralization was decreased (Fig. 2a). In addition, when cultured with GM, the intensities of ALP and ARS staining were also slightly decreased in the CDR1as knockdown groups (Fig. 2a). The results of qRT-PCR revealed that the mRNA expression of ALP and RUNX2 decreased in the CDR1as knockdown groups (Fig. 2b). Additionally, the Western blotting results showed that protein expression of RUNX2 was also downregulated when CDR1as was knocked down in PDLSCs (Fig. 2c).

Fig. 2
figure 2

Knockdown of CDR1as inhibited osteogenic differentiation of PDLSCs. a Images of alkaline phosphatase (ALP) staining in the small interfering RNA negative control (si-NC), small interfering RNAs targeting CDR1as (si-CDR1as)-1, and si-CDR1as-2 groups (top row). Cells were cultured in growth medium (GM) or osteogenic medium (OM) for 7 days. Images of alizarin red S (ARS) staining for mineralized matrix in the si-NC, si-CDR1as-1, and si-CDR1as-2 groups (bottom row). Cells were cultured in GM or OM for 14 days. Histograms show ALP activity (left panel) and quantification of ARS staining (right panel) by spectrophotometry (normalized to the si-NC groups). b The efficiency of transient transduction of si-CDR1as-1 and si-CDR1as-2 (left panel) and relative mRNA expression of ALP (middle panel) and RUNX2 (right panel) measured by qRT-PCR at day 3 of osteogenic induction. GAPDH was used for normalization relative to si-NC groups. c Western blot analysis of protein expression of RUNX2 and the internal control β-ACTIN at day 3 of osteogenic induction. Histograms show the quantification of band intensities. β-ACTIN was used for normalization relative to si-NC groups. *p < 0.05, **p < 0.001. RUNX2 Runt-related transcription factor 2

miR-7 inhibits osteoblast differentiation of PDLSCs

To determine whether miR-7 was involved in the regulation of osteoblast differentiation in PDLSCs, transfection was conducted to overexpress or knock down miR-7. qRT-PCR showed that miR-7 was increased by approximately twofold in the miR-7 overexpression group, and reduced by approximately 70% in the miR-7 knockdown group (Fig. 3b). ALP staining and activity were decreased in the miR-7 overexpression group and increased in the miR-7 knockdown group after 7 days of induction (Fig. 3a). A similar trend in ARS staining was detected after induction for 14 days (Fig. 3a). qRT-PCR analysis confirmed that transfection with miR-7 mimic decreased the mRNA expression of ALP and RUNX2, whereas knockdown of miR-7 increased the gene expression of these osteogenic markers (Fig. 3b). The protein level of RUNX2 was also increased, as revealed by Western blot analysis (Fig. 3c).

Fig. 3
figure 3

miR-7 inhibited osteogenic differentiation of PDLSCs. a Images of alkaline phosphatase (ALP) staining in the miR-NC, miR-7 mimics, and miR-7 inhibitor groups (top row). Cells were cultured in growth medium (GM) or osteogenic medium (OM) for 7 days. Images of alizarin red S (ARS) staining for mineralized matrix in the miR-NC, miR-7 mimics, and miR-7 inhibitor groups (bottom row). Cells were cultured in GM or OM for 14 days. Histograms show ALP activity (left panel) and quantification of ARS staining (right panel) by spectrophotometry (normalized to the microRNA normal control (miR-NC) groups). b The efficiency of transient transduction of miR-7 mimics and miR-7 inhibitor (left panel), and relative mRNA expression of ALP (middle panel) and RUNX2 (right panel) measured by qRT-PCR at day 3 of osteogenic induction. GAPDH was used for normalization relative to miR-NC groups. c Western blot analysis of protein expression of RUNX2 and the internal control β-ACTIN at day 3 of osteogenic induction. The histogram shows the quantification of band intensities. β-ACTIN was used for normalization relative to miR-NC groups. *p < 0.05, **p < 0.001. RUNX2 Runt-related transcription factor 2

GDF5 is regulated by CDR1as/miR-7

To uncover the downstream molecular mechanism of CDR1as/miR-7 regulating the osteoblast lineage differentiation of PDLSCs, we used a target prediction algorithm (RNA22 software) to search for potential targets of miR-7. Notably, we found that the 3′ untranslated region (UTR) of the osteogenic-related gene GDF5 contained miR-7 binding sites (Fig. 4a). Next, two types of luciferase reporters for GDF5 were constructed. The wild-type (WT) reporter contained a wild 3’-UTR of GDF5. The mutant-type (MU) reporter contained a mutant 3’-UTR designed with mutated sequences of the miR-7 binding site. The results indicated that overexpression of miR-7 markedly decreased the luciferase activity in the WT group, whereas this phenomenon was not observed in the MU group (Fig. 4b).

Fig. 4
figure 4

GDF5 was regulated by CDR1as/miR-7. a The seed sequences of miR-7 match the 3’ untranslated region (UTR) of growth differentiation factor 5 (GDF5). b The 293 T cells were cotransfected with 100 nM microRNA normal control (miR-NC), miR-7 mimics, or miR-7 inhibitor and the luciferase constructs carrying 40 ng GDF5 luciferase reporter plasmids. Luciferase assays were performed 24 h after transfection. c Quantification of mRNA expression of GDF5 measured by qRT-PCR in miR-7 inhibitor, miR-7 mimic, small interfering RNAs targeting CDR1as (si-CDR1as), si-CDR1as/miR-7 inhibitor, and si-CDR1as/miR-7 mimic transfected PDLSCs relative to the control groups. GAPDH mRNA levels were used for normalization. d Western blot analysis of protein expression of GDF5 and the internal control β-ACTIN in miR-7 inhibitor, miR-7 mimic, si-CDR1as, si-CDR1as/miR-7 inhibitor, si-CDR1as/miR-7 mimic, and normal control (NC) groups. The histogram shows the quantification of band intensities. *p < 0.05, **p < 0.001

To validate whether GDF5 acts as a downstream target of CDR1as/miR-7, we transfected PDLSCs with miR-7 mimic, miR-7 inhibitor, si-CDR1as, and the corresponding mi-NC and si-NC. qRT-PCR and Western blotting were used to detect the mRNA and protein expression levels of GDF5, respectively. Knockdown of CDR1as significantly repressed both the mRNA and protein levels of GDF5. Consistently, GDF5 was significantly downregulated in cells transfected with miR-7 mimic, whereas it was significantly upregulated in cells transfected with miR-7 inhibitor (Fig. 4c, d). Furthermore, we cotransfected si-CDR1as and miR-7 mimic or inhibitor into PDLSCs and detected the expression of GDF5. Compared with the si-CDR1as group, the mRNA and protein levels of GDF5 were decreased in cells cotransfected with si-CDR1as and miR-7 mimic, whereas GDF5 levels were increased in cells cotransfected with si-CDR1as and miR-7 inhibitor (Fig. 4c, d). These results indicated that CDR1as acts upstream of miR-7 and inhibits the effect of miR-7 in targeting GDF5.

Downregulation of GDF5 inhibits osteogenic differentiation of PDLSCs and partially blocks the osteogenesis-promoting function of CDR1as

To verify whether GDF5 participated in regulating the osteogenic differentiation of PDLSCs, GDF5 was transiently knocked down. qRT-PCR showed that the expression of GDF5 was approximately 70–80% lower in the GDF5 knockdown groups compared with WT (Fig. 5b). Following osteogenic induction for 7 and 14 days, ALP and ARS staining were markedly decreased in the GDF5 knockdown groups (Fig. 5a). Consistently, qRT-PCR showed that the osteogenic markers ALP and RUNX2 were also decreased in the GDF5 knockdown groups (Fig. 5b).

Fig. 5
figure 5

Knockdown of GDF5 inhibited osteogenic differentiation of PDLSCs and partially blocked the pro-osteogenic function of miR-7 inhibitor. a Images of alkaline phosphatase (ALP) staining (top row) in the small interfering RNA normal control (si-NC) and three small interfering RNAs targeting growth differentiation factor 5 (si-GDF5) groups. Cells were cultured in growth medium (GM) or osteogenic medium (OM) for 7 days. Images of alizarin red S (ARS) staining (bottom row) for mineralized matrix in the si-NC and three si-GDF5 groups. Cells were cultured in GM or OM for 14 days. Histograms show ALP activity and quantification of ARS staining by spectrophotometry (normalized to the si-NC groups). b The efficiency of transient transduction of si-GDF5 (left panel) and relative mRNA expression of ALP (middle panel) and RUNX2 (right panel) measured by qRT-PCR at day 3 of osteogenic induction. GAPDH was used for normalization relative to si-NC groups. c Images of ALP staining (top row) in the normal control (NC), miR-7 inhibitor, si-GDF5, and miR-7 inhibitor + si-GDF5 groups. Cells were cultured in GM or OM for 7 days. Images of ARS staining (bottom row) for mineralized matrix in the NC, miR-7 inhibitor, si-GDF5, and miR-7 inhibitor + si-GDF5 groups. Cells were cultured in GM or OM for 14 days. Histograms show ALP activity and quantification of ARS staining by spectrophotometry normalized to the si-NC groups. *p < 0.05, **p < 0.001. RUNX2 Runt-related transcription factor 2

To further confirm whether CDR1as/miR-7 mediate osteoblast differentiation of PDLSCs by altering GDF5 expression, we cotransfected PDLSCs with miR-7 inhibitor and GDF5 siRNA. The cells were then induced to osteoblast lineage using OM. ALP staining and activity were increased in the miR-7 inhibitor group after 7 days of induction. Furthermore, ALP staining and activity were decreased in cells cotransfected with GDF5 siRNA compared with the miR-7 inhibitor group (Fig. 5c). Cotransfection with GDF5 siRNA partially blocked the promotion of matrix mineralization in miR-7 inhibitor-treated PDLSCs (Fig. 5c).

CDR1as/miR-7/GDF5 induces osteogenic differentiation of PDLSCs partly via enhanced Smad1/5/8 and p38 MAPK phosphorylation

GDF5, also known as BMP14, is a member of the transforming growth factor (TGF)-β superfamily. According to previous reports, the MAPK and TGF-β/Smad signaling pathways are involved in the osteogenic process mediated by BMP [28,29,30]. To investigate whether CDR1as/GDF5 induces osteogenic differentiation of PDLSCs via accelerated Smad1/5/8 phosphorylation, Western blot analysis was used to measure pSmad1/5/8 expression in si-CDR1as-, miR-7 mimic-, miR-7 inhibitor-, and si-GDF5-treated PDLSCs. The results indicated that pSmad1/5/8 was markedly decreased in GDF5 knockdown, CDR1as knockdown, and miR-7 mimic groups, while increased in the miR-7 inhibitor group (Fig. 6a, b, and c, respectively). Furthermore, knockdown of GDF5 or CDR1as also inhibited the phosphorylation of p38 MAPK (Fig. 6d and e, respectively). Consistently, p-p38 MAPK was significantly decreased in the miR-7 overexpression group, while it was increased in the miR-7 inhibitor group (Fig. 6f).

Fig. 6
figure 6

CDR1as/miR-7/GDF5 affected the phosphorylation of Smad1/5/8 and p38 MAPK in PDLSCs. Western blot analyses (left) and the relative quantification (right) of phosphorylated Smad1/5/8 (pSmad1/5/8), growth differentiation factor 5 (GDF5), and the internal control β-ACTIN in si-GDF5 (a), si-CDR1as (b), and miR-7 mimic- or inhibitor-treated PDLSCs (c) are shown. Western blot analysis (left) and the relative quantification (right) of phosphorylated p38 mitogen-activated protein kinase (p-p38), total p38 mitogen-activated protein kinase (p38), and the internal control β-ACTIN in si-GDF5 (d), si-CDR1as (e), and miR-7 mimic- or inhibitor-treated PDLSCs (f) are shown. *p < 0.05, **p < 0.001. NC normal control, si small interfering RNA

CDR1as knockdown inhibits bone formation in vivo

To further verify the role of CDR1as in osteoblast differentiation of PDLSCs, we conducted animal experiments in vivo. CDR1as siRNA- and negative control siRNA-treated PDLSCs were loaded on scaffold materials. The seeded scaffolds were then gently implanted in the calvarial defect area of nude mice (five mice per group). After 8 weeks, the mice were sacrificed and skull specimens were harvested for further studies. Three-dimensional reconstructed micro-CT images were used to visualize the repair of bone defects. Less bone formation and a larger defect area were observed in the CDR1as knockdown group compared with the control group (Fig. 7a). The BMD and BV/TV were significantly decreased in the CDR1as knockdown group (Fig. 7a). The histological observations were consistent with the findings from micro-CT analysis. H&E staining indicated new bone formation around the bone defect areas in both groups. However, the amount of new bone formation in the CDR1as knockdown group was significantly lower than that in the control group. In the control group, the newly formed bone tissue fused with the edge of the bone defect, but in the CDR1as knockdown group the ends of the bone were blunted and new bone was scarcely observed (Fig. 7b). Moreover, immunofluorescence staining showed that the bone tissue was positive for OCN staining, whereas the staining intensity in the CDR1as knockdown group was less than that in the control group (Fig. 7c).

Fig. 7
figure 7

Knockdown of CDR1as inhibited bone formation in vivo. a Reconstructed three-dimensional micro-CT images of the calvarial defect area of nude mice contained normal control (NC) and small interfering RNA (si)-CDR1as treated PDLSCs. The histograms show the bone mineral density (BMD; left panel) and the ratio of new bone volume to existing tissue volume (BV/TV; right panel) of the two groups. b Hematoxylin and eosin (H&E) staining in NC and si-CDR1as groups. Bone formation (B) around the scaffold (S) and the original cranial bones (C) were identified. Scale bar = 100 μm. c Confocal microscopy of osteocalcin (OCN) with DAPI counterstaining in NC and si-CDR1as groups in the calvarial defect area. Scale bars = 50 μm. d Schematic illustration showed the signaling of CDR1as/miR-7/GDF5/Smad and p38 MAPK in PDLSCs. CDR1as relieved the negative regulation of miR-7 on growth differentiation factor 5 (GDF5) by acting as an miR-7 sponge and subsequently activated the pSmad1/5/8 and p38 MAPK (p38) pathway

Discussion

In this study, circRNA CDR1as was demonstrated to promote the osteogenic differentiation of PDLSCs. PDLSCs are an excellent cell source with potential for multilineage differentiation, and are considered ideal candidates for use in periodontal tissue and bone regenerative applications [3, 31]. Therefore, it is essential to understand the mechanisms underlying their multilineage differentiation potential. Although numerous studies have investigated mRNAs, the regulatory role of noncoding RNAs in the osteogenesis process of PDLSCs is relatively unknown. CircRNAs, a novel type of noncoding RNAs, were initially identified as a product formed by aberrant gene splicing and were long considered nonfunctional transcriptional noise [32, 33]. However, studies have increasingly revealed that circRNAs participate in a variety of cellular activities and are involved in many human diseases [11, 34,35,36]. CDR1as is a well-identified circRNA [15]. The importance of CDR1as in embryonic development and many disease processes is still under investigation [11, 17,18,19,20], and its role in the osteogenic differentiation process of stem cells has not been defined. Using RNA sequencing, we previously revealed that the expression pattern of CDR1as is significantly altered during osteogenic lineage induction of PDLSCs [14]. In this study, ALP and ARS staining indicated that knockdown of CDR1as inhibits osteoblast differentiation of PDLSCs in vitro. At the same time, qRT-PCR and Western blot analysis confirmed this at the RNA and protein levels. In addition, a critical-sized mouse calvarial defect model was used to further confirm the role of CDR1as in the osteogenesis process in vivo. The results of micro-CT scanning and histological examination were consistent with the in-vitro experiments. These results indicate that CDR1as has a critical function in the osteogenic process. Furthermore, previous studies have shown that CDR1as impacted biological functions by acting as an miRNA sponge or inhibitor of miR-7 [15]. Moreover, miR-7 has been reported to be involved in osteogenesis of human skeletal muscle-derived progenitor cells [37]. We also observed that miR-7 knockdown promoted osteogenic differentiation, whereas osteogenic differentiation was inhibited in the presence of miR-7 mimics. These findings suggest that CDR1as promotes the osteogenic differentiation of PDLSCs by negatively regulating miR-7.

The CDR1as/miR-7 pathway promotes the osteogenic differentiation of PDLSCs partially through the regulation of GDF5. Bioinformatics analysis showed that the 3’-UTR region of GDF5 contained miR-7 binding sites. Overexpression of miR-7 significantly reduced the luciferase reporter activity compared with the WT group. When the sequence of the binding site was mutated, this inhibitory effect was effectively reversed. Moreover, CDR1as knockdown and overexpression of miR-7 markedly decreased the mRNA and protein expression levels of GDF5, whereas knockdown of miR-7 increased GDF5 mRNA and protein expression. Furthermore, cotransfection with miR-7 mimic and si-CDR1as led to a more significant decrease of GDF5 compared with the si-CDR1as alone, whereas cotransfection with miR-7 inhibitor and si-CDR1as partially blocked the effects of si-CDR1as. These results confirmed that CDR1as acts upstream of miR-7 and inhibits the effect of miR-7 in targeting GDF5. GDF5, also known as BMP14, is a member of the TGF-β superfamily and has been demonstrated to be closely involved in the osteogenesis process [38, 39]. Previous studies have shown that GDF5 plays critical roles in skeletal development and encourages osteogenic differentiation [40, 41]. We also found that GDF5 knockdown significantly inhibited the osteogenic differentiation of PDLSCs, and the effects were similar to those observed with CDR1as siRNA or miR-7 mimic. Furthermore, GDF5 knockdown partially blocked the osteogenesis-promoting effects of an miR-7 inhibitor. These finding indicate that the CDR1as/miR-7 pathway promotes the osteogenic differentiation of PDLSCs partially by regulating GDF5. However, knockdown of GDF5 did not completely block the osteogenic effects of miR-7 knockdown. These results suggest that other downstream signaling pathways may participate in the CDR1as regulation of osteogenic differentiation of PDLSCs.

CDR1as/GDF5 induced the osteogenic differentiation of PDLSCs partially via enhanced Smad1/5/8 and p38 MAPK phosphorylation. Previous studies have shown that GDF5 can bind to BMP type I and II receptors located on the cell membrane, transmit the signal into the cytoplasm, and then activate the phosphorylation of Smad1/5/8. pSmad1/5/8 further forms complexes with Smad4 and translocates into the nucleus and then binds to target genes to regulate transcription [42, 43]. In this study, we observed that knockdown of both GDF5 and CDR1as decreased the protein expression of pSmad1/5/8. Many studies have reported a relationship between pSmad1/5/8 and the osteogenic differentiation process [44, 45]. Besides, GDF5 can also mediates activation of the MAPK signaling transduction pathway [46]. The p38 MAPK pathway has been reported to participate in the osteogenic differentiation of mesenchymal stem cells [29, 47, 48]. Here, we also observed that knockdown of GDF5 or CDR1as decreased the protein expression of p-p38 MAPK, and p-p38 MAPK was significantly decreased in the miR-7 overexpression group while it was increased in the miR-7 inhibitor group. Therefore, we hypothesize that the circRNA CDR1as/GDF5-related promotion of osteogenic differentiation of PDLSCs is modulated through a Smad-dependent and p38 MAPK signaling pathway (Fig. 7d).

Conclusions

In conclusion, our results suggest that circRNA CDR1as inhibits miR-7 expression, relieving the negative regulation of GDF5 by miR-7, and further activates the pSmad1/5/8 and p38 MAPK pathway, thus promoting osteogenesis in PDLSCs.

Abbreviations

ALP:

Alkaline phosphatase

ARS:

Alizarin red S

BMD:

Bone mineral density

BMP2:

Bone morphogenetic protein 2

BV/TV:

Ratio of new bone volume to existing tissue volume

CDR1as:

Antisense to the cerebellar degeneration-related protein 1 transcript

circRNA:

Circular RNA

ciRS-7:

Circular RNA sponge for miR-7

CT:

Computed tomography

EDTA:

Ethylenediaminetetraacetic acid

GAPDH:

Glyceraldehyde 3-phosphate dehydrogenase

GDF:

Growth differentiation factor

GM:

Growth medium

H&E:

Hematoxylin and eosin

MEM:

Modified Eagle’s medium

miRNA:

MicroRNA

miR-NC:

MicroRNA control

MU:

Mutant-type

OCN:

Osteocalcin

OM:

Osteogenic induction medium

PBS:

Phosphate-buffered saline

PDLSC:

Periodontal ligament stem cell

PLGA:

Polylactic-co-glycolic acid

p38 MAPK:

p38 mitogen-activated protein kinase

p-p38 MAPK:

Phosphorylated p38 mitogen-activated protein kinase

pSmad1/5/8:

Phosphorylated Smad1/5/8

PVDF:

Polyvinylidene difluoride

qRT-PCR:

Quantitative reverse-transcription polymerase chain reaction

RIPA:

Radio immunoprecipitation assay

RUNX2:

Runt-related transcription 2

si-CDR1as:

Small interfering RNAs targeting CDR1as

si-GDF5:

Small interfering RNAs targeting GDF5

si-NC:

Small interfering RNA control

siRNA:

Small interfering RNA

TGF:

Transforming growth factor

UTR:

Untranslated region

WT:

Wild-type

References

  1. Beertsen W, McCulloch CA, Sodek J. The periodontal ligament: a unique, multifunctional connective tissue. Periodontol 2000. 1997;13:20–40.

    Article  PubMed  CAS  Google Scholar 

  2. Shimono M, Ishikawa T, Ishikawa H, Matsuzaki H, Hashimoto S, Muramatsu T, et al. Regulatory mechanisms of periodontal regeneration. Microsc Res Tech. 2003;60:491–502.

    Article  PubMed  CAS  Google Scholar 

  3. Seo BM, Miura M, Gronthos S, Bartold PM, Batouli S, Brahim J, et al. Investigation of multipotent postnatal stem cells from human periodontal ligament. Lancet. 2004;364:149–55. https://doi.org/10.1016/S0140-6736(04)16627-0.

    Article  PubMed  CAS  Google Scholar 

  4. Park JY, Jeon SH, Choung PH. Efficacy of periodontal stem cell transplantation in the treatment of advanced periodontitis. Cell Transplant. 2011;20:271–85. https://doi.org/10.3727/096368910X519292.

    Article  PubMed  Google Scholar 

  5. Menicanin D, Mrozik KM, Wada N, Marino V, Shi S, Bartold PM, et al. Periodontal-ligament-derived stem cells exhibit the capacity for long-term survival, self-renewal, and regeneration of multiple tissue types in vivo. Stem Cells Dev. 2014;23:1001–11. https://doi.org/10.1089/scd.2013.0490.

    Article  PubMed  CAS  Google Scholar 

  6. Yu BH, Zhou Q, Wang ZL. Periodontal ligament versus bone marrow mesenchymal stem cells in combination with bio-Oss scaffolds for ectopic and in situ bone formation: a comparative study in the rat. J Biomater Appl. 2014;29:243–53. https://doi.org/10.1177/0885328214521846.

    Article  PubMed  CAS  Google Scholar 

  7. Tsumanuma Y, Iwata T, Washio K, Yoshida T, Yamada A, Takagi R, et al. Comparison of different tissue-derived stem cell sheets for periodontal regeneration in a canine 1-wall defect model. Biomaterials. 2011;32:5819–25. https://doi.org/10.1016/j.biomaterials.2011.04.071.

    Article  PubMed  CAS  Google Scholar 

  8. Sanger HL, Klotz G, Riesner D, Gross HJ, Kleinschmidt AK. Viroids are single-stranded covalently closed circular RNA molecules existing as highly base-paired rod-like structures. Proc Natl Acad Sci U S A. 1976;73:3852–6.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  9. Danan M, Schwartz S, Edelheit S, Sorek R. Transcriptome-wide discovery of circular RNAs in archaea. Nucleic Acids Res. 2012;40:3131–42. https://doi.org/10.1093/nar/gkr1009.

    Article  PubMed  CAS  Google Scholar 

  10. Chen LL. The biogenesis and emerging roles of circular RNAs. Nat Rev Mol Cell Biol. 2016;17:205–11. https://doi.org/10.1038/nrm.2015.32.

    Article  PubMed  CAS  Google Scholar 

  11. 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. https://doi.org/10.1038/nature11928.

    Article  PubMed  CAS  Google Scholar 

  12. Hansen TB, Kjems J, Damgaard CK. Circular RNA and miR-7 in cancer. Cancer Res. 2013;73:5609–12. https://doi.org/10.1158/0008-5472.CAN-13-1568.

    Article  PubMed  CAS  Google Scholar 

  13. Guarnerio J, Bezzi M, Jeong JC, Paffenholz SV, Berry K, Naldini MM, et al. Oncogenic role of fusion-circRNAs derived from cancer-associated chromosomal translocations. Cell. 2016;166:1055–6. https://doi.org/10.1016/j.cell.2016.03.020.

    Article  PubMed  CAS  Google Scholar 

  14. Zheng Y, Li X, Huang Y, Jia L, Li W. The circular RNA landscape of periodontal ligament stem cells during osteogenesis. J Periodontol. 2017;88:906–14. https://doi.org/10.1902/jop.2017.170078.

    Article  PubMed  Google Scholar 

  15. Hansen TB, Jensen TI, Clausen BH, Bramsen JB, Finsen B, Damgaard CK, et al. Natural RNA circles function as efficient microRNA sponges. Nature. 2013;495:384–8. https://doi.org/10.1038/nature11993.

    Article  PubMed  CAS  Google Scholar 

  16. Hansen TB, Wiklund ED, Bramsen JB, Villadsen SB, Statham AL, Clark SJ, et al. miRNA-dependent gene silencing involving Ago2-mediated cleavage of a circular antisense RNA. EMBO J. 2011;30:4414–22. https://doi.org/10.1038/emboj.2011.359.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  17. Xu H, Guo S, Li W, Yu P. The circular RNA Cdr1as, via miR-7 and its targets, regulates insulin transcription and secretion in islet cells. Sci Rep. 2015;5:12453. https://doi.org/10.1038/srep12453.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Tang W, Ji M, He G, Yang L, Niu Z, Jian M, et al. Silencing CDR1as inhibits colorectal cancer progression through regulating microRNA-7. Onco Targets Ther. 2017;10:2045–56. https://doi.org/10.2147/OTT.S131597.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  19. Yang X, Xiong Q, Wu Y, Li S, Ge F. Quantitative proteomics reveals the regulatory networks of circular RNA CDR1as in hepatocellular carcinoma cells. J Proteome Res. 2017;16:3891–902. https://doi.org/10.1021/acs.jproteome.7b00519.

    Article  PubMed  CAS  Google Scholar 

  20. Geng HH, Li R, Su YM, Xiao J, Pan M, Cai XX, et al. The circular RNA Cdr1as promotes myocardial infarction by mediating the regulation of miR-7a on its target genes expression. PLoS One. 2016;11:e0151753. https://doi.org/10.1371/journal.pone.0151753.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  21. Xu L, Zhang M, Zheng X, Yi P, Lan C, Xu M. The circular RNA ciRS-7 (Cdr1as) acts as a risk factor of hepatic microvascular invasion in hepatocellular carcinoma. J Cancer Res Clin Oncol. 2017;143:17–27. https://doi.org/10.1007/s00432-016-2256-7.

    Article  PubMed  CAS  Google Scholar 

  22. Huang Y, Zheng Y, Jia L, Li W. Long noncoding RNA H19 promotes osteoblast differentiation via TGF-beta1/Smad3/HDAC signaling pathway by deriving miR-675. Stem Cells. 2015;33:3481–92. https://doi.org/10.1002/stem.2225.

    Article  PubMed  CAS  Google Scholar 

  23. Huang Y, Jin C, Zheng Y, Li X, Zhang S, Zhang Y, et al. Knockdown of lncRNA MIR31HG inhibits adipocyte differentiation of human adipose-derived stem cells via histone modification of FABP4. Sci Rep. 2017;7:8080. https://doi.org/10.1038/s41598-017-08131-6.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Huang Y, Zheng Y, Jin C, Li X, Jia L, Li W. Long non-coding RNA H19 inhibits adipocyte differentiation of bone marrow mesenchymal stem cells through epigenetic modulation of histone deacetylases. Sci Rep. 2016;6:28897. https://doi.org/10.1038/srep28897.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Krebsbach PH, Mankani MH, Satomura K, Kuznetsov SA, Robey PG. Repair of craniotomy defects using bone marrow stromal cells. Transplantation. 1998;66:1272–8.

    Article  PubMed  CAS  Google Scholar 

  26. Cooper GM, Mooney MP, Gosain AK, Campbell PG, Losee JE, Huard J. Testing the critical size in calvarial bone defects: revisiting the concept of a critical-size defect. Plast Reconstr Surg. 2010;125:1685–92. https://doi.org/10.1097/PRS.0b013e3181cb63a3.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Herberg S, Susin C, Pelaez M, Howie RN, Moreno de Freitas R, Lee J, et al. Low-dose bone morphogenetic protein-2/stromal cell-derived factor-1beta cotherapy induces bone regeneration in critical-size rat calvarial defects. Tissue Eng Part A. 2014;20:1444–53. https://doi.org/10.1089/ten. TEA.2013.0442

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Wu Y, Xia L, Zhou Y, Xu Y, Jiang X. Icariin induces osteogenic differentiation of bone mesenchymal stem cells in a MAPK-dependent manner. Cell Prolif. 2015;48:375–84. https://doi.org/10.1111/cpr.12185.

    Article  PubMed  CAS  Google Scholar 

  29. Kim HY, Park SY, Choung SY. Enhancing effects of myricetin on the osteogenic differentiation of human periodontal ligament stem cells via BMP-2/Smad and ERK/JNK/p38 mitogen-activated protein kinase signaling pathway. Eur J Pharmacol. 2018;834:84–91. https://doi.org/10.1016/j.ejphar.2018.07.012.

    Article  PubMed  CAS  Google Scholar 

  30. Zhu JH, Liao YP, Li FS, Hu Y, Li Q, Ma Y, et al. Wnt11 promotes BMP9-induced osteogenic differentiation through BMPs/Smads and p38 MAPK in mesenchymal stem cells. J Cell Biochem. 2018; https://doi.org/10.1002/jcb.27262.

  31. Tassi SA, Sergio NZ, Misawa MYO, Villar CC. Efficacy of stem cells on periodontal regeneration: systematic review of pre-clinical studies. J Periodontal Res. 2017;52:793–812. https://doi.org/10.1111/jre.12455.

    Article  PubMed  CAS  Google Scholar 

  32. Hsu MT, Coca-Prados M. Electron microscopic evidence for the circular form of RNA in the cytoplasm of eukaryotic cells. Nature. 1979;280:339–40.

    Article  PubMed  CAS  Google Scholar 

  33. Cocquerelle C, Daubersies P, Majerus MA, Kerckaert JP, Bailleul B. Splicing with inverted order of exons occurs proximal to large introns. EMBO J. 1992;11:1095–8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  34. Jeck WR, Sorrentino JA, Wang K, Slevin MK, Burd CE, Liu J, et al. Circular RNAs are abundant, conserved, and associated with ALU repeats. RNA. 2013;19:141–57. https://doi.org/10.1261/rna.035667.112.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  35. Zhang XO, Wang HB, Zhang Y, Lu X, Chen LL, Yang L. Complementary sequence-mediated exon circularization. Cell. 2014;159:134–47. https://doi.org/10.1016/j.cell.2014.09.001.

    Article  PubMed  CAS  Google Scholar 

  36. Zhao ZJ, Shen J. Circular RNA participates in the carcinogenesis and the malignant behavior of cancer. RNA Biol. 2017;14:514–21. https://doi.org/10.1080/15476286.2015.1122162.

    Article  PubMed  Google Scholar 

  37. Oishi T, Uezumi A, Kanaji A, Yamamoto N, Yamaguchi A, Yamada H, et al. Osteogenic differentiation capacity of human skeletal muscle-derived progenitor cells. PLoS One. 2013;8:e56641. https://doi.org/10.1371/journal.pone.0056641.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  38. Baur ST, Mai JJ, Dymecki SM. Combinatorial signaling through BMP receptor IB and GDF5: shaping of the distal mouse limb and the genetics of distal limb diversity. Development. 2000;127:605–19.

    PubMed  CAS  Google Scholar 

  39. Nickel J, Kotzsch A, Sebald W, Mueller TD. A single residue of GDF-5 defines binding specificity to BMP receptor IB. J Mol Biol. 2005;349:933–47. https://doi.org/10.1016/j.jmb.2005.04.015.

    Article  PubMed  CAS  Google Scholar 

  40. Le Maitre CL, Freemont AJ, Hoyland JA. Expression of cartilage-derived morphogenetic protein in human intervertebral discs and its effect on matrix synthesis in degenerate human nucleus pulposus cells. Arthritis Res Ther. 2009;11:R137. https://doi.org/10.1186/ar2808.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  41. Qu X, Chen Z, Fan D, Sun C, Zeng Y. MiR-132-3p regulates the osteogenic differentiation of thoracic ligamentum flavum cells by inhibiting multiple osteogenesis-related genes. Int J Mol Sci. 2016;17:8. https://doi.org/10.3390/ijms17081370.

    Article  CAS  Google Scholar 

  42. Sieber C, Ploger F, Schwappacher R, Bechtold R, Hanke M, Kawai S, et al. Monomeric and dimeric GDF-5 show equal type I receptor binding and oligomerization capability and have the same biological activity. Biol Chem. 2006;387:451–60. https://doi.org/10.1515/BC.2006.060.

    Article  PubMed  CAS  Google Scholar 

  43. Heinecke K, Seher A, Schmitz W, Mueller TD, Sebald W, Nickel J. Receptor oligomerization and beyond: a case study in bone morphogenetic proteins. BMC Biol. 2009;7:59. https://doi.org/10.1186/1741-7007-7-59.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  44. Prall WC, Haasters F, Heggebo J, Polzer H, Schwarz C, Gassner C, et al. Mesenchymal stem cells from osteoporotic patients feature impaired signal transduction but sustained osteoinduction in response to BMP-2 stimulation. Biochem Biophys Res Commun. 2013;440:617–22. https://doi.org/10.1016/j.bbrc.2013.09.114.

    Article  PubMed  CAS  Google Scholar 

  45. Kim DY, Kim GW, Chung SH. Nectandrin a enhances the BMP-induced osteoblastic differentiation and mineralization by activation of p38 MAPK-Smad signaling pathway. Korean J Physiol Pharmacol. 2013;17:447–53. https://doi.org/10.4196/kjpp.2013.17.5.447.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  46. Daniels J, Binch AA, Le Maitre CL. Inhibiting IL-1 signaling pathways to inhibit catabolic processes in disc degeneration. J Orthop Res. 2017;35:74–85. https://doi.org/10.1002/jor.23363.

    Article  PubMed  CAS  Google Scholar 

  47. Lu Y, Zhao Q, Liu Y, Zhang L, Li D, Zhu Z, et al. Vibration loading promotes osteogenic differentiation of bone marrow-derived mesenchymal stem cells via p38 MAPK signaling pathway. J Biomech. 2018;71:67–75. https://doi.org/10.1016/j.jbiomech.2018.01.039.

    Article  PubMed  Google Scholar 

  48. Li S, Wang J, Han Y, Li X, Liu C, Lv Z, et al. Carbenoxolone inhibits mechanical stress-induced osteogenic differentiation of mesenchymal stem cells by regulating p38 MAPK phosphorylation. Exp Ther Med. 2018;15:2798–803. https://doi.org/10.3892/etm.2018.5757.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Funding

This study was financially supported by grants from the National Natural Science Foundation of China (nos. 81670957, 81772876, and 81700938), the Beijing Natural Science Foundation (no. 7172239), and sponsored by the Fund for Fostering Young Scholars of Peking University Health Science Center (BMU2018PY022).

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Author information

Authors and Affiliations

Authors

Contributions

XL was responsible for performing the experiments, collecting and analyzing the data, and drafting and revising the manuscript. YuZ was responsible for the study design, analyzing the data, and drafting and revising the manuscript. YaZ and YiZ were responsible for collecting the data. YH was responsible for analyzing the data. LJ was responsible for the study design, data analysis, and manuscript revision. WL was responsible for the study design, administrative support, and financial support. All authors approved the final version of the manuscript. XL and YuZ contributed equally to this work.

Corresponding authors

Correspondence to Lingfei Jia or Weiran Li.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Committee of Peking University School of Stomatology (PKUSSIRB-201837096).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1:

Table S1. The sequences of the RNA oligoribonucleotide used in this study. Table S2. The primers used in this study. (DOCX 22 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, X., Zheng, Y., Zheng, Y. et al. Circular RNA CDR1as regulates osteoblastic differentiation of periodontal ligament stem cells via the miR-7/GDF5/SMAD and p38 MAPK signaling pathway. Stem Cell Res Ther 9, 232 (2018). https://doi.org/10.1186/s13287-018-0976-0

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13287-018-0976-0

Keywords