Abstract

A novel noninvasive technique was applied to measure velocity within slow conducting axons in the distal extreme of the sciatic nerve (i.e., digital nerve) in a rat model. The technique is based on the extraction of rectified multiple unit activity (MUA) from in vivo whole nerve compound responses. This method reliably identifies compound action potentials in thinly myelinated fibers conducting at a range of 9–18 m/s ( axons), as well as in a subgroup of unmyelinated C fibers conducting at approximately 1-2 m/s. The sensitivity of the method to C-fiber conduction was confirmed by the progressive decrement of the responses in the 1-2 m/s range over a 20-day period following the topical application of capsaicin (ANOVA ). Increasing the frequency of applied repetitive stimulation over a range of 0.75 Hz to 6.0 Hz produced slowing of conduction and a significant decrease in the magnitude of the compound C-fiber response (ANOVA ). This technique offers a unique opportunity for the noninvasive, repeatable, and quantitative assessment of velocity in the subsets of and C fibers in parallel with the evaluation of fast nerve conduction.

1. Introduction

The evaluation of peripheral nerve conduction velocity (NCV) (CAP: capsaicin, HFRS: high frequency repetitive stimulation, MUA: multiple unit activity, NCV: nerve conduction velocity, NFD: nerve fiber density, and SFN: small fiber neuropathy) is a well-established and highly utilized procedure for the assessment of the pattern, symmetry and extent of neuropathy in both pre-clinical and clinical studies [14]. When properly performed, NCV provides a sensitive, objective, and specific index of the onset and progression of multiple forms of nerve damage. C fibers make up the majority of fibers in some nerve segments; however, standard NCV measures are only sensitive to activity in fast conducting, heavily myelinated axons. The in vivo assessment of conduction in small diameter axons is a formidable task due to the requirement of high-threshold stimulation, asynchronous, slow conduction and phase cancellation with distance [5]. The inability to assess conduction, in thinly myelinated Aδ axons or in unmyelinated C fibers is a severe limitation of NCV procedures. There is a growing recognition that many forms of idiopathic small fiber neuropathies (SFNs) may have been poorly diagnosed and their prevalence underestimated due, in part, to a lack of sensitive, non-invasive measures of this condition [510]. SFN can result in pain, paresthesias, abnormal thermal sensitivity, and various autonomic dysfunctions. Emerging data suggests that SFNs have unique etiologies and distinct genetic contributions and that they potentially differ in their response to therapeutic approaches [1124].

The clinical diagnosis of SFN is generally based on symptoms or the presence of specific sensory deficits revealed by questionnaires and/or physical examination (e.g., diminished perception of “sharp”). These measures are subjective and often highly variable. The addition of epidermal fiber density has added sensitivity to the assessment of SFNs; however, the interpretation of these findings is often limited by wide ranges of “normal values” and differences in fiber counts as a function of location [2537].

The situation is even more difficult in preclinical models of SFN, which generally are forced to rely on behavioral assessment as a measure of small fiber dysfunction [3845]. The behaviors in question (e.g., thermal withdrawal reaction) are often complex; they involve both CNS and PNS functions, and they can be influenced by learning and motivation. Although neuropathologic evaluation can be applied to animal models, the accurate evaluation of structural deficits in C fibers requires invasive, time-consuming, teased fiber techniques [11, 14, 46, 47]. Epidermal fiber density measures have also been applied to preclinical studies, but the sensitivity of the procedures has often been limited by the same high variability issues seen in the clinic [4852]. Single unit recording provides the most complete and specific evaluation of the function of slow conducting fibers. These studies are critical for defining cellular and molecular mechanisms of SFN, but they require complicated, invasive procedures often restricted to the proximal segments of the nerves and have limited applications in assessment of the disease progression and/or in screening of potential therapeutics [5363]. Microneurographic studies have reported that painful SFNs are often associated with an increase in spontaneous activity and changes in excitability of C fibers [6467]. The introduction of microneurography to animal modeling confirmed similar patterns of C-fiber activity in rats and humans in both normal and pathological conditions; however, the accurate tracing of C-fiber activity in animal models required exposing the sciatic nerve, precluding longitudinal assessments [68, 69].

