key: cord-0733638-y6ui1wwd
authors: Antonazzo, Ippazio Cosimo; Fornari, Carla; Maumus-Robert, Sandy; Cei, Eleonora; Paoletti, Olga; Ferrara, Pietro; Conti, Sara; Cortesi, Paolo Angelo; Mantovani, Lorenzo Giovanni; Gini, Rosa; Mazzaglia, Giampiero
title: Antidepressants Drug Use during COVID-19 Waves in the Tuscan General Population: An Interrupted Time-Series Analysis
date: 2022-01-28
journal: J Pers Med
DOI: 10.3390/jpm12020178
sha: 3de5ea52d45730895dc4c80b0911ff2b2674c930
doc_id: 733638
cord_uid: y6ui1wwd

In Italy, during the COVID-19 waves two lockdowns were implemented to prevent virus diffusion in the general population. Data on antidepressant (AD) use in these periods are still scarce. This study aimed at exploring the impact of COVID-19 lockdowns on prevalence and incidence of antidepressant drug use in the general population. A population-based study using the healthcare administrative database of Tuscany was performed. We selected a dynamic cohort of subjects with at least one ADs dispensing from 1 January 2018 to 27 December 2020. The weekly prevalence and incidence of drug use were estimated across different segments: pre-lockdown (1 January 2018–8 March 2020), first lockdown (9 March 2020–15 June 2020), post-first lockdown (16 June 2020–15 November 2020) and second lockdown (16 November 2020–27 December 2020). An interrupted time-series analysis was used to assess the effect of lockdowns on the observed outcomes. Compared to the pre-lockdown we observed an abrupt reduction of ADs incidence (Incidence-Ratio: 0.82; 95% Confidence-Intervals: 0.74–0.91) and a slight weekly decrease of prevalence (Prevalence-Ratio: 0.997; 0.996–0.999). During the post-first lockdown AD use increased, with higher incidence- and similar prevalence values compared with those expected in the absence of the outbreak. This pandemic has impacted AD drug use in the general population with potential rebound effects during the period between waves. This calls for future studies aimed at exploring the mid–long term effects of this phenomenon.

The COVID-19 pandemic caused by the SARS-CoV-2 virus has resulted in more than 228 million verified cases of virus infections and over 4.6 million COVID-19 related deaths globally [1] . Italy was one of the first countries in Europe to be hit by the new virus [2] [3] [4] . The emergency situation, characterized by an exponential growth of cases with high fatality rates, forced the Italian government to impose the first national lockdown from 9 March 2020 to 15 June 2020 to stop the virus spread in the general population [5] . In the following summer months, restrictions have been eased and people were allowed to work and travel. This resulted in a second large COVID-19 wave with a higher number of cases and consequent implementation of new lockdown restrictions from November 2020 to the end of December 2020 [5] . The measures adopted during the two lockdowns encompassed social distancing, smart-working, quarantine, people movement restriction, closure of non-essential activities, and limitation of the physical interactions.

The impact of stressful events on mental health has been already documented in the literature [6, 7] . In particular, during the first COVID-19 wave an increased number of psychiatric diagnoses such as post-traumatic stress disorder (PTSD), depression, and anxiety were reported [7] [8] [9] [10] [11] [12] [13] [14] . Previous studies already investigated the impact of lockdown on antidepressant drug (AD) use. In particular, some studies reported an increased use of these medications during the first wave [15, 16] that in some case persist also during the subsequent period of 2020 [17] [18] [19] [20] . On the contrary, other studies found an initial reduction of AD use during the peak of the first wave followed by an increased use in the subsequent periods [21] . Finally, some authors reported a non-significant change in psychotropic medication dispensing, including ADs, during the outbreak [22] . To the best of our knowledge, evidence on the impact of lockdowns on AD use in Italy is still scant and deserves careful evaluation.

This study, therefore, aimed at assessing whether the lockdown measures implemented in Italy during COVID-19 waves influenced prevalence and incidence of AD use in the general population.

