Europe PMC

This website requires cookies, and the limited processing of your personal data in order to function. By using the site you are agreeing to this as outlined in our privacy notice and cookie policy.

Abstract 


ABSTRACT

Using longitudinal health records from 45.7 million adults in England followed for a year, our study compared the incidence of thrombotic and cardiovascular complications after first, second and booster doses of brands and combinations of COVID-19 vaccines used during the first two years of the UK vaccination program with the incidence before or without the corresponding vaccination. The incidence of common arterial thrombotic events (mainly acute myocardial infarction and ischaemic stroke) was generally lower after each vaccine dose, brand and combination. Similarly, the incidence of common venous thrombotic events, (mainly pulmonary embolism and lower limb deep venous thrombosis) was lower after vaccination. There was a higher incidence of previously reported rare harms after vaccination: vaccine-induced thrombotic thrombocytopenia after first ChAdOx1 vaccination, and myocarditis and pericarditis after first, second and transiently after booster mRNA vaccination (BNT-162b2 and mRNA-1273). These findings support the wide uptake of future COVID-19 vaccination programs.

Free full text 


medRxiv

PPRID: PPR805294
EMSID: EMS194133
medRxiv preprint, version 2, posted 2024 February 13
https://doi.org/10.1101/2024.02.12.24302698

Cohort study of cardiovascular safety of different COVID-19 vaccination doses among 46 million English adults

Affiliations

  1. 1.British Heart Foundation Cardiovascular Epidemiology Unit, Department of Public Health and Primary Care, University of Cambridge, UK
  2. 2.Centre for Cancer Genetic Epidemiology, University of Cambridge, UK
  3. 3.Victor Phillip Dahdaleh Heart and Lung Research Institute, University of Cambridge, Cambridge UK
  4. 4.Department of Population Health Sciences, Bristol Medical School, University of Bristol, UK
  5. 5.MRC Integrative Epidemiology Unit, Bristol, UK
  6. 6.Population Data Science, Swansea University Medical School, Faculty of Medicine, Health, and Life Science, Swansea University, Wales, UK
  7. 7.NIHR Bristol Biomedical Research Centre, UK
  8. 8.Health Data Research UK South-West, Bristol
  9. 9.Health Data Research UK (S.D.), London
  10. 10.Institute of Health Informatics, University College London, UK
  11. 11.University College London Hospitals Biomedical Research Centre, University College London, UK
  12. 12.BHF Accelerator, London, UK
  13. 13.Diabetes Research Centre, University of Leicester, UK
  14. 14.British Heart Foundation Data Science Centre, Health Data Research UK, London, UK
  15. 15.Centre for Clinical Brain Sciences, University of Edinburgh, UK
  16. 16.Department of Public Health and Primary Care, NIHR Blood and Transplant Research Unit in Donor Health and Genomics, University of Cambridge, UK
  17. 17.British Heart Foundation Centre of Research Excellence, University of Cambridge, UK
  18. 18.Wellcome Genome Campus, Health Data Research UK Cambridge
  19. 19.NIHR Cambridge Biomedical Research Centre, UK
  20. 20.Cambridge Centre for AI in Medicine, UK

Copyright and license information

This article is a preprint. It may not have been peer reviewed.
A preprint is a complete scientific manuscript that an author uploads on a public server for free viewing. Initially it is posted without peer review, but may acquire feedback or reviews as a preprint, and may eventually be published in a peer-reviewed journal. The posting of preprints on public servers allows almost immediate dissemination and scientific feedback early in the 'publication' process.

Abstract

Using longitudinal health records from 45.7 million adults in England followed for a year, our study compared the incidence of thrombotic and cardiovascular complications after first, second and booster doses of brands and combinations of COVID-19 vaccines used during the first two years of the UK vaccination program with the incidence before or without the corresponding vaccination. The incidence of common arterial thrombotic events (mainly acute myocardial infarction and ischaemic stroke) was generally lower after each vaccine dose, brand and combination. Similarly, the incidence of common venous thrombotic events, (mainly pulmonary embolism and lower limb deep venous thrombosis) was lower after vaccination. There was a higher incidence of previously reported rare harms after vaccination: vaccine-induced thrombotic thrombocytopenia after first ChAdOx1 vaccination, and myocarditis and pericarditis after first, second and transiently after booster mRNA vaccination (BNT-162b2 and mRNA-1273). These findings support the wide uptake of future COVID-19 vaccination programs.

Introduction

SARS-CoV-2 vaccination prevented 14.4 million deaths from COVID-19 worldwide in the first year of the pandemic 1. In England, which entered its third COVID-19 vaccine season in autumn 2023 2, around 90% of people aged ≥12 years have been vaccinated at least once 3. COVID-19 vaccines are associated with rare cardiovascular complications: mRNA-based brands (e.g. BNT-162b2, mRNA1273) with myocarditis and adenovirus-based brands (e.g. ChAdOx1) with vaccine-induced thrombotic thrombocytopenia (VITT), leading to intracranial venous thrombosis (ICVT) and thrombocytopenia 4,5. It is important to understand the risk of thrombotic and cardiovascular complications arising from second and subsequent doses, within the general population and subpopulations 6.

In this work, we use whole population longitudinal electronic health records from 45.7 million adults in England to quantify associations of first, second, and booster doses of mRNA and non-mRNA COVID-19 vaccine brands with subsequent thrombotic and cardiovascular events. Data were accessed within the NHS England Secure Data Environment (NHSE SDE) 7, encompassing primary care, hospital admissions, COVID-19 testing and vaccination data, dispensed medication records in primary care data, and the Office of National Statistics death registrations. We use Cox regression to estimate adjusted hazard ratios (aHRs) and corresponding 95% confidence intervals (95% CIs) in time intervals since vaccination, adjusted for a wide range of co-morbidities, age, sex, and prior COVID-19. We show that the incidence of thrombotic and cardiovascular complications was generally lower after each dose of each vaccine brand, except for previously recognised rare complications of the ChAdOx1 vaccine and the mRNA vaccines.

Results

Characteristics of the study populations for first, second and booster vaccination

