Behavior modulates effective connectivity between cortex and striatum

Dynamics Learning Rate Bias in Pigeons: Insights from Reinforcement Learning and Neural Correlates

Research in reinforcement learning indicates that animals respond differently to positive and negative reward prediction errors, which can be calculated by assuming learning rate bias. Many studies have shown that humans and other animals have learning rate bias during learning, but it is unclear whether and how the bias changes throughout the entire learning process. Here, we recorded the behavior data and the local field potentials (LFPs) in the striatum of five pigeons performing a probabilistic learning task. Reinforcement learning models with and without learning rate biases were used to dynamically fit the pigeons' choice behavior and estimate the option values. Furthemore, the correlation between the striatal LFPs power and the model-estimated option values was explored. We found that the pigeons' learning rate bias shifted from negative to positive during the learning process, and the striatal Gamma (31 to 80 Hz) power correlated with the option values modulated by dynamic learning rate bias. In conclusion, our results support the hypothesis that pigeons employ a dynamic learning strategy in the learning process from both behavioral and neural aspects, providing valuable insights into reinforcement learning mechanisms of non-human animals.

Altered Effective Connectivity Measured by Resting-State Functional Magnetic Resonance Imaging in Posterior Parietal-Frontal-Striatum Circuit in Patients With Disorder of Consciousness

Objective

Disorder of consciousness (DoC) resulting from severe brain injury is characterized by cortical and subcortical dysconnectivity. However, research on seed-based effective connectivity (EC) of DoC might be questioned as to the heterogeneity of prior assumptions.

Methods

Functional MRI data of 16 DoC patients and 16 demographically matched healthy individuals were analyzed. Revised coma recovery scale (CRS-R) scores of patients were acquired. Seed-based d mapping permutation of subject images (SDM-PSI) of meta-analysis was performed to quantitatively synthesize results from neuroimaging studies that evaluated resting-state functional activity in DoC patients. Spectral dynamic causal modeling (spDCM) was used to assess how EC altered between brain regions in DoC patients compared to healthy individuals.

Results

We found increased effective connectivity in left striatum and decreased effective connectivity in bilateral precuneus (preCUN)/posterior cingulate cortex (PCC), bilateral midcingulate cortex and left middle frontal gyrus in DoC compared with the healthy controls. The resulting pattern of interaction in DoC indicated disrupted connection and disturbance of posterior parietal-frontal-striatum, and reduced self-inhibition of preCUN/PCC. The strength of self-inhibition of preCUN/PCC was negatively correlated with the total score of CRS-R.

Conclusion

This impaired EC in DoC may underlie disruption in the posterior parietal-frontal-striatum circuit, particularly damage to the cortico-striatal connection and possible loss of preCUN/PCC function as the main regulatory hub.

Mapping functional gradients of the striatal circuit using simultaneous microelectric stimulation and ultrahigh-field fMRI in non-human primates

Advances in functional magnetic resonance imaging (fMRI) have significantly enhanced our understanding of the striatal system of both humans and non-human primates (NHP) over the last few decades. However, its circuit-level functional anatomy remains poorly understood, partly because in-vivo fMRI cannot directly perturb a brain system and map its casual input-output relationship. Also, routine 3T fMRI has an insufficient spatial resolution. We performed electrical microstimulation (EM) of the striatum in lightly-anesthetized NHPs while simultaneously mapping whole-brain activation, using contrast-enhanced fMRI at ultra-high-field 7T. By stimulating multiple positions along the striatum's main (dorsal-to-ventral) axis, we revealed its complex functional circuit concerning mutually connected subsystems in both cortical and subcortical areas. Indeed, within the striatum, there were distinct brain activation patterns across different stimulation sites. Specifically, dorsal stimulation revealed a medial-to-lateral elongated shape of activation in upper caudate and putamen areas, whereas ventral stimulation evoked areas confined to the medial and lower caudate. Such dorsoventral gradients also appeared in neocortical and thalamic activations, indicating consistent embedding profiles of the striatal system across the whole brain. These findings reflect different forms of within-circuit and inter-regional neuronal connectivity between the dorsal and ventromedial striatum. These patterns both shared and contrasted with previous anatomical tract-tracing and in-vivo resting-state fMRI studies. Our approach of combining microstimulation and whole-brain fMRI mapping in NHPs provides a unique opportunity to integrate our understanding of a targeted brain area's meso- and macro-scale functional systems.

