Multimodal validation of focal enhancement in intracranial aneurysms as a surrogate marker for aneurysm instability

Purpose Circumferential enhancement on MR vessel wall imaging has been proposed as a biomarker of a higher risk of rupture in intracranial aneurysms. Focal enhancement is frequently encountered in unruptured aneurysms, but its implication for risk stratification and patient management remains unclear. This study investigates the association of focal wall enhancement with hemodynamic and morphological risk factors and histologic markers of wall inflammation and degeneration. Methods Patients with an unruptured middle cerebral artery aneurysm who underwent 3D rotational angiography and 3T MR vessel wall imaging showing focal wall enhancement were included. Hemodynamic parameters were calculated based on flow simulations and compared between enhanced regions and the entire aneurysm surface. Morphological parameters were semiautomatically extracted and quantitatively associated with wall enhancement. Histological analysis included detection of vasa vasorum, CD34, and myeloperoxidase staining in a subset of patients. Results Twenty-two aneurysms were analyzed. Enhanced regions were significantly associated with lower AWSS, lower maxOSI, and increased LSA. In multivariate analysis, higher ellipticity index was an independent predictor of wall enhancement. Histologic signs of inflammation and degeneration and higher PHASES score were significantly associated with focal enhancement. Conclusion Focal wall enhancement is colocalized with hemodynamic factors that have been related to a higher rupture risk. It is correlated with morphological factors linked to rupture risk, higher PHASES score, and histologic markers of wall destabilization. The results support the hypothesis that focal enhancement could serve as a surrogate marker for aneurysm instability.


Introduction
Intracranial aneurysms have a prevalence of up to 3% and generally harbor a low risk of rupture [1,2]. Still, subarachnoid hemorrhage following rupture of an intradural aneurysm is associated with considerable mortality and morbidity [3]. Therefore, risk stratification of patients diagnosed with an unruptured intradural aneurysm is crucial, but optimal management remains controversial. Recently, published data suggest an association of wall enhancement in intracranial saccular aneurysms on MR vessel wall imaging with a higher risk for rupture [4][5][6][7]. Specifically, thick circumferential wall enhancement has been linked to a higher risk of rupture [4]. Moreover, wall enhancement probably identifies ruptured aneurysms in patients with subarachnoid hemorrhage and multiple aneurysms [8], and in AVMs, respectively [9]. In a prior study, inflammatory and degenerative wall changes were Naomi Larsen and Charlotte Flüh contributed equally to this work. detected in a histologic analysis of unruptured aneurysms showing wall enhancement [10].
Focal enhancement in unruptured aneurysms is frequently encountered on MR vessel wall imaging; still, its implication for patient risk stratification and management remains unclear. Hence, the goal of this study was to determine the correlation of focal enhancement with previously described risk factors in order to assess its suitability as a surrogate marker for a higher risk of rupture. To achieve this goal, we chose a multimodal approach: We investigated whether focal wall enhancement is colocalized with specific hemodynamic conditions related to a higher risk of rupture. Furthermore, the association with morphological parameters linked to rupture risk was assessed. Histologic data were analyzed to associate pathologic processes encountered during aneurysm evolution and wall remodeling with wall enhancement.