Over the study period from 8th December 2020 to 23rd January 2022, 45.7 million individuals met the eligibility criteria for our first-dose analyses (Table 1, Supplementary Figure 1). Among these, 37.3 million people received a first ChAdOx1, BNT-162b2 or mRNA1273 vaccination and were eligible for the second dose analyses (Supplementary Table 1). We refer to the first and second dose vaccinations as the "primary course", distinguishing them from booster vaccine doses. Third dose vaccination, which is distinct from the booster dose and is administered as part of an extended primary course, was not considered. Of those who received a first dose, 35.9 million people received ChAdOx1, BNT-162b2 or mRNA1273 primary course vaccinations and were eligible for analyses of booster vaccination (Supplementary Table 2).

Table 1. Characteristics (number [%]) of the population under study - dose 1 cohort. CKD: Chronic Kidney Disease, COCP: Combined Oral Contraceptive Pill, COPD: Chronic Obstructive Pulmonary Disease, HRT: Hormone Replacement Therapy
CharacteristicAll eligible for dose 1 analysisVaccination status during dose 1 follow-up
UnvaccinatedChAdOx1BNT162b2mRNA-1273
TotalAll45673965141713151921465511840635447360
SexMale22607555 (49.5%)7714835 (54.4%)9361990 (48.7%)5285235 (44.6%)245495 (54.9%)
Female23066410 (50.5%)6456475 (45.6%)9852665 (51.3%)6555400 (55.4%)201870 (45.1%)
Age (years)18-298847920 (19.4%)6627225 (46.8%)1017765 (5.30%)1169505 (9.88%)33425 (7.47%)
30-398484105 (18.6%)3765585 (26.6%)1621590 (8.44%)2863315 (24.2%)233615 (52.2%)
40-497540600 (16.5%)1798740 (12.7%)4267675 (22.2%)1301565 (11.0%)172620 (38.6%)
50-597700005 (16.9%)1006920 (7.11%)5262530 (27.4%)1424845 (12.0%)5710 (1.28%)
60-695871400 (12.9%)522695 (3.69%)3754280 (19.5%)1592895 (13.4%)1530 (0.34%)
70-794594010 (10.1%)264575 (1.87%)2526135 (13.2%)1802920 (15.2%)380 (0.08%)
80-892157220 (4.72%)136385 (0.96%)574395 (2.99%)1446365 (12.2%)75 (0.02%)
90+478710 (1.05%)49195 (0.35%)190285 (0.99%)239225 (2.02%)5 (0.00%)
Ethnic groupAsian3740455 (8.19%)1491890 (10.5%)1302850 (6.78%)917205 (7.75%)28510 (6.37%)
Black1675070 (3.67%)932490 (6.58%)477935 (2.49%)255595 (2.16%)9050 (2.02%)
Mixed717695 (1.57%)385870 (2.72%)193815 (1.01%)130360 (1.10%)7650 (1.71%)
Other1448035 (3.17%)839310 (5.92%)356930 (1.86%)236285 (2.00%)15510 (3.47%)
White36479155 (79.9%)9565930 (67.5%)16457960 (85.7%)10082615 (85.2%)372650 (83.3%)
Unknown1613555 (3.53%)955820 (6.74%)425170 (2.21%)218580 (1.85%)13985 (3.13%)
Index of multiple deprivation (1 most, 10 least deprived)1-28651540 (18.9%)3461120 (24.4%)3224430 (16.8%)1899160 (16.0%)66830 (14.9%)
3-49015905 (19.7%)3254940 (23.0%)3511375 (18.3%)2167490 (18.3%)82100 (18.4%)
5-68815855 (19.3%)2585605 (18.2%)3786185 (19.7%)2360150 (19.9%)83915 (18.8%)
7-88580765 (18.8%)2126380 (15.0%)3935320 (20.5%)2434215 (20.6%)84850 (19.0%)
9-108300280 (18.2%)1756675 (12.4%)4026770 (21.0%)2422325 (20.5%)94510 (21.1%)
Missing2309620 (5.06%)986590 (6.96%)730570 (3.80%)557300 (4.71%)35160 (7.86%)
Smoking statusCurrent7788670 (17.1%)3238785 (22.9%)2970250 (15.5%)1506115 (12.7%)73520 (16.4%)
Former10414755 (22.8%)1881795 (13.3%)5236765 (27.2%)3202600 (27.1%)93595 (20.9%)
Never25332265 (55.5%)7480385 (52.8%)10706425 (55.7%)6881715 (58.1%)263740 (59.0%)
Missing2138275 (4.68%)1570345 (11.1%)301215 (1.57%)250210 (2.11%)16505 (3.69%)
Medical historyAMI1161760 (2.54%)95090 (0.67%)565960 (2.95%)500385 (4.23%)325 (0.07%)
Diabetes3907740 (8.56%)371355 (2.62%)1992030 (10.4%)1541315 (13.0%)3040 (0.68%)
Depression9426550 (20.6%)2196140 (15.5%)4576145 (23.8%)2578470 (21.8%)75795 (16.9%)
Obesity4384350 (9.60%)636440 (4.49%)2290445 (11.9%)1435515 (12.1%)21950 (4.91%)
Cancer6465720 (14.2%)1025045 (7.23%)3113295 (16.2%)2269930 (19.2%)57450 (12.8%)
COPD1549765 (3.39%)133930 (0.95%)795500 (4.14%)619575 (5.23%)760 (0.17%)
Liver disease225630 (0.49%)39465 (0.28%)116685 (0.61%)69215 (0.58%)265 (0.06%)
CKD2867785 (6.28%)254935 (1.80%)1342385 (6.99%)1268080 (10.7%)2385 (0.53%)
Dementia510600 (1.12%)51920 (0.37%)267380 (1.39%)191240 (1.62%)60 (0.01%)
All stroke835305 (1.83%)75680 (0.53%)421650 (2.19%)337615 (2.85%)360 (0.08%)
All VTE642035 (1.41%)80740 (0.57%)328255 (1.71%)232215 (1.96%)825 (0.18%)
Thrombophilia43440 (0.10%)7665 (0.05%)20745 (0.11%)14765 (0.12%)265 (0.06%)
Major surgery in the last year3438335 (9.23%)3899655 (8.54%)732470 (5.17%)1813645 (9.44%)1330225 (11.2%)
No. of unique diseases in the last year036130080 (79.1%)12473515 (88.0%)14609555 (76.0%)8651765 (73.1%)395245 (88.3%)
1-59412710 (20.6%)1683845 (11.9%)4536390 (23.6%)3140420 (26.5%)52055 (11.6%)
6+131170 (0.29%)13955 (0.10%)68705 (0.36%)48450 (0.41%)60 (0.01%)
RegionNorthwest5658370 (12.4%)1609700 (11.4%)2412950 (12.6%)1581275 (13.3%)54445 (12.2%)
Southeast6574915 (14.4%)1772660 (12.5%)2989610 (15.6%)1735340 (14.7%)77305 (17.3%)
London7210945 (15.8%)3203600 (22.6%)2329760 (12.1%)1605725 (13.6%)71860 (16.1%)
East of England4040560 (8.85%)1118725 (7.89%)1823935 (9.49%)1062575 (8.97%)35325 (7.90%)
Southwest3574255 (7.83%)928270 (6.55%)1610965 (8.38%)984180 (8.31%)50840 (11.4%)
Yorkshire/Humber5096465 (11.2%)1678780 (11.8%)2036080 (10.6%)1318885 (11.1%)62720 (14.0%)
East Midlands4245940 (9.30%)1200815 (8.47%)1893155 (9.85%)1100865 (9.30%)51105 (11.4%)
West Midlands3125510 (6.84%)831225 (5.87%)1427785 (7.43%)862795 (7.29%)3705 (0.83%)
Northeast4462635 (9.77%)1360330 (9.60%)1932690 (10.1%)1142780 (9.65%)26835 (6.00%)
Missing1684370 (3.69%)467225 (3.30%)757715 (3.94%)446210 (3.77%)13220 (2.96%)
Prior COVID-19 at index date2496910 (5.5%)1101915 (2.41%)307405 (2.17%)462060 (2.40%)323250 (2.73%)
Medications taken in the last 3 monthsAntiplatelet2445805 (5.35%)142340 (1.00%)1212305 (6.31%)1090350 (9.21%)810 (0.18%)
Antihypertensive8390090 (18.4%)449315 (3.17%)4509545 (23.5%)3421830 (28.9%)9400 (2.10%)
Lipid lowering6651900 (14.6%)289075 (2.04%)3544825 (18.4%)2814740 (23.8%)3260 (0.73%)
Anticoagulant1300380 (2.85%)70965 (0.50%)612260 (3.19%)616910 (5.21%)245 (0.05%)
COCP612760 (1.34%)283510 (2.00%)142915 (0.74%)177015 (1.49%)9320 (2.08%)
HRT530515 (1.16%)28765 (0.20%)349435 (1.82%)149790 (1.27%)2525 (0.56%)
Clinically vulnerableNo35461085 (77.6%)13080935 (92.3%)14037335 (73.1%)7904020 (66.8%)438795 (98.1%)
Vulnerable8084085 (17.7%)897230 (6.33%)4079130 (21.2%)3100120 (26.2%)7605 (1.70%)
Extremely2128795 (4.66%)193140 (1.36%)1098190 (5.72%)836500 (7.06%)965 (0.22%)

