Automated detection of early signs of irreversible ischemic change on CTA source images in patients with large vessel occlusion.

PURPOSE: To create and validate an automated pipeline for detection of early signs of irreversible ischemic change from admission CTA in patients with large vessel occlusion (LVO) stroke. METHODS: We retrospectively included 368 patients for training and 143 for external validation. All patients had anterior circulation LVO stroke, endovascular therapy with successful reperfusion, and follow-up diffusion-weighted imaging (DWI). We devised a pipeline to automatically segment Alberta Stroke Program Early CT Score (ASPECTS) regions and extracted their relative Hounsfield unit (rHU) values. We determined the optimal rHU cut points for prediction of final infarction in each ASPECT region, performed 10-fold cross-validation in the training set, and measured the performance via external validation in patients from another institute. We compared the model with an expert neuroradiologist for prediction of final infarct volume and poor functional outcome. RESULTS: We achieved a mean area under the receiver operating characteristic curve (AUC), accuracy, sensitivity, and specificity of 0.69±0.13, 0.69±0.09, 0.61±0.23, and 0.72±0.11 across all regions and folds in cross-validation. In the external validation cohort, we achieved a median [interquartile] AUC, accuracy, sensitivity, and specificity of 0.71 [0.68-0.72], 0.70 [0.68-0.73], 0.55 [0.50-0.63], and 0.74 [0.73-0.77], respectively. The rHU-based ASPECTS showed significant correlation with DWI-based ASPECTS (rS = 0.39, p<0.001) and final infarct volume (rS = -0.36, p<0.001). The AUC for predicting poor functional outcome was 0.66 (95%CI: 0.57-0.75). The predictive capabilities of rHU-based ASPECTS were not significantly different from the neuroradiologist's visual ASPECTS for either final infarct volume or functional outcome. CONCLUSIONS: Our study demonstrates the feasibility of an automated pipeline and predictive model based on relative HU attenuation of ASPECTS regions on baseline CTA and its non-inferior performance in predicting final infarction on post-stroke DWI compared to an expert human reader.

Radiomics-Based Prediction of Collateral Status from CT Angiography of Patients Following a Large Vessel Occlusion Stroke.

BACKGROUND: A major driver of individual variation in long-term outcomes following a large vessel occlusion (LVO) stroke is the degree of collateral arterial circulation. We aimed to develop and evaluate machine-learning models that quantify LVO collateral status using admission computed tomography angiography (CTA) radiomics. METHODS: We extracted 1116 radiomic features from the anterior circulation territories from admission CTAs of 600 patients experiencing an acute LVO stroke. We trained and validated multiple machine-learning models for the prediction of collateral status based on consensus from two neuroradiologists as ground truth. Models were first trained to predict (1) good vs. intermediate or poor, or (2) good vs. intermediate or poor collateral status. Then, model predictions were combined to determine a three-tier collateral score (good, intermediate, or poor). We used the receiver operating characteristics area under the curve (AUC) to evaluate prediction accuracy. RESULTS: We included 499 patients in training and 101 in an independent test cohort. The best-performing models achieved an averaged cross-validation AUC of 0.80 ± 0.05 for poor vs. intermediate/good collateral and 0.69 ± 0.05 for good vs. intermediate/poor, and AUC = 0.77 (0.67-0.87) and AUC = 0.78 (0.70-0.90) in the independent test cohort, respectively. The collateral scores predicted by the radiomics model were correlated with (rho = 0.45, p = 0.002) and were independent predictors of 3-month clinical outcome (p = 0.018) in the independent test cohort. CONCLUSIONS: Automated tools for the assessment of collateral status from admission CTA-such as the radiomics models described here-can generate clinically relevant and reproducible collateral scores to facilitate a timely treatment triage in patients experiencing an acute LVO stroke.

A deep learning analysis of stroke onset time prediction and comparison to DWI-FLAIR mismatch.

INTRODUCTION: When time since stroke onset is unknown, DWI-FLAIR mismatch rating is an established technique for patient stratification. A visible DWI lesion without corresponding parenchymal hyperintensity on FLAIR suggests time since onset of under 4.5 h and thus a potential benefit from intravenous thrombolysis. To improve accuracy and availability of the mismatch concept, deep learning might be able to augment human rating and support decision-making in these cases. METHODS: We used unprocessed DWI and coregistered FLAIR imaging data to train a deep learning model to predict dichotomized time since ischemic stroke onset. We analyzed the performance of Group Convolutional Neural Networks compared to other deep learning methods. Unlabeled imaging data was used for pre-training. Prediction performance of the best deep learning model was compared to the performance of four independent junior and senior raters. Additionally, in cases deemed indeterminable by human raters, model ratings were used to augment human performance. Post-hoc gradient-based explanations were analyzed to gain insights into model predictions. RESULTS: Our best predictive model performed comparably to human raters. Using model ratings in cases deemed indeterminable by human raters improved rating accuracy and interrater agreement for junior and senior ratings. Post-hoc explainability analyses showed that the model localized stroke lesions to derive predictions. DISCUSSION: Our analysis shows that deep learning based clinical decision support has the potential to improve the accessibility of the DWI-FLAIR mismatch concept by supporting patient stratification.