Impairment in the goal-directed corticostriatal learning system as a biomarker for obsessive-compulsive disorder

BACKGROUND

Compulsive behaviors in obsessive-compulsive disorder (OCD) have been related to impairment within the associative cortical-striatal system connecting the caudate and prefrontal cortex that underlies consciously-controlled goal-directed learning and behavior. However, little is known whether this impairment may serve as a biomarker for vulnerability to OCD.

METHODS

Using resting-state functional magnetic resonance imaging (fMRI), we employed Granger causality analysis (GCA) to measure effective connectivity (EC) in previously validated striatal sub-regions, including the caudate, putamen, and the nucleus accumbens, in 35 OCD patients, 35 unaffected first-degree relatives and 35 matched healthy controls.

RESULTS

Both OCD patients and their first-degree relatives showed greater EC than controls between the left caudate and the orbital frontal cortex (OFC). Both OCD patients and their first-degree relatives showed lower EC than controls between the left caudate and lateral prefrontal cortex. These results are consistent with findings from task-related fMRI studies which found impairment in the goal-directed system in OCD patients.

CONCLUSIONS

The same changes in EC were present in both OCD patients and their unaffected first-degree relatives suggest that impairment in the goal-directed learning system may be a biomarker for OCD.

Lack of mutant huntingtin in cortical efferents improves behavioral inflexibility and corticostriatal dynamics in Huntington's disease mice

Abnormal communication between cerebral cortex and striatum plays a major role in the motor symptoms of Huntington's disease (HD), a neurodegenerative disorder caused by a mutation of the huntingtin gene (mHTT). Because cortex is the main driver of striatal processing, we recorded local field potential (LFP) activity simultaneously in primary motor cortex (M1) and dorsal striatum (DS) in BACHD mice, a full-length HD gene model, and in a conditional BACHD/Emx-1 Cre (BE) model in which mHTT is suppressed in cortical efferents, while mice freely explored a plus-shaped maze beginning at 20 wk of age. Relative to wild-type (WT) controls, BACHD mice were just as active across >40 wk of testing but became progressively less likely to turn into a perpendicular arm as they approached the choice point of the maze, a sign of HD motor inflexibility. BE mice, in contrast, turned as freely as WT throughout testing. Although BE mice did not exactly match WT in LFP activity, the reduction in alpha (8-13 Hz), beta (13-30 Hz), and low-gamma (30-50 Hz) power that occurred in M1 of turning-impaired BACHD mice was reversed. No reversal occurred in DS. In fact, BE mice showed further reductions in DS theta (4-8 Hz), beta, and low-gamma power relative to the BACHD model. Coherence analysis indicated a dysregulation of corticostriatal information flow in both BACHD and BE mice. Collectively, our results suggest that mHTT in cortical outputs drives the dysregulation of select cortical frequencies that accompany the loss of behavioral flexibility in HD.NEW & NOTEWORTHY BACHD mice, a full-length genetic model of Huntington's disease (HD), express aberrant local field potential (LFP) activity in primary motor cortex (M1) along with decreased probability of turning into a perpendicular arm of a plus-shaped maze, a motor inflexibility phenotype. Suppression of the mutant huntingtin gene in cortical output neurons prevents decline in turning and improves alpha, beta, and low-gamma activity in M1. Our results implicate cortical networks in the search for therapeutic strategies to alleviate HD motor signs.

Translating striatal activity from brain slice to whole animal neurophysiology: A guide for neuroscience research integrating diverse levels of analysis

An important goal of this review is highlighting research in neuroscience as examples of multilevel functional and anatomical analyses addressing basic science issues and applying results to the understanding of diverse disorders. The research of Dr. Michael Levine, a leader in neuroscience, exemplifies this approach by uncovering fundamental properties of basal ganglia function and translating these findings to clinical applications. The review focuses on neurophysiological research connecting results from in vitro and in vivo recordings. A second goal is to utilize these research connections to produce novel, accurate descriptions for corticostriatal processing involved in varied, complex functions. Medium spiny neurons in striatum act as integrators combining input with baseline activity creating motivational "events." Basic research on corticostriatal synapses is described and links developed to issues with clinical relevance such as inhibitory gating, self-injurious behavior, and relative reward valuation. Work is highlighted on dopamine-glutamate interactions. Individual medium spiny neurons express both D₁ and D₂ receptors and encode information in a bivalent manner depending upon the mix of receptors involved. Current work on neurophysiology of reward processing has taken advantage of these basic approaches at the cellular and molecular levels. Future directions in studying physiology of reward processing and action sequencing could profit by incorporating the divergent ways dopamine modulates incoming neurochemical signals. Primary investigators leading research teams should mirror Mike Levine's efforts in "climbing the mountain" of scientific inquiry by performing analyses at different levels of inquiry, integrating the findings, and building comprehensive answers to problems unsolvable without this bold approach.

A Tutorial for Information Theory in Neuroscience

Understanding how neural systems integrate, encode, and compute information is central to understanding brain function. Frequently, data from neuroscience experiments are multivariate, the interactions between the variables are nonlinear, and the landscape of hypothesized or possible interactions between variables is extremely broad. Information theory is well suited to address these types of data, as it possesses multivariate analysis tools, it can be applied to many different types of data, it can capture nonlinear interactions, and it does not require assumptions about the structure of the underlying data (i.e., it is model independent). In this article, we walk through the mathematics of information theory along with common logistical problems associated with data type, data binning, data quantity requirements, bias, and significance testing. Next, we analyze models inspired by canonical neuroscience experiments to improve understanding and demonstrate the strengths of information theory analyses. To facilitate the use of information theory analyses, and an understanding of how these analyses are implemented, we also provide a free MATLAB software package that can be applied to a wide range of data from neuroscience experiments, as well as from other fields of study.

Corticostriatal network dysfunction in Huntington's disease: Deficits in neural processing, glutamate transport, and ascorbate release

AIMS

This review summarizes evidence for dysfunctional connectivity between cortical and striatal neurons in Huntington's disease (HD), a fatal neurodegenerative condition caused by a single gene mutation. The focus is on data derived from recording of electrophysiological signals in behaving transgenic mouse models.

DISCUSSIONS

Firing patterns of individual neurons and the frequency oscillations of local field potentials indicate a disruption in corticostriatal processing driven, in large part, by interactions between cells that contain the mutant gene rather than the mutant gene alone. Dysregulation of glutamate, an excitatory amino acid released by cortical afferents, plays a key role in the breakdown of corticostriatal communication, a process modulated by ascorbate, an antioxidant vitamin found in high concentration in striatum. Up-regulation of glutamate transport by drug administration or viral-vector delivery improves ascorbate homeostasis and neurobehavioral processing in HD mice. Further analysis of electrophysiological data, including the use of sophisticated computational strategies, is required to discern how behavioral demands modulate the flow of corticostriatal information and its disruption by HD.

CONCLUSIONS

Long before massive cell loss occurs, HD impairs the mechanisms by which cortical and striatal neurons communicate. A key problem identified in transgenic animal models is dysregulation of the dynamic changes in extracellular glutamate and ascorbate. Improved understanding of how these neurochemical systems impact corticostriatal communication is necessary before an effective therapeutic strategy can emerge.

