Skip to content
Frontiers
Search

Original Research

Experimental Biology and Medicine ·

Volume 251 - 2026 | https://doi.org/10.3389/ebm.2026.10948

Genetic evidence for causal relationships between brain functional networks and domain-specific recovery after nondisabling ischemic stroke

  • Huan Cai1
  • Zhenchun Huang2
  • Jialin Liang3
  • Hao Zhang4*corresponding author
  • Zhonghua Liu1*corresponding author
  1. 1Department of Rehabilitation, Zhongshan City People’s Hospital, Zhongshan, Guangdong, China
  2. 2Department of Cardiology, Shantou Central Hospital, Shantou, Guangdong, China
  3. 3Department of Endocrinology and Metabolism, Zhongshan City People’s Hospital, Zhongshan, Guangdong, China
  4. 4Department of Neurology, Affiliated Hangzhou First People’s Hospital, Westlake University School of Medicine, Hangzhou, Zhejiang, China
Article metrics0 viewsView details
PDFCiteMetrics

Abstract

Intrinsic brain networks are crucial for post-stroke recovery, but the causal relationships between specific networks and domain-specific recovery outcomes, as well as the role of lipid metabolism, remain unclear. This study leveraged Mendelian randomization (MR) to evaluate 191 resting-state functional magnetic resonance imaging (rs-fMRI) BOLD-derived phenotypes in relation to post-stroke recovery after nondisabling ischemic stroke. Genetic instruments for rs-fMRI phenotypes were derived from a UK Biobank genome-wide association study (n = 34,691). Outcomes included motor, cognitive, and global recovery after nondisabling ischemic stroke, assessed using longitudinal National Institutes of Health Stroke Scale subscales over 2 years (n = 1,270). Primary analyses used the multiplicative random-effects inverse-variance weighted method. A two-step MR analysis investigated whether brain networks mediate the effects of lipids on post-stroke outcomes. Higher BOLD-derived functional connectivity within the triple network (default mode network, central executive network, and salience network) was associated with better motor and cognitive outcomes. Higher genetically predicted orbitofrontal node amplitude in the limbic network correlated with better motor recovery, while stronger parieto-frontal connectivity was associated with cognitive recovery. Genetically proxied higher low-density lipoprotein cholesterol (LDL-C) was associated with poorer cognitive recovery, with evidence suggesting partial mediation through differences in BOLD-derived triple-network connectivity. This MR study supports a potential causal role of BOLD-derived functional network phenotypes, particularly the triple network, in motor and cognitive recovery, and further suggests that differences in triple-network connectivity act as a partial mediator linking elevated LDL-C liability to impaired cognitive recovery. These findings provide hypothesis-generating evidence for future mechanistic studies and for exploring whether specific brain network-targeted interventions could have a role in stroke recovery.

Keywords

Impact statement

While correlations between functional connectivity and specific aspects of stroke recovery have been reported, the causal nature of these associations warrants clarification. This study provides genetic evidence supporting potential causal contributions of specific intrinsic brain networks in post-stroke recovery, moving beyond prior correlational findings. It identifies genetically predicted higher BOLD-derived functional connectivity within the triple network as relevant to both motor and cognitive recovery and suggests that elevated low-density lipoprotein cholesterol was associated with poorer cognitive outcomes partially through differences in this network. These findings offer a preliminary and hypothesis generating rationale for targeting triple network connectivity as a potential biomarker and candidate intervention-relevant network feature in post-stroke recovery, and for testing whether lipid management is associated with preservation of intrinsic network function and cognitive recovery in prospective clinical studies.

Introduction

Ischemic stroke is a leading cause of death and disability worldwide, highlighting the need to understand its risk factors and post-stroke recovery mechanisms. Despite advances in understanding stroke pathogenesis, the mechanisms driving post-stroke functional recovery remain poorly understood [1].

Functional connectivity (FC) refers to temporal dependencies in neural activity across distinct brain regions and is fundamental to neural plasticity [2]. Resting-state functional magnetic resonance imaging (rs-fMRI) non-invasively maps these connectivity patterns by measuring spontaneous low-frequency blood oxygen level–dependent (BOLD) signal fluctuations, which indirectly index neural activity through the hemodynamic response and are influenced by cerebral blood flow, vascular reactivity, oxygen metabolism, and endothelial function [2]. Given the importance of neural plasticity in post-stroke recovery, FC was prioritized as a key network-level phenotype for investigating the brain’s capacity to reorganize and support recovery after stroke-related injury.

Growing evidence links altered intrinsic brain functional networks to post-stroke recovery outcomes. Disrupted FC after a stroke correlates with functional impairment, while its reestablishment facilitates recovery [3–5]. Although FC correlates with specific aspects of stroke recovery [6], such as motor function [7, 8] and cognition [9], causality remains unclear.

Mendelian randomization (MR) infers causality using genetic variants as instrumental variables (IVs). A recent bidirectional MR study did not establish robust causal links between 191 rs-fMRI phenotypes and post-stroke functional recovery [10], potentially attributable to the use of the modified Rankin Scale (mRS), which combines varied neurological impairments into a single measure and may obscure domain-specific recovery. No MR study has yet examined causal relationships between FC and domain-specific stroke recovery.

