Observational study of haloperidol in hospitalized patients with COVID-19

About the Authors:

Nicolas Hoertel

Roles Conceptualization, Formal analysis, Methodology, Writing – original draft, Writing – review & editing

Affiliations Département de Psychiatrie, Hôpital Corentin-Celton, Centre Université de Paris, AP-HP, Issy-les-Moulineaux, France, INSERM, Institut de Psychiatrie et Neurosciences de Paris, Paris, France, Faculté de Santé, UFR de Médecine, Université de Paris, Paris, France

Marina Sánchez-Rico

Roles Conceptualization, Data curation, Formal analysis, Methodology, Writing – review & editing

* E-mail: [email protected]

Affiliations Département de Psychiatrie, Hôpital Corentin-Celton, Centre Université de Paris, AP-HP, Issy-les-Moulineaux, France, Faculty of Psychology, Department of Psychobiology & Behavioural Sciences Methods, Universidad Complutense de Madrid, Madrid, Spain

ORCID logo https://orcid.org/0000-0002-1121-8641

Raphaël Vernet

Roles Data curation, Formal analysis

Affiliation: Biostatistics and Public Health Department, Hôpital Européen Georges Pompidou, Medical Informatics, AP-HP, Centre-Université de Paris, Paris, France

Anne-Sophie Jannot

Roles Conceptualization, Data curation, Writing – review & editing

Affiliations Faculté de Santé, UFR de Médecine, Université de Paris, Paris, France, Biostatistics and Public Health Department, Hôpital Européen Georges Pompidou, Medical Informatics, AP-HP, Centre-Université de Paris, Paris, France, INSERM, UMR_S 1138, Cordeliers Research Center, Université de Paris, Paris, France

Antoine Neuraz

Roles Data curation

Affiliations INSERM, UMR_S 1138, Cordeliers Research Center, Université de Paris, Paris, France, Department of Medical Informatics, Necker-Enfants Malades Hospital, AP-HP, Centre-Université de Paris, Paris, France

Carlos Blanco

Roles Writing – review & editing

Affiliation: National Institute on Drug Abuse, Bethesda, MD, United States of America

Cédric Lemogne

Roles Writing – review & editing

Affiliations Département de Psychiatrie, Hôpital Corentin-Celton, Centre Université de Paris, AP-HP, Issy-les-Moulineaux, France, INSERM, Institut de Psychiatrie et Neurosciences de Paris, Paris, France, Faculté de Santé, UFR de Médecine, Université de Paris, Paris, France

ORCID logo https://orcid.org/0000-0002-3487-4721

Guillaume Airagnes

Roles Writing – review & editing

Affiliations Département de Psychiatrie, Hôpital Corentin-Celton, Centre Université de Paris, AP-HP, Issy-les-Moulineaux, France, INSERM, Institut de Psychiatrie et Neurosciences de Paris, Paris, France, Faculté de Santé, UFR de Médecine, Université de Paris, Paris, France

Nicolas Paris

Roles Data curation, Writing – review & editing

Affiliations AP-HP, DSI-WIND, Paris, France, LIMSI, CNRS, Université Paris-Sud, Université Paris-Saclay, Orsay, France

Christel Daniel

Roles Data curation, Writing – review & editing

Affiliations AP-HP, DSI-WIND, Paris, France, Sorbonne University, University Paris 13, Sorbonne Paris Cité, INSERM UMR_S 1142, Paris, France

ORCID logo https://orcid.org/0000-0001-7938-8156

Alexandre Gramfort

Roles Data curation, Writing – review & editing

Affiliation: Université Paris-Saclay, Inria, CEA, Palaiseau, France

Guillaume Lemaitre

Roles Data curation, Writing – review & editing

Affiliation: Université Paris-Saclay, Inria, CEA, Palaiseau, France

Mélodie Bernaux

Roles Data curation, Writing – review & editing

Affiliation: Direction de la stratégie et de la transformation, AP-HP, Paris, France

Ali Bellamine

Roles Writing – review & editing

Affiliation: Unité de Recherche clinique, Hôpital Cochin, AP-HP, Centre-Université de Paris, Paris, France

Nathanaël Beeker

Roles Conceptualization, Data curation, Writing – review & editing

Affiliation: Unité de Recherche clinique, Hôpital Cochin, AP-HP, Centre-Université de Paris, Paris, France

Frédéric Limosin

Roles Conceptualization, Writing – review & editing

