Modeling fluctuations in default-mode brain network using a spiking neural network

Recently, numerous attempts have been made to understand the dynamic behavior of complex brain systems using neural network models. The fluctuations in blood-oxygen-level-dependent (BOLD) brain signals at less than 0.1 Hz have been observed by functional magnetic resonance imaging (fMRI) for subjects in a resting state. This phenomenon is referred to as a "default-mode brain network." In this study, we model the default-mode brain network by functionally connecting neural communities composed of spiking neurons in a complex network. Through computational simulations of the model, including transmission delays and complex connectivity, the network dynamics of the neural system and its behavior are discussed. The results show that the power spectrum of the modeled fluctuations in the neuron firing patterns is consistent with the default-mode brain network's BOLD signals when transmission delays, a characteristic property of the brain, have finite values in a given range.     

The brainstem reticular formation is a small-world, not scale-free, network

Recently, it has been demonstrated that several complex systems may have simple graph-theoretic characterizations as so-called 'small-world' and 'scale-free' networks. These networks have also been applied to the gross neural connectivity between primate cortical areas and the nervous system of Caenorhabditis elegans. Here, we extend this work to a specific neural circuit of the vertebrate brain--the medial reticular formation (RF) of the brainstem--and, in doing so, we have made three key contributions. First, this work constitutes the first model (and quantitative review) of this important brain structure for over three decades. Second, we have developed the first graph-theoretic analysis of vertebrate brain connectivity at the neural network level. Third, we propose simple metrics to quantitatively assess the extent to which the networks studied are small-world or scale-free. We conclude that the medial RF is configured to create small-world (implying coherent rapid-processing capabilities), but not scale-free, type networks under assumptions which are amenable to quantitative measurement.     

Brain networks of explicit and implicit learning

Are explicit versus implicit learning mechanisms reflected in the brain as distinct neural structures, as previous research indicates, or are they distinguished by brain networks that involve overlapping systems with differential connectivity? In this functional MRI study we examined the neural correlates of explicit and implicit learning of artificial grammar sequences. Using effective connectivity analyses we found that brain networks of different connectivity underlie the two types of learning: while both processes involve activation in a set of cortical and subcortical structures, explicit learners engage a network that uses the insula as a key mediator whereas implicit learners evoke a direct frontal-striatal network. Individual differences in working memory also differentially impact the two types of sequence learning.     

Integration of sleep homeostasis and navigation in Drosophila

During sleep, the brain undergoes dynamic and structural changes. In Drosophila, such changes have been observed in the central complex, a brain area important for sleep control and navigation. The connectivity of the central complex raises the question about how navigation, and specifically the head direction system, can operate in the face of sleep related plasticity. To address this question, we develop a model that integrates sleep homeostasis and head direction. We show that by introducing plasticity, the head direction system can function in a stable way by balancing plasticity in connected circuits that encode sleep pressure. With increasing sleep pressure, the head direction system nevertheless becomes unstable and a sleep phase with a different plasticity mechanism is introduced to reset network connectivity. The proposed integration of sleep homeostasis and head direction circuits captures features of their neural dynamics observed in flies and mice.     

A Neural Population Model Incorporating Dopaminergic Neurotransmission during Complex Voluntary Behaviors

Assessing brain activity during complex voluntary motor behaviors that require the recruitment of multiple neural sites is a field of active research. Our current knowledge is primarily based on human brain imaging studies that have clear limitations in terms of temporal and spatial resolution. We developed a physiologically informed non-linear multi-compartment stochastic neural model to simulate functional brain activity coupled with neurotransmitter release during complex voluntary behavior, such as speech production. Due to its state-dependent modulation of neural firing, dopaminergic neurotransmission plays a key role in the organization of functional brain circuits controlling speech and language and thus has been incorporated in our neural population model. A rigorous mathematical proof establishing existence and uniqueness of solutions to the proposed model as well as a computationally efficient strategy to numerically approximate these solutions are presented. Simulated brain activity during the resting state and sentence production was analyzed using functional network connectivity, and graph theoretical techniques were employed to highlight differences between the two conditions. We demonstrate that our model successfully reproduces characteristic changes seen in empirical data between the resting state and speech production, and dopaminergic neurotransmission evokes pronounced changes in modeled functional connectivity by acting on the underlying biological stochastic neural model. Specifically, model and data networks in both speech and rest conditions share task-specific network features: both the simulated and empirical functional connectivity networks show an increase in nodal influence and segregation in speech over the resting state. These commonalities confirm that dopamine is a key neuromodulator of the functional connectome of speech control. Based on reproducible characteristic aspects of empirical data, we suggest a number of extensions of the proposed methodology building upon the current model.     

