Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Temporal patterns in electrical nerve stimulation: Burst gap code shapes tactile frequency perception

  • Kevin K. W. Ng ,

    Roles Conceptualization, Formal analysis, Investigation, Visualization, Writing – original draft, Writing – review & editing

    kevin.ng@unsw.edu.au

    Affiliations School of Medical Sciences, UNSW Sydney, Sydney, Australia, Neuroscience Research Australia, Sydney, Australia

  • Christoffer Olausson,

    Roles Formal analysis, Investigation, Visualization, Writing – review & editing

    Affiliations Neuroscience Research Australia, Sydney, Australia, Faculty of Medicine and Health Sciences, Linköping University, Linköping, Sweden

  • Richard M. Vickery ,

    Roles Conceptualization, Funding acquisition, Resources, Supervision, Writing – review & editing

    ‡ These authors are joint senior authors on this work.

    Affiliations School of Medical Sciences, UNSW Sydney, Sydney, Australia, Neuroscience Research Australia, Sydney, Australia

  • Ingvars Birznieks

    Roles Conceptualization, Funding acquisition, Resources, Supervision, Writing – review & editing

    ‡ These authors are joint senior authors on this work.

    Affiliations School of Medical Sciences, UNSW Sydney, Sydney, Australia, Neuroscience Research Australia, Sydney, Australia

Abstract

We have previously described a novel temporal encoding mechanism in the somatosensory system, where mechanical pulses grouped into periodic bursts create a perceived tactile frequency based on the duration of the silent gap between bursts, rather than the mean rate or the periodicity. This coding strategy may offer new opportunities for transmitting information to the brain using various sensory neural prostheses and haptic interfaces. However, it was not known whether the same coding mechanisms apply when using electrical stimulation, which recruits a different spectrum of afferents. Here, we demonstrate that the predictions of the burst gap coding model for frequency perception apply to burst stimuli delivered with electrical pulses, re-emphasising the importance of the temporal structure of spike patterns in neural processing and perception of tactile stimuli. Reciprocally, the electrical stimulation data confirm that the results observed with mechanical stimulation do indeed depend on neural processing mechanisms in the central nervous system, and are not due to skin mechanical factors and resulting patterns of afferent activation.

Introduction

Temporal features of neural spike activity play a major role in encoding tactile information [14]. In a previous study, we used pulsatile mechanical stimulation to investigate the neural basis of frequency perception in the flutter range. When we grouped pulsatile stimulation patterns into periodic bursts, we found that perceived frequency was determined by the duration of the silent gap between bursts of spikes in peripheral afferents (termed the burst gap code), irrespective of the number of spikes within a burst or the burst rate [5]. This stands in contrast to the mean spike rate or periodicity code that are often proposed to be neural codes for frequency perception in higher neural centres [69]. Electrical burst stimulation has been suggested to be useful in restoring the sense of touch in amputees using a prosthesis [1012] or in the treatment of movement disorders through deep brain stimulation [1315]. These studies were generally concerned with the efficacy of inducing neural activity, rather than the effect of burst stimulation on information encoding. We believe that careful consideration of neural coding related to the temporal aspects of these burst patterns could potentially improve the performance of such stimulation strategies.

The discovery of the burst gap code for perceived frequency depended on the use of brief mechanical pulses that each only evoke a single spike in responding afferents, as had been verified using microneurography, allowing us to precisely control the spiking pattern of responding tactile afferents [5, 16, 17]. Electrical stimulation activates peripheral axons in a manner different than that of mechanical stimulation [18], by bypassing the specialised mechanoreceptors in the skin and directly stimulating the tactile afferent axons [19, 20]. Electrocutaneous stimulation non-selectively activates the different types of tactile afferents (associated with different skin receptors) [21, 22], which contrasts with mechanical stimulation where the type of afferent recruited varies depending on the sensitivity of the fibre type to stimulus characteristics such as the vibration frequency of a mechanical sinusoid [23, 24].

A single tap on a finger generates a complex propagating wave on the fingertip skin [25, 26] that will evoke complex and potentially variable population spike patterns due to sporadically responding afferents [27]. In particular, such a stimulus would likely reliably activate afferents whose receptive fields are directly under the probe, but this activation would become stochastic at some distance from the locus of stimulation [24, 28]. This raises the possibility that the burst gap code that we observed with mechanical stimulation is in some way dependent on these spatial activation patterns and/or predominant activation of certain afferent types. As an alternative, electrical stimulation is a robust method to reliably activate axons, and a simple technique to implement in neuroprosthetic interfaces [29] or haptic displays [30]. Thus, in this study we aimed to use electrical stimuli to directly evoke spike trains in tactile afferent fibres while avoiding complex mechanical skin phenomena and establish the burst gap model for frequency perception.

Materials and methods

Subjects

Two separate experiments were conducted. Twelve subjects (aged 21–26, 5 females) participated in Experiment 1. Eight subjects (aged 21–25, 5 females) participated in Experiment 2. Three of the subjects participated in both experiments. All subjects were healthy and without any known history of altered tactile function. The experimental protocols were approved by the Human Research Ethics Committee (HC16245) of UNSW Sydney and written consent was obtained from subjects before experimentation began.

