Spike train analysis

"The great mixing machine": multisensory integration and brain-breath coupling in the cerebral cortex

It is common to distinguish between "holist" and "reductionist" views of brain function, where the former envisions the brain as functioning as an indivisible unit and the latter as a collection of distinct units that serve different functions. Opposing reductionism, a number of researchers have pointed out that cortical network architecture does not respect functional boundaries, and the neuroanatomist V. Braitenberg proposed to understand the cerebral cortex as a "great mixing machine" of neuronal activity from sensory inputs, motor commands, and intrinsically generated processes. In this paper, we offer a contextualization of Braitenberg's point, and we review evidence for the interactions of neuronal activity from multiple sensory inputs and intrinsic neuronal processes in the cerebral cortex. We focus on new insights from studies on audiovisual interactions and on the influence of respiration on brain functions, which do not seem to align well with "reductionist" views of areal functional boundaries. Instead, they indicate that functional boundaries are fuzzy and context dependent. In addition, we discuss the relevance of the influence of sensory, proprioceptive, and interoceptive signals on cortical activity for understanding brain-body interactions, highlight some of the consequences of these new insights for debates on embodied cognition, and offer some suggestions for future studies.

Study on the mechanism of high-frequency stimulation inhibiting low-Mg2+-induced epileptiform discharges in juvenile rat hippocampal slices

Study on the mechanism of high-frequency stimulation inhibiting low-Mg²⁺-induced epileptiform discharges in juvenile rat hippocampal slices High-frequency stimulation (HFS) has been demonstrated to be an effective treatment for inhibiting epilepsy in some clinical and laboratory studies. However, the mechanisms underlying the therapeutic effects of HFS are not yet fully understood. In our present study, epileptiform discharges (EDs) in acutely isolated hippocampal slices of male Sprague-Dawley (SD) juvenile rats induced by low-Mg²⁺ artificial cerebrospinal fluid (ACSF), and electrical stimulation (square wave, 900 pulses, 50 % duty-cycle, 130 Hz) was performed on the CA3 using concentric bipolar electrodes. EDs of neurons in hippocampal were recorded by multi-electrode arrays (MEA). After stable EDs events had been recorded for at least 20 min, HFS was added, followed by 10 μmol/L gamma-aminobutyric acid type A (GABA_A) receptors blocker bicuculline (BIC). The results show that the HFS can increase the discharges frequency of inter-ictal discharges (IIDs) and decrease the duration of ictal discharges (IDs). However, the HFS had no effect on the slices with 10 μmol/L BIC. These results indicated that the GABA_A receptors are activated when HFS inhibited EDs, thereby achieving the inhibition of low-Mg²⁺-induced EDs in slices.

Subcortical neuronal ensembles: an analysis of motor task association, tremor, oscillations, and synchrony in human patients

Deep brain stimulation (DBS) has expanded as an effective treatment for motor disorders, providing a valuable opportunity for intraoperative recording of the spiking activity of subcortical neurons. The properties of these neurons and their potential utility in neuroprosthetic applications are not completely understood. During DBS surgeries in 25 human patients with either essential tremor or Parkinson's disease, we acutely recorded the single-unit activity of 274 ventral intermediate/ventral oralis posterior motor thalamus (Vim/Vop) neurons and 123 subthalamic nucleus (STN) neurons. These subcortical neuronal ensembles (up to 23 neurons sampled simultaneously) were recorded while the patients performed a target-tracking motor task using a cursor controlled by a haptic glove. We observed that modulations in firing rate of a substantial number of neurons in both Vim/Vop and STN represented target onset, movement onset/direction, and hand tremor. Neurons in both areas exhibited rhythmic oscillations and pairwise synchrony. Notably, all tremor-associated neurons exhibited synchrony within the ensemble. The data further indicate that oscillatory (likely pathological) neurons and behaviorally tuned neurons are not distinct but rather form overlapping sets. Whereas previous studies have reported a linear relationship between power spectra of neuronal oscillations and hand tremor, we report a nonlinear relationship suggestive of complex encoding schemes. Even in the presence of this pathological activity, linear models were able to extract motor parameters from ensemble discharges. Based on these findings, we propose that chronic multielectrode recordings from Vim/Vop and STN could prove useful for further studying, monitoring, and even treating motor disorders.

