Oral Administration of Quercetin or Its Derivatives Inhibit Bone Loss in Animal Model of Osteoporosis

Objectives. Quercetin (Q) and its derivatives are the major members of the naturally occurring flavonoid family, which possess beneficial effects on disease prevention including osteoporosis. The present study is aimed at further investigating the efficacy of the Q and its derivatives on bone pathology, bone-related parameters under imageology, bone maximum load, and serum bone metabolism indexes in animal model of osteoporosis. Potential mechanisms of Q and its derivatives in the treatment of osteoporosis as well as the existing problems regarding the modeling method and limitations of researches in this area were also summarized. Eight databases were searched from their inception dates to February 2020. Nineteen eligible studies containing 21 comparisons were identified ultimately. The risk of bias and data on outcome measures were analyzed by the CAMARADES 10-item checklist and Rev-Man 5.3 software separately. The results displayed the number of criteria met varied from 3/10 to 7/10 with an average of 5.05. The present study provided the preliminary preclinical evidence that oral administration of Q or its derivatives was capable of improving bone pathology, bone-related parameters under imageology and bone maximum load, increasing serum osteocalcin, alkaline phosphatase, and estradiol, and reducing serum c-terminal cross-linked telopeptide of type I collagen (P < 0.05). No statistical difference was seen in survival rate, index of liver, or kidney function (P > 0.05). Q and its derivatives partially reverse osteopenia probably via antioxidant, anti-inflammatory, promoting osteogenesis, inhibiting osteoclasts, and its estrogen-like effect. The findings reveal the possibility of developing Q or its derivatives as a drug or an ingredient in diet for clinical treatment of osteoporosis.


Introduction
Osteoporosis, as a systemic skeletal disease characterized by loss of bone mass and bone microarchitectural deterioration, causes enhanced bone fragility and a greater risk of fractures [1]. Some risk factors are bound up with osteopathic fracture, containing age, race, cigarette smoking, low physical activity, hormone-related factors, the use of drugs (e.g., glucocorticoids), low calcium and vitamin D levels, and prior history of fracture [2]. It gradually becomes a major public health issue affecting exceeded 200 million people worldwide [3], and the burden of society is continuously increasing with the aging of the world population [4]. Thus, how to effec-tively prevent and treat osteoporosis has attracted more and more attention of the world.
Physical activity [5], supplementing with calcium and vitamin D, [6,7] and pharmacological therapies [8] (including estrogen [9], bisphosphonates [10], denosumab [11], or teriparatide [12]) are the most commonly used clinical approaches for the prevention or treatment of osteoporosis. However, there is no definite solution for this disease because there are still some problems in the above treatment methods [13]. For example, although bisphosphonates are widely used globally as the main treatment for osteoporosis, it does not include proven efficacy against nonvertebral fractures, and the evidence for benefit beyond 5 years in intervention studies is limited [14]. Besides, there are concerns that atypical femur fractures and osteonecrosis of the jaw may be caused by long-term use of bisphosphonates [14]. Highimpact exercise (e.g., running) is conducive to producing changes that reflect distal tibial bone adaptation [5]. However, a well-designed training intervention is difficult to achieve in normal people much less in patients with osteoporosis or fracture. The benefits of supplementing with vitamin D and calcium for fracture prevention and holistic fall remain uncertain [15]. There are inconsistent findings for adverse effects on the cardiovascular system of calcium supplements with or without vitamin D [7,16,17], and the adverse effects on the skeleton and musculature of highdose vitamin D were reported by recent evidence [18,19]. Despite its excellent antiosteoporosis effect for postmenopausal osteoporosis, long-term use of estrogen for osteoporosis caused many side effects, containing an increased incidence of cardiovascular accidents, endometrial cancer [20][21][22], as well as the loss of mesenteric ganglion neurons and dominant ovarian follicles [23]. Therefore, obtaining a new effective drug without side effects is needed urgently in the treatment of osteoporosis.
Quercetin (Q, C 15 H 10 O 7 , Figure 1(a)) and its derivatives (Figures 1(b) and 1(c)) are the major members of the naturally occurring flavonoid family. Flavonoid is well known for its antioxidant and anti-inflammatory properties. Q and its derivatives have been widely used as a natural antioxidant in the treatment of cardiovascular diseases, tumor, and some metabolic diseases including osteoporosis for many years in China [24,25]. Recent studies [26,27] demonstrate that Q and its derivatives are potent natural osteogenic agents by multiple pharmacological activities including the function of antioxidant, anti-inflammatory, and estrogen-like effect in vivo and vitro studies. However, the scattered evidence and uncertain mechanisms limited the application of Q and its derivatives in the clinic. Systematic review and metaanalysis of animal studies are considered to be a valuable tool to provide important insights into the validity of animal studies, improve the precision of estimated effects, and support further generalization to human clinical trials [28]. Thus, the present study is aimed at investigating the preclinical evidence and possible mechanisms of Q and its derivatives in animal models of osteoporosis.