The present study was designed to define a reliable experimental technique for the assessment of slow conducting axons. The approach differs from the well-established use of single axon recording in several important ways: (1) it assesses whole nerve compound responses, allowing quantitative evaluation of the effects of experimental treatment or disease on both fast and slow conducting axons in parallel; (2) it is noninvasive, permitting repeat measurement of the slow nerve conduction as deficits progress and possibly recover; (3) it is a relatively simple procedure that could be implemented, in most cases, with existing electrophysiologic equipment. Application of this technique could augment the specificity of microneurography with longitudinal and quantitative assessment of activity in the whole intact nerve.

2. Materials and Methods

2.1. Animals

Electrophysiological studies were performed on 15 Sprague-Dawley female rats (Charles River Lab., Wilmington, MA) that were 5-6 months old at the start of the study. All animals were given free access to standard laboratory chow and water throughout the study. Five rats were treated with the topical application of capsaicin (CAP). The magnitude of C-fiber conduction in these rats were compared to that determined in 5 age-matched control rats at baseline (prior to the treatment), 10, 20, and 30–45 days after the administration of CAP. The remaining rats ( ) were used to help optimize the recording procedures and analysis. The use of all animals and procedures was in accordance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals and approved by the Animal Institute Care and Use Committee at Albert Einstein College of Medicine.

2.2. Electrophysiological Procedures

All experiments were done under general anesthesia (isoflurane 1.5–2.5%/O2), with the rats positioned on a heating pad and monitored for respiration and rectal temperature. Evoked activity was recorded in the digital nerve using subcutaneous platinum needle electrodes (Grass Astro-Med, Inc., Warwick, RI) for both recording and stimulation. The digital nerve was stimulated at the dorsal-superficial surface of the third cuneiform/proximal end of metatarsal (Figure 1, S1). A distal set of recording electrodes was positioned on the lateral side of the third toe (Figure 1, R1) 20 mm distal from the stimulation point. Another set of recording electrodes was positioned at the sciatic notch (75–85 mm proximally to the stimulation point) to monitor nerve conduction over the proximal segment of the nerve (Figure 1, R2). After evaluation of the proximal nerve conduction, percutaneous injection of 0.1–0.3 mL of 2% lidocaine (SPARHAWK, Laboratories, Inc., Lenexa, KS) was applied between the recording sites to block nerve conduction [70] and prevent the masking of slow conducting responses by central reflexes (e.g., dorsal root reflexes). The effective blocking of conduction in fast fibers was verified by the absence of neural response at the proximal site and monitored throughout the experiment. Responses were evoked by constant current electrical stimulation (Grass stimulator-11, isolation unit SIU-7, Warwick, RI), which was systematically increased in intensity (0.1–15 mA and duration 0.01–3.0 ms) until the magnitude of the C-fiber signal (i.e., velocities <2.0 m/s) was maximal.

The effects of high frequency repetitive stimulation (HFRS) on C-fiber activity were explored using stimulation rates of 0.75, 1.5, 3.0, and 6.0 Hz. The use of this procedure has been previously described for the evaluation of fast conducting fibers [71]. HFRS was conducted in step-like order from the lowest (0.75 Hz) to the highest (6.0 Hz) frequency. Prior to each escalation in frequency, a 0.75 Hz signal was evaluated to test the return of the slow fiber neural activity to the baseline level. Neuroelectric signals were impedance-matched and differentially amplified with a gain of 20,000 and a frequency band of 20 Hz–3 kHz using a BIOPAC MP100 and AcqKnowledge software version 3.7.2 (BIOPAC Systems Inc., Goleta, CA). Data were scanned for artifacts, digitized at the rate of 20 kHz, and averaged across 50–100 stimuli for slow nerve conduction and 5–20 stimuli for the fast nerve conduction. Single sweeps of slow nerve conduction were stored in a digital format (MDI DT-1600) for further off-line analysis.

2.3. Multiple Unit Activity (MUA)

Measures of MUA were extracted from single sweep recording by first high-pass filtering of the response above 450 Hz and then full-wave rectifying of the signal prior to averaging [7274]. Rectification minimizes phase cancellation of the asynchronous neural signals of slow conduction during averaging. MUA data represent composite voltages of the averaged whole nerve responses digitized at a constant sampling rate (i.e., voltage per time bin). At the specific latencies over a fixed distance from the site of stimulation, these responses were integrated across a specific time period to reflect the signal strength in a specific range of velocities (i.e., 1 to 2 m/s). To balance the difference between the time periods of faster and slower velocities (e.g., 3 ms for 1.5–2 m/s and 20 ms for 0.5–1 m/s), the response integrated over each range of NCV was divided by the number of time bins (voltage samplings) within this range. The integrated MUA signal was expressed as a magnitude of response over the range of 0.5–2 m/s with respect to the level of spontaneous activity in the prestimulation epoch [75, 76].

