On imputing function to structure from the behavioural effects of brain lesions

Flexibility of in vitro cortical circuits influences resilience from microtrauma

Background

Small clusters comprising hundreds to thousands of neurons are an important level of brain architecture that correlates single neuronal properties to fulfill brain function, but the specific mechanisms through which this scaling occurs are not well understood. In this study, we developed an in vitro experimental platform of small neuronal circuits (islands) to probe the importance of structural properties for their development, physiology, and response to microtrauma.

Methods

Primary cortical neurons were plated on a substrate patterned to promote attachment in clusters of hundreds of cells (islands), transduced with GCaMP6f, allowed to mature until 10-13 days in vitro (DIV), and monitored with Ca²⁺ as a non-invasive proxy for electrical activity. We adjusted two structural factors-island size and cellular density-to evaluate their role in guiding spontaneous activity and network formation in neuronal islands.

Results

We found cellular density, but not island size, regulates of circuit activity and network function in this system. Low cellular density islands can achieve many states of activity, while high cellular density biases islands towards a limited regime characterized by low rates of activity and high synchronization, a property we summarized as "flexibility." The injury severity required for an island to lose activity in 50% of its population was significantly higher in low-density, high flexibility islands.

Conclusion

Together, these studies demonstrate flexible living cortical circuits are more resilient to microtrauma, providing the first evidence that initial circuit state may be a key factor to consider when evaluating the consequences of trauma to the cortex.

Systematic perturbation of an artificial neural network: A step towards quantifying causal contributions in the brain

Lesion inference analysis is a fundamental approach for characterizing the causal contributions of neural elements to brain function. This approach has gained new prominence through the arrival of modern perturbation techniques with unprecedented levels of spatiotemporal precision. While inferences drawn from brain perturbations are conceptually powerful, they face methodological difficulties. Particularly, they are challenged to disentangle the true causal contributions of the involved elements, since often functions arise from coalitions of distributed, interacting elements, and localized perturbations have unknown global consequences. To elucidate these limitations, we systematically and exhaustively lesioned a small artificial neural network (ANN) playing a classic arcade game. We determined the functional contributions of all nodes and links, contrasting results from sequential single-element perturbations with simultaneous perturbations of multiple elements. We found that lesioning individual elements, one at a time, produced biased results. By contrast, multi-site lesion analysis captured crucial details that were missed by single-site lesions. We conclude that even small and seemingly simple ANNs show surprising complexity that needs to be addressed by multi-lesioning for a coherent causal characterization.

The motto “No causation without manipulation” is canonical to scientific endeavors. In particular, neuroscience seeks to identify which brain elements are causally involved in cognition and behavior, by perturbing them. However, due to multi-dimensional interactions among the elements, this goal has remained challenging. Here, we used an Artificial Neural Network as a ground-truth model to compare the inferential capacities of two principal approaches, lesioning a system one element at a time versus sampling from the set of all possible combinations of lesions. We show that lesioning one element at a time provides misleading results. Hence, we argue for employing exhaustive perturbation regimes. We further advocate using simulation experiments and ground-truth models to verify the assumptions and limitations of current approaches for brain mapping by perturbation.

Multiple timescales of sensory-evidence accumulation across the dorsal cortex

Cortical areas seem to form a hierarchy of intrinsic timescales, but the relevance of this organization for cognitive behavior remains unknown. In particular, decisions requiring the gradual accrual of sensory evidence over time recruit widespread areas across this hierarchy. Here, we tested the hypothesis that this recruitment is related to the intrinsic integration timescales of these widespread areas. We trained mice to accumulate evidence over seconds while navigating in virtual reality and optogenetically silenced the activity of many cortical areas during different brief trial epochs. We found that the inactivation of all tested areas affected the evidence-accumulation computation. Specifically, we observed distinct changes in the weighting of sensory evidence occurring during and before silencing, such that frontal inactivations led to stronger deficits on long timescales than posterior cortical ones. Inactivation of a subset of frontal areas also led to moderate effects on behavioral processes beyond evidence accumulation. Moreover, large-scale cortical Ca2+ activity during task performance displayed different temporal integration windows. Our findings suggest that the intrinsic timescale hierarchy of distributed cortical areas is an important component of evidence-accumulation mechanisms.

Preparing for a Second Attack: A Lesion Simulation Study on Network Resilience After Stroke

Background:

Does the brain become more resilient after a first stroke to reduce the consequences of a new lesion? Although recurrent strokes are a major clinical issue, whether and how the brain prepares for a second attack is unknown. This is due to the difficulties to obtain an appropriate dataset of stroke patients with comparable lesions, imaged at the same interval after onset. Furthermore, timing of the recurrent event remains unpredictable.

Methods:

Here, we used a novel clinical lesion simulation approach to test the hypothesis that resilience in brain networks increases during stroke recovery. Sixteen highly selected patients with a lesion restricted to the primary motor cortex were recruited. At 3 time points of the index event (10 days, 3 weeks, 3 months), we mimicked recurrent infarcts by deletion of nodes in brain networks (resting-state functional magnetic resonance imaging). Graph measures were applied to determine resilience (global efficiency after attack) and wiring cost (mean degree) of the network.

Results:

At 10 days and 3 weeks after stroke, resilience was similar in patients and controls. However, at 3 months, although motor function had fully recovered, resilience to clinically representative simulated lesions was higher compared to controls (cortical lesion P =0.012; subcortical: P =0.009; cortico-subcortical: P =0.009). Similar results were found after random ( P =0.012) and targeted ( P =0.015) attacks.

Conclusions:

Our results suggest that, in this highly selected cohort of patients with lesions restricted to the primary motor cortex, brain networks reconfigure to increase resilience to future insults. Lesion simulation is an innovative approach, which may have major implications for stroke therapy. Individualized neuromodulation strategies could be developed to foster resilient network reconfigurations after a first stroke to limit the consequences of future attacks.

The Resilience and Vulnerability of Human Brain Networks Across the Lifespan