Efficient communication dynamics on macro-connectome, and the propagation speed

Global communication dynamics in the brain can be captured using fMRI, MEG, or electrocorticography (ECoG), and the global slow dynamics often represent anatomical constraints. Complementary single-/multi-unit recordings have described local fast temporal dynamics. However, global fast temporal dynamics remain incompletely understood with considering of anatomical constraints. Therefore, we compared temporal aspects of cross-area propagations of single-unit recordings and ECoG, and investigated their anatomical bases. First, we demonstrated how both evoked and spontaneous ECoGs can accurately predict latencies of single-unit recordings. Next, we estimated the propagation velocity (1.0–1.5 m/s) from brain-wide data and found that it was fairly stable among different conscious levels. We also found that the shortest paths in anatomical topology strongly predicted the latencies. Finally, we demonstrated that Communicability, a novel graph-theoretic measure, is able to quantify that more than 90% of paths should use shortest paths and the remaining are non-shortest walks. These results revealed that macro-connectome is efficiently wired for detailed communication dynamics in the brain.

Functional network stability and average minimal distance - A framework to rapidly assess dynamics of functional network representations

BACKGROUND

Recent advances in neurophysiological recording techniques have increased both the spatial and temporal resolution of data. New methodologies are required that can handle large data sets in an efficient manner as well as to make quantifiable, and realistic, predictions about the global modality of the brain from under-sampled recordings.

NEW METHOD

To rectify both problems, we first propose an analytical modification to an existing functional connectivity algorithm, Average Minimal Distance (AMD), to rapidly capture functional network connectivity. We then complement this algorithm by introducing Functional Network Stability (FuNS), a metric that can be used to quickly assess the global network dynamic changes over time, without being constrained by the activities of a specific set of neurons.

RESULTS

We systematically test the performance of AMD and FuNS (1) on artificial spiking data with different statistical characteristics, (2) from spiking data generated using a neural network model, and (3) using in vivo data recorded from mouse hippocampus during fear learning. Our results show that AMD and FuNS are able to monitor the change in network dynamics during memory consolidation.

COMPARISON WITH OTHER METHODS

AMD outperforms traditional bootstrapping and cross-correlation (CC) methods in both significance and computation time. Simultaneously, FuNS provides a reliable way to establish a link between local structural network changes, global dynamics of network-wide representations activity, and behavior.

CONCLUSIONS

The AMD-FuNS framework should be universally useful in linking long time-scale, global network dynamics and cognitive behavior.

Local or Not Local: Investigating the Nature of Striatal Theta Oscillations in Behaving Rats

In the cortex and hippocampus, neuronal oscillations of different frequencies can be observed in local field potentials (LFPs). LFPs oscillations in the theta band (6–10 Hz) have also been observed in the dorsolateral striatum (DLS) of rodents, mostly during locomotion, and have been proposed to mediate behaviorally-relevant interactions between striatum and cortex (or between striatum and hippocampus). However, it is unclear if these theta oscillations are generated in the striatum. To address this issue, we recorded LFPs and spiking activity in the DLS of rats performing a running sequence on a motorized treadmill. We observed an increase in rhythmical activity of the LFP in the theta-band during run compared to rest periods. However, several observations suggest that these oscillations are mainly generated outside of the striatum. First, theta oscillations disappeared when LFPs were rereferenced against a striatal recording electrode and the imaginary coherence between LFPs recorded at different locations within the striatum was null. Second, 8% of the recorded neurons had their spiking activity phase-locked to the theta rhythm. Third, Granger causality analyses between LFPs simultaneously recorded in the cortex and the striatum revealed that the interdependence between these two signals in the theta range was mostly accounted for by a common external source. The most parsimonious interpretation of these results is that theta oscillations observed in striatal LFPs are largely contaminated by volume-conducted signals. We propose that striatal LFPs are not optimal proxies of network dynamics in the striatum and should be interpreted with caution.

