Physiological and pathological oscillatory networks in the human motor system.

Human brain functions are heavily contingent on neural interactions both at the single neuron and the neural population or system level. Accumulating evidence from neurophysiological studies strongly suggests that coupling of oscillatory neural activity provides an important mechanism to establish neural interactions. With the availability of whole-head magnetoencephalography (MEG) macroscopic oscillatory activity can be measured non-invasively from the human brain with high temporal and spatial resolution. To localise, quantify and map oscillatory activity and interactions onto individual brain anatomy we have developed the 'dynamic imaging of coherent sources' (DICS) method which allows to identify and analyse cerebral oscillatory networks from MEG recordings. Using this approach we have characterized physiological and pathological oscillatory networks in the human sensorimotor system. Coherent 8 Hz oscillations emerge from a cerebello-thalamo-premotor-motor cortical network and exert an 8 Hz oscillatory drive on the spinal motor neurons which can be observed as a physiological tremulousness of the movement termed movement discontinuities. This network represents the neurophysiological substrate of a discrete mode of motor control. In parkinsonian resting tremor we have identified an extensive cerebral network consisting of primary motor and lateral premotor cortex, supplementary motor cortex, thalamus/basal ganglia, posterior parietal cortex and secondary somatosensory cortex, which are entrained in the tremor or twice the tremor rhythm. This low frequency entrapment of motor areas likely plays an important role in the pathophysiology of parkinsonian motor symptoms. Finally, studies on patients with postural tremor in hepatic encephalopathy revealed that this type of tremor results from a pathologically slow thalamocortical and cortico-muscular coupling during isometric hold tasks. In conclusion, the analysis of oscillatory cerebral networks provides new insights into physiological mechanisms of motor control and pathophysiological mechanisms of tremor disorders.     

Spontaneous Local Gamma Oscillation Selectively Enhances Neural Network Responsiveness

Synchronized oscillation is very commonly observed in many neuronal systems and might play an important role in the response properties of the system. We have studied how the spontaneous oscillatory activity affects the responsiveness of a neuronal network, using a neural network model of the visual cortex built from Hodgkin-Huxley type excitatory (E-) and inhibitory (I-) neurons. When the isotropic local E-I and I-E synaptic connections were sufficiently strong, the network commonly generated gamma frequency oscillatory firing patterns in response to random feed-forward (FF) input spikes. This spontaneous oscillatory network activity injects a periodic local current that could amplify a weak synaptic input and enhance the network's responsiveness. When E-E connections were added, we found that the strength of oscillation can be modulated by varying the FF input strength without any changes in single neuron properties or interneuron connectivity. The response modulation is proportional to the oscillation strength, which leads to self-regulation such that the cortical network selectively amplifies various FF inputs according to its strength, without requiring any adaptation mechanism. We show that this selective cortical amplification is controlled by E-E cell interactions. We also found that this response amplification is spatially localized, which suggests that the responsiveness modulation may also be spatially selective. This suggests a generalized mechanism by which neural oscillatory activity can enhance the selectivity of a neural network to FF inputs.     

Synaptic interactions mediating synchrony and oscillations in primate sensorimotor cortex.

The appearance of oscillatory modes of 'gamma' activity in many cortical areas of different species has generated interest in understanding their underlying mechanisms and possible functions. This paper reviews evidence from studies on primate motor cortex showing that oscillatory activity entrains many neurons during periods of exploratory manipulative behavior. These oscillatory episodes synchronize widely spread neurons in sensorimotor cortex bilaterally, including descending corticospinal neurons, as evidenced by correlated modulations in EMG activity. The resulting neural synchronization involves task-related and -unrelated neurons similarly, suggesting that it is more likely to play some global role in attention than mediating any obvious interactions involved in coordinating movements. Intracellular recordings have elucidated the strength and types of synaptic interactions between motor cortical neurons that are involved in both normal and oscillatory activity. Spike-triggered averages (STAs) of intracellular membrane potentials have revealed serial connections in the form of unitary excitatory and inhibitory post-synaptic potentials (EPSPs and IPSPs). More commonly, STAs showed large synchronous excitatory or inhibitory potentials (ASEPs and ASIPs) beginning before the trigger spike and composed of multiple unitary events. ASEPs involved synchronous activity in a larger and more widespread group of presynaptic neurons than ASIPs. During oscillatory episodes synchronized excitatory and inhibitory synaptic potentials occurred in varying proportions. EPSPs evoked by stimulating neighboring cortical sites during the depolarizing phase of spontaneous oscillations showed evidence of transient potentiation. These observations are consistent with several functional hypotheses, but fit best with a possible role in attention or arousal.     

