
==== Front
Brain Commun
Brain Commun
braincomms
Brain Communications
2632-1297
Oxford University Press UK

10.1093/braincomms/fcae129
fcae129
Original Article
AcademicSubjects/MED00310
AcademicSubjects/SCI01870
Bayesian modelling disentangles language versus executive control disruption in stroke
https://orcid.org/0000-0002-8084-1330
Hartwigsen Gesa Wilhelm Wundt Institute for Psychology, Leipzig University, 04109 Leipzig, Germany
Research Group Cognition and Plasticity, Max Planck Institute for Human Cognitive and Brain Sciences, 04103 Leipzig, Germany

Lim Jae-Sung Department of Neurology, Asan Medical Center, University of Ulsan College of Medicine, Seoul, 05505, South Korea

Bae Hee-Joon Department of Neurology, Seoul National University Bundang Hospital, Seoul National University College of Medicine, Seoul, 13620, South Korea

Yu Kyung-Ho Department of Neurology, Hallym University Sacred Heart Hospital, Hallym University College of Medicine, Anyang, 14068, Republic of Korea

Kuijf Hugo J Image Sciences Institute, University Medical Center Utrecht, 3584 CX Utrecht, The Netherlands

Weaver Nick A Department of Neurology and Neurosurgery, Utrecht Brain Center, University Medical Center Utrecht, 3584 CX Utrecht, The Netherlands

https://orcid.org/0000-0001-7017-2148
Biesbroek J Matthijs Department of Neurology and Neurosurgery, Utrecht Brain Center, University Medical Center Utrecht, 3584 CX Utrecht, The Netherlands
Department of Neurology, Diakonessenhuis Hospital, 3582 KE Utrecht, The Netherlands

Kopal Jakub Department of Biomedical Engineering, Faculty of Medicine, McConnell Brain Imaging Centre, Montreal Neurological Institute, McGill University, Montreal, Quebec H3A 2BA, Canada
Mila—Quebec Artificial Intelligence Institute, Montreal, Quebec H2S 3H1, Canada

Bzdok Danilo Department of Biomedical Engineering, Faculty of Medicine, McConnell Brain Imaging Centre, Montreal Neurological Institute, McGill University, Montreal, Quebec H3A 2BA, Canada
Mila—Quebec Artificial Intelligence Institute, Montreal, Quebec H2S 3H1, Canada

Correspondence to: Gesa Hartwigsen Wilhelm Wundt Institute for Psychology, Leipzig University, Neumarkt 9-19, 04109 Leipzig, Saxony, Germany E-mail: hartwigsen@cbs.mpg.de
Correspondence may also be addressed to: Jakub Kopal Department of Biomedical Engineering, Faculty of Medicine, McConnell Brain Imaging Centre, Montreal Neurological Institute McGill University, 3775 Rue University, Montréal, QC, H3A 2B4 Canada E-mail: kopal.kuba@gmail.com
Correspondence may also be addressed to: Danilo Bzdok E-mail: danilo.bzdok@mcgill.ca
Jakub Kopal and Danilo Bzdok contributed equally to this work.

2024
10 4 2024
10 4 2024
6 3 fcae12917 10 2023
06 2 2024
09 4 2024
04 5 2024
© The Author(s) 2024. Published by Oxford University Press on behalf of the Guarantors of Brain.
2024
https://creativecommons.org/licenses/by/4.0/ This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Stroke is the leading cause of long-term disability worldwide. Incurred brain damage can disrupt cognition, often with persisting deficits in language and executive capacities. Yet, despite their clinical relevance, the commonalities and differences between language versus executive control impairments remain under-specified. To fill this gap, we tailored a Bayesian hierarchical modelling solution in a largest-of-its-kind cohort (1080 patients with stroke) to deconvolve language and executive control with respect to the stroke topology. Cognitive function was assessed with a rich neuropsychological test battery including global cognitive function (tested with the Mini-Mental State Exam), language (assessed with a picture naming task), executive speech function (tested with verbal fluency tasks), executive control functions (Trail Making Test and Digit Symbol Coding Task), visuospatial functioning (Rey Complex Figure), as well as verbal learning and memory function (Soul Verbal Learning). Bayesian modelling predicted interindividual differences in eight cognitive outcome scores three months after stroke based on specific tissue lesion topologies. A multivariate factor analysis extracted four distinct cognitive factors that distinguish left- and right-hemispheric contributions to ischaemic tissue lesions. These factors were labelled according to the neuropsychological tests that had the strongest factor loadings: One factor delineated language and general cognitive performance and was mainly associated with damage to left-hemispheric brain regions in the frontal and temporal cortex. A factor for executive control summarized mental flexibility, task switching and visual-constructional abilities. This factor was strongly related to right-hemispheric brain damage of posterior regions in the occipital cortex. The interplay of language and executive control was reflected in two distinct factors that were labelled as executive speech functions and verbal memory. Impairments on both factors were mainly linked to left-hemispheric lesions. These findings shed light onto the causal implications of hemispheric specialization for cognition; and make steps towards subgroup-specific treatment protocols after stroke.

Hartwigsen et al. report different cognitive factors that underlie cognitive impairments after stroke in a large patient cohort. Combining Bayesian modelling with comprehensive behavioural testing revealed four cognitive factors that differentiated left- and right-hemispheric contributions to brain lesions. These findings elucidate the causal implications of hemispheric specialization for cognition.

Graphical Abstract

Graphical Abstract

speech
control
lateralization
domain-general
lesion
Brain Canada Foundation 10.13039/100009408 Health Canada, National Institutes of Health NIH R01 AG068563A NIH R01 R01DA053301-01A1 Canadian Institute of Health Research 10.13039/501100000024 CIHR 438531 CIHR 470425 Healthy Brains Healthy Lives initiative CIFAR Artificial Intelligence Chairs program Lise Meitner Excellence funding from the Max Planck Society European Research Council 10.13039/100010663 ERC-2021-COG 101043747 German Research Foundation 10.13039/501100001659 HA 6314/3-1 HA 6314/4-2 HA 6314/9-1 Open Access Publishing Fund of Leipzig University
==== Body
pmcIntroduction

In our rapidly aging societies, stroke is now the leading cause of long-term disability worldwide, with 12.2 million new cases each year.1 Globally, one in four people will be affected by stroke in their lifetime. Stroke often severely affects cognition and can cause loss of language and executive functions.2-4 These cognitive faculties are crucial for interpersonal interaction in everyday life. Language is a key ability for communication uniquely developed in humans, including production and comprehension abilities.5 Aside from specific linguistic operations, efficient communication also requires controlled planning, focusing and flexible thinking. These mental skills are subsumed as executive functions and include inhibitory control, working memory and cognitive flexibility.6 Brain damage can severely affect executive functions, including the ability to complete basic and complex activities of daily living and participate in work, social and leisure activities.7 Yet, despite the wide-ranging impact of stroke-induced dysfunction on the individual patient’s cognitive abilities and our society at large, it remains unclear how the neural networks for cognitive functions recover from tissue damage.8 A deeper look into the commonalities and differences of stroke-induced tissue dysfunctions leading to language versus executive control deficits would identify biologically valid subgroups, thereby paving the way for more accurate outcome predictions and better-targeted therapeutics in the future of precision medicine.9,10

Indeed, today’s cognitive therapy after stroke routinely ignores the individual topography of the original brain lesion.11 Such one-fits-all treatment is probably a culprit for the observation that there is substantial variability across patients in treatment success, leading to overall small effect sizes.12,13 The critical gap in our knowledge about the specificity of stroke-induced impairments on different cognitive functions may partly result from the fact that most previous studies focused on impairments in single cognitive domains derived from small patient cohorts.14-16 Such studies suggest that language impairments result from damage to key language regions in the left hemisphere, while executive deficits have been preferentially linked to lesions in the right hemisphere.8,16 This distinction ignores the functional interplay of both domains, especially under challenging conditions. Indeed, recent work provides clear clues that successful language recovery after stroke may extend beyond classically studied language regions: neuroimaging studies suggest that domain-specific recovery of language functions also includes recruitment of domain-general networks for cognitive control, such as the multiple-demand network.17-20 This implies that some patients with deficits in the language domain may particularly benefit from training in executive functions. For example, in the case of prevailing verbal executive deficits, treatment may be most successful if centered on both executive and speech functions.

However, the contributions of domain-general regions to language processing remain largely obscure. If domain-general regions essentially contribute to language, one may expect overlap in the lesion-deficit associations for language and executive control, especially for language operations with higher executive demands. This hypothesis is based on the observation that under challenging conditions (e.g. noisy environments or cognitive decline), language-related activity engages both language-related as well as executive control regions. e.g.21-23

Yet, previous studies neglected a comprehensive characterization across several key cognitive domains in a single clean analysis. Indeed, most existing stroke patient cohorts are under-phenotyped and do not typically lend insight into cognitive measures that subsume various domains derived from the same cohort. Consequently, commonalities and differences between cognitive domains are unclear. Moreover, only a few studies to date take a network perspective and consider the role of both hemispheres across various cognitive domains in patients with stroke. e.g.15,16,24,25 While the few existing studies lend initial evidence for substantial lateralization differences between language and executive control, they do not address the functional interplay of domains.

