Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Effect of common maintenance drugs on the risk and severity of COVID-19 in elderly patients

  • Kin Wah Fung ,

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

    kfung@mail.nih.gov

    Affiliation Lister Hill National Center for Biomedical Communications, National Library of Medicine, National Institutes of Health, Bethesda, Maryland, United States of America

  • Seo H. Baik,

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

    Affiliation Lister Hill National Center for Biomedical Communications, National Library of Medicine, National Institutes of Health, Bethesda, Maryland, United States of America

  • Fitsum Baye,

    Roles Formal analysis, Writing – review & editing

    Affiliation Lister Hill National Center for Biomedical Communications, National Library of Medicine, National Institutes of Health, Bethesda, Maryland, United States of America

  • Zhaonian Zheng,

    Roles Formal analysis, Writing – review & editing

    Affiliation Lister Hill National Center for Biomedical Communications, National Library of Medicine, National Institutes of Health, Bethesda, Maryland, United States of America

  • Vojtech Huser,

    Roles Investigation, Writing – review & editing

    Affiliation Lister Hill National Center for Biomedical Communications, National Library of Medicine, National Institutes of Health, Bethesda, Maryland, United States of America

  • Clement J. McDonald

    Roles Conceptualization, Methodology, Writing – review & editing

    Affiliation Lister Hill National Center for Biomedical Communications, National Library of Medicine, National Institutes of Health, Bethesda, Maryland, United States of America

Abstract

Background

Maintenance drugs are used to treat chronic conditions. Several classes of maintenance drugs have attracted attention because of their potential to affect susceptibility to and severity of COVID-19.

Methods

Using claims data on 20% random sample of Part D Medicare enrollees from April to December 2020, we identified patients diagnosed with COVID-19. Using a nested case-control design, non-COVID-19 controls were identified by 1:5 matching on age, race, sex, dual-eligibility status, and geographical region. We identified usage of angiotensin-converting enzyme inhibitors (ACEI), angiotensin-receptor blockers (ARB), statins, warfarin, direct factor Xa inhibitors, P2Y12 inhibitors, famotidine and hydroxychloroquine based on Medicare prescription claims data. Using extended Cox regression models with time-varying propensity score adjustment we examined the independent effect of each study drug on contracting COVID-19. For severity of COVID-19, we performed extended Cox regressions on all COVID-19 patients, using COVID-19-related hospitalization and all-cause mortality as outcomes. Covariates included gender, age, race, geographic region, low-income indicator, and co-morbidities. To compensate for indication bias related to the use of hydroxychloroquine for the prophylaxis or treatment of COVID-19, we censored patients who only started on hydroxychloroquine in 2020.

Results

Up to December 2020, our sample contained 374,229 Medicare patients over 65 who were diagnosed with COVID-19. Among the COVID-19 patients, 278,912 (74.6%) were on at least one study drug. The three most common study drugs among COVID-19 patients were statins 187,374 (50.1%), ACEI 97,843 (26.2%) and ARB 83,290 (22.3%). For all three outcomes (diagnosis, hospitalization and death), current users of ACEI, ARB, statins, warfarin, direct factor Xa inhibitors and P2Y12 inhibitors were associated with reduced risks, compared to never users. Famotidine did not show consistent significant effects. Hydroxychloroquine did not show significant effects after censoring of recent starters.

Conclusion

Maintenance use of ACEI, ARB, warfarin, statins, direct factor Xa inhibitors and P2Y12 inhibitors was associated with reduction in risk of acquiring COVID-19 and dying from it.

1. Background and significance

Maintenance drugs are indicated for chronic conditions, taken indefinitely on regular, usually daily basis. Reports have attributed protective or aggravating effects of some maintenance drugs on the likelihood and/or severity of SARS-CoV-2 infection. These include angiotensin-converting enzyme inhibitors (ACEIs), angiotensin-receptor blockers (ARBs), statins, antithrombotic agents, hydroxychloroquine, and famotidine.

ACEI and ARB attracted early attention because they increase the expression of angiotensin-converting enzyme 2 (ACE2) [1, 2], the cellular doorway for SARS-CoV-2 [3, 4]. There are concerns that their use could increase the likelihood, or the severity, of a COVID-19 infection. However, most studies published to-date show no such increases [59], and some actually suggest that they might decrease COVID-19 severity [10, 11]. Due to the potential beneficial effect, several ongoing prospective trials are testing whether initiating ACEI or ARB treatment for patients after they are diagnosed with COVID-19 would confer any benefits [1214].

Statins are among the commonest maintenance drugs in the elderly population. Attention was drawn to statins because of their potential effect on virus entry, replication or degradation, and their well-known anti-inflammatory properties [15, 16]. Several subsequent studies have shown the beneficial effects of statins in COVID-19 patients [17, 18].