Bulbocortical interplay in olfactory information processing via synchronous oscillations

Emergence of synchronous oscillatory activity is an inherent feature of the olfactory systems of insects, mollusks and mammals. A class of simple computational models of the mammalian olfactory system consisting of olfactory bulb and olfactory cortex is constructed to explore possible roles of the related neural circuitry in olfactory information processing via synchronous oscillations. In the models, the bulbar neural circuitry is represented by a chain of oscillators and that of cortex is analogous to an associative memory network with horizontal synaptic connections. The models incorporate the backprojection from cortical units to the bulbar oscillators in particular ways. They exhibit rapid and robust synchronous oscillations in the presence of odorant stimuli, while they show either nonoscillatory states or propagating waves in the absence of stimuli, depending on the values of model parameters. In both models, the backprojection is shown to enhance the establishment of large-scale synchrony. The results suggest that the modulation of neural activity through centrifugal inputs may play an important role at the early stage of cortical information processing.     

Development of neural circuitry for precise temporal sequences through spontaneous activity, axon remodeling, and synaptic plasticity

Temporally precise sequences of neuronal spikes that span hundreds of milliseconds are observed in many brain areas, including songbird premotor nucleus, cat visual cortex, and primary motor cortex. Synfire chains-networks in which groups of neurons are connected via excitatory synapses into a unidirectional chain-are thought to underlie the generation of such sequences. It is unknown, however, how synfire chains can form in local neural circuits, especially for long chains. Here, we show through computer simulation that long synfire chains can develop through spike-time dependent synaptic plasticity and axon remodeling-the pruning of prolific weak connections that follows the emergence of a finite number of strong connections. The formation process begins with a random network. A subset of neurons, called training neurons, intermittently receive superthreshold external input. Gradually, a synfire chain emerges through a recruiting process, in which neurons within the network connect to the tail of the chain started by the training neurons. The model is robust to varying parameters, as well as natural events like neuronal turnover and massive lesions. Our model suggests that long synfire chain can form during the development through self-organization, and axon remodeling, ubiquitous in developing neural circuits, is essential in the process.     

Bilaterally propagating waves of spontaneous activity arising from discrete pacemakers in the neonatal mouse cerebral cortex

Spontaneous electrical activity that moves in synchronized waves across large populations of neurons plays widespread and important roles in nervous system development. The propagation patterns of such waves can encode the spatial location of neurons to their downstream targets and strengthen synaptic connections in coherent spatial patterns. Such waves can arise as an emergent property of mutually excitatory neural networks, or can be driven by a discrete pacemaker. In the mouse cerebral cortex, spontaneous synchronized activity occurs for approximately 72 h of development centered on the day of birth. It is not known whether this activity is driven by a discrete pacemaker or occurs as an emergent network property. Here we show that this activity propagates as a wave that is initiated at either of two homologous pacemakers in the temporal region, and then propagates rapidly across both sides of the brain. When these regions of origin are surgically isolated, waves do not occur. Therefore, this cortical spontaneous activity is a bilateral wave that originates from a discrete subset of pacemaker neurons. © 2009 Wiley Periodicals, Inc. Develop Neurobiol, 2009     

Spatially structured oscillations in a two-dimensional excitatory neuronal network with synaptic depression

We study the spatiotemporal dynamics of a two-dimensional excitatory neuronal network with synaptic depression. Coupling between populations of neurons is taken to be nonlocal, while depression is taken to be local and presynaptic. We show that the network supports a wide range of spatially structured oscillations, which are suggestive of phenomena seen in cortical slice experiments and in vivo. The particular form of the oscillations depends on initial conditions and the level of background noise. Given an initial, spatially localized stimulus, activity evolves to a spatially localized oscillating core that periodically emits target waves. Low levels of noise can spontaneously generate several pockets of oscillatory activity that interact via their target patterns. Periodic activity in space can also organize into spiral waves, provided that there is some source of rotational symmetry breaking due to external stimuli or noise. In the high gain limit, no oscillatory behavior exists, but a transient stimulus can lead to a single, outward propagating target wave.     

