Statistical limits for detecting change in the cumulative sum derivative of the peristimulus time histogram

Substantia Nigra Dopamine Neuronal Responses to Habenular Stimulation and Foot Shock Are Altered by Lesions of the Rostromedial Tegmental Nucleus

Dopamine (DA) neurons of the substantia nigra (SN) and ventral tegmental area generally respond to aversive stimuli or the absence of expected rewards with transient inhibition of firing rates, which can be recapitulated with activation of the lateral habenula (LHb) and eliminated by lesioning the intermediating rostromedial tegmental nucleus (RMTg). However, a minority of DA neurons respond to aversive stimuli, such as foot shock, with a transient increase in firing rate, an outcome that rarely occurs with LHb stimulation. The degree to which individual neurons respond to these two stimulation modalities with the same response phenotype and the role of the RMTg is not known. Here, we record responses from single SN DA neurons to alternating activation of the LHb and foot shock in male rats. Lesions of the RMTg resulted in a shift away from inhibition to no response during both foot shock and LHb stimulation. Furthermore, lesions unmasked an excitatory response during LHb stimulation. The response correspondence within the same neuron between the two activation sources was no different from chance in sham controls, suggesting that external inputs rather than intrinsic DA neuronal properties are more important to response outcome. These findings contribute to a literature that shows a complex neurocircuitry underlies the regulation of DA activity and, by extension, behaviors related to learning, anhedonia, and cognition.

Breathing-driven prefrontal oscillations regulate maintenance of conditioned-fear evoked freezing independently of initiation

Brain-body interactions are thought to be essential in emotions but their physiological basis remains poorly understood. In mice, regular 4 Hz breathing appears during freezing after cue-fear conditioning. Here we show that the olfactory bulb (OB) transmits this rhythm to the dorsomedial prefrontal cortex (dmPFC) where it organizes neural activity. Reduction of the respiratory-related 4 Hz oscillation, via bulbectomy or optogenetic perturbation of the OB, reduces freezing. Behavioural modelling shows that this is due to a specific reduction in freezing maintenance without impacting its initiation, thus dissociating these two phenomena. dmPFC LFP and firing patterns support the region's specific function in freezing maintenance. In particular, population analysis reveals that network activity tracks 4 Hz power dynamics during freezing and reaches a stable state at 4 Hz peak that lasts until freezing termination. These results provide a potential mechanism and a functional role for bodily feedback in emotions and therefore shed light on the historical James-Cannon debate.

Pauses in cholinergic interneuron firing exert an inhibitory control on striatal output in vivo

The cholinergic interneurons (CINs) of the striatum are crucial for normal motor and behavioral functions of the basal ganglia. Striatal CINs exhibit tonic firing punctuated by distinct pauses. Pauses occur in response to motivationally significant events, but their function is unknown. Here we investigated the effects of pauses in CIN firing on spiny projection neurons (SPNs) – the output neurons of the striatum – using in vivo whole cell and juxtacellular recordings in mice. We found that optogenetically-induced pauses in CIN firing inhibited subthreshold membrane potential activity and decreased firing of SPNs. During pauses, SPN membrane potential fluctuations became more hyperpolarized and UP state durations became shorter. In addition, short-term plasticity of corticostriatal inputs was decreased during pauses. Our results indicate that, in vivo, the net effect of the pause in CIN firing on SPNs activity is inhibition and provide a novel mechanism for cholinergic control of striatal output.

Nerve cells or neurons communicate with one another using electrical impulses and chemical messengers called neurotransmitters. Additional molecules known as neuromodulators regulate the communication process. In contrast to neurotransmitters, neuromodulators do not send messages directly from one neuron to the next. Instead they change the way that neurons respond to neurotransmitters.

For example, the neuromodulator acetylcholine is most abundant in a region called the striatum. Located deep within the brain, the striatum contributes to learning and memory, motivation, and movement. Studies in rodents show that neurons within the striatum called cholinergic interneurons are almost continuously active. Each time these cells fire, they release acetylcholine. But whenever an animal experiences something unusual or important, the interneurons temporarily stop firing. Zucca et al. wanted to know whether these pauses in firing also act as a signal within the striatum.

To find out, Zucca et al. inserted a light-sensitive ion channel into cholinergic interneurons in the mouse striatum. Activating the ion channels with a laser beam stopped the interneurons from firing. Zucca et al. showed that these pauses in firing reduced the activity of another group of neurons, the spiny projection neurons. These are the major output neurons of the striatum. They send messages from the striatum to other parts of the brain. The results thus suggest that cholinergic interneurons signal notable events by temporarily blocking output from the striatum.

Understanding how cholinergic interneurons work will help reveal how the striatum drives behavior. It may also lead to treatments for diseases caused by cholinergic system dysfunction. Many patients with Parkinson’s disease or schizophrenia take medicines to block the effects of acetylcholine. Understanding how acetylcholine affects the striatum may help clarify how these treatments work.

A nociceptive stress model of adolescent physical abuse induces contextual fear and cingulate nociceptive neuroplasticities

Adolescent physical abuse impairs emotional development and evokes cingulate pathologies, but its neuronal and circuit substrates are unknown. Conditioning adolescent rabbits with noxious colorectal distension for only 2 h over 3 weeks simulated the human child abuse in amplitude, frequency, and duration. Thermal withdrawal thresholds were unchanged suggesting that sensitized spinal mechanisms may not be operable. Unchanged weight, stools, colorectal histology, and no evidence of abdominal pain argue against tissue injury or irritable bowel syndrome. Contextual fear was amplified as they avoided the site of their abuse. Conditioning impacted anterior cingulate and anterior midcingulate (ACC, aMCC) neuron excitability: (1) more neurons responded to cutaneous and visceral (VNox) noxious stimuli than controls engaging latent nociception (present but not manifest in controls). (2) Rear paw stimulation increased responses over forepaws with shorter onsets and longer durations, while forepaw responses were of higher amplitude. (3) There were more VNox responses with two excitatory phases and longer durations. (4) Some had unique three-phase excitatory responses. (5) Long-duration VNox stimuli did not inhibit neurons as in controls, suggesting the release of an inhibitory circuit. (6) aMCC changes in cutaneous but not visceral nociception confirmed its role in cutaneous nociception. For the first time, we report neuroplasticities that may be evoked by adolescent physical abuse and reflect psychogenic pain: i.e., no ongoing peripheral pain and altered ACC nociception. These limbic responses may be a cognitive trace of abuse and may shed light on impaired human emotional development and sexual function.

Single and Multiple Change Point Detection in Spike Trains: Comparison of Different CUSUM Methods

In a natural environment, sensory systems are faced with ever-changing stimuli that can occur, disappear or change their properties at any time. For the animal to react adequately the sensory systems must be able to detect changes in external stimuli based on its neuronal responses. Since the nervous system has no prior knowledge of the stimulus timing, changes in stimulus need to be inferred from the changes in neuronal activity, in particular increase or decrease of the spike rate, its variability, and shifted response latencies. From a mathematical point of view, this problem can be rephrased as detecting changes of statistical properties in a time series. In neuroscience, the CUSUM (cumulative sum) method has been applied to recorded neuronal responses for detecting a single stimulus change. Here, we investigate the applicability of the CUSUM approach for detecting single as well as multiple stimulus changes that induce increases or decreases in neuronal activity. Like the nervous system, our algorithm relies exclusively on previous neuronal population activities, without using knowledge about the timing or number of external stimulus changes. We apply our change point detection methods to experimental data obtained by multi-electrode recordings from turtle retinal ganglion cells, which react to changes in light stimulation with a range of typical neuronal activity patterns. We systematically examine how variations of mathematical assumptions (Poisson, Gaussian, and Gamma distributions) used for the algorithms may affect the detection of an unknown number of stimulus changes in our data and compare these CUSUM methods with the standard Rate Change method. Our results suggest which versions of the CUSUM algorithm could be useful for different types of specific data sets.

Estimating reflex responses in large populations of motor units by decomposition of the high-density surface electromyogram

KEY POINTS

