Neural Plasticity in Human Brain Connectivity: The Effects of Long Term Deep Brain Stimulation of the Subthalamic Nucleus in Parkinson’s Disease

Background

Positive clinical outcomes are now well established for deep brain stimulation, but little is known about the effects of long-term deep brain stimulation on brain structural and functional connectivity. Here, we used the rare opportunity to acquire pre- and postoperative diffusion tensor imaging in a patient undergoing deep brain stimulation in bilateral subthalamic nuclei for Parkinson’s Disease. This allowed us to analyse the differences in structural connectivity before and after deep brain stimulation. Further, a computational model of spontaneous brain activity was used to estimate the changes in functional connectivity arising from the specific changes in structural connectivity.

Results

We found significant localised structural changes as a result of long-term deep brain stimulation. These changes were found in sensory-motor, prefrontal/limbic, and olfactory brain regions which are known to be affected in Parkinson’s Disease. The nature of these changes was an increase of nodal efficiency in most areas and a decrease of nodal efficiency in the precentral sensory-motor area. Importantly, the computational model clearly shows the impact of deep brain stimulation-induced structural alterations on functional brain changes, which is to shift the neural dynamics back towards a healthy regime. The results demonstrate that deep brain stimulation in Parkinson’s Disease leads to a topological reorganisation towards healthy bifurcation of the functional networks measured in controls, which suggests a potential neural mechanism for the alleviation of symptoms.

Conclusions

The findings suggest that long-term deep brain stimulation has not only restorative effects on the structural connectivity, but also affects the functional connectivity at a global level. Overall, our results support causal changes in human neural plasticity after long-term deep brain stimulation and may help to identify the underlying mechanisms of deep brain stimulation.     

Encoding time in neural dynamic regimes with distinct computational tradeoffs

Converging evidence suggests the brain encodes time in dynamic patterns of neural activity, including neural sequences, ramping activity, and complex dynamics. Most temporal tasks, however, require more than just encoding time, and can have distinct computational requirements including the need to exhibit temporal scaling, generalize to novel contexts, or robustness to noise. It is not known how neural circuits can encode time and satisfy distinct computational requirements, nor is it known whether similar patterns of neural activity at the population level can exhibit dramatically different computational or generalization properties. To begin to answer these questions, we trained RNNs on two timing tasks based on behavioral studies. The tasks had different input structures but required producing identically timed output patterns. Using a novel framework we quantified whether RNNs encoded two intervals using either of three different timing strategies: scaling, absolute, or stimulus-specific dynamics. We found that similar neural dynamic patterns at the level of single intervals, could exhibit fundamentally different properties, including, generalization, the connectivity structure of the trained networks, and the contribution of excitatory and inhibitory neurons. Critically, depending on the task structure RNNs were better suited for generalization or robustness to noise. Further analysis revealed different connection patterns underlying the different regimes. Our results predict that apparently similar neural dynamic patterns at the population level (e.g., neural sequences) can exhibit fundamentally different computational properties in regards to their ability to generalize to novel stimuli and their robustness to noise—and that these differences are associated with differences in network connectivity and distinct contributions of excitatory and inhibitory neurons. We also predict that the task structure used in different experimental studies accounts for some of the experimentally observed variability in how networks encode time.     

Leader neurons in leaky integrate and fire neural network simulations

