Action Potential Spike Research Articles

Characterization of Early Cortical Neural Network Development in Multiwell Microelectrode Array Plates

We examined neural network ontogeny using microelectrode array (MEA) recordings made in multiwell MEA (mwMEA) plates over the first 12 days in vitro (DIV). In primary cortical cultures, action potential spiking activity developed rapidly between DIV 5 and 12. Spiking was sporadic and unorganized at early DIV, and became progressively more organized with time, with bursting parameters, synchrony, and network bursting increasing between DIV 5 and 12. We selected 12 features to describe network activity; principal components analysis using these features demonstrated segregation of data by age at both the well and plate levels. Using random forest classifiers and support vector machines, we demonstrated that four features (coefficient of variation [CV] of within-burst interspike interval, CV of interburst interval, network spike rate, and burst rate) could predict the age of each well recording with >65% accuracy. When restricting the classification to a binary decision, accuracy improved to as high as 95%. Further, we present a novel resampling approach to determine the number of wells needed for comparing different treatments. Overall, these results demonstrate that network development on mwMEA plates is similar to development in single-well MEAs. The increased throughput of mwMEAs will facilitate screening drugs, chemicals, or disease states for effects on neurodevelopment.

Genetic alteration of the metal/redox modulation of Cav3.2 T‐type calcium channel reveals its role in neuronal excitability

In this study, we describe a new knock-in (KI) mouse model that allows the study of the H191-dependent regulation of T-type Cav3.2 channels. Sensitivity to zinc, nickel and ascorbate of native Cav3.2 channels is significantly impeded in the dorsal root ganglion (DRG) neurons of this KI mouse. Importantly, we describe that this H191-dependent regulation has discrete but significant effects on the excitability properties of D-hair (down-hair) cells, a sub-population of DRG neurons in which Cav3.2 currents prominently regulate excitability. Overall, this study reveals that the native H191-dependent regulation of Cav3.2 channels plays a role in the excitability of Cav3.2-expressing neurons. This animal model will be valuable in addressing the potential in vivo roles of the trace metal and redox modulation of Cav3.2 T-type channels in a wide range of physiological and pathological conditions. Cav3.2 channels are T-type voltage-gated calcium channels that play important roles in controlling neuronal excitability, particularly in dorsal root ganglion (DRG) neurons where they are involved in touch and pain signalling. Cav3.2 channels are modulated by low concentrations of metal ions (nickel, zinc) and redox agents, which involves the histidine 191 (H191) in the channel's extracellular IS3-IS4 loop. It is hypothesized that this metal/redox modulation would contribute to the tuning of the excitability properties of DRG neurons. However, the precise role of this H191-dependent modulation of Cav3.2 channel remains unresolved. Towards this goal, we have generated a knock-in (KI) mouse carrying the mutation H191Q in the Cav3.2 protein. Electrophysiological studies were performed on a subpopulation of DRG neurons, the D-hair cells, which express large Cav3.2 currents. We describe an impaired sensitivity to zinc, nickel and ascorbate of the T-type current in D-hair neurons from KI mice. Analysis of the action potential and low-threshold calcium spike (LTCS) properties revealed that, contrary to that observed in WT D-hair neurons, a low concentration of zinc and nickel is unable to modulate (1) the rheobase threshold current, (2) the afterdepolarization amplitude, (3) the threshold potential necessary to trigger an LTCS or (4) the LTCS amplitude in D-hair neurons from KI mice. Together, our data demonstrate that this H191-dependent metal/redox regulation of Cav3.2 channels can tune neuronal excitability. This study validates the use of this Cav3.2-H191Q mouse model for further investigations of the physiological roles thought to rely on this Cav3.2 modulation.

Central Nervous System–Toxic Lidocaine Concentrations Unmask L-Type Ca2+ Current–Mediated Action Potentials in Rat Thalamocortical Neurons: An In Vitro Mechanism of Action Study

