Suppression without inhibition: how retinal computation contributes to saccadic suppression.

Visual perception remains stable across saccadic eye movements, despite the concurrent strongly disruptive visual flow. This stability is partially associated with a reduction in visual sensitivity, known as saccadic suppression, which already starts in the retina with reduced ganglion cell sensitivity. However, the retinal circuit mechanisms giving rise to such suppression remain unknown. Here, we describe these mechanisms using electrophysiology in mouse, pig, and macaque retina, 2-photon calcium imaging, computational modeling, and human psychophysics. We find that sequential stimuli, like those that naturally occur during saccades, trigger three independent suppressive mechanisms in the retina. The main mechanism is triggered by contrast-reversing sequential stimuli and originates within the receptive field center of ganglion cells. It does not involve inhibition or other known suppressive mechanisms like saturation or adaptation. Instead, it relies on temporal filtering of the inherently slow response of cone photoreceptors coupled with downstream nonlinearities. Two further mechanisms of suppression are present predominantly in ON ganglion cells and originate in the receptive field surround, highlighting another disparity between ON and OFF ganglion cells. The mechanisms uncovered here likely play a role in shaping the retinal output following eye movements and other natural viewing conditions where sequential stimulation is ubiquitous.

Perceptual saccadic suppression starts in the retina.

Visual sensitivity, probed through perceptual detectability of very brief visual stimuli, is strongly impaired around the time of rapid eye movements. This robust perceptual phenomenon, called saccadic suppression, is frequently attributed to active suppressive signals that are directly derived from eye movement commands. Here we show instead that visual-only mechanisms, activated by saccade-induced image shifts, can account for all perceptual properties of saccadic suppression that we have investigated. Such mechanisms start at, but are not necessarily exclusive to, the very first stage of visual processing in the brain, the retina. Critically, neural suppression originating in the retina outlasts perceptual suppression around the time of saccades, suggesting that extra-retinal movement-related signals, rather than causing suppression, may instead act to shorten it. Our results demonstrate a far-reaching contribution of visual processing mechanisms to perceptual saccadic suppression, starting in the retina, without the need to invoke explicit motor-based suppression commands.

Massively parallel microwire arrays integrated with CMOS chips for neural recording.

Multi-channel electrical recordings of neural activity in the brain is an increasingly powerful method revealing new aspects of neural communication, computation, and prosthetics. However, while planar silicon-based CMOS devices in conventional electronics scale rapidly, neural interface devices have not kept pace. Here, we present a new strategy to interface silicon-based chips with three-dimensional microwire arrays, providing the link between rapidly-developing electronics and high density neural interfaces. The system consists of a bundle of microwires mated to large-scale microelectrode arrays, such as camera chips. This system has excellent recording performance, demonstrated via single unit and local-field potential recordings in isolated retina and in the motor cortex or striatum of awake moving mice. The modular design enables a variety of microwire types and sizes to be integrated with different types of pixel arrays, connecting the rapid progress of commercial multiplexing, digitisation and data acquisition hardware together with a three-dimensional neural interface.

Medullary Respiratory Circuit Is Reorganized by a Seasonally-Induced Program in Preparation for Hibernation.

Deep hibernators go through several cycles of profound drops in body temperature during the winter season, with core temperatures sometimes reaching near freezing. Yet unlike non-hibernating mammals, they can sustain breathing rhythms. The physiological processes that make this possible are still not understood. In this study, we focused on the medullary Ventral Respiratory Column of a facultative hibernator, the Syrian hamster. Using shortened day-lengths, we induced a "winter-adapted" physiological state, which is a prerequisite for hibernation. When recording electrophysiological signals from acute slices in the winter-adapted pre-Bötzinger complex (preBötC), spike trains showed higher spike rates, amplitudes, complexity, as well as higher temperature sensitivity, suggesting an increase in connectivity and/or synaptic strength during the winter season. We further examined action potential waveforms and found that the depolarization integral, as measured by the area under the curve, is selectively enhanced in winter-adapted animals. This suggests that a shift in the ion handling kinetics is also being induced by the winter-adaptation program. RNA sequencing of respiratory pre-motor neurons, followed by gene set enrichment analysis, revealed differential regulation and splicing in structural, synaptic, and ion handling genes. Splice junction analysis suggested that differential exon usage is occurring in a select subset of ion handling subunits (ATP1A3, KCNC3, SCN1B), and synaptic structure genes (SNCB, SNCG, RAB3A). Our findings show that the hamster respiratory center undergoes a seasonally-cued alteration in electrophysiological properties, likely protecting against respiratory failure at low temperatures.

Optimal Electrode Size for Multi-Scale Extracellular-Potential Recording From Neuronal Assemblies.

Advances in microfabrication technology have enabled the production of devices containing arrays of thousands of closely spaced recording electrodes, which afford subcellular resolution of electrical signals in neurons and neuronal networks. Rationalizing the electrode size and configuration in such arrays demands consideration of application-specific requirements and inherent features of the electrodes. Tradeoffs among size, spatial density, sensitivity, noise, attenuation, and other factors are inevitable. Although recording extracellular signals from neurons with planar metal electrodes is fairly well established, the effects of the electrode characteristics on the quality and utility of recorded signals, especially for small, densely packed electrodes, have yet to be fully characterized. Here, we present a combined experimental and computational approach to elucidating how electrode size, and size-dependent parameters, such as impedance, baseline noise, and transmission characteristics, influence recorded neuronal signals. Using arrays containing platinum electrodes of different sizes, we experimentally evaluated the electrode performance in the recording of local field potentials (LFPs) and extracellular action potentials (EAPs) from the following cell preparations: acute brain slices, dissociated cell cultures, and organotypic slice cultures. Moreover, we simulated the potential spatial decay of point-current sources to investigate signal averaging using known signal sources. We demonstrated that the noise and signal attenuation depend more on the electrode impedance than on electrode size, per se, especially for electrodes <10 µm in width or diameter to achieve high-spatial-resolution readout. By minimizing electrode impedance of small electrodes (<10 µm) via surface modification, we could maximize the signal-to-noise ratio to electrically visualize the propagation of axonal EAPs and to isolate single-unit spikes. Due to the large amplitude of LFP signals, recording quality was high and nearly independent of electrode size. These findings should be of value in configuring in vitro and in vivo microelectrode arrays for extracellular recordings with high spatial resolution in various applications.

Technologies to Study Action Potential Propagation With a Focus on HD-MEAs.

Axons convey information in neuronal circuits via reliable conduction of action potentials (APs) from the axon initial segment (AIS) to the presynaptic terminals. Recent experimental findings increasingly evidence that the axonal function is not limited to the simple transmission of APs. Advances in subcellular-resolution recording techniques have shown that axons display activity-dependent modulation in spike shape and conduction velocity, which influence synaptic strength and latency. We briefly review here, how recent methodological developments facilitate the understanding of the axon physiology. We included the three most common methods, i.e., genetically encoded voltage imaging (GEVI), subcellular patch-clamp and high-density microelectrode arrays (HD-MEAs). We then describe the potential of using HD-MEAs in studying axonal physiology in more detail. Due to their robustness, amenability to high-throughput and high spatiotemporal resolution, HD-MEAs can provide a direct functional electrical readout of single cells and cellular ensembles at subcellular resolution. HD-MEAs can, therefore, be employed in investigating axonal pathologies, the effects of large-scale genomic interventions (e.g., with RNAi or CRISPR) or in compound screenings. A combination of extracellular microelectrode arrays (MEAs), intracellular microelectrodes and optical imaging may potentially reveal yet unexplored repertoires of axonal functions.

Automatic spike sorting for high-density microelectrode arrays.

High-density microelectrode arrays can be used to record extracellular action potentials from hundreds to thousands of neurons simultaneously. Efficient spike sorters must be developed to cope with such large data volumes. Most existing spike sorting methods for single electrodes or small multielectrodes, however, suffer from the "curse of dimensionality" and cannot be directly applied to recordings with hundreds of electrodes. This holds particularly true for the standard reference spike sorting algorithm, principal component analysis-based feature extraction, followed by k-means or expectation maximization clustering, against which most spike sorters are evaluated. We present a spike sorting algorithm that circumvents the dimensionality problem by sorting local groups of electrodes independently with classical spike sorting approaches. It is scalable to any number of recording electrodes and well suited for parallel computing. The combination of data prewhitening before the principal component analysis-based extraction and a parameter-free clustering algorithm obviated the need for parameter adjustments. We evaluated its performance using surrogate data in which we systematically varied spike amplitudes and spike rates and that were generated by inserting template spikes into the voltage traces of real recordings. In a direct comparison, our algorithm could compete with existing state-of-the-art spike sorters in terms of sensitivity and precision, while parameter adjustment or manual cluster curation was not required. NEW & NOTEWORTHY We present an automatic spike sorting algorithm that combines three strategies to scale classical spike sorting techniques for high-density microelectrode arrays: 1) splitting the recording electrodes into small groups and sorting them independently; 2) clustering a subset of spikes and classifying the rest to limit computation time; and 3) prewhitening the spike waveforms to enable the use of parameter-free clustering. Finally, we combined these strategies into an automatic spike sorter that is competitive with state-of-the-art spike sorters.