Intrinsic Functional Brain Architecture Derived from Graph Theoretical Analysis in the Human Fetus

The human brain undergoes dramatic maturational changes during late stages of fetal and early postnatal life. The importance of this period to the establishment of healthy neural connectivity is apparent in the high incidence of neural injury in preterm infants, in whom untimely exposure to ex-uterine factors interrupts neural connectivity. Though the relevance of this period to human neuroscience is apparent, little is known about functional neural networks in human fetal life. Here, we apply graph theoretical analysis to examine human fetal brain connectivity. Utilizing resting state functional magnetic resonance imaging (fMRI) data from 33 healthy human fetuses, 19 to 39 weeks gestational age (GA), our analyses reveal that the human fetal brain has modular organization and modules overlap functional systems observed postnatally. Age-related differences between younger (GA <31 weeks) and older (GA≥31 weeks) fetuses demonstrate that brain modularity decreases, and connectivity of the posterior cingulate to other brain networks becomes more negative, with advancing GA. By mimicking functional principles observed postnatally, these results support early emerging capacity for information processing in the human fetal brain. Current technical limitations, as well as the potential for fetal fMRI to one day produce major discoveries about fetal origins or antecedents of neural injury or disease are discussed.     

Dynamical correlation patterns and corresponding community structure in neural spontaneous activity at criticality

It has been considered that the state in the vicinity of a critical point, which is the point between ordered and disordered states, can underlie and facilitate information processing of the brain in various aspects. In this research, we numerically study the influence of criticality on one aspect of brain information processing, i.e., the community structure, which is an important characteristic of complex networks. We examine community structure of the functional connectivity in simulated brain spontaneous activity, which is based on dynamical correlations between neural activity patterns at different positions. The brain spontaneous activity is simulated by a neural field model whose parameter covers subcritical, critical, and supercritical regions. Then, the corresponding dynamical correlation patterns and community structure are compared. In the critical region, we found some distinctive properties, namely high correlation and correlation switching, high modularity and a low number of modules, high stability of the dynamical functional connectivity, and moderate flexibility of the community structure across temporal scales. We also discuss how these characteristics might improve information processing of the brain.     

3D multicolor super-resolution imaging offers improved accuracy in neuron tracing

The connectivity among neurons holds the key to understanding brain function. Mapping neural connectivity in brain circuits requires imaging techniques with high spatial resolution to facilitate neuron tracing and high molecular specificity to mark different cellular and molecular populations. Here, we tested a three-dimensional (3D), multicolor super-resolution imaging method, stochastic optical reconstruction microscopy (STORM), for tracing neural connectivity using cultured hippocampal neurons obtained from wild-type neonatal rat embryos as a model system. Using a membrane specific labeling approach that improves labeling density compared to cytoplasmic labeling, we imaged neural processes at 44 nm 2D and 116 nm 3D resolution as determined by considering both the localization precision of the fluorescent probes and the Nyquist criterion based on label density. Comparison with confocal images showed that, with the currently achieved resolution, we could distinguish and trace substantially more neuronal processes in the super-resolution images. The accuracy of tracing was further improved by using multicolor super-resolution imaging. The resolution obtained here was largely limited by the label density and not by the localization precision of the fluorescent probes. Therefore, higher image resolution, and thus higher tracing accuracy, can in principle be achieved by further improving the label density.     

Resting state functional connectivity in patients with chronic hallucinations

