Cross-correlation assessment of synaptic strength of single Ia fibre connections with triceps surae motoneurones in cats

Extensive Cortical Convergence to Primate Reticulospinal Pathways

Early evolution of the motor cortex included development of connections to brainstem reticulospinal neurons; these projections persist in primates. In this study, we examined the organization of corticoreticular connections in five macaque monkeys (one male) using both intracellular and extracellular recordings from reticular formation neurons, including identified reticulospinal cells. Synaptic responses to stimulation of different parts of primary motor cortex (M1) and supplementary motor area (SMA) bilaterally were assessed. Widespread short latency excitation, compatible with monosynaptic transmission over fast-conducting pathways, was observed, as well as longer latency responses likely reflecting a mixture of slower monosynaptic and oligosynaptic pathways. There was a high degree of convergence: 56% of reticulospinal cells with input from M1 received projections from M1 in both hemispheres; for SMA, the equivalent figure was even higher (70%). Of reticulospinal neurons with input from the cortex, 78% received projections from both M1 and SMA (regardless of hemisphere); 83% of reticulospinal cells with input from M1 received projections from more than one of the tested M1 sites. This convergence at the single cell level allows reticulospinal neurons to integrate information from across the motor areas of the cortex, taking account of the bilateral motor context. Reticulospinal connections are known to strengthen following damage to the corticospinal tract, such as after stroke, partially contributing to functional recovery. Extensive corticoreticular convergence provides redundancy of control, which may allow the cortex to continue to exploit this descending pathway even after damage to one area.SIGNIFICANCE STATEMENT The reticulospinal tract (RST) provides a parallel pathway for motor control in primates, alongside the more sophisticated corticospinal system. We found extensive convergent inputs to primate reticulospinal cells from primary and supplementary motor cortex bilaterally. These redundant connections could maintain transmission of voluntary commands to the spinal cord after damage (e.g., after stroke or spinal cord injury), possibly assisting recovery of function.

Plastic changes in primate motor cortex following paired peripheral nerve stimulation

Repeated paired stimulation of two peripheral nerves can produce lasting changes in motor cortical excitability, but little is known of the underlying neuronal basis. Here, we trained two macaque monkeys to perform selective thumb and index finger abduction movements. Neural activity was recorded from the contralateral primary motor cortex during task performance, and following stimulation of the ulnar and median nerves, and the nerve supplying the extensor digitorum communis (EDC) muscle. Responses were compared before and after 1 h of synchronous or asynchronous paired ulnar/median nerve stimulation. Task performance was significantly enhanced after asynchronous and impaired after synchronous stimulation. The amplitude of short latency neural responses to median and ulnar nerve stimulation was increased after asynchronous stimulation; later components were reduced after synchronous stimulation. Synchronous stimulation increased neural activity during thumb movement and decreased it during index finger movement; asynchronous stimulation decreased activity during both movements. To assess how well neural activity could separate behavioral or sensory conditions, linear discriminant analysis was used to decode which nerve was stimulated, or which digit moved. Decoding accuracy for nerve stimulation was decreased after synchronous and increased after asynchronous paired stimulation. Decoding accuracy for task performance was decreased after synchronous but was unchanged after asynchronous paired stimulation. Paired stimulation produces changes in motor cortical circuits that outlast the stimulation. Some of these changes depend on precise stimulus timing.

NEW & NOTEWORTHY Paired stimulation of peripheral nerves for 1 h induced lasting changes in neural responses within the motor cortex to nerve stimulation and to performance of a behavioral task. These changes were sufficient to alter the efficiency with which activity could encode stimulus type. Stimuli that can be easily applied noninvasively in human subjects can alter central motor circuits.

Quantitative input-output relationships between human soleus muscle spindle afferents and motoneurons

A method is described that, for the first time, allows instantaneous estimation of the Ia fiber input to human soleus motoneurons following electrical stimulation of the tibial nerve. The basis of the method is to determine the thresholds of the most and least excitable 1a fibers to electrical stimulation, and to treat the intervening thresholds as having a normal distribution about the mean; the validity of this approach is discussed. It was found that, for the same Ia fiber input, the percentage of soleus motoneurons contributing to the H (Hoffmann)-reflex differed considerably among subjects; when the results were pooled, however, there was an approximately linear relationship between Ia input and motoneuron output. Weak extension of the great toe diminished the soleus motoneuron reflex discharge in all but 2 of 16 subjects; the results for weak ankle plantarflexion were less consistent, but overall, there was a reduction in soleus motoneuron output also. The methodology should provide new insights into disorders of movement and tone, especially as it permits estimates of motoneuron depolarization to be made. NEW & NOTEWORTHY Assuming a normal distribution of Ia fiber thresholds to electrical stimulation and using the H-reflex, we determined for the first time an Ia input-α-motoneuron output relationship for the human soleus muscle. The relationship varies greatly among subjects but, overall, is approximately linear. Minimal contraction of a toe muscle alters the relationship dramatically, probably due to presynaptic inhibition of Ia fibers. Drawing on the literature, we can calculate changes in α-motoneuron membrane potential.

Corticospinal gating during action preparation and movement in the primate motor cortex

During everyday actions there is a need to be able to withhold movements until the most appropriate time. This motor inhibition is likely to rely on multiple cortical and subcortical areas, but the primary motor cortex (M1) is a critical component of this process. However, the mechanisms behind this inhibition are unclear, particularly the role of the corticospinal system, which is most often associated with driving muscles and movement. To address this, recordings were made from identified corticospinal (PTN, n = 94) and corticomotoneuronal (CM, n = 16) cells from M1 during an instructed delay reach-to-grasp task. The task involved the animals withholding action for ~2 s until a GO cue, after which they were allowed to reach and perform the task for a food reward. Analysis of the firing of cells in M1 during the delay period revealed that, as a population, non-CM PTNs showed significant suppression in their activity during the cue and instructed delay periods, while CM cells instead showed a facilitation during the preparatory delay. Analysis of cell activity during movement also revealed that a substantial minority of PTNs (27%) showed suppressed activity during movement, a response pattern more suited to cells involved in withholding rather than driving movement. These results demonstrate the potential contributions of the M1 corticospinal system to withholding of actions and highlight that suppression of activity in M1 during movement preparation is not evenly distributed across different neural populations.

NEW & NOTEWORTHY Recordings were made from identified corticospinal (PTN) and corticomotoneuronal (CM) cells during an instructed delay task. Activity of PTNs as a population was suppressed during the delay, in contrast to CM cells, which were facilitated. A minority of PTNs showed a rate profile that might be expected from inhibitory cells and could suggest that they play an active role in action suppression, most likely through downstream inhibitory circuits.

Different phase delays of peripheral input to primate motor cortex and spinal cord promote cancellation at physiological tremor frequencies

Neurons in the spinal cord and motor cortex (M1) are partially phase-locked to cycles of physiological tremor, but with opposite phases. Convergence of spinal and cortical activity onto motoneurons may thus produce phase cancellation and a reduction in tremor amplitude. The mechanisms underlying this phase difference are unknown. We investigated coherence between spinal and M1 activity with sensory input. In two anesthetized monkeys, we electrically stimulated the medial, ulnar, deep radial, and superficial radial nerves; stimuli were timed as independent Poisson processes (rate 10 Hz). Single units were recorded from M1 (147 cells) or cervical spinal cord (61 cells). Ninety M1 cells were antidromically identified as pyramidal tract neurons (PTNs); M1 neurons were additionally classified according to M1 subdivision (rostral/caudal, M1r/c). Spike-stimulus coherence analysis revealed significant coupling over a broad range of frequencies, with the strongest coherence at <50 Hz. Delays implied by the slope of the coherence phase-frequency relationship were greater than the response onset latency, reflecting the importance of late response components for the transmission of oscillatory inputs. The spike-stimulus coherence phase over the 6–13 Hz physiological tremor band differed significantly between M1 and spinal cells (phase differences relative to the cord of 2.72 ± 0.29 and 1.72 ± 0.37 radians for PTNs from M1c and M1r, respectively). We conclude that different phases of the response to peripheral input could partially underlie antiphase M1 and spinal cord activity during motor behavior. The coordinated action of spinal and cortical feedback will act to reduce tremulous oscillations, possibly improving the overall stability and precision of motor control.

THE USE OF NONPARAMETRIC METHODS OF STATIONARY POINT PROCESSES IN THE STUDY OF COMPLEX INTERACTIONS IN THE NEUROMUSCULAR SYSTEM

In this paper we study the complex interactions involved in the incoming stimulus, from a gamma (γ) and/or an alpha (α) motoneuron, and the outgoing response from the muscle spindle transmitted by the Ia sensory afferent neuron to the spinal cord. The most interesting case is the γ and α coactivation to the function of the muscle spindle, while the effect from a single (γ or α) motoneuron is also presented as a comparison. The mathematical background of this analysis is based on the theory of stationary point processes. A kernel type method of estimating second- and third-order conditional densities is used. Under certain conditions the asymptotic distributions of these estimates are Normal and the construction of 95% approximate confidence intervals is feasible. The application of these asymptotic results to the system of the muscle spindle enables us to detect and interpret its excitatory and/or inhibitory behavior.

Assessment of across-muscle coherence using multi-unit vs. single-unit recordings

Coherence between electromyographic (EMG) signals has been used to identify correlated neural inputs to motor units (MUs) innervating different muscles. Simulations using a motor-unit model (Fuglevand et al. 1992) were performed to determine the ability of coherence between two multi-unit EMGs (mEMG) to detect correlated MU activity and the range of correlation strengths in which mEMG coherence can be usefully employed. Coherence between motor-unit and mEMG activities in two muscles was determined as we varied the strength of a 30-Hz periodic common input, the number of correlated MU pairs and variability of MU discharge relative to the common input. Pooled and mEMG coherence amplitudes positively and negatively accelerated, respectively, toward the strongest and most widespread correlating inputs. Furthermore, the relation between pooled and mEMG coherence was also nonlinear and was essentially the same whether correlation strength varied by changing common input strength or its distribution. However, the most important finding is that while the mEMG coherence saturates at the strongest common input strengths, this occurs at common input strengths greater than found in most physiological studies. Thus, we conclude that mEMG coherence would be a useful measure in many experimental conditions and our simulation results suggest further guidelines for using and interpreting coherence between mEMG signals.

Permanent reorganization of Ia afferent synapses on motoneurons after peripheral nerve injuries

After peripheral nerve injuries to a motor nerve, the axons of motoneurons and proprioceptors are disconnected from the periphery and monosynaptic connections from group I afferents and motoneurons become diminished in the spinal cord. Following successful reinnervation in the periphery, motor strength, proprioceptive sensory encoding, and Ia afferent synaptic transmission on motoneurons partially recover. Muscle stretch reflexes, however, never recover and motor behaviors remain uncoordinated. In this review, we summarize recent findings that suggest that lingering motor dysfunction might be in part related to decreased connectivity of Ia afferents centrally. First, sensory afferent synapses retract from lamina IX, causing a permanent relocation of the inputs to more distal locations and significant disconnection from motoneurons. Second, peripheral reconnection between proprioceptive afferents and muscle spindles is imperfect. As a result, a proportion of sensory afferents that retain central connections with motoneurons might not reconnect appropriately in the periphery. A hypothetical model is proposed in which the combined effect of peripheral and central reconnection deficits might explain the failure of muscle stretch to initiate or modulate firing of many homonymous motoneurons.

Deciphering the contribution of intrinsic and synaptic currents to the effects of transient synaptic inputs on human motor unit discharge

The amplitude and time course of synaptic potentials in human motoneurons can be estimated in tonically discharging motor units by measuring stimulus-evoked changes in the rate and probability of motor unit action potentials. However, in spite of the fact that some of these techniques have been used for over 30 years, there is still no consensus on the best way to estimate the characteristics of synaptic potentials or on the accuracy of these estimates. In this review, we compare different techniques for estimating synaptic potentials from human motor unit discharge and also discuss relevant animal models in which estimated synaptic potentials can be compared to those directly measured from intracellular recordings. We also review the experimental evidence on how synaptic noise and intrinsic motoneuron properties influence their responses to synaptic inputs. Finally, we consider to what extent recordings of single motor unit discharge in humans can be used to distinguish the contribution of changes in synaptic inputs versus changes in intrinsic motoneuron properties to altered motoneuron responses following CNS injury.

Synaptic interactions between forelimb-related motor cortex neurons in behaving primates

We investigated the synaptic interactions between neighboring motor cortex cells in monkeys generating isometric ramp-and-hold torques about the wrist. For pairs of cortical cells the response patterns were determined in response-aligned averages and their synaptic interactions were identified by cross-correlation histograms. Cross-correlograms were compiled for 215 cell pairs and 84 (39%) showed significant features. The most frequently found feature (65/84 = 77%) was a central peak, straddling the origin and representing a source of common synaptic input to both cells. One third of these also had superimposed lagged peaks, indicative of a serial excitatory connection. Pure lagged peaks and lagged troughs, indicative of serial excitatory or inhibitory linkages, respectively, both occurred in 5% of the correlograms with features. A central trough appeared in 13% of the correlograms. The magnitude of the synaptic linkage was measured as the normalized area of the correlogram feature. Plotting the strength of synaptic interaction against response similarity during alternating wrist torques revealed a positive relationship for the correlated cell pairs. A linear fit yielded a positive slope: the pairs with excitatory interactions tended to covary more often than countervary. This linear fit had a positive offset, reflecting a tendency for both covarying and countervarying cells to have excitatory common input. Plotting the cortical location of the cell pairs showed that the strongest interactions occurred between cells separated by <400 microns. The correlational linkages between cells of different cortical layers showed a large proportion of common input to cells in layer V.

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.

Responses of human motoneurons to high-frequency stimulation of Ia afferents

This study was designed to extend to humans the findings of classical studies on anesthetized cats, which have examined the discharge of spinal motoneurons in response to high-frequency stimulus trains delivered to Ia afferents. Experiments were conducted on the monosynaptic pathway in the flexor carpi radialis (FCR) and soleus muscles. Subjects maintained a rhythmic discharge of a single motor unit (SMU) in either the FCR or soleus while homonymous Ia afferents were stimulated with either a single- or multipulse train. An n@IPI stimulus train had n pulses (n = 2-4) and an interpulse interval (IPI) of 1-8 ms. For each condition and motor unit, surface electromyographic (EMG) activity was averaged, and peristimulus-time histograms (PSTHs) were constructed for the SMU. The magnitude of the EMG was high for IPI = 1 ms, low for IPI = 2-3 ms, and high for IPI = 4-8 ms. SMU responses showed a similar pattern, which indicated that the increased EMG response was due to the presence of multiple peaks in a PSTH. The key results indicate that: (1) a short, high-frequency stimulus train enhances the discharge probability of a motoneuron above that observed with a single pulse; and (2) the increased motoneuron responses are significantly greater for the FCR than for the soleus muscle.

Functional relations between primate motor cortex cells and muscles: fixed and flexible

In behaving monkeys the effects of motor cortex cells on muscles are inferred from two quite different types of 'correlational' evidence: their coactivation and cross-correlation. Many precentral cells are coactivated with limb muscles, suggesting that they make a proportional contribution to muscle activity; however, such coactivation is typically quite flexible, and can be changed by operantly conditioning the dissociation of cell and muscle activity. Cross-correlating cells and muscles by spike-triggered averaging of the electromyogram (EMG) shows that certain cells produce short-latency post-spike facilitation of EMG; this correlational linkage is relatively fixed under different behavioural conditions and its time course suggests it is mediated by a corticomotoneuronal (CM) synaptic connection. CM cells typically facilitate a set of coactivated agonist muscles, and some also inhibit their antagonists. The firing patterns of CM cells can differ significantly from those of their target muscles. During ramp-and-hold wrist responses most CM cells discharge a phasic burst that precedes target muscle onset and that contributes to changes in muscle activity. At low force levels many CM cells are activated without their target motor units. Conversely, many CM cells are paradoxically inactive during rapid forceful movements that vigorously activate their target muscles; they appear to be preferentially active during finely controlled movements. Thus CM cells, with a fixed correlational linkage to their target muscles, may be recruited without their target muscles, and vice versa.

Spike timing and reliability in cortical pyramidal neurons: effects of EPSC kinetics, input synchronization and background noise on spike timing

