Skip to main content
Download PDF

Machine learning-based segmentation of ischemic penumbra by using diffusion tensor metrics in a rat model

Journal of Biomedical Science volume 27, Article number: 80 (2020) Cite this article

Abstract

Background

Recent trials have shown promise in intra-arterial thrombectomy after the first 6–24 h of stroke onset. Quick and precise identification of the salvageable tissue is essential for successful stroke management. In this study, we examined the feasibility of machine learning (ML) approaches for differentiating the ischemic penumbra (IP) from the infarct core (IC) by using diffusion tensor imaging (DTI)-derived metrics.

Methods

Fourteen male rats subjected to permanent middle cerebral artery occlusion (pMCAO) were included in this study. Using a 7 T magnetic resonance imaging, DTI metrics such as fractional anisotropy, pure anisotropy, diffusion magnitude, mean diffusivity (MD), axial diffusivity, and radial diffusivity were derived. The MD and relative cerebral blood flow maps were coregistered to define the IP and IC at 0.5 h after pMCAO. A 2-level classifier was proposed based on DTI-derived metrics to classify stroke hemispheres into the IP, IC, and normal tissue (NT). The classification performance was evaluated using leave-one-out cross validation.

Results

The IC and non-IC can be accurately segmented by the proposed 2-level classifier with an area under the receiver operating characteristic curve (AUC) between 0.99 and 1.00, and with accuracies between 96.3 and 96.7%. For the training dataset, the non-IC can be further classified into the IP and NT with an AUC between 0.96 and 0.98, and with accuracies between 95.0 and 95.9%. For the testing dataset, the classification accuracy for IC and non-IC was 96.0 ± 2.3% whereas for IP and NT, it was 80.1 ± 8.0%. Overall, we achieved the accuracy of 88.1 ± 6.7% for classifying three tissue subtypes (IP, IC, and NT) in the stroke hemisphere and the estimated lesion volumes were not significantly different from those of the ground truth (p = .56, .94, and .78, respectively).

Conclusions

Our method achieved comparable results to the conventional approach using perfusion–diffusion mismatch. We suggest that a single DTI sequence along with ML algorithms is capable of dichotomizing ischemic tissue into the IC and IP.

Background

Stroke is one of the major causes of long-term disability and death, and nearly 80% of stroke cases are ischemic [1]. Treatment options for acute ischemic stroke (AIS) are rapid recanalization of the occluded large vessels by using intravenous (IV) thrombolysis with tissue plasminogen activator (tPA) and intra-arterial (IA) thrombectomy to mechanically disrupt or remove the thrombus. In either treatment, identifying a substantial and salvageable ischemic penumbra (IP) is essential for a patient to be eligible for therapy [2,3,4]. In the DAWN and DEFUSE 3 trials, which included acute stroke patients within 6–24 h of onset, obtaining perfusion imaging computed tomography (CT) perfusion or magnetic resonance imaging (MRI) perfusion-weighted imaging (PWI) or an MRI with a diffusion-weighted imaging (DWI) sequence was recommended to help determine whether the patient is a candidate for mechanical thrombectomy [2, 5]. In the acute setting, an infarct core (IC) can be identified through DWI and combined with the hypoperfusion area depicted by PWI, which allows for the specific definition of the salvageable IP and IC by using the concept of perfusion–diffusion mismatch (PDM). However, quick and accurate delineation of IP is demanded by clinicians for AIS management.

Diffusion tensor imaging (DTI) has been used in clinical applications for measuring cerebral microstructural changes induced by neurological diseases [6]. In AIS, DTI-derived metrics, such as fractional anisotropy (FA), mean diffusivity (MD), pure anisotropic diffusion (q), and diffusion magnitude (L), have demonstrated the feasibility of DTI in assessing the damage of ischemic brain tissue [7], determining the onset time of AIS in hours in an animal model [8] or in humans [9, 10], estimating the salvageable tissue [8], and microstructurally discriminating benign oligemia from the “true” penumbral tissue [11]. Based on the previous studies, DTI may provide comprehensive characterization of the pathophysiological process of cerebral ischemia.

Machine learning (ML) has become a useful aid for physicians in the diagnosis of, treatment of, and prediction of complications and patient outcomes for numerous diseases. While extracting meaningful and discriminative imaging features that exhibits the characteristics of lesion part [12, 13], ML-based algorithms can establish predictive models for various clinical applications [12, 14]. In the current study, we attempted to develop a 2-level ML classifier based on DTI-derived metrics for characterizing ischemic tissue subtypes. We aimed to determine whether the DTI sequence along with ML algorithms could classify the stroke hemisphere into the IP, IC, and unaffected normal tissue (NT) during the AIS stage.

Materials and methods

Animals

Fourteen male Sprague-Dawley rats (weight, 270–350 g; Taipei Medical University Animal Center, Taiwan) were used in this study. The rats were housed in a humidity- and temperature-controlled environment and placed under a 12:12-h light–dark cycle, with free access to sterile food and water. All the rats underwent permanent middle cerebral artery occlusion (pMCAO) through an intraluminal suture method based on the modified Zea Longa approach [15]. All animal experiments were approved and performed in accordance with guidelines and regulations of the Institutional Animal Care and Use Committee of Taipei Medical University (IACUC approval No: LAC-2015-0033).

Image acquisition

