Skip to main content

Predictors and outcomes of acute pulmonary embolism in COVID-19; insights from US National COVID cohort collaborative

Abstract

Objectives

To investigate whether COVID-19 patients with pulmonary embolism had higher mortality and assess the utility of d-dimer in predicting acute pulmonary embolism.

Patients and methods

Using the National Collaborative COVID-19 retrospective cohort, a cohort of hospitalized COVID-19 patients was studied to compare 90-day mortality and intubation outcomes in patients with and without pulmonary embolism in a multivariable cox regression analysis. The secondary measured outcomes in 1:4 propensity score-matched analysis included length of stay, chest pain incidence, heart rate, history of pulmonary embolism or DVT, and admission laboratory parameters.

Results

Among 31,500 hospitalized COVID-19 patients, 1117 (3.5%) patients were diagnosed with acute pulmonary embolism. Patients with acute pulmonary embolism were noted to have higher mortality (23.6% vs.12.8%; adjusted Hazard Ratio (aHR) = 1.36, 95% CI [1.20–1.55]), and intubation rates (17.6% vs. 9.3%, aHR = 1.38[1.18–1.61]). Pulmonary embolism patients had higher admission D-dimer FEU (Odds Ratio(OR) = 1.13; 95%CI [1.1–1.15]). As the d-dimer value increased, the specificity, positive predictive value, and accuracy of the test increased; however, sensitivity decreased (AUC 0.70). At cut-off d-dimer FEU 1.8 mcg/ml, the test had clinical utility (accuracy 70%) in predicting pulmonary embolism. Patients with acute pulmonary embolism had a higher incidence of chest pain and history of pulmonary embolism or deep vein thrombosis.

Conclusions

Acute pulmonary embolism is associated with worse mortality and morbidity outcomes in COVID-19. We present d-dimer as a predictive risk tool in the form of a clinical calculator for the diagnosis of acute pulmonary embolism in COVID-19.

Introduction

As of Feb 2023, over 6 million deaths have been attributed to COVID-19 worldwide [1]. COVID-19 associated coagulopathy confers an increased risk of microvascular and macrovascular thrombosis [2]. Histological examination of pulmonary vessels in COVID-19 is distinguished from other viruses by widespread micro-thrombosis [3]. In the meta-analysis, high incidence rates of pulmonary embolism, 16.5% (95% CI: 11.6–22.9), have been noted in COVID-19 [4]. Relatively small studies have reported worse outcomes associated with PE in COVID-19; however, to our knowledge, no clear association of increased mortality associated with PE in COVID-19 has yet been reported [5, 6]. We wished to assess whether COVID-19 patients with diagnosed PE have increased mortality when compared to COVID-19 patients without PE.

The screening of acute PE in COVID-19 patients remains an important clinical question to date. Our traditional clinical predictors of screening for acute PE such as tachycardia, dyspnea, or hypoxia are confounded by the clinical presentation of COVID-19 pneumonia. The conventional scoring systems, such as Well’s criteria, have been noted to be non-discriminatory (AUC 0.54) in predicting the risk of acute PE in COVID-19 patients [7, 8]. Previous studies have pointed towards the possible role of d-dimer in the assessment of COVID-19 patients for screening of PE; however, no clear prediction scheme has been devised [5, 9, 10]. d-dimer is generated as a result of the lysis of blood clots as a part of the normal healing process. d-dimer is released when the cross-linked fibrin undergoes lysis by the plasmin-mediated degradation process [11]. d-dimer would be released in the blood during the breakdown process of a blood clot, as would be the case with pulmonary embolism. Using a national database, the National Collaborative COVID-19 cohort (N3C), we intended to develop an easy-to-use clinical calculator using d-dimer for predicting the risk of acute PE in hospitalized COVID-19 patients [12]. We also studied the commonly used clinical laboratory parameters to formulate the pathogenesis of acute PE in COVID-19.

Methods

Data source

The National Collaborative COVID-19 cohort (N3C) is the largest US database registry for COVID-19 patients and their representative controls [13, 14]. The N3C Data Enclave is approved under the authority of the NIH Institutional Review Board. An Institutional data-use agreement between the University of Kentucky and N3C data Enclave provided access to the de-identified database. N3C harmonizes its dataset based on four heterogeneous sources (ACT Network, TriNetX, OHDSI, PCORNet) and provides a unified platform within the N3C enclave for statistical analysis. The N3C cohort represents a rich, highly granular, ethnically, and geographically diverse dataset from clinical sites throughout the US (including Southeast, Midatlantic, and Midwest USA). The N3C cohort includes patients with any encounter after Jan 1, 2020, based on one of a set of SARS-CoV-2 laboratory tests or diagnostic codes. Further details about the N3C cohort are listed in the data source Additional file 1.

Measures, definition, and outcomes

The N3C database was accessed on Jan 7, 2021 (the time period of encounters from Jan 1, 2020, to Jan 7, 2021), which provided data from 29 clinical sites. We performed a retrospective observational cohort analysis using N3C, including individuals aged ≥ 18 years, with a positive SARS-CoV-2 PCR test within 2 weeks of hospitalization (Fig. 1). Based on the diagnosis of acute PE during the admission encounter, hospitalized COVID-19 patients were divided into PE and non-PE groups. Hospitalized COVID-19 patients are not universally screened with CT pulmonary embolism in the clinical setting. Hence, the small possibility of undiagnosed pulmonary embolism exists in this unscreened population—designated as non-PE in our study, based on the absence of acute pulmonary embolism diagnosis in the database. Baseline characteristics included demographics (age, sex, and race), comorbidities (atherosclerosis, congestive heart failure, diabetes, hypertension, atrial fibrillation, chronic kidney injury, acute kidney injury, chronic obstructive pulmonary disease, obstructive sleep apnea, hypothyroidism, obesity, anemia, nicotine dependence, malignancy, and sepsis) and inpatient medications use (steroids, remdesivir, tocilizumab, hydroxychloroquine). The continuous variables, if normally distributed, were reported as means ± SD or as medians and interquartile ranges, if not normally distributed. The primary outcomes were 90-day mortality rate and 90-day intubation rate from the time of hospitalization. Secondary analyzed parameters included length of stay, chest pain incidence, heart rate, history of deep vein thrombosis or PE, and comparison of admission laboratory parameters (ferritin, C-reactive protein [CRP], d-dimer fibrinogen equivalent units [FEU], and d-dimer units [DDU], fibrinogen, lactate dehydrogenase [LDH], and lymphocyte count). The secondary analyzed parameters included a history of deep vein thrombosis or PE, as we wanted to evaluate whether it is an associated risk factor for PE in COVID-19. The data accessed on Jan 7, 2021, was used for both primary and secondary outcomes analysis. In the clinical setting, two different d-dimer assays (d-dimer FEU and d-dimer DDU) are widely utilized, depending on the lab preference of the clinical site. The Fibrinogen Equivalent Unit (FEU) reporting of d-dimer levels is based on the molecular weight of fibrinogen (340 kDa) and the d-Dimer Unit (DDU) reporting of d-dimer levels is based on the molecular weight of d-dimer (195 kDa), which is about half that of fibrinogen. We performed a separate analysis for each assay.

Fig. 1
figure 1

Flow chart describing the schema of the study population. Primary outcomes included mortality and intubation analysis. Secondary outcomes included length of stay, chest pain incidence, heart rate, history of deep vein thrombosis or PE, and comparison of admission laboratory parameters (ferritin, C-reactive protein [CRP], d-dimer fibrinogen equivalent units [FEU], and d-dimer units [DDU], fibrinogen, lactate dehydrogenase [LDH], and lymphocyte count). Separate propensity score matched analysis was performed for each lab parameter

Statistical analysis

For the primary outcomes study, a multivariable time-dependent cox proportional hazard regression model was adjusted for demographics, comorbidities, and medication use as described earlier, to compare the mortality and intubation outcomes in the acute PE and non-PE hospitalized COVID-19 patient groups [15]. Patients were followed up until the time of an event or censored at 90 days for mortality or intubation outcomes. For the intubation outcomes analysis, in the PE group, patients who had PE diagnosed after the day of intubation were excluded. Overall survival and hazard for intubation were evaluated with the log-rank test. The Kaplan Meier and hazard curves were constructed for mortality and intubation outcomes, respectively. The covariates in the cox-proportional hazard model were carefully chosen based on existing evidence of their association with COVID-19 or all-cause mortality [16]. Further details about the assumption of proportionality of hazards in cox regression are listed in the Additional file 1 biostatistical methods. For the secondary outcomes study—which required comparison of absolute values post-matching—a separate adjusted 1:4 propensity score matching with the logistic regression model, greedy nearest neighbor, without replacement and caliper of 0.2 (Additional file 1: methods) was used to construct PE and non-PE groups (Fig. 1). The propensity score matching was preferentially used for the analysis of secondary parameters so that post-matching, the secondary parameters can be visually compared which can signify clinical significance (Fig. 4) and ROC curves can be constructed for d-dimer on the matched data (Additional file 1: Figs. S4, S5). 1:4 propensity score matching was performed for a larger number of non-PE patients than PE patients in order to reduce data wasting and hence, improve study power [17,18,19]. The covariates in the propensity score-matched model were either related to the outcome due to their association with COVID-19 or treatment (PE in COVID-19) [19, 20]. The covariates in the propensity score-matched analysis were also demographics (age, sex, and race), comorbidities (atherosclerosis, congestive heart failure, diabetes, hypertension, atrial fibrillation, chronic kidney injury, acute kidney injury, chronic obstructive pulmonary disease, obstructive sleep apnea, hypothyroidism, obesity, anemia, nicotine dependence, malignancy, and sepsis) and inpatient medications use (steroids, remdesivir, tocilizumab, hydroxychloroquine), similar to the covariates used in the primary outcomes study in cox-regression analysis. The standardized mean deviations were calculated to compare the balance between the two matched groups (Additional file 1: Fig. S1). They were then analyzed to compare the length of stay, chest pain incidence, heart rate (first reading recorded), and deep vein thrombosis or PE history. The propensity score matched results were analyzed via univariate logistic regression to calculate odds ratio. Since the admission labs were not present in all the patients, a separate adjusted propensity score-matched analysis was performed for each lab parameter (Additional file 1: Table S1). ROC curves were constructed to determine the characteristics of admission d-dimer FEU and d-dimer DDU in predicting acute PE in hospitalized COVID-19 patients. The sensitivity, specificity, positive predictive value, and accuracy of d-Dimer FEU and d-Dimer DDU to predict pulmonary embolism at different threshold levels are recorded in Additional file 1: Tables S2 and S3, which were subsequently exported to construct an easy-to-use calculator which was published online by our group [21]. The COVID-19 d-dimer calculator thus provides sensitivity, specificity, positive predictive value and accuracy of d-dimer FEU and d-dimer DDU to predict PE for various threshold levels [21]. The ROC optimal cut off value was selected based on the accuracy (accuracy ≥ 70%), as has been suggested previously by Korevaar et al. [22] Missing data (age 1.25%, race 12.9%, gender 0.5%) was imputed with the multivariate imputation by chained equations (MICE) throughout the data analysis. Significance was defined as p < 0.05 when using a two-tailed test. The odds ratio for categorical variables was calculated using Altman’s method on MedCalc software 20.013 [23, 24]. Odds ratio for continuous variables was calculated via univariate logistic regression. Further details about the statistical methods are provided in the Additional file 1. The analysis was performed within the Palantir Foundry, hosted within the cloud-based N3C enclave using Python 3.6 (Python Software Foundation) and R 3.5.1(R project for statistical computing).

