Acute pulmonary embolism in hospitalized patients with SARS-CoV-2-related pneumonia: multicentric experience from Italian endemic area.

PURPOSE: To analyze pulmonary embolism (PE) on chest computed tomography pulmonary angiography (CTPA) in hospitalized patients affected by SARS-CoV-2, according to the severity of lung disease based both on temporal CT features changes and on CT-severity lung involvement (CT-severity score), along with the support of clinical and laboratory findings.
METHODS: We retrospectively enrolled a total of 170 patients with confirmed SARS-CoV-2 infection who underwent CTPA examination for PE suspicion. Pulmonary arteries diameters, right ventricle/left ventricle (RV/LV) ratio, presence, absence, and distribution of PE, pulmonary artery obstructive index (PAO index), and lobe involvement were recorded. All CT scans were reviewed to assess temporal CT changes and the COVID CT-severity score.
RESULTS: A total of 76 out of 170 patients (44.7%) developed PE without having any major risk factors for venous thromboembolism. The most severe pulmonary arteries involvement, expressed in terms of PAO Index, occurred in those patients with markedly elevated D-dimer and C-reactive protein (CRP) values and those patients with an advanced temporal stage of lung disease. The majority PE-positive patients were hospitalized in non-intensive wards. PE-positive patients showed a slightly higher hospitalization time in comparison with PE-negative ones. In the three months of study, overall 85.9% of patients were discharged while 14.1% died, of whom 13 PE-positive (54.2%).
CONCLUSIONS: Patients hospitalized for SARS-CoV-2 infection present a higher cumulative incidence of PE compared to the general population of hospitalized patients, regardless of the severity of lung inflammation or the temporal stage of the disease.

Introduction

Since December 2019, several cases of viral pneumonia caused by a novel coronavirus (COVID-19), identified as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), have been reported worldwide after the initial outbreak in Wuhan, China.

As of August 26, 2020, there were about 24 million confirmed cases of SARS-CoV-2 all over the world and 820,000 patients died with a mortality rate of 5.1%. Due to the high contagious rate, laboratory detection has become time-consuming and it may not be available for all people with suspected infection owing to the shortage of test kits for SARS-CoV-2. According to the current experience, pulmonary findings can be detected by imaging techniques in the initial phase of infection, earlier than clinical symptoms, and therefore imaging plays an important role in the preclinical screening, primary diagnosis, and management of SARS-CoV-2-related pneumonia [1-4].

According to Fleischner Society guidelines and European Society of Radiology advice paper, computed tomography (CT) can be considered a useful tool in patients with respiratory symptoms and not as a screening [5, 6].

CT features of SARS-CoV-2-related pneumonia have now been extensively described [7-10] and as reported in the recent literature, coagulation disorders are common and thromboembolic complications are increasingly being reported [11-14]. In this new SARS-CoV-2 infection, recent studies have shown a higher risk of developing thrombotic events [11-13], emphasizing the need to evaluate the D-Dimer value and other coagulation tests during the management, to regularly perform risk assessment. Abnormal coagulation results, especially markedly elevated D-dimer is associated with poor prognosis [15]. Infection, immobilization, respiratory failure, and mechanical ventilation are further considered potential thromboembolic-inducing factors in SARS-CoV-2 patients. Moreover, they might also be prone to develop venous and arterial thromboembolic disease due to excessive inflammation, hypoxia, and disseminated intravascular coagulation (DIC) [16-19].

On these bases, this study aims to analyze the pulmonary embolism detected on computed tomography pulmonary angiography (CTPA) in hospitalized patients affected by SARS-CoV-2, according to the disease severity based both on temporal CT features changes and on semiquantitative CT-severity lung involvement, along with the support of clinical and laboratory findings.

Methods

The study was performed under the Declaration of Helsinki; clinical and radiological data was anonymized during data collection. Given the retrospective nature of the study and in the presence of technical difficulties during the pandemic, informed consent will be waived, following article 89 of the GDPR EU Regulation 2016/679.

Patient population

From March 4th, 2020 until April 24th, 2020 we retrospectively enrolled all hospitalized patients with SARS-CoV-2 infection, confirmed through a pan-coronavirus conventional polymerase chain reaction (PCR) assay.

Data were extracted from the electronic medical records database of three different hospitals and tertiary care referral centers in Lombardy, a Northern Italian region.

Inclusion criteria were: (1) PCR-confirmed diagnosis of SARS-CoV-2 infection, (2) laboratory test, including D-dimer value, white blood cell (WBC) differential, C-reactive protein (CRP) values, (3) complete arterial blood gas (ABG) test. The exclusion criterion was the presence of severe breathing-induced artifacts on the CTPA scan.

All patients with clinical suspicion of PE including chest pain, worsening of respiratory symptoms, irregular or new-onset rapid heartbeat, worsening of fever, aggravation of ABG parameters, and a marked increase over time of D-dimer (normal range 262–62.000 ng/ml) and/or fibrinogen values (normal range 246–983 mg/dl) despite the anticoagulant therapy, were addressed to CTPA.

Patients characteristics

Patient history and major predisposing risk factors, including coagulation disorders, oral contraceptives, recent surgical procedures, pregnancy, cancer, fractures of the lower limb, and previous deep vein thrombosis, were recorded as well.

According to the guidelines implemented by the hospitals taking part in the study, patients with D-dimer score over 5000 µg/mL underwent a complete anticoagulation therapy (Enoxaparin sodium, 100 UI/kg twice a day), while patients with D-dimer score under 5000 µg/mL underwent prophylactic anticoagulation (Enoxaparin sodium, 100 UI/kg, once a day).

Laboratory findings, ventilation, and outcome

For each patient, we collected the following laboratory test results, performed on the same day of the CTPA: (1) D-dimer value, (2) white blood cell count (WBC), (3) neutrophils, (4) lymphocytes, (5) platelets (PLT), (6) CRP value, (7) fibrinogen, (8) antithrombin III (ATIII), and (9) international normalized ratio (INR). For each patient, ABG data, including pO2, pCO2, HCO3-, pH, SaO2, was recorded. P/F ratio, expressing the ratio between arterial pO2 and FiO2, was calculated and reported as well.

Information about ventilation, in particular noninvasive ventilation (NIV) and invasive ventilation by endotracheal tube (ET), were recorded.

We registered information about overall hospitalization time, in particular, dividing patients according to the addressing ward (non-intensive ward and intensive-care unit), time-to-discharge, and time-to-death in the period between March 4, and June 4, 2020.

CTPA technique

All CTPA examinations were performed using a 256-slice scanner (iCT Elite—Philips, Eindhoven), during contrast media injection (50 ml, Iobitridol 350—Xenetix 350 Guerbet, France; flow rate 3.5 ml/s) followed by saline flush (50 ml, flow rate 3.5 ml/s) by using bolus tracking technique, with a trigger level of 120 Hounsfield unit (HU) and a delay time of 8 s. The trigger area was manually placed in the main pulmonary artery. The parameters for CTPA scans were tube voltage 80 kV, automated tube current modulation, collimation 128 × 0.625, pitch 0.617, rotation time 0.27 s, field of view (FOV) 350, matrix 512 × 512, thickness 0.8 mm, and increment 0.4 mm. CT images were acquired during a single breath-hold. Multiplanar reconstruction on the sagittal and the coronal plane, and the maximum intensity projection (MIP) reconstruction images were obtained.

Image analysis

Image analysis was performed using the institutional digital database system PACS (Enterprise imaging, AGFA Healthcare, Belgium) by a radiologist with 15 years of experience in chest imaging and a resident radiologist with 4 years of experience, who independently reviewed the images. Final scores were determined by consensus. For each patient, the CTPA scan evaluation included: (1) main, right, and left pulmonary artery (MPA, RPA, LPA) diameters, (2) right ventricle/left ventricle (RV/LV) ratio, (3) presence of PE, (4) pulmonary artery obstructive index (PAO Index), (5) PE distribution, central o peripheral, (s6) the temporal stage of lung disease and (7) the extent of lung involvement.

