Objects guide human gaze behavior in dynamic real-world scenes.

The complexity of natural scenes makes it challenging to experimentally study the mechanisms behind human gaze behavior when viewing dynamic environments. Historically, eye movements were believed to be driven primarily by space-based attention towards locations with salient features. Increasing evidence suggests, however, that visual attention does not select locations with high saliency but operates on attentional units given by the objects in the scene. We present a new computational framework to investigate the importance of objects for attentional guidance. This framework is designed to simulate realistic scanpaths for dynamic real-world scenes, including saccade timing and smooth pursuit behavior. Individual model components are based on psychophysically uncovered mechanisms of visual attention and saccadic decision-making. All mechanisms are implemented in a modular fashion with a small number of well-interpretable parameters. To systematically analyze the importance of objects in guiding gaze behavior, we implemented five different models within this framework: two purely spatial models, where one is based on low-level saliency and one on high-level saliency, two object-based models, with one incorporating low-level saliency for each object and the other one not using any saliency information, and a mixed model with object-based attention and selection but space-based inhibition of return. We optimized each model's parameters to reproduce the saccade amplitude and fixation duration distributions of human scanpaths using evolutionary algorithms. We compared model performance with respect to spatial and temporal fixation behavior, including the proportion of fixations exploring the background, as well as detecting, inspecting, and returning to objects. A model with object-based attention and inhibition, which uses saliency information to prioritize between objects for saccadic selection, leads to scanpath statistics with the highest similarity to the human data. This demonstrates that scanpath models benefit from object-based attention and selection, suggesting that object-level attentional units play an important role in guiding attentional processing.

Empirically Identifying and Computationally Modeling the Brain-Behavior Relationship for Human Scene Categorization.

Humans effortlessly make quick and accurate perceptual decisions about the nature of their immediate visual environment, such as the category of the scene they face. Previous research has revealed a rich set of cortical representations potentially underlying this feat. However, it remains unknown which of these representations are suitably formatted for decision-making. Here, we approached this question empirically and computationally, using neuroimaging and computational modeling. For the empirical part, we collected EEG data and RTs from human participants during a scene categorization task (natural vs. man-made). We then related EEG data to behavior to behavior using a multivariate extension of signal detection theory. We observed a correlation between neural data and behavior specifically between ∼100 msec and ∼200 msec after stimulus onset, suggesting that the neural scene representations in this time period are suitably formatted for decision-making. For the computational part, we evaluated a recurrent convolutional neural network (RCNN) as a model of brain and behavior. Unifying our previous observations in an image-computable model, the RCNN predicted well the neural representations, the behavioral scene categorization data, as well as the relationship between them. Our results identify and computationally characterize the neural and behavioral correlates of scene categorization in humans.

Interrelations and functional roles of key oscillatory activities during daytime sleep in older adults.

Certain neurophysiological characteristics of sleep, in particular slow oscillations (SOs), sleep spindles, and their temporal coupling, have been well characterised and associated with human memory abilities. Delta waves, which are somewhat higher in frequency and lower in amplitude compared to SOs, and their interaction with spindles have only recently been found to play a critical role in memory processing of rodents, through a competitive interaction between SO-spindle and delta-spindle coupling. However, human studies that comprehensively address delta wave interactions with spindles and SOs, as well as their functional role for memory are still lacking. Electroencephalographic data were acquired across three naps of 33 healthy older human participants (17 female) to investigate delta-spindle coupling and the interplay between delta- and SO-related activity. Additionally, we determined intra-individual stability of coupling measures and their potential link to the ability to form novel memories in a verbal memory task. Our results revealed weaker delta-spindle compared to SO-spindle coupling. Contrary to our initial hypothesis, we found no evidence for an opposing dependency between SO- and delta-related activities during non-rapid eye movement sleep. Moreover, the ratio between SO- and delta-nested spindles rather than SO-spindle and delta-spindle coupling measures by themselves predicted the ability to form novel memories best. In conclusion, our results do not confirm previous findings in rodents on competitive interactions between delta activity and SO-spindle coupling in older adults. However, they support the hypothesis that SO, delta wave, and spindle activity should be jointly considered when aiming to link sleep physiology and memory formation.