Addressing these shortcomings, the present investigation was designed to tease apart the overlap and diverging consequences of lesion topologies for language and executive control in a richly phenotyped cohort of >1000 patients with stroke. To draw a complete picture of major domains for human cognition, we examined multiple verbal and non-verbal assessments with varying executive demands. To assess language function, we used a picture naming task which is a measure of basic verbal functions and language skills that probes word retrieval abilities26 and is frequently used for language assessment in clinical populations. Executive control was tested with several tasks requiring, among other functions, mental flexibility, task switching and visuospatial constructional abilities. We also tested executive speech functions and verbal memory (see next section for details about the specific tests). This comprehensive characterization, in combination with multivariate factor analysis, allowed us to extract key dimensions of human cognition and identify the impact of stroke lesions on different regions of functional capacity.

Our first hypothesis was that lesions to the left-hemispheric language network (e.g. left prefrontal, anterior and posterior temporal cortex55,27) would preferentially track language impairments. In contrast, disturbed executive control functions should index damage of the distributed (potentially right-dominant) multiple-demand network.19,28,29 Second, we expected a considerable degree of overlap of executive functions with the whole-brain lesion distribution for language assessments that require high control demands. These assessments should draw on both language-specific and domain-general control regions, with a potential strongest overlap in bilateral prefrontal control regions.19,20,30,31

To foreshadow the main results of our population-scale lesion-network investigation, our factor analysis identified four factors that define unique left- and right-hemispheric contributions to different cognitive domains. One factor delineated language and general cognitive performance and was mainly associated with damaged regions in the left hemisphere. An executive control factor summarized control and visual-constructional abilities and was strongly related to right-hemispheric brain damage. The interplay of language and executive control was reflected in two different factors, delineating executive speech functions and verbal memory. Impairments on both were mainly linked to left-hemispheric lesions. Collectively, our findings formulate new insight into causal elements of structure-function underpinnings and hemispheric specialization in an approach that cuts across cognitive domains usually studied in isolation.

Material and methods

Characteristics of the participant sample and neuropsychological tests

1080 patients from a prospective South-Korean stroke registry32 were included (614 males, 464 females, mean age 67.4 ± 0.7). Details on participant recruitment can be found in the SI Materials and Methods. MRI acquisition was performed within one week after the stroke to predict cognitive outcomes three months after the stroke. Patient inclusion was not restricted to specific lesion topographies. Most of the patients had lesions in either the left or right vascular territory of the middle cerebral artery and there was no significant difference in the number of lesioned voxels per patient between the left and right hemisphere (two-sided t-test: P = 0.81; see24 for details on the spatial coverage). The maximum of lesioned tissue was localized in subcortical zones. Details about structural brain image assessment and resulting segmented lesion patterns were described elsewhere.24,25

Patients underwent a rich battery of neuropsychological tests ∼three months after the acute onset of stroke (median time post-stroke: 98 days, Standard Deviation: 71.72 days).32 Neuropsychological and demographic data were obtained from the Korean-Vascular Cognitive Impairment Harmonization Standards-Neuropsychology Protocol.32,33 Specifically, each patient was characterized by eight key assessments of post-stroke cognitive performance: global cognitive function (Korean version of the Mini-Mental State Exam, MMSE),34 language (Korean short version of the Boston Naming Test, BNT),35 executive speech function (Korean version of phonemic and semantic fluency tests),32 executive control functions (Trail Making Test, TMT, version A and B and Digit Symbol Coding Task, DSCT),36,37 visuospatial functioning (Rey Complex Figure, RCF, including both the copying phase and the delayed recall phase)38 as well as verbal learning and memory function (Soul Verbal Learning, SVL).32

Performance in the MMSE reflects global cognition, including the orientation to time and place, as well as calculation or language performance. The BNT (Korean short version) is a standardized clinical test that measures word retrieval of patients by asking them to name 15 pictured nouns (short version). Semantic and phonemic fluency tests are prototypical measures of verbal fluency, probing executive functioning, speed and attention, and access to the mental lexicon. In semantic fluency tests, subjects are required to generate words that belong to a specific category for a limited time window. Phonemic fluency tests require the subject to generate words starting with a given letter. Since we averaged across semantic and phonemic tests, the term ‘fluency’ will henceforth be used to describe both. The TMT probes visual attention and task switching in two parts in which the subject is instructed to connect a set of dots or numbers and letters as quickly as possible while maintaining accuracy. The test provides information about visual search speed, scanning and processing speed, mental flexibility and executive functioning. The DSCT was designed to measure processing speed, working memory, visuospatial processing and attention. In this test, subjects learn a code in which each digit is represented by a symbol. Subjects have to substitute the correct symbols for a series of digits as quickly and accurately as possible. In the RCF, patients are first asked to copy a complex line drawing (copy phase) and then draw it from memory (delayed recall phase). The test measures visuospatial constructional abilities and visual memory. The SVL, in turn, examines episodic memory performance requires the auditory learning of a word list, and tests its memorization by an immediate recall task.

Finally, the Informant Questionnaire on Cognitive Decline in the Elderly (IQCODE) was used to capture the cognitive performance before the ischaemic event.39 IQCODE is an established retrospective assessment relying on health-proxy reports that provide a widely used and validated measurement instrument of cognitive decline in the ten years before stroke onset. All tests and questionnaires were conducted by a trained clinical psychometrician blinded to clinical or neuroradiological patient information.40 Additionally, we collected available sociodemographic and clinical information, such as age, sex, time since stroke onset (in days), years of education and lesion volume (grey and white matter). Each continuous (non-binary) variable was z-scored across all subjects to ensure comparability.

Statistical analysis

For details on the preprocessing, and the neuroimaging protocols, please refer to the Supplementary Information, Materials and Methods. We note that our analysis framework has two phases. In the model development phase, the lesion patterns and cognitive factors are identified. In the application phase, each individual patient has a specific combination of lesion patterns (i.e. a weighted combination of each of the derived 10 lesion patterns) and cognitive factors (identified by the multivariate factor analysis, see below).

Atlas-based extraction of target lesion load signals

The structural brain scan for each subject was characterized by a total of 435 642 grey matter voxels of 1 mm3. To provide a more generalizable and interpretable form, we parsed each patient’s lesion fingerprint by summarizing the lesion load within 54 parcels (108 for both hemispheres) based on the Harvard-Oxford cortical atlas with 47 regions and subcortical atlas with 7 regions in each hemisphere.41 We first counted the number of voxels affected per atlas-defined brain region. In doing so, we obtained 54 regional measures of lesion load per hemisphere in each participant. We then log-transformed and concatenated the ensuing lesion load measures for the left and right hemispheres.

Data-driven deconvolution of lesion patterns

After summarizing the lesion load across 54 anatomical brain regions for each subject, we sought to isolate coherent topographical patterns that may be hidden in these regional measures. To this end, we used non-negative matrix factorization (NMF)42 as a multivariate encoding strategy to identify ‘templates’ of the most common lesion locations across the entire brain. NMF decomposed the input lesion map X of n subjects and m regions into two low-rank non-negative matrices W and H, such that:

X≈WH,

where X is the [m × n] lesion load matrix, W is a [m × k] non-negative matrix for k latent patterns and H is a non-negative matrix of size [k × n].

The matrix W represents a set of non-negative basis vectors (i.e. latent factor representations), which denotes region-wise implications in each of the k patterns. The latent pattern expression matrix H indicated how relevant each emerging lesion pattern is to describe the constituent parts of an individual patient’s overall spatial lesion distribution in the brain. Specifically, our latent factorization deconvolved the actual lesion constellation into 10 unique combinations of spatially distributed region damage that will henceforth be labelled as ‘lesion atoms’ (see Fig. 1). The optimal number of 10 lesion atoms per hemisphere was validated in previous works on the identical dataset.24,25 This global decomposition of local lesion load indicators strikes a balance between capturing a substantial amount of lesion variability, on the one hand, and keeping the number of quantities low for neuroscientific domain interpretation, on the other hand.

Figure 1 Lesion atoms of stroke patterns reveal unique lesion topologies across the whole brain. Voxel-wise information on stroke-induced lesions from >1000 patients was summarized based on 108 anatomical region definitions from a reference atlas (54 per hemisphere) and derived by dimensionality-reducing pattern learning. Region-wise lesion measures were then compressed into 10 essential lesion-pattern ‘prototypes’ in each hemisphere by capitalizing on non-negative matrix factorization (NMF). A. Our derived 10 brain lesion atoms projected on the cortical surface. The resulting lesion atoms capture biologically plausible lesion pattern topographies. B. Relevance of specific brain regions within each of the 10 lesion atoms (quantified as NMF loading) shows whole-brain coverage with distributed lesions in frontal, temporal and parietal as well as subcortical regions. Lesion atoms 1 and 2 implicate the insular cortex and thalamus, respectively. Lesion atom 3 covers the prefrontal cortex (including the inferior frontal gyrus). Occipital regions are covered by lesion atom 4 (inferior occipital cortex and adjacent inferior temporal and parietal regions) and 5 (superior occipital cortex and heteromodal association regions of the inferior parietal lobe, angular gyrus, supramarginal gyrus). Temporal regions are covered by lesion atoms 6 (superior, middle and inferior temporal cortex) and 9 (fusiform gyrus, hippocampus and parahippocampal regions, as well as precuneus and cuneus in the parietal lobe). Precentral and postcentral regions are included in lesion atoms 7 and 10. Lesion atom 8 includes the basal ganglia. C. Z-scored behavioural performance. Raincloud plots show the performance of each subject for the eight cognitive assessments. Please note that the Figure visualizes the input data. We did not perform any statistical tests on these data.