How Diverse Retinal Functions Arise from Feedback at the First Visual Synapse.

Many brain regions contain local interneurons of distinct types. How does an interneuron type contribute to the input-output transformations of a given brain region? We addressed this question in the mouse retina by chemogenetically perturbing horizontal cells, an interneuron type providing feedback at the first visual synapse, while monitoring the light-driven spiking activity in thousands of ganglion cells, the retinal output neurons. We uncovered six reversible perturbation-induced effects in the response dynamics and response range of ganglion cells. The effects were enhancing or suppressive, occurred in different response epochs, and depended on the ganglion cell type. A computational model of the retinal circuitry reproduced all perturbation-induced effects and led us to assign specific functions to horizontal cells with respect to different ganglion cell types. Our combined experimental and theoretical work reveals how a single interneuron type can differentially shape the dynamical properties of distinct output channels of a brain region.

Acquisition of Bioelectrical Signals with Small Electrodes.

Although the mechanisms of recording bioelectrical signals from different types of electrogenic cells (neurons, cardiac cells etc.) by means of planar metal electrodes have been extensively studied, the recording characteristics and conditions for very small electrode sizes are not yet established. Here, we present a combined experimental and computational approach to elucidate, how the electrode size influences the recorded signals, and how inherent properties of the electrode, such as impedance, noise, and transmission characteristics shape the signal. We demonstrate that good quality recordings can be achieved with electrode diameters of less than 10 µm, provided that impedance reduction measures have been implemented and provided that a set of requirements for signal amplification has been met.

Tracking individual action potentials throughout mammalian axonal arbors.

Axons are neuronal processes specialized for conduction of action potentials (APs). The timing and temporal precision of APs when they reach each of the synapses are fundamentally important for information processing in the brain. Due to small diameters of axons, direct recording of single AP transmission is challenging. Consequently, most knowledge about axonal conductance derives from modeling studies or indirect measurements. We demonstrate a method to noninvasively and directly record individual APs propagating along millimeter-length axonal arbors in cortical cultures with hundreds of microelectrodes at microsecond temporal resolution. We find that cortical axons conduct single APs with high temporal precision (~100 µs arrival time jitter per mm length) and reliability: in more than 8,000,000 recorded APs, we did not observe any conduction or branch-point failures. Upon high-frequency stimulation at 100 Hz, successive became slower, and their arrival time precision decreased by 20% and 12% for the 100th AP, respectively.

Combination of High-density Microelectrode Array and Patch Clamp Recordings to Enable Studies of Multisynaptic Integration.

We present a novel, all-electric approach to record and to precisely control the activity of tens of individual presynaptic neurons. The method allows for parallel mapping of the efficacy of multiple synapses and of the resulting dynamics of postsynaptic neurons in a cortical culture. For the measurements, we combine an extracellular high-density microelectrode array, featuring 11'000 electrodes for extracellular recording and stimulation, with intracellular patch-clamp recording. We are able to identify the contributions of individual presynaptic neurons - including inhibitory and excitatory synaptic inputs - to postsynaptic potentials, which enables us to study dendritic integration. Since the electrical stimuli can be controlled at microsecond resolution, our method enables to evoke action potentials at tens of presynaptic cells in precisely orchestrated sequences of high reliability and minimum jitter. We demonstrate the potential of this method by evoking short- and long-term synaptic plasticity through manipulation of multiple synaptic inputs to a specific neuron.

Structures of Neural Correlation and How They Favor Coding.

The neural representation of information suffers from "noise"-the trial-to-trial variability in the response of neurons. The impact of correlated noise upon population coding has been debated, but a direct connection between theory and experiment remains tenuous. Here, we substantiate this connection and propose a refined theoretical picture. Using simultaneous recordings from a population of direction-selective retinal ganglion cells, we demonstrate that coding benefits from noise correlations. The effect is appreciable already in small populations, yet it is a collective phenomenon. Furthermore, the stimulus-dependent structure of correlation is key. We develop simple functional models that capture the stimulus-dependent statistics. We then use them to quantify the performance of population coding, which depends upon interplays of feature sensitivities and noise correlations in the population. Because favorable structures of correlation emerge robustly in circuits with noisy, nonlinear elements, they will arise and benefit coding beyond the confines of retina.

Congenital Nystagmus Gene FRMD7 Is Necessary for Establishing a Neuronal Circuit Asymmetry for Direction Selectivity.

