Moving Beyond a Single Perfusion Threshold to Define Penumbra
A Novel Probabilistic Mismatch Definition

Abstract
Background and Purpose—The mismatch lesion volumes defined by perfusion-weighted imaging exceeding diffusion-weighted imaging have been used as a marker of ischemic penumbral tissue. Defining the perfusion lesion by thresholding has shown promise as a practical tool; several positron emission tomography studies have indicated a more probabilistic relationship between perfusion and infarction. Here, we used a randomized controlled trial dataset of tissue-type plasminogen activator 3 to 6 hours after stroke to: (1) quantify the relationship between severity of hypoperfusion (measured by Tmax) and risk of infarction; (2) exploit this relationship to present a novel definition of mismatch based on infarct probabilities rather than dichotomies; and (3) examine the treatment response in the subgroup of patients with mismatch by the new definition.
Methods—Patients from the Echoplanar Imaging Thrombolytic Evaluation Trial (EPITHET) were included. Baseline perfusion-weighted imaging and 90-day T2-weighted imaging were coregistered. Perfusion-weighted imaging lesion volumes were divided into 10 Tmax delay strata, and infarct risk was defined as the fraction of the tissue at a given Tmax strata that progressed to infarction by day 90.
Results—Sixty-two patients were studied. Infarct risk was an increasing function of Tmax for all subgroups, including the whole cohort. The probabilistic approach outperformed all Tmax thresholds, with exception of the Tmax ≥10 threshold, for which it was only favored by a trend.
Conclusions—Infarct risk and treatment effect increased with severity of perfusion abnormalities. This suggests that a severity-weighted mismatch definition may define penumbral tissue more accurately.
- magnetic resonance imaging
- mismatch
- penumbra
- perfusion-weighted imaging
- tissue-type plasminogen activator
Introduction
Ischemic penumbral tissue is at risk of infarction but has the potential to be salvaged by reperfusion and/or other strategies.1 With the advent of functional MRI techniques, the concept of mismatch defined by perfusion-weighted imaging (PWI) exceeding diffusion-weighted imaging (DWI) lesion volumes was postulated as a penumbral marker.2 Patients with PWI/DWI mismatch are likely to benefit from timely reperfusion and may be eligible for thrombolytic therapy beyond the current therapeutic time window.3 A recent meta-analysis showed that MRI-selected patients treated with thrombolysis beyond 3 hours had improved clinical outcomes with reperfusion/recanalization.4
In acute stroke, PWI is used to estimate the tissue regions that will infarct if hypoperfusion persists and reperfusion is not achieved. Positron emission tomography studies have focused on determining thresholds for the separation between benign oligemia and penumbra, the so-called penumbral thresholds. Based on these findings, several MR studies have aimed at identifying parameters and thresholds that correspond best with this penumbral threshold.5–7 In practice, the hypoperfused at-risk region is most commonly defined as a region above a certain threshold value of mean transit time or time to peak of the residue function (Tmax).8 It is clear that lower threshold values may include significant amounts of benign oligemia,9 and that no single threshold is capable of accurately predicting tissue fate, even in homogenous groups of patients without reperfusion.10 Therefore, the actual choice of threshold is inherently a compromise between sensitivity and specificity for infarction. Because the threshold-based approach to penumbra definition effectively dichotomizes the perfusion values as normal or abnormal, it masks the distribution of perfusion severities in the individual patient and oversimplifies any probabilistic relationship between hypoperfusion severity and risk of infarction. Figure 1 shows as an example how 2 similar threshold-based volumes (similar Tmax ≥2-second volumes) can mask 2 different patterns of perfusion severities, as seen by the amount of tissue affected by higher Tmax values.
Illustrative example of different Tmax lesion profiles. Patients A and B had comparable Tmax≥2 lesion volumes (229 and 255 mL, respectively); however, when examining the Tmax severity profile, patient A had larger amounts of very prolonged Tmax (see color codes) than patient B. Patient A had an acute diffusion-weighted imaging (DWI) lesion of 145 mL expanding to 188 mL, whereas patient B presented with a 3-mL DWI lesion expanding to 5 mL on day 90. The severity profile is for all slices, but only 1 representative slice for each patient is shown here.
In this study, we used data from a phase II randomized controlled trial of tissue-type plasminogen activator with MRI-based outcome measures to: (1) quantify the relationship between severity of hypoperfusion (as measured by Tmax) and risk of infarction; (2) exploit this relationship to present a novel definition of mismatch based on infarct probabilities rather than dichotomies; and (3) examine the treatment response in the subgroup of patients with mismatch by the new definition.
Materials and Methods
Patients and Imaging Protocol
The study design of the Echoplanar Imaging Thrombolytic Evaluation Trial (EPITHET), including patient eligibility, treatment allocation, and imaging protocol, has been previously reported.11 Briefly, patients with acute hemispheric ischemic stroke who presented 3 to 6 hours after symptom onset were randomized to intravenous alteplase or placebo, and had DWI, PWI, and MRA sequences with 1.5-Tesla echoplanar-equipped MRI scanners before treatment and repeated on days 3 to 5. On day 90, follow-up T2-weighted images were obtained.
Quantification of Infarct Risk Versus Tmax Severity
To quantify the relationship between severity of hypoperfusion and risk of infarction, we coregistered the original EPITHET trial day 90 final infarct regions of interest to the acute PWI images. Patients who had parenchymal hemorrhage, defined by the European Cooperative Acute Stroke Study (ECASS) algorithm,12,13 and those who had inadequate imaging quality for coregistration were excluded from this analysis.
Automatic registration (autoreg software; The McConnell Brain Imaging Centre) was used to register the day 90 T2-weighted imaging to the baseline PWI. If automatic registration failed, then a manual landmark-based registration was used to initialize another automatic registration attempt. Final acceptance of the coregistration quality was subject to a consensus-based quality control by 3 investigators (Y.N., S.C., and T.O.) by interactively blending the coregistered images to verify correct superposition.
Infarct Risks in Each Tmax Strata
The Tmax maps were stratified into 9 2-second strata as well as a ≥20-second strata, as depicted in Figure 1. In the individual patient, the infarct risk for each Tmax strata was then defined as the fraction of the strata with infarction on follow-up, as determined from the coregistered day 90 infarct regions of interest. This analysis was restricted to a manually edited region of interest containing only tissue with Tmax ≥2 seconds and excluding obvious false-positive regions. Regions that were DWI-positive were included in the analysis.
We investigated the infarct risks in all patients as well as in subgroups defined by treatment versus placebo and by recanalization versus no recanalization. As per the EPITHET protocol, recanalization was defined as improvement from baseline to day 3 to 5 arterial obstruction by ≥2 points, based on an adaptation of the Thrombolysis In Myocardial Infarction grading on MRA (0=complete occlusion, 1=severe stenosis, 2=mild to moderate stenosis, and 3=normal arterial caliber).14
Prediction of Final Lesion Volumes
We used the median infarct risks from nonrecanalized patients as a severity-weighted predictive model for the relationship between Tmax severity and final infarct risk. In each patient, the predicted infarct volume was defined as the sum of the volume of each Tmax strata (VolumeTmax) multiplied by the median infarct risks of that strata (RiskTmax):
We used a jack-knifing technique to avoid using the same data for model estimation as for performance validation. The predicted infarct volume for each patient was estimated using a model trained on the 10 remaining patients. Cycling through the patients in this way, we estimated these coefficients of RiskTmax 11 times, each time applying the model to the patient left out of the coefficient estimation. The severity-weighted infarct volume predictions were compared with those obtained by Tmax thresholding (Tmax 2 to Tmax 14 seconds). Analysis was restricted to nonrecanalized patients.
Severity-Weighted Mismatch Definition
Finally, we applied the Tmax severity versus risk model to the entire EPITHET cohort of 101 patients and classified mismatch patients using this new definition. A severity-weighted mismatch was defined as predicted infarct volume/DWI >1.2 and absolute difference >10 mL. We examined whether infarct growth was attenuated by alteplase in patients with mismatch by the new definition.
Statistical Analysis
Statistical analyses were performed using Stata/IC 10 (StataCorp). The 2-sample Wilcoxon rank-sum test was used for continuous variables, and Fisher exact test was used for categorical variables between patients with alteplase and placebo and between recanalization and no recanalization. The relationship between Tmax severity and infarct risk was quantified by a multilevel random-effect linear model with patients as a grouping variable. To compare infarct risks by alteplase versus placebo and by recanalization versus no recanalization, multilevel random-effect linear models with appropriate interaction terms and patients as a grouping variable were utilized. As per the original EPITHET analysis protocol, the primary outcome measure (infarct growth attenuation in mismatch patients between alteplase and placebo assessed by a t test for ratio of the geometric means) as well as secondary outcome measures (relative growth, absolute growth, and the difference of cube-root–transformed volume change) were assessed.
Results
Out of 101 patients enrolled in EPITHET, 99 patients had baseline PWI and 72 had both baseline PWI and day 90 T2-weighted imaging (day 92 interquartile range, 90–95). Six patients with parenchymal hematoma and 4 patients with severely distorted EPI imaging were excluded from the assessment of Tmax stratified infarct risk. For the remaining 62 patients (29 received alteplase and 33 received placebo), baseline PWI and day 90 T2-weighted imaging were used for this analysis. Recanalization status could be assessed in 29 patients with adequate pretreatment and posttreatment MRA scans. Of these, recanalization was present in 18 patients. The Table shows baseline characteristics, and there was no significant difference between the alteplase and placebo groups. Also, the baseline risk factors and imaging characteristics did not differ between recanalization and no recanalization groups, except for prevalence of current or past smoking and median DWI lesion volumes.
Baseline Characteristics
Infarct Risks in Each Tmax Strata
Figure 2A shows the volume of tissue within specified Tmax strata on baseline PWI and the volume of infarcted tissue within specified Tmax strata on day 90 T2-weighted imaging for all 62 patients. Infarct risks in each Tmax strata increased progressively with the Tmax values (Figure 2B). A strong relationship between the infarct risks and severity of Tmax values was observed (increase of 2.2% per unit of Tmax; 95% CI, 2.0%–2.4%; P<0.001; intraclass correlation coefficient=0.65).
Baseline tissue volume and infarcted tissue volume within specified Tmax strata (A) and infarct risks in each Tmax strata (B) for all patients (n=62). Boxes indicate median in A and interquatile ranges in B. Bars indicate maximum and minimum values, and lines in the boxes and numbers indicate median values.
The positive relationship between infarct risks and Tmax severity was maintained in subgroups of treatment (Figure 3A) and recanalization status (Figure 3B). Corresponding increases in risk differed significantly between alteplase and placebo subgroups (interaction Tmax by treatment, P<0.001) and between recanalization and no recanalization subgroups (interaction Tmax by recanalization, P<0.001). Patients without recanalization and those treated by placebo had high infarct risks, particularly in severe Tmax strata.
A, Infarct risks by treatment arm (tissue-type plasminogen activator=29; placebo=33). B, Infarct risks by recanalization status (recanalization=18; no recanalization=11). Bars indicate maximum and minimum values, boxes indicate interquatile ranges, and lines and numbers in the boxes indicate median values.
Prediction of Final Lesion Volumes
Figure 4 shows the Tmax versus severity relationships in the 11 patients with no recanalization. To assess the accuracy of infarct volume prediction, we calculated the differences (residuals) between baseline PWI volumes and 90-day infarct volumes for conventional Tmax thresholds in 11 patients with no recanalization (Figure 5A). Low Tmax thresholds, such as ≥2, ≥4, and ≥6 seconds, overestimated the 90-day infarct volumes, whereas high thresholds, such as ≥12 and ≥14, underestimated the same volumes. Tmax thresholds of ≥8 and ≥10 seconds provided the most accurate estimation of the 90-day infarct volumes (median differences of 31 and −35 mL, respectively). The Tmax ≥10 estimates did not differ significantly from the severity-weighted estimate (P=0.2), but all other Tmax estimates showed significantly worse estimation of the final infarct with a maximum P=0.0068. For comparison, we also analyzed the DWI performance, which was outperformed by the severity-weighted prediction (P=0.032).
Infarct risks in each Tmax strata for the 11 patients with no recanalization. The plot highlights intrapatient variability. The infarct risk generally is increasing with Tmax in the individual, although often with a plateau beginning at approximately the 12 to 14 strata. The plot also highlights interpatient variability. The severity–risk relationship in the individual is subject to considerable variability across patients. It is clear from this plot that a threshold approach may not be the best way to model the continuous relationship between Tmax and risk of infarction.
A, Difference between the perfusion-weighted imaging (PWI) volumes and the 90-day infarct volumes for patients with no recanalization (n=11). Bars indicate maximum and minimum values, boxes indicate interquartile ranges, and lines in the boxes indicate median values. B, Prevalence of mismatch patients by various PWI definitions in all EPITHET patients (n=99). Mismatch: PWI ÷ DWI volume >1.2 and PWI − DWI volume ≥10 mL.
Severity-Weighted Mismatch Definition
When mismatch was assessed in all EPITHET patients using the severity-weighted predicted infarct volumes generated by the infarct risks in the nonrecanalized patients, the prevalence of the severity-weighted mismatch was 62% (61/99; Figure 5B). For the EPITHET primary outcome measure, the patients with severity-weighted mismatch showed trends to infarct growth attenuation with alteplase compared with placebo, such as geometric mean (P=0.2681), relative growth (P=0.0716), absolute growth (P=0.0873), and difference in cube-root lesion volume (P=0.0563), although it did not reach statistical significance, perhaps because of the lower number of patients. In comparison, the original EPITHET mismatch definition using the threshold approach with Tmax ≥2 seconds that classified 86% (85/99) of patients with mismatch was also nonsignificant for all analytical methods of infarct growth, such as geometric mean (P=0.239), relative growth (P=0.054), absolute growth (P=0.126), and difference in cube-root lesion volume (P=0.069).
Discussion
The present study showed that the risk of infarction increased with the severity of perfusion deficits instead of showing the expected threshold, although across patients Tmax 8 to 10 seconds did exhibit a sharp slope (Figure 3B). The source of this probabilistic relationship appeared to be caused by both intrapatient and interpatient variability in the Tmax–infarct risk relationship. This suggests that no single PWI threshold valid across patients can perfectly distinguish between salvageable tissue and areas destined to infarct (Figure 4). Hence, the mismatch between single PWI thresholds and DWI lesion volumes may provide an oversimplified penumbra estimate. Even so, the calculation of the predicted infarct volume using the severity-weighted probabilistic approach is no more complex than the quantification of a threshold-defined mismatch.
Differentiation of true penumbra from tissue experiencing benign oligemia has been the major challenge to the PWI/DWI mismatch concept.9 In general, the unthresholded visual PWI/DWI mismatch volume is greater than the volume of penumbra defined by positron emission tomography.15–18 Obviously, higher PWI thresholds provide smaller PWI lesion volumes, so using a higher PWI threshold increases specificity, but at the cost of sensitivity. This has been shown in a number of studies, and these findings are consistent with the increasing severity–risk relationship.5,6,19 In our study, Tmax thresholds of ≥8 and ≥10 seconds appeared to best-predict the final infarct volume; however, the best thresholds for final infarction in terms of spatial sensitivity and specificity are likely lower than these numbers.19
A model to predict final infarction must be based on patients with no reperfusion. The clinically relevant question is what will happen to the tissue if there is no reperfusion, and this needs to be estimated in patients with no reperfusion. It has to be acknowledged that many patients have reperfusion spontaneously and will never have infarction to an extent predicted by the model. It is likely that further research can help identify this subset of patients by means of occlusion site or other reperfusion determining factors and, in this way, can improve the model for these patients. In the present study, we based the predictive model on patients without recanalization as a proxy for no reperfusion because the best reperfusion definition for this type of study is unknown. Interestingly, the infarct risks did not reach 100%, even for the highest Tmax range in the nonrecanalized group. We speculate that factors such as collateral supply in some cases inflate Tmax and do not cause infarction, which is consistent with recent findings that Tmax reflects a combination of hemodynamic effects.8 Visual inspection indicated that regions close to sinuses and edges were often associated with high Tmax values and absence of infarction. Lacking automatic ways to distinguish artifactual areas, and given that they often are of modest volume, we did not correct for these artifacts; however, future efforts to detect these artifactual regions are needed. Patients who died before the day 90 follow-up imaging represent a possible data bias. Another potential bias is tissue shrinkage by day 90, causing underestimates of the volume of infarction and thereby underestimation of risk for a given Tmax value. Future studies may be able to address this with earlier imaging follow-up, limiting the influence of atrophy.
The probabilistic approach to infarct prediction used in this study is not novel and positron emission tomography studies have demonstrated a probabilistic, rather than threshold-based, relationship between acute perfusion and tissue infarction.20 Several methods have been proposed that integrate information from several modalities and use a variety of transformation techniques to estimate the infarct probability.21–23 However, these techniques still have not found widespread use despite their theoretical appeal. The method presented here only uses DWI and PWI to classify a mismatch and is simple to implement. It does require exclusion of false-positive areas, in our case performed by a manual exclusion of false-positives, which is relatively quick to perform (3–5 minutes). The actual calculation of the predicted volume is then a few algebraic operations of microseconds in duration. To the best of our knowledge, a severity-weighted mismatch definition has not been presented before. In addition, a map of infarct risk can be generated by substituting each Tmax strata with the associated probability to produce a visual topographical representation of infarct risk in the individual. In its current form, the method does not require coregistration to calculate a volumetric mismatch; however, as fully automated registration is becoming increasingly available, it is likely that registering the acute DWI to the PWI could provide more accurate mismatch estimates that account for DWI regions already reperfused. With coregistered data, it would also be possible to attribute higher infarct risks to tissue already DWI-positive, irrespective of the concurrent Tmax value.
There are some limitations to this study. First, the analyzed cohort was small and we had a limited number of patients without recanalization. The estimated probability of infarction in each Tmax strata could be improved by a larger cohort with acute PWI and follow-up images, and it would enhance the predicted infarct volume by the severity-weighted approach. Second, regarding all perfusion parameters, Tmax is sensitive to several hemodynamic effects and this is bound to introduce some variability in the risk estimates.8 Third, our estimates of risk will be affected by fluctuating perfusion levels. Ideally, this analysis could be improved by limiting it to nonreperfusing tissue as determined by coregistration of acute and follow-up perfusion imaging. Unfortunately, they day 3 perfusion reperfusion assessment in these data makes this impossible for patients with edema, which peaks at day 3 to 5 and is largest in patients without reperfusion, which is exactly the group of interest. Earlier reperfusion assessments, at the end of infarct evolution, but before severe edema development, may partially address this issue. Limiting the at-risk region to nonreperfused tissue, albeit limited by coregistration errors, did not significantly improve prediction of final infarct volume over those based on Tmax ≥8-second estimates. Fourth, the mismatch calculation in this study assumes that DWI is a reasonable representation of a core matter that remains controversial, but this topic previously has been systematically examined in this data set and has been found to have negligible impact on mismatch classifications.24
Despite these limitations, this study demonstrated that the probability of infarction is strongly tied to the severity of PWI deficits. A penumbral definition based on Tmax ≥2 is now considered inadequate because of a penumbral definition and our results corroborating this (showing that Tmax 2 seconds to Tmax 4 seconds have a median risk of only 8%); also, in the absence of more severe areas of hypoperfusion, a penumbral definition of 2 seconds would grossly overestimate the area of final infarction.25
In summary, we have shown that infarct risks increased with the severity of perfusion deficits and that in a small population of 11 patients a probabilistic method for infarct volume prediction outperforms all Tmax thresholds statistically, with exception of the Tmax ≥10 threshold, for which it was only favored by a trend. Thus, although we do not provide formal evidence that the proposed method performs better than using a Tmax ≥10 threshold, the data suggest it is of potential interest and deserves further investigation in larger samples. It would seem ideal to perform both the reperfusion and final infarct assessments at an earlier time point to minimize the confounding effects of late futile reperfusion and infarct shrinkage.
Sources of Funding
Supported by National Health and Medical Research Council, Australia; National Stroke Foundation, Australia; and Heart Foundation of Australia.
Disclosures
None.
- Received November 3, 2011.
- Accepted January 18, 2012.
- © 2012 American Heart Association, Inc.
References
- 1.↵
- Donnan GA,
- Baron JC,
- Davis SM,
- Sharp FR
- 2.↵
- 3.↵
- Albers GW
- 4.↵
- Mishra NK,
- Albers GW,
- Davis SM,
- Donnan GA,
- Furlan AJ,
- Hacke W,
- et al
- 5.↵
- Shih LC,
- Saver JL,
- Alger JR,
- Starkman S,
- Leary MC,
- Vinuela F,
- et al
- 6.↵
- Butcher KS,
- Parsons M,
- MacGregor L,
- Barber PA,
- Chalk J,
- Bladin C,
- et al
- 7.↵
- 8.↵
- Calamante F,
- Christensen S,
- Desmond PM,
- Ostergaard L,
- Davis SM,
- Connelly A
- 9.↵
- Kidwell CS,
- Alger JR,
- Saver JL
- 10.↵
- Christensen S,
- Mouridsen K,
- Wu O,
- Hjort N,
- Karstoft H,
- Thomalla G,
- et al
- 11.↵
- 12.↵
- Berger C,
- Fiorelli M,
- Steiner T,
- Schäbitz WR,
- Bozzao L,
- Bluhmki E,
- et al
- 13.↵
- Larrue V,
- von Kummer RR,
- Müller A,
- Bluhmki E
- 14.↵
- Hacke W,
- Albers G,
- Al-Rawi Y,
- Bogousslavsky J,
- Davalos A,
- Eliasziw M,
- et al
- 15.↵
- Neumann-Haefelin T,
- Wittsack HJ,
- Wenserski F,
- Siebler M,
- Seitz RJ,
- Modder U,
- et al
- 16.↵
- Parsons MW,
- Yang Q,
- Barber PA,
- Darby DG,
- Desmond PM,
- Gerraty RP,
- et al
- 17.↵
- Thijs VN,
- Adami A,
- Neumann-Haefelin T,
- Moseley ME,
- Marks MP,
- Albers GW
- 18.↵
- Sobesky J,
- Zaro Weber O,
- Lehnhardt FG,
- Hesselmann V,
- Neveling M,
- Jacobs A,
- et al
- 19.↵
- Olivot JM,
- Mlynash M,
- Thijs VN,
- Kemp S,
- Lansberg MG,
- Wechsler L,
- et al
- 20.↵
- 21.↵
- 22.↵
- Mitsias PD,
- Jacobs MA,
- Hammoud R,
- Pasnoor M,
- Santhakumar S,
- Papamitsakis NI,
- et al
- 23.↵
- Wu O,
- Christensen S,
- Hjort N,
- Dijkhuizen RM,
- Kucinski T,
- Fiehler J,
- et al
- 24.↵
- Chemmanam T,
- Campbell BC,
- Christensen S,
- Nagakane Y,
- Desmond PM,
- Bladin CF,
- et al
- 25.↵
- Takasawa M,
- Jones PS,
- Guadagno JV,
- Christensen S,
- Fryer TD,
- Harding S,
- et al
Jump to
This Issue
Article Tools
- Moving Beyond a Single Perfusion Threshold to Define PenumbraYoshinari Nagakane, Soren Christensen, Toshiyasu Ogata, Leonid Churilov, Henry Ma, Mark W. Parsons, Patricia M. Desmond, Christopher R. Levi, Kenneth S. Butcher, Stephen M. Davis and Geoffrey A. Donnan for the EPITHET InvestigatorsStroke. 2012;43:1548-1555, originally published May 25, 2012https://doi.org/10.1161/STROKEAHA.111.643932
Citation Manager Formats
Share this Article
- Moving Beyond a Single Perfusion Threshold to Define PenumbraYoshinari Nagakane, Soren Christensen, Toshiyasu Ogata, Leonid Churilov, Henry Ma, Mark W. Parsons, Patricia M. Desmond, Christopher R. Levi, Kenneth S. Butcher, Stephen M. Davis and Geoffrey A. Donnan for the EPITHET InvestigatorsStroke. 2012;43:1548-1555, originally published May 25, 2012https://doi.org/10.1161/STROKEAHA.111.643932












