Next Article in Journal
Early Chronic Kidney Disease (G1-G3a) in Combination with Steatosis as a Predictor of Incident Ischemic Heart Disease: A Longitudinal Study in Non-Diabetic Koreans
Next Article in Special Issue
Overview of Acute Ischemic Stroke Evaluation and Management
Previous Article in Journal
Cardiovascular Benefits from Gliflozins: Effects on Endothelial Function
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biomedicines
Volume 9
Issue 10
10.3390/biomedicines9101357
biomedicines-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Treatment Efficacy Analysis in Acute Ischemic Stroke Patients Using In Silico Modeling Based on Machine Learning: A Proof-of-Principle

by
Anthony Winder
1,*,
Matthias Wilms
1,2,
Jens Fiehler
3 and
Nils D. Forkert
1,2,4,5
1
Department of Radiology, University of Calgary, Calgary, AB T2N 2T9, Canada
2
Hotchkiss Brain Institute, University of Calgary, Calgary, AB T2N 4N1, Canada
3
Department of Diagnostic and Interventional Neuroradiology, University Medical Center Hamburg-Eppendorf, 20251 Hamburg, Germany
4
Department of Clinical Neuroscience, University of Calgary, Calgary, AB T2N 2T9, Canada
5
Alberta Children’s Hospital Research Institute, University of Calgary, Calgary, AB T3B 6A8, Canada
*
Author to whom correspondence should be addressed.
Biomedicines 2021, 9(10), 1357; https://doi.org/10.3390/biomedicines9101357
Submission received: 20 July 2021 / Revised: 17 September 2021 / Accepted: 26 September 2021 / Published: 29 September 2021
(This article belongs to the Special Issue Stroke—Pathophysiology and New Therapeutic Strategies)
Browse Figures
Versions Notes

Abstract

:
Interventional neuroradiology is characterized by engineering- and experience-driven device development with design improvements every few months. However, clinical validation of these new devices requires lengthy and expensive randomized controlled trials. This contribution proposes a machine learning-based in silico study design to evaluate new devices more quickly with a small sample size. Acute diffusion- and perfusion-weighted MRI, segmented one-week follow-up imaging, and clinical variables were available for 90 acute ischemic stroke patients. Three treatment option-specific random forest models were trained to predict the one-week follow-up lesion segmentation for (1) patients successfully recanalized using intra-arterial mechanical thrombectomy, (2) patients successfully recanalized using intravenous thrombolysis, and (3) non-recanalizing patients as an analogue for conservative treatment for each patient in the sample, independent of the true group membership. A repeated-measures analysis of the three predicted follow-up lesions for each patient revealed significantly larger lesions for the non-recanalizing group compared to the successful intravenous thrombolysis treatment group, which in turn showed significantly larger lesions compared to the successful mechanical thrombectomy treatment group (p < 0.001). A groupwise comparison of the true follow-up lesions for the three treatment options showed the same trend but did not reach statistical significance (p = 0.19). We conclude that the proposed machine learning-based in silico trial design leads to clinically feasible results and can support new efficacy studies by providing additional power and potential early intermediate results.

1. Introduction

Stroke is a major cause of death and disability in the world. For example, the annual cost for stroke patients amounts to approximately CAD 2.8 billion in the Canadian health care system [1], USD 28 billion for the US health care system [2], and EUR 60 billion in the European health care systems [3]. Even a small improvement in clinical stroke outcomes could improve the daily life of many patients and result in significant cost savings in the health care system. The majority of all strokes are caused by an artery blockage due to a blood clot, the so-called ischemic stroke [4]. Such a blockage leads to a partial or complete restriction of blood flow to the brain tissue supplied by this artery and may result in irreversible damage to the brain cells in that region.
For many years, intravenous tissue plasminogen activator (IV tPA) therapy was the only acute ischemic stroke patient treatment together with stroke unit care. However, IV tPA is only effective in 30–40% of treated patients to achieve early recanalization [5]. Multiple randomized controlled trials have shown an overwhelming efficacy of intra-arterial therapy for treatment of large vessel occlusions with considerably higher recanalization rates compared to IV tPA alone [6]. Based on the results of these recent randomized controlled trials, new medical devices for intra-arterial therapy are currently being developed with the aim to increase recanalization rates. However, the actual clinical benefit of new devices must be shown in randomized control trials with large patient samples. For example, 500 patients (233 intra-arterial therapy, 267 IV tPA controls) were recruited in the MR CLEAN study [7], 315 patients (165 intra-arterial therapy, 150 IV tPA controls) in the ESCAPE study [8], and 196 patients (98 intra-arterial therapy, 98 IV tPA controls) in the SWIFT PRIME study [9], all of which compared intra-arterial therapy to IV tPA alone. Unlike pharmacological therapy trials, interventional neuroradiology is characterized by several iterative steps in engineering- and experience-driven device development with tentative design improvements every few months. Once a randomized controlled trial is finished and published, the next device generation might be already available. Additionally, the benefit of other novel acute ischemic stroke treatment approaches, including neuroprotectant drugs [10] or hypothermia [11], also needs to be tested in randomized controlled trials.
However, the smaller the expected effect size of a new treatment approach compared to the current standard of care is, the more patients need to be recruited into fully powered randomized controlled trials. This is not only delaying patients’ access to promising new technology or treatments but also increases the costs for future medical devices or drugs. Furthermore, a considerable number of patients in a randomized controlled trial might be treated with a potentially less effective or even harmful treatment compared to current standard of care. To overcome these limitations, in silico evaluation of new treatment devices or drugs has gained increasing interest in some fields outside of acute ischemic stroke for efficacy testing [12,13].
The term “in silico” refers to the usage of computer methods for understanding biological processes in the living organism. In silico validation of medical devices is typically performed using a computational biomechanical model [14]. However, no such model exists for the evaluation of stroke treatments due to the high complexity of the cerebrovascular system and interventional procedure, which renders an in silico validation, for example using computational fluid dynamics simulations, infeasible. Instead, we propose the use of machine learning models, which have been used previously to model the evolution of acute ischemic stroke under particular treatment conditions [15,16,17,18,19,20,21,22,23]. In these studies, real patient data were used to empirically and automatically optimize (train) the machine learning model for the purpose of predicting patients’ final voxel-wise tissue outcome given their baseline imaging and clinical information. In case such a model is trained using data from patients treated exclusively with a particular drug or device, it can be assumed that it practically models the effect of the given treatment approach on the evolution of the ischemic lesion between the patients’ baseline and follow-up imaging. Thus, a collection of models, each trained to capture the effect of a different treatment, could be used to simulate hypothetical tissue outcomes for multiple treatment approaches in a single patient. Comparing patients’ predicted tissue outcomes between different treatment approaches would then provide a means of in silico efficacy testing for medical devices and drugs in the context of acute ischemic stroke. This study design is presented graphically in Figure 1.
Compared to previous randomized controlled trials, such an in silico approach could significantly reduce the burden of patient recruitment. For instance, treatment-specific models can be trained using a few datasets of patients treated with the methods to be compared. After training, those models can be applied to historical datasets that meet the trial’s inclusion criteria, i.e., the same criteria that were used to select the few patients for the actual treatment. While state-of-the-art deep machine learning approaches for tissue outcome prediction have been shown to lead to good results, they are generally believed to require a large sample size for training [24] so that these models are of reduced interest for the proposed in silico trial design. Compared to this, random decision forests (RDF) have been shown to have relatively small sample size requirements for tissue outcome prediction [25] while still producing competitive tissue outcome prediction results compared to other conventional machine learning techniques, such as generalized linear models and k-nearest neighbors [21]. However, the general framework proposed in this paper can be used with any machine learning model. In addition, the ability to compare treatment-contingent predicted tissue outcomes within, rather than between, patients could reduce the effects of inter-patient variability on the analysis, increasing its statistical power. Overall, introducing in silico evaluation methods to test new acute ischemic stroke treatment approaches is likely to increase the speed and power of these trials while decreasing the cost. The aim of this study is to evaluate if random decision forests for tissue outcome prediction can be used for the proposed in silico trial design, and whether they convey the aforementioned benefits.

