Incidence, risk factors, and thrombotic load of pulmonary embolism in patients hospitalized for COVID-19 infection.

OBJECTIVE: To determine the incidence, characteristics, and risk factors of pulmonary embolism (PE) among patients hospitalized for COVID-19. PATIENTS AND METHODS: We performed a prospective observational study of a randomly selected cohort of consecutive patients hospitalized for COVID-19 infection between March 8, 2020 through April 25, 2020. All eligible patients underwent a computed tomography pulmonary angiography independently of their PE clinical suspicion and were pre-screened for a baseline elevated D-dimer level.
RESULTS: 119 patients were randomly selected from the 372 admitted to one tertiary hospital in Valencia (Spain) for COVID-19 infection during the period of study. Seventy-three patients fulfilled both the inclusion criteria and none of the exclusion criteria and were finally included in the study. Despite a high level of pharmacological thromboprophylaxis (89%), the incidence of PE was 35.6% (95% confidence interval [CI], 29.6 to 41.6%), mostly with a peripheral location and low thrombotic load (Qanadli score 18.5%). Multivariate analysis showed that heart rate (Hazard Ratio [HR], 1.04), room-air oxygen saturation (spO2) (HR, 0.87), D-dimer (HR, 1.02), and C-reactive protein (CRP) levels (HR, 1.01) at the time of admission were independent predictors of incident PE during hospitalization. A risk score was constructed with these four variables showing a high predictive value of incident PE (AUC-ROC: 0.86; 95% CI: 0.80 to 0.93).
CONCLUSIONS: Our findings confirmed a high incidence of PE in hospitalized COVID-19 patients. Heart rate, spO2, D-dimer, and CRP levels at admission were associated with higher rates of PE during hospitalization.

Introduction

Since its first outbreak in Wuhan in late December 2019, the coronavirus infection known as COVID-19 has spread rapidly around the globe, with more than 55 million people infected and nearly 1350,000 deceased (Centers for Disease Control and Prevention [CDC], WHO; November 18, 2020). This virus, classified as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), seems to enter cells by endocytosis, using angiotensin converting enzyme-2 (ACE2) protein receptors found in the lungs, heart, kidneys, and gastrointestinal tract, as well as in the blood vessels. The abundant expression of ACE2 in type II alveolar cells may cause diffuse alveolar damage and hyperinflammation resulting in pneumonia and acute respiratory distress syndrome, as well as several extrapulmonary manifestations, including cardiovascular, hematologic, and thrombotic complications due to direct and indirect effects of the viral illness.3, 4, 5, 6

As regards the coagulopathy observed in association with SARS-CoV-2 infection, the mechanisms that activate coagulation have been hypothesized as being linked to immune responses, through the release of pro-inflammatory mediators (such as interleukin-1, interleukin-6 [IL6] and tumor necrosis factor-α) that interact with platelets, stimulate the expression of tissue factor (initiating extrinsic coagulation cascade), induce an upregulation of plasminogen activator inhibitor-1, suppress the fibrinolytic system, lead to endothelial dysfunction and complement pathway activation, triggering thrombogenesis.7, 8, 9, 10, 11

One of the most relevant and consistent features of this process is the increase in D-dimer levels, which has been associated in COVID-19 patients with a greater need for mechanical ventilation and admission to an intensive care unit (ICU), as well as death.12, 13, 14, 15, 16, 17, 18, 19

In the last few months, much scientific literature has been devoted to descriptions of different aspects of COVID-19 infection associated with thrombotic complications, particularly pulmonary embolism (PE). The results have been mixed, probably because most studies to date have been performed only on patients attended in an ICU and have not included any systematic or comprehensive investigation of prospective protocols.20, 21, 22, 23, 24, 25, 26, 27, 28 As a consequence, the incidence, characteristics, and risk factors of hemodynamically stable PE associated with COVID-19 infection still remain unknown.

Therefore, to address this critical gap in our knowledge, the main objective of the current study was to assess the incidence of PE in a representative series of patients randomly selected from all those hospitalized for COVID-19 infection. We also aimed to determine the risk factors and thrombotic load of PE in these patients, and to construct a predictive score for PE during hospitalization with easy-to-use variables obtained at the time of admission.

