Front. Pain Res.Frontiers in Pain ResearchFront. Pain Res.2673-561XFrontiers Media S.A.10.3389/fpain.2022.1112614Pain ResearchBrief Research ReportInvestigations into an overlooked early component of painful nociceptive withdrawal reflex responses in humansThorellOumie12YdreforsJohannes2SvantessonMats2GerdleBjörn3OlaussonHåkan2MahnsDavid A.1NagiSaad S.12*School of Medicine, Western Sydney University, Sydney, NSW, AustraliaDepartmentof Biomedical and Clinical Sciences, Linköping University, Linköping, SwedenPain and Rehabilitation Centre, and Department of Health, Medicine and Caring Sciences, Linköping University, Linköping, Sweden
Edited by: Craig T. Hartrick, Oakland University, United States
Reviewed by: Giancarlo carli, University of Siena, Italy Antti Pertovaara, University of Helsinki, Finland
Correspondence: Saad S. Nagi saad.nagi@liu.se
Specialty Section: This article was submitted to Pain Mechanisms, a section of the journal Frontiers in Pain Research
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Introduction
The role of pain as a warning system necessitates a rapid transmission of information from the periphery for the execution of appropriate motor responses. The nociceptive withdrawal reflex (NWR) is a physiological response to protect the limb from a painful stimulus and is often considered an objective measure of spinal nociceptive excitability. The NWR is commonly defined by its latency in the presumed Aδ-fiber range consistent with the canonical view that “fast pain” is signaled by Aδ nociceptors. We recently demonstrated that human skin is equipped with ultrafast (Aβ range) nociceptors. Here, we investigated the short-latency component of the reflex and explored the relationship between reflex latency and pain perception.
Methods
We revisited our earlier work on NWR measurements in which, following convention, only reflex responses in the presumed Aδ range were considered. In our current analysis, we expanded the time window to search for shorter latency responses and compared those with pain ratings.
Results
In both cohorts, we found an abundance of recordings with short-latency reflex responses. In nearly 90% of successful recordings, only single reflex responses (not dual) were seen which allowed us to compare pain ratings based on reflex latencies. We found that shorter latency reflexes were just as painful as those in the conventional latency range.
Conclusion
We found a preponderance of short-latency painful reflex responses. Based on this finding, we suggest that short-latency responses must be considered in future studies. Whether these are signaled by the ultrafast nociceptors remains to be determined.
aβ fiberfibromyalgianociceptive withdrawal reflexnociceptorheightened painRIIRIIIEMGThis work was supported by the Swedish Research Council (SSN and BG), Knut and Alice Wallenberg Foundation (HO), ALF Grants, Region Östergötland (SSN), Svenska Läkaresällskapet (SSN), and Western Sydney University (DAM).Introduction
The nociceptive withdrawal reflex (NWR) is a physiological response of the limb away from a painful stimulus. It has been investigated both as a tool to probe spinal nociceptive excitability and because of its close association with subjective pain thresholds (1–7). The NWR response is measured by electromyography and considered to have two latency components: the first component, referred to as RII, is mediated by Aβ or Group II fibers; the second component, referred to as RIII, is mediated by Aδ or Group III fibers (8). The conventional view is that the first or short-latency response is exclusively tactile, and the second or long-latency response is nociceptive (but also see Willer et al. 1978 (9)). While the RII-RIII latency cut-off varies across studies, the exclusion of short-latency responses is common practice (1, 10–12).
Using microneurography, we recently showed that humans, akin to other mammals, are equipped with ultrafast (Aβ range) nociceptors in the skin (13). Considering this finding, we revisited our earlier work on NWR measurements (14): in that study, following convention, NWR responses were only selected if they occurred ≥90 ms, consistent with the presumed Aδ-fiber range. Here we expanded the time window to search for shorter latency responses with the hypothesis that those are nociceptive, corresponding to painful sensations.
We found an abundance of short-latency reflex responses, and these were just as painful as those in the conventional latency range, suggesting that by discarding shorter latencies, we may be overlooking valuable quantitative measures of pain processing.
Material and methodsParticipants
NWR responses and pain ratings were successfully extracted for 20 fibromyalgia patients (FM: 24–56 years, all female) and 10 healthy controls (HC: 24–53 years, all female). For details on patient eligibility criteria, refer to Ydrefors et al. (14). Raw data were unavailable for 10 HC and therefore another 10 HC were recruited (19–39 years, all female). The new participants were not age-matched: these data were collected during the pandemic, and it was considered an unnecessary risk to recruit older participants. Additional data collection was approved by the Swedish Ethical Review Authority (Dnr: 2020–04207), and the study procedures complied with the revised Declaration of Helsinki. All participants gave their written informed consent.
Testing procedure and NWR determination
For full details on the testing procedure, refer to Ydrefors et al. (14). Briefly, electrical stimuli were delivered to the surface of the foot sole using a constant current stimulator generating a train of 5 square wave pulses (1 ms, 200 Hz), and electromyographic (EMG) responses were recorded from the ipsilateral tibialis anterior muscle. In Ydrefors et al. (14), reflex responses with Z-scores ≥12 were detected using an automated approach. The maximum amplitude (peak amplitude) in the time window of 90 to 150 ms after stimulus onset and the mean amplitude in the −65 to −5 ms pre-stimulus onset (baseline) were determined on a trial-to-trial basis. To determine the Z-score, the difference between peak amplitude and mean baseline amplitude was divided by the baseline standard deviation.
All participants rated the intensity of the sensation, immediately after receiving the electrical stimulus, on a descriptive numeric scale from 0 to 10. Zero corresponded to “no feeling”, 1 to a “slight feeling”, 2 to a “distinct feeling”, 3 to “unpleasantness”, 4 to “just noticeable pain”, 5 to “slight pain”, 6 to “distinct pain”, 7 to “moderately intense pain”, 8 to “intense pain”, 9 to “very intense pain” and 10 to the “worst imaginable pain” (14).
New data analysis
The data were pseudonymized and information on latency, Z-scores, NWR thresholds, and age were stored in a relational database. NWR responses were converted from text files (.txt) into graphs (.png), using a script made in Python Distribution (v3.7.4, Python Software Foundation, Beaverton, United States) and latencies were visually inspected by the author (OT) in LabChart (v8.1.16 ADInstruments, Dunedin, New Zealand) and in MATLAB (r2021b, MathWorks Inc, Natick, Massachusetts). Z-scores were calculated for the early time window of 50 to 89 ms after stimulus onset, using the same MATLAB algorithm that was used for the 90 to 150 ms time window. Careful visual inspection of the data allowed us to extract reflex responses with Z-scores ≥6.
Statistical analysis
Statistical analysis was done in GraphPad Prism (v. 9.1.2, GraphPad Software, San Diego, United States). QQ plots, means, standard deviations, and skewness were assessed to determine the normal distribution of the data. Where possible, parametric tests were used. If assumptions for parametric tests were violated, non-parametric tests were conducted.
To compare independent differences between HC and FM, unpaired two-tailed t-test was used, while small sample size data were analysed using two-tailed Mann-Whitney U test. Two-way ANOVA was used to compare multiple independent groups with Tukey's test as a multiple comparison (post-hoc) test. Only main effects were analyzed due to uneven sample sizes. A statistical value of p < 0.05 was considered statistically significant.
Effect sizes were calculated as Hedges’ g for the unpaired t-tests, due to different sample sizes, and partial eta square (ƞ2p) for the two-way ANOVA. Common language effect sizes (CLES) are shown for statistically significant ANOVA results. When using non-parametric tests, CLES is shown to compare effects to the parametric results. Effect sizes were calculated in statistical calculators (15, 16). Numbers are presented as mean and standard deviation for parametric tests and median and interquartile range for non-parametric tests.
Results
Three hundred and eighty-two painful reflex responses with Z-scores ≥6 were successfully extracted from 340 EMG recordings: 166 NWR responses in 20 HC and 216 NWR responses in 20 FM (Figures 1A,B). The Z-scores ranged from 6.1 to 726.4 with rectified amplitudes of 6 to 698 mV. These reflex responses corresponded to ratings of 4 (“just noticeable pain”) and higher. 63 reflex responses (14.2% of total (382 + 63)) corresponded to ratings below 4 (i.e., nonpainful).
(A). Examples of reflex recordings with RII, RIII, and dual responses superimposed. The five peaks at the beginning of the graph represent the electrical stimulus (5 square pulses). (B). Latency spread of painful NWR responses. Healthy controls (HC) had a bimodal distribution while fibromyalgia patients (FM) had a more even distribution throughout the time analysis window. The y-axis shows the number of NWR responses, and the x-axis shows reflex latencies. (C). Latencies of all NWR responses from HC and FM. Latencies did not differ between the two groups (HC: 95.7 ± 23.5 ms, FM: 95.7 ± 24.0 ms, t (380) = 0.043, p = 0.965, 95% CI [-4.722, 4.936], Hedges’ g = 0.004, CLES = 50.1%). (D). Stimulus intensities of all NWR responses in HC and FM. Stimulus intensities were not different between HC and FM (HC: 15.0 ± 5.0 mA, FM: 14.7 ± 6.4 mA, t (380) = 0.568, p = 0.570, 95% CI [-1.524,.840], Hedges’ g = 0.059, CLES = 51.7%) but were significantly higher for eliciting a painful reflex than just pain regardless of subject type (at least p < 0.05, Wilcoxon test). (E). Duration of silent EMG period intervening a dual reflex response occasionally seen. No statistical difference was found in the duration of the silent period between HC and FM (HC: 50.0 (16.0) ms, FM: 47.0 (15.6) ms, U = 203, p = 0.667, CLES = 52.1%).
Reflex latencies and stimulus intensities did not differ between HC and FM groups (Figures 1C,D). In 42 out of 340 EMG recordings (21 each in HC and FM groups, 12.4% of the total), two reflex responses were seen (84 reflex responses). These dual responses were separated by a silent EMG period (SP) with a duration of 23 to 67 ms (mean 49.1 ms). The SP duration was not different between HC and FM groups (Figure 1E).
In terms of RII-III prevalence, 192 (50.3%) reflex responses were identified in the 90- to 150-ms latency range corresponding to RIII, and 190 reflexes (49.7%) were identified in the 50- to 90 ms latency range corresponding to RII. The RII data are new: the pre-set 90 ms latency cut-off implemented in Ydrefors et al. (14) resulted in an automatic discounting of shorter latency responses. To compare pain with NWR responses, only those reflex responses that were painful (at least a 4 rating on a 0–10 scale) were included in the main analysis.
Comparison of RII and RIII responses between FM and HC groups
Eighteen out of 20 HC (90%) and 13 out of 20 FM (65%) had an RII. No differences were found in stimulus intensities between RII and RIII (p = .717) (Figure 2A). To determine the relationship between subjective pain rating and reflex latency, only single reflex responses were considered (298 reflex responses) (Figure 2B). FM had higher pain ratings than HC for both RII and RIII responses. Within each group (FM/HC), when pain ratings were compared between RII and RIII, they were not different, suggesting that RII was just as painful as RIII (Figure 2C). A small proportion of the NWR responses were non-painful: 58 in the HC group and 5 in the FM group (Figure 2D).
(A). Stimulus intensities required to evoke RII and RII in HC and FM. A pre-set cut-off of 90 ms was implemented to separate RII and RIII responses. The stimulus intensities required to evoke RII and RIII responses were not different (F (1, 283) = 0.131, p = 0.717, CI [-0.292, 0.424], η2p < 0.000, CLES = 50.0%). There was a main effect of subject type (HC or FM) (F (1, 283) = 79.9, p < 0.000, η2p = 0.022, CLES = 77.3%) but post hoc test indicated no differences in stimulus intensities for subject or reflex type. (B). Proportion of single and dual NWR EMG recordings. The dual recordings (84 reflex responses) were excluded from subsequent perception analysis. (C). Pain ratings corresponding to RII and RIII response in HC and FM. Simple main effects indicated that the reflex type had no effect on pain ratings (F (1, 295) = 0.011, p = 0.916, CI [-0.333, 0.370], η2p < 0.001, CLES = 50.0%). Subject type (HC or FM) did have a large effect on pain ratings (F (1, 295) = 82.6, p < 0.001, CI [-2.001, −1.288], η2p = 0.218, CLES = 77.2%). (D). Latency spread of non-painful NWR responses. Only a few NWR responses (14.6%) were reported as non-painful (pain rating <4), almost entirely by HC (HC: 82.0 (15.5) ms, n = 58. FM: 96.0 (25.0) ms, n = 5). These responses were not included in our analysis. The y-axis shows the number of NWR responses, and the x-axis shows reflex latencies.
As an alternative to a pre-set latency cut-off for RII/III (90 ms), we used the data from dual reflex responses to distinguish between RII and RIII latencies. No dual responses occurred before 93 ms or after 99 ms, therefore we took an in-between value of 96 ms to separate RII and RIII responses (Figure 3A). Predictably, this increased the proportion of RII responses: in the HC group, the increase was 13.3% (22 additional responses) and in the FM group, the increase was 13.9% (an additional 30 responses) (Figure 3B,C). However, the overall results did not change (Figure 3D,E).
(A). Latencies of all dual reflex responses. No dual responses were observed before 93 or after 99 ms (dotted lines) so an in-between value of 96 ms was taken to separate RII and RIII responses (solid black line). (B-C). Proportion of RII and RIII with a divider set at 96 ms. The pie chart on the left shows NWR responses with a pre-set cut-off at 90 ms, and the pie chart on the right shows NWR responses with a divider set at 96 ms, a cut-off derived from dual responses. Predictably, in both HC and FM groups, the prevalence of RII increased and the prevalence of RIII decreased when a 96 ms cutoff was implemented. (D-E). Stimulus intensities and pain ratings for RII and RIII in HC and FM using 96 ms as the divider. Current intensities had an effect on reflex type (F (1, 379) = 5.262, p = 0.022, CI [-2.630, −0.202], η2p = 0.013, CLES = 56.5%), but not on subject type (F (1,379) = 0.587, p = 0.444, CI [-0.720, 1.640], η2p < 0.001, CLES = 50.0%). Post hoc test indicated no difference in stimulus intensities for subject or reflex type. Subject type had an effect on pain ratings (F (1, 295) = 84.8, p < 0.001, CI [-2.022, −1.310], η2p = .223, CLES = 77.6%) but not reflex type (F (1, 295) = 1.125, p = 0.290, CI [-0.173, 0.577], η2p = .003, CLES = 53.1%).
Z-scores
In Ydrefors et al. (14), an automated method was used to detect reflex responses with Z-scores of ≥12. Here, visual inspection of the data allowed us to include responses with Z-scores ≥6. For comparison, we performed the analysis implementing the original Z-score (≥12) condition. A total of 234 NWR responses had a Z-score of 12 or higher (Supplementary Figure S1A). The proportion of dual responses increased from 12.4% to 17.6% (Supplementary Figure S1B). Reflex latencies were shorter in the FM group compared to the HC group (Supplementary Figure S1C). Stimulus intensities did not differ between HC and FM (Supplementary Figure S1D). FM had higher pain ratings than HC, but pain ratings and stimulus intensities were not different between RII and RIII for either group (Supplementary Figure S1E-F).
Discussion
In this study, we reanalyzed the reflex data from Ydrefors et al. (14): expanding the time window revealed an abundance of RII NWR responses at stimulus intensities deemed painful. Remarkably, most recordings contained a single reflex response, contrary to the notion that the reflex usually consists of a double burst of EMG activity. We excluded the dual responses from the perception analysis and compared reflex latency (all single responses) with corresponding pain ratings. We found that RII responses were just as painful as RIII responses. The canonical view is that the short-latency component of the NWR response is purely tactile (i.e., nonpainful) and signaled by Aβ low-threshold mechanoreceptors. However, our data show a preponderance of painful reflexes with short latencies.
In the literature, the proposed RII-RIII latency cut-off can be anywhere between 60 and 115 ms post-stimulus onset (for references, see 17). Here, we compared the data based on a pre-set latency cut-off of 90 ms with a post hoc approach using dual responses to set the RII-RIII cut-off; this did not change the overall results. In the early reflex work showing double-burst EMG activity, the first EMG response had a lower electrical threshold (18, 19) and was less painful or non-painful (8, 1). In our study, nonpainful reflex responses were rarely observed. Further, our focus on short-latency limb reflexes avoided the confound of startle responses which are of longer latency (20, 11). Using dual nerve stimulation, Willer and colleagues were able to evoke an NWR response while the small-fiber inputs were blocked by anesthesia (9). Indeed, in our single-unit microneurography study, we had confirmed that Aβ nociceptive fibers not only respond to and encode nociceptive stimuli, but also evoke a painful percept when selectively activated during intraneural electrical stimulation (13).
In the FM group, pain ratings were higher and nearly all reflex responses were painful. We observed no latency differences between HC and FM groups except when only responses with Z-scores ≥12 were considered; in that case, FM had shorter latencies than HC. Higher pain ratings and shorter latencies could be attributed to peripheral and/or central sensitization in the patient population (21). We did not detect differences in the SP duration between the two groups. The SPs arise due to postsynaptic inhibition in the motor neurons following a strong electrical stimulus of a muscle or cutaneous nerve. Prolonged SP duration has been previously observed in FM patients (22) and is thought to reflect spinal dysregulation.
In Ydrefors et al. (14), a Z-score of ≥12 was considered a successful muscle response. In that study, a fully automated method was used for NWR detection, therefore the Z-score had to be large enough to ensure that noise would not be interpreted as muscle response. In the current study, we visually inspected all data and were able to reliably detect responses with Z-scores ≥ 6, increasing our sample size by over a third. The conclusions were the same regardless of the Z-score threshold.
Limitations
Raw data were not available from 10 HC in the original sample, therefore new participants had to be recruited. The recruitment was done during the pandemic, and it was considered an unnecessary risk to recruit older participants, thus the new sample is not age matched. Other than that, care was taken to ensure that the experimental protocol was as similar as possible to the original study. A comparison between the two HC samples did not reveal any differences in reflex thresholds or pain ratings.
We did not calculate the conduction velocity (CV) of the NWR in the current study. One study based on afferent CVs from single painful shocks with near-nerve stimulation of the tibial nerve reported velocities of 18.5 ± 1.3 m/s with onset latencies between 100 and 200 ms (23). Another study based on a train of 5 pulses reported conduction velocities of 49.0 ± 11.3 m/s with onset latencies between 50 and 100 ms (24). Classification based on conduction velocity into Aβ and Aδ groups is not clear-cut in humans. In animal studies, the D-hair units are considered a benchmark for the Aδ velocity range (25), however, in humans no detailed account of D-hair units exists.
Conclusions
We found a great many short-latency NWR responses that were as prevalent and as painful as the conventional longer latency NWR responses. Reflex responses that were not painful rarely occurred. Only a minority of NWR recordings consisted of dual reflexes. Pain ratings were similar across all latencies, suggesting that the short-latency component is not tactile but nociceptive. Whether this fast signaling involves Aβ nociceptors remains to be determined.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by Swedish Ethical Review Authority. The patients/participants provided their written informed consent to participate in this study.
Author contributions
OT, DAM, SSN and HO contributed to the conception and design of this study. MS optimized data analysis. OT wrote the first draft of the manuscript. DAM, SSN, HO, JY, BG contributed to subsequent revisions of the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Swedish Research Council (SSN and BG), Knut and Alice Wallenberg Foundation (HO), ALF Grants, Region Östergötland (SSN), Svenska Läkaresällskapet (SSN), and Western Sydney University (DAM).
Acknowledgments
We thank Magnus Kronander for his helpful contribution to data analysis.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpain.2022.1112614/full#supplementary-material.
Referencesde WillerJC. Comparative study of perceived pain and nociceptive flexion reflex in man. Pain. (1977) 3(1):69–80. 10.1016/0304-3959(77)90036-7876668ChanCWDallaireM. Subjective pain sensation is linearly correlated with the flexion reflex in man. Brain Res. (1989) 479(1):145–50. 10.1016/0006-8993(89)91344-92924143GuieuRBlinOPougetJSerratriceG. Nociceptive threshold and physical activity. Can J Neurol Sci. (1992) 19(1):69–71. 10.1017/S03171671000425661562910SandriniGArrigoABonoGNappiG. The nociceptive flexion reflex as a tool for exploring pain control systems in headache and other pain syndromes. Cephalalgia. (1993) 13(1):21–7. 10.1046/j.1468-2982.1993.1301021.x8448783Biurrun ManresaJANeziriAYCuratoloMArendt-NielsenLAndersenOK. Test-retest reliability of the nociceptive withdrawal reflex and electrical pain thresholds after single and repeated stimulation in patients with chronic low back pain. Eur J Appl Physiol. (2011) 111(1):83–92. 10.1007/s00421-010-1634-020814801LimECSterlingMPedlerACoombesBKVicenzinoB. Evidence of spinal cord hyperexcitability as measured with nociceptive flexion reflex (NFR) threshold in chronic lateral epicondylalgia with or without a positive neurodynamic test. J Pain. (2012) 13(7):676–84. 10.1016/j.jpain.2012.04.00522687426LimECWSterlingMStoneAVicenzinoB. Central hyperexcitability as measured with nociceptive flexor reflex threshold in chronic musculoskeletal pain: a systematic review. Pain. (2011) 152(8):1811–20. 10.1016/j.pain.2011.03.03321524852HugonM. Exteroceptíve reflexes to stimulation of the sural nerve in normal man. New Development in Electromyography and Clin Neurophysiol. (1973) 3:713–29. 10.1159/000394186WillerJCBoureauFAlbe-FessardD. Role of large diameter cutaneous afferents in transmission of nociceptive messages: electrophysiological study in man. Brain Res. (1978) 152(2):358–64. 10.1016/0006-8993(78)90264-0209873DowmanR. Spinal and supraspinal correlates of nociception in man. Pain. (1991) 45(3):269–81. 10.1016/0304-3959(91)90051-X1876436DowmanR. Possible startle response contamination of the spinal nociceptive withdrawal reflex. Pain. (1992) 49(2):187–97. 10.1016/0304-3959(92)90142-X1608645RhudyJLFranceCR. Defining the nociceptive flexion reflex (NFR) threshold in human participants: a comparison of different scoring criteria. Pain. (2007) 128(3):244–53. 10.1016/j.pain.2006.09.02417070999NagiSSMarshallAGMakdaniAJarockaELiljencrantzJRidderstromMAn ultrafast system for signaling mechanical pain in human skin. Sci Adv. (2019) 5:7. 10.1126/sciadv.aaw1297YdreforsJKarlssonTWentzel OlaussonUGhafouriBJohanssonACOlaussonHAutomated nociceptive withdrawal reflex measurements reveal normal reflex thresholds and augmented pain ratings in patients with fibromyalgia. J Clin Med. (2020) 9:6. 10.3390/jcm9061992LakensD. Calculating and reporting effect sizes to facilitate cumulative science: a practical primer for t-tests and ANOVAs. Front Psychol. (2013) 4:863. 10.3389/fpsyg.2013.0086324324449LenhardWLA (2016). Computation of effect sizes.https://www.psychometrica.de/effect_size.html[Accessed July 05, 2022].AndersenOK. Studies of the organization of the human nociceptive withdrawal reflex. Focus on sensory convergence and stimulation site dependency. Acta Physiol (Oxf). (2007) 189(Suppl 654):1–35. 10.1111/j.1748-1716.2007.01706.x17439638HagbarthKE. Spinal withdrawal reflexes in the human lower limbs. J Neurol Neurosurg Psychiatr. (1960) 23:222–7. 10.1136/jnnp.23.3.222ShahaniBTYoungRR. Human flexor reflexes. J Neurol Neurosurg Psychiatr. (1971) 34(5):616–27. 10.1136/jnnp.34.5.616DavisM. (1984). “The mammalian startle response,” in neural mechanisms of startle behavior, ed.EatonR. C. (New York: Plenum Press), 287–351.BoureauFLuuMDoubrereJF. Study of experimental pain measures and nociceptive reflex in chronic pain patients and normal subjects. Pain. (1991) 44(2):131–8. 10.1016/0304-3959(91)90126-I2052379BaekSHSeokHYKooYSKimBJ. Lengthened cutaneous silent period in fibromyalgia suggesting central sensitization as a pathogenesis. PLoS ONE. (2016) 11(2):e0149248. 10.1371/journal.pone.014924826871583ErtekinCErtekinNKarciogluM. Conduction velocity along human nociceptive reflex afferent nerve fibres. J Neurol Neurosurg Psychiatr. (1975) 38(10):959–65. 10.1136/jnnp.38.10.959EllrichJTreedeRD. Convergence of nociceptive and non-nociceptive inputs onto spinal reflex pathways to the tibialis anterior muscle in humans. Acta Physiol Scand. (1998) 163:391–401. 10.1046/j.1365-201X.1998.t01-1-00392.x9789583DjouhriLLawsonSN. Abeta-fiber nociceptive primary afferent neurons: a review of incidence and properties in relation to other afferent A-fiber neurons in mammals. Brain Res Brain Res Rev. (2004) 46(2):131–45. 10.1016/j.brainresrev.2004.07.01515464202