Article Contributions

The main contributions of this paper can be summarized as follows:
  • An experimental framework is proposed for predicting the differential treatment effects of current gold-standard medical interventions using machine learning methods.
  • The proposed in silico study setup is validated using a large database of real stroke patient data.
  • It is demonstrated that treatment efficacy analysis using the proposed in silico framework benefits from an increase in statistical power.

2. Materials and Methods

2.1. Patient Data and Image Acquisition

The patient data used in this retrospective study, including information about image acquisition protocols and preprocessing methods utilized, have been described in detail in a previous publication [21]. These are summarized in the following for the reader’s ease of reference. Briefly described, this study utilized the anonymized data of 90 acute ischemic stroke patients with a major vessel occlusion of the ICA or MCA M1 segment treated at the University Medical Center Hamburg-Eppendorf, Germany. All patients used in this study were selected so that they fulfill typical guidelines for intra-arterial treatment. The patients were treated with an intravenous thrombolytic agent either exclusively (IV, n = 48) or in conjunction with intra-arterial mechanical thrombectomy (IA, n = 42). All study protocols and procedures were conducted in accordance with the ethical guidelines (Ethics committee of the Hamburg Chamber of Physicians, Hamburg, Germany) and in compliance with the Declaration of Helsinki. As this was a retrospective analysis, the requirement of informed consent was waived.
In all cases, the datasets available for each patient included diffusion-weighted (DWI) and perfusion-weighted (PWI) magnetic resonance imaging (MRI) acquired at admission, follow-up MRI or CT imaging acquired 5–7 days after stroke onset, and the following clinical variables: NIHSS score and age at admission, sex, time from symptom onset to imaging, and a binarized measure of recanalization status (TICI ≥ 2b). PWI and DWI datasets were acquired using a 1.5T Sonata or Avanto scanner (Siemens, Erlangen, Germany). Datasets were first divided by binarized recanalization status and then, for successfully recanalizing datasets, by treatment type to yield the following three patient groups: IA-recanalized (IAR, n = 33), IV-recanalized (IVR, n = 23), and non-recanalized (NR, n = 34). For this feasibility analysis, we used the non-recanalized patient group as an analogue for conservative treatment; recanalization status, being unavailable in the acute stage of stroke, is neither used for training nor prediction in the models.

2.2. Summary of Training Features

Tissue viability in acute ischemic stroke is a complex phenomenon that has a diverse set of contributing factors related to cerebral hemodynamics, patient-specific risks, differences in tissue vulnerability, and medical intervention. We selected features for model training with consideration for each of these factors as well as the design of prior tissue outcome prediction studies. Additional factors affecting stroke outcome are constantly being discovered [26]. However, for this study, we limited the data used for model training to that generally available for all patients in a retrospectively accessed database of clinical stroke patients.

2.2.1. Voxel-Wise Medical Imaging Features

Biologically relevant hemodynamic changes in the human brain due to a stroke are typically observed and measured in the clinical setting using medical imaging techniques such as multi-modal MRI or CT to quantify their combined effects on a voxel-level scale. These voxel-level imaging parameters have historically been used to determine tissue viability using a threshold-based approach that segments the affected tissue into three compartments: (1) the ischemic core, which is already infarcted and cannot be salvaged, (2) the ischemic penumbra, which will eventually develop infarction without reperfusion, and (3) the benign oligemia, which is hypoperfused but not at risk of infarction [27]. These imaging parameters are considered sufficiently reliable that stroke treatment decisions are often based on the volume of the ischemic core and relative volume of the ischemic penumbra [28]. However, there is a growing acknowledgement that these imaging parameters are an oversimplification of the complex tissue perfusion and do not allow the making of true personalized outcome predictions [29].
Methods using machine learning to better interpret imaging parameters for stroke outcome prediction have already been proposed in a number of independent studies [15,16,17,18,19,20,21,22,23] and international challenges [30]. These methods typically consider the full range of imaging features previously used in clinical practice as salient predictors of tissue outcome. For perfusion-weighted imaging, these parameters are cerebral blood flow (CBF), cerebral blood volume (CBV), mean transit time (MTT), and time-to-maximum (Tmax). For diffusion-weighted MRI, this is the parameter apparent diffusion coefficient (ADC). All of these voxel-wise imaging parameters are included in this work in the machine learning model.

2.2.2. Patient-Level Clinical Features

The extent to which imaging parameters predict poor tissue viability is likely mediated by several patient-specific factors. Infarction in the late hyperacute stage, for example, has been shown to depend not only on the magnitude of the hypoperfusion as indicated on imaging, but also the duration of the hypoperfusion [18]. Specifically, in patients who are rapidly recanalized, criteria based solely on imaging information may substantially overpredict the extent of the final infarct [31]. For this reason, we included the time from symptom onset to imaging as a training feature. Similarly, the mean values of perfusion and diffusion imaging parameters in healthy tissue are known to vary with respect to patient age [32,33] and sex [32]. Thus, both parameters are included as additional training features. Finally, we also included the NIHSS score at admission as it is a strong independent predictor of tissue fate in similar stroke outcome prediction models [17,18].

2.2.3. Tissue-Specific and Spatial Features

Even within an individual brain, the infarction thresholds for perfusion and diffusion imaging parameters are known to vary between white and gray brain matter [34], between different brain regions [35] and, due to the spatial organization of stroke, likely with proximity to the ischemic core. Therefore, we included atlas-based white matter probabilities and brain region labels as training features. Because our model processes individual voxels without spatial context, we included the Euclidean distance of each voxel to the closest point in the ischemic core, as determined by ADC thresholding, as our final training feature.

2.2.4. Medical Intervention

Although the treatment administered to a patient has a profound effect on their tissue outcome, we opted not to include the treatment as a training feature and to instead train a different model for each successful treatment of interest, as well as unsuccessful treatment as a surrogate for conservative treatment. The advantage of doing so is twofold. First, it enables the use of different training datasets for each model so that a model meant to learn one treatment effect does not learn from patients administered other treatments. Second, a model for a new treatment can be added to the ensemble or an existing model can be tweaked without having to retrain any of the other models.

2.3. Image Preprocessing and Feature Extraction

