1 Introduction

Nowadays, chronic psychological stress, a critical public health issue, is getting more and more attention. It is responsible for various systemic diseases, one of which is the disruption of bone metabolism. Several studies had shown that chronic stress and even depression can lead to osteoporosis(Calarge et al., 2014; Furlan et al., 2005; Yirmiya et al., 2006). Whereas, how does psychological stress signal the brain but affect peripheral organs? Former research suggested that psychological stress can stimulate the autonomic nerves or the central nervous system to disrupt the activation of sympathetic nerves(Kondo & Togari, 2003; Ma et al., 2023) and the endocrine homeostasis(Raison & Miller, 2003) respectively. They alters the levels of neurotransmitters(Chen & Wang, 2017; Foertsch et al., 2017; Guo et al., 2023; Khosla et al., 2018; Teong et al., 2017), neuropeptides(Baldock et al., 2014; Ng & Chin, 2021), hormones(Hasan et al., 2012; Henneicke et al., 2017), related biomolecules, and their corresponding receptors in vivo, so as to affect the activations of several signaling pathways. Therefore, it contributes to the increased osteoclast activity, decreased osteoblast activity, and promoted bone matrix degradation, ultimately leading to bone loss. However, treatment of antidepressant and molecules targeting the above-mentioned receptors can only partially alleviate the bone loss(Lam et al., 2022; Marenzana et al., 2007), which suggest that there may be other unanticipated molecular mechanisms involved in psychological stress-induced osteoporosis and efficient for the clinical treatment. Besides, these mentioned neurotransmitters or hormones and their related receptors have broad biological effects, and their effects in vivo are highly dose-dependent, which usually give rise to an unpredictable side-effect when using hormonotherapy or neuropeptide therapy for the patients with psychological stress-related osteoporosis(S. Kim et al., 2021; Marenzana et al., 2007; Teong et al., 2017). Therefore, it is full of necessity to further explore alternative mechanisms that implicated in psychological stress-related osteoporosis, since any of these new discoveries will be possible for attracting an ideal clinical therapy with higher biosafety and efficiency.

MicroRNAs (miRNAs) carried by extracellular vesicles (EVs) or proteins in vivo are highly involved in biocommunication between multiple organs. They usually have a high regulatory efficiency for targeting multiple transcripts, and a wide regulatory range for controlling the expression of genes coding for protein. Particularly, due to the protection of EVs, they are not only stable in the extracellular environment, such as blood(Chen et al., 2008), but also free to pass through the blood-brain barrier, indicating that miRNAs may work as the linkers between brain and bone. Consistently, a previous study reported that miR-483-5p, stored in brain-derived EVs, could promote bone-fat imbalance in Alzheimer’s disease(Liu et al., 2023). Take the consideration that several endogenic miRNAs in the brain(Olejniczak et al., 2018) and blood(Makrygianni & Chrousos, 2023) were reported to change their expression levels after chronic psychological stress, it’s deducible that psychological stress may interfere the bone metabolism through brain-derived miRNAs. Therefore, it is significative and valuable to explore the biofunctions of brain-derived miRNAs during the pathogenesis of psychological stress-related osteoporosis, which can also preliminarily discover the roles of brain-derived miRNAs during the central-peripheral linkage. What’s more, as miRNAs are ideal targets for precision therapy due to their effective biofunctions, easy synthesis and longer half-life, exploring their regulatory role in the pathogenesis of the psychological stress-related osteoporosis is crucial for building the theoretical foundation, based on which, a more efficient and much safer clinical treatment can be developed.

