Skip to main content
Download PDF

Impact of a publicly-funded pharmacare program policy on benzodiazepine dispensing among children and youth: a population-based natural experiment

BMC Pediatrics volume 23, Article number: 519 (2023) Cite this article

Abstract

Background

In January 2018, the Government of Ontario, Canada, initiated a universal pharmacare program (OHIP+) for all individuals aged 24 years and younger. In April 2019, the program was amended to cover only children and youth without private insurance. Because benzodiazepines are commonly prescribed to children and youth despite their potential hazards, we examined whether changes in publicly-funded drug coverage influenced benzodiazepine dispensing trends in this demographic.

Methods

We conducted a population-based natural experiment study of benzodiazepine dispensing to children and youth in Ontario between January 2013 and March 2020. We used interventional autoregressive integrated moving average models to estimate the impact of OHIP + and its subsequent modification on these trends.

Results

The implementation of OHIP + was associated with an immediate increase in the monthly rate of benzodiazepine dispensing of 12.9 individuals per 100,000 population (95% confidence interval [CI]; 7.5 to 18.3 per 100,000). Benzodiazepine dispensing rates rose from 214.2 to 241.5 per 100,000 from December 2017 to March 2019, a 12.8% (95% CI 9.6–16.0%) increase. In stratified analyses, increases were most pronounced among females, children and youth living in the lowest income neighbourhoods and individuals aged 20 to 24. The April 2019 modification to OHIP + was not associated with changes in monthly benzodiazepine dispensing trends (0.39 individuals per 100,000; 95% CI -1.3 to 2.1 per 100,000). However, rates remained elevated relative to the period preceding OHIP + implementation.

Conclusions

Implementation of a publicly-funded pharmacare program resulted in more children and youth being prescribed benzodiazepines.

Peer Review reports

Introduction

The use of psychotropic medications in children and youth has increased worldwide [1]. However, in contrast to stimulants, antipsychotics, and antidepressants, little is known about patterns of benzodiazepine use in children and youth [1]. Studies from several jurisdictions in Europe and North America with publicly funded healthcare or combined private and public healthcare systems have found that the prevalence of benzodiazepine dispensing in children and youth ranges from 0.2 to 9.8%, with use increasing over time [1,2,3,4,5,6,7,8]. In a study from Manitoba, Canada, the incidence of benzodiazepine use among children and adolescents aged 17 and under increased from 2.2 to 3.95 per 1000 population between 1996/97 and 2011/12 [7]. In addition, findings from Ontario, Canada indicate that increased use in children and youth may be associated with harm. Specifically, between 2013 and 2020, rates of benzodiazepine-related toxicity healthcare encounters increased from 11.1 to 16.0 per 100,000 and 39.9 to 66.6 per 100,000 among children and adolescents under the age of 18 and youth aged 19 to 24, respectively [9].

While increased benzodiazepine use may reflect improved recognition and diagnosis of anxiety disorders in children and youth, there is a lack of research supporting the effectiveness and long-term safety of benzodiazepines in this population [10, 11]. Consequently, benzodiazepines are not approved as anxiolytics for children in North America or Europe, with clinical practice guidelines preferentially advocating behavioural therapies over medication [12, 13]. Moreover, although these drugs are sometimes used short-term in combination with antidepressants approved for anxiety disorder, [14] this practice is not recommended in individuals below the age of 18, and benzodiazepines are not endorsed as a therapeutic option for any psychiatric disorder in this population [12, 13]. In light of the increasing off-label use of benzodiazepines in children and youth and concerns about the possibilities of dependence, diversion and nonmedical misuse, [15,16,17,18,19] research examining patterns of pediatric benzodiazepine use and the influence of specific interventions on the use of these drugs is needed to inform clinical practice and policy.

Although several studies have examined the impact of policies and guidance intended to curb benzodiazepine use among adults, [20,21,22] less is known about the influence on benzodiazepine use of interventions intended to increase access to prescription medications more generally. In January 2018, the Ontario government implemented a publicly-funded pharmacare program known as OHIP + to provide all Ontarians aged 24 and younger prescription medications listed on the Ontario Drug Benefit formulary at no cost [23]. Coverage was automatic, with no deductibles or copayments. The program was subsequently modified in March 2019 to only cover medications for children and youth without private insurance [24]. Because OHIP + eliminated out-of-pocket costs and disparities in drug insurance coverage, we postulated that dispensing of benzodiazepines to children and youth may have increased immediately following implementation. Additional changes in benzodiazepine use may have also occurred with the subsequent modification of OHIP + in April 2019, which restricted universal drug coverage to children and youth without private insurance. Therefore, the implementation and modification of OHIP + offered a natural experiment for evaluating the influence of universal pharmacare on the use of benzodiazepines in children and youth and evaluating the differential impact of these policy changes on various populations of children and youth. Accordingly, we studied the impact of OHIP + and its modification on benzodiazepine dispensing among the entire population of eligible individuals aged 24 and younger in Ontario, home to approximately 40% of Canadian children and youth [25]. We speculated that the implementation of OHIP + would be associated with an increase in benzodiazepine use and that a more gradual change would be observed following its modification to cover only those children without private insurance.

Methods

Setting

We conducted a population-based study of all residents in Ontario aged 0 to 24 years between January 1, 2013, and March 31, 2020.

Data sources