2.4. Application of Capsaicin

Local skin denervation was induced by topical application of 0.1% CAP crème (Capzasin-HP, Chattem, Inc., Chattanooga, TN) to the third toe of rat hind paw. CAP was applied under continuous light isoflurane anesthesia (0.5%–1.0% in 0.3 L/min flow of O2) for 10 hours/day over a four-day period. This procedure minimized discomfort of the animal during CAP application and assured consistency of treatment across animals. After each daily application, the residual crème was cleaned from the skin. The first post-CAP electrophysiologic measures occurred 10 days after discontinuation of CAP treatment; at that time, there was no evidence of inflammation or skin redness.

2.5. Statistics

Comparisons across conditions utilized standard inferential procedures (i.e., ANOVA); alpha levels were set at 0.05. If warranted, a Tukey’s Multiple Comparison Test was used for cross-sectional evaluations. Experimental data are presented as mean SD, unless otherwise indicated.

3. Results

3.1. Distal Digital Nerve Electrophysiology

The patterns of digital nerve activity varied as a function of stimulus intensity and distal-to-proximal recording site. The lowest intensity evaluated (5 mA, 0.01 ms) uniquely activated low threshold axons and resulted in a composite response dominated by velocities over 55 m/s in the proximal nerve segments and between 26 and 42 m/s at the distal segments (Figure 2(a)). Increasing the stimulation intensity (10–15 mA-0.01–0.05 ms) added the activation of slower conduction fibers with velocities less than 20 m/s (semilate responses, Figure 2(a)). In the distal segment of digital nerve, maximal NCV was  m/s with an amplitude of μV; the NCV of semilate responses was evaluated by peak latency and varied between 9.0 and 18.2 m/s (first peak mean NCV—  m/s). Simultaneously recorded neural response at the sciatic notch (proximal segment) demonstrated maximal NCV of  m/s reflecting conduction in Aαβ fibers. These measurements are in agreement with previously published evaluations of maximal NCV and single unit recordings [62, 7779] as well as histological evaluations of the morphometric fiber composition of the digital nerve as a distal branch of superficial peroneal nerve [80].

The activation of C fibers with associated conduction velocities less than 2.0 m/s required the use of high intensity stimulation (15 mA, 0.5−2 ms, Figure 2(b)). These stimulation currents and durations produced artifactual contamination of the earliest responses, and therefore the analysis of C-fiber activity was limited to the later MUA signals. The maximum magnitude of the slow NCV was observed in the range of 1-2 m/s, with the peak of activity in a range of 1.5–1.7 m/s (Figure 3). As expected, C-fiber responses extracted from whole nerve recordings were small and variable across animals (Figure 3(a)). Data recorded at a fixed distance from the stimulating cathode were averaged across subjects to improve the signal to noise ratio for the low amplitude, slow conducting responses. MUA was integrated across fixed time windows corresponding to velocities in the ranges of 2.0–1.6 m/s, 1.5–1.1 m/s, and 1.0–0.5 m/s and compared to the strength of the signal prior to stimulation. As a result, each compound C-fiber response was evaluated for a total magnitude of activity and its distribution in the range of velocities from 0.5 to 2 m/s (Figures 3(b) and 3(c)).

3.2. Effects of Topical CAP

The topical application of CAP resulted in a progressive decrease in the magnitude of the C-fiber compound response in the digital nerve over first 20 days after the end of treatment. Figure 4 illustrates the effects of CAP-induced denervation for responses with velocities in the C-fiber range. The effects of CAP are clearly seen in MUA averaged across 5 rats at four time points (Figure 4(a)). Evaluation of total activity of C fibers in the range of slow NCV (0.5–2 m/s) across time points confirmed a statistically significant decrease in activity (Figure 4(b), ANOVA ). In addition to a decrease in total magnitude of the response, C fibers also demonstrated redistribution of MUA with a shift towards a slower range of velocities at the early stage of denervation. Untreated age-matched controls demonstrated no change in the same measures of C-fiber activity. The recovery of slow NCV was observed in a period of 30–45 days after drug application.

3.3. Effects of HFRS

