Skip to main content

Association of depressive symptoms with incidence and mortality rates of COVID-19 over 2 years among healthcare workers in 20 countries: multi-country serial cross-sectional study

Abstract

Background

Long-term deterioration in the mental health of healthcare workers (HCWs) has been reported during and after the COVID-19 pandemic. Determining the impact of COVID-19 incidence and mortality rates on the mental health of HCWs is essential to prepare for potential new pandemics. This study aimed to investigate the association of COVID-19 incidence and mortality rates with depressive symptoms over 2 years among HCWs in 20 countries during and after the COVID-19 pandemic.

Methods

This was a multi-country serial cross-sectional study using data from the first and second survey waves of the COVID-19 HEalth caRe wOrkErS (HEROES) global study. The HEROES study prospectively collected data from HCWs at various health facilities. The target population included HCWs with both clinical and non-clinical roles. In most countries, healthcare centers were recruited based on convenience sampling.

As an independent variable, daily COVID-19 incidence and mortality rates were calculated using confirmed cases and deaths reported by Johns Hopkins University. These rates represent the average for the 7 days preceding the participants’ response date. The primary outcome was depressive symptoms, assessed by the Patient Health Questionnaire-9. A multilevel linear mixed model (LMM) was conducted to investigate the association of depressive symptoms with the average incidence and mortality rates.

Results

A total of 32,223 responses from the participants who responded to all measures used in this study on either the first or second survey, and on both the first and second surveys in 20 countries were included in the analysis. The mean age was 40.1 (SD = 11.1), and 23,619 responses (73.3%) were from females. The 9323 responses (28.9%) were nurses and 9119 (28.3%) were physicians. LMM showed that the incidence rate was significantly and positively associated with depressive symptoms (coefficient = 0.008, standard error 0.003, p = 0.003). The mortality rate was significantly and positively associated with depressive symptoms (coefficient = 0.049, se = 0.020, p = 0.017).

Conclusions

This is the first study to show an association between COVID-19 incidence and mortality rates with depressive symptoms among HCWs during the first 2 years of the outbreak in multiple countries. This study’s findings indicate that additional mental health support for HCWs was needed when the COVID-19 incidence and mortality rates increase during and after the early phase of the pandemic, and these findings may apply to future pandemics.

Trial registration

Clinicaltrials.gov, NCT04352634.

Peer Review reports

Manuscript

Background

Deterioration in the mental health of healthcare workers (HCWs) has been reported not only in the early stages of the novel coronavirus disease (COVID-19) pandemic [1,2,3,4,5,6,7], but also in the long term during and after the COVID-19 pandemic [8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25]. Mental health problems have been associated with absenteeism, job performance, and intention to leave [26,27,28]. Prevention and countermeasures for mental health problems among HCWs are essential for sustaining the healthcare system. Few studies have examined long-term mental health outcomes among HCWs during and after the COVID-19 pandemic, and the long-term factors related to mental health among HCWs have not been adequately elucidated [13, 29,30,31,32,33]. Therefore, it is necessary to investigate the factors associated with the long-term mental health deterioration among HCWs during and after the COVID-19 pandemic.

The severity of the COVID-19 outbreak, measured by the incidence and mortality rates of COVID-19, has been reported in many cross-sectional studies as a salient factor associated with increased mental health issues [30, 34,35,36,37,38]. The incidence and mortality rates are factors that increase and decrease over time, which warrant data-recording and analytical approaches adequate to model such time-varying nature. In the context of high COVID-19 incidence rates, HCWs may have experienced deterioration in their mental health due to increased workload, changes in work activities, increased risk of infection, worsening of their fear of infection, and limited access to adequate PPE, among other reasons [39]. Similarly, an increase in the COVID-19 related mortality rate may have led to a deterioration in the mental health of HCWs due to traumatic experiences such as the death of a patient or fear of their own death [40, 41]. Determining the impact of COVID-19 incidence and mortality rates on the mental health of HCWs is essential to prepare for future pandemics.

There are some essential gaps in prior studies of COVID-19 incidence and mortality rates and long-term mental health outcomes [42,43,44,45]. First, no study has investigated the long-term and time-related association between the country's mortality rate and HCWs' mental health during and after the COVID-19 pandemic. Second, long-term associations between the incidence rate and mental health have been reported in single-region or country-specific studies but not in a multi-country study. It is unclear whether the long-term association between the incidence or mortality rate and mental health varies by country or between high-income and low-and-middle-income countries. Third, to our knowledge, there was only one study with depressive symptoms as an outcome, which examined the long-term associations between the incidence rate and depressive symptoms. In the previous study in New York City, daily changes in the local and national incidence rate were modestly correlated with changes in daily mean depressive symptoms among HCWs from May 2020 to January 2021 [44]. Fourth, the timeframe of previous studies was only one year at the most during the COVID-19 pandemic. It is not clear whether the incidence and mortality rates were associated with depressive symptoms among HCWs only during the early phase of the COVID-19 pandemic or also long-term during the COVID-19 pandemic. Therefore, further research using long-term and multi-country data is needed to elucidate the direction and strength of the association of COVID-19 incidence and mortality rates with depressive symptoms. This research may reveal changes in depressive symptoms over time and identify when mental health support for HCWs is needed. Therefore, if these research gaps are studied and show long-term associations over two years between incidence rate and mortality rate and depressive symptoms, long-term factors associated with depressive symptoms can be clarified, and prevention and improvement of long-term depressive symptoms from the early stages of an outbreak can be expected. The results of this study may not only help identify when depressive symptoms become severe and the need to take countermeasures, but also help prevent deterioration by identifying when depressive symptoms gradually worsen from their usual state and take preventive measures.

Accordingly, the present study aimed to investigate the association of depressive symptoms with the prior week’s incidence and mortality rates of COVID-19 over 2 years among HCWs in 20 countries during and after the COVID-19 outbreak. In addition, this study aimed to investigate whether the association of depressive symptoms with the incidence and mortality rates of COVID-19 varied by factors such as individual characteristics and whether the participants lived in a high-income or low-and-middle-income country.

Methods

Study design

The design of this study was a multi-country serial cross-sectional study using multiple waves of survey data. Data were retrieved from the first and second surveys of the international COVID-19 HEalth caRe wOrkErS (HEROES) study. The HEROES study aimed to assess the impact of the COVID-19 pandemic on HCWs’ mental health [38, 46, 47]. While the HEROES study used a prospective cohort design that included HCWs from pre-selected healthcare facilities in the participating countries, the present analysis captures the data from the first and second surveys, which occurred during a 2-year timeframe, thus allowing to examine the effect of time-varying stressors linked to the evolving, dynamic nature of the pandemic. The reason for using data from the first and second surveys was to examine the long-term association between the incidence and mortality rates were associated with depressive symptoms from the early stages of the COVID-19 pandemic, as well as because fewer countries conducted surveys after the second survey. The target population included HCWs with both clinical and non-clinical roles (e.g., physicians, nurses, psychologists, dentists, managers, administrative staff, security, cleaning staff) employed in a wide range of healthcare centers (e.g., hospitals, primary care centers, mental health facilities, elderly homes, rehabilitation centers, emergency medical services) in the 24 participating countries. Healthcare centers were recruited in most countries based on convenience sampling; however, facilities were randomly selected in Brazil, Colombia, Japan, and Lebanon. Sites recruited participants by sharing a link to the study’s digital platform via their work email address or the healthcare center’s internal communication system. The recruitment period for the first survey of the HEROES study varied by country, with the first survey taking place after March 2020 and the second survey taking place 6 months after the first survey in the countries where the second survey was conducted. In most countries, recruitment of the first survey took place during the first wave of the pandemic, at a time when vaccines were scarcely available. The reason why the timing of the first survey differed from country to country was that the survey was started in countries that were ready to begin the survey in order to obtain data from the early stages of the COVID-19 pandemic, which was unpredictable and sudden. The second survey was intended to take place 6 months after the first survey, but in some cases the survey took place after 6 months. Participants generally participated in the longitudinal survey, although some participants began participating in the second survey. An online questionnaire based on standardized measures and ad hoc items was employed. The study used a secure platform based on the REDCap model [48, 49], which was designed to ensure data security and quality. This system was created ad hoc to facilitate translation into different languages, making the survey available in all participants’ languages.

