Correlation of spatial event plots with simulated population responses of mechanoreceptive fibers

The experimental setup for generating spatial event plots (SEPs) from single mechanoreceptive fibers of the skin was computationally simulated. The generic fibers used in the simulations were similar to the rapidly adapting fibers (RAs), and had variable refractoriness and receptive-field size. The speed, lateral shift, and the contact width of the drum scanned across the receptive field of the fiber are adjustable parameters. The stimulus patterns used on the drum mimicked stimuli used by several other investigators. These were dot patterns, grating patterns, and the letter "E". First, the effects of simulation parameters on the SEPs were studied. The simulation output confirms the results of physiological experiments that SEPs contain information on the spatiotemporal resolution of the fiber. The next series of simulations involved generating SEPs of fibers obtained from the same or varying spatial distributions of receptive fields. Three hypothetical distributions were used: homogeneous rectangular, uniformly random, and Gaussian. The momentary population response at each case was found using the technique by Johansson and Vallbo (Brain Res 184: 353-366, 1980). The population responses were not isomorphic images of the stimulus patterns due to the variations in field sizes and locations. However, every fiber, no matter which distribution it came from, generated almost identical SEPs given similar response properties. Furthermore, the SEPs looked like the outline of the stimulus. These observations show that SEPs do not contain information about the population response. Therefore, reconstructing the population response using SEPs can result in misleading conclusions on central-nervous-system processing and should be viewed cautiously when formulating psychophysical/physiological linking hypotheses.     

Asymmetric response properties of rapidly adapting mechanoreceptive fibers in the rat glabrous skin

Previous histological and neurophysiological studies have shown that the innervation density of rapidly adapting (RA) mechanoreceptive fibers increases towards the fingertip. Since the psychophysical detection threshold depends on the contribution of several RA fibers, a high innervation density would imply lower thresholds. However, our previous human study showed that psychophysical detection thresholds for the Non-Pacinian I channel mediated by RA fibers do not improve towards the fingertip. By recording single-unit spike activity from rat RA fibers, here we tested the hypothesis that the responsiveness of RA fibers is asymmetric in the proximo-distal axis which may counterbalance the effects of innervation density. RA fibers (n?=?32) innervating the digital glabrous skin of rat hind paw were stimulated with 40-Hz sinusoidal mechanical bursts at five different stimulus locations relative to the receptive field (RF) center (two distal, one RF center, two proximal). Different contactor sizes (area: 0.39, 1.63, 2.96?mm²) were used. Rate-intensity functions were constructed based on average firing rates, and the absolute spike threshold and the entrainment threshold were obtained for each RA fiber. Thresholds for proximal stimulus locations were found to be significantly higher than those for distal stimulus locations, which suggests that the mechanical stimulus is transmitted better towards the proximal direction. The effect of contactor size was not significant. Mechanical impedance of the rat digital glabrous skin was further measured and a lumped-parameter model was proposed to interpret the relationship between the asymmetric response properties of RA fibers and the mechanical properties of the skin.     

Distribution of the intensity-characteristic parameters of cat rapidly adapting mechanoreceptive fibers

Modeling population responses of nerve fibers requires statistical characterization of fiber-response properties. The rate/intensity characteristics of cat rapidly adapting (RA) fibers were fitted by four-parameter, piece-wise linear functions using nonlinear regression (n = 14; R2 > 0.958). The parameters were tested against the null hypothesis that they are log normally distributed. The test fail to reject this hypothesis (Kolmogorov-Smirnov p>0.477). However, a significant statistical difference was found between the specific lognormal distributions obtained from monkey (Johnson, J Neurophysiol 37: 48-72, 1974) and cat for all four parameters (Kolmogorov-Smirnov, p<0.0075, p<0.05, p<0.0001, p<0.00007). Although the stimulus contactor size was not the same in monkey and cat studies, the differences between monkey and cat fibers are attributed to anatomical differences in the glabrous sin of both species. Modeling studies suggest that the absolute firing thresholds of RA fibers have a right-skewed distribution because of the anatomical constraints present in both species' skin. Meissner corpuscles, which are the sensory end-organs of RA fibers, are likely to be found deeper in the skin within dermal papilla, therefore, the thresholds can be elevated. However, the thresholds are bounded at lower end, probably due to the epidermal junction that acts as a superficial mechanical barrier for these corpuscles.     