Affiliations Département de Psychiatrie, Hôpital Corentin-Celton, Centre Université de Paris, AP-HP, Issy-les-Moulineaux, France, INSERM, Institut de Psychiatrie et Neurosciences de Paris, Paris, France, Faculté de Santé, UFR de Médecine, Université de Paris, Paris, France

on behalf of the AP-HP/Universities/INSERM Covid-19 research collaboration and AP-HP Covid CDR Initiative

¶Membership of the AP-HP/Universities/INSERM Covid-19 research collaboration and AP-HP Covid CDR Initiative is provided in the Acknowledgments.

Introduction

The novel coronavirus SARS-CoV-2, the causative agent of coronavirus disease 2019 (COVID-19), has caused worldwide health, social and economic disruption [1,2]. In the absence of antiviral medications with proven clinical efficacy [3,4], the search for an effective treatment for patients with COVID-19 among all available medications is urgently needed [4,5].

Based on advances in the knowledge of molecular details of SARS-CoV-2 infection [4], it has been suggested that two sets of pharmacological agents that show in-vitro antiviral activity should be prioritized in that search: the inhibitors of mRNA translation and the predicted regulators of the Sigma1 and Sigma2 receptors [4]. Molecules that target Sigma receptors may reduce virus infectivity through different mechanisms, including lipid remodeling and endoplasmic reticulum stress response [4,6].

Haloperidol, a butyrophenone-derivative antipsychotic widely used in the treatment of psychoses and delirium, has been suggested as potentially useful for patients with COVID-19 on the grounds of its in-vitro antiviral effects against SARS-CoV-2, possibly through sigma-1 receptor antagonist effect [4,7].

Short-term use of haloperidol is generally well tolerated [8], although side effects can occur, including extrapyramidal symptoms and QT interval prolongation [9].

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

To this end, 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 associations of haloperidol use with the risk of intubation or death and the time to discharge home among adult patients who have been admitted to AP-HP hospitals for COVID-19. We hypothesized that haloperidol use would be associated with lower risk of a composite endpoint of intubation or death, and with shorter time from hospital admission to discharge home in time-to-event analyses adjusting for patient characteristics and other psychotropic medications.

Methods

Setting

We conducted a multicenter observational retrospective study at AP-HP, which includes 39 hospitals, 23 of which are acute, 20 are adult and 3 are pediatric hospitals. We included all adults aged 18 years or over who have been admitted for COVID-19 to these medical centers from the beginning of the epidemic in France, i.e. January 24th, until May 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 non-interventional retrospective 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). AP-HP clinical Data Warehouse initiative ensures patients’ information and 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). Participants who did not consent to participate in the study were excluded prior to the construction of the database. 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 AP-HP hospital. The data obtained included patient demographic characteristics, 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: sex; age (binarized at the median value observed in the full sample); hospital, which was categorized into 2 classes following the administrative clustering of AP-HP hospitals in Paris and its suburbs based on their geographical location (i.e., AP-HP Centre–Paris University, Henri Mondor University Hospitals and at home hospitalization; and AP-HP Nord and Hôpitaux Universitaires Paris Seine-Saint-Denis, Paris Saclay University, and Sorbonne University); obesity, defined as having a body-mass index higher than 30 kg/m2 or an International Statistical Classification of Diseases and Related Health Problems (ICD-10) diagnosis code for obesity (E66.0, E66.1, E66.2, E66.8, E66.9); self-reported current smoking status; any medical conditions associated with increased risk of severe COVID-19 [10–14] based on ICD-10 diagnosis codes, including diabetes mellitus (E11), diseases of the circulatory system (I00-I99), diseases of the respiratory system (J00-J99), neoplasms (C00-C96), and diseases of the blood and blood-forming organs and certain disorders involving the immune mechanism (D5-D8); and any medication prescribed according to compassionate use or as part of a clinical trial (e.g. hydroxychloroquine, azithromycin, remdesivir, tocilizumab, sarilumab, or dexamethasone). To take into account possible confounding by indication bias for haloperidol, we recorded whether patients had any current current diagnosis of psychiatric disorders, including delirium (F00-F99 and R41.0), any prescribed antipsychotic other than haloperidol, and any other prescribed psychotropic medication (i.e., antidepressants, benzodiazepines, Z-drugs, and mood stabilizers, including lithium and antiepileptic medications with mood stabilizing effects).

