Multiple Midfrontal Thetas Revealed by Source Separation of Simultaneous MEG and EEG.

Theta-band (∼6 Hz) rhythmic activity within and over the medial PFC ("midfrontal theta") has been identified as a distinctive signature of "response conflict," the competition between multiple actions when only one action is goal-relevant. Midfrontal theta is traditionally conceptualized and analyzed under the assumption that it is a unitary signature of conflict that can be uniquely identified at one electrode (typically FCz). Here we recorded simultaneous MEG and EEG (total of 328 sensors) in 9 human subjects (7 female) and applied a feature-guided multivariate source-separation decomposition to determine whether conflict-related midfrontal theta is a unitary or multidimensional feature of the data. For each subject, a generalized eigendecomposition yielded spatial filters (components) that maximized the ratio between theta and broadband activity. Components were retained based on significance thresholding and midfrontal EEG topography. All of the subjects individually exhibited multiple (mean 5.89, SD 2.47) midfrontal components that contributed to sensor-level midfrontal theta power during the task. Component signals were temporally uncorrelated and asynchronous, suggesting that each midfrontal theta component was unique. Our findings call into question the dominant notion that midfrontal theta represents a unitary process. Instead, we suggest that midfrontal theta spans a multidimensional space, indicating multiple origins, but can manifest as a single feature at the sensor level because of signal mixing.SIGNIFICANCE STATEMENT "Midfrontal theta" is a rhythmic electrophysiological signature of the competition between multiple response options. Midfrontal theta is traditionally considered to reflect a single process. However, this assumption could be erroneous because of "mixing" (multiple sources contributing to the activity recorded at a single electrode). We investigated the dimensionality of midfrontal theta by applying advanced multivariate analysis methods to a multimodal MEG/EEG dataset. We identified multiple topographically overlapping neural sources that drove response conflict-related midfrontal theta. Midfrontal theta thus reflects multiple uncorrelated signals that manifest with similar EEG scalp projections. In addition to contributing to the cognitive control literature, we demonstrate both the feasibility and the necessity of signal demixing to understand the narrowband neural dynamics underlying cognitive processes.

Flexible Frequency Switching in Adult Mouse Visual Cortex Is Mediated by Competition Between Parvalbumin and Somatostatin Expressing Interneurons.

Neuronal networks in rodent primary visual cortex (V1) can generate oscillations in different frequency bands depending on the network state and the level of visual stimulation. High-frequency gamma rhythms, for example, dominate the network's spontaneous activity in adult mice but are attenuated upon visual stimulation, during which the network switches to the beta band instead. The spontaneous local field potential (LFP) of juvenile mouse V1, however, mainly contains beta rhythms and presenting a stimulus does not elicit drastic changes in network oscillations. We study, in a spiking neuron network model, the mechanism in adult mice allowing for flexible switches between multiple frequency bands and contrast this to the network structure in juvenile mice that lack this flexibility. The model comprises excitatory pyramidal cells (PCs) and two types of interneurons: the parvalbumin-expressing (PV) and the somatostatinexpressing (SOM) interneuron. In accordance with experimental findings, the pyramidal-PV and pyramidal-SOM cell subnetworks are associated with gamma and beta oscillations, respectively. In our model, they are both generated via a pyramidal-interneuron gamma (PING) mechanism, wherein the PCs drive the oscillations. Furthermore, we demonstrate that large but not small visual stimulation activates SOM cells, which shift the frequency of resting-state gamma oscillations produced by the pyramidal-PV cell subnetwork so that beta rhythms emerge. Finally, we show that this behavior is obtained for only a subset of PV and SOM interneuron projection strengths, indicating that their influence on the PCs should be balanced so that they can compete for oscillatory control of the PCs. In sum, we propose a mechanism by which visual beta rhythms can emerge from spontaneous gamma oscillations in a network model of the mouse V1; for this mechanism to reproduce V1 dynamics in adult mice, balance between the effective strengths of PV and SOM cells is required.

Comprehensive characterization of oscillatory signatures in a model circuit with PV- and SOM-expressing interneurons.

