Phenotype identification of DPSCs
The morphology of primary generation DPSCs were fibroblast- or spindle-like (Figure S1A). To identify the phenotype and qualification of extracted DPSCs, the multipotency, including chondrogenic, adipogenic and osteogenic differentiation, was tested [23]. Flow cytometry results demonstrated that the isolated DPSCs were negative for hematopoietic markers (CD34, CD45) (Figure S1B), but positive for MSC markers (CD29, CD90, CD73 and CD105) (Figure S1C). Tri‐lineage differentiation of DPSCs was firstly confirmed (Figure S1D). Meanwhile, immunofluorescence staining results showed that DPSCs were positive for the MSC surface molecule STRO-1 (Figure S1D). The above results all verified the stem cell characteristics of isolated DPSCs.
Identification of the circular structure
Head-to-tail splicing of circSIPA1L1, and its genome size and sequences were confirmed by Sanger sequencing (Figure 1A). Moreover, divergent and convergent primers were used and it is found that circ-SIPA1L1, but not linear SIPA1L1 was resistant to RNase R digestion (Figure 1B). To exclude the possibility that head-to-tail splicing product of circSIPA1L1 comes from genomic rearrangement or trans-splicing, its cDNA and gDNA of 293T cells either with RNase R or not were detected. Supplemental expression level of reverse splicing or canonical form of SIPA1L1 was shown (Figure 1C). Subsequently, FISH identified that circSIPA1L1 was mainly distributed in the cytoplasm of DPSCs with 18S as an internal control (Figure 1D). We hypothesized that circSIPA1L1 regulates the biological characteristics of DPSCs via the ceRNA mechanism. In summary, circSIPA1L1 was identified as a stable circRNA and deserved further exploration.
CircSIPA1L1 is upregulated and miR-617 is downregulated during osteogenesis of DPSCs
Dynamically expressed circSIPA1L1 and miR-617 during osteogenesis in DPSCs were detected. Three circSIPA1L1 small interfering RNAs (siRNAs) specifically targeting the backsplice junction sequences at different binding sites in circSIPA1L1 were designed (Figure 1E). Small interfering RNA transfection efficiency was detected by RT-PCR. The results showed that si-circSIPA1L1-1 and si-circSIPA1L1-3 could effectively knock down the expression of circSIPA1L1(Figure 1F). Meanwhile, the expression of circSIPA1L1 between NC and circSIPA1L1 group weredetected by RT-PCR, the results showed that circSIPA1L1 could effectively increase the expression of circSIPA1L1(Figure 1G).
CircSIPA1L1 was time-dependently upregulated and miR-617 was downregulated in osteogenic DPSCs. Moreover, mRNA levels of osteogenesis markers ALP, OSX and RUNX2 were remarkably upregulated during the process of osteogenesis (Figure 1H), demonstrating the successful induction of osteogenesis.
CircSIPA1L1 have no effect on DPSCs proliferation
To elucidate the role of circSIPA1L1 in DPSCs proliferation, CCK-8, flow cytometry and EdU assay were conducted. FCM analysis did not show significant differences in the proliferation index (PI = G2M ± S) between the NC group (9.15%) and the circSIPA1L1 group (8.15%, P> 0.05, Figure S2A). Similarly, no significant difference was found in the proliferation index between the si-NC group (4.72%), the si-circ-SIPA1L1-1 group (5.23%), and the si-circ-SIPA1L1-3 group (5.32%, P>0.05, Figure S2A). In addition, the results of the EdU assay showed no significant difference between the NC group and the circSIPA1L1 group (P > 0.05, Figure S2B, C) or between the si-NC, si-circ-SIPA1L1-1 and si-circ-SIPA1L1-3 groups (Figure S2B,D). The CCK-8 assay showed no significant difference in proliferation rates between the NC group and the circSIPA1L1 group (Figure S2E) or between the si-NC, si-circ-SIPA1L1-1 and si-circ-SIPA1L1-3 groups from 0 days to 9 days (P> 0.05) (Figure S2F). Taken together, the data demonstrated that circSIPA1L1 does not affect the proliferation of DPSCs.
CircSIPA1L1 stimulates DPSCs osteogenesis
To further analyze the effect of circSIPA1L1 on osteogenic differentiation of DPSCs, protein and mRNA levels of ALP, OSX and RUNX2 were detected by Western blot and RT-PCR in osteogenic DPSCs. Western blot results showed that protein expression of ALP, OSX and RUNX2 were up-regulated in the overexpression group of circSIPA1L1 (Figure 2A).The results of RT-PCR indicated that circSIPA1L1 overexpression increased ALP, OSX and RUNX2 (Figure 2E).Whereas the expression of protein level was down-regulated when circSIPA1L1 was knocked down in DPSCs (Figure. 2B), and the results of RT-PCR indicated that the level of ALP, OSX and RUNX2 were decreased in circSIPA1L1 knockdown of DPSCs (Figure 2H).After 7 days of osteogenesis, ALP staining showed decreased ALP activity after knockdown of circSIPA1L1 and obviously upregulated by circSIPA1L1 overexpression (Figure 2C). After 14 days of induction, alizarin red staining showed reduced matrix mineralization in DPSCs with circSIPA1L1 knockdown whereas circSIPA1L1 overexpression obtained the opposite effects (Figure 2C, F, I). Identically, positive expressions of ALP and OSX were downregulated by circSIPA1L1 knockdown in DPSCs as immunofluorescence revealed (Figure 4D). These results indicated that circSIPA1L1 stimulated DPSCs osteogenesis.
DPSCs stably down expressing circSIPA1L1 and controls were loaded on Bio-Oss Collagen scaffolds, and implanted in the subcutaneous tissues of nude mice for 8 weeks growth.Both H&E and Masson staining showed less bone-like structures and collagen deposit in DPSCs of the circSIPA1L1-downexpressing group than control group (Figure 2G).
CircSIPA1L1 sponges miR-617
CircRNAs are able to regulate downstream gene expressions and functions by sponging corresponding miRNAs. Transfection efficiency of miR-617 mimics and inhibitor was verified by qRT-PCR (Figure 3A). It is shown that circSIPA1L1 expression was negatively regulated by miR-617 (Figure 3B). To further validate their interaction, FISH analysis was conducted in DPSCs, and the results revealed that miR-617 colocalized with circSIPA1L1 in the cytoplasm (Figure 3C, D).Through analyses on miRDB, miRTarBase, and TargetScan database, a binding site in 3’UTR of miR-617 and circSIPA1L1 was discovered (Figure 3E). Subsequently, dual-luciferase reporter assay was conducted to test the interaction between circSIPA1L1 and miR-617. 293T cells were co-transfected with miR-617 mimics/negative control and wild-type/mutant-type circSIPA1L1, respectively. Overexpression of miR-617 markedly quenched luciferase activity in wild-type circSIPA1L1 compared with controls, verifying the direction interaction between circSIPA1L1 and miR-617 (Figure 3F). These observations indicated that circSIPA1L1 and miR-617 coexisted in the cytoplasm, and circSIPA1L1 acted as a miRNA sponge for miR-617 in DPSCs.
MiR-617 inhibits DPSCs osteogenesis
Next, the potential influence of miR-617 on DPSCs osteogenesis was explored. Western blot revealed that the protein levels of ALP, OSX and RUNX2 increased in the miR-617 knockdown group and the opposite effect was observed in the miR-617 overexpression group (Figure 4A, B).RT-PCR analysis confirmed that the expression of osteogenic related genes ALP, OSX and RUNX2 were significantly lower in the miR-617 overexpressing group than in the control group, while knockdown of miR-617 increased the gene expression of these osteogenic markers (Figure 4C). After 7 days of osteogenesis, ALP staining showed that miR-617 negatively regulated ALP activity in DPSCs (Figure 4D). After 14 days of osteogenesis, alizarin red staining showed that the formation of mineralized nodules in DPSCs was negatively mediated by miR-617 as well (Figure 4E, F). Immunofluorescence staining analysis showed that positive expressions of ALP and OSX were upregulated in DPSCs with miR-617 knockdown, which were downregulated in those overexpressing miR-617 (Figure 4G).In conclusion, the above findings demonstrated that miR-617 was a negative regulator in DPSCs osteogenesis.
MiR-617 directly targets Smad3
Similarly, downstream genes binding miR-617 were predicted using miRDB, miRTarBase, miRWalk and TargetScan algorithms (Figure 5A). A total of 10,461 potential target genes of miR-617 were obtained (see supplement file1). GO and KEGG pathway analysis indicated that these target genes were mainly involved in intracellular activities (Figure 5B, C). Interestingly, Smad3 was a shared gene predicted in miRDB, miRWalk and TargetScan databases. As an intracellular protein, Smad3 induces nuclear transportation of extracellular transforming growth factor β ligands, thereafter activating transcription of downstream genes. Binding sequences in 3’UTR of Smad3 and miR-617 were shown (Figure 5D), and the complementary regions between these different species were also highly conserved.
CircSIPA1L1/miR-617/Smad3 axis is responsible for DPSCs osteogenesis
To test the interaction between miR-617 and Smad3, pciCHECK2-Smad3 and psiCHECK2-mut-Smad3 were constructed. Dual-luciferase reporter assay uncovered decreased luciferase activity after co-transfection of miR-617 mimics and pciCHECK2-Smad3, confirming the direct interaction between miR-617 and Smad3 (Figure 5D). Interestingly, Smad3 level was positively regulated by circSIPA1L1, but negatively regulated by miR-617 (Figure 6A-C). Immunofluorescence assay revealed a similar result (Figure 6D).
MiR-617 reversed regulatory effect of circSIPA1L1 on DPSCs osteogenesis
Rescue experiments were conducted to clarify the involvement of miR-617/Smad3 in circSIPA1L1-mediated osteogenesis. Western blot results showed that downregulated RUNX2, ALP, OSX and Smad3 in osteogenic DPSCs with circSIPA1L1 knockdown were partially reversed by co-silence of miR-617 (Figure 6E, F).