In vivo studies have shown that neurons in the neocortex can generate action potentials at high temporal precision. The mechanisms controlling timing and reliability of action potential generation in neocortical neurons, however, are still poorly understood. Here we investigated the temporal precision and reliability of spike firing in cortical layer V pyramidal cells at near-threshold membrane potentials. Timing and reliability of spike responses were a function of EPSC kinetics, temporal jitter of population excitatory inputs, and of background synaptic noise. We used somatic current injection to mimic population synaptic input events and measured spike probability and spike time precision (STP), the latter defined as the time window (Δt) holding 80% of response spikes. EPSC rise and decay times were varied over the known physiological spectrum. At spike threshold level, EPSC decay time had a stronger influence on STP than rise time. Generally, STP was highest (≤2.45 ms) in response to synchronous compounds of EPSCs with fast rise and decay kinetics. Compounds with slow EPSC kinetics (decay time constants>6 ms) triggered spikes at lower temporal precision (≥6.58 ms). We found an overall linear relationship between STP and spike delay. The difference in STP between fast and slow compound EPSCs could be reduced by incrementing the amplitude of slow compound EPSCs. The introduction of a temporal jitter to compound EPSCs had a comparatively small effect on STP, with a tenfold increase in jitter resulting in only a five fold decrease in STP. In the presence of simulated synaptic background activity, precisely timed spikes could still be induced by fast EPSCs, but not by slow EPSCs.

A simple connectivity scheme for sparse coding in an olfactory system

Recent studies, using unbiased sampling of neuronal activity in vivo, indicate the existence of sparse codes in the brain. These codes are characterized by highly specific, associative (i.e., dependent on combinations of features) and often invariant neuronal responses. Sparse representations present many advantages for memory storage and are, thus, of wide interest in sensory physiology. Here, we study the statistics of connectivity in an olfactory network that contributes to the generation of such codes: Kenyon cells (KCs), the intrinsic neurons of the mushroom body (a structure involved in learning and memory in insects) receive inputs from a small population of broadly tuned principal neurons; from these inputs, KCs generate exquisitely selective responses and, thus, sparse representations. We find, surprisingly, that KCs are on average each connected to about 50% of their input population. Simple analysis indicates that such connectivity indeed maximizes the difference between input vectors to KCs and helps to explain their high specificity.

Neural spike train synchronization indices: definitions, interpretations, and applications

A comparison of previously defined spike train syncrhonization indices is undertaken within a stochastic point process framework. The second-order cumulant density (covariance density) is shown to be common to all the indices. Simulation studies were used to investigate the sampling variability of a single index based on the second-order cumulant. The simulations used a paired motoneurone model and a paired regular spiking cortical neurone model. The sampling variability of spike trains generated under identical conditions from the paired motoneurone model varied from 50% to 160% of the estimated value. On theoretical grounds, and on the basis of simulated data a rate dependence is present in all synchronization indices. The application of coherence and pooled coherence estimates to the issue of synchronization indices is considered. This alternative frequency domain approach allows an arbitrary number of spike train pairs to be evaluated for statistically significant differences, and combined into a single population measure. The pooled coherence framework allows pooled time domain measures to be derived, application of this to the simulated data is illustrated. Data from the cortical neurone model is generated over a wide range of firing rates (1-250 spikes/s). The pooled coherence framework correctly characterizes the sampling variability as not significant over this wide operating range. The broader applicability of this approach to multielectrode array data is briefly discussed.

Dynamic-clamp analysis of the effects of convergence on spike timing. I. Many synaptic inputs

Precise action potential timing is crucial in sensory acuity and motor control. Convergence of many synaptic inputs is thought to provide a means of decreasing spike-timing variability ("jitter"), but its effectiveness has never been tested in real neurons. We used the dynamic-clamp technique in mouse auditory brain stem slices to examine how convergence controls spike timing. We tested the roles of several synaptic properties that are influenced by ongoing activity in vivo: the number of active inputs (N), their total synaptic conductance (G(tot)), and their timing, which can resemble an alpha or a Gaussian distribution. Jitter was reduced most with large N, up to a factor of over 20. Variability in N is likely to occur in vivo, but this added little jitter. Jitter reduction also depended on the timing of inputs: alpha-distributed inputs were more effective than Gaussian-distributed inputs. Furthermore, the two distributions differed in their sensitivity to synaptic conductance: for Gaussian-distributed inputs, jitter was most reduced when G(tot) was 2-3 times threshold, whereas alpha-distributed inputs showed continued jitter reduction with higher G(tot). However, very high G(tot) caused the postsynaptic cell to fire multiple times, particularly when the input jitter outlasted the cell's refractory period, which interfered with jitter reduction. G(tot) also greatly affected the response latency, which could influence downstream computations. Changes in G(tot) are likely to arise in vivo through activity-dependent changes in synaptic strength. High rates of postsynaptic activity increased the number of synaptic inputs required to evoke a postsynaptic response. Despite this, jitter was still effectively reduced, particularly in cases when this increased threshold eliminated secondary spikes. Thus these studies provide insight into how specific features of converging inputs control spike timing.

Factors determining the precision of the correlated firing generated by a monosynaptic connection in the cat visual pathway

Across the visual pathway, strong monosynaptic connections generate a precise correlated firing between presynaptic and postsynaptic neurons. The precision of this correlated firing is not the same within thalamus and visual cortex. While retinogeniculate connections generate a very narrow peak in the correlogram (peak width < 1 ms), the peaks generated by geniculocortical and corticocortical connections have usually a time course of several milliseconds. Several factors could explain these differences in timing precision such as the amplitude of the monosynaptic EPSP (excitatory postsynaptic potential), its time course or the contribution of polysynaptic inputs. While it is difficult to isolate the contribution of each factor in physiological experiments, a first approximation can be done in modelling studies. Here, we simulated two monosynaptically connected neurons to measure changes in their correlated firing as we independently modified different parameters of the connection. Our results suggest that the precision of the correlated firing generated by strong monosynaptic connections is mostly determined by the EPSP time course of the connection and much less by other factors. In addition, we show that a polysynaptic pathway is unlikely to emulate the correlated firing generated by a monosynaptic connection unless it generates EPSPs with very small latency jitter.

Background synaptic conductance and precision of EPSP-spike coupling at pyramidal cells

The temporal precision of converting excitatory postsynaptic potentials (EPSPs) into spikes at pyramidal cells is critical for the coding of information in the cortex. Several in vitro studies have shown that voltage-dependent conductances in pyramidal cells can prolong the EPSP time course resulting in an imprecise EPSP-spike coupling. We have used dynamic-clamp techniques to mimic the in vivo background synaptic conductance in cortical slices and investigated how the ongoing synaptic activity may affect the EPSP time course near threshold and the EPSP spike coupling. We report here that background synaptic conductance dramatically diminished the depolarization related prolongation of the EPSPs in pyramidal cells and improved the precision of spike timing. Furthermore, we found that background synaptic conductance can affect the interaction among action potentials in a spike train. Thus the level of ongoing synaptic activity in the cortex may regulate the capacity of pyramidal cells to process temporal information.

Estimation of postsynaptic potentials in rat hypoglossal motoneurones: insights for human work

Classical techniques for estimating postsynaptic potentials in motoneurones include spike-triggered averages of rectified surface and multiunit electromyographic recordings (SEMG and MU-EMG), as well as the compilation of peristimulus time histograms (PSTH) based on the discharge of single motor units (SMU). These techniques rely on the probability of spike occurrence in relation to the stimulus and can be contaminated by count- and synchronization-related errors, arising from post-spike refractoriness and the discharge statistics of motoneurones. On the other hand, since these probability-based techniques are easy to use and require only inexpensive equipment, it is very likely that they will continue to be used in clinical and laboratory settings for the foreseeable future. One aim of the present study was to develop a modification of these probability-based analyses in order to provide a better estimate of the initial phase of postsynaptic potentials. An additional aim was to combine probability-based analyses with frequency-based analyses to provide a more reliable estimate of later phases of postsynaptic potentials. To achieve these aims, we have injected simple as well as complex current transients into regularly discharging hypoglossal motoneurones recorded in vitro from rat brainstem slices. We examined the discharge output of these cells using both probability- and frequency-based analyses to identify which of the two represented the profile of the postsynaptic potential more closely. This protocol was designed to obtain PSTHs of the responses of single motor units to repeated application of the same afferent input. We have also simulated multiunit responses to afferent input by replacing the times of spike occurrence in individual trials with a representation of either an intramuscular or surface-recording single motor unit waveform and summing many of these trials to obtain either a simulated SEMG or MU-EMG. We found that in a regularly discharging motoneurone, the rising phase of an EPSP moves the occurrence of spikes forward and hence induces a substantial peak in all probability-based records. This peak is followed immediately by a period of reduced activity ('silent period') due to the phase advancement of spikes that were to occur at this period. Similarly, the falling phase of an IPSP delays spikes so that they occur during the rising phase of the IPSP. During the delay, the probability-based analyses display gaps and during the occurrence of the delayed spikes they generate peaks. We found that all the probability-based analyses (SEMG, MU-EMG and PSTH) can be made useful for illustrating the underlying initial PSP by a special use of the cumulative sum (CUSUM) calculation. We have illustrated that, in most cases, the CUSUM of probability-based analyses can overcome the delay- or advance-related (i.e. the count-related) errors of the classical methods associated with the first PSP only. The probability-based records also induce secondary and tertiary peaks and troughs due to synchronization of the spikes in relation to the stimulus (i.e. the synchronization-related errors) by the first PSP to occur at fixed times from the stimulus. Special CUSUM analyses cannot overcome these synchronization-related errors. Frequency-based analysis (PSFreq) of individual and summed trials gave comparable and often better indications of the underlying PSPs than the probability-based analyses. When used in combination, these analyses compliment each other so that a more accurate estimation of the underlying PSP is possible. Since the correct identification of the connections in the central nervous system is of utmost importance in order to understand the operation of the system, we suggest that as well as the using the special CUSUM approach on probability-based records, researchers should seriously consider the use of frequency-based analyses in their indirect estimation of stimulus-induced compound synaptic potentials in human motoneurones.