Alterations in Functional Cortical Hierarchy in Hemiparkinsonian Rats

Parkinson's disease and experimentally induced hemiparkinsonism are characterized by increased beta synchronization between cortical and subcortical areas. This change in beta connectivity might reflect either a symmetric increase in interareal influences or asymmetric changes in directed influences among brain areas. We assessed patterns of functional and directed connectivity within and between striatum and six cortical sites in each hemisphere of the hemiparkinsonian rat model. LFPs were recorded in resting and walking states, before and after unilateral 6-hydroxydopamine lesion. The hemiparkinsonian state was characterized by increased oscillatory activity in the 20-40 Hz range in resting and walking states, and increased interhemispheric coupling (phase lag index) that was more widespread at rest than during walking. Spectral Granger-causality analysis revealed that the change in symmetric functional connectivity comprised profound reorganization of hierarchical organization and directed influence patterns. First, in the lesioned hemisphere, the more anterior, nonprimary motor areas located at the top of the cortical hierarchy (i.e., receiving many directed influences) tended to increase their directed influence onto the posterior primary motor and somatosensory areas. This enhanced influence of "higher" areas may be related to the loss of motor control due to the 6-OHDA lesion. Second, the drive from the nonlesioned toward the lesioned hemisphere (in particular to striatum) increased, most prominently during walking. The nature of these adaptations (disturbed signaling or compensation) is discussed. The present study demonstrates that hemiparkinsonism is associated with a profound reorganization of the hierarchical organization of directed influence patterns among brain areas, perhaps reflecting compensatory processes.SIGNIFICANCE STATEMENT Parkinson's disease classically first becomes manifest in one hemibody before affecting both sides, suggesting that degeneration is asymmetrical. Our results suggest that asymmetrical degeneration of the dopaminergic system induces an increased drive from the nonlesioned toward the lesioned hemisphere and a profound reorganization of functional cortical hierarchical organization, leading to a stronger directed influence of hierarchically higher placed cortical areas over primary motor and somatosensory cortices. These changes may represent a compensatory mechanism for loss of motor control as a consequence of dopamine depletion.

Dysregulated corticostriatal activity in open-field behavior and the head-twitch response induced by the hallucinogen 2,5-dimethoxy-4-iodoamphetamine

The substituted amphetamine, 2,5-dimethoxy-4-iodoamphetamine (DOI), is a hallucinogen that has been used to model a variety of psychiatric conditions. Here, we studied the effect of DOI on neural activity recorded simultaneously in the primary motor cortex (M1) and dorsal striatum of freely behaving FvB/N mice. DOI significantly decreased the firing rate of individually isolated neurons in M1 and dorsal striatum relative to pre-drug baseline. It also induced a bursting pattern of activity by increasing both the number of spikes within a burst and burst duration. In addition, DOI increased coincident firing between simultaneously recorded neuron pairs within the striatum and between M1 and dorsal striatum. Local field potential (LFP) activity also increased in coherence between M1 and dorsal striatum after DOI in the low frequency gamma band (30-50 Hz), while corticostriatal coherence in delta, theta, alpha, and beta activity decreased. We also assessed corticostriatal LFP activity in relation to the DOI-induced head-twitch response (HTR), a readily identifiable behavior used to assess potential treatments for the conditions it models. The HTR was associated with increased delta and decreased theta power in both M1 and dorsal striatum. Together, our results suggest that DOI dysregulates corticostriatal communication and that the HTR is associated with this dysregulation.

A novel cross-frequency coupling detection method using the generalized Morse wavelets

BACKGROUND

Cross-frequency coupling (CFC) occurs when non-identical frequency components entrain one another. A ubiquitous example from neuroscience is low frequency phase to high frequency amplitude coupling in electrophysiological signals. Seminal work by Canolty revealed CFC in human ECoG data. Established methods band-pass the data into component frequencies then convert the band-passed signals into the analytic representation, from which we infer the instantaneous amplitude and phase of each component. Though powerful, such methods resolve signals with respect to time and frequency without addressing the multiresolution problem.