Patients
Ethics approval for this study was obtained from the local ethics committee. Written informed consent was obtained. We selected all patients from our institutional MR vessel wall imaging database with an unruptured middle cerebral artery aneurysm, who received high-resolution 3T MR vessel wall imaging showing focal wall enhancement and threedimensional rotational angiography (3D-RA) from 7/2016 to 3/2019. Only cases with MRI and 3D-RA images of sufficient quality were included. Demographic data as age, gender, intake of acetylsalicylic acid and nicotine, hypertension, diabetes, and PHASES score were extracted from patient records. Histologic data were collected from a subset of patients who underwent microsurgical clipping of the aneurysm. In our institution, the neurosurgeon and interventional neuroradiologist determine in an interdisciplinary consensus-based decision-making process whether to recommend treatment of the aneurysm and the preferred method.
Contrast enhancement of the aneurysm wall was visually assessed on multiplanar reformatted postcontrast 3D T1 black blood images.
Focal wall enhancement was defined as a hyperintense signal adjacent to the black blood lumen or wall-adherent thrombus with at least one clearly discernible unenhanced wall segment of any extent. The presence of a double-layer appearance indicative of thrombus, defined as a non-enhancing, walladjacent structure bordered by an enhancing layer on both the luminal and the wall-adjacent surface, was documented.
Two neuroradiologists (N.L., F.W.), each with at least 9 years of experience in neurovascular imaging, who were blinded for clinical and hemodynamic data, independently assessed the aneurysm wall. Disagreement between the readers was solved by consensus.

Digital subtraction angiography
Digital subtraction angiography (DSA) was performed under local anesthesia on a biplane flat panel DSA unit (Allura Xper FD 20/10, Philips, Best, The Netherlands) via a femoral access. 3D-RA was performed with a selective contrast injection of 20 ml Imeron 300 (Iomeprol, Bracco Imaging, Milan, Italy) in the ipsilateral internal carotid artery at a flow rate of 2 ml/s. The dataset was transferred to a dedicated workstation (XtraVision, Philips, Best, The Netherlands) and reconstructed with a voxel size of 0.27 × 0.27 × 0.27 mm 3 .

Vessel segmentation and image co-registration
Segmentation of vessel surfaces from 3D-RA datasets and coregistration to MRI was carried out using MeVisLab 2.7 (MeVis Medical Solutions AG, Bremen, Germany). A threshold-based segmentation approach was applied to the 3D-RA images, followed by marching cubes to obtain a triangulated surface mesh. Enhanced regions were manually delineated in black blood MR images using the contour segmentation objects library of MeVisLab. The enhancement mask was converted into a triangle surface mesh as well. Finally, coregistration utilizing a rigid multi-resolution registration [25] with a Quasi-Newton optimizer within MeVisLab yielded transformation matrices for each patient's dataset. Applying the transformation matrices to the manually segmented enhancement regions allowed for the combination with the vessel surface (see Fig. 1). The percentage of the aneurysm surface area covered by enhanced regions was determined for each aneurysm (enhancement area EA) and utilized as a quantitative measure for the extent of enhancement.

Hemodynamic simulations
Based on the segmented aneurysm surfaces, hemodynamic simulations derived from 3D-RA datasets were carried out in all cases using STAR-CCM+ 13.04 (Siemens PLM Software Inc., Plano, TX, USA). For spatial discretization a base size of Δx = 0.08 mm was chosen leading to 1.2 to 3.1 million polyhedral and prism elements in total depending on the considered vasculature. At each inlet, a time-dependent flow waveform was applied obtained from PC-MRI measurements in a healthy volunteer [26]. Flow-splitting at all outlet cross-sections was defined according to the recent model by Chnafa et al. [27]. Blood was treated as an incompressible (ρ = 1055 kg/m 3 ) and Newtonian (η = 4 mPa s) fluid, while the flow was assumed to be laminar. Three cardiac cycles were simulated, while only the last one was considered for analysis.
Quantification of shear-related parameters was conducted for enhanced and non-enhanced wall segments separately, including time-averaged wall shear stress (AWSS), normalized wall shear stress (nWSS), maximum oscillatory shear index (maxOSI), and low shear area (LSA). Mean neck inflow rate (meanNIR), inflow concentration index (ICI), and shear concentration index (SCI) were determined for entire aneurysms [28].

Semi-automatic extraction of morphological parameters
In each aneurysm, a semiautomatic neck curve reconstruction was performed. First, four points in the neck region of the aneurysm were determined following a manual one-click selection of the aneurysm. Based on these points, a neck curve was created and the ostium plane derived. For a detailed description of the computational process, see Saalfeld et al. [29]. Utilizing the neck curve and ostium plane, the following morphological parameters were automatically extracted: aneurysm surface area (A), aneurysm volume (V), ostium area1 (OA1), ostium area2 (OA2), maximum diameter (Dmax), maximum height (Hmax), maximum width perpendicular to Hmax (Wmax), maximum height perpendicular to the ostium plane (Hortho), maximum width parallel to the ostium plane (Wortho), maximum neck curve diameter (Nmax), average neck curve diameter (Navg), aspect ratio 1 (AR1), aspect ratio 2 (AR2), volume of the convex hull (Vch), surface area of the convex hull (Ach), ellipticity index (EI), non-sphericity index (NSI), undulation index (UI), and aneurysm tilt angle γ [29].

Aneurysm clipping and histologic analysis
Microsurgical clipping of the aneurysms was performed via a pterional-transsylvian approach.
Preoperative, postoperative, and intraoperative management of the patients was performed according to institutional standards, including preoperative and postoperative DSA of the intracranial vessels, and, in most cases, intraoperative angiography with indocyanine green. After clipping of the aneurysm, the aneurysm sac was removed.

Statistical analysis
Statistical analysis was performed in RStudio (RStudio v1.2.1335 with R v3.5.1, RStudio Inc., Boston, USA). Normally distributed variables are presented as mean (range), non-normally distributed continuous variables as median (range) or median (standard deviation), and categorical variables in the form of percentage distributions. A Shapiro-Wilk test was employed to determine, if a continuous variable was normally distributed. Non-normally distributed data were nonlinearly transformed before they were included in univariate analysis. Hemodynamic parameters were analyzed using the Wilcoxon signed-rank test for comparison of enhanced regions and the entire aneurysm surface, and the Mann-Whitney U test for comparison of the enhancement area EA in aneurysms with and without histologic markers of inflammation and wall degeneration. Pearson's correlation test was used to assess the relationship of hemodynamic and morphological parameters with the enhancement area EA. Only parameters that were significantly associated with EA in univariate analysis were entered into a linear regression model in order to restrict the number of variables included in multivariate analysis. Furthermore, since some of the morphological variables have been reported to show a strong degree of correlation [29], only the variable with the lowest P value was selected for multivariate analysis. A P < 0.05 was considered significant in univariate and multivariate analyses.

Results
Twenty-one patients with 22 aneurysms were included; 71% were located on the right side. Clinical data were available for all included patients. The interrater agreement concerning focal enhancement was rated as substantial (Cohen's κ = 0.742). Demographic information of the patients is summarized in Table 1. The PHASES score was positively correlated with EA (P = 0.0275). An association of gender, age, hypertension, diabetes, ASA and statin intake, smoking, and aneurysm location with EA could not be demonstrated.
Microsurgical clipping with subsequent histologic analysis of the resected aneurysm sac was conducted in 9 patients. The surgical indication was based on the patient's preference in 5 cases, and on aneurysm morphology (wide-necked aneurysm) in 4. Features of all aneurysms alongside histologic data are listed in Table 2.

Wall enhancement and blood flow evaluation
Flow simulation and calculation of hemodynamic variables were successful in all patients. A qualitative observation of the blood flow revealed that flow structures perpendicular to the main flow direction occurred in the vicinity of the enhancement zones. In cases with an increased aspect ratio, it was observed that these flows move sideward through the aneurysm sac instead of following the global flow direction. Additionally, it is important to highlight that enhancement appeared predominantly at the sides of the wall and was only present at the aneurysm dome in one case. Thus, the entering flow-jet and the corresponding impingement zone (which causes increased shear stresses on the luminal surface) were not associated with the enhanced wall segments.
To further quantify the flow behavior along the aneurysmal lumen comparing enhanced and non-enhanced regions, three shear-related parameters were assessed.
There was no significant correlation of inflow-related hemodynamic parameters or SCI with EA. There was a trend towards higher LSA in aneurysms with a larger extent of enhancement (P = 0.08222). For a summary of the analysis of hemodynamic and morphologic parameters, see Table 3.

Wall enhancement and morphologic parameters
OA2, Dmax, Hmax, Wortho, Vch, Ach, EI, NSI, and UI were significantly correlated with EA in univariate analysis (Table 3). EI showed the lowest P value and was entered alongside Dmax into multiple regression analysis, the remaining variables were not considered due to collinearity as described before [29]. EI remained significantly associated with EA (P = 0.0395).

Histologic features
Four aneurysms stained positive for MPO-expressing cells. In three aneurysms, vasa vasorum were present. CD34-positive neovessel endothelium could not be detected in any of the In this subset of patients, the presence of one or two of these histologic markers was significantly associated with a larger extent of wall enhancement (P = 0.03175). Size (as manually measured on 2D angiographic images), enhancement area, and histologic features of all aneurysm. Histologic specimens were available from aneurysms 1-9 Wall-adherent thrombus was detected on MR vessel wall imaging in two aneurysms, and intraluminal wall adherent thrombus was confirmed by histologic analysis. Enhancement was visible both on the luminal surface of the thrombus and the adjacent wall, leading to a double-layer appearance (Fig. 2c).
Extracellular lipid deposition was detected in two aneurysms. Hyalinosis, fibrosis, and calcification were common features in most aneurysms.

Focal enhancement, hemodynamic parameters, and morphology
In this study, focal enhancement was colocalized with low AWSS, low maxOSI, and high LSA in agreement with results recently published by Xiao et al. and Khan et al. [13,31], confirming the assumption that enhanced regions possibly indicate the presence and localization of low-flow conditions in an aneurysm.
This observation adds to the assumption that the extent of focal enhancement depends on the prevalence of low-flow conditions. Furthermore, it seems likely that there is a continuous spectrum of the extent of wall enhancement correlated to a hemodynamic environment promoting wall destabilization that has to be considered when assessing rupture risk, rather than an either-or presumption that is implied when dichotomizing unruptured aneurysms into either enhancing (and thus at higher rupture risk) or non-enhancing. There was no significant correlation of AWSS, maxOSI, and LSA with EA when computed for whole aneurysms. Still, there was a trend towards higher LSA in aneurysms with a larger extent of enhancement (P = 0.08222). This could be explained by the fact that the percentage of the aneurysm surface that was covered by enhancement was relatively low in the aneurysms included in this study (median 27%; see Table 3). Future longitudinal studies focusing on quantification of aneurysm wall enhancement on MR vessel wall imaging and correlated clinical parameters are warranted.
Numerous studies have described morphologic characteristics associated with a ruptured state, such as larger aneurysm size and greater aspect ratio [15,17,22,32,33]. Yet, a retrospective longitudinal analysis conducted by Leemans et al. yielded no significant differences in hemodynamic and morphological parameters at baseline imaging in aneurysms that had grown versus stable aneurysms [23]. Morphological and hemodynamic analysis in cases where longitudinal data are not available might therefore not be able to sufficiently discriminate between stable and unstable aneurysms. In our patient cohort, aneurysm enhancement was positively correlated with hemodynamic, morphologic, and histologic markers that have been linked to an unstable state, which supports the assumption that wall enhancement might reflect the complex interactions and multifactorial pathologic processes in aneurysm evolution leading to wall destabilization, while a direct causal relation to either of these factors cannot readily be assumed.

Low-flow conditions, thrombi, and histologic findings
In a comprehensive study linking macroscopic wall appearances to specific hemodynamic conditions, Cebral et al. found low AWSS to be associated with macroscopically hyperplastic and atherosclerotic wall segments [34]. Concordantly, we observed histologic markers of inflammatory and remodeling processes in aneurysms with a greater extent of focal wall enhancement, suggesting that these pathologic conditions are more abundant in conjunction with low-flow conditions. Moreover, rupture sites of intracranial aneurysms have been associated with low wall shear stress [14,[34][35][36] and focal enhancement [11,12]. Still, a direct association of lowflow areas and possible future rupture site cannot directly be deduced from the results of the present study, because local matching of histologic findings and enhancements areas was not conducted.
Recently, Sato et al. published their results from a 7T vessel wall MRI study describing the double-layer appearance in thrombosed aneurysms and linking them to histologic signs of inflammation, in accordance with our results [37].
Matsushige et al. described focal enhancement as possibly attributable to the presence of thrombotic material at rupture sites [11]. Another study attempted to categorize aneurysm wall types and found thrombosis as a feature of categories with a higher proportion of ruptured aneurysms [38]. These observations add to the assumption that slow blood flow and thrombus formation might play a substantial role in wall destabilization and can be visualized with MR vessel wall imaging.
Still, the observation that focal enhancement, albeit to a lesser extent, is also encountered in aneurysms lacking histologic signs of wall inflammation, possibly reflects the proposed theory that there are different pathways of wall destabilization in intracranial aneurysms [39]. Substantial wall enhancement might primarily reflect the inflammatory pathway of wall remodeling, while the sensitivity for the mural-cellmediated pathway could be diminished.
Cornelissen et al. assumed that wall enhancement is not causally linked with inflammatory processes but rather attributable to method-inherent failure of blood signal suppression in low-flow conditions [40,41]. The results of our study suggest that inflammatory changes probably add to the extent of wall enhancement, and that even without convincing proof of a causal link, wall enhancement indicates the presence of a hemodynamic environment promoting inflammatory and degenerative processes and might therefore still serve as a surrogate marker for aneurysm wall instability.
The extent of enhancement seems to depend on the sequence employed in the vessel wall MR protocol, with the T1 TSE black blood sequence like the one used in our study possibly being more sensitive for the detection of slow blood flow compared with motion-sensitized driven equilibrium (MSDE)-prepared sequences [42].
This study has several limitations. The retrospective study design might imply a selection bias towards patients who underwent DSA and were potentially pre-selected for therapy due to their risk profiles. Additionally, histologic analysis was exclusively performed in patients who underwent clipping, possibly inducing bias towards aneurysms with morphologies unfavorable for endovascular therapy. Histological data were only available in nine aneurysms. Microsurgical clipping provides only specimens of the aneurysm dome; neck segments under the clip cannot be retrieved. The small sample size precluded subgroup analysis of histologic features with adequate power. We chose to restrict our analysis to middle cerebral artery aneurysms in an effort to exclude a possible heterogeneity in aneurysm phenotype caused by different pathogenetic mechanisms potentially underlying the evolution and progression of aneurysms in different locations [43]. Additionally, an overestimation of the neck size using 3D imaging techniques can occur compared with 2D images, but careful segmentation using advanced algorithms was applied [44,45]. Finally, no patient-specific inflow conditions were available for the blood flow simulations. However, a comparison between enhanced regions and the whole aneurysm remains feasible.

Conclusion
Focal wall enhancement in unruptured intracranial aneurysms is colocalized with hemodynamic factors that have been related to a higher risk of rupture. It is correlated with morphological factors linked to rupture risk, a higher PHASES score, and histologic markers of wall destabilization. The results support the hypothesis that focal enhancement could serve as a surrogate marker for aneurysm instability.
Data availability The datasets analyzed in the current study are available from the corresponding author on reasonable request.

Compliance with ethical standards
Conflict of interest The authors declare that they have no conflict of interest.
Ethical approval All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.
Approval was obtained from the Ethics Committee of the Christian-Albrechts-University Kiel, Germany.
This article does not contain any studies with animals performed by any of the authors.
Informed consent Informed consent was obtained from all individual participants included in the study.
Code availability Not applicable.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.