We used Ontario’s administrative health databases. These datasets were linked using unique encoded identifiers and analyzed at ICES (formerly known as the Institute for Clinical Evaluative Sciences). We identified prescriptions for benzodiazepines using the Narcotics Monitoring System database, which contains comprehensive records of all prescriptions for controlled substances dispensed from community pharmacies in Ontario, regardless of payer. Because we were focusing on the use of benzodiazepines for mental health conditions, we excluded prescriptions for clobazam, which is used primarily for seizure disorders. We used the Registered Persons Database, a registry of all individuals eligible for the publicly-funded Ontario Health Insurance Plan, to ascertain demographic characteristics for all children and youth dispensed benzodiazepines over the study period. We determined the proportions of individuals diagnosed with anxiety, mood, or seizure disorders within 30 days prior to or on the prescription dispensing date using outpatient physician claims data from the Ontario Health Insurance Plan database, emergency department data from the Canadian Institute for Health Information National Ambulatory Care Reporting System, and hospitalization data from the Canadian Institute for Health Information Discharge Abstract Database and Ontario Mental Health Reporting System (see supplemental Table 1 for diagnostic codes). The OHIP database includes all physician claims (primary care and specialist) for medical services covered under the provincial health insurance plan. We obtained prescriber specialty information using the ICES Corporate Provider Database. The use of data in this project was authorized under Sect. 45 of Ontario’s Personal Health Information Protection Act, which does not require review by a Research Ethics Board, or individual patient consent. ICES approved the study and all analyses of raw data were conducted at ICES. No other administrative permissions were required.

Study Population and Outcomes

For each month in the study period, we defined our study population as all Ontario residents aged 0 to 24 who were alive on the first day of the month. Our primary outcome was the monthly rate of benzodiazepine use per 100,000 children and youth, defined as the number of individuals dispensed a benzodiazepine in a given month divided by the population of children and youth aged 0 to 24 for that period. To determine whether characteristics of children and youth dispensed a benzodiazepine changed following the implementation and modification of OHIP+, we compared demographic characteristics and prescriber type among individuals receiving benzodiazepines during the period immediately preceding the implementation of OHIP+ (January 1, 2017, to December 31, 2017), during OHIP+ (January 2018 to March 2019) and the period in which OHIP + was available only for children and youth without private insurance (April 2019 to March 2020).

Statistical analysis

We used interventional autoregressive integrated moving average (ARIMA) models to examine the impact of the implementation of OHIP + and its subsequent modification on benzodiazepine dispensing rates among children and youth [26, 27]. We used the Dickey-Fuller test to determine the stationarity of the time series and applied first order and seasonal differencing to arrive at a stationary series if needed [27, 28]. The Dickey-Fuller statistic tests the null hypothesis a unit root is present in a time series and that the time series is non-stationary (i.e. does not exhibit constant variance over time) against the alternative hypothesis of stationarity. A p-value of greater than 0.05 and a failure to reject the null hypothesis indicates the presence of a unit root and that differencing is required to render the time series stationary. We used the autocorrelation function and partial autocorrelation function to identify autoregressive and/or moving average components in each time series and correct for any autocorrelation remaining after differencing, and selected the best models using goodness of fit tests [26, 27]. We used residual plots and the Portmanteau statistic to confirm that residuals from specified ARIMA models were a white noise process [29]. The null hypothesis of the Portmanteau test is that there is no residual autocorrelation and that the residuals are a white noise process. A p-value of greater than 0.05 is therefore required to infer the lack of residual autocorrelation before the ARIMA model can be used for inference and forecasting. Finally, once the ARIMA models were specified, we used a step intervention function to test for a change in the rate of benzodiazepine dispensing during the period in which OHIP + provided universal coverage of prescription medication to Ontario children and youth (January 1, 2018 to March 31, 2019), and compared predicted values of dispensing during this period with the observed rates [26, 27]. The ARIMA model also included a ramp intervention function to determine if benzodiazepine dispensing rates changed following modification to the program in April 2019, providing coverage only for children and youth with no private insurance. To explore heterogeneity in the impact of the OHIP + program, we stratified our analyses by sex, age category (0 to 9 years, 10 to 14 years, 15 to 19 years, 20 to 24 years), neighbourhood income quintile and urban versus rural residence, defined on the first day of the month of interest. We used standardized differences to compare demographic characteristics between individuals receiving a benzodiazepine during the pre-OHIP+, OHIP + and modified OHIP + periods, with differences greater than 0.1 representing an imbalance between the two groups [30]. The analyses used SAS Enterprise Guide, version 7.1 (SAS Institute Inc., Cary, NC, USA).

Results

There were 218,299 children and youth aged 0 to 24 that received a benzodiazepine prescription between January 1, 2013 and March 31, 2020. The median age of benzodiazepine-treated children and youth was 20 years (interquartile range 17 to 22 years), and the majority (n = 136,091; 62.3%) were female (Table 1). General practitioners accounted for most benzodiazepine prescribing (n = 135,218; 61.9%), with fewer children and youth receiving these drugs from pediatricians (n = 10,470; 4.8%) and psychiatrists (n = 31,863; 14.6%) (Table 1). The majority (n = 148,892; 68.3%) of individuals were dispensed a short-term (i.e. 14 days or less) supply of benzodiazepines, with nearly 1 in 5 receiving a supply exceeding 30 days in duration. The number of children with a diagnosis of anxiety, mood and seizure disorder in the 30 days prior to receiving a benzodiazepine prescription was 102,397 (46.9%), 26,599 (12.2%) and 8,817 (4.0%), respectively. Relative to individuals dispensed less than a 30 days’ supply, those dispensed more than a 30 days’ supply of were more likely to be between the ages of 20 and 24 (55.9% vs. 51.1%; standardized difference [SD] = 0.10), receive their prescription from a psychiatrist (30.9% vs. 11.4%; SD = 0.49) and have a mood disorder in the 30 days prior to receiving a benzodiazepine (18.1% vs. 11.0%; SD = 0.34) (supplemental Table 2). Overall, the demographic characteristics of children and youth receiving a benzodiazepine did not change appreciably between the pre-OHIP+, OHIP + and modified-OHIP + periods (Table 1).