Images were acquired using a 7 T MRI scanner (PharmaScan 70/16; Bruker Biospin, Ettlingen, Germany). During image acquisition, the rats were placed under anesthesia by using 1.5–2% isoflurane and the rectal temperature was maintained at approximately 37 °C using a warm water bath with continuous circulation through a water-bath temperature controller set outside the magnet. DTI was performed with 6 noncollinear diffusion encoding gradients with a b factor of 1200 s/mm2 and 1 b = 0 s/mm2. Multishot echoplanar imaging (repetition time [TR] = 3000 ms, echo time [TE] = 37 ms, number of excitations = 6) with the navigator-echo correction technique was used as the signal readout module. To obtain the largest IP area and avoid fast diminishing, PWI was performed once at 0.5 h after pMCAO by using a dynamic susceptibility contrast technique. A series of gradient-echo echoplanar coronal images with a TR/TE value of 1000/20 ms and 300 repetitions were acquired. A bolus of the susceptibility contrast agent gadolinium-diethylenetriamine penta-acetic acid (0.25 mmol/kg; Magnevist, Bayer Schering Pharma, Berlin, Germany) was injected manually through the rat tail vein approximately 30 s after the start of image acquisition. All the images acquired from DTI and PWI were with a field-of-view of 20 mm × 20 mm and a matrix of 64 × 64, which were subsequently zero-filled to 128 × 128 with a resolution of 0.16 mm × 0.16 mm for further analyses.

Data analysis

Calculation of the relative cerebral blood flow and DTI metrics

The relative cerebral blood flow (rCBF) and DTI metrics were calculated using in-house algorithms in MATLAB (MathWorks, Natick, MA, USA). First, the relative cerebral blood volume (rCBV) and relative mean transit time (rMTT) were determined using the integral and normalized first moment of gamma variate fitting, respectively. Next, rCBF was derived as the quotient of rCBV divided by rMTT by using the central volume principle [16]. For DTI metrics, the eigenvalues of each image voxel were computed and then applied to derive the MD, FA, q, L, axial diffusivity (AD), and radial diffusivity (RD). To avoid noise or artifacts, all the DTI-derived and rCBF maps were computed through a convolution with a Gaussian kernel using the weighted mean intensity value. All the smoothed maps were then normalized linearly to the range [0, 1] for inter-rat comparisons.

Delineation of the IP, IC, and NT

With supervised learning procedures in the classifiers, the labels of the IP, IC, and NT should be determined in advance. According to previous research, abnormal MD (i.e., IC) is defined using a reduction of 30% of the contralateral hemisphere with the exclusion of the ventricles (Fig. 1a) [17]. Perfusion deficit is defined with a lower CBF threshold of 46% reduction of the contralateral hemisphere (Fig. 1b) [18]. The rCBF map was coregistered to the MD maps to delineate the perfusion–diffusion mismatch (i.e., the IP). Regions without CBF deficit within the ipsilateral brain were defined as the NT. Contiguity correction was performed to remove “misclassified” pixels (Fig. 1c) [19, 20]. Finally, as shown in Fig. 1d, the regions of the IP, IC and NT were depicted and the corresponding voxels were labeled.

Fig. 1
figure 1

Definitions of the Ischemic Penumbra (IP), Infarct core (IC), and Normal Tissue (NT) in a Rat Subjected to Permanent Middle Cerebral Artery Occlusion (pMCAO). IC was defined as the blue area in the mean diffusivity (MD) map (a) and perfusion deficit at 0.5 h after pMCAO is shown in (b). Perfusion–diffusion mismatch is illustrated in (c) and (d), where the red region indicates the IC and the green region indicates the IP. The NT was defined as the region in the ipsilateral hemisphere except for the IP and IC [white in (d)]. d Indicates the “label” for the IC, IP, and NT

Full size image

Feature extraction

Each voxel was characterized by 110 features and was segmented into the IP, IC, or NT according to the definitions in the previous section. Three types of features (the relative DTI-derived metrics; Mahalanobis distance; normalized histogram described by kurtosis, skewness, and bin counts) were extracted from regions of interest in the voxel-located slices and adjacent slices. These features represented the spatial patterns around the voxel as well as the spatial relationship between the voxel and its neighboring voxels.

Relative DTI-derived metrics

Six relative DTI-derived metrics (rMD, rAD, rRD, rFA, rL, and rq) were obtained for each pixel. Once the regions of the IP, IC and NT were depicted, the relative DTI metrics were derived on a pixel-by-pixel basis in relation to the contralateral homologous tissue as follows: rX = (Xipsilateral − Xcontralateral) / Xcontralateral, where X indicates the particular DTI index. In consideration of the lesion spatial contiguity, metrics of two vertical adjacent voxels in the adjacent slices were also used. Thus, a total of 18 DTI-derived features were computed for each voxel.

Mahalanobis distance

Although previous studies have adopted the spatial Euclidean distance from the lesion location as features [21, 22], various ischemic distributions and multifocal ischemia make stroke lesion segmentation difficult [21]. In this study, once the IC and non-IC were classified by the first-level classifier (Fig. 2a), the Mahalanobis distances [19] between each non-IC feature vector and the distribution of all IC feature vectors were calculated as discriminative features for the second-level classifier.

Fig. 2
figure 2

Strategy of 2-level Classification and Validation. a The proposed 2-level strategy for voxel-wise classification to classify every voxel in the hemisphere into a tissue subtype (i.e., the IC, IP, or NT). b The validation methods include 5-fold cross validation in training phase and leave-one-out cross validation for the final prediction

Full size image

Normalized histogram (kurtosis, skewness, and bin counts)

The normalized histogram-specified 11 bins of the relative DTI metrics in an area of the 7 × 7 window in the coronal planes were calculated [21]. In addition, we reconstructed the axial projection maps from the coronal slices to calculate the normalized histogram in an area of the 3 × 3 window in the axial planes. The kurtosis, skewness, and bin counts of each histogram distribution were calculated as features. A total of 90 features were obtained for each voxel.

