Ultra-high density electrodes improve detection, yield, and cell type identification in neuronal recordings

To understand the neural basis of behavior, it is essential to sensitively and accurately measure neural activity at single neuron and single spike resolution. Extracellular electrophysiology delivers this, but it has biases in the neurons it detects and it imperfectly resolves their action potentials. To minimize these limitations, we developed a silicon probe with much smaller and denser recording sites than previous designs, called Neuropixels Ultra (NP Ultra). This device samples neuronal activity at ultra-high spatial density (~10 times higher than previous probes) with low noise levels, while trading off recording span. NP Ultra is effectively an implantable voltage-sensing camera that captures a planar image of a neuron's electrical field. We use a spike sorting algorithm optimized for these probes to demonstrate that the yield of visually-responsive neurons in recordings from mouse visual cortex improves up to ~3-fold. We show that NP Ultra can record from small neuronal structures including axons and dendrites. Recordings across multiple brain regions and four species revealed a subset of extracellular action potentials with unexpectedly small spatial spread and axon-like features. We share a large-scale dataset of these brain-wide recordings in mice as a resource for studies of neuronal biophysics. Finally, using ground-truth identification of three major inhibitory cortical cell types, we found that these cell types were discriminable with approximately 75% success, a significant improvement over lower-resolution recordings. NP Ultra improves spike sorting performance, detection of subcellular compartments, and cell type classification to enable more powerful dissection of neural circuit activity during behavior.

Lightning Pose: improved animal pose estimation via semi-supervised learning, Bayesian ensembling, and cloud-native open-source tools

Contemporary pose estimation methods enable precise measurements of behavior via supervised deep learning with hand-labeled video frames. Although effective in many cases, the supervised approach requires extensive labeling and often produces outputs that are unreliable for downstream analyses. Here, we introduce "Lightning Pose," an efficient pose estimation package with three algorithmic contributions. First, in addition to training on a few labeled video frames, we use many unlabeled videos and penalize the network whenever its predictions violate motion continuity, multiple-view geometry, and posture plausibility (semi-supervised learning). Second, we introduce a network architecture that resolves occlusions by predicting pose on any given frame using surrounding unlabeled frames. Third, we refine the pose predictions post-hoc by combining ensembling and Kalman smoothing. Together, these components render pose trajectories more accurate and scientifically usable. We release a cloud application that allows users to label data, train networks, and predict new videos directly from the browser.

Identifying Interpretable Latent Factors with Sparse Component Analysis

In many neural populations, the computationally relevant signals are posited to be a set of 'latent factors' - signals shared across many individual neurons. Understanding the relationship between neural activity and behavior requires the identification of factors that reflect distinct computational roles. Methods for identifying such factors typically require supervision, which can be suboptimal if one is unsure how (or whether) factors can be grouped into distinct, meaningful sets. Here, we introduce Sparse Component Analysis (SCA), an unsupervised method that identifies interpretable latent factors. SCA seeks factors that are sparse in time and occupy orthogonal dimensions. With these simple constraints, SCA facilitates surprisingly clear parcellations of neural activity across a range of behaviors. We applied SCA to motor cortex activity from reaching and cycling monkeys, single-trial imaging data from C. elegans, and activity from a multitask artificial network. SCA consistently identified sets of factors that were useful in describing network computations.

DREDge: robust motion correction for high-density extracellular recordings across species

High-density microelectrode arrays (MEAs) have opened new possibilities for systems neuroscience in human and non-human animals, but brain tissue motion relative to the array poses a challenge for downstream analyses, particularly in human recordings. We introduce DREDge (Decentralized Registration of Electrophysiology Data), a robust algorithm which is well suited for the registration of noisy, nonstationary extracellular electrophysiology recordings. In addition to estimating motion from spikes in the action potential (AP) frequency band, DREDge enables automated tracking of motion at high temporal resolution in the local field potential (LFP) frequency band. In human intraoperative recordings, which often feature fast (period <1s) motion, DREDge correction in the LFP band enabled reliable recovery of evoked potentials, and significantly reduced single-unit spike shape variability and spike sorting error. Applying DREDge to recordings made during deep probe insertions in nonhuman primates demonstrated the possibility of tracking probe motion of centimeters across several brain regions while simultaneously mapping single unit electrophysiological features. DREDge reliably delivered improved motion correction in acute mouse recordings, especially in those made with an recent ultra-high density probe. We also implemented a procedure for applying DREDge to recordings made across tens of days in chronic implantations in mice, reliably yielding stable motion tracking despite changes in neural activity across experimental sessions. Together, these advances enable automated, scalable registration of electrophysiological data across multiple species, probe types, and drift cases, providing a stable foundation for downstream scientific analyses of these rich datasets.

Bayesian target optimisation for high-precision holographic optogenetics

Two-photon optogenetics has transformed our ability to probe the structure and function of neural circuits. However, achieving precise optogenetic control of neural ensemble activity has remained fundamentally constrained by the problem of off-target stimulation (OTS): the inadvertent activation of nearby non-target neurons due to imperfect confinement of light onto target neurons. Here we propose a novel computational approach to this problem called Bayesian target optimisation. Our approach uses nonparametric Bayesian inference to model neural responses to optogenetic stimulation, and then optimises the laser powers and optical target locations needed to achieve a desired activity pattern with minimal OTS. We validate our approach in simulations and using data from in vitro experiments, showing that Bayesian target optimisation considerably reduces OTS across all conditions we test. Together, these results establish our ability to overcome OTS, enabling optogenetic stimulation with substantially improved precision.

Bypassing spike sorting: Density-based decoding using spike localization from dense multielectrode probes

Neural decoding and its applications to brain computer interfaces (BCI) are essential for understanding the association between neural activity and behavior. A prerequisite for many decoding approaches is spike sorting, the assignment of action potentials (spikes) to individual neurons. Current spike sorting algorithms, however, can be inaccurate and do not properly model uncertainty of spike assignments, therefore discarding information that could potentially improve decoding performance. Recent advances in high-density probes (e.g., Neuropixels) and computational methods now allow for extracting a rich set of spike features from unsorted data; these features can in turn be used to directly decode behavioral correlates. To this end, we propose a spike sorting-free decoding method that directly models the distribution of extracted spike features using a mixture of Gaussians (MoG) encoding the uncertainty of spike assignments, without aiming to solve the spike clustering problem explicitly. We allow the mixing proportion of the MoG to change over time in response to the behavior and develop variational inference methods to fit the resulting model and to perform decoding. We benchmark our method with an extensive suite of recordings from different animals and probe geometries, demonstrating that our proposed decoder can consistently outperform current methods based on thresholding (i.e. multi-unit activity) and spike sorting. Open source code is available at https://github.com/yzhang511/density_decoding.

maskNMF: A denoise-sparsen-detect approach for extracting neural signals from dense imaging data

A number of calcium imaging methods have been developed to monitor the activity of large populations of neurons. One particularly promising approach, Bessel imaging, captures neural activity from a volume by projecting within the imaged volume onto a single imaging plane, therefore effectively mixing signals and increasing the number of neurons imaged per pixel. These signals must then be computationally demixed to recover the desired neural activity. Unfortunately, currently-available demixing methods can perform poorly in the regime of high imaging density (i.e., many neurons per pixel). In this work we introduce a new pipeline (maskNMF) for demixing dense calcium imaging data. The main idea is to first denoise and temporally sparsen the observed video; this enhances signal strength and reduces spatial overlap significantly. Next we detect neurons in the sparsened video using a neural network trained on a library of neural shapes. These shapes are derived from segmented electron microscopy images input into a Bessel imaging model; therefore no manual selection of "good" neural shapes from the functional data is required here. After cells are detected, we use a constrained non-negative matrix factorization approach to demix the activity, using the detected cells' shapes to initialize the factorization. We test the resulting pipeline on both simulated and real datasets and find that it is able to achieve accurate demixing on denser data than was previously feasible, therefore enabling faithful imaging of larger neural populations. The method also provides good results on more "standard" two-photon imaging data. Finally, because much of the pipeline operates on a significantly compressed version of the raw data and is highly parallelizable, the algorithm is fast, processing large datasets faster than real time.

The spatial and temporal structure of neural activity across the fly brain

What are the spatial and temporal scales of brainwide neuronal activity? We used swept, confocally-aligned planar excitation (SCAPE) microscopy to image all cells in a large volume of the brain of adult Drosophila with high spatiotemporal resolution while flies engaged in a variety of spontaneous behaviors. This revealed neural representations of behavior on multiple spatial and temporal scales. The activity of most neurons correlated (or anticorrelated) with running and flailing over timescales that ranged from seconds to a minute. Grooming elicited a weaker global response. Significant residual activity not directly correlated with behavior was high dimensional and reflected the activity of small clusters of spatially organized neurons that may correspond to genetically defined cell types. These clusters participate in the global dynamics, indicating that neural activity reflects a combination of local and broadly distributed components. This suggests that microcircuits with highly specified functions are provided with knowledge of the larger context in which they operate.

MULTIMODAL MICROSCOPY IMAGE ALIGNMENT USING SPATIAL AND SHAPE INFORMATION AND A BRANCH-AND-BOUND ALGORITHM

Multimodal microscopy experiments that image the same population of cells under different experimental conditions have become a widely used approach in systems and molecular neuroscience. The main obstacle is to align the different imaging modalities to obtain complementary information about the observed cell population (e.g., gene expression and calcium signal). Traditional image registration methods perform poorly when only a small subset of cells are present in both images, as is common in multimodal experiments. We cast multimodal microscopy alignment as a cell subset matching problem. To solve this non-convex problem, we introduce an efficient and globally optimal branch-and-bound algorithm to find subsets of point clouds that are in rotational alignment with each other. In addition, we use complementary information about cell shape and location to compute the matching likelihood of cell pairs in two imaging modalities to further prune the optimization search tree. Finally, we use the maximal set of cells in rigid rotational alignment to seed image deformation fields to obtain a final registration result. Our framework performs better than the state-of-the-art histology alignment approaches regarding matching quality and is faster than manual alignment, providing a viable solution to improve the throughput of multimodal microscopy experiments.

Learning Probabilistic Piecewise Rigid Atlases of Model Organisms via Generative Deep Networks

Atlases are crucial to imaging statistics as they enable the standardization of inter-subject and inter-population analyses. While existing atlas estimation methods based on fluid/elastic/diffusion registration yield high-quality results for the human brain, these deformation models do not extend to a variety of other challenging areas of neuroscience such as the anatomy of C. elegans worms and fruit flies. To this end, this work presents a general probabilistic deep network-based framework for atlas estimation and registration which can flexibly incorporate various deformation models and levels of keypoint supervision that can be applied to a wide class of model organisms. Of particular relevance, it also develops a deformable piecewise rigid atlas model which is regularized to preserve inter-observation distances between neighbors. These modeling considerations are shown to improve atlas construction and key-point alignment across a diversity of datasets with small sample sizes including neuron positions in C. elegans hermaphrodites, fluorescence microscopy of male C. elegans, and images of fruit fly wings. Code is accessible at https://github.com/amin-nejat/Deformable-Atlas.

ROBUST ONLINE MULTIBAND DRIFT ESTIMATION IN ELECTROPHYSIOLOGY DATA

High-density electrophysiology probes have opened new possibilities for systems neuroscience in human and non-human animals, but probe motion poses a challenge for downstream analyses, particularly in human recordings. We improve on the state of the art for tracking this motion with four major contributions. First, we extend previous decentralized methods to use multiband information, leveraging the local field potential (LFP) in addition to spikes. Second, we show that the LFP-based approach enables registration at sub-second temporal resolution. Third, we introduce an efficient online motion tracking algorithm, enabling the method to scale up to longer and higher-resolution recordings, and possibly facilitating real-time applications. Finally, we improve the robustness of the approach by introducing a structure-aware objective and simple methods for adaptive parameter selection. Together, these advances enable fully automated scalable registration of challenging datasets from human and mouse.

Neuroscience Cloud Analysis As a Service: An open-source platform for scalable, reproducible data analysis

A key aspect of neuroscience research is the development of powerful, general-purpose data analyses that process large datasets. Unfortunately, modern data analyses have a hidden dependence upon complex computing infrastructure (e.g., software and hardware), which acts as an unaddressed deterrent to analysis users. Although existing analyses are increasingly shared as open-source software, the infrastructure and knowledge needed to deploy these analyses efficiently still pose significant barriers to use. In this work, we develop Neuroscience Cloud Analysis As a Service (NeuroCAAS): a fully automated open-source analysis platform offering automatic infrastructure reproducibility for any data analysis. We show how NeuroCAAS supports the design of simpler, more powerful data analyses and that many popular data analysis tools offered through NeuroCAAS outperform counterparts on typical infrastructure. Pairing rigorous infrastructure management with cloud resources, NeuroCAAS dramatically accelerates the dissemination and use of new data analyses for neuroscientific discovery.

Blind demixing methods for recovering dense neuronal morphology from barcode imaging data

Cellular barcoding methods offer the exciting possibility of 'infinite-pseudocolor' anatomical reconstruction-i.e., assigning each neuron its own random unique barcoded 'pseudocolor,' and then using these pseudocolors to trace the microanatomy of each neuron. Here we use simulations, based on densely-reconstructed electron microscopy microanatomy, with signal structure matched to real barcoding data, to quantify the feasibility of this procedure. We develop a new blind demixing approach to recover the barcodes that label each neuron, and validate this method on real data with known barcodes. We also develop a neural network which uses the recovered barcodes to reconstruct the neuronal morphology from the observed fluorescence imaging data, 'connecting the dots' between discontiguous barcode amplicon signals. We find that accurate recovery should be feasible, provided that the barcode signal density is sufficiently high. This study suggests the possibility of mapping the morphology and projection pattern of many individual neurons simultaneously, at high resolution and at large scale, via conventional light microscopy.

Reconstruction of neocortex: Organelles, compartments, cells, circuits, and activity

We assembled a semi-automated reconstruction of L2/3 mouse primary visual cortex from ∼250 × 140 × 90 μm3 of electron microscopic images, including pyramidal and non-pyramidal neurons, astrocytes, microglia, oligodendrocytes and precursors, pericytes, vasculature, nuclei, mitochondria, and synapses. Visual responses of a subset of pyramidal cells are included. The data are publicly available, along with tools for programmatic and three-dimensional interactive access. Brief vignettes illustrate the breadth of potential applications relating structure to function in cortical circuits and neuronal cell biology. Mitochondria and synapse organization are characterized as a function of path length from the soma. Pyramidal connectivity motif frequencies are predicted accurately using a configuration model of random graphs. Pyramidal cells receiving more connections from nearby cells exhibit stronger and more reliable visual responses. Sample code shows data access and analysis.

Visualizing the organization and differentiation of the male-specific nervous system of C. elegans

Sex differences in the brain are prevalent throughout the animal kingdom and particularly well appreciated in the nematode Caenorhabditis elegans, where male animals contain a little-studied set of 93 male-specific neurons. To make these neurons amenable for future study, we describe here how a multicolor reporter transgene, NeuroPAL, is capable of visualizing the distinct identities of all male-specific neurons. We used NeuroPAL to visualize and characterize a number of features of the male-specific nervous system. We provide several proofs of concept for using NeuroPAL to identify the sites of expression of gfp-tagged reporter genes and for cellular fate analysis by analyzing the effect of removal of several developmental patterning genes on neuronal identity acquisition. We use NeuroPAL and its intrinsic cohort of more than 40 distinct differentiation markers to show that, even though male-specific neurons are generated throughout all four larval stages, they execute their terminal differentiation program in a coordinated manner in the fourth larval stage. This coordinated wave of differentiation, which we call 'just-in-time' differentiation, couples neuronal maturation programs with the appearance of sexual organs.

Partitioning variability in animal behavioral videos using semi-supervised variational autoencoders

Recent neuroscience studies demonstrate that a deeper understanding of brain function requires a deeper understanding of behavior. Detailed behavioral measurements are now often collected using video cameras, resulting in an increased need for computer vision algorithms that extract useful information from video data. Here we introduce a new video analysis tool that combines the output of supervised pose estimation algorithms (e.g. DeepLabCut) with unsupervised dimensionality reduction methods to produce interpretable, low-dimensional representations of behavioral videos that extract more information than pose estimates alone. We demonstrate this tool by extracting interpretable behavioral features from videos of three different head-fixed mouse preparations, as well as a freely moving mouse in an open field arena, and show how these interpretable features can facilitate downstream behavioral and neural analyses. We also show how the behavioral features produced by our model improve the precision and interpretation of these downstream analyses compared to using the outputs of either fully supervised or fully unsupervised methods alone.

Non-parametric Vignetting Correction for Sparse Spatial Transcriptomics Images

Spatial transcriptomics techniques such as STARmap [15] enable the subcellular detection of RNA transcripts within complex tissue sections. The data from these techniques are impacted by optical microscopy limitations, such as shading or vignetting effects from uneven illumination during image capture. Downstream analysis of these sparse spatially resolved transcripts is dependent upon the correction of these artefacts. This paper introduces a novel non-parametric vignetting correction tool for spatial transcriptomic images, which estimates the illumination field and background using an efficient iterative sliced histogram normalization routine. We show that our method outperforms the state-of-the-art shading correction techniques both in terms of illumination and background field estimation and requires fewer input images to perform the estimation adequately. We further demonstrate an important downstream application of our technique, showing that spatial transcriptomic volumes corrected by our method yield a higher and more uniform gene expression spot-calling in the rodent hippocampus. Python code and a demo file to reproduce our results are provided in the supplementary material and at this github page: https://github.com/BoveyRao/Non-parametric-vc-for-sparse-st.

Chronic, cortex-wide imaging of specific cell populations during behavior

Measurements of neuronal activity across brain areas are important for understanding the neural correlates of cognitive and motor processes such as attention, decision-making and action selection. However, techniques that allow cellular resolution measurements are expensive and require a high degree of technical expertise, which limits their broad use. Wide-field imaging of genetically encoded indicators is a high-throughput, cost-effective and flexible approach to measure activity of specific cell populations with high temporal resolution and a cortex-wide field of view. Here we outline our protocol for assembling a wide-field macroscope setup, performing surgery to prepare the intact skull and imaging neural activity chronically in behaving, transgenic mice. Further, we highlight a processing pipeline that leverages novel, cloud-based methods to analyze large-scale imaging datasets. The protocol targets laboratories that are seeking to build macroscopes, optimize surgical procedures for long-term chronic imaging and/or analyze cortex-wide neuronal recordings. The entire protocol, including steps for assembly and calibration of the macroscope, surgical preparation, imaging and data analysis, requires a total of 8 h. It is designed to be accessible to laboratories with limited expertise in imaging methods or interest in high-throughput imaging during behavior.

Nonlinear Decoding of Natural Images From Large-Scale Primate Retinal Ganglion Recordings

Decoding sensory stimuli from neural activity can provide insight into how the nervous system might interpret the physical environment, and facilitates the development of brain-machine interfaces. Nevertheless, the neural decoding problem remains a significant open challenge. Here, we present an efficient nonlinear decoding approach for inferring natural scene stimuli from the spiking activities of retinal ganglion cells (RGCs). Our approach uses neural networks to improve on existing decoders in both accuracy and scalability. Trained and validated on real retinal spike data from more than 1000 simultaneously recorded macaque RGC units, the decoder demonstrates the necessity of nonlinear computations for accurate decoding of the fine structures of visual stimuli. Specifically, high-pass spatial features of natural images can only be decoded using nonlinear techniques, while low-pass features can be extracted equally well by linear and nonlinear methods. Together, these results advance the state of the art in decoding natural stimuli from large populations of neurons.

High-fidelity estimates of spikes and subthreshold waveforms from 1-photon voltage imaging in vivo