Table 1 Demographic and clinical characteristics of individuals aged 0 to 24 dispensed a benzodiazepine, January 2013 to March 2020*
Full size table

The most frequently prescribed benzodiazepines during the study period were clonazepam and lorazepam, with variation observed according to age and neighbourhood income quintile. Specifically, lorazepam and clonazepam represented 56.5% and 28.0% of prescriptions dispensed to children under the age of 12, while respective estimates for those 12 years and older were 42.7% and 42.5% (supplemental Table 3). This pattern was consistent across the pre-OHIP+, OHIP + and modified OHIP + periods (supplemental Table 3). In terms of socioeconomic status, clonazepam was the most frequently dispensed benzodiazpine to children and youth in the lowest income neighbourhoods, representing 45.9% of benzodiazepine prescriptions during the study period (supplemental table 4a). The proportion of prescriptions that were clonazepam declined as neighbourhood income quintile increased, representing 38.9% of benzodiazepine prescriptions in the highest income neighbourhoods. The opposite pattern was observed for lorazepam, which accounted for 47.5% and 38.0% of benzodiazepine prescriptions in the highest and lowest income neighbourhoods, respectively (supplemental table 4a). Other benzodiazepines for which income gradients were observed included alprazolam, which was more common in the highest- relative to the lowest-income neighbourhoods (4.3% versus 3.2%), and diazepam, which was more frequently dispensed among individuals in the lowest relative to the highest income neighbourhoods (7.3% versus 4.7%) (supplemental table 4a). These patterns were consistent across the three study periods (supplemental tables 4b, 4c and 4d).

Change in benzodiazepine dispensing rates following implementation of OHIP+

The average monthly percent change in benzodiazepine dispensing in the period preceding OHIP + was 0.3% (0.3–0.4%). We observed a modest relative percent increase in benzodiazepine dispensing to children and youth following the implementation of OHIP+, with rates increasing 12.8% (95% confidence interval [CI] 9.6–16.0%) between December 2017 and March 2019 (214.2 vs. 241.5 per 100,000 population, respectively) (Table 2). In stratified analyses, the increase was more pronounced in females (255.5 vs. 295.5 individuals per 100,000 population) than males (174.9 vs. 190.3 individuals per 100,000 population), with relative percent increases of 15.7% (95% CI 11.4–20.0%) and 8.8% (95% CI 4.0–13.7%), respectively, between December 2017 and March 2019. In addition, individuals in the lowest income neighbourhoods had the greatest relative percent increase in benzodiazepine dispensing, with rates increasing 15.5% (95% CI 8.4–22.9%) (223.0 vs. 257.6 individuals per 100,000 population) and 19.6% (95% CI 11.9–27.6%) (215.8 vs. 258.2 individuals per 100,000 population) in the lowest and second-lowest income quintile areas, respectively. Conversely, rates increased only 5.2% (95% CI -1.3–11.8%) in the highest income neighbourhoods (224.8 vs. 236.4 individuals per 100,000 population). Furthermore, individuals aged 20 to 24 (581.6 vs. 675.9 individuals per 100,000 population) and 15 to 19 (303.1 vs. 326.2 individuals per 100,000 population) had the greatest increase in benzodiazepine dispensing following the implementation of OHIP+, corresponding to relative percent increases of 16.2% (95% CI 12.0–20.5%) and 7.6% (95% CI 1.9–13.9%), respectively, between December 2017 and March 2019.

Table 2 Changes in benzodiazepine dispensing to children and youth following the introduction of the OHIP+ pharmacare program in January 2018
Full size table

Following ARIMA modelling, implementation of OHIP + was associated with a significant immediate increase in the monthly rate of benzodiazepine dispensing of 12.9 per 100,000 population (95% CI; 7.5 to 18.3 per 100,000). (Table 2; Fig. 1). The monthly rate of increase in benzodiazepine dispensing prior to OHIP + was 0.77 per 100,000 population (95% CI; 0.59 to 0.94). We observed similar results in stratified analyses (supplemental Fig. 1 to 4), with the largest immediate increases among females (16.5 per 100,000 population; 95% CI 8.3 to 24.7), individuals aged 15 to 19 (16.5 per 100,000 population; 95% CI 7.0 to 26.0) and 20 to 24 (41.1 per 100,000 population; 95% CI 26.0 to 56.1), and children and youth living in rural areas (14.6 per 100,000 population; 95% CI 4.6 to 24.6) (Table 2).

Fig. 1
figure 1

Impact of OHIP + implementation (January 2018) and modification (April 2019) on monthly rates of benzodiazepine dispensing among Ontario residents between the ages of 0 and 24, January 2013 to March 2020

Full size image

Observed benzodiazepine dispensing rates in the first 12 months of OHIP + were higher than those predicted by the ARIMA model in the absence of this program (Table 3; Fig. 2). The largest difference was observed in October 2018, corresponding to an additional 1,119 children and youth dispensed a benzodiazepine (Table 3).

Table 3 Projected and Actual Benzodiazepine Dispensing in Children and Youth, January 2018 to December 2018
Full size table
Fig. 2
figure 2

Actual versus predicted rates of benzodiazepine dispensing during OHIP+ (January 2018 to December 2018)