Personalised simulation of hemodynamics in cerebrovascular disease: lessons learned from a study of diagnostic accuracy.

Intracranial atherosclerotic disease (ICAD) poses a significant risk of subsequent stroke but current prevention strategies are limited. Mechanistic simulations of brain hemodynamics offer an alternative precision medicine approach by utilising individual patient characteristics. For clinical use, however, current simulation frameworks have insufficient validation. In this study, we performed the first quantitative validation of a simulation-based precision medicine framework to assess cerebral hemodynamics in patients with ICAD against clinical standard perfusion imaging. In a retrospective analysis, we used a 0-dimensional simulation model to detect brain areas that are hemodynamically vulnerable to subsequent stroke. The main outcome measures were sensitivity, specificity, and area under the receiver operating characteristics curve (ROC AUC) of the simulation to identify brain areas vulnerable to subsequent stroke as defined by quantitative measurements of relative mean transit time (relMTT) from dynamic susceptibility contrast MRI (DSC-MRI). In 68 subjects with unilateral stenosis >70% of the internal carotid artery (ICA) or middle cerebral artery (MCA), the sensitivity and specificity of the simulation were 0.65 and 0.67, respectively. The ROC AUC was 0.68. The low-to-moderate accuracy of the simulation may be attributed to assumptions of Newtonian blood flow, rigid vessel walls, and the use of time-of-flight MRI for geometric representation of subject vasculature. Future simulation approaches should focus on integrating additional patient data, increasing accessibility of precision medicine tools to clinicians, addressing disease burden disparities amongst different populations, and quantifying patient benefit. Our results underscore the need for further improvement of mechanistic simulations of brain hemodynamics to foster the translation of the technology to clinical practice.

Anatomical labeling of intracranial arteries with deep learning in patients with cerebrovascular disease.

Brain arteries are routinely imaged in the clinical setting by various modalities, e.g., time-of-flight magnetic resonance angiography (TOF-MRA). These imaging techniques have great potential for the diagnosis of cerebrovascular disease, disease progression, and response to treatment. Currently, however, only qualitative assessment is implemented in clinical applications, relying on visual inspection. While manual or semi-automated approaches for quantification exist, such solutions are impractical in the clinical setting as they are time-consuming, involve too many processing steps, and/or neglect image intensity information. In this study, we present a deep learning-based solution for the anatomical labeling of intracranial arteries that utilizes complete information from 3D TOF-MRA images. We adapted and trained a state-of-the-art multi-scale Unet architecture using imaging data of 242 patients with cerebrovascular disease to distinguish 24 arterial segments. The proposed model utilizes vessel-specific information as well as raw image intensity information, and can thus take tissue characteristics into account. Our method yielded a performance of 0.89 macro F1 and 0.90 balanced class accuracy (bAcc) in labeling aggregated segments and 0.80 macro F1 and 0.83 bAcc in labeling detailed arterial segments on average. In particular, a higher F1 score than 0.75 for most arteries of clinical interest for cerebrovascular disease was achieved, with higher than 0.90 F1 scores in the larger, main arteries. Due to minimal pre-processing, simple usability, and fast predictions, our method could be highly applicable in the clinical setting.

Dataset on acute stroke risk stratification from CT angiographic radiomics.

With advances in high-throughput image processing technologies and increasing availability of medical mega-data, the growing field of radiomics opened the door for quantitative analysis of medical images for prediction of clinically relevant information. One clinical area in which radiomics have proven useful is stroke neuroimaging, where rapid treatment triage is vital for patient outcomes and automated decision assistance tools have potential for significant clinical impact. Recent research, for example, has applied radiomics features extracted from CT angiography (CTA) images and a machine learning framework to facilitate risk-stratification in acute stroke. We here provide methodological guidelines and radiomics data supporting the referenced article "CT angiographic radiomics signature for risk-stratification in anterior large vessel occlusion stroke." The data were extracted from the stroke center registry at Yale New Haven Hospital between 1/1/2014 and 10/31/2020; and Geisinger Medical Center between 1/1/2016 and 12/31/2019. It includes detailed radiomics features of the anterior circulation territories on admission CTA scans in stroke patients with large vessel occlusion stroke who underwent thrombectomy. We also provide the methodological details of the analysis framework utilized for training, optimization, validation and external testing of the machine learning and feature selection algorithms. With the goal of advancing the feasibility and quality of radiomics-based analyses to improve patient care within and beyond the field of stroke, the provided data and methodological support can serve as a baseline for future studies applying radiomics algorithms to machine-learning frameworks, and allow for analysis and utilization of radiomics features extracted in this study.