Methods

Patients and study design

We conducted a prospective observational study of consecutive hospitalized patients admitted to the respiratory and internal medicine departments with confirmed COVID-19 infection in one tertiary university hospital in Spain (Polytechnic and University La Fe Hospital. Valencia) between March 8, 2020 and April 25, 2020.

Each COVID-19 patient admitted into hospital was randomly assigned to one pulmonologist or one internal medicine specialist. Around 50% of the specialists involved (seven pulmonologists and three internal medicine specialists) were randomly selected to participate in the study (and all accepted to participate) and all their patients were initially included in the study, independently of their pre-test clinical probability of PE (n = 119). D-dimer levels were assessed before final inclusion on the basis that normal levels render acute PE unlikely and D-dimer is more frequently elevated in hospitalized patients, in infection or inflammatory disease, and in COVID-19 patients. ,

The inclusion criteria were a) patient age≥18 years and b) elevated D-dimer levels, defined as>500 ng/mL, except in patients older than 50 years, where an age-adjusted cut-off was used (age x 10 ng/mL). ,

The exclusion criteria were a) pregnancy; b) hemodynamic instability, defined as cardiac arrest or obstructive shock (systolic blood pressure < 90 mmHg or vasopressors required to achieve blood pressure ≥ 90 mmHg; c) contraindicated exposure to radiological contrast due to severe renal failure [creatinine clearance < 30 mL/min/1.73 m2]; d) known iodine allergy; e) supine intolerance, and/or f) any other indication of full-dose anticoagulation (i.e., atrial fibrillation, prosthetic heart valve).

Patients were treated at their attending physician's discretion, according to our internal protocol and their clinical judgement. All the eligible patients (n = 73) underwent a computed tomography pulmonary angiography (CTPA) to assess the presence of PE, in accordance with the research protocol.

The research protocol was approved by the ethics committee of the Hospital Universitari i Politècnic La Fe (Valencia, Spain), with the register number 2020–197–1. In accordance with this committee's requirements, all the patients provided informed consent to participate in the study orally, because of the biological risk involved in obtaining written consent. We followed the guidelines for observational studies of the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement.

Diagnosis of COVID-19 infection and pneumonia

Diagnosis of COVID-19 was confirmed either by a positive result in a real-time polymerase chain reaction (PCR) test of a nasopharyngeal or sputum sample or by a positive result from serological testing, COVID-19 suggestive opacities on chest X-ray and clinically compatible presentation (fever, cough, dyspnea, asthenia, myalgia, or diarrhea). , 34, 35, 36

The extent of the effects on the lung in a chest X-ray (CXR) was assessed by an expert thoracic radiologist, following our internal protocol. Each lung was divided on CXR in upper (suprahilar), medium (hilar) and inferior field (infrahilar) with a total of 6 lung fields. The affected lung extent was graded according to the fields involved: I: no lung opacities; II: mild extent (< 1 field); III: moderate extent (1–2 fields); IV: severe extent (3–4 fields); V: very severe extent (5–6 fields).

Clinical and analytical measurements

A complete set of clinical variables, including clinical history, clinical picture, onset of symptoms, previous treatments, comorbidities, physical examination, and daily evolution, was carefully recorded from all patients at admission and during hospitalization. Peripheral blood samples were collected at admission to analyze a set of acute phase reactants, including lactate dehydrogenase (LDH), total and differential leukocyte count, platelet count, D-dimer levels, IL6 levels, high sensitivity C-reactive protein (CRP), and albumin. Treatments used for the COVID-19 infection were also carefully recorded (Supplemental Table 1). The local D-dimer laboratory test was HemosILⓇ D-dimer HS500, a quantitative D-dimer assay with a threshold of 500 ng/mL. , The clinical risk of patients with confirmed PE was determined by the simplified Pulmonary Embolism Severity Index (sPESI). On this basis, the patients were divided into two groups: low-risk (0 points) and high-risk (1 point or more) (Supplemental Table 2).39, 40

