Localization of neurons from extracellular footprints.

BACKGROUND: High density microelectrode arrays (HD-MEAs) are now widely used for both in-vitro and in-vivo recordings, as they allow spikes from hundreds of neurons to be recorded simultaneously. Since extracellular recordings do not allow visualization of the recorded neurons, algorithms are needed to estimate their physical positions, especially to track their movements when the are drifting away from recording devices. NEW METHOD: The objective of this study was to evaluate the performance of multiple algorithms for neuron localization solely from extracellular traces (MEA recordings), either artificial or obtained from mouse retina. The algorithms compared included center-of-mass, monopolar, and grid-based algorithms. The first method is a barycenter calculation. The second algorithm infers the position of the cell using triangulation with the assumption that the neuron behaves as a monopole. Finally, grid-based methods rely on comparing the recorded spike with a projection of spikes of hypothetical neurons with different positions. RESULTS: The Grid-Based algorithm yielded the most satisfactory outcomes. The center-of-mass exhibited a minimal computational cost, yet its average localization was suboptimal. Monopolar algorithms gave cell localizations with an average error of less than 10µm, but they had considerable variability and a high computational cost. For the grid-based method, the variability was smaller, with satisfactory performance and low computational cost. COMPARISON WITH EXISTING METHOD(S): The accuracy of the different localization methods benchmarked in this article had not been properly tested with ground-truth recordings before. CONCLUSION: The objective of this article is to provide guidance to researchers on the selection of optimal methods for localizing neurons based on MEA recordings.

Altered Functional Responses of the Retina in B6 Albino Tyrc/c Mice.

Purpose: Mammals with albinism present low visual discrimination ability and different proportions of certain retinal cell subtypes. As the spatial resolution of the retina depends on the visual field sampling by retinal ganglion cells (RGCs) based on the convergence of upstream cell inputs, it could be affected in albinism and thus modify the RGC function. Methods: We used the Tyrc/c line, a mouse model of oculocutaneous albinism type 1 (OCA1), carrying a tyrosinase mutation, and previously characterized by a total absence of pigment and severe visual deficits. To assess their retinal function, we recorded the light responses of hundreds of RGCs ex vivo using multi-electrode array (MEA). We estimated the receptive field (RF)-center diameter of Tyr+/c and Tyrc/c RGCs using a checkerboard stimulation before simultaneously stimulating the center and surround of RGC RFs with full-field flashes. Results: Following checkerboard stimulation, the RF-center diameters of RGCs were indistinguishable between Tyrc/c and Tyr+/c retinas. Nevertheless, RGCs from Tyrc/c retinas presented more OFF responses to full-field flashes than RGCs from Tyr+/c retinas. Unlike Tyr+/c retinas, very few OFF-center RGCs switched polarity to ON or ON-OFF responses after full-field flashes in Tyrc/c retinas, suggesting a different surround suppression in these retinas. Conclusions: The retinal output signal is affected in Tyrc/c retinas, despite intact RF-center diameters of their RGCs. Adaptive mechanisms during development are probably responsible for this change in RGC responses, related to the absence of ocular pigments.

Temporal pattern recognition in retinal ganglion cells is mediated by dynamical inhibitory synapses.

A fundamental task for the brain is to generate predictions of future sensory inputs, and signal errors in these predictions. Many neurons have been shown to signal omitted stimuli during periodic stimulation, even in the retina. However, the mechanisms of this error signaling are unclear. Here we show that depressing inhibitory synapses shape the timing of the response to an omitted stimulus in the retina. While ganglion cells, the retinal output, responded to an omitted flash with a constant latency over many frequencies of the flash sequence, we found that this was not the case once inhibition was blocked. We built a simple circuit model and showed that depressing inhibitory synapses were a necessary component to reproduce our experimental findings. A new prediction of our model is that the accuracy of the constant latency requires a sufficient amount of flashes in the stimulus, which we could confirm experimentally. Depressing inhibitory synapses could thus be a key component to generate the predictive responses observed in the retina, and potentially in many brain areas.

Encoding surprise by retinal ganglion cells.