Full size image

Change in benzodiazepine dispensing rates following modification of OHIP+

Following modification of the OHIP + program, there was a small relative percent decrease in benzodiazepine dispensing, with rates decreasing 4.0% (95% CI -6.7% to -1.3%) between March 2019 and March 2020 (241.5 vs. 231.9 individuals per 100,000 population) (Table 4). However, rates remained elevated relative to the period preceding OHIP + implementation. Following ARIMA modelling, there was no significant change in monthly benzodiazepine dispensing trends following the April 2019 modification to OHIP+ (0.39 per 100,000; 95% CI -1.3 to 2.1 per 100,000) (Table 3; Fig. 1), with similar results obtained in stratified analyses (Table 4, supplemental Fig. 1 to 4).

Table 4 Changes in benzodiazepine dispensing to children and youth following the April 2019 modification of the OHIP+ pharmacare program restricting coverage only to individuals without private insurance
Full size table

Discussion

In our population-based study, we observed an increase in the rate of benzodiazepine dispensing among children and youth following the implementation of a publicly-funded pharmacare program covering all prescription drug costs for Ontario residents between the ages of 0 and 24. Increased benzodiazepine dispensing was most pronounced among females, children and youth living in lower-income neighbourhoods, and individuals aged 15 to 19 and 20 to 24 years. Modifying the program to maintain coverage only for individuals lacking private insurance was not associated with a significant change in benzodiazepine dispensing trends.

Our findings add to previous research. As in other jurisdictions, we observed increased benzodiazepine dispensing over time, with greater use in females and older youth relative to younger children [4,5,6,7]. These findings reflect known patterns in the diagnosis of anxiety disorder in children and youth, with a higher prevalence among females and a typical age of diagnosis between early adolescence and young adulthood [31,32,33]. However, our study is the first to estimate the influence of a newly-implemented pharmacare program on benzodiazepine dispensing. Our finding of greater than expected benzodiazepine dispensing in the 12 months following the implementation of OHIP + is especially salient, highlighting the relative increase in benzodiazepene use attributable to this program. Moreover, we were able to study the differential impact of such a program on benzodiazepine use in specific sub-populations of children and youth. Therefore, our work extends the study of trends in benzodiazepine use among children and youth to a setting with no financial barriers to drug therapy.

Our study has important policy and practice implications for children and youth with mental health conditions. The largest increase in dispensing was in 20 to 24-year-olds, the group most commonly prescribed these drugs prior to the implementation of OHIP+. Although considered in the definition of youth, individuals aged to 20 to 24 can arguably be regarded as young adults for whom benzodiazepine medications need to be dispensed with all the cautions given to older adults. This age group includes individuals enrolled at post-secondary institutions with school insurance plans that typically cover only a portion of drug costs or place caps on annual prescription drug benefits. However, the increased dispensing of these drugs to adolescents aged 15 to 19 years immediately following the implementation of OHIP + is of concern, given the absence of evidence supporting long-term effectiveness for anxiety or sleep disorders and the potential for serious adverse effects, including decreased alertness, dependence, withdrawal, diversion and injury [34,35,36]. Moreover, there is concern that increased use of these drugs in children and youth could lead to a public health crisis similar to that observed with opioids [4, 37, 38]. For these reasons, benzodiazepines are not recommended for treating pediatric anxiety or sleep disorders and are considered second-line or short-term adjunctive therapies for adults [11, 12, 39,40,41,42]. Furthermore, although 4.0% of children and youth dispensed a benzodiazepine had a diagnosis of seizure disorder, the finding that nearly 1 in 5 individuals were dispensed a supply of benzodiazepines exceeding 30 days suggests long-term use for mental health conditions and symptoms could also be occurring. This practice is inconsistent with clinical guidelines advising short-term therapy of mental health symptoms until evidence-based therapies have been accessed [12]. Notably, most benzodiazepine prescriptions were written by general practitioners who may be less familiar with guidelines or have limited access to colleagues with expertise in treating pediatric mental health conditions. Increased use of telemedical health programs facilitating access to pediatric mental health specialists represents a promising mechanism for supporting general practitioners caring for children and youth with mental health conditions and optimizing psychotropic prescribing in this population [43, 44].

In addition, we observed a greater relative increase in benzodiazepine dispensing among children and youth in lower-income neighbourhoods following the implementation of OHIP+. This finding may reflect inequities created by long wait times for publicly-funded mental health and behavioural programs and out-of-pocket costs associated with accessing privately-funded interventions not covered by Ontario’s publicly-funded healthcare program, such as clinical psychologists and social workers [45,46,47]. A similar phenomenon may account for the greater immediate increase in benzodiazepine dispensing to children and youth in rural communities relative to urban centres following the implementation of OHIP+, where delays are incurred because of local service gaps and the travel costs associated with accessing care [45]. In contrast, there are no delays associated with accessing medications, and prescription drugs were available at no cost for all children and youth immediately following the implementation of OHIP+. It is therefore possible that systemic disparities in access to non-drug therapies create conditions that promote disproportionate reliance on benzodiazepines for symptom management among lower-income children and youth and individuals living in rural communities.

Strengths of our study include complete benzodiazepine dispensing data for all children and youth in Ontario, regardless of insurance status. However, our study has some limitations. First, we cannot reliably ascertain the appropriateness of benzodiazepine use or indication for therapy. This limitation is common to all population-based studies using claims-based data to study benzodiazepine use. Second, we could not determine whether benzodiazepines were dispensed for as-needed or ongoing use, or the proportion of children and youth progressing to long-term treatment with these drugs. Similarly, it is likely that a portion of children and youth dispensed a short supply of benzodiazepines were receiving these drugs for procedural anxiety, a practice which is generally well tolerated and does not involve continued use [48]. Third, it is possible that the period following the modification of OHIP + was too brief to detect a significant change. Fourth, misclassification of individuals with certain mental health diagnoses is possible because we used previously derived administrative algorithms that were not validated for our study [49]. Fifth, we did not longitudinally examine whether OHIP + changed the prevalence of long-term (e.g., 180 days) benzodiazepine use. Finally, our study was conducted in a single Canadian province, potentially limiting the generalizability of our findings. However, ours is the first population-based natural experiment study estimating the impact of a newly introduced pharmacare program on benzodiazepine dispensing among children and youth. Our findings suggest that additional measures may be needed in jurisdictions contemplating similar interventions to mitigate increases in benzodiazepine use, including improved access to behavioural and mental health interventions and measures to facilitate collaboration between primary care providers and pediatric mental health specialists.

In conclusion, implementing a publicly-funded pharmacare program was associated with increased benzodiazepine dispensing to children and youth and greater than expected use of these drugs. Further, we found disproportionate increases in the rate of benzodiazepine dispensing among females, adolescents and young adults and low-income children and youth. Although these trends align with the patterns of anxiety diagnosis in adolescents and young adults, socioeconomic differences in dispensing suggest that removing financial barriers to benzodiazepines may have increased the use of these drugs in a population with structural barriers to accessing guideline-recommended non-pharmacologic treatments. Moreover, past research in adults has found an association between increasing use of benzodiazepines, long-term use, misuse and benzodiazepine-related toxicity with low-income status or residence in low-income neighbourhoods [9, 50, 51, 52]. Further research is needed to determine if these associations extend to children and youth and to examine the role of alternative healthcare delivery models, such as telehealth and programs training primary care physicians to provide specialist care, is warranted to determine whether such interventions can improve prescribing. Moreover, interventions to support increased and equitable access to non-pharmacologic mental health expertise are required to minimize reliance on benzodiazepines for the treatment of mental health conditions in children and youth.

Data Availability

The datasets generated and/or analysed during the current study are not publicly available due to data sharing agreements which prohibit ICES from making the data set publicly available. The data set from this study is held securely in coded form at ICES, and access may be granted to those who meet pre-specified criteria for confidential access, available at www.ices.on.ca/DAS. The full data set creation plan and underlying analytic code are available from the authors upon request (contact Dr. Tony Antoniou), understanding that the programs may rely upon coding templates or macros that are unique to ICES.

Abbreviations

ARIMA:

Autoregressive integrated moving average

Cis:

Confidence intervals

IQR:

Interquartile range

SD:

Standardized difference

