Exosomal hsa_circ_0006859 is a potential biomarker for postmenopausal osteoporosis and enhances adipogenic versus osteogenic differentiation in human bone marrow mesenchymal stem cells by sponging miR-431-5p
Stem Cell Research & Therapy volume 12, Article number: 157 (2021)
As one of the most common chronic diseases in the world, osteoporosis occurs especially in postmenopausal women. Circular RNAs (circRNAs) are emerging as major drivers in human disease. The aim of the present study was to analyse circRNA expression profiles in osteoporosis and to explore the clinical significance and the regulatory molecular mechanism of hsa_circ_0006859 during osteoporosis.
Exosomes were isolated from clinically collected serum samples. A circRNA microarray was performed to screen differentially expressed circRNAs. Quantitative real-time PCR (qRT-PCR) and western blot were performed to analyse target gene mRNA expression and protein expression. Alizarin red staining (ARS) was performed to evaluate the mineralization ability of human bone marrow mesenchymal stem cells (hBMSCs). Oil Red O staining was performed to evaluate the lipid droplet formation ability of hBMSCs. Bioinformatics analysis and the luciferase reporter assay were performed to investigate the interaction between two genes.
Hsa_circ_0006859 was identified as one of the most upregulated circRNAs in the microarray analysis. Hsa_circ_0006859 in exosomes was upregulated in osteoporosis patients compared to healthy controls. Hsa_circ_0006859 differentiated osteopenia or osteoporosis patients from healthy controls with high sensitivity and specificity. Hsa_circ_0006859 suppressed osteoblastic differentiation and promoted adipogenic differentiation of hBMSCs. Hsa_circ_0006859 directly bound to miR-431-5p, and ROCK1 was identified as a novel target gene of miR-431-5p. Hsa_circ_0006859 is a competing endogenous RNA (ceRNA) of miR-431-5p that promotes ROCK1 expression. Hsa_circ_0006859 suppressed osteogenesis and promoted adipogenesis by sponging miR-431-5p to upregulate ROCK1.
Exosomal hsa_circ_0006859 is a potential biomarker for postmenopausal osteoporosis and controls the balance between osteogenesis and adipogenesis in hBMSCs by sponging miR-431-5p.
Osteoporosis is a prevalent age-related disease characterized by bone mass reduction and tissue microstructure degeneration, which lead to a high risk of bone fracture, especially in postmenopausal women . Osteoporosis is widely believed to be the consequence of an imbalance between bone formation, which is regulated by osteoblasts, and bone resorption, which is regulated by osteoclasts . Though osteoporosis management strategies have made some strides in the past few years, the complex molecular mechanisms causing osteoporosis and the lack of potential therapeutic targets hinder improvements in its prevention and treatment . Thus, exploring novel targets for the treatment of osteoporosis is imperatively needed.
Circular RNAs (circRNAs) are a novel class of endogenous covalently linked RNA molecules without 5′ to 3′ polarity or a poly A tail . CircRNAs are evolutionarily conserved and are abundant in eukaryotes, suggesting their importance in various physiological conditions . CircRNAs are stable and enriched in exosomes, indicating the potential of circRNAs as promising biomarkers for complex diseases [6, 7]. CircRNAs have multiple functions, such as modulating gene transcription, regulating alternative splicing, modulating translation, interacting with RNA-binding proteins (RBPs), and functioning as microRNA (miRNA) sponges . MicroRNAs (miRNAs) are evolutionarily conserved small non-coding RNAs approximately 22 nucleotides long that inhibit gene expression by base-pairing with the 3′-untranslated region (3′-UTR) or 5′-UTR of mRNAs . Due to the special circular structure, circRNAs are highly tolerant to exonucleases in body fluids, which makes them feasible as special molecular markers in diseases . CircRNAs are also involved in diverse physiological programmes, and the dysregulation of circRNAs is associated with a wide range of diseases, including osteoporosis . For example, hsa_circ_0002060 in serum/plasma was significantly upregulated in osteoporosis patients and was associated with low bone mineral density (BMD). Furthermore, hsa_circ_0002060 showed good potential for diagnosing osteoporosis with a high sensitivity and specificity . Hsa_circ_0006393 promotes osteogenesis by sponging miR-145-5p and upregulating FOXO1 in osteoporosis . Hsa_circ_0016624 prevents osteoporosis and promotes BMP2 expression by sponging miR-98 . Hsa_circ_0076906 promotes osteogenesis and alleviates osteoporosis by acting as a sponge for miR-1305 and upregulating osteoglycin . With the development of different screening methods for circRNAs, such as circRNA microarrays, circRNA sequencing, Northern blotting, and PCR-based analyses, increasing circRNAs are being identified as involved in the progression of various diseases. However, only a few circRNAs have been investigated in osteoporosis. Therefore, it is imperative to investigate the association between circRNA expression and osteoporosis development.
In the present study, we compared circRNA expression profiles of exosomes isolated from the serum of osteoporosis patients and matched healthy controls using circRNA microarray technology. Hsa_circ_0006859, one of the most upregulated circRNAs as revealed by microarray analysis, was investigated to explore its clinical significance and potential regulatory mechanism in osteoporosis. Hsa_circ_0006859 was reported to be upregulated during chondrogenic differentiation of human mesenchymal stromal cells by interacting with five possible miRNAs . Hsa_circ_0006859 was reported to be upregulated in colorectal cancer . Hsa_circ_0006859 was also reported to be downregulated in bladder cancer . However, these reports barely discuss the role of hsa_circ_0006859 on miRNAs and the related cell signalling in their research on human diseases. In the present study, our results indicated that exosomal circRNAs were differentially expressed between osteoporosis patients and healthy controls. Therefore, Hsa_circ_0006859 is a potential biomarker for postmenopausal osteoporosis. Furthermore, hsa_circ_0006859 controls the balance between osteogenesis and adipogenesis by sponging miR-431-5p to regulate ROCK1.
Patients and samples
This study was approved by the Ethics Committee of the Third Affiliated Hospital of Soochow University (No.CZYY2020032). Each participant signed an informed consent document. All participants were postmenopausal women and were consecutively recruited at Osteoporosis Outpatient Clinic in the Third Affiliated Hospital of Soochow University from January 2017 to June 2019. Every participant signed the informed consent. Lumbar spine L1-L4 BMD was evaluated by dual energy X-ray absorptiometry (DXA) (Hologic Discovery Wi, Hologic, USA) according to standard operation instructions. T-score is the standard deviation of the BMD value from the mean. Recruited participants were classified as the osteoporosis group, defined by lumbar spine (L1-L4) BMD T-score ≤ − 2.5; the osteopenia group, defined by lumbar spine (L1-L4) BMD T-score ≤ − 1.0 and > − 2.5; and the control group, defined by lumbar spine (L1-L4) BMD T-score ≥ − 1.0. All recruited participants were newly diagnosed with osteoporosis or osteopenia. Participants with other relevant chronic diseases were excluded. The detailed characteristics of study subjects are summarized in Table 1. Thirty of the 58 participants in the osteoporosis group received anti-osteoporotic treatment using alendronate supplemented with calcium and vitamin D for 9 months. Peripheral blood was collected from each participant at initial diagnosis and from each participant who received anti-osteoporotic treatment at 9 months. Serum was centrifuged and collected after the blood was left for 1 h at room temperature to allow clotting. Aliquots were stored at − 80 °C until use.
Assays for parathyroid hormone (PTH), 25-hydroxyvitamin D (25(OH)D), β-isomerized C-terminal telopeptides (β-CTx), total N-terminal procollagen of type l collagen (tP1NP), osteocalcin, calcium, and inorganic phosphorus in serum were measured by the Department of Laboratory Medicine in the hospital. PTH was analysed on UniCel DxI 800 (Beckman Coulter, USA). 25(OH) D, β-CTx, tP1NP, and osteocalcin were analysed on Cobas E601 (Roche Diagnostics, Mannheim, Germany). Calcium and inorganic phosphorus were analysed on AU5831 (Beckman Coulter, USA).
Exosomes were isolated from serum using the miRCURY Exosome Serum/Plasma Kit (Thermo Fisher, USA) according to the manufacturer’s instructions. Then, 100 ng synthetic cel-miR-39-3p (Thermo Fisher, USA) was added as a spike-in control similar to previous reports [17, 18]. Briefly, 500 μl serum was incubated with the precipitation buffer for 1 h at 4 °C. The mixture was fully vortexed and then centrifuged at 1500g for 30 min at room temperature, and the supernatant was discarded to obtain the exosome precipitate. Purified exosomes were resuspended by adding resuspension buffer.
The circRNA microarray assay was performed on Human CircRNA Array V2 (CapitalBio Technology, Beijing, China) by CapitalBio (CapitalBio Technology, Beijing, China) as previously reported . Briefly, total RNA was extracted from the exosomes, which were isolated from the serum of 3 osteoporosis patients and 3 matched controls. Extracted RNAs were digested, amplified, labelled with Cy3-dCTP, and hybridized to the microarray, according to the manufacturer’s instructions. Data were summarized, normalized, and analysed using the GeneSpring software V13.0 (Agilent, USA). The data were deposited in NCBI’s Gene Expression Omnibus under the accession number (GSE161361).
Human bone marrow mesenchymal stem cells (hBMSCs) (HUXMA-01001) were obtained from Cyagen Biosciences (Soochow, China). hBMSCs were cultured in Alpha Modified Eagle’s Medium (α-MEM) (Thermo Fisher, USA) containing 10% FBS (Thermo Fisher, USA), 100 U/mL penicillin (Thermo Fisher, USA), and 100 mg/mL streptomycin (Thermo Fisher, USA) in the incubator with 5% CO2 and 95% humidity at 37 °C. Cells were trypsinized and passaged at 85% confluence using 0.25% trypsin-EDTA (Thermo Fisher, USA). To induce osteogenic differentiation, cells were seeded into six-well plates until they reached 70% confluence and were then cultured in osteogenic induction medium (OIM) (HUXMA-90021, Cyagen Biosciences, China) for 3 days. The medium was completely refreshed every 3 days. Cells were washed and collected after induction for 2–3 weeks. To induce adipogenic differentiation, cells were seeded into six-well plates until they reached 100% confluence and were then cultured in adipogenic induction medium (AIM) Solution A (HUXMA-90031, Cyagen Biosciences, China) for 3 days. The medium was completely replaced with AIM Solution B (Cyagen Biosciences, China) for 24 h. After three cycles of Solution A/Solution B, cells were incubated in Solution B for an additional 5 days.
Small interfering RNA (siRNA) against hsa_circ_0006859 (si-circ), siRNA against ROCK1 (si-ROCK1), siRNA negative control siRNA (si-NC), miR-431-5p mimics (miR-431 M), miR-431-5p inhibitor (miR-431I), miRNA negative control (miR-NC), pcDNA 3.1 vector containing hsa_circ_0006859 (OE-circ), pcDNA 3.1 vector containing ROCK1 without 3′-UTR (OE-ROCK1), and empty vector control (EV) were custom-designed and synthesized by Beierbo (Beierbo, China). Briefly, hBMSCs were seeded into a 6-well plate (1 × 106 cells/well) and were cultured until they reached 85% confluence. Cells were then transfected with different treatments using Lipofectamine 2000 reagent (Thermo Fisher, USA) according to the manufacturer’s instructions.
RNA extraction and quantitative real-time PCR (qRT-PCR)
Total RNA was extracted from exosomes or cells using TRIzol reagent (Thermofisher, USA). The purity and concentration of RNA was spectrophotometrically analysed using a Nanodrop One (Thermo Fisher, USA). To determine expression of circRNAs, RNase R was added to digest linear transcripts. cDNA was transcribed using the PrimeScript RT reagent kit (TaKaRa, Japan). To determine expression of miRNAs, cDNA was transcribed using the TaqMan MicroRNA Reverse Transcription Kit (Thermo Fisher, USA). To determine expression of mRNAs, cDNA was transcribed using the First Strand cDNA Synthesis Kit (Beyotime, China). PCR was performed on an ABI 7500 Real-Time PCR system (Thermo Fisher, USA) using the TB Green Premix Ex Taq II (Taraka, Japan). Cel-miR-39 was the reference gene for circRNA in the exosome isolates, U6 was the reference gene for miRNA in the cell, and GAPDH was the reference gene for mRNA or circRNA in the cell. The primers used in the study were as follows: hsa_circ_0006859, Forward: 5′-ACTGTCCCGTACTCGCTGAT-3′, Reverse: 5′-CGGTTCTCATGCACACTGAC-3′; miR-431-5p, Forward: 5′-ACGCGTGTCTTGCAGGCCGT-3′, Reverse: 5′-ATCCAGTGCAGGGTCCGAGG-3′; ROCK1, Forward: 5′-ACCTGTAACCCAAGGAGATGTG-3′, Reverse: 5′-CACAATTGGCAGGAAAGTGG-3′; GAPDH, Forward: 5′-TGACCACAGTCCATGCCATCAC-3′; Reverse: 5′-GCCTGCTTCACCACCTTCTTGA-3′. U6, Forward: 5′-CTCGCTTCGGCAGCACATATACTA-3′; Reverse: 5′-ACGAATTTGCGTGTCATCCTTGC-3′. Each experiment was repeated independently at least three times.
Alkaline phosphatase (ALP) activity assay
ALP activity was detected using the Alkaline Phosphatase Assay Kit (Jiancheng, China) according to standard operation protocol. Absorbance was determined at 520 nm by BioTek Synergy 2 (BioTek, USA). Protein concentrations were analysed using the BCA Protein Assay Kit (Beyotime, China). The final value was calculated after normalization to protein content. Each experiment was repeated independently at least three times.
Alizarin red staining (ARS) assay
Cells were fixed in standard solution (PBS+4% paraformaldehyde) at room temperature and were stained with 0.1% Alizarin red (pH = 4.2, Solarbio, China) according to the manufacturer’s procedure. Cells were washed with PBS three times to remove unbound dye. The cells were observed and imaged under the IX71 microscope (Olympus, Japan). To quantify mineralization, the bound dye was dissolved in the eluent containing 10% cetylpyridinium chloride and then transferred to a 96-well plate. Absorbance was measured at 570 nm using a BioTek Synergy 2 (BioTek, USA). Each experiment was independently repeated at least three times.
Oil Red O staining
Cells were washed twice with PBS and then fixed in standard solution for 30 min at room temperature. Cells were then incubated with Oil Red O solution (Jiancheng, China) for 1 h at room temperature. Then, stained cells were washed with PBS three times and visualized under an IX71 microscope (Olympus, Japan).
Bioinformatics analysis and dual-luciferase assay
The bioinformatics database Circular RNA Interactome was used to search for potential miRNAs that may be targeted by hsa_circ_0006859 . Seventy-one miRNAs were predicted to be the targets of hsa_circ_0006859. As miR-431-5p was found to be correlated with adipogenic differentiation of hBMSCs , miR-431-5p was selected as one of the possible targets of hsa_circ_0006859. The two predicted binding sites on miR-431-5p targeted by hsa_circ_0006859 were at position 1668–1674 (#1) and position 401–407 (#2). Target genes of miR-431-5p were predicted by two independent databases, TargetScan  and DIANA Tools . Though IRS2 was proven to be one of the direct target genes of miR-431-5p , genes directly targeted by miR-431-5p were screened since one miRNA may target many genes. Target genes predicted by both computational algorithms were selected, and a literature review was performed to further select genes associated with osteoporosis. ROCK1 was selected as one of the miR-431-5p target genes. The prediction was further verified in vitro by luciferase assay. Two hsa_circ_0006859 fragments and the ROCK1 3′-UTR possessing wild type (Wt) or mutant (Mut) seed region were synthesized and cloned into psiCHECK-2 vectors (Promega, USA), termed circ-#1-Wt, circ-#1-Mut, circ-#2-Wt, circ-#2-Mut, ROCK1-Wt, and ROCK1-Mut. After hBMSCs reached 85% confluence, these constructed reporter plasmids were co-transfected with miR-431-5p mimics using Lipofectamine 2000 reagent (Thermo Fisher, USA) according to the manufacturer’s protocol. After transfection for 24 h, the fluorescence in each well was detected using the Dual Luciferase Reporter Assay System (Thermo Fisher, USA).
RNA pull-down assay
Biotinylated miR-NC, miR-431-5p-Wt, and miR-431-5p-Mut were transfected into hBMSCs with Lipofectamine 2000 reagent (Thermo Fisher, USA). After sonication, cell lysates were collected and incubated with Dynabeads M-280 Streptavidin (Thermo Fisher, USA) according to standard operation. Bound RNAs were purified using TRIzol reagent (Thermo Fisher, USA), and hsa_circ_0006859 was detected by qRT-PCR.
Total protein was extracted from hBMSCs after different treatments using RIPA lysis buffer (Beyotime, China) supplemented with protease inhibitors (Thermo Fisher, USA), and protein concentration was quantified using the BCA Protein Assay Kit (Beyotime, China). SDS-PAGE electrophoresis was performed to separate proteins. The gel was then transferred onto a PVDF membrane (Thermo Fisher, USA). The membrane was blocked in 5% non-fat milk for 1 h at room temperature and was subsequently incubated with specific primary antibody at 4 °C overnight. The membrane was rinsed three times with TBST and then incubated with specific secondary antibody for 1 h at room temperature. The membrane was washed, incubated with an enhanced chemiluminescence reagent (Thermo Fisher, USA), and observed on the ChemiDoc Touch Imaging System (Bio-Rad, USA). All experiments were independently repeated three times.
Each experiment was repeated three times. All data are expressed as the mean ± SEM. Statistical analysis was performed using GraphPad Prism 7.0. Student’s t test was used to test the difference between two groups. The one-way analysis of variance (ANOVA) was used to test the difference between three or more groups. Receiver operating characteristic (ROC) analysis was performed to evaluate the diagnostic ability of gene expression for disease prediction. Differences were considered statistically significant when P < 0.05.
CircRNAs in exosome are differentially expressed in osteoporosis patients
The circRNA expression profiles of exosomes isolated from the serum of osteoporosis patients (N = 3) and matched healthy controls (N = 3) were examined by circRNA microarray. To select potential differentially expressed circRNAs, the selection criteria were as follows: (1) circRNA was detected among all samples, (2) fold change > 4 or < 0.25, and (3) P < 0.05 according to the t test. The results showed that circRNAs in exosomes were differentially expressed between the two groups (Fig. 1a). Five hundred eighty-nine circRNAs were differentially expressed between the two groups, of which 376 circRNAs were upregulated and 213 were downregulated in osteoporosis patients. The five most upregulated or downregulated circRNAs in osteoporosis patients are listed in Supplementary Table 1. Though hsa_circ_006859 was the second most upregulated circRNA, hsa_circ_006859 was selected for further study due to hsa_circ_0006859 being found to exist in human serum , cells , and tissues , and it was upregulated in colorectal cancer tissues compared to the adjacent non-tumour tissues , while there are no reports concerning the role of hsa_circ_0047341. Thus, we chose to examine hsa_circ_006859. In total, these results indicated that circRNAs in exosome are differentially expressed in osteoporosis patients, and hsa_circ_006859 was selected for further experiments.
hsa_circ_0006859 in exosomes is upregulated in osteoporosis patients
The demographic characteristics of participants included in this study are shown in Table 1. A total of 168 postmenopausal women from the Third Affiliated Hospital of Soochow University were recruited, including 60 participants with normal bone mass as healthy controls, 50 participants with osteopenia, and 58 participants with osteoporosis. No significant difference was observed in age, body mass index (BMI), PTH, 25(OH) D, calcium, or inorganic phosphorus among the three groups. BMD, T-score, Z-score, β-CTx, tP1NP, and osteocalcin were significantly different among the osteoporosis, osteopenia, and healthy groups, indicating that participants with osteopenia or osteoporosis exhibited worsening bone parameters compared to the control group, especially the osteoporosis group (Table 1). The 3 osteoporosis patients and 3 healthy controls for circRNA microarray were not included. To determine whether hsa_circ_0006859 is involved in the development of osteoporosis, expression levels of hsa_circ_0006859 in exosomes were detected using qRT-PCR. Expression levels of hsa_circ_0006859 in exosomes were significantly higher in osteoporosis and osteopenia patients compared to non-osteoporotic controls (p < 0.0001, Fig. 1b). Expression levels of hsa_circ_0006859 were also significantly higher in osteoporosis patients than in osteopenia patients (p < 0.0001, Fig. 1b). We next analysed the relationship between expression levels of hsa_circ_0006859 in exosome and clinical parameters. Expression levels of hsa_circ_0006859 in exosome showed a strong negative correlation with L1-L4 spine T-score (Fig. 1c) and L1-L4 spine Z-score (Fig. 1d) across the entire cohort. To assess the diagnostic value of exosomal hsa_circ_0006859 to discriminate osteoporosis or osteopenia patients from healthy controls, ROC curve analysis was conducted. Hsa_circ_0006859 differentiated osteopenia patients from healthy controls with an AUC of 0.913 (95% CI 0.8617–0.9643, p < 0.0001), 94.00% sensitivity, and 76.67% specificity (Fig. 1e). Hsa_circ_0006859 also significantly differentiated osteoporosis patients from healthy controls with an AUC = 0.8974 (95% CI 0.8248–0.97, p < 0.0001), 93.1% sensitivity, and 93.33% specificity (Fig. 1f). Furthermore, hsa_circ_0006859 significantly differentiated osteoporosis patients from osteopenia patients with an AUC = 0.8873 (95% CI 0.8209–0.9536, p < 0.0001), 75.0% sensitivity, and 93.3% specificity (Fig. 1g). Expression levels of hsa_circ_0006859 in exosomes extracted from the serum of OPO patients after anti-osteoporotic treatment for 9 months were also significantly downregulated compared to those at diagnosis (Fig. 1h). These results suggested that hsa_circ_0006859 in exosomes was upregulated in patients with osteopenia or osteoporosis, and it may distinguish these patients from healthy controls, suggesting its potential for the diagnosis of osteoporosis and its critical role in the pathophysiology of osteoporosis.
hsa_circ_0006859 suppresses osteoblastic differentiation
To investigate the function of hsa_circ_0006859 in osteoporosis, expression of hsa_circ_0006859 during osteoblastic differentiation of hBMSCs induced by OIM was examined by qRT-PCR at indicated days (0, 7, 10, and 14 days). Hsa_circ_0006859 in cell was gradually decreased during osteoblastic differentiation in hBMSCs (Fig. 2a). The OE-circ or si-circ was transfected into hBMSCs to overexpress or silence hsa_circ_0006859, respectively. Expression of hsa_circ_0006859 was significantly increased by approximately 8-fold after OE-circ transfection and was decreased about half-fold after si-circ transfection, indicating that the cell transfection was successful (Fig. 2b). Osteocalcin (OCN) and alkaline phosphatase (ALP) are commonly used osteoblast differentiation markers. Relative mRNA expression levels of OCN (Fig. 2c) and ALP (Fig. 2d) were significantly decreased in hBMSCs transfected with OE-circ, while relative mRNA expression levels of OCN (Fig. 2c) and ALP (Fig. 2d) were significantly increased in hBMSCs transfected with si-circ after induction of osteoblastic differentiation for 10 days. Hsa_circ_0006859 overexpression significantly decreased protein levels of OCN and ALP in hBMSCs, while hsa_circ_0006859 inhibition had the opposite effect, which was illustrated by western blot analysis after hBMSCs were cultured in OIM for 10 days (Fig. 2e). In addition, ALP activity and ARS staining assays were performed to evaluate the effects of hsa_circ_0006859 on the osteogenic process. ALP activity in hBMSCs transfected with OE-circ was significantly decreased, while ALP activity was significantly enhanced in hBMSCs transfected with si-circ (Fig. 2f). Calcified nodules were markedly decreased following hsa_circ_0006859 overexpression, while calcified nodules were significantly increased following hsa_circ_0006859 suppression (Fig. 2g). These results imply that hsa_circ_0006859 suppresses osteoblastic differentiation of hBMSCs.
Hsa_circ_0006859 promotes adipogenic differentiation
We next focused on the biological function of hsa_circ_0006859 in adipogenic differentiation. As shown in Fig. 3a, hsa_circ_0006859 gradually increased during adipogenic differentiation in hBMSCs. The mRNA expression levels of adipogenic markers C/EBPα and PPARγ were examined by qRT-PCR after adipogenic differentiation for 12 days. hBMSCs transfected with OE-circ exhibited higher mRNA expression levels of C/EBPα (Fig. 3b) and PPARγ (Fig. 3c) compared to EV, while hBMSCs transfected with si-circ presented lower mRNA expression levels of C/EBPα (Fig. 3b) and PPARγ (Fig. 3c) compared to si-NC. Protein levels of C/EBPα and PPARγ were examined by western blot analysis, and their changes were similar to the mRNA changes (Fig. 3d). Furthermore, overexpression of hsa_circ_0006859 by OE-circ enhanced the formation of lipid droplets, while suppression of hsa_circ_0006859 by si-circ attenuated adipogenic differentiation of hBMSCs as indicated by Oil red O staining (Fig. 3e). Therefore, these results indicate that hsa_circ_0006859 enhances adipogenic differentiation of hBMSCs.
hsa_circ_0006859 directly binds to miR-431-5p
Given that circRNAs act as miRNA sponges to negatively regulate miRNA activity , the regulatory function of hsa_circ_0006859 at the post-transcriptional level was explored. Based on the bioinformatics analysis and literature review, miR-431-5p may directly bind to hsa_circ_0006859. Two potential miR-431-5p binding sites were identified in the hsa_circ_0006859 sequence at position 1668–1674 and position 401–407 (Fig. 4a). The minimum free energy values of the two hybrids were − 18.9 kcal/mol and − 15.6 kcal/mol (Fig. 4a). Direct interactions between hsa_circ_0006859 and miR-431-5p were verified by dual luciferase reporter assays. circ-#1-Wt, circ-#1-Mut, circ-#2-Wt, or circ-#2-Mut was co-transfected with miR-431 M into hBMSCs. The luciferase activity of circ-#1-Wt (Fig. 4b) or circ-#2-Wt (Fig. 4c) was significantly inhibited by miR-431 M, while the luciferase activity of circ-#1-Mut (Fig. 4b) or circ-#2-Mut (Fig. 4c) was not significantly changed by miR-431 M (Fig. 4c). Moreover, expression levels of miR-431-5p in hBMSCs transfected with hsa_circ_0006859 overexpression vector or hsa_circ_0006859 siRNA were determined by qRT-PCR. The expression of miR-431-5p was significantly downregulated by OE-circ, while it was significantly upregulated by si-circ (Fig. 4d). The RNA pull-down assay was performed to further verify the interaction between hsa_circ_0006859 and miR-431-5p. Relative hsa_circ_0006859 enrichment was significantly higher in the Bio-miR-431-5p-Wt group than in the Bio-miR-431-5p-Mut group or the Bio-miR-NC group (Fig. 4e). Next, expression of miR-431-5p in exosomes was detected using qRT-PCR. Expression of miR-431-5p in exosomes was significantly lower in osteoporosis patients compared to osteopenia patients and non-osteoporotic controls (Fig. 4f). The correlation between hsa_circ_0006859 expression and miR-431-5p expression in the exosomes of clinical samples was further illustrated using Pearson’s correlation analysis. There was an inverse correlation between hsa_circ_0006859 expression and miR-431-5p expression in exosomes (Fig. 4g). These findings support the conclusion that hsa_circ_0006859 directly sponges miR-431-5p.
hsa_circ_0006859 acts as a competing endogenous RNA (ceRNA) to promote ROCK1 expression
In this experiment, two biological algorithms, TargetScan and DIANA, were used to predict potential target genes that may share the binding sites of miR-431-5p matched in hsa_circ_0006859. ROCK1 was identified as a potential target of miR-431-5p (Fig. 5a). To identify whether miR-431-5p directly targets ROCK1 3′-UTR, Wt-ROCK1 or Mut-ROCK1 was co-transfected with miR-431-5p mimics into hBMSCs. The dual-luciferase assay showed that the luciferase activity of ROCK1-Wt was significantly inhibited by miR-431 M, while ROCK1-Mut relieved the effect induced by miR-431 M (Fig. 5b). Moreover, mRNA expression of ROCK1 and protein expression of ROCK1 in hBMSCs transfected with miR-431 M or miR-431I were determined by qRT-PCR and western blot analysis. The qRT-PCR results demonstrated that mRNA expression of ROCK1 was significantly upregulated by miR-431I and significantly downregulated by miR-431 M (Fig. 5c). Western blot results showed that protein expression of ROCK1 was significantly downregulated by miR-431 M and significantly upregulated by miR-431I (Fig. 5d). These findings support the conclusion that ROCK1 is a direct target gene of miR-431-5p. Moreover, the role of hsa_circ_0006859 in ROCK1 regulation was investigated. The dual-luciferase assay demonstrated that the luciferase activity of ROCK1-Wt was decreased by hsa_circ_0006859 siRNA, while the reduced luciferase activity was rescued by miR-431I (Fig. 5e). However, the luciferase activity of ROCK1-Mut was barely changed by si-circ or si-circ+miR-431I (Fig. 5e). mRNA expression levels of ROCK1 (Fig. 5f) and protein expression levels of ROCK1 (Fig. 5g) were significantly decreased by hsa_circ_0006859 knockdown and were rescued by miR-431I. Taken together, these results suggest that hsa_circ_0006859 promotes ROCK1 expression by competitively binding to miR-431-5p in hBMSCs.
hsa_circ_0006859 suppresses osteogenesis and promotes adipogenesis by sponging miR-431-5p to upregulate ROCK1
To investigate the function of miR-431-5p in osteoporosis, miR-431 M or miR-431I was transfected into hBMSCs to overexpress or silence miR-431-5p, respectively. Expression of miR-431-5p was significantly increased by approximately 8-fold after miR-431 M transfection and by about half-fold decreased after miR-431I transfection, indicating that the cell transfection was successful (Fig. 6a). ALP activity and ARS staining assays were performed to evaluate the effects of miR-431-5p on the osteogenic process. ALP activity in hBMSCs transfected with miR-431 M was significantly enhanced, while the ALP activity was significantly suppressed in hBMSCs transfected with miR-431I (Fig. 6b). Calcified nodules were increased following miR-431-5p overexpression and decreased following miR-431-5p suppression (Fig. 6c). We next focused on the biological function of miR-431-5p in adipogenic differentiation. Overexpression of miR-431-5p by miR-431 M weakened the formation of lipid droplets, while suppression of miR-431-5p by miR-431I strengthened the adipogenic differentiation of hBMSCs as indicated by the Oil red O staining (Fig. 6d). Therefore, these results indicate that miR-431-5p promotes osteoblastic differentiation and inhibits adipogenic differentiation of hBMSCs.
Next, we investigated the functional role of ROCK1 in the regulation of osteoblastic differentiation and adipogenic differentiation. OE-ROCK1 was transfected into hBMSCs to overexpress ROCK1, while si-ROCK1 was transfected into hBMSCs to inhibit ROCK1. qRT-PCR analysis showed that OE-ROCK1 significantly enhanced ROCK1 mRNA expression, while si-ROCK1 #2 and #3 significantly decreased ROCK1 mRNA expression in hBMSCs (Fig. 6e). Western blot analysis demonstrated that ROCK1 protein expression was significantly enhanced by ROCK1 OE transfection, while ROCK1 protein expression was significantly decreased by ROCK1 si-ROCK1 #2 and #3 in hBMSCs (Fig. 6f). As si-ROCK1 #3 was very effective in knocking down ROCK1, both at the mRNA and protein levels, it was selected in the following experiments. ROCK1 overexpression significantly decreased the mineralized nodule formation, while ROCK1 knockdown significantly enhanced mineralized nodule formation compared to their respective controls (Fig. 6g). Furthermore, ROCK1 overexpression significantly enhanced the formation of lipid droplets, while ROCK1 knockdown significantly decreased the formation of lipid droplets compared to their respective controls (Fig. 6h). These results indicate that ROCK1 suppresses osteoblastic differentiation and promotes adipogenic differentiation of hBMSCs.
Because hsa_circ_0006859 was found to directly bind with miR-431-5p and ROCK1 was directly targeted by miR-431-5p, we speculated that the regulatory effects of hsa_circ_0006859 on hBMSCs osteogenesis and adipogenesis were mediated by the miR-431-5p/ROCK1 axis. ARS staining revealed that the reduced osteoblastic differentiation caused by hsa_circ_0006859 overexpression was partially reversed by ROCK1 knockdown (Fig. 6i). Oil red O staining revealed that the enhanced adipogenic differentiation caused by hsa_circ_0006859 overexpression was partially reversed by ROCK1 knockdown (Fig. 6j). Collectively, hsa_circ_0006859 suppressed osteogenesis and promoted adipogenesis via sponging miR-431-5p to upregulate ROCK1.
Osteoporosis is a skeletal disorder characterized by impaired bone density and quality, which increase the risk of bone fractures, especially in postmenopausal women . Many studies have shown that circRNAs are involved in the progression of a variety of pathological conditions and are likely to be new potential clinical biomarkers or therapeutic targets in human diseases [5, 29]. In the present study, we found that circRNAs in exosome were differentially expressed in osteoporosis patients compared to healthy controls. Hsa_circ_0006859 in exosomes was significantly upregulated in osteoporosis patients and had a good discrimination to distinguish osteoporosis patients from healthy controls. Furthermore, hsa_circ_0006859 inhibited osteoblastic differentiation and promoted adipogenic differentiation of hBMSCs by sponging miR-431-5p to upregulate ROCK1 in vitro. Taken together, our results demonstrated that hsa_circ_0006859 may be a potential biomarker for detection and disease progression of osteoporosis and may represent a promising therapeutic target for the treatment of osteoporosis in postmenopausal women.
CircRNAs are abundant in exosomes, and more than one thousand exosomal circRNAs have been identified in human serum . CircRNAs are ideal diagnostic biomarkers in human diseases due to their universality (prevailing expression), conservatism (evolutionally conserved), stability (highly resistant to RNase R), and specificity (tissue- and stage-specific) . For example, exosomal hsa_circ_0004771 was upregulated in the serum of colorectal cancer patients and was downregulated in the serum of postoperative colorectal cancer patients. Exosomal hsa_circ_0004771 was able to differentiate colorectal cancer patients from healthy controls with high sensitivity and specificity . Several serum exosomal circRNAs, such as hsa_circ_0082333, were able to distinguish patients with colon cancer from healthy controls . Exosomal hsa_circ_0130810 was decreased in gastric cancer patients, and its reduction was correlated with shorter overall survival time, indicating that exosomal hsa_circ_0130810 may serve as a novel biomarker for gastric cancer . Currently, the preferred technology for diagnosing osteoporosis is the DXA (dual-energy X-ray absorptiometry) scan, which measures BMD in the spine, hip, or forearm and predicts skeletal fracture risk . Although DXA is valuable, it can only reveal already established changes in bone structure, which may take a long time to become detectable and can only provide limited information on disease progression and treatment efficacy . To date, several circRNAs have been proposed as potentially valuable biomarkers for osteoporosis. Hsa_circ_0001275 in peripheral blood mononuclear cells was higher in osteoporosis patients when compared to healthy controls, and its expression was negatively correlated with T-score. ROC curve analysis showed that hsa_circ_0001275 efficiently differentiated osteoporosis patients from healthy controls . Furthermore, expression of serum hsa_circ_0002060 was upregulated in osteoporosis patients compared to control groups, and its expression was associated with low BMD in osteoporosis patients. ROC analysis indicated that hsa_circ_0002060 shows potential diagnostic value for osteoporosis . In the present study, exosomal hsa_circ_0006859 was upregulated in postmenopausal women with osteoporosis. Its expression in exosomes was negatively correlated with L1-L4 spine T-score and Z-score. Furthermore, hsa_circ_0006859 significantly differentiated osteoporosis patients from healthy controls. Taken together, our results suggest the clinical potential of exosomal hsa_circ_0006859 in the diagnosis of osteoporosis. To our knowledge, this is the first study in which exosomal hsa_circ_0006859 in serum was successfully identified as a potential biomarker for postmenopausal osteoporosis.
An increasing number of studies have demonstrated that circRNAs play critical roles in multiple diseases through different mechanisms . Some circRNAs interact with miRNAs, which releases miRNA-mediated mRNA inhibition in osteoporosis. Plasma circRNA_0016624 is downregulated in osteoporosis. It was found to prevent osteoporosis and promote BMP2 expression by sponging miR-98 . Hsa_circ_0006393 promotes osteogenesis by sponging miR-145-5p and upregulating FOXO1 in osteoporosis . Hsa_circ_0028313 knockdown suppressed osteoclast differentiation and prevented bone loss via the miR-195a/CSF1 axis . Hsa_circ_0076906 promoted osteogenesis and alleviated osteoporosis by acting as a sponge for miR-1305 and upregulating osteoglycin . In the present study, hsa_circ_0006859 was found to suppress osteogenesis and promote adipogenesis by sponging miR-431-5p, which released the inhibitory effect of miR-431-5p on its direct target, ROCK1. This is the first study demonstrating the important role of the hsa_circ_0006859/miR-431-5p/ROCK1 axis during osteoporosis.
Multiple RNA populations have been identified in exosomes, including mRNAs and many types of non-coding RNAs, such as circRNAs, miRNAs, and lncRNAs . These miRNAs in exosomes are categorized into circulating miRNAs, which are involved in many aspects of osteoporosis. Some circulating miRNAs may be applied at the prediction stage for osteoporosis. Upregulation of serum miR-21-5p and miR-23a-3p and downregulation of serum miR-125b-5p were correlated with bone loss and muscle loss in postmenopausal women [36, 37]. Circulating miR-194-5p was upregulated in the blood of postmenopausal patients with osteoporosis and was negatively correlated with T-scores . Some circulating miRNAs may be applied at the diagnostic stage for osteoporosis. Plasma miR-19b-3p was significantly lower in osteoporosis patients and differentiated osteoporosis patients with vertebral fracture from healthy controls with relative high accuracy . Serum miR-21-5p clearly distinguished healthy controls, osteopenia and osteoporosis patients with relative high sensitivity and specificity in diagnosis [37, 40]. Serum miR-30b-5p was significantly downregulated in osteopenia or osteoporosis patients and was positively correlated with BMD . Some circulating miRNAs may be applied at the treatment stage for osteoporosis. Serum miR-29b-3p and miR-324-p were reduced in patients undergoing anti-resorptive therapy . Serum miR-33-3p and miR-133a-3p were decreased during teriparatide treatment . In our study, miR-431-5p was also found in osteoporosis serum exosomes. Furthermore, expression levels of miR-431-5p were negatively correlated with expression levels of hsa_circ_0006859 in exosomes. Our results indicate the important role of circulating miR-431-5p in osteoporosis.
Increasing evidence has suggested that miRNAs play a pivotal role in bone homeostasis. During ageing, hBMSCs reduce the ability to differentiate into osteoblasts but increase the ability to differentiate into adipocytes, which causes bone loss, resulting in osteoporosis. For example, miR-199a-5p targeting TET2 , miR-218-5p targeting COL1A1 , and miR-291a-3p targeting DKK1  promote osteogenic differentiation of hBMSCs and promote bone formation to prevent bone loss, while miR-23a-3p targeting PGC-1α , miR-31-5p targeting Frizzled-3 , and miR-128 targeting SIRT6  inhibit osteoblast differentiation. Some miRNAs are involved in positive or negative regulation of hBMSC adipogenic differentiation. miR-7b-5p targeting IRS2  suppresses adipogenic differentiation of hBMSCs, while miR-214-5p targeting COL4A1  may promote adipogenic differentiation of hBMSCs. In the present study, miR-431-5p was proven to promote osteogenic differentiation and inhibit adipogenic differentiation of hBMSCs by targeting ROCK1, indicating the potential role of miR-431-5p in controlling the osteoblastic and adipogenic fate. Furthermore, reduced osteoblastic differentiation and enhanced adipogenic differentiation caused by hsa_circ_0006859 overexpression was partially reversed by ROCK1 knockdown, indicating that hsa_circ_0006859 suppresses osteogenesis and promotes adipogenesis via sponging miR-431-5p to upregulate ROCK1.
In conclusion, exosomal hsa_circ_0006859 is upregulated in osteoporosis patients. Hsa_circ_0006859 inhibits osteoblastic differentiation and promotes adipogenic differentiation of hBMSCs by sponging miR-431-5p to upregulate ROCK1. Our findings suggest a role of hsa_circ_0006859 as a novel circulating biomarker for osteoporosis, and the mechanism of hsa_circ_0006859 induces imbalance between osteogenesis and adipogenesis, representing a novel target for osteoporosis treatment.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Quantitative real-time PCR
Alizarin red staining
Human bone marrow mesenchymal stem cells
Competing endogenous RNA
Bone mineral density
Dual energy X-ray absorptiometry
β-Isomerized C-terminal telopeptides
N-Terminal procollagen of type l collagen
Alpha modified eagle’s medium
Osteogenic induction medium
Adipogenic induction medium
Small interfering RNA
siRNA against hsa_circ_0006859
siRNA against ROCK1
siRNA negative control siRNA
- miR-431 M:
miRNA negative control
pcDNA 3.1 vector containing hsa_circ_0006859
pcDNA 3.1 vector containing ROCK1 without 3′-UTR
Empty vector control
Alizarin red staining
Receiver operating characteristic
Black DM, Rosen CJ. Postmenopausal osteoporosis. N Engl J Med. 2016;374(3):254–62.
Yang T-L, Shen H, Liu A, et al. A road map for understanding molecular and genetic determinants of osteoporosis. Nature Reviews Endocrinol. 2019;16(2):91–103.
Tu KN, Lie JD, Wan CKV, et al. Osteoporosis: a review of treatment options. P & T : a Peer-Reviewed J Formulary Management 2018;43(2):92–104.
Memczak S, Jens M, Elefsinioti A, et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013;495(7441):333–8.
Kristensen LS, Andersen MS, Stagsted LVW, et al. The biogenesis, biology and characterization of circular RNAs. Nat Rev Genet. 2019;20(11):675–91.
Li Y, Zheng Q, Bao C, et al. Circular RNA is enriched and stable in exosomes: a promising biomarker for cancer diagnosis [in eng]. Cell Res. 2015;25(8):981–4.
Fanale D, Taverna S, Russo A, et al. Circular RNA in Exosomes. Adv Exp Med Biol. 2018;1087:109–17.
Gebert LFR, MacRae IJ. Regulation of microRNA function in animals. Nat Rev Mol Cell Biol. 2019;20(1):21–37.
Aufiero S, Reckman YJ, Pinto YM, et al. Circular RNAs open a new chapter in cardiovascular biology. Nat Rev Cardiol. 2019;16(8):503–14.
Huang Y, Xie J, Li E. Comprehensive circular RNA profiling reveals circ_0002060 as a potential diagnostic biomarkers for osteoporosis. J Cell Biochem. 2019;120(9):15688–94.
Wang XB, Li PB, Guo SF, et al. circRNA_0006393 promotes osteogenesis in glucocorticoidinduced osteoporosis by sponging miR1455p and upregulating FOXO1. Mol Med Rep. 2019;20(3):2851–8.
Yu L, Liu Y. circRNA_0016624 could sponge miR-98 to regulate BMP2 expression in postmenopausal osteoporosis. Biochem Biophys Res Commun. 2019;516(2):546–50.
Wen J, Guan Z, Yu B, et al. Circular RNA hsa_circ_0076906 competes with OGN for miR-1305 biding site to alleviate the progression of osteoporosis. Int J Biochem Cell Biol. 2020;122:105719.
Della Bella E, Menzel U, Basoli V, et al. Differential regulation of circRNA, miRNA, and piRNA during early osteogenic and chondrogenic differentiation of human mesenchymal stromal cells. Cells. 2020;9(2):398.
Yuan W, Peng S, Wang J, et al. Identification and characterization of circRNAs as competing endogenous RNAs for miRNA-mRNA in colorectal cancer. PeerJ. 2019;7:e7602.
Zhuang C, Huang X, Yu J, et al. Circular RNA hsa_circ_0075828 promotes bladder cancer cell proliferation through activation of CREB1. BMB Rep. 2020;53(2):82–7.
Vilades D, Martinez-Camblor P, Ferrero-Gregori A, et al. Plasma circular RNA hsa_circ_0001445 and coronary artery disease: performance as a biomarker. FASEB J. 2020;34(3):4403–14.
Hang D, Zhou J, Qin N, et al. A novel plasma circular RNA circFARSA is a potential biomarker for non-small cell lung cancer. Cancer Med. 2018;7(6):2783–91.
Zhao Z, Li X, Jian D, et al. Hsa_circ_0054633 in peripheral blood can be used as a diagnostic biomarker of pre-diabetes and type 2 diabetes mellitus. Acta Diabetol. 2017;54(3):237–45.
Dudekula DB, Panda AC, Grammatikakis I, et al. CircInteractome: a web tool for exploring circular RNAs and their interacting proteins and microRNAs [in eng]. RNA Biol. 2016;13(1):34–42.
Wang Y, Yang L, Liu X, et al. miR-431 inhibits adipogenic differentiation of human bone marrow-derived mesenchymal stem cells via targeting insulin receptor substance 2. Stem Cell Res Therapy. 2018;9(1):231.
Lewis BP, Shih IH, Jones-Rhoades MW, et al. Prediction of mammalian microRNA targets. Cell. 2003;115(7):787–98.
Paraskevopoulou MD, Georgakilas G, Kostoulas N, et al. DIANA-microT web server v5.0: service integration into miRNA functional analysis workflows. Nucleic Acids Res. 2013;41(Web Server issue):W169–73.
Maass PG, Glažar P, Memczak S, et al. A map of human circular RNAs in clinically relevant tissues [in eng]. J Molecular Med. 2017;95(11):1179–89.
Della Bella E, Menzel U, Basoli V, et al. Differential regulation of circRNA, miRNA, and piRNA during early osteogenic and chondrogenic differentiation of human mesenchymal stromal cells. Cells 2020;9(2):398.
Yuan W, Peng S, Wang J, et al. Identification and characterization of circRNAs as competing endogenous RNAs for miRNA-mRNA in colorectal cancer [in eng]. PeerJ. 2019;7:e7602.
Hansen TB, Jensen TI, Clausen BH, et al. Natural RNA circles function as efficient microRNA sponges. Nature. 2013;495(7441):384–8.
Ensrud KE, Crandall CJ. Bisphosphonates for postmenopausal osteoporosis. JAMA. 2019. https://doi.org/10.1001/jama.2019.15781.
Zhang Z, Yang T, Xiao J. Circular RNAs: promising biomarkers for human diseases [in eng]. EBioMedicine. 2018;34:267–74.
Pan B, Qin J, Liu X, et al. Identification of serum exosomal hsa-circ-0004771 as a novel diagnostic biomarker of colorectal cancer. Front Genet. 2019;10:1096.
Tang W, Fu K, Sun H, et al. CircRNA microarray profiling identifies a novel circulating biomarker for detection of gastric cancer. Mol Cancer. 2018;17(1):137.
Choksi P, Jepsen KJ, Clines GA. The challenges of diagnosing osteoporosis and the limitations of currently available tools. Clinical Diabetes Endocrinol. 2018;4:12.
Bottani M, Banfi G, Lombardi G. Perspectives on miRNAs as epigenetic markers in osteoporosis and bone fracture risk: a step forward in personalized diagnosis. Front Genet. 2019;10:1044.
Zhao K, Zhao Q, Guo Z, et al. Hsa_Circ_0001275: a potential novel diagnostic biomarker for postmenopausal osteoporosis. Cell Physiol Biochem. 2018;46(6):2508–16.
Chen X, Ouyang Z, Shen Y, et al. CircRNA_28313/miR-195a/CSF1 axis modulates osteoclast differentiation to affect OVX-induced bone absorption in mice. RNA Biol. 2019;16(9):1249–62.
Chen Z, Bemben MG, Bemben DA. Bone and muscle specific circulating microRNAs in postmenopausal women based on osteoporosis and sarcopenia status. Bone. 2019;120:271–8.
Kelch S, Balmayor ER, Seeliger C, et al. miRNAs in bone tissue correlate to bone mineral density and circulating miRNAs are gender independent in osteoporotic patients. Sci Rep. 2017;7(1):15861.
Meng J, Zhang D, Pan N, et al. Identification of miR-194-5p as a potential biomarker for postmenopausal osteoporosis. PeerJ. 2015;3:e971.
Sun M, Hu L, Wang S, et al. Circulating microRNA-19b identified from osteoporotic vertebral compression fracture patients increases bone formation. J Bone Mineral Res. 2019;35(2):306–16.
Yavropoulou MP, Anastasilakis AD, Makras P, et al. Expression of microRNAs that regulate bone turnover in the serum of postmenopausal women with low bone mass and vertebral fractures. Eur J Endocrinol. 2017;176(2):169–76.
Chen J, Li K, Pang Q, et al. Identification of suitable reference gene and biomarkers of serum miRNAs for osteoporosis. Sci Rep. 2016;6:36347.
Feichtinger X, Muschitz C, Heimel P, et al. Bone-related circulating microRNAs miR-29b-3p, miR-550a-3p, and miR-324-3p and their association to bone microstructure and histomorphometry. Sci Rep. 2018;8(1):4867.
Anastasilakis AD, Makras P, Pikilidou M, et al. Changes of circulating microRNAs in response to treatment with teriparatide or denosumab in postmenopausal osteoporosis. J Clin Endocrinol Metab. 2018;103(3):1206–13.
Qi X-B, Jia B, Wang W, et al. Role of miR-199a-5p in osteoblast differentiation by targeting TET2. Gene. 2020;726:144193.
Kou J, Zheng X, Guo J, et al. MicroRNA-218-5p relieves postmenopausal osteoporosis through promoting the osteoblast differentiation of bone marrow mesenchymal stem cells. J Cell Biochem. 2020;121(2):1216–26.
Li ZH, Hu H, Zhang XY, et al. MiR-291a-3p regulates the BMSCs differentiation via targeting DKK1 in dexamethasone-induced osteoporosis. Kaohsiung J Med Sci. 2020;36(1):35–42.
Dai Y, Zheng C, Li H. Inhibition of miR-23a-3p promotes osteoblast proliferation and differentiation. J Cell Biochem. 2019. https://doi.org/10.1002/jcb.29497.
Weilner S, Schraml E, Wieser M, et al. Secreted microvesicular miR-31 inhibits osteogenic differentiation of mesenchymal stem cells. Aging Cell. 2016;15(4):744–54.
Zhao J, Liu S, Zhang W, et al. MiR-128 inhibits the osteogenic differentiation in osteoporosis by down-regulating SIRT6 expression. Biosci Rep. 2019;39(9):BSR20191405.
Li QQ, Wei Q, Zhai XC, et al. MiRNA-7b-5p attenuates the progression of osteoporosis by inhibiting adipose differentiation of hMSCs via regulating IRS2. Eur Rev Med Pharmacol Sci. 2019;23(21):9207–14.
Qiu J, Huang G, Na N, et al. MicroRNA-214-5p/TGF-beta/Smad2 signaling alters adipogenic differentiation of bone marrow stem cells in postmenopausal osteoporosis. Mol Med Rep. 2018;17(5):6301–10.
This study was funded by the Changzhou High-Level Medical Talents Training Project (2016CZBJ006) and Major Science and Technology Project of Changzhou Health Commission (ZD202008).
Ethics approval and consent to participate
The study was approved by the Ethics Committee of the Third Affiliated Hospital of Soochow University (No.CZYY2020032). Each participant signed an informed consent document.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhi, F., Ding, Y., Wang, R. et al. Exosomal hsa_circ_0006859 is a potential biomarker for postmenopausal osteoporosis and enhances adipogenic versus osteogenic differentiation in human bone marrow mesenchymal stem cells by sponging miR-431-5p. Stem Cell Res Ther 12, 157 (2021). https://doi.org/10.1186/s13287-021-02214-y