Apparatus

Stimulus patterns were generated via a CED Power1401 mk II data acquisition system using Spike2 (Cambridge Electronic Design, Cambridge, UK) and MATLAB (Mathworks, Natick, MA, USA) software. These patterns triggered constant current electrical pulses from a DS7AH stimulator (Digitimer, Welwyn Garden City, UK). The electrical stimuli were delivered to the right index finger of subjects using Kendall 200 Series foam electrodes (Covidien, Mansfield, MA, USA), with one electrode placed on the distal interphalangeal joint and another on the proximal phalanx, to stimulate the digital nerve. A button box was used by subjects to indicate their responses to the psychophysical task. Compound action potentials were recorded using a PowerLab data acquisition system and LabChart software (ADInstruments, Sydney, Australia).

Stimulation patterns

Experiment 1.

The stimulus timing patterns used were based on that of our previous study [5]. In our first experiment, the aim was to determine whether burst gap coding would best describe the perception of frequency in subjects. The test stimuli consisted of one of four trains with periodic bursts of 2–4 pulses (Fig 1A), such as would be generated in the afferent response to high-amplitude vibration [24, 31]. The individual pulses within the burst were spaced either 4.35 ms (stimuli a-c) or 8.7 ms (stimulus d) apart. The period between repeating bursts was fixed at 43.5 ms, resulting in a burst rate of 23 Hz. However, the mean pulse rate and the inter-burst interval (tL) varied noticeably between the four stimuli, ranging from 46 to 92 pulses/s and 26.1–39.15 ms, respectively.

thumbnail
Fig 1. Schematic of stimulus patterns and experimental protocol, each red vertical line represents an electrical pulse.

(A) In Experiment 1, test stimuli were pulse trains with periodic bursts of 2–4 pulses. The stimuli have the same periodicity of 43.5 ms (burst rate of 23 Hz), but different inter-burst intervals (tL) and overall pulse rates. (B) In Experiment 2, test stimuli consisted of a short interval doublet (tS), followed by an 87 ms long inter-pulse interval (tL). The short inter-pulse interval was varied in steps of 4.35 ms, from 4.35 ms up to 56.55 ms. (C) Two interval forced choice procedure used to determine perceived frequency of test stimuli.

https://doi.org/10.1371/journal.pone.0237440.g001

Experiment 2.

In our second experiment, we measured the critical value of the inter-pulse interval within a burst that determined how the pulses contribute to perceived frequency. The test stimuli were one of 13 trains of a 2-pulse burst (doublet) with a short interval (tS) and a long interval (tL) (Fig 1B). The short inter-pulse interval (tS) ranged from 4.35 ms to 56.55 ms, increasing in steps of 4.35 ms. The long interval (tL) was fixed at 87 ms. The burst gap code predicts that when two pulses are close enough in time, they will be treated as a burst and accordingly, the perceived frequency will be determined by the burst gap (approximately 11.5 Hz in this experiment).

Electrical stimulation parameters

Electrical stimuli were 0.1 ms wide square pulses. Prior to the start of each experiment, the electrical current level was optimised for the subject so that the electrical pulses were perceived clearly as a series of taps without discomfort. Current levels used ranged from 6.5 to 11.5 mA across subjects and were kept at a fixed level within a session.

Verifying the stable activation of a population of afferents

To verify that the afferent axons fired reliably to each electrical pulse within a burst in our closely spaced stimulation patterns, sensory nerve action potentials (SNAPs) were measured in two subjects by stimulating the sensory nerves distally (as described above) and recording the compound action potential more proximally [32]. The recording electrodes were located just proximal to the wrist crease over the median nerve, with the ground electrode placed on the thenar eminence. SNAPs were recorded from 200 repetitions of the burst patterns used in Experiment 1 (Fig 1A).

Measuring the point of subjective equality for frequency

To determine the perceived frequency of each test stimulus, a two-interval forced choice paradigm was used as in our previous study [5]. Test stimuli were compared against 6–8 different comparison stimuli, which were regular trains of pulses with repetition rates ranging from 14–54 Hz for Experiment 1 and 7–21 Hz for Experiment 2. For each trial, a test stimulus and a comparison stimulus were presented for 1 s each, in a random order and separated by 500 ms (Fig 1C). The subject indicated by button press which of the two stimuli had the higher frequency. Following optimisation of the electrical current, a practice round was conducted which comprised a selection of test and comparison stimuli to familiarise subjects with the task, before the actual test procedure began. To obtain psychometric curves, each test condition was compared 20 times against each of the regular comparison stimuli, giving rise to 120–160 trials per test stimulus.

At each comparison frequency, we determined the proportion of times the subject indicated that the comparison stimulus had a higher frequency than the test stimulus (PH). To obtain a linear psychometric function, the logit transformation ln(PH/(1−PH)) was applied to the data [33]. The value of the comparison frequency that is equally likely to be judged as higher or lower than the test stimulus was taken as the point of subjective equality (PSE). This value corresponds to the frequency where the zero crossing of the logit axis occurs for the regression line fitted to the logit transformed data (example in S1 Fig).