Computed tomography pulmonary angiography

The images were acquired with a 256-slice multidetector CT (Philips Brilliance iCT) after intravenous injection of iodinated contrast agent (Iomeprol-400) at 4 ml/s, triggered on the main pulmonary artery.

Each CTPA was read by two independent experienced thoracic radiologists, blinded to patient status as well as to clinical and biological features. PE was defined by the presence of at least one intravascular contrast filling defect, and the closest pulmonary vessel affected was used to characterize the PE location as being in the main, lobar, segmental, or subsegmental pulmonary arteries. If any images proved difficult to interpret, a third independent radiologist examined them to provide a consensus.

The report also included the thrombotic load, measured via the Qanadli Score (Supplemental Figure 1), the pulmonary artery trunk diameter (PA) -at the level of the pulmonary artery bifurcation-, the ratio of the PA/Aortic diameter (Ao), the right ventricle (RV) diameter and the RV/left ventricle (LV) ratio. The extent of lung involvement was visually assessed as the percentage of lungs with COVID-19 suggestive ground-glass opacities or consolidations on CTPA. Five categories of involvement were described: I: no lung opacities; II: mild extent (≤ 25%); III: moderate extent (26–50%); IV: severe extent (51–75%); V: very severe extent (≥76%).

Statistical analysis

Continuous variables were reported as mean ± standard deviation (SD), and qualitative or dichotomous variables as percentages. Normal distribution of the variables was assessed by the Kolmogorov-Smirnov test. To analyze the characteristics of patients according to the presence or not of PE, a Student's t-test (or non-parametric test when appropriate) and Chi-Square test (or Fisher exact when appropriate) were used to compare the groups.

A sample size was calculated to analyze the estimated incidence of PE (main variable of the study) in the whole group of 372 patients hospitalized because of COVID-19 infection. According to the current literature, the expected incidence of PE in hemodynamically stable (mainly non-ICU) patients with elevated D-dimer and routine thrombosis prophylaxis (as similar as possible to our patients) is 10%. Assuming a calculation error of 6%, a confidence level of 95%, 73 patients with a valid reading of CTPA scan were needed. Agreement between radiologists on the diagnosis of PE was assessed via the Kappa index. The study was completed when the sample size of valid patients was reached.

Multivariable logistic regression was used to determine the optimal independent variables associated with the diagnosis of PE at admission. The variables included in the logistic model were those considered clinically relevant by the authors after an assessment of the current literature.

In order to construct a predictive score for PE in our patients, those quantitative variables related to the presence of PE in the logistic regression analysis were dichotomized. The cut-off point for each variable was selected using the best area under curve ROC (AUC-ROC) for the prognosis of PE via the Youden index. After dichotomization, these variables were again included in a logistic regression model. The beta coefficient value of each variable (rounded to the closest whole value) was selected to assess the relative weight of each variable in the model. The sum of the relative weight of all the variables included in the score corresponded to its total value. The score was divided into terciles and a post-test probability of suffering a PE during hospitalization was calculated for each tercile. Finally, a ROC-curve for the prognostic value of the score was constructed.

A Spearman test was used to assess the correlation between the Qanadli score and the analytical and CTPA measurements.

A two-sided p<0.05 was considered to be statistically significant and all the calculations were made with the corresponding 95% confidence interval (CI). Statistical analysis was carried out using SPSS software, version 24 (SPSS, Inc., Chicago, IL, USA).

Results

Characteristics of patients at baseline and treatment

In order to achieve 73 valid patients, we randomly screened a total cohort of 119 (32% patients from the total of 372 hospitalized from March 8, 2020 to April 25, 2020) with a confirmed diagnosis of COVID-19 infection. Although our initial protocol included the possibility of a serological test for inclusion in the study, in the end all the patients were included after a positive PCR test. Forty-six patients were excluded, mainly for a D-dimer value below the threshold level (25 patients) (Fig. 1 ).

Fig. 1

Flow chart of the study.

Abbreviations: CT, computed tomography; CTPA, computed tomography pulmonary angiography; DD, D-dimer; PE, pulmonary embolism.

