Inferring neuronal network connectivity from spike data: A temporal data mining approach

Understanding the functioning of a neural system in terms of its underlying circuitry is an important problem in neuroscience. Recent d evelopments in electrophysiology and imaging allow one to simultaneously record activities of hundreds of neurons. Inferring the underlying neuronal connectivity patterns from such multi-neuronal spike train data streams is a challenging statistical and computational problem. This task involves finding significant temporal patterns from vast amounts of symbolic time series data. In this paper we show that the frequent episode mining methods from the field of temporal data mining can be very useful in this context. In the frequent episode discovery framework, the data is viewed as a sequence of events, each of which is characterized by an event type and its time of occurrence and episodes are certain types of temporal patterns in such data. Here we show that, using the set of discovered frequent episodes from multi-neuronal data, one can infer different types of connectivity patterns in the neural system that generated it. For this purpose, we introduce the notion of mining for frequent episodes under certain temporal constraints; the structure of these temporal constraints is motivated by the application. We present algorithms for discovering serial and parallel episodes under these temporal constraints. Through extensive simulation studies we demonstrate that these methods are useful for unearthing patterns of neuronal network connectivity.

Analysis of connectivity map: Control to glutamate injured and phenobarbital treated neuronal network

We study the responses of a cultured neural network when it is exposed to epileptogenesis glutamate injury causing epilepsy and subsequent treatment with phenobarbital by constructing connectivity map of neurons using correlation matrix. This study is particularly useful in understanding the pharmaceutical drug induced changes in the neuronal network properties with insights into changes at the systems biology level. (C) 2010 American Institute of Physics. [doi:10.1063/1.3398025]

Phase transitions and nonlinear phenomena in neuronal network models

Communication and cooperation between billions of neurons underlie the power of the brain. How do complex functions of the brain arise from its cellular constituents? How do groups of neurons self-organize into patterns of activity? These are crucial questions in neuroscience. In order to answer them, it is necessary to have solid theoretical understanding of how single neurons communicate at the microscopic level, and how cooperative activity emerges. In this thesis we aim to understand how complex collective phenomena can arise in a simple model of neuronal networks. We use a model with balanced excitation and inhibition and complex network architecture, and we develop analytical and numerical methods for describing its neuronal dynamics. We study how interaction between neurons generates various collective phenomena, such as spontaneous appearance of network oscillations and seizures, and early warnings of these transitions in neuronal networks. Within our model, we show that phase transitions separate various dynamical regimes, and we investigate the corresponding bifurcations and critical phenomena. It permits us to suggest a qualitative explanation of the Berger effect, and to investigate phenomena such as avalanches, band-pass filter, and stochastic resonance. The role of modular structure in the detection of weak signals is also discussed. Moreover, we find nonlinear excitations that can describe paroxysmal spikes observed in electroencephalograms from epileptic brains. It allows us to propose a method to predict epileptic seizures. Memory and learning are key functions of the brain. There are evidences that these processes result from dynamical changes in the structure of the brain. At the microscopic level, synaptic connections are plastic and are modified according to the dynamics of neurons. Thus, we generalize our cortical model to take into account synaptic plasticity and we show that the repertoire of dynamical regimes becomes richer. In particular, we find mixed-mode oscillations and a chaotic regime in neuronal network dynamics.

Architecture for Living Neuronal Network Control of a mobile robot

The intelligent controlling mechanism of a typical mobile robot is usually a computer system. Research is however now ongoing in which biological neural networks are being cultured and trained to act as the brain of an interactive real world robot – thereby either completely replacing or operating in a cooperative fashion with a computer system. Studying such neural systems can give a distinct insight into biological neural structures and therefore such research has immediate medical implications. The principal aims of the present research are to assess the computational and learning capacity of dissociated cultured neuronal networks with a view to advancing network level processing of artificial neural networks. This will be approached by the creation of an artificial hybrid system (animat) involving closed loop control of a mobile robot by a dissociated culture of rat neurons. This paper details the components of the overall animat closed loop system architecture and reports on the evaluation of the results from preliminary real-life and simulated robot experiments.