Cortical Resonance Frequencies Emerge from Network Size and Connectivity

Neural oscillations occur within a wide frequency range with different brain regions exhibiting resonance-like characteristics at specific points in the spectrum. At the microscopic scale, single neurons possess intrinsic oscillatory properties, such that is not yet known whether cortical resonance is consequential to neural oscillations or an emergent property of the networks that interconnect them. Using a network model of loosely-coupled Wilson-Cowan oscillators to simulate a patch of cortical sheet, we demonstrate that the size of the activated network is inversely related to its resonance frequency. Further analysis of the parameter space indicated that the number of excitatory and inhibitory connections, as well as the average transmission delay between units, determined the resonance frequency. The model predicted that if an activated network within the visual cortex increased in size, the resonance frequency of the network would decrease. We tested this prediction experimentally using the steady-state visual evoked potential where we stimulated the visual cortex with different size stimuli at a range of driving frequencies. We demonstrate that the frequency corresponding to peak steady-state response inversely correlated with the size of the network. We conclude that although individual neurons possess resonance properties, oscillatory activity at the macroscopic level is strongly influenced by network interactions, and that the steady-state response can be used to investigate functional networks.     

A thalamo-cortical neural mass model for the simulation of brain rhythms during sleep

Cortico-thalamic interactions are known to play a pivotal role in many brain phenomena, including sleep, attention, memory consolidation and rhythm generation. Hence, simple mathematical models that can simulate the dialogue between the cortex and the thalamus, at a mesoscopic level, have a great cognitive value. In the present work we describe a neural mass model of a cortico-thalamic module, based on neurophysiological mechanisms. The model includes two thalamic populations (a thalamo-cortical relay cell population, TCR, and its related thalamic reticular nucleus, TRN), and a cortical column consisting of four connected populations (pyramidal neurons, excitatory interneurons, inhibitory interneurons with slow and fast kinetics). Moreover, thalamic neurons exhibit two firing modes: bursting and tonic. Finally, cortical synapses among pyramidal neurons incorporate a disfacilitation mechanism following prolonged activity. Simulations show that the model is able to mimic the different patterns of rhythmic activity in cortical and thalamic neurons (beta and alpha waves, spindles, delta waves, K-complexes, slow sleep waves) and their progressive changes from wakefulness to deep sleep, by just acting on modulatory inputs. Moreover, simulations performed by providing short sensory inputs to the TCR show that brain rhythms during sleep preserve the cortex from external perturbations, still allowing a high cortical activity necessary to drive synaptic plasticity and memory consolidation. In perspective, the present model may be used within larger cortico-thalamic networks, to gain a deeper understanding of mechanisms beneath synaptic changes during sleep, to investigate the specific role of brain rhythms, and to explore cortical synchronization achieved via thalamic influences.     

Generating oscillatory bursts from a network of regular spiking neurons without inhibition

Avian nucleus isthmi pars parvocellularis (Ipc) neurons are reciprocally connected with the layer 10 (L10) neurons in the optic tectum and respond with oscillatory bursts to visual stimulation. Our in vitro experiments show that both neuron types respond with regular spiking to somatic current injection and that the feedforward and feedback synaptic connections are excitatory, but of different strength and time course. To elucidate mechanisms of oscillatory bursting in this network of regularly spiking neurons, we investigated an experimentally constrained model of coupled leaky integrate-and-fire neurons with spike-rate adaptation. The model reproduces the observed Ipc oscillatory bursting in response to simulated visual stimulation. A scan through the model parameter volume reveals that Ipc oscillatory burst generation can be caused by strong and brief feedforward synaptic conductance changes. The mechanism is sensitive to the parameter values of spike-rate adaptation. In conclusion, we show that a network of regular-spiking neurons with feedforward excitation and spike-rate adaptation can generate oscillatory bursting in response to a constant input.     

