Bursting in cerebellar stellate cells induced by pharmacological agents: Non-sequential spike adding

ElecFeX is a user-friendly toolbox for efficient feature extraction from single-cell electrophysiological recordings

Characterizing neurons by their electrophysiological phenotypes is essential for understanding the neural basis of behavioral and cognitive functions. Technological developments have enabled the collection of hundreds of neural recordings; this calls for new tools capable of performing feature extraction efficiently. To address the urgent need for a powerful and accessible tool, we developed ElecFeX, an open-source MATLAB-based toolbox that (1) has an intuitive graphical user interface, (2) provides customizable measurements for a wide range of electrophysiological features, (3) processes large-size datasets effortlessly via batch analysis, and (4) yields formatted output for further analysis. We implemented ElecFeX on a diverse set of neural recordings; demonstrated its functionality, versatility, and efficiency in capturing electrical features; and established its significance in distinguishing neuronal subgroups across brain regions and species. ElecFeX is thus presented as a user-friendly toolbox to benefit the neuroscience community by minimizing the time required for extracting features from their electrophysiological datasets.

Classification of bursting patterns: A tale of two ducks

Bursting is one of the fundamental rhythms that excitable cells can generate either in response to incoming stimuli or intrinsically. It has been a topic of intense research in computational biology for several decades. The classification of bursting oscillations in excitable systems has been the subject of active research since the early 1980s and is still ongoing. As a by-product, it establishes analytical and numerical foundations for studying complex temporal behaviors in multiple timescale models of cellular activity. In this review, we first present the seminal works of Rinzel and Izhikevich in classifying bursting patterns of excitable systems. We recall a complementary mathematical classification approach by Bertram and colleagues, and then by Golubitsky and colleagues, which, together with the Rinzel-Izhikevich proposals, provide the state-of-the-art foundations to these classifications. Beyond classical approaches, we review a recent bursting example that falls outside the previous classification systems. Generalizing this example leads us to propose an extended classification, which requires the analysis of both fast and slow subsystems of an underlying slow-fast model and allows the dissection of a larger class of bursters. Namely, we provide a general framework for bursting systems with both subthreshold and superthreshold oscillations. A new class of bursters with at least 2 slow variables is then added, which we denote folded-node bursters, to convey the idea that the bursts are initiated or annihilated via a folded-node singularity. Key to this mechanism are so-called canard or duck orbits, organizing the underpinning excitability structure. We describe the 2 main families of folded-node bursters, depending upon the phase (active/spiking or silent/nonspiking) of the bursting cycle during which folded-node dynamics occurs. We classify both families and give examples of minimal systems displaying these novel bursting patterns. Finally, we provide a biophysical example by reinterpreting a generic conductance-based episodic burster as a folded-node burster, showing that the associated framework can explain its subthreshold oscillations over a larger parameter region than the fast subsystem approach.

Bursting is ubiquitous in cellular excitable rhythms and comes in a plethora of patterns, both experimentally recorded and reproduced through models. As these different patterns may reflect different coding or information properties, it is therefore crucial to develop modeling frameworks that can both capture them and understand their characteristics. In this review, we propose a comprehensive account of the main bursting classification systems that have been developed over the past 40 years, together with recent developments allowing us to extend these classifications. Based upon bifurcation theory and heavily reliant on timescale separation, these schemes take full advantage of the fast subsystem analysis, obtained when slow variables are frozen and considered as bifurcation parameters. We complement this classical view by showing that nontrivial slow subsystem may also encode key informations important to classify bursting rhythms, due to the presence of so-called folded-node singularities. We provide minimal idealized models as well as one generic conductance-based example displaying bursting oscillations that require our extended classification in order to be fully characterized. We also highlight examples of biological data that could be suitably revisited with the lenses of this extended classifications and could lead to new models of complex cellular activity.

Bursting in cerebellar stellate cells induced by pharmacological agents: Non-sequential spike adding

Cerebellar stellate cells (CSCs) are spontaneously active, tonically firing (5-30 Hz), inhibitory interneurons that synapse onto Purkinje cells. We previously analyzed the excitability properties of CSCs, focusing on four key features: type I excitability, non-monotonic first-spike latency, switching in responsiveness and runup (i.e., temporal increase in excitability during whole-cell configuration). In this study, we extend this analysis by using whole-cell configuration to show that these neurons can also burst when treated with certain pharmacological agents separately or jointly. Indeed, treatment with 4-Aminopyridine (4-AP), a partial blocker of delayed rectifier and A-type K⁺ channels, at low doses induces a bursting profile in CSCs significantly different than that produced at high doses or when it is applied at low doses but with cadmium (Cd²⁺), a blocker of high voltage-activated (HVA) Ca²⁺ channels. By expanding a previously revised Hodgkin–Huxley type model, through the inclusion of Ca²⁺-activated K⁺ (K(Ca)) and HVA currents, we explain how these bursts are generated and what their underlying dynamics are. Specifically, we demonstrate that the expanded model preserves the four excitability features of CSCs, as well as captures their bursting patterns induced by 4-AP and Cd²⁺. Model investigation reveals that 4-AP is potentiating HVA, inducing square-wave bursting at low doses and pseudo-plateau bursting at high doses, whereas Cd²⁺ is potentiating K(Ca), inducing pseudo-plateau bursting when applied in combination with low doses of 4-AP. Using bifurcation analysis, we show that spike adding in square-wave bursts is non-sequential when gradually changing HVA and K(Ca) maximum conductances, delayed Hopf is responsible for generating the plateau segment within the active phase of pseudo-plateau bursts, and bursting can become “chaotic” when HVA and K(Ca) maximum conductances are made low and high, respectively. These results highlight the secondary effects of the drugs applied and suggest that CSCs have all the ingredients needed for bursting.

Excitable cells, including neurons, fire action potentials (APs) in their membrane voltage that allow them to communicate with each other and to serve certain physiological purposes. They do so either tonically by firing APs periodically, or episodically by repeatedly firing clusters of APs (called bursts) separated by quiescent periods. Each one of those firing patterns can be neuron-specific and dependent on synaptic inputs and/or their physiological environment. Cerebellar stellate cells (CSCs) that synapse onto Purkinje cells, the sole output of the cerebellum responsible for motor control, are spontaneously active inhibitory interneurons that fire APs tonically. We previously studied the excitability properties of these neurons and showed that they possess several important key features, including type I excitability, runup, non-monotonic first spike latency and switching in responsiveness. In this study, we show that CSCs can also exhibit two modes of burst firing, called square-wave and pseudo-plateau, when treated with certain pharmacological agents. Using bifurcation theory, we demonstrate that spike adding in the square-wave burst is non-sequential, changing by several spikes when certain conductances are altered gradually. This study thus sheds lights onto the overall effects of the pharmacological agents and highlights the ability of CSCs to burst in certain biological conditions.