This study was approved by accredited Institutional Review Boards (IRBs) in each participating country. This study was reported in accordance with the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement [50, 51].

Participants

The target population included HCWs with both clinical and non-clinical roles employed in a wide range of healthcare centers in the participating countries. The inclusion criteria for potential participants included being of legal age, working in one of the preselected health facilities, working in a health facility that provides care to confirmed or suspected cases of COVID-19, and having an internet connection to answer the online questionnaire.

Measurements

Outcome

Depressive symptoms were measured by the Patient Health Questionnaire-9 (PHQ-9) and used as a continuous score in this study. The PHQ-9 is a self-report questionnaire developed to assess depressive symptoms by rating one’s condition for the previous 2 weeks [52]. PHQ-9 was answered by participants at the first and second surveys. The validity and reliability of the various language versions have been verified [52, 53]. The PHQ-9 is recommended as an international measure of depressive symptoms by the Common Measures in Mental Health Science Initiative [54, 55].

Independent variables

The incidence and mortality rates of COVID-19 were calculated as a moving average of the 7 days prior to each participant’s response date based on the number of daily confirmed COVID-19 cases and deaths in each country as reported by the COVID-19 Data Repository at Johns Hopkins University [56]. The incidence rate was calculated as the number of COVID-19 cases per 100,000 population, and the mortality rate was calculated by dividing the number of recorded deaths by the number of COVID-19 cases and then multiplying by 100 to obtain a percentage [56].

In addition, four researchers (HA, KW, YM, and DN) discussed the time period for calculating the incidence and mortality rates of COVID-19 with reference to the previous studies [38, 44]. There was a possibility of errors in government statistics for any particular day of the week, and the participant might not know the incidence and mortality rates of COVID-19 of their response date. As a result of the discussion, it was decided that the average incidence and mortality rates of COVID-19 from 7 days to 1 day prior to the response date would be more appropriate than the incidence and mortality rates on the response date. Thus, the incidence and mortality rates of COVID-19 in this study were defined as the average incidence and mortality rates of COVID-19 over the 7 days preceding the date the participants responded (Additional file 1: Figure S1 and Figure S2).

Demographic variables

Questions about age, gender, occupation, experience of contact with patients who had suspected or confirmed cases of COVID-19, and history of mental health problems were initially developed in the HEROES study [47]. The occupation was categorized into five groups based on a previous study: physicians, nurses, other HCWs (e.g., psychologist, social worker, radiology technician, laboratory technician, physical therapist, respiratory therapist, dietician, dentist, occupational therapist, speech therapist, first responder, clinical manager/director, social-community mental health worker), ancillary workers (e.g., non-clinical manager, administrator/secretary/admission, patient transportation, food/hospitality, cleaning staff, maintenance staff, security staff), and other [38]. The question about the experience of contact with patients who were suspected or confirmed cases of COVID-19 was “During the past week, have you been close to patients who were suspected or confirmed cases of COVID-19?” and answered with three choices (yes/ no/ I don’t know). The question about the history of mental health problems was “Before the pandemic, did you have a mental health diagnosis?” and answered with three choices (yes/ no/ I don’t know). Based on the DAC List of ODA Recipients for reporting on aid in 2022 and 2023 published by the Organization for Economic Co-operation and Development, participants’ countries of residence were classified as high-income countries (Chile, Czech Republic, Germany, Italy, Japan, Netherlands, Saudi Arabia, Spain, and Uruguay) or low-and-middle-income countries (Argentina, Armenia, Bolivia, Brazil, Colombia, Guatemala, Mexico, Nigeria, Peru, Tunisia, and Venezuela) [57].

Statistical analysis

We analyzed the dataset of responses of the participants who responded to all measures used in this study on either the first or second survey, and on both the first and second surveys. Missing values were considered missing not at random because some variables were not measured in some countries, so all participants in those countries had missing values for those variables, and participants with severe depressive symptoms may have had missing values because they had difficulty responding to the PHQ-9. We did not adopt any imputation methods because missing in the present study should occur not at random and imputations were problematic. Most of the excluded participants in this study responded only once and did not respond to all of the measures used. The participants with missing values might have had systematic differences from the completers. Instead, we compared and reported statistical differences in demographic characteristics. Mann–Whitney U tests for continuous variables, and Pearson’s chi-squared tests for categorical variables were used to compare variables between the analyzed and excluded responses. The average incidence and mortality rates of COVID-19 and PHQ-9 were calculated for the entire period of this study, each year, and each month based on the participants’ response dates, regardless of the first and second surveys. The intraclass correlation coefficients (ICC) within countries for PHQ-9 were calculated for all responses of the analyzed participants.

Multilevel linear mixed models (LMM) were conducted to investigate the association of depressive symptoms with the incidence and mortality rates [58, 59]. The data were conceptualized as a two-level model, consisting of individuals at the first level and countries at the second level. We used restricted maximum likelihood (ReML) estimation for all LMM analyses [60]. Three steps of multilevel LMM were conducted. The first step of the analysis was performed without explanatory variables (null model: Model 1). This model contains a random effect between countries of the intercept for the depressive symptoms and the fixed effect of time. The time effect was coded as a daily-level continuous variable based on the survey period, with January 1, 2020, as 0. The outcome, the PHQ-9, was used at the daily-level based on the participant’s response date. The incidence or mortality rates were entered into Model 2 (crude model). The average incidence and mortality rates of COVID-19 from 7 days to 1 day prior to the response date were also used at the daily-level based on the participant’s response date. The demographic and work-related variables were entered into Model 3 (adjusted model) as the covariates. Akaike information criteria (AIC) was calculated to check the model fit [61,62,63]. In model 3 of the multilevel linear mixed model (LMM), the grand mean centered on age established a useful zero point, and dummy variables were created for gender, occupation, contact with patients who had suspected or confirmed cases of COVID-19, and history of mental health problems. The equation for Model 3 could be explained as follows;

  • Level 1 (individual level)

    $$Y{(depressive\ symptoms:PHQ-9)}_{ijt}={\beta }_{0j}+{\beta }_{1}*{\left(age\right)}_{ijt}+{\beta }_{2}*{\left(gender: female\right)}_{ijt}+{\beta 3*\left(gender: other\,gender\right)}_{ijt}+{\beta }_{4}*{\left(Occupation: nurses\right)}_{ijt}+{\beta }_{5}*{\left(Occupation: other\,healthcare\,workers\right)}_{ijt}+{\beta }_{6}*{\left(Occupation:\,ancillary\,workers\right)}_{ijt}+{\beta }_{7}*{\left(Occupation:\,other\right)}_{ijt}+{\beta }_{8}*{\left(Contact\,with\,patients\,who\,had\,suspected\,or\,confirmed\,cases\,of\,COVID-19:\,yes\right)}_{ijt}+{\beta }_{9}*{\left(Contact\,with\,patients\,who\,had\,suspected\,or\,confirmed\,cases\,of\,COVID-19:\,I\,do\,not\,know\right)}_{ijt}+{\beta }_{10}*{\left(History\,of\,mental\,health\,problems:\,yes\right)}_{ijt}+{\beta }_{11}*{\left(History\,of\,mental\,health\,problems:\,I\,do\,not\,know\right)}_{ijt}+{e}_{ijt}$$
  • Level 2 (country level)

    $${\beta }_{0jt}=\gamma 00*{\left(incidence\,rate\,or\,mortality\,rate\,of\,COVID-19\right)}_{jt}+{\varepsilon }_{0jt}$$
    $${\varepsilon }_{0jt}\sim \left(0,{\sigma }^{2}\right)$$

Note: Yijt means the score of depressive symptoms of individual healthcare worker j who lived in a country i at Time t. ε0jt means a random effect, and eijt means a residual.

Three similar steps of multilevel LMMs were conducted in each of the subgroups classified by gender (men and women), occupation (physicians, nurses, other HCWs, ancillary workers, and others), contact with patients who had suspected or confirmed cases of COVID-19 (yes, and no), history of mental health problems (yes, and no), and high-income and low-and-middle-income countries. As a sensitivity analysis, similar multilevel LMMs were conducted on 30,520 responses, including the 28,817 responses which were responded to either the first or second survey, and only the responses to the first survey from the 1703 participants who responded to both the first and second surveys.

The statistical significance level was set at a p value of less than 0.05. All analyses were conducted by using SPSS version 29.0 J for Windows (SPSS, Tokyo, Japan).

Results

Participants flow

The 37,105 participants from 24 countries (Argentina, Armenia, Belgium, Bolivia, Brazil, Chile, Colombia, Czech Republic, Germany, Guatemala, Italy, Japan, Lebanon, Mexico, Netherlands, Nigeria, Peru, Puerto Rico, Saudi Arabia, South Africa, Spain, Uruguay, Tunisia, and Venezuela) participated in the first survey. The 15,117 participants from 14 countries (Argentina, Armenia, Chile, Colombia, Czech Republic, Guatemala, Japan, Mexico, Netherlands, Nigeria, Peru, Puerto Rico, and Spain, Tunisia) participated in the second survey. From the first and second surveys, 52,222 responses were collected from 24 countries (Argentina, Armenia, Belgium, Bolivia, Brazil, Chile, Colombia, Czech Republic, Germany, Guatemala, Italy, Japan, Lebanon, Mexico, Netherlands, Nigeria, Peru, Puerto Rico, Saudi Arabia, Tunisia, Venezuela, South Africa, Spain, and Uruguay). The 19,999 responses that did not respond to the variables used in this study were excluded. All participants in 4 countries (Belgium, Lebanon, Puerto Rico, and South Africa) were excluded because they did not respond to the variables used in this study, such as PHQ-9. Thus, 32,223 responses from 30,520 participants in 20 countries (Argentina, Armenia, Bolivia, Brazil, Chile, Colombia, Czech Republic, Germany, Guatemala, Italy, Japan, Mexico, Netherlands, Nigeria, Peru, Saudi Arabia, Spain, Tunisia, Uruguay, Venezuela) were included in the analysis. The 28,817 responses were responded to either the first or second survey, and the 1703 participants who responded to both the first and second surveys (Fig. 1).

Fig. 1
figure 1

Participant flow. The design of this study was a multi-country serial cross-sectional study, not a longitudinal study, although data were retrieved from the first and second surveys of the international COVID-19 HEalth caRe wOrkErS (HEROES) study. We analyzed the dataset of responses of the participants who responded to all measures used in this study

Participant characteristics

The mean age was 40.1 years (SD = 11.1); 23,619 responses (73.3%) were from women, 9323 (28.9%) were nurses, and 9119 (28.3%) were physicians (Table 1). The mean PHQ-9 score was 6.0 (SD = 5.5). The ICC within countries for the PHQ-9 was 0.07 among all analyzed participants’ responses.

Table 1 Characteristics of the participants (N = 32,223)

The mean incidence rate of COVID-19 for the entire period of this study was 11.90 (SD = 16.46), 6.77 (SD = 7.78) in 2020, 20.87 (SD = 22.73) in 2021, and 12.00 (SD = 8.80) in 2022 (Table 2). The mean mortality rate of COVID-19 for the entire period of this study was 6.47 (SD = 4.96), 7.70 (SD = 5.32) in 2020, 4.42 (SD = 3.37) in 2021, and 2.28 (SD = 0.01) in 2022. The number of responses from the analyzed participants each month for each country is shown in Table S1 (Additional file 1: Table S1). The earliest response date was April 24, 2020, and the latest was March 5, 2022, and there were participant responses every month during this period.

Table 2 The annual and monthly means of the incidence and mortality rates of COVID-19, and PHQ-9

The results of Mann–Whitney U tests and Pearson’s chi-square tests, gender (p < 0.01), age (the analyzed responses 40.1 years, the excluded responses 38.6 years; p < 0.01), occupation (p < 0.01), the experience of contact with patients who had suspected or confirmed cases of COVID-19 (p < 0.01), history of mental health problems (p < 0.01), country (p < 0.01), the incidence rate of COVID-19 (the analyzed responses 11.9, the excluded responses 22.6; p < 0.01) and the mortality rate of COVID-19 (the analyzed responses 6.5, the excluded responses 5.7; p < 0.01) were significantly different between the analyzed responses (n = 32,223) and the excluded responses (n = 19,999) (Table 3). There was no significant difference in PHQ-9 scores between the analyzed and excluded responses.

Table 3 Demographic and outcome variables in the analyzed responses and the excluded responses

Association between incidence or mortality rate of COVID-19 and depressive symptoms

Multilevel LMM showed that time was significantly negatively associated with depressive symptoms in model 1 (coefficient = − 0.003, se = 0.000, p < 0.001). The incidence rate was significantly positively associated with depressive symptoms in model 2 (coefficient = 0.011, se = 0.003, p < 0.001) and model 3 (coefficient = 0.008, se = 0.003, p = 0.003) (Table 4). In this multilevel LMM, model 3 had the best model fit (AIC = 194,434.531).

Table 4 Association between incidence rate of COVID-19a and depressive symptoms: results of multilevel linear mixed model (N = 32,223)

The mortality rate was significantly positively associated with depressive symptoms in model 2 (coefficient = 0.076, se = 0.022, p < 0.001) and model 3 (coefficient = 0.049, se = 0.020, p = 0.017) (Table 5). In this multilevel LMM, model 3 had the best model fit (AIC = 194,433.259).

Table 5 Association between mortality rate of COVID-19a and depressive symptoms: results of multilevel linear mixed model (N = 32,223)

In subgroup analyses by gender, there was a significant positive association between the incidence rate of COVID-19 and depressive symptoms in models 2 and 3 in both men and women (Additional file 1: Table S2 and Table S3). There was a significant association in women between the mortality rate of COVID-19 and depressive symptoms in model 2 (coefficient = 0.074, se = 0.026, p = 0.004) and model 3 (coefficient = 0.050, se = 0.024, p = 0.035), but no significant association in men (Additional file 1: Table S4 and Table S5).

In subgroup analyses by occupation, there was a significant association between the incidence rate of COVID-19 and depressive symptoms in model 2 (coefficient = 0.013, se = 0.005, p = 0.018) among physicians, model 2 (coefficient = 0.017, se = 0.005, p = 0.002) and model 3 (coefficient = 0.014, se = 0.005, p = 0.006) among nurses, and model 2 (coefficient = 0.014, se = 0.006, p = 0.037) among other HCWs, but not among ancillary workers and other occupations (Additional file 1: Table S6 to Table S10). There was no significant association between the mortality rate of COVID-19 and depressive symptoms among any of the occupations (Additional file 1: Table S11 to Table S15).

In subgroup analyses by the experience of contact with patients who had suspected or confirmed cases of COVID-19, there was a significant association between the incidence rate of COVID-19 and depressive symptoms in model 2 (coefficient = 0.011, se = 0.005, p = 0.032) among participants without the experience of contact with patients who had suspected or confirmed cases of COVID-19, model 2 (coefficient = 0.009, se = 0.004, p = 0.015) and model 3 (coefficient = 0.009, se = 0.004, p = 0.010) among participants with the experience of contact with patients who had suspected or confirmed cases of COVID-19 (Additional file 1: Table S16 and Table S17). There was a significant association between the mortality rate of COVID-19 and depressive symptoms in model 2 (coefficient = 0.087, se = 0.031, p = 0.005) among participants with the experience of contact with patients who had suspected or confirmed cases of COVID-19, but no significant association among participants without the experience of contact with patients who had suspected or confirmed cases of COVID-19 (Additional file 1: Table S18 and Table S19).

In subgroup analyses by the history of mental health problems, there was a significant association between the incidence rate of COVID-19 and depressive symptoms in model 2 (coefficient = 0.012, se = 0.003, p < 0.001) and model 3 (coefficient = 0.010, se = 0.003, p < 0.001) among participants without a history of mental health problems, but no significant association among participants with a history of mental health problems (Additional file 1: Table S20 and Table S21). There was a significant association between the mortality rate of COVID-19 and depressive symptoms in model 2 (coefficient = 0.119, se = 0.052, p = 0.024) among participants with a history of mental health problems, but no significant association among participants without a history of mental health problems (Additional file 1: Table S22 and Table S23).

In subgroup analyses in high-income or low-and-middle-income countries, there was a significant association between the incidence rate of COVID-19 and depressive symptoms in model 2 (coefficient = 0.015, se = 0.006, p = 0.014) in low-and-middle-income countries, but no significant association in high-income countries (Additional file 1: Table S24 and Table S25). There was a significant positive association between the mortality rate of COVID-19 and depressive symptoms in model 2 and model 3 in both high-income and low-and-middle-income countries (Additional file 1: Table S26 and Table S27).

In summary, subgroup analyses showed that as the incidence rate increased, depressive symptoms worsened among men, women, nurses, those who had contact with patients who had suspected or confirmed cases of COVID-19, those without a history of mental health problems, and those living in low-and-middle-income countries. Subgroup analyses also showed that as the mortality rate increased, depressive symptoms worsened among women, those who had contact with patients who had suspected or confirmed cases of COVID-19, those with a history of mental health problems and those living in both high-income and low-and-middle-income countries.

In a sensitivity analysis of 30,520 responses, there was a similarly significant positive association of incidence and mortality rates with depressive symptoms in the multilevel LMM (Additional file 1: Table S28 and Table S29).

Discussion

This study showed the association of prior week incidence and mortality rates with depressive symptoms over 2 years among HCWs in 20 countries during and after the COVID-19 outbreak. In addition, the association of the incidence and mortality rates with depressive symptoms varied by the country and individual characteristics of the participants.

This study found that the increased daily incidence rates were associated with worsened depressive symptoms among HCWs in 20 countries, including both high-income and low-and-middle-income countries, regardless of the country of residence. Table 3 shows that in the best-fitting model 3, an increase of 1 per 100,000 in the mean incidence of COVID-19 from 7 days to 1 day before the survey date in the healthcare worker’s country of residence increased the value of PHQ-9 of the participant by 0.008. In the Netherlands on February 12, 2022, which had the highest incidence rate of COVID-19 of 731.31 during the study period in 20 countries, the PHQ-9 score was calculated to have increased by 5.85 compared to when the incidence rate was 0. These results suggest that although the coefficient was small at 0.008, a high incidence rate would result in a large increase in the PHQ-9 score, which may exceed the cutoff score of 10 for major depressive disorder, and that the incidence rate was a serious associated factor for depressive symptoms among healthcare workers during the COVID-19 pandemic. The association between incidence rate and depressive symptoms among HCWs during the COVID-19 pandemic was consistent with the previous longitudinal study conducted in New York [44]. This study adds to the findings of previous studies that the incidence rate of COVID-19 is associated with depressive symptoms among HCWs not only in a local region but also in 20 countries, including both high-income and low-and-middle-income countries, and that the association persists for as long as 2 years, not just in the early stages of the COVID-19 pandemic. Reasons for the association over 2 years include a higher workload with increased incidence of COVID-19 and exhaustion with a prolonged COVID-19 pandemic, combined with the sense of no longer being valued by organizations or even by patients in an understaffed environment during and after the COVID-19 pandemic [15].

This study showed an association between depressive symptoms and the mortality rate among multinational HCWs during the first 2 years of the COVID-19 pandemic in 20 countries. Table 4 shows that in the best-fitting model 3, an increase of 1% in the mean mortality rate from 7 days to 1 day before the survey date in the participant’s country of residence increased the PHQ-9 score of the participant by 0.049. Some reasons for observing the association even 2 years after the start of the COVID-19 pandemic include moral injury, difficult clinical decisions regarding patients with COVID-19, and repeated traumatic experiences of patient deaths throughout the pandemic [64]. Therefore, this study suggested that mortality rate might be a long-term associated factor in the mental health of healthcare workers internationally during the COVID-19 pandemic.

This study showed a long-term association between the incidence and mortality rates and depressive symptoms and revealed depressive symptoms over time and when mental health support was needed for HCWs during and after the COVID-19 pandemic. The results of Model 1 in Table 3 and Table 4 show that the PHQ-9 score of a participant decreased by 0.002 from 6.444 with each passing day from January 1, 2020. Thus, considering only time, depressive symptoms among HCWs were worst in the early stages of the COVID-19 pandemic and gradually improved over time. In addition, this study showed that the increased daily incidence and mortality rates were associated with worsening depressive symptoms over 2 years. As it has been reported that future trajectories of the incidence and mortality rates can be simulated [65,66,67], this study’s results are important in identifying periods of worsening depressive symptoms among HCWs during and after the COVID-19 pandemic. This study suggested that countermeasures for depressive symptoms were needed in any country at the beginning of the pandemic and at the times of an increase in the incidence and mortality rates during the 2 years from the beginning of the COVID-19 pandemic. Training on infectious diseases and training on mental health, such as psychological first aid, may be an effective countermeasure for healthcare workers before a pandemic begins or in the very early stages [68, 69]. Psychoeducational interventions (e.g., online stress reduction information, behavioral changes including exercise) and psychosocial interventions (e.g., “talk therapy,” crisis lines), which have been reported to be effective in previous studies [70,71,72,73], may be beneficial to the mental health of health care workers as a countermeasure during a pandemic, especially when infection and mortality rates increase. It may also be important for the government and affiliated organizations such as hospitals to plan and implement programs and countermeasures for the mental health of healthcare workers during a pandemic. This study indicated the need for additional mental health care for healthcare workers when the incidence rate and mortality rate of COVID-19 increased during the pandemic and suggests that specific measures may be needed for those whose mental health is particularly affected by the increases. Based on the results of the subgroup analysis, when the incidence and mortality rates increased, mental health countermeasures were especially needed for HCWs who were women, nurses, those who had contact with suspected or confirmed cases of COVID-19, and those who lived in low-and-middle-income countries.

Furthermore, this study’s results may be applicable during future outbreaks of emerging infectious diseases other than COVID-19. It may be helpful to implement mental health countermeasures for HCWs during potential new global health emergency crises, both at the beginning of an epidemic and at the time of an increase in the incidence and mortality rates of emerging infectious diseases. This study also suggested that during outbreaks of emerging infectious diseases other than COVID-19, it may be necessary to identify the characteristics of healthcare workers who are particularly affected in their mental health by the increases in the incidence rate and mortality rate and to implement specific countermeasures for these healthcare workers. In addition, there might have been benefits to the participants in this study, including the fact that they were able to assess their own mental health and the possibility that the results of this study could prevent or improve their mental health in the future.

Limitations