Furthermore, emerging evidence implicates a complex interplay between lipid profiles and brain networks. Imbalances in lipid metabolism may disrupt FC, especially in regions associated with cognition [11, 12], which are critical for post-stroke recovery. However, whether brain networks can mediate the effect of blood lipids on recovery remains unknown.

Thus, this study employed MR to evaluate relationships between intrinsic rs-fMRI phenotypes and domain-specific stroke recovery, as well as whether these network phenotypes mediate associations between lifelong genetic liability to lipid traits and stroke recovery.

Materials and methods

All summary-level genomic data are publicly available (Supplementary Table 1 in the electronic Supplementary Material). The study design is shown in Figure 1. All data sources were approved by relevant institutional review boards in the original studies, and all participants provided informed consent. This study follows the STROBE-MR statement [13].

FIGURE 1. Study design and three core assumptions of the MR analysis. IV, instrumental variable; MR, Mendelian randomization; SNP, single-nucleotide polymorphism; rs-fMRI, resting-state functional magnetic resonance imaging; IVW, inverse-variance weighted; cML-MA, Constrained Maximum Likelihood Mixed-Model Association; MR-PRESSO, MR-Pleiotropy Residual Sum and Outlier; BWMR, Bayesian Weighted MR.

Genome-wide association study (GWAS) of 191 rs-fMRI phenotypes

Summary GWAS data for rs-fMRI phenotypes were sourced from the UK Biobank (n = 34,691) [14]. Fluctuation amplitude was defined as the temporal standard deviation of node features within a given component. Pairwise FC was defined as temporal correlation between BOLD signal time series for each pair of nodes. Global connectivity was defined by extracting the top six principal components from the node adjacency matrix. At the genome-wide significance level, 191 phenotypes were linked to genetic variation, comprising 75 regional BOLD fluctuation amplitude traits, 111 pairwise FC traits, and 5 global FC traits (Supplementary Table 2).

GWAS of domain-specific stroke recovery

Genetic associations for post-stroke motor, cognitive, and global recovery were obtained from a GWAS of 1,270 nondisabling ischemic stroke patients, defined as modified Rankin Scale (mRS) ≤ 3, from the Vitamin Intervention for Stroke Prevention (VISP) trial with a National Institutes of Health Stroke Scale (NIHSS) score greater than zero at randomization and no recurrent stroke [15]. All participants were followed up every 3 months for 2 years.

Three recovery phenotypes were defined by changes in motor, cognitive, and global impairments at each time point. Motor impairment was measured as the change in the sum of the NIHSS motor scores for Facial Palsy, Dysarthria, Motor Arm, and Motor Leg, termed NIH Motor-6. Similarly, cognitive impairment was defined as the change in the sum of the NIHSS cognitive scores, termed NIHSS Cog-4, for Level of Consciousness (LOC) Questions, LOC Commands, Aphasia, and Inattention/Neglect. Global impairment was the change in the total NIHSS score. A linear mixed-effects model estimated the changes in NIHSS scores for each phenotype (delta = initial score − follow-up score) [15].

Genetic instrument selection

We first selected independent single-nucleotide polymorphisms (SNPs) (r *2* < 0.001, window size = 1,000 kb) with genome-wide significance (p < 5 × 10−8). SNPs with an F-statistic < 10 were excluded. We then excluded SNPs linked to confounders (e.g., depression [16], migraine [17], insomnia [16]) affecting both brain networks and post-stroke functional outcomes [18–20]. We used the GWAS Catalog database to identify and exclude these SNPs (Supplementary Table 3). Third, we applied Steiger filtering to reduce the likelihood of reverse-direction instruments, and removed networks with fewer than three IVs required for some sensitivity analyses.

Negligible sample overlap between exposure and outcome datasets was confirmed. We used proxy SNPs (r *2* > 0.8) from 1000 Genomes European reference data if needed, and removed ambiguous palindromic SNPs.

Two-sample MR analyses

Multiplicative random-effects inverse-variance weighted (IVW) MR was the primary method. Supplementary analyses included weighted median [21], MR-Egger [22], Constrained Maximum Likelihood Mixed-Model Association [23], debiased IVW [24], and Bayesian Weighted Mendelian Randomization [25]. MR-Pleiotropy Residual Sum and Outlier (MR-PRESSO) was used to identify potential outliers [26]. These methods collectively enhance robustness to IV assumption violations.

A causal relationship was inferred if all MR methods showed directionally consistent estimates and the IVW p < 8.73 × 10−5 (Bonferroni correction for 191 phenotypes and 3 outcomes).

Sensitivity analyses

We tested SNP heterogeneity with Cochran’s Q statistic and used leave-one-out plots to check if single SNPs drove the results. In addition, we repeated our analysis using a relaxed IV selection threshold of p < 1 × 10−6 to check for consistency in the direction with those obtained from the genome-wide significant threshold. To assess the possibility of reverse causality, we performed reverse MR analyses in significant rs-fMRI phenotypes from forward MR using IVs with a suggestive significance threshold (p < 5 × 10−6). MR-PRESSO global test and MR-Egger intercept were used to detect horizontal pleiotropy.