Neural circuits contain a wide variety of interneuron types, which differ in their biophysical properties and connectivity patterns. The two most common interneuron types, parvalbumin-expressing and somatostatin-expressing cells, have been shown to be differentially involved in many cognitive functions. These cell types also show different relationships with the power and phase of oscillations in local field potentials. The mechanisms that underlie the emergence of different oscillatory rhythms in neural circuits with more than one interneuron subtype, and the roles specific interneurons play in those mechanisms, are not fully understood. Here, we present a comprehensive analysis of all possible circuit motifs and input regimes that can be achieved in circuits comprised of excitatory cells, PV-like fast-spiking interneurons and SOM-like low-threshold spiking interneurons. We identify 18 unique motifs and simulate their dynamics over a range of input strengths. Using several characteristics, such as oscillation frequency, firing rates, phase of firing and burst fraction, we cluster the resulting circuit dynamics across motifs in order to identify patterns of activity and compare these patterns to behaviors that were generated in circuits with one interneuron type. In addition to the well-known PING and ING gamma oscillations and an asynchronous state, our analysis identified three oscillatory behaviors that were generated by the three-cell-type motifs only: theta-nested gamma oscillations, stable beta oscillations and theta-locked bursting behavior, which have also been observed in experiments. Our characterization provides a map to interpret experimental activity patterns and suggests pharmacological manipulations or optogenetics approaches to validate these conclusions.

Sensitivity to Stimulus Irregularity Is Inherent in Neural Networks.

Behavior is controlled by complex neural networks in which neurons process thousands of inputs. However, even short spike trains evoked in a single cortical neuron were demonstrated to be sufficient to influence behavior in vivo. Specifically, irregular sequences of interspike intervals (ISIs) had a more reliable influence on behavior despite their resemblance to stochastic activity. Similarly, irregular tactile stimulation led to higher rates of behavioral responses. In this study, we identify the mechanisms enabling this sensitivity to stimulus irregularity (SSI) on the neuronal and network levels using simulated spiking neural networks. Matching in vivo experiments, we find that irregular stimulation elicits more detectable network events (bursts) than regular stimulation. Dissecting the stimuli, we identify short ISIs-occurring more frequently in irregular stimulations-as the main drivers of SSI rather than complex irregularity per se. In addition, we find that short-term plasticity modulates SSI. We subsequently eliminate the different mechanisms in turn to assess their role in generating SSI. Removing inhibitory interneurons, we find that SSI is retained, suggesting that SSI is not dependent on inhibition. Removing recurrency, we find that SSI is retained due to the ability of individual neurons to integrate activity over short timescales ("cell memory"). Removing single-neuron dynamics, we find that SSI is retained based on the short-term retention of activity within the recurrent network structure ("network memory"). Finally, using a further simplified probabilistic model, we find that local network structure is not required for SSI. Hence, SSI is identified as a general property that we hypothesize to be ubiquitous in neural networks with different structures and biophysical properties. Irregular sequences contain shorter ISIs, which are the main drivers underlying SSI. The experimentally observed SSI should thus generalize to other systems, suggesting a functional role for irregular activity in cortex.

Characterization of network structure in stereoEEG data using consensus-based partial coherence.

Coherence is a widely used measure to determine the frequency-resolved functional connectivity between pairs of recording sites, but this measure is confounded by shared inputs to the pair. To remove shared inputs, the 'partial coherence' can be computed by conditioning the spectral matrices of the pair on all other recorded channels, which involves the calculation of a matrix (pseudo-) inverse. It has so far remained a challenge to use the time-resolved partial coherence to analyze intracranial recordings with a large number of recording sites. For instance, calculating the partial coherence using a pseudoinverse method produces a high number of false positives when it is applied to a large number of channels. To address this challenge, we developed a new method that randomly aggregated channels into a smaller number of effective channels on which the calculation of partial coherence was based. We obtained a 'consensus' partial coherence (cPCOH) by repeating this approach for several random aggregations of channels (permutations) and only accepting those activations in time and frequency with a high enough consensus. Using model data we show that the cPCOH method effectively filters out the effect of shared inputs and performs substantially better than the pseudo-inverse. We successfully applied the cPCOH procedure to human stereotactic EEG data and demonstrated three key advantages of this method relative to alternative procedures. First, it reduces the number of false positives relative to the pseudo-inverse method. Second, it allows for titration of the amount of false positives relative to the false negatives by adjusting the consensus threshold, thus allowing the data-analyst to prioritize one over the other to meet specific analysis demands. Third, it substantially reduced the number of identified interactions compared to coherence, providing a sparser network of connections from which clear spatial patterns emerged. These patterns can serve as a starting point of further analyses that provide insight into network dynamics during cognitive processes. These advantages likely generalize to other modalities in which shared inputs introduce confounds, such as electroencephalography (EEG) and magneto-encephalography (MEG).

The missing link: Predicting connectomes from noisy and partially observed tract tracing data.

Our understanding of the wiring map of the brain, known as the connectome, has increased greatly in the last decade, mostly due to technological advancements in neuroimaging techniques and improvements in computational tools to interpret the vast amount of available data. Despite this, with the exception of the C. elegans roundworm, no definitive connectome has been established for any species. In order to obtain this, tracer studies are particularly appealing, as these have proven highly reliable. The downside of tract tracing is that it is costly to perform, and can only be applied ex vivo. In this paper, we suggest that instead of probing all possible connections, hitherto unknown connections may be predicted from the data that is already available. Our approach uses a 'latent space model' that embeds the connectivity in an abstract physical space. Regions that are close in the latent space have a high chance of being connected, while regions far apart are most likely disconnected in the connectome. After learning the latent embedding from the connections that we did observe, the latent space allows us to predict connections that have not been probed previously. We apply the methodology to two connectivity data sets of the macaque, where we demonstrate that the latent space model is successful in predicting unobserved connectivity, outperforming two baselines and an alternative model in nearly all cases. Furthermore, we show how the latent spatial embedding may be used to integrate multimodal observations (i.e. anterograde and retrograde tracers) for the mouse neocortex. Finally, our probabilistic approach enables us to make explicit which connections are easy to predict and which prove difficult, allowing for informed follow-up studies.

Correction: The missing link: Predicting connectomes from noisy and partially observed tract tracing data.

[This corrects the article DOI: 10.1371/journal.pcbi.1005374.].

A Developmental Switch for Hebbian Plasticity.

Hebbian forms of synaptic plasticity are required for the orderly development of sensory circuits in the brain and are powerful modulators of learning and memory in adulthood. During development, emergence of Hebbian plasticity leads to formation of functional circuits. By modeling the dynamics of neurotransmitter release during early postnatal cortical development we show that a developmentally regulated switch in vesicle exocytosis mode triggers associative (i.e. Hebbian) plasticity. Early in development spontaneous vesicle exocytosis (SVE), often considered as 'synaptic noise', is important for homogenization of synaptic weights and maintenance of synaptic weights in the appropriate dynamic range. Our results demonstrate that SVE has a permissive, whereas subsequent evoked vesicle exocytosis (EVE) has an instructive role in the expression of Hebbian plasticity. A timed onset for Hebbian plasticity can be achieved by switching from SVE to EVE and the balance between SVE and EVE can control the effective rate of Hebbian plasticity. We further show that this developmental switch in neurotransmitter release mode enables maturation of spike-timing dependent plasticity. A mis-timed or inadequate SVE to EVE switch may lead to malformation of brain networks thereby contributing to the etiology of neurodevelopmental disorders.

Translating single-neuron axonal reconstructions into meso-scale connectivity statistics in the mouse somatosensory thalamus.

Characterizing the connectomic and morphological diversity of thalamic neurons is key for better understanding how the thalamus relays sensory inputs to the cortex. The recent public release of complete single-neuron morphological reconstructions enables the analysis of previously inaccessible connectivity patterns from individual neurons. Here we focus on the Ventral Posteromedial (VPM) nucleus and characterize the full diversity of 257 VPM neurons, obtained by combining data from the MouseLight and Braintell projects. Neurons were clustered according to their most dominantly targeted cortical area and further subdivided by their jointly targeted areas. We obtained a 2D embedding of morphological diversity using the dissimilarity between all pairs of axonal trees. The curved shape of the embedding allowed us to characterize neurons by a 1-dimensional coordinate. The coordinate values were aligned both with the progression of soma position along the dorsal-ventral and lateral-medial axes and with that of axonal terminals along the posterior-anterior and medial-lateral axes, as well as with an increase in the number of branching points, distance from soma and branching width. Taken together, we have developed a novel workflow for linking three challenging aspects of connectomics, namely the topography, higher order connectivity patterns and morphological diversity, with VPM as a test-case. The workflow is linked to a unified access portal that contains the morphologies and integrated with 2D cortical flatmap and subcortical visualization tools. The workflow and resulting processed data have been made available in Python, and can thus be used for modeling and experimentally validating new hypotheses on thalamocortical connectivity.