Reprogramming movements: extraction of motor intentions from cortical ensemble activity when movement goals change

The ability to inhibit unwanted movements and change motor plans is essential for behaviors of advanced organisms. The neural mechanisms by which the primate motor system rejects undesired actions have received much attention during the last decade, but it is not well understood how this neural function could be utilized to improve the efficiency of brain-machine interfaces (BMIs). Here we employed linear discriminant analysis (LDA) and a Wiener filter to extract motor plan transitions from the activity of ensembles of sensorimotor cortex neurons. Two rhesus monkeys, chronically implanted with multielectrode arrays in primary motor (M1) and primary sensory (S1) cortices, were overtrained to produce reaching movements with a joystick toward visual targets upon their presentation. Then, the behavioral task was modified to include a distracting target that flashed for 50, 150, or 250 ms (25% of trials each) followed by the true target that appeared at a different screen location. In the remaining 25% of trials, the initial target stayed on the screen and was the target to be approached. M1 and S1 neuronal activity represented both the true and distracting targets, even for the shortest duration of the distracting event. This dual representation persisted both when the monkey initiated movements toward the distracting target and then made corrections and when they moved directly toward the second, true target. The Wiener filter effectively decoded the location of the true target, whereas the LDA classifier extracted the location of both targets from ensembles of 50–250 neurons. Based on these results, we suggest developing real-time BMIs that inhibit unwanted movements represented by brain activity while enacting the desired motor outcome concomitantly.

Analysis of motoneuron responses to composite synaptic volleys (computer simulation study)

This paper deals with the analysis of changes in motoneuron (MN) firing evoked by repetitively applied stimuli aimed toward extracting information about the underlying synaptic volleys. Spike trains were obtained from computer simulations based on a threshold-crossing model of tonically firing MN, subjected to stimulation producing postsynaptic potentials (PSPs) of various parameters. These trains were analyzed as experimental results, using the output measures that were previously shown to be most effective for this purpose: peristimulus time histogram, raster plot and peristimulus time intervalgram. The analysis started from the effects of single excitatory and inhibitory PSPs (EPSPs and IPSPs). The conclusions drawn from this analysis allowed the explanation of the results of more complex synaptic volleys, i.e., combinations of EPSPs and IPSPs, and the formulation of directions for decoding the results of human neurophysiological experiments in which the responses of tonically firing MNs to nerve stimulation are analyzed.

Cortical correlates of fitts' law

Fitts’ law describes the fundamental trade-off between movement accuracy and speed: it states that the duration of reaching movements is a function of target size (TS) and distance. While Fitts’ law has been extensively studied in ergonomics and has guided the design of human–computer interfaces, there have been few studies on its neuronal correlates. To elucidate sensorimotor cortical activity underlying Fitts’ law, we implanted two monkeys with multielectrode arrays in the primary motor (M1) and primary somatosensory (S1) cortices. The monkeys performed reaches with a joystick-controlled cursor toward targets of different size. The reaction time (RT), movement time, and movement velocity changed with TS, and M1 and S1 activity reflected these changes. Moreover, modifications of cortical activity could not be explained by changes of movement parameters alone, but required TS as an additional parameter. Neuronal representation of TS was especially prominent during the early RT period where it influenced the slope of the firing rate rise preceding movement initiation. During the movement period, cortical activity was correlated with movement velocity. Neural decoders were applied to simultaneously decode TS and motor parameters from cortical modulations. We suggest that sensorimotor cortex activity reflects the characteristics of both the movement and the target. Classifiers that extract these parameters from cortical ensembles could improve neuroprosthetic control.

Voltage fluctuations in neurons: signal or noise?

