NaV1.1 is essential for proprioceptive signaling and motor behaviors
Abstract
The voltage-gated sodium channel (NaV), NaV1.1, is well-studied in the central nervous system; conversely, its contribution to peripheral sensory neuron function is more enigmatic. Here, we identify a new role for NaV1.1 in mammalian proprioception. RNAscope analysis and in vitro patch-clamp recordings in genetically identified mouse proprioceptors show ubiquitous channel expression and significant contributions to intrinsic excitability. Notably, genetic deletion of NaV1.1 in sensory neurons caused profound and visible motor coordination deficits in conditional knockout mice of both sexes, similar to conditional Piezo2-knockout animals, suggesting that this channel is a major contributor to sensory proprioceptive transmission. Ex vivo muscle afferent recordings from conditional knockout mice found that loss of NaV1.1 leads to inconsistent and unreliable proprioceptor firing characterized by action potential failures during static muscle stretch; conversely, afferent responses to dynamic vibrations were unaffected. This suggests that while a combination of Piezo2 and other NaV isoforms is sufficient to elicit activity in response to transient stimuli, NaV1.1 is required for transmission of receptor potentials generated during sustained muscle stretch. Impressively, recordings from afferents of heterozygous conditional knockout animals were similarly impaired, and heterozygous conditional knockout mice also exhibited motor behavioral deficits. Thus, NaV1.1 haploinsufficiency in sensory neurons impairs both proprioceptor function and motor behaviors. Importantly, human patients harboring NaV1.1 loss-of-function mutations often present with motor delays and ataxia; therefore, our data suggest that sensory neuron dysfunction contributes to the clinical manifestations of neurological disorders in which NaV1.1 function is compromised. Collectively, we present the first evidence that NaV1.1 is essential for mammalian proprioceptive signaling and behaviors.
Editor's evaluation
This article provides insight into the importance of a voltage-gated sodium channel in proprioceptors, a group of mechanosensory neurons that target muscle. Using pharmacology, gene knockout, behavior, and histology in mice, the authors show quite convincingly that NaV1.1 in sensory neurons is essential for normal motor behavior and contributes to proprioceptor excitability. The work has interesting implications for human subjects with inherited variants of Nav1.1.
https://doi.org/10.7554/eLife.79917.sa0Introduction
Voltage-gated sodium channels (NaVs) are critical mediators of neuronal excitability and are responsible for action potential generation and propagation (Ahern et al., 2016; Bean, 2007). In the mammalian nervous system, there are nine isoforms (NaV1.1–1.9), each with unique biophysical properties, as well as distinguishing cellular expression and subcellular localization patterns (Bennett et al., 2019; Catterall, 2017). Of these different subtypes, NaV1.1 is notable for its role in brain disease (Escayg and Goldin, 2010; Mulley et al., 2005; Ogiwara et al., 2007). Indeed, Scn1a, the gene that encodes NaV1.1, is referred to as a ‘super culprit’ gene, with over 1000 associated mutations that lead to abnormal brain function, resulting in brain disorders such as epilepsy and migraine, as well as neurodivergent phenotypes, such as autism spectrum disorder (Ding et al., 2021; Lossin, 2009). Homozygous Scn1a-/- global knockout mice are ataxic and die by postnatal day (P) 15, while heterozygous Scn1a+/- animals develop seizures and begin to die sporadically starting at P21 (Yu et al., 2006). In addition to the central nervous system (CNS), NaV1.1 is also expressed in the peripheral nervous system (PNS) (Sharma et al., 2020; Usoskin et al., 2015); yet, the prominent role this channel plays in brain function has left its physiological roles in sensory neuron populations understudied.
Peripheral sensory neurons of the dorsal root ganglia (DRG) and trigeminal ganglia (TG) are tasked with encoding somatic sensations, such as touch, temperature, pain, and proprioception, and are anatomically and functionally heterogenous (Kupari et al., 2021; Nguyen et al., 2021; Oliver et al., 2021; Wu et al., 2021). Scn1a transcript and protein have been observed primarily in myelinated mechanosensory DRG and TG neurons (Fukuoka et al., 2008; Ho and O’Leary, 2011; Osteen et al., 2016). Indeed, subcutaneous injection of the NaV1.1 activator, Hm1a, into mouse hind paw causes noninflammatory mechanical pain and spontaneous pain behaviors (Osteen et al., 2016). Interestingly, pharmacological inhibition of NaV1.1 does not affect mechanical thresholds in uninjured mice but does reduce mechanical pain in a spared-nerve injury model (Salvatierra et al., 2018), suggesting that NaV1.1 may have a more prominent role in mechanical pain as opposed to normal touch sensing. NaV1.1 in TG neurons has also been reported to mediate mechanical pain in an orbitofacial chronic constriction injury model (Pineda-Farias et al., 2021). In addition to somatosensory neurons, NaV1.1 is found in colon-innervating vagal neurons, where it contributes to firing of colonic mechano-nociceptors and is upregulated in a mouse model of chronic visceral hypersensitivity (Osteen et al., 2016; Salvatierra et al., 2018). Lastly, NaV1.1 contributes to action potential firing in a subset of DRG neurons that express the cold-sensitive ion channel, transient receptor potential melastatin 8 (TRPM8), suggesting that the channel may also contribute to thermosensory transmission (Griffith et al., 2019). While most data support a role for NaV1.1 in pain, the limited number of studies that have investigated NaV1.1 function in sensory neurons has left gaps in our knowledge regarding other potential roles this channel may play in somatosensation.
Given the relatively underexplored role of NaV1.1 in the PNS, we set out to determine what other somatosensory modalities rely on NaV1.1 expression in sensory neurons. Here, we show that 100% of proprioceptors express Scn1a mRNA, where it makes notable contributions to the somal whole-cell sodium current and intrinsic excitability. A functional role for NaV1.1 in proprioceptive signaling was also supported by ex vivo electrophysiological recordings from functionally identified muscle spindle afferents. Importantly, mice lacking NaV1.1 in all sensory neurons display visible and profound motor deficits and ataxic-like behavior, which were quantified in rotarod and open-field assays. Surprisingly, we found that NaV1.1 is haploinsufficient for normal proprioceptor function and behavior, in ex vivo recordings and the open-field assay, respectively. Collectively, our data provide the first evidence that peripherally expressed NaV1.1 is critical for sensory proprioceptive signaling and motor coordination.
Results
Most studies have localized NaV1.1 expression primarily to myelinated sensory neurons that transmit mechanical signals (Fukuoka et al., 2008; Ho and O’Leary, 2011; Osteen et al., 2016; Wang et al., 2011). In line with prior work, RNAscope analysis of DRG sections from adult mice showed that 93% of myelinated neurons, as determined by neurofilament heavy chain (NFH) labeling, express Scn1a transcripts (Figure 1A). RNA-sequencing datasets have consistently identified NaV1.1 expression in proprioceptors; thus, we next analyzed NaV1.1 expression in these cells using a ParvalbuminCre;Rosa26Ai14 reporter line (PvalbAi14) and found that 100% of genetically identified proprioceptors were positive for Scn1a message (Figure 1B). This contrasted with low expression of Scn1a mRNA in both calcitonin gene-related peptide (CGRP)-expressing neurons, which represent peptidergic nociceptors, and non-peptidergic polymodal Mrgprd-expressing nociceptors (24 and 3%, respectively, Figure 1C and D). Frequency and cumulative distribution plots show the spread of integrated fluorescence density measurements obtained for Scn1a transcripts in proprioceptors (Figure 1E and F).
In addition to NaV1.1, NaV1.6 and NaV1.7 expression has also been reported in proprioceptors (Carrasco et al., 2017). As with Scn1a, mRNA for Scn8a and Scn9a is also found in 100% of genetically identified proprioceptors (Figure 2A and B). Cumulative distribution plots of Scn8a and Scn9a integrated fluorescence density measurements showed higher variability as compared to NaV1.1 (Figures 1E and 2C and D). This was quantified using the coefficient of variation (CV), a relative measure of the extent of variations within data. The CV for Scn1a transcript expression was calculated to be 75.6, whereas this value increased to 97.3 and 88.1 for Scn8a and Scn9a, respectively. This indicates that while all three isoforms are ubiquitously expressed in proprioceptors, the relative levels differ, with NaV1.1 having the most consistent level of expression across neurons analyzed. Furthermore, the average integrated density of the NaV1.1 signal for a given proprioceptive DRG neuron was significantly higher than the average integrated densities calculated for both NaV1.6 and NaV1.7 (Figure 2E, p<0.0001).
Due to the ubiquitous expression of NaV1.1, NaV1.6, and NaV1.7 in proprioceptors, we sought to determine the relative contributions of these isoforms to the proprioceptor whole-cell sodium current (INa). We performed in vitro voltage-clamp experiments on TdTomato+ neurons harvested from thoracic spinal levels of adult PvalbAi14 mice (de Nooij et al., 2013) and used serial application of selective NaV channel antagonists to determine the specific contributions of NaV1.1, NaV1.6, and NaV1.7 (Figure 2F). We first applied the selective NaV1.1 antagonist ICA 1214314 (ICA, 500 nM), followed by 9-anhydroustetrodoxin (AH-TTX, 300 nM), which is a selective NaV1.6 blocker but also partially blocks NaV1.1 (Denomme et al., 2020; Griffith et al., 2019). We reasoned that by first blocking the NaV1.1-mediated current, the effect of AH-TTX should largely be due to inhibition of NaV1.6. Finally, we blocked NaV1.7 channels using the antagonist PF-05089771 (25 nM), followed by tetrodotoxin (TTX, 300 nM) to block any residual current, as proprioceptors do not express TTX-resistant NaVs. On average, 7.8% of the current remained unblocked following serial application of these antagonists due to incomplete block by the drugs used. We found that the ICA-sensitive component was 44.8% of INa. Conversely, 14.5 and 32.9% of INa was sensitive to AH-TTX and PF-05089771, respectively (Figure 2G). No significant effect of the 0.1% DMSO vehicle solution on INa amplitude was observed (Figure 2—figure supplement 1). Collectively, these data suggest that in proprioceptors NaV1.1 is a dominant functional NaV subtype.
We next determined the biophysical features of the whole-cell sodium current (INa) in proprioceptors, which has not been previously reported (Figure 3A–D). The current–voltage relationship shows the first detectable current appeared at voltages near –50 mV and was maximal at voltages near –30 mV when evoked from a holding potential of –90 mV (Figure 3B). Voltage dependence of peak conductance was best fit to a single Boltzmann function and the voltage for half-maximal activation was –38.7 mV (Figure 3C). The voltage dependence of inactivation was determined with 40 ms prepulse steps ranging from –120 mV to +10 mV. The midpoint of the inactivation curve was –64.5 mV and was best fit to a single Boltzmann function (Figure 3D). To analyze recovery from fast inactivation, TdTomato+ neurons were depolarized to –20 mV, followed by a series of recovery periods ranging from 0.5 ms to 10 ms before a second test step to –20 mV was given to assess sodium channel availability. INa recovery was rapid (τ = 0.54ms), with greater than 50% of INa recovered after 0.5 ms (Figure 3E). Finally, entry into slow inactivation was determined. Cells were held at 0 mV during conditioning voltage steps ranging from 10 ms to 2000 ms, separated by two 2 ms pulses to –20 mV to compare channel availability before and after the conditioning pulse (Figure 3F). The tau for entry into slow inactivation was 928.6 ms, with more than 50% of channels available after a 2000 ms conditioning pulse.
We next asked how blocking NaV1.1 channels affects proprioceptor function in vitro. Similar to serial pharmacological experiments, INa density in proprioceptors was significantly reduced from ~–140 pA/pF to ~–75 pA/pF following ICA application (Figure 4B, p=0.0003). The proprioceptor INa had an average tau of 0.6 ms, which was significantly slowed to 1.0 ms following application of ICA (Figure 4C, p=0.0007), in line with loss of a fast gating channel. Blocking NaV1.1 did not change INa rise time (Figure 4D, p=0.1611). Quantification of the ICA-sensitive component found an average tau of 0.4 ms and an average current density of –96 pA/pF (Figure 4E). Of note, there was a wide distribution of current densities for the ICA-sensitive component, ranging from ~–28 pA/pF to ~–263 pA/pF, suggesting some variability in the contribution of NaV1.1 to proprioceptor excitability that may be proprioceptor subtype dependent. We next used current-clamp experiments to determine the effect of ICA on proprioceptor intrinsic excitability (Figure 4F). Pharmacological inhibition of NaV1.1 significantly reduced the number of evoked action potentials in most genetically identified proprioceptors (Figure 4G); however, five of the cells recorded had low firing rates that were not inhibited by ICA. This further suggests that NaV1.1 is important for repetitive firing in most proprioceptors, but some subtypes with lower intrinsic excitability instead rely on a combination of NaV1.6 and NaV1.7. Action potential amplitude (Figure 4H, p=0.0420) and action potential threshold (Figure 4I, p=0.0186) were also significantly reduced following ICA application, and action potential full-width half-max was significantly increased following ICA application (Figure 4K, p=0.0068), in line with loss of the fast NaV1.1-mediated current.
It is important to note that ICA 121431 also blocks NaV1.3 channels, which could be upregulated in our cultured DRG neuron preparations (Wangzhou et al., 2020). Indeed, a small but significant decrease in INa was observed in recordings from large-diameter DRG neurons harvested from PirtCre;Scn1a-floxed mice (Scn1a-cKO), which lack NaV1.1 in all sensory neurons. Thus, inhibition of upregulated NaV1.3 channels could contribute to the effect of on ICA on the proprioceptor INa(Figure 4—figure supplement 1).
To clarify the importance of NaV1.1 to proprioceptor function and avoid the caveats associated with in vitro pharmacological studies, we took an in vivo approach and analyzed motor behaviors in Scn1a-cKO mice of both sexes. We were precluded from using a Pvalbcre driver line to directly interrogate a role for NaV1.1 in proprioceptors as loss of NaV1.1 in Pvalb-expressing brain interneurons produces an epilepsy phenotype that prevents behavioral analyses in adult animals (Ogiwara et al., 2007). Consistent with in vitro data, Scn1a-cKO animals of both sexes displayed profound and visible motor abnormalities. These abnormalities include ataxic-like tremors when suspended in the air (Video 1), abnormal limb positioning (Videos 2 and 3), and paw clasping, which are absent in Scn1a-floxed littermate controls and heterozygous animals (PirtCre;Scn1afl/+, Scn1a-Het, respectively, Figure 5A). We first ran animals in the open-field test for 10 min each to quantify spontaneous locomotor behaviors (Figure 5B). We found that Scn1a-cKO animals traveled significantly less (Figure 5C) and slower (Figure 5D) than Scn1a-floxed littermate controls (p=0.0077 and p=0.0057, respectively). Surprisingly, Scn1a-Het mice also displayed motor abnormalities in the open-field test, performing similarly to Scn1a-cKO animals (Figure 5B–D), demonstrating NaV1.1 haploinsufficiency in sensory neurons for motor behaviors. No genotype-dependent differences were observed in the amount of time spent moving, suggesting gross motor function was intact (Figure 5E). Additionally, the amount of time spent in the center of the open-field chamber was also independent of genotype (Figure 5F). We next used the rotarod assay to investigate differences in motor coordination. Mice were assayed on three consecutive days and latency to fall and revolutions per minute (RPM) were quantified. Unlike in the open-field assay, both Scn1a-floxed and Scn1a-Het mice performed at similar levels during the 3-day period (Figure 5G and H). Conversely, Scn1a-cKO animals performed significantly worse. By day 3, on average they were only able to maintain their position on the rotarod for 41 s, falling over 50% faster Scn1a-floxed and Scn1a-Het mice. We did not observe any sex-dependent differences in performance in the open-field or rotarod tests (Figure 5—figure supplement 1). We confirmed that our mouse model selectively targeted sensory neurons by crossing a PirtCre driver with a fluorescent reporter line (PirtCre;Rosa26Ai14). We observed little-to-no neuronal expression of TdTomato in both dorsal and ventral spinal cord (Figure 5—figure supplement 2). In contrast, DRG somata and axons showed strong labeling. Collectively, our behavioral data provide evidence for a new in vivo role of NaV1.1 in sensory neurons in mammalian proprioception.
Could the motor deficits observed in Scn1a-cKO mice be due to abnormal proprioceptor development? To address this question, we performed RNAscope analysis of DRG sections from Scn1a-floxed, Scn1a-Het, and Scn1a-cKO mice. We quantified the number of neurons per DRG section that were positive for both Runx3 and Pvalb transcript, the molecular signature of mature proprioceptors (Figure 6, Oliver et al., 2021). We found no significant genotype-dependent differences in the number of proprioceptors in Scn1a-Het and Scn1a-cKO mice compared to Scn1a-floxed controls (p=0.3824 and p=0.1665, respectively), indicating that the behavioral deficits observed in Scn1a-cKO mice are not the result of a developmental loss of proprioceptors. We also analyzed muscle spindle morphology to determine whether aberrant sensory end organ development may contribute to the observed motor abnormalities. Similar to conditional Piezo2-knockout animals (Woo et al., 2015), no qualitative differences were observed between genotypes (Figure 6—figure supplement 1). Thus, abnormal proprioceptor development does not contribute to the overall phenotype of Scn1a-cKO mice.
We next asked whether the motor deficits of Scn1a-cKO mice are due to altered synaptic connectivity between proprioceptive axons and motor neurons in the ventral spinal cord. Spinal cord sections were harvested from Scn1a-floxed, Scn1a-Het, and Scn1a-cKO mice and stained with antibodies against vesicular glutamate transporter 1 (VGLUT1) to label proprioceptor axons, and choline acetyltransferase (ChAT). ChAT primarily labels α- and γ-motor neurons in the ventral horn, which can be distinguished based on size. We analyzed the number of VGLUT1 puncta on the somata and proximal dendrites of individual cholinergic neurons greater than 400 μm2 to bias our quantification towards α-motorneurons. We found no significant decrease in VGLUT puncta density per ChAT+ neuron between Scn1a-Het or Scn1a-cKO when compared to Scn1a-floxed littermate controls (Figure 7, Scn1a-Het, p>0.9999, Scn1a-cKO, p=0.4573). This suggests that general deficits in proprioceptor innervation of motor neurons do not contribute to the phenotype of Scn1a-cKO mice.
We next asked whether proprioceptor electrical signaling is altered in Scn1a-cKO mice. While in vitro patch-clamp electrophysiology can assess NaV function at DRG somata and provide insight as to how they contribute to intrinsic excitability, the physiological contributions of ion channels in DRG soma to somatosensory transmission in vivo are not well understood. Thus, to directly investigate how NaV1.1 shapes action potential propagation down proprioceptor axons, we used an ex vivo preparation to record muscle afferent activity during ramp-and-hold stretch and sinusoidal vibration. Afferents from both Scn1a-Het and Scn1a-cKO mice exhibited impaired static stretch sensitivity as evidenced by a decreased likelihood of firing during rest as compared to Scn1a-floxed mice, as well as an inability to maintain firing throughout the entire 4 s stretch. Almost all afferents from Scn1a-floxed mice could fire consistently throughout the entire 4 s hold phase (Figure 8A), but loss of one or both copies of NaV1.1 led to either firing only near the beginning of stretch or inconsistent firing in a high percentage of afferents lacking NaV1.1 (Figure 8B and C). We quantified this inconsistent firing by determining the CV of the interspike interval (ISI) during the plateau phase of stretch (1.5–3.5 s into the hold phase) across different stretch lengths and found a significant effect of genotype, with the knockout afferents both having higher ISI CV than the Scn1a-floxed afferents (Figure 8D; 0.074 ± 0.06, 0.313 ± 0.456, 0.497 ± 0.831, at 7.5% Lo, Scn1a-floxed, Scn1a-Het, and Scn1a-cKO afferents, respectively, two-way ANOVA, main effect of genotype, p=0.015). In contrast to the clear deficits in static sensitivity in afferents lacking NaV1.1, dynamic sensitivity was not significantly impaired. The maximum firing frequency during the rampup phase (Dynamic Peak) was independent of genotype, and even trended slightly higher in afferents lacking NaV1.1 (Figure 8E; two-way ANOVA, effect of genotype p=0.0633).
We next examined the requirement of NaV1.1 for proprioceptor afferent responses to sinusoidal vibration, which is a measure of dynamic sensitivity, and found no differences with loss of NaV1.1 (Figure 9A–C, Tables 1–3). We characterized a unit as having entrained to vibration if it fired at approximately the same time every cycle of the 9 s vibration. In most cases, afferents lacking NaV1.1 were equally likely to entrain to vibration than Scn1a-floxed afferents (Figure 9D–F). Indeed, Scn1a-cKO afferents were able to maintain firing during the entire 9 s sinusoidal vibration, in contrast to their inability to maintain consistent firing during 4 s of static stretch. There were no significant differences in firing rate during vibration between Scn1a-floxed, Scn1a-Het, and Scn1a-cKO afferents (Figure 9D–F). Taken together, our ex vivo recordings suggest that behavioral deficits in Scn1a-cKO result from abnormal proprioceptor responses to static muscle movement, whereas afferent responsiveness to dynamic stimuli is NaV1.1-independent. Furthermore, recordings from Scn1a-Het animals support the notion that NaV1.1 is haploinsufficient for proprioceptor function at the cellular level.
Discussion
The critical role for NaV1.1 in various brain disorders has overshadowed the potential contributions of this channel in peripheral signaling. The results presented in this study are the first to provide functional evidence that NaV1.1 in peripheral sensory neurons is required for normal proprioception. We found that mice lacking NaV1.1 in sensory neurons exhibit visible motor deficits and ataxic-like behaviors, which we propose is largely attributed to loss of NaV1.1 in proprioceptors. Indeed, RNAscope analysis showed expression of NaV1.1 mRNA in 100% of genetically identified proprioceptors. While NaV1.6 and NaV1.7 were also ubiquitously expressed in proprioceptors, NaV1.1 displayed higher expression levels, consistent with a previous RNA-sequencing study (Zheng et al., 2019). There are anywhere from 5 to 8 different proprioceptor molecular subclasses (Oliver et al., 2021; Wu et al., 2021), however, and it is possible that these distinct classes rely on different combinations of these channels for function. Nevertheless, our functional in vitro patch-clamp experiments found NaV1.1 to be the dominant subtype across recorded neurons, comprising nearly half of the proprioceptor INa. In line with this, pharmacological inhibition of NaV1.1 is sufficient to significantly attenuate action potential firing in most proprioceptors; yet, it should be noted that 25% of the neurons we recorded fired action potentials that were insensitive to ICA application. This suggests that some proprioceptor subtypes rely more heavily on NaV1.6 and NaV1.7 for electrical activity. Interestingly, the proprioceptors that were insensitive to ICA application also were less intrinsic excitable, firing only 1–2 action potentials in response to current injection, as opposed to the other 75% on neurons recorded, which fired repetitively during the 100 ms injection protocol. This further supports a role for NaV1.1 in maintaining action potential firing in response to sustained stimulation. These results, however, should be interpreted with the caveat that in these experiments ICA may also be acting on NaV1.3 channels that are upregulated during DRG neuron culturing (Wangzhou et al., 2020). Nevertheless, at the afferent level Scn1a-cKO and Scn1a-Het animals show clear deficits in static stretch sensitivity, but not dynamic sensitivity, and could even entrain to vibrations as fast as 100 Hz, suggesting a specific role for NaV1.1 in proprioceptors responses to static muscle movement. Finally, we found that loss of NaV1.1 in sensory neurons had no effect on proprioceptor development, muscle spindle morphology, or proprioceptive afferent innervation of ChAT+ motor neurons in the spinal cord. Thus, the observed motor behavioral deficits are likely due to reduced static sensitivity of proprioceptor afferents.
Our model proposes that NaV1.1 is tasked with maintaining consistent firing in spindle afferents during static muscle stretch for normal motor behaviors, whereby activation of the mechanotransduction channel Piezo2 initiates electrical signaling, which in turn activates a complement of TTX-sensitive NaV channels (Carrasco et al., 2017; Florez-Paz et al., 2016; Woo et al., 2015; Figure 10). During dynamic or vibratory stimuli, Piezo2, and likely a combination of other molecular mediators, including NaV1.6 and NaV1.7, are sufficient to elicit normal electrical activity. Conversely, during prolonged muscle stretch when Piezo2 channels presumably inactivate, NaV1.1 is required for regular and reliable firing. While other signaling molecules and channels, such as vesicle-released glutamate (Bewick et al., 2005; Than et al., 2021), and mechanosensitive ASIC channels (Lin et al., 2016) and ENaC channels (Bewick and Banks, 2015), also contribute to mammalian proprioception, our data suggests that NaV1.1 is a critical for muscle spindle afferent mechanotransduction, given the overt behavioral deficits observed in Scn1a-cKO mice. Due to the ubiquitous expression of NaV1.1 in all proprioceptors and the importance of Golgi tendon organ (GTO) feedback to motor control, alterations in function in those proprioceptors likely contribute to the behavioral deficits we observed; however, we did not directly measure their function.
To date, our knowledge of the functional contributions of NaV1.1 in the PNS is limited. Most studies have identified roles for this channel in mechanical pain signaling in DRG, TG, and vagal sensory neurons. Intraplantar pharmacological activation of NaV1.1 induces spontaneous pain behaviors and mechanical pain, which are absent in mice lacking NaV1.1 in small- and medium-diameter sensory neurons (Osteen et al., 2016). Inhibition of NaV1.1 prevented the development of mechanical pain in several preclinical models, including spared nerve injury (Salvatierra et al., 2018), an irritable bowel syndrome mouse model (Salvatierra et al., 2018), and infraorbital nerve chronic constriction injury (Pineda-Farias et al., 2021). Additionally, blocking NaV1.1 channels inhibited firing in TRPM8-expressing neurons in vitro, suggesting a potential role for this channel in thermosensation (Griffith et al., 2019). No prior studies, however, have reported a functional role for NaV1.1, or NaVs in general, in proprioception.
The loss of consistent firing we observed during static stretch in Scn1a-cKO and Scn1a-Het animals is functionally similar to deletion of NaV1.1 in other brain cell types. Indeed, loss of a single copy of NaV1.1 is sufficient to attenuate sustained action potential firing in parvalbumin-positive hippocampal interneurons (Ogiwara et al., 2007; Yu et al., 2006) and cerebellar Purkinje neurons (Yu et al., 2006). NaV1.1 has been associated with persistent sodium current (INaP) and resurgent sodium current (INaR), both of which promote repetitive firing in a wide variety of cell types in the CNS and PNS (Barbosa et al., 2015; Kalume et al., 2007; Khaliq et al., 2003). NaV1.1 also recovers rapidly from fast inactivation compared to other channel subtypes (Herzog et al., 2003; Patel et al., 2015) and has been shown to be refractory to entry into slow inactivation in TRPM8-expressing DRG neurons (Griffith et al., 2019). We observed similar characteristics when analyzing the proprioceptor INa. Future studies will determine whether proprioceptors rely on these features of NaV1.1 for reliable and consistent encoding of static muscle stretch.
Loss of NaV1.1 notably impacted proprioceptor afferent static sensitivity during ramp-and-hold stretch, but not dynamic sensitivity as measured by entrainment to sinusoidal vibrations using ex vivo muscle-nerve recordings. Afferents from Scn1a-cKO animals were more likely to have action potential failures and thus were largely unable to fire consistently throughout the 4 s of stretch, which was accompanied by a higher coefficient of variability in the ISI. This indicates that NaV1.1 has a critical role in transmitting static stretch information to the CNS. Interestingly, however, dynamic sensitivity in these afferents was unimpaired. Both Scn1a-cKO and Scn1a-Het afferents were able to entrain throughout the entire 9 s vibration; therefore, NaV1.1 does not appear to have a generalized role in maintaining high-frequency firing, but a more specific contribution to muscle-afferent static sensitivity. NaV1.1 has been localized to muscle spindle afferent endings and has been hypothesized to help amplify receptor current (Carrasco et al., 2017). Our results support a model whereby current from Piezo2 and potentially other mechanically sensitive ion channels at the start of stretch produces a sufficient receptor potential to generate firing at the heminode, but that amplification of the receptor potential by NaV1.1 is necessary to maintain firing during held stretch.
A similar deficit in static but not dynamic sensitivity was seen following loss of synaptic-like vesicle released glutamate from afferent endings (Than et al., 2021); however, in those afferents, firing only occurred at the beginning of stretch and patchy firing was never observed. This may indicate that glutamate plays a more general role maintaining excitability, whereas NaV1.1 is required for reliable action potential generation at heminodes during static stimuli. Alternatively, or in addition, NaV1.1 expressed along the axon could be essential for sustained static firing. A detailed examination of NaV1.1 subcellular localization along proprioceptor afferents could shed light on how this channel contributes to signal propagation. The lack of an effect on dynamic sensitivity could suggest the upregulation of other NaV subtypes or other molecules as a compensatory mechanism to counteract the loss of NaV1.1. Indeed, our in vitro electrophysiological experiments found a more pronounced effect of acute NaV1.1 inhibition on proprioceptor excitability. This could be due, however, to artificially upregulated NaV1.3 channel activity in culturing conditions, or conversely, a higher density of NaV1.1 expression in proprioceptor somata. Future studies using temporally controlled deletion of NaV1.1 in sensory neurons could tease this apart. Nevertheless, as static sensitivity is still very much impaired in both Scn1a-cKO and Scn1a-Het afferents, NaV1.1 may play a potentially unique role in maintaining afferent firing during the sustained stretch.
Loss of NaV1.1 in sensory neurons did not impact proprioceptor development as the number of proprioceptors in DRG sections was unchanged between genotypes, and muscle spindles developed normally in our model. While we did not directly examine GTOs, we do not anticipate that abnormal end organ development would be restricted to that proprioceptor subtype. Additionally, we did not observe a general decrease in α-motor neuron innervation, which is consistent with the findings of Mendelsohn et al., 2015, who reported loss of proprioceptor activity does not generally reduce proprioceptive input into the spinal cord. Interestingly, however, the authors did observe changes in heteronymous sensory-motor connectivity when proprioceptor transmission was blocked, whereby proprioceptor innervation of motor neurons that project to antagonistic muscles was increased. Whether loss of NaV1.1 in proprioceptive afferents causes them to ‘mis-wire’ in our model to make connections with inappropriate motor neuron pools is unclear. Future studies will determine whether this contributes to the observed behavioral deficits.
We found effects of both pharmacological inhibition and genetic deletion of NaV1.1 in in vitro and ex vivo electrophysiological experiments, respectively (Figures 2, 4, 8 and 9); however, in our mouse model NaV1.1 is deleted in all sensory neurons. Thus, we cannot rule out that loss of NaV1.1 in other mechanosenory neuron populations, such as touch receptors, contributes to the motor deficits observed. Deletion of NaV1.1 in small- and medium-diameter DRG neurons using a peripherin-Cre driver did not produce visible motor deficits (Osteen et al., 2016), indicating that sensory neuron populations in those categories are not involved. We did observe NaV1.1 transcripts in the vast majority of myelinated DRG neurons (a combination of large- and medium-diameter DRG neurons), consistent with its presence in different subclasses of tactile sensory neurons (Zheng et al., 2019). However, the severe motor phenotype of Scn1a-cKO mice precludes mechanical threshold analysis using von Frey or tactile sensitivity using tape test. Notably, baseline mechanical thresholds were unchanged following intraplantar injection of a selective NaV1.1 inhibitor (Salvatierra et al., 2018). This suggests that while NaV1.1 mRNA is expressed in most tactile sensory neurons, functional protein may only be upregulated in these populations during pathological states.
Despite this limitation, one noteworthy and intriguing finding from our study was the haploinsufficiency of NaV1.1 in sensory neurons for proprioceptor function and normal motor behavior in the open-field test. At the afferent level, heterozygous and homozygous loss of NaV1.1 produced similar deficits in static firing, suggesting that loss of less than a quarter of the proprioceptor INa is sufficient to impair proprioceptor responsiveness to muscle stretch. NaV1.1 is haploinsufficient in several brain neuron cell types for normal excitability and function (Ogiwara et al., 2007; Yu et al., 2006), suggesting the contributions of NaV1.1 to neuronal function are highly sensitive to genetic perturbations. At the behavioral level, Scn1a-Het mice had an identical phenotype to Scn1a-cKO mice, moving more slowly and less than controls, despite not having the more severe and visible motor coordination deficits. Indeed, their performance on the rotarod was identical to that of Scn1a-floxed controls (Figure 5). How these behavioral differences arise given the similar transmission deficits in Scn1a-cKO and Scn1a-Het afferents is unclear. One possibility is a presynaptic role for NaV1.1 in proprioceptive terminals that is unveiled when both copies of NaV1.1 are lost. For example, loss of presynaptic NaV1.7 channels in the spinal cord reduced glutamate release from nociceptive afferents onto dorsal horn neurons (MacDonald et al., 2021). If a similar mechanism is at play for NaV1.1 in proprioceptors, reduced neurotransmitter release from Scn1a-Het afferent terminals could be sufficient to produce quantifiable, albeit more subtle, motor deficits. Future studies are required to test this possibility.
Notably, Scn1a is a super culprit gene with over 1000 associated disease-causing mutations, most of which are linked to different forms of epilepsy. Many epilepsy patients with hemizygous NaV1.1 loss of function display ataxia and motor delays and deficiencies (Claes et al., 2001; Fujiwara et al., 2003), which has traditionally been attributed to loss of NaV1.1 function in the brain, namely, the cerebellum (Kalume et al., 2007). Our findings suggest that some of the clinical manifestations associated with epilepsy are not solely due to NaV1.1 loss of function in the brain, but also may manifest in part as a result of unreliable coding by peripheral proprioceptors.
Data presented in this study provide new evidence of a role for peripherally expressed NaV1.1 in motor coordination. We show that NaV1.1 is ubiquitously and strongly expressed in proprioceptors, contributes to proprioceptor excitability in vitro and ex vivo, and is haploinsufficient in sensory neurons for normal motor behaviors. Collectively, this work identifies a new role for NaV1.1 in mammalian proprioception.
Materials and methods
Key Resources Table contains a list of key resources and supplies used for this study.
Animals
Pirtcre mice were a kind gift from Dr. Xinzhong Dong (Johns Hopkins University). Rosa26Ai14 (stock #007914; Madisen et al., 2010) and Pvalbcre (stock # 008069) were obtained from Jackson Laboratories. Scn1a-floxed (stock #041829-UCD) mice were purchased from the UC Davis MMRRC. Mrgprd mice were originally published in Zheng et al., 2019. All mice used were on a mixed C57BL/6 background (non-congenic). Genotyping was outsourced to Transnetyx. Animal use was conducted according to guidelines from the National Institutes of Health’s Guide for the Care and Use of Laboratory Animals and was approved by the Institutional Animal Care and Use Committee of Rutgers University-Newark (PROTO201900161), UC Davis (#21947 and #22438), and San José State University (#990, ex vivo muscle recordings). Mice were maintained on a 12 hr light/dark cycle, and food and water was provided ad libitum.
Rotarod
Request a detailed protocolTo assess motor coordination, a rotarod machine (IITC Life Sciences, Woodland Hills, CA) that has an accelerating rotating cylinder was used. 8–10-week-old mice of both sexes were acclimated to the behavior room for 2 hr prior to testing. Mice were assayed on the rotarod for three consecutive days, with three trials per day and an intertrial interval of at least 15 min. The average of the three trials per day was used. The experimenter was blind to genotype.
Open-field test
Request a detailed protocol8–10-week-old mice of both sexes were acclimated to the behavior room for 2 hr prior to testing. The open-field apparatus consisted of a black square sound attenuating box of dimensions 40.6 cm × 40.6 cm. A camera suspended above the arena was connected to a computer running Ethovision XT software, which tracked animal movement and velocity. An animal was placed in the center of the arena and allowed to freely explore for a 10-min trial. The experimenter was blind to genotype.
Tissue processing
Request a detailed protocolFor spinal cord immunolabeling experiments, whole spinal columns from adult (8–15 weeks) PirtCre;Rosa26Ai14 and PirtCre;Scn1a-floxed animals of both sexes were harvested on ice. For Tdtomato immunohistochemistry, spinal columns were fixed overnight at 4°C in 4% paraformaldehyde. For vesicular glutamate transporter 1 (VGLUT1) and choline acetyltransferase (ChAT) co-labeling experiments, spinal columns were fixed in 4% paraformaldehyde for 1 hr on ice. Tissue was then placed in 30% sucrose solution overnight at 4°C. Following cryoprotection, tissue was embedded in optimal cutting temperature compound (OCT, Tissue-Tek Sakura) and stored at –80°C until sectioning. DRG from adult (8–15 weeks) animals of both sexes were harvested from thoracic spinal levels and fixed in 4% formaldehyde for 15 min at 4°C and were then incubated in 30% sucrose for 2–4 hr at 4°C. DRG were embedded in OCT and stored at –80°C until sectioning.
Immunohistochemistry
Request a detailed protocolImmunohistochemistry of spinal cord cryostat sections (30 μm) was performed using the following primary antibodies: rabbit anti-DsRed (1:3000, Takara Bio, 632496), guinea pig anti-VGLUT1 (1:8000, Zuckerman Institute, 1705), and rabbit anti-ChAT (1:10,000, Zuckerman Institute, 1574). Secondary antibodies used were as follows: anti-rabbit 594 (1:1000, Thermo Fisher, A32740), anti-guinea pig 488 (1:1000, Thermo Fisher, A11073), and anti-chicken 647 (Thermo Fisher, A32733). Specimens were mounted with Fluoromount-G with DAPI (SouthernBiotech, 0100-20). EDL muscles used in ex vivo muscle afferent recordings were placed in ice-cold 4% paraformaldehyde for 1 hr followed by ice-cold methanol for 15 min. Muscles were incubated in blocking solution (0.3% PBS-T and 1% BSA) followed by incubation in primary antibodies (guinea pig anti-VGLUT1 1:800 and chicken anti-NFH 1:300, Thermo Fisher ab4680) for 3–6 days at 4°C. After primary antibody treatment, tissue was washed in blocking solution and treated with secondary antibody (anti-guinea pig 488 1:50 and anti-chicken 594 1:300, Invitrogen, WA316328) for 2–3 days. Specimens were mounted with VECTASHIELD with DAPI (H-2000, Vector Laboratories). All specimens were imaged in three dimensions on a Zeiss LSM880 Airyscan confocal microscope. Images were analyzed using ImageJ software.
Multiplex in situ hybridization
Request a detailed protocolFixed-frozen DRG tissue from 8- to 15-week-old mice of both sexes was cut into 25 μm sections and placed on electrostatically coated slides. Sections were processed for RNA in situ detection using a modified version of the manufacturer’s protocol (Griffith et al., 2019, Advanced Cell Diagnostics) and the following probes: Pvalb (421931C1, mouse), Runx3 (451271-C3, mouse), Scn1a (556181-C2, mouse), Scn8a (313341-C2, mouse), and Scn9a (434191-C2, mouse). Following in situ hybridization, sections were incubated in blocking solution (5% normal goat serum, 0.1% PBS-T) for 1 hr at room temperature (RT). Tissue was then incubated in primary antibodies overnight at 4°C. The following antibodies were used: rabbit DsRed (1:3000, Takara Bio, 632496), rabbit β3-tubulin (1:3000, Abcam, ab18207), chicken β3-tubulin (1:500, Abcam, ab41489), rabbit CGRP (1:1000, ImmunoStar, 24112), chicken GFP (1:3000, Abcam, ab13970), and chicken NFH (1:3000, Abcam, ab4680). Tissue was treated with the following secondary antibodies for 45 min at RT: anti-rabbit 448 (1:1000, Invitrogen, A32731), 594 (1:1000, Invitrogen, A11037) and 647 (1:1000, Invitrogen, A32733), anti-chicken 488 (1:1000, Invitrogen, A32931) and 594 (1:1000, Invitrogen, A32740). Sections were washed and mounted with Fluoromount-G with DAPI and imaged in three dimensions (2 μm axial steps) on an Olympus confocal (LV3000) using a ×40 0.90 NA water objective lens. Images were autothresholded and the probe signal integrated density for individual neurons was analyzed using ImageJ software. The coefficients of variation for Scn1a, Scn8a, and Scn9a integrated densities were calculated in Prism 9.0 (GraphPad Software) using the following formula: CV = μ/σ * 100.
DRG culture preparation
Request a detailed protocolDRGs were harvested from thoracic spinal levels of adult Pvalbcre;Rosa26Ai14 and Pirtcre;Scn1a-floxed (6–16 weeks) mice of both sexes and transferred to Ca2+-free and Mg2+-free HBSS solution (Invitrogen, 14170-112). Upon isolation, processes were trimmed, and ganglia were transferred into collagenase (1.5 mg/mL; Type P, Sigma-Aldrich, 11213865001) in HBSS for 20 min at 37°C followed by TrypLE Express (Thermo Fisher, 12605-010) for 3 min with gentle rotation. TrypLE was neutralized with 10% horse serum (heat-inactivated; Invitrogen, 26050-070) and supplemented with culture media (MEM with l-glutamine, Phenol Red, without sodium pyruvate, Thermo Fisher, 11095-080), containing 10,000 U/mL penicillin-streptomycin (Thermo Fisher, 15140-122), MEM Vitamin Solution (Invitrogen, 11120-052), and B-27 supplement (Thermo Fisher, 17504-044). Serum-containing media was decanted and cells were triturated using a fire-polished Pasteur pipette in the MEM culture media described above. Cells were resuspended and triturated using a plastic pipette tip. Cells were plated on glass coverslips that had been washed in 2 M NaOH for at least 4 hr, rinsed with 70% ethanol, UV-sterilized, and treated with laminin (0.05 mg/mL, Sigma-Aldrich, L2020-1MG) for 1 hr prior to plating. Cells were then incubated at 37°C in 5% CO2. Cells were used for electrophysiology experiments 14–36 hr post-plating.
In vitro electrophysiology
Request a detailed protocolWhole-cell voltage-clamp recordings were made from dissociated DRG neurons using patch pipettes pulled from Model P-1000 (Sutter Instruments). Patch pipettes had a resistance of 3–5 MΩ when filled with an internal solution containing the following (in mM): 140 CsF, 10 NaCl, 1.1 EGTA,.1 CaCl2, 10 HEPES, and 2.5 MgATP, pH with CsOH to 7.2. Seals and whole-cell configuration were obtained in an external solution containing the following (in mM): 145 NaCl, 5 KCl, 10 HEPES, 10 glucose, 2 CaCl2, 2 MgCl2, pH 7.3 with NaOH, osmolarity ~320 mOsm. Series resistance ranged from 6 to 11 MΩ and was compensated by 70–80%. Voltage errors were not directly assessed in current-clamp recordings. To isolate whole-cell sodium currents during voltage-clamp experiments, a modified external solution was applied containing the following (in mM): 15 NaCl, 130 TEA-Cl, 10 HEPES, 2 BaCl2, 13 glucose, 0.03 CdCl2, pH 7.3 with NaOH, osmolarity ~320 mOsm. Voltage-clamp recordings were performed at RT and current-clamp recordings were conducted at 37°C. Bath temperature was controlled and monitored using CL-100 (Warner Instruments).
Ex vivo electrophysiology
Request a detailed protocolThe effect of the loss of NaV1.1 on muscle spindle afferent firing rates during muscle stretch and sinusoidal vibration was determined using an isolated muscle nerve preparation. The extensor digitorum longus muscle and innervating peroneal branch of the sciatic nerve were dissected from adult (2–4-month-old) mice of both sexes. Muscles were held at optimal length (Lo), or the length of the muscle that maximal force of twitch contraction occurred. A series of nine 4 s ramp-and-hold stretches were given to three different stretch lengths repeated three times each (2.5, 5, and 7.5% Lo; ramp speed 40% Lo/s). A series of twelve 9 s sinusoidal vibrations were given (25, 50, and 100 µm amplitude; 10, 25, 50, and 100 Hz frequency). A 1-min rest was given between each length change. Firing rates during a 10 s baseline before stretch (resting discharge or RD) and the maximal firing rate during the rampup phase of stretch (dynamic peak or DP) were calculated for all animals. We determined whether the response to static stretch was maintained consistently throughout the 4 s stretch, as well as the coefficient of variability of the ISI during the plateau phase of stretch (CV = Std Dev/Mean of ISI over the time period of 1.5–3.5 s after end of rampup). Average firing rate during the 9 s of vibration and whether the unit could entrain in a 1:1 fashion to vibration was also determined. Detailed methods can be found in Wilkinson et al., 2012.
Pharmacology
Request a detailed protocolICA 121431 (#5066), 4,9-anhydrotetrodotoxin (AH-TTX, #6159), and PF-05089771 (#5931) were purchased from Tocris Bioscience. Tetrodotoxin (TTX, ab120054) was purchased from Abcam. All other chemicals were from Sigma-Aldrich and Fisher Chemical.
Data acquisition and analysis
Request a detailed protocolCurrents and voltages were acquired using pClamp software v11.2 (Molecular Devices). Recordings were obtained using an AxoPatch 200b patch-clamp amplifier and a Digidata 1550B and filtered at 5 kHz and digitized at 10 kHz. For biophysical analysis of whole-cell sodium currents, conductance (G) was calculated as G = I/(V – ENa), in which I is the peak current, V is the voltage step, and ENa is the reversal potential for sodium calculated from Nernst equation based on the intracellular and extracellular sodium concentrations in our recording solutions (10.38 mV). Conductance data were normalized by the maximum conductance value, Gmax, and data was fit with the Boltzmann equation: Fraction available = Minimum + ([Maximum-Minimum]/[1 + exp(V50 Vm)/k]), where V50 denotes the membrane potential at which half the channels are inactivated and k denotes the Boltzmann constant/slope factor. For voltage dependence of steady-state inactivation, peak current data were normalized based on the maximum current, Imax. Analysis of action potential amplitude, full-width half max, and threshold was performed on the first action potential elicited in response to a 100 pA current injection (100 ms). Action potential threshold was calculated as the membrane potential at which the first derivative of the somatic membrane potential (dV/dT) reached 10 mV ms–1 (Griffith et al., 2019; Kress et al., 2008). Tau values were calculated from 20 ms voltage steps from –90 mV to –30 mV and analyzed with single-exponential curve fits. Voltage-clamp and current-clamp experiments were analyzed with Clampfit software v11.2 (Molecular Devices) and custom MATLAB Scripts. Ex vivo recordings were obtained using an A-M Systems Model 1800 extracellular amplifier with headstage and digitized using an ADInstruments PowerLab. Data was analyzed using ADInstruments LabChart software using the Spike Histogram function.
Experimental design and statistical analysis
Request a detailed protocolSummary data are presented as mean ± SEM, from n cells, afferents, or sections, or N animals. For quantitative analysis of in situ hybridization data, at least three biological replicates per condition were used and the investigator was blinded to genotype for analysis. Behavioral experiments and analysis were also performed genotype-blind. Statistical differences were determined using parametric tests for normally distributed data and nonparametric tests for data that did not conform to Gaussian distributions or had different variances. Statistical tests are listed in ‘Results’ and/or figure legends. Statistical significance in each case is denoted as follows: *p<0.05, **p<0.01, ***p<0.001, and ****p<0.0001. Statistical tests and curve fits were performed using Prism 9.0 (GraphPad Software). All data generated or analyzed during this study are included in the article. Source data files have been uploaded to Mendeley for all figures (https://data.mendeley.com/datasets/kt23th75v9). Code has been uploaded to GitHub (https://github.com/darikoneil/PropAnalysisScripts; O’Neil, 2022a). A Key Resources Table with specific organism and reagent information has been included in the ‘Materials and methods’ section.
Data availability
Source Data files have been uploaded to Mendeley for all figures (https://data.mendeley.com/datasets/kt23th75v9). Code has been uploaded to GitHub (https://github.com/doctheagrif/Current-Clamp-Matlab-Code_O-Neil-DA, copy archived at swh:1:rev:12667fd74eb68a8cd7ec858f33d9fec6e5c46772). A key resources table with specific organism and reagent information has been included in the methods section.
-
Mendeley DataNav1.1 is essential for proprioceptive signaling and motor behaviors.https://doi.org/10.17632/kt23th75v9.3
References
-
The hitchhiker’s guide to the voltage-gated sodium channel galaxyThe Journal of General Physiology 147:1–24.https://doi.org/10.1085/jgp.201511492
-
The action potential in mammalian central neuronsNature Reviews. Neuroscience 8:451–465.https://doi.org/10.1038/nrn2148
-
The role of voltage-gated sodium channels in pain signalingPhysiological Reviews 99:1079–1151.https://doi.org/10.1152/physrev.00052.2017
-
Mechanotransduction in the muscle spindlePflugers Archiv 467:175–190.https://doi.org/10.1007/s00424-014-1536-9
-
Distribution of TTX-sensitive voltage-gated sodium channels in primary sensory endings of mammalian muscle spindlesJournal of Neurophysiology 117:1690–1701.https://doi.org/10.1152/jn.00889.2016
-
Forty years of sodium channels: structure, function, pharmacology, and epilepsyNeurochemical Research 42:2495–2504.https://doi.org/10.1007/s11064-017-2314-9
-
De novo mutations in the sodium-channel gene SCN1A cause severe myoclonic epilepsy of infancyAmerican Journal of Human Genetics 68:1327–1332.https://doi.org/10.1086/320609
-
Comparative study of the distribution of the alpha-subunits of voltage-gated sodium channels in normal and axotomized rat dorsal root ganglion neuronsThe Journal of Comparative Neurology 510:188–206.https://doi.org/10.1002/cne.21786
-
Single-Cell analysis of sodium channel expression in dorsal root ganglion neuronsMolecular and Cellular Neurosciences 46:159–166.https://doi.org/10.1016/j.mcn.2010.08.017
-
A catalog of SCN1A variantsBrain & Development 31:114–130.https://doi.org/10.1016/j.braindev.2008.07.011
-
Molecular correlates of muscle spindle and golgi tendon organ afferentsNature Communications 12:1451.https://doi.org/10.1038/s41467-021-21880-3
-
SoftwareCurrent-clamp-matlab-code_o-neil-DA, version swh:1:rev:12667fd74eb68a8cd7ec858f33d9fec6e5c46772Software Heritage.
-
NIH image to imagej: 25 years of image analysisNature Methods 9:671–675.https://doi.org/10.1038/nmeth.2089
-
Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencingNature Neuroscience 18:145–153.https://doi.org/10.1038/nn.3881
-
Piezo2 is the principal mechanotransduction channel for proprioceptionNature Neuroscience 18:1756–1762.https://doi.org/10.1038/nn.4162
Article and author information
Author details
Funding
Burroughs Wellcome Fund (#1017423)
- Theanne N Griffith
National Institute of General Medical Sciences (5T32GM099608-10)
- Cyrrus M Espino
National Institute of General Medical Sciences (1T32GM1144303-01A1)
- Cyrrus M Espino
National Institute of General Medical Sciences (5SC3GM127195)
- Katherine A Wilkinson
National Institute of General Medical Sciences (5R25GM71381)
- Serena Ortiz
National Institute of Neurological Disorders and Stroke (R25NS063307)
- Kaylee M Wells
- Darik A O'Neil
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This research was supported by a Postdoctoral Enrichment Program Award from the Burroughs Welcome Fund (TNG) and 5T32GM099608-10 and 1T32GM1144303-01A1 (CME). Work at San José State was supported by NIGMS 5SC3GM127195 (KAW) and NIGMS 5R25GM71381 (SO). Core facilities were supported by P30 EY12576. Part of this project was carried out at the Marine Biological Laboratory Neurobiology Course, with support from NINDS R25NS063307. Thanks to Drs. Jon Sack and Xinzhong Dong for sharing mouse lines, Drs. Jorge Contreras and Ioana Carcea for sharing behavioral equipment, Miguel Gonzalez Fernandez for writing code to automate ex vivo data analysis, and members of the Griffith laboratory for helpful discussions.
Ethics
Animal use was conducted according to guidelines from the National Institutes of Health's Guide for the Care and Use of Laboratory Animals and was approved by the Institutional Animal Care and Use Committee of Rutgers University-Newark (PROTO201900161), UC Davis (#21947 and #22438) and San José State University (#990, ex vivo muscle recordings).
Copyright
© 2022, Espino et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,836
- views
-
- 377
- downloads
-
- 18
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
- Neuroscience
The blood-brain barrier (BBB) controls the movement of molecules into and out of the central nervous system (CNS). Since a functional BBB forms by mouse embryonic day E15.5, we reasoned that gene cohorts expressed in CNS endothelial cells (EC) at E13.5 contribute to BBB formation. In contrast, adult gene signatures reflect BBB maintenance mechanisms. Supporting this hypothesis, transcriptomic analysis revealed distinct cohorts of EC genes involved in BBB formation and maintenance. Here, we demonstrate that epigenetic regulator’s histone deacetylase 2 (HDAC2) and polycomb repressive complex 2 (PRC2) control EC gene expression for BBB development and prevent Wnt/β-catenin (Wnt) target genes from being expressed in adult CNS ECs. Low Wnt activity during development modifies BBB genes epigenetically for the formation of functional BBB. As a Class-I HDAC inhibitor induces adult CNS ECs to regain Wnt activity and BBB genetic signatures that support BBB formation, our results inform strategies to promote BBB repair.
-
- Neuroscience
Specialized chemosensory signals elicit innate social behaviors in individuals of several vertebrate species, a process that is mediated via the accessory olfactory system (AOS). The AOS comprising the peripheral sensory vomeronasal organ has evolved elaborate molecular and cellular mechanisms to detect chemo signals. To gain insight into the cell types, developmental gene expression patterns, and functional differences amongst neurons, we performed single-cell transcriptomics of the mouse vomeronasal sensory epithelium. Our analysis reveals diverse cell types with gene expression patterns specific to each, which we made available as a searchable web resource accessed from https://www.scvnoexplorer.com. Pseudo-time developmental analysis indicates that neurons originating from common progenitors diverge in their gene expression during maturation with transient and persistent transcription factor expression at critical branch points. Comparative analysis across two of the major neuronal subtypes that express divergent GPCR families and the G-protein subunits Gnai2 or Gnao1, reveals significantly higher expression of endoplasmic reticulum (ER) associated genes within Gnao1 neurons. In addition, differences in ER content and prevalence of cubic membrane ER ultrastructure revealed by electron microscopy, indicate fundamental differences in ER function.