Supplemental Table 3 shows that there were no differences between the randomly selected (n = 119) and non-selected (n = 253) groups of patients.

Table 1 shows the main characteristics of the 73 patients finally included in the study. Their mean age was 65.4 (SD 16) years (range: 28 to 92) and 52 (71%) were men. The mean body mass index (BMI) was 29.3 kg/m2 (SD 5.8). The most common comorbidities were arterial hypertension (50%), diabetes (18%), asthma (10%), and obstructive sleep apnea syndrome (10%). The most common symptoms at the time of admission were cough (76%) and dyspnea (45%). Most patients presented multi-lobe pneumonia (77%) and one third were moved to ICU (34%). During the hospitalization, patients received hydroxychloroquine or chloroquine (93%), azithromycin (92%), corticosteroids (46%), lopinavir/ritonavir (34%), tocilizumab (22%), and baricitinib (13%). Furthermore, 86.6% received hydroxychloroquine or chloroquine plus azithromycin, and 43% triple combination plus systemic corticosteroids. Significantly, most of the patients received pharmacological thromboprophylaxis (93%), mainly with enoxaparin 40 mg per day (70%) or 60 mg per day (19%). The mean hospitalization period was 27.3 days (SD 16.1).

Table 1

Baseline comparative characteristics of patients with and without pulmonary embolism. Data are expressed as mean (SD) or n (%).

Variable	All patients (n = 73)	With PE (n = 26)	Without PE (n = 47)	p
Age, yr	65.4 (16)	65.3 (14.5)	65.5 (16.9)	0.93
Gender, male/female%	71%/29%	69%/31%	72%/28%	0.78
BMI, Kg/m2	29.3 (5.8)	29.9 (6.2)	28.6 (5.5)	0.49
Symptoms,%-Fever (> 37 °C)-Dyspnea-Cough-Chest pain-Hemoptysis-Abdominal pain-Neurological disorder-Diarrhea-Limbs edema-DVT signs	77%45%76%17%4%6%1%30%3%3%	77%54%83%19%8%8%0%27%8%6%	78%40%72%16%2%13%2%31%0%0%	0.930.260.360.690.290.470.450.740.060.19
Treatment during hospitalization,%-Systemic steroids-Baricitinib-Lopinavir/ritonavir-Tocilizumab-Azithromycin-Chloroquine (or hydroxychloroquine)-Thromboprophylaxis	46%13%34%22%92%93%94%	48%14%42%39%88%93%88%	44%11%29%10%93%93%98%	0.870.370.260.0010.480.480.14
Cardiovascular disease,%-Arterial hypertension-Diabetes-Dyslipidemia-Stroke-Ischemic heart disease-Heart failure-Atrial fibrillation-Valvular disease	50%18%34%6%10%3%4%3%	50%8%35%4%8%4%4%0%	51%24%33%7%11%2%4%4%	0.920.080.910.620.640.690.900.29
Respiratory disease,%-Asthma-COPD-OSA-Smoke habit (current or former)	10%1%10%30%	15%0%12%38%	7%2%9%19%	0.240.450.720.13
Other risk factors,%-NSAID-Ongoing cancer-Thrombophilia-Previous VTE-Chronic venous disease-Previous surgery-Antiaggregant therapy	3%5%1%1%8%3%15%	3%8%0%0%15%4%16%	3%2%2%2%4%2%13%	0.890.390.450.450.110.690.71
At admission-Temperature, °C-Heart rate, bpm-Respiratory rate, rpm-Systolic BP, mmHg-Diastolic BP, mmHg-spO2/FiO2 ratio-room-air sp02,%	37.3 (0.9)94.7 (17)21.5 (5.5)133 (19.3)79.7 (9.8)429 (58)91.5 (8.6)	37.4 (0.8)102 (19)26 (15)133.6 (18.7)81 (10.5)411 (79)89.1 (9.8)	37.2 (1)88.4 (12.3)20 (5)132.3 (19.3)78.3 (9.6)440 (38)92.5 (7.9)	0.610.0010.010.820.400.0430.001
ICU,%	34%	42%	28%	0.25
Orotracheal intubation,%	27%	38%	19%	0.09
Onset of symptoms, days	7.3 (6.2)	6.3 (4.9)	7.8 (6.8)	0.31
Length of stay, days	27.6 (16.1)	27.8 (17.9)	27.4 (14.3)	0.92
CXR pneumonia score-Mild-Moderate-Severe-Very severe	23.3%27.4%20.528.7%	23.1%23.1%23.1%34.6%	23.4%29.8%19.1%27.7%	0.39
Analytical at admission-D-dimer, ng/mL-GPT, IU/mL-GOT, IU, mL-Ferritin, pg/mL-IL-6, pg/mL-LDH, IU/mL-CRP, mg/L-Leukocytes, cell/uL-Hemoglobin,%-Neutrophils, cell/uL-Lymphocytes, cell/uL-Platelets, (x103)-proBNP, mg/L-HS troponin T, mg/mL-INR-Quick Index,%-Fibrinogen, mg/dL-aPTT, secs	4327 (8321)40.9 (72)58.9 (77.6)704.5 (533)185 (379)328 (124)105 (107)8130 (3.3)14.7 (13.4)6245 (3.01)1289 (1.1)224 (87)672 (966)58.9 (142)1.1 (0.1)84.7 (14.1)651 (91)31.6 (15)	6270 (13,814)57 (110)79 (115)915 (618)226 (456)379 (147)154 (128)8991 (3.6)12.7 (1.8)6913 (3.2)1330 (1.6)236 (99)605 (1155)62 (23)1.1 (0.2)83.8 (13.7)690 (18)28.9 (3.7)	2384 (6134)31 (35)46 (44)572 (442)161 (326)297 (127)97 (96)7758 (3.1)15.8 (16.7)5923 (2.9)1120 (0.6)217 (80)686 (259)117 (53) 1.%2 (0.1) 85 (14.4) 630 (107)33.2 (18.3)	0.0060.160.270.110.040.040.0010.140.340.210.450.370.740.580.730.760.020.31