Conduction in C-fiber activity was clearly altered by increasing frequency rates. The initial decline in the magnitude of the C-fiber activity was evident with an HFRS of 1.5 Hz; at a rate of 6.0 Hz, the C-fiber response had been reduced by approximately 50% (ANOVA ) compared to responses to low rates of stimulation (Figure 5(a)). In parallel with changes in magnitude, the distribution of velocities in the C-fiber range was also shifted to slower response as a function of increase in stimulation frequency (Figure 5(c)). Surprisingly, C-fiber response recorded at the frequency of 0.75 Hz right after HFRS demonstrated minimal or no decline in total MUA (Figure 5(b), 0 min), while a recording after a 5 min rest period showed a significant reduction in magnitude (Figure 5(b), 5 min, ANOVA ). This paradox appears more prominent at the highest rates of stimulation (e.g., 3 and 6 Hz) and is most likely related to the slowing of Aδ fibers into the range of C-fiber activity under the pressure of the elevated rates of stimulation. Immediately after termination of HFRS, the responses of some Aδ fibers may still be delayed and cause an increase in compound response within the range of C-fiber conduction, while after 5 min rest, they recover and shift back to the shorter latencies, leaving the compound C-fiber response at the level of the magnitude comparable to the one during HFRS. Figure 5(d) provides an example of changes in representative traces of the distal Aβ and Aδ responses under conditions of HFRS. The slowest subset of Aδ fibers (NCV range of 3-4 m/s) diminished the amplitude of responses during HFRS of 3–6 Hz and, therefore, may slow down into the range of conduction less than 2 m/s. As expected, the population of faster Aδ fibers (e.g., NCV range of 7–12 m/s) demonstrated better resistance to the elevation of frequency and only slightly slowed conduction at the frequency of 6 Hz, while response in Aβ fibers (peak NCV—17 m/s) did not change at any frequencies applied.

4. Discussion

The application of MUA techniques to in vivo electrophysiological measures provides an opportunity to measure the magnitude and range of NCV in C fibers in parallel with the assessment of conduction in the faster Aβ and Aδ fibers. This approach can identify velocities less than 2 m/s and can be applied repeatedly to the vulnerable distal segments of specific nerves. In the present study, up to 85% of C-fiber activity was recorded in the range of NCV between 1 and 2 m/s. The results are comparable to those observed in microneurographic recordings in humans [66, 8183]. However, single unit studies, especially with the use of teased fiber preparation, can identify C-fiber velocities less than 1 m/s [54, 56, 60, 62, 63, 77, 84]. In the noninvasive, whole nerve approach used in the present study, the temporal dispersion and phase cancellation of activity in the slowest subset of the C fibers may simply be too great for the summed response to significantly exceed prestimulus baseline levels. In addition, the use of stimulation of 0.75 Hz and above in the present study may be too rapid for activation of the slowest fibers, which in the single fiber studies are driven at the rate of 0.25–0.5 Hz [55, 63, 68, 82]. The increased sensitivity to activity in the subset of neurons conducting slower than 1 m/s in the single fiber techniques is gained at the expense of the sampling of a limited number of axons and/or the need for mechanical manipulation of the tissue (e.g., teased fiber preparation) which is generally performed at the more proximal nerve segments.

CAP selectively affects distal endings of Aδ and C fibers [85]. Prolonged topical application of CAP leads to local skin desensitization and depletion of peripheral endings of small diameter fibers evident by immunohistochemical markers [8691]. We employed CAP to confirm that MUA in the range of slow NCV originates from the activity in C fibers and to explore the ability of new electrophysiological evaluation to register small fiber degeneration and regeneration. The compound C-fiber response was significantly reduced during the first 20 days after CAP application and recovered in 30–45 days. These functional measures were consistent with the results observed in clinic for the same population of axons evaluated by epidermal and subepidermal nerve fiber density (NFD) and quantitative sensory testing [90, 91]. In humans, a reduction of NFD in the skin was observed in 2–7 days after topical CAP application, while initial reinnervation of CAP-treated skin was marked between 20 and 40 days after application [91]. The longer time period needed for electrophysiological registration of denervations is most likely related to the greater extent of axonal degeneration which is required for the elimination of nerve compound response compared to the reduction of density of nerve endings in the skin and subcutaneous tissue. At the early stage of degeneration (10 days), C fibers demonstrated a redistribution of MUA toward a slower range of NCV, in addition to a decrease in total magnitude of the response. Similar signs of dysfunction were observed in distal diabetic autonomic neuropathy [76].