Data Sources and Search Strategies.
A comprehensive literature search about animal experimental studies of Q or its derivatives for osteoporosis was conducted in the Chinese Science and Technology Journal Database, WanFang, China National Knowledge Infrastructure, Chinese Biomedical Database, PubMed, EMBASE, Cochrane library, and Web of Science database from their inception dates to February 2020. The following search terms were used in PubMed and were modified to suit other databases: "Quercetin" AND "Osteoporosis OR Bone Loss". In addition, reference lists from the resulting publications and reviews were also searched carefully for the eligible studies. Wang) selected the studies separately by browsing the abstracts and full texts via the eligibility criteria. The study was included if it met the following criteria: (1) controlled studies assessing the administration of Q and its derivatives for osteoporosis or bone loss animal models established by various ways; (2) the treatment group received Q or its derivatives as monotherapy with unrestricted dosage, medicament type, route of administration, and time for the medicine application. Blank treatment or isasteric placebo was received in the control group; (3) the present study received bone pathology and/or bone mineral density (including femur bone mineral density (F-BMD), lumbar spine bone mineral density (L-BMD)) and/or bone histomorphometric parameters under micro-CT (trabecular number (Tb.N), trabecular thickness (Tb.Th)) and/or bone maximum load and/or bone turnover markers (serum alkaline phosphatase (ALP), C-terminal cross-linked telopeptide of type I collagen (CTX) and osteocalcin (OC)) and/or serum estradiol and/or uterine weight and/or indicators of adverse reactions as the primary outcome measures, while the antiosteoporosis mechanisms of Q or its derivatives was selected as the second outcome measures. Exclusion criteria was as below: (1) other types of studies (in vitro studies, case reports, clinical trials, reviews, abstracts, comments, and editorials);