Two-level classification

We propose a 2-level classification model composed of two binary classifiers to hierarchically classify the stroke hemisphere into three tissue subtypes, as displayed in Fig. 2a, and the classification performance was evaluated through the leave-one-out cross validation (LOOCV) method (i.e., 13 rats were used for training and a remaining rat for testing each time, Fig. 2b). A total of 141,806 samples (voxels) were obtained from 14 rats within 71 slices. During the training phase, five-fold cross validation (CV) was used to prevent the possible bias of overfitting. We trained the first-level classifiers by using the 18 DTI features to classify hemisphere into the IC and non-IC. We then used all 110 features for the second-level classifiers. For comparison, a single-level classification model using either support vector machine (SVM), k-nearest neighbors (KNN), or decision tree algorithms was constructed to classify three tissue subtypes simultaneously based on 18 DTI features. The classification methods were implemented using the Statistics and Machine Learning Toolbox in the MATLAB.

Validation and statistical analysis

The performance of the classification was evaluated by calculating the accuracy, sensitivity, specificity, and areas under the receiver operating characteristic curve (AUC). The slice-to-slice correspondence analysis between the estimated lesion volume and the volume defined by PDM was also conducted. To evaluate the applicability and potential of our proposed approach for the whole ischemic lesion volume, the Mann–Whitney U-test was also used for comparison between the average estimated lesion volumes and those defined by PDM. Statistical tests were performed using SPSS® (Version 19.0; SPSS Inc., Chicago, IL, USA). All the group data are reported as mean ± SD, and the significance level was defined at a p value < .05.

Results

Figure 3a illustrates the maps of DTI metrics at 0.5 h post-pMCAO. The q, L, MD, AD, and RD maps demonstrate initial hypointensity changes in the ischemic areas while the FA map exhibits symmetrical signal intensity. Figure 3b shows six relative DTI-derived metrics (rMD, rAD, rRD, rFA, rL, and rq) obtained for each voxel.

Fig. 3
figure 3

Maps of Diffusion Tensor Imaging (DTI) Metrics Measured at 0.5 Hours After pMCAO. a Significant hypointensities on the ischemic lesion can be observed from MD, L, q, AD and RD maps but not from FA map. b The relative DTI metrics are shown. The intensity of the map represents the quantitative decreases or increases of the DTI metrics compared with the corresponding contralateral homologous tissue

Full size image

For the 2-level classifiers, the SVM algorithm revealed the best overall discrimination performance. Its first-level classifier exhibited a high performance for discriminating the IC from the non-IC in the training dataset (AUC: 0.99; accuracy: 96.3%; sensitivity: 0.95; specificity: 0.97). The discrimination performance remained high for the second-level classifier (IP vs. NT), with AUCs, accuracies, and specificities up to 0.96, 95.0% and 0.97, respectively. However, the sensitivities were relatively low (0.85 ~ 0.86). The detailed classification performance using the LOOCV is listed in Table 1. On the other hand, the single-level classifiers exhibited poor sensitivities for IP detection (SVM: 29.0%; KNN: 46.8%; decision tree: 4.9%), as shown in Table 2.

Table 1 Performances of the 2-level classifiers for the training dataset
Full size table
Table 2 Performances of the single-level classifiers for the training dataset
Full size table

The trained SVM classifiers were then applied to the testing dataset, and the results are presented in Table 3. The classification accuracy of the IC and non-IC segmentations was 96.0 ± 2.3%. For further subclassification within the non-IC, the classification accuracy of IP and NT segmentation was 80.1 ± 8.0%. Overall, the classification accuracy of the segmentation of the 3 tissue subtypes was 88.1 ± 6.7%, and the median was 90.0% within the cerebral hemisphere of stroke (ranging from 69.5 to 96.9%). We observed a favorable classification outcome for the testing rats.

Table 3 Performances of the 2-level SVM Classifiers for the Testing Dataset
Full size table

The slice-to-slice correspondence analysis between the classifier-estimated volume and the corresponding volume defined by PDM for the 3 tissue subtypes in all 71 slices is illustrated in Fig. 4. The correlation coefficients between the volume estimated from one slice and the volume defined by PDM were 0.613, 0.999, and 0.932 (with all p < .001) for the IP, IC, and NT, respectively. The estimated IC and NT volumes exhibited excellent correlations with the conventional PDM measures, whereas the IP exhibited a moderate correlation.

Fig. 4
figure 4

Results of Slice-to-Slice Correspondence Analysis Between the Classifier-Estimated Volume and the Perfusion-Diffusion Mismatch -Defined Volume for the IP (a), IC (b), and NT (c)

Full size image

Figure 5 depicts a comparison of the average of the classifier-estimated lesion volumes with the average of the volumes defined by PDM across 14 rats for the 3 tissue subtypes. The volumes of the 3 tissue subtypes were calculated and summed up across slices. Overall, the classifier-estimated lesion volume was comparable with those defined by PDM-defined lesion volume in all tissue subtypes. In the Mann–Whitney U-test, the estimated IP (17.8 ± 10.4 mm3), IC (104.8 ± 72.6 mm3), and NT (110.7 ± 55.9 mm3) volumes were not significantly different from those defined by PDM for the IP (20.9 ± 10.5 mm3, p = .56), IC (104.3 ± 73.2 mm3, p = .94), and NT (108.7 ± 56.8 mm3, p = .78) volumes.

Fig. 5
figure 5

