- Open Access
Circular RNA CDR1as regulates osteoblastic differentiation of periodontal ligament stem cells via the miR-7/GDF5/SMAD and p38 MAPK signaling pathway
Stem Cell Research & Therapyvolume 9, Article number: 232 (2018)
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.
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.
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.
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.
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 . 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 . 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 , 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 . 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 . CDR1as was also found to regulate the transcription and secretion of insulin by blocking miR-7 . 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.
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.
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 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 . 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 . 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 . 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 . 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 . 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 . 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).
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.
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).
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).
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).
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).
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).
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).
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).
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 . 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 . 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 . Moreover, miR-7 has been reported to be involved in osteogenesis of human skeletal muscle-derived progenitor cells . 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 . 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).
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.
Alizarin red S
Bone mineral density
Bone morphogenetic protein 2
Ratio of new bone volume to existing tissue volume
Antisense to the cerebellar degeneration-related protein 1 transcript
Circular RNA sponge for miR-7
Glyceraldehyde 3-phosphate dehydrogenase
Growth differentiation factor
Hematoxylin and eosin
Modified Eagle’s medium
Osteogenic induction medium
Periodontal ligament stem cell
- p38 MAPK:
p38 mitogen-activated protein kinase
- p-p38 MAPK:
Phosphorylated p38 mitogen-activated protein kinase
Quantitative reverse-transcription polymerase chain reaction
Radio immunoprecipitation assay
Runt-related transcription 2
Small interfering RNAs targeting CDR1as
Small interfering RNAs targeting GDF5
Small interfering RNA control
Small interfering RNA
Transforming growth factor
Beertsen W, McCulloch CA, Sodek J. The periodontal ligament: a unique, multifunctional connective tissue. Periodontol 2000. 1997;13:20–40.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Krebsbach PH, Mankani MH, Satomura K, Kuznetsov SA, Robey PG. Repair of craniotomy defects using bone marrow stromal cells. Transplantation. 1998;66:1272–8.
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.
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
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. The sequences of the RNA oligoribonucleotide used in this study. Table S2. The primers used in this study. (DOCX 22 kb)