Causality and synchronisation in complex systems with applications to neuroscience

This thesis presents an investigation, of synchronisation and causality, motivated by problems in computational neuroscience. The thesis addresses both theoretical and practical signal processing issues regarding the estimation of interdependence from a set of multivariate data generated by a complex underlying dynamical system. This topic is driven by a series of problems in neuroscience, which represents the principal background motive behind the material in this work. The underlying system is the human brain and the generative process of the data is based on modern electromagnetic neuroimaging methods . In this thesis, the underlying functional of the brain mechanisms are derived from the recent mathematical formalism of dynamical systems in complex networks. This is justified principally on the grounds of the complex hierarchical and multiscale nature of the brain and it offers new methods of analysis to model its emergent phenomena. A fundamental approach to study the neural activity is to investigate the connectivity pattern developed by the brain’s complex network. Three types of connectivity are important to study: 1) anatomical connectivity refering to the physical links forming the topology of the brain network; 2) effective connectivity concerning with the way the neural elements communicate with each other using the brain’s anatomical structure, through phenomena of synchronisation and information transfer; 3) functional connectivity, presenting an epistemic concept which alludes to the interdependence between data measured from the brain network. The main contribution of this thesis is to present, apply and discuss novel algorithms of functional connectivities, which are designed to extract different specific aspects of interaction between the underlying generators of the data. Firstly, a univariate statistic is developed to allow for indirect assessment of synchronisation in the local network from a single time series. This approach is useful in inferring the coupling as in a local cortical area as observed by a single measurement electrode. Secondly, different existing methods of phase synchronisation are considered from the perspective of experimental data analysis and inference of coupling from observed data. These methods are designed to address the estimation of medium to long range connectivity and their differences are particularly relevant in the context of volume conduction, that is known to produce spurious detections of connectivity. Finally, an asymmetric temporal metric is introduced in order to detect the direction of the coupling between different regions of the brain. The method developed in this thesis is based on a machine learning extensions of the well known concept of Granger causality. The thesis discussion is developed alongside examples of synthetic and experimental real data. The synthetic data are simulations of complex dynamical systems with the intention to mimic the behaviour of simple cortical neural assemblies. They are helpful to test the techniques developed in this thesis. The real datasets are provided to illustrate the problem of brain connectivity in the case of important neurological disorders such as Epilepsy and Parkinson’s disease. The methods of functional connectivity in this thesis are applied to intracranial EEG recordings in order to extract features, which characterize underlying spatiotemporal dynamics before during and after an epileptic seizure and predict seizure location and onset prior to conventional electrographic signs. The methodology is also applied to a MEG dataset containing healthy, Parkinson’s and dementia subjects with the scope of distinguishing patterns of pathological from physiological connectivity.

Disconnected brains:what is the role of fMRI in connectivity research?

In this paper we consider how functional Magnetic Resonance Imaging (fMRI) has been used to study cortical connectivity in autism and autistic spectrum disorders (ASD). We discuss those studies that have contributed to the evidence supporting a model of disordered cortical connectivity in autism and (ASD), with a focusing emphasis on the application to research into the underconnectivity model. We note that the analytical techniques employed are limited and do not allow interpretation in terms of effective, or directional connectivity, nor do they provide information about the temporal or spectral characteristics of the functional networks being studied. We highlight how currently the features of neural generators that are being assessed by functional connectivity in fMRI are unclear. In addition, we note the importance in clinical studies of considering the consequences of task choice for the nature of the imaging data that can be collected and also of individual differences in participant state and trait characteristics for the accurate interpretation of imaging data. We discuss how alternative techniques such as EEG/MEG may address the limitations of fMRI in assessing brain connectivity, and additionally consider the potential of multimodal approaches. We conclude that fMRI has made significant contributions towards our understanding of the brain in terms of neural systems but that the conclusions drawn from its application in the sphere of autism research need to be approached with caution. It is important in research of this kind that we are aware of the need to examine the methodological and analytical techniques closely when applying findings in clinical populations, not only when they are used to support the development of theoretical models but also to inform diagnostic or treatment decisions.