Noise and variability are fundamental companions to ion channels and synapses and thus inescapable elements of brain function. The overriding unresolved issue is to what extent noise distorts and limits signaling on one hand and at the same time constitutes a crucial and fundamental enrichment that allows and facilitates complex adaptive behavior in an unpredictable world. Here we review the growing experimental evidence that functional network activity is associated with intense fluctuations in membrane potential and spike timing. We trace origins and consequences of noise and variability. Finally, we discuss noise-free neuronal signaling and detrimental and beneficial forms of noise in large-scale functional neural networks. Evidence that noise and variability in some cases go hand in hand with behavioral variability and increase behavioral choice, richness, and adaptability opens new avenues for future studies.

iRaster: a novel information visualization tool to explore spatiotemporal patterns in multiple spike trains

Over the last few years, simultaneous recordings of multiple spike trains have become widely used by neuroscientists. Therefore, it is important to develop new tools for analysing multiple spike trains in order to gain new insight into the function of neural systems. This paper describes how techniques from the field of visual analytics can be used to reveal specific patterns of neural activity. An interactive raster plot called iRaster has been developed. This software incorporates a selection of statistical procedures for visualization and flexible manipulations with multiple spike trains. For example, there are several procedures for the re-ordering of spike trains which can be used to unmask activity propagation, spiking synchronization, and many other important features of multiple spike train activity. Additionally, iRaster includes a rate representation of neural activity, a combined representation of rate and spikes, spike train removal and time interval removal. Furthermore, it provides multiple coordinated views, time and spike train zooming windows, a fisheye lens distortion, and dissemination facilities. iRaster is a user friendly, interactive, flexible tool which supports a broad range of visual representations. This tool has been successfully used to analyse both synthetic and experimentally recorded datasets. In this paper, the main features of iRaster are described and its performance and effectiveness are demonstrated using various types of data including experimental multi-electrode array recordings from the ganglion cell layer in mouse retina. iRaster is part of an ongoing research project called VISA (Visualization of Inter-Spike Associations) at the Visualization Lab in the University of Plymouth. The overall aim of the VISA project is to provide neuroscientists with the ability to freely explore and analyse their data. The software is freely available from the Visualization Lab website (see www.plymouth.ac.uk/infovis).

Computer simulation study of the relationship between the profile of excitatory postsynaptic potential and stimulus-correlated motoneuron firing

This paper shows the results of computer simulation of changes in motoneuron (MN) firing evoked by a repetitively applied synaptic volley that consists of a single excitatory postsynaptic potential (EPSP). Spike trains produced by the threshold-crossing MN model were analyzed as experimental results. Various output functions were applied for analysis; the most useful was a peristimulus time histogram, a special modification of a raster plot and a peristimulus time frequencygram (PSTF). It has been shown that all functions complement each other in distinguishing between the genuine results evoked by the excitatory volley and the secondary results of the EPSP-evoked synchronization. The EPSP rising edge was best reproduced by the PSTF. However, whereas the EPSP rise time could be estimated quite accurately, especially for high EPSP amplitudes at high MN firing rates, the EPSP amplitude estimate was also influenced by factors unrelated to the synaptic volley, such as the afterhyperpolarization duration of the MN or the amplitude of synaptic noise, which cannot be directly assessed in human experiments. Thus, the attempts to scale any estimate of the EPSP amplitude in millivolts appear to be useless. The decaying phase of the EPSP cannot be reproduced accurately by any of the functions. For the short EPSPs, it is extinguished by the generation of an action potential and a subsequent decrease in the MN excitability. For longer EPSPs, it is inseparable from the secondary effects of synchronization. Thus, the methods aimed at extracting information about long-lasting and complex postsynaptic potentials from stimulus-correlated MN firing, should be refined, and the theoretical considerations checked in computer simulations.

Interspike interval analysis in a patient with peripheral nerve hyperexcitability and potassium channel antibodies

Neuromyotonia or Isaacs' syndrome is a rare peripheral nerve hyperexcitability disorder caused by antibodies against potassium channels of myelinated axons. We present the high-density surface electromyographic (EMG) recordings of a patient with fasciculations and cramps due to neuromyotonia. To characterize the time course of hyperexcitability, we analyzed the interspike intervals (ISIs) between fasciculation potentials, doublet, and multiplet discharges. ISI duration increased within each burst. The ISI histograms found can be explained by the recovery cycle of the myelinated axon and its dependency on the slow potassium conductance. We conclude that ISI analysis is a useful tool to understand the membrane dynamics underlying abnormal motor unit activity.

Discrete stimulus estimation from neural responses in the turtle retina

In this paper, we investigate the decoding of flashed, full-field visual stimuli while recording from a population of retinal ganglion cells. We present a direct statistical method for determining the likelihood that a response was evoked by a particular stimulus, and use this method to estimate stimuli based on microelectrode array recordings in the turtle retina. This method uses the well-known time-varying Poisson model of neural firing, along with extensions to accommodate neural refractory periods. Unlike other approaches commonly used for Poisson processes, the specific formulation presented here is bin free and requires few user-specified parameters. Statistical dependency issues and the effects of stationarity on the estimation method are also discussed.

Comparing neuronal spike trains with inhomogeneous Poisson distribution: evaluation procedure and experimental application in cases of cyclic activity

A procedure is proposed to compare single-unit spiking activity elicited in repetitive cycles with an inhomogeneous Poisson process (IPP). Each spike sequence in a cycle is discretized and represented as a point process on a circle. The interspike interval probability density predicted for an IPP is computed on the basis of the experimental firing probability density; differences from the experimental interval distribution are assessed. This procedure was applied to spike trains which were repetitively induced by opening-closing movements of the distal article of a lobster leg. As expected, the density of short interspike intervals, less than 20-40 ms in length, was found to lie greatly below the level predicted for an IPP, reflecting the occurrence of the refractory period. Conversely, longer intervals, ranging from 20-40 to 100-120 ms, were markedly more abundant than expected; this provided evidence for a time window of increased tendency to fire again after a spike. Less consistently, a weak depression of spike generation was observed for longer intervals. A Monte Carlo procedure, implemented for comparison, produced quite similar results, but was slightly less precise and more demanding as concerns computation time.

The correlation grid: analysis of synchronous spiking in multi-dimensional spike train data and identification of feasible connection architectures

This paper presents a visualization technique specifically designed to support the analysis of synchronous firings in multiple, simultaneously recorded, spike trains. This technique, called the correlation grid, enables investigators to identify groups of spike trains, where each pair of spike trains has a high probability of generating spikes approximately simultaneously or within a constant time shift. Moreover, the correlation grid was developed to help solve the following reverse problem: identification of the connection architecture between spike train generating units, which may produce a spike train dataset similar to the one under analysis. To demonstrate the efficacy of this approach, results are presented from a study of three simulated, noisy, spike train datasets. The parameters of the simulated neurons were chosen to reflect the typical characteristics of cortical pyramidal neurons. The schemes of neuronal connections were not known to the analysts. Nevertheless, the correlation grid enabled the analysts to find the correct connection architecture for each of these three data sets.

Strength and voluntary activation of the quadriceps femoris muscle at different severities of osteoarthritic knee joint damage

OBJECTIVE

Improvements of quadriceps motor deficits represent a major therapeutical target in knee osteoarthritis (OA). In the present study, we investigated changes in quadriceps function at different stages of osteoarthritic cartilage damage.

METHODS

Measurements of quadriceps voluntary activation (VA) and maximum voluntary contraction (MVC) were performed by a twitch interpolation technique and the total muscular capacity (TMC) was calculated as the ratio of MVC and VA. We assessed 68 patients (56.7+/-9.5 years) with stage II and 154 patients (65.6+/-6.0 years) with stage IV chondropathy. As controls, we used 85 age related healthy subjects (58.1+/-8.7 years).

RESULTS

While TMC was significantly lower in stage IV (90.6+/-43.7 Nm) than in stage II chondropathy (109.6+/-51.0) there were no differences in the MVCs between both groups. Quadriceps VA was even higher in stage IV (77.2+/-13.2%) than in stage II chondropathy (70.8+/-16.0%). In the controls, MVC, VA and TMC were significantly higher than in both OA groups.

CONCLUSION

We assume that a decrease of TMC might occur within the course of OA and, in consequence, VA increases to maintain quadriceps MVC.