Statistical analysis

For Experiment 1, linear regression was performed to determine if the slope of the line fitted to the data would be better explained by the burst gap (i.e. the inter-burst interval) or the burst rate (periodicity), which was constant at 23 Hz. A two-way ANOVA (repeated measures within the same subject for the different stimulation methods) was used to determine if there were differences between the data collected here using electrical stimulation and the mechanical data collected in our previous study. For Experiment 2, a joinpoint regression model created using Joinpoint Trend Analysis Software (National Cancer Institute, Bethesda, MD, USA) was used to determine the trend that best fitted the psychophysical data. An alpha level of 0.05 was used for all statistical tests.

Results

Stable activation of the afferent population

Electrical stimulation with patterns used in Experiment 1 (Fig 1A) resulted in SNAPs that were consistent within a burst irrespective of inter-pulse intervals used or the number of pulses within a burst (Fig 2, Fig 3A). To examine pulse-to-pulse variability within individual bursts, we normalised the amplitude of subsequent SNAPs to that of the first SNAP for each burst; this showed that the nerve response to pulses occurring later in a burst were not negatively affected by preceding pulses (Fig 3B), which suggests that the stimulus burst patterns were faithfully reproduced in the nerve firing patterns.

thumbnail
Fig 2. Sensory Nerve Action Potential (SNAP) recordings.

Average of 200 compound action potentials recorded from the median nerve in one subject, in response to Experiment 1 stimulation patterns illustrated in Fig 1A. The red circles indicate SNAPs, while the clipped waveform occurring before each SNAP represents the electrical stimulation artefact.

https://doi.org/10.1371/journal.pone.0237440.g002

thumbnail
Fig 3. SNAP amplitudes within a burst.

(A) The mean peak-to-peak amplitude (n = 200) of evoked SNAPs in one subject. Error bars denote standard deviation. (B) Geometric mean of SNAP amplitudes after normalising to the first SNAP within each burst response. Error bars denote 95% confidence intervals.

https://doi.org/10.1371/journal.pone.0237440.g003

How is perceived frequency encoded?

In Experiment 1, the aim was to determine which coding mechanism could best describe the perception of frequency in subjects when the electrical pulses were grouped into bursts. The mean data from the 12 subjects are shown in Fig 4. The observed perceived frequencies were 2- to 4-fold lower than that predicted by the mean pulse rate (Fig 4A green triangles), indicating that a simple rate code could not explain perceived frequency in subjects. Linear regression was performed between the reciprocal of the burst gap (tL) and the perceived frequency for each stimulus (Fig 4B). The regression line had a slope that was non-zero (1.2, 95% CI 1.0–1.4; R2 = 0.76; p < 0.0001), which does not correspond to that of the burst rate / periodicity (expected slope 0, Fig 4A pink triangles), but is a good match to the reciprocal of the burst gap (tL). The slope here closely resembles that of the mechanical data from our previous study (slope = 1.3, 95% CI 1.0–1.5; R2 = 0.68; p < 0.0001), which also best matched 1/tL [5].

thumbnail
Fig 4. Possible coding schemes for frequency perception of electrical burst stimuli.

(A) Boxplots with whiskers represent the subjective perceived frequency of the electrical stimuli used in Experiment 1 (n = 12). Green and pink arrowheads indicate predicted PSE values if the frequency judgement had been based on the mean pulse rate and the burst rate (periodicity), respectively. (B) Fitted regression line (black line) indicating that the mean electrical PSE data (black circles; n = 12) for perceived frequency is a close match to the reciprocal of the burst gap tL. The regression line fitted to mechanical data from our previous study [5] is shown for comparison (orange line). The dotted line is the line of identity, where the PSE would equal the reciprocal of the burst gap tL. Error bars denote 95% confidence intervals.

https://doi.org/10.1371/journal.pone.0237440.g004

A two-way ANOVA indicated that there was a statistically significant main effect by the stimulus pattern (p < 0.0001, F2,44 = 143, 71% of total variation) and the method of stimulation (p = 0.0084, F1,22 = 8.4, 4.7% of total variation). However, there was no interaction effect (p = 0.36, F3,66 = 1.1), suggesting that varying the stimulation pattern had the same effect on frequency perception regardless of whether mechanical or electrical stimulation was used.

What defines an inter-pulse interval as part of a burst?

For Experiment 2, we varied the duration of the short interval to determine which inter-pulse intervals (tS) would be treated as a burst, so that the perceived frequency depended only on the longer interval (tL). The results from the 8 subjects are shown in Fig 5 which plots the apparent frequency of each stimulus. A joinpoint regression model was fitted, which indicated that the data was best explained by three linear segments. The first segment extended from 4.35 to 13.05 ms, with a slope that was not significantly different from zero (p = 0.64). The next two segments had slopes that were significantly non-zero. The second segment spanned 13.05 to 34.8 ms and had a slope of 0.2 (p = 0.0002). The last segment had a range between 34.8 and 56.55 ms, with a slope of −0.08 (p = 0.0002).

thumbnail
Fig 5. Relationship between the short inter-pulse interval and perceived frequency.

