key: cord-0841040-01qy7hc4
authors: Lund, Lars Christian; Kristensen, Kasper Bruun; Reilev, Mette; Christensen, Steffen; Thomsen, Reimar Wernich; Christiansen, Christian Fynbo; Støvring, Henrik; Johansen, Nanna Borup; Brun, Nikolai Constantin; Hallas, Jesper; Pottegård, Anton
title: Adverse outcomes and mortality in users of non-steroidal anti-inflammatory drugs who tested positive for SARS-CoV-2: A Danish nationwide cohort study
date: 2020-09-08
journal: PLoS Med
DOI: 10.1371/journal.pmed.1003308
sha: e81d87ea35be397b6cafe6853d9dff21b5b7a1d1
doc_id: 841040
cord_uid: 01qy7hc4

BACKGROUND: Concerns over the safety of non-steroidal anti-inflammatory drug (NSAID) use during severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection have been raised. We studied whether use of NSAIDs was associated with adverse outcomes and mortality during SARS-CoV-2 infection. METHODS AND FINDINGS: We conducted a population-based cohort study using Danish administrative and health registries. We included individuals who tested positive for SARS-CoV-2 during the period 27 February 2020 to 29 April 2020. NSAID users (defined as individuals having filled a prescription for NSAIDs up to 30 days before the SARS-CoV-2 test) were matched to up to 4 non-users on calendar week of the test date and propensity scores based on age, sex, relevant comorbidities, and use of selected prescription drugs. The main outcome was 30-day mortality, and NSAID users were compared to non-users using risk ratios (RRs) and risk differences (RDs). Secondary outcomes included hospitalization, intensive care unit (ICU) admission, mechanical ventilation, and acute renal replacement therapy. A total of 9,236 SARS-CoV-2 PCR-positive individuals were eligible for inclusion. The median age in the study cohort was 50 years, and 58% were female. Of these, 248 (2.7%) had filled a prescription for NSAIDs, and 535 (5.8%) died within 30 days. In the matched analyses, treatment with NSAIDs was not associated with 30-day mortality (RR 1.02, 95% CI 0.57 to 1.82, p = 0.95; RD 0.1%, 95% CI −3.5% to 3.7%, p = 0.95), risk of hospitalization (RR 1.16, 95% CI 0.87 to 1.53, p = 0.31; RD 3.3%, 95% CI −3.4% to 10%, p = 0.33), ICU admission (RR 1.04, 95% CI 0.54 to 2.02, p = 0.90; RD 0.2%, 95% CI −3.0% to 3.4%, p = 0.90), mechanical ventilation (RR 1.14, 95% CI 0.56 to 2.30, p = 0.72; RD 0.5%, 95% CI −2.5% to 3.6%, p = 0.73), or renal replacement therapy (RR 0.86, 95% CI 0.24 to 3.09, p = 0.81; RD −0.2%, 95% CI −2.0% to 1.6%, p = 0.81). The main limitations of the study are possible exposure misclassification, as not all individuals who fill an NSAID prescription use the drug continuously, and possible residual confounding by indication, as NSAIDs may generally be prescribed to healthier individuals due to their side effects, but on the other hand may also be prescribed for early symptoms of severe COVID-19. CONCLUSIONS: Use of NSAIDs was not associated with 30-day mortality, hospitalization, ICU admission, mechanical ventilation, or renal replacement therapy in Danish individuals who tested positive for SARS-CoV-2. TRIAL REGISTRATION: The European Union electronic Register of Post-Authorisation Studies EUPAS34734

a1111111111 a1111111111 a1111111111 a1111111111 a1111111111

As severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) began to transmit and spread in Europe and the United States, a letter suggesting that ibuprofen could influence the prognosis of coronavirus disease 2019 (COVID-19) through upregulation of angiotensin converting enzyme 2 receptors was circulated widely [1] . Because of this letter and case reports of otherwise healthy young patients with severe COVID-19 who had used NSAIDs in the early stage of disease [2] , concerns regarding the safety of NSAID use during the COVID-19 pandemic were widely circulated, including warnings against NSAID use in COVID-19 from the French health minister [2] and the World Health Organization. However, the European Medicines Agency stated that evidence to support warnings against use of NSAIDs in COVID-19 was lacking and stressed the need for further evidence on any effects of NSAIDs on disease prognosis in COVID-19 [3] .