Results

As of Jan 7, 2021, 390,681 patients were identified as COVID-19 positive based on PCR testing. Out of these, 31,500 COVID-19 adult patients were hospitalized. 1,117(3.5%) of these patients were diagnosed with acute PE during their admission. Primary propensity score matching yielded 1,117 PE patients and 4,468 non-PE patients (Fig. 1). 128 patients had PE diagnosed after the day of intubation and were excluded from the intubation outcomes analysis. Acute PE was diagnosed in 302 (9.7%) of 3,116 intubated patients (including those who were diagnosed post-intubation).

Baseline characteristics

Acute PE diagnosis in COVID-19 patients was recorded at median 0, interquartile range (IQR) 0–3 days post-admission. In our univariable analysis of COVID-19 patients, older patients (age 62 vs. 57 years; odds ratio (OR) 1.014, 95% CI [1.011–1.018]), and males (58% vs. 51%; OR 1.35 [1.20–1.52]) were more likely to have acute PE. A higher incidence of acute PE was noted in African Americans (OR 1.34 [1.18–1.52]) and Asian Americans (OR 1.45 [1.04–2.2]). Patients with acute PE had a higher prevalence of comorbidities, notably atrial fibrillation, hypertension, diabetes, congestive heart failure, obstructive sleep apnea, obesity, anemia, chronic kidney disease, malignancy. Higher incidence of DVTs (9% vs. 0.2%; OR 46.16 [33.65–63.3]), cerebrovascular diseases (5.8% vs. 3.5%; OR 1.66 [1.28–2.18]), acute encephalopathy (11.2% vs. 6%; OR 1.98 [1.63–2.39]), acute kidney injury (42.2% vs. 27.8%; OR 1.89 [1.67–2.14]) and sepsis (34.5 vs. 22.1%; OR 1.85 [1.63–2.1]) was also noted in acute PE patients. Patients with acute PE were more likely to receive treatments such as steroids, remdesivir, and tocilizumab (Table 1).

Table 1 Baseline characteristics of patients with and without acute pulmonary embolism

Primary outcomes