Average of the Estimated Volume and Perfusion-Diffusion Mismatch Defined Volume for the IP, IC, and NT for the 14 rats. No significant differences were observed between the classifier-estimated volume and the perfusion-diffusion mismatch defined volume in the IP (P = .56), IC (P = .94) and NT (P = .78)

Full size image

Figure 6 illustrates the comparison of the classifier-estimated IC and IP with the corresponding PDM-defined IC and IP for a rat. In the suture-occlusion model, the IP is relatively small (even sparse) in areas at the margin of a large IC. Nevertheless, the proposed classification model can successfully segment the ischemic regions into the IP and IC through visual inspection.

Fig. 6
figure 6

Demonstration of a Predicted Mismatch for 2 Rats. The conventional perfusion–diffusion mismatch and estimated mismatch are illustrated, where the red region indicates the IC and the green region indicates the IP. The NT is displayed in grayscale

Full size image

Discussion

In the present study, we developed a 2-level classification model with an overall accuracy of 88.1 ± 6.7% for discriminating the stroke hemisphere into the IC, IP, and NT regions on a voxel-wise basis in a pMCAO model. According to the analysis results, we suggest that a single DTI sequence combined with ML algorithms is capable of dichotomizing ischemic tissue into the IC and IP, which are comparable to the conventional PDM and DTI studies. This classifier system could assist clinical AIS management for acute triage and patient selections without the use of intravenous gadolinium-based contrast agent. Thus, the costs for contrast agent, the time for IV access, and the possibility of nephrogenic systemic fibrosis in renal-insufficient patients can be avoid [23]. Most previous imaging studies using ML tools have focused on diagnostic accuracy and prediction of prognosis after stroke [14, 24]. The current study moved a step further on the development of a quick and well-validated classifier to depict the salvageable tissue during an AIS phase.

Currently, the therapeutic time window for patients with AIS who are eligible for intravenous thrombolysis is within 4.5 h of onset [25,26,27]. In the DEFUSE 3 and DAWN trials, the time window for intra-arterial thrombectomy was extended up to 16–24 h [2, 5, 28]. These penumbra-based stroke trials have shown that patients with favorable clinical response generally have a relatively small IC (volume < 70 cm3) and large penumbra (PWI volume/DWI volume ratio > 1.8). Therefore, a quick and direct measurement of the IP at the acute phase through single imaging acquisition could be considerably helpful in stroke trials. Recent studies employing DTI metrics to characterize the IP after pMCAO have shown a persistent IP up to 6.5 h. DTI metrics exhibit promise in differentiating the IP from the oligemia and IC and even in determining the time of stroke onset [8]. DTI has been applied to measure the cerebral microstructural integrity of cell membranes after ischemia by characterizing their tensor magnitude, orientation, and anisotropy such that different tissue injuries can be stratified. The proposed 2-level binary classification is based on the same principle that DTI metrics could provide information about the varying degree of microstructural damages following ischemic changes. The combination of DTI metrics with ML algorithms could thus enable the precise stratification of tissue subtypes, which would lead to quick estimation of the salvageable tissue in the acute phase of a stroke.

In our study, the classifier exhibited a strong correlation (r = 0.999, p < .001; Fig. 4b) between the estimated IC volume and abnormal MD-defined IC volume. However, only a moderate correlation could be achieved between the estimated IP volume and the PDM-defined volume (r = 0.613, p < .001; Fig. 4a). Several factors may have caused the moderate correlation. One factor is that the intraluminal suture MCAO model may produce variable lesion distributions and the CBF reduction can be highly heterogeneous [29,30,31], which presumably leads to CBF-based viability threshold deviation in the range of ± 11% [18]. In other words, the deviation can result in overlapping measurements at the borders between the IP and NT areas, which leads to inappropriate preassigned class for ML. Another factor is the minor changes of DTI metrics in the IP (decrease approximately 10%) during the hyperacute stage. Such changes have been reported in humans [32, 33] and in rodent models [8, 11]. The small difference in DTI metrics between the IP and NT (an approximately 40% difference in DTI metrics between the IC and NT) that are challenging to separate in the original or even feature space may have contributed to the relatively low accuracy of discrimination between the IP and the NT (80.1 ± 8.0%; Table 3). Nonetheless, our findings showed that the classifiers can assist in estimating the IP and IC volumes, which are comparable with the lesion volumes calculated from conventional PDM (all P > .556, Fig. 5). The proposed classifiers may also be used for the 14–27% of patients who have unknown onset time (such as in wake-up strokes) [34] but still have substantial salvageable tissue volumes. Nonetheless, IP gradually diminish and eventually progress into IC tissue without recanalization. Previous study, however, has shown 100% normobaric hyperoxia (NBO) was helpful to “freeze” the IP and prolong the treatment time window [35]. Combined the proposed classifiers may help to select a candidate for reperfusion and evaluate the treatment efficacy repeatedly.