The available evidence stems mainly from studies on community-acquired bacterial pneumonia and shows that use of NSAIDs is associated with bacterial complications, specifically empyema and lung abscesses [4] [5] [6] [7] . For viral illness, use of NSAIDs was not associated with mortality in intensive care unit (ICU) patients with influenza H1N1 during the 2009 pandemic [8] , and a recent study found that use of NSAIDs was not associated with mortality in patients hospitalized for influenza [9] . As use of ibuprofen and other NSAIDs is widespread, data on their safety are urgently needed to guide clinicians and patients on how to use NSAIDs during the COVID-19 pandemic. We therefore examined whether use of NSAIDs prior to infection with SARS-CoV-2 was associated with increased risk of hospitalization, ICU admission, and mortality compared to non-use of NSAIDs.

We conducted a Danish nationwide registry-based cohort study investigating the association between NSAID use and 30-day mortality and other adverse outcomes, specified as hospitalization, ICU admission, mechanical ventilation, and acute renal replacement therapy, in all patients who tested positive for SARS-CoV-2. For a graphical representation of the study design, see Fig 1. All analyses followed the publicly registered protocol [10] , except for a change in the matching algorithm and a post hoc analysis of test-negative individuals (both detailed below). This study is reported as per the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guideline (S1 STROBE Checklist). The institutional data protection board at the University of Southern Denmark and the Danish Health Data Authority approved the research project. According to Danish law, studies based entirely on registry data do not require approval from an ethics review board [11] . All source code used to conduct this study is freely available at https://source.coderefinery.org/lcl/nsaid-covid19.

Data on all Danish residents with a positive test for SARS-CoV-2 were obtained from Danish health and administrative registries as described elsewhere [12] . In brief, identification of the study population was based on prospectively collected data on all Danish residents receiving a polymerase chain reaction (PCR) test for SARS-CoV-2 from the Danish Microbiology Database [13] . Data were linked to the Danish Civil Registration System [14] , the Danish National Prescription Registry [15] , the Danish National Patient Registry [16] , and the Danish Register of Causes of Death [17] by means of the unique personal identifier assigned to all Danish residents at birth or immigration. Data on ICU treatment, mechanical ventilation, and renal replacement therapy from the Danish National Patient Registry were supplemented with daily reports on admitted patients from the 5 Danish regions [16, 18] . For more details regarding individual registries, see S1 Appendix.

All Danish residents who had a positive PCR test for SARS-CoV-2 during the period 27 February 2020 to 29 April 2020 were included in the study. To ensure complete information on exposure and covariates prior to cohort entry, individuals with less than 1 year of residence in Denmark prior to the positive test for SARS-CoV-2 were excluded. For all secondary outcomes, individuals with an outcome during 30 days to 1 day prior to cohort entry were excluded, partly to ensure that outcomes were incident and plausibly occurring due to 

COVID-19 and partly because in-hospital drug use was not available from the Danish National Prescription Registry.

We conducted a post hoc supplementary analysis where we examined the same association in a cohort of all Danish patients who tested negative for SARS-CoV-2 in the study period.

The exposure of interest was current use of any NSAID prior to a positive SARS-CoV-2 test. Current use was defined as having filled a prescription for any NSAID in the 30 days prior to the positive test. Filled NSAID prescriptions were identified from the Danish National Prescription Registry, with information on all dispensed prescriptions at community pharmacies in Denmark since 1995 [19] . Users of NSAIDs were compared to individuals without NSAID use in the corresponding time window. In Denmark, NSAIDs are sold by prescription except for low-dose (200 mg) ibuprofen sold over-the-counter in pack sizes of no more than 20 tablets. In 2018, over-the-counter purchases of ibuprofen constituted 15% of total ibuprofen sales and a smaller proportion of total NSAID sales [20] . Hence, the potential to identify NSAID use from the Danish National Prescription Registry is high compared to many other countries where over-the-counter use of NSAIDs is common [19, 21] .

The main outcome of interest was death within 30 days of a positive test for SARS-CoV-2.