References

  1. Piovani D, Clavenna A, Bonati M. Prescription prevalence of psychotropic drugs in children and adolescents: an analysis of international data. Eur J Clin Pharmacol. 2019;75:1333–46.

    Article  PubMed  Google Scholar 

  2. O’Sullivan K, Reulbach U, Boland F, Motterlini N, Kelly D, Bennett K, et al. Benzodiazepine prescribing in children under 15 years of age receiving free medical care on the General Medical Services scheme in Ireland. BMJ Open. 2015;5:e007070.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Bénard-Laribière A, Noize P, Pambrun E, Bazin F, Verdoux H, Tournier M, et al. Trends in incident use of benzodiazepines and Z-drugs in France from 2006 to 2012: a population-based study. Pharmacoepidemiol Drug Saf. 2017;26:162–9.

    Article  PubMed  Google Scholar 

  4. Sidorchuk A, Isomura K, Molero Y, Hellner C, Lichtenstein P, Chang Z, Franck J, de la Fernández L, Mataix-Cols D. Benzodiazepine prescribing for children, adolescents, and young adults from 2006 through 2013: a total population register-linkage study. PLoS Med. 2018;15:e1002635.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Hugtenburg JG, Heerdink ER, Egberts AC. Increased psychotropic drug consumption by children in the Netherlands during 1995–2001 is caused by increased use of methylphenidate by boys. Eur J Clin Pharmacol. 2004;60:377–9.

    Article  CAS  PubMed  Google Scholar 

  6. Huerta C, Abbing-Karahagopian V, Requena G, Oliva B, Alvarez Y, Gardarsdottir H, et al. Exposure to benzodiazepines (anxiolytics, hypnotics and related drugs) in seven european electronic healthcare databases: a cross-national descriptive study from the PROTECT-EU Project. Pharmacoepidemiol Drug Saf. 2016;25(Suppl 1):56–65.

    Article  CAS  PubMed  Google Scholar 

  7. Alessi-Severini S, Bolton JM, Enns MW, Dahl M, Collins DM, Chateau D, et al. Use of benzodiazepines and related drugs in Manitoba: a population-based study. CMAJ Open. 2014;2:E208–16.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Nielsen ES, Rasmussen L, Hellfritzsch M, Thomsen PH, Nørgaard M, Laursen T. Trends in off-label prescribing of sedatives, hypnotics and antidepressants among children and adolescents - a danish, Nationwide Register-Based study. Basic Clin Pharmacol Toxicol. 2017;120:360–7.

    Article  CAS  PubMed  Google Scholar 

  9. Campbell TJ, Men S, Shearer D, Ebejer T, Joosse M, Quercia J et al. The epidemiology of benzodiazepine-related toxicity in Ontario, Canada: a population-based descriptive study. Can J Public Health 2023 Jun 15:1–11. https://doi.org/10.17269/s41997-023-00784-3. Epub ahead of print.

  10. Patel DR, Feucht C, Brown K, Ramsay J. Pharmacological treatment of anxiety disorders in children and adolescents: a review for practitioners. Transl Pediatr. 2018;7:23–35.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Witek MW, Rojas V, Alonso C, Minami H, Silva RR. Review of benzodiazepine use in children and adolescents. Psychiatr Q. 2005;76:283–96.

    Article  PubMed  Google Scholar 

  12. Walter HJ, Bukstein OG, Abright AR, Keable H, Ramtekkar U, Ripperger-Suhler J, et al. Clinical practice Guideline for the Assessment and treatment of children and adolescents with anxiety Disorders. J Am Acad Child Adolesc Psychiatry. 2020;59:1107–24.

    Article  PubMed  Google Scholar 

  13. Katzman MA, Bleau P, Blier P, Chokka P, Kjernisted K, Van Ameringen M, et al. Canadian clinical practice guidelines for the management of anxiety, posttraumatic stress and obsessive-compulsive disorders. BMC Psychiatry. 2014;14(Suppl 1):1. https://doi.org/10.1186/1471-244X-14-S1-S1.

    Article  Google Scholar 

  14. Bushnell GA, Compton SN, Dusetzina SB, Gaynes BN, Brookhart MA, Walkup JT, et al. Treating Pediatric anxiety: initial use of SSRIs and other Antianxiety prescription medications. J Clin Psychiatry. 2018;79:16m11415.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Braüner JV, Johansen LM, Roesbjerg T, Pagsberg AK. Off-label prescription of Psychopharmacological Drugs in Child and Adolescent Psychiatry. J Clin Psychopharmacol. 2016;36:500–7.

    Article  PubMed  Google Scholar 

  16. Smogur M, Onesanu A, Plessen KJ, Eap CB, Ansermot N. Psychotropic drug prescription in children and adolescents: approved Medications in European Countries and the United States. J Child Adolesc Psychopharmacol. 2022;32:80–8.

    Article  PubMed  Google Scholar 

  17. McCabe SE, West BT. Medical and nonmedical use of prescription benzodiazepine anxiolytics among U.S. high school seniors. Addict Behav. 2014;39:959–64.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Boyd CJ, Austic E, Epstein-Ngo Q, Veliz PT, McCabe SE. A prospective study of adolescents’ nonmedical use of anxiolytic and sleep medication. Psychol Addict Behav. 2015;29:184–91.

    Article  PubMed  Google Scholar 

  19. Votaw VR, Geyer R, Rieselbach MM, McHugh RK. The epidemiology of benzodiazepine misuse: a systematic review. Drug Alcohol Depend. 2019;200:95–114.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Ross-Degnan D, Simoni-Wastila L, Brown JS, Gao X, Mah C, Cosler LE, et al. A controlled study of the effects of state surveillance on indicators of problematic and non-problematic benzodiazepine use in a Medicaid population. Int J Psychiatry Med. 2004;34:103–23.

    Article  PubMed  Google Scholar 

  21. Bijlsma MJ, Hak E, Bos J, De Jong-van den Berg LT, Janssen F. Assessing the effect of a guideline change on drug use prevalence by including the birth cohort dimension: the case of benzodiazepines. Pharmacoepidemiol Drug Saf. 2013;22:933–41.

    Article  CAS  PubMed  Google Scholar 

  22. Pearson SA, Soumerai S, Mah C, Zhang F, Simoni-Wastila L, Salzman C, et al. Racial disparities in access after regulatory surveillance of benzodiazepines. Arch Intern Med. 2006;166:572–9.

    Article  PubMed  Google Scholar 

  23. Ministry of Health and Long-Term Care. OHIP+: Children and Youth Pharmacare Key Facts for Prescribers. Available at: https://www.afhto.ca/wp-content/uploads/ohip_prescriber-002.pdf. Accessed October 2, 2022.

  24. Ministry of Health and Long-Term Care. OHIP + Redesign Update – FAQs for Patients. Available at: https://www.health.gov.on.ca/en/pro/programs/drugs/opdp_eo/notices/fq_pat_exec_office_20190311.pdf. Accessed October2, 2022.

  25. Statistics Canada. Focus on Geography Series., 2016 Census. Available at: https://www12.statcan.gc.ca/census-recensement/2016/as-sa/fogs-spg/Index-eng.cfm. Accessed October 4, 2022.

  26. Helfenstein U. The use of transfer function models, intervention analysis and related time series methods in epidemiology. Int J Epidemiol. 1991;20:808–15.

    Article  CAS  PubMed  Google Scholar 

  27. Schaffer AL, Dobbins TA, Pearson SA. Interrupted time series analysis using autoregressive integrated moving average (ARIMA) models: a guide for evaluating large-scale health interventions. BMC Med Res Methodol. 2021;21:58.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Dickey DA, Fuller WA. Distribution of the estimators for autoregressive time series with a unit root. J Am Stat Assoc. 1979;74:427–31.

    Google Scholar 

  29. Ljung GM, Box GEP. On a measure of lack of fit in time series models. Biometrika. 1978;65:297–303.

    Article  Google Scholar 

  30. Austin PC, Grootendorst P, Anderson GM. A comparison of the ability of different propensity score models to balance measured variables between treated and untreated subjects: a Monte Carlo study. Stat Med. 2007;26:734–53.

    Article  PubMed  Google Scholar 

  31. Solmi M, Radua J, Olivola M, Croce E, Soardo L, Salazar de Pablo G, et al. Age at onset of mental disorders worldwide: large-scale meta-analysis of 192 epidemiological studies. Mol Psychiatry. 2022;27:281–95.

    Article  CAS  PubMed  Google Scholar 

  32. Lijster JM, Dierckx B, Utens EM, Verhulst FC, Zieldorff C, Dieleman GC, et al. The age of onset of anxiety Disorders. Can J Psychiatry. 2017;62:237–46.

    Article  PubMed  Google Scholar 

  33. Merikangas KR, He JP, Burstein M, Swanson SA, Avenevoli S, Cui L, et al. Lifetime prevalence of mental disorders in U.S. adolescents: results from the National Comorbidity Survey replication–adolescent supplement (NCS-A). J Am Acad Child Adolesc Psychiatry. 2010;49:980–9.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Edinoff AN, Nix CA, Hollier J, Sagrera CE, Delacroix BM, Abubakar T, et al. Benzodiazepines: uses, Dangers, and clinical considerations. Neurol Int. 2021;13:594–607.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Bushnell GA, Gerhard T, Crystal S, Olfson M. Benzodiazepine Treatment and Fracture Risk in Young persons with anxiety Disorders. Pediatrics. 2020;146:e20193478.

    Article  PubMed  Google Scholar 

  36. Guina J, Merrill B, Benzodiazepines I. Upping the care on Downers: the evidence of risks, benefits and Alternatives. J Clin Med. 2018;7:17.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Chen J. Are Benzodiazepines the New Opioids? Yale Medicine Reviews, December 11, 2019. Available at: https://www.yalemedicine.org/news/benzodiazepine-epidemic. Accessed October 5, 2022.

  38. Friedrich JM, Sun C, Geng X, Calello DP, Gillam M, Medeiros KL, et al. Child and adolescent benzodiazepine exposure and overdose in the United States: 16 years of poison center data. Clin Toxicol (Phila). 2020;58:725–31.

    Article  CAS  PubMed  Google Scholar 

  39. Kansagra S. Sleep Disorders in adolescents. Pediatrics. 2020;145(Suppl 2):204–S209.

    Article  Google Scholar 

  40. Ophoff D, Slaats MA, Boudewyns A, Glazemakers I, Van Hoorenbeeck K, Verhulst SL. Sleep disorders during childhood: a practical review. Eur J Pediatr. 2018;177:641–8.

    Article  CAS  PubMed  Google Scholar 

  41. Ipser JC, Stein DJ, Hawkridge S, Hoppe L. Pharmacotherapy for anxiety disorders in children and adolescents. Cochrane Database Syst Rev. 2009;3:CD005170.

    Google Scholar 

  42. Sateia MJ, Buysse DJ, Krystal AD, Neubauer DN, Heald JL. Clinical practice guideline for the pharmacologic treatment of chronic insomnia in adults: an American Academy of Sleep Medicine clinical practice guideline. J Clin Sleep Med. 2017;13:307–49.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Myers KM, Valentine JM, Melzer SM. Child and adolescent telepsychiatry: utilization and satisfaction. Telemed J E Health. 2008;14:131–7.

    Article  PubMed  Google Scholar 

  44. Gloff NE, LeNoue SR, Novins DK, Myers K. Telemental health for children and adolescents. Int Rev Psychiatry. 2015;27:513–24.

    Article  PubMed  Google Scholar 

  45. Children’s Mental Health Ontario. Kids can’t wait: 2020 report on wait lists and wait times for child and youth mental health care in Ontario. Available at: https://cmho.org/wp-content/uploads/CMHO-Report-WaitTimes-2020.pdf. Accessed October 5, 2022.

  46. Gratzer D, Goldbloom D. Making evidence-based psychotherapy more accessible in Canada. Can J Psychiatry. 2016;61:618–23.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Mental Health Commission of Canada. Options for improving access to counselling, psychotherapy and psychological services for mental health problems and illnesses. Available at: https://www.mentalhealthcommission.ca/wp-content/uploads/drupal/2017-07/Options_for_improving_access_to_counselling_psychotherapy_and_psychological_services_eng.pdf. Accessed October 5, 2022.

  48. Kuang H, Johnson JA, Mulqueen JM, Bloch MH. The efficacy of benzodiazepines as acute anxiolytics in children: a meta-analysis. Depress Anxiety. 2017;34:888–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. The Mental Health and Addictions Scorecard and Evaluation Framework Research Team. The mental health of children and youth in Ontario: 2017 scorecard. Toronto, ON: Institute for Clinical Evaluative Sciences; 2017.

    Google Scholar 

  50. Fride Tvete I, Bjørner T, Skomedal T. Risk factors for excessive benzodiazepine use in a working age population: a nationwide 5-year survey in Norway. Scand J Prim Health Care. 2015;33:252–9.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Sjöstedt C, Ohlsson H, Li X, Sundquist K. Socio-demographic factors and long-term use of benzodiazepines in patients with depression, anxiety or insomnia. Psychiatry Res. 2017;249:221–5.

    Article  PubMed  Google Scholar 

  52. Taipale H, Särkilä H, Tanskanen A, Kurko T, Taiminen T, Tiihonen J, et al. Incidence of and characteristics Associated with Long-term Benzodiazepine Use in Finland. JAMA Netw Open. 2020;3:e2019029.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

