Vascular tortuosity of the internal carotid artery is related to the RNF213 c.14429G > A variant in moyamoya disease

Recent studies have implicated RNF213 mutations in the pathogenesis of moyamoya disease (MMD). However, the underlying mechanism of disease development is not fully elucidated. Nonetheless, a possible relationship between vascular morphology and hemodynamics related with MMD has been proposed. Here, we aimed to investigate the relationship between a variant of RNF213 and the morphology of the internal carotid artery (ICA). We enrolled bilateral MMD patients who had undergone genetic testing for RNF213. Patients were divided into mutant and wild-type groups. Six anatomy-specific three-dimensional coordinates were collected using magnetic-resonance angiography. From these, five vectors between two adjacent points and four angles between two adjacent vectors were calculated. The tortuosity was defined as the ratio between the actual and the linear length of the ICAs. Among 58 patients, 44 and 14 belonged to the mutant and wild-type groups, respectively. The tortuosity of ICAs was significantly lower in the mutant group (p = 0.010). The change in blood flow direction was more prominent in the wild-type group (p = 0.002). The tortuosity was significantly lower in MMD patients than normal controls (p < 0.001). Our results indicate that RNF213 could play a role in the lower tortuosity observed in patients with RNF213 mutation.


Methods patients. This study was approved by the institutional review board of Seoul National University Bundang
Hospital (B-1902/522-112), and informed consent from patients was waived. We screened our institution's electronic medical records and identified patients who (1) had been diagnosed with MMD based on digital subtraction angiography, (2) had undergone genetic examination for RNF213 from November 2017 to October 2018, and (3) had undergone time-of-flight magnetic resonance angiography (TOF-MRA). Patients with unilateral MMD or a history of previous revascularization surgery were excluded because these factors could alter the flow through the ICA 11 . The same number of controls matched by age and sex were screened among those with normal MRA readings. Patients were divided into a mutant group (patients with the RNF213 c.14429G > A mutation, regardless of its homozygous or heterozygous state) and a wild-type group. Two ICAs from each subject were analyzed as independent vessels. three-dimensional measurement. The source images used in this study were derived from TOF-MRA.
The original Digital Imaging and Communications in Medicine (DICOM) files were imported into a 3D image analysis tool (Slicer 4.10.1; open source software, www.slicer.org/). Six anatomy-specific 3D coordinates under the RAS axes system (Right side, Anterior side, and Superior side as corresponding for X, Y, and Z coordinates) within each intracranial ICA were defined as follows: P1, the beginning of the vertical petrous ICA at the level of the carotid canal; P2, the inflection point between the vertical petrous ICA and the horizontal petrous ICA; P3, the inflection point between the horizontal petrous ICA and the vertical cavernous ICA; P4, the inflection point between the vertical cavernous ICA and the horizontal cavernous ICA; P5, the anteriormost point of the carotid siphon; and P6, the point of emergence of the anterior choroidal artery. Each point was marked at the center of the vascular lumen in the corresponding position.
Using 3D coordinates, five vectors and angles were defined. . As a person has two ICAs, each measured value from the right and left sides was averaged as an individual value. The more advanced Suzuki grade between the right and left side was used as the representative grade.  www.nature.com/scientificreports www.nature.com/scientificreports/ range; IQR) because some variables showed non-normal distribution by the Kolmogorov-Smirnov test. Those variables were analyzed by the Kruskal-Wallis test and Wilcoxon's rank-sum test. Bonferroni's correction was also applied for significant variables in the Kruskal-Wallis test. All length variables were measured in millimeters, whereas angle variables were in degrees. A p-value of < 0.05 was considered statistically significant.

Results patient characteristics.
In one year (from November 2017 to October 2018), a total of 208 patients had undergone genetic tests for the RNF213 gene (c.14429G > A mutation). Among them, 70 and 45 patients were excluded due to the lack of cerebral angiography and MRA data, respectively. Among the remaining 93 patients, 27 patients with unilateral MMD or those who underwent any prior revascularization surgery were excluded. The resultant 66 patients were divided into two groups. The mutant group comprised 52 (78.8%) patients with heterozygotic (GA) mutation.; the wild-type (without mutation) group comprised 14 (21.2%) patients. There was no patient with homozygotic mutation (AA). All the patients included in this study were Koreans. When medical records were reviewed, there was a single MMD patient (1.5%) with neurofibromatosis type 1 in mutant group. Eight patients with Suzuki grade 1 or 5 were excluded in the main analysis, because all wild-type group patients had Suzuki grade of either 2, 3, or 4.
The [27/44]) was not significantly different from that in the wild-type group (57.1% [8/14]; p = 1.000). More patients with advanced Suzuki grade were in the mutant group; however, this was not a statistically significant finding (p = 0.056). The age and sex ratios were exactly the same between MMD patients and control subjects. tortuosity analysis between mutant and wild-type RNF213. As to segmental lengths, → V23 was significantly shorter (p = 0.001) in the mutant than in the wild-type group. Other segmental lengths were similar. As to the sum of segmental lengths, ALEN was significantly shorter in the mutant group (53.80 [50.38-58.41]) than in the wild-type group (60.31 [56.20-61.00], p = 0.023). Considering the similar length of LLEN between the two groups (p = 0.518), the tortuosity of ICAs was calculated to be significantly lower in the mutant group (p = 0.010) (Fig. 2). Measured angles, ∠123 and ∠234, showed that the angle change of the ICA was significantly lower in the mutant group. Other angles, ∠345 and ∠456, were similar in each comparison. The total change in blood flow direction was more prominent in the ICAs of wild-type group patients (p = 0.002, Fig. 3). tortuosity analysis among MMD patients and controls. In pair-wise comparisons, mutant MMD patients had a greater angle change than control individuals (p < 0.001) and wild-type MMD patients (p = 0.005). However, the total change in angle was similar between wild-type MMD patients and control individuals (p = 0.654). In short, although the tortuosity of wild-type MMD patients was lower than that of control subjects, the difference was more pronounced in the mutant MMD group.    .1]; p = 0.002). As for Suzuki grade, the age distribution was similar among grade groups. The overall tortuosity and change in angle tended to decrease with a higher grade, although the differences between adjacent grades were not statistically significant except between grade 3 and 4 (Fig. 4).