The secondary outcomes included hospitalization, ICU admission, mechanical ventilation, and acute renal replacement therapy within 14 days of a positive test for SARS-CoV-2.

Eligible individuals were followed until the end of follow up (30 days for the main outcome, 14 days for secondary outcomes) or the event of interest.

We used propensity score (PS) matching to increase comparability between NSAID users and non-users. The PS is the estimated probability of receiving the treatment of interest given a set of characteristics [22] . PSs were estimated on the day of the positive SARS-CoV-2 test using logistic regression. Independent variables in the PS model were age, sex, relevant comorbidities, use of selected prescription drugs, and phase of the outbreak. For details of these independent variables, see S2 Appendix. A separate PS was estimated for each exposure definition in the main and supplementary analyses. To evaluate the appropriateness of the model, PS distributions were plotted separately for each cohort and overlap assessed visually. To reduce unmeasured confounding, individuals in the tails of the PS distribution were trimmed asymmetrically [23] . Up to 4 non-users were matched to each NSAID user using a nearest neighbor algorithm. Non-users could be matched to multiple NSAID users, and the maximum allowed difference in the PS between matches was 0.05 [24] . The Danish SARS-CoV-2 test strategy was subject to marked changes from a limited capacity setting in the beginning of the study period to a setting where widespread testing was available at the end of the study period. To account for this [25] , we included calendar week of the test date as a forced matching variable. This decision was made post hoc, i.e., not included in the protocol. Covariate balance before and after matching was assessed using standardized mean differences [26] .

Descriptive statistics were used to describe NSAID users and non-users at baseline. Continuous variables were reported as medians and interquartile ranges. Dichotomous variables were reported as frequencies and percentages.

Risks and risk differences (RDs) were estimated using generalized linear models with a binomial distribution and an identity link. Risk ratios (RRs) were estimated similarly but using a log link. Matched analyses were implemented using frequency weighting, i.e., NSAID users were assigned a weight of 1, and non-users' weights were assigned according to each individual user's number of matches. Robust 95% confidence intervals were calculated using the sandwich estimator of variance where the assumption regarding independence of observations was relaxed in the matched analyses. Data management and statistical analyses were performed using Stata 16 MP. The codes used to define exposures, covariates, and outcomes are available in S2 Appendix.

To explore treatment effect heterogeneity, we repeated the main analyses stratifying by age (<65 years, 65+ years), sex, and history of cardiovascular disease. To examine whether widespread testing of healthcare workers influenced the findings, we repeated the main analyses excluding healthcare workers from the study population. We used the same PS as estimated in the main analyses for the subgroup analyses [27] .

We conducted the following supplementary analyses. (1) We relaxed the exposure definition by using an extended NSAID exposure assessment window of 60 days prior to the positive test and repeated the main analyses with this exposure definition. (2) To explore whether reverse causation may have influenced the findings, effect estimates were obtained using an exposure assessment window of 60 days to 14 days before cohort entry (i.e., disregarding NSAID prescriptions filled during the 14 days immediately prior to cohort entry). (3) To evaluate the robustness of the findings with regards to the chosen outcome assessment windows, we obtained 60-day risk estimates for mortality and 30-day risk estimates for secondary outcomes. (4) To examine the potential for residual confounding, we conducted a negative control analysis by repeating the main analyses within the test-negative population, i.e. individuals who tested negative for SARS-CoV-2 (and did not later test positive). If an individual was tested more than once, the first test date was used as the cohort entry date. This post hoc analysis was not specified in the protocol.

We identified 9,370 individuals who tested positive for SARS-CoV-2 during the study period. Of these, 134 were excluded due to migration within 1 year prior to cohort entry, resulting in an eligible population of 9,236 individuals followed for a total of 705 person-years. The median age in the study cohort was 50 years, and 58% were female. Overall, 535 individuals (5.8%) died within 30 days, 1,512 (16%) were hospitalized within 14 days, 290 (3.1%) were admitted to the ICU, 235 (2.5%) received mechanical ventilation, and 61 (0.7%) received acute renal replacement therapy. In total, 248 (2.7%) patients had filled a prescription for an NSAID within 30 days before the test date. Compared to non-users, NSAID users were older (median age 55 versus 49 years) and more likely to have a history of hospital-diagnosed overweight or obesity (13% versus 9%), to have medical indications for NSAID use such as osteoarthritis (19% versus 12%) or rheumatoid arthritis (7% versus 3%), and to have been prescribed opioids the year before sampling date (24% versus 11%). After matching, covariates were well balanced, with standardized mean differences � 0.1 ( Table 1) . Use of opioids was strongly associated with use of NSAIDs and 30-day mortality, while cardiovascular disease and dementia was negatively associated with use of NSAIDs and positively associated with death (S1 Table) .

