Statistical signs of synaptic interaction in neurons

Visual perception of ambiguous figures: synchronization based neural models

We develop and study two neural network models of perceptual alternations. Both models have a star-like architecture of connections with a central element connected to a set of peripheral elements. A particular perception is simulated in terms of partial synchronization between the central element and some sub-group of peripheral elements. The first model is constructed from phase oscillators and the mechanism of perceptual alternations is based on chaotic intermittency under fixed parameter values. Similar to experimental evidence, the distribution of times between perceptual alternations is represented by the gamma distribution. The second model is built of spiking neurons of the Hodgkin-Huxley type. The mechanism of perceptual alternations is based on plasticity of inhibitory synapses which increases the inhibition from the central unit to the neural assembly representing the current percept. As a result another perception is formed. Simulations show that the second model is in good agreement with behavioural data on switching times between percepts of ambiguous figures and with experimental results on binocular rivalry of two and four percepts.

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

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

Organization of interneuronal connections in the nucleus accumbens in "impulsive" and "self-controlled" behavior in cats

In behavioral experiments, cats placed in a situation of choosing between a high-value time-delayed and a low-value rapid food reinforcement elected to wait for the preferred reward (they demonstrated "self-control") or to obtain the worse reward quickly (they demonstrated impulsive behavior). On the basis of the selected behavioral strategy, the cats were divided into three groups - "impulsive," "ambivalent," and "self-controlled." Cross-correlation analysis was used to assess the linked activity of cells in the nucleus accumbens, which reflects the nature of interactions between close-lying neurons. In cats with self-control, interneuronal interactions appeared in a significantly larger proportion of cases than in impulsive cats. In combinations resulting in long-latency reactions, cats with self-controlled and impulsive behavior showed no significant difference in the occurrence frequency of interneuronal interactions. The numbers of interneuronal interactions were greater during erroneous responses as compared with correctly performed reactions in animals of the different groups. These data indicate a key role for the interrelated activity of nucleus accumbens neurons in organizing the pattern of long-latency responses typical of self-controlled behavior.

Pontine-ventral respiratory column interactions through raphe circuits detected using multi-array spike train recordings

Recently, Segers et al. identified functional connectivity between the ventrolateral respiratory column (VRC) and the pontine respiratory group (PRG). The apparent sparseness of detected paucisynaptic interactions motivated consideration of other potential functional pathways between these two regions. We report here evidence for "indirect" serial functional linkages between the PRG and VRC via intermediary brain stem midline raphé neurons. Arrays of microelectrodes were used to record sets of spike trains from a total of 145 PRG, 282 VRC, and 340 midline neurons in 11 decerebrate, vagotomized, neuromuscularly blocked, ventilated cats. Spike trains of 13,843 pairs of neurons that included at least one raphé cell were screened for respiratory modulation and short-time scale correlations. Significant correlogram features were detected in 7.2% of raphé-raphé (291/4,021), 4.3% of VRC-raphé (292/6,755), and 4.0% of the PRG-raphé (124/3,067) neuron pairs. Central peaks indicative of shared influences were the most common feature in correlations between pairs of raphé neurons, whereas correlated raphé-PRG and raphé-VRC neuron pairs displayed predominantly offset peaks and troughs, features suggesting a paucisynaptic influence of one neuron on the other. Overall, offset correlogram features provided evidence for 33 VRC-to-raphé-to-PRG and 45 PRG-to-raphé-to-VRC correlational linkage chains with one or two intermediate raphé neurons. The results support a respiratory network architecture with parallel VRC-to-PRG and PRG-to-VRC links operating through intervening midline circuits, and suggest that raphé neurons contribute to the respiratory modulation of PRG neurons and shape the respiratory motor pattern through coordinated divergent actions on both the PRG and VRC.

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

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

Spatial and temporal scales of neuronal correlation in primary visual cortex

The spiking activity of cortical neurons is correlated. For instance, trial-to-trial fluctuations in response strength are shared between neurons, and spikes often occur synchronously. Understanding the properties and mechanisms that generate these forms of correlation is critical for determining their role in cortical processing. We therefore investigated the spatial extent and functional specificity of correlated spontaneous and evoked activity. Because feedforward, recurrent, and feedback pathways have distinct extents and specificity, we reasoned that these measurements could elucidate the contribution of each type of input. We recorded single unit activity with microelectrode arrays which allowed us to measure correlation in many hundreds of pairings, across a large range of spatial scales. Our data show that correlated evoked activity is generated by two mechanisms that link neurons with similar orientation preferences on different spatial scales: one with high temporal precision and a limited spatial extent (approximately 3 mm), and a second that gives rise to correlation on a slow time scale and extends as far as we were able to measure (10 mm). The former is consistent with common input provided by horizontal connections; the latter likely involves feedback from extrastriate cortex. Spontaneous activity was correlated over a similar spatial extent, but approximately twice as strongly as evoked activity. Visual stimuli thus caused a substantial decrease in correlation, particularly at response onset. These properties and the circuit mechanism they imply provide new constraints on the functional role that correlation may play in visual processing.

Bayesian inference of functional connectivity and network structure from spikes

Current multielectrode techniques enable the simultaneous recording of spikes from hundreds of neurons. To study neural plasticity and network structure it is desirable to infer the underlying functional connectivity between the recorded neurons. Functional connectivity is defined by a large number of parameters, which characterize how each neuron influences the other neurons. A Bayesian approach that combines information from the recorded spikes (likelihood) with prior beliefs about functional connectivity (prior) can improve inference of these parameters and reduce overfitting. Recent studies have used likelihood functions based on the statistics of point-processes and a prior that captures the sparseness of neural connections. Here we include a prior that captures the empirical finding that interactions tend to vary smoothly in time. We show that this method can successfully infer connectivity patterns in simulated data and apply the algorithm to spike data recorded from primary motor (M1) and premotor (PMd) cortices of a monkey. Finally, we present a new approach to studying structure in inferred connections based on a Bayesian clustering algorithm. Groups of neurons in M1 and PMd show common patterns of input and output that may correspond to functional assemblies.

Inferring functional connections between neurons