Vibrotactile thresholds of the Non-Pacinian I channel: II. Predicting the effects of contactor location on the phalanx

The firing-rate-based population model for rapidly-adapting (RA) mechanoreceptive fibers by Güçlü and Bolanowski (Güçlü B, Bolanowski SJ. 2002. Modeling population responses of rapidly-adapting mechanoreceptive fibers. Journal of Computational Neuroscience 12:201–218.) is extended by including temporal-response properties of RA fibers. This representation allows for the generation of action-potential (spike) times for each fiber when a sinusoidal, steady-state stimulus is applied onto the skin. Signal detection theory was used to predict human psychophysical thresholds. Specifically, the effects of sensorineural innervation pattern, stimulus-contactor location and selected decision rules on the model predictions were studied. The predicted thresholds were lowest when the decision rule was one spike and highest when many active fibers were required for detection. These predictions were empirically tested by measuring vibrotactile thresholds of the Non-Pacinian I (NP I) channel, which required the special techniques discussed in the preceding article. Although the model predicted thresholds to decrease distally due to the known innervation density which is higher distally, the thresholds of the NP I psychophysical channel were found to be approximately constant (20–25?dB re 1?µm peak amplitude) from the proximal site on the terminal phalanx to the most distal portion. Interestingly, the mechanical impedance of the skin was found not to be constant along the proximo-distal axis. This latter result implies that the space-invariant mechanical attenuation function used in the model may not be valid at every location on the fingertip. Because of this, the discrepancy between the model's predictions and the psychophysical results may be reconciled.     

Vibrotactile thresholds of the Non-Pacinian I channel: II. Predicting the effects of contactor location on the phalanx

The firing-rate-based population model for rapidly-adapting (RA) mechanoreceptive fibers by Gü?lü and Bolanowski is extended by including temporal-response properties of RA fibers. This representation allows for the generation of action-potential (spike) times for each fiber when a sinusoidal, steady-state stimulus is applied onto the skin. Signal detection theory was used to predict human psychophysical thresholds. Specifically, the effects of sensorineural innervation pattern, stimulus-contactor location and selected decision rules on the model predictions were studied. The predicted thresholds were lowest when the decision rule was one spike and highest when many active fibers were required for detection. These predictions were empirically tested by measuring vibrotactile thresholds of the Non-Pacinian I (NP I) channel, which required the special techniques discussed in the preceding article. Although the model predicted thresholds to decrease distally due to the known innervation density which is higher distally, the thresholds of the NP I psychophysical channel were found to be approximately constant (20-25 dB re 1 microm peak amplitude) from the proximal site on the terminal phalanx to the most distal portion. Interestingly, the mechanical impedance of the skin was found not to be constant along the proximo-distal axis. This latter result implies that the space-invariant mechanical attenuation function used in the model may not be valid at every location on the fingertip. Because of this, the discrepancy between the model's predictions and the psychophysical results may be reconciled.     

Sensory Innervation of the Raccoon Forepaw: 1. Receptor Types in Glabrous and Hairy Skin and Deep Tissue

Electrophysiological recordings were made from the median, ulnar, radial, and dorsal ulnar nerves to determine the types of mechanosensory receptors serving glabrous and hairy skin surfaces of the raccoon forepaw. In addition to the cutaneous mechanoreceptors, fibers innervating deep tissues were also recorded from each of these nerves. These included sensory fibers innervating muscles, joints, claws, and the subcutaneous pulp.