To address collider bias risk when analyses are conditioned on ischemic stroke survivors, we repeated MR analyses using genetic effects on recovery outcomes adjusted via corrected weighted least squares (CWLS) [27], based on summary ischemic stroke data from the GIGASTROKE consortium of European ancestry [28].

Mediation MR

We used two-step MR to examine whether genetically proxied lipid traits were linked to post-stroke outcomes through BOLD-derived rs-fMRI network phenotypes. Summary statistics of three blood lipid measures, including triglycerides, low-density lipoprotein cholesterol (LDL-C), and high-density lipoprotein cholesterol (HDL-C), were obtained from the most recent circulating metabolite GWAS analysis of 136,016 participants [29]. We first estimated the causal effect of lipid traits on intrinsic functional networks (coefficient βXM) and on post-stroke recovery patterns (coefficient βXY), then assessed the causal effect of rs-fMRI phenotypes on specific subdomains of post-ischemic stroke recovery (coefficient βMY). We calculated the mediation effect of BOLD-derived network phenotypes using the product of coefficients method [30]. The mediation effect (βXM × βMY) was then divided by the total effect (βXY) to calculate the proportion of the mediation effect. Standard errors for the mediation effect were calculated with the delta method. Notably, only rs-fMRI phenotypes showing robust evidence in the primary univariable MR were taken forward to mediation analyses.

Lipid and rs-fMRI traits were measured continuously and standardized to standard deviations. Analyses were performed in R (v4.2.1) with TwoSampleMR (v0.6.8) [31], RMediation (v1.2.2) [32], and MRPRESSO (v1.0) [26] packages.

Results

Genetically predicted rs-fMRI phenotypes and specific domains of post-stroke recovery

The full primary MR results are shown in Supplementary Tables 4–6. Twenty-four rs-fMRI phenotypes were nominally associated with domain-specific stroke recovery (p < 0.05, Figure 2). We identified four BOLD-derived brain functional network phenotypes showing robust MR evidence of association with two distinct recovery phenotypes (motor and cognitive recovery) after Bonferroni correction in the primary IVW analysis (Figure 3). The genetic variants used as instruments for significant intrinsic brain activity phenotypes are presented in Supplementary Table 7. However, no significant association survived multiple testing correction for post-stroke global recovery.

FIGURE 2. Circular heatmap for the IVW MR estimates of associations between rs-fMRI phenotypes and domain-specific stroke recovery.

FIGURE 3. Causal relationships that remained significant after Bonferroni correction in MR analysis results. Left: forest plot illustrating the significant causal estimates using six MR methods. Right: pattern diagram showing the brain anatomical regions of the corresponding rs-fMRI phenotypes. CI, confidence interval.

Two rs-fMRI phenotypes showed genetically supported evidence of association with motor improvement (Figure 3). Pheno593 reflects connectivity in the precuneus, angular gyrus, or cingulate, and parietal regions. The BOLD-derived signals of these brain regions involve FC in the triple network (default mode network (DMN)/central executive network (CEN) and CEN/salience network (SN)). Higher genetically predicted FC in these networks was associated with greater motor improvement (IVW β = 0.399, 95% confidence interval [CI] = 0.302–0.496, p = 8.33 × 10−16). In the limbic network, increased BOLD fluctuation amplitude of the orbitofrontal region was linked to a better motor outcome (IVW β = 0.241, 95% CI = 0.150–0.331, p = 7.78 × 10−6).

Three rs-fMRI phenotypes showed evidence for a potential causal contribution to post-stroke cognitive recovery (Figure 3). Pheno674 involves FC between parietal and frontal regions. Higher genetically predicted CEN-SN connectivity was associated with a better cognitive outcome (IVW β = 0.245, 95% CI = 0.176–0.315, p = 5.34 × 10−12). For Pheno593, connectivity in the precuneus, angular gyrus, or cingulate area, and parietal area positively correlated with cognitive recovery (IVW β = 0.232, 95% CI = 0.138–0.327, p = 1.52 × 10−6). For Pheno1276, resting frontal BOLD-derived signal features were linked to cognitive changes after stroke, with higher FC between CEN and SN or DMN or CEN associated with greater cognitive improvement (IVW β = 0.188, 95% CI = 0.106–0.270, p = 7.78 × 10−6).

These findings were supported by other MR analyses (Figure 3; Supplementary Tables 8,9). MR-PRESSO detected no significant outliers in exposure-outcome pairs for which the test was applicable.

Sensitivity analysis

Cochran’s Q test showed no statistical evidence of heterogeneity among individual SNPs (Table 1). The scatterplots, forest plots, and leave-one-out plots suggested that no single SNP drove the associations (Supplementary Figures 1–6 in the electronic Supplementary Material). No statistical evidence of horizontal pleiotropy was detected by MR-PRESSO global test and MR-Egger intercept (all p > 0.05, Table 1). Results using the relaxed threshold were still statistically significant and directionally consistent with the primary analysis (Supplementary Table 10). Reverse MR analyses (Supplementary Tables 11, 12) showed no evidence of reverse causality. Associations remained stable after adjusting for collider bias (Supplementary Tables 8,9).