Multivariate autoregressive-based neuronal network flow analysis for in-vitro recorded bursts

Neuroscientific studies of in vitro neuron cell cultures has attracted paramount attention to investigate the behaviour of neuronal networks in response to different environmental conditions and external stimuli such as drugs, optical and electrical stimulations. Microelec trodearray (MEA) technology has been widely adopted as a tool for this investigation. In this work, we present a new approach to estimate interconnectivity of neural spikes using multivariate autoregressive (MVAR) analysis and Partial Directed Coherence (PDC). The proposed approach has the potential to discover hidden intra-burst causal connectivity patterns and to help understand the spatiotemporal communication patterns within bursts, pre and post stimulations.

The influence of connectivity on the firing rate in a neuronal network with electrical and chemical synapses

We study the firing rate properties of a cellular automaton model for a neuronal network with chemical synapses. We propose a simple mechanism in which the nonlocal connections are included, through electrical and chemical synapses. In the latter case, we introduce a time delay which produces self-sustained activity. Nonlocal connections, or shortcuts, are randomly introduced according to a specified connection probability. There is a range of connection probabilities for which neuron firing occurs, as well as a critical probability for which the firing ceases in the absence of time delay. The critical probability for nonlocal shortcuts depends on the network size according to a power-law. We also compute the firing rate amplification factor by varying both the connection probability and the time delay for different network sizes. (C) 2011 Elsevier B.V. All rights reserved.

Tactile Perception - Perception of tactile distance on a single skin surface changes with stimulus orientation: a neuronal network modelling study

