Properties of human motoneurones and their synaptic noise deduced from motor unit recordings with the aid of computer modelling.

This paper reviews two new facets of the behaviour of human motoneurones; these were demonstrated by modelling combined with analysis of long periods of low-frequency tonic motor unit firing (sub-primary range). 1) A novel transformation of the interval histogram has shown that the effective part of the membrane's post-spike voltage trajectory is a segment of an exponential (rather than linear), with most spikes being triggered by synaptic noise before the mean potential reaches threshold. The curvature of the motoneurone's trajectory affects virtually all measures of its behaviour and response to stimulation. The 'trajectory' is measured from threshold, and so includes any changes in threshold during the interspike interval. 2) A novel rhythmic stimulus (amplitude-modulated pulsed vibration) has been used to show that the motoneurone produces appreciable phase-advance during sinusoidal excitation. At low frequencies, the advance increases with rising stimulus frequency but then, slightly below the motoneurones mean firing rate, it suddenly becomes smaller. The gain has a maximum for stimuli at the mean firing rate (the 'carrier'). Such behaviour is functionally important since it affects the motoneurone's response to any rhythmic input, whether generated peripherally by the receptors (as in tremor) or by the CNS (as with cortical oscillations). Low mean firing rates favour tremor, since the high gain and reduced phase advance at the 'carrier' reduce the stability of the stretch reflex.     

Spike discharges of skeletomotor neurons during random noise modulated transmembrane current stimulation and muscle stretch

 Spike discharges of skeletomotor neurons innervating triceps surae muscles elicited by white noise modulated transmembrane current stimulation and muscle stretch were studied in decerebrated cats. The white noise modulated current intensity ranged from 4.3 to 63.2 nA peak-to-peak, while muscle stretches ranged from 100 μm to 4.26 mm peak-to-peak. The neuronal responses were studied by averaging the muscle length records centered at the skeletomotor action potentials (peri-spike average, PSA) and by Wiener analysis. Skeletomotor spikes appeared after a sharp peak in PSA of the injected current, preceded by a longer-lasting smaller wavelet of either depolarizing or hyperpolarizing direction. The PSA amplitude was not related to the injected current amplitude nor showed any differences related to the motor unit type. The PSA amplitudes were virtually independent of the stretching amplitude σ, after an initial increase with stretching amplitudes in the range of 15–40 μm (S.D.), or 100–270 μm peak-to-peak.Analyses of cross-spectra indicated a small or absent increase in gain with frequency in response to injected current, but about 20 dB/decade in the range 10–100 Hz in response to muscle stretch. The peaks of both Wiener kernels in response to current injection appear to decrease with the amplitude of injected current, but this decrease was not statistically significant. The narrow first-order kernels suggest that the transfer function between the current input and spike discharge is lowpass with a wide passband, i.e. there is very little change in dynamics. The values of the second-order kernels appear to be nonzero only along the main diagonal. This is characteristic of a simple Hammerstein type cascade, i.e. a zero memory nonlinearity followed by a linear system. Small values of second-order kernels away from the origin and narrow first-order kernels suggest that the linear cascade contributes very little to the overall dynamic response.In contrast to Wiener kernels found in response to current injection, the Wiener kernels in response to stretch showed a decreasing trend with stretch amplitude. The size of the second-order kernels decreased to a somewhat larger extent with input amplitude than that of the first-order kernels, indicating an amplitude-dependent nonlinearity. Overall, the transformation between length and spike output was described as an LNNL cascade with second-order nonlinearities. Received: 1 April 1993/Accepted in revised form: 24 March 1994     

White noise analysis of graded response in a wind-sensitive,nonspiking interneuron of the cockroach

