Nodakenin Ameliorates Ovariectomy-Induced Bone Loss by Regulating Gut Microbiota

Disordered gut microbiota (GM) structure and function may contribute to osteoporosis (OP). Nodakenin has been shown to ameliorate osteoporosis; however, its anti-osteoporotic mechanism is unknown. This study aimed to further reveal the mechanism of the anti-osteoporotic action of nodakenin from the perspective of the microbiome and metabolome. An osteoporosis model was induced in mice through ovariectomy (OVX), with bone mass and microstructure assessed using μCT. Subsequently, ELISA and histologic examination were used to detect biochemical indicators of bone conversion and intestinal morphology. Using metabolomics and 16S rRNA sequencing, it was possible to determine the composition and abundance of the gut microbiota in feces. The results revealed that nodakenin treatment improved the bone microstructure and serum levels of bone turnover markers, and increased the intestinal mucosal integrity. 16S rRNA sequencing analysis revealed that nodakenin treatment decreased the relative abundance of Firmicutes and Patescibacteria, as well as the F/B ratio, and elevated the relative abundance of Bacteroidetes in OVX mice. In addition, nodakenin enhanced the relative abundance of Muribaculaceae and Allobaculum, among others, at the genus level. Moreover, metabolomics analysis revealed that nodakenin treatment significantly altered the changes in 113 metabolites, including calcitriol. A correlation analysis revealed substantial associations between various gut microbiota taxa and both the osteoporosis phenotype and metabolites. In summary, nodakenin treatment alleviated OVX-induced osteoporosis by modulating the gut microbiota and intestinal barrier.


Introduction
Osteoporosis (OP) is a condition characterized by systemic disruptions in bone metabolism, leading to decreased bone mineral density (BMD) and structural deterioration [1].The main population that suffers from osteoporosis is postmenopausal women, where metabolic disturbances in bone homeostasis are triggered by aging or estrogen deficiency [2].Due to the growing prevalence of aging worldwide, osteoporosis has emerged as a global public health issue.Currently, the prevention and treatment of OP mainly depend on lifestyle modifications and various types of pharmacologic interventions [3,4].Nevertheless, osteoporosis treatment is sometimes hindered by the financial burden and adverse reactions associated with drugs [5], so there is an urgent need to develop more efficient and safer treatment options.
Gut microbiota can interact with various organs and systems of the body and are essential influencers of metabolism-related diseases [6].Almost 10 trillion bacteria in the human gastrointestinal system can interact with the host and play a crucial role in influencing metabolic illnesses like obesity, diabetes, and osteoporosis [7][8][9].Clinical studies have shown that gut microbiota can influence bone quality and strength [10][11][12].In addition, recent studies have shown that estrogen deficiency can induce imbalances in the gut microbiota (GM), impair intestinal mucosal barrier function, and increase reactivity in the inflammatory system [13][14][15].However, their complex mechanisms have not been completely clarified.Hence, the mechanism underlying the gut/bone axis deserves further investigation, and it would be interesting to search for an effective drug candidate that exerts anti-osteoporotic effects by modulating the gut microbiota.
Nodakenin, a furanocoumarin glycoside isolated from the roots of Angelica sinensis, has anti-inflammatory, antimicrobial, antioxidant, and antitumor effects [16][17][18][19][20].It has recently been found that nodakenin can increase the subchondral bone volume and thereby reduce osteoarthritis [21].We also demonstrated in our previous study that nodakenin alleviated osteoporosis by regulating bone formation and bone resorption [22].Nodakenin exhibits significantly low oral bioavailability yet is efficacious.This indicates that nodakenin may exercise its pharmacological effects by modulating the gut microbiota and the resulting metabolites [23,24].However, whether the mechanism by which nodakenin inhibits bone loss is related to changes in the gut microbiota has not been reported to date.
In this study, we used an OVX mouse model to mimic postmenopausal estrogen deficiency-induced osteoporosis.Furthermore, we utilized fecal microbiome and metabolome analysis to investigate how nodakenin alters the gut microbiota and its metabolite profile.We found that nodakenin treatment may ameliorate OVX-induced osteoporosis in mice by modulating the gut microbiota and repairing the gut barrier.

Nodakenin Alleviates Osteoporosis in OVX Mice
As demonstrated, we administered nodakenin medication for a duration of 6 weeks following an 8-week period of ovariectomy in mice.This was carried out to assess the impact of nodakenin on the bone loss generated by ovariectomy in mice with OP (Figure 1A).OVX mice gained significant body weight compared with the Sham group, while nodakenin treatment reduced the OVX-induced weight increase.There were no significant differences in the weights of the liver, kidney, and spleen (Supplementary Figure S1).µCT confirmed a significant decrease in BMD, BV/TV, and Tb.N and a significant increase in SMI in OVX mice compared to the Sham group, suggesting the successful modeling of osteoporosis in mice.BMD and bone microstructural parameters were significantly improved in the mice after nodakenin treatment (Figure 1B,C).It is suggested that nodakenin can alleviate OVX-induced bone loss.

Nodakenin Regulates Estrogen Deficiency-Induced Bone Turnover
In order to evaluate the impact of nodakenin on osteoporosis, we examined the changes in serum bone turnover marker levels via ELISA.The serum levels of the B-ALP and BGP bone formation markers were significantly elevated in nodakenin-treated mice compared to OVX mice, while TRAP and RANKL bone resorption markers were decreased (Figure 2A-D).In addition, we used TRAP staining to assess the femoral histomorphometry in order to investigate the impact of nodakenin on estrogen shortage-induced bone loss.The findings demonstrated that a decrease in estrogen levels led to a notable rise in osteoclasts, while nodakenin exhibited the ability to impede osteoclastogenesis in mice (Figure 2E,F).

Nodakenin Regulates Estrogen Deficiency-Induced Bone Turnover
In order to evaluate the impact of nodakenin on osteoporosis, we examined the changes in serum bone turnover marker levels via ELISA.The serum levels of the B-ALP and BGP bone formation markers were significantly elevated in nodakenin-treated mice compared to OVX mice, while TRAP and RANKL bone resorption markers were decreased (Figure 2A-D).In addition, we used TRAP staining to assess the femoral histomorphometry in order to investigate the impact of nodakenin on estrogen shortage-induced bone loss.The findings demonstrated that a decrease in estrogen levels led to a notable rise in osteoclasts, while nodakenin exhibited the ability to impede osteoclastogenesis in mice (Figure 2E,F).