The ability to probe the membrane potential of multiple genetically defined neurons simultaneously would have a profound impact on neuroscience research. Genetically encoded voltage indicators are a promising tool for this purpose, and recent developments have achieved a high signal-to-noise ratio in vivo with 1-photon fluorescence imaging. However, these recordings exhibit several sources of noise and signal extraction remains a challenge. We present an improved signal extraction pipeline, spike-guided penalized matrix decomposition-nonnegative matrix factorization (SGPMD-NMF), which resolves supra- and subthreshold voltages in vivo. The method incorporates biophysical and optical constraints. We validate the pipeline with simultaneous patch-clamp and optical recordings from mouse layer 1 in vivo and with simulated and composite datasets with realistic noise. We demonstrate applications to mouse hippocampus expressing paQuasAr3-s or SomArchon1, mouse cortex expressing SomArchon1 or Voltron, and zebrafish spines expressing zArchon1.

BARcode DEmixing through Non-negative Spatial Regression (BarDensr)

Modern spatial transcriptomics methods can target thousands of different types of RNA transcripts in a single slice of tissue. Many biological applications demand a high spatial density of transcripts relative to the imaging resolution, leading to partial mixing of transcript rolonies in many voxels; unfortunately, current analysis methods do not perform robustly in this highly-mixed setting. Here we develop a new analysis approach, BARcode DEmixing through Non-negative Spatial Regression (BarDensr): we start with a generative model of the physical process that leads to the observed image data and then apply sparse convex optimization methods to estimate the underlying (demixed) rolony densities. We apply BarDensr to simulated and real data and find that it achieves state of the art signal recovery, particularly in densely-labeled regions or data with low spatial resolution. Finally, BarDensr is fast and parallelizable. We provide open-source code as well as an implementation for the 'NeuroCAAS' cloud platform.

NeuroPAL: A Multicolor Atlas for Whole-Brain Neuronal Identification in C. elegans

Comprehensively resolving neuronal identities in whole-brain images is a major challenge. We achieve this in C. elegans by engineering a multicolor transgene called NeuroPAL (a neuronal polychromatic atlas of landmarks). NeuroPAL worms share a stereotypical multicolor fluorescence map for the entire hermaphrodite nervous system that resolves all neuronal identities. Neurons labeled with NeuroPAL do not exhibit fluorescence in the green, cyan, or yellow emission channels, allowing the transgene to be used with numerous reporters of gene expression or neuronal dynamics. We showcase three applications that leverage NeuroPAL for nervous-system-wide neuronal identification. First, we determine the brainwide expression patterns of all metabotropic receptors for acetylcholine, GABA, and glutamate, completing a map of this communication network. Second, we uncover changes in cell fate caused by transcription factor mutations. Third, we record brainwide activity in response to attractive and repulsive chemosensory cues, characterizing multimodal coding for these stimuli.

Reinforcement Learning Recruits Somata and Apical Dendrites across Layers of Primary Sensory Cortex

The mammalian brain can form associations between behaviorally relevant stimuli in an animal’s environment. While such learning is thought to primarily involve high-order association cortex, even primary sensory areas receive long-range connections carrying information that could contribute to high-level representations. Here, we imaged layer 1 apical dendrites in the barrel cortex of mice performing a whisker-based operant behavior. In addition to sensory-motor events, calcium signals in apical dendrites of layers 2/3 and 5 neurons and in layer 2/3 somata track the delivery of rewards, both choice related and randomly administered. Reward-related tuft-wide dendritic spikes emerge gradually with training and are task specific. Learning recruits cells whose intrinsic activity coincides with the time of reinforcement. Layer 4 largely lacked reward-related signals, suggesting a source other than the primary thalamus. Our results demonstrate that a sensory cortex can acquire a set of associations outside its immediate sensory modality and linked to salient behavioral events.

Previously, the only known triggers of apical dendritic spikes were “bottom-up”events, such as appropriate sensory stimuli or an animal’s location in space. Lacefield et al. show that reinforced associations are powerful triggers of apical dendrite activity and that reward can manipulate perceptions at their earliest stages of cortical processing.

Localized semi-nonnegative matrix factorization (LocaNMF) of widefield calcium imaging data

Widefield calcium imaging enables recording of large-scale neural activity across the mouse dorsal cortex. In order to examine the relationship of these neural signals to the resulting behavior, it is critical to demix the recordings into meaningful spatial and temporal components that can be mapped onto well-defined brain regions. However, no current tools satisfactorily extract the activity of the different brain regions in individual mice in a data-driven manner, while taking into account mouse-specific and preparation-specific differences. Here, we introduce Localized semi-Nonnegative Matrix Factorization (LocaNMF), a method that efficiently decomposes widefield video data and allows us to directly compare activity across multiple mice by outputting mouse-specific localized functional regions that are significantly more interpretable than more traditional decomposition techniques. Moreover, it provides a natural subspace to directly compare correlation maps and neural dynamics across different behaviors, mice, and experimental conditions, and enables identification of task- and movement-related brain regions.

Bright and photostable chemigenetic indicators for extended in vivo voltage imaging

Genetically encoded voltage indicators (GEVIs) enable monitoring of neuronal activity at high spatial and temporal resolution. However, the utility of existing GEVIs has been limited by the brightness and photostability of fluorescent proteins and rhodopsins. We engineered a GEVI, called Voltron, that uses bright and photostable synthetic dyes instead of protein-based fluorophores, thereby extending the number of neurons imaged simultaneously in vivo by a factor of 10 and enabling imaging for significantly longer durations relative to existing GEVIs. We used Voltron for in vivo voltage imaging in mice, zebrafish, and fruit flies. In the mouse cortex, Voltron allowed single-trial recording of spikes and subthreshold voltage signals from dozens of neurons simultaneously over a 15-minute period of continuous imaging. In larval zebrafish, Voltron enabled the precise correlation of spike timing with behavior.

Complementary networks of cortical somatostatin interneurons enforce layer specific control

The neocortex is functionally organized into layers. Layer four receives the densest bottom up sensory inputs, while layers 2/3 and 5 receive top down inputs that may convey predictive information. A subset of cortical somatostatin (SST) neurons, the Martinotti cells, gate top down input by inhibiting the apical dendrites of pyramidal cells in layers 2/3 and 5, but it is unknown whether an analogous inhibitory mechanism controls activity in layer 4. Using high precision circuit mapping, in vivo optogenetic perturbations, and single cell transcriptional profiling, we reveal complementary circuits in the mouse barrel cortex involving genetically distinct SST subtypes that specifically and reciprocally interconnect with excitatory cells in different layers: Martinotti cells connect with layers 2/3 and 5, whereas non-Martinotti cells connect with layer 4. By enforcing layer-specific inhibition, these parallel SST subnetworks could independently regulate the balance between bottom up and top down input.

Neural data science: accelerating the experiment-analysis-theory cycle in large-scale neuroscience

Modern large-scale multineuronal recording methodologies, including multielectrode arrays, calcium imaging, and optogenetic techniques, produce single-neuron resolution data of a magnitude and precision that were the realm of science fiction twenty years ago. The major bottlenecks in systems and circuit neuroscience no longer lie in simply collecting data from large neural populations, but also in understanding this data: developing novel scientific questions, with corresponding analysis techniques and experimental designs to fully harness these new capabilities and meaningfully interrogate these questions. Advances in methods for signal processing, network analysis, dimensionality reduction, and optimal control-developed in lockstep with advances in experimental neurotechnology-promise major breakthroughs in multiple fundamental neuroscience problems. These trends are clear in a broad array of subfields of modern neuroscience; this review focuses on recent advances in methods for analyzing neural time-series data with single-neuronal precision.

Community-based benchmarking improves spike rate inference from two-photon calcium imaging data

In recent years, two-photon calcium imaging has become a standard tool to probe the function of neural circuits and to study computations in neuronal populations. However, the acquired signal is only an indirect measurement of neural activity due to the comparatively slow dynamics of fluorescent calcium indicators. Different algorithms for estimating spike rates from noisy calcium measurements have been proposed in the past, but it is an open question how far performance can be improved. Here, we report the results of the spikefinder challenge, launched to catalyze the development of new spike rate inference algorithms through crowd-sourcing. We present ten of the submitted algorithms which show improved performance compared to previously evaluated methods. Interestingly, the top-performing algorithms are based on a wide range of principles from deep neural networks to generative models, yet provide highly correlated estimates of the neural activity. The competition shows that benchmark challenges can drive algorithmic developments in neuroscience.

Efficient and accurate extraction of in vivo calcium signals from microendoscopic video data

In vivo calcium imaging through microendoscopic lenses enables imaging of previously inaccessible neuronal populations deep within the brains of freely moving animals. However, it is computationally challenging to extract single-neuronal activity from microendoscopic data, because of the very large background fluctuations and high spatial overlaps intrinsic to this recording modality. Here, we describe a new constrained matrix factorization approach to accurately separate the background and then demix and denoise the neuronal signals of interest. We compared the proposed method against previous independent components analysis and constrained nonnegative matrix factorization approaches. On both simulated and experimental data recorded from mice, our method substantially improved the quality of extracted cellular signals and detected more well-isolated neural signals, especially in noisy data regimes. These advances can in turn significantly enhance the statistical power of downstream analyses, and ultimately improve scientific conclusions derived from microendoscopic data.

Anxiety Cells in a Hippocampal-Hypothalamic Circuit

The hippocampus is traditionally thought to transmit contextual information to limbic structures where it acquires valence. Using freely moving calcium imaging and optogenetics, we show that while the dorsal CA1 subregion of the hippocampus is enriched in place cells, ventral CA1 (vCA1) is enriched in anxiety cells that are activated by anxiogenic environments and required for avoidance behavior. Imaging cells defined by their projection target revealed that anxiety cells were enriched in the vCA1 population projecting to the lateral hypothalamic area (LHA) but not to the basal amygdala (BA). Consistent with this selectivity, optogenetic activation of vCA1 terminals in LHA but not BA increased anxiety and avoidance, while activation of terminals in BA but not LHA impaired contextual fear memory. Thus, the hippocampus encodes not only neutral but also valence-related contextual information, and the vCA1-LHA pathway is a direct route by which the hippocampus can rapidly influence innate anxiety behavior.