Assessment of activity-dependent changes is an important experimental tool for the characterization of C-fiber function. While fast conducting fibers can easily follow frequencies up to 100 Hz, slow conducting fibers are known for significant limitations in the ability to respond to rapid stimulation [53, 54, 57, 61, 63, 64, 66, 8184, 9296]. Using repetitive stimulation, conduction of C fibers is altered at rates of 1 Hz and, in some circumstances, blocked at frequencies greater than 2 Hz [81, 92, 97]. In the present study, HFRS altered C-fiber compound response in two ways: higher stimulation rates reduced the overall response magnitude of the C-fiber response and shifted the distribution of NCV to slower values. The initial drop in magnitude of the compound response was observed at 1.5 Hz, while stimulation at 3 Hz increased the slowing of conduction. At a stimulation rate of 6 Hz, the C-fiber response was reduced by 60% compared to responses at the rate of 0.75 Hz, and the range of NCV has shifted from 1.5–2 m/s to less than 1.5 m/s.

The complex effects of stimulation rates in the present study support the contribution of a variety of unmyelinated fiber types to the compound C-fiber response. Single unit recordings confirm that high rates of repetitive stimulation are associated with progressive slowing, partial failure of conduction, and an increase in temporal dispersion in C fibers [68, 81, 92, 9799]. The nature and severity of these effects also depend on the class/modality of activated C fibers [54, 84, 99]. For instance, mechano-insensitive C fibers are very sensitive to stimulation rates, while sympathetic efferent fibers demonstrated minimal changes in conduction under the same conditions. When measured as part of a composite whole nerve response, the effects on the various sub-types of C fibers can create complex patterns reflecting both slowing (e.g., redistribution of NCV) and blocking of slow conduction (e.g., reduction in magnitude of response), as it was observed in the compound responses of the current study.

Another difficulty in the assessment of slow conducting axons is the possible interaction of activity from overlapping populations of the slowest subset of Aδ fibers and the fastest of the C fibers. With greater distances, these velocities would be clearly segregated, but when the examining patterns are isolated from the distal extremes of the whole nerve over a few centimeters, activity in these populations can be complex. Aδ fibers demonstrate slowing and partial failure of conduction at much lower frequencies than Aα and Aβ fibers; activity in these fibers comes very close to the margins of C-fiber conduction in skin innervation [77, 92, 100, 101]. The subcomponents of the late responses can be explored using different time periods after HFRS. Relatively fast conducting subpopulations demonstrate minimal suppression following HFRS, while activity in slower C fibers can continue to be diminished for a period of up to 5 min after stimulation. Under the pressure of HFRS, some Aδ fibers may slow down to the C-fiber range of velocities [100]; however, these larger diameter axons can be identified by their faster rate of recovery.

4.1. Advantages and Limitations of the Novel Approach

The advantages of the use of noninvasive measures of C-fiber activity are obvious. The outlined technique can reliably quantify activity in axons conducting at 1-2 m/s, the measures are relatively rapid, they can be repeated in the same nerve segment in longitudinal studies, and the procedure can measure conduction in fast and slow axons in the same experiment. As expected, the compound measure of slow conduction is sensitive to the effects of CAP, and subcomponents of slow conducting fibers can be differentiated by their response to stimulation rate and their fatigue. In contrast to the single unit assessments, the obtained results reflect the activity of the population of C fibers within intact nerves. The measures are parametric and subject to statistical comparison across time and experimental conditions. The limitations of this technique are also clear. The activation of high-threshold C fibers requires intense stimulation which may limit this technique to anesthetized preparation. The composite nature of the C-fiber response confounds the identification of activity in specific subpopulations, although this can be facilitated by use of HFRS and recovery profiles. As measured in this study, the compound response is insensitive to conduction below 1 m/s, but this may be altered by the use of lower stimulation rates and a higher resolution averaging.

The relatively simple evaluation of C-fiber conduction outlined in the study offers a novel tool for examining unique electrophysiological aspects of SFN and may complement structural measures of skin biopsy with functional evaluation of deficits in transmembrane ion distribution, nodal integrity, and energy utilization (mitochondrial function) acquired from the whole nerve response of small diameter nerve fibers. The longitudinal measure of conduction in Aδ and C fibers will bridge the gap between discoveries on the axonal level and their preclinical evaluation by gross behavior (e.g., withdrawal reaction).

Conflict of Interests

The authors report no conflict of interests.

Acknowledgments

The authors would like to thank Shirley Seto and Linda O’Donnell for their assistance in preparation of the paper. This work was supported by NIH Grant R21 RR025896.