Eligibility Criteria and Data Extraction. Two authors (Yue-Yue Huang and Zi-Hao
(2) combination with other compounds; (3) compared with other traditional Chinese medicine; (4) no any primary outcome indicator were involved or incomplete date; (5) inconsistent of graphic and textual data; (6) no control group; (7) duplicate publications; (8) not osteoporosis or bone loss model.
The details were extracted from included studies by two independent authors (Yue-Yue Huang and Zi-Hao Wang) using a predefined form. The information included the authors and years of publication; information of animals; modeling method; the use of anesthetics anaesthetize in the course of the experiment; the therapeutic regimen of treatment and control group; and primary and/or secondary outcomes and its intergroup differences. Only data from the osteoporosis group and Q+osteoporosis groups were included for analysis when a study is involved in multiple intervention groups. When the outcomes were displayed through gradient doses of drug therapy or determined at different times, only the data of the highest dose and the final measurement was included for analysis.

Risk of Bias in Individual
Studies. CAMARADES 10-item quality checklist [29] with minor modification was used to assess the study quality by two independent authors (Yue-Yue Huang and Li-Hui Deng). The modification is listed as follows: D: blinded induction of model (group randomly after modeling); F: use of anesthetic without significant bony protective activity or nephrotoxicity; G: appropriate animal model with complications or risk factors (including age, hyperlipemia, diabetes, or hypertensive). Disagreements in the process of selecting studies, extracting data, and assessing the quality of studies were resolved by consensus or arbitration by the correspondence author (Qun Zheng).

2
Oxidative Medicine and Cellular Longevity 2.4. Statistical Analysis. The RevMan 5.3 software was used for data analysis where possible; otherwise, comparison between groups was performed. The bar graphs were drawn via Prism 6. In meta-analysis, standardized mean differences (SMDs) and 95% confidence intervals (95% CIs) were calculated to estimate the combined overall effect sizes when the outcomes were determined in various ways or the unit of measurement is different. Heterogeneity was assessed using the Cochrane Q-statistic test (P < 0:05 was considered statistically significant) and the I 2 -statistic test (I 2 < 50% was considered homogeneous). Random (I 2 > 50%) or fixed-effects model (I 2 < 50%) was selected according to the results of I 2 .
In order to ensure the reliability of results, a sensitivity analysis was performed, and potential publication bias was assessed by the visual inspection of the funnel plot and asymmetry test. Moreover, in order to explore the impact of potential confounding factors on the estimates of combined effect size, subgroup analyses were conducted in this study. The significance level was set at P < 0:05.

Study Quality.
Detailed results of methodological quality are presented in Table 3. The number of criteria met in each study varied from 3/10 to 7/10 with the average of 5.05. Only 2 studies [30,39] were not a peer-reviewed publication, and 3 studies [32,44,45] did not mention control of temperature. Two [31,46] of the included studies did not declare randomization. The ways of blinding induction of model were reported in 5 studies [4,31,35,36,43], and all of them reported the animals were grouped randomly after modeling. No study mentioned the calculation of sample size and none used a blinding method during outcome assessment and appropriate animal model.  [31,38,41] utilized bone pathology as a primary outcome measure. Among them, 1 study [31] reported that osteoporotic rats treated with Q showed marked the improvement of the structure of femoral cortical bone compared with osteoporotic rats induced by 13cRA, which showed thickness was nearly similar to that of the control group although a few small intracortical cavities were still present. One study [38] reported that Q treatment was observed to prevent trabecular fracture and osteoblast apoptosis and maintain normal distribution of trabecular. Another study [41] reported that the administration of 150 mg/kg Q increased femoral trabecular and cortical thickness by 36% and 22%, respectively, compared with the osteoporotic rats induced by methylprednisolone sodium succinate.         Figure 6). In consideration of high heterogeneity, sensitivity analyses of the above indicators were carried out, and the result showed that the heterogeneity did not substantially alter after removing any 1 study.

Serum ALP, OC, and CTX.
Compared with the control group, Q and its derivatives were reported that they existed Note: studies fulfilling the criteria of: A: peer-reviewed publication; B: control of temperature; C: random allocation to treatment or control; D: blinded induction of model (group randomly after modeling); E: blinded assessment of outcome; F: use of anesthetic without significant protective and toxic effects on bones; G: appropriate animal model (aged, hyperlipemia, hypertensive, or diabetes); H: sample size calculation; I: compliance with animal welfare regulations (including three or more of the following points: preoperative anaesthesia, postoperative analgesia, nutrition, disinfection, environment temperature, environment humidity, circadian rhythm, and euthanasia); J: statement of potential conflict of interests.  11 Oxidative Medicine and Cellular Longevity significant effect for increasing serum ALP (P < 0:05) in 7 studies [30,[33][34][35][36]42], no significant effect for serum ALP (P > 0:05) in 1 study [31], and significant effect for reducing serum ALP (P < 0:05) in 1 study [45]. Besides, 5 studies [30,35,37,41,42] and 3 studies [34,35,41] reported separately that Q could increase the serum OC (P < 0:05) and/or reduce the serum CTX (P < 0:05).

Indicators of Adverse
Reactions. Four studies [4,33,38,45] utilized survival rate as primary outcome measure, and meta-analysis of 3 studies showed no statistical difference of Q or its derivatives on survival rate (n = 94, RR 1.00, 95% CI (0.89, 1.13), P = 0:94; heterogeneity: χ 2 = 2:34, I 2 = 0%, Figure 9). BUN was measured in 2 studies [31,34] and SCr in 1 study [34] to assess the adverse effect of Q to the renal function. In addition, AST and ALT were measured in 1 study [31] to assess the adverse effect of Q to the liver, and   (b) Figure 5: (a) The forest plot: effects of Q or its derivatives for increasing Tb.Th compared with the control group; (b) The forest plot: effects of Q or its derivatives for increasing Tb.N compared with the control group.

