- 1Department of Vascular Surgery, The First Affiliated Hospital of Kunming Medical University, Kunming, China
- 2Department of Orthopedics, The First Affiliated Hospital of Kunming Medical University, Kunming, China
- 3Imaging Intervention Center, The First Affiliated Hospital of Kunming Medical University, Kunming, China
Objective: This study aims to establish a radiomics-based machine learning model that predicts the risk of transient ischemic attack in patients with mild carotid stenosis (30–50% North American Symptomatic Carotid Endarterectomy Trial) using extracted computed tomography radiomics features and clinical information.
Methods: A total of 179 patients underwent carotid computed tomography angiography (CTA), and 219 carotid arteries with a plaque at the carotid bifurcation or proximal to the internal carotid artery were selected. The patients were divided into two groups; patients with symptoms of transient ischemic attack after CTA and patients without symptoms of transient ischemic attack after CTA. Then we performed random sampling methods stratified by the predictive outcome to obtain the training set (N = 165) and testing set (N = 66). 3D Slicer was employed to select the site of plaque on the computed tomography image as the volume of interest. An open-source package PyRadiomics in Python was used to extract radiomics features from the volume of interests. The random forest and logistic regression models were used to screen feature variables, and five classification algorithms were used, including random forest, eXtreme Gradient Boosting, logistic regression, support vector machine, and k-nearest neighbors. Data on radiomic feature information, clinical information, and the combination of these pieces of information were used to generate the model that predicts the risk of transient ischemic attack in patients with mild carotid artery stenosis (30–50% North American Symptomatic Carotid Endarterectomy Trial).
Results: The random forest model that was built based on the radiomics and clinical feature information had the highest accuracy (area under curve = 0.879; 95% confidence interval, 0.787–0.979). The combined model outperformed the clinical model, whereas the combined model showed no significant difference from the radiomics model.
Conclusion: The random forest model constructed with both radiomics and clinical information can accurately predict and improve discriminative power of computed tomography angiography in identifying ischemic symptoms in patients with carotid atherosclerosis. This model can aid in guiding the follow-up treatment of patients at high risk.
Introduction
In the United States, 690,000 patients experience ischemic stroke yearly. Ischemic brain injury is caused by large-arterial atherosclerosis, cardioembolism, small-vessel disease, or cryptogenic cirrhosis. Carotid bifurcation plaques in atherosclerosis are prone to transient ischemic attack (TIA) and stroke. Among all patients with ischemic stroke, 173,000 (approximately 25%) suffer from atherosclerotic carotid artery disease (1). Endarterectomy and stenting are used to treat patients with 70–99% carotid stenosis (2–4), but moderate stenosis is more common and the main culprit of the plaques (5).
In addition, <50% of patients with stenosis are more difficult to treat and suffer from recurrent neurological disorders such as stroke and TIA (3). According to the 2017 European Society for Vascular Surgery (ESVS) clinical practice guidelines, patients with carotid stenosis who had cerebral infarction or TIA in the past 6 months were divided into symptomatic carotid stenosis and asymptomatic carotid stenosis groups. Best medical therapy (BMT) is recommended for patients with asymptomatic carotid stenosis with <60% stenosis and those with symptomatic carotid stenosis with <50% stenosis. In the past 5 years, the risk of stroke in 1,429 patients with 30–49% stenoses randomized to carotid endarterectomy (CEA) was 22.8%, compared with 25.5% on BMT. If symptoms persist, despite BMT, it is recommended to undergo carotid endarterectomy or carotid artery stenting (CEA/CAS) treatment (3).
Carotid computed tomography angiography (CTA) can measure carotid lumen stenosis and provide information about arterial wall calcification. While the stenosis degree provides crucial information about the disease, it does not determine the underlying plaque stability or inflammation degree and the probability of the patient experiencing a second event (6). Radiomics is a multi-step process that transforms medical images into high-dimensional structures to comprehensively analyze the regions of interest and correlate them with clinical, diagnostic, and prognostic information (7, 8). Standardized radiomics analysis comprises image acquisition, reconstruction, image preprocessing and processing, feature extraction, selection, and classification/regression modeling (6, 9, 10). Radiomics and machine learning are used to diagnose disease diagnosis and prognosis, limited to dermatology (11), oncology (12–14), and cardiac diseases (15). Radiological studies that evaluate carotid disease, especially in mildly stenotic carotid CTA imaging, are limited. Therefore, this study aims to identify patients' CTA region of interest and clinical information (6). Machine learning methods were used to establish a model that predicts the risk of TIA events with carotid artery stenosis by 30–50% and guides patients' follow-up treatment plans.
Materials and methods
Patients
This study was approved by the ethics committee of The First Affiliated Hospital of Kunming Medical University (No. 2022-274). Written informed consent from all subjects (patients) was waived by the Ethics Committee of The First Affiliated Hospital of Kunming Medical University, because of the retrospective nature of the study. The flowchart of the current research protocol was presented in Figure 1. Nearly 615 patients with carotid artery stenosis were selected from October 2016 to March 2022. The inclusion criteria comprised patients undergoing carotid CTA and carotid artery stenosis of 30–50% with clinical laboratory tests (uric acid, triglycerides, low-density lipoprotein, homocysteine, and fibrinogen). The patients were excluded if they had (1) carotid artery dissection and aneurysm; (2) intracranial vascular disease (e.g., intracranial atherosclerosis with stenosis <50%, vasculitis, aneurysm); (3) posterior circulation stroke; (4) intracerebral hemorrhage; (5) other causes of hemorrhagic stroke (e.g., cardioembolic source and chest embolism); or (6) patients with other neurological diseases such as brain tumors and demyelinating diseases.
Figure 1. Flowchart for building a predictive model. (A) The carotid CTA image of a patient. (B) The manual VOI segmentation was performed on an image. (C) The radiomic features extracted from each VOI. (D) The feature data were selected. (E) The establishment of prediction models. CTA, computed tomography angiography; VOI, volume of interest; ROC, receiver operating characteristic; AUC, area under curve.
After the screening process, 179 patients were identified and 219 carotid arteries were included in this study. Follow-up continued until the outcome occurred or November 2022, whichever came first. Through follow-up, 34 patients (including 39 carotid arteries) developed TIA. There were no recurrences of stroke. The patients were divided into two groups: the symptomatic group: patients with symptoms of transient ischemic attack after CTA and the asymptomatic group: patients without symptoms of transient ischemic attack after CTA (Table 1).
Stenosis was measured according to the North American Symptomatic Carotid Endarterectomy Trial (NASCET) criteria (3). The symptoms were defined as patients who were considered symptomatic after a TIA or ischemic stroke in the ipsilateral cerebral hemisphere. TIA is defined as brief (24-h) episodes of neurological dysfunction, such as hemiplegia, sensory disturbances, dysarthria, speech disturbances, or monocular blindness.
Image acquisition
Carotid computed tomography angiography imaging (Siemens SOMATOM Definition Flash DSCT) was used to scan the area from the aorta to the atlas. About 60–80 mL of non-ionic contrast agent Iopromide injection (360 mg/mL) and 50 mL of normal saline were injected from the cubital vein at a flow rate of 5 mL/s. Smart prep technology was used to observe the density of the targeted region in the aortic arch, and the scan was administered when the density reached 150 HU. The scanning parameters were tube voltage (120 kV), automatic tube current modulation technology, layer thickness (0.625 mm), layer spacing (0.625 mm), pitch (0.938:1), and rotation speed (0.5 r/s).
Region of interest segmentation
Soft plaques have attenuation values of <50 HU. A mixed plaque is one that has attenuation values between 50 and 119 HU. Calcified plaques have attenuation values of >120 HU (16). During the image examination, window/level (window 900–1,200 HU, level 300–500 HU) to view calcified and non-calcified plaque. The calcified plaque was easily distinguished from non-calcified plaque by manually and visually manipulating the window/level. Vascular luminal contrast enhancement was easily separated visually from calcified plaque by the manual manipulation of window/level settings which allowed for optimal plaque component distinction (17). The manual volume of interest (VOI) segmentation was performed on all images using the 3D Slicer software package (version 4.11) (18). The carotid artery with arterial plaque was manually segmented into the targeted region. The proximal and distal carotid arteries were divided into normal vessels at least 2 cm away from the plaque. Figure 2 shows the sample patient images. The VOI of optimal size was accurately segmented to establish a precise relationship between the radiomics of carotid plaque and the transient ischemic attack incident.
Figure 2. The procedure of radiomics segmentation. The VOI in each image was segmented with the green arterial lumen and yellow plaques, and a 3D model was synthesized. (A) The edges of atherosclerotic plaques were delineated. (B) Atherosclerotic plaques were segmented. (C) The edges of arteries were delineated. (D) Arteries were segmented. (E) A 3D model was built based on the segmentation of each layer image. (F) Carotid artery CTA corresponding image of the patient. VOI, volume of interest; CTA, computed tomography angiography.
Radiomics feature extraction
An open-source package PyRadiomics (version 2.4, https://pyradiomics.readthedocs.io/) was used to extract the radiomic features. The radiomic features extracted from each volume were automatically calculated using Pyradiomics in Python. They were divided into eight groups: (a) first-order features, (b) shape features (3D), (c) shape features (2D), (d) grayscale co-occurrence matrix features, (e) grayscale size regions matrix function, (f) grayscale run-length matrix feature, (g) adjacent grayscale difference matrix feature, and (h) grayscale dependency matrix feature. A total of 129 radiomics features were extracted from each patient's VOI (19).
Radiomics feature selection
The primary features that significantly impact the occurrence of TIA were identified by reducing the number of extracted features. About 129 radiomic features were converted into digital form through feature data. Seventeen features that could not be converted into digital form were excluded, and the remaining 112 features were included in the feature set. A feature selection algorithm based on the RF (20) method was used as the basic tool, incorporating the features into the importance order from high to low. Each time a feature was included, the optimal number of features was determined using the logistic regression model, and the classification area under curve (AUC) value was calculated. A total of 112 AUC values were calculated, of which the numbers of features corresponding to the top three AUC values were 73, 3, and 9, respectively. The prediction models were constructed with the features, respectively. The prediction performance of each model was optimal when 73 features were included, but it was not parsimonious. When only three features were included, the prediction performance of the subsequent prediction model was poor, thus, we chose to include the first nine features.
Clinical features and combined features selection
The clinical features included four components, such as routine blood test data (triglyceride, low-density lipoprotein, homocysteine, uric acid, and fibrinogen (21–25), which were demonstrated to be closely related to carotid plaque formation according to previous studies), demographic data [age, gender and body mass index (BMI), high blood pressure, diabetes, and smoking (26–31)], stenosis of the carotid artery, and medication data [antiplatelet and lipid-lowering drugs (32)]. We performed the same process as radiomics feature selection. The final result was a model with better prediction performance when the first two features were included.
Similarly, as earlier, the radiomics and clinical feature data were merged into the same dataset, and feature screening was performed. As a result, the first three important features were selected to build the model (Figure 3).
Figure 3. Feature importance in the random forest model. Radiomics features and combined features show only the top 20 features. Clinical features show all features. (A) The random forest method was used to rank clinical features from high to low importance. (B) The random forest method was used to rank radiomics features from high to low importance. (C) The random forest method was used to rank combined features from high to low importance. BMI, body mass index; UA, uric acid; TG, triglycerides; LDL, low-density lipoprotein; Fib, fibrinogen; HCY, homocysteine.
Predictive model building and model validation
Five classification algorithms were used, namely, RF, eXtreme Gradient Boosting (XGBoost), logistic regression (LR), support vector machines (SVMs), and k-nearest neighbors (KNNs). The hyperparameters for each classifier were tuned via a grid search process to maximize the model's performance. The two groups of patients were randomly assigned; 70% to the training set (N = 165) and 30% to the testing set (N = 66). Due to the limited sample size of the dataset, the 5-fold cross-validation method was used to train the training set and the best threshold corresponding to the maximum Youden's index was selected in the training set and applied in the testing set (33, 34). The models were validated in the testing set (33, 34). The 95% confidence intervals of the training and testing sets were obtained using bootstrap methods (repeated 1,000 times) (35). For each of these five different algorithms, the screened radiomics, clinical, and combined features information mentioned above were used to construct models, respectively.
Statistical analysis
Continuous variables were recorded as mean ± standard deviation, while categorical variables were recorded as percentages. The results were compared using a t-test for continuous variables and a χ2 test for categorical variables. The percentages of missing values were as followed: fibrinogen (18.7%) and homocysteine (21.9 %). Missing values were filled with the mean values.
The model prediction performance was determined based on the accuracy, precision, specificity, and sensitivity test. The receiver operating characteristic curve (ROC) was constructed, and the area under the ROC curve (AUC) was presented as the model's predictive ability. The Delong test checked the different feature sets in similar models (6). All statistical analysis procedures were performed using Python (version 3.7.6). The P < 0.05 was considered to be statistically significant.
Results
Clinical characteristics
Baseline characteristics and imaging CTA of 615 patients were initially obtained. Under the specified exclusion conditions, 112 who had carotid dissection and aneurysm, 214 who had intracranial vascular disease, 68 who had posterior circulation stroke, 26 who had a cerebral hemorrhage, and 16 who had an ischemic stroke caused by other etiologies were all excluded. The final study included 179 patients (mean age, 65.4 years; 69.9% males), and a total of 219 carotidarteries were included in the final analysis, of which 39 had a TIA (Figure 4).
Figure 4. Patient selection. The flowchart shows the process of enrolling patients in this study. The original population consisted of 6,314 subjects who had undergone carotid CTA studies, and we first screened the 615 subjects with mild carotid stenosis. According to the inclusion and exclusion criteria, we finally enrolled 179 subjects and included patients with TIA after CTA in the symptomatic group. CTA, computed tomography angiography; TIA, Transient ischemic attack.
Analysis of clinical features, radiomics features, and combined features
After we extract and select the feature information, the first nine features selected for the radiomics mode were maximum 3D diameter, major axis length, maximum 2D diameter (column), maximum 2D diameter (row), inverse variance, maximal correlation coefficient, gray level variance, strength, and voxel volume. The first two features selected for the clinical model were LDL and UA. The first three features selected for the combined model were maximum 3D diameter, LDL, and UA.
Predictive model evaluation
Overall, the RF model showed good predictive ability, followed by the XGBoost model, while the remaining three models had poor predictive performance (Figure 5). The top nine important features of radiomics data were selected, and the models were constructed and evaluated in the training and testing sets, respectively, which showed good diagnostic performance in the training and testing sets. The best model was RF (training set: AUC = 0.982, 95%CI, 0.973–0.999, ACC = 0.973, precision = 0.878, sensitivity = 0.968, and specificity = 0.865; testing set: AUC = 0.746, 95%CI, 0.761–0.979, ACC = 0.787, precision = 0.556, sensitivity = 0.417, and specificity = 0.926), and the performance of other models is shown in Table 2. In the machine learning models based on clinical features, the first two important features were LDL and UA, and the models were constructed and evaluated in the training and testing sets, respectively. The best model was XGBoost (training set: AUC = 0.948, 95%CI, 0.921–0.981, ACC = 0.936, precision = 0.828, sensitivity = 0.762, and specificity = 0.841; testing set: AUC = 0.765, 95%CI, 0.574–0.847, ACC = 0.727, precision = 0.429, sensitivity = 0.500, and specificity = 0.852). In the machine learning model based on combined features, the first three features were maximum 3D diameter, LDL, and UA, respectively, and the best model was RF (training set: AUC = 0.983, 95%CI, 0.977–0.998, ACC = 0.988, precision = 0.882, sensitivity = 0.952, and specificity = 0.873; testing set: AUC = 0.879, 95%CI, 0.787–0.979, ACC = 0.863, precision = 0.778, sensitivity = 0.583, and specificity = 0.963). In conclusion, the RF model constructed based on radiomics features and clinical features achieved the best predictive performance. We used the Delong test (goodness of fit test) to check whether different data were statistically different on the same model. In the RF model, the clinical information-based model was not statistically different from the radiomics-based model (AUC = 0.721 vs. AUC = 0.746, p = 0.813). However, the combined model outperformed the model based on clinical information (AUC = 0.879 vs. AUC = 0.721, p = 0.004) (Figure 6). In the XGBoost model, the models established by the three data types were not statistically different from each other (Figure 7).
Figure 5. ROC curves of various models. The area under the ROC curves (AUCs) of five different machine learning models constructed from radiomics features, clinical features, and combined features, respectively. (A) ROC curves of five machine learning models built with clinical features. (B) ROC curves of five machine learning models built with radiomics features. (C) ROC curves of five machine learning models built with combined features. ROC, receiver operating characteristic; AUC, area under curve; RF, random forest; XGBoost, eXtreme gradient boosting; LR, logistic regression; SVM, support vector machines; KNN, k-nearest neighbors.
Figure 6. Comparative analysis of ROC in the random forest model. Among the five models, the random forest model showed the best predictive performance. The random forest model was constructed with radiomics features, clinical features, and combined features, respectively, and the DeLong test was performed. The combined model was significantly better than the clinical model. (A) A DeLong test was used to compare the radiomics model and clinical model. P = 0.813. (B) A DeLong test was used to compare the radiomics model and the combined model. P = 0.098. (C) The DeLong test was used to compare the clinical model and the combined model. P = 0.004. ROC, receiver operating characteristic; RF, random forest.
Figure 7. Comparative analysis of ROC in the XGBoost model. In the XGBoost model, there were no statistical differences in the models based on clinical features, radiomics features, and combined features. (A) A DeLong test was used to compare the clinical model and the radiomics model. P = 0.598. (B) A DeLong test was used to compare the radiomics model and the combined model. P = 0.676. (C) The DeLong test was used to compare the clinical model and the combined model. P = 0.950. ROC, receiver operating characteristic; RF, random forest.
Discussion
In this study, patients with mild carotid stenosis (30%−50% NASCET) had a cumulative incidence of 17.8% of TIA on the ipsilateral side of the carotid artery after follow-up. This is slightly higher than previously reported results (10–15%) (3), possibly because living above 3,500 meters appears to be associated with a significantly increased risk of ischemic stroke and TIA, which may be related to polycythemia and other related factors such as increased blood viscosity (36).
Radiomics refers to the extraction of information from imaging images, including images that are difficult to identify by the human eye. Radiomics can capture tissue and lesion characteristics, thus, radiomics was first used to analyze tumor aggressiveness. One study predicted axillary lymph node metastasis in breast cancer. The XGBoost algorithm was used to establish a prediction model, and its AUC was 0.890 (20). Radiomics can also be used to predict clinical endpoints such as survival and treatment response. Nazari et al. predicted the mortality risk of patients with clear cell renal cell carcinoma within 5 years. They combined radiomics features with clinical data to establish a predictive model using four machine learning algorithms, of which the XGBoost model showed the best performance (AUC was in the range of 0.95~0.98) (37). Radiomics signatures can be mined, and in sufficiently large data, radiomics can be used to discover previously unknown markers and models of disease development, progression, and response to treatment, as well as new imaging studies for known diseases' classification. Liu et al. classified patients with rectal cancer and predicted the clinical outcomes of patients (38). Recent studies on radiomics and machine learning primarily focus on tumor classification or the prediction of clinical endpoints. According to the literature, there is no relevant research on the use of machine learning to build a model to predict the complications of patients with carotid artery stenosis of 30–50% based on CT imaging and clinical laboratory results; hence, the motivation for this study apply CT radiomics features combined with clinical features to predict the risk of TIA in patients with mild carotid stenosis.
In this study, we trained RF, XGBoost, LR, SVM, and KNN five radiomics group models on carotid artery CTA images to predict the occurrence of cerebrovascular symptoms in participants. When we evaluated our results in the testing set, an RF model based on radiomics and clinical information identified symptomatic participants with optimal diagnostic accuracy. Furthermore, we demonstrate that, although the radiomics-based RF model was not statistically significantly different when compared with the assessment of clinical characteristics, the combined information-based RF model could well distinguish between the high-risk and low-risk groups from participating (AUC = 0.879 vs. AUC = 0.721, p = 0.004).
In addition to predicting TIA in carotid artery stenosis, this study can also provide a valuable reference for formulating treatment plans through prediction. The ESVS issued a set of guidelines proposing medical therapy as an intervention for asymptomatic patients with stenosis <60% and symptomatic patients with stenosis <50%, and subsequently indicated that if symptoms persist, despite the absence of BMT, it may be reasonable to consider CEA/CAS (3). Therefore, there remains uncertainty in the treatment plan of patients with stenosis of 30–49%. If the patient predicts the risk of TIA in the future through the model, under the premise of BMT, CEA/CAS surgery can be given a higher priority to avoid complications. Therefore, we performed radiomics feature analysis of carotid arteries and their plaques, which have the potential ability to identify plaques leading to the risk of cerebrovascular events and have instructive value in guiding and optimizing the management of patients with mild stenosis. We aim to establish a model with optimal performance in predicting complications, combining both radiomics and clinical features. The models built together have the best predictive performance and are statistically different from the clinical features model. Traditional CTA imaging can only provide limited information on carotid plaque characteristics, and we need more objective methods that rely less on medical expertise while increasing prediction accuracy. Radiomics has the potential to be a tool that facilitates accurate phenotyping of abnormalities based on radiological images (39). A previous study explored that CT texture features could distinguish between symptomatic and asymptomatic patients, with AUCs ranging from 0.68 to 0.81 (40). The study used radiomics to acquire extensive data on stenotic carotid arteries and utilized machine learning (ML) models to improve diagnostic performance. The results showed that ML is beneficial for carotid CTA plaque analysis. However, the assessment of the carotid artery was not performed using clinical information added to the model. Another study on CTA used a machine learning model to identify symptomatic participants using radiology and routine assessment. Routine assessment refers to the analysis of carotid artery stenosis, plaque length, plaque thickness, and plaque ulceration, combined with imaging features, to identify symptomatic participants. The AUC was 0.85 (39). ML models based on radiomics and clinical features can be an ideal tool to identify which specific CTA plaque features are associated with increased TIA risk. Our study uses radiomics features and clinical laboratory tests (triglycerides, low-density lipoprotein, homocysteine, etc.) to assess the risk of complications based on machine learning models.
Although the findings are promising, some limitations must be addressed. The cases in this study were collected from a single-center. We included a small number of cases, and the incidence of cerebral ischemic events in patients with mild carotid stenosis was relatively low. The number of negative patients was much larger than that of positive patients, which may bias our results. To overcome overfitting our model, we calculated all diagnostic scores using 5-fold cross-validation. Prospective studies with larger sample sizes are still needed to assess long-term accuracy and stability. Manual segmentation may also affect parameter values, while automatic segmentation may reduce inter-observer variability and improve the feasibility of applying these methods on larger datasets.
Conclusion
In conclusion, we established five models for predicting ischemic cerebrovascular events in patients with mild carotid stenosis using radiomics, clinical data, and combined data based on machine learning models. After comparison, the random forest model of the combined data showed the best accuracy and maybe a more comprehensive and specific prediction model than previously reported methods. According to the predicted results, patients can have improved guided treatment plans using the proposed model.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by the Ethics Committee of The First Affiliated Hospital of Kunming Medical University (No. 2022-274). Written informed consent from the patients/participants or patients/participants legal guardian/next of kin was not required to participate in this study in accordance with the national legislation and the institutional requirements.
Author contributions
HX and LY: conception and design. YF: administrative support. WZ and JH: provision of study materials or patients. HX, YL, and YB: collection and assembly of data. HX, CZ, LY, and LZ: data analysis and interpretation. All authors: manuscript writing and final approval of manuscript.
Acknowledgments
We would like to thank all the clinicians for participating in this study.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Abbreviations
VOI, Volume of interest; NASCET, North American Symptomatic Carotid Endarterectomy Trial; RF, Random forest; XGBoost, EXtreme Gradient Boosting; LR, Logistic regression; SVM, Support vector machine; KNN, k-nearest neighbors; TIA, Transient ischemic attack; BMT, Best medical therapy; ESVS, European Society for Vascular Surgery; CEA, Carotid endarterectomy; CAS, Carotid artery stenting; ML, Machine learning.
References
1. Scheffler M, Pellaton A, Boto J, Barnaure I, Delattre BM, Remuinan J, et al. Hemorrhagic plaques in mild carotid stenosis: the risk of stroke. Can J Neurol Sci. (2021) 48:218–25. doi: 10.1017/cjn.2020.177
2. North American Symptomatic Carotid Endarterectomy Trial C, Barnett HJM, Taylor DW, Haynes RB, Sackett DL, Peerless SJ, et al. Beneficial Effect of Carotid Endarterectomy in Symptomatic Patients with High-Grade Carotid Stenosis. N Engl J Med. (1991) 325:445–53. doi: 10.1056/NEJM199108153250701
3. Naylor AR, Ricco JB, de Borst GJ, Debus S, de Haro J, Halliday A, et al. Editor's choice - management of atherosclerotic carotid and vertebral artery disease: 2017 clinical practice guidelines of the European society for vascular surgery (Esvs). Eur J Vasc Endovasc Surg. (2018) 55:3–81. doi: 10.1016/j.ejvs.2017.06.021
4. Warlow C. Mrc European carotid surgery trial: interim results for symptomatic patients with severe (70-99%) or with mild (0-29%) carotid stenosis. Lancet. (1991) 337:1235–43. doi: 10.1016/0140-6736(91)92916-P
5. Barnett HJ, Taylor DW, Eliasziw M, Fox AJ, Ferguson GG, Haynes RB, et al. Benefit of Carotid Endarterectomy in Patients with Symptomatic Moderate or Severe Stenosis. North American Symptomatic Carotid Endarterectomy Trial Collaborators. N Engl J Med. (1998) 339:1415–25. doi: 10.1056/NEJM199811123392002
6. Le EPV, Rundo L, Tarkin JM, Evans NR, Chowdhury MM, Coughlin PA, et al. Assessing robustness of carotid artery CT angiography radiomics in the identification of culprit lesions in cerebrovascular events. Sci Rep. (2021) 11:3499. doi: 10.1038/s41598-021-82760-w
7. Yan KK, Wang X, Lam WWT, Vardhanabhuti V, Lee AWM, Pang HH. Radiomics analysis using stability selection supervised component analysis for right-censored survival data. Comput Biol Med. (2020) 124:103959. doi: 10.1016/j.compbiomed.2020.103959
8. Gillies RJ, Kinahan PE, Hricak H. Radiomics: images are more than pictures, they are data. Radiology. (2016) 278:563–77. doi: 10.1148/radiol.2015151169
9. Lambin P, Leijenaar RTH, Deist TM, Peerlings J, de Jong EEC, van Timmeren J, et al. Radiomics: the bridge between medical imaging and personalized medicine. Nat Rev Clin Oncol. (2017) 14:749–62. doi: 10.1038/nrclinonc.2017.141
10. Vallieres M, Zwanenburg A, Badic B, Cheze Le Rest C, Visvikis D, Hatt M. Responsible radiomics research for faster clinical translation. J Nucl Med. (2018) 59:189–93. doi: 10.2967/jnumed.117.200501
11. Esteva A, Kuprel B, Novoa RA, Ko J, Swetter SM, Blau HM, et al. Dermatologist-level classification of skin cancer with deep neural networks. Nature. (2017) 542:115–8. doi: 10.1038/nature21056
12. Traverso A, Wee L, Dekker A, Gillies R. Repeatability and reproducibility of radiomic features: a systematic review. Int J Radiat Oncol Biol Phys. (2018) 102:1143–58. doi: 10.1016/j.ijrobp.2018.05.053
13. Schwier M, van Griethuysen J, Vangel MG, Pieper S, Peled S, Tempany C, et al. Repeatability of multiparametric prostate mri radiomics features. Sci Rep. (2019) 9:9441. doi: 10.1038/s41598-019-45766-z
14. Zhao B, Tan Y, Tsai WY, Qi J, Xie C, Lu L, et al. Reproducibility of radiomics for deciphering tumor phenotype with imaging. Sci Rep. (2016) 6:23428. doi: 10.1038/srep23428
15. Cikes M, Sanchez-Martinez S, Claggett B, Duchateau N, Piella G, Butakoff C, et al. Machine learning-based phenogrouping in heart failure to identify responders to cardiac resynchronization therapy. Eur J Heart Fail. (2019) 21:74–85. doi: 10.1002/ejhf.1333
16. Saba L, Mallarin G. Window settings for the study of calcified carotid plaques with multidetector CT angiography. AJNR Am J Neuroradiol. (2009) 30:1445–50. doi: 10.3174/ajnr.A1509
17. McKinney AM, Casey SO, Teksam M, Lucato LT, Smith M, Truwit CL, et al. Carotid bifurcation calcium and correlation with percent stenosis of the internal carotid artery on CT angiography. Neuroradiology. (2005) 47:1–9. doi: 10.1007/s00234-004-1301-4
18. Fedorov A, Beichel R, Kalpathy-Cramer J, Finet J, Fillion-Robin JC, Pujol S, et al. 3d slicer as an image computing platform for the quantitative imaging network. Magn Reson Imaging. (2012) 30:1323–41. doi: 10.1016/j.mri.2012.05.001
19. van Griethuysen JJM, Fedorov A, Parmar C, Hosny A, Aucoin N, Narayan V, et al. Computational radiomics system to decode the radiographic phenotype. Cancer Res. (2017) 77:e104–e7. doi: 10.1158/0008-5472.CAN-17-0339
20. Song BI. A Machine learning-based radiomics model for the prediction of axillary lymph-node metastasis in breast cancer. Breast Cancer. (2021) 28:664–71. doi: 10.1007/s12282-020-01202-z
21. Kamei K, Konta T, Hirayama A, Ichikawa K, Kubota I, Fujimoto S, et al. Associations between serum uric acid levels and the incidence of nonfatal stroke: a nationwide community-based cohort study. Clin Exp Nephrol. (2017) 21:497–503. doi: 10.1007/s10157-016-1311-7
22. Homocysteine Studies C. Homocysteine and risk of ischemic heart disease and stroke: a meta-analysis. JAMA. (2002) 288:2015–22. doi: 10.1001/jama.288.16.2015
23. Brunzell JD, Davidson M, Furberg CD, Goldberg RB, Howard BV, Stein JH, et al. Lipoprotein management in patients with cardiometabolic risk: consensus statement from the American Diabetes Association and the American College of Cardiology Foundation. Diabetes Care. (2008) 31:811–22. doi: 10.2337/dc08-9018
24. Drexel H, Amann FW, Beran J, Rentsch K, Candinas R, Muntwyler J, et al. Plasma triglycerides and three lipoprotein cholesterol fractions are independent predictors of the extent of coronary atherosclerosis. Circulation. (1994) 90:2230–5. doi: 10.1161/01.CIR.90.5.2230
25. Sabeti S, Exner M, Mlekusch W, Amighi J, Quehenberger P, Rumpold H, et al. Prognostic impact of fibrinogen in carotid atherosclerosis: nonspecific indicator of inflammation or independent predictor of disease progression? Stroke. (2005) 36:1400–4. doi: 10.1161/01.STR.0000169931.96670.fc
26. Casalnuovo G, Gerdts E, de Simone G, Izzo R, De Marco M, Giudice R, et al. Arterial stiffness is associated with carotid atherosclerosis in hypertensive patients (the Campania Salute Network). Am J Hypertens. (2012) 25:739–45. doi: 10.1038/ajh.2012.38
27. Kim DE, Lee KB, Jang IM, Roh H, Ahn MY, Lee J. Associations of cigarette smoking with intracranial atherosclerosis in the patients with acute ischemic stroke. Clin Neurol Neurosurg. (2012) 114:1243–7. doi: 10.1016/j.clineuro.2012.03.012
28. Alloubani A, Saleh A, Abdelhafiz I. Hypertension and diabetes mellitus as a predictive risk factors for stroke. Diabetes Metab Syndr. (2018) 12:577–84. doi: 10.1016/j.dsx.2018.03.009
29. Liu X, Zhang D, Liu Y, Sun X, Hou Y, Wang B, et al. A J-shaped relation of Bmi and stroke: systematic review and dose-response meta-analysis of 443 million participants. Nutr Metab Cardiovasc Dis. (2018) 28:1092–9. doi: 10.1016/j.numecd.2018.07.004
30. Kragsterman B, Bjorck M, Lindback J, Bergqvist D, Parsson H, Swedish Vascular R. Long-term survival after carotid endarterectomy for asymptomatic stenosis. Stroke. (2006) 37:2886–91. doi: 10.1161/01.STR.0000248967.44015.88
31. Iemolo F, Martiniuk A, Steinman DA, Spence JD. Sex differences in carotid plaque and stenosis. Stroke. (2004) 35:477–81. doi: 10.1161/01.STR.0000110981.96204.64
32. Herder M, Arntzen KA, Johnsen SH, Eggen AE, Mathiesen EB. Long-term use of lipid-lowering drugs slows progression of carotid atherosclerosis: the tromso study 1994 to 2008. Arterioscler Thromb Vasc Biol. (2013) 33:858–62. doi: 10.1161/ATVBAHA.112.300767
33. Erickson BJ, Korfiatis P, Akkus Z, Kline TL. Machine learning for medical imaging. Radiographics. (2017) 37:505–15. doi: 10.1148/rg.2017160130
34. Guo W, Li H, Zhu Y, Lan L, Yang S, Drukker K, et al. Prediction of clinical phenotypes in invasive breast carcinomas from the integration of radiomics and genomics data. J Med Imaging. (2015) 2:041007. doi: 10.1117/1.JMI.2.4.041007
35. Buckland ST, Efron B, Tibshirani RJ. An introduction to the bootstrap. Biometrics. (1994) 23:49–54. doi: 10.2307/2532810
36. Liu M, Yan M, Guo Y, Xie Z, Li R, Li J, et al. Acute ischemic stroke at high altitudes in China: early onset and severe manifestations. Cells. (2021) 10:809. doi: 10.3390/cells10040809
37. Nazari M, Shiri I, Zaidi H. Radiomics-based machine learning model to predict risk of death within 5-years in clear cell renal cell carcinoma patients. Comput Biol Med. (2021) 129:104135. doi: 10.1016/j.compbiomed.2020.104135
38. Liu Z, Zhang XY, Shi YJ, Wang L, Zhu HT, Tang Z, et al. Radiomics analysis for evaluation of pathological complete response to neoadjuvant chemoradiotherapy in locally advanced rectal cancer. Clin Cancer Res. (2017) 23:7253–62. doi: 10.1158/1078-0432.CCR-17-1038
39. Dong Z, Zhou C, Li H, Shi J, Liu J, Liu Q, et al. Radiomics vs. conventional assessment to identify symptomatic participants at carotid computed tomography angiography. Cerebrovasc Dis. (2022) 51:647–54. doi: 10.1159/000522058
Keywords: CT angiography, carotid artery, machine learning, prediction model, transient ischemic attack
Citation: Xia H, Yuan L, Zhao W, Zhang C, Zhao L, Hou J, Luan Y, Bi Y and Feng Y (2023) Predicting transient ischemic attack risk in patients with mild carotid stenosis using machine learning and CT radiomics. Front. Neurol. 14:1105616. doi: 10.3389/fneur.2023.1105616
Received: 22 November 2022; Accepted: 18 January 2023;
Published: 08 February 2023.
Edited by:
Siti Balqis Samdin, Xiamen University, MalaysiaReviewed by:
Siti Anom Ahmad, Putra Malaysia University, MalaysiaKazuya Matsuo, Japanese Red Cross Kobe Hospital, Japan
Copyright © 2023 Xia, Yuan, Zhao, Zhang, Zhao, Hou, Luan, Bi and Feng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yaoyu Feng, Znlhb3l1JiN4MDAwNDA7MTI2LmNvbQ==
†These authors have contributed equally to this work and share first authorship