Reflex responses of single motor units have been used for the study of spinal circuitries but the methods employed are invasive and limited to the assessment of a relatively small number of motor units. We propose a new approach to investigate reflexes on individual motor units based on high-density surface electromyography (HDsEMG) decomposition. The decomposition of HDsEMG has been previously validated in voluntary isometric contractions but never during reflex activities. The use of HDsEMG decomposition for reflex studies at the individual motor unit level, during constant force contractions, with excitatory and inhibitory stimuli, was validated here by the comparison of results with concurrently recorded intramuscular EMG signals. The validation results showed that HDsEMG decomposition allows an accurate quantification of reflex responses for a large number of individual motor units non-invasively, for both excitatory and inhibitory stimuli.

ABSTRACT

We propose and validate a non-invasive method that enables accurate detection of the discharge times of a relatively large number of motor units during excitatory and inhibitory reflex stimulations. High-density surface electromyography (HDsEMG) and intramuscular EMG (iEMG) were recorded from the tibialis anterior muscle during ankle dorsiflexions performed at 5%, 10% and 20% of the maximum voluntary contraction (MVC) force, in nine healthy subjects. The tibial nerve (inhibitory reflex) and the peroneal nerve (excitatory reflex) were stimulated with constant current stimuli. In total, 416 motor units were identified from the automatic decomposition of the HDsEMG. The iEMG was decomposed using a state-of-the-art decomposition tool and provided 84 motor units (average of two recording sites). The reflex responses of the detected motor units were analysed using the peri-stimulus time histogram (PSTH) and the peri-stimulus frequencygram (PSF). The reflex responses of the common motor units identified concurrently from the HDsEMG and the iEMG signals showed an average disagreement (the difference between number of observed spikes in each bin relative to the mean) of 8.2 ± 2.2% (5% MVC), 6.8 ± 1.0% (10% MVC) and 7.5 ± 2.2% (20% MVC), for reflex inhibition, and 6.5 ± 4.1%, 12.0 ± 1.8% and 13.9 ± 2.4%, for reflex excitation. There was no significant difference between the characteristics of the reflex responses, such as latency, amplitude and duration, for the motor units identified by both techniques. Finally, reflex responses could be identified at higher force (4 of the 9 subjects performed contraction up to 50% MVC) using HDsEMG but not iEMG, because of the difficulty in decomposing the iEMG at high forces. In conclusion, single motor unit reflex responses can be estimated accurately and non-invasively in relatively large populations of motor units using HDsEMG. This non-invasive approach may enable a more thorough investigation of the synaptic input distribution on active motor units at various force levels.

Quantification and classification of neuronal responses in kernel-smoothed peristimulus time histograms

Peristimulus time histograms are a widespread form of visualizing neuronal responses. Kernel convolution methods transform these histograms into a smooth, continuous probability density function. This provides an improved estimate of a neuron's actual response envelope. We here develop a classifier, called the h-coefficient, to determine whether time-locked fluctuations in the firing rate of a neuron should be classified as a response or as random noise. Unlike previous approaches, the h-coefficient takes advantage of the more precise response envelope estimation provided by the kernel convolution method. The h-coefficient quantizes the smoothed response envelope and calculates the probability of a response of a given shape to occur by chance. We tested the efficacy of the h-coefficient in a large data set of Monte Carlo simulated smoothed peristimulus time histograms with varying response amplitudes, response durations, trial numbers, and baseline firing rates. Across all these conditions, the h-coefficient significantly outperformed more classical classifiers, with a mean false alarm rate of 0.004 and a mean hit rate of 0.494. We also tested the h-coefficient's performance in a set of neuronal responses recorded in humans. The algorithm behind the h-coefficient provides various opportunities for further adaptation and the flexibility to target specific parameters in a given data set. Our findings confirm that the h-coefficient can provide a conservative and powerful tool for the analysis of peristimulus time histograms with great potential for future development.

Visualization of plasticity in fear-evoked calcium signals in midbrain dopamine neurons

Dopamine is broadly implicated in fear-related processes, yet we know very little about signaling dynamics in these neurons during active fear conditioning. We describe the direct imaging of calcium signals of dopamine neurons during Pavlovian fear conditioning using fiber-optic confocal microscopy coupled with the genetically encoded calcium indicator GCaMP3. We observed calcium transients in a subset of dopamine neurons to an unconditioned fear stimulus on the first day of Pavlovian fear conditioning. On the second day, calcium transients occurred in response to conditioned and unconditioned stimuli. These results demonstrate plasticity in dopamine neuron calcium signals and the occurrence of activity-dependent processes in these neurons during fear conditioning.

Influence of common synaptic input to motor neurons on the neural drive to muscle in essential tremor

Tremor in essential tremor (ET) is generated by pathological oscillations at 4-12 Hz, likely originating at cerebello-thalamo-cortical pathways. However, the way in which tremor is represented in the output of the spinal cord circuitries is largely unknown because of the difficulties in identifying the behavior of individual motor units from tremulous muscles. By using novel methods for the decomposition of multichannel surface EMG, we provide a systematic analysis of the discharge properties of motor units in nine ET patients, with concurrent recordings of EEG activity. This analysis allowed us to infer the contribution of common synaptic inputs to motor neurons in ET. Motor unit short-term synchronization was significantly greater in ET patients than in healthy subjects. Furthermore, the strong association between the degree of synchronization and the peak in coherence between motor unit spike trains at the tremor frequency indicated that the high synchronization levels were generated mainly by common synaptic inputs specifically at the tremor frequency. The coherence between EEG and motor unit spike trains demonstrated the presence of common cortical input to the motor neurons at the tremor frequency. Nonetheless, the strength of this input was uncorrelated to the net common synaptic input at the tremor frequency, suggesting a contribution of spinal afferents or secondary supraspinal pathways in projecting common input at the tremor frequency. These results provide the first systematic analysis of the neural drive to the muscle in ET and elucidate some of its characteristics that determine pathological tremulous muscle activity.

Functional connectivity in the brain estimated by analysis of gamma events

It is known that gamma activity is generated by local networks. In this paper we introduced a new approach for estimation of functional connectivity between neuronal networks by measuring temporal relations between peaks of gamma event amplitudes. We have shown in freely moving rats that gamma events recorded between electrodes 1.5 mm apart in the majority of cases, are generated by different neuronal modules interfering with each other. The map of functional connectivity between brain areas during the resting state, created based on gamma event temporal relationships is in agreement with anatomical connections and with maps described by fMRI methods during the resting state. The transition from the resting state to exploratory activity is accompanied by decreased functional connectivity between most brain areas. Our data suggest that functional connectivity between interhemispheric areas depends on GABAergic transmission, while intrahemispheric functional connectivity is kainate receptor dependent. This approach presents opportunities for merging electrographic and fMRI data on brain functional connectivity in normal and pathological conditions.

Human stretch reflex pathways reexamined

Reflex responses of tibialis anterior motor units to stretch stimuli were investigated in human subjects. Three types of stretch stimuli were applied (tap-like, ramp-and-hold, and half-sine stretch). Stimulus-induced responses in single motor units were analyzed using the classical technique, which involved building average surface electromyogram (SEMG) and peristimulus time histograms (PSTH) from the discharge times of motor units and peristimulus frequencygrams (PSF) from the instantaneous discharge rates of single motor units. With the use of SEMG and PSTH, the tap-like stretch stimulus induced five separate reflex responses, on average. With the same single motor unit data, the PSF technique indicated that the tap stimulus induced only three reflex responses. Similar to the finding using the tap-like stretch stimuli, ramp-and-hold stimuli induced several peaks and troughs in the SEMG and PSTH. The PSF analyses displayed genuine increases in discharge rates underlying the peaks but not underlying the troughs. Half-sine stretch stimuli induced a long-lasting excitation followed by a long-lasting silent period in SEMG and PSTH. The increase in the discharge rate, however, lasted for the entire duration of the stimulus and continued during the silent period. The results are discussed in the light of the fact that the discharge rate of a motoneuron has a strong positive linear association with the effective synaptic current it receives and hence represents changes in the membrane potential more directly and accurately than the other indirect measures. This study suggests that the neuronal pathway of the human stretch reflex does not include inhibitory pathways.

Functional connectivity between brain areas estimated by analysis of gamma waves