This study was supported by ICES, which is funded by an annual grant from the Ontario Ministry of Health (MOH) and the Ministry of Long-Term Care (MLTC). Parts of this material are based on data and/or information compiled and provided by CIHI and the Ontario Ministry of Health. The analyses, conclusions, opinions and statements expressed herein are solely those of the authors and do not reflect those of the funding or data sources; no endorsement is intended or should be inferred. This document used data adapted from the Statistics Canada Postal CodeOM Conversion File, which is based on data licensed from Canada Post Corporation, and/or data adapted from the Ontario Ministry of Health Postal Code Conversion File, which contains data copied under license from ©Canada Post Corporation and Statistics Canada. We thank IQVIA Solutions Canada Inc. for use of their Drug Information File.

Funding

This study was funded by the Canadian Institutes of Health Research (funding reference number 166149).

Author information

Authors and Affiliations

  1. Li Ka Shing Knowledge Institute, St. Michael’s Hospital, Toronto, ON, Canada

    Tony Antoniou, Muhammad Mamdani & Tara Gomes

  2. ICES, Toronto, ON, Canada

    Tony Antoniou, Daniel McCormack, Sophie Kitchen, William Gardner, Yona Lunsky, Mina Tadrous, Muhammad Mamdani, David N. Juurlink & Tara Gomes

  3. Department of Family and Community Medicine, University of Toronto, Toronto, ON, Canada

    Tony Antoniou

  4. Department of Family and Community Medicine, St. Michael’s Hospital, Toronto, ON, Canada

    Tony Antoniou

  5. Children’s Hospital of Eastern Ontario Research Institute, Ottawa, ON, Canada

    Kathleen Pajer

  6. Department of Psychiatry, University of Ottawa, Ottawa, ON, Canada

    Kathleen Pajer

  7. School of Epidemiology and Public Health, University of Ottawa, Ottawa, ON, Canada

    William Gardner

  8. Azrieli Adult Neurodevelopmental Centre, Centre for Addiction and Mental Health, Toronto, Canada

    Yona Lunsky

  9. Department of Psychiatry, University of Toronto, Toronto, ON, Canada

    Yona Lunsky

  10. Autism Research Centre, Bloorview Research Institute, Holland Bloorview Kids Rehabilitation Hospital, Toronto, Canada

    Melanie Penner

  11. Department of Pediatrics, University of Toronto, Toronto Ontario, Canada

    Melanie Penner & David N. Juurlink

  12. Leslie Dan Faculty of Pharmacy, University of Toronto, Toronto, ON, Canada

    Mina Tadrous, Muhammad Mamdani & Tara Gomes

  13. Li Ka Shing Centre for Healthcare Analytics Research & Training, Unity Health, Toronto, ON, Canada

    Muhammad Mamdani

  14. Temerty Faculty of Medicine, University of Toronto, Toronto, ON, Canada

    Muhammad Mamdani

  15. Institute of Health Policy, Management, and Evaluation (Mamdani), University of Toronto, Toronto, ON, Canada

    Muhammad Mamdani & Tara Gomes

  16. Department of Medicine, University of Toronto, Toronto, ON, Canada

    David N. Juurlink