TABLE 1. Pleiotropy and heterogeneity tests for significant exposure-outcome pairs in forward Mendelian randomization.

ExposureOutcomePleiotropyHeterogeneity (IVW)
MR-PRESSO global p valueMR-Egger interceptp valueCochran’s Qp value
Pheno593Post-stroke motor recoveryNA−0.0180.8510.0570.972
Pheno71Post-stroke motor recovery0.9940.0180.9110.0640.996
Pheno593Post-stroke cognitive recovery0.909−0.0090.7760.4790.923
Pheno1276Post-stroke cognitive recovery0.974−0.0080.8760.2260.973
Pheno674Post-stroke cognitive recoveryNA−0.0220.7940.1200.942

Abbreviations: MR-PRESSO, Mendelian randomization-pleiotropy residual sum and outlier; IVW, inverse-variance weighted; NA, not applicable.

Mediation analysis

We found that genetically predicted higher LDL-C levels, but not triglycerides or HDL-C, were associated with poorer cognitive recovery (IVW β = −0.057; 95% CI = −0.105, −0.008; p = 0.022, Supplementary Figures 7,8). We found no MR evidence supporting associations between lipid traits and motor recovery. We then conducted MR analysis of LDL-C and its subfractions, including total cholesterol in large LDL, total cholesterol in medium LDL (M-LDL-C), and total cholesterol in small LDL (S-LDL-C), to examine whether they were associated with three imaging phenotypes that showed a potential causal contribution to post-stroke cognitive recovery. We found that LDL-C and its subfractions were associated with Pheno593, but not Pheno674 or Pheno1276 (Supplementary Figures 8,9). Finally, for Pheno593, our exploratory mediation analysis suggested that BOLD-derived connectivity in the precuneus, angular gyrus, or cingulate, and parietal regions, as well as the triple network (DMN/CEN and CEN/SN), may partially mediate the associations between genetic liability to LDL-C, M-LDL-C, and S-LDL-C and cognitive recovery (Figure 4). Specifically, the mediation effects were estimated to be −0.008 (95% CI -0.016, −0.001; p = 0.033) for LDL-C, −0.008 (95% CI -0.016, −0.001; p = 0.032) for M-LDL-C, and −0.009 (95% CI -0.017, −0.001; p = 0.025) for S-LDL-C, with mediation proportions of 14.46% (95% CI 0.96%–27.96%), 14.60% (95% CI 1.01%–28.18%), and 15.09% (95% CI 1.70%–28.48%), respectively. No statistical evidence of horizontal pleiotropy or heterogeneity was detected, and the Steiger directionality tests supported the inferred causal direction (Supplementary Table 13).

FIGURE 4. Associations between genetically predicted LDL-C, M-LDL-C, S-LDL-C and post-stroke cognitive recovery mediated by brain networks. LDL-C: low-density lipoprotein cholesterol; M-LDL-C, total cholesterol in medium LDL; S-LDL-C, total cholesterol in small LDL.

Discussion

Using MR, we found that specific rs-fMRI phenotypes were associated with motor and cognitive recovery, and that BOLD-derived connectivity may partly mediate associations between genetic liability to lipid traits and cognitive recovery. These findings may help inform future research on connectivity profiles related to poor recovery.

Our results support a genetically informed association between triple-network (DMN-CEN-SN) connectivity and both motor and cognitive recovery after stroke. Higher genetically predicted BOLD-derived connectivity within this large-scale network system was associated with better recovery outcomes (Pheno593), suggesting that post-stroke recovery may depend not only on restoration of domain-specific circuits, but also on the integrity of domain-general networks that coordinate attention, executive control, self-monitoring, and adaptive learning [33, 34]. Within the triple-network framework, the SN, anchored in the anterior insula and dorsal anterior cingulate cortex, detects behaviorally relevant stimuli and facilitates dynamic switching between the DMN and the CEN. The CEN supports goal-directed control, working memory, performance monitoring, and error correction, whereas the DMN contributes to memory-based prediction, self-referential processing, and internally guided cognition [33]. These processes are directly relevant to cognitive recovery after stroke, as impairments in attention allocation, executive control, memory retrieval, and self-monitoring are common contributors to post-stroke cognitive dysfunction [33]. Consistent with this view, previous studies have shown that alterations in triple-network organization are associated with cognitive impairment [35], and that preservation or restoration of FC within and between regions of the DMN is related to better cognitive outcomes after stroke [6]. Importantly, the same network architecture may also support motor recovery. Although motor restoration is strongly influenced by the integrity of sensorimotor pathways and interhemispheric motor connectivity [6], successful motor rehabilitation also requires attention to task-relevant cues, selection of compensatory strategies, feedback-based error correction, and motor relearning [36]. These functions depend heavily on interactions between the SN and CEN and may interact with DMN-related memory and self-monitoring processes [37]. Thus, more favorable triple-network connectivity may provide a domain-general network-level scaffold through which cognitive-control systems support adaptive motor performance, particularly when primary sensorimotor pathways are disrupted [37], although the rs-fMRI phenotype may also partly index vascular support for network function.