Motifs in health and disease: the promise of circuit interrogation by optogenetics.

Identifying the dominant dynamical motifs in cortical circuits and determining their functional relevance is of the utmost importance to understand the underlying mechanisms of psychiatric diseases and to develop effective therapies. Optogenetics can be used to interrogate cortical circuits to determine the dominant motif and thereby identify the relevant biophysical time scales that set the oscillation frequency. We review how computational models of cortical networks can help guide optogenetics experiments. We focus our attention on the pyramidal interneuron gamma motif, which is comprised of reciprocally connected excitatory and inhibitory neurons, and determine how the different biophysical time scales of the circuit components are reflected in the resonance of the power in the local field potential at the frequency of stimulation as a function of that frequency. Cardin et al. [J.A. Cardin et al. (2009)Nature, 459, 663-667] find that periodic stimulation of inhibitory cells leads to a resonance at gamma frequencies (30-80 Hz), but that stimulation of excitatory cells does not lead to a resonance. We can account for these results when the pyramidal cells are endowed with an intrinsic frequency preference due to a slow hyperpolarizing current. Furthermore, when fast α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA)-mediated excitatory currents are replaced by slow N-methyl-d-aspartate (NMDA)-mediated ones in inhibitory cells, the gamma frequency resonance is reduced; however, when the same replacement is made in excitatory cells, gamma oscillations are enhanced. The results are relevant to schizophrenia, because there is evidence that NMDA receptors on parvalbumin-positive cells are primarily affected and that the regulation of gamma oscillations is impaired.

Hedonic and nucleus accumbens neural responses to a natural reward are regulated by aversive conditioning.

The nucleus accumbens (NAc) plays a role in hedonic reactivity to taste stimuli. Learning can alter the hedonic valence of a given stimulus, and it remains unclear how the NAc encodes this shift. The present study examined whether the population response of NAc neurons to a taste stimulus is plastic using a conditioned taste aversion (CTA) paradigm. Electrophysiological and electromyographic (EMG) responses to intraoral infusions of a sucrose (0.3 M) solution were made in naïve rats (Day 1). Immediately following the session, half of the rats (n = 6; Paired) received an injection of lithium chloride (0.15 M; i.p.) to induce malaise and establish a CTA while the other half (n = 6; Unpaired) received a saline injection. Days later (Day 5), NAc recordings during infusions of sucrose were again made. Electrophysiological and EMG responses to sucrose did not differ between groups on Day 1. For both groups, the majority of sucrose responsive neurons exhibited a decrease in firing rate (77% and 71% for Paired and Unpaired, respectively). Following conditioning, in Paired rats, EMG responses were indicative of aversion. Moreover, the majority of responsive NAc neurons now exhibited an increase in firing rate (69%). Responses in Unpaired rats were unchanged by the experience. Thus, the NAc differentially encodes the hedonic value of the same stimulus based on learned associations.

Uncovering Statistical Links Between Gene Expression and Structural Connectivity Patterns in the Mouse Brain.

Finding links between genes and structural connectivity is of the utmost importance for unravelling the underlying mechanism of the brain connectome. In this study we identify links between the gene expression and the axonal projection density in the mouse brain, by applying a modified version of the Linked ICA method to volumetric data from the Allen Institute for Brain Science for identifying independent sources of information that link both modalities at the voxel level. We performed separate analyses on sets of projections from the visual cortex, the caudoputamen and the midbrain reticular nucleus, and we determined those brain areas, injections and genes that were most involved in independent components that link both gene expression and projection density data, while we validated their biological context through enrichment analysis. We identified representative and literature-validated cortico-midbrain and cortico-striatal projections, whose gene subsets were enriched with annotations for neuronal and synaptic function and related developmental and metabolic processes. The results were highly reproducible when including all available projections, as well as consistent with factorisations obtained using the Dictionary Learning and Sparse Coding technique. Hence, Linked ICA yielded reproducible independent components that were preserved under increasing data variance. Taken together, we have developed and validated a novel paradigm for linking gene expression and structural projection patterns in the mouse mesoconnectome, which can power future studies aiming to relate genes to brain function.

Progress towards a cellularly resolved mouse mesoconnectome is empowered by data fusion and new neuroanatomy techniques.

Over the past decade there has been a rapid improvement in techniques for obtaining large-scale cellular level data related to the mouse brain connectome. However, a detailed mapping of cell-type-specific projection patterns is lacking, which would, for instance, allow us to study the role of circuit motifs in cognitive processes. In this work, we review advanced neuroanatomical and data fusion techniques within the context of a proposed Multimodal Connectomic Integration Framework for augmenting the cellularly resolved mouse mesoconnectome. First, we emphasize the importance of registering data modalities to a common reference atlas. We then review a number of novel experimental techniques that can provide data for characterizing cell-types in the mouse brain. Furthermore, we examine a number of data integration strategies, which involve fine-grained cell-type classification, spatial inference of cell densities, latent variable models for the mesoconnectome and multi-modal factorisation. Finally, we discuss a number of use cases which depend on connectome augmentation techniques, such as model simulations of functional connectivity and generating mechanistic hypotheses for animal disease models.

Where is Cingulate Cortex? A Cross-Species View.

To compare findings across species, neuroscience relies on cross-species homologies, particularly in terms of brain areas. For cingulate cortex, a structure implicated in behavioural adaptation and control, a homologous definition across mammals is available - but currently not employed by most rodent researchers. The standard partitioning of rodent cingulate cortex is inconsistent with that in any other model species, including humans. Reviewing the existing literature, we show that the homologous definition better aligns results of rodent studies with those of other species, and reveals a clearer structural and functional organisation within rodent cingulate cortex itself. Based on these insights, we call for widespread adoption of the homologous nomenclature, and reinterpretation of previous studies originally based on the nonhomologous partitioning of rodent cingulate cortex.

Human stereoEEG recordings reveal network dynamics of decision-making in a rule-switching task.

The processing steps that lead up to a decision, i.e., the transformation of sensory evidence into motor output, are not fully understood. Here, we combine stereoEEG recordings from the human cortex, with single-lead and time-resolved decoding, using a wide range of temporal frequencies, to characterize decision processing during a rule-switching task. Our data reveal the contribution of rostral inferior parietal lobule (IPL) regions, in particular PFt, and the parietal opercular regions in decision processing and demonstrate that the network representing the decision is common to both task rules. We reconstruct the sequence in which regions engage in decision processing on single trials, thereby providing a detailed picture of the network dynamics involved in decision-making. The reconstructed timeline suggests that the supramarginal gyrus in IPL links decision regions in prefrontal cortex with premotor regions, where the motor plan for the response is elaborated.

Connectomic Analysis of Brain Networks: Novel Techniques and Future Directions.

Brain networks, localized or brain-wide, exist only at the cellular level, i.e., between specific pre- and post-synaptic neurons, which are connected through functionally diverse synapses located at specific points of their cell membranes. "Connectomics" is the emerging subfield of neuroanatomy explicitly aimed at elucidating the wiring of brain networks with cellular resolution and a quantified accuracy. Such data are indispensable for realistic modeling of brain circuitry and function. A connectomic analysis, therefore, needs to identify and measure the soma, dendrites, axonal path, and branching patterns together with the synapses and gap junctions of the neurons involved in any given brain circuit or network. However, because of the submicron caliber, 3D complexity, and high packing density of most such structures, as well as the fact that axons frequently extend over long distances to make synapses in remote brain regions, creating connectomic maps is technically challenging and requires multi-scale approaches, Such approaches involve the combination of the most sensitive cell labeling and analysis methods available, as well as the development of new ones able to resolve individual cells and synapses with increasing high-throughput. In this review, we provide an overview of recently introduced high-resolution methods, which researchers wanting to enter the field of connectomics may consider. It includes several molecular labeling tools, some of which specifically label synapses, and covers a number of novel imaging tools such as brain clearing protocols and microscopy approaches. Apart from describing the tools, we also provide an assessment of their qualities. The criteria we use assess the qualities that tools need in order to contribute to deciphering the key levels of circuit organization. We conclude with a brief future outlook for neuroanatomic research, computational methods, and network modeling, where we also point out several outstanding issues like structure-function relations and the complexity of neural models.