1. A novel approach using a Gaussian white noise as stimulus is described which allowed quantitative analysis of neuronal responses in the cercal system of the cockroach, Periplaneta americana. Cerci were stimulated by air displacement which was modulated by a sinusoidal and a white noise signal. During the stimulation, intracellular recordings were made from a uniquely identifiable, nonspiking, local interneuron which locates within the terminal abdominal ganglion. The white noise stimulation was cross-correlated with the evoked response to compute first- and second-order kernels that could define the cell's response dynamics. 2. The interneuron, cell 101, has an exceptionally large transverse neurite that connects two asymmetrical dendritic arborizations located on both sides of the ganglion. 3. The first-order Wiener kernels in cell 101 were biphasic (differentiating). The waveforms of the kernels produced by the ipsilateral and contralateral stimulations were roughly mirror images of each other: the kernels produced by wind stimuli on the side ipsilateral to the cell body of the interneuron are initially depolarized and then hyperpolarized, whereas those on the other side are initially hyperpolarized. The polarity reversal occurred along the midline of the animal's body, and no well-defined kernel was produced by a stimulus directed head on or from the tail. 4. Mean square error (MSE) between the actual response and the model prediction suggests that the linear component in cell 101 comprises half of the cell's total response (MSEs for the linear models were about 50% at preferred directions), whereas the second-order, non-linear component is insignificant. The linear component of the wind-evoked response was bandpass with the preferred frequency of 70-90 Hz. 5. Accounting for a noise, we reasonably assumed that at high frequencies the graded response in cell 101 is linearly related to a modulation of the air displacement and sensitive to the rate of change of the signal (i.e., wind velocity) and the direction of its source. It is suggested that the dynamics of the first-order kernel simply reflect the dynamics of sensory receptors that respond linearly to wind stimulation.     

White-noise analysis of nonlinear behavior in an insect sensory neuron: Kernel and cascade approaches

A functional expansion was used to model the relationship between a Gaussian white noise stimulus current and the resulting action potential output in the single sensory neuron of the cockroach femoral tactile spine. A new precise procedure was used to measure the kernels of the functional expansion. Very similar kernel estimates were obtained from separate sections of the data produced by the same neuron with the same input noise power level, although some small time-varying effects were detectable in moving through the data. Similar kernel estimates were measured using different input noise power levels for a given cell, or when comparing different cells under similar stimulus conditions. The kernels were used to identify a model for sensory encoding in the neuron, comprising a cascade of dynamic linear, static nonlinear, and dynamic linear elements. Only a single slice of the estimated experimental second-order kernel was used in identifying the cascade model. However, the complete second-order kernel of the cascade model closely resembled the estimated experimental kernel. Moreover, the model could closely predict the experimental action potential train obtained with novel white noise inputs.     

Nonlinear analysis of the human visual evoked response

Using time-domain correlation techniques, the first- and second-order Wiener kernels have been calculated for the system mediating the human visual evoked response. The first-order kernels indicate the linear element is a resonant one, with a natural frequency near 20 Hz, and a memory of approximately 250 ms. The transport delay associated with this element is approximately 56 ms. The second-order kernels indicate a quadratic nonlinear element with a memory less than 20 ms. The analytic form of this element can be approximated by a parabola shifted to the right of the origin. A close correspondance between the spectrum of the first-order kernel and the spectrum of the main diagonal of the second-order kernel suggests the nonlinear element preceeds the linear one. Tests of reproducibility on the first-order kernel and the main diagonal of the second-order kernel suggest they are reliable describing functions for the system mediating the human visual evoked response.     

White noise analysis of a chromatic type horizontal cell in the Xenopus retina

The dynamics of color-coded signal transmission in the light-adapted Xenopus retina were studied by a combination of white noise (Wiener) analysis and simultaneous recordings from two types of horizontal cells: chromatic-type horizontal cells (C-HCs) are hyperpolarized by blue light and depolarized by red light, whereas luminosity-type horizontal cells (L-HCs) are hyperpolarized by all wave-lengths. The retina was stimulated by two superimposed fields of red and blue light modulated by two independent white noise signals, and the resulting intracellular responses were decomposed into red and blue components (first-order kernels). The first-order kernels predict the intracellular responses with a small degree of error (3.5-9.5% in terms of mean square error) under conditions where modulated responses exceeded 30 mV in amplitude peak-to-peak, thus demonstrating that both red and blue modulation responses are linear. Moreover, there is little or no interaction between the red- and blue-evoked responses; i.e., nearly identical first-order kernels were obtained for one color whether the other color was modulated or not. In C-HCs (but not L-HCs), there were consistent differences in the dynamics of the red and blue responses. In the C-HC, the cutoff frequency of the red response was higher than for the blue (approximately 12 vs 5 Hz), and the red kernel was more bandpass than the blue. In the L-HC, kernel waveform and cutoff frequencies were similar for both colors (approximately 12 Hz or greater), and the time-to-peak of the L-HC kernel was always shorter than either the red or blue C-HC kernel. These results have implications for the mechanisms underlying color coding in the distal retina, and they further suggest that nonlinear phenomena, such as voltage-dependent conductances in HCs, do not contribute to the generation of modulation responses under the experimental conditions used here.     