COVID-19 patients with acute PE had higher 90-day mortality (264 (23.6%) vs. 3918(12.8%), adjusted hazard ratio (aHR) = 1.36 [1.20–1.55] and 90-day intubation rates (174 (17.6%) vs. 2874 (9.3%), aHR = 1.38 [1.18–1.61] when compared to non-PE COVID-19 patients (p < 0.001 for both by log-rank test) (Figs. 2, 3). The unadjusted hazard ratio (HR) for mortality outcome was 1.92 [1.7–2.18] and 1.94 [1.66–2.26] for intubation outcome (p < 0.001 for both). Time from hospitalization to intubation was recorded at a median of 1 [IQR 0–4] days for both groups (OR 0.99[0.98–1.01]). For PE patients, who were intubated, the time difference between PE diagnosis and intubation since hospitalization was noted at 0 [IQR − 12 to 1.75] days (Additional file 1: Fig. S2). Hospital length of stay for patients who died was 14 [IQR 6.5–26] days for PE patients vs. 11 [IQR 5–20] days for non-PE patients (OR 0.99 [0.99–1]).

Fig. 2
figure 2

Kaplan Meier estimates of 90-day Survival in the COVID-19 patients with and without acute pulmonary embolism. Acute pulmonary embolism is associated with increased mortality in COVID-19 patients

Fig. 3
figure 3

Hazard curve of mechanical ventilation in the COVID-19 patients with and without Acute Pulmonary Embolism. Acute pulmonary embolism is associated with increased mechanical ventilation rates in COVID-19 patients

Secondary parameters

Hospital length of stay was longer for acute PE patients (median 9 [IQR 5–18] vs. 6 [3–11] days; OR 1.02 [1.02–1.03]. A history of DVT or PE (8.2% vs. 1.7%; OR 5.11 [3.75–6.98]) and incidence of chest pain (8.3% vs. 5.5%; OR 1.55 [1.21–1.99]) were more common in patients with acute PE (Table 2).

Table 2 Outcomes of COVID-19 patients with and without acute pulmonary embolism

For the laboratory/monitoring parameters, PE COVID-19 patients had higher CRP (OR 1.001 [1.0009–1.002]) and LDH values (1.0003 [1.00001–1.0006]); however, no appreciable clinical difference was noted, as the values were mostly overlapping (Additional file 1: Fig. S3). The two groups did not significantly differ in ferritin (OR 1 [0.99–1]) or lymphocyte count (OR 1.04 [0.99–1.09]). The PE patients had decreased fibrinogen (OR 0.99 [0.998–0.99]) compared to the non-PE patients. Heart rate was slightly increased in PE patients (OR 1.009 [1.004–1.013]).

The admission d-dimer FEU and DDU measurements were recorded at median 0 [IQR0-1] days since admission. Statistically and clinically meaningful differences were noted in d-dimer in the PE patients, compared to non-PE patients (Fig. 4). PE patients had higher d-dimer FEU (OR 1.13 [1.1–1.15]) and d-dimer DDU (OR 1.13 [1.06–1.21]) values. In ROC curve analyses, d-dimer FEU and d-dimer DDU were noted to have reasonable AUCs (0.70 for both) (Fig. 5). A clinical calculator indicating the sensitivity, specificity, positive predictive value, and accuracy of d-dimer in predicting acute PE in COVID-19 was devised and published online by our group [21]. For d-Dimer FEU, at 1.8 mcg/ml, the test reached the sensitivity 55%, specificity 74%, positive predictive value 34%, and accuracy 70%. At 0.8 mcg/ml cut off, d-dimer DDU reached the sensitivity 54%, specificity 75%, positive predictive value 37%, and accuracy 70%. At higher cut-offs, 7 mcg/ml d-Dimer FEU and 3.5mcg/ml d-dimer DDU, specificity and accuracy improved to 95% and 81%, respectively, at the expense of poor sensitivity (≤ 25%) (Additional file 1: Tables S2, S3).

Fig. 4
figure 4

Box plot comparison of d-dimer FEU, d-Dimer DDU, Heart Rate and Lymphocytes Count in COVID-19 with and without Acute Pulmonary Embolism. Difference in d-dimer (FEU and DDU) in patients with and without pulmonary embolism can be clearly appreciated. Very slight difference is noted in heart rate. No difference is noted in the lymphocyte count. Out of 31,500 COVID-19 hospitalized patients, 1117(3.5%) of these patients were diagnosed with acute PE in this study

Fig. 5
figure 5

ROC curve for d-Dimer FEU and D d-Dimer DDU predicting acute pulmonary embolism in COVID-19 patients. Area under curve for both d-Dimer FEU and d-dimer DDU was 0.70

Discussion

In the large population-based national cohort in the US, we report worse mortality and morbidity outcomes associated with the diagnosis of acute PE in hospitalized COVID-19 patients. D-dimer was noted to have utility in predicting pulmonary embolism starting at the cut-off levels 1.8 mcg/ml D-dimer FEU and 0.8 mcg/ml d-dimer DDU. We also present an admission d-dimer based calculator as a prediction tool for diagnosing acute PE in hospitalized COVID-19 patients.

The incidence of acute PE in hospitalized COVID-19 patients was noted at 3.5% in our study, which is lower than has been recorded in previous studies. In the French multicenter study of 1240 patients, the incidence of acute PE was noted at 8.7% in hospitalized COVID-19 patients; however, patients without CTPE (computerized axial tomography pulmonary embolism) imaging were excluded [6]. Similarly, the overall PE incidence of 9.5% in mechanically ventilated patients in our study is lower than has been previously reported for critically ill COVID-19 patients [25]. The lower incidence of PE in COVID-19 patients in our study represents the real-life clinical data, where—unlike protocol-driven studies—not all the patients are screened for acute PE. This points out the possibility of underdiagnosis of venous thromboembolism in this difficult-to-diagnose population. In our univariable analysis, COVID-19 patients with acute PE were more likely to have traditional risk factors (older age, obesity, history of malignancy) than patients without PE [26, 27]. Patients with PE were more likely to have cardiovascular risk factors; hypertension, diabetes, atrial fibrillation, and congestive heart failure [28, 29]. PE patients were sicker with a higher incidence of cerebrovascular diseases, acute encephalopathy, and sepsis. A high proportion of COVID-19 patients had acute kidney injury and chronic kidney injury (Table 1), more so in PE patients, which reflects upon the real-life dilemma that clinicians have to face when balancing the need for frequent image screening of PE against the high likelihood of contrast-induced kidney injury in this population.

Previous studies have noted increased mechanical ventilation rates consistent with the results of our study, but they have been underpowered to detect mortality differences [5, 6]. We report a clear association of mortality with PE diagnosis in COVID-19 patients. PE was mostly diagnosed earlier on during the admission and at the time of clinical deterioration during the peri-intubation period (Additional file 1: Fig. S2). Increased hospital length of stay was associated with the PE diagnosis. This underscores the importance of diagnosing acute PE in COVID-19 patients, which has an impact on important clinical outcomes.

Higher admission D-dimers of statistical and clinical significance were noted in PE patients than non-PE COVID-19 patients. Previous studies have reported an association of high D-dimer levels with increased mortality in COVID-19, which may be partly explained by the increased incidence of acute PE in patients with high d-dimer as noted in our study [30, 31].The close relationship of d-dimer levels with acute PE has been previously reported; however, it has not been established as a predictive risk tool for the evaluation of PE in COVID-19 [5, 32, 33]. Earlier smaller reports have reported variable d-dimer cut-off’s (ranging from 1 mcg/ml to 7.5 mcg/ml) to predict PE. However, the small size of the studies limits the ability to ascertain accurate prediction analysis of d-dimer [34,35,36,37,38]. In our study, at cut off values 1.8 mcg/ml d-dimer FEU and 0.8 mcg/ml d-dimer DDU, the test began to have clinical utility (accuracy 70%) for predicting acute PE in COVID-19 with specificity recorded at 74% for d-dimer FEU and 75% for d-dimer DDU. However, since the sensitivity was low (55% for d-dimer FEU and 54% for d-dimer DDU) at these cut-off values, the possibility of acute PE below the cut-off values cannot be ruled out. As the d-dimer values further increased, d-dimer had better ability to discriminate PE—specificity, positive predictive value, and accuracy improved—from the expected elevation in the setting of COVID-19 hospitalization; however, sensitivity continued to drop (Additional file 1: Tables S2, S3). We thus found d-dimer to be helpful in predicting pulmonary embolism in COVID-19 at the cut-off values 1.8 mcg/ml d-dimer FEU and 0.8 mcg/ml d-dimer DDU. We devised a clinical calculator (listing sensitivity, specificity, positive predictive values, and accuracy) using the threshold values of d-dimers to predict the probability of acute PE in hospitalized COVID-19 patients [21].

Patients with PE were more likely to have a history of PE or deep vein thrombosis with a high odds ratio of 5.11 [3.75–6.98], suggesting a lower threshold of screening for PE in this patient population. Chest pain was also more frequent in patients with PE. CRP and LDH were significantly increased in PE patients with COVID-19. This may still point towards the slightly higher inflammation noted in PE patients, supporting the hypothesis of ‘thromboinflammation’ for PE development [39, 40]. Slightly lower fibrinogen levels were noted in PE patients, as has been described before [41]. The heart rate differences pointed out between the two groups were not of clinical significance or usefulness (Fig. 4).

Our study has several limitations. As would be the case with any observational electronic health database study, the data retrieved via chart analysis may be subject to inaccuracies such as documentation errors. The number of patients that were screened with CT PE could not be retrieved due to the limitations of the database. The data may include patients who have not had the time to age into the outcome measurement yet, as the data is being regularly retrieved from the original sites. We adjusted for known confounders in our study, but the possibility of residual confounding exists, although sensitivity analysis indicated that addition of further covariates will not significantly affect the results (Additional file 1: Methods, Table S4).

There are several strengths of our study. We report by far the largest cohort analyzed for acute pulmonary embolism in COVID-19 patients from a geographically dispersed and demographically diverse multicenter population in the US. This allowed us sufficient power to delineate the effect of PE on COVID-19 mortality. The diagnosis of acute PE in COVID-19 is fraught with difficulty. For the first time, we provide a d-dimer based clinical calculator for predicting the probability of the diagnosis. The predictive risk tool underscores the utility and limitations of d-dimer at different threshold levels. We further defined optimal cut-offs where the utility of d-dimers is observed in predicting pulmonary embolism in hospitalized COVID-19 patients.

In conclusion, we report increased mortality and mechanical ventilation rates with acute PE in COVID-19 patients. d-dimer had potential utility in predicting pulmonary embolism in hospitalized COVID-19 patients. Further prospective studies should be performed to validate the pulmonary embolism risk model presented in our study.

Availability of data and materials

The data repository is located at the N3C Data Enclave. An institutional database agreement needs to be in place prior to joining the N3C for level II access which is required for data-sharing. Will individual participant data be available (including data dictionaries)?, All data has been de-identified in the N3C database. What data in particular will be shared?, Access to the de-identified database can be shared upon joining the N3C. What other documents will be available?, Not applicable. When will data be available (start and end dates)?, Immediately following the publication to one year post-publication (the project is hosted on the N3C enclave platform for one year). With whom?, Investigator who has institutional data use agreement with N3C. For what types of analyses?, Any purpose. By what mechanism will data be made available?, Data is accessible on the N3C enclave platform with the institutional database agreement.

Abbreviations

PE:

Pulmonary embolism

CTPE:

Computerized axial tomography pulmonary embolism

N3C:

National Collaborative COVID-19 cohort

LDH:

Lactate dehydrogenase

CRP:

C-reactive protein

ROC:

Receiver operative characteristic

AUC:

Area under curve

FEU:

Fibrinogen equivalent unit

DDU:

D-dimer units

aHR:

Adjusted hazard ratio

HR:

Hazard ratio

OR:

Odds ratio

IQR:

Interquartile range

References

  1. Johns Hopkins Corona Virus Resource Center Global Cases. Available from: https://coronavirus.jhu.edu/map.html. Accessed 24 Jan 2023.

  2. Morrissey SM, Geller AE, Hu X, et al. A specific low-density neutrophil population correlates with hypercoagulation and disease severity in hospitalized COVID-19 patients. JCI Insight. 2021;6(9).

  3. Ackermann M, Verleden SE, Kuehnel M, et al. Pulmonary vascular endothelialitis, thrombosis, and angiogenesis in Covid-19. N Engl J Med. 2020;383(2):120–8. https://doi.org/10.1056/NEJMoa2015432.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Suh YJ, Hong H, Ohana M, et al. Pulmonary embolism and deep vein thrombosis in COVID-19: a systematic review and meta-analysis. Radiology. 2021;298(2):E70–80. https://doi.org/10.1148/radiol.2020203557.

    Article  PubMed  Google Scholar 

  5. Ameri P, Inciardi RM, DiPasquale M, et al. Pulmonary embolism in patients with COVID-19: characteristics and outcomes in the Cardio-COVID Italy multicenter study. Clin Res Cardiol. 2020;41:1–9.

    Google Scholar 

  6. Fauvel C, Weizman O, Trimaille A, et al. Pulmonary embolism in COVID-19 patients: a French multicentre cohort study. Eur Heart J. 2020;41(32):3058–68.

    Article  CAS  PubMed  Google Scholar 

  7. Tritschler T, Kraaijpoel N, Le Gal G, et al. Venous thromboembolism: advances in diagnosis and treatment. JAMA. 2018;320(15):1583–94.

    Article  PubMed  Google Scholar 

  8. Kirsch B, Aziz M, Kumar S, et al. Wells score to predict pulmonary embolism in patients with coronavirus disease 2019. Am J Med. 2021;134(5):688–90. https://doi.org/10.1016/j.amjmed.2020.10.044.

    Article  CAS  PubMed  Google Scholar 

  9. Kwee RM, Adams HJA, Kwee TC. Pulmonary embolism in patients with COVID-19 and value of d-dimer assessment: a meta-analysis. Eur Radiol. 2021;31:9168. https://doi.org/10.1007/s00330-021-08003-8.

    Article  CAS  Google Scholar 

  10. Riyahi S, Dev H, Behzadi A, et al. Pulmonary embolism in hospitalized patients with COVID-19: a multicenter study. Radiology. 2021;301:426. https://doi.org/10.1148/radiol.2021210777.

    Article  Google Scholar 

  11. Weitz JI, Fredenburgh JC, Eikelboom JW. A test in context: d-Dimer. J Am Coll Cardiol. 2017;70(19):2411–20. https://doi.org/10.1016/j.jacc.2017.09.024.

    Article  CAS  PubMed  Google Scholar 

  12. NCATS. National Collaborative COVID-19 cohort (N3C) 2021. Available from: https://ncats.nih.gov/n3c. Accessed 17 Sept 2021.

  13. Haendel MA, Chute CG, Bennett TD, et al. The National COVID Cohort Collaborative (N3C): rationale, design, infrastructure, and deployment. J Am Med Inform Assoc. 2020;28(3):427–43. https://doi.org/10.1093/jamia/ocaa196.

    Article  PubMed Central  Google Scholar 

  14. Github. Phenotype Data Acquisition. Available from: https://github.com/National-COVID-Cohort-Collaborative/Phenotype_Data_Acquisition. Accessed 27 Jun 2021.

  15. Brazauskas R, Logan BR. Observational studies: matching or regression? Biol Blood Marrow Transplant. 2016;22(3):557–63. https://doi.org/10.1016/j.bbmt.2015.12.005.

    Article  PubMed  Google Scholar 

  16. Tuty Kuswardhani RA, Henrina J, Pranata R, et al. Charlson comorbidity index and a composite of poor outcomes in COVID-19 patients: a systematic review and meta-analysis. Diabetes Metab Syndr. 2020;14(6):2103–9. https://doi.org/10.1016/j.dsx.2020.10.022.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Baek S, Park SH, Won E, et al. Propensity score matching: a conceptual review for radiology researchers. Korean J Radiol. 2015;16(2):286–96. https://doi.org/10.3348/kjr.2015.16.2.286.

    Article  PubMed  PubMed Central  Google Scholar 

  18. McAfee AT, Ming EE, Seeger JD, et al. The comparative safety of rosuvastatin: a retrospective matched cohort study in over 48 000 initiators of statin therapy. Pharmacoepidemiol Drug Saf. 2006;15(7):444–53.

    Article  CAS  PubMed  Google Scholar 

  19. Ali MS, Prieto-Alhambra D, Lopes LC, et al. Propensity score methods in health technology assessment: principles, extended applications, and recent advances. Front Pharmacol. 2019;10:973. https://doi.org/10.3389/fphar.2019.00973.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Zöller D, Wockner LF, Binder H. Automatic variable selection for exposure-driven propensity score matching with unmeasured confounders. Biom J. 2020;62(3):868–84. https://doi.org/10.1002/bimj.201800190.

    Article  PubMed  Google Scholar 

  21. Gul M. Hamdan's PE Calculator 2021. Available from: https://www.covid19ddimer.com/. Accessed 18 Sept 2021.

  22. Korevaar DA, Josiean van ES. Pulmonary embolism in COVID-19: d-dimer threshold selection should not be based on maximising Youden's index. Eur Respir J. 2021;57(2). https://doi.org/10.1183/13993003.04279-2020

  23. MedCalc Software Ltd. Odds ratio calculator. Available from: https://www.medcalc.org/calc/odds_ratio.php (Version 20.013). Accessed 25 Sept 2021.

  24. Altman DG. Practical statistics for medical research. London, New York: Chapman and Hall; 1991.

    Google Scholar 

  25. Klok F, Kruip M, Van der Meer N, et al. Incidence of thrombotic complications in critically ill ICU patients with COVID-19. Thromb Res. 2020;191:145–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Eichinger S, Hron G, Bialonczyk C, et al. Overweight, obesity, and the risk of recurrent venous thromboembolism. Arch Intern Med. 2008;168(15):1678–83.

    Article  PubMed  Google Scholar 

  27. Stein PD, Huang HI, Afzal A, et al. Incidence of acute pulmonary embolism in a general hospital: relation to age, sex, and race. Chest. 1999;116(4):909–13. https://doi.org/10.1378/chest.116.4.909.

    Article  CAS  PubMed  Google Scholar 

  28. Rodríguez-Núñez N, Ruano-Raviña A, Lama A, et al. Impact of cardiovascular risk factors on the clinical presentation and survival of pulmonary embolism without identifiable risk factor. J Thorac Dis. 2020;12(10):5411–9. https://doi.org/10.21037/jtd-20-1634.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Sørensen HT, Horvath-Puho E, Lash TL, et al. Heart disease may be a risk factor for pulmonary embolism without peripheral deep venous thrombosis. Circulation. 2011;124(13):1435–41.

    Article  PubMed  Google Scholar 

  30. Soni M, Gopalakrishnan R, Vaishya R, et al. d-dimer level is a useful predictor for mortality in patients with COVID-19: analysis of 483 cases. Diabetes Metab Syndr. 2020;14(6):2245–9. https://doi.org/10.1016/j.dsx.2020.11.007[publishedOnlineFirst:2020/11/17].

    Article  PubMed  PubMed Central  Google Scholar 

  31. Yao Y, Cao J, Wang Q, et al. d-dimer as a biomarker for disease severity and mortality in COVID-19 patients: a case control study. J Intensive Care. 2020;8(1):49. https://doi.org/10.1186/s40560-020-00466-z.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Poyiadji N, Cormier P, Patel PY, et al. Acute pulmonary embolism and COVID-19. Radiology. 2020;297(3):E335–8. https://doi.org/10.1148/radiol.2020201955.

    Article  PubMed  Google Scholar 

  33. Valle C, Bonaffini P, Dal Corso M, et al. Association between pulmonary embolism and COVID-19 severe pneumonia: experience from two centers in the core of the infection Italian peak. Eur J Radiol. 2021;137: 109613.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Mestre-Gómez B, Lorente-Ramos RM, Rogado J, et al. Incidence of pulmonary embolism in non-critically ill COVID-19 patients. Predicting factors for a challenging diagnosis. J Thromb Thromboly. 2021;51(1):40–6. https://doi.org/10.1007/s11239-020-02190-9.

    Article  CAS  Google Scholar 

  35. Léonard-Lorant I, Delabranche X, Séverac F, et al. Acute pulmonary embolism in patients with COVID-19 at CT angiography and relationship to d-dimer levels. Radiology. 2020;296(3):E189–91.

    Article  PubMed  Google Scholar 

  36. Choi JJ, Wehmeyer GT, Li HA, et al. D-dimer cut-off points and risk of venous thromboembolism in adult hospitalized patients with COVID-19. Thromb Res. 2020;196:318–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Rodriguez-Sevilla JJ, Rodó-Pin A, Espallargas I, et al. Pulmonary embolism in patients with Covid-19 pneumonia: the utility of d-dimer. Arch Bronconeumol. 2020;56(11):758.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Mouhat B, Besutti M, Bouiller K, et al. Elevated d-dimers and lack of anticoagulation predict PE in severe COVID-19 patients. Eur Respir J. 2020;56(4):2001811. https://doi.org/10.1183/13993003.01811-2020.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Connors JM, Levy JH. Thromboinflammation and the hypercoagulability of COVID-19. J Thromb Haemost. 2020;18(7):1559–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Gu SX, Tyagi T, Jain K, et al. Thrombocytopathy and endotheliopathy: crucial contributors to COVID-19 thromboinflammation. Nat Rev Cardiol. 2021;18(3):194–209. https://doi.org/10.1038/s41569-020-00469-1.

    Article  CAS  PubMed  Google Scholar 

  41. Kucher N, Kohler H-P, Dornhöfer T, et al. Accuracy of d-dimer/fibrinogen ratio to predict pulmonary embolism: a prospective diagnostic study. J Thromb Haemost. 2003;1(4):708–13. https://doi.org/10.1046/j.1538-7836.2003.00145.x.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We acknowledge Dr. Philip Kern, James Aaron, and the CCTS University of Kentucky Department (funded by the NIH National Center for Advancing Translational Sciences through grant number UL1TR001998) for the help provided with the N3C database. We gratefully acknowledge the following core contributors to N3C: Adam B. Wilcox, Adam M. Lee, Alexis Graves, Alfred (Jerrod) Anzalone, Amin Manna, Amit Saha, Amy Olex, Andrea Zhou, Andrew E. Williams, Andrew Southerland, Andrew T. Girvin, Anita Walden, Anjali A. Sharathkumar, Benjamin Amor, Benjamin Bates, Brian Hendricks, Brijesh Patel, Caleb Alexander, Carolyn Bramante, Cavin Ward-Caviness, Charisse Madlock-Brown, Christine Suver, Christopher Chute, Christopher Dillon, Chunlei Wu, Clare Schmitt, Cliff Takemoto, Dan Housman, Davera Gabriel, David A. Eichmann, Diego Mazzotti, Don Brown, Eilis Boudreau, Elaine Hill, Elizabeth Zampino, Emily Carlson Marti, Emily R. Pfaff, Evan French, Farrukh M Koraishy, Federico Mariona, Fred Prior, George Sokos, Greg Martin, Harold Lehmann, Heidi Spratt, Hemalkumar Mehta, Hongfang Liu, Hythem Sidky, J.W. Awori Hayanga, Jami Pincavitch, Jaylyn Clark, Jeremy Richard Harper, Jessica Islam, Jin Ge, Joel Gagnier, Joel H. Saltz, Joel Saltz, Johanna Loomba, John Buse, Jomol Mathew, Joni L. Rutter, Julie A. McMurry, Justin Guinney, Justin Starren, Karen Crowley, Katie Rebecca Bradwell, Kellie M. Walters, Ken Wilkins, Kenneth R. Gersing, Kenrick Dwain Cato, Kimberly Murray, Kristin Kostka, Lavance Northington, Lee Allan Pyles, Leonie Misquitta, Lesley Cottrell, Lili Portilla, Mariam Deacy, Mark M. Bissell, Marshall Clark, Mary Emmett, Mary Morrison Saltz, Matvey B. Palchuk, Melissa A. Haendel, Meredith Adams, Meredith Temple-O'Connor, Michael G. Kurilla, Michele Morris, Nabeel Qureshi, Nasia Safdar, Nicole Garbarini, Noha Sharafeldin, Ofer Sadan, Patricia A. Francis, Penny Wung Burgoon, Peter Robinson, Philip R.O. Payne, Rafael Fuentes, Randeep Jawa, Rebecca Erwin-Cohen, Rena Patel, Richard A. Moffitt, Richard L. Zhu, Rishi Kamaleswaran, Robert Hurley, Robert T. Miller, Saiju Pyarajan, Sam G. Michael, Samuel Bozzette, Sandeep Mallipattu, Satyanarayana Vedula, Scott Chapman, Shawn T. O'Neil, Soko Setoguchi, Stephanie S. Hong, Steve Johnson, Tellen D. Bennett, Tiffany Callahan, Umit Topaloglu, Usman Sheikh, Valery Gordon, Vignesh Subbian, Warren A. Kibbe, Wenndy Hernandez, Will Beasley, Will Cooper, William Hillegass, Xiaohan Tanner Zhang. Details of contributions available at covid.cd2h.org/core-contributors. The analyses described in this publication were conducted with data or tools accessed through the NCATS N3C Data Enclave covid.cd2h.org/enclave and supported by NCATS U24 TR002306. The data-release partners are listed in Additional file 2: Appendix B.

Funding

The study did not receive any funding directly from a grant. We acknowledge Dr. Philip Kern, James Aaron, and the CCTS University of Kentucky Department (funded by the NIH National Center for Advancing Translational Sciences through grant number UL1TR001998) for the help provided with the N3C database. The N3C enclave is supported by supported by NCATS U24 TR002306. This specific study did not receive any funding.

Author information

Authors and Affiliations

Authors

Contributions

MHG, ZMH, JW, MS, MI conceived the idea. MHG analyzed the data and verified by SA and MV. All authors contributed towards interpreting the data, drafting, and critical revision of the work. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Muhammad H. Gul.

Ethics declarations

Ethics approval and consent to participate

The N3C Data Enclave is approved under the authority of the NIH Institutional Review Board for Protocol 000082 associated with NIH iRIS reference number: 546652 entitled: NCATS National COVID-19 Cohort Collaborative (N3C) Data Enclave Repository. An Institutional data-use agreement between the University of Kentucky and N3C data Enclave provided access to the de-identified database. This study did not need the consent to participate as the study was conducted on a de-identified database. The database was approved for use under the IRB aforementioned.

Consent for publication

Not applicable.

Competing interests

JW acknowledges the grant from National Heart Lung and Blood Institute (HL 129193) Pfizer Inc (Investigator-Initiated Grant) without any relationship to the current manuscript. All the other authors did not have any conflicts of interest to disclose.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

: Fig S1. Standardized Mean Differences of covariates before and after propensity score matching in the main analysis; Dot-plot Graph. The clinical conditions and the treatments are well-aligned. Fig S2. Temporal Relationship of Pulmonary Embolism diagnosis to Intubation. Fig S3. Box-plot comparison of CRP, Ferritin, Fibrinogen and LDH Count in COVID-19 with and without Acute Pulmonary Embolism. Table S1. Patient distribution in main and subgroup 1:4 propensity score-matched analysis. Table S2. Testing D-dimer FEU as a diagnostic test for Acute Pulmonary Embolism in COVID-19 for different cut off values. Table S3. Testing D dimer DDU as a diagnostic test for Acute Pulmonary Embolism in COVID-19 for different cut off value. Table S4: Sensitivity Analysis, outcomes of COVID-19 patients with and with-out acute pulmonary embolism after addition of three covariates; acute encephalopathy, stroke, and DVT. Table S5. Recurrent Mortality Analysis with Variable Censoring Times (Primary Analysis Data). Table S6. Recurrent Mortality Analysis with Variable Censoring Times (Sensitivity Analysis Data). Table S7. Recurrent Intubation analysis with Variable Censoring Times (Primary Analysis Data). Table S8. Recurrent Intubation Analysis with Variable Censoring Times (Sensitivity Analysis Data). Table S9. Diagnostic Codes for Clinical Conditions.

Additional file 2

: Appendix B.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gul, M.H., Htun, Z.M., de Jesus Perez, V. et al. Predictors and outcomes of acute pulmonary embolism in COVID-19; insights from US National COVID cohort collaborative. Respir Res 24, 59 (2023). https://doi.org/10.1186/s12931-023-02369-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12931-023-02369-7

Keywords