Discovering spike patterns in neuronal responses.

When a cortical neuron is repeatedly injected with the same fluctuating current stimulus (frozen noise) the timing of the spikes is highly precise from trial to trial and the spike pattern appears to be unique. We show here that the same repeated stimulus can produce more than one reliable temporal pattern of spikes. A new method is introduced to find these patterns in raw multitrial data and is tested on surrogate data sets. Using it, multiple coexisting spike patterns were discovered in pyramidal cells recorded from rat prefrontal cortex in vitro, in data obtained in vivo from the middle temporal area of the monkey (Buracas et al., 1998) and from the cat lateral geniculate nucleus (Reinagel and Reid, 2002). The spike patterns lasted from a few tens of milliseconds in vitro to several seconds in vivo. We conclude that the prestimulus history of a neuron may influence the precise timing of the spikes in response to a stimulus over a wide range of time scales.

Circuit to construct mapping: a mathematical tool for assisting the diagnosis and treatment in major depressive disorder.

Major depressive disorder (MDD) is a serious condition with a lifetime prevalence exceeding 16% worldwide. MDD is a heterogeneous disorder that involves multiple behavioral symptoms on the one hand and multiple neuronal circuits on the other hand. In this review, we integrate the literature on cognitive and physiological biomarkers of MDD with the insights derived from mathematical models of brain networks, especially models that can be used for fMRI datasets. We refer to the recent NIH research domain criteria initiative, in which a concept of "constructs" as functional units of mental disorders is introduced. Constructs are biomarkers present at multiple levels of brain functioning - cognition, genetics, brain anatomy, and neurophysiology. In this review, we propose a new approach which we called circuit to construct mapping (CCM), which aims to characterize causal relations between the underlying network dynamics (as the cause) and the constructs referring to the clinical symptoms of MDD (as the effect). CCM involves extracting diagnostic categories from behavioral data, linking circuits that are causal to these categories with use of clinical neuroimaging data, and modeling the dynamics of the emerging circuits with attractor dynamics in order to provide new, neuroimaging-related biomarkers for MDD. The CCM approach optimizes the clinical diagnosis and patient stratification. It also addresses the recent demand for linking circuits to behavior, and provides a new insight into clinical treatment by investigating the dynamics of neuronal circuits underneath cognitive dimensions of MDD. CCM can serve as a new regime toward personalized medicine, assisting the diagnosis and treatment of MDD.

Reliability and bifurcation in neurons driven by multiple sinusoids.

The brain produces dynamical rhythms at many frequencies that shift in amplitude and phase. To understand the functional consequences of mixtures of oscillations at the single cell level, we recorded the spike trains from single rat cortical neurons in vitro in response to two mixed sine wave currents. The reliability of spike timing was measured as a function of the relative power, phase and frequencies of the sine wave mixture. Peaks in the reliability were observed at a preferred phase difference, frequency and relative power. These results have a natural interpretation in terms of spike train attractors and bifurcations.

Synchronization as a mechanism for attentional gain modulation.

Responses of neurons in monkey visual cortex are modulated when attention is directed into the receptive field of the neuron: the gain or sensitivity of the response is increased or the synchronization of the spikes to the local field potential (LFP) is increased. We investigated, using model simulations, whether the synchrony of inhibitory networks could link these observations. We found that, indeed, an increase in inhibitory synchrony could enhance the coherence of the model neurons with the simulated LFP, and could have different effects on the firing rate. When the firing rate vs. current (f-I) response curves saturated at high I, attention yielded a shift in sensitivity; alternatively, when the f-I curves were non-saturating, the most significant effect was on the gain of the response. This suggests that attention may act through changes in the synchrony of inhibitory networks.