Differences in characteristics of people in the three vaccination cohorts reflect known vaccination uptake profiles. Considering the strong similarity in characteristics of the second and booster dose cohorts, we compared characteristics of the second and booster dose cohorts with the first dose cohort. Compared to the first dose cohort (Table 1), people in the second and booster vaccine dose cohorts (Supplementary Tables 1-2), who had received at least one COVID-19 vaccine dose, were older (<40 years: 30% vs 38%), and less likely to have a non-White ethnicity (16% vs 20%), or be deprived (IMD deciles 1-4, 36% vs 39%). They were modestly more likely to have cancer (18-19% vs 14%), to have contracted COVID-19 previously (7-8% vs 6%), to take blood pressure (22% vs 18%) or lipid-lowering medication (18% vs 15%), or to be clinically vulnerable (25% vs 22%).

Characteristics of people receiving different vaccine brands reflect vaccine availability during the UK vaccine rollout, which was in order of priority groups specified by the Joint Committee on Vaccination and Immunisation (JCVI). People younger than 40 years old and not in a priority group were offered only mRNA (BNT-162b2 or mRNA1273) vaccine brands 8.

Incidence of different cardiovascular events

From the start of vaccine rollout (8th December 2020) up to first vaccination, during approximately 21 million person-years there were 75655 arterial and 21230 venous incident thrombotic events (Table 2 and Supplementary Table 4). Arterial thromboses included acute myocardial infarction (AMI; 37915) and ischaemic stroke (36720). Venous thromboses included pulmonary embolism (PE; 11835) and lower limb deep venous thrombosis (DVT; 9075), intracranial venous thrombosis (ICVT; 370) and portal vein thrombosis (PVT; 185). Other incident cardiovascular events included subarachnoid haemorrhage and haemorrhagic stroke (SAH; 5235), mesenteric thrombus (1515), thrombocytopenia (1885), myocarditis (590) and pericarditis (455). The numbers of events and incidence rates after first, second and booster vaccinations reflected the older age of vaccinated, compared with unvaccinated, people (Supplementary Tables 4 to 7).