There are several limitations to this study. First, because this study was based on convenience sampling in most countries and recruitment methods differed somewhat from country to country, it is necessary to interpret the findings of this study with attention to the representativeness of the participants. Selection bias due to the methods used to recruit participants may have limited the generalizability of the results. Examples of selection bias could include the fact that the recruitment methods created biases in the medical specialties of the participants and the hospitals in which they worked, resulting in characteristics of the participants that differed from those of the general healthcare workforce during the COVID-19 pandemic. HCWs who experienced a severe worsening of depressive symptoms during the pandemic may not have been able to participate in this study because of severe depressive symptoms or because they had retired from their jobs. The mean PHQ-9 score for participants in this study was 6.0, and 14.2% had a history of mental health problems, indicating some inclusion of HCWs with worsening depressive symptoms, but selection bias may have occurred if those with severe depressive symptoms were unable to participate. Future studies are needed using recruitment methods that do not cause selection bias, such as random sampling. Second, many participants were excluded due to missing values in their responses. As most participants in this study responded only once and missing values in the present study should occur not at random, the implementation of missing value imputation, such as multiple imputation, was problematic and missing values could not be imputed. The bias in the characteristics of the excluded responses, such as incidence rate, mortality rate, demographic variables, and work-related variables, compared to the included responses in this study may have affected the results and could affect the generalizability of the findings. The exclusion of responses from all participants in all four countries due to missing values could affect the generalizability of the findings. To address this limitation, this study conducted several subgroup analyses to investigate the association between incidence and mortality rates and depressive symptoms by demographic variables such as gender, occupation, and high-income and low-and-middle-income countries. Third, because the study applied the incidence and mortality rates of COVID-19 for each country, the rates in the participant’s region might have differed from those in their country if there were significant regional differences. In future studies, examining associations using national and regional level incidence and mortality rates of COVID-19 may provide further understanding of the association with mental health. Fourth, some factors that might have been related to the incidence and mortality rates of COVID-19 and depressive symptoms among HCWs during the COVID-19 pandemic could not be included in the model in the analysis. Although the incidence and mortality rates of COVID-19 were country-level factors, and thus confounding factors were unlikely in their association with depressive symptoms, the model in the analysis might have been more appropriate if variables such as vaccination against COVID-19, which could not be measured because the study was conducted before the development of the vaccine against COVID-19, could have been included [74]. Fifth, participants’ response dates were scattered across approximately all days of the period, although some days had no one who had responded. Although we addressed this limitation by conducting the LMM analyses, it may have affected the results of this study. Sixth, this study examined the association between the incidence rate and mortality rate of COVID-19 and depressive symptoms with time by skillfully creating and analyzing the time variable. Still, the causality cannot be clarified due to the cross-sectional study design, and caution should be needed in interpreting the results. It would be useful to conduct a longitudinal study to investigate the long-term association between the incidence rate and mortality rate and the mental health of healthcare workers during a pandemic in the future. Finally, the ICC for the PHQ-9 was low among all responses of the analyzed participants, which means only 7% of the total variance in depressive symptoms was explained at the country level. The multilevel analysis of country and individual was appropriate because country differences in the random effect were significant in each of the LMM models. Still, factors related to the individual may explain more about HCWs’ mental health during the COVID-19 pandemic than factors related to the country. However, it should be remembered that factors related to individuals also include social determinants of health.

Conclusions

This is the first study to show an association between depressive symptoms and the incidence and mortality rates of COVID-19 among multi-country HCWs during the first 2 years of the pandemic. This study indicates that additional mental health care for HCWs was needed when the incidence and mortality rates of COVID-19 were increasing during and after the early phase of the pandemic. In addition, the association of depressive symptoms with the incidence and mortality rates of COVID-19 varied by the country and individual characteristics of the participants. The study’s findings may apply to future pandemics.

Availability of data and materials

The dataset and statistical/analytical code that support the findings of this study are available from the corresponding author, D.N., upon reasonable request. The dataset is not publicly available due to privacy reasons.

Data Availability

No datasets were generated or analysed during the current study.

Abbreviations

AIC:

Akaike Information Criteria

COVID-19:

Novel coronavirus disease

HCWs:

Healthcare workers

HEROES study:

COVID-19 HEalth caRe wOrkErS study

ICC:

Intraclass correlation coefficients

IRBs:

Institutional Review Boards

LMM:

Linear mixed models

PHQ-9:

Patient Health Questionnaire-9

REDCap:

Research Electronic Data Capture

ReML:

Restricted maximum likelihood

SD:

Standard deviation

STROBE:

Strengthening the Reporting of Observational Studies in Epidemiology

References

  1. Andhavarapu S, Yardi I, Bzhilyanskaya V, et al. Post-traumatic stress in healthcare workers during the COVID-19 pandemic: a systematic review and meta-analysis. Psychiatry Res. 2022;317:114890. https://doi.org/10.1016/j.psychres.2022.114890.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Asaoka H, Koido Y, Kawashima Y, Ikeda M, Miyamoto Y, Nishi D. Post-traumatic stress symptoms among medical rescue workers exposed to COVID-19 in Japan. Psychiatry Clin Neurosci. 2020;74(9):503–5. https://doi.org/10.1111/pcn.13092.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Kisely S, Warren N, McMahon L, Dalais C, Henry I, Siskind D. Occurrence, prevention, and management of the psychological effects of emerging virus outbreaks on healthcare workers: rapid review and meta-analysis. BMJ. 2020;369:m1642. https://doi.org/10.1136/bmj.m1642.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Lai J, Ma S, Wang Y, et al. Factors associated with mental health outcomes among health care workers exposed to coronavirus disease 2019. JAMA Netw Open. 2020;3(3):e203976. https://doi.org/10.1001/jamanetworkopen.2020.3976.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Marvaldi M, Mallet J, Dubertret C, Moro MR, Guessoum SB. Anxiety, depression, trauma-related, and sleep disorders among healthcare workers during the COVID-19 pandemic: a systematic review and meta-analysis. Neurosci Biobehav Rev. 2021;126:252–64. https://doi.org/10.1016/j.neubiorev.2021.03.024.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Pappa S, Ntella V, Giannakas T, Giannakoulis VG, Papoutsi E, Katsaounou P. Prevalence of depression, anxiety, and insomnia among healthcare workers during the COVID-19 pandemic: a systematic review and meta-analysis [published correction appears in Brain Behav Immun. 2021 Feb;92:247]. Brain Behav Immun. 2020;88:901–7. https://doi.org/10.1016/j.bbi.2020.05.026.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Tong J, Zhang J, Zhu N, et al. Effects of COVID-19 pandemic on mental health among frontline healthcare workers: A systematic review and meta-analysis. Front Psychol. 2023;13:1096857. https://doi.org/10.3389/fpsyg.2022.1096857.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Duden GS, Reiter J, Paswerg A, Weibelzahl S. Mental health of healthcare professionals during the ongoing COVID-19 pandemic: a comparative investigation from the first and second pandemic years. BMJ Open. 2023;13(3):e067244. https://doi.org/10.1136/bmjopen-2022-067244.

    Article  PubMed  Google Scholar 

  9. Elliott TR, Perrin PB, Bell AS, Powers MB, Warren AM. Resilience, coping, and distress among healthcare service personnel during the COVID-19 pandemic. BMC Psychiatry. 2021;21(1):489. https://doi.org/10.1186/s12888-021-03506-6. Published 2021 Oct 6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Kapil V, Collett G, Godec T, et al. Longitudinal comparisons of mental health, burnout and well-being in patient-facing, non-patient-facing healthcare professionals and non-healthcare professionals during the COVID-19 pandemic: findings from the CoPE-HCP study. BJPsych Open. 2022;8(5):e173. https://doi.org/10.1192/bjo.2022.579.

    Article  PubMed  Google Scholar 

  11. Kok N, van Gurp J, Teerenstra S, et al. Coronavirus disease 2019 immediately increases burnout symptoms in ICU professionals: a longitudinal cohort study. Crit Care Med. 2021;49(3):419–27. https://doi.org/10.1097/CCM.0000000000004865.

    Article  CAS  PubMed  Google Scholar 

  12. Kotera Y, Liu G, Colman R, et al. A longitudinal study of mental health in healthcare workers in Japan during the initial phase of COVID-19 pandemic: comparison with the general population [published online ahead of print, 2023 Mar 10]. Curr Psychol. 2023;1–4. https://doi.org/10.1007/s12144-023-04444-0

  13. Lee BEC, Ling M, Boyd L, Olsson C, Sheen J. The prevalence of probable mental health disorders among hospital healthcare workers during COVID-19: A systematic review and meta-analysis. J Affect Disord. 2023;330:329–45. https://doi.org/10.1016/j.jad.2023.03.012.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Li W, Frank E, Zhao Z, et al. Mental Health of Young Physicians in China During the Novel Coronavirus Disease 2019 Outbreak. JAMA Netw Open. 2020;3(6):e2010705. https://doi.org/10.1001/jamanetworkopen.2020.10705. Published 2020 Jun 1.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Linzer M, Jin JO, Shah P, et al. Trends in clinician burnout with associated mitigating and aggravating factors during the COVID-19 pandemic [published correction appears in JAMA Health Forum 2023 Feb 3;4(2):e230002]. JAMA Health Forum. 2022;3(11):e224163. https://doi.org/10.1001/jamahealthforum.2022.4163.

    Article  PubMed  PubMed Central  Google Scholar 

  16. López Steinmetz LC, Herrera CR, Fong SB, Godoy JC. A longitudinal study on the changes in mental health of healthcare workers during the COVID-19 pandemic. Psychiatry. 2022;85(1):56–71. https://doi.org/10.1080/00332747.2021.1940469.

    Article  PubMed  Google Scholar 

  17. McGuinness SL, Eades O, Grantham KL, et al. Mental health and wellbeing of health and aged care workers in Australia, May 2021 - June 2022: a longitudinal cohort study. Med J Aust. 2023;218(8):361–7. https://doi.org/10.5694/mja2.51918.

    Article  PubMed  Google Scholar 

  18. Müller MM, Baillès E, Blanch J, et al. Burnout among hospital staff during the COVID-19 pandemic: longitudinal results from the international Cope-Corona survey study. J Psychosom Res. 2023;164:111102. https://doi.org/10.1016/j.jpsychores.2022.111102.

    Article  PubMed  Google Scholar 

  19. Roberts T, Daniels J, Hulme W, et al. Psychological distress and trauma in doctors providing frontline care during the COVID-19 pandemic in the United Kingdom and Ireland: a prospective longitudinal survey cohort study. BMJ Open. 2021;11(7):e049680. https://doi.org/10.1136/bmjopen-2021-049680.

    Article  PubMed  Google Scholar 

  20. Sasaki N, Asaoka H, Kuroda R, Tsuno K, Imamura K, Kawakami N. Sustained poor mental health among healthcare workers in COVID-19 pandemic: a longitudinal analysis of the four-wave panel survey over 8 months in Japan. J Occup Health. 2021;63(1):e12227. https://doi.org/10.1002/1348-9585.12227.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Sasaki N, Kuroda R, Tsuno K, Kawakami N. The deterioration of mental health among healthcare workers during the COVID-19 outbreak: a population-based cohort study of workers in Japan. Scand J Work Environ Health. 2020;46(6):639–44. https://doi.org/10.5271/sjweh.3922.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Serpa ALO, Pinto ALB, Diaz AP, et al. The mental health of Brazilian healthcare professionals during the COVID-19 pandemic: a longitudinal study. Braz J Psychiatry. 2022;44(4):401–8. https://doi.org/10.47626/1516-4446-2021-2347.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Th’ng F, Rao KA, Ge L, et al. A one-year longitudinal study: changes in depression and anxiety in frontline emergency department healthcare workers in the COVID-19 pandemic. Int J Environ Res Public Health. 2021;18(21):11228. https://doi.org/10.3390/ijerph182111228. Published 2021 Oct 26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Umbetkulova S, Kanderzhanova A, Foster F, Stolyarova V, Cobb-Zygadlo D. Mental Health Changes in Healthcare Workers DURING COVID-19 Pandemic: A Systematic REVIEW of Longitudinal Studies [published online ahead of print, 2023 May 4]. Eval Health Prof. 2023;1632787231165076. https://doi.org/10.1177/01632787231165076

  25. Xiong NN, Fan TT, Leonhart R, et al. Workplace factors can predict the stress levels of healthcare workers during the COVID-19 pandemic: First interim results of a multicenter follow-up study. Front Public Health. 2022;10:1002927. https://doi.org/10.3389/fpubh.2022.1002927.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Johnston DA, Harvey SB, Glozier N, Calvo RA, Christensen H, Deady M. The relationship between depression symptoms, absenteeism and presenteeism. J Affect Disord. 2019;256:536–40. https://doi.org/10.1016/j.jad.2019.06.041.

    Article  CAS  PubMed  Google Scholar 

  27. Tawfik DS, Scheid A, Profit J, Shanafelt T, Trockel M, Adair KC, Sexton JB, Ioannidis JPA. Evidence relating health care provider burnout and quality of care: a systematic review and meta-analysis. Ann Intern Med. 2019;171(8):555–67. https://doi.org/10.7326/M19-1152.

    Article  PubMed  PubMed Central  Google Scholar 

  28. van Rijn RM, Robroek SJ, Brouwer S, Burdorf A. Influence of poor health on exit from paid employment: a systematic review. Occup Environ Med. 2014;71(4):295–301. https://doi.org/10.1136/oemed-2013-101591.

    Article  PubMed  Google Scholar 

  29. Kunzler AM, Röthke N, Günthner L, et al. Mental burden and its risk and protective factors during the early phase of the SARS-CoV-2 pandemic: systematic review and meta-analyses. Global Health. 2021;17(1):34. https://doi.org/10.1186/s12992-021-00670-y.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Tian H, Qiao T, Teng J, et al. Factors associated with depression among healthcare workers during the COVID-19 pandemic: a systematic review and meta-analysis [published online ahead of print, 2023 Sep 15]. Psychol Med. 2023;1–10. https://doi.org/10.1017/S0033291723002271

  31. Manchia M, Gathier AW, Yapici-Eser H, et al. The impact of the prolonged COVID-19 pandemic on stress resilience and mental health: a critical review across waves. Eur Neuropsychopharmacol. 2022;55:22–83. https://doi.org/10.1016/j.euroneuro.2021.10.864.

    Article  CAS  PubMed  Google Scholar 

  32. Muller AE, Hafstad EV, Himmels JPW, et al. The mental health impact of the covid-19 pandemic on healthcare workers, and interventions to help them: A rapid systematic review. Psychiatry Res. 2020;293:113441. https://doi.org/10.1016/j.psychres.2020.113441.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Witteveen AB, Young SY, Cuijpers P, et al. COVID-19 and common mental health symptoms in the early phase of the pandemic: An umbrella review of the evidence. PLoS Med. 2023;20(4):e1004206. https://doi.org/10.1371/journal.pmed.1004206. Published 2023 Apr 25.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Abdessater M, Rouprêt M, Misrai V, et al. COVID19 pandemic impacts on anxiety of French urologist in training: outcomes from a national survey. Prog Urol. 2020;30(8–9):448–55. https://doi.org/10.1016/j.purol.2020.04.015.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Khusid JA, Weinstein CS, Becerra AZ, et al. Well-being and education of urology residents during the COVID-19 pandemic: results of an American national survey. Int J Clin Pract. 2020;74(9):e13559. https://doi.org/10.1111/ijcp.13559.

    Article  CAS  PubMed  Google Scholar 

  36. Prati G, Mancini AD. The psychological impact of COVID-19 pandemic lockdowns: a review and meta-analysis of longitudinal studies and natural experiments. Psychol Med. 2021;51(2):201–11. https://doi.org/10.1017/S0033291721000015.

    Article  PubMed  Google Scholar 

  37. Xu J, Xu QH, Wang CM, Wang J. Psychological status of surgical staff during the COVID-19 outbreak. Psychiatry Res. 2020;288:112955. https://doi.org/10.1016/j.psychres.2020.112955.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Czepiel D, McCormack C, da Silva ATC, et al. Inequality on the frontline: A multi-country study on gender differences in mental health among healthcare workers during the COVID-19 pandemic. Glob Ment Health (Camb). 2024;11:e34. https://doi.org/10.1017/gmh.2024.18.

    Article  PubMed  Google Scholar 

  39. Balai MK, Avasthi RD, Va R, Jonwal A. Psychological impacts among health care personnel during COVID-19 pandemic: a systematic review. J Caring Sci. 2022;11(2):118–25. https://doi.org/10.34172/jcs.2022.14.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Asaoka H, Koido Y, Kawashima Y, Ikeda M, Miyamoto Y, Nishi D. Association between clinical decision for patients with COVID-19 and post-traumatic stress symptoms among healthcare professionals during the COVID-19 pandemic. Environ Occup Health Pract. 2022;425;(1). https://doi.org/10.1539/eohp.2022-0018-OA

  41. Greenberg N, Docherty M, Gnanapragasam S, Wessely S. Managing mental health challenges faced by healthcare workers during covid-19 pandemic. BMJ. 2020;368:m1211. https://doi.org/10.1136/bmj.m1211.

    Article  PubMed  Google Scholar 

  42. Rosenström T, Tuisku K, Suvisaari J, et al. Healthcare workers’ heterogeneous mental-health responses to prolonging COVID-19 pandemic: a full year of monthly follow up in Finland. BMC Psychiatry. 2022;22(1):724. https://doi.org/10.1186/s12888-022-04389-x.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Maunder RG, Heeney ND, Hunter JJ, et al. Trends in burnout and psychological distress in hospital staff over 12 months of the COVID-19 pandemic: a prospective longitudinal survey. J Occup Med Toxicol. 2022;17(1):11. https://doi.org/10.1186/s12995-022-00352-4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Zebley B, Wolk D, McAllister M, Lynch CJ, Mikofsky R, Liston C. Individual differences in the affective response to pandemic-related stressors in COVID-19 health care workers. Biol Psychiatry Glob Open Sci. 2021;1(4):336–44. https://doi.org/10.1016/j.bpsgos.2021.08.008.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Pan CX, Crupi R, August P, et al. An intensive longitudinal assessment approach to surveilling trajectories of burnout over the first year of the COVID pandemic. Int J Environ Res Public Health. 2023;20(4):2930. https://doi.org/10.3390/ijerph20042930.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Mascayano F, van der Ven E, Moro MF, et al. The impact of the COVID-19 pandemic on the mental health of healthcare workers: study protocol for the COVID-19 HEalth caRe wOrkErS (HEROES) study. Soc Psychiatry Psychiatr Epidemiol. 2022;57(3):633–45. https://doi.org/10.1007/s00127-021-02211-9.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Rivera-Segarra E, Mascayano F, Alnasser L, et al. Global mental health research and practice: a decolonial approach. Lancet Psychiatry. 2022;9(7):595–600. https://doi.org/10.1016/S2215-0366(22)00043-8.

    Article  PubMed  Google Scholar 

  48. Harris PA, Taylor R, Minor BL, et al. The REDCap consortium: Building an international community of software platform partners. J Biomed Inform. 2019;95:103208. https://doi.org/10.1016/j.jbi.2019.103208.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Harris PA, Taylor R, Thielke R, Payne J, Gonzalez N, Conde JG. Research electronic data capture (REDCap)–a metadata-driven methodology and workflow process for providing translational research informatics support. J Biomed Inform. 2009;42(2):377–81. https://doi.org/10.1016/j.jbi.2008.08.010.

    Article  PubMed  Google Scholar 

  50. Ghaferi AA, Schwartz TA, Pawlik TM. STROBE reporting guidelines for observational studies. JAMA Surg. 2021;156(6):577–8. https://doi.org/10.1001/jamasurg.2021.0528.

    Article  PubMed  Google Scholar 

  51. Vandenbroucke JP, von Elm E, Altman DG, et al. Strengthening the Reporting of Observational Studies in Epidemiology (STROBE): explanation and elaboration. PLoS Med. 2007;4(10):e297. https://doi.org/10.1371/journal.pmed.0040297.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Kroenke K, Spitzer RL, Williams JB. The PHQ-9: validity of a brief depression severity measure. J Gen Intern Med. 2001;16(9):606–13. https://doi.org/10.1046/j.1525-1497.2001.016009606.x.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Muramatsu K, Miyaoka H, Kamijima K, et al. Performance of the Japanese version of the Patient Health Questionnaire-9 (J-PHQ-9) for depression in primary care. Gen Hosp Psychiatry. 2018;52:64–9. https://doi.org/10.1016/j.genhosppsych.2018.03.007.

    Article  PubMed  Google Scholar 

  54. Farber GK, Gage S, Kemmer D. A Collaborative effort to establish common metrics for use in mental health [published online ahead of print, 2023 Aug 2]. JAMA Psychiatry. 2023. https://doi.org/10.1001/jamapsychiatry.2023.2282.

  55. Farber GK, Gage S, Kemmer D, White R. Common measures in mental health: a joint initiative by funders and journals. Lancet Psychiatry. 2023;10(6):465–70. https://doi.org/10.1016/S2215-0366(23)00139-6.

    Article  PubMed  Google Scholar 

  56. Dong E, Du H, Gardner L. An interactive web-based dashboard to track COVID-19 in real time [published correction appears in Lancet Infect Dis. 2020 Sep;20(9):e215]. Lancet Infect Dis. 2020;20(5):533–4. https://doi.org/10.1016/S1473-3099(20)30120-1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. The Organisation for Economic Co-operation and Development, the DAC List of ODA Recipients for reporting on aid in 2022 and 2023, 2023. https://www.oecd.org/dac/financing-sustainable-development/development-finance-standards/DAC-List-of-ODA-Recipients-for-reporting-2022-23-flows.pdf. Accessed 17 Apr 2024.

  58. Laird NM, Ware JH. Random-effects models for longitudinal data. Biometrics. 1982;38(4):963–74.

    Article  CAS  PubMed  Google Scholar 

  59. Ryoo JH. Model selection with the linear mixed model for longitudinal data. Multivariate Behav Res. 2011;46(4):598–624. https://doi.org/10.1080/00273171.2011.589264.

    Article  PubMed  Google Scholar 

  60. Tak S, Ye JC. Statistical analysis of fNIRS data: a comprehensive review. Neuroimage. 2014;85(Pt 1):72–91. https://doi.org/10.1016/j.neuroimage.2013.06.016.

    Article  PubMed  Google Scholar 

  61. Akaike H. Information theory and an extension of the maximum likelihood principle. In: Petrov BN and Csaki F, eds. International Symposium on Information Theory, 1973;267–281.

  62. Vaida F, Blanchard S. Conditional Akaike information for mixed-effects models. Biometrika. 2005;92:351–70.

    Article  Google Scholar 

  63. Wenren C, Shang J, Pan J. Marginal conceptual predictive statistic for mixed model selection. Open J Stat. 2016;6:239–53. https://doi.org/10.4236/ojs.2016.62021.

    Article  Google Scholar 

  64. Portillo-Van Diest A, Vilagut G, Alayo I, et al. Traumatic stress symptoms among Spanish healthcare workers during the COVID-19 pandemic: a prospective study. Epidemiol Psychiatr Sci. 2023;32:e50. https://doi.org/10.1017/S2045796023000628.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Kissler SM, Tedijanto C, Goldstein E, Grad YH, Lipsitch M. Projecting the transmission dynamics of SARS-CoV-2 through the postpandemic period. Science. 2020;368(6493):860–8. https://doi.org/10.1126/science.abb5793.

    Article  CAS  PubMed  Google Scholar 

  66. Singh DE, Marinescu MC, Guzmán-Merino M, et al. Simulation of COVID-19 propagation scenarios in the Madrid metropolitan area [published correction appears in Front Public Health. 2023 Mar 16;11:1180932]. Front Public Health. 2021;9:636023. https://doi.org/10.3389/fpubh.2021.636023. Published 2021 Mar 16.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Yajada M, Karimi Moridani M, Rasouli S. Mathematical model to predict COVID-19 mortality rate. Infect Dis Model. 2022;7(4):761–76. https://doi.org/10.1016/j.idm.2022.11.005.

    Article  PubMed  PubMed Central  Google Scholar 

  68. Asaoka H, Koido Y, Kawashima Y, Ikeda M, Miyamoto Y, Nishi D. Longitudinal change of psychological distress among healthcare professionals with and without psychological first aid training experience during the COVID-19 pandemic. Int J Environ Res Public Health. 2021;18(23):12474. https://doi.org/10.3390/ijerph182312474. Published 2021 Nov 26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. De Kock JH, Latham HA, Leslie SJ, Grindle M, Munoz SA, Ellis L, Polson R, O’Malley CM. A rapid review of the impact of COVID-19 on the mental health of healthcare workers: implications for supporting psychological well-being. BMC Public Health. 2021;21(1):104. https://doi.org/10.1186/s12889-020-10070-3. Published 2021 Jan 9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Boucher VG, Haight BL, Hives BA, Zumbo BD, Merali-Dewji A, Hutton S, Liu Y, Nguyen S, Beauchamp MR, Black AT, Puterman E. Effects of 12 weeks of at-home, application-based exercise on health care workers' depressive symptoms, burnout, and absenteeism: a randomized clinical trial [published online ahead of print, 2023 Aug 9]. JAMA Psychiatry. 2023;e232706. https://doi.org/10.1001/jamapsychiatry.2023.2706

  71. Buselli R, Corsi M, Veltri A, Baldanzi S, Chiumiento M, Lupo ED, Marino R, Necciari G, Caldi F, Foddis R, Guglielmi G, Cristaudo A. Mental health of Health Care Workers (HCWs): a review of organizational interventions put in place by local institutions to cope with new psychosocial challenges resulting from COVID-19. Psychiatry Res. 2021;299:113847. https://doi.org/10.1016/j.psychres.2021.113847.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Otte C, Gold SM, Penninx BW, Pariante CM, Etkin A, Fava M, Mohr DC, Schatzberg AF. Major depressive disorder. Nat Rev Dis Primers. 2016;2:16065. https://doi.org/10.1038/nrdp.2016.65. Published 2016 Sep 15.

    Article  PubMed  Google Scholar 

  73. World Health Organization. WHO guidelines on mental health at work. September 28, 2022. https://www.who.int/publications/i/item/9789240053052. Accessed 9 Jul 2024.

  74. Asaoka H, Koido Y, Kawashima Y, Ikeda M, Miyamoto Y, Nishi D. Longitudinal change in depressive symptoms among healthcare professionals with and without COVID-19 vaccine hesitancy from October 2020 to June 2021 in Japan. Ind Health. 2022;60(4):387–94. https://doi.org/10.2486/indhealth.2021-016.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

