- Open Access
Analysis of changes in trends in the consumption rates of benzodiazepines and benzodiazepine-related drugs
Journal of Pharmaceutical Policy and Practicevolume 11, Article number: 1 (2018)
To analyse trends in the rates of consumption of benzodiazepine (BZD) anxiolytics, BZD hypnotics and non-BZD hypnotics and the association with contextual factors.
Descriptive time series study. Units of analysis were monthly dose per inhabitant per day (DID) and dose per medical card per day(DCD) of benzodiazepine(BZD anxiolytics(BZD-A), BZD hypnotics(BZD-H) and non-BZD hypnotics(Non-BZD-H) between January 2006-December 2015. We analysed 6 primary healthcare districts(PHD) and used defined daily doses (DDDs) to calculate the monthly DIDs(overall and by ATC group). Trends and monthly percentage change (MPC) were analysed through joinpoint regression.
The annual DID increased by 26% overall, the trend was different across ATC groups. Consumption in BZD-A and BZD-H increased (27.1%,61.9%), consumption in Non-BZD-H decreased by 35%. There was high variability in DCD across the PHD, with an overall increase of 10.2%(5.7%-22.9%). By ATC, DCD increased by 10.4% in BZD-A(4.2%-22.2%) and by 44.2% in BZD-H(33.2%-76.5%). The overall DCD in the Non-BZD-H decreased by 42.1%(19.7%-50.8%). We found an initial upward trend in consumption of BZD-A until April/2008(monthly percentage change –MPC- +0.5%), followed by a slightly slower increase (+0.1%). No changes in trend were detected in BZD-H. In Non-BZD-H, we observed an upward trend until February/2013(+0.1%), followed by a sharp decrease until August/2013(−6.3%), and finally a slight decrease(−0.3%).
BZD consumption has increased in the last decade, with variability across areas. The changes in trends do not coincide with the financial crisis, introduction of prescriptions by active ingredient, electronic prescriptions or copayment. The only decrease in the Non-BZD-H may be linked to an intervention.
Benzodiazepines (BZDs) are used to treat insomnia, anxiety and chronic back pain. BZDs have significant adverse effects, including memory loss, increased risk of accidents and falls, and dependency [1,2,3]. Although guidelines recommend psychological treatments as first-line therapy before pharmacological treatment , and that BZDs only be used in the short term (2-4 weeks) when symptoms are incapacitating [5,6,7], numerous studies have shown that BZDs are being over-prescribed for prolonged periods in many countries [8,9,10].
For several years, various initiatives have promoted the rational use of diagnostic and therapeutic resources to avoid unnecessary interventions and/or those with potential risks for patients in specific situations. These include the Choosing Wisely initiative (2009)  and the NICE ‘do not do’ recommendations (2011) . In 2013, a new initiative, called the Commitment to Quality of Spanish Scientific Societies, was launched in Spain. As part of this programme, the Spanish Society of Family and Community Medicine (SEMFYC) published a series of recommendations that included “Do not prescribe benzodiazepines (or non-benzodiazepine hypnotics) for long-term use in patients with insomnia” .
BZDs were first used in clinical practice in the early 1960s , and since then have been one of the most-prescribed drug groups in most developed countries . In Spain, BZD consumption has undergone sustained growth, from a defined daily dose (DDD) per 1000 inhabitants per day (DID) of 32.7 in 1992, to 56.7 in 2000, and 89.3 in 2012 [16, 17]. This 54.7% increase between 2000 and 2012 is beyond comparison with trends in any other countries in our setting with the exception of Portugal, where BZD consumption increased by 24% between 2003 and 2010 . In all other countries, BZD consumption rates remained stable or even decreased, with a more marked reduction in anxiolytics than hypnotics [18,19,20,21,22].
According to the latest Annual Report of the Spanish National Health Service (2015), the 15 most-prescribed active ingredients include 3 benzodiazepines: 2 anxiolytics (lorazepam and alprazolam, 16.4 and 11.8 million bottles, respectively) and a sedative hypnotic (lormetazepam, 9.3 million bottles). Prescriptions for these 3 active ingredients alone had a retail value of 79.3 million euros during 2015 .
As evidence-based clinical guidelines  on BZD use are not being applied in clinical practice, we saw the need for an analysis of how prescribing trends have changed over recent years. The aim of this study is therefore to study trends in BZD consumption rates in the province of Seville and its primary healthcare districts (PHDs) and healthcare management areas (HMAs) between 2006 and 2015. We set out to quantify those trends, delimit the resulting trend periods and correlate changes to relevant events during the period, such as the start of the financial crisis, the introduction of prescriptions by active ingredient, electronic prescriptions and copayment, or safety alerts issued by regulatory agencies (FDA, AEMPS) [24, 25].
Descriptive time series study.
Units of analysis
Total and disaggregated monthly dose per inhabitant per day (DID) and dose per medical card per day (DCD) of medicines derived from and related to benzodiazepine, in 3 therapeutic groups (BZD anxiolytics – N05BA (group BZD-A); BZD hypnotics – N05CD (group BZD-H); and non-BZD hypnotics – N05CF (group Non-BZD-H)) between January 2006 and December 2015.
Four primary healthcare districts (Seville (DS), Aljarafe (DA), Seville North (DSN) and Seville South (DSS)) and 2 healthcare management areas (Osuna (HMAO) and Seville South (HMASS)) in the province of Seville, Spain. Those organizations serve 8.2% of the total population of the region of Andalusia and 1,5% of the whole of Spain. All inhabitants are covered by the Public Health System.
Defined daily dose (DDD)
The average maintenance dose for a drug used for its main indication via one route of administration in adults . The DDD is not the same as the recommended or prescribed dose. DDDs are established by the World Health Organization (WHO) and available on the website of the WHO Collaborating Centre for Drug Statistics.
Dose per card per day (DCD)
Defined as the DDD per 1000 medical cards per day. DCDs can be used to compare PHDs and HMAs with different population densities.
Dose per inhabitant per day (DID)
Defined as the DDD per 1000 inhabitants per day. The DID provides information about a population’s exposure to a particular drug or group of drugs. The annual DID can be calculated based on the number of packs/bottles dispensed using the following formula:
U = units (packs/bottles); PF = number of pharmaceutical forms per unit; Q = amount of active ingredient in each pharmaceutical form. The monthly DID can be calculated using the same formula, replacing the 365 days of the year with the number of days in the month studied. Data regarding prescriptions and treatment costs per day in the primary healthcare district of Seville were obtained from the Andalusian Public Health Service.
DDDs and DCDs by ATC group were obtained based on data in the dispensing database of the Andalusian Health Service, which records all prescriptions dispensed by community pharmacies throughout Andalusia. We calculated DIDs using data from Seville’s provincial census, which gives the number of inhabitants per year on 1 January of that year .
The database includes the prescriptions detailed by age, sex, active drug principle, dose, dose schedule and duration of treatment. The quality of the database is acceptable, since the payment to community pharmacies depends in part on this.
We used the DDDs to calculate the 120 monthly DIDs (overall and broken down by ATC group ) in the 10-year study period. We created graphs showing monthly trends and the differences between consumption rates in the first and last years, or in the first and last months of each period. We used joinpoint regression to explore trends using Joinpoint Trend Analysis Software. This statistical analysis program was developed by the Surveillance, Epidemiology, and End Results (SEER) Program of the National Cancer Institute to analyse trend joinpoints. We set the maximum number of joinpoints (changes in trend) for each analysis to 5.
Consumption rates over time
The annual DID increased by 26% overall. However, the trend was not equal across all ATC therapeutic groups: consumption rates of drugs in groups BZD-A and BZD-H increased by 27.1% and 61.9%, respectively, but consumption of drugs in the Non-BZD-H group decreased by 35%. The results were similar when we compared the variation in overall DID between January 2006 and 2015. There was an overall increase of 21.8% (29.3% in BZD- A, 64.5% in BZD-H). However, DID in the Non-BZD-H group decreased by 35.5% (Table 1 and Fig. 1).
Between December 2015 and January 2006, there was high variability in DCD by therapeutic group across the different PHDs and HMAs. Overall DCD increased by 10.2%, but the increase ranged from 5.7% in DS to 22.9% in HMAO. By ATC group, overall DCD increased by 10.4% in group BZD-A (ranging from 4.2% in DS to 22.2% in HMAO), while in group BZD-H overall DCD increased by 44.2% (33.2% in DSN to 76.5% in HMAO). Finally, the overall DCD in the Non-BZD-H group decreased by 42.1% (19.7% in DSN to 50.8% in DSS) (Fig. 2).
When we analysed the 120 monthly DIDs by therapeutic group, we found an upward trend in group BZD-A, with 2 distinct segments separated by a joinpoint. Between January 2006 and April 2008, the monthly percentage change (MPC) was 0.5%, while between April 2008 and December 2015 the MPC was slightly lower at 0.1%. No changes in trend were observed in the BZD-H group, where there was a sustained increase (0.4%) in DID. In the Non-BZD-H group we found 3 trend segments separated by 2 joinpoints. DID increased initially (MPC 0.1%) between January 2006 and February 2013. There was then a sharp decrease (−6.3%) between February 2013 and August 2013. Finally, there was a slight decrease (−0.3%) between August 2013 and December 2015 (Table 2).
When we examined all groups together, the overall trend was upwards, with a single joinpoint in March 2009. There was an increase (MPC 0.4%) in the first segment, followed by moderate growth (0.1%) until December 2015. When we analysed the overall DCD in the different PHDs and HMAs, we found 2 segments divided by 1 joinpoint. In the first segment, DCD increased with an MPC of 0.2% (in DS and DA) or 0.3% (DSN, DSS and HMAO). In the second segment, DID decreased with an MPC of −0.1% (DA, DSN and HMAO) or −0.2% (DS and DSS), with the joinpoint between June 2011 (DS) and March 2012 (DA) (Table 2).
We observed an overall decrease in the Non-BZD-H group, and therefore carried out a detailed analysis of the monthly DCD trends in this therapeutic group by PHD/HMA. In the DS, DA and DSN districts we found 3 segments separated by 2 joinpoints. In the first segment, there was a slight upward trend with an MPC of 0.1% (DS and DA) or 0.2% (DSN) until January (DS and DSN) or February (DA) 2013. In the second segment, there was a decrease in DCD with an MPC of −1.9% (DSN), −6.2% (DA) and −8.4% (DS) until August 2013 (DS), September 2013 (DA) or June 2014 (DSN). There was a smaller decrease of −0.2% (DS and DSN) or −0.3% (DA) in the final segment. We should highlight that the first segment for DSS and HMAO is almost flat, i.e. the MPC is close to or equal to zero and not statistically significant. The first joinpoint for DSS is in September 2011, followed by a second, slight downward segment with an MPC of −0.6% until February 2013. The second joinpoint in DSS, in February 2013, marks the start of a sharp decrease in DCD (−10.7%). This increase lasts just 4 months, until June 2013. The decrease then becomes less extreme (−1.2%) until August 2014. The last segment, up to May 2016, is relatively flat and not statistically significant. In HMAO, after an initial flat segment, the second segment beginning in March 2013 shows a marked decrease (−15.2%), although it lasts for just 3 months. Its third segment starts in 2013 with a slight decrease (−0.7%) until October 2015. There is a fourth segment between October 2015 and May 2016, with an MPC of 0.2% (Table 2, Fig. 2).
This study shows that BZD consumption in the province of Seville has increased by more than 20% in 10 years. This finding is consistent with the results of a national study by the Spanish Ministry of Health, which found an annual increase of 2.5-3% between 2000 and 2011 with some spikes . We conducted our study in the province of Seville, home to 23.1% of the population of the region of Andalusia . This is a large sample, and our study is the first to break down consumption rates by primary healthcare district and healthcare management area, enabling a detailed analysis of trends. The increase detected in the rate of consumption of BZDs cannot be down solely to an ageing or growing population, as the increase in consumption is almost 5 times the increase in the number of inhabitants .
This study has some limitations which should be considered when interpreting the results. Firstly, we used 2 indicators for the analysis of trends: DID for provincial trends and DCD for the healthcare districts and healthcare management areas. Although they are very similar, DCD may have been underestimated as there may be users with medical cards who are not included in the census. When we analysed the increase by level of care, we found that the increase in primary care (overall DCD in PHDs/HMAs) was 10.2%, while the overall provincial increase (DID, primary and hospital care) was 21.8%. We can therefore assume that approximately 11.6% of the increase is attributable to hospital prescriptions. Because of the significant impact of such prescriptions (>50% of the increase), we recommend that studies of prescriptions in this level of care be conducted.
When we compare the overall DIDs for the province of Seville with those for Spain as a whole for the last available period (2006-12) , it appears that the growth has been slower (a provincial increase of 20.2% vs. a national increase of 26.4%), although the initial rates were higher (DID of 74.1 vs. 70.6). Both DIDs hit around 89 in 2012, and continued to increase in all subgroups. A study in Asturias had almost identical results . Because we do not have the updated national data, we do not know if the Non-BZD-H group followed the same downward trend nationally as found in Seville from 2013 onwards. BDZs are one of the most commonly prescribed classes of psychotropic medications in developed countries  but rates of consumption in the last years remains high and stable. For example, in 2012, there were approximately 85 million BDZs prescriptions written in the United States to outpatients, which was not significantly changed from the 90 million written in 2001 In the Dutch population aged 55–64, overall b BDZs use remained stable from 1992 to 2012, with a high proportion of long-term users, despite the effort to reduce BDZs use and the renewal of the guidelines .
The relative importance of drugs in group BZD-A (>80% of prescriptions) means that the increase in the overall rate (26%) is almost exactly the same as the increase in that group alone (27.1%). This group cancels out the increase in the BZD-H group and decrease in the Non-BZD-H group. This can be seen in the trend graphs, where the trend for group BZD-A is very similar to the overall trend, although the latter is slightly modified by the decrease in Non-BZD-H consumption.
One noteworthy result is the significant decrease in consumption of drugs in the Non-BZD-H group. Zolpidem is the most-prescribed drug in this group, accounting for 90% of all prescriptions. In January 2013, the FDA  warned of the risk of using zolpidem at high doses. As a result, primary care pharmacists designed a training intervention for primary care doctors in the province of Seville. The aim was to raise awareness of the translated FDA alert, providing a list of patients affected and proposed interventions. Given the effectiveness of this initiative in the province of Seville , interventions of this sort should be extended to all healthcare professionals in both primary and hospital care settings, whose actions largely determine the success or failure of any healthcare intervention.
When we began our analysis, we expected the start of the global financial crisis (2009)  to cause an increase in the rates of consumption of these drugs. However, we observed no changes in trend consistent with this hypothesis. This result is consistent with those of studies in Paterna  and Asturias . In fact, the only change observed was in March 2009, when the MPC of overall DID decreased from 0.4 to 0.1%. We also hypothesised that the introduction of copayment of medicines (July 2012)  would lead to a decrease in consumption. However, as in a previous study in Murcia , there were small decreases in the first few months but the rates bounced back quickly in the following months, with no overall change in trend. Electronic prescriptions and prescriptions by active ingredient were both introduced in 2005 and are now in widespread use. We expected the introduction of these measures to cause a change in BZD consumption. However, we observed no changes in trend consistent with this hypothesis.
The use of BZDs has increased in the last decade, with wide variability in the rates of BZD consumption within the province of Seville. The variability between neighbouring and homogenous small areas makes us hypothesize that the increase in use of BZDs is not only due to the increased prevalence of insomnia and due to other associated conditions.
The changes in trends do not coincide with the start of the financial crisis or the introduction of prescriptions by active ingredient, electronic prescriptions or co-payment. The only decrease in the Non-BZD-H group could be linked to a training intervention conducted with doctors, run in response to an alert issued by the US Food and Drug Administration (FDA) . The multiple causes of the increase in BZD consumption need to be studied in detail, and further training interventions of this sort must be run.
BZD consumption has increased in the last decade, with variability across geographical areas. The changes in trends do not coincide with the financial crisis, introduction of prescriptions by active ingredient, electronic prescriptions or copayment. The only decrease in the Non-BZD-H may be linked to an educational intervention for physicians.
Healthcare management areas
Primary healthcare districts
Barker MJ, Greenwood KM, Jackson M, Crowe SF. Persistence of cognitive effects after withdrawal from long-term benzodiazepine use: a meta-analysis. Arch Clin Neuropsychol. 2004 Apr;19(3):437–54.
Lader MH. Limitations on the use of benzodiazepines in anxiety and insomnia: are they justified? Eur Neuropsychopharmacol. 1999 Dec;9(Suppl 6):S399–405.
Centro Andaluz de Información del Medicamento (CADIME). Benzodiazepines: risks and strategies for their withdrawal. Granada; Available from: https://doi.org/10.11119/BTA2014-29-02. Accessed 21 Dec 2017.
NICE. Guidance on the use of zaleplon, zolpidem and zopiclone for the short-term management of insomnia. National Institute Centre for Excellence: London; 2004.
Straand J, Rokstad K. General practitioners’ prescribing patterns of benzodiazepine hypnotics: are elderly patients at particular risk for overprescribing? A report from the more & Romsdal prescription study. Scand J Prim Health Care. 1997;15(1):16–21.
Millar A, Hughes C, Ryan C. Evaluating the prevalence of potentially inappropriate prescribing in older adults in intermediate care facilities: a cross-sectional observational study. Int J Clin Pharmacol. 2017;39(3):527–35.
Siriwardena AN, Qureshi Z, Gibson S, Collier S, Latham M. GPs’ attitudes to benzodiazepine and “Z-drug” prescribing: a barrier to implementation of evidence and guidance on hypnotics. Br J Gen Pract. 2006;56(533):964–7.
Siriwardena AN, Qureshi MZ, Dyas JV, Middleton H, Orner R. Magic bullets for insomnia? Patients’ use and experiences of newer (Z drugs) versus older (benzodiazepine) hypnotics for sleep problems in primary care. Br J Gen Pract. 2008 Jun;58(551):417–22.
Cook JM, Marshall R, Masci C, Coyne JC. Physicians’ perspectives on prescribing benzodiazepines for older adults: a qualitative study. J Gen Intern Med. 2007;22(3):303–7.
Peritogiannis V, Manthopoulou T, Mavreas V. Long-term benzodiazepine treatment in patients with psychotic disorders attending a mental health Service in Rural Greece. J Neurosci Rural Pract. 2016;7(Suppl 1):S26–30.
ABIM Foundation. Choosing Wisely. Available from: http://www.choosingwisely.org/. Accessed 21 Dec 2017.
Olano-Espinosa E, Minue-Lorenzo C. “Do not do” also as regards tobacco. Aten Primaria. 2016;48(7):493–9.
Ministerio de Sanidad Servicios Sociales e Igualdad. Compromiso por la Calidad de las Sociedades Científicas para consensuar recomendaciones «no hacer» basadas en la evidencia científica. 2013. Available from: http://www.msssi.gob.es/organizacion/sns/planCalidadSNS/cal_sscc.htm.
Correas Lauffer J, Braquehais Conesa D, Barbudo Del Cura E, Ochoa ME. Abuse, tolerance and dependence of zolpidem: three case reports. Actas Esp Psiquiatr. 2002;30(4):259–62.
Hollingworth SA, Siskind DJ. Anxiolytic, hypnotic and sedative medication use in Australia. Pharmacoepidemiol Drug Saf. 2010;19(3):280–8.
Carrasco-Garrido P, Hernández-Barrera V, Jiménez-Trujillo I, et al. Time Trend in Psychotropic Medication Use in Spain: A Nationwide Population-Based Study. Tchounwou PB, ed. Int J Env Res Public Health. 2016;13(12):1177.
Vicente Sanchez MP, Macias Saint-Gerons D, de la Fuente HC, Gonzalez Bermejo D, Montero Corominas D, Catala-Lopez F. Trends of use of anxiolytics and hypnotics in Spain from 2000 to 2011. Rev Esp Salud Publica. 2013;87(3):247–55.
Ciuna A, Andretta M, Corbari L, Levi D, Mirandola M, Sorio A, et al. Are we going to increase the use of antidepressants up to that of benzodiazepines? Eur J Clin Pharmacol. 2004;60(9):629–34.
Sonnenberg CM, Bierman EJM, Deeg DJH, Comijs HC, van Tilburg W, Beekman ATF. Ten-year trends in benzodiazepine use in the Dutch population. Soc Psychiatry Psychiatr Epidemiol. 2012;47(2):293–301.
Simo MJ. Use of prescription drugs in Spain and Europe. Aten Primaria. 2012;44(6):335–47.
Stephenson CP, Karanges E, McGregor IS. Trends in the utilisation of psychotropic medications in Australia from 2000 to 2011. Aust N Z J Psychiatry. 2013;47(1):74–87.
Ministerio de Sanidad servicios sociales e igualdad. Use of anxiolytic and hypnotic drugs in Spain during the period 2000-2012. Drug utilization report. Madrid; 2015. Available from: https://www.aemps.gob.es/medicamentosUsoHumano/observatorio/docs/ansioliticos_hipnoticos-2000-2012.pdf. Accessed 21 Dec 2017.
Lillo Fernández de Cuevas JM, Rodríguez Blas M del C. Estadística de Gasto Sanitario Público 2015. Minist Sanidad, Serv Soc e Igual. 2015;1–18.
Agencia Española del Medicamento y Productos sanitarios (AEMPS). Zolpidem. 2015. Available from: https://www.aemps.gob.es/en/informa/notasInformativas/medicamentosUsoHumano/seguridad/2014/NI-MUH_FV_05-2014-zolpidem.html.
Gale CK, Millichamp J. Generalised anxiety disorder. BMJ Clin Evid 2011;2011.
Ministerio de Sanidad servicios sociales e igualdad. Methods to mesure drugs use. Observatorio de uso de medicamentos. 2014. Available from: https://www.aemps.gob.es/medicamentosUsoHumano/observatorio/metodologia.html.
National Institute of statistics. Statistics of the continuous register as of January 1, 2015. 2015 [cited 2017 Nov 21]. Available from: http://www.ine.es/jaxi/Datos.htm?path=/t20/e245/p04/a2015/l0/&file=00041002.px.
WHO Collaborating Centre for Drug Statistics Methodology. Guidelines for ATC classification and DDD assignment. Oslo; 2013. https://www.whocc.no/filearchive/publications/1_2013guidelines.pdf.
Nicieza-Garcia ML, Alonso-Lorenzo JC, Suarez-Gil P, Rilla-Villar N. Effect of the economic crisis on consumption of psychotropic drugs in Asturias (Spain). Gac Sanit. 2016;30(6):464–7.
Cunningham CM, Hanley GE, Morgan S. Patterns in the use of benzodiazepines in British Columbia: examining the impact of increasing research and guideline cautions against long-term use. Health Policy (New York). 2010;97(2–3):122–9. [cited 2017 Nov 22] Available from: http://www.sciencedirect.com/science/article/pii/S0168851010001028?via%3Dihub
Ferrer López I, Segura del Real F, Benavente Cantalejo SR, Fernández Urrusuno R SPA. Impact of a pharmaceutical-medical intervention to minimize patients at risk due to high doses of zolpidem. Farm comunitarios. 2015;6(Sup. 1).
Alvarez-Galvez J, Salinas-Perez JA, Rodero-Cosano ML, Salvador-Carulla L. Methodological barriers to studying the association between the economic crisis and suicide in Spain. BMC Public Health. 2017;17(1):694.
Navarro Vicente M, Sempere Verdú E, Pereira García A, MM SSA MJ. Crisis and mental health in Paterna (Valencia). FML. 2013;17(16):4p.
Ministerio de Sanidad servicios sociales e igualdad. Real Decreto-ley 16/2012, de 20 de abril, de medidas urgentes para garantizar la sostenibilidad del Sistema Nacional de Salud y mejorar la calidad y seguridad de sus prestaciones. Available from: https://www.boe.es/diario_boe/txt.php?id=BOE-A-2012-5403. Accessed 21 Dec 2017.
Sanchez DP, Guillen JJ, Torres AM, Arense JJ, Lopez A, Sanchez FI. Prescription drug consumption recovery following the co-payment change: evidence from a regional health service. Aten Primaria. 2015;47(7):411–8.
The study was funded by the CIBERESP (Ciber of Epidemiology and Public Health), Spain. The funding body had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Availability of data and materials
These data are available from the Seville Primary Care Pharmacy Clinical Management Unit, but restrictions limit access to these data, which were used under a data use agreement for the current study and therefore are not publicly available. Data are however available from the authors upon reasonable request and with permission of the Primary Care Pharmacy Clinical Management Unit.
Ethics approval and consent to participate
Ethical approval was obtained from the research committee of the Andalusian School of Public Health.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.