The impact of the COVID-19 pandemic on diagnoses of common mental health disorders in adults in Catalonia, Spain

Objectives: To investigate how incidence trends of anxiety and depressive disorders have been affected by the COVID-19 pandemic. Design: Population-based cohort study. Setting: Observational cohort study from 2018 to 2021 using the Information System for Research in Primary Care (SIDIAP) database in Catalonia, Spain. Participants: 4,255,847 individuals aged 18 or older in SIDIAP on 1 March, 2018 with no prior history of anxiety and depressive disorders. Primary and secondary outcomes measures: Incidence of anxiety and depressive disorders prior to COVID-19 (March, 2018 to February, 2020), during the COVID-19 lockdown (March to June, 2020) and post-lockdown periods (from July, 2020 to March, 2021) were calculated. Forecasted rates over COVID-19 periods were estimated using negative binomial regression models based on previous data. The percentage reduction was estimated by comparing forecasted versus observed events, overall and by age, sex and socioeconomic status. Results: The incidence rates per 100,000 person-months of anxiety and depressive disorders were 171.0 (95%CI: 170.2-171.8) and 46.6 (46.2-47.0), respectively, during the pre-lockdown period. We observed an increase of 39.7% (95%PI: 26.5 to 53.3) in incident anxiety diagnoses compared to the expected in March, 2020, followed by a reduction of 16.9% (8.6 to 24.5) during the post-lockdown periods. A reduction of incident depressive disorders occurred during the lockdown and post-lockdown periods (46.6% [38.9 to 53.1] and 23.2% [12.0 to 32.7], respectively). Reductions were higher among adults aged 18 to 34 and individuals living in most deprived areas. Conclusions: The COVID-19 pandemic in Catalonia was associated with an initial increase in anxiety disorders diagnosed in primary care, but a reduction in cases as the pandemic continued. Diagnoses of depressive disorders were lower than expected throughout the pandemic.


Summary box:
What is already known on this topic -While previous self-reported studies have provided evidence of increased mental health burden during the initial phase of the COVID-19 pandemic, a number of studies observed that fewer diagnoses were made in primary care settings than would have been expected during the initial stages of the pandemic.
-Population data that examine the impact of COVID-19 on temporal trends of incident cases of common mental health disorders are lacking in Catalonia, Spain.

What this study adds
-This study has quantified the impact of the COVID-19 pandemic on trends of incidence of anxiety and depressive disorders among adults living in Catalonia.
-Reductions in incident cases of anxiety and depressive disorders were higher for young adults and people living in most deprived areas.
-Incident diagnoses of anxiety and depressive disorders have not been fully recovered to what would have been expected.

INTRODUCTION
The outbreak of the COVID-19 pandemic and the associated control measures have impacted many aspects of people's lives. In Spain, a national lockdown was implemented on 14 March. 1 Strict control measures were adopted in order to contain the spread of the virus, limiting mobility with few exceptions such as grocery shopping, health emergencies or essential work.
The restrictions were gradually lifted following different de-escalation phases, which started on 28 May, 2020 and ended on 21 June, 2020. 2 Due to a resurgence of COVID-19 infections, new nationwide measures were implemented on 25 October, 2020, including mobility restrictions and curfew hours, which were extended up to 9 May, 2021. 3 Studies based on self-reported surveys have provided evidence of elevated rates of anxiety, depression and stress in the initial stages of the pandemic. [4][5][6][7] In Catalonia (a northern region of Spain), self-reported surveys performed in April-May 2020 found increased levels of anxiety and depression during COVID-19 lockdown, especially among women and young adults. 8,9 Conversely, health care contacts related to mental conditions were substantially reduced in primary care, emergency departments and hospital settings after the lockdown announcement in March, 2020. 10 Serious concerns have been raised about the long-term and far-reaching mental health impact of the pandemic. 11 A discrepancy between increasing levels of mental health disorders and reductions in primary care diagnoses during the COVID-19 pandemic may well result in a substantial burden of ill-health, as long delays in diagnosis are associated with negative health outcomes among adults with mood and/or anxiety disorders. 12,13 To date, most evidence regarding the mental health impact of COVID-19 is from self-report surveys performed during the initial stages of the pandemic. Yet, we do not know whether survey results are mirrored in the rates of recorded incident diagnoses of common mental disorders in primary care, and how they have evolved during the pandemic in Catalonia.
. CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted August 9, 2021. ;https://doi.org/10.1101https://doi.org/10. /2021 The use of longitudinal studies to monitor rates of mental disorders and to identify gaps in mental health care have been described as an urgent research priority in the response to  In this study, we aim to investigate how incidence trends of anxiety and depressive disorders have been affected by the various stages of the COVID-19 pandemic through the analysis of a large primary care longitudinal dataset representative of the population living in Catalonia, Spain.

