Quantification of ischemic brain edema after mechanical thrombectomy using dual-energy computed tomography in patients with ischemic stroke

Net water uptake (NWU) is a quantitative imaging biomarker used to assess cerebral edema resulting from ischemia via Computed Tomography (CT)-densitometry. It serves as a strong predictor of clinical outcome. Nevertheless, NWU measurements on follow-up CT scans after mechanical thrombectomy (MT) can be affected by contrast staining. To improve the accuracy of edema estimation, virtual non-contrast images (VNC-I) from dual-energy CT scans (DECT) were compared to conventional polychromatic CT images (CP-I) in this study. We examined NWU measurements derived from VNC-I and CP-I to assess their agreement and predictive value in clinical outcome. 88 consecutive patients who received DECT as follow-up after MT were included. NWU was quantified on CP-I (cNWU) and VNC-I (vNWU). The clinical endpoint was functional independence at discharge. cNWU and vNWU were highly correlated (r = 0.71, p < 0.0001). The median difference between cNWU and vNWU was 8.7% (IQR: 4.5–14.1%), associated with successful vessel recanalization (mTICI2b-3) (ß: 11.6%, 95% CI 2.9–23.0%, p = 0.04), and age (ß: 4.2%, 95% CI 1.3–7.0%, p = 0.005). The diagnostic accuracy to classify outcome between cNWU and vNWU was similar (AUC:0.78 versus 0.77). Although there was an 8.7% median difference, indicating potential edema underestimation on CP-I, it did not have short-term clinical implications.


Study design and study population
The local ethics committee (of University Medical Hospital Hamburg-Eppendorf, Hamburg; Number WF 04/13) approved the study and waived the requirement to obtain informed consent.The design and procedures of the study were carried out in accordance with the principles of the Declaration of Helsinki.This retrospective singlecenter study investigated whether dual-energy CT scans using VNC-I are superior for the assessment of NWU in acute stroke patients for clinical outcome prediction compared to traditional polychromatic CT scans.The local ethics committee approved the study and waived the requirement to obtain informed consent.We anonymized and retrospectively analyzed the data of 88 consecutive patients referred to our hospital between June 2020 and December 2021 with LVO of the anterior circulation who received MT.Inclusion criteria were defined as follows: (1) acute ischemic stroke in the anterior circulation due to an isolated occlusion of the intracranial carotid bifurcation or of the M1 or M2 segment of the middle cerebral artery confirmed on admission CT; (2) MT was performed according to the ESO/ESMINT-Guidelines 4 ; (3) dual-energy follow-up cranial CT scan after MT; (4) visually hypoattenuated and circumscribable infarct lesion on follow-up imaging.

Image acquisition and reconstruction
All CT scans were performed on a dual-energy dual-layer CT (IQon spectral CT, Philips Healthcare, USA) with the following imaging parameters: nonenhanced CT with 120 kVp, 230 mAs, and 5.0 mm slice reconstruction.As post-processing software IntelliSpace (Version 11.1, Philips Healthcare, USA) was used to generate CP-I and VNC-I from spectral based datasets.The CP-I were reconstructed using iterative model-based reconstructions (level 1, filter UB), whereas the VNC-I were generated using a spectral reconstruction mode (level 2) 14 .

Image analysis
The data was anonymized and the rater was blinded for all patient related data and clinical information.The rating was done by a neuroradiologist (PS) with > 6 years of experience in stroke image analysis.Results were confirmed by a second neuroradiologist (GB, 7 years of experience).Discrepancies were settled in a consensus reading.The rating was done using a commercially available software (Analyze 11.0, Biomedical Imaging Resource, Mayo Clinic, Rochester, MN).A standardized procedure to quantify the proportion of ischemic edema due to NWU was used as previously published and extended to account for CP-I and VNC-I comparison 7,23 .In summary, ischemic edema was quantified by placing a region of interest (ROI) on a sligle slice on visually evident edematous hypoattenuation on the conventional CT-scans for density measurements with knowledge of the core ischemic lesion from CT perfusion images.A symmetric ROI was mirrored within unaffected brain tissue on the contralateral hemisphere.Sulci were excluded manually.ROI histograms were sampled between 20 and 80 Hounsfield units to exclude voxels belonging to cerebrospinal fluid or calcifications.Both density measurements were used to calculate the proportion of edema (cNWU) as the difference between density of the infarcted area to the unaffected hemisphere.The ROIs from the conventional images were saved and reused for the same person's VNC-I (for vNWU) to ensure that both, CP-I and VNC-I density measurements were performed at identical

