Efficacy and safety of medical therapy for low bone mineral density in patients with Crohn disease: A systematic review with network meta-analysis.

BACKGROUND: Low bone mineral density (BMD) is a frequent complication of inflammatory bowel disease (IBD), particularly in patients with Crohn disease (CD). The aim of our study is to determine the efficacy and safety of different drugs used to treat low BMD in patients with CD.
METHODS: PUBMED/MEDLINE, EMBASE, and Cochrane Central Register of Controlled Trials were searched for eligible studies. A random-effects model within a Bayesian framework was applied to compare treatment effects as standardized mean difference (SMD) with their corresponding 95% credible interval (CrI), while odds ratio (OR) was applied to compare adverse events with 95% CrI. The surface under the cumulative ranking area (SUCRA) was calculated to make the ranking of the treatments for outcomes.
RESULTS: Twelve randomized controlled trials (RCTs) were eligible. Compared with placebo, zoledronate (SMDs 2.74, 95% CrI 1.36-4.11) and sodium-fluoride (SMDs 1.23, 95% CrI 0.19-2.26) revealed statistical significance in increasing lumbar spine BMD (LSBMD). According to SUCRA ranking, zoledronate (SUCRA = 2.5%) might have the highest probability to be the best treatment for increasing LSBMD in CD patients among all agents, followed by sodium-fluoride (27%). For safety assessment, the incidence of adverse events (AEs) demonstrated no statistical difference between agents and placebo. The corresponding SUCRA values indicated that risedronate (SUCRA = 77%) might be the most safe medicine for low BMD in CD patients and alendronate ranked the worst (SUCRA = 16%).
CONCLUSIONS: Zoledronate might have the highest probability to be the best therapeutic strategy for increasing LSBMD. For the safety assessment, risedronate showed the greatest trend to decrease the risk of AEs. In the future, more RCTs with higher qualities are needed to make head-to-head comparison between 2 or more treatments.

Introduction