Optimal control of a Wilson-Cowan model of neural population dynamics.

Nonlinear dynamical systems describe neural activity at various scales and are frequently used to study brain functions and the impact of external perturbations. Here, we explore methods from optimal control theory (OCT) to study efficient, stimulating "control" signals designed to make the neural activity match desired targets. Efficiency is quantified by a cost functional, which trades control strength against closeness to the target activity. Pontryagin's principle then enables to compute the cost-minimizing control signal. We then apply OCT to a Wilson-Cowan model of coupled excitatory and inhibitory neural populations. The model exhibits an oscillatory regime, low- and high-activity fixed points, and a bistable regime where low- and high-activity states coexist. We compute an optimal control for a state-switching (bistable regime) and a phase-shifting task (oscillatory regime) and allow for a finite transition period before penalizing the deviation from the target state. For the state-switching task, pulses of limited input strength push the activity minimally into the target basin of attraction. Pulse shapes do not change qualitatively when varying the duration of the transition period. For the phase-shifting task, periodic control signals cover the whole transition period. Amplitudes decrease when transition periods are extended, and their shapes are related to the phase sensitivity profile of the model to pulsed perturbations. Penalizing control strength via the integrated 1-norm yields control inputs targeting only one population for both tasks. Whether control inputs drive the excitatory or inhibitory population depends on the state-space location.

Towards Optimization of Oscillatory Stimulation During Sleep.

BACKGROUND: Oscillatory rhythms during sleep, such as slow oscillations (SOs) and spindles and, most importantly, their coupling, are thought to underlie processes of memory consolidation. External slow oscillatory transcranial direct current stimulation (so-tDCS) with a frequency of 0.75 Hz has been shown to improve this coupling and memory consolidation; however, effects varied quite markedly between individuals, studies, and species. In this study, we aimed to determine how precisely the frequency of stimulation must match the naturally occurring SO frequency in individuals to best improve SO-spindle coupling. Moreover, we systematically tested stimulation durations necessary to induce changes. MATERIALS AND METHODS: We addressed these questions by comparing so-tDCS with individualized frequency to standardized frequency of 0.75 Hz in a within-subject design with 28 older participants during napping while stimulation train durations were systematically varied between 30 seconds, 2 minutes, and 5 minutes. RESULTS: Stimulation trains as short as 30 seconds were sufficient to modulate the coupling between SOs and spindle activity. Contrary to our expectations, so-tDCS with standardized frequency indicated stronger aftereffects regarding SO-spindle coupling than individualized frequency. Angle and variance of spindle maxima occurrence during the SO cycle were similarly modulated. CONCLUSIONS: In sum, short stimulation trains were sufficient to induce significant changes in sleep physiology, allowing for more trains of stimulation, which provides methodological advantages and possibly even larger behavioral effects in future studies. Regarding individualized stimulation frequency, further options of optimization need to be investigated, such as closed-loop stimulation, to calibrate stimulation frequency to the SO frequency at the time of stimulation onset. CLINICAL TRIAL REGISTRATION: The Clinicaltrials.gov registration number for the study is NCT04714879.

Correction: Biophysically grounded mean-field models of neural populations under electrical stimulation.

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

Biophysically grounded mean-field models of neural populations under electrical stimulation.

