Thrombosis is a frequently underdiagnosed condition associated with high mortality in neglected cases. Many factors, including geoheliophysical and biochemical ones, are responsible for thrombosis modulation. Routine investigations may sometimes be inconsistent and, thus, unreliable in a clinical setting.
Data were collected from patients treated in the Department of Vascular Surgery at the ‘Timofei Moșneaga’ Republican Clinical Hospital, Chișinău, Republic of Moldova. A total of 1,865 patients were initially included in the study. After applying rigorous inclusion and exclusion criteria, 263 eligible patients were identified, and their complete blood counts and biochemical reports were retrospectively analyzed.
The analysis revealed increased mean values for absolute polymorphonuclear neutrophils, absolute monocytes, erythrocyte sedimentation rate (ESR), and glucose. The median values of these indicators, except for absolute polymorphonuclear neutrophils and ESR in female patients, were also elevated above normal ranges. Significant Pearson and Spearman correlations were identified among the analyzed indicators, and a binary logistic regression model was constructed using the most statistically significant variables.
Usual mathematical models that outline thrombosis consider deep vein thrombosis without a sustainable arterial assessment. The sensitivity of our model is lower than that of the D-dimer, while the specificity is almost the same. Platelets and clotting tests are well-known, reliable indicators; however, novel contemporary augmentations to these may, in turn, increase the predictive capability of our model if applied. This study has its limitations due to the lack of variance in the variance inflation factors (VIF), preventing the evaluation of multicollinearity among the included biomarkers.
The mathematical model developed in this study shows potential for further clinical application; however, additional research, validation, and the incorporation of non-biochemical indicators may be necessary to enhance its predictive accuracy.
There are no mathematical models that can accurately predict thrombosis risk in a clinical setting; thus, the development of one is a mainstay task.
There are certain biomarkers that may be used to develop a mathematical model capable of predicting thrombosis risk.
The development of a new mathematical model that may become more reliable and sustainable over time.
Thrombosis is defined as the formation of a blood clot that can obstruct the normal permeability of the circulatory system [1]. The mortality burden associated with thromboembolic events is substantial (25% fatal outcomes) [2].
Differences in thrombosis rates across regions are thought to be influenced by variations in plasma rheology and clotting factors, such as fibrinogen, and factor VII, which could potentially overshadow the impact of seasonality [3]. The interactions of the endocrine system and its array of pathologies are of utmost importance when considering the etiology of thrombosis. Conditions like hyperthyroidism, hypercortisolism (endogenous or exogenous), growth hormone excess, hyperprolactinemia, and pregnancy, which involve direct or indirect involvement of the hormonal system, are known to have a higher incidence of thrombosis [4].
In the past, erythrocytes were considered to play a passive role in thrombotic events. However, they are now recognized as active contributors. Erythrocytes can enhance blood viscosity, interact with vessel walls, and influence the spatial organization of platelets during aggregation [5]. This aggregation, a critical process in thrombosis, is driven by elevated levels of adenosine diphosphate (ADP) [6]. Leukocytes, especially monocytes, are renowned for their capacity to promote clot formation in an environment that has the components required for thrombosis to emerge [7]. Serum bilirubin has remarkable antioxidative and cytoprotective effects [8], which in turn warrants that its low concentrations will lead to increased venous thromboembolism events in an acute respiratory infection setting [9-10].
Routine hemoleukogram, which measures parameters such as red blood cell count, white blood cell count, and hemoglobin levels, remains a valuable tool in the assessment of deep vein thrombosis (DVT). These routine measurements are essential for identifying patients at risk of pre-thrombotic or thrombotic conditions [11]. While advanced biomarkers, including D-dimers, selectins, microparticles, and inflammatory cytokines, demonstrate high accuracy in thrombosis evaluation, their utilization is often restricted to private healthcare settings due to their high cost. This financial barrier poses significant challenges for many healthcare systems globally, limiting the accessibility of these biomarkers in routine clinical practice [12].
Current mathematical models for assessing thrombosis risk, such as the Wells score, demonstrate limited performance. Their efficiency varies significantly, ranging from 0.5% to 23.4%, depending on the expertise of the specialist applying the model [13]. Additionally, these models perform poorly in critically ill patients, further limiting their reliability in high-risk populations [14].
The aim of this study was to develop a mathematical model for thrombosis diagnosis using serum biomarkers, particularly in patients lacking imaging data.
Study population. The study was conducted in the Department of Vascular Surgery at the Timofei Moșneaga Republican Clinical Hospital (RCH) in Chișinău, Republic of Moldova. In 2019, a total of 1,865 patients were admitted to the department for various vascular conditions, including atherosclerosis, atherothrombosis, thrombosis, and thromboembolism. The department primarily relies on a combination of biochemical laboratory investigations and imaging studies as the cornerstone of its clinical research and diagnostic methodologies.
Inclusion and exclusion criteria. To minimize systemic bias, reduce the risk of incomplete data inputs, and ensure consistency in data collection, stringent inclusion and exclusion criteria were established for the study population. Only patients who underwent a complete blood count (CBC) and biochemical investigation during the inpatient stage (on the first day of hospital admission) and tested negative for methicillin-resistant Staphylococcus (MRS) antibodies were included. A comprehensive range of comorbidities and conditions were considered exclusion criteria to ensure the reliability and objectivity of the biochemical results. The excluded conditions included hepatic and splenic disorders, oncological and infectious conditions, hematological disorders, autoimmune and dermatological conditions, vascular and cardiac conditions, gastrointestinal and endocrine disorders, and additional exclusions such as trauma, open wounds, systemic diseases requiring artificial grafts or stoma implants, post-thrombotic or post-inflammatory syndromes, chronic obstructive pulmonary disease (COPD), and patients receiving anticoagulant therapy. The exclusion criteria were carefully designed to reduce confounding variables and ensure high-quality biochemical data. It should be mentioned that the controls were confirmed to be patients who did not have thrombosis, using the same investigations as the patients with confirmed thrombosis (imaging and laboratory methods).
Statistical Data. A total of 263 patients included in the final study cohort were systematically documented and organized using Microsoft® Excel® 2013 (15.0.4569.1504). The data were collected retrospectively. The dataset encompassed a broad range of demographic, clinical, and laboratory parameters, including: demographic information (age, gender); clinical features (affected limb region: upper or lower limbs, type of thrombosis: arterial or venous, and season of occurrence); hematological parameters (white blood cell count [WBC], red blood cell count [RBC], hemoglobin [HGB], hematocrit [HCT], mean corpuscular volume [MCV], mean corpuscular hemoglobin [MCH], mean corpuscular hemoglobin concentration [MCHC], platelet count [PLT], red cell distribution width – coefficient of variation [RDW-CV], and differential counts [absolute and relative values] of neutrophils, lymphocytes, monocytes, eosinophils, and basophils); inflammatory markers (erythrocyte sedimentation rate [ESR]; biochemical parameters (urea, aspartate aminotransferase [ASAT], alanine aminotransferase [ALAT], creatinine, total bilirubin, conjugated bilirubin, unconjugated bilirubin, glucose, albumin, and total protein); coagulation profiles (Quick prothrombin, international normalized ratio [INR], fibrinogen levels, activated partial thromboplastin time [aPTT], and thrombin time [TT]. This comprehensive collection of variables ensured a robust dataset for the subsequent analysis, enabling a thorough investigation of potential associations between the data.
Statistical Methods. The data were analyzed using IBM SPSS Statistics (version 26.0.0.0). The Shapiro-Wilk and Kolmogorov-Smirnov tests were applied to assess the normality of the data distribution. For normally distributed data, a one-sample t-test was performed with a t-value of 1.97 and a 95% confidence interval (CI). For non-normally distributed data, the one-sample Wilcoxon signed-rank test was utilized. Descriptive statistics, including mean values, standard deviations, medians, and minimum and maximum values, were calculated to summarize the data. Pearson and Spearman correlation coefficients were determined to evaluate relationships between the studied variables. Binary logistic regression analysis was conducted to identify, using mathematical methods, biochemical markers with potential predictive value for thrombosis events. An algorithm was developed using the statistically significant biomarkers identified in the analysis. Sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) were calculated.
The mean age of the patients was 59.21±13.50 years while the median was 62 years, ranging in the interval of 23-90 years. Based on gender distribution, the study included 164 male patients (62.4%) and 99 female patients (37.6%). Among these, 74 patients (28.1%) presented during the winter season, 57 (21.7%) during spring, 72 (27.4%) during autumn, and 60 (22.8%) during summer. Of the total participants, 143 patients had thrombosis, while 120 served as controls.
In 125 cases (87.41%), the lower limbs were affected, whereas the upper limbs were affected in 18 cases (12.59%). Regarding vascular involvement, the arterial system was implicated in 63 cases (44.06%), while the venous system was affected in 80 cases (55.94%).
The Shapiro-Wilk and Kolmogorov-Smirnov tests demonstrated that the data followed a normal distribution, except for hematocrit, relative polymorphonucleated neutrophils, and relative lymphocytes.
The one-sample t-test indicated that all parameters, except for relative eosinophils and conjugated bilirubin, were statistically significant, with n = 263, t = 1.97, and a 95% confidence interval (CI). The one-sample Wilcoxon signed-rank test was applied to relative eosinophils, where the sample median (MDN = 1.40) differed from the hypothesized population median (2.5), n = 263, Z = -7.818, 95% CI, r = 0.482. Similarly, the test for conjugated bilirubin demonstrated a sample median (MDN = 0) differing from the hypothesized population median (3.0), n = 263, Z = -8.432, 95% CI, r = 0.519, confirming their statistical significance.
Mean values for absolute polymorphonucleated neutrophils, absolute monocytes, erythrocyte sedimentation rate (ESR) in both male and female patients, and glucose were elevated above normal values (Table 1). Additionally, the median values for absolute monocytes, ESR in male patients, and glucose were also above normal values (Table 1).
The strongest Pearson and Spearman correlations (Table 2) are detailed. It is important to note that these values do not necessarily reflect the patient’s condition during thrombosis events or in controls. Instead, they aim to identify potential bias factors in the diagnostic process. While such elevations are commonly associated with pathological states, they do not appear to directly correlate with the severity of the condition.
Table 1. Indicators above normal statistical values | |||||
Indicator | Mean±SD | Median value | Minimum | Maximum | Normal range |
NEUT# | 5.35±2.43 | 4.88 | 1.45 | 20.11 | 2-5 x 103/µL |
MONO# | 0.69±0.52 | 0.63 | 0.17 | 7.7 | |
Table 2. Pearson and Spearman correlations between the biochemical markers | |||
Pearson correlations | Value | Spearman correlations | Value |
NEUT% x LYMPH% | -0.950** | NEUT% x LYMPH% | -0.952** |
NEUT# x LYMPH% | -0.644** | QP x INR | -0.923** |
QP x INR | -0.617** | NEUT# x LYMPH% | -0.671** |
Considering the statistical significance of the b-values derived from the binary logistic regression model and the predictive capability for thrombosis events, as indicated by a Nagelkerke R² value of 0.203, a mathematical model has been developed. This model includes a constant and the following statistically significant variables: platelets (PLT), mean platelet volume (MPV), and thrombin time (TT). The following mathematical model may be applied: b = 4.854 – 0.009 x PLT – 0.461 x MPV + 0.096 x TT. Once the b-values are determined, the thrombosis chance (TC) is calculated using the following formula: TC=ⅇ^b/(1+ⅇ^b ). The specificity of this model was 55.8%, while its sensitivity was 76.9%. The negative predictive value (NPV) was 67%, and the positive predictive value (PPV) was 67.48%.
Previous models for thrombosis evaluation were primarily developed to assess deep vein thrombosis (affecting the upper or lower limbs) and pulmonary embolism, without adequately addressing arterial thrombosis [15-17]. In this research, a new model has been proposed that integrates both arterial and venous thrombosis, providing a comprehensive framework for the empirical management of patients. The proposed model relies on widely available and inexpensive biomarkers, making it particularly suitable for resource-limited settings where advanced imaging or specialized tests are not accessible. However, its specificity (55.8%) limits its standalone utility in clinical practice.
The widely used screening test of D-dimers has a high incidence of false-positive results due to its low specificity. As a result, its application may introduce a significant range of systematic errors, despite its clinical utility [18]. Our score is generally less sensitive than D-dimers; however, given their similar specificity, it could serve as an additional confirmatory tool for patients who cannot undergo this test.
Elevated absolute polymorphonucleated basophil values have been associated with increased mortality in patients with coronary artery disease (CAD), suggesting their potential role as a pro-thrombotic marker [19].
Platelets, the primary circulating elements responsible for clot formation [20], are expected to fluctuate during thrombosis events, reflecting their critical role in the coagulation process.
Although platelets are the key players in primary hemostasis, they can also contribute to secondary hemostasis through mediators found in their granules, including alpha, dense, and lysosomal granules. Additionally, cytosolic factors can sustain the coagulation cascade, further enhancing biochemical interactions at this level [21]. These molecules may serve as potential predictive biomarkers for future algorithm development [22-23].
Mean platelet volume (MPV) has been shown to reduce thrombosis risk in multiple regression analysis, consistent with the findings of this study [24]. However, this contrasts with classical high-evidence studies that identify elevated mean platelet volume as a risk factor for thrombosis [25]. MPV has been found to have low significance in predicting thrombosis in COVID-19 patients [26], while demonstrating good sensitivity and specificity in assessing thrombosis recurrence in conditions such as antiphospholipid syndrome (APS) [27]. However, in a retrospective study, MPV was found to be less effective compared to D-dimers [28].
Thrombin time (TT) reflects the ability of fibrinogen to convert into fibrin. An elevated TT may indicate fibrinogen deficiency, dysfunction, or disruption of other coagulation factors [29]. A study conducted in Japan states that TT, when enhanced by clot waveform analysis (CWA), may provide valuable insights into hemostatic abnormalities in patients without deficiencies in other clotting components [30]. This novel global coagulation assessment tool may also mitigate the influence of certain drugs due to its biophysical nature [31-32]. However, both activated partial thromboplastin time (aPTT) and TT are known for their inherent variability and limited diagnostic utility in thrombotic events. Individually assessed biomarkers are often deemed unreliable in thrombosis diagnosis [33]. A significant portion of false-negative results may be attributed to the depletion of key thrombogenic elements–platelets, fibrinogen/fibrin, protein C, and protein S–which directly impact Quick prothrombin (PT) and aPTT [34]. Notably, a widely recognized scoring system requires platelet count, D-dimer levels, PT, and fibrinogen values to predict the onset of disseminated intravascular coagulation (DIC) [35].
It can be hypothesized that enhancing current standard methods for thrombosis assessment may improve the specificity and sensitivity of the proposed mathematical model. However, this study did not evaluate variance inflation factors (VIF) to assess multicollinearity among the included biomarkers.
This study proposes a simple and cost-effective model for thrombosis risk prediction based on platelets, MPV, and thrombin time. While the model shows promise for clinical application, further validation, refinement, and the integration of advanced biomarkers are needed to enhance its predictive accuracy and clinical applicability.
None declared.
DC designed the study, collected, and analyzed the data. IT critically revised the manuscript and analyzed the data. EP critically evaluated the results and assessed their applicability. OA designed the protocol for patient’ eligibility. EC critically revised the manuscript.
No approval was required for this study.
The authors declare no external funding.
Not commissioned, externally peer-reviewed.
Dan Croitoru – https://orcid.org/0000-0002-8915-0157
Iurie Trohin – https://orcid.org/0009-0001-8680-5402
Ecaterina Pavlovschi – https://orcid.org/0000-0003-0385-4805
Oleg Arnaut – https://orcid.org/0000-0002-5483-8672
Eugen Cerevan – https://orcid.org/0000-0002-3221-7584
0.09-0.6 x 103/µL
ESRtotal | 16.87±28.89 | 12 | 1 | 445 | 2-10 mm/h |
ESRmale gender | 17.89±35.45 | 12 | 1 | 445 | 2-10 mm/h |
ESRfemale gender | 15.19±11.72 | 12 | 2 | 65 | 3-12 mm/h |
Glucose | 5.98±2.07 | 5.6 | 1.47 | 20.71 | 3.5-5.5 mM/L |
Note: NEUT# – absolute polymorphonucleated neutrophils; MONO# – absolute monocytes; ESR – erythrocyte sedimentation rate; SD – standard deviation. |
MCH x RDW-CV
-0.540** |
NEUT# x BASO% |
-0.542** |
NEUT% x LYMPH# | -0.476** | NEUT% x BASO% | -0.516** |
BASO# x BASO% | 0.766* | HGB x HCT | 0.831** |
RBC x HCT | 0.803** | RBC x HCT | 0.838** |
EO# x EO% | 0.883** | EO# x EO% | 0.908** |
TotalBil x UnconjBil | 0.909** | TotalBil x UnconjBil | 0.919** |
WBC x NEUT# | 0.942** | WBC x NEUT# | 0.935** |
Note: * Statistical significance at p < 0.05; ** Statistical significance at p<0.01; NEUT% – relative polymorphonucleated neutrophils; LYMPH% – relative lymphocytes; BASO% – relative polymorphonucleated basophils; EO% – relative polymorphonucleated eosinophils; NEUT# – absolute polymorphonucleated neutrophils; LYMPH# – absolute lymphocytes; BASO# – absolute polymorphonucleated basophils; EO# – absolute polymorphonucleated eosinophils; RBC – red blood cells; WBC – white blood cells; HGB – hemoglobin; HCT – hematocrit; QP – Quick prothrombin; INR – international normalized ratio; MCH – mean corpuscular hemoglobin; RDW-CV – red cell distribution width (coefficient of variation); TotalBil – total bilirubin; UnconjBil – unconjugated bilirubin |