The goal of this study is to investigate functional connectivity between different brain regions by analyzing the temporal relationship of the maxima of gamma waves recorded in multiple brain areas. Local field potentials were recorded from motor cortex, hippocampus, entorhinal cortex and piriform cortex of rats. Gamma activity was filtered and separated into two bands; high (65-90Hz) and low (30-55Hz) gamma. Maxima for gamma activity waves were detected and functional connectivity between different brain regions was determined using Shannon entropy for perievent histograms for each pair channels. Significant Shannon entropy values were reported as connectivity factors. We defined a connectivity matrix based the connectivity factors between different regions. We found that maxima of low and high frequency gamma occur in strong temporal relationship between some brain areas, indicating the existence of functional connections between these areas. The spatial pattern of functional connections between brain areas was different for slow wave sleep and waking states. However for each behavioral state in the same animal the pattern of functional connections was stable over time within 30min of continuous analysis and over a 5 day period. With the same electrode montage the pattern of functional connectivity varied from one subject to another. Analysis of the temporal relationship of maxima of gamma waves between various brain areas could be a useful tool for investigation of functional connections between these brain areas. This approach could be applied for analysis of functional alterations occurring in these connections during different behavioral tasks and during processes related to learning and memory. The specificity in the connectivity pattern from one subject to another can be explained by the existence of unique functional networks for each subject.

Trial time warping to discriminate stimulus-related from movement-related neural activity

In tasks where different sensory, cognitive, and motor events are mixed in a sequence it is difficult to determine whether neural activity is related to any behavioral parameter. Here, we consider the case in which two alternative trial-alignment schemes correspond to two different neural representations, one stimulus-related and the other movement-related, using both simulations of neural activity and real recordings in the medial premotor areas during a multiple-interval tapping task called synchronization-continuation task (SCT). To discover whether neural responses are better aligned to sensory or motor events we introduce a family of trial-alignment time-warping functions indexed by a single parameter such that when the parameter takes the value 0 the trials are aligned to the stimulus and when the parameter takes the value 1 they are aligned to the movement. We then characterize neurons by the best-fitting alignment scheme (in the sense of maximum likelihood) under the assumption that the correct alignment would produce homogeneous trials without excess trial-to-trial variation. We use Bayes factors to determine the evidence in favor of sensory or motor neural alignments. The simulations revealed that the variability in neural responses and sequential motor outputs are key parameters to obtain appropriate warping results. In addition, the analysis on the activity of 500 neurons in the medial premotor areas of monkeys executing the SCT showed that most of the neural responses (54.2%) were aligned to the tapping movements instead of the stimuli used to drive the temporal behavior.

Recruitment of motor units in two fascicles of the semispinalis cervicis muscle

This study investigated the behavior of motor units in the semispinalis cervicis muscle. Intramuscular EMG recordings were obtained unilaterally at levels C2 and C5 in 15 healthy volunteers (8 men, 7 women) who performed isometric neck extensions at 5%, 10%, and 20% of the maximal force [maximum voluntary contraction (MVC)] for 2 min each and linearly increasing force contractions from 0 to 30% MVC over 3 s. Individual motor unit action potentials were identified. The discharge rate and interspike interval variability of the motor units in the two locations did not differ. However, the recruitment threshold of motor units detected at C2 (n = 16, mean ± SD: 10.3 ± 6.0% MVC) was greater than that of motor units detected at C5 (n = 92, 6.9 ± 4.3% MVC) (P < 0.01). A significant level of short-term synchronization was identified in 246 of 307 motor unit pairs when computed within one spinal level but only in 28 of 110 pairs of motor units between the two levels. The common input strength, which quantifies motor unit synchronization, was greater for pairs within one level (0.47 ± 0.32) compared with pairs between levels (0.09 ± 0.07) (P < 0.05). In a second experiment on eight healthy subjects, interference EMG was recorded from the same locations during a linearly increasing force contraction from 0 to 40% MVC and showed significantly greater EMG amplitude at C5 than at C2. In conclusion, synaptic input is distributed partly independently and nonuniformly to different fascicles of the semispinalis cervicis muscle.

A modelling study on transmission of the central oscillator in tremor by a motor neuron pool

In spite of decades of intense research, pathological tremors still constitute unknown disorders. This study addresses, based on a multi-scale model, the behavior of an entire pool of motor neurons in tremor, under the hypothesis that tremor is an oscillation of central origin commonly projected to all motor neurons that innervate a muscle. Our results show that under such conditions both paired discharges and enhanced motor neuron synchronization, two of the characteristic landmarks of tremor, emerge. Moreover, coherence and correlation analyses suggest that the central tremor oscillator is transmitted linearly by the motor neuron pool given that a small set (7 or 8) of motor neurons are sampled.

Low-frequency oscillations of the neural drive to the muscle are increased with experimental muscle pain

We investigated the influence of nociceptive stimulation on the accuracy of task execution and motor unit spike trains during low-force isometric contractions. Muscle pain was induced by infusion of hypertonic saline into the abductor digiti minimi muscle of 11 healthy men. Intramuscular EMG signals were recorded from the same muscle during four isometric contractions of 60-s duration at 10% of the maximal force [maximal voluntary contraction (MVC)] performed before injection (baseline), after injection of isotonic (control) or hypertonic saline (pain), and 15 min after pain was no longer reported. Each contraction was preceded by three 3-s ramp contractions from 0% to 10% MVC. The low-frequency oscillations of motor unit spike trains were analyzed by the first principal component of the low-pass filtered spike trains [first common component (FCC)], which represents the effective neural drive to the muscle. Pain decreased the accuracy of task performance [coefficient of variation (CoV) for force: baseline, 2.8 ± 1.8%, pain, 3.9 ± 1.8%; P < 0.05] and reduced motor unit discharge rates [11.6 ± 2.3 pulses per second (pps) vs. 10.7 ± 1.7 pps; P < 0.05]. Motor unit recruitment thresholds (2.2 ± 1.2% MVC vs. 2.4 ± 1.6% MVC), interspike interval variability (18.4 ± 4.9% vs. 19.1 ± 5.4%), strength of motor unit short-term synchronization [common input strength (CIS) 1.02 ± 0.44 vs. 0.83 ± 0.22], and strength of common drive (0.47 ± 0.08 vs. 0.47 ± 0.06) did not change across conditions. The FCC signal was correlated with force (R = 0.45 ± 0.06), and the CoV for FCC increased in the painful condition (5.69 ± 1.29% vs. 7.83 ± 2.61%; P < 0.05). These results indicate that nociceptive stimulation increased the low-frequency variability in synaptic input to motoneurons.

Central chemoreceptor modulation of breathing via multipath tuning in medullary ventrolateral respiratory column circuits

Ventrolateral respiratory column (VRC) circuits that modulate breathing in response to changes in central chemoreceptor drive are incompletely understood. We employed multielectrode arrays and spike train correlation methods to test predictions of the hypothesis that pre-Bötzinger complex (pre-BötC) and retrotrapezoid nucleus/parafacial (RTN-pF) circuits cooperate in chemoreceptor-evoked tuning of ventral respiratory group (VRG) inspiratory neurons. Central chemoreceptors were selectively stimulated by injections of CO(2)-saturated saline into the vertebral artery in seven decerebrate, vagotomized, neuromuscularly blocked, and artificially ventilated cats. Among sampled neurons in the Bötzinger complex (BötC)-to-VRG region, 70% (161 of 231) had a significant change in firing rate after chemoreceptor stimulation, as did 70% (101 of 144) of the RTN-pF neurons. Other responsive neurons (24 BötC-VRG; 11 RTN-pF) had a change in the depth of respiratory modulation without a significant change in average firing rate. Seventy BötC-VRG chemoresponsive neurons triggered 189 offset-feature correlograms (96 peaks; 93 troughs) with at least one responsive BötC-VRG cell. Functional input from at least one RTN-pF cell could be inferred for 45 BötC-VRG neurons (19%). Eleven RTN-pF cells were correlated with more than one BötC-VRG target neuron, providing evidence for divergent connectivity. Thirty-seven RTN-pF neurons, 24 of which were chemoresponsive, were correlated with at least one chemoresponsive BötC-VRG neuron. Correlation linkage maps and spike-triggered averages of phrenic nerve signals suggest transmission of chemoreceptor drive via a multipath network architecture: RTN-pF modulation of pre-BötC-VRG rostral-to-caudal excitatory inspiratory neuron chains is tuned by feedforward and recurrent inhibition from other inspiratory neurons and from "tonic" expiratory neurons.

Decorrelation of cortical inputs and motoneuron output