Study design
We conducted a cohort study based on primary care records from 1 March 2018 to 31 March 2021 in Catalonia, Spain. Time trends were assessed in three different periods: 1) the prelockdown-19 period from 1 March 2018 to 29 February 2020; 2) the lockdown period from 1 March 2020 to 30 June 2020, and 3) the post-lockdown period, from 1 July 2020 to 31 March 2021. The lockdown period coincides with the home confinement and different phases of deescalation in Spain. The post-lockdown period was divided into three trimesters, which were broadly aligned with the different stages of the pandemic in Catalonia: 1) easing of restrictions from 1 July, 2020 to 30 September, 2020; 2) implementation of new control measures from 1 October, 2020 to 31 December, 2020; 3) extension of control measures and start of the vaccination campaign from 1 January, 2021 to 31 March, 2021.

Setting
The Catalan Health Care System dispenses services for 7.5 inhabitants, providing universal health coverage to residents of Catalonia through a tax-based system. 15 Primary care services in Spain are the first point of contact with healthcare services for the population, 16 and they carry the major burden of detection, management and treatment of mild mental health problems. 17

Study participants and data source
. CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted August 9, 2021. Individuals with less than one year of prior history available were excluded so that study participants had sufficient prior observation time to truly detect incident cases of anxiety and depressive disorders. Individuals with episodes of anxiety or depressive disorders prior to index date (i.e., date of start of the cohort) were excluded. Individuals were observed until the event of interest (anxiety or depressive disorders), until they were transferred or died, or until the study period ended.

Variables
Individuals' age and sex were extracted. Information on socioeconomic status (SES) was available through the Mortalidad en áreas pequeñas Españolas y Desigualdades Socioeconómicas y Ambientales (MEDEA) deprivation index, linked to each residential census area of the population. 20 The deprivation index was only available for urban areas, defined as municipalities with more than 10,000 inhabitants and a population density greater than 150 habitants/km 2 ; remaining areas were considered rural areas. The deprivation index is categorized in quintiles, in which the first and fifth quintiles are the least and most deprived, respectively.
The outcomes of the study were incidence of anxiety or depressive disorders. Conditions were identified on the basis of the International Statistical Classification of Diseases and Related Health Problems, Tenth Revision, Clinical Modification (ICD-10-CM) codes of interest. The ICD-10-CM codes included for depressive disorders were F32 (depressive episode) and F33 . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted August 9, 2021. ;https://doi.org/10.1101https://doi.org/10. /2021 (recurrent depressive disorder). All descendants were included except the ones referring to episodes in partial or complete remission (F32. 4, F32.5 and F33.4). For anxiety disorder the F41 (other anxiety disorder) code and all its descendants were included.

Statistical methods
We first summarised the socio-demographics characteristics of individuals included in the study, with counts and percentages for categorical variables and median and interquartile ranges (IQR) for continuous variables. We structured data in a monthly time-series format with incident and person-months at risk aggregated with stratification by sex (males and females), age group (18 to 34, 35 to 64, >65 years) and SES (quintiles of MEDEA deprivation index).
Incidence rates with 95% CI were calculated for each month and study period by dividing the number of first recorded cases of anxiety and depression by 100,000 person-months at risk, overall and stratified by sex, age group and SES. Incidence Rate Ratios (IRR) with 95% CI were calculated to compare the differences in incidence of each strata of the population during the lockdown period and post-lockdown period (divided into trimesters) compared to the prelockdown period.
We used negative binomial regression models to estimate expected monthly incident cases from March, 2020 to March, 2021, using data collected in the pre-lockdown period. To account for possible seasonality and linear trends, we fitted calendar month as a categorical variable and time as a continuous variable. The number of months since the start of the study was considered as the unit of measurement for time.
The estimated number of underdiagnoses were calculated by subtracting the number of expected from the observed diagnoses. Reduction in diagnoses were calculated dividing the number of underdiagnoses by the expected diagnoses, and were estimated monthly and by study period (lockdown and post-lockdown, overall and divided by trimesters). Reduction in diagnoses are expressed as percentages with 95% prediction intervals (PI), 21 as previously reported in a similar study. 22 . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted August 9, 2021. ; https://doi.org/10. 1101 For each month during the COVID-19 study period, observed and expected incident counts were converted to rates using the observed person-month denominator. We plotted monthly expected rates and corresponding 95%PI against the observed rates.
To validate our modelling approach, we developed a negative binomial model based on two years of prior history (from March, 2017 to February, 2019). We forecasted our expected values from March, 2019 to February, 2020, and checked if values predicted from the model fell within the calculated 95%PI. We validated our approach for overall diagnoses and stratifying by sex, age group and SES.
All analyses were performed in R version 4.0.4.