Thalamocortical control of feed-forward inhibition in awake somatosensory 'barrel' cortex

Intracortical inhibition plays a role in shaping sensory cortical receptive fields and is mediated by both feed-forward and feedback mechanisms. Feed-forward inhibition is the faster of the two processes, being generated by inhibitory interneurons driven by monosynaptic thalamocortical (TC) input. In principle, feed-forward inhibition can prevent targeted cortical neurons from ever reaching threshold when TC input is weak. To do so, however, inhibitory interneurons must respond to TC input at low thresholds and generate spikes very quickly. A powerful feed-forward inhibition would sharpen the tuning characteristics of targeted cortical neurons, and interneurons with sensitive and broadly tuned receptive fields could mediate this process. Suspected inhibitory interneurons (SINs) with precisely these properties are found in layer 4 of the somatosensory (S1) 'barrel' cortex of rodents and rabbits. These interneurons lack the directional selectivity seen in most cortical spiny neurons and in ventrobasal TC afferents, but are much more sensitive than cortical spiny neurons to low-amplitude whisker displacements. This paper is concerned with the activation of S1 SINs by TC impulses, and with the consequences of this activation. Multiple TC neurons and multiple S1 SINs were simultaneously studied in awake rabbits, and cross-correlation methods were used to examine functional connectivity. The results demonstrate a potent, temporally precise, dynamic and highly convergent/divergent functional input from ventrobasal TC neurons to SINs of the topographically aligned S1 barrel. Whereas the extensive pooling of convergent TC inputs onto SINs generates sensitive and broadly tuned inhibitory receptive fields, the potent TC divergence onto many SINs generates sharply synchronous activity among these elements. This TC feed-forward inhibitory network is well suited to provide a fast, potent, sensitive and broadly tuned inhibition of targeted spiny neurons that will suppress spike generation following all but the most optimal feed-forward excitatory inputs.

Firing rate modulation of motoneurons activated by cutaneous and muscle receptor afferents in the decerebrate cat

The aim of this study was to investigate whether activation of spinal motoneurons by sensory afferents of the caudal cutaneous sural (CCS) nerve evokes an atypical motor control scheme. In this scheme, motor units that contract fast and forcefully are driven by CCS afferents to fire faster than motor units that contract more slowly and weakly. This is the opposite of the scheme described by the size principle. Earlier studies from this lab do not support the atypical scheme and instead demonstrate that both CCS and muscle stretch recruit motor units according to the size principle. The latter finding may indicate that CCS and muscle-stretch inputs have similar functional organizations or that comparison of recruitment sequence was simply unable to resolve a difference. In the present experiments, we examine this issue using rate modulation as a more sensitive index of motoneuron activation than recruitment. Quantification of the firing output generated by these two inputs in the same pairs of motoneurons enabled direct comparison of the functional arrangements of CCS versus muscle-stretch inputs across the pool of medial gastrocnemius (MG) motoneurons. No systematic difference was observed in the rate modulation produced by CCS versus muscle-stretch inputs for 35 pairs of MG motoneurons. For the subset of 24 motoneuron pairs exhibiting linear co-modulation of firing rate (r > 0.5) in response to both CCS and muscle inputs, the slopes of the regression lines were statistically indistinguishable between the two inputs. For individual motoneuron pairs, small differences in slope between inputs were not related to differences in conduction velocity (CV), recruitment order, or, for a small sample, differences in motor unit force. We conclude that an atypical motor control scheme involving selective activation of typically less excitable motoneurons, if it does occur during normal movement, is not an obligatory consequence of activation by sural nerve afferents. On average and for both muscle-stretch and skin-pinch inputs, the motoneuron with the faster CV in the pair tended to be driven to fire at slightly but significantly faster firing rates. Computer simulations based in part on frequency-current relations measured directly from motoneurons revealed that properties intrinsic to motoneurons are sufficient to account for the higher firing rates of the faster CV motoneuron in a pair.

Input-output functions of mammalian motoneurons

Our intent in this review was to consider the relationship between the biophysical properties of motoneurons and the mechanisms by which they transduce the synaptic inputs they receive into changes in their firing rates. Our emphasis has been on experimental results obtained over the past twenty years, which have shown that motoneurons are just as complex and interesting as other central neurons. This work has shown that motoneurons are endowed with a rich complement of active dendritic conductances, and flexible control of both somatic and dendritic channels by endogenous neuromodulators. Although this new information requires some revision of the simple view of motoneuron input-output properties that was prevalent in the early 1980's (see sections 2.3 and 2.10), the basic aspects of synaptic transduction by motoneurons can still be captured by a relatively simple input-output model (see section 2.3, equations 1-3). It remains valid to describe motoneuron recruitment as a product of the total synaptic current delivered to the soma, the effective input resistance of the motoneuron and the somatic voltage threshold for spike initiation (equations 1 and 2). However, because of the presence of active channels activated in the subthreshold range, both the delivery of synaptic current and the effective input resistance depend upon membrane potential. In addition, activation of metabotropic receptors by achetylcholine, glutamate, noradrenaline, serotonin, substance P and thyrotropin releasing factor (TRH) can alter the properties of various voltage- and calcium-sensitive channels and thereby affect synaptic current delivery and input resistance. Once motoneurons are activated, their steady-state rate of repetitive discharge is linearly related to the amount of injected or synaptic current reaching the soma (equation 3). However, the slope of this relation, the minimum discharge rate and the threshold current for repetitive discharge are all subject to neuromodulatory control. There are still a number of unresolved issues concerning the control of motoneuron discharge by synaptic inputs. Under dynamic conditions, when synaptic input is rapidly changing, time- and activity-dependent changes in the state of ionic channels will alter both synaptic current delivery to the spike-generating conductances and the relation between synaptic current and discharge rate. There is at present no general quantitative expression for motoneuron input-output properties under dynamic conditions. Even under steady-state conditions, the biophysical mechanisms underlying the transfer of synaptic current from the dendrites to the soma are not well understood, due to the paucity of direct recordings from motoneuron dendrites. It seems likely that resolving these important issues will keep motoneuron afficiandoes well occupied during the next twenty years.

Spike transmission and synchrony detection in networks of GABAergic interneurons

The temporal pattern and relative timing of action potentials among neocortical neurons may carry important information. However, how cortical circuits detect or generate coherent activity remains unclear. Using paired recordings in rat neocortical slices, we found that the firing of fast-spiking cells can reflect the spiking pattern of single-axon pyramidal inputs. Moreover, this property allowed groups of fast-spiking cells interconnected by electrical and gamma-aminobutyric acid (GABA)-releasing (GABAergic) synapses to detect the relative timing of their excitatory inputs. These results indicate that networks of fast-spiking cells may play a role in the detection and promotion of synchronous activity within the neocortex.

Precise burst synchrony in the superior colliculus of the awake cat during moving stimulus presentation