Conectividade estrutural do cérebro

Mestrado em Radiações Aplicadas às Tecnologias da Saúde. Área de especialização: Ressonância Magnética

Arquitetura distribuída para análises multimodais de conectividade cerebral

Dissertação de mestrado integrado em Engenharia Biomédica (área de especialização em Informática Médica)

A Monte Carlo-Based Fiber Tracking Algorithm using Diffusion Tensor MRI

Diffusion tensor magnetic resonance imaging, which measures directional information of water diffusion in the brain, has emerged as a powerful tool for human brain studies. In this paper, we introduce a new Monte Carlo-based fiber tracking approach to estimate brain connectivity. One of the main characteristics of this approach is that all parameters of the algorithm are automatically determined at each point using the entropy of the eigenvalues of the diffusion tensor. Experimental results show the good performance of the proposed approach

The connectome mapper: an open-source processing pipeline to map connectomes with MRI.

Researchers working in the field of global connectivity analysis using diffusion magnetic resonance imaging (MRI) can count on a wide selection of software packages for processing their data, with methods ranging from the reconstruction of the local intra-voxel axonal structure to the estimation of the trajectories of the underlying fibre tracts. However, each package is generally task-specific and uses its own conventions and file formats. In this article we present the Connectome Mapper, a software pipeline aimed at helping researchers through the tedious process of organising, processing and analysing diffusion MRI data to perform global brain connectivity analyses. Our pipeline is written in Python and is freely available as open-source at www.cmtk.org.

Glutathione Precursor, N-Acetyl-Cysteine, Modulates EEG Synchronization in Schizophrenia Patients

Background: Glutathione (GSH) dysregulation at the gene, protein and functional levels observed in schizophrenia patients, and schizophrenia-like anomalies in GSH deficit experimental models, suggest that genetic glutathione synthesis impairments represent one major risk factor for the disease (Do et al., 2009). In a randomized, double blind, placebo controlled, add-on clinical trial of 140 patients, the GSH precursor N-Acetyl-Cysteine (NAC, 2g/day, 6 months) significantly improved the negative symptoms and reduced sideeffects due to antipsychotics (Berk et al., 2008). In a subset of patients (n=7), NAC (2g/day, 2 months, cross-over design) also improved auditory evoked potentials, the NMDA-dependent mismatch negativity (Lavoie et al, 2008). Methods: To determine whether increased GSH levels would modulate the topography of functional brain connectivity, we applied a multivariate phase synchronization (MPS) estimator (Knyazeva et al, 2008) to dense-array EEGs recorded during rest with eyes closed at the protocol onset, the point of crossover, and at its end. Results: The whole-head imaging revealed a specific synchronization landscape in NAC compared to placebo condition. In particular, NAC increased MPS over frontal and left temporal regions in a frequency-specific manner. The topography and direction of MPS changes were similar and robust in all 7 patients. Moreover, these changes correlated with the changes in the Liddle's score of disorganization, thus linking EEG synchronization to the improvement of the clinical picture. Conclusions: The data suggest an important pathway towards new therapeutic strategies that target GSH dysregulation in schizophrenia. They also show the utility of MPS mapping as a marker of treatment efficacy.

RECHERCHE TRANSLATIONNELLE EN NEUROSCIENCES PSYCHIATRIQUES: DE L'ANIMAL MODELE DE SCHIZOPHRENIE A LA DETECTION D'UNE DESORIENTATION SPATIALE, PRATIQUE OU SYMBOLIQUE CHEZ DES PATIENTS