The central amygdala controls learning in the lateral amygdala

Experience-driven synaptic plasticity in the lateral amygdala is thought to underlie the formation of associations between sensory stimuli and an ensuing threat. However, how the central amygdala participates in such a learning process remains unclear. Here we show that PKC-δ-expressing central amygdala neurons are essential for the synaptic plasticity underlying learning in the lateral amygdala, as they convey information about the unconditioned stimulus to lateral amygdala neurons during fear conditioning.

Electrical stimulus artifact cancellation and neural spike detection on large multi-electrode arrays

Simultaneous electrical stimulation and recording using multi-electrode arrays can provide a valuable technique for studying circuit connectivity and engineering neural interfaces. However, interpreting these measurements is challenging because the spike sorting process (identifying and segregating action potentials arising from different neurons) is greatly complicated by electrical stimulation artifacts across the array, which can exhibit complex and nonlinear waveforms, and overlap temporarily with evoked spikes. Here we develop a scalable algorithm based on a structured Gaussian Process model to estimate the artifact and identify evoked spikes. The effectiveness of our methods is demonstrated in both real and simulated 512-electrode recordings in the peripheral primate retina with single-electrode and several types of multi-electrode stimulation. We establish small error rates in the identification of evoked spikes, with a computational complexity that is compatible with real-time data analysis. This technology may be helpful in the design of future high-resolution sensory prostheses based on tailored stimulation (e.g., retinal prostheses), and for closed-loop neural stimulation at a much larger scale than currently possible.

The Spatiotemporal Organization of the Striatum Encodes Action Space

Activity in striatal direct- and indirect-pathway spiny projection neurons (SPNs) is critical for proper movement. However, little is known about the spatiotemporal organization of this activity. We investigated the spatiotemporal organization of SPN ensemble activity in mice during self-paced, natural movements using microendoscopic imaging. Activity in both pathways showed predominantly local but also some long-range correlations. Using a novel approach to cluster and quantify behaviors based on continuous accelerometer and video data, we found that SPN ensembles active during specific actions were spatially closer and more correlated overall. Furthermore, similarity between different actions corresponded to the similarity between SPN ensemble patterns, irrespective of movement speed. Consistently, the accuracy of decoding behavior from SPN ensemble patterns was directly related to the dissimilarity between behavioral clusters. These results identify a predominantly local, but not spatially compact, organization of direct- and indirect-pathway SPN activity that maps action space independently of movement speed.

•

Direct- and indirect-pathway SPNs show locally biased spatiotemporal organization

•

SPNs active during a particular behavior are more correlated and spatially closer

•

SPN ensembles encode action identity independently of movement speed

•

Distance between behaviors corresponds to distance between SPN ensemble patterns

Multi-scale approaches for high-speed imaging and analysis of large neural populations

Progress in modern neuroscience critically depends on our ability to observe the activity of large neuronal populations with cellular spatial and high temporal resolution. However, two bottlenecks constrain efforts towards fast imaging of large populations. First, the resulting large video data is challenging to analyze. Second, there is an explicit tradeoff between imaging speed, signal-to-noise, and field of view: with current recording technology we cannot image very large neuronal populations with simultaneously high spatial and temporal resolution. Here we describe multi-scale approaches for alleviating both of these bottlenecks. First, we show that spatial and temporal decimation techniques based on simple local averaging provide order-of-magnitude speedups in spatiotemporally demixing calcium video data into estimates of single-cell neural activity. Second, once the shapes of individual neurons have been identified at fine scale (e.g., after an initial phase of conventional imaging with standard temporal and spatial resolution), we find that the spatial/temporal resolution tradeoff shifts dramatically: after demixing we can accurately recover denoised fluorescence traces and deconvolved neural activity of each individual neuron from coarse scale data that has been spatially decimated by an order of magnitude. This offers a cheap method for compressing this large video data, and also implies that it is possible to either speed up imaging significantly, or to "zoom out" by a corresponding factor to image order-of-magnitude larger neuronal populations with minimal loss in accuracy or temporal resolution.

Population-Level Representation of a Temporal Sequence Underlying Song Production in the Zebra Finch

The zebra finch brain features a set of clearly defined and hierarchically arranged motor nuclei that are selectively responsible for producing singing behavior. One of these regions, a critical forebrain structure called HVC, contains premotor neurons that are active at precise time points during song production. However, the neural representation of this behavior at a population level remains elusive. We used two-photon microscopy to monitor ensemble activity during singing, integrating across multiple trials by adopting a Bayesian inference approach to more precisely estimate burst timing. Additionally, we examined spiking and motor-related synaptic inputs using intracellular recordings during singing. With both experimental approaches, we find that premotor events do not occur preferentially at the onsets or offsets of song syllables or at specific subsyllabic motor landmarks. These results strongly support the notion that HVC projection neurons collectively exhibit a temporal sequence during singing that is uncoupled from ongoing movements.

Cerebellar granule cells acquire a widespread predictive feedback signal during motor learning

Cerebellar granule cells, which constitute half the brain's neurons, supply Purkinje cells with contextual information necessary for motor learning, but how they encode this information is unknown. Here we show, using two-photon microscopy to track neural activity over multiple days of cerebellum-dependent eyeblink conditioning in mice, that granule cell populations acquire a dense representation of the anticipatory eyelid movement. Initially, granule cells responded to neutral visual and somatosensory stimuli as well as periorbital airpuffs used for training. As learning progressed, two-thirds of monitored granule cells acquired a conditional response whose timing matched or preceded the learned eyelid movements. Granule cell activity covaried trial by trial to form a redundant code. Many granule cells were also active during movements of nearby body structures. Thus, a predictive signal about the upcoming movement is widely available at the input stage of the cerebellar cortex, as required by forward models of cerebellar control.

25th Annual Computational Neuroscience Meeting: CNS-2016

A1 Functional advantages of cell-type heterogeneity in neural circuits

Tatyana O. Sharpee

A2 Mesoscopic modeling of propagating waves in visual cortex

Alain Destexhe

A3 Dynamics and biomarkers of mental disorders

Mitsuo Kawato

F1 Precise recruitment of spiking output at theta frequencies requires dendritic h-channels in multi-compartment models of oriens-lacunosum/moleculare hippocampal interneurons

Vladislav Sekulić, Frances K. Skinner

F2 Kernel methods in reconstruction of current sources from extracellular potentials for single cells and the whole brains

Daniel K. Wójcik, Chaitanya Chintaluri, Dorottya Cserpán, Zoltán Somogyvári

F3 The synchronized periods depend on intracellular transcriptional repression mechanisms in circadian clocks.

Jae Kyoung Kim, Zachary P. Kilpatrick, Matthew R. Bennett, Kresimir Josić

O1 Assessing irregularity and coordination of spiking-bursting rhythms in central pattern generators

Irene Elices, David Arroyo, Rafael Levi, Francisco B. Rodriguez, Pablo Varona

O2 Regulation of top-down processing by cortically-projecting parvalbumin positive neurons in basal forebrain

Eunjin Hwang, Bowon Kim, Hio-Been Han, Tae Kim, James T. McKenna, Ritchie E. Brown, Robert W. McCarley, Jee Hyun Choi

O3 Modeling auditory stream segregation, build-up and bistability

James Rankin, Pamela Osborn Popp, John Rinzel

O4 Strong competition between tonotopic neural ensembles explains pitch-related dynamics of auditory cortex evoked fields

Alejandro Tabas, André Rupp, Emili Balaguer-Ballester

O5 A simple model of retinal response to multi-electrode stimulation

Matias I. Maturana, David B. Grayden, Shaun L. Cloherty, Tatiana Kameneva, Michael R. Ibbotson, Hamish Meffin

O6 Noise correlations in V4 area correlate with behavioral performance in visual discrimination task

Veronika Koren, Timm Lochmann, Valentin Dragoi, Klaus Obermayer

O7 Input-location dependent gain modulation in cerebellar nucleus neurons

Maria Psarrou, Maria Schilstra, Neil Davey, Benjamin Torben-Nielsen, Volker Steuber

O8 Analytic solution of cable energy function for cortical axons and dendrites

Huiwen Ju, Jiao Yu, Michael L. Hines, Liang Chen, Yuguo Yu

O9 C. elegans interactome: interactive visualization of Caenorhabditis elegans worm neuronal network

Jimin Kim, Will Leahy, Eli Shlizerman

O10 Is the model any good? Objective criteria for computational neuroscience model selection

Justas Birgiolas, Richard C. Gerkin, Sharon M. Crook

O11 Cooperation and competition of gamma oscillation mechanisms

Atthaphon Viriyopase, Raoul-Martin Memmesheimer, Stan Gielen

O12 A discrete structure of the brain waves

Yuri Dabaghian, Justin DeVito, Luca Perotti

O13 Direction-specific silencing of the Drosophila gaze stabilization system

Anmo J. Kim, Lisa M. Fenk, Cheng Lyu, Gaby Maimon

O14 What does the fruit fly think about values? A model of olfactory associative learning

Chang Zhao, Yves Widmer, Simon Sprecher,Walter Senn

O15 Effects of ionic diffusion on power spectra of local field potentials (LFP)

Geir Halnes, Tuomo Mäki-Marttunen, Daniel Keller, Klas H. Pettersen,Ole A. Andreassen, Gaute T. Einevoll

O16 Large-scale cortical models towards understanding relationship between brain structure abnormalities and cognitive deficits

Yasunori Yamada

O17 Spatial coarse-graining the brain: origin of minicolumns

Moira L. Steyn-Ross, D. Alistair Steyn-Ross

O18 Modeling large-scale cortical networks with laminar structure

Jorge F. Mejias, John D. Murray, Henry Kennedy, Xiao-Jing Wang

O19 Information filtering by partial synchronous spikes in a neural population

Alexandra Kruscha, Jan Grewe, Jan Benda, Benjamin Lindner

O20 Decoding context-dependent olfactory valence in Drosophila

Laurent Badel, Kazumi Ohta, Yoshiko Tsuchimoto, Hokto Kazama

P1 Neural network as a scale-free network: the role of a hub