High systemic lidocaine concentrations exert well-known toxic effects on the central nervous system (CNS), including seizures, coma, and death. The underlying mechanisms are still largely obscure, and the actions of lidocaine on supraspinal neurons have received comparatively little study. We recently found that lidocaine at clinically neurotoxic concentrations increases excitability mediated by Na-independent, high-threshold (HT) action potential spikes in rat thalamocortical neurons. Our goal in this study was to characterize these spikes and test the hypothesis that they are generated by HT Ca currents, previously implicated in neurotoxicity. We also sought to identify and isolate the specific underlying subtype of Ca current. We investigated the actions of lidocaine in the CNS-toxic concentration range (100 μM-1 mM) on ventrobasal thalamocortical neurons in rat brain slices in vitro, using whole-cell patch-clamp recordings aided by differential interference contrast infrared videomicroscopy. Drugs were bath applied; action potentials were generated using current clamp protocols, and underlying currents were identified and isolated with ion channel blockers and electrolyte substitution. Lidocaine (100 μM-1 mM) abolished Na-dependent tonic firing in all neurons tested (n = 46). However, in 39 of 46 (85%) neurons, lidocaine unmasked evoked HT action potentials with lower amplitudes and rates of de-/repolarization compared with control. These HT action potentials remained during the application of tetrodotoxin (600 nM), were blocked by Cd (50 μM), and disappeared after superfusion with an extracellular solution deprived of Ca. These features implied that the unmasked potentials were generated by high-voltage-activated Ca channels and not by Na channels. Application of the L-type Ca channel blocker, nifedipine (5 μM), completely blocked the HT potentials, whereas the N-type Ca channel blocker, ω-conotoxin GVIA (1 μM), had little effect. At clinically CNS-toxic concentrations, lidocaine unmasked in thalamocortical neurons evoked HT action potentials mediated by the L-type Ca current while substantially suppressing Na-dependent excitability. On the basis of the known role of an increase in intracellular Ca in the pathogenesis of local anesthetic neurotoxicity, this novel action represents a plausible contributing candidate mechanism for lidocaine's CNS toxicity in vivo.

Power-Law Dynamics of Membrane Conductances Increase Spiking Diversity in a Hodgkin-Huxley Model.

We studied the effects of non-Markovian power-law voltage dependent conductances on the generation of action potentials and spiking patterns in a Hodgkin-Huxley model. To implement slow-adapting power-law dynamics of the gating variables of the potassium, n, and sodium, m and h, conductances we used fractional derivatives of order η≤1. The fractional derivatives were used to solve the kinetic equations of each gate. We systematically classified the properties of each gate as a function of η. We then tested if the full model could generate action potentials with the different power-law behaving gates. Finally, we studied the patterns of action potential that emerged in each case. Our results show the model produces a wide range of action potential shapes and spiking patterns in response to constant current stimulation as a function of η. In comparison with the classical model, the action potential shapes for power-law behaving potassium conductance (n gate) showed a longer peak and shallow hyperpolarization; for power-law activation of the sodium conductance (m gate), the action potentials had a sharp rise time; and for power-law inactivation of the sodium conductance (h gate) the spikes had wider peak that for low values of η replicated pituitary- and cardiac-type action potentials. With all physiological parameters fixed a wide range of spiking patterns emerged as a function of the value of the constant input current and η, such as square wave bursting, mixed mode oscillations, and pseudo-plateau potentials. Our analyses show that the intrinsic memory trace of the fractional derivative provides a negative feedback mechanism between the voltage trace and the activity of the power-law behaving gate variable. As a consequence, power-law behaving conductances result in an increase in the number of spiking patterns a neuron can generate and, we propose, expand the computational capacity of the neuron.

In-depth Characterization of Silicon Nanowire Field-Effect Transistor (SiNW-FET) for Neural Recording and Direct Performance Comparison with Passive MEA

In-depth Characterization of Silicon Nanowire Field-Effect Transistor (SiNW-FET) for Neural Recording and Direct Performance Comparison with Passive MEA

Spike-count distribution in a neuronal population under weak common stimulation.

We study the probability distribution of the number of synchronous action potentials (spike count) in a model network consisting of a homogeneous neural population that is driven by a common time-dependent stimulus. We derive two analytical approximations for the count statistics, which are based on linear response theory and hold true for weak input correlations. Comparison to numerical simulations of populations of integrate-and-fire neurons in different parameter regimes reveals that our theory correctly predicts how much a weak common stimulus increases the probability of common firing and of common silence in the neural population.

Chronic Intermittent Ethanol Exposure Enhances the Excitability and Synaptic Plasticity of Lateral Orbitofrontal Cortex Neurons and Induces a Tolerance to the Acute Inhibitory Actions of Ethanol.

Alcoholism is associated with changes in brain reward and control systems, including the prefrontal cortex. In prefrontal areas, the orbitofrontal cortex (OFC) has been suggested to have an important role in the development of alcohol-abuse disorders and studies from this laboratory demonstrate that OFC-mediated behaviors are impaired in alcohol-dependent animals. However, it is not known whether chronic alcohol (ethanol) exposure alters the fundamental properties of OFC neurons. In this study, mice were exposed to repeated cycles of chronic intermittent ethanol (CIE) exposure to induce dependence and whole-cell patch-clamp electrophysiology was used to examine the effects of CIE treatment on lateral OFC (lOFC) neuron excitability, synaptic transmission, and plasticity. Repeated cycles of CIE exposure and withdrawal enhanced current-evoked action potential (AP) spiking and this was accompanied by a reduction in the after-hyperpolarization and a decrease in the functional activity of SK channels. CIE mice also showed an increase in the AMPA/NMDA ratio, and this was associated with an increase in GluA1/GluA2 AMPA receptor expression and a decrease in GluN2B NMDA receptor subunits. Following CIE treatment, lOFC neurons displayed a persistent long-term potentiation of glutamatergic synaptic transmission following a spike-timing-dependent protocol. Lastly, CIE treatment diminished the inhibitory effect of acute ethanol on AP spiking of lOFC neurons and reduced expression of the GlyT1 transporter. Taken together, these results suggest that chronic exposure to ethanol leads to enhanced intrinsic excitability and glutamatergic synaptic signaling of lOFC neurons. These alterations may contribute to the impairment of OFC-dependent behaviors in alcohol-dependent individuals.

Spike rate of multi-unit muscle sympathetic nerve fibers after catheter-based renal nerve ablation

Patients with treatment-resistant arterial hypertension exhibited profound reductions in single sympathetic vasoconstrictor fiber firing rates after renal nerve ablation. In contrast, integrated multi-unit muscle sympathetic nerve activity (MSNA) changed little or not at all. We hypothesized that conventional MSNA analysis may have missed single fiber discharges, thus, obscuring sympathetic inhibition after renal denervation. We studied patients with difficult-to-control arterial hypertension (age 45–74 years) before, 6 (n = 11), and 12 months (n = 8) after renal nerve ablation. Electrocardiogram, respiration, brachial, and finger arterial blood pressure (BP), as well as the MSNA and raw MSNA signals were analyzed. We detected MSNA action-potential spikes using 2 stage kurtosis wavelet denoising techniques to assess mean, median, and maximum spike rates for each beat-to-beat interval. Supine heart rate and systolic BP did not change at 6 (ΔHR: −2 ± 3 bpm; ΔSBP: 2 ± 9 mm Hg) or at 12 months (ΔHR: −1 ± 3 mm Hg, ΔSBP: −1 ± 9 mm Hg) after renal nerve ablation. Mean burst frequency and mean spike frequency at baseline were 34 ± 3 bursts per minute and 8 ± 1 spikes per second. Both measurements did not change at 6 months (−1.4 ± 3.6 bursts/minute; −0.6 ± 1.4 spikes/second) or at 12 months (−2.5 ± 4.0 bursts/minute; −2.0 ± 1.6 spikes/second) after renal nerve ablation. After renal nerve ablation, BP decreased in 3 of 11 patients. BP and MSNA spike frequency changes were not correlated (slope = −0.06; P = .369). Spike rate analysis of multi-unit MSNA neurograms further suggests that profound sympathetic inhibition is not a consistent finding after renal nerve ablation.

