Abstract
Background The novel coronavirus SARS-Cov2 uses the angiotensin-converting enzyme 2 (ACE2) receptor as an entry point to the cell. Cardiovascular disease (CVD) is a risk factor for the novel coronavirus disease (Covid-19) with poor outcomes. We hypothesized that the rate of ACE inhibitor (ACEI) and angiotensin receptor blocker (ARB) use is associated with the rate of Covid-19 confirmed cases and deaths.
Methods We conducted a geospatial study using publicly available county-level data. The Medicare ACEI and ARB prescription rate was exposure. The Covid-19 confirmed case and death rates were outcomes. Spatial autoregression models were adjusted for the percentage of Black residents, children, residents with at least some college degree, median household income, air quality index, CVD hospitalization rate in Medicare beneficiaries, and CVD death rate in a total county population.
Findings After adjustment for confounders, the ACEI use rate did not associate with Covid-19 confirmed case rate (direct county-own effect +0.11 %; 95%CI -0.31 to 0.53; P=0.600, and indirect spillover effect -0.53 %; 95%CI -3.89 to 2.84; P=0.760). The ARB use rate was associated with increased Covid-19 confirmed case rate (direct county-owned effect +0.12 %; 95%CI 0.05-0.19; P=0.002, and indirect spillover effect -0.33 %; 95%CI -2.11 to 1.44; P=0.714). Sensitivity analysis indicated an absence of significant reverse causality bias for analyses with Covid-19 confirmed case rate, but not death rate outcome.
Interpretation Our results highlight the safety of ACEI use for patients with clinical indications for ACEI use. However, an increase in ARB use by 1% was associated with a 0.12 % increase in Covid-19 confirmed cases. The use of ARB, due to known ACE2 upregulation, may facilitate SARS-CoV-2 entry into target cells and increase infectivity. Cluster-randomized controlled trial is warranted to answer the question of whether the replacement of ARB by ACEI may reduce the Covid-19 confirmed case rate.
Funding LGT was supported in part by the National Institute of Health (HL118277).
Introduction
The novel coronavirus disease 2019 (Covid-19) caused by SARS-Cov2 virus was named a pandemic officially by the World Health Organization on March 11th, 2020.1 It has been shown that SARS-Cov2 uses the angiotensin-converting enzyme 2 (ACE2) receptor as an entry point into a cell.2
With the ACE2 receptor acting as a binding site for SARS-Cov2, the Renin-Angiotensin- Aldosterone System (RAAS), and the medications affecting it become important points of discussion.3 Angiotensin-converting enzyme inhibitors (ACEI) and angiotensin receptor blockers (ARB) are two classes of medications widely used in patients with hypertension, diabetes, cardiovascular disease (CVD), and latent or manifest left ventricular dysfunction. Hypertension, diabetes, and CVD emerged as risk factors for severe Covid-19 cases and deaths.4 Previous experiments showed that ACE2 expression is associated with susceptibility to SARS-Cov infection.5
Notably, available clinical studies consider ACEI and ARB together. However, the effects of ACEI and ARB on ACE2 levels and activity differ.3 While there is strong evidence that ARBs increase ACE2 expression6-11 and augment ACE2 activity6,12, only Ferrario et al13 showed that ACEIs increase ACE2 expression, whereas other studies showed no effect of ACEI on ACE2 activity.14-16
The consensus is reached by all international cardiac societies about the importance of the continuation of ACEI and ARB use in Covid-19 pandemic. However, it remains unknown whether clinically indicated use of ACEI and ARB improve or worsen infectivity or the course of Covid-19. In this rapidly growing pandemic, time is of an essence. To address an urgent need, we conducted a geospatial study. We hypothesized that in the geospatial analysis, the rate of ACEI and ARB use are associated with the number of confirmed Covid-19 cases and deaths in the United States (US).
Methods
We conducted a geospatial disease mapping study using publicly available county-level data. The study was reviewed by the Oregon Health & Science University Institutional Review Board and assigned a determination of Not Human Research. We provided the study dataset and STATA (StataCorp, College Station, Texas) code at https://github.com/Tereshchenkolab/geospatial, allowing future replication and update of the study results as COVID-19 pandemic is unfolding.
Geographical framework
An individual county in the US was an observation unit in this study. We used a Federal Information Processing Standard (FIPS) county code to link the data. Data with missing FIPS codes were excluded from the study. Geographic information about each county was obtained from the cartographic boundary files (shapefiles) provided by the US Census Bureau’s MAF/TIGER geographic database.17
Exposure: rate of ACEI and ARB use by the Medicare Part D beneficiaries
We used the 2017 Centers for Medicare & Medicaid Services (CMS) public dataset, the Medicare Provider Utilization and Payment Data: Part D Prescriber Public Use File, with information on prescription drugs prescribed by individual physicians and other health care providers and paid for under the Medicare Part D Prescription Drug Program in 2017.18 The dataset included the total number of prescriptions that were dispensed (total day supply), which include original prescriptions and any refills, and, therefore, reflects ACEI and ARB usage.
The Medicare dataset only includes city and state information and not county information. In order to map the Medicare prescription data to their corresponding FIPS Code, we used the Google Geocoding API.19 We then loaded the Medicare data and geocoded data into the SQLite database to produce the final datasets with prescription counts per county. Prescriptions from county-equivalents (independent cities) were manually matched with their corresponding FIPS code. Medicare prescriptions with misspelled cities or prescriptions that lacked valid city and state descriptions, if unable to be assigned, were excluded. Excluded prescriptions accounted for <0.01% of the data.
We calculated a drug class use rate as a sum of total day supply in a county for all drugs comprising a particular class (Table 1), normalized by the total county population estimate. We used the US Census Annual Resident Population Estimates for July 1st, 2019.
Outcomes: COVID-19 confirmed cases and deaths
The primary outcome was confirmed Covid-19 cases. The secondary outcome was Covid-19 deaths. We imported the raw COVID-19 data from the Johns Hopkins GitHub repository.20 The number of confirmed COVID-19 cases and deaths in each county as reported for June 11th, 2020 was divided by the total population in each county (2019 county population estimate) and multiplied by 100,000 to convert to cases and deaths per 100,000 population.
Covariates: Population characteristics
Demographic and socioeconomic characteristics
We used the US Census County Population Estimates, released in March 2020, and included reported deaths and births in period July 1st, 2018 to June 30th, 2019.21 Due to known negative impact of Covid-19 on the population of nursing homes and prisons/jails, we included July 1st, 2019 Group Quarters total population estimate. Group Quarters Facilities include correctional facilities for adults, nursing homes, college/university student housing, military quarters, and group homes. Group Quarters data was gathered from an estimated 20,000 randomly selected facilities. Data was then collected through resident interviews of these selected facilities using the American Community Survey conducted by the US Census Bureau.21 The total 2019 county population estimate21 normalized all demographic characteristics.
To characterize socioeconomic characteristics, we used the 2018 median household income expressed as a percent of the state total, and percent of the total population in poverty, as reported by the Economic Research Service of the US Department of Agriculture.22 We also used the data compiled by the County Health Rankings & Roadmaps program, which is a collaboration between the Robert Wood Johnson Foundation and the University of Wisconsin Population Health Institute.23
Cardiovascular disease prevalence and severity
To characterize CVD prevalence and severity, we used the Centers for Disease Control and Prevention (CDC) estimates24 of total CVD death rate per 100,000 population (2016-18), total CVD hospitalizations (2015-17) per 1,000 Medicare beneficiaries, heart failure (HF) death rate per 100,000 population (2016-18), HF hospitalization rate per 1,000 Medicare beneficiaries (2015-17), coronary heart disease (CHD) death rate per 100,000 population (2016-18), CHD hospitalization rate per 1,000 Medicare beneficiaries (2015-17), and age-adjusted diabetes percentage in adults (age > 20y). This data was obtained from the Interactive Atlas of Heart Disease and Stroke, published by the CDC. Within this atlas, death rates were gathered from the Deaths National Vital Statistics Program, hospitalization rates were gathered from the Centers for Medicare and Medicaid Services Medicare Provider Analysis and Review file Part A, and diabetes percentages were collected from the Division of Diabetes Translation.24
To characterize the use of cardiovascular medications, we calculated the rate of cardiovascular medications use, which included original prescriptions and any refills (total day supply), as reported in the 2017 CMS Part D Medicare Prescriber Public Use File.18 We considered the total day supply data for 20 medication groups (Table 1): ACEI, ARB, beta- blockers, alpha-and-beta-blockers, alpha-blockers, class I, III, and V antiarrhythmic medications, dihydropyridine, and non-dihydropyridine calcium channel blockers, aldosterone antagonists, central acting antihypertensive medications, vasodilators, diuretics, lipid-lowering drugs, insulins, and oral hypoglycemic agents, anticoagulants and antiplatelet medications. We normalized the cardiovascular medications day supply for each county by the 2019 county population estimate.21
Statistical analyses
Data are summarized as mean ± standard deviation or as the median and interquartile range (IQR) if non-normally distributed. A paired t-test was used to compare an average rate of ACEI and ARB use.
As Covid-19 is a contagious disease, confirmed cases and deaths in neighboring counties are spatially correlated. Therefore, we used spatial autoregression model25 that allows modeling the spatial dependence among the outcomes, covariates, and among unobserved errors.26 The spatial autoregression model used the generalized spatial two-stage (method-of-moment), least-squares estimator.27 The model included spatial lags of the outcome variable, spatial lags of covariates, and spatially autoregressive errors. The lag operator was a spatial weighting (inverse-distance) matrix, which summarized spatial relationships between counties, based on the distance between county centroids. Weighting matrix was scaled so that its largest eigenvalue is 1, which guarantees nonsingularity in the model estimation. We constructed cross-sectional spatial autoregressive models. The estimator treated the errors as heteroskedastic, thus relaxing the assumption that errors represent identically distributed disturbance.
We conducted the Moran test to determine whether exposure, outcome, and covariate variables are spatially dependent.
First, we constructed unadjusted spatial autoregression models, to investigate a geospatial association of the county population characteristics with the ACEI and ARB use rate, calculated as follows:
ACEI use rate = ln (1+total ACEI day supply/county population).
ARB use rate = ln (1+total ARB day supply/county population).
Each model included spatial lags of the outcome variable (ACEI or ARB use rate, one-by- one), spatial lags of the test’s predictor variable (demographic, socioeconomic, and CVD prevalence characteristics, one-by-one), and spatially autoregressive errors (reflecting unobserved factors).
Next, we constructed unadjusted and bivariate spatial autoregression models, to evaluate a geospatial association of the county population characteristics with the rate of Covid-19 confirmed cases and deaths. To normalize the distribution of the outcome variables, and to improve the interpretability of models, we transformed outcome variables as follows:
Covid-19 confirmed case rate = ln (1+ confirmed Covid-19 cases/100,000 population); Covid-19 death rate = ln (1+ confirmed Covid-19 deaths/100,000 population).
Each model included spatial lags of the outcome variable (Covid-19 confirmed case rate and death rate, one-by-one), spatial lags of the tested predictor variable (demographic and socioeconomic characteristics, cardiovascular risk factors and CVD prevalence, and the rate of
ACEI or ARB use, one-by-one), and spatially autoregressive errors (unobserved impacts). The bivariate model included both predictors (ACEI and ARB use) together.
Finally, we constructed adjusted spatial autoregression models to answer the main study question: whether there is an independent association of ACEI and ARB use rate with Covid-19 confirmed case and death rate. The selection of covariates for adjustment was guided by confounders observed in this study and model fit. We adjusted for demographic and socioeconomic characteristics of a county population (percentage of Black non-Hispanic county residents, percentage of a county residents younger than 18 years of age, percentage of residents with at least some college degree, median household income as a percent of state total), air quality index, as well as CDC-reported CVD hospitalization rate in Medicare beneficiaries, and CVD death rate in a total county population. Each model included spatial lags of the outcome variable (Covid-19 confirmed cases and deaths, one-by-one), spatial lags of the tested predictor variable (rate of ACEIs or ARBs usage, one-by-one), spatial lags of the selected (listed above) 7 covariates (altogether) and spatially autoregressive errors (unobserved influences).
Sensitivity analyses
Cross-sectional geospatial analysis is susceptible to reverse causality bias. It is well- documented that patients with CVD and cardiovascular risk factors (hypertension, diabetes) have a higher rate of Covid-19 confirmed cases and deaths. The rate of ACEIs and ARBs use indirectly indicates CVD prevalence and severity. While we adjusted our models for the broad range of confounders, including CVD mortality in a total county population, and CVD hospitalization rate among Medicare beneficiaries, reverse causality remained of concern. To assess the possibility and extent of reverse causality bias, we constructed the aforementioned spatial autoregression models for the use rate of other cardiovascular medications, for each class of drugs separately, one-by-one.
As it is known that the angiotensin receptor–neprilysin inhibitors are used as second-line therapy for the treatment of advanced HF, we conducted a sensitivity analysis with pure ARB class of medications to minimize the reverse causality bias.
Results
Rate of ACEI and ARB use, and their association with population characteristics
We analyzed the data of 3,141 counties and county-equivalents in the 50 states and the District of Columbia. ACEI were the second most ubiquitous medications, surpassed only by lipid-lowering drugs (Table 1). Average county characteristics are reported in Table 2. Figure 1 shows the ACEI and ARB total day supply rates across the US. On average, the total day supply rate was higher for ACEI than for ARB (10.3±6.5 vs. 6.0±4.6 days; P<0.0001), as shown in Figure 1C. The Moran test indicated that the rates of ACEI and ARB use were spatially dependent (P<0.0001).
In unadjusted spatial autoregression analysis (Table 3), as expected, CVD prevalence, general demographic characteristics, uninsured rate, and air quality were associated with the use of both ACEI and ARB. A higher percentage of adults above 65 y of age was associated with higher use of ACEI, but not ARB. A higher percentage of Asians was associated with the use of ARB, but not ACEI.
Covid-19 confirmed case and death rate
Covid-19 confirmed case and death rate (Figure 2) had similar geographic distributions. As of June 11th, 2020, in an average county, there were medians of 168.3 (IQR 67.6 – 410.9) confirmed cases and 2.95 (IQR 0 – 13.7) deaths per 100,000 population.
In unadjusted spatial autoregression analysis (Table 4), higher CVD and HF hospitalization rates among Medicare beneficiaries, higher percentage of Blacks, Hispanics, Asians, and Pacific Islanders, greater proportions of children (age < 18y), higher percentage of residents admitting excessive drinking, higher percentage of households with high housing costs, and worse air quality were associated with a higher rate of confirmed Covid-19 cases. In contrast, a larger proportion of county residents admitting physical inactivity, with some college education, greater percentage of adults above 65 years of age and non-Hispanic Whites, and higher CVD, HF, and total mortality across all ages were associated with a lower rate of confirmed Covid-19 cases (Table 4). The same factors that affected Covid-19 confirmed case rate also affected Covid-19 death rate, but the strength of the impact on deaths was lesser. As expected, we observed significant indirect (spillover) effects of socioeconomic factors coming from neighboring counties on confirmed Covid-19 case and death rate to a given county.
Association of ACEI and ARB use rate with confirmed Covid-19 case and death rate
In unadjusted analysis, as expected, both ACEI and ARB use rate had a direct and indirect impact on confirmed Covid-19 case and death rate (Table 5). In bivariate analysis, only ARB but not ACEI use was associated with higher rate of confirmed cases. The association of ACEI use with Covid-19 was fully explained by confounders (Table 5). In adjusted analysis, the ACEI use rate had no direct or indirect impact on Covid-19 confirmed case rate, whereas the ARB use was associated with higher rate of confirmed Covid-19 cases. On average, an increase in ARB use rate in a given county by 1% was associated with an increase in Covid-19 confirmed case rate by 0.12%.
In unadjusted analysis, both ACEI and ARB use rate was associated with Covid-19 death rate. In the bivariate model, ACEI and ARB use rates demonstrated the opposite direct (county- own) impact on the death rate. Higher ACEI use rate was associated with a lower death rate, whereas higher ARB use rate was associated with higher death rate (Table 5). After adjustment for demographic and socioeconomic characteristics, air quality index, and CVD prevalence and mortality, the rates of ACEI and ARB use remained associated with Covid-19 death rate. On average, an increase in ARB use rate in a given county by 1% was associated with a 0.06% increase in Covid-19 deaths in that county (direct, county-own effect).
Of note, adjustment fully explained the indirect (spillover) impact of ACEI and ARB use rate coming from all neighboring counties on Covid-19 confirmed case and death rate to a given county, while all models confirmed strong spatial dependence (Moran test of spatial terms for all models P<0.00001). The marginal analysis showed no significant differences in outcomes if ACEI and ARB use rate would change in the same direction (either increase or decrease) in all counties (Supplemental Figure 1).
Sensitivity analysis
In unadjusted analysis, the use of cardiovascular medications reflects CVD prevalence. As expected, the use of all types of cardiovascular medications was associated with Covid-19 confirmed case rate (Supplemental Table 1). Overall, the association with the Covid-19 death rate was similar but had a smaller effect size. The effect size for ACEI and ARB use rate was comparable to that of other cardiovascular medications.
In adjusted analysis (Supplemental Table 2), the rate of use of the vast majority of cardiovascular medications had no association with Covid-19 confirmed case rate, with very few exceptions. This finding suggests an absence of significant reverse causality bias and allows meaningful interpretation of the main study results for the primary outcome.
Similarly to the main results, after adjustment, most of the cardiovascular mediations remained significantly associated with Covid-19 death rate (Supplemental Table 2). Therefore, reverse causality bias was likely present in our analyses with the Covid-19 death rate outcome, which limited meaningful interpretation of the main results for the secondary outcome.
Results of the analyses with pure ARB class of medications did not differ from the main study results (Table 5).
Discussion
Our study confirmed that the rate of ACEI use does not impact Covid-19 confirmed case rate. However, we observed that ARB use associated with slightly increased Covid-19 confirmed case rates. After adjustment for demographic and socioeconomic confounders and CVD prevalence and severity, an increase in ARB use by 1% was associated with a 0.12 % increase in Covid-19 confirmed cases. These findings suggest that long-term use of ARB, due to known ACE2 upregulation, may facilitate SARS-CoV-2 entry into target cells and increase infectivity. ACEI and ARB are frequently prescribed interchangeably for the same clinical indications. Cluster- randomized controlled trial is warranted to answer the question of whether the replacement of ARB by ACEI may reduce the Covid-19 confirmed case rate. Our results highlight the safety and indicate possible benefits of ACEI use for patients with clinical indications for ACEI in the Covid-19 era, consistently with several other studies of ACEI in Covid-19.28,29
Importantly, in this observational geospatial study, residual confounding and reverse causality bias cannot be completely ruled out. The use of ARB may indicate a subgroup of CVD patients who are especially vulnerable to the virus. In such a case, the use of ARBs is not a cause, but a marker of risk. Reverse causality bias likely explains an association of ACEI and ARB use with the Covid-19 death rate, which was similar to the vast majority of other cardiovascular medications.
SARS-CoV-2 virus may preferentially infect individuals taking ARB, but not ACEI
There is strong evidence that the entry of the SARS-CoV-2 virus into the host cell depends on the SARS-CoV receptor ACE2. ACE2 is a type I integral membrane glycoprotein expressed mainly in the respiratory tract, heart, kidneys, and gastrointestinal tract. ACE2 tissue expression facilitates the virus entry in target cells.5
Clinical indications for ACEI and ARB are similar, and in many clinical studies, these two classes of drugs considered together under the common name “RAAS inhibitors.” However, the effects of ACEI and ARB on ACE2 levels and activity are different.3 Our findings support the hypothesis that long-term use of ARB, but not ACEI may facilitate SARS-CoV-2 entry and increase infectivity.
Several unmeasured confounders could be responsible for our findings. First, ARB are indicated for patients who are intolerant to ACEI, usually because of bradykinin-mediated cough. Second, both ACEI and ARB can be used together in patients with advanced HF, and thus indicate a high-risk patient population. For those patient categories, switching from ARB to ACEI is not an option. Nevertheless, regardless of whether ARB indeed increases infectivity or simply indicates a high-risk patient population, it would be wise to reinforce effective Covid-19 prevention strategies, to minimize the risks of exposure to the virus.
Our findings support the notion of ACEI and ARB playing a double-edged sword role. The strength of associations of ACEI and ARB use with the Covid-19 death rate was twice less than with confirmed case rate, suggesting no adverse effect on Covid-19 disease severity.
Limitations
Although the Medicare Part D Prescriber Public Use File has a wealth of information, the dataset has several limitations. The data may not be representative of a physician’s entire practice or all of Medicare as it only includes information on beneficiaries enrolled in the Medicare Part D prescription drug program (approximately two-thirds of all Medicare beneficiaries). Besides, available data were for the year 2017 and did not reflect the most recent use of medications in 2020. Nevertheless, we measured exposure before the outcome, which is essential for the interpretation of the study findings. Furthermore, we did not adjust for adherence to medications. Nevertheless, a recent geospatial study of ACEI/ARB adherence30 showed a relatively consistent geographic distribution of ACEI/ARB adherence across the US. An observational cross-sectional geospatial study is susceptible to reverse causality bias. To address this limitation, we performed a rigorous analysis of all other classes of cardiovascular medications. The results of sensitivity analyses helped to identify the reverse causality bias in analyses with the secondary outcome. Finally, unobserved confounding was likely present in this observational study. The most apparent missing data included the rate of Covid-19 testing. Therefore, observed effect sizes have to be interpreted with caution. However, unobserved confounding would unlikely to affect a relative comparison of ACEI and ARB.
Data Availability
We provide the study dataset and STATA (StataCorp, College Station, Texas) code at https://github.com/Tereshchenkolab/geospatial.
Funding Sources
LGT was supported in part by the National Institute of Health (HL118277).
Disclosures
None.
Acknowledgments
Authors acknowledge numerous scientists and institutions that provided open data.