All medical notes and prescriptions are computerized in Greater Paris University hospitals. Medications including their dose, frequency, date, and mode of administration were identified from medication administration data or scanned hand-written medical prescriptions, through two deep learning models based on BERT contextual embeddings [15], one for the medications and one for their mode of administration. The model was trained on the APmed corpus [16], a previously annotated dataset for this task. Extracted medications names were then normalized to the Anatomical Therapeutic Chemical (ATC) terminology using approximate string matching.

Haloperidol use

Study baseline was defined as the date of hospital admission. Haloperidol use was defined as receiving haloperidol within the first 48 hours 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, or the treatment may not have been recorded in the computerized medication administration data the first day of admission. In this observational study, no specific clinical guidelines were given to practitioners to administer haloperidol.

Endpoints

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. The secondary outcome was the time from study baseline to discharge home among survivors. Patients without an end-point event had their data censored on May 1st, 2020.

Statistical analysis

We calculated frequencies and means (± standard deviations (SD)) of each variable described above in patients receiving or not receiving haloperidol and compared them using standardized mean differences (SMD). A SMD higher than 0.1 was considered to reflect substantial imbalance [17].

To examine the association of haloperidol use with the primary composite endpoint of intubation or death and the secondary endpoint of discharge home among survivors, we performed Cox proportional-hazards regression models [18]. To help account for the nonrandomized prescription of haloperidol and reduce the effects of confounding, the primary analysis used a univariate Cox regression model in a matched analytic sample for each outcome. We selected four controls for each exposed case, based on age, sex, hospital, obesity, smoking status, any medical condition, any medication prescribed according to compassionate use or as part of a clinical trial, any current diagnosis of psychiatric disorders, any prescribed antipsychotic other than haloperidol, and any other prescribed psychotropic medication [19–22]. To reduce the effects of confounding, optimal matching was used in order to obtain the smallest average absolute distance across all these characteristics between each exposed patient and its corresponding non-exposed matched controls [23]. Weighted Cox regression models were used when proportional hazards assumption was not met [24]. Kaplan-Meier curves were performed [25] and their 95% pointwise confidence intervals were estimated using the nonparametric bootstrap method [26].

We conducted six sensitivity analyses. First, we performed propensity score analyses with inverse probability weighting (IPW) [27,28]. The individual propensities for haloperidol prescription were estimated by a multivariable logistic regression model that included as covariates the same variables used in the primary analyses. The predicted probabilities from the propensity-score model were used to calculate the stabilized inverse-probability-weighting weights [27]. Associations between haloperidol use and the two outcomes were then estimated using IPW Cox regression models. In cases of non-balanced covariates, IPW multivariable Cox regression models adjusting for these non-balanced covariates were also performed. Kaplan-Meier curves were performed using the inverse-probability-weighting weights [25,26]. Second, we performed multivariable Cox regression models including as covariates the same variables as in the primary analyses and the inverse-probability-weighted analyses. Third, to address a potential immortality bias in the exposed group due to a treatment initiation after hospital admission, we performed multivariable cox regression models while considering haloperidol use as a time dependent variable [18], including all participants who received haloperidol at any time from hospital admission until the end of the index hospitalization or intubation or death. In this type of analysis, patients who received haloperidol after study baseline were allowed to come into the analysis risk-sets at the time of actual first initiation of haloperidol. Fourth, we examined whether our findings were similar in models imputing missing data using multiple imputation [29] instead of excluding patients with any missing data as done in the main analyses. Fifth, in order to account for potential latent effects of the variable hospital, we examined whether our findings were similar while considering this variable as a random effects covariate in the main analyses [18]. Finally, because psychotropic medications other than haloperidol could have been prescribed to patients after they received haloperidol, we examined whether the results were similar when not including these variables as covariates in the main analyses.

We also performed additional analyses. First, to increase our confidence that the results might not be due to unmeasured confounding or indication bias, we examined whether the two endpoints differed between patients receiving haloperidol only in the 3 months before hospital admission and those who received it only during the visit. Second, we examined a potential dose-effect relationship by testing the association between the daily dose received (dichotomized at the median value) with the two endpoints among patients who received haloperidol.