Abbreviations: aPTT, activated partial thromboplastin time; BMI, body mass index; BP, blood pressure; COPD, chronic obstructive pulmonary disease; CRP, c-reactive protein; CXR, chest X-ray; DVT, deep venous thrombosis; fiO2, inspiratory oxygen concentration; GOT, glutamate-oxaloacetate transaminase; GPT, glutamate pyruvate transaminase; HS, high-sensivity; ICU, intensive care unit; IL6, interleukin-6; INR, international normalized ratio; LDH, lactate dehydrogenase; NSAID, nonsteroidal anti-inflammatory drug; O2, oxygen; OSA, obstructive sleep apnoea; PE, pulmonary embolism; SpO2, oxygen saturation; VTE, venous thromboembolism.

Radiological findings and PE location

Twenty-six out of 73 patients included in the study had a confirmed diagnosis of PE with an incidence of 35.6% (95% CI, 29.6 to 41.6); 51.4% of PE had a sPESI=0 points and 48.6% had a sPESI>0 points. Supplemental Figure 2 shows the distribution of the observed PE. Most of the intraluminal filling defects were detected in peripheral (segmental and subsegmental) branches (45.1% of the right PE and 42.4% of the left PE of the total lumen contrast filling defects seen). No patient had a main artery trunk PE, 5 (19.2%) presented PE in the right or left artery trunk, and 18 (69.2%) presented at least one lobar PE. An example of patient with a lobar and subsegmental PE can be seen in Fig. 2 . The Kappa index between radiologists for the diagnosis of PE was 81% (100% for main, lobar and segmental pulmonary arteries and 69% for subsegmental pulmonary arteries). The mean time lapse between the onset of symptoms and the application of CTPA was 18.2 (SD 13.2) days.

Fig. 2