Electrical stimulation of neural systems is a key tool for understanding neural dynamics and ultimately for developing clinical treatments. Many applications of electrical stimulation affect large populations of neurons. However, computational models of large networks of spiking neurons are inherently hard to simulate and analyze. We evaluate a reduced mean-field model of excitatory and inhibitory adaptive exponential integrate-and-fire (AdEx) neurons which can be used to efficiently study the effects of electrical stimulation on large neural populations. The rich dynamical properties of this basic cortical model are described in detail and validated using large network simulations. Bifurcation diagrams reflecting the network's state reveal asynchronous up- and down-states, bistable regimes, and oscillatory regions corresponding to fast excitation-inhibition and slow excitation-adaptation feedback loops. The biophysical parameters of the AdEx neuron can be coupled to an electric field with realistic field strengths which then can be propagated up to the population description. We show how on the edge of bifurcation, direct electrical inputs cause network state transitions, such as turning on and off oscillations of the population rate. Oscillatory input can frequency-entrain and phase-lock endogenous oscillations. Relatively weak electric field strengths on the order of 1 V/m are able to produce these effects, indicating that field effects are strongly amplified in the network. The effects of time-varying external stimulation are well-predicted by the mean-field model, further underpinning the utility of low-dimensional neural mass models.

Weak electric fields promote resonance in neuronal spiking activity: Analytical results from two-compartment cell and network models.

Transcranial brain stimulation and evidence of ephaptic coupling have sparked strong interests in understanding the effects of weak electric fields on the dynamics of neuronal populations. While their influence on the subthreshold membrane voltage can be biophysically well explained using spatially extended neuron models, mechanistic analyses of neuronal spiking and network activity have remained a methodological challenge. More generally, this challenge applies to phenomena for which single-compartment (point) neuron models are oversimplified. Here we employ a pyramidal neuron model that comprises two compartments, allowing to distinguish basal-somatic from apical dendritic inputs and accounting for an extracellular field in a biophysically minimalistic way. Using an analytical approach we fit its parameters to reproduce the response properties of a canonical, spatial model neuron and dissect the stochastic spiking dynamics of single cells and large networks. We show that oscillatory weak fields effectively mimic anti-correlated inputs at the soma and dendrite and strongly modulate neuronal spiking activity in a rather narrow frequency band. This effect carries over to coupled populations of pyramidal cells and inhibitory interneurons, boosting network-induced resonance in the beta and gamma frequency bands. Our work contributes a useful theoretical framework for mechanistic analyses of population dynamics going beyond point neuron models, and provides insights on modulation effects of extracellular fields due to the morphology of pyramidal cells.

Differential polarization of cortical pyramidal neuron dendrites through weak extracellular fields.

The rise of transcranial current stimulation (tCS) techniques have sparked an increasing interest in the effects of weak extracellular electric fields on neural activity. These fields modulate ongoing neural activity through polarization of the neuronal membrane. While the somatic polarization has been investigated experimentally, the frequency-dependent polarization of the dendritic trees in the presence of alternating (AC) fields has received little attention yet. Using a biophysically detailed model with experimentally constrained active conductances, we analyze the subthreshold response of cortical pyramidal cells to weak AC fields, as induced during tCS. We observe a strong frequency resonance around 10-20 Hz in the apical dendrites sensitivity to polarize in response to electric fields but not in the basal dendrites nor the soma. To disentangle the relative roles of the cell morphology and active and passive membrane properties in this resonance, we perform a thorough analysis using simplified models, e.g. a passive pyramidal neuron model, simple passive cables and reconstructed cell model with simplified ion channels. We attribute the origin of the resonance in the apical dendrites to (i) a locally increased sensitivity due to the morphology and to (ii) the high density of h-type channels. Our systematic study provides an improved understanding of the subthreshold response of cortical cells to weak electric fields and, importantly, allows for an improved design of tCS stimuli.

Risk profiles for heavy drinking in adolescence: differential effects of gender.