Limited studies have applied ML algorithms for the discrimination of the IC and salvageable tissue. A study using 2 convolutional neural networks demonstrated a 94% lesion detection rate by training the segmented IC in DWI from 741 patients [36]. Despite the difference in ML algorithms, the results of the aforementioned study were similar (approximately 95–96% for the IC; Table 1) to the results obtained with the support vector machine (SVM)-based model in this study. The high IC detection rates obtained are probably attributable to the inherent high sensitivity of DWI in depicting acute infarct [37,38,39,40]. Our classifier can achieve not only a high IC detection rate but can also separate the non-IC into the IP and NT through DTI alone. An animal study by Huang et al. [41] used an SVM-based model to predict acute ischemic tissue fate through the CBF and apparent diffusion coefficient mismatch. The AUCs in the 30-min, 60-min, and pMCAO groups were 88 ± 2 .9%, 94 ± 0.8%, and 97 ± 0.9%, respectively. The marginally lower performance of our 2-level classification model was because we did not use the perfusion parametric maps as features in our training process. Multimodal MR features from different imaging sequences, such as contrast enhancement, T2 fluid attenuation inversion recovery, and T1/T2-weighted imaging, have been shown to enhance the performance of ML-based approaches, either in lesion segmentation or the prediction of tissue outcome [21, 22, 41, 42]. However, adding multiple MR sequences may result in long scanning times, which may not be feasible for patients with AIS simply because “time is brain.”

The performance of classifiers highly depends on the features extracted from the images. Although the end-to-end deep learning model, which covers both feature extraction and classification procedures, has exhibited suitable capability for classifying natural images, the small samples size of medical images limits its reliability in clinical applications. Radiomics is a recently developed computational pipeline for extracting quantitative features from medical images [43, 44]. These features include first-order statistical features, shape-based features, and the most used texture features in oncology studies [43,44,45]. In the current experimental model, we used first-order statistical features (i.e., skewness and kurtosis), shape-based features (i.e., distance measures for similarity [21, 22]), and texture features for the spatial information of the lesions. Due to the relatively small number of rats and imbalanced distribution of the IC, IP, and NT, we adopted the LOOCV in the training phase to obtain a larger number of training data than that obtained using the K-fold partition scheme [46]. This approach also provides an unbiased evaluation of the performance [47]. The small number of rats made the usage of a convolutional neural network a challenge because the training procedure was based on the entire image. However, the pixel-wise analysis along with SVM provided sufficient samples (N = 141,806) for model training as well as a reasonable ratio of samples to data dimensions (d = 110). To avoid the possible problem of overfitting, we adopted a 5-fold CV for evaluating the proposed methods. In addition, estimation of volumetric IP and localization would be implemented by reconstruction of 3D from 2D results to precisely know how the spatially location relative to normal tissue is, which may help physicians’ assessment before reperfusion surgery.

This study has some limitations. First, correctly labeled data is essential in supervised learning such as SVM. Our research used a PDM overlay for DTI segmentation. The conventional PDM region may exceed the true penumbral area and is usually much larger than the final size of the infarct [48,49,50]. The perfusion deficit overestimates the region at risk, including the penumbral tissue and benign oligemia [51]. Other modalities, such as positron emission tomography [52, 53] or MRI using oxygen challenge with T2* signal change [54], could provide detailed insight into the complex pathophysiological changes of the brain after ischemia and may be used to additionally define the penumbral tissue. In our study, we adopted 46% reduction of CBF to define penumbra by referencing Meng’s result [18]. With this criterion, we also observed the consistency with the final infarct regions in the T2-weighted images. Second, possible redundant features may be eliminated through advanced feature selection or dimensionality reduction methods to improve the accuracy and reduce the computational time. Also, because features are essential in classification, further efforts should be made to explore stroke-related features. Finally, the region-wise histopathological correlation, which is difficult to perform because of the dynamic changes within the penumbra, may be implemented to provide further proof of the value of our proposed classifiers.

Conclusions

Our results suggested that single DTI combined with ML algorithms could provide a noninvasive, quick, and reliable method of assessing the salvageable tissue, thus accelerating the management of patients with AIS. To the best of our knowledge, this is the first ML-based study to demonstrate the potential of using a single DTI sequence for substituting the conventional approach of PDM, offering a practical workflow for clinical decision-making and stroke trials.

Availability of data and materials

The data generated or analyzed during the current study is not publicly available due to restrictions in the ethical permit, but may be available from the corresponding author on request.

