{'fraction_non_alphanumeric': 0.03855503246253731, 'fraction_numerical': 0.022062217298993814, 'mean_word_length': 5.376299784186776, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 2, 'https://': 0, 'lorem ipsum': 0, 'www.': 2, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Meta-Analysis and Systematic Review of Coagulation Disbalances in COVID-19: 41 Studies and 17,601 Patients

Introduction Coagulation parameters are important determinants for COVID-19 infection. We conducted meta-analysis to assess the association between early hemostatic parameters and infection severity. Methods Electronic search was made for papers that addressed clinical characteristics of COVID-19 patients and disease severity. Results were filtered using exclusion and inclusion criteria and then pooled into a meta-analysis to estimate the standardized mean difference (SMD) with 95% confidence interval (CI) for D-dimers, fibrinogen, prothrombin time, platelet count (PLT), activated partial thromboplastin time. To explore the heterogeneity and robustness of our fundings, sensitivity and subgroup analyses were conducted. Publication bias was assessed with contour-enhanced funnel plots and Egger's test by linear regression. Coagulation parameters data from retrospective cohort study of 451 patients with COVID-19 at National Research Center for Cardiac Surgery were included in meta-analysis of published studies. Results Overall, 41 original studies (17,601 patients) on SARS-CoV-2 were included. For the two groups of patients, stratified by severity, we identified that D-dimers, fibrinogen, activated partial thromboplastin time, and prothrombin time were significantly higher in the severe group [SMD 0.6985 with 95%CI (0.5155; 0.8815); SMD 0.661 with 95%CI (0.3387; 0.9833); SMD 0.2683 with 95%CI (0.1357; 0.4009); SMD 0.284 with 95%CI (0.1472; 0.4208)]. In contrast, PLT was significantly lower in patients with more severe cases of COVID-19 [SMD −0.1684 with 95%CI (−0.2826; −0.0542)]. Neither the analysis by the leave-one-out method nor the influence diagnostic have identified studies that solely cause significant change in the effect size estimates. Subgroup analysis showed no significant difference between articles originated from different countries but revealed that severity assessment criteria might have influence over estimated effect sizes for platelets and D-dimers. Contour-enhanced funnel plots and the Egger's test for D-dimers and fibrinogen revealed significant asymmetry that might be a sign of publication bias. Conclusions The hemostatic laboratory parameters, with exception of platelets, are significantly elevated in patients with severe COVID-19. The two variables with strongest association to disease severity were D-dimers and fibrinogen levels. Future research should aim outside conventional coagulation tests and include analysis of clotting formation and platelet/platelet progenitors characteristics.

There is a need for early detection of elevated coagulation biomarkers to optimize risk stratification of patients with COVID-19. Importantly, multiple findings suggest that thromboembolic complications should be considered as a cause of clinical deterioration in severe COVID-19 infection and potentially beyond hospital discharge (27)(28)(29). However, due to the limitations of early pandemic research with small study sizes and highly heterogeneous datasets, it is still debated.
The aim of the current study was to evaluate the validity of stratification based on early coagulation parameters provided by different clinical studies worldwide. We also aimed to use subgroup analysis to explore whether differences in patient subpopulations partly explained heterogeneity in results. The coagulation parameters of our own cohort were included in a systematic meta-analysis of published studies. Furthermore, we also discuss a complex mechanism of hemostatic dysregulation that fuels the coagulation cycle.

MATERIALS AND METHODS
In this study, we aimed to identify the relationship between coagulation biomarkers taken at admission and severity of COVID-19 in adult (above 18) patients by estimating the effect size of five laboratory coagulation tests: D-dimers (DD), platelet count (PLT), fibrinogen (FIB), activated partial thromboplastin time (APTT), and prothrombin time (PT).

Search Strategy, Exclusion, and Inclusion Criteria
To broaden our search, we included Ovid, PubMed, Web of Sciences, and Google Scholar in the study for the period ending September 30th, 2021. They were thoroughly scanned using the following keywords: SARS-CoV-2, COVID-19, coagulation, severity, characteristics, features, D-dimers, platelets, fibrinogen, cohort, observational, retrospective. First, satisfying search results were exported to an Excel table, and duplicates were eliminated. Second, we conducted a preliminary review of retrieved articles' abstracts. Papers that focused on pregnant women, children, or specific age groups were avoided. Similarly, works that focused on the effect of medical treatment or patients with particular commodities were not included in the metaanalysis. We also excluded studies that compared mortality among groups of COVID-19 patients. While diagnostic and classification of COVID-19 severity follow official guidelines, the differentiation between whether COVID-19 caused deaths might be ambiguous. In order to avoid such bias, we only considered studies that grouped patients according to the severity of their symptoms. Then, each publication was checked according to the inclusion criteria: sample size, severity assessment criteria, presence of at least one coagulation marker of interest expressed as a continuous variable. To be included, a study should have had more than 15 patients and divided them according to the disease severity-WHO guidelines, intensive care unit (ICU) admission, disease aggravation, or need for oxygen therapy. Next, a list of biomarkers should have included at least some of the coagulation parameters: PLT, DD, FIB, APTT, PT, thrombin time (TT), Activated thrombin (AT), Factor VIII (FVIII), and VWF. Finally, the quality of selected articles was assessed with the NIH Quality Assessment Tool for Observational Cohort and Cross-Sectional Studies (30) (Supplementary Table 1). Questions 5, 6, 12, 13 were not applicable to selected studies due to their observational/retrospective nature; Question 8 was not applicable to studies with dichotomous exposure values, like ICU admission and need for oxygen therapy. Hence, studies with a quality score of <5 were excluded from the analysis. Search and quality check were conducted independently by two authors. Disagreements were resolved during a joint discussion with a third author.

Data Extraction
The primary goal was to retrieve data from each article about two groups of patients differing in disease severity by: group size, stratification criteria, and results of the laboratory tests. Apart from these, we kept a record of the country, hospital, and admission time to identify and exclude duplicate patients.

Participants and Data Collection
From a large cohort of 560 patients consecutively admitted from June 2020 to August 2021 due to coronavirus infection to the cardiology department of the National Research Center for Cardiac Surgery, Nur-Sultan, Kazakhstan, we selected 451 patients with a confirmed diagnosis of COVID-19, aged ≥18 years, and those with the onset of the disease ≤21 days, all of whom were hospitalized for at least 24 h with COVID-19. All demographic, clinical, laboratory data were extracted from the electronic records of the NRCSC. Laboratory tests included PLT, PT, international normalized ratio, APTT, FIB, DD. Blood samples were collected at admission by the clinical team. Clotting tests were performed according to standard methods. For coagulation tests (PT, international normalized ratio, APTT, FIB), we used AVATUBE R with sodium citrate, and it was measured on Sysmex CS-2500 (Sysmex, Japan). Blood for DD was collected in AVATUBE R with sodium citrate, and D-dimers were quantified on Cobas R 6000. PLT counts were measured on Sysmex XS-500 (Sysmex, Japan) with blood collection performed in AVATUBE R with K2 EDTA. The confirmation of SARS-CoV-2 infection was done by real-time quantitative reversetranscription polymerase chain reaction (RT-PCR) assay on nose/throat swab or sputum samples using CFX96 R Real-Time System (Bio-Rad Laboratories, Inc., USA).

Clinical Assessment
The COVID-19 population was initially divided into four groups according to disease severity following local guidelines "Diagnosis and treatment of Covid-19 in Adults" (http:// www.rcrz.kz/index.php/ru/2017-03-12-10-51-13/klinicheskieprotokoly) (Supplementary Data 1). The study protocol was approved by the medical ethics committee of the National Research Center for Cardiac Surgery, Nur-Sultan, Kazakhstan waiving the need for written informed consent due to the retrospective design.

Data Transformation
All data analysis was performed in the RStudio 1.4.1717 with R version 4.1.0 (Integrated Development Environment for R. RStudio, PBC, USA). In most articles, results were reported as the first, second (median), and third quartiles, reflecting that original data might be abnormally distributed. Since the authors of the primary studies do not report the underlying distributions, and there are no established methods to analyze the difference of medians, we performed data transformation according to the method described by McGrath et al. (31). The workflow takes advantage of the formulas proposed by Luo et al. (32) and Wan et al. (33) for estimation of mean and standard deviation, respectively. Both methods assume the input data to be normally distributed, so the package also includes the preceding transformation step by the Box-Cox method. In the scope of this meta-analysis, we compare two groups: more and less severe. Patients classified as severe, critical, admitted to ICU, with observed aggravation of a disease, or in need of oxygen therapy were assigned to the first group. The second group included mild, moderate, outpatients, general ward patients, recovered and stable patients (Supplementary Table 2). If, for example, a primary study had more than two groups-mild and moderate, severe and critical-then mild patients were combined with moderate, and severe were combined with critical.

Meta-Analysis
The "meta" and "metafor" packages were used for the metaanalysis. Heterogeneity was assessed using the I 2 , tau-squared, and prediction intervals. Following the recommendations of Veroniki et al. (34), Paule-Mandel (PM) tau-squared estimator and Q-Profile (QP) tau-squared confidence interval estimator were used in this meta-analysis. To account for different measurement scales, the standardized mean difference (SMD), or Hedges g was estimated. Before pooling point estimates, each study was assigned a weight based on inverse variance approach-smaller variance gives higher precision, thus yielding greater weight for a study. For datasets with high heterogeneity, the Knapp-Hartung adjusted random-effects model is implemented. In case of low heterogeneity, the fixedeffects model is applied.

Sensitivity and Subgroup Analysis
We also performed several sensitivity analyses using the "dmetar" package. First, we recalculated the effect size using the leaveone-out method. Second, basic outliers, studies whose confidence intervals do not overlap with the pooled effect confidence interval, were removed. Third, we ran the diagnostic for influential cases with Graphic Display of Study Heterogeneity (GOSH) plots and excluded them from the meta-analysis. To further explore the heterogeneity of the data, we conducted subgroup analyses: based on grouping criteria and location.

Publication Bias
We used the counter-enhanced Funnel plots and Egger's test by linear regression to assess the possibility for publication bias in selected studies.

Results of the Web Search
The stepwise process of selection is depicted in Figure 1. Initially, there were 4,056 papers that matched search parameters. After the removal of duplicates and studies that did not meet the exclusion criteria, only 195 unique studies remained. A more thorough review based on the inclusion criteria has proven 125 studies to be unfitting. In addition to 70 that were left, 15 studies were extracted from the references lists of previously published meta-analyses and systematic reviews. Finally, we performed a quality check and removed duplicate patients, which resulted in 40 studies being included in the meta-analysis. Besides original articles, we have included data from the National Research Center for Cardiac Surgery located in Nur-Sultan, Kazakhstan. This data set included data from 451 patients with acute COVID-19.

Description of Selected Studies
Overall, this meta-analysis comprises 17,601 patients from 40 research papers and the NRCSC. Of the total, 24 studies originated from China, while the rest included data from France, Germany, Italy, Mexico, Singapore, South Korea, Turkey, the UK, and the USA. Most of the articles reported the division of patients according to the WHO guidelines. Out of the rest, 12 articles made groups based on the ICU admission, three focused on disease aggravation, and three on the need for oxygen therapy. This information is depicted in Figure 2 in more detail.

Pooled Effect Sizes
We collected data for five coagulation parameters: PLT, DD, FIB, APTT, and PT. In all cases, the heterogeneity was extremely high (I 2 > 80%), so the Knapp-Hartung adjusted random-effects model was applied in all parts of the meta-analysis. DD, FIB, The sizes of red squares are proportional to weights of respective studies, which are adjusted random-effects weights: reciprocal of the sum of squares of study's standard error and the variance of the distribution of true effect sizes. Pooled effect size and its confidence intervals are shown by blue diamond at the bottom. Detailed output of the meta-analysis is given in Table 1.

Sensitivity and Subgroup Analysis
Neither the analysis by the leave-one-out method nor the influence diagnostic has identified studies that solely cause a significant change in the effect size estimates. Removal of several outliers based on the pooled and individual confidence intervals allowed significant decrease in heterogeneity, yet it did not affect the interpretation of the results. Adjusted effect sizes and heterogeneity measurements are summarized in Table 2. All influence diagnostic plots, such as Baujat and leave-one-out plots, Cook's distance, Covariance Ratio, etc., are provided in the Supplementary Figures 1-16 along with resulting Forest plots.
GOSH plots were built for all five models. K-means, DBSCAN, and Gaussian Mixture Model were implemented as clustering methods. This diagnostic was not applied on PLT and DD since they only have one cluster each (  Figures 6C-E). GOSH plots of the prothrombin time effect size identified four studies, three of which (36,39,40) pull effect size to the left (Figures 6F-H), and one (41) pulls effect size to the right ( Figure 6I). We have excluded these studies and recalculated the pooled effect size for each parameter. Results are summarized in Table 2.
To explore the possible sources of heterogeneity and test the robustness of our data, we have conducted subgroup analyses based on articles' country of origin and criteria for severity stratification. Results are summarized in Table 3. Subgroup analysis showed no significant difference between articles originated from different countries but revealed that severity assessment criteria might have an influence on overestimated effect sizes for PLT and DD. High heterogeneity levels seem to be persistent oven different subgroups. Being the only exception, values of APTT obtained from non-Chinese articles were less heterogenous-47.1% compared to 84.7% in studies originated in China.

Publication Bias
Contour-enhanced funnel plots for five coagulation parameters are depicted in Figure 7, and the results of Egger's test are summarized in Table 4. Based on the linear regression results, significant asymmetry is observed in published studies that included results for PLT, DD [Intercept (p-value): −1.297 (0.0336), 3.182 (0.0013)], and somewhat significant   . While it is impossible to predict whether asymmetry is coming solely from publication bias, contour-enhanced funnel plots for DD and FIB (Figures 6B,C) show that the number of papers with insignificant findings is very small (p > 0.1). This may advocate for the presence of some publication bias. In contrast, Figure 6A illustrates that even though the funnel plot for PLT has some asymmetry, there are plenty of publications with insignificant SMDs, thus decreasing the possibility of publication bias.
To further explore the robustness of our findings, we have conducted several additional tests. As reported by Zwetsloot et al. (42), plotting SMD against its standard error to assess the publication bias can result in distorted funnel plots, leading to false-positive conclusions. Since SMD is necessary for our analysis and cannot be substituted, we used a sample sizebased precision estimate instead of the standard error-1/ √ n. The resulting funnel plots (Supplementary Figure 16) showed no significant distortion. Similar corrections were made to Egger's test: since standardized mean difference and standard error are not independent, we can expect the test to yield some false-positive results. In our calculations, we used the corrected SE suggested by Pustejovsky and Rodgers (43). Results are given in the Supplementary Table 3. The only model that was affected by this correction was FIB (intercept 2.142 with a p-value of 0.0526 compared to the initial 1.872 and 0.0493), indicating that previously concluded asymmetry might have been a falsepositive. Interpretation of results for other parameters remained the same.

DISCUSSION
There is a need to identify early indicators of the risk of critical patients progressing with COVID-19 to ICU. Abnormal coagulation parameters, including DD and FIB levels, and changes in PLT and PT, are observed in many patients with COVID-19 infection at admission [rev. (8,13)]. However, as discussed by different authors, these indicators provide contradictory information regarding risk stratification and prognosis (44)(45)(46). Several previous meta-analysis reviews are researching abnormal coagulation and the severity of COVID-19. The majority of analyzed studies were conducted at the beginning of the pandemic and examined the Chinese population of patients [before April 2020- (14,(47)(48)(49); June 2020- (50)].
Our meta-analysis study found that DD and FIB levels, APTT, and prothrombin time were significantly higher in the severe group of COVID-19 patients. In contrast, PLT was significantly lower in these patients. There was no evidence of effects across different geographic subgroups, except for APTT values.

Platelet Changes
In this meta-analysis, we systematically analyze evidence for the utility of lowered PLT levels on admission in relation to the severity of infection. Platelets play a crucial role in the maintenance of hemostasis, and contribute to thromboinflammatory processes (51). At different stages of Frontiers in Cardiovascular Medicine | www.frontiersin.org viral infection, changes in platelet production or destruction may result in coagulation imbalances leading to pro-thrombotic events or in platelet dysfunction and bleeding risks (52) and increased activation of platelets (19). Thrombocytopenia (platelet count <150 × 10 9 /L) seems to represent an important marker of COVID-19 severity (8,14,18,53) and mortality outcome (14,18). The reasons for thrombocytopenia are multi-factorial may include early suppression of platelets progenitors' production and damage in bone marrow and lungs, hemophagocytosis and other reasons. Significant thrombocytopenia in COVID-19 is less frequent than it was reported for the largest SARS-1 cohort (45%) (54). SARS-1 induced platelet depletion was related to direct viral infection of megakaryocytes and hematopoietic progenitors [rev. by (4)]. However, almost half of patients developed thrombocytosis following initial thrombocytopenia during the SARS epidemics in 2003 (55). It is possible that opposite PLT changes may reflect different phases of COVID-19 infection (56)-acute and convalescent phases, though the explanation has not been established yet. The reactive thrombocytosis that may occur during the convalescent phase requires time observation, including a period preceding patient discharge. The increase of immature reticulated platelets reported recently in COVID-19 patients by different groups (57,58) may play a role in the limiting effects of anti-platelet therapies (59)(60)(61). Change in platelets parameters (increased mean volume) (62) and reactivity (63,64) are suggested to be associated with a severe COVID-19 infection (65). In addition to platelet changes, there are an increasing number of studies related to the dysfunction of megakaryocytes (platelets progenitors) in COVID-19 patients. Thus, elevated and abnormal megakaryocytes (MKs) were found in post-mortem and autopsy studies in the lungs (66)(67)(68)(69), heart (68), brain (70,71), bone marrow (72,73). Recently, the concept that megakaryopoiesis occurs in the lung parenchyma with significant contribution of the lungs in platelet biogenesis has emerged (74)(75)(76). There is a close association between platelet biogenesis and alveolarization of the lung (77,78), and hypoxia-induced thrombocytopenia associated with reduction of lung megakaryocytes and impaired generation of platelets in the lungs (79). The SARS-CoV-2 virus elicited up-regulation of megakaryocyte progenitors and elevated circulating megakaryocytes in blood in severe COVID-19 (80,81), hypothetically through infection of early megakaryocyte progenitors. Recently, Fortmann et al. (80) detected infected calprotectin+ spike+ megakaryocytes in the blood of patients with severe COVID-19. We hypothesize that viral RNA next may be transferred from megakaryocytes to platelets and circulating plasma EVs. The strong supporting argument is a recent finding of Brunet-Ratnasingam et al. (82) that plasma SARS-CoV-2 RNA is strongly associated with increased mortality. There is a current gap of knowledge regarding platelet non-canonical function and regeneration, and further research is crucial for our understanding of COVID-19 pathogenesis and stratification of increased thromboembolic risk.

DD and FIB Levels and Thromboembolic Complications
Concordant with other studies, we found DD and FIB levels significantly elevated in severe COVID-19 infection. However, DD levels results may suffer from some publication bias as revealed by meta-regression analysis. It may be related to the variability of cut-offs for varying DD measurement assays used in the different studies [(83); rev. (84)]. The usefulness of hemostasis tests, particularly, D-dimers, is limited by interlaboratory variability reflecting differences in methodologies, antibody sources and specificity, different detection methods, calibrators, and diagnostic thresholds (85,86). As confirmed by autopsy results, the elevated DD levels can be associated with fibrin deposits within pulmonary extravascular space and alveoli (67), but they may be non-specific to intravascular fibrin formation (87). High DD levels, together with elevated neutrophils, were reported to be predictive factors of pulmonary embolism in COVID-19 patients (88). DD and FIB levels are elevated in both COVID-19 and thromboembolism, and therefore as single tests are unspecific and unhelpful in the differentiation of these conditions. In part, it may be explained by the involvement in abnormal coagulation in COVID-19 of different systems, including endothelial cells, complement activation, and hypofibrinolysis resulting in changes undetected by routine tests (39,88). A significantly higher amount of neutrophil extracellular traps [and related markers such as cell-free DNA and myeloperoxidase-DNA (MPO-DNA)] was found when compared with non-COVID-19 thrombi (89). Neutrophil extracellular traps (NETs) markers, and cell-free DNA was most closely aligned with D-dimers, with inflammatory markers (C-reactive protein, lactate dehydrogenase) and more specific markers of NETs, MPO-DNA (90)(91)(92). The presence and high concentrations of NETs at COVID-19 patients' tracheal aspirates, blood samples (93,94), and different organs (73) are associated with COVID-19 severity. The formation of NETs in the lungs and bloodstream can be critically associated with thrombosis (95). Platelet activation and endothelial cell damage (96,97) in COVID-19 patients are also resulting in elevated levels of extracellular vesicles (EVs) (98-100) expressing tissue factor (101,102) and associated with severity of disease (99,101,103).

PT, APTT, FIB, and Clot Formation Parameters
PT has been included in the International Society on Thrombosis and Hemostasis (ISTH) criteria for diagnosing disseminated intravascular coagulation (DIC) (104), prolonged at COVID-19 with acute distress respiratory syndrome (ARDS) (10) and is significantly elevated in severe COVID-19 patients (105). Progressive prolongation of PT is considered a predictor of fatal outcome.
The first reports from China described a shortened APTT which is considered to be a predictive marker of hypercoagulation. Later studies reported prolonged APTT (45,106). A prolongation of APTT may indicate a deficiency of clotting factors (II, V, VIII, IX, X, XI, or XII) or the presence of inhibitor (heparin therapy). Another likely explanation is the presence of antiphospholipid antibodies (aPL) or lupus anticoagulant (LA) observed in 45% of COVID-19 patients in the Mulhouse series (107). In Bowles's study (106), most of the patients with prolonged APTT were positive for LA (91%). However, increasing of heparin therapy usage, especially in severe COVID-19 cases, may also lead to prolonged APTT.
Discordant reports exist concerning the fibrinogen levels and severity of COVID-19 disease. There are reports that COVID-19 is associated with high fibrinogen levels (16,108) and alternative findings (109). Some authors reported that severe COVID-19 infection induces fibrinolysis shutdown (26,45), which authors consider partly responsible for high DD and FIB levels. However, there are also reports of an increase in plasminantiplasmin complexes in severe COVID-19 and findings of mild consumption coagulopathy seen in a proportion of COVID-19 patients (110). Elevated tissue-type plasminogen activator was also found to be increased in the patients with severe COVID-19 (101).
Fibrinogen is a critical molecule that links coagulation, complement system, and inflammation, and at COVID-19, a significant association between FIB levels and interleukin-6 was described (16). Recent post-mortem studies of COVID-19 patients identified microvascular thrombi in multiple organs, including kidney (67), which may lead to increased DD even in the situation of reduced fibrinolysis. Depression of fibrinolysis was reported not only in ICU patients with severe infection and clinical signs of thromboembolism (26) but also in the general unselected ICU cohort (108,109). Viscoelastic testing ([ROTEM, tissue-type plasminogen activator (tPA) ROTEM, TEG] and evaluation of clot formation by Clot Waveworm analysis have demonstrated a hypercoagulable state, characterized by increased clot stiffness and severely impaired fibrinolysis (17,108,111,112). However, no consistent association between abnormal VET (viscoelastometric testing) pattern at COVID-19 and clinical outcome have been demonstrated [rev. (113)].
Elevated routine coagulation parameters taken at admission may help with the initial differentiation of severe and nonsevere COVID-19 patient groups. The hypothetical mechanisms of fibrinolysis shutdown implicate lowering of fibrinolysis factors (plasminogen) and elevated amounts of inhibitors of fibrinolysis (such as α2-antiplasmin and plasminogen activator inhibitor PAI-1).
Longitudinal observational studies of COVID-19 infection are particularly important for understanding infection dynamics and disease outcomes (114)(115)(116)(117)(118)(119). Hardy et al. (117) reported daily changes of relevant parameters of hemostasis of severe COVID-19 patients, including standard hemostatic tests (D-dimers, PT) as well as functional integrative tests for thrombin generation and fibrinolysis. Importantly, the day of hospital and/or ICU admission does not represent the same timepoint in the COVID-19 disease course that increases variability and complicates interpretation of studies results.
Although our findings were consistent across different geographic groups, they were tempered by significant heterogeneity. The DD, FIB, and PLT elevation trends speak about the necessity of investigating the potential of additional specific biomarkers for the stratification of COVID-19 patients. We suggest that thorough characterization of platelets and their progenitors (megakaryocytes), as well as clot formation testing could aid in the assessment of risk forthromboembolic events in COVID-19 patient groups.

LIMITATIONS OF THE STUDY
Our findings should be interpreted with some limitations. Firstly, most of the included studies were either retrospective and observational by design, thus prone to recall or misclassification bias. Second, such studies often imply recruitment by convenience sampling, in which case the representativeness of a study is questioned. In the context of this meta-analysis, several different endpoints were allowed: severity guidelines, ICU admissions, oxygen therapy requirement, and disease aggravation. In an attempt to overcome this issue, we have conducted a subgroup analysis. Third, uncontrolled variables might also pose a limitation-both the disease severity and coagulation parameters can be affected by several potential confounding variables, such as comorbidities, age, etc. Fourth, data transformation was applied to those studies that reported their results in quartiles or maximum and minimum values. This could potentially introduce some degree of error in the pooled estimate. Fifth, the random-effects model was used due to the high heterogeneity of research. Although sensitivity and subgroup analyses were conducted, some residual heterogeneity may affect the interpretation of results. Finally, the underlying mechanisms of coagulation disbalances at COVID-19 still need to be investigated.

CONCLUSIONS
Our findings support using a combination of coagulation parameters for risk stratification of patients with COVID-19 infection at the time of admission, in particular DD and FB levels. However, additional evaluation using clotting formation methods or combination with specific biomarkers of production and activation of platelet and platelet progenitors is required. This information may help physicians triage patients with thromboembolism. An important consideration is an intermediate and long follow-up of COVID-19 patients.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

ETHICS STATEMENT
The studies involving human participants were reviewed and approved by Ethics Committee of the National Research Center for Cardiac Surgery (NCJSC), Nur-Sultan, Kazakhstan. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

AUTHOR CONTRIBUTIONS
PL: conceptualization, methodology, formal analysis, investigation, data curation, and writing original draft. GI, ZS, and AK: methodology, data curation, formal analysis, and investigation. ATT, MS, AD, MA, and MB: methodology, formal analysis, investigation, data curation of NRCSC cohort, and writing a draft. EP and AT: methodology, analysis, and funding. NB: conceptualization, methodology, formal analysis, investigation, data curation, finding, and writing a draft. MB and NB: supervision. All authors read and edit a final draft. All authors contributed to the article and approved the submitted version.

FUNDING
NB and AT were supported by FDCRGP SSH2020028 and NB was funded by OPCRP2020018 grants from Nazarbayev University.

ACKNOWLEDGMENTS
We are very grateful to Sadyk Khamitov and Aliya Sailybayeva from the Research Department of NCJSC National Research Center for Cardiac Surgery for their excellent organizational help and to Alexander Kovaldgi from 2nd Mathematical Liceum for helpful discussion.