The efficient coding hypothesis posits that early sensory neurons transmit maximal information about sensory stimuli, given internal constraints. A central prediction of this theory is that neurons should preferentially encode stimuli that are most surprising. Previous studies suggest this may be the case in early visual areas, where many neurons respond strongly to rare or surprising stimuli. For example, previous research showed that when presented with a rhythmic sequence of full-field flashes, many retinal ganglion cells (RGCs) respond strongly at the instance the flash sequence stops, and when another flash would be expected. This phenomenon is called the 'omitted stimulus response'. However, it is not known whether the responses of these cells varies in a graded way depending on the level of stimulus surprise. To investigate this, we presented retinal neurons with extended sequences of stochastic flashes. With this stimulus, the surprise associated with a particular flash/silence, could be quantified analytically, and varied in a graded manner depending on the previous sequences of flashes and silences. Interestingly, we found that RGC responses could be well explained by a simple normative model, which described how they optimally combined their prior expectations and recent stimulus history, so as to encode surprise. Further, much of the diversity in RGC responses could be explained by the model, due to the different prior expectations that different neurons had about the stimulus statistics. These results suggest that even as early as the retina many cells encode surprise, relative to their own, internally generated expectations.

Optogenetic targeting of AII amacrine cells restores retinal computations performed by the inner retina.

Most inherited retinal dystrophies display progressive photoreceptor cell degeneration leading to severe visual impairment. Optogenetic reactivation of inner retinal neurons is a promising avenue to restore vision in retinas having lost their photoreceptors. Expression of optogenetic proteins in surviving ganglion cells, the retinal output, allows them to take on the lost photoreceptive function. Nonetheless, this creates an exclusively ON retina by expression of depolarizing optogenetic proteins in all classes of ganglion cells, whereas a normal retina extracts several features from the visual scene, with different ganglion cells detecting light increase (ON) and light decrease (OFF). Refinement of this therapeutic strategy should thus aim at restoring these computations. Here we used a vector that targets gene expression to a specific interneuron of the retina called the AII amacrine cell. AII amacrine cells simultaneously activate the ON pathway and inhibit the OFF pathway. We show that the optogenetic stimulation of AII amacrine cells allows restoration of both ON and OFF responses in the retina, but also mediates other types of retinal processing such as sustained and transient responses. Targeting amacrine cells with optogenetics is thus a promising avenue to restore better retinal function and visual perception in patients suffering from retinal degeneration.

Small-correlation expansion to quantify information in noisy sensory systems.

Neural networks encode information through their collective spiking activity in response to external stimuli. This population response is noisy and strongly correlated, with a complex interplay between correlations induced by the stimulus, and correlations caused by shared noise. Understanding how these correlations affect information transmission has so far been limited to pairs or small groups of neurons, because the curse of dimensionality impedes the evaluation of mutual information in larger populations. Here, we develop a small-correlation expansion to compute the stimulus information carried by a large population of neurons, yielding interpretable analytical expressions in terms of the neurons' firing rates and pairwise correlations. We validate the approximation on synthetic data and demonstrate its applicability to electrophysiological recordings in the vertebrate retina, allowing us to quantify the effects of noise correlations between neurons and of memory in single neurons.

Differences in nonlinearities determine retinal cell types.

Classifying neurons in different types is still an open challenge. In the retina, recent works have taken advantage of the ability to record from a large number of cells to classify ganglion cells into different types based on functional information. Although the first attempts in this direction used the receptive field properties of each cell to classify them, more recent approaches have proposed to cluster ganglion cells directly based on their response to stimuli. These two approaches have not been compared directly. Here, we recorded the responses of a large number of ganglion cells and compared two methods for classifying them into functional groups, one based on the receptive field properties, and the other one using directly their responses to stimuli with various temporal frequencies. We show that the response-based approach allows separation of more types than the receptive field-based method, leading to a better classification. This better granularity is due to the fact that the response-based method takes into account not only the linear part of ganglion cell function but also some of the nonlinearities. A careful characterization of nonlinear processing is thus key to allowing functional classification of sensory neurons.NEW & NOTEWORTHY In the retina, ganglion cells can be classified based on their response to visual stimuli. Although some methods are based on the modeling of receptive fields, others rely on responses to characteristic stimuli. We compared these two classes of methods and show that the latter provides a higher discrimination performance. We also show that this gain arises from the ability to account for the nonlinear behavior of neurons.

Stimulus invariant aspects of the retinal code drive discriminability of natural scenes.

Everything that the brain sees must first be encoded by the retina, which maintains a reliable representation of the visual world in many different, complex natural scenes while also adapting to stimulus changes. Decomposing the population code into independent and cell-cell interactions reveals how broad scene structure is encoded in the adapted retinal output. By recording from the same retina while presenting many different natural movies, we see that the population structure, characterized by strong interactions, is consistent across both natural and synthetic stimuli. We show that these interactions contribute to encoding scene identity. We also demonstrate that this structure likely arises in part from shared bipolar cell input as well as from gap junctions between retinal ganglion cells and amacrine cells.