This concurrent association with motor and cognitive recovery may reflect large-scale network integration, shared plasticity mechanisms, and vascular support for metabolically demanding recovery processes. Stroke recovery is thought to rely on experience-dependent plasticity, including unmasking of latent synaptic connections, long-term potentiation-like mechanisms, and functional reorganization of spared cortical and subcortical circuits [1]. The hubs of the SN and CEN are anatomically and functionally connected with association cortices and secondary motor regions, including the premotor cortex and supplementary motor area, through cortico-cortical and cortico-striatal-thalamic pathways [38, 39]. These pathways may serve as convergence zones where salience detection, executive control, memory-based prediction, and motor planning are integrated during rehabilitation-related learning. Therefore, stronger genetically predicted triple-network connectivity may reflect an inherited propensity for more efficient large-scale network communication or more favorable neurovascular coupling and vascular reactivity, which could facilitate plastic reorganization and coordinated recovery across both cognitive and motor domains.

Within the limbic network, we found that genetically predicted higher node amplitude in the orbitofrontal region was linked to better motor outcomes (Pheno71). The limbic network’s involvement in motor recovery may suggest an interplay between emotional/motivational factors and physical rehabilitation, expanding traditional models focused solely on sensorimotor pathways.

The parieto-frontal network is integral to various cognitive processes, including attention, working memory, and executive functions [40]. Our MR analysis showed that the parieto-frontal network (Pheno674) exhibited a positive association with better cognitive outcomes, consistent with the previous findings of reduced global and local efficiency in the parieto-frontal network in subjects with cognitive impairment [41].

In conjunction with a prior study [12], our results provide genetic evidence consistent with the possibility that FC mediates the relationship between atherogenic lipid profiles and cognition, and extend this hypothesis to stroke recovery. While LDL-C is a well-established risk factor for ischemic stroke, its role in recovery is less understood, particularly because post-stroke LDL-C levels are influenced by treatment and clinical factors. Prior research indicated that higher LDL-C levels correlated with greater cognitive impairment in patients with cerebral small vessel disease, suggesting that LDL-C may contribute to cognitive decline through vascular mechanisms [42]. Our findings suggest that genetic liability to atherogenic lipid (LDL-C, M-LDL-C, S-LDL-C) traits may be associated with poorer cognitive recovery after stroke partly through genetically proxied differences in BOLD-derived triple-network connectivity, possibly through effects on neuroinflammation [43, 44], vascular integrity, endothelial function, cerebral perfusion, neurovascular coupling, amyloidogenic factors [45], neural repair mechanisms, or shared genetic influences. Future studies should test whether optimized lipid management can preserve brain networks or mitigate cognitive deficits after stroke.

Several limitations should be considered. First, because the VISP outcome cohort consisted predominantly of patients with nondisabling ischemic stroke of European ancestry, the modest sample limits generalizability to more clinically heterogeneous populations and to other ancestry groups. Second, resting-state BOLD response reflects spontaneous low-frequency neural activity but is also influenced by neurovascular coupling, vascular reactivity, and metabolic integrity. Therefore, this phenotype may capture variation in both functional network organization and the vascular–metabolic support systems. In the context of lipid metabolism, this is particularly important because atherogenic lipid profiles may affect endothelial function, cerebral perfusion, and vascular integrity, which could in turn alter rs-fMRI connectivity estimates independently of, or in combination with, neural communication. Thus, mediation by rs-fMRI phenotypes should be interpreted as reflecting BOLD-derived neurovascular or hemodynamic pathways, not purely neural connectivity. Third, recovery phenotypes based on NIHSS subscales are relatively coarse compared to detailed assessments. Fourth, MR estimates reflect the lifelong genetic liability to BOLD-derived brain functional phenotypes, with genetic instruments for the 191 rs-fMRI traits derived from healthy UK Biobank participants. Therefore, these findings do not prove that modifying FC after stroke would improve recovery, especially given lesion-related reorganization and disconnection. Fifth, our summary-level MR design could not account for lesion localization. Lesions situated within or adjacent to the key networks examined here may themselves alter the measured FC and could modulate the observed network–recovery relationships. Because the outcome GWAS did not provide lesion-location information, we were unable to stratify or adjust for this anatomical confounder. Finally, limited sample size of the outcome GWAS and instrument availability may reduce power and precision of MR estimates, and increase susceptibility to winner’s curse or overestimation. Although Steiger directionality tests and pleiotropy sensitivity analyses supported the robustness of the mediation-related associations, shared genetic architecture and unmeasured pleiotropic pathways cannot be ruled out. The reported effect sizes for the detected associations should be interpreted as preliminary rather than definitive estimates of magnitude, especially for the mediation analysis, because the lipid instruments represent lifelong genetic predisposition rather than longitudinally measured or treated LDL-C levels. Despite these limitations, our findings suggest a biologically plausible and testable pathway linking lipid genetic liability, brain network organization, and recovery.

Conclusion

