Changes in Aβ non-nociceptive primary sensory neurons in a rat model of osteoarthritis pain
© Wu and Henry; licensee BioMed Central Ltd. 2010
Received: 10 May 2010
Accepted: 1 July 2010
Published: 1 July 2010
Pain is a major debilitating factor in osteoarthritis (OA), yet few mechanism-based therapies are available. To address the need to understand underlying mechanisms the aim of the present study was to determine changes in sensory neurons in an animal model of OA pain.
The model displayed typical osteoarthritis pathology characterized by cartilage degeneration in the knee joint and also manifested knee pathophysiology (edema and increased vasculature permeability of the joint) and altered nociception of the affected limb (hind paw tenderness and knee articulation-evoked reduction in the tail flick latency). Neurons included in this report innervated regions throughout the entire hind limb. Aβ-fiber low threshold mechanoreceptors exhibited a slowing of the dynamics of action potential (AP) genesis, including wider AP duration and slower maximum rising rate, and muscle spindle neurons were the most affected subgroup. Only minor AP configuration changes were observed in either C- or Aδ-fiber nociceptors.
Thus, at one month after induction of the OA model Aβ-fiber low threshold mechanoreceptors but not C- or Aδ-fiber nociceptors had undergone changes in electrophysiological properties. If these changes reflect a change in functional role of these neurons in primary afferent sensory processing, then Aβ-fiber non-nociceptive primary sensory neurons may be involved in the pathogenesis of OA pain. Further, it is important to point out that the patterns of the changes we observed are consistent with observations in models of peripheral neuropathy but not models of peripheral inflammation.
Osteoarthritis (OA) is the most common form of arthritis and it is a major health issue. Pain is the major complaint of patients with OA, and pain is reported to be more disabling than loss of joint function . Various therapeutic approaches are available to treat OA pain, yet the effectiveness of existing drug therapies for OA pain is poor, with only moderate relief. A number of proposed mechanisms for OA pain have focussed on changes in knee joint nociceptors and their sensitivity. These mechanisms include activation of sensitized nociceptors in the knee by local inflammation [2–4], bone marrow lesions or micro fractures and increased intra-osseous pressure . However, clinical data on OA suggest widespread changes in the properties of primary and secondary sensory neurons that might include low threshold knee joint afferents as well as neurons innervating tissues outside the joint. For example, most OA patients experience pain from the arthritic joint as well as referred pain from areas remote from the arthritic joint . Further, even months following total hip replacement joint pain still persists in approximately 12% of patients . Moreover, most OA patients also experience loss of proprioception [8–11] and loss of vibrational sense , which are mediated by Aβ-fiber low threshold neurons. This evidence suggests a general change in sensory neurons beyond simply changes in nociceptive neurons.
Changes in the functional properties of primary afferent neurons may be able to initiate these changes but they are largely overlooked as possible origins of the pain of OA, even though such changes have been suggested in other models of chronic pain, including models of inflammatory pain [13–15] and neuropathic pain [15–18]. Inflammatory pain models are associated with changes only in small dorsal root ganglion (DRG) neurons, possibly C- and Aδ-fiber neurons [13, 19], which are usually associated with pain, or nociception. On the other hand, in neuropathic models all neuronal populations appear to be changed, including large Aβ-fiber neurons [20–27].
Therefore, the aim of the current study was to look beyond knee joint nociceptors and to investigate changes in sensory neurons within dorsal root ganglia (DRG) that innervate not only the knee but also neighboring DRG neurons that innervate other areas of the leg in an animal model of OA pain. In particular, we determined the electrophysiological properties of C- and Aδ-fiber nociceptive neurons as well as those of Aβ-fiber low threshold mechanoreceptors (LTMs), comparing these properties in naïve control animals and in OA animals one month after model induction. We report here that significant changes in action potential (AP) configuration were observed only in Aβ-fiber LTMs at a time during model development when the knee joint histopathology, knee pathophysiology and nociceptive responses of the affected limb confirmed that this was an animal model of OA that also exhibited altered sensory function. Only minor changes were observed in the small diameter nociceptive neurons. If changes in physiological properties reflect changes in functional role of these neurons in sensory processing, the present findings do not demonstrate changes in C- or Aδ-fiber nociceptors and therefore possible participation in the pathogenesis of pain in this surgically-induced rat model of OA.
Materials and methods
This study was carried out on female Sprague Dawley rats (180-225 g) obtained from Charles River Inc. (Saint Constant, QC, Canada). All experimental procedures were approved by the McMaster University Animal Review Ethics Board and conform to the Guide to the Care and Use of Laboratory Animals of the Canadian Council of Animal Care, Vols.1 and 2, and experiments adhered to the guidelines of the Committee for Research and Ethical Issues of IASP published in PAIN®, 16(1983) 109-110.
Following model induction, animals were housed for one month before the acute electrophysiological experiment. At one month animals were tested for changes in nociceptive scores, and some animals were selected for histological and further physiological studies. After the end of the acute electrophysiological experiment each animal was euthanized by an overdose of anesthetic.
Induction of the model of OA
Details of the surgical procedure to establish the surgical knee derangement model of OA have been reported previously . In brief, animals were anesthetized with a ketamine-based anesthetic (ketamine, 100 mg/ml; xylazine, 20 mg/ml; and acepromazine, 10 mg/ml). The right side medial meniscus was removed, and the right side anterior cruciate ligament was cut to induce the unilateral OA. After surgery, the animals were given 0.05 ml of the antibiotic Trimel (sulfamethoxazole plus trimethoprim; Novopharm, Toronto, ON, Canada) once per day for 3 consecutive days, and the analgesic buprenorphine hydrochloride (Temgesic, Schering-Plough, Kenilworth, NJ, USA) twice per day for 2 consecutive days. Animals were allowed to survive for one month, as previous work in our group has suggested that typical signs of OA are entrenched by one month after surgery .
Knee joint histopathology
Knee joints from the model animals were processed for histopathological evaluations in order to confirm that the model successfully mimicked OA. The knee joints were harvested and decalcified in 5% formic acid. Histological processing and assessment of tissues were done by Bolder BioPATH Inc. (Boulder, CO, USA). Briefly, knees were trimmed into two approximately equal frontal halves, processed through graded alcohols, and embedded in paraffin. An initial section was cut, and two additional step sections were cut at 150 μm for a total of three sections, which were stained with toluidine blue and evaluated microscopically for cartilage damage, osteophyte formation and the degree of joint instability. Cartilage degeneration in the tibia and femur was scored none to severe using the following criteria described by Janusz et al.: 1 = minimal superficial zone only; 2 = mild extends into the upper middle zone; 3 = moderate well into the middle zone; 4 = marked into the deep zone but not to tidemark; 5 = severe full thickness degeneration to tidemark . Cartilage degeneration scores were measured for the medial and lateral tibia and femur, and all values from all three slides were summed to provide a total joint cartilage degeneration sum.
Osteophytes were scored 1, 2, 3, 4 or 5 for small (< 299 μm), moderate (300-399 μm), large (400-499 μm), very large (500-599 μm), or extremely large (> 600 μm) depending on the size using an ocular micrometer. Medial and lateral osteophyte scores were added to the total joint cartilage sum to derive a total joint score.
The degree of joint stability was subjectively determined based on fibroproliferative and chondrogenic changes in the synovium/collateral ligament as well as transected cruciate area: mild - some fibroplasia with minor proteoglycan deposition, small foci of chondrogenesis, marginal zone proliferation and/or small osteophytes; moderate - diffuse fibroplasia with proteoglycan deposition and larger foci of chondrogenesis, small to medium osteophytes; severe - diffuse severe thickening of synovium and ligaments with proteoglycan deposition and chondrogenesis, generally large osteophytes. The scoring system was as follows: mild = 1, moderate = 2, and severe = 3.
Tissue edema and plasma extravasation of the knee joint
One group of rats was used to study pathophysiological changes associated with the joint capsule in this model, in particular to assay edema and vascular permeability. Briefly, knee joints were dissected and dried in an oven at 60°C for 24 h. Edema in the knee joint was determined as the weight difference of each knee joint after this dehydration procedure.
To evaluate vascular permeability the degree of plasma extravasation was determined. Evans blue dye (VWR, Mississauga, ON, Canada) was injected through the jugular vein. Twenty minutes following Evans blue dye injection, animals were perfused intracardially with 500 ml of physiological saline. Knee joints were harvested and dried as described above, and then were placed in vials each with 3 ml of formamide (Fisher Scientific, Ottawa, ON, Canada) overnight in an oven at 60°C. Twenty-four hours later, fluids in the vials were filtered and evaluated by the absorbance measured by color spectrophotometry (Biochrom Ltd., Cambridge, UK), compared to pure formamide at wavelength 620 nm. The optical density of Evans blue dye in the joint capsule tissue indicated plasma extravasation, and was calculated as follows: ipsilateral absorbance/ipsilateral weight.
von Frey test to determine hind paw mechanical withdrawal threshold
To determine changes in nocifensive behaviors in OA animals von Frey test was conducted. Animals were placed in the testing chamber and allowed to acclimatize for 30 min prior to testing. von Frey filaments from Stoelting (Wood Dale, IL, USA) were applied to the soft tissue of the plantar surface of the hind paw to determine the withdrawal threshold . The response pattern described by Chaplan was adopted to calculate 50% response threshold . The maximum score possible was 15 grams, and the minimum was 0.25 grams.
Effects of repeated flexion and extension of the OA knee on tail flick latency
Previously, we have established that noxious peripheral stimuli alter reaction time in the tail flick test [33–35]. Thus, we applied a similar approach to determine whether repeated flexion and extension of the deranged knee would be a noxious stimulus and therefore alter tail withdrawal latency.
Rats were gently wrapped in clean surgical drapes that covered the entire body to the base of the tail. They were acclimatized to the wrapping for 20-25 min, twice each day over a two-day period prior to surgery. Tail withdrawal latency was then determined on Model 33 tail flick Analgesia meter (IITC, Woodland Hills, CA, USA) at a point 10 cm from the tip that was blackened prior to the test. The intensity of light beam was set so that a baseline reaction time of 8-10 sec was achieved. Once stable baseline readings had been taken, the deranged knee was then articulated with a full extension-flexion mode through the normal plane of motion 20 times over a 30 sec period. Readings in the tail-flick test were then taken again 3 and 6 min after the articulation. The mean of the three baseline responses was taken as 100%. All subsequent responses were normalized as a percentage of the baseline value.
Animal preparation for acute electrophysiological recordings
Full details of the animal preparation and intracellular recordings have been reported previously . Recordings were made from the L4 DRG partly because it is one of the DRGs that contain knee joint afferents. In brief, the L4 dorsal root was cut close to the spinal cord to allow a 12-15 mm length for electrical stimulation. One pair of bipolar platinum stimulating electrodes (FHC, Bowdoinham, ME, USA) was placed beneath the L4 dorsal root.
The animals were ventilated to achieve an end-tidal CO2 concentration around 40 mmHg. Rectal temperature was maintained at approximately 37°C using an in-house servo-controlled infrared heating lamp. Immediately before the start of recording an initial 1 mg/kg dose of pancuronium was given to eliminate muscle tone. The effect of pancuronium was allowed to wear off hourly in order to confirm a surgical level of anesthesia by observing the pupil for dilation and testing for reflex withdrawal from a pinch to a forelimb. Throughout the experiments, supplemental pentobarbital was added hourly to maintain a surgical level of anesthesia.
Recordings were made intracellularly using sharp glass micropipettes with DC resistance around 40-70MΩ. The microelectrode was advanced with an EXFO IW-800 micromanipulator (Montreal, QC, Canada) until a resting membrane potential of at least -40 mV suddenly occurred and an AP could be evoked by stimulation of the dorsal root. Once this occurred the recording was allowed to stabilize over a five min period. Then, the stimulating electrode was used to deliver a single electrical pulse to evoke an AP for analysis. Recordings were made with a Multiclamp 700B amplifier (Molecular Devices, Union City, CA, USA) and digitized on-line via a Digidata 1322A interface (Molecular Devices) with pClamp 9.2 software (Molecular Devices). The first evoked AP in each neuron was used to determine any differences in configuration between control and OA animals. Measurements of the electrophysiological parameters have been demonstrated . These included conduction velocity, resting membrane potential, AP duration at base, AP half width, AP amplitude, AP rise time, AP fall time, maximum rising rate, maximum falling rate, afterhyperpolarization amplitude, 50% afterhyperpolarization recovery time and 80% afterhyperpolarization recovery time. Analysis was done off-line using pClamp 9.2 software.
Classification of DRG neurons
Response properties of neurons to natural stimuli of peripheral receptive fields were identified by various mechanical stimuli, and classified as previously described . The criterion for the classification of C-, Aδ- and Aβ-fiber DRG neurons was mainly based on dorsal root conduction velocities: ≤ 0.8 m/s for C-fiber neurons, 1.5-6.5 m/s for Aδ-fiber neurons and ≥ 6.5 m/s for Aβ-fiber neurons .
The differentiation of nociceptive and non-nociceptive neurons was based specifically on their sensory properties identified during receptive field searching. Those responding to high intensity, potentially noxious, stimuli were classified as nociceptive neurons, whereas those responding to low intensity, innocuous stimuli were classified as non-nociceptive neurons.
Three major factors were considered in grouping Aβ-fiber LTMs: the threshold of activation, the depth of the receptive field and the pattern of adaptation. Non-nociceptive Aβ-fiber neurons were identified as low threshold mechanoreceptors using a soft brush, light pressure with a blunt probe and light manual tap. These neurons included various subtypes, such as guard hair, field hair, Pacinian, glabrous rapidly adapting, slowly adapting types I and II, and muscle spindle types I and II. Guard and field hair neurons were both rapidly adapting cutaneous hair units and are included together. Pacinian and glabrous neurons were both rapidly adapting non-hair neurons, and were named rapidly adapting neurons. Slowly adapting neurons adapted slowly to light tactile stimuli to the cutaneous receptive fields. Muscle spindle neurons were slowly adapting neurons with subcutaneous receptive fields
For C- and Aδ-fiber DRG neurons, only high threshold or unresponsive C- and Aδ-fiber neurons were recorded and included in the current electrophysiological study. High threshold neurons were those that were activated only by high intensity stimuli such as pinch and squeeze applied with a fine forceps, a coarse-toothed forceps or a sharp object such as a syringe needle. Unresponsive neurons were those not excited by any of the non-noxious or noxious mechanical stimuli listed above, and as defined by Lawson et al .
Normality of electrophysiological data was done with the D'Agostino and Pearson omnibus test. Wherever appropriate, Student's t-test or the Mann-Whitney U-test was used for comparisons between OA and control animals in various neuronal subtypes and for various parameters. All statistical tests and graphing were done using Prism 4 software (GraphPad, La Jolla, CA, USA). P-values are indicated in the graphs and P < 0.05 was considered to indicate a statistically significant difference.
A total of 81 animals was used for different purposes in this study, 50 for electrophysiology (26 control and 24 OA animals), 17 for behavioral studies and knee joint pathology (7 control and 10 OA animals) and 14 for knee pathophysiological evaluations (7 control and 7 OA animals).
Histopathological changes in the knee joint
Pathophysiological changes in the knee joint
Extravasation of Evans blue dye, usually taken as a measurement of vascular permeability, was greater in OA rats. The optical density of Evans blue dye in knees from control rats was 0.1 ± 0.01 (N = 7), and was 0.2 ± 0.02 in OA rats (N = 7; P = 0.029; Fig. 2B), suggesting a loss of vascular integrity in the knees from OA animals.
Nocifensive behaviors of the OA model
The latency to withdrawal of the tail in the tail-flick test was also determined in these animals at one month after model induction. Repeated flexion and extension of the knee had no effect on the latency of the tail flick reflex in control rats (97.1 ± 5.18% of the baseline value). However, in OA rats the same manipulation significantly decreased the latency of the tail flick reflex to 64.7 ± 4.91% of the baseline reading (P < 0.001; Fig. 3B).
AP configurations in C- and Aδ-fiber neurons
Acute electrophysiological experiments were run at one month after model induction.
The C-fiber pool was comprised of 24 neurons (5 neurons with an identifiable receptive field) from 14 OA rats and 32 neurons (19 neurons with an identifiable receptive field) from 21 control rats. The Aδ-neuron pool was comprised of 15 neurons (5 neurons with an identifiable receptive field) from 10 OA rats and 18 neurons (9 neurons with an identifiable receptive field) from 15 control rats.
The resting membrane potential in C-fiber neurons was similar in control rats (-58.7 ± 1.76 mV; N = 29) and in OA model rats (-56.7 ± 2.46 mV; N = 25; P = 0.497). However, resting membrane potential in Aδ-fibers in control rats (-65.1 ± 2.04 mV, N = 17) was less depolarized than in OA model rats (-58.7 ± 2.11 mV; N = 15; P = 0.038; Fig. 4B).
AP amplitude was similar in both control and OA rats (86.6 ± 1.49 mV, N = 33 in control C-fiber neurons vs. 84.1 ± 2.64 mV, N = 25 in OA C-fiber neurons; P = 0.387; 74.9 ± 2.48 mV, N = 18 in control Aδ-fiber neurons vs. 77.7 ± 2.85 mV, N = 15 in OA Aδ-fiber neurons; P = 0.474; Fig. 4D).
AP rise time reflects the duration of the depolarization phase of the AP. No significant difference between OA and control animals was found in the AP rise time in either C-fiber neurons or Aδ-fiber neurons. AP rise time in C-fiber neurons was 1.7 ± 0.14 ms in control (N = 33) vs. 2.1 ± 0.41 ms in OA (N = 25; P = 0.47); AP rise time in Aδ-fiber neurons was 0.9 ± 0.05 ms in control (N = 18) vs. 1.1 ± 0.09 ms in OA (N = 15; P = 0.095). Maximum rising rate was derived from the differentiated conversion of the AP curve, which reflects the rate of depolarization over time. The maximum rising rate in C-fiber neurons was similar in control and OA rats (154.5 ± 6.51 mV/ms in control, N = 33 vs. 158.1 ± 9.95 mV/ms in OA, N = 25; P = 0.796). Similar is maximum rising rate in Aδ-fiber neurons, 178.1 ± 11.08 mV/ms in control, N = 18 vs. 171.7 ± 10.97 mV/ms in OA, N = 15 (P = 0.689). The data are shown in (Fig. 4E).
A similar rationale was adopted to determine the dynamics of repolarization, where AP fall time and maximum falling rate were used to measure the dynamics of the repolarization phase. Repolarization of the AP in either C-fibers or Aδ-fibers in OA animals was not different from that of control animals. AP fall time in C-fiber neurons in control rats (2.6 ± 0.19 ms, N = 33) was similar as that in OA rats (3.3 ± 0.62 ms, N = 25; P = 0.594), as was AP fall time in Aδ-fiber neurons (1.3 ± 0.08 ms, N = 18 in control vs. 1.6 ± 0.14 ms, N = 15 in OA animals; P = 0.084). Maximum falling rate in C-fiber neurons in OA animals (65.1 ± 5.27 mV/ms (N = 24) was similar to that in control animals (67.2 ± 4.84 mV/ms; N = 33; P = 0.734). Maximum falling rate in Aδ-fiber neurons in OA animals (90.5 ± 6.35 mV/ms; N = 15) was similar to that in control animals (98.3 ± 5.54 mV/ms; N = 18; P = 0.364). The data are shown in (Fig. 4F).
Nociceptors typically have a longer afterhyperpolarization period than non-nociceptors . Therefore, measurements of the afterhyperpolarization associated parameters, particularly, 80% afterhyperpolarization recovery time, are liable to be compromised by noise signals during recording. In C-fiber neurons, a total of 16 out of 56 neurons (9 in OA and 7 in control groups) had at least one missing value for afterhyperpolarization, 50% afterhyperpolarization recovery time, or 80% afterhyperpolarization recovery time. A similar problem was observed in Aδ-fiber neuron recordings, where 12 out of 33 neurons (6 in OA and 6 in control groups) lacked the full complement of afterhyperpolarization associated readings. Nonetheless, examination of these afterhyperpolarization associated parameters revealed no difference between OA and control rats, irrespective of the duration or the amplitude. Afterhyperpolarization amplitude was 11.6 ± 0.82 mV (N = 30) and 11.9 ± 1.08 mV (N = 22) in C-fiber control and OA neurons, respectively (P = 0.831), and was 12.1 ± 0.84 mV (N = 16) and 10.3 ± 0.99 mV (N = 12) in Aδ-fiber control and OA neurons, respectively (P = 0.191). The 50% afterhyperpolarization recovery time was 11.4 ± 0.93 ms (N = 30 in control) vs. 12.8 ± 1.51 ms (N = 21 in OA) in C-fiber neuron group (P = 0.414), and was 8.8 ± 1.54 ms (N = 16 in control) vs. 10.7 ± 2.28 ms (N = 12 in OA) in Aδ-fiber neuron group (P = 0.491). The 80% afterhyperpolarization recovery time was 27.5 ± 2.37 ms in control C-fiber neurons (N = 24) vs. 30.4 ± 3.16 ms in OA C-fiber neurons (N = 18; P = 0.451), and was 20.4 ± 2.66 ms in control Aδ-fiber neurons (N = 12) vs. (21.1 ± 6.97 ms in OA Aδ-fiber neurons (N = 8; P = 0.512).
AP configurations in Aβ-fiber LTMs
For comparison of Aβ-fiber LTMs, 83 such neurons were recorded from 25 naïve control rats and 79 were recorded from 22 OA rats. In terms of the breakdown of different types of Aβ-fiber LTMs, both groups of animals yielded comparable numbers of each neuronal subtype. For example, guard/field hair neurons were recorded from 14 rats in the control group and 15 rats in the OA group. Similarly, muscle spindle neurons (slowly adapting with subcutaneous receptive field) were recorded from 15 control rats and 20 OA rats.
Representative electrophysiological parameters of control A-fiber LTMs, such as resting membrane potential, AP duration at base, AP amplitude, maximum rising rate and maximum falling rate, were comparable to what has been reported in vivo [13, 14, 37], and also similar to what Ma defined in the low threshold mechanoreceptor category in an in vitro recording that allowed activation of peripheral receptive fields . Electrophysiological parameters of the control neurons were also in the range of what has been reported in medium to large size neurons in vitro [25, 26, 38, 39].
Locations of receptive fields of Aβ-fiber low threshold mechanoreceptors recorded in both the OA rats and the naive control rats
G/F CTL ( N )
G/F OA ( N )
RA CTL ( N )
RA OA ( N )
MS CTL ( N )
MS OA ( N )
SA CTL ( N )
SA OA ( N )
Aβ LTM CTL ( N )
Aβ LTM OA ( N )
In general, the dynamics of AP genesis were slower in the OA animals, particularly in the depolarization phase of the AP. The duration of the AP was longer in A-fiber LTMs in animals following knee derangement. Compared with the control group (1.0 ± 0.03 ms; N = 83), the AP duration at base was significantly wider in the OA group (1.1 ± 0.02 ms; N = 79; P = 0.028; Fig. 5A). AP half width was significantly longer in neurons in OA animals (0.4 ± 0.02 ms, N = 83 in control vs. 0.5 ± 0.01 ms, N = 79 in OA; P = 0.046). In contrast to the control group (0.4 ± 0.01 ms; N = 83), AP rise time was significantly longer in the OA group (0.5 ± 0.01 ms; N = 79; P < 0.001; Fig. 5B). Maximum rising rate was 291.9 ± 9.41 mV/ms in the control group (N = 83), which was significantly faster than in the OA group (254.3 ± 6.22 mV/ms, N = 79; P = 0.001; Fig. 5E).
However, the AP fall time was not significantly different between the control group and the OA group (P = 0.262); readings were 0.6 ± 0.02 ms (N = 83) in the control group and 0.6 ± 0.02 ms (N = 79) in the OA group (Fig. 5C). Maximum falling rate was similar in the control and the OA group, at 173.3 ± 8.77 mV/ms (N = 82) vs. 145.1 ± 4.39 mV/ms (N = 79), respectively (P = 0.114; Fig. 5F). The AP amplitude was not different in the control vs. the OA group (58.9 ± 1.10 mV, N = 83 vs. 57.6 ± 0.98 mV, N = 79, respectively; P = 0.527; Fig. 5D). The remaining parameters, including conduction velocity, resting membrane potential, 50% afterhyperpolarization recovery time and 80% afterhyperpolarization recovery time, were also not different between the two groups (data not shown).
Changes in AP configuration in subgroups of Aβ-fiber LTMs
Further comparison was made between the OA group and the control group for each subset of Aβ-fiber LTMs based on the 4 subsets described above: guard/field hair, rapidly adapting, slowly adapting and muscle spindle neurons. Muscle spindle neurons were the most affected, followed by guard/field hair neurons. Surprisingly, no significant difference was identified between control and OA groups in either the rapidly adapting neurons or the slowly adapting neurons. The relatively small number of slowly adapting neurons may have contributed to the lack of a significant difference between the OA neurons and the control neurons.
In muscle spindle neurons, the slower dynamics of the AP was the most obvious of all of the parameters studied. Compared with 0.8 ± 0.06 ms (N = 23) in the control group, AP duration at base was significantly wider in the OA group (0.9 ± 0.04 ms, N = 24; P = 0.04). The AP rise time was 0.3 ± 0.02 ms in control (N = 23), which is significantly shorter than that in the OA group (0.4 ± 0.02 ms, N = 24; P = 0.002). Correspondingly, the maximum rising rate was significantly decreased in OA (307.2 ± 19.73 mV/ms, N = 23 in control vs. 246.8 ± 11.34 mV/ms, N = 24 in OA; P = 0.01). The remaining parameters were not different between the OA and the control groups.
In guard/field hair neurons, the slowing of the AP rise time in OA rats was the only statistically significant change that related to the duration of the AP (0.4 ± 0.02 ms, N = 24 in control vs. 0.5 ± 0.02 ms, N = 20 in OA; P = 0.038). No other significant difference was identified in the remaining parameters related to the duration of AP, such as AP duration at base, AP half width and AP fall time. The resting membrane potential of neurons in the control neurons (-66.3 ± 1.55 mV, N = 23) was significantly less depolarized than that of the OA group (-61.1 ± 2.15 mV, N = 19; P = 0.047).
In the present in vivo study using intracellular recording techniques in rat DRG neurons, the electrophysiological properties of Aβ-fiber LTMs (non-nociceptive) and C- and Aδ-fiber nociceptive primary sensory neurons were systematically evaluated in a rat model of OA at one month following model induction. This model was confirmed with fully established osteoarthritis characteristics, namely characteristic cartilage degeneration within the knee joint, edema, increased permeability of the knee vasculature and tactile hypersensitivity of the affected lower limb.
Several important observations were made in the electrophysiological studies. There were prominent changes in electrophysiological properties of Aβ-fiber LTMs suggesting a slowing of the dynamics of AP generation, including a wider duration of the AP and a slower maximum rising rate. Importantly, even Aβ-fiber LTMs innervating non-articular structures were affected by the injury initiated in the knee joint. It is not clear what is driving the changes in non-nociceptive neurons or how these neurons are preferentially affected, but the coincidence of changes in nociception and the selective changes in these neurons might imply a role of Aβ-fiber LTMs in the pathogenesis of OA pain.
In contrast, no changes were observed in Aδ- or C-fiber neurons, except a more depolarized Vm in Aδ neurons. This lack of change in the functional properties was surprising, given the changes in these neurons in animal models of inflammatory pain, as discussed below. However, caution should be reminded to consider the alone-standing depolarized Vm in Aδ neurons as a proof of lowered activation threshold in these neurons and consequently as a mechanism of joint pain, as additional evidence of altered AP genesis in these neurons could not be found.
Joint vs. non-joint afferents
Accumulating clinical data and our own observations suggest sensory neuron changes beyond simply changes in knee joint nociceptors, and these may relate to OA pain pathogenesis. This study was initially designed to investigate changes in DRG neurons that can be activated by stimulating knee joint structures, such as knee joint ligaments, muscle attached to the joint and skin covering the joint. Unexpectedly, during pilot studies other sensory neurons within the same DRG, with receptive fields far beyond the knee joint such as foot, seemed to have undergone changes. This observation is consistent with pain referred to other areas beyond the joint, as reported in OA patients [6, 40]. Thus, we felt compelled to change our initial experimental design and to include in this study all sensory neurons regardless of their receptive fields within L4 DRG, which is known as having the largest number of knee joint afferents . One might argue in favor of recording from articular nerve to achieve the highest yield of knee joint nociceptors in order to comment on the role of these neurons. However, mounting evidence has shown that C- or Aδ-fiber joint nociceptors are not necessarily the cause of OA pain. For example, there is limited correlation of the severity of joint pathology with the severity of joint pain [42–44], suggesting that non-articular factors may give rise to OA pain pathogenesis. Further, in approximately 12% of patients, joint pain is not relieved by total joint replacement , and therefore this post-replacement pain cannot be maintained by sensitization or activation of nociceptors in the joint. Our observations are in line with these clinical presentations in OA patients. We report here that the most significant change in properties was in non-knee joint afferents and between Aβ-fiber low threshold non-nociceptive neurons and small diameter C- or Aδ-fiber "pain" neurons.
The neuropathic pattern of affected neuronal types in OA
The pattern of prominent changes in large Aβ-fiber neurons and the lack of change in small C-fiber neurons is commonly reported in neuropathic models of chronic pain [20–26], and is rarely seen in inflammation models of chronic pain [13, 19]. Therefore, we propose that the electrophysiological changes in our OA rats may be associated with a neuropathic etiology that follows model induction. OA rats exhibited signs of movement-evoked pain behavior as suggested by the tail flick results in these animals. Pain behaviors in our OA animals and characteristic movement-evoked pain in OA patients resemble the movement-evoked pain that has been reported in bone cancer pain animals, a complex pain condition that is believed to involve neuropathic mechanisms and large diameter Aβ-fiber neurons [45, 46]. Further, our results and interpretation are consistent with the suggestion of Ivanavicius et al. (2007), who referred to OA pain as having a neuropathic component in view of the mild immune cell infiltration and the poor efficacy of non-steroidal anti-inflammatory drugs , and are consistent with the clinical picture of OA patients who commonly describe their pain with terms typically associated with neuropathic pain processes .
This prompts a comparison of our data with data from other groups investigating changes in DRG neurons in models of peripheral neuropathic pain vs. models of peripheral inflammatory pain. It has been suggested that inflammation and neuropathic etiologies likely affect distinct populations of DRG neurons in various chronic pain models. In superficial inflammation models, for example as induced by injecting complete Freund's adjuvant subcutaneously [13, 19], only Aδ-fiber neurons and C-fiber neurons undergo significant changes in electrophysiological properties, with changes in C-fiber neurons being more severe. In adjuvant-induced joint inflammation, to our knowledge no studies on the properties of DRG neurons are available. However, indirect evidence has suggested that the effects of joint inflammation do not influence the large, non-nociceptive A-fiber neurons, as indicated by the lack of expression of two pain-related peptides, calcitonin gene-related peptide and substance P in those neuronal types [48, 49]. On the contrary, in classic neuropathic models, such as the complete sciatic nerve transection model , the partial sciatic nerve transection model , and the lumbar spinal nerve transection model [21, 23–26], changes in AP configuration in A-fiber neurons characterize changes in primary sensory system, and are common. Although in some studies on neuropathic models [20, 21, 24] changes in C-fiber neurons have been reported, such changes are less prominent than those in A-fiber neurons.
Functional changes in Aβ-fiber LTMs
What prompts us to question the participation of C- or Aδ-fiber neurons in OA pain is the lack of correlation of the various changes in nociception with changes in the function of these neurons. It is obvious that some other mechanisms should account for the changes in nociception, including lowered activation threshold of the hind paw and painful flexion and extension of the affected knee joint. These changes in nociception, mainly mechanical sensitivity, occur along with changes in functional properties of Aβ-fiber LTMs. Studies from other research groups also suggest a possible role of Aβ-fiber LTMs in sensory deficits, such as allodynia [20, 21, 23–25, 50], although detailed mechanisms have not been identified. One possible explanation is that some Aβ-fiber non-nociceptive neurons take up a new role in nociception and begin to convey signals along novel pathways leading to nociception/pain after model induction .
The observed changes in AP configuration in Aβ-fiber non-nociceptor neurons, including wider AP duration, longer AP rise time and slower maximum rising rate, reflect slowed dynamics of depolarization and therefore suggest a change in sodium currents in these neurons, either a functional change or a change in expression. However, the specific ionic mechanisms remain unknown, partly because details of the specific sodium channel composition has not yet been identified in functionally classified sensory neuron subtypes, such as hair, Pacinian, glabrous rapidly adapting or muscle spindle neurons (all examples of Aβ-fiber LTMs). According to a recent paper by Fukuoka et al., large A-fiber neurons are thought to express both TTX-sensitive sodium channels (Nav 1.1, Nav 1.6 and Nav 1.7) and TTX-resistant sodium channels (Nav 1.8, Nav 1.9) . Moreover, after axotomy, 75% of A-fiber neurons re-express the embryonic TTX-sensitive Nav 1.3 channel . Therefore, the significance of our data should also be considered in the context of the possible summation effect of various changes in sodium channels Nav 1.1, Nav 1.3, Nav 1.6, Nav 1.7, Nav 1.8, and Nav 1.9.
The patterns of the changes in the electrophysiological properties of Aβ-fiber LTMs but not in C- or Aδ-fiber neurons are consistent with observations from other laboratories in models of peripheral neuropathy but not models of peripheral inflammation. These changes might reflect a change in functional role of primary afferent sensory processing, which might then constitute a novel mechanism in the pathogenesis of pain at the early phase of OA.
The authors are grateful to Dr. Sally N. Lawson for expert guidance during the earliest stages of this project. We thank Dr. Kiran Yashpal and Ms. Sheila Bouseh for the excellent technical support in behavior tests. We thank Mrs. Chang Ye for help in data input and statistical analysis. This work was generously supported by the Canadian Arthritis Network (CAN), the Canadian Institutes of Health Research (CIHR) and McMaster University. Qi Wu was a CIHR Strategic Training Fellow in Pain: Molecules to Community, and was also supported by CIHR, CAN and Canadian Pain Society.
- Jordan J, Luta G, Renner J, Dragomir A, Hochberg M, Fryer J: Knee pain and knee osteoarthritis severity in self-reported task specific disability: the Johnston County Osteoarthritis Project. J Rheumatol 1997, 24: 1344–1349.PubMedGoogle Scholar
- McDougall JJ, Watkins L, Li Z: Vasoactive intestinal peptide (VIP) is a modulator of joint pain in a rat model of osteoarthritis. Pain 2006, 123: 98–105. 10.1016/j.pain.2006.02.015PubMedView ArticleGoogle Scholar
- Schaible HG, Schmidt RF: Effects of an experimental arthritis on the sensory properties of fine articular afferent units. J Neurophysiol 1985, 54: 1109–1122.PubMedGoogle Scholar
- Schuelert N, McDougall JJ: Electrophysiological evidence that the vasoactive intestinal peptide receptor antagonist VIP6–28 reduces nociception in an animal model of osteoarthritis. Osteoarthritis Cartilage 2006, 14: 1155–1162. 10.1016/j.joca.2006.04.016PubMedView ArticleGoogle Scholar
- Mach DB, Rogers SD, Sabino MC, Luger NM, Schwei MJ, Pomonis JD, Keyser CP, Clohisy DR, Adams DJ, O'Leary P, Mantyh PW: Origins of skeletal pain: sensory and sympathetic innervation of the mouse femur. Neuroscience 2002, 113: 155–166. 10.1016/S0306-4522(02)00165-3PubMedView ArticleGoogle Scholar
- Kean WF, Kean R, Buchanan WW: Osteoarthritis: symptoms, signs and source of pain. Inflammopharmacology 2004, 12: 3–31. 10.1163/156856004773121347PubMedView ArticleGoogle Scholar
- Nikolajsen L, Brandsborg B, Lucht U, Jensen TS, Kehlet H: Chronic pain following total hip arthroplasty: a nationwide questionnaire study. Acta Anaesthesiol Scand 2006, 50: 495–500. 10.1111/j.1399-6576.2006.00976.xPubMedView ArticleGoogle Scholar
- Barrett DS, Cobb AG, Bentley G: Joint proprioception in normal, osteoarthritic and replaced knees. J Bone Joint Surg Br 1991, 73: 53–56.PubMedGoogle Scholar
- Hurley MV, Scott DL, Rees J, Newham DJ: Sensorimotor changes and functional performance in patients with knee osteoarthritis. Ann Rheum Dis 1997, 56: 641–648. 10.1136/ard.56.11.641PubMed CentralPubMedView ArticleGoogle Scholar
- Katz WA: Osteoarthritis: clinical presentations. In Osteoarthritis: diagnosis and medical/surgical management. 3rd edition. Edited by: Moskowitz RW, Howell DS, Altman RD, Buckwalter JA, Goldberg V. W.B. Saunders; 2001:231–238.Google Scholar
- Sharma L: Proprioceptive impairment in knee osteoarthritis. Rheum Dis Clin North Am 1999, 25: 299–314. 10.1016/S0889-857X(05)70069-7PubMedView ArticleGoogle Scholar
- Shakoor N, Agrawal A, Block JA: Reduced lower extremity vibratory perception in osteoarthritis of the knee. Arthritis Rheum 2008, 59: 117–121. 10.1002/art.23241PubMed CentralPubMedView ArticleGoogle Scholar
- Djouhri L, Lawson SN: Changes in somatic action potential shape in guinea-pig nociceptive primary afferent neurones during inflammation in vivo. J Physiol 1999, 520: 565–576. 10.1111/j.1469-7793.1999.t01-1-00565.xPubMed CentralPubMedView ArticleGoogle Scholar
- Djouhri L, Lawson SN: Differences in the size of the somatic action potential overshoot between nociceptive and non-nociceptive dorsal root ganglion neurones in the guinea-pig. Neuroscience 2001, 108: 479–491. 10.1016/S0306-4522(01)00423-7PubMedView ArticleGoogle Scholar
- Djouhri L, Koutsikou S, Fang X, McMullan S, Lawson SN: Spontaneous pain, both neuropathic and inflammatory, is related to frequency of spontaneous firing in intact C-fiber nociceptors. J Neurosci 2006, 26: 1281–1292. 10.1523/JNEUROSCI.3388-05.2006PubMedView ArticleGoogle Scholar
- Michaelis M, Liu X, Janig W: Axotomized and intact muscle afferents but no skin afferents develop ongoing discharges of dorsal root ganglion origin after peripheral nerve lesion. J Neurosci 2000, 20: 2742–2748.PubMedGoogle Scholar
- Newton RA, Bingham S, Davey PD, Medhurst AD, Piercy V, Raval P, Parsons AA, Sanger GJ, Case CP, Lawson SN: Identification of differentially expressed genes in dorsal root ganglia following partial sciatic nerve injury. Neuroscience 2000, 95: 1111–1120. 10.1016/S0306-4522(99)00515-1PubMedView ArticleGoogle Scholar
- Rausch O, Newton RA, Bingham S, Macdonald R, Case CP, Sanger GJ, Lawson SN, Reith AD: Nerve injury-associated kinase: a sterile 20-like protein kinase up-regulated in dorsal root ganglia in a rat model of neuropathic pain. Neuroscience 2000, 101: 767–777. 10.1016/S0306-4522(00)00392-4PubMedView ArticleGoogle Scholar
- Xu GY, Huang LY, Zhao ZQ: Activation of silent mechanoreceptive cat C and Adelta sensory neurons and their substance P expression following peripheral inflammation. J Physiol 2000, 528: 339–348. 10.1111/j.1469-7793.2000.00339.xPubMed CentralPubMedView ArticleGoogle Scholar
- Abdulla FA, Smith PA: Axotomy- and autotomy-induced changes in the excitability of rat dorsal root ganglion neurons. J Neurophysiol 2001, 85: 630–643.PubMedGoogle Scholar
- Kim YI, Na HS, Kim SH, Han HC, Yoon YW, Sung B, Nam HJ, Shin SL, Hong SK: Cell type-specific changes of the membrane properties of peripherally-axotomized dorsal root ganglion neurons in a rat model of neuropathic pain. Neuroscience 1998, 86: 301–309. 10.1016/S0306-4522(98)00022-0PubMedView ArticleGoogle Scholar
- Liu B, Eisenach JC: Hyperexcitability of axotomized and neighboring unaxotomized sensory neurons is reduced days after perineural clonidine at the site of injury. J Neurophysiol 2005, 94: 3159–3167. 10.1152/jn.00623.2005PubMedView ArticleGoogle Scholar
- Liu CN, Wall PD, Ben-Dor E, Michaelis M, Amir R, Devor M: Tactile allodynia in the absence of C-fiber activation: altered firing properties of DRG neurons following spinal nerve injury. Pain 2000, 85: 503–521. 10.1016/S0304-3959(00)00251-7PubMedView ArticleGoogle Scholar
- Ma C, Shu Y, Zheng Z, Chen Y, Yao H, Greenquist KW, White FA, LaMotte RH: Similar electrophysiological changes in axotomized and neighboring intact dorsal root ganglion neurons. J Neurophysiol 2003, 89: 1588–1602. 10.1152/jn.00855.2002PubMedView ArticleGoogle Scholar
- Sapunar D, Ljubkovic M, Lirk P, McCallum JB, Hogan QH: Distinct membrane effects of spinal nerve ligation on injured and adjacent dorsal root ganglion neurons in rats. Anesthesiology 2005, 103: 360–376. 10.1097/00000542-200508000-00020PubMedView ArticleGoogle Scholar
- Stebbing MJ, Eschenfelder S, Habler HJ, Acosta MC, Janig W, McLachlan EM: Changes in the action potential in sensory neurones after peripheral axotomy in vivo. Neuroreport 1999, 10: 201–206. 10.1097/00001756-199902050-00001PubMedView ArticleGoogle Scholar
- Zhang JM, Song XJ, LaMotte RH: Enhanced excitability of sensory neurons in rats with cutaneous hyperalgesia produced by chronic compression of the dorsal root ganglion. J Neurophysiol 1999, 82: 3359–3366.PubMedGoogle Scholar
- Wu Q, Henry JL: Delayed onset of changes in soma action potential genesis in nociceptive A-beta DRG neurons in vivo in a rat model of osteoarthritis. Mol Pain 2009, 5: 57. 10.1186/1744-8069-5-57PubMed CentralPubMedView ArticleGoogle Scholar
- Appleton CT, McErlain DD, Henry JL, Holdsworth DW, Beier F: Molecular and histological analysis of a new rat model of experimental knee osteoarthritis. Ann N Y Acad Sci 2007, 1117: 165–174. 10.1196/annals.1402.022PubMedView ArticleGoogle Scholar
- Janusz MJ, Bendele AM, Brown KK, Taiwo YO, Hsieh L, Heitmeyer SA: Induction of osteoarthritis in the rat by surgical tear of the meniscus: Inhibition of joint damage by a matrix metalloproteinase inhibitor. Osteoarthritis Cartilage 2002, 10: 785–791. 10.1053/joca.2002.0823PubMedView ArticleGoogle Scholar
- Pitcher GM, Ritchie J, Henry JL: Paw withdrawal threshold in the von Frey hair test is influenced by the surface on which the rat stands. J Neurosci Methods 1999, 87: 185–193. 10.1016/S0165-0270(99)00004-7PubMedView ArticleGoogle Scholar
- Chaplan SR, Bach FW, Pogrel JW, Chung JM, Yaksh TL: Quantitative assessment of tactile allodynia in the rat paw. J Neurosci Methods 1994, 53: 55–63. 10.1016/0165-0270(94)90144-9PubMedView ArticleGoogle Scholar
- Cridland RA, Henry JL: Facilitation of the tail-flick reflex by noxious cutaneous stimulation in the rat: antagonism by a substance P analogue. Brain Res 1988, 462: 15–21. 10.1016/0006-8993(88)90579-3PubMedView ArticleGoogle Scholar
- Pitcher GM, Yashpal K, Coderre TJ, Henry JL: Mechanisms underlying antinociception provoked by heterosegmental noxious stimulation in the rat tail-flick test. Neuroscience 1995, 65: 273–281. 10.1016/0306-4522(94)00477-MPubMedView ArticleGoogle Scholar
- Romita VV, Henry JL: Intense peripheral electrical stimulation differentially inhibits tail vs. limb withdrawal reflexes in the rat. Brain Res 1996, 720: 45–53. 10.1016/0006-8993(96)00105-9PubMedView ArticleGoogle Scholar
- Lawson SN, Crepps BA, Perl ER: Relationship of substance P to afferent characteristics of dorsal root ganglion neurones in guinea-pig. J Physiol 1997, 505: 177–191. 10.1111/j.1469-7793.1997.00177.xPubMed CentralPubMedView ArticleGoogle Scholar
- Fang X, McMullan S, Lawson SN, Djouhri L: Electrophysiological differences between nociceptive and non-nociceptive dorsal root ganglion neurones in the rat in vivo. J Physiol 2005, 565: 927–943. 10.1113/jphysiol.2005.086199PubMed CentralPubMedView ArticleGoogle Scholar
- Villiere V, McLachlan EM: Electrophysiological properties of neurons in intact rat dorsal root ganglia classified by conduction velocity and action potential duration. J Neurophysiol 1996, 76: 1924–1941.PubMedGoogle Scholar
- Xu ZQ, Zhang X, Grillner S, Hokfelt T: Electrophysiological studies on rat dorsal root ganglion neurons after peripheral axotomy: changes in responses to neuropeptides. Proc Natl Acad Sci USA 1997, 94: 13262–13266. 10.1073/pnas.94.24.13262PubMed CentralPubMedView ArticleGoogle Scholar
- Hawker GA, Stewart L, French MR, Cibere J, Jordan JM, March L, Suarez-Almazor M, Gooberman-Hill R: Understanding the pain experience in hip and knee osteoarthritis--an OARSI/OMERACT initiative. Osteoarthritis Cartilage 2008, 16: 415–422. 10.1016/j.joca.2007.12.017PubMedView ArticleGoogle Scholar
- Salo PT, Theriault E: Number, distribution and neuropeptide content of rat knee joint afferents. J Anat 1997, 190: 515–522. 10.1046/j.1469-7580.1997.19040515.xPubMed CentralPubMedView ArticleGoogle Scholar
- Hannan MT, Felson DT, Pincus T: Analysis of the discordance between radiographic changes and knee pain in osteoarthritis of the knee. J Rheumatol 2000, 27: 1513–1517.PubMedGoogle Scholar
- Hochberg MC, Lawrence RC, Everett DF, Cornoni-Huntley J: Epidemiologic associations of pain in osteoarthritis of the knee: data from the National Health and Nutrition Examination Survey and the National Health and Nutrition Examination-I Epidemiologic Follow-up Survey. Semin Arthritis Rheum 1989, 18: 4–9. 10.1016/0049-0172(89)90008-5PubMedView ArticleGoogle Scholar
- Lethbridge-Cejku M, Scott WW Jr, Reichle R, Ettinger WH, Zonderman A, Costa P, Plato CC, Tobin JD, Hochberg MC: Association of radiographic features of osteoarthritis of the knee with knee pain: data from the Baltimore Longitudinal Study of Aging. Arthritis Care Res 1995, 8: 182–188. 10.1002/art.1790080311PubMedView ArticleGoogle Scholar
- Mantyh PW: Cancer pain and its impact on diagnosis, survival and quality of life. Nat Rev Neurosci 2006, 7: 797–809. 10.1038/nrn1914PubMedView ArticleGoogle Scholar
- Peters CM, Ghilardi JR, Keyser CP, Kubota K, Lindsay TH, Luger NM, Mach DB, Schwei MJ, Sevcik MA, Mantyh PW: Tumor-induced injury of primary afferent sensory nerve fibers in bone cancer pain. Exp Neurol 2005, 193: 85–100. 10.1016/j.expneurol.2004.11.028PubMedView ArticleGoogle Scholar
- Ivanavicius SP, Ball AD, Heapy CG, Westwood FR, Murray F, Read SJ: Structural pathology in a rodent model of osteoarthritis is associated with neuropathic pain: increased expression of ATF-3 and pharmacological characterisation. Pain 2007, 128: 272–282. 10.1016/j.pain.2006.12.022PubMedView ArticleGoogle Scholar
- Bulling DG, Kelly D, Bond S, McQueen DS, Seckl JR: Adjuvant-induced joint inflammation causes very rapid transcription of beta-preprotachykinin and alpha-CGRP genes in innervating sensory ganglia. J Neurochem 2001, 77: 372–382. 10.1046/j.1471-4159.2001.00175.xPubMedView ArticleGoogle Scholar
- Hanesch U, Pfrommer U, Grubb BD, Heppelmann B, Schaible HG: The proportion of CGRP-immunoreactive and SP-mRNA containing dorsal root ganglion cells is increased by a unilateral inflammation of the ankle joint of the rat. Regul Pept 1993, 46: 202–203. 10.1016/0167-0115(93)90033-5PubMedView ArticleGoogle Scholar
- Devor M: Sodium channels and mechanisms of neuropathic pain. J Pain 2006, 7: S3-S12. 10.1016/j.jpain.2005.09.006PubMedView ArticleGoogle Scholar
- Devor M: Ectopic discharge in Abeta afferents as a source of neuropathic pain. Exp Brain Res 2009, 196: 115–128. 10.1007/s00221-009-1724-6PubMedView ArticleGoogle Scholar
- Fukuoka T, Kobayashi K, Yamanaka H, Obata K, Dai Y, Noguchi K: Comparative study of the distribution of the alpha-subunits of voltage-gated sodium channels in normal and axotomized rat dorsal root ganglion neurons. J Comp Neurol 2008, 510: 188–206. 10.1002/cne.21786PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.