In this paper, we highlight the topological properties of leader neurons whose existence is an experimental fact. Several experimental studies show the existence of leader neurons in population bursts of activity in 2D living neural networks (Eytan and Marom, J Neurosci 26(33):8465–8476, 2006; Eckmann et al., New J Phys 10(015011), 2008). A leader neuron is defined as a neuron which fires at the beginning of a burst (respectively network spike) more often than we expect by chance considering its mean firing rate. This means that leader neurons have some burst triggering power beyond a chance-level statistical effect. In this study, we characterize these leader neuron properties. This naturally leads us to simulate neural 2D networks. To build our simulations, we choose the leaky integrate and fire (lIF) neuron model (Gerstner and Kistler 2002; Cessac, J Math Biol 56(3):311–345, 2008), which allows fast simulations (Izhikevich, IEEE Trans Neural Netw 15(5):1063–1070, 2004; Gerstner and Naud, Science 326:379–380, 2009). The dynamics of our lIF model has got stable leader neurons in the burst population that we simulate. These leader neurons are excitatory neurons and have a low membrane potential firing threshold. Except for these two first properties, the conditions required for a neuron to be a leader neuron are difficult to identify and seem to depend on several parameters involved in the simulations themselves. However, a detailed linear analysis shows a trend of the properties required for a neuron to be a leader neuron. Our main finding is: A leader neuron sends signals to many excitatory neurons as well as to few inhibitory neurons and a leader neuron receives only signals from few other excitatory neurons. Our linear analysis exhibits five essential properties of leader neurons each with different relative importance. This means that considering a given neural network with a fixed mean number of connections per neuron, our analysis gives us a way of predicting which neuron is a good leader neuron and which is not. Our prediction formula correctly assesses leadership for at least ninety percent of neurons.     

Emergent functional properties of neuronal networks with controlled topology

The interplay between anatomical connectivity and dynamics in neural networks plays a key role in the functional properties of the brain and in the associated connectivity changes induced by neural diseases. However, a detailed experimental investigation of this interplay at both cellular and population scales in the living brain is limited by accessibility. Alternatively, to investigate the basic operational principles with morphological, electrophysiological and computational methods, the activity emerging from large in vitro networks of primary neurons organized with imposed topologies can be studied. Here, we validated the use of a new bio-printing approach, which effectively maintains the topology of hippocampal cultures in vitro and investigated, by patch-clamp and MEA electrophysiology, the emerging functional properties of these grid-confined networks. In spite of differences in the organization of physical connectivity, our bio-patterned grid networks retained the key properties of synaptic transmission, short-term plasticity and overall network activity with respect to random networks. Interestingly, the imposed grid topology resulted in a reinforcement of functional connections along orthogonal directions, shorter connectivity links and a greatly increased spiking probability in response to focal stimulation. These results clearly demonstrate that reliable functional studies can nowadays be performed on large neuronal networks in the presence of sustained changes in the physical network connectivity.     

The impact of sparsity in low-rank recurrent neural networks

Neural population dynamics are often highly coordinated, allowing task-related computations to be understood as neural trajectories through low-dimensional subspaces. How the network connectivity and input structure give rise to such activity can be investigated with the aid of low-rank recurrent neural networks, a recently-developed class of computational models which offer a rich theoretical framework linking the underlying connectivity structure to emergent low-dimensional dynamics. This framework has so far relied on the assumption of all-to-all connectivity, yet cortical networks are known to be highly sparse. Here we investigate the dynamics of low-rank recurrent networks in which the connections are randomly sparsified, which makes the network connectivity formally full-rank. We first analyse the impact of sparsity on the eigenvalue spectrum of low-rank connectivity matrices, and use this to examine the implications for the dynamics. We find that in the presence of sparsity, the eigenspectra in the complex plane consist of a continuous bulk and isolated outliers, a form analogous to the eigenspectra of connectivity matrices composed of a low-rank and a full-rank random component. This analogy allows us to characterise distinct dynamical regimes of the sparsified low-rank network as a function of key network parameters. Altogether, we find that the low-dimensional dynamics induced by low-rank connectivity structure are preserved even at high levels of sparsity, and can therefore support rich and robust computations even in networks sparsified to a biologically-realistic extent.     

Brain rhythms reveal a hierarchical network organization

Recordings of ongoing neural activity with EEG and MEG exhibit oscillations of specific frequencies over a non-oscillatory background. The oscillations appear in the power spectrum as a collection of frequency bands that are evenly spaced on a logarithmic scale, thereby preventing mutual entrainment and cross-talk. Over the last few years, experimental, computational and theoretical studies have made substantial progress on our understanding of the biophysical mechanisms underlying the generation of network oscillations and their interactions, with emphasis on the role of neuronal synchronization. In this paper we ask a very different question. Rather than investigating how brain rhythms emerge, or whether they are necessary for neural function, we focus on what they tell us about functional brain connectivity. We hypothesized that if we were able to construct abstract networks, or "virtual brains", whose dynamics were similar to EEG/MEG recordings, those networks would share structural features among themselves, and also with real brains. Applying mathematical techniques for inverse problems, we have reverse-engineered network architectures that generate characteristic dynamics of actual brains, including spindles and sharp waves, which appear in the power spectrum as frequency bands superimposed on a non-oscillatory background dominated by low frequencies. We show that all reconstructed networks display similar topological features (e.g. structural motifs) and dynamics. We have also reverse-engineered putative diseased brains (epileptic and schizophrenic), in which the oscillatory activity is altered in different ways, as reported in clinical studies. These reconstructed networks show consistent alterations of functional connectivity and dynamics. In particular, we show that the complexity of the network, quantified as proposed by Tononi, Sporns and Edelman, is a good indicator of brain fitness, since virtual brains modeling diseased states display lower complexity than virtual brains modeling normal neural function. We finally discuss the implications of our results for the neurobiology of health and disease.     

Cooperation prevails when individuals adjust their social ties

Conventional evolutionary game theory predicts that natural selection favours the selfish and strong even though cooperative interactions thrive at all levels of organization in living systems. Recent investigations demonstrated that a limiting factor for the evolution of cooperative interactions is the way in which they are organized, cooperators becoming evolutionarily competitive whenever individuals are constrained to interact with few others along the edges of networks with low average connectivity. Despite this insight, the conundrum of cooperation remains since recent empirical data shows that real networks exhibit typically high average connectivity and associated single-to-broad–scale heterogeneity. Here, a computational model is constructed in which individuals are able to self-organize both their strategy and their social ties throughout evolution, based exclusively on their self-interest. We show that the entangled evolution of individual strategy and network structure constitutes a key mechanism for the sustainability of cooperation in social networks. For a given average connectivity of the population, there is a critical value for the ratio W between the time scales associated with the evolution of strategy and of structure above which cooperators wipe out defectors. Moreover, the emerging social networks exhibit an overall heterogeneity that accounts very well for the diversity of patterns recently found in acquired data on social networks. Finally, heterogeneity is found to become maximal when W reaches its critical value. These results show that simple topological dynamics reflecting the individual capacity for self-organization of social ties can produce realistic networks of high average connectivity with associated single-to-broad–scale heterogeneity. On the other hand, they show that cooperation cannot evolve as a result of “social viscosity” alone in heterogeneous networks with high average connectivity, requiring the additional mechanism of topological co-evolution to ensure the survival of cooperative behaviour.     

Neural networks for perceptual processing: from simulation tools to theories

Neural networks are modelling tools that are, in principle, able to capture the input-output behaviour of arbitrary systems that may include the dynamics of animal populations or brain circuits. While a neural network model is useful if it captures phenomenologically the behaviour of the target system in this way, its utility is amplified if key mechanisms of the model can be discovered, and identified with those of the underlying system. In this review, we first describe, at a fairly high level with minimal mathematics, some of the tools used in constructing neural network models. We then go on to discuss the implications of network models for our understanding of the system they are supposed to describe, paying special attention to those models that deal with neural circuits and brain systems. We propose that neural nets are useful for brain modelling if they are viewed in a wider computational framework originally devised by Marr. Here, neural networks are viewed as an intermediate mechanistic abstraction between 'algorithm' and 'implementation', which can provide insights into biological neural representations and their putative supporting architectures.     

Dynamic changes in brain lateralization correlate with human cognitive performance

Hemispheric lateralization constitutes a core architectural principle of human brain organization underlying cognition, often argued to represent a stable, trait-like feature. However, emerging evidence underlines the inherently dynamic nature of brain networks, in which time-resolved alterations in functional lateralization remain uncharted. Integrating dynamic network approaches with the concept of hemispheric laterality, we map the spatiotemporal architecture of whole-brain lateralization in a large sample of high-quality resting-state fMRI data (N = 991, Human Connectome Project). We reveal distinct laterality dynamics across lower-order sensorimotor systems and higher-order associative networks. Specifically, we expose 2 aspects of the laterality dynamics: laterality fluctuations (LF), defined as the standard deviation of laterality time series, and laterality reversal (LR), referring to the number of zero crossings in laterality time series. These 2 measures are associated with moderate and extreme changes in laterality over time, respectively. While LF depict positive association with language function and cognitive flexibility, LR shows a negative association with the same cognitive abilities. These opposing interactions indicate a dynamic balance between intra and interhemispheric communication, i.e., segregation and integration of information across hemispheres. Furthermore, in their time-resolved laterality index, the default mode and language networks correlate negatively with visual/sensorimotor and attention networks, which are linked to better cognitive abilities. Finally, the laterality dynamics are associated with functional connectivity changes of higher-order brain networks and correlate with regional metabolism and structural connectivity. Our results provide insights into the adaptive nature of the lateralized brain and new perspectives for future studies of human cognition, genetics, and brain disorders.