Table 2. Number of events, person-years and incidence rates for composite arterial and venous thrombotic events for first, second and booster vaccinations 2.
Vaccination doseBrandEventEvents/100K person-yearsIncidence rate
Dose 1No vaccinationAfter first doseNo vaccinationAfter first dose
ChAdOx1Arterial75655/205.0073330/ 95.63369.05766.83
Venous21230/205.0017525/ 96.09103.56182.39
BNT-162b2Arterial75655/205.0057050/ 80.99369.05704.38
Venous21230/205.0011860/ 81.00103.56146.42
mRNA-1273Arterial75655/205.00255/ 5.11369.0549.89
Venous21230/205.00190/ 5.10103.5637.25
Doses 1 and 2After first doseAfter second doseAfter first doseAfter second dose
ChAdOx1Arterial39425/41.2370485/93.36956.17755.01
Venous8685/41.2716120/93.60210.46172.22
BNT-162b2Arterial26975/33.6355620/72.44802.22767.8
Venous5260/33.6711400/72.69156.22156.83
mRNA-1273Arterial160/2.17200/4.0673.6649.24
Venous90/2.18135/4.0641.3733.29
BoosterPrimary course onlyAfter booster vaccinationPrimary course onlyAfter booster vaccination
Primary course of ChAdOx1 or BNT-162b2 or mRNA-1273BNT-162b2Arterial157235/182.2139735/38.72862.931026.14
Venous32545/182.547855/ 38.97178.29201.54
mRNA-1273Arterial157235/182.223200/ 7.46862.88429.07
Venous32545/182.56760/ 7.47178.27101.79
BNT-162b2 or mRNA-1273Arterial157235/182.2142935/ 46.20862.93929.36
Venous32545/182.588615/ 46.37178.25185.79
Primary course of BNT-162b2 or mRNA-1273BNT-162b2Arterial70215/80.3121425/19.04874.281125.34
Venous13830/80.453975/19.15171.91207.62
mRNA-1273Arterial70215/80.29870/ 2.23874.52389.78
Venous13830/80.45185/ 2.24171.9182.76
BNT-162b2 or mRNA-1273Arterial22300/21.25874.731049.3322300/21.25
Venous4165/21.39171.88194.724165/21.39
Primary course of ChAdOx1BNT-162b2Arterial87020/101.9318310/ 19.76853.73926.59
Venous18715/102.133880/ 19.81183.25195.83
mRNA-1273Arterial87020/101.952325/ 5.21853.59446.61
Venous18715/102.12575/ 5.24183.27109.82
BNT-162b2 or mRNA-1273Arterial87020/101.9520635/ 24.94853.54827.32
Venous18715/102.084455/ 25.08183.33177.61
Primary course of BNT-162b2BNT-162b2Arterial70005/76.6021415/18.78913.881140.44
Venous13700/76.753965/18.91178.51209.65
mRNA-1273Arterial70005/76.59860/ 2.02914426.6
Venous13700/76.80180/ 2.04178.3888.15
BNT-162b2 or mRNA-1273Arterial70005/76.5922275/20.79913.991071.28
Venous13700/76.734145/20.96178.54197.8

COVID-19 vaccination and arterial, venous, and other thrombotic events

We used Cox models to estimate adjusted hazard ratios (aHRs) and corresponding 95% CIs, comparing the incidence of thrombotic and cardiovascular events after first, second and booster vaccine doses with the incidence before or without the corresponding vaccine dose, adjusting for a wide range of potential confounding factors (Supplementary Table 8-25). Associations with second dose were estimated in people who had a first dose, and associations with booster vaccination in people who had both first and second vaccinations. The aHRs after first doses of mRNA-1273 were imprecisely estimated, because this brand was received by fewer people who were generally younger, so were reported in supplementary tables but omitted from the figures.

The incidence of composite arterial thrombotic events (AMI, ischaemic stroke and other arterial embolism) was similar or lower after first, second and booster doses of ChAdOx1 and BNT-162b2 vaccines and booster doses of mRNA-1273, compared to follow-up before or without the corresponding vaccine dose in eligible people (Figure 1, Supplementary Tables 8-25). After second and booster vaccine doses, the reduction in incidence of composite arterial thrombosis was greater than after first vaccination. For example, aHRs for arterial thrombotic events 13-24 weeks after first vaccine dose were 0.99 (95%CI 0.97-1.02) after ChAdOx1 and 0.90 (0.88-0.93) after BNT-162b2. Corresponding aHRs after second doses were 0.73 (0.70-0.76) and 0.80 (0.77-0.83) respectively. aHRs for mRNA booster vaccination after primary course of ChAdOx1 were 0.71 (0.66-0.76) 13-24 weeks after BNT-162b2 and 0.67 (0.62-0.72) 5-24 weeks after mRNA-1273 respectively. After a primary course of BNT-162b2, aHRs 13-24 weeks after booster vaccination were 0.73 (0.69-0.77) and 1.21 (0.38-3.86) after BNT-162b2 and mRNA-1273 respectively. For all vaccine brands and doses, aHRs in the first few weeks were lower than in later weeks. The aHR profiles for AMI and ischaemic stroke were similar to those of composite arterial thrombosis, for all vaccine brands and doses.

Figure 1
Open in new tabFigure 1: Adjusted hazard ratios (aHRs) and 95% confidence intervals (95% CIs) for arterial thromboses following COVID-19 vaccination, by dose and brand. Vertical lines depict 95% CIs; these are not visible when they are very narrow.1

Similar to arterial events, the incidence of composite venous thrombotic events (PE, DVT, ICVT and PVT) was generally lower after first, second and booster doses vaccination, compared to follow-up before or without the corresponding vaccine dose (Figure 2, Supplementary Tables 8-25). After second and booster vaccine doses, the reduction in incidence of composite venous thrombotic events was greater than after first vaccination. For example, aHRs 13-24 weeks after first vaccine doses were 0.94 (95%CI 0.90-0.98) after ChAdOx1 and 0.85 (0.81-0.88) after BNT-162b2. Corresponding aHRs after second doses were 0.68 (0.63-0.73) and 0.77 (0.72-0.83) respectively. aHRs for mRNA booster vaccination after primary course of ChAdOx1 were 0.63 (0.54-0.74) 13-24 weeks after BNT-162b2 and 0.55 (0.47-0.65) 5-24 weeks after mRNA-1273. After a primary course of BNT-162b2, aHRs for booster vaccination were 0.56 (0.49-0.63) 13-24 weeks after BNT-162b2 and 0.58 (0.45-0.74) 5-24 weeks after mRNA-1273. The aHR profiles for common venous events PE and DVT were similar to those of composite venous thrombosis for all vaccine brands and doses. In contrast, there was a higher incidence of the known rare complication ICVT immediately after first dose of ChAdOx1 with greatest aHR 2 weeks after vaccination (5.92; 95%CI 4.07-8.63). There was no increase in incidence of ICVT after second dose of ChAdOx1 (Supplementary Table 8), or after any other vaccine brand.