Ectopic expression of a mechanosensitive channel confers spatiotemporal resolution to ultrasound stimulations of neurons for visual restoration.

Remote and precisely controlled activation of the brain is a fundamental challenge in the development of brain-machine interfaces for neurological treatments. Low-frequency ultrasound stimulation can be used to modulate neuronal activity deep in the brain, especially after expressing ultrasound-sensitive proteins. But so far, no study has described an ultrasound-mediated activation strategy whose spatiotemporal resolution and acoustic intensity are compatible with the mandatory needs of brain-machine interfaces, particularly for visual restoration. Here we combined the expression of large-conductance mechanosensitive ion channels with uncustomary high-frequency ultrasonic stimulation to activate retinal or cortical neurons over millisecond durations at a spatiotemporal resolution and acoustic energy deposit compatible with vision restoration. The in vivo sonogenetic activation of the visual cortex generated a behaviour associated with light perception. Our findings demonstrate that sonogenetics can deliver millisecond pattern presentations via an approach less invasive than current brain-machine interfaces for visual restoration.

All-optical inter-layers functional connectivity investigation in the mouse retina.

We developed a multi-unit microscope for all-optical inter-layers circuits interrogation. The system performs two-photon (2P) functional imaging and 2P multiplexed holographic optogenetics at axially distinct planes. We demonstrated the capability of the system to map, in the mouse retina, the functional connectivity between rod bipolar cells (RBCs) and ganglion cells (GCs) by activating single or defined groups of RBCs while recording the evoked response in the GC layer with cell-type specificity and single-cell resolution. We then used a logistic model to probe the functional connectivity between cell types by deriving the "cellular receptive field" describing how RBCs impact each GC type. With the capability to simultaneously image and control neuronal activity at axially distinct planes, the system enables a precise interrogation of multi-layered circuits. Understanding this information transfer is a promising avenue to dissect complex neural circuits and understand the neural basis of computations.

Context-dependent selectivity to natural images in the retina.

Retina ganglion cells extract specific features from natural scenes and send this information to the brain. In particular, they respond to local light increase (ON responses), and/or decrease (OFF). However, it is unclear if this ON-OFF selectivity, characterized with synthetic stimuli, is maintained under natural scene stimulation. Here we recorded ganglion cell responses to natural images slightly perturbed by random noise patterns to determine their selectivity during natural stimulation. The ON-OFF selectivity strongly depended on the specific image. A single ganglion cell can signal luminance increase for one image, and luminance decrease for another. Modeling and experiments showed that this resulted from the non-linear combination of different retinal pathways. Despite the versatility of the ON-OFF selectivity, a systematic analysis demonstrated that contrast was reliably encoded in these responses. Our perturbative approach uncovered the selectivity of retinal ganglion cells to more complex features than initially thought.

Learning low-dimensional generalizable natural features from retina using a U-net.

Much of sensory neuroscience focuses on presenting stimuli that are chosen by the experimenter because they are parametric and easy to sample and are thought to be behaviorally relevant to the organism. However, it is not generally known what these relevant features are in complex, natural scenes. This work focuses on using the retinal encoding of natural movies to determine the presumably behaviorally-relevant features that the brain represents. It is prohibitive to parameterize a natural movie and its respective retinal encoding fully. We use time within a natural movie as a proxy for the whole suite of features evolving across the scene. We then use a task-agnostic deep architecture, an encoder-decoder, to model the retinal encoding process and characterize its representation of "time in the natural scene" in a compressed latent space. In our end-to-end training, an encoder learns a compressed latent representation from a large population of salamander retinal ganglion cells responding to natural movies, while a decoder samples from this compressed latent space to generate the appropriate future movie frame. By comparing latent representations of retinal activity from three movies, we find that the retina has a generalizable encoding for time in the natural scene: the precise, low-dimensional representation of time learned from one movie can be used to represent time in a different movie, with up to 17 ms resolution. We then show that static textures and velocity features of a natural movie are synergistic. The retina simultaneously encodes both to establishes a generalizable, low-dimensional representation of time in the natural scene.

Mice Lacking Gpr179 with Complete Congenital Stationary Night Blindness Are a Good Model for Myopia.