The RV/LV ratio was calculated from the value of the minor axes of heart right and left ventricles, measured in the transverse plane at their widest points between the inner surface of the free wall and the surface of the interventricular septum. An RV/LV ratio lesser than 1.0 was considered normal.

The PAO index was calculated according to Qanadli et al. [20], as the product of N*D, where N is the value of the proximal clot site equal to the number of segmental branches arising distally, and D is the degree of obstruction, defined as 1 for non-occlusive vessels (D1) and 2 for totally occluded vessels (D2). Thus, an occluded embolus in a segmental artery was assigned a value of 2. A clot in a proximal artery that gives rise to segmental arteries was scored as though each segmental artery had an embolus with the same degree of obstruction as the proximal embolus. Subsegmental emboli were scored as a partial obstruction of a segmental artery. The maximum obstruction score value was 40 and then converted into a percentage.

The temporal CT changes which indicate the disease progression was evaluated using the temporal CT staging system according to recent literature findings reported by Pan et al. [21]: early stage (0–4 days) as normal CT or ground-glass opacities (GGO); progressive stage (5–8 days) as increased GGO and crazy-paving pattern; peak stage (9–13 days) as consolidation; absorption stage (> 14 days) as fibrous stripes.

Moreover, the lung CT involvement was assessed using a semiquantitative CT-severity score, previously reported by Liu et al. [22]: each of the five lung lobes was visually scored from 0 to 5 as follows: 0) no involvement, (1) less than 5% involvement, (2) 6–25% involvement, (3) 26–49% involvement, (4) 50–75% involvement, or (5) more than 75% involvement. The total CT score was the sum of the five individual lobar scores and ranged from 0 (no involvement) to 25 (maximum involvement).

Statistical analysis

Continuous variables were expressed as means and standard deviations and compared by using the Mann–Whitney test. Categorical variables were expressed as median value and interquartile range (IQR) and compared by using the Chi-square test or Fisher exact test, as appropriate. Relative risk (RR) and odds ratios (OR) were computed using crosstabs, with a 95% confidence interval (CI) as reference. A RR value of 1 implies that the presence of that variable poses no additional risk, RR greater than 1 or less than 1 implies additional or reduction in risk, respectively. OR greater than 1 represents an association between two variables. Correlations were computed with the Pearson or Spearman correlation coefficients, as appropriate. All tests were two-sided, and p < 0.05 was considered statistically significant. All statistical analyses were performed using the SPSS statistical package software (SPSS v26.0—Chicago).

Results

Demographic, clinical, and laboratory data of the entire cohort

By applying inclusion and exclusion criteria, a total of 170 patients were enrolled (n = 24 from Center 1, n = 30 from Center 2, and n = 116 from Center 3). The majority of patients were male, with a mean age of 63 years (± 12), with no major risk factors nor a history of PE (n = 152, 89.4%). All laboratory data are summarized in Table 1.

Table 1

CT and laboratory findings, medical and ventilation therapy in the cohort and between the two groups, according to the presence of PE. p values in bold denote a significant difference

N = 170	All	PE positive (n = 76)	PE negative (n = 94)	p value
Age (yo ± SD)	63 ± 12	62 ± 15	60 ± 14	0.519
Sex male (n; %)	116; 68.2	58; 76.31	58; 61.70	0.061
MPA (mm ± SD)	27.17 ± 3.42	27.74 ± 3.79	26.72 ± 3.06	0.238
RPA (mm ± SD)	21.87 ± 3.82	22.63 ± 3.54	21.26 ± 2.96	0.039
LPA (mm ± SD)	21.57 ± 3.47	22.68 ± 3.43	20.68 ± 3.29	0.006
RV/LV ratio (± SD)	0.99 ± 0.18	1.14 ± 0.17	0.87 ± 0.11	< 0.0001
D-Dimer (µg/ml ± SD)	6352 ± 12,362.41	32,231.26 ± 5223.97	10,051.80 ± 17,289.01	0.019
WBC (× 103/mm3 ± SD)	8.71 ± 4.06	9.75 ± 4.29	8.08 ± 3.85	0.114
Neutrophils (× 103/mm3 ± SD)	7.07 ± 3.65	7.81 ± 4.18	6.71 ± 3.39	0.492
Lymphocytes (× 103/mm3 ± SD)	1.21 ± 0.91	1.51 ± 1.37	1.07 ± 0.54	0.806
PLT (× 103/mm3 ± SD)	266.66 ± 128.82	270.42 ± 126.82	254.21 ± 132.29	0.621
CRP (mg/l ± SD)	10.38 ± 11.75	14.63 ± 14.74	14.68 ± 13.38	0.582
Fibrinogen (mg/dl ± SD)	586.06 ± 195.41	620.18 ± 195.13	558.36 ± 194.82	0.325
ATIII (% ± SD)	88.39 ± 19.39	85.17 ± 20.55	94.20 ± 16.53	0.089
INR (± SD)	1.24 ± 0.51	1.19 ± 0.28	1.29 ± 0.67	0.934
SaO2 (%)	92.96 ± 7.90	94 ± 8	92 ± 8	0.256
pO2 (mmHg ± SD)	95.15 ± 55.45	91.4 ± 49.5	97.4 ± 59.6	0.930
pCO2 (mmHg ± SD)	37.99 ± 14.01	40.0 ± 20.6	36.7 ± 8.0	0.807
HCO3− (mEq/l ± SD)	24.68 ± 5.71	24.9 ± 4.0	24.4 ± 6.6	0.902
pH (± SD)	7.0964 ± 1.5213	7.02 ± 1.72	7.15 ± 1.41	0.243
P/F ratio (± SD)	220.91 ± 140.56	217.52 ± 165.50	233.49 ± 122.50	0.774
CT-staging system (median; IQR)	12.5; 7–15.25	14; 9–17	11; 7–14	0.103
Medical therapy (n; %)	68; 40.0	32; 18.9	36; 21.2	0.518
NIV (n; %)	76; 44.7	34; 20.0	42; 24.7	0.520
ET (n; %)	26; 15.3	10; 5.9	16; 9.4	0.400

MPA, main pulmonary artery; RPA, right pulmonary artery, LPA, left pulmonary artery; RV/LV, right ventricle/left ventricle ratio; WBC, white blood cell; PLT, platelets; CRP, C-reactive protein; ATIII, antithrombin III; INR, international normalized ratio; NIV, noninvasive ventilation; ET, endotracheal tube

Imaging findings

The majority of patients were negative for PE at CT scan (n = 94, 55.3%). Patients with a CT scan positive for PE (44.7%) were more frequently male (M/F = 58/18, p = 0.061 ) with a mean age of 63 years (± 14). The majority of PE-positive patients showed a peripheral (n = 44, 57.9%) and bilateral (n = 42, 55.3%) distribution of the embolic event (Figs. 1 and 2), and a higher number of RV/LV enlargement (n = 68, 89.5%, p < 0.0001, OR = 71.40).

Fig. 1

CTPA images on the axial planes a, b and MPR c, d— Temporal CT-staging system Stage 1. a window width and level for the evaluation of mediastinum allow to correctly evaluate the presence of clot (arrows) inside the RPA and in the lobar branch for the left lower lobe. The clot appears as a hypoattenuating finding inside the vessel. b window width and level for the evaluation of lung parenchyma allow us to correctly evaluate the presence of multiple GGOs, distributed bilaterally, especially in the subpleural space. Moreover, reticulations are diffuse bilaterally. c MPR image with window width and level for the evaluation of mediastinum confirms the presence of the clot and its extent. The clot inside the RPA lumen is considered completely occlusive, while in the left side sub-occlusive. d MPR image with window width and level for the evaluation of lung parenchyma confirms the presence of GGOs, especially in both upper lobes in the subpleural spaces. No lung effusion was present

Fig. 2

CTPA images on the axial planes—CT-staging system Stage 3. a window width and level for the evaluation of mediastinum allow us to correctly evaluate the presence of clot (arrows) inside the lobar branch for the right upper lobe. b window width and level for the evaluation of lung parenchyma allow us to correctly evaluate the presence of diffuse and confluent GGOs with air bronchogram associated. c window width and level for the evaluation of mediastinum allow us to correctly evaluate the presence of small clot (arrow) inside the lobar branch of RPA. d window width and level for the evaluation of lung parenchyma allow us to correctly evaluate the presence of diffuse and confluent GGOs with air bronchogram associated. Moreover, it is possible to appreciate the presence of bilateral consolidations, with a subpleural distribution

Overall, the mean MPA diameter was 27.17 mm (± 3.42), RPA 21.87 mm (± 3.82), LPA 21.57 (± 3.47), and RV/LV ratio 0.99 (± 0.18). The mean RPA and LPA diameters and RV/LV ratio were significantly higher in the PE-positive group (p = 0.039, p = 0.006, and p < 0.0001, respectively). No significant difference was found in the MPA diameter between PE-positive and negative groups (p = 0.238). According to laboratory data, only D-Dimer showed a statistically significant difference between the two groups (p = 0.019), with a higher value in those patients positive for PE.

PAO index

The median value of the PAO index was 12.5 (2.5–50.0). Only six patients (7.9%) presented a PAO index > 40%. According to laboratory data, significant correlations with D-Dimer (r = 0.584, p = 0.002),

CRP (r = 0.384, p = 0.013), pH (r = 0.339, p = 0.029) and P/F ratio (r =  − 0.463, p = 0.031) values were found.

Temporal CT staging system

No significant difference according to the temporal CT-staging system, between the two groups, was found (p = 0.273).

CT scans, both negative and positive for PE, were more frequently characterized by a score of 3 (n = 54, 20.0% and n = 32, 18.8% respectively). The presence of PE had an OR = 2.26 to obtain a score ≥ 3 and an OR = 2.12 for a score equal to 4.

Positive correlations between temporal CT-staging system and MPA (r = 0.229, p = 0.018), RPA (r = 0.187, p = 0.047), LPA (r = 0.186, r = 0.045) diameters, RV/LV ratio (r = 0.415, p < 0.0001), were found.

Moreover, a positive correlation between the temporal CT-staging system and the PAO index was found (r = 0.291, p = 0.036) (Figs. 3 and 4).

Fig. 3

CTPA images on the axial planes. CT-staging system Stage 4. a window width and level for the evaluation of mediastinum allow us to correctly evaluate the presence of clot (arrows) inside the lobar branch for the left lower lobe. b window width and level for the evaluation of lung parenchyma allow us to correctly evaluate the presence of diffuse GGOs, in particular in the subpleural space of the left lung, associated with more evident reticulations. In the right lower lobe, it is possible to appreciate the presence of a consolidation

Fig. 4

CTPA images on the axial planes (a, b), and coronal plane (c, d) - CT-staging system Stage 2. (a, c) Window width and level for the evaluation of mediastinum allow us to correctly evaluate the presence of clot (arrows) inside the segmental and subsegmental branches of the RUL. (b, d) Window width and level for the evaluation of lung parenchyma allow us to correctly evaluate the presence of multiple and confluent GGOs in the RUL. (c)

CT-severity Score

The overall median value of the CT-severity score was 12.5 (IQR 7–15.25). No significant difference regarding CT-severity score values was obtained between the two groups [median value 11 (7–14) and 14 (9–17) respectively, p = 0.103].

Significant correlations between the CT-severity score and D-Dimer (r = 0.615, p=0.001), fibrinogen (r = 0.317, p = 0.015), ATIII (r =   − 0.572, p = 0.001), pCO2 (r = 0.199, p = 0.050), HCO3-(r = 0.220, p = 0.039), and FiO2 (r = 0.436, p = 0.001), were found.

Ventilation

Overall, 68 patients (40.0%) underwent medical therapy while 102 (60.0%) ventilation, of whom 76 (74.5%) NIV and 26 (25.5%) with ET. From a total of 76 PE-positive patients, 32 (42.1%) underwent medical therapy, while 44 (57.9%) ventilation [n = 34 (44.7%) NIV, and n = 10 (13.2%) invasive ventilation by ET]. The median PAO index value was significantly higher in patients who underwent ventilation, both NIV and ET (p = 0.010). According to laboratory data, CRP values and P/F ratio showed significant differences in patients who underwent ventilation (p = 0.002 and p = 0.013 respectively) (supplementary table).

The temporal CT-staging system between patients who underwent ventilation and those who did notshowed a statistically significant difference, both in terms of NIV and ET (p = 0.004), in particular for patients included in classes 3 and 4 (p < 0.0001). CT-severity score values were significantly higher in patients who underwent ventilation (p = 0.007), but without a significant difference in terms of hospitalization ward (p = 0.277).

Hospitalization and outcome

From a total of 76 PE-positive patients, 66 (87%) were hospitalized in non-intensive wards while only 10 (13%) were hospitalized in ICUs. PE-positive patients showed a slightly higher hospitalization time in comparison with PE-negative ones [mean 30 (± 18) vs 24 (± 19), median 26 (17–46) vs. 19 (7–40) days respectively], without statistically significant difference (p = 0.177). Moreover, hospitalization time did not show a statistically significant difference between ICU and non-intensive wards [mean 29 (± 22) vs. 26 (± 18), median 19 (15–36) vs. 22 (12–42) days respectively, p = 0.668] (supplementary table). In the 3-month period of study, 146 patients (85.9%) were discharged, while 24 patients (14.1%) died, of whom 13 PE-positive (54.2%) and 9 PE-negative (45.8%) (p = 0.365, OR = 1.54—95% CI 0.430–5.495). No statistically significant differences among all CT findings nor laboratory data and patient’s death were found (all p > 0.05).

Discussion

Along with other established risk factors, the presence of acute infections in both hospitalized and non-hospitalized patients is deeply associated with the increased risk of venous thromboembolism, which includes deep vein thrombosis and PE [21, 22].

The increased risk of PE after acute infection, especially the respiratory tract, might be caused by local inflammation leading to local activation of coagulation and local vasoconstriction induced by the hypoxic lung environment in the infected area [23].

Therefore, accurate knowledge of thrombotic complications, such as acute PE, in patients infected by SARS-CoV-2 is crucial for monitoring disease evolution and depicting worsening or other complications.

Among all patients enrolled in our study, 44.7% developed pulmonary embolism during hospitalization despite any major risk factors for venous thromboembolism and having undergone, according to their D-dimer score, either prophylactic or complete anticoagulation therapy before CTPA examination, according to the guidelines implemented by the hospitals taking part in the study. The causal relationship between PE and severe infection in our cohort is made clear by the absence of predisposing thrombo-embolic factors.

Some authors supported the existence of a significant correlation between SARS-CoV-2 and PE due to endothelial dysfunction and to a hypercoagulation state with the activation of platelets and leukocytes which result in thrombin generation and widespread fibrin deposition [17, 18]. This condition is further aggravated by hypoxia, which augments thrombosis by both increasing blood viscosity and hypoxia-inducible transcription factor-dependent signaling pathway [23-25].

More specifically, a recent study found a pro-coagulation pattern characterized by increased clot strength in SARS-CoV-2 patients, which significantly responded to thromboprophylaxis [25].

However, the reason for PE occurrence in SARS-CoV-2 patients is still unclear.