Given the thrombotic propensity of COVID-19 [1923], some think that antithrombotics could reduce the severity of COVID-19, and accordingly, anticoagulants are now recommended for hospitalized and high risk patients [24], though results of studies of maintenance antithrombotics have been mixed [2529]. Famotidine has attracted attention because it was shown to have potential anti- SARS-CoV-2 activity based on structural homology modeling [30], and it improved COVID-19 outcomes in small studies [3133]. Hydroxychloroquine was considered as treatment option for COVID-19 because of its immune modulation and anti-SARS-CoV-2 activity in vitro [3436], and reported benefit in one clinical study [37].

Most of the early primary studies on the effect of these maintenance drugs have been small (tens or hundreds of COVID-19 patients who were also on the drugs) and would most likely lack statistical power to see important associations.

In the U.S., the Virtual Research Data Center (VRDC) [38] of the Centers for Medicare and Medicaid Services (CMS) carries de-identified data which include complete drug prescription, diagnoses, encounters as well as geographic, socio-demographic and vital status information for most (>93%) of the 65 and older U.S. residents. Over 80% of COVID-19 deaths occur in this age group, according to the U.S. Centers for Disease Control and Prevention (CDC). Therefore, this would be a propitious population in which to explore the association of maintenance drugs with COVID-19 risks.

In this study, we used the VRDC data and extended Cox regression analysis to examine associations between the above eight drug types and the occurrence of three outcomes: 1) COVID-19 infection, 2) COVID-19 hospitalization, and 3) death after a COVID19 diagnosis.

2. Materials and methods

2.1 Study population and case definition

VRDC provided us with a 20% random sample of all Medicare Part D enrollees. This set included 374,299 patients 65 or above, who had at least one record of the COVID-19 specific ICD10-CM diagnosis code of U07.1 between April 1 and December 31, 2020. We used both inpatient and outpatient claims and Medicare’s vital status to identify COVID-19 cases and the outcome events. We only counted COVID-19 cases occurring on or after April 1, 2020, when the specific code for COVID-19 first became available. We stopped accruing cases after December 31, 2020, because of potential incompleteness of data due to the time lag (at least 3–4 months) between data capture and availability through VRDC. This study was declared not human subject research by the Office of Human Research Protection at the National Institutes of Health and by the CMS’s Privacy Board.

2.2 Drugs, exposure definition and comorbidities

We identified all study drug preparations that were available in the U.S. market. For drug classes we identified the class members through Anatomical Therapeutic Chemical (ATC) classification codes and used generic drug names to find them in the CMS data. Our eight study drugs/drug classes were ACEI, ARB, statins, warfarin, direct factor Xa inhibitors, P2Y12 inhibitors, famotidine, and hydroxychloroquine (see S1 Table for full list of drugs and prescription frequencies). H2 blockers other than famotidine were not included because the potential beneficial effect was based on the chemical structure of famotidine and not the pharmacological class. Note that hydroxychloroquine was a class which also included two other aminoquinolines (chloroquine and primaquine), but hydroxychloroquine accounted for 99.7% of prescriptions.

We assumed that patients were on a given study drug during the window from the prescription dispensing date to 30 days after the end-of-supply day (we called this period the “current use period”). We added the 30-day buffer after the end-of-supply day because of the common behavior of drug stockpiling—patients maintaining a stock of drugs so that they will not run out immediately in case refill is interrupted. Drug use status was treated as a time-varying covariate in the Cox regression model. We defined a patient to be a current user if an outcome event fell within a current use period, former user if the event fell outside the current use period, and never user if they never had a prescription for that drug. Our primary analysis compared current users with never users as control. In a supplementary analysis, we used former users as control to see if that would give different results.

In the first few months of the pandemic, there were media reports of the potential beneficial or harmful effects of some of the study drugs on COVID-19. Publicity about the good effects of hydroxychloroquine was particularly rife, which could have led people to start taking it because of symptoms, or fear, of COVID -19. We looked at the trend of the usage of the study drugs, starting from January 2019, to identify abnormal patterns that could be related to COVID-19. To study the potential effect of COVID-19 affecting the use of drugs (“reverse causality”), we did separate analyses with special treatment of patients who were only recently started on a study drug (see section 2.3 below).

Medicare specifies the onset of 67 chronic conditions by algorithm and we followed the algorithm to define the occurrence and onset date of each condition [39]. To adjust for illness burden in our analysis we included 57 Medicare chronic conditions with >1% prevalence in the Master Beneficiary Summary File.

2.3 Statistical analysis and covariates

We excluded patients who were not enrolled in Parts A (hospital) and B (medical) to ensure complete capture of hospitalization and encounter diagnoses. We considered the effect of the study drugs on three outcomes: 1) the risk of acquiring a COVID-19 diagnosis, 2) the risk of COVID-19 hospitalization, and 3) the risk of death after being diagnosed with COVID-19. For the first outcome (diagnosis of COVID-19), we used a nested case-control analysis [40], where each index COVID-19 patient was matched to patients with no COVID-19 diagnosis up to the date COVID-19 was diagnosed in the index patient.

For the first outcome (COVID-19 diagnosis), each index case was matched to five controls on age in years, race, sex, dual-eligibility status and five regions of residence (Northeast, Midwest, South, West, and others) at the time of the diagnosis of COVID-19. All cases and control were followed from January 1, 2020, until COVID-19 diagnosis, death, disenrollment from Medicare Parts A/B/D or December 31, 2020, whichever came first.

For the analyses of the second and third outcomes (hospitalization and death), we included only patients diagnosed with COVID-19. For the hospitalization outcome, we followed all patients from COVID-19 diagnosis until COVID-19-related hospitalization, or the other censoring points (except COVID-19 diagnosis) as described for the first outcome. For the death outcome, we did the same, swapping hospitalization for death.

We approximated patient’s income level using the monthly indicators of dual-eligibility and low-income subsidy (LIS), which divided patients into three income groups: 1) dual-eligible: income ≤ 135% federal poverty line (FPL), 2) non-dual LIS: income > 135% and ≤ 150% FPL, and 3) non-dual no LIS: income > 150% FPL. To explore the effect of each study drug, we employed extended Cox regression analysis with days-on-study as time scale. We included age, gender, race, regions of residence, Medicare insurance type (Advantage or fee-for-service) and the degree of LIS as covariates. We also included binary flags for each of the 57 Medicare chronic conditions to control for the effect of co-morbidities. In order to protect against the immortal time bias and violation of proportional hazard assumption [41], we treated the use of each of the study drugs and almost all covariates as time-varying covariates. Only age, sex, race, and regions of residence were time-fixed, and their values were defined as of their values on January 1, 2020. Since the disease covariates were chronic diseases, we considered them as always present after their onset and always absent before then. The values of all time-varying covariates were reset at the time of each event in the Cox regression. We used Efron’s adjustment for tied events.

To mitigate the potential of COVID-19 influencing the use of drugs, we did a separate Cox regression in which patients who only started a study drug in 2020 (they had never been on the drug before January 1,2020) were censored when they first started the drug (“recent starter censoring”). We compared the results with and without such censoring.

In order to mitigate selection bias toward use of the study drugs, we developed time-varying propensity scores (PS) [42] separately for each study drug using logistic regressions. The PS was the likelihood of receiving a study drug, conditional on patient’s characteristics (demographics, socioeconomics, and presence of the chronic conditions). We iteratively estimated the PSs every month among the patients who remained in follow-up, considering all covariates that preceded the end of a given monthly cycle [43], and ran all Cox regression analyses with time-varying PSs as additional adjustments.

3. Results

3.1 Study population

In our 20% random sample of Part D enrollees, 374,299 Medicare beneficiaries aged 65 and above were diagnosed with COVID-19 between April 1 and December 31, 2020, among which 65,108 (17.4%) died (Table 1 and S1 Fig). Over the same period, CDC registered a national total of 19,852,636 COVID-19 cases and 350,510 deaths [44]. In the CDC statistics, even though patients ≥65 accounted for only 13.9% of cases, they accounted for 80.5% of mortalities. Therefore, our study population represents patients most at risk for COVID-19 death. Projecting from the CDC statistics, a 20% random sample of patients 65 and above would have 551,903 COVID-19 cases and 56,432 mortalities (Table 1). Therefore, our study population captured 67.8% of COVID-19 cases nationally. Our mortality number is 15% higher than the projected national count. Our death numbers were higher because Medicare data did not distinguish between causes of death. It is possible that some mortalities were not related to COVID-19. However, most of the mortalities occurred within a short time after the COVID-19 diagnosis (median 14 days, inter-quartile range 6–36 days). Compared with national statistics, the Medicare population was over-represented in females, whites and under-represented in Hispanics. Geographically, there was over-representation of the North-East region.

thumbnail
Table 1. Study population compared to U.S. national statistics, April to December 2020 (SRS—Simple random sample; CDC—Centers for Disease Control and Prevention).

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

3.2 Drug exposure and matching

Table 2 shows the breakdown of drug usage (all patients who had prescription for the drug during the follow up period, starting on January 1, 2020) among the COVID-19 patients and controls. Overall, 278,912 (74.6%) cases and 1,351,244 (72.4%) controls had prescriptions for at least one study drug. The three most common study drugs among COVID-19 patients were statins 187,374 (50.1%), ACEI 97,872 843 (26.2%) and ARB 83,290 (22.3%). For all study drugs except ACEI and ARB, the usage was significantly higher among COVID-19 patients than controls. The unmatched rate (less than five controls found) was 0.4%.

thumbnail
Table 2. Characteristics of Covid-19 patients and matched controls.

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

As for abnormal patterns of drug usage related to COVID-19, we detected a sharp rise in the use of hydroxychloroquine in March and April 2020 (S2 Fig), which was not seen in other study drugs. FDA granted emergency use authorization for hydroxychloroquine for COVID-19 treatment on March 28, and the percentage of COVID-19 patients on hydroxychloroquine surged to over 0.8% compared to the baseline of 0.3% before the pandemic. There was a similar but smaller rise in non-COVID-19 patients, showing that some patients might be taking hydroxychloroquine for prophylaxis or suspicion of COVID-19. After FDA revoked the emergency use authorization on June 15, the usage of hydroxychloroquine began to drop, but remained slightly above the baseline before the pandemic.