B. Kahng

P2 Hemodynamic responses to emotions and decisions using near-infrared spectroscopy optical imaging

Nicoladie D. Tam

P3 Phase space analysis of hemodynamic responses to intentional movement directions using functional near-infrared spectroscopy (fNIRS) optical imaging technique

Nicoladie D.Tam, Luca Pollonini, George Zouridakis

P4 Modeling jamming avoidance of weakly electric fish

Jaehyun Soh, DaeEun Kim

P5 Synergy and redundancy of retinal ganglion cells in prediction

Minsu Yoo, S. E. Palmer

P6 A neural field model with a third dimension representing cortical depth

Viviana Culmone, Ingo Bojak

P7 Network analysis of a probabilistic connectivity model of the Xenopus tadpole spinal cord

Andrea Ferrario, Robert Merrison-Hort, Roman Borisyuk

P8 The recognition dynamics in the brain

Chang Sub Kim

P9 Multivariate spike train analysis using a positive definite kernel

Taro Tezuka

P10 Synchronization of burst periods may govern slow brain dynamics during general anesthesia

Pangyu Joo

P11 The ionic basis of heterogeneity affects stochastic synchrony

Young-Ah Rho, Shawn D. Burton, G. Bard Ermentrout, Jaeseung Jeong, Nathaniel N. Urban

P12 Circular statistics of noise in spike trains with a periodic component

Petr Marsalek

P14 Representations of directions in EEG-BCI using Gaussian readouts

Hoon-Hee Kim, Seok-hyun Moon, Do-won Lee, Sung-beom Lee, Ji-yong Lee, Jaeseung Jeong

P15 Action selection and reinforcement learning in basal ganglia during reaching movements

Yaroslav I. Molkov, Khaldoun Hamade, Wondimu Teka, William H. Barnett, Taegyo Kim, Sergey Markin, Ilya A. Rybak

P17 Axon guidance: modeling axonal growth in T-Junction assay

Csaba Forro, Harald Dermutz, László Demkó, János Vörös

P19 Transient cell assembly networks encode persistent spatial memories

Yuri Dabaghian, Andrey Babichev

P20 Theory of population coupling and applications to describe high order correlations in large populations of interacting neurons

Haiping Huang

P21 Design of biologically-realistic simulations for motor control

Sergio Verduzco-Flores

P22 Towards understanding the functional impact of the behavioural variability of neurons

Filipa Dos Santos, Peter Andras

P23 Different oscillatory dynamics underlying gamma entrainment deficits in schizophrenia

Christoph Metzner, Achim Schweikard, Bartosz Zurowski

P24 Memory recall and spike frequency adaptation

James P. Roach, Leonard M. Sander, Michal R. Zochowski

P25 Stability of neural networks and memory consolidation preferentially occur near criticality

Quinton M. Skilling, Nicolette Ognjanovski, Sara J. Aton, Michal Zochowski

P26 Stochastic Oscillation in Self-Organized Critical States of Small Systems: Sensitive Resting State in Neural Systems

Sheng-Jun Wang, Guang Ouyang, Jing Guang, Mingsha Zhang, K. Y. Michael Wong, Changsong Zhou

P27 Neurofield: a C++ library for fast simulation of 2D neural field models

Peter A. Robinson, Paula Sanz-Leon, Peter M. Drysdale, Felix Fung, Romesh G. Abeysuriya, Chris J. Rennie, Xuelong Zhao

P28 Action-based grounding: Beyond encoding/decoding in neural code

Yoonsuck Choe, Huei-Fang Yang

P29 Neural computation in a dynamical system with multiple time scales

Yuanyuan Mi, Xiaohan Lin, Si Wu

P30 Maximum entropy models for 3D layouts of orientation selectivity

Joscha Liedtke, Manuel Schottdorf, Fred Wolf

P31 A behavioral assay for probing computations underlying curiosity in rodents

Yoriko Yamamura, Jeffery R. Wickens

P32 Using statistical sampling to balance error function contributions to optimization of conductance-based models

Timothy Rumbell, Julia Ramsey, Amy Reyes, Danel Draguljić, Patrick R. Hof, Jennifer Luebke, Christina M. Weaver

P33 Exploration and implementation of a self-growing and self-organizing neuron network building algorithm

Hu He, Xu Yang, Hailin Ma, Zhiheng Xu, Yuzhe Wang

P34 Disrupted resting state brain network in obese subjects: a data-driven graph theory analysis

Kwangyeol Baek, Laurel S. Morris, Prantik Kundu, Valerie Voon

P35 Dynamics of cooperative excitatory and inhibitory plasticity

Everton J. Agnes, Tim P. Vogels

P36 Frequency-dependent oscillatory signal gating in feed-forward networks of integrate-and-fire neurons

William F. Podlaski, Tim P. Vogels

P37 Phenomenological neural model for adaptation of neurons in area IT

Martin Giese, Pradeep Kuravi, Rufin Vogels

P38 ICGenealogy: towards a common topology of neuronal ion channel function and genealogy in model and experiment

Alexander Seeholzer, William Podlaski, Rajnish Ranjan, Tim Vogels

P39 Temporal input discrimination from the interaction between dynamic synapses and neural subthreshold oscillations

Joaquin J. Torres, Fabiano Baroni, Roberto Latorre, Pablo Varona

P40 Different roles for transient and sustained activity during active visual processing

Bart Gips, Eric Lowet, Mark J. Roberts, Peter de Weerd, Ole Jensen, Jan van der Eerden

P41 Scale-free functional networks of 2D Ising model are highly robust against structural defects: neuroscience implications

Abdorreza Goodarzinick, Mohammad D. Niry, Alireza Valizadeh

P42 High frequency neuron can facilitate propagation of signal in neural networks

Aref Pariz, Shervin S. Parsi, Alireza Valizadeh

P43 Investigating the effect of Alzheimer’s disease related amyloidopathy on gamma oscillations in the CA1 region of the hippocampus

Julia M. Warburton, Lucia Marucci, Francesco Tamagnini, Jon Brown, Krasimira Tsaneva-Atanasova

P44 Long-tailed distributions of inhibitory and excitatory weights in a balanced network with eSTDP and iSTDP

Florence I. Kleberg, Jochen Triesch

P45 Simulation of EMG recording from hand muscle due to TMS of motor cortex

Bahar Moezzi, Nicolangelo Iannella, Natalie Schaworonkow, Lukas Plogmacher, Mitchell R. Goldsworthy, Brenton Hordacre, Mark D. McDonnell, Michael C. Ridding, Jochen Triesch

P46 Structure and dynamics of axon network formed in primary cell culture

Martin Zapotocky, Daniel Smit, Coralie Fouquet, Alain Trembleau

P47 Efficient signal processing and sampling in random networks that generate variability

Sakyasingha Dasgupta, Isao Nishikawa, Kazuyuki Aihara, Taro Toyoizumi

P48 Modeling the effect of riluzole on bursting in respiratory neural networks

Daniel T. Robb, Nick Mellen, Natalia Toporikova

P49 Mapping relaxation training using effective connectivity analysis

Rongxiang Tang, Yi-Yuan Tang

P50 Modeling neuron oscillation of implicit sequence learning

Guangsheng Liang, Seth A. Kiser, James H. Howard, Jr., Yi-Yuan Tang

P51 The role of cerebellar short-term synaptic plasticity in the pathology and medication of downbeat nystagmus

Julia Goncharenko, Neil Davey, Maria Schilstra, Volker Steuber

P52 Nonlinear response of noisy neurons

Sergej O. Voronenko, Benjamin Lindner

P53 Behavioral embedding suggests multiple chaotic dimensions underlie C. elegans locomotion

Tosif Ahamed, Greg Stephens

P54 Fast and scalable spike sorting for large and dense multi-electrodes recordings

Pierre Yger, Baptiste Lefebvre, Giulia Lia Beatrice Spampinato, Elric Esposito, Marcel Stimberg et Olivier Marre

P55 Sufficient sampling rates for fast hand motion tracking

Hansol Choi, Min-Ho Song

P56 Linear readout of object manifolds

SueYeon Chung, Dan D. Lee, Haim Sompolinsky

P57 Differentiating models of intrinsic bursting and rhythm generation of the respiratory pre-Bötzinger complex using phase response curves

Ryan S. Phillips, Jeffrey Smith

P58 The effect of inhibitory cell network interactions during theta rhythms on extracellular field potentials in CA1 hippocampus

Alexandra Pierri Chatzikalymniou, Katie Ferguson, Frances K. Skinner

P59 Expansion recoding through sparse sampling in the cerebellar input layer speeds learning

N. Alex Cayco Gajic, Claudia Clopath, R. Angus Silver

P60 A set of curated cortical models at multiple scales on Open Source Brain

Padraig Gleeson, Boris Marin, Sadra Sadeh, Adrian Quintana, Matteo Cantarelli, Salvador Dura-Bernal, William W. Lytton, Andrew Davison, R. Angus Silver

P61 A synaptic story of dynamical information encoding in neural adaptation

Luozheng Li, Wenhao Zhang, Yuanyuan Mi, Dahui Wang, Si Wu

P62 Physical modeling of rule-observant rodent behavior

Youngjo Song, Sol Park, Ilhwan Choi, Jaeseung Jeong, Hee-sup Shin

P64 Predictive coding in area V4 and prefrontal cortex explains dynamic discrimination of partially occluded shapes

Hannah Choi, Anitha Pasupathy, Eric Shea-Brown

P65 Stability of FORCE learning on spiking and rate-based networks

Dongsung Huh, Terrence J. Sejnowski

P66 Stabilising STDP in striatal neurons for reliable fast state recognition in noisy environments

Simon M. Vogt, Arvind Kumar, Robert Schmidt

P67 Electrodiffusion in one- and two-compartment neuron models for characterizing cellular effects of electrical stimulation

Stephen Van Wert, Steven J. Schiff

P68 STDP improves speech recognition capabilities in spiking recurrent circuits parameterized via differential evolution Markov Chain Monte Carlo

Richard Veale, Matthias Scheutz