In this context, we primarily established the CUMS model mice, characterized its osteoporosis phenotype, and analyzed the activities of osteoblasts and osteoclasts. Subsequently, using high-throughput miRNA-seq analysis of bone, we identified the potential miRNA candidate that is responsible for bone loss under psychological stress, as well as, investigated the molecular mechanism of this miRNA. Our findings show that the abnormal osteoclast activity is much more responsible for the psychological stress induced osteoporosis. Moreover, miR-335-3p, which can inhibit Fos translation and its integration with Nftac1, is significantly reduced under psychological stress and sequentially enhances the activation of Nftac1 signaling pathway, thereby upregulating the osteoclast differentiation and maturation. These results imply that miR-335-3p plays an important role in the communication between central nervous system and peripheral bone tissue, in addition, provide a new perspective on neuroregulatory mechanisms of the brain that participated in bone metabolism when suffering psychological stress.

2 Results

2.1 CUMS induces behavioral and physiological changes

The mental state of mice were examined by body weight, behavioral and serum parameters. It is showed that the CUMS mice had significantly lower body weights (Figure 1B) compared to the control group. They had considerably lower distances (Figure 1C) and times (Figure 1D) within the central area of the open field but did not differ from the control group in terms of the total distance moved in the open field (Figure 1E). Besides, the tail-hanging experiment showed that the immobilization time of mice in the CUMS group to give up struggling was significantly increased (Figure 1F). At the same time, they had lower sucrose preference index (Figure 1G). In addition, serum corticosterone (Figure 1H) and norepinephrine (Figure 1I) levels were significantly elevated in CUMS mice. Taken together, these results all indicated that CUMS modeling caused the mice to develop anxiety and depression-like behaviors.

Figure 1.

2.2 CUMS leads to bone loss

Next, mouse femurs from different groups were dissected for micro-CT scanning (Figure 2A-D) and H&E staining (Figure 2E). Then quantitatively analyzed histomorphometry, bone mineral density, and bone microarchitecture to assess the bone remodeling influenced by the psychological stress. It was found that the femurs of CUMS mice had sparse cancellous bone and thinner cortical bone compared to control mice, indicating that psychological stress harmed both cancellous and cortical bone.

Figure 2.

2.3 CUMS promotes bone formation slightly but enhance osteoclast activity significantly

To further explore the cause of psychological stress induced bone loss, the osteogenic activity of mice were then examined in the CUMS group. Firstly, Masson staining visualized that more osteoid was present in cortical and cancellous bone in CUMS mice (Figure 3A). In addition, the distance between two fluorochromes calcein lines was found wider in CUMS mice, suggesting the mineral deposition rate in CUMS mice was faster than control ones (Figure 3D).Then RNA were isolated to examine osteoblast-related genes, and found that the expression of Osterix (Osx), Osteocalcin (Ocn), osteopontin (Opn), and dentin matrix protein 1 (Dmp1) were elevated at the mRNA level in CUMS mice, whereas no significant difference was observed in runt-related transcription factor 2 (Runx2) and alkaline phosphatase (Alp) (Figure 3E). Further, immunohistochemical staining of OSX and OCN displayed that the protein expression of OSX was higher in CUMS mice (Figure 3B), but no considerable difference could be found in OCN (Figure 3C). Taken together, these data suggested that psychological stress promoted osteoblast functions in CUMS mice.

Figure 3.

Since the enhanced osteogenesis could not explain the reduced bone mass under psychological stress, we further examined the osteoclastic activity of CUMS mice. Tartrate-resistant acid phosphatase (TRAP) staining showed a significant increase in the number and area of TRAP-positive osteoclasts under CUMS (Figure 3G-H). Then, RT-PCR analysis of mouse femurs revealed that several osteoclastic genes, such as Ca2, Mmp9, Nfatc1, and Acp5, were significantly increased in CUMS mice (Figure 3I, Figure 3— figure supplement 1). Consistently, immunohistochemistry analyses also showed an elevated MMP9 expression in CUMS mice (Figure 3J). Thus, these data demonstrate that osteoclastic activity is elevated in vivo under psychological stress, which may be a more direct and vital cause of the decline in bone mass.

2.4 miR-335-3p is apparently downregulated in CUMS mice