CT angiographic radiomics signature for risk stratification in anterior large vessel occlusion stroke.

BACKGROUND AND PURPOSE: As "time is brain" in acute stroke triage, the need for automated prognostication tools continues to increase, particularly in rapidly expanding tele-stroke settings. We aimed to create an automated prognostication tool for anterior circulation large vessel occlusion (LVO) stroke based on admission CTA radiomics. METHODS: We automatically extracted 1116 radiomics features from the anterior circulation territory on admission CTAs of 829 acute LVO stroke patients who underwent mechanical thrombectomy in two academic centers. We trained, optimized, validated, and compared different machine-learning models to predict favorable outcome (modified Rankin Scale ≤ 2) at discharge and 3-month follow-up using four different input sets: "Radiomics", "Radiomics + Treatment" (radiomics, post-thrombectomy reperfusion grade, and intravenous thrombolysis), "Clinical + Treatment" (baseline clinical variables and treatment), and "Combined" (radiomics, treatment, and baseline clinical variables). RESULTS: For discharge outcome prediction, models were optimized/trained on n = 494 and tested on an independent cohort of n = 100 patients from Yale. Receiver operating characteristic analysis of the independent cohort showed no significant difference between best-performing Combined input models (area under the curve, AUC = 0.77) versus Radiomics + Treatment (AUC = 0.78, p = 0.78), Radiomics (AUC = 0.78, p = 0.55), or Clinical + Treatment (AUC = 0.77, p = 0.87) models. For 3-month outcome prediction, models were optimized/trained on n = 373 and tested on an independent cohort from Yale (n = 72), and an external cohort from Geisinger Medical Center (n = 232). In the independent cohort, there was no significant difference between Combined input models (AUC = 0.76) versus Radiomics + Treatment (AUC = 0.72, p = 0.39), Radiomics (AUC = 0.72, p = 0.39), or Clinical + Treatment (AUC = 76, p = 0.90) models; however, in the external cohort, the Combined model (AUC = 0.74) outperformed Radiomics + Treatment (AUC = 0.66, p < 0.001) and Radiomics (AUC = 0.68, p = 0.005) models for 3-month prediction. CONCLUSION: Machine-learning signatures of admission CTA radiomics can provide prognostic information in acute LVO stroke candidates for mechanical thrombectomy. Such objective and time-sensitive risk stratification can guide treatment decisions and facilitate tele-stroke assessment of patients. Particularly in the absence of reliable clinical information at the time of admission, models solely using radiomics features can provide a useful prognostication tool.

Similar admission NIHSS may represent larger tissue-at-risk in patients with right-sided versus left-sided large vessel occlusion.

BACKGROUND: We investigated the effects of the side of large vessel occlusion (LVO) on post-thrombectomy infarct volume and clinical outcome with regard to admission National Institutes of Health Stroke Scale (NIHSS) score. METHODS: We retrospectively identified patients with anterior LVO who received endovascular thrombectomy and follow-up MRI. Applying voxel-wise general linear models and multivariate analysis, we assessed the effects of occlusion side, admission NIHSS, and post-thrombectomy reperfusion (modified Thrombolysis in Cerebral Infarction, mTICI) on final infarct distribution and volume as well as discharge modified Rankin Scale (mRS) score. RESULTS: We included 469 patients, 254 with left-sided and 215 with right-sided LVO. Admission NIHSS was higher in those with left-sided LVO (median (IQR) 16 (10-22)) than in those with right-sided LVO (14 (8-16), p>0.001). In voxel-wise analysis, worse post-thrombectomy reperfusion, lower admission NIHSS score, and poor discharge outcome were associated with right-hemispheric infarct lesions. In multivariate analysis, right-sided LVO was an independent predictor of larger final infarct volume (p=0.003). There was a significant three-way interaction between admission stroke severity (based on NIHSS), LVO side, and mTICI with regard to final infarct volume (p=0.041). Specifically, in patients with moderate stroke (NIHSS 6-15), incomplete reperfusion (mTICI 0-2b) was associated with larger final infarct volume (p<0.001) and worse discharge outcome (p=0.02) in right-sided compared with left-sided LVO. CONCLUSIONS: When adjusted for admission NIHSS, worse post-thrombectomy reperfusion is associated with larger infarct volume and worse discharge outcome in right-sided versus left-sided LVO. This may represent larger tissue-at-risk in patients with right-sided LVO when applying admission NIHSS as a clinical biomarker for penumbra.

