Directional data is ubiquitious in science. Due to its circular nature such data cannot be analyzed with commonly used statistical techniques. Despite the rapid development of specialized methods for ...directional statistics over the last fifty years, there is only little software available that makes such methods easy to use for practioners. Most importantly, one of the most commonly used programming languages in biosciences, MATLAB, is currently not supporting directional statistics. To remedy this situation, we have implemented the CircStat toolbox for MATLAB which provides methods for the descriptive and inferential statistical analysis of directional data. We cover the statistical background of the available methods and describe how to apply them to data. Finally, we analyze a dataset from neurophysiology to demonstrate the capabilities of the CircStat toolbox.
Single-cell transcriptomics yields ever growing data sets containing RNA expression levels for thousands of genes from up to millions of cells. Common data analysis pipelines include a dimensionality ...reduction step for visualising the data in two dimensions, most frequently performed using t-distributed stochastic neighbour embedding (t-SNE). It excels at revealing local structure in high-dimensional data, but naive applications often suffer from severe shortcomings, e.g. the global structure of the data is not represented accurately. Here we describe how to circumvent such pitfalls, and develop a protocol for creating more faithful t-SNE visualisations. It includes PCA initialisation, a high learning rate, and multi-scale similarity kernels; for very large data sets, we additionally use exaggeration and downsampling-based initialisation. We use published single-cell RNA-seq data sets to demonstrate that this protocol yields superior results compared to the naive application of t-SNE.
Deep learning (DL) has revolutionized the field of computer vision and image processing. In medical imaging, algorithmic solutions based on DL have been shown to achieve high performance on tasks ...that previously required medical experts. However, DL-based solutions for disease detection have been proposed without methods to quantify and control their uncertainty in a decision. In contrast, a physician knows whether she is uncertain about a case and will consult more experienced colleagues if needed. Here we evaluate drop-out based Bayesian uncertainty measures for DL in diagnosing diabetic retinopathy (DR) from fundus images and show that it captures uncertainty better than straightforward alternatives. Furthermore, we show that uncertainty informed decision referral can improve diagnostic performance. Experiments across different networks, tasks and datasets show robust generalization. Depending on network capacity and task/dataset difficulty, we surpass 85% sensitivity and 80% specificity as recommended by the NHS when referring 0-20% of the most uncertain decisions for further inspection. We analyse causes of uncertainty by relating intuitions from 2D visualizations to the high-dimensional image space. While uncertainty is sensitive to clinically relevant cases, sensitivity to unfamiliar data samples is task dependent, but can be rendered more robust.
Standard preprocessing of single-cell RNA-seq UMI data includes normalization by sequencing depth to remove this technical variability, and nonlinear transformation to stabilize the variance across ...genes with different expression levels. Instead, two recent papers propose to use statistical count models for these tasks: Hafemeister and Satija (Genome Biol 20:296, 2019) recommend using Pearson residuals from negative binomial regression, while Townes et al. (Genome Biol 20:295, 2019) recommend fitting a generalized PCA model. Here, we investigate the connection between these approaches theoretically and empirically, and compare their effects on downstream processing.
We show that the model of Hafemeister and Satija produces noisy parameter estimates because it is overspecified, which is why the original paper employs post hoc smoothing. When specified more parsimoniously, it has a simple analytic solution equivalent to the rank-one Poisson GLM-PCA of Townes et al. Further, our analysis indicates that per-gene overdispersion estimates in Hafemeister and Satija are biased, and that the data are in fact consistent with the overdispersion parameter being independent of gene expression. We then use negative control data without biological variability to estimate the technical overdispersion of UMI counts, and find that across several different experimental protocols, the data are close to Poisson and suggest very moderate overdispersion. Finally, we perform a benchmark to compare the performance of Pearson residuals, variance-stabilizing transformations, and GLM-PCA on scRNA-seq datasets with known ground truth.
We demonstrate that analytic Pearson residuals strongly outperform other methods for identifying biologically variable genes, and capture more of the biologically meaningful variation when used for dimensionality reduction.
In the vertebrate visual system, all output of the retina is carried by retinal ganglion cells. Each type encodes distinct visual features in parallel for transmission to the brain. How many such ...'output channels' exist and what each encodes are areas of intense debate. In the mouse, anatomical estimates range from 15 to 20 channels, and only a handful are functionally understood. By combining two-photon calcium imaging to obtain dense retinal recordings and unsupervised clustering of the resulting sample of more than 11,000 cells, here we show that the mouse retina harbours substantially more than 30 functional output channels. These include all known and several new ganglion cell types, as verified by genetic and anatomical criteria. Therefore, information channels from the mouse eye to the mouse brain are considerably more diverse than shown thus far by anatomical studies, suggesting an encoding strategy resembling that used in state-of-the-art artificial vision systems.
In the mouse retina, three different types of photoreceptors provide input to 14 bipolar cell (BC) types. Classically, most BC types are thought to contact all cones within their dendritic field; ...ON-BCs would contact cones exclusively via so-called invaginating synapses, while OFF-BCs would form basal synapses. By mining publically available electron microscopy data, we discovered interesting violations of these rules of outer retinal connectivity: ON-BC type X contacted only ~20% of the cones in its dendritic field and made mostly atypical non-invaginating contacts. Types 5T, 5O and 8 also contacted fewer cones than expected. In addition, we found that rod BCs received input from cones, providing anatomical evidence that rod and cone pathways are interconnected in both directions. This suggests that the organization of the outer plexiform layer is more complex than classically thought.
In the dentate gyrus - a key component of spatial memory circuits - granule cells (GCs) are known to be morphologically diverse and to display heterogeneous activity profiles during behavior. To ...resolve structure-function relationships, we juxtacellularly recorded and labeled single GCs in freely moving rats. We found that the vast majority of neurons were silent during exploration. Most active GCs displayed a characteristic spike waveform, fired at low rates and showed spatial activity. Primary dendritic parameters were sufficient for classifying neurons as active or silent with high accuracy. Our data thus support a sparse coding scheme in the dentate gyrus and provide a possible link between structural and functional heterogeneity among the GC population.
In the eye, the function of same-type photoreceptors must be regionally adjusted to process a highly asymmetrical natural visual world. Here, we show that UV cones in the larval zebrafish area ...temporalis are specifically tuned for UV-bright prey capture in their upper frontal visual field, which may use the signal from a single cone at a time. For this, UV-photon detection probability is regionally boosted more than 10-fold. Next, in vivo two-photon imaging, transcriptomics, and computational modeling reveal that these cones use an elevated baseline of synaptic calcium to facilitate the encoding of bright objects, which in turn results from expressional tuning of phototransduction genes. Moreover, the light-driven synaptic calcium signal is regionally slowed by interactions with horizontal cells and later accentuated at the level of glutamate release driving retinal networks. These regional differences tally with variations between peripheral and foveal cones in primates and hint at a common mechanistic origin.
Display omitted
•Larval zebrafish prey detection uses specialized single UV cones in the acute zone•High-gain UV cones in this region boost detection of UV-bright contrasts•Cellular and molecular mechanisms of this tuning tally with those of the primate fovea•Further optimization occurs at the level of UV cone glutamate release
Yoshimatsu et al. show that larval zebrafish rely on single UV cones at a time to support visual prey capture. For this, zebrafish combine molecular, cellular, and circuit tuning to regionally boost detectability of prey in their acute zone. The mechanisms of this specialization tally with those of the primate fovea.
Abstract
Color vision is essential for an animal’s survival. It starts in the retina, where signals from different photoreceptor types are locally compared by neural circuits. Mice, like most ...mammals, are dichromatic with two cone types. They can discriminate colors only in their upper visual field. In the corresponding ventral retina, however, most cones display the same spectral preference, thereby presumably impairing spectral comparisons. In this study, we systematically investigated the retinal circuits underlying mouse color vision by recording light responses from cones, bipolar and ganglion cells. Surprisingly, most color-opponent cells are located in the ventral retina, with rod photoreceptors likely being involved. Here, the complexity of chromatic processing increases from cones towards the retinal output, where non-linear center-surround interactions create specific color-opponent output channels to the brain. This suggests that neural circuits in the mouse retina are tuned to extract color from the upper visual field, aiding robust detection of predators and ensuring the animal’s survival.
Neural computation relies on the integration of synaptic inputs across a neuron's dendritic arbour. However, it is far from understood how different cell types tune this process to establish ...cell-type specific computations. Here, using two-photon imaging of dendritic Ca
signals, electrical recordings of somatic voltage and biophysical modelling, we demonstrate that four morphologically distinct types of mouse retinal ganglion cells with overlapping excitatory synaptic input (transient Off alpha, transient Off mini, sustained Off, and F-mini Off) exhibit type-specific dendritic integration profiles: in contrast to the other types, dendrites of transient Off alpha cells were spatially independent, with little receptive field overlap. The temporal correlation of dendritic signals varied also extensively, with the highest and lowest correlation in transient Off mini and transient Off alpha cells, respectively. We show that differences between cell types can likely be explained by differences in backpropagation efficiency, arising from the specific combinations of dendritic morphology and ion channel densities.