Patient
This study was approved by the Ethics Committee at Beijing Tiantan Hospital. A 31-year-old male patient presented with sudden posterior headache without any apparent cause. The first HR-MRI examination was performed and revealed a superior basilar aneurysm with Grade-1 of the aneurysmal wall enhancement. Considering the relatively high risk of endovascular treatment for this case, conservative treatment was applied. After seven months, the patient came for re-examination with progressive headache, and Grade-2 of AWE with obvious IA enlargement were found. Unfortunately, in one month, the aneurysm ruptured. All of the imaging examinations were performed via a 3.0-T MR scanner (Achieva TX, Philips, Netherlands) . The legal representatives of the participant gave their informed consent.
Reconstruction of Aneurysms and Enhancement Sections and Model co-registration
The patient-specific aneurysm models were reconstructed based on the TOF MRI images and the enhancement regions were extracted based on the T1-Vista sequence, via a semi-automatic image processing tool, MIMICS (Materialize, Leuven, Belgium). The lumen profiles of the aneurysm and the segmented enhancement regions were confirmed by two experienced neurointerventional experts. Aneurysm/enhancement models at the first and follow-up examinations, denoted as A-T1/E-T1 and A-T2/E-T2 models, were thus confirmed (Figure 1A). The longitudinal models were registried via the reference point system alignment method. Optimal fit alignment of two models was achieved by a three-dimensional transformation matrix shown in Figure 1B. In this matrix, Aij, i, j=1,2,3 controls the rotation of the 3D model, [Px Py Pz]T is the perspective transformation matrix, [tx ty tz] denotes the translation and S is the overall scale factor with a default value of 1. Furthermore, due to the various scanning resolution, medical images showing enhancement regions contain 6 slices and 34 slices for E-T1 and E-T2 model, respectively (Figure 1C and D).
Morphological Features and Growth Quantification
Imaging follow-ups facilitate morphological analysis on growing aneurysms. However, the hemodynamic and pathological mechanism that promotes aneurysm growth is still unclear. To precisely quantify the aneurysm growth by digital representation, two quantitative indices, Displacement Index (DI) and Area Stretch Ratio (ASR)[21], were calculated. DI represented the Eulerian distance between two aneurysm models in the direction of the normal vector at each voxel on A-T1 model. The ASR provided a more sophisticated descriptor of aneurysm growth compared to DI and the growth was considered as the stretching of the A-T1 surface over time. This parameter was quantified as the ratio of the A-T2 surface area to that at A-T1, which was a function of the displacement between the A-T1 and A-T2 surfaces.
Hemodynamics Modeling
The reconstructed aneurysm geometries were meshed in ICEM CFD with tetrahedron elements in the core region surrounded by 5-layer prismatic elements to reinforce the grid resolution near the vessel wall. Hemodynamic computations were conducted on A-T1 and A-T2 models after 3D registration. The blood flow was assumed as laminar, incompressible and Newtonian flow with a constant viscosity of 0.0035 Pa.s and constant density of 1066 kg/m-3. The property of arterial wall was set as rigid and non-slip boundary. The unstructured mesh was imported into ANSYS CFX 19.2 (ANSYS Inc. Canonsburg, PA, USA) to solve the flow governing Naviér-Strokes equation. This specific boundary condition was extracted from the software @neufuse which possesses the database of a 1D model of human vascular tree[22]. Computations were conducted for three cardiac cycles to guarantee the numerical stability and the results of last cardiac cycle were applied for further data analysis.
Signal Intensity Standardization of Enhancement
Jorge A. Roa et al. demonstrated three measurement methods of the signal intensity (SI) normalization and revealed that the SI value of aneurysm enhancement scandalized by pituitary stalk achieved higher sensitivity and specificity, and allowed enhancement measurements standardization in HR-MRI [23]. In the currently study, the SI of E-T1 and E-T2 models were extracted and compared with the pituitary stalk to standardize the SI of enhancement area, which were defined in Equation 1, where SI-ave-pi denoted the average SI of pituitary stalk.
stanSI =SI-en / SI-ave-pi (1)
The E-T1 and E-T2 models were mapped with stanSI to obtain the enhanced 3D model with SI for subsequent visualizations.
Thickness and SI Quantification of Enhancement on Single Layer
To explore the correlation between hemodynamic and enhancement characteristics (Aim 1), i.e. SI and thickness, hemodynamic results on aneurysm lumen and enhancement features on each slice were extracted (40 slices). Furthermore, six new slices of enhancement in the post-registered E-T2 model were extracted aligned with the enhancement layers in the E-T1 model, which aims to investigate the effective markers for aneurysm local growth (Aim 2). Therefore, a total of 46 slices were enrolled for Aim 1 and only 6 slices were included for Aim 2 due to the limited layers in the E-T1 model.
To satisfy the analysis requirements of above two aims, the key procedure is the quantification of the thickness and SI of enhancement regions and the establishment of voxel-based match of hemodynamics and enhancement. Firstly, the 3D coordinates of enhancement region on slice 17 were selected to determine the origin of the local coordinate system because of its medial position. By using the principal component analysis (PCA) method, two largest principal components that were mutually orthogonal were calculated as the X/C1-axis and Y/C2-axis of the local coordinate system (Figure 2A). For each single slice, the localized origin O(n) was determined as the mean value of the coordinates on the outline of the aneurysm lumen. Taking the positive Y-axis as the starting direction, the outline of aneurysm lumen was clockwise discretized into 200 intersection points with an interval of = /100. Hemodynamic information, i.e. time averaged wall shear stress (TAWSS), OSI, gradient oscillatory number (GON), aneurysm formation index at middle diastole (AFImd), pressure at systole (Ps), at all points were extracted. Along with each direction, the enhancement thickness, maximum gray value and average gray value were also calculated.
Statistical Analysis
Significance and correlation test for all statistical data were conducted using SPSS (IBM, Armonk, New York, USA). Mann-Whitney test was performed to analyze the differences of hemodynamic variables between AWE patterns. P <0.05 was considered statistically significant. Shapiro-Wilk test was performed to verify whether the parameters conform to normal distribution. For nonnormally distributed parameters, Spearman correlation coefficient was performed to analyze the correlation between hemodynamic variables and AWE.