|
|
||||||||
1Krieger Mind/Brain Institute, 2Department of Neuroscience, and 3Department of Biomedical Engineering, Johns Hopkins University, Baltimore, Maryland
Submitted 3 January 2005; accepted in final form 11 July 2005
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The time-course along which afferent adaptation and recovery operate is largely unknown, in part because the phenomenon itself has been difficult to study for reasons discussed in the companion paper (Bensmaïa et al. 2005
). Recording from RA afferents in monkeys and cats, Whitsel et al. (2000)
observed a slight drop in mean spike rate within the first few seconds of a vibratory stimulus. The observed decline in responsivity was accompanied by a slight increase in the phase angle of the entrained response within the first 100500 ms of the stimulus. Recording the responses from PC fibers in cat mesentery, O'Mara et al. (1988)
found that the spike rates recovered to their preadaptation levels along approximately exponential time-courses, with time constants <30 s. In this study, we investigated the time-course of adaptation in slowly adapting type I (SA1), RA, and pacinian (PC) fibers using two indices of afferent sensitivity: 1) absolute and entrainment thresholds measured using a novel tracking algorithm described in Bensmaïa et al. (2005)
and 2) the phase at which afferents produce a spike within each cycle of the adapting stimulus. We show that the observed rates of adaptation and recovery are compatible with predictions derived from a simple integrate-and-fire model. Finally, we compare the time-courses of afferent adaptation and recovery with previously reported time-courses of psychophysical adaptation and recovery and conclude that central factors must also play a role in the psychophysical phenomenon of vibratory adaptation.
| METHODS |
|---|
|
|
|---|
Threshold shifts
Absolute and entrainment thresholds of the afferents were measured on-line using a tracking algorithm, described in detail in the companion paper. The averages of the last five I0 and I1 estimates from the adaptation period were used to estimate the adapted thresholds. Similarly, the averages of the last five I0 and I1 measurements from the recovery period provided estimates of the recovered thresholds. The time-course of vibrotactile adaptation and recovery was characterized post hoc as a series of exponential functions. Preadaptation thresholds were fitted with a horizontal line, because thresholds theoretically remained constant during this period. Thresholds from the adaptation period were fitted with the method of least squares to an exponential function of the form
![]() | (1) |
is the threshold shift, ta is the time of onset of the adapting stimulus, and
a is the adaptation time constant. Similarly, recovering thresholds were fit to an exponential of the form
![]() | (2) |
is the degree of threshold recovery, and
r is the recovery time constant. (Note that
a,
r,
, and
are free parameters and
=
if thresholds recover to preadaptation levels.) As can be seen from plots of thresholds versus time (Figs. 1
|
|
|
Phase shift
Along with increases in threshold, adaptation also causes a progressive shift in the phase at which the afferent produces a spike within each stimulus cycle: as the sensitivity of the receptor declines, it takes longer for the receptor potential to reach threshold and the impulse phase progressively increases. Note that accurate estimates of phase shift require near-perfect entrainment of the neural response as impulse phase is also affected by the immediate spiking history of the afferent (Freeman and Johnson 1982b
). Specifically, the impulse occurs later in a stimulus cycle if the previous cycle evoked an impulse than if it did not (Freeman and Johnson 1982b
). If every stimulus cycle evokes the same number of impulses, systematic shifts in impulse phase are caused by changes in the excitability of the neuron, i.e., to adaptation, and thus can yield an additional estimate of the time course of adaptation.
We measured the phase adaptation, i.e., the progressive shift in impulse phase during the adaptation period, for stimulus runs on which afferent responses were entrained to the adapting stimulus (i.e., produced the same number of spikes on each stimulus cycle). The criterion for entrainment was that the afferent fire a given number of spikes (typically 1, but sometimes 2 or even 3) on 80% or more of adapting stimulus cycles over the course of the adaptation period. For each run meeting this criterion, the phase of the first spike produced on each cycle of the adapting stimulus was measured, along with its time of occurrence within the adapting period. The impulse phase was the time of the first spike produced in a given cycle, normalized by the period of the cycle (i.e., 1/Fa) and multiplied by 360°. The phases of the spikes evoked during the first few cycles of the stimulus in each adapting interval were not included in the analysis because impulse phase shifts rapidly before reaching steady state (i.e., a state changing over the span of seconds) for reasons unrelated to adaptation (but that are beyond the scope of this study). We show here that phase adaptation accompanies threshold shift and, like threshold shift, progresses along an exponential time-course. Thus time constants of adaptation could be independently estimated using the procedure described above then compared with those obtained using threshold shifts.
| RESULTS |
|---|
|
|
|---|
Time-course of the threshold adaptation
Figures 13 show time-courses of adaptation and recovery obtained from typical SA1, RA, and PC afferents. SA1 and RA thresholds shifted rapidly upon onset of the adapting stimulus and remained at their adapted levels throughout the remainder of the adapting period. On the other hand, PC thresholds often had not stabilized by the end of adaptation periods lasting 2, 4, and sometimes even 8 min. SA1 afferents tended to recover most rapidly, whereas RA and PC afferents recovered at similar, slower rates. Differences in the rates of recovery across fiber types were less pronounced than differences in the rates of adaptation.
PC fibers are exquisitely sensitive to mechanical stimulation, often exhibiting submicron absolute thresholds. Furthermore, PC fibers are relatively unsusceptible to adaptation (Bensmaïa et al. 2005
). As a result, the recovered thresholds obtained from these afferents were relatively noisy. More reliable estimates of the time-course of PC adaptation were obtained by estimating the time-course of the shift in impulse phase during the adaptation period.
Time-course of adaptation
The left panels of Fig. 4 show the distributions of threshold adaptation time constants for the three types of afferents. SA1 fibers tended to adapt most quickly, with mean time constants of 11 and 9 s for I0 and I1. RA fibers had somewhat longer adaptation time constants (
17.9 s;
: 13.2 s). PC fibers were the slowest to adapt (
: 41.8 s;
: 28.1 s; Fig. 4, left). SA1 I0 and I1 adaptation time constants were not significantly different (paired t-test, P > 0.05). RA and PC I0 adaptation time constants were 50% longer than their I1 counterparts, but that difference was only significant for RA afferents (paired t-test, P < 0.01); the difference was not significant for PC fibers, likely because of the high degree of variability in the time constants obtained from these afferents.
|
The right panels of Fig. 4 show the distributions of recovery time constants for the three types of afferents. SA1 fibers recovered most rapidly, with mean time constants of 9.4 and 17.3 s for I0 and I1 thresholds, respectively. Recovery time constants for RA fibers (
I0: 27 s;
I1: 26.6 s), and PC afferents (
I0: 20.9 s;
I1: 29.6 s) were somewhat slower (Fig. 4, right). I0 and I1 recovery time constants were not significantly different from one another for both RA and PC fibers (paired t-test, P > 0.05). SA1 I0 recovery time constants, on the other hand, were almost twice as fast as I1 recovery time constants (paired t-test, P < 0.001).
Entrainment thresholds measured at the end of the recovery periods did not differ systematically from their preadaptation levels (Fig. 5). Preadaptation and recovered absolute thresholds exhibited greater variability than their I1 counterparts, likely because of the fact that I0 tended to be much smaller than I1 and thus more susceptible to measurement noise. Overall, these results suggest that all three afferent types generally returned to their resting responsiveness within the allotted recovery periods (which ranged from 2 to 16 min for SA1 and RA fibers and from 4 to 16 min for PC fibers).
|
SA1 I1 recovery was slower than I1 adaptation (t-test, P < 0.001), but the I0 rates were not significantly different from one another. RA recovery tended to be slower than adaptation (Fig. 6) for both I0 and I1 (t-test, P < 0.001). Finally, PC adaptation and recovery time constants were not significantly different for both I0 and I1 (again, perhaps because of the variability of the PC time constants).
|
In the previous analyses, we assumed adaptation and recovery time constants to be independent of the stimulus parameters. To test this assumption, we performed an ANOVA on each type of time constant (adaptation and recovery, absolute, and entrainment threshold) with Ft, Fa, and Aa as covariates (with an additional factor to control for differences across fibers of a given type). Time constants were found to be largely independent of stimulus parameters, i.e., the effects of test and adapting frequency and adapting amplitude were generally nonsignificant (
= 0.05). In the few cases in which there was a significant effect of a stimulus parameter on
, the effect accounted for a negligible fraction of the variance in the time constants (
p2
0.1). These results suggest that the time constants are an inherent property of the afferent fibers themselves and are not stimulus-dependent. Given the inconsequentiality of these effects, they are not reported in detail here.
Time-course of phase adaptation
On stimulus runs in which the afferent response was entrained with the adapting stimulus, each cycle of the adapting stimulus contributed a data point to the estimate of the impulse phase adaptation time constant. Figure 7 shows the time-course of phase adaptation for one receptor of each type. For all three afferent types, impulse phase increased exponentially with time. Figure 8 shows the cumulative histograms of phase adaptation time constants for the three fiber types. SA phase adaptation time constants were faster than their RA and PC counterparts (t = 2.95 and 3.98 for RA and PC fibers, respectively; P < 0.01). On the other hand, RA and PC phase adaptation time constants were not significantly different (t = 1.77, P > 0.05). Furthermore, the variability of these time constant estimates was much lower than that of their threshold shift counterparts, likely because of the fact that impulse phase could be estimated on each cycle of the adapting stimulus, whereas only one estimate of threshold was obtained on each adaptation trial. The median number of observations used to compute phase adaptation time constants was
4,500 (mean,
7,000), two orders of magnitude higher than the median number of adaptation trials per run (
30).
|
|
| DISCUSSION |
|---|
|
|
|---|
The time-course of adaptation was measured using two methods: one using thresholds and the other using impulse phase in entrained responses. We found the rates of adaptation measured using the two methods to be different, particularly for PC afferents (Table 1). Specifically, phase adaptation time constants were faster than their threshold shift counterparts for PC fibers. The discrepancies between the time constants measured using the two methods could be caused by measurement error or by differences in the mechanisms underlying the measured quantities.
|
T is the residual increment in spiking threshold on a stimulus cycle after a cycle on which a spike was produced (reflecting the relative refractoriness of the spiking mechanism). When the afferent response is entrained with the stimulus, the impulse phase is proportional to the instantaneous threshold of the afferent and inversely proportional to the transducer gain. Thus
![]() | (3) |
![]() | (4) |
On the other hand, PC threshold adaptation time constants were substantially longer than their phase adaptation counterparts. This discrepancy is likely because of measurement error in the threshold shift time constants. Indeed, PC thresholds were low, sometimes submicron, as were threshold shifts. Measurements of instantaneous PC thresholds were thus error prone. (Note that estimates of steady-state threshold shift were much less susceptible to measurement error because these rely on five instantaneous estimates of threshold.) In contrast, the time-course of PC phase adaptation can be reliably measured because each adaptation run yielded many estimates of the impulse phase, approximately one for each cycle of the adapting stimulus. Furthermore, PC afferents produce very regular entrained responses to high-frequency stimuli (Freeman and Johnson 1982b
) so there was little variation (
0.10.2 ms) in impulse phase over small time intervals (Fig. 7). In summary, the time-course of phase adaptation for RA and SA1 afferents was comparable with that of I1 adaptation, as predicted by a simple integrate-and-fire model. We propose that the discrepancy in the time constants obtained for PC threshold and phase adaptation is caused by measurement error in the former.
Comparing psychophysics and neurophysiology
Time-course of adaptation.
The time-course of vibrotactile adaptation has been studied in a number of psychophysical studies, either by measuring absolute thresholds at various intervals after the onset of the adapting stimulus (Hahn 1968a
,b
; Hollins et al. 1990
, 1991
) or by tracking the subjective intensity of a suprathreshold stimulus over time (Berglund and Berglund 1970
; Hahn 1968b
). Because psychophysical absolute thresholds depend on a small group of peripheral fibers, perhaps even a single one (Johansson and Vallbo 1979
), the time-course of adaptation and recovery of (psychophysical) absolute thresholds is more comparable with that of afferent thresholds than is the time-course of adaptation and recovery of sensation magnitude, which is mediated by a larger population of peripheral fibers. To compare the time-courses of psychophysical and peripheral adaptation, we fitted exponential functions (analogous to that shown in Eq. 1) to the time-courses of psychophysical adaptation obtained from a pair of adaptation studies (Hollins et al. 1990
, 1991
).
Psychophysical adaptation time constants for adapting frequencies of 10 and 50 Hz varied from
62 to 212 s (Hollins et al. 1990
, 1991
). Given that the RA channel is more sensitive than the SA1 channel at all but the lowest frequencies (< approximately 5 Hz; Bolanowski et al. 1988
), and that test stimulus frequencies fell within the most sensitive portion of the RA range (10 and 50 Hz), it is the time-course of adaptation within this sensory channel that was likely measured in the two threshold-tracking studies, although the PC channel may mediate thresholds at 50 Hz (Hollins et al. 1990
, 1991
).
In this study, RA adaptation time constants and their PC counterparts estimated from phase data ranged from 1 s to about 40 s (Figs. 4 and 8). If vibrotactile adaptation were exclusively a peripheral phenomenon, psychophysical and afferent thresholds would adapt within a similar time frame. However, the time-course of psychophysical adaptation was much longer than that predicted from our neurophysiological measurements.
Time-course of recovery from adaptation.
To quantify the rate of recovery from psychophysical adaptation, we fitted exponential functions (analogous to that shown in Eq. 2) to psychophysical (absolute) thresholds as they recovered after the termination of the adapting stimulus. These psychophysical recovery time constants ranged from about 26 to 316 s (Gescheider and Wright 1968
; Hahn 1968b
; Hollins et al. 1990
, 1991
). Recovery time constants derived from RA fibers ranged from 2 to 45 s. PC recovery time constants ranged from 18 to 92 s. Thus RA and PC time constants were somewhat faster than their psychophysical counterparts but fell within an overlapping range. That psychophysical adaptation and recovery operate on a slower time scale than their counterparts at the sensory periphery suggests central factors contribute to the psychophysical phenomenon.
Neural mechanisms
In the companion paper (Bensmaïa et al. 2005
), we propose that adaptation is produced by an influx of ions, probably calcium, through mechanosensitive channels. The increase in concentration of this ion in turn causes a change in the spiking threshold of the afferent or possibly in the transducer gain of the receptor. Because of its dependence on the electrochemical gradient across the receptor membrane, the magnitude of the change in the intracellular concentration of this ion is in part determined by its instantaneous intracellular concentration. The greater the intracellular concentration, the lower the driving force on the ion, and the lower the change in intracellular concentration, and vice versa. The dependence of the change in concentration on the instantaneous concentration is the type of relation that gives rise to an exponential time-course. Thus the ionic hypothesis of adaptation is consistent with its time-course.
The question remains why the time-course of adaptation for entrainment thresholds is different from its absolute threshold counterpart for RA afferents. Similarly, it is unclear why the time-course of I1 recovery is faster than that of I0 recovery for SA1 fibers. Within the simple integrate-and-fire framework sketched above, the absolute threshold is proportional to the spiking threshold and inversely proportional to the transducer gain (Bensmaïa et al. 2005
)
![]() | (5) |
T (when the afferent has reached steady state) is equal to its preadaptation level in RA afferents (Bensmaïa et al. 2005
and
is that the maintenance of
T as T changes is a process with a faster time-course than the shift in T itself. Similarly, differences in the rates of I0 and I1 recovery in SA1 afferents may be caused by differences in the rates of recovery of T and
T. To evaluate these hypotheses, a more detailed biophysical model of mechanoreceptive transduction is required. What is adaptation?
It is paradoxical that SA1 fibers, named for their slowly declining firing rates in response to a constant indentation, adapt most rapidly to a vibratory stimulus, whereas their rapidly adapting counterparts, RA and PC fibers, adapt to the same stimulus more slowly. Although they unfortunately share the same name, these two kinds of adaptation clearly rely on very different mechanisms. Rapid adaptation to constant indentation is measured in milliseconds and is, at least in the case of PC fibers, determined to a large extent by the mechanical properties of the Pacinian corpuscle (Loewenstein and Skalak 1966
). Adaptation to constant vibratory stimulation in the three types of fibers is measured in tens of seconds, which points to a completely different mechanism. For SA1 fibers, vibratory adaptation may be a manifestation of the same mechanism that underlies the slow decline in the response to a constant indentation that inspired their designation; in both cases, ion flow into the receptor may cause the receptor to become desensitized. According to this hypothesis, the time-course of the two types of adaptation should be identical. Figure 9 shows the firing rate evoked in an SA1 afferent by a 1-mm indentation lasting 30 s. The spike rate declines along an approximately exponential time-course with a time constant of 8.4 s, matching its vibratory adaptation counterpart. We recorded the responses of two other SA1 afferents to a step indentation and observed similar time courses (data not shown). We tentatively conclude that, for SA1 afferents, vibratory adaptation and adaptation to a constant indentation reflect the same underlying phenomenon, although more data will be needed to establish this conclusively.
|
| GRANTS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
Address for reprint requests and other correspondence: S. J. Bensmaïa, 3400 N. Charles St., Krieger 338, Baltimore, MD 21218 (E-mail: sliman{at}jhu.edu)
| REFERENCES |
|---|
|
|
|---|
Chubbuck JG. Small motion biological stimulator. Johns Hopkins APL Tech Digest 5: 1823, 1966.
Freeman AW and Johnson KO. A model accounting for effects of vibratory amplitude on responses of cutaneous mechanoreceptors in macaque monkey. J Physiol 323: 4364, 1982a.
Freeman AW and Johnson KO. Cutaneous mechanoreceptors in macaque monkey: temporal discharge patterns evoked by vibration, and a receptor model. J Physiol 323: 2141, 1982b.
Gescheider GA and Wright JH. Effects of sensory adaptation on the form of the psychophysical magnitude function for cutaneous vibration. J Exp Psychol 77: 308313, 1968.[CrossRef][ISI][Medline]
Hahn JF. Low-frequency vibrotactile adaptation. J Exp Psychol 78: 655659, 1968a.
Hahn JF. Tactile adaptation. In: The Skin Senses, edited by Kenshalo DR. Springfield, IL, Charles C Thomas, 1968b, p 322330.
Hollins M, Delemos KA, and Goble AK. Vibrotactile adaptation on the face. Percept Psychophys 49: 2130, 1991.[ISI][Medline]
Hollins M, Goble AK, Whitsel BL, and Tommerdahl M. Time course and action spectrum of vibrotactile adaptation. Somatosens Mot Res 7: 205221, 1990.[ISI][Medline]
Johansson RS and Vallbo ÅB Detection of tactile stimuli. Thresholds of afferent units related to psychophysical thresholds in the human hand. J Physiol 297: 405422, 1979.
Loewenstein WR and Skalak R. Mechanical transmission in a Pacinian corpuscle. An analysis and a theory. J Physiol 182: 346378, 1966.
O'Mara S, Rowe MJ, and Tarvin RP. Neural mechanisms in vibrotactile adaptation. J Neurophysiol 59: 607622, 1988.
Talbot WH, Darian-Smith I, Kornhuber HH, and Mountcastle VB. The sense of flutter-vibration: comparison of the human capacity with response patterns of mechanoreceptive afferents from the monkey hand. J Neurophysiol 31: 301334, 1968.
Whitsel BL, Kelly EF, Delemos KA, Xu M, and Quibrera PM. Stability of rapidly adapting afferent entrainment vs responsivity. Somatosens Mot Res 17: 1331, 2000.[CrossRef][ISI][Medline]
Bensmaïa SJ, Leung YY, Hsiao SS, and Johnson KO.Vibratory adaptation of cutaneous mechanoreceptive afferents. J Neurophysiol 94: 30243037, 2005.
Bolanowski SJ, Gescheider GA, Verrillo RT, and Checkosky CM. Four channels mediate the mechanical aspects of touch. J Acoust Soc Am 64: 16801694, 1988.
This article has been cited by other articles:
![]() |
M. A. Muniak, S. Ray, S. S. Hsiao, J. F. Dammann, and S. J. Bensmaia The Neural Coding of Stimulus Intensity: Linking the Population Response of Mechanoreceptive Afferents with Psychophysical Behavior J. Neurosci., October 24, 2007; 27(43): 11687 - 11699. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. S. Beauchamp, N. E. Yasar, N. Kishan, and T. Ro Human MST But Not MT Responds to Tactile Stimulation J. Neurosci., August 1, 2007; 27(31): 8261 - 8267. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. J. Berryman, J. M. Yau, and S. S. Hsiao Representation of Object Size in the Somatosensory System J Neurophysiol, July 1, 2006; 96(1): 27 - 39. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Hollins, F. Lorenz, and D. Harper Somatosensory coding of roughness: the effect of texture adaptation in direct and indirect touch. J. Neurosci., May 17, 2006; 26(20): 5582 - 5588. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. J. Bensmaia, J. C. Craig, and K. O. Johnson Temporal Factors in Tactile Spatial Acuity: Evidence for RA Interference in Fine Spatial Processing J Neurophysiol, March 1, 2006; 95(3): 1783 - 1791. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. W. Goodwin Paradoxes in Tactile Adaptation. Focus on "Vibratory Adaptation in Cutaneous Mechanoreceptive Afferents" and "Time-Course of Vibratory Adaptation and Recovery in Cutaneous Mechanoreceptive Afferents" J Neurophysiol, November 1, 2005; 94(5): 2995 - 2996. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |