Introduction

Global spread of the novel coronavirus SARS-CoV-2, the causative agent of coronavirus disease 2019 (COVID-19), has created an unprecedented infectious disease crisis worldwide [1, 2]. In the current absence of a vaccine or curative treatment, the search for an effective treatment for patients with COVID-19 among all available medications is urgently needed [3, 4].

Both COVID-19 and severe acute respiratory syndrome (SARS) are characterized by an overexuberant inflammatory response [5] and, for COVID-19, viral load is associated with the worsening of symptoms [6].

Accumulating evidence suggests that severe COVID-19 is associated with an increased plasma level of inflammatory mediators including cytokines and chemokines such as interleukin (IL)-2, IL-6, IL-7, IL-10, tumor necrosis factor alpha (TNF-α), monocyte chemoattractant protein-1 (MCP1; also known as CCL2), macrophage inflammatory protein 1 alpha (MIP1α; also known as CCL3), C-reactive protein, ferritin, and D-dimers [7]. A recent meta-analysis [8] of studies conducted in individuals with major depressive disorder following antidepressant treatment, mostly including selective serotonin reuptake inhibitors (SSRIs), supports that, overall, antidepressants may be associated with decreased plasma levels of 4 of 16 tested inflammatory mediators, including IL-10, TNF-α, and CCL-2, which are associated with COVID-19 severity [7], as well as IL-6, which is highly correlated with disease mortality [7, 9].

Furthermore, a recent study [10] suggests that functional inhibitors of acid sphingomyelinase activity (FIASMA), including several SSRI and non-SSRI antidepressants, may prevent the infection of epithelial cells with SARS-CoV-2 [11].

A recent in-vitro study [12] also supports antiviral effects of the SSRI fluoxetine on SARS-CoV-2, although this effect was not observed for other SSRIs, including paroxetine and escitalopram.

Finally, a recent randomized trial suggests that fluvoxamine, a SSRI, sigma-1 receptor (S1R) agonist and FIASMA, may prevent clinical deterioration in outpatients with acute COVID-19 compared to placebo [13].

In this context, we hypothesized that antidepressants could be potentially useful in reducing the risk of intubation or death in patients with COVID-19. Short-term use of low to moderate doses of antidepressants, and particularly of SSRIs, is generally well tolerated [14], and notably in older adults [15,16,17].

To our knowledge, no clinical study has examined to date the potential usefulness of antidepressants in patients hospitalized for COVID-19. Observational studies of patients with COVID-19 taking medications for other indications can help decide which treatment should be prioritized for randomized clinical trials and minimize the risk for patients of being exposed to potentially harmful and ineffective treatments.

We took advantage of the Assistance Publique-Hôpitaux de Paris (AP-HP) Health Data Warehouse, which includes data on all patients with COVID-19 who had been consecutively admitted to any of the 39 AP-HP Greater Paris University hospitals.

In this report, we examined the association between antidepressant use and the risk of intubation or death among adult patients who have been admitted to these medical centers with COVID-19. If a significant protective association were found, we sought to perform exploratory analyses to examine whether this association may be specific to certain classes of antidepressants (i.e. SSRIs, serotonin–norepinephrine reuptake inhibitors (SNRIs), tricyclic, tetracyclic and α2-antagonist antidepressants) or certain individual antidepressant medications, in order to help guide future studies. Our main hypothesis was that antidepressant use would be associated with reduced risk of intubation or death in time-to-event analyses adjusting for patient characteristics, clinical and biological markers of disease severity, and other psychotropic medications.

Methods

Setting

We conducted a multicenter observational retrospective study at AP-HP, which includes 39 hospitals of which 23 are acute, 20 are adult and 3 are pediatric hospitals. We included all adults aged 18 years or over who have been admitted with COVID-19 to these medical centers from the beginning of the epidemic in France, i.e. January 24th, until April 1st. COVID-19 was ascertained by a positive reverse-transcriptase–polymerase-chain-reaction (RT-PCR) test from analysis of nasopharyngeal or oropharyngeal swab specimens. This observational study using routinely collected data received approval from the Institutional Review Board of the AP-HP clinical data warehouse (decision CSE-20-20_COVID19, IRB00011591, April 8th, 2020). AP-HP clinical Data Warehouse initiative ensures patient information and informed consent regarding the different approved studies through a transparency portal in accordance with European Regulation on data protection and authorization n°1980120 from National Commission for Information Technology and Civil Liberties (CNIL). All procedures related to this work adhered to the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008.