NSAID use was not associated with 30-day mortality in the crude (unmatched) analyses or adjusted (matched) analyses (Table 2 ). In the adjusted analyses, the 30-day mortality rate was 6.3% (95% CI 3.1% to 9.4%) in NSAID users and 6.1% (95% CI 4.4% to 7.8%) in non-users, 

In the crude analyses, use of NSAIDs was associated with an increased risk of hospitalization (RR 1.42, 95% CI 1.13 to 1.79, p < 0.01) but was not associated with an increased risk of ICU 

The subgroup analyses were limited by low power and wide confidence intervals (Table 3 ). In individuals below 65 years of age, the adjusted RR for death associated with NSAID use was NSAID use was defined as having an NSAID prescription filled within 30 days prior to the date of cohort entry. � Patients with a secondary outcome occurring during the exclusion assessment window were excluded, resulting in exclusion of n = 594 patients for hospitalization, n = 33 for ICU admission, n = 18 for mechanical ventilation, and n = 6 for renal replacement therapy in unmatched cohorts, and n = 90, n = 8, n = 5, and n < 5, respectively, in matched cohorts. �� Censored to preserve anonymity for counts n < 5.

ICU, intensive care unit; NSAID, non-steroidal anti-inflammatory drug.

https://doi.org/10.1371/journal.pmed.1003308.t002 Table 3 . 

Increasing the duration of follow-up to 60 days for death and 30 days for hospitalization, ICU admission, mechanical ventilation, and acute renal replacement therapy did not influence the findings (S2 Table) . When extending the definition of current NSAID use to a filled prescription up to 60 days before the test date, we observed similar null findings (S3 Table) . Likewise, defining exposure as a filled prescription 60 days to 14 days before the test date yielded comparable results (S4 Table) .

We identified 204,920 individuals with a negative SARS-CoV-2 PCR test in the study period. Of these, 3,506 were excluded due to migration within 1 year prior to cohort entry, resulting in a population of 201,414 individuals followed up for a total of 15,840 person-years. Use of NSAIDs was associated with a decreased risk of death (RR 0.64, 95% CI 0.49 to 0.84, p < 0.01) and increased risk of hospitalization (RR 1.18, 95% CI 1.08 to 1.28, p < 0.001) and ICU admission (RR 1.39, 95% CI 1.00-1.95, p = 0.05) in the adjusted analyses (Table 4 ).

We examined whether use of NSAIDs was associated with 30-day mortality and adverse outcomes in a nationwide population of SARS-CoV-2-positive individuals. Use of NSAIDs was not associated with increased 30-day mortality, a finding that was robust in a range of supplementary analyses. Likewise, use of NSAIDs was not associated with an increased risk of hospitalization, ICU admission, mechanical ventilation, or renal replacement therapy in the adjusted analyses. 

The Danish nationwide registries allowed for identification of all individuals who had been tested for SARS-CoV-2 in Denmark and allowed for obtaining data on prescription drug use, medical history, migration, hospital admissions, and death through individual-level linkage between health and administrative registries. Thus, we were able to include SARS-CoV-2-positive individuals managed in the community, unlike other data sources that mainly record hospitalized patients. Further, the SARS-CoV-2 test-negative population allowed us to conduct a negative control analysis.

The main limitations of our study are potential misclassification of non-users as NSAID users; potentially unmeasured confounding due to a lack of information on study participants' bodyweight, an important risk factor of severe COVID-19; and confounding by indication due to NSAIDs possibly being prescribed due to early symptoms of severe COVID-19.

