- 1Division of Rheumatology, Department of Medicine, University of Washington, Seattle, WA, United States
- 2Rheumatology Section, Veterans Affairs (VA) Puget Sound Healthcare System, Seattle, WA, United States
- 3Department of Biostatistics, University of Washington, Seattle, WA, United States
- 4Advanced Cancer Research Group, Kirkland, WA, United States
- 5Public Health Sciences and Clinical Research Divisions, Fred Hutchinson Cancer Center, Seattle, WA, United States
- 6Schools of Medicine, Public Health and Pharmacy, University of Washington, Seattle, WA, United States
- 7Department of Rheumatology, Martina Hansens Hospital, Oslo, Norway
- 8Department of Clinical Sciences Lund, Rheumatology, Lund University, Lund, Sweden
Giant cell arteritis (GCA) that affects older patients is an independent risk factor for thromboembolic events. The objective of this study was to identify predictive factors for thromboembolic events in patients with GCA and develop quantitative predictive tools (prognostic nomograms) for pulmonary embolism (PE) and deep venous thrombosis (DVT). A total of 13,029 patients with a GCA diagnosis were included in this retrospective study. We investigated potential predictors of PE and DVT using univariable and multivariable Cox regression models. Nomograms were then constructed based on the results of our Cox models. We also assessed the accuracy and predictive ability of our models by using calibration curves and cross-validation concordance index. Age, inpatient status at the time of initial diagnosis of GCA, number of admissions before diagnosis of GCA, and Charlson comorbidity index were each found to be independent predictive factors of thromboembolic events. Prognostic nomograms were then prepared based on these predictors with promising prognostic ability. The probability of developing thromboembolic events over an observation period of 5 years was estimated by with time-to-event analysis using the method of Kaplan and Meier, after stratifying patients based on predicted risk. The concordance index of the time-to-event analysis for both PE and DVT was > 0.61, indicating a good predictive performance. The proposed nomograms, based on specific predictive factors, can accurately estimate the probability of developing PE or DVT among patients with GCA.
Introduction
Giant cell arteritis (GCA), a form of large vessel vasculitis, is a chronic systemic inflammatory disease of unknown etiology that primarily affects older patients (1). GCA can be complicated by thromboembolic events such as pulmonary embolism (PE) and deep venous thrombosis (DVT) (2).
In recent years, terms such as immunothrombosis and thromboinflammation have been emerged. In the absence of infection, uncontrolled inflammation may induce a pro-coagulant state leading to the development of thrombotic events (3, 4). Recent studies demonstrate that patients with GCA have a higher incidence of PE and DVT within the first year of their diagnosis (5) and an increased risk of cerebrovascular events (6), especially in the presence of atherosclerotic risk factors at the time of diagnosis of GCA (7). Atherosclerotic risk factors are also known to be significantly associated with venous thrombotic events (8). In another report, an increased risk of thromboembolism was related to the period of higher disease activity at the onset of GCA (9). However, in a different study the incidence of both thromboembolic and cerebrovascular events did not differ between biopsy confirmed GCA and non-GCA patients (10). Interestingly, hospitalizations in patients with GCA were found to increase the risk of thromboembolic events compared to hospitalizations in the general population (11).
It is also known that thrombocytosis is associated with GCA (12). Leukocyte and platelet activation may contribute to vessel inflammation and possibly to thromboembolic events in patients with GCA (13). In a small study, 13 out of 19 patients with biopsy proven GCA had evidence of thrombocytosis several months to a year before the diagnosis and treatment of GCA (14). Medications that can be “thrombogenic” include diuretics that may cause hemoconcentration and lead to increased blood viscosity via changes in fibrinogen concentrations and coagulation time (15). Of note a recent study demonstrated that diuretic use might be weakly associated with a lower risk of GCA (16).
Predictive factors for thromboembolic events in patients with GCA have not been carefully investigated. To the best of our knowledge, there are no models that provide quantitative prediction of thromboembolic events among patients diagnosed with GCA. One of the main objectives of this study was to propose specific predictive factors that will help physicians to better assess patient risk for thromboembolic events and develop quantitative tools (such as nomograms) for prediction in patients with a diagnosis of GCA among a Veteran-based population. The secondary objective of our study was to perform a time-to-event analysis to estimate the probability of a thromboembolic event over an observation period of 5 years.
Methods
Data sources and study population
In this institutional review board (IRB)-approved retrospective study (#01854) patients with GCA were included between January 1st 1999 and January 13th 2022. The Veterans Affairs Clinical Data Warehouse was used for data extraction and collection. Thromboembolic events were identified, by diagnostic codes of hospital admissions and outpatient visits within the Veterans Health Administration system as previously conducted by our group (5).
Study design and data collection
Our study cohort included patients with a diagnosis of GCA without a history or prior diagnosis of GCA over the prior 2 years. The study cohort included patients: (a) aged ≥50 years; and (b) classified with at least one GCA clinical modification code from the International Classification of Diseases Ninth or Tenth Revision (ICD-9-CM or ICD-10-CM) (Supplementary Table 1) by a rheumatologist or other health provider.
The following data were collected: patient age, race/ethnicity, sex, baseline body mass index (BMI), smoking status prior to the diagnosis of GCA, use of aspirin, anticoagulants (heparin, warfarin, low molecular weight heparin, rivaroxaban, apixaban) or diuretics (thiazide, loop, potassium sparing, and osmotic diuretics, carbonic anhydrase inhibitors) within 5 years prior to diagnosis of GCA. The Charlson comorbidity index, which is frequently used as a measure of general health based on 19 medical conditions (17, 18), was also generated using data up to 5 years prior to the index date. Thrombocytosis was defined as platelet count ≥ 450x109/L (19) between a two-week period before and after the date of diagnosis of GCA. Data about the initial diagnosis of GCA in the outpatient or inpatient setting, and number of hospital admissions within 5-years for any cause before diagnosis of GCA, were also recorded. These clinical and laboratory covariates were then used as potential predictive factors.
Study outcomes of interest
The main outcomes of interest were the incidence of PE or DVT after the diagnosis of GCA. We defined PE and DVT events using ICD‐9 and ICD-10 procedure codes (Supplementary Table 1).
Statistical analysis
Baseline demographic characteristics were detailed based on the stratification with or without thromboembolic events after the diagnosis of GCA. We reported mean ± standard deviation for quantitative variables. For categorical variables counts and proportions (based on all non-missing observations) are presented. Covariates specified a priori age, race, sex, BMI, smoking status, Charlson comorbidity index, diagnosis of GCA in the outpatient versus inpatient setting, number of hospital admissions for any reason before diagnosis of GCA, use of anticoagulation, aspirin and diuretics as well as hypertension within 5 years prior to the diagnosis of GCA, and thrombocytosis to be assessed as potential predictors in our study.
Univariable and multivariable Cox proportional hazards regression models were then applied to explore potential predictors of PE and DVT events for patients with GCA. The censoring events in the time-to-event analysis were: 1) death, 2) end of study (January 13th 2022), and 3) development of malignancies. The baseline covariates we used for these models are reported in Table 1. We report the hazard ratio (HR) with 95% confidence interval (CI) for each potential predictive factor, together with the associated p-value in Tables 2, 3. For all our analyses, p-values < 0.05, were considered statistically significant.
Table 1 Baseline demographic characteristics of GCA patients without and with thromboembolic events.
Using the R package rms, nomograms were created for each type of thromboembolic events with covariates selected by Cox models. Calibration curves were then plotted to evaluate the calibration ability of our nomograms. Calibration was also based on cross validation.
Stratified Kaplan-Meier estimates of the outcomes are also presented in our study, where the grouping is based on the predicted risk score of multivariable Cox models. Subjects having a risk score higher (or lower) than the sample median were classified as high (or low) risk respectively. Two well- separated Kaplan-Meier plots would indicate a better predictive ability of the Cox models. The log-rank test was used to assess comparisons. Moreover, to quantify the predictive power of the multivariable Cox model, we also report the testing concordance index, which is calculated using 5-fold sample splitting (20). All statistical analyses were done in R (version 3.6.3, http://www.r-project.org/).
Results
Baseline demographic characteristics of our study population
A total of 13,029 patients with a diagnosis of GCA were included in our analysis. Among them, 626 were diagnosed with thromboembolic events (PE and/or DVT) after the diagnosis of GCA, whereas 12,403 were not. Specifically, 261 patients developed PE and 440 patients DVT, whereas 75 patients developed both PE and DVT after the diagnosis of GCA. Our stratified summary of baseline covariates is presented in Table 1. The majority patients were white male. The mean age of patients within our two study groups was above 72 years old.
The initial diagnosis of GCA was placed at outpatient visits in more than 80% of the patients and the mean BMI was approximately 29 for both study groups with and without thromboembolic events. The 5-year Charlson comorbidity index was slightly higher in the group with thromboembolic events reflecting a worse state of health, compared to the group without thromboembolic events. Less than one third of the patients with GCA in both groups were recorded as smokers. Admission data in Table 1 shows that the mean number of hospital admissions for any cause was slightly higher in patients with GCA with thromboembolic events compared to patients without.
Only 2% of the patients in both groups were on anticoagulation treatment within 5 years prior to diagnosis of GCA. Interestingly, a higher percentage of patients with GCA that developed thromboembolic events had thrombocytosis, compared to patients without thromboembolic events. Similar percentages of both groups were on aspirin or diuretic therapy within 5 years prior to the diagnosis of GCA. The frequency of hypertension, which is an independent risk factor for thromboembolism (21), was similar in patients within both groups.
Risk factors of thromboembolic events among patients with GCA
We explored potential predictors of developing PE and DVT during an observation period of five years after the diagnosis of GCA using Cox regression models (Tables 2, 3). In the univariable Cox regression analysis, for the outcome of PE, the HR was 1.13 (p<0.0001) for Charlson comorbidity index, 1.05 (p=0.006) for hospital admissions, 0.65 (p=0.03) for diagnosis of GCA in the outpatient setting (left column, Table 2). In the multivariable analysis, only Charlson index (HR 1.13, p<0.0001) was significantly associated with the outcome of PE (right column, Table 2).
On the other hand, age (HR 1.01, p=0.02), Charlson comorbidity index (HR 1.12, p<0.0001), diagnosis of GCA in the inpatient setting, and hospital admissions before the diagnosis of GCA (HR 1.06, p<0.0001), were strong independent predictors of DVT in our multivariable model (right column, Table 3). Results were different for the univariable analysis, where thrombocytosis was also found to be a predictor of DVT (HR 3.09, p=0.0009), as depicted in the left column of Table 3. In conclusion, Charlson comorbidity index was a statistically significant independent risk factor for both PE and DVT. Other covariates may also have predictive ability as we show later in this study.
Evaluation of predictive models
We present the stratified Kaplan-Meier curves for outcomes of PE and DVT in Figures 1, 2 respectively. As demonstrated in these figures, patients with GCA, corresponding to a higher risk score (according to the multivariable Cox model), have distinct time-to-event incidence of PE or DVT events compared to those with a lower score. The x-axis represents time in months after the diagnosis of GCA and the y-axis depicts the probability of patients with GCA surviving without experiencing PE or DVT. The blue and red lines represent Kaplan-Meier plots of the low and high-risk groups of GCA for PE or DVT. For example, at 40 months after the diagnosis of GCA, the predicted probability of experiencing PE is almost 1% for the low-risk group, whereas the patients in the predicted high-risk group have approximately 2.5-fold greater chance of developing PE (Figure 1).
Figure 1 Kaplan-Meier plot for the time to PE from the time of diagnosis of GCA stratified by predicted risk (blue and red lines) respectively.
Figure 2 Kaplan-Meier plot for the time to DVT from the time of diagnosis of GCA stratified by predicted risk (blue and red lines) respectively.
For PE the testing concordance index (C-index) is 0.62 ± 0.02, and for DVT it is 0.61 ± 0.01, indicating a good predictive discrimination. The standard error was estimated using bootstrap. These results suggest that if we know the baseline covariates (as listed in Tables 2, 3) of two patients with GCA there is around 60% chance that we can correctly predict which one will develop PE or DVT events earlier, based on our multivariable regression models.
Nomograms of thromboembolic events for patients with GCA
The univariable and multivariable Cox regression analyses indicate that age, Charlson comorbidity index, diagnosis of GCA in the inpatient setting, hospital admissions for any cause before the diagnosis of GCA, and thrombocytosis are independent factors that are associated with the occurrence of thromboembolic events in patients with GCA. However, these models do not give a simple or direct quantity that physicians can use readily in practice. For this purpose, we created nomograms based on some of these predictive variables to estimate the 1-year probability of developing PE or DVT in patients with GCA (Figures 3, 4). We also present the calibration curves for the nomograms to evaluate their accuracy (Figures 5, 6). For both PE and DVT outcomes, the predicted risk is in agreement with the apparent proportion for low to moderate risk patients, and this statistical regime includes most of the patients. However, for the PE event the nomogram would overestimate the risk for patients having a higher risk >0.03 that may be due to the linearity assumption of logistic regression. On the other hand, for the DVT event the nomogram works well for the whole range of risk.
Figure 3 Nomogram for predicting the 1-year probability of PE in patients with GCA. NumAdmissions, number of admissions; Outpt, outpatient; Inpt, inpatient.
Figure 4 Nomogram for predicting the 1-year probability of DVT in patients with GCA. NumAdmissions, number of admissions; Outpt, outpatient; Inpt, inpatient.
Figure 5 Calibration plot of the nomogram for predicting 1-year overall probability for PE in patients with GCA. Nomogram-predicted probability for PE is plotted on the x-axis while actual probability is plotted on the y-axis; The grey line represents the ideal line. .
Figure 6 Calibration plot of the nomogram for predicting 1-year overall probability for DVT in patients with GCA. Nomogram-predicted probability for DVT is plotted on the x-axis while actual probability is plotted on the y-axis; The grey line represents the ideal line.
To read the nomogram, for each variable draw a straight line upward to the axis labeled as Points at the same vertical position, and then calculate the total points. Afterwards, draw a straight line downward to the axis labeled as 1-year probability PE or DVT at the same vertical position to find the patient’s likelihood of 1-year probability of PE or DVT. For instance, a 65 year-old patient with Charlson comorbidity index of 6 is admitted to the hospital 15 times before the diagnosis of GCA. The points for each of the variables are 18 (65 year-old), 30 (Charlson score of 6), and 20 (15 hospital admissions before diagnosis of GCA); the patient has a total score of 68 points, and the 1-year probability of developing DVT turns out to be 0.02 (Figure 4).
Discussion
In this retrospective study, we propose four independent predictive factors for thromboembolic events among veterans diagnosed with GCA that include age, Charlson comorbidity index, initial diagnosis of GCA in the inpatient setting, and number of hospital admissions before the diagnosis of GCA in our multivariable model. We then constructed a clinical prediction model based on these predictive factors to predict the 1-year thromboembolic event probability of patients with GCA at any time after the diagnosis of GCA, which could guide clinicians in their practice and advance the health of patients with GCA by creating thromboembolic risk stratification models. Future prospective multi-center studies are required to validate the accuracy and effectiveness of our prediction models/nomograms.
Most of the predictive factors for thromboembolic events that are proposed in this study are consistent with previous studies. Similar to our finding, Charlson comorbidity index has been identified as a predictor of VTE among high-risk cancer patients undergoing chemotherapy in the ambulatory setting (22). In the Khorana risk score, increased platelet count was also found to be an important predictor for VTE among patients with cancer (23, 24). Although thrombocytosis was significantly associated with the outcome of DVT, it was not included in our predictive models because of missing data, but its validity as a potential risk factor for VTE can be assessed in prospective multicenter studies. Older age has also been recognized as a predictor for VTE (25). In another study there was moderate evidence that there is a probable association between risk of VTE and increasing age (≥60 years old) in admitted medical patients (26).
Our finding that hospital admissions is a strong independent predictor of thromboembolic events among patients with GCA is supported by existing literature that shows that hospitalized GCA patients are at a greater risk of VTE (11). Our overall findings lay the foundation for future studies investigating the complexity of interactions between thrombo-inflammation, aging, multi-morbidities, and thromboembolic risk in patients with GCA, and/or other inflammatory conditions.
One of the limitations of our study is that potential confounders such as, sepsis, trauma or surgery, were unable to be taken into consideration, and residual confounding is a possibility. Another limitation is possible selection bias because of the retrospective nature of this study. It is also uncertain whether our results can be applied to the general population as the military population is represented heavily by men. It is well known in the literature that males with GCA are at higher risk of developing VTE compared to females with GCA (27). Finally, histologic evidence of temporal arteritis, representing the gold standard for the diagnosis of GCA, was not available for patients limiting diagnostic accuracy.
There is also a possibility of misclassification, as we did not include in the definition of GCA, glucocorticoid therapy for at least 6 months to avoid survival bias. For example, patients with GCA that developed fatal PE may not have survived for up to 6 months to receive corticosteroid therapy for the full period. Thus, glucocorticoid therapy was not part of our inclusion criteria, as our main goal was to identify high-risk patients at the time of GCA diagnosis. Additionally, inflammatory markers, clinical subtypes of GCA, advanced imaging studies such as magnetic resonance angiography and fluorodeoxyglucose (FDG)-positron emission tomography, and disease activity status, were not recorded for this study.
Our study has several strengths. This is the first large study that proposes independent predictors and risk model for thromboembolic events among patients with GCA diagnosis in a US military population. Our findings may help physicians to accurately calculate with readily available clinical risk factors the probability of developing PE or DVT among veterans with GCA over a period of one year, integrating specific predictive variables that can easily be extracted from the medical chart for most hospitals and academic centers, by using our nomogram model as a reference.
The results of our study may also assist physicians to stratify patients with GCA at low and high risk for VTE and predict the probability of a thromboembolic event over an observation period of five years in a male population with GCA. Additionally, this stratification strategy may be useful for developing strategies for screening of lower extremities with Doppler ultrasonography among high-risk GCA patients for VTE in an attempt to identify asymptomatic events of DVT (28, 29).
In conclusion, use of thrombosis predictive models to predict the risk of VTE may help with risk stratification of patients at increased risk according to the probability of a thromboembolic event. Thus, there is an unmet need for further validation of the clinical applicability of these predictive models by large multicenter, prospective clinical studies.
Data availability statement
The datasets presented in this article are not readily available due to ethical/privacy restrictions. Requests to access the datasets should be directed to BN, Bernard.Ng@va.gov.
Ethics statement
The Research and Development Committee from the Department of Veterans Affairs approved the study (MIRB#01854). Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author contributions
DM, TZ, NK, GL, AD, PS, and BN contributed to the conception and design of this study. DM, TZ, and BN participated in data collection and data analysis. Only DM and BN had access to raw data. All authors contributed to data interpretation, critically reviewed and revised the manuscript and approved the final submission.
Funding
DM is supported by the T32 NIH grant from Nutrition, Obesity and Atherosclerosis (#5T32HL007028-44). The funder body was not involved in the study design, collection, analysis and interpretation of the data, the writing of this article or the decision to submit it for publication.
Conflict of interest
DM received Advisory Board fees from ChemoCentryx and a research grant from Pfizer US pharmaceuticals. NK reports consulting fees from Beyond Spring, Bristol Myers Squibb, G1 Therapeutics, Sandoz, Seattle Genetics, Janssen, Pfizer, and Spectrum and research funding from Amgen to her institution (family member). GL reports research grant support from Amgen to the Fred Hutchinson Cancer Center and consulting fees from Beyond Spring, G1 Therapeutics, Partner Therapeutics, Samsung Bioepis, Merck, Jazz, TEVA, Squibb, Sandoz, Seattle Genetics, and Fresenius Kabi.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022.997347/full#supplementary-material
References
1. Pugh D, Karabayas M, Basu N, Cid MC, Goel R, Goodyear CS, et al. Large Vessel vasculitis. Nat Rev Dis Primers (2022) 7(1):93. doi: 10.1038/s41572-021-00327-5
2. Unizony S, Lu N, Tomasson G, Zhang Y, Merkel PA, Stone JH, et al. Temporal trends of venous thromboembolism risk before and after diagnosis of giant cell arteritis. Arthritis Rheumatol (2017) 69(1):176–84. doi: 10.1002/art.3984
3. Stark K, Massberg S. Interplay between inflammation and thrombosis in cardiovascular pathology. Nat Rev Cardiol (2021) 18(9):666–82. doi: 10.1038/s41569-021-00552-1
4. Aksu K, Donmez A, Keser G. Inflammation-induced thrombosis: mechanisms, disease associations and management. Curr Pharm Des (2012) 18(11):1478–93. doi: 10.2174/138161212799504731
5. Michailidou D, Zhang T, Stamatis P, Ng B. Risk of thromboembolism in giant cell arteritis and/or polymyalgia rheumatica: A veterans health administration population-based study. J Intern Med (2022) 291(5):665–75. doi: 10.1111/joim.13446
6. Ungprasert P, Wijarnpreecha K, Koster MJ, Thongprayoon C, Warrington K. Cerebrovascular accident in patients with giant cell arteritis: A systematic review and meta-analysis of cohort studies. Semin Arthritis Rheumatol (2016) 46(3):361–6. doi: 10.1016/j.semarthrit.2016.07.005
7. Gonzalez-Gay MA, Piñeiro A, Gomez-Gigirey A, Garcia-Porrua C, Pego-Reigosa R, Dierssen-Sotos T, et al. Influence of traditional risk factors of atherosclerosis in the development of severe ischemic complications in giant cell arteritis. Med (Baltimore) (2004) 83(6):342–7. doi: 10.1097/01.md.0000145369.25558.b5
8. Mi Y, Yan S, Lu Y, Liang Y, Li C. Venous thromboembolism has the same risk factors as atherosclerosis: A PRISMA-compliant systemic review and meta-analysis. Med (Baltimore) (2016) 95(32):e4495. doi: 10.1097/MD.0000000000004495
9. Azhar SS, Tang RA, Dorotheo EU. Giant cell arteritis: diagnosing and treating inflammatory disease in older adults. Geriatrics (2005) 60(8):26–30.
10. Lo Gullo A, Koster MJ, Crowson CS, Makol A, Ytterberg SR, Saitta A, et al. Venous thromboembolism and cerebrovascular events in patients with giant cell arteritis: A populationbased retrospective cohort study. PLoS One (2016) 11(2):e0149579. doi: 10.1371/journal.pone.0149579
11. Unizony S, Menendez ME, Rastalsky N, Stone JH. Inpatient complications in patients with giant cell arteritis: Decreased mortality and increased risk of thromboembolism, delirium and adrenal insufficiency. Rheumatol (Oxford) (2015) 54(8):1360–8. doi: 10.1093/rheumatology/keu483
12. Price N, Clearkin LG. Thrombocytosis and giant cell arteritis (let-ter). Lancet (1994) 343:672. doi: 10.1016/s0140-6736(94)92665-4
13. Maugeri N, Baldini M, Rovere-Querini P, Maseri A, Sabbadini MG, Manfredi AA. Leukocyte and platelet activation in patients with giant cell arteritis and polymyalgia rheumatica: A clue to thromboembolic risks? Autoimmunity (2009) 42(4):386–8. doi: 10.1080/08916930902832629
14. Lincoff NS, Erlich PD, Brass LS. Thrombocytosis in temporal arteritis rising platelet counts: a red flag for giant cell arteritis. J Neuroophthalmol (2000) 20(2):67–72.
15. Marvel RJ, Shullenberger WA. Thromboembolic phenomena associated with rapid diuresis in the treatment of congestive heart failure. Am Heart J (1951) 42:194. doi: 10.1016/0002-8703(51)90240-2
16. Turkiewicz A, Stamatis P, Mohammad AJ. Cardiovascular drug treatment, statins and biopsy-confirmed giant cell arteritis: A population-based case-control study. RMD Open (2020) 6(2):e001285. doi: 10.1136/rmdopen-2020-001285
17. Glasheen WP, Cordier T, Gumpina R, Haugh G, Davis J, Renda A, et al. Charlson comorbidity index: ICD-9 update and ICD-10 translation. Am Health Drug Benefits (2019) 12(4):188–97.
18. Bannay A, Chaignot C, Blotière PO, Basson M, Weill A, Ricordeau P, et al. The best use of the charlson comorbidity index with electronic health care database to predict mortality. Med Care (2016) 54(2):188–94. doi: 10.1097/MLR.0000000000000471
19. Bleeker JS, Hogan WJ. Thrombocytosis:diagnostic evaluation, thrombotic risk stratification, and risk-based management strategies. Thrombosis (2011) 2011:536062. doi: 10.1155/2011/536062
20. Hastie T, Tibshirani R, Friedman J. The elements of statistical learning: Data mining, inference, and prediction Vol. 2 New York: springer (2009).
21. Lowe GD. Common risk factors for both arterial and venous thrombosis. Br J Haematol (2008) 140(5):488–95. doi: 10.1111/j.1365-2141.2007.06973.x
22. Khorana AA, Dalal M, Lin J, Connolly GC. Incidence and predictors of venous thromboembolism (VTE) among ambulatory high-risk cancer patients undergoing chemotherapy in the United States. Cancer (2013) 119(3):648–55. doi: 10.1002/cncr.27772
23. Khorana AA, Kuderer NM, Culakova E, Lyman GH, Francis CW. Development and validation of a predictive model for chemotherapy-associated thrombosis. Blood (2008) 111(10):4902–7. doi: 10.1182/blood-2007-10-116327
24. Khorana AA, Soff GA, Kakkar AK, Vadhan-Raj S, Riess H, Wun T, et al. Rivaroxaban for thromboprophylaxis in high risk ambulatory patients with cancer. N Engl J Med (2019) 380(8):720–8. doi: 10.1056/NEJMoa1814630
25. Brand JS, Hedayati E, Bhoo-Pathy N, Bergh J, Hall P, Humphreys K, et al. Time-dependent risk and predictors of venous thromboembolism in breast cancer patients: A population-based cohort study. Cancer (2017) 123(3):468–75. doi: 10.1002/cncr.30364
26. Darzi AJ, Karam SG, Charide R, Etxeandia-Ikobaltzeta I, Cushman M, Gould MK, et al. Prognostic factors for VTE and bleeding in hospitalized medical patients: a systematic review and meta-analysis. Blood (2020) 135(20):1788–810. doi: 10.1182/blood.2019003603
27. Aviña-Zubieta JA, Bhole VM, Amiri N, Sayre EC, Choi HK. The risk of deep venous thrombosis and pulmonary embolism in giant cell arteritis: A general population-based study. Ann Rheum Dis (2016) 75(1):148–54. doi: 10.1136/annrheumdis-2014-205665
28. Khorana AA, Vadhan-Raj S, Kuderer NM, Wun T, Liebman H, Soff G, et al. Rivaroxaban for preventing venous thromboembolism in high-risk ambulatory patients with cancer: Rationale and design of the CASSINI trial. rationale and design of the CASSINI trial. Thromb Haemost (2017) 117(11):2135–45. doi: 10.1160/TH17-03-0171
Keywords: giant cell arteritis, thromboinflammation, pulmonary embolism, deep venous thrombosis, thromboembolic events, nomograms, thrombocytosis, predictors
Citation: Michailidou D, Zhang T, Kuderer NM, Lyman GH, Diamantopoulos AP, Stamatis P and Ng B (2022) Predictive models for thromboembolic events in giant cell arteritis: A US veterans health administration population-based study. Front. Immunol. 13:997347. doi: 10.3389/fimmu.2022.997347
Received: 18 July 2022; Accepted: 17 August 2022;
Published: 09 November 2022.
Edited by:
Angelo A. Manfredi, Vita-Salute San Raffaele University, ItalyReviewed by:
Georgios Sogkas, Hannover Medical School, GermanyAlberto Lo Gullo, Garibaldi Hospital, Italy
Copyright © 2022 Michailidou, Zhang, Kuderer, Lyman, Diamantopoulos, Stamatis and Ng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Despina Michailidou, despina1@uw.edu
†These authors have contributed equally to this work