Computed tomography (CT) pulmonary angiography: A and B: a 66-year old female patient hospitalized for COVID-19 pneumonia; axial lung tissue setting (A) and soft tissue setting coronal (B) CT image reconstructions showing characteristic COVID-19 pulmonary lesions with ground-glass opacities and consolidations involving both the lung parenchyma in predominantly peripheral distribution (A) and a subsegmental contrast filling defect (white arrow, B). C, D and E: axial CT image reconstructions of a 38-year old male patient who had fever, cough, dyspnea, and left-side chest pain and was admitted with COVID-19 pneumonia; upper lung levels showing patchy ground-glass opacities in the left lung (C), and left lower lung with a triangular subpleural consolidation (D) corresponding to a pulmonary infarction (red arrow) related to the presence of pulmonary thrombus in the left lower lobe, visible in a soft tissue setting (white arrow, E).

Risk factors for PE

The comparative characteristics of patients with PE versus those without PE are shown in Table 1. At the time of admission, patients with PE presented lower room-air oxygen saturation (spO2), higher respiratory (RR) and heart rates, and higher values of systemic inflammation (D-dimer, LDH, CRP, and IL6 concentrations). We failed to find any differences in age, gender, body mass index, blood pressure, or extent of affected parenchyma lung involvement extent at admission. Furthermore, patients with PE were more frequently critically ill, transferred to ICU (42% vs 28%; p = 0.06), and in need of orotracheal intubation (38% vs 19%; p = 0.07).

Supplemental Table 4 shows the comparative CTPA characteristics of those patients with and without PE. There were no differences between groups as regards vascular and heart dimensions, or as regards the extent of pneumonic infiltration.

The following variables assessed at admission were included in a logistic regression model examining their potential clinical relevance: heart rate, RR, LDH, D-dimer and CRP concentrations, BMI, age, and room-air spO2. Of these, only room-air spO2, heart rate, and D-dimer and CRP concentrations at admission were independently associated with the occurrence of PE ( Table 2 ).

Table 2

Logistic regression. Baseline (hospital admission) quantitative variables independently associated with pulmonary embolism.

Variable	HR	95% CI	P value
Heart rate, bpm	1.04	1.01–1.09	0.036
Baseline room-air sp02,%	0.87	0.76–0.98	0.041
D-dimer, ng/mL	1.02	1.01–1.04	0.022
CRP, mg/L	1.01	1.01–1.05	0.037

Abbreviations: CI, confidence interval; CRP, C-reactive protein; HR, hazard ratio; spO2, oxygen saturation.

The chod score

The AUC-ROC method was used to determine the diagnostic value of each selected quantitative variable. The Youden index was used to dichotomize the selected variables and assess the cut-off point with the best prognostic value, and the corresponding AUC-ROC were calculated (Supplemental Table 5).

The dichotomized variables were included in another multivariable logistic regression in order to construct the score. The beta coefficients were rounded to establish the relative weight of each variable in the score (Supplemental Table 6).

Table 3 shows the final score that was constructed. It was named CHOD, an acronym of (CRP concentration + Heart rate + Oxygen saturation + D-Dimer levels).

Table 3

The CHOD score.

Variables (at hospital admission)	Values
C-Reactive protein<50 mg/L≥ 50 mg/L	01
Heart Rate< 90 bpm≥ 90 bpm	02
Oxygen Saturation (room air)> 92%≤92% or less	02
D-dimer< 956 ng/mL≥ 956 ng/mL	02
TOTAL	Range (0–7) points

The AUC-ROC of the CHOD score was 0.86 (95% CI: 0.80 to 0.93) ( Fig. 3 ). The score ranged from 0 to 7 points. The probability of incident PE during the hospitalization was low (4.5%) at 0–2 points, moderate (36.8%) at 3–5 points, and high (100%) at 6–7 points.

Fig. 3

Area under curve ROC of the CHOD score.

Factors associated with the thrombotic load

The mean value of the thrombotic load assessed with the Qanadli score was 18.5% (SD 15.9). Table 4 shows the significant correlations between the Qanadli index and various biomarkers and CTPA findings. A greater thrombotic load was seen in those patients with higher numbers of platelets and higher D-dimer values. The higher the thrombotic load, the larger the PA diameter, PA/Ao ratio, and RV/LV ratio.