Discussion
The protein RNF213, encoded by the RNF213 gene, is known to be associated with angiogenesis 12 . RNF213 deficiency in the endothelium leads to the overexpression of several matrix metalloproteinase genes. This deficiency has deleterious effects on endothelial cells in the maintenance of vascular structure 13 . Sonobe et al. revealed that mice lacking RNF213 showed suppressed vascular remodeling after common carotid artery ligation 14 . The ischemic insult caused increased vascular expression of matrix metalloproteinase-9 and thinning of the vascular wall distal to the ligation. Similarly, RNF213 knockout mice showed enhanced angiogenesis following hindlimb ischemia. This enhancement in angiogenesis corresponds to the abnormal cerebral vascular network in the chronic ischemic state. However, after cerebral ischemia, no significant difference in brain angiogenesis was found compared with the wild-type 15 .
However, vasculopathies that are strongly related to genetic mutation tend to involve multisystemic vessels. Since MMD inarguably involves a specific part of the vessel, some efforts have been made to study the RNF213 c.14429G > A mutation and its relationship with intracranial vascular morphology 16,17 . However, the link between the vasculopathies and the specific location of the MMD culprit, i.e., the distal ICA, remains unsolved. Computational fluid dynamics demonstrated that lower ICA tortuosity results in higher wall shear stress, especially at the terminal ICA 6,7 . Thus, the hypothesis derived here was that the RNF213 c.14429G > A mutation is related to distal ICA stenosis/occlusion of MMD through the intracranial ICA phenotype involving tortuosity.
We showed the morphological differences of the ICA according to the mutation status of RNF213 c.14429G > A. Although the most accurate method for calculating the tortuosity of vessels is the extraction of centerlines from 3D-rendered vascular models produced by manual segmentation of a source image, its analysis requires excessive time and manpower 7 . Thus, in this study, we adopted a method for defining six landmark points and calculating the straight lines connecting the adjacent points, thereby approximating the trajectory of the actual vessels. By plotting 3D coordinates using the patients' MRA, we could obtain segmental length, angle, and tortuosity of the vessel close to real values, such that an acceptable data comparison could be made. Through the analysis, we could show the lower tortuosity and flow angle change of the ICA in MMD patients with mutated RNF213 c.14429G > A compared with those in wild-type MMD patients. Considering the relationship among tortuosity, RNF213 mutation, and Suzuki grade, it could be inferred that patients with RNF213 mutation would be diagnosed with MMD in more advanced status, and that differences in one of the phenotypes would result in vascular tortuosity.
Numerous studies have been carried out to investigate the pathogenesis of MMD, and complex regulation of angiogenesis or vasculopathy have been revealed, with some contradictory results between studies 18 . If gross vascular morphological characteristics related with the RNF213 c.14429G > A mutation can be demonstrated, various angiogenic or vasculopathy properties of MMD can be understood as a gross phenomenon. Our study successfully revealed the relationship between the RNF213 c.14429G > A mutation and greater ICA tortuosity.
A crucial limitation of this study was related to the design of the study. Since we performed this investigation as a retrospective case-control study at a certain time, we could only show correlations between the morphological differences of the ICA and the RNF213 c.14429G > A mutation and were unable to prove causality. Serial tortuosity of the ICA and cerebral perfusion need to be measured throughout the disease course for validating the causality. Another limitation was the modest sample size. Although over 200 subjects underwent genetic testing for the RNF213 c.14429G > A mutation, almost 68% were excluded to secure a homogenous study population (bilateral MMD, availability of cerebral angiography and MRA). As only 14 subjects were allocated to the wild-type group with Suzuki grade 2-4, additional eight subjects were also excluded from mutant group who presented as Suzuki grade 1 or 5. The distributions of sex and age were statistically similar between the two groups. Although the distribution of Suzuki grade showed statistical homogeneity, the p-value was marginal. Thus, we www.nature.com/scientificreports www.nature.com/scientificreports/ performed additional analysis among the tortuosity, angle change, and Suzuki grade which showed that ICA tortuosity tended lower in earlier Suzuki stage patients.
Another potential weakness of this study was related to the simplified methodology of measurement. A previous tortuosity measurement was made by actual measurement after 3D centerline extraction following segmentation and rendering 7 . Although 3D segmentation and centerline extraction is the most accurate method to measure vascular morphology, it requires considerable effort. Our methodology used a simplified method approximating real measurement and required the measurement of only six fiducial points, which could be easily obtained using an open source software. Although the accuracy might be slightly decreased, the range of measured tortuosity seemed to be acceptable compared with that of a previous report, considering the longer segment investigated (to the point of emergence of the anterior choroidal artery) than that in the earlier report (to the carotid siphon) 7 .
Together, our results indicate that RNF213 might be involved in the shaping of the ICA, with RNF213 deficiency resulting in terminal ICA stenosis/occlusion through focal high wall shear stress from the lower tortuosity of conduits 6,7 . Further serial observational studies using a larger disease cohort are warranted for further clarification of the causal relationship between the morphology of the ICA and MMD.

Data Availability
The data are not available for public access because of patient privacy concerns but are available from the corresponding author on reasonable request approved by the institutional review boards of Seoul National University Bundang Hospital.