Abstract
Background Substantial evidence suggests that severe Covid-19 leads to an increased risk of Venous Thromboembolism (VTE). We aimed to quantify the risk of VTE associated with ambulatory Covid-19, study the potential protective role of vaccination, and establish key clinical and genetic determinants of post-Covid VTE.
Methods We analyzed a cohort of ambulatory Covid-19 patients from UK Biobank, and compared their 30-day VTE risk with propensity-score-matched non-infected participants. We fitted multivariable models to study the associations between age, sex, ethnicity, socio-economic status, obesity, vaccination status and inherited thrombophilia with post-Covid VTE.
Results Overall, VTE risk was nearly 20-fold higher in Covid-19 vs matched non-infected participants (hazard ratio [HR] 19.49, 95% confidence interval [CI] 11.50 to 33.05). However, the risk was substantially attenuated amongst the vaccinated (HR: 2.79, 95% CI 0.82 to 9.54). Older age, male sex, and obesity were independently associated with higher risk, with adjusted HRs of 2.00 (1.61 to 2.47) per 10 years, 1.66 (1.28 to 2.15), and 1.85 (1.29 to 2.64), respectively. Further, inherited thrombophilia led to an HR 2.05, 95% CI 1.15 to 3.66.
Conclusions Ambulatory Covid-19 was associated with a striking 20-fold increase in incident VTE, but no elevated risk after breakthrough infection in the fully vaccinated. Older age, male sex, and obesity were clinical determinants of Covid-19-related VTE. Additionally, inherited thrombophilia doubled risk further, comparable to the effect of 10-year ageing. These findings reinforce the need for vaccination, and call for targeted strategies to prevent VTE during outpatient care of Covid-19.
Introduction
Numerous hospital-based studies and case series have demonstrated a high risk of venous thromboembolism (VTE) in severe Covid-19 patients. A recent meta-analysis reported a pooled VTE rate of 14.7% and 23.2% among those admitted to hospital and in intensive care units, respectively.1 Additionally, emerging randomized controlled trials have shown the benefit-risk of anticoagulation for Covid-19 patients at different stages of the diseases,2 and clinical societies have recommended initiating routine antithrombotic therapy during hospital admission.3
In contrast, fewer clinical interventions have been implemented to prevent VTE amongst ambulatory patients with SARS-CoV-2 infection, partially due to conflicting findings on the association between infection and VTE risk, and mixed evidence of benefit from oral anticoagulation.4 Given the ongoing global rollout of vaccines, relaxation of public health restrictions, and widespread highly transmissible Omicron variant, the absolute number of milder Covid-19 cases managed in ambulatory settings continues to rise worldwide.5,6 These collectively suggest that prophylaxis, including timing and dosing regimens, requires further refinement, particularly in the outpatient setting.2,7–10 Moreover,11 a lack of insight into the impact of clinical, socioeconomic, and genetic risk factors for infection-related VTE persists.
This study aimed to (1) quantify the magnitude of acute VTE risk attributable to SARS-CoV-2 infection identified in ambulatory settings; and (2) investigate the clinical and genetic determinants of VTE risk after SARS-CoV-2 infection.
Methods
We included UK Biobank (UKBB) participants from England who were alive on 01 March 2020. All participants provided written informed consent at the UKBB cohort recruitment. This study received ethical approval from UKBB Ethics Advisory Committee (EAC) and was performed under the application of 65397.
Data sources and study cohorts
We obtained data from UK Biobank comprising multiple linked sources, including baseline surveys conducted between 2006-2010, individual genetic data, primary care electronic medical records, hospital inpatient data from Hospital Episode Statistics, diagnostic Covid-19 test data from Public Health England and death records from the national death registry (Office of National Statistics).
We curated an infected cohort by enrolling individuals with positive polymerase-chain-reaction SARS-CoV-2 results confirmed between 1 March 2020 and 30 September 2021. Participants who were never tested or only had negative test results were classified into the non-infected cohort. The index date was the date of the first positive specimen sample for the infected cohort, and a random date during the same calendar period for the matched non-infected. Participants with historical VTE or who used antithrombotic drugs one year before index dates were. Additionally, we excluded those in the infected cohort who tested positive in hospital settings.
Inherited thrombophilia
Information on genotyping and imputation procedures in UK Biobank has been detailed in previous studies.12 Briefly, genome-wide single nucleotide polymorphisms (SNPs) were genotyped using two closely related purpose-designed arrays (the UK BiLEVE Axiom array and UK Biobank Axiom array). We defined inherited thrombophilia carriers as having any of two risk SNP variants in Factor V Leiden (rs6025) or Prothrombin G20210A (rs1799963). We also defined a positive genetic control exposure by calculating a 297-SNPs polygenic risk score (PRS) for VTE that did not include these two variants13 (see Supplementary Appendix for details).
Covariates
We pre-specified a list of covariates for adjustment, including demographics (age, sex, ethnicity), socio-economic status measured by the Index of Multiple Deprivation (a summary deprivation measurement used in England containing seven aspects in crime, education, employment, health, housing, income and living environment),14 body mass index (BMI), medications prescribed within one year before the index date, and comorbidities recorded in primary care records (Table 1). The number of hospital admissions in the past one year (proxy of healthcare utilization) and vaccination status (no, partial 1-dose, or full 2-dose vaccination) were also studied.
Outcomes
Incident VTE, comprising either deep vein thrombosis or pulmonary embolism, was identified using ICD-10 codes based on hospital records. Eligible participants were followed up for up to 30 days after the index date, given that VTE occurring after 30 days were much less likely to be due to a direct consequence of SARS-CoV-2 infection.
Statistical analyses
We used propensity score (PS) matching to minimize confounding in the study of the association between SARS-CoV-2 infection and VTE. We fitted multivariable logistic regression models to estimate PS as the probability of infection based on the predefined covariates. We then matched infected with non-infected individuals with a ratio of 1:5 based on PS values with a caliper width of up to 0.2 standard deviations of the logit of the PS, with exact-matching on index dates.15,16 We assessed covariate balance between the cohorts before and after matching using absolute standardized mean differences (SMD) and specified SMD>0.1 as relevant imbalances.17 Cause-specific Cox survival models were applied to estimate hazard ratios (HRs) for VTE according to exposure, in which death was considered a competing risk.18
To study clinical determinants, we fitted multivariable Cox models of 30-day VTE in the ambulatory Covid-19 cohort, including age, sex, ethnicity, socio-economic status, obesity (BMI less than vs equal or more than 30 kg/m2) and vaccination status. For the analysis of inherited thrombophilia’ effect on post Covid-19 VTE, we adjusted Cox models for age, sex and genetic ancestry (quantified by the first three principal components), assuming that genetic variants were independent of baseline characteristics. We introduced a positive exposure control (PRS for VTE) and a negative outcome (arterial thromboembolism [ATE]) experiment to detect residual confounding and potential unresolved bias.19
A sensitivity analysis was performed repeating the modelling of clinical determinants with VTE in the non-infected cohort, where no association is expected between vaccination and Covid-19-unrelated VTE.
All statistical tests were 2-sided, where a P = 0.05 or a 95% CI that did not contain unity were considered statistically significant for the primary analyses. All analyses and data visualizations were conducted using R statistical software. Genetic data management and quality controls were performed using Plink 1.9.20
Results
Baseline characteristics
Out of 407,311 UKBB participants, 26,210 were infected with SARS-CoV-2 between 1 March 2020 and 30 September 2021. After applying exclusion criteria, 21,724/ 26,210 (83%) infected and 318,754/ 407,311 (78.23%) non-infected were eligible for analyses (Supplementary Results 1). For all infections, 2,877 (13 %) and 18,847 (87 %) tested in hospital and outpatient settings, respectively. Only the latter were included for subsequent analyses.
Baseline characteristics by infection status are summarized in Table 1. Overall, ambulatory Covid-19 participants were younger than the uninfected (mean age [sd] 64.32 [8.03] vs 67.85 [8.04]), more likely male (43.8% vs 41.5%), from non-White ethnicity (12.0 % vs 6.4%), and more socioeconomically deprived and obese. After PS matching, almost all ambulatory Covid-19 cases (98.7%) were successfully matched to at least one non-infected counterpart, and all covariates became balanced. In addition, index dates and calendar periods were accurately aligned between the cohorts, as depicted in Supplementary Results 2.
Association of SARS-CoV-2 infection with incident VTE
Figure 1. depicts cumulative incidence of VTE according to infection status, showing an early separation of the matched cohorts, with continued divergence over time. A total of 73 and 17 VTE events were seen within 30 days amongst the ambulatory Covid-19 patients and matched uninfected, corresponding to incidence rates of 51.00 and 2.61 per 1,000 person-years, respectively. Survival analyses (Table 2) suggested that SARS-CoV-2 infection was associated with a striking increase in VTE risk: HR: 19.49, 95% CI 11.50 to 33.05. The observed risk was more pronounced in the unvaccinated patients (HR: 29.37, 95% CI 15.11 to 57.08) and substantially mitigated in the fully vaccinated (HR: 2.79, 95% CI 0.82 to 9.54).
Clinical determinants of post-Covid-19 VTE
The associations between socio-demographic and clinical factors (including vaccination status) and risk of post-Covid VTE are shown in Figure 2. Older participants had a higher risk, with an approximate doubling of risk per each 10-year increase in age (adjusted HR: 2.00, 95% CI 1.61 to 2.47). Men were at higher risk than women (adjusted HR: 1.66, 95% CI 1.28 to 2.15), and people with obesity at a higher risk than non-obese (adjusted HR: 1.85, 95% 1.29 to 2.64). These associations were equally seen for Covid-19-unrelated VTE in both direction and magnitude (Figure 2; Supplementary Results 3).
Additionally, 2-dose vaccination was associated with a substantially lower risk of Covid-19-related VTE (adjusted HR: 0.13, 95% 0.06 to 0.28). A sensitivity analysis demonstrated no association between vaccination status and Covid-19-unrelated VTE (Figure 2; Supplementary Results 3).
Inherited thrombophilia and post-Covid-19 VTE
Among 21,055 infected participants with complete genetic data, 1,287 (6.11%) had inherited thrombophilia, with 909 (4.32%) and 392 (1.86%) carrying risk variant/s of Factor V Leiden and Prothrombin G20210A, respectively (Table 3). The frequency of these genetic variations in the infected cohort was like that in the overall UKBB cohort and consistent with reports from previous literature (Supplementary Results 4). No differences in any of the measured covariates (except for ethnicity), including socio-demographics, medications or comorbidities, were observed when comparing those with vs without inherited thrombophilia (Supplementary Results 5).
Participants with inherited thrombophilia had a higher risk of VTE following SARS-CoV-2 infection than those without (adjusted HR: 2.05, 95% CI 1.15 to 3.66). For each risk variant, the adjusted HR was 2.17 (95% CI 1.13 to 4.15) for Factor V Leiden carriers and 1.52 (95% CI 0.48 to 4.79) for Prothrombin G20210A carriers. Also, individuals with higher PRS values had greater VTE risk: adjusted HR per 1-SD increase of PRS 1.33, 95% 1.11 to 1.59 (Table 3). As expected, no associations were observed between inherited thrombophilia and the negative control outcome ATE, with adjusted HR ranging from 0.91 (95% CI 0.29 to 2.84) to 0.99 (95% CI 0.49 to 2.01).
Discussion
In this community-based cohort of UKBB participants, including 26,210 with PCR-confirmed ambulatory Covid-19 cases and 381,063 matched contemporary non-infected controls, we found SARS-CoV-2 infection was associated with an approximately 20-fold increase in VTE risk within 30 days of a positive test. However, this risk was largely attenuated in the fully vaccinated participants who then suffered a breakthrough infection.
Known clinical determinants of VTE, including older age, male sex, and obesity, applied to post-Covid-19 VTE. Complete vaccination was associated with an almost 90% relative risk reduction of Covid-19-related VTE. As expected, vaccination did not impact VTE risk in the uninfected peers. Finally, inherited thrombophilia carriers had an additional double risk of post-Covid-19 VTE compared to non-carriers, equivalent to the excess risk attributable to an increase of 10 years of age.
Our finding of a substantially higher incidence of VTE following ambulatory Covid-19 disagreed with a previous meta-analysis of seven heterogeneous small Covid-19 cohorts, which had suggested that mild Covid-19 was not a risk factor for VTE 21. However, our data are in line with larger self-controlled case series studies22,23 that better accounted for within-person confounding and consistently showed orders of magnitude increases in VTE risk after SARS-CoV-2 infection.
Public interest and concerns have been placed predominantly on vaccine-related rare thromboembolic events 24, which has led to vaccine hesitancy and restrictions on its use.25 Our study found that the vaccination can offset SARS-CoV-2-induced VTE risk even if people get a breakthrough infection. This evident benefit should not be ignored in the ongoing global vaccination campaigns.
For the first time, we showed that inherited thrombophilia conferred a double risk of Covid-19-related VTE, echoing clinical findings of elevated factor V activity in severe Covid-19 patients.26–28 Although a relatively high proportion of congenital thrombophilia was previously detected in a small pilot study of 87 Covid-19 patients, the minimal sample size precluded any robust inference.29
We analyzed rich linked data combining extensive community SARS-CoV-2 testing, well-recorded vaccination status, ambulatory and hospital-based clinical outcomes, as well as large scale genotyping data readily available for UKBB participants. The results of this analysis have many noteworthy implications. First, VTE risk management needs re-evaluation for milder ambulatory Covid-19 patients. With emerging evidence and guidelines focusing on VTE prophylaxis for hospitalized Covid-19 patients, further work is now necessary to mitigate risk in the community. Second, although the aetiology of post-Covid VTE is complex and multifaceted, our findings elucidated the role of Factor V and possibly Prothrombin proteins as contributing factors. Third, although genetic testing of inherited thrombophilia for VTE prevention has been previously discussed in many clinical scenarios,30,31 this newly identified effect on Covid-19-related VTE, comparable to a 10-year ageing effect, supports the potential value of targeted genetic screening for thrombophilia in the infected older adults. Fourth, traditional known clinical risk factors of VTE are still informative to identify subgroups of people at high risk of Covid-19-related VTE. Finally, our data demonstrate the significant impact of vaccination to minimize the risk of Covid-19 VTE.
This study also has some limitations. Residual confounding cannot be ruled out in this observational study, although robust statistical approaches for causal inference was applied, including PS matching and knowledge-driven negative control analyses. Participants recruited in UKBB might not fully represent the general population and were likely healthier. Although those Covid-19 participants were from non-hospital settings, they were likely to be symptomatic and more severe cases. The extent to which asymptomatic infection is associated with VTE risk warrants further investigation. Finally, the estimates from our analyses were an average and mixed effect of several SARS-CoV-2 strains, which should be cautiously extrapolated to novel variants, such as Omicron.
Community-based ambulatory Covid-19 is associated with a striking 20-fold excess risk of VTE. This risk was much higher among unvaccinated individuals, and increased with older age, in men, and in patients with obesity. Inherited thrombophilia further doubled VTE risk, comparable to a 10-year ageing effect. Our findings call for targeted prevention strategies for post-Covid VTE in outpatient settings, and suggest an aetiological role of inherited thrombophilia. The clinical utility of genetic testing to stratify infected patients and inform personalized thromboprophylaxis warrants further research.
Data Availability
Bonafide researchers can apply to use the UK Biobank dataset by registering and applying at http://ukbiobank.ac.uk/register-apply/. The datasets generated and/or analysed during the current study are available from the corresponding author on reasonable request
Declaration of interests
Funding
Mr Xie is funded through Jardine-Oxford Graduate Scholarship and a titular Clarendon Fund Scholarship. DG is supported by the British Heart Foundation Research Centre of Excellence (RE/18/4/34215) at Imperial College London and by a National Institute for Health Research Clinical Lectureship (CL-2020-16-001) at St. George’s, University of London. Prof Prieto-Alhambra is funded through an NIHR Senior Research Fellowship (grant SRF-2018-11-ST2-004), and received partial support from the Oxford NIHR Biomedical Research Centre. APU has received funding from the Medical Research Council (MRC) [MR/K501256/1, MR/N013468/1].
Conflicts of interest
Prof Prieto-Alhambra that his research group has received grants and advisory or speaker fees from Amgen, Astellas, AstraZeneca, Chiesi-Taylor, Johnson & Johnson, and UCB; and that Janssen, on behalf of Innovative Medicines Initiative–funded European Health Data Evidence Network and European Medical Information Framework consortiums and Synapse Management Partners, have supported training programs, open to external participants, organized by his department. DG is employed part-time by Novo Nordisk.
Roger Paredes has participated in advisory boards for Gilead, MSD, ViiV Healthcare, Theratechnologies and Lilly. His institution has received research support from Gilead, MSD, and ViiV Healthcare.
Ethical approval
All participants provided written informed consent at the UKBB cohort recruitment. This study received ethical approval from UKBB Ethics Advisory Committee (EAC) and was performed under the application of 65397.
Data sharing
Bonafide researchers can apply to use the UK Biobank dataset by registering and applying at http://ukbiobank.ac.uk/register-apply/.The datasets generated and/or analysed during the current study are available from the corresponding author on reasonable request.
Transparency declaration
The lead author affirms that this manuscript is an honest, accurate, and transparent account of the study being reported; that no important aspects of the study have been omitted.
Contributors
D.P.A., J.Q.X., D.G, R.P were responsible for the study design. J.Q.X. did the data analyses, and A.P.U. checked the statistical codes. J.Q.X. and D.P.A. drafted the manuscript, and all co-authors reviewed and approved it for submission.