Abstract
BACKGROUND AND PURPOSE: Previous studies have reported that MCA bifurcation aneurysms usually emerge on inclined bifurcations; however, the reason is unclear. We designed this study to explore hemodynamic mechanisms that correlate with the initiation of MCA bifurcation aneurysms.
MATERIALS AND METHODS: Fifty-four patients with unilateral MCA bifurcation aneurysms and 54 control patients were enrolled in this study after propensity score matching, and their clinical and CTA data were collected. We extracted the morphologic features of aneurysmal MCA bifurcations to build a simplified MCA bifurcation model and performed a computational fluid dynamics analysis.
RESULTS: The presence of MCA aneurysms correlated with smaller parent-daughter angles of MCA bifurcations (P < .001). Aneurysmal MCA bifurcations usually presented with inclined shapes. The computational fluid dynamics analysis demonstrated that when arterial bifurcations became inclined, the high-pressure regions and low wall shear stress regions shifted from the apexes of the arterial bifurcations to the inclined daughter arteries, while the initial sites of MCA bifurcation aneurysms often overlapped with the shifted high-pressure regions and low wall shear stress regions.
CONCLUSIONS: Our results suggest that the initiation of MCA bifurcation aneurysms may correlate with shifts of high-pressure regions and low wall shear stress regions that occur on inclined MCA bifurcations.
ABBREVIATIONS:
- CFD
- computational fluid dynamics
- HPR
- high-pressure region
- LWSS
- low wall shear stress
- LWSSR
- low wall shear stress region
- PSM
- propensity score matching
- RD
- ratio of diameter
- ROC
- receiver operating characteristic
The occurrence of intracranial aneurysms is generally thought to be due to arterial wall weakening and/or the influence of hemodynamics.1⇓⇓–4 Most studies on aneurysm etiology have focused on hemodynamics. Because the geometric shapes of intracranial arteries are diverse, their hemodynamics are also complex. Arterial bifurcations are common sites for intracranial aneurysms. The bifurcations are generally exposed to extreme hemodynamic stress; those with special morphologic features that significantly divert blood flow from the direction of flow of the parent vessels are a risk factor for aneurysm formation.1,5,6 Previous studies have reported that the inclined MCA bifurcations with a widening angle are likely to harbor aneurysms,7,8 but mechanisms by which the aneurysms develop on MCA bifurcations are controversial and need to be elucidated.8⇓⇓-11 The present study was designed to extract the morphologic features of aneurysmal MCA bifurcations from clinical data, build a simplified MCA bifurcation model based on these extracted features, perform computational fluid dynamics (CFD) analysis on the aneurysmal MCA bifurcations and simplified models, and finally, to explore the hemodynamic mechanisms that could trigger aneurysm initiation.
MATERIALS AND METHODS
Case Selection and Matching
Patients with unilateral MCA bifurcation aneurysms confirmed by CTA from July 2016 to March 2019 were enrolled in our study. Exclusion criteria were as follows: 1) fusiform or dissecting aneurysms, 2) trifurcation types of MCAs, and 3) aneurysms of >10 mm. In addition, we reviewed healthy subjects without intracranial aneurysms from the general population during the same time period. Of the 162 patients enrolled in our study, 54 with unilateral MCA aneurysms were assigned to the aneurysm group and 108 without intracranial aneurysms, to the control group. General data including age, sex, smoking, drinking, hypertension, hyperlipidemia, and diabetes were collected. Propensity score matching (PSM) was used for selecting the cases in the 2 groups, and the variables listed above were used as matching parameters. The matched cases included 54 with unilateral MCA aneurysms and 54 without intracranial aneurysms. A flow diagram of our study is shown in Fig 1. All procedures in this retrospective study that involved human participants were approved by the ethics committee of our hospital and were performed in accordance with the institutional ethics standards, the 1964 Declaration of Helsinki and its later amendments, or comparable ethical standards. Informed consent was obtained from all individual participants included in the study.
Morphologic Features
3D CTA data were reconstructed using Materialise Mimics software (Version 21.0; https://www.materialise.com/en/medical/mimics-innovation-suite/mimics), and morphologic data were measured using the Materialise 3-matic research software (Version 13.0; https://www.materialise.com/en/software/3-matic). Two certified neuroradiologists (F.Y.Z. and C.Q.C.), who were blinded to the clinical data, independently evaluated the morphologic data and resolved disagreements by consensus. The parent artery was defined as the M1 trunk of the MCA, and the daughter artery, as the M2 branch of the MCA.
The measuring methods of 6 angles in the control group are shown in Fig 2A. The measuring methods of diameters are depicted in Fig 2D. Four ratios of diameters (RDs) were used for analysis. They were separately defined as RDL1, the ratio of the left upper daughter-parent artery diameters; RDL2, the ratio of left lower daughter-parent artery diameters; RDR1, the ratio of right upper daughter-parent artery diameters; and RDR2, the ratio of right lower daughter-parent artery diameters.
The measuring methods of 6 angles in the aneurysm group are described in Fig 2B, and the diameter measuring methods are shown in Fig 2D. The ratio of parent-daughter artery diameters was separately termed as RDA on the aneurysm side and RDNA on the nonaneurysmal side. For the contralateral side, the ratio of upper daughter-parent artery diameters was defined as RDC1, and the ratio of lower daughter-parent artery diameters was RDC2. The method for deciding the location of an MCA bifurcation aneurysm is described in Fig 2C.12 The long axis of an aneurysm was measured as aneurysm size. Aneurysm neck width was also measured.
CFD Analysis
We reconstructed 3D MCA bifurcation models in Mimics research software (Materialise NV) for all 54 cases with MCA aneurysms and 20 cases in the control group. Then, the aneurysmal MCA bifurcation models were virtually removed using the 3-matic research software, and their initial status was simulated. All 3D models were exported to Fluent solver (Version 19.1; ANSYS) for CFD analysis.
After extracting some hemodynamic features from the above CFD analyses, we constructed a group of simplified 3D MCA bifurcations models with variable single parent-daughter angles. For these models, parent-daughter angle of the MCA bifurcations was set to 120° (based on the mean value of ϕNA, 118.81° ± 21.35°), and another was varied from 120° to 30° with intervals of 10° (based on the range of ϕA, 31.95°–124.89°). The diameters of the parent and daughter arteries were set to 4 and 3 mm, respectively. For all 3D models, blood flow was modeled as a laminar Newtonian fluid, with a density of 1050 kg/m3 and dynamic viscosity of 0.0032 Pa × s. The inlet boundary of the parent artery was defined as the mass-flow inlet at 0.0034 kg/s, and the outlet boundaries of 2 daughter arteries were defined as zero-static pressure.8,13,14
Statistical Analysis
Data were analyzed using R Studio (Version 1.1.383; http://rstudio.org/download/desktop). Student t tests and ANOVA analyses were used to compare continuous variables, and the Pearson χ2 test was used for comparing categoric variables. Pearson correlation tests were used to analyze correlations between 2 continuous variables. Multiple variables were analyzed using logistic regression. To predict which daughter artery was prone to developing an MCA aneurysm, receiver operating characteristic (ROC) curve analysis was performed. P values < .05 were considered statistically significant.
RESULTS
Demographics
The clinical data of both groups before and after PSM are summarized in Table 1. Of 54 subjects in the control group, 28 were women and 26 men, with a mean age of 57.39 ± 8.79 years. Of 54 patients with aneurysms, 36 were women and 18 men, with a mean age of 60.35 ± 9.87 years. Twenty-two (40.74%) patients were confirmed to have ruptured aneurysms, and 17 (31.48%), to have multiple aneurysms. Logistic regression analysis based on morphologic parameters of MCA bifurcations showed that the presence of MCA aneurysms correlated with the parent-daughter angle (Table 2).
Morphologic Features of Nonaneurysmal and Aneurysmal MCA Bifurcations
There was no statistically significant difference between groups among the 7 parent-daughter angles (ϕL1, ϕL2, ϕR1, ϕR2, ϕNA, ϕC1, ϕC2; Table 3 and Fig 3A). The daughter-daughter angles (γA and γC) in patients with aneurysms were significantly larger than those (γL and γR) in the control group (Table 3 and Fig 3B).
For the aneurysm group, parent-daughter angles on the inclining side of the aneurysm (ϕA) were significantly smaller than the other 3 parent-daughter angles (ϕNA, ϕC1, ϕC2) (Fig 3A), whereas the daughter-daughter angles on the aneurysm side (γA) were significantly larger than the angle of the contralateral side (γC) (Fig 3B). The mean aneurysm size was 5.05 ± 1.98 mm, and the mean aneurysm neck diameter was 4.68 ± 1.54 mm. The optimal parent-daughter angle threshold to predict the daughter artery that was likely to harbor an MCA aneurysm was determined by ROC analysis. The ROC curve showed that the optimal threshold was 100.06 (area under the curve, 0.903), with 84.1% sensitivity and 81.5% specificity (Fig 3C).
Influence of Aneurysm Rupture on MCA Bifurcation Morphology
There was no significant difference between the parent-daughter and daughter-daughter angles or the RDs on aneurysmal and nonaneurysmal sides of the ruptured and unruptured aneurysmal MCA bifurcations (Table 4).
Influence of Aneurysm Size and Neck Width on Parent-Daughter Angles of MCA Bifurcations
Correlation analysis showed that the sizes and neck widths of MCA aneurysms were not significantly correlated with parent-daughter angles (On-line Figure).
Hemodynamic Features of Aneurysmal MCA Bifurcations and Simplified MCA Bifurcation Models
The CFD analysis performed on all MCA aneurysms showed that after the aneurysms were virtually removed, compared with control MCA bifurcations (On-line Digital Content 1; https://yaleedu-my.sharepoint.com/:b:/g/personal/fenghua_chen_yale_edu/EV5HuHzoJrRBsp1nO86pKgQBxPemPz4zSsXWBH2Ap-eR3w?e=u79Us0), their initial sites presented high-pressure regions (HPRs) and low-wall shear stress regions (LWSSRs) (Fig 4A and On-line Digital Content 2; https://yaleedu-my.sharepoint.com/:b:/g/personal/fenghua_chen_yale_edu/EV5HuHzoJrRBsp1nO86pKgQBxPemPz4zSsXWBH2Ap-eR3w?e=u79Us0). The CFD analysis of simplified 3D MCA bifurcation models revealed that when a single parent-daughter angle decreased from 120° to 30°, the HPR and LWSSR shifted from the apex of the bifurcation to the daughter artery wall on the same side, and the maximum shifted distance was 2.03 mm (Fig 4B–D and On-line Digital Content 3; https://yaleedu-my.sharepoint.com/:b:/g/personal/fenghua_chen_yale_edu/EV5HuHzoJrRBsp1nO86pKgQBxPemPz4zSsXWBH2Ap-eR3w?e=u79Us0).
DISCUSSION
The mechanism of intracranial aneurysm initiation remains unclear. Intracranial arterial bifurcations frequently harbor aneurysms. When blood flow enters an arterial bifurcation, it exhibits complex hemodynamic features, some of which may correlate with aneurysm formation.2,5,15,16 In this study, we analyzed a series of MCA bifurcations with and without aneurysms and explored the hemodynamic factors involved in aneurysm development. We found that MCA aneurysms were commonly located on the inclined MCA bifurcations, consistent with previous studies.8,10 Our results revealed that inclined MCA aneurysmal bifurcations are usually caused by smaller single parent-daughter angles, but other parent-daughter angles maintain sizes similar to the 2 parent-daughter angles on the contralateral “normal” MCA bifurcations. CFD analysis based on the above morphologic features showed that when MCA bifurcations become inclined, the HPRs and LWSSRs shift from the apexes of the arterial bifurcations to the ipsilateral daughter arteries, and this hemodynamic feature may correlate with the initiation of MCA bifurcation aneurysms.
Considering that the parent-daughter angle may be affected by the aneurysm size, we enrolled only patients with small MCA aneurysms. We analyzed the correlation between the parent-daughter angle of the MCA bifurcation and aneurysm size and neck width. The results suggested that the size of the parent-daughter angle on the aneurysm side might not be affected by aneurysm growth. Our data also showed that aneurysm rupturing had little influence on MCA bifurcation morphology.8 Unlike some previous groups that based their measurements on 2D projecting images of MCA bifurcations,8,17 we used 3D measuring tools available in Mimics research software, so we did not adopt all the parameters used in those studies. We considered that by using the 3D coordinate, our data would be more representative of the real conditions in MCA bifurcations.
Some studies have reported that the daughter-daughter angle is correlated with the presence of an intracranial aneurysm,8,10 which is consistent with our results. We found that of all the parent-daughter angles in the aneurysm group, the only difference was observed in the parent-daughter angle (ϕA) on the aneurysm inclining side. This suggests that the difference between the daughter-daughter angle (γA) on the aneurysmal side and the contralateral daughter-daughter angle (γc) is largely attributable to the parent-daughter angle (ϕA). Moreover, ROC curve analysis showed that the parent-daughter angle (ϕA) was useful for predicting which daughter artery is likely to harbor an MCA aneurysm. These results suggest that a single smaller parent-daughter angle may play an important role in MCA aneurysm occurrence. Prospective studies should be performed to test the predictive value in subjects with smaller parent-daughter angles of MCA bifurcations.
To explore the possible hemodynamic mechanism behind this unique shape, we further performed CFD analysis on 3D aneurysmal MCA bifurcations reconstructed from CTA data. By comparing the differences before and after virtual MCA aneurysm removal, we observed that the initiation sites of aneurysms overlapped with HPRs and LWSSRs, and these areas inclined to the daughter artery that had a smaller parent-daughter angle. Because hemodynamic data of CFD analysis on aneurysmal bifurcations are difficult to compare, we constructed a group of simplified bifurcation models to analyze hemodynamic changes on varying a single parent-daughter angle. The CFD analysis of these simplified models reconfirmed the previous finding in the 3D models of actual aneurysmal MCA bifurcations—that is, when one of the parent-daughter angles becomes smaller, the HPR and LWSSR shift from the apex of the arterial bifurcation to the daughter artery. We also observed that the HPR range widened and the distal high-pressure gradient became steeper, but the pressure and wall shear stress were not dramatically altered. Combining these results with the fact that the inclined HPR and LWSSR overlap with an MCA aneurysm site, we inferred that the inclined HPRs and LWSSRs may correlate with the initiation of the aneurysm.
According to previous studies, high or low wall shear stress (LWSS) may cause aneurysm formation.1,3,4,18 The wall shear stress curve data supported the hypothesis that LWSS was a possible reason. Similar LWSS values also exist on the apexes of normal MCA bifurcations, however, the normal distribution of LWSS seldom causes aneurysms. Rather, the abnormal shift of LWSS may be the real reason. Why can the HPR and LWSSR on the daughter arterial wall form an aneurysm but cannot when they are located on the apex? Previous studies offer some clues.3,18⇓⇓⇓-22 Specifically, unlike the arterial wall, the apex of a cerebral arterial bifurcation lacks a tunica media and forms a “gap” that is filled with fibers from the tunica adventitia, namely the “apical ridge.” This structure is composed of collagen fibers and can, therefore, bear stronger impingement of blood flow than the arterial wall.
On the basis of the above inferences, we propose the following mechanism of MCA aneurysm initiation. Compared with symmetric MCA bifurcations in which high-pressure and LWSS fall on the “strong” apical ridge, the high-pressure and LWSS of asymmetric MCA bifurcations fall on the relatively “weak” arterial wall. Because the arterial wall cannot bear this stress, arterial structures—especially the elastic laminas—will be damaged under the long-term impact of blood flow and an aneurysm will emerge (Fig 5). Although this hypothesis considers the influence of the MCA bifurcation structure on the formation of an intracranial aneurysm to some extent, it also emphasizes the core role of hemodynamics. Whether these events are universal is unknown, and our theory still needs to be examined and confirmed by further investigations based on biologic models.
Our study has some limitations that need to be considered. First, we examined cross-sectional data. Although our analysis showed no significant correlation between parent-daughter angles on the aneurysmal bifurcations and aneurysm size or neck width, the local geometry of the arteries may still be affected, which would reduce the accuracy of hemodynamic analysis. To minimize the effect of aneurysm size and represent the initial status of aneurysm formation, only small MCA aneurysms of <10 mm were included in our study; therefore, our results may not be generalizable to all MCA aneurysms. Second, our findings revealed a stronger correlation between the parent-daughter angle of the MCA bifurcation and the presence of aneurysms; however, it is unclear whether this is a cause-and-effect relationship. Future studies should be performed to help identify high-risk individuals who are likely to develop aneurysms.
CONCLUSIONS
An inclined MCA bifurcation caused by a single smaller parent-daughter angle is prone to develop an aneurysm. When a single parent-daughter angle of an MCA bifurcation decreases, the HPR and LWSSR shift from the apex of the arterial bifurcation to the daughter artery wall. This hemodynamic change may be associated with the initiation of an MCA bifurcation aneurysm.
ACKNOWLEDGMENTS
The authors thank Jun Zhou, PhD (College of Chemistry and Chemical Engineering, Hunan University), for assistance with CFD analysis.
Footnotes
This work was supported by National Natural Science Foundation of China (grant No. 81873756).
Disclosures: Zheng Huang—RELATED: Grant: National Natural Science Foundation of China (grant No. 81873756), Comments: The grant was used to buy some software for this study.* Feiyue Y. Zeng—RELATED: Grant: National Science Foundation of China (grant No.81873756).* Longbo B. Zhang—RELATED: Grant: National Natural Science Foundation of China (grant No.81873756).* Fenghua H. Chen—RELATED: Grant: National Natural Science Foundation of China (grant No. 81873756).* *Money paid to the institution.
Indicates open access to non-subscribers at www.ajnr.org
References
- 1.
- 2.
- 3.
- 4.
- 5.
- 6.
- 7.
- 8.
- 9.
- 10.
- 11.
- 12.
- 13.
- 14.
- 15.
- 16.
- 17.
- 18.
- 19.
- 20.
- 21.
- 22.
- Received November 29, 2019.
- Accepted after revision February 22, 2020.
- © 2020 by American Journal of Neuroradiology