All imaging features described above were extracted for the patients available in this study using the AnToNIa software tool [36]. Briefly described, DWI acquired with diffusion weightings of b = 0 and b = 1000 s/mm2 were used to calculate a quantitative apparent diffusion coefficient (ADC) map, which was then used to segment the brain tissue, cerebrospinal fluid (ADC > 1200 × 10−6 mm²/s), and the ischemic core (ADC < 550 × 10−6 mm2/s). Subsequently, the ischemic core segmentation was used to compute a map of the shortest Euclidean distance from each voxel to the ischemic core. For each PWI dataset, maps of cerebral blood volume (CBV), cerebral blood flow (CBF), mean transit time (MTT), and time-to-maximum of the residual curve (Tmax) were calculated using a block-circulant singular value decomposition with a threshold of 0.15 [37]. The arterial input function needed for this was automatically identified using an atlas-based approach. For each follow-up image, an experienced medical expert used AnToNIa to manually segment the final lesion, yielding a binary mask of the real tissue outcome. Finally, image registration was used to transfer each perfusion parameter map (CBV, CBF, MTT, Tmax), the tissue type information (voxel-wise white and gray matter probabilities) and region labels from the MNI atlas, and the follow-up real tissue outcome to the patient’s ADC space. Clinical variables were associated with each voxel of the ADC image by simple repetition. As a result, each voxel was associated with a total of 12 features: ADC, distance to ischemic core, tissue type, anatomical location, CBV, MTT, Tmax, CBF, NIHSS, age, sex, and time from symptom onset, whereas the real tissue outcome was used as the outcome variable for training and testing of the machine learning models.

2.4. Machine Learning Model

Based on the preprocessed voxel-wise data, separate optimized random forest classifier models were trained to predict the voxel-wise treatment outcomes of successful IA therapy, successful IV therapy, and unsuccessful reperfusion therapy using datasets from the IAR, IVR, and NR groups, respectively. Training a separate prediction model for each treatment group, rather than a single combined model using the cohort information as a patient-level variable, guarantees that the effect of each treatment is learned independently of the other treatments being evaluated. In a combined model, for example, datasets from the IAR group would influence the prediction of treatment outcomes under IVR conditions, which could introduce significant bias.
Each classifier comprised 100 decision trees and was implemented using ALGLIB (www.alglib.net, accessed on 15 August 2021). For each classifier, the single hyperparameter r, representing the proportion of the training set used for the construction of individual trees, was set to 0.5. Stratified random undersampling was used to account for class imbalances resulting from the relatively small volume of the final infarct compared to the whole ipsilateral hemisphere. For each dataset, this entailed sampling every voxel that had a real infarct tissue outcome as well as an equal number of randomly selected non-infarct voxels from the ipsilateral brain parenchyma. A more granular description of this machine learning model has been previously published [21].
To virtually increase the sample size following the in silico principle described in the introduction, all three tissue outcomes (IAR, IVR, NR) were predicted for every patient available for this feasibility study. In all cases where the predicted treatment outcome matched the true treatment, which is referred to as ‘standard prediction’ in the following, the predicted and true tissue outcomes were compared using the Dice metric as a means of model validation. In all other cases, which are referred to as ‘cross-predictions’ in the following, the predicted tissue outcome represented an estimate of how the tissue outcome would look like if the patient had received the different treatment. In case of the standard predictions, a leave-one-out approach was used to maximize the quantity of datasets available. Essentially, each model was trained from scratch before every standard prediction using all datasets from the corresponding group except for the dataset for which the outcome was being predicted. For all cross-predictions, each model was trained only once using all data from one group without any omissions, since the dataset being predicted was not in the data partition used for training. For both the standard and cross-predictions, the predicted maps were postprocessed using morphological closing with a 3 × 3 × 3 kernel followed by connected-component analysis excluding components less than ten voxels in size as a means of noise reduction.

2.5. Statistical Evaluation

To test for significant differences between the predicted tissue outcomes simulating successful IA therapy, successful IV therapy, and unsuccessful reperfusion therapy, the predicted lesion volumes for all 90 patients were compared using a one-way repeated-measures ANOVA. In case of significant results, planned contrasts were performed between each pair of models using a Bonferroni-adjusted p-value to achieve a family-wise error rate of 5%. To demonstrate the difference between this analysis and one which does not use tissue outcome prediction models, a one-way ANOVA was conducted comparing the ground-truth lesion volumes between the corresponding patient groups (IAR, IVR, NR). A log transformation was applied to both the true and the predicted lesion volumes prior to conducting the respective ANOVAs to obtain variables that are normally distributed, as assessed by the Shapiro–Wilk test.
Confusion matrices and ROC analyses for the models were computed after applying the brain tissue masks derived from the ADC maps to exclude background and CSF voxels. Therefore, true positives were counted as the number of voxels positively labelled in all three of the brain tissue, predicted infarct, and ground-truth infarct segmentations. Conversely, true negatives were counted as the number of voxels positively labelled in the brain mask, but negatively labelled in both the predicted and ground-truth infarct segmentations. False positives and negatives were counted according to similar definitions, only considering voxels where the labels on the ground-truth infarct segmentation differed from the predicted infarct. Accuracy, sensitivity, specificity, and precision were derived from these four quantities according to standard definitions. Separate statistics were computed for each patient, then averaged together to avoid a heavy bias in the reported values toward larger lesions. AUC values were derived from the ROC curves using Simpson’s method.
The Dice metric, which is commonly used to quantify the overlap between two binary segmentations, is reported as a measure of the agreement between each model’s standard predictions and the corresponding ground-truth lesion segmentations. These Dice values were compared using a one-way ANOVA to determine if any model was significantly more reliable than another. Finally, the volume error of each predicted lesion was computed as a proportion of the true lesion volume and statistically compared between each pair of models using independent samples t-tests. All statistical analyses were performed in SPSS25.

3. Results

An in-depth description of the patient and imaging characteristics has been previously published [21]. Briefly summarized, there was no significant difference between any of the IAR, IVR, and NR patient groups with respect to the acute stroke lesion volume, penumbral tissue volume, follow-up lesion volume, age, sex distribution, or NIHSS scores at admission. However, the NR patients showed a significantly longer onset-to-imaging time compared to the IVR patients (p = 0.03).
The repeated-measures analysis, made possible by the proposed in silico trial design, revealed that the three sets of predicted lesion volumes were significantly different (p < 0.001). Planned contrasts indicated that the NR model predicted on average significantly larger lesion volumes (M = 68.75 mL, SD = 66.83 mL) compared to the IVR model (M = 49.67 mL, SD = 59.57 mL, p < 0.001), which in turn produced significantly larger lesion volumes than the IAR model (M = 34.42 mL, SD = 49.09, p < 0.001). Accordingly, the differences between the NR model and the IAR model (p < 0.001) were also significant. Exemplary tissue outcome predictions applying the three different models to the same patients are shown in Figure 2.
In contrast to the findings comparing the predicted lesion volumes, there was no significant difference comparing the real final lesion volumes between the three groups (p = 0.19). However, there was a non-significant trend between the three groups that followed the same pattern identified for the predicted lesion volumes, showing the smallest lesion volumes for IAR, slightly larger lesion volumes for IVR, and the largest lesion volumes for the NR group. The distribution of the real and predicted lesion volumes, and the similarities therein, can be seen in Figure 3.
The mean Dice values comparing the standard predictions to their corresponding ground-truths were not significantly different between the IAR, IVR, and NR machine learning models (p = 0.90). All three models produced mean Dice values > 0.44, which is comparable to similar tissue outcome prediction studies in the field of acute ischemic stroke [15,16,17,18,19,20,21,22,23]. Additionally, the proportional error, which quantifies the difference between the predicted and true lesion volume as a proportion of the true lesion volume, did not differ significantly between models (p = 0.60). These statistics can be found in Table 1.
Additional statistics including the mean accuracy, sensitivity, specificity, and precision of the models operating at their optimal thresholds are reported in Table 2, while Figure 4 shows receiver operating characteristic (ROC) curves derived from each model operating at thresholds varying from 0 to 1 by intervals of 0.01. The area under the ROC curve (AUC) for the IAR, IVR, and NR models were 0.979, 0.976, and 0.957, respectively. While these additional metrics are provided in keeping with the conventions of other machine learning studies in the field of stroke medicine [16,18,19,22,23,38], it should be noted that their interpretation is complicated by the fact that voxel-based tissue outcome prediction is a highly imbalanced problem. That is to say, the vast majority of brain tissue voxels will not proceed to infarction and are easy to identify as such, leading to especially high values for accuracy, specificity, and the area under the ROC curve (AUC), which may not intuitively reflect the ability of the model to discriminate between healthy and ischemic voxels.