Embedding lesion load across the brain into 10 lesion atoms using NMF provided at least two key advantages to alternative dimensionality reduction tools. First, in contrast to clustering approaches that consider the effect of each location only once, each brain location could belong to several latent lesion components of W to varying degrees. In this way, each location could contribute to the prediction of cognitive scores through relative contributions of multiple components, each of which reflected extracted lesion archetypes distributed across the whole brain. As a second key advantage, the non-negativity of the segmented brain lesion information and the non-negativity constraint of the NMF model allowed for intuitive neurobiologically meaningful interpretations. That is, each latent component Wk represented a unique and directly interpretable aspect of the overall topographical lesion pattern variation. The neurobiologically interpretable sum-of-parts representation enabled by NMF contrasts with latent representations learned by alternative matrix factorization algorithms. For example, in principal component analysis, individual lesions would be recovered through convoluted additions and subtractions of several components with positive and negative weights. For this reason, the overall effect of all principal components, yet not the effect of ensuring individual components, would have been as easily and intuitively interpretable to draw neuroscientific conclusions.

Latent factors driving cognitive outcomes

We explored the possibility that the observed inter-correlations between our eight target cognitive scores are influenced by one or more underlying factors that are not directly observable. This assumption is in direct analogy to the q-factor in intelligence research as well as the big-5 personality model, where factor analysis was applied to personality survey data to summarize aspects of personality into five broad dimensions.43 The factor-analysis-derived overarching domains contain and subsume the most known personality traits and are assumed to represent the basic structure behind all personality traits.44 Hence, we turned to factor analysis to uncover coherent latent patterns explaining the interrelationships among the set of our cognitive scores (Supplementary Information, Materials and Methods for details).

Predicting cognitive outcomes at 3 months post-stroke

The NMF-derived expressions of topographical lesion atoms provided the neurobiological input into our Bayesian hierarchical models.45 Carefully tailored to two different analytical strategies; we built two distinct classes of Bayesian hierarchical models. First, we designed a single multivariate, multivariable Bayesian hierarchical model to explain interindividual differences in eight cognitive outcome scores. The multivariable aspect allowed us to jointly estimate population means, variances and covariance of interrelated post-stroke outcomes. This model is labelled as MIMO (multiple-input multiple-output) throughout the manuscript. In this step, we derived a set of four multivariate Bayesian hierarchical models dedicated to the four factors representing investigated cognitive domains. In other words, each model distilled knowledge from multiple inputs to provide predictions for a single output (i.e. factor). This setup is referred to as MISO (multiple-input single-output) throughout the manuscript. Finally, all models took into account several covariates, including age, age,2 sex, time since stroke onset, education years, pre-morbid cognitive performance and total lesion volume following previous research.24,25 Details on the Bayesian model specification and posterior Bayesian estimates for all parameters and models are available in the SI Materials and Methods section.

Our Bayesian hierarchical approach facilitated the careful dissection of predictive relevance allocated to different levels of the model. For all outcome models, we first evaluated lateralization effects inferred from the left and right hemisphere posterior dispersion distributions. Subsequently, we considered the lateralization effects of specific lesion atoms and, lastly, reverted back the predictive relevance of lesion atoms46 to the level of the anatomical brain regions for each modelled outcome. Finally, we included age of education and pre-stroke cognitive decline (estimated by the IQCODE39) as covariates to our Bayesian model to quantify their impact on Lesion atom-outcome associations for each factor in both hemispheres.

Quantifying the similarity with language and multiple-demand network

We explored the similarity of our factor-specific brain maps of beta coefficients with maps of the language network and domain-general multiple-demand network.47,48 Both maps are freely available online (https://evlab.mit.edu/funcloc/). As a first step, for each factor-specific Bayesian model, we inverse-projected beta coefficients from the latent space of non-negative matrix factorization to the space of the Harvard-Oxford atlas. Thus, we obtained a coefficient for each region of this atlas. Then, we projected the absolute values of region-specific coefficients on the cortex. Finally, to quantify the similarity, we computed Pearson’s r using neuromaps toolbox49 between the brain map of coefficients pertinent to a given factor and brain maps of language as well as multiple-demand networks.

Results

We first derived a whole-brain representation of distributed lesion patterns from expert-segmented post-stroke lesion maps. Next, we devised and deployed a principled Bayesian framework. A fully probabilistic framework is well suited to trace how brain tissue impairments impact interindividual differences in key cognitive domains after acute ischaemic stroke24,25 and quantify the extent of overlap of these brain-behaviour correspondences. The combination of a comprehensive cognitive test battery with factor analysis in a Bayesian modelling framework should capture the impact of stroke-induced dysfunction in distributed brain regions for key faculties of human cognition.

Stroke lesion atoms isolate a typology of unique lesion constellations

We first summarized the high-dimensional lesion information at the voxel level in 108 cortical and subcortical brain regions (54 per hemisphere, based on the Harvard-Oxford cortical atlas). To uncover coherent hidden patterns of lesion topography from the brain scan data, we used a previously established set of 10 distinct topographical lesion configurations, so-called ‘lesion atoms’ (lesion patterns) based on non-negative matrix factorization (NMF).24,25,42 Note that the definition of lesion patterns was symmetrically identical in both hemispheres. The spatial definition of the atomic lesion patterns corresponded to biologically plausible components of stroke-induced brain lesions in both hemispheres and their co-occurrences across patients (Fig. 1A). Accordingly, lesion atoms were reminiscent of territories of arterial blood supply via the anterior (anterior part of lesion atom 3, midline structures of lesion atom 7), middle (lesion atoms 1, 4, 5, 6, 7, 8, 10, posterior part of lesion atom 3) and posterior (lesion atoms 2 and 9) cerebral artery territories. Important to our present analyses, our collection of lesion atoms showed whole-brain coverage with distributed lesions in frontal, temporal and parietal as well as subcortical regions (Fig. 1B for details).

Having parsed the lesion-segmented brain scans by means of the recurring stroke topographies, we investigated the overlap versus dissociation of lesion patterns on cognitive outcomes after stroke (Fig. 1C for the distribution of cognitive scores).

Bayesian modelling predicts atom-impairment-links in a hemisphere-aware fashion

Next, we performed a Bayesian analysis to test how the patients’ configurations of lesion atoms explain interindividual differences in post-stroke outcomes. Since these post-stroke outcomes are expected to be interrelated, we derived a single multiple-input multiple-output (MIMO) model. The MIMO model specification allowed us to jointly estimate population means, variances and covariance of our collection of eight outcome measures. This analysis can infer accurate uncertainties for the contribution of each hemisphere and their dependent lesion atom to the successful endpoint prediction at the single-patient level. Our brain measurements were modelled by taking into account potential sources of variation outside of primary scientific interest: total lesion volume, age, age,2 sex, education in years, and premorbid cognitive performance.24,25 As a litmus test of model quality, the posterior predictive checks confirmed the adequate fit of our estimated MIMO prediction model. The variation was described by how individual lesion atom expressions explained differences in cognitive performance outcomes in patients,50 albeit with different achieved performances across the target cognitive scores (coefficient of determination R2 ranged from 0.11 to 0.49 for phonemic fluency and MMSE, respectively, Fig. 2A). This empirical verification of the MIMO model’s ability for brain-based predictions of post-stroke capacities is a well-recognized proxy of external validation given the actual data at hand.51 The estimated standard deviation posteriors captured the overall left and right hemispheric relevance with good certainty (i.e. the width of the posterior variance parameter distributions) and dissociated hemisphere-specific predictive contributions tiled across all candidate lesion atoms (Fig. 2B).

Figure 2 Bayesian hemisphere-aware analysis can robustly predict individual clinical outcomes. A. Model performance of the inferred MIMO (multiple-input multiple-output) Bayesian analytical solution using brain lesion atoms in predicting cognitive impairments. Posterior predictive checks were performed for the Bayesian models that were estimated to predict interindividual differences in eight cognitive measures (on z-scale). Model-based simulations (y-axis) were compared to the observed data (x-axis) to compute the overall explained variance (coefficient of determination, R2). B. Posterior probabilities for the hemisphere-specific lateralization effects for each clinical endpoint. Lesions to the left hemisphere were more relevant for single-patient predictions of cognitive outcomes in the Mini Mental State Exam, BNT, Seoul Verbal Learning, verbal fluency tests, and DSCT. In contrast, damage in the right hemisphere was more explanatory for performance in the RCF and TMT. The posterior model parameters correspond to the upper hemisphere level of our Bayesian multilevel modeling strategy that uncovered the hemisphere-specific model certainty for each cognitive performance dimension. Surface projections: blue colours: lesions are associated with relatively stronger impairments, red colours: lesions are associated with relatively weaker impairments. Std., standard deviation. C. The covariance matrix quantifies the co-occurrences between cognitive outcome measures across patients. Covariance is high (>0.5) for Mini Mental State Exam, BNT and Seoul Verbal Learning, as well as for verbal fluency measures and Digit Symbol Coding Test. In comparison to empirical covariance, the model covariance matrix quantifies the intrinsic relationship between cognitive outcomes given brain lesion conditions and other covariates. D. Similarity of brain-behaviour representations: Hierarchical clustering of model-specific coefficients obtained per clinical endpoint shows distances (similarities) among cognitive measures based on Ward’s method. A first cluster summarizes both fluency assessments and Digit Symbol Coding Test. Based on distances, RCF and TMT can also be grouped. Finally, Seoul Verbal Learning, Mini Mental State Exam and BNT show low distances, suggesting high similarity. Overall, our hemisphere-aware Bayesian model revealed differences in cognitive outcome predictions that follow the distinction between measures of language and executive function.

As a central feature of our findings, lesions in both hemispheres contributed in specific amounts to prediction success across the examined cognitive endpoints. Lesions in the left hemisphere were overall more relevant for single-patient predictions of most cognitive outcomes. More specifically, the estimated model parameters for hemisphere-specific dispersion showed relative left-biased lateralization effects for MMSE, BNT, verbal fluency tests, Seoul Verbal Learning and DSCT; but relative right-biased lateralization for RCF and TMT. These findings were in line with the empirical and model-derived covariance matrices between cognitive outcome measures which also revealed strong similarity between MMSE, BNT and SVL, as well as between both fluency measures and DSCT (Fig. 2C). Finally, we used hierarchical clustering of the derived beta estimates (corresponding to lesion atom involvements from the estimated MISO prediction model) to explore similarities between cognitive outcome tests (Fig. 2D). This post-hoc analysis showed similarities for both fluency assessments and DSCT. Rey Complex Figure and TMT were also close to each other in this clustering space. The last cluster based on the Ward similarity measure was formed by SVL, MMSE and BNT. Collectively, our Bayesian approach revealed hemisphere-specific differences between cognitive measures based on predictions from lesion topographies. Our observations suggested that some of the cognitive measures can be considered as representative cognitive key dimensions.

Four impairment dimensions reveal hemisphere-specific effects on language and executive functions

Our exploratory factor analysis extracted principled representations of cognitive domains that summarized different outcome measures. Based on the screen test criterion, we favoured a four-factor solution that characterized the behavioural variation across eight clinical endpoint assessments (Fig. 3A). We labelled the factors based on the task functions they represent. Factor 1 was labelled as executive speech functions because it explained the combined variation of phonemic and semantic fluency and Digit Symbol Coding Test in our patients. Factor 2 tracked the overlap of naming and general cognitive performance by explaining a notable variation of MMSE, BNT and with smaller contributions estimated for RCF and Seoul Verbal Learning. The overlap may reflect both language abilities as well as perceptual and general task processing components. We labelled this factor as language and general cognitive performance.

Figure 3 Four factors trace out the similarity and dissimilarity between language-related and executive function-related clinical endpoints. Low-dimensional representation of clinical endpoints dissociates between language-related and executive control-related factors. A. Scree plot illustrating the cumulative amount of explained variance (cum. expl. variance) in outcome measures and associated eigenvalues. According to the screen test criterion, a four-factor solution is an effective low-dimensional data representation that explains a substantial portion of the variation. B. Factor loadings for the included outcome measures. C. Factor biplot depicting participant scores and loading vectors for the first three factors. Loading vectors represent how strongly each characteristic influences the resulting factor. The angle between a pair of vectors corresponds to the correlation between the given characteristics. Participant factor scores are displayed as points in the plane formed by three principal components. Dot colours code for factor scores (bright colour: high score). D. Performance of the inferred MISO (multiple-input single-output) Bayesian analytical solutions in predicting cognitive impairments for each factor. Posterior predictive checks are shown for the Bayesian models that were estimated to predict interindividual differences in (z-scored) factors. Model-based simulations (y-axis) were compared to the observed data (x-axis) to compute the overall explained variance (coefficient of determination, R2). Right panels: Relevance of the posterior parameter distributions (std = standard deviation) of the left and right hemispheres obtained through separate Bayesian hierarchical models dedicated to predicting each factor. Lesions to the left hemisphere were more relevant for single-patient predictions of cognitive outcomes for factors 1, 2 and 4. In contrast, lesions in the right hemisphere were more relevant for outcome predictions for factor 3. The posterior model parameters correspond to the upper hemisphere level of our Bayesian multilevel modeling strategy that uncovered the hemisphere-specific model certainty for each factor. The four-factor solution, in combination with four Bayesian models, successfully untangled the unique contributions of cognitive outcomes along the language versus executive control axis.

Factor 3 captured executive control, as shown by the highest loadings of TMT and lower contributions of RCF. Finally, factor 4 represented verbal memory with high loadings of Seoul Verbal Learning. The four-factor solution (Fig. 3B and C) supported and complemented the results from the hierarchical clustering approach (Fig. 2D). Specifically, the four-factor solution revealed that some impairment dimensions of language and executive functions are clearly distinguishable, while there is a shared core in the dimension of executive speech functions and verbal memory.

To explore the brain lesion implications of the derived four intrinsic cognitive factors, we estimated a set of four Bayesian hierarchical models. This multiple-input single-output (MISO) model class was carried out so that each model was dedicated to a single factor. The probabilistic prediction of patient-specific factor expressions across the four MISO models achieved 9–45% explained variance (R2, coefficient of determination; factor 1: 9%, factor 2: 45%, factor 3: 27%, factor 4: 16%). In line with the modelling results for the single outcome measures, factors 1, 2 and 4 were left-lateralized, whereas factor 3 showed right-hemispheric lateralization (Fig. 3D).

Cognitive factors associated with distinct lesion patterns in specific brain regions

Next, we explored how the derived four MISO models tied into four cognitive factors (i.e. cognitive symptom sets that vary across patients) are associated with specific atlas regions (Fig. 4) and particular lesion atoms (Fig. 5, see next section). Individual subject scores across lesion atoms and cognitive factors are provided in SI Fig. 1. For cognitive factor 1 (indexing primarily verbal fluency and executive control), we found that lesions to the left putamen, precentral gyrus, caudate nucleus, posterior superior and middle temporal gyrus (STG and MTG), and supramarginal gyrus were driving deficits in executive speech functions. In contrast, lesion-deficit associations for cognitive factor 2 (particularly reflecting word retrieval and overall cognitive abilities) were located more posterior, with the strongest impairments in language and general cognitive performance caused by lesions to the left lateral occipital regions, lingual gyrus, hippocampus and occipital-fusiform gyrus. Complementing the pattern for factor 2, cognitive impairments in factor 3 (indexing primarily executive control) showed prominent associations of executive deficits with right occipital regions, as well as the postcentral gyrus, superior parietal lobe and occipital-fusiform gyrus. For factor 4 (indexing primarily verbal learning and memory), we found relevant associations between cognitive deficits in verbal memory and lesions located in the left post- and precentral somatomotor regions, while damage to the left and right occipital regions was associated with preserved functions. Overall, the derived cognitive impairment dimensions revealed distinct ties to brain lesion patterns in both hemispheres, following a left-versus-right and anterior-posterior axis: impairments in executive speech functions were associated with left-hemispheric damage to subcortical, precentral and more anterior temporal regions. Deficits in general cognitive and language abilities were related to more posterior damage, including in the left occipital and temporal regions. Impaired executive functions were associated with lesions in the right occipital and temporal regions. Finally, deficits in verbal memory were related to lesions in the left postcentral regions. The overall distribution of lesion patterns across factors showed that language and executive functions mainly dissociate in posterior regions (including occipital, posterior temporal and subcortical regions), with the language being located relatively more anterior. In contrast, the overlap of both functions was found in more anterior regions, including pre- and postcentral somatomotor gyri of the cerebral cortex.

Figure 4 Four factor-specific models trace out distinct stroke-induced lesion patterns for single regions across different cognitive dimensions. Left side: Lesion-deficit prediction for each of the four factors. Coloured brains reflect associations of brain regions with lost (blue colour) or preserved (red colour) cognitive functions. Note that factors 1, 2 and 4 show stronger left-hemispheric lateralization, while factor 3 is right lateralized. Right side: 54 parcels per hemisphere were included based on the Harvard-Oxford cortical and subcortical atlas. The predictions of four factors with differential contributions to language versus executive functions are subserved by distinct brain patterns. β, beta value.

Figure 5 Distinct associations of lesion atoms and cognitive factors disentangle language and executive deficits. Results for the four-factor solution. A. Lesion-deficit associations for the different factors and lesion atoms. Factors 1, 2 and 4 show the most prominent brain lesion effects for atom 6 in the left hemisphere (including temporal regions), while factor 3 shows strong implications for atom 4 (summarizing occipital regions) in the right hemisphere. B. Pearson’s correlations (r) between model-specific beta values dedicated to each factor. Left-dominant factors show stronger correlations with each other. Asterisk denotes False Discovery Rate (FDR)-corrected P-values obtained from a spin-permutation test across the whole brain. C. Distinct relationships between strongest lesion atoms and cognitive impairments. Left: lesion atom 4 in the left hemisphere, with the overall strongest load of the left lateral occipital cortex, is the strongest contributor to factor 3. Middle: Overlap and distinct contributions of the four factors, where each factor is represented by the highest absolute atom loading. Right: lesion atom 6 is characterized by strong left-hemispheric contributions from factors 1, 2, 4, with the strongest load of the middle temporal gyrus. L, R, left, right; NMF, non-negative matrix factorization. D. Correlations of factors with maps of the language network and multiple-demand network47,48 show a relatively stronger overlap of executive speech functions and language with the language network. Executive functions and verbal memory show relatively stronger associations with the multiple-demand network. Please note that we illustrate an exploratory correlation analysis but none of the values reached significance after correction for multiple comparisons (using False Discovery Rate (FDR)-corrected thresholds of P < 0.05). E. Hemispheric relevance depends on key covariates. The plot depicts the interrelation between the relevance of lesion load in each hemisphere (y-axis, left hemisphere: light colours, right hemisphere: dark colours) and marginal posterior parameters of two key covariates (x-axis) in predicting cognitive performance. Years of education had positive effects on all four factors. The strongest effect was on factor 2. Conversely, an increase on the Informant Questionnaire on Cognitive Decline in the Elderly (IQCODE) scale,39 i.e. a higher pre-stroke cognitive decline, predicted a decrease in scores of factors 1, 2 and 4. Lesion atom-outcome associations uncovered factor-specific language and executive control deficits: a prominent contribution of the lesion atom covering left-hemispheric superior and middle temporal regions to language, executive speech functions and verbal memory and strong implications of right occipital brain regions for executive functions. IQ, Informant Questionnaire scale. β, beta value.

Lesion atom-outcome associations uncover factor-specific language and executive control deficits

Finally, we complemented region-based lesion-deficit associations with pattern-level results based on the identified lesion atoms. Factors 1, 2 and 4 showed the most prominent brain lesion effects for atom 6 (including inferior, middle and superior temporal regions) in the left hemisphere (Fig. 5A). These factors also shared lesion effects in lesion atom 8 in the left hemisphere (covering the basal ganglia). In contrast, factor 3 featured the strongest brain effects for atom 4 (including occipital and adjacent inferior temporal and parietal regions) in the right hemisphere (Fig. 5C). Moreover, factor 3 showed an overlap of lesion effects with factor 4 for lesion atoms 4, 6 and 7 in the right hemisphere (covering occipital, temporal, inferior parietal, precentral and middle to superior frontal areas).

Additionally, factor 1 showed a robust implication of tissue lesions falling into the territories of lesion atoms 5 and 7 in the left hemisphere (covering the superior occipital cortex and adjacent inferior parietal regions, precentral cortex and adjacent middle as well as inferior frontal regions). Factor 2 showed additional effects for atoms 2, 3 and 4 in the left hemisphere (covering thalamus, (pre-)frontal cortex and occipital cortex) and atom 1 in the right hemisphere (insular cortex). Factor 3 had additional implications of atom 1 in the left hemisphere (insular cortex) and atoms 9 and 10 in the right hemisphere (covering temporal regions, (para-)hippocampus and adjacent fusiform gyrus, thalamus, occipital pole as well as pre- and postcentral cortex, inferior and middle frontal regions and the superior parietal cortex). Finally, factor 4 was additionally associated with brain lesion effects captured by all atoms except atom 9 in the left hemisphere and atom 2 in the right hemisphere (thalamus). Overall, these results showed that stroke-induced language deficits could be mainly characterized by left-hemispheric lesions to prefrontal and temporal regions, while executive deficits are related to right-hemispheric damage in occipital and postcentral regions. Impairments in the overlap of both domains, summarized as executive speech functions and verbal memory, were associated with more distributed lesions in the left precentral and prefrontal regions as well as posterior temporal and parietal regions.

Notably, we found correlations between factors with left-hemispheric dominant lesion-outcome associations (factors 1, 2, 4, Fig. 5B). Factors 3 and 4 were also correlated.

Correlating our factors with maps of the language network and domain-general multiple-demand network43,44 revealed a relatively stronger overlap with the language network for factors 1 and 2 (executive speech functions as well as language and general cognitive performance). In contrast, factors 3 (executive functions) and 4 (verbal memory) showed more overlap with the multiple-demand network (Fig. 5D). We also explored interrelations between important covariates (i.e. age and years of education, as well as pre-stroke cognitive decline and number of lesions) in the prediction of cognitive outcomes for each factor (Fig. 5E). While we observed a positive impact of education on cognition-related factor scores, the inverse was true for cognitive decline.

Collectively, our results showed unique hemisphere-specific contributions of language and executive functions. These specific contributions dissociated at the level of brain regions (including prominent contributions of temporal, (pre-)frontal and occipital regions) and networks. Aggregated cognitive factors derived from the factor analysis showed stronger associations of executive speech functions (factor 1) as well as language and general cognitive abilities (factor 2) with the previously described language network.42 In contrast, executive control (factor 3) and verbal memory (factor 4) were relatively more strongly associated with the multiple-demand network.49 While factors 1, 2 and 4 showed prominent associations with left-hemispheric lesion atoms, factor 4 showed additional brain effects in right-hemispheric lesion atoms, and factor 3 revealed a prominent association with right-hemispheric lesion atoms. The observed constellations of specific and overlapping associations helped to uncover unique and general contributions and argue for a complex interplay of specialized and domain-general networks and brain regions across different key domains of human cognition.

Discussion

Here, we shed light on the causal relevance of distributed brain regions for language versus executive control functions. We combined Bayesian hierarchical modelling with a rich battery of neuropsychological outcome measures in a large cohort of >1000 patients with stroke. Seizing this opportunity, we were able to robustly predict individual clinical outcomes based on specific tissue lesion topologies. We identified four distinct cognitive factors that characterized the similarity and dissimilarity between language-related and executive control-related clinical outcomes. Our results elucidate specific lesion topologies associated with deficits in each of the two cognitive domains and detail their functional interplay.

As a core finding, our analytical framework could disentangle anterior-posterior and left-right principles for impairment of language and executive functions due to distributed stroke insults. Language and general cognitive performance impairments were primarily linked to lesions of prefrontal and temporal regions in the left hemisphere, while executive control was preferentially affected by damage to occipital and postcentral regions in the right hemisphere. As a second main finding, the interplay of both cognitive dimensions was reflected in two distinct factors, tracking executive speech functions and verbal memory capacities. Contrasting our expectations, impairments on these two cognitive factors were mainly linked to left-hemispheric lesion atoms, with only minor contributions of right-hemispheric regions. These identified left-hemispheric regions covered pre- and postcentral regions associated with cognitive control as well as motor cortices. Considering specific lesion topologies in different areas across the whole brain provided a fine-grained picture of impairment dimensions associated with distinct lesion atom patterns and specific brain regions for each identified cognitive factor.

Overall, our findings illustrate that stroke-induced lesions in the left hemisphere were relatively more relevant for single-patient predictions of most examined cognitive outcomes. Accordingly, factors 1, 2 and 4 showed left-hemispheric dominance of lesion-outcome associations with the most prominent brain lesion effects for atom 6 in the left hemisphere. This lesion atom tracked grey matter lesions across the superior, middle, and inferior temporal brain regions, which are associated with language operations and global cognitive functions in studies of healthy and lesioned brains.5,52-55 The observed dominance of the left hemisphere in our findings may reflect our selection of cognitive outcome measures with excellent coverage of verbal memory as well as speech and language-related faculties (i.e. MMSE, BNT, verbal fluency, Seoul Verbal Learning and DSCT). Accordingly, the verbal component of the tasks contributing to factors 1 and 4 is a likely explanation for the unexpected left-lateralization of both factors. Based on previous results, we initially expected that verbal tasks with a strong executive load (particularly phonemic fluency) should be more bilaterally distributed and engage the right inferior frontal cortex and insula (e.g.;56 see57 for a meta-analysis). Likewise, the Seoul Verbal Learning test has been associated with bilateral frontal activity, which was linked to the memory component of the task.58 Accordingly, verbal learning impairments have been related to stroke-induced lesions of the right hemisphere before.59 Indeed, for factor 4, we found associations with right-hemispheric lesion atoms, including prefrontal, temporal and parietal regions. However, verbal fluency selectively revealed associations with left-hemispheric lesion atoms.

Nevertheless, regarding the lesion patterns for the executive control-related factor 3, we were also able to isolate specific contributions of right-hemispheric regions. We can only speculate how the inclusion of more tests for spatial attention and executive functions might have influenced our results. In case of overlap with our executive control-related factor, we would expect strong loadings of such tests on factor 3. Alternatively, one might expect a more fine-grained picture, leading to an additional factor with right-hemispheric coverage that might, for example, reflect visuospatial attention.

Moving beyond incumbent narratives from previous stroke research, our results further show that the shared substrates of language and executive control functions, that is, executive speech functions and verbal memory, mainly suffer from distributed left-hemispheric damage. Associated lesion atoms included both language regions as well as pre- and postcentral motor and control regions. In this context, we wish to highlight that relative to classical voxel-wise lesion-symptom mapping approaches, our patient-specific predictions are better suited to forecasting clinical endpoints.24

Accordingly, hierarchical clustering and exploratory factor analyses quantitatively dissected the underlying cognitive dimensions in both hemispheres. The resulting factors showed links to specific lesion-deficit associations, establishing an anterior-to-posterior axis for executive speech versus language operations. The first factor summarized executive speech functions. Impairment of this cognitive factor was linked to lesions in distributed left-hemispheric brain regions. These included posterior and middle temporal as well as inferior parietal regions, which are frequently associated with phonological and verbal working memory processes (superior temporal gyrus and supramarginal gyrus60-63) or semantic control processes (middle temporal gyrus64,65). These neural processes are known to be relevant for verbal fluency. Other regions comprised left precentral and subcortical regions in the basal ganglia (putamen and caudate nucleus), linked to sequencing operations during speech processing and often reported in fluency tasks in the healthy brain.66,67 While the contribution of subcortical regions beyond the thalamus to clinical language impairments is still debated,68,69 our in-depth quantifications add evidence for a key contribution of the basal ganglia in the left hemisphere to fluent speech operations. Notably, our findings in single brain regions were supported by complementary information from lesion atom-deficit associations, providing evidence for the relevance of the left precentral cortex and adjacent middle and inferior frontal regions, which are key regions for executive semantic processes31,64,70 and implicated in executive speech control.71 These results open a new window into how executive control is necessary for fluent language performance. Strikingly, the overlap in the lesion patterns shows that rather than recruiting additional control regions in the opposite hemisphere, fluent speech processing mainly draws on distinct left hemispheric regions, with a strong emphasis on the prefrontal and temporal cortex. This insight paves the way for future individualized lesion- and deficit-oriented therapeutic approaches that should include training of domain-general executive control elements for patients with impaired speech fluency.

The second factor particularly prioritized interindividual differences in language and general cognitive functions. Impairments were linked to lesions located in more posterior temporal regions of the left hemisphere which are associated with general cognitive abilities, as well as language and semantic memory.71-73 These regions included lateral occipital regions, lingual gyrus, hippocampus and occipital-fusiform gyrus. Results from lesion atoms supported these findings and emphasized the role of left temporal and frontal regions for language and general cognitive abilities. (see also24) Collectively, the results for factor 2 highlight a strong contribution of lateral temporo-occipital regions to general cognitive and language functions associated with word retrieval. The additional association of factor 2 with lesion atom 2 (covering the right insula) likely reflects more general control-related processes of the Mini Mental State Exam and picture naming task. Indeed, the right insula is often engaged in production tasks, especially with an increased task load (e.g.74).

The third factor carefully separated the stroke symptom component indexing central executive control, summarizing two measures (TMT and RCF) that are commonly considered to reflect executive control and visual-constructional abilities.37,75-78 A shared underlying cognitive dimension for both measures is supported by the observed similarities in our hierarchical clustering analysis. Cognitive impairments in factor 3 showed prominent associations of deficits in executive control functions with right occipital regions, as well as the postcentral gyrus, superior parietal lobe and occipital-fusiform gyrus. Accordingly, factor 3 showed strong implications for lesion atom 4, encompassing occipital and adjacent parietal regions in the right hemisphere. This lesion-deficit pattern likely reflects the strong reliance of the underlying tasks (TMT and RCF) on visual scanning and visual-constructional abilities.76,79 We note that the observed differences in the lesion-deficits patterns between factors 1 and 3 offer a fine-grained parcellation of the neural proxies of executive control. This implies that visuo-constructional abilities draw more strongly on posterior right hemispheric regions while verbal executive control is primarily related to left prefrontal and temporal regions. We wish to emphasize that our results do not stand against the contribution of more anterior (right-hemispheric frontal) regions to executive functions which may have been better covered by other cognitive control tasks.

Collectively, our results move beyond a classic ‘left hemisphere = language, right hemisphere = executive functions’ distinction, arguing for a key contribution of left hemispheric regions to verbal executive functions. As a limitation of the present study, handedness as a surrogate for the dominant hemisphere was not recorded. We note that the two tests with the highest factor loadings for factor 3 (TMT and RCF) require hand-motor function. However, we could not assess the number of patients who performed these tests with their non-dominant hand due to stroke-induced motor impairment, although this may have influenced task performance and lateralization.

The fourth factor described verbal learning and memory performance in our patients, emphasizing the Seoul Verbal Learning task. Cognitive impairments were tied to damage in the left post- and precentral motor regions. Complementary evidence from lesion atom-deficit associations showed contributions of distributed regions in both hemispheres, likely reflecting the domain-general aspect of the SVL, drawing on verbal learning and memory functions.32,80 We note a close relationship of factor 4 with the neuropsychological assessments characterizing factor 2, as reflected in relatively high loadings of the SVL on factor 2, the strong positive correlation between factor 2 and 4, and the similarities between the SVL and the MMSE and BNT in the hierarchical clustering analysis. These results show that verbal memory is a unique cognitive faculty that quantifies the memory-related intersection of language and executive control, which can be clearly distinguished from executive speech functions as reflected in factor 1.

Finally, correlations with well-described brain networks for language and domain-general cognitive control47,48 provided complementary insight into lesion-deficit topologies. The two speech and language-related factors 1 and 2 showed positive correlations with the language network and negative correlations with the domain-general MDN. While this result was expected for the language-related factor 2, the specific overlap for factor 1 with language rather than MDN regions is surprising. This finding implies that executive speech functions may draw more strongly on language-related control processes rather than domain-general control functions associated with the MDN. This interpretation is supported by previous work arguing for a strong contribution of left prefrontal and posterior middle temporal regions to semantic control processes,31,64,80 which are also necessary for accurate semantic fluency performance.81,82 However, we note that the employed language localization atlas48 was largely based on sentence reading and not optimized for covering overt production and speech motor processes. These processes are associated with premotor and medial frontal areas, including the pre-supplementary motor area.55 The limited coverage of medial frontal areas the in the language localization atlas may explain the relatively low correlation for the two factors that engage overt speech and language functions (factors 1 and 2).

The remaining factors 3 and 4 showed weak positive correlations with the MDN. The association of the right-dominant executive control factor 3 with this domain-general network was expected.8,16,24 In contrast, a correlation of the left-dominant factor 4 with the MDN may be surprising at first glance but likely reflects the learning and memory-related component of the underlying cognitive task. This association is likely explained by the contribution of the bilateral parietal cortex. The overlap of factor 4 with general control functions is further supported by the relatively strong correlation of factors 3 and 4.32,83

It should be borne in mind that lesions are usually not equally distributed across cortical areas. Most stroke-induced brain lesions are caused by middle cerebral artery infarctions and often affect frontal areas. However, midline structures such as the pre-supplementary motor area are usually spared by MCA infarctions because they are supplied by the anterior cerebral artery. This may have contributed to the relatively low prediction accuracy of our model for the fluency measures (and thus, factor 1) since fluency is known to be strongly associated with medial frontal areas, especially the pre-supplementary area (e.g.81).

In summary, our results provide insight into common and distinct lesion-deficit patterns for language and executive control. The identified four factors of cognitive facets may inform future personalized cognitive therapy approaches optimized for individual lesion-deficit profile patterns. Individual deficit profile patterns may identify patients at risk for specific cognitive deficits early after stroke to enable targeted cognitive testing and therapy (see84). While previous studies in patients with post-stroke aphasia suggested that recruitment of domain-general regions for cognitive control is linked to favourable language recovery,17-20 none of these studies provided insight into the overlap of language and executive control functions. Our results fill this gap by showing that language-related control is associated with left-hemispheric regions while executive control is associated with the right hemisphere. Strikingly, the overlap of both functions is also mainly located in left-hemispheric regions. Yet, to date, individual lesion patterns are rarely considered for optimizing cognitive therapy. Based on our findings, we may speculate that lesion topologies associated with factor 2 may benefit most from language-specific treatment. In contrast, lesion topologies associated with factors 1 and 4 would additionally require the inclusion of more executive therapy elements.

Supplementary Material

fcae129_Supplementary_Data

Acknowledgements

The authors are grateful to Angelina K. Kancheva and Goezdem Arikan for their help with performing the manual infarct segmentations.

Supplementary material

Supplementary material is available at Brain Communications online.

Funding

DB was supported by the Brain Canada Foundation, through the Canada Brain Research Fund, with the financial support of Health Canada, National Institutes of Health (NIH R01 AG068563A, NIH R01 R01DA053301-01A1), the Canadian Institute of Health Research (CIHR 438531, CIHR 470425), the Healthy Brains Healthy Lives initiative (Canada First Research Excellence fund), Google (Research Award, Teaching Award), and by the CIFAR Artificial Intelligence Chairs program (Canada Institute for Advanced Research). GH was supported by Lise Meitner Excellence funding from the Max Planck Society, the European Research Council (ERC-2021-COG 101043747) and the German Research Foundation (HA 6314/3-1, HA 6314/4-2, HA 6314/9-1). Publication of this article was funded by the Open Access Publishing Fund of Leipzig University supported by the German Research Foundation within the program Open Access Publication Funding.

Competing interests

The authors report no competing interests.

Data availability

Anonymized data that support the findings of this study are available from the corresponding author DB upon reasonable request. Analyses were conducted in a Python 3.7 environment and predominantly relied on the packages nilearn and pymc3. Full code is accessible to and open for reuse by the reader here: https://github.com/jakubkopal/bayesian_stroke.
==== Refs
References

1 GBD 2019 Stroke Collaborators . Global, regional, and national burden of stroke and its risk factors, 1990–2019: A systematic analysis for the global burden of disease study 2019. Lancet Neurol. 2021;20 (10 ):795–820.34487721
2 Tarantino  V, Burgio  F, Toffano  R, et al  Efficacy of a training on executive functions in potentiating rehabilitation effects in stroke patients. Brain Sci. 2021;11 (8 ):1002.34439621
3 Glosser  G, Goodglass  H. Disorders in executive control functions among aphasic and other brain-damaged patients. J Clin Exp Neuropsychol. 1990;12 (4 ):485–501.1698809
4 Chung  C, Pollock  A, Campbell  T, Durward  B, Hagen  S. Cognitive rehabilitation for executive dysfunction in adults with stroke or other adult nonprogressive acquired brain damage. Stroke. 2013;44 (7 ):e77–e78.23798559
5 Friederici  AD . The brain basis of language processing: From structure to function. Physiol Rev. 2011;91 (4 ):1357–1392.22013214
6 Miyake  A, Friedman  NP, Emerson  MJ, Witzki  AH, Howerter  A, Wager  TD. The unity and diversity of executive functions and their contributions to complex “frontal lobe” tasks: A latent variable analysis. Cogn Psychol.  2000;41 (1 ):49–100.10945922
7 Stolwyk  RJ, Mihaljcic  T, Wong  DK, Chapman  JE, Rogers  JM. Poststroke cognitive impairment negatively impacts activity and participation outcomes: A systematic review and meta-analysis. Stroke. 2021;52 (2 ):748–760.33493048
8 Schumacher  R, Halai  AD, Ralph  L, A  M. Assessing and mapping language, attention and executive multidimensional deficits in stroke aphasia. Brain. 2019;142 (10 ):3202–3216.31504247
9 Bzdok  D, Ioannidis  JPA. Exploration, inference, and prediction in neuroscience and biomedicine. Trends Neurosci. 2019;42 (4 ):251–262.30808574
10 Bzdok  D, Meyer-Lindenberg  A. Machine learning for precision psychiatry: Opportunities and challenges. Biol Psychiatry Cogn Neurosci Neuroimaging. 2018;3 (3 ):223–230.29486863
11 Hamzei  F, Erath  G, Kücking  U, Weiller  C, Rijntjes  M. Anatomy of brain lesions after stroke predicts effectiveness of mirror therapy. Eur J Neurosci. 2020;52 (6 ):3628–3641.32031282
12 Crinion  JT, Leff  AP. Using functional imaging to understand therapeutic effects in poststroke aphasia. Curr Opin Neurol. 2015;28 (4 ):330–337.26110806
13 Hartwigsen  G, Saur  D. Neuroimaging of stroke recovery from aphasia—Insights into plasticity of the human language network. Neuroimage. 2019;190 :14–31.29175498
14 Turkeltaub  PE, Messing  S, Norise  C, Hamilton  RH. Are networks for residual language function and recovery consistent across aphasic patients?  Neurology. 2011;76 (20 ):1726–1734.21576689
15 Siegel  JS, Ramsey  LE, Snyder  AZ, et al  Disruptions of network connectivity predict impairment in multiple behavioral domains after stroke. Proc Natl Acad Sci U S A. 2016;113 (30 ):E4367–E4376.27402738
16 Siegel  JS, Shulman  GL, Corbetta  M. Mapping correlated neurological deficits after stroke to distributed brain networks. Brain Struct Funct. 2022;227 (9 ):3173–3187.35881254
17 Geranmayeh  F, Brownsett  SL, Wise  RJ. Task-induced brain activity in aphasic stroke patients: What is driving recovery?  Brain. 2014;137 (Pt 10 ):2632–2648.24974382
18 Geranmayeh  F, Chau  TW, Wise  RJS, Leech  R, Hampshire  A. Domain-general subregions of the medial prefrontal cortex contribute to recovery of language after stroke. Brain. 2017;140 (7 ):1947–1958.29177494
19 Brownsett  SL, Warren  JE, Geranmayeh  F, Woodhead  Z, Leech  R, Wise  RJ. Cognitive control and its impact on recovery from aphasic stroke. Brain. 2014;137 (Pt 1 ):242–254.24163248
20 Stockert  A, Wawrzyniak  M, Klingbeil  J, et al  Dynamics of language reorganization after left temporo-parietal and frontal stroke. Brain. 2020;143 (3 ):844–861.32068789
21 Vaden  KI  Jr, Kuchinsky  SE, Cute  SL, Ahlstrom  JB, Dubno  JR, Eckert  MA. The cingulo-opercular network provides word-recognition benefit. J Neurosci. 2013;33 (48 ):18979–18986.24285902
22 Hartwigsen  G, Bzdok  D, Klein  M, et al  Rapid short-term reorganization in the language network. Elife. 2017;6 :e25964.28537558
23 Martin  S, Williams  KA, Saur  D, Hartwigsen  G. Age-related reorganization of functional network architecture in semantic cognition. Cereb Cortex. 2023;33 (8 ):4886–4903.36190445
24 Bonkhoff  AK, Lim  JS, Bae  HJ, et al  Generative lesion pattern decomposition of cognitive impairment after stroke. Brain Commun. 2021;3 (2 ):fcab110.34189457
25 Kernbach  JM, Hartwigsen  G, Lim  JS, et al  Bayesian stroke modeling details sex biases in the white matter substrates of aphasia. Commun. Biol. 2023;6 (1 ):354.37002267
26 Polito  C, Conca  F, Santi  GC, et al  Comparing two picture naming tasks in primary progressive aphasia: Insights from behavioural and neural results. Cortex. 2023;166 :1–18.37295234
27 Hickok  G, Poeppel  D. The cortical organization of speech processing. Nat Rev. 2007;8 (5 ):393–402.
28 Camilleri  JA, Muller  VI, Fox  P, et al  Definition and characterization of an extended multiple-demand network. Neuroimage. 2018;165 :138–147.29030105
29 Duncan  J . The multiple-demand (MD) system of the primate brain: Mental programs for intelligent behaviour. Trends Cogn Sci. 2010;14 (4 ):172–179.20171926
30 Sliwinska  MW, Violante  IR, Wise  RJS, et al  Stimulating multiple-demand cortex enhances vocabulary learning. J Neurosci. 2017;37 (32 ):7606–7618.28676576
31 Jackson  RL . The neural correlates of semantic control revisited. Neuroimage. 2021;224 :117444.33059049
32 Yu  KH, Cho  SJ, Oh  MS, et al  Cognitive impairment evaluated with vascular cognitive impairment harmonization standards in a multicenter prospective stroke cohort in Korea. Stroke. 2013;44 (3 ):786–788.23271507
33 Kim  BJ, Park  JM, Kang  K, et al  Case characteristics, hyperacute treatment, and outcome information from the clinical research center for stroke-fifth division registry in South Korea. J Stroke. 2015;17 (1 ):38–53.25692106
34 Folstein  MF, Folstein  SE, McHugh  PR. “Mini-mental state”. A practical method for grading the cognitive state of patients for the clinician. J Psychiatr Res. 1975;12 (3 ):189–198.1202204
35 Kaplan  E, Goodglass  H, Weintraub  S. Boston naming test. Lea & Febiger; 1983.
36 Wechsler  D.  Manual for the Wechsler intelligence scale. Archives of Clinical Neuropsychology. Psychological Corp; 1955.
37 Reitan  RM . Validity of the trail making test as an indicator of organic brain damage. Percept Mot Skills.  1958;8 :271–276.
38 Shine  JM, Breakspear  M, Bell  PT, et al  Human cognition involves the dynamic integration of neural activity and neuromodulatory systems. Nat Neurosci. 2019;22 (2 ):289–296.30664771
39 Lee  DW, Lee  JY, Ryu  SG, et al  Validity of the Korean version of informant questionnaire on the cognitive decline in the elderly (IQCODE). J Korean Geriatr Soc. 2005;9 (3 ):196–202.
40 Lim  JS, Kim  N, Jang  MU, et al  Cortical hubs and subcortical cholinergic pathways as neural substrates of poststroke dementia. Stroke. 2014;45 (4 ):1069–1076.24603067
41 Desikan  RS, Ségonne  F, Fischl  B, et al  An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. NeuroImage. 2006;31 (3 ):968–980.16530430
42 Lee  DD, Seung  HS. Learning the parts of objects by non-negative matrix factorization. Nature. 1999;401 (6755 ):788–791.10548103
43 Goldberg  LR . The structure of phenotypic personality traits. Am Psychol. 1993;48 :26–34.8427480
44 O’Connor  BP . A quantitative review of the comprehensiveness of the five-factor model in relation to popular personality inventories. Assessment. 2002;9 (2 ):188–203.12066834
45 Gelman  A, Hill  J. Data analysis using regression and multilevel/hierarchical models. Higher Education from Cambridge University Press; 2006.
46 Bzdok  D, Varoquaux  G, Grisel  O, Eickenberg  M, Poupon  C, Thirion  B. Formal models of the network co-occurrence underlying mental operations. PLoS Comput Biol.  2016;12 (6 ):e1004994.27310288
47 Fedorenko  E, Behr  MK, Kanwisher  N. Functional specificity for high-level linguistic processing in the human brain. Proc Natl Acad Sci USA. 2011;108 (39 ):16428–16433.21885736
48 Lipkin  B, Tuckute  G, Affourtit  J, et al  Probabilistic atlas for the language network based on precision fMRI data from >800 individuals. Sci Data. 2022;9 (1 ):529.36038572
49 Markello  RD, Hansen  JY, Liu  ZQ, et al  Neuromaps: Structural and functional interpretation of brain maps. Nat Methods. 2022;19 (11 ):1472–1479.36203018
50 Kiesow  H, Uddin  LQ, Bernhardt  BC, Kable  J, Bzdok  D. Dissecting the midlife crisis: Disentangling social, personality and demographic determinants in social brain anatomy. Commun Biol. 2021;4 (1 ):728.34140617
51 Kruschke  J.  Doing Bayesian data analysis: A tutorial with R, JAGS, and Stan. Elsevier, Academic Press: 2014.
52 Baldo  JV, Arévalo  A, Patterson  JP, Dronkers  NF. Grey and white matter correlates of picture naming: Evidence from a voxel-based lesion analysis of the Boston naming test. Cortex. 2013;49 (3 ):658–667.22482693
53 Munsch  F, Sagnier  S, Asselineau  J, et al  Stroke location is an independent predictor of cognitive outcome. Stroke. 2016;47 (1 ):66–73.26585396
54 Zhao  L, Biesbroek  JM, Shi  L, et al  Strategic infarct location for post-stroke cognitive impairment: A multivariate lesion-symptom mapping study. J Cereb Blood Flow Metab. 2018;38 (8 ):1299–1311.28895445
55 Turker  S, Kuhnke  P, Eickhoff  SB, Caspers  S, Hartwigsen  G. Cortical, subcortical, and cerebellar contributions to language processing: A meta-analytic review of 403 neuroimaging experiments. Psychol Bull. 2023;149 (11-12 ):699–723.
56 Scheuringer  A, Harris  TA, Pletzer  B. Recruiting the right hemisphere: Sex differences in inter-hemispheric communication during semantic verbal fluency. Brain Lang. 2020;207 :104814.32502896
57 Wagner  S, Sebastian  A, Lieb  K, Tüscher  O, Tadić  A. A coordinate-based ALE functional MRI meta-analysis of brain activation during verbal fluency tasks in healthy control subjects. BMC Neurosci. 2014;15 :19.24456150
58 Ryu  HJ, Yang  DW. The Seoul neuropsychological screening battery (SNSB) for comprehensive neuropsychological assessment. Dement Neurocogn Disord. 2023;22 (1 ):1–15.36814700
59 Welte  PO . Indices of verbal learning and memory deficits after right hemisphere stroke. Arch Phys Med Rehabil. 1993;74 (6 ):631–636.8503754
60 Bhaya-Grossman  I, Chang  EF. Speech computations of the human superior temporal gyrus. Annu Rev Psychol. 2022;73 :79–102.34672685
61 Deschamps  I, Baum  SR, Gracco  VL. On the role of the supramarginal gyrus in phonological processing and verbal working memory: Evidence from rTMS studies. Neuropsychologia. 2014;53 :39–46.24184438
62 Hartwigsen  G, Baumgaertner  A, Price  CJ, Koehnke  M, Ulmer  S, Siebner  HR. Phonological decisions require both the left and right supramarginal gyri. Proc Natl Acad Sci U S A.  2010;107 (38 ):16494–16499.20807747
63 Grossman  M, Peelle  JE, Smith  EE, et al  Category-specific semantic memory: Converging evidence from bold fMRI and Alzheimer’s disease. Neuroimage. 2013;68 :263–274.23220494
64 Whitney  C, Kirk  M, O’Sullivan  J, Lambon Ralph  MA, Jefferies  E. The neural organization of semantic control: TMS evidence for a distributed network in left inferior frontal and posterior middle temporal gyrus. Cereb Cortex. 2011;21 (5 ):1066–1075.20851853
65 Noonan  KA, Jefferies  E, Visser  M, Ralph  L, A  M. Going beyond inferior prefrontal involvement in semantic control: Evidence for the additional contribution of dorsal angular gyrus and posterior middle temporal cortex. J Cogn Neurosci.  2013;25 (11 ):1824–1850.23859646
66 Camerino  I, Ferreira  J, Vonk  JM, et al  Systematic review and meta-analyses of word production abilities in dysfunction of the basal ganglia: Stroke, small vessel disease, Parkinson’s disease, and Huntington’s disease. Neuropsychol Rev. 2024;34 (1 ):1–26.36564612
67 Biesbroek  JM, Lim  JS, Weaver  NA, et al  Anatomy of phonemic and semantic fluency: A lesion and disconnectome study in 1231 stroke patients. Cortex. 2021;143 :148–163.34450565
68 Radanovic  M, Mansur  LL. Aphasia in vascular lesions of the basal ganglia: A comprehensive review. Brain Lang. 2017;173 :20–32.28570947
69 Nadeau  SE, Crosson  B. Subcortical aphasia. Brain Lang. 1997;58 (3 ):355–402. discussion 418–423.9222518
70 Thompson-Schill  SL, Esposito  D, Aguirre  M, Farah  GK, J  M. Role of left inferior prefrontal cortex in retrieval of semantic knowledge: A reevaluation. Proc Natl Acad Sci U S A.  1997;94 (26 ):14792–14797.9405692
71 Baron  JC, Chételat  G, Desgranges  B, et al  In vivo mapping of gray matter loss with voxel-based morphometry in mild Alzheimer’s disease. Neuroimage. 2001;14 (2 ):298–309.11467904
72 Chételat  G, Desgranges  B, Landeau  B, et al  Direct voxel-based comparison between grey matter hypometabolism and atrophy in Alzheimer’s disease. Brain. 2008;131 (Pt 1 ):60–71.18063588
73 Chen  Y, Huang  L, Chen  K, et al  White matter basis for the hub-and-spoke semantic representation: Evidence from semantic dementia. Brain. 2020;143 (4 ):1206–1219.32155237
74 Fiori  V, Kunz  L, Kuhnke  P, Marangolo  P, Hartwigsen  G. Transcranial direct current stimulation (tDCS) facilitates verb learning by altering effective connectivity in the healthy brain. Neuroimage. 2018;181 :550–559.30030198
75 Grossman  M, Carvell  S, Peltzer  L, Stern  MB, Gollomp  S, Hurtig  HI. Visual construction impairments in Parkinson’s disease. Neuropsychology. 1993;7 :536–547.
76 Weber  RC, Riccio  CA, Cohen  MJ. Does Rey complex figure copy performance measure executive function in children?  Appl Neuropsychol Child. 2013;2 (1 ):6–12.23427771
77 Leyhe  T, Saur  R, Eschweiler  GW, Milian  M. Clock test deficits are associated with semantic memory impairment in Alzheimer disease. J Geriatr Psychiatry Neurol. 2009;22 (4 ):235–245.19433861
78 McMorris  T . Exercise-cognition interaction: Neuroscience perspectives. Elsevier; 2016.
79 Avers  D, Wong  RA. Guccione’s geriatric physical therapy. 4th Edn. Elsevier; 2020.
80 Davey  J, Cornelissen  PL, Thompson  HE, et al  Automatic and controlled semantic retrieval: TMS reveals distinct contributions of posterior middle temporal gyrus and angular gyrus. J Neurosci. 2015;35 (46 ):15230–15239.26586812
81 Martin  S, Saur  D, Hartwigsen  G. Age-dependent contribution of domain-general networks to semantic cognition. Cereb Cortex. 2022;32 (4 ):870–890.34464442
82 Kircher  T, Nagels  A, Kirner-Veselinovic  A, Krach  S. Neural correlates of rhyming vs. Lexical and semantic fluency. Brain Res.  2011;1391 :71–80.21447325
83 Delis  DC, Kramer  JH, Kaplan  E, Ober  BA. California verbal learning test. 2nd edn. Psychological Corporation; 1987.
84 Seghier  ML, Patel  E, Prejawa  S, et al  The PLORAS database: A data repository for predicting language outcome and recovery after stroke. Neuroimage. 2016;124 (Pt B ):1208–1212.25882753