This MR study provides preliminary genetic evidence supporting potential causal contributions of four BOLD-derived rs-fMRI phenotypes to motor and cognitive recovery after nondisabling ischemic stroke, highlighting the possible relevance of the triple network (DMN, CEN, SN). Additionally, BOLD-derived connectivity features involving DMN, CEN, and SN networks may in part mediate the association between genetically proxied lipid traits and cognitive outcomes. However, this potential pathway should be interpreted cautiously and requires validation in additional cohorts and larger GWAS datasets with more available instruments, ideally including detailed lesion localization, serial post-stroke lipid measurements, medication data, and post-stroke rs-fMRI measures. Future research could use these hypothesis-generating findings to identify candidate network targets for mechanistic studies and to evaluate whether interventions such as noninvasive brain stimulation can modulate relevant connectivity patterns and improve stroke recovery.

Data availability

Publicly available datasets were analyzed in this study. The summary level data for rs-fMRI phenotypes are available at https://doi.org/10.5281/zenodo.5775047. The summary level data for post-stroke recovery phenotypes and GIGASTROKE data are available at the Cerebrovascular Disease Knowledge Portal. The summary level data for blood lipid measures are available through the NHGRI-EBI GWAS Catalogue.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.ebm-journal.org/articles/10.3389/ebm.2026.10948/full#supplementary-material

Supplementary material

Supplementary material

Statements

Conflict of interest

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This study was supported by the Social Public Welfare and Basic Research Project of Zhongshan City (No. 2023B3025) and the Research Project of the Traditional Chinese Medicine Bureau of Guangdong Province (No. 20251440).

Ethics approval

The studies involving humans were approved by relevant institutional review boards in the original studies. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

HC: Writing – original draft, Supervision, Conceptualization. ZH: Supervision, Writing – review and editing, Software. JL: Methodology, Formal Analysis, Data curation. HZ: Writing – review and editing, Validation, Supervision. ZL: Writing – review and editing, Supervision, Methodology, Investigation, Data curation. All authors contributed to the article and approved the submitted version.

Acknowledgements

We gratefully thank the UK Biobank, International Stroke Genetics Consortium (ISGC) Cerebrovascular Disease Knowledge Portal, and GIGASTROKE consortium for providing summary statistics data.