Auditory verbal hallucinations (AVH) are not only among the most common but also one of the most distressing symptoms of schizophrenia. Despite elaborate research, the underlying brain mechanisms are as yet elusive. Functional MRI studies have associated the experience of AVH with activation of bilateral language-related areas, in particular the right inferior frontal gyrus (rIFG) and the left superior temporal gyrus (lSTG). While these findings helped to understand the neural underpinnings of hearing voices, they provide little information about possible brain mechanisms that predispose a person to experience AVH, i.e. the traits to hallucinate. In this study, we compared resting state connectivity between 49 psychotic patients with chronic AVH and 49 matched controls using the rIFG and the lSTG as seed regions, to identify functional brain systems underlying the predisposition to hallucinate. The right parahippocampal gyrus showed increased connectivity with the rIFG in patients as compared to controls. Reduced connectivity with the rIFG in patients was found for the right dorsolateral prefrontal cortex. Reduced connectivity with the lSTG in patients was identified in the left frontal operculum as well as the parietal opercular area. Connectivity between the lSTG and the left hippocampus was also reduced in patients and showed a negative correlation with the severity of hallucinations. Concluding, we found aberrant connectivity between the seed regions and medial temporal lobe structures which have a prominent role in memory retrieval. Moreover, we found decreased connectivity between language-related areas, indicating aberrant integration in this system potentially including corollary discharge mechanisms.     

Asymmetry in neural fields: a spatiotemporal encoding mechanism

Neural field models have been successfully applied to model diverse brain mechanisms like visual attention, motor control, and memory. Most theoretical and modeling works have focused on the study of the dynamics of such systems under variations in neural connectivity, mainly symmetric connectivity among neurons. However, less attention has been given to the emerging properties of neuron populations when neural connectivity is asymmetric, although asymmetric activity propagation has been observed in cortical tissue. Here we explore the dynamics of neural fields with asymmetric connectivity and show, in the case of front propagation, that it can bias the population to follow a certain trajectory with higher activation. We find that asymmetry relates linearly to the input speed when the input is spatially localized, and this relation holds for different kernels and input shapes. To illustrate the behavior of asymmetric connectivity, we present an application: standard video sequences of human motion were encoded using the asymmetric neural field and compared to computer vision techniques. Overall, our results indicate that asymmetric neural fields are a competitive approach for spatiotemporal encoding with two main advantages: online classification and distributed operation.     

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.     

Dynamic effective connectivity of inter-areal brain circuits

Anatomic connections between brain areas affect information flow between neuronal circuits and the synchronization of neuronal activity. However, such structural connectivity does not coincide with effective connectivity (or, more precisely, causal connectivity), related to the elusive question “Which areas cause the present activity of which others?”. Effective connectivity is directed and depends flexibly on contexts and tasks. Here we show that dynamic effective connectivity can emerge from transitions in the collective organization of coherent neural activity. Integrating simulation and semi-analytic approaches, we study mesoscale network motifs of interacting cortical areas, modeled as large random networks of spiking neurons or as simple rate units. Through a causal analysis of time-series of model neural activity, we show that different dynamical states generated by a same structural connectivity motif correspond to distinct effective connectivity motifs. Such effective motifs can display a dominant directionality, due to spontaneous symmetry breaking and effective entrainment between local brain rhythms, although all connections in the considered structural motifs are reciprocal. We show then that transitions between effective connectivity configurations (like, for instance, reversal in the direction of inter-areal interactions) can be triggered reliably by brief perturbation inputs, properly timed with respect to an ongoing local oscillation, without the need for plastic synaptic changes. Finally, we analyze how the information encoded in spiking patterns of a local neuronal population is propagated across a fixed structural connectivity motif, demonstrating that changes in the active effective connectivity regulate both the efficiency and the directionality of information transfer. Previous studies stressed the role played by coherent oscillations in establishing efficient communication between distant areas. Going beyond these early proposals, we advance here that dynamic interactions between brain rhythms provide as well the basis for the self-organized control of this “communication-through-coherence”, making thus possible a fast “on-demand” reconfiguration of global information routing modalities.     