References

  1. Mozaffarian D, Benjamin EJ, Go AS, Arnett DK, Blaha MJ, Cushman M, et al. Executive summary: heart disease and stroke statistics—2015 update: a report from the American Heart Association. Circulation. 2015;131(4):434–41.

    Google Scholar 

  2. Albers GW, Marks MP, Kemp S, Christensen S, Tsai JP, Ortega-Gutierrez S, et al. Thrombectomy for stroke at 6 to 16 hours with selection by perfusion imaging. N Engl J Med. 2018;378(8):708–18.

    PubMed  PubMed Central  Google Scholar 

  3. Thomalla G, Simonsen CZ, Boutitie F, Andersen G, Berthezene Y, Cheng B, et al. MRI-guided thrombolysis for stroke with unknown time of onset. N Engl J Med. 2018;379(7):611–22.

    PubMed  Google Scholar 

  4. Albers GW, Goyal M, Jahan R, Bonafe A, Diener HC, Levy EI, et al. Ischemic core and hypoperfusion volumes predict infarct size in SWIFT PRIME. Ann Neurol. 2016;79(1):76–89.

    PubMed  Google Scholar 

  5. Hacke W. A new DAWN for imaging-based selection in the treatment of acute stroke. N Engl J Med. 2018;378(1):81–3.

    PubMed  Google Scholar 

  6. Sundgren P, Dong Q, Gomez-Hassan D, Mukherji S, Maly P, Welsh R. Diffusion tensor imaging of the brain: review of clinical applications. Neuroradiology. 2004;46(5):339–50.

    CAS  PubMed  Google Scholar 

  7. Urbanski M, De Schotten MT, Rodrigo S, Oppenheim C, Touzé E, Méder J-F, et al. DTI-MR tractography of white matter damage in stroke patients with neglect. Exp Brain Res. 2011;208(4):491–505.

    CAS  PubMed  Google Scholar 

  8. Kuo D-P, Lu C-F, Liou M, Chen Y-C, Chung H-W, Chen C-Y. Differentiation of the infarct core from ischemic penumbra within the first 4.5 hours, using diffusion tensor imaging-derived metrics: a rat model. Korean J Radiol. 2017;18(2):269–78.

    PubMed  PubMed Central  Google Scholar 

  9. Sakai K, Yamada K, Nagakane Y, Mori S, Nakagawa M, Nishimura T. Diffusion tensor imaging may help the determination of time at onset in cerebral ischaemia. J Neurol Neurosurg Psychiatry. 2009;80(9):986–90.

    CAS  PubMed  Google Scholar 

  10. Puig J, Blasco G, Daunis-I-Estadella J, Thomalla G, Castellanos M, Soria G, et al. Increased corticospinal tract fractional anisotropy can discriminate stroke onset within the first 4.5 hours. Stroke. 2013;44(4):1162–5.

    PubMed  Google Scholar 

  11. Chiu F-Y, Kuo D-P, Chen Y-C, Kao Y-C, Chung H-W, Chen C-Y. Diffusion tensor-derived properties of benign oligemia, true “at risk” penumbra, and infarct core during the first three hours of stroke onset: a rat model. Korean J Radiol. 2018;19(6):1161–71.

    PubMed  PubMed Central  Google Scholar 

  12. Erickson BJ, Korfiatis P, Akkus Z, Kline TL. Machine learning for medical imaging. Radiographics. 2017;37(2):505–15.

    PubMed  Google Scholar 

  13. Nixon M, Aguado A. Feature extraction and image processing for computer vision: Academic press; 2019.

    Google Scholar 

  14. Lee E-J, Kim Y-H, Kim N, Kang D-W. Deep into the brain: artificial intelligence in stroke imaging. J Stroke. 2017;19(3):277.

    PubMed  PubMed Central  Google Scholar 

  15. Longa EZ, Weinstein PR, Carlson S, Cummins R. Reversible middle cerebral artery occlusion without craniectomy in rats. Stroke. 1989;20(1):84–91.

    CAS  PubMed  Google Scholar 

  16. Calamante F, Gadian D, Connelly A. Quantification of perfusion using bolus tracking magnetic resonance imaging in stroke: assumptions, limitations, and potential implications for clinical use. Stroke. 2002;33(4):1146–51.

    CAS  PubMed  Google Scholar 

  17. Shen Q, Fisher M, Sotak CH, Duong TQ. Effects of reperfusion on ADC and CBF pixel-by-pixel dynamics in stroke: characterizing tissue fates using quantitative diffusion and perfusion imaging. J Cereb Blood Flow Metab. 2004;24(3):280–90.

    PubMed  Google Scholar 

  18. Meng X, Fisher M, Shen Q, Sotak CH, Duong TQ. Characterizing the diffusion/perfusion mismatch in experimental focal cerebral ischemia. Ann Neurol. 2004;55(2):207–12.

    PubMed  PubMed Central  Google Scholar 

  19. Shen Q, Ren H, Fisher M, Bouley J, Duong TQ. Dynamic tracking of acute ischemic tissue fates using improved unsupervised ISODATA analysis of high-resolution quantitative perfusion and diffusion data. J Cereb Blood Flow Metab. 2004;24(8):887–97.

    PubMed  Google Scholar 

  20. Huang S, Shen Q, Duong TQ. Artificial neural network prediction of ischemic tissue fate in acute stroke imaging. J Cereb Blood Flow Metab. 2010;30(9):1661–70.

    PubMed  PubMed Central  Google Scholar 

  21. Maier O, Wilms M, von der Gablentz J, Krämer UM, Münte TF, Handels H. Extra tree forests for sub-acute ischemic stroke lesion segmentation in MR sequences. J Neurosci Methods. 2015;240:89–100.

    PubMed  Google Scholar 

  22. Maier O, Schröder C, Forkert ND, Martinetz T, Handels H. Classifiers for ischemic stroke lesion segmentation: a comparison study. PLoS One. 2015;10(12):e0145118.

    PubMed  PubMed Central  Google Scholar 

  23. Cowper SE, Robin HS, Steinberg SM, Su LD, Gupta S, LeBoit PE. Scleromyxoedema-like cutaneous diseases in renal-dialysis patients. Lancet. 2000;356(9234):1000–1.

    CAS  PubMed  Google Scholar 

  24. Kamal H, Lopez V, Sheth SA. Machine learning in acute ischemic stroke neuroimaging. Front Neurol. 2018;9:945.

    PubMed  PubMed Central  Google Scholar 

  25. Wahlgren N, Ahmed N, Dávalos A, Hacke W, Millán M, Muir K, et al. Thrombolysis with alteplase 3–4· 5 h after acute ischaemic stroke (SITS-ISTR): an observational study. Lancet. 2008;372(9646):1303–9.

    CAS  PubMed  Google Scholar 

  26. Ahmed N, Wahlgren N, Grond M, Hennerici M, Lees KR, Mikulik R, et al. Implementation and outcome of thrombolysis with alteplase 3–4· 5 h after an acute stroke: an updated analysis from SITS-ISTR. Lancet Neurol. 2010;9(9):866–74.

    CAS  PubMed  Google Scholar 

  27. Hacke W, Kaste M, Bluhmki E, Brozman M, Dávalos A, Guidetti D, et al. Thrombolysis with alteplase 3 to 4.5 hours after acute ischemic stroke. N Engl J Med. 2008;359(13):1317–29.

    CAS  PubMed  Google Scholar 

  28. Powers WJ, Rabinstein AA, Ackerson T, Adeoye OM, Bambakidis NC, Becker K, et al. 2018 guidelines for the early management of patients with acute ischemic stroke: a guideline for healthcare professionals from the American Heart Association/American Stroke Association. Stroke. 2018;49(3):e46–99.

    PubMed  Google Scholar 

  29. Li F, Han S, Tatlisumak T, Carano RA, Irie K, Sotak CH, et al. A new method to improve in-bore middle cerebral artery occlusion in rats. Stroke. 1998;29:1715–20.

    CAS  PubMed  Google Scholar 

  30. Roussel S, Van Bruggen N, King M, Gadian D. Identification of collaterally perfused areas following focal cerebral ischemia in the rat by comparison of gradient echo and diffusion-weighted MRI. J Cereb Blood Flow Metab. 1995;15(4):578–86.

    CAS  PubMed  Google Scholar 

  31. Lythgoe MF, Williams SR, Busza AL, Wiebe L, McEwan AJ, Gadian DG, et al. The relationship between magnetic resonance diffusion imaging and autoradiographic markers of cerebral blood flow and hypoxia in an animal stroke model. Magn Reson Med. 1999;41(4):706–14.

    CAS  PubMed  Google Scholar 

  32. Schaefer PW, Ozsunar Y, He J, Hamberg LM, Hunter GJ, Sorensen AG, et al. Assessing tissue viability with MR diffusion and perfusion imaging. Am J Neuroradiol. 2003;24(3):436–43.

    PubMed  PubMed Central  Google Scholar 

  33. Røhl L, Østergaard L, Simonsen CZ, Vestergaard-Poulsen P, Andersen G, Sakoh M, et al. Viability thresholds of ischemic penumbra of hyperacute stroke defined by perfusion-weighted MRI and apparent diffusion coefficient. Stroke. 2001;32(5):1140–6.

    PubMed  Google Scholar 

  34. Kang D-W, Kwon JY, Kwon SU, Kim JS. Wake-up or unclear-onset strokes: are they waking up to the world of thrombolysis therapy? Int J Stroke. 2012;7(4):311–20.

    PubMed  Google Scholar 

  35. Baron J-C. Protecting the ischaemic penumbra as an adjunct to thrombectomy for acute stroke. Nat Rev Neurol. 2018;14(6):325.

    CAS  PubMed  Google Scholar 

  36. Chen L, Bentley P, Rueckert D. Fully automatic acute ischemic lesion segmentation in DWI using convolutional neural networks. Neuroimage. 2017;15:633–43.

    PubMed  PubMed Central  Google Scholar 

  37. Smajlović D, Sinanović O. Sensitivity of the neuroimaging techniques in ischemic stroke. Med Arh. 2004;58(5):282–4.

    PubMed  Google Scholar 

  38. Muir KW, Buchan A, von Kummer R, Rother J, Baron J-C. Imaging of acute stroke. Lancet Neurol. 2006;5(9):755–68.

    PubMed  Google Scholar 

  39. Moseley M, Kucharczyk J, Mintorovitch J, Cohen Y, Kurhanewicz J, Derugin N, et al. Diffusion-weighted MR imaging of acute stroke: correlation with T2-weighted and magnetic susceptibility-enhanced MR imaging in cats. Am J Neuroradiol. 1990;11(3):423–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  40. Schellinger P, Bryan R, Caplan L, Detre J, Edelman R, Jaigobin C, et al. Evidence-based guideline: the role of diffusion and perfusion MRI for the diagnosis of acute ischemic stroke: report of the Therapeutics AND Technology Assessment Subcommittee of the American Academy of Neurology. Neurology. 2010;75(2):177–85.

    CAS  PubMed  PubMed Central  Google Scholar 

  41. Huang S, Shen Q, Duong TQ. Quantitative prediction of acute ischemic tissue fate using support vector machine. Brain Res. 2011;1405:77–84.

    CAS  PubMed  PubMed Central  Google Scholar 

  42. Bouts MJ, Tiebosch IA, Van Der Toorn A, Viergever MA, Wu O, Dijkhuizen RM. Early identification of potentially salvageable tissue with MRI-based predictive algorithms after experimental ischemic stroke. J Cereb Blood Flow Metab. 2013;33(7):1075–82.

    PubMed  PubMed Central  Google Scholar 

  43. Lambin P, Rios-Velazquez E, Leijenaar R, Carvalho S, Van Stiphout RG, Granton P, et al. Radiomics: extracting more information from medical images using advanced feature analysis. Eur J Cancer. 2012;48(4):441–6.

    PubMed  PubMed Central  Google Scholar 

  44. Aerts HJ, Velazquez ER, Leijenaar RT, Parmar C, Grossmann P, Carvalho S, et al. Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach. Nat Commun. 2014;5:4006.

    CAS  PubMed  Google Scholar 

  45. Gillies RJ, Kinahan PE, Hricak H. Radiomics: images are more than pictures, they are data. Radiology. 2015;278(2):563–77.

    PubMed  Google Scholar 

  46. Bishop CM. Pattern recognition and machine learning: springer; 2006.

    Google Scholar 

  47. Cawley GC, Talbot NL. Efficient approximate leave-one-out cross-validation for kernel logistic regression. Mach Learn. 2008;71(2–3):243–64.

    Google Scholar 

  48. Schlaug G, Benfield A, Baird A, Siewert B, Lövblad K, Parker R, et al. The ischemic penumbra: operationally defined by diffusion and perfusion MRI. Neurology. 1999;53(7):1528.

    CAS  PubMed  Google Scholar 

  49. Sorensen AG, Copen WA, Østergaard L, Buonanno FS, Gonzalez RG, Rordorf G, et al. Hyperacute stroke: simultaneous measurement of relative cerebral blood volume, relative cerebral blood flow, and mean tissue transit time. Radiology. 1999;210(2):519–27.

    CAS  PubMed  Google Scholar 

  50. Neumann-Haefelin T, Wittsack H Jr, Wenserski F, Siebler M, Seitz RJ, Mödder U, et al. Diffusion-and perfusion-weighted MRI: the DWI/PWI mismatch region in acute stroke. Stroke. 1999;30(8):1591–7.

    CAS  PubMed  Google Scholar 

  51. Kidwell CS, Alger JR, Saver JL. Beyond mismatch: evolving paradigms in imaging the ischemic penumbra with multimodal magnetic resonance imaging. Stroke. 2003;34(11):2729–35.

    PubMed  Google Scholar 

  52. Sobesky J. Refining the mismatch concept in acute stroke: lessons learned from PET and MRI. J Cereb Blood Flow Metab. 2012;32(7):1416–25.

    PubMed  PubMed Central  Google Scholar 

  53. Sobesky J, Weber OZ, Lehnhardt F-G, Hesselmann V, Neveling M, Jacobs A, et al. Does the mismatch match the penumbra? Magnetic resonance imaging and positron emission tomography in early ischemic stroke. Stroke. 2005;36(5):980–5.

    PubMed  Google Scholar 

  54. Santosh C, Brennan D, McCabe C, Macrae IM, Holmes WM, Graham DI, et al. Potential use of oxygen as a metabolic biosensor in combination with T2*-weighted MRI to define the ischemic penumbra. J Cereb Blood Flow Metab. 2008;28(10):1742–53.

    CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