Mutations in GPR179 are one of the most common causes of autosomal recessive complete congenital stationary night blindness (cCSNB). This retinal disease is characterized in patients by impaired dim and night vision, associated with other ocular symptoms, including high myopia. cCSNB is caused by a complete loss of signal transmission from photoreceptors to ON-bipolar cells. In this study, we hypothesized that the lack of Gpr179 and the subsequent impaired ON-pathway could lead to myopic features in a mouse model of cCSNB. Using ultra performance liquid chromatography, we show that adult Gpr179-/- mice have a significant decrease in both retinal dopamine and 3,4-dihydroxyphenylacetic acid, compared to Gpr179+/+ mice. This alteration of the dopaminergic system is thought to be correlated with an increased susceptibility to lens-induced myopia but does not affect the natural refractive development. Altogether, our data added a novel myopia model, which could be used to identify therapeutic interventions.

Inferring the function performed by a recurrent neural network.

A central goal in systems neuroscience is to understand the functions performed by neural circuits. Previous top-down models addressed this question by comparing the behaviour of an ideal model circuit, optimised to perform a given function, with neural recordings. However, this requires guessing in advance what function is being performed, which may not be possible for many neural systems. To address this, we propose an inverse reinforcement learning (RL) framework for inferring the function performed by a neural network from data. We assume that the responses of each neuron in a network are optimised so as to drive the network towards 'rewarded' states, that are desirable for performing a given function. We then show how one can use inverse RL to infer the reward function optimised by the network from observing its responses. This inferred reward function can be used to predict how the neural network should adapt its dynamics to perform the same function when the external environment or network structure changes. This could lead to theoretical predictions about how neural network dynamics adapt to deal with cell death and/or varying sensory stimulus statistics.

Predicting synchronous firing of large neural populations from sequential recordings.

A major goal in neuroscience is to understand how populations of neurons code for stimuli or actions. While the number of neurons that can be recorded simultaneously is increasing at a fast pace, in most cases these recordings cannot access a complete population: some neurons that carry relevant information remain unrecorded. In particular, it is hard to simultaneously record all the neurons of the same type in a given area. Recent progress have made possible to profile each recorded neuron in a given area thanks to genetic and physiological tools, and to pool together recordings from neurons of the same type across different experimental sessions. However, it is unclear how to infer the activity of a full population of neurons of the same type from these sequential recordings. Neural networks exhibit collective behaviour, e.g. noise correlations and synchronous activity, that are not directly captured by a conditionally-independent model that would just put together the spike trains from sequential recordings. Here we show that we can infer the activity of a full population of retina ganglion cells from sequential recordings, using a novel method based on copula distributions and maximum entropy modeling. From just the spiking response of each ganglion cell to a repeated stimulus, and a few pairwise recordings, we could predict the noise correlations using copulas, and then the full activity of a large population of ganglion cells of the same type using maximum entropy modeling. Remarkably, we could generalize to predict the population responses to different stimuli with similar light conditions and even to different experiments. We could therefore use our method to construct a very large population merging cells' responses from different experiments. We predicted that synchronous activity in ganglion cell populations saturates only for patches larger than 1.5mm in radius, beyond what is today experimentally accessible.

Optogenetic therapy: high spatiotemporal resolution and pattern discrimination compatible with vision restoration in non-human primates.

Vision restoration is an ideal medical application for optogenetics, because the eye provides direct optical access to the retina for stimulation. Optogenetic therapy could be used for diseases involving photoreceptor degeneration, such as retinitis pigmentosa or age-related macular degeneration. We describe here the selection, in non-human primates, of a specific optogenetic construct currently tested in a clinical trial. We used the microbial opsin ChrimsonR, and showed that the AAV2.7m8 vector had a higher transfection efficiency than AAV2 in retinal ganglion cells (RGCs) and that ChrimsonR fused to tdTomato (ChR-tdT) was expressed more efficiently than ChrimsonR. Light at 600 nm activated RGCs transfected with AAV2.7m8 ChR-tdT, from an irradiance of 1015 photons.cm-2.s-1. Vector doses of 5 × 1010 and 5 × 1011 vg/eye transfected up to 7000 RGCs/mm2 in the perifovea, with no significant immune reaction. We recorded RGC responses from a stimulus duration of 1 ms upwards. When using the recorded activity to decode stimulus information, we obtained an estimated visual acuity of 20/249, above the level of legal blindness (20/400). These results lay the groundwork for the ongoing clinical trial with the AAV2.7m8 - ChR-tdT vector for vision restoration in patients with retinitis pigmentosa.

Towards optogenetic vision restoration with high resolution.