Data sources

We used data from the AP-HP Health Data Warehouse (‘Entrepôt de Données de Santé (EDS)’). This warehouse contains all the clinical data available on all inpatient visits for COVID-19 to any of the 39 AP-HP Greater Paris University hospitals. The data obtained included patient demographic characteristics, vital signs, laboratory test and RT-PCR test results, medication administration data, medication lists during current and past hospitalizations in AP-HP hospitals, current diagnoses, discharge disposition, ventilator use data, and death certificates.

Variables assessed

We obtained the following data for each patient at the time of the hospitalization through electronic health records:[17,18,19] sex, age, hospital, obesity, self-reported current smoking status, any medical condition associated with increased risk of severe COVID-19 [17, 19,20,21], any medication prescribed according to compassionate use or as part of a clinical trial [22], clinical and biological markers of disease severity at admission [23, 24], any current mood or anxiety disorder, any other current psychiatric disorder, and any psychotropic medication, including any antidepressant, benzodiazepine or Z-drug, mood stabilizer (i.e. lithium or antiepileptic medications with mood stabilizing effects), and antipsychotic medication [25, 26]. These variables are detailed in Supplementary Text 1.

Antidepressant use

Study baseline was defined as the date of hospital admission. Antidepressant use was defined as receiving any antidepressant within the first 48 h of hospital admission and before the end of the index hospitalization or intubation or death. We used this delay because we considered that, in a context of overwhelming of all hospital units during the COVID-19 peak incidence, patients may not have received or been prescribed the treatment the first day of their admission.

Primary endpoint

The primary endpoint was the time from study baseline to intubation or death. For patients who died after intubation, the timing of the primary endpoint was defined as the time of intubation. Patients without an end-point event had their data censored on May 20th, 2020.

Statistical analysis

We calculated frequencies and means (±standard deviations (SD)) of each baseline characteristic described above in patients receiving or not receiving antidepressants and compared them using chi-square tests or Welch’s t-tests.

To examine the associations between the use of any antidepressant, each class of antidepressants, and each individual antidepressant with the composite endpoint of intubation or death, we performed Cox proportional-hazards regression models. To help account for the nonrandomized prescription of antidepressants and reduce the effects of confounding, the primary analyses used a propensity score analysis with inverse probability weighting [27, 28]. The individual propensities for exposures were estimated by multivariable logistic regression models that included sex, age, hospital, obesity, smoking status, any medical condition, any medication prescribed according to compassionate use or as part of a clinical trial, clinical and biological markers of severity of COVID-19, the presence of mood or anxiety or other current psychiatric disorders, and any prescribed benzodiazepine or Z-drug, mood stabilizer, and antipsychotic medication. In the inverse-probability-weighted analyses, the predicted probabilities from the propensity-score models were used to calculate the stabilized inverse-probability-weighting weights [27]. Associations of any antidepressant, each class of antidepressants, and each individual antidepressant with the primary endpoint were then estimated using multivariable Cox regression models including the inverse-probability-weighting weights. Kaplan–Meier curves were performed using the inverse-probability-weighting weights [29], and their pointwise 95% confidence intervals were estimated using the nonparametric bootstrap method [30].

We conducted sensitivity analyses, including multivariable Cox regression models comprising as covariates the same variables as the inverse-probability-weighted analyses, and univariate Cox regression models in matched analytic samples. For these latter analyses, we decided a priori to select one control for each exposed case for exposures to any antidepressant and each class of antidepressants, and two controls for each exposed case for individual antidepressant medications, based on the same variables used for both the inverse-probability-weighted and the multivariable Cox regression analyses. Weighted Cox regression models were used when proportional hazards assumption was not met. To reduce the effects of confounding, optimal matching was used in order to obtain the smallest average absolute distance across all clinical characteristics between exposed patients and non-exposed matched controls.