Oscillations in the primary motor cortex are transmitted through the corticospinal tract to the motoneuron pool. This pathway is believed to produce an effective and direct command from the motor cortex to the spinal motoneurons for the modulation of the force output. In this study, we used a computational model of a population of motoneurons to investigate the factors that can influence the transmission of the cortical input to the output of motoneurons, since it can be quantified by coherence analysis. The simulations demonstrated that, despite the nonlinearity of the motoneurons, oscillations present in the cortical input are transmitted to the output of the motoneuron pool at the same frequency. However, the interference introduced by the nonlinearity of the system increases the variability of the oscillations in output, introducing spectral lines whose frequency depends on the input frequencies and the motoneuron discharge rates. Moreover, an additional source of synaptic input common to all motoneurons but independent from the corticospinal component decorrelates the cortical input and motoneuron output and, thus, decreases the magnitude of the estimated coherence, even if the effective cortical drive does not change. These results indicate that the corticospinal input can effectively be sampled by a small population of motoneurons. However, the transmission of a corticospinal drive to the motoneuron pool is influenced by the nonlinearity of the spiking processes of the active motoneurons and by synaptic inputs common to the motoneuron population but independent from the cortical input.

Multiple extra-synaptic spillover mechanisms regulate prolonged activity in cerebellar Golgi cell-granule cell loops

Despite a wealth of in vitro and modelling studies it remains unclear how neuronal populations in the cerebellum interact in vivo. We address the issue of how the cerebellar input layer processes sensory information, with particular focus on the granule cells (input relays) and their counterpart inhibitory interneurones, Golgi cells. Based on the textbook view, granule cells excite Golgi cells via glutamate forming a negative feedback loop. However, Golgi cells express inhibitory mGluR2 receptors suggesting an inhibitory role for glutamate. We set out to test this glutamatergic paradox in Golgi cells. Here we show that granule cells and Golgi cells interact through extra-synaptic signalling mechanisms during sensory information processing, as well as synaptic mechanisms. We demonstrate that such interactions depend on granule cell-derived glutamate acting via inhibitory mGluR2 receptors leading causally to the suppression of Golgi cell activity for several hundreds of milliseconds. We further show that granule cell-derived inhibition of Golgi cell activity is regulated by GABA-dependent extra-synaptic Golgi cell inhibition of granule cells, identifying a regulatory loop in which glutamate and GABA may be critical regulators of Golgi cell–granule cell functional activity. Thus, granule cells may promote their own prolonged activity via paradoxical feed-forward inhibition of Golgi cells, thereby enabling information processing over long timescales.

Ventrolateral medullary functional connectivity and the respiratory and central chemoreceptor-evoked modulation of retrotrapezoid-parafacial neurons

The medullary ventral respiratory column (VRC) of neurons is essential for respiratory motor pattern generation; however, the functional connections among these cells are not well understood. A rostral extension of the VRC, including the retrotrapezoid nucleus/parafacial region (RTN-pF), contains neurons responsive to local perturbations of CO(2)/pH. We addressed the hypothesis that both local RTN-pF interactions and functional connections from more caudal VRC compartments--extending from the Bötzinger and pre-Bötzinger complexes to the ventral respiratory group (Böt-VRG)--influence the respiratory modulation of RTN-pF neurons and their responses to central chemoreceptor and baroreflex activation. Spike trains from 294 RTN-pF and 490 Böt-VRG neurons were monitored with multielectrode arrays along with phrenic nerve activity in 14 decerebrate, vagotomized cats. Overall, 214 RTN-pF and 398 Böt-VRG neurons were respiratory modulated; 124 and 95, respectively, were cardiac modulated. Subsets of these neurons were tested with sequential, selective, transient stimulation of central chemoreceptors and arterial baroreceptors; each cell's response was evaluated and categorized according to the change in firing rate (if any) following the stimulus. Cross-correlation analysis was applied to 2,884 RTN-pF↔RTN-pF and 8,490 Böt-VRG↔RTN-pF neuron pairs. In total, 174 RTN-pF neurons (59.5%) had significant features in short-time scale correlations with other RTN-pF neurons. Of these, 49 neurons triggered cross-correlograms with offset peaks or troughs (n = 99) indicative of paucisynaptic excitation or inhibition of the target. Forty-nine Böt-VRG neurons (10.0%) were triggers in 74 Böt-VRG→RTN-pF correlograms with offset features, suggesting that Böt-VRG trigger neurons influence RTN-pF target neurons. The results support the hypothesis that local RTN-pF neuron interactions and inputs from Böt-VRG neurons jointly contribute to respiratory modulation of RTN-pF neuronal discharge patterns and promotion or limitation of their responses to central chemoreceptor and baroreceptor stimulation.

Cerebellar Golgi cells in the rat receive convergent peripheral inputs via a lateral reticular nucleus relay

Golgi cells are important players in the function of the cerebellar cortex, controlling the flow of incoming information from mossy fibres to the granule cells, which excite other cortical neurons. We recently showed that in anaesthetized rats most Golgi cells respond to stimulation of afferents from a very wide peripheral receptive field with a long-lasting depression of firing. These responses are mediated via a crossed ascending afferent pathway but the supraspinal part of this pathway is unknown. Here we have examined the hypothesis that the lateral reticular nucleus, a brainstem nucleus with known broad afferent convergence that projects mossy fibres to much of the cerebellum, is involved. First, we showed that single-pulse electrical microstimulation within the lateral reticular nucleus can elicit long-lasting depressions in Golgi cells, which are qualitatively similar to those evoked by peripheral afferent stimulation. Second, we showed that the amplitude of the depressions of Golgi cell firing evoked by peripheral stimulation can be reduced by pharmacological manipulation of the lateral reticular nucleus, either ipsilateral or contralateral to the stimulus site, with local injections of either the GABA(A) receptor agonist muscimol or the AMPA receptor blocker 6-cyano-7-nitroquinoxaline-2,3-dione. This evidence suggests that the lateral reticular nucleus is a relay nucleus in the brainstem for peripheral afferent information in a pathway that generates Golgi cell long-lasting depression responses.

Fluctuations in isometric muscle force can be described by one linear projection of low-frequency components of motor unit discharge rates

The aim of the study was to investigate the relation between linear transformations of motor unit discharge rates and muscle force. Intramuscular (wire electrodes) and high-density surface EMG (13 x 5 electrode grid) were recorded from the abductor digiti minimi muscle of eight healthy men during 60 s contractions at 5%, 7.5% and 10% of the maximal force. Spike trains of a total of 222 motor units were identified from the EMG recordings with decomposition algorithms. Principal component analysis of the smoothed motor unit discharge rates indicated that one component (first common component, FCC) described 44.2 +/- 7.5% of the total variability of the smoothed discharge rates when computed over the entire contraction interval and 64.3 +/- 10.2% of the variability when computed over 5 s intervals. When the FCC was computed from four or more motor units per contraction, it correlated with the force produced by the muscle (62.7 +/- 10.1%) by a greater degree (P < 0.001) than the smoothed discharge rates of individual motor units (41.4 +/- 7.8%). The correlation between FCC and the force signal increased up to 71.8 +/- 13.1% when the duration and the shape of the smoothing window for discharge rates were similar to the average motor unit twitch force. Moreover, the coefficients of variation (CoV) for the force and for the FCC signal were correlated in all subjects (R(2) range = 0.14-0.56; P < 0.05) whereas the CoV for force was correlated to the interspike interval variability in only one subject (R(2) = 0.12; P < 0.05). Similar results were further obtained from measures on the tibialis anterior muscle of an additional eight subjects during contractions at forces up to 20% of the maximal force (e.g. FCC explained 59.8 +/- 11.0% of variability of the smoothed discharge rates). In conclusion, one signal captures most of the underlying variability of the low-frequency components of motor unit discharge rates and explains large part of the fluctuations in the motor output during isometric contractions.

Cardiorespiratory response to spontaneous cortical arousals during stage 2 and rapid eye movement sleep in healthy children