4. Discussion

The main finding of this study was that the proposed in silico trial design leads to clinically feasible results, with successful recanalization following IA therapy resulting in significantly smaller predicted lesion volumes than successful recanalization following IV therapy, which in turn resulted in significantly smaller predicted lesion volumes than non-successful or conservative treatment. This trend is already well-known from previous studies comparing intra-arterial therapy [6] and intravenous thrombolysis [39], which confirms that the RDF models were able to train successfully on the relatively small dataset used for this study and make plausible treatment-specific tissue outcome predictions. By contrast, a direct comparison of the ground-truth lesion volumes for the three treatment outcomes yielded the same trend at a statistically insignificant level. This difference in statistical significance demonstrates that the proposed in silico trial design can be used to reduce the minimum sample size of new trials or provide early intermediate results that might be helpful to identify potentially unsuccessful trials at an earlier time point.
A major contributing factor to the significance of our results was likely the paired and repeated-measures statistical tests made possible by the in silico trial design, which offer several advantages over a conventional independent samples test. First, the paired and repeated-measures tests operate on each patient’s multiple predicted tissue outcomes, (as opposed to their single ground-truth tissue outcome), which represents a virtual expansion of the study sample size. Second, these tests compare the multiple predicted tissue outcomes within, rather than between, individual patients, which serves to control for the between-subjects variability and further increase the statistical power of the analysis.
It should also be noted that all patients used in this study were selected so that they fulfilled typical guidelines for intra-arterial treatment. The patients in the NR group only showed a significantly longer symptom onset-to-imaging time than patients in the IVR group while all other clinical characteristics were similar. This suggests that the ability to predict treatment effects was learned by the machine learning models and not simply inherited from bias in the training or testing data. The absence of significant differences in the proportional volume error of the models further suggests that the results are not biased by a systematic underprediction of lesion volumes by the IAR model nor an overprediction of lesion volumes by the NR model. The quality of the standard predictions from the three models, as evaluated by the Dice coefficient, was also found to be not significantly different. Taken together, it is reasonable to conclude that the random forest models were successfully able to learn the differential treatment effects of successful IA therapy, successful IV therapy, and unsuccessful therapy used as a proxy for conservative treatment. Within this context, it should be noted that a much higher reperfusion/recanalization rate is typically observed for IAR compared to IVR treatment. Since the IAR and IVR treatment machine learning models developed in this work both assume successful treatment in all cases (as only those patients were used for training), the effects identified are rather an underestimation of the true effects that would be found when comparing the two clinical treatment cohorts without taking the recanalization success into account, which gives even more credibility to the ability of the suggested in silico modeling approach to being able to determine significant treatment effects even with a small sample.
The proposed in silico evaluation scheme might be desirable for research and clinical studies in general, where patient recruitment can be costly and labor-intensive, but is especially relevant in situations where it is difficult to recruit a large number of patients or to match patients between groups. For example, clinical trials are warranted to examine the prospective benefit of mechanical thrombectomy in less-common patient demographics such as those with occlusion of the posterior circulation [40,41], as well as children and adolescents [42,43]. For these demographics in particular, the limited availability of patient datasets may hinder efforts to construct a sufficiently powered regular study. The proposed method would allow researchers to acquire only a few datasets for model training, then apply these trained models to the large number of historical datasets available in most primary stroke care centers for the in silico evaluation. In silico outcome prediction may also have clinical utility in the field of computer-aided decision making. For example, comparing the predicted lesions for new stroke patients under different treatment conditions could assist clinicians in determining the most appropriate treatment on an individual basis, though this idea remains to be evaluated in future studies.
It should be mentioned that this is not the first study to use predictive modeling approaches for an in silico evaluation [44]. Fiehler et al. used a similar machine learning model based on CTP datasets to evaluate the effectiveness of a new recanalization device together with a specific intermediate catheter compared to intravenous tPA [15]. The major difference is that, in the cited clinical study, only a single model was trained (using data from patients treated with tPA). This trained machine learning model was then applied to the patients treated with the new recanalization device and the volume of saved tissue was quantified by comparing the true lesion volume in these patients with the predicted tPA treatment lesion volume. Compared to this, the proposed method of training multiple machine learning models has several theoretical benefits. First, one can make full use of many historical datasets available for the actual in silico evaluation as the real lesion outcome is not required. Second, it is statistically more sound to compare outcomes that were generated using the same method. Comparing a true lesion outcome with the predicted lesion outcome might be biased, for example, due to regularization used within the machine learning models. Additionally, blood clots can move downstream or break up into smaller clots because of the treatment so that the anatomical lesion location varies between baseline and follow-up imaging, which cannot be accounted for when comparing true and virtual lesion outcomes. While this is less of a problem for the direct comparison of two lesion predictions, it would be important to exclude these cases from the model training.
There is also a much larger body of knowledge regarding predictive modeling for stroke outcome prediction outside the context of in silico treatment evaluation. Many of the earliest machine learning models trained for stroke tissue outcome prediction were based on linear [17,22] or logistic [18] regression. While these methods were useful for quantifying the contribution of particular variables to the overall infarct probabilities [17], more recent publications comparing regression models and decision forests unanimously agree that decision forests facilitate more accurate tissue outcome predictions [16,19,21]. This was the motivating factor behind using decision forests in the current work. However, it should be noted that the purpose of this investigation was not to propose the next state-of-the-art model for tissue outcome prediction, but rather to propose an in silico framework for treatment evaluation that can be used with any tissue outcome prediction model. More recently, there are more works using deep learning models to predict stroke tissue outcomes from imaging and clinical data showing high accuracy [23]. Following trends in other medical domains [24], it seems likely that deep learning models will significantly outperform conventional machine learning models for stroke tissue outcome prediction. However, this remains to be shown via comparative studies. In this context, the results of the current study have a couple of important implications. First, we show that deep learning is not strictly necessary for the proposed in silico approach. This may be especially significant for studies where data are scarce or difficult to obtain, considering the relatively small sample size requirements of decision forests [25]. Conversely, deep learning methods are uniquely suited to leverage the large amount of data present in electronic health records, otherwise known as clinical big data [45]. Previous predictive modeling studies have used a few additional variables including age [15,17], sex [15,17], admission NIHSS [15,17,18], and tissue priors [16,17] to enrich the training data of their stroke tissue outcome models. The current study includes all of these as training features. However, if more salient variables are discovered, it may become less feasible to use conventional (non-deep) classifiers. In this eventuality, the proposed in silico method could still be used alongside deep learning classifiers to take advantage of big medical data. Second, we provide a benchmark method for the evaluation of future deep learning classifiers. Demonstrating the prospective benefit of deep learning in the context of in silico treatment evaluation requires a standardized, accessible, and well-validated comparison method. Random forests are already implemented in a number of freely available software packages and have been repeatedly validated as high-performing stroke tissue outcomes in the existing literature [15,16,19,21]. As novel deep learning models (and other future methods) are explored as a way to improve the performance of the proposed in silico method, we expect this study to provide a valuable point of reference.
There are some limitations of this work that have to be mentioned. First, large clinical stroke trials typically use patients’ disability and functional independence score on the modified Rankin Scale (mRS) as the primary endpoint because the mRS can be evaluated consistently across trials and conveys no additional risk to the patient [46]. In comparison, the current study analyzes patients’ follow-up lesion volume, which is a common secondary endpoint in clinical stroke trials. However, there is a strong inverse correlation between the likelihood of a good functional outcome (mRS 0–2) and the volume of the patient’s follow-up lesion volume [47,48]. It has been reported that this correlation also holds true for follow-up lesions predicted from acute PWI data [49], raising the possibility that mRS scores could be predicted directly from acute image data and clinical information in future and used in the same manner as the prediction lesion volume for an in silico evaluation. As mRS is used as both a clinical and a functional outcome, an imaging-based approach to predicting mRS would likely benefit from the availability of information about the infarct location, which adds valuable context when interpreting the functional severity of stroke from the lesion volume [50], as well as the potential to add salient covariates [51] and additional neurological data [38] to the model. More recently, it has been suggested that the post-treatment NIHSS score may be a suitable alternative primary outcome measure to mRS for clinical trials of acute ischemic stroke [52]. One compelling argument for the incorporation of post-treatment NIHSS into future predictive models is that especially high or low post-treatment NIHSS scores are associated with discrepancies between patients’ post-treatment infarct volume and 90-day mRS [53]. Beyond using post-treatment NIHSS as a surrogate for mRS, there is the potential to use the true or predicted post-treatment NIHSS alongside the voxel-wise tissue fate as each appears to provide context for the other in determining the patient’s true functional outcome. Generally, it would be possible to include any kind of clinical or functional outcome score into the framework as long as it can be predicted with reasonable accuracy from the baseline imaging data or predicted lesion outcome.
A second limitation of this study is that it employed a leave-one-out cross-validation approach to a single-center dataset, making a full discussion of the study’s external validity difficult. Within this context, it should be noted that machine learning models that were trained using datasets from a single institution should not be applied to a multi-center database of historical datasets without caution, as differences due to image acquisition or patient selection might bias the machine learning models. In this case, it would be more appropriate to also train the machine learning models using datasets from the same centers where the test datasets were acquired. More generally, care must be taken to ensure that the training data for a model is representative of the data that it is applied to. An insufficient quantity of training data or the presence of significant groupwise differences between the training and test datasets could bias the results of an in silico experiment. Investigating a novel treatment device, patient demographic, or complication, for example, would require that a new model be trained using a sample of the novel data. Finally, it should be noted that the proposed framework still needs be validated in more detail using prospectively collected trial data, such as a multi-center, double-blind, randomized controlled ESCAPE-NA1 trial [54] or other recent studies.