Kinesthetic sense of the shoulder in patients with impingement syndrome

A proprioceptive deficit is an important determinant of disability in various shoulder disorders, such as instability and osteoarthrosis. In 15 patients with impingement syndrome stage II (Neer 1983), who were treated by arthroscopic subacromial decompression, we measured movement sense by determining threshold levels for the perception of motion of the shoulder. The patients were placed in a specially designed chair allowing continuous passive motion of the shoulder joint, while avoiding cutaneous, auditory and visual stimuli. To assess movement detection thresholds, passive abduction movements of the shoulder were performed at a starting angle of 60 degrees, an amplitude of 10 degrees and an angular velocity of 1.3 degrees/s. Before surgery, all patients had higher threshold levels for the perception of motion in their affected shoulders then in the other side. After decompression, proprioception had improved on the decompressed side, but was unchanged on the other side.

Weakly modulated spike trains: significance, precision, and correction for sample size

Many single-unit electrophysiological studies of visual cortex have investigated strong evoked responses to simple stimuli such as oriented gratings. Experiments involving other types of stimuli, such as natural scenes, higher-order features, and surface brightness, produce single-unit responses that are more difficult to interpret. Experiments with brightness, in particular, evoke single-unit responses that are typically weakly modulated. When the brightness is generated by a visual illusion such as the Cornsweet illusion, statistical tests are often necessary to distinguish true responses from baseline fluctuations. Here, using data collected from cat Areas 17 and 18 in response to real and illusory brightness stimuli, we provide a method for detecting and quantifying weak but significant periodic responses. By randomizing spike trains (via bootstrap methods), we provide confidence levels for response significance, permitting the evaluation of both weak and strong responses. We show that because of a strong dependence on total spike number, response significance can only be appropriately determined with randomized spike trains of similar spike number. Such randomizations can be performed for both stimulus-elicited and spontaneously occurring spike trains. By developing a method for generating randomized modulated spike trains (phase-restricted randomization) from actual recordings, we calculate upper and lower confidence limits of modulated spike trains and describe how measurement precision varies as a function of total spike count. Finally, using this randomization method, we describe how a correction function can be determined to correct for measurement bias introduced at low spike counts. These methods may also be useful in the study of small but potentially significant responses in other systems.

Evaluation of quadriceps strength and voluntary activation after unicompartmental arthroplasty for medial osteoarthritis of the knee

INTRODUCTION

In early and moderate stages of osteoarthritis (OA) of the knee, arthrogenous muscle inhibition (AMI) is an important factor for the initiation and the progression of the disease. Although AMI has been shown to be reduced after physiotherapeutical exercises resulting in significant improvements in disability, implantation of unicondylar knee arthroplasties is much provided in these stages of OA. Therefore, in the present study we investigate changes in quadriceps muscle after implantation of such prostheses as compared to physiotherapeutical treatment, alone.

METHODS

In eighteen patients with bilateral moderate knee OA, who were treated with unicondylar knee arthroplasty we investigated voluntary activation (VA) and maximum voluntary contraction (MVC) of the quadriceps femoris muscle. There were 7 males and 11 females, the mean age at time of operation was 67 years (range 58-76 years). Measurements on both sides were performed preoperatively and 18 months postoperatively using the twitch-interpolation technique.

RESULTS

Follow-up assessment revealed a significant VA and MVC increase in both the surgically treated knees and in the contralateral knees treated by physiotherapy alone. However, VA and MVC improvements were significantly higher in the operated on knees than in those treated by physiotherapy alone.

DISCUSSION

Both physiotherapeutical exercise and unicondylar knee replacements lead to an improvement of quadriceps motor function in knee OA. The greater improvement in knees with both knee replacement and physiotherapy might be related to the intraoperative removal of arthritic tissue in these knees.

Magnetic stimulation-induced modulations of motor unit firings extracted from multi-channel surface EMG