In Denmark, NSAIDs can only be obtained via prescription, except for low-dose (200 mg) ibuprofen. With the limited availability and use of NSAIDs over-the-counter in Denmark, misclassification from over-the-counter purchases of NSAIDs is of less concern [19] . However, exposure misclassification may still be present since information on adherence, intended duration, and dose was not available. Some patients do not use NSAIDs continuously after filling a prescription, and consequently some of these patients will have been incorrectly classified as NSAID users at the time of the positive test for SARS-CoV-2. The fact that a patient had filled a prescription for an NSAID may be considered an indicator of availability of NSAIDs rather than of actual use of NSAIDs. As information on exposure status was collected prior to our outcomes of interest, the misclassification is nondifferential and may introduce a bias towards the null. This could lead to wrongly dismissing a causal detrimental effect of NSAIDs on the prognosis of COVID-19 if many NSAID users truly were non-users. The time of cohort entry was defined by the SARS-CoV-2 test date because information on time of symptom onset was not available. Thus, the timing of NSAID use relative to cohort entry will not necessarily reflect NSAID use in the early course of COVID-19 disease.

Users of NSAIDs were more likely to be overweight or obese than non-users, both among individuals who tested positive for SARS-CoV-2 and those who tested negative. This is possibly explained by the fact that a diagnosis of overweight or obesity in the Danish National Patient Registry is dependent on a hospital admission or outpatient clinic visit. This is also a likely explanation for the relatively low prevalence of overweight or obesity in this study. The positive predictive value of these diagnoses in the registries is, however, high [28] .

Use of NSAIDs has been associated with lower mortality among elderly individuals, presumably because of confounding by indication [29] , i.e., NSAIDs are preferentially prescribed to younger, less frail patients because of the established renal, gastrointestinal, and cardiovascular adverse effects [30] . Considerable media attention to the hypothesized risks of use of NSAIDs in COVID-19 may also have influenced how physicians prescribed NSAIDs to selected patients. Finally, severe symptoms early in the course of COVID-19 disease, before the patient is known to be infected with SARS-CoV-2, may increase the likelihood of being prescribed NSAIDs, which would bias the effect estimates towards an increased risk of severe disease associated with NSAIDs.

The secondary outcomes reflecting in-hospital treatment decisions may be more prone to confounding by indication because of the clinical selection of patients to be hospitalized or admitted to the ICU. For example, besides disease severity, factors such as age, comorbidity, and expected outcomes are involved in the ICU triage decision.

Ideally, confounding by indication would be mitigated using an active comparator; however, a suitable active comparator does not exist for ibuprofen. In a previous study, users of paracetamol differed more from NSAID users than did non-users of NSAIDs [9] .

The lack of association between use of NSAIDs and adverse outcomes in individuals who tested positive for SARS-CoV-2 in our study could be explained by multiple factors. First, NSAIDs may not increase angiotensin converting enzyme 2 in humans. The original hypothesis stems from experiments conducted in diabetic rats [31] , and the findings may not transfer between organisms. Second, increased angiotensin converting enzyme 2 expression may not affect the risk of severe COVID-19. Studies on angiotensin converting enzyme inhibitors and angiotensin receptor blockers and the risk of contracting COVID-19 or experiencing a severe course of disease have not shown any association [32] [33] [34] . Third, the adverse effects of NSAIDs on the course of pneumonia may be specific to bacterial infections. A recent study on the risk of adverse outcomes in users of NSAIDs hospitalized for influenza found no association between use of NSAIDs and ICU admission or death [9] , similar to the findings in this study.

Considering the available evidence, there is no reason to withdraw well-indicated use of NSAIDs during the SARS-CoV-2 pandemic. However, the well-established adverse effects of NSAIDs, particularly their renal, gastrointestinal, and cardiovascular effects, should always be considered, and NSAIDs should be used in the lowest possible dose for the shortest possible duration for all patients [30] .