3.3 Risk of being diagnosed with COVID-19

The total number of patients followed (cases and controls) was about 2.2 million, and around 350,000 ended up with the diagnosis of COVID-19 (Table 3). We did separate analyses with and without recent starter censoring for all study drugs. Only hydroxychloroquine exhibited significantly different results with censoring. We report the results with recent starter censoring for hydroxychloroquine as the main results in Table 3. The results for hydroxychloroquine without such censoring are shown for comparison.

thumbnail
Table 3. Effect of drug use on Covid-19 diagnosis, Covid-19 hospitalization and death (other non-drug covariates are listed in Table 4).

https://doi.org/10.1371/journal.pone.0266922.t003

Compared to patients who never used the drug, current users of ACEI, ARB, statin, warfarin, direct factor Xa inhibitors and P2Y12 inhibitors were associated with decreased risk of contracting COVID-19, ranging from 1% reduction (direct factor Xa inhibitors) to 12% reduction (warfarin). For famotidine, there was a 12% increase in risk. For hydroxychloroquine, with the proper censoring of recent starters, there was no effect on the risk of being diagnosed with COVID-19. However, without censoring, current users of hydroxychloroquine would appear to be associated with an increased risk of getting COVID-19 compared to never users (63% increase).

The results for the non-drug covariates are shown in Table 4. Compared with the youngest group of patients in our population (65–69), older age was associated with decreased risk of COVID-19 diagnosis. Female sex was associated with a 11% risk reduction compared to male. The risk of COVID-19 diagnosis was higher in all race groups compared to whites, 13% higher for Asians, 54% for Blacks and 85% for Hispanics. The comorbidities associated with the greatest increase in risk were dementia (132% increase), hypertension (125% increase) and schizophrenia (68% increase).

thumbnail
Table 4. Effect of non-drug covariates on Covid-19 diagnosis, Covid-19 hospitalization and death (drug-related covariates are listed in Table 3).

https://doi.org/10.1371/journal.pone.0266922.t004

3.4 Risk of COVID-19 hospitalization

Overall, 142,004 (39.6%) COVID-19 patients were hospitalized for COVID-19. Current users of ACEI, ARB, statin, warfarin, direct factor Xa inhibitors, P2Y12 inhibitors and famotidine were associated with 2–11% decreased risk of COVID-19 hospitalization when compared with never users (Table 3). Hydroxychloroquine did not show any effect when recent starters were censored. The risk of COVID-19 hospitalization increased monotonically with age (Table 4). Compared to the 65–69 age group, the risk increased by 10%, 40%, 63% and 95% for groups 70–74, 75–79, 80–84 and ≥85 respectively. Female sex was associated with a 20% reduced risk compared to male. The risk of COVID-19 hospitalization for Hispanics was 19% higher than whites. Poorer patients were associated with higher risks of hospitalization with COVID-19. The three co-morbidities associated with largest increase in risk were hypertension (134%), chronic kidney disease (34%) and dementia (32%).

3.5 Risk of mortality

Overall, 61,778 (17.2%) COVID-19 patients died. Current users of ACEI, ARB, statin, warfarin, direct Xa inhibitors and P2Y12 inhibitors among COVID-19 patients were all associated with lower risks of mortality, from 6% less than never users (P2Y12 inhibitors) to 20% less (direct factor Xa inhibitors) (Table 3). Famotidine and hydroxychloroquine showed no significant effect. Mortality risk increased monotonically with age (Table 4). Compared to patients 65–69, all other age groups exhibited increased risk of death: 12%, 55%, 103% and 233% increased risk among 70–74, 75–79, 80–84 and ≥85 respectively. Female sex was associated with a 35% reduced risk of death compared to male. All race groups were associated with notably greater mortality risk than whites. Unlike analyses of other outcomes, the poorest patients (dual-eligible) were associated with 36% decreased risk of death, while the second poorest patients (non-dual LIS) had 206% increased risk of death compared to patients not dually eligible and not receiving LIS. The three co-morbidities associated with the greatest mortality risk were dementia (127%) chronic kidney disease (56%) and lung cancer (53%).

4. Discussion

The most important finding of our study is that after controlling for age, gender, race, socio-economic, geographic factors and co-morbidities, there was an associated decline in mortality risk of 12% or greater among current users for ACEIs, ARB, statins, and anticoagulants, and by a lesser amount (6%) among P2Y12 inhibitors users. To the best of our knowledge, ours is first single source study to show such beneficial associations, though it has been hinted at in meta-analyses [10, 11]. While more study is needed to identify the exact reasons for severity and mortality modification of the study drugs in COVID-19 patients, some possible mechanisms have been proposed. One predominant hypothesis is that SARS-CoV-2 down-regulates ACE2 expression, resulting in unabated angiotensin II activity that may be responsible for organ damage in COVID-19 [6]. ACEI and ARB reduce angiotensin II activity and so are protective [45]. Our findings provide strong support to the advice from professional societies and the WHO that ACEI and ARB be continued in COVID-19 patients [46, 47]. The benefits of antithrombotics can be explained by reduction in the thromboembolism observed in COVID-19 patients [21, 23, 48]. These results strongly support the current advice about the use of anti-clotting drugs in COVID-19 patients. The benefits of statins can be associated with their anti-viral and anti-inflammatory properties [15, 16]. In addition, our study gives justification for the need of clinical trials to initiate treatment with drugs that are potentially beneficial (such as ACEI, ARB and statins) in patients diagnosed with COVID-19.