Table 4

Correlations between thrombotic load and peripheral biomarkers and computed tomography pulmonary angiography (CTPA) findings.

Variables	Spearman correlation
At admission
Platelets	r = 0.40; p = 0.045
D-dimer	r = 0.47; p = 0.016
CT findings
PA/Ao ratio	r = 0.40; p = 0.046
PA diameter	r = 0.42, p = 0.034
RV/LV ratio	r = 0.35; p = 0.04

Abbreviations: Ao: Aortic diameter; LV; Left Ventricle; PA, Pulmonary artery; RV, Right ventricle.

Discussion

Our data, obtained from a randomly selected cohort of hospitalized patients with confirmed COVID-19 infection and D-dimer value above the threshold level, reveal that the cumulative incidence of PE was higher than expected (35.6%), thus suggesting that SARS-CoV-2 infection has a thrombogenic element. At admission, the risk factors for developing PE during hospitalization were a high heart rate, high D-dimer and CRP levels, and lower room-air oxygen saturation. These variables were used to construct an easy-to-use score that showed an excellent predictive value. The overall thrombotic load of PE was low, predominantly peripheral, and associated with various inflammatory biomarkers and CTPA findings. In our cohort, however, we failed to find any relationship between the extension of radiological pneumonia and the PE diagnosis. We suggest that physicians attending hospitalized patients with COVID-19 infection should be aware of their increased risk of PE, particularly in a severe form with systemic inflammation and hypoxemia.

The data currently available on thromboembolic risk in COVID-19 patients are still scarce and mainly refer to patients admitted to an ICU. . In ICU patients, a high cumulative incidence of 18–50% of venous thromboembolism (VTE) has been reported, depending on the research methodology used. , , 45, 46, 47, 48, 49, 50 As for non-ICU patients, Middeldorp et al. have observed that the cumulative incidence of VTE in those patients attended in an ICU (48%, 95% CI, 33 to 61%) was significantly greater than that seen in non-ICU patients (10%; 95% CI, 3 to 24%). Similarly, Bompard et al. retrospectively analyzed the CTPA performed on those COVID-19 patients with clinical suspicion of PE or elevated D-dimer levels. Again, PE incidence was greater in the ICU patients (50%, 95% CI, 30 to 70%), compared to non-ICU patients (18%, 95% CI, 12 to 27%). More recently, a retrospective single-center study found that the incidence of symptomatic PE in patients with COVID-19 was lower (5.4%), but CTPA was only performed on 14% (214/1477) of the subjects and the diagnostic yield of CTPA in these patients was similar to our findings (37%). However, these studies had several limitations, such as a lack of sample size calculation to estimate the incidence in the whole population of hospitalized patients, as well as a lack of analysis of the risk factors or thrombotic burden of PE in patients infected by COVID-19.

To the best of our knowledge, this is the first study to use a representative and random selection of patients to calculate the cumulative incidence of PE in the entire population of patients hospitalized with COVID-19 infection. The resulting incidence of 35.6% (95% CI, 29.6 to 41.6) is greater than any previously reported in non-ICU patients. Interestingly, this reported rate is similar to those previously reported in ICU patients with influenza A (H1N1) virus infection. Moreover, this is also the first study to analyze the main risk factors for suffering a PE during hospitalization, using clinical and analytical variables easily measured at the point of admission into hospital: increased heart rate, CRP, and D-dimer levels, and decreased room-air oxygen saturation. Significantly, all these variables could be altered by either PE of any origin or COVID-19 infection, but it seems that any alterations are more pronounced when both diseases appear simultaneously, or when COVID-19 infection produces a state of hypercoagulability that could lead to PE. In fact, current evidence suggests that D-dimer values are increased in COVID-19 patients, and this elevation correlates with the severity and prognosis of the disease, independently of the presence of PE.13, 14, 15 , , Furthermore, autopsy specimen studies have revealed that most lung samples present thrombotic features, especially in peripheral vessels.52, 53, 54, 55