Image-to-image generative adversarial networks for synthesizing perfusion parameter maps from DSC-MR images in cerebrovascular disease.

Stroke is a major cause of death or disability. As imaging-based patient stratification improves acute stroke therapy, dynamic susceptibility contrast magnetic resonance imaging (DSC-MRI) is of major interest in image brain perfusion. However, expert-level perfusion maps require a manual or semi-manual post-processing by a medical expert making the procedure time-consuming and less-standardized. Modern machine learning methods such as generative adversarial networks (GANs) have the potential to automate the perfusion map generation on an expert level without manual validation. We propose a modified pix2pix GAN with a temporal component (temp-pix2pix-GAN) that generates perfusion maps in an end-to-end fashion. We train our model on perfusion maps infused with expert knowledge to encode it into the GANs. The performance was trained and evaluated using the structural similarity index measure (SSIM) on two datasets including patients with acute stroke and the steno-occlusive disease. Our temp-pix2pix architecture showed high performance on the acute stroke dataset for all perfusion maps (mean SSIM 0.92-0.99) and good performance on data including patients with the steno-occlusive disease (mean SSIM 0.84-0.99). While clinical validation is still necessary for future studies, our results mark an important step toward automated expert-level perfusion maps and thus fast patient stratification.

A precision medicine framework for personalized simulation of hemodynamics in cerebrovascular disease.

BACKGROUND: Cerebrovascular disease, in particular stroke, is a major public health challenge. An important biomarker is cerebral hemodynamics. To measure and quantify cerebral hemodynamics, however, only invasive, potentially harmful or time-to-treatment prolonging methods are available. RESULTS: We present a simulation-based approach which allows calculation of cerebral hemodynamics based on the patient-individual vessel configuration derived from structural vessel imaging. For this, we implemented a framework allowing segmentation and annotation of brain vessels from structural imaging followed by 0-dimensional lumped simulation modeling of cerebral hemodynamics. For annotation, a 3D-graphical user interface was implemented. For 0D-simulation, we used a modified nodal analysis, which was adapted for easy implementation by code. The simulation enables identification of areas vulnerable to stroke and simulation of changes due to different systemic blood pressures. Moreover, sensitivity analysis was implemented allowing the live simulation of changes to simulate procedures and disease progression. Beyond presentation of the framework, we demonstrated in an exploratory analysis in 67 patients that the simulation has a high specificity and low-to-moderate sensitivity to detect perfusion changes in classic perfusion imaging. CONCLUSIONS: The presented precision medicine approach using novel biomarkers has the potential to make the application of harmful and complex perfusion methods obsolete.

BRAVE-NET: Fully Automated Arterial Brain Vessel Segmentation in Patients With Cerebrovascular Disease.

Introduction: Arterial brain vessel assessment is crucial for the diagnostic process in patients with cerebrovascular disease. Non-invasive neuroimaging techniques, such as time-of-flight (TOF) magnetic resonance angiography (MRA) imaging are applied in the clinical routine to depict arteries. They are, however, only visually assessed. Fully automated vessel segmentation integrated into the clinical routine could facilitate the time-critical diagnosis of vessel abnormalities and might facilitate the identification of valuable biomarkers for cerebrovascular events. In the present work, we developed and validated a new deep learning model for vessel segmentation, coined BRAVE-NET, on a large aggregated dataset of patients with cerebrovascular diseases. Methods: BRAVE-NET is a multiscale 3-D convolutional neural network (CNN) model developed on a dataset of 264 patients from three different studies enrolling patients with cerebrovascular diseases. A context path, dually capturing high- and low-resolution volumes, and deep supervision were implemented. The BRAVE-NET model was compared to a baseline Unet model and variants with only context paths and deep supervision, respectively. The models were developed and validated using high-quality manual labels as ground truth. Next to precision and recall, the performance was assessed quantitatively by Dice coefficient (DSC); average Hausdorff distance (AVD); 95-percentile Hausdorff distance (95HD); and via visual qualitative rating. Results: The BRAVE-NET performance surpassed the other models for arterial brain vessel segmentation with a DSC = 0.931, AVD = 0.165, and 95HD = 29.153. The BRAVE-NET model was also the most resistant toward false labelings as revealed by the visual analysis. The performance improvement is primarily attributed to the integration of the multiscaling context path into the 3-D Unet and to a lesser extent to the deep supervision architectural component. Discussion: We present a new state-of-the-art of arterial brain vessel segmentation tailored to cerebrovascular pathology. We provide an extensive experimental validation of the model using a large aggregated dataset encompassing a large variability of cerebrovascular disease and an external set of healthy volunteers. The framework provides the technological foundation for improving the clinical workflow and can serve as a biomarker extraction tool in cerebrovascular diseases.