P69 Bidirectional transformation between dominant cortical neural activities and phase difference distributions

Sang Wan Lee

P70 Maturation of sensory networks through homeostatic structural plasticity

Júlia Gallinaro, Stefan Rotter

P71 Corticothalamic dynamics: structure, number of solutions and stability of steady-state solutions in the space of synaptic couplings

Paula Sanz-Leon, Peter A. Robinson

P72 Optogenetic versus electrical stimulation of the parkinsonian basal ganglia. Computational study

Leonid L. Rubchinsky, Chung Ching Cheung, Shivakeshavan Ratnadurai-Giridharan

P73 Exact spike-timing distribution reveals higher-order interactions of neurons

Safura Rashid Shomali, Majid Nili Ahmadabadi, Hideaki Shimazaki, S. Nader Rasuli

P74 Neural mechanism of visual perceptual learning using a multi-layered neural network

Xiaochen Zhao, Malte J. Rasch

P75 Inferring collective spiking dynamics from mostly unobserved systems

Jens Wilting, Viola Priesemann

P76 How to infer distributions in the brain from subsampled observations

Anna Levina, Viola Priesemann

P77 Influences of embedding and estimation strategies on the inferred memory of single spiking neurons

Lucas Rudelt, Joseph T. Lizier, Viola Priesemann

P78 A nearest-neighbours based estimator for transfer entropy between spike trains

Joseph T. Lizier, Richard E. Spinney, Mikail Rubinov, Michael Wibral, Viola Priesemann

P79 Active learning of psychometric functions with multinomial logistic models

Ji Hyun Bak, Jonathan Pillow

P81 Inferring low-dimensional network dynamics with variational latent Gaussian process

Yuan Zaho, Il Memming Park

P82 Computational investigation of energy landscapes in the resting state subcortical brain network

Jiyoung Kang, Hae-Jeong Park

P83 Local repulsive interaction between retinal ganglion cells can generate a consistent spatial periodicity of orientation map

Jaeson Jang, Se-Bum Paik

P84 Phase duration of bistable perception reveals intrinsic time scale of perceptual decision under noisy condition

Woochul Choi, Se-Bum Paik

P85 Feedforward convergence between retina and primary visual cortex can determine the structure of orientation map

Changju Lee, Jaeson Jang, Se-Bum Paik

P86 Computational method classifying neural network activity patterns for imaging data

Min Song, Hyeonsu Lee, Se-Bum Paik

P87 Symmetry of spike-timing-dependent-plasticity kernels regulates volatility of memory

Youngjin Park, Woochul Choi^, Se-Bum Paik

P88 Effects of time-periodic coupling strength on the first-spike latency dynamics of a scale-free network of stochastic Hodgkin-Huxley neurons

Ergin Yilmaz, Veli Baysal, Mahmut Ozer

P89 Spectral properties of spiking responses in V1 and V4 change within the trial and are highly relevant for behavioral performance

Veronika Koren, Klaus Obermayer

P90 Methods for building accurate models of individual neurons

Daniel Saska, Thomas Nowotny

P91 A full size mathematical model of the early olfactory system of honeybees

Ho Ka Chan, Alan Diamond, Thomas Nowotny

P92 Stimulation-induced tuning of ongoing oscillations in spiking neural networks

Christoph S. Herrmann, Micah M. Murray, Silvio Ionta, Axel Hutt, Jérémie Lefebvre

P93 Decision-specific sequences of neural activity in balanced random networks driven by structured sensory input

Philipp Weidel, Renato Duarte, Abigail Morrison

P94 Modulation of tuning induced by abrupt reduction of SST cell activity

Jung H. Lee, Ramakrishnan Iyer, Stefan Mihalas

P95 The functional role of VIP cell activation during locomotion

Jung H. Lee, Ramakrishnan Iyer, Christof Koch, Stefan Mihalas

P96 Stochastic inference with spiking neural networks

Mihai A. Petrovici, Luziwei Leng, Oliver Breitwieser, David Stöckel, Ilja Bytschok, Roman Martel, Johannes Bill, Johannes Schemmel, Karlheinz Meier

P97 Modeling orientation-selective electrical stimulation with retinal prostheses

Timothy B. Esler, Anthony N. Burkitt, David B. Grayden, Robert R. Kerr, Bahman Tahayori, Hamish Meffin

P98 Ion channel noise can explain firing correlation in auditory nerves

Bahar Moezzi, Nicolangelo Iannella, Mark D. McDonnell

P99 Limits of temporal encoding of thalamocortical inputs in a neocortical microcircuit

Max Nolte, Michael W. Reimann, Eilif Muller, Henry Markram

P100 On the representation of arm reaching movements: a computational model

Antonio Parziale, Rosa Senatore, Angelo Marcelli

P101 A computational model for investigating the role of cerebellum in acquisition and retention of motor behavior

Rosa Senatore, Antonio Parziale, Angelo Marcelli

P102 The emergence of semantic categories from a large-scale brain network of semantic knowledge

K. Skiker, M. Maouene

P103 Multiscale modeling of M1 multitarget pharmacotherapy for dystonia

Samuel A. Neymotin, Salvador Dura-Bernal, Alexandra Seidenstein, Peter Lakatos, Terence D. Sanger, William W. Lytton

P104 Effect of network size on computational capacity

Salvador Dura-Bernal, Rosemary J. Menzies, Campbell McLauchlan, Sacha J. van Albada, David J. Kedziora, Samuel Neymotin, William W. Lytton, Cliff C. Kerr

P105 NetPyNE: a Python package for NEURON to facilitate development and parallel simulation of biological neuronal networks

Salvador Dura-Bernal, Benjamin A. Suter, Samuel A. Neymotin, Cliff C. Kerr, Adrian Quintana, Padraig Gleeson, Gordon M. G. Shepherd, William W. Lytton

P107 Inter-areal and inter-regional inhomogeneity in co-axial anisotropy of Cortical Point Spread in human visual areas

Juhyoung Ryu, Sang-Hun Lee

P108 Two bayesian quanta of uncertainty explain the temporal dynamics of cortical activity in the non-sensory areas during bistable perception

Joonwon Lee, Sang-Hun Lee

P109 Optimal and suboptimal integration of sensory and value information in perceptual decision making

Hyang Jung Lee, Sang-Hun Lee

P110 A Bayesian algorithm for phoneme Perception and its neural implementation

Daeseob Lim, Sang-Hun Lee

P111 Complexity of EEG signals is reduced during unconsciousness induced by ketamine and propofol

Jisung Wang, Heonsoo Lee

P112 Self-organized criticality of neural avalanche in a neural model on complex networks

Nam Jung, Le Anh Quang, Seung Eun Maeng, Tae Ho Lee, Jae Woo Lee

P113 Dynamic alterations in connection topology of the hippocampal network during ictal-like epileptiform activity in an in vitro rat model

Chang-hyun Park, Sora Ahn, Jangsup Moon, Yun Seo Choi, Juhee Kim, Sang Beom Jun, Seungjun Lee, Hyang Woon Lee

P114 Computational model to replicate seizure suppression effect by electrical stimulation

Sora Ahn, Sumin Jo, Eunji Jun, Suin Yu, Hyang Woon Lee, Sang Beom Jun, Seungjun Lee

P115 Identifying excitatory and inhibitory synapses in neuronal networks from spike trains using sorted local transfer entropy

Felix Goetze, Pik-Yin Lai

P116 Neural network model for obstacle avoidance based on neuromorphic computational model of boundary vector cell and head direction cell

Seonghyun Kim, Jeehyun Kwag

P117 Dynamic gating of spike pattern propagation by Hebbian and anti-Hebbian spike timing-dependent plasticity in excitatory feedforward network model

Hyun Jae Jang, Jeehyun Kwag

P118 Inferring characteristics of input correlations of cells exhibiting up-down state transitions in the rat striatum

Marko Filipović, Ramon Reig, Ad Aertsen, Gilad Silberberg, Arvind Kumar

P119 Graph properties of the functional connected brain under the influence of Alzheimer’s disease

Claudia Bachmann, Simone Buttler, Heidi Jacobs, Kim Dillen, Gereon R. Fink, Juraj Kukolja, Abigail Morrison

P120 Learning sparse representations in the olfactory bulb

Daniel Kepple, Hamza Giaffar, Dima Rinberg, Steven Shea, Alex Koulakov

P121 Functional classification of homologous basal-ganglia networks

Jyotika Bahuguna,Tom Tetzlaff, Abigail Morrison, Arvind Kumar, Jeanette Hellgren Kotaleski

P122 Short term memory based on multistability

Tim Kunze, Andre Peterson, Thomas Knösche

P123 A physiologically plausible, computationally efficient model and simulation software for mammalian motor units

Minjung Kim, Hojeong Kim

P125 Decoding laser-induced somatosensory information from EEG

Ji Sung Park, Ji Won Yeon, Sung-Phil Kim

P126 Phase synchronization of alpha activity for EEG-based personal authentication

Jae-Hwan Kang, Chungho Lee, Sung-Phil Kim

P129 Investigating phase-lags in sEEG data using spatially distributed time delays in a large-scale brain network model

Andreas Spiegler, Spase Petkoski, Matias J. Palva, Viktor K. Jirsa

P130 Epileptic seizures in the unfolding of a codimension-3 singularity

Maria L. Saggio, Silvan F. Siep, Andreas Spiegler, William C. Stacey, Christophe Bernard, Viktor K. Jirsa

P131 Incremental dimensional exploratory reasoning under multi-dimensional environment

Oh-hyeon Choung, Yong Jeong

P132 A low-cost model of eye movements and memory in personal visual cognition

Yong-il Lee, Jaeseung Jeong

P133 Complex network analysis of structural connectome of autism spectrum disorder patients

Su Hyun Kim, Mir Jeong, Jaeseung Jeong

P134 Cognitive motives and the neural correlates underlying human social information transmission, gossip

Jeungmin Lee, Jaehyung Kwon, Jerald D. Kralik, Jaeseung Jeong

P135 EEG hyperscanning detects neural oscillation for the social interaction during the economic decision-making