L’interazione che abbiamo con l’ambiente che ci circonda dipende sia da diverse tipologie di stimoli esterni che percepiamo (tattili, visivi, acustici, ecc.) sia dalla loro elaborazione per opera del nostro sistema nervoso. A volte però, l’integrazione e l’elaborazione di tali input possono causare effetti d’illusione. Ciò si presenta, ad esempio, nella percezione tattile. Infatti, la percezione di distanze tattili varia al variare della regione corporea considerata. Il concetto che distanze sulla cute siano frequentemente erroneamente percepite, è stato scoperto circa un secolo fa da Weber. In particolare, una determinata distanza fisica, è percepita maggiore su parti del corpo che presentano una più alta densità di meccanocettori rispetto a distanze applicate su parti del corpo con inferiore densità. Oltre a questa illusione, un importante fenomeno osservato in vivo è rappresentato dal fatto che la percezione della distanza tattile dipende dall’orientazione degli stimoli applicati sulla cute. In sostanza, la distanza percepita su una regione cutanea varia al variare dell’orientazione degli stimoli applicati. Recentemente, Longo e Haggard (Longo & Haggard, J.Exp.Psychol. Hum Percept Perform 37: 720-726, 2011), allo scopo di investigare come sia rappresentato il nostro corpo all’interno del nostro cervello, hanno messo a confronto distanze tattili a diverse orientazioni sulla mano deducendo che la distanza fra due stimoli puntuali è percepita maggiore se applicata trasversalmente sulla mano anziché longitudinalmente. Tale illusione è nota con il nome di Illusione Tattile Orientazione-Dipendente e diversi risultati riportati in letteratura dimostrano che tale illusione dipende dalla distanza che intercorre fra i due stimoli puntuali sulla cute. Infatti, Green riporta in un suo articolo (Green, Percpept Pshycophys 31, 315-323, 1982) il fatto che maggiore sia la distanza applicata e maggiore risulterà l’effetto illusivo che si presenta. L’illusione di Weber e l’illusione tattile orientazione-dipendente sono spiegate in letteratura considerando differenze riguardanti la densità di recettori, gli effetti di magnificazione corticale a livello della corteccia primaria somatosensoriale (regioni della corteccia somatosensoriale, di dimensioni differenti, sono adibite a diverse regioni corporee) e differenze nella dimensione e forma dei campi recettivi. Tuttavia tali effetti di illusione risultano molto meno rilevanti rispetto a quelli che ci si aspetta semplicemente considerando i meccanismi fisiologici, elencati in precedenza, che li causano. Ciò suggerisce che l’informazione tattile elaborata a livello della corteccia primaria somatosensoriale, riceva successivi step di elaborazione in aree corticali di più alto livello. Esse agiscono allo scopo di ridurre il divario fra distanza percepita trasversalmente e distanza percepita longitudinalmente, rendendole più simili tra loro. Tale processo assume il nome di “Rescaling Process”. I meccanismi neurali che operano nel cervello allo scopo di garantire Rescaling Process restano ancora largamente sconosciuti. Perciò, lo scopo del mio progetto di tesi è stato quello di realizzare un modello di rete neurale che simulasse gli aspetti riguardanti la percezione tattile, l’illusione orientazione-dipendente e il processo di rescaling avanzando possibili ipotesi circa i meccanismi neurali che concorrono alla loro realizzazione. Il modello computazionale si compone di due diversi layers neurali che processano l’informazione tattile. Uno di questi rappresenta un’area corticale di più basso livello (chiamata Area1) nella quale una prima e distorta rappresentazione tattile è realizzata. Per questo, tale layer potrebbe rappresentare un’area della corteccia primaria somatosensoriale, dove la rappresentazione della distanza tattile è significativamente distorta a causa dell’anisotropia dei campi recettivi e della magnificazione corticale. Il secondo layer (chiamato Area2) rappresenta un’area di più alto livello che riceve le informazioni tattili dal primo e ne riduce la loro distorsione mediante Rescaling Process. Questo layer potrebbe rappresentare aree corticali superiori (ad esempio la corteccia parietale o quella temporale) adibite anch’esse alla percezione di distanze tattili ed implicate nel Rescaling Process. Nel modello, i neuroni in Area1 ricevono informazioni dagli stimoli esterni (applicati sulla cute) inviando quindi informazioni ai neuroni in Area2 mediante sinapsi Feed-forward eccitatorie. Di fatto, neuroni appartenenti ad uno stesso layer comunicano fra loro attraverso sinapsi laterali aventi una forma a cappello Messicano. E’ importante affermare che la rete neurale implementata è principalmente un modello concettuale che non si preme di fornire un’accurata riproduzione delle strutture fisiologiche ed anatomiche. Per questo occorre considerare un livello astratto di implementazione senza specificare un’esatta corrispondenza tra layers nel modello e regioni anatomiche presenti nel cervello. Tuttavia, i meccanismi inclusi nel modello sono biologicamente plausibili. Dunque la rete neurale può essere utile per una migliore comprensione dei molteplici meccanismi agenti nel nostro cervello, allo scopo di elaborare diversi input tattili. Infatti, il modello è in grado di riprodurre diversi risultati riportati negli articoli di Green e Longo & Haggard.

A hierarchical neuronal network for planning behavior

Planning a goal-directed sequence of behavior is a higher function of the human brain that relies on the integrity of prefrontal cortical areas. In the Tower of London test, a puzzle in which beads sliding on pegs must be moved to match a designated goal configuration, patients with lesioned prefrontal cortex show deficits in planning a goal-directed sequence of moves. We propose a neuronal network model of sequence planning that passes this test and, when lesioned, fails in a way that mimics prefrontal patients’ behavior. Our model comprises a descending planning system with hierarchically organized plan, operation, and gesture levels, and an ascending evaluative system that analyzes the problem and computes internal reward signals that index the correct/erroneous status of the plan. Multiple parallel pathways connecting the evaluative and planning systems amend the plan and adapt it to the current problem. The model illustrates how specialized hierarchically organized neuronal assemblies may collectively emulate central executive or supervisory functions of the human brain.