The method used in this study have been already described in other publications aimed at assessing the impact of lockdown on other drug classes in the Tuscan general population [23, 24] . In this observational retrospective study, the healthcare administrative database (HAD) of Tuscany was used, covering approximately 3.6 million inhabitants. The HAD contains data on healthcare services accesses in the regional area, reimbursed by the Italian Healthcare System (NHS) to all regional citizens. The HADs consulted for this study included: (i) the archive of demographic and administrative data of individuals living in the catchment area who receive national healthcare service assistance; and (ii) the pharmacy claim registry that provides information (i.e., dispensing date, number of packages, substance name, anatomical therapeutic chemical code [ATC] and defined daily dose [DDD]) on all community prescription reimbursed by the NHS [25] . Approval for the study was obtained from the "Agenzia Regionale di Sanità della Toscana" Internal Governance Board.

All individuals with at least one dispensing of ADs (ATC: N06A*) between 1 January 2018 and 27 December 2020 were selected. The date of the first AD dispensing in the period was considered as the index date in the study. Subjects registered in the HAD less than 1 year before the index date were excluded from the cohort. Each patient was followed from the index date until 27 December 2020 or death/emigration, whichever came first.

During follow-up, patients were considered exposed to ADs when they claimed an AD dispensing of the same drug class (selective serotonin reuptake inhibitors-SSRIs, other antidepressants and tricyclic antidepressants). They were considered exposed until the end of duration of their last consecutive AD dispensing (treatment episode). We considered a dispensing as consecutive if it was claimed within the duration of the prior dispensing plus a 30-day grace period, which was added to account for irregular refill patterns and minor non-compliance. The duration of each dispensing was calculated by dividing the total quantity of active substance dispensed by the relevant defined daily dose (DDD) [26, 27] . The number of treatment episodes was used to estimate the number of patients with incidental or prevalent AD use in each week of the study period. Patients switching therapy or using concomitantly more than one AD in the week, were counted as just one user.

The outcomes of interest were weekly prevalence and incidence of AD use during the observation period. In particular, the weekly prevalence of AD use was calculated by dividing the weekly number of patients under AD treatment by the number of inhabitants living in the region at 1 January of each corresponding calendar year as reference population [28] . Similarly, the weekly incidence of AD use was estimated by including only patients who started an incident treatment episode at the numerator and the same reference population. For each treatment episode, we classified the patient as incident if he/she did not have an AD dispensed in the year prior the starting date of the treatment episode. During the study period, patients might have had more than one incident treatment episode.

First, we estimated the weekly prevalence and incidence of ADs use in each segment of the study period: pre-lockdown period from 1 January 2018 to 8 March 2020 (114 weeks), the first lockdown from 9 March 2020 to 15 June 2020 (14 weeks), the post-First lockdown from 16 June 2020 to 15 November 2020 (22 weeks) and the initial period of second lockdown from 16 November 2020 to 27 December 2020 (6 weeks). Second, the impact of lockdown restrictions on the study outcomes was assessed by performing an interrupted time series (ITS) analysis with a quasi-Poisson generalized additive model [29] [30] [31] . The model included the weekly count of the observed outcome as response variable (Y) and the reference population as offset variable to transform the count outcome in incidence or prevalence. The fitted model was of the following form:

The model included a non-linear function of the week (f(week i ), spline function) and a dummy holiday indicator (0 = no, 1 = yes) to account for time trend and seasonality. Model coefficients were: β 0 , which represents the baseline level of outcome during the pre-lockdown period; β 2 , β 4 and β 6 which estimate the level change during first, post and second lockdown, respectively and β 3 , β 5 and β 7 which estimate the trends/slopes of the time series during the aforementioned periods [31, 32] . Negative β i coefficients indicate reduction over the corresponding period. In the model, a level change means an abrupt effect of the intervention whereas a change in trend/slope represents a gradual change in the estimated outcome [31] . In this specific case, a level change indicated an immediate and sustained effect on drug use, whereas a trend/slope change imply that the use of the antidepressants changed gradually in the studied segment. Next, we also investigated a possible delayed effect of lockdown restrictions implementation or cancellation, using the delayed or "lagged" level and slope indicators for all segments.

Finally, we plotted graphs of the outcomes over time to illustrate trends. Statistical significance of the parameters and the goodness of fit of the model were used to choose the best model. Significance was defined as a p-value less than 0.05. The models were implemented separately for each outcome.

In this study both the data processing and data analysis were performed using the R studio software (version 4.0.2 Rstudio, PBC: Boston, MA, USA).

The research was conceived and performed by public institutions. All authors are employed or collaborate with public institutions.

As reported in Table 1 and Figure 1A , a variation of the ADs prevalence was observed during the COVID-19 pandemic. In particular, the ITS analysis suggested a slight increase of weekly prevalence of ADs use in the first week of the first lockdown (Prevalence Ratio: 1.01; 95% CI: 1.00-1.02) followed by a weekly reduction (0.997; 0.996-0.999) for the entire first lockdown segment. During the post-First lockdown segment we observed a level change indicating a rapid increase of AD prevalence (1.02; 1.01-1.03) in the first week and slope change (1.002; 1.001-1.003) indicating a weekly increase of the prevalence for the rest of the period (Table 1 and Figure 1A ). During the initial phase of second lockdown, no significant variation in AD prevalence was observed (Table 1 and Figure 1A ). 

In the study period, a variation in the trend of AD incidence was observed in the Tuscan general population across the different segments ( Figure 1B ). As shown in Table 2 , individuals were mainly female and aged 40 years and older. Among the study ADs drug classes, the most used was the selective serotonin reuptake inhibitors (SSRIs), whereas the less used were tricyclic ADs ( Figure A1 ). As reported in Figure 1B , a level change, indicating an abrupt reduction of the AD incidence was observed during the first week of the first lockdown segment (Incidence Ratio: 0.82; 95% CI: 0.74-0.91). At post-First lockdown a level change, which indicates a rapid increase of AD incidence (1.30; 1. 16-1.45) in the first week, was shown. During the initial phase of second lockdown, no significant variation was observed ( Figure 1B and Table 1 ). 

The findings of this observational descriptive study showed that an abrupt reduction of AD incidence and a gradual reduction of AD prevalence occurred in the general population during the first lockdown. This initial reduction was followed by new increase of both prevalence and incidence during the subsequent periods (post-first lockdown and second lockdown). In particular, the AD prevalence reverted to values similar to those registered in the pre-lockdown period following the secular trend. On the contrary, after the firstlockdown the AD incidence reached higher values compared with expected values in case of no lockdown measures implementation.

Our data on AD use during the first lockdown appeared to be consistent with those reported in previous studies from Italy and other European countries that analysed mental health services access during the first lockdown [33] [34] [35] [36] . The first COVID-19 outbreak has disrupted these services in most countries, especially during the first weeks of lockdown implementation. This might have caused an abrupt reduction of new diagnosis of psychiatric diseases with consequent reduction of new treatment cycle initiation. In parallel, the slight increase of the AD prevalence observed around the first week of the first lockdown and its subsequent reduction during the following weeks might result by combination of two phenomena. In fact, it is possible to speculate that during the early phase of virus spread (end of February) stress, worries and fear may have been the precursor of clinical exacerbation of existing anxiety/depression symptoms with consequent re-initiation of AD treatment in patients who had previously used ADs. On the contrary, during the subsequent weeks of lockdown patients might have been more reluctant to seek treatment, following the pattern of behaviour already described in other medical fields [34] . It is also possible that during the COVID-19 outbreak physicians were less prone to prescribe ADs to these patients by preferring other alternatives such as benzodiazepines. This hypothesis is supported by findings from Balestri and colleagues that documented an increased use of benzodiazepines in the psychiatric setting during the first lockdown [37] .

