Objectives To describe the risk of venous thromboembolism (VTE), and risk factors for VTE, in people with immune-mediated inflammatory diseases (IMID) (ulcerative colitis, Crohn’s disease (CD), rheumatoid arthritis (RA) and psoriatic arthritis (PsA)), compared with a matched control population.
Methods A total of 53 378 people with an IMID were identified over 1999–2019 in the UK Royal College of General Practitioners (RCGP) Research and Surveillance Centre (RSC) primary care database and were matched to 213 512 people without an IMID. The association between the presence of any IMID, and each IMID separately, and risk of VTE was estimated using unadjusted and multivariable-adjusted Cox proportional hazards models. The prevalence of VTE risk factors, and associations between VTE risk factors and risk of VTE, were estimated in people with and without an IMID.
Results People with an IMID were at increased risk of VTE (adjusted HR [aHR] 1.46, 95% CI 1.36,1.56), compared with matched controls. When assessing individual diseases, risk was increased for CD (aHR 1.74, 95% CI 1.45 to 2.08), ulcerative colitis (aHR 1.27, 95% CI 1.10 to 1.45) and RA (aHR 1.54, 95% CI 1.40 to 1.70) but there was no evidence of an association for PsA (aHR 1.21, 95% CI 0.96 to 1.52). In people with an IMID, independent risk factors for VTE included male sex, overweight/obese body mass index, current smoking, history of fracture, and, across study follow-up, abnormal platelet count.
Conclusions VTE risk is increased in people with IMIDs. Routinely available clinical information may be helpful to identify individuals with an IMID at increased future risk of VTE.
Observational study registration number Clinicaltrials.gov (NCT03835780).
- Autoimmune Diseases
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.
Statistics from Altmetric.com
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.
Venous thromboembolism (VTE), comprising pulmonary embolism (PE) and deep vein thrombosis (DVT), is relatively common, with an incidence in the general population of around 3 cases per 1000 patient years.1 It is associated with significant morbidity and mortality.2 3
Inflammation increases the risk of VTE,4 and observational data demonstrate higher VTE rates in individuals with immune-mediated inflammatory diseases (IMID) including ulcerative colitis (UC), Crohn’s disease (CD) and rheumatoid arthritis (RA) compared with the general population.5–9 Evidence for VTE risk in other inflammatory diseases, including psoriatic arthritis (PsA) is more limited.6 Risk factors for VTE have been well described in the general population, and include obesity, fractures, surgery, use of oral corticosteroids and hormone therapies.10 11 and high platelet count which has been reported to be a risk factor for VTE in hospital inpatients,12 and is recognised as a marker of inflammation in inflammatory bowel disease and RA.13 14 There has however been little systematic interrogation of whether VTE risk factors convey the same risk in individuals with and without an IMID.15
In this study, we set out to use a large UK primary care database to establish the excess risk of VTE in people with an IMID (UC, CD, RA and PsA) compared with a control population without any of these conditions. We then compared the prevalence of traditional VTE risk factors in people with and without an IMID, and the associations between these features and future risk of VTE.
We performed a cohort study using matched populations to compare VTE risk in adults with an IMID (UC, CD, RA and PsA) and controls between 1999 and 2018 inclusive, using UK population-based primary care data.
Data were sourced from the Royal College of General Practitioners Research (RCGP) and Surveillance Centre (RSC) database. RCGP RCS derives data from a representative network of general practices distributed across England, currently covering a registered population of 2 million people.16 RCGP RSC contains information on demographics, clinical features and diagnoses, laboratory tests and prescriptions, and studies using RCGP RSC data have been published across a range of chronic diseases.17–20
Adults (aged ≥18) were eligible for inclusion if registered with a general practice between January 1, 1999 and December 31, 2018, with at least one consultation over that period (to minimise the impact of ‘ghost’ patients), and no history of VTE.
Definition of the exposed cohort with IMID
The exposed cohort was defined as all individuals with an existing or incident diagnosis of UC, CD, RA or PsA in the RCGP RSC database over the study period. UC, CD and RA were identified using Read diagnostic codes and algorithms previously validated by review of individual patient records or collection of questionnaires from general practitioners in UK primary care studies.21–25 In the absence of a validated method to identify the presence of PsA from UK primary care data, this was identified using a Read code list generated in accordance with published guidance.26 27 The index date for start of follow-up for exposed individuals began on the latest of the date of diagnosis indicated by first diagnostic code, January 1, 1999, or 180 days after practice registration.
Definition of the matched unexposed cohort
People with an IMID were matched at their index date with four unexposed individuals at general practice level by current age (per year), sex and years since practice registration (nearest neighbour matching, with replacement). The eligible pool of unexposed individuals at each index date comprised individuals registered at that date with no history of an IMID and at least 1 year of follow-up in RCGP RCS (to minimise the risk they had a non-recorded existing IMID diagnosis). Follow-up for each matched individual started on the index date of their matched case. Individuals with an incident diagnosis of an IMID during the study period were included in the pool of eligible unexposed individuals, but if matched were censored on the date of their diagnosis of an IMID; that is, these individuals were eligible to contribute to unexposed person time before their diagnosis of an IMID. Follow-up for each individual ended at the earliest of the study end-date (December 31, 2018), the date an individual was transferred from an included practice, date of death or the date an individual developed an outcome of interest.
The primary outcome was a diagnosis of VTE (a composite of PE or DVT). The secondary outcomes were individual diagnoses of PE and DVT. When both PE and DVT occurred on the same date this was classified as PE. Outcomes were identified using updated Read code lists previously validated by review of patient records and provision of general practitioner questionnaires.28 Risk of each outcome was compared between individuals with an IMID and the matched control population, and between individuals with UC, CD, RA and PsA and their matched counterparts.
Recorded characteristics and VTE risk factors
Baseline features comprised sociodemographic characteristics, clinical VTE risk factors, comorbidities and medication use. VTE risk factors were selected based on existing literature demonstrating an established association with VTE6 10 and clinical expertise. Clinical VTE risk factors were body mass index (BMI), smoking status, alcohol use, evidence of reduced mobility, thrombophilia, fracture of the lower limb and family history of VTE. Socioeconomic status was defined using index of multiple deprivation (IMD), the official national measure of socioeconomic status in the UK.29 Ethnicity was extracted from the primary care record and grouped into major UK ethnic groups: white, black, Asian, mixed and others.30 BMI, smoking status and alcohol use were defined using the most recently recorded data prior to the index date. Diagnostic codes were used to define the following baseline comorbidities: hypertension, hyperlipidaemia, type 2 diabetes, peripheral vascular disease, cardiovascular disease (atrial fibrillation, angina, myocardial infarction, congestive heart failure), stroke, malignancy, chronic obstructive pulmonary disease (COPD), chronic kidney disease (CKD) (stages 3–5), liver disease and thrombophilia. Type 2 diabetes was identified using an algorithm developed for use within RCGP RSC.31 Read codes used to describe cardiovascular disease within RCGP RSC have been previously reported.30 32 Platelet count measures were extracted at baseline (the most recent value up to 2 years prior to the index date) and across study follow-up.
We examined the following medications commonly used for the management of IMIDs: non-steroidal anti-inflammatory drugs (NSAIDs), oral corticosteroids, non-biologic immunosuppressant medications and biologic therapies recorded in primary care. We also examined hormone therapy (hormonal contraceptives, hormone replacement therapy (HRT)), antiplatelet agents (aspirin or ADP receptor inhibitors), warfarin, direct oral anticoagulants (DOACs), and statins. Hormonal contraceptives comprised only combined oestrogen and progestogen preparations; progesterone only contraceptives were not included as these preparations are not associated with VTE.33 HRT comprised systemic oestrogen only preparations. Active prescribing was defined as an issued prescription in the 3 months preceding and/or 1 month after the index date.
We estimated the risk of VTE, the primary outcome, using unadjusted Cox proportional hazards models, stratified by matched set (exposed cohort vs unexposed cohort), to provide overall HRs with 95% CI for the association. Models were subsequently adjusted for all sociodemographic, clinical and VTE risk factors, as described above, in multivariable analysis. We then repeated the same analyses for PE and DVT as separate endpoints and each condition (UC, CD, RA and PsA) separately. Proportional hazards assumptions for each model were checked graphically by plotting Schoenfeld residuals.
VTE risk factors
We used multivariable Cox models to examine the influence of baseline-recorded characteristics and VTE risk factors on risk of VTE. Models were run separately in cohorts with and without an IMID. If baseline platelet count was significant in the model, we proposed to explore the impact of changing platelet count over study follow-up on VTE risk in time-updated analysis, by including platelet count as a time-updated exposure in unadjusted and multivariable-adjusted Cox models. Platelet count was both categorised as low (<150×109/L), normal (150–400×109/L) or high (>400×109/L), and analysed continuously using a restricted cubic spline pre-specified with 3 knots.
Sensitivity to the introduction of the Quality and Outcomes Framework (QOF), an incentivised programme to monitor clinical and health improvement indicators for general practice that rewards completeness of electronic coding,34 was tested by repeating the main analyses with the study follow-up beginning on January 1, 2004. All statistical analyses used R (version 3.4.1).
A total of 53 378 people with an IMID were included, of whom 14 182 (26%) had a first diagnosis of UC, 9489 (18%) CD, 23 410 (44%) RA and 6297 (12%) PsA (table 1). Matched controls comprised 213 512 people without an IMID of interest. Average study follow-up was 8.2 (SD 6.2) years.
People with an IMID were similar in characteristics to their matched counterparts (table 1). Several comorbidities were more common in the exposed group including type 2 diabetes, COPD and chronic liver disease. BMI was similar although differences were observed between individuals with an IMID; more people with PsA were obese (32.9%) than people with UC (16.7%) or CD (14.6%), and more people with CD were underweight (5.6%) compared with other IMIDs (range 1.0–2.6%). Use of NSAIDs, corticosteroids and immunosuppressive medications were, as expected, considerably higher in the IMID group.
Risk of VTE
Unadjusted VTE event rates were higher in the IMID group (34.9 [95% CI 33.2 to 36.7] per 10 000 person-years) compared with controls (21.7 [95% CI 21.0 to 22.4] per 10 000 person-years, p<0.001) (figure 1); 1532 (2.9%) people with an IMID developed VTE compared with 3804 (1.8%) controls. Table 2 reported study follow-up and outcome events for the primary VTE outcome and the secondary outcomes of PE and DVT.
In the primary outcome analysis, adjusted models demonstrated an association between UC, CD and RA and the development of VTE, with the strength of association greatest for people with CD. Associations were consistent in analyses of separate PE and DVT endpoints (table 2, online supplemental figure 1). For PsA, a significant increase in risk was seen only for the DVT endpoint (table 2). Sensitivity analysis exploring the impact of QOF demonstrated primary results were consistent with follow-up beginning in 2004 (online supplemental table 1).
Table 3 shows associations between record characteristics and risk of developing VTE in people with an IMID and matched controls. Increasing age at entry, being overweight/obese and thrombophilia history were associated with an increased risk of VTE in both groups. Associations between VTE risk factors and risk of VTE differed, with male sex, history of fracture, current smoking and alcohol abstinence associated with an increased VTE risk only in the IMID group. Reduced mobility (increased risk) and Asian ethnicity (decreased risk) were risk factors only in controls. COPD, chronic liver disease, peripheral vascular disease and family history of VTE were all associated with an increased risk of VTE only in controls. CKD was associated with increased risk of VTE only in people with an IMID. For medication, warfarin, DOACs, corticosteroids, and, in females, use of oral contraceptives, were associated with an increased risk of VTE in both groups. Statins (decreased risk) and NSAIDs (increased risk) were risk factors only in the IMID group. Risk of VTE was increased in people with an IMID and low number of platelets (<150×109/L), and in people without an IMID and high number of platelets (>400×109/L). In the IMID group, relative to people with UC, risk of VTE was increased in CD only.
Association of platelet count across study follow-up with risk of VTE
To further interrogate the relationship between baseline platelets and VTE in each group, a time-updated analysis was undertaken. Individuals with at least one platelet count were included in the analysis (96% of those with an IMID and 75% of those without an IMID). Platelet count across study follow-up was initially categorised as low, normal or high (table 4). High and low platelet counts were more common in individuals with an IMID (proportion of individuals with 1+ one high platelet count 29.2%, low count 10.9%) compared with those without an IMID (high count 11.5%, low count 7.4%).
Higher time-varying platelet counts were associated with an increased risk of VTE in individuals with and without an IMID (table 4, figure 1). Figure 2 confirms the association between time-updated lower and higher platelet count and higher risk of VTE in both groups when modelling platelet count as a non-linear continuous variable; a positive association was also seen for platelet counts <200×109/L.
Our study shows that VTE is more common in people with UC, CD, RA, and PsA compared with people without these IMIDs. UC, CD and RA were independently associated to a similar degree with increased risk of VTE, while the wide CI for PSA suggests we lacked statistical power to detect a difference in this group. Established risk factors for VTE had a similar prevalence in people with an IMID compared with the wider population, and similar strengths of association were observed in people with and without an IMID for higher age, being overweight or obese, thrombophilia, malignancy, and corticosteroid use. Notable differences were also observed; only in people with an IMID did we find evidence that male sex, current smoking, CKD, and history of fracture were independent risk factors for VTE. Across study follow-up, abnormal platelet counts were found to be independently associated with risk of VTE in both groups but were substantially more common in people with an IMID.
Our study for the first time assesses the risk of developing VTE across four common IMIDs using the same study design. VTE incidence in this study was similar to that previously reported in the UK,1 and a similar difference in incidence rates between people with PsA, RA and psoriasis and matched controls was recently reported using UK primary care data.6 Results are in keeping with previous studies that have consistently found people with RA to be at increased risk of VTE.5 6 9 35 Ogdie et al conducted the first observational study of VTE risk in people with PsA and similar to our study, likely lacked power to detect a difference for overall VTE risk,6 suggesting further evaluation of VTE risk in patients with PsA in even larger cohorts is an important area for future research.
We also demonstrate an interesting u-shaped association between platelet count and VTE risk, with both high and low platelet count demonstrated to be markers of increased risk compared with normal platelet count in people with IMIDs managed in primary care. Given initial positive associations with baseline platelet count, and the recognised interaction between inflammatory cytokines and platelet function,36 we explored this association in depth using time-updated platelet counts across study follow-up to further delineate thresholds of risk/association with time to VTE. Time-updated high and low platelets were independently associated with risk of VTE in both people with and without an IMID; however, high and low platelet counts were much more common in people with an IMID, suggesting particular clinical utility in this group. Although the direction of effects was the same for the baseline and time-varying platelet analysis, differences in statistical significance and effect size may relate to the increased power and greater predictive ability gained from incorporating time-updated platelet measures. Our findings for high platelet count are in keeping with previous studies that have demonstrated thrombocytosis to be both a risk factor for VTE in inpatient populations,12 and to be associated with increased mortality risk in population-based cohorts.37 To our knowledge, the association between low platelet count and increased VTE risk is novel, with one possible explanation that clumping of platelets occurs with platelet activation and could cause an artificially low platelet count.
Our evaluation of VTE risk factors is in keeping with other less comprehensive previous studies, which have demonstrated the influence of obesity, fractures, smoking, BMI and medications including oral corticosteroids and oral contraceptives.10 38 We were able to explore these and other risk factors with adjustment for other patient characteristics. Results highlight an interesting absence of association with VTE for traditional cardiovascular disease risk factors including hypertension and hyperlipidaemia.
Strengths of our large, long-term population-based study include the comprehensive capture of VTE risk factors and patient characteristics, allowing interrogation not only of VTE risk across multiple diseases in adjusted analysis but also assessment of independent risk factors for VTE. Exposures and outcomes were defined using algorithms previously validated in primary care. Interpretation of coefficients for individual risk factors may be limited by the potential of confounding, and these estimates do not provide a causal interpretation.39 A further limitation of the study, similar to all studies using routine data, include the potential of unmeasured confounding and selection bias. Findings may not be generalisable to more ethnically diverse populations than the UK. Despite the use of validated algorithms to classify CD, UC, RA and the use of published guidance to define PsA, the lack of medical record review and use of clinical criteria to classify these IMIDs is a further limitation of the study, since diagnoses were recorded in primary care and may not have been made by specialists. When evaluating VTE risk factors, chance findings offer a potential explanation for differences in the groups with and without an IMID due to the number of associations tested. Family history of VTE is poorly captured in primary care data, and a resultant lack of power offers the most likely explanation for the observation that family history of VTE was not a significant risk factor in the IMID cohort. Similarly, this study will have systematically under captured biologic medication prescribing as, in the UK, these are prescribed by specialists and not captured in primary care. Secondary care data were not available to evaluate risk associated with surgery, an established major VTE risk factor. For analysis of time-updated platelet count, we used a complete-case approach, and for other missing information including BMI and ethnicity, we used the missing indicator variable method, as data are likely to be missing not at random meaning multiple imputation may lack validity.40
Our study is timely and of particular relevance in the context of the clinical interest in VTE in people with IMID. Our data provide an understanding of the contextual risk in IMID populations, and suggests considerable potential to update or augment existing VTE risk stratification decision aids such as the Wells Score with more refined multivariable prediction models incorporating routinely measured clinical patient characteristics. Another interesting direction for future research would be to use time-updated risk models to evaluate the temporal association between measures of IMID disease extent and severity, and by use of medication, in particular immunosuppressive treatment, and risk of VTE. This would provide further important information for clinicians responsible for monitoring patients with IMIDs in primary care.
In summary, VTE is more common in people with UC, CD, RA, and PsA compared with those without these diseases, highlighting the need for increased awareness among clinicians. Although associations do not have a causal interpretation, this study refines our understanding of classical VTE risk factors in people with an IMID compared with the wider population. Our data provide an initial platform for the risk assessment of individual patients with an IMID, and support active monitoring and strategies to mitigate VTE risk in people with an IMID.
What is already known about this subject?
Risk of venous thromboembolism (VTE) is increased in people with immune-mediated inflammatory diseases (IMIDs; ulcerative colitis, Crohn’s disease, rheumatoid arthritis and psoriatic arthritis) compared with the general population, but differences in VTE risk have not been systematically compared across these conditions.
The magnitude and relevance of VTE risk from traditional VTE risk factors (such as obesity, fractures, and use of specific medications) in IMIDs is unknown.
What does this study add?
In over 266 890 people, risk of VTE was increased to a similar degree in people with ulcerative colitis, Crohn’s disease and rheumatoid arthritis. For psoriatic arthritis, risk was not significantly increased, likely due to lack of statistical power.
Risk factors identified in people with IMIDs include male sex, overweight/obese BMI, smoking, fractures, use of corticosteroids and oral contraceptives, and abnormal platelet count.
How might this impact on clinical practice?
Knowledge of specific risk factors in people with immune-mediated inflammatory diseases can help identify those susceptible to developing VTE.
Patients and practices who are members of the Royal College of General Practitioners (RCGP) Research and Surveillance Centre (RSC) network, who allow their data to be shared for surveillance, research, quality improvement and education. The collaboration of primary care computerised medical record system providers EMIS, TPP, InPractice Systems and Apollo in facilitating the RCGP RSC data.
Contributors JG, KB, PI, KK, MN, RN, SdeL and MHB designed the study, supervised the data analysis, provided the interpretation of results, and contributed to the drafting and critical review of the manuscript. All authors approved the final draft. Medical writing and statistical support, funded by Pfizer, was provided by John Dennis (PhD), Andrew McGovern (MD) and Anita Lynam (PhD) [Momentum Data], with project management support from Filipa Ferreira (PhD) [University of Surrey]. As corresponding author, MHB attests that all listed authors meet authorship criteria and that no others meeting the criteria have been omitted. MHB is the guarantor and accepts full responsibility for the conduct of the study, had access to the data, and controlled the decision to publish.
Funding This study, and medical writing support, was funded by Pfizer. As study authors, representatives of Pfizer contributed to the design and conduct of the study, interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Competing interests All authors have completed the ICMJE uniform disclosure form at URL: www.icmje.org/coi_disclosure.pdf and declare: JG has received honoraria and/or sponsorships for conferences from AbbVie, Celgene, Janssen, Pfizer and UCB. KB has received honoraria from Tillots, Thermo Fisher Scientific, Boehringer Ingelheim, Pfizer, and Yakult. PI has received lecture fees from AbbVie, Celgene, Falk Pharma, Ferring, MSD, Janssen, Pfizer, Takeda, Tillotts, Sapphire Medical, Sandoz, Shire, and Warner Chilcott; financial support for research from MSD, Pfizer, and Takeda; advisory fees from AbbVie, Arena, Genentech, Gilead, Hospira, Janssen, Lilly, MSD, Pfizer, Pharmacosmos, Prometheus, Roche, Sandoz, Samsung Bioepis, Takeda, Topivert, VH2, Vifor Pharma, and Warner Chilcott. KK, MN and RN are employees of Pfizer. SdeL is Director of the RCGP RSC, he has received funding through his universities from Eli Lilly, GSK, Astra Zeneca, MSD, Sanofi, Seqirus, and Takeda. MHB has provided expert advice and received consultant fees and/or sponsorship for conference from AbbVie, Boehringer ingelheim, Eli Lilly, EMD Serono, Gilead, Pfizer, Roche, and Sanofi and has received research grants paid to her employer from Pfizer, Roche and UCB.
Patient consent for publication This research was done without patient involvement. Patients were not invited to comment on the study design and were not consulted to develop patient relevant outcomes or interpret the results. Patients were not invited to contribute to the writing or editing of this document for readability or accuracy.
Ethics approval Study approval was granted by the Research Committee of the RCGP RSC. The study did not meet the requirements for formal ethics board review as defined using the NHS Health Research Authority research decision tool (URL: http://www.hra-decisiontools.org.uk/research/).
Provenance and peer review Not commissioned; externally peer reviewed.
Data availability statement No additional data are available from the authors although Royal College of General Practitioners (RCGP) Research and Surveillance Centre (RSC) United Kingdom primary care data are available by application to the Research Committee of the RCGP RS.
Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.