5. Conclusions

The results of this study suggest that machine learning models can be used for an efficient in silico trial design that leads to clinically plausible results and could be used for evaluating the efficacy of new treatment devices or approaches for acute ischemic stroke patients. Even though the proposed in silico study design has the potential to significantly accelerate acute stroke studies and reduce sample sizes, the proposed method should be rather considered as a valuable tool to support clinical trials (e.g., early intermediate efficacy analyses) rather than replace them entirely at this point as this approach still needs to be evaluated and validated in more detail.

Author Contributions

Conceptualization, A.W., M.W., and N.D.F.; methodology, A.W. and N.D.F.; software, N.D.F.; investigation, A.W. and N.D.F.; data curation, J.F.; writing—original draft preparation, A.W. and N.D.F.; writing—review and editing, A.W., M.W., J.F., and N.D.F.; supervision, N.D.F. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the Heart and Stroke Foundation of Canada Grant in aid (G-17-0018368), the Canada Research Chairs program, as well as by the River Fund at Calgary Foundation.

Institutional Review Board Statement

All study protocols and procedures were conducted in accordance with the ethical guidelines (Ethics committee of the Hamburg Chamber of Physicians, Hamburg, Germany) and in compliance with the Declaration of Helsinki. As this was a retrospective analysis, the requirement of informed consent was waived.

Informed Consent Statement

As this was a retrospective analysis, the requirement of informed consent was waived.

Data Availability Statement