Resilience, the ability for a system to maintain its basic functionality when suffering from lesions, is a critical property for human brain, especially in the brain aging process. This study adopted a novel metric of network resilience, the Resilience Index (RI), to assess human brain resilience with three different lifespan datasets. Based on the structural brain networks constructed from diffusion tensor imaging (DTI), we observed an inverted-U relationship between RI and age, that is, RI increased during development and early adulthood, reached a peak at about 35 years old, and then decreased during aging, which suggested that brain resilience could be quantified by RI. Furthermore, we studied brain network vulnerability by the decreases in RI when virtual lesions occurred to nodes (i.e., brain regions) or edges (i.e., structural brain connectivity). We found that the strong edges were markedly vulnerable, and the homotopic edges were the most prominent representatives of vulnerable edges. In other words, an arbitrary attack on homotopic edges would have a high probability to degrade brain network resilience. These findings suggest the change of human brain resilience across the lifespan and provide a new perspective for exploring human brain vulnerability.

Investigating the network consequences of focal brain lesions through comparisons of real and simulated lesions

Given the increased interest in the functional human connectome, a number of computer simulation studies have sought to develop a better quantitative understanding of the effects of focal lesions on the brain’s functional network organization. However, there has been little work evaluating the predictions of this simulation work vis a vis real lesioned connectomes. One of the few relevant studies reported findings from real chronic focal lesions that only partially confirmed simulation predictions. We hypothesize that these discrepancies arose because although the effects of focal lesions likely consist of two components: short-term node subtraction and long-term network re-organization, previous simulation studies have primarily modeled only the short-term consequences of the subtraction of lesioned nodes and their connections. To evaluate this hypothesis, we compared network properties (modularity, participation coefficient, within-module degree) between real functional connectomes obtained from chronic stroke participants and “pseudo-lesioned” functional connectomes generated by subtracting the same sets of lesioned nodes/connections from healthy control connectomes. We found that, as we hypothesized, the network properties of real-lesioned connectomes in chronic stroke differed from those of the pseudo-lesioned connectomes which instantiated only the short-term consequences of node subtraction. Reflecting the long-term consequences of focal lesions, we found re-organization of the neurotopography of global and local hubs in the real but not the pseudo-lesioned connectomes. We conclude that the long-term network re-organization that occurs in response to focal lesions involves changes in functional connectivity within the remaining intact neural tissue that go well beyond the short-term consequences of node subtraction.

Behavioral, Physiological and EEG Activities Associated with Conditioned Fear as Sensors for Fear and Anxiety

Anxiety disorders impose substantial costs upon public health and productivity in the USA and worldwide. At present, these conditions are quantified by self-report questionnaires that only apply to behaviors that are accessible to consciousness, or by the timing of responses to fear- and anxiety-related words that are indirect since they do not produce fear, e.g., Dot Probe Test and emotional Stroop. We now review the conditioned responses (CRs) to fear produced by a neutral stimulus (conditioned stimulus CS+) when it cues a painful laser unconditioned stimulus (US). These CRs include autonomic (Skin Conductance Response) and ratings of the CS+ unpleasantness, ability to command attention, and the recognition of the association of CS+ with US (expectancy). These CRs are directly related to fear, and some measure behaviors that are minimally accessible to consciousness e.g., economic scales. Fear-related CRs include non-phase-locked phase changes in oscillatory EEG power defined by frequency and time post-stimulus over baseline, and changes in phase-locked visual and laser evoked responses both of which include late potentials reflecting attention or expectancy, like the P300, or contingent negative variation. Increases (ERS) and decreases (ERD) in oscillatory power post-stimulus may be generalizable given their consistency across healthy subjects. ERS and ERD are related to the ratings above as well as to anxious personalities and clinical anxiety and can resolve activity over short time intervals like those for some moods and emotions. These results could be incorporated into an objective instrumented test that measures EEG and CRs of autonomic activity and psychological ratings related to conditioned fear, some of which are subliminal. As in the case of instrumented tests of vigilance, these results could be useful for the direct, objective measurement of multiple aspects of the risk, diagnosis, and monitoring of therapies for anxiety disorders and anxious personalities.

Task-Dependent Changes in the Large-Scale Dynamics and Necessity of Cortical Regions

Neural activity throughout the cortex is correlated with perceptual decisions, but inactivation studies suggest that only a small number of areas are necessary for these behaviors. Here we show that the number of required cortical areas and their dynamics vary across related tasks with different cognitive computations. In a visually guided virtual T-maze task, bilateral inactivation of only a few dorsal cortical regions impaired performance. In contrast, in tasks requiring evidence accumulation and/or post-stimulus memory, performance was impaired by inactivation of widespread cortical areas with diverse patterns of behavioral deficits across areas and tasks. Wide-field imaging revealed widespread ramps of Ca2+ activity during the accumulation and visually guided tasks. Additionally, during accumulation, different regions had more diverse activity profiles, leading to reduced inter-area correlations. Using a modular recurrent neural network model trained to perform analogous tasks, we argue that differences in computational strategies alone could explain these findings.

Network curvature as a hallmark of brain structural connectivity

Although brain functionality is often remarkably robust to lesions and other insults, it may be fragile when these take place in specific locations. Previous attempts to quantify robustness and fragility sought to understand how the functional connectivity of brain networks is affected by structural changes, using either model-based predictions or empirical studies of the effects of lesions. We advance a geometric viewpoint relying on a notion of network curvature, the so-called Ollivier-Ricci curvature. This approach has been proposed to assess financial market robustness and to differentiate biological networks of cancer cells from healthy ones. Here, we apply curvature-based measures to brain structural networks to identify robust and fragile brain regions in healthy subjects. We show that curvature can also be used to track changes in brain connectivity related to age and autism spectrum disorder (ASD), and we obtain results that are in agreement with previous MRI studies.

The brain can often continue to function despite lesions in many areas, but damage to particular locations may have serious effects. Here, the authors use the concept of Ollivier-Ricci curvature to investigate the robustness of brain networks.

The impact of ischemic stroke on connectivity gradients

The functional organization of the brain can be represented as a low-dimensional space that reflects its macroscale hierarchy. The dimensions of this space, described as connectivity gradients, capture the similarity of areas' connections along a continuous space. Studying how pathological perturbations with known effects on functional connectivity affect these connectivity gradients provides support for their biological relevance. Previous work has shown that localized lesions cause widespread functional connectivity alterations in structurally intact areas, affecting a network of interconnected regions. By using acute stroke as a model of the effects of focal lesions on the connectome, we apply the connectivity gradient framework to depict how functional reorganization occurs throughout the brain, unrestricted by traditional definitions of functional network boundaries. We define a three-dimensional connectivity space template based on functional connectivity data from healthy controls. By projecting lesion locations into this space, we demonstrate that ischemic strokes result in dimension-specific alterations in functional connectivity over the first week after symptom onset. Specifically, changes in functional connectivity were captured along connectivity Gradients 1 and 3. The degree of functional connectivity change was associated with the distance from the lesion along these connectivity gradients (a measure of functional similarity) regardless of the anatomical distance from the lesion. Together, these results provide support for the biological validity of connectivity gradients and suggest a novel framework to characterize connectivity alterations after stroke.

Perturbation-driven paradoxical facilitation of visuo-spatial function: Revisiting the 'Sprague effect'

The 'Sprague Effect' described in the seminal paper of James Sprague (Science 153:1544-1547, 1966a) is an unexpected paradoxical effect in which a second brain lesion reversed functional deficits induced by an earlier lesion. It was observed initially in the cat where severe and permanent contralateral visually guided attentional deficits generated by the ablation of large areas of the visual cortex were reversed by the subsequent removal of the superior colliculus (SC) opposite to the cortical lesion or by the splitting of the collicular commissure. Physiologically, this effect has been explained in several ways-most notably by the reduction of the functional inhibition of the ipsilateral SC by the contralateral SC, and the restoration of normal interactions between cortical and midbrain structures after ablation. In the present review, we aim at reappraising the 'Sprague Effect' by critically analyzing studies that have been conducted in the feline and human brain. Moreover, we assess applications of the 'Sprague Effect' in the rehabilitation of visually guided attentional impairments by using non-invasive therapeutic approaches such as transcranial magnetic stimulation (TMS) and transcranial direct-current stimulation (tDCS). We also review theoretical models of the effect that emphasize the inhibition and balancing between the two hemispheres and show implications for lesion inference approaches. Last, we critically review whether the resulting inter-hemispheric rivalry theories lead toward an efficient rehabilitation of stroke in humans. We conclude by emphasizing key challenges in the field of 'Sprague Effect' applications in order to design better therapies for brain-damaged patients.

Structural connectivity in multiple sclerosis and modeling of disconnection

BACKGROUND

Multiple sclerosis (MS) is characterized by focal white matter damage, and when the brain is modeled as a network, lesions can be treated as disconnection events.

OBJECTIVE

To evaluate whether modeling disconnection caused by lesions helps explain motor and cognitive impairment in MS.

METHODS

Pathways connecting 116 cortical regions were reconstructed with magnetic resonance imaging (MRI) tractography from diffusion tensors averaged across healthy controls (HCs); maps of pathways were applied to 227 relapse-onset MS patients and 50 HCs to derive structural connectivity. Then, the likelihood of individual connections passing through lesions was used to model disconnection. Patients were grouped according to clinical phenotype (113 relapsing-remitting multiple sclerosis (RRMS), 69 secondary progressive multiple sclerosis (SPMS), 45 benign MS), and then network metrics were compared between groups (analysis of variance (ANOVA)) and correlated with motor and cognitive scores (linear regression).

RESULTS

Global metrics differentiated RRMS from SPMS and benign MS patients, but not benign from SPMS patients. Nodal connectivity strength replicated global results. After disconnection, few nodes were significantly different between benign MS and RRMS patients. Correlations revealed nodes pertinent to motor and cognitive dysfunctions; these became slightly stronger after disconnection.

CONCLUSION

Connectivity did not change greatly after modeled disconnection, suggesting that the brain network is robust against damage caused by MS lesions.

The search of "canonical" explanations for the cerebral cortex

This paper addresses a fundamental line of research in neuroscience: the identification of a putative neural processing core of the cerebral cortex, often claimed to be "canonical". This "canonical" core would be shared by the entire cortex, and would explain why it is so powerful and diversified in tasks and functions, yet so uniform in architecture. The purpose of this paper is to analyze the search for canonical explanations over the past 40 years, discussing the theoretical frameworks informing this research. It will highlight a bias that, in my opinion, has limited the success of this research project, that of overlooking the dimension of cortical development. The earliest explanation of the cerebral cortex as canonical was attempted by David Marr, deriving putative cortical circuits from general mathematical laws, loosely following a deductive-nomological account. Although Marr's theory turned out to be incorrect, one of its merits was to have put the issue of cortical circuit development at the top of his agenda. This aspect has been largely neglected in much of the research on canonical models that has followed. Models proposed in the 1980s were conceived as mechanistic. They identified a small number of components that interacted as a basic circuit, with each component defined as a function. More recent models have been presented as idealized canonical computations, distinct from mechanistic explanations, due to the lack of identifiable cortical components. Currently, the entire enterprise of coming up with a single canonical explanation has been criticized as being misguided, and the premise of the uniformity of the cortex has been strongly challenged. This debate is analyzed here. The legacy of the canonical circuit concept is reflected in both positive and negative ways in recent large-scale brain projects, such as the Human Brain Project. One positive aspect is that these projects might achieve the aim of producing detailed simulations of cortical electrical activity, a negative one regards whether they will be able to find ways of simulating how circuits actually develop.

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.

Brain networks under attack: robustness properties and the impact of lesions

A growing number of studies approach the brain as a complex network, the so-called 'connectome'. Adopting this framework, we examine what types or extent of damage the brain can withstand-referred to as network 'robustness'-and conversely, which kind of distortions can be expected after brain lesions. To this end, we review computational lesion studies and empirical studies investigating network alterations in brain tumour, stroke and traumatic brain injury patients. Common to these three types of focal injury is that there is no unequivocal relationship between the anatomical lesion site and its topological characteristics within the brain network. Furthermore, large-scale network effects of these focal lesions are compared to those of a widely studied multifocal neurodegenerative disorder, Alzheimer's disease, in which central parts of the connectome are preferentially affected. Results indicate that human brain networks are remarkably resilient to different types of lesions, compared to other types of complex networks such as random or scale-free networks. However, lesion effects have been found to depend critically on the topological position of the lesion. In particular, damage to network hub regions-and especially those connecting different subnetworks-was found to cause the largest disturbances in network organization. Regardless of lesion location, evidence from empirical and computational lesion studies shows that lesions cause significant alterations in global network topology. The direction of these changes though remains to be elucidated. Encouragingly, both empirical and modelling studies have indicated that after focal damage, the connectome carries the potential to recover at least to some extent, with normalization of graph metrics being related to improved behavioural and cognitive functioning. To conclude, we highlight possible clinical implications of these findings, point out several methodological limitations that pertain to the study of brain diseases adopting a network approach, and provide suggestions for future research.

Technical considerations of a game-theoretical approach for lesion symptom mapping

Background

Various strategies have been used for inferring brain functions from stroke lesions. We explored a new mathematical approach based on game theory, the so-called multi-perturbation Shapley value analysis (MSA), to assess causal function localizations and interactions from multiple perturbation data. We applied MSA to a dataset composed of lesion patterns of 148 acute stroke patients and their National Institutes of Health Stroke Scale (NIHSS) scores, to systematically investigate the influence of different parameter settings on the outcomes of the approach. Specifically, we investigated aspects of MSA methodology including the choice of the predictor algorithm (typology and kernel functions), training dataset (original versus binary), as well as the influence of lesion thresholds. We assessed the suitability of MSA for processing real clinical lesion data and established the central parameters for this analysis.

Results

We derived general recommendations for the analysis of clinical datasets by MSA and showed that, for the studied dataset, the best approach was to use a linear-kernel support vector machine predictor, trained with a binary training dataset, where the binarization was implemented through a median threshold of lesion size for each region. We demonstrated that the results obtained with different MSA variants lead to almost identical results as the basic MSA.

Conclusions

MSA is a feasible approach for the multivariate lesion analysis of clinical stroke data. Informed choices need to be made to set parameters that may affect the analysis outcome.

Connectivity-based parcellation: Critique and implications

Regional specialization and functional integration are often viewed as two fundamental principles of human brain organization. They are closely intertwined because each functionally specialized brain region is probably characterized by a distinct set of long-range connections. This notion has prompted the quickly developing family of connectivity-based parcellation (CBP) methods in neuroimaging research. CBP assumes that there is a latent structure of parcels in a region of interest (ROI). First, connectivity strengths are computed to other parts of the brain for each voxel/vertex within the ROI. These features are then used to identify functionally distinct groups of ROI voxels/vertices. CBP enjoys increasing popularity for the in-vivo mapping of regional specialization in the human brain. Due to the requirements of different applications and datasets, CBP has diverged into a heterogeneous family of methods. This broad overview critically discusses the current state as well as the commonalities and idiosyncrasies of the main CBP methods. We target frequent concerns faced by novices and veterans to provide a reference for the investigation and review of CBP studies.

Causal functional contributions and interactions in the attention network of the brain: an objective multi-perturbation analysis

Spatial attention is a prime example for the distributed network functions of the brain. Lesion studies in animal models have been used to investigate intact attentional mechanisms as well as perspectives for rehabilitation in the injured brain. Here, we systematically analyzed behavioral data from cooling deactivation and permanent lesion experiments in the cat, where unilateral deactivation of the posterior parietal cortex (in the vicinity of the posterior middle suprasylvian cortex, pMS) or the superior colliculus (SC) cause a severe neglect in the contralateral hemifield. Counterintuitively, additional deactivation of structures in the opposite hemisphere reverses the deficit. Using such lesion data, we employed a game-theoretical approach, multi-perturbation Shapley value analysis (MSA), for inferring functional contributions and network interactions of bilateral pMS and SC from behavioral performance in visual attention studies. The approach provides an objective theoretical strategy for lesion inferences and allows a unique quantitative characterization of regional functional contributions and interactions on the basis of multi-perturbations. The quantitative analysis demonstrated that right posterior parietal cortex and superior colliculus made the strongest positive contributions to left-field orienting, while left brain regions had negative contributions, implying that their perturbation may reverse the effects of contralateral lesions or improve normal function. An analysis of functional modulations and interactions among the regions revealed redundant interactions (implying functional overlap) between regions within each hemisphere, and synergistic interactions between bilateral regions. To assess the reliability of the MSA method in the face of variable and incomplete input data, we performed a sensitivity analysis, investigating how much the contribution values of the four regions depended on the performance of specific configurations and on the prediction of unknown performances. The results suggest that the MSA approach is sensitive to categorical, but insensitive to gradual changes in the input data. Finally, we created a basic network model that was based on the known anatomical interactions among cortical-tectal regions and reproduced the experimentally observed behavior in visual orienting. We discuss the structural organization of the network model relative to the causal modulations identified by MSA, to aid a mechanistic understanding of the attention network of the brain.

Why is it difficult to predict language impairment and outcome in patients with aphasia after stroke?

One of the most devastating consequences of stroke is aphasia. Communication problems after stroke can severely impair the patient's quality of life and make even simple everyday tasks challenging. Despite intense research in the field of aphasiology, the type of language impairment has not yet been localized and correlated with brain damage, making it difficult to predict the language outcome for stroke patients with aphasia. Our primary objective is to present the available evidence that highlights the difficulties of predicting language impairment after stroke. The different levels of complexity involved in predicting the lesion site from language impairment and ultimately predicting the long-term outcome in stroke patients with aphasia were explored. Future directions and potential implications for research and clinical practice are highlighted.

Comparative analysis of the macroscale structural connectivity in the macaque and human brain

The macaque brain serves as a model for the human brain, but its suitability is challenged by unique human features, including connectivity reconfigurations, which emerged during primate evolution. We perform a quantitative comparative analysis of the whole brain macroscale structural connectivity of the two species. Our findings suggest that the human and macaque brain as a whole are similarly wired. A region-wise analysis reveals many interspecies similarities of connectivity patterns, but also lack thereof, primarily involving cingulate regions. We unravel a common structural backbone in both species involving a highly overlapping set of regions. This structural backbone, important for mediating information across the brain, seems to constitute a feature of the primate brain persevering evolution. Our findings illustrate novel evolutionary aspects at the macroscale connectivity level and offer a quantitative translational bridge between macaque and human research.

What are the commonalities and differences of human brains when compared to the brains of other primates? The brain can be conceived as a complex network. Its topological properties constrain its function. Ethical and technical reasons necessitate the use of animal brains, like the macaque monkey, as models for the human brain. However, evolutionary changes, including “brain rewiring”, might result in unique human features. Hence, a detailed and quantitative comparative analysis of the connectivity of the brains of the two species is needed. Here, we undertake this task by adopting techniques analogous to those used in comparative studies in other scientific fields. Our approach reveals converging but also diverging wiring patterns. The brain of the two species as a whole is similarly wired. The majority of the brain regions appear to have evolutionary conserved connectivity patterns while for certain regions this appears not to be the case. We also uncover an evolutionary conserved “structural backbone” in the brain of the two species. Our findings highlight common and unique “wiring properties” of the brains of these two primate species and offer a quantitative basis for translating findings from macaque research to human research.

Methods of association and dissociation for establishing selective brain-behavior relations

Methods for identifying and understanding brain structure-function relations have evolved over the past century, from astute observations of selective impairments associated with focal brain damage to dissociations measured by combining quantitative neuropsychologic assessment and brain imaging. Enhanced spatial and temporal resolution in brain imaging modalities has led to refined visualization and quantification of the brain's substructures, microstructural integrity, and functional connectivity of neural networks. The double dissociation model has been a gold standard used to demonstrate that a particular cognitive, emotional, sensory, or motor process is selectively related to a particular brain region or neural network and not to others. This model has provided a fruitful means for testing hypotheses of functional localization and enabled examination and establishment of component processes contributing to complex cognitive and motor functions, parsing multifactorial behaviors and identifying brain regions, and networks subserving these complex abilities. In this chapter we discuss the evolution of the dissociation model and highlight how the modifications of this model are used presently to establish selective brain-behavior relationships in disorders such as chronic alcoholism with a neuropathologic signature but no localizable, space-occupying lesion.

Default Mode Network Connectivity in Stroke Patients

The pathophysiology of episodic memory dysfunction after infarction is not completely understood. It has been suggested that infarctions located anywhere in the brain can induce widespread effects causing disruption of functional networks of the cortical regions. The default mode network, which includes the medial temporal lobe, is a functional network that is associated with episodic memory processing. We investigated whether the default mode network activity is reduced in stroke patients compared to healthy control subjects in the resting state condition. We assessed the whole brain network properties during resting state functional MRI in 21 control subjects and 20 ‘first-ever’ stroke patients. Patients were scanned 9–12 weeks after stroke onset. Stroke lesions were located in various parts of the brain. Independent component analyses were conducted to identify the default mode network and to compare the group differences of the default mode network. Furthermore, region-of-interest based analysis was performed to explore the functional connectivity between the regions of the default mode network. Stroke patients performed significantly worse than control subjects on the delayed recall score on California verbal learning test. We found decreased functional connectivity in the left medial temporal lobe, posterior cingulate and medial prefrontal cortical areas within the default mode network and reduced functional connectivity between these regions in stroke patients compared with controls. There were no significant volumetric differences between the groups. These results demonstrate that connectivity within the default mode network is reduced in ‘first-ever’ stroke patients compared to control subjects. This phenomenon might explain the occurrence of post-stroke cognitive dysfunction in stroke patients.

Context-dependent interactions of left posterior inferior frontal gyrus in a local visual search task unrelated to language

The Embedded Figures Task (EFT) involves search for a target hidden in a complex geometric pattern. Even though the EFT is designed to probe local visual search functions, not language-related processes, neuropsychological studies have demonstrated a strong association between aphasia and impairment on this task. A potential explanation for this relationship was offered by a recent functional MRI study (Manjaly et al., 2003), which demonstrated that a part of the left posterior inferior frontal gyrus (pIFG), overlapping with Broca's region, is crucially involved in the execution of the EFT. This result suggested that pIFG, an area strongly associated with language-related functions, is also part of a network subserving cognitive functions unrelated to language. In this study, we tested this conjecture by analysing the data of Manjaly et al. for context-dependent functional interactions of the pIFG during execution of the EFT. The results showed that during EFT, compared to a similar visual matching task with minimal local search components, pIFG changed its interactions with areas commonly involved in visuospatial processing: Increased contributions to neural activity in left posterior parietal cortex, cerebellar vermis, and extrastriate areas bilaterally, as well as decreased contributions to bilateral temporo-parietal cortex, posterior cingulate cortex, and left dorsal premotor cortex were found. These findings demonstrate that left pIFG can be involved in nonlanguage processes. More generally, however, they provide a concrete example of the notion that there is no general one-to-one mapping between cognitive functions and the activations of individual areas. Instead, it is the spatiotemporal pattern of functional interactions between areas that is linked to a particular cognitive context.

Neuroimaging in aphasia treatment research: quantifying brain lesions after stroke

New structural and functional neuroimaging methods continue to rapidly develop, offering promising tools for cognitive neuroscientists. In the last 20 years, advanced magnetic resonance imaging (MRI) techniques have provided invaluable insights into how language is represented and processed in the brain and how it can be disrupted by damage to, or dysfunction of, various parts of the brain. Current functional MRI (fMRI) approaches have also allowed researchers to purposefully investigate how individuals recover language after stroke. This paper presents recommendations for quantification of brain lesions derived from discussions among international researchers at the Neuroimaging in Aphasia Treatment Research Workshop held at Northwestern University (Evanston, Illinois, USA). Methods for detailing and characterizing the brain damage that can influence results of fMRI studies in chronic aphasic stroke patients are discussed. Moreover, we aimed to provide the reader with a set of general practical guidelines and references to facilitate choosing adequate structural imaging strategies that facilitate fMRI studies in aphasia treatment research.

Focal brain lesions to critical locations cause widespread disruption of the modular organization of the brain

Although it is generally assumed that brain damage predominantly affects only the function of the damaged region, here we show that focal damage to critical locations causes disruption of network organization throughout the brain. Using resting state fMRI, we assessed whole-brain network structure in patients with focal brain lesions. Only damage to those brain regions important for communication between subnetworks (e.g., "connectors")--but not to those brain regions important for communication within sub-networks (e.g., "hubs")--led to decreases in modularity, a measure of the integrity of network organization. Critically, this network dysfunction extended into the structurally intact hemisphere. Thus, focal brain damage can have a widespread, nonlocal impact on brain network organization when there is damage to regions important for the communication between networks. These findings fundamentally revise our understanding of the remote effects of focal brain damage and may explain numerous puzzling cases of functional deficits that are observed following brain injury.

Predicting language outcome and recovery after stroke: the PLORAS system

The ability to comprehend and produce speech after stroke depends on whether the areas of the brain that support language have been damaged. Here, we review two different ways to predict language outcome after stroke. The first depends on understanding the neural circuits that support language. This model-based approach is a challenging endeavor because language is a complex cognitive function that involves the interaction of many different brain areas. The second approach, by contrast, does not require an understanding of why a lesion impairs language; instead, predictions are made on the basis of the recovery of previous patients with the same lesion. This approach requires a database that records the speech and language capabilities of a large population of patients who have, collectively, incurred a comprehensive range of focal brain lesions. In addition, a system is required that converts an MRI scan from a new patient into a three-dimensional description of the lesion and compares this lesion against all others on the database. The outputs of this system are the longitudinal language outcomes of corresponding patients in the database. This approach will provide the patient with a range of probable recovery patterns over a variety of language measures.

Cortical hubs form a module for multisensory integration on top of the hierarchy of cortical networks

Sensory stimuli entering the nervous system follow particular paths of processing, typically separated (segregated) from the paths of other modal information. However, sensory perception, awareness and cognition emerge from the combination of information (integration). The corticocortical networks of cats and macaque monkeys display three prominent characteristics: (i) modular organisation (facilitating the segregation), (ii) abundant alternative processing paths and (iii) the presence of highly connected hubs. Here, we study in detail the organisation and potential function of the cortical hubs by graph analysis and information theoretical methods. We find that the cortical hubs form a spatially delocalised, but topologically central module with the capacity to integrate multisensory information in a collaborative manner. With this, we resolve the underlying anatomical substrate that supports the simultaneous capacity of the cortex to segregate and to integrate multisensory information.

Can structure predict function in the human brain?

Over the past decade, scientific interest in the properties of large-scale spontaneous neural dynamics has intensified. Concurrently, novel technologies have been developed for characterizing the connective anatomy of intra-regional circuits and inter-regional fiber pathways. It will soon be possible to build computational models that incorporate these newly detailed structural network measurements to make predictions of neural dynamics at multiple scales. Here, we review the practicality and the value of these efforts, while at the same time considering in which cases and to what extent structure does determine neural function. Studies of the healthy brain, of neural development, and of pathology all yield examples of direct correspondences between structural linkage and dynamical correlation. Theoretical arguments further support the notion that brain network topology and spatial embedding should strongly influence network dynamics. Although future models will need to be tested more quantitatively and against a wider range of empirical neurodynamic features, our present large-scale models can already predict the macroscopic pattern of dynamic correlation across the brain. We conclude that as neuroscience grapples with datasets of increasing completeness and complexity, and attempts to relate the structural and functional architectures discovered at different neural scales, the value of computational modeling will continue to grow.

Modeling the impact of lesions in the human brain

Lesions of anatomical brain networks result in functional disturbances of brain systems and behavior which depend sensitively, often unpredictably, on the lesion site. The availability of whole-brain maps of structural connections within the human cerebrum and our increased understanding of the physiology and large-scale dynamics of cortical networks allow us to investigate the functional consequences of focal brain lesions in a computational model. We simulate the dynamic effects of lesions placed in different regions of the cerebral cortex by recording changes in the pattern of endogenous ("resting-state") neural activity. We find that lesions produce specific patterns of altered functional connectivity among distant regions of cortex, often affecting both cortical hemispheres. The magnitude of these dynamic effects depends on the lesion location and is partly predicted by structural network properties of the lesion site. In the model, lesions along the cortical midline and in the vicinity of the temporo-parietal junction result in large and widely distributed changes in functional connectivity, while lesions of primary sensory or motor regions remain more localized. The model suggests that dynamic lesion effects can be predicted on the basis of specific network measures of structural brain networks and that these effects may be related to known behavioral and cognitive consequences of brain lesions.

Organization of excitable dynamics in hierarchical biological networks

This study investigates the contributions of network topology features to the dynamic behavior of hierarchically organized excitable networks. Representatives of different types of hierarchical networks as well as two biological neural networks are explored with a three-state model of node activation for systematically varying levels of random background network stimulation. The results demonstrate that two principal topological aspects of hierarchical networks, node centrality and network modularity, correlate with the network activity patterns at different levels of spontaneous network activation. The approach also shows that the dynamic behavior of the cerebral cortical systems network in the cat is dominated by the network's modular organization, while the activation behavior of the cellular neuronal network of Caenorhabditis elegans is strongly influenced by hub nodes. These findings indicate the interaction of multiple topological features and dynamic states in the function of complex biological networks.

Dynamical consequences of lesions in cortical networks

To understand the effects of a cortical lesion it is necessary to consider not only the loss of local neural function, but also the lesion-induced changes in the larger network of endogenous oscillatory interactions in the brain. To investigate how network embedding influences a region's functional role, and the consequences of its being damaged, we implement two models of oscillatory cortical interactions, both of which inherit their coupling architecture from the available anatomical connection data for macaque cerebral cortex. In the first model, node dynamics are governed by Kuramoto phase oscillator equations, and we investigate the sequence in which areas entrain one another in the transition to global synchrony. In the second model, node dynamics are governed by a more realistic neural mass model, and we assess long-run inter-regional interactions using a measure of directed information flow. Highly connected parietal and frontal areas are found to synchronize most rapidly, more so than equally highly connected visual and somatosensory areas, and this difference can be explained in terms of the network's clustered architecture. For both models, lesion effects extend beyond the immediate neighbors of the lesioned site, and the amplitude and dispersal of nonlocal effects are again influenced by cluster patterns in the network. Although the consequences of in vivo lesions will always depend on circuitry local to the damaged site, we conclude that lesions of parietal regions (especially areas 5 and 7a) and frontal regions (especially areas 46 and FEF) have the greatest potential to disrupt the integrative aspects of neocortical function.

On the relationship between emotion and cognition

The current view of brain organization supports the notion that there is a considerable degree of functional specialization and that many regions can be conceptualized as either 'affective' or 'cognitive'. Popular examples are the amygdala in the domain of emotion and the lateral prefrontal cortex in the case of cognition. This prevalent view is problematic for a number of reasons. Here, I will argue that complex cognitive-emotional behaviours have their basis in dynamic coalitions of networks of brain areas, none of which should be conceptualized as specifically affective or cognitive. Central to cognitive-emotional interactions are brain areas with a high degree of connectivity, called hubs, which are critical for regulating the flow and integration of information between regions.

Motor imagery: a window into the mechanisms and alterations of the motor system

Motor imagery is a widely used paradigm for the study of cognitive aspects of action control, both in the healthy and the pathological brain. In this paper we review how motor imagery research has advanced our knowledge of behavioral and neural aspects of action control, both in healthy subjects and clinical populations. Furthermore, we will illustrate how motor imagery can provide new insights in a poorly understood psychopathological condition: conversion paralysis (CP). We measured behavioral and cerebral responses with functional magnetic resonance imaging (fMRI) in seven CP patients with a lateralized paresis of the arm as they imagined moving the affected or the unaffected hand. Imagined actions were either implicitly induced by the task requirements, or explicitly instructed through verbal instructions. We previously showed that implicitly induced motor imagery of the affected limb leads to larger ventromedial prefrontal responses compared to motor imagery of the unaffected limb. We interpreted this effect in terms of greater self-monitoring of actions during motor imagery of the affected limb. Here, we report new data in support of this interpretation: inducing self-monitoring of actions of both the affected and the unaffected limb (by means of explicitly cued motor imagery) abolishes the activation difference between the affected and the unaffected hand in the ventromedial prefrontal cortex. Our results show that although implicit and explicit motor imagery both entail motor simulations, they differ in terms of the amount of action monitoring they induce. The increased self-monitoring evoked by explicit motor imagery can have profound cerebral consequences in a psychopathological condition.

Mechanisms of hemispheric specialization: insights from analyses of connectivity

Traditionally, anatomical and physiological descriptions of hemispheric specialization have focused on hemispheric asymmetries of local brain structure or local functional properties, respectively. This article reviews the current state of an alternative approach that aims at unraveling the causes and functional principles of hemispheric specialization in terms of asymmetries in connectivity. Starting with an overview of the historical origins of the concept of lateralization, we briefly review recent evidence from anatomical and developmental studies that asymmetries in structural connectivity may be a critical factor shaping hemispheric specialization. These differences in anatomical connectivity, which are found both at the intra- and inter-regional level, are likely to form the structural substrate of different functional principles of information processing in the two hemispheres. The main goal of this article is to describe how these functional principles can be characterized using functional neuroimaging in combination with models of functional and effective connectivity. We discuss the methodology of established models of connectivity which are applicable to data from positron emission tomography and functional magnetic resonance imaging and review published studies that have applied these approaches to characterize asymmetries of connectivity during lateralized tasks. Adopting a model-based approach enables functional imaging to proceed from mere descriptions of asymmetric activation patterns to mechanistic accounts of how these asymmetries are caused.

The neural bases of prosopagnosia and pure alexia: recent insights from functional neuroimaging

PURPOSE OF REVIEW

To discuss whether recent functional neuroimaging results can account for clinical phenomenology in visual associative agnosias.

RECENT FINDINGS

Functional neuroimaging studies in healthy human subjects have identified only two regions of ventral occipitotemporal cortex that invariantly respond to individual faces and visual words, respectively. The signature of face identity coding in the fusiform neural response was shown to be missing in a patient with prosopagnosia. Another case study established that a surgical lesion close to the region sensitive to visual words can result in pure alexia.

SUMMARY

Evidence is increasing that functional specialization for processing face identity and visual word forms is restricted to two specialized sensory modules in the occipitotemporal cortex. A structural or functional lesion to face-sensitive and word-sensitive regions in the ventral occipitotemporal cortex can provide the most parsimonious account for the clinical syndromes of prosopagnosia and agnosic alexia. This review suggests that functional specialization should be considered in terms of whether exclusively one brain region (instead of many) underpins a defined function and not as whether this brain region underpins exclusively one cognitive function. Such functional specialization seems to exist for at least two higher-order visual perceptual functions, face and word identification.

Operational principles of neurocognitive networks

Large-scale neural networks are thought to be an essential substrate for the implementation of cognitive function by the brain. If so, then a thorough understanding of cognition is not possible without knowledge of how the large-scale neural networks of cognition (neurocognitive networks) operate. Of necessity, such understanding requires insight into structural, functional, and dynamical aspects of network operation, the intimate interweaving of which may be responsible for the intricacies of cognition. Knowledge of anatomical structure is basic to understanding how neurocognitive networks operate. Phylogenetically and ontogenetically determined patterns of synaptic connectivity form a structural network of brain areas, allowing communication between widely distributed collections of areas. The function of neurocognitive networks depends on selective activation of anatomically linked cortical and subcortical areas in a wide variety of configurations. Large-scale functional networks provide the cooperative processing which gives expression to cognitive function. The dynamics of neurocognitive network function relates to the evolving patterns of interacting brain areas that express cognitive function in real time. This article considers the proposition that a basic similarity of the structural, functional, and dynamical features of all neurocognitive networks in the brain causes them to function according to common operational principles. The formation of neural context through the coordinated mutual constraint of multiple interacting cortical areas, is considered as a guiding principle underlying all cognitive functions. Increasing knowledge of the operational principles of neurocognitive networks is likely to promote the advancement of cognitive theories, and to seed strategies for the enhancement of cognitive abilities.

The Cingulate as a Catalyst Region for Global Dysfunction: a Dynamical Modelling Paradigm

The anterior cingulate (AC) often exhibits both structural and functional abnormalities in affective disorders. Neither the cause for this vulnerability nor its effect on behaviour is known. Due to its extensive connectivity, minor output changes from the AC may exert widespread consequences. A causal model describing coupling coefficients (effective connectivity) among several brain regions in healthy subjects performing a memory task inspired our work. This stationary causal analysis provides a theoretical framework for our nonlinear dynamical models. We tested the effects of global and local perturbations upon stability of a systems-level neural network of interconnected brain regions. Interactions between regions, represented by path coefficients, were modelled using connectivity matrices. We found that both characteristic behaviour and response to perturbation differed in networks representing perceptual matching and long-delay conditions. Owing to the highly interconnected character of the networks, activation of a few areas was sufficient to trigger characteristic patterns of behaviour. However, only perturbation of key regions resulted in global dysfunction. Likewise, recovery of function was possible by increasing output from some, but not all, regions. We suggest for this recovery to be context specific, conditional on the task, integrity of other regions and global properties such as neuronal excitability.

On the role of general system theory for functional neuroimaging

One of the most important goals of neuroscience is to establish precise structure-function relationships in the brain. Since the 19th century, a major scientific endeavour has been to associate structurally distinct cortical regions with specific cognitive functions. This was traditionally accomplished by correlating microstructurally defined areas with lesion sites found in patients with specific neuropsychological symptoms. Modern neuroimaging techniques with high spatial resolution have promised an alternative approach, enabling non-invasive measurements of regionally specific changes of brain activity that are correlated with certain components of a cognitive process. Reviewing classic approaches towards brain structure-function relationships that are based on correlational approaches, this article argues that these approaches are not sufficient to provide an understanding of the operational principles of a dynamic system such as the brain but must be complemented by models based on general system theory. These models reflect the connectional structure of the system under investigation and emphasize context-dependent couplings between the system elements in terms of effective connectivity. The usefulness of system models whose parameters are fitted to measured functional imaging data for testing hypotheses about structure-function relationships in the brain and their potential for clinical applications is demonstrated by several empirical examples.

Network participation indices: characterizing component roles for information processing in neural networks

We propose a set of indices that characterize-on the basis of connectivity data-how a network node participates in a larger network and what roles it may take given the specific sub-network of interest. These Network Participation Indices are derived from simple graph theoretic measures and have the interesting property of linking local features of individual network components to distributed properties that arise within the network as a whole. We use connectivity data on large-scale cortical networks to demonstrate the virtues of this approach and highlight some interesting features that had not been brought up in previously published material. Some implications of our approach for defining network characteristics relevant to functional segregation and functional integration, for example, from functional imaging studies are discussed.

Cerebral localization, then and now

We review some of the progress made in understanding the nature of functional specialization in the human brain, beginning with the anatomical claim that all mental faculties have their own distinct material substrate in different regions of the brain and the psychological claim that each mental faculty is characterized by the content domain with which it deals. This conceptual framework led behavioral neurologists to show how discrete brain lesions provoked different types of language, praxic, gnostic, spatial, and memory disorders. The simplest way of interpreting these anatomoclinical associations was to conjecture that the normal function (now impaired by brain damage) was localized within that lesioned region. It was also realized that cognitive impairments could arise from lesions that spared the functional centers themselves but disconnected them from other centers. Nonetheless, many neuroscientists remained skeptical of the entire paradigm. Accordingly, in the late 19th century functional localization began to be studied in the intact human brain by such techniques as measuring the temperature of different brain regions when different cognitive tasks were performed. During the 20th century these crude techniques gave way to positron emission tomography, functional magnetic resonance imaging, and magnetoencephalography. The relatively precise spatial and temporal resolution of modern methods now raises a crucial question: Do the functional localizations obtained by the anatomoclinical method converge with those implied by the functional neuroimaging of cognition in healthy volunteers? We then conclude with some recent suggestions that functional specialization is not such a fixed property of brain regions as previously supposed.

Controlled analysis of neurocontrollers with informational lesioning

How does one aim to understand neural information processing? One of the difficult first challenges is to identify the roles of the network's elements. To this end a functional contribution analysis (FCA) method has been developed and applied for studying the neurocontrollers of evolutionary autonomous agents (EAAs). The FCA processes data composed of multiple lesion experiments and the corresponding performance levels that the agent obtains under these lesions. It calculates the contribution values (CVs) of the network's elements such that the ability to predict the agent's performance under new, unseen lesions is maximized. Previous analysis has found a strong dependence of the CVs and the prediction error on the specific type of lesioning method used, i.e. on the way in which the activity of lesioned neurons is disrupted. We present a new, informational lesioning method (ILM), which views a lesion as a noisy channel and applies a controlled lesion to the network by varying the lesioning level from large to arbitrarily small magnitudes. Studying the ILM within the FCA framework, our main results are threefold: first, that lower lesioning levels permit more accurate FCA predictions; second, that the usage of minute ILM lesioning levels can uncover the long-term effects of elements on the network's functioning; and third, that as the lesioning level decreases, the CVs tend to approach limit values, reflecting the importance of these elements in the intact, normal-functioning neurocontroller.

Localization of function via lesion analysis

This article presents a general approach for employing lesion analysis to address the fundamental challenge of localizing functions in a neural system. We describe functional contribution analysis (FCA), which assigns contribution values to the elements of the network such that the ability to predict the network's performance in response to multilesions is maximized. The approach is thoroughly examined on neurocontroller networks of evolved autonomous agents. The FCA portrays a stable set of neuronal contributions and accurate multilesion predictions that are significantly better than those obtained based on the classical single lesion approach. It is also used for a detailed synaptic analysis of the neurocontroller connectivity network, delineating its main functional backbone. The FCA provides a quantitative way of measuring how the network functions are localized and distributed among its elements. Our results question the adequacy of the classical single lesion analysis traditionally used in neuroscience and show that using lesioning experiments to decipher even simple neuronal systems requires a more rigorous multilesion analysis.

Evolutionary coherence of the mammalian amygdala

Despite great interest in the role of the amygdala in animal and human behaviour, its very existence as a structurally and functionally unified brain component has been questioned, on the grounds that cell groups within it display divergent pharmacological and connectional characteristics. We argue that the question of whether particular brain nuclei constitute a valid structural and functional unit is inherently an evolutionary question, and we present a method for answering it. The method involves phylogenetic analysis of comparative data to determine whether or not separate regions of the putative brain structure show statistically correlated evolution. We find that, in three separate groups of mammals (primates and two groups of insectivores), evolutionary changes in the volumes of amygdala components are strongly correlated, even after controlling for volumetric change in a wide range of limbic and other brain structures. This allows us to reject the strong claim that the amygdala is neither a structural nor a functional unit, and demonstrates the importance of evolutionary analysis in resolving such issues in systems neuroscience.