The Absence of Sensory Axon Bifurcation Affects Nociception and Termination Fields of Afferents in the Spinal Cord

Introduction

The primary sensory representation of the body within the central nervous system is based on the intricate innervation patterns of dorsal root ganglion (DRG) neurons into the spinal cord. This projection represents an attractive system to study the branching of axons and has enabled the characterization of a cGMP-dependent signaling cascade essential for the bifurcation of sensory axons. This specific form of neuronal branching proceeds during embryonic development at the so-called dorsal root entry zone (DREZ) where axons of DRG neurons split into ascending and descending stem axons that grow along the lateral margin of the spinal cord (Brown, 1981). Previous investigations showed that the cGMP signaling cascade composed of C-type natriuretic peptide (CNP), the receptor guanylyl cyclase Npr2 (also designated GC-B or Npr-B), and the cGMP-dependent protein kinase I (cGKI, also known as PKGI) are essential for the bifurcation of axons from DRG as well as cranial sensory ganglion neurons (Schmidt et al., 2002, 2007, 2009; Zhao and Ma, 2009; Zhao et al., 2009; Ter-Avetisyan et al., 2014). In the absence of any one of these components sensory axons no longer bifurcate and instead turn either in a rostral or caudal direction. Consistent with these observations is the timing and pattern of localization of CNP in the dorsal spinal cord and Npr2 and cGKI in sensory neurons (Schmidt et al., 2009). A critical missing link of the Npr2-mediated cGMP signaling pathway is the characterization of phosphorylation targets of cGKIα in sensory growth cones that mediate axon bifurcation. Such data might provide mechanistic insights into the machinery for bifurcation. The nitric oxide-sensitive guanylyl cyclases (NO-GCs) are not expressed in embryonic DRG neurons and thus not implicated in sensory axon branching in vivo (Schmidt et al., 2007, 2009). Also, variations in cGMP levels caused by the absence of phosphodiesterase 2A do not interfere with proper bifurcation of sensory axons (Schmidt et al., 2016). Notably, collateral sprouting originating from the stem axons as well as branch formation in the periphery of the body are not affected by the absence of Npr2-mediated cGMP signaling (Schmidt et al., 2007, 2009; Ter-Avetisyan et al., 2014).

Additionally, CNP and Npr2 are also involved in the process of endochondral ossification which is essential for long bone growth. Consequently, biallelic loss-of-function mutations including missense, nonsense, frame-shift mutations, insertions and deletions, and splice site mutations in the human Npr2 gene result in acromesomelic dysplasia type Maroteaux (AMDM; OMIM602875), a skeletal dysplasia with an extremely short and disproportionate stature (Bartels et al., 2004; Potter, 2011; Kuhn, 2016). Whether AMDM patients also reveal bifurcation errors of sensory axons when entering the spinal cord is currently not known and unfortunately, neurological qualities have so far not been characterized on patients with mutations in the Npr2 gene. Similarly to human patients, constitutive Npr2-deficient mice show dwarfism (Chusho et al., 2001; Tamura et al., 2004; Tsuji and Kunieda, 2005). Due to their decreased survival rate at post-weaning stages, Npr2 global mouse knockouts are of limited use for further anatomical, physiological, and behavioral studies. Furthermore, the disproportionate bone growth in Npr2-deficient mice might cause the vertebrate column to squeeze on spinal or cranial nerves which in turn might affect sensation. Unlike Npr2 and its ligand CNP, cGKI is not implicated in long bone growth. However, constitutive cGKI-deficient mice have a number of other deficits including gastrointestinal and cardiovascular impairments which limits their use for investigations on the impact of disturbed axon bifurcation on sensory information processing (Hofmann et al., 2006).

To study the functional consequences of the absence of axon bifurcation in the spinal cord in the absence of other phenotypes that may complicate the interpretation of results, we have generated floxed alleles of Npr2 and cGKI for conditional inactivation of Npr2 or cGKI in DRG neurons at early stages. As in constitutive Npr2 or cGKI knockout mice, DRG axons from mice with a conditional inactivation of Npr2 or cGKI in the neural crest completely lack bifurcations which was associated with altered spinal termination patterns. Behavioral testing indicated that noxious heat sensation as well as nociception induced by chemical irritants was impaired, whereas the behavior driven by mechanical stimulation or balance of body position and motor coordination remained unchanged in mice without sensory axon bifurcations.

Materials and Methods

Mice

Generation of Npr2^fl/fl Mice

Based on the bacterial artificial chromosome clone bMQ331a20 (BioScience) two loxP sites were inserted by standard procedures into the intronic sequences flanking exons 17 and 18 of the murine Npr2 gene to generate the target construct. R1 embryonic stem cells (ESCs) (129X1x129S1) were electroporated and clones that had incorporated the targeting vector into their genome were selected by G418 and analyzed for homologous recombination after digestion with EcoRV or BlnI by Southern blotting using either the 3′ or 5′ probe as indicated in Figure 1A. ESCs were injected into blastocysts from C57Bl/6 mice and chimera that transmitted the floxed Npr2 allele were identified by Southern blotting and by PCR genotyping using the following primer sequences: P1 5-GCCACTTTTGCACCCGGATG-3, P2 5-GTGACGCTGTCGAAGGCCTC-3, P3 5-CCTGCTTTGATGCCATTATCG-3, and P4 5-CTGCAACAACCAAAGCTCAG-3. Crossbreeding with Flpe-deleter mice induced the excision of the neomycin cassette flanked by FRT sites (B6; 129S7-Npr2^tm4(flox)Fgr).

FIGURE 1

Figure 1. Generation of a conditional Npr2-deficient mouse mutant and analysis of sensory axon bifurcation at the DREZ. (A) Targeting strategy for the generation of the floxed Npr2 allele and inactivation by Wnt1-Cre. (B) Southern blotting of the wild type, heterozygotes, and homozygous floxed allele using genomic DNA isolated from liver tissue. (C) PCR genotyping of wild type, heterozygotes and homozygous allele. Upper panel primer P1 and P2 flanking the loxP1 site, lower panel primer P3 and P4 flanking loxP2. (D) Western blotting of membrane fractions from E13.5 DRGs demonstrating the absence of Npr2 protein after recombination with Wnt1-Cre. The heavy chain of clathrin served as loading control. (E) Depiction of wild type, Npr2^LacZ/lacZ, Npr2^fl/fl or Npr^2fl/fl;Wnt1^Cre mice of mature stages (see also Figure S1). (F) Survival rate of Npr2^fl/fl or Npr2^fl/fl;Wnt1^Cre mice at different ages. Numbers of animals inspected are given in parenthesis. (G–J) Naso-anal length and body weight of Npr2^fl/fl or Npr2^fl/fl;Wnt1^Cre males and females at different ages. (K) Analysis of DRG axon bifurcation by DiI tracing in Npr2^fl/fl or Npr2^fl/fl;Wnt1^Cre E12.5 embryos (see also Figure S2). Arrow heads mark bifurcations or turns of sensory axons in wild type or mutants, respectively. Bar, 25 μm. (L) Quantification of bifurcation errors demonstrated the complete absence of DRG axon bifurcation in Npr2^fl/fl;Wnt1^Cre E12.5 embryos. The numbers of counted axons and embryos are indicated below the columns. (M) Bifurcation errors persisted at mature stages (P15) as indicated by the reporter Thy1-YFP-H in the Npr2^fl/fl;Wnt1^Cre;Thy1^YFP−H mouse. Arrow heads mark bifurcations or turns of sensory axons in wild type or mutants, respectively. Bar, 250 μm. (N) Analysis of grip strength revealed indistinguishable maximal muscle strength (mean ± SEM) of the forelimbs between mature Npr2^fl/fl (n = 6 males and 6 females) and Npr2^fl/fl;Wnt1^Cre (n = 8 males and 7 females) of both sexes (paired t-test; p = 0.90 and 0.72 in males and females, respectively).