This study aimed to characterize the synchrony that occurs between cell discharges in the superior colliculus of the awake cat. We trained cats to perform a visual fixation in the presence of a visual moving stimulus and then recorded 686 pairs of neighboring cells in the superior colliculus during task performance. A new method to assess the significance of precise discharge synchronization is described, which permits analysis of nonstationary data. Of 181 pairs with sufficient data for quantitative analysis, 125 showed a cross-correlation histogram (CCH) with features assessed as significant using this approach. CCHs frequently showed an isolated central peak (41 of 125) or a peak flanked by one or two troughs (68 of 125), and in a few cases an oscillatory pattern of approximately 65 Hz (16 of 125). This is in contrast to the oscillation frequency reported for the visual cortex and shows that oscillations in the superior colliculus probably arise from a cortex-independent mechanism. Our method also permits direct quantification of the correlation shift predictors, assessing precise time locking of spikes to the stimulus. Only 1 of 125 cross-correlation shift predictors had a significant central peak, meaning that most of the CCH features were not related to cell discharges time-locked to the stimulus presentation. Further investigation using a burst-jittering method showed that synchrony in the superior colliculus is attributable to precise synchronization of short bursts of spikes. Such synchrony could be related to the network dynamics and the common inhibitory feedback from local interneurons, which would act as temporal selectors of the cells with greatest or fastest response.

Precise burst synchrony in the superior colliculus of the awake cat during moving stimulus presentation

This study aimed to characterize the synchrony that occurs between cell discharges in the superior colliculus of the awake cat. We trained cats to perform a visual fixation in the presence of a visual moving stimulus and then recorded 686 pairs of neighboring cells in the superior colliculus during task performance. A new method to assess the significance of precise discharge synchronization is described, which permits analysis of nonstationary data. Of 181 pairs with sufficient data for quantitative analysis, 125 showed a cross-correlation histogram (CCH) with features assessed as significant using this approach. CCHs frequently showed an isolated central peak (41 of 125) or a peak flanked by one or two troughs (68 of 125), and in a few cases an oscillatory pattern of approximately 65 Hz (16 of 125). This is in contrast to the oscillation frequency reported for the visual cortex and shows that oscillations in the superior colliculus probably arise from a cortex-independent mechanism. Our method also permits direct quantification of the correlation shift predictors, assessing precise time locking of spikes to the stimulus. Only 1 of 125 cross-correlation shift predictors had a significant central peak, meaning that most of the CCH features were not related to cell discharges time-locked to the stimulus presentation. Further investigation using a burst-jittering method showed that synchrony in the superior colliculus is attributable to precise synchronization of short bursts of spikes. Such synchrony could be related to the network dynamics and the common inhibitory feedback from local interneurons, which would act as temporal selectors of the cells with greatest or fastest response.

Amplification and linear summation of synaptic effects on motoneuron firing rate

The aim of this study was to measure the effects of synaptic input on motoneuron firing rate in an unanesthetized cat preparation, where activation of voltage-sensitive dendritic conductances may influence synaptic integration and repetitive firing. In anesthetized cats, the change in firing rate produced by a steady synaptic input is approximately equal to the product of the effective synaptic current measured at the resting potential (I(N)) and the slope of the linear relation between somatically injected current and motoneuron discharge rate (f-I slope). However, previous studies in the unanesthetized decerebrate cat indicate that firing rate modulation may be strongly influenced by voltage-dependent dendritic conductances. To quantify the effects of these conductances on motoneuron firing behavior, we injected suprathreshold current steps into medial gastrocnemius motoneurons of decerebrate cats and measured the changes in firing rate produced by superimposed excitatory synaptic input. In the same cells, we measured I(N) and the f-I slope to determine the predicted change in firing rate (Delta F = I(N) * f-I slope). In contrast to previous results in anesthetized cats, synaptically induced changes in motoneuron firing rate were greater-than-predicted. This enhanced effect indicates that additional inward current was present during repetitive firing. This additional inward current amplified the effective synaptic currents produced by two different excitatory sources, group Ia muscle spindle afferents and caudal cutaneous sural nerve afferents. There was a trend toward more prevalent amplification of the Ia input (14/16 cells) than the sural input (11/16 cells). However, in those cells where both inputs were amplified (10/16 cells), amplification was similar in magnitude for each source. When these two synaptic inputs were simultaneously activated, their combined effect was generally very close to the linear sum of their amplified individual effects. Linear summation is also observed in medial gastrocnemius motoneurons of anesthetized cats, where amplification is not present. This similarity suggests that amplification does not disturb the processes of synaptic integration. Linear summation of amplified input was evident for the two segmental inputs studied here. If these phenomena also hold for other synaptic sources, then the presence of active dendritic conductances underlying amplification might enable motoneurons to integrate multiple synaptic inputs and drive motoneuron firing rates throughout the entire physiological range in a relatively simple fashion.

The influence of single VB thalamocortical impulses on barrel columns of rabbit somatosensory cortex

Extracellular recordings were obtained from single neurons in ventrobasal (VB) thalamus of awake rabbits while field potentials were recorded at various depths within topographically aligned and nonaligned barrel columns of somatosensory cortex (S1). Spike-triggered averages of cortical field potentials were obtained following action potentials in thalamic neurons. Action potentials in a VB neuron elicited a cortical response within layer 4 with three distinct components. 1) A biphasic, initially positive response (latency <1 ms) was interpreted to reflect activation of the VB axon terminals (the AxTP). This response was not affected by infusion of an alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA)/kainate receptor antagonist within the barrel. In contrast, later components of the response were completely eliminated and were interpreted to reflect focal synaptic potentials. 2) A negative potential [focal synaptic negativity (FSN)] occurred at a mean latency of 1.65 ms and lasted approximately 4 ms. This response had a rapid rise time ( approximately 0.7 ms) and was interpreted to reflect monosynaptic excitation. 3) The third component was a positive potential (the FSP), with a slow rise time and a half-amplitude duration of approximately 30 ms. The FSP showed a weak reversal in superficial cortical layers and was interpreted to reflect di/polysynaptic inhibition. The amplitudes of the AxTP, the FSN, and the FSP reached a peak near layer 4 and were highly attenuated in both superficial and deep cortical layers. All components were attenuated or absent when the cortical electrode was missaligned from the thalamic electrode by a single cortical barrel. Deconvolution procedures revealed that the autocorrelogram of the presynaptic VB neuron had very little influence on either the amplitude or duration of the AxTP or the FSN, and only a minor influence (mean, 11%) on the amplitude of the FSP. We conclude that individual VB thalamic impulses entering a cortical barrel engage both monosynaptic excitatory and di/polysynaptic inhibitory mechanisms. Putative inhibitory interneurons of an S1 barrel receive a highly divergent/convergent monosynaptic input from the topographically aligned VB barreloid, and this results in sharp synchrony among these interneurons. We suggest that single-fiber access to disynaptic inhibition is facilitated by this sharp synchrony, and that the FSP reflects a consequent synchronous wave of feed-forward inhibition within the S1 barrel.

Correlations between neurograms and locomotor drive potentials in motoneurons during fictive locomotion: implications for the organization of locomotor commands

The patterns of correlation found between motoneuron pools during fictive locomotion are the same whether the coherence functions used to detect the correlations are determined using pairs of rectified ENGs or motoneuron LDPs and rectified ENGs. This finding suggests that the higher frequencies in rectified ENGs (and, perhaps, EMGs) contain information about the synaptic input to motoneurons. Nevertheless, differences between the coherence functions of rectified ENG pairs and those of LDPs and rectified ENGs suggests that this information is distorted by harmonics introduced by rectification. The activities of many motoneuron pools are correlated during the flexor or extensor phase of fictive locomotion, indicating that they receive common synaptic input from branched presynaptic axons or from pools of interneurons whose activities are synchronized. Similar findings were reported by Bayev (1978), based on temporal correlations. Our results indicate that the investigated motor nuclei, which innervate muscles with actions at the hip, knee and ankle, are subject to a set of common locomotor commands. These commands are also received by inhibitory interneurons that project to the motor nuclei of antagonists, as indicated by the correlations between the hyperpolarizing phase of LDPs and activity in the rectified ENGs of antagonists. This last result is consistent with a modular organization for the spinal locomotor generator, in which one set of interneurons drives a motor pool and the inhibitory interneurons that project to the motor pool's antagonist (Jordan, 1991). However, these results also suggest that the spinal modules for locomotion may not be separable into independent unit-burst generators that produce commands for control of each joint as Grillner (1981) has suggested. Our results are more consistent with a model in which a generator distributes flexor and extensor commands to many motor pools (like the half-center model) with as yet unidentified spinal mechanisms that determine differences in the initiation and termination of activity of individual motor nuclei. Alternatively, the correlations between motor pools that we have observed could be explained by spinal mechanisms that synchronize the activity of unit-burst type generators. Despite the distribution of common locomotor commands to many functionally diverse motor nuclei, the spinal locomotor pattern generator is differentiated to the extent that some motor nuclei, like EDL and FDL, receive separate locomotor commands. This conclusion is consistent with other observations. EDL and FDL display distinctive, individualized patterns of locomotor activity that may vary in a facultative manner or in different forms of locomotion (O'Donovan et al., 1980; Trank et al., 1996). A recent study has shown that during fictive locomotion EDL and FDL motoneurons receive input from different sets of last-order interneurons than those which project to other motor pools (Degtyarenko et al., 1998). These results suggest that spinal locomotor generators are differentiated for the individualized control of some digit muscles, like FDL and EDL.

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.