Arousal from sleep is associated with transient and abrupt cardiorespiratory changes, and elevated arousals associated with sleep disorders may trigger adverse cardiovascular sequela. In this paper, we provide the first data in children on cardiorespiratory responses to cortical arousal. Heart rate and ventilatory responses to arousal from stage 2 and rapid eye movement (REM) sleep were investigated in 40 normal, healthy Caucasian children (age: 7.7 +/- 2.6 years; body mass index z-score: 0.30 +/- 0.8). All children underwent overnight polysomnography studies. Cortical arousals were scored according to standard criteria. Heart rate changes were assessed over 30 s, starting 15 s prior to cortical arousal onset. Breathing rates were quantified three breaths before and after arousal onset. Arousals from stage 2 as well as REM sleep resulted in an R-R interval shortening of about 15%, independent of age and gender. The R-R interval shortening initiated at least 3 s before the cortical arousal onset. The breathing interval immediately after cortical arousal onset was significantly shortened (P < 0.001). In conclusion, cortical arousals in children are associated with an increase in breathing rate and significant heart rate accelerations, which typically precede the cortical arousal onset.

Arousal in obstructive sleep apnoea patients is associated with ECG RR and QT interval shortening and PR interval lengthening

Sudden cardiac death appears to be more prevalent during the normal sleeping hours in obstructive sleep apnoea (OSA) patients compared with the general population as well as to cardiovascular disease patients. The reasons for this remain unclear, but there are three likely main contributors to nocturnal death in OSA patients; cardiac arrhythmias, stroke/ruptured cerebral aneurism and myocardial infarction. Particularly marked cardiovascular system activation with arousal may play a role in initiating sudden adverse cardiovascular events in OSA. The purpose of this study was to investigate cardiac RR, QT and PR interval changes in the electrocardiogram (ECG) associated with spontaneous and respiratory-related arousals in OSA patients. A detailed observational study of ECG records obtained during conventional diagnostic sleep study with no further interventions was carried out in 20 patients (12 males, age 42.8 +/- 2.1 years, body mass index 35.1 +/- 1.9 kg m(-2), and respiratory disturbance index 51.8 +/- 6.4 events/hour). RR and QT intervals showed significant shortening during arousals. RR interval shortening was found to be greater during respiratory arousals when compared to spontaneous arousals. PR interval showed a trend toward a greater prolongation during respiratory arousal. QT interval shortening was weakly correlated with arterial oxygen saturation levels preceding arousal. In conclusion, these data suggest that despite greater cardiac acceleration following respiratory versus spontaneous arousals from sleep, QT shortening and PR prolongation responses are similar independent of arousal type. These data support that arousals produce quite marked and differential cardiac conduction system activation in OSA and that the degree and pattern of activation may be partly influenced by the presence and severity of preceding respiratory events.

Comparison between the degree of motor unit short-term synchronization and recurrence quantification analysis of the surface EMG in two human muscles

OBJECTIVE

To verify if non-linear recurrence analysis of the surface EMG is a suitable tool for assessing motor unit short-term synchronization.

METHODS

Surface and intramuscular EMG signals were recorded from the abductor digiti minimi and vastus medialis muscles of 12 and 10 healthy men, respectively, during isometric contractions. In the abductor digiti minimi, EMG signals were additionally recorded after a contraction sustained for 1min at 50% of the maximal force. In both muscles, percent of determinism (%DET) was estimated from the surface EMG and common input strength (CIS) index was computed from motor unit recordings.

RESULTS

For both muscles, CIS did not correlate with %DET (abductor digiti minimi: R(2)=0.11, P=0.12; vastus medialis: R(2)=0.04, P=0.56). Although the values of CIS for the vastus medialis were lower than those of the abductor digiti minimi (P<0.001), the %DET values did not differ between the two muscles (71.6+/-5.5% vs 66.9+/-8.7%; P=0.12).

CONCLUSION

The variable %DET extracted from the surface EMG is a poor indicator of the degree of motor unit short-term synchronization.

SIGNIFICANCE

The study provides a systematic evaluation of a technique previously proposed for the estimation of a clinically relevant characteristic of motor unit behavior.

Central and peripheral chemoreceptors evoke distinct responses in simultaneously recorded neurons of the raphé-pontomedullary respiratory network

The brainstem network for generating and modulating the respiratory motor pattern includes neurons of the medullary ventrolateral respiratory column (VRC), dorsolateral pons (PRG) and raphé nuclei. Midline raphé neurons are proposed to be elements of a distributed brainstem system of central chemoreceptors, as well as modulators of central chemoreceptors at other sites, including the retrotrapezoid nucleus. Stimulation of the raphé system or peripheral chemoreceptors can induce a long-term facilitation of phrenic nerve activity; central chemoreceptor stimulation does not. The network mechanisms through which each class of chemoreceptor differentially influences breathing are poorly understood. Microelectrode arrays were used to monitor sets of spike trains from 114 PRG, 198 VRC and 166 midline neurons in six decerebrate vagotomized cats; 356 were recorded during sequential stimulation of both receptor classes via brief CO(2)-saturated saline injections in vertebral (central) and carotid arteries (peripheral). Seventy neurons responded to both stimuli. More neurons were responsive only to peripheral challenges than those responsive only to central chemoreceptor stimulation (PRG, 20 : 4; VRC, 41 : 10; midline, 25 : 13). Of 16 474 pairs of neurons evaluated for short-time scale correlations, similar percentages of reference neurons in each brain region had correlation features indicative of a specific interaction with at least one target neuron: PRG (59.6%), VRC (51.0%) and raphé nuclei (45.8%). The results suggest a brainstem network architecture with connectivity that shapes the respiratory motor pattern via overlapping circuits that modulate central and peripheral chemoreceptor-mediated influences on breathing.

Functional connectivity in the pontomedullary respiratory network

Current models propose that a neuronal network in the ventrolateral medulla generates the basic respiratory rhythm and that this ventrolateral respiratory column (VRC) is profoundly influenced by the neurons of the pontine respiratory group (PRG). However, functional connectivity among PRG and VRC neurons is poorly understood. This study addressed four model-based hypotheses: 1) the respiratory modulation of PRG neuron populations reflects paucisynaptic actions of multiple VRC populations; 2) functional connections among PRG neurons shape and coordinate their respiratory-modulated activities; 3) the PRG acts on multiple VRC populations, contributing to phase-switching; and 4) neurons with no respiratory modulation located in close proximity to the VRC and PRG have widely distributed actions on respiratory-modulated cells. Two arrays of microelectrodes with individual depth adjustment were used to record sets of spike trains from a total of 145 PRG and 282 VRC neurons in 10 decerebrate, vagotomized, neuromuscularly blocked, ventilated cats. Data were evaluated for respiratory modulation with respect to efferent phrenic motoneuron activity and short-timescale correlations indicative of paucisynaptic functional connectivity using cross-correlation analysis and the "gravity" method. Correlogram features were found for 109 (3%) of the 3,218 pairs composed of a PRG and a VRC neuron, 126 (12%) of the 1,043 PRG-PRG pairs, and 319 (7%) of the 4,340 VRC-VRC neuron pairs evaluated. Correlation linkage maps generated for the data support our four motivating hypotheses and suggest network mechanisms for proposed modulatory functions of the PRG.

Climbing fibre-dependent changes in Golgi cell responses to peripheral stimulation

Golgi cells are important elements of the cerebellar cortex, controlling the flow of mossy fibre information to other cells via granule cells. Several anatomical reports suggest that climbing fibre afferents contact Golgi cells, and electrophysiological studies suggest that they depress Golgi cell firing. We reinvestigated this issue and, given that climbing fibres mediate synaptic plasticity in the cerebellar cortex, we have examined the effects of conjunctive stimulation of peripheral afferents and climbing fibres on Golgi cell responses. The results confirm that climbing fibre stimulation depresses Golgi cell firing at short latency. Golgi cells responded to stimulation of peripheral afferents with longer latency depressions of firing and after conjunctive stimulation with climbing fibres these were significantly reduced. The reductions developed progressively over 20 min of conjunctive stimulation and were persistent (up to 84 min). Temporal conjunction of the inputs was important because non-synchronous stimulation of climbing fibres and peripheral afferents failed to alter the peripheral afferent-evoked response in Golgi cells. In control experiments using either the same climbing fibre stimulation alone, or peripheral afferent stimulation paired with brainstem stimulation that did not activate climbing fibres, responses were not depressed. The results thus show that conjunctive stimulation of climbing fibres with other inputs to Golgi cells can induce long-term changes in Golgi cell responses in vivo. This raises the possibility that changes in Golgi cell peripheral responses mediated by climbing fibres can potentially contribute to cerebellar motor learning.

Motor unit synchronization in FDI and biceps brachii muscles of strength-trained males

Motor unit (MU) synchronization is the simultaneous or near-simultaneous firing of two MUs which occurs more often than would be expected by chance. The present study sought to investigate the effects of exercise training, muscle group, and force level, by comparing the magnitude of synchronization in the biceps brachii (BB) and first dorsal interosseous (FDI) muscles of untrained and strength-trained college-aged males at two force levels, 30% of maximal voluntary contraction (MVC) and 80% MVC. MU action potentials were recorded directly via an intramuscular needle electrode. The magnitude of synchronization was assessed using previously-reported synchronization indices: k', E, and CIS. Synchronization was significantly higher in the FDI than in the BB. Greater synchronization was observed in the strength-trained group with CIS, but not with E or k'. Also, synchronization was significantly greater at 80% MVC than at 30% MVC with E, but only moderately greater with CIS and there was no force difference with k'. Synchronization prevalence was found to be greater in the BB (80.1%) than in the FDI (71.5%). Thus, although the evidence is a bit equivocal, it appears that MU synchronization is greater at higher forces, and greater in strength-trained individuals than in untrained subjects.

Amplitude cancellation of motor-unit action potentials in the surface electromyogram can be estimated with spike-triggered averaging

The study presents analytical, simulation, and experimental analyses of amplitude cancellation of motor-unit action potentials (APs) in the interference electromyogram (EMG) and its relation to the size of the spike-triggered average (STA) EMG. The amount of cancellation of motor-unit APs decreases monotonically as a function of the ratio between the root mean square (RMS) of the motor-unit AP and the RMS of the interference EMG signal. The theoretical derivation of this association indicates a method to measure cancellation in individual motor units by STA of the interference and squared EMGs. The theoretical relation was examined in both simulated EMG signals generated by populations of 200 motor units and experimental recordings of 492 and 184 motor-unit APs in the vastus medialis and abductor digiti minimi muscles, respectively. Although the theoretical relation predicted (R2 = 0.95; P < 0.001) the amount of cancellation in the simulated EMGs, the presence of motor-unit synchronization decreased the strength of the association for small APs. The decrease in size of the STA obtained from the squared EMG relative to that extracted from the interference EMG was predicted by the experimental measure of cancellation (R2 = 0.65; P < 0.001, for vastus medialis; R2 = 0.26; P < 0.05, for abductor digiti minimi). The results indicate that cancellation of APs in the interference EMG can be analytically predicted and experimentally measured with STA from the discharge times of the motor units into the surface EMG.

Distribution and properties of visceral nociceptive neurons in rabbit cingulate cortex

Human imaging localizes most visceral nociceptive responses to anterior cingulate cortex (ACC), however, imaging in conscious subjects cannot completely control anticipatory and reflexive activity or resolve neuron activity. This study overcame these shortcomings by recording individual neuron responses in 12 anesthetized and paralyzed rabbits to define the visceronociceptive response pattern by region and layer. Balloon distension was applied to the colon at innocuous (15 mmHg) or noxious (60 mmHg) intensities, and innocuous and noxious mechanical, thermal and electrical stimuli were applied to the skin. Simultaneous recording from multiple regions assured differences were not due to anesthesia and neuron responses were resolved by spike sorting using principal components analysis. Of the total 346 neurons, 48% were nociceptive; responding to noxious levels of visceral or cutaneous stimulation, or both. Visceronociceptive neurons were most frequent in ACC (39%) and midcingulate cortex (MCC, 36%) and infrequent in retrosplenial cortex (RSC, 12%). In contrast, cutaneous nociceptive units were higher in MCC (MCC, 43%; ACC, 32%; RSC, 23%). Visceral-specific neurons were proportionately more frequent in ACC (37%), while cutaneous-specific units predominated in RSC (62.5%). Visceral nociceptive response durations were longer than those for cutaneous responses. Postmortem analysis of electrode tracks confirmed regional designations, and laminar analysis found inhibitory responses mainly in superficial layers and excitatory in deep layers. Thus, cingulate visceral nociception extends beyond ACC, this is the first report of nociceptive activity in RSC including nociceptive cutaneous responses, and these regional differences require a new model of cingulate nociceptive processing.

Lateral habenula stimulation inhibits rat midbrain dopamine neurons through a GABA(A) receptor-mediated mechanism

Transient changes in the activity of midbrain dopamine neurons encode an error signal that contributes to associative learning. Although considerable attention has been devoted to the mechanisms contributing to phasic increases in dopamine activity, less is known about the origin of the transient cessation in firing accompanying the unexpected loss of a predicted reward. Recent studies suggesting that the lateral habenula (LHb) may contribute to this type of signaling in humans prompted us to evaluate the effects of LHb stimulation on the activity of dopamine and non-dopamine neurons of the anesthetized rat. Single-pulse stimulation of the LHb (0.5 mA, 100 micros) transiently suppressed the activity of 97% of the dopamine neurons recorded in the substantia nigra and ventral tegmental area. The duration of the cessation averaged approximately 85 ms and did not differ between the two regions. Identical stimuli transiently excited 52% of the non-dopamine neurons in the ventral midbrain. Electrolytic lesions of the fasciculus retroflexus blocked the effects of LHb stimulation on dopamine neurons. Local application of bicuculline but not the SK-channel blocker apamin attenuated the effects of LHb stimulation on dopamine cells, indicating that the response is mediated by GABA(A) receptors. These data suggest that LHb-induced suppression of dopamine cell activity is mediated indirectly by orthodromic activation of putative GABAergic neurons in the ventral midbrain. The habenulomesencephalic pathway, which is capable of transiently suppressing the activity of dopamine neurons at a population level, may represent an important component of the circuitry involved in encoding reward expectancy.

Differential control of the scapulothoracic muscles in humans

The control of the scapulothoracic muscles trapezius (Tr) and serratus anterior (SA) has been examined in normal human subjects. Electromyographic recordings were made from the SA and Tr muscles (upper trapezius UTr, lower trapezius LTr) using surface electrodes placed bilaterally. Magnetic stimulation of the motor cortex and electrical stimulation of peripheral nerves were used to examine their descending and reflex control. The average optimal site of cortical stimulation was found to be the same for SA, UTr and LTr (an approximate centre of gravity of -0.6 cm, 3.7 cm where the centre of gravity is expressed as the mean anterio-posterior position, the mean medio-lateral position). Some asymmetry in the cortical representation of UTr was found in each individual tested. Magnetic stimulation evoked bilateral MEPs in Tr (latency contralateral (c) UTr 8.5 +/- 1.6 ms, ipsilateral (i) UTr 19.0 +/- 2.7 ms) but only contralateral responses were evoked in SA (11.2 +/- 2.6 ms). Electrical stimulation of the long thoracic nerve at two sites was used to examine homonymous and heteronymous reflexes of SA, while electrical stimulation of cervical nerve of C3/4 was used to examine the heteronymous reflexes of Tr. Ipsilateral SA H reflexes were evoked at a latency of 9.9 +/- 0.8 ms (proximal site) and 10.8 +/- 1.2 ms (distal site). No group I reflexes were evoked from SA to its contralateral homologue. No group I reflexes were evoked between Tr and SA. Finally, cross-correlation of activity from the Tr muscle pairs and the SA muscle pair revealed that the motoneurones of the Tr muscles share some common presynaptic input whereas there was no detectable common presynaptic input to the SA muscle pair. This study extends and consolidates knowledge regarding the neural control of trapezius and for the first time explores the neural control of SA. The study demonstrates a contrasting bilateral control of Tr and SA. These patterns of connections are discussed in relation to the contrasting bilateral functional roles of these muscles.

The amplitude of lower leg motor evoked potentials is a reliable measure when controlled for torque and motor task

OBJECTIVES

Motor evoked potential (MEP) amplitudes have the disadvantage of a high variability when repeatedly assessed. This affects the reliability of MEP amplitude measurements taken during the course of motor incomplete spinal cord injury (iSCI). The study investigated the reliability of anterior tibial (TA) MEP measures controlled for dorsal flexion torque and motor task.

METHODS

TA MEPs were recorded at 10, 20, 40 and 60% of maximal voluntary contraction (MVC) during a static and dynamic (isometric increase of dorsal flexion torque) motor task. To determine reliability, 20 healthy and five chronic iSCI subjects were tested twice (> or =7 days) by the same investigator. Intraclass correlation coefficients (ICCs) were calculated. MEP amplitudes and latencies were compared between 20 healthy and 29 iSCI subjects.

RESULTS

The reliability of MEP amplitude was in general good (ICC > or = 0.52) and was highest during the static task at 40% MVC (ICC = 0.77). The increased facilitation by the dynamic motor task showed the best reliability at 20% MVC (ICC = 0.48). The reliability was good to excellent for MEP latency (0.46 < or = ICC < or = 0.81), MVC (ICC > or = 0.90) and for the TMS threshold required to evoke a MEP response (ICC > or = 0.77). The torque generated by the MEP response ()0.02 < or = ICC < or = 0.55) and the duration of the silent period (0.07 < or = ICC < or = 0.50) were not reliable. Both MEP amplitudes and latencies differed significantly between healthy and iSCI subjects.

CONCLUSIONS

Controlling for torque generation and motor task establishes a reliability of TA MEP amplitudes that is sufficient for longitudinal assessments in motor incomplete SCI.

Cardiac changes during arousals from non-REM sleep in healthy volunteers

Our aim was to evaluate cardiac changes evoked by spontaneous and sound-induced arousals from sleep. Cardiac responses to spontaneous and auditory-induced arousals were recorded during overnight sleep studies in 28 young healthy subjects (14 males, 14 females) during non-rapid eye movement sleep. Computerized analysis was applied to assess beat-to-beat changes in heart rate, atrio-ventricular conductance, and ventricular repolarization from 30 s before to 60 s after the auditory tone. During both types of arousals, the most consistent change was the increase in the heart rate (in 62% of spontaneous and in 89% of sound-induced arousals). This was accompanied by an increase or no change in PR interval and by a decrease or no change in QT interval. The magnitude of all cardiac changes was significantly higher for tone-induced vs. spontaneous arousals (mean +/- SD for heart rate: +9 +/- 8 vs. +13 +/- 9 beats per min; for PR prolongation: 14 +/- 16 vs. 24 +/- 22 ms; for QT shortening: -12 +/- 6 vs. -20 +/- 9 ms). The prevalence of transient tachycardia and PR prolongation was also significantly higher for tone-induced vs. spontaneous arousals (tachycardia: 85% vs. 57% of arousals, P < 0.001; PR prolongation: 51% vs. 25% of arousals, P < 0.001). All cardiac responses were short-lasting (10-15 s). We conclude that cardiac pacemaker region, conducting system, and ventricular myocardium may be under independent neural control. Prolongation of atrio-ventricular delay may serve to increase ventricular filling during arousal from sleep. Whether prolonged atrio-ventricular conductance associated with increased sympathetic outflow to the ventricular myocardium contributes to arrhythmogenesis during sudden arousal from sleep remains to be evaluated.

Different responses of rat cerebellar Purkinje cells and Golgi cells evoked by widespread convergent sensory inputs

While the synaptic properties of Golgi cell-mediated inhibition of granule cells are well studied, less is known of the afferent inputs to Golgi cells so their role in information processing remains unclear. We investigated the responses of cerebellar cortical Golgi cells and Purkinje cells in Crus I and II of the posterior lobe cerebellar hemisphere to activation of peripheral afferents in vivo, using anaesthetized rats. Recordings were made from 70 Golgi cells and 76 Purkinje cells. Purkinje cells were identified by the presence of climbing fibre responses. Golgi cells were identified by both spontaneous firing pattern and response properties, and identification was confirmed using juxtacellular labelling of single neurones (n = 16). Purkinje cells in Crus II showed continuous firing at relatively high rates (25-60 Hz) and stimulation of peripheral afferents rarely evoked substantial responses. The most common response was a modest, long-latency, long-lasting increase in simple spike output. By comparison, the most common response evoked in Golgi cells by the same stimuli was a long-latency, long-lasting depression of firing, found in approximately 70% of the Golgi cells tested. The onsets of Golgi cell depressions had shorter latencies than the Purkinje cell excitations. Brief, short-latency excitations and reductions in firing were also evoked in some Golgi cells, and rarely in Purkinje cells, but in most cases long-lasting depressions were the only significant change in spike firing. Golgi cell responses could be evoked using air puff or tactile stimuli and under four different anaesthetic regimens. Long-lasting responses in both neurone types could be evoked from wide receptive fields, in many cases including distal afferents from all four limbs, as well as from trigeminal afferents. These Golgi cell responses are not consistent with the conventional feedback inhibition or 'gain control' models of Golgi cell function. They suggest instead that cerebellar cortical activity can be powerfully modulated by the general level of peripheral afferent activation from much of the body. On this basis, Golgi cells may act as a context-specific gate on transmission through the mossy fibre-granule cell pathway.

Detecting neuronal activity changes using an interspike interval algorithm compared with using visual inspection

A universally accepted method for efficiently detecting neuronal activity changes (NACs) in neurophysiological studies has not been established. Visual inspection is still considered to be one of the most reliable methods, although it is limited when it is used for analyzing large quantities of data. In this study, an algorithm that considers interspike intervals (ISIs) was developed to define the onset of NACs. Two criteria, involving the mean and the standard deviation (S.D.) of the ISIs during a control period, were used in the ISI algorithm to evaluate the NACs that occurred during a detection period. The first, an ISI decrease of more than 1 S.D. from the mean ISI of the control period, proved to be an effective criterion for qualifying the increased NACs (firing rate increases). The second, an ISI increase greater than 3 S.D.s, efficiently demarcated periods of decreased NACs (firing rate decreases). Statistically significant correlations between the detection of NAC onset times by the ISI algorithm and the detection of those times by visual inspections were observed after offline analyses of recorded neuronal activity. The present results suggest that this ISI algorithm is a reliable and efficient way of defining the onset of NACs.

Estimation of cortical silent period following transcranial magnetic stimulation using a computerised cumulative sum method

The cortical silent period (CSP) following transcranial magnetic stimulation (TMS) of the motor cortex can be used to measure intra-cortical inhibition and changes in a number of important pathologies affecting the central nervous system. The main drawback of this technique has been the difficulty in accurately identifying the onset and offset of the cortical silent period leading to inter-observer variability. We developed an automated method based on the cumulative sum (Cusum) technique to improve the determination of the duration and area of the cortical silent period. This was compared with experienced raters and two other automated methods. We showed that the automated Cusum method reliably correlated with the experienced raters for both duration and area of CSP. Compared with the automated methods, the Cusum also showed the strongest correlation with the experienced raters. Our results show the Cusum method to be a simple, graphical and powerful method of detecting low-intensity CSP that can be easily automated using standard software.

Excitatory synaptic potentials in spastic human motoneurons have a short rise-time

This study assessed whether changes in size or time-course of excitatory postsynaptic potentials (EPSPs) in motoneurons innervating spastic muscle could induce a greater synaptic response, and thereby contribute to reflex hyperexcitability. We compared motor unit (MU) firing patterns elicited by tendon taps applied to both spastic and contralateral (nonspastic) biceps brachii muscle in hemiparetic stroke subjects. Based on recordings of 115 MUs, significantly shortened EPSP rise times were present on the spastic side, but with no significant differences in estimated EPSP amplitude. These changes may contribute to hyperexcitable reflex responses at short latency, but the EPSP amplitude changes appear insufficient to account for global differences in reflex excitability.

Motor unit synchronization between medial and lateral vasti muscles

OBJECTIVE

Accurate neuromuscular control of the patellofemoral joint is important in knee joint mechanics. Strategies to coordinate the vasti muscles, such as motor unit synchronization, may simplify control of patellar tracking. This study investigated motor unit synchronization between vastus medialis (VM) and lateralis (VL).

METHODS

Electromyographic (EMG) recordings of single motor unit action potentials (MUAPs) were made from VM and single- and multi-unit recordings were made from VL. Synchronization was quantified from peaks in the cross-correlogram generated from single MUAP pairs in VL and VM. The proportion of motor units in VM with synchronized firing in VL was also quantified from peaks in averages of multiunit VL EMG triggered from the VM MUAP.

RESULTS

A high degree of synchronization of motor unit firing between VM and VL was identified. Results were similar for cross-correlation ( approximately 45% of cases) and triggered averages (approximately 41% of cases).

CONCLUSIONS

The data suggest that synchronization between VM and VL is higher than expected. Agreement between traditional cross-correlation and triggered averaging methods suggest that this new technique may provide a more clinically viable method to quantify synchronization.

SIGNIFICANCE

High synchronization between VM and VL may provide a solution to simplify control of the mechanically unstable patellofemoral joint.

Role in the inspiratory off-switch of vagal inputs to rostral pontine inspiratory-modulated neurons

Neurons of the pontine respiratory group (PRG) in the region of the nucleus parabrachialis medialis and the Kolliker-Fuse nucleus are believed to play an important role in promoting the inspiratory (I) off-switch (IOS). In decerebrate gallamine-paralyzed cats ventilated with a cycle-triggered pump system (lung inflation during the neural I phase), we studied the effects of vagal (V) afferent inputs on firing of I-modulated neurons (the most numerous population in PRG) and on I duration. The predominant V effect on unit activity was inhibitory, as shown by two types of test: (a) withholding of inflation during an I phase, which produced increase of unit firing and of its respiratory modulation (58/66 neurons in 14 cats), indicating disinhibition due to removal of phasic V input; (b) delivery of afferent V stimulus trains during a no-inflation I phase, which produced decrease of unit firing and of its respiratory modulation (20 neurons). In addition, application of V input during the neural expiratory (E) phase, which lengthened E phase duration, produced no effect on the neurons' firing, suggesting that the inhibition during I was presynaptic in origin. The results may be interpreted by the hypothesis that the medullary late-I presumptive IOS neurons receive excitatory inputs from the PRG I-modulated neurons as well as from V afferents.. With relatively strong V input, this pontine excitatory output is reduced by inhibition, whereas with relatively weak V input that excitatory output is increased due to reduction of inhibition. Thus the pontine and the V influences on the IOS can operate in a complementary manner dependent on state.

Integration-plot test for peri-stimulus time histograms in human motor units

A statistical test is proposed for peri-stimulus time histograms in human motor units for the case where a test stimulus is delivered at a constant interval after a previous discharge. This mathematically described test included the notion of the multiple comparison and thus achieved higher sensitivity than the previously proposed method. With regard to data acquisition, the interval from a sham stimulus to the next discharge was acquired as a control, and the total number of samples was set to be four times as large as that in the test situation to reduce the statistical scattering noise. A newly defined statistical object, the integration plot (timewise accumulation of the test histogram without control subtraction) was used for this statistical test. The integration plot had less noise than the cumulative sum (CUSUM) plot (1/square root(2) in theory) and thus represented the neural effect. To compare this integration-plot test with that of the CUSUM, a simulation experiment that compared two sample histograms (one of which had a faint structural change from 20 ms) was conducted. As a result, the present test succeeded in detecting the onset of the change point earlier (23 ms on average) than the CUSUM test (27 ms on average), and the detection probability was also higher (9 out of 10) than the CUSUM (6 out of 10). It was therefore confirmed experimentally that the present statistical test had higher sensitivity than that of the CUSUM proposed previously.

Formula-based approach of statistical tests for peri-stimulus time histograms

We introduce formula-based approaches for statistical tests regarding peri-stimulus time histograms (PSTHs) in cases when a conditioning stimulus was triggered by a motor-unit discharge with a constant delay time. The individual bin test and the cumulative sum test have recently been reported as methods for the quantitative analysis of PSTHs in the evaluation of human neural projections. These tests calculate confidence intervals using simulated point processes on the discharge of a motor unit, but require a significant amount of calculation time for the point process simulation. In order to overcome such a disadvantage in practical use, we provide a statistical formulization of the two above-mentioned tests using a combination theory. In general, the time required for calculating confidence intervals using these formula-based approaches was 2-13 times faster than when using the previous approaches. Unlike the previous ones, these formula-based approaches do not need to judge that a simulation process converges to the substationary state, and calculate an ideal distribution of statistical noise on the histogram, thereby providing high accuracy. We conclude that the formula-based approaches increase reliability and are sufficiently sophisticated for practical use.

Reassessment of parameters for applying motor-unit triggered stimuli in peri-stimulus time histograms

We reassessed the response properties of peri-stimulus time histograms (PSTHs) in cases when a test stimulus was triggered by a motor-unit discharge with a constant delay time. In this experiment, single motor unit action potentials were recorded from the right tibialis anterior (TA) muscle of five healthy persons. A test stimulus of the common peroneal nerve with low intensity to activate only Ia afferents of the TA was applied through a bipolar stimulating electrode placed distal to the neck of the fibula. We obtained several PSTHs with various delay times and stimulus intensities in the same recording session for maintaining the background property as the same among the test situations. As a result, we confirmed three characteristics of PSTHs from observed data: (1) given the same delay time (the same background firing properties), a weaker stimulus intensity evokes a lessened effect on PSTHs, naturally; (2) delay time alters the induction balance of direct and indirect effects on PSTHs even if the stimulus intensity is the same because the background firing properties are different; and (3) response probabilities do not correspond directly to stimulus intensities when background firing properties are different; it is possible for a relatively strong intensity stimulus to produce a weaker effect than a weak stimulus. We concluded that comparisons of effects taken at different phases in the control distribution (and also effects taken from different control distributions or different motor units) can be misleading. Therefore, such comparisons should only be made within data obtained from the same phase in the same control distribution of the same motor unit.

Reflex inhibition of human inspiratory muscles in response to contralateral phrenic nerve stimulation

In animals, high-intensity unilateral stimulation of the phrenic nerve results in short-latency inhibition of phrenic and intercostal nerve activity bilaterally. This study provides the first demonstration in human subjects of a short-latency inhibitory response in the contralateral scalene, parasternal intercostal and diaphragm muscles to single stimuli delivered at cervical level to the phrenic nerve. Electromyographic (EMG) responses were recorded with intramuscular and surface electrodes. An inhibitory response with an onset latency of approximately 35 ms followed by a long-latency excitatory response at approximately 100 ms were observed in the three inspiratory muscles. The inhibition was evident in single trials, averaged EMG, histograms of the discharge of single motor units, and even when the phrenic nerve stimulus intensity was relatively low. Thus, the inhibition may be mediated by large-diameter muscle afferents. The latency of this potent inhibitory response to contralateral phrenic nerve stimulation is too long to be mediated via a simple spinal circuit and may involve a brainstem projection.

Electrocardiographic changes associated with the nasopharyngeal reflex in conscious rabbits: vago-sympathetic co-activation

Electrocardiographic responses were assessed in conscious rabbits when the nasopharyngeal reflex was elicited by inhalation of formaldehyde vapour. There was a profound fall in heart rate (224+/-5 to 64+/-4 beats per min (bpm)) associated with abnormal or absent P-waves. There were no changes in the QRS complex. The R-T interval (control value 118+/-4 ms) was initially shortened to 107+/-3 ms and then prolonged to 130+/-4 ms. Heart rate and P-wave changes were prevented by muscarinic cholinergic blockade with methylscopolamine. The R-T shortening was reduced by 79+/-4% by beta-adrenergic blockade with propranolol. Methylscopolamine also unmasked small tachycardic responses (5-25 bpm) in 5/7 animals. This tachycardia was prevented by propranolol. Thus both parasympathetic vagal cardiac nerves and sympathetic cardiac nerves are activated during the nasopharyngeal reflex, with increased vagal effects in the sino-atrial node, and increased sympathetic effects in the ventricular myocardium.

Single-unit activity related to bimanual arm movements in the primary and supplementary motor cortices

Single units were recorded from the primary motor (MI) and supplementary motor (SMA) areas of Rhesus monkeys performing one-arm (unimanual) and two-arm (bimanual) proximal reaching tasks. During execution of the bimanual movements, the task related activity of about one-half the neurons in each area (MI: 129/232, SMA: 107/206) differed from the activity during similar displacements of one arm while the other was stationary. The bulk of this "bimanual-related" activity could not be explained by any linear combination of activities during unimanual reaching or by differences in kinematics or recorded EMG activity. The bimanual-related activity was relatively insensitive to trial-to-trial variations in muscular activity or arm kinematics. For example, trials where bimanual arm movements differed the most from their unimanual controls did not correspond to the ones where the largest bimanual neural effects were observed. Cortical localization established by using a mixture of surface landmarks, electromyographic recordings, microstimulation, and sensory testing suggests that the recorded neurons were not limited to areas specifically involved with postural muscles. By rejecting this range of alternative explanations, we conclude that neural activity in MI as well as SMA can reflect specialized cortical processing associated with bimanual movements.