Figure 2
Open in new tabFigure 2: Adjusted hazard ratios (aHRs) and 95% confidence intervals (95% CIs) for venous thromboses following COVID-19 vaccination, by dose and brand. Vertical lines depict 95% CIs; these are not visible when they are very narrow.

There was a higher incidence of thrombocytopenia after first dose of ChAdOx1 (Figure 3, Supplementary Table 8), compared with no vaccination, with greatest aHR 2 weeks after vaccination (2.07; 95%CI 1.67-2.58) but no increase in incidence after second dose of ChAdOx1 or after first or second dose of BNT-162b2, or after booster dose of mRNA-1273 following primary course of ChAdOx1 (Supplementary Table 9, 11, 12). There was a higher incidence of thrombocytopenia 13-24 weeks after a booster dose of BNT-162b2 following primary course of ChAdOx1 (aHR 2.16; 95%CI 1.26-3.69) (Figure 3, Supplementary Table 15).

Figure 3
Open in new tabFigure 3: Adjusted hazard ratios (aHRs) and 95% confidence intervals (95% CIs) for other cardiovascular events following COVID-19 vaccination, by dose and brand. Vertical lines depict 95% CIs; these are not visible when they are very narrow.

The incidence of SAH and mesenteric thrombus was similar or lower after all doses and brands of vaccination (Supplementary Tables 8-25), compared with before or without the corresponding dose. The incidence of mesenteric thrombus and SAH were markedly lower beyond 4 weeks after second dose of ChAdOx1 or BNT-162b2 vaccines. For example, aHRs for mesenteric thrombus 13-24 weeks after second dose were 0.64 (95%CI 0.49-0.83) for ChAdOx1 and 0.64 (0.49-0.83) for BNT-162b2. Corresponding aHRs for SAH 13-24 weeks after second dose were 0.64 (0.56-0.74) for ChAdOx1 and 0.84 (0.73-0.98) for BNT-162b2.

The incidence of myocarditis was higher after first dose of BNT-162b2 vaccine, with greatest aHR 1 week after vaccination (2.05; 95%CI 1.28-3.29) (Figure 3; Supplementary Table 9), higher one week after second dose of BNT-162b2 (aHR 3.14; 95%CI 2.04-4.85) (Supplementary Table 12) and higher after some mRNA booster vaccinations (for example, the aHR 1 week after booster BNT-162b2 vaccination was 1.65; 95%CI 1.07-2.57 following primary course of any of ChAdOx1, BNT-162b2 or mRNA-1273) (Supplementary Tables 19,21,22,24,25). Otherwise, the incidence of myocarditis after vaccination was similar to or lower than that before or without vaccination.

The incidence of pericarditis was higher after first dose of ChAdOx1, with greatest aHR 2 weeks after vaccination (1.74; 95%CI 1.04-2.91), after first dose of BNT-162b2 (for example, the aHR 13-24 weeks after vaccination was 1.50; 95%CI 1.17-1.92) (Supplementary Tables 8, 9), after second dose of BNT-162b2 (aHR 3-4 weeks after vaccination 2.42; 95%CI 1.62-3.62) (Supplementary Table 12), and after mRNA-based booster vaccination (for example, the aHR 2 weeks after booster dose of BNT-162b2 vaccine following primary course of any vaccine brand was 1.73; 95%CI 1.05-2.83) (Supplementary Tables 22, 25). Otherwise, the incidence of pericarditis after vaccination was similar to, or lower than, that before or without vaccination, with a halving in incidence beyond 4 weeks after second dose of ChAdOx1 (for example, the aHR after 13-24 weeks was 0.32; 95%CI 0.19-0.51) (Supplementary Table 11).

When follow-up was censored at the time of the US Centers for Disease Control (CDC) public announcement on myocarditis and pericarditis on 17th May 2021, the incidence of myocarditis and pericarditis after both first and second doses of ChAdOx1 and BNT-162b2 was similar to that before or without such vaccinations. One exception was a higher incidence of pericarditis 5-24 weeks after first dose BNT-162b2 vaccination (aHR 1.89; 95%CI 1.13-3.16) (Supplementary Figure 8; Supplementary Table 30). Note that the smaller number of events in this sensitivity analyses meant that aHRs were less precisely estimated than in the main analyses.

COVID-19 vaccination and thrombotic events within population subgroups

Subgroup analyses by age group, ethnic group, previous history of COVID-19, history of the event of interest and sex were conducted for composite arterial and composite venous outcomes (Supplementary Figures 2-7; Supplementary Tables 26-29). Associations between vaccination and thrombotic events were generally similar across subgroups, with the following exceptions. For composite arterial events after first doses of ChAdOx1 or BNT-162b2, aHRs in people with unknown ethnicity were higher than in people with known ethnicity (Supplementary Figures 2,4; Supplementary Table 26-27). In general, aHRs for composite arterial and composite venous events after first and second doses of ChAdOx1 or BNT-162b2 were higher in males than females (Supplementary Figures 2-5).

Discussion

This study used whole population longitudinal health records from over 45.7 million adults in England to quantify associations of first, second and booster doses of COVID-19 vaccine brands used during the first two years of the UK vaccine rollout with the incidence of arterial and venous thromboses, thrombocytopenia and myocarditis. Estimated hazard ratios were adjusted for a wide range of potential confounders. The incidence of thrombotic and cardiovascular complications was generally lower after each dose of each vaccine brand. Exceptions, consistent with previous findings that have been recognised by medicines regulators, included rare complications of the ChAdOx1 vaccine (ICVT and thrombocytopenia, due to vaccine-induced immune thrombocytopenia and thrombosis) and the mRNA vaccines (myocarditis and pericarditis). There were few differences between subgroups defined by demographic and clinical characteristics. These findings, in conjunction with the long-term higher risk of severe cardiovascular and other complications associated with COVID-19, offer compelling evidence supporting the net cardiovascular benefit of COVID vaccination.

