Bayesian hypothesis testing and experimental design for two-photon imaging data.

Variability, stochastic or otherwise, is a central feature of neural activity. Yet the means by which estimates of variation and uncertainty are derived from noisy observations of neural activity is often heuristic, with more weight given to numerical convenience than statistical rigour. For two-photon imaging data, composed of fundamentally probabilistic streams of photon detections, the problem is particularly acute. Here, we present a statistical pipeline for the inference and analysis of neural activity using Gaussian Process regression, applied to two-photon recordings of light-driven activity in ex vivo mouse retina. We demonstrate the flexibility and extensibility of these models, considering cases with non-stationary statistics, driven by complex parametric stimuli, in signal discrimination, hierarchical clustering and other inference tasks. Sparse approximation methods allow these models to be fitted rapidly, permitting them to actively guide the design of light stimulation in the midst of ongoing two-photon experiments.

Correction: Bayesian hypothesis testing and experimental design for two-photon imaging data.

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

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.

Connectomics of synaptic microcircuits: lessons from the outer retina.

Photoreceptors form a sophisticated synaptic complex with bipolar and horizontal cells, transmitting the signals generated by the phototransduction cascade to downstream retinal circuitry. The cone photoreceptor synapse shows several characteristic anatomical connectivity motifs that shape signal transfer: typically, ON-cone bipolar cells receive photoreceptor input through invaginating synapses; OFF-cone bipolar cells form basal synapses with photoreceptors. Both ON- and OFF-cone bipolar cells are believed to sample from all cone photoreceptors within their dendritic span. Electron microscopy and immunolabelling studies have established the robustness of these motifs, but have been limited by trade-offs in sample size and spatial resolution, respectively, constraining precise quantitative investigation to a few individual cells. 3D-serial electron microscopy overcomes these limitations and has permitted complete sets of neurons to be reconstructed over a comparatively large section of retinal tissue. Although the published mouse dataset lacks labels for synaptic structures, the characteristic anatomical motifs at the photoreceptor synapse can be exploited to identify putative synaptic contacts, which has enabled the development of a quantitative description of outer retinal connectivity. This revealed unexpected exceptions to classical motifs, including substantial interaction between rod and cone pathways at the photoreceptor synapse, sparse photoreceptor sampling and atypical contacts. Here, we summarize what was learned from this study in a more general context: we consider both the implications and limitations of the study and identify promising avenues for future research.

Redefining the role of Ca2+-permeable channels in photoreceptor degeneration using diltiazem.

Hereditary degeneration of photoreceptors has been linked to over-activation of Ca2+-permeable channels, excessive Ca2+-influx, and downstream activation of Ca2+-dependent calpain-type proteases. Unfortunately, after more than 20 years of pertinent research, unequivocal evidence proving significant and reproducible photoreceptor protection with Ca2+-channel blockers is still lacking. Here, we show that both D- and L-cis enantiomers of the anti-hypertensive drug diltiazem were very effective at blocking photoreceptor Ca2+-influx, most probably by blocking the pore of Ca2+-permeable channels. Yet, unexpectedly, this block neither reduced the activity of calpain-type proteases, nor did it result in photoreceptor protection. Remarkably, application of the L-cis enantiomer of diltiazem even led to a strong increase in photoreceptor cell death. These findings shed doubt on the previously proposed links between Ca2+ and retinal degeneration and are highly relevant for future therapy development as they may serve to refocus research efforts towards alternative, Ca2+-independent degenerative mechanisms.

Systematic spatiotemporal mapping reveals divergent cell death pathways in three mouse models of hereditary retinal degeneration.

Calcium (Ca2+ ) dysregulation has been linked to neuronal cell death, including in hereditary retinal degeneration. Ca2+ dysregulation is thought to cause rod and cone photoreceptor cell death. Spatial and temporal heterogeneities in retinal disease models have hampered validation of this hypothesis. We examined the role of Ca2+ in photoreceptor degeneration, assessing the activation pattern of Ca2+ -dependent calpain proteases, generating spatiotemporal maps of the entire retina in the cpfl1 mouse model for primary cone degeneration, and in the rd1 and rd10 models for primary rod degeneration. We used Gaussian process models to distinguish the temporal sequences of degenerative molecular processes from other variability sources.In the rd1 and rd10 models, spatiotemporal pattern of increased calpain activity matched the progression of primary rod degeneration. High calpain activity coincided with activation of the calpain-2 isoform but not with calpain-1, suggesting differential roles for both calpain isoforms. Primary rod loss was linked to upregulation of apoptosis-inducing factor, although only a minute fraction of cells showed activity of the apoptotic marker caspase-3. After primary rod degeneration concluded, caspase-3 activation appeared in cones, suggesting apoptosis as the dominant mechanism for secondary cone loss. Gaussian process models highlighted calpain activity as a key event during primary rod photoreceptor cell death. Our data suggest a causal link between Ca2+ dysregulation and primary, nonapoptotic degeneration of photoreceptors and a role for apoptosis in secondary degeneration of cones, highlighting the importance of the spatial and temporal location of key molecular events, which may guide the evaluation of new therapies.

The temporal structure of the inner retina at a single glance.

The retina decomposes visual stimuli into parallel channels that encode different features of the visual environment. Central to this computation is the synaptic processing in a dense layer of neuropil, the so-called inner plexiform layer (IPL). Here, different types of bipolar cells stratifying at distinct depths relay the excitatory feedforward drive from photoreceptors to amacrine and ganglion cells. Current experimental techniques for studying processing in the IPL do not allow imaging the entire IPL simultaneously in the intact tissue. Here, we extend a two-photon microscope with an electrically tunable lens allowing us to obtain optical vertical slices of the IPL, which provide a complete picture of the response diversity of bipolar cells at a "single glance". The nature of these axial recordings additionally allowed us to isolate and investigate batch effects, i.e. inter-experimental variations resulting in systematic differences in response speed. As a proof of principle, we developed a simple model that disentangles biological from experimental causes of variability and allowed us to recover the characteristic gradient of response speeds across the IPL with higher precision than before. Our new framework will make it possible to study the computations performed in the central synaptic layer of the retina more efficiently.

Local Signals in Mouse Horizontal Cell Dendrites.

The mouse retina contains a single type of horizontal cell, a GABAergic interneuron that samples from all cone photoreceptors within reach and modulates their glutamatergic output via parallel feedback mechanisms. Because horizontal cells form an electrically coupled network, they have been implicated in global signal processing, such as large-scale contrast enhancement. Recently, it has been proposed that horizontal cells can also act locally at the level of individual cone photoreceptors. To test this possibility physiologically, we used two-photon microscopy to record light stimulus-evoked Ca2+ signals in cone axon terminals and horizontal cell dendrites as well as glutamate release in the outer plexiform layer. By selectively stimulating the two mouse cone opsins with green and UV light, we assessed whether signals from individual cones remain isolated within horizontal cell dendritic tips or whether they spread across the dendritic arbor. Consistent with the mouse's opsin expression gradient, we found that the Ca2+ signals recorded from dendrites of dorsal horizontal cells were dominated by M-opsin and those of ventral horizontal cells by S-opsin activation. The signals measured in neighboring horizontal cell dendritic tips varied markedly in their chromatic preference, arguing against global processing. Rather, our experimental data and results from biophysically realistic modeling support the idea that horizontal cells can process cone input locally, extending the classical view of horizontal cell function. Pharmacologically removing horizontal cells from the circuitry reduced the sensitivity of the cone signal to low frequencies, suggesting that local horizontal cell feedback shapes the temporal properties of cone output.