References

  1. 1. Grefkes C, Fink GR. Recovery from stroke: current concepts and future perspectives. Neurol Res Pract (2020) 2: 17. doi:10.1186/s42466-020-00060-6
  2. 2. van den Heuvel MP, Hulshoff Pol HE. Exploring the brain network: a review on resting-state fMRI functional connectivity. Eur Neuropsychopharmacol (2010) 20: 519 – 34. doi:10.1016/j.euroneuro.2010.03.008
  3. 3. van Meer MP, van der Marel K, Wang K, Otte WM, El Bouazati S, Roeling TA Recovery of sensorimotor function after experimental stroke correlates with restoration of resting-state interhemispheric functional connectivity. J Neurosci (2010) 30: 3964 – 72. doi:10.1523/JNEUROSCI.5709-09.2010
  4. 4. Arun KM, Smitha KA, Sylaja PN, Kesavadas C. Identifying resting-state functional connectivity changes in the motor cortex using fNIRS during recovery from stroke. Brain Topogr (2020) 33: 710 – 9. doi:10.1007/s10548-020-00785-2
  5. 5. Golestani AM, Tymchuk S, Demchuk A, Goodyear BG, Group V-S. Longitudinal evaluation of resting-state FMRI after acute stroke with hemiparesis. Neurorehabil Neural Repair (2013) 27: 153 – 63. doi:10.1177/1545968312457827
  6. 6. Ismail UN, Yahya N, Manan HA. Investigating functional connectivity related to stroke recovery: a systematic review. Brain Res (2024) 1840: 149023. doi:10.1016/j.brainres.2024.149023
  7. 7. Park CH, Chang WH, Ohn SH, Kim ST, Bang OY, Pascual-Leone A Longitudinal changes of resting-state functional connectivity during motor recovery after stroke. Stroke (2011) 42: 1357 – 62. doi:10.1161/STROKEAHA.110.596155
  8. 8. Wu CW, Lin SN, Hsu LM, Yeh SC, Guu SF, Lee SH Synchrony between default-mode and sensorimotor networks facilitates motor function in stroke rehabilitation: a pilot fMRI study. Front Neurosci (2020) 14: 548. doi:10.3389/fnins.2020.00548
  9. 9. Vicentini JE, Weiler M, Casseb RF, Almeida SR, Valler L, de Campos BM Subacute functional connectivity correlates with cognitive recovery six months after stroke. Neuroimage Clin (2021) 29: 102538. doi:10.1016/j.nicl.2020.102538
  10. 10. Li MZ, Shi YL, He XJ, Wang SC, Liu J, Wang Z Resting-state fMRI and post-ischemic stroke functional recovery: unraveling causality and predicting therapeutic targets. Int J Mol Sci (2025) 26: 3608. doi:10.3390/ijms26083608
  11. 11. Xia W, Zhang B, Yang Y, Wang P, Yang Y, Wang S. Poorly controlled cholesterol is associated with cognitive impairment in T2DM: a resting-state fMRI study. Lipids Health Dis (2015) 14: 47. doi:10.1186/s12944-015-0046-x
  12. 12. Zhang W, Fu L, Bi Y, Liu J, Li X, Zhang X Large-scale functional network connectivity mediates the associations between lipids metabolism and cognition in type 2 diabetes. J Cereb Blood Flow Metab (2024) 44: 384 – 96. doi:10.1177/0271678X231204426
  13. 13. Skrivankova VW, Richmond RC, Woolf BAR, Yarmolinsky J, Davies NM, Swanson SA Strengthening the reporting of observational studies in epidemiology using mendelian randomization: the STROBE-MR statement. JAMA (2021) 326: 1614 – 21. doi:10.1001/jama.2021.18236
  14. 14. Zhao B, Li T, Smith SM, Xiong D, Wang X, Yang Y Common variants contribute to intrinsic human brain functional networks. Nat Genet (2022) 54: 508 – 17. doi:10.1038/s41588-022-01039-6
  15. 15. Aldridge CM, Braun R, Lohse K, de Havenon A, Cole JW, Cramer SC Genome-wide association studies of 3 distinct recovery phenotypes in mild ischemic stroke. Neurology (2024) 102: e208011. doi:10.1212/WNL.0000000000208011
  16. 16. Mu C, Dang X, Luo XJ. Mendelian randomization analyses reveal causal relationships between brain functional networks and risk of psychiatric disorders. Nat Hum Behav (2024) 8: 1417 – 28. doi:10.1038/s41562-024-01879-8
  17. 17. Lu Q, Jia Z, Gu H. Association between brain resting-state functional activities and migraine: a bidirectional mendelian randomization study. Sci Rep (2024) 14: 23901. doi:10.1038/s41598-024-74745-2
  18. 18. Zhang Z, Wang M, Gill D, Liu X, Zhu W. Association of genetically predicted anxiety and depression with functional outcome after ischemic stroke: a mendelian randomization study. Neurology (2024) 103: e209776. doi:10.1212/WNL.0000000000209776
  19. 19. Zhang Z, Wang M, Gill D, Zhu W, Liu X. Genetically predicted sleep traits and functional outcome after ischemic stroke: a mendelian randomization study. Neurology (2023) 100: e1159 – e65. doi:10.1212/WNL.0000000000206745
  20. 20. Wang M, Daghlas I, Zhang Z, Ye D, Li S, Liu D. Genetic liability to migraine and functional outcome after ischemic stroke. Eur Stroke J (2023) 8: 517 – 21. doi:10.1177/23969873231164728
  21. 21. Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol (2016) 40: 304 – 14. doi:10.1002/gepi.21965
  22. 22. Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through egger regression. Int J Epidemiol (2015) 44: 512 – 25. doi:10.1093/ije/dyv080
  23. 23. Xue H, Shen X, Pan W. Constrained maximum likelihood-based mendelian randomization robust to both correlated and uncorrelated pleiotropic effects. Am J Hum Genet (2021) 108: 1251 – 69. doi:10.1016/j.ajhg.2021.05.014
  24. 24. Su Y, Xu S, Ma Y, Yin P, Hao X, Zhou J A modified debiased inverse-variance weighted estimator in two-sample summary-data Mendelian randomization. Stat Med (2024) 43: 5484 – 96. doi:10.1002/sim.10245
  25. 25. Zhao J, Ming J, Hu X, Chen G, Liu J, Yang C. Bayesian weighted Mendelian randomization for causal inference based on summary statistics. Bioinformatics (2020) 36: 1501 – 8. doi:10.1093/bioinformatics/btz749
  26. 26. Verbanck M, Chen CY, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet (2018) 50: 693 – 8. doi:10.1038/s41588-018-0099-7
  27. 27. Cai S, Hartley A, Mahmoud O, Tilling K, Dudbridge F. Adjusting for collider bias in genetic association studies using instrumental variable methods. Genet Epidemiol (2022) 46: 303 – 16. doi:10.1002/gepi.22455
  28. 28. Mishra A, Malik R, Hachiya T, Jurgenson T, Namba S, Posner DC Stroke genetics informs drug discovery and risk prediction across ancestries. Nature (2022) 611: 115 – 23. doi:10.1038/s41586-022-05165-3
  29. 29. Karjalainen MK, Karthikeyan S, Oliver-Williams C, Sliz E, Allara E, Fung WT Genome-wide characterization of circulating metabolic biomarkers. Nature (2024) 628: 130 – 8. doi:10.1038/s41586-024-07148-y
  30. 30. Carter AR, Sanderson E, Hammerton G, Richmond RC, Davey Smith G, Heron J Mendelian randomisation for mediation analysis: current methods and challenges for implementation. Eur J Epidemiol (2021) 36: 465 – 78. doi:10.1007/s10654-021-00757-1
  31. 31. Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D The MR-Base platform supports systematic causal inference across the human phenome. Elife (2018) 7: e34408. doi:10.7554/eLife.34408
  32. 32. Tofighi D, MacKinnon DP. RMediation: an R package for mediation analysis confidence intervals. Behav Res Methods (2011) 43: 692 – 700. doi:10.3758/s13428-011-0076-x
  33. 33. Menon V. Large-scale brain networks and psychopathology: a unifying triple network model. Trends Cogn Sci (2011) 15: 483 – 506. doi:10.1016/j.tics.2011.08.003
  34. 34. Siegel JS, Ramsey LE, Snyder AZ, Metcalf NV, Chacko RV, Weinberger K Disruptions of network connectivity predict impairment in multiple behavioral domains after stroke. Proc Natl Acad Sci U S A (2016) 113: E4367 – 76. doi:10.1073/pnas.1521083113
  35. 35. Wei HL, Wei C, Yu YS, Yu X, Chen Y, Li J Dysfunction of the triple-network model is associated with cognitive impairment in patients with cerebral small vessel disease. Heliyon (2024) 10: e24701. doi:10.1016/j.heliyon.2024.e24701
  36. 36. Kwakkel G, Stinear C, Essers B, Munoz-Novoa M, Branscheidt M, Cabanas-Valdes R Motor rehabilitation after stroke: European stroke organisation (ESO) consensus-based definition and guiding framework. Eur Stroke J (2023) 8: 880 – 94. doi:10.1177/23969873231191304
  37. 37. Lee J, Kim YH. Does a cognitive network contribute to motor recovery after ischemic stroke? Neurorehabil Neural Repair (2023) 37: 458 – 65. doi:10.1177/15459683231177604
  38. 38. Menon V, Uddin LQ. Saliency, switching, attention and control: a network model of insula function. Brain Struct Funct (2010) 214: 655 – 67. doi:10.1007/s00429-010-0262-0
  39. 39. Peters SK, Dunlop K, Downar J. Cortico-striatal-thalamic loop circuits of the salience network: a central pathway in psychiatric disease and treatment. Front Syst Neurosci (2016) 10: 104. doi:10.3389/fnsys.2016.00104
  40. 40. Frontzkowski L, Fehring F, Frey BM, Wrobel PP, Reibelt A, Higgen F Frontoparietal structural network disconnections correlate with outcome after a severe stroke. Hum Brain Mapp (2024) 45: e70060. doi:10.1002/hbm.70060
  41. 41. Chen H, Huang L, Yang D, Ye Q, Guo M, Qin R Nodal global efficiency in front-parietal lobe mediated periventricular white matter hyperintensity (PWMH)-related cognitive impairment. Front Aging Neurosci (2019) 11: 347. doi:10.3389/fnagi.2019.00347
  42. 42. Cheng Y, Li L, Lv Y, Zhang L, Chen W, Xu G. Association between cerebral small vessel disease and plasma levels of LDL cholesterol and homocysteine: implications for cognitive function. J Med Biochem (2024) 43: 696 – 703. doi:10.5937/jomb0-50100
  43. 43. Sohrabi Y, Schwarz D, Reinecke H. LDL-C augments whereas HDL-C prevents inflammatory innate immune memory. Trends Mol Med (2022) 28: 1 – 4. doi:10.1016/j.molmed.2021.11.003
  44. 44. Izumida T, Nakamura Y, Hino Y, Ishikawa S. Combined effect of small dense low-density lipoprotein cholesterol (sdLDL-C) and remnant-like particle cholesterol (RLP-C) on low-grade inflammation. J Atheroscler Thromb (2020) 27: 319 – 30. doi:10.5551/jat.49528
  45. 45. Wingo AP, Vattathil SM, Liu J, Fan W, Cutler DJ, Levey AI LDL cholesterol is associated with higher AD neuropathology burden independent of APOE. J Neurol Neurosurg Psychiatry (2022) 93: 930 – 8. doi:10.1136/jnnp-2021-328164

Citation

Cai, H., et al. (2026). Genetic evidence for causal relationships between brain functional networks and domain-specific recovery after nondisabling ischemic stroke. Experimental Biology and Medicine, Volume 251. https://doi.org/
Received
Accepted
Published
Volume
251 - 2026

Copyright

Copyright © 2026 Cai, Huang, Liang, Zhang and Liu.

License: https://creativecommons.org/licenses/by/4.0/

*Correspondence: Hao Zhang; Zhonghua Liu

Article metrics

Views
0
Downloads
0
Citations
0
Shares
0

Related articles

ReviewApplication of natural products combined with plant-derived extracellular vesicles in cancer therapyRuyue Wang, Menglin Ren, Yuanjian Xu, Yan Sun, Zenghui XuOriginal ResearchRefined protocol for quantifying hippocampal neurogenesis in organotypic slice culturesIrem Akyel, Graham K. SheridanReviewMitochondrial metabolic plasticity in cancer: hybrid bioenergetic states and translational targeting in surgical oncologyDebdeep Chattopadhyay, Ramya Lakshmi Rajendran, Nida Mubin, Durva Tendulkar, Thoudam Debraj Singh, Subhrojyoti Ghosh, Saurabh Kumar Jha, Prakash Gangadaran, Byeong-Cheol Ahn