Electromechanical coupling and synchronous firing of single wrist extensor motor units in sporadic amyotrophic lateral sclerosis

Electrical and contractile properties of motor units (MU) were studied in the extensor carpi radialis muscles during voluntary contraction. The discharge of 234 single MUs was recorded in 11 patients with sporadic amyotrophic lateral sclerosis (ALS) and compared with that of the 260 MUs recorded in 12 healthy control subjects. Characteristics of the MU twitches and of the macro-potentials, the electromechanical coupling and the synchronization of the motor neurone discharges, were compared. In 5 patients (population ALS1), the twitch contraction force and macro-MUP area values were much larger than those of the controls. In the 6 other patients (population ALS2), the twitch force was considerably depressed, whereas the macro-MUP area was slightly, but significantly, increased. In ALS1, as well as in ALS2, the electromechanical coupling was much weaker than in the controls, and the fast-contracting MUs were more severely affected than the slowly contracting MUs. The motoneuronal synchronization was assessed by performing cross-correlation analysis on MUs discharges, and was used as an index to the strength of the common motoneuronal inputs. The rate of occurrence of synchronous firing was conspicuously lower in both populations of patients than in the control group. This might reflect the loss of corticospinal projections that occurs in ALS. The data are discussed in terms of the time course of motor neurone axonal sprouting, and in terms of the neuronal and muscular dysfunction possibly involved in ALS disease.

Isodirectional tuning of adjacent interneurons and pyramidal cells during working memory: evidence for microcolumnar organization in PFC

Studies on the cellular mechanisms of working memory demonstrated that neurons in dorsolateral prefrontal cortex (dPFC) exhibit directionally tuned activity during an oculomotor delayed response. To determine the particular contributions of pyramidal cells and interneurons to spatial tuning in dPFC, we examined both individually and in pairs the tuning properties of regular-spiking (RS) and fast-spiking (FS) units that represent putative pyramidal cells and interneurons, respectively. Our main finding is that FS units possess spatially tuned sensory, motor, and delay activity (i. e., "memory fields") similar to those found in RS units. Furthermore, when recorded simultaneously at the same site, the majority of neighboring neurons, whether FS or RS, displayed isodirectional tuning, i.e., they shared very similar tuning angles for the sensory and delay phases of the task. As the trial entered the response phase of the task, many FS units shifted their direction of tuning and became cross-directional to adjacent RS units by the end of the trial. These results establish that a large part of inhibition in prefrontal cortex is spatially oriented rather than being untuned and simply regulating the threshold response of pyramidal cell output. Moreover, the isodirectional tuning between adjacent neurons supports a functional microcolumnar organization in dPFC for spatial memory fields similar to that found in other areas of cortex for sensory receptive fields.

A comparison of human motoneuron data to simulated data using cat motoneuron models

The response of repetitively firing human motoneurons to a composite excitatory input was evaluated. It was clearly shown that the response of the motoneurons to the transient input decreased with an increase in the background firing rate of the cell. The current model of repetitively firing human motoneurons could not account for this experimental result. Therefore, a compartmental modelling approach was used to simulate the repetitive firing properties of anaesthetised cat motoneurons under current clamp conditions. The modelled motoneurons were used in simulations similar to the experimental paradigms where the response to a composite excitatory input was evaluated at different background firing rates. The motoneuron models also showed a decrease in response to the excitatory input at faster background firing rates. The results suggest that human motoneurons are more comparable to motoneurons in the anaesthetised cat preparation than formerly thought. The results also demonstrate that the apparent efficacy of a synaptic input may be modulated by changes in background firing rate of the postsynaptic neuron.

Functional identification of the input-output transforms of mammalian motoneurones

We studied the responses of rat hypoglossal and cat lumbar motoneurones to a variety of excitatory and inhibitory injected current transients during repetitive discharge. The amplitudes and time courses of the transients were comparable to those of the synaptic currents underlying postsynaptic potentials (PSPs) recorded in these cells. Poisson trains of these current transients were combined with an additional independent, high frequency random waveform to approximate band-limited white noise. The composite, white noise waveform was then superimposed on long duration suprathreshold current steps. We used the responses of the motoneurones to the white noise stimulus to derive zero-, first- and second-order Wiener kernels, which provide a quantitative description of the relation between injected current and discharge probability. The convolution integral computed for an injected current waveform and the first-order Wiener kernel provides the best linear prediction of the associated peristimulus time histogram (PSTH). This linear model provided good matches to most of the PSTHs compiled between the times of occurrence of individual current transients and motoneurone discharges. However, for the largest amplitude current transients, a significant improvement in the PSTH match was often achieved by expanding the model to include the convolution of the second-order Wiener kernel with the input. The overall transformation of current inputs into firing rate could be approximated by a second-order Wiener Model, i.e., a cascade of a dynamic, linear filter followed by a static non-linearity. At a given mean firing rate, the non-linear component of the motoneurone's response could be described by the square of the linear component multiplied by a constant coefficient. The amplitude of the response of the linear component increased with the average firing rate, whereas the value of the multiplicative coefficient in the nonlinear component decreased. As a result, the overall transform could be predicted from the mean firing rate and the linear impulse response, yielding a relatively simple, general description of the motoneurone's input-output function.

Activity of spinal interneurons and their effects on forearm muscles during voluntary wrist movements in the monkey

We studied the activity of 577 neurons in the C6-T1 spinal cord of three awake macaque monkeys while they generated visually guided, isometric flexion/extension torques about the wrist. Spike-triggered averaging of electromyographic activity (EMG) identified the units' correlational linkages with </=12 forearm muscles. One hundred interneurons produced changes in the level of average postspike EMG with onset latencies consistent with mono- or oligosynaptic connections to motoneurons; these were classified as premotor interneurons (PreM-INs). Most PreM-INs (82%) produced postspike facilitations in forearm muscles. Earlier spike-related features, often beginning before the trigger spike, were seen in spike-triggered averages from 72 neurons. Postspike effects were present in one muscle for 64% of the PreM-INs. Neurons with divergent linkages to larger "muscle fields" usually generated postspike effects in synergistic muscles. Fifty-eight percent of the PreM-INs had postspike effects in flexor muscles only and 29% in extensor muscles only. Postspike effects were distributed relatively evenly among the primary flexor and extensor muscles studied. The mean percent change in EMG level from baseline and the mean onset latencies for postspike facilitations and postspike suppressions were similar. PreM-INs exhibited a variety of response patterns during the generation of isometric wrist torque. The response patterns and output effects of 24% of the PreM-INs were consistent with a strict reciprocal organization of flexor and extensor muscle control. For another 60% of the PreM-INs, there was a congruent relation between activity and output effects for only one direction of torque production. These neurons were active for both flexion and extension torques, including 37 neurons that exhibited bidirectional increases in discharge rate. The relatively small number of postspike suppressions observed suggests that inhibitory interneurons were silent when their target muscles were recruited. Compared with premotor neurons in the motor cortex, the red nucleus and the C8-T1 dorsal root ganglia, spinal PreM-INs affected flexor muscles in greater proportions and had smaller muscle fields. The magnitudes of postspike facilitations were similar in all premotor populations. Bidirectional activity, common for PreM-INs, was rare for corticomotoneuronal and premotor dorsal root ganglion cells, which discharge only for torques in their preferred direction.

Frequency-dependent synaptic depression modifies postsynaptic firing probability in cats

1. The influence of stimulus trains applied to single I a axons on the firing behaviour of single motoneurones was assessed in anaesthetized cats. The change in motoneurone firing probability associated with a single I a afferent spike was measured from short-latency peaks in peristimulus time histograms or cross-correlograms. Some synapses showed frequency-dependent depression of the short-latency peak, which is consonant with the frequency-dependent depression reported for the I a-motoneurone excitatory postsynaptic potential (EPSP). 2. Where they could be measured, EPSPs superimposed on the depolarizing ramps of potential recorded from motoneurones as they fired repetitively showed frequency-dependent changes in amplitude that parallelled those of the simultaneously recorded histograms. 3. Thus it appears that at synapses with small EPSPs, which are typical in the mammalian CNS, modulation of the EPSP should result in similar modulation of cell firing.