Use of NSAIDs was not associated with an increased risk of 30-day mortality or adverse outcomes in patients infected with SARS-CoV-2 in this cohort of all Danish residents who tested positive for SARS-CoV-2. Table. Association between current NSAID use and mortality, hospitalization, ICU admission, mechanical ventilation, and renal replacement therapy using an outcome assessment window of 60 days for mortality and 30 days for secondary outcomes. (DOCX) S3 Table. Association between use of NSAIDs within 60 days before cohort entry and 30-day mortality, hospitalization, ICU admission, mechanical ventilation, and renal replacement therapy. (DOCX) S4 Table. Association between NSAID prescription fills within the period 60 days to 14 days before cohort entry and 30-day mortality, hospitalization, ICU admission, mechanical ventilation, and renal replacement therapy. 

Are patients with hypertension and diabetes mellitus at increased risk for COVID-19 infection?

Covid-19: ibuprofen should not be used for managing symptoms, say doctors and scientists

EMA gives advice on the use of non-steroidal anti-inflammatories for COVID-19. EMA/136850/2020. Amsterdam: European Medicines Agency

Nonsteroidal antiinflammatory drug use and clinical outcomes of community-acquired pneumonia

Risks related to the use of nonsteroidal anti-inflammatory drugs in community-acquired pneumonia in adult and pediatric patients

Predicting complicated parapneumonic effusion in community acquired pneumonia: hospital based case-control study

Risks of nonsteroidal antiinflammatory drugs in undiagnosed intensive care unit pneumococcal pneumonia: younger and more severely affected patients

Nonsteroidal anti-inflammatory drug and aspirin use, and mortality among critically ill pandemic H1N1 influenza patients: an exploratory analysis

Association of nonsteroidal anti-inflammatory drug use and adverse outcomes among patients hospitalized with influenza

Use of non-steroidal anti-inflammatory drugs and clinical outcome of COVID-19: a Danish nationwide cohort study. EUPAS34734. Amsterdam: European Network of Centres for Pharmacoepidemiology and Pharmacovigilance

Introduction to Danish (nationwide) registers on health and social issues: structure, access, legislation, and archiving

Existing data sources in clinical epidemiology: the Danish COVID-19 cohort

The Danish Microbiology Database (MiBa) 2010 to 2013

The Danish health care system and epidemiological research: from health care contacts to database records

Data resource profile: the Danish National Prescription Registry

The Danish National Patient Registry: a review of content, data quality, and research potential

The Danish Register of Causes of Death

Validity of the coding for intensive care admission, mechanical ventilation, and acute dialysis in the Danish National Patient Registry: a short report

Potential of prescription registries to capture individual-level use of aspirin and other nonsteroidal anti-inflammatory drugs in Denmark: trends in utilization 1999-2012

Data resource profile: Danish online drug use statistics (MED-STAT)

Nonaspirin nonsteroidal antiinflammatory drug use in the Nordic countries from a cardiovascular risk perspective, 2000-2016: a drug utilization study

The central role of the propensity score in observational studies for causal effects

Treatment effects in the presence of unmeasured confounding: dealing with observations in the tails of the propensity score distribution-a simulation study

PSMATCH2: Stata module to perform full Mahalanobis and propensity score matching, common support graphing, and covariate imbalance testing. Boston: Boston College Department of Economics

Considerations for pharmacoepidemiological analyses in the SARS-CoV-2 pandemic

Metrics for covariate balance in cohort studies of causal effects: metrics for covariate balance in cohort studies of causal effects

Applying propensity scores estimated in a full cohort to adjust for confounding in subgroup analyses

Validity of ICD-10 diagnoses of overweight and obesity in Danish hospitals

Paradoxical relations of drug treatment with mortality in older persons

Cardiovascular safety of non-aspirin non-steroidal anti-inflammatory drugs: review and position paper by the working group for Cardiovascular Pharmacotherapy of the European Society of Cardiology

Ibuprofen attenuates cardiac fibrosis in streptozotocin-induced diabetic rats

Association of use of angiotensinconverting enzyme inhibitors and angiotensin II receptor blockers with testing positive for coronavirus disease 2019 (COVID-19)

Renin-angiotensin-aldosterone system blockers and the risk of Covid-19

Renin-angiotensin-aldosterone system inhibitors and risk of Covid-19

The departments of clinical microbiology throughout Denmark are acknowledged for contributing to the national infectious disease surveillance. The Danish Health Data Authority and Statens Serum Institut are acknowledged for valuable support with preparation and linkage of data.