Introduction

A frequent and often lethal complication of a venous thromboembolism is a pulmonary embolism (PE). The in-cidence in the US is about 112 per 100,000, with rates of about 1/5th of this in Asian countries [1], and overall increasing rates in recent years [2,3]. A PE causes partial or complete blockage of arterial blood flow to the lungs [4]. Patients present with variable symptoms, and the specific diagnosis of PE is difficult to reach based on clinical information alone [4]. Conventional computed tomographic pulmonary angiography (CTPA) is a non-invasive and highly accurate method for early detection of PE, leading to early treatment [1,5]. A previous meta-analysis reported that the overall sensitivity and specificity of CTPA were 82% and 94.9%, respectively [6]. Dual-energy CTPA has recently been proposed as a technique to aid in the detection of PE [7-12]. The addition of lung perfusion information from a postprocessing iodine distribution map can enhance the reading results, without excess radiation or acquisition time [1,8-10,12-14].

Currently, functional imaging of the lung is a standard component of nuclear medicine using lung scintigraphy. Many studies have found that a dual-energy CTPA with an iodine distribution map increased the sensitivity of the PE detection from 67% to 89%, as well as providing more information than lung scintigraphy [1,12,15-17]. Dual-energy CTPA demonstrates anatomy, thrombus, other pulmonary diseases if present, and perfusion information, while lung scintigraphy provides perfusion information alone [18-20]. Moreover, dual-energy CT uses a single iodinated contrast injection. Thus, there is no need for subsequent lung scintigraphy with an additional radiopharmaceutical injection, which may lead to more allergic reactions to the radiotracer and injection site reaction than a single iodinated contrast injection [20-23].

Many studies found correlations between perfusion defects in PE patients and adverse clinical outcomes [18,24,25]. Perfusion defects revealed by lung scintigraphy were associated with dyspnoea, shorter distances covered in the 6-minute walk test, higher pulmonary vascular resistance, and higher mean pulmonary artery pressure [25,26]. The extension of perfusion changes can be used as a predictor of haemodynamic changes in acute PE [26]. Knowing a perfusion defect is present can help localize a PE clot and its severity, which is beneficial for treatment planning. Thus, the quantification of perfusion changes in PE patients could become an integral part of risk assessment [18].

Nowadays, conventional CTPA studies are increasingly used by clinicians, but they cannot provide perfusion information as can dual-energy CTPA. Although perfusion studies have many advantages, the required technologies are not available in most middle to lower income countries. Thus, it would be useful if there was a way to use information received from a conventional CT to predict the extent of perfusion defect.

To our knowledge, only a few studies have used their own CT scoring system on conventional CTPA to evaluate the degree of pulmonary vascular obstruction, assess the severity of PE, and predict perfusion defects on lung scintigraphy [15,25,27]. In the current study, we correlated the PE details from conventional CTPA against perfusion details from iodine distribution images. The objective was to predict the perfusion defects by clot size and pulmonary vascular stenotic degree on conventional CTPA images. This measurement method is simple and easy to use for every level of radiologist, whether with few or many years of experience. The data are highly quantitative, and thus consistent interpretations are easier to reach among radiologists of different levels of experience.

Material and methods

Study population

Our retrospective study was approved by the Human Research Ethics Committee (HREC). A total of 186 adults aged 15 years or above were enrolled in the study, all suspected to have PE because of major symptoms such as dyspnoea, hypoxia, chest pain, or leg swelling. All of them underwent a CTPA using a dual-energy CT machine between 15 November 2018 and 15 September 2020. Their baseline demographics were retrieved from the Hospital Information System (HIS). A picture archiving and communication system (PACS) inquiry was performed by the principal investigator (P.C.) to retrieve all dual-energy CTPA examinations. We used qualitative quality scores for qualitative control [14]. Fifteen CTPA studies were excluded due to poor imaging quality (2 studies), missing workstation data (8 studies), and the presence of artificial embolized mate-rials (5 studies), leaving 171 PE studies for analysis.

Imaging protocol, imaging reconstruction, and post-processing images

Conventional computed tomographic pulmonary angiography

All examinations were performed as a fundamental CTPA chest study protocol on a single-source CT scanner with a single detector layer (Revolution CT 1.0; GE Healthcare, Milwaukee, Wis). The dual-energy CT technique was applied only in arterial phase images, while the venous phase was reconstructed as a conventional image.