To search the key miRNAs involved in the abnormal bone metabolism correlated to chronic stress, femurs (dissected off the growth plate) were cellected from different groups for sRNA-seq, respectively. The heat map and Volcano plot analysis displayed that a number of miRNAs were significantly altered in CUMS mice (Figure 4A-B). The most significant variations including miR-335-3p, miR-133a-3p, miR-1298-5p, miR-144-5p, miR-1b-5p, mi5-582-3p, based on criteria of |Log2 Fold Change (FC)| > 1.5.

Figure 4.

In addition, based on these different expressed miRNAs, miRDB and miRTarBase databases were used to predicted potential target genes, which then performed bioinformatics analysis. In total, 184 GO terms were enriched in target genes according to the criteria FDR < 0.05. The results identified alteration in the DNA binding, transcription activity, cell differentiation, and nervous system development, with most predicted genes being enriched (Figure 4C-D).

miR-335-3p, miR-133a-3p were verified to have significant down-regulation Later, a 3.4-fold and 31-fold decrease in miR-335-3p, miR-133a-3p were further verified by qRT-PCR (Figure 4E). Since upregulation of miR133a has been previously reported to play an important role in postpartum osteoporosis(Li et al., 2018), and is function seems opposite to that in CUMS mice. Together, these results suggest that the decline of miR-335-3p in CUMS mice may be of more concern to us.