The transducer and encoder of frog muscle spindles are essentially nonlinear. Physiological conclusions from a white-noise analysis

1. Nonlinear second order white-noise analysis has been applied to the isolated frog muscle spindle. Power (²) of the Gaussian white noise (GWN) and the average prestretch level L were varied and the response of both the isolated receptor potential (transducer) and the action potential (encoder) level were analysed. 2. The standard white-noise method is briefly presented. Particular emphasis, however, is put on the limitations in the range of validity of the method and, consequently, on the use and interpretation of the kernels as a Wiener model. Conclusions in the present paper are within this frame and are mainly of qualitative nature. 3. The analysis reveals that the nonlinear contributions of the model are essential for approximating physiological results, thus ruling out purely linear modelling for this receptor organ. 4. The dependence of the transducer kernels on are compatible with the behaviour of a rectifier. Rectification is represented by the lack of hyperpolarization within the isolated receptor potential and is enhanced by the substantial memory in the linear and nonlinear kernels as demonstrated by their extent in time. This is equivalent to low power in high frequencies of the response. Obviously, the hyperpolarizing potentials following each spike counteract the long transducer memory. 5. At the encoder level the memory of the system is strongly reduced. This is achieved by using predominantly high frequency components of the receptor potential for triggering the process of impulse generation, and by the precise coupling and high frequency content of the impulses. This coupling precision is possible because of the sensitivity of the spike-generating mechanism to steep rising transients of the receptor potential and also owing to the reduction in transducer memory by the hyperpolarizing afferpotentials. 6. The preference given to the high frequency components is also read from the structure of the second order transducer kernel and from both the linear and the second order encoder kernels, which allows the most effective input waveform for triggering action potentials to be determined. 6. When the operating point is changed to higher prestretch values, kernel heights increase strongly implying higher response strength of the muscle spindle. The kernel structure is changed as well in the direction of reducing the effective memory already at the level of the receptor potentials, probably a means to prevent too high depolarization values.     

Membrane potential changes of skeletomotor neurons in response to random stretches of the triceps surae muscles in decerebrate cats

 The properties of membrane potential changes of skeletomotor neurons (S, FR, and FF) innervating triceps surae muscles during pseudorandom stretching of these muscles were studied in decerebrate cats. Peak amplitudes of pseudorandom muscle stretches ranged from 119 μm to 4.15 mm peak-to-peak. Sequences of ten identical stretching periods were applied for averaging. Shapes of membrane potential changes and probability density distribution of amplitudes of the input and output signals and power spectra suggest that the skeletomotor neuron membrane has nonlinear properties. First- and second-order Wiener kernels were determined by applying the cross-correlation (Lee-Schetzen) method. The results suggest that the transfer function between muscle stretches and subthreshold membrane potentials is a Wiener-type cascade. This cascade is consistent with a linear, second-order, underdamped transfer function followed by a simple quadratic nonlinearity [linear (L) system followed by nonlinear (N) system, or LN cascade]. Including the nonlinear component calculated from the second-order Wiener kernel improved the model significantly over its linear counterpart, especially in S-type motoneurons. Qualitatively similar results were obtained with all types of motoneurons studied. Received: 1 April 1993/Accepted in revised form: 24 March 1994     

Intensity Coding in the Frog Retina : Quantitative relations between impulse and graded activity

