Dear Sir:
Acute ischemic stroke (AIS) is a leading cause of disability and mortality, warranting detailed investigation into its underlying mechanisms. A major subtype of AIS is atherothrombotic brain infarction, which often results from intracranial atherosclerotic disease (ICAD). In addition to traditional atherosclerotic risk factors, emerging evidence suggests a role for genetic contributors [1]. The p.R4810K variant of the ring finger protein 213 gene (RNF213, mysterin), a susceptibility gene for Moyamoya disease, has also been associated with ICAD, particularly in the anterior circulation [2]. This variant is relatively prevalent among East Asian populations, especially in Japan, where approximately 1%-2% of the general population carries the RNF213 p.R4810K allele [2]. Consequently, the “two-hit hypothesis” has been proposed, suggesting that additional genetic or environmental factors are necessary to trigger vascular stenosis. The search for such secondary factors is ongoing [3]. This study aimed to evaluate the association between RNF213 p.R4810K and ICAD and to investigate, in line with the two-hit hypothesis, effect modification by traditional atherosclerotic risk factors, using data from a multicenter stroke registry.
We analyzed 14,471 patients with AIS or transient ischemic attack enrolled in the Fukuoka Stroke Registry [4]. Details of patient selection are provided in the Supplementary Methods and in the Supplementary Figure 1. The study design was approved by the Institutional Review Boards of all participating hospitals. Written informed consent was obtained from all patients or their legal representatives. Genetic testing for RNF213 p.R4810K (rs112735431) was performed using the Melt Analysis of Mismatch Amplification Mutation Assays [5]. Patients were classified as wild-type (G/G) or variant carriers (G/A or A/A). ICAD was defined as ≥50% stenosis or occlusion of major intracranial arteries on magnetic resonance angiography or computed tomography angiography [6], as initially interpreted independently at each participating hospital by at least one radiologist and one stroke neurologist. For the sensitivity analysis, intracranial atherosclerotic stenosis (ICAS) was defined as ≥50% stenosis without occlusion.
Logistic regression analysis was used to estimate odds ratios (ORs) and 95% confidence intervals for each outcome. Exploratory interaction analyses were also conducted using multiplicative terms to assess potential effect modification by traditional atherosclerotic risk factors. Body mass index (BMI) was categorized into four groups (<18.5, 18.5-24.9, 25.0-29.9, and ≥30.0 kg/m2) [7]. Heterogeneity was further explored by stratifying the analyses by BMI. Patients were divided into eight subgroups based on RNF213 p.R4810K carrier status and BMI category. Interaction effects were tested by including multiplicative interaction terms in the corresponding logistic regression models. Two-tailed P values <0.05 were considered statistically significant. All statistical analyses were performed using STATA version 16.0 (StataCorp., College Station, TX, USA). Further details on the ethics statement, genotyping procedures, ICAD definitions, clinical variables, and statistical analyses are provided in the Supplementary Methods.
The mean age of the 14,471 patients was 73 (±12) years, and 8,593 (59.4%) were men. Of all participants, 14,100 (97.4%) were non-carriers and 371 (2.6%) were carriers of the p.R4810K variant (Supplementary Table 1). A total of 3,347 patients (23.1%) had ICAD, including 2,358 (16.3%) with anterior circulation lesions (anterior-intracranial atherosclerotic disease [a-ICAD]) and 1,686 (11.7%) with posterior circulation lesions (posterior-intracranial atherosclerotic disease [p-ICAD]). Among those with ICAD, 150 were variant carriers and 3,197 were non-carriers (Supplementary Table 2). The RNF213 p.R4810K variant was significantly associated with an increased risk of ICAD, including both a-ICAD and p-ICAD, even after adjusting for age, sex, and vascular risk factors (Supplementary Table 3). When stratified by covariates, significant heterogeneity was observed in the ORs for ICAD among variant carriers by age (P for heterogeneity <0.001) and BMI categories (P=0.027) (Supplementary Table 4). Table 1 presents baseline characteristics by BMI category, including the distribution of RNF213 variant carriers and ICAD prevalence. Both RNF213 variant carrier and ICAD prevalence tended to increase with higher BMI. In the subgroup analysis across the four BMI categories, the risk of ICAD among variant carriers increased progressively with increasing BMI. Significant heterogeneity was observed for overall ICAD and a-ICAD but not for p-ICAD (Figure 1). When the ICAD risk was analyzed across variant and BMI categories, an increasing trend was observed only among variant carriers, with no such trend among non-carriers (P for interaction=0.014) (Table 2). These findings were consistent with the sensitivity analyses using ICAS as the outcome (Supplementary Tables 5 and 6).
Based on our findings, the RNF213 p.R4810K variant was significantly associated with an increased risk of both a-ICAD and p-ICAD. This association appeared particularly strong for a-ICAD in younger individuals and those with a higher BMI. Our results suggest that an elevated BMI, that is, obesity, may increase the risk of a-ICAD among RNF213 p.R4810K carriers (Figure 1 and Table 2). RNF213 encodes an E3 ubiquitin ligase that modulates signaling pathways involving nuclear factor of activated T cells and hypoxia-inducible factor, both of which are implicated in obesity, angiogenesis, and atherosclerosis [8,9]. Reduced RNF213 function due to the variant [10] may therefore increase susceptibility to vascular disease under conditions of metabolic stress. Although the mechanism by which these factors selectively promote a-ICAD remains unclear, regional vascular characteristics, such as differences in shear stress, may contribute. Further investigation should focus on identifying obesity-related genetic or environmental factors—such as insulin resistance, chronic inflammation, and dysregulated lipid metabolism—as potential second-hit modifiers in RNF213-related a-ICAD.
Although the RNF213 variant precedes clinical events and may support causal inference, this study was cross-sectional and conducted in a single region of Japan, which may limit the temporal and geographic generalizability of our findings. Although exploratory, the interaction analysis revealed patterns that warrant further investigation. These findings suggest that RNF213 p.R4810K carriers may be particularly susceptible to ICAD in the context of obesity, helping to clarify the heterogeneous vascular phenotypes associated with RNF213 and guiding future research on modifiable cofactors. Further investigation into obesity-related genetic and environmental factors may elucidate the mechanisms underlying RNF213-related vasculopathy and support the development of novel therapeutic strategies for RNF213 p.R4810K carriers.