PSE values obtained using Experiment 2 stimuli, indicated by black circles (Electrical PSE; n = 8). The black line signifies the three fitted linear segments as determined by the joinpoint regression model. The green and blue lines represent predicted PSE values if frequency judgement had been based on the mean pulse rate and the burst gap code, respectively. The orange line (Mechanical PSE) is the joinpoint model based on the mechanical data from our previous study [5]. Error bars denote 95% confidence intervals.

https://doi.org/10.1371/journal.pone.0237440.g005

The slope for the segments here were similar to those reported in our previous study involving mechanical stimulation (Fig 5, Mechanical PSE), though the joinpoint between the second and third segments occurred at 30.45 ms in that study [5], thus at a slightly shorter inter-pulse interval than in the current study (34.8 ms). Nevertheless, the same basic findings were observed in both studies, as shown by the confidence intervals overlapping with the joinpoint model based on the mechanical data and an identical first jointpoint. This suggests that the same definition of a burst can be applied to both mechanical and electrical stimulation in the flutter range.

Discussion

The main objective of the present study was to verify whether the burst gap coding model for frequency perception could be observed with electrical stimulation in the flutter range. This is an important step to demonstrate due to the potential value of this coding strategy for implementation in bionic systems for providing somatosensory information. Our results show that the burst gap model, based on the duration of the silent inter-burst interval, could reliably predict perceived frequency of complex bursting spike patterns when the spikes were elicited by electrical pulse stimulation of axons. This means that the same overall spike rate in an axon can produce a different frequency percept, depending on how those spikes are arranged in time. This reinforces the importance of the temporal structure of spiking patterns in neural processing and sensory perception.

Encoding perceived frequency

Electrical stimulation bypasses the complex viscoelastic properties of the skin that may have affected previous work using mechanical stimulation, and so provides strong independent confirmation of the burst gap code. We verified that the efficiency of the electrical stimulus was not affected by the repetition rate by measuring SNAP recordings, which showed that even the shortest inter-pulse interval (4.35 ms) resulted in reliable and consistent activation of tactile afferents. As we used a reasonably strong stimulation current at a supra-threshold level, it was expected that a majority, if not all, of the large myelinated afferents of the digital nerve would be reliably activated regardless of the inter-pulse interval chosen. The ability for the peripheral sensory nerve to faithfully relay the presented burst stimulation patterns, as evident in the SNAP recordings, suggests that perceived frequency depends on the processing of these signals within the central nervous system and not peripheral factors.

While the results in our first experiment, where we used electrical pulses to create bursting trains of the same periodicity but different number of spikes, were best explained by the burst gap code rather than a periodicity or rate code, the mean PSEs were lower than that would be predicted simply from the inter-burst duration. Though the difference in PSEs between electrical and mechanical stimulation was statistically significant, it is unknown whether this difference indicates any physiological importance, as it is well within the Weber fraction of ~0.2–0.3 that has been previously reported in the literature [16, 3438]. One possible explanation for this bias towards lower frequencies may be that electrical stimulation recruits all afferent types non-selectively, so that the slowly adapting (SA) afferents would be recruited in addition to the fast adapting (FA) afferents, whereas the mechanical stimulus predominantly activated FA afferents [5]. Although previous research has found that SA afferents and their corresponding cortical sensory neurons do not contribute to vibrotactile frequency perception [3941], their activation and the corresponding sense of pressure they convey [23, 39] may bias subject judgements in psychophysics tests.

These findings are consistent with recent findings from our laboratory which revealed that the spiking pattern, rather than activated afferent type, determined the perceived frequency of repetitive mechanical stimuli [16]. In that study, we found that low-frequency spike trains in Pacinian (FA2) afferents could induce a vibratory percept, even though Meissner (FA1) afferents had been traditionally thought to be solely responsible for perception of low frequency flutter vibrations (<60 Hz), whereas higher frequency activation of Pacinian (FA2) afferents would evoke a different percept of vibratory hum [24]. A generalised frequency processing system that operates with inputs from a variety of receptor type inputs would explain these findings, as well as the constancy of vibrotactile frequency perception across different skin regions innervated by different afferent types [34, 42], e.g. hairy skin lacks Meissner corpuscles [43, 44]. The congruence of the present results using electrical stimulation whereby all types of afferents are recruited non-selectively, with the previous results using mechanical stimuli that activated predominantly FA afferents, supports the suggestion of a generalised frequency processing mechanism.

Natural stimuli encountered during object manipulation and surface exploration generally evoke complex discharge patterns in tactile afferents [2, 4547]. Additionally, numerous studies have implicated the importance of frequency content in skin vibrations for encoding tactile information such as texture [45, 48, 49]. Accordingly, the frequency perception mechanism substantiated here may represent a more universal method of encoding complex naturalistic patterns of vibration, which may not necessarily have a simple fixed periodicity. For instance, we previously showed that this sensory coding strategy also applied to stimuli that are irregular or aperiodic in nature, as it relies on the processing of individual inter-pulse intervals by weighting their perceptual contribution to perceived frequency depending on their duration, rather than detection of periodicity or spike counting [50].

The value of bursts

The second experiment, with doublets of different inter-pulse durations, clearly indicate that when the inter-pulse interval for the doublet is sufficiently brief (<15 ms), the pulses are interpreted as belonging to a single burst event and the burst gap determines the perceived frequency. As the separation of the two pulses is increased beyond 15 ms, the propensity for two consecutive pulses to be regarded as a burst decreases and the within-doublet interval appears to have a partial weighted contribution to frequency perception [50]. Beyond 30–35 ms, each individual pulse is treated as a discrete event and the perceived frequency approaches the mean pulse rate.

Burst firing is a neural feature common in various sensory systems [51]. It has been proposed that bursts play a major role in the reliable signalling of neural information, due to its ability to induce long-term synaptic plasticity and to encode more information as compared to single spikes [52, 53]. Moreover, bursting inputs may be required to cause the firing of postsynaptic neurons through interactions at the preferred oscillation or resonance frequency of the particular synapse and may represent a way to increase transmission security of weak inputs [54]. The burst duration reported in such studies ranges between 10–25 ms, which agrees with the 15 ms value reported in our current study.

Burst stimulation has been used as a strategy for delivering sensory information in studies investigating tactile displays [55, 56] and neural prosthetics [1012]. Using stimulus bursts may provide additional parameters that can be varied in controlling the perception of tactile features such as intensity sensation. Previous studies have varied the width of the electrical pulse, which varies the charge delivered and changes the number of afferents recruited, in turn modulating perceived intensity [37]. However, altering the charge delivered by a pulse risks generating discomfort in the subject, so a method which can extend the range of intensity sensation without changing individual electrical impulses may be of value [57]. Whilst modulating the pulse rate at a fixed pulse width can influence perceived intensity [37], it would also change perceived frequency [57]. The benefit of the burst stimulation method we describe here is that we can control perceived frequency, which is determined by the burst gap, regardless of burst features such as the number of pulses within a burst [5]. Hence, in haptic devices and brain-machine interfaces, this coding strategy may provide a way to modulate frequency perception by means of changing the inter-burst interval and possibly other perceptual features of tactile stimuli encoded by spiking activity within individual bursts [58].

The present study used external electrodes to stimulate the digital nerve and confirm the burst gap phenomenon for frequency encoding, albeit with a small frequency shift compared with mechanical activation. The use of intrafascicular stimulation with multi-electrode arrays could activate nerve fascicles which may have a more restricted range of afferent types, in addition to offering improved spatial resolution [29, 59]. It would be interesting to see whether predominantly FA-containing fascicles show frequency responses that are a closer match to the results with mechanical stimulation in our previous study.

Conclusions

When electrical stimulus pulses are grouped into periodic bursts, perceived frequency is best explained by the inter-burst interval, i.e. the duration of the silent gap between bursts, as opposed to the burst rate (periodicity) or mean pulse rate. This finding, consistent with what we have shown with mechanical pulsatile stimulation, suggests that this frequency encoding arises from central processing rather than factors such as skin mechanics. The method of burst stimulation described here may represent a way of providing sensory feedback information such as frequency for tactile displays and sensory prostheses.

Supporting information

S1 Fig. Example psychometric function.

An example of a logistic regression fit applied to psychometric data from a single subject in one test condition (Experiment 1, stimulus d). PSE (point of subjective equality) is determined as the frequency that gives logit = 0 from the fitted regression line. In this case, PSE was estimated to be 38.1 Hz.

https://doi.org/10.1371/journal.pone.0237440.s001

(EPS)

S1 Dataset. Individual subject PSE values for each experiment.

https://doi.org/10.1371/journal.pone.0237440.s002

(XLSX)

References

  1. 1. Saal HP, Wang X, Bensmaia SJ. Importance of spike timing in touch: an analogy with hearing? Curr Opin Neurobiol. 2016;40:142–9. Epub 2016/08/10. pmid:27504741; PubMed Central PMCID: PMC5315566.
  2. 2. Mackevicius EL, Best MD, Saal HP, Bensmaia SJ. Millisecond precision spike timing shapes tactile perception. J Neurosci. 2012;32(44):15309–17. pmid:23115169; PubMed Central PMCID: PMC3752122.
  3. 3. Johansson RS, Birznieks I. First spikes in ensembles of human tactile afferents code complex spatial fingertip events. Nat Neurosci. 2004;7(2):170–7. Epub 2004/01/20. pmid:14730306.
  4. 4. Saal HP, Vijayakumar S, Johansson RS. Information about complex fingertip parameters in individual human tactile afferent neurons. J Neurosci. 2009;29(25):8022–31. Epub 2009/06/26. pmid:19553442; PubMed Central PMCID: PMC6666033.
  5. 5. Birznieks I, Vickery RM. Spike timing matters in novel neuronal code involved in vibrotactile frequency perception. Curr Biol. 2017;27(10):1485–90.e2. Epub 2017/05/10. pmid:28479322.
  6. 6. Salinas E, Hernandez A, Zainos A, Romo R. Periodicity and firing rate as candidate neural codes for the frequency of vibrotactile stimuli. J Neurosci. 2000;20(14):5503–15. Epub 2000/07/07. pmid:10884334; PubMed Central PMCID: PMC6772326.
  7. 7. Mountcastle VB, Talbot WH, Sakata H, Hyvärinen J. Cortical neuronal mechanisms in flutter-vibration studied in unanesthetized monkeys. Neuronal periodicity and frequency discrimination. J Neurophysiol. 1969;32(3):452–84. Epub 1969/05/01. pmid:4977839.
  8. 8. Luna R, Hernandez A, Brody CD, Romo R. Neural codes for perceptual discrimination in primary somatosensory cortex. Nat Neurosci. 2005;8(9):1210–9. Epub 2005/08/02. pmid:16056223.
  9. 9. Romo R, Salinas E. Flutter discrimination: neural codes, perception, memory and decision making. Nat Rev Neurosci. 2003;4(3):203–18. Epub 2003/03/04. pmid:12612633.
  10. 10. McLeod KJ, Lovely DF, Scott RN. A biphasic pulse burst generator for afferent nerve stimulation. Med Biol Eng Comput. 1987;25(1):77–80. Epub 1987/01/01. pmid:3695608.
  11. 11. Pena AE, Rincon-Gonzalez L, Abbas JJ, Jung R. Effects of vibrotactile feedback and grasp interface compliance on perception and control of a sensorized myoelectric hand. PLoS One. 2019;14(1):e0210956. Epub 2019/01/17. pmid:30650161; PubMed Central PMCID: PMC6334959.
  12. 12. Szeto AY, Saunders FA. Electrocutaneous stimulation for sensory communication in rehabilitation engineering. IEEE Trans Biomed Eng. 1982;29(4):300–8. Epub 1982/04/01. pmid:7068167.
  13. 13. Grill WM. Temporal pattern of electrical stimulation is a new dimension of therapeutic innovation. Curr Opin Biomed Eng. 2018;8:1–6. Epub 2019/03/25. pmid:30906909; PubMed Central PMCID: PMC6426311.
  14. 14. Brocker DT, Swan BD, Turner DA, Gross RE, Tatter SB, Koop MM, et al. Improved efficacy of temporally non-regular deep brain stimulation in Parkinson's disease. Exp Neurol. 2013;239:60–7. Epub 2012/10/02. pmid:23022917; PubMed Central PMCID: PMC3547657.
  15. 15. Baker KB, Zhang J, Vitek JL. Pallidal stimulation: effect of pattern and rate on bradykinesia in the non-human primate model of Parkinson's disease. Exp Neurol. 2011;231(2):309–13. Epub 2011/07/20. pmid:21767534; PubMed Central PMCID: PMC3536492.
  16. 16. Birznieks I, McIntyre S, Nilsson HM, Nagi SS, Macefield VG, Mahns DA, et al. Tactile sensory channels over-ruled by frequency decoding system that utilizes spike pattern regardless of receptor type. eLife. 2019;8. Epub 2019/08/07. pmid:31383258; PubMed Central PMCID: PMC6684274.
  17. 17. Gardner EP, Palmer CI. Simulation of motion on the skin. I. Receptive fields and temporal frequency coding by cutaneous mechanoreceptors of OPTACON pulses delivered to the hand. J Neurophysiol. 1989;62(6):1410–36. Epub 1989/12/01. pmid:2600632.
  18. 18. Kaczmarek KA, Tyler ME, Brisben AJ, Johnson KO. The afferent neural response to electrotactile stimuli: preliminary results. IEEE Trans Rehabil Eng. 2000;8(2):268–70. Epub 2000/07/15. pmid:10896199.
  19. 19. Rollman GB. Electrocutaneous stimulation: Psychometric functions and temporal integration. Percept Psychophys. 1969;5(5):289–93.
  20. 20. Rollman GB. Electrocutaneous stimulation. In: Geldard FA, editor. Conference on Cutaneous Communication Systems and Devices. Austin, Texas: Psychonomic Society; 1974. p. 38–51.
  21. 21. Steenbergen P, Buitenweg JR, Trojan J, van der Heide EM, van den Heuvel T, Flor H, et al. A system for inducing concurrent tactile and nociceptive sensations at the same site using electrocutaneous stimulation. Behav Res Methods. 2012;44(4):924–33. Epub 2012/07/19. pmid:22806702; PubMed Central PMCID: PMC3509324.
  22. 22. van der Heide EM, Buitenweg JR, Marani E, Rutten WL. Single pulse and pulse train modulation of cutaneous electrical stimulation: a comparison of methods. J Clin Neurophysiol. 2009;26(1):54–60. Epub 2009/01/20. pmid:19151619.
  23. 23. Johansson RS, Landström U, Lundström R. Responses of mechanoreceptive afferent units in the glabrous skin of the human hand to sinusoidal skin displacements. Brain Res. 1982;244(1):17–25. Epub 1982/07/22. pmid:6288178.
  24. 24. Talbot WH, Darian-Smith I, Kornhuber HH, Mountcastle VB. The sense of flutter-vibration: comparison of the human capacity with response patterns of mechanoreceptive afferents from the monkey hand. J Neurophysiol. 1968;31(2):301–34. Epub 1968/03/01. pmid:4972033.
  25. 25. Manfredi LR, Baker AT, Elias DO, Dammann JF 3rd, Zielinski MC, Polashock VS, et al. The effect of surface wave propagation on neural responses to vibration in primate glabrous skin. PLoS One. 2012;7(2):e31203. Epub 2012/02/22. pmid:22348055; PubMed Central PMCID: PMC3278420.
  26. 26. Dandu B, Shao Y, Stanley A, Visell Y. Spatiotemporal Haptic Effects from a Single Actuator via Spectral Control of Cutaneous Wave Propagation. 2019 IEEE World Haptics Conference (WHC); 2019 9–12 July 2019; Tokyo, Japan. https://doi.org/10.1109/WHC.2019.8816149
  27. 27. Johnson KO. Reconstruction of population response to a vibratory stimulus in quickly adapting mechanoreceptive afferent fiber population innervating glabrous skin of the monkey. J Neurophysiol. 1974;37(1):48–72. Epub 1974/01/01. pmid:4204567.
  28. 28. Bensmaia SJ. Tactile intensity and population codes. Behav Brain Res. 2008;190(2):165–73. Epub 2008/04/19. pmid:18420286; PubMed Central PMCID: PMC2659538.
  29. 29. Yildiz KA, Shin AY, Kaufman KR. Interfaces with the peripheral nervous system for the control of a neuroprosthetic limb: a review. J Neuroeng Rehabil. 2020;17(1):43. Epub 2020/03/11. pmid:32151268; PubMed Central PMCID: PMC7063740.
  30. 30. Kajimoto H. Electro-tactile Display: Principle and Hardware. In: Kajimoto H, Saga S, Konyo M, editors. Pervasive Haptics: Science, Design, and Application. Tokyo: Springer Japan; 2016. p. 79–96. https://doi.org/10.1007/978-4-431-55772-2_5.
  31. 31. von Békésy G. Synchronism of neural discharges and their demultiplication in pitch perception on the skin and in hearing. J Acoust Soc Am. 1959;31(3):338–49.
  32. 32. Mallik A, Weir AI. Nerve conduction studies: essentials and pitfalls in practice. J Neurol Neurosurg Psychiatry. 2005;76(Suppl 2):ii23–ii31. Epub 2005/06/18. pmid:15961865; PubMed Central PMCID: PMC1765692.
  33. 33. Vickery RM, Morley JW, Rowe MJ. The role of single touch domes in tactile perception. Exp Brain Res. 1993;93(2):332–4. Epub 1993/01/01. pmid:8491272.
  34. 34. Mahns DA, Perkins NM, Sahai V, Robinson L, Rowe MJ. Vibrotactile frequency discrimination in human hairy skin. J Neurophysiol. 2006;95(3):1442–50. Epub 2005/12/02. pmid:16319219.
  35. 35. Rothenberg M, Verrillo RT, Zahorian SA, Brachman ML, Bolanowski SJ Jr., Vibrotactile frequency for encoding a speech parameter. J Acoust Soc Am. 1977;62(4):1003–12. Epub 1977/10/01. pmid:908786.
  36. 36. Goff GD. Differential discrimination of frequency of cutaneous mechanical vibration. J Exp Psychol. 1967;74(2):294–9. Epub 1967/06/01. pmid:6048472.
  37. 37. Graczyk EL, Schiefer MA, Saal HP, Delhaye BP, Bensmaia SJ, Tyler DJ. The neural basis of perceived intensity in natural and artificial touch. Sci Transl Med. 2016;8(362):362ra142. Epub 2016/11/01. pmid:27797958; PubMed Central PMCID: PMC5713478.
  38. 38. Li M, Zhang D, Chen Y, Chai X, He L, Chen Y, et al. Discrimination and recognition of phantom finger sensation through transcutaneous electrical nerve stimulation. Front Neurosci. 2018;12:283. Epub 2018/05/16. pmid:29760647; PubMed Central PMCID: PMC5937010.
  39. 39. Vallbo ÅB, Olsson KÅ, Westberg KG, Clark FJ. Microstimulation of single tactile afferents from the human hand. Sensory attributes related to unit type and properties of receptive fields. Brain. 1984;107(3):727–49. Epub 1984/09/01. pmid:6478176.
  40. 40. Romo R, Hernandez A, Zainos A, Brody CD, Lemus L. Sensing without touching: psychophysical performance based on cortical microstimulation. Neuron. 2000;26(1):273–8. Epub 2000/05/08. pmid:10798410.
  41. 41. Merzenich MM, Harrington TH. The sense of flutter-vibration evoked by stimulation of the hairy skin of primates: comparison of human sensory capacity with the responses of mechanoreceptive afferents innervating the hairy skin of monkeys. Exp Brain Res. 1969;9(3):236–60. Epub 1969/01/01. pmid:4984464.
  42. 42. McIntyre S, Birznieks I, Andersson R, Dicander G, Breen PP, Vickery RM. Temporal Integration of Tactile Inputs from Multiple Sites. In: Bello F, Kajimoto H, Visell Y, editors. Haptics: Perception, Devices, Control, and Applications. Eurohaptics 2016. Lecture Notes in Computer Science. 9774. Cham: Springer International Publishing; 2016. p. 204–13. https://doi.org/10.1007/978-3-319-42321-0_19
  43. 43. Vallbo ÅB, Olausson H, Wessberg J, Kakuda N. Receptive field characteristics of tactile units with myelinated afferents in hairy skin of human subjects. J Physiol. 1995;483(3):783–95. Epub 1995/03/15. pmid:7776258; PubMed Central PMCID: PMC1157818.
  44. 44. Bolton CF, Winkelmann RK, Dyck PJ. A quantitative study of Meissner's corpuscles in man. Neurology. 1966;16(1):1–9. Epub 1966/01/01. pmid:5948004.
  45. 45. Manfredi LR, Saal HP, Brown KJ, Zielinski MC, Dammann JF 3rd, Polashock VS, et al. Natural scenes in tactile texture. J Neurophysiol. 2014;111(9):1792–802. Epub 2014/02/14. pmid:24523522.
  46. 46. Hollins M, Bensmaia SJ, Roy EA. Vibrotaction and texture perception. Behav Brain Res. 2002;135(1–2):51–6. Epub 2002/10/03. pmid:12356433.
  47. 47. Johansson RS, Flanagan JR. Coding and use of tactile signals from the fingertips in object manipulation tasks. Nat Rev Neurosci. 2009;10(5):345–59. Epub 2009/04/09. pmid:19352402.
  48. 48. Delhaye B, Hayward V, Lefevre P, Thonnard JL. Texture-induced vibrations in the forearm during tactile exploration. Front Behav Neurosci. 2012;6:37. Epub 2012/07/12. pmid:22783177; PubMed Central PMCID: PMC3390558.
  49. 49. Weber AI, Saal HP, Lieber JD, Cheng JW, Manfredi LR, Dammann JF 3rd, et al. Spatial and temporal codes mediate the tactile perception of natural textures. Proc Natl Acad Sci U S A. 2013;110(42):17107–12. Epub 2013/10/02. pmid:24082087; PubMed Central PMCID: PMC3800989.
  50. 50. Ng KKW, Birznieks I, Tse ITH, Andersen J, Nilsson S, Vickery RM. Perceived Frequency of Aperiodic Vibrotactile Stimuli Depends on Temporal Encoding. In: Prattichizzo D, Shinoda H, Tan HZ, Ruffaldi E, Frisoli A, editors. Haptics: Science, Technology, and Applications. EuroHaptics 2018. Lecture Notes in Computer Science. 10893. Cham: Springer; 2018. p. 199–208. https://doi.org/10.1007/978-3-319-93445-7_18
  51. 51. Krahe R, Gabbiani F. Burst firing in sensory systems. Nat Rev Neurosci. 2004;5(1):13–23. Epub 2003/12/09. pmid:14661065.
  52. 52. Lisman JE. Bursts as a unit of neural information: making unreliable synapses reliable. Trends Neurosci. 1997;20(1):38–43. Epub 1997/01/01. pmid:9004418.
  53. 53. Kepecs A, Lisman J. Information encoding and computation with spikes and bursts. Network. 2003;14(1):103–18. Epub 2003/03/05. pmid:12613553.
  54. 54. Izhikevich EM, Desai NS, Walcott EC, Hoppensteadt FC. Bursts as a unit of neural information: selective communication via resonance. Trends Neurosci. 2003;26(3):161–7. Epub 2003/02/20. pmid:12591219.
  55. 55. Kaczmarek KA, Haase SJ. Pattern identification and perceived stimulus quality as a function of stimulation waveform on a fingertip-scanned electrotactile display. IEEE Trans Neural Syst Rehabil Eng. 2003;11(1):9–16. Epub 2003/06/12. pmid:12797720.
  56. 56. Perez CA, Holzmann CA, Jaeschke HE. Two-point vibrotactile discrimination related to parameters of pulse burst stimulus. Med Biol Eng Comput. 2000;38(1):74–9. Epub 2000/06/01. pmid:10829394.
  57. 57. Kaczmarek KA, Tyler ME, Okpara UO, Haase SJ. Interaction of Perceived Frequency and Intensity in Fingertip Electrotactile Stimulation: Dissimilarity Ratings and Multidimensional Scaling. IEEE Trans Neural Syst Rehabil Eng. 2017;25(11):2067–74. Epub 2017/05/16. pmid:28504942; PubMed Central PMCID: PMC5681405.
  58. 58. Vickery RM, Ng KKW, Potas JR, Shivdasani MN, McIntyre S, Nagi SS, et al. Tapping into the language of touch: using non-invasive stimulation to specify tactile afferent firing patterns. Front Neurosci. 2020;14:500. Epub 2020/06/09. pmid:32508581; PubMed Central PMCID: PMC7248323.
  59. 59. Rijnbeek EH, Eleveld N, Olthuis W. Update on peripheral nerve electrodes for closed-loop neuroprosthetics. Front Neurosci. 2018;12:350. Epub 2018/06/19. pmid:29910705; PubMed Central PMCID: PMC5992394.