NEW METHOD

We build upon the ground-breaking work of Canolty and others and derive a wavelet-based CFC detection algorithm that efficiently searches a range of frequencies using a sequence of filters with optimal trade-off between time and frequency resolution. We validate our method using simulated data and analyze CFC within and between the primary motor cortex and dorsal striatum of rats under ketamine-xylazine anesthesia.

RESULTS

Our method detects the correct CFC in simulated data and reveals CFC between frequency bands that were previously shown to participate in corticostriatal effective connectivity.

COMPARISON WITH EXISTING METHODS

Other CFC detection methods address the need to increase bandwidth when analyzing high frequency components but none to date permit rigorous bandwidth selection with no a priori knowledge of underlying CFC. Our method is thus particularly useful for exploratory studies.

CONCLUSIONS

The method developed here permits rigorous and efficient exploration of a hypothesis space and is particularly useful when the frequencies participating in CFC are unknown.

High-Degree Neurons Feed Cortical Computations

Recent work has shown that functional connectivity among cortical neurons is highly varied, with a small percentage of neurons having many more connections than others. Also, recent theoretical developments now make it possible to quantify how neurons modify information from the connections they receive. Therefore, it is now possible to investigate how information modification, or computation, depends on the number of connections a neuron receives (in-degree) or sends out (out-degree). To do this, we recorded the simultaneous spiking activity of hundreds of neurons in cortico-hippocampal slice cultures using a high-density 512-electrode array. This preparation and recording method combination produced large numbers of neurons recorded at temporal and spatial resolutions that are not currently available in any in vivo recording system. We utilized transfer entropy (a well-established method for detecting linear and nonlinear interactions in time series) and the partial information decomposition (a powerful, recently developed tool for dissecting multivariate information processing into distinct parts) to quantify computation between neurons where information flows converged. We found that computations did not occur equally in all neurons throughout the networks. Surprisingly, neurons that computed large amounts of information tended to receive connections from high out-degree neurons. However, the in-degree of a neuron was not related to the amount of information it computed. To gain insight into these findings, we developed a simple feedforward network model. We found that a degree-modified Hebbian wiring rule best reproduced the pattern of computation and degree correlation results seen in the real data. Interestingly, this rule also maximized signal propagation in the presence of network-wide correlations, suggesting a mechanism by which cortex could deal with common random background input. These are the first results to show that the extent to which a neuron modifies incoming information streams depends on its topological location in the surrounding functional network.

We recorded the electrical activity of hundreds of neurons simultaneously in brain tissue from mice and we analyzed these signals using state-of-the-art tools from information theory. These tools allowed us to ascertain which neurons were transmitting information to other neurons and to characterize the computations performed by neurons using the inputs they received from two or more other neurons. We found that computations did not occur equally in all neurons throughout the networks. Surprisingly, neurons that computed large amounts of information tended to be recipients of information from neurons with a large number of outgoing connections. Interestingly, the number of incoming connections to a neuron was not related to the amount of information that neuron computed. To better understand these results, we built a network model to match the data. Unexpectedly, the model also maximized information transfer in the presence of network-wide correlations. This suggested a way that networks of cortical neurons could deal with common random background input. These results are the first to show that the amount of information computed by a neuron depends on where it is located in the surrounding network.

Untangling Cortico-Striatal Connectivity and Cross-Frequency Coupling in L-DOPA-Induced Dyskinesia