The datasets generated during and/or analyzed during the current study are not publicly available due to containing information that could compromise the privacy of research participants but are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mittmann, N.; Seung, S.J.; Hill, M.D.; Phillips, S.J.; Hachinski, V.; Coté, R.; Buck, B.H.; Mackey, A.; Gladstone, D.J.; Howse, D.C.; et al. Impact of disability status on ischemic stroke costs in Canada in the first year. Can. J. Neurol. Sci. 2012, 39, 793–800. [Google Scholar] [CrossRef] [Green Version]
  2. Benjamin, E.J.; Muntner, P.; Alonso, A.; Bittencourt, M.S.; Callaway, C.W.; Carson, A.P.; Chamberlain, A.M.; Chang, A.R.; Cheng, S.; Das, S.R.; et al. Heart disease and stroke statistics—2019 update: A report from the American heart association. Circulation 2019, 139, e56–e528. [Google Scholar] [CrossRef] [PubMed]
  3. Luengo-Fernandez, R.; Violato, M.; Candio, P.; Leal, J. Economic burden of stroke across Europe: A population-based cost analysis. Eur. Stroke J. 2020, 5, 17–25. [Google Scholar] [CrossRef] [PubMed]
  4. Donnan, G.A.; Fisher, M.; Macleod, M.; Davis, S.M. Stroke. Lancet 2008, 371, 1612–1623. [Google Scholar] [CrossRef]
  5. Mijajlovic, M.D. Thrombolytic or endovascular therapy for acute ischemic stroke: Time is brain. J. Neurosci. Rural. Pr. 2014, 5, 3–5. [Google Scholar] [CrossRef]
  6. Goyal, M.; Menon, B.K.; van Zwam, W.H.; Dippel, D.W.J.; Mitchell, P.J.; Demchuk, A.M.; Dávalos, A.; Majoie, C.B.L.M.; van der Lugt, A.; de Miquel, M.A.; et al. Endovascular thrombectomy after large-vessel ischaemic stroke: A meta-analysis of individual patient data from five randomised trials. Lancet 2016, 387, 1723–1731. [Google Scholar] [CrossRef]
  7. Berkhemer, O.A.; Fransen, P.S.S.; Beumer, D.; van den Berg, L.A.; Lingsma, H.F.; Yoo, A.J.; Schonewille, W.J.; Vos, J.A.; Nederkoorn, P.J.; Wermer, M.J.H.; et al. A Randomized Trial of Intraarterial Treatment for Acute Ischemic Stroke. N. Engl. J. Med. 2015, 372, 11–20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Goyal, M.; Demchuk, A.M.; Menon, B.K.; Eesa, M.; Rempel, J.L.; Thornton, J.; Roy, D.; Jovin, T.G.; Willinsky, R.A.; Sapkota, B.L.; et al. Randomized Assessment of Rapid Endovascular Treatment of Ischemic Stroke. N. Engl. J. Med. 2015, 372, 1019–1030. [Google Scholar] [CrossRef] [PubMed]
  9. Saver, J.L.; Goyal, M.; Bonafe, A.; Diener, H.C.; Levy, E.I.; Pereira, V.M.; Albers, G.W.; Cognard, C.; Cohen, D.J.; Hacke, W.; et al. Stent-Retriever Thrombectomy after Intravenous t-PA vs. t-PA Alone in Stroke. N. Engl. J. Med. 2015, 372, 2285–2295. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Minnerup, J.; Sutherland, B.A.; Buchan, A.M.; Kleinschnitz, C. Neuroprotection for Stroke: Current Status and Future Perspectives. Int. J. Mol. Sci. 2012, 13, 11753–11772. [Google Scholar] [CrossRef]
  11. Yenari, M.A.; Hemmen, T.M. Therapeutic Hypothermia for Brain Ischemia: Where have we come and where do we go? Stroke 2010, 41, S72–S74. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Boyle, C.J.; Lennon, A.B.; Prendergast, P.J. In Silico Prediction of the Mechanobiological Response of Arterial Tissue: Application to Angioplasty and Stenting. J. Biomech. Eng. 2011, 133, 081001. [Google Scholar] [CrossRef] [PubMed]
  13. Damiano, R.J.; Ma, D.; Xiang, J.; Siddiqui, A.H.; Snyder, K.V.; Meng, H. Finite element modeling of endovascular coiling and flow diversion enables hemodynamic prediction of complex treatment strategies for intracranial aneurysm. J. Biomech. 2015, 48, 3332–3340. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Carniel, E.L.; Toniolo, I.; Fontanella, C.G. Computational Biomechanics: In-Silico Tools for the Investigation of Surgical Procedures and Devices. Bioengineering 2020, 7, 48. [Google Scholar] [CrossRef]
  15. Fiehler, J.; Thomalla, G.; Bernhardt, M.; Kniep, H.; Berlis, A.; Dorn, F.; Eckert, B.; Kemmling, A.; Langner, S.; Remonda, L.; et al. ERASER. Stroke 2019, 50, 1275–1278. [Google Scholar] [CrossRef]
  16. Grosser, M.; Gellißen, S.; Borchert, P.; Sedlacik, J.; Nawabi, J.; Fiehler, J.; Forkert, N.D. Improved multi-parametric prediction of tissue outcome in acute ischemic stroke patients using spatial features. PLoS ONE 2020, 15, e0228113. [Google Scholar] [CrossRef]
  17. Kemmling, A.; Flottmann, F.; Forkert, N.D.; Minnerup, J.; Heindel, W.; Thomalla, G.; Eckert, B.; Knauth, M.; Psychogios, M.; Langner, S.; et al. Multivariate Dynamic Prediction of Ischemic Infarction and Tissue Salvage as a Function of Time and Degree of Recanalization. J. Cereb. Blood Flow Metab. 2015, 35, 1397–1405. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Kidwell, C.S.; Wintermark, M.; De Silva, D.A.; Schaewe, T.J.; Jahan, R.; Starkman, S.; Jovin, T.; Hom, J.; Jumaa, M.; Schreier, J.; et al. Multiparametric MRI and CT Models of Infarct Core and Favorable Penumbral Imaging Patterns in Acute Ischemic Stroke. Stroke 2013, 44, 73–79. [Google Scholar] [CrossRef] [Green Version]
  19. Livne, M.; Boldsen, J.K.; Mikkelsen, I.K.; Fiebach, J.B.; Sobesky, J.; Mouridsen, K. Boosted Tree Model Reforms Multimodal Magnetic Resonance Imaging Infarct Prediction in Acute Stroke. Stroke 2018, 49, 912–918. [Google Scholar] [CrossRef]
  20. Modrau, B.; Winder, A.; Hjort, N.; Johansen, M.N.; Andersen, G.; Fiehler, J.; Vorum, H.; Forkert, N.D. Machine Learning-Based Prediction of Brain Tissue Infarction in Patients With Acute Ischemic Stroke Treated With Theophylline as an Add-On to Thrombolytic Therapy: A Randomized Clinical Trial Subgroup Analysis. Front. Neurol. 2021, 12, 613029. [Google Scholar] [CrossRef]
  21. Winder, A.J.; Siemonsen, S.; Flottmann, F.; Thomalla, G.; Fiehler, J.; Forkert, N.D. Technical considerations of multi-parametric tissue outcome prediction methods in acute ischemic stroke patients. Sci. Rep. 2019, 9, 1–12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Wu, O.; Koroshetz, W.J.; Østergaard, L.; Buonanno, F.S.; Copen, W.A.; Gonzalez, R.G.; Rordorf, G.; Rosen, B.R.; Schwamm, L.H.; Weisskoff, R.M.; et al. Predicting Tissue Outcome in Acute Human Cerebral Ischemia Using Combined Diffusion- and Perfusion-Weighted MR Imaging. Stroke 2001, 32, 933–942. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Yu, Y.; Xie, Y.; Thamm, T.; Gong, E.; Ouyang, J.; Huang, C.; Christensen, S.; Marks, M.P.; Lansberg, M.G.; Albers, G.W.; et al. Use of Deep Learning to Predict Final Ischemic Stroke Lesions From Initial Magnetic Resonance Imaging. JAMA Netw. Open 2020, 3, e200772. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Vercio, L.L.; Amador, K.; Bannister, J.J.; Crites, S.; Gutierrez, A.; Macdonald, M.E.; Moore, J.; Mouches, P.; Rajashekar, D.; Schimert, S.; et al. Supervised machine learning tools: A tutorial for clinicians. J. Neural Eng. 2020, 17, 062001. [Google Scholar] [CrossRef]
  25. Forkert, N.D.; Fiehler, J. Effect of sample size on multi-parametric prediction of tissue outcome in acute ischemic stroke using a random forest classifier. In Proceedings of the SPIE Medical Imaging 2015: Biomedical Applications in Molecular, Structural, and Functional Imaging, Orlando, FL, USA, 21–26 February 2015; Volume 9417. [Google Scholar] [CrossRef]
  26. Zang, N.; Lin, Z.; Huang, K.; Pan, Y.; Wu, Y.; Wu, Y.; Wang, S.; Wang, D.; Ji, Z.; Pan, S. Biomarkers of Unfavorable Outcome in Acute Ischemic Stroke Patients with Successful Recanalization by Endovascular Thrombectomy. Cerebrovasc. Dis. 2020, 49, 583–592. [Google Scholar] [CrossRef] [PubMed]
  27. Demeestere, J.; Wouters, A.; Christensen, S.; Lemmens, R.; Lansberg, M.G. Review of Perfusion Imaging in Acute Ischemic Stroke. Stroke 2020, 51, 1017–1024. [Google Scholar] [CrossRef]
  28. Boulanger, J.M.; Lindsay, M.P.; Gubitz, G.; Smith, E.E.; Stotts, G.; Foley, N.; Bhogal, S.; Boyle, K.; Braun, L.; Goddard, T.; et al. Canadian Stroke Best Practice Recommendations for Acute Stroke Management: Prehospital, Emergency Department, and Acute Inpatient Stroke Care, 6th Edition, Update 2018. Int. J. Stroke 2018, 13, 949–984. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Goyal, M.; Ospel, J.M.; Menon, B.; Almekhlafi, M.; Jayaraman, M.; Fiehler, J.; Psychogios, M.; Chapot, R.; van der Lugt, A.; Liu, J.; et al. Challenging the Ischemic Core Concept in Acute Ischemic Stroke Imaging. Stroke 2020, 51, 3147–3155. [Google Scholar] [CrossRef] [PubMed]
  30. Winzeck, S.; Hakim, A.; McKinley, R.; Pinto, J.A.A.D.S.R.; Alves, V.; Silva, C.; Pisov, M.; Krivov, E.; Belyaev, M.; Monteiro, M.; et al. ISLES 2016 and 2017-Benchmarking Ischemic Stroke Lesion Outcome Prediction Based on Multispectral MRI. Front. Neurol. 2018, 9, 679. [Google Scholar] [CrossRef]
  31. Martins, N.; Aires, A.; Mendez, B.; Boned, S.; Rubiera, M.; Tomasello, A.; Coscojuela, P.; Hernandez, D.; Muchada, M.; Rodriguez-Luna, D.; et al. Ghost Infarct Core and Admission Computed Tomography Perfusion: Redefining the Role of Neuroimaging in Acute Ischemic Stroke. Interv. Neurol. 2018, 7, 513–521. [Google Scholar] [CrossRef]
  32. Shin, W.; Horowitz, S.; Ragin, A.; Chen, Y.; Walker, M.; Carroll, T.J. Quantitative cerebral perfusion using dynamic susceptibility contrast MRI: Evaluation of reproducibility and age- and gender-dependence with fully automatic image postprocessing algorithm. Magn. Reson. Med. 2007, 58, 1232–1241. [Google Scholar] [CrossRef] [PubMed]
  33. Klimas, A.; Drzazga, Z.; Kluczewska, E.; Hartel, M. Regional ADC measurements during normal brain aging in the clinical range of b values: A DWI study. Clin. Imaging 2013, 37, 637–644. [Google Scholar] [CrossRef] [PubMed]
  34. Arakawa, S.; Wright, P.M.; Koga, M.; Phan, T.G.; Reutens, D.C.; Lim, I.; Gunawan, M.R.; Ma, H.; Perera, N.; Ly, J.; et al. Ischemic Thresholds for Gray and White Matter. Stroke 2006, 37, 1211–1216. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Payabvash, S.; Souza, L.C.S.; Wang, Y.; Schaefer, P.W.; Furie, K.L.; Halpern, E.F.; Gonzalez, R.G.; Lev, M.H. Regional Ischemic Vulnerability of the Brain to Hypoperfusion. Stroke 2011, 42, 1255–1260. [Google Scholar] [CrossRef] [Green Version]
  36. Forkert, N.D.; Cheng, B.; Kemmling, A.; Thomalla, G.; Fiehler, J. ANTONIA Perfusion and Stroke. Methods Inf. Med. 2014, 53, 469–481. [Google Scholar] [CrossRef]
  37. Wu, O.; Østergaard, L.; Weisskoff, R.M.; Benner, T.; Rosen, B.R.; Sorensen, A.G. Tracer arrival timing-insensitive technique for estimating flow in MR perfusion-weighted imaging using singular value decomposition with a block-circulant deconvolution matrix. Magn. Reson. Med. 2003, 50, 164–174. [Google Scholar] [CrossRef]
  38. Hussain, I.; Park, S.J. Quantitative Evaluation of Task-Induced Neurological Outcome after Stroke. Brain Sci. 2021, 11, 900. [Google Scholar] [CrossRef]
  39. National Institute of Neurological Disorders and Stroke rt-PA Stroke Study Group. Tissue Plasminogen Activator for Acute Ischemic Stroke. N. Engl. J. Med. 1995, 333, 1581–1588. [Google Scholar] [CrossRef]
  40. Maus, V.; Styczen, H.; Liman, J.; Maier, I.; Brehm, A.; Tsogkas, I.; Psychogios, M.N. Intracranial mechanical thrombectomy of large vessel occlusions in the posterior circulation using SAVE. BMC Neurol. 2019, 19, 197. [Google Scholar] [CrossRef] [Green Version]
  41. Weber, R.; Minnerup, J.; Nordmeyer, H.; Eyding, J.; Krogias, C.; Hadisurya, J.; Berger, K.; The REVASK Investigators. Thrombectomy in posterior circulation stroke: Differences in procedures and outcome compared to anterior circulation stroke in the prospective multicentre REVASK registry. Eur. J. Neurol. 2019, 26, 299–305. [Google Scholar] [CrossRef]
  42. Bhogal, P.; Hellstern, V.; AlMatter, M.; Ganslandt, O.; Bäzner, H.; Pérez, M.A.; Henkes, H. Mechanical thrombectomy in children and adolescents: Report of five cases and literature review. Stroke Vasc. Neurol. 2018, 3, 245–252. [Google Scholar] [CrossRef] [PubMed]
  43. Sun, L.R.; Felling, R.J.; Pearl, M.S. Endovascular mechanical thrombectomy for acute stroke in young children. J. NeuroInterv. Surg. 2019, 11, 554–558. [Google Scholar] [CrossRef]
  44. Yedavalli, V.S.; Tong, E.; Martin, D.; Yeom, K.W.; Forkert, N.D. Artificial intelligence in stroke imaging: Current and future perspectives. Clin. Imaging 2021, 69, 246–254. [Google Scholar] [CrossRef] [PubMed]
  45. Yu, Y.; Li, M.; Liu, L.; Li, Y.; Wang, J. Clinical big data and deep learning: Applications, challenges, and future outlooks. Big Data Min. Anal. 2019, 2, 288–305. [Google Scholar] [CrossRef] [Green Version]
  46. Hicks, K.A.; Mahaffey, K.W.; Mehran, R.; Nissen, S.E.; Wiviott, S.D.; Dunn, B.; Solomon, S.D.; Marler, J.R.; Teerlink, J.R.; Farb, A.; et al. 2017 Cardiovascular and Stroke Endpoint Definitions for Clinical Trials. Circulation 2018, 137, 961–972. [Google Scholar] [CrossRef]
  47. Forkert, N.D.; Verleger, T.; Cheng, B.; Thomalla, G.; Hilgetag, C.C.; Fiehler, J. Multiclass Support Vector Machine-Based Lesion Mapping Predicts Functional Outcome in Ischemic Stroke Patients. PLoS ONE 2015, 10, e0129569. [Google Scholar] [CrossRef]
  48. Kim, S.M.; Kwon, S.U.; Kim, J.S.; Kang, D.W. Early infarct growth predicts long-term clinical outcome in ischemic stroke. J. Neurol. Sci. 2014, 347, 205–209. [Google Scholar] [CrossRef]
  49. Ma, L.; Gao, P.Y.; Hu, Q.M.; Lin, Y.; Jing, L.N.; Xue, J.; Chen, Z.J.; Wang, Y.J.; Liu, M.L.; Cai, Y.F. Effect of baseline magnetic resonance imaging (MRI) apparent diffusion coefficient lesion volume on functional outcome in ischemic stroke. Neurol. Res. 2011, 33, 494–502. [Google Scholar] [CrossRef]
  50. Menezes, N.M.; Ay, H.; Zhu, M.W.; Lopez, C.J.; Singhal, A.B.; Karonen, J.O.; Aronen, H.J.; Liu, Y.; Nuutinen, J.; Koroshetz, W.J.; et al. The Real Estate Factor. Stroke 2007, 38, 194–197. [Google Scholar] [CrossRef] [Green Version]
  51. Rajashekar, D.; Wilms, M.; Hecker, K.G.; Hill, M.D.; Dukelow, S.; Fiehler, J.; Forkert, N.D. The Impact of Covariates in Voxel-Wise Lesion-Symptom Mapping. Front. Neurol. 2020, 11, 854. [Google Scholar] [CrossRef] [PubMed]
  52. Chalos, V.; van der Ende, N.A.M.; Lingsma, H.F.; Mulder, M.J.H.L.; Venema, E.; Dijkland, S.A.; Berkhemer, O.A.; Yoo, A.J.; Broderick, J.P.; Palesch, Y.Y.; et al. National Institutes of Health Stroke Scale An Alternative Primary Outcome Measure for Trials of Acute Treatment for Ischemic Stroke. Stroke 2020, 51, 282–290. [Google Scholar] [CrossRef] [PubMed]
  53. Ganesh, A.; Menon, B.K.; Assis, Z.A.; Demchuk, A.M.; Al-Ajlan, F.S.; Al-Mekhlafi, M.A.; Rempel, J.L.; Shuaib, A.; Baxter, B.W.; Devlin, T.; et al. Discrepancy between post-treatment infarct volume and 90-day outcome in the ESCAPE randomized controlled trial. Int. J. Stroke 2021, 16, 593–601. [Google Scholar] [CrossRef] [PubMed]
  54. Hill, M.D.; Goyal, M.; Menon, B.K.; Nogueira, R.G.; McTaggart, R.A.; Demchuk, A.M.; Poppe, A.Y.; Buck, B.H.; Field, T.S.; Dowlatshahi, D.; et al. Efficacy and safety of nerinetide for the treatment of acute ischaemic stroke (ESCAPE-NA1): A multicentre, double-blind, randomised controlled trial. Lancet 2020, 395, 878–887. [Google Scholar] [CrossRef]