Modulation of neuronal network activity in the primary motor cortex

In the present study I investigated the mechanisms of modulation of neuronal network activity in rat primary motor cortex using pharmacological manipulations employing the in vitro brain slice technique. Preparation of the brain slice in sucrose-based aCSF produced slices with low viability. Introducing the neuroprotectants N-acetyl-cysteine, taurine and aminoguanidine to the preparatory method saw viability of slices increase significantly. Co-application of low dose kainic acid and carbachol consistently generated beta oscillatory activity in M1. Analyses indicated that network activity in M1 relied on the involvement of GABAA receptors. Dose-response experiments performed in M1 showed that beta activity can be modulated by benzodiazepine site ligands. Low doses of positive allosteric modulators consistently desynchronised beta oscillatory activity, a mechanism that may be driven by a1-subunit containing GABAA receptors. Higher doses increased the power of beta oscillatory activity. Whole-cell recordings in M1 uncovered three interneuronal subtypes regularly encountered in M1; Fast-spiking, regular-spiking non-Pyramidal and low threshold spiking. With the paradoxical effects of positive allosteric modulators in mind, subsequent voltage-clamp recordings in FS cells revealed a constitutively active tonic inhibitory current that could be modulated by zolpidem in two different ways. Low dose zolpidem increased the tonic inhibitory current in FS cells, consistent with the desynchronisation of network oscillatory activity seen at this concentration. High dose zolpidem decreased the inhibitory tonic current seen in FS cells, coinciding with an increase in oscillatory power. These studies indicate a fundamental role for a tonic inhibitory current in the modulation of network activity. Furthermore, desynchronisation of beta activity in M1 decreased as viability of the in vitro brain slice increased, suggesting that the extent of desynchronisation is dependent upon the pathophysiological state of the network. This indicates that low dose zolpidem could be used as a therapeutic agent specifically for the desynchronisation of pathological oscillations in oscillopathies such as Parkinson’s disease.

Neuronal network pharmacodynamics of GABAergic modulation in the human cortex determined using pharmaco-magnetoencephalography

Neuronal network oscillations are a unifying phenomenon in neuroscience research, with comparable measurements across scales and species. Cortical oscillations are of central importance in the characterization of neuronal network function in health and disease and are influential in effective drug development. Whilst animal in vitro and in vivo electrophysiology is able to characterize pharmacologically induced modulations in neuronal activity, present human counterparts have spatial and temporal limitations. Consequently, the potential applications for a human equivalent are extensive. Here, we demonstrate a novel implementation of contemporary neuroimaging methods called pharmaco-magnetoencephalography. This approach determines the spatial profile of neuronal network oscillatory power change across the cortex following drug administration and reconstructs the time course of these modulations at focal regions of interest. As a proof of concept, we characterize the nonspecific GABAergic modulator diazepam, which has a broad range of therapeutic applications. We demonstrate that diazepam variously modulates ? (4–7 Hz), a (7–14 Hz), ß (15–25 Hz), and ? (30–80 Hz) frequency oscillations in specific regions of the cortex, with a pharmacodynamic profile consistent with that of drug uptake. We examine the relevance of these results with regard to the spatial and temporal observations from other modalities and the various therapeutic consequences of diazepam and discuss the potential applications of such an approach in terms of drug development and translational neuroscience.

The role of GABAergic modulation in motor function related neuronal network activity