The study was supported by the Ministry of Science and Technology, Taiwan (MOST105-2221-E-038-007-MY3) and by Taipei Medical University Hospital, Taipei, Taiwan (106TMUH-SP-04 and 107TMUH-SP-02).

Author information

Authors and Affiliations

  1. Department of Medical Imaging, Taipei Medical University Hospital, No.250, Wu-Hsing St, Taipei, 11031, Taiwan

    Duen-Pang Kuo, Yung-Chieh Chen & Cheng-Yu Chen

  2. Department of Radiology, Taoyuan Armed Forces General Hospital, Taoyuan, Taiwan

    Duen-Pang Kuo

  3. Institute for Medical Engineering and Science, Massachusetts Institute of Technology, Cambridge, MA, USA

    Po-Chih Kuo

  4. Department of Biomedical Imaging and Radiological Sciences, National Yang-Ming University, No.155, Sec.2, Linong St, Taipei, 11221, Taiwan

    Yu-Chieh Jill Kao & Cheng-Yu Chen

  5. Department of Radiology, School of Medicine, College of Medicine, Taipei Medical University, No.250, Wu-Hsing St, Taipei, 11031, Taiwan

    Cheng-Yu Chen

  6. TMU Center for Big Data and Artificial Intelligence in Medical Imaging, Taipei Medical University Hospital, Taipei, Taiwan

    Ching-Yen Lee

  7. TMU Research Center for Artificial Intelligence in Medicine, Taipei Medical University Hospital, Taipei, Taiwan

    Ching-Yen Lee

  8. Graduate Institute of Biomedical Electrics and Bioinformatics, National Taiwan University, Taipei, Taiwan

    Hsiao-Wen Chung

  9. Radiogenomic Research Center, Taipei Medical University Hospital, No.250, Wu-Hsing St, Taipei, 11031, Taiwan

    Cheng-Yu Chen

  10. Center for Artificial Intelligence in Medicine, Taipei Medical University, No.250, Wu-Hsing St, Taipei, 11031, Taiwan

    Cheng-Yu Chen

  11. Department of Radiology, National Defense Medical Center, No.250, Wu-Hsing St, Taipei, 11031, Taiwan

    Cheng-Yu Chen