All patients received intravenous iopromide (Ultravist 370 mg) infused via peripheral venous catheter. A dual-bolus intravenous contrast technique was used to perform both the pulmonary arterial (PA) and venous phases, in line with the current literature [28-30].

The region of interest (ROI) was at the level of the main pulmonary artery with a trigger threshold of 100 Hounsfield units (HU). The scanning length covered from the aortic arch to the dome of the diaphragm.

Iodine distribution map

The iodine distribution maps were generated on a workstation computer with specific dual-energy post-processing software (GE Healthcare with AW Server version 3.2) [31-33]. For the iodine maps, a slice thickness of 2 mm with a 2-mm increment of the pulmonary arterial phase image at 70 KEV was generated.

Imaging interpretation

Conventional computed tomographic pulmonary angiography

All images were independently reviewed by 2 thoracic radiologists (13 years each of attending thoracic radiology experience; S.N. and N.K.), who were completely blinded. The CTPA images were interpreted on the pulmonary arterial phase with soft tissue reconstruction and mediastinal window (window width [WW] 350-500 HU and window-level [WL] 40-50 HU) [34,35].

The measurements for all affected luminal diameters were obtained by P.C. from the axial images and made from inner wall to inner wall and perpendicular to the wall at least twice. The largest diameter on the transverse axial image was used.

Each PE lesion was categorized by type, location, size, stenotic degree, and other pulmonary findings. Each clot size was classified either as large branch (main, lobar, or segmental PE) or small branch (subsegmental PE). The degree of stenosis was assessed on the cross-sectional view of the affected vessel. The expected percentage of obstruction was calculated from the value of the occluded area divided by the total luminal area using the area demarcated manually by the radiologist (Figure 1). The highest degree of stenosis of each size was used for analysis. Each radio-logist reviewed the images on 2 occasions separated by at least 2 days, then classified the degree of vascular stenosis into 5 categories by percentage reduction in vascular diameter as less than 25%, 25-49%, 50-74%, 75-99%, and 100% [27,36]. Any discrepancy between the 2 radiologists was reconciled by discussion until consensus was reached.

Figure 1

Measurement of stenotic degree at the affected pulmonary artery shows a reduction of 66%. A cross-sectional view of the highest degree of stenosis was reconstructed, then the area was demarcated manually by the radiologist

https://www.polradiol.com/f/fulltexts/146512/PJR-87-47891-g001_min.jpg

Iodine distribution map

Axial and coronal images were used for interpretation. A perfusion defect was defined as an area of low iodine distribution within the lung parenchyma demonstrating hypoattenuation on the iodine map. Patterns of perfusion defects were distinguished and classified into presence (triangular-shaped defect and peripheral location) or absence (normal lung perfusion, irregularly shaped defect, or band-like artifact) (Figure 2) [10]. The lowest density in each perfusion defect area by visualization was chosen for measurement of iodine concentration.

Figure 2

Pulmonary embolism on conventional computed tomographic pulmonary angiography (CTPA) and iodine distribution map. A CTPA with a mediastinal window shows a thrombus in the large branch of the right lower lobe (A, arrow), consistent with a triangular-shaped perfusion defect with peripheral location in the iodine distribution map (B and D, white frame). A band-like artifact with an irregular shape is also seen in the left upper lobe on the same image (B, dotted frame). The iodine concentrations using ROI analysis were measured (C). The value at the defect area (ROI 3) was 2.53 mg/cm3 and at the normal area (ROI 4) it was 26.77 mg/cm3. The iodine concentration ratio was 0.09

https://www.polradiol.com/f/fulltexts/146512/PJR-87-47891-g002_min.jpg

The iodine concentration was measured in milligrams per cubic centimetre (mg/cm3) by placing a 50 mm2 ROI at the same anatomical level (Figure 2), in order to assess the comparability between the scan with and without a perfusion defect. The iodine concentration ratio was calculated from the value at the defect area divided by the normal area.

Statistical analysis

Quantitative parameters are presented as mean ± standard deviation (SD) for continuous variables with normal distribution and median ± interquartile range (IQR) for those without normal distribution. The categorical variables were summarized using frequencies and percentages. The study analysis was computed by R program version 4.0.2 with RStudio server version 1.3.1073. Demographic and clinical characteristics were compared between the perfusion defect groups. The c2 test was used for categorical variables and t-test was used for continuous variables. Significant variables were noted as possible confounders. For further main hypothesis testing, lesions with and without perfusion defects were compared. Finally, GEE logistic regression was used to predict perfusion defects detected at each lesion, with the subject ID number being the repeating unit.

Results

From a total of 186 CTPA studies, only acute or chronic PE studies were enrolled, classified into 52 PE patients and 144 PE locations. PEs were detected in 16 patients and not detected in 36 patients. Of the PE locations, 45 lesions were rated as having a perfusion defect and 99 lesions were rated as no perfusion defect (Figure 3).

Figure 3

Study flow chart

https://www.polradiol.com/f/fulltexts/146512/PJR-87-47891-g003_min.jpg

Table 1 shows the patient characteristics compared between the group with and without perfusion defects using univariate analysis. There was no evidence of association between gender, age, major underlying disease, presenting symptoms, timing of clinical presentation, blood pressure, heart rate, respiratory rate, oxygen saturation, body weight, choice of treatment, and clinical outcome within 60 days of diagnosis.

Table 1

Demographic and clinical characteristics by perfusion defect (n = 52 subjects)

VariablePerfusion defect correlated with pulmonary embolismp-value
Present (n = 16)Absent (n = 36)
Gender
Male4 (28.6)10 (71.4)1.00
Female12 (31.6)26 (68.4)
Age at exam date (years)59.3 (13.4)60.4 (17.0)0.82
Major underlying disease
Cancer6 (23.1)20 (76.9)0.39
Hypercoagulative disease1 (25.0)3 (75.0)
Trauma2 (66.7)1 (33.3)
Other7 (36.8)12 (63.2)
Presenting symptoms
Hypoxia5 (38.5)8 (61.5)0.65
Dyspnoea4 (21.1)15 (78.9)
Chest pain1 (50.0)1 (50.0)
Leg swelling3 (42.9)4 (57.1)
Other3 (27.3)8 (72.7)
Timing of presentation (days)§1 (0, 1)1 (0.8, 5)0.14
Blood pressure (mmHg)
Systolic121.5 (25.2)126.1 (22.6)0.51
Diastolic73.1 (11.1)72.8 (12.2)0.93
Heart rate (beats/minute)94.3 (15.0)100.2 (20.6)0.31
Respiratory rate (breaths/minute)22.6 (2.7)24.5 (4.9)0.16
Oxygen saturation (%)§93.5 (90.8, 96.2)95 (90, 98)0.57
Body weight (kilograms)58.8 (12.0)55.4 (9.4)0.32
Choice of treatment
Medication14 (32.6)29 (67.4)1.00
Inferior vena cava filter1 (25.0)3 (75.0)
Thrombectomy0 (0)1 (100.0)
None1 (25.0)3 (75.0)
Clinical outcome within 60 days
Improved8 (25.0)24 (75.0)0.28
Not improved0 (0)2 (100)
Death8 (44.4)10 (55.6)

[i] Data are number with percentage in parentheses. Data are mean with SD in parentheses. §Data are median with interquartile range in parentheses.

Table 2 shows that there were no statistically significant differences in volume of contrast media, injection rate of contrast media, and quantitative quality score of the CT studies among the PE patients with and without a perfusion defect. Forty-five perfusion defect areas were measured for their iodine concentration. The mean and median concentrations were lower in the defect areas than in the normal areas. The iodine concentration ratios had a median of 0.11 and ranged from 0.01 to 0.61.

Table 2

Technical data and iodine concentration measurement

VariablePerfusion defect correlated with pulmonary embolismTotal (n = 52)p-value
Present (n = 16)Absent (n = 36)
Contrast media
Volume (ml)§40 (38.8, 40)40 (40, 50)40 (40, 50)0.21
Injection rate (ml/s)§4 (4, 4)4 (4, 4)4 (4, 4)0.61
Quantitative quality score
33 (18.8)4 (11.1)7 (13.5)0.77
48 (50)18 (50)26 (50)
55 (31.2)14 (38.9)19 (36.5)
Measured iodine concentrationMean (SD)MedianRange
Defect area (mg/cm3)2.48 (2.50)1.810.07-10.78
Normal area (mg/cm3)18.98 (8.28)17.215.71-41.73
Defect to normal area ratio0.14 (0.14)0.110.01-0.61

[i] Data are number with percentage in parentheses. §Data are median with interquartile range in parentheses.

The inter-observer agreement calculations are shown in Table 3. All assessments of PE including detection, type, size, location, stenotic degree, and detection of perfusion defect had excellent agreement.

Table 3

Interobserver variability

Imaging parameterObserved agreement (%)κ-valuep-value
Detection of pulmonary embolism100.001.00< 0.001
Type of pulmonary embolism93.060.85< 0.001
Size of pulmonary embolism82.640.72< 0.001
Location of pulmonary embolism100.001.00< 0.001
Stenotic degree81.250.74< 0.001
Detection of perfusion defect95.140.88< 0.001

A comparison of the CT findings between the PE patients with and without perfusion defect by PE location is presented in Table 4. The size of PE and the degree of vascular stenosis were significantly different between the 2 groups (p = 0.001 and p < 0.001, respectively). Lesions with higher degrees of vascular stenosis have higher percentages of perfusion defects. Other variables did not show statistically significant differences.

Table 4

Comparison of computed tomography (CT) findings by perfusion defect and pulmonary embolism (PE) location (n = 144 locations)

CT findingsPerfusion defect correlated with PEp-value
Present (n = 45)Absent (n = 99)
Type of pulmonary embolism
Acute32 (34)62 (66)0.422
Chronic13 (26)37 (74)
Location of pulmonary embolism
Right upper lobe12 (42.9)16 (57.1)0.523
Right middle lobe4 (20.0)16 (80.0)
Right lower lobe13 (32.5)27 (67.5)
Left upper lobe7 (26.9)19 (73.1)
Left lower lobe7 (30.0)23 (70.0)
Size of pulmonary embolism
Large and small branches27 (44.3)34 (55.7)0.001
Small branch18 (21.7)65 (78.3)
Degree of vascular stenosis
Less than 25%2 (10.5)17 (89.5)< 0.001
25-49%1 (7.1)13 (92.9)
50-74%6 (40.0)9 (60.0)
75-99%17 (26.6)47 (73.4)
100%19 (59.4)13 (40.6)
Other pulmonary findings
Consolidation5 (29.4)12 (70.6)0.573
Pleural effusion16 (29.6)38 (70.4)
Consolidation with pleural effusion6 (54.5)5 (45.5)
Others1 (20)4 (80)
Absence17 (29.8)40 (70.2)
Right heart strain
Presence13 (33.3)26 (66.7)0.899
Absence32 (30.5)73 (69.5)
Right/left ventricular ratio§1 (0.9, 1.6)0.9 (0.8, 1.2)0.192
Size of pulmonary trunk§2.9 (2.7, 3.1)2.9 (2.7, 3.2)0.978
Size of right pulmonary artery§2.1 (1.9, 2.2)2.1 (1.8, 2.4)0.479
Size of left pulmonary artery§2 (1.9, 2.2)2 (1.8, 2.2)0.977

[i] Data are number with percentage in parentheses. §Data are median with interquartile range in parentheses.

Table 5 summarizes the results of the generalized estimating equation (GEE) logistic regression model predicting the presence of perfusion defects. This multivariate analysis result confirms those from univariate analysis shown in Table 4.

Table 5

Generalized estimating equation (GEE) for perfusion defect

VariableOdds ratio95% CI
Type of pulmonary embolism
Acute (ref.)1.00
Chronic1.040.38, 2.88
Location of pulmonary embolism
Right upper lobe (ref.)1.00
Right middle lobe0.310.08, 1.20
Right lower lobe0.670.23, 1.96
Left upper lobe0.380.11, 1.22
Left lower lobe0.660.20, 2.20
Size of pulmonary embolism
Large and some small branches (ref.)1.00
Small branches0.700.49, 0.98
Degree of vascular stenosis
Less than 25% (ref.)1.00
25-49%0.530.04, 6.86
50-74%5.410.79, 37.24
75-99%2.890.56, 15.87
100%11.272.06, 61.81

[i] CI – confidence interval, ref. – reference

Discussion

Our study found no correlation between patient demographics, clinical characteristics, and anatomical data between the group of PE patients with and without perfusion defects. Univariate and multivariate analyses both found that size of PE and vascular stenotic degree could predict perfusion defect on conventional CTPA. The 2 significant predictors of perfusion defect were occluded vessels in both small and large branches together, or complete occlusion. Also, we found that lesions with higher degrees of vascular stenosis also have higher percentages of perfusion defect.