We simultaneously recorded local field potentials (LFPs) in the primary motor cortex and sensorimotor striatum in awake, freely behaving, 6-OHDA lesioned hemi-parkinsonian rats in order to study the features directly related to pathological states such as parkinsonian state and levodopa-induced dyskinesia. We analyzed the spectral characteristics of the obtained signals and observed that during dyskinesia the most prominent feature was a relative power increase in the high gamma frequency range at around 80 Hz, while for the parkinsonian state it was in the beta frequency range. Here we show that during both pathological states effective connectivity in terms of Granger causality is bidirectional with an accent on the striatal influence on the cortex. In the case of dyskinesia, we also found a high increase in effective connectivity at 80 Hz. In order to further understand the 80-Hz phenomenon, we performed cross-frequency analysis and observed characteristic patterns in the case of dyskinesia but not in the case of the parkinsonian state or the control state. We noted a large decrease in the modulation of the amplitude at 80 Hz by the phase of low frequency oscillations (up to ~10 Hz) across both structures in the case of dyskinesia. This may suggest a lack of coupling between the low frequency activity of the recorded network and the group of neurons active at ~80 Hz.

Rich-Club Organization in Effective Connectivity among Cortical Neurons

The performance of complex networks, like the brain, depends on how effectively their elements communicate. Despite the importance of communication, it is virtually unknown how information is transferred in local cortical networks, consisting of hundreds of closely spaced neurons. To address this, it is important to record simultaneously from hundreds of neurons at a spacing that matches typical axonal connection distances, and at a temporal resolution that matches synaptic delays. We used a 512-electrode array (60 μm spacing) to record spontaneous activity at 20 kHz from up to 500 neurons simultaneously in slice cultures of mouse somatosensory cortex for 1 h at a time. We applied a previously validated version of transfer entropy to quantify information transfer. Similar to in vivo reports, we found an approximately lognormal distribution of firing rates. Pairwise information transfer strengths also were nearly lognormally distributed, similar to reports of synaptic strengths. Some neurons transferred and received much more information than others, which is consistent with previous predictions. Neurons with the highest outgoing and incoming information transfer were more strongly connected to each other than chance, thus forming a "rich club." We found similar results in networks recorded in vivo from rodent cortex, suggesting the generality of these findings. A rich-club structure has been found previously in large-scale human brain networks and is thought to facilitate communication between cortical regions. The discovery of a small, but information-rich, subset of neurons within cortical regions suggests that this population will play a vital role in communication, learning, and memory. Significance statement: Many studies have focused on communication networks between cortical brain regions. In contrast, very few studies have examined communication networks within a cortical region. This is the first study to combine such a large number of neurons (several hundred at a time) with such high temporal resolution (so we can know the direction of communication between neurons) for mapping networks within cortex. We found that information was not transferred equally through all neurons. Instead, ∼70% of the information passed through only 20% of the neurons. Network models suggest that this highly concentrated pattern of information transfer would be both efficient and robust to damage. Therefore, this work may help in understanding how the cortex processes information and responds to neurodegenerative diseases.

Functional Clusters, Hubs, and Communities in the Cortical Microconnectome

Although relationships between networks of different scales have been observed in macroscopic brain studies, relationships between structures of different scales in networks of neurons are unknown. To address this, we recorded from up to 500 neurons simultaneously from slice cultures of rodent somatosensory cortex. We then measured directed effective networks with transfer entropy, previously validated in simulated cortical networks. These effective networks enabled us to evaluate distinctive nonrandom structures of connectivity at 2 different scales. We have 4 main findings. First, at the scale of 3-6 neurons (clusters), we found that high numbers of connections occurred significantly more often than expected by chance. Second, the distribution of the number of connections per neuron (degree distribution) had a long tail, indicating that the network contained distinctively high-degree neurons, or hubs. Third, at the scale of tens to hundreds of neurons, we typically found 2-3 significantly large communities. Finally, we demonstrated that communities were relatively more robust than clusters against shuffling of connections. We conclude the microconnectome of the cortex has specific organization at different scales, as revealed by differences in robustness. We suggest that this information will help us to understand how the microconnectome is robust against damage.

Brain activity mapping at multiple scales with silicon microprobes containing 1,024 electrodes