Emergence of coherent traveling waves controlling quadruped gaits in a two-dimensional spinal cord model

 The concerted and self-organizing behavior of spinal cord segments in generating locomotor patterns is modulated by afferent sensory information and controlled by descending pathways from the brainstem, cerebellum, or cortex. The purpose of this study was to define a minimal set of parameters that could control a similar self-organizing behavior in a two-dimensional neural network. When we implemented synaptic depression and active membrane repolarization as two properties of the neurons, the two-dimensional neural network generated traveling waves. Their wavelength and angle of propagation could be independently controlled by two parameters that modulated excitatory premotor neurons and inhibitory commissural neurons. It is further demonstrated that the selection of wave parameters corresponds to the selection of quadruped gaits. Received: 30 July 2001 / Accepted in revised form: 17 April 2002 Correspondence to: A. Kaske (e-mail: alexander.kaske@mtc.ki.se, alexander.kaske@vglab.com)     

Computational Modeling of Seizure Dynamics Using Coupled Neuronal Networks: Factors Shaping Epileptiform Activity

Epileptic seizure dynamics span multiple scales in space and time. Understanding seizure mechanisms requires identifying the relations between seizure components within and across these scales, together with the analysis of their dynamical repertoire. Mathematical models have been developed to reproduce seizure dynamics across scales ranging from the single neuron to the neural population. In this study, we develop a network model of spiking neurons and systematically investigate the conditions, under which the network displays the emergent dynamic behaviors known from the Epileptor, which is a well-investigated abstract model of epileptic neural activity. This approach allows us to study the biophysical parameters and variables leading to epileptiform discharges at cellular and network levels. Our network model is composed of two neuronal populations, characterized by fast excitatory bursting neurons and regular spiking inhibitory neurons, embedded in a common extracellular environment represented by a slow variable. By systematically analyzing the parameter landscape offered by the simulation framework, we reproduce typical sequences of neural activity observed during status epilepticus. We find that exogenous fluctuations from extracellular environment and electro-tonic couplings play a major role in the progression of the seizure, which supports previous studies and further validates our model. We also investigate the influence of chemical synaptic coupling in the generation of spontaneous seizure-like events. Our results argue towards a temporal shift of typical spike waves with fast discharges as synaptic strengths are varied. We demonstrate that spike waves, including interictal spikes, are generated primarily by inhibitory neurons, whereas fast discharges during the wave part are due to excitatory neurons. Simulated traces are compared with in vivo experimental data from rodents at different stages of the disorder. We draw the conclusion that slow variations of global excitability, due to exogenous fluctuations from extracellular environment, and gap junction communication push the system into paroxysmal regimes. We discuss potential mechanisms underlying such machinery and the relevance of our approach, supporting previous detailed modeling studies and reflecting on the limitations of our methodology.     

Synconset Waves and Chains: Spiking Onsets in Synchronous Populations Predict and Are Predicted by Network Structure

Synfire waves are propagating spike packets in synfire chains, which are feedforward chains embedded in random networks. Although synfire waves have proved to be effective quantification for network activity with clear relations to network structure, their utilities are largely limited to feedforward networks with low background activity. To overcome these shortcomings, we describe a novel generalisation of synfire waves, and define ‘synconset wave’ as a cascade of first spikes within a synchronisation event. Synconset waves would occur in ‘synconset chains’, which are feedforward chains embedded in possibly heavily recurrent networks with heavy background activity. We probed the utility of synconset waves using simulation of single compartment neuron network models with biophysically realistic conductances, and demonstrated that the spread of synconset waves directly follows from the network connectivity matrix and is modulated by top-down inputs and the resultant oscillations. Such synconset profiles lend intuitive insights into network organisation in terms of connection probabilities between various network regions rather than an adjacency matrix. To test this intuition, we develop a Bayesian likelihood function that quantifies the probability that an observed synfire wave was caused by a given network. Further, we demonstrate it''s utility in the inverse problem of identifying the network that caused a given synfire wave. This method was effective even in highly subsampled networks where only a small subset of neurons were accessible, thus showing it''s utility in experimental estimation of connectomes in real neuronal-networks. Together, we propose synconset chains/waves as an effective framework for understanding the impact of network structure on function, and as a step towards developing physiology-driven network identification methods. Finally, as synconset chains extend the utilities of synfire chains to arbitrary networks, we suggest utilities of our framework to several aspects of network physiology including cell assemblies, population codes, and oscillatory synchrony.     

Peripheral and central inputs shape network dynamics in the developing visual cortex in vivo

