Skip to main content

The hemorrhagic transformation index score: a prediction tool in middle cerebral artery ischemic stroke



We aimed to develop a tool, the hemorrhagic transformation (HT) index (HTI), to predict any HT within 14 days after middle cerebral artery (MCA) stroke onset regardless of the intravenous recombinant tissue plasminogen activator (IV rtPA) use. That is especially important in the light of missing evidence-based data concerning the timing of anticoagulant resumption after stroke in patients with atrial fibrillation (AF).


We retrospectively analyzed 783 consecutive MCA stroke patients. Clinical and brain imaging data at admission were recorded. A follow-up period was 2 weeks after admission. The patients were divided into derivation (DC) and validation (VC) cohorts by generating Bernoulli variates with probability parameter 0.7. Univariate/multivariate logistic regression, and factor analysis were used to extract independent predictors. Validation was performed with internal consistency reliability and receiver operating characteristic (ROC) analysis. Bootstrapping was used to reduce bias.


The HTI was composed of 4 items: Alberta Stroke Program Early CT score (ASPECTS), National Institutes of Health Stroke Scale (NIHSS), hyperdense MCA (HMCA) sign, and AF on electrocardiogram (ECG) at admission. According to the predicted probability (PP) range, scores were allocated to ASPECTS as follows: 10–7 = 0; 6–5 = 1; 4–3 = 2; 2–0 = 3; to NIHSS: 0–11 = 0; 12–17 = 1; 18–23 = 2; >23 = 3; to HMCA sign: yes = 1; to AF on ECG: yes = 1. The HTI score varied from 0 to 8. For each score, adjusted PP of any HT with 95% confidence intervals (CI) was as follows: 0 = 0.027 (0.011–0.042); 1 = 0.07 (0.043–0.098); 2 = 0.169 (0.125–0.213); 3 = 0.346 (0.275–0.417); 4 = 0.571 (0.474–0.668); 5 = 0.768 (0.676–0.861); 6 = 0.893 (0.829–0.957); 7 = 0.956 (0.92–0.992); 8 = 0.983 (0.965–1.0). The optimal cutpoint score to differentiate between HT-positive and negative groups was 2 (95% normal-based CI, 1–3) for the DC and VC alike. ROC area/sensitivity/specificity with 95% normal-based CI for the DC and VC were 0.85 (0.82–0.89)/0.82 (0.73–0.9)/0.89 (0.8–0.97) and 0.83 (0.78–0.88)/0.8 (0.66–0.94)/0.87 (0.73–1.0) respectively. McDonald’s categorical omega with 95% bias-corrected and accelerated CI for the DC and VC was 0.81 (0.77–0.84) and 0.82 (0.76–0.86) respectively.


The HTI is a simple yet reliable tool to predict any HT within 2 weeks after MCA stroke onset regardless of the IV rtPA use.

Peer Review reports


Hemorrhagic transformation (HT), either asymptomatic (AHT) or symptomatic (SHT), is considered to be a notorious complication of acute ischemic stroke (AIS), associated with limited treatment options and long-term adverse outcomes [1]. It seems reasonable that efforts should be directed towards preventing HT before it occurs. Fortunately, it is more predictable than other types of intracranial hemorrhage.

In AIS patients, the incidence of HT induced by intravenous recombinant tissue plasminogen activator (IV rtPA) is reported to be 4.5–39.6% for AHT and 5.2–7.3% for SHT. In contrast, the rate of spontaneous AHT and SHT ranges from 13% to 43% and from 0.6% to 20% respectively [2, 3]. Although the proportion of AIS patients treated with IV rtPA is relatively small (4.7–21.4%) [4], the majority of authors have focused on searching HT predictors coupled with IV rtPA over the past decade. As a result, a variety of predictive clinical scores have emerged [5,6,7,8,9].

On the other hand, there is a lack of tools for making an accurate HT prediction in AIS patients who are not eligible for IV rtPA. That is especially important in the light of missing evidence-based data concerning the timing of anticoagulant resumption after AIS in patients with atrial fibrillation (AF). Recommendations on the initiation of anticoagulation are currently based on consensus opinion, in what is known as the “1–3–6–12 day rule” [10]. Therefore, the two-week timeframe following the AIS onset is the most critical for developing HT. In this instance, an accurate prediction of HT could make a difference in decision making to reinstitute anticoagulation. The middle cerebral artery (MCA) is by far the largest cerebral artery and is the vessel most commonly affected by cerebrovascular accident.

Given the background, we aimed to develop a simple and yet reliable instrument called the hemorrhagic transformation index (HTI) to predict any HT within 14 days after AIS onset in the MCA territory regardless of the use of IV rtPA.



Using prospectively collected clinical and radiological databases, we retrospectively identified 783 consecutive patients with AIS in the MCA territory who were admitted to the stroke unit of the Interregional Clinical Diagnostic Center, Kazan, Russia, within 12 h after onset between January 2013 and May 2016. The exclusion criteria were: involvement of other vascular territories; AIS following any surgery or endovascular procedure within 1 month; brain ischemic lesions due to an intracranial tumor, infection, cerebral venous thrombosis, subarachnoid hemorrhage, and arteriovenous malformation/fistula. In total, 1361 AIS patients were admitted over the specified period. The sample was drawn from the local Caucasian population.