Authors
  1. Duen-Pang Kuo
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Po-Chih Kuo
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Yung-Chieh Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Yu-Chieh Jill Kao
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Ching-Yen Lee
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Hsiao-Wen Chung
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Cheng-Yu Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Duen-Pang Kuo participated in the design of the study, carried out the data collection and data interpretation, and drafted the manuscript. Po-Chih Kuo participated in the design of the study, carried out the data analysis and statistical analysis, and drafted the manuscript. Yung-Chieh Chen participated in data interpretation and revised the manuscript. Yu-Chieh Jill Kao participated in building the pMCAO model. Ching-Yen Lee assisted the data interpretation and revised the manuscript. Hsiao-Wen Chung participated in the coordination of the study and assisted in data interpretation. Cheng-Yu Chen participated in the design and coordination of the study, assisted in data interpretation, revised the manuscript, and supervised the study. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Cheng-Yu Chen.

Ethics declarations

Ethics approval

All animal experiments were approved and performed in accordance with guidelines and regulations of the Institutional Animal Care and Use Committee of Taipei Medical University (IACUC approval No: LAC-2015-0033).

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kuo, DP., Kuo, PC., Chen, YC. et al. Machine learning-based segmentation of ischemic penumbra by using diffusion tensor metrics in a rat model. J Biomed Sci 27, 80 (2020). https://doi.org/10.1186/s12929-020-00672-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12929-020-00672-9

Keywords

Journal of Biomedical Science

ISSN: 1423-0127