A central question in neuroscience is how interactions between neurons give rise to behavior. In many electrophysiological experiments, the activity of a set of neurons is recorded while sensory stimuli or movement tasks are varied. Tools that aim to reveal underlying interactions between neurons from such data can be extremely useful. Traditionally, neuroscientists have studied these interactions using purely descriptive statistics (cross-correlograms or joint peri-stimulus time histograms). However, the interpretation of such data is often difficult, particularly as the number of recorded neurons grows. Recent research suggests that model-based, maximum likelihood methods can improve these analyses. In addition to estimating neural interactions, application of these techniques has improved decoding of external variables, created novel interpretations of existing electrophysiological data, and may provide new insight into how the brain represents information.

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.

Brain dynamics for perception of tactile allodynia (touch-induced pain) in postherpetic neuralgia

Postherpetic neuralgia (PHN) is a debilitating chronic pain condition often accompanied by a sensation of pain when the affected region is touched (tactile allodynia). Here we identify brain regions involved in stimulus-induced touch-evoked pain (dynamical mechanical allodynia, DMA), compare brain activity between DMA and spontaneous pain (described earlier for the same patients in [Geha PY, Baliki MN, Chialvo DR, Harden RN, Paice JA, Apkarian AV. Brain activity for spontaneous pain of postherpetic neuralgia and its modulation by lidocaine patch therapy. Pain 2007;128:88-100]), delineate regions that specifically code the magnitude of perceived allodynia, and show the transformation of allodynia-related information in the brain as a time-evolving network. Eleven PHN patients were studied for DMA and its modulation with Lidoderm therapy (patches of 5% lidocaine applied to the PHN affected body part). Continuous ratings of pain while the affected body part was brushed during fMRI were contrasted with non-painful touch when brushing was applied to an equivalent opposite body site, and with fluctuations of a bar observed during scanning, at three sessions relative to Lidoderm treatment. Lidoderm treatment did not decrease DMA ratings but did decrease spontaneous pain. Multiple brain areas showed preferential activity for allodynia. However, mainly responses in the bilateral putamen and left medial temporal gyrus were related to the magnitude of allodynia. Both DMA and spontaneous pain perceptions were best represented within the same sub-cortical structures but with minimal overlap, implying that PHN pain modulates behavioral learning and hedonics. These results have important clinical implications regarding adequate therapy.

Interhemisphere asymmetry of the neocortex and hippocampus during orientational-investigative behavior and freezing in rabbits

Cross-correlation histograms were used to study the discharges of cortical neurons in symmetrical leads in both hemispheres (visual and parietal areas) and in hippocampal field CA1 on the right and left sides in rabbits in conditions of free behavior during exposure to emotionally significant stimuli. During active orientational-investigative responses to stimuli, as compared with baseline, the neocortex showed increases in left-sided influences on cells in the right hemisphere, with delays of up to 100 msec, which led to the appearance of asymmetry in the interhemisphere interaction, with left-sided dominance. On freezing, the left-sided influence weakened and the right hemisphere became dominant. Interhemisphere asymmetry in hippocampal neuron activity was seen, and was reciprocal to the asymmetry observed in the neocortex. The active investigative response increased right-sided influences in the hippocampus with delays of up to 200 msec, leading to right-sided dominance. Freezing was associated with increases in left-sided influences, such that the left side was dominant. The interaction of cells in the hippocampus was largely at the frequencies of the theta rhythm during active movement and in the delta range on freezing. These data lead to the conclusion that the active or passive nature of behavioral movement reactions to emotionally significant stimuli correlates with changes in the asymmetry of interhemisphere neuron interaction at the levels of the cerebral cortex and hippocampus.

Auditory nerve inputs to cochlear nucleus neurons studied with cross-correlation

The strength of synapses between auditory nerve (AN) fibers and ventral cochlear nucleus (VCN) neurons is an important factor in determining the nature of neural integration in VCN neurons of different response types. Synaptic strength was analyzed using cross-correlation of spike trains recorded simultaneously from an AN fiber and a VCN neuron in anesthetized cats. VCN neurons were classified as chopper, primarylike, and onset using previously defined criteria, although onset neurons usually were not analyzed because of their low discharge rates. The correlograms showed an excitatory peak (EP), consistent with monosynaptic excitation, in AN-VCN pairs with similar best frequencies (49% 24/49 of pairs with best frequencies within +/-5%). Chopper and primarylike neurons showed similar EPs, except that the primarylike neurons had shorter latencies and shorter-duration EPs. Large EPs consistent with end bulb terminals on spherical bushy cells were not observed, probably because of the low probability of recording from one. The small EPs observed in primarylike neurons, presumably spherical bushy cells, could be derived from small terminals that accompany end bulbs on these cells. EPs on chopper or primarylike-with-notch neurons were consistent with the smaller synaptic terminals on multipolar and globular bushy cells. Unexpectedly, EPs were observed only at sound levels within about 20 dB of threshold, showing that VCN responses to steady tones shift from a 1:1 relationship between AN and VCN spikes at low sound levels to a more autonomous mode of firing at high levels. In the high level mode, the pattern of output spikes seems to be determined by the properties of the postsynaptic spike generator rather than the input spike patterns. The EP amplitudes did not change significantly when the presynaptic spike was preceded by either a short or long interspike interval, suggesting that synaptic depression and facilitation have little effect under the conditions studied here.

Interactions between neurons in the frontal cortex and hippocampus in cats trained to select reinforcements of different value in conditions of cholinergic deficiency

An operant food-related conditioned reflex was developed in six cats by the "active choice" protocol: short-latency pedal presses were followed by presentation of low-quality reinforcement (bread-meat mix), while long-latency pedal presses were followed by presentation of high-quality reinforcement (meat). Animals differed in terms of their food-procuring strategies, displaying "self-control," "ambivalence," or "impulsivity." Multineuron activity was recorded from the frontal cortex and hippocampus (field CA3). Cross-correlation analysis of interneuronal interactions within (local networks) and between (distributed networks) study structures showed that the numbers of interneuronal interactions in both local and distributed networks were maximal in animals with "self-control." On the background of systemic administration of the muscarinic cholinoreceptor blockers scopolamine and trihexyphenidyl, the numbers of interneuronal interactions decreased, while "common source" influences increased. This correlated with impairment of the reproduction of the selected strategy, primarily affecting the animals' self-controlled behavior. These results show that the "self-control" strategy is determined by the organization of local and distributed networks in the frontal cortex and hippocampus.

Generation of synthetic spike trains with defined pairwise correlations

Recent technological advances as well as progress in theoretical understanding of neural systems have created a need for synthetic spike trains with controlled mean rate and pairwise cross-correlation. This report introduces and analyzes a novel algorithm for the generation of discretized spike trains with arbitrary mean rates and controlled cross correlation. Pairs of spike trains with any pairwise correlation can be generated, and higher-order correlations are compatible with common synaptic input. Relations between allowable mean rates and correlations within a population are discussed. The algorithm is highly efficient, its complexity increasing linearly with the number of spike trains generated and therefore inversely with the number of cross-correlated pairs.

A simple connectivity scheme for sparse coding in an olfactory system

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

Differential taste coding of salt and acid by correlative activities between taste-sensitive neuron types in rat gustatory cortex

Using a multi-electrode recording technique, the present study aimed to elucidate the role of broadly-tuned taste-sensitive neurons in the rat gustatory cortex in discriminating between salt and acid. A majority of taste-sensitive neurons (94/119 neurons; 78%) were classified as NaCl (N)-, HCl (H)- or NaCl and HCl (NH)-best neurons. Of 63 neuron pairs (94 neurons), 31 showed significant peaks and/or troughs in their cross-correlograms (CCs) during taste stimulation periods. During NaCl stimulation, the incidence of significant correlation and the mean frequency of correlated discharges (FC) in the N/N and NH/NH pairs were higher than those in the other best-taste pairs. In contrast, during HCl stimulation both indices in the N/N or H/H pairs were very low, while those in the NH/NH pairs were high. These results suggest that (1) correlated activities between N-best neurons and those between NH-best neurons play a significant role in taste quality coding of salt, and that (2) correlated activities between NH-best neurons may be important for sour taste coding as well. Peak formation in CCs tended to be more frequent in the homo-types (N/N, H/H and NH/NH pairs) than in the hetero-types (N/NH, N/H and H/NH pairs). In contrast, troughs were observed mostly in the hetero-types. Inhibitory interaction in hetero-type pairs together with coactivation in homo-type pairs may enhance taste discrimination by taste-sensitive neuron populations.

Organisation of common inputs to motoneuron pools of human masticatory muscles

OBJECTIVE

To determine the pattern of organization of common inputs to the motoneuron pools of individual muscles in the masticatory system.

METHODS

Six subjects bit on a rubber-coated wooden splint placed between the upper and lower incisor teeth. We recorded the surface electromyogram (EMG) of co-contracting masseter, temporalis and digastric muscles bilaterally during isometric jaw closing at 5%, 10%, 20% and 40% of maximal voluntary masseter EMG.

RESULTS

The cross-correlograms of the EMGs of homologous muscle pairs indicate that there are common synaptic inputs to the motoneuron pools of the left and right masseter, and left and right digastric muscles, but not to left and right temporalis. The amplitude of the central peak in masseter and digastric correlograms increased with bite force. When the activity of ipsilateral muscle pairs was cross-correlated, central peaks were prominent for masseter-digastric and masseter-temporalis muscle pairs, and the peak amplitudes increased significantly with bite force. In contrast, no significant central peak was observed for temporalis-digastric muscle pairs at any level of voluntary biting.

CONCLUSIONS

We conclude that there is synchronous modulation of input bilaterally to the masseter muscles and to the digastric muscles but not to the temporalis muscles. There is synchronous modulation of input to ipsilateral masseter-digastric and masseter-temporalis muscle pairs but not to temporalis and digastric muscles.

SIGNIFICANCE

The extent of common input to motoneuron pools of muscles acting around a common joint varies for different muscle pairs, and is not simply a function of whether the muscles of the pair are synergists or antagonists.

Synchronous Activity in Cat Visual Cortex Encodes Collinear and Cocircular Contours

We explored how contour information in primary visual cortex might be embedded in the simultaneous activity of multiple cells recorded with a 100-electrode array. Synchronous activity in cat visual cortex was more selective and predictable in discriminating between drifting grating and concentric ring stimuli than changes in firing rate. Synchrony was found even between cells with wholly different orientation preferences when their receptive fields were circularly aligned, and membership in synchronous groups was orientation and curvature dependent. The existence of synchrony between cocircular cells reinforces its role as a general mechanism for contour integration and shape detection as predicted by association field concepts. Our data suggest that cortical synchrony results from common and synchronous input from earlier visual areas and that it could serve to shape extrastriate response selectivity.

Linkage of Neuron Spike Activity in the Right and Left Amygdala in Food Motivation and Emotional Tension

Cross-correlation histograms were plotted to study the linkage of spike activity in simultaneously recorded neurons in the central nucleus of the right and left amygdala in rabbits in calm waking, after 24 hours of food deprivation, in satiation, and in emotional tension (on presentation of dogs). Histograms showed peaks displaced from zero in 50-67% of cases. In hunger, many cases (52%) showed pairs in which the first neuron to discharge was in the left amygdala, this being followed by a neuron in the right amygdala (peaks from 10-50 and 130-150 msec). Firing in the opposite order was seen more rarely (36%). On presentation of dogs, there was an increase in the number of cases showing a common input to neurons, along with an increase in the leading role of neurons in the right amygdala (57%), due to increases in inhibitory influences from this area on cells in the left amygdala. The interaction of amygdalar neurons in these states was in most cases at frequencies in the delta range (74%), mainly at 2-4 Hz.

Correlations Between Purkinje Cell Single-Unit Activity and Simultaneously Recorded Field Potentials in the Immediately Underlying Granule Cell Layer

Evidence from both anatomical and physiological studies suggests that the ascending segment of the granule cell axon provides a large, driving input to overlying Purkinje cells. In the current experiments, we used dual recording electrodes to simultaneously record spike activity of Purkinje cells and multiunit field potential activity in the directly underlying granule cell layer. These dual recordings were performed both during periods of spontaneous (“background”) firing and also after peripheral tactile stimulation. The results demonstrate that in the large majority of cases, there is a strong positive correlation between spontaneous Purkinje cell simple spikes and spontaneous activity in the immediately underlying granule cell layer. The strength of this correlation was dependent on both the firing rate of the Purkinje cell as well as on the rate of granule cell layer multiunit activity. In addition, for any given pair of recordings, the correlation seen during spontaneous activity accurately predicted the magnitude and time course of responses evoked by peripheral tactile stimulation. These results provide additional evidence that the synapses associated with the ascending segment of the granule cell axon have a substantial influence on Purkinje cell output. This relationship is considered in the context of our ongoing reevaluation of the physiological relationship between cerebellar granule and Purkinje cells.