Computer simulation of post-spike facilitation in spike-triggered averages of rectified EMG

When the spikes of a motor cortical cell are used to compile a spike-triggered average (STA) of rectified electromyographic (EMG) activity, a post-spike facilitation (PSF) is sometimes seen. This is generally thought to be indicative of direct corticomotoneuronal (CM) connections. However, it has been claimed that a PSF could be caused by synchronization between CM and non-CM cells. This study investigates the generation of PSF using a computer model. A population of cortical cells was simulated, some of which made CM connections to a pool of 103 motoneurons. Motoneurons were simulated using a biophysically realistic model. A subpopulation of the cortical cells was synchronized together. After a motoneuron discharge, a motor unit action potential was generated; these were summed to produce an EMG output. Realistic values were used for the corticospinal and peripheral nerve conduction velocity distribution, for slowing of impulse conduction in CM terminal axons, and for the amount of cortical synchrony. STA of the rectified EMG from all cortical neurons showed PSF; however, these were qualitatively different for CM versus non-CM cells. Using an epoch analysis to determine reliability in a quantitative manner, it was shown that the onset latency of PSF did not distinguish the two classes of cells after 10,000 spikes because of high noise in the averages. The time of the PSF peak and the peak width at half-maximum (PWHM) could separate CM from synchrony effects. However, only PWHM was robust against changes in motor unit action-potential shape and duration and against changes in the width of cortical synchrony. The amplitude of PSF from a CM cell could be doubled by the presence of synchrony. It is proposed that, if a PSF has PWHM < 7 ms, this reliably indicates that the trigger is a CM cell projecting to the muscle whose EMG is averaged. In an analysis of experimental data where macaque motor cortical cells facilitated hand and forearm muscle EMG, 74% of PSFs fulfilled this criterion. The PWHM criterion could be applied to other STA studies in which it is important to exclude the effects of synchrony.

Identification of patterns of neuronal connectivity--partial spectra, partial coherence, and neuronal interactions

The cross-correlation histogram has provided the primary tool for inferring the structure of common inputs to pairs of neurones. While this technique has produced useful results it not clear how it may be extended to complex networks. In this report we introduce a linear model for point process systems. The finite Fourier transform of this model leads to a regression type analysis of the relations between spike trains. An advantage of this approach is that the full range of techniques for multivariate regression analyses becomes available for spike train analysis. The two main parameters used for the identification of neural networks are the coherence and partial coherences. The coherence defines a bounded measure of association between two spike trains and plays the role of a squared correlation coefficient defined at each frequency lambda. The partial coherences, analogous to the partial correlations of multiple regression analysis, allow an assessment of how any number of putative input processes may influence the relation between any two output processes. In many cases analytic solutions may be found for coherences and partial coherences for simple neural networks, and in combination with simulations may be used to test hypotheses concerning proposed networks inferred from spike train analyses.

Functional identification of the input-output transforms of motoneurones in the rat and cat

1. We studied the responses of rat hypoglossal and cat lumbar motoneurones to a variety of excitatory and inhibitory injected current transients during repetitive discharge. The amplitudes and time courses of the transients were comparable to those of the synaptic currents underlying unitary and small compound postsynaptic potentials (PSPs) recorded in these cells. Poisson trains of ten of these excitatory and ten inhibitory current transients were combined with an additional independent, high-frequency random waveform to approximate band limited white noise. The white noise waveform was then superimposed on long duration (39 s) suprathreshold current steps. 2. We measured the effects of each of the current transients on motoneurone discharge by compiling peristimulus time histograms (PSTHs) between the times of occurrence of individual current transients and motoneurone discharges. We estimated the changes in membrane potential associated with each current transient by approximating the passive response of the motoneurone with a simple resistance-capacitance circuit. The relations between the features of these simulated PSPs and those of the PSTHs were similar to those reported previously for real PSPs: the short-latency PSTH peak (or trough) was generally longer than the initial phase of the PSP derivative, but shorter than the time course of the PSP itself. Linear models of the PSP to PSTH transform based on the PSP time course, the time derivative of the PSP, or a linear combination of the two parameters could not reproduce the full range of PSTH profiles observed. 3. We also used the responses of the motoneurones to the white noise stimulus to derive zero-, first- and second-order Wiener kernels, which provide a quantitative description of the relation between injected current and discharge probability. The convolution integral computed for an injected current waveform and the first-order Wiener kernel should provide the best linear prediction of the associated PSTH. This linear model provided good matches to the PSTHs associated with a wide range of current transients. However, for the largest amplitude current transients, a significant improvement in the PSTH match was often achieved by expanding the model to include the convolution of the second-order Wiener kernel with the input. 4. The overall transformation of current inputs into firing rate could be approximated by a second-order Wiener model, i.e. a cascade of a dynamic, linear filter followed by a static non-linearity. At a given mean firing rate, the non-linear component of the response of the motoneurone could be described by the square of the linear component multiplied by a constant coefficient. The amplitude of the response of the linear component increased with the average firing rate, whereas the value of the multiplicative coefficient in the non-linear component decreased. As a result, the overall transform could be predicted from the mean firing rate and the linear impulse response, yielding a relatively simple, general description of the motoneurone input-output function.

Evidence from motoneurone synchronization for disynaptic pathways in the control of inspiratory motoneurones in the cat

1. Motoneurone synchronization was measured by cross-correlation between paired inspiratory discharges in external and internal intercostal nerves or their intramuscular branches (T3 to T8) or in the phrenic nerve (C5 root or both C5 and C6 roots independently) in anaesthetized, paralysed cats. 2. All cross-correlation histograms showed central peaks, for which the durations at half-amplitude (half-widths) from internal nerve pairs in adjacent segments were all less than for external nerve pairs in adjacent segments or within a segment (means, 1.6 ms vs. 3.4 ms for adjacent segments). Values for external-internal pairs covered the ranges for both these two. Lowest values came from two phrenic pairs (1.2 and 1.4 ms). 3. The peaks from ipsisegmental external-internal pairs were usually asymmetric and the maximum of the peak was often displaced to a lag of about -1 ms (external nerve providing the reference spikes), whereas peaks from external-external pairs were always symmetrical and centred on zero. Phrenic-internal peaks gave maxima with lags about 1 ms less than for phrenic-external peaks from the same segments. 4. Two explanations were considered possible for the differences in duration and timing: an extra synapse on the pathway to the external nerve motoneurones, or a correlation kernel for a monosynaptic connection to the external nerve motoneurones that had a slower time course than that for the internal or phrenic nerve motoneurones. Computer simulations, assuming the extra synapse, gave a good fit to the observed time courses of the correlation peaks for all categories of nerve pairs using single values of parameters (e.g. EPSP rise time) consistent with those in the literature. This could not be achieved with the different correlation kernel model. The timing of high-frequency oscillation (HFO), which was sometimes present in the correlations, was also better predicted with the extra synapse model. 5. It is concluded that most of the synchronization between external nerve motoneurones is derived from disynaptic common inputs and that any motoneurone synchronization peak with a half-width greater than about 2.2 ms should be assumed to be likely to contain di- or oligosynaptically derived components.

Synaptic interactions between primate precentral cortex neurons revealed by spike-triggered averaging of intracellular membrane potentials in vivo

To document synaptic interactions between neurons in the precentral cortex of macaque monkeys, we recorded in vivo the intracellular (IC) membrane potentials of cortical neurons simultaneously with extracellular (EC) action potentials of neighboring cells. The synaptic potentials correlated with EC spikes were obtained by spike-triggered averages (STA) of the IC membrane potentials for 373 cell pairs recorded in anesthetized and awake behaving monkeys. Sixty-three STAs (17%) showed excitatory postsynaptic potentials (EPSPs), beginning after the trigger spike. Pure EPSPs had onset latencies of 0.9 +/- 0.7 msec (mean +/- SD) and amplitudes of 226 +/- 130 microV. Sixteen STAs (4%) showed postspike inhibitory postsynaptic potentials (IPSPs), with onset latencies of 0.4 +/- 0.4 msec and amplitudes of -274 +/- 188 microV. The most common waveform, observed in 82% of the STAs with features, was a broad depolarization straddling the trigger spikes, reflecting synchronized synaptic input to both IC and EC neurons. These average synchronous excitation potentials (ASEPs) began 14.3 +/- 6.6 msec before the trigger spike and had amplitudes of 1064 +/- 867 microV. Twenty-three STAs (6%) showed an average synchronous inhibitory potential (ASIP): a hyperpolarization beginning before the trigger spike and reflecting IPSPs produced by a group of local inhibitory cells synchronized with the trigger cell. ASIPs had an onset latency of -5.5 +/- 2.7 msec and amplitude of -589 +/- 502 microV. Combinations of synchronous and postspike potentials were also observed. Successive recordings provided examples of convergent and divergent connections between EC and IC cells. Neuron pairs with depolarizing postsynaptic potentials (PSPs) in the STA yielded peaks in the cross-correlograms of the IC and EC action potentials; the peak area was proportional to the amplitude of the PSP. These data suggest that a significantly larger proportion of cortical neurons interact through synchronous activity than through simple serial interactions; moreover, synchronous excitation affected more widely separated cell pairs than EPSPs and IPSPs, which were seen most often among the closest cells.