The impulse discharge of single on-off neurons and a graded field potential, the proximal negative response (PNR), were simultaneously recorded with an extracellular microelectrode in the inner frog retina. Normalized amplitude-intensity functions for the on-response of the PNR and the neuron's post-stimulus time histogram (PSTH) were nearly coincident and typically showed a dynamic range spanning approximately 2 log units of intensity. Thus a nearly linear relation is found between the amplitude of the PNR and the neuron's PSTH. A neuron's PSTH amplitude and maximum instantaneous frequency of discharge were usually highly correlated, but occasional marked disparities indicate that temporal jitter of the first spike latency is an additional, relatively independent variable influencing PSTH amplitude. It typically changes by a factor of 20–30 over the intensity range. These and other findings have implications for the functional significance of the PNR and the PSTH, for a possible linear link between amacrine and on-off ganglion cells, and for a mechanism of intensity coding in which temporal jitter of latency exerts a major role.     

Retinal responses to white-noise modulated current stimuli

Electric current was injected into a rabbit's eye with white-noise modulations of the current amplitude. A variable D.C. bias was added to the whitenoise stimulus to study the effects of stimulus bias. For each bias level, the ERG response to the electrical stimulus was cross-correlated with the random stimulus to estimate first-and second-order Wiener kernels. The kernels indicated both linear and nonlinear characteristics of the Electrical ERG. The results for the zero biased stimulus are particularly relevant for clinical testing because the root mean square (RMS) level of the stimulus was less than 0.2 mA.This work was supported by NIH Grant Number EY03022     

Single unit activity in the guinea-pig cochlear nucleus during sleep and wakefulness.

The effects of waking and sleep on the response properties of auditory units in the ventral cochlear nucleus (CN) were explored by using extracellular recordings in chronic guinea-pigs. Significant increases and decreases in firing rate were detected in two neuronal groups, a) the "sound-responding" and b) the "spontaneous" (units that do not show responses to any acoustic stimuli controlled by the experimenter). The "spontaneous" may be considered as belonging to the auditory system because the corresponding units showed a suppression of their discharge when the receptor was destroyed. The auditory CN units were characterized by their PSTH in response to tones at their characteristic frequency and also by the changes in firing rate and probability of discharge evaluated during periods of waking, slow wave and paradoxical sleep. The CNS performs functions dependent on sensory inputs during wakefulness and sleep phases. By studying the auditory input at the level of the ventral CN with constant sound stimuli, it was shown that, in addition to the firing rate shifts, some units presented changes in the temporal probability of discharge, implying central actions on the corresponding neurons. The mean latency of the responses, however, did not show significant changes throughout the sleep-waking cycle. The auditory efferent pathways are postulated to modulate the auditory input at CN level during different animal states. The probability of firing and the changes in the temporal pattern, as shown by the PSTH, are thus dependent on both the auditory input and the functional brain state related to the sleep-waking cycle.     

System analysis of Phycomyces light-growth response: madC, madG, and madH mutants.

The light-growth response of Phycomyces has been studied further with the sum-of-sinusoids method in the framework of the Wiener theory of nonlinear system identification. The response was treated as a black box with the logarithm of light intensity as the input and elongation rate as the output. The nonlinear input-output relation of the light-growth response can be represented mathematically by a set of weighting functions called kernels, which appear in the Wiener intergral series. The linear (first-order) kernels of wild type, and of single and double mutants affected in genes madA to madG were determined previously with Gaussian white noise test stimuli, and were used to investigate the interactions among the products of these genes (R.C. Poe, P. Pratap, and E.D. Lipson. 1986. Biol. Cybern. 55:105.). We have used the more precise sum-of-sinusoids method to extend the interaction studies, including both the first- and second-order kernels. Specifically, we have investigated interactions of the madH ("hypertropic") gene product with the madC ("night blind") and madG ("stiff") gene products. Experiments were performed on the Phycomyces tracking machine. The log-mean intensity of the stimulus was 6 x 10(-2) W m-2 and the wavelength was 477 nm. The first- and second-order kernels were analyzed in terms of nonlinear kinetic models.(ABSTRACT TRUNCATED AT 250 WORDS)     

Excitability of single firing human motoneurones to single and repetitive stimulation (experiment and model)

The activity of single motoneurones of m. flexor carpi ulnaris (FCU) was investigated by recording their motor unit (MU) action potentials during weak and moderate voluntary muscle contractions. The MU firing rate range was 4.5-15 imp/s. The excitability of motoneurones was tested with a number of single stimuli eliciting a monosynaptic H-reflex of low amplitude. Two different indices were defined which relate to motoneuronal excitability: the response index--the ratio of the number of responses of a motoneurone to the total number of stimuli, and the response time--the time after the last background MU discharge at which motoneurone is ready to respond to the excitatory volley. Both the response index and the response time were determined for single motoneurones at different levels of background activity. In the lower range of firing rates, the response index for all motoneurones decreased when increasing the firing rate, but it remained constant in the higher rate range. This kind of response seems to be a typical motoneuronal response to the stimulation with single stimuli. The data on the response time were used to study the excitability of the same single motoneurones to computer simulated repetitive stimulation (stimulation rate 40-100 imp/s). In this case, the excitability of each motoneurone was determined as an increment of its firing rate in response to the stimulation. For the lower firing rate range, the excitability for all motoneurones also decreased when the firing rates increased whereas a variety of slopes was obtained in the higher rate range.(ABSTRACT TRUNCATED AT 250 WORDS)     

Modelling of the control of heart rate by breathing using a kernel method

The process of the breathing (input) to the heart rate (output) of man is considered for system identification by the input-output relationship, using a mathematical model expressed as integral equations. The integral equation is considered and fixed so that the identification method reduces to the determination of the values within the integral, called kernels, resulting in an integral equation whose input-output behaviour is nearly identical to that of the system. This paper uses an algorithm of kernel identification of the Volterra series which greatly reduces the computational burden and eliminates the restriction of using white Gaussian input as a test signal. A second-order model is the most appropriate for a good estimate of the system dynamics. The model contains the linear part (first-order kernel) and quadratic part (second-order kernel) in parallel, and so allows for the possibility of separation between the linear and non-linear elements of the process. The response of the linear term exhibits the oscillatory input and underdamped nature of the system. The application of breathing as input to the system produces an oscillatory term which may be attributed to the nature of sinus node of the heart being sensitive to the modulating signal the breathing wave. The negative-on diagonal seems to cause the dynamic asymmetry of the total response of the system which opposes the oscillatory nature of the first kernel related to the restraining force present in the respiratory heart rate system. The presence of the positive-off diagonal of the second-order kernel of respiratory control of heart rate is an indication of an escape-like phenomenon in the system.     

Synaptic integration in spinal motoneurones.

Spinal motoneurones receive thousands of presynaptic excitatory and inhibitory synaptic contacts distributed throughout their dendritic trees. Despite this extensive convergence, there have been very few studies of how synaptic inputs interact in mammalian motoneurones when they are activated concurrently. In the experiments reported here, we measured the effective synaptic currents and the changes in firing rate evoked in cat spinal motoneurones by concurrent repetitive activation of two separate sets of presynaptic neurons. We compared these effects to those predicted by a linear sum of the effects produced by activating each set of presynaptic neurons separately. We generally found that when two inputs were activated concurrently, both the effective synaptic currents and the synaptically-evoked changes in firing rate they produced in motoneurones were generally linear, or slightly less than the linear sum of the effects produced by activating each input alone. The results suggest that the spatial distribution synaptic terminals on the dendritic trees of motoneurones may help isolate synapses from one another, minimizing non-linear interactions.     

On-line identification of sensory systems using pseudorandom binary noise perturbations.

A technique of on-line identification of linear system characteristics from sensory systems with spike train or analog voltage outputs was developed and applied to the semicircular canal. A pseudorandom binary white noise input was cross-correlated with the system's output to produce estimates of linear system unit impulse responses (UIRs), which were then corrected for response errors of the input transducers. The effects of variability in the system response characteristics and sensitivity were studied by employing the technique with known linear analog circuits. First-order unit afferent responses from the guitarfish horizontal semicircular canal were cross-correlated with white noise rotational acceleration inputs to produce non-parametric UIR models. In addition, the UIRs were fitted by nonlinear regression to truncated exponential series to produce parametric models in the form of low-order linear system equations. The experimental responses to the white noise input were then compared with those predicted from the UIR models linear convolution, and the differences were expressed as a percent mean-square-error (%MSE). The average difference found from a population of 62 semicircular canal afferents was relatively low mean and standard deviation of 10.2 +/- 5.9 SD%MSE, respectively. This suggests that relatively accurate inferences can be made concerning the physiology of the semicircular canal from the linear characteristics of afferent responses.