The main strength of this study is the size of the study population and relative completeness of demographic, prescription, and co-morbidity data. The U.S. has the largest collection of COVID-19 cases in the world. The elderly represents the most at-risk population for severe COVID-19. Our study population covered over 370,000 COVID-19 cases and 65,000 mortalities in U.S. patients over 65. According to CDC statistics, over 80% of COVID-19 mortalities occurred in patients over 65. A large U.K. study based on NHS records registered 10,926 COVID-19 deaths, our study surpassed that number by almost six-fold [49]. In terms of study drug usage among COVID-19 patients, our numbers were also much higher. For instance, one of the biggest early original ACEI studies was from Italy, covering 1,502 ACEI users among COVID-19 patients [7]. We had 65 times that number. Recently published studies do have larger patient numbers. One such study supported by the American Heart Association’s Rapid Response Grant COVID-19 is by An et al. [50], with 4,652 patients on ACEI and 2,546 on ARB. However, this is still lower than the number of patients observed in our study, 97,843 on ACEI and 83,290 on ARB. Sample sizes of this magnitude are not seen in other studies, including meta-analyses [8, 10, 11]. Moreover, other studies tended to be restricted to institutional or regional populations. Our Medicare population includes most US elderly individuals. By one estimate 93% of all US adults over 65 are enrolled to Medicare [51]. Unlike meta-analysis that pools multiple, not necessarily comparable data sources, we used a single data source that has relatively complete and longitudinal medication information with well-documented data capture procedures.

The emergence of hydroxychloroquine as a potential “game-changer” in COVID-19 had been accompanied by the dramatic rise in its use, not only among COVID-19 patients but in the general population as well. Even after the emergency use authorization was revoked by FDA, the use of hydroxychloroquine remained somewhat higher than the baseline level before the pandemic. This shows that media claims, even those that are eventually shown to be unsubstantiated or disproved, may have lasting impact on the public psyche. In our observational study, the use of hydroxychloroquine is a classic example of indication bias—the indication for the exposure is directly related to the outcome being observed [52]. Since we have longitudinal prescription data (starting from 2019), we were able to compensate for this by censoring patients who were first started on hydroxychloroquine during the pandemic. Indeed, our results show that, without censoring, the use of hydroxychloroquine would have appeared to be associated with an increase in the risk of COVID-19 diagnosis and death, but the effects disappeared with proper censoring.

Regarding the risk of catching COVID-19, we found that ACEI, ARB, statins and antithrombotics were associated with a reduction in the risk of getting a COVID-19 diagnosis. One possible explanation is that if these drugs blunted the effects of COVID-19 infection, infected patients might have milder symptoms and would be more likely to not seek medical care and remain undiagnosed. The apparent reduction in risk of catching COVID-19 with advancing age could be related to the reduced level of social activity at extreme age, thus lowering the risk of exposure.

Apart from drug usage, we found the following risk factors for severe COVID-19—advanced age, male, non-white, and co-morbidities including chronic kidney diseases, dementia, hypertension, heart diseases, chronic obstructive pulmonary disease, chronic liver diseases and some malignancies. These findings concur with similar studies and are well-documented. In our study, the group with the lowest income (dual-eligible) was associated with a reduced risk of death. One possible explanation is that this group has better access to healthcare, since patients are eligible to both Medicaid and Medicare. However, this explanation only applies to the Medicaid and Medicare patients and not low income patients in general.