We also performed several additional analyses. First, to increase our confidence that the results might not be due to unmeasured confounding or indication bias, we examined whether the primary endpoint differed between patients receiving an antidepressant only in the 3 months before hospital admission (during a prior hospitalization in a AP-HP hospital) and those who received it during the visit only. Second, to examine a potential immortal bias in the exposed group, we performed additional Cox regression analyses to compare the potential effect of antidepressant use to that of an active comparator, i.e., the benzodiazepine diazepam. Third, we performed multivariable Cox regression models including interaction terms to examine whether the association between antidepressant use and the endpoint significantly differed across subgroups defined by baseline characteristics. Fourth, we examined a potential dose-effect relationship by testing the association between the daily dose received (converted into fluoxetine-equivalent dose [31] and dichotomized at the median value) with the endpoint among patients receiving antidepressants. Fifth, we examined whether exposure to a combination of antidepressants was associated with a significantly different risk of intubation or death than exposure to only one antidepressant. Finally, we reproduced the main analyses (i) among patients with critical COVID-19 hospitalized in intensive care units (ICUs) and (ii) using death and intubation as separate endpoints.

For all associations, we performed residual analyses to assess the fit of the data, check assumptions, including proportional hazards assumptions, and examined the potential influence of outliers. We followed the recommendations of The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) Initiative [32]. Our main analysis focused on the association between antidepressant use and the primary endpoint. Statistical significance was fixed a priori at two-sided p value < 0.05. Only if a significant protective association were found, we planned to perform exploratory analyses and examine the associations between each class of antidepressants and each individual antidepressant medication with the endpoint. All analyses were conducted in R software version 2.4.3 (R Project for Statistical Computing).

Results

Characteristics of the cohort

Of the 9509 patients consecutively admitted to 33 participating AP-HP Greater Paris University hospitals with COVID-19, 2164 patients (22.8%) were excluded because of missing data or their young age (i.e. less than 18 years of age). In addition, 115 patients who were prescribed an antidepressant were excluded because they received it more than 48 h from hospital admission or after intubation. Of the remaining 7230 adult inpatients, 345 patients (4.8%) received an antidepressant within 48 h of hospital admission at a mean fluoxetine-equivalent dose of 21.6 mg (SD = 14.1) per day (Fig. 1). Doses of each antidepressant medication are described in Supplementary Table 1. Among patients exposed to antidepressants, 298 (86.4%) were exposed to only one antidepressant, and 47 (15.8%) were exposed to at least two antidepressants.

Fig. 1: Study cohort.
figure 1

Antidepressant use was defined as receiving any antidepressant during the first 48 h of hospital admission and before the end of the index hospitalization or intubation or death. SSRIs selective serotonin reuptake inhibitors, SNRIs serotonin–norepinephrine reuptake inhibitors.

First positive COVID-19 RT-PCR tests were obtained after a median delay of 1 day (SD = 12.0) from study baseline. This delay was not significantly different between patients receiving or not receiving antidepressants [median in the exposed group = 1 day (SD = 14.1); median in the non-exposed group = 1 day (SD = 11.9); Mood’s median test Chi-square = 0.21, p = 0.650)].

Over a mean follow-up of 18.4 days (SD = 27.2; median = 5 days; range: 1–117 days), 1272 patients (18.1%) had a primary end-point event at the time of data cutoff on May 20th. Times to follow-up by antidepressant medication are shown in Supplementary Table 2.

All baseline characteristics, when examined independently, were significantly associated with the endpoint (Table 1). In the multivariable analysis, these associations remained significant, except for current smoking, any medication prescribed according to compassionate use or as part of a clinical trial, any current mood or anxiety disorder, any other current psychiatric disorder, any mood stabilizer medication, and any antipsychotic medication (Table 1).

Table 1 Associations of baseline clinical characteristics with the composite endpoint of intubation or death in the cohort of adult patients who had been admitted to the hospital for COVID-19 (n = 7230).

The distributions of patient characteristics according to antidepressant use, and exposure by class of antidepressants and by individual medication are shown in Supplementary Tables 313. In the full sample, antidepressant use and exposure to each antidepressant class and each individual medication significantly differed according to most baseline characteristics. The direction of these associations indicated older age [mean age of patients exposed to antidepressants was 73.7 (SD = 15.8) years, whereas it was 56.8 (SD = 19.3) years in those who were not] and overall greater medical severity of patients who received antidepressants than those who did not. After applying the propensity score weights, these differences were substantially reduced (Supplementary Tables 313). In the matched analytic samples, there were no significant differences in patient characteristics across the different exposures, except for exposures to any antidepressant, non-SSRI antidepressants and SNRI antidepressants, for which sex ratio significantly differed (Supplementary Tables 3–13).