The eligible patients received diagnostic tests and treatment in accordance with current national stroke guidelines. The permissible hospital length of stay was at least 14 days, which was determined by the state mandatory medical insurance standard for AIS patients.

Clinical baseline variables, including age, sex, risk factors, pre-admission medication, stroke subtype according to the Trial of ORG 10172 in Acute Stroke Treatment classification, NIHSS score, vital signs, blood tests, electrocardiogram (ECG), echocardiogram, and chest X-ray findings at admission were extracted from the medical charts. The NIHSS score was routinely and systematically assessed by neurologists. The time of AIS onset was documented as described by the patient or witness; if unknown, it was considered to be the last time the patient was seen well. In-hospital antithrombotic medication was logged for 14 consecutive days; the log was withdrawn earlier if HT occurred.

Imaging protocol

Brain non-contrast computed tomography (CT) was performed using a multidetector CT scanner (Aquilion 64; Toshiba Medical Systems, Otawara, Japan). All CT scans were obtained with 0.5 mm slice thickness; the technical parameters were as follows: 120 kVp, 300 mA, rotation time 0.75 s, matrix size 512 × 512, helical scan mode, total scan time 9.7 s, reconstruction interval 5 mm. Window levels and widths were optimized for gray/white matter distinction. The Alberta Stroke Program Early CT score (ASPECTS), hyperdense MCA (HMCA) sign, and leukoaraiosis were routinely and systematically recorded at admission by radiologists. The HMCA sign was assessed by measurements of absolute attenuation of the affected and normal vessels. Absolute density of the affected MCA of >43 Hounsfield units and the MCA ratio of >1.2 on a non-contrast CT scan were regarded as the positive HMCA sign [11]. Diffuse hypodense areas involving the periventricular and/or centrum semiovale white matter were considered as leukoaraiosis. A follow-up CT scan was routinely done on hospitalization day 7 and 14 or at any time if required by a treating neurologist. All patients had at least one follow-up CT scan.

Outcome measures

The outcome was retrospectively revised based on prospectively collected data. Any HT on a follow-up CT scan within 14 days after AIS onset was taken into account. A hemorrhage was considered symptomatic if it was not seen on a previous CT scan and there had subsequently been either a suspicion of hemorrhage or any decline in neurologic status [12]. According to the ECASS I trial [13], HT was further classified into hemorrhagic infarction type 1 (HI-1), or type 2 (HI-2), or parenchymal hematoma type 1 (PH-1), or type 2 (PH-2).

Statistical analysis

Multinomial logistic regression with relative risk ratio (RRR) estimation was used to highlight the AHT and SHT association with poor outcomes (death, malignant edema, and dependency defined as the modified Rankin scale of >2 at discharge); the baseline category was the HT-negative group.

The intraclass correlation coefficient (ICC) was computed to assess inter-rater agreement for the ASPECTS and NIHSS. In order to calculate the ICC for the NIHSS, the same set of video files with NIHSS examination from six patients with different stroke severity was demonstrated to seven neurologists who regularly admitted patients to our stroke unit. The ICC for the ASPECTS was obtained in the similar manner: the identical pool of brain non-contrast CT scans from 33 patients with different ischemic burden was presented to four radiologists, who regularly evaluated brain CT scans at admission. Each doctor was evaluated separately and independently; the NIHSS and ASPECTS reference manuals were available on request. He or she had a chance to make any corrections during the evaluation process, but was not allowed to do so after his or her assessment had been completed.

A preliminary data analysis showed that 14 variables in 39.21% of observations were missing (Table 1).

Table 1 Variables with Missing Values

The data were not missing completely at random (Little’s test: χ2(1706), 2336.69; p < 0.001). However, the omissions did occur accidentally because some tests were not available at the time of patient’s admission or the results were lost. Moreover, the missing variables correlated with other collected data. Therefore, it was reasonable to assume that the data were missing at random and multiple imputation (MI) was an appropriate technique to manage the absent values (Table 2).

Table 2 MI Specification

After obtaining the imputed data, the observations were divided into derivation (DC) and validation (VC) cohorts by generating Bernoulli variates with probability parameter 0.7.

Descriptive statistics included median with interquartile range (IQR) and percentage for continuous (the distribution was not normal) and categorical data respectively. The NIHSS and ASPECTS were treated as continuous variables because of multiple categories. Lists of univariate and multivariate HT predictors were obtained by fitting a binary logistic regression (BLR) model. Variables with univariate p-value ≤0.25 were further included in multivariate analysis, whereas only items with p-value <0.05 were kept in the multivariate BLR equation. Once the list of HT predictors was obtained by fitting a multivariate BLR model, we dropped the MI dataset because the included variables had no missing values in the source data.

In order to proceed with exploratory (EFA) and confirmatory (CFA) factor analysis, the Bartlett’s test of sphericity, Kaiser-Meyer-Olkin measure, and Doornik-Hansen test were carried out to check for patterned relationships between the HTI items, data sufficiency, and multivariate normality respectively. EFA by means of principal factor (PF) and principal component factor (PCF) techniques was performed with varimax and promax rotations to assess dimensionality of the HTI items and to extract variables with shared variance; an eigenvalue cut-off was 1.0. CFA with maximum likelihood estimation was applied to select the final model; goodness of fit was assessed with the Satorra-Bentler scaled χ2 test to adjust for data non-normality.

HTI internal consistency reliability (ICR) was evaluated with the ordinal α, Guttman λ2 and λ4 bounds, Raykov’s ρ, McDonald’s ω, and greatest lower bound [14, 15]. The values ≥0.7 were considered to be reliable. As the model included ordinal and dichotomous variables, a polychoric correlation matrix was used for EFA and ICR analysis except for McDonald’s categorical ω. The latter was computed by using the Green and Yang method [16].

The DC and VC were compared by using the Mann-Whitney U and Pearson χ2 tests for continuous and categorical variables respectively. The equality of kernel density estimate (KDE) for predicted probability (PP) of any HT between the multivariate BLR model and HTI score as well as for HTI scores between the DC and VC was evaluated with the two-sample Kolmogorov-Smirnov test.

Receiver operating characteristic (ROC) analysis was conducted to assess prognostic performance. The optimal cutpoint score to distinguish between HT-positive and negative groups was defined with the Youden index. Based on the VC appraisal, the predictive ability of the HTI was compared with several alternative tools by testing the area under the ROC curve (AUC) of each score against the HTI one. For each comparison, the Šidák-adjusted p-value was reported. The AUC equality was evaluated by using the DeLong algorithm [17].

Whenever possible, bootstrapping was performed with 1000 samples and computing either adjusted for ties bias-corrected and accelerated (BCa) or normal-based (NB) confidence intervals (CI) to reduce sampling bias, overfitting, and prediction errors.


Overall, HT occurred in 186 (23.8%) out of 783 cases, whereas SHT was determined in 98 (12.5%). HI-1 was established in 26 (3.3%); HI-2, 120 (15.3%); PH-1, 13 (1.7%); and PH-2, 27 (3.4%) observations. Patients with any HT were more likely to have a poor outcome: death (AHT: RRR, 4.8; 95% CI, 2.2–10.6; p < 0.001; SHT: RRR, 11.7; 95% CI, 6–22.6; p < 0.001), malignant edema (AHT: RRR, 16.7; 95% CI, 5.6–49.2; p < 0.001; SHT: RRR, 52.1; 95% CI, 19.6–138.9; p < 0.001), and dependency (AHT: RRR, 4.4; 95% CI, 2.5–7.7; p < 0.001; SHT: RRR, 30.4; 95% CI, 9.5–97.1; p < 0.001).

There was no difference between AHT and SHT groups with regard to timing (p = 0.08): the former, median day 2 (IQR, 1–4); the latter, median day 2 (IQR, 1–3). The median AIS onset time was 6 h (IQR, 2–11). Just under half of all cases (346; 44.2%) were admitted within the 4.5-h therapeutic window; in total, 67 (8.6%) patients were treated with IV rtPA.

The ICC was 0.95 (95% CI, 0.85–0.99; p < 0.001) and 0.78 (95% CI, 0.59–0.88; p < 0.001) for the NIHSS and ASPECTS respectively, which indicated excellent inter-rater agreement. There was no difference between the DC and VC (Table 3).

Table 3 Baseline Characteristics in the DC and VC

Univariate analysis was summarized in Table 4.

Table 4 Univariate Analysis in the DC Using the MI Dataset

Multivariate analysis

Although univariate p-values for leukoaraiosis and international normalized ratio (INR) were above our acceptable threshold, we included them in the multivariate analysis because some authors had proposed them as risk factors [18, 19]. As a result of fitting a multivariate BLR model, seven variables were kept in the final equation (Table 5).

Table 5 Multivariate Analysis (DC, n = 535)

Swapping AF on ECG for the AF history variable increased the Bayesian and Akaike information criteria by 1.18, which slightly favored the initial model. Overall, the multivariate BLR model was statistically significant (Wald test: χ2(7), 87.76; p < 0.001; −2log-likelihood, 271.93; Cox-Snell pseudo-R2, 0.44; Nagelkerke pseudo-R2, 0.67). It explained variance of 78.4% (variance of latent variable, y* = 11.13; error, ε = 3.29) and fitted the data well (Hosmer-Lemeshow goodness of fit test for 10 groups: χ2(8), 6.87; p = 0.551). Neither significant interactions nor polynomial terms were found.

There was no specification error (Pregibon’s link test: linear predicted value, p < 0.001; linear predicted value squared, p = 0.54). Assumption of linearity between independent variables and log odds was confirmed by the LOWESS graph. Multicollinearity was not an issue: the extrema of the variance inflation factor were 1.01 and 2.29.

Although standardized Pearson and deviance residuals exceeded 2 in a few observations, their leverage and Pregibon’s influential statistics (dbeta) turned out to be very small. Moreover, removing those observations did not significantly change the equation coefficients. Influence of each individual observation on the coefficient estimate (not adjusted for the covariate pattern), dfbeta, was not strong. However, the most sensitive was the INR variable (Fig. 1).

Fig. 1
figure 1

Multivariate BLR Model Diagnostics. a LOWESS graph, bandwidth = 0.8. Each hollow circle is an individual observation. b Influence of each individual observation on the coefficient estimate (dfbeta). c Standardized Pearson residuals. d Deviance residuals. e Leverage. The dashed line corresponds to the 3-fold leverage mean. f Pregibon’s influential statistics (dbeta)

Model sensitivity, specificity, positive and negative predictive values were 76.2%, 95.6%, 84.2%, and 92.9% respectively. The model accurately classified 91% of the observations, whereas the equation without any independent variables classified correctly only 76.5% of the cases.

Based on PP tables and plots, we divided each predictor into categories and allocated them points according to the PP range in order to draw the HTI score. Given the Doornik-Hansen test (χ2(14), 3665.64; p < 0.001), the distribution of the newly derived HTI items was not multivariate normal (Table 6; Figs. 2 and 3).

Table 6 Derivation of the HTI Score
Fig. 2
figure 2

Derivation of HTI Score. a ASPECTS. b NIHSS. c INR. d HR on ECG. e HMCA sign. f AF on ECG. g Sex. Grey areas and capped lines represent 95% CI. Numbers at the top indicate HTI scores; at the bottom, units of variables. The confounders are held constant at their observed values

Fig. 3
figure 3

KDE for HTI Items. Kernel = Epanechnikov for all items. a ASPECTS; bandwidth = 0.1899. b NIHSS; bandwidth = 0.1899. c INR; bandwidth = 0.116. d HR on ECG; bandwidth = 0.1528. e Sex; bandwidth = 0.1278. f HMCA sign; bandwidth = 0.113. g AF on ECG; bandwidth = 0.1182

Factor analysis

The Bartlett’s test of sphericity (χ2(21), 721.74; p < 0.001) and the Kaiser-Meyer-Olkin measure of 0.71 demonstrated that the HTI items did have patterned relationships and were sufficient for EFA. PF EFA established a unidimensional scale, i.e. there was only one factor that explained a cumulative variance of 82.8%. By means of PCF EFA, the factor was discovered to consist of cerebral (ASPECTS, NIHSS, and HMCA sign) and extracerebral (resting heart rate (HR) on ECG, AF on ECG, sex, and INR) components as we called them. However, the resting HR on ECG, INR, and sex variables showed high uniqueness values and low factor loadings; therefore, we had to drop them. The final 4-item HTI was strongly supported by subsequent CFA (Tables 7 and 8; Fig. 4).

Table 7 Factor/Component Loadings and Uniqueness in EFA of the HTI Items
Table 8 CFA. Goodness of Fit Statistics
Fig. 4
figure 4

Scree Plot. EFA by using (a) PF and (b) PCF methods

Once the HTI was definitively established, crude PP of any HT was computed for each score by using BLR. Although the dropped items were no longer a part of the HTI, we put them into the BLR equation for confounding adjustment. Unsurprisingly, the confounders exerted only a minor influence on the overall HTI OR by increasing it up to 12% and had no significant effect on the OR of the separate HTI items (Tables 9, 10 and 11; Fig. 5).

Table 9 BLR Coefficients, Bias, and Bootstrap SE for the HTI (DC, n = 535)
Table 10 OR for the HTI (DC, n = 535)
Table 11 PP of Any HT for Each HTI Score
Fig. 5
figure 5

Calibration Plot. The scatterplots display PP of any HT for individual observations obtained from multivariate analysis (Table 5) and arranged by the HTI score: (a) The DC, (b) The VC. The superimposed connected line graph drawn from the DC demonstrates PP of any HT for each HTI score with 95% CI adjusted for sex, INR, and HR on ECG (Table 11). The confounders are held constant at their observed values

KDE for PP of any HT was equal between the multivariate BLR model and HTI score (D = 0.184; p = 0.371). Thus, the HTI score was considered as a surrogate for the multivariate BLR model (Fig. 6).

Fig. 6
figure 6

KDE for PP of Any HT. Kernel = Epanechnikov, bandwidth = 0.0623. The solid line represents the multivariate BLR model (Table 5); the dashed line, the HTI score (Table 11)

ICR and ROC analysis

Given multiple reliability tests, HTI ICR was considered to be fair enough for the DC and VC alike. There was also no difference in the AUC (χ2(1), 0.01; p = 0.93) and KDE for the HTI scores (D = 0.02; p = 1.0) between both cohorts (Table 12; Fig. 7).

Table 12 HTI ICR and ROC Analysis
Fig. 7
figure 7

KDE for HTI Scores. Kernel = Epanechnikov, bandwidth = 0.5697. The solid line represents the DC; the dashed line, the VC

Taking into account that alternative scores had been developed in AIS patients with slightly different clinical settings, the HTI prognostic performance was considered to be at least non-inferior to the competitors (Tables 13 and 14).

Table 13 Comparison of Different Composite Scores for Predicting HT in AIS Patients
Table 14 AUC Comparison Between the HTI and Alternative Scores (VC, n = 248)


The incidence of HT showed in our study echoes the rate reported in literature. We have also reaffirmed the concept that AHT is not clinically innocuous. The study unequivocally reiterates infarct size, stroke severity, large-artery occlusion, and cardioembolism defined by ASPECTS, NIHSS, HMCA sign, and AF respectively are well-established independent HT predictors [20, 21]. In our HTI score, we use the presence of AF on ECG at admission rather than AF history for the reason discussed in the Results. However, other known predictors – platelet count, cholesterol level, age, hypertension, renal failure, hyperglycemia, and leukoaraiosis – have shown no independent association with any HT in our cohort; similar results were obtained by other authors [22,23,24,25,26]. The mechanism of this association needs to be explored.

Although women tend to be more likely HT-positive in our univariate BLR model, multivariate analysis reveals the opposite. It seems there is still controversy about the sex propensity for developing HT [27].

Accelerated HR at rest is known to be associated with an increased risk of stroke especially recurrent [28]. Since AIS commonly induces change in cardiovascular responses, post-stroke HR at admission could be a potential marker to identify patients at risk for short-term deterioration and long-term poor outcomes [29]. To the best of our knowledge, we have not found any literature, concerning HR correlations with HT. Here, we report that the more the HR is accelerated on ECG (but not the pulse rate), the more likely HT can occur.

Hypercoagulability at AIS onset measured by INR and other tests is known to be associated with an increased thrombotic tendency. As long as the hypercoagulable state persists, both arterial and venous thromboembolic recurrences can be expected. The association of these coagulation abnormalities with HT is not always clear [30]. As we have shown here, the less the INR, the higher the risk of HT.

Having identified seven independent variables in our multivariate analysis, we applied the PP range followed by factor analysis to assigned HTI scores and to regroup HTI variables into a limited set of clusters based on shared variance. We assume probability is more intuitive for interpreting than OR. In contrast, authors of other predictive tools allocated scores based on OR changes only [5,6,7,8,9, 31,32,33]. While BLR analyzes effects of each individual predictor on the dependent variable, factor analysis isolating constructs and concepts treats the model as a whole [34]. Thus, it helped us to avoid overfitting.

Among the compared predictive tools, the SPAN-100 and HeRS were derived from the cohorts, which were very similar in clinical settings to ours (Table 13). As to the former, we have reaffirmed that it is far inferior in prognostic performance to other scores [35, 36]. Regarding the latter, the infarct size was measured in different ways: we chose the ASPECTS estimation on non-contrast CT, whereas the HeRS scored it in milliliters on magnetic resonance diffusion-weighted imaging (DWI MRI). Although both approaches are widely acceptable in hyperacute stroke settings, the ASPECTS is more suitable for non-contrast CT assessment, while the lesion volume can easily be quantified on DWI [37,38,39]. Furthermore, the ASPECTS, as well as the NIHSS, correlates strongly with the infarct volume (ASPECTS: Spearman ρ, −0.88; p < 0.001; NIHSS: Spearman ρ, 0.71; p < 0.001; VC, n = 248), but moderately with each other (Spearman ρ, −0.66; p < 0.001; VC, n = 248). If the stroke volume variable had been added to our multivariate BLR equation, a multicollinearity issue would have occurred. We also suppose that a combination of clinical and imaging features is more reliable than the imaging data alone. To make our HTI score as much easy-to-use as possible, we have purposely refrained from MRI since CT is readily available in the most hospitals. Moreover, the HeRS score is computationally complicated; therefore, it seems less attractive from the practical point of view.

There are some important peculiarities between posterior and anterior circulation stroke. The differences include the value of screening instruments, optimum diagnostic modalities, clinical features, and outcomes [40, 41]. For instance, patients with vertebrobasilar infarction have lower NIHSS score and HT rates, less often AF, higher blood glucose level and rates of false-negative DWI findings, more WBC counts, and a better long-term outcome than those with carotid stroke [42,43,44]. However, all aforementioned tools predict HT regardless of the vascular basin (Table 13). Moreover, scores with imaging modalities, like the HAT and SEDAN, include CT signs of MCA stroke only. Meanwhile, a scoring system, the pc-ASPECTS, has been developed and validated for posterior circulation stroke [45]. Thus, the accuracy of predictive tools could be further improved by distinguishing the infarcted vascular basins; therefore, we have decided to restrict our study to the MCA territory.

There are a few limitations in our study. A relatively small, but sufficient for statistical inferences, sample size and lack of ethnic and racial diversity could be a source of potential bias. Almost all patients came from our local community, which was populated with Russian, Tatar, and Jewish ethnic groups; there were no patients of African, Asian or Hispanic origin. Furthermore, AIS patients following endovascular interventions were excluded from the analysis due to a small number of observations. Although some clinical and imaging data were collected prospectively, the research was retrospective in nature. As a result, we were not blinded to the outcome. The study was also confined to a single clinical center; to cope with that bias, we used bootstrapping. Finally, prospective multicenter external validation would be desirable.


The HTI is a four-item tool composed of ASPECTS, NIHSS, HMCA sign, and presence of AF on ECG at admission. The total score ranges from zero to eight. The higher the score, the more likely HT can occur. Knowing probability of any HT in advance could exert a significant influence on decision making to reinstitute anticoagulation in AIS patients with AF. It is a simple yet reliable instrument to predict any HT within 2 weeks after onset of AIS in the MCA territory regardless of the use of IV rtPA.



Atrial fibrillation


Asymptomatic hemorrhagic transformation


Acute ischemic stroke


Alanine transaminase


Activated partial thromboplastin time


Alberta Stroke Program Early CT score


Aspartate transaminase


Area under the ROC curve


Bias-corrected and accelerated (adjusted for ties)


Binary logistic regression


Confirmatory factor analysis


Confidence interval


Computed tomography


Diastolic blood pressure


Derivation cohort


Diffusion-weighted imaging




Exploratory factor analysis


Fraction of missing information


Hemorrhagic infarction type 1


Hemorrhagic infarction type 2


Hyperdense middle cerebral artery


Heart rate


Hemorrhagic transformation


Hemorrhagic transformation index


Intraclass correlation coefficient


Internal consistency reliability


International normalized ratio


Interquartile range

IV rtPA:

Intravenous recombinant tissue plasminogen activator


Kernel density estimate


Locally weighted scatterplot smoothing


Left ventricular ejection fraction (echocardiography, Simpson method)


Middle cerebral artery


Markov chain Monte Carlo


Multiple imputation


Magnetic resonance imaging




National Institutes of Health Stroke Scale


Odds ratio


Principal component factor


Principal factor


Parenchymal hematoma type 1


Parenchymal hematoma type 2


Predicted probability


Prothrombin time


Red blood cells


Relative efficiency


Relative increase in variance


Receiver operating characteristic


Relative risk ratio


Systolic blood pressure


Standard error


Symptomatic hemorrhagic transformation


Transient ischemic attack


Validation cohort


White blood cells


  1. Lei C, Wu B, Liu M, Chen Y. Asymptomatic hemorrhagic transformation after acute ischemic stroke: is it clinically innocuous? J Stroke Cerebrovasc Dis. 2014;23:2767–72.

    Article  PubMed  Google Scholar 

  2. Sussman ES, Connolly ES Jr. Hemorrhagic transformation: a review of the rate of hemorrhage in the major clinical trials of acute ischemic stroke. Front Neurol. 2013;4:69.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Zhang J, Yang Y, Sun H, Xing Y. Hemorrhagic transformation after cerebral infarction: current concepts and challenges. Ann Transl Med. 2014;2:81.

    PubMed  PubMed Central  Google Scholar 

  4. Quain DA, Parsons MW, Loudfoot AR, Spratt NJ, Evans MK, Russell ML, et al. Improving access to acute stroke therapies: a controlled trial of organised pre-hospital and emergency care. Med J Aust. 2008;189:429–33.

    PubMed  Google Scholar 

  5. Lou M, Safdar A, Mehdiratta M, Kumar S, Schlaug G, Caplan L, et al. The HAT score: a simple grading scale for predicting hemorrhage after thrombolysis. Neurology. 2008;71:1417–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Cucchiara B, Tanne D, Levine SR, Demchuk AM, Kasner S. A risk score to predict intracranial hemorrhage after recombinant tissue plasminogen activator for acute ischemic stroke. J Stroke Cerebrovasc Dis. 2008;17:331–3.

    Article  PubMed  Google Scholar 

  7. Mazya M, Egido JA, Ford GA, Lees KR, Mikulik R, Toni D, et al. Predicting the risk of symptomatic intracerebral hemorrhage in ischemic stroke treated with intravenous alteplase: Safe Implementation of Treatments in Stroke (SITS) symptomatic intracerebral hemorrhage risk score. Stroke 2012;43:1524–31.

    Article  CAS  PubMed  Google Scholar 

  8. Strbian D, Engelter S, Michel P, Meretoja A, Sekoranja L, Ahlhelm FJ, et al. Symptomatic intracranial hemorrhage after stroke thrombolysis: the SEDAN score. Ann Neurol. 2012;71:634–41.

    Article  PubMed  Google Scholar 

  9. Menon BK, Saver JL, Prabhakaran S, Reeves M, Liang L, Olson DM, et al. Risk score for intracranial hemorrhage in patients with acute ischemic stroke treated with intravenous tissue-type plasminogen activator. Stroke. 2012;43:2293–9.

    Article  CAS  PubMed  Google Scholar 

  10. Heidbuchel H, Verhamme P, Alings M, Antz M, Diener HC, Hacke W, et al. Updated European Heart Rhythm Association Practical Guide on the use of non-vitamin K antagonist anticoagulants in patients with non-valvular atrial fibrillation. Europace. 2015;17:1467–507.

    Article  PubMed  Google Scholar 

  11. Koo CK, Teasdale E, Muir KW. What constitutes a true hyperdense middle cerebral artery sign? Cerebrovasc Dis. 2000;10:419–23.

    Article  CAS  PubMed  Google Scholar 

  12. 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–7.

    Article  Google Scholar 

  13. Fiorelli M, Bastianello S, von Kummer R, del Zoppo GJ, Larrue V, Lesaffre E, et al. Hemorrhagic transformation within 36 hours of a cerebral infarct: relationships with early clinical deterioration and 3-month outcome in the European Cooperative Acute Stroke Study I (ECASS I) cohort. Stroke. 1999;30:2280–4.

    Article  CAS  PubMed  Google Scholar 

  14. Gadermann AM, Guhn M, Bruno DZ. Estimating ordinal reliability for Likert-type and ordinal item response data: a conceptual, empirical, and practical guide. Practical Assessment, Research & Evaluation. 2012;17:3.

    Google Scholar 

  15. Trizano-Hermosilla I, Alvarado JM. Best alternatives to Cronbach's alpha reliability in realistic conditions: congeneric and asymmetrical measurements. Front Psychol. 2016;7:769.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Green SB, Yang Y. Reliability of summed item scores using structural equation modeling: an alternative to coefficient alpha. Psychometrika. 2009;74:155–67.

    Article  Google Scholar 

  17. DeLong ER, DeLong DM, Clarke-Pearson DL. Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biometrics. 1988;44:837–45.

    Article  CAS  PubMed  Google Scholar 

  18. Neumann-Haefelin T, Hoelig S, Berkefeld J, Fiehler J, Gass A, Humpich M, et al. Leukoaraiosis is a risk factor for symptomatic intracerebral hemorrhage after thrombolysis for acute stroke. Stroke. 2006;37:2463–6.

    Article  PubMed  Google Scholar 

  19. Marsh EB, Llinas RH, Schneider AL, Hillis AE, Lawrence E, Dziedzic P, et al. Predicting hemorrhagic transformation of acute ischemic stroke: prospective validation of the HeRS score. Medicine (Baltimore). 2016;95:e2430.

    Article  CAS  Google Scholar 

  20. Tan S, Wang D, Liu M, Zhang S, Wu B, Liu B. Frequency and predictors of spontaneous hemorrhagic transformation in ischemic stroke and its association with prognosis. J Neurol. 2014;261:905–12.

    Article  PubMed  Google Scholar 

  21. Zou M, Churilov L, He A, Campbell B, Davis SM, Yan B. Hyperdense middle cerebral artery sign is associated with increased risk of hemorrhagic transformation after intravenous thrombolysis for patients with acute ischaemic stroke. J Clin Neurosci. 2013;20:984–7.

    Article  PubMed  Google Scholar 

  22. Terruso V, D’Amelio M, Di Benedetto N, Lupo I, Saia V, Famoso G, et al. Frequency and determinants for hemorrhagic transformation of cerebral infarction. Neuroepidemiology. 2009;33:261–5.

    Article  PubMed  Google Scholar 

  23. Paciaroni M, Agnelli G, Corea F, Ageno W, Alberti A, Lanari A, et al. Early hemorrhagic transformation of brain infarction: rate, predictive factors, and influence on clinical outcome: results of a prospective multicenter study. Stroke. 2008;39:2249–56.

    Article  PubMed  Google Scholar 

  24. Pundik S, McWilliams-Dunnigan L, Blackham KL, Kirchner HL, Sundararajan S, Sunshine JL, et al. Older age does not increase risk of hemorrhagic complications after intravenous and/or intra-arterial thrombolysis for acute stroke. J Stroke Cerebrovasc Dis. 2008;17:266–72.

    Article  PubMed  Google Scholar 

  25. Demchuk AM, Khan F, Hill MD, Barber PA, Silver B, Patel S, et al. Importance of leukoaraiosis on CT for tissue plasminogen activator decision making: evaluation of the NINDS rt-PA stroke study. Cerebrovasc Dis. 2008;26(2):120–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Nardi K, Leys D, Eusebi P, Cordonnier C, Gautier S, Hénon H, et al. Influence of lipid profiles on the risk of hemorrhagic transformation after ischemic stroke: systematic review. Cerebrovasc Dis Extra. 2011;1:130–41.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Kaesmacher J, Kaesmacher M, Maegerlein C, Zimmer C, Gersing AS, Wunderlich S, et al. Hemorrhagic transformations after thrombectomy: risk factors and clinical relevance. Cerebrovasc Dis. 2017;43:294–304.

    Article  CAS  PubMed  Google Scholar 

  28. Sandset EC, Berge E, Kjeldsen SE, Julius S, Holzhauer B, Krarup LH, et al. Heart rate as a predictor of stroke in high-risk, hypertensive patients with previous stroke or transient ischemic attack. J Stroke Cerebrovasc Dis. 2014;23:2814–8.

    Article  PubMed  Google Scholar 

  29. Chen CH, Huang PW, Tang SC, Shieh JS, Lai DM, Wu AY, et al. Complexity of heart rate variability can predict stroke-in-evolution in acute ischemic stroke patients. Sci Rep. 2015;5:17552.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Tatlisumak T, Fisher M. Hematologic disorders associated with ischemic stroke. J Neurol Sci. 1996;140:1–11.

    Article  CAS  PubMed  Google Scholar 

  31. Marsh EB, Llinas RH, Hillis AE, Gottesman RF. Hemorrhagic transformation in patients with acute ischaemic stroke and an indication for anticoagulation. Eur J Neurol. 2013;20:962–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Saposnik G, Fang J, Kapral MK, Tu JV, Mamdani M, Austin P, et al. The iScore predicts effectiveness of thrombolytic therapy for acute ischemic stroke. Stroke. 2012;43:1315–22.

    Article  PubMed  Google Scholar 

  33. Saposnik G, Guzik AK, Reeves M, Ovbiagele B, Johnston SC. Stroke prognostication using age and NIH Stroke Scale: SPAN-100. Neurology. 2013;80:21–8.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Yong AG, Pearce S. A Beginner’s guide to factor analysis: focusing on exploratory factor analysis. Tutor Quant Methods Psychol. 2013;9:79–94.

    Article  Google Scholar 

  35. Asuzu D, Nystrom K, Amin H, Schindler J, Wira C, Greer D, et al. Comparison of 8 scores for predicting symptomatic intracerebral hemorrhage after IV thrombolysis. Neurocrit Care. 2015;22:229–33.

    Article  CAS  PubMed  Google Scholar 

  36. Sung SF, Chen SC, Lin HJ, Chen YW, Tseng MC, Chen CH. Comparison of risk-scoring systems in predicting symptomatic intracerebral hemorrhage after intravenous thrombolysis. Stroke. 2013;44:1561–6.

    Article  CAS  PubMed  Google Scholar 

  37. Demeestere J, Garcia-Esperon C, Garcia-Bermejo P, Ombelet F, McElduff P, Bivard A, et al. Evaluation of hyperacute infarct volume using ASPECTS and brain CT perfusion core volume. Neurology. 2017;88:2248–53.

    Article  PubMed  PubMed Central  Google Scholar 

  38. González RG, Schaefer PW, Buonanno FS, Schwamm LH, Budzik RF, Rordorf G, et al. Diffusion-weighted MR imaging: diagnostic accuracy in patients imaged within 6 hours of stroke symptom onset. Radiology. 1999;210:155–62.

    Article  PubMed  Google Scholar 

  39. Schröder J, Cheng B, Ebinger M, Köhrmann M, Wu O, Kang DW, et al. Validity of acute stroke lesion volume estimation by diffusion-weighted imaging-Alberta Stroke Program Early Computed Tomographic Score depends on lesion location in 496 patients with middle cerebral artery stroke. Stroke. 2014;45:3583–8.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Merwick Á, Werring D. Posterior circulation ischaemic stroke. BMJ. 2014;348:g3175.

    Article  PubMed  Google Scholar 

  41. Sarikaya H, Arnold M, Engelter ST, Lyrer PA, Mattle HP, Georgiadis D, et al. Outcomes of intravenous thrombolysis in posterior versus anterior circulation stroke. Stroke. 2011;42:2498–502.

    Article  PubMed  Google Scholar 

  42. Libman RB, Kwiatkowski TG, Hansen MD, Clarke WR, Woolson RF, Adams HP. Differences between anterior and posterior circulation stroke in TOAST. Cerebrovasc Dis. 2001;11:311–6.

    Article  CAS  PubMed  Google Scholar 

  43. Tong X, Liao X, Pan Y, Cao Y, Wang C, Liu L, et al. Intravenous thrombolysis is more safe and effective for posterior circulation stroke: data from the Thrombolysis Implementation and Monitor of Acute Ischemic Stroke in China (TIMS-China). Medicine (Baltimore). 2016;95:e3848.

  44. Oppenheim C, Stanescu R, Dormont D, Crozier S, Marro B, Samson Y, et al. False-negative diffusion-weighted MR findings in acute ischemic stroke. AJNR Am J Neuroradiol. 2000;21:1434–40.

    CAS  PubMed  Google Scholar 

  45. Puetz V, Sylaja PN, Coutts SB, Hill MD, Dzialowski I, Mueller P, et al. Extent of hypoattenuation on CT angiography source images predicts functional outcome in patients with basilar artery occlusion. Stroke. 2008;39:2485–90.

    Article  PubMed  Google Scholar 

  46. Pisters R, Lane DA, Nieuwlaat R, de Vos CB, Crijns HJ, Lip GY. A novel user-friendly score (HAS-BLED) to assess 1-year risk of major bleeding in patients with atrial fibrillation: the Euro Heart Survey. Chest. 2010;138:1093–100.

    Article  PubMed  Google Scholar 

Download references


Not applicable.


No funding was received.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Author information

Authors and Affiliations



M.N.K. was involved in the design, data collection, analysis and interpretation, and reporting of this project. D.R.K. and M.M.I. helped with the design, data collection, data analysis, and revision of the manuscript. All authors have read and approved the final manuscript for publication.

Corresponding author

Correspondence to Mikhail N. Kalinin.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Local Ethics Committee of the Kazan State Medical University, Kazan, Russia. Informed consent was not required since the study was observational and retrospective in nature.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

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 distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kalinin, M.N., Khasanova, D.R. & Ibatullin, M.M. The hemorrhagic transformation index score: a prediction tool in middle cerebral artery ischemic stroke. BMC Neurol 17, 177 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: