What does the mean mean? A simple test for neuroscience.

Trial-averaged metrics, e.g. tuning curves or population response vectors, are a ubiquitous way of characterizing neuronal activity. But how relevant are such trial-averaged responses to neuronal computation itself? Here we present a simple test to estimate whether average responses reflect aspects of neuronal activity that contribute to neuronal processing. The test probes two assumptions implicitly made whenever average metrics are treated as meaningful representations of neuronal activity: Reliability: Neuronal responses repeat consistently enough across trials that they convey a recognizable reflection of the average response to downstream regions.Behavioural relevance: If a single-trial response is more similar to the average template, it is more likely to evoke correct behavioural responses. We apply this test to two data sets: (1) Two-photon recordings in primary somatosensory cortices (S1 and S2) of mice trained to detect optogenetic stimulation in S1; and (2) Electrophysiological recordings from 71 brain areas in mice performing a contrast discrimination task. Under the highly controlled settings of Data set 1, both assumptions were largely fulfilled. In contrast, the less restrictive paradigm of Data set 2 met neither assumption. Simulations predict that the larger diversity of neuronal response preferences, rather than higher cross-trial reliability, drives the better performance of Data set 1. We conclude that when behaviour is less tightly restricted, average responses do not seem particularly relevant to neuronal computation, potentially because information is encoded more dynamically. Most importantly, we encourage researchers to apply this simple test of computational relevance whenever using trial-averaged neuronal metrics, in order to gauge how representative cross-trial averages are in a given context.

Analyzing Thalamocortical Tract-Tracing Experiments in a Common Reference Space.

Current mesoscale connectivity atlases provide limited information about the organization of thalamocortical projections in the mouse brain. Labeling the projections of spatially restricted neuron populations in thalamus can provide a functionally relevant level of connectomic analysis, but these need to be integrated within the same common reference space. Here, we present a pipeline for the segmentation, registration, integration and analysis of multiple tract-tracing experiments. The key difference with other workflows is that the data is transformed to fit the reference template. As a test-case, we investigated the axonal projections and intranuclear arrangement of seven neuronal populations of the ventral posteromedial nucleus of the thalamus (VPM), which we labeled with an anterograde tracer. Their soma positions corresponded, from dorsal to ventral, to cortical representations of the whiskers, nose and mouth. They strongly targeted layer 4, with the majority exclusively targeting one cortical area and the ones in ventrolateral VPM branching to multiple somatosensory areas. We found that our experiments were more topographically precise than similar experiments from the Allen Institute and projections to the primary somatosensory area were in agreement with single-neuron morphological reconstructions from publicly available databases. This pilot study sets the basis for a shared virtual connectivity atlas that could be enriched with additional data for studying the topographical organization of different thalamic nuclei. The pipeline is accessible with only minimal programming skills via a Jupyter Notebook, and offers multiple visualization tools such as cortical flatmaps, subcortical plots and 3D renderings and can be used with custom anatomical delineations.

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.

A Standards Organization for Open and FAIR Neuroscience: the International Neuroinformatics Coordinating Facility.

There is great need for coordination around standards and best practices in neuroscience to support efforts to make neuroscience a data-centric discipline. Major brain initiatives launched around the world are poised to generate huge stores of neuroscience data. At the same time, neuroscience, like many domains in biomedicine, is confronting the issues of transparency, rigor, and reproducibility. Widely used, validated standards and best practices are key to addressing the challenges in both big and small data science, as they are essential for integrating diverse data and for developing a robust, effective, and sustainable infrastructure to support open and reproducible neuroscience. However, developing community standards and gaining their adoption is difficult. The current landscape is characterized both by a lack of robust, validated standards and a plethora of overlapping, underdeveloped, untested and underutilized standards and best practices. The International Neuroinformatics Coordinating Facility (INCF), an independent organization dedicated to promoting data sharing through the coordination of infrastructure and standards, has recently implemented a formal procedure for evaluating and endorsing community standards and best practices in support of the FAIR principles. By formally serving as a standards organization dedicated to open and FAIR neuroscience, INCF helps evaluate, promulgate, and coordinate standards and best practices across neuroscience. Here, we provide an overview of the process and discuss how neuroscience can benefit from having a dedicated standards body.

Geodesic-based distance reveals nonlinear topological features in neural activity from mouse visual cortex.

An increasingly popular approach to the analysis of neural data is to treat activity patterns as being constrained to and sampled from a manifold, which can be characterized by its topology. The persistent homology method identifies the type and number of holes in the manifold, thereby yielding functional information about the coding and dynamic properties of the underlying neural network. In this work, we give examples of highly nonlinear manifolds in which the persistent homology algorithm fails when it uses the Euclidean distance because it does not always yield a good approximation to the true distance distribution of a point cloud sampled from a manifold. To deal with this issue, we instead estimate the geodesic distance which is a better approximation of the true distance distribution and can therefore be used to successfully identify highly nonlinear features with persistent homology. To document the utility of the method, we utilize a toy model comprised of a circular manifold, built from orthogonal sinusoidal coordinate functions and show how the choice of metric determines the performance of the persistent homology algorithm. Furthermore, we explore the robustness of the method across different manifold properties, like the number of samples, curvature and amount of added noise. We point out strategies for interpreting its results as well as some possible pitfalls of its application. Subsequently, we apply this analysis to neural data coming from the Visual Coding-Neuropixels dataset recorded at the Allen Institute in mouse visual cortex in response to stimulation with drifting gratings. We find that different manifolds with a non-trivial topology can be seen across regions and stimulus properties. Finally, we interpret how these changes in manifold topology along with stimulus parameters and cortical region inform how the brain performs visual computation.