Jaehwan Jahng, Dong-Uk Hwang, Jaeseung Jeong

P136 Detecting purchase decision based on hyperfrontality of the EEG

Jae-Hyung Kwon, Sang-Min Park, Jaeseung Jeong

P137 Vulnerability-based critical neurons, synapses, and pathways in the Caenorhabditis elegans connectome

Seongkyun Kim, Hyoungkyu Kim, Jerald D. Kralik, Jaeseung Jeong

P138 Motif analysis reveals functionally asymmetrical neurons in C. elegans

Pyeong Soo Kim, Seongkyun Kim, Hyoungkyu Kim, Jaeseung Jeong

P139 Computational approach to preference-based serial decision dynamics: do temporal discounting and working memory affect it?

Sangsup Yoon, Jaehyung Kwon, Sewoong Lim, Jaeseung Jeong

P141 Social stress induced neural network reconfiguration affects decision making and learning in zebrafish

Choongseok Park, Thomas Miller, Katie Clements, Sungwoo Ahn, Eoon Hye Ji, Fadi A. Issa

P142 Descriptive, generative, and hybrid approaches for neural connectivity inference from neural activity data

JeongHun Baek, Shigeyuki Oba, Junichiro Yoshimoto, Kenji Doya, Shin Ishii

P145 Divergent-convergent synaptic connectivities accelerate coding in multilayered sensory systems

Thiago S. Mosqueiro, Martin F. Strube-Bloss, Brian Smith, Ramon Huerta

P146 Swinging networks

Michal Hadrava, Jaroslav Hlinka

P147 Inferring dynamically relevant motifs from oscillatory stimuli: challenges, pitfalls, and solutions

Hannah Bos, Moritz Helias

P148 Spatiotemporal mapping of brain network dynamics during cognitive tasks using magnetoencephalography and deep learning

Charles M. Welzig, Zachary J. Harper

P149 Multiscale complexity analysis for the segmentation of MRI images

Won Sup Kim, In-Seob Shin, Hyeon-Man Baek, Seung Kee Han

P150 A neuro-computational model of emotional attention

René Richter, Julien Vitay, Frederick Beuth, Fred H. Hamker

P151 Multi-site delayed feedback stimulation in parkinsonian networks

Kelly Toppin, Yixin Guo

P152 Bistability in Hodgkin–Huxley-type equations

Tatiana Kameneva, Hamish Meffin, Anthony N. Burkitt, David B. Grayden

P153 Phase changes in postsynaptic spiking due to synaptic connectivity and short term plasticity: mathematical analysis of frequency dependency

Mark D. McDonnell, Bruce P. Graham

P154 Quantifying resilience patterns in brain networks: the importance of directionality

Penelope J. Kale, Leonardo L. Gollo

P155 Dynamics of rate-model networks with separate excitatory and inhibitory populations

Merav Stern, L. F. Abbott

P156 A model for multi-stable dynamics in action recognition modulated by integration of silhouette and shading cues

Leonid A. Fedorov, Martin A. Giese

P157 Spiking model for the interaction between action recognition and action execution

Mohammad Hovaidi Ardestani, Martin Giese

P158 Surprise-modulated belief update: how to learn within changing environments?

Mohammad Javad Faraji, Kerstin Preuschoff, Wulfram Gerstner

P159 A fast, stochastic and adaptive model of auditory nerve responses to cochlear implant stimulation

Margriet J. van Gendt, Jeroen J. Briaire, Randy K. Kalkman, Johan H. M. Frijns

P160 Quantitative comparison of graph theoretical measures of simulated and empirical functional brain networks

Won Hee Lee, Sophia Frangou

P161 Determining discriminative properties of fMRI signals in schizophrenia using highly comparative time-series analysis

Ben D. Fulcher, Patricia H. P. Tran, Alex Fornito

P162 Emergence of narrowband LFP oscillations from completely asynchronous activity during seizures and high-frequency oscillations

Stephen V. Gliske, William C. Stacey, Eugene Lim, Katherine A. Holman, Christian G. Fink

P163 Neuronal diversity in structure and function: cross-validation of anatomical and physiological classification of retinal ganglion cells in the mouse

Jinseop S. Kim, Shang Mu, Kevin L. Briggman, H. Sebastian Seung, the EyeWirers

P164 Analysis and modelling of transient firing rate changes in area MT in response to rapid stimulus feature changes

Detlef Wegener, Lisa Bohnenkamp, Udo A. Ernst

P165 Step-wise model fitting accounting for high-resolution spatial measurements: construction of a layer V pyramidal cell model with reduced morphology

Tuomo Mäki-Marttunen, Geir Halnes, Anna Devor, Christoph Metzner, Anders M. Dale, Ole A. Andreassen, Gaute T. Einevoll

P166 Contributions of schizophrenia-associated genes to neuron firing and cardiac pacemaking: a polygenic modeling approach

Tuomo Mäki-Marttunen, Glenn T. Lines, Andy Edwards, Aslak Tveito, Anders M. Dale, Gaute T. Einevoll, Ole A. Andreassen

P167 Local field potentials in a 4 × 4 mm² multi-layered network model

Espen Hagen, Johanna Senk, Sacha J. van Albada, Markus Diesmann

P168 A spiking network model explains multi-scale properties of cortical dynamics

Maximilian Schmidt, Rembrandt Bakker, Kelly Shen, Gleb Bezgin, Claus-Christian Hilgetag, Markus Diesmann, Sacha Jennifer van Albada

P169 Using joint weight-delay spike-timing dependent plasticity to find polychronous neuronal groups

Haoqi Sun, Olga Sourina, Guang-Bin Huang, Felix Klanner, Cornelia Denk

P170 Tensor decomposition reveals RSNs in simulated resting state fMRI

Katharina Glomb, Adrián Ponce-Alvarez, Matthieu Gilson, Petra Ritter, Gustavo Deco

P171 Getting in the groove: testing a new model-based method for comparing task-evoked vs resting-state activity in fMRI data on music listening

Matthieu Gilson, Maria AG Witek, Eric F. Clarke, Mads Hansen, Mikkel Wallentin, Gustavo Deco, Morten L. Kringelbach^, Peter Vuust

P172 STochastic engine for pathway simulation (STEPS) on massively parallel processors

Guido Klingbeil, Erik De Schutter

P173 Toolkit support for complex parallel spatial stochastic reaction–diffusion simulation in STEPS

Weiliang Chen, Erik De Schutter

P174 Modeling the generation and propagation of Purkinje cell dendritic spikes caused by parallel fiber synaptic input

Yunliang Zang, Erik De Schutter

P175 Dendritic morphology determines how dendrites are organized into functional subunits

Sungho Hong, Akira Takashima, Erik De Schutter

P176 A model of Ca²⁺/calmodulin-dependent protein kinase II activity in long term depression at Purkinje cells

Criseida Zamora, Andrew R. Gallimore, Erik De Schutter

P177 Reward-modulated learning of population-encoded vectors for insect-like navigation in embodied agents

Dennis Goldschmidt, Poramate Manoonpong, Sakyasingha Dasgupta

P178 Data-driven neural models part II: connectivity patterns of human seizures

Philippa J. Karoly, Dean R. Freestone, Daniel Soundry, Levin Kuhlmann, Liam Paninski, Mark Cook

P179 Data-driven neural models part I: state and parameter estimation

Dean R. Freestone, Philippa J. Karoly, Daniel Soundry, Levin Kuhlmann, Mark Cook

P180 Spectral and spatial information processing in human auditory streaming

Jaejin Lee, Yonatan I. Fishman, Yale E. Cohen

P181 A tuning curve for the global effects of local perturbations in neural activity: Mapping the systems-level susceptibility of the brain

Leonardo L. Gollo, James A. Roberts, Luca Cocchi

P182 Diverse homeostatic responses to visual deprivation mediated by neural ensembles

Yann Sweeney, Claudia Clopath

P183 Opto-EEG: a novel method for investigating functional connectome in mouse brain based on optogenetics and high density electroencephalography

Soohyun Lee, Woo-Sung Jung, Jee Hyun Choi

P184 Biphasic responses of frontal gamma network to repetitive sleep deprivation during REM sleep

Bowon Kim, Youngsoo Kim, Eunjin Hwang, Jee Hyun Choi

P185 Brain-state correlate and cortical connectivity for frontal gamma oscillations in top-down fashion assessed by auditory steady-state response

Younginha Jung, Eunjin Hwang, Yoon-Kyu Song, Jee Hyun Choi

P186 Neural field model of localized orientation selective activation in V1

James Rankin, Frédéric Chavane

P187 An oscillatory network model of Head direction and Grid cells using locomotor inputs

Karthik Soman, Vignesh Muralidharan, V. Srinivasa Chakravarthy

P188 A computational model of hippocampus inspired by the functional architecture of basal ganglia

Karthik Soman, Vignesh Muralidharan, V. Srinivasa Chakravarthy

P189 A computational architecture to model the microanatomy of the striatum and its functional properties

Sabyasachi Shivkumar, Vignesh Muralidharan, V. Srinivasa Chakravarthy

P190 A scalable cortico-basal ganglia model to understand the neural dynamics of targeted reaching

Vignesh Muralidharan, Alekhya Mandali, B. Pragathi Priyadharsini, Hima Mehta, V. Srinivasa Chakravarthy

P191 Emergence of radial orientation selectivity from synaptic plasticity

Catherine E. Davey, David B. Grayden, Anthony N. Burkitt

P192 How do hidden units shape effective connections between neurons?

Braden A. W. Brinkman, Tyler Kekona, Fred Rieke, Eric Shea-Brown, Michael Buice

P193 Characterization of neural firing in the presence of astrocyte-synapse signaling

Maurizio De Pittà, Hugues Berry, Nicolas Brunel

P194 Metastability of spatiotemporal patterns in a large-scale network model of brain dynamics

James A. Roberts, Leonardo L. Gollo, Michael Breakspear

P195 Comparison of three methods to quantify detection and discrimination capacity estimated from neural population recordings

Gary Marsat, Jordan Drew, Phillip D. Chapman, Kevin C. Daly, Samual P. Bradley