Hemispheric lateralization constitutes a core architectural principle of human brain organization, often argued to represent a stable, trait-like feature, but how does this fit with our increasing appreciation of the inherently dynamic nature of brain networks? This neuroimaging study reveals the dynamic nature of functional brain lateralization at resting-state and its relationship with language function and cognitive flexibility.     

Dynamic large-scale connectivity of intrinsic cortical oscillations supports adaptive listening in challenging conditions

In multi-talker situations, individuals adapt behaviorally to this listening challenge mostly with ease, but how do brain neural networks shape this adaptation? We here establish a long-sought link between large-scale neural communications in electrophysiology and behavioral success in the control of attention in difficult listening situations. In an age-varying sample of N = 154 individuals, we find that connectivity between intrinsic neural oscillations extracted from source-reconstructed electroencephalography is regulated according to the listener’s goal during a challenging dual-talker task. These dynamics occur as spatially organized modulations in power-envelope correlations of alpha and low-beta neural oscillations during approximately 2-s intervals most critical for listening behavior relative to resting-state baseline. First, left frontoparietal low-beta connectivity (16 to 24 Hz) increased during anticipation and processing of a spatial-attention cue before speech presentation. Second, posterior alpha connectivity (7 to 11 Hz) decreased during comprehension of competing speech, particularly around target-word presentation. Connectivity dynamics of these networks were predictive of individual differences in the speed and accuracy of target-word identification, respectively, but proved unconfounded by changes in neural oscillatory activity strength. Successful adaptation to a listening challenge thus latches onto two distinct yet complementary neural systems: a beta-tuned frontoparietal network enabling the flexible adaptation to attentive listening state and an alpha-tuned posterior network supporting attention to speech.

This study investigates how intrinsic neural oscillations, acting in concert, tune into attentive listening. Using electroencephalography signals collected from people in a dual-talker listening task, the authors find that network connectivity of frontoparietal beta and posterior alpha oscillations is regulated according to the listener’s goal.     

Intrinsic Brain Network Abnormalities in Migraines without Aura Revealed in Resting-State fMRI

Background

Previous studies have defined low-frequency, spatially consistent intrinsic connectivity networks (ICN) in resting functional magnetic resonance imaging (fMRI) data which reflect functional interactions among distinct brain areas. We sought to explore whether and how repeated migraine attacks influence intrinsic brain connectivity, as well as how activity in these networks correlates with clinical indicators of migraine.

Methods/Principal Findings

Resting-state fMRI data in twenty-three patients with migraines without aura (MwoA) and 23 age- and gender-matched healthy controls (HC) were analyzed using independent component analysis (ICA), in combination with a “dual-regression” technique to identify the group differences of three important pain-related networks [default mode network (DMN), bilateral central executive network (CEN), salience network (SN)] between the MwoA patients and HC. Compared with the HC, MwoA patients showed aberrant intrinsic connectivity within the bilateral CEN and SN, and greater connectivity between both the DMN and right CEN (rCEN) and the insula cortex - a critical region involving in pain processing. Furthermore, greater connectivity between both the DMN and rCEN and the insula correlated with duration of migraine.

Conclusions

Our findings may provide new insights into the characterization of migraine as a condition affecting brain activity in intrinsic connectivity networks. Moreover, the abnormalities may be the consequence of a persistent central neural system dysfunction, reflecting cumulative brain insults due to frequent ongoing migraine attacks.     

Shaping bursting by electrical coupling and noise