In our series, the most diffuse pulmonary arteries involvement occurred in those patients with markedly elevated D-dimer and CRP values, despite the high dose anticoagulation regimen, and in those patients with a low value of P/F. These results suggested that SARS-CoV-2 patients with higher levels of inflammation and higher values of D-dimer are more prone to develop pulmonary embolism underlining the importance of monitoring these biochemical values to promptly identify the patients at higher risk of PE, and quickly perform CTPA.

Although our cohort did not show a direct correlation between the development of an embolic event and the CT staging system, our results showed an increment of the occurrence of a severe thromboembolic event, according to PAO index values, in the advanced stages of lung involvement Interestingly, from pathophysiologic view, the positive correlations between temporal CT-staging system and MPA (r = 0.229, p = 0.018), RPA (r = 0.187, p = 0.047), LPA (r = 0.186, r = 0.045) diameters, RV/LV ratio (r = 0.415, p < 0.0001), indicate the increased resistance of pressure of the arterial pulmonary flow. These results demonstrate that SARS-CoV-2 infection, in combination with mechanical ventilation, can cause an increased RV afterload, with the consequent possibility to develop PE.

As for the CT-severity score, in line with the result of the recent paper of Gervaise et al. [26], our study did not show any significant difference concerning the extent of lung involvement in COVID-19 patients with or without PE. This supports the evidence that the increased risk of developing pulmonary embolism is not directly correlated to the severity of lung disease.

In our series, most of the registered thromboembolic events were bilateral and non-massive with a median PAO index value of 12.50. Only 6 out of 76 patients (7.8%) presented a PAO Index > 40%, indicating a massive embolic event. A slight majority of patients (57.9%) developed a peripheral segmental or subsegmental embolism.

Although we did not find a significant difference in terms of therapy, both medical and mechanical ventilation (noninvasive ventilation and endotracheal intubation), between patients who developed PE and those did not, our results indicates that patients with higher PAO Index values underwent more frequently ventilation, both NIV and ET. This may suggest that patients who developed a more severe embolic event are more prone to be hospitalized in ICU which underlines again the importance to evaluate the presence of pulmonary embolism and in particular the extent of the embolic event.

Furthermore, invasive and noninvasive ventilation was frequently found in association with the more advanced stages of parenchymal involvement.

Moreover, our results highlight that PE-positive patients presented the worst general health condition as they needed a slightly longer recovery than patients without PE.

In line with a recent paper [27], we did not find a significant difference in the ICU admission and the need for intubation between the two groups (PE positive and PE negative).

Indeed 87% of patients positive for pulmonary embolism were not hospitalized in ICU suggesting that every patient affected by COVID-19, even if they do not have severe health conditions to qualify for ICU care, can develop an acute pulmonary embolism.

A low mortality incidence of patients with PE and SARS-CoV-2-related pneumonia was registered in our series, maybe because of the prompt and rapid anticoagulant administration, and the strict observation regimen in those hospitalized patients.

In our series, we confirmed that PE should be considered a common complication in patients with SARS-CoV-2-related pneumonia, in line with recently published studies [13, 28, 29].

This study presents some limitations. Firstly, its retrospective nature. Secondly, the CTPAs were performed only in those patients with clinical and laboratory data suspicious for PE, and therefore the overall correct incidence and prevalence of PE in SARS-CoV-2 patients cannot be determined. Thirdly, CTPAs were usually performed when most of the hospitalized patients already presented advanced pneumonia, thus making an evaluation of the association between PE and CT-stage system difficult.

In conclusion, patients hospitalized for SARS-CoV-2 infection present a higher cumulative incidence of PE compared to the general population of hospitalized patients, regardless of the severity of lung inflammation or the temporal stage of the disease.

Therefore, it is required constant monitoring of D-dimer and CRP values along with the support of CTPA, for the evaluation of both parenchymal involvement and pulmonary artery obstruction, useful to address the correct clinical management of patients.

Supplementary file1 (DOCX 16 kb)