Kv3.4 channel function and dysfunction in nociceptors

Recently, we reported the isolation of the Kv3.4 current in dorsal root ganglion (DRG) neurons and described dysregulation of this current in a spinal cord injury (SCI) model of chronic pain. These studies strongly suggest that rat Kv3.4 channels are major regulators of excitability in DRG neurons from pups and adult females, where they help determine action potential (AP) repolarization and spiking properties. Here, we characterized the Kv3.4 current in rat DRG neurons from adult males and show that it transfers 40–70% of the total repolarizing charge during the AP across all ages and sexes. Following SCI, we also found remodeling of the repolarizing currents during the AP. In the light of these studies, homomeric Kv3.4 channels expressed in DRG nociceptors are emerging novel targets that may help develop new approaches to treat neuropathic pain.

Phase synchronization of coupled bursting neurons and the generalized Kuramoto model

Bursting neurons fire rapid sequences of action potential spikes followed by a quiescent period. The basic dynamical mechanism of bursting is the slow currents that modulate a fast spiking activity caused by rapid ionic currents. Minimal models of bursting neurons must include both effects. We considered one of these models and its relation with a generalized Kuramoto model, thanks to the definition of a geometrical phase for bursting and a corresponding frequency. We considered neuronal networks with different connection topologies and investigated the transition from a non-synchronized to a partially phase-synchronized state as the coupling strength is varied. The numerically determined critical coupling strength value for this transition to occur is compared with theoretical results valid for the generalized Kuramoto model.

Synaptic Conductance Estimates of the Connection Between Local Inhibitor Interneurons and Pyramidal Neurons in Layer 2/3 of a Cortical Column.

Stimulation of a principal whisker yields sparse action potential (AP) spiking in layer 2/3 (L2/3) pyramidal neurons in a cortical column of rat barrel cortex. The low AP rates in pyramidal neurons could be explained by activation of interneurons in L2/3 providing inhibition onto L2/3 pyramidal neurons. L2/3 interneurons classified as local inhibitors based on their axonal projection in the same column were reported to receive strong excitatory input from spiny neurons in L4, which are also the main source of the excitatory input to L2/3 pyramidal neurons. Here, we investigated the remaining synaptic connection in this intracolumnar microcircuit. We found strong and reliable inhibitory synaptic transmission between intracolumnar L2/3 local-inhibitor-to-L2/3 pyramidal neuron pairs [inhibitory postsynaptic potential (IPSP) amplitude −0.88 ± 0.67 mV]. On average, 6.2 ± 2 synaptic contacts were made by L2/3 local inhibitors onto L2/3 pyramidal neurons at 107 ± 64 µm path distance from the pyramidal neuron soma, thus overlapping with the distribution of synaptic contacts from L4 spiny neurons onto L2/3 pyramidal neurons (67 ± 34 µm). Finally, using compartmental simulations, we determined the synaptic conductance per synaptic contact to be 0.77 ± 0.4 nS. We conclude that the synaptic circuit from L4 to L2/3 can provide efficient shunting inhibition that is temporally and spatially aligned with the excitatory input from L4 to L2/3.

A practice of caution: spontaneous action potentials or artifactual spikes?

