Skip to main content

Non-steroidal anti-inflammatory drugs and bone healing in animal models—a systematic review and meta-analysis

Abstract

Background

Non-steroidal anti-inflammatory drugs (NSAID) have excellent anti-inflammatory and analgesic properties and are extensively used to treat post-traumatic or surgical musculoskeletal pain. Although an extensive literature exists on the administration of NSAID on animal bone healing, no systematic review and meta-analysis of animal studies that investigate the effect of NSAID administration on bone fracture healing. Objective of this study is to conduct a systematic review and meta-analysis to estimate the effect of NSAIDs administration on bone healing biomechanical and histomorphometric measurements in different animal models after bone fracture surgery.

Methods

We performed a systematic review and meta-analysis of animal studies to estimate the effect of NSAID administration after bone fracture on healing outcomes. We searched eight databases without limiting the search to starting date up to 1 February 2021 for articles on fractured bone healing in animal models in which NSAID were administered.

Results

Out of 6732 articles screened, 47 were included and 3 common bone healing outcomes were analysed: biomechanical properties (maximum force to break, stiffness, and work-to-failure), micro-computed tomography (μ-CT), and histomorphometric measurements. The studies were generally of low-quality scores because crucial information, especially concerning randomization, blinding, and allocation concealment, was poorly reported. Our results show that the negative effects of NSAID after bone fracture on certain biomechanical properties of the healing bones was not statistically significant in mice compared with other animals, in females compared with males, and in younger compared with older animals.

Conclusion

The findings demonstrated that NSAIDs administration decreased the biomechanical properties of healing bones after fracture surgery in comparison to the control group. Moreover, different effect on certain outcomes was detected among different sites, sex of the animals, and the time of assessment.

Trial registration

Protocol published and registered in SYstematic Review Center for Laboratory animal Experimentation (SYRCLE) in 2017, https://www.radboudumc.nl/getmedia/757ec408-7a9e-4635-8233-ae951effea54/Non-Steroidal-Anti-inflammatory-Drugs-and-bone-healing-in-animal-Models-Systematic-Review-and-Meta-Analysis.aspx

Peer Review reports

Background

Non-steroidal anti-inflammatory drugs (NSAID) have been extensively used to treat post-traumatic and surgical musculoskeletal pain because of the excellent anti-inflammatory and analgesic properties of cyclooxygenase (COX) inhibitors [1,2,3,4,5]. However, clinical and biological data suggest that COX inhibition has a negative impact on bone tissue repair. Results from several animal trials have suggested that decreased bone healing, including biomechanical and histomorphometric properties, was associated with NSAID administration after bone fracture [6,7,8,9,10,11,12]. While few studies have investigated the effect of NSAIDs on bone healing in humans, their results are contradictory, possibly reflecting methodological limitations [13,14,15,16]. Indeed, authors have argued for well-designed, large, multicentric randomized controlled trials with appropriately defined endpoints [17]. To the best of our knowledge, there are no clinical studies to accurately determine the effect of NSAID on bone healing.

Several animal models including different animal strains, sex, and fracture techniques have been used to test the effect of types and administration durations of NSAID on bone healing outcomes [10, 18,19,20,21,22,23,24]. The assessments of bone healing in these studies include integrated multilevel measurements at the organ (e.g., biomechanical to tissue levels using micro-computed tomography (μ-CT) analysis), cellular (e.g., histology, histomorphometric analysis), and gene (e.g., mRNA microarray to explain the associated healing pathway) levels [25, 26]. Overall, the techniques used to assess bone healing can be divided into four main categories: imaging analyses, biomechanical tests, detection of serologic markers, and clinical examinations [25]. Among these categories, the biomechanical confirmed by histomorphometric analysis is the best way to determine the success of fracture healing [27, 28]. These tests are widely used in animal studies at different time points of fracture healing, which are selected according to the biological, physiological, and pathological changes of the fractured bone [28,29,30,31].

Although an extensive literature exists on the administration of NSAID on animal bone healing, no systematic review and meta-analysis have yet been conducted on the subject. Such work is important because it can identify the key histomorphometric and biomechanics characteristics during the healing process as well as it provides information on different factors that may affect the healing process after NSAID administration (e.g., NSAID type, animal type and strain, sex of the animal, fracture type, time of assessment). This information may help to design future experimental trials and facilitate knowledge translation. Therefore, we performed a systematic review and meta-analysis of all identified and available animal studies on the effect of NSAID administration after bone fracture on healing outcomes.

Specifically, our objectives were to estimate the extent to which the effect of NSAID administration after bone fracture using animal models: (i) results in less favorable bone-biomechanical and morphometric healing measurements; and (ii) differs by animal type, age and sex, type of NSAID, and length of follow-up.

Methods

The systematic review methodology was specified and documented in advance using SYstematic Review Center for Laboratory animal Experimentation (SYRCLE)’s protocol template for animal studies and registered in the SYRCLE database (Supplementary material Table S1) [32, 33]. In our initial protocol, we had proposed to carry out a sensitivity analysis to assess whether the methodological quality of the studies included in the meta-analysis greatly influences the findings of the review. However, most of the studies in our systematic reviewed received low-quality scores; we therefore amended the review protocol to remove this part of the analysis. In addition, we mentioned using Holm-Bonferroni method in our protocol, but we did not apply it in our analysis.

Search strategy and selection of studies

We searched eight databases (Embase, Scopus, MEDLINE, CINAHL, BIOSIS, Cochrane, Central, and DARE) for original articles concerning the effects of NSAID on fractured bone healing in animal models. We considered the period from the starting date of each database to 1 February 2021. The main terms used in the search strategy, developed with the help of the Liaison Librarian for Life Sciences at McGill University, were “anti-inflammatory agents,” “non-steroidal,” “bone,” and “animals” (the complete search strategy is provided in Supplementary material Table S2). We used a search filter to detect animal studies and exclude human studies. Although we did not impose language restrictions while searching, only English articles were reviewed. Also, we did not include conference abstracts because they do not provide sufficient data to allow for an evaluation. No other restrictions were used. Additionally, we hand searched the reference lists of eligible articles, and screened review articles for relevant references. The first selection was performed based on independent selection by two reviewers (HA, AP) using RAYYAN software (httpp://rayyan.org, Doha, State of Qatar) [34]. Any disagreement was solved by discussion or by including the third reviewer (BN). Full-text articles were screened and included when they met the following pre-specified criteria: (1) a controlled NSAID interventional design after bone fracture; and (2) a description of outcome measures related to bone healing (biomechanical characteristics, μ-CT scan measures, radiographic bone assessment, and/or histomorphometric and histological based grading). Papers were excluded if they fulfilled one of the following criteria: not an original article (e.g., review or letter), use of a bone graft or other material, included only outcomes that were not biomechanical or histomorphometric, not NSAID used, not experimental animal studies, and duplicate studies. A list of articles excluded as well as the reasons for exclusion are available from the authors upon request.

Data extraction

The final set of articles was assessed independently by two reviewers who extracted the data using DistillerSR (Evidence Partners, Ottawa, Canada) following a piloted data extraction method. Disagreements were resolved by consensus.

Data retrieved from the articles included characteristics of the animals (e.g., animal model used, weight and sex), test methods (bone fracture [e.g., site, type, number and technique used to perform the fracture, whether there was fixation or not and type of fixation method]), use of medication (e.g., opioid, antibiotics, and NSAID [e.g., type, name, duration of use and route of administration]), and outcome data (time points of outcome measures collection, type of outcome). We grouped the outcomes into main classes: (i) biomechanical (e.g., maximum force or ultimate force load, stiffness, work-to-failure), (ii) μ-CT assessment of healing (volume and density), and (iii) histomorphometric characterization of fracture callus (bone, cartilage, and mineralized tissue).

Raw data or group averages (mean, median), standard deviation (SD), standard error (SE), or ranges and number of animals per group (n) were extracted for all continuous outcome measures. We contacted the authors to obtain original data if results were presented graphically or were incomplete. If data could not be retrieved, the study was excluded from further analysis.

Risk of bias assessment

Two blinded reviewers (HA, AM) assessed the internal validity of the included studies using SYRCLE’s risk of bias tool (Table S3) [35]. The tool, an adaptation of the Cochrane risk-of-bias tool, considers aspects of bias specific to animal studies. It contains 10 entries related to 6 types of bias (selection, performance, detection, attrition, reporting, and other bias). The score (yes) indicates a low risk of bias, (no) indicates a high risk of bias, and (?) indicates an unclear risk of bias. We were concerned that many items would be rated as having an unclear risk of bias because of the known poor reporting of experimental designs [36]. To overcome this problem, we added four entries to the tool, pertaining to randomization, blinding, sample size calculation, and time of day of the NSAID administration or time of day at which surgery was performed [35]. For these items, ‘yes’ and ‘no’ indicates reported and unreported, respectively.

Data analysis and synthesis

The meta-analysis was conducted according to Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines, and using Comprehensive Meta-Analysis software (version 2.2.064, Biostat Inc., Englewood) when five or more independent comparisons from at least three different studies per outcome category were included (provided that outcome measure assessments were sufficiently comparable). We calculated the standardized mean differences (SMD) through Hedges g effect sizes [37]. The calculation was (SMD = the mean of the NSAIDs group minus the mean of the control vehicle group divided by the pooled standard deviations of the two groups) to account for the differences in the units of measurements. We used Hedges g effect to calculate the SMD, Hedge’s g (which is based on Cohen’s D but includes a correction factor for small sample size bias) [38, 39]. Again, the calculations need to take into account the direction of effect.

Despite the anticipated heterogeneity, the individual effect sizes were subsequently pooled to obtain an overall SMD and 95% confidence interval (95%CI) [38, 40]. We used a random-effects model [40], which takes into account the precision of individual studies and the variation among them, and weights each study accordingly. If multiple independent experimental groups were compared to the same control group within the meta-analysis, the number of animals in the control group was corrected by dividing it by the number of experimental groups.

Rather than computing a single summary measure, an important objective of meta-analysis is to explore the sources of heterogeneity [41], a measure of the degree of variability in study results, and assess which variables influence the effect of NSAID on bone healing outcomes. We conducted subgroup analyses according to sex and animal species (mice, rats, and rabbits), type of NSAID (non-selective/COX-2 selective), type of fracture, and period of outcome measurement or data collection (early healing less than 21 days, 21 to 48 days, and more than 48 days). We present below the results for subgroups containing at least 10 comparisons. A minimum of three independent comparisons per subgroup was needed to record the subgroup characteristics. The interpretation of differences between subgroups should be used mainly to construct new hypotheses rather than drawing definite conclusions. Heterogeneity for subgroup analyses was assessed using I2 and the Q statistic.

The subgroup analysis was done and generated through three methods for comparing the effect size across subgroups. One method was to use a Z-test to compare the effect sizes directly. Another method was to use a Q-test to partition the variance and test the between-subgroups portion of the variance. A third one was to use a Q-test to assess the dispersion of the summary effect about the combined effect. All the methods used by the software assess the variance across subgroups effects relative to the variance within subgroups. The row tables for each subgroup analysis generated by the comprehensive meta-analysis software were added in the supplementary information file for detailed results based on the previous described methods (Tables S5 to S27). We assessed publication bias for different healing outcomes by evaluating the possible asymmetry using funnel plots Trim and fill and “Egger’ regression test” if the analysis contained at least 20 comparisons” [38]. All statistical analysis was conducted using Comprehensive Meta-Analysis Version 3 Software (NJ, USA) [42]. A P value less than 0.05 was considered to be statistically significant.

Results

Study selection and characteristics

After the full-text assessment, 47 publications were included in the systematic review (see Fig. 1 for the PRISMA flowchart). The authors were contacted when we could not retrieve or understand the data and only two out of ten responded to the request and sent the raw data. Table S4 presents characteristics of the included studies (the complete list is available in Supplementary file 1). Overall, study characteristics varied considerably; most studies were performed in rats (31 studies; 66%), 5 in mice (10.6), 9 in rabbits (19.1%), and 2 in dogs. Eight studies did not report the sex of the animal, while 26 (55.3) and 13 (27.7%) used only male or female animals, respectively, and none used both sexes. Eleven studies (23.4%) used a selective COX-2 NSAID as an experimental intervention drug, 22 (46.8%) and 14 (29.8%) used a non-selective or both NSAID types, respectively. There was a great variability on the primary outcomes for biomechanical characteristics; 37 (78.7%) studies reported 1 or more biomechanical characteristics, and 10 (21.2%) studies did not report any of these characteristics.

Fig. 1
figure 1

PRISMA flow diagram of the selection of the studies for the systematic review and meta-analysis included

Risk bias and quality of the studies

The assessment results for risk of bias and quality of reporting related to randomization, blinding, sample size calculation, and time of day for NSAID administration or surgery are summarized in Figs. 2 and 3. Scores for each study are presented in Supplementary file 2. Among the 47 included studies, 30 (63.8%) mentioned the term “randomization” at any step in the study, but no article provided details on the method used. Only 14 (29.8%) studies reported blinding which for most of them was on the histological outcome assessment. Among all included studies, only 6 (12.8%) reported a sample size calculation; 1 article specified the time of day at which NSAID was administered or the time that surgery was performed (day or night). Due to poor reporting, many items evaluating the risk of bias on the assessment tool showed an unclear score. For example, “selective outcome reporting bias” was assessed as unclear for all studies because none reported using a research protocol defining primary and secondary outcomes.

Fig. 2
figure 2

Risk of bias assessment using SYRCLE risk of bias assessment tool

Fig. 3
figure 3

Quality of reporting assessment of the included studies

For more details regarding each included study (n = 47) in the systematic review and the risk of bias assessment for each parameter is presented in the Supplementary file 2 including how many studies with high, low, or unclear for each parameter of the tool.

Meta-analysis of NSAID administration during fractured bone healing

Due to missing information in 11 studies regarding outcome data, or not suitable outcome measurement, or the intervention is not fracture, or the drug was not NSAID or was given with other analgesic such as opioid or steroidal anti-inflammatory for the same group during the course of the study [6, 43,44,45,46,47,48,49,50], we included 36 studies in the meta-analysis. Thirty-two studies compared the effect of administration of one or more NSAID on biomechanical characteristics (e.g., maximum force (MF) to fracture, stiffness, and work-to-failure) to a control group. For three-point mechanical bending properties, the analysis includes 186 experiments covering different animal models, NSAID types, and measurement time points. Four and seven studies were included in the analysis of the effect of NSAID administration on the μ-CT and histological assessment healing outcomes, respectively. The average timing of data collection after bone fracture to assess the mechanical bending maximum force of healing bones was an average of 29.6 days (minimum, 5 days; maximum, 84 days).

Biomechanical assessment

Results from 30 studies including 94 comparisons showed that the maximum force to fracture was significantly decreased, indicating bone healing delay, in animals that received an NSAID after bone fracture compared to the control group (SMD − 0.58, 95%CI [− 0.74, − 0.42]; Table 1). Heterogeneity was moderate (I2, 55.04%). Similarly, animals that received NSAID had an overall decrease in bone stiffness and work-to-failure properties (SMD − 0.56 [− 0.76, − 0.37] and SMD − 0.58 [− 0.95, − 0.20]) respectively compared to controls (Fig. 4; Table 1). Between-study heterogeneity was moderate for both stiffness (I2, 60.41%) and work to failure outcomes (I2, 56.29%).

Table 1 Meta-analysis based on several subgroups showing the effect of NSAID administration on three-points mechanical bending measurements (maximum force, stiffness, and work to failure) of healing bones after fracture
Fig. 4
figure 4

Forest plot of the included studies (experimental groups), which used three-point mechanical bending a maximum force (MF), b stiffness, c work-to-failure. The forest plot displays the standard mean differences (SMDs) Hedges’ g, 95% confidence interval. The diamond indicates the overall estimation and its 95% confidence interval

We explored the sources of heterogeneity by examining the effect sizes in predefined subgroups: animal sex, age and species, time of bone collection, and type of fractured bone. While animal age and type of bone were source of heterogeneity for the maximum force to break, time of sample collection and animal sex, age, and species were for the stiffness analysis. Moreover, sex and time of sample collection were sources of heterogeneity in the work to failure analysis (Table 1).

Table 1 shows the subgroup analysis for three-point mechanical bending measurements. For maximum force measurement, NSAID administration did not delay bone healing among mice (SMD − 0.28 [− 0.68, 0.10]) but did it in other animals. In addition, we observed a difference in this measurement for the subgroup analysis of bone model, while femur (SMD − 0.68 [− 0.88, − 0.48]) showed a significant difference between NSAID and control, tibia did not (SMD − 0.19 [− 0.49, 0.10]). Moreover, when comparing SMD across bone models, the effect of NSAID administration was significantly larger in femur compared to tibia (P = 0.007; Fig. 5b).

Fig. 5
figure 5

Effect of a animal models’ characteristics, b type of fractured bone, c time of collection healing bones, d type of NSAID, and e age of the animals on biomechanical bending measurements (maximum force to fracture) (MF) after administration of NSAID compared to the administration of a control vehicle. The columns indicate the effect estimate with the 95% confidence interval of the subgroups. SMD, standard mean difference–Hedges’ g mean. NM, not mentioned; NS-COX, non-selective cyclooxygenase inhibitor; COX2, selective-cyclooxygenase2 inhibitor

Bone stiffness among mice (SMD − 0.07 [− 0.55, 0.40]) and animals older than 16 weeks (SMD − 0.31 [− 0.82, 0.18]) did not differ between NSAID and control groups (Table 1). However, compared to controls, bone healing was better in mice taking NSAIDs than in rabbits (P = 0.01; Fig. 6a). The effect of NSAID administration was significantly different when the bone samples were harvested before 21 days compared to other time points between 21 and 48 days after surgery (P = 0.03; Fig. 6c).

Fig. 6
figure 6

Effect of a animal models’ characteristics, b type of fractured bone, c time of collection healing bones, d type of NSAID, and e age of the animals on biomechanical bending measurements (stiffness) after administration of NSAID compared to the administration of a control vehicle. The columns indicate the effect estimate with the 95% confidence interval of the subgroups. SMD, standard mean difference–Hedges’ g mean. NM, not mentioned; NS-COX, non-selective-cyclooxygenase; COX2, selective-cyclooxygenase 2 inhibitor

Regarding work to failure, there was no significant effect of NSAID administration on bone healing in the groups of female animals (SMD − 0.09 [− 0.55, 0.36]), those that received selective-cyclooxygenase2 NSAID (SMD, − 0.58 [− 1.26, 0.09]), and for the femur bone model fracture (SMD, − 0.38 [− 0.83, 0.05]) (Table 1; Fig. 7).

Fig. 7
figure 7

Effect of a animal models’ characteristics, b type of fractured bone, c time of collection healing bones, d type of NSAID, and e age of the animals on biomechanical bending measurements (work to failure) after administration of NSAID compared to the administration of a control vehicle. The columns indicate the effect estimate with the 95% confidence interval of the subgroups. SMD, standard mean difference–Hedges’ g mean. NM, not mentioned; NS-COX, non-selective-cyclooxygenase; COX2, selective-cyclooxygenase 2 inhibitor

Micro-computed tomography assessment (bone assessment)

We included five comparisons from four studies that measured healing bone using a μ-CT scan in the meta-analysis. The average time of bone collection after animal euthanasia was 19.5 days (range, 17–21 days). Figure 8 and Table 2 show the distribution of the data. Although the subgroup analyses were not performed because the number of comparisons was small, the overall analysis shows a significant difference in bone volume measurements for animals that received NSAID compared to controls (SMD, − 1.63 [− 2.87, − 0.39]), but this was associated with high heterogeneity among the studies (I2 83.32, P < 0.001).

Fig. 8
figure 8

Forest plot of the included studies (experimental groups), which used micro CT analysis bone volume and bone density measurements. The forest plot displays the standard mean differences (SMDs) Hedges’ g, 95% confidence interval. The diamond indicates the global (overall) estimation and its 95% confidence interval

Table 2 Meta-analysis based on several subgroups showing the effect of NSAID administration on bone volume and density measurements (μ CT assessment) of healing bones after fracture

Histomorphometric assessment

Seven studies including 33 experimental comparisons between NSAID administration and a control group showed no significant difference in all three (callus size, cartilage, and bone tissue) histomorphometric measurements (SMD, − 0.16 [− 0.49, 0.17], I2 = 54.64). Animal models and types of fractured bones were sources of heterogeneity for histomorphometric measurements of healing bones among studies (Fig. 9; Table 3). Interestingly, no mouse model was used to study the histomorphometric measurements related to bone, cartilage, or callus size. Rat models and histomorphometric evaluation at less than 21 days showed that bone healing was delayed in the NSAID group compared to controls (Table 3).

Fig. 9
figure 9

Forest plot of the included studies (experimental groups) that used histomorphometric analysis on microscopic image including callus size, cartilage tissue, bone volume, and bone area measurements. A forest plot displays the standards mean differences (SMDs) Hedges’ g, 95% confidence. The diamond indicates the overall estimation and its 95% confidence interval

Table 3 Meta-analysis based on several subgroups showing the effect of NSAID administration on bone volume measurements (histomorphometric analysis) of healing bones after fracture

Moreover, when comparing SMD across animal species, bone models, and time of collection, the effect of NSAID administration was significantly larger in rats compared to rabbits (P = 0.01; Fig. 10a), in femur compared to fibula (P = 0.02; Fig. 10b), and in the groups of bone samples that have been harvested less than 21 days (P = 0.03; Fig. 10c) after surgery.

Fig. 10
figure 10

Effect of a animal models’ characteristics, b type of fractured bone, c time of collection healing bones, d type of NSAID, and e age of the animals on histomorphometric measurements (callus, bone) after administration of NSAID compared to the administration of a control vehicle. The columns indicate the effect estimate with the 95% confidence interval of the subgroups. SMD, standard mean difference–Hedges’ g mean. NM, not mentioned; NS-COX, non-selective-cyclooxygenase; COX2, selective-cyclooxygenase 2 inhibitor

Publication bias

The possible presence of publication bias was observed when assessing the biomechanical bending and histomorphometric outcome measurements. The inspection of the funnel plot suggested asymmetry resulting from the underrepresentation of studies regarding the effect of NSAID (Fig. 11 a, and b). For biomechanical outcome, trim and fill analysis resulted in data points, indicating the presence of publication bias (Duval and Tweedi’s trim and fill for adjusted values; point estimate = − 0.36; 95%CI − 0.48 to − 0.23; Q value 630.69). The Egger’s regression test intercept result (95% CI − 2.33 to − 0.67); P < 0.001. For the histomorphometric outcome; trim and fill analysis resulted in data points, indicating the presence of publication bias (Duval and Tweedi’s trim and fill for adjusted values; point estimate = − 0.248; 95%CI − 0.469 to − 0.027; Q value 70.55). The Egger’s regression test intercept result (95% CI 1.67 to 5.8); P < 0.001.

Fig. 11
figure 11

Funnel plot for a biomechanical and b histomorphometric measurements of bone in the included studies

Discussion

This unique systematic review and meta-analysis was designed to answer a specific research question regarding the effect of different types of NSAID on bone fracture healing in animal models (in vivo studies). Three important outcomes commonly used to assess bone healing were analyzed: biomechanical properties (maximum force to break, stiffness, and work to failure), micro CT, and histomorphometric measurements.

NSAID administration had a negative effect on the biomechanical properties in different animal models of the included studies [10, 12, 18, 19, 48, 51,52,53,54,55,56,57,58,59]. However, the results for histomorphometric assessments did not show a difference (Table 3).

Depending on the type of NSAID, they are known to inhibit both COX isoforms. Both non-selective NSAID and COX-2-selective drugs decrease prostaglandin production, which plays an essential regulatory role in all phases of bone healing, especially the inflammatory phase [60,61,62,63]. Thus, it is reasonable to expect that NSAID administration after bone fracture may delay or impair healing outcomes [12, 53, 62].

As early as the 1970s, many animal trials strongly emphasized the negative effect of NSAID on bone healing [7, 9, 12, 19, 24, 45, 53, 59, 64]. Most studies that used rodents and rabbits demonstrated that non-selective and selective COX inhibitors impair the bone healing process [7, 9,10,11,12, 19, 24, 45, 53, 59, 64]. Conversely, only a few studies indicated that NSAID has little or no effect on fracture healing outcomes [20, 56, 57, 65]. However, these studies have significant limitations, because they tested only one time point, did not measure clear bone healing outcomes that include mechanical or histomorphometric analysis, or used very low NSAID doses, and some did not perform a proper statistical analysis [43, 66, 67]. Recently, a review by Huss et al. in 2019 stated the need for more evidence based and weighting the risk and benefit regarding administration of NSAID for orthopedic treatment and the majority of animal experimental data support short term and perioperative administration of NSAID after bone fracture [68]. Clinically, there are few retrospective studies and even fewer prospective clinical trials [13,14,15,16]. The results from the retrospective studies were contradictory. Some of them have confirmed the negative effect of NSAID on bone formation and healing following hip or femoral neck fracture as well as hip arthroplasty [2, 13], while others have shown no effect on bone healing. One of the few prospective clinical trials reported beneficial effects of NSAID on bone healing in humans, which was among Colles’ fracture patients who were treated with casting and reduction [14, 16]. Borgeat et al. 2018, demonstrated in their systematic review that results from the available human trials did not show strong evidence that NSAID administration is related to increase non-union after bone fracture. In addition, they emphasised on the need for further randomized clinical trials to support or refuse this hypothesis [69]. It is widely accepted that trying to understand the effect of NSAID administration on bone healing is extremely challenging in a clinical setting especially from a methodological perspective. For example, controlling the many confounding factors (e.g., smoking, diabetes, obesity) in a prospective manner requires considerable time and planning [17]. Therefore, it is imperative to translate the evidence from the available in vivo experimental studies. In fact, animal studies helped in understanding the physiological process of bone healing and can also provide important insights on the effect of NSAID on bone healing.

Some methodological issues that might hamper the interpretation of the experimental animal data and their subsequent translation to the clinical setting should be discussed. First, there was substantial heterogeneity among the various animal studies. We performed subgroup analyses to investigate factors that may modify the effect of NSAID on bone healing outcomes (e.g., animal species, sex, type of bone, age of the animal, type of NSAID, and time at which the outcome was measured. More studies are required, especially in mice, because contrary to other models such as rats, mice show no negative effect of NSAID on the stiffness of the harvested bone compared to the control group. Because the mouse genetic map is similar to that of humans, the mouse may be the best available model to study the effect of NSAID on bone healing in different human genetic conditions [70]. Overall, the pharmacokinetic variations between species, sexes, and ages should be considered, especially regarding drug absorption. Our meta-analysis showed non-significant negative effect of NSAID administration after bone fracture not only in mice compared to other animals, but also in females compared to males, and in younger compared to older animals.

Our results should be interpreted with the limitation of the included studies regarding higher risk of bias, using healthy animal model and the sample size calculation. Experimental studies that compare different sexes and ages within the same experiment are needed for stratification and comparison.

With the above-mentioned limitations, we observed that the negative effect of NSAID administration on biomechanical properties differs between animal species and this may suggest that rodent models (mice and rats) may be more sensitive than others for bone healing outcomes. Additionally, this finding suggests that the negative effect of NSAID on bone healing is species related for certain outcomes, and this need to be taken into consideration for knowledge translation. Other factors should be taking into consideration in selecting animal model for experimental bone regeneration studies which are highlighted in the 2018 systemic review by Peric et al. such as skeletal features of the selected animals, the model that can mimic the clinical scenarios including relevant doses and statistically supported sample size [71]. The timing of healing outcome measurements also seems to modify the results because they do not show a significant effect of NSAID administration on μ-CT and histomorphometric outcomes early in the healing process compared to the control groups, but results differ when measurements are taken after 28 days. It is important to consider this information during the design of further experimental protocols, and it may be crucial in managing research efforts and reducing the unnecessary use of animals.

Finally, we did not find any studies that used animals of both sexes in their experimental design, which seems to be important for future studies because results differed between males and females for the work-to-failure biomechanical outcome.

Methodological quality of the studies

Our study quality checklist assessed aspects of both internal and external validity, and we observed many studies were generally of low-quality scores and tended to overstate the effect size. The overall quality score accounted for a significant proportion of between-study heterogeneity; however, the correlation between the aggregate quality score and the effect was not clear. The reporting and risk of bias assessments indicate the need for protocol registration or publication, and for the reporting of the elements of randomization, blinding, and allocation concealment [35]. This concern is shared and addressed by others [33, 35, 38]. It is crucial that future animal studies improve the reporting of study procedures, allowing others to replicate and build on previously published work. With better reporting, systematic reviews of higher quality will also become feasible.

Limitations

Several limitations of this work should be considered. One important limitation of this study was the risk of bias analysis which showed most of the included animal studies are poorly reported. This lack of reporting important methodological details might lead to increase bias. This limitation needs to be considered during interpretation of our conclusion from the included animal studies. In addition, we observed a level of heterogeneity present among studies reviewed for certain measured outcomes. This is can be expected as animal studies are explorative and heterogenous with respect species, sex, design, and drug administration protocol compared to clinical trials [35]. Exploring heterogeneity among animal studies is one of the objectives of performing meta-analysis of these studies that might help to inform future study design. The use of random-effects models as we did in this study can help to account for expecting heterogeneity [35]; however, appropriate caution still needs to be counted when interpreting the results. A fourth limitation is that the small sample size in animal studies which may exaggerate biases like publication bias which can affect the reliability and the validity of the study outcome. To overcome this limitation, one of the approaches we used was to calculate the standardized mean differences (SMD) through Hedges’ g effect sizes which is based on Cohen’s D but includes a correction factor for small sample size bias [35]. More specific limitation for this study was regarding the inclusion criteria of the systematic review, and we did not include animal models of disease or pathology, as we want to estimate the effect on healthy animal model, but this can be investigated and tested in future systematic review and meta-analysis with specific research question on the effect of NSAID administration on bone healing in diseased animal models. Other limitation is the meta-analysis of previous experimental studies did not include healing outcomes other than bone fracture (e.g., pain behavior or levels of inflammatory mediators). The measurements of pain behavior are subjective and depend on the interpretation of the examiner, which may lead to high heterogeneity among studies.

Clinical implications

Drug pharmacokinetics vary between species, and this must be considered when extrapolating data from animals to humans. Therefore, it is important to investigate NSAID doses that are equivalent to those used in humans after bone fracture (through interspecies allometric scaling for dose conversion from animal to human studies) and to evaluate effects using different types of animal and bone fracture models. The increasing amount of evidence from animal studies and the results from this systematic review and meta-analysis indicate that caution should be exercised when using NSAID after bone fractures or with specific orthopedic surgical procedures until prospective human clinical studies indicate otherwise.

Conclusions

Our findings provide some guidance for future laboratory and clinical research. First, it is important to test different hypotheses of bone fracture healing in small animals, and mice especially because mice models provide opportunities to examine genetics and create knock-out species. Our results also indicate the need for studies that compare the effect of NSAID administration on bone healing outcomes between male and female animal models. Overall, choosing appropriate animal model to test the effect of NSAID on fracture bone healing should take into consideration the animals’ species, type of the bone, and age of the animal.

Moreover, our results demonstrate it is important to choose the suitable time of sample collection based on what healing outcome to be measured. Histomorphometric outcome measurements require more than 21 days to show results that are comparable to the control group. Second, improvements in internal (study quality) and external (publication bias) validity might provide more information for the translation of the data to clinical trials, and more robust exploration of the efficacy limits in such studies could inform inclusion and exclusion criteria for these trials.

It is increasingly clear that the function of COX and their products are critical for bone healing. Most animal and human studies support the conclusion that NSAID administration can delay or impair bone fracture healing. However, the anti-inflammatory and analgesic effect of NSAID seems to be beneficial in treating post-traumatic pain and edema. The need for a new approach to using NSAID that preserves its pain control properties without affecting bone healing is justifiable and important.

Availability of data and materials

All data generated or analyzed during this study are included in this published article and its Supplementary files 1 and 2.

Abbreviations

NSAID:

Non-steroidal anti-inflammatory drugs

μ-CT:

Micro-computed tomography

SYRCLE:

SYstematic Review Center for Laboratory animal Experimentation

COX:

Cyclooxygenase

SD:

Standard deviation

SE:

Standard error

SMD:

Standardized mean differences

MF:

Maximum force

References

  1. Bergese S, Castellon-Larios K. The effectiveness of a single dose of oral ibuprofen plus caffeine in acute postoperative pain in adults. Evid-Based Med. 2016;21(1):24. https://0-doi-org.brum.beds.ac.uk/10.1136/ebmed-2015-110278.

    Article  PubMed  Google Scholar 

  2. Bhattacharyya T, Levin R, Vrahas MS, Solomon DH. Nonsteroidal antiinflammatory drugs and nonunion of humeral shaft fractures. Arthritis Rheum. 2005;53(3):364–7. https://0-doi-org.brum.beds.ac.uk/10.1002/art.21170.

    Article  CAS  PubMed  Google Scholar 

  3. Moore RA, Derry S, McQuay HJ, Wiffen PJ. Single dose oral analgesics for acute postoperative pain in adults. Cochrane Database Syst Rev. 2011;(9):Cd008659.

  4. Pozzi A, Gallelli L. Pain management for dentists: the role of ibuprofen. Annali di stomatologia. 2011;2(3-4 Suppl):3–24.

    PubMed  Google Scholar 

  5. Vuolteenaho K, Moilanen T, Moilanen E. Non-steroidal anti-inflammatory drugs, cyclooxygenase-2 and the bone healing process. Basic Clin Pharmacol Toxicol. 2008;102(1):10–4. https://0-doi-org.brum.beds.ac.uk/10.1111/j.1742-7843.2007.00149.x.

    Article  CAS  PubMed  Google Scholar 

  6. Altman RD, Latta LL, Keer R, Renfree K, Hornicek FJ, Banovac K. Effect of nonsteroidal antiinflammatory drugs on fracture healing: a laboratory study in rats. Journal of Orthopaedic Trauma. 1995;9(5):392–400. https://0-doi-org.brum.beds.ac.uk/10.1097/00005131-199505000-00006.

    Article  CAS  PubMed  Google Scholar 

  7. Beck A, Krischak G, Sorg T, Augat P, Farker K, Merkel U, et al. Influence of diclofenac (group of nonsteroidal anti-inflammatory drugs) on fracture healing. Arch Orthop Trauma Surg. 2003;123(7):327–32. https://0-doi-org.brum.beds.ac.uk/10.1007/s00402-003-0537-5.

    Article  CAS  PubMed  Google Scholar 

  8. Bo J, Sudmann E, Marton PF. Effect of indomethacin on fracture healing in rats. Acta Orthop Scand. 1976;47(6):588–99. https://0-doi-org.brum.beds.ac.uk/10.3109/17453677608988744.

    Article  CAS  PubMed  Google Scholar 

  9. Cappello T, Nuelle JAV, Katsantonis N, Nauer RK, Lauing KL, Jagodzinski JE, et al. Ketorolac administration does not delay early fracture healing in a juvenile rat model: a pilot study. J Pediatric Orthop. 2013;33(4):415–21. https://0-doi-org.brum.beds.ac.uk/10.1097/BPO.0b013e318288b46f.

    Article  Google Scholar 

  10. Gerstenfeld LC, Al-Ghawas M, Alkhiary YM, Cullinane DM, Krall EA, Fitch JL, et al. Selective and nonselective cyclooxygenase-2 inhibitors and experimental fracture-healing. Reversibility of effects after short-term treatment. J Bone Joint Surg Am Vol. 2007;89(1):114–25. https://0-doi-org.brum.beds.ac.uk/10.2106/00004623-200701000-00016.

    Article  CAS  Google Scholar 

  11. Sandberg O, Aspenberg P. Different effects of indomethacin on healing of shaft and metaphyseal fractures. Acta Orthopaedica. 2015;86(2):243–7. https://0-doi-org.brum.beds.ac.uk/10.3109/17453674.2014.973328.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Simon AM, O'Connor JP, Simon AM, O'Connor JP. Dose and time-dependent effects of cyclooxygenase-2 inhibition on fracture-healing. J Bone Joint Surg Am Vol. 2007;89(3):500–11. https://0-doi-org.brum.beds.ac.uk/10.2106/JBJS.F.00127.

    Article  Google Scholar 

  13. Zura R, Xiong Z, Einhorn T, Watson JT, Ostrum RF, Prayson MJ, et al. Epidemiology of fracture nonunion in 18 human bones. JAMA Surgery. 2016;151(11):e162775. https://0-doi-org.brum.beds.ac.uk/10.1001/jamasurg.2016.2775.

    Article  PubMed  Google Scholar 

  14. Davis TR, Ackroyd CE. Non-steroidal anti-inflammatory agents in the management of Colles' fractures. Brit J Clin Pract. 1988;42(5):184–9.

    CAS  PubMed  Google Scholar 

  15. Adolphson P, Abbaszadegan H, Jonsson U, Dalen N, Sjoberg HE, Kalen S. No effects of piroxicam on osteopenia and recovery after Colles' fracture. A randomized, double-blind, placebo-controlled, prospective trial. Arch Orthop Trauma Surg. 1993;112(3):127–30. https://0-doi-org.brum.beds.ac.uk/10.1007/BF00449987.

    Article  CAS  PubMed  Google Scholar 

  16. Sagi HC, Jordan CJ, Barei DP, Serrano-Riera R, Steverson B. Indomethacin prophylaxis for heterotopic ossification after acetabular fracture surgery increases the risk for nonunion of the posterior wall. J Orthop Trauma. 2014;28(7):377–83. https://0-doi-org.brum.beds.ac.uk/10.1097/BOT.0000000000000049.

    Article  PubMed  Google Scholar 

  17. Marquez-Lara A, Hutchinson ID, Nunez F, Jr., Smith TL, Miller AN. Nonsteroidal anti-inflammatory drugs and bone-healing: a systematic review of research quality. JBJS Rev. 2016;4(3):1–14.

  18. Dimmen S, Nordsletten L, Madsen JE. Parecoxib and indomethacin delay early fracture healing: a study in rats. Clin Orthop. 2009;467(8):1992–9. https://0-doi-org.brum.beds.ac.uk/10.1007/s11999-009-0783-0.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Endo K, Sairyo K, Komatsubara S, Sasa T, Egawa H, Ogawa T, et al. Cyclooxygenase-2 inhibitor delays fracture healing in rats. Acta Orthopaedica. 2005;76(4):470–4. https://0-doi-org.brum.beds.ac.uk/10.1080/17453670510041439.

    Article  PubMed  Google Scholar 

  20. Inal S, Kabay S, Cayci MK, Kuru HI, Altikat S, Akkas G, et al. Comparison of the effects of dexketoprofen trometamol, meloxicam and diclofenac sodium on fibular fracture healing, kidney and liver: an experimental rat model. Injury. 2014;45(3):494–500. https://0-doi-org.brum.beds.ac.uk/10.1016/j.injury.2013.10.002.

    Article  PubMed  Google Scholar 

  21. Chappard D, Basle MF, Aguado E, Goyenvalle E, Nyangoga H. A non-steroidal anti-inflammatory drug (ketoprofen) does not delay beta-TCP bone graft healing. Acta Biomaterialia. 2010;6(8):3310–7.

    Article  Google Scholar 

  22. Donovan M, George M, Petersen E, Fredericks D, Femino JE, Lack WD, et al. Effect of aspirin on bone healing in a rabbit ulnar osteotomy model. J Bone Joint Surg – Am Vol. 2013;95(6):488–96.

    Article  Google Scholar 

  23. Goodman SB, Ma T, Miyanishi K, Smith RL, Oh K, Wadsworth S, et al. Effects of a p38 MAP kinase inhibitor on bone ingrowth and tissue differentiation in rabbit chambers. J Biomed Materials Res Part A. 2007;81(2):310–6.

    Article  CAS  Google Scholar 

  24. Singh A, Shekhar S, Saraf SK, Garbyal RS. Effect of etoricoxib on fracture healing - an experimental study. Biomed Res. 2011;22(1):52–6.

    CAS  Google Scholar 

  25. Morgan EF, De Giacomo A, Gerstenfeld LC. Overview of skeletal repair (fracture healing and its assessment). Methods Mol Biol. 2014;1130:13–31. https://0-doi-org.brum.beds.ac.uk/10.1007/978-1-62703-989-5_2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Histing T, Garcia P, Holstein JH, Klein M, Matthys R, Nuetzi R, et al. Small animal bone healing models: standards, tips, and pitfalls results of a consensus meeting. Bone. 2011;49(4):591–9. https://0-doi-org.brum.beds.ac.uk/10.1016/j.bone.2011.07.007.

    Article  CAS  PubMed  Google Scholar 

  27. Bissinger O, Kreutzer K, Gotz C, Hapfelmeier A, Pautke C, Vogt S, et al. A biomechanical, micro-computertomographic and histological analysis of the influence of diclofenac and prednisolone on fracture healing in vivo. BMC Musculoskelet Disord. 2016;17(1):383. https://0-doi-org.brum.beds.ac.uk/10.1186/s12891-016-1241-2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Foux A, Black RC, Uhthoff HK. Quantitative measures for fracture healing: an in-vitro biomechanical study. J Biomech Eng. 1990;112(4):401–6. https://0-doi-org.brum.beds.ac.uk/10.1115/1.2891203.

    Article  CAS  PubMed  Google Scholar 

  29. O'Connor AM, Sargeant JM. Critical appraisal of studies using laboratory animal models. ILAR J. 2014;55(3):405–17. https://0-doi-org.brum.beds.ac.uk/10.1093/ilar/ilu038.

    Article  CAS  PubMed  Google Scholar 

  30. Gerstenfeld LC, Wronski TJ, Hollinger JO, Einhorn TA. Application of histomorphometric methods to the study of bone repair. J Bone Mineral Res. 2005;20(10):1715–22. https://0-doi-org.brum.beds.ac.uk/10.1359/JBMR.050702.

    Article  Google Scholar 

  31. Hacker SA, Healey RM, Yoshioka M, Coutts RD. A methodology for the quantitative assessment of articular cartilage histomorphometry. Osteoarthritis Cartilage. 1997;5(5):343–55. https://0-doi-org.brum.beds.ac.uk/10.1016/S1063-4584(97)80038-6.

    Article  CAS  PubMed  Google Scholar 

  32. Al-Waeli H, Ana Paula Reboucas MM, Mansour A, Tamimi F, Nicolau B. Non steroidal anti inflammatory drugs and bone healing in animal models systematic review and meta analysis; 2017.

    Google Scholar 

  33. de Vries RBM, Hooijmans CR, Langendam MW, van Luijk J, Leenaars M, Ritskes-Hoitinga M, et al. A protocol format for the preparation, registration and publication of systematic reviews of animal intervention studies. Evid-Based Preclinical Med. 2015;2(1):1–9.

    Article  Google Scholar 

  34. Ouzzani M, Hammady H, Fedorowicz Z, Elmagarmid A. Rayyan-a web and mobile app for systematic reviews. Systematic reviews. 2016;5(1):210. https://0-doi-org.brum.beds.ac.uk/10.1186/s13643-016-0384-4.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Hooijmans CR, Rovers MM, de Vries RB, Leenaars M, Ritskes-Hoitinga M, Langendam MW. SYRCLE's risk of bias tool for animal studies. BMC Med Res Methodol. 2014;14(1):43. https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2288-14-43.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Kilkenny C, Parsons N, Kadyszewski E, Festing MF, Cuthill IC, Fry D, et al. Survey of the quality of experimental design, statistical analysis and reporting of research using animals. PLoS One. 2009;4(11):e7824. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pone.0007824.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Chalmers I, Hedges LV, Cooper H. A brief history of research synthesis. Eval Health Prof. 2002;25(1):12–37. https://0-doi-org.brum.beds.ac.uk/10.1177/0163278702025001003.

    Article  PubMed  Google Scholar 

  38. Vesterinen HM, Sena ES, Egan KJ, Hirst TC, Churolov L, Currie GL, et al. Meta-analysis of data from animal studies: a practical guide. J Neurosci Methods. 2014;221:92–102. https://0-doi-org.brum.beds.ac.uk/10.1016/j.jneumeth.2013.09.010.

    Article  CAS  PubMed  Google Scholar 

  39. Hedges LV, Olkin I. Statistical methods for meta-analysis. Academic Press; 2014.

  40. Hooijmans CR, IntHout J, Ritskes-Hoitinga M, Rovers MM. Meta-analyses of animal studies: an introduction of a valuable instrument to further improve healthcare. ILAR J. 2014;55(3):418–26. https://0-doi-org.brum.beds.ac.uk/10.1093/ilar/ilu042.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Petitti DB. Approaches to heterogeneity in meta-analysis. Stat Med. 2001;20(23):3625–33. https://0-doi-org.brum.beds.ac.uk/10.1002/sim.1091.

    Article  CAS  PubMed  Google Scholar 

  42. Pierce CA. Software Review: Borenstein, M., Hedges, L. V., Higgins, J. P. T., & Rothstein, H. R. (2006). Comprehensive meta-analysis (Version 2.2.027) [Computer software]. Englewood, NJ: Biostat. Organ Res Methods. 2008;11(1):188–91. https://0-doi-org.brum.beds.ac.uk/10.1177/1094428106296641.

    Article  Google Scholar 

  43. Akman Ş, Gögüs A, Sener N, Bilgiç B, Aksoy B, Seckin F. Effect of diclofenac sodium on union of tibial fractures in rats. Adv Therapy. 2002;19(3):119–25. https://0-doi-org.brum.beds.ac.uk/10.1007/BF02850267.

    Article  CAS  Google Scholar 

  44. Allen HL, Wase A, Bear WT. Indomethacin and aspirin: effect of nonsteroidal anti-inflammatory agents on the rate of fracture repair in the rat. Acta Orthop Scand. 1980;51(4):595–600. https://0-doi-org.brum.beds.ac.uk/10.3109/17453678008990848.

    Article  CAS  PubMed  Google Scholar 

  45. Bergenstock M, Min W, Simon AM, Sabatino C, O'Connor JP. A comparison between the effects of acetaminophen and celecoxib on bone fracture healing rats. J Orthop Trauma. 2005;19(10):717–23. https://0-doi-org.brum.beds.ac.uk/10.1097/01.bot.0000184144.98071.5d.

    Article  PubMed  Google Scholar 

  46. Brown KM, Saunders MM, Kirsch T, Donahue HJ, Reid JS. Effect of COX-2-specific inhibition on fracture-healing in the rat femur. J Bone Joint Surg - Am Vol. 2004;86(1):116–23. https://0-doi-org.brum.beds.ac.uk/10.2106/00004623-200401000-00017.

    Article  Google Scholar 

  47. Giordano V, Giordano M, Knackfuss IG, Apfel MIR, Gomes RDC. Effect of tenoxicam on fracture healing in rat tibiae. Injury. 2003;34(2):85–94. https://0-doi-org.brum.beds.ac.uk/10.1016/S0020-1383(02)00199-7.

    Article  PubMed  Google Scholar 

  48. Krischak GD, Augat P, Blakytny R, Claes L, Kinzl L, Beck A. The non-steroidal anti-inflammatory drug diclofenac reduces appearance of osteoblasts in bone defect healing in rats. Arch Orthop Trauma Surg. 2007;127(6):453–8. https://0-doi-org.brum.beds.ac.uk/10.1007/s00402-007-0288-9.

    Article  CAS  PubMed  Google Scholar 

  49. Li K-H, Cheng L, Zhu Y, Deng G-B, Long H-T. Effects of a selective cyclooxygenase-2 inhibitor (celecoxib) on fracture healing in rats. Indian J Orthop. 2013;47(4):395–401. https://0-doi-org.brum.beds.ac.uk/10.4103/0019-5413.114930.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Sudmann E, Bang G. Indomethacin-induced inhibition of haversian remodelling in rabbits. Acta Orthop Scand. 1979;50(6 Pt 1):621–7. https://0-doi-org.brum.beds.ac.uk/10.3109/17453677908991283.

    Article  CAS  PubMed  Google Scholar 

  51. Utvag SE, Fuskevag OM, Shegarfi H, Reikeras O. Short-term treatment with COX-2 inhibitors does not impair fracture healing. J Invest Surg. 2010;23(5):257–61. https://0-doi-org.brum.beds.ac.uk/10.3109/08941939.2010.481009.

    Article  PubMed  Google Scholar 

  52. Tiseo BC, Namur GN, De Paula EJL, Mattar R Jr, De Oliveira CRGCM. Experimental study of the action of COX-2 selective nonsteroidal anti-inflammatory drugs and traditional anti-inflammatory drugs in bone regeneration. Clinics. 2006;61(3):223–30. https://0-doi-org.brum.beds.ac.uk/10.1590/S1807-59322006000300007.

    Article  PubMed  Google Scholar 

  53. Simon AM, Manigrasso MB, O'Connor JP. Cyclo-oxygenase 2 function is essential for bone fracture healing. J Bone Mineral Res. 2002;17(6):963–76. https://0-doi-org.brum.beds.ac.uk/10.1359/jbmr.2002.17.6.963.

    Article  CAS  Google Scholar 

  54. Sassioto MC, Inouye CM, Aydos RD, Figueiredo AS. Bone repair in rats treated with sodic diclofenac and calcitonin. Acta Cirurgica Brasileira. 2006;21(suppl 4):40–4. https://0-doi-org.brum.beds.ac.uk/10.1590/S0102-86502006001000010.

    Article  PubMed  Google Scholar 

  55. Reikeraas O, Engebretsen L. Effects of ketoralac tromethamine and indomethacin on primary and secondary bone healing. An experimental study in rats. Arch Orthop Trauma Surg. 1998;118(1):50–2. https://0-doi-org.brum.beds.ac.uk/10.1007/s004020050310.

    Article  CAS  PubMed  Google Scholar 

  56. Mullis BH, Copland ST, Weinhold PS, Miclau T, Lester GE, Bos GD. Effect of COX-2 inhibitors and non-steroidal anti-inflammatory drugs on a mouse fracture model. Injury. 2006;37(9):827–37. https://0-doi-org.brum.beds.ac.uk/10.1016/j.injury.2005.12.018.

    Article  PubMed  Google Scholar 

  57. Matsumoto MA, De Oliveira A, Ribeiro Junior PD, Nary Filho H, Ribeiro DA. Short-term administration of non-selective and selective COX-2 NSAIDs do not interfere with bone repair in rats. J Mol Histol. 2008;39(4):381–7. https://0-doi-org.brum.beds.ac.uk/10.1007/s10735-008-9176-z.

    Article  CAS  PubMed  Google Scholar 

  58. Gerstenfeld LC, Thiede M, Siebert K, Mielke C, Phippard D, Svagr B, et al. Differential inhibition of fracture healing by non-selective and cyclooxygenase-2 selective non-steroidal anti-inflammatory drugs. J Orthop Res. 2003;21(4):670–5. https://0-doi-org.brum.beds.ac.uk/10.1016/S0736-0266(03)00003-2.

    Article  CAS  PubMed  Google Scholar 

  59. Dimmen S, Nordsletten L, Engebretsen L, Steen H, Madsen JE. Negative effect of parecoxib on bone mineral during fracture healing in rats. Acta Orthop. 2008;79(3):438–44. https://0-doi-org.brum.beds.ac.uk/10.1080/17453670710015373.

    Article  PubMed  Google Scholar 

  60. Lin HN, Cottrell J, O'Connor JP. Variation in lipid mediator and cytokine levels during mouse femur fracture healing. J Orthopaedic Res. 2016;34(11):1883–93. https://0-doi-org.brum.beds.ac.uk/10.1002/jor.23213.

    Article  CAS  Google Scholar 

  61. O'Connor JP, Manigrasso MB, Kim BD, Subramanian S. Fracture healing and lipid mediators. BoneKEy Reports. 2014;3:517.

    PubMed  PubMed Central  Google Scholar 

  62. Su B, O'Connor JP. NSAID therapy effects on healing of bone, tendon, and the enthesis. J Appl Physiol. 2013;115(6):892–9. https://0-doi-org.brum.beds.ac.uk/10.1152/japplphysiol.00053.2013.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. O'Connor JP, Simon AM, Manigrasso MB. Cyclooxygenase-2 function is essential for bone fracture healing. Inflammopharmacol. 2003;11(4):541–3.

    Google Scholar 

  64. Spiro AS, Beil FT, Baranowsky A, Barvencik F, Schilling AF, Nguyen K, et al. BMP-7-induced ectopic bone formation and fracture healing is impaired by systemic NSAID application in C57BL/6-mice. J Orthop Res. 2010;28(6):785–91. https://0-doi-org.brum.beds.ac.uk/10.1002/jor.21044.

    Article  CAS  PubMed  Google Scholar 

  65. Sevimli R, Üzel M, Sayar H, Kalender AM, Dökmeci O. The effect of dexketoprofen trometamol on the healing of diaphysis fractures of rat tibia. Acta Orthopaedica et Traumatologica Turcica. 2013;47(6):423–9. https://0-doi-org.brum.beds.ac.uk/10.3944/AOTT.2013.3093.

    Article  PubMed  Google Scholar 

  66. Dorey FJ, Kody MH, Kabo JM, More RC, Meals RA. The effects of two nonsteroidal antiinflammatory drugs on limb swelling, joint stiffness, and bone torsional strength following fracture in a rabbit model. Clin Orthop. 1989;0(247):306-312.

  67. Huo MH, Troiano NW, Pelker RR, Gundberg CM, Friedlaender GE. The influence of ibuprofen on fracture repair: biomechanical, biochemical, histologic, and histomorphometric parameters in rats. J Orthop Res. 1991;9(3):383–90. https://0-doi-org.brum.beds.ac.uk/10.1002/jor.1100090310.

    Article  CAS  PubMed  Google Scholar 

  68. Huss MK, Felt SA, Pacharinsak C. Influence of pain and analgesia on orthopedic and wound-healing models in rats and mice. Comp Med. 2019;69(6):535–45. https://0-doi-org.brum.beds.ac.uk/10.30802/AALAS-CM-19-000013.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Borgeat A, Ofner C, Saporito A, Farshad M, Aguirre J. The effect of nonsteroidal anti-inflammatory drugs on bone healing in humans: A qualitative, systematic review. J Clin Anesth. 2018;49:92–100. https://0-doi-org.brum.beds.ac.uk/10.1016/j.jclinane.2018.06.020.

    Article  CAS  PubMed  Google Scholar 

  70. Haffner-Luntzer M, Kovtun A, Rapp AE, Ignatius A. Mouse Models in Bone Fracture Healing Research. Curr Mol Biol Reports. 2016;2(2):101–11. https://0-doi-org.brum.beds.ac.uk/10.1007/s40610-016-0037-3.

    Article  Google Scholar 

  71. Peric M, Dumic-Cule I, Grcevic D, Matijasic M, Verbanac D, Paul R, et al. The rational use of animal models in the evaluation of novel bone regenerative therapies. Bone. 2015;70:73–86. https://0-doi-org.brum.beds.ac.uk/10.1016/j.bone.2014.07.010.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Thank you to Dr. Genevieve Castonguay for her help in the editing of the manuscript writing.

Funding

No applicable

Author information

Authors and Affiliations

Authors

Contributions

HA conceived the topic of the manuscript and object of the analyses, performed all steps of the systematic review (searching the literature, selection of studies, data extraction, and quality assessment), statistical analyses, and wrote the manuscript. AR participated in the search process, screening, and selection of eligible studies. AM participated in the data extraction and the quality assessment of the included studies. MM designed the search strategies of eight databases and participated in the writing of the methodology. FT supervised the conception of the manuscript, and all steps of the systematic review and meta-analyses including the statistical analyses and writing the manuscript. BN supervised the conception of the manuscript, and all steps of the systematic review and meta-analyses including the statistical analyses and writing the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Haider Al-Waeli.

Ethics declarations

Ethics approval and consent to participate

Not applicable

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Systematic review and meta-analysis registered protocol in SYRCLE website: https://www.radboudumc.nl/en/research/technology-centers/animal-research-facility/systematic-review-center-for-laboratory-animal-experimentation/protocols NSAIDs effect on bone healing after bone fracture in animal models. Table S2. Search strategy. Table S3. SYRCLE’s tool for assessing the risk of bias (Hooijman et al 2014) (1). Table S4. Study characteristics. Table S5. Subgroup analysis by species regarding the effect of NSAID administration vs control on maximum force to fracture outcome. Table S6. Subgroup analysis by sex regarding the effect of NSAID administration vs control on maximum force to fracture outcome. Table S7. Subgroup analysis by age regarding the effect of NSAID administration vs control on maximum force to fracture outcome (1= <8 wks, 2=8-16wks, >16wks, 4=not mentioned). Table S8. Subgroup analysis by type pf NSAID regarding the effect of NSAID administration vs control on maximum force to fracture outcome. Table S9. Subgroup analysis by time point regarding the effect of NSAID administration vs control on maximum force to fracture outcome (1=<21days, 2=21-48days, 3=>48days). Table S10. Subgroup analysis by bone fracture site regarding the effect of NSAID administration vs control on maximum force to fracture outcome. Table S11. Subgroup analysis by species regarding the effect of NSAID administration vs control on stiffness to fracture outcome. Table S12. Subgroup analysis by sex regarding the effect of NSAID administration vs control on stiffness to fracture outcome. Table S13. Subgroup analysis by age regarding the effect of NSAID administration vs control on stiffness to fracture outcome (1= <8 wks, 2=8-16wks, >16wks, 4=not mentioned). Table S14. Subgroup analysis by type of NSAID regarding the effect of NSAID administration vs control on stiffness to fracture outcome. Table S15. Subgroup analysis by type of time point regarding the effect of NSAID administration vs control on stiffness to fracture outcome (1=<21days, 2=21-48days, 3=>48days). Table S16. Subgroup analysis by bone fracture site regarding the effect of NSAID administration vs control on stiffness to fracture outcome. Table S17. Subgroup analysis by species regarding the effect of NSAID administration vs control on work to failure outcome. Table S18. Subgroup analysis by sex regarding the effect of NSAID administration vs control on work to failure outcome. Table S19. Subgroup analysis by age regarding the effect of NSAID administration vs control on work to failure outcome (1= <8 wks, 2=8-16wks, >16wks, 4=not mentioned). Table S20. Subgroup analysis by type of NSAID regarding the effect of NSAID administration vs control on work to failure outcome. Table S21. Subgroup analysis by time point regarding the effect of NSAID administration vs control on work to failure outcome (1=<21days, 2=21-48days, 3=>48days). Table S22. Subgroup analysis by bone fracture site the effect of NSAID administration vs control on work to failure outcome. Table S23. Subgroup analysis by species the effect of NSAID administration vs control on histomorphometric outcome. Table S24. Subgroup analysis by sex the effect of NSAID administration vs control on histomorphometric outcome. Table S25. Subgroup analysis by age the effect of NSAID administration vs control on histomorphometric outcome (1= <8 wks, 2=8-16wks, >16wks, 4=not mentioned). Table S26. Subgroup analysis by type of NSAID the effect of NSAID administration vs control on histomorphometric outcome. Table S27. Subgroup analysis by time point the effect of NSAID administration vs control on histomorphometric outcome (1=<21days, 2=21-48days, 3=>48days). Table S28. Subgroup analysis by bone fracture site the effect of NSAID administration vs control on histomorphometric outcome.

Additional file 2.

Additional file 3.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Al-Waeli, H., Reboucas, A.P., Mansour, A. et al. Non-steroidal anti-inflammatory drugs and bone healing in animal models—a systematic review and meta-analysis. Syst Rev 10, 201 (2021). https://0-doi-org.brum.beds.ac.uk/10.1186/s13643-021-01690-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s13643-021-01690-w

Keywords