Genotyping of the floxed cGKI mouse (B6.129-Prkg1^tm2Naw) (Wegener et al., 2002) (RRID:MGI:2668654), Thy1-YFP-H (Tg(Thy1-YFP)H^Jrs/J) (Feng et al., 2000) (RRID:IMSR_JAX:003782), Wnt1-Cre (Tg(Wnt1-Cre)11Rth) (Danielian et al., 1998) (RRID:IMSR_JAX:003829), FLPE-deleter strain (B6.129S4-Gt(ROSA)26Sor^tm1(FLP1)Dym/RainJ; The Jackson Laboratories) (Farley et al., 2000) (RRID:IMSR_JAX:009086), TrkA-Cre mice (B6;129S4-Ntrk1^tm1(cre)Lfr/Mmucd; Mutant Mouse Regional Resource Centers, https://www.mmrrc.org/) (Quina et al., 2015) (RRID:MMRRC_015500-UCD), TrkC-Cre mice (B6.129X1-Ntrk3^tm1(cre)Lfr/Mmucd, originating in the laboratory of Louis Reichardt, provided by Karina Gültig, University of Tübingen) (Funfschilling et al., 2004) (RRID:MMRRC_000364-UCD), Npr2-LacZ (B6.129P2-Npr2^{tm1.1(nlslacZ)/Fgr}) (Ter-Avetisyan et al., 2014) (RRID:MGI:5568090), Npr2-CreERT2 (B6.129S7-Npr2^{tm1.2(CreERT2)/Fgr}) (Ter-Avetisyan et al., 2014) or Rosa26-tdTomato (B6.Cg-Gt(ROSA)26Sor^{tm14(CAG−TdTomato)Hze}/J) (Madisen et al., 2010) (RRID:IMSR_JAX:007914) was performed by PCR as described.

Animals were housed on a 12/12 h light/dark cycle with free access to food. The animal procedures were performed according to the guidelines from directive 2010/63/EU of the European Parliament on the protection of animals used for scientific purposes. All experiments were approved by the local authorities of Berlin(LaGeSO) (numbers T0313/97, 0143/07, G0370/13, G0239/11, X9014/15, and G0222/14), the Landesamt für Natur, Umwelt und Verbraucherschutz Nordrhein-Westfahlen (number 84-02.04.2012.A422), and Regierungspräsidium Darmstadt (number FU/1102). All experiments were also approved by the local Ethics Committee for Animal Research and adhered to the guidelines of the Committee for Research and Ethical Issues of the International Association for the Study of Pain.

Axon Tracing by DiI, Immunohistochemistry, and Western Blotting

DiI tracing and immunofluorescent staining of cryostat sections (16 μm) from paraformaldehyde-fixed tissue were performed as described previously (Schmidt et al., 2007; Schmidt and Rathjen, 2011).

The following primary and secondary antibodies using the indicated concentrations or dilutions were applied to tissue sections or Western blotting: Guinea pig antiserum to the extracellular domain of Npr2 (dilution: 1:5,000) (Ter-Avetisyan et al., 2014), guinea pig antiserum to amino acids 2-89 of cGKIα (1:25,000 for immunofluorescence staining), affinity-purified rabbit antibody to amino acids 2-89 of cGKIα (0.25 μg/ml for Western blot analysis), rabbit antibody to full length cGKIα expressed in Sf9 cells (1 μg/ml of the IgG fraction for Western blot analysis), chicken anti-β-galactosidase (1:5,000; Abcam, ab 9361; RRID:AB_307210), mouse anti-clathrin heavy chain (0.05 μg/ml; BD Biosciences, 610499; RRID:AB_397865), mouse anti-GAPDH (1:7,500; Novus Biologicals, NB300-221; RRID:AB_10077627), rabbit anti-PGP9.5 (1:1,000, Dako, Z5116, RRID:AB_2622233), rabbit anti-trkA (1:1,000; Millipore, 06-574; RRID:AB_310180), rabbit anti-parvalbumin (1:1,000; Swant, PV28; RRID:AB_2315235), rabbit anti-CGRP (1:1,000; Chemicon, AB1971; RRID:AB_2313629), rabbit anti-red fluorescent protein (1:2,500, ABIN129578; RRID:AB_10781500), isolectin GS-IB4-Alexa-488 (1:50; Life Technologies; RRID:AB_2314662), donkey anti-chick-IgY-Cy3 (1:1,000; Dianova), goat-anti-rabbit-Alexa488 (1:1,000; Dianova), goat anti-guinea pig-Alexa488 (1:1,000; Dianova), donkey anti-guinea pig-HRP (Dianova), goat anti-rabbit-HRP, and goat anti-mouse-HRP (all 1:20,000; Dianova,). 0.1 mg per g body weight of tamoxifen was applied by oral gavaging as described (Ter-Avetisyan et al., 2014).

Microscopic images were obtained at room temperature by confocal imaging using a Carl Zeiss LSM 710 NLO Laser Scanning Microscope equipped with ZEN 2010 software and the following lenses: a Plan-Neofluar 10x/0.30 NA objective, a Plan-Achromat 40x/1.40 NA oil objective, or a Plan-Achromat 63x/1.40 NA oil objective (all from Carl Zeiss MicroImaging, GmbH). Images were imported into Photoshop CS5 (Adobe) for uniform adjustment of contrast and brightness. Figures were assembled using Illustrator CS5 (Adobe).

Transganglionic Labeling Using Cholera Toxin Subunit B, Optical Clearing, and Imaging of Fixed Spinal Cord Tissue

Five to six-week old mice of either sex were anesthetized by an intraperitoneal injection of ketamine (100 mg/kg) and xylazine (10 mg/kg). 0.2 μl of 1.5% cholera toxin subunit B conjugated with Alexa Fluor 594 (CTB-AL594) (Wan et al., 1982; Robertson and Arvidsson, 1985; Conte et al., 2009) in 0.1 M phosphate buffered saline (PBS) were injected subcutaneously into the plantar surface of the left second hindpaw digit using a pulled glass capillary attached to a Hamilton microliter syringe. The glass capillary was inserted into the most distal interphalangeal crease and advanced under the skin toward the next proximal crease where the tracer was slowly injected. Five days post-injection, allowing for transganglionic transport of the tracer, mice were transcardially perfused with PBS and ice-cold 4% paraformaldehyde (PFA). Subsequently, the lumbar spinal cord was dissected out and post-fixed overnight in 4% PFA at 4°C. Fixed spinal cords were washed three times with 0.1 M PBS for 10 min each at RT. Subsequently, the tissue was immersed in ascending concentration series of 2,2'-thiodiethanol (TDE) (Staudt et al., 2007; Kloepper et al., 2010; Aoyagi et al., 2015; Costantini et al., 2015) for 24 h each at RT. The applied concentrations were 10, 25, 50, and 97% TDE diluted with 0.1 M PBS. During all incubation steps, the samples were kept on a vibrating table in the dark. For two-photon imaging, the cleared spinal cords were mounted on glass slides in 97% TDE which has a refractive index of 1.52 (Staudt et al., 2007) using press-to-seal silicone isolators with the dorsal surface facing up.

Two-photon imaging was performed using a laser scanning microscope equipped with a tunable Ti:sapphire laser. Two channels were recorded sequentially to collect Alexa Fluor 594 fluorescence (excitation wavelength: 810 nm; emission range: 600–735 nm) and tissue autofluorescence (excitation wavelength: 810 nm; emission range: 505–575 nm). A 25x multi-immersion objective (0.8 numerical aperture) was used with immersion oil for cleared sample imaging. Tiled stacks were taken through the spinal cord dorsal horn (pixel size: 0.5 × 0.5 μm; z step size: 1.5 μm). All images were processed using ImageJ (Schneider et al., 2012). Tiled stacks were stitched using either the imaging software ZEN 2010 or the ImageJ plugin “Stitching 2D/3D” (Preibisch et al., 2009). Subsequently, images were cropped to the same size and reduced to the same slice number. Next, background fluorescence was reduced by subtracting the autofluorescence channel from the CTB-AL594 channel. Using stack histogram-based thresholding, the image stacks were binarized. The threshold was set as the mean gray value plus three times the standard deviation. Finally, single pixels were removed to reduce noise, e.g., hot pixels. In order to enable comparative analyses of spinal terminal fields of cutaneous myelinated afferents, the three-dimensional centers of mass of the voxel clouds representing CTB-AL594 labeled fiber terminals were determined using the ImageJ plugin “3D ImageJ Suite” (Ollion et al., 2013). All images were aligned to the center of mass of the voxel cloud, i.e., images were cropped to the same size and reduced to the same slice number around the respective centers of mass. Summed dorsoventral, rostrocaudal, and/or mediolateral projections of the binary image stacks were constructed to enable two-dimensional visualization of terminal fields.

Morphometric Analysis of Spinal Terminal Fields

Mediolateral, rostrocaudal, and dorsoventral spans of the terminal fields were measured in summed dorsoventral and rostrocaudal projections of binary image stacks. Summed projections were thresholded with the threshold being set as the mean gray value plus one standard deviation. Subsequently, the dimensions of the bounding rectangles enclosing all pixels representing CTB-AL594 labeled fiber terminals were measured.

Density Analysis of Spinal Terminal Fields

Areal densities (as voxels per area) of spinal terminal fields were calculated in summed dorsoventral, mediolateral, and rostrocaudal projections of binary image stacks. Using the ImageJ plugin “3D ImageJ Suite” (Ollion et al., 2013), the total number of voxels representing CTB-AL594 labeled fiber terminals was determined in binary image stacks. Subsequently, the number of voxels was divided by the area (in μm²) that was occupied by positive pixels in summed dorsoventral, mediolateral, and rostrocaudal projections, respectively.

Ex Vivo Skin Nerve Preparation

Ex vivo skin nerve preparations were performed as described (Moshourab et al., 2013). Briefly, mice were sacrificed and the glabrous skin on the left hindlimb was removed. The saphenous nerve was exposed dissected free along the lower leg up. Subsequently, the skin was carefully removed from the musculoskeletal and the connective tissue of the paw. The skin-nerve preparation was placed in an organ bath filled with oxygenated 32°C warm synthetic interstitial fluid (SIF) consisting of 123 mM NaCl, 3.5 mM KCl, 0.7 mM MgSO₄, 1.7 mM NaH₂PO₄, 2.0 mM CaCl₂, 9.5 mM sodium gluconate, 5.5 mM glucose, 7.5 mM sucrose and 10 mM HEPES, at a pH of 7.4. Using insect needles, the skin was mounted in the organ bath with its epidermis facing the bottom of the chamber, exposing the dermis to the solution. The nerve was pulled through a hole into the adjacent recording chamber which was filled with mineral oil. Finally, using fine forceps the nerve was desheathed by removing its epineurium and small filaments were teased of the nerve. Throughout the whole experiment the skin was superfused with oxygenated SIF at a flow rate of 15 ml/min.

Single-Unit Recordings

Single receptor units were recorded as previously described (Moshourab et al., 2013). Teased filaments were attached to a recording electrode and the receptive fields (RF) of individual units were identified by manually probing the skin with blunt forceps. To probe for heat-responsive units, 1 ml of hot (48°C) SIF buffer was washed over the surface of the skin. For immediate visual identification of single units, whole action potential waveforms were resolved on an oscilloscope. Data was acquired using a PowerLab 4/30 system which was controlled with the software LabChart 7.1.The conduction velocities (CVs) of single fibers were determined by evoking a local action potential with a platinum iridium electrode. The electrical impulse was conducted nearly instantaneously through the solution whereas the triggered action potential conducted by the fiber was delayed depending on the fiber type and the distance of the fiber's receptive field from the electrode. Hence, the distance between the RF of a unit to the electrode was measured and the CV was calculated as distance divided by time delay. Fibers with CVs below 1.3 m/s were classified as C-fibers.

Stimulation Protocols

A Peltier device (custom device built by the Yale School of medicine Instrument Repair and Design) was used to apply warm or cool stimulation parameters to the RFs of cutaneous afferent fibers. Two heat ramps were applied to skin: firstly a ramp and hold stimulus from 32 to 48°C with a 0.5 s ramp, followed by a 3 s hold time and a further 0.5 s ramp; and secondly a slow 15 s ramp from 32 to 48°C. Two cold ramps were also applied: first a ramp and hold stimulus from 32 to 12°C with a 0.5 s ramp, followed by a 3 s hold time, and a further 0.5 s ramp; and secondly a slow 15 s ramp from 32 to 12°C.

A computer-controlled Nanomotor (Kleindieck, Germany) was used to apply controlled mechanical stimuli of known amplitude and velocity to the RF of afferent fibers. Four 3 s ramp and hold stimuli were applied to the RF of each afferent fiber, each increasing in amplitude (40, 75, 150, 220 mN).

Analysis of Extracellular Afferent Recordings

Thermosensitive C-fibers were searched for in recordings and for further analysis. Spike sorting was performed with a spike analysis plug-in of the LabChart software (AD Instruments). C-mechanoheat (C-MH), C-mechanoheatcold (C-MHC), and C-mechanocold (C-MC) fibers were pooled in each experimental group (Milenkovic et al., 2014). Thermal thresholds represent the temperature required to cause the first action potential spike in slow (15 s) heat and cold ramp protocols. Mechanical thresholds represent the average force required to cause the first action potential spike over the four ramp and hold stimuli. Mean group data were compared using t-tests and two-way repeated measures ANOVA tests.

Behavioral Testing

Littermate mice of either sex were used in all behavioral tests (8 to 10-weeks old). The Npr2^fl/fl;Wnt1^Cre mice were on a mixed genetic background (Bl6/SV129) and cGKI^KO/fl;Wnt1^Cre mice were on a Bl6 background. Animals were habituated to the experimental room and were investigated by observers blinded for the genotype.

Hot Plate Test

Mice were placed into a Plexiglas cylinder on a metal surface maintained at 50, 52, or 54°C (Hot Plate; Ugo Basile, Comerio Italy). Cut-off times were 60, 40, and 20 s, respectively, to prevent tissue damage. The time between placement and shaking or licking of the hindpaws or jumping off the plate was recorded.

Acetone Test

For acetone-evoked evaporative cooling, animals were habituated in a glass enclosure on a mesh floor, a drop of acetone (50 μl) was applied to the plantar side of a hindpaw and responses were observed for 1 min. The responses were scored as follows: 0, no response (sniffing or walking was not considered nociceptive); 0.5, a licking response; 1, flinching and brushing of the paw; 2, strong flinching (Caspani et al., 2009).

Capsaicin Test

Capsaicin (5 μg dissolved in 10 μl of a solution consisting of 90% sterile saline, 5% ethanol, and 5% Tween 80) was injected subcutaneously (s.c.) into the dorsal surface of a hindpaw. The time spent licking the capsaicin-injected paw and the number of shakes of the injected paw was recorded for 5 min.

Formalin Test

Formalin (20 μl of a 0.5% formaldehyde solution) was s.c. injected into the dorsal surface of a hindpaw (Braz and Basbaum, 2010). The time spent licking the formalin-injected paw was recorded in 5 min intervals up to 60 min after formalin injection.

Hargreaves Test

The Hargreaves test was performed with an idle value of 5% and a stimulus value of 50% of the maximum intensity as described (Hargreaves et al., 1988). The withdrawal reaction induced by heat stimulation was measured 25 times per animal over period of 3 consecutive days using model 400 Heated Base (IITC Life Science Inc.).

Dynamic Plantar Test

The mechanical sensitivity of the plantar side of a hindpaw was assessed with an automated von Frey-type testing device (Dynamic Plantar Aesthesiometer; Ugo Basile). This device pushes a thin probe (0.5 mm diameter) with increasing force through a wire-grated floor against the plantar surface of the paw from beneath, and it automatically stops and records the latency time, after which the animal withdraws the paw. The force increased from 0 to 5 g within 10 s (0.5 g/s ramp) and was then held at 5 g for an additional 10 s (Lu et al., 2015). The paw withdrawal latency was calculated as the mean of five to six consecutive trials with at least 20 s in-between.

Balance and Motor Coordination in the Absence of Npr2

Grip Strength Assay

Motor functions were assessed using the grip strength assay that measures maximal muscle strength of the forelimbs (TSE Grip Strength Meter). An animal holding onto the meter handle was pulled backwards steadily in horizontal orientation, and the mean of three maximum force measurements was recorded.

Balance and coordination were examined using the beam walking test, rotarod performance test, staircase test, adhesion removal test, food grasping and reaching test, and walking track (MouseWalker) which provide readouts for general motor performance.

The Beam Walking Test

The beam walking test assesses motor coordination and balance from the ability of mice to traverse a series of narrow, horizontal beams graded in diameter. Wooden beams were 1 m × 10 or 15 or 20 mm diameter at 50 cm height with one starting end (with a 60 W aversive light) and the other end attached to an escape box. Mice were trained for 2 days with four trials per day, motivated by a reward in the escape box until traversing took < 20 s. On the test day (day 3) mice received two consecutive trials per beam, progressing from wide to narrow, allowing up to 60 s for a traversal. A fall off counted as incomplete trial. Durations of beam traversals und instances of hind feet slipping off the beam were recorded (Carter et al., 2001; Luong et al., 2011).

The Staircase Test

The staircase test (Campden Instruments) allows measurement of lateralized effects on skilled motor function required for reaching and grasping a food pellet from variable distance. Food was removed from subjects' home cages 4–6 h before first training and prior to later experiments animals were fed one pellet of regular feed. Animals were habituated to sucrose pellet rewards at least three times prior to testing to avoid food neophobia. Mice were trained for 3 days in the staircase test box with the double, eight-step staircase baited with two pellets per step (16 on each side, total 32). Testing was conducted with one pellet per step and sessions lasted for 15 min. We recorded collected pellets, and maximum distance reached.

The Adhesion Removal Test

The adhesion removal test assesses sensory and motor deficits related to paw and mouth coordination. It was performed in an empty mouse cage where a mouse received two adhesive tape strips (3 × 4 mm) onto the hairless part of the forepaws. The time to first contact reaction was recorded for both paws, as well as the time to tape removal.

The Reaching and Grasping Test

The reaching and grasping test was performed with an automated forepaw reaching chamber (Campden Instruments) providing sensitive measures of fine motor coordination (Marques and Olsson, 2010). Prior to testing subjects were habituated to sucrose pellets. Four days before testing, animals were trained daily to reach for pellets, each training lasting until either 10 pellets were retrieved or 10 min had passed. During the test one pellet at a time was placed in the tray. Reaching attempts were counted until the pellet was removed. Performance parameters were latency to first reach, reaching accuracy (pellets retrieved per reach attempt), latency to first pellet retrieval, interval to retrieve five pellets after the first pellet had been retrieved.

The Rotarod Test

The Rotarod test was performed using Rota-Rod 47600 (Ugo Basile) with an acceleration of the rotating cylinder from 5 to 40 rpm in 360 s. Animals were given 2 days of habituation to the Rotarod (same protocol as in test). Four recordings per day were done.

The mouse's gait was tested using a home-made MouseWalker as described in Mendes et al. (2015). This system is based on the total internal reflection of light on transparent plates. Paw contacts of the mice disrupt this reflection and generate scattered light which is then detected by a high speed digital camera. Parameters such as foot print clustering and stance linearity were evaluated.

Experimental Design and Statistical Analysis

Sample size for behavioral experiments and for analysis of axon branching was deduced from previously published studies (Schmidt et al., 2007) or in publications listed above in the sub-section on behavioral testing. No further statistical methods were used to predetermine sample size. Sample size are given directly in the figures or in the legends. Experiments were done blind with respect to genotype. Littermates were used as control. Statistical analysis was performed with SPSS software (RRID:SCR_002865) using the Student's t-test for paired comparisons or Mann-Whitney U-test. Power analysis of statistical data was performed with the software “nQuery + nTerim 4.0.” Depending on the equality or inequality of variances and on the equality and inequality of the number of measurements t-test or Satterthwaite t-test were applied.

Results

Conditional Inactivation of Npr2 in DRG Neurons Interrupted Sensory Axon Bifurcation

To generate a conditional allele of Npr2, two loxP sites were introduced into the introns flanking exons 17 and 18 of the mouse Npr2 gene which encode the guanylyl cyclase domain of Npr2. The neo selection cassette flanked by two FRT sites was removed by crossing this mouse line with a FLPE-deleter strain. The resulting floxed Npr2 allele contains two loxP sites surrounding exon 17 and 18 (Figure 1A). The correct integration of the loxP sites was verified by Southern blotting using (Figure 1B) a probe as indicated in Figure 1A and PCR was applied for genotyping using two distinct primer pairs that amplify genomic regions including loxP1 or loxP2, respectively (Figure 1C). Conditional inactivation of Npr2 in DRG neurons at early embryonic stages when sensory axons enter the spinal cord was obtained by crossbreeding with the Wnt1-Cre transgenic mouse line. At early embryonic stages Wnt1 is restricted to the midbrain, dorsal spinal cord and DRGs as detected by the Wnt1-lacZ reporter mouse (Echelard et al., 1994). Cre-mediated excision of the exon 17 and exon 18 generated a frameshift in exon 19 which leads to a premature stop codon. This resulted in nearly undetectable expression of Npr2 protein in DRGs as indicated by Western blotting using antibodies to the extracellular domain of Npr2 (Figure 1D).

About 50% of these Wnt1-Cre inactivated Npr2 mutants developed normally as demonstrated by a nasal to anal body length and a body weight that were indistinguishable from wild type or the homozygous floxed strain (Figures 1E–J). Furthermore, grip strength in these mutant mice appeared normal (Figure 1N). However, about 50% of Npr2^fl/fl;Wnt1^Cre mice were retarded in their development which revealed also a reduced survival rate at post-weaning stages (Figure 1F). Premature death of these mutants was accompanied by a severe malocclusion of the upper and lower incisors and intestinal distension which both affected normal nutrition (Figure S1). These mice were excluded from further testing. A malocclusion in conditional mutants of the Npr2 ligand CNP (Nakao et al., 2015) and milk retention in the stomach and intestinal distension have also been observed in the Npr2 mutant slw/slw (Sogawa et al., 2010).

Analysis of sensory axon bifurcation by DiI tracing at embryonic day 12.5 or by the reporter line Thy1-YFP-H at postnatal day P15 indicated a complete absence of axon bifurcation at the dorsal root entry zone of the spinal cord in the conditional mutant but not in the homozygous floxed allele (Figures 1K–M). DRG axons turn either in a rostral or caudal direction at the lateral margin of the cord. In contrast, the formation of collaterals from stem sensory axons in the dorsal funiculus is not impaired in the absence of Npr2 (data not shown) indicating that this type of branch formation is not dependent on Npr2 signaling as described previously (Schmidt et al., 2007). In contrast to conditional inactivation of Npr2 by Wnt1-Cre the use of other Cre-lines such as trkA-Cre or trkC-Cre did not induce bifurcation errors most likely due to later initiation of Cre expression with respect to Npr2 in embryonic DRG neurons (Figure S2). Overall, our data indicate that the Npr2^fl/fl;Wnt1^Cre mutant is a suitable model to study the consequences of impaired axon bifurcation in adult mice.

Receptive Fields in the Spinal Cord Are Altered in the Absence of Sensory Axon Bifurcation

Sensory information from a large number of afferent axons converges within the spinal cord in nociceptive, mechanoreceptive or proprioceptive fields by the generation of collaterals that terminate in specific layers. Bifurcation of sensory axons at the DREZ contributes to this overall representational organization (Brown, 1981). The normal body length of the conditional Npr2 mutant allowed us to clarify whether the absence of bifurcation affects the size or shape of termination fields of afferents from the skin. A lack of bifurcation at the DREZ predicts altered receptive fields and less overlap between different sensory modalities. Therefore, we mapped the spinal terminations of cutaneous myelinated afferents from the hind paw in wild type and Npr2^fl/fl;Wnt1^Cre mice. The topography of primary afferent termination was visualized by transganglionic transport of fluorescently labeled cholera toxin B (CTB) subcutaneously injected into the second digit of the left hind paw of 5-week old mice (Figure 2A). CTB binds to the ganglioside GM1 which is highly expressed on axons (Kobbert et al., 2000). Injection in digit two resulted in a characteristic fork-like afferent fiber termination field in control Npr2^fl/fl spinal cords (Figure 2B). Total number of positive voxels representing CTB-labeled fibers as well as spans and area densities in dorso-ventral, medio-lateral, and rostro-caudal planes were calculated (Figures 2C–I). The total number of voxels was reduced to 64% for the 2nd digit (Figure 2C) indicating that the number of incoming sensory axons and collaterals were strongly reduced in comparison to wild type - surprisingly not to 50% as one might have been expected from the total lack of bifurcation at the DREZ. This observation might be explained by an increase in the density of collaterals protruding from the stem axons or by an increase in the terminal branching of collaterals in the mutant. Interestingly, the dorsoventral span of fluorescence intensity was increased to about 131% of control for the 2nd digit in the mutant whereas the medial-to-lateral distribution of terminal field was decreased by 25%. No difference was observed in the rostral-to-caudal span of the projection (Figures 2D–F). The density, the voxels per area, decreased to 66 and 58% in the rostrocaudal and in the mediolateral direction, respectively, whereas the dorsoventral direction decreased only to 76% (Figures 2G–I) for digit two. In summary, the absence of bifurcation caused decreased projection intensities together with a reduction in the mediolateral span as well as an increase in the dorsoventral extension of sensory afferent termination fields in the spinal cord.

FIGURE 2

Figure 2. Receptive fields of spinal terminations of tactile afferents from the hindpaw in wild type and Npr2-deficient conditional mice. (A) Subcutaneous injection of 0.3 μl of CTB conjugated to AL594 into the second digit of the left hind paw in order to label the terminal fields of cutaneous myelinated afferents in the spinal cord dorsal horn. A fluorescent micrograph of the left hind paw is shown. Scale bar, 1 mm. (B) Heat maps of averaged summed projections of spinal terminal fields in Npr2^fl/fl;Wnt1^Cre and Npr2^fl/fl control mice. Terminal fields of fibers innervating the second digit of the left hind paw appear different in Npr2^fl/fl;Wnt1^Cre mice (n = 6) than in Npr2^fl/fl mice (n = 6) in dorsoventral, mediolateral, and rostrocaudal projections. The ImageJ color lookup table “16 Colors” was applied (see color map). The color code in each pixel denotes the number of voxels found in corresponding positions along different axes averaged across mice. Scale bars, 100 μm. (C) In Npr2^fl/fl;Wnt1^Cre mice, the number of voxels representing CTB-labeled fibers is significantly lower compared to those in Npr2^fl/fl mice. (D–F) Morphometrics of spinal terminal fields in Npr2^fl/fl;Wnt1^Cre and Npr2^fl/fl control mice. Spinal terminal fields of fibers innervating the second digit of the left hind paw span similarly far in the (D) RC dimension but in a significantly different range in the (E) ML and further in the (F) DV dimension in Npr2^fl/fl;Wnt1^Cre (red, n = 6) as compared to Npr2^fl/fl (blue, n = 6) mice. (G–I) Areal densities of spinal terminal fields in Npr2^fl/fl and Npr2^fl/fl;Wnt1^Cre mice. The areal densities of terminal fields innervating the second digit of the left hind paw digit were significantly lower in (G) the dorsoventral projection, (H) the mediolateral projection, and (I) the rostrocaudal projection in Npr2^fl/fl;Wnt1^Cre (red, n = 6) as compared to Npr2^fl/fl (blue, n = 6) mice. Individual data points and mean values ± SEM are shown. Each data set was compared using a two-tailed unpaired t-test (*p < 0.05; **p < 0.01; ***p < 0.001). CTB, cholera toxin subunit B; WT, wild type; M, medial; V, ventral; C, caudal; D, dorsal; RC, rostrocaudal; ML, mediolateral; DV, dorsoventral.

Balance and Motor Coordination Is Not Impaired in the Absence of Sensory Axon Bifurcation

Since axon bifurcation is compromised in all DRG neurons of Npr2-deficient mice we asked whether proprioception might be impaired. Therefore, balance and motor coordination in Npr2^fl/fl;Wnt1^Cre mice was examined using the balance beam test, rotarod, staircase assay, food grasping and reaching assay, and walking track analysis (MouseWalker). Balance beam and rotarod measure the ability to perform complex and coordinated movements. Npr2-deficient mice did not show a reduced latency to fall off or to slip off the beam (Figures 3A,B) or to fall off the rotating rod (Figure 3C). Similarly, no deficits were observed in the staircase and reaching and grasping test (Figures 3E–I), indicating that motor coordination was not altered in the absence of sensory axon bifurcation. To assess sensory and motor abilities related to the forepaw and mouth coordination the adhesion removal test was performed. Adhesive tape was fixed gently to the hairless parts of the forepaws and the time to first contact and the removal of both adhesive strips was recorded. Control mice and Npr2 mutants made regular attempts to remove the tape (Figure 3D). To ask whether the absence of sensory axon bifurcation interferes with the control of locomotion we applied the MouseWalker system which allows a comprehensive and quantitative description of parameters of freely walking Npr2-deficient mice (Figures 3J–L).

FIGURE 3

Figure 3. Analysis of balance and motor coordination in the absence of Npr2^fl/fl;Wnt1^Cre. Balance and motor coordination in Npr2-deficient mice was not altered in the absence of DRG axon bifurcation using the (A) beam walking test (n = 13 Npr2^fl/fl and 14 Npr2^fl/fl;Wnt1^Cre), passage time and (B) number of slips, (C) rotarod performance test (n = 5), (D) tape removal test (n = 13 Npr2^fl/fl and 14 Npr2^fl/fl;Wnt1^Cre), (E) staircase test, no. of pellets and (F) reaching distance (n = 13 Npr2^fl/fl and 14 Npr2^fl/fl;Wnt1^Cre), (G) food grasping and reaching, number of attempts, (H) number of pellets and (I) successful attempts (n = 10 Npr2^fl/fl and 15 Npr2^fl/fl;Wnt1^Cre), (J) walking track analysis using the MouseWalker, (K) Plots of stand traces of free walking animals and (L) summary of paw positioning in Npr2^fl/fl (n = 4) and Npr2^fl/fl;Wnt1^Cre (n = 4) are shown. AEP, anterior extreme position; PEP, posterior extreme postion; L, left; R, right; F, forepaw; H, hind paw; error bar: SEM. No significant differences were observed between Npr2^fl/fl and Npr2^fl/fl;Wnt1^Cre mice in any of these tests.

In summary, no deficits were observed in these behavioral assays in conditional Npr2 mutants (Figure 3). Npr2-deficient mutants possess considerable coordination capabilities.

Heat Sensation Is Impaired in the Absence of Sensory Axon Bifurcation in Npr2^fl/fl;Wnt1^Cre Mice

To assess the consequences of impaired axonal bifurcation of sensory neurons for pain processing in vivo, we analyzed the nociceptive behavior of Npr2^fl/fl;Wnt1^Cre and respective control mice in various animal models of pain. We tested both male and female mice, but no significant effects of sex were detected in any assay.

We first determined the acute pain thresholds for noxious thermal stimuli in Npr2^fl/fl;Wnt1^Cre and littermate control mice using the hot plate test (50–54°C). Npr2^fl/fl;Wnt1^Cre mice displayed typical nocifensive behaviors (shaking or licking of the hindpaws or jumping off the plate), but at all tested temperatures they did so at consistently longer latencies than control mice (Figure 4A). Paw withdrawal responses to radiant heat (Hargreaves method) were also significantly increased in Npr2^fl/fl;Wnt1^Cre mice as compared to control mice (Figure 4B). In contrast to the impaired heat sensitivity, responses to acetone-evoked evaporative cooling were normal in Npr2^fl/fl;Wnt1^Cre mice in the group tested (Figure 4C; p = 0.295). Moreover, mechanical sensitivity of the hindpaw assessed using a dynamic plantar aesthesiometer was similar in Npr2^fl/fl;Wnt1^Cre and control mice (Figure 4D). Together, these data suggest that the loss of axonal bifurcation impairs the rapid response to avoid noxious heat, whereas behavioral thresholds and response latencies to cold or mechanical stimuli were not affected.

FIGURE 4

Figure 4. Pain behavior in Npr2^fl/fl;Wnt1^Cre and littermate control mice. (A) Thermal heat pain was assessed on a 50, 52, and 54°C hot plate (n = 15–16; *p < 0.05; power 97, 99, and 53% for 50, 52, and 54°C, respectively). (B) Thermosensation by the Hargreaves method (Plantar test) performed at 3 consecutive days (n = 8; *p < 0.05 and **p < 0.01; power 70, 56, and 85% for day 1, 2, and 3, respectively). (C) Cold pain was investigated using acetone-evoked evaporative cooling (n = 15–16, power 17%). (D) Mechanical sensitivity was measured using a dynamic plantar test (n = 11, power 11%). (E,F) The nocifensive response to chemical pain stimuli was determined by injection of 5 μg capsaicin (n = 9, power 36 and 55% for licking and paw shakes, respectively) or 0.5% formalin (n = 5, power 5 and 84% for phase 1 and 2, respectively) into a hindpaw. *p < 0.05.

We next examined the tonic pain behavior in Npr2^fl/fl;Wnt1^Cre and control mice after injection of capsaicin or formalin into a hindpaw. The capsaicin injection (5 μg) resulted in licking and shaking the injected paw in both genotypes. Whereas, the licking behavior was not statistically different between groups (p = 0.104), the number of paw shakes was significantly reduced in Npr2^fl/fl;Wnt1^Cre mice compared to control mice in the test group (Figure 4E). Injection of 0.5% formalin into a hindpaw resulted in the typical biphasic paw licking. In the first phase (1–10 min), which results from acute activation of nociceptors, the licking behavior of Npr2^fl/fl;Wnt1^Cre and control littermates was comparable (Figure 4F). However, the formalin-induced paw licking behavior was significantly reduced in Npr2^fl/fl;Wnt1^Cre mice in the second phase (11–60 min), which involves a combination of ongoing afferent firing and central sensitization (Shields et al., 2010). These data indicate that nociception induced by the chemical irritants capsaicin or formalin are impaired by the loss of axonal bifurcation. Consistent with these results are our previously published electrophysiological recordings on global Npr2 and CNP knockouts which revealed a reduced number of neurons in the dorsal spinal cord responding to capsaicin (Schmidt et al., 2007, 2009).

Npr2 Expression in DRGs Declines Postnatally

In addition to its role in sensory axon bifurcation Npr2 could potentially have an acute function in the processing of sensory information for example in noxious heat perception or nociception induced by chemical irritants. Furthermore, the second messenger cGMP and cGKI contribute to the processing of inflammatory and neuropathic pain in the spinal cord (Tegeder et al., 2004; Schmidtko et al., 2008, 2009; Luo et al., 2012). To further explore whether the measured deficits in pain behaviors in Npr2 mutants are caused by the absence of sensory axon bifurcation or by changes in primary acute pain perception, we studied the developmental expression profile of Npr2 in DRGs by histochemical determination of β-galactosidase activity in cryostat sections from the Npr2^wt/lacZ reporter mouse and by Western blot analysis of tissue extracts of DRGs (Ter-Avetisyan et al., 2014). Consistent with our previous in situ hybridization data we observed that all DRG neurons strongly express Npr2 at early embryonic stages when sensory axons enter the spinal cord (Figures 5A, C–E). At postnatal stages, however, the level and the number of Npr2 expressing cells gradually decreases in DRGs and at post-natal day 75 only 10% were positive for Npr2. Similar results were obtained by analyzing red fluorescent protein expression in DRGs at P75 which was induced by crossing the Rosa26-tdTomato reporter mouse line with Npr2-CreERT2 and tamoxifen application at 3 consecutive days (P71, P72, and P73) (data not shown). Further analysis revealed that Npr2 is not confined to a specific subset of DRG neurons. It is expressed in a small group of trkA-, parvalbumin-, CGRP-, or IB4-expressing DRG neurons. The percentage of Npr2-positive cells in these subpopulations - except for parvalbumin - decreased during maturation (Figures 5F–I).

FIGURE 5

Figure 5. Expression of Npr2 and cGKIα in DRGs at different developmental stages. (A,B) Western blot of membrane or cytosolic fractions of DRGs using antibodies to Npr2 or to cGKIα, respectively. The heavy chain of clathrin or GAPDH served as loading control. (C–E) Sections through thoracic DRGs at different stages. Localization of anti-βGal staining indicating Npr2-expression is demonstrated in the Npr2^wt/LacZ-mouse reporter and localization of cGKIα by a polyclonal antibody to mouse cGKIα. Scale bars in (C₁,D₁), 50 μm; in (C₂-C₅,D₂-D₅), 100 μm; in (E) 25 μm. (F) Quantification of Npr2-positive and cGKIα-positive neurons in DRGs counterstained with an antibody to PGP9.5 at different stages. (G–I) A small proportion of Npr2-positive DRG neurons at P45 and P75 express trkA, parvalbumin, CGRP or IB4. Scale bars in (H,I), 20 μm.

Cutaneous C-Fiber Thermal and Mechanical Sensitivity Is Normal in Npr2 Mutants

Co-stainings indicated that a small percentage of Npr2-positive DRG neurons at P75 express parvalbumin or trkA; therefore, it is possible that changes in peripheral nociceptor function may underlie altered heat sensitivity in Npr2^fl/fl;Wnt1^Cre mice. To test this, extracellular recordings from single C-fiber units in the hindpaw skin were performed using the ex vivo skin nerve preparation. C-fiber responses to heat, cold and mechanical stimulation of receptive fields were characterized in a total of 24 thermal-responsive C-fibers in Npr2^fl/fl;Wnt1^Cre mice, and 33 units in littermate control mice (Table S1). Example heat-evoked firing activity of single C-fibers in Npr2^fl/fl;Wnt1^Cre mice and littermate controls are shown in Figure 6A.

FIGURE 6

Figure 6. Responses of cutaneous hindpaw C-fibers to mechanical and thermal stimuli. Hindpaw cutaneous C-fibers with receptive fields on hairy foot skin were recorded from Npr2^fl/fl;Wnt1^Cre mice and littermate controls using the ex vivo skin nerve preparation. (A) Example traces of C-fiber firing activity recorded from Npr2^fl/fl;Wnt1^Cre mice and littermate controls in response to heat stimulation. (B) Heat thresholds of C-fibers were investigated using a 32–48°C 15 s heat ramp protocol. (C) The total number of spikes were quantified over the 4 s heat stimulus shown in (A). (D) Cold thresholds were investigated using a 32–12°C 15 s heat ramp protocol. (E) The total numbers of spikes were quantified over the 4 s heat stimulus shown in (A). (F) Mechanical thresholds of C-fibers. (G) Mechanical stimulus-evoked firing activity of C-fibers was quantified in response to four ramp and hold stimuli and compared between groups (see also Table S1).

C-fiber heat thresholds, defined as the temperature required to evoke the first action potential, did not differ between Npr2^fl/fl;Wnt1^Cre mice and littermate controls (Figure 6B). Similarly, there was no significant difference in the number of action potentials in C-fibers during heat stimulation between Npr2^fl/fl;Wnt1^Cre mice and littermate controls in response to suprathreshold heat stimulation (Figure 6C, exemplified in Figure 6A). Cold thresholds of C-fibers were also not different between Npr2^fl/fl;Wnt1^Cre mice and littermate controls (Figure 6D). Mechanical thresholds and suprathreshold mechanical stimulus-evoked firing activity of C-fibers also did not differ between Npr2^fl/fl;Wnt1^Cre mice and littermate controls (Figures 6E,F). Example mechanical-evoked firing activity of a C-fibers recorded from a Npr2^fl/fl;Wnt1^Cre mouse is shown in Figure 6G. These data indicate that thermal and mechanical coding properties of cutaneous C-fibers are not functionally impaired in Npr2 mutant mice.

Conditional Inactivation of cGKI Abolishes Sensory Axon Bifurcation and Partially Mimics the Phenotype of Npr2 Mutants in Pain-Related Behavior

cGKI is an essential signaling component downstream of Npr2 to induce sensory axon bifurcation as shown previously (Schmidt et al., 2002, 2007; Zhao et al., 2009). However, the developmental expression profile of cGKI is distinct from Npr2. At early embryonic stages when sensory axons are bifurcating in the spinal cord a complete overlap between Npr2 and cGKIα in all sensory neurons was observed (Figures 5B,D,E). In contrast to Npr2, cGKIα remains expressed in all sensory neurons also at mature stages suggesting that its function might be regulated by other cGMP synthesizing enzymes than by Npr2 in adult mice.

The absence of cGKI in a constitutive knockout does not interfere with bone growth but shows deficits in the regulation of smooth muscle contraction with vascular and intestinal dysfunctions which cause pre-mature death (Hofmann et al., 2006). Behavioral testing therefore also required conditional inactivation of cGKI in DRG neurons, which was obtained by a floxed allele of the cGKI gene (Wegener et al., 2002) and using Wnt1-Cre as described for Npr2. This resulted in a complete loss of cGKI protein in DRGs as indicated by Western blotting using antibodies to cGKI (Figure 7A). Consequently, sensory axons failed to form bifurcational branches at the DREZ (Figures 7B,C). In contrast to conditional Npr2 mutants, no decrease in the survival rate or any problems on the health status of 23 inspected cGKI^KO/fl;Wnt1^Cre animals in which cGKI was conditionally inactivated was observed. In some tissues Npr2-induced cGMP signals are mediated by the cGKII instead of cGKI which might explain these phenotypic differences on the survival rate between cGKI and Npr2.

FIGURE 7

Figure 7. Conditional inactivation of cGKI. (A) Conditional inactivation by Wnt1-Cre results in the complete absence of cGKI protein in embryonic DRGs as revealed by Western blot detection of cGKI in cytosolic fractions of E13.5 DRG from cGKI^wt/wt, cGKI^KO/fl, and cGKI^KO/fl;Wnt1^Cre mice. Equal loading was demonstrated by reprobing of the blot with an antibody to GAPDH. (B,C) Complete absence of sensory axon bifurcation in cGKI^KO/fl;Wnt1^Cre and quantification. The numbers of counted axons and embryos are indicated below the columns. Scale bar, 25 μm.

To estimate whether the observed Npr2-dependent changes in pain processing are mediated by a Npr2/cGMP/cGKI signaling pathway, we analyzed the nociceptive behavior of cGKI^KO/fl;Wnt1^Cre and littermate control mice (cGKI^KO/fl). In the hot plate test, cGKI^KO/fl;Wnt1^Cre mice demonstrated reduced noxious heat sensation (Figure 8A), albeit with significant differences compared to control mice only at 50 and 54°C. In the acetone test (Figure 8B) and dynamic plantar test (Figure 8C), the responses to cold and mechanical stimuli were normal in cGKI^KO/fl;Wnt1^Cre mice. After capsaicin injection into a hindpaw, the paw licking behavior of cGKI^KO/fl;Wnt1^Cre mice was reduced as compared to cGKI^KO/fl control mice (Figure 8D). Altogether, both Npr2^fl/fl;Wnt1^Cre and cGKI^KO/fl;Wnt1^Cre mice demonstrated impaired responses to noxious heat and chemical stimuli but normal responses to cold and mechanical stimuli, suggesting that the loss of bifurcation causes specific deficits in nociceptive processing.

FIGURE 8

Figure 8. Pain behavior in cGKI^KO/fl;Wnt1^Cre and littermate control mice. (A) Thermal heat pain was assessed on a 50, 52, and 54°C hot plate (n = 16; *p < 0.05; power 99, 47, and 90% for 50, 52, and 54°C, respectively). (B) Cold pain was investigated using acetone-evoked evaporative cooling (n = 16, power 11%). (C) Mechanical sensitivity was measured using a dynamic plantar test (n = 12, power 20%). (D) The nocifensive response to a chemical pain stimulus was determined by injection of 5 μg capsaicin (n = 9, power 61% for licking and 12% for paw shakes) into a hindpaw. *p < 0.05.

Discussion

Previous studies using knockout mice and axon tracing methods have implicated a cGMP signaling cascade composed of the ligand CNP, the receptor guanylyl cyclase Npr2 and the cGMP-dependent kinase I in bifurcation of DRG and cranial sensory ganglia (Schmidt and Rathjen, 2010; Gibson and Ma, 2011; Ter-Avetisyan et al., 2014). Since constitutive knockouts of Npr2 and cGKI are not suitable for further anatomical and behavioral studies on the consequences of the primary defect, the branching error at the DREZ, we generated mouse models that allowed an inactivation of Npr2 or cGKI specifically in embryonic DRGs at very early developmental stages before sensory axons enter the spinal cord. The Wnt1-Cre line was suitable to completely inactivate Npr2 or cGKI resulting in a complete lack of bifurcation for all DRG axons as previously described for the constitutive knockouts. Other Cre lines such as trkA-Cre or trkC-Cre did not induce bifurcation errors most likely by mediating recombination at relatively late stages of development and therefore cannot be used to assess a role of axonal bifurcation. The conditional inactivation of Npr2 by Wnt1-Cre allowed us to study the behavioral consequences and additional anatomical changes in the termination fields in the spinal cord caused by the disturbed axon bifurcation at the DREZ.

Due to the lack of axon bifurcation a reduction in the number of collaterals originating from primary afferents of DRG neurons within the spinal cord is not unexpected. However, the observed reduction to 64% of control values and not to 50% when total voxels numbers were summed up points to compensatory mechanisms such as an increase in the density of collaterals on stem axons, increased terminal branching of collaterals in spinal target areas or increased cell number in DRGs. The latter can be excluded since numbers were unchanged in the absence of cGMP signaling (Schmidt et al., 2002). Topographic representation of the sensory periphery within the spinal cord is of fundamental importance for information processing (Brown et al., 1991; Florence et al., 1991; Levinsson et al., 2002; Schouenborg, 2004; Granmo et al., 2008). Our work reported here showed that not only a quantitative reduction of incoming fibers but also a change in the pattern in the termination fields was detected in the absence of bifurcation. Npr2-deficiency caused an increase of the dorsoventral span to 131% of the termination field of digit two of the hindpaw whereas the mediolateral extension was narrowed by 25%. In principle pre- and postsynaptic mechanisms contribute to the formation of terminal fields. In the Npr2 mutants the balance between these interacting structures is disordered in the spinal cord. Our observations suggest that the incoming sensory axons and most likely less the postsynaptic neurons within the spinal cord might be the causative determinants on the development of these modified fields. Less sensory collaterals compete for postsynaptic neurons which might result in an increase of terminal branching in dorsoventral direction. In other words, in Npr2 mutants collaterals project not just to the normal zone of responding postsynaptic cells. The opposite is found for the mediolateral extension and no change was observed for the rostrocaudal span. The latter was expected since the longitudinal growth of the stem axons was not diminished in the absence of cGMP signaling (Schmidt et al., 2002).

In behavioral assays responses to noxious heat and chemical stimuli were impaired by the loss of axonal bifurcation either by the absence of Npr2 or cGKI which is in line with our previously published patch clamp recordings in Npr2- or CNP-deficient slices of the dorsal horn of the spinal cord. A reduced number of responding cells upon capsaicin treatment was measured in both mutants in comparison to wild type tissue, whereas synaptic transmission per se was not affected (Schmidt et al., 2007, 2009). In addition, in a recent study using global Npr2 mutants - in which, however, bone growth is disturbed - deficits in the auditory system have been described (Lu et al., 2014). Although the downstream effector cGKI shows a widespread pattern of localization in DRGs that is distinct from the more restricted expression of Npr2 in mature stages, similar - however not fully identical - deficits in nociception were observed in the behavioral tests upon the conditional ablation of cGKI. Its widespread localization in mature DRG neurons, its multiple phosphorylation targets (Hofmann and Wegener, 2013) and its involvement in synaptic potentiation in spinal neurons (Luo et al., 2012) might contribute to the differences in the licking time vs. paw shakes that were observed upon capsaicin treatment in Npr2 and cGKI conditional mutants. Our interpretation that the deficits in noxious heat perception in Npr2 and cGKI mutated mice were caused by the absence of axon bifurcation are in line with studies on a mouse mutant in which cGKI was inactivated by a Cre recombinase under control of the Na_v1.8 promoter (SNS-Cre). These mice do not show axon bifurcation defects at the DREZ due to the expression of the Cre recombinase after the ingrowth of sensory axons. They responded normally to noxious heat as measured by paw withdrawal to radiant heat (Hargreaves method) (Luo et al., 2012). In contrast to nociception, balance and motor performance appear normal in the conditional absence of Npr2-mediated bifurcation suggesting that pain processing might be more sensitive to a lack of bifurcation than proprioception. Moreover, nociceptors and other chemoreceptors have small-diameter axons which are either unmyelinated or thinly myelinated and therefore conduct action potentials more slowly. In contrast, proprioception, which is mediated by myelinated axons with a large diameter, might be more flexible due to fast conductance of action potentials and might therefore adjust more rapidly. It is also conceivable that elaborate compensatory mechanisms are implemented in the proprioceptive system to reorganize neuronal circuits in the absence of bifurcation. These compensational mechanisms might be less active in the nociceptive system.

It also cannot be excluded that the significant differences observed between nociception and proprioception in the absence of bifurcation might be related to the experimental systems that were applied. Pain was always analyzed by paw withdrawal reactions and a time delay in responses was measured (Weng and Schouenborg, 1998). In such a system pain perception needs to be transferred to muscle activity and this might require axon bifurcation at the DREZ to be well-coordinated and effectively executed. Further tactile feedback might be important for the nociceptive withdrawal reflex system. In addition accurate pain signaling might involve the integration of sensory information from several afferent subtypes and might require a complex crosstalk between spinal interneurons which might not be fully developed in the absence of axon bifurcation.

An acute or a morphological role of Npr2 on interneurons in the spinal cord or in the brain which might contribute to the measured behavioral deficits cannot currently be fully excluded. However, we never observed any changes in the arrangement of specific cell layers in the spinal cord, in the overall growth of proprioceptive collaterals or any synaptic deficits in the absence of Npr2 or CNP (Schmidt et al., 2007, 2009). Bifurcation errors of sensory axons were the only changes in the absence of the CNP/Npr2/cGKI signaling cascade we observed in the nervous system so far. Our electrophysiological recordings on the C-fibers in the nerve-skin preparation make it unlikely that Npr2 has an acute function in the conduction of sensory information. Nevertheless, a detailed pattern of expression of Npr2 in the developing and mature brain and spinal cord remains to be established including investigations on a co-localization with cGKI and Wnt1-Cre recombinase. In the future selective inactivation of Npr2 in spinal interneurons but not in DRG neurons might further substantiate or relativize our interpretations on the deficits caused by axonal bifurcation.

In humans, loss-of-function mutations in the NPR2 gene cause a skeletal dysplasia, termed acromesomelic dysplasia type Maroteaux (AMDM) with an extremely short and disproportionate stature (Bartels et al., 2004). AMDM is a rare autosomal-recessive genetic disorder with a prevalence of one out of 1,000,000. Point mutations leading to single amino acid exchanges were found throughout the NPR2 protein including the ligand binding, kinase homology or guanylyl cyclase domain (for a synopsis of human NPR2 mutations see Vasques et al., 2014). A number of these NPR2-missense mutations resulted in retention of the protein in the endoplasmic reticulum and poor targeting of the protein to the plasma membrane (Hume et al., 2009; Vasques et al., 2013; Amano et al., 2014; Wang et al., 2016) while others reached the cell surface (Dickey et al., 2016). Radiographic images demonstrated abnormal growth plates and short bones in the limbs detectable by two years of age. Carrier parents of individuals with AMDM and heterozygous mutations in NPR2 are associated with a slight reduction in height in comparison to the population average (Olney et al., 2006; Vasques et al., 2013; Hisado-Oliva et al., 2015) which was also found for mice heterozygous for Npr2 (Tamura et al., 2004). Heterozygosity in Npr2 or cGKI does not affect axonal bifurcation in mice (Schmidt et al., 2007). In contrast, NPR2 gain-of-function mutations result in an overgrowth syndrome (Miura et al., 2014). Currently, it can only be hypothesized whether the absence of the Npr2-mediated cGMP signaling axis in DRG neurons causes branching errors of sensory axons within the spinal cord in these patients. In addition, detailed neurological tests which might reveal neurological deficits have not been performed in AMDM patients so far. However, studies on mutations in rodents on axonal guidance factors have established a good correlation between animal models and human diseases (Jen et al., 2004; Engle, 2010; Srour et al., 2010; Depienne et al., 2011; Nugent et al., 2012; Chedotal, 2014). Our data using mouse genetics indicated deficits in noxious heat perception and nociception induced by the chemical irritants capsaicin or formalin. These observations might provide a framework for future studies to characterize neurological qualities of human patients with mutations in the NPR2 gene.

Author Contributions

PT, JH, JP, OD, FS, GT-A, SP, SF, FR, and HS performed experiments. Conception, drafting, approval, and agreement to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved: PT, JH, JP, OD, AS, FS, GL, GT-A, YW, SP, SF, RF, FR, and HS.

Funding

This work was supported by the MDC and DFG grants SFB 665 (grant B2 and B6 to FR and GL, respectively) and FOR 2060 (FE 438/5-1 and FE 438/6-1 to RF, and SCHM 2371/1 to HS, and SCHM 2629/3-1 to AS).

Conflict of Interest Statement

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.

Acknowledgments

The technical help of Madlen Driesner, Mechthild Henning, and Karola Bach is greatly acknowledged. We thank Melissa Long from the Animal Outcome Core Facility (Charité Universitätsmedizin Berlin) for her dedicated support. We thank Dr. Christina Eichhorn (MDC) for help on statistical evaluations of data, Dr. César S. Mendes (Universida de Nova, Lisboa) for discussions on the MouseWalker and Drs. Carmen Birchmeier (MDC, Berlin) and Marina Gültig (University of Tübingen) for mice.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnmol.2018.00019/full#supplementary-material

Figure S1. Severe malocclusions of P22 Npr2^fl/fl;Wnt1^Cre mice which prevents normal nutrition and correlated with premature death. Such animals have been excluded from behavioral testing.

Figure S2. No bifurcation errors of DRG axons were observed when the Npr2-flox line was crossbred with trkA-Cre or trkC-Cre. (A) Example of a DiI-labeled axon from E12.5 Npr2^fl/fl;trkA^Cre. Scale bar, 50 μm. (B) Quantification of bifurcation behavior of DRG axons in Npr2^fl/fl and Npr2^fl/fl;trkA^Cre embryos at E12.5. Numbers below the columns indicate counted axons and embryos. (C) Dorsolateral view of DRG axons at the dorsal root entry zone of a Npr2^fl/fl;trkC^Cre;Thy1^YFP−H mouse at P21. Bar, 100 μm.

Table S1. C-fiber ex vivo skin nerve preparation recorded from the hindpaw hairy skin. C-mechanoheat (CMH), C-mechanoheatcold (CMHC), and C-mechanocold (CMC) fibers were recorded from Npr2^fl/fl;Wnt1^Cre mice and littermate controls. Mean CVs are shown ± SEM.

References

Amano, N., Mukai, T., Ito, Y., Narumi, S., Tanaka, T., Yokoya, S., et al. (2014). Identification and functional characterization of two novel NPR2 mutations in Japanese patients with short stature. J. Clin. Endocrinol. Metab. 99, E713–E718. doi: 10.1210/jc.2013-3525

PubMed Abstract | CrossRef Full Text | Google Scholar

Aoyagi, Y., Kawakami, R., Osanai, H., Hibi, T., and Nemoto, T. (2015). A rapid optical clearing protocol using 2,2'-thiodiethanol for microscopic observation of fixed mouse brain. PLoS ONE 10:e0116280. doi: 10.1371/journal.pone.0116280

PubMed Abstract | CrossRef Full Text | Google Scholar

Bartels, C. F., Bükülmez, H., Padayatti, P., Rhee, D. K., van Ravenswaaij-Arts, C., Pauli, R. M., et al. (2004). Mutations in the transmembrane natriuretic peptide receptor NPR-B impair skeletal growth and cause acromesomelic dysplasia, type Maroteaux. Am. J. Hum. Genet. 75, 27–34. doi: 10.1086/422013

PubMed Abstract | CrossRef Full Text | Google Scholar

Braz, J. M., and Basbaum, A. I. (2010). Differential ATF3 expression in dorsal root ganglion neurons reveals the profile of primary afferents engaged by diverse noxious chemical stimuli. Pain 150, 290–301. doi: 10.1016/j.pain.2010.05.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Brown, A. G. (1981). Organization in the Spinal Cord. Berlin; Heidelberg; New York, NY: Springer.

Google Scholar

Brown, P. B., Gladfelter, W. E., Culberson, J. C., Covalt-Dunning, D., Sonty, R. V., Pubols, L. M., et al. (1991). Somatotopic organization of single primary afferent axon projections to cat spinal cord dorsal horn. J. Neurosci. 11, 298–309.

PubMed Abstract | Google Scholar

Carter, R. J., Morton, J., and Dunnett, S. B. (2001). Motor coordination and balance in rodents. Curr. Protoc. Neurosci. Chapter 8:Unit. doi: 10.1002/0471142301.ns0812s15

PubMed Abstract | CrossRef Full Text | Google Scholar

Caspani, O., Zurborg, S., Labuz, D., and Heppenstall, P. A. (2009). The contribution of TRPM8 and TRPA1 channels to cold allodynia and neuropathic pain. PLoS ONE 4:e7383. doi: 10.1371/journal.pone.0007383

PubMed Abstract | CrossRef Full Text | Google Scholar

Chedotal, A. (2014). Development and plasticity of commissural circuits: from locomotion to brain repair. Trends Neurosci. 37, 551–562. doi: 10.1016/j.tins.2014.08.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Chusho, H., Tamura, N., Ogawa, Y., Yasoda, A., Suda, M., Miyazawa, T., et al. (2001). Dwarfism and early death in mice lacking C-type natriuretic peptide. Proc. Natl. Acad. Sci. U.S.A. 98, 4016–4021. doi: 10.1073/pnas.071389098

PubMed Abstract | CrossRef Full Text | Google Scholar

Conte, W. L., Kamishina, H., and Reep, R. L. (2009). Multiple neuroanatomical tract-tracing using fluorescent Alexa Fluor conjugates of cholera toxin subunit B in rats. Nat. Protoc. 4, 1157–1166. doi: 10.1038/nprot.2009.93

PubMed Abstract | CrossRef Full Text | Google Scholar

Costantini, I., Ghobril, J. P., Di Giovanna, A. P., Allegra Mascaro, A. L., Silvestri, L., Mullenbroich, M. C., et al. (2015). A versatile clearing agent for multi-modal brain imaging. Sci. Rep. 5:9808. doi: 10.1038/srep09808

PubMed Abstract | CrossRef Full Text | Google Scholar

Danielian, P. S., Muccino, D., Rowitch, D. H., Michael, S. K., and McMahon, A. P. (1998). Modification of gene activity in mouse embryos in utero by a tamoxifen-inducible form of Cre recombinase. Curr. Biol. 8, 1323–1326. doi: 10.1016/S0960-9822(07)00562-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Depienne, C., Cincotta, M., Billot, S., Bouteiller, D., Groppa, S., Brochard, V., et al. (2011). A novel DCC mutation and genetic heterogeneity in congenital mirror movements. Neurology 76, 260–264. doi: 10.1212/WNL.0b013e318207b1e0

PubMed Abstract | CrossRef Full Text | Google Scholar

Dickey, D. M., Edmund, A. B., Otto, N. M., Chaffee, T. S., Robinson, J. W., and Potter, L. R. (2016). Catalytically active guanylyl cyclase b requires endoplasmic reticulum-mediated glycosylation, and mutations that inhibit this process cause dwarfism. J. Biol. Chem. 291, 11385–11393. doi: 10.1074/jbc.M115.704015

PubMed Abstract | CrossRef Full Text | Google Scholar

Echelard, Y., Vassileva, G., and McMahon, A. P. (1994). Cis-acting regulatory sequences governing Wnt-1 expression in the developing mouse CNS. Development 120, 2213–2224.

PubMed Abstract | Google Scholar

Engle, E. C. (2010). Human genetic disorders of axon guidance. Cold Spring Harb. Perspect. Biol. 2:a001784. doi: 10.1101/cshperspect.a001784

PubMed Abstract | CrossRef Full Text | Google Scholar

Farley, F. W., Soriano, P., Steffen, L. S., and Dymecki, S. M. (2000). Widespread recombinase expression using FLPeR (flipper) mice. Genesis 28, 106–110. doi: 10.1002/1526-968X(200011/12)28:3/4<106::AID-GENE30>3.0.CO;2-T

PubMed Abstract | CrossRef Full Text | Google Scholar

Feng, G., Mellor, R. H., Bernstein, M., Keller-Peck, C., Nguyen, Q. T., Wallace, M., et al. (2000). Imaging neuronal subsets in transgenic mice expressing multiple spectral variants of GFP. Neuron 28, 41–51. doi: 10.1016/S0896-6273(00)00084-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Florence, S. L., Wall, J. T., and Kaas, J. H. (1991). Central projections from the skin of the hand in squirrel monkeys. J. Comp. Neurol. 311, 563–578. doi: 10.1002/cne.903110410

PubMed Abstract | CrossRef Full Text | Google Scholar

Funfschilling, U., Ng, Y. G., Zang, K., Miyazaki, J., Reichardt, L. F., and Rice, F. L. (2004). TrkC kinase expression in distinct subsets of cutaneous trigeminal innervation and nonneuronal cells. J. Comp. Neurol. 480, 392–414. doi: 10.1002/cne.20359

PubMed Abstract | CrossRef Full Text | Google Scholar

Gibson, D. A., and Ma, L. (2011). Developmental regulation of axon branching in the vertebrate nervous system. Development 138, 183–195. doi: 10.1242/dev.046441

PubMed Abstract | CrossRef Full Text | Google Scholar

Granmo, M., Petersson, P., and Schouenborg, J. (2008). Action-based body maps in the spinal cord emerge from a transitory floating organization. J. Neurosci. 28, 5494–5503. doi: 10.1523/JNEUROSCI.0651-08.2008

PubMed Abstract | CrossRef Full Text | Google Scholar

Hargreaves, K., Dubner, R., Brown, F., Flores, C., and Joris, J. (1988). A new and sensitive method for measuring thermal nociception in cutaneous hyperalgesia. Pain 32, 77–88.

PubMed Abstract | Google Scholar

Hisado-Oliva, A., Garre-Vazquez, A. I., Santaolalla-Caballero, F., Belinchon, A., Barreda-Bonis, A. C., Vasques, G. A., et al. (2015). Heterozygous NPR2 mutations cause disproportionate short stature, similar to leri-weill dyschondrosteosis. J. Clin. Endocrinol. Metab. 100, E1133–E1142. doi: 10.1210/jc.2015-1612

PubMed Abstract | CrossRef Full Text | Google Scholar

Hofmann, F., and Wegener, J. W. (2013). cGMP-dependent protein kinases (cGK). Methods Mol. Biol. 1020, 17–50. doi: 10.1007/978-1-62703-459-3_2

PubMed Abstract | CrossRef Full Text | Google Scholar

Hofmann, F., Feil, R., Kleppisch, T., and Schlossmann, J. (2006). Function of cGMP-dependent protein kinases as revealed by gene deletion. Physiol. Rev. 86, 1–23. doi: 10.1152/physrev.00015.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Hume, A. N., Buttgereit, J., Al-Awadhi, A. M., Al-Suwaidi, S. S., John, A., Bader, M., et al. (2009). Defective cellular trafficking of missense NPR-B mutants is the major mechanism underlying acromesomelic dysplasia-type Maroteaux. Hum. Mol. Genet. 18, 267–277. doi: 10.1093/hmg/ddn354

PubMed Abstract | CrossRef Full Text | Google Scholar

Jen, J. C., Chan, W. M., Bosley, T. M., Wan, J., Carr, J. R., Rüb, U., et al. (2004). Mutations in a human ROBO gene disrupt hindbrain axon pathway crossing and morphogenesis. Science 304, 1509–1513. doi: 10.1126/science.1096437

PubMed Abstract | CrossRef Full Text | Google Scholar

Kloepper, J. E., Biro, T., Paus, R., and Cseresnyes, Z. (2010). Point scanning confocal microscopy facilitates 3D human hair follicle imaging in tissue sections. Exp. Dermatol. 19, 691–694. doi: 10.1111/j.1600-0625.2010.01110.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kobbert, C., Apps, R., Bechmann, I., Lanciego, J. L., Mey, J., and Thanos, S. (2000). Current concepts in neuroanatomical tracing. Prog. Neurobiol. 62, 327–351. doi: 10.1016/S0301-0082(00)00019-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuhn, M. (2016). Molecular physiology of membrane guanylyl cyclase receptors. Physiol. Rev. 96, 751–804. doi: 10.1152/physrev.00022.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

Levinsson, A., Holmberg, H., Broman, J., Zhang, M., and Schouenborg, J. (2002). Spinal sensorimotor transformation: relation between cutaneous somatotopy and a reflex network. J. Neurosci. 22, 8170–8182.

PubMed Abstract | Google Scholar

Lu, C. C., Cao, X. J., Wright, S., Ma, L., Oertel, D., and Goodrich, L. V. (2014). Mutation of Npr2 leads to blurred tonotopic organization of central auditory circuits in mice. PLoS Genet. 10:e1004823. doi: 10.1371/journal.pgen.1004823

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, R., Bausch, A. E., Kallenborn-Gerhardt, W., Stoetzer, C., Debruin, N., Ruth, P., et al. (2015). Slack channels expressed in sensory neurons control neuropathic pain in mice. J. Neurosci. 35, 1125–1135. doi: 10.1523/JNEUROSCI.2423-14.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, C., Gangadharan, V., Bali, K. K., Xie, R. G., Agarwal, N., Kurejova, M., et al. (2012). Presynaptically localized cyclic GMP-dependent protein kinase 1 is a key determinant of spinal synaptic potentiation and pain hypersensitivity. PLoS Biol. 10:e1001283. doi: 10.1371/journal.pbio.1001283

PubMed Abstract | CrossRef Full Text | Google Scholar

Luong, T. N., Carlisle, H. J., Southwell, A., and Patterson, P. H. (2011). Assessment of motor balance and coordination in mice using the balance beam. J. Vis. Exp. e2376 doi: 10.3791/2376

PubMed Abstract | CrossRef Full Text | Google Scholar

Madisen, L., Zwingman, T. A., Sunkin, S. M., Oh, S. W., Zariwala, H. A., Gu, H., et al. (2010). A robust and high-throughput Cre reporting and characterization system for the whole mouse brain. Nat. Neurosci. 13, 133–140. doi: 10.1038/nn.2467

PubMed Abstract | CrossRef Full Text | Google Scholar

Marques, J. M., and Olsson, I. A. (2010). Performance of juvenile mice in a reach-to-grasp task. J. Neurosci. Methods 193, 82–85. doi: 10.1016/j.jneumeth.2010.08.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Mendes, C. S., Bartos, I., Marka, Z., Akay, T., Marka, S., and Mann, R. S. (2015). Quantification of gait parameters in freely walking rodents. BMC Biol. 13:50. doi: 10.1186/s12915-015-0154-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Milenkovic, N., Zhao, W. J., Walcher, J., Albert, T., Siemens, J., Lewin, G. R., et al. (2014). A somatosensory circuit for cooling perception in mice. Nat. Neurosci. 17, 1560–1566. doi: 10.1038/nn.3828

PubMed Abstract | CrossRef Full Text | Google Scholar

Miura, K., Kim, O. H., Lee, H. R., Namba, N., Michigami, T., Yoo, W. J., et al. (2014). Overgrowth syndrome associated with a gain-of-function mutation of the natriuretic peptide receptor 2 (NPR2) gene. Am. J. Med. Genet. A 164A, 156–163. doi: 10.1002/ajmg.a.36218

PubMed Abstract | CrossRef Full Text | Google Scholar

Moshourab, R. A., Wetzel, C., Martinez-Salgado, C., and Lewin, G. R. (2013). Stomatin-domain protein interactions with acid-sensing ion channels modulate nociceptor mechanosensitivity. J. Physiol. 591, 5555–5574. doi: 10.1113/jphysiol.2013.261180

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakao, K., Osawa, K., Yasoda, A., Yamanaka, S., Fujii, T., Kondo, E., et al. (2015). The Local CNP/GC-B system in growth plate is responsible for physiological endochondral bone growth. Sci. Rep. 5:10554. doi: 10.1038/srep10554

PubMed Abstract | CrossRef Full Text | Google Scholar

Nugent, A. A., Kolpak, A. L., and Engle, E. C. (2012). Human disorders of axon guidance. Curr. Opin. Neurobiol. 22, 837–843. doi: 10.1016/j.conb.2012.02.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Ollion, J., Cochennec, J., Loll, F., Escude, C., and Boudier, T. (2013). TANGO: a generic tool for high-throughput 3D image analysis for studying nuclear organization. Bioinformatics 29, 1840–1841. doi: 10.1093/bioinformatics/btt276

PubMed Abstract | CrossRef Full Text | Google Scholar

Olney, R. C., Bukulmez, H., Bartels, C. F., Prickett, T. C., Espiner, E. A., Potter, L. R., et al. (2006). Heterozygous mutations in natriuretic peptide receptor-B (NPR2) are associated with short stature. J. Clin. Endocrinol. Metab. 91, 1229–1232. doi: 10.1210/jc.2005-1949

PubMed Abstract | CrossRef Full Text | Google Scholar

Potter, L. R. (2011). Guanylyl cyclase structure, function and regulation. Cell. Signal. 23, 1921–1926. doi: 10.1016/j.cellsig.2011.09.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Preibisch, S., Saalfeld, S., and Tomancak, P. (2009). Globally optimal stitching of tiled 3D microscopic image acquisitions. Bioinformatics 25, 1463–1465. doi: 10.1093/bioinformatics/btp184

PubMed Abstract | CrossRef Full Text | Google Scholar

Quina, L. A., Tempest, L., Ng, L., Harris, J. A., Ferguson, S., Jhou, T. C., et al. (2015). Efferent pathways of the mouse lateral habenula. J. Comp. Neurol. 523, 32–60. doi: 10.1002/cne.23662

PubMed Abstract | CrossRef Full Text | Google Scholar

Robertson, B., and Arvidsson, J. (1985). Transganglionic transport of wheat germ agglutinin-HRP and choleragenoid-HRP in rat trigeminal primary sensory neurons. Brain Res. 348, 44–51. doi: 10.1016/0006-8993(85)90357-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmidt, H., and Rathjen, F. G. (2010). Signalling mechanisms regulating axonal branching in vivo. Bioessays 32, 977–985. doi: 10.1002/bies.201000054

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmidt, H., and Rathjen, F. G. (2011). DiI-labeling of DRG neurons to study axonal branching in a whole mount preparation of mouse embryonic spinal cord. J. Vis. Exp. e3667. doi: 10.3791/3667

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmidt, H., Peters, S., Frank, K., Wen, L., Feil, R., and Rathjen, F. G. (2016). Dorsal root ganglion axon bifurcation tolerates increased cyclic GMP levels: the role of phosphodiesterase 2A and scavenger receptor Npr3. Eur. J. Neurosci. 44, 2991–3000. doi: 10.1111/ejn.13434

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmidt, H., Stonkute, A., Juttner, R., Koesling, D., Friebe, A., and Rathjen, F. G. (2009). C-type natriuretic peptide (CNP) is a bifurcation factor for sensory neurons. Proc. Natl. Acad. Sci. U.S.A. 106, 16847–16852. doi: 10.1073/pnas.0906571106

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmidt, H., Stonkute, A., Juttner, R., Schaffer, S., Buttgereit, J., Feil, R., et al. (2007). The receptor guanylyl cyclase Npr2 is essential for sensory axon bifurcation within the spinal cord. J. Cell Biol. 179, 331–340. doi: 10.1083/jcb.200707176

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmidt, H., Werner, M., Heppenstall, P. A., Henning, M., More, M. I., Kuhbandner, S., et al. (2002). cGMP-mediated signaling via cGKIalpha is required for the guidance and connectivity of sensory axons. J. Cell Biol. 159, 489–498. doi: 10.1083/jcb.200207058

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmidtko, A., Gao, W., Konig, P., Heine, S., Motterlini, R., Ruth, P., et al. (2008). cGMP produced by NO-sensitive guanylyl cyclase essentially contributes to inflammatory and neuropathic pain by using targets different from cGMP-dependent protein kinase I. J. Neurosci. 28, 8568–8576. doi: 10.1523/JNEUROSCI.2128-08.2008

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmidtko, A., Tegeder, I., and Geisslinger, G. (2009). No NO, no pain? The role of nitric oxide and cGMP in spinal pain processing. Trends Neurosci. 32, 339–346. doi: 10.1016/j.tins.2009.01.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Schneider, C. A., Rasband, W. S., and Eliceiri, K. W. (2012). NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671–675. doi: 10.1038/nmeth.2089

PubMed Abstract | CrossRef Full Text | Google Scholar

Schouenborg, J. (2004). Learning in sensorimotor circuits. Curr. Opin. Neurobiol. 14, 693–697. doi: 10.1016/j.conb.2004.10.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Shields, S. D., Cavanaugh, D. J., Lee, H., Anderson, D. J., and Basbaum, A. I. (2010). Pain behavior in the formalin test persists after ablation of the great majority of C-fiber nociceptors. Pain 151, 422–429. doi: 10.1016/j.pain.2010.08.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Sogawa, C., Abe, A., Tsuji, T., Koizumi, M., Saga, T., and Kunieda, T. (2010). Gastrointestinal tract disorder in natriuretic peptide receptor B gene mutant mice. Am. J. Pathol. 177, 822–828. doi: 10.2353/ajpath.2010.091278

PubMed Abstract | CrossRef Full Text | Google Scholar

Srour, M., Riviere, J. B., Pham, J. M., Dube, M. P., Girard, S., Morin, S., et al. (2010). Mutations in DCC cause congenital mirror movements. Science 328:592. doi: 10.1126/science.1186463

PubMed Abstract | CrossRef Full Text | Google Scholar

Staudt, T., Lang, M. C., Medda, R., Engelhardt, J., and Hell, S. W. (2007). 2,2′-thiodiethanol: a new water soluble mounting medium for high resolution optical microscopy. Microsc. Res. Technol. 70, 1–9. doi: 10.1002/jemt.20396

PubMed Abstract | CrossRef Full Text | Google Scholar

Tamura, N., Doolittle, L. K., Hammer, R. E., Shelton, J. M., Richardson, J. A., and Garbers, D. L. (2004). Critical roles of the guanylyl cyclase B receptor in endochondral ossification and development of female reproductive organs. Proc. Natl. Acad. Sci. U.S.A. 101, 17300–17305. doi: 10.1073/pnas.0407894101

PubMed Abstract | CrossRef Full Text | Google Scholar

Tegeder, I., Del, T. D., Schmidtko, A., Sausbier, M., Feil, R., Hofmann, F., et al. (2004). Reduced inflammatory hyperalgesia with preservation of acute thermal nociception in mice lacking cGMP-dependent protein kinase I. Proc. Natl. Acad. Sci. U.S A. 101, 3253–3257. doi: 10.1073/pnas.0304076101

PubMed Abstract | CrossRef Full Text | Google Scholar

Ter-Avetisyan, G., Rathjen, F. G., and Schmidt, H. (2014). Bifurcation of axons from cranial sensory neurons is disabled in the absence of Npr2-induced cGMP signaling. J. Neurosci. 34, 737–747. doi: 10.1523/JNEUROSCI.4183-13.2014

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsuji, T., and Kunieda, T. (2005). A loss-of-function mutation in natriuretic peptide receptor 2 (Npr2) gene is responsible for disproportionate dwarfism in cn/cn mouse. J. Biol. Chem. 280, 14288–14292. doi: 10.1074/jbc.C500024200

PubMed Abstract | CrossRef Full Text | Google Scholar

Vasques, G. A., Amano, N., Docko, A. J., Funari, M. F., Quedas, E. P., Nishi, M. Y., et al. (2013). Heterozygous mutations in natriuretic peptide receptor-B (NPR2) gene as a cause of short stature in patients initially classified as idiopathic short stature. J. Clin. Endocrinol. Metab. 98, E1636–E1644. doi: 10.1210/jc.2013-2142

PubMed Abstract | CrossRef Full Text | Google Scholar

Vasques, G. A., Arnhold, I. J., and Jorge, A. A. (2014). Role of the natriuretic peptide system in normal growth and growth disorders. Horm. Res. Paediatr. 82, 222–229. doi: 10.1159/000365049

PubMed Abstract | CrossRef Full Text | Google Scholar

Wan, X. C., Trojanowski, J. Q., and Gonatas, J. O. (1982). Cholera toxin and wheat germ agglutinin conjugates as neuroanatomical probes: their uptake and clearance, transganglionic and retrograde transport and sensitivity. Brain Res. 243, 215–224. doi: 10.1016/0006-8993(82)90244-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, W., Song, M. H., Miura, K., Fujiwara, M., Nawa, N., Ohata, Y., et al. (2016). Acromesomelic dysplasia, type maroteaux caused by novel loss-of-function mutations of the NPR2 gene: Three case reports. Am. J. Med. Genet. A 170A, 426–434. doi: 10.1002/ajmg.a.37463

PubMed Abstract | CrossRef Full Text | Google Scholar

Wegener, J. W., Nawrath, H., Wolfsgruber, W., Kuhbandner, S., Werner, C., Hofmann, F., et al. (2002). cGMP-dependent protein kinase I mediates the negative inotropic effect of cGMP in the murine myocardium. Circ. Res. 90, 18–20. doi: 10.1161/hh0102.103222

PubMed Abstract | CrossRef Full Text | Google Scholar

Weng, H. R., and Schouenborg, J. (1998). On the cutaneous receptors contributing to withdrawal reflex pathways in the decerebrate spinal rat. Exp. Brain Res. 118, 71–77. doi: 10.1007/s002210050256

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Z., and Ma, L. (2009). Regulation of axonal development by natriuretic peptide hormones. Proc. Natl. Acad. Sci. U.S.A. 106, 18016–18021. doi: 10.1073/pnas.0906880106

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Z., Wang, Z., Gu, Y., Feil, R., Hofmann, F., and Ma, L. (2009). Regulate axon branching by the cyclic GMP pathway via inhibition of glycogen synthase kinase 3 in dorsal root ganglion sensory neurons. J. Neurosci. 29, 1350–1360. doi: 10.1523/JNEUROSCI.3770-08.2009

PubMed Abstract | CrossRef Full Text | Google Scholar