Gap-junctional coupling is an important way of communication between neurons and other excitable cells. Strong electrical coupling synchronizes activity across cell ensembles. Surprisingly, in the presence of noise synchronous oscillations generated by an electrically coupled network may differ qualitatively from the oscillations produced by uncoupled individual cells forming the network. A prominent example of such behavior is the synchronized bursting in islets of Langerhans formed by pancreatic β-cells, which in isolation are known to exhibit irregular spiking (Sherman and Rinzel, Biophys J 54:411–425, 1988; Sherman and Rinzel, Biophys J 59:547–559, 1991). At the heart of this intriguing phenomenon lies denoising, a remarkable ability of electrical coupling to diminish the effects of noise acting on individual cells. In this paper, building on an earlier analysis of denoising in networks of integrate-and-fire neurons (Medvedev, Neural Comput 21 (11):3057–3078, 2009) and our recent study of spontaneous activity in a closely related model of the Locus Coeruleus network (Medvedev and Zhuravytska, The geometry of spontaneous spiking in neuronal networks, submitted, 2012), we derive quantitative estimates characterizing denoising in electrically coupled networks of conductance-based models of square wave bursting cells. Our analysis reveals the interplay of the intrinsic properties of the individual cells and network topology and their respective contributions to this important effect. In particular, we show that networks on graphs with large algebraic connectivity (Fiedler, Czech Math J 23(98):298–305, 1973) or small total effective resistance (Bollobas, Modern graph theory, Graduate Texts in Mathematics, vol. 184, Springer, New York, 1998) are better equipped for implementing denoising. As a by-product of the analysis of denoising, we analytically estimate the rate with which trajectories converge to the synchronization subspace and the stability of the latter to random perturbations. These estimates reveal the role of the network topology in synchronization. The analysis is complemented by numerical simulations of electrically coupled conductance-based networks. Taken together, these results explain the mechanisms underlying synchronization and denoising in an important class of biological models.     

Induction and propagation of transient synchronous activity in neural networks endowed with short-term plasticity

Transient, task related synchronous activity within neural populations has been recognized as the substrate of temporal coding in the brain. The mechanisms underlying inducing and propagation of transient synchronous activity are still unknown, and we propose that short-term plasticity (STP) of neural circuits may serve as a supplemental mechanism therein. By computational modeling, we showed that short-term facilitation greatly increases the reactivation rate of population spikes and decreases the latency of response to reactivation stimuli in local recurrent neural networks. Meanwhile, the timing of population spike reactivation is controlled by the memory effect of STP, and it is mediated primarily by the facilitation time constant. Furthermore, we demonstrated that synaptic facilitation dramatically enhances synchrony propagation in feedforward neural networks and that response timing mediated by synaptic facilitation offers a scheme for information routing. In addition, we verified that synaptic strengthening of intralayer or interlayer coupling enhances synchrony propagation, and we verified that other factors such as the delay of synaptic transmission and the mode of synaptic connectivity are also involved in regulating synchronous activity propagation. Overall, our results highlight the functional role of STP in regulating the inducing and propagation of transient synchronous activity, and they may inspire testable hypotheses for future experimental studies.     

Functional Connectivity in MRI Is Driven by Spontaneous BOLD Events

Functional brain signals are frequently decomposed into a relatively small set of large scale, distributed cortical networks that are associated with different cognitive functions. It is generally assumed that the connectivity of these networks is static in time and constant over the whole network, although there is increasing evidence that this view is too simplistic. This work proposes novel techniques to investigate the contribution of spontaneous BOLD events to the temporal dynamics of functional connectivity as assessed by ultra-high field functional magnetic resonance imaging (fMRI). The results show that: 1) spontaneous events in recognised brain networks contribute significantly to network connectivity estimates; 2) these spontaneous events do not necessarily involve whole networks or nodes, but clusters of voxels which act in concert, forming transiently synchronising sub-networks and 3) a task can significantly alter the number of localised spontaneous events that are detected within a single network. These findings support the notion that spontaneous events are the main driver of the large scale networks that are commonly detected by seed-based correlation and ICA. Furthermore, we found that large scale networks are manifestations of smaller, transiently synchronising sub-networks acting dynamically in concert, corresponding to spontaneous events, and which do not necessarily involve all voxels within the network nodes oscillating in unison.     