The array of receptors serving raccoon glabrous skin was the same as found in monkeys and humans: Rapidly adapting (RA), slowly adapting (SA), and Pacinian (Pc) fibers were characterized. Pacinian fibers have been rarely described in previous physiological studies of the raccoon peripheral nerves, but in the present study they composed between 14% and 18% of the glabrous skin mechanoreceptors recorded. A distal-proximal gradient in the density of skin innervation was evident for all three types of receptors.

Receptors characterized in the hairy skin of the dorsal paw were similar to those described in other mammals, and included both down and guard hair afferents, non-hair-associated RA fibers, and SA I and SA II fibers. The relative proportions of these fibers differed from those generally reported for the hairy skin of other mammals. SA hair-associated afferent fibers, which have been reported previously only in primate hairy skin, were also found in large numbers in the raccoon. Similarities and differences in the frequency and types of receptors innervating the raccoon forepaw, the forepaws of other mammals, and the hands of primates (including humans) are discussed.     

Modulatory effects of parallel fiber and molecular layer interneuron synaptic activity on purkinje cell responses to ascending segment input: a modeling study

Based on anatomical, physiological, and model-based studies, it has been proposed that synapses associated with the ascending segment of granule cell axons provide the principle excitatory drive on Purkinje cells which is then modulated by the more numerous parallel fiber synapses. In this study we have evaluated this idea using a detailed compartmental model of a cerebellar Purkinje cell by providing identical ascending segment synaptic inputs during different levels of random parallel fiber and molecular interneuron input. Results suggest that background inputs from parallel fibers and molecular layer interneurons can have a substantial effect on the response of Purkinje cells to ascending segment inputs. Interestingly, these effects are not reflected in the average firing rate of the Purkinje cell and are thus entirely dendritic in effect. These results are considered in the context of the known segregated spatial distribution of the parallel fibers and ascending segment synapses and a new hypothesis concerning the functional organization of cerebellar cortical circuitry.     

Modeling the response of a population of olfactory receptor neurons to an odorant

We modeled the firing rate of populations of olfactory receptor neurons (ORNs) responding to an odorant at different concentrations. Two cases were considered: a population of ORNs that all express the same olfactory receptor (OR), and a population that expresses many different ORs. To take into account ORN variability, we replaced single parameter values in a biophysical ORN model with values drawn from statistical distributions, chosen to correspond to experimental data. For ORNs expressing the same OR, we found that the distributions of firing frequencies are Gaussian at all concentrations, with larger mean and standard deviation at higher concentrations. For a population expressing different ORs, the distribution of firing frequencies can be described as the superposition of a Gaussian distribution and a lognormal distribution. Distributions of maximum value and dynamic range of spiking frequencies in the simulated ORN population were similar to experimental results.     

Assessment of the non-Gaussianity and non-linearity levels of simulated sEMG signals on stationary segments

The purpose of this paper was to evaluate the effects of the longitudinal single differential (LSD), the longitudinal double differential (LDD) and the normal double differential (NDD) spatial filters, the electrode shape, the inter-electrode distance (IED) on non-Gaussianity and non-linearity levels of simulated surface EMG (sEMG) signals when the maximum voluntary contraction (MVC) varied from 10% to 100% by a step of 10%. The effects of recruitment range thresholds (RR), the firing rate (FR) strategy and the peak firing rate (PFR) of motor units were also considered.A cylindrical multilayer model of the volume conductor and a model of motor unit (MU) recruitment and firing rate were used to simulate sEMG signals in a pool of 120 MUs for 5 s. Firstly, the stationarity of sEMG signals was tested by the runs, the reverse arrangements (RA) and the modified reverse arrangements (MRA) tests. Then the non-Gaussianity was characterised with bicoherence and kurtosis, and non-linearity levels was evaluated with linearity test.The kurtosis analysis showed that the sEMG signals detected by the LSD filter were the most Gaussian and those detected by the NDD filter were the least Gaussian. In addition, the sEMG signals detected by the LSD filter were the most linear. For a given filter, the sEMG signals detected by using rectangular electrodes were more Gaussian and more linear than that detected with circular electrodes. Moreover, the sEMG signals are less non-Gaussian and more linear with reverse onion-skin firing rate strategy than those with onion-skin strategy. The levels of sEMG signal Gaussianity and linearity increased with the increase of the IED, RR and PFR.     

Patterns in the discharge of simple and complex visual cortical cells

The activity of visual cortical neurons (area 17) was recorded in anaesthetized cats in response to sinusoidal drifting gratings. The statistical structure of the discharge of simple and complex cells has been studied as a function of the various parameters of a drifting grating: spatial frequency, orientation, drifting velocity and contrast. For simple cells it has been found that the interspike interval distributions in response to drifting gratings of various spatial frequencies differ only by a time scale factor. They can be reduced to a unique distribution by a linear time transformation. Variations in the spatial frequency of the grating induce variations in the mean firing rate of the cell but leave unchanged the statistical structure of the discharge. On the contrary, the statistical structure of the simple cell activity changes when the contrast or the velocity of the stimulus is varied. For complex cells it has been found that the invariance property described above for simple cells is not valid. Complex cells present in their activity in response to visual stimuli two different firing patterns: spikes organized in clusters and spikes that do not show this organization ('isolated spikes'). The clustered component is the only component of the complex cell discharge that is tuned for spatial frequency and orientation, while the isolated spike component is correlated with the contrast of the stimulus.     

Frequency responses of cat rapidly adapting mechanoreceptive fibers

The frequency responses of 11 rapidly adapting (RA) fibers in cat were studied by representing the average firing rate as a function of sinusoidal stimulus amplitude and stimulus frequency. Specifically, rate-intensity functions at different stimulation frequencies were fitted by four-parameter (a0, a1, a2, a3), piece-wise linear functions using nonlinear regression (n = 59; R2 > 0.877). Rate-intensity functions at intermediate frequencies were found by linear interpolation. The result of this analysis is rate-amplitude-frequency functions plotted as two-dimensional surfaces. The surfaces consist of five regions separated and sufficiently defined by four space curves. At 14 different frequencies, the statistical distribution of each rate-intensity-function parameter could be approximated by a particular lognormal distribution (n = 56; R2 > 0.796). The Kolmogorov-Smirnov test fails to reject this hypothesis for each combination of frequency and parameter (56 tests; p > 0.39). Therefore, at a given frequency, the variation of the parameters can be represented by lognormal distributions with specific means and standard deviations. Responses of six RA fibers, which are different from the data-set used for modeling, were compared with the stochastic model at different frequencies. The parameters of those fibers were tested against the null hypotheses that they were sampled from the particular parameter distributions dictated by the model. The Kolmogorov-Smirnov test fails to reject all the hypotheses at the alpha = 0.05 level (44 tests). At the alpha = 0.10 level, only a few test parameters were found to be departing from the model (a0 and a1 at 5 Hz; a2 at 20 Hz; a2 and a3 at 50 Hz). The remaining test parameters could be accurately described by the model. Having confirmed the validity of the model, the logarithmic means and the logarithmic standard deviations of the lognormally distributed rate-intensity-function parameters were estimated in the frequency range of 4-200 Hz. The rate-amplitude-frequency surfaces sampled from the established stochastic model completely characterize the rate responses of RA fibers to sinusoidal stimuli and are superior to tuning curves which require selecting criterion responses. The current rate-response model is promising for future computational work, especially on population modeling.     

Mesenteric and tactile Pacinian corpuscles are anatomically and physiologically comparable

Pacinian corpuscles (PCs) in cat mesentery have been studied extensively to help determine the structural and functional bases of tactile mechanotransduction. Although we, like many other investigators, have found that the mesenteric receptors are anatomically very similar to those found in mammalian skin, few physiological characteristics of the mesenteric PCs and those of the skin have been compared. Action-potential rate-amplitude and frequency characteristics (10 Hz–1 KHz), as well as interval (IH) and peri-stimulus-time (PSTH) histograms in response to sinusoidal displacements were obtained from nerve fibers innervating mesenteric PCs and from PC fibers innervating cat glabrous skin. The intensity characteristics obtained on both preparations showed similar response profiles, including equal slopes for low stimulus intensities (approximately 10, with impulse ratios/20 dB displacement) and one and two impulse/cycle entrainment. The frequency characteristics of both groups were U-shaped with similar low-frequency slopes (?12.5 dB/octave) and bandwidths (Q_3dB = 1.4). The best frequency for both the tactile PCs' and mesenteric PCs was 250 Hz, which is in the expected range. The IHs showed entrainment and the PSTHs showed neither transient responses nor adaptation to steady-state sinusoidal stimuli. The functional similarity between mesenteric PCs' nerve responses and those of tactile PC afferents, as well as the receptors’ anatomical similarity, lead us to suggest that the mesenteric PC can act as a model for those in the skin. Furthermore, since the frequency characteristics of the two PC types are similar, it is concluded that the skin, while attenuating stimulus intensity, does not impart temporal filtering of vibratory stimuli.     

Localization of Merkel cells in the monkey skin: an anatomical model

The Merkel cell-neurite complex is considered to be one class of mechanoreceptors in the skin. Merkel cells are innervated by slowly adapting type I (SAI) tactile nerve fibers. In this paper, the detailed distribution of Merkel cells is studied by immunohistochemical labeling of the monkey (Macaca fascicularis) digital glabrous skin. Specific morphometric variables (density of intermediate epidermal ridges and Merkel cells, distance between skin surface and ridge tips and bases, maximum and average cell counts per ridge, distance between cells and ridges) were measured by a combination of light/fluorescence microscopy and computer-image analysis. The morphometric results were similar for each digit of the monkey's hand. Next, the anatomical data were used to form a three-dimensional reconstruction of the Merkel-cell distribution in the fingertip skin. A patch of the distal-pad surface was then computationally flattened to obtain the two-dimensional distribution of Merkel cells. Based on previous anatomical and physiological data, SAI fibers were simulated to innervate clusters of Merkel cells in the distal-pad surface. On average, 28 cells were innervated by a single fiber. The resulting anatomical model may be used to estimate the population response of SAI fibers by incorporating spike generation.     

Organization of Meissner corpuscles in the glabrous skin of monkey and cat

Reconstructing neural-population responses in the form of spatial event plots assumes that the receptors are organized in a dense linear array. We have found that this assumption is not valid by determining the spatial organization of Meissner corpuscles (MCs) in the glabrous skin of both cat and monkey. The tissue was excised from animals that had been cardiac perfused with 4% paraformaldehyde. One-micrometer plastic sections revealed that the morphology of these receptors is different in the two animal species. However, in both species, they reside in approximately the same place in the dermal pegs of the skin, between the epidermal ridges, and electrophysiologically they both respond to ramp-and-hold stimuli with a rapidly adapting firing pattern. Thus, in this study we will refer to the receptors of the cat as "Meissner-like". In monkey, MCs are located in the dermal papillae between the epidermal limiting and intermediate ridges, forming orderly rows, the contours of which follow the overlying fingerprint. Although the average density of MCs is 45/mm2, they are distributed along the dermal pegs in such a manner as to give rise to three significantly ( p < 0.017) different average distances between corpuscles. We note that "fingerprints" vary in topography across the hand and this is also reflected in the underlying MC arrays. In the cat, these "Meissner-like" receptors display no specific pattern and have a density much lower than in the monkey. Cat glabrous skin does not have "fingerprints". The results emphasize that the spatial organization of tactile receptors must be taken into account when interpreting reconstructed population responses.     

An analysis of a sample of afferent discharges provoked in frog dorsal root by a passive movement

Summary Discharges of individual fibers from the IX-th frog dorsal root, under stretching of the hind leg by falling weights, have been recorded. About 1/3 of the recorded fibers are leading off from rapidly adapting receptors. The fibers coming from slowly adapting receptors and those not responding to the applied stimulus are analyzed with some details.The number of pulse per second as a function of time and load is examined. It has been found that a group of fibers increases the rate of firing regularly with the load, a second group gives the greatest response in intermediate zones of the loads scale and a third group shows a higher threshold.Some remarks on the aggregate message reaching the spinal cord are developed.     

Directionality of auditory nerve fiber responses to pure tone stimuli in the grassfrog, Rana temporaria. II. Spike timing

We studied the directionality of spike timing in the responses of single auditory nerve fibers of the grass frog, Rana temporaria, to tone burst stimulation. Both the latency of the first spike after stimulus onset and the preferred firing phase during the stimulus were studied. In addition, the directionality of the phase of eardrum vibrations was measured. The response latency showed systematic and statistically significant changes with sound direction at both low and high frequencies. The latency changes were correlated with response strength (spike rate) changes and were probably the result of directional changes in effective stimulus intensity. Systematic changes in the preferred firing phase were seen in all fibers that showed phaselocking (i.e., at frequencies below 500–700 Hz). The mean phase lead for stimulation from the contralateral side was approximately 140° at 200 Hz and decreased to approximately 100° at 700 Hz. These phaseshifts correspond to differences in spike timing of approximately 2 ms and 0.4 ms respectively. The phaseshifts were nearly independent of stimulus intensity. The phase directionality of eardrum vibrations was smaller than that of the nerve fibers. Hence, the strong directional phaseshifts shown by the nerve fibers probably reflect the directional characteristics of extratympanic pathways. Accepted: 23 November 1996     

The frequency selectivity of information-processing channels in the tactile sensory system

The frequency selectivity of the P, NP I, and NP II channels of the four-channel model of mechanoreception for glabrous skin was measured psychophysically by an adaptation tuning curve procedure. The results substantially extend the frequency range over which the frequency selectivity of these channels is known and further confirm the hypothesis that the input stage of each of these channels consists of specific sensory nerve fibers and associated receptors. Specifically, the frequency characteristics of Pacinian nerve fibers, rapidly adapting (RA) nerve fibers, and slowly adapting Type II (SA II) nerve fibers were found to be the peripheral neurophysiological correlates of the P, NP I, and NP II channels, respectively. The finding that the tuning characteristic for a test stimulus of 250 Hz delivered through a small (0.008 cm2) contactor depended dramatically on the duration of the test stimulus whereas the detection threshold did not, provides new evidence in support of the hypothesis that separate NP II and P channels exist.     

Mesenteric and tactile Pacinian corpuscles are anatomically and physiologically comparable

Pacinian corpuscles (PCs) in cat mesentery have been studied extensively to help determine the structural and functional bases of tactile mechanotransduction. Although we, like many other investigators, have found that the mesenteric receptors are anatomically very similar to those found in mammalian skin, few physiological characteristics of the mesenteric PCs and those of the skin have been compared. Action-potential rate-amplitude and frequency characteristics (10 Hz-1 KHz), as well as interval (IH) and peri-stimulus-time (PSTH) histograms in response to sinusoidal displacements were obtained from nerve fibers innervating mesenteric PCs and from PC fibers innervating cat glabrous skin. The intensity characteristics obtained on both preparations showed similar response profiles, including equal slopes for low stimulus intensities (approximately 10, with impulse ratios/20 dB displacement) and one and two impulse/cycle entrainment. The frequency characteristics of both groups were U-shaped with similar low-frequency slopes (-12.5 dB/octave) and bandwidths (Q(3dB) = 1.4). The best frequency for both the tactile PCs' and mesenteric PCs was 250 Hz, which is in the expected range. The IHs showed entrainment and the PSTHs showed neither transient responses nor adaptation to steady-state sinusoidal stimuli. The functional similarity between mesenteric PCs' nerve responses and those of tactile PC afferents, as well as the receptors' anatomical similarity, lead us to suggest that the mesenteric PC can act as a model for those in the skin. Furthermore, since the frequency characteristics of the two PC types are similar, it is concluded that the skin, while attenuating stimulus intensity, does not impart temporal filtering of vibratory stimuli.