BackgroundHigh density surface electromyogram (EMG) techniques with electrode arrays have been used to record spontaneous muscle activity, which is important, both for supporting the diagnosis of neuromuscular diseases and for laboratory based neurophysiological investigations. This short report addresses a practical issue we have experienced during recording of spontaneous muscle activity using electrode arrays from subjects with major neuromuscular disorders.FindingsWe show that recording artifacts can appear similar to spontaneous action potential spikes. Moreover, a causal filter may induce asymmetric distortions of an artifact and thus confuse it with a real action potential spike. As a consequence, for a single channel surface EMG recording, it might be difficult to judge whether a voltage transient is a real action potential or an artifact. Further investigation of the signal distributions among other channels of the array can be used to reach a more confident judgment.ConclusionsDuring examination of spontaneous muscle activity using electrode arrays, caution is required for differentiation of physiological signals from artifactual spikes, which is important for accurate extraction of diagnostic or investigatory information.

(±)3,4-Methylenedioxyamphetamine inhibits the TEA-sensitive K+ current in the hippocampal neuron and the Kv2.1 current expressed in H1355 cells

The whole-cell patch clamp method was used to study the effects of (±)3,4-methylenedioxyamphetamine (MDA) in hippocampal CA1 neurons from neonatal rats and in lung epithelial H1355 cells expressing Kv2.1. Extracellular application of MDA (30 μM) induced bursts of action potentials in hippocampal CA1 neurons exhibiting single spike action potentials without a bursting firing pattern, and promoted action potential bursts in hippocampal CA1 neurons exhibiting bursting firing of action potentials. Whereas MDA (30 and 100 μM) markedly decreased the delayed outward current in hippocampal CA1 neurons, MDA (100 μM) only slightly decreased the fast-inactivating K+ current (IA) current. Furthermore, MDA (100 μM) substantially decreased the delayed outward current in the presence of 4-aminopyridine (4-AP; 3 mM), but did not significantly decrease the delayed outward current in the presence of tetraethylammonium (TEA; 30 mM). MDA (100 μM) also inhibited the current in H1355 cells expressing Kv2.1. Our results have shown that MDA inhibits the TEA-sensitive K+ current in the hippocampus and the Kv2.1 current expressed in H1355 cells. These effects may contribute to the pharmacological and toxicological effects of MDA.

Spike detection and sorting using PARAFAC2 method.

In this contribution we introduce the Parallel Factor 2 (PARAFAC2) analysis as a novel method for the simultaneous detection and classification of neural action potentials. In order to measure these action potentials (spike signals), stem cell derived neuronal cells are cultivated on the surface of a Micro Electrode Array (MEA). Here, the neuronal cells produce ion currents, which can be measured as extracellular electric potentials. Whenever a cell or a group of cells produces ion currents, either spontaneously or evoked by a stimulus, a spike signal can be measured by the electrodes of the MEA. Stimulated cells produce spikes and groups of spikes (bursts) which propagate in space over the MEA. In the recorded data, different source types (e.g., cells which respond directly to external stimuli and cells which are triggered by other neural cells) are characterized by different spike shapes. The proposed PARAFAC2 method is able to separate these spike shapes (sources) in time, frequency and space (channels) enabling an improved performance in noisy scenarios. Furthermore, PARAFAC2 allows for a causality analysis on the measured spike signals (i.e. the identification of different signal paths). Thereby, the PARAFAC2 decomposition is able to exploit the multi-dimensional structure of the MEA data.

Matching of feedback inhibition with excitation ensures fidelity of information flow in the anterior piriform cortex

Odor-evoked responses in mitral cells of the olfactory bulb are characterized by prolonged patterns of action potential (spike) activity. If downstream neurons are to respond to each spike in these patterns, the duration of the excitatory response to one spike should be limited, enabling cells to respond to subsequent spikes. To test for such mechanisms, we performed patch-clamp recordings in slices of the mouse anterior piriform cortex. Mitral cell axons in the lateral olfactory tract (LOT) were stimulated electrically at different intensities and with various frequency patterns to mimic changing input conditions that the piriform cortex likely encounters in vivo. We found with cell-attached measurements that superficial pyramidal (SP) cells in layer 2 consistently responded to LOT stimulation across conditions with a limited number (1–2) of spikes per stimulus pulse. The key synaptic feature accounting for the limited spike number appeared to be somatic inhibition derived from layer 3 fast-spiking cells. This inhibition tracked the timing of the first spike in SP cells across conditions, which naturally limited the spike number to 1–2. These response features to LOT stimulation were, moreover, not unique to SP cells, also occurring in a population of fluorescently labeled interneurons in glutamic acid decarboxylase 65-eGFP mice. That these different cortical cells respond to incoming inputs with 1–2 spikes per stimulus may be especially critical for relaying bulbar information contained in synchronized oscillations at beta (15–30Hz) or gamma (30–80Hz) frequencies.