At rest, the primary motor cortex (M1) exhibits spontaneous neuronal network oscillations in the beta (15–30 Hz) frequency range, mediated by inhibitory interneuron drive via GABA-A receptors. However, questions remain regarding the neuropharmacological basis of movement related oscillatory phenomena, such as movement related beta desynchronisation (MRBD), post-movement beta rebound (PMBR) and movement related gamma synchronisation (MRGS). To address this, we used magnetoencephalography (MEG) to study the movement related oscillatory changes in M1 cortex of eight healthy participants, following administration of the GABA-A modulator diazepam. Results demonstrate that, contrary to initial hypotheses, neither MRGS nor PMBR appear to be GABA-A dependent, whilst the MRBD is facilitated by increased GABAergic drive. These data demonstrate that while movement-related beta changes appear to be dependent upon spontaneous beta oscillations, they occur independently of one other. Crucially, MRBD is a GABA-A mediated process, offering a possible mechanism by which motor function may be modulated. However, in contrast, the transient increase in synchronous power observed in PMBR and MRGS appears to be generated by a non-GABA-A receptor mediated process; the elucidation of which may offer important insights into motor processes.

Beta frequency neuronal network activity in the primary motor cortex

In this study I investigated the mechanisms of neuronal network oscillatory activity in rat M1 using pharmacological manipulations and electrical stimulation protocols, employing the in vitro brain slice technique in rat and magnetoencephalography (MEG) in man. Co-application of kainic acid and carbachol generated in vitro beta oscillatory activity in all layers in M1. Analyses indicated that oscillations originated from deep layers and indicated significant involvement of GABAA receptors and gap junctions. A modulatory role of GABAB, NMDA, and dopamine receptors was also evident. Intracellular recordings from fast-spiking (FS) GABAergic inhibitory cells revealed phase-locked action potentials (APs) on every beta cycle. Glutamatergic excitatory regular-spiking (RS) and intrinsically-bursting (IB) cells both received phase locked inhibitory postsynaptic potentials, but did not fire APs on every cycle, suggesting the dynamic involvement of different pools of neurones in the overall population oscillations. Stimulation evoked activity at high frequency (HFS; 125Hz) evoked gamma oscillations and reduced ongoing beta activity. 20Hz stimulation promoted theta or gamma oscillations whilst 4Hz stimulation enhanced beta power at theta frequency. I also investigated the modulation of pathological slow wave (theta and beta) oscillatory activity using magnetoencephalography. Abnormal activity was suppressed by sub-sedative doses of GABAA receptor modulator zolpidem and the observed desynchronising effect correlated well with improved sensorimotor function. These studies indicate a fundamental role for inhibitory neuronal networks in the patterning beta activity and suggest that cortical HFS in PD re-patterns abnormally enhanced M1 network activity by modulating the activity of FS cells. Furthermore, pathological oscillation may be common to many neuropathologies and may be an important future therapeutic target.

Functional astrocyte-neuron lactate shuttle in a human stem cell-derived neuronal network

The NT2.D1 cell line is one of the most well-documented embryocarcinoma cell lines, and can be differentiated into neurons and astrocytes. Great focus has also been placed on defining the electrophysiological properties of the neuronal cells, and more recently we have investigated the functional properties of their associated astrocytes. We now show for the first time that human stem cell-derived astrocytes produce glycogen and that co-cultures of these cells demonstrate a functional astrocyte-neuron lactate shuttle (ANLS). The ANLS hypothesis proposes that during neuronal activity, glutamate released into the synaptic cleft is taken up by astrocytes and triggers glucose uptake, which is converted into lactate and released via monocarboxylate transporters for neuronal use. Using mixed cultures of NT2-derived neurons and astrocytes, we have shown that these cells modulate their glucose uptake in response to glutamate. Additionally, we demonstrate that in response to increased neuronal activity and under hypoglycaemic conditions, co-cultures modulate glycogen turnover and increase lactate production. Similar results were also shown after treatment with glutamate, potassium, isoproterenol, and dbcAMP. Together, these results demonstrate for the first time a functional ANLS in a human stem cell-derived co-culture. © 2013 ISCBFM.

Neuronal network dynamics during epileptogenesis in the medial temporal lobe