In many cases of inherited retinal degenerations, ganglion cells are spared despite photoreceptor cell death, making it possible to stimulate them to restore visual function. Several studies have shown that it is possible to express an optogenetic protein in ganglion cells and make them light sensitive, a promising strategy to restore vision. However the spatial resolution of optogenetically-reactivated retinas has rarely been measured, especially in the primate. Since the optogenetic protein is also expressed in axons, it is unclear if these neurons will only be sensitive to the stimulation of a small region covering their somas and dendrites, or if they will also respond to any stimulation overlapping with their axon, dramatically impairing spatial resolution. Here we recorded responses of mouse and macaque retinas to random checkerboard patterns following an in vivo optogenetic therapy. We show that optogenetically activated ganglion cells are each sensitive to a small region of visual space. A simple model based on this small receptive field predicted accurately their responses to complex stimuli. From this model, we simulated how the entire population of light sensitive ganglion cells would respond to letters of different sizes. We then estimated the maximal acuity expected by a patient, assuming it could make an optimal use of the information delivered by this reactivated retina. The obtained acuity is above the limit of legal blindness. Our model also makes interesting predictions on how acuity might vary upon changing the therapeutic strategy, assuming an optimal use of the information present in the retinal activity. Optogenetic therapy could thus potentially lead to high resolution vision, under conditions that our model helps to determinine.

A biophysical model explains the spontaneous bursting behavior in the developing retina.

During early development, waves of activity propagate across the retina and play a key role in the proper wiring of the early visual system. During a particular phase of the retina development (stage II) these waves are triggered by a transient network of neurons, called Starburst Amacrine Cells (SACs), showing a bursting activity which disappears upon further maturation. The underlying mechanisms of the spontaneous bursting and the transient excitability of immature SACs are not completely clear yet. While several models have attempted to reproduce retinal waves, none of them is able to mimic the rhythmic autonomous bursting of individual SACs and reveal how these cells change their intrinsic properties during development. Here, we introduce a mathematical model, grounded on biophysics, which enables us to reproduce the bursting activity of SACs and to propose a plausible, generic and robust, mechanism that generates it. The core parameters controlling repetitive firing are fast depolarizing V-gated calcium channels and hyperpolarizing V-gated potassium channels. The quiescent phase of bursting is controlled by a slow after hyperpolarization (sAHP), mediated by calcium-dependent potassium channels. Based on a bifurcation analysis we show how biophysical parameters, regulating calcium and potassium activity, control the spontaneously occurring fast oscillatory activity followed by long refractory periods in individual SACs. We make a testable experimental prediction on the role of voltage-dependent potassium channels on the excitability properties of SACs and on the evolution of this excitability along development. We also propose an explanation on how SACs can exhibit a large variability in their bursting periods, as observed experimentally within a SACs network as well as across different species, yet based on a simple, unique, mechanism. As we discuss, these observations at the cellular level have a deep impact on the retinal waves description.

Modeling the Correlated Activity of Neural Populations: A Review.

The principles of neural encoding and computations are inherently collective and usually involve large populations of interacting neurons with highly correlated activities. While theories of neural function have long recognized the importance of collective effects in populations of neurons, only in the past two decades has it become possible to record from many cells simultaneously using advanced experimental techniques with single-spike resolution and to relate these correlations to function and behavior. This review focuses on the modeling and inference approaches that have been recently developed to describe the correlated spiking activity of populations of neurons. We cover a variety of models describing correlations between pairs of neurons, as well as between larger groups, synchronous or delayed in time, with or without the explicit influence of the stimulus, and including or not latent variables. We discuss the advantages and drawbacks or each method, as well as the computational challenges related to their application to recordings of ever larger populations.

A Simple Model for Low Variability in Neural Spike Trains.

Neural noise sets a limit to information transmission in sensory systems. In several areas, the spiking response (to a repeated stimulus) has shown a higher degree of regularity than predicted by a Poisson process. However, a simple model to explain this low variability is still lacking. Here we introduce a new model, with a correction to Poisson statistics, that can accurately predict the regularity of neural spike trains in response to a repeated stimulus. The model has only two parameters but can reproduce the observed variability in retinal recordings in various conditions. We show analytically why this approximation can work. In a model of the spike-emitting process where a refractory period is assumed, we derive that our simple correction can well approximate the spike train statistics over a broad range of firing rates. Our model can be easily plugged to stimulus processing models, like a linear-nonlinear model or its generalizations, to replace the Poisson spike train hypothesis that is commonly assumed. It estimates the amount of information transmitted much more accurately than Poisson models in retinal recordings. Thanks to its simplicity, this model has the potential to explain low variability in other areas.