The coordinated activity of neural ensembles across multiple interconnected regions has been challenging to study in the mammalian brain with cellular resolution using conventional recording tools. For instance, neural systems regulating learned behaviors often encompass multiple distinct structures that span the brain. To address this challenge we developed a three-dimensional (3D) silicon microprobe capable of simultaneously measuring extracellular spike and local field potential activity from 1,024 electrodes. The microprobe geometry can be precisely configured during assembly to target virtually any combination of four spatially distinct neuroanatomical planes. Here we report on the operation of such a device built for high-throughput monitoring of neural signals in the orbitofrontal cortex and several nuclei in the basal ganglia. We perform analysis on systems-level dynamics and correlations during periods of conditioned behavioral responding and rest, demonstrating the technology's ability to reveal functional organization at multiple scales in parallel in the mouse brain.

Human subthalamic nucleus in movement error detection and its evaluation during visuomotor adaptation

Monitoring and evaluating movement errors to guide subsequent movements is a critical feature of normal motor control. Previously, we showed that the postmovement increase in electroencephalographic (EEG) beta power over the sensorimotor cortex reflects neural processes that evaluate motor errors consistent with Bayesian inference (Tan et al., 2014). Whether such neural processes are limited to this cortical region or involve the basal ganglia is unclear. Here, we recorded EEG over the cortex and local field potential (LFP) activity in the subthalamic nucleus (STN) from electrodes implanted in patients with Parkinson's disease, while they moved a joystick-controlled cursor to visual targets displayed on a computer screen. After movement offsets, we found increased beta activity in both local STN LFP and sensorimotor cortical EEG and in the coupling between the two, which was affected by both error magnitude and its contextual saliency. The postmovement increase in the coupling between STN and cortex was dominated by information flow from sensorimotor cortex to STN. However, an information drive appeared from STN to sensorimotor cortex in the first phase of the adaptation, when a constant rotation was applied between joystick inputs and cursor outputs. The strength of the STN to cortex drive correlated with the degree of adaption achieved across subjects. These results suggest that oscillatory activity in the beta band may dynamically couple the sensorimotor cortex and basal ganglia after movements. In particular, beta activity driven from the STN to cortex indicates task-relevant movement errors, information that may be important in modifying subsequent motor responses.

Multiplex networks of cortical and hippocampal neurons revealed at different timescales

Recent studies have emphasized the importance of multiplex networks--interdependent networks with shared nodes and different types of connections--in systems primarily outside of neuroscience. Though the multiplex properties of networks are frequently not considered, most networks are actually multiplex networks and the multiplex specific features of networks can greatly affect network behavior (e.g. fault tolerance). Thus, the study of networks of neurons could potentially be greatly enhanced using a multiplex perspective. Given the wide range of temporally dependent rhythms and phenomena present in neural systems, we chose to examine multiplex networks of individual neurons with time scale dependent connections. To study these networks, we used transfer entropy--an information theoretic quantity that can be used to measure linear and nonlinear interactions--to systematically measure the connectivity between individual neurons at different time scales in cortical and hippocampal slice cultures. We recorded the spiking activity of almost 12,000 neurons across 60 tissue samples using a 512-electrode array with 60 micrometer inter-electrode spacing and 50 microsecond temporal resolution. To the best of our knowledge, this preparation and recording method represents a superior combination of number of recorded neurons and temporal and spatial recording resolutions to any currently available in vivo system. We found that highly connected neurons ("hubs") were localized to certain time scales, which, we hypothesize, increases the fault tolerance of the network. Conversely, a large proportion of non-hub neurons were not localized to certain time scales. In addition, we found that long and short time scale connectivity was uncorrelated. Finally, we found that long time scale networks were significantly less modular and more disassortative than short time scale networks in both tissue types. As far as we are aware, this analysis represents the first systematic study of temporally dependent multiplex networks among individual neurons.