Statistical analysis
The difference between cNWU and vNWU was defined es edema underestimation with higher values indicating an increase in edema underestimation.The median NWU difference between cNWU and vNWU was used to divide the study cohort into two groups of patients (i.e.low edema underestimation versus high edema underestimation).Non-parametrical Kruskal-Wallis/Mann Whitney U-Test was used for continuous variables and were reported as median and interquartile range (IQR).Pearson's chi-squared test was used for categorical/binary variables and are reported as counts and percentages.cNWU and vNWU were compared using paired sample t-tests with means and 95% confidence intervals.Pairwise Correlation between cNWU and vNWU was tested.Multivariable linear regression analysis was performed with the absolute difference of cNWU and vNWU as dependent variable implementing variables with established impact on edema formation (time from onset, age, ASPECTS, and status of vessel recanalization).Finally, the potential clinical relevance of edema underestimation was tested by comparing two multivariable logistic regression models with mRS 0-2 at discharge as dependent variable.Variables were selected based on established associations to outcome (Age, NIHSS, ASPECTS, and degree of recanalization).For Model 1, cNWU was added, and for Model 2, vNWU was implemented.The area under the curve (AUC) for both models were calculated.Analyses were performed using Stata

Patient Characteristics and stratification by low and high edema underestimation
88 consecutive patients were assessed for eligibility, of which 70 patients met the inclusion criteria.18 patients were excluded due to an absence of any infarct demarcation on follow up imaging (14 patients), parenchymal hemorrhage type 2 according to Heidelberg Bleeding Classification 24 (2), strong artifacts from extracranial metal (1), and reinfarction with thrombectomy before follow-up CT scan (1).Among the included patients, mean Age was 76.1 years (range: 47-93) with 54.3% being female.The median NIHSS score on admission was 13 points (IQR: 8-18) and the baseline ASPECTS was 8 (IQR 7-9).Successful recanalization defined as mTICI 2b, 2c or 3 was achieved in 58 cases (83%) of patients.
Stratified by low (< 8.7%) and high (≥ 8.7%) edema underestimation, no significant differences between groups were observed regarding clinical parameters (NIHSS and mRS) on admission or discharge.No difference on admission ASPECTS was seen.Time between symptom onset to flow restoration differed significantly between the two study groups with 289 min in patients with high edema underestimation compared to 372 min in patients with low edema underestimation (p = 0.019).Successful recanalization was achieved in 26 (74.3%) patients with low edema underestimation and in 32 (91.4%) patients with high edema underestimation showing a trend but no significant difference (p = 0.057).Table 1 shows baseline, treatment, and outcome characteristics for the two study groups.

Independent variables associated with edema underestimation
Multivariable linear regression analysis was used to test the association between baseline and treatment variables with the difference between cNWU and vNWU as a measure of edema underestimation (Table 2).Successful recanalization (ß: 11.8%, 95% CI 0.27-23.4,p = 0.04), and age (ß: 0.04%, 95% CI 0.01-0.07,p = 0.005) were significantly associated with edema underestimation.In multivariable linear regression, time from symptom onset to flow restoration did not show significant difference between groups.The relationship between age and successful recanalization as well edema underestimation is shown in Fig. 3.

Functional outcome depending on NWU and the degree of underestimation
The rate of good functional outcome (mRS 0-2 at discharge) was 40% in patients with high edema underestimation compared to 45% good outcome in patients with low edema underestimation (Fig. 4).The median NIHSS at discharge was 5 in both groups (IQR 1.00-16.50for patients with low edema underestimation and IQR 2.00-14.00for patients with high edema underestimation), which was not significantly different p = 0.92.Similarly, the median mRS at discharge was also not different between groups (p = 0.87) (Table 1).

Discussion
The main purpose of this study was to investigate whether the use of conventional CT scans to measure ischemic edema may result in underestimation due to the presence of contrast staining or hemorrhage.In this study, the viability of VNC-I as a superior alternative to CP-I for NWU measurements is explored to mitigate the negative effects of contrast staining on NWU and reduce the underestimation of ischemic edema.While previous research has highlighted the clinical value of NWU as an imaging biomarker 10,25,26 , it is uncertain whether VNC-I has any benefits over CP-I for quantifying ischemic edema 27,28 .A pilot study investigating the use of VNC-I demonstrated that DECT may be advantageous in 20% of patients compared to conventional CT 14 , which is consistent with another study that reported the presence of hyperattenuating lesions after MT in approximately 20% of patients on 24-h follow-up CT scans 13 .It is important to acknowledge that the occurrence of contrast staining is believed to arise from a complex of multiple factors, leading to increased blood-brain-barrier permeability in infarcted brain parenchyma 29 .While Hemorrhage remains in the extravascular space for days and should be hyperdense beyond 48 h, contrast staining is time sensitive 30 and is likely be absorbed within 48 h 12 .Prolonged contrast staining on follow-up DECT CT 48-72 h after thrombectomy seems to be associated with poor outcome, independently of recanalization status 31 .Given that vNWU is measured using VNC-I, it is likely that the differences between cNWU and vNWU are attributable to contrast staining.For the current study, no difference in time (flow restoration to follow-up CT) was seen between groups of low and high edema underestimation.
The main findings of this study were that the median underestimation of edema (difference between cNWU and vNWU) was 8.7%, and that age and successful recanalization influenced edema underestimation.Overall, correlation between cNWU and vNWU was strong.Median cNWU was 19.7% which is consistent with Table 2. Multivariable linear regression analysis to determine independent variables associated with edema underestimation.n = 70 patients were included in multivariable regression analysis.Age, Blood Glucose, Admission ASPECTS and Time from symptom onset to flow restoration was calculated per 1 respective unit.Successful recanalization is defined as mTICI 2b or higher.ASPECTS Alberta Stroke Programme Early CT Score.* Indicates significant difference between groups.previously published data stating NWU values between 17.6% for successfully recanalized patients (mTICI 2b or better) and 24.9% for persistent LVO (mTICI 0-1) on follow-up CT scans 8 .In multivariate analysis, two parameters remained significantly associated with edema underestimation: Successful recanalization and age.Recent studies showed that, successful recanalization is associated with lower ischemic edema formation 8,32 and decreased incidence of symptomatic intracranial hemorrhage (SICH) 33 .While SICH is easily identifiable on CP-I and should always be excluded from the chosen ROI for NWU measurements, other hyperdense lesions such as contrast staining and/or hemorrhagic transformation can be obscured on conventional CT techniques.Contrast staining tends to have a lower density than hemorrhagic transformation (HU < 50 vs.HU > 90) on immediate postprocedural CT scans 34 , making it especially difficult to differentiate from vital brain parenchyma.In the current study, successful recanalization was associated with increased edema underestimation, suggestive of reperfusion injury with contrast staining.Yet, the short-term clinical outcome was not different between low and high edema underestimation, questioning this conclusion.An extensive review addressing contrast extravasation after cerebral revascularization, concluded a lack of consensus on how to best interpret these radiographic findings with some papers describing contrast staining as a predictor for poor clinical outcome while others did not find any correlation to clinical outcome 29 .

Beta coefficient 95% confidence interval p value
In elderly patients, the vulnerability of the vessels is increased due to lower vessel compliance 35 , resulting in more frequent hemorrhagic transformation and contrast staining 36 .Both increase the Hounsfield Units in CP-I.Especially contrast staining leads to a stronger difference to VNC-I, and thus to an increased underestimation of edema.Moreover, the combination of age-related small vessel disease and reduction of brain volume may directly impair the accuracy of densitometric assessments due to an increased number of voxels that do not represent brain parenchyma.
In the past, it has been observed that NWU measured in follow-up CT is a predictor of outcome based on quantification on CP-I.In this study, no significant difference was seen in the association between outcome and degree of edema underestimation, further strengthening the impression that both cNWU and vNWU are reasonable imaging techniques for NWU quantification.One reason for the high correlation lies within the ROIs chosen for NWU measurements.Wherever possible, the ROI was placed so that visually apparent hemorrhage or staining was not within the area used for measurement (as depicted in Fig. 1).This is an established method and important to reduce the effect of contrast staining and hemorrhage on NWU quantification 10,26,37 .This is a retrospective observational, single-center study with several limitations.Firstly, a relatively small number of patients were included owing to inclusion requirements such as dual-energy CT as follow-up examination, which was only recently adopted for follow-up imaging at our clinic.Moreover, the quantification of NWU using DECT imaging has not been fully established, and caution should be exercised when interpreting the results for treatment analysis.Studies have demonstrated that VNC-I and iodine overlay maps can effectively distinguish between contrast agents and blood products.However, it is important to note that VNC-I cannot be directly compared to standard non-enhanced CT scans due to differences in tissue contrast.For example, www.nature.com/scientificreports/ the differentiation between gray and white matter may be compromised and the impact of a secondary (nonstandardized) calculation of the virtual unenhanced non-contrast image on the quantification of NWU is not yet fully understood.Secondly, both patients with good prognosis (no identifiable demarcation) and patients with poor prognosis (PH2) were excluded.This might lead to outcome biases and impair generalization of the results.Thirdly, the routine use of DECT as a follow-up in acute ischemic stroke patients could increase patient radiation exposure, raising concerns in daily practice.Therefore, its application should always be carefully considered in light of the benefits.Lastly, short-term clinical outcome was evaluated as mRS at discharge whereas mRS at 90 days is often used in literature but was not obtained for this study.However, mRS at discharge strongly predicts 90 days mRS 38 .
Besides the technical implications of improving NWU quantification, the impact on outcome prediction may be of particular relevance.In this study, NWU measurements in conventional CP-I and VNC-I were highly correlated.A median edema underestimation of 8.7% was measured for cNWU, but it had no short-term clinical implication as edema underestimation was not associated with outcome at discharge using both NIHSS and mRS scores.VNC-I may have its use in patients with low edema formation as well as in patients with diffuse hyperattenuations, which might be hard to identify on CP-I.Future research is required to investigate the association of www.nature.com/scientificreports/edema underestimation with long-term functional outcome and to assess whether vNWU quantification may be used as a tool to improve the comparability of infarct imaging endpoints in ischemic stroke studies.

Conclusion
NWU measurements in conventional follow-up CT was highly correlated to NWU measurements in DECT-based VNC-I maps.A median edema underestimation of 8.7% was observed but without short-term clinical implications.The results confirm that CP-I as well as VNC-I can be used to quantify NWU reliably.Further research is needed for long-term clinical implications.

Figure 1 .
Figure 1.Examples for region of interest (ROI) placements in polychromatic (1st and 3rd images from left) as well as corresponding virtual non-contrast (2nd and 4th images from left) images for net water uptake measurements.The ischemic infarction (Red) was compared to the corresponding unaffected contralateral region (green).Top row: Follow-up dual-energy CT within 16 h in an 84-year-old female with right middle cerebral artery (MCA) occlusion (initial NHISS 20; ASPECTS 7) after thrombectomy (mTICI 0).Example to demonstrate how visually apparent hyperattenuations in conventional polychromatic images (CP-I) were spared during ROI placement.Edema underestimation was 12.8%.Bottom row: Follow-up dual-layer dual-energy CT within 16 h in a 69-year-old man with right MCA occlusion (initial NHISS 17; ASPECTS 4) after thrombectomy (mTICI 2b).Example of ischemic ROI placement if no visually apparent hyperattenuations were seen on CP-I.Edema underestimation was 4.5%.

Figure 3 .
Figure 3. Relationship of successful recanalization and age on the degree of edema underestimation in %.Elderly patients show an increase in edema underestimation, irrespectively of successful or unsuccessful recanalization.Overall, successful recanalization was associated with higher edema underestimation.mTICI modified Thrombolysis in Cerebral Infarction.

Figure 4 .
Figure 4. Distribution of modified Rankin Scale at admission and at discharge stratified by low (< 8.7%, LowEU) and high (≥ 8.7%, HighEU) edema underestimation.Shown is the distribution of modified Rankin Scale (mRS) which range from 0 (no symptoms) to 6 (death).The rate of good functional outcome (mRS 0-2 at discharge) was 40% for patients with high edema underestimation compared to 45% good outcome for patients with low edema underestimation (p = 0.54).

Table 1 .
Baseline, treatment, and outcome characteristics stratified by < 8.7% and ≥ 8.7% edema quantification difference.Patient characteristic of subgroups.Successful recanalisation is defined as mTICI 2b or higher.Time data is given in minutes.N number of patients, SD standard deviation, IQR interquartile range, NIHSS National Institutes of Health Stroke Scale, mRS modified Rankin Scale, ASPECTS Alberta Stroke Programme Early CT Score.P-value < 0.05 indicates clinical significance.