Computer simulation of the responses of human motoneurons to composite 1A EPSPS: effects of background firing rate

Two compartmental models of spinal alpha motoneurons were constructed to explore the relationship between background firing rate and response to an excitatory input. The results of these simulations were compared with previous results obtained from human motoneurons and discussed in relation to the current model for repetitively firing human motoneurons. The morphologies and cable parameters of the models were based on two type-identified cat motoneurons previously reported in the literature. Each model included five voltage-dependent channels that were modeled using Hodgkin-Huxley formalism. These included fast Na+ and K+ channels in the initial segment and fast Na+ and K+ channels as well as a slow K+ channel in the soma compartment. The density and rate factors for the slow K+ channel were varied until the models could reproduce single spike AHP parameters for type-identified motoneurons in the cat. Excitatory synaptic conductances were distributed along the equivalent dendrites with the same density described for la synapses from muscle spindles to type-identified cat motoneurons. Simultaneous activation of all synapses on the dendrite resulted in a large compound excitatory postsynaptic potential (EPSP). Brief depolarizing pulses injected into a compartment of the equivalent dendrite resulted in pulse potentials (PPs), which resembled the compound EPSPs. The effects of compound EPSPs and PPs on firing probability of the two motoneuron models were examined during rhythmic firing. Peristimulus time histograms, constructed between the stimulus and the spikes of the model motoneuron, showed excitatory peaks whose integrated time course approximated the time course of the underlying EPSP or PP as has been shown in cat motoneurons. The excitatory peaks were quantified in terms of response probability, and the relationship between background firing rate and response probability was explored. As in real human motoneurons, the models exhibited an inverse relationship between response probability and background firing rate. The biophysical properties responsible for the relationship between response probability and firing rate included the shapes of the membrane voltage trajectories between spikes and nonlinear changes in PP amplitude during the interspike interval at different firing rates. The results from these simulations suggest that the relationship between response probability and background firing rate is an intrinsic feature of motoneurons. The similarity of the results from the models, which were based on the properties of cat motoneurons, and those from human motoneurons suggests that the biophysical properties governing rhythmic firing in human motoneurons are similar to those of the cat.

Synaptic interactions between primate precentral cortex neurons revealed by spike-triggered averaging of intracellular membrane potentials in vivo

To document synaptic interactions between neurons in the precentral cortex of macaque monkeys, we recorded in vivo the intracellular (IC) membrane potentials of cortical neurons simultaneously with extracellular (EC) action potentials of neighboring cells. The synaptic potentials correlated with EC spikes were obtained by spike-triggered averages (STA) of the IC membrane potentials for 373 cell pairs recorded in anesthetized and awake behaving monkeys. Sixty-three STAs (17%) showed excitatory postsynaptic potentials (EPSPs), beginning after the trigger spike. Pure EPSPs had onset latencies of 0.9 +/- 0.7 msec (mean +/- SD) and amplitudes of 226 +/- 130 microV. Sixteen STAs (4%) showed postspike inhibitory postsynaptic potentials (IPSPs), with onset latencies of 0.4 +/- 0.4 msec and amplitudes of -274 +/- 188 microV. The most common waveform, observed in 82% of the STAs with features, was a broad depolarization straddling the trigger spikes, reflecting synchronized synaptic input to both IC and EC neurons. These average synchronous excitation potentials (ASEPs) began 14.3 +/- 6.6 msec before the trigger spike and had amplitudes of 1064 +/- 867 microV. Twenty-three STAs (6%) showed an average synchronous inhibitory potential (ASIP): a hyperpolarization beginning before the trigger spike and reflecting IPSPs produced by a group of local inhibitory cells synchronized with the trigger cell. ASIPs had an onset latency of -5.5 +/- 2.7 msec and amplitude of -589 +/- 502 microV. Combinations of synchronous and postspike potentials were also observed. Successive recordings provided examples of convergent and divergent connections between EC and IC cells. Neuron pairs with depolarizing postsynaptic potentials (PSPs) in the STA yielded peaks in the cross-correlograms of the IC and EC action potentials; the peak area was proportional to the amplitude of the PSP. These data suggest that a significantly larger proportion of cortical neurons interact through synchronous activity than through simple serial interactions; moreover, synchronous excitation affected more widely separated cell pairs than EPSPs and IPSPs, which were seen most often among the closest cells.

Effects of background noise on the response of rat and cat motoneurones to excitatory current transients

1. We studied the responses of rat hypoglossal motoneurones to excitatory current transients (ECTs) using a brainstem slice preparation. Steady, repetitive discharge at rates of 12-25 impulses s-1 was elicited from the motoneurones by injecting long (40 s) steps of constant current. Poisson trains of the ECTs were superimposed on these steps. The effects of additional synaptic noise was simulated by adding a zero-mean random process to the stimuli. 2. We measured the effects of the ECTs on motoneurone discharge probability by compiling peristimulus time histograms (PSTHs) between the times of occurrence of the ECTs and the motoneurone spikes. The ECTs produced modulation of motoneurone discharge similar to that produced by excitatory postsynaptic currents. 3. The addition of noise altered the pattern of the motoneurone response to the current transients: both the amplitude and the area of the PSTH peaks decreased as the power of the superimposed noise was increased. Noise tended to reduce the efficacy of the ECTs, particularly when the motoneurones were firing at lower frequencies. Although noise also increased the firing frequency of the motoneurones slightly, the effects of noise on ECT efficacy did not simply result from noise-induced changes in mean firing rate. 4. A modified version of the experimental protocol was performed in lumbar motoneurones of intact, pentobarbitone-anaesthetized cats. These recordings yielded results similar to those obtained in rat hypoglossal motoneurones in vitro. 5. Our results suggest that the presence of concurrent synaptic inputs reduces the efficacy of any one input. The implications of this change in efficacy and the possible underlying mechanisms are discussed.

Role of reticular activation in the modulation of intracortical synchronization

During aroused states of the brain, electroencephalographic activity is characterized by fast, irregular fluctuations of low amplitude, which are thought to reflect desynchronization of neuronal activity. This phenomenon seems at odds with the proposal that synchronization of cortical responses may play an important role in the processing of sensory signals. Here, activation of the mesencephalic reticular formation (MRF), an effective way to "desynchronize the electroencephalogram," was shown to facilitate oscillatory activity in the gamma frequency range and to enhance the stimulus-specific synchronization of neuronal spike responses in the visual cortex of cats.

Abnormal excitability of the corticospinal pathway in patients with amyotrophic lateral sclerosis: a single motor unit study using transcranial magnetic stimulation

The pathophysiology of corticospinal tract degeneration in amyotrophic lateral sclerosis (ALS) was investigated by studying the effect of transcranial magnetic stimulation on discharge characteristics of single motor units during voluntary activation. The motor units were recorded from the first dorsal interosseus muscles of 12 patients with ALS, 14 healthy subjects, 12 patients with upper motor neuron lesions and 9 with pure lower motor neuron diseases. More than 100 magnetic stimuli were delivered over the scalp during minimal muscle contraction. The occurrence of motor unit discharges was plotted in a peristimulus time histogram. An increase in discharge probability at latencies of 20-30 msec, that represents monosynaptic activation (primary peak) was found in normal units. Motor units from ALS patients with short disease durations had significantly increased discharge probabilities in the primary peak (P < 0.001). Motor units from 4 ALS patients with upper motor neuron signs showed double primary peaks: an initial synchronized peak followed by a dispersed peak. The latter was ascribed to a slow corticospinal pathway, which remains undetected or is functionally insignificant in healthy subjects. We conclude that the excitabilities of the surviving corticospinal tract pathways are abnormally increased in ALS, especially in the early stage.