In addition, it is stress signals received by the brain that induced local miR-335-3p alterations in bone tissue. miR-335-3p is highly expressed in brain, according to mouse small noncoding RNA tissue atlas (https://www.ccb.uni-saarland.de/asra/) (Figure 4— figure supplement 1). So the expression of miR-335-3p were detected subsequently in five critical brain region, including medial prefrontal cortex (mPFC), NAC, amygdala, hypothalamus, hippocampus, which play key role in receiving and processing chronic stress. Only miR-335-3p was significantly decreased in NAC (Figure 4G). Consistently, the miR-335-3p level in the serum of CUMS mice was also reduced (Figure 4H). These results imply that brain-derived and femur-expressed miR-335-3p maintain bone homeostasis under normal conditions. Psychological stress leads to a decrease in miR-335-3p secreted by NAC, which in turn reduces levels in the femur.

2.5 miR-335-3p attenuates osteoclast activity

To further confirm whether miR-335-3p affected osteoclast activity, RAW264.7 cells were transfected with mimic-NC, mimic-miR-335-3p, inhibitor-NC or inhibitor-miR-335-3p, respectively (Figure 5—figure supplement 1). Cytochemical TRAP staining showed that mimic-miR-335-3p impaired the activity of osteoclasts, while inhibitor-miR-335-3p significantly promotes the activity of osteoclasts (Figure 5A-D). Next, the mRNA levels of Acp5, Nfatc1, and Mmp9 in the indicated RAW264.7 cells were determined by using qRT-PCR. And a significantly decreased expression in these osteoclast-related genes was observed in the mimic-miR-335-3p group, but an increased expression happened in the inhibitor-miR-335-3p group (Figure 5F). Additionally, in the inhibitor-miR-335-3p group, more intranuclear NFATC1, FOS and co-localization signals were observed by immunofluorescence staining, and this response was inhibited by transfection of mimic (Figure 5E).

Figure 5.

2.6 miR-335-3p alters osteoclast activity by targeting Fos

To seek the downstream molecular mechanisms that mediate the effects of miR-335-3p on the osteoclast, the ingenuity pathways analysis (IPA) was made based on miR-335-3p and its predicted target gene in miRDB and miRTarBase database. The results show that Fos, one of mRNA that miR-335-3p could target with a high prediction score and have a tight connection with osteoclast differentiation and maturation and osteoporosis (Figure 6A). FOS has been previously described as a regulator of transcription factors. It could induce NFATC1 translocation and finally increase the expression of osteogenic markers such as Acp5, Ctsk, Mmp9, and so on(Deng et al., 2022). Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis based on the above predicted target genes resulted in the same higher enrichment for transcription factors activity and DNA binding (Figure 6B).

Figure 6.

Figure 7.

Next, the dual-luciferase, co-transfection of miR-335-3p, suppressed the luciferase activity of the reporter containing wildtype Fos 3’UTR sequence by dual-luciferase reporter assay, which further strongly verifies the direct target relationship between miR-335-3p and Fos. But it failed to completely return to normal luciferase activity of two predicted target sites (Figure 6D) deleted construct, which means other binding sites still exist (Figure 6C). Furthermore, at the protein level, there was stronger signals that FOS co-localized with CTSK in the femur of CUMS mice, whose miR-335-3p level decreased (Figure 6E), although no significant difference at mRNA level (Figure 6—figure supplement 1). These results suggest that the enhanced osteoclastic activity under psychological stress may be due to the reduction of miR-335-3p in the femur which in turn slows down its inhibitory effect on Fos in osteoclast.

Moreover, the effect promoting osteoclastic activity of inhibitor-miR-335-3p could be reversed via T-5224, an inhibitor of FOS, measured via TRAP staining and counting (Figure 6F-H). The opposing effects on the mRNA mentioned above levels were also exerted via T-5224 (Figure 6I). These results further demonstrated that miR-335-3p regulates osteoclast activity by targeting Fos.

3 Discussion

Here, we found that psychological stress disturbs bone metabolism in vivo, and ultimately leads to bone loss in CUMS mice. We reported for the first time that miR-335-3p could target Fos to interfere with its translation and inhibit the activity of NFATC1 signaling pathway, besides, it was significantly reduced in mouse NAC, serum, and bone tissues under psychological stress. Therefore, the osteoclast activity is apparently promoted in CUMS mice, leading to the psychological stress induced osteoporosis. Our study demonstrated that miR-335-3p is a key mediator for nerve-regulated bone metabolism under psychological stress, and it might hold great potential to be a novel therapeutic target.

Previous studies have attributed the psychological stress-induced osteoporosis to the altered neurotransmitters, neuropeptides, and hormones that circulate and localize in tissues via neural and neuroendocrine signaling. In recent years, there has been an increasing focus on the regulatory role of miRNAs in the interaction between organs or systems, since miRNAs that protected by EV can be stably stored and trans, Ported to nearby or even isolated organs and system in vivo. Especially their ability to overcome the blood-brain barrier, broaden new horizons for the studies on the interaction between central nervous system and peripheral tissues, for example, bone. Currently, it has been reported that the expression profile of miRNAs was significantly altered in patients with neuropsychiatric disorders(Rao et al., 2013; Wang, 2021; Wu et al., 2022; Żurawek & Turecki, 2021). And some of these miRNAs contained in EVs play a vital role in central regulation in bone metaboliasm. miR-328a-3p and miR-150-5p enriched in EVs from hippocampus after traumatic brain injury accelerate bone healing(Xia et al., 2021). miR-384-5p in brain-derived EVs induces a bone-fat imbalance and bone loss in Alzheimer’s disease (Liu et al., 2023). Psychological stress, one of the most common type of neuropsychiatric disorders, is well established as a promoter of osteoporosis, so could miRNAs act as messengers between it and bone homeostasis?

Psychological stress correlated to several brain regions, such as the mPFC, the NAC, the amygdala, the hypothalamus, the hippocampus, etc. Expression levels of several miRNAs related to bone metabolism altered in brain under psychological stress. Hypotheses have been developed on miRNAs that may regulate bone loss under chronic stress(He et al., 2021). However, no study has yet reported whether altered miRNAs in the brain under stress can influence bone metabolism through blood circulation. In this study, for the identification of the key brain region that involved in chronic psychological stress-related bone metabolism, several miRNAs, which changed their expression levels in the CUMS femur comparing to the control ones, were examined in above mentioned brain regions, such as the mPFC, NAC, amygdala, hypothalamus. It was found that miR-335-3p, which is originally highly expressed in the brain and closely associated with bone metabolism, was significantly down-regulated in the NAC, as well as the serum of CUMS mice. This result is in accordance with the miRNAs sequencing results in the NAC of the CUMS mice reported by Ma K et al(Ma et al., 2019). Taken together, its indicated that the NAC may be one of the important brain regions what is responsible for the psychological stress induced osteoporosis via its secretion of miR-335-3p. In another word, these also suggest that the NAC derived and blood-trans, Ported miR-335-3p may have a protective effect on bone homeostasis under normal conditions, while in contrast, chronic stress diminishes the protective effect of miR-335-3p and leads to bone loss.

miR-335-3p is one of the maturation products after shear processing of its precursor mir-335. The other is miR-335-5p, also previously known as miR-335. They may be able to synergistically regulate the progression of diseases (Wilson et al., 2023). Previous psychiatric disease-related studies have mainly focused on miR-335. It play a vital role in the development of depression, both in brain and serum(Bocchio-Chiavetto et al., 2013; Fan et al., 2022; Li et al., 2015; Smalheiser et al., 2012).Therefore, miR-335-3p, which also be present in brain-derived EVs(Tian et al., 2022), could be another significant miRNA involved in psychiatric disorders. At the same time, miR-335-3p also plays a crucial role in the skeletal system. One study reported a significant downregulation of miR-335-3p in the mouse bone tissue of hindlimb unloading osteoporosis model(Huai et al., 2020). Our study further revealed the mechanism of miR-335-3p: It targets Fos so as to influence the activity of FOS-NFATC1 signaling pathway. As a result, the down-regulation of miR-335-3p under chronic psychological stress promotes FOS expression and osteoclast activity, leading to osteoporosis.

However, although it’s well verified the consistency of the trend of miR-335-3p in the NAC, serum, and femur, respectively, suggesting miR-335-3p may derived from NAC and targets bone via blood circulation system. Further direct validation is still needed, such as using miRNA tracer mice to demonstrate the importance of brain-derived miRNA. Meanwhile, the potential impact of exogenous miR-335-3p on other bone cells, such as osteoblasts and mesenchymal stem cells, remains to be investigated. Moreover, in addition to the well demonstrated protective effect of miR-335-3p on bone mass in vitro, in vivo evidences need to be established considering the further therapeutic strategies for psychological stress-induced osteoporosis. Nevertheless, our study confirmed that miR-335-3p may works as a protective effector for bone metabolism during normal conditions and is a key regulator during the pathogenesis of psychological stress-induced osteoporosis.

In conclusion, our study suggests that enhanced osteoclast activity is a major cause of psychological stress-induced osteoporosis. It is found that miR-335-3p regulates osteoclast activity and protects bone metabolic balance by suppressing FOS expression and FOS-NAFCT1 signaling pathway. However, this protective effect is diminished by chronic psychological stress. These findings provide new insights into the mechanisms of osteoporosis in patients with chronic psychological stress and may aid in the development of new therapeutic strategies.

4 Materials and Methods

4.1 Mice

Male C57BL/6J mice, aged 6 weeks, purchased from Shanghai Model Organisms Center, Inc., were raised at Hubei Laboratory Animal Center, Tongji University, with a 12-hour light/12-hour dark cycle, a room temperature of 22-24 degrees, humidity of 45-55% and free access to food and water. All experimental procedures were approved by and performed in accordance with the standards of the Animal Welfare Committees of Tongji University in Shanghai, China (Number: TJ-HB-LAC-2023-39).

4.2 Establishment of CUMS mouse models

After 1 week of acclimation, 18 mice were divided into the control group and CUMS model Group. CUMS was carried out as described previously(Antoniuk et al., 2019; Willner, 2017) with slight modifications, and the stimulations included 24-hour food deprivation, 24-hour water deprivation, 1 min tail pinch, 7-hour 45° cage tilt, 24-hour wet cage, 5 minutes forced swimming, 5-hour physical restraint, 24-hour day/night reversal. These stressors were randomized to CUMS mice and continued for 6 weeks, while control mice were handled daily in the housing room. Body weight was monitored once every week at 15:30-16:30. Figure 1 A shows the schedule of this study.

4.3 Behavior Analysis

All behavioral tests were done in a dimly lit and quiet room. Animals were acclimatized in the room for at least 2 hours before formal testing. The testing period did not span the diurnal period of the mice. The behavioral performance was monitored and recorded using a digital camera connected to a computer running the Tracking Master software.

4.3.1 Open Field Test (OFT)

The mice were placed in a white open field (40 cm long, 40 cm wide, and 40cm high). The total horizontal distance traveled in 10 minutes, as well as the time and distance traveled in the center area (20 cm* 20 cm) were recorded and calculated. 75% alcohol was used to clean the box during the tests of different mice.

4.3.2 Tail suspension Test (TST)

The mice were suspended individually, with one end of an adhesive tape fixed about 1cm from the tip of the tail, and the other end stocked to the TST apparatus. The mice were held at a 25-30 cm distance between the tip of their nose and the ground. Their behavior was monitored for 6 min, and the immobility time in the latter 4 min was calculated.

4.3.3 Sucrose preference test (SPT)

Mice were single housed and provided with two identical bottles of drinking water. First, Mice were acclimatized to two bottles of normal drinking water for 24 hours. Then, they were given one bottle of normal water and one bottle of 1% sucrose water for the next 24 hours and switched the positions of the two water bottles at the 12th hour. Mice were water and food deprived for 24 hours before the test. During the test, they had free access to normal water and 1% sucrose water, and the positions of the bottles were switched in the middle of the test. The liquid consumption in each bottle over 12 hours was measured, and the sucrose preference index was calculated by dividing the consumption of sucrose solution by the total amount of fluid (water + sucrose).

4.4 Serum Elisa

The concentrations of corticosterone and norepinephrine in serum were measured using elisa kits (Jiangsu MEIMIAN Industrial Co. Ltd., Yancheng, China).

4.5 Micro–CT Analysis

The left femur was dissected from each group of mice, fixed in 4% PFA for 48 hours, and stored in 0.4% PFA until scanning. The femurs were analyzed by micro–CT 50 (Scanco Medical, Zurich, Switzerland) with a scan resolution of a 10-μm slice, a voltage of 70 kV, and a current of 200 μA. To quantify trabecular bone parameters, 60 CT slices immediately proximal to the distal femoral growth plate were analyzed. Bone volume (BV), bone volume over tissue volume (BV/TV), bone mineral density (BMD), trabecular separation (Tb.Sp), trabecular number (Tb.N), trabecular thickness (Tb.Th) were quantified according to the standard procedures.(Y. Kim et al., 2021)

4.6 Histologic Analysis

The right femur samples were fixed with 4% PFA for 48 hours, decalcified in 0.5 mol/L ethylene diamine tetraacetic acid (EDTA, pH 7.4) for 3 weeks, and paraffin-embedded. Next, the samples were sectioned into 4um-thick slices. Hematoxylin & eosin (H&E) (Beyotime, China; C0105S) was conducted to evaluate histological morphology. Osteogenic activities were detected by masson (Servicebio, China; G1006) and immunohistochemical staining for OSX (Abcam, Cambridge, MA, USA; ab22552) and OCN (Affinity, China; DF12303). For double fluorochrome labeling, a total of 2 intraperitoneal injections of calcein (5 mg/kg; Sigma-Aldrich, St. Louis, MO, USA) were given at 7-day intervals at the end of CUMS modeling. Osteoclastic activities were detected by immunohistochemical staining for MMP9 (Abclonal, China; A0289) and TRAP staining (Sigma-Aldrich, St. Louis, MO, USA; 387A). The numbers of osteoclasts (N. OCs/BS/mm) and the TRAP positive surface of osteoclast (OC.S/B) were calculated. At the level of mechanism validation, immunohistochemical staining for FOS (Cell Signaling, USA; #2250), CTSK (Santa Cruz Biotechnology, CA, USA; sc-48353), NFATC1 (Santa Cruz Biotechnology, CA, USA; sc-7294), HOCHEST (ThermoFisher, Waltham, USA; H3570) was used to observe the location and abundance of protein expression

4.7 Quantitative RT-PCR Analysis

Total RNA was extracted from tissues or cells using RNAiso Plus (Takara Biotechnology, Japan; 9109) according to the manufacturer’s guidelines. RNA concentration was assessed by NanoDrop One. For mRNAs, cDNA was synthesized using PrimeScript™ RT Master Mix (Takara Biotechnology, Japan; RR036A),and a reaction mix was prepared using SYBR green master mix (Yeasen, Shanghai, China; 11201ES08). For miRNAs, miDETECT A TrackTM miRNA qRT-PCR Starter Kit (RIBOBIO, Guangzhou, China; C10712-1) was used. The expression levels of the target gene were calculated by the comparative Ct (2^−ΔΔCT) method using GAPDH or U6 for normalization, respectively, using CFX Opus 96 (Bio-Rad Laboratories, Hercules, CA, USA). The amplifications of the miRNAs were carried out by using miR-335-3p-F, miR-133a-3p-F, miR-1298-5p-F, miR-144-5p-F, miR-1b-5p-F, miR-582-3p-F, miR-141-3p-F, and U6-F as forward primer and unified reverse primer (URP) as a reverse primer to amplify miR-335-3p, miR-133a-3p, miR-193b-3p, miR-1298-5p, miR-144-5p, miR-1b-5p, miR-582-3p, miR-141-3p and U6, respectively. The primer sequences were as follows:

4.8 microRNA sequencing (sRNA-Seq)

Total RNA was extracted from M-exos and HM-exos using the Trizol reagent (Takara Biotechnology, Japan; 9109), and RNA integrity was assessed using a 2100 bioanalyzer (Agilent Technologies, CA, USA). An Illumina TruSeq Small RNA kit (Illumina, San Diego, CA, USA) was used to construct the library, and a high-throughput sequencing platform was used to sequence the enriched 18–32 nt Small RNA fragments. The differentially expressed miR NAs in the two groups of exosomes were screened (p-value < 0.05), and miRDB (http://www.mirdb.org/) and DIANA (https://dianalab.e-ce.uth.gr/) databases were used to predict the target genes of the miRNAs, and KEGG enrichment analysis was performed for the identified target genes.

4.9 3’UTR luciferase reporter assay

A fragment of the Fos 3′UTR containing wildtype (WT) or mutant (Mut) predicted two binding sites from TargetScan (http://www.targetscan.org) database for miR-335-3p was inserted into the pMIR-REPORT vector. HEK293 cells were cotransfected with either Fos-WT or Fos-Mut vector and either miR-335-3p mimic or miRNA mimic negative control (mimic NC) using Lipofectamine 2000 (Invitrogen, USA; 11668019). Luciferase activity was determined at 48 h post-transfection using the Dual-Glo Luciferase Assay System (Promega, Madison, WI; E2940) following the manufacturer’s instructions. Data were normalized by dividing firefly luciferase activity with that of Renilla luciferase enzyme activity.

4.10 Cell culture and transfection

The Raw 264.7 cells (gifts from Dr. Shengbing Yang, Shanghai Key Laboratory of Orthopaedic Implants, Department of Orthopaedic Surgery, Shanghai Ninth People’s Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, 200125, China) were used for osteoclastic differentiation induction via RANKL (R&D Systems, USA; 462-TEC) (50nM). The entire inducing process lasted 4 days, and medium refreshed every three days. In addition, BMSCs were transfected with miR-335-3p mimic, miR-335-3p inhibitor and NC that synthesized by RiboBio (Guangzhou, China) via Fugene HD (Promega, Madison, WI, USA; E2311) when cells confluency reached 50%.

4.11 Statistical analysis

Data are represented as the mean±SD. Differences between 2 groups were determined by Student’s t test. Differences among 3 groups were determined by one-way ANOVA with Tukey’s post hoc test. All experiments were repeated at least three times. A value of P <0.05 was considered significant.

5 Data availability

All data generated or analysed during this study are included in the manuscript and supporting file.

7 Article and author information

Jiayao Zhang

Department of Prosthodontics, Stomatological Hospital and Dental School of Tongji University, Shanghai Engineering Research Center of Tooth Restoration and Regeneration, Shanghai, China.

Contribution: Data curation, Investigation, Writing - original draft, Writing - review and editing

For correspondence: michellezhang_07@163.com

Competing interests: No competing interests declared

Juan Li

Department of Prosthodontics, Stomatological Hospital and Dental School of Tongji University, Shanghai Engineering Research Center of Tooth Restoration and Regeneration, Shanghai, China.

Contribution: Data curation, Investigation, Validation; Formal analysis

For correspondence: 15025002369@163.com

Competing interests: No competing interests declared

Jiehong Huang

Department of Neurology and Neurological Rehabilitation, Shanghai Yangzhi Rehabilitation Hospital, School of Medicine, Tongji University, Shanghai, China Contribution: Methodology; Software; Validation; Formal analysis

For correspondence: jh_huang1996@163.com

Competing interests: No competing interests declared

Xuerui Xiang

Department of Neurology and Neurological Rehabilitation, Shanghai Yangzhi Rehabilitation Hospital, School of Medicine, Tongji University, Shanghai, China Contribution: Methodology; Software; Validation; Formal analysis

For correspondence: sherry199503@163.com

Competing interests: No competing interests declared

Ruoyu Li

Department of Neurology and Neurological Rehabilitation, Shanghai Yangzhi Rehabilitation Hospital, School of Medicine, Tongji University, Shanghai, China Contribution: Methodology; Software; Validation; Formal analysis

For correspondence: liry_1996@163.com

Competing interests: No competing interests declared

Yun Zhai

Department of Prosthodontics, Stomatological Hospital and Dental School of Tongji University, Shanghai Engineering Research Center of Tooth Restoration and Regeneration, Shanghai, China.

Contribution: Methodology; Data curation

For correspondence: 904890219@qq.com

Competing interests: No competing interests declared

Shuxian Lin

Department of Prosthodontics, Stomatological Hospital and Dental School of Tongji University, Shanghai Engineering Research Center of Tooth Restoration and Regeneration, Shanghai, China.

Contribution: Conceptualization; Funding acquisition; Supervision; Writing - review and editing

For correspondence: shuxian.lin@hotmail.com

Contributed equally with: Weicai Liu

Competing interests: No competing interests declared

Weicai Liu

Department of Prosthodontics, Stomatological Hospital and Dental School of Tongji University, Shanghai Engineering Research Center of Tooth Restoration and Regeneration, Shanghai, China.

Contribution: Conceptualization; Funding acquisition; Supervision; Writing - review and editing

For correspondence: vogi@163.com

Contributed equally with: Shuxian Lin

Competing interests: No competing interests declared

8 Funding

National Natural Science Foundation of China (81873709)

• Shuxian Lin

Shanghai Health Committee (20204Y0278)

• Shuxian Lin

National Natural Science Foundation of China (8227032547)

• Weicai Liu

Natural Science Foundation of Shanghai (20ZR1463000)

• Weicai Liu