During the post-first lockdown period, an increased use of ADs was registered. In particular, our results showed that the AD prevalence returned to pre-lockdown levels, whereas the number of new AD treatments were higher during this period compared with expected values in case of no lockdown measures implementation. In this regard, it should be noted that during the last phase of the first-lockdown and subsequent periods, several digital solutions were implemented to guarantee medical assistance to psychiatry patients [38] . This might have represented a crucial event for medical/therapeutic assistance to patients who experienced anxiety/depression symptoms onset during the last part of the first lockdown and those who had symptoms exacerbation with consequent mental health problems [39, 40] . Additionally, data on incidence of drug use might also suggest a potential rebound/mid-long term effect of lockdown. In fact, it is possible to speculate that the emergency situation might have contribute to development of mental problems in individuals without pre-existing mental health problems who were resilient during the early phase of lockdown [37, 41] . In this regard, stressors such as fear of contagion, long quarantine, isolation and movement restriction, lack of information, financial loss, inadequate supplies and stigma, and loss of relatives due to COVID-19 might have contributed to the deterioration in psychological outcomes with consequent onset of anxiety, depression and related symptoms [42] [43] [44] .

The mid/long term consequences of the under-diagnosis and under-treatment of anxiety and depression during the first lockdown are still under-investigation. Previous studies suggested that an increase in psychiatric disorders and suicidal ideations and attempts might appear several months after stressful events even if a temporary reduction might be observed during the emergency situation [45] [46] [47] [48] [49] . This highlights the need of continuous monitoring of patients with psychiatry problems under exceptional circumstances like future COVID-19 waves. From clinical point of view, our findings suggest the necessity to encourage new tools such as telemedicine in order to guarantee medical assistance to psychiatric patients. This could be crucial to reduce the risk of adverse events due to non-optimal drug treatment. Additionally, telemedicine could increase also the effect of non-pharmacological management of these patients contributing to the aforementioned events improvements.

Future research should be performed to investigate whether new users of ADs during the pandemic outbreak had different demographic and clinical characteristics compared with those who started an AD treatment in the pre-lockdown period. In addition, new studies are warranted to investigate the potential consequence associated with the observed variations in AD use in the study population.

An important strength of this study is that data related to an unselected population from one of the largest Italian regions, which accounts for over over 3.5 million of inhabitants. This study has some limitations. First, in Italy both ADs and anxiolytic drugs (i.e., benzodiazepines) are indicated for the treatment of the mental health symptoms after traumatic events. However, in Italy only ADs are fully reimbursed by the NHS and thus captured in the HAD, while the cost of anxiolytics is totally charged to patients. As a consequence, data about benzodiazepines use is generally not available in the Italian HAD. Moreover, the study did not include information about drug indications because in Italy this information is not available in the HAD. Additionally, for the second lockdown we included only six data points, which is slightly lower than the minimum number of points (eight) required to have robust estimates, therefore the results on the second lockdown should be considered with caution. Finally, as per study design we were not able to ensure a causal correlation between COVID-19 waves and studied outcomes, although the observed changes in the ADs use occurred around the studied events.

In conclusion, we found a change in AD use during pandemic lockdowns in the general population. There was a reduction of both prevalence and incidence of AD use during the first lockdown that reversed during the post-first lockdown reaching higher value compared with the expected results observed in the counterfactual scenario. The potential undertreatment occurred during the pandemic waves might be associated with increased risk of negative events (i.e., suicidal attempts) occurrence in the general population. These findings might be used to implement strategies to ensure access to adequate mental healthcare services for vulnerable population during extreme circumstances such as lockdown periods. Institutional Review Board Statement: This is an observational, non-interventional database study, reutilising administrative data for the purpose of addressing a research question. The Agenzia regionale di Sanità della Toscana has general legal and ethical frameworks that allow conducting research by making secondary use of administrative data, upon approval of an internal governance board. Permission to contribute to this study was granted by the governance boards of the Agenzia regionale di Sanità della Toscana.

