You are here

  • Home
  • Archive
  • Volume 10, Issue 11
  • Evaluating the intended and unintended consequences of opioid-prescribing interventions on primary care in British Columbia, Canada: protocol for a retrospective population-based cohort study

Article Text

Article menu

Download PDFPDF

XML
Addiction
Protocol
Evaluating the intended and unintended consequences of opioid-prescribing interventions on primary care in British Columbia, Canada: protocol for a retrospective population-based cohort study
  1. Dimitra Panagiotoglou1,
  2. Rita McCracken2,
  3. M Ruth Lavergne3,4,
  4. Erin C Strumpf1,5,
  5. Tara Gomes6,7,
  6. Benedikt Fischer8,9,10,
  7. Austyn Brackett11,
  8. Cheyenne Johnson12,13,
  9. Perry Kendall12
  1. 1Department of Epidemiology, Biostatistics and Occupational Health, McGill University, Montreal, Québec, Canada
  2. 2Department of Family Practice, University of British Columbia Faculty of Medicine, Vancouver, British Columbia, Canada
  3. 3Centre for Applied Research in Mental Health and Addiction, Simon Fraser University, Burnaby, British Columbia, Canada
  4. 4Faculty of Health Sciences, Simon Fraser University, Burnaby, British Columbia, Canada
  5. 5Department of Economics, McGill University, Montreal, Québec, Canada
  6. 6Li Ka Shing Knowledge Institute, St. Michael's Hospital, Toronto, Ontario, Canada
  7. 7Institute for Clinical Evaluative Sciences, University of Toronto, Toronto, Ontario, Canada
  8. 8Faculty of Medical and Health Sciences, The University of Auckland, Auckland, New Zealand
  9. 9Institute for Mental Health Policy Research, Centre for Addiction and Mental, University of Toronto, Toronto, Ontario, Canada
  10. 10Department of Psychiatry, University of Toronto, Toronto, Ontario, Canada
  11. 11Patient Voices Network, Vancouver, British Columbia, Canada
  12. 12British Columbia Centre on Substance Use, Vancouver, British Columbia, Canada
  13. 13School of Nursing, University of British Columbia, Vancouver, British Columbia, Canada
  1. Correspondence to Dr Dimitra Panagiotoglou; dimitra.panagiotoglou{at}mcgill.ca

Abstract

Introduction Between 2015 and 2018, there were over 40 000 opioid-related overdose events and 4551 deaths among residents in British Columbia (BC). During this time the province mobilised a variety of policy levers to encourage physicians to expand access to opioid agonist treatment and the College of Physicians and Surgeons of British Columbia (CPSBC) released a practice standard establishing legally enforceable minimum thresholds of professional behaviour in the hopes of curtailing overdose events. Our goal is to conduct a comprehensive investigation of the intended and unintended consequences of these policy changes. Specifically, we aim to understand the effects of these measures on physician prescribing behaviours, identify physician characteristics associated with uptake of the new measures, and measure the effects of the policy changes on patients’ access to quality primary care.

Methods and analysis This is a population-level, retrospective cohort study of all BC primary care physicians who prescribed any opioid medication for opioid-use disorder or chronic non-cancer pain during the study period, and their patients. The study period is 1 January 2013–31 December 2018, with a 1-year wash-in period (1 January 2012–31 December 2012) to exclude patients who initiated long-term opioid treatment prior to our study period or whose pain type (ie, ‘chronic non-cancer’, ‘acute’, ‘cancer or palliative’, or ‘other’) cannot be confirmed. The project combines five administrative health datasets under the authority of the BC Ministry of Health, with the CPSBC’s Physician Registry, BC Cancer Agency’s Cancer Registry and Vital Statistics’ Mortality data. We will create measures of prescribing concordance, access, continuity, and comprehensiveness to assess primary care delivery and quality at both the physician and patient level. We will use generalised estimating equations, interrupted time series, mixed effects models, and funnel plots to identify factors related to changes in prescribing and evaluate the impact of the changes to prescribing policies. Results will be reported using appropriate Enhancing the QUAlity and Transparency Of health Research guidelines (eg, STrengthening the Reporting of OBservational studies in Epidemiology).