Abnormalities across different domains of neuropsychological functioning may constitute a risk factor for heavy drinking during adolescence and for developing alcohol use disorders later in life. However, the exact nature of such multi-domain risk profiles is unclear, and it is further unclear whether these risk profiles differ between genders. We combined longitudinal and cross-sectional analyses on the large IMAGEN sample (N ≈ 1000) to predict heavy drinking at age 19 from gray matter volume as well as from psychosocial data at age 14 and 19-for males and females separately. Heavy drinking was associated with reduced gray matter volume in 19-year-olds' bilateral ACC, MPFC, thalamus, middle, medial and superior OFC as well as left amygdala and anterior insula and right inferior OFC. Notably, this lower gray matter volume associated with heavy drinking was stronger in females than in males. In both genders, we observed that impulsivity and facets of novelty seeking at the age of 14 and 19, as well as hopelessness at the age of 14, are risk factors for heavy drinking at the age of 19. Stressful life events with internal (but not external) locus of control were associated with heavy drinking only at age 19. Personality and stress assessment in adolescents may help to better target counseling and prevention programs. This might reduce heavy drinking in adolescents and hence reduce the risk of early brain atrophy, especially in females. In turn, this could additionally reduce the risk of developing alcohol use disorders later in adulthood.

Approximate Inference for Time-Varying Interactions and Macroscopic Dynamics of Neural Populations.

The models in statistical physics such as an Ising model offer a convenient way to characterize stationary activity of neural populations. Such stationary activity of neurons may be expected for recordings from in vitro slices or anesthetized animals. However, modeling activity of cortical circuitries of awake animals has been more challenging because both spike-rates and interactions can change according to sensory stimulation, behavior, or an internal state of the brain. Previous approaches modeling the dynamics of neural interactions suffer from computational cost; therefore, its application was limited to only a dozen neurons. Here by introducing multiple analytic approximation methods to a state-space model of neural population activity, we make it possible to estimate dynamic pairwise interactions of up to 60 neurons. More specifically, we applied the pseudolikelihood approximation to the state-space model, and combined it with the Bethe or TAP mean-field approximation to make the sequential Bayesian estimation of the model parameters possible. The large-scale analysis allows us to investigate dynamics of macroscopic properties of neural circuitries underlying stimulus processing and behavior. We show that the model accurately estimates dynamics of network properties such as sparseness, entropy, and heat capacity by simulated data, and demonstrate utilities of these measures by analyzing activity of monkey V4 neurons as well as a simulated balanced network of spiking neurons.

Spatial separation of two different pathways accounting for the generation of calcium signals in astrocytes.

Astrocytes integrate and process synaptic information and exhibit calcium (Ca2+) signals in response to incoming information from neighboring synapses. The generation of Ca2+ signals is mostly attributed to Ca2+ release from internal Ca2+ stores evoked by an elevated metabotropic glutamate receptor (mGluR) activity. Different experimental results associated the generation of Ca2+ signals to the activity of the glutamate transporter (GluT). The GluT itself does not influence the intracellular Ca2+ concentration, but it indirectly activates Ca2+ entry over the membrane. A closer look into Ca2+ signaling in different astrocytic compartments revealed a spatial separation of those two pathways. Ca2+ signals in the soma are mainly generated by Ca2+ release from internal Ca2+ stores (mGluR-dependent pathway). In astrocytic compartments close to the synapse most Ca2+ signals are evoked by Ca2+ entry over the plasma membrane (GluT-dependent pathway). This assumption is supported by the finding, that the volume ratio between the internal Ca2+ store and the intracellular space decreases from the soma towards the synapse. We extended a model for mGluR-dependent Ca2+ signals in astrocytes with the GluT-dependent pathway. Additionally, we included the volume ratio between the internal Ca2+ store and the intracellular compartment into the model in order to analyze Ca2+ signals either in the soma or close to the synapse. Our model results confirm the spatial separation of the mGluR- and GluT-dependent pathways along the astrocytic process. The model allows to study the binary Ca2+ response during a block of either of both pathways. Moreover, the model contributes to a better understanding of the impact of channel densities on the interaction of both pathways and on the Ca2+ signal.

Low-dimensional spike rate models derived from networks of adaptive integrate-and-fire neurons: Comparison and implementation.

