Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

DWI Intensity Values Predict FLAIR Lesions in Acute Ischemic Stroke

  • Vince I. Madai,

    Affiliations Center for Stroke Research Berlin (CSB), Charité-Universitätsmedizin, Berlin, Germany, Department of Neurology, Charité-Universtitätsmedizin, Berlin, Germany

  • Ivana Galinovic,

    Affiliation Center for Stroke Research Berlin (CSB), Charité-Universitätsmedizin, Berlin, Germany

  • Ulrike Grittner,

    Affiliations Center for Stroke Research Berlin (CSB), Charité-Universitätsmedizin, Berlin, Germany, Department for Biostatistics and Clinical Epidemiology, Charité-Universitätsmedizin, Berlin, Germany

  • Olivier Zaro-Weber,

    Affiliations Center for Stroke Research Berlin (CSB), Charité-Universitätsmedizin, Berlin, Germany, Max-Planck-Institute for Neurological Research, Cologne, Germany

  • Alice Schneider,

    Affiliations Center for Stroke Research Berlin (CSB), Charité-Universitätsmedizin, Berlin, Germany, Department for Biostatistics and Clinical Epidemiology, Charité-Universitätsmedizin, Berlin, Germany

  • Steve Z. Martin,

    Affiliation Center for Stroke Research Berlin (CSB), Charité-Universitätsmedizin, Berlin, Germany

  • Federico C. v. Samson-Himmelstjerna,

    Affiliations Center for Stroke Research Berlin (CSB), Charité-Universitätsmedizin, Berlin, Germany, Fraunhofer MEVIS, Bremen, Germany

  • Katharina L. Stengl,

    Affiliations Center for Stroke Research Berlin (CSB), Charité-Universitätsmedizin, Berlin, Germany, Department of Neurology, Charité-Universtitätsmedizin, Berlin, Germany

  • Matthias A. Mutke,

    Affiliations Center for Stroke Research Berlin (CSB), Charité-Universitätsmedizin, Berlin, Germany, Department of Neurology, Charité-Universtitätsmedizin, Berlin, Germany

  • Walter Moeller-Hartmann,

    Affiliation Department of Radiology, Krankenhaus Ludmillenstift, Meppen, Germany

  • Martin Ebinger,

    Affiliations Center for Stroke Research Berlin (CSB), Charité-Universitätsmedizin, Berlin, Germany, Department of Neurology, Charité-Universtitätsmedizin, Berlin, Germany

  • Jochen B. Fiebach,

    Affiliation Center for Stroke Research Berlin (CSB), Charité-Universitätsmedizin, Berlin, Germany

  • Jan Sobesky

    jan.sobesky@charite.de

    Affiliations Center for Stroke Research Berlin (CSB), Charité-Universitätsmedizin, Berlin, Germany, Department of Neurology, Charité-Universtitätsmedizin, Berlin, Germany

Abstract

Background and Purpose

In acute stroke, the DWI-FLAIR mismatch allows for the allocation of patients to the thrombolysis window (<4.5 hours). FLAIR-lesions, however, may be challenging to assess. In comparison, DWI may be a useful bio-marker owing to high lesion contrast. We investigated the performance of a relative DWI signal intensity (rSI) threshold to predict the presence of FLAIR-lesions in acute stroke and analyzed its association with time-from-stroke-onset.

Methods

In a retrospective, dual-center MR-imaging study we included patients with acute stroke and time-from-stroke-onset ≤12 hours (group A: n = 49, 1.5T; group B: n = 48, 3T). DW- and FLAIR-images were coregistered. The largest lesion extent in DWI defined the slice for further analysis. FLAIR-lesions were identified by 3 raters, delineated as regions-of-interest (ROIs) and copied on the DW-images. Circular ROIs were placed within the DWI-lesion and labeled according to the FLAIR-pattern (FLAIR+ or FLAIR−). ROI-values were normalized to the unaffected hemisphere. Adjusted and nonadjusted receiver-operating-characteristics (ROC) curve analysis on patient level was performed to analyze the ability of a DWI- and ADC-rSI threshold to predict the presence of FLAIR-lesions. Spearman correlation and adjusted linear regression analysis was performed to assess the relationship between DWI-intensity and time-from-stroke-onset.

Results

DWI-rSI performed well in predicting lesions in FLAIR-imaging (mean area under the curve (AUC): group A: 0.84; group B: 0.85). An optimal mean DWI-rSI threshold was identified (A: 162%; B: 161%). ADC-maps performed worse (mean AUC: A: 0.58; B: 0.77). Adjusted regression models confirmed the superior performance of DWI-rSI. Correlation coefficents and linear regression showed a good association with time-from-stroke-onset for DWI-rSI, but not for ADC-rSI.

Conclusion

An easily assessable DWI-rSI threshold identifies the presence of lesions in FLAIR-imaging with good accuracy and is associated with time-from-stroke-onset in acute stroke. This finding underlines the potential of a DWI-rSI threshold as a marker of lesion age.

Introduction

In patients with acute ischemic stroke, the combination of a hyperintense lesion in diffusion weighted imaging (DWI) and the absence of a corresponding lesion in T2-weighted fluid-attenuated inversion recovery (FLAIR) imaging, the so called DWI-FLAIR mismatch, can predict the time from stroke onset <4.5 h [1]. This can be attributed to the time-dependent appearance of FLAIR-lesions within the first hours after stroke onset [2][4]. This finding is promising, as identification of patients eligible for thrombolysis with unknown stroke onset, e.g. in wake-up stroke, may be facilitated by specific imaging markers. However, the visual assessment of FLAIR lesions may be difficult [5], [6] and the automated analysis of FLAIR images is challenging owing to low contrast and partial volume effects [7]. Lesions on DW-images, on the other hand, show high contrast and can be easily delineated by automated software solutions [8] making DWI suitable for clinical stroke trials. In the present work, we hypothesized that DWI signal-intensity increases with the time-from-stroke-onset. To test this hypothesis we investigated, whether relative DWI intensity (DWI-rSI) is associated a) with the presence of hyperintensities in FLAIR-imaging and b) with time-from-stroke-onset.

Materials and Methods

Ethics Statement

All patients gave informed written consent prior to the study. The study was conducted according to the principles expressed in the Declaration of Helsinki and was approved by the authorized institutional review boards (IRB) of the University of Cologne and the Charité-Universitätsmedizin Berlin.

Study Design

We performed a dual center retrospective observational imaging study. Imaging data including DW- and FLAIR-images were acquired from two stroke imaging databases: Group A, 1.5 T MR-imaging, University of Cologne, neurological imaging data base. Stroke patients available for the analysis were imaged consecutively between 2/2002 and 5/2004, in total 430 patients; Group B, 3 T MR-imaging, Charité-Universitätsmedizin Berlin, stroke imaging data base. Stroke patients available for the analysis were imaged consecutively between 3/2008 and 8/2010, in total 347 patients. Databases were screened and patients were included according to the following criteria: 1) clinically proven stroke, 2) confirmed symptom onset <12 h, 3) confirmed unilateral stroke lesion in DW-imaging, 4) available FLAIR imaging. Exclusion criteria were: 1) insufficient image quality, 2) incomplete clinical data, 3) punctate lesions and 4) brainstem infarctions. For a flow chart depicting the exclusions from the individual databases, see figure 1.

thumbnail
Figure 1. Database screening results and final study inclusion rate.

In both databases the exclusion criterion with the highest exclusion rate was a stroke-to-imaging-time higher than 12 hours. In the 1.5 T database, the number of patients, which had to be excluded due to insufficient image quality (mainly of FLAIR images), was much higher (17.1%) than in the 3 T database (2.0%).

https://doi.org/10.1371/journal.pone.0092295.g001

Magnetic Resonance Imaging Hardware

MR-imaging was performed at 1.5 T on a Philips Gyroscan Intera Master whole-body system (Philips Medical Systems, Best, The Netherlands). At 3 T, a Magnetom Tim Trio whole-body system (Siemens Healthcare, Erlangen, Germany) was used.

Magnetic Resonance Imaging Parameters

DW- and FLAIR imaging parameters were:

  • At 1.5 T:
    1. DWI: single shot SE-EPI (TE: 96 ms, TR: 3560 ms, flip angle: 90°, matrix: 256×256, FoV: 230×230, b:0 and b:1000, pixel size: 0.9×0.9 mm2, slice thickness: 6 mm, interslice gap: 0.6 mm)
    2. T2-weighted FLAIR: (TE: 100 ms, TR: 6000 ms, TI: 2000 ms, flip angle: 90°, matrix: 256×256, FoV: 220×220, pixel size 0.9×0.9 mm2, slice thickness: 6 mm, interslice gap: 0.6 mm)
  • At 3 T:
    1. DWI: single shot SE-EPI (TE: 93 ms, TR: 7600 ms, flip angle: 90°, matrix: 192×192, FoV: 230×230, b:0 and b:1000, pixel size 1.2×1.2 mm2, slice thickness: 2.5 mm)
    2. T2-weighted FLAIR: (TE: 100 ms, TR: 8000 ms, TI: 2370 ms, flip angle: 130°, matrix: 256×256, FoV: 220×220, pixel size 0.9×0.9 mm2, slice thickness: 5 mm, interslice gap: 0.5 mm).

Data Postprocessing and Image Analysis

Co-registration and post-processing of DW- and FLAIR images was performed with VINCI, Version 2.63 (Max-Planck-Institute for Neurological Research, Cologne, Germany) [9]. DWI lesion volumina were assessed using MRIcron (Chris Rorden, http://www.mccauslandcenter.sc.edu/mricro/). At 3 T, DW images were resized in the z-axis to match FLAIR images.

In DWI, the slice with the largest lesion extent was identified visually and used for the complete further analysis. DW-images were co-registered to FLAIR images and the absence or presence of FLAIR lesions was assessed by three raters blinded to clinical data and DW-images (experience in stroke imaging is indicated for each rater; Rater 1, VIM: 3 years, Rater 2, ME: 5 years, Rater 3, JS: 10 years). Prior to the rating, raters were encouraged to look for subtle intensity changes by adjusting contrast and brightness of the images and to compare intensities of potentially hyperintense regions with the healthy contralateral hemisphere. Such subtle intensity changes were also rated as a FLAIR hyperintensity. The area of the FLAIR lesion was individually delineated by each rater and copied on the DW-images. Then, 6 mm regions of interest (ROIs) were placed within the whole DWI-lesion. Each ROI was labeled according to its position in regard to the FLAIR-ROI. If it was located inside the FLAIR lesion, it was labeled FLAIR+ (positive), if it was located outside of the FLAIR lesion, it was labeled as FLAIR- (negative). If a FLAIR lesion was absent, all DWI ROIs of this patient were labeled as FLAIR-. For a graphical overview of the analysis see figure 2. ROI-values were normalized as a ratio: [100% x (mean ROI value/mean value of the unaffected hemisphere)], taken from a slice at the height of the lateral ventricles and above the putamen encompassing the corona radiata. In cases of cerebellar infarction, a ROI of the contralateral cerebellar hemisphere was used to normalize the ROI values. Above steps were performed equally for apparent diffusion coefficient (ADC)-maps.

thumbnail
Figure 2. Labeling of lesions in DWI imaging according to the FLAIR pattern.

A) In case of a FLAIR lesion encompassing the whole DWI-lesion, the hyperintensity was delineated as a region of interest (ROI) (A, second row). The FLAIR-ROI was then copied on the DWI and filled with 6 mm circular ROIs (A, third row)). These ROIs were classified as FLAIR+. B) In cases, in which the FLAIR-ROI did not completely match the DWI-lesion (B, second row), ROIs inside the FLAIR-ROI were classified as FLAIR+, and those outside as FLAIR- (B, third row). C) If no FLAIR lesion was identified (C, second row), the whole DWI lesion was filled with circular ROIs, which were classified as FLAIR- (C, third row). These steps were performed equally in ADC-maps. D), E) and F) show examples in analogy to the scheme, D) showing a patient, where the delineated FLAIR-ROI encompasses the whole DWI-lesion, E) depicting a patient, where the FLAIR-lesion only partially covers the DWI-lesion. Lastly, F) shows a patient, where all DWI-ROIs were labelled as FLAIR- in the absence of a visible FLAIR lesion.

https://doi.org/10.1371/journal.pone.0092295.g002

Statistical Analysis

Owing to skewed distribution of some variables, results are presented as median and interquartile range (IQR) if not indicated otherwise. Differences in clinical data between groups were assessed using the Mann-Whitney U rank sum test.

Agreement between raters for the identification of FLAIR-hyperintensities was analyzed using free-marginal kappa [10], [11]. Kappa values were evaluated as suggested by Landis and Koch [12].

ROI-analysis was performed on patient level and separately for the two centres. We used mean ROI-values per patient and per rater for positive ROIs and negative ROIs separately:

a) If all raters had some ROIs of a patient classified as having a positive FLAIR status, we used only the mean of the positive values and classified the patient as having a positive flair status.

b) If not all raters found positive ROIs for a patient, but all raters had ROIs classified as negative, we used only the mean of the negative ROI-values and classified the patient as negative. *

c) Only for some patients one of the raters classified all ROIs in another category as the other raters. For them we used the classification of the two corresponding raters and set the ROI value for the rater not corresponding to missing.

In the next step, the ability of a relative DWI-intensity threshold to predict the presence of corresponding FLAIR-hyperintensities was analyzed using an unadjusted receiver operating characteristics (ROC) curve analysis. The area under the curve (AUC) and the 95% confidence limits for the raters are reported. To get the optimal threshold, the Youden Index was used ( [13]. Sensitivity, specificity, and predictive values for the optimal thresholds are also reported.

To adjust for possible confounders, a multiple logistic regression model with the dependent variable “FLAIR-status” and independent variables “lesion volume”, “sex”, “thrombolysis” and “NIHSS” was used as a basic model (m0). In the additional model 1 (m1), “age” was added. Finally, different models were compared with regard to their ability to discriminate individuals in their FLAIR-status:

i) In model 2 (m2), we added “time-from-stroke-onset” to the m1 model.

ii) In model 3 (m3), we added the ROI-intensity for each rater separately to the m1 model. Paired sample statistical techniques were used for the comparison of two models. The method exploits the mathematical equivalence of the AUC to the Mann-Whitney U-statistic [14]. The ROC curves were calculated using SPSS Statistics 21, Release Version 21.0.0.0 (SPSS, Inc., 2012, Chicago, IL, www.spss.com). The comparisons of ROC curves and the linear mixed models were done using SAS software, Version 9.3 of the SAS System for Windows. (2010 SAS Institute Inc., Cary, NC, USA).

For analyzing the association between DWI-rSI and time-from-stroke-onset we calculated the mean DWI intensity over all ROIs and raters for every patient and used unadjusted and adjusted correlation analysis (Spearman's rank correlation) and a multiple linear regression analysis adjusted for “age”, “lesion volume” and “thrombolysis”. We calculated multiple linear regressions with (log-transformed) “mean DWI-value” as dependent variable and “age”, “thrombolysis” and (log-transformed) “lesion volume” as independent variables. Mean DWI intensities and lesion volume values were log-transformed to overcome the skewness in the distribution of the values. We analyzed the adjusted association between mean DWI-rSI and time-from-stroke-onset by analyzing the association of the residuals from the regression analysis with time-from-stroke-onset.

Above steps were performed equally for ADC-maps.

Results

In group A (1.5 T), 49 patients (16 females) and in group B (3 T) 48 patients (22 females) were included in the analysis. Median values for clinical data were (Group A/Group B): Time-from-stroke-onset (h) was 2.4/2.0; age (years) was 62/74; NIHSS (points) was 8/5 and the lesion volume (ml) was 22.9/6.8. The two groups differed significantly in age, NIHSS and lesion volume, but not in the time-from-stroke-onset. Five patients had cerebellar infarction (1 in group A and 4 in group B). Detailed patient data are shown in Table 1.

thumbnail
Table 1. Clinical data, imaging data and comparison of patient groups.

https://doi.org/10.1371/journal.pone.0092295.t001

Interrater agreement for the rating of FLAIR images as positive or negative for hyperintensities was substantial for both group A and group B with a kappa value of 0.62/0.69 (overall agreement was 81%/85%).

In the unadjusted ROC curve analysis, relative DWI-intensities performed well in discriminating hyperintensities in FLAIR imaging in both groups and for all 3 raters (results for groups and 3 raters; Group A: AUC 0.84, 0.91, 0.76 [mean: 0.84]; Group B: AUC 0.87, 0.86, 0.83 [mean: 0.85]). The Youden-Index identified comparable optimal relative DWI-intensity thresholds for both groups for all 3 raters (in %; Group A: 162, 158, 167, mean: 162; Group B: 163, 161, 159, mean: 161). In contrast, relative ADC-intensity values performed worse (group A: AUC 0.56, 0.55, 0.64; group B: AUC 0.74, 0.80, 0.77). Detailed data including sensitivity, specificity, positive and negative predictive value for the identified thresholds are listed in Table 2.

thumbnail
Table 2. Detailed results of the unadjusted ROC analysis for all 3 raters at 1.5 and 3 T.

https://doi.org/10.1371/journal.pone.0092295.t002

In the adjusted ROC-analysis for DWI, the basic model m0 including “lesion volume”, “sex”, “thrombolysis” and “NIHSS” had only a weak discrimination value for FLAIR status (AUC; Group A = 0.65; Group B = 0.59). The m1 model, where information on “age” was added, had a higher discrimination value (AUC; Group A = 0.73; Group B = 0.78). Adding “time-from-stroke-onset” to model m1 led to a further (significant) increase of the discrimination (model 2; AUC; Group A = 0.87; Group B = 0.91). Adding the ROI values for each rater as a variable in the m1-model, the increase was also significant in comparison with the m1 model and the discrimination value was comparable to the model 2 (model 3; mean AUC; Group A: 0.93; Group B: 0.93). Detailed data including p-values for the model comparison are listed in table 3. ROC-curves for model m0, m1, m2 and m3 for each rater are shown in figure 3.

thumbnail
Figure 3. Adjusted ROC curves for the detection of presence of FLAIR-lesions by a relative DWI- and ADC-threshold.

ROC-curves belonging to the detailed data presented in table 3 and 4 (please see legends of table 3 and 4 for further details). DWI-models for Group A (1.5 T) and B (3 T) (A,B) and ADC-models for Group A and B (C,D).

https://doi.org/10.1371/journal.pone.0092295.g003

thumbnail
Table 3. Detailed results of the adjusted ROC analysis for DWI and all 3 raters at 1.5 and 3 T.

https://doi.org/10.1371/journal.pone.0092295.t003

For ADC in contrast, adding of ROI values in the m3 model did not increase the discrimination value. On the contrary, model 3 performed even worse than model 2, which was based on “time-from-stroke-onset”. Detailed data including p-values for the model comparison are listed in table 4. ROC-curves for model m0, m1, m2 and m3 for each rater are shown in figure 3.

thumbnail
Table 4. Detailed results of the adjusted ROC analysis for ADC and all 3 raters at 1.5 and 3 T.

https://doi.org/10.1371/journal.pone.0092295.t004

In the unadjusted correlation analysis, a significant moderate to good correlation between mean relative DWI intensity and time-from-stroke-onset was found (Group A: rs = 0.54 (p<0.001); Group B: rs = 0.73 (p<0.011). The adjusted correlation confirmed a moderate to good correlation (Group A: rs = 0.45 (p<0.001;, Group B: rs = 0.69 (p<0.011) with a significant moderate fit in the linear regression analysis (see figure 4).

thumbnail
Figure 4. Adjusted linear regression analysis to evaluate the association of relative DWI-intensity and time-from-stroke-onset.

Adjusted linear regression analysis was performed to identify a possible association of relative (A) DWI-intensity and (B) ADC-intensity (y-axis,) and time-from-stroke-onset (x-axis) at 1.5 T (blue circles) and 3 T (green circle). At both field strengths, a significant association was found for DWI (A) with moderate adjusted Rsquare values (1.5 T: 0.28; 3 T: 0.44). Adjusted correlation (Spearman's rank correlation) was: 1.5 T = 0.45 (p<0.001), 3 T = 0.69 (p<0.001). In contrast, no association was found for ADC-maps (B) with adjusted Rsquare values near zero (1.5 T: 0.04; 3 T: 0.01) and weak to no adjusted correlation (1.5 T = −0.22, 3 T = 0.05). Plots are shown in logarithmic scale.

https://doi.org/10.1371/journal.pone.0092295.g004

For ADC in contrast, no correlation was found between ADC and time-frome-stroke-onset (unadjusted analysis: Group A rs = −0.25, Group B rs = 0.07; adjusted analysis: Group A rs = −0.22, Group B rs = 0.05) and no fit was present in the linear regression analysis (see figure 4).

Discussion

We report on the ability of a relative DWI-intensity threshold to discriminate with good accuracy between absence or presence of hyperintensities in corresponding FLAIR-images at both 1.5 and 3 T. The presence of FLAIR-hyperintensities was determined by 3 readers, who showed substantial interrater agreement. At both field strengths, similar DWI-intensity thresholds were identified. In addition, DWI intensity showed a significant association with time-from-stroke-onset.

To date, patients with unknown time from stroke onset are excluded from intravenous thrombolysis [15], [16]. As stroke incidence rates are higher in the morning hours compared to the rest of the day [17], patients in theory eligible for thrombolysis are in practice excluded from thrombolysis if time of stroke onset is unknown. Strategies to identify patients eligible for thrombolysis by MRI have been a major focus of interest in stroke research [18]. While ischemic lesions are visible in DWI as early as several minutes after stroke [19], lesions in T2-weighted FLAIR imaging show a later appearance, where the majority of patients displays FLAIR lesions only after several hours of stroke [2],[3],[4]. Thus, it is not surprising that the DWI-FLAIR mismatch allows for the allocation of patients to the current thrombolysis time window (i.e. <4.5 h after stroke) with a high specificity and a high positive predictive value [1]. The use of the DWI-FLAIR mismatch, however, mainly relies on the visual assessment of FLAIR hyperintensities. On one hand, visual qualitative assessment of FLAIR-images is challenging [6]. On the other hand, an automated analysis and the delineation of lesions in FLAIR imaging may be difficult owing to low contrast and partial volume effects [7]. This explains, why the use of quantitative relative FLAIR-intensity values has led to heterogenous results [3],[4],[20],[21]. We therefore hypothesized that DWI-maps might yield a surrogate of FLAIR imaging and we were able to show that a DWI based threshold predicts the presence of a FLAIR lesion with good performance across two different MR field strengths and three independent raters. In addition, mean DWI intensity of the lesion showed significant association with time-from-stroke-onset. These results strongly indicate two important points: First, DWI-intensity might exhibit a time-dependent increase after onset of ischemia in the acute phase. This specifies previous findings by Petkova et al. that described different DWI intensity values in patient samples stratified according to time-from-stroke-onset [4]. In that work, DWI-rSI was able to allocate patients to the thrombolysis time window <3 h (AUC: 0.75), but did not perform better than ADC-rSI (AUC: 0.74). It should be noted, however, that the authors projected the DWI lesion as a mask on ADC-maps and did not delineate the lesion in ADC-maps individually to derive rSI values. The reported threshold for DWI for the allocation of patients to the thrombolysis window was 19% using the formula “(DWI-lesion - contralateral value)/(DWI-lesion + contralateral value) × 100”. Recalculating our threshold of 160% rSI according to that formula leads to a similar threshold of 23% corroborating the pivotal association between FLAIR-appearance, DWI-rSI and time-from-stroke-onset.

Second, our results indicate that DWI-intensity is associated with tissue fate. The finding that a DWI-intensity threshold allows the prediction of FLAIR hyperintensity suggests that DWI-intensity follows a pathophysiologically driven time course similar to the manifestation of FLAIR lesions rather than simply increasing linearly with time. This is corroborated by previous results showing that a certain DWI-intensity threshold (118%) was able to predict permanent infarction as shown in an acute stroke sample imaged by MRI and comparative positron emission tomography [22]. This is important, as the DWI/FLAIR mismatch is limited by a high percentage of patients showing FLAIR positive lesion early after stroke [1]. Recently, it has been reported that this is even more pronounced at 3T, where a high percentage of patients showed FLAIR hyperintensities within the thrombolysis time window (44.5%) [23]. In this context, biomarkers are of interest, which are linked to tissue fate [3] and DWI intensity might serve this purpose.

It should be noted, that our analysis is an explorative and hypothesis-generating imaging study, in which only one representative slice per patient was analyzed. Based on our findings, it will be of major clinical interest to establish the relationship of DWI- intensity values and time-from-stroke-onset and tissue fate by analyzing the whole DWI lesion volume in future studies. Using such an approach, also the allocation to the 4.5 hours thrombolysis time window by DWI intensity can be investigated.

Our results rely in part on the assessment of FLAIR images as “positive” or “negative”. The raters were encouraged to look for even subtle changes by adjusting contrast and brightness and by comparing the putative lesion intensity with the intensity of the contralateral hemisphere. Following this predefined algorithm, a substantial agreement as measured by interrater kappa could be achieved. However, in this clinically relevant technique, overall agreement did not exceed 85%. This finding corroborates that the assessment of FLAIR imaging is prone to a subjective bias, especially in a clinical setting, where less rigid algorithms are applied to image rating. Despite of these findings, FLAIR imaging is an important clinical tool for the stratification of acute stroke patients and may be supported by other MRI parameters in clinical decision making in the future.

In contrast to DWI, ADC maps were not able to predict the time-dependent appearance of corresponding FLAIR lesions in acute stroke and showed no association with time-from-stroke-onset. Therefore, our results suggest that DWI intensity follows a time-dependent increase in intensity, while ADC-values do not. Several studies have described serial changes of ADC-values and DWI-intensity in acute human stroke, but focused on changes between the (hyper)acute, subacute or chronic stage. ADC values were reported to decrease in the acute stage and to increase again in the subacute and chronic stages [24][29]. DWI values were reported to increase between the acute and subacute stage [27] and to decrease in the chronic stage [30]. There is a substantial lack of knowledge regarding the evolution of DWI and ADC value changes within the acute phase of stroke. Based on these considerations, we see the need to characterize the evolution of DWI and FLAIR intensity values in humans by serial multiparametric MRI within the (hyper)acute phase of stroke.

DWI is a composite parameter of diffusion imaging and T2-imaging and thus a surrogate of very early restricted diffusion as well as following edema. ADC, on the other hand, is a quantification of diffusion alone. In our study, only DWI-intensity thresholds were able to predict the presence of FLAIR-lesions. A possible reason might be the additional information from T2 imaging present in DWI. In this respect, future studies evaluating the time-dependency of DWI from stroke-onset should shed light on the role of diffusion-weighting and T2-weighting as factors of the time-dependent intensity increase.

Importantly, we found that patient age was significantly associated with FLAIR-hyperintensities in both groups. Hence, patient age might be a confounder in the assessment of FLAIR-images for patient stratification in the acute stroke setting. This finding needs investigation to further define the predicition value of qualitative or quantitative FLAIR-imaging.

Our study has several limitations. First, a ROI based approach was chosen instead of a voxel-based analysis. A voxel-based analysis might be more accurate, but ROI-based approaches are less sensitive to spatial distortions occurring in echo planar imaging (EPI). Second, as evidenced by descriptive statistics (s. Table 1), lesion size was heterogenous in both groups, which might include a bias. Third, stroke volumes and stroke severity (based on NIHSS) were only moderate. It remains to be shown, whether our results are also applicable to patient samples with larger mean infarct volumes. Fourth, the patient groups were rather small. Our results should be validated in larger patient samples in future studies. Fifth, comparison of the results obtained at 1.5 T and 3 T is limited by a large time span between the measurements of the study groups. This led to a higher exclusion rate of images due to inferior image quality in the 1.5 T group as a results of advances in MRI techniques. This could have affected the study results. Sixth, owing to mild NIHSS values and small lesion volumes the cohort measured at 3T is not representative for patients eligible for thrombolysis. Our results must therefore be validated for a more diverse patient sample in future studies.

Summary

In conclusion, a relative DWI-intensity threshold predicted the presence of hyperintensities in FLAIR imaging at both 1.5 and 3 T with good accuracy in a retrospective sample. Moreover, DWI-intensity values were associated with time-from-stroke-onset. These findings suggest a time-dependent increase of DWI-intensity in the hyperacute phase of stroke. Future studies should investigate the value of DWI-intensity measurement as an easily accessable estimate of lesion-age.

Acknowledgments

We express our gratitude to the MRI staff of both the Department of Neuroradiology, University of Cologne, and the Academic Neuroradiology, Center for Stroke Research Berlin, Charite Universitätsmedizin Berlin.

Author Contributions

Conceived and designed the experiments: VIM IG WM ME JBF JS. Performed the experiments: VIM IG OZW JBF JS. Analyzed the data: VIM IG UG AS SZM FCvS KLS MAM. Wrote the paper: VIM IG OZW FCvS WM ME JBF JS UG.

References

  1. 1. Thomalla G, Cheng B, Ebinger M, Hao Q, Tourdias T, et al. (2011) DWI-FLAIR mismatch for the identification of patients with acute ischaemic stroke within 4·5 h of symptom onset (PRE-FLAIR): a multicentre observational study. Lancet Neurol 10: 978–986
  2. 2. Aoki J, Kimura K, Iguchi Y, Shibazaki K, Sakai K, et al. (2010) FLAIR can estimate the onset time in acute ischemic stroke patients. J Neurol Sci 293: 39–44
  3. 3. Ebinger M, Galinovic I, Rozanski M, Brunecker P, Endres M, et al. (2010) Fluid-Attenuated Inversion Recovery Evolution Within 12 Hours From Stroke Onset: A Reliable Tissue Clock? Stroke 41: 250–255
  4. 4. Petkova M, Rodrigo S, Lamy C, Oppenheim G, Touzé E, et al. (2010) MR imaging helps predict time from symptom onset in patients with acute stroke: implications for patients with unknown onset time. Radiology 257: 782–792
  5. 5. Galinovic I, Puig J, Neeb L, Guibernau Lisitano J, Kemmling A, et al.. (2014) Visual and Region of Interest-based Inter-Rater Agreement in the Assessment of the Diffusion-Weighted Imaging-Fluid-Attenuated Inversion Recovery Mismatch. Stroke, 2014 Feb 20 [epub ahead of print]
  6. 6. Ziegler A, Ebinger M, Fiebach JB, Audebert HJ, Leistner S (2011) Judgment of FLAIR signal change in DWI–FLAIR mismatch determination is a challenge to clinicians.J Neurol. Available: http://www.springerlink.com/content/623w4k071r6681g2/. Accessed 7 November 2011.
  7. 7. Khademi A, Venetsanopoulos A, Moody AR (2012) Robust White Matter Lesion Segmentation in FLAIR MRI. IEEE Trans Biomed Eng 59: 860–871
  8. 8. Lansberg MG, Lee J, Christensen S, Straka M, Silva DAD, et al. (2011) RAPID Automated Patient Selection for Reperfusion Therapy A Pooled Analysis of the Echoplanar Imaging Thrombolytic Evaluation Trial (EPITHET) and the Diffusion and Perfusion Imaging Evaluation for Understanding Stroke Evolution (DEFUSE) Study. Stroke 42: 1608–1614
  9. 9. Cízek J, Herholz K, Vollmar S, Schrader R, Klein J, et al. (2004) Fast and robust registration of PET and MR images of human brain. NeuroImage 22: 434–442
  10. 10. Randolph JJ (2005) Free-marginal multirater kappa: An alternative to Fleiss' fixed-marginal multirater kappa. Paper presented at the Joensuu University Learning and Instruction Symposium 2005, Joensuu, Finland. Available: http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.59.8776&rep=rep1&type=pdf.
  11. 11. Brennan RL, Prediger DJ (1981) Coefficient Kappa: Some Uses, Misuses, and Alternatives. Educ Psychol Meas 41: 687–699
  12. 12. Landis JR, Koch GG (1977) The measurement of observer agreement for categorical data. Biometrics 33: 159–174.
  13. 13. YOUDEN WJ (1950) Index for rating diagnostic tests. Cancer 3: 32–35.
  14. 14. DeLong ER, DeLong DM, Clarke-Pearson DL (1988) Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biometrics 44: 837–845.
  15. 15. Guidelines for Management of Ischaemic Stroke and Transient Ischaemic Attack 2008 (2008) Cerebrovasc Dis. 25: 457–507
  16. 16. Adams HP, Zoppo G del, Alberts MJ, Bhatt DL, Brass L, et al. (2007) Guidelines for the Early Management of Adults With Ischemic Stroke A Guideline From the American Heart Association/American Stroke Association Stroke Council, Clinical Cardiology Council, Cardiovascular Radiology and Intervention Council, and the Atherosclerotic Peripheral Vascular Disease and Quality of Care Outcomes in Research Interdisciplinary Working Groups: The American Academy of Neurology affirms the value of this guideline as an educational tool for neurologists. Stroke 38: 1655–1711
  17. 17. Elliott WJ (1998) Circadian Variation in the Timing of Stroke Onset A Meta-analysis. Stroke 29: 992–996
  18. 18. Sobesky J (2012) Refining the mismatch concept in acute stroke: lessons learned from PET and MRI. J Cereb Blood Flow Metab 32: 1416–1425
  19. 19. Schaefer PW, Grant PE, Gonzalez RG (2000) Diffusion-weighted MR Imaging of the Brain1. Radiology 217: 331–345.
  20. 20. Cheng B, Brinkmann M, Forkert ND, Treszl A, Ebinger M, et al.. (2012) Quantitative measurements of relative fluid-attenuated inversion recovery (FLAIR) signal intensities in acute stroke for the prediction of time from symptom onset. J Cereb Blood Flow Metab. Available: http://www.nature.com/jcbfm/journal/vaop/ncurrent/full/jcbfm2012129a.html. Accessed 28 November 2012.
  21. 21. Song SS, Latour LL, Ritter CH, Wu O, Tighiouart M, et al. (2012) A Pragmatic Approach Using Magnetic Resonance Imaging to Treat Ischemic Strokes of Unknown Onset Time in a Thrombolytic Trial. Stroke 43: 2331–2335
  22. 22. Heiss W-D, Sobesky J, Smekal U v, Kracht LW, Lehnhardt F-G, et al. (2004) Probability of Cortical Infarction Predicted by Flumazenil Binding and Diffusion-Weighted Imaging Signal Intensity: A Comparative Positron Emission Tomography/Magnetic Resonance Imaging Study in Early Ischemic Stroke. Stroke 35: 1892–1898
  23. 23. Emeriau S, Serre I, Toubas O, Pombourcq F, Oppenheim C, et al. (2013) Can Diffusion-Weighted Imaging–Fluid-Attenuated Inversion Recovery Mismatch (Positive Diffusion-Weighted Imaging/Negative Fluid-Attenuated Inversion Recovery) at 3 Tesla Identify Patients With Stroke at <4.5 Hours? Stroke 44: 1647–1651
  24. 24. Warach S, Gaa J, Siewert B, Wielopolski P, Edelman RR (1995) Acute human stroke studied by whole brain echo planar diffusion-weighted magnetic resonance imaging. Ann Neurol 37: 231–241
  25. 25. Schlaug G, Siewert B, Benfield A, Edelman RR, Warach S (1997) Time course of the apparent diffusion coefficient (ADC) abnormality in human stroke. Neurology 49: 113–119
  26. 26. Schwamm LH, Koroshetz WJ, Sorensen AG, Wang B, Copen WA, et al. (1998) Time Course of Lesion Development in Patients With Acute Stroke Serial Diffusion- and Hemodynamic-Weighted Magnetic Resonance Imaging. Stroke 29: 2268–2276
  27. 27. Lansberg MG, Thijs VN, O'Brien MW, Ali JO, Crespigny AJ de, et al. (2001) Evolution of Apparent Diffusion Coefficient, Diffusion-weighted, and T2-weighted Signal Intensity of Acute Stroke. Am J Neuroradiol 22: 637–644.
  28. 28. Fiebach JB, Jansen O, Schellinger PD, Heiland S, Hacke W, et al. (2002) Serial analysis of the apparent diffusion coefficient time course in human stroke. Neuroradiology 44: 294–298
  29. 29. Srivastava A, Mehrotra G, Bhargava S, Agarwal S, Tripathi R (2008) Studies on the time course of apparent diffusion coefficient and signal intensities on T2- and diffusion-weighted MR Imaging in acute cerebral ischemic stroke. J Med Phys 33: 162
  30. 30. Eastwood JD, Engelter ST, MacFall JF, Delong DM, Provenzale JM (2003) Quantitative Assessment of the Time Course of Infarct Signal Intensity on Diffusion-Weighted Images. AJNR Am J Neuroradiol 24: 680–687.