Informed Consent Statement: Not applicable. The Agenzia regionale di Sanità della Toscana has general legal and ethical frameworks that allow conducting research by making secondary use of administrative data, upon approval of an internal governance board.

Conflicts of Interest: I.C.A., C.F., S.C., P.F., P.C. and G.M. have no disclosure to declare. L.G.M. reported receiving grants from Bayer, Daiiki-Sankyo, and Boehringer Ingelheim outside the submitted work and speaker fees from Pfizer and Bayer. O.P. and R.G. are employed by ARS, a public health agency that conducts or participates in pharmacoepidemiology studies compliant with the ENCePP Code of Conduct. The budget of ARS is partially sustained by such studies. S.M.-R. is part of DRUGS-SAFER centre of pharmacoemidemiology, which is granted by the Agence Nationale de Sécurité du Médicament et des Produits de Santé (French Medicine Agency). The centre is conducting drug utilization studies on the impact of COVID-19, in which SMR is not involved.

WH Coronavirus (COVID-19) Dashboard

The COVID-19 Pandemic: A Comprehensive Review of Taxonomy, Genetics, Epidemiology, Diagnosis, Treatment, and Control

Clinical characteristics and differential clinical diagnosis of novel coronavirus disease 2019 (COVID-19)

Magnitude and time-course of excess mortality during COVID-19 outbreak: Population-based empirical evidence from highly impacted provinces in northern Italy

Ulteriori Disposizioni Attuative del Decreto-Legge 23 Febbraio 2020, n. 6, Recante Misure Urgenti in Materia di Contenimento e Gestione Dell'emergenza Epidemiologica da COVID-19

The psychological impact of quarantine and how to reduce it: Rapid review of the evidence

Global prevalence and burden of depressive and anxiety disorders in 204 countries and territories in 2020 due to the COVID-19 pandemic

Impact of the coronavirus 2019 (COVID-References 1. World Health Organization (WHO). WH Coronavirus (COVID-19) Dashboard

The COVID-19 Pandemic: A Comprehensive Review of Taxonomy, Genetics, Epidemiology, Diagnosis, Treatment, and Control

Clinical characteristics and differential clinical diagnosis of novel coronavirus disease 2019 (COVID-19)

Magnitude and time-course of excess mortality during COVID-19 outbreak: Population-based empirical evidence from highly impacted provinces in northern Italy

Ulteriori Disposizioni Attuative del Decreto-Legge 23 Febbraio 2020, n. 6, Recante Misure Urgenti in Materia di Contenimento e Gestione Dell'emergenza Epidemiologica da COVID-19, Applicabili Sull'intero Territorio Nazionale

The psychological impact of quarantine and how to reduce it: Rapid review of the evidence

Global prevalence and burden of depressive and anxiety disorders in 204 countries and territories in 2020 due to the COVID-19 pandemic

Impact of the coronavirus 2019 (COVID-19) pandemic on anxiety diagnosis in general practices in Germany

Increased generalized anxiety, depression and distress during the COVID-19 pandemic: A cross-sectional study in Germany

Depression, Anxiety, and Stress in the time of COVID-19 Pandemic in Kurdistan Region

A Nationwide Survey of Psychological Distress among Italian People during the COVID-19 Pandemic: Immediate Psychological Responses and Associated Factors

Assessing the anxiety level of Iranian general population during COVID-19 outbreak

Psychological effects of the COVID-19 outbreak and lockdown among students and workers of a Spanish university

Depression, Anxiety and stress during COVID-19: Associations with changes in physical activity, sleep, tobacco and alcohol use in australian adults

Comparing actual and forecasted numbers of unique patients dispensed select medications for opioid use disorder, opioid overdose reversal, and mental health, during the COVID-19 pandemic, United States

The effects of COVID-19 on the dispensing rates of antidepressants and benzodiazepines in Canada

Antidepressants, primary care, and adult mental health services in England during COVID-19

Trends in the Use of Benzodiazepines, Z-Hypnotics, and Serotonergic Drugs among US Women and Men Before and During the COVID-19 Pandemic

Assessment of Psychotropic Drug Prescribing Among Nursing Home Residents in Ontario, Canada, During the COVID-19 Pandemic

Changes in antidepressant use in Australia: A nationwide analysis prior to and during the COVID-19 pandemic

Effects of the COVID-19 pandemic on primary care-recorded mental illness and self-harm episodes in the UK: A population-based cohort study

Impact of the COVID-19 Pandemic on the General Mental Health in Sweden: No Observed Changes in the Dispensed Amount of Common Psychotropic Medications in the Region of Scania. Front

COVID-19 Outbreak Impact on Anticoagulants Utilization: An Interrupted Time-Series Analysis Using Health Care Administrative Databases

Impact of COVID-19 Lockdown, during the Two Waves, on Drug Use and Emergency Department Access in People with Epilepsy: An Interrupted Time-Series Analysis

The Role of European Healthcare Databases for Post-Marketing Drug Effectiveness, Safety and Value Evaluation: Where Does Italy Stand? Drug Saf

Comparison of dose standard units for drug utilisation studies

Methods for constructing treatment episodes and impact on exposure-outcome associations

Italian Office of National Statistics-Istituto Nazionale di Statistica (ISTAT)

Use of Interrupted Time Series Analysis in Evaluating Health Care Quality Improvements

Regression based quasi-experimental approach when randomisation is not an option: Interrupted time series analysis

Segmented regression analysis of interrupted time series studies in medication use research

Time Series Analysis: The Cross-Correlation Function

Psychiatric Admissions, Referrals, and Suicidal Behavior Before and During the COVID-19 Pandemic in Denmark: A Time-Trend Study

The Early Impact of the COVID-19 Pandemic on Acute Care Mental Health Services

Psychiatric hospitalization rates in Italy before and during COVID-19: Did they change? An analysis of register data

The effect of the March 2020 COVID-19 lockdown on national psychiatric contacts in Denmark: An interrupted time series analysis

Emergency Psychiatric Consultations During and After the COVID-19 Lockdown in Italy

Adoption of Digital Technologies in Health Care During the COVID-19 Pandemic: Systematic Review of Early Scientific Literature

The emotional impact of Coronavirus 2019-nCoV (new Coronavirus disease)

Patients with mental health disorders in the COVID-19 epidemic

A psychiatrist's perspective from a COVID-19 epicentre: A personal account

The anatomy of panic buying related to the current COVID-19 pandemic

Coping behaviors associated with decreased anxiety and depressive symptoms during the COVID-19 pandemic and lockdown

The COVID-19 Global Pandemic: Implications for People with Schizophrenia and Related Disorders

Effect of 11 September 2001 terrorist attacks in the USA on suicide in areas surrounding the crash sites

Effect of the First World War on suicide rates in Ireland: An investigation of the 1864-1921 suicide trends

Psychosocial characteristics as potential predictors of suicide in adults: An overview of the evidence with new results from prospective cohort studies

Patients at high risk of suicide before and during a COVID-19 lockdown: Ecological momentary assessment study

Consequences of the COVID-19 pandemic on admissions to general hospital psychiatric wards in Italy: Reduced psychiatric hospitalizations and increased suicidality

J. Pers. Med. 2022, 12, x FOR PEER REVIEW 9 of 11 Figure A1 . Incidence of different antidepressant drugs (ADs) classes use in the general population in the study period. Legend: Orange dots: estimated incidence of Selective Serotonin Reuptake inhibitors (SSRIs); blue dots: estimated incidence of Others antidepressants and red dots: estimated incidence of tricyclic antidepressants. Figure A1 . Incidence of different antidepressant drugs (ADs) classes use in the general population in the study period. Legend: Orange dots: estimated incidence of Selective Serotonin Reuptake inhibitors (SSRIs); blue dots: estimated incidence of Others antidepressants and red dots: estimated incidence of tricyclic antidepressants.