Analgesic medicines for adults with low back pain: protocol for a systematic review and network meta-analysis
Systematic Reviews volume 9, Article number: 255 (2020)
There is limited evidence for the comparative effectiveness of analgesic medicines for adults with low back pain. This systematic review and network meta-analysis aims to determine the analgesic effect, safety, acceptability, effect on function, and relative rank according to analgesic effect, safety, acceptability, and effect on function of a single course of [an] analgesic medicine(s) or combination of these medicines for people with low back pain.
We will include published and unpublished randomised trials written in any language that compare an analgesic medicine to either another medicine, placebo/sham, or no intervention in adults with low back pain, grouped according to pain duration: acute (fewer than 6 weeks), sub-acute (6 to 12 weeks), and chronic (greater than 12 weeks). The co-primary outcomes are pain intensity following treatment and safety (adverse events). The secondary outcomes are function and acceptability (all-cause dropouts). We will perform a network meta-analysis to compare and rank analgesic medicines. We will form judgements of confidence in the results using the Confidence in Network Meta-Analysis (CINeMA) methodology.
This network meta-analysis will establish which medicine, or combination of medicines, is most effective for reducing pain and safest for adults with low back pain.
Systematic review registration
Low back pain (LBP) is a global health problem; it is the leading cause of disability in 126 out of 195 countries according to the Global Burden of Disease Study 2017 . Non-specific LBP represents over 90% of cases, in which a pathoanatomical source of pain cannot be reliably identified [2, 3], and is classified by the duration of pain: acute (fewer than 6 weeks), sub-acute (6 to 12 weeks), and chronic (greater than 12 weeks) [4, 5]. The estimated lifetime prevalence of LBP is up to 80%, meaning many adults will experience an episode of LBP at least once [6, 7]. Following onset, most people experience a reduction in pain within the first 6 weeks . However, only 60% are fully recovered at 12 weeks and 70% will experience a recurrence of LBP within 12 months of recovery [9, 10]. Approximately 20% of the global population live with chronic LBP at any time .
Despite recommendations for non-pharmacological interventions to play a greater role in the management of LBP , analgesic medicines are the most commonly prescribed intervention across a range of primary care settings (e.g. general practice, emergency department) [12,13,14,15,16,17,18]. Analgesic medicines include non-steroidal anti-inflammatory drugs (NSAIDs), opioids, paracetamol (acetaminophen), anti-convulsants, anti-depressants, muscle relaxants, and corticosteroids. These medicines are used in a range of countries [13,14,15,16, 19, 20], but there appear to be no common patterns to prescribing . For example, muscle relaxant medicines are commonly prescribed in the USA and Italy [14, 15, 19], but seldom in Australia .
People with LBP (and clinicians) want to know the most effective medicine for their condition [21,22,23,24]. This requires information about comparative effectiveness—the effect of a medicine compared to other medicines—to be available for clinical decision making. Comparative effectiveness has been insufficiently described in syntheses of the literature to date [25,26,27,28,29,30,31,32,33,34,35]. This is understandable, as most of these systematic reviews [28, 29, 31, 32] investigated a single comparison, usually efficacy, the effect of a medicine compared to sham. Several of these reviews [25, 35] additionally examined a limited number of effectiveness comparisons. No quantitative synthesis was made of these data, which was appropriate because methods for single comparisons should not be used across multiple comparisons. Comparative effectiveness data that are limited and lack synthesis have low utility for clinical decision-making.
Network meta-analysis (NMA) can be used to synthesise data across multiple comparisons. NMA provides valid estimates of comparative effectiveness by fitting a single statistical model to a connected network of interventions when there is confidence in the assumptions of transitivity and coherence (see the ‘Assumption of transitivity’ and ‘Assessment of network heterogeneity and coherence’ sections) [36,37,38,39]. The results are relative effect estimates for each comparison between medicines of interest. These data are applicable to decision scenarios if adults with LBP, clinicians, and policy makers are determining which of several medicines should be used. The relative effect estimates can be used to rank interventions based on their effect on relevant outcomes, which may also assist clinical decision-making. Therefore, this NMA will evaluate the comparative effectiveness of analgesic medicines for adults with LBP.
The objectives of this review are to:
Determine the effect on pain intensity, safety, acceptability, and effect on function of a single course of [an] analgesic medicine(s) or combination of these medicines.
Determine a relative rank for each intervention according to its effect on pain intensity, safety, acceptability, and effect on function.
This protocol follows the Preferred Reporting Items for Systematic Reviews and Meta-Analysis Protocols (PRISMA-P) guidelines [40, 41], provided as Additional file 1. The review was registered with the International Prospective Register of Systematic Reviews on 31 October 2019 (PROSPERO ID: CRD42019145257).
Eligibility criteria for inclusion in the review
We will include randomised trials (RCTs) that provide at least one comparison between two interventions of interest (see the ‘Interventions’ section). There are no restrictions on language or publication status—we will include unpublished data because it may meaningfully alter the relative effectiveness of a medicine [42, 43]. We will include data from parallel group trials. We will also include data from the first phase of crossover RCTs because of the lack of relative stability in recent-onset LBP and known carry over effects in different classes of analgesic medicines (e.g. anti-depressants) [44, 45]. This approach has been used in published NMAs [36, 46, 47]. We will exclude cluster RCTs.
We will exclude enriched-enrolment RCTs. Although the single-arm run-in phase of enriched-enrolment trials might provide information on adverse effects, the NMA will not include data from uncontrolled trials. Additionally, the double-blind/randomised phase of enriched-enrolment RCTs has questionable external validity  and is limited in detecting adverse effects [48, 49].
We will include studies that randomised adults with non-specific LBP, defined as a primary area of pain between the twelfth rib and gluteal fold, with or without associated leg pain [4, 5]. We will consider three different pain durations, which we will analyse separately: acute (fewer than 6 weeks), sub-acute (6 to 12 weeks), and chronic (greater than 12 weeks) . We will include studies that randomised participants with heterogeneous pain conditions if separate data is obtainable for the participants with LBP. Participants may be experienced or naïve to the trial intervention, which we will assess for the evaluation of transitivity (see the ‘Assumption of transitivity’ section). We will exclude interventional (surgical/operative) settings and studies where greater than 20% of participants have leg pain that meets the definition of sciatica according to Koes et al. , where patients have unilateral leg pain greater than LBP radiating to the foot or toes with associated neurological indications, or where LBP is attributable to a specific pathology, such as infection, neoplasm, metastasis, inflammatory disease, or fractures.
We will include medicines from the following classes: NSAIDs, paracetamol (acetaminophen), opioid analgesics, anti-convulsants, anti-depressants, muscle relaxant medicines, and corticosteroids. These medicines must be listed on the World Health Organization (WHO) Anatomical Therapeutic Chemical (ATC) system and licensed for current use by at least one of the following agencies: US Food and Drug Administration (FDA), the Australian Therapeutic Goods Administration (TGA), the European Medicines Agency (EMA), or the Medicine and Healthcare Products Regulatory Agency (MHRA). These interventions of interest are listed in Additional file 2 with their ATC codes and licensing status. We will include additional analgesic medicines from these classes, identified during the review, provided they are licensed by one of the above agencies. Medicines may be delivered as mono or combination therapy via any systemic route of administration (e.g. oral, intravenous, intramuscular, buccal, sublingual). We will exclude non-systemic administration (such as topical, intraarticular, or epidural administration). We will not exclude studies that assign non-pharmacological co-interventions to one or more of the intervention arms. We will consider these studies in the evaluation of transitivity .
We will represent each intervention of interest as a separate node in a network of all possible comparisons between interventions (not shown due to visual complexity). We define the placebo/sham intervention node as any drug intervention that does not contain an active analgesic ingredient (including ‘active’ placebo arms). We consider that no-treatment includes continuation of usual care or being placed on a waitlist. We will further define nodes according to route of administration, which means that a single drug may be represented by multiple nodes in the network (e.g. oral diclofenac and intravenous diclofenac are separate nodes).
We will classify the dose of each medicine as either (i) standard dosing range (SDR), (ii) below the SDR, or (iii) above the SDR. We will source the relevant SDR for each medicine using the following hierarchy: Prescriber’s Digital Reference , MIMS , or the Australian Medicines Handbook . If the SDR is not provided by any of these sources, we will use the medicine’s licensed dosing range (LDR). Typically, the LDR is equivalent to the SDR used in clinical practice, although the SDR may be lower than the LDR. We will identify a medicine’s LDR using the following hierarchy: FDA, MHRA, EMA, or TGA. Where a study contains two or more intervention arms within the same dosing range, we will combine these arms using formulae in the Cochrane Handbook for Systematic Reviews of Interventions .
The primary outcomes for this review are pain and safety.
Pain is defined as pain intensity, measured at the time point closest to the end of treatment. Pain intensity may be measured with a continuous self-report scale (e.g. visual analogue scale (VAS) or numeric rating scale (NRS)), a rating scale within a composite measure of pain (e.g. McGill Pain Questionnaire), or an ordinal scale (we will consider such ordinal scales to exhibit continuous properties). We will not exclude studies that use other measurement tools.
Safety is the number of participants who experience an adverse effect during the treatment period. We define adverse effects according to the US FDA, as ‘any untoward medical occurrence associated with the use of a drug in humans, whether or not considered drug related’ is considered an adverse effect (55). We will consider ‘adverse event’, ‘adverse drug reaction’, ‘side effect’, ‘toxic effect’, or ‘complication’ as indicating adverse effect. No change or an increase in pain intensity is not considered an adverse effect.
The secondary outcomes are function, serious adverse events, and acceptability:
Function, defined as low back specific function, measured at the end of treatment, may be measured with a continuous, self-report scale (e.g. Roland Morris Disability Questionnaire (RMDQ) or Oswestry Disability Index (ODI)), a rating scale within a composite measure (e.g. Short Form-36), or an ordinal scale (we will consider such ordinal scales to exhibit continuous properties). We will not exclude studies that use other measurement tools.
We will also compare the number of participants who experience a serious adverse effect during the treatment period using the US FDA classification [55, 56], where a serious adverse effect: ‘results in death; is life threatening; requires inpatient hospitalisation or causes prolongation of existing hospitalisation; results in persistent or significant disability/incapacity; may have caused a congenital anomaly/birth defect; or requires intervention to prevent permanent impairment or damage’. A consistent definition of serious adverse effect ensures the different medicines can be compared across the network.
Acceptability is defined as the number of participants who leave the trial for any reason before the end of treatment .
Search strategy and study selection
We will search the following electronic databases from inception to current:
MEDLINE (Ovid) (1946 to current)
EMBASE (Ovid) (1980 to current)
CINAHL (EBSCO) (1982 to current)
Cochrane Central Register of Controlled Trials (CENTRAL) in the Cochrane Library, current issue
EU Clinical Trials Register (eudract.ema.europa.eu)
WHO International Clinical Trial Registry Platform (apps.who.int/trialsearch/Default.aspx)
Our search strategies incorporate the recommended strategies from the Cochrane Back and Neck Group (4) to identify randomised trials of low back pain and terms for the interventions of interest . The search strategy for MEDLINE is listed in Additional file 3. We will also search previous systematic reviews and the reference lists of included studies to identify any additional trials. Records identified through all searches will be downloaded and managed in a custom relational database.
We will conduct record screening in Covidence systematic review software . We will conduct two stages of screening: (i) title and abstract and (ii) full text. Two reviewers will independently screen studies for eligibility at each stage. Disagreements will be resolved through discussion, with arbitration from a third author (JHM) if required. We will contact a study’s corresponding author up to three times to obtain additional information to determine eligibility, and if no reply is received, we will exclude the study from this iteration of the review. Studies in languages other than English will be translated. We will summarise the literature search using an adapted PRISMA flow diagram .
We will manage the included records in the relational database. We will conduct record linkage to establish unique studies for data extraction, which may consist of multiple records. We will search for the protocols and trial registrations of all included trials. We will use the following hierarchy to prioritise records for data extraction: (i) primary report (typically the journal article reporting the results of the primary analysis of the trial), (ii) secondary report (secondary analysis of the trial), (iii) conference abstract (a report of a secondary analysis), (iv) trial registration, (v) other secondary records, and (vi) other conference abstracts.
Two reviewers will independently extract and enter data from included trials into standardised spreadsheets. Review authors will not extract data from any trial in which they have had any involvement. Data will be taken from previous reviews conducted by the authors when possible. Disagreements between reviewers will be resolved through discussion, with arbitration from a third author (JHM) if required. We will not extract data from interventions that do not meet the eligibility criteria for this review.
We will extract data on:
Trial characteristics: country, setting, and number of trial sites; sample size; and study duration.
Participants: diagnosis, duration of LBP, age, male/female ratio, arm-level pain intensity at baseline (as mean (standard deviation [SD])), experience or naivety with the trial intervention, and co-morbidities, including alternate sites of pain.
Interventions: medicine(s) tested, control; duration of intervention; dosage regimen; routes of administration; and usage of rescue medication.
Outcomes: type and dimensions of the scale/measure used to assess pain or function and the time from randomisation at which the end of treatment data were obtained in individual trials. We will extract the definition of ‘adverse effect’ and ‘serious adverse effect’ used in each study. We will extract data on study results including participant allocation to each intervention group; compliance to the intervention (including the definition of compliance); the number of participants who discontinued due to an adverse event; the event rate and descriptions of all reported adverse effects; and pain intensity and function at the completion of treatment.
If studies report more than one measure for pain, we will prioritise extraction in the following order: 100 mm VAS, 10 cm VAS, 11-point NRS, rating scale for pain intensity from a composite measure of pain (e.g. McGill Pain Questionnaire), and ordinal scale. We will preferentially extract the outcome score and measure of variance at the end of treatment (or closest time point) for each group, followed by the change from baseline and measure of variance. If data are not available for each trial arm, we will extract the between-group statistics at the end of treatment.
If studies report more than one measure for function, we will prioritise extraction in the following order: ODI, RMDQ, rating scale for functional ability from a composite measure, and ordinal scale. We will preferentially extract the outcome score and measure of variance at the end of treatment (or closest time point) for each group, followed by the change from baseline and measure of variance. If data are not available for each trial arm, we will extract the between-group statistics at the end of treatment.
We will contact a trial’s corresponding author up to three times via email to request missing data, which will be considered unobtainable if no reply is received within 6 weeks. If data for outcomes of pain and function are not presented in an appropriate form for meta-analysis (such as median and range instead of SDs, standard errors, t-statistics, or p values), we will attempt to impute these using established methods [54, 59]. We will conduct sensitivity analyses for pain at the end of treatment and safety if we impute missing data for either of these outcomes.
Risk of bias
We will appraise each trial’s risk of bias using the Cochrane ‘Risk of bias’ tool, version 5.1  and recommendations by Furlan et al. . Two reviewers will independently appraise trial-level risk of bias for 13 items across the domains of selection, performance, attrition, detection, reporting, and other sources of bias. If an item is typically rated at outcome level, which may differ between our two primary outcomes (pain intensity and safety), we will use the more conservative rating (e.g. using high risk over unclear risk). Review authors will not appraise risk of bias for any trial in which they have had any involvement (e.g. trial investigator). Risk of bias assessments will be taken from previous reviews of analgesic medicines conducted by our author team, where the same approach was used.
We will determine an overall risk of bias for each trial by adapting the process from Furukawa et al. : low overall risk is determined when three or fewer items are rated ‘unclear’ risk and no domains are rated ‘high’; moderate overall risk is determined if a single item is rated as ‘high’ risk of bias, or no item is rated as ‘high’ risk but four or more are rated as ‘unclear’; and high overall risk otherwise.
We will perform separate analyses for the three classifications of pain duration: acute (fewer than 6 weeks), sub-acute (6 to 12 weeks), and chronic (greater than 12 weeks).
Summary of the network
Within each classification, we will present descriptive statistics for each included trial, including the comparison(s) and a clinical/methodological summary (e.g. year of publication, sponsorship, clinical setting). We will represent the network of trials in a network graph; the size of the node will reflect the total number of participants, the width of each edge will reflect the number of studies presenting direct evidence for the comparison, and the colour of each edge will represent the overall risk of bias (see the ‘Risk of bias’ section).
We will synthesise the data for each comparison using pairwise random-effects meta-analysis in R . We will compare the effects of competing interventions on pain and function using mean differences (MD) with 95% confidence intervals (CIs) and on safety and acceptability using odds ratios (OR) with 95% CIs. For pain intensity and function, we will convert outcome data to common 0- to 100-point scales (mean (SD)), which has been used in reviews of analgesic medicines to enable greater clinical translation of results [25, 29, 33, 36]. This approach is based on evidence that measurement scales within the constructs of pain intensity and function are highly correlated [61, 62] and enables different types of study data to be pooled (e.g. endpoint, change from baseline).
We assume that the heterogeneity variance is different for each comparison in pairwise meta-analyses and will estimate this parameter (τ2) for each comparison. We will also test for the presence of statistical heterogeneity within each comparison using the Q statistic. We will calculate 95% prediction intervals and consider intervals spanning greater than 15 points (on a 0- to 100-point scale) on either side to indicate important heterogeneity . We will visually inspect the distribution of effect sizes in the forest plots for each comparison and consider an I2 value greater than 50% indicative of important variability across studies that is not due to sampling error [63, 64].
Assumption of transitivity
Transitivity is the key assumption underlying the valid estimation of effects for indirect comparisons in NMA [38, 39]. Transitivity is the assumption that the distributions of effect modifiers (covariates associated with intervention effects) are balanced across comparisons in the network [39, 65]. Given the lack of evidence for robust effect modifiers in LBP trials , we have used clinical and methodological experience to identify the following potential effect modifiers:
Baseline mean pain intensity (continuous variable)
Assigned co-interventions (dichotomous variable), categorised as (i) yes or (ii) no
Small sample size  (dichotomous variable, categorised as (i) total sample fewer than 50 participants and (ii) total sample greater than 50 participants)
Experience with test medicine (dichotomous variable), categorised as (i) yes or (ii) no
Naivety to test medicine (dichotomous variable), categorised as (i) yes or (ii) no
Dose of medicine (trichotomous variable), categorised as (i) SDR, (ii) below the SDR, or (iii) above the SDR
We will represent the distribution of these effect modifiers in a range of covariate contribution plots . The authors will visually assess the distributions of effect modifiers across all treatment comparisons in the network and determine by consensus whether there is sufficient dissimilarity between comparisons in the network to threaten the assumption of transitivity. We will explore the influence of effect modifiers that demonstrate dissimilarity on incoherence/heterogeneity using network meta-regression or subgroup analyses (or both). We will consider not proceeding with NMA, or altering the network structure, if we observe considerable dissimilarity. We anticipate that insufficient reporting of effect modifiers and pairwise comparisons containing few studies will limit the assessment of transitivity .
We will perform random-effects NMA within an electrical network and graph theory framework using the netmeta package in R [69, 70]. We will account for multi-arm trials using the weighting method based on back-calculating variances (using the Laplacian matrix and its pseudoinverse) . We will use MDs on a common 0 to 100-point scale for pain and function and ORs for safety and acceptability. We will present the results for each intervention compared to placebo in NMA forest plots for each outcome. We will consider a 10-point difference to constitute the minimal clinically important difference for pain intensity and function .
We will rank the effect of all interventions on pain intensity and safety using P-scores, a method for estimating treatment rankings that does not require re-sampling methods . We will present a contributions matrix to indicate the weighting of direct evidence contributions to each NMA effect size, which will also be used to evaluate the confidence in the overall evidence (see the ‘Confidence in cumulative evidence’ section) . We will illustrate the contribution of each design to an NMA effect size using net heat plots .
Assessment of network heterogeneity and coherence
We will assume a common heterogeneity variance across the network . We will present the estimate for this parameter (τ2network) from the NMA models and the estimated proportion of variability across the entire network that is not due to sampling error (I2network). We will estimate the Q statistics for total network heterogeneity (Qtotal), heterogeneity within designs (Qwithin), and heterogeneity between designs (Qbetween), where designs constitute the individual elements of the set of trial designs.
Coherence is a property of closed loops of evidence, whereby it reflects the agreement between direct and indirect treatment effects . We will evaluate coherence across the entire network using the Q statistics (above); the decomposed Qwithin and Qbetween, an alternative estimate for Qbetween using the ‘design-by-treatment’ interaction model [75, 76]; and the Separating Indirect from Direct Evidence (SIDE; aka node-splitting) approach . We will illustrate the extent of incoherence across the network using net heat plots . We will illustrate local coherence estimates using forest plots grouped into direct and indirect estimates for each available comparison. We will form judgements about important incoherence using all the measures of global and local heterogeneity and coherence.
If we encounter important incoherence, we will examine the dataset for data extraction errors and explore the observed incoherence using pre-specified covariates in network meta-regression and subgroup analyses, provided sufficient studies are available. We may consider not proceeding with NMA if important unexplained incoherence remains. This judgement will involve the clinical and methodological evaluation of transitivity, the approaches to identify incoherence, and the knowledge that small amounts of incoherence may be due to chance .
We aim to perform random-effects network meta-regression within a Bayesian hierarchical framework using the gemtc package in R . The gemtc package will automatically determine uninformative prior distributions for all parameters in our model , which are commonly applied in NMA [47, 81]. We will run the Markov Chain Monte Carlo simulation with four chains for each model, using 100,000 iterations, a burn-in of 5000 iterations and extraction of every 10th value. We will assess convergence with the Gelman-Rubin-Brooks plot and Potential Scale Reduction Factor (a threshold of < 1.05 indicates adequate convergence).
We will investigate baseline pain intensity (continuous variable) and sample size (dichotomous variable) as possible sources of incoherence or heterogeneity by default. We will specify our assumptions (common or exchangeable covariate-comparison interaction) for each network meta-regression model once the data are available, so as to make best use of the available data [82,83,84]. We assume coherent relative treatment effects estimated at the covariate value 0 and coherent regression coefficients for the treatment effect by covariate interaction [84, 85]. We hypothesise that:
Increasing baseline pain intensity increases the effect size between intervention and comparator.
Increasing sample size reduces the effect size between intervention and comparator.
We may also investigate clinical or methodological factors identified during the review process that may threaten transitivity as sources of incoherence or heterogeneity, or both .
We will conduct sensitivity analyses on pain and safety excluding studies at high risk of bias, provided the original network structure remains the same. We will also conduct sensitivity analyses on pain and safety excluding doses above or below the SDR, provided the original network structure remains the same. If sufficient data are not available for network meta-regression baseline pain intensity or sample size, we will conduct sensitivity analyses by removing trials with baseline mean pain intensity higher than 70/100 (VAS) or sample size less than 50, respectively. We will also conduct sensitivity analyses for pain at end of treatment and safety if we impute missing data for either of these outcomes. We will consider the effects in the sensitivity analyses to be important when their interpretation differs compared to the primary analysis, for example, if a statistically significant effect becomes non-significant.
We will assess small study effects in pairwise comparisons using comparison-adjusted  and contour-enhanced  funnel plots when there are at least 10 studies available. Such plots assist interpretation of asymmetry that is due to publication bias rather than other factors, such as lesser methodological quality. We will interpret an absence of studies in areas of non-significance as suggestive of publication bias for that pairwise comparison.
Confidence in cumulative evidence
We will form judgements of confidence in the effect estimates and rankings for pain, safety, function, and acceptability using the Confidence in Network Meta-Analysis (CINeMA) web application [88, 89]. CINeMA considers six domains: within-study bias, reporting bias, indirectness, imprecision, heterogeneity, and incoherence (Additional file 4) . Initially, judgements will be rated as ‘high’ because all included trials will be RCTs.
Summary of findings
We will present the results of the NMA in adapted ‘Summary of findings’ tables for pain and safety [90, 91]. The tables will contain details of the clinical question, a network geometry plot, relative and absolute effect estimates, certainty of evidence, ranking of treatments, and interpretation of findings.
Availability of data and materials
Data sharing is not applicable to this article as no datasets were generated or analysed during the current study.
Anatomical Therapeutic Chemical
Confidence in Network Meta-Analysis
European Medicines Agency
Licensed dosing range
Low back pain
Medicine and Healthcare Products Regulatory Agency
Non-steroidal anti-inflammatory drug
Numeric rating scale
Oswestry Disability Index
Preferred Reporting Items for Systematic Reviews and Meta-Analysis
Roland Morris Disability Questionnaire
Standard dosing range
Therapeutic Goods Administration
US Food and Drug Administration
Visual analogue scale
World Health Organization
James SLAD, Abate KH, Abay SM, Abbafati C, Abbasi N, et al. Global, regional, and national incidence, prevalence, and years lived with disability for 354 diseases and injuries for 195 countries and territories, 1990-2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet. 2018;392(10159):1789–858.
Koes BW, van Tulder MW, Thomas S. Diagnosis and treatment of low back pain. BMJ. 2006;332(7555):1430–4.
Maher C, Underwood M, Buchbinder R. Non-specific low back pain. Lancet. 2017;389(10070):736–47.
Furlan AD, Malmivaara A, Chou R, Maher CG, Deyo RA, Schoene M, et al. 2015 Updated Method Guideline for Systematic Reviews in the Cochrane Back and Neck Group. Spine. 2015;40(21):1660–73.
Treede R-D, Rief W, Barke A, Aziz Q, Bennett MI, Benoliel R, et al. A classification of chronic pain for ICD-11. Pain. 2015;156(6):1003–7.
Rubin DI. Epidemiology and risk factors for spine pain. Neurol Clin. 2007;25(2):353–71.
Hoy D, Bain C, Williams G, March L, Brooks P, Blyth F, et al. A systematic review of the global prevalence of low back pain. Arthritis Rheumatol. 2012;64(6):2028–37.
da C Menezes Costa L, Maher CG, Hancock MJ, McAuley JH, Herbert RD, Costa LOP. The prognosis of acute and persistent low-back pain: a meta-analysis. Can Med Assoc J. 2012;184(11):E613–E24.
Henschke N, Maher CG, Refshauge KM, Herbert RD, Cumming RG, Bleasel J, et al. Prognosis in patients with recent onset low back pain in Australian primary care: inception cohort study. BMJ. 2008;337:a171.
da Silva T, Mills K, Brown BT, Pocovi N, de Campos T, Maher C, et al. Recurrence of low back pain is common: a prospective inception cohort study. J Physiother. 2019;65(3):159–65.
Oliveira CB, Maher CG, Pinto RZ, Traeger AC, Lin C-WC, Chenot J-F, et al. Clinical practice guidelines for the management of non-specific low back pain in primary care: an updated overview. Eur Spine J. 2018;27(11):2791–803.
Mathieson S, Valenti L, Maher CG, Britt H, Li Q, McLachlan AJ, et al. Worsening trends in analgesics recommended for spinal pain in primary care. Eur Spine J. 2018;27(5):1136–45.
Bishop PB, Wing PC. Compliance with clinical practice guidelines in family physicians managing worker's compensation board patients with acute lower back pain. Spine J. 2003;3(6):442–50.
Ivanova JI, Birnbaum HG, Schiller M, Kantor E, Johnstone BM, Swindle RW. Real-world practice patterns, health-care utilization, and costs in patients with low back pain: the long road to guideline-concordant care. Spine J. 2011;11(7):622–32.
Piccoliori G, Engl A, Gatterer D, Sessa E, in der Schmitten J, Abholz HH. Management of low back pain in general practice - is it of acceptable quality: an observational study among 25 general practices in South Tyrol (Italy). BMC Fam Pract. 2013;14:148.
Williams CM, Maher CG, Hancock MJ, McAuley JH, McLachlan AJ, Britt H, et al. Low back pain and best practice care: a survey of general practice physicians. Arch Intern Med. 2010;170(3):271–7.
Ferreira GE, Machado GC, Abdel Shaheed C, Lin C-WC, Needs C, Edwards J, et al. Management of low back pain in Australian emergency departments. BMJ Qual Saf. 2019;28(10):826.
Friedman BW, Chilstrom M, Bijur PE, Gallagher EJ. Diagnostic testing and treatment of low back pain in United States emergency departments: a national perspective. Spine. 2010;35(24):E1406–E11.
Gore M, Tai KS, Sadosky A, Leslie D, Stacey BR. Use and costs of prescription medications and alternative treatments in patients with osteoarthritis and chronic low back pain in community-based settings. Pain Pract. 2012;12(7):550–60.
Hart OR, Uden RM, McMullan JE, Ritchie MS, Williams TD, Smith BH. A study of National Health Service management of chronic osteoarthritis and low back pain. Prim Health Care Res Dev. 2015;16(2):157–66.
Del Fiol G, Workman TE, Gorman PN. Clinical questions raised by clinicians at the point of care: a systematic review. JAMA Int Med. 2014;174(5):710–8.
Chou L, Ranger TA, Peiris W, Cicuttini FM, Urquhart DM, Sullivan K, et al. Patients’ perceived needs of health care providers for low back pain management: a systematic scoping review. Spine J. 2018;18(4):691–711.
Chou L, Ranger TA, Peiris W, Cicuttini FM, Urquhart DM, Sullivan K, et al. Patients’ perceived needs for medical services for non-specific low back pain: a systematic scoping review. PLoS One. 2018;13(11):e0204885.
Lim YZ, Chou L, Au RTM, Seneviwickrama KLMD, Cicuttini FM, Briggs AM, et al. People with low back pain want clear, consistent and personalised information on prognosis, treatment options and self-management strategies: a systematic review. J Physiother. 2019;65(3):124–35.
Abdel Shaheed C, Maher CG, Williams KA, McLachlan AJ. Efficacy and tolerability of muscle relaxants for low back pain: systematic review and meta-analysis. Eur J Pain. 2017;21(2):228–37.
Chou R, Deyo R, Friedly J, Skelly A, Weimer M, Fu R, et al. Systemic pharmacologic therapies for low back pain: a systematic review for an American College of Physicians clinical practice guideline. Ann Int Med. 2017;166(7):480–92.
Deyo RA, Von Korff M, Duhrkoop D. Opioids for low back pain. BMJ. 2015;350:g6380.
Enke O, New HA, New CH, Mathieson S, McLachlan AJ, Latimer J, et al. Anticonvulsants in the treatment of low back pain and lumbar radicular pain: a systematic review and meta-analysis. Can Med Assoc J. 2018;190(26):E786–E93.
Machado GC, Maher CG, Ferreira PH, Day RO, Pinheiro MB, Ferreira ML. Non-steroidal anti-inflammatory drugs for spinal pain: a systematic review and meta-analysis. Ann Rheum Dis. 2017;76(7):1269–78.
Qaseem A, Wilt TJ, McLean RM, Forciea MA. Noninvasive treatments for acute, subacute, and chronic low back pain: a clinical practice guideline from the American College of Physicians. Ann Intern Med. 2017;166(7):514–30.
Saragiotto BT, Machado GC, Ferreira ML, Pinheiro MB, Abdel Shaheed C, Maher CG. Paracetamol for low back pain. Cochrane Database Syst Rev. 2016(6):CD012230.
Urquhart DM, Hoving JL, Assendelft WJJ, Roland M, van Tulder MW. Antidepressants for non-specific low back pain. Cochrane Database Syst Rev. 2008(1):CD001703.
Abdel Shaheed C, Maher CG, Williams KA, Day R, McLachlan AJ. Efficacy, tolerability, and dose-dependent effects of opioid analgesics for low back pain: a systematic review and meta-analysis. JAMA Int Med. 2016;176(7):958–68.
Seidel S, Aigner M, Ossege M, Pernicka E, Wildner B, Sycha T. Antipsychotics for acute and chronic pain in adults. Cochrane Database Syst Rev. 2013(8):CD004844.
Tucker H-R, Scaff K, McCloud T, Carlomagno K, Daly K, Garcia A, et al. Harms and benefits of opioids for management of non-surgical acute and chronic low back pain: a systematic review. Br J Sports Med. 2020;54:664.
Bagg MK, McLachlan AJ, Maher CG, Kamper SJ, Williams CM, Henschke N, et al. Paracetamol, NSAIDS and opioid analgesics for chronic low back pain: a network meta-analysis. Cochrane Database Syst Rev. 2018(6):CD013045.
Bagg MK, Salanti G, McAuley JH. Comparing interventions with network meta-analysis. J Physiother. 2018;64(2):128–32.
Cipriani A, Higgins JP, Geddes JR, Salanti G. Conceptual and technical challenges in network meta-analysis. Ann Intern Med. 2013;159(2):130–7.
Salanti G. Indirect and mixed-treatment comparison, network, or multiple-treatments meta-analysis: many names, many benefits, many concerns for the next generation evidence synthesis tool. Res Synth Methods. 2012;3(2):80–97.
Moher D, Shamseer L, Clarke M, Ghersi D, Liberati A, Petticrew M, et al. Preferred reporting items for systematic review and meta-analysis protocols (PRISMA-P) 2015 statement. Syst Rev. 2015;4(1):1.
Shamseer L, Moher D, Clarke M, Ghersi D, Liberati A, Petticrew M, et al. Preferred reporting items for systematic review and meta-analysis protocols (PRISMA-P) 2015: elaboration and explanation. BMJ. 2015;349:g7647.
Bagg MK, O’Hagan E, Zahara P, Wand BM, Hübscher M, Moseley GL, et al. Reviews may overestimate the effectiveness of medicines for back pain: systematic review and meta-analysis. J Clin Epidemiol. 2019;124:149–59.
Turner EH, Matthews AM, Linardatos E, Tell RA, Rosenthal R. Selective publication of antidepressant trials and its influence on apparent efficacy. N Engl J Med. 2008;358(3):252–60.
Spineli LM, Higgins JP, Cipriani A, Leucht S, Salanti G. Evaluating the impact of imputations for missing participant outcome data in a network meta-analysis. Clin Trials. 2013;10(3):378–88.
Elbourne DR, Altman DG, Higgins JP, Curtin F, Worthington HV, Vail A. Meta-analyses involving cross-over trials: methodological issues. Int J Epidemiol. 2002;31(1):140–9.
Cipriani A, Furukawa TA, Salanti G, Chaimani A, Atkinson LZ, Ogawa Y, et al. Comparative efficacy and acceptability of 21 antidepressant drugs for the acute treatment of adults with major depressive disorder: a systematic review and network meta-analysis. Lancet. 2018;391(10128):1357–66.
Furukawa TA, Salanti G, Atkinson LZ, Leucht S, Ruhe HG, Turner EH, et al. Comparative efficacy and acceptability of first-generation and second-generation antidepressants in the acute treatment of major depression: protocol for a network meta-analysis. BMJ Open. 2016;6(7):e010919.
Furlan A, Chaparro LE, Irvin E, Mailis-Gagnon A. A comparison between enriched and nonenriched enrollment randomized withdrawal trials of opioids for chronic noncancer pain. Pain Res Manag. 2011;16(5):337–51.
Yamato TP, Maher CG, Saragiotto BT, Shaheed CA, Moseley AM, Lin C-WC, et al. Comparison of effect sizes between enriched and nonenriched trials of analgesics for chronic musculoskeletal pain: a systematic review. Br J Clin Pharmacol. 2017;83(11):2347–55.
Koes BW, van Tulder MW, Peul WC. Diagnosis and treatment of sciatica. BMJ. 2007;334(7607):1313–7.
PDR. Prescriber’s Digital Reference. Available from https://www.pdr.net.
MIMS. MIMS. Available from https://www.mims.co.uk.
AMH. Australian Medicines Handbook. Available from https://amhonline.amh.net.au.
Higgins JPT, Thomas J, Chandler J, Cumpston M, Li T, Page MJ, et al. Cochrane Handbook for Systematic Reviews of Interventions Version 6.0 (updated July 2019). Cochrane; 2011.
US Food and Drug Administration. Code of Federal Regulations Title 21. 2018.
Tomlinson A, Efthimiou O, Boaden K, New E, Mather S, Salanti G, et al. Side effect profile and comparative tolerability of 21 antidepressants in the acute treatment of major depression in adults: protocol for a network meta-analysis. Evid Based Ment Health. 2019;22(2):61–6.
Innovation VH. Covidence systematic review software. Melbourne, Australia: Veritas Health Innovation. Available at www.covidence.org.
Liberati A, Altman DG, Tetzlaff J, Mulrow C, Gotzsche PC, Ioannidis JP, et al. The PRISMA statement for reporting systematic reviews and meta-analyses of studies that evaluate healthcare interventions: explanation and elaboration. BMJ. 2009;339:b2700.
Wan X, Wang W, Liu J, Tong T. Estimating the sample mean and standard deviation from the sample size, median, range and/or interquartile range. BMC Med Res Methodol. 2014;14:135.
R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2019.
Hjermstad MJ, Fayers PM, Haugen DF, Caraceni A, Hanks GW, Loge JH, et al. Studies comparing Numerical Rating Scales, Verbal Rating Scales, and Visual Analogue Scales for assessment of pain intensity in adults: a systematic literature review. J Pain Symptom Manage. 2011;41(6):1073–93.
Smeets R, Köke A, Lin C-W, Ferreira M, Demoulin C. Measures of function in low back pain/disorders: Low Back Pain Rating Scale (LBPRS), Oswestry Disability Index (ODI), Progressive Isoinertial Lifting Evaluation (PILE), Quebec Back Pain Disability Scale (QBPDS), and Roland-Morris Disability Questionnaire (RDQ). Arthritis Care Res (Hoboken). 2011;63(S11):S158–S73.
Borenstein M, Hedges LV, Higgins JP, Rothstein HR. Introduction to Meta-Analysis. Chichester: Wiley; 2009.
Borenstein M, Higgins JPT, Hedges LV, Rothstein HR. Basics of meta-analysis: I2 is not an absolute measure of heterogeneity. Res Synth Methods. 2017;8(1):5–18.
Jansen JP, Naci H. Is network meta-analysis as valid as standard pairwise meta-analysis? It all depends on the distribution of effect modifiers. BMC Med. 2013;11(1):159.
Saragiotto BT, Maher CG, Moseley AM, Yamato TP, Koes BW, Sun X, et al. A systematic review reveals that the credibility of subgroup claims in low back pain trials was low. J Clin Epidemiol. 2016;79:3–9.
Dechartres A, Trinquart L, Faber T, Ravaud P. Empirical evaluation of which trial characteristics are associated with treatment effect estimates. J Clin Epidemiol. 2016;77:24–37.
Donegan S, Dias S, Tudur-Smith C, Marinho V, Welton NJ. Graphs of study contributions and covariate distributions for network meta-regression. Res Synth Methods. 2018;9(2):243–60.
Rücker G, Krahn U, König J, Efthimiou O, Schwarzer G. netmeta: Network meta-analysis using frequentist methods. R package version 1.1-0. 2019. https://CRAN.R-project.org/package=netmeta.
Rücker G. Network meta-analysis, electrical networks and graph theory. Res Synth Methods. 2012;3(4):312–24.
Rücker G, Schwarzer G. Ranking treatments in frequentist network meta-analysis works without resampling methods. BMC Med Res Methodol. 2015;15(1):58.
Salanti G, Del Giovane C, Chaimani A, Caldwell DM, Higgins JPT. Evaluating the quality of evidence from a network meta-analysis. PLoS One. 2014;9(7):e99682.
Schwarzer G, Carpenter JR, Rücker G. Network Meta-Analysis. In: Schwarzer G, Carpenter JR, Rücker G, editors. Meta-Analysis with R. Cham: Springer International Publishing; 2015. p. 187–216.
Higgins JPT, Whitehead A. Borrowing strength from external trials in a meta-analysis. Stat Med. 1996;15(24):2733–49.
Higgins JPT, Jackson D, Barrett JK, Lu G, Ades AE, White IR. Consistency and inconsistency in network meta-analysis: concepts and models for multi-arm studies. Res Synth Methods. 2012;3(2):98–110.
White IR, Barrett JK, Jackson D, Higgins JPT. Consistency and inconsistency in network meta-analysis: model estimation using multivariate meta-regression. Res Synth Methods. 2012;3(2):111–25.
Dias S, Welton NJ, Caldwell DM, Ades AE. Checking consistency in mixed treatment comparison meta-analysis. Stat Med. 2010;29(7-8):932–44.
Veroniki AA, Vasiliadis HS, Higgins JPT, Salanti G. Evaluation of inconsistency in networks of interventions. Int J Epidemiol. 2013;42(1):332–45.
van Valkenhoef G, Kuiper J. Gemtc: Network Meta-Analysis Using Bayesian Methods; 2016. https://cran.r-project.org/web/packages/gemtc/index.html.
van Valkenhoef G, Lu G, de Brock B, Hillege H, Ades AE, Welton NJ. Automating network meta-analysis. Res Synth Methods. 2012;3(4):285–99.
Huhn M, Nikolakopoulou A, Schneider-Thoma J, Krause M, Samara M, Peter N, et al. Comparative efficacy and tolerability of 32 oral antipsychotics for the acute treatment of adults with multi-episode schizophrenia: a systematic review and network meta-analysis. Lancet. 2019;394(10202):939–51.
Chaimani A, Salanti G. Using network meta-analysis to evaluate the existence of small-study effects in a network of interventions. Res Synth Methods. 2012;3(2):161–76.
Efthimiou O, Debray TPA, van Valkenhoef G, Trelle S, Panayidou K, Moons KGM, et al. GetReal in network meta-analysis: a review of the methodology. Res Synth Methods. 2016;7(3):236–63.
Donegan S, Welton NJ, Tudur Smith C, D'Alessandro U, Dias S. Network meta-analysis including treatment by covariate interactions: consistency can vary across covariate values. Res Synth Methods. 2017;8(4):485–95.
Donegan S, Dias S, Welton NJ. Assessing the consistency assumptions underlying network meta-regression using aggregate data. Res Synth Methods. 2019;10(2):207–24.
Jansen JP, Schmid CH, Salanti G. Directed acyclic graphs can help understand bias in indirect and mixed treatment comparisons. J Clin Epidemiol. 2012;65(7):798–807.
Peters JL, Sutton AJ, Jones DR, Abrams KR, Rushton L. Contour-enhanced meta-analysis funnel plots help distinguish publication bias from other causes of asymmetry. J Clin Epidemiol. 2008;61(10):991–6.
ISPM. CINeMA: confidence in network meta-analysis. cinema.ispm.unibe.ch: Institute of Social and Preventive Medicine, University of Bern; 2017.
Nikolakopoulou A, Higgins JPT, Papakonstantinou T, Chaimani A, Del Giovane C, Egger M, et al. CINeMA: An approach for assessing confidence in the results of a network meta-analysis. PLoS Med. 2020:17(4):e1003082.
Brignardello-Petersen R, Bonner A, Alexander PE, Siemieniuk RA, Furukawa TA, Rochwerg B, et al. Advances in the GRADE approach to rate the certainty in estimates from a network meta-analysis. J Clin Epidemiol. 2018;93:36–44.
Yepes-Nuñez JJ, Li S-A, Guyatt G, Jack SM, Brozek JL, Beyene J, et al. Development of the summary of findings table for network meta-analysis. J Clin Epidemiol. 2019;115:1–13.
The ANiMALIA investigators: Aidan G Cashin (Centre for Pain IMPACT, Neuroscience Research Australia/Prince of Wales Clinical School, University of New South Wales); Richard O Day (Clinical Pharmacology & Toxicology, St. Vincent’s Hospital/ St. Vincent’s Clinical School, University of New South Wales); Michael C Ferraro (Centre for Pain IMPACT, Neuroscience Research Australia/School of Medical Sciences, University of New South Wales); Thiago Folly (Centre for Pain IMPACT, Neuroscience Research Australia); Sylvia M Gustin (Centre for Pain IMPACT, Neuroscience Research Australia/ School of Psychology, University of New South Wales); Amanda D Hagstrom (School of Medical Sciences, University of New South Wales); Hayley B Leake (IIMPACT in Health, University of South Australia); Colleen K Loo (School of Psychiatry, University of New South Wales/ Black Dog Institute); Andrew J McLachlan (The University of Sydney School of Pharmacy); Edel T O’Hagan (Centre for Pain IMPACT, Neuroscience Research Australia/ Prince of Wales Clinical School, University of New South Wales); Rodrigo RN Rizzo (Centre for Pain IMPACT, Neuroscience Research Australia/School of Medical Sciences, University of New South Wales); Siobhan M Schabrun (Centre for Pain IMPACT, Neuroscience Research Australia).
This project did not receive funding.
MAW was supported by a University Postgraduate Award and School of Medical Sciences Top-Up Scholarship from the University of New South Wales and a Postgraduate Scholarship from the National Health and Medical Research Council of Australia. MKB was supported by a NeuRA PhD Candidature Scholarship and Supplementary Scholarship. AGC was supported by the University of New South Wales Prince of Wales Clinical School Postgraduate Research Scholarship and a NeuRA PhD Candidature Supplementary Scholarship. MCF was supported by an Australian Medical Research Future Fund Grant GNTID1170205. HBL was supported by Australian Government Research Training Program Scholarship. ETO was supported by an Australian Government Research Training Program Scholarship and a NeuRA PhD Candidature Supplementary Scholarship. RRNR was supported by the University of New South Wales School of Medical Sciences Postgraduate Research Scholarship and a NeuRA PhD Candidature Supplementary Scholarship.
Ethics approval and consent to participate
Ethics approval is not required for this study.
Consent for publication
MAW, MDJ, and JHM declare no competing interests.
MKB received travel support from Memorial University of Newfoundland to speak about engagement with research evidence, including the effects of medicines.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Wewege, M.A., Bagg, M.K., Jones, M.D. et al. Analgesic medicines for adults with low back pain: protocol for a systematic review and network meta-analysis. Syst Rev 9, 255 (2020). https://doi.org/10.1186/s13643-020-01506-3
- Analgesics [MeSH]
- Low back pain [MeSH]
- Meta-analysis as topic [MeSH]
- Network meta-analysis [MeSH]
- Systematic reviews as topic [MeSH]