How do neural connectivity and time delays influence bimanual coordination?

Multilevel crosstalk as a neural basis for motor control has been widely discussed in the literature. Since no natural process is instantaneous, any crosstalk model should incorporate time delays, which are known to induce temporal coupling between functional elements and stabilize or destabilize a particular mode of coordination. In this article, we systematically study the dynamics of rhythmic bimanual coordination under the influence of varying connection topology as realized by callosal fibers, cortico-thalamic projections, and crossing peripheral fibers. Such connectivity contributes to various degrees of neural crosstalk between the effectors which we continuously parameterize in a mathematical model. We identify the stability regimes of bimanual coordination as a function of the degree of neural crosstalk, movement amplitude and the time delays involved due to signal processing. Prominent examples include explanations of the decreased stability of the antiphase mode of coordination in split brain patients and the role of coupling in mediating bimanual coordination.     

Neurovesicles in Brain Development

Long before the nervous system is organized into electrically active neural circuits, connectivity emerges between cells of the developing brain through extracellular signals. Extracellular vesicles that shuttle RNA, proteins, and lipids from donor cells to recipient cells are candidates for mediating connectivity in the brain. Despite the abundance of extracellular vesicles during brain development, evidence for their physiological functions is only beginning to materialize. Here, we review evidence of the existence, content, and functions of extracellular vesicles in brain development.     

Neural pattern dynamics in an oscillator model of the thalamo-reticular system.

One of development issues for information processing with synchronous oscillations in the brain is how new information is coded and how a comparison with already existing information is performed. In the present work we study a simple neural network model of the thalamo-reticular system based on the Wilson-Cowan model of neuronal oscillatory behavior. Our results show that both cortical control over the thalamus and external sensory input are essential in coordinating and generating spatio-temporal patterns of synchronous activity. A main finding of the numerical simulations is that the network connectivity and the intrinsic oscillatory properties of the neurons result in distinct collective behaviors within the network. By varying the connectivity schemes comparable with lesionated or damaged brain regions our results are in good agreement with in vivo experimental results. Suppressing the sensory input results in temporal oscillatory activity in the beta and gamma range and a strong spatial dependence of the network activity.     

Gene expression in the rodent brain is associated with its regional connectivity

The putative link between gene expression of brain regions and their neural connectivity patterns is a fundamental question in neuroscience. Here this question is addressed in the first large scale study of a prototypical mammalian rodent brain, using a combination of rat brain regional connectivity data with gene expression of the mouse brain. Remarkably, even though this study uses data from two different rodent species (due to the data limitations), we still find that the connectivity of the majority of brain regions is highly predictable from their gene expression levels-the outgoing (incoming) connectivity is successfully predicted for 73% (56%) of brain regions, with an overall fairly marked accuracy level of 0.79 (0.83). Many genes are found to play a part in predicting both the incoming and outgoing connectivity (241 out of the 500 top selected genes, p-value<1e-5). Reassuringly, the genes previously known from the literature to be involved in axon guidance do carry significant information about regional brain connectivity. Surveying the genes known to be associated with the pathogenesis of several brain disorders, we find that those associated with schizophrenia, autism and attention deficit disorder are the most highly enriched in the connectivity-related genes identified here. Finally, we find that the profile of functional annotation groups that are associated with regional connectivity in the rodent is significantly correlated with the annotation profile of genes previously found to determine neural connectivity in C. elegans (Pearson correlation of 0.24, p<1e-6 for the outgoing connections and 0.27, p<1e-5 for the incoming). Overall, the association between connectivity and gene expression in a specific extant rodent species' brain is likely to be even stronger than found here, given the limitations of current data.     

Estimation of Directed Effective Connectivity from fMRI Functional Connectivity Hints at Asymmetries of Cortical Connectome