Ethics and dissemination This study has been approved by McGill University’s Institutional Review Board (#A11-M55-19A), and the University of British Columbia’s Research Ethics Board (#H19-03537). We will disseminate results via a combination of open access peer-reviewed journal publications, conferences, lay summaries and OpEds.

  • protocols & guidelines
  • quality in health care
  • health policy
http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

https://doi.org/10.1136/bmjopen-2020-038724

Statistics from Altmetric.com

    Request Permissions

    If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

    Strengths and limitations of this study

    • This is a population-level study of all primary care physicians in British Columbia (approximately 6000, 1200 of which prescribed at least one opioid agonist treatment during the study period) and their patients (>4.8 million) over 6 years (2013–2018).

    • By working with a comprehensive dataset, master drug lists and validated coding algorithms in a context where the majority of health services are provided at no cost and all prescription drug dispensations are recorded, we have minimised the potential for misclassification and capture as many patients as possible outside of a strict longitudinal cohort study.

    • We use a combination of methods and explore the effects of the prescribing policies from both the physician and patient perspectives.

    • We are unable to operationalise and evaluate adherence to all stipulations in the practice standard (ie, documentation of discussions for non-pharmaceutical alternatives or take-home naloxone kits, pill counts or urine tests).

    Introduction

    Between 2015 and 2018, there were over 40 000 opioid-related overdose events and 4551 deaths among residents in British Columbia (BC).1–4 The annual provincial overdose mortality rate rose from 4.7 per 100 000 population in 2009 to 30.8 per 100 000 population in 2018, and for Vancouver’s Downtown Eastside exceeded 100 deaths per 100 000 population.5 6 In an attempt to curtail the number of overdose events, the province launched new policies to encourage physicians and nurse practitioners to expand access to opioid agonist treatment (OAT) (OAT refers to a set of pharmacological antagonists (eg, methadone, buprenorphine or naltrexone) for addiction to opioids such as heroin, fentanyl, hydromorphone and oxycodone) including open listing buprenorphine/naloxone in the PharmaCare and First Nations Health Benefits programmes.7 Further, the province amended existing fees and introduced new ones to better compensate physicians who provide treatment for opioid-use disorder (OUD).8 9 Concurrently, the BC Centre on Substance Use (BCCSU) developed provincial guidelines and began offering comprehensive education and training for prescribers of OAT.10

    Meanwhile, the College of Physicians and Surgeons of British Columbia (CPSBC) released the Safe Prescribing of Drugs with Potential for Misuse/Diversion practice standard (box 1) to ‘prevent an increasing toll of prescription drug misuse and overdose deaths’.11 Recently renamed Safe Prescribing of Opioids and Sedatives, the CPSBC’s practice standard evolved from the Centers for Disease Control and Prevention Guideline for Prescribing Opioids for Chronic Pain–USA, 2016, the College’s Prescribing Principles (2012), and the National Opioid Use Guideline Group’s Safe and Effective Use of Opioids for Chronic Non-Cancer Pain (2010). The practice standard reflects the findings that opioids compared with placebo modestly improve physical functioning and quality of life for people living with chronic non-cancer pain (lasting longer than 3 months); but that the risks of adverse health effects outweigh the benefits.12–14 Unlike its predecessors which recommended courses of action and allowed physicians to ‘exercise reasonable discretion in their decision to act on guidance provided’, the practice standard established legally enforceable minimum thresholds of professional behaviour.15 Non-compliant physicians can be disciplined or fined under the Health Professions Act, RSBC 1996, c.183 and College Bylaws.15 16

    Box 1 Safe prescribing of drugs with potential for misuse/diversion practice standard, June 2016

    Physicians must:

    • Review patients’ current medications before prescribing opioids, sedatives or stimulants.

    • Base long-term treatment with medications with known risks, including opioids, sedatives and stimulants, on clinical diagnosis and objective evidence. Continuing to prescribe medication solely on the basis that they have been previously prescribed is not acceptable.

    • Document discussion with patients that non-pharmacological therapy and non-opioid analgesics are preferred for chronic non-cancer pain, and that the potential benefit of long-term opioid treatment (LTOT) is modest and risk significant.

    • Advise patients that LTOT is not indicated for certain medical conditions including headache disorders, fibromyalgia and axial low back pain.

    • Always prescribe the lowest effective dosage of opioid medication. Doses >50 morphine milligram equivalents (MME) per day warrant careful reassessment and documentation. Doses >90 MME per day warrant substantive evidence of exceptional need and benefit. (This advice excludes treatment with methadone.)

    • When treating patients with acute pain conditions, prescribe only immediate release opioids in quantities that the patient will need before community follow-up will be resumed (3–7 days is often adequate).

    • When discharging patients from acute-care settings, or postoperatively, prescribe only the quantities of opioids, sedatives or stimulants that the patient will need before community follow-up will be resumed.

    • Base decisions to prescribe long-term psychoactive medications, including LTOT, on well-documented, comprehensive initial assessments and frequent (at least every 3 months) reassessments. These assessments and reassessments must include documented history and physical examination of the patient. There must also be documentation that the patient has been screened regularly for the presence or emergence of mental health and substance-use disorders and risk factors, and advised about safety-sensitive occupational risks, child care responsibilities and driving.

    • Document the offer of a take-home naloxone prescription to all patients who are at risk of respiratory depression as a consequence of receiving opioid medications.

    • Document having directed and regularly reminded patients for whom they are prescribing LTOT to abstain from alcohol and non-prescription sedatives.

    • Order at least annual random urine drug testing and/or random pill counts for all adult patients on long-term opioids, benzodiazepines, sedative hypnotics or stimulants.

    Further, physicians must not:

    • Prescribe benzodiazepines or sedative hypnotics to patients on LTOT, other than as a documented taper.

    • Prescribe combinations of opioids with benzodiazepines and/or sedative hypnotics.

    • Provide prescriptions allowing dispenses of opioids, sedatives and stimulants, which exceed a 3-month supply or 250 tablets, whichever is less.

    • Initiate treatment with drugs with a high-risk profile such as methadone and fentanyl without relevant training and experience.

    People with OUD are consistently marginalised and stigmatised from mainstream healthcare delivery systems,17 18 relying on expensive and fragmented care from emergency departments and walk-in clinics to address their acute and chronic medical conditions.19 20 Policies aimed at expanding the provision of OAT may therefore have positive spillover effects beyond the benefits of the treatment itself. An expanded pool of prescribers, mostly family medicine physicians, may enable access to primary care services; improving continuity of care and reducing the use of walk-in clinics and emergency departments.21–24

    Conversely, the CPSBC’s practice standard may have inadvertently restricted access to quality primary care for patients with chronic non-cancer pain.25 Across Canada and the USA there is growing anecdotal evidence that opioid-prescribing guidelines have negatively affected people who could benefit from opioids. Incorrect interpretations of the standard can result in aggressive weaning without consent, cause debilitating pain and serious withdrawal symptoms, strain patient–physician relationships and increase risk of overdose for patients who self-medicate for pain relief.26–32

    Our goal is to conduct a comprehensive investigation of the intended and unintended consequences of the changes to OAT prescribing and the CPSBC’s practice standard for patients and primary care physicians in BC. To meet this objective, we will:

    1. Examine the uptake of OAT prescribing and the adoption of the CPSBC’s practice standard among primary care physicians based on the following research questions (RQs):

      • 1.1 How do physicians who begin prescribing OAT between 2016 and 2018 differ from those who do not, and with providers already offering OAT prior to recent prescribing changes?

      • 1.2 To what extent do physicians vary in adopting the CPSBC’s practice standard?

      • 1.3 What are the effects of the practice standard on prescribing patterns of long-term opioid treatment (LTOT)?

      • 1.4 What physician characteristics are associated with ceasing prescription of controlled substances, and with terminating primary care for patients on LTOT?

    2. Determine the effects of the changes to opioid prescribing on primary care for patients based on the following RQs:

      • 2.1 What are the characteristics of patients who are newly prescribed OAT (2016–2018)?

      • 2.2 What are the characteristics of patients who experience rapid tapering and/or termination with their primary care provider following the release of the CPSBC’s practice standard?

      • 2.3 Do patients who begin OAT experience changes in primary care access, continuity or comprehensiveness?

      • 2.4 Do patients treated with LTOT experience changes in primary care access, continuity or comprehensiveness following implementation of the new practice standard?

    Methods and analysis

    This project is set to start September 2020 and end August 2024.

    Study setting

    BC is Canada’s most western province, and has a population of approximately 4.8 million residents. There are over 1200 OAT prescribers33 among approximately 6000 primary care physicians and addiction specialists.34 The province provides single-payer coverage of inpatient and outpatient health services through its Medical Services Plan (MSP). Residents excluded from the insurance programme include newly landed immigrants and people covered under federal insurance programmes including refugees, asylum seekers, military personnel and First Nations’ members (representing less than 4% of the population).35

    In 2012, fentanyl was first detected in the illicit drug supply, and 4% of the province’s 270 overdose deaths were fentanyl-related (illicit overdose events include indication of street drugs (controlled and illegal: heroin, cocaine, MDMA, methamphetamine, illicit fentanyl), and medications not prescribed to the decedent but obtained/purchased on the street, from unknown means, or where origin of drug not known). By 2018, fentanyl was detected in over 80% of the province’s 1535 drug overdose deaths.36 Dr Perry Kendall, the provincial medical health officer at the time, declared the opioid overdose epidemic a public health emergency in April 2016.37 The declaration enabled the BC Centre for Disease Control, Ministry of Health, Regional Health Authorities, BC Coroner Services and related stakeholders to quickly expand surveillance efforts and adopt new harm reduction programmes (eg, overdose prevention sites, take-home naloxone programmes), treatment and recovery interventions. Meanwhile, the Ministry of Health and the College of Physicians focused on supply side interventions to prevent overdose deaths and reduce harms.37 A variety of studies are currently underway to evaluate the effects of these interventions.

    Data sources and linkages

    This project combines five administrative health datasets under the authority of the BC Ministry of Health, with the CPSBC’s Physician Registry, BC Cancer Agency’s Cancer Registry and Vital Statistics’ Mortality data (table 1). Population Data BC, a multi-university data resource, will provide secure access to individual level, linked and de-identified data for our research purposes.

    View this table:
    Table 1

    Administrative datasets used to build cohort

    Patient and public involvement

    As part of the project’s development, we recruited a patient-partner with lived experience from the Patient Voices Network. Through a series of meetings, the patient-partner has informed the research aims of the project and outcome measures of interest. This patient-partner will remain an integral member of the team assisting in the interpretation of patient-level results and knowledge dissemination efforts. As this project uses secondary administrative data, no additional patients were recruited for the conduct of the study.

    Study design and study population

    This is a population-level, retrospective cohort study of all BC primary care physicians who prescribed any opioid medication for OUD or chronic non-cancer pain during the study period, and their patients. The study period is 1 January 2013–31 December 2018, with a 1-year wash-in period (1 January 2012–31 December 2012) to exclude patients who initiated LTOT before 2012 or whose pain type (ie, ‘chronic non-cancer’, ‘acute’, ‘cancer or palliative’, or ‘other’) cannot be confirmed (figure 1).

    Figure 1
    Figure 1

    Primary use decision tree. BC CDC, British Columbia Centre for Disease Control; OAT, opioid agonist treatment; Rx, prescription; Tx, treatment.

    The cohort will include all primary care physicians and their patients using the CPSBC’s Physician Registry and MSP billings. We will use PharmaNet data to identify all physicians (including specialists) who prescribed any opioid during the study period so we can determine the physician who initiated each opioid prescription. For all patients, MSP, National Ambulatory Care Reporting System (NACRS) and Discharge Abstract Database (DAD) will provide complete health service use history during the wash-in and study periods; indicate the types of services the prescription initiating physician provides and enable us to control for patient-level comorbidities. The Cancer Registry and PharmaNet will allow us to identify and exclude patients being treated for cancer and palliative care. DAD and mortality data will be used to identify patient subpopulations who are disproportionately harmed by the changes. All records in the MSP, NACRS, DAD and PharmaNet files include de-identified physician and referring practitioner numbers for referral pattern purposes.

    Operationalised measures

    We will operationalise a series of variables using our linked dataset (tables 2 and 3).

    View this table:
    Table 2

    Prescribing measures

    View this table:
    Table 3

    Quality of primary care measures

    Figure 2
    Figure 2

    Long-term opioid treatment (LTOT) concordant fill decision tree. *Excluding benzodiazepine and sedative hypnotic tapers. CNCP, chronic non-cancer pain; MME, morphine milligram equivalents; Rx, prescription.

    Analysis plan for each objective

    Data will be prepared using SAS (V.9.4) and statistical analyses will be performed using R (V.3.2.5). In general, missing and incomplete data will be excluded from analyses and the number of observations omitted from analyses due to missing data will be documented. If any imputation is used, the method and extent will be reported. We will report p values<0.05 and 95% confidence intervals. Results will be reported using appropriate Enhancing the QUAlity and Transparency Of health Research guidelines (eg, STrengthening the Reporting of OBservational studies in Epidemiology).

    Physicians’ uptake of OAT prescribing and the adoption of the CPSBC’s practice standard

    To examine the differences between primary care physicians who begin prescribing OAT to patients with OUD and those who do not (RQ1.1), we will restrict our analysis to OAT-naive physicians (ie, never prescribed OAT during the 3-year wash-in period: 2012–2014, inclusive), and saw OAT-naive patients with indication of an OUD between 2015 and 2018, inclusive. This will allow us to examine the cohort of physicians who were susceptible to OAT expansion interventions. We will use generalised estimating equations (GEE) for logistic regression to identify physician-level characteristics (eg, age, sex, years since training completion, geography, practice type and prescribing history) associated with OAT prescribing. These non-parametric models allow us to account for the repeated measures and hierarchical structure of our data by specifying joint distribution in their random effects terms, and are well suited to identifying population average differences.38 39 We will repeat the analysis to compare physician characteristics associated with new OAT prescribing (ie, first prescription post-2015) compared with early adopters (ie, primary care physicians prescribing OAT pre-2015).

    To understand the extent to which physicians vary in adopting the CPSBC’s practice standard (RQ1.2), we will restrict analysis to physicians who prescribed LTOT between June 2014 and May 2016 and again during the effective period of the CPSBC’s practice standard (June 2016–May 2018). We will use funnel plots to quantify the extent of deviation in prescription concordance (concordance prior to the implementation of the CPSBC’s practice standard (June 2014–May 2016) will be an artificial measure used to identify physician outliers susceptible to the effects of the practice standard post-implementation) (figure 2) at the physician-level before and after the implementation of the practice standard. Physicians whose observed proportion of non-concordant prescriptions remain above the upper 95% control limits of the expected proportion of non-concordant fills (given the number of prescriptions they prescribed, and controlling for patient and geographical differences) following the implementation of the standard, will be compared with primary care peers whose prescribing becomes concordant. We will use mixed-effects models to identify physician characteristics associated with non-concordance.

    To estimate the effects of the practice standard (implemented June 2016) on usual prescribers of care’s LTOT prescriptions (RQ1.3), we will use interrupted time series (ITS) analysis. ITS is a quasi-experimental study design that estimates the effects of service or policy interventions before and after implementation in contexts where randomised controlled trials are not feasible or ethical.40 41 The advantage of this before–after comparison with a single population is that selection bias and confounding due to between-group differences are limited; and within-group characteristics that change slowly over time (eg, physician characteristics), secular changes, random fluctuations from one time point to the next and regression to the mean are also controlled.42 The primary assumption for ITS is that without the intervention (here the introduction of the CPSBC’s practice standard), the observed pre-intervention outcome trends would continue unchanged into the post-intervention period. This assumption is supported by Crabtree et al’s43 work which found no change in trends of defined daily dose of opioids prescribed pre-implementation versus post-implementation of the practice standard. Additionally, to our knowledge there were no other interventions that would affect physician prescribing of LTOT implemented during our study period. Using this method, we will look at the effects of the practice standard on the number of LTOT prescriptions filled with: (a) a daily dose greater than 90 morphine milligram equivalents; (b) a benzodiazepine co-prescription; and (c) a supply that exceeds 3 months or 250 tablets (whichever is less). We will also look at (d) the number of physicians who terminate any controlled substance and (e) the number of LTOT prescriptions terminated. We selected these outcomes because the CPSBC’s primary aim for establishing the practice standard was to reduce prescription drug misuse. However, we would also like to quantify the extent of inappropriate treatment cessation, as suggested by mounting anecdotal evidence. Given the potential rarity of the outcomes, we will aggregate counts to eight quarterly periods pre-implementation and post-implementation (pre-intervention: June 2014–May 2016; post-intervention: June 2016–May 2018). We will use the two-sided Durbin-Watson test, plot of residuals, and autocorrelation plots to identify and adjust for autocorrelation and moving averages where necessary. We will test the hypothesis that the practice standard had no effect on these outcomes using ordinary least squares and segmented regression.

    Lastly, we will use GEE for logistic regression, and GEE for Poisson or negative binomial distribution to identify physician characteristics associated with any controlled substances cessation, and with counts of treatment termination between 2016 and 2018, respectively (RQ1.4). Negative binomial models are similar to Poisson models with the exception that the mean and variance of the count data do not have to be the same. These models include an additional parameter to handle overdispersion in the data, and are well suited for zero inflation and unobserved heterogeneity in the data.44 45 To determine which of the two model types is best suited to the data, we will use Pearson Χ2 dispersion statistics and residual plots.45 46

    Effects of the changes to opioid prescribing on primary care for patients

    We will use GEE for logistic regression to identify patient-level characteristics (eg, age, sex, geography, comorbidities, contraindications) associated with OAT initiation between 2015 and 2018 among OAT-naive patients (no OAT prescription filled between 2012 and 2014, inclusive) with indication of OUD (RQ 2.1). For patients on LTOT as of June 2016, we will use GEE to identify patient-level characteristics associated with inappropriate treatment termination (including rapid tapering) and regression models to estimate risk of mortality (RQ2.2).

    For RQ 2.3, we will use a variety of multiple linear (or linearised, where appropriate) regression models to measure the population-level association between expanded provision of OAT defined as the number of primary care physicians newly prescribing OAT as of 2015 (OAT prescribing naive 2012–2014) and changes to patients’ access, continuity and comprehensiveness of primary care (independently), while controlling for patient-level characteristics (eg, age, sex, geography, indication of OUD, comorbidities and contraindications). Models will include a lag for the potential delayed effects between expanded provision of OAT and our outcomes. Lags between the number of physicians prescribing OAT and primary care outcomes will be estimated using the weighted cumulative approach47–49 and informed with expert input from our knowledge users/stakeholders.50 51 We anticipate at least a 1-year lag between expanded OAT provision and changes in access, continuity and comprehensiveness of care given how these measures will be constructed.

    To measure the effects of the practice standard on patients’ access, continuity and comprehensiveness of primary care (RQ2.4), we will use controlled ITS analysis. This modified method allows us to account for time-varying confounders that may have influenced the delivery of primary care. Patients on LTOT just before the implementation of the practice standard will be matched with patients with diabetes not on opioids for pain management (eg, for neuropathic pain)52 or OAT whose pharmacotherapy is overseen by the same physician, and on age, sex, Charlson Comorbidity Index and treatment initiation month. Patients with diabetes were selected as our negative control group because:

    1. The majority of diabetes care is provided in primary care settings.

    2. There were no contemporaneous policies affecting diabetes pharmacotherapy prescribing behaviours during the study period.

    3. The opioid-prescribing standard is not expected to affect diabetes pharmacotherapy prescriptions.

    4. There are clear practice guidelines for physicians treating patients with diabetes, including frequency of physician oversight, which are similar to the guidelines for patients on LTOT.53

    5. Prescription disruptions (eg, termination, change in prescriber) are unusual for this patient population.

    We will use 24 monthly intervals pre-implementation and post-implementation for a total of 48 time periods between June 2014 and May 2018 (inclusive). At each interval, the outcome measure will be assessed for each patient dispensed a prescription during that month (LTOT for exposed, diabetes drug for matched controls). We will test our hypothesis of no change in access, continuity or comprehensiveness of care using ordinary least squares and segmented regression, by fitting the following regression model, per outcome:

    Embedded Image

    where j is the intervention, t is the study time in monthly intervals pre-intervention (negative time) and post-intervention (positive time), and k distinguishes between intervention and control group. Significant values for coefficients β6 and β7 will indicate an effect of the practice standard on patients’ access, continuity and comprehensiveness of primary care after accounting for level and trend changes among controls with diabetes.

    Ethics and dissemination

    This study has been approved by McGill University’s Institutional Review Board (#A11-M55-19A), and the University of British Columbia’s Research Ethics Board (#H19-03537).

    All data used in this project will be linked using de-identified personal health numbers or physician practice numbers. The data will be stored in Population Data BC’s Secure Research Environment, a central server for data storage and analysis, including encrypted backups, software and other services to ensure compliance with data access requirements.54 All members of the research team who will have access to the data have had the necessary tri-council privacy training and will complete privacy training provided by Population Data BC. Study results will be screened by data stewards prior to publication to ensure privacy and confidentiality requirements are maintained, there is no gross misuse of the data and data are appropriately referenced. Linked data will remain within Population Data BC’s Red Zone (terminals with no external connection) for up to 7 years after project completion before being destroyed by authorised data personnel.

    Discussion

    Our objective is to evaluate the effects of interventions aimed to improve access to OAT and limit overprescribing of opioid analgesics on physician prescribing behaviour and patient access to comprehensive primary care. Given the potential for patient harm, it is important to understand the effects of the interventions to prevent over/underprescribing of opioids, and to mitigate the effects of exacerbating or introducing new inequities in access to healthcare. We will use a variety of health service delivery outcomes and address limitations of existing research. Results from aim 1 will be useful for isolating the effects of the interventions on physician prescribing behaviour, and to identify physician characteristics associated with high-error prescribing in light of the new standard. Results from aim 2 will estimate the effects of the interventions on patients’ recent primary care experiences; and enable policy makers and physicians to identify potential subgroups harmed by the recent changes in prescribing. Together, these results will provide invaluable information on the effects of recent opioid-prescribing policies at both the patient and prescriber level, and equip policymakers and regulatory colleges with the much-needed information to understand the intended and unintended consequences of opioid-prescribing interventions to fine-tune their policies.

    The primary challenge of this project relates to the use of administrative health data to create quantifiable measures of access, continuity and comprehensiveness of care. In settings without comprehensive records of service use or where access to services and treatment is restricted, such methods can lead to under-ascertainment and estimation of key patient groups (eg, people with OUD) and undercounting primary care service use/provision.55 For example, reliance on administrative data to accurately identify patient health needs (eg, opioid use disorder, chronic pain) may lead to misclassification. However, by working with the comprehensive dataset described above in a context where the majority of health services are provided at no cost, and by applying modern surveillance methods, we minimise the potential for misclassification and capture as many patients as possible outside of a strict longitudinal cohort study.56 Of note, while most provinces (eg, Quebec) lack access to complete drug dispensation records for their populations, BC is unique in that it includes all outpatient prescription fills for all residents, irrespective of payer, within PharmaNet’s data file. These complete records, along with the BC Cancer Agency’s Cancer Registry enable the cross-validation of patient pain type identified using algorithms which rely on inpatient and outpatient records only, and allow us to work around the issue of inadequate specificity in International Classification of Diseases coding described elsewhere. Further, recognising the limitations of working with administrative data, we have intentionally recruited primary care physicians, members of the BCCSU and CPSBC, and a patient-partner living with chronic pain in BC to inform all stages of the project including the operationalisation of performance measures, their analysis and the interpretation of results. For all health system performance measures, we will also conduct extensive sensitivity analyses.

    Knowledge translation is integrated throughout the proposed study. Our team includes a primary care provider in active clinical practice, experts in substance-use treatment, policymakers and a patient with pertinent lived experience. The inclusion of diverse stakeholders within the project team has helped to confirm relevance of research aims and refine specific research questions. As the project progresses, we will meet with stakeholder groups to share interim findings and identify emergent policy areas where information from our project may inform decision-making. In addition to communicating results of research to stakeholders and policymakers, and through traditional academic channels (publications and conference presentations), we will also work to ensure end-of-grant research findings are accessible to people prescribed opioids in BC and to the public. We will disseminate results via a combination of open access peer-reviewed publications in high-impact journals, conferences, lay summaries, OpEds and ongoing meetings with our stakeholders/knowledge users.

    The overdose epidemic in BC is unique to that observed elsewhere in North America. In the USA, the opioid epidemic is described as a triple wave of overdose deaths starting with prescription opioids, followed by heroin, and more recently, fentanyl.57 58 Other regions in Canada demonstrate a similar epidemiological transition from prescription opioids to illicit substances. For BC, contamination of the illicit drug supply seems to have driven the epidemic since the beginning.59 Although the underlying drivers of BC’s overdose epidemic differ from those described in other parts of Canada and the USA, many of the responses to the epidemic parallel those observed elsewhere. The removal of the federal section 56 exemption from the Controlled Drugs and Substances Act for physicians to prescribe methadone, and subsequent provincial requirements for buprenorphine, along with coverage on public insurance programmes, expanded access to OAT for residents in Alberta, Ontario and Quebec.60 61 Similarly, while enforceable practice standards for the prescribing of opioid analgesics have not been implemented in other jurisdictions, provinces and states are struggling with purported effects of changes in prescribing guidelines. Given these response similarities, we expect findings will be of international relevance.

    References

    1. ↵
      1. Special Advisory Committee on the Epidemic of Opioid Overdoses
      . National report: apparent opioid-related deaths in Canada (January 2016 to December 2018). Ottawa: Public Health Agency of Canada, 2019.
    2. ↵
      1. Belzak L,
      2. Halverson J
      . The opioid crisis in Canada: a national perspective. Health Promot Chronic Dis Prev Can 2018;38:224–33.doi:10.24095/hpcdp.38.6.02pmid:http://www.ncbi.nlm.nih.gov/pubmed/29911818
      OpenUrlCrossRefPubMed
    3. ↵
      1. Ye X,
      2. Sutherland J,
      3. Henry B, et al
      . At-a-glance - impact of drug overdose-related deaths on life expectancy at birth in British Columbia. Health Promot Chronic Dis Prev Can 2018;38:248–51.doi:10.24095/hpcdp.38.6.05pmid:http://www.ncbi.nlm.nih.gov/pubmed/29911821
      OpenUrlCrossRefPubMed
    4. ↵
      1. Health P
      . Chronic disease prevention in Canada editorial. correction: at-a-glance - Impact of drug overdose-related deaths on life expectancy at birth in British Columbia. Health Promot Chronic Dis Prev Can 2019;39:64.
    5. ↵
      1. Vancouver Police Department
      . The opioid crisis: the need for treatment on demand. Vancouver Police Department, 2017.
    6. ↵
      1. Vancouver Coastal Health
      . Response to the opioid overdose crisis in Vancouver coastal health. Vancouver, BC: Vancouver Coastal Health, 2019.
    7. ↵
      1. Adams E, et al
      . A guideline for the clinical management: opioid use disorder. Vancouver, BC: BC Centre for Substance Use, 2017.
    8. ↵
      1. Dr. Bill
      . 5’addiction medicine fee codes’ available in BC, 2018.
    9. ↵
      1. Doctors of BC
      . New fees, updated guidelines support treatment of opioid use disorder. Doctors of BC, 2018.
    10. ↵
      1. British Columbia Centre on Substance Use
      . Provincial opioid addiction treatment support program. Vancouver, BC: British Columbia Centre on Substance Use, 2016.
    11. ↵
      1. College of Physicians and Surgeons of British Columbia
      . Safe prescribing of drugs with potential for Misuse/Diversion. College of Physicians and Surgeons BC, 2016.
    12. ↵
      1. Colvin LA,
      2. Bull F,
      3. Hales TG
      . Perioperative opioid analgesia-when is enough too much? A review of opioid-induced tolerance and hyperalgesia. Lancet 2019;393:1558–68.doi:10.1016/S0140-6736(19)30430-1pmid:http://www.ncbi.nlm.nih.gov/pubmed/30983591
      OpenUrlCrossRefPubMed
    13. ↵
      1. Busse JW, et al
      . Opioids for chronic noncancer pain: a systematic review and meta-analysis opioids for chronic noncancer pain. JAMA 2018;320:2448–60.
      OpenUrlCrossRefPubMed
    14. ↵
      1. Chou R,
      2. Deyo R,
      3. Devine B, et al
      . The effectiveness and risks of long-term opioid treatment of chronic pain. Evid Rep Technol Assess 2014:1–219.doi:10.23970/AHRQEPCERTA218pmid:http://www.ncbi.nlm.nih.gov/pubmed/30313000
    15. ↵
      1. Oetter HM
      . Clarification from the College of physicians and surgeons of bc on commentary about limitations of the CDC guideline for prescribing opioids. CMAJ 2017;189:E508. doi:10.1503/cmaj.732899pmid:http://www.ncbi.nlm.nih.gov/pubmed/28385898
      OpenUrlFREE Full Text
    16. ↵
      1. Health Professions Act
      . [RSBC 1996] Chapter 183. Victoria, BC, 2019.
    17. ↵
      1. Pollack HA,
      2. Khoshnood K,
      3. Blankenship KM, et al
      . The impact of needle exchange-based health services on emergency department use. J Gen Intern Med 2002;17:341–8.doi:10.1007/s11606-002-0037-2pmid:http://www.ncbi.nlm.nih.gov/pubmed/12047730
      OpenUrlPubMedWeb of Science
    18. ↵
      1. Wang L,
      2. Panagiotoglou D,
      3. Min JE, et al
      . Inability to access health and social services associated with mental health among people who inject drugs in a Canadian setting. Drug Alcohol Depend 2016;168:22–9.doi:10.1016/j.drugalcdep.2016.08.631pmid:http://www.ncbi.nlm.nih.gov/pubmed/27610937
      OpenUrlPubMed
    19. ↵
      1. French MT,
      2. McGeary KA,
      3. Chitwood DD, et al
      . Chronic illicit drug use, health services utilization and the cost of medical care. Soc Sci Med 2000;50:1703–13.doi:10.1016/S0277-9536(99)00411-6pmid:http://www.ncbi.nlm.nih.gov/pubmed/10798326
      OpenUrlCrossRefPubMedWeb of Science
    20. ↵
      1. Stein MD,
      2. Anderson B
      . Injection frequency mediates health service use among persons with a history of drug injection. Drug Alcohol Depend 2003;70:159–68.doi:10.1016/S0376-8716(02)00344-7pmid:http://www.ncbi.nlm.nih.gov/pubmed/12732409
      OpenUrlCrossRefPubMed
    21. ↵
      1. Schwarz R,
      2. Zelenev A,
      3. Bruce RD, et al
      . Retention on buprenorphine treatment reduces emergency department utilization, but not hospitalization, among treatment-seeking patients with opioid dependence. J Subst Abuse Treat 2012;43:451–7.doi:10.1016/j.jsat.2012.03.008pmid:http://www.ncbi.nlm.nih.gov/pubmed/22534003
      OpenUrlPubMed
    22. ↵
      1. Friedmann PD,
      2. Hendrickson JC,
      3. Gerstein DR, et al
      . Do mechanisms that link addiction treatment patients to primary care influence subsequent utilization of emergency and hospital care? Med Care 2006;44:8–15.doi:10.1097/01.mlr.0000188913.50489.77pmid:http://www.ncbi.nlm.nih.gov/pubmed/16365607
      OpenUrlCrossRefPubMedWeb of Science
    23. ↵
      1. Laine C,
      2. Lin Y-T,
      3. Hauck WW, et al
      . Availability of medical care services in drug treatment clinics associated with lower repeated emergency department use. Med Care 2005;43:985–95.doi:10.1097/01.mlr.0000178198.79329.89pmid:http://www.ncbi.nlm.nih.gov/pubmed/16166868
      OpenUrlCrossRefPubMedWeb of Science
    24. ↵
      1. Spithoff S,
      2. Kiran T,
      3. Khuu W, et al
      . Quality of primary care among individuals receiving treatment for opioid use disorder. Can Fam Physician 2019;65:343–51.pmid:http://www.ncbi.nlm.nih.gov/pubmed/31088874
      OpenUrlAbstract/FREE Full Text
    25. ↵
      1. Volkow ND,
      2. McLellan AT
      . Opioid abuse in chronic pain-misconceptions and mitigation strategies. N Engl J Med 2016;374:1253–63.doi:10.1056/NEJMra1507771pmid:http://www.ncbi.nlm.nih.gov/pubmed/27028915
      OpenUrlCrossRefPubMed
    26. ↵
      1. Kroenke K,
      2. Alford DP,
      3. Argoff C, et al
      . Challenges with implementing the centers for disease control and prevention opioid guideline: a consensus panel report. Pain Med 2019;20:724–35.doi:10.1093/pm/pny307pmid:http://www.ncbi.nlm.nih.gov/pubmed/30690556
      OpenUrlPubMed
    27. ↵
      1. Voon P,
      2. Callon C,
      3. Nguyen P, et al
      . Denial of prescription analgesia among people who inject drugs in a Canadian setting. Drug Alcohol Rev 2015;34:221–8.doi:10.1111/dar.12226pmid:http://www.ncbi.nlm.nih.gov/pubmed/25521168
      OpenUrlPubMed
    28. ↵
      1. Busse JW,
      2. Craigie S,
      3. Juurlink DN, et al
      . Guideline for opioid therapy and chronic noncancer pain. CMAJ 2017;189:E659–66.doi:10.1503/cmaj.170363pmid:http://www.ncbi.nlm.nih.gov/pubmed/28483845
      OpenUrlFREE Full Text
    29. ↵
      1. Antoniou T,
      2. Ala-Leppilampi K,
      3. Shearer D, et al
      . ‘Like being put on an ice floe and shoved away’: a qualitative study of the impacts of opioid-related policy changes on people who take opioids. Int J Drug Policy 2019;66:15–22.doi:10.1016/j.drugpo.2019.01.015pmid:http://www.ncbi.nlm.nih.gov/pubmed/30685650
      OpenUrlPubMed
    30. ↵
      1. Dassieu L,
      2. Kaboré J-L,
      3. Choinière M, et al
      . Chronic pain management among people who use drugs: a health policy challenge in the context of the opioid crisis. Int J Drug Policy 2019;71:150–6.doi:10.1016/j.drugpo.2019.03.023pmid:http://www.ncbi.nlm.nih.gov/pubmed/30975594
      OpenUrlPubMed
    31. ↵
      1. Darnall BD,
      2. Juurlink D,
      3. Kerns RD, et al
      . International Stakeholder community of pain experts and leaders call for an urgent action on forced opioid tapering. Pain Medicine 2019;20:429–33.doi:10.1093/pm/pny228
      OpenUrl
    32. ↵
      1. Lupick T
      . Opioid refugees: How the fentanyl crisis led to a backlash against doctors that’s leaving people in pain. In: The georgia straight. Vancouver, BC, 2018.
    33. ↵
      1. BC Centre for Disease Control
      . Overdose response indicators. In: Data & reports. Vancouver, BC: BC Centre for Disease Control, 2019.
    34. ↵
      1. Hedden L,
      2. Barer ML,
      3. McGrail K, et al
      . In British Columbia, the supply of primary care physicians Grew, but their rate of clinical activity declined. Health Aff 2017;36:1904–11.doi:10.1377/hlthaff.2017.0014pmid:http://www.ncbi.nlm.nih.gov/pubmed/29137511
      OpenUrlPubMed
    35. ↵
      1. Morrow RL,
      2. Bassett K,
      3. Wright JM, et al
      . Influence of opioid prescribing standards on drug use among patients with long-term opioid use: a longitudinal cohort study. CMAJ Open 2019;7:E484–91.doi:10.9778/cmajo.20190003pmid:http://www.ncbi.nlm.nih.gov/pubmed/31345786
      OpenUrlAbstract/FREE Full Text
    36. ↵
      1. BC Coroners Service
      . Illicit drug overdose deaths in bc: findings of coroners’ investigations. Ministry of Public Safety and Solicitor General, 2018.
    37. ↵
      1. BC Centre for Disease Control
      . Public health emergency in BC. In Stories. Vancouver, BC: BC Centre for Disease Control, 2017.
    38. ↵
      1. Gebregziabher M,
      2. Zhao Y,
      3. Dismuke CE, et al
      . Joint modeling of multiple longitudinal cost outcomes using multivariate generalized linear mixed models. Health Serv Outcomes Res Methodol 2013;13:39–57.doi:10.1007/s10742-012-0103-0pmid:http://www.ncbi.nlm.nih.gov/pubmed/30555270
      OpenUrlPubMed
    39. ↵
      1. Hubbard AE,
      2. Ahern J,
      3. Fleischer NL, et al
      . To GEE or not to GEE: comparing population average and mixed models for estimating the associations between neighborhood risk factors and health. Epidemiology 2010;21:467–74.doi:10.1097/EDE.0b013e3181caeb90pmid:http://www.ncbi.nlm.nih.gov/pubmed/20220526
      OpenUrlCrossRefPubMedWeb of Science
    40. ↵
      1. Cook TD
      . Quasi-experimentation: design analysis issues for field settings. Boston, MA: Houghton Mifflin Company, 1979.
    41. ↵
      1. Zhang F,
      2. Wagner AK,
      3. Ross-Degnan D
      . Simulation-based power calculation for designing interrupted time series analyses of health policy interventions. J Clin Epidemiol 2011;64:1252–61.doi:10.1016/j.jclinepi.2011.02.007pmid:http://www.ncbi.nlm.nih.gov/pubmed/21640554
      OpenUrlCrossRefPubMed
    42. ↵
      1. Lopez Bernal J,
      2. Cummins S,
      3. Gasparrini A
      . The use of controls in interrupted time series studies of public health interventions. Int J Epidemiol 2018;47:2082–93.doi:10.1093/ije/dyy135pmid:http://www.ncbi.nlm.nih.gov/pubmed/29982445
      OpenUrlCrossRefPubMed
    43. ↵
      1. Crabtree A,
      2. Rose C,
      3. Chong M, et al
      . Effects of the new prescribing standards in British Columbia on consumption of opioids and benzodiazepines and z drugs. Can Fam Physician 2019;65:e231–7.pmid:http://www.ncbi.nlm.nih.gov/pubmed/31088889
      OpenUrlAbstract/FREE Full Text
    44. ↵
      1. Pedan A
      . Analysis of count data using the SAS system. SAS Institute Statistics, Data analysis, and data mining.
    45. ↵
      1. Liu W,
      2. Cela J
      . Count Data in SAS. SAS Institute statistics and data analysis 2008.
    46. ↵
      1. Meyer J
      . Poisson or negative binomial? Using count model diagnostics to select model. In the analysis factor, 2018.
    47. ↵
      1. Sylvestre M-P,
      2. Abrahamowicz M
      . Flexible modeling of the cumulative effects of time-dependent exposures on the hazard. Stat Med 2009;28:3437–53.doi:10.1002/sim.3701pmid:http://www.ncbi.nlm.nih.gov/pubmed/19708037
      OpenUrlCrossRefPubMedWeb of Science
    48. ↵
      1. Danieli C,
      2. Abrahamowicz M
      . Competing risks modeling of cumulative effects of time-varying drug exposures. Stat Methods Med Res 2019;28:248–62.doi:10.1177/0962280217720947pmid:http://www.ncbi.nlm.nih.gov/pubmed/28882094
      OpenUrlPubMed
    49. ↵
      1. Xiao Y,
      2. Abrahamowicz M,
      3. Moodie EEM, et al
      . Flexible marginal structural models for estimating the cumulative effect of a time-dependent treatment on the hazard: reassessing the cardiovascular risks of didanosine treatment in the Swiss HIV cohort study. J Am Stat Assoc 2014;109:455–64.doi:10.1080/01621459.2013.872650
      OpenUrlCrossRef
    50. ↵
      1. Soyiri IN,
      2. Reidpath DD,
      3. Sarran C
      . Forecasting asthma-related hospital admissions in London using negative binomial models. Chron Respir Dis 2013;10:85–94.doi:10.1177/1479972313482847pmid:http://www.ncbi.nlm.nih.gov/pubmed/23620439
      OpenUrlCrossRefPubMed
    51. ↵
      1. Briët OJT,
      2. Amerasinghe PH,
      3. Vounatsou P
      . Generalized seasonal autoregressive integrated moving average models for count data with application to malaria time series with low case numbers. PLoS One 2013;8:e65761. doi:10.1371/journal.pone.0065761pmid:http://www.ncbi.nlm.nih.gov/pubmed/23785448
      OpenUrlCrossRefPubMed
    52. ↵
      1. Mu A,
      2. Weinberg E,
      3. Moulin DE, et al
      . Pharmacologic management of chronic neuropathic pain: review of the Canadian pain Society consensus statement. Can Fam Physician 2017;63:844–52.pmid:http://www.ncbi.nlm.nih.gov/pubmed/29138154
      OpenUrlAbstract/FREE Full Text
    53. ↵
      1. Diabetes Canada Clinical Practice Guidelines Expert Committee
      . Diabetes Canada 2018 clinical practice guidelines for the prevention and management of diabetes in Canada. Can J Diabetes 2018;42:S1–325.
      OpenUrlCrossRefPubMed
    54. ↵
      1. PopulationDataBC
      . Data access, 2018. Available: https://www.popdata.bc.ca/dataaccess [Accessed 25 Oct 2018].
    55. ↵
      1. Unick GJ,
      2. Rosenblum D,
      3. Mars S, et al
      . Intertwined epidemics: national demographic trends in hospitalizations for heroin- and opioid-related overdoses, 1993-2009. PLoS One 2013;8:e54496. doi:10.1371/journal.pone.0054496pmid:http://www.ncbi.nlm.nih.gov/pubmed/23405084
      OpenUrlCrossRefPubMed
    56. ↵
      1. Haut ER,
      2. Pronovost PJ,
      3. Schneider EB
      . Limitations of administrative databases. JAMA 2012;307:2589–90.doi:10.1001/jama.2012.6626pmid:http://www.ncbi.nlm.nih.gov/pubmed/22735421
      OpenUrlPubMed
    57. ↵
      1. Ciccarone D
      . The triple wave epidemic: supply and demand drivers of the US opioid overdose crisis. Int J Drug Policy 2019;71:183–8.doi:10.1016/j.drugpo.2019.01.010pmid:http://www.ncbi.nlm.nih.gov/pubmed/30718120
      OpenUrlPubMed
    58. ↵
      1. Schuchat A,
      2. Houry D,
      3. Guy GP
      . New data on opioid use and prescribing in the United StatesNew data on opioid use and prescribing, 2006-2015. JAMA 2017;318:425–6.
      OpenUrlPubMed
    59. ↵
      1. British Columbia Coroners Service
      . Illicit drug toxicity deaths in BC: January 1, 2010-June 30, 2020. Victoria, BC: BC Coroners Service, 2020.
    60. ↵
      1. Régie de l’assurance maladie du Québec
      . Prescription drug insurance, 2018.
    61. ↵
      1. Régie de l’assurance maladie du Québec
      . List of medications, 2018.
      1. Lacasse A,
      2. Ware MA,
      3. Dorais M, et al
      . Is the Quebec provincial administrative database a valid source for research on chronic non-cancer pain? Pharmacoepidemiol Drug Saf 2015;24:980–90.doi:10.1002/pds.3820pmid:http://www.ncbi.nlm.nih.gov/pubmed/26105572
      OpenUrlPubMed
      1. Clarke H,
      2. Soneji N,
      3. Ko DT, et al
      . Rates and risk factors for prolonged opioid use after major surgery: population based cohort study. BMJ 2014;348:g1251 doi:10.1136/bmj.g1251pmid:24519537
      OpenUrlAbstract/FREE Full Text
      1. College of Pharmacists of British Columbia
      . New safety measures for codeine liquid preparations coming soon. Victoria, BC, 2019.
      1. Liu Y,
      2. Logan JE,
      3. Paulozzi LJ, et al
      . Potential misuse and inappropriate prescription practices involving opioid analgesics. Am J Manag Care 2013;19:648–58.pmid:http://www.ncbi.nlm.nih.gov/pubmed/24304213
      OpenUrlPubMed
      1. Kantor ED,
      2. Rehm CD,
      3. Haas JS, et al
      . Trends in prescription drug use among adults in the United States from 1999-2012. JAMA 2015;314:1818–31.doi:10.1001/jama.2015.13766pmid:http://www.ncbi.nlm.nih.gov/pubmed/26529160
      OpenUrlCrossRefPubMed
      1. College of Pharmacists of British Columbia
      . Controlled prescription program. Vol. 2018.3. Vancouver, BC, 2018.
      1. Schultz SE,
      2. Glazier RH
      . Identification of physicians providing comprehensive primary care in Ontario: a retrospective analysis using linked administrative data. CMAJ Open 2017;5:E856–63.doi:10.9778/cmajo.20170083pmid:http://www.ncbi.nlm.nih.gov/pubmed/29259018
      OpenUrlAbstract/FREE Full Text

    Footnotes

    • Twitter @DrRitaMc

    • Contributors DP and RM are the co-principal Investigators of the study. DP, RM and MRL conceived the study and overall design, with critical revisions from TG, ECS, BF, AB, CJ and PK. All coauthors developed the analysis plan, with specific methodology input from TG and ECS. DP drafted the manuscript. All authors provided feedback on the manuscript, approved the final version to be published and agree to be accountable to all aspects of the work ensuring the project’s accuracy and integrity.

    • Funding This proposal is currently under review for funding. To date, this research has received no specific grant from any funding agency in the public, commercial or not-for-profit sectors. It is supported by new hire start-up funds from McGill University.

    • Competing interests None declared.

    • Patient and public involvement Patients and/or the public were involved in the design, or conduct, or reporting, or dissemination plans of this research. Refer to the Methods section for further details.

    • Patient consent for publication Not required.

    • Provenance and peer review Not commissioned; externally peer reviewed.