Neuronal circuit asymmetries are important components of brain circuits, but the molecular pathways leading to their establishment remain unknown. Here we found that the mutation of FRMD7, a gene that is defective in human congenital nystagmus, leads to the selective loss of the horizontal optokinetic reflex in mice, as it does in humans. This is accompanied by the selective loss of horizontal direction selectivity in retinal ganglion cells and the transition from asymmetric to symmetric inhibitory input to horizontal direction-selective ganglion cells. In wild-type retinas, we found FRMD7 specifically expressed in starburst amacrine cells, the interneuron type that provides asymmetric inhibition to direction-selective retinal ganglion cells. This work identifies FRMD7 as a key regulator in establishing a neuronal circuit asymmetry, and it suggests the involvement of a specific inhibitory neuron type in the pathophysiology of a neurological disease.

Time-lapse electrical impedance spectroscopy for monitoring the cell cycle of single immobilized S. pombe cells.

As a complement and alternative to optical methods, wide-band electrical impedance spectroscopy (EIS) enables multi-parameter, label-free and real-time detection of cellular and subcellular features. We report on a microfluidics-based system designed to reliably capture single rod-shaped Schizosaccharomyces pombe cells by applying suction through orifices in a channel wall. The system enables subsequent culturing of immobilized cells in an upright position, while dynamic changes in cell-cycle state and morphology were continuously monitored through EIS over a broad frequency range. Besides measuring cell growth, clear impedance signals for nuclear division have been obtained. The EIS system has been characterized with respect to sensitivity and detection limits. The spatial resolution in measuring cell length was 0.25 µm, which corresponds to approximately a 5-min interval of cell growth under standard conditions. The comprehensive impedance data sets were also used to determine the occurrence of nuclear division and cytokinesis. The obtained results have been validated through concurrent confocal imaging and plausibilized through comparison with finite-element modeling data. The possibility to monitor cellular and intracellular features of single S. pombe cells during the cell cycle at high spatiotemporal resolution renders the presented microfluidics-based EIS system a suitable tool for dynamic single-cell investigations.

A method for electrophysiological characterization of hamster retinal ganglion cells using a high-density CMOS microelectrode array.

Knowledge of neuronal cell types in the mammalian retina is important for the understanding of human retinal disease and the advancement of sight-restoring technology, such as retinal prosthetic devices. A somewhat less utilized animal model for retinal research is the hamster, which has a visual system that is characterized by an area centralis and a wide visual field with a broad binocular component. The hamster retina is optimally suited for recording on the microelectrode array (MEA), because it intrinsically lies flat on the MEA surface and yields robust, large-amplitude signals. However, information in the literature about hamster retinal ganglion cell functional types is scarce. The goal of our work is to develop a method featuring a high-density (HD) complementary metal-oxide-semiconductor (CMOS) MEA technology along with a sequence of standardized visual stimuli in order to categorize ganglion cells in isolated Syrian Hamster (Mesocricetus auratus) retina. Since the HD-MEA is capable of recording at a higher spatial resolution than most MEA systems (17.5 µm electrode pitch), we were able to record from a large proportion of RGCs within a selected region. Secondly, we chose our stimuli so that they could be run during the experiment without intervention or computation steps. The visual stimulus set was designed to activate the receptive fields of most ganglion cells in parallel and to incorporate various visual features to which different cell types respond uniquely. Based on the ganglion cell responses, basic cell properties were determined: direction selectivity, speed tuning, width tuning, transience, and latency. These properties were clustered to identify ganglion cell types in the hamster retina. Ultimately, we recorded up to a cell density of 2780 cells/mm(2) at 2 mm (42°) from the optic nerve head. Using five parameters extracted from the responses to visual stimuli, we obtained seven ganglion cell types.

Spike sorting of synchronous spikes from local neuron ensembles.

Synchronous spike discharge of cortical neurons is thought to be a fingerprint of neuronal cooperativity. Because neighboring neurons are more densely connected to one another than neurons that are located further apart, near-synchronous spike discharge can be expected to be prevalent and it might provide an important basis for cortical computations. Using microelectrodes to record local groups of neurons does not allow for the reliable separation of synchronous spikes from different cells, because available spike sorting algorithms cannot correctly resolve the temporally overlapping waveforms. We show that high spike sorting performance of in vivo recordings, including overlapping spikes, can be achieved with a recently developed filter-based template matching procedure. Using tetrodes with a three-dimensional structure, we demonstrate with simulated data and ground truth in vitro data, obtained by dual intracellular recording of two neurons located next to a tetrode, that the spike sorting of synchronous spikes can be as successful as the spike sorting of nonoverlapping spikes and that the spatial information provided by multielectrodes greatly reduces the error rates. We apply the method to tetrode recordings from the prefrontal cortex of behaving primates, and we show that overlapping spikes can be identified and assigned to individual neurons to study synchronous activity in local groups of neurons.

