This retrospective study was approved by our institutional review board; patient informed consent was waived.
Patients
Between April 2015 and August 2019, 154 consecutive patients with suspected thoracic aortic aneurysms (TAA) underwent both contrast-enhanced CTA and TEVAR. Contrast-enhanced CTA combined with delayed-phase imaging performed one year later identified 14 patients with- (group 1) and 131 without endoleaks (group 2). Excluded from further study were 9 patients (renal dysfunction, n = 6; severe calcification, n = 3). As shown in Table 1, there was no significant inter-group difference in the age, sex, and body habitus.
Table 1
| Endoleak group | Non-endoleak group | p-value |
Age (year) | 79 (44–90) | 76 (43–93) | 0.12 |
Sex (M / F) | 11 / 3 | 104 / 27 | 0.27 |
Height (cm) | 162.8 (140.0–174.3) | 163.7 (130.0–185.0) | 0.39 |
Weight (kg) | 51.8 (35.5–73.2) | 58.5 (37.7–97.2) | 0.1 |
Thoracic CTA- and contrast medium injection protocol
All patients were scanned on a 64-detector-row CT scanner (Lightspeed VCT; GE Healthcare, Milwaukee, WI). Thoracic CTA was from the top of the lung apex to the lower end of the diaphragm in the craniocaudal direction. The scan parameters were 0.5-sec rotation scan, 1.25-mm detector row width, 0.984 helical pitch (beam pitch), 78.75-mm table movement, and 50-cm scan field of view. The tube voltage was 100 kVp; the tube current was changed from 200 to 770 mA to maintain the image quality (noise index = 9) using automatic tube current modulation. The scan length along the z axis ranged from 30 to 40 cm depending on the patient body size. The scanning time varied from 7.0 to 13.0 sec. The contrast material (Omnipaque-300; Daiichi-Sankyo, Tokyo, Japan) was injected through a 20-gauge catheter into the antecubital vein with a power injector (Dual Shot; Nemoto-Kyorindo, Tokyo, Japan). We used smart prep to obtain the scan delay for thoracic CTA. The injection volume was the total body weight in kg x 1.8 ml, administered during 30 sec.
TEVAR
All interventions were performed under aseptic conditions. The same angiography system (Alphenix Hybrid, Canon, Japan) was used by a multidisciplinary team that included interventional radiologists and vascular surgeons, each with at least 5 years of experience in the endovascular treatment of TAAs. Based on the vascular morphology, the stent-graft s were TAG (W. L. Gore & Associates, Flagstaff, AZ, USA), Zenith TX2 (William Cook Europe, ApS, Bjaeverskov, Denmark), RELAY PLUS (ABS Bolton Medical, Barcelona, Spain), VALIANT (Medtronic, Santa Rosa, CA, USA), and ZenithAlpha (Cook Medical, Bloomington, IN, USA).
The contrast material was Iopamidol 370 (BRACCO Imaging, Constance, Germany). The stent-graft was deployed in the descending aorta via the right femoral artery. The condition of the stent-graft was checked on angiograms; the entry tear exclusion was confirmed, as was the absence of endoleaks and of backflow from the distal re-entry site.
Machine-learning
One year post-TEVAR we compared the 145 patients with- (group 1, n = 14) and without endoleaks (group 2, n = 131) with respect to their age, sex, height, and weight, and the 22 vessel measurements shown in Table 2. We used data in their electronic medical records (Fujitsu, Japan) and in the Picture Archiving and Communications System (PACS). Vascular data were obtained on thoracic CTA images.
Table 2
Pre-procedure vessel measurements obtained in 14 patients with- and 131 patients without post-TEVAR endoleaks
| Endoleak group | Non-endoleak group | p-value |
Angle (º) | 19.0 (15.0–25.0) | 37.0 (7.0–47.0) | < 0.01 |
Distance between valsalva and brachiocephalic artery | 84.0 (68.0–113.0) | 84.0 (53.0–135.0) | 0.76 |
Short axis of brachiocephalic artery | 33.0 (23.0–41.0) | 32.0 (24.0–41.0) | 0.83 |
Short axis of vertebral artery | 31.5 (23.0–42.0) | 31.0 (24.0–46.0) | 0.76 |
Short axis of left subclavian artery | 31.5 (22.0–40.0) | 31.0 (24.0–45.0) | 0.58 |
Diameter of subclavian artery | 11.0 (9.0–14.0) | 11.0 (8.0–17.0) | 0.84 |
Proximal short axis of landing zone 1 | 3.0 ( 23.0–54.0) | 30.0 (18.0–42.0) | 0.22 |
Proximal short axis of landing zone 2 | 3.0 ( 23.0–54.0) | 30.0 (17.0–41.0) | 0.21 |
Proximal short axis of landing zone 3 | 3.0 ( 23.0–54.0) | 30.0 (18.0–42.0) | 0.23 |
Proximal long axis of aneurysm | 30.5 (24.5–55.0) | 31.0 (19.0–42.0) | 0.73 |
Proximal short axis of aneurysm | 30.5 (22.4–54.0) | 30.0 (18.0–41.0) | 0.44 |
Long axis of inner diameter of aneurysm | 47.0 (27.0–61.0) | 42.0 (16.0–75.0) | 0.65 |
Short axis of inner diameter of aneurysm | 39.5 (22.0–57.0) | 36.0 (11.0–65.0) | 0.12 |
Long axis of outer diameter of aneurysm | 49.0 ( 35.0–88.0) | 48.0 (31.0–95.0) | 0.44 |
Short axis of outer diameter of aneurysm | 42.0 (34.7–75.0) | 40.0 (29.0–94.0) | 0.24 |
Distal long axis of aneurysm | 30.0 ( 22.0–52.0) | 29.0 (21.0–44.0) | 0.69 |
Distal short axis of aneurysm | 29.5 ( 21.0–50.0) | 29.0 (21.0–43.0) | 0.87 |
Distal short axis of landing zone 1 | 28.5 (24.0–35.0) | 29.0 (22.0–37.0) | 0.6 |
Distal short axis of landing zone 2 | 29.0 (23.0–35.0) | 29.0 (22.0–37.0) | 0.87 |
Distal short axis of landing zone 3 | 28.5 (24.0–35.0) | 29.0 (22.0–37.0) | 0.71 |
Diameter between subclavian artery and aneurysm | 17.0 ( 10.0–54.0) | 32.0 (9.2–189.0) | < 0.01 |
Long axis of aneurysm | 115.0 (35.0–220.0) | 121.0 (35.0–241.0) | 0.65 |
TEVAR: Thoracic endovascular aortic repair |
Unless otherwise indicated, measurements are the mean and shown in mm |
We created machine-learning classifiers that applied extreme Gradient Boosting (XGBoost) (version 1.6.2; https://pypi.org/project/xgboost/) using Python (version 3.5; https://www.python.org/) and scikit-learn (version 0.18.1, http://scikit-learn.org/stable/) to discriminate between lumens with/without post-EVAR endoleaks. Scikit-learn is a python library and XGBoost is a library that can work with the same python. We divided the145 thoracic CTA patients into training data and test data in an 70/30 ratio so that the ratio of endoleak, without endoleak was as comparable as possible. As a result, we obtained 102 cases of training data (10 endoleak and 91 without endoleak) and 43 cases of test data (4 endoleak and 40 without endoleak). For parameter tuning we used the K-fold cross-validation function (K = 10), and applied oversampling for imbalanced data (synthetic minority over-sampling, SMOTE) [9], and Bayesian optimization (Hyperopt) [10]. The parameters for Bayesian optimization were: {'x_colsample_bytree': 0.8510707298078505, 'x_gamma': 0.45545493021965733, 'x_max_depth': 5.0, 'x_min_child': 2.0, 'x_reg_lambda': 0.6139857890064586, 'x_subsample': 0.9242267670478332}. Cross-validation was applied twice: the first time between training and validation data for tuning hyperparameters; the second time between training and test data after obtaining tuning data for model evaluation.
Machine-learning classifiers
For the post-TEVAR endoleak and non-endoleak classification we recorded the implanted stent-graft device (TAG, Zenith TX2, RELAY PLUS, VALIANT, ZenithAlpha), the patient age, sex, height, and weight (Table 1), and the 22 vessel standard set of measurements obtained on pre-operative CTA images (Fig. 1, Table 2). As the classification endoleak/non-endoleak with XGBoost using multiple items was unsatisfactory, we applied SMOTE for correcting imbalance data. Simple normalization was performed by first calculating the mean and standard deviation (SD) for each item. We then subtracted the mean from each item and divided the result by the SD. We trained the machine-learning classifiers on endoleaks and non-endoleaks by applying K-fold cross-validation to separate training sessions performed during testing of the discriminant using the classifiers (K = 10). The cross-validation function of scikit-learn automatically separated the population into K-fold cross-validation (K = 10). The parameters of XGBoost were selected by Bayesian optimization; to evaluate the importance of each item we applied the feature-importance algorithm.
Statistical analysis
Statistical analyses were performed with free “R” (R, version 3.2.2; The R Project for Statistical Computing; http://www.r-project.org/). Vessel measurements were performed between with the endoleak and without the endoleak group to identify the factors that affect for the endoleak. We compared the endoleak/non-endoleak groups with respect to 26 items using the Mann-Whitney U-test. The importance of individual features in the development of endoleaks was determined using XGBoost. We calculated their importance by applying XGBoost to machine learning and compared our findings between with those of conventional vessel measurement-based methods such as the 22 vascular features by using the Pearson correlation coefficients. To determine the strength of association, we took ‘r’ as the absolute value of the normal data distribution; r = 0–0.19 was recorded as very weak-, 0.2–0.39 as weak-, 0.40–0.59 as moderate, 0.6–0.79 as strong-, and 0.8–1.0 as very strong correlation. We compared the correlation coefficients using the method of Meng et al. [11].