Synchronized activity of renal neurons and their pattern of discharge in rabbits

Mechanisms of synchronization of renal neurons were studied by correlating its primary features (peak width and peak height of correlation histograms) with parameters derived from interspike-interval histograms of their resting activity. In anaesthetized rabbits the synchronous correlograms had a peak width of 210+/-8.8 ms (x+/-SEM; n=156) and their peak height was 4.3+/-2 spikes/s. Following parameters were calculated from interval-histograms of single renal units: the shortest and preferred interspike-intervals (their values were 9+/-0.5 ms and 11+/-1 ms, respectively), and longest interspike-interval and spread of histogram (amounting to 5.79+/-0.45 s and 5.76+/-0.37 s, respectively). Peak width of cross-correlogram was significantly correlated both with the longest interspike-interval (r=-0.426) and histogram spread (r=-0.431) while the strength of relationship with the shortest and preferred interspike-intervals was non-significant. On the other hand, peak height of correlograms was significantly correlated with the longest and preferred interspike-intervals (r=-0.179 and r=-0.191) and histogram spread (r=0.191). These data suggest various properties of both primary features of synchronized firing. They also support the concept of DiBona that various intrarenal effectors may be activated by single renal sympathetic neuron due to information encoded in its discharge pattern.

Analyzing functional connectivity using a network likelihood model of ensemble neural spiking activity

Analyzing the dependencies between spike trains is an important step in understanding how neurons work in concert to represent biological signals. Usually this is done for pairs of neurons at a time using correlation-based techniques. Chornoboy, Schramm, and Karr (1988) proposed maximum likelihood methods for the simultaneous analysis of multiple pair-wise interactions among an ensemble of neurons. One of these methods is an iterative, continuous-time estimation algorithm for a network likelihood model formulated in terms of multiplicative conditional intensity functions. We devised a discrete-time version of this algorithm that includes a new, efficient computational strategy, a principled method to compute starting values, and a principled stopping criterion. In an analysis of simulated neural spike trains from ensembles of interacting neurons, the algorithm recovered the correct connectivity matrices and interaction parameters. In the analysis of spike trains from an ensemble of rat hippocampal place cells, the algorithm identified a connectivity matrix and interaction parameters consistent with the pattern of conjoined firing predicted by the overlap of the neurons' spatial receptive fields. These results suggest that the network likelihood model can be an efficient tool for the analysis of ensemble spiking activity.

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

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

Stimulus dependence of neuronal correlation in primary visual cortex of the macaque

Nearby cortical neurons often have correlated trial-to-trial response variability, and a significant fraction of their spikes occur synchronously. These two forms of correlation are both believed to arise from common synaptic input, but the origin of this input is unclear. We investigated the source of correlated responsivity by recording from pairs of single neurons in primary visual cortex of anesthetized macaque monkeys and comparing correlated variability and synchrony for spontaneous activity and activity evoked by stimuli of different orientations and contrasts. These two stimulus manipulations would be expected to have different effects on the cortical pool providing input to the recorded pair: changing stimulus orientation should recruit different populations of cells, whereas changing stimulus contrast affects primarily the relative strength of sensory drive and ongoing cortical activity. Consistent with this predicted difference, we found that correlation was affected by these stimulus manipulations in different ways. Synchrony was significantly stronger for orientations that drove both neurons well than for those that did not, but correlation on longer time scales was orientation independent. Reducing stimulus contrast resulted in a decrease in the temporal precision of synchronous firing and an enhancement of correlated response variability on longer time scales. Our results thus suggest that correlated responsivity arises from mechanisms operating at two distinct timescales: one that is orientation tuned and that determines the strength of temporally precise synchrony, and a second that is contrast sensitive, of low temporal frequency, and present in ongoing cortical activity.

Periodic modulation of motor-unit activity in extrinsic hand muscles during multidigit grasping

We recently examined the extent to which motor units of digit flexor muscles receive common input during multidigit grasping. This task elicited moderate to strong motor-unit synchrony (common input strength, CIS) across muscles (flexor digitorum profundus, FDP, and flexor pollicis longus, FPL) and across FDP muscle compartments, although the strength of this common input was not uniform across digit pairs. To further characterize the neural mechanisms underlying the control of multidigit grasping, we analyzed the relationship between firing of single motor units from these hand muscles in the frequency domain by computing coherence. We report three primary findings. First, in contrast to what has been reported in intrinsic hand muscles, motor units belonging to different muscles and muscle compartments of extrinsic digit flexors exhibited significant coherence in the 0- to 5- and 5- to 10-Hz frequency ranges and much weaker coherence in the higher 10-20 Hz range (maximum 0.0025 and 0.0008, respectively, pooled across all FDP compartment pairs). Second, the strength and incidence of coherence differed considerably across digit pairs. Third, contrary to what has been reported in the literature, across-muscle coherence can be stronger and more prevalent than within-muscle coherence, as FPL-FDP2 (thumb-index digit pair) exhibited the strongest and most prevalent coherence in our data (0.010 and 43% at 3 Hz, respectively). The heterogeneous organization of common input to these muscles and muscle compartments is discussed in relation to the functional role of individual digit pairs in the coordination of multiple digit forces in grasping.

Pooled spike trains of correlated presynaptic inputs as realizations of cluster point processes

The pooled spike trains of correlated presynaptic terminals acting synchronously upon a single neuron are realizations of cluster point processes: the notions of spikes synchronizing in bursts and of points bunching in clusters are conceptually identical. The primary processes constituent specifies the timing of the cluster series; subsidiary processes and poolings specify burst structure and tightness. This representation and the Poisson process representation of independent terminals complete the formal approach to pooled trains. The notion's usefulness was illustrated by expressing physiological questions in terms of those constituents, each possessing a clear biological embodiment; constituents provided the control variables in simulations using leaky integrate-and-fire postsynaptic neurons excited by multiple weak terminals. Regular or irregular primary processes and bursts series determined low or high postsynaptic dispersions. When convergent set synchrony increased, its postsynaptic consequences approached those of single powerful synapses; concomitantly, output spike trains approached periodic, quasiperiodic, or aperiodic behaviors. The sequence in which terminals fired within bursts affected the predictee and predictor roles of presynaptic and postsynaptic spikes; when inhibition was added, EPSP and IPSP delays and order were influential (summation was noncommutative). Outputs to different correlations were heterogeneous; heterogeneity was accentuated by conditioning by variables such as DC biases.