For all associations, we performed residual analyses to assess the fit of the data, check assumptions, including proportional hazards assumption using proportional hazards tests and diagnostics based on weighted residuals [30], and examined the potential presence and influence of outliers [31]. Statistical significance was fixed a priori at p<0.05. All analyses were conducted in R software version 2.4.3. Statistical code used for the main analyses has been deposited in a recognized public source code repository (GitHub, https://github.com/mlsrico/haloperidol_and_covid19).

Results

Characteristics of the cohort

Of the 17,076 hospitalized adult patients with a positive COVID-19 RT-PCR test, 1,908 patients (11.2%) were excluded because of missing data or their young age (i.e. less than 18 years of age). Of the 86 adult patients who received haloperidol at any time during the visit, 47 (54.7%) patients were excluded because they received it more than 48 hours from hospital admission (N = 47, 54.7%) or after intubation (N = 16, 18.6%). Of the remaining 15,121 inpatients, 39 patients (0.3%) received haloperidol at baseline (i.e., within the first 48 hours of hospital admission) at a median daily dose of 3.0 mg per day (SD = 5.2; mean = 4.5; first quartile = 1.0; third quartile = 5.0; range = 0.5–20.0 mg) for a median duration of 7.0 days (SD = 7.2; mean = 8.4; first quartile = 2.5; third quartile = 12.0; range = 1–26). Of these 39 patients, 4 patients (10.2%) had a medication administration by intramuscular injection. Median delay between study baseline and haloperidol initiation was lower than 0 day (SD = 1.0; mean = 0.6; first quartile = 0.0; third quartile = 1.0; range = 0–2 days) (S1 Fig).

First positive COVID-19 RT-PCR tests were obtained after a median delay of 1.2 days (SD = 12.7) from study baseline. This delay was significantly but not substantially different between patients receiving or not receiving haloperidol [median in the exposed group = 1.0 day (SD = 11.2); median in the non-exposed group = 1.2 days (SD = 12.8); Mood’s median test Chi-square = 3.76, p = 0.001].

Over a median follow-up of 7 days (SD = 17.9; mean = 13.8; first quartile = 1.0; third quartile = 22.0; range: 1–98 days), 2,024 patients (13.4%) had a primary end-point event and 10,179 patients (77.6%) were discharged home at the time of study end on May 1st. Patients receiving haloperidol had a median follow-up of 7 days (SD = 14.0; mean = 10.9; first quartile = 4; third quartile = 22; range: 1–80 days), while the non-exposed group had a median follow-up of 7 days (SD = 17.9; mean = 13.8; first quartile = 1; third quartile = 22; range: 1–98 days) [Mood’s median test Chi-square = 0.41, p = 0.679].

All baseline characteristics, when examined independently, were significantly associated with both endpoints (S1 and S2 Tables). In the multivariable analysis, these associations remained significant for the primary outcome, except for hospital type, smoking status, any medication according to compassionate use or as part of a clinical trial, and any antipsychotic other than haloperidol, as they did for the secondary outcome, except for smoking status.

The distribution of the patient characteristics according to haloperidol use is shown in Table 1. In the full sample, haloperidol use substantially differed (i.e., SMD>0.1) across all characteristics except for any medication according to compassionate use or as part of a clinical trial. After applying the propensity score weights, all these differences became not substantial (i.e., SMD<0.1), except for any antipsychotic other than haloperidol and any other psychotropic medication (Table 1). In the matched analytic sample comprising 195 patients (i.e., 39 patients who received haloperidol at baseline and 156 patients who did not receive haloperidol during the visit from the matched group), there were no substantial differences in patient characteristics according to haloperidol use (all SMD<0.1) (Table 1).

[Figure omitted. See PDF.]

Table 1. Characteristics of hospitalized patients with COVID-19 receiving or not receiving haloperidol in the full sample and in the matched analytic sample.

https://doi.org/10.1371/journal.pone.0247122.t001

Study endpoints

The primary endpoint event of intubation or death respectively occurred in 9 patients (23.1%) who received haloperidol and 2,015 patients (13.4%) who did not (Table 2). In both the crude, unadjusted analysis and the primary analysis using a univariate Cox regression model in the matched analytic sample, there were no significant associations between haloperidol use and the primary endpoint (hazard ratio (HR), 1.68; 95% CI, 0.87 to 3.23; p = 0.120; and HR, 0.80; 95% CI, 0.39 to 1.62; p = 0.531, respectively) (Fig 1; Table 2).

[Figure omitted. See PDF.]

Fig 1.

Kaplan-Meier curves for intubation or death in the full sample crude analysis (N = 15,121) (A), in the full sample analysis with inverse probability weighting (N = 15,121) (B) and in the matched analytic sample using a 1:4 ratio (N = 195) (C) of patients who had been hospitalized for COVID-19, according to haloperidol use. The shaded areas represent pointwise 95% confidence intervals.

https://doi.org/10.1371/journal.pone.0247122.g001

[Figure omitted. See PDF.]

Table 2. Associations between haloperidol use and the endpoints of intubation or death and discharge home among survivors, in the full sample and in the matched analytic sample of patients hospitalized for COVID-19.

https://doi.org/10.1371/journal.pone.0247122.t002

Among survivors, the secondary endpoint of discharge home occurred in 16 patients (61.5%) who were prescribed haloperidol and 9,907 patients (85.8%) who were not. Haloperidol use was significantly and negatively associated with the secondary endpoint in the crude, unadjusted analysis (HR, 0.24; 95% CI, 0.13 to 0.44; p<0.001), but this association was not significant in the primary analysis using a univariate Cox regression model in the matched analytic sample (HR, 1.30; 95% CI, 0.74 to 2.28; p = 0.355) (Fig 2; Table 2).

[Figure omitted. See PDF.]

Fig 2.

Kaplan-Meier curves for discharge home among survivors in the full sample crude analysis (N = 11,572) (A), in the full sample analysis with inverse probability weighting (N = 11,572) (B) and in the matched analytic sample using a 1:4 ratio (N = 130) (C) of patients who had been hospitalized for COVID-19, according to haloperidol use. The shaded areas represent pointwise 95% confidence intervals.

https://doi.org/10.1371/journal.pone.0247122.g002

Sensitivity analyses, including multivariable Cox regression models and propensity score analyses with inverse probability weighting in the full sample yielded similar non-significant results for the two endpoints (Figs 1 and 2; Table 2). Similar non-significant results were found in the inverse probability weighting analyses adjusting for the two unbalanced covariates (i.e. any antipsychotic medication other than haloperidol and any other psychotropic medication) (i.e., HR, 1.32; 95% CI, 0.65 to 2.71; p = 0.441 for the primary outcome; and HR, 1.18; 95% CI, 0.71 to 1.95; p = 0.528 for the secondary outcome), as well as in the analyses including all participants who received haloperidol during the hospitalization until the end of the index hospitalization or intubation or death and considering haloperidol use as a time dependent variable (S3 Table). Findings were also similar when considering hospital as a random effects variable (S4 Table). Finally, models imputing missing data using multiple imputation yielded very similar results as in the main analyses (S5 Table), as did models not including psychotropic medications other than haloperidol as covariates (S6 Table).

Additional analyses indicated that risks for both endpoints were not significantly different between patients who were prescribed haloperidol only in the three months before the hospitalization and those who received it only during the visit (S7 Table). Exposure to higher rather than lower doses of haloperidol was not significantly associated with the primary or secondary endpoints (S8 and S9 Tables).

A post-hoc analysis indicated that in the full sample, we had 80% power to detect unweighted and unadjusted hazard ratios of at least 0.20/2.41 for the primary endpoint and 0.27/5.57 for the secondary endpoint, while we had 80% power to detect unweighted and unadjusted hazard ratios of at least 0.39/2.18 for the primary endpoint and 0.39/3.07 for the secondary endpoint in the matched analytic sample.

Discussion

In this multicenter retrospective observational study involving a large number of adult patients hospitalized for COVID-19, the risk of intubation or death and the time to discharge home among survivors were not significantly different between patients who received haloperidol and those who did not. Although these findings should be interpreted with caution due to the observational design, the relatively wide confidence intervals for estimates, and the fact that this is, to our knowledge, the first study examining these associations in patients with COVID-19, they suggest that haloperidol prescribed at a mean daily dose of 4.5 mg per day (SD = 5.2) for a mean duration of 8.4 days (SD = 7.2) is not associated with risk of intubation or death, or with time to discharge home, among patients hospitalized for COVID-19.

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 a univariate Cox regression models in matched analytic samples and multivariable Cox regression models with inverse probability weighting to minimize the effects of confounding by indication [27,28]. Second, although some amount of unmeasured confounding may remain, our analyses adjusted for numerous potential confounders, including age, sex, hospital, obesity, current smoking status, any medical condition, any medication prescribed according to compassionate use or as part of a clinical trial, any current diagnosis of psychiatric disorders, any prescribed antipsychotic other than haloperidol, and any other psychotropic medication. Finally, the lack of significant associations between the daily dose of haloperidol and the two endpoints also supports our conclusion.

Other limitations include missing data for some variables (i.e., 11.2%) and potential for inaccuracies in the electronic health records, such as the possible lack of documentation of illnesses or medications, or the misidentification of treatment mode of administration (e.g., dose), especially for hand-written medical prescriptions. However, results remained unchanged when using multiple imputation to account for missing data. Second, given the limited number of patients who received haloperidol, our analyses were powered to detect only substantial effect sizes (i.e., 80% power to detect unweighted and unadjusted hazard ratios of at least 0.20/2.41 for the primary endpoint and 0.27/5.57 for the secondary endpoint in the full sample). In addition, the fact that the inverse probability weighting analyses in the full sample were not successful in balancing two covariates between the two groups (i.e. exposed to haloperidol vs. non-exposed) likely has led to reduced statistical power. Third, haloperidol was prescribed at a relatively low dose, i.e., at a mean daily dose of 4.5 mg per day (SD = 5.2), and its antiviral properties might be observable at higher doses. Fourth, associations reported in our study do not imply causal relationships [32]. Finally, despite the multicenter design, our results may not be generalizable to other settings, e.g. outpatients, or regions.

Conclusion

Findings form this multicenter observational study suggest that haloperidol prescribed at a mean daily dose of 4.5 mg per day (SD = 5.2) for a mean duration of 8.4 days (SD = 7.2) may not be associated with risk of intubation or death, or with time to discharge home, among patients hospitalized for COVID-19.

Supporting information

S1 Fig. Study cohort.

https://doi.org/10.1371/journal.pone.0247122.s001

(DOCX)

S1 Table. Associations of baseline clinical characteristics with the composite endpoint of intubation or death in the cohort of adult patients hospitalized for COVID-19 (N = 15,121).

https://doi.org/10.1371/journal.pone.0247122.s002

(DOCX)

S2 Table. Associations of baseline clinical characteristics with the endpoint of discharge home in the cohort of adult patients hospitalized for COVID-19 who survived (N = 11,572).

https://doi.org/10.1371/journal.pone.0247122.s003

(DOCX)

S3 Table. Associations between haloperidol use and the endpoints of intubation or death and discharge home among survivors, in the full sample and in the matched analytic sample of patients hospitalized for COVID-19, while including all patients who received haloperidol at any time during the visit and before the end of the index hospitalization or intubation or death and considering haloperidol use as a time dependent variable.

https://doi.org/10.1371/journal.pone.0247122.s004

(DOCX)

S4 Table. Associations between haloperidol use and the endpoints of intubation or death and discharge home among survivors, in the full sample and in the matched analytic sample of patients hospitalized for COVID-19, when considering the variable hospital as a random effects variable.

https://doi.org/10.1371/journal.pone.0247122.s005

(DOCX)

S5 Table. Associations between haloperidol use and the endpoints of intubation or death and discharge home among survivors, in the full sample and in the matched analytic sample of patients hospitalized for COVID-19, following imputation of missing data with multiple imputation.

https://doi.org/10.1371/journal.pone.0247122.s006

(DOCX)

S6 Table. Associations between haloperidol use and the endpoints of intubation or death and discharge home among survivors, in the full sample and in the matched analytic sample of patients hospitalized for COVID-19, in multivariable models not including psychotropic medications other than haloperidol as covariates.

https://doi.org/10.1371/journal.pone.0247122.s007

(DOCX)

S7 Table. Comparing haloperidol use only during the three months before hospitalization versus only during the visit for the endpoints of intubation or death and discharge home among survivors.

https://doi.org/10.1371/journal.pone.0247122.s008

(DOCX)

S8 Table. Association between haloperidol dose and the endpoint of intubation or death.

https://doi.org/10.1371/journal.pone.0247122.s009

(DOCX)

S9 Table. Association between haloperidol dose and the endpoint of discharge home among survivors.

https://doi.org/10.1371/journal.pone.0247122.s010

(DOCX)

Acknowledgments

The authors warmly thank the EDS APHP Covid consortium integrating the APHP Health Data Warehouse team as well as all the APHP staff and volunteers who contributed to the implementation of the EDS-Covid database and operating solutions for this database.

Collaborators EDS APHP Covid consortium: Pierre-Yves ANCEL, Alain BAUCHET, Nathanaël BEEKER, Vincent BENOIT, Mélodie BERNAUX, Ali BELLAMINE, Romain BEY, Aurélie BOURMAUD, Stéphane BREANT, Anita BURGUN, Fabrice CARRAT, Charlotte CAUCHETEUX, Julien CHAMP, Sylvie CORMONT, Christel DANIEL, Julien DUBIEL, Catherine DUCLOAS, Loic ESTEVE, Marie FRANK, Nicolas GARCELON, Alexandre GRAMFORT, Nicolas GRIFFON, Olivier GRISEL, Martin GUILBAUD, Claire HASSEN-KHODJA, François HEMERY, Martin HILKA, Anne Sophie JANNOT, Jerome LAMBERT, Richard LAYESE, Judith LEBLANC, Léo LEBOUTER, Guillaume LEMAITRE, Damien LEPROVOST, Ivan LERNER, Kankoe LEVI SALLAH, Aurélien MAIRE, Marie-France MAMZER, Patricia MARTEL, Arthur MENSCH, Thomas MOREAU, Antoine NEURAZ, Nina ORLOVA, Nicolas PARIS, Bastien RANCE, Hélène RAVERA, Antoine ROZES, Elisa SALAMANCA, Arnaud SANDRIN, Patricia SERRE, Xavier TANNIER, Jean-Marc TRELUYER, Damien VAN GYSEL, Gaël VAROQUAUX, Jill Jen VIE, Maxime WACK, Perceval WAJSBURT, Demian WASSERMANN, Eric ZAPLETAL. Membership of the AP-HP/Universities/INSERM Covid-19 research collaboration and AP-HP Covid CDR Initiative is provided at https://www.aphp.fr/, https://u-paris.fr/, and https://www.inserm.fr/.

Disclaimer

The views and opinions expressed in this report are those of the authors and should not be construed to represent the views of any of the sponsoring organizations, agencies, or the US government.

Citation: Hoertel N, Sánchez-Rico M, Vernet R, Jannot A-S, Neuraz A, Blanco C, et al. (2021) Observational study of haloperidol in hospitalized patients with COVID-19. PLoS ONE 16(2): e0247122. https://doi.org/10.1371/journal.pone.0247122

References

1. Hoertel N, Blachier M, Blanco C, Olfson M, Massetti M, Sanchez-Rico M, et al. A stochastic agent-based model of the SARS-CoV-2 epidemic in France. Nature Medicine. 2020. pp. 1417–1421. pmid:32665655

2. Hoertel N, Blachier M, Blanco C, Olfson M, Massetti M, Limosin F, et al. Facing the COVID-19 epidemic in NYC: a stochastic agent-based model of various intervention strategies. medRxiv. 2020. pmid:32511467

3. Anderson RM, Heesterbeek H, Klinkenberg D, Hollingsworth TD. How will country-based mitigation measures influence the course of the COVID-19 epidemic? Lancet. 2020;395: 931–934. pmid:32164834

4. Gordon DE, Jang GM, Bouhaddou M, Xu J, Obernier K, White KM, et al. A SARS-CoV-2 protein interaction map reveals targets for drug repurposing. Nature. 2020;583: 459–468. pmid:32353859

5. Chevance A, Gourion D, Hoertel N, Llorca P-M, Thomas P, Bocher R, et al. Ensuring mental health care during the SARS-CoV-2 epidemic in France: A narrative review. L’Encephale. 2020;46: 193–201. pmid:32370982

6. Mitsuda T, Omi T, Tanimukai H, Sakagami Y, Tagami S, Okochi M, et al. Sigma-1Rs are upregulated via PERK/eIF2α/ATF4 pathway and execute protective function in ER stress. Biochem Biophys Res Commun. 2011;415: 519–525. pmid:22079628

7. Dalwadi DA, Kim S, Schetz JA. Activation of the sigma-1 receptor by haloperidol metabolites facilitates brain-derived neurotrophic factor secretion from human astroglia. Neurochem Int. 2017;105: 21–31. pmid:28188803

8. Zayed Y, Barbarawi M, Kheiri B, Banifadel M, Haykal T, Chahine A, et al. Haloperidol for the management of delirium in adult intensive care unit patients: a systematic review and meta-analysis of randomized controlled trials. J Crit Care. 2019;50: 280–286. pmid:30665181

9. Xiong GL, Pinkhasov A, Mangal J, Huang H, Rado J, Gagliardi J, et al. QTc monitoring in adults with medical and psychiatric comorbidities: Expert consensus from the Association of Medicine and Psychiatry. J Psychosom Res. 2020; 110138. pmid:32442893

10. Zhou F, Yu T, Du R, Fan G, Liu Y, Liu Z, et al. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study. The lancet. 2020. pmid:32171076

11. Ruan Q, Yang K, Wang W, Jiang L, Song J. Clinical predictors of mortality due to COVID-19 based on an analysis of data of 150 patients from Wuhan, China. Intensive Care Med. 2020;46: 846–848. pmid:32125452

12. Salje H, Kiem CT, Lefrancq N, Courtejoie N, Bosetti P, Paireau J, et al. Estimating the burden of SARS-CoV-2 in France. Science. 2020;369: 208–211. pmid:32404476

13. Hur K, Price CP, Gray EL, Gulati RK, Maksimoski M, Racette SD, et al. Factors Associated With Intubation and Prolonged Intubation in Hospitalized Patients With COVID-19. Otolaryngol Neck Surg. 2020. pmid:32423368

14. Williamson E, Walker AJ, Bhaskaran KJ, Bacon S, Bates C, Morton CE, et al. OpenSAFELY: factors associated with COVID-19-related hospital death in the linked electronic health records of 17 million adult NHS patients. MedRxiv. 2020. https://doi.org/10.1101/2020.05.06.20092999.

15. Devlin J, Chang M-W, Lee K, Toutanova K. Bert: Pre-training of deep bidirectional transformers for language understanding. arXiv. 2018. 1810.04805.

16. Jouffroy J, Feldman SF, Lerner I, Rance B, Neuraz A, Burgun A. MedExt: combining expert knowledge and deep learning for medication extraction from French clinical texts. 2020.

17. Stuart EA, Lee BK, Leacy FP. Prognostic score–based balance measures can be a useful diagnostic for propensity score methods in comparative effectiveness research. J Clin Epidemiol. 2013;66: S84–S90. pmid:23849158

18. Therneau Terry M., Grambsch Patricia M. Modeling Survival Data: Extending the Cox Model. New York: Springer; 2000.

19. Hoertel N, Sánchez M, Vernet R, Beeker N, Neuraz A, Alvarado J, et al. Dexamethasone use and Mortality in Hospitalized Patients with Coronavirus Disease 2019: a Multicenter Retrospective Observational Study. medRxiv. 2020. https://doi.org/10.1101/2020.10.23.20218172.

20. Hoertel N, Sanchez Rico M, Vernet R, Jannot A-S, Neuraz A, Blanco C, et al. Observational Study of Chlorpromazine in Hospitalized Patients with Covid-19. medRxiv. 2020.

21. Hoertel N, Sánchez M, Vernet R, Beeker N, Neuraz A, Blanco C, et al. Association between Hydroxyzine Use and Reduced Mortality in Patients Hospitalized for Coronavirus Disease 2019: Results from a multicenter observational study. medRxiv. 2020.

22. Hoertel N, Sánchez-Rico M, Vernet R, Beeker N, Jannot A-S, Neuraz A, et al. Association between antidepressant use and reduced risk of intubation or death in hospitalized patients with COVID-19: results from an observational study. Mol Psychiatry. 2021. pmid:33536545

23. Hansen BB, Klopfer SO. Optimal full matching and related designs via network flows. J Comput Graph Stat. 2006;15: 609–627.

24. Dunkler D, Ploner M, Schemper M, Heinze G. Weighted Cox Regression Using the R Package coxphw. J Stat Softw. 2018;84: 1–26. pmid:30450020

25. Kassambara A, Kosinski M, Biecek P. survminer: Drawing Survival Curves using “ggplot2.” 2020. Available: https://CRAN.R-project.org/package=survminer.

26. Efron B. Nonparametric standard errors and confidence intervals. Can J Stat. 1981;9: 139–158.

27. Robins JM, Hernán MÁ, Brumback B. Marginal Structural Models and Causal Inference in Epidemiology. Epidemiology. 2000;11: 550–560. pmid:10955408

28. Geleris J, Sun Y, Platt J, Zucker J, Baldwin M, Hripcsak G, et al. Observational Study of Hydroxychloroquine in Hospitalized Patients with Covid-19. N Engl J Med. 2020;382: 2411–2418. pmid:32379955

29. Stekhoven DJ, Buehlmann P. MissForest—non-parametric missing value imputation for mixed-type data. Bioinformatics. 2012;28: 112–118. pmid:22039212

30. Grambsch PM, Therneau TM. Proportional hazards tests and diagnostics based on weighted residuals. Biometrika. 1994;81: 515–526.

31. Von Elm E, Altman DG, Egger M, Pocock SJ, Gøtzsche PC, Vandenbroucke JP. The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement: guidelines for reporting observational studies. Ann Intern Med. 2007;147: 573–577. pmid:17938396

32. Le Strat Y, Hoertel N. Correlation is no causation: gymnasium proliferation and the risk of obesity. Addict Abingdon Engl. 2011;106: 1871–1872. pmid:21917042