Biomedicines 09 01357 g001 550
Figure 1. Graphical abstract illustrating the proposed in silico method.
Figure 1. Graphical abstract illustrating the proposed in silico method.
Biomedicines 09 01357 g001
Biomedicines 09 01357 g002 550
Figure 2. Selected follow-up lesion segmentations generated for each combination of outcome prediction model and patient, alongside the patients’ corresponding ground-truth follow-up lesion segmentation. Recanalization outcomes are abbreviated IAR (intra-arterial therapy, recanalizing), IVR (intravenous tPA, recanalizing), and NR (non-recanalizing).
Figure 2. Selected follow-up lesion segmentations generated for each combination of outcome prediction model and patient, alongside the patients’ corresponding ground-truth follow-up lesion segmentation. Recanalization outcomes are abbreviated IAR (intra-arterial therapy, recanalizing), IVR (intravenous tPA, recanalizing), and NR (non-recanalizing).
Biomedicines 09 01357 g002
Biomedicines 09 01357 g003 550
Figure 3. Boxplots representing the distribution of lesion volumes calculated for (a) the expert-segmented follow-up lesions, separated by the treatment outcome of the patient, (b) the predicted follow-up lesion, separated by predicted treatment outcomes. Recanalization outcomes are abbreviated IAR (intra-arterial therapy, recanalizing), IVR (intravenous tPA, recanalizing), and NR (non-recanalizing). Significance markers (***) represent p < 0.001 for post-hoc tests following the one-way ANOVA performed in (a) and the one-way repeated-measures ANOVA performed in (b). Circular markers represent outliers with values more than 1.5 times the interquartile range.
Figure 3. Boxplots representing the distribution of lesion volumes calculated for (a) the expert-segmented follow-up lesions, separated by the treatment outcome of the patient, (b) the predicted follow-up lesion, separated by predicted treatment outcomes. Recanalization outcomes are abbreviated IAR (intra-arterial therapy, recanalizing), IVR (intravenous tPA, recanalizing), and NR (non-recanalizing). Significance markers (***) represent p < 0.001 for post-hoc tests following the one-way ANOVA performed in (a) and the one-way repeated-measures ANOVA performed in (b). Circular markers represent outliers with values more than 1.5 times the interquartile range.
Biomedicines 09 01357 g003
Biomedicines 09 01357 g004 550
Figure 4. Receiver operating characteristic (ROC) curves for the standard predictions of each model. Recanalization outcomes for the models are (a) IAR (intra-arterial therapy, recanalizing), (b) IVR (intravenous tPA, recanalizing), and (c) NR (non-recanalizing). AUC values are (a) 0.979, (b) 0.976, and (c) 0.957.
Figure 4. Receiver operating characteristic (ROC) curves for the standard predictions of each model. Recanalization outcomes for the models are (a) IAR (intra-arterial therapy, recanalizing), (b) IVR (intravenous tPA, recanalizing), and (c) NR (non-recanalizing). AUC values are (a) 0.979, (b) 0.976, and (c) 0.957.
Biomedicines 09 01357 g004
Table
Table 1. Descriptive statistics comparing each model’s standard predictions and the corresponding follow-up lesion segmentations. Recanalization outcomes are abbreviated IAR (intra-arterial therapy, recanalizing), IVR (intravenous tPA, recanalizing), and NR (non-recanalizing).
Table 1. Descriptive statistics comparing each model’s standard predictions and the corresponding follow-up lesion segmentations. Recanalization outcomes are abbreviated IAR (intra-arterial therapy, recanalizing), IVR (intravenous tPA, recanalizing), and NR (non-recanalizing).
Dice CoefficientProportional Volume Error
MeanStandard Dev.MedianInterquartile Range
IAR Model0.4420.2310.1821.80
IVR Model0.4480.255−0.1392.10
NR Model0.4730.2450.0541.59
Table
Table 2. Additional statistics derived from the confusion matrices comparing each model’s standard predictions and the corresponding follow-up lesion segmentations. Recanalization outcomes are abbreviated IAR (intra-arterial therapy, recanalizing), IVR (intravenous tPA, recanalizing), and NR (non-recanalizing).
Table 2. Additional statistics derived from the confusion matrices comparing each model’s standard predictions and the corresponding follow-up lesion segmentations. Recanalization outcomes are abbreviated IAR (intra-arterial therapy, recanalizing), IVR (intravenous tPA, recanalizing), and NR (non-recanalizing).
AccuracySensitivitySpecificityPrecision
MeanStd. Dev.MeanStd. Dev.MeanStd. Dev.MeanStd. Dev.
IAR Model0.9780.0260.5990.2760.9900.0130.5070.308
IVR Model0.9860.0090.5670.2630.9930.0070.4740.298
NR Model0.9710.0220.6190.2950.9840.0170.4890.279
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Winder, A.; Wilms, M.; Fiehler, J.; Forkert, N.D. Treatment Efficacy Analysis in Acute Ischemic Stroke Patients Using In Silico Modeling Based on Machine Learning: A Proof-of-Principle. Biomedicines 2021, 9, 1357. https://doi.org/10.3390/biomedicines9101357

AMA Style

Winder A, Wilms M, Fiehler J, Forkert ND. Treatment Efficacy Analysis in Acute Ischemic Stroke Patients Using In Silico Modeling Based on Machine Learning: A Proof-of-Principle. Biomedicines. 2021; 9(10):1357. https://doi.org/10.3390/biomedicines9101357

Chicago/Turabian Style

Winder, Anthony, Matthias Wilms, Jens Fiehler, and Nils D. Forkert. 2021. "Treatment Efficacy Analysis in Acute Ischemic Stroke Patients Using In Silico Modeling Based on Machine Learning: A Proof-of-Principle" Biomedicines 9, no. 10: 1357. https://doi.org/10.3390/biomedicines9101357

APA Style

Winder, A., Wilms, M., Fiehler, J., & Forkert, N. D. (2021). Treatment Efficacy Analysis in Acute Ischemic Stroke Patients Using In Silico Modeling Based on Machine Learning: A Proof-of-Principle. Biomedicines, 9(10), 1357. https://doi.org/10.3390/biomedicines9101357

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop