- Open Access
Metabolomics and physiological analysis of the effect of calcium supplements on reducing bone loss in ovariectomized rats by increasing estradiol levels
Nutrition & Metabolism volume 18, Article number: 76 (2021)
Data from the 2010–2012 Chinese National Nutrition and Health Survey showed that the vast majority of postmenopausal women in China had dual deficiencies in calcium and estrogen.
This study aimed to clarify whether calcium supplementation alleviated bone loss caused by calcium restriction combined with estrogen deficiency in rats.
Forty-eight female rats aged 9 weeks were assigned to 4 groups and fed a low-calcium diet: sham-operated (SHAM-LC), ovariectomized (OVX-LC), and ovariectomized rats treated with 750 mg/kg (OVX-LC-M) or 2800 mg/kg CaCO3 (OVX-LC-H). CaCO3 or distilled water was administered orally for 13 weeks. Bone mineral density (BMD) and histomorphometry of the femur, serum biochemical parameters, and serum metabolites were analyzed.
The OVX-LC rats showed a significant increase in body weight and serum levels of lipid markers, a significant decrease in serum estradiol, calcium, phosphorus, and 25(OH)D levels, and deterioration of the femur. At 750 mg/kg and 2800 mg/kg, CaCO3 reduced the deterioration of trabecular bone and increased the trabecular area percentage (Tb.Ar %) and BMD of the femur. Serum estradiol levels increased in a dose-dependent manner after CaCO3 supplementation (p < 0.01). The administration of 2800 mg/kg CaCO3 decreased serum triglyceride and high-density lipoprotein levels (p < 0.05) and decreased the levels of the bone turnover markers osteocalcin, N-telopeptide of type I collagen and β-crosslaps. The results of the metabolomics analysis showed that the glycerophospholipid metabolism pathway was closely related to calcium supplementation, and more DG (44:6 n3), LysoPC (22:2) and PE (P-34:3) and less Cer (d43:0) and PE-NMe2 (46:3) were produced.
The results clearly indicated that calcium supplementation was beneficial for decreasing bone loss in OVX-LC rats. The present study is the first to show that calcium supplementation increased the estradiol content in OVX-LC rats, and the effect of calcium on bone loss may be partially attributed to the increase in the estrogen level that subsequently induced the changes in metabolite levels, eventually increasing the bone mineral density to a relatively higher level to reduce bone deterioration.
Osteoporosis is a disease characterized by a decrease in bone mineral density and the destruction of the bone microstructure, which eventually leads to increased bone fragility and fracture [1,2,3]. With the aging of the population, osteoporosis has become a common disease that endangers the health of elderly individuals worldwide and has imposed heavy economic and social burdens . According to an epidemiological survey of osteoporosis in China in 2018, the incidence of osteoporosis in people over 50 years old was 19.2%, and the incidence was even higher at 30% and 50% for females over 50 and 65 years old, respectively .
Many factors contribute to osteoporosis, such as the insufficient intake of calcium (Ca) and protein, a low level of vitamin D, endocrine factors, heredity, and living habits [6,7,8,9,10,11]. For postmenopausal women, osteoporosis is mainly caused by estrogen deficiency. The 2010–2012 Chinese National Nutrition and Health Survey showed that the daily Ca intake of Chinese residents was 366.1 mg/d, which was only half of the recommended Ca intake, and the Ca intake of 96.6% of Chinese residents was lower than the estimated average requirement for Ca . Therefore, the vast majority of postmenopausal women in China are present with dual deficiencies in Ca and estrogen.
Ca is an important component of bone, and > 99% of Ca in the body is deposited in bone and teeth . Ca supplementation has been widely used to prevent osteoporosis and subsequent fractures in postmenopausal women based on the hypothesis that an adequate Ca level is crucial for maintaining bone health [14, 15]. This hypothesis has also been supported by many studies. A meta-analysis showed that Ca intake effectively postpones the decreasing trend in the BMD and reduces the risk of fractures in postmenopausal women [16,17,18,19]. However, the effect of Ca on preventing postmenopausal osteoporosis was not consistent in recent data. Reid reported that the effect of calcium supplementation on fractures in healthy older women remained uncertain , and some meta-analyses showed that the evidence of the effectiveness of Ca treatment for fractures in postmenopausal women remained limited .
Many explanations for the controversial results of meta-analyses of population studies have been proposed, such as different dietary Ca intake levels, national Ca recommendations, vitamin D intake statuses, Ca supplement dosages and durations, and poor long-term compliance. A rat bone loss model was induced by performing OVX and administering a low Ca diet to clarify whether Ca supplementation is beneficial for the prevention of osteoporosis induced by estrogen and Ca deficiencies. Then, different concentrations of calcium carbonate (CaCO3) were administered to rats to observe the effects of Ca supplementation on bone loss and explore the possible mechanism.
Metabolomics can be used to collect information on metabolites, including lipids, amino acids, sugars and vitamins, in blood, urine and tissues, and changes in metabolites are sensitive indicators of nutrition and metabolism [22,23,24]. Metabolomics has been widely used to identify biomarkers and explore the mechanism of osteoporosis [25,26,27]. Here, the mechanism underlying the effect of Ca on osteoporosis was explored with metabolomics after rats were supplemented with different concentrations of Ca for 13 weeks.
CaCO3 was obtained from the Harbin Pharmaceutical Group, China (Batch No. 20190227). Forty-eight 7- to 8-week-old female Sprague–Dawley rats with body weights of 250–300 g were purchased from Beijing WeitongLihua Experimental Animal Technology Co. Ltd., license number: SCXK (Peking) 2016-0011. The rats were maintained on a standard 12 h light/12 h dark illumination cycle with water and chow provided ad libitum.
After 7 days of acclimation, all rats were anesthetized with an intraperitoneal injection of pentobarbital sodium (30 mg/kg BW). All rats underwent either bilateral OVX (n = 36) or sham operation (SO, n = 12). For the sham-operated rats, bilateral abdominal incisions and sutures were performed without oophorectomy . Uterine atrophy (observed during dissection) indicated a successful operation in all OVX rats. After 7 days of convalescence, all rats were assigned to 4 groups and fed a low Ca diet: sham-operated rats fed a low Ca diet (SHAM-LC, n = 12), OVX rats fed a low Ca diet (OVX-LC, n = 12), OVX rats fed a low Ca diet and treated with 750 mg/kg CaCO3 (OVX-LC-M, n = 12) or OVX rats fed a low Ca diet and treated with 2800 mg/kg CaCO3 (OVX-LC-H, n = 12). The rats in the SHAM-LC and OVX-LC groups were administered distilled water. CaCO3 or distilled water was administered to rats by oral gavage for 13 weeks. The experimental protocol was approved by the Animal Ethics Committee of National Institute of Nutrition and Health, Chinese Center for Diseases Control and Prevention.
Preparation of the low Ca diet
The low Ca diet was prepared according to the “Test method for improving bone mineral density function” inspection and assessment standard for health food issued by the Ministry of Health, People’s Republic of China. The diet was composed of 32% corn starch, 30% sucrose, 23% casein, 5% fiber, 5% corn oil, 3.5% mineral mixture, 1% vitamin mixture, 0.3% DL-methionine, and 0.2% choline tartrate. The following mineral mixture (per kilogram diet) was used: MnSO4 110 mg, CuSO4 0.8 mg, FeSO4 1.2 mg, KI 18.0 mg, ZnSO4 2960 mg, CaHPO4 2890 mg, and MgSO4 12.5 g. The following vitamin mixture was used (per kilogram diet): vitamin A 1.4 × 104 IU, vitamin D 1500 IU, vitamin E 120 mg, vitamin K 3 mg, vitamin B1 12 mg, vitamin B2 20 mg, vitamin B6 12 mg, vitamin B12 0.03 mg, nicotinic acid 60 mg, pantothenic acid 24 mg, folate 6 mg, and biotin 0.54 mg. The Ca content in the low Ca diet was 2.2 g/kg diet after examination.
Analyses of serum parameters
At the termination of the study, all rats were fasted overnight, and abdominal aortic blood was collected. The levels of Ca, phosphorus (P), glucose (GLU), total cholesterol (TC), triglycerides (TGs), high-density lipoprotein (HDL), low-density lipoprotein (LDL), and nonesterified fatty acids (NEFAs) in serum were examined with an automatic biochemistry analyzer (Hitachi 7600, Japan). The content of 25(OH)D in serum was determined with an ultra-performance liquid chromatography-tandem mass spectrometer (UPLC-MS/MS, Shimadzu 8060, Japan). Serum estradiol (E2) concentrations were assayed with an automatic immune analyzer (DXI800, Beckman, USA). Bone turnover markers (BTMs), including procollagen I N-terminal peptide (PINP), osteocalcin (OC), N-telopeptide of type I collagen (NTX), and β-crosslaps (β-CTX), were detected with rat ELISA kits (Cusabio Biotech Co., China).
Femoral BMD analysis
The left femur was extracted and examined with a dual-energy X-ray absorptiometry system (DXA, Hologic, USA). The BMD of the whole femur and distal end of femur was analyzed. The placement positions of each femur were consistent. The region of interest (ROI) was defined manually after the scouting scan. The densitometer was calibrated using small animal quality assurance phantoms provided by the manufacturer and performed using established procedures before the series of measurements was collected.
Histomorphometric analysis of the femur
After fixation with 10% buffered formalin for 5 days, the right femurs were decalcified with 10% EDTA for 3 days. Then, the distal metaphyses of femurs were dehydrated in ethanol, defatted in xylene, embedded in paraffin, and sliced into longitudinal sections (5-μm thick). H&E staining was performed on the sections. The morphology of the rat trabecular bone in the distal femoral metaphyses was observed, and the static parameter of Tb.Ar % within 3 mm under the epiphyseal plate were calculated with Image-Pro Plus 6.0 software.
The metabolomics method described in previous studies was used [29, 30]. Serum (100 μL) was mixed with 0.9 mL of 80% methanol containing 0.1% formic acid (FA). After vortexing for 30 s and ultrasonication for 20 min, all samples were frozen for 1 h at − 20 °C for protein precipitation. Then, the samples were centrifuged at 12,000 × g for 10 min at 4 °C, and the supernatant (800 μL) of each sample was collected and transferred to a sample vial. Quality control (QC) samples were prepared by pooling aliquots of all serum samples for the serum metabolomics analysis. Blank samples (80% methanol containing 0.1% formic acid) and QC samples were repeated every ten samples during data acquisition. The sample vials were stored at − 20 °C until detection.
The UPLC-QTOF MS analysis was performed using a UPLC system (ACQUITY UPLC I-Class, Waters) coupled to an electrospray ionization quadruple time-of-flight mass spectrometer (ESI-QTOF MS) (SYNAPT G2-Si HDMS, Waters). A Waters ACQUITY BEH C18 column (1.7 μm; 100 mm × 2.1 mm) was used for LC separation, and the column temperature was maintained at 40 °C. The flow rate was 0.4 mL/min, and the sample injection volume was 10 μL. Mobile phase A was 0.1% FA in water, and mobile phase B was 0.1% FA in ACN. The following linear gradient was set: Initial to 1 min: 10% B, 0–3 min 10% B to 80% B, 3–8 min: 95% B to 95% B, and 8.1–10 min: 10% B.
High-accuracy MS data were recorded using MassLynx 4.1 software. The capillary voltage was 2.5 kV for positive mode and 2 kV for negative mode, whereas the cone voltage was 30 V for both modes. The source temperature was set to 120 °C with a cone gas flow of 50 L/h, and the desolvation temperature was set to 500 °C with a desolvation gas flow of 800 L/h. Leucine-enkephalin (LE) was used as the lock mass, generating a reference ion at m/z 556.2771 in positive mode and m/z 554.2615 in negative mode, which was introduced by a lockspray at a rate of 10 μL/min for data calibration. The MSE data were acquired in continuum mode using ramp collision energy in two scan functions. For low-energy mode, a scan range of 50–1200 Da, scan time of 0.2 s, and collision energy of 6 V were used. For high-energy mode, a scan range of 50–1200 Da, scan time of 0.2 s, and collision energy ramp of 15–45 V were used.
Raw data were imported into the commercial software Progenesis QI (Version 2.4, Waters) for data processing, which included peak selection, peak alignment and acquiring compound-associated information such as the m/z, retention time and intensity. Next, data filtering was performed to delete low-quality data. Ions with a relative standard deviation (RSD) of more than 30% in QC samples were filtered. These filtered ions fluctuated substantially among samples and were not included in further analyses. PLS-DA (partial least squares discriminant analysis) was performed, and VIP (variable importance in projection) was calculated using MetaboAnalyst 4.0 (https://www.metaboanalyst.ca/)  and R project software. R project software was also applied in further data processing and statistical analyses. The pathway analysis was performed using MetaboAnalyst 4.0 software. The column graph was drawn using GraphPad Prism 8.0 software.
All data are presented as the means ± SD. The differences between groups were analyzed using the T-test or one-way ANOVA with SPSS 16.0 software. Differences were considered statistically significant at p < 0.05.
Effect of CaCO3 on body weight
As shown in Fig. 1, all rats in the four groups had similar initial body weights. The body weight of rats in the OVX-LC group increased significantly after surgery compared to that of the SHAM-LC group (p < 0.01). A tendency of reduced body weight after CaCO3 treatment was observed, although a significant difference was not observed among the OVX-LC, OVX-LC-M and OVX-LC-H groups. The results indicated that Ca supplementation slightly inhibited the weight gain induced by OVX-LC.
CaCO3 significantly increased the E2 content
The induction of postmenopausal conditions by OVX resulted in a dramatic reduction in E2 levels. Compared to the SHAM-LC group, serum E2 levels in the OVX-LC group were reduced significantly (p < 0.01) by approximately threefold. Interestingly, serum E2 levels increased significantly in the OVX-LC-M and OVX-LC-H groups compared to the OVX-LC group (p < 0.01), and a significant positive dose–response relationship between the dose of CaCO3 and serum E2 level was observed (Fig. 2). Ca supplementation increased the E2 concentration in OVX-LC rats, which might be an important factor to alleviate osteoporosis.
Effect of CaCO3 on biochemical parameters
As shown in Table 1, compared to the SHAM-LC group, the serum Ca, P, and 25(OH)D concentrations decreased significantly, and the serum TC, HDL and LDL concentrations increased significantly in the OVX-LC group (p < 0.05 or p < 0.01). The results indicated that the induction of postmenopausal conditions by OVX led to significant disorders of bone metabolism and lipid metabolism parameters. Compared to the OVX-LC group, the TG and HDL contents in the OVX-LC-H group decreased significantly (p < 0.05), and an obvious dose–response relationship was observed among the three groups. No obvious differences were observed in the other parameters examined in serum among the OVX-LC, OVX-LC-M and OVX-LC-H groups. Based on the results, Ca supplementation partially counteracted the disorder of lipid metabolism induced by OVX-LC.
Effect of CaCO3 on BTMs
As shown in Fig. 3, compared to the SHAM-LC group, no obvious differences were observed in the four examined BTMs in the OVX-LC group. Compared to the OVX-LC group, the content of OC in the OVX-LC-M and OVX-LC-H groups decreased significantly (p < 0.01), and the contents of NTX and β-CTX in the OVX-LC-H group decreased significantly (p < 0.05).
Effect of CaCO3 on the BMD of the femur
As shown in Fig. 4, compared to the SHAM-LC group, OVX caused a significant decrease in the BMDs of the whole femur and distal end of the femur in the OVX-LC group (p < 0.01). The BMDs of the whole femur and distal end of the femur in the OVX-LC-M and OVX-LC-H groups increased significantly compared to those of the OVX-LC group (p < 0.05). These results confirmed that the rat bone loss model was successfully established by the OVX operation, while Ca supplementation partially reversed the decrease in BMD in the OVX-LC rats.
Effect of CaCO3 on the histomorphometry of the femur
As shown in Figs. 5 and 6, the structure of trabecular bone in the OVX-LC group became sparse, slender and fractured compared to that of the SHAM-LC group. The changes induced by OVX were accompanied by a significant decrease in the Tb.Ar % (p < 0.01). After 13 weeks of CaCO3 treatment, improved continuity, integrity and numbers of trabecular bone were observed in the OVX-LC-M and OVX-LC-H groups. Tb.Ar % also increased significantly in the OVX-LC-M and OVX-LC-H groups compared to the OVX-LC group (p < 0.05 or p < 0.01). Therefore, Ca supplementation improved the structure of the bulk trabecular bone, thereby exerting a protective effect on the bones of OVX-LC rats.
Effect of CaCO3 on serum metabolites
PLS-DA model analysis
Nine hundred sixty seven metabolites in rat serum were identified using metabolomics, including glycerophospholipids, triglycerides, diglycerides, ceramides, organic acids, amino acids, fatty acids, and vitamins. According to the PLS-DA multivariate statistical model analysis, obvious differences were discovered on the score graphs among the OVX-LC, OVX-LC-M and OVX-LC-H groups, yet little difference was detected between the OVX-LC-M and OVX-LC-H groups (Fig. 7). The results indicated that calcium supplementation induced changes in metabolites, but no significant difference was observed after the administration of higher concentrations of the Ca supplement.
We used MetaboAnalyst 4.0 software to analyze and compare the metabolic pathways among the OVX-LC, OVX-LC-M and OVX-LC-H groups, and found that the glycerophospholipid and butanoate metabolic pathways play a major role in the effect of Ca supplementation on OVX rats. However, no significant difference was observed in the butanoate metabolic pathway (p values were 0.76, 0.23, and 0.8 for comparisons between the OVX-LC and OVX-LC-M groups, OVX-LC and OVX-LC-H groups, and OVX-LC-M and OVX-LC-H groups, respectively), but a significant difference in the glycerophospholipid metabolic pathway was identified (p values were 5.12E−7, 4.08E−8, and 0.26 for comparisons between the OVX-LC and OVX-LC-M groups, OVX-LC and OVX-LC-H groups, and OVX-LC-M and OVX-LC-H groups, respectively). The main differences in glycerophospholipid metabolism were identified in phosphatidylcholines, phosphatidylethanolamines and ceramides (Fig. 8).
Cross analysis of differentially altered metabolites
Through pairwise comparisons, numerous differentially altered metabolites were identified, including 239 metabolites between the OVX-LC and OVX-LC-M groups, 365 metabolites between the OVX-LC and OVX-LC-H groups, 86 metabolites between the OVX-LC-M and OVX-LC-H groups, and 46 metabolites in common among the three groups, as shown in the center of Fig. 9A. The results of the heatmap analysis of the 46 metabolites showed that 6 of them were decreased significantly and 40 of them were increased significantly after Ca supplementation (shown in Fig. 9B).
Correlation analysis of Ca and estrogen with differentially altered metabolites
The Pearson method was used to analyze the correlations of 46 differentially altered metabolites with Ca and E2 levels. Among the 46 metabolites, 10 were significantly correlated with Ca and E2 levels, and the correlation coefficient was greater than 0.6. The correlation diagram showed that 5 metabolites were positively correlated with Ca supplementation, including DG (44:6n3), LysoPC (22:2) and PE (P-34:3), while the other 5 metabolites were negatively correlated with Ca supplementation, including Cer (d43:0) and PE-NMe2 (46:3), as shown in Fig. 10. This result indicated that Ca supplementation had a close correlation with changes in the levels of some metabolites, especially changes in glycerophospholipid levels.
Ca is an important micronutrient that regulates various physiological activities in the human body, such as the formation of bone and teeth, and maintains the normal function of cells . Sufficient dietary Ca intake is important for maintaining bone health. When Ca intake is insufficient, the organism is forced to increase the process of osteolysis to maintain Ca homeostasis . In the present study, a low Ca diet (22% Ca recommended dose) was used to induce a severe Ca deficiency. Two doses of CaCO3 were administered to rats to produce a moderate deficiency and a sufficient Ca intake level to observe the effect of Ca on bone loss induced by OVX-LC.
Estrogen deficiency is the major factor contributing to postmenopausal osteoporosis. Many similarities have been observed between OVX-induced bone loss in rats and postmenopausal osteoporosis in humans. Thus, OVX rats have been widely used as an animal model to study postmenopausal osteoporosis [34, 35]. In the present study, a rat bone loss model was established by OVX-LC. A BMD measurement is considered the standard method for diagnosing osteoporosis in humans, and DXA has been used to measure BMD in rats in many studies [36, 37]. As expected, a significant decrease in the BMD of the femur was observed in the OVX-LC group compared to the SHAM-LC group, which confirmed that estrogen is important for maintaining bone density. The BMD of the whole femur and distal end of the femur in the OVX-LC-M and OVX-LC-H groups increased significantly after 13 weeks of Ca supplementation compared to the values of the OVX-LC group. The changes in the BMD of the femur indicated that Ca supplementation is beneficial to counteract the decrease in bone loss induced by OVX-LC. Gao also reported that a low Ca diet (0.1% calcium) accelerates bone loss in OVX rats , consistent with the results of our study.
Deterioration of trabecular bone has been implicated in the incidence of osteoporosis in humans. Thus, an assessment of trabecular bone structure is important to evaluate the quality of the femur in addition to BMD. In this study, a histomorphometric analysis of the femur was performed at the end of the study. The trabecular structure of the femur deteriorated obviously, and the Tb.Ar % decreased significantly compared to that of the SHAM-LC group. The deterioration of trabecular bone was induced by estrogen deficiency, which has also been reported in other studies. After 13 weeks of Ca supplementation, the deterioration of trabecular bone of the femur was ameliorated, and the Tb.Ar % increased significantly in the OVX-LC-M and OVX-LC-H groups. Based on these findings, Ca supplementation contributed to the restoration of trabecular bone.
Decreases in the BMD and deterioration of trabecular bone of the femur were linked to an imbalance in bone remodeling. BTMs are enzymes associated with osteoblasts and osteoclasts [39, 40]. They are intermediate products in the process of bone remodeling. The contents of BTMs in serum or urine reflect changes in bone remodeling before changes in BMD, and BTMs have been suggested as independent risk factors for osteoporotic fractures . Many studies have reported a close association between an increase in the levels some BTMs and an increased risk of hip fracture in the population [42, 43]. For postmenopausal women, the BMD decreases markedly and BTMs increase significantly upon the withdrawal of estrogen [44,45,46]. In our study, BTMs, including ALP, PINP, OC, β-CTX and NTX, were detected to understand the process of bone formation and bone resorption. Compared to the SHAM-LC group, no significant changes were observed in the OVX-LC group. The changes in BTMs induced by OVX in rats were not consistent with those documented in previous studies [47, 48]. Yan Zhang  reported increased urinary deoxypyridinoline (one biomarker of bone resorption) levels at the 4th week after OVX followed by a decrease at the 18th week, and serum ALP (one biomarker of bone formation) levels decreased early and increased late after OVX, indicating that the changes in BTMs after OVX were time-dependent. On the other hand, 2–3-month-old OVX rats were used in the present study. N Patlas reported that rats at 1 and 3 months were more suitable for research on bone histomorphometric parameters . However, continuous physiological bone growth occurs in rats of this age, which might lead to an increase in BTMs. In the present study, serum OC, NTX and β-CTX concentrations in the OVX-LC-H group decreased significantly compared to those in the OVX-LC group, indicating that Ca supplementation slowed the bone remodeling process.
The results of the metabolomics analysis revealed 46 differentially altered metabolites shared among the three groups after Ca supplementation, 6 of which decreased significantly and 40 of which increased significantly. A strong correlation was observed between Ca supplementation and 10 metabolites (correlation coefficient > 0.6), 5 of which had a positive correlation and 5 of which had a negative correlation. These metabolites were mainly glycerophospholipids, including Cer (d43:0), PE-NMe2 (46:3), DG (44:6n3), LysoPC (22:2), and PE (P-34:3). Jiaqi W also reported that one anti-osteoporosis medicine alleviates osteoporosis by regulating glycerophospholipids , which is similar to our results showing that glycerophospholipids play an important role in osteoporosis. In addition, studies have also reported that changes in glycerophospholipids are regulated by human hormones, such as estradiol, and glycerophospholipid metabolism is the main potential target pathway of E2 .
Estrogen plays an important role in maintaining the balance of bone absorption and bone formation, and estrogen deficiency is known to induce osteoporosis in a variety of animals and in humans . In humans, estrogen replacement appears to be the most efficient method to alleviate postmenopausal osteoporosis regardless of its side effects . In animal studies, estrogen supplementation can also successfully prevent the reduction of BMD in OVX rats [55, 56]. In our study, the serum E2 content in rats increased significantly after Ca supplementation, and an obvious dose–response relationship between the E2 level and Ca content was observed among the OVX-LC, OVX-LC-M and OVX-LC-H groups. Researchers also reported that electroacupuncture and traditional Chinese medicine induce an obvious increase in E2 levels in OVX rats [57, 58]. Since the rats used in our study were OVX rats, E2 was unable to be produced by the ovary. Evidence from recent studies revealed that E2 is synthesized by many nongonadal sites, such as bones, neurons, pancreas, adipose tissue, mesenteric lymph nodes and Peyer’s patches [12, 59,60,61,62]. For OVX rats, the capacity of nongonadal sites to synthesize estradiol may be the determining factor in the eventual increase in the BMD.
The results of biochemical analyses in our study showed significant decreases in TG and HDL levels after Ca supplementation. For postmenopausal women, notable changes in blood lipid profiles have been observed upon the withdrawal of estrogen. The results of a meta-analysis showed higher serum HDL and TC concentrations in postmenopausal women with osteoporosis than in normal persons . The results from our study indicated that Ca supplementation alleviated blood lipid disorders in OVX-LC rats.
We are the first to show that Ca supplementation partially restored the estradiol level in OVX-LC rats. According to the results of previous studies, we speculate that the effect of Ca on osteoporosis may be partially mediated by an increase in E2 levels, which then induces changes in metabolites, especially changes in glycerophospholipid and serum lipid levels, and alters the contents of BMTs, eventually increasing the BMD to a relatively higher level to reduce the degree of osteoporosis.
There are two limitations of the study. First, as mentioned in the introduction, the 2010–2012 Chinese National Nutrition and Health Survey showed that the daily Ca intake of Chinese residents was less than half of the recommended intake, leading to Ca deficiency throughout life and not only after menopause. The OVX rats used in the study only mimics postmenopausal Ca and estradiol deficiency. It does not represent the premenopausal Ca deficiency. Second, the rats used in this study are equivalent to teenagers in humans. Old age is not just about menopause but a number of other factors that change during aging. The Ca supplementation proposed in the study may be effective in teenage equivalent OVX rats, however, until proven, it may not produce the same results when tested in old age rats.
Our results clearly confirm that Ca supplementation is beneficial to reduce bone loss in OVX-LC rats. The main explanation may be that Ca supplementation partially restores the level of estradiol, alters lipid metabolism, induces the production of higher levels of DG (44:6n3), LysoPC (22:2) and PE (P-34:3), and lower levels of Cer (d43:0) and PE-NMe2 (46:3), eventually decreasing the loss of trabecular bone and increasing the BMD of the femur. Originally, researchers proposed that Ca is directly deposited in bone to increase the BMD. A new mechanism by which Ca supplementation potentially increases BMD is by increasing estradiol levels and altering lipid metabolism was proposed. In the future, more studies will be performed to study the effect of Ca on estradiol, and multiomics technology will be used to analyze the upstream and downstream relationship to provide a deeper understanding of the relationship between Ca and osteoporosis.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author upon reasonable request.
Analysis of variance
Bone mineral density
Bone turnover markers
Dual-energy X-ray absorptiometry system
Nonesterified fatty acid
N-Telopeptide of type I collagen
Ovariectomized rats fed a low calcium diet
Ovariectomized rats fed a low calcium diet and treated with 2800 mg/kg CaCO3
Ovariectomized rats fed a low calcium diet and treated with 750 mg/kg CaCO3
Procollagen I N-terminal peptide
Partial least squares discriminant analysis
Sham-operated rats fed a low calcium diet
- Tb.Ar %:
Trabecular area percentage
Ultra-performance liquid chromatography–tandem mass spectrometry
- UPLC–QTOF MS:
Ultra-performance liquid chromatography–quadrupole time-of-flight mass spectrometry
Variable importance in projection
Rachner TD, Khosla S, Hofbauer LC. Osteoporosis: now and the future. Lancet. 2011;377(9773):1276–87.
Solmi M, Veronese N, Correll CU, Favaro A, Santonastaso P, Caregaro L, Vancampfort D, Luchini C, De Hert M, Stubbs B. Bone mineral density, osteoporosis, and fractures among people with eating disorders: a systematic review and meta-analysis. Acta Psychiatr Scand. 2016;133(5):341–51.
Kanis JA. Osteoporosis III: diagnosis of osteoporosis and assessment of fracture risk. Lancet. 2002;359(9321):1929–36.
Cummings SR, Melton LJ. Epidemiology and outcomes of osteoporotic fractures. Lancet. 2002;359(9319):1761–7.
Osteoporosis and bone mineral Diseases chapter of Chinese Medical Association. Results release of Epidemiological survey of osteoporosis in China and “healthy bones” special action. Chinese J Osteoporos Bone Miner Res. 2019;12(4):317–8.
Lips P. Vitamin D deficiency and secondary hyperparathyroidism in the elderly: consequences for bone loss and fractures and therapeutic implications. Endocr Rev. 2001;22(2):477–501.
Ralston SH. Genetic control of susceptibility to osteoporosis. J Clin Endocrinol Metab. 2002;87(6):2460–6.
Stewart TL, Ralston SH. Role of genetic factors in the pathogenesis of osteoporosis. J Endocrinol. 2000;166(2):235–45.
Stathopoulou MG, Grigoriou E, Dedoussis GVZ. Calcium and vitamin D intake interactions with genetic variants on bone phenotype. Curr Nutr Rep. 2012;1:169–74.
Prentice A. Diet, nutrition and the prevention of osteoporosis. Public Health Nutr. 2004;7:227–43.
Bhattarai HK, Shrestha S, Rokka K, Shakya R. Vitamin D, calcium, parathyroid hormone, and sex steroids in bone health and effects of aging. J Osteoporos. 2020;2020:9324505.
Chang J, Wang Y. Report on Chinese Residents’ nutrition and health status survey (2010–2013). Beijing: Peking University Medical Press; 2016. p. 38.
Cashman KD. Calcium intake, calcium bioavailability and bone health. Br J Nutr. 2002;87(Suppl 2):S169–77.
Dawson-Hughes B, Harris SS, Krall EA, Dallal GE. Effect of calcium and vitamin D supplementation on bone density in men and women 65 years of age or older. N Engl J Med. 1997;337(10):670–6.
Power ML, Heaney RP, Kalkwarf HJ, Pitkin RM, Repke JT, Tsang RC, Schulkin J. The role of calcium in health and disease. Am J Obstet Gynecol. 1999;181(6):1560–9.
Nordin BE. The effect of calcium supplementation on bone loss in 32 controlled trials in postmenopausal women. Osteoporos Int. 2009;20(12):2135–43.
Weaver CM, Alexander DD, Boushey CJ, Dawson-Hughes B, Lappe JM, LeBoff MS, Liu S, Looker AC, Wallace TC, Wang DD. Calcium plus vitamin D supplementation and risk of fractures: an updated meta-analysis from the National Osteoporosis Foundation. Osteoporos Int. 2016;27(1):367–76.
Wu J, Xu L, Lv Y, Dong L, Zheng Q, Li L. Quantitative analysis of efficacy and associated factors of calcium intake on bone mineral density in postmenopausal women. Osteoporos Int. 2017;28(6):2003–10.
Thulkar J, Singh S, Sharma S, Thulkar T. Preventable risk factors for osteoporosis in postmenopausal women: systematic review and meta-analysis. J Midlife Health. 2016;7(3):108–13.
Reid IR, Mason B, Horne A, Ames R, Reid HE, Bava U, Bolland MJ, Gamble GD. Randomized controlled trial of calcium in healthy older women. Am J Med. 2006;119(9):777–85.
Barrionuevo P, Kapoor E, Asi N, Alahdab F, Mohammed K, Benkhadra K, Almasri J, Farah W, Sarigianni M, Muthusamy K, Al Nofal A, Haydour Q, Wang Z, Murad MH. Efficacy of pharmacological therapies for the prevention of fractures in postmenopausal women: a network meta-analysis. J Clin Endocrinol Metab. 2019;104(5):1623–30.
Kaddurah-Daouk R, Kristal BS, Weinshilboum RM. Metabolomics: a global biochemical approach to drug response and disease. Annu Rev Pharmacol Toxicol. 2008;48:653–83.
Newgard CB. Metabolomics and metabolic diseases: where do we stand? Cell Metab. 2017;25(1):43–56.
Johnson CH, Ivanisevic J, Siuzdak G. Metabolomics: beyond biomarkers and towards mechanisms. Nat Rev Mol Cell Biol. 2016;17(7):451–9.
Lv H, Jiang F, Guan D, Lu C, Guo B, Chan C, Peng S, Liu B, Guo W, Zhu H, Xu X, Lu A, Zhang G. Metabolomics and its application in the development of discovering biomarkers for osteoporosis research. Int J Mol Sci. 2016;17(12):2018.
Si Z, Zhou S, Shen Z, Luan F. High-throughput metabolomics discovers metabolic biomarkers and pathways to evaluating the efficacy and exploring potential mechanisms of osthole against osteoporosis based on UPLC/Q-TOF-MS coupled with multivariate data analysis. Front Pharmacol. 2020;11:741.
Zhao H, Li X, Zhang D, Chen H, Chao Y, Wu K, Dong X, Su J. Integrative bone metabolomics–lipidomics strategy for pathological mechanism of postmenopausal osteoporosis mouse model. Sci Rep. 2018;8(1):16456.
Cui L, Wu T, Liu YY, Deng YF, Ai CM, Chen HQ. Tanshinone prevents cancellous bone loss induced by ovariectomy in rats. Acta Pharmacol Sin. 2004;25(5):678–84 (PMID: 15132837).
Liu XY, Zhang AH, Fang H, Li MX, Song Q, Su J, Yu MD, Yang L, Wang XJ. Serum metabolomics strategy for understanding the therapeutic effects of Yin-Chen-Hao-Tang against Yanghuang syndrome. RSC Adv R Soc Chem. 2018;8(14):7403–13.
Chu H, Zhang A, Han Y, Lu S, Kong L, Han J, Liu Z, Sun H, Wang X. Metabolomics approach to explore the effects of Kai-Xin-San on Alzheimer’s disease using UPLC/ESI-Q-TOF mass spectrometry. J Chromatogr B Analyt Technol Biomed Life Sci. 2016;1015–1016:50–61.
Chong J, Soufan O, Li C, Caraus I, Li S, Bourque G, Wishart DS, Xia J. MetaboAnalyst 4.0: towards more transparent and integrative metabolomics analysis. Nucleic Acids Res. 2018;46(W1):W486–94.
Zhu K, Prince RL. Calcium and bone. Clin Biochem. 2012;45(12):936–42.
Karpouzos A, Diamantis E, Farmaki P, Savvanis S, Troupis T. Nutritional aspects of bone health and fracture healing. J Osteoporos. 2017;2017:4218472.
Yamazaki I, Yamaguchi H. Characteristics of an ovariectomized osteopenic rat model. J Bone Miner Res. 1989;4(1):13–22.
Omi N, Ezawa I. The effect of ovariectomy on bone metabolism in rats. Bone. 1995;17(4 Suppl):163S-168S.
Kanis JA, Melton LJ 3rd, Christiansen C, Johnston CC, Khaltaev N. The diagnosis of osteoporosis. J Bone Miner Res. 1994;9(8):1137–41.
Rico H, Valencia MJ, Villa LF, Hernández ER, Seco C, Sánchez-Atrio A, Revilla M. Calcitonin versus clodronate in the prevention of ovariectomy-induced osteopenia in rats. Clin Rheumatol. 2000;19(1):47–50.
Gao X, Ma W, Dong H, Yong Z, Su R. Establishing a rapid animal model of osteoporosis with ovariectomy plus low calcium diet in rats. Int J Clin Exp Pathol. 2014;7(8):5123–8.
Seibel MJ. Biochemical markers of bone turnover: part I: biochemistry and variability. Clin Biochem Rev. 2005;26(4):97–122.
Cox G, Einhorn TA, Tzioupis C, Giannoudis PV. Bone-turnover markers in fracture healing. J Bone Jt Surg Br. 2010;92(3):329–34.
Devogelaer JP, Boutsen Y, Gruson D, Manicourt D. Is there a place for bone turnover markers in the assessment of osteoporosis and its treatment? Rheum Dis Clin North Am. 2011;37(3):365–86.
Bauer DC, Garnero P, Harrison SL, Cauley JA, Eastell R, Ensrud KE, Orwoll E, Osteoporotic Fractures in Men (MrOS) Research Group. Biochemical markers of bone turnover, hip bone loss, and fracture in older men: the MrOS study. J Bone Miner Res. 2009;24(12):2032–8.
Dai Z, Wang R, Ang LW, Yuan JM, Koh WP. Bone turnover biomarkers and risk of osteoporotic hip fracture in an Asian population. Bone. 2016;83:171–7.
Eastell R, Szulc P. Use of bone turnover markers in postmenopausal osteoporosis. Lancet Diabetes Endocrinol. 2017;5(11):908–23.
Weisman SM, Matkovic V. Potential use of biochemical markers of bone turnover for assessing the effect of calcium supplementation and predicting fracture risk. Clin Ther. 2005;27(3):299–308.
Garnero P, Sornay-Rendu E, Claustrat B, Delmas PD. Biochemical markers of bone turnover, endogenous hormones and the risk of fractures in postmenopausal women: the OFELY study. J Bone Miner Res. 2000;15(8):1526–36.
Cui G, Leng H, Wang K, Wang J, Zhu S, Jia J, Chen X, Zhang W, Qin L, Bai W. Effects of remifemin treatment on bone integrity and remodeling in rats with ovariectomy-induced osteoporosis. PLoS ONE. 2013;8(12):e82815.
Han J, Wang W. Effects of tanshinol on markers of bone turnover in ovariectomized rats and osteoblast cultures. PLoS ONE. 2017;12(7):e0181175.
Zhang Y, Lai WP, Leung PC, Wu CF, Wong MS. Short- to mid-term effects of ovariectomy on bone turnover, bone mass and bone strength in rats. Biol Pharm Bull. 2007;30(5):898–903.
Patlas N, Zadik Y, Yaffe P, Schwartz Z, Ornoy A. Oophorectomy-induced osteopenia in rats in relation to age and time postoophorectomy. Cells Tissues Organs. 2000;166(3):267–74.
Lee MY, Kim HY, Singh D, Yeo SH, Baek SY, Park YK, Lee CH. Metabolite profiling reveals the effect of dietary Rubus coreanus vinegar on ovariectomy-induced osteoporosis in a rat model. Molecules. 2016;21(2):149.
Liu X, Liu Y, Cheng M, Zhang X, Xiao H. A metabolomics study of the inhibitory effect of 17-beta-estradiol on osteoclast proliferation and differentiation. Mol Biosyst. 2015;11(2):635–46. https://doi.org/10.1039/c4mb00528g.
Ebrahimpour P, Fakhrzadeh H, Heshmat R, Ghodsi M, Bandarian F, Larijani B. Metabolic syndrome and menopause: a population-based study. Diabetes Metab Syndr Clin Res Rev. 2010;4(1):5–9.
Levin VA, Jiang X, Kagan R. Estrogen therapy for osteoporosis in the modern era. Osteoporos Int. 2018;29(5):1049–55.
Rajfer RA, Flores M, Abraham A, Garcia E, Hinojosa N, Desai M, Artaza JN, Ferrini MG. Prevention of osteoporosis in the ovariectomized rat by oral administration of a nutraceutical combination that stimulates nitric oxide production. J Osteoporos Hindawi. 2019;2019:1–11.
Geng Q, Gao H, Yang R, Guo K, Miao D. Pyrroloquinoline quinone prevents estrogen deficiency-induced osteoporosis by inhibiting oxidative stress and osteocyte senescence. Int J Biol Sci. 2019;15:58–68.
Wu C, Chen J, Cheng J, Liang L, Zhu S. Effect of Shoutai Pill on estrogen content and bone morphometry in ovariectomized rats. Trad Chinese Med Clin Pharmacol. 2014;25(3):295–8.
Shen J, Shen M, Li Z. Effect of electroacupuncture on estrogen and endometrium in ovariectomized rats. Liaoning J Tradit Chinese Med. 2009;36(12):2167–8.
Terasawa E, Kenealy BP. Neuroestrogen, rapid action of estradiol, and GnRH neurons. Front Neuroendocrinol. 2012;33(4):364–75.
Nelson LR, Bulun SE. Estrogen production and action. J Am Acad Dermatol. 2001;45(3 Suppl):S116–24.
Oakley OR, Kim KJ, Lin PC, Barakat R, Cacioppo JA, Li Z, Whitaker A, Chung KC, Mei W, Ko C. Estradiol synthesis in gut-associated lymphoid tissue: leukocyte regulation by a sexually monomorphic system. Endocrinology. 2016;157(12):4579–87.
Eyre LJ, Bland R, Bujalska IJ, Sheppard MC, Stewart PM, Hewison M. Characterization of aromatase and 17 beta-hydroxysteroid dehydrogenase expression in rat osteoblastic cells. J Bone Miner Res. 1998;13(6):996–1004.
Chen YY, Wang WW, Yang L, Chen WW, Zhang HX. Association between lipid profiles and osteoporosis in postmenopausal women: a meta-analysis. Eur Rev Med Pharmacol Sci. 2018;22(1):1–9.
This study was supported by grants from the Center for Food Evaluation, State Food and Drug Administration (Grant No. ZBW-2017-BJSP-01).
This animal study was performed according to institutional guidelines with the approval of the Animal Ethics Committee of National Institute of Nutrition and Health, Chinese Center for Diseases Control and Prevention.
Consent for publication
All authors have declared that no conflicts of interest exist.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Mao, H., Wang, W., Shi, L. et al. Metabolomics and physiological analysis of the effect of calcium supplements on reducing bone loss in ovariectomized rats by increasing estradiol levels. Nutr Metab (Lond) 18, 76 (2021). https://doi.org/10.1186/s12986-021-00602-y
- Bone mineral density