Crohn disease (CD), one of major phenotype of inflammatory bowel disease (IBD), is characterized by chronic relapsing inflammatory disorder of the gastrointestinal tract.[ In the past 50 years, the incidence and prevalence of CD have increased up to 6 to 15/100,000 and 50 to 200/100,000 persons, respectively, in the west.[ Various extra-intestinal damages have been reported in CD patients, of these, low bone mineral density (BMD) is considered a common extra-intestinal manifestation.[ Low BMD is often referred to as osteopenia and osteoporosis in the clinic. The prevalence of osteopenia in patients with IBD was up to 62% while osteoporosis occurred in 38%.[ Although some studies showed that ulcerative colitis was also associated with low BMD,[ a majority of studies demonstrated that the lower bone density was more likely to be involved in CD.[ In patients with CD, the risk factors contributing to the development of low BMD mainly include: sex, malabsorption of vitamin D, calcium, and other nutrients, the use of corticosteroids, inflammatory cytokines, and bowel resection.[ Moreover, it has been reported that low BMD was tightly associated with osteoporosis-related fractures,[ thus, it is necessary for CD patients to receive the treatment against osteopenia and osteoporosis.

Some randomized controlled trials (RCTs) have demonstrated the effectiveness and safety of several drugs, such as the antiresorptive and bone-anabolic drugs. These clinical trials would benefit the decision-making strategies in clinical practices.[ Recently, a meta-analysis concluded that both bisphosphonates (e.g., alendronate, risedronate, ibandronate, and zoledronate) and sodium fluoride were effective in improvement of low BMD in patients with IBD.[ However, it was only conducted on the basis of pairwise evidence by comparing various medical interventions with placebo. Thus, the optimal selection of drugs still remains controversial because it was unlikely to acquire comparisons between each available treatment options.

Given the lack of head-to-head RCTs between antiresorptive and bone-anabolic drugs for low BMD in CD patients, we conducted a systematic review with network meta-analysis, which permits the integration of direct and indirect comparisons and allows us to simultaneously compare multiple treatments.[ In order to enhance statistical power and improve estimates of effect size, we analyzed and combined results of previous RCTs to assess the effectiveness and safety of antiresorptive and bone-anabolic drugs for low BMD in CD patients.

Methods

Search strategy

Electronic databases of PUBMED, MEDLINE, EMBASE, and Cochrane Central Register of Controlled Trials were searched from inception to 7 April 2016 for RCTs. Medical subject heading (MeSH) terms and keywords as follows: “inflammatory bowel disease,” “IBD,” “Crohn disease,” “CD,” “bone mineral density,” “bone density,” “osteoporosis,” “osteopenia,” “osteoporotic,” “OP,” and “bone loss” were used in the search strategy without language or date restrictions. We also scanned previous systematic reviews and meta-analysis to identify potential relevant studies. Titles and abstracts were carefully screened. Any disagreement between the two authors was resolved by further discussion. If there was no consensus, a third reviewer (corresponding author) was consulted. All analyses were based on previous published studies, thus ethical approval and patient consent were not required.

Selection criteria

We included (RCTs) comparing three interventions categories: antiresorptive drugs, bone-anabolic drugs, or placebo. Participants should be CD patients with osteoporosis or osteopenia. The primary outcome should include the change of lumbar spine bone mineral density (LSBMD), as measured by dual energy x-ray absorptiometry (DXA). We excluded (i) non-RCTs; (ii) review and meeting abstract; (iii) non-medicine intervention and mixed with ulcerative colitis (UC), or indeterminate colitis (IC) patients. There was no restriction on language, age, or sex.

Data extraction

Two authors independently reviewed abstracts and full texts of all eligible studies. We extracted information from all eligible studies on author names, publication year, number of patients, mean age, interventions, follow-up time. Any discrepancies in data extraction were resolved by the corresponding author.

Outcome measures

Primary outcome was improving LSBMD at 12 months (if the definite BMD value at 12 months was not available, the last time-point assessment in the trial would be taken); secondary outcome was the change of total hip BMD (THBMD) at 12 months (if the direct results of THBMD was not available, the value of total hip components would be taken, such as femoral neck, femoral trochanter) and incidence of adverse events (AEs).

Statistical analysis

We used WinBUGS (version 1.4.3, MRC Biostatistics Unit, Cambridge, UK) with a random effects model proposed by Chaimani (obtained from www.mtm.uoi.gr) within a Bayesian framework to perform network meta-analyses. The posterior parameters were calculated by Markov chain Monte Carlo methods in the network meta-analysis.[ Non-informative uniform and normal prior distributions were performed and an automatically generated starting value was utilized to fit the model.[ After an initial burn-in of 50,000, we conducted another 250,000 iterations. For dichotomous variables and continuous variables, treatment effects were summarized as odds ratio (OR) and standardized mean difference (SMD) with their corresponding 95% confidence interval (CI) or credible interval (CrI) (95% CI for direct evidences and CrI for indirect evidences), respectively. To make the ranking of the treatments for each outcomes, we calculated the surface under the cumulative ranking area (SUCRA).[ For the primary outcomes, lower SUCRA value indicated better rank of the treatment, whereas it was opposite regarding the secondary outcomes.

We compared the pooled OR and SMD from network meta-analysis with the corresponding OR or SMD from traditional pair-wise meta-analysis to evaluate whether there was inconsistency between direct and indirect evidences. The traditional pair-wise meta-analysis was conducted for trials that directly compared different interventions using Stata (version 13.0, StataCorp, College Station, TX). To account for heterogeneity, the method of DerSimonian and Laird random effects model was used.[ The heterogeneity between eligible studies was explored by Chi-square based Q-test in traditional pair-wise meta-analysis and the presence of heterogeneity was considered significant if P < 0.1.[ Nevertheless, inconsistency was also quantitatively assessed by calculating inconsistency factors and their 95% CI in closed loops. The goodness-of-fit of the model was examined by calculating the posterior mean residual deviance, and the model was considered to fit the data well when the posterior mean residual deviance approximated the number of data points in the present study.[

To detect the small study effects of the data, we conducted comparison-adjusted funnel plots.[ Sensitivity analysis was performed to authenticate the stability of the results according to the quality of study (excluding studies with a high risk of bias) and follow-up time (excluding studies with a short-term follow-up less than 12 months). This study was conducted and reported in accordance with the PRISMA guidelines.[

Assessment of risk of bias

We independently assessed the methodologic quality of included trials using the Cochrane Collaboration tool,[ which included the following items: (i) random sequence generation; (ii) allocation concealment; (iii) blinding of patients and personnel; (iv) blinding of outcome assessment; (v) incomplete outcome data; (vi) selective reporting; and (vii) other bias.

Quality of evidence

A four-step approach was used to rate the quality of therapeutic effect estimates from network meta-analysis based on the grading of recommendation, assessment, development, and evaluation (GRADE).[ Evidence evaluation includes direct, indirect as well as network estimates and it is rated as high, moderate, low and very low quality. At the beginning of the assessment, the quality of direct evidence was considered high and it may be rated down for the following reasons: risk of bias; inconsistency; indirectness; imprecision; publication bias. The rating for indirect evidence from the lower rating of the quality of direct evidence would be further rated down due to imprecision and indirectness. We used the higher rating of direct and indirect evidence as the quality rating for the network estimates.

Results

Characteristics of included studies

The flow diagram of study selection is provided in Fig. 1. First, 1682 literature citations were identified. After duplicates’ removal, a total of 1314 citations were yielded. After reading titles and abstracts, 1277 ineligible records were excluded, leaving 37 articles. Furthermore, 25 studies were excluded for the following reasons: non-RCTs, review article, meeting abstract, non-medicine intervention, and mixed with UC or IC. Finally, 12 eligible studies were included with a total of 920 subjects who received 1 of the 8 agents or placebo in this network meta-analysis. These studies come from different countries, 5[ were reported in Germany, 3[ from Canada, and 4[ from Netherlands, Ireland, Danmark, and UK, respectively. The mean age of participants in trials ranged from 14.6 to 45.1 years.

Figure 1

Flow diagram of study selection. A total of 12 studies were included in this network meta-analysis after literature search and selection.

Among these 12 eligible studies, there were 2-arm (n = 10) or 3-arm (n = 2) RCTs. Four studies compared ibandronate and sodium-fluoride with each other or against placebo (head to head trials: n = 3). Two studies compared risedronate with placebo with same dosage schedule. Likewise, 2 studies compared zoledronate with placebo. The remaining 4 articles were compared etidronate, pamidronate, alendronate, phylloquinone with placebo, respectively (Table 1).

Table 1

Characteristics of included studies.

Risk of bias in included studies

The risk of bias in all included studies is illustrated in Supplementary Figure S1. Nine studies provide sufficient details of randomization. Three studies revealed inadequate in terms of allocation concealment. Two studies did not blind to participants or study personnel and 1 trial had a high risk of bias due to unblinding of outcome assessment. Considering incomplete outcome data domain, 9 studies were cited as a low risk of bias. Eight studies had a low risk of reporting bias.

Outcomes

Efficacy: The improvement of LSBMD

The comparisons on changes of LSBMD with various medical therapies in CD patients by network meta-analysis are shown in Supplementary Figure S2 (A). The 12 RCTs were included in our study enrolled 920 patients (Table 1), Of these patients, 353 participants were assigned to placebo therapy,[ 149 to sodium-fluoride therapy,[ 122 to ibandronate therapy,[ 109 to risedronate therapy,[ 72 to etidronate therapy,[ 37 to pamidronate therapy,[ 36 to phylloquinone therapy,[ 27 to zoledronate therapy,[ 15 to alendronate therapy.[

Results of comparisons on the change of LSBMD in various treatment strategies in our network meta-analysis are shown in Fig. 2(A). Compared with placebo, zoledronate (SMDs 2.74, 95% CrI 1.36–4.11) reached statistical significance in improving LSBMD of CD patients, followed by sodium-fluoride (SMDs 1.23, 95% CrI 0.19–2.26). Similarly, zoledronate (SMDs 2.49, 95% CrI 0.21–4.78 and SMDs 2.68, 95% CrI 0.42–4.96) showed significant difference with pamidronate and etidronate, respectively. Additionally, alendronate (SMDs 1.23, 95% CrI −0.66 to 3.14), risedronate (SMDs 0.15, 95% CrI −1.13 to 1.44), ibandronate (SMDs 0.79, 95% CrI −0.37 to 1.95), etidronate (SMDs 0.05, 95% CrI −1.77 to 1.87), and pamidronate (SMDs 0.24, 95% CrI –1.59 to 2.07) except phylloquinone (SMDs −0.16, 95% CrI −1.99 to 1.67) all had a tendency to improve LSBMD in CD patients compared with the placebo group, but there were no statistical significance. In addition, zoledronate was significantly superior to risedronate (SMDs −2.59, 95% CrI −4.45 to −0.72) and phylloquinone (SMDs −2.89, 95% CrI −5.17 to −0.60) in improving LSBMD of CD patients.

Figure 2

Efficacy and safety of agents in Crohn disease patients with low bone mineral density (BMD): the efficacy was estimated in 2 lower triangle comparing column-defining with row-defining treatments (A, the change of lumbar spine BMD (LSBMD); B, the change of total hip BMD [THBMD]). The efficacy was estimated in upper triangle (A) comparing row-defining with column-defining treatments. For efficacy and safety assessment, treatment effects were summarized as standardized mean difference (SMD) and odds ratio (OR) with their corresponding 95% credible intervals, respectively. For change of BMD, SMD lower than 0 favor the row-defining treatment while for adverse effects, ORs lower than 1 favor the column-defining treatment.

SUCRA expressed as percentages (Fig. 3) were: 2.5% for zoledronate; 27% for sodium-fluoride; 32% for alendronate; 44% for ibandronate; 62% for pamidronate; 66% for risedronate; 68% for etidronate; 74% for phylloquinone; 75% for placebo. According to the synthetical data from SUCRA values, zoledronate (SUCRA = 2.5%) might have the highest probability to be the best treatment for improving LSBMD in CD patients and phylloquinone seemed to be the worst one within all available interventions (74%).

Figure 3

Surface under the cumulative ranking curve expressed as percentages ranking therapeutic effects of treatments for Crohn disease patients with low bone mineral density (BMD). The surface under the cumulative ranking area of each bar represents the probability size of each treatment. The color of the bar represents the specific therapeutic schedule. The abscissa represents the therapeutic probability of various drugs and the ordinate denotes the outcomes about the change of BMD and the incidence of adverse events. For the change of BMD, the pharmacological agent with the longest bar means the most efficacious treatment while agent with the shortest bar indicates the worest one. Oppositely, for the incidence of adverse events, the pharmacological agent with the shortest bar means the most safety treatment while agent with the longest bar indicates the worest one. LSBMD = lumbar spine bone mineral density, THBMD = total hip bone mineral density.

Efficacy: The change of THBMD

There were 7 trials introduced the improvement of THBMD as an assessment of drug efficacy. The network of included medical therapies comparisons is shown in Supplementary Figure S2 (B). A total of 280 CD patients with low BMD were assigned to placebo group,[ 15 to alendronate group,[ 109 to Risedronate group,[ 36 to phylloquinone group,[ 36 to sodium-fluoride group,[ 35 to ibandronate,[ 72 to etidronate group,[ 37 to pamidronate group[ (Table 1).

Network results regarding the change of THBMD are shown in Fig. 2(B). Most of the included drugs revealed an increase of THBMD when compared with placebo, unfortunately, no obviously significance could be identified (alendronate: SMD 0.49, 95%CrI −0.63 to 1.62; risedronate: SMD 0.27, 95%CrI −0.44 to 1.00; sodium-fluoride: SMD 0.03, 95%CrI −1.07 to 1.12; ibandronate: SMD 0.07, 95%CrI −1.02 to 1.15; pamidronate: SMD 0.21, 95%CrI −0.82 to 1.24).

The value of SUCRA as shown in Fig. 3 was partially consistent with the results of primary outcome, and alendronate (SUCRA = 24%) indicated the highest probability of being the most efficacious agent for improving THBMD, nevertheless, phylloquinone seemed to be the worst one (SUCRA = 72%).

Safety: Adverse events (AEs)

The comparisons on the AEs in our network are shown in Supplementary Figure S2 (C). Nine 2-arm and 2 3-arm studies comparing sodium-fluoride (n = 4), ibandronate (n = 3), etidronate (n = 1), risedronate (n = 1), pamidronate (n = 1), phylloquinone (n = 1), zoledronate (n = 2), or alendronate (n = 1) with each other or against placebo were included in the network meta-analysis. A total of 262 CD patients with low BMD were assigned to placebo group,[ 134 to sodium-fluoride group,[ 117 to ibandronate group,[ 49 to etidronate group,[ 45 to risedronate group,[ 37 to pamidronate group,[ 36 to phylloquinone group,[ 27 to zoledronate group[, and 15 to alendronate group[ (Table 1).

Results of comparisons on the incidence of AEs in various treatment strategies in this study are shown in Fig. 2(A). Although the occurrence of AEs in the RCTs demonstrated no statistical difference between agents and placebo, the incidence of AEs in alendronate (OR 0.13, 95% CrI 0.01–1.51), risedronate (OR 1.68, 95% CrI 0.27–10.16), phylloquinone (OR 0.30, 95% CrI 0.02–3.59), zoledronate (OR 0.66, 95%CrI 0.15–2.77), sodium-fluoride (OR 0.81, 95% CrI 0.36–1.79), ibandronate (OR 1.06, 95% CrI 0.39–2.78), etidronate (OR 0.35, 95% CrI 0.04–2.18), and pamidronate (OR 0.70, 95% CrI 0.11–4.60) compared with placebo, which had a certain trend to be higher than placebo.

Simultaneously, the corresponding SUCRA values (provided the optium treatment decision for clinicians) is depicted in Fig. 3. The values were 77% for risedronate, 70% for ibandronate, 68% for placebo, 54% for sodium-fluoride, 52% for pamidronate, 51% for zoledronate, 32% for phylloquinone, 31% for etidronate, 16% for alendronate, respectively. According to the synthetical evaluation of SUCRA values, risedronate (SUCRA = 77%) had the highest probability of being the safest medicine for increasing LSBMD of CD patients and alendronate ranked the worst (16%).

Evaluation of inconsistency and fit of the models

The results of the pair-wise and corresponding Bayesian network meta-analysis are shown in Fig. 4, Supplementary Figure S3 and Figure S4. The effect size and relevant CI or CrI delivered no obvious discrepancy between the 2 different types of comparisons, indicating that there were no inconsistencies. Moreover, the consistency was also confirmed by the quantitative assessment in closed loops (Table 2). The result of the model test indicated that the posterior mean residual deviance approximated the data points in both primary and secondary outcomes (Table 2), namely, the present model fit the data well.

Figure 4

Forest plot with 95% confidence interval (CI) or credible intervals (CrI) (95% CI for direct evidences and CrI for indirect evidences) for continuous variable (change of bone mineral density in lumbar spine) by traditional meta-analysis and Bayesian network meta-analysis. Ale = alendronate, Ris = risedronate, Phy = phylloquinone, Zol = zoledronate, Sod = sodium-fluoride, Iba = ibandronate, Eti = etidronate, Pam = pamidronate, Pla = placebo.

Table 2

Heterogeneity and evaluation of model fit in included.

The GRADE approach was applied to the primary outcomes about the change of BMD in lumbar spine (LS) at 12 months (Table 3). The quality of direct and indirect evidence was very low, low, or moderate for all comparisons. Likewise, the quality for network meta-analysis was also very low, low, or moderate.

Table 3

Estimates of effects and quality ratings for comparison of drugs for low BMD in CD patients.

Publication bias and sensitivity analysis

The result of the comparison-adjusted funnel plots did not uncover any evidence of apparent asymmetry (Supplementary Figure S5). Sensitivity analysis regarding the quality of study and follow-up time did not significantly alter the results of 3 outcomes (Supplementary Table S1 and Table S2).

Discussion

In this network meta-analysis, we included 12 RCTs comparing the efficacy and safety of 8 agents for the management of low BMD in CD patients: 3 studies used ibandronate, 2 studies used zoledronate, 2 studies used risedronate, 4 studies used sodium-fluoride, the rest 4 studies used etidronate, pamidronate, alendronate, and phylloquinone, respectively. The efficacy and safety of these agents were analyzed as the primary and secondary outcomes, respectively. We provided hierarchies for clinicians in the treatment process. Our results suggested that zoledronate might have the highest probability to be the best treatment for increasing LSBMD in CD patients among all included agents. The comprehensive analysis of the secondary efficacy assessment concerning the change of THBMD revealed no obvious significance between all agents and placebo, and alendronate yielded a trend to be the best agent for improving THBMD. For the safety assessment, risedronate showed the greatest power to decrease the risk of AEs while alendronate performed the worst, but all of them were not statistically sound, possibly due to small sample size.

Previous studies indicated that bisphosphonates had beneficial effect for low BMD among postmenopausal women or patients without CD who were undertaken glucocorticoid administration.[ Stokkers et al[ reported that 49 patients with IBD (44 CD patients) received injection of pamidronate, the BMD of lumbar spine was increased. Similarly, Bartram et al[ also confirmed that intravenous pamidronate had a therapeutic effect in CD patients with low BMD. The study by Sbrocchi et al[ demonstrated that zoledronate lead to a significant increase in LSBMD at 6 and 12 months with a well-tolerated infusion. Although the above studies showed that several drugs were effective in CD patients with low BMD, there still lacks study to show which agent specifically targets on bone health of CD patients. Therefore, we performed this network meta-analysis to evaluate the efficacy and safety of different agents of antiresorptive drugs (bisphosphonates), bone-anabolic drugs (sodium-fluoride and phylloquinone) for treatment of low BMD in CD patients.

We found that a single dose of intravenous zoledronate might have the highest probability to be the most beneficial strategy for increasing LSBMD of CD patients, followed by sodium-fluoride and both reached statistical significance. The outcome was consistent with the confirmed beneficial effect of a single dose zoledronate infusion in a study conducted by Reid et al.[ Actually, it is generally believed that oral bisphosphonate is associated with low absorption and poor patients’ compliance due to its side effects in an empty stomach in CD patients with low BMD whose digestive tracts have been compromised,[ and our results were consistent with the statement that zoledronate was better than other bisphosphonate. Interestingly, only zoledronate had statistical beneficial effect of therapy for lumbar spine bone density in this network meta-analysis while other types of bisphosphonate rendered no significant difference when compared with placebo. Several possible reasons to explain this result as follows. First, small cohorts in our study may not have enough powerful effect size to show a statistical difference between the treatment and placebo group. Second, few articles were published to demonstrate the efficacy of etidronate, pamidronate, alendronate, and phylloquinone. Third, patients in different RCTs have different severity of low bone density.

The results obtained by the network meta-analysis demonstrated that sodium-fluoride had a beneficial effect in increasing bone density of CD patients, which was consistent with previous studies on postmenopausal women reported in a traditional meta-analysis.[ This traditional meta-analysis confirmed that the incidence of new nonvertebral fractures and gastrointestinal adverse events were higher than those in sodium-fluoride group especially with high doses and in a non-slow-release form than placebo group at the forth year and reached statistical difference, thus, the use of sodium-fluoride was not recommended because the higher dose of sodium-fluoride, the higher incidence of nonvertebral fractures, and gastrointestinal tract side effects. Also, it had no effects on the reduction of the incidence of vertebral fracture. In view of our results, the long-term use of sodium-fluoride in clinical practice should be with a bit more cautions.

Besides the spine, a markedly increased rates of hip fracture could also be observed among patients with CD or other types of IBD.[ Thus, additionally, we induced an assessment of BMD change within total hip to further detect the effect of included treatments. However, the final results were unconspicuous. To determine the potential role of etidronate in the improvement of BMD, Siffledeen et al[ conducted a RCTs recruiting 154 CD patients with a follow-up of 24 months. They proposed an interesting result that both the etidronate and the non-etidronate treated groups showed a statistical trend of BMD improvement compared with the baseline state at all involved bones, except the total hip. Although the improvement was probably attributed to the utilization of supplemental calcium and vitamin D, the total hip might render an insensitive to these exogenous agents and other supplements. This could be partially interpreted by the difference of surface-to-volume ratio and vascularity between trabecular bone and cortical bone.[

The incidence of AEs in bisphosphonates, sodium-fluoride, and phylloquinone concomitant has a trend to be higher than placebo, but no significant differences were identified. Several factors could interpret this result: first, the overall 12 RCTs in our network meta-analysis was not powerful enough to demonstrate statistical significance difference between intervention and placebo group; second, few articles were published to demonstrate the safety of etidronate, pamidronate, alendronate, and phylloquinone.

Some strengths in our article: (i) it is the first and the most comprehensive network meta-analysis to investigate different intervention in treating low BMD in patients with CD; (ii) design and outcome assessments were consistent in all trials included in our study; (iii) the latest guideline of GRADE approach was applied to rate the quality of evidence of eight agents for network meta-analysis. However, our article also has limitations: (i) due to the lack of data to yield outcomes (such as fracture rate) in most trails included in our article, we could only extract the change of LSBMD value to evaluate the efficiency of various treatment options; (ii) due to the small sample sizes in our study, most network comparisons were in low quality by GRADE assessment.

In conclusion, our network meta-analysis suggests that zoledronate might have the highest probability to be the best treatment for increasing LSBMD in CD patients among all other agents. For the safety assessment, risedronate showed the greatest power to decrease the risk of AEs while alendronate performed the worst, but all of them were not statistically sound. Future, more high quality RCTs (larger number of CD patients with low BMD and longer follow-up time) are needed to make head-to-head comparisons of treatments.