12
Oxidative Medicine and Cellular Longevity the results showed that there was no statistical difference in renal and liver function between Q group and control group (P > 0:05).

Subgroup
Analysis. F-BMD was reported to be improved greatly in 15 comparisons [4, 27, 30-34, 39, 42-46]. The potential confounding factors which may increase the heterogeneity of F-BMD were explored via subgroup analysis. First, we divide the 15 comparisons into the ovariectomized model group and nonovariectomized model group according to different modeling methods. As the result, no significant difference was observed in the effect size of two groups (SMD = 2:00 ± 0:36 versus SMD = 1:93 ± 0:58, P = 0:85, Figure 10(a)) and heterogeneity of both groups did not decrease obviously. On the basis of the result of the subgroup analysis above, we analyzed the F-BMD in different subgroups stratified according to the following variables in the ovariectomized model group: different animal species, different laboratory drugs, and the duration of treatment. In the subgroup analysis of these factors, the mice group showed better effect size than the rat group (SMD = 5:55 ± 2:64 versus SMD = 1:93 ± 0:37, P = 0:008, Figure 10(b)) with significantly reduced heterogeneity of both groups. Another subgroup analysis indicated that the effect of Q derivatives was better than Q with better evaluation effect size (SMD = 3:28 ± 0:78 versus SMD = 1:64 ± 0:41, P = 0:0003, Figure 10(c)), and the heterogeneity experienced a marked decline in Q derivatives treatment group. In addition, the shorter period of Q or its derivatives treatment showed batter effect size than the longer treatment (SMD = 5:55 ± 2:64 versus SMD = 1:69 ± 0:43 versus SMD = 2:16 ± 0:64, P = 0:01, Figure 10(d)).

Discussion
4.1. Summary of Evidence. The first-ever preclinical systematic review included a batch of studies of acceptable quality to estimate the efficacy and mechanisms of Q and its derivatives in animal models of osteoporosis. The findings revealed Q and its derivatives are potential antiosteoporosis drug via multiple mechanisms.

4.2.
Limitations. Some limitations of the meta-analysis and the system evaluation were listed as follows: (1) there may still be a certain degree of selective bias due to the lack of negative studies and the studies from other databases or in other languages; (2) defects in aspects of blinding assessment of outcome and sample size calculation may affect the accuracy of findings [47]; (3)   13 Oxidative Medicine and Cellular Longevity between groups are carried out rather than systematic evaluations in those studies; (4) no study utilized animals with relevant complication.