La modélisation, chez l'animal, de maladies psychiatriques telles que la schizophrénie repose sur différentes démarches visant à induire des perturbations cérébrales similaires à celles observées dans la maladie. Nous avons cherché à étudier chez le rat les effets d'une diminution (50%) transitoire en glutathion (GSH) durant le développement (PND 5 à PND 16) à partir de l'implication, chez des adultes, des conséquences de cette perturbation dans des mécanismes fondamentaux de traitement de l'information sensorielle. Cette thèse évalue et documente les déficits de compétences de navigation spatiale dans ce modèle. Nous avons mis en évidence des effets comportementaux à partir de l'identification de différences particulières dans des tâches d'orientation: des difficultés, chez les rats ayant subi un déficit en GSH, à élaborer une représentation globale de l'environnement dans lequel ils se déplacent, difficultés compensées par une attention particulière aux détails visuels le composant. Cette stratégie réactive compensatoire est efficace lorsque les conditions permettent un ajustement continu aux repères visuels environnementaux. Elle ne permet cependant pas des prédictions et des attentes sur ce qui devrait être rencontré et perçu dans une certaine direction, dès qu'une partie des informations visuelles familières disparaît. Il faudrait pour cela une capacité fondée sur une représentation abstraite, à distance des modalités sensorielles qui en ont permis son élaboration. Notre thèse soutient que les déficits, supposés participer à l'émergence de certains symptômes de la maladie, auraient également des conséquences sur l'élaboration de la représentation spatiale nécessaire à des capacités d'orientation effectives et symboliques. - The study of a psychiatric disease such as schizophrenia in an animal model relies on different approaches attempting to replicate brain perturbations similar to those observed in the illness. In the present work, behavioural consequences of a functional deficit in brain connectivity and coordination were assessed in rats with a transitory glutathione (GSH) deficit induced during the postnatal development (PND 5-PND 16) with daily injections of BSO (1- buthionine-(S,R)- sulfoximine). We searched for a theoretical syndrome associating ecologically relevant behavioural adaptive deficits and resulting from the weakening of sensory integration processes. Our results revealed significant and specific deficit of BSO treated rats in spatial orientation tasks designed to test for cognitive mapping abilities. Treated rats behaved as if impaired in the proactive strategies supported by an abstract representation such as a cognitive map. In contrast their performances were preserved whenever the environmental conditions allowed for adaptative reactive strategies, an equivalent of the visual affordances described by Gibson (1958). This supports our thesis that BSO treated rats expressed difficulties in elaborating a global representation of the environment. This deficit was completely - or - partially compensated by the development of an increased attention to the environment's visual details. This compensatory reactive strategy requires a rich environment allowing for continuous adjustment to visual cues. However, such adjustment doesn't allow to predictions and expectancies about what should be met and perceived in a certain direction, when familiar visual spatial cues are missing. Such competencies require orientation based on the use of an abstract spatial representation, independent from the specific sensory modalities that have participated to its elaboration. The impairment of BSO rats such spatial representation could result from a deficit in the integration and organization of perceptual information. Our model leads to the hypothesis that these fundamental deficits might account for certain symptoms of schizophrenia. They would also interfere with in the capacity to elaborate spatial representation necessary for optimal orientation in natural, artificial or symbolic environment.

A la recherche de symptômes de désorientation dépendant d'un déficit d'intégration multisensorielle chez l'animal : un éclairage sur les mécanismes fondamentaux de la schizophrénie

Abstract Fundamental research in psychiatric neurosciences assumes that psychiatric disorders are associated with neurobiological factors. Identification of these factors would provide therapeutic targets as well as a better understanding of the relationship between- brain and behaviour in pathological processes. We conducted experiments in an animal model of schizophrenia. Several behavioural tasks were used to evaluate spatial and working memory in these animals. The model is based on glutathione deficit during cerebral development. Indeed, a 50% decrease of glutathione has been reported in prefrontal cortex of patients with schizophrenia. Glutathione is a major antioxidant in the brain and its deficit could lead to abnormal brain connectivity. The glutathione deficit was induced in rats by perinatal (PS-P16) subcutaneous injections with Lbuthionine-(S,R)-sullfoximine (BSO), an inhibitor of glutathione synthesis. This treatment leads to a transitory 50% glutathione levels during brain development. In parallel, we conducted behavioural testing in rats with a medial prefrontal cortex lesion. This allowed us to compare early damage induced by BSO treatment with a focal lesion in adults of a brain area known to present anomalies in schizophrenia. Finally, we conducted a series of experiments in senescent rats to evaluate if cognitive deficits could be related to neurobiological changes. Our results show that an early glutathione deficit provokes cognitive deficits in adulthood. These spatial and working memory deficits resemble the cognitive deficits observed in schizophrenia. The comparison with prefrontal rats revealed that the early brain glutathione deficit provoked more severe cognitive deficits than the prefrontal lesion in adult rats. Moreover, in both cases, we observed a dissociation in memory deficits depending on the type of locomotion that was used in behavioural experiments. Indeed, BSO treated rats as well as prefrontal rats showed place learning or working memory deficits in tasks conducted on dry surfaces where they had to walk. In contrast, they showed no deficit when the same cognitive functions were tested in the water maze. This dissociation might be sustained by a difference in requirement of sensory integration between walking and swimming tasks. Résumé La recherche fondamentale en neurosciences psychiatriques repose sur le présupposé selon lequel les symptômes manifestés dans les troubles psychiatriques auraient des concomitants neurobiologiques. Ceux-ci, une fois identifiés, offriraient des cibles pour une démarche thérapeutique ainsi que des modèles permettant de mieux comprendre les soubassements biologiques du comportement et des activités mentales. Nos expériences s'articulent autour de la question de la modélisation de la schizophrénie chez l'animal. Nous avons recherché chez ces animaux des troubles cognitifs et sensoriels associés à la schizophrénie. En effet, chez l'homme comme chez l'animal, la mémoire spatiale et la mémoire de travail dépendent fortement de la capacité d'intégration et d'organisation des informations sensorielles. Les premières expériences ont été menées suite à une perturbation périnatale du développement cérébral. Celle-ci visait à reproduire une diminution du taux de glutathion dans le cerveau, des recherches précédentes ayant observé une diminution de 50% du taux de glutathion dans le cortex préfrontal de patients schizophrènes. Le glutathion étant un antioxydant majeur dans le cerveau, son déficit pourrait conduire à des perturbations de la circuiterie cérébrale. Nous avons reproduit ce déficit chez le rat, par injection de Lbuthionine-(S,R)-sullfoximine (BSO), un inhibiteur de la synthèse du glutathion... Ce traitement a été administré pendant la période périnatale (du jour postnatal 5 au jour 16) provoquant une diminution de 50% du taux de glutathion. Nous avons ensuite évalué lës répercussions de cette atteinte précoce sur le comportement des rats à l'âge adulte. Ce modèle s'inscrit donc dans l'hypothèse neurodéveloppementale qui associe la schizophrénie à une atteinte du développement cérébral normal. Nous avons ensuite conduit des expériences similaires chez des rats ayant subi une lésion du cortex préfrontal pour comparer les répercussions du traitement périnatal avec une lésion, à l'âge adulte, d'une aire cérébrale connue pour présenter des anomalies chez les patients. Finalement, nous avons évalué si les processus sensoriels et cognitifs précédemment étudiés pouvaient également être affectés lors du vieillissement normal en recherchant des corrélats biologiques des déficits de mémoire liés à l'âge avancé. Nos résultats montrent que ce déficit précoce en glutathion peut avoir des répercussions surale comportement à l'âge adulte. On a relevé une similarité avec les déficits cognitifs associés.à la schizophrénie, incluant des déficits de mémoire de travail ainsi que des déficits de mémoire spatiale. Ces déficits étaient fortement liés au type de locomotion utilisée et n'ont été observés que dans les tâches où les animaux devaient rejoindre un but en marchant mais pas dans lés tests dans lesquels ils devaient localiser une cible en nageant. Les déficits induits par la lésion préfrontale chez l'adulte étaient beaucoup plus légers que ceux découlant de l'atteinte périnatale mais présentaient une dissociation analogue en fonction du type de locomotion. De plus, des tests similaires menés au cours du vieillissement confirment que la mémoire de travail peut être affectée sélectivement par le vieillissement dans une tâche où les animaux doivent marcher, tout en restant intacte dans le bassin de Morris. Les déficits cognitifs liés au vieillissement étaient significativement corrélés à des différences de niveaux des protéines post-synaptiques PSD95 (postsynaptic density 95). L'ensemble des résultats montre que les tests qui sont fréquemment utilisés pour évaluer la mémoire chez l'animal pourraient faire appel à des processus différents. Cette différence pourrait notamment tenir au niveau d'intégration sensorielle requis pour résoudre la tâche, qui est particulièrement sollicitée au cours d'une locomotion intermittente.

Seizure semiology: an overview of the 'inverse problem'.

In clinical practice, a classification of seizures based on clinical signs and symptoms leads to an improved understanding of epilepsy-related issues and therefore strongly contributes to a better patient care. The inverse problem involves inferring the anatomical brain localization of a seizure from the scalp surface EEG, a concept we apply here to correlate seizure origin with seizure semiology. The spheres of sensorium, motor features, consciousness changes and autonomic alterations during ictal and postictal manifestations are reviewed, including several subdivisions used to better categorize particular features. Particular attention is given to behavioral features, as well as to features occurring in idiopathic generalized epileptic syndromes and psychogenic nonepileptic spells.

Comparing connectomes across subjects and populations at different scales.

Brain connectivity can be represented by a network that enables the comparison of the different patterns of structural and functional connectivity among individuals. In the literature, two levels of statistical analysis have been considered in comparing brain connectivity across groups and subjects: 1) the global comparison where a single measure that summarizes the information of each brain is used in a statistical test; 2) the local analysis where a single test is performed either for each node/connection which implies a multiplicity correction, or for each group of nodes/connections where each subset is summarized by one single test in order to reduce the number of tests to avoid a penalizing multiplicity correction. We comment on the different levels of analysis and present some methods that have been proposed at each scale. We highlight as well the possible factors that could influence the statistical results and the questions that have to be addressed in such an analysis.

Increased eeg synchronization in schizophrenia patients by the glutathione precursor, n-acetyl-cysteine

Background: Glutathione (GSH) dysregulation at the gene, protein and functional levels observed in schizophrenia patients, and schizophrenia-like anomalies in GSH deficit experimental models, suggest that genetic glutathione synthesis impairments represent one major risk factor for the disease (Do et al., 2009). In a randomized, double blind, placebo controlled, add-on clinical trial of 140 patients, the GSH precursor N-Acetyl-Cysteine (NAC, 2 g/day, 6 months) significantly improved the negative symptoms and reduced side-effects due to antipsychotics (Berk et al., 2008). In a subset of patients (n=7), NAC (2 g/day, 2 months, cross-over design) also improved auditory evoked potentials, the NMDAdependent mismatch negativity (Lavoie et al, 2008). Methods: To determine whether increased GSH levels would modulate the topography of functional brain connectivity, we applied a multivariate phase synchronization (MPS) estimator (Knyazeva et al, 2008) to dense-array EEGs recorded during rest with eyes closed at the protocol onset, the point of crossover, and at its end. Phase synchronization phenomena are appealing because they can be associated to synchronized phases while the amplitudes stay uncorrelated. MPS measures the degree of interactions among the recorded neuronal oscillators by quantifiying to what extent they behave like a macro-oscillator (i.e. the oscillators are phase synchronous). To assess the whole-head synchronization topography, we computed the MPS sensor-wise over the cluster of locations defined by the sensor itself and he surrounding ones belonging to its second-order neighborhood (Carmeli et al, 2005). Such a cluster spans about 12 cm on average. Abstracts 245 Results: The whole-head imaging revealed a specific synchronization landscape in NAC compared to placebo condition. In particular, NAC increased MPS over frontal and left temporal regions in a frequency-specific manner. Importantly, the topography and direction of MPS changes were similar and robust in all 7 patients. Moreover, these changes correlated with the changes in the Liddle's score of disorganization (Liddle, 1987) thus linking EEG synchronization to the improvement of clinical picture. Discussion: The data suggest an important pathway towards new therapeutic strategies that target GSH dysregulation in schizophrenia. They also show the utility of MPS mapping as a marker of treatment efficacy.

ConnectomeViewer - Multi-Modal Multi-Level Network Visualization and Analysis.

Introduction: The field of Connectomic research is growing rapidly, resulting from methodological advances in structural neuroimaging on many spatial scales. Especially progress in Diffusion MRI data acquisition and processing made available macroscopic structural connectivity maps in vivo through Connectome Mapping Pipelines (Hagmann et al, 2008) into so-called Connectomes (Hagmann 2005, Sporns et al, 2005). They exhibit both spatial and topological information that constrain functional imaging studies and are relevant in their interpretation. The need for a special-purpose software tool for both clinical researchers and neuroscientists to support investigations of such connectome data has grown. Methods: We developed the ConnectomeViewer, a powerful, extensible software tool for visualization and analysis in connectomic research. It uses the novel defined container-like Connectome File Format, specifying networks (GraphML), surfaces (Gifti), volumes (Nifti), track data (TrackVis) and metadata. Usage of Python as programming language allows it to by cross-platform and have access to a multitude of scientific libraries. Results: Using a flexible plugin architecture, it is possible to enhance functionality for specific purposes easily. Following features are already implemented: * Ready usage of libraries, e.g. for complex network analysis (NetworkX) and data plotting (Matplotlib). More brain connectivity measures will be implemented in a future release (Rubinov et al, 2009). * 3D View of networks with node positioning based on corresponding ROI surface patch. Other layouts possible. * Picking functionality to select nodes, select edges, get more node information (ConnectomeWiki), toggle surface representations * Interactive thresholding and modality selection of edge properties using filters * Arbitrary metadata can be stored for networks, thereby allowing e.g. group-based analysis or meta-analysis. * Python Shell for scripting. Application data is exposed and can be modified or used for further post-processing. * Visualization pipelines using filters and modules can be composed with Mayavi (Ramachandran et al, 2008). * Interface to TrackVis to visualize track data. Selected nodes are converted to ROIs for fiber filtering The Connectome Mapping Pipeline (Hagmann et al, 2008) processed 20 healthy subjects into an average Connectome dataset. The Figures show the ConnectomeViewer user interface using this dataset. Connections are shown that occur in all 20 subjects. The dataset is freely available from the homepage (connectomeviewer.org). Conclusions: The ConnectomeViewer is a cross-platform, open-source software tool that provides extensive visualization and analysis capabilities for connectomic research. It has a modular architecture, integrates relevant datatypes and is completely scriptable. Visit www.connectomics.org to get involved as user or developer.

Mapping the human connectome at multiple scales with diffusion spectrum MRI.

The global structural connectivity of the brain, the human connectome, is now accessible at millimeter scale with the use of MRI. In this paper, we describe an approach to map the connectome by constructing normalized whole-brain structural connection matrices derived from diffusion MRI tractography at 5 different scales. Using a template-based approach to match cortical landmarks of different subjects, we propose a robust method that allows (a) the selection of identical cortical regions of interest of desired size and location in different subjects with identification of the associated fiber tracts (b) straightforward construction and interpretation of anatomically organized whole-brain connection matrices and (c) statistical inter-subject comparison of brain connectivity at various scales. The fully automated post-processing steps necessary to build such matrices are detailed in this paper. Extensive validation tests are performed to assess the reproducibility of the method in a group of 5 healthy subjects and its reliability is as well considerably discussed in a group of 20 healthy subjects.