Nodakenin Alters OVX-Induced Gut Microbiota Imbalance in Mice
Recent research has demonstrated that the gut microbiota significantly influences bone metabolism [25,26].Therefore, we analyzed the composition of the gut microbiota via the systematic sequencing of 16S rRNA gene amplicons in fecal samples.First, the Sobs and Shannon indices flattened the dilution curves, suggesting that the sequencing quantities could adequately capture the majority of the microbial diversity present in the samples (Supplementary

Nodakenin Alters OVX-Induced Gut Microbiota Imbalance in Mice
Recent research has demonstrated that the gut microbiota significantly influences bone metabolism [25,26].Therefore, we analyzed the composition of the gut microbiota via the systematic sequencing of 16S rRNA gene amplicons in fecal samples.First, the Sobs and Shannon indices flattened the dilution curves, suggesting that the sequencing quantities could adequately capture the majority of the microbial diversity present in the samples (Supplementary Figure S2).Chao 1 analysis showed no significant differences in the microbial community richness between groups (Figure 3A).The Shannon index showed an increased diversity of gut microbial species after nodakenin treatment compared to OVX mice (Figure 3B).Both UniFrac-based principal coordinates analysis (PCoA) and hierarchical clustering showed significant differences in the gut microbiota of OVX mice compared to nodakenin-treated mice (Figure 3C,D).The findings indicate that nodakenin may ameliorate OVX-induced osteoporosis by regulating the imbalance of the gut microbiota.To analyze particular alterations in the gut microbiota, we assessed the relative abundance of the top 10 or 20 major taxa.At the phylum level, Bacteroides and Firmicutes were the most dominant phyla, accounting for more than 90% of the total microbial composition (Figure 4A).The research showed that there was a higher proportion of Firmicutes and Patescibacteria and a lower proportion of Bacteroidetes in OVX mice compared to the Sham mice.Nodakenin treatment decreased the relative abundance of Firmicutes and Patescibacteria and increased the relative abundance of Bacteroidetes (Figure 4B).In addition, the ratio of Firmicutes and Bacteroidetes (F/B) is often used as an assessment of the impact of the gut microbe on numerous disease species.Our study found a notable rise in the F/B ratio in OVX mice compared to the Sham mice, while nodakenin treatment decreased the F/B ratio (Figure 4C).At the genus level, the relative abundances of Muribaculaceae, Parasutterella, and Coriobacteriaceae_UCG-002 were significantly reduced, and the relative abundances of Muribaculum and Roseburia were increased in OVX mice.Meanwhile, nodakenin treatment restored the relative abundances of Muribaculaceae, Parasutterella, Faecalibaculum, Allobaculum, and Coriobacteriaceae_UCG-002, and decreased the relative abundances of Roseburia, Blautia, and Lachnoclostridium (Figure 4D,E).Additionally, to To analyze particular alterations in the gut microbiota, we assessed the relative abundance of the top 10 or 20 major taxa.At the phylum level, Bacteroides and Firmicutes were the most dominant phyla, accounting for more than 90% of the total microbial composition (Figure 4A).The research showed that there was a higher proportion of Firmicutes and Patescibacteria and a lower proportion of Bacteroidetes in OVX mice compared to the Sham mice.Nodakenin treatment decreased the relative abundance of Firmicutes and Patescibacteria and increased the relative abundance of Bacteroidetes (Figure 4B).In addition, the ratio of Firmicutes and Bacteroidetes (F/B) is often used as an assessment of the impact of the gut microbe on numerous disease species.Our study found a notable rise in the F/B ratio in OVX mice compared to the Sham mice, while nodakenin treatment decreased the F/B ratio (Figure 4C).At the genus level, the relative abundances of Muribaculaceae, Parasutterella, and Coriobacteriaceae_UCG-002 were significantly reduced, and the relative abundances of Muribaculum and Roseburia were increased in OVX mice.Meanwhile, nodakenin treatment restored the relative abundances of Muribaculaceae, Parasutterella, Faecalibaculum, Allobaculum, and Coriobacteriaceae_UCG-002, and decreased the relative abundances of Roseburia, Blautia, and Lachnoclostridium (Figure 4D,E).Additionally, to find bacterial taxonomic markers linked to OVX-induced OP, we conducted a further LEfSe analysis (Figure 4F) and created branching diagrams (Figure 4G).Significant differences in bacterial taxonomic composition were observed in different groups.In the nodakenin group, g__Allobaculum, s__gut_metagenome, f__Atopobiaceae, g__Coriobacteriaceae_UCG_002, s__Kribbella_sp, and g__Kribbella were identified (Figure 4G).Taken together, our data suggest that nodakenin treatment modulates OVX-induced GM imbalance in OP mice.Spearman correlation analysis was used to evaluate the associations between alterations in the gut microbiota composition and the markers associated with osteoporosis.At the phylum level, Patescibacteria and Firmicutes exhibited a positive correlation with RANKL, TRAP, and SMI, while displaying a negative correlation with Tb.N, BMD, BV/TV, BGP, and B-ALP.Bacteroidota showed a positive correlation with Tb.N, BMD, BV/TV, BGP, and B-ALP and a negative correlation with RANKL, TRAP, and SMI.Fusobacteriota showed a positive association with Tb.N, BMD, BV/TV, and BGP, while it exhibited a negative association with RANKL and TRAP.Proteobacteria showed a positive correlation with BMD and BV/TV (Figure 5A).At the genus level, Muribaculaceae, Faecalibaculum, and Coriobacteriaceae_UCG-002 were positively correlated with Tb.N, BMD, BV/TV, and bone formation-associated BTMs, and negatively correlated with TRAP, SMI, and RANKL.Parasutterella was positively correlated with Tb.N, BMD, BV/TV, and B-ALP.Muribaculum and Bacteroides were positively correlated with RANKL, TRAP, and SMI.Muribaculum was negatively correlated with BV/TV and B-ALP.Bacteroides was negatively correlated with BV/TV, Tb.N, and BGP.Roseburia was negatively correlated with B-ALP and positively correlated with TRAP (Figure 5B).Muribaculum was negatively correlated with BV/TV and B-ALP.Bacteroides had a negative correlation with BV/TV, Tb.N, and BGP.Roseburia had a negative correlation with B-ALP and a substantial positive correlation with TRAP (Figure 5B).

Metabolite Profiles Change Caused by Nodakenin Treatment
Alterations in the composition of gut microbes can result in changes in intestinal metabolites.To further explore the effect of nodakenin on intestinal metabolites, an untargeted metabolomic analysis of feces was performed using LC-MS.The PLS-DA analysis showed that the sample points in the OVX and OVX-NK groups did not overlap at all.This means that the mixtures of metabolites were very different (Figure 6A).The variable significances of VIP > 1 and p < 0.05 were used to screen for differential metabolites between the two groups.Out of a total of 2080 metabolites, 113 showed significant differences.Among them, 48 exhibited considerably higher relative levels, while 65 had significantly lower relative levels (Figure 6B).The top 20 differential metabolites with VIP values were selected for presentation, including LPG 22:6, traumatic acid, and calcitriol (Figure 6C).Notably, the group that was administered nodakenin had much higher levels of calcitriol, a metabolite that helps bones form (Figure 6C).Further metabolite enrichment analysis showed that nodakenin treatment had significant effects on multiple metabolic pathways, involving endocrine and other factor-regulated calcium reabsorption, tuberculosis, retrograde endocannabinoid signaling, alpha-Linolenic acid metabolism, thermogenesis, riboflavin metabolism, ether lipid metabolism, mineral absorption, etc., (Figure 6D).

Metabolite Profiles Change Caused by Nodakenin Treatment
Alterations in the composition of gut microbes can result in changes in intestinal metabolites.To further explore the effect of nodakenin on intestinal metabolites, an untargeted metabolomic analysis of feces was performed using LC-MS.The PLS-DA analysis showed that the sample points in the OVX and OVX-NK groups did not overlap at all.This means that the mixtures of metabolites were very different (Figure 6A).The variable significances of VIP > 1 and p < 0.05 were used to screen for differential metabolites between the two groups.Out of a total of 2080 metabolites, 113 showed significant differences.Among them, 48 exhibited considerably higher relative levels, while 65 had significantly lower relative levels (Figure 6B).The top 20 differential metabolites with VIP values were selected for presentation, including LPG 22:6, traumatic acid, and calcitriol (Figure 6C).Notably, the group that was administered nodakenin had much higher levels of calcitriol, a metabolite that helps bones form (Figure 6C).Further metabolite enrichment analysis showed that nodakenin treatment had significant effects on multiple metabolic pathways, involving endocrine and other factor-regulated calcium reabsorption, tuberculosis, retrograde endocannabinoid signaling, alpha-Linolenic acid metabolism, thermogenesis, riboflavin metabolism, ether lipid metabolism, mineral absorption, etc., (Figure 6D).Consistent with our metabolomics findings, we observed a notable increase in serum levels of calcitriol, a pivotal metabolite, in the nodakenin-treated group compared to OVX mice (Figure 7A).Calcitriol functions by binding to the vitamin D receptor (VDR).Remarkably, we found a significant elevation in VDR expression within the bone tissue of nodakenin-treated mice (Figure 7B).Furthermore, the mRNA levels of key osteogenic marker genes, including alkaline phosphatase (ALP), collagen type I (Col-1), and osteocalcin (OCN), exhibited significant increases (Figure 7C-E).
In addition, Spearman correlation analysis was performed on intestinal metabolites (VIP top 20), significantly different gut microbiota, and osteoporosis indicators (Figure 8).We discovered substantial associations among the aforementioned metrics, indicating that the beneficial impact of nodakenin on bone loss in mice produced by OVX involves the regulation and recovery of gut microbiota and metabolite profiles.Nevertheless, additional verification is required to establish the potential pathways that influence the gut microbiota and metabolites in OP therapy.Consistent with our metabolomics findings, we observed a notable increase in serum levels of calcitriol, a pivotal metabolite, in the nodakenin-treated group compared to OVX mice (Figure 7A).Calcitriol functions by binding to the vitamin D receptor (VDR).Remarkably, we found a significant elevation in VDR expression within the bone tissue of nodakenin-treated mice (Figure 7B).Furthermore, the mRNA levels of key osteogenic marker genes, including alkaline phosphatase (ALP), collagen type I (Col-1), and osteocalcin (OCN), exhibited significant increases (Figure 7C-E).
In addition, Spearman correlation analysis was performed on intestinal metabolites (VIP top 20), significantly different gut microbiota, and osteoporosis indicators (Figure 8).We discovered substantial associations among the aforementioned metrics, indicating that the beneficial impact of nodakenin on bone loss in mice produced by OVX involves the regulation and recovery of gut microbiota and metabolite profiles.Nevertheless, additional verification is required to establish the potential pathways that influence the gut microbiota and metabolites in OP therapy.

Nodakenin Ameliorates Intestinal Barrier Impairment in OVX Mice
To explore the effect of nodakenin on intestinal barrier integrity in OVX-induced OP mice, we analyzed the intestinal barrier via an intestinal villus morphology assessment (Figure 9A).Compared with that in the Sham group, the height of the intestinal villi was decreased in OVX mice (Figure 9B), while the depth of intestinal crypts was increased (Figure 9C) and the V/C ratio (Figure 9D) was significantly decreased.It was also found that the OVX mice had much lower levels of intestinal Occludin (Figure 9E) and Zonulae-1 (ZO-1) (Figure 9F) mRNA, while they had much higher levels of IL-1β and TNFα mRNA (Figure 9G,H).Nodakenin treatment elevated the intestinal villus height, and Occludin and ZO-1 mRNA levels, and decreased the depth of crypts and the expression of inflammatory factors in the intestines of OVX mice (Figure 9B-H).It is suggested that nodakenin can improve OVX-induced intestinal barrier impairment.

Nodakenin Ameliorates Intestinal Barrier Impairment in OVX Mice
To explore the effect of nodakenin on intestinal barrier integrity in OVX-induced OP mice, we analyzed the intestinal barrier via an intestinal villus morphology assessment (Figure 9A).Compared with that in the Sham group, the height of the intestinal villi was decreased in OVX mice (Figure 9B), while the depth of intestinal crypts was increased (Figure 9C) and the V/C ratio (Figure 9D) was significantly decreased.It was also found that the OVX mice had much lower levels of intestinal Occludin (Figure 9E) and Zonulae-1 (ZO-1) (Figure 9F) mRNA, while they had much higher levels of IL-1β and TNFα mRNA (Figure 9G,H).Nodakenin treatment elevated the intestinal villus height, and Occludin and ZO-1 mRNA levels, and decreased the depth of crypts and the expression of inflammatory factors in the intestines of OVX mice (Figure 9B-H).It is suggested that nodakenin can improve OVX-induced intestinal barrier impairment.

Discussion
The burden of OP and OP-induced fractures is increasing with global population aging [27].In the present study, our data demonstrated that nodakenin treatment elevated bone mass, modulated serum bone turnover indexes, and decreased osteoclast activity.Generally, the glycosides contained in herbal medicines are poorly absorbed in the intestine.
Due to their low bioavailability, glycosides are retained in the intestine for a considerable period.After being metabolized by the gut microbiota, they are mainly absorbed in the form of glycosides, although it has been shown that the oral bioavailability of nodakenin is extremely low [23], which greatly limits its absorption.However, the oral administration of nodakenin still significantly elevated bone mineral density in OVX mice.It is suggested that it may exert its pharmacological effects by affecting the gut microbiota and the produced metabolites.There is increasing evidence that GM imbalance is associated with the pathogenesis and clinical manifestations of OP [28].Therefore, the present study focused on the possible mechanism of action of nodakenin in suppressing OP in terms of the gut/bone axis.
The gastrointestinal tract is acknowledged as a distinct organ system that has a significant impact on the maintenance of bone health.Recent research has indicated that there are connections between the stomach, bone signaling pathways, and microbiota that play a role in the regulation of bone health [29][30][31][32].Firmicutes can activate osteoclasts and exacerbate inflammation [33].The abundance and proportion (F/B ratio) of Firmicutes and Bacteroides are important factors indicating gut microbiota disorders, which may reflect the process of disease development.In this study, Firmicutes and Bacteroides were the two dominant species, accounting for more than 90% of the gut microbiota at the phylum level.In the OVX group, the relative abundance of Firmicutes was significantly increased and significantly negatively correlated with BMD, BV/TV, etc., and the F/B ratio was significantly increased.The increase in the abundance of Firmicutes and the increase in the F/B ratio were reversed after nodakenin treatment.Consistent with the present study, it has been found that the ratio of Firmicutes to Bacteroides was negatively correlated with bone volume [34,35].Furthermore, Bacteroides was significantly reduced in OVX mice and nodakenin treatment restored its relative abundance.A positive correlation between Bacteroides and calcium uptake has been reported [36].At the genus level, nodakenin treatment restored OVX-induced reductions in the relative abundances of Muribaculaceae, Parasutterella, Faecalibaculum, Allobaculum, and Coriobacteriaceae_UCG-002 compared to the OVX group.Muribaculaceae and Parasutterella were positively correlated with BMD, BV/TV, Tb.N, and B-ALP.Palmatine and intermittent parathyroid hormone (PTH) have been reported to increase the abundance of Muribaculaceae and increase bone mass in OVX rats [37,38].Parasutterella is involved in Kefir Peptides' prevention of estrogen deficiencyinduced bone loss and GM disorders [13].Allobaculum, a known butyrate producer is involved in the anti-osteoporotic effects of several active compounds [39,40].Faecalibaculum and Coriobacteriaceae_UCG-002, although not yet reported to be associated with increased bone volume, play an important role in the maintenance of intestinal homeostasis and in preventing intestinal damage [41][42][43].In conclusion, these flora may be involved in the correction of nodakenin anti-OVX-induced GM ecological imbalances.
Metabolomics is a technique used to investigate and analyze small-molecule metabolites in order to understand and describe pathological alterations.The use of metabolomics allows for the screening of metabolic biomarkers, which can help to elucidate the pathogenic mechanisms of osteoporosis and develop new drug targets [44].Our study showed that nodakenin treatment significantly altered the fecal metabolite profile in OVX mice.During nodakenin alleviation of OP, metabolomic screening yielded a significant upregulation of calcitriol, a differential metabolite closely related to osteoporosis, and four signaling pathways associated with calcitriol were identified by KEGG enrichment analysis.Therefore, calcitriol might be a key metabolite regulated by nodakenin.Calcitriol has been reported as a drug used in the clinical treatment of osteoporosis, and it is also one of the most important active metabolites of vitamin D3, which enhances intestinal calcium absorption [45,46].It is also a potent transcriptional activator of osteogenesis-related genes encoding type I collagen (Col-I), alkaline phosphatase (ALP), and osteocalcin (OCN) in osteoblasts, which regulate calcium homeostasis through stimulating osteoblast activity, thereby combating osteoporosis [47,48].Similarly, a significant increase in calcitriol was detected in the fecal metabolites of nodakenin-treated mice.Also, a significant increase in the levels of calcitriol was confirmed in the serum of nodakenin-treated mice, which was accompanied by an increase in the expression of the osteogenic marker genes ALP, OCN, and Col-1 in bone tissue.This suggests that calcitriol is a key metabolite of nodakenin in alleviating OP.In addition, both epidemiological investigations and basic studies have shown that there is a link between obesity and osteoporosis, with multiple important and intersecting molecular and cellular signaling pathways [49].Estrogen deficiency-induced weight gain has also been induced in mouse models of osteoporosis [50].LPG 22:6 has been found to be significantly elevated in high-fat-induced obese mice [51].Our study showed that LPG 22:6 was elevated in OVX mice, whereas nodakenin treatment decreased the level of LPG 22:6.Traumatic acid is a potent wound-healing agent in plants with antioxidant properties.It can also control lipid accumulation in adipocytes and reduce obesity caused by a high-fat diet [52].Our study showed that traumatic acid was decreased in OVX mice and increased under nodakenin treatment.It is suggested that the alteration in these obesity-related metabolites may be the key to nodakenin's ability to inhibit OVX-induced weight gain.In addition, KEGG enrichment analysis revealed that of the pathways in which nodakenin alleviated the significant changes induced by OP, four of them were closely associated with osteoporosis, including endocrine and other factor-regulated calcium reabsorption; parathyroid hormone synthesis, secretion, and action; thyroid hormone synthesis and mineral absorption.Moreover, calcitriol, as a hub of metabolites is closely related to these four pathways.It plays a crucial role in bone regeneration by binding to the vitamin D receptor (VDR) to activate various downstream transcription factors involved in signaling pathways such as calcium reabsorption, hormone synthesis and secretion [53][54][55][56].Consistent with this, our experiments showed that VDR expression was elevated in the bone tissue of nodakenin-treated mice.
Estrogen deficiency is tightly linked to the disruption of the intestinal mucosal barrier, which adversely affects skeletal properties through the modulation of the GM-bone axis [57,58].PMOP is strongly associated with systemic chronic inflammation, which is often accompanied by the augmentation of pro-osteoclastogenic factors such as TNF-α and IL-1β [59,60].Interestingly, nodakenin-treated mice exhibit reduced levels of expression of the pro-osteoclastogenic factors TNF-α and IL-1β in the intestine, which is likely an important reason for their inhibition of osteoclast formation.Increased intestinal permeability has been shown to play a key role in the regulation of inflammation and is inextricably linked to the colonization of Firmicutes as well as the diversity of gut microbes [33].Our study provides strong support for this view, as nondakenin treatment significantly reduced intestinal permeability in OVX mice, accompanied by a decrease in the release of proinflammatory cytokines IL-1β and TNF-α.This phenomenon is likely related to the role of nondakenin in enhancing gut microbial diversity and reducing the relative abundance of Firmicutes.Notably, calcitriol protects the barrier function of endothelial and epithelial cells and effectively inhibits intestinal inflammation [61].In our study, nodakenin-treated mice exhibited a significant rise in calcitriol, suggesting that the reduction in intestinal permeability and expression of pro-inflammatory factors by nodakenin may be related to calcitriol.This finding not only deepens our understanding of the role of nodakenin in maintaining intestinal health but also provides new ideas for future therapies targeting bone-related diseases.

Animal Experiment and Design
We used C57BL/6J mice (6-8 weeks); one week after acclimatization, six mice were randomly selected as controls (Sham) for cesarean section and the rest of the mice received bilateral ovariectomy.In detail, after anesthesia with 1% sodium pentobarbital, the dorsal test fur was shaved and disinfected, and then, the bilateral ovaries were removed and the oviducts were ligated.After surgery, OVX mice were randomly divided into two groups (n = 6 in each group) and the vehicle (intragastric administration) (1% DMSO + 5% tween 80 + saline) and nodakenin (Chemfaces, Wuhan, China) (20 mg/kg nodakenin + 5% tween 80 + saline) were administered once daily for 6 weeks.Body weights were measured at the end of the experiment, and serum, feces, femurs, and viscera were collected for subsequent experimental analysis.The Qingdao Agricultural University granted approval for all the animal experimentation processes.

Micro-CT (µCT) Analysis
Mouse femurs were dissected and then analyzed using a SkyScan-1276 X-ray microtomograph (Micro-CT) (Bruker, Billerica, MA, Germany).Using the bottom-most portion of the lateral growth plate of the femoral knee as the baseline, the region of the bone marrow cavity, with a thickness of 1 mm underneath was selected and designated as the region of interest (ROI) for 3D reconstruction; 3D images were reconstructed using N-Recon software and analyzed in three dimensions using CT-AN software.

Histomorphometric Analysis
The femurs of mice, with the exception of muscle tissue were immersed in a 4% paraformaldehyde solution for 24 h.Subsequently, they were transferred to an EDTA decalcification solution for an additional 4 weeks.Paraffin-embedded sections were cut into 4 µm thick sections and stained with TRAP.The sections were observed under a microscope, and images were captured.
Sections of paraffin-embedded colon tissue (5 µm) were stained with hematoxylin and eosin to assess intestinal barrier integrity.The morphological characteristics of intestinal villi, including the intestinal villus height (µm), intestinal crypt depth (µm), and the ratio of villus height to crypt depth (V/C) were assessed using Image-Pro Plus.

Quantitative Real-Time PCR
Mouse colon RNA was prepared by using TriZol reagent (Invitrogen, Waltham, MA, USA) and then quantified using a Nano Drop 2000 (Thermo Fisher, Waltham, MA, USA).The reverse transcription of cDNA was performed using the PrimeScript RT Reagent Kit with gDNA Eraser (Takara, Kyoto, Japan).A real-time polymerization chain reaction was performed using a CFX96 Real-Time Fluorescence Quantitative PCR System (Bio-Rad, Hercules, CA, USA).GAPDH was used as an internal reference, and the primers used are shown in Supplementary Table S1.The 2 −∆∆CT method was used for analysis.

16S rRNA Analysis
Fresh mouse feces were collected and frozen at −80 • C (n = 6).Total genomic DNA was extracted from the fecal samples and assayed for quality using a Nanodrop 2000 and gel electrophoresis.16rRNA gene V3-V4 region amplification was performed using the Bacterial Universal Priming Kit (Omega Bio-tek, Norcross, GA, USA), and paired-end sequencing (2 × 250 bp) was carried out with the Illumina MiSeq platform (Illumina, San Diego, CA, USA) by Beijing Allwegene Tech, Ltd. (Beijing, China).Based on the OTUs and richness results, the Alpha Diversity Index and Beta Diversity Distance Matrix were calculated using QIIME v1.8.0 software [62].The Bray-Curtis distance was measured using the R package (R version 3.6.0)to compare the differences in microbial composition between groups.Linear discriminant analysis (LDA), combined with effect size (LEfSe, Blair, WI, USA, Galaxy version 1.0) was used to identify key OTUs that responded to different interventions.

Metabolic Analysis of Feces by LC-MS/MS
The samples were removed from the freezer at −80 • C and thawed slowly at 4 • C.An appropriate amount of feces (n = 6) was added into pre-cooled MeOH:ACN:H 2 O (v:v:v = 2:2:1) solution and ground and sonicated in order to homogenize the sample; the sample was allowed to stand at −20 • C for 1 h and then centrifuged at 13,000 rpm for 15 min at 4 • C. The supernatant was then freeze-dried.For mass spectrometry analysis, an appropriate amount of ACN:H 2 O solution (v:v = 1:1) was added for vortexing and sonication to make it fully re-dissolve, and then centrifuged at 13,000 rpm 4 • C for 15 min; the supernatant was aspirated into the injection bottle and then analyzed by LC-MS/MS.
Metabolomics was used for separation using a Waters ACQUITY UPLC HSS T3 column (2.1 mm × 100 mm, 1.8 µm, Milford, MA, USA).The mobile phases were A (0.1% formic acid in water) and B (0.1% formic acid in acetonitrile) for gradient elution: 95% B, 0.5 min; 95-75% B, 0.5-2.5 min; 70-100% B, 2.5-7.5 min; 100% B, 7.5-9 min; 100-5% B, 9-9.5 min and 5% B, 9-12 min.The injection volume was 2 µL, the flow rate was 0.3 mL/min, and the column temperature was 50 • C. The sample was collected in the positive and negative ion modes.Sample metabolic analytes exiting the column were collected once in positive ion mode and once in negative ion mode using a high-resolution mass spectrometer, Triple TOF 5600+(SCIEX, Washington, DC, USA).The details of the parameters are as follows: Ion Source Gas: 150 psi; Ion Source Gas2: 50 psi; Curtain Gas: 35 psi; Source Temperature: 500 • C; IonSapary Voltage Floating:5500 V and −4500 V (positive and negative); Declustering Potential (DP): ±80 V (positive and negative); TOF MS scan m/z range:60-1200 Da; Product Ion Scan m/z Range: 25-1200 Da; TOF MS Scan Accumulation Time: 0.25 s/spectra; Product Ion Scan Accumulation Time: 0.03 s/spectra.The secondary mass spectra were obtained by Information Dependent Acquisition (IDA), and the extracted data were analyzed for metabolite structure identification and data analysis using the High Sensitivity mode, CE: 30 V ± 15, to screen for different metabolites in each group.
Peak identification, peak extraction, peak alignment, and integration were performed using an R package written in-house (kernel XCMS), followed by metabolite identification and screening using the in-house MS2 database (Allwegene, Beijing, China).The differential metabolic pathways were screened against authoritative metabolite databases such as KEGG and PubChem.

Statistics
All data are presented as the mean ± standard deviation.t-tests were employed for comparisons between two groups, and one-way ANOVA was used for multiple comparisons.Spearman's correlation analysis was used to analyze the relationship between the gut microbiota, differential metabolites, and bone metabolic parameters.Heat map correlation statistics were analyzed using SPSS 25.0 (IBM).* p < 0.05; ** p < 0.01.

Conclusions
In summary, the present study confirmed that nodakenin modulates gut microbiota composition and gut metabolite disorders and restores intestinal permeability.Nodakenin treatment decreased the F/B ratio, as well as produced some specific gut microbiota changes at the genus level.Meanwhile, nodakenin may affect hormone synthesis and calcium reabsorption by regulating metabolites such as calcitriol, thereby ameliorating osteoporosis.The present study provides new perspectives for exploring nodakenin's anti-osteoporotic effects and lays the foundation for considering gut microbiota as a drug target.

Figure 3 .
Figure 3. Clustering and diversity analysis of gut microbiota in different groups of mice.Changes in fecal α-diversity of gut microorganisms: Chao 1 estimator (A) and Shannon index (B).(C) Plot-based UniFrac PCoA.(D) Cluster analysis of gut microbiota among samples.Mean ± SEM (n = 6).* p < 0.05, ns means no significant difference.

Figure 4 .
Figure 4. Relative abundance analysis of the gut microbiota at the multispecies level for different groups.(A) Composition of gut microbiota at the phylum level.(B) Gut microbiota are significantly different at the phylum level.(C) Abundance ratio of F/B.(D) Composition of gut microbiota at the phylum level.(E) Gut microbiota significantly differ at the genus level.(F) Histogram of the distribution of LDA values.(G) Evolutionary branching diagram of LEfSe analysis.Mean ± SEM (n = 6).* p < 0.05; ** p < 0.01.Spearman correlation analysis was used to evaluate the associations between alterations in the gut microbiota composition and the markers associated with osteoporosis.At the phylum level, Patescibacteria and Firmicutes exhibited a positive correlation with RANKL, TRAP, and SMI, while displaying a negative correlation with Tb.N, BMD, BV/TV, BGP, and B-ALP.Bacteroidota showed a positive correlation with Tb.N, BMD,

Figure 4 .
Figure 4. Relative abundance analysis of the gut microbiota at the multispecies level for different groups.(A) Composition of gut microbiota at the phylum level.(B) Gut microbiota are significantly different at the phylum level.(C) Abundance ratio of F/B.(D) Composition of gut microbiota at the phylum level.(E) Gut microbiota significantly differ at the genus level.(F) Histogram of the distribution of LDA values.(G) Evolutionary branching diagram of LEfSe analysis.Mean ± SEM (n = 6).* p < 0.05; ** p < 0.01.

Molecules 2024 ,
29,  x FOR PEER REVIEW 7 of 18 BV/TV, and bone formation-associated BTMs, and negatively correlated with TRAP, SMI, and RANKL.Parasutterella was positively correlated with Tb.N, BMD, BV/TV, and B-ALP.Muribaculum and Bacteroides were positively correlated with RANKL, TRAP, and SMI.Muribaculum was negatively correlated with BV/TV and B-ALP.Bacteroides was negatively correlated with BV/TV, Tb.N, and BGP.Roseburia was negatively correlated with B-ALP and positively correlated with TRAP (Figure5B).Muribaculum was negatively correlated with BV/TV and B-ALP.Bacteroides had a negative correlation with BV/TV, Tb.N, and BGP.Roseburia had a negative correlation with B-ALP and a substantial positive correlation with TRAP (Figure5B).

Figure 5 .
Figure 5. Analysis of correlations between gut microbiota and bone turnover indicators.(A) Heat map of Spearman's r correlation between gut microbiota and bone biological parameters at the phylum level.(B) Heat map of Spearman's r correlation between gut microbiota and bone biological parameters at the genus level.*p < 0.05; ** p < 0.01.

Figure 5 .
Figure 5. Analysis of correlations between gut microbiota and bone turnover indicators.(A) Heat map of Spearman's r correlation between gut microbiota and bone biological parameters at the phylum level.(B) Heat map of Spearman's r correlation between gut microbiota and bone biological parameters at the genus level.*p < 0.05; ** p < 0.01.

Molecules 2024 , 18 Figure 6 .
Figure 6.Composition and differences in gut metabolites and analysis of metabolic pathways.(A) OPLS-DA analysis of the OVX group versus the nodakenin-treated group.(B) Volcano plot of differential metabolites.(C) VIP value plot of differential metabolites (top 20 VIP values).(D) KEGG differential metabolic pathway plot (top 20 metabolic pathway hits).

Figure 6 .
Figure 6.Composition and differences in gut metabolites and analysis of metabolic pathways.(A) OPLS-DA analysis of the OVX group versus the nodakenin-treated group.(B) Volcano plot of differential metabolites.(C) VIP value plot of differential metabolites (top 20 VIP values).(D) KEGG differential metabolic pathway plot (top 20 metabolic pathway hits).

Figure S2 :
Rarefaction curves and Shannon-Wiener curves; Table S1: Primer sequence for qRT-PCR.Author Contributions: Conceptualization, G.L. and C.L.; methodology, G.L. and C.L.; software, J.C. and Z.W.; validation, C.L., Z.W. and Y.L.; formal analysis, C.L., J.C. and Z.W.; investigation, C.L., J.C., Z.W. and Y.Z.; writing-original draft preparation, C.L.; writing-review and editing, G.L. and C.L.; project administration, G.L. and C.L.; funding acquisition, G.L. All authors have read and agreed to the published version of the manuscript.Funding: This study was funded by the Start-up Fund for Scientific Research of High-level Talents of Qingdao Agricultural University (1121021) and Shandong Modern Agricultural Technology & Industry System (SDAIT-21-01) to Guangyu Li.Institutional Review Board Statement: All animal care and experimental procedures in this study were approved by the Board of Animal Study of Qingdao Agricultural University (protocol code 2023032701, approval date: 27 March 2023).