Spontaneous network activity constitutes a central theme during the development of neuronal circuitry [1, 2]. Before the onset of vision, retinal neurons generate waves of spontaneous activity that are relayed along the ascending visual pathway [3, 4] and shape activity patterns in these regions [5, 6]. The spatiotemporal nature of retinal waves is required to establish precise functional maps in higher visual areas, and their disruption results in enlarged axonal projection areas (e.g., [7-10]). However, how retinal inputs shape network dynamics in the visual cortex on the cellular level is unknown. Using in vivo two-photon calcium imaging, we identified two independently occurring patterns of network activity in the mouse primary visual cortex (V1) before and at the onset of vision. Acute manipulations of spontaneous retinal activity revealed that one type of network activity largely originated in the retina and was characterized by low synchronicity (L-) events. In addition, we identified a type of high synchronicity (H-) events that required gap junction signaling but were independent of retinal input. Moreover, the patterns differed in wave progression and developmental profile. Our data suggest that different activity patterns have complementary functions during the formation of synaptic circuits in the developing visual cortex.     

Wave-Processing of Long-Scale Information by Neuronal Chains

Investigation of mechanisms of information handling in neural assemblies involved in computational and cognitive tasks is a challenging problem. Synergetic cooperation of neurons in time domain, through synchronization of firing of multiple spatially distant neurons, has been widely spread as the main paradigm. Complementary, the brain may also employ information coding and processing in spatial dimension. Then, the result of computation depends also on the spatial distribution of long-scale information. The latter bi-dimensional alternative is notably less explored in the literature. Here, we propose and theoretically illustrate a concept of spatiotemporal representation and processing of long-scale information in laminar neural structures. We argue that relevant information may be hidden in self-sustained traveling waves of neuronal activity and then their nonlinear interaction yields efficient wave-processing of spatiotemporal information. Using as a testbed a chain of FitzHugh-Nagumo neurons, we show that the wave-processing can be achieved by incorporating into the single-neuron dynamics an additional voltage-gated membrane current. This local mechanism provides a chain of such neurons with new emergent network properties. In particular, nonlinear waves as a carrier of long-scale information exhibit a variety of functionally different regimes of interaction: from complete or asymmetric annihilation to transparent crossing. Thus neuronal chains can work as computational units performing different operations over spatiotemporal information. Exploiting complexity resonance these composite units can discard stimuli of too high or too low frequencies, while selectively compress those in the natural frequency range. We also show how neuronal chains can contextually interpret raw wave information. The same stimulus can be processed differently or identically according to the context set by a periodic wave train injected at the opposite end of the chain.     

Formation and Dynamics of Waves in a Cortical Model of Cholinergic Modulation

Acetylcholine (ACh) is a regulator of neural excitability and one of the neurochemical substrates of sleep. Amongst the cellular effects induced by cholinergic modulation are a reduction in spike-frequency adaptation (SFA) and a shift in the phase response curve (PRC). We demonstrate in a biophysical model how changes in neural excitability and network structure interact to create three distinct functional regimes: localized asynchronous, traveling asynchronous, and traveling synchronous. Our results qualitatively match those observed experimentally. Cortical activity during slow wave sleep (SWS) differs from that during REM sleep or waking states. During SWS there are traveling patterns of activity in the cortex; in other states stationary patterns occur. Our model is a network composed of Hodgkin-Huxley type neurons with a M-current regulated by ACh. Regulation of ACh level can account for dynamical changes between functional regimes. Reduction of the magnitude of this current recreates the reduction in SFA the shift from a type 2 to a type 1 PRC observed in the presence of ACh. When SFA is minimal (in waking or REM sleep state, high ACh) patterns of activity are localized and easily pinned by network inhomogeneities. When SFA is present (decreasing ACh), traveling waves of activity naturally arise. A further decrease in ACh leads to a high degree of synchrony within traveling waves. We also show that the level of ACh determines how sensitive network activity is to synaptic heterogeneity. These regimes may have a profound functional significance as stationary patterns may play a role in the proper encoding of external input as memory and traveling waves could lead to synaptic regularization, giving unique insights into the role and significance of ACh in determining patterns of cortical activity and functional differences arising from the patterns.     

In vivo and in vitro patch-clamp recording analysis of the process of sensory transmission in the spinal cord and sensory cortex

Following the integration and modification of the sensory inputs in the spinal cord, the information is transmitted to the primary sensory cortex where the integrated information is further processed and perceived. Processing of the sensory information in the spinal cord has been intensively investigated. However, the mechanisms of how the inputs are processed in the cortex are still unclear. To know the correlation of the sensory processing in the dorsal horn and cortex, in vivo and in vitro patch-clamp recordings were made from rat dorsal horn and sensory cortex. Although dorsal horn neurons showed spontaneous and evoked EPSCs by noxious and non-noxious stimuli, most somatosensory neurons located at 100 to 1000 microm from the surface of the cortex exhibited an oscillatory activity and received synaptic inputs from non-noxious but not noxious receptors. These observations suggest that the synaptic responses in cortical neurons are processed in a more complex manner; and this may be due to the reciprocal synaptic connection between thalamus and cortex.     

Neural network simulations of coupled locomotor oscillators in the lamprey spinal cord

The segmental locomotor network in the lamprey spinal cord was simulated on a computer using a connectionist-type neural network. The cells of the network were identical except for their excitatory levels and their synaptic connections. The synaptic connections used were based on previous experimental work. It was demonstrated that the connectivity of the circuit is capable of generating oscillatory activity with the appropriate phase relations among the cells. Intersegmental coordination was explored by coupling two identical segmental networks using only the cells of the network. Each of the possible couplings of a bilateral pair of cells in one oscillator with a bilateral pair of cells in the other oscillator produced stable phase locking of the two oscillators. The degree of phase difference was dependent upon synaptic weight, and the operating range of synaptic weights varied among the pairs of connections. The coupling was tested using several criteria from experimental work on the lamprey spinal cord. Coupling schemes involving several pairs of connecting cells were found which 1) achieved steadystate phase locking within a single cycle, 2) exhibited constant phase differences over a wide range of cycle periods, and 3) maintained stable phase locking in spite of large differences in the intrinsic frequencies of the two oscillators. It is concluded that the synaptic connectivity plays a large role in producing oscillations in this network and that it is not necessary to postulate a separate set of coordinating neurons between oscillators in order to achieve appropriate phase coupling.     

Modeling of Rhythmogenesis in an Ensemble of Cortical Neurons Connected with Electrical Synapses

Based on our own data on generation of spindle-like field electrical activity in neuronal barrels of the rat somatic cortex and also on the published data on the properties of voltage-dependent channels in the membranes of cortical cells, we developed a model of the ensemble (simple network) of neurons connected by electrical synapses. Such connections were found earlier in neurophysiological and ultramicroscopic studies. Model neurons with membranes having sodium, potassium, and calcium channels described in the literature were capable of generating bursting rhythmic impulse activity under conditions of switching off of synaptic connections between cells (isolation). With switching on of electrical synapses, spiking generated by separate neurons, which initially was nonsynchronous, became synchronized in time. Ipso facto, we demonstrated the ability of pacemaker oscillatory activity to be electrotonically synchronized in ensembles of neurons connected with electrical synapses.     

M-channels modulate the intrinsic excitability and synaptic responses of layer 2/3 pyramidal neurons in auditory cortex

Neurons in the auditory cortex are believed to utilize temporal patterns of neural activity to accurately process auditory information but the intrinsic neuronal mechanism underlying the control of auditory neural activity is not known. The slowly activating, persistent K⁺ channel, also called M-channel that belongs to the Kv7 family, is already known to be important in regulating subthreshold neural excitability and synaptic summation in neocortical and hippocampal pyramidal neurons. However, its functional role in the primary auditory cortex (A1) has never been characterized. In this study, we investigated the roles of M-channels on neuronal excitability, short-term plasticity, and synaptic summation of A1 layer 2/3 regular spiking pyramidal neurons with whole-cell current-clamp recordings in vitro. We found that blocking M-channels with a selective M-channel blocker, XE991, significantly increased neural excitability of A1 layer 2/3 pyramidal neurons. Furthermore, M-channels controled synaptic responses of intralaminar-evoked excitatory postsynaptic potentials (EPSPs); XE991 significantly increased EPSP amplitude, decreased the rate of short-term depression, and increased the synaptic summation. These results suggest that M-channels are involved in controlling spike output patterns and synaptic responses of A1 layer 2/3 pyramidal neurons, which would have important implications in auditory information processing.