Study endpoint

Among patients receiving any antidepressant, SSRIs, and non-SSRI antidepressants, the primary endpoint of intubation or death occurred respectively in 84 patients (24.3%), 45 patients (23.1%) and 39 patients (26.0%), while 1188 non-exposed patients (17.3%) had this outcome (Table 2). Unadjusted hazard ratio age-stratified estimates of the association between antidepressant use and the endpoint were non-significant (all p > 0.060) (Supplementary Table 14). Following adjustments, the primary multivariable analyses with inverse probability weighting showed significant associations between use of any antidepressant (HR, 0.56; 95% CI, 0.43–0.73, p < 0.001) (Fig. 2), SSRI antidepressants (HR, 0.51; 95% CI, 0.316–0.72, p < 0.001), and non-SSRI antidepressants (HR, 0.65; 95% CI, 0.45–0.93, p = 0.018) (Fig. 3) with reduced risk of intubation or death (Table 2). Among non-SSRI antidepressants, there were no significant associations between the endpoint and exposures to each class of antidepressants, except for α2-antagonist antidepressants (HR, 0.22; 95% CI, 0.05–0.93, p = 0.040) (Table 2; Fig. 3). Finally, exposure to escitalopram, fluoxetine, paroxetine, venlafaxine or mirtazapine was significantly associated with reduced risk of intubation or death (all p < 0.05). There were no significant differences in this risk according to exposure to any other molecule (Table 2; Supplementary Figs. 15).

Table 2 Associations of the use of any antidepressant, each antidepressant class and each individual treatment with the composite endpoint of intubation or death in the full sample and in the matched analytic samples.
Fig. 2: Association between antidepressant use and the composite endpoint of intubation or death.
figure 2

Kaplan–Meier curves of the composite endpoint of intubation or death in the full sample (n = 7230) (A) and in the matched analytic sample using a 1:1 ratio (n = 648) (B) of patients who had been admitted with Covid-19, according to antidepressant exposure. The shaded areas represent pointwise 95% confidence intervals.

Fig. 3: Associations of SSRIs and other antidepressants with the composite endpoint of intubation or death.
figure 3

Kaplan–Meier curves of the composite endpoint of intubation or death in the full sample (A, C and E) and in the matched analytic samples using a 1:1 ratio (B, D and F) of patients who had been admitted with Covid-19, according to exposure to antidepressant classes. The shaded areas represent pointwise 95% confidence intervals.

In sensitivity analyses, multivariable Cox regression models in the full sample yielded similar results, as did univariate Cox regression models in the matched analytic samples, except for the association between SNRI antidepressants and the endpoint, which was significant in both sensitivity analyses, and that between α2-antagonist antidepressants and the outcome, which was not significant in the matched analytic sample (Table 2; Supplementary Figs. 15).

Additional analyses indicated that the association between antidepressant use and the endpoint was significant in both genders (Supplementary Table 15) and in a subsample of patients excluding those receiving an antipsychotic medication, who might be at increased risk for immune dysregulation linked to a psychotic disorder (Supplementary Table 16), supporting the robustness of our results. This association did not significantly differ across subgroups defined by baseline characteristics (Supplementary Table 17). Patients who were prescribed an antidepressant in the three months before but not during the visit were significantly at higher risk of intubation or death than those who received an antidepressant only during the visit (Supplementary Table 18). When considering diazepam as an active comparator, we found that antidepressant use was significantly associated with reduced risk of intubation or death in the crude (HR, 0.50; 95% CI, 0.32–0.77, p = 0.002) and in the main (HR, 0.47; 95% CI, 0.25–0.88, p = 0.018) analyses, as well as in sensitivity analyses (Supplementary Table 19). Participants receiving diazepam did not have a significantly different risk of intubation or death compared to those not receiving either medication. Exposure to higher rather than lower doses of antidepressants was not significantly associated with the endpoint, except for tetracyclic antidepressants, for which higher doses were associated with an increased risk (Supplementary Table 20). Exposure to a combination of antidepressants was not significantly associated with the endpoint compared to exposure to only one antidepressant (Supplementary Table 21). When using intubation and death as separate endpoints, we found that antidepressant use was significantly associated with lower risks of death and intubation (Supplementary Table 22). Finally, among the 577 patients admitted to ICUs, antidepressant use was significantly associated with reduced risk of intubation or death (Supplementary Table 23). A post-hoc analysis of statistical power is presented in Supplementary Text 2.