Visual coding with a population of direction-selective neurons.

The brain decodes the visual scene from the action potentials of ∼20 retinal ganglion cell types. Among the retinal ganglion cells, direction-selective ganglion cells (DSGCs) encode motion direction. Several studies have focused on the encoding or decoding of motion direction by recording multiunit activity, mainly in the visual cortex. In this study, we simultaneously recorded from all four types of ON-OFF DSGCs of the rabbit retina using a microelectronics-based high-density microelectrode array (HDMEA) and decoded their concerted activity using probabilistic and linear decoders. Furthermore, we investigated how the modification of stimulus parameters (velocity, size, angle of moving object) and the use of different tuning curve fits influenced decoding precision. Finally, we simulated ON-OFF DSGC activity, based on real data, in order to understand how tuning curve widths and the angular distribution of the cells' preferred directions influence decoding performance. We found that probabilistic decoding strategies outperformed, on average, linear methods and that decoding precision was robust to changes in stimulus parameters such as velocity. The removal of noise correlations among cells, by random shuffling trials, caused a drop in decoding precision. Moreover, we found that tuning curves are broad in order to minimize large errors at the expense of a higher average error, and that the retinal direction-selective system would not substantially benefit, on average, from having more than four types of ON-OFF DSGCs or from a perfect alignment of the cells' preferred directions.

ViSAPy: a Python tool for biophysics-based generation of virtual spiking activity for evaluation of spike-sorting algorithms.

BACKGROUND: New, silicon-based multielectrodes comprising hundreds or more electrode contacts offer the possibility to record spike trains from thousands of neurons simultaneously. This potential cannot be realized unless accurate, reliable automated methods for spike sorting are developed, in turn requiring benchmarking data sets with known ground-truth spike times. NEW METHOD: We here present a general simulation tool for computing benchmarking data for evaluation of spike-sorting algorithms entitled ViSAPy (Virtual Spiking Activity in Python). The tool is based on a well-established biophysical forward-modeling scheme and is implemented as a Python package built on top of the neuronal simulator NEURON and the Python tool LFPy. RESULTS: ViSAPy allows for arbitrary combinations of multicompartmental neuron models and geometries of recording multielectrodes. Three example benchmarking data sets are generated, i.e., tetrode and polytrode data mimicking in vivo cortical recordings and microelectrode array (MEA) recordings of in vitro activity in salamander retinas. The synthesized example benchmarking data mimics salient features of typical experimental recordings, for example, spike waveforms depending on interspike interval. COMPARISON WITH EXISTING METHODS: ViSAPy goes beyond existing methods as it includes biologically realistic model noise, synaptic activation by recurrent spiking networks, finite-sized electrode contacts, and allows for inhomogeneous electrical conductivities. ViSAPy is optimized to allow for generation of long time series of benchmarking data, spanning minutes of biological time, by parallel execution on multi-core computers. CONCLUSION: ViSAPy is an open-ended tool as it can be generalized to produce benchmarking data or arbitrary recording-electrode geometries and with various levels of complexity.

Bayes optimal template matching for spike sorting - combining fisher discriminant analysis with optimal filtering.

Spike sorting, i.e., the separation of the firing activity of different neurons from extracellular measurements, is a crucial but often error-prone step in the analysis of neuronal responses. Usually, three different problems have to be solved: the detection of spikes in the extracellular recordings, the estimation of the number of neurons and their prototypical (template) spike waveforms, and the assignment of individual spikes to those putative neurons. If the template spike waveforms are known, template matching can be used to solve the detection and classification problem. Here, we show that for the colored Gaussian noise case the optimal template matching is given by a form of linear filtering, which can be derived via linear discriminant analysis. This provides a Bayesian interpretation for the well-known matched filter output. Moreover, with this approach it is possible to compute a spike detection threshold analytically. The method can be implemented by a linear filter bank derived from the templates, and can be used for online spike sorting of multielectrode recordings. It may also be applicable to detection and classification problems of transient signals in general. Its application significantly decreases the error rate on two publicly available spike-sorting benchmark data sets in comparison to state-of-the-art template matching procedures. Finally, we explore the possibility to resolve overlapping spikes using the template matching outputs and show that they can be resolved with high accuracy.