Authors
  1. Tony Antoniou
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Daniel McCormack
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Sophie Kitchen
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Kathleen Pajer
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. William Gardner
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Yona Lunsky
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Melanie Penner
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Mina Tadrous
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Muhammad Mamdani
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. David N. Juurlink
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Tara Gomes
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

TA, DM, SK, TG, MP, MT, KP, WG, MM, YL, and DJ conceptualized and designed the study, and were involved in the interpretation of the data. DM and SK conducted the analyses. TA drafted the initial manuscript. All authors critically reviewed the manuscript for intellectual content and revised the manuscript. All authors approved the final manuscript as submitted and agree to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.

Corresponding author

Correspondence to Tony Antoniou.

Ethics declarations

Ethics approval and consent to participate

ICES is an independent, non-profit research institute whose legal status under Ontario’s health information privacy law allows it to collect and analyze health care and demographic data, without consent, for health system evaluation and improvement. ICES is a prescribed entity under Sect. 45 of Ontario’s Personal Health Information Protection Act. Section 45 authorizes ICES to collect personal health information, without consent, for the purpose of analysis or compiling statistical information with respect to the management of, evaluation or monitoring of, the allocation of resources to or planning for all or part of the health system. Projects conducted under Sect. 45, by definition, are exempt from and do not require review by a Research Ethics Board. This is confirmed in a letter from the REB of Sunnybrook Health Sciences Centre, ICES’ Research Ethics Board of Record (See Appendix A). This project was conducted under Sect. 45, and approved by ICES’ Privacy and Compliance Office.

All methods were carried out in accordance with relevant guidelines and regulations.

Consent for publication

Not applicable.

Competing interests

Melanie Penner has received consulting fees for unrelated work from Addis & Associates/Roche and from the Government of Nova Scotia. Mina Tadrous has received consulting fees for unrelated work from Green Shield Canada and the Canadian Agency for Drugs and Technologies in Health. Tara Gomes has received funding from the Ontario MOH for unrelated work. Tony Antoniou has no competing interests. Daniel McCormack has no competing interests. Sophie Kitchen has no competing interests. Kathleen Pajer has no competing interests. William Gardner has no competing interests. Muhammad Mamdani has no competing interests. Yona Lunsky has no competing interests. David Juurlink has no competing interests.

Additional information

Publisher’s Note

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

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Supplementary Material 2

Supplementary Material 3

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Antoniou, T., McCormack, D., Kitchen, S. et al. Impact of a publicly-funded pharmacare program policy on benzodiazepine dispensing among children and youth: a population-based natural experiment. BMC Pediatr 23, 519 (2023). https://doi.org/10.1186/s12887-023-04331-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12887-023-04331-4

Keywords

BMC Pediatrics

ISSN: 1471-2431

Contact us