We recognize the following limitations. In our retrospective observation study, drug exposure was implied from prescription records and there was no verification that the drugs were dispensed or taken. Drug data from Medicare claims are incomplete for hospitals and nursing homes. As a result, the continuation of medications for inpatients could not be ascertained. For certain drugs (e.g., statins), it has been reported that rebound effects were associated with drug discontinuation which could potentially affect the outcomes of COVID-19 infection [53]. Some of the study drugs are often prescribed in combination with other drugs. For example, statins are often used together with anti-hypertensive medications. In our Medicare population, statins are taken together with ACEI (in 35% of patients on statins), ARB (28.6%), beta blockers (44.3%), calcium channel blockers (37.7%) and diuretics (39.8%). (S2 Table) This can be a potential confounder if not properly addressed. In our multiple variable analysis, the use of ACEI, ARB and statins are included in the same unified regression model as time-varying covariates. This ensures that the observed effects of each individual study drug are already adjusted for the co-prescription of other drugs included in the model [54]. It is true that the combined use of statins with other anti-hypertensive drugs (e.g., diuretics, beta blockers) has not been adjusted for. However, up to the start of the study, there was no strong evidence that diuretics, calcium channel blockers or beta blockers significantly affect the susceptibility to and severity of COVID-19. It is therefore unlikely that the observed effects of the study drugs on COVID-19 are caused by these other drugs. Furthermore, there could also be other residual confounding factors outside the scope of our regression models. We did not include COVID-19 patients before April because the COVID-19 specific diagnosis code was not in use. Compared with national data from CDC, we could be missing about one-third of COVID-19 cases. It is likely that milder cases would be missed since the patients did not seek medical care. This would also explain the relatively high hospitalization and mortality rate in our population compared to some other studies [55]. We used all-cause mortality because of the lack of an accurate cause of death in CMS data. Because of this, it is possible that some reduction of mortality is due to the general protective effect of the study drugs and not related to COVID-19 per se (e.g., statins and hypotensive agents can reduce cardiovascular mortality). However, since the mortality rate was 17% in our COVID-19 patients overall, 29% among those hospitalized, and death occurred within a median of 14 days of COVID-19 diagnosis, it is highly likely that COVID-19 was the main contributor to mortality. Another potential confounding factor is that some patients could be diagnosed with COVID-19 because of more frequent routine testing during an encounter for another reason (e.g., preparation for coronary bypass), which makes COVID-19 an incidental finding rather than the trigger event. However, overall, only 3.2% of our Medicare population were diagnosed with COVID-19, so the chance of a purely incidental COVID-19 diagnosis is small.

5. Conclusion

Analysis of more than 370,000 Medicare enrollees over 65 diagnosed with COVID-19 showed that the use of ACEI, ARB, statins, warfarin, direct factor Xa inhibitors and P2Y12 inhibitors was associated with a reduction in the risk of catching COVID-19 and developing severe disease. Hydroxychloroquine and famotidine were not associated with significant effects in these outcomes.

Supporting information

S2 Fig. Trends of drug usage in 2019 and 2020.

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

(PDF)

S1 Table. Drug classes, clinical drugs and usage frequencies among all our study population.

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

(DOCX)

S2 Table. The combination use of anti-hypertensive drugs with statins in Medicare patients for 2019–2020.

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

(DOCX)

S3 Table. Data values for drug drug usage trend graphs in S2 Fig.

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

(XLSX)

S4 Table. Detailed statistical data of the Cox regression models for death, hospitalization and acquiring COVID-19.

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

(XLSX)

Acknowledgments

We would like to thank the diligent and helpful staff at the VRDC and Research Data Assistance Center (ResDAC), without them this study would not be possible.

References

  1. 1. Ferrario C.M., et al., Effect of angiotensin-converting enzyme inhibition and angiotensin II receptor blockers on cardiac angiotensin-converting enzyme 2. Circulation, 2005. 111(20): p. 2605–10. pmid:15897343
  2. 2. Soler M.J., et al., Pharmacologic modulation of ACE2 expression. Curr Hypertens Rep, 2008. 10(5): p. 410–4. pmid:18775121
  3. 3. Li W., et al., Angiotensin-converting enzyme 2 is a functional receptor for the SARS coronavirus. Nature, 2003. 426(6965): p. 450–4. pmid:14647384
  4. 4. Hoffmann M., et al., SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor. Cell, 2020. 181(2): p. 271–280 e8. pmid:32142651
  5. 5. Li J., et al., Association of Renin-Angiotensin System Inhibitors With Severity or Risk of Death in Patients With Hypertension Hospitalized for Coronavirus Disease 2019 (COVID-19) Infection in Wuhan, China. JAMA Cardiol, 2020. 5(7): p. 825–830. pmid:32324209
  6. 6. Vaduganathan M., et al., Renin-Angiotensin-Aldosterone System Inhibitors in Patients with Covid-19. N Engl J Med, 2020. 382(17): p. 1653–1659. pmid:32227760
  7. 7. Mancia G., et al., Renin-Angiotensin-Aldosterone System Blockers and the Risk of Covid-19. N Engl J Med, 2020. 382(25): p. 2431–2440. pmid:32356627
  8. 8. Fosbol E.L., et al., Association of Angiotensin-Converting Enzyme Inhibitor or Angiotensin Receptor Blocker Use With COVID-19 Diagnosis and Mortality. JAMA, 2020. 324(2): p. 168–177. pmid:32558877
  9. 9. Reynolds H.R., et al., Renin-Angiotensin-Aldosterone System Inhibitors and Risk of Covid-19. N Engl J Med, 2020. 382(25): p. 2441–2448. pmid:32356628
  10. 10. Guo X., Zhu Y., and Hong Y., Decreased Mortality of COVID-19 With Renin-Angiotensin-Aldosterone System Inhibitors Therapy in Patients With Hypertension: A Meta-Analysis. Hypertension, 2020. 76(2): p. e13–e14. pmid:32458694
  11. 11. Barochiner J. and Martinez R., Use of inhibitors of the renin-angiotensin system in hypertensive patients and COVID-19 severity: A systematic review and meta-analysis. J Clin Pharm Ther, 2020. pmid:32767823
  12. 12. Efficacy of Captopril in Covid-19 Patients With Severe Acute Respiratory Syndrome (SARS) CoV-2 Pneumonia (CAPTOCOVID) [clinical trial]. https://clinicaltrials.gov/ct2/show/NCT04355429.
  13. 13. Angiotensin Converting Enzyme Inhibitors in Treatment of Covid 19 [clinical trial]. https://clinicaltrials.gov/ct2/show/NCT04345406.
  14. 14. Valsartan for Prevention of Acute Respiratory Distress Syndrome in Hospitalized Patients With SARS-COV-2 (COVID-19) Infection Disease [clinical trial]. https://clinicaltrials.gov/ct2/show/NCT04335786.
  15. 15. Rodrigues-Diez R.R., et al., Statins: Could an old friend help in the fight against COVID-19? Br J Pharmacol, 2020. 177(21): p. 4873–4886. pmid:32562276
  16. 16. Pawlos A., et al., COVID-19: Direct and Indirect Mechanisms of Statins. Int J Mol Sci, 2021. 22(8). pmid:33920709
  17. 17. Permana H., et al., In-hospital use of statins is associated with a reduced risk of mortality in coronavirus-2019 (COVID-19): systematic review and meta-analysis. Pharmacol Rep, 2021. 73(3): p. 769–780. pmid:33608850
  18. 18. Pal R., et al., Statin use and clinical outcomes in patients with COVID-19: An updated systematic review and meta-analysis. Postgrad Med J, 2021. pmid:33541927
  19. 19. Tang N., et al., Abnormal coagulation parameters are associated with poor prognosis in patients with novel coronavirus pneumonia. J Thromb Haemost, 2020. 18(4): p. 844–847. pmid:32073213
  20. 20. Bikdeli B., et al., COVID-19 and Thrombotic or Thromboembolic Disease: Implications for Prevention, Antithrombotic Therapy, and Follow-Up: JACC State-of-the-Art Review. J Am Coll Cardiol, 2020. 75(23): p. 2950–2973. pmid:32311448
  21. 21. Zhang Y., et al., Coagulopathy and Antiphospholipid Antibodies in Patients with Covid-19. N Engl J Med, 2020. 382(17): p. e38. pmid:32268022
  22. 22. Cui S., et al., Prevalence of venous thromboembolism in patients with severe novel coronavirus pneumonia. J Thromb Haemost, 2020. 18(6): p. 1421–1424. pmid:32271988
  23. 23. Klok F.A., et al., Confirmation of the high cumulative incidence of thrombotic complications in critically ill ICU patients with COVID-19: An updated analysis. Thromb Res, 2020. 191: p. 148–150. pmid:32381264
  24. 24. Rico-Mesa J.S., et al., The Role of Anticoagulation in COVID-19-Induced Hypercoagulability. Curr Cardiol Rep, 2020. 22(7): p. 53. pmid:32556892
  25. 25. Tremblay D., et al., Impact of anticoagulation prior to COVID-19 infection: a propensity score-matched cohort study. Blood, 2020. 136(1): p. 144–147. pmid:32462179
  26. 26. Russo V., et al., Clinical impact of pre-admission antithrombotic therapy in hospitalized patients with COVID-19: A multicenter observational study. Pharmacol Res, 2020. 159: p. 104965. pmid:32474087
  27. 27. Sivaloganathan H., Ladikou E.E., and Chevassut T., COVID-19 mortality in patients on anticoagulants and antiplatelet agents. Br J Haematol, 2020. pmid:32584423
  28. 28. Inama G., et al., Coronavirus disease 2019 infection in patients with recent cardiac surgery: does chronic anticoagulant therapy have a protective effect? J Cardiovasc Med (Hagerstown), 2020. 21(10): p. 765–771. pmid:32890069
  29. 29. Rossi R., et al., Protective role of chronic treatment with direct oral anticoagulants in elderly patients affected by interstitial pneumonia in COVID-19 era. Eur J Intern Med, 2020. 77: p. 158–160. pmid:32564905
  30. 30. Wu C., et al., Analysis of therapeutic targets for SARS-CoV-2 and discovery of potential drugs by computational methods. Acta Pharm Sin B, 2020. 10(5): p. 766–788. pmid:32292689
  31. 31. Bourinbaiar A.S. and Fruhstorfer E.C., The effect of histamine type 2 receptor antagonists on human immunodeficiency virus (HIV) replication: identification of a new class of antiviral agents. Life Sci, 1996. 59(23): p. PL 365–70. pmid:8950301
  32. 32. Janowitz T., et al., Famotidine use and quantitative symptom tracking for COVID-19 in non-hospitalised patients: a case series. Gut, 2020. 69(9): p. 1592–1597. pmid:32499303
  33. 33. Freedberg D.E., et al., Famotidine Use Is Associated With Improved Clinical Outcomes in Hospitalized COVID-19 Patients: A Propensity Score Matched Retrospective Cohort Study. Gastroenterology, 2020. 159(3): p. 1129–1131 e3. pmid:32446698
  34. 34. Keyaerts E., et al., In vitro inhibition of severe acute respiratory syndrome coronavirus by chloroquine. Biochem Biophys Res Commun, 2004. 323(1): p. 264–8. pmid:15351731
  35. 35. Vincent M.J., et al., Chloroquine is a potent inhibitor of SARS coronavirus infection and spread. Virol J, 2005. 2: p. 69. pmid:16115318
  36. 36. Sinha N. and Balayla G., Hydroxychloroquine and COVID-19. Postgrad Med J, 2020. 96(1139): p. 550–555. pmid:32295814
  37. 37. Arshad S., et al., Treatment with hydroxychloroquine, azithromycin, and combination in patients hospitalized with COVID-19. Int J Infect Dis, 2020. 97: p. 396–403. pmid:32623082
  38. 38. Research Data Assistance Center (ResDAC). CMS Virtual Research Data Center (VRDC). https://www.resdac.org/cms-virtual-research-data-center-vrdc.
  39. 39. The Chronic Condition Warehouse. Chronic Conditions Data Warehouse: CCW Chronic Condition Algorithms. https://www.ccwdata.org/web/guest/condition-categories.
  40. 40. Ernster V.L., Nested case-control studies. Prev Med, 1994. 23(5): p. 587–90. pmid:7845919
  41. 41. Azzato E.M., et al., Prevalent cases in observational studies of cancer survival: do they bias hazard ratio estimates? Br J Cancer, 2009. 100(11): p. 1806–11. pmid:19401693
  42. 42. Rosenbaum P.R. and Rubin D.B., The Central Role of the Propensity Score in Observational Studies for Causal Effects. Biometrika, 1983. 70: p. 41–55.
  43. 43. Wyss R., et al., Use of Time-Dependent Propensity Scores to Adjust Hazard Ratio Estimates in Cohort Studies with Differential Depletion of Susceptibles. Epidemiology, 2020. 31(1): p. 82–89. pmid:31569120
  44. 44. Centers for Disease Control and Prevention. COVID-19 Response. COVID-19 Case Surveillance Public Data Access, Summary, and Limitations. 2020; https://data.cdc.gov/Case-Surveillance/COVID-19-Case-Surveillance-Public-Use-Data/vbim-akqf.
  45. 45. Elshafei A., et al., RAAS, ACE2 and COVID-19; a mechanistic review. Saudi J Biol Sci, 2021. 28(11): p. 6465–6470. pmid:34305426
  46. 46. Bozkurt B., Kovacs R., and Harrington B., Joint HFSA/ACC/AHA Statement Addresses Concerns Re: Using RAAS Antagonists in COVID-19. J Card Fail, 2020. 26(5): p. 370. pmid:32439095
  47. 47. World Health Organization. COVID-19 and the use of angiotensin-converting enzyme inhibitors and receptor blockers—Scientific Brief. 2020; https://www.who.int/news-room/commentaries/detail/covid-19-and-the-use-of-angiotensin-converting-enzyme-inhibitors-and-receptor-blockers.
  48. 48. Fontelo P., et al., A review of thromboembolic events in hospitalized COVID-19 patients. Thromb J, 2021. 19(1): p. 47. pmid:34187490
  49. 49. Williamson E.J., et al., Factors associated with COVID-19-related death using OpenSAFELY. Nature, 2020. 584(7821): p. 430–436. pmid:32640463
  50. 50. An J., et al., COVID-19 morbidity and mortality associated with angiotensin-converting enzyme inhibitors or angiotensin receptor blockers use among 14,129 patients with hypertension from a US integrated healthcare system. Int J Cardiol Hypertens, 2021. 9: p. 100088. pmid:34155486
  51. 51. Administration on Aging—U.S. Department of Health and Human Services. 2018 Profile of Older Americans. https://acl.gov/aging-and-disability-in-america/data-and-research/profile-older-americans.
  52. 52. Csizmadi I., Collet J.-P., and Boivin J.-F., Bias and Confounding in Pharmacoepidemiology, in Pharmacoepidemiology. 2006. p. 791–809.
  53. 53. Cubeddu L.X. and Seamon M.J., Statin withdrawal: clinical implications and molecular mechanisms. Pharmacotherapy, 2006. 26(9): p. 1288–96. pmid:16945051
  54. 54. Pourhoseingholi M.A., Baghestani A.R., and Vahedi M., How to control confounding effects by statistical analysis. Gastroenterol Hepatol Bed Bench, 2012. 5(2): p. 79–83. pmid:24834204
  55. 55. Verity R., et al., Estimates of the severity of coronavirus disease 2019: a model-based analysis. Lancet Infect Dis, 2020. 20(6): p. 669–677. pmid:32240634