The spiking activity of single neurons can be well described by a nonlinear integrate-and-fire model that includes somatic adaptation. When exposed to fluctuating inputs sparsely coupled populations of these model neurons exhibit stochastic collective dynamics that can be effectively characterized using the Fokker-Planck equation. This approach, however, leads to a model with an infinite-dimensional state space and non-standard boundary conditions. Here we derive from that description four simple models for the spike rate dynamics in terms of low-dimensional ordinary differential equations using two different reduction techniques: one uses the spectral decomposition of the Fokker-Planck operator, the other is based on a cascade of two linear filters and a nonlinearity, which are determined from the Fokker-Planck equation and semi-analytically approximated. We evaluate the reduced models for a wide range of biologically plausible input statistics and find that both approximation approaches lead to spike rate models that accurately reproduce the spiking behavior of the underlying adaptive integrate-and-fire population. Particularly the cascade-based models are overall most accurate and robust, especially in the sensitive region of rapidly changing input. For the mean-driven regime, when input fluctuations are not too strong and fast, however, the best performing model is based on the spectral decomposition. The low-dimensional models also well reproduce stable oscillatory spike rate dynamics that are generated either by recurrent synaptic excitation and neuronal adaptation or through delayed inhibitory synaptic feedback. The computational demands of the reduced models are very low but the implementation complexity differs between the different model variants. Therefore we have made available implementations that allow to numerically integrate the low-dimensional spike rate models as well as the Fokker-Planck partial differential equation in efficient ways for arbitrary model parametrizations as open source software. The derived spike rate descriptions retain a direct link to the properties of single neurons, allow for convenient mathematical analyses of network states, and are well suited for application in neural mass/mean-field based brain network models.

Interaction of Instrumental and Goal-Directed Learning Modulates Prediction Error Representations in the Ventral Striatum.

Goal-directed and instrumental learning are both important controllers of human behavior. Learning about which stimulus event occurs in the environment and the reward associated with them allows humans to seek out the most valuable stimulus and move through the environment in a goal-directed manner. Stimulus-response associations are characteristic of instrumental learning, whereas response-outcome associations are the hallmark of goal-directed learning. Here we provide behavioral, computational, and neuroimaging results from a novel task in which stimulus-response and response-outcome associations are learned simultaneously but dominate behavior at different stages of the experiment. We found that prediction error representations in the ventral striatum depend on which type of learning dominates. Furthermore, the amygdala tracks the time-dependent weighting of stimulus-response versus response-outcome learning. Our findings suggest that the goal-directed and instrumental controllers dynamically engage the ventral striatum in representing prediction errors whenever one of them is dominating choice behavior. SIGNIFICANCE STATEMENT: Converging evidence in human neuroimaging studies has shown that the reward prediction errors are correlated with activity in the ventral striatum. Our results demonstrate that this region is simultaneously correlated with a stimulus prediction error. Furthermore, the learning system that is currently dominating behavioral choice dynamically engages the ventral striatum for computing its prediction errors. This demonstrates that the prediction error representations are highly dynamic and influenced by various experimental context. This finding points to a general role of the ventral striatum in detecting expectancy violations and encoding error signals regardless of the specific nature of the reinforcer itself.

Extending Integrate-and-Fire Model Neurons to Account for the Effects of Weak Electric Fields and Input Filtering Mediated by the Dendrite.

Transcranial brain stimulation and evidence of ephaptic coupling have recently sparked strong interests in understanding the effects of weak electric fields on the dynamics of brain networks and of coupled populations of neurons. The collective dynamics of large neuronal populations can be efficiently studied using single-compartment (point) model neurons of the integrate-and-fire (IF) type as their elements. These models, however, lack the dendritic morphology required to biophysically describe the effect of an extracellular electric field on the neuronal membrane voltage. Here, we extend the IF point neuron models to accurately reflect morphology dependent electric field effects extracted from a canonical spatial "ball-and-stick" (BS) neuron model. Even in the absence of an extracellular field, neuronal morphology by itself strongly affects the cellular response properties. We, therefore, derive additional components for leaky and nonlinear IF neuron models to reproduce the subthreshold voltage and spiking dynamics of the BS model exposed to both fluctuating somatic and dendritic inputs and an extracellular electric field. We show that an oscillatory electric field causes spike rate resonance, or equivalently, pronounced spike to field coherence. Its resonance frequency depends on the location of the synaptic background inputs. For somatic inputs the resonance appears in the beta and gamma frequency range, whereas for distal dendritic inputs it is shifted to even higher frequencies. Irrespective of an external electric field, the presence of a dendritic cable attenuates the subthreshold response at the soma to slowly-varying somatic inputs while implementing a low-pass filter for distal dendritic inputs. Our point neuron model extension is straightforward to implement and is computationally much more efficient compared to the original BS model. It is well suited for studying the dynamics of large populations of neurons with heterogeneous dendritic morphology with (and without) the influence of weak external electric fields.

Altered behavioral and neural responsiveness to counterfactual gains in the elderly.

Counterfactual information processing refers to the consideration of events that did not occur in comparison to those actually experienced, in order to determine optimal actions, and can be formulated as computational learning signals, referred to as fictive prediction errors. Decision making and the neural circuitry for counterfactual processing are altered in healthy elderly adults. This experiment investigated age differences in neural systems for decision making with knowledge of counterfactual outcomes. Two groups of healthy adult participants, young (N = 30; ages 19-30 years) and elderly (N = 19; ages 65-80 years), were scanned with fMRI during 240 trials of a strategic sequential investment task in which a particular strategy of differentially weighting counterfactual gains and losses during valuation is associated with more optimal performance. Elderly participants earned significantly less than young adults, differently weighted counterfactual consequences and exploited task knowledge, and exhibited altered activity in a fronto-striatal circuit while making choices, compared to young adults. The degree to which task knowledge was exploited was positively correlated with modulation of neural activity by expected value in the vmPFC for young adults, but not in the elderly. These findings demonstrate that elderly participants' poor task performance may be related to different counterfactual processing.

Integrator or Coincidence Detector: A Novel Measure Based on the Discrete Reverse Correlation to Determine a Neuron's Operational Mode.

In this letter, we propose a definition of the operational mode of a neuron, that is, whether a neuron integrates over its input or detects coincidences. We complete the range of possible operational modes by a new mode we call gap detection, which means that a neuron responds to gaps in its stimulus. We propose a measure consisting of two scalar values, both ranging from -1 to +1: the neural drive, which indicates whether its stimulus excites the neuron, serves as background noise, or inhibits it; the neural mode, which indicates whether the neuron's response is the result of integration over its input, of coincidence detection, or of gap detection; with all three modes possible for all neural drive values. This is a pure spike-based measure and can be applied to measure the influence of either all or subset of a neuron's stimulus. We derive the measure by decomposing the reverse correlation, test it in several artificial and biological settings, and compare it to other measures, finding little or no correlation between them. We relate the results of the measure to neural parameters and investigate the effect of time delay during spike generation. Our results suggest that a neuron can use several different modes simultaneously on different subsets of its stimulus to enable it to respond to its stimulus in a complex manner.

A hippocampal signature of perceptual learning in object recognition.

Perceptual learning is the improvement in perceptual performance through training or exposure. Here, we used fMRI before and after extensive behavioral training to investigate the effects of perceptual learning on the recognition of objects under challenging viewing conditions. Objects belonged either to trained or untrained categories. Trained categories were further subdivided into trained and untrained exemplars and were coupled with high or low monetary rewards during training. After a 3-day training, object recognition was markedly improved. Although there was a considerable transfer of learning to untrained exemplars within categories, an enhancing effect of reward reinforcement was specific to trained exemplars. fMRI showed that hippocampus responses to both trained and untrained exemplars of trained categories were enhanced by perceptual learning and correlated with the effect of reward reinforcement. Our results suggest a key role of hippocampus in object recognition after perceptual learning.

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.

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.