Deep Variational Lesion-Deficit Mapping

Causal mapping of the functional organisation of the human brain requires evidence of \textit{necessity} available at adequate scale only from pathological lesions of natural origin. This demands inferential models with sufficient flexibility to capture both the observable distribution of pathological damage and the unobserved distribution of the neural substrate. Current model frameworks -- both mass-univariate and multivariate -- either ignore distributed lesion-deficit relations or do not model them explicitly, relying on featurization incidental to a predictive task. Here we initiate the application of deep generative neural network architectures to the task of lesion-deficit inference, formulating it as the estimation of an expressive hierarchical model of the joint lesion and deficit distributions conditioned on a latent neural substrate. We implement such deep lesion deficit inference with variational convolutional volumetric auto-encoders. We introduce a comprehensive framework for lesion-deficit model comparison, incorporating diverse candidate substrates, forms of substrate interactions, sample sizes, noise corruption, and population heterogeneity. Drawing on 5500 volume images of ischaemic stroke, we show that our model outperforms established methods by a substantial margin across all simulation scenarios, including comparatively small-scale and noisy data regimes. Our analysis justifies the widespread adoption of this approach, for which we provide an open source implementation: https://github.com/guilherme-pombo/vae_lesion_defici

The dimensionalities of lesion-deficit mapping

Lesion-deficit mapping remains the most powerful method for localising function in the human brain. As the highest court of appeal where competing theories of cerebral function conflict, it ought to be held to the most stringent inferential standards. Though at first sight elegantly transferable, the mass-univariate statistical framework popularized by functional imaging is demonstrably ill-suited to the task, both theoretically and empirically. The critical difficulty lies with the handling of the data's intrinsically high dimensionality. Conceptual opacity and computational complexity lead lesion-deficit mappers to neglect two distinct sets of anatomical interactions: those between areas unified by function, and those between areas unified by the natural pattern of pathological damage. Though both are soluble through high-dimensional multivariate analysis, the consequences of ignoring them are radically different. The former will bleach and coarsen a picture of the functional anatomy that is nonetheless broadly faithful to reality; the latter may alter it beyond all recognition. That the field continues to cling to mass-univariate methods suggests the latter problem is misidentified with the former, and that their distinction is in need of elaboration. We further argue that the vicious effects of lesion-driven interactions are not limited to anatomical localisation but will inevitably degrade purely predictive models of function such as those conceived for clinical prognostic use. Finally, we suggest there is a great deal to be learnt about lesion-mapping by simulation-based modelling of lesion data, for the fundamental problems lie upstream of the experimental data themselves

Shared lesion correlates of semantic and letter fluency in post-stroke aphasia

Lesion–symptom mapping studies have reported a temporal versus frontal dissociation between semantic and letter fluency, and mixed evidence regarding the role of white matter. Mass-univariate and multivariate lesion–symptom mapping was used to identify regions associated with semantic and letter fluency deficits in post-stroke aphasia. Multivariate LSM revealed broad networks including underlying white matter, and substantial overlap between both types of fluency, suggesting that semantic fluency and letter fluency largely rely on the same neural system. All data are available on OSF

Decoding post-stroke motor function from structural brain imaging

AbstractClinical research based on neuroimaging data has benefited from machine learning methods, which have the ability to provide individualized predictions and to account for the interaction among units of information in the brain. Application of machine learning in structural imaging to investigate diseases that involve brain injury presents an additional challenge, especially in conditions like stroke, due to the high variability across patients regarding characteristics of the lesions. Extracting data from anatomical images in a way that translates brain damage information into features to be used as input to learning algorithms is still an open question. One of the most common approaches to capture regional information from brain injury is to obtain the lesion load per region (i.e. the proportion of voxels in anatomical structures that are considered to be damaged). However, no systematic evaluation has yet been performed to compare this approach with using patterns of voxels (i.e. considering each voxel as a single feature). In this paper we compared both approaches applying Gaussian Process Regression to decode motor scores in 50 chronic stroke patients based solely on data derived from structural MRI. For both approaches we compared different ways to delimit anatomical areas: regions of interest from an anatomical atlas, the corticospinal tract, a mask obtained from fMRI analysis with a motor task in healthy controls and regions selected using lesion-symptom mapping. Our analysis showed that extracting features through patterns of voxels that represent lesion probability produced better results than quantifying the lesion load per region. In particular, from the different ways to delimit anatomical areas compared, the best performance was obtained with a combination of a range of cortical and subcortical motor areas as well as the corticospinal tract. These results will inform the appropriate methodology for predicting long term motor outcomes from early post-stroke structural brain imaging

Evaluating the granularity and statistical structure of lesions and behaviour in post-stroke aphasia

The pursuit of relating the location of neural damage to the pattern of acquired language and general cognitive deficits post-stroke stems back to 19th century behavioural neurology. Whilst spatial specificity has improved dramatically over time, from the large areas of damage specified by post-mortem investigation to the millimetre precision of modern MRI, there is an underlying issue that is rarely addressed, which relates to the fact that damage to a given area of the brain is not random but constrained by the brain’s vasculature. Accordingly, the aim of this study was to uncover the statistical structure underlying the lesion profile in chronic aphasia post-stroke. By applying varimax-rotated principal component analysis to the lesions of 70 patients with chronic post-stroke aphasia, we identified 17 interpretable clusters, largely reflecting the vascular supply of middle cerebral artery sub-branches and other sources of individual variation in vascular supply as shown in classical angiography studies. This vascular parcellation produced smaller displacement error in simulated lesion-symptom analysis compared with individual voxels and Brodmann regions. A second principal component analysis of the patients’ detailed neuropsychological data revealed a four factor solution reflecting phonological, semantic, executive-demand and speech fluency abilities. As a preliminary exploration, stepwise regression was used to relate behavioural factor scores to the lesion principal components. Phonological ability was related to two components, which covered the posterior temporal region including the posterior segment of the arcuate fasciculus, and the inferior frontal gyrus. Three components were linked to semantic ability and were located in the white matter underlying the anterior temporal lobe, the supramarginal gyrus and angular gyrus. Executive-demand related to two components covering the dorsal edge of the middle cerebral artery territory, while speech fluency was linked to two components that were located in the middle frontal gyrus, precentral gyrus, and subcortical regions (putamen and thalamus). Future studies can explore in formal terms the utility of these PCA-derived lesion components for relating post-stroke lesions and symptoms

Large-scale inference in the focally damaged human brain

Clinical outcomes in focal brain injury reflect the interactions between two distinct anatomically distributed patterns: the functional organisation of the brain and the structural distribution of injury. The challenge of understanding the functional architecture of the brain is familiar; that of understanding the lesion architecture is barely acknowledged. Yet, models of the functional consequences of focal injury are critically dependent on our knowledge of both. The studies described in this thesis seek to show how machine learning-enabled high-dimensional multivariate analysis powered by large-scale data can enhance our ability to model the relation between focal brain injury and clinical outcomes across an array of modelling applications. All studies are conducted on internationally the largest available set of MR imaging data of focal brain injury in the context of acute stroke (N=1333) and employ kernel machines at the principal modelling architecture. First, I examine lesion-deficit prediction, quantifying the ceiling on achievable predictive fidelity for high-dimensional and low-dimensional models, demonstrating the former to be substantially higher than the latter. Second, I determine the marginal value of adding unlabelled imaging data to predictive models within a semi-supervised framework, quantifying the benefit of assembling unlabelled collections of clinical imaging. Third, I compare high- and low-dimensional approaches to modelling response to therapy in two contexts: quantifying the effect of treatment at the population level (therapeutic inference) and predicting the optimal treatment in an individual patient (prescriptive inference). I demonstrate the superiority of the high-dimensional approach in both settings

Assessing and mapping language, attention and executive multidimensional deficits in stroke aphasia.

There is growing awareness that aphasia following a stroke can include deficits in other cognitive functions and that these are predictive of certain aspects of language function, recovery and rehabilitation. However, data on attentional and executive (dys)functions in individuals with stroke aphasia are still scarce and the relationship to underlying lesions is rarely explored. Accordingly in this investigation, an extensive selection of standardized non-verbal neuropsychological tests was administered to 38 individuals with chronic post-stroke aphasia, in addition to detailed language testing and MRI. To establish the core components underlying the variable patients' performance, behavioural data were explored with rotated principal component analyses, first separately for the non-verbal and language tests, then in a combined analysis including all tests. Three orthogonal components for the non-verbal tests were extracted, which were interpreted as shift-update, inhibit-generate and speed. Three components were also extracted for the language tests, representing phonology, semantics and speech quanta. Individual continuous scores on each component were then included in a voxel-based correlational methodology analysis, yielding significant clusters for all components. The shift-update component was associated with a posterior left temporo-occipital and bilateral medial parietal cluster, the inhibit-generate component was mainly associated with left frontal and bilateral medial frontal regions, and the speed component with several small right-sided fronto-parieto-occipital clusters. Two complementary multivariate brain-behaviour mapping methods were also used, which showed converging results. Together the results suggest that a range of brain regions are involved in attention and executive functioning, and that these non-language domains play a role in the abilities of patients with chronic aphasia. In conclusion, our findings confirm and extend our understanding of the multidimensionality of stroke aphasia, emphasize the importance of assessing non-verbal cognition in this patient group and provide directions for future research and clinical practice. We also briefly compare and discuss univariate and multivariate methods for brain-behaviour mapping

Metabolic lesion-deficit mapping of human cognition

In theory the most powerful technique for functional localization in cognitive neuroscience, lesion-deficit mapping is in practice distorted by unmodelled network disconnections and strong ‘parasitic’ dependencies between collaterally damaged ischaemic areas. High-dimensional multivariate modelling can overcome these defects, but only at the cost of commonly impracticable data scales. Here we develop lesion-deficit mapping with metabolic lesions—discrete areas of hypometabolism typically seen on interictal 18F-fluorodeoxyglucose PET imaging in patients with focal epilepsy—that inherently capture disconnection effects, and whose structural dependence patterns are sufficiently benign to allow the derivation of robust functional anatomical maps with modest data. In this cross-sectional study of 159 patients with widely distributed focal cortical impairments, we derive lesion-deficit maps of a broad range of psychological subdomains underlying affect and cognition. We demonstrate the potential clinical utility of the approach in guiding therapeutic resection for focal epilepsy or other neurosurgical indications by applying high-dimensional modelling to predict out-of-sample verbal IQ and depression from cortical metabolism alone