Juxtasomal Biocytin Labeling to Study the Structure-function Relationship of Individual Cortical Neurons

The cerebral cortex is characterized by multiple layers and many distinct cell-types that together as a network are responsible for many higher cognitive functions including decision making, sensory-guided behavior or memory. To understand how such intricate neuronal networks perform such tasks, a crucial step is to determine the function (or electrical activity) of individual cell types within the network, preferentially when the animal is performing a relevant cognitive task. Additionally, it is equally important to determine the anatomical structure of the network and the morphological architecture of the individual neurons to allow reverse engineering the cortical network. Technical breakthroughs available today allow recording cellular activity in awake, behaving animals with the valuable option of posthoc identifying the recorded neurons. Here, we demonstrate the juxtasomal biocytin labeling technique, which involves recording action potential spiking in the extracellular (or loose-patch) configuration using conventional patch pipettes. The juxtasomal recording configuration is relatively stable and applicable across behavioral conditions, including anesthetized, sedated, awake head-fixed, and even in the freely moving animal. Thus, this method allows linking cell-type specific action potential spiking during animal behavior to reconstruction of the individual neurons and ultimately, the entire cortical microcircuit. In this video manuscript, we show how individual neurons in the juxtasomal configuration can be labeled with biocytin in the urethane-anaesthetized rat for posthoc identification and morphological reconstruction.

A Clinically Applicable Approach for Detecting Spontaneous Action Potential Spikes in Amyotrophic Lateral Sclerosis With a Linear Electrode Array

Examination of spontaneous muscle activity is an important part of the routine electromyogram (EMG) in assessing neuromuscular diseases. The EMG is specifically valuable as a diagnostic test in supporting the diagnosis of amyotrophic lateral sclerosis. High-density surface EMG is a relatively new technique that has until now been used in research but has the potential for clinical application. This study presents a simple high-density surface EMG method for automatic detection of spontaneous action potentials from surface electrode array recordings of patients with amyotrophic lateral sclerosis. To reduce computational complexity while maintaining useful information from the electrode array recording, the multichannel high-density surface EMG was transferred to single-dimensional data by calculating the maximum difference across all channels of the electrode array. A spike detection threshold was then set in the single-dimensional domain to identify the firing times of each spontaneous action potential spike, whereas a spike extraction threshold was used to define the onset and offset of the spontaneous spikes. These data were used to extract the spontaneous spike waveforms from the electrode array EMG. A database of detected spontaneous spikes was thus obtained, including their waveforms, on all channels along with their corresponding firing times. This newly developed method makes use of the information from different channels of the electrode array EMG recording. It also has the primary feature of being simple and fast in implementation, with convenient parameter adjustment and user-computer interaction. Hence, it has good possibilities for clinical application.

Glial Regulation of Extrasynaptic NMDA Receptor-Mediated Excitation of Supraoptic Nucleus Neurones During Dehydration