P196 Quantifying the constraints for independent evoked and spontaneous NMDA receptor mediated synaptic transmission at individual synapses

Sat Byul Seo, Jianzhong Su, Ege T. Kavalali, Justin Blackwell

P199 Gamma oscillation via adaptive exponential integrate-and-fire neurons

LieJune Shiau, Laure Buhry, Kanishka Basnayake

P200 Visual face representations during memory retrieval compared to perception

Sue-Hyun Lee, Brandon A. Levy, Chris I. Baker

P201 Top-down modulation of sequential activity within packets modeled using avalanche dynamics

Timothée Leleu, Kazuyuki Aihara

Q28 An auto-encoder network realizes sparse features under the influence of desynchronized vascular dynamics

Ryan T. Philips, Karishma Chhabria, V. Srinivasa Chakravarthy

Neuroprosthetic Decoder Training as Imitation Learning

Neuroprosthetic brain-computer interfaces function via an algorithm which decodes neural activity of the user into movements of an end effector, such as a cursor or robotic arm. In practice, the decoder is often learned by updating its parameters while the user performs a task. When the user's intention is not directly observable, recent methods have demonstrated value in training the decoder against a surrogate for the user's intended movement. Here we show that training a decoder in this way is a novel variant of an imitation learning problem, where an oracle or expert is employed for supervised training in lieu of direct observations, which are not available. Specifically, we describe how a generic imitation learning meta-algorithm, dataset aggregation (DAgger), can be adapted to train a generic brain-computer interface. By deriving existing learning algorithms for brain-computer interfaces in this framework, we provide a novel analysis of regret (an important metric of learning efficacy) for brain-computer interfaces. This analysis allows us to characterize the space of algorithmic variants and bounds on their regret rates. Existing approaches for decoder learning have been performed in the cursor control setting, but the available design principles for these decoders are such that it has been impossible to scale them to naturalistic settings. Leveraging our findings, we then offer an algorithm that combines imitation learning with optimal control, which should allow for training of arbitrary effectors for which optimal control can generate goal-oriented control. We demonstrate this novel and general BCI algorithm with simulated neuroprosthetic control of a 26 degree-of-freedom model of an arm, a sophisticated and realistic end effector.

Bayesian Sparse Regression Analysis Documents the Diversity of Spinal Inhibitory Interneurons

Documenting the extent of cellular diversity is a critical step in defining the functional organization of tissues and organs. To infer cell-type diversity from partial or incomplete transcription factor expression data, we devised a sparse Bayesian framework that is able to handle estimation uncertainty and can incorporate diverse cellular characteristics to optimize experimental design. Focusing on spinal V1 inhibitory interneurons, for which the spatial expression of 19 transcription factors has been mapped, we infer the existence of ~50 candidate V1 neuronal types, many of which localize in compact spatial domains in the ventral spinal cord. We have validated the existence of inferred cell types by direct experimental measurement, establishing this Bayesian framework as an effective platform for cell-type characterization in the nervous system and elsewhere.

Mapping nonlinear receptive field structure in primate retina at single cone resolution

The function of a neural circuit is shaped by the computations performed by its interneurons, which in many cases are not easily accessible to experimental investigation. Here, we elucidate the transformation of visual signals flowing from the input to the output of the primate retina, using a combination of large-scale multi-electrode recordings from an identified ganglion cell type, visual stimulation targeted at individual cone photoreceptors, and a hierarchical computational model. The results reveal nonlinear subunits in the circuity of OFF midget ganglion cells, which subserve high-resolution vision. The model explains light responses to a variety of stimuli more accurately than a linear model, including stimuli targeted to cones within and across subunits. The recovered model components are consistent with known anatomical organization of midget bipolar interneurons. These results reveal the spatial structure of linear and nonlinear encoding, at the resolution of single cells and at the scale of complete circuits.

DOI:http://dx.doi.org/10.7554/eLife.05241.001

Light that enters the eye begins the process of vision by activating two types of photoreceptors: rods, which support vision under low light levels, and cones, which are responsible for fine detail and color vision. Activation of either type of photoreceptor triggers responses in bipolar cells, which activate the ganglion cells that transmit visual signals to the brain.

Bipolar cells therefore belong to a class of cells called interneurons, which relay information from certain cell types to others. Interneurons play an important role in information processing throughout the brain, but directly accessing them or characterizing their role in neural computation is often difficult.

To address this problem, Freeman, Field et al. have developed a combined computational and experimental approach to describe the flow of sensory signals through the circuits within the retina of primates. Large arrays of electrodes were used to record the responses of many retinal ganglion cells in response to the activation or de-activation of pairs of cones. These experiments revealed that the responses of ganglion cells are not simply the sum of the inputs that they receive from cones; specifically, the activation of one cone is not cancelled by the deactivation of another. Instead, the data suggest that bipolar cells add cone inputs together and then pass on the total activation (but not deactivation) to ganglion cells.

By analyzing the responses of ganglion cells to numerous random patterns of cone activation, Freeman, Field et al. were able to estimate the locations and arrangements of bipolar cells that connect to them. These predicted patterns of connectivity agreed with observations from anatomical studies. This work provides detailed insights into how the primate retina works. It also suggests that similar approaches may be used to characterize how signals flow across other brain networks in which large-scale recordings are now possible.

DOI:http://dx.doi.org/10.7554/eLife.05241.002

Efficient "Shotgun" Inference of Neural Connectivity from Highly Sub-sampled Activity Data

Inferring connectivity in neuronal networks remains a key challenge in statistical neuroscience. The “common input” problem presents a major roadblock: it is difficult to reliably distinguish causal connections between pairs of observed neurons versus correlations induced by common input from unobserved neurons. Available techniques allow us to simultaneously record, with sufficient temporal resolution, only a small fraction of the network. Consequently, naive connectivity estimators that neglect these common input effects are highly biased. This work proposes a “shotgun” experimental design, in which we observe multiple sub-networks briefly, in a serial manner. Thus, while the full network cannot be observed simultaneously at any given time, we may be able to observe much larger subsets of the network over the course of the entire experiment, thus ameliorating the common input problem. Using a generalized linear model for a spiking recurrent neural network, we develop a scalable approximate expected loglikelihood-based Bayesian method to perform network inference given this type of data, in which only a small fraction of the network is observed in each time bin. We demonstrate in simulation that the shotgun experimental design can eliminate the biases induced by common input effects. Networks with thousands of neurons, in which only a small fraction of the neurons is observed in each time bin, can be quickly and accurately estimated, achieving orders of magnitude speed up over previous approaches.

Optical imaging of the activity in a neuronal network is limited by the scanning speed of the imaging device. Therefore, typically, only a small fixed part of the network is observed during the entire experiment. However, in such an experiment, it can be hard to infer from the observed activity patterns whether (1) a neuron A directly affects neuron B, or (2) another, unobserved neuron C affects both A and B. To deal with this issue, we propose a “shotgun” observation scheme, in which, at each time point, we observe a small changing subset of the neurons from the network. Consequently, many fewer neurons remain completely unobserved during the entire experiment, enabling us to eventually distinguish between cases (1) and (2) given sufficiently long experiments. Since previous inference algorithms cannot efficiently handle so many missing observations, we develop a scalable algorithm for data acquired using the shotgun observation scheme, in which only a small fraction of the neurons are observed in each time bin. Using this kind of simulated data, we show the algorithm is able to quickly infer connectivity in spiking recurrent networks with thousands of neurons.

Primacy of Flexor Locomotor Pattern Revealed by Ancestral Reversion of Motor Neuron Identity

Spinal circuits can generate locomotor output in the absence of sensory or descending input, but the principles of locomotor circuit organization remain unclear. We sought insight into these principles by considering the elaboration of locomotor circuits across evolution. The identity of limb-innervating motor neurons was reverted to a state resembling that of motor neurons that direct undulatory swimming in primitive aquatic vertebrates, permitting assessment of the role of motor neuron identity in determining locomotor pattern. Two-photon imaging was coupled with spike inference to measure locomotor firing in hundreds of motor neurons in isolated mouse spinal cords. In wild-type preparations, we observed sequential recruitment of motor neurons innervating flexor muscles controlling progressively more distal joints. Strikingly, after reversion of motor neuron identity, virtually all firing patterns became distinctly flexor like. Our findings show that motor neuron identity directs locomotor circuit wiring and indicate the evolutionary primacy of flexor pattern generation.

Encoder-decoder optimization for brain-computer interfaces

Neuroprosthetic brain-computer interfaces are systems that decode neural activity into useful control signals for effectors, such as a cursor on a computer screen. It has long been recognized that both the user and decoding system can adapt to increase the accuracy of the end effector. Co-adaptation is the process whereby a user learns to control the system in conjunction with the decoder adapting to learn the user's neural patterns. We provide a mathematical framework for co-adaptation and relate co-adaptation to the joint optimization of the user's control scheme ("encoding model") and the decoding algorithm's parameters. When the assumptions of that framework are respected, co-adaptation cannot yield better performance than that obtainable by an optimal initial choice of fixed decoder, coupled with optimal user learning. For a specific case, we provide numerical methods to obtain such an optimized decoder. We demonstrate our approach in a model brain-computer interface system using an online prosthesis simulator, a simple human-in-the-loop pyschophysics setup which provides a non-invasive simulation of the BCI setting. These experiments support two claims: that users can learn encoders matched to fixed, optimal decoders and that, once learned, our approach yields expected performance advantages.

On quadrature methods for refractory point process likelihoods

Parametric models of the conditional intensity of a point process (e.g., generalized linear models) are popular in statistical neuroscience, as they allow us to characterize the variability in neural responses in terms of stimuli and spiking history. Parameter estimation in these models relies heavily on accurate evaluations of the log likelihood and its derivatives. Classical approaches use a discretized time version of the spiking process, and recent work has exploited the existence of a refractory period (during which the conditional intensity is zero following a spike) to obtain more accurate estimates of the likelihood. In this brief letter, we demonstrate that this method can be improved significantly by applying classical quadrature methods directly to the resulting continuous-time integral.