Spatio-Temporal Correlation Tensors Reveal Functional Structure in Human Brain

Resting state functional magnetic resonance imaging (fMRI) has been commonly used to measure functional connectivity between cortical regions, while diffusion tensor imaging (DTI) can be used to characterize structural connectivity of white matter tracts. In principle combining resting state fMRI and DTI data could allow characterization of structure-function relations of distributed neural networks. However, due to differences in the biophysical origins of their signals and in the tissues to which they apply, there has been no direct integration of these techniques to date. We demonstrate that MRI signal variations and power spectra in a resting state are largely comparable between gray matter and white matter, that there are temporal correlations of fMRI signals that persist over long distances within distinct white matter structures, and that neighboring intervoxel correlations of low frequency resting state signals showed distinct anisotropy in many regions. These observations suggest that MRI signal variations from within white matter in a resting state may convey similar information as their corresponding fluctuations of MRI signals in gray matter. We thus derive a local spatio-temporal correlation tensor which captures directional variations of resting-state correlations and which reveals distinct structures in both white and gray matter. This novel concept is illustrated with in vivo experiments in a resting state, which demonstrate the potential of the technique for mapping the functional structure of neural networks and for direct integration of structure-function relations in the human brain.     

The effect of negative feedback loops on the dynamics of boolean networks

Feedback loops play an important role in determining the dynamics of biological networks. To study the role of negative feedback loops, this article introduces the notion of distance-to-positive-feedback which, in essence, captures the number of independent negative feedback loops in the network, a property inherent in the network topology. Through a computational study using Boolean networks, it is shown that distance-to-positive-feedback has a strong influence on network dynamics and correlates very well with the number and length of limit cycles in the phase space of the network. To be precise, it is shown that, as the number of independent negative feedback loops increases, the number (length) of limit cycles tends to decrease (increase). These conclusions are consistent with the fact that certain natural biological networks exhibit generally regular behavior and have fewer negative feedback loops than randomized networks with the same number of nodes and same connectivity.     

Modeling brain dynamics using computational neurogenetic approach

The paper introduces a novel computational approach to brain dynamics modeling that integrates dynamic gene–protein regulatory networks with a neural network model. Interaction of genes and proteins in neurons affects the dynamics of the whole neural network. Through tuning the gene–protein interaction network and the initial gene/protein expression values, different states of the neural network dynamics can be achieved. A generic computational neurogenetic model is introduced that implements this approach. It is illustrated by means of a simple neurogenetic model of a spiking neural network of the generation of local field potential. Our approach allows for investigation of how deleted or mutated genes can alter the dynamics of a model neural network. We conclude with the proposal how to extend this approach to model cognitive neurodynamics.

Motor Imagery Learning Modulates Functional Connectivity of Multiple Brain Systems in Resting State

Background

Learning motor skills involves subsequent modulation of resting-state functional connectivity in the sensory-motor system. This idea was mostly derived from the investigations on motor execution learning which mainly recruits the processing of sensory-motor information. Behavioral evidences demonstrated that motor skills in our daily lives could be learned through imagery procedures. However, it remains unclear whether the modulation of resting-state functional connectivity also exists in the sensory-motor system after motor imagery learning.

Methodology/Principal Findings

We performed a fMRI investigation on motor imagery learning from resting state. Based on previous studies, we identified eight sensory and cognitive resting-state networks (RSNs) corresponding to the brain systems and further explored the functional connectivity of these RSNs through the assessments, connectivity and network strengths before and after the two-week consecutive learning. Two intriguing results were revealed: (1) The sensory RSNs, specifically sensory-motor and lateral visual networks exhibited greater connectivity strengths in precuneus and fusiform gyrus after learning; (2) Decreased network strength induced by learning was proved in the default mode network, a cognitive RSN.

Conclusions/Significance

These results indicated that resting-state functional connectivity could be modulated by motor imagery learning in multiple brain systems, and such modulation displayed in the sensory-motor, visual and default brain systems may be associated with the establishment of motor schema and the regulation of introspective thought. These findings further revealed the neural substrates underlying motor skill learning and potentially provided new insights into the therapeutic benefits of motor imagery learning.