Magnocellular neurosecretory cells (MNCs) of the supraoptic nucleus (SON) project to the posterior pituitary gland where they release the hormones, vasopressin and oxytocin into the circulation to maintain plasma osmolality. Hormone release is proportionate to SON MNC action potential (spike) firing rate. When activated by ambient extracellular glutamate, extrasynaptic NMDA receptors (eNMDARs) mediate a tonic (persistent) depolarisation to increase the probability of action potential firing. In the present study, in vivo single-unit electrophysiological recordings were made from urethane-anaesthetised female Sprague-Dawley rats to investigate the impact of tonic eNMDAR activation on MNC activity. Water deprivation (for up to 48 h) caused an increase in the firing rate of SON MNCs that was associated with a general increase in post-spike excitability. To determine whether eNMDAR activation contributes to the increased MNC excitability during water deprivation, memantine, which preferentially blocks eNMDARs, was administered locally into the SON by microdialysis. Memantine significantly decreased the firing rate of MNCs recorded from 48-h water-deprived rats but had no effect on MNCs recorded from euhydrated rats. In the presence of the glial glutamate transporter-1 (GLT-1) blocker, dihydrokainate, memantine also reduced the MNC firing rate in euhydrated rats. Taken together, these observations suggest that GLT-1 clears extracellular glutamate to prevent the activation of eNDMARs under basal conditions and that, during dehydration, eNMDAR activation contributes to the increased firing rate of MNCs.

Intraspinal transplantation and modulation of donor neuron electrophysiological activity

Rat fetal spinal cord (FSC) tissue, naturally enriched with interneuronal progenitors, was introduced into high cervical, hemi-resection (Hx) lesions. Electrophysiological analyses were conducted to determine if such grafts exhibit physiologically-patterned neuronal activity and if stimuli which increase respiratory motor output also alter donor neuron bursting. Three months following transplantation, the bursting activity of FSC neurons and the contralateral phrenic nerve were recorded in anesthetized rats during a normoxic baseline period and brief respiratory challenges. Spontaneous neuronal activity was detected in 80% of the FSC transplants, and autocorrelation of action potential spikes revealed distinct correlogram peaks in 87% of neurons. At baseline, the average discharge frequency of graft neurons was 13.0±1.7Hz, and discharge frequency increased during a hypoxic respiratory challenge (p<0.001). Parallel studies in unanesthetized rats showed that FSC tissue recipients had larger inspiratory tidal volumes during brief hypoxic exposures (p<0.05 vs. C2Hx rats). Anatomical connectivity was explored in additional graft recipients by injecting a transsynaptic retrograde viral tracer (pseudorabies virus, PRV) directly into matured transplants. Neuronal labeling occurred throughout graft tissues and also in the host spinal cord and brainstem nuclei, including those associated with respiratory control. These results underscore the neuroplastic potential of host–graft interactions and training approaches to enhance functional integration within targeted spinal circuitry.

A finite rate of innovation algorithm for fast and accurate spike detection from two-photon calcium imaging

Objective. Inferring the times of sequences of action potentials (APs) (spike trains) from neurophysiological data is a key problem in computational neuroscience. The detection of APs from two-photon imaging of calcium signals offers certain advantages over traditional electrophysiological approaches, as up to thousands of spatially and immunohistochemically defined neurons can be recorded simultaneously. However, due to noise, dye buffering and the limited sampling rates in common microscopy configurations, accurate detection of APs from calcium time series has proved to be a difficult problem. Approach. Here we introduce a novel approach to the problem making use of finite rate of innovation (FRI) theory (Vetterli et al 2002 IEEE Trans. Signal Process. 50 1417–28). For calcium transients well fit by a single exponential, the problem is reduced to reconstructing a stream of decaying exponentials. Signals made of a combination of exponentially decaying functions with different onset times are a subclass of FRI signals, for which much theory has recently been developed by the signal processing community. Main results. We demonstrate for the first time the use of FRI theory to retrieve the timing of APs from calcium transient time series. The final algorithm is fast, non-iterative and parallelizable. Spike inference can be performed in real-time for a population of neurons and does not require any training phase or learning to initialize parameters. Significance. The algorithm has been tested with both real data (obtained by simultaneous electrophysiology and multiphoton imaging of calcium signals in cerebellar Purkinje cell dendrites), and surrogate data, and outperforms several recently proposed methods for spike train inference from calcium imaging data.