ABSTRACT
Objectives To evaluate antihypertensive medications and COVID-19 diagnosis and mortality, accounting for healthcare seeking behaviour.
Design A population-based case control study with additional cohort analysis.
Setting Primary care patients from the UK Clinical Practice Research Datalink (CPRD).
Participants 16 866 patients with COVID-19 events in the CPRD from 29th January to June 25th 2020 and 70 137 matched controls.
Main outcome measures We explored associations between COVID-19 diagnosis and prescriptions for angiotensin converting enzyme inhibitors (ACEIs), angiotensin receptor blockers (ARBs), beta-blockers (B), calcium-channel blockers (C), thiazide diuretics (D) and other antihypertensive drugs (O). We evaluated all-cause mortality among COVID-19 cases. Analyses were adjusted for covariates and consultation frequency.
Results In covariate adjusted analyses, ACEIs were associated with lower odds of COVID-19 diagnosis (0.82, 95% confidence interval 0.77 to 0.88) as were ARBs, 0.87 (0.80 to 0.95) with little attenuation from adjustment for consultation frequency. In fully adjusted analyses, C and D were also associated with lower odds of COVID-19. Increased odds of COVID-19 for B (1.19, 1.12 to 1.26), were attenuated after adjustment for consultation frequency (1.01, 0.95 to 1.08). In adjusted analyses, patients treated with ACEIs or ARBs had similar mortality to patients treated with classes B, C, D or O (1.00, 0.83 to 1.20) or patients receiving no antihypertensive therapy (0.99, 0.83 to 1.18).
Conclusions Associations were sensitive to adjustment for confounding and healthcare seeking, but there was no evidence that antihypertensive therapy is associated with increased risk of COVID-19 diagnosis or mortality; most classes of antihypertensive therapy showed negative associations with COVID-19 diagnosis.
INTRODUCTION
Coronavirus disease 2019 (COVID-19) caused by the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has become a global pandemic. The World Health Organization’s 29 June 2020 situation report documented 10 million cases and nearly 500 000 deaths of COVID-19 globally.1 The United Kingdom is among the countries with the highest COVID-19 incidence and death rates in the world, with figures from 14 September 2020 indicating that among those first testing positive for the virus, in excess of 41 637 died within 28 days.2
Biomedical scientific evidence suggests a role of the renin-angiotensin-aldosterone system (RAAS) in COVID-19. SARS-CoV-2 enters host cells via interaction with the angiotensin-converting enzyme 2 (ACE2) receptor which is part of the renin-angiotensin-aldosterone system (RAAS).3-6 Angiotensin-converting enzyme inhibitors (ACEIs) and angiotensin receptor blockers (ARBs) modulate the RAAS system and treatment with ACEIs and ARBs may enhance ACE2 activity thereby increasing SARS-CoV-2 susceptibility and COVID-19 severity.7,8 Conversely, increased ACE2 might have a protective effect by competitive inhibition of SARS-CoV-2 entry into respiratory epithelium or via negative regulation of the RAAS for anti-inflammatory, antioxidative and vasodilatory effects.7,9
Antihypertensive drugs represent the most frequently prescribed medicines in the UK, used by 15% of adults with long term conditions including diabetes, hypertension and cardiovascular diseases.10 These pre-existing conditions are frequent among those receiving healthcare for COVID-196,11 and have also been associated with high COVID-19 case fatality rates,12-14 raising concerns about possible associations with AHT. There have been conflicting results from several observational studies exploring the relationship between AHT and COVID-19 susceptibility and severity since the start of the pandemic.15,16
A living systematic review, updated to 3rd August 2020, included three studies with 8,766 COVID-19 patients and found no evidence of association between either ACEI or ARB treatment and positive COVID-19 test results.15 Two of the studies drew on patients attending hospital in the US,17,18 while the third study was population-based, including both primary care and hospital attendances in the Lombardy region in Italy.19 A large international study which was initially included in the review has since been retracted by the journal in which it was published. A retrospective cohort study using Danish national disease registries which found no evidence of association of prior treatment with ACEIs/ARBs with COVID-19 diagnosis, severity or death.20 In contrast, a systematic review and meta-analysis of 28,872 patients found that treatment with RAAS inhibitors (RAASi) was associated with lower risk of death or critical events.16 A primary care database study in the UK reported that prior use of ACEIs and ARBs were associated with lower risk of COVID-19 diagnosis, but there was no association with intensive care unit admission.21
Non-randomised studies to evaluate the association of RAASi drugs on COVID-19 susceptibility and severity face several methodological challenges. Firstly, there may be substantial confounding. Several previous studies have only presented ‘fully adjusted’ estimates, making it difficult to assess the possibility of residual confounding. Secondly, analysis of clinical data from electronic health records may be associated with misclassification of confounders, which might bias estimates in either positive or negative directions. Thirdly, and perhaps most importantly, cases of COVID-19 documented into electronic health records represent only a minority of cases occurring in the community during the first wave of the pandemic. If hypertension and COVID-19 are both associated with increased health care utilisation, this might introduce spurious associations between antihypertensive therapy and COVID-19 through collider bias.22,23 Selection pressures may have biased the sample toward those with greater symptom severity24 or who have comorbidities that make contact with health services more likely.
This study aimed to add to the evidence concerning antihypertensive therapy and COVID-19 diagnosis and mortality through analysis of a large population-based sample registered in UK primary care. We aimed to evaluate the association between AHT classes including ACEIs, ARBs, calcium-channel blockers, beta-blockers and thiazide diuretics and COVID-19 diagnosis and mortality. We aimed to explore the effects of covariate adjustment including patients’ consultation frequency in primary care.
METHODS
Overview
We undertook a case-control analysis of CPRD GOLD to estimate associations between antihypertensive treatments and COVID-19 susceptibility. We examined associations between antihypertensive treatments and COVID-19 mortality using date of death in the CPRD within 30 days as the outcome. Analyses were adjusted for covariates and consultation frequency.
Case and control selection
The CPRD GOLD is among the world’s largest databases of anonymised electronic health records from primary care. It is estimated that 7% of UK general practices contribute to CPRD GOLD enabling the database to have extensive coverage and good sociodemographic and geographic representativeness of the UK population.25 The high quality of the CPRD GOLD data has been verified by several studies.26 The protocol was approved by the CPRD Independent Scientific Advisory Committee (ISAC protocol 20_081RA).
Patients recorded with COVID-19 were identified from the July 2020 release of CPRD GOLD. COVID-19 events were identified from Read codes recorded into patients’ clinical, referral and test records: ‘Suspected disease caused by 2019-nCoV (novel coronavirus)’, 39%; ‘Suspected coronavirus infection’, 19%; ‘Telephone consultation for suspected 2019-nCoV (novel coronavirus)’, 16%; ‘2019-nCoV (novel coronavirus) detected’, 12%; ‘Disease caused by 2019-nCoV (novel coronavirus)’, 7%; ‘Coronavirus infection’, 5%; ‘Coronavirus nucleic acid detection’, 3%; ‘[X] Coronavirus infection, unspecified’ and ‘Coronavirus as cause of disease classified to other chapters’, less than 1%. We excluded the small number of COVID-19 diagnosis and mortality events dated on or before 29th January 2020, the official date of the UK’s first confirmed COVID-19 case.
COVID-19 cases were compared with a sample of matched control patients who did not have a COVID-19 diagnosis. Up to five control patients for each case were randomly sampled from the entire registered population of CPRD GOLD, matching on age, gender, index date and general practice.
Exposures
The exposure was defined as prescription of antihypertensive treatments (AHT) in the six months before the index date of the following classes: ACEIs, ARBs, beta blockers (B), calcium channel blockers (C), thiazide diuretics (D) and other antihypertensive drugs (O). Each class of AHT was coded as either prescribed or not prescribed.
Covariates
Covariates were defined using data recorded in the five-year period before date of diagnosis. These included smoking status (non-smoker, current smoker, ex-smoker), body mass index (BMI, underweight, normal weight, overweight, obese, and missing), systolic blood pressure (SBP) and diastolic blood pressure (DBP) in categories of 10 mm Hg. Ethnicity was classified as ‘white’, ‘black’, ‘Asian’, ‘mixed’, ‘other’ and ‘not known’. We used Clegg’s electronic Frailty Index (eFI) to evaluate frailty into categories of non-frail, mild, moderate and severe frailty 27. We also evaluated the 15 comorbidities of the Charlson comorbidity index as present or absent, following the recommendations of Khan et al.28,29 We evaluated healthcare seeking behaviour by calculating the rate of events in each patient’s clinical record file, including consultations and other contacts with the practice, in the year preceding the index date. This rate per patient year was entered as a continuous variable.
Analysis
A conditional logistic regression model was employed for the case control analysis. Unadjusted odds ratios were compared with covariate adjusted odds ratios including ethnicity, BMI, blood pressure, smoking status, frailty level, comorbidities and treatment with each class of AHT. Finally, additional adjustment was made for consultation frequency. As a sensitivity analysis we repeated the analysis only including patients with confirmed COVID-19 diagnoses, excluding patients recorded to have ‘suspected’ COVID-19. We also explored evidence of an interaction between age group and each class of AHT, conducting a sub-group analysis to examine age as an effect modifier.
A Cox proportional hazards model was employed to examine the association of AHT with patient survival. Survival time in days from the date of first COVID-19 diagnosis to date of death or end of record was the outcome. Patients who were still alive at the end of the study period were censored. Covariates were month of COVID-19 diagnosis, gender, age (zero to four, five to nine and 10 to 14 and then 10 years age groups up to 85 years and over), and smoking status, BMI, blood pressure, ethnic group, eFI category comorbidities and treatment with each class of AHT. A secondary analysis also adjusted for region of practice. We did not adjust for consultation frequency because all patients in this analysis had already accessed care from their general practice. We also evaluated ACEIs and ARBs combined, as renin-angiotensin-aldosterone system inhibiting drugs (RAASi), making comparisons with other classes of AHT drugs (BCDO) or no antihypertensive treatment. The proportional hazards assumption was checked graphically and with a test for proportional hazards. All data were analysed in R, version 4.0.2.
RESULTS
Case-control analysis: COVID-19 diagnosis
There were 16 866 COVID-19 cases identified in CPRD GOLD from 29 January to 20 June 2020 (Supplementary Figure 1), which were compared to 70 137 matched controls. COVID-19 cases had consulted their general practice more frequently over the preceding year than controls, with a mean rate of clinical record events in the preceding year of 28.0 per person year, compared to 15.4 per person year among controls.
Table 1 and Supplementary Table 1 show descriptive data for COVID-19 cases and controls as well as unadjusted and covariate adjusted odds ratios and 95% confidence intervals.
COVID-19 diagnosis was associated with current smoking, underweight and obese BMI, frailty and black and minority ethnicity. Both low and high systolic or diastolic blood pressure values were associated with increased odds of COVID-19 diagnosis when compared with intermediate blood pressure values. Comorbidities associated with COVID-19 diagnosis in covariate adjusted analyses included chronic pulmonary disease, metastatic tumour, congestive heart disease, cerebrovascular disease and dementia. Diabetes and renal disease were associated with COVID-19 diagnosis in unadjusted but not covariate adjusted analyses.
There were 4 852 (29%) COVID-19 cases prescribed antihypertensive treatment (AHT) in the year preceding diagnosis, compared to 17 978 (26%) of controls. Positive associations of ACEI and ARB with COVID-19 diagnosis in unadjusted analyses (Supplementary Table 1) were attenuated in the covariate adjusted analysis (Table 1), with both ACEIs and ARBs being negatively associated with COVID-19 diagnosis (ACEI: odds ratio 0.82, 95% CI 0.77 to 0.88; ARB: 0.87, 0.80 to 0.95). Further adjustment for consultation frequency had little effect on estimated associations (ACEI: 0.80, 0.75 to 0.86; ARB: 0.84, 0.77 to 0.92) (Figure 1). In the fully adjusted model, calcium channel blockers and thiazide diuretics were also associated with lower odds of COVID-19 diagnosis (Figure 1). In covariate adjusted analyses, beta-blocker treatment was positively associated with COVID-19 diagnosis (1.19, 1.12 – 1.26) but this was attenuated after further adjustment for consultation frequency (1.01, 0.95 to 1.08) (Figure 1). When the sample was restricted to only those cases with confirmed COVID-19 diagnoses and their matched controls, there were 4 233 cases matched to 17 700 controls, the pattern of association was unaltered though estimates were necessarily less precise (Supplementary Figure 3).
Adding an interaction term for age within logistic regression models for each class of AHT improved the model goodness of fit, therefore we conducted a subgroup analysis for cases and controls by AHT treatment and age (Supplementary Table 3). Whilst lower age categories had very few observations, AHT therapy was positively associated with higher odds of COVID-19 diagnosis across all classes ABCDO within the mid-range age categories, such as those aged 45-54 years on ACEI (1.47, 1.27 to 1.69). These associations were attenuated after adjusting for all covariates including consultation frequency.
Primary cohort analysis: COVID-19 mortality
Among the 16 866 COVID-19 cases, 921 (5%) died within 30 days. Patients treated with antihypertensive drugs were more highly represented among deceased patients than the overall sample and this was true for each class of AHT drugs (Table 2 and Supplementary Table 2). Table 2 and Supplementary Table 2 provide estimates for covariates which show that male gender, older age, black and ethnic minority status, diabetes, metastatic tumour, dementia and mild liver disease were associated with greater mortality. In covariate adjusted analyses, there was no evidence that any of the classes of AHT drugs might be associated with mortality after COVID-19 diagnosis (Table 2). Additional adjustment for practice region did not influence associations (Figure 2). There was no evidence that treatment with RAASi drugs, including ACEI and ARB drugs, might be associated with higher mortality than other classes of AHT drugs (BCDO) with hazard ratio 1.00 (0.83 to 1.20). Similarly there was no evidence that treatment with RAASi drugs was associated with greater risk than no antihypertensive treatment (0.99, 0.83 to 1.18), nor was treatment with BCD and O classes associated with higher mortality than no AHT treatment (0.99, 0.84 to 1.18).
DISCUSSION
Statement of principle findings
This large population-based study included patients diagnosed with suspected or confirmed COVID-19 in UK primary care. During the first wave of the pandemic, there was limited testing capacity and most suspected cases remained unconfirmed; however, restricting our analysis to confirmed cases yielded similar results. After adjusting for covariates that characterise case-mix, we found no evidence that treatment with ACEIs, ARBs, calcium channel blockers, thiazide diuretics or other antihypertensive drugs might be associated with greater risk of COVID-19 diagnosis. There was evidence that beta-blocker treatment was associated with greater odds of COVID-19 diagnosis but additional adjustment for patients’ underlying consultation frequency removed this association, which might be attributable to collider bias. After adjusting for covariates, including blood pressure, and consultation frequency, ACEIs, ARBs, calcium channel blockers and thiazide diuretics were associated with lower odds of COVID-19 diagnosis. This might suggest that at a given level of blood pressure, patients treated with antihypertensive drugs may be at lower risk, but this pattern of association does not appear to be specific to any class of antihypertensive drug.
Associations of AHT with COVID-19 diagnosis may be modified by age, but the effects of adjusting for covariates, including comorbidities, suggest this attributable to confounding. In this cohort analysis of patients diagnosed in UK primary care with suspected or confirmed COVID-19, recorded mortality was considerably higher than the overall reported infection fatality ratio for this condition but there was no evidence that any class of antihypertensive drug might be associated with increased mortality risk, and the adjusted hazard for drugs acting on the renin-angiotensin-aldosterone system was similar to that for patients treated with other classes of AHT drugs or patients not receiving AHT therapy.
Strengths and weaknesses of the study
This study drew on a large, longitudinal population-based data resource that enabled us to conduct a matched case-control analysis of risk factors for clinical COVID-19 diagnosis, as well as a cohort study of risk factors for COVID-19 mortality. As noted above, most cases were diagnosed clinically as suspected cases because of the limited capacity for testing in the early stage of the pandemic. General practice systems capture comprehensive data for all prescriptions issued by the general practice, so we can be confident that this exposure was accurately recorded. However, prescription utilisation may not be universal and non-adherence to prescribed medicines may be widespread. Data for covariates was not always completely recorded. Data on ethnicity were missing for almost half the sample as reported by others;25 an important limitation given the disproportionate impact COVID-19 has had on ethnic minority populations in the US and UK.14,30,31 COVID-19 susceptibility and severity have also been associated with measures of deprivation.14 Matching on practice may have accounted for differences in area-level deprivation to an extent, but deprivation based on participants’ home postcodes or individual-level deprivation measures might improve precision. It is known that observational studies on COVID-19 outcomes may be susceptible to collider bias.22,23 This can lead to spurious associations if both antihypertensive therapy and COVID-19 are associated with greater likelihood of general practice consultations.
Selection pressures may have biased the sample toward those with increased symptom severity24 or who are otherwise more likely to recognise COVID-19 symptoms and make contact with health services. We accounted for the effects of health-seeking behaviour by introducing into the analysis the rate of events in each patient’s clinical record in the year preceding the index date. There were substantial differences in consultation frequency between cases and controls and adjusting for this metric nullified an apparent association between beta-blocker use and COVID-19 diagnosis. However, randomisation allocation will be preferred to provide a higher level of evidence.
Strengths and weaknesses in relation to other studies
A substantial number of observational studies have previously evaluated whether drugs acting on the renin-angiotensin system, including ACEIs and ARBs, might be associated with susceptibility, severity and mortality from COVID-19, but few studies have considered other classes of antihypertensive medications. Several systematic reviews have summarised the main findings, which generally find no evidence for worse COVID-19 outcomes in patients treated with these medications.32-35 Some studies suggest possible benefits from ACEI treatment.21 However, systematic reviews also highlight limitations of the evidence presented to date.35 A high proportion of studies has been based on data from patients admitted to hospital with COVID-19, sometimes with small samples. Many studies failed to include adequate adjustment for confounding. Few, if any, studies have directly addressed the question of collider bias, which distort associations when data are gathered from patients conditional upon their attendance for health care. Recently, COVID-19 has been evaluated in an ongoing randomised controlled trial of ramipril treatment.36 The analysis found no evidence for an effect of ramipril on COVID-19 incidence or severity but the analysis comprised 102 patients with 11 cases of COVID-19.36 Further trials are ongoing.15
Conclusions
Drawing on data for a large population-based sample and using rigorous analytical methods, this study adds to the evidence that antihypertensive therapy may be safely continued during the SARS-CoV-2 pandemic. While previous studies have largely evaluated drugs acting on the renin-angiotensin system, the study found no evidence that any class of antihypertensive therapy might be associated with greater risk of COVID-19 diagnosis or mortality. There was evidence that, after adjusting for covariates including blood pressure, several classes of antihypertensive therapy might be associated with lower risk of a clinical COVID-19 diagnosis, but this pattern of association was not apparent for COVID-19 mortality. While this might be interpreted as evidence of a protective effect, in this observational study it is not possible to exclude the possibility that this pattern of association may be caused by bias.
Data Availability
Requests for access to data from the study should be addressed to martin.gulliford@kcl.ac.uk. All proposals requesting data access will need to specify planned uses with approval of the study team and CPRD before data release.
Data sources
The study is based in part on data from the Clinical Practice Research Datalink obtained under license from the UK Medicines and Healthcare products Regulatory Agency. However, the interpretation and conclusions contained in this report are those of the authors alone.
Data sharing
Requests for access to data from the study should be addressed to martin.gulliford{at}kcl.ac.uk. All proposals requesting data access will need to specify planned uses with approval of the study team and CPRD before data release.
Funding
MG was supported by the NIHR Biomedical Research Centre at Guy’s and St Thomas’ Hospitals. The views expressed are those of the authors and not necessarily those of the NHS, the NIHR, or the Department of Health. The funder of the study had no role in study design, data collection, data analysis, data interpretation, or writing of the report. The authors had full access to all the data in the study and all authors shared final responsibility for the decision to submit for publication.
Conflict of Interest
Authors declare no conflicts of interest.
Author Contributions
ERP designed the study with advice from MG, AD and PC; ERP analysed the data and drafted the paper; MG, AD and PC reviewed the analysis and contributed to drafts of the paper. All authors approved the final draft. The corresponding author attests that all listed authors meet authorship criteria and that no others meeting the criteria have been omitted.