Epilepsy is one of the most common neurological disorders, a large fraction of which is resistant to pharmacotherapy. In this light, understanding the mechanisms of epilepsy and its intractable forms in particular could create new targets for pharmacotherapeutic intervention. The current project explores the dynamic changes in neuronal network function in the chronic temporal lobe epilepsy (TLE) in rat and human brain in vitro. I focused on the process of establishment of epilepsy (epileptogenesis) in the temporal lobe. Rhythmic behaviour of the hippocampal neuronal networks in healthy animals was explored using spontaneous oscillations in the gamma frequency band (SγO). The use of an improved brain slice preparation technique resulted in the natural occurence (in the absence of pharmacological stimulation) of rhythmic activity, which was then pharmacologically characterised and compared to other models of gamma oscillations (KA- and CCh-induced oscillations) using local field potential recording technique. The results showed that SγO differed from pharmacologically driven models, suggesting higher physiological relevance of SγO. Network activity was also explored in the medial entorhinal cortex (mEC), where spontaneous slow wave oscillations (SWO) were detected. To investigate the course of chronic TLE establishment, a refined Li-pilocarpine-based model of epilepsy (RISE) was developed. The model significantly reduced animal mortality and demonstrated reduced intensity, yet high morbidy with almost 70% mean success rate of developing spontaneous recurrent seizures. We used SγO to characterize changes in the hippocampal neuronal networks throughout the epileptogenesis. The results showed that the network remained largely intact, demonstrating the subtle nature of the RISE model. Despite this, a reduction in network activity was detected during the so-called latent (no seizure) period, which was hypothesized to occur due to network fragmentation and an abnormal function of kainate receptors (KAr). We therefore explored the function of KAr by challenging SγO with kainic acid (KA). The results demonstrated a remarkable decrease in KAr response during the latent period, suggesting KAr dysfunction or altered expression, which will be further investigated using a variety of electrophysiological and immunocytochemical methods. The entorhinal cortex, together with the hippocampus, is known to play an important role in the TLE. Considering this, we investigated neuronal network function of the mEC during epileptogenesis using SWO. The results demonstrated a striking difference in AMPAr function, with possible receptor upregulation or abnormal composition in the early development of epilepsy. Alterations in receptor function inevitably lead to changes in the network function, which may play an important role in the development of epilepsy. Preliminary investigations were made using slices of human brain tissue taken following surgery for intratctable epilepsy. Initial results showed that oscillogenesis could be induced in human brain slices and that such network activity was pharmacologically similar to that observed in rodent brain. Overall, our findings suggest that excitatory glutamatergic transmission is heavily involved in the process of epileptogenesis. Together with other types of receptors, KAr and AMPAr contribute to epilepsy establishment and may be the key to uncovering its mechanism.

Functional astrocyte-neuron lactate shuttle in a human stem cell derived neuronal network

The development of stem cell-derived neuronal networks will promote experimental system development for drug screening, toxicological testing and disease modelling, providing that they mirror closely the functional competencies of their in vivo counterparts. The NT2 cell line is one of the best documented embryocarcinoma cell lines, and can be differentiated into neurons and astrocytes. Great focus has also been placed on defining the electrophysiological properties of these cells, and more recently we have investigated the functional properties of their associated astrocytes. We now show for the first time in a human stem cell derived co-culture model that these cultures are also metabolically competent and demonstrate a functional astrocyte neuron lactate shuttle (ANLS). The ANLS hypothesis proposes that during neuronal activity, glutamate released into the synaptic cleft is taken up by astrocytes and triggers glucose uptake which is converted into lactate and released via monocarboxylate transporters for neuronal use. Using mixed cultures of NT2 derived neurons and astrocytes we have shown that these cells modulate their glucose uptake in response to glutamate, an effect that was blocked by cytochalasin B and ouabain. Additionally we demonstrate that in response to increased neuronal activity and under hypoglycaemic conditions, co-cultures modulate glycogen turnover and increase lactate production. Similar results were also shown following treatment with glutamate, potassium, Isoproterenol and dbcAMP. Together these results demonstrate for the first time a functional ANLS in a human stem cell derived co-culture.