The strengths of this study lie in the representativeness of the whole population data, offering an overview of thrombotic events after vaccination, as well as the comprehensive analyses of different vaccine dose and brand combinations in the general population. Consequently, the findings should apply to nations with comparable demographics and healthcare systems. The very large sample size facilitated estimation of associations with rare outcomes, within time periods after vaccination, and within population subgroups. The extensive longitudinal coverage of the health records also allowed examination of events after the first, second, and booster vaccinations. We addressed potential confounding by adjusting for a wide-range of demographic factors and prior diagnoses available in primary and secondary care records, defined using clinician-validated code lists that are accessible via our GitHub repository (available at https://github.com/BHFDSC/CCU002_06). Our analysis adhered to a pre-specified protocol with one deviation: we censored all analyses 26 weeks after vaccination to avoid interference from subsequent vaccinations 9.

This study has several limitations. First, residual confounding, including that linked to delayed vaccination in high-risk individuals, may persist despite extensive adjustments for available covariates. We were able to identify some, but not all people who were clinically vulnerable (and hence might have been eligible for earlier vaccination): for example, younger adults in long-stay settings could not be reliably identified. Second, we did not adjust for potential confounding by time-varying post-baseline factors that may have influenced receipt of vaccination and the outcomes of interest: for example, development of respiratory symptoms or being admitted into hospital leading to postponement of vaccination. Such confounding may explain estimated lower hazard ratios soon after vaccination 10. Third, ascertainment of some outcomes may have been influenced by public announcements from regulatory agencies, such as the European Medicines Agency Pharmacovigilance Risk Assessment Committee announcement 11 or the CDC announcement on myocarditis 12,13. This was addressed in sensitivity analyses for myocarditis and pericarditis, censoring follow-up at the time of public announcements of these adverse effects of vaccination, although the shorter follow-up times and corresponding smaller numbers of events in the restricted analyses meant that aHRs were estimated with reduced precision. Fourth, outcomes may be underreported, particularly from people in nursing homes or among those with severe health conditions, due to diagnostic challenges; also, routine electronic health records, not intended for research, may under-ascertain less severe, non-hospitalised events. Both forms of potential underreporting, however, are expected to be uncommon for hospitalised thrombotic events 14. Fifth, we restricted follow-up to 26 weeks after vaccination to prevent an influence of subsequent vaccinations on estimated associations and limit the impact of delayed vaccination on our findings. Horne et al. demonstrated selection bias in estimated HRs for non-COVID-19 death arising from deferred next-dose vaccination in people with a recent confirmed COVID-19 diagnosis or in poor health 9. Sixth, we did not address long-term safety of vaccination, or the impact of subsequent booster doses.

The incidence of arterial and venous thrombotic events was generally lower after COVID-19 vaccination than before or without vaccination. The higher incidence of cardiovascular events after COVID-19 is well-established 1417 and a plausible explanation for reductions in these events after vaccination is that vaccination prevents COVID-19, particularly severe COVID-19 1826. However, quantifying the mediating role of COVID-19 in cardioprotective effects of vaccination is beyond the scope of this paper.

The cohort design of our study allows estimation of incidence rates (see Table 2) and hence facilitates an understanding of the population-level impact of vaccination on common and rare cardiovascular events. However, our findings may be biased by unmeasured confounding, as discussed above. The self-controlled case series (SCCS) method aims to avoid the need to control for time-invariant confounders, but entails assumptions such as event-independent observation periods and medically-informed risk period definitions, and is particularly suited to short-term outcomes. Our results thus complement, and may be triangulated with, previous publications that examined cardiovascular implications of COVID-19 vaccination using the SCCS method 2729.

This extensive England-wide study offers reassurance regarding the cardiovascular safety of COVID-19 vaccines, with lower incidence of common cardiovascular events outweighing the higher incidence of their known rare cardiovascular complications. We found no novel cardiovascular complications or new associations with subsequent doses. Our findings support the wide uptake of future COVID-19 vaccination programs. We hope this evidence addresses public concerns, supporting continued trust and participation in vaccination programs and adherence to public health guidelines.

Methods

This analysis was performed according to a pre-specified analysis plan published on GitHub, along with the code lists to define variables and analysis code: https://github.com/BHFDSC/CCU002_06.

Study population

We analysed pseudonymised data, made available for the BHF Data Science Centre’s CVD-COVID-UK/COVID-IMPACT Consortium within the NHS England Secure Data Environment (7,27), which is a secure, privacy-protecting platform. This data consists of linked datasets including General Practice Extraction Service Extract for Pandemic Planning and Research (GDPPR), hospital admission data from Secondary Uses Service (SUS), Hospital Episode Statistics for admitted patient care (HES-APC), national laboratory COVID-19 testing data from the UK Health Security Agency (UKHSA) Second Generation Surveillance System (SGSS), Office for National Statistics (ONS) Civil Registration of Deaths (ONS deaths registry), medicines dispensed in primary care data and COVID-19 vaccination data.

The primary course of vaccination consists of the first and second vaccinations and, for certain groups such as people with severe immunosuppression, a third vaccination: this is distinct from booster vaccination that is given some time after the primary course 31.

People were included in the study if they were alive on 8th December 2020; aged 18-110 years inclusive; had a record in the GDPPR; recorded as male or female; and living in England. People with missing Lower-layer Super Output Area (LSOA) data were assumed to live in England. People were excluded if (1) they were vaccinated before 8th December 2020; (2) they were recorded as having a second dose and/or a booster or third dose, before or without records of first and second dose vaccinations respectively; (3) the interval between their first and second vaccination was less than 21 days 32; (4) they had mixed first and second vaccine brands where the second dose was given on or before 7th May 2021 3335; (5) the interval between their second and booster vaccinations was less than 90 days 36; (6) they had conflicting vaccination records or a situation code attached to any vaccination indicating that the vaccination was not given 37. We applied general quality checks, including removing people from the analysis who had nonsensical dates of birth or death (for details see https://github.com/BHFDSC/CCU002_06).

All eligible people were considered for the first dose analysis. Those who received a first dose of ChAdOx1, BNT-162b2, or mRNA1 were included in the second dose analyses. People who received the same vaccine brand for their first and second doses were included in the booster vaccination analyses.

Data for the study were extracted on 12th May 2022._The study spanned records from 8th December 2020, the start of the UK's vaccine rollout, to 23rd January 2022, the latest available date to ensure completeness of records from across all datasets at data extraction. Follow-up for first dose vaccination analyses began on 8th December 2020 for all people. For the second and booster vaccination analyses, follow-up commenced on the date of the preceding-dose vaccination. The term "index date" will henceforth denote the start of follow-up for each individual in each analysis.

Exposures

For the first and second dose analyses, the brands ChAdOx1, BNT-162b2 and mRNA-1273 were analysed separately. For the booster dose analyses, we analysed all combinations of primary courses (ChAdOx1, BNT-162b2, BNT-162b2/mRNA-1273) with booster vaccines (BNT-162b2, mRNA-1273, and BNT-162b2/mRNA-1273). (Supplementary Table 3).

Confounders

We considered the following confounders: age, sex, ethnic group, Index of Multiple Deprivation (2010 IMD deciles grouped into Deciles 1-4, 5-6, 7-10/missing), smoking status (never/ever, with missing classified as “never”), medical history (including acute myocardial infarction (AMI), diabetes, depression, obesity, cancer, chronic obstructive pulmonary disease (COPD), liver disease, chronic kidney disease, dementia, all stroke, all venous thromboembolic events and thrombophilia), major surgery in the last year, number of unique medical conditions in the last year, prior COVID-19 at index date, medications taken in the last 90 days (including antiplatelets, blood pressure lowering, lipid-lowering, oral anticoagulants, combined oral contraceptives (COCP) and hormone replacement therapy (HRT)) and clinical vulnerability (clinically extremely vulnerable/clinically vulnerable/neither). The medical history covariates were defined as a diagnosis of the condition before the index date except for diabetes which was additionally defined as a record of diabetic medication in the GDPPR data in the 90 days before the index date. Clinical vulnerability was defined on 8th December 2020. People were flagged as “clinically extremely vulnerable” using the SNOMED code 1300561000000107 38, and “clinically vulnerable” by identifying component conditions as applied in Table 3 of the COVID-19 chapter of the Green Book 31. Except for sex, clinical vulnerability and ethnic group, all other covariates were updated at individual-specific index dates for the dose 1, dose 2 and booster analyses. History of confirmed COVID-19 diagnosis was ascertained using established algorithms that combine information from SGSS, HES-APC, SUS, ONS deaths registry 39.

Outcomes

Eleven cardiovascular outcomes were analysed: AMI, ischaemic stroke, lower limb deep venous thrombosis (DVT), pulmonary embolism (PE), intracranial venous thrombosis (ICVT), mesenteric thrombus, portal vein thrombosis (PVT), any thrombocytopenia, subarachnoid haemorrhage & haemorrhagic stroke (SAH & HS), myocarditis and pericarditis. In addition, two composite outcomes were analysed: composite arterial (AMI, ischaemic stroke and other arterial embolism) and composite venous (PE, DVT, ICVT and PVT). We selected the earliest date of outcome event on or after index date from GDPPR, SUS, HES-APC and ONS deaths registry. We considered only the first/primary position from HES-APC and SUS and used the underlying cause from the death data, to differentiate acute, new events, from prevalent conditions. Further, we had previously found that aHRs for outcomes recorded as primary or secondary reason for admission or death were consistent with those from analyses of outcomes in the primary position5.

Statistical Analyses

Eligibility criteria and index dates for each vaccination course are detailed in Supplementary Table 3. We censored follow-up at the earliest of death, outcome event of interest, receipt of another vaccine brand, 26 weeks since the vaccination under consideration 9, or the study end date. Baseline demographic and clinical characteristics were detailed for each study cohort, and the number of outcome events and person-years of follow-up were quantified both before and after each vaccination, with incidence rates expressed per 100,000 person-years. Following the NHSE statistical disclosure control process, any counts less than 10 are presented as “10” and all numbers above 10 are rounded to the nearest 5. Percentages and incidence rates were calculated using rounded counts. Therefore, incidence rates corresponding to event counts less than 10 should be regarded as upper bounds and interpreted cautiously.

We analysed the time since vaccination to the first event for each outcome by fitting Cox models with a calendar time scale, with 8th December 2020 as time zero for the first dose analysis and the date of the previous dose as time zero for the second and booster dose analyses. We estimated aHRs and corresponding 95% CIs comparing follow-up after first, second and booster vaccine doses with follow-up before or without the corresponding vaccine dose for time intervals since vaccination: 1 weeks, 2 weeks, 3-4 weeks, 5-12 weeks, 13-24 weeks, and 25-26 weeks after vaccination. We expect that on average, and in the absence of bias, 5% of 95% CIs will exclude the true value of the aHR. Each comparison group included only individuals eligible to receive the vaccine brand and dose under consideration. For example, booster doses can only be received by individuals who have received a primary vaccine course. In the absence of any events in one of these intervals, we consolidated the periods into 1-4 weeks and 5-26 weeks after vaccination. For categorical confounders where fewer than two people in any category experienced an event, we combined categories where feasible and removed the confounder if not feasible. Furthermore, we stratified all models by region to account for between-region variation.

For computational feasibility, our analysis datasets comprised people who experienced the outcome (‘cases’) during follow-up and a randomly selected subsample with size twenty times the number of cases of those who did not experience the outcome. We applied inverse probability weights to account for this sampling method and used robust standard errors to compute confidence intervals. For each outcome and vaccination course, we estimated hazard ratios adjusting for (i) age and sex; and (ii) all measured confounders (maximally-adjusted).

We performed subgroup analyses for composite arterial and venous outcomes, by age group, ethnic group, prior history of the event (arterial: AMI, ischaemic stroke and other arterial embolism; venous: PE, DVT, ICVT, PVT and other DVT), prior history of confirmed COVID-19 diagnosis, and sex. We repeated the maximally-adjusted analyses for myocarditis and pericarditis, censoring follow-up on 17th May 2021, the date of the CDC‘s public announcement of their potential associations with COVID-19 mRNA vaccines.

Software

Analyses used SQL and Python (in Databricks, version 3.68), and RStudio (Professional) Version 1.3.1093.1 driven by R Version 4.0.3 (10th October 2020).

Ethical approval and information governance

The North East – Newcastle and North Tyneside 2 research ethics committee provided ethical approval for the CVD-COVID-UK/COVID-IMPACT research programme (REC No 20/NE/0161) to access, within secure trusted research environments, unconsented, whole-population, de-identified data from electronic health records collected as part of patients’ routine healthcare.

Supplementary Material

Acknowledgements

The British Heart Foundation Data Science Centre (grant No SP/19/3/34678, awarded to Health Data Research (HDR) UK) funded co-development (with NHS England) of the Secure Data Environment service for England, provision of linked datasets, data access, user software licenses, computational usage, and data management and wrangling support, with additional contributions from the HDR UK Data and Connectivity component of the UK Government Chief Scientific Adviser’s National Core Studies programme to coordinate national COVID-19 priority research. Consortium partner organisations funded the time of contributing data analysts, biostatisticians, epidemiologists, and clinicians. This work was supported by the Longitudinal Health and Wellbeing COVID-19 National Core Study (UKRI Medical Research Council MC_PC_20030 and MC_PC_20059).

SI was funded by the International Alliance for Cancer Early Detection, a partnership between Cancer Research UK C18081/A31373, Canary Center at Stanford University, the University of Cambridge, OHSU Knight Cancer Institute, University College London and the University of Manchester. SI and YL are supported by Cancer Research UK EDDPMA-May22\100062. AB is chief investigator of the NIHR-funded STIMULATE-ICP study (COV-LT2-0043). VW is supported by the COVID-19 Longitudinal Health and Wellbeing National Core Study, which is funded by the Medical Research Council (MC_PC_20059) and NIHR (COV-LT-0009), and the Medical Research Council Integrative Epidemiology Unit at the University of Bristol [MC_UU_00032/03].

Author Information

ku.ca.mac@02iyh

Author Contributions

Author contributions are reported below in line with the Contributor Roles Taxonomy (CRediT).

Conceptualization: J. A. C. S., A. M. W., V. W., W. N. W., S. I., T.-L. N.

Methodology: J. A. C. S., A. M. W., V. W., W. N. W., S. I., T.-L. N., A. A., E. H.

Software: S. I., T.-L. N., E. H., S. K.

Validation: S. I., T.-L. N., F. T.

Formal analysis: S. I., T.-L. N.

Investigation: S. I., T.-L. N.

Resources: C. S.

Data curation: S. I., T.-L. N.

Writing - Original Draft: S. I., T.-L. N., Y. L., W. N. W., J. A. C. S., A. M. W., V. W.

Writing - Review & Editing: S. I., T.-L. N., F. T., Y. L., H. A., A. A., E. H., R. D., S. K., S. D., A. B., K. K., C. S., W. N. W., J. A. C. S., A. M. W., V. W.

Visualization: S. I., T.-L. N., Y. L., W. N. W., J. A. C. S., A. M. W., V. W.

Project administration: S. I., T.-L. N., J. A. C. S., A. M. W., V. W.

Funding acquisition: J. A. C. S., A. M. W., C. S.

# The following authors contributed equally:

  • Samantha Ip
  • Teri-Louise North

These authors jointly supervised this work: Jonathan A C Sterne’ Angela M Wood, Venexia Walker

Competing Interests

KK was chair of the ethnicity subgroup of the UK Scientific Advisory Group for Emergencies (SAGE) and was a member of SAGE. CS is Director of the BHF Data Science Centre (whose main funding support comes from the British Heart Foundation) and Chief Scientist and Deputy Director at Health Data Research UK, the UK’s Institute for Health Data Science. She played a key role in co-developing NHS England’s national secure data environment and leads the CVD-COVID UK/COVID-IMPACT Consortium, which enabled this work. WW is supported by the Chief Scientists Office (CAF/01/17) and Stroke Association (SA CV 20100018). WW has given expert testimony to UK courts. The remaining authors declare no competing interests.

Data Availability

Data used in this study are available in NHS England’s Secure Data Environment service for England for England, but as restrictions apply are not publicly available (https://digital.nhs.uk/services/secure-data-environment-service). The CVD-COVID-UK/COVID-IMPACT programme led by the BHF Data Science Centre (https://bhfdatasciencecentre.org/) received approval to access data in the NHS England’s SDE service for England from the Independent Group Advising on the Release of Data (IGARD) (https://digital.nhs.uk/about-nhs-digital/corporate-information-and-documents/independent-group-advising-on-the-release-of-data) via an application made in the Data Access Request Service (DARS) Online system (reference: DARS-NIC-381078-Y9C5K; https://digital.nhs.uk/services/data-access-request-service-dars/dars-products-and-services). The CVD-COVID-UK/COVID-IMPACT Approvals & Oversight Board ( https://bhfdatasciencecentre.org/areas/cvd-covid-uk-covid-impact/) subsequently granted approval to this project to access the data within NHS England’s SDE service for England. The de-identified data used in this study were made available to accredited researchers only. Those wishing to access the data should contact ku.ca.kurdh@csdfhb.

Code Availability

All code and code lists are shared openly for review and re-use under an MIT open license (https://github.com/BHFDSC/CCU002_06/).

Notes

1 There were no AMI events during weeks 13-26 after mRNA-1273 booster vaccination, so follow-up is grouped as 1-4 and 5-26 weeks post-vaccination.

2 All counts rounded to the nearest 5 and counts less than 10 displayed as ‘10’

References

History

  • Posted February 13, 2024.

Citations & impact 


This article has not been cited yet.

Impact metrics

Alternative metrics

Altmetric item for https://www.altmetric.com/details/159498767
Altmetric
Discover the attention surrounding your research
https://www.altmetric.com/details/159498767

Data 


Data behind the article

This data has been text mined from the article, or deposited into data resources.

Funding 


Funders who supported this work.

British Heart Foundation (1)