The brain exhibits complex spatio-temporal patterns of activity. This phenomenon is governed by an interplay between the internal neural dynamics of cortical areas and their connectivity. Uncovering this complex relationship has raised much interest, both for theory and the interpretation of experimental data (e.g., fMRI recordings) using dynamical models. Here we focus on the so-called inverse problem: the inference of network parameters in a cortical model to reproduce empirically observed activity. Although it has received a lot of interest, recovering directed connectivity for large networks has been rather unsuccessful so far. The present study specifically addresses this point for a noise-diffusion network model. We develop a Lyapunov optimization that iteratively tunes the network connectivity in order to reproduce second-order moments of the node activity, or functional connectivity. We show theoretically and numerically that the use of covariances with both zero and non-zero time shifts is the key to infer directed connectivity. The first main theoretical finding is that an accurate estimation of the underlying network connectivity requires that the time shift for covariances is matched with the time constant of the dynamical system. In addition to the network connectivity, we also adjust the intrinsic noise received by each network node. The framework is applied to experimental fMRI data recorded for subjects at rest. Diffusion-weighted MRI data provide an estimate of anatomical connections, which is incorporated to constrain the cortical model. The empirical covariance structure is reproduced faithfully, especially its temporal component (i.e., time-shifted covariances) in addition to the spatial component that is usually the focus of studies. We find that the cortical interactions, referred to as effective connectivity, in the tuned model are not reciprocal. In particular, hubs are either receptors or feeders: they do not exhibit both strong incoming and outgoing connections. Our results sets a quantitative ground to explore the propagation of activity in the cortex.     

Defining functional gene‐circuit interfaces in the mouse nervous system

Complexity in the nervous system is established by developmental genetic programs, maintained by differential genetic profiles and sculpted by experiential and environmental influence over gene expression. Determining how specific genes define neuronal phenotypes, shape circuit connectivity and regulate circuit function is essential for understanding how the brain processes information, directs behavior and adapts to changing environments. Mouse genetics has contributed greatly to current percepts of gene‐circuit interfaces in behavior, but considerable work remains. Large‐scale initiatives to map gene expression and connectivity in the brain, together with advanced techniques in molecular genetics, now allow detailed exploration of the genetic basis of nervous system function at the level of specific circuit connections. In this review, we highlight several key advances for defining the function of specific genes within a neural network .     

The Rich Get Richer: Brain Injury Elicits Hyperconnectivity in Core Subnetworks

There remains much unknown about how large-scale neural networks accommodate neurological disruption, such as moderate and severe traumatic brain injury (TBI). A primary goal in this study was to examine the alterations in network topology occurring during the first year of recovery following TBI. To do so we examined 21 individuals with moderate and severe TBI at 3 and 6 months after resolution of posttraumatic amnesia and 15 age- and education-matched healthy adults using functional MRI and graph theoretical analyses. There were two central hypotheses in this study: 1) physical disruption results in increased functional connectivity, or hyperconnectivity, and 2) hyperconnectivity occurs in regions typically observed to be the most highly connected cortical hubs, or the “rich club”. The current findings generally support the hyperconnectivity hypothesis showing that during the first year of recovery after TBI, neural networks show increased connectivity, and this change is disproportionately represented in brain regions belonging to the brain''s core subnetworks. The selective increases in connectivity observed here are consistent with the preferential attachment model underlying scale-free network development. This study is the largest of its kind and provides the unique opportunity to examine how neural systems adapt to significant neurological disruption during the first year after injury.     

Identification of neural activity based on fMRI data: A simulation study

This study proposes a technique for determining effective connectivity among brain regions which operates at the level of neuronal dynamics. We propose an alternative time-variant dynamic causal model (TV-DCM) where neuronal dynamic activity evolves based on generalized dynamic neural networks (GDNNs). The identification of brain architecture connectivity is carried out based on a least squares criterion and on a global search technique. Computer simulations carried out in the paper show that TV-DCM may provide multiple solutions, i.e. a set of different architectures all of which approximate the data equally well. Numerical comparisons between TV-DCM and DCM are also given. In order to determine the unique causal structure of brain regions, we apply an additional criterion, i.e. an estimation of generalization error, known from the theory of neural networks. Computer simulations also confirm the validity of our techniques.