Another novel finding of our study is the construction of a score to predict PE during hospitalization using four easy-to-measure variables (CRP, heart rate, oxygen saturation, and D-dimer) at admission. This CHOD score was constructed after a dichotomization of the variables, choosing the cut-off points with the best prognostic value (>50 UI/mL, >90 bpm, <92% and >956 mg/dl, respectively), and it has a range of values from 0 to 7 points. The prognostic value of the CHOD score was excellent (AUC-ROC: 0.86 [CI 95%, 0.80 to 0.93]). A CHOD score of 0 to 2 points at admission showed a low probability of PE during hospitalization (4.5%), compared with a very high probability (100%) in those patients with 6 or 7 points (and a moderate probability, 36.8%, in those with 3 to 5 points).

Thus, this score could be a useful tool for identifying those COVID-19 patients who, independently of the clinical suspicion of the attending physician, could benefit from a CTPA scan at admission because of a pro-thrombotic state caused by SARS-CoV-2. The predictive value of the CHOD score remained consistent in different series of patients with even very high cut-off points of D-Dimer levels (>2000 ng/ml) and different baseline characteristics.

It is worth noting that, despite the high incidence of PE, the vast majority of patients had no identifiable clinical symptoms suggestive of PE. This is probably because there was an overlap of symptoms of COVID-19 infection and PE, or because most of the PE observed were segmentary and subsegmentary asymptomatic thrombus. It is also striking that there is no association between the extension of radiological pneumonia (measured by both CXR and CTPA) and the confirmation of PE, even though a possible correlation could have been expected as both pneumonia extension and PE are related to a hyper-inflammatory state. These two situations could, however, correspond to different pathophysiological pathways: PE to a state of hypercoagulability with systemic inflammation and endothelial dysfunction, and the extension of pneumonia to a more local pulmonary inflammation. This possible divergence needs to be corroborated by further studies. Finally, it is also remarkable that PE had such a high incidence of PE in spite of the massive use of pharmacological thromboprophylaxis at correct doses (at least 40 mg of enoxaparin) in more than 80% of patients, reinforcing the theory of high pro-thrombotic risk from COVID-19 infection.

Finally, we observed that the thrombotic load, as assessed by the Qanadli index, was low, due to the high percentage of patients who presented effects on only segmentary and subsegmentary branches. Our results concur with those found by van Dam et al., who reported that PE patients with COVID-19 have a more peripheral and less extensive diagnosis than those with no COVID-19 infection, suggesting that PE in COVID-19 patients could be a combination of thromboembolic disease and in situ thrombosis. Not surprisingly, as in the case of PE patients with no COVID-19 infection, the thrombotic load of those with both diseases was associated with both more biomarkers of coagulation (number of platelets and D-dimer levels) and CTPA evidence of vascular repercussions, such as increased rates of PA/Ao and RV/LV.

The main strength of our study is that is the first to analyze the incidence and risk factors of PE in a consecutive hospitalized cohort of patients with COVID-19 infection, independently of the clinical probability of PE. However, some limitations also need to be considered: Our study is observational in nature, and therefore no causal association could be established. Our results are only applicable to patients with elevated D-dimer levels (although in the overall group of hospitalized patients in our study less than 20% presented D-dimer levels below the threshold).

Moreover, we did not screen for the presence of deep venous thrombosis, so the incidence of VTE could be even higher. On the other hand, by excluding patients with shock and cardiac arrest a number of PEs might have been lost and anticoagulation is not 100% protective against PE, so this could be a selection bias. Antifactor Xa levels were not assessed to monitor enoxaparin thromboprophylaxis dosing. Finally, a larger-scale external validation of our findings is still needed before our results can be applied to clinical practice.

In summary, our study confirmed the presence of a high cumulative incidence of PE in hemodynamically stable hospitalized patients with both a confirmed diagnosis of COVID-19 infection and increased D-dimer levels, despite appropriate thromboprophylaxis and independently of any clinical suspicion of PE. At admission to hospital, increased heart rate, CRP, and D-dimer levels and low oxygen saturation while breathing room air were the main risk factors for developing PE during hospitalization.

Declaration of Competing Interest

The authors declare no conflict of interest related to this study.