The members of the HEROES study group were listed in Table S30. We also acknowledge all HCWs, all hospitals, and healthcare institutions that participated in the survey.

Funding

This work was supported by the Japan Society for the Promotion of Science KAKENHI Grant Number JP23KJ0461 and The Health Care Science Institute Research Grant to Hiroki Asaoka, by Health and Labor Sciences Research Grants (19IA2014) to Daisuke Nishi, and by the Ministry of Health of the Czech Republic (NU22J-09–00064) to Dominika Seblova. The COVID-19 HEROES study in Italy was funded by the Italian Ministry of Education, Universities and Research (Bando FISR 2020IP_05308) and Fondazione di Sardegna (Bando 2020). The Colombia Heroes project is funded by the MInciencias (grant #: CT 860-2020), and by the UK Research and Innovation (UKRI) (grant #: ES/V013157/1). The funding organizations had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.

Author information

Authors and Affiliations

Authors

Consortia

Contributions

HA had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. Concept and design: HA, KW, YM, DN. Acquisition, analysis, or interpretation of data: HA, KW, YM, AR, EV, MF, LA, OA, AB1. AB2, SD, ME, EF, FF, LG, HH, RE, JL, HM, GM, RM, CM, JN, UO, AB3, EC, JR, AM, DS, AT, LV, OG, DB, MG, AI, AJ, JS, MT, RA, ES, FM, DN. Drafting of the manuscript: HA, KW, YM, MF, DN. Critical review of the manuscript for important intellectual content: HA, KW, YM, AR, EV, MF, LA, OA, AB1. AB2, SD, ME, EF, FF, LG, HH, RE, JL, HM, GM, RM, CM, JN, UO, AB3, EC, JR, AM, DS, AT, LV, OG, DB, MG, AI, AJ, JS, MT, RA, ES, FM, DN. Statistical analysis: HA, KW, DN. Obtained funding: HA, DN. Administrative, technical, or material support: HA, DN. Supervision: DN. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Daisuke Nishi.

Ethics declarations

Ethics approval and consent to participate

This study was approved by accredited Institutional Review Boards (IRBs) in each participating country, including the IRB of the Faculty of Medicine, University of Chile; Columbia University Medical Center IRB; Hospital La Paz Ethics Committee; Research Ethics Committee of the Graduate School of Medicine and Faculty of Medicine at the University of Tokyo (No. 2019164NI-(1)(2)(3)(4)) and the research ethics committee of the National Hospital Organization Disaster Medical Center (No. 2019–19) in Japan; National Committee of Ethics in Research (CONEP/CNS/MS); Ethics Committee of the Azienda Mista Ospedaliero-Universitaria of Cagliari; University of Ibadan/University College Hospital Ethics Review Committee; Instituto Jalisciense de Salud Mental (SALME); Comité de Ética de Investigación en Seres Humanos, de la Facultad de Medicina, Universidad de Chile; IRB University of Applied Sciences Emden / Leer; Comité de Bioética, de la Facultad de Medicina, Universidad Mayor de San Simón; Ponce Health Sciences University; Ethics Committee of the Czech Ministry of Health, Ethics Committee at the Second Faculty of Medicine, Charles University in Prague; Institutional Ethics Committee at the Institute of Molecular Biology of National Academy of Sciences of Armenia; King Abdullah International Medical Research Center (KAIMRC) IRB committee; Comité de Bioética en Investigación en Salud, Universidad de San Carlos de Guatemala; Bioethics Committee of the National School of Public Health, University of Antioquia and Pan-American Health Organization Ethics Review Committee. In the Netherlands, ethical approval was waived by the Medical Ethical boards of the Maastricht University Medical Center and the Amsterdam Medical Center, because the participants were not regarded as patients or identifiable persons sharing sensitive information following the Medical Research Involving Human Subjects Act (WMO). All participants provided informed consent to participate.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

12916_2024_3585_MOESM1_ESM.docx

Additional file 1: FigureS1, S2- [Trends in the mean of the incidence and mortality rates in each country]. Table S1- [Number of responses in each month for each country]. Table S2-S29- [Results of multilevel LMMs were conducted in each of the subgroups classified by gender, occupation, contact with patients who had suspected or confirmed cases of COVID-19, history of mental health problems, and high-income and low-and-middle-income countries]. Table S30- [The members of the HEROES study group].

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, 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 you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. 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-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Asaoka, H., Watanabe, K., Miyamoto, Y. et al. Association of depressive symptoms with incidence and mortality rates of COVID-19 over 2 years among healthcare workers in 20 countries: multi-country serial cross-sectional study. BMC Med 22, 386 (2024). https://doi.org/10.1186/s12916-024-03585-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12916-024-03585-8

Keywords