Structure of dependent relationships between neurons in the sensorimotor cortex of the left and right hemispheres in rabbits in immobilization catatonia

Dependence in the activity of sensorimotor cortex neurons recorded simultaneously in the left and right hemispheres was detected in rabbits in baseline conditions, during the state of immobilization ("animal hypnosis"), and recovery of animals from this state. In baseline conditions, the total percentage of dependent relationships between close-lying (within 50 microm) neurons in the left hemisphere was significantly smaller than in the right hemisphere and did not change either in the state of immobilization or on recovery from it. The total percentage of dependent relationships between close-lying neurons in the right hemisphere decreased significantly during immobilization and returned to baseline levels on recovery from this state. The percentage of dependent relationships between distant (500 microm) neurons in immobilization, conversely, showed no change in the cortex of the right hemisphere, though it changed significantly in the cortex of the left hemisphere, returning to baseline values when the rabbits recovered from this state. Further analysis showed that this cortical interhemisphere asymmetry was based on the asymmetrical activity of individual neurons and small neuronal populations. Thus, changes in the structure of dependent relationships between neurons in microareas of the cortex of the left and macroareas of the cortex of the right hemisphere could be in different directions, while changes in microareas of the right hemisphere and macroareas of the left hemisphere were synergistic. Thus, asymmetry was detected at different levels of neuronal combinations (neuron pairs, micro- and macrogroups of neurons), which suggests mosaicism in neuron structure, which ultimately leads to overall functional asymmetry in "animal hypnosis." Some changes in the structure of dependent relationships between sensorimotor cortex neurons arising in "animal hypnosis" persisted or even became more marked after recovery of animals from this state.

EMG, force and discharge rate analysis of human jaw reflexes in response to axial stimulation of the incisor

Reflex studies utilising controlled stimulation along the long axis of human incisors are relatively new, and the effects that various stimulus parameters have on the elicited reflexes are not fully understood. Twelve subjects were recruited to determine the effects that contraction level, stimulus force and amount of constant force applied between stimuli have on the reflex response of the masseter muscle. Multi-unit intramuscular electromyogram (EMG) was recorded alongside surface EMG to determine whether any differences existed between the two. Furthermore, cumulative peri-stimulus "dischargegrams" were constructed to determine whether events seen in the EMG corresponded to changes in the discharge rates of the underlying motor units. Axial stimulation of the incisor induced a response in the EMG comprising of peak-trough-peak, with the trough being the most dominant. The bite force record showed only a reduction (relaxation) in response to the stimulation. The most significant experimental factor affecting the reflex occurrence and strength was the stimulus force. Although the latency, duration and occurrence rates were not significantly different, the strength of the responses was greater in intramuscular recordings compared with the surface recordings. Discharge rate analysis showed that approximately two-thirds of the late peaks detected in the EMG did not correspond to an increase in the discharge rates of the underlying units; hence they were due to the clustering of action potentials following the trough and not to a change in the membrane potential of the motoneurone. It was also found that the duration of the trough, as seen by the reduced cumulative discharge rate of the underlying units, was longer than indicated by the EMG.

Presumed Inhibitory Neurons in the Macaque Inferior Temporal Cortex: Visual Response Properties and Functional Interactions With Adjacent Neurons

Neurons in area TE of the monkey inferior temporal cortex respond selectively to images of particular objects or their characteristic visual features. The mechanism of generation of the stimulus selectivity, however, is largely unknown. This study addresses the role of inhibitory TE neurons in this process by examining their visual response properties and interactions with adjacent target neurons. We applied cross-correlation analysis to spike trains simultaneously recorded from pairs of adjacent neurons in anesthetized macaques. Neurons whose activity preceded a decrease in activity from their partner were presumed to be inhibitory neurons. Excitatory neurons were also identified as the source neuron of excitatory linkage as evidenced by a sharp peak displaced from the 0-ms bin in cross-correlograms. Most inhibitory neurons responded to a variety of visual stimuli in our stimulus set, which consisted of several dozen geometrical figures and photographs of objects, with a clear stimulus preference. On average, 10% of the stimuli increased firing rates of the inhibitory neurons. Both excitatory and inhibitory neurons exhibited a similar degree of stimulus selectivity. Although inhibitory neurons occasionally shared the most preferred stimuli with their target neurons, overall stimulus preferences were less similar between adjacent neurons with inhibitory linkages than adjacent neurons with common inputs and/or excitatory linkages. These results suggest that inhibitory neurons in area TE are activated selectively and exert stimulus-specific inhibition on adjacent neurons, contributing to shaping of stimulus selectivity of TE neurons.

A multi-channel correlation method detects traveling gamma-waves in monkey visual cortex

Correlations among simultaneously recorded signals are mostly analyzed pairwise and include temporal averaging. However, pairwise methods are not suitable for characterizing relationships among multiple channels for signals which vary temporally in an unpredictable way. Here we develop a time-resolved spatio-temporal correlation (STC) measure among simultaneously recorded signals. We demonstrate the capabilities of the method with artificial data sets and with multiple-channel recordings from striate cortex of awake monkeys. We concentrate on correlations in the gamma-frequency range (gamma: 30-90 Hz) because they were prominent in the analyzed recordings and gained high interest in the recent years due to their assumed role in associative processing, including perceptual binding. Former analyses of gamma-activities in visual cortex, using pairwise correlation methods, mostly revealed zero-delay correlation, indicating synchrony. In cat and monkey visual cortex this gamma-synchrony is restricted to 1.5-3.0 mm (half-height decline). However, our spatio-temporal correlation (STC)-method demonstrates for striate cortex from awake monkeys that gamma-synchrony is a local phenomenon of more global traveling plane waves that appear stimulus-induced at randomly varying orientations. These gamma-waves are coupled over much larger cortical distances (approximately 7 mm half-height decline) than the gamma-synchrony ranges obtained by pairwise correlation analyses from the same data. Our STC-method therefore suggests that the previously reported results of short-range and zero-delay correlations were often due to temporal averaging of traveling gamma-waves.

Functional connectivity of disparity-tuned neurons in the visual cortex

Different mechanisms have been proposed concerning how disparity-tuned neurons might be connected to produce the signals for depth perception. Here we present neurophysiological evidence providing insight on this issue. We have recorded simultaneously from pairs of disparity-tuned neurons in the cat's striate cortex. The purpose was to determine the relationships between disparity tuning and functional connectivity revealed through neural cross-correlograms. Monosynaptic connections tend to be stronger between pairs of cells with similar disparity tuning. Pairs of complex cells tend to have either similar tuning or nearly opposite tuning with an absence of quadrature relations. Pairs with at least one simple cell do have some nearly quadrature relationships when they are recorded from the same electrode. Coarse-to-fine connections (i.e., the presynaptic cell has lower disparity frequency and larger disparity range) tend to be stronger but less frequent than those of a fine-to-coarse nature. Our results are consistent with a system that produces weighted averaging across cells that are tuned to similar disparities but different disparity scales to reduce false matches.

Organization of neural networks in the neocortex

Implanted semimicroelectrodes were used in conscious cats to record spike discharges from groups of close-lying neurons, i.e., multineuron activity, in the deep layers of the frontal and motor areas of the cortex at different levels of food motivation. Spike activity was extracted from 4-7 neurons and interneuronal interactions were studied by cross-correlation analysis between neighboring neurons in each zone (local networks) and between neurons in two zones (distributed networks) with analysis epochs of 0-100 msec. The results showed that neurons in local networks can be divided into two subgroups: neurons with high-amplitude spikes and a predominance of output (divergent) connections and neurons with low-amplitude spikes and a predominance of input (convergent) connections. Local networks are based on powerful monosynaptic connections (with delays of up to 2 msec) between large and small neurons. Most connections in distributed networks were between small neurons in local networks of the frontal cortex and large neurons in local networks in the motor cortex. Food deprivation for 24 h mainly affected late (with delays of 2-100 msec) cross-correlation interneuronal relationships in both local and distributed networks.

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

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

Organization of the central control of muscles of facial expression in man

Surface EMGs were recorded simultaneously from ipsilateral pairs of facial muscles while subjects made three different common facial expressions: the smile, a sad expression and an expression of horror, and three contrived facial expressions. Central peaks were found in the cross-correlograms of EMG activity recorded from the orbicularis oculi and zygomaticus major during smiling, the corrugator and depressor anguli oris during the sad look and the frontalis and mentalis during the horror look. The size of the central peak was significantly greater between the orbicularis oculi and zygomaticus major during smiling. It is concluded that co-contraction of facial muscles during some facial expressions are accompanied by the presence of common synaptic drive to the motoneurones supplying the muscles involved. Central peaks were found in the cross-correlograms of EMG activity recorded from the frontalis and depressor anguli oris during a contrived expression. However, no central peaks were found in the cross-correlograms of EMG activity recorded from the frontalis and orbicularis oculi or from the frontalis and zygomaticus major during the other two contrived expressions. It is concluded that a common synaptic drive is not present between all possible facial muscle pairs and suggests a functional role for the synergy. The origin of the common drive is discussed. It is concluded that activity in branches of common stem last-order presynaptic input fibres to motoneurones innervating the different facial muscles and presynaptic synchronization of input activity to the different motoneurone pools is involved. The former probably contributes more to the drive to the orbicularis oculi and zygomaticus major during smiling, while the latter is probably more prevalent in the corrugator and depressor anguli oris during the sad look, the frontalis and mentalis during the horror look and the frontalis and depressor anguli oris during one of the contrived expressions. The strength of common synaptic drive is inversely related to the degree of separate control that can be exhibited by the facial muscles involved.

Connections between respiratory neurones in the neonatal rat transverse medullary slice studied with cross-correlation

In the transverse medullary slice prepared from neonatal rats the hypoglossal nerve rootlets exhibit a bursting 'respiratory' rhythm as do neurones in the pre-Bötzinger complex (PBC). We used cross-correlation analysis of the rhythmic multiunit discharges recorded from hypoglossal nerve rootlets, hypoglossal nucleus neurones and PBC neurones to investigate the connections between these groups. All cross-correlograms computed between left and right hypoglossal nerves, and between hypoglossal neurones and contralateral hypoglossal nerves, displayed central peaks with broad half-amplitude widths (mean +/- S.D. of 29.6 +/- 10.4 and 37.3 +/- 6.0 ms, respectively), which we interpreted as evidence for activation from a common source. Five of the 18 cross-correlograms computed between left and right PBC neurones displayed peaks either side of time zero with narrower half-amplitude widths (mean +/- S.D. of 9.3 +/- 1.9 ms) superimposed on broader central peaks, which we interpreted as evidence for mutual excitation and common activation, respectively. Cross-correlograms computed between PBC neurones and contralateral hypoglossal neurones or nerves did not display consistent features, but some of those computed between PBC and ipsilateral hypoglossal neurones (two of eight) or nerves (two of five) displayed peaks with broad half-amplitude widths (mean +/- S.D. of 36.8 +/- 6.9 ms), offset from time zero by 6 ms (except for one at 18 ms), which we interpreted as evidence for excitation of hypoglossal neurones and motoneurones by PBC neurones. We concluded that rhythm is synchronised between left and right sides by mutual excitatory connections between left and right PBC neurones. The rhythm is transmitted to ipsilateral hypoglossal neurones by a paucisynaptic pathway. Both hypoglossal neurones and PBC neurones receive a common activation from as yet unidentified sources.

No parallel fiber volleys in the cerebellar cortex: evidence from cross-correlation analysis between Purkinje cells in a computer model and in recordings from anesthetized rats

Purkinje cells aligned on the medio-lateral axis share a large proportion of their approximately 175,000 parallel fiber inputs. This arrangement has led to the hypothesis that movement timing is coded in the cerebellum by beams of synchronously active parallel fibers. In computer simulations I show that such synchronous activation leads to a narrow spike cross-correlation between pairs of Purkinje cells. This peak was completely absent when shared parallel fiber input was active in an asynchronous mode. To determine the presence of synchronous parallel fiber beams in vivo I recorded from pairs of Purkinje cells in crus IIa of anesthetized rats. I found a complete absence of precise spike synchronization, even when both cells were strongly modulated in their spike rate by trains of air-puff stimuli to the face. These results indicate that Purkinje cell spiking is not controlled by volleys of synchronous parallel fiber inputs in the conditions examined. Instead, the data support a model by which granule cells primarily control Purkinje cell spiking via dynamic population rate changes.

Pharmacologically induced enhancement of recurrent inhibition in humans: effects on motoneurone discharge patterns

The aim of the present study was to investigate the effects of spinal recurrent inhibition on human motoneurone discharge patterns. The tonic discharge activity of motor unit pairs was recorded in the extensor carpi radialis (ECR) and abductor digiti minimi (ADM) muscles during voluntary isometric contraction. While undergoing continuous intravenous saline (NaCl 0.9 %) perfusion, the subjects were given a short lasting injection of L-acetylcarnitine (L-Ac), which has been found to potentiate recurrent inhibition in humans. The variability, synchronization and coherence of the motor unit discharges were analysed during four successive test periods (lasting 2-3 min each). A significant decrease in the inter-spike interval (ISI) coefficient of variation was observed in the discharge patterns of the motor units tested in the ECR and not in the ADM, which were not accompanied by any consistent changes in the mean ISIs of the motor unit activity in either muscle. The L-Ac injection also led to a significant increase in the synchronization in half of the motor unit pairs tested in the ECR muscle (n = 29), whereas no consistent changes were observed with the ADM motor units (n = 25). However, coherence analysis failed to reveal any consistent differences in the incidence of significant values of coherence spectrum between the pre-injection and injection periods among the motor unit pairs tested with either saline or L-Ac injections, in either the ECR or ADM muscles. The contrasting effects on the variability and the synchronization of the motor unit discharges observed with ECR motoneurones known to undergo recurrent inhibition and with ADM motoneurones known to lack recurrent inhibition suggest that the drug may have specific effects which are mediated by an enhancement of the Renshaw cell activity. The decrease in the ISI variability is in line with the hypothesis that recurrent inhibition may contribute along with the post-spike after-hyperpolarization to limiting the influence of the synaptic noise on the firing times of steadily discharging motoneurones. The present data, which suggest that recurrent inhibition plays a synchronizing rather than a desynchronizing role, are in keeping with the fact that the Renshaw cells may provide an important source of common inhibitory inputs.

Stereoscopic depth processing in the visual cortex: a coarse-to-fine mechanism

For binocular animals viewing a three-dimensional scene, the left and right eyes receive slightly different information, and the brain uses this 'binocular disparity' to interpret stereoscopic depth. An important theoretical conjecture in this mechanism is that coarse processing precedes and constrains finely detailed processing. We present three types of neurophysiological data from the cat's visual cortex that are consistent with a temporal coarse-to-fine tuning of disparity information. First, the disparity tuning of cortical cells generally sharpened during the time course of response. Second, cells responsive to large and small spatial scale had relatively shorter and longer temporal latencies, respectively. Third, cross-correlation analysis between simultaneously recorded pairs of cortical cells showed that connections between disparity-tuned neurons were generally stronger for coarse-to-fine processing than for fine-to-coarse processing. These results are consistent with theoretical and behavioral studies and suggest that rapid, coarse percepts are refined over time in stereoscopic depth perception.

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

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

Firing modes of midbrain dopamine cells in the freely moving rat

There is a large body of data on the firing properties of dopamine cells in anaesthetised rats or rat brain slices. However, the extent to which these data relate to more natural conditions is uncertain, as there is little quantitative information available on the firing properties of these cells in freely moving rats. We examined this by recording from the midbrain dopamine cell fields using chronically implanted microwire electrodes. (1) In most cases, slowly firing cells with broad action potentials were profoundly inhibited by the dopamine agonist apomorphine, consistent with previously accepted criteria. However, a small group of cells was found that were difficult to classify because of ambiguous combinations of properties. (2) Presumed dopamine cells could be divided into low and high bursting (>40% of their spikes in bursts) groups, with the majority having low bursting rates. The distribution of burst incidence was similar to that previously reported with chloral hydrate anaesthesia, but the average intraburst frequency was higher in the conscious animal at rest and was higher again in bursts triggered by salient stimuli. (3) There was no evidence for spike frequency adaptation within bursts on average, consistent with the hypothesis that afterhyperpolarisation currents may be disabled during behaviourally induced bursting. (4) Presumed dopamine cells responded to reward-related stimuli with increased bursting rates and significantly higher intraburst frequencies compared to bursts emitted outside task context, indicating that modulation of afferent activity might not only trigger bursting, but may also regulate burst intensity. (5) In addition to the irregular single spike and bursting modes we found that extremely regular (clock-like) firing, previously only described for dopamine cells in reduced preparations, can also be expressed in the freely moving animal. (6) Cross-correlation analysis of activity recorded from simultaneously recorded neurones revealed coordinated activity in a quarter of dopamine cell pairs consistent with at least "functional" connectivity. On the other hand, most dopamine cell pairs showed no correlation, leaving open the possibility of functional sub-groupings within the dopamine cell fields. Taken together, the data suggest that the basic firing modes described for dopamine cells in reduced or anaesthetised preparations do reflect natural patterns of activity for these neurones, but also that the details of this activity are dependent upon modulation of afferent inputs by behavioural stimuli.

Functional connectivity between neurons generating resting discharge in renal sympathetic neurons in the rabbit

In anaesthetised rabbits, we analysed the cross-correlations of resting discharges in pairs of simultaneously recorded renal neurons. The study of 428 significant cross-correlations showed three pure types of neuronal co-ordination and two combined effects. Pure shared input was observed in 202 pairs of neurons (47.2%). The width of cross-correlogram peak was 222.9 +/- 7.2 ms. The peak height of shared input was 4.8 +/- 0.2 spikes s(-1) and was significantly related with correlogram width. Pure excitatory connection was found in 24 neuronal pairs (5.6%). Its peak width was 7.9 +/- 1.9 ms and it was shifted from time zero by 13.3 +/- 3.4 ms. The amplitude of the peak was 19.6 +/- 4.6 spikes s(-1). In 165 pairs of neurons (38.6%), excitatory connection was combined with shared input component. Pure reciprocal activation was seen in only 5 pairs of neurons (1.1%) while in 32 pairs (7.5%) it was accompanied by shared input co-ordination. The distance between narrow peaks measured in the combined type of co-ordination was 32.5 +/- 3.5 ms and the mean peak height was 4.1 +/- 0.7 spikes s(-1). Shared input in pure form and that accompanying both excitatory connection and reciprocal activation were significantly related to the frequency of discharge of neurons contributing to the cross-correlogram. The incidence and basic properties of the encountered types of neural coordination may indicate a pattern of interconnections between cells generating vasomotor tone in renal sympathetic neurons.

Linked spike activity of neurons in the right and left lateral hypothalamus in conditions of food motivation

Cross- and autocorrelation histograms were constructed to investigate the linkage of spike activity from neurons in the right and left lateral hypothalamus, recorded in conditions of restful waking, after 24-h food deprivation, and after satiation. In conditions of hunger, there was a preferred order in which neurons fired in these two hypothalamic regions: a significant majority of cases (59%) showed initial firing by a neuron in the left hypothalamus followed by a neuron in the right, with delays of up to 200 msec and peaks in the delay plot at 30 and 160 msec. The opposite firing order was seen less frequently (21% of cases). In conditions of hunger, linked spike activity showed an increase in the probability of detecting frequencies in the theta range (from 11% to 29%).

Respiratory pre-motor control of hypoglossal motoneurons in the rat

The goal of this study was to determine the origin and transmission pathway of respiratory drive to hypoglossal motoneurons. First we recorded intracellularly from 28 antidromically activated inspiratory hypoglossal motoneurons (resting membrane potential, -50+/-3 mV), and found that injection of chloride ions had no discernible effect on the shape of their membrane potential trajectories. We concluded that the membrane potential trajectories of these hypoglossal motoneurons were determined primarily by inspiratory excitation. To determine the origin of this excitation we cross-correlated the extracellular discharge of medullary inspiratory neurons, including those in the hypoglossal motor nucleus, with the hypoglossal nerve discharge. We found 27 inspiratory neurons within the hypoglossal motor nucleus that were not antidromically activated from the ipsilateral hypoglossal nerve; their cross-correlograms featured either central peaks (1.7+/-0.2 ms) alone (n=14; 39%), or central peaks (1.3+/-0.2 ms) followed by troughs (1.3+/-0.1 ms) at short latencies (1.1+/-0.4 ms) (n=13; 36%), and suggest that these neurons are hypoglossal interneurons. We recorded from 238 inspiratory neurons throughout the rest of the medulla; the cross-correlograms of 19 neurons (8%), located mostly in the lateral tegmental field, displayed narrow half-amplitude peaks (1.0+/-0.1 ms) at short latencies (0.9+/-0.1 ms), which we interpreted as evidence for monosynaptic excitation of hypoglossal motoneurons.We conclude that the respiratory control of hypoglossal motoneurons originates from inspiratory premotor neurons scattered throughout the lateral tegmental field and interneurons within the hypoglossal motor nucleus.

Synchronized neuronal discharge in the basal ganglia of parkinsonian patients is limited to oscillatory activity

It has been proposed that an increase in synchronization between neurons in the basal ganglia contributes to the clinical features of Parkinson's disease (PD). To examine this hypothesis, we looked for correlations in the discharge activity of pairs of neurons in the globus pallidus internus (GPi), globus pallidus externus (GPe), and the substantia nigra pars reticulata (SNr). Recordings were performed in PD patients undergoing functional stereotactic mapping for pallidotomy (eight patients) or subthalamic nucleus deep brain stimulation (four patients). A double-microelectrode setup was used to simultaneously record from neurons separated by distances as small as 250 microm. In the five pallidotomy patients without limb tremor during the procedure, none of the 73 GPi pairs and 15 GPe pairs displayed synchronous activity. In the three pallidotomy patients with limb tremor, 6 of 21 GPi pairs and 5 of 29 GPe pairs displayed oscillatory synchronization in the frequency range of the ongoing limb tremor (3-6 Hz) or at higher frequencies (15-30 Hz). Synchronized activity was not observed in the SNr (10 pairs). The findings indicate that oscillatory synchronization between pairs of GPi or GPe neurons is found in patients with limb tremor. These results also suggest that overt neuronal synchronization, which may be attributable to an increase in direct synaptic connections or common collateral afferent inputs, is not present in the basal ganglia of patients with PD.

Noise and the PSTH response to current transients: II. Integrate-and-fire model with slow recovery and application to motoneuron data

A generalized version of the integrate-and-fire model is presented that qualitatively reproduces firing rates and membrane trajectories of motoneurons. The description is based on the spike-response model and includes three different time constants: the passive membrane time constant, a recovery time of the input conductance after each spike, and a time constant of the spike afterpotential. The effect of stochastic background input on the peristimulus time histogram (PSTH) response to spike input is calculated analytically. Model results are compared with the experimental data of Poliakov et al. (1996). The linearized theory shows that the PSTH response to an input spike is proportional to a filtered version of the postsynaptic potential generated by the input spike. The shape of the filter depends on the background activity. The full nonlinear theory is in close agreement with simulated PSTH data.

Effects of common excitatory and inhibitory inputs on motoneuron synchronization

We compared the effects of common excitatory and inhibitory inputs on motoneuron synchronization by simulating synaptic inputs with injected current transients. We elicited repetitive discharge in hypoglossal motoneurons recorded in slices of rat brain stem using a combination of a suprathreshold injected current step with superimposed noise to mimic the synaptic drive likely to occur during physiological activation. The effects of common inputs to motoneurons were simulated by the addition of a waveform composed of from 6 to 300 trains of current transients designed to mimic excitatory and/or inhibitory synaptic currents. We compared the discharge records obtained in several trials in which the same "common input" waveform was applied repeatedly in the presence of different background noise waveforms. The effects of the common input on motoneuron discharge probability and discharge rate were determined by compiling a cross-correlation histogram (CCHist) and a perispike frequencygram (PSFreq) between the discharges of the same cell at different times. Both excitatory and inhibitory common inputs induced synchronous discharge that was evident by a large central peak in the CCHist. The CCHists produced by common excitatory inputs were characterized by larger and narrower central peaks than those generated by common inhibitory inputs. The PSFreqs produced by common excitatory inputs indicated an increase in the discharge rate of motoneurons around time 0 that coincided with the narrow and large central peak in the CCHist. On the other hand, inhibitory inputs often generated very little, if any, change in the discharge rate around time 0 corresponding with the small and wide central peak in the CCHist. These results suggest that the CCHist indicates the effective strength of the net common input but not its sign. Although correlated changes in discharge rate are often quite different for net excitatory and inhibitory common input, except in some restricted conditions, the PSFreq analysis also cannot be used to unambiguously distinguish net excitation from net inhibition.