The clinical backgrounds and characteristics of our patients were similar to other studies, such as those of Sanchez et al. and Lee et al. [2,25]. The confirmed PE cases were predominantly found in women in the 6th and 7th decades of life. These elderly patients had higher prevalence of cardiac or pulmonary comorbidities due to a higher incidence of conditions associated with venous thromboembolic disease [37,38]. The women with PE were also older than the men, and they were more likely to have immobilization and congestive heart failure, which have been reported as significant risk factors for a thromboembolic event [39-42]. The common comorbidity in our study was malignancy, but there were no statistically significant differences in presenting symptoms, timing of clinical presentation, or choice of treatment. Even though our study conducted indirect assessments of the patients’ haemodynamic status from vital signs instead of clinical findings, as in a study by Sanchez et al., there were no significant differences in the vital signs between the groups with perfusion defects in the 2 studies [25]. This could be due to the populations of the studies having similar clinical backgrounds and characteristics.

The volume and rate of contrast medium injection were not associated with perfusion defects in our study. We had initially proposed that perfusion defects might differ by differing volumes and rates of contrast medium injection; however, the study findings falsified our hypothesis because these 2 variables were not significant predictors of perfusion defects. The main reason for this negative finding was probably due to the use of a standard protocol as the volume and rate of contrast injection were calculated based on each patient’s body weight [28-30]. To date there have been only a few studies reporting on quantitatively measured iodine concentrations [10,43]. The measurement was a feasible and reliable method to evaluate regional perfusion based on individual physiological variations. For these reasons, we created an iodine concentration ratio to reduce the effect of physiologic variation.

Interobserver agreement in our study was good, similarly to many prior studies, such as the one by Chae et al., which used subjective visualization for correlation of perfusion defects with the degree of pulmonary occlusion and location [24]. We used a simple technique for measurement and interpretation, including a vascular stenosis scoring system. With such high agreement, CTPA for diagnosis of PE should be reliable for general use by trained radiologists.

In this study, there was no significant association between perfusion defects and type of PE in the pulmonary arterial phase, a finding similar to that of a study by Hong et al. [44]. The Hong study also showed more enhancement of iodine attenuation at centrally located thrombi from the systemic collateral supply, which was found only in the delayed phase [44-46]. We did not perform delayed-phase dual-energy CT, which may have resulted in no association between perfusion defect and location of PE in our study. However, we chose a single pulmonary arterial phase in our routine practice to save time for the post-processing steps in our workflow. Unfortunately, we did not have a delayed-phase dual-energy CT during the study period, so we could not confirm whether it would have been useful.

Our study found that total occlusion of the pulmonary artery in conventional CTPA was associated with perfusion defects, a finding consistent with studies by Pontana et al. and Geyer et al. We found that pulmonary arterial occlusion involving both the large and small branches together was associated with perfusion defect rather than those small branches alone [8,10]. This result can be explained by noting that obstructed large vessels cause a large blood flow obstruction and decrease blood flow to the capillary bed, reflecting in increased perfusion defects in PE involving both large and small branches [19], while small vessel obstruction would promote systemic collateral supply in the normal lung via capillaries, resulting in fewer perfusion defects in PE involving only the small branches [47]. The evaluation of systemic collateral formation is important because it is a favourable prognostic factor for postsurgical outcome. Thus, the involvement of PE in only the small vessels may be an indicator of treatment success [44].

Lung scintigraphy has been used as a second-line diagnostic test because it has inferior spatial resolution to CTPA studies, inefficient for evaluation of anatomical data, and it takes a lot of time to perform [9]. Thus, an indirect method using conventional CT for assessment of pulmonary perfusion would be useful to provide additional information to improve diagnostic accuracy and predict clinical outcomes [18,24,25].

We found that perfusion defects could be predicted using conventional CTPA studies without any additional perfusion studies such as lung scintigraphy or iodine distribution map images. However, the imaging quality of the conventional CTPA study should be controlled, especially in small pulmonary arterial branches, including PE conspicuity, degree of pulmonary arterial enhancement, motion, and noise [14], because these common pitfalls in routine practice might cause limitations in interpreting the type and stenotic degree, leading to incorrect prediction of perfusion defects.

Our study had several limitations that should be noted. First, there were only data from a limited number of patients with pulmonary embolism, and no negative PE patients were enrolled from whom we could calculate the sensitivity and specificity. In addition, we did not perform delayed-phase dual-energy CT to evaluate lung perfusion. Further studies are recommended in both acute and chronic PE cases.

Conclusions

The measurement criteria on conventional CTPA images to predict perfusion defects were found with obstructions of both small and large vessels together or total occlusion of the pulmonary artery.