Learning at Variable Attentional Load Requires Cooperation of Working Memory, Meta-learning, and Attention-augmented Reinforcement Learning.

Flexible learning of changing reward contingencies can be realized with different strategies. A fast learning strategy involves using working memory of recently rewarded objects to guide choices. A slower learning strategy uses prediction errors to gradually update value expectations to improve choices. How the fast and slow strategies work together in scenarios with real-world stimulus complexity is not well known. Here, we aim to disentangle their relative contributions in rhesus monkeys while they learned the relevance of object features at variable attentional load. We found that learning behavior across six monkeys is consistently best predicted with a model combining (i) fast working memory and (ii) slower reinforcement learning from differently weighted positive and negative prediction errors as well as (iii) selective suppression of nonchosen feature values and (iv) a meta-learning mechanism that enhances exploration rates based on a memory trace of recent errors. The optimal model parameter settings suggest that these mechanisms cooperate differently at low and high attentional loads. Whereas working memory was essential for efficient learning at lower attentional loads, enhanced weighting of negative prediction errors and meta-learning were essential for efficient learning at higher attentional loads. Together, these findings pinpoint a canonical set of learning mechanisms and suggest how they may cooperate when subjects flexibly adjust to environments with variable real-world attentional demands.

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.

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.

Interneuron-specific gamma synchronization indexes cue uncertainty and prediction errors in lateral prefrontal and anterior cingulate cortex.

Inhibitory interneurons are believed to realize critical gating functions in cortical circuits, but it has been difficult to ascertain the content of gated information for well-characterized interneurons in primate cortex. Here, we address this question by characterizing putative interneurons in primate prefrontal and anterior cingulate cortex while monkeys engaged in attention demanding reversal learning. We find that subclasses of narrow spiking neurons have a relative suppressive effect on the local circuit indicating they are inhibitory interneurons. One of these interneuron subclasses showed prominent firing rate modulations and (35-45 Hz) gamma synchronous spiking during periods of uncertainty in both, lateral prefrontal cortex (LPFC) and anterior cingulate cortex (ACC). In LPFC, this interneuron subclass activated when the uncertainty of attention cues was resolved during flexible learning, whereas in ACC it fired and gamma-synchronized when outcomes were uncertain and prediction errors were high during learning. Computational modeling of this interneuron-specific gamma band activity in simple circuit motifs suggests it could reflect a soft winner-take-all gating of information having high degree of uncertainty. Together, these findings elucidate an electrophysiologically characterized interneuron subclass in the primate, that forms gamma synchronous networks in two different areas when resolving uncertainty during adaptive goal-directed behavior.

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.

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.

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.

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.

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.

Prediction of a Cell-Class-Specific Mouse Mesoconnectome Using Gene Expression Data.

Reconstructing brain connectivity at sufficient resolution for computational models designed to study the biophysical mechanisms underlying cognitive processes is extremely challenging. For such a purpose, a mesoconnectome that includes laminar and cell-class specificity would be a major step forward. We analyzed the ability of gene expression patterns to predict cell-class and layer-specific projection patterns and assessed the functional annotations of the most predictive groups of genes. To achieve our goal we used publicly available volumetric gene expression and connectivity data and we trained computational models to learn and predict cell-class and layer-specific axonal projections using gene expression data. Predictions were done in two ways, namely predicting projection strengths using the expression of individual genes and using the co-expression of genes organized in spatial modules, as well as predicting binary forms of projection. For predicting the strength of projections, we found that ridge (L2-regularized) regression had the highest cross-validated accuracy with a median r2 score of 0.54 which corresponded for binarized predictions to a median area under the ROC value of 0.89. Next, we identified 200 spatial gene modules using a dictionary learning and sparse coding approach. We found that these modules yielded predictions of comparable accuracy, with a median r2 score of 0.51. Finally, a gene ontology enrichment analysis of the most predictive gene groups resulted in significant annotations related to postsynaptic function. Taken together, we have demonstrated a prediction workflow that can be used to perform multimodal data integration to improve the accuracy of the predicted mesoconnectome and support other neuroscience use cases.

Adaptive spike-artifact removal from local field potentials uncovers prominent beta and gamma band neuronal synchronization.

BACKGROUND: Many neurons synchronize their action potentials to the phase of local field potential (LFP) fluctuations in one or more frequency bands. Analyzing this spike-to-LFP synchronization is challenging, however, when neural spikes and LFP are generated in the same local circuit, because the spike's action potential waveform leak into the LFP and distort phase synchrony estimates. Existing approaches to address this spike bleed-through artifact relied on removing the average action potential waveforms of neurons, but this leaves artifacts in the LFP and distorts synchrony estimates. NEW METHOD: We describe a spike-removal method that surpasses these limitations by decomposing individual action potentials into their frequency components before their removal from the LFP. The adaptively estimated frequency components allow for variable spread, strength and temporal variation of the spike artifact. RESULTS: This adaptive approach effectively removes spike bleed-through artifacts in synthetic data with known ground truth, and in single neuron and LFP recordings in nonhuman primate striatum. For a large population of neurons with both narrow and broad action potential waveforms, the use of adaptive artifact removal uncovered 20-35 Hz beta and 35-45 Hz gamma band spike-LFP synchronization that would have remained contaminated otherwise. COMPARISON WITH EXISTING METHODS: We demonstrate that adaptive spike-artifact removal cleans LFP data that remained contaminated when applying existing Bayesian and non-Bayesian methods of average spike-artifact removal. CONCLUSIONS: Applying adaptive spike-removal from field potentials allows to estimate the phase at which neurons synchronize and the consistency of their phase-locked firing for both beta and low gamma frequencies. These metrics may prove essential to understand cell-to-circuit neuronal interactions in multiple brain systems.

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.