The noninvasive assessment of motor unit (MU) firing patterns on the basis of topographical information from 128-channel high-density surface electromyography (SEMG) is reported. First, multi-channel MU action potential (MUAP) templates are obtained by clustering detected firing events according to the surface topography of the MUAP. Second, a template-matching algorithm is used to find all firings of a MU, including the superimpositions of MUAPs. From a single recording, the firing pattern of up to five MUs could be derived. The modulation of MU firing by transcranial magnetic stimulation was analyzed in peri-stimulus time histograms. The results are similar to previous results of transcranial magnetic stimulation (TMS) obtained by needle electromyographic (EMG) recordings. The method can be used to investigate MU firing patterns in patients with central motor disorders. An additional advantage of the technique, apart from its noninvasiveness, is the structural and functional information that it provides on the MUs, which is not obtained by needle EMG.

Effects of large excitatory and inhibitory inputs on motoneuron discharge rate and probability

We elicited repetitive discharge in hypoglossal motoneurons recorded in slices of rat brain stem using a combination of a suprathreshold injected current step with superimposed noise to mimic the synaptic drive likely to occur during physiological activation. The effects of repetitive en mass stimulation of afferent nerves were simulated by the further addition of trains of injected current transients of varying shapes and sizes. The effects of a given current transient on motoneuron discharge timing and discharge rate were measured by calculating a peristimulus time histogram (PSTH) and a peristimulus frequencygram (PSF). The amplitude and time course of the simulated postsynaptic potentials (PSPs) produced by the current transients were calculated by convolving the current transient with an estimate of the passive impulse response of the motoneuron. We then compared the shape of the injected current transient and the simulated PSP to the profiles of the PSTH and the PSF records. The PSTHs produced by excitatory PSPs (EPSPs) were characterized by a large, short-latency increase in firing probability that lasted slightly longer than the rising phase of the EPSP, followed by a reduced discharge probability during the falling phase of the EPSP. In contrast, the PSF analysis revealed a proportionate increase in discharge rate over the entire profile of the EPSP, even though relatively few spikes occurred during the falling phase. The PSTHs associated with inhibitory PSPs (IPSPs) indicated a reduction in discharge probability during the initial, hyperpolarizing phase of the IPSP, followed by an increase in the discharge probability during its subsequent repolarizing phase. Using the PSF analysis, the initial phase of the IPSP appeared as a large hole in the record where a very small number or no discharges occurred. The subsequent phase of the IPSP was associated with frequency values that were lower than the background values. The primary features of both PSTHs and PSFs can be used to estimate the relative amplitudes of the underlying EPSPs and IPSPs. However, PSTHs contain secondary peaks and troughs that are not directly related to the underlying PSP but instead reflect the regular recurrence of spikes following those affected by the PSP. The PSF analysis is more useful for indicating the total duration and the profile of the underlying PSP. The shape of the underlying PSP can be obtained directly from the PSF records because the discharge frequency of the spikes follow the PSPs very closely, especially for EPSPs.

Assessing changes in presynaptic inhibition of Ia afferents during movement in humans

Different methods, based on different principles, have been proposed to estimate changes in presynaptic inhibition of Ia terminals (accompanied by primary afferent depolarization, (PAD)) during voluntary contraction in humans. (i) A discrepancy between the H-reflex amplitude, at an equal level of EMG activity, in two situations (e.g., walking and standing) may be taken as suggesting a different control of PAD interneurones in the two cases. (ii) A conditioning stimulation (vibration or electrical stimulation) is used to activate PAD interneurones and to evoke presynaptic inhibition of the afferent volley of the test reflex. The resulting long-lasting depression of the reflex depends on the excitability of PAD interneurones, but can be contaminated by long-lasting post-synaptic effects. (iii) The amount of reflex facilitation evoked by a purely monosynaptic Ia volley varies inversely with the on-going presynaptic inhibition of Ia afferents mediating the conditioning volley, and can be used to assess this on-going presynaptic inhibition. None of these methods can provide by itself unequivocal evidence for a change in presynaptic inhibition of Ia terminals, but reasonably reliable interpretations may be proposed when congruent results are obtained with different methods. Thus it has been shown that, during selective voluntary contraction, presynaptic inhibition is decreased on Ia afferents projecting on motoneurones of the contracting muscle and increased on Ia afferents projecting on motor nuclei not involved in the contraction.