- Research article
- Open Access
A circRNA-miRNA-mRNA network plays a role in the protective effect of diosgenin on alveolar bone loss in ovariectomized rats
BMC Complementary Medicine and Therapies volume 20, Article number: 220 (2020)
The present study aimed to assess the perturbation in circular RNA (circRNA)/mRNA expression profiles and a circRNA-miRNA-mRNA coexpression network involved in the potential protective effect of diosgenin (DIO) on alveolar bone loss in rats subjected to ovariectomy (OVX).
The Wistar rats (female) manipulated with sham operation were classified as the SHAM group and the grouping of OVX rats administered with DIO, estradiol valerate or vehicle for 12 weeks was DIO group, EV group and OVX group respectively. Following treatments, the plasmatic levels of osteocalcin and tumor necrosis factor-alpha and the microstructure of alveolar bone were assayed. Based on microarray analyses, we identified differentially expressed (DE) circRNAs and mRNAs in alveolar bone of rats in both OVX and DIO group. The DE circRNAs and DE mRNAs involved in the bone metabolism pathway validated by RT-qPCR were considered key circRNAs/mRNAs. On the basis of these key circRNAs/mRNAs, we predicted the overlapping relative miRNAs of key circRNAs/mRNAs, and a circRNA-miRNA-mRNA network was built.
DIO showed an anti-osteopenic effect on the rat alveolar bone loss induced by OVX. In total, we found 10 DE circRNAs (6 downregulated and 4 upregulated) and 614 DE mRNAs (314 downregulated and 300 upregulated) in samples of the DIO group compared with those of the OVX group. However, only one circRNA (rno_circRNA_016717) and seven mRNAs (Sfrp1, Csf1, Il1rl1, Nfatc4, Tnfrsf1a, Pik3c2g, and Wnt9b) were validated by qRT-PCR and therefore considered key circRNA/mRNAs. According to these key circRNA/mRNAs and overlapping predicted miRNAs, a coexpression network was constructed. After network analysis, one circRNA-miRNA-mRNA axis (circRNA_016717/miR-501-5p/Sfrp1) was identified.
The mechanism of DIO inhibiting alveolar bone loss after OVX is possibly relevant to the simultaneous inhibition of osteogenesis and osteoclastogenesis by mediating the expression of important molecules in the Wnt, PI3K, RANK/RANKL or osteoclastogenic cytokine pathways. The circRNA_016717/miR-501-5p/Sfrp1 axis may play important roles in these processes.
Tooth and periodontal diseases associated with alveolar bone loss are high risk complications for and frequently occur in postmenopausal women suffering from osteoporosis because of estrogen deficiency [1, 2]. Unlike other bones, the biological metabolic cycle of an alveolar bone is much shorter , but this bone is more tolerant to estrogen deficiency than appendicular long bones . Some studies have suggested that estrogen can protect against alveolar bone loss not only in humans  but also in rodents .
Many active constituents in plants, especially steroid saponins, are thought to have estrogen-like effects and are called phytoestrogens . Diosgenin (DIO), a phytosteroid sapogenin, is commonly categorized as a phytoestrogen [8,9,10]. The results regarding the estrogen-like effect of DIO have been controversial. A study suggested that DIO (20–200 mg/kg) had no effect on the epithelium height, uterine weight, volume density of endometrium, endometrial/vaginal epithelia or endometrial glands in rats . Nonetheless, other researchers have found that DIO had adverse effects. For example, in ovariectomized (OVX) mice, DIO may stimulate the growth of mammary glands . Regarding bone metabolism, our previous studies showed that decreasing the ratio of receptor activator of nuclear factor-kappa B ligand (RANKL)/osteoprotegerin (OPG) of the tibia  or regulating the expression of long noncoding RNAs in alveolar bone  was the reason for the anti-bone loss effect of DIO in OVX rats.
Compared with linear RNA (e.g., mRNA, miRNA, lncRNA) characterized by 3′ tale and 5′ cap structures, circular RNA (circRNA) features a covalently closed continuous loop generated by a unique splicing strategy . CircRNAs are thought to participate in various diseases because of their special “miRNA sponge” function, whereby they bind miRNAs to efficiently suppress miRNA transcription, and this inhibitory effect on miRNAs can further regulate the expression of downstream mRNAs . In bone metabolism, circRNAs play an important role in modulating mRNA expression via the circRNA-miRNA-mRNA axis not only in a murine model  but also in humans . Thus, we raise the question of whether the regulatory effect of DIO on bone metabolism was mediated via circRNA-miRNA-mRNA interactions. To answer this question, we carried out the present study to explore the action of DIO on gene profiles or the circRNA-miRNA-mRNA network in the alveolar bone of OVX rats.
Experimental animals and treatments
Six-month-old female rats undergoing ovariectomy have been commonly used as a model to study osteoporosis in postmenopausal women [18, 19]. This study included 48 female Wistar rats (six months old). We obtained the rats (average weight 300 ± 20.0 g) from the National Institutes for Food and Drug Control of China. All rats were kept under a regular 12 h/12 h light/dark cycle and at constant room temperature (22 ± 1 °C). The rats underwent sham operation (SHAM, n = 12) or bilateral OVX (n = 36). The ventral approach was used in bilateral ovariectomy . Briefly, at first, the anesthetized rats were fixed ventrally and the abdominal skin in the mid-lower back was preparated and disinfected. A 0.5 cm incision closed the spine under the costal margin was made and skin was separated. Secondly, dorsal muscle was cut and then ovary was exposed adequately. Thirdly, the fat tissue coating the ovary was separated by blunt dissection and the blood vessels and oviduct were ligated. Fourthly, ovariectomy was completed and the preserved tissue or organs were put back into the abdominal cavity, followed by suturing incision. Ovariectomy on the other side was in the same way. The SHAM group rats underwent same operation except that the ovaries were preserved.
The model rats subjected to OVX were classified into 3 categories randomly, namely, OVX group (n = 12), estradiol valerate (EV) group (n = 12) and DIO group (n = 12). According to studies published previously [21,22,23], EV group rats were administered oral gavage of EV (0.1 mg/kg body weight, Bayer Health Care, Guangzhou, China) every day, while DIO group rats were treated by oral gavage with DIO (100 mg/kg body weight, purity≥93%, Sigma-Aldrich, Saint Louis, MO, USA) every day. For the SHAM or OVX group, rats were administered distilled water (of equal volume) by oral gavage. A standard chow was used to feed all rats during our present experiment. The treatments started one week after the operation and lasted 12 weeks. No rats died during the whole treatment period.
At the end of treatment, all experimental rats were anesthetized via intraperitoneal injection by ketamine (80 mg/kg) and xylazine (12 mg/kg) and were subsequently exsanguinated for sacrifice. Under anesthesia, we punctured the abdominal aortae to collect blood (8–10 mL) which was subsequently transferred into tubes added with heparin. The rats subjected to haemospasia were palpated for 5 min to ensure asystole and were considered as dead after confirming asystole, respiration cease and corectasis. Next, the plasma was separated from blood by centrifugation (12 min, 2500 g, 4 °C) and stored (− 80 °C) for the following experiments. To observe and determine the alveolar bone microstructure of the rats, the right mandibles were first excised and kept under the temperature of − 20 °C and then scanned using microcomputed tomography (micro-CT). Thereafter, the right mandibles also served as the specimens for histological examination. For microarray assays and reverse transcription-quantitative polymerase chain reaction (RT-qPCR), we used osseous tissue between the first molar and incisor in left mandibles.
Enzyme-linked immunosorbent assays (ELISA)
We used ELISA to determine plasma concentrations of osteocalcin (OCN) and tumor necrosis factor-alpha (TNF-alpha). A rat OCN ELISA kit (Novus Biologicals, Littleton, CO, USA) and a rat TNF-alpha ELISA kit (Abcam, Cambridge, MA, USA) were used to determine the rat plasma levels of TNF-αand OCN from all groups. The absorption value at 450 nm was detected by a BioTek ELISA reader (BioTek Instruments., Winooski, VT, USA).
The micro-CT scanning method has been previously described . The right mandibles of the rats were scanned by using a high-resolution micro-CT system (Skyscan 1174, Bruker Corporation, Ettlingen, Germany) without any sample processing. The resolution of the scan was 9.8 μm, and analuminum filter (0.5 mm) was used to remove image noise. A global threshold (upper gray threshold: 255, lower gray threshold: 55) was used for the quantity of the parameters about trabecular bone.
The image capture conditions were 800 μA and 50 keV. A cube (0.5 mm × 0.5 mm × 0.5 mm) termed as “volume of interest” (VOI) was reconstructed starting from 0.5 mm beneath the base of crown in the first molar. The trabecular bone morphological characteristics within the VOI were measured applying the Skyscan software package. Furthermore, a 3-D analysis was performed to determine eight key parameters, namely, the trabecular bone volume fraction (BV/TV), the bone surface (BS), the trabecular separation (Tb.Sp), the trabecular number (Tb.N), the trabecular thickness (Tb.Th), the trabecular pattern factor (Tb.Pf), the degree of anisotropy (DA) and the structural model index (SMI) of the identical VOL .
Observations of histology
The prepared solution containing formalin (10%) was used to fix the right mandibles. Next, the mandibles were decalcified with ethylenediaminetetraacetic acid (EDTA, 14%) and embedded in paraffin. After a regular microtome cutting process, the glass slides were used to affix sections, which then were stained with hematoxylin and eosin.
CircRNA and mRNA microarray analysis
Twelve samples of alveolar bones in the DIO group (n = 6) and OVX group (n = 6) were chosen for microarray studies randomly. The microarray assay was completed by KangChen Biotech Inc. (Shanghai, China).
For our circRNA study, we used an Arraystar Rat circRNA Array (Arraystar, Rockville, MD, USA). We used RNase R (Epicentre, Madison, WI, USA) to digested total RNA, and further remove linear RNA and concentrate circular RNA. Subsequently, we amplified and transcribed the concentrated circRNA into fluorescent RNA. This step was performed using an Arraystar Super RNA Labeling Kit (Arraystar, Rockville, MD, USA), by a random priming method. Thereafter, the labeled cRNA was subjected to hybridization, washing and scanning on Arraystar mouse circRNA Array v1.0 using DNA Microarray Scanner (Agilent Technologies, Santa Clara, CA, USA).
For mRNAs, we used a rat 4 × 44 K Gene Expression Array (Agilent Technology, Santa Clara, CA, USA). In line with the protocol, samples were labeled and arrays were hybridized. Then, the obtained array images were evaluated by the Agilent Feature Extraction software (version 22.214.171.124). The data was normalized and further processed using an R software package.
CircRNAs/mRNAs with fold change≥1.5 and P-value< 0.05 were regarded as differentially expressed (DE) circRNAs/mRNAs.
Analysis of ingenuity pathway analysis (IPA)
All DE mRNAs obtained from our array analyses were input into the IPA system, a system that predicts canonical pathways, global functions and biological networks of a particular gene dataset and is referred to as the Ingenuity Pathways Knowledge Base (IPKB). In IPKB, “Role of Osteoblasts, Osteoclasts and Chondrocytes in Rheumatoid Arthritis (ROOCRA)” pathway is a exclusive pathway because this pathway includes nearly all known pivotal molecules and signaling pathways that are closely related to the development, differentiation, degradation, mineralization and apoptosis of osteoblasts and osteoclasts. As this pathway is essential in bone metabolism, the DE mRNAs that belong to the ROOCRA pathway were therefore regarded as key DE mRNAs.
Building of a circRNA-miRNA-mRNA coexpression network
Potential target miRNAs of circRNAs were predicted by miRanda (v3.3a)  and TargetScan (version 7.2) . Software developed by Arraystar was applied to predict the interactive relationship between circRNAs and miRNAs, while the predicted interactions of mRNAs and miRNAs were obtained using miRanda (version 3.3a). After Pearson correlation analyses of the predicted mRNA-miRNA and circRNA-miRNA coexpression networks, the circRNA-miRNA-mRNA coexpression network were constructed. Cytoscape (version 2.8.2)  was applied to visualize the whole circRNA-miRNA-mRNA network.
Validations of DE circRNAs, key DE mRNAs and predicted miRNAs using RT-qPCR
qRT-PCR assays used for detecting the expression levels of the predicted miRNAs were all performed in an ABI 7500 system (Applied Biosystems, Foster City, CA, USA) using SYBR RT-PCR kits (Takara, Dalian, Liaoning, China). We used the method of 2-ΔΔCt cycle threshold, and the expression of U6 served as the internal normalization control. For circRNAs/mRNAs expression data, the internal control was Gapdh expression.
The mean ± standard deviation was used to express all values. We utilized SPSS 19.0 (SPSS Inc., Chicago, IL, USA) for analyzing data statistically. The difference of the assessed parameters between rats from two groups or four groups was tested using the the t-test or the one-way analysis of variance (ANOVA) and subsequently the least significant difference (LSD) test, respectively. Kolmogorov-Smirnov statistics was used to test the normality of data from all groups. Statistical significance was set at p < 0.05.
Regulation on plasmic OCN and TNF-alpha by DIO
Figure 1A and B showed the plasmic concentrations of OCN and TNF-alpha in the rats that were subjected to 12-week long treatment. After the treatment, the OCN and TNF-alpha contents of rats in the OVX group were remarkably superior to those of rats in the SHAM group (p < 0.01). Contrastively, the DIO and EV group rats exhibited significantly lower plasmic OCN and TNF-alpha levels than the OVX group rats (p < 0.05).
Effects of DIO on the microstructure of trabecular bones
Our micro-CT results suggested that in the OVX group, morphological parameters of bone, such as BV/TV, BS, Tb.Th, and Tb.N were substantially reduced, yet Tb.Pf, Tb.Sp and SMI were raised in comparison with the SHAM group (Fig. 2A-H). As shown in Fig. 2, significant changes in the bone morphological parameters were found when the rats were treated with EV or DIO. Furthermore, trabecular impairment due to ovariectomy was lessened by DIO or EV treatment (Fig. 3A-D).
Effect of DIO on histomorphological change of alveolar bone
To understand the effect of DIO from a histological perspective, we observed the alveolar bone beneath the first molar. Representative histomorphological images from four group rats are shown in Fig. 4. It is clear that the SHAM group rats had thicker trabeculae in alveolar bone and scant and smaller medullary cavity (Fig. 4A) in comparison to the OVX group rats. It appears that OVX decreases alveolar bone volume and increases the medullary cavity size (Fig. 4B). Furthermore, bone volume beneath the first molar was significantly enlarged after DIO or EV treatment (Fig. 4C and D), and EV appeared to have a stronger effect in this regard. With these data, our results from histological observations and micro-CT were consistent with each other.
Regulation of the expression profiles of circRNA and mRNAs by DIO
Our microarray results revealed that there were, in total, 10 DE exonic circRNAs (Table 1) and 614 DE mRNAs (Table S4) in the alveolar bone between DIO and OVX group rats. Four and 6 circRNAs were up-and downregulated, respectively, while almost half of the mRNAs were increased and half were decreased (300 were upregulated and 314 downregulated).
Regulation of the bone metabolic pathway by DIO
Eight key DE mRNAs detected from the groups of OVX and DIO were assigned to the pathway of ROOCRA in IPA (Table 2). Figures 5 and 6 illustrate the potential effects of these key mRNAs on osteoclasts and osteoblasts, respectively, on the basis of the ROOCRA pathway.
Validation of DE circRNAs and key DE mRNAs using RT-qPCR
The 10 circRNA expressions are listed in Tables 1 and 8 key DE mRNA expressions listed in Table 2 were measured. As a result, rno_circRNA_016717 was the only circRNA whose RT-qPCR data were in agreement with the microarray results (Fig. 7). As for mRNAs, RT-qPCR data of all mRNAs except Birc3 agreed with the microarray results (Fig. 8). Tables 3-4 listed the primers of DE circRNAs and key DE mRNAs for qRT-PCR experiments.
A circRNA-miRNA-mRNA coexpression network showing the modulatory effect of DIO
Based on the validated DE circRNA (rno_circRNA_016717) and 7 DE key mRNAs (Sfrp1, Pik3c2g, Wnt9b, Csf1, Il1rl1, Nfatc4, and Tnfrsf1a), 60 potential miRNA targets were predicted using miRanda and TargetScan. Thereafter, a circRNA-miRNA-mRNA coexpression network was established that included 380 nodes and 2173 edges (Fig. S1).
According to the interactions among circRNA, miRNAs and mRNAs, we searched for overlapping miRNAs that were downstream molecules of the DE circRNA (rno_circRNA_016717) and were upstream molecules of one DE upregulated mRNA (Sfrp1). Finally, only one circRNA-miRNA-mRNA axis, circRNA_016717/miR-501-5p/Sfrp1 (Fig. 9), was explored. The roles of this axis in the regulatory effect of DIO require further intensive study.
Validation of the predicted miRNA by RT-qPCR
MiR-501-5p, the only predicted miRNA, was validated by RT-qPCR. As predicted, miR-501-5p expression in alveolar bone of DIO group rats was decreased compared to that in OVX rats (Fig. 10). Tables 5 listed the primers of predicted miRNA for qRT-PCR experiments.
DIO is a phytoestrogen that has been shown to have an anti-osteopenic function [29, 30]; however, the mechanism of action remains unclear. Since Hansen TB et al. demonstrated that circRNAs can play a role as “miRNA sponge” , the important regulatory role of circRNAs on miRNAs has received increasing attention. The circRNA profile and the circRNA-miRNA-mRNA network in the protective effect of DIO on alveolar bone in OVX rats has never been reported, so we conducted a study combining microarray analysis and bioinformatics to fill the gap.
In our present study, we chose a rat model induced by OVX to mimic postmenopausal women with alveolar high-turnover bone loss. These patients usually feature high bone resorption and bone formation processes [32, 33]. Our 3-D bone microstructure modeling revealed many parameters, e.g., Tb.Th, BV/TV, BS, Tb.N, Tb.Pf, Tb.Sp, as well as SMI were significantly changed. The observation suggested that DIO and EV group rats exhibited less bone loss in alveolar bones in comparison to the model rats in the OVX group. It is worth noting that the anti-bone loss action of DIO on alveolar bone was less impactful than that of EV (Fig. 2). Our histomorphological findings in alveolar bone (Fig. 4) were largely in line with the results from micro-CT (Fig. 3).
Both TNF-alpha and OCN were upregulated in the OVX group compared to the SHAM group. However, after treatment with EV or DIO, the increased levels of the TNF-alpha and OCN were reduced (Fig. 1A and B). Our findings on EV are consistent with those in previous reports [34, 35] and strongly suggested that DIO, or other estrogens, may have a negative effect on both bone resorption and bone formation in OVX rats; thus, we confirmed that DIO can act as an estrogen-like chemical.
It is interesting to note that our histological observations showed that DIO had an anti-bone loss influence on rat alveolar bones induced by OVX, and this effect was confirmed by results from micro-CT and assay of TNF-alpha and OCN.
We conducted microarray assays to understand the role of perturbations of circRNA and mRNA profiles in the anti-bone loss action of DIO. Due to poor understanding about the function of DE circRNAs, we focused on the DE mRNAs involved in the anti-osteopenic effects of DIO. Seven validated DE mRNAs (Sfrp1, Csf1, Il1rl1, Nfatc4, Tnfrsf1a, Pik3c2g and Wnt9b) were associated with the ROOCRA pathway. It is clear that DIO has a multitarget inhibitory effect on the signaling pathways of both bone resorption and bone formation (Figs. 5 and 6).
The Wnt pathway is acknowledged to be one of the key signaling pathways that mediate the osteogenic differentiation not only mesenchymal stem cells but preosteoblasts [36, 37], whereas RANKL/the receptor activator of nuclear factor-kappa B (RANK) pathway plays an essential role in osteoclastogenic differentiation and bone resorption [38, 39]. Interestingly, sFRPs have dual inhibitory effects on the two pathways above. sFRPs can function as Wnt inhibitors to resemble the ligand-binding cysteine rich domain (CRD) of the Frizzled family of Wnt receptors and inhibit both canonical Wnt/beta-catenin signaling and noncanonical Wnt/planar cell polarity (PCP) signaling [40, 41]. In osteoclasts, sFRP1 can bind to RANKL directly and further block the interaction of RANKL/RANK and osteoclastogenesis . DIO has been reported to inhibit breast cancer stem-like cells by deregulating the activation of Wnt/beta-catenin signaling via sFRP4 . We speculated that sFRP1 could be a key target of DIO in its anti-bone loss effect.
Phosphatidylinositol 3-kinase (PI3K) controls numerous cellular functions like motility after ligand activation and cell proliferation . In osteoblasts, PI3K has been well studied. PI3K inhibits osteoblast apoptosis through activating Akt (protein kinase B) and further activates PI3K/Akt by Wnt3a, and heparin promotes osteoblast differentiation . For osteoclastogenesis, the PI3K/Akt pathway also plays a fundamental role. Some researchers recently reported that the activation of Akt can limit osteoclast differentiation through activating the glycogen synthase kinase 3 beta (GSK 3beta)/nuclear factor of activated T cells (NFAT) c1 signaling cascade . PI3K-gamma (Pik3c2g) was also characterized recently. It has been shown that this protein is associated with the promotion of osteoclastogenesis and bone mass reduction in mice and therefore may be a potential target in osteoporosis . In our present study, we inferred that DIO could attenuate osteogenesis and osteoclastogenesis by decreasing the expression of PI3K in osteoblasts and osteoclasts at the mRNA level.
Our data also showed that DIO could slow down the process of bone resorption by decreasing the signaling of three potent stimulator pathways in osteoclastogenesis (Fig. 6). DIO attenuated the ligands or receptors of macrophage colony-stimulating factor (M-CSF), interleukin-1 receptor (IL-1R) and tumor necrosis factor receptor (TNF-R) [48,49,50]. In particular, the results on TNF-R1 and IL-1R from the microarray assay were consistent with those of our previous study . In addition, NFATC is regarded as a master transcription factor indispensable for the osteoclastogenesis induced by RANKL. The auto amplification and activation of NFATc1 give rise to a rapid upregulation of osteoclast-specific genes [51, 52] and NFATc4 was necessary for the osteoclastogenic effect of NFATc1 . In this study, we surmised that TNF-R1, IL-1R and NFATc4 could be other target molecules of DIO in the reduction of bone resorption.
To sum up, the findings from the gene chip and pathway analyses suggested that the anti-bone loss action of diosgenin on alveolar bone was ascribed to inhibition of osteogenesis and osteoclastogenesis synchronously by mediating the expression of important molecules in the Wnt, PI3K, RANK/RANKL or osteoclastogenic cytokine pathways (e.g., Sfrp1, Pik3c2g, Wnt9b, Csf1, Il1rl1, Nfatc4, and Tnfrsf1a).
We built a circRNA-miRNA-mRNA coexpression network involved in the inhibitory action of DIO on alveolar bone loss. This network is beneficial for identifying pivotal miRNAs that are linked to the regulatory effect of rno_circRNA_016717 on the 7 key mRNAs mentioned above.
There were 380 nodes and 2173 edges in the coexpression network (Fig. S1). However, we only found one circRNA-miRNA-mRNA axis, circRNA_016717/miR-501-5p/Sfrp1 (Fig. 9). In 2016, a study reported that miR-501-5p was considerably upregulated in human gastric cancer tissues and cell lines. The authors found that miR-501-5p may directly bind and suppress several important repressors of the Wnt/beta-catenin signaling cascade [such as GSK3 beta, dickkopf-related protein 1 (DKK1) and naked cuticle 1 (NKD1)], which lead to hyperactivated signaling in gastric cancer cells . No study has reported the relationship between miR-501-5p and Sfrp1 in vivo or in vitro, but we hypothesize that this relationship likely exists considering that Sfrp1 is also a repressor of Wnt/beta-catenin signaling. Further efforts are needed with the aim of revealing the potential mechanism of the circRNA_016717/miR-501-5p/Sfrp1 axis in the anti-osteoporotic effect of DIO.
Our study proved that diosgenin had a protective effect on rat alveolar bone loss induced by ovariectomized. The potential mechanism of this protective effect from diosgenin was possibly associated with circRNA_016717/miR-501-5p/Sfrp1 axis that could inhibit osteogenesis and osteoclastogenesis simultaneously by regulating the expression of important molecules in the Wnt, PI3K, RANK/RANKL or osteoclastogenic cytokine pathways.
Availability of data and materials
The current publication contains all data produced or analyzed in the study. The corresponding author will provide associated raw data which was used and/or analyzed throughout this study if requested reasonably.
One-way analysis of variance
Trabecular bone volume fraction
Cysteine rich domain
Degree of anisotropy
Dickkopf-related protein 1
Enzyme-linked immunosorbent assays
- GSK 3beta:
Glycogen synthase kinase 3 beta
Ingenuity Pathways Knowledge Base
Least significant difference
Macrophage colony-stimulating factor
Nuclear factor of activated T cells
Naked cuticle 1
Planar cell polarity
Receptor activator of nuclear factor-kappa B ligand
Receptor activator of nuclear factor-kappa B
Role of Osteoblasts, Osteoclasts and Chondrocytes in Rheumatoid Arthritis
Reverse transcription-quantitative polymerase chain reaction
Structural model index
Trabecular pattern factor
Tumor necrosis factor-alpha
Tumor necrosis factor receptor
Volume of interest
Taguchi A, Tanimoto K, Suei Y, Wada T. Tooth loss and mandibular osteopenia. Oral surgery, oral medicine, oral pathology, oral radiology, and endodontics. 1995;79:127–32.
Krall EA, Garcia RI, Dawson-Hughes B. Increased risk of tooth loss is related to bone loss at the whole body, hip, and spine. Calcif Tissue Int. 1996;59:433–7.
Ikeo T, Goda S, Domae E. Metabolism of alveolar bone. Clinical calcium. 2006;16:117–21.
Du Z, Steck R, Doan N, Woodruff MA, Ivanovski S, Xiao Y. Estrogen deficiency-associated bone loss in the maxilla: a methodology to quantify the changes in the maxillary intra-radicular alveolar bone in an Ovariectomized rat osteoporosis model. Tissue engineering Part C, Methods. 2015;21:458–66.
Payne JB, Zachs NR, Reinhardt RA, Nummikoski PV, Patil K. The association between estrogen status and alveolar bone density changes in postmenopausal women with a history of periodontitis. J Periodontol. 1997;68:24–31.
Duarte PM, Goncalves P, Casati MZ, de Toledo S, Sallum EA, Nociti FH Jr. Estrogen and alendronate therapies may prevent the influence of estrogen deficiency on the tooth-supporting alveolar bone: a histometric study in rats. J Periodontal Res. 2006;41:541–6.
Turner JV, Agatonovic-Kustrin S, Glass BD. Molecular aspects of phytoestrogen selective binding at estrogen receptors. J Pharm Sci. 2007;96:1879–85.
Sirotkin AV, Alexa R, Alwasel S, Harrath AH. The phytoestrogen, diosgenin, directly stimulates ovarian cell functions in two farm animal species. Domest Anim Endocrinol. 2019;69:35–41.
Shen M, Qi C, Kuang YP, Yang Y, Lyu QF, Long H, Yan ZG, Lu YY. Observation of the influences of diosgenin on aging ovarian reserve and function in a mouse model. Eur J Med Res. 2017;22:42.
Zhang Z, Chen Y, Xiang L, Wang Z, Xiao GG, Ju D. Diosgenin protects against alveolar bone loss in ovariectomized rats via regulating long non-coding RNAs. Experimental and therapeutic medicine. 2018;16:3939–50.
Medigovic I, Ristic N, Zivanovic J, Sosic-Jurjevic B, Filipovic B, Milosevic V, Nestorovic N. Diosgenin does not express estrogenic activity: a uterotrophic assay. Can J Physiol Pharmacol. 2014;92:292–8.
Aradhana, Rao AR, Kale RK. Diosgenin--a growth stimulator of mammary gland of ovariectomized mouse. Indian journal of experimental biology. 1992;30:367–70.
Zhang Z, Song C, Fu X, Liu M, Li Y, Pan J, Liu H, Wang S, Xiang L, Xiao GG, et al. High-dose diosgenin reduces bone loss in ovariectomized rats via attenuation of the RANKL/OPG ratio. Int J Mol Sci. 2014;15:17130–47.
Salzman J, Gawad C, Wang PL, Lacayo N, Brown PO. Circular RNAs are the predominant transcript isoform from hundreds of human genes in diverse cell types. PLoS One. 2012;7:e30733.
Wilusz JE, Sharp PA. Molecular biology. A circuitous route to noncoding RNA. Science. 2013;340:440–1.
Chen X, Ouyang Z, Shen Y, Liu B, Zhang Q, Wan L, Yin Z, Zhu W, Li S, Peng D. CircRNA_28313/miR-195a/CSF1 axis modulates osteoclast differentiation to affect OVX-induced bone absorption in mice. RNA Biol. 2019:1–14.
Yu L, Liu Y. circRNA_0016624 could sponge miR-98 to regulate BMP2 expression in postmenopausal osteoporosis. Biochemical and biophysical research communications. 2019.
Li CM, Dong XL, Fan XD, Wu JH, Wang QH, Tian XL, Guo DJ, Wong MS, Qiu TQ, Chan SW. Aqueous extract of danshen (Salvia miltiorrhiza Bunge) protects ovariectomized rats fed with high-fat diet from endothelial dysfunction. Menopause. 2013;20:100–9.
Maimoun L, Brennan-Speranza TC, Rizzoli R, Ammann P. Effects of ovariectomy on the changes in microarchitecture and material level properties in response to hind leg disuse in female rats. Bone. 2012;51:586–91.
Lane NE, Yao W, Kinney JH, Modin G, Balooch M, Wronski TJ. Both hPTH(1-34) and bFGF increase trabecular bone mass in osteopenic rats but they have different effects on trabecular bone architecture. Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research. 2003;18:2105–15.
Hsu KH, Chang CC, Tsai HD, Tsai FJ, Hsieh YY. Effects of yam and diosgenin on calpain systems in skeletal muscle of ovariectomized rats. Taiwanese journal of obstetrics & gynecology. 2008;47:180–6.
Ma Y, Zeng R, Hu QQ, Yan HX, Yang LX, Dong Y, Qu Y. Preventive effects of Polygonum orientale L. on ovariectomy-induced osteoporosis in rats. Climacteric : the journal of the International Menopause Society. 2020;1–9.
Tang YQ, Li C, Sun XJ, Liu Y, Wang XT, Guo YB, Wang LL, Ma RF, Niu JZ, Fu M, et al. Fructus Ligustri Lucidi modulates estrogen receptor expression with no uterotrophic effect in ovariectomized rats. BMC Complement Altern Med. 2018;18:118.
Yang J, Pham SM, Crabbe DL. High-resolution micro-CT evaluation of mid- to long-term effects of estrogen deficiency on rat trabecular bone. Acad Radiol. 2003;10:1153–8.
Bouxsein ML, Boyd SK, Christiansen BA, Guldberg RE, Jepsen KJ, Muller R. Guidelines for assessment of bone microstructure in rodents using micro-computed tomography. Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research. 2010;25:1468–86.
Betel D, Wilson M, Gabow A, Marks DS, Sander C. The microRNA.org resource: targets and expression. Nucleic Acids Res. 2008;36:D149–53.
Agarwal V, Bell GW, Nam JW, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. eLife. 2015;4.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.
Hung YT, Tikhonova MA, Ding SJ, Kao PF, Lan HH, Liao JM, Chen JH, Amstislavskaya TG, Ho YJ. Effects of chronic treatment with diosgenin on bone loss in a d-galactose-induced aging rat model. The Chinese journal of physiology. 2014;57:121–7.
Zhao S, Niu F, Xu CY, Liu Y, Ye L, Bi GB, Chen L, Tian G, Nie TH. Diosgenin prevents bone loss on retinoic acid-induced osteoporosis in rats. Ir J Med Sci. 2016;185:581–7.
Hansen TB, Jensen TI, Clausen BH, Bramsen JB, Finsen B, Damgaard CK, Kjems J. Natural RNA circles function as efficient microRNA sponges. Nature. 2013;495:384–8.
Alikhani M, Alikhani M, Alansari S, Almansour A, Hamidaddin MA, Khoo E, Lopez JA, Nervina JM, Nho JY, Oliveira SM, et al. Therapeutic effect of localized vibration on alveolar bone of osteoporotic rats. PLoS One. 2019;14:e0211004.
Wang L, Cheng L, Zhang B, Wang N, Wang F. Tanshinone prevents alveolar bone loss in ovariectomized osteoporosis rats by up-regulating phosphoglycerate dehydrogenase. Toxicol Appl Pharmacol. 2019;376:9–16.
Rao YQ, Li J, Wang WJ. Effects of Gengnianchun on learning and memory ability, neurotransmitter, cytokines, and leptin in ovariectomized rats. Int J Clin Exp Med. 2015;8:8648–60.
Cheng M, Wang Q, Fan Y, Liu X, Wang L, Xie R, Ho CC, Sun W. A traditional Chinese herbal preparation, Er-Zhi-Wan, prevent ovariectomy-induced osteoporosis in rats. J Ethnopharmacol. 2011;138:279–85.
Marcellini S, Henriquez JP, Bertin A. Control of osteogenesis by the canonical Wnt and BMP pathways in vivo: cooperation and antagonism between the canonical Wnt and BMP pathways as cells differentiate from osteochondroprogenitors to osteoblasts and osteocytes. BioEssays : news and reviews in molecular, cellular and developmental biology. 2012;34:953–62.
Lin GL, Hankenson KD. Integration of BMP, Wnt, and notch signaling pathways in osteoblast differentiation. J Cell Biochem. 2011;112:3491–501.
Amin N, Boccardi V, Taghizadeh M, Jafarnejad S. Probiotics and bone disorders: the role of RANKL/RANK/OPG pathway. Aging Clin Exp Res. 2019.
Hanada R, Hanada T, Penninger JM. Physiology and pathophysiology of the RANKL/RANK system. Biol Chem. 2010;391:1365–70.
Li Y, Rankin SA, Sinner D, Kenny AP, Krieg PA, Zorn AM. Sfrp5 coordinates foregut specification and morphogenesis by antagonizing both canonical and noncanonical Wnt11 signaling. Genes Dev. 2008;22:3050–63.
Satoh W, Matsuyama M, Takemura H, Aizawa S, Shimono A. Sfrp1, Sfrp2, and Sfrp5 regulate the Wnt/beta-catenin and the planar cell polarity pathways during early trunk formation in mouse. Genesis. 2008;46:92–103.
Hausler KD, Horwood NJ, Chuman Y, Fisher JL, Ellis J, Martin TJ, Rubin JS, Gillespie MT. Secreted frizzled-related protein-1 inhibits RANKL-dependent osteoclast formation. Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research. 2004;19:1873–81.
Bhuvanalakshmi G, Basappa, Rangappa KS, Dharmarajan A, Sethi G, Kumar AP, Warrier S. Breast Cancer Stem-Like Cells Are Inhibited by Diosgenin, a Steroidal Saponin, by the Attenuation of the Wnt beta-Catenin Signaling via the Wnt Antagonist Secreted Frizzled Related Protein-4. Frontiers in pharmacology. 2017;8:124.
Cantley LC. The phosphoinositide 3-kinase pathway. Science. 2002;296:1655–7.
Ling L, Dombrowski C, Foong KM, Haupt LM, Stein GS, Nurcombe V, van Wijnen AJ, Cool SM. Synergism between Wnt3a and heparin enhances osteogenesis via a phosphoinositide 3-kinase/Akt/RUNX2 pathway. J Biol Chem. 2010;285:26233–44.
Moon JB, Kim JH, Kim K, Youn BU, Ko A, Lee SY, Kim N. Akt induces osteoclast differentiation through regulating the GSK3beta/NFATc1 signaling cascade. J Immunol. 2012;188:163–9.
Kang H, Chang W, Hurley M, Vignery A, Wu D. Important roles of PI3Kgamma in osteoclastogenesis and bone homeostasis. Proc Natl Acad Sci U S A. 2010;107:12901–6.
Hattersley G, Owens J, Flanagan AM, Chambers TJ. Macrophage colony stimulating factor (M-CSF) is essential for osteoclast formation in vitro. Biochem Biophys Res Commun. 1991;177:526–31.
Kudo O, Fujikawa Y, Itonaga I, Sabokbar A, Torisu T, Athanasou NA. Proinflammatory cytokine (TNFalpha/IL-1alpha) induction of human osteoclast formation. J Pathol. 2002;198:220–7.
Kwan Tat S, Padrines M, Theoleyre S, Heymann D, Fortun Y. IL-6, RANKL, TNF-alpha/IL-1: interrelations in bone resorption pathophysiology. Cytokine Growth Factor Rev. 2004;15:49–60.
Takayanagi H. Mechanistic insight into osteoclast differentiation in osteoimmunology. J Mol Med (Berl). 2005;83:170–9.
Sitara D, Aliprantis AO. Transcriptional regulation of bone and joint remodeling by NFAT. Immunol Rev. 2010;233:286–300.
Kuroda Y, Hisatsune C, Mizutani A, Ogawa N, Matsuo K, Mikoshiba K. Cot kinase promotes Ca2+ oscillation/calcineurin-independent osteoclastogenesis by stabilizing NFATc1 protein. Mol Cell Biol. 2012;32:2954–63.
Fan D, Ren B, Yang X, Liu J, Zhang Z. Upregulation of miR-501-5p activates the wnt/beta-catenin signaling pathway and enhances stem cell-like phenotype in gastric cancer. Journal of experimental & clinical cancer research : CR. 2016;35:177.
This research was funded by the National Natural Science Foundation of China (No. 81473450, 81302960) in the design of the study and collection, analysis, and interpretation of data, the Beijing Natural Science Foundation (No. 7172148) in collection, analysis, and interpretation of data, the Top Young Scientist Funds of Beijing University of Chinese Medicine (No. BUCM-2019-JCRC007) in collection of data, as well as the Fundamental Research Funds for the Central Public Welfare Research Institutes (No. YZ-1780, YZ-202021) in writing the manuscript and data collection.
Ethics approval and consent to participate
The ethics committee of Institute of Basic Theory, China Academy of Chinese Medical Sciences provided the approval to the entire protocol of this study which adhered to the recommendations mentioned in the Guidelines for Laboratory Animal Care and Use from the Chinese Ministry of Science and Technology (2006) (Approval number: 2017–015).
Consent for publication
No conflict of interest was declared by the authors.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The differential expressions of mRNAs. In the current file, 614 mRNAs with differential expressions (p ≤ 0.05, fold changes ≥1.5) extracted from samples of alveolar bones from both DIO and OVX group were presented.
A co-expression network of circRNA-miRNA-mRNA related to the roles exerted by DIO on alveolar bone loss among OVX rats. The green nodes represented key mRNAs, the blue nodes were the related mRNAs introduced into the network, the brown node represented differentially expressed circRNAs, and the red nodes represented predicted miRNAs. The lines indicate a correlation between mRNAs, circRNA and miRNAs.
About this article
Cite this article
Zhang, Z., Yue, L., Wang, Y. et al. A circRNA-miRNA-mRNA network plays a role in the protective effect of diosgenin on alveolar bone loss in ovariectomized rats. BMC Complement Med Ther 20, 220 (2020). https://doi.org/10.1186/s12906-020-03009-z
- Alveolar bone
- Osteoprotective effect
- Coexpression network