Patient and public involvement
This research was done without patient involvement. Patients were not invited to comment on the study design and were not consulted to develop patient relevant outcomes or interpret the results. Patients were not invited to contribute to the writing or editing of this document for readability or accuracy.

RESULTS
Our dataset included 4,255,847 individuals, of whom 2,105,171 (49.5%) were female. The median age of individuals on the index date was 47 years (IQR 35 to 62). Detailed demographic data of the study population is included in the Supplementary Table S1.
During the lockdown period, the incidence rate per 100,000 person-months of anxiety disorders increased up to 201.6 (199.4 to 203.9) whereas rates for depressive disorders sharply decreased down to 29.5 per 100,000 person-months (28.6 to 30.3). Incidence rates of anxiety disorders were substantially reduced from July to September, 2020 (142.0 per 100,000 person-months [139.9 to 144.2]), but subsequently increased and exceeded the figures . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted August 9, 2021. from January to March, 2021) ( Figure 1). Detailed data of incidence rates stratified by age, gender and SES during the lockdown and post-lockdown periods is included in the Supplementary Materials.
Incidence Rate Ratios (IRR) comparing the differences in incidence of each strata of the population during the lockdown and post-lockdown periods compared to the pre-lockdown period are shown in Table 1. IRR of anxiety disorders were significantly higher during the lockdown period for all age groups and deprivation quintiles (IRR range 1.06 to 1.37) except for females aged between 18 to 34 years (IRR 1.00 [0.97 to 1.03]). Significant reduction in diagnoses of anxiety disorders were found for all strata during the following two trimesters of 2020. Individuals aged 18 to 34 years or individuals living in rural areas obtained significantly higher IRR during the first trimester of 2021 compared to the pre-lockdown period (ie, IRR is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted August 9, 2021. ; https://doi.org/10.1101/2021.08.06.21261709 doi: medRxiv preprint Incident diagnoses of anxiety disorders showed a sudden peak at the onset of the lockdown measures in Catalonia, and accounted for an increase of 39.7% (95%PI 26.5 to 53.3) new diagnoses compared to the expected in March, 2020 ( Figure 1A). This peak was particularly pronounced among individuals aged 18 to 64 years and individuals living in most deprived areas ( Figure 2A and Figure 3A). Compared to the younger groups, males aged 65 or older showed a lower increase of incident cases of anxiety disorders from April to May, 2020, which was much more subtle among their female peer group ( Figure 2A).
The estimated number of underdiagnoses and the percentage of reduction are presented stratified by age, sex and SES in Table 2. An increase of diagnoses of anxiety disorders was found during the lockdown period for all stratas except from individuals aged 18 to 34 years, as new diagnoses of anxiety disorders among this subgroup were lower than expected from April, 2020 to February, 2021 ( Figure 2A; Table 2). The increase of diagnoses was followed by a decrease of 19.4% (95%PI 11.4 to 26.6) compared to the expected from July to August, 2020. This percentage was progressively reduced during the following two trimesters (17.3% [9.3 to 24.8] from October to December, 2020; 14.3 % [5.6 to 22.6] from January to March, 2021). The reduction of incident cases of post-lockdown period was 16.9% (8.6 to 24.5) ( Table   2).
Reductions in incident cases of anxiety disorders were greater in females during the lockdown period and were slightly higher in males during the post-lockdown period (lockdown: 2.2% [-9.3 to 12.0] in females vs -5.7% [-17.4 to 4.5] in males; post-lockdown: 16.8% [6.7 to 25.5] in females vs 17.5% [8.0 to 25.8] in males; data not shown). Individuals aged 18 to 34 years had greater reductions of incident diagnoses of anxiety disorders during lockdown and postlockdown periods. All age groups showed a decrease in the percentage of reduction comparing the first trimester of 2021 with the trimester following the lockdown end, except from males aged 65 or over ( Table 2). is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted August 9, 2021. ; https://doi.org/10.1101/2021.08.06.21261709 doi: medRxiv preprint New cases of depressive disorders were substantially reduced for all age groups compared to the expected at the onset of the lockdown measures, accounting for a reduction of 67.2% (95%PI 62.5 to 71.3) compared to the expected in April, 2020 ( Figure 2B). Overall monthly incidence rates remained under the 95% PI from March, 2020 to February, 2021 ( Figure 1B).
We estimated that the difference between expected and observed diagnoses accounted for a reduction of 46.6% (38.7 to 53.1) during the lockdown period compared to the expected. This reduction was substantially decreased to 25.3% (14.6 to 34.5) from July, 2020 to September, 2020, and was progressively diminished during the following two trimesters ( in males, data not shown). All age groups showed a marked decline in the percentage of reduction of incident cases of depressive disorders during the trimester following the lockdown end ( Table 2). Reductions were progressively diminished over the following trimester for all age groups with a few exceptions. Males aged 35 to 64 years and 65 years or above experienced a subsequent increase in the reduction of diagnosis of depressive disorders. This increase was limited to the last trimester of 2020 among males aged 35 to 64 years, whereas progressively increased among males aged 65 or above, reaching a 10.6 % absolute increase comparing the first quarter of 2021 with the trimester following the lockdown end. Reductions among females were progressively diminished over time, but slightly increased among females aged 18 to 34 years from January to March 2021. While reductions among individuals aged 18 to 34 and 35 to 64 years were similar from July, 2020 to March, 2021, individuals aged 18 to 34 years were the age group with higher reductions during the lockdown and postlockdown periods ( Table 2). is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted August 9, 2021. ; https://doi.org/10.1101/2021.08.06.21261709 doi: medRxiv preprint Most deprived urban areas accounted for higher reductions of incident diagnoses of anxiety disorders (during the lockdown period) and depressive disorders (during the lockdown and post-lockdown periods) ( Table 2). Reductions were lower in rural areas compared to urban areas for both periods, but were particularly high for anxiety disorders from July to September, 2020 (23.7 % [95%PI 14.8 to 31.4]

DISCUSSION
In this large population-based cohort study, we observed large variations in the incidence of anxiety and depression disorders diagnosed in primary care in Catalonia by the different periods of the COVID-19 pandemic. Incidence of primary care-coded anxiety disorders among individuals aged 18 to 64 years showed a sudden peak compared to expected rates in March, 2020, and was followed by reduction of diagnoses during the next months. Contrary to anxiety disorders, the incidence of depressive disorders sharply decreased among all sociodemographic groups reaching its minimum value in April, 2020. Based on observed data, incidence of anxiety and depressive disorders seemed to reach pre-lockdown levels in the first quarter of 2021. According to our model, which captures the underlying upward trend of diagnoses over time, the incidence of anxiety and depressive disorders did not fully reach what would have been expected. These results suggest that the effects of the measures . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted August 9, 2021. ; https://doi.org/10.1101/2021.08.06.21261709 doi: medRxiv preprint implemented to control the spread of the pandemic are still negatively affecting the diagnoses of these mental health conditions in primary care.
Although the negative impact on the incidence of anxiety and depressive disorders appeared to occur across all ages, sex and SES, this was mostly pronounced among young adults or individuals living in the most deprived areas. Greater reductions were also found among females except for incident diagnoses of anxiety disorders during the post-lockdown period.
These findings are particularly concerning, as these groups have been previously identified as experiencing worse mental health burden due to the pandemic. 4 Incidence rate ratios provided by our study suggest that individuals aged 18 to 34 years were the only age group with significant increases in cases of anxiety and depressive disorders during the first trimester of 2021 compared to the pre-lockdown period. A possible explanation could be the increased impact of the pandemic on daily life, which has been heavily disrupted among the youth. Yet, they experienced the highest percentage of reductions according to our model during and after the lockdown period.
Reductions in common mental health diagnoses might be explained by the disruption of the normal functioning of healthcare systems due to the COVID-19 pandemic. In Catalonia, nonessential health activities were interrupted during lockdown in order to prioritize COVID-19 services, and health authorities advised against going to health centres except in the event of serious illness or urgent situations. While it is certainly true that healthcare systems reached capacity constraints, continuous dissemination of mass media messages reporting the strain of healthcare systems might have influenced health-care seeking behaviour. People might have delayed seeking care in order to avoid overburdening health systems or out of fear of contracting the disease. 23 In such an extraordinary situation, low perceived legitimacy of mental issues might have also prevented those psychologically suffering to discern their emotional distress as an appropriate reason to seek care. In this line, we believe that a large number of people experiencing psychological distress might have developed coping strategies without seeking help in the public healthcare system. Prior evidence suggests that clinicians . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted August 9, 2021. ; https://doi.org/10.1101/2021.08.06.21261709 doi: medRxiv preprint perceive the emotional distress of their patients, but have difficulties in specifying standardized psychiatric diagnoses. 17,24 Consultations with a recorded diagnosis of a psychological problem take longer, 25 and might require re-assessment of suspected cases to increase accuracy. 26 It seems likely that the diagnoses detected in primary care during the pandemic might account for more severe cases of mental health issues. In addition, it is important to emphasize the increasing use of telemedicine in Catalonia, where the share of telemedicine-based visits has increased to 56.2%, from 15.4% in the pre-pandemic year. 28 Since the pandemic, telemedicine tools serve a wider and more representative cross-section of the population in Catalonia, including underrepresented groups such as lower income individuals and those living in rural areas. 27 While telemedicine has ensured the continuity of care of many health care processes during lockdown, it might face limitations generating diagnoses that require more assessments. 28

A recent study performed in Central Catalonia (a health administrative region of Catalonia)
found an average decline of 31.1% in new diagnoses compared to 2019. 28 This study found an association between the intensity of the pandemic and the increased monthly declines in new diagnoses across ICD-10 groups, highlighting the strain on capacity of the healthcare system to identify and address the healthcare needs of the population during the pandemic.
Another study performed in Catalonia has found a reduction in primary care visits associated with diagnoses related to chronic pathologies, while visits associated with diagnoses related to socio-economic and housing problems have increased. 29 Both studies found moderate reductions in diagnoses of mental health conditions, which were lower than the ones observed for other diagnoses groups. 28,29 Similar studies performed in the UK have used the same modelling approach to forecast incidence of anxiety and depressive disorders during lockdown and the first months of its aftermath. 22,30 Larger reductions were also found for adults aged 18 to 64 years and for patients registered at practices in more deprived areas. 30 In contrast to our findings, none of these studies have found an increase in new cases of anxiety disorders during lockdown. Comparison between studies performed in other countries render difficulties . CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted August 9, 2021. ; https://doi.org/10.1101/2021.08.06.21261709 doi: medRxiv preprint given the different strategies adopted to control de COVID-19 pandemic, which could have also played a role in the incidence of common mental disorders.
The reduced number of new diagnoses observed compared with the expected numbers obtained in this study are most likely to represent a large number of disorders that have gone undiagnosed and untreated. Previous studies suggest that diagnostic delays of anxiety and depressive disorders are associated with poorer outcomes. 12,13 Possible consequences of this unmet need may include increased demand for mental health services, increased use of emergency departments for mental issues and heightened risk of suicide. The economic impacts of the pandemic will likely increase the risk of mental health problems and exacerbate health-care disparities. 31 Monitoring patterns in primary care recording of common mental disorders, will provide crucial information to ensure health services can meet future demand.
Data on psychotropic drug prescriptions, primary care contacts and secondary care referrals can provide valuable information on the management of anxiety and depressive disorders during the COVID-19 pandemic. Further research should also investigate the evolution of conditions over time using linkages to data from hospital emergency departments and outpatient mental health centers.
The main strengths of this study are the sample size and the real-world nature of the data.
Our study provided data from a broadly representative setting, and included more than 4 million people registered in primary care in Catalonia. The data used were obtained from primary care electronic health records, which have been shown to be a useful tool for research in many areas, including 33 In order to strengthen our results, we validated our modelling approach forecasting expected diagnoses from one year prior to the pandemic, is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted August 9, 2021. ;https://doi.org/10.1101https://doi.org/10. /2021 We acknowledge several limitations to our study. First, our data is limited to recorded diagnoses in primary care, which will have led to omit people with whom anxiety and depressive disorders were discussed or people who were treated for these conditions, but not entered as a diagnosis onto their clinical record. It is therefore likely that this is an underestimate of episodes of anxiety and depressive disorders. Second, reductions in diagnoses of anxiety and depressive disorders may have been underestimated as they were based on pre-pandemic data, and did not take into consideration an increase of incident cases of psychological distress due to the pandemic.
In conclusion, the sudden peak in primary-care recorded incident cases of anxiety disorders  is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted August 9, 2021. ;https://doi.org/10.1101https://doi.org/10. /2021 Contributors BR, EB and TDS conceived and designed the study. SFB, MA, EB, and TDS mapped source data to the OMOP CDM. BR, EB and AP performed the statistical analysis. BR drafted the initial version of the manuscript. All authors (BR, AP, CR, SFB, MA, AB, CJA, LMP, EB, TDS) interpreted the results, critically reviewed the manuscript and approved the final version for submission.

Ethics approval
This project was approved by the Clinical Research Ethics Committee of the IDIAPJGol (project code: 21/052-PCV).

Transparency statement
The lead author affirms that the manuscript is an honest, accurate, and transparent account of the study being reported; that no important aspects of the study have been omitted; and that any discrepancies from the study as originally planned have been explained.

Funding
This research received no specific grant from any funding agency in the public, commercial or not-for-profit sectors.

Data availability statement
In accordance with current European and national law, the data used in this study is only available for the researchers participating in this study. Thus, we are not allowed to distribute or make publicly available the data to other parties. However, researchers from public institutions can request data from SIDIAP if they comply with certain requirements. Further information is available online (https://www.sidiap.org/index.php/menusolicitudesen/application-proccedure) or by contacting SIDIAP (sidiap@idiapjgol.org).
. CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted August 9, 2021. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted August 9, 2021. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted August 9, 2021. ;https://doi.org/10.1101https://doi.org/10. /2021  is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted August 9, 2021. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted August 9, 2021. ; https://doi.org/10.1101/2021.08.06.21261709 doi: medRxiv preprint Incidence Rate Ratios (IRR) were calculated by dividing the incidence rate obtained for each strata during the COVID-19 periods (lockdown and post-lockdown periods) by the incidence rate obtained for the same strata during the pre-lockdown period. The MEDEA deprivation index is calculated at the census tract level in urban areas of Catalonia and is categorised in quintiles of deprivation.

Figure 1. Expected and observed incidence rates of anxiety and depressive disorders in primary care in Catalonia (March, 2018-March, 2021).
Points indicate monthly incidence rates of anxiety and depressive disorders. Number of expected cases (95% PI) were estimated with negative binomial models, using data from 1 March, 2018 to 29 February, 2020. Shaded areas in blue represent 95 % PI. Vertical lines show 1 March, 2020.

Figure 2. Expected and observed incidence rates of anxiety and depressive disorders stratified by age group and gender in primary care in Catalonia (March, 2019-March, 2021).
Points indicate monthly incidence rates of anxiety and depressive disorders. Number of expected cases (95% PI) were estimated with negative binomial models, using data from 1 March, 2018 to 29 February, 2020. Shaded areas in blue represent 95 % PI. Vertical lines show 1 March, 2020.

Figure 3. Expected and observed incidence rates of anxiety and depressive disorders stratified by MEDEA deprivation index in primary care in Catalonia (March, 2019-March, 2021).
Points indicate monthly incidence rates of anxiety and depressive disorders. Number of expected cases (95% PI) were estimated with negative binomial models, using data from 1 March, 2018 to 29 February, 2020. Shaded areas in blue represent 95 % PI. Vertical lines show 1 March, 2020. The MEDEA deprivation index is calculated at the census tract level in urban areas of Catalonia and is categorised in quintiles of deprivation. It also includes a rural category for individuals living in rural areas. Demographic characteristics of individuals in the study population. Individuals were aged 18 years or older and had at least one year of prior history before index date (1st March 2018) in the database. The MEDEA deprivation index is calculated at the census tract level in urban areas of Catalonia and is Supplementary Figure 2. Expected and observed incidence of anxiety and depressive disorders stratified by age group and gender in primary care in Catalonia -February, 2020.

Supplementary materials
Validation of the modelling approach for incidence rates of anxiety and depressive disorders stratified by age and gender. Number of expected cases (95% PI) were estimated with negative binomial models, using data from 1 March, 2017 to 28 February, 2019. Vertical lines show 1 March, 2019. Shaded areas in blue represent 95 % PI.