Discussion

In this multicenter retrospective observational study involving a large sample of patients hospitalized for COVID-19, we found that antidepressant use, at a mean dosage of 21.6 (SD = 14.1) fluoxetine-equivalent milligrams, was significantly and substantially associated with reduced risk of intubation or death, independently of patient characteristics, clinical and biological markers of disease severity, and other psychotropic medications. This association remained significant in multiple sensitivity analyses. Exploratory analyses suggest that SSRI and non-SSRI antidepressants, and specifically the SSRIs escitalopram, fluoxetine, and paroxetine, and the SNRI venlafaxine, were significantly associated with reduced risk of intubation or death. These findings should be interpreted with caution due to the observational design of the study. However, our findings provide support for conducting controlled randomized clinical trials of antidepressant medications for COVID-19.

The association observed between antidepressant use and reduced risk of intubation or death might be explained by several mechanisms. First, prior studies suggest that several antidepressants may inhibit acid sphingomyelinase activity [10], which may prevent the infection of epithelial cells with SARS-CoV-2 [11]. Second, numerous antidepressants are S1R agonists [13], which might prevent the cytokine storm observed in severe COVID-19 [3]. Third, antidepressant use may be associated with reduced plasma levels of several inflammatory mediators associated with severe COVID-19 [7], including IL-6, IL-10, TNF-α, and CCL-2 [8]. Finally, certain antidepressants may exert antiviral effects on SARS-CoV-2, such as fluoxetine [12]. If this link were confirmed, studies would be needed to determine the exact mechanisms underlying this association.

We found that individuals with psychiatric disorders and hospitalized for COVID-19 would be at higher risk of intubation or death, in line with findings of recent studies [33, 34]. Our results suggest that this association may be explained by the greater prevalence of clinical characteristics associated with severe COVID-19 in individuals with psychiatric disorders than their counterparts, including medical comorbidity, and clinical and biological markers of disease severity. Nevertheless, future studies focusing specifically on this vulnerable population are needed to examine this issue.

Our study has several limitations. First, there are two possible major inherent biases in observational studies: unmeasured confounding and confounding by indication. We tried to minimize the effects of confounding in several different ways. First, we used multivariable regression models with inverse probability weighting to minimize the effects of confounding by indication [27, 28]. Second, we performed sensitivity analyses that showed similar results. Third, although some amount of unmeasured confounding may remain, our analyses adjusted for numerous potential confounders. Furthermore, our findings indicate that the association between antidepressant use and the endpoint was only observed for patients receiving an antidepressant during the visit and not in those who received it only within the 3 months before hospital admission. Fourth, this association remained significant when using death and intubation as separate endpoints. Fifth, the low number of patients with antidepressants admitted to ICUs in our study (12 of 577 patients, i.e. 2.1%) and the significant association between antidepressant use and the endpoint in these patients further give strength to our conclusion. Finally, when considering diazepam as an active comparator, we found similar results.

Other limitations include missing data for some baseline characteristic variables, including baseline clinical and biological markers of severity of COVID-19, which may be explained by the overwhelming of all hospital units during the COVID-19 peak incidence, and potential for inaccuracies in electronic health records in this context, especially for hand-written medical prescriptions. However, most associations observed between baseline characteristics and mortality or intubation are in line with prior epidemiological studies [17, 19]. Second, our study cannot establish a causal relationship between antidepressant use and reduced risk of intubation or death [35]. Third, inflation of type I error might have occurred in exploratory analyses due to multiple testing. Fourth, information on the reason for prescribing the antidepressant, the duration and adherence to its prescription, the prescription record of patients before the admission, and the date that COVID-19 symptoms appeared, was not available. Fifth, the criteria used to determine whether a patient with COVID-19 should be hospitalized or not might have varied during the study period. However, its length was relatively limited (i.e., from January 24th until April 1st) and our main results remained unchanged following additional adjustment for it (Supplementary Tables 2426). Finally, despite the multicenter design, our results may not be generalizable to outpatients and other countries.

In this multicenter observational retrospective study, antidepressant use at usual antidepressant doses was significantly and substantially associated with lower risk of intubation or death among adult patients hospitalized for COVID-19. Double-blind controlled randomized clinical trials of antidepressant medications for COVID-19 are needed.