Implication
4.3.1. Animal Model Selection. Using different animal models at different research stages of disease is crucial to study it pathophysiology and treatments [48]. Factors that need to be considered include pathogenesis of model, availability of the animals, technical requirements, and cost and ethical considerations [49]. According to the pathogenesis, animal models of osteoporosis can be divided into two types: models with increased bone resorption as the dominant mechanism (such as ovariectomized osteoporosis model, disused osteoporosis model, retinoic acid induction model, nutritional osteoporosis model, and glucocorticoid model) and models with reduced bone formation as the dominant mechanism (such as senile osteoporosis model and n-ZnO induction model) [50]. The present study comprehensively contains the ovariectomized osteoporosis model, retinoic acid or n-ZnO induction model, diabetic osteoporosis, and glucocorticoid model to estimate the efficacy and mechanisms of Q and its derivatives for osteoporosis. And the results of metaanalysis and subgroup analysis reflected that Q and its deriv-atives could play a role in both two mechanisms of osteoporosis (SMD = 2:00 ± 0:36 versus SMD = 1:93 ± 0:58, P = 0:85, Figure 10(a)). However, some key points to establish incorporating models still deserve attention. The ovariectomized rats, as the best recognized postmenopausal osteoporosis model, are the most adopted model in eligible studies. Animals with mature skeleton are obligatory for osteoporosis researches [51]. In the present study, some immature rats (<12 weeks) whose bone mass was below its peak [32,33,43] were used to establish models, causing that the confounding factors are introduced to animals that are still accruing bone [51]. Thus, the application of animals with mature bone needs to be emphasized for future osteoporosis experiments. In addition, the dosage of glucocorticoids should be grasped well to avoid no change in bone mass in low dosage or animal death due to overimmunization in high dosage. Preexperiment as Derakhshanian et al. [41] did to explore the appropriate dosage is a recommended practice.
The subgroup analysis showed better effect size in the mice group than the rat group (SMD = 5:55 ± 2:64 versus SMD = 1:93 ± 0:37, P = 0:008, Figure 10(b)) with significantly reduced heterogeneity of both groups, suggesting that different animals may be one of the main sources of heterogeneity. The commonly used animals in osteoporosis experiments are rodent, rabbit, dog, sheep, primates, and so on     14 Oxidative Medicine and Cellular Longevity [52,53]. Rodents such as rats and mice possess the advantages of being repeatable, cheap, and convenient to be bred and anesthetized. It also takes less time to form a new balance of bone remodeling. After ovariectomy in rodents, the bone mass of cancellous bone decrease and the bone turnover rate increase, which resembles to that of osteoporosis in postmenopausal women and estrogen replacement therapy could alleviate bone mass loss [54]. Therefore, female rodents are widely used in the study of postmenopausal osteoporosis. It was regarded as a preferred animal for small-or mediumsized laboratories to investigate the efficacy and mechanisms of drugs. However, it is not suitable for the study of the bone cortex because of the absence of haversian system in the bone cortex of rats or mice. Additionally, difficulties arose when the arm of the study was to implant fixation or prosthetic devices as well as that studies which need several the collection of high blood volumes or surgical operation or several biopsies due to its small size especially in mice [55]. Compared to rodents, rabbits and dogs have the convenience of appropriate cost in terms of purchase and maintenance and reasonable anatomical size which was advantageous to biopsies, blood collection repeatedly, and surgical treatments such as bone implant and bone-implant interface. On the other hand, it is suitable for the study of the effect of cortical bone because of the obvious Harvard reconstruction. However, no significant change was observed in bone mineral density after ovariectomy in rabbits and dogs [56]. The combination of surgery and subsequent glucocorticoid treatment is the best solution at present to obtain bone mineral reduction consistently in the short term in rabbits [48]. However, it cannot be achieved in dogs now which may be related to the low estrogen level in female dogs. Apart from the advantages of rabbits and dogs, sheep are ideal models for the study of vertebral osteoporosis [57] and they also have similar toxic effects on osteocytes by fluoride to humans [58]. Unfortunately, they are not suitable to be adopted to study the efficacy of oral administration of drugs given that the sheep are ruminants. Additionally, bone mineral density, blood biochemical parameters, and bone histomorphology of sheep all fluctuated seasonally [59] which may influence the accuracy of experimental results. From both a physiological and anatomical standpoint, the characteristics of the skeleton of primates are most close to humans than any other type of animal. Nevertheless, obtaining licenses to use them as experimental animal become increasingly difficult due to ethical considerations [60] and potential epidemic animal-borne diseases in primates [61]. Meanwhile, the high cost to purchase and maintain primates restricts their use in experimentation [49]. In consideration of that, the arm of the present studies is to explore the preliminary effect of and mechanisms of Q in animal models of osteoporosis, thus the use of rodents is acceptable at this stage. We suggest the advanced animals (sheep, primates, etc.) or animal models (transgenic or knockout rodents, etc.) for osteoporosis should be chosen in the future on the basis of the experimental purpose and permissible conditions.

Other Subgroup Analyses.
A better evaluation effect size in a group of Q derivatives (SMD = 3:28 ± 0:78 versus SMD = 1:64 ± 0:41, P = 0:0003, Figure 10(c)) was indicated in   17 Oxidative Medicine and Cellular Longevity subgroup analysis based on contrast of the overall effect in different experiments. Given the existence of mixed factors in different experiments, all included studies have been perused again for these studies were designed to contrast the different therapeutic effects between Q and its derivatives under the same experimental conditions. Among, QCG (quercetin-6-C-A-D-glucopyranoside) [27] and rutin (quercetin-3-O-B-rutinoside) [33] were reported to improve bone biomechanical quality more effectively than Q via positive modifications of bone microarchitecture and bone mineral density without hyperplastic effect on the uterus, which possibly is related to that the synthesis of modified groups attached to Q improves solubility and bioavailability [62]. In this context, all the existing derivatives of Q were collected and listed in Table 4, which are recommended as potential antiosteoporosis drugs in the future researches. Another subgroup analysis showed the shorter period of Q or its derivative treatment showed better effect size than the longer treatment (SMD = 5:55 ± 2:64 versus SMD = 1:69 ± 0:43 versus SMD = 2:16 ± 0:64, P = 0:01, Figure 10(d)), suggesting that the duration of treatment may be a source of high heterogeneity. For the reason, we attribute it to that osteoporosis is a progressive and irreversible disease when pathogeny persist, extending the treatment time of Q and its derivatives is merely conducive to delay the progression of osteoporosis rather than reverse it.

Possible Mechanisms.
Systemic review of preclinical studies is conducive to understand comprehensively pathological mechanisms of disease and pharmacological effects of drugs [63]. We have summarized the possible mechanisms of Q and its derivatives mediated bone protection from cur-rent findings and listed them as follows: (1) Q alleviated oxidative damage by decreasing NO [36] and increasing GSH [31,40], SOD [39,45], and CAT [45] to reduce the release of MDA [31,39] in the ovaries and bone tissue. A decreased MDA level in the ovaries directly increased estrogen activity [31] which has been shown to have antioxidant properties [64,65]. In addition, Gsk-3β/Nrf2 signal pathway was reported to participate in the regulation of the abovementioned antioxidant process [39]. (2) Q and quercetin-3-Orutinose have the similar effect of phytoestrogen on inhibiting bone resorption by participating in the binding of estrogen receptor (ER) [33,45,66], especially in ERβ mainly expressed in bone [67]. However, another study [43] showed that Q did not appear to have this activity through either ERα or ERβ, which suggested that Q might affect bone metabolism through ERs independent pathway. (3) Q alleviated inflammatory reaction by inhibiting the expression of TNFα [4,32,36], IL-6 [4,36], INF-γ [4], and CRP [36]. And then TNF-α activated NF-κB, increased the expression of NF-κB protein, and promoted the degradation of β-catenin protein [31]. (4) Q promoted bone synthesis by enhancing the expression of osteogenic protein (FoxO1, Bglap2, Collal, Osterix, and Runx2) [30,39] via PI3K/Akt/Fox O1/NF-κB signal pathway [30] and BMP2/smad4 signaling pathway [35]. (5) QCG inhibited the expression of osteoclast markers including RANK and c-fos in bone marrow cells (BMCs) cultured in the presence of RANK ligand and M-CSF [27]. Q and quercetin-3-O-beta-D-glucuronide inhibited RANKinduced osteoclast formation in a dose-dependent manner in RAW264.7 cells, and the RANK ligand-stimulated expression of osteoclast related genes including NFATc1 was inhibited by Q [43]. In addition, as the most abundant p38  Figure 11: A schematic representation of osteoprotective mechanisms of Q and its derivatives for osteoporosis. 18 Oxidative Medicine and Cellular Longevity member in the bone and bone marrow [68], the ablation of p38 MAPK signaling in osteoblast lineage cells protects the mice from bone loss induced by estrogen deficiency [34,69] found that Q could attenuate osteoporosis by downregulating MAPK signaling pathways. (6) Prostaglandins (PGs) played a role in IL-1-induced bone resorption [70]. Q was observed to reduce the production of PGs by inhibiting cyclooxygenase and phosphoesterase A2 [71]. The mechanism diagram is summarized in Figure 11.

Conclusion
The present study provided the preliminary preclinical evidence that oral administration of Q and its derivatives was capable of partially reversing osteopenia in animal models probably via antioxidant, anti-inflammatory, promoting osteogenesis, inhibiting osteoclasts, and its estrogen-like effect. The findings reveal the possibility of developing Q and its derivatives as a drug or an ingredient in diet for the clinical treatment of osteoporosis.