# The Human SCN10AG1662S Point Mutation Established in Mice Impacts on Mechanical, Heat, and Cool Sensitivity

Céleste Chidiac, Yaping Xue, Maria del Mar Muniz Moreno, Ameer Abu Bakr Rasheed, Romain Lorentz, Marie-Christine Birling, Claire Gaveriaux-Ruff, Yann Herault

► To cite this version:

Céleste Chidiac, Yaping Xue, Maria del Mar Muniz Moreno, Ameer Abu Bakr Rasheed, Romain Lorentz, et al.. The Human SCN10AG1662S Point Mutation Established in Mice Impacts on Mechanical, Heat, and Cool Sensitivity. Frontiers in Pharmacology, 2021, 12, pp.780132. 10.3389/fphar.2021.780132. [hal-03706653](https://hal.science/hal-03706653)

HAL Id: hal-03706653

https://hal.science/hal-03706653

Submitted on 9 Dec 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# The Human SCN10A<sup>G1662S</sup> Point Mutation Established in Mice Impacts on Mechanical, Heat, and Cool Sensitivity

Celeste Chidiac<fnref n="1" />, Yaping Xue<fnref n="1" />, Maria del Mar Muniz Moreno<fnref n="1" />, Ameer Abu Bakr Rasheed<fnref n="1" />, Romain Lorentz<fnref n="2" />, Marie-Christine Birling<fnref n="2" />, Claire Gaveriaux-Ruff<fnref n="1" />† and Yann Herault<fnref n="1" />,<fnref n="2" />*†

The voltage-gated sodium channel NAV1.8 is expressed in primary nociceptive neurons and is involved in pain transmission. Mutations in the SCN10A gene (encoding NAV1.8 channel) have been identified in patients with idiopathic painful small fiber neuropathy (SFN) including the SCN10A<sup>G1662S</sup> gain-of-function mutation. However, the role of this mutation in pain sensation remains unknown. We have generated the first mouse model for the G1662S mutation by using homologous recombination in embryonic stem cells. The corresponding Scn10a<sup>G1663S</sup> mouse line has been analyzed for Scn10a expression, intraepidermal nerve fiber density (IENFD), and nociception using behavioral tests for thermal and mechanical sensitivity. The Scn10a<sup>G1663S</sup> mutants had a similar Scn10a expression level in dorsal root ganglia (DRG) to their wild-type littermates and showed normal IENFD in hindpaw skin. Mutant mice were more sensitive to touch than wild types in the von Frey test. In addition, sexual dimorphism was observed for several pain tests, pointing to the relevance of performing the phenotypical assessment in both sexes. Female homozygous mutants tended to be more sensitive to cooling stimuli in the acetone test. For heat sensitivity, male homozygous mutants showed shorter latencies to radiant heat in the Hargreaves test while homozygous females had longer latencies in the tail flick test. In addition, mutant males displayed a shorter reaction latency on the 54°C hot plate. Collectively, Scn10a<sup>G1663S</sup> mutant mice show a moderate but consistent increased sensitivity in behavioral tests of nociception. This altered nociception found in Scn10a<sup>G1663S</sup> mice demonstrates that the corresponding G1662 mutation of SCN10A found in SFN patients with pain contributes to their pain symptoms.

**Keywords:** SCN10A, sodium channel, small fiber neuropathy, nociception, neuropathic pain, statistical modelling

## INTRODUCTION

Neuropathic pain “arises as a direct consequence of a lesion or diseases affecting the somatosensory system” as defined by IASP (Treede et al., 2008; Woolf, 2010; Brouwer et al., 2015). It affects 7% of the general population and 40% of patients with neurological diseases (Bouhassira, 2019). Neuropathic pain is generally characterized by abnormal sensations such as hyperalgesia (increased response to painful stimuli) and allodynia (pain in response to a stimulus that does not normally provoke pain)

**OPEN ACCESS**

Edited by:

Jing Yao,
Wuhan University, China

Reviewed by:

Sara Gonzalez-Rodriguez,
University of Oviedo, Spain
Michael E. O'Leary,
Cooper Medical School of Rowan
University, United States

*Correspondence:

Yann Herault
herault@igbmc.fr

†These authors have contributed
equally to this work

**Specialty section:**

This article was submitted to
Pharmacology of Ion Channels and
Channelopathies,
a section of the journal
Frontiers in Pharmacology

Received: 20 September 2021

Accepted: 01 November 2021

Published: 01 December 2021

**Citation:**

Chidiac C, Xue Y, Muniz Moreno MdM, Bakr Rasheed AA, Lorentz R, Birling M-C, Gaveriaux-Ruff C and Herault Y (2021) The Human SCN10A<sup>G1662S</sup> Point Mutation Established in Mice Impacts on Mechanical, Heat, and Cool Sensitivity. *Front. Pharmacol.* 12:780132
doi: 10.3389/fphar.2021.780132

<footnote n="1">¹CNRS, INSERM Institut de Génétique et de Biologie Moléculaire et Cellulaire (IGBMC), Université de Strasbourg, Illkirch, France,</footnote>

<footnote n="2">²CNRS, INSERM, PHENOMIN-Institut Clinique de la Souris, Université de Strasbourg, Illkirch, France</footnote>
(Farrar and Portenoy, 2001) and is related to functional changes in primary afferents and sensitization within the central nervous system (CNS) (Brouwer et al., 2015; Bennett et al., 2019). Primary sensory neurons express unique repertoires of voltage-gated sodium channels. These channels are critical in initiating and propagating the all-or-none action potential to the CNS (Bennett et al., 2019). Tissue damage or injury can lead to changes in expression of sodium channel subunits in the sensory neurons, leading to an imbalance between excitatory and inhibitory signaling. These modifications will affect pain thresholds and activity resulting in hyper- or hyposensitivity to pain (Nassar et al., 2005; Colloca et al., 2017).

Inherited disorders of ion channels are known to account for a wide spectrum of human diseases called channelopathies. More than 1,000 disease-related mutations have been identified in the family of NAV channels (NAV1.1-NAV1.9) (Huang et al., 2017). Nowadays, about 800 variants have been described in the SCN10A gene encoding the NAV1.8 channel in the ClinVar NCBI database (https://www.ncbi.nlm.nih.gov/clinvar queried the August 26, 2021), and potentially pathogenic variants of the SCN10A gene were found in 3.7% of patients with small fiber neuropathy (SFN) (Eijkenboom et al., 2019). The missense point mutation c.4984G > A (p.Gly1662Ser) was identified in three SFN patients with years of pain history (Faber et al., 20121944; Han et al., 2014; Eijkenboom et al., 2019). The first patient was 24 years old and suffered from continuous tingling and cramps and reported more pain when exposed to cold temperature or exercising. She had abnormal sensitivity thresholds for warmth and cold stimuli as assessed by quantitative sensory testing (QST). The second patient was 62 years old, experienced continuous burning and stabbing pain, and was very intolerant to sheet. She showed hypoesthesia and reduced vibration sense in physical examination and abnormal warmth sensation levels revealed by QST (Han et al., 2014). The third patient was recently reported. She was 28 years old and showed abnormal temperature sensitivity thresholds (Eijkenboom et al., 2019). In functional studies, the SCN10A<sup>G1662S</sup> mutation was shown to be the first human disease-causing mutation that enhanced the generation of tetrodotoxin-resistant (TTX-R) resurgent currents (Xiao et al., 2019), impaired channel inactivation, and accelerated recovery from inactivation, which result in hyperexcitability of DRG neurons (Han et al., 2014).

With the aim to study the causality of the SNC10A G1662S mutation in pain, we have generated the corresponding mouse model. The comparison of human and mouse SCN10A protein sequences showed 82% homology with a resultant of one additional amino acid found in the mouse vs. human sequence before the G1662 residue, leading to a change in the position of the human G1662 to G1663 in mice (see **Supplementary Figure S1**). We have created the Scn10a<sup>G1663S</sup> mouse model corresponding to the G1662S mutation in the human gene and have investigated the p.G1662S genotype-phenotype link and in order to relate this mutation to SFN pain symptoms. The Scn10a transcripts are well expressed in the DRG of mutant mice as compared to wild-type controls. The mutant animals showed no alteration of growth, survival, and health state. The

Scn10a<sup>G1663S</sup> mutant mice of both sexes were evaluated for pain behavioral responses to noxious cold, heat, and mechanical stimuli. We found that Scn10a<sup>G1663S</sup> mutant mice displayed a moderate pain phenotype and identified the behaviors contributing the most to the pain phenotype using a novel statistical modeling analysis. This suggests that the SCN10A<sup>G1662S</sup> mutation identified in SFN patients with pain favors their pain symptoms.

# MATERIALS AND METHODS

## Animals

### Experimental Subjects and Ethical Approval

Animal experiments were performed after approval from the local ethical committee (Com'Eth, "Comité d'Ethique pour l'Expérimentation Animale IGBMC-ICS") with the agreement number 14691. The number of mice used for the experiments was adapted according to the 3R principles. Mouse breeding and behavioral experiments were conducted in SPF conditions at the ICS animal facility. Animals were housed under a 12-h/12-h light/ dark cycle (lights on at 7 a.m.) and 21 ± 1°C, 55 ± 10% humidity condition. Standard Chow (D04 for maintenance and D03 for breeding, SAFE, Les Oncins, France) and autoclaved tap water were provided ad libitum with two to five mice housed per cage. Prior to beginning behavioral experiments, all mice were habituated to the environmental conditions and to handling.

## Establishment of the Genetic Animal Model

The mouse line creation was done on the C57BL/6NCrl genetic background. Mutant mice carrying the G1663S mutation in the Scn10a gene corresponding to the G1662S mutation in the human gene (Supplementary Figure S1) were created by homologous recombination in ES cells. A targeted allele was engineered bearing two homology arms of the original Scn10a endogenous gene and a floxed Neo resistance auto-selection cassette. A 1.2-kb fragment encompassing Scn10a exon 28 (ENSMUSE00000582971) and 5' intronic sequence was amplified by PCR on ES cell C57BL/6N genomic DNA in two steps to allow the introduction of the GGC > AGC point mutation (leading to the G1663S mutation) and subcloned in an ICS proprietary vector. The 5' arm of homology (1.1 kb) which consisted of part of intron 27 was amplified by PCR and subcloned in a step1 plasmid to generate the final targeting construct. An RNA guide, gR68 (guide RNA sequence CTA ATAAGATATTCTGGGTA), targeting the intronic sequence at the site of insertion of the NeoR selection cassette was cloned in pX330 from Addgene (CRISPR/Cas9 plasmid). This vector was used to increase the efficiency of homologous recombination at the Scn10a locus.

The homologous recombination strategy is shown in **Figure 1A**. Both plasmids targeting plasmid and CRISPR/Cas9 plasmid were electroporated in C57BL/6N mouse ES cells (ICS proprietary line). After G418 selection, targeted clones were identified by long-range PCR and further confirmed by Southern blot with an internal (Neo) probe and a 5' external probe.
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

**FIGURE 1** | Generation and validation of the Scn10a<sup>G1663S</sup> mouse model. (A) Targeting strategy by homologous recombination. The scheme shows the Scn10a WT allele, the targeted allele, and the G1663S knock-in allele. The targeted allele was engineered bearing two homology arms of the original Scn10a endogenous gene and an auto-selection cassette. The 3' arm of homology (1.2 kb) consisted of a part of intron 28 and of a part of exon 28 bearing the S1663 mutation. The 3' arm was cloned into a targeting vector containing the loxP-flanked ("floxed") neo cassette. The 5' arm of homology (1.1 kb) consisted of part of intron 28 and was cloned into the vector containing the auto-selection cassette and the 3' arm. The final mutant allele was obtained after excision of the selection cassette. Light blue box: Exon 28 containing the locus G1663/S1663; the PM was indicated with a brown star; blue double arrow lines: indicate size of the arms; red triangles: LoxP sites; light gray box: neomycin selection cassette; green pentagon: Cre cassette. Different-colored lightning block arrows show the primers used for genotyping. The green rectangle indicates the primer located in the Cre cassette (Fcre) used for DNA sequencing. (B) Mutant allele sequence characterized through Sanger sequencing on a positive ES clone using a specific primer located in the Cre cassette (Fcre) amplifying only the mutant allele. The mutated nucleotide is indicated in red. (C) Example of genotyping results on electrophoresis gel showing three different genotypes. For all samples, PCR3 and 4 gave no bands indicating that the selection marker was excised. For sample #7106: PCR1 and 2 gave one band each at 363 and 236 bp, respectively. PCR5 gave one band of 194 bp. These reactions confirm that one LoxP site is still there and the Scn10a<sup>G1663S/G1663S</sup> genotype. For sample #7107: PCR1 and 2 gave two bands each at 363 and 250 bp and 236 and 123 bp, respectively. One band at 194 bp for PCR5. These PCR reactions confirm that one LoxP site is still there and the Scn10a<sup>G1663S/G1663S</sup> genotype. For sample #7108: PCR1 and 2 gave one band each at 250 and 123 bp, respectively. No band for PCR5, confirming the wild-type genotype. (D) Wild-type and G1663S allele mRNA expression in DRG of Scn10a<sup>G1663S</sup> heterozygous and homozygous mutant mice. Scn10a mRNA expression was normalized to Hprt expression. Homozygous mutant mice expressed Scn10a<sup>G1663S</sup> transcript at levels comparable to wt transcript expression, and heterozygous mice express half level of each transcript type. Females, *n* = 7–9/group; males, *n* = 7–9/group. One-way ANOVA analysis was followed by Sidak's multiple-comparison tests. The bars and *p* values indicate significant differences (*p* < 0.05) between wt and heterozygous mutant mice, and between heterozygous and homozygous mutant mice in Sidak's multiple-comparison tests. See **Supplementary Table S3** for detailed statistical analysis. (E) The G1663S mutant mice showed comparable IENFD compared to wt. Data are expressed as % of the mean of IENFD of each animal to the mean of wt of the same sex. Circles are used for females and squares for males. *n* = 4 mice per genotype and sex. The two-way ANOVA for the effects of genotype and sex as well as one-way ANOVA on grouped sexes showed no sex or genotype effect. Data are expressed as means ± SEM in all panels. See **Supplementary Table S4** for detailed statistical analysis.
We confirmed Neo cassette insertion by long-range PCR and verified that the construct was not integrated randomly in the genome. Clones with the correct genotype were injected into BALB/CN blastocysts, and resulting chimeric males having the correct targeted genotype and a color score >85% were bred with C57BL/6NCrl WT females for germline transmission. Heterozygous Scn10a<sup>+/G1663S</sup> male and female mice were crossed to generate the heterozygous Scn10<sup>+/+</sup> (= wt), heterozygous Scn10<sup>+/G1663S</sup> (het), and homozygous Scn10a<sup>G1663S/G1663S</sup> (homo) mice.

## Determination of Genotype

For genotype determination, 0.5 cm of mouse tails was harvested. Genomic DNA was extracted by DirectPCR Lysis Reagent Tail (Viagen Biotech, Los Angeles, CA, USA; Cat # 101-T) according to the manufacturer's instructions. Five hundred nanograms of genomic DNA extracted from WT or mutant mouse was PCR-amplified by 0.2 µl Taq DNA polymerase (Roche), 1.5 µl 10× PCR buffer with MgCl₂, 0.3 µl dNTP PCR-grade nucleotide (dATP, dCTP, dGTP, dTTP at 10 mM, Thermo Scientific, Waltham, MA, USA), 0.5 µl each primer at 10 µM, and H₂O in a total volume of 20 ml. Using a T100 thermocycler (Bio-Rad, Hercules, CA, USA), PCRs were subject to the following thermal conditions: a denaturation step at 94°C for 2 min followed by 40 cycles of 94°C for 30 s, a gradient of annealing temperatures between 55°C and 65°C for 45 s and 72°C for 1 min/kb, and a final elongation step for 5 min at 72°C. The PCR outcome was analyzed on a 1.5% agarose gel. To confirm the G1663S mutation, genomic DNA from mutated mice was PCR-amplified and confirmed by DNA sequencing. Genotyping of G1663S mutant mice consisted of five PCR reactions: PCR1 and PCR2 to confirm the excision of the selection marker using Ef/Er2 and Ef2/Er3 primers and giving bands of 4,484 and 4,357 bp for the recombinant allele, 363 and 236 bp for the PM allele, and 250 and 123 bp for the WT allele, respectively. PCR3 and PCR4 targeting the 5' and 3' sides of the selection marker, respectively, using Ef/Mq1r and Mq1f/Er primers giving respectively bands of 286 and 461 bp for the targeted allele and no bands for PM and WT alleles. The PCR5 specific for the LoxP site using Ef/Lxr primers gave a band of 194 bp for the targeted and PM alleles. Primer sequences and PCR reaction information are shown in **Supplementary Tables S1A,B**, and primers cutting sites are shown in the homologous recombination strategy in **Figure 1A**.

## Determination of Transcript Expression by Real-Time Droplet Digital PCR

Dorsal root ganglia were collected from WT and mutant mice and flash frozen in liquid nitrogen. All tissues were disrupted using Precellys® CK14 Lysing Kit in TRIzol Reagent and total RNA purified using RNeasy Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions. cDNA synthesis was performed using the SuperScript™ VILO™ cDNA Synthesis Kit (Invitrogen, Carlsbad, CA, USA). The Scn10a transcription level was determined by using ddPCR and Hprt as a standard gene. ddPCR was performed in 20-µl reactions containing 1× ddPCR Supermix for Probes (No dUTP), 125 nM of each probe,

250 nM specific primers, and 50 ng DNA according to the manufacturer's recommendations (PCR conditions: 10 min at 95°C, 40 cycles of 20 s at 95°C 30 s at 59.2°C, 2 min at 72°C; and 10 min at 98°C). Droplet generation and absolute droplet quantification were performed in a QX200 Droplet Digital PCR System (Bio-Rad) and analyzed by QuantaSoft Software (Bio-Rad). Probe and primer sequences are shown in **Supplementary Table S2**.

## Histological Analysis

### Tissue Preparation

Naïve adult mice (aged 23–28 weeks) were used for the histological analysis. Mice were anesthetized with ketamine/xylazine and perfuse intracardially with 40 ml of phosphate-buffered saline 0.1 M pH 7.4 (PBS) following 4% paraformaldehyde (PFA) in 0.1 M pH 7.4 PBS. Glabrous hindpaw skin samples (5 mm long) were carefully dissected out and postfixed for overnight at 4°C in 4% PFA in PBS, then tissues are removed from PFA and cryoprotected at 4°C in 30% sucrose in phosphate buffer (PB) for 3 days, embedded in optimal cutting temperature medium, frozen, and kept at -80°C. Glabrous skin longitudinal sections (30 µm thick) were cut with a cryostat and kept at -20°C.

## Immunohistochemistry

IHC of skin sections was performed with free-floating sections. Light antigen retrieval was performed in 0.1 M PB pH 7.3, 0.2% hydrogen peroxide, and 0.05% Triton X-100 for 25 min at room temperature (RT). Sections were then washed two times in PBS and incubated in blocking solution of Tris-buffered saline (TBS) with 0.05% Tween 20, 2% bovine serum albumin (BSA), and 2% normal donkey serum (NDS) for 1 h before applying a primary antibody overnight at 4°C. The following primary antibody (diluted in the blocking solution) was used: rabbit anti-PGP9.5 (1:500, ab108986, Abcam). For detection of the primary antibody, a secondary antibody raised in donkey and conjugated with Alexa-647 fluorophore was used (1:500, Molecular Probes, Thermo Fisher Scientific) for 1 h at RT. DAPI staining (1:2,000, Molecular Probes) was performed as the same time as the secondary antibody. Sections were then washed two times with PBS then placed on superfrost glass slides. One drop of Immu-Mount is added over the tissue sections, and a cover glass is placed over the slide.

## Image Acquisition and Analysis

N = 4 mice per genotype and sex were used for IENFD determination. Both hindpaw skin and four randomly selected sections per paw were analyzed. Images were visualized with Leica SP8 UV/Visible Laser Confocal Microscope using a ×20 oil objective. The ×20 resolution was achieved with a digital zoom factor of 3. Image acquisitions in the sequential mode (excitation intervals: 415–550 nm and 643–750 nm) were used for marker co-localization to avoid potential cross talk between the different fluorescence emissions. Images were acquired with the LCS (Leica, Wetzlar, Germany) software. The whole length of the skin sections was acquired as a tile scan, and an overview of the epidermis-dermis was reconstructed by acquisition of z-stacks
with a total depth of 30 µm. IENF were counted, and the length of the epidermal surface is measured using ImageJ software on digitized confocal images (Schneider et al., 2012). Counting on blinded samples was performed manually on five randomly chosen sections per paw per animal. Afferents crossing the dermal-epidermal junction were counted. Nerves that branch after crossing the junction were counted as a single unit. Nerves that split below the junction or appear to branch within the junction were counted as two units (Lauria et al., 2005). IENFD was determined by dividing the number of IENF by the total length of skin analyzed per section. We defined for each animal the IENFD as the mean of the IENFD calculated for the four sections per animal. The % of IENFD in mutants was calculated by dividing the mean of IENFD of each mouse by the mean IENFD of wt mice of the same sex.

## Behavioral Assessment of the Scn10a<sup>G1663S</sup> Mouse Model

A pipeline for behavioral assessment was defined and conducted on mutants and control littermates to investigate nociception and proprioception. All behavioral tests were conducted between 9:00 a.m. and 5:00 p.m. Animals were transferred to the experimental room 30 min before each experimental test. Body weights of wt and mutant mice were measured at the age of 7-11 weeks. The behavioral tests were carried on from the least noxious to the most noxious test in the following order: string test, crenellated bar, von Frey, Hargreaves plantar, acetone test, tail flick, tail pressure, cold plate, and hot plate. At least 2 days were kept between two consecutive tests. Females and males of the different genotypes were tested at 2 months of age. The experimenter was blind to mouse genotype.

The grip string test (in-house) was used to measure muscle strength. It was built as a horizontal wire placed at 40 cm above a table. The mouse was hanged by the forelimb to the wire, and the latency to hindlimb traction was measured. Three consecutive trials were done by 5-min intervals.

The notched/crenellated bar test was used to evaluate motor coordination and balance. The method described by Carter et al. was used (Carter et al., 2001). Briefly, mice were put on an elevated crenellated bar and had to traverse the distance from the beginning of the bar until the end where the home cage is placed. The time required to traverse the whole bar distance was recorded. In addition, the number of mistakes was defined as the times when mice were touching the gaps of the crenellated bar by their hindpaws while passing over the bar.

Tactile sensitivity was analyzed with the Von Frey test. One day prior to the test, mice were placed in small transparent plastic boxes over a mesh floor and were left there to acclimatize for 45 min. On the test day, mice were habituated to the box for 30 min before starting the test. Mechanical sensitivity was determined using different von Frey filaments of specific forces from 0.008 to 2 g (Bioseb, Vitrolles, France) and the “up-down” von Frey method (Weibel et al., 2013). Briefly, filaments were applied to the mid-plantar hindpaws starting with the filament of 0.4 g. Each trial was done by applying the filaments onto one hindpaw of mice placed consecutively and

then applying the filaments onto the other hindpaw of each mouse. The 50% mechanical threshold was calculated for each animal. The test was reproduced on two consecutive days with two trials per day.

The “Rodent pincher” test (Bioseb, France) was used to determine sensitivity to noxious mechanical stimuli, by applying a gradually increasing pressure to mice tail. This test allows calibrated forceps to induce quantifiable mechanical stimulation on a linear scale with a cutoff of 500 g. Two days prior to the Rodent pincher test, animals were habituated to be restrained in a 50-ml tube, with the whole tail exposed. Three measures were performed for each mouse, recording on three different locations of the tail. The mean of tail withdrawal threshold was calculated.

The acetone test was used to measure sensitivity to cool stimuli. Mice were placed and habituated in the same manner as for the von Frey test. Ten microliters of acetone was applied to the center of the plantar surface of each hindpaw, using a 50 µl Hamilton Gastight Syringe (Model 1705 TLL). Acetone was applied in three successive testing sessions for each paw. The interval between each application was at least 5 min. The duration of hindpaw withdrawal and flicking was measured during 30 s following acetone application. The values of hindpaw withdrawal and flicking duration displayed in the six trials were added.

The cold plate assay (Bioseb, France) was used to assess sensitivity to noxious cold temperatures. One day prior to the test, mice were acclimatized to the plate at RT, then on the test day, mice were tested on a 5°C cold plate. To avoid tissue damage, a maximum testing time of 5 min was set up. The mean numbers of hindpaw lifts and jumps were measured.

The Hargreaves test consisted of applying a radiant infrared heat stimulus to the hindpaws of mice to determine heat thresholds. One day prior to the test, mice were habituated to the apparatus by placing them in small transparent plastic boxes over a glass plate. Withdrawal latency was recorded on two consecutive days, with two trials per day. A trial was done by applying the heat stimulus on one hindpaw for all the mice, then applying it to the second hindpaw for all the mice, leaving at least 5 min between each hindpaw recording of the same mice. The mean of all measurements of both wt paws was calculated. The radiant heat was applied using a temperature ramp speed allowing to get a paw withdrawal latency of 7 s in wt mice.

Tail flick consisted in applying a focused beam of light to the mouse tail, to determine heat thresholds. Two days prior to the test, mice are habituated to be restrained in a 50-ml tube, with the whole tail exposed. The heat stimuli were applied at three different places of the tail tip, and the mean of the three withdrawal latencies was calculated.

A hot plate test was used as an integrated test to determine heat sensitivity. One day before the test, mice were habituated to the plate at RT. In our study, the mice were tested at three different temperatures on 3 separate days with specific cutoffs to avoid tissue damage: 3 min at 47°C, 1 min at 50°C, and 30 s 54°C. At each temperature, the latency to display first a reaction was recorded. In addition, we also measured coping reactions, a second line behavioral response to pain in addition to reflex responses, as reported recently by Huang et al. (Huang et al.,
2019). Mutant mice with ablated spinal preprotachykinin-
positive neurons showed withdrawal latencies comparable to
their wt littermates on the hot plate while they showed
reduced number of licking episodes compared to wts (Huang
et al., 2019), revealing the importance of scoring coping behaviors
as an additional information of heat pain sensitivity in this test. In
our study, the coping behaviors included lifting, flicking, licking
of the hindpaws, and jumping. The number of total reactions was
calculated per minute of test at each temperature.

**Statistical Analysis**

Results are expressed as mean ± SEM. Statistical analyses were
performed by using GraphPad Prism 9 software. All data sets
were tested for normality using the D'Agostino and Pearson
normality test. One-way ANOVA analysis was done when data
were found normally distributed, and Kruskal–Wallis analysis
was used when data had non-normal distribution. For transcript
expression level analysis, the genotype effect in each sex was
analyzed using one-way ANOVA analysis followed by Šidák's
multiple-comparison test when appropriate. The analysis of
behavioral tests results was performed by one-way ANOVA or
Kruskal–Wallis analyses to test for genotype difference on
grouped sexes and separated sexes, followed by *post-hoc*
analysis when appropriate. A *p*-value < 0.05 was considered
statistically significant. Supplementary statistical analyses were
done to identify the variables contributing the most to the
genotype or sex discrimination using the Gdaphen R pipeline.
The principle of Gdaphen analysis is described in the
**Supplementary Material**.

RESULTS

Generation and Characterization of the
Scn10a<sup>G1663S</sup> Mouse Model

Initially, we tried to create the Scn10a<sup>G1663S</sup> mouse model with the CRISPR-Cas9 innovative technology. We tried two different strategies that were unsuccessful due to the high similarity between Scn10a and other Scn genes, especially the Scn5a gene that shares the highest similarity with Scn10a and is expressed on the same chromosome. We moved then to homologous recombination in ES cells to generate the Scn10a<sup>G1663S</sup> mouse line (Figure 1A). Briefly, the G1663S-targeting vector was electroporated into C57BL/6N ES cells, together with the plasmid vector containing a guide RNA named gR68. The guide RNA gR68 in association with the Cas9 protein specifically generates a double-strand break in the intronic region close to the targeted region. The intronic sequence, in contrast to the exonic sequence, does not show any similarity to intronic sequences of any other Scn genes. The CRISPR/Cas9 double-strand break driven by the guide RNA gR68 forced the ES cells to repair, and the presence of the targeting construct with the selection pressure (NeoR) increased the homologous recombination rate (even if the size of the homology arms is much less important than with classical ES cell targeting). Among positive ES clones, one clone was selected. The genomic DNA from the mutant allele of this

clone was PCR-amplified with a forward primer located in
the Cre cassette (Fcre) (Figure 1A) and a 3' external reverse
primer and was then sequenced. Sequencing results
demonstrated the presence of the mutation in this positive
ES clone (Figure 1B). Heterozygous Scn10a<sup>+/G1663S</sup> mice were
crossed to generate the Scn10a<sup>G1663S/G1663S</sup> mice and their
control wt mice. An example of genotyping results is shown in
Figure 1C. The Scn10a transcript expression level was
evaluated in the DRG of G1663S and control mice by RT-
ddPCR. One-way ANOVA revealed a genotype effect for each
allele-derived (wt or mutant G1663S) transcript expression in
grouped or separate sexes. When ratios of Scn10a<sup>+</sup> and
Scn10a<sup>G1663S</sup> transcripts were summed, no effect on
genotype or sex was detected by the two-way ANOVA.
Both transcripts were equally expressed in the DRG of wt
or G1663S/G1663S mutant mice respectively and in both
sexes, indicating that the mutation did not alter the Scn10a
expression. Heterozygous mice expressed about 50% of each
allelic transcript in DRG (Figure 1D and Supplementary
Table S3). SFN patients with heterozygous c.4984G > A,
p.G1662S mutation showed a normal IENFD. We have
quantified IENFD in adult control and mutant mice. No
genotype or sex effect was detected in the two-way ANOVA
test. The one-way ANOVA on grouped sexes also showed no
difference between mutant and wt animals (Figure 1E and
Supplementary Table S4).

Normal Body Weight and Proprioception in
the Scn10a<sup>G1663S</sup> Mice

The Scn10a<sup>G1663S</sup> mutation did not affect the body weight of the
mutant mice. Heterozygous +/G1663S and homozygous
G1663S/G1663S mice showed a normal body weight
compared to wt littermates (Figure 2A and Supplementary
Table S5). To evaluate whether the mutation alters the
proprioception of mutant mice, string test and crenelated bar
were performed. String test latency analysis revealed a sex effect
(two-way ANOVA F<sub>1, 105</sub> = 10.76, *p* = 0.0014) but no genotype
effect. Similarly, no genotype or sex effect was detected in the
crenelated bar test. Heterozygous and homozygous mutant mice
behaved in the same manner than control mice in these two tests
(Figure 2B-D and Supplementary Table S5). We concluded
that the presence of one or two alleles with the G1663S mutation
did not change body weight and proprioception in the
mutant mice.

Enhanced Sensitivity to Mechanical Stimuli
in the Scn10a<sup>G1663S</sup> Mice

The oldest SFN patient with heterozygous G1662S mutation developed sensitive skin and intolerance to sheets over her feet with continuous pain (Han et al., 2014; Eijkenboom et al., 2019). We used the von Frey test to assess touch allodynia and the rodent pincher test to evaluate sensitivity to noxious mechanical stimuli in the Scn10a<sup>G1663S</sup> mutant mice. In the von Frey test, a genotype effect was present (two-way ANOVA F<sub>2, 93</sub> = 3.461, *p* = 0.035). The one-way ANOVA showed a genotype effect on grouped sexes
<figure><img src="image_9.png" /></figure>

**FIGURE 2** | *Scn10a*<sup>G1663S</sup> mice show normal healthy conditions and proprioception capacities. (A) Body weight of *Scn10a*<sup>G1663S</sup> mice measured at the age of 7–11 weeks. Females, *n* = 13–15/group; males, *n* = 9–16/group. The one-way ANOVA for genotype on separated sexes showed no genotype effect. (B) Muscle strength was evaluated by string test. Females, *n* = 18–20/group; males, *n* = 15–21/group. A sex effect was detected in the string test. The Kruskal–Wallis test on separated sexes showed no genotype effect. (C–D) Motor coordination and balance were evaluated by the latency and the number of mistakes in the crenellated bar test. Females, *n* = 17–21/group; males, *n* = 16–21/group. The Kruskal–Wallis tests on separated or grouped sexes showed no sex or genotype effect. Data are expressed as means ± SEM in all panels. See **Supplementary Table S5** for detailed statistical analysis.

where heterozygote and homozygote mice showed a higher sensitivity to mechanical stimuli compared to their wt littermates (Figure 3A and Supplementary Table S6). This did not reach significance when males and females were analyzed separately. In the rodent pincher test, two-way ANOVA analysis resulted in a sex effect (two-way ANOVA F₁, ₁₀₆ = 4.186, *p* = 0.043). In addition, there was a tendency for higher sensitivity in the mutants that was above the significance threshold (Figure 3B and Supplementary Table S6). Therefore, *Scn10a*<sup>G1663S</sup> mutant mice developed mechanical allodynia.

## Sensitivity to Cooling and Cold Stimuli in the *Scn10a*<sup>G1663S</sup> Mice

One of the SFN patients with G1662S mutation reported that cold temperature aggravated her complaints and QST revealed abnormal thresholds for cold sensation (Han et al., 2014). We performed the acetone test to investigate the sensitivity of the *Scn10a*<sup>G1663S</sup> mice to cooling stimuli (12°C–15°C) and cold plate at 5°C for cold stimuli. In the acetone test, the duration of withdrawal and flicking reactions revealed an interaction between genotype and sex (two-way ANOVA F<sub>2, 107</sub> = 4.336, *p* = 0.015). The Kruskal–Wallis analysis showed a

higher behavioral response to the acetone stimulus in mutant females while mutant males reacted comparably to wt males (Figure 3C and Supplementary Table S6). In the cold plate test, the number of hindpaw lifts and jumps showed a sex effect (two-way ANOVA F<sub>1, 107</sub> = 21.61, *p* < 0.0001) with fewer reactions in males. Mutant and control mice of both sexes had comparable reactions on the cold plate (Figure 3D and Supplementary Table S6). Therefore, *Scn10a*<sup>G1663S</sup> mutant mice showed a higher sensitivity to cool, but not cold temperatures.

## Enhanced Heat Pain Sensitivity in the *Scn10a*<sup>G1663S</sup> Mice

An abnormal threshold for warmth sensation was also an important clinical feature of SFN patients with the G1662S mutation (Han et al., 2014; Eijkenboom et al., 2019). *Scn10a*<sup>G1663S</sup> mutant mice were tested for their sensitivity to radiant heat stimuli in the Hargreaves and tail flick tests (Figure 4). The statistical analysis of the latencies in the Hargreaves and tail flick tests identified an interaction between sex and genotype (two-way ANOVA F<sub>2, 113</sub> = 3.808, *p* = 0.025 and F<sub>2, 108</sub> = 3.345, *p* = 0.039 respectively). In the Hargreaves test, the
<figure><img src="image_10.png" /></figure>

**FIGURE 3** | Scn10a<sup>G1663S</sup> mice show enhanced pain sensitivity to mechanical and cool stimuli. (A) G1663S mutant mice displayed hypersensitivity to von Frey filaments. Females, *n* = 17–18/group; males *n* = 14–17/group. The one-way ANOVA showed a general genotype effect in sex-grouped animals. Unpaired two-tailed t-tests showed a genotype effect for wt vs. heterozygous and wt vs. homozygous mutant sex-grouped mice. The 2 bars and *p* values are for genotype difference in grouped females and males, wt vs. heterozygous and wt vs. homozygous mutants. (B) Noxious mechanical sensitivity was analyzed by the tail pressure test. Females *n* = 18–20/group; males *n* = 16–20/group. Kruskal–Wallis analysis indicated a trend for genotype effect in sex-grouped animals. The bar and *p* value represent the overall genotype effect in grouped females and males. (C) Sensitivity to cool temperature was analyzed by the acetone test. Females *n* = 19–21/group; males *n* = 16–20/group. Mutant females tended to display longer duration of paw reactions following acetone application. Kruskal–Wallis analysis. The bar and *p* value represent the overall genotype effect in females. (D) G1663S mutant mice showed no phenotype for cold response on the 5°C cold plate. Females *n* = 18–21/group; males *n* = 17–20/group. Kruskal–Wallis analysis showed no genotype effect in either sex-grouped or sex separated animals. Data are expressed as means ± SEM in all panels. See **Supplementary Table S6** for detailed statistical analysis.

mutation lowered response latency, with homozygous males withdrawing their hindpaw earlier than wt males (Figure 4A and Supplementary Table S7). In the Tail Flick assay, female mutant mice showed a general genotype effect (Figure 4B and Supplementary Table S7).

The hot plate test is classically used to determine heat sensitivity. Latencies to first hindpaw reaction and to first jump reaction are important variables in this assay. Huang and others introduced the scoring of the number and duration of licking and jumping as outcomes of coping behaviors (Huang et al., 2019). Accordingly, we scored both the latency to first hindpaw reaction at 47°C, 50°C, and 54°C and the number of coping reactions (Figures 4C–E). The latency to first hindpaw reaction on the hot plate showed a sex effect at 54°C (two-way ANOVA F₁, 108 = 5.703, *p* = 0.019), and a genotype effect was detected at 54°C for the homozygous males (Figure 4E and Supplementary Table S7). For the coping reactions, a tendency for genotype effect was recorded at 47°C although it did not reach significance (Figure 4C and Supplementary Table

S7). Therefore, Scn10a<sup>G1663S</sup> mutant male mice showed an increased sensitivity to heat.

## Gdaphen Analysis for the Identification of the Variables Contributing the Most to the Genotype or Sex Discrimination

We sought to identify which variables contribute most to genotype and sex discrimination, by assessing how well each variable could discriminate between the three genotypes (wt, heterozygous, and homozygous) and between the two sexes. By performing the Gdaphen analysis, we considered 14 variables in total, genotype, sex, and 12 behavioral variables (see Supplementary Method). Among these 14, 9 variables were detected as the most relevant (>30%) to discriminate between the three mouse genotypes. These variables are sex, coping reactions on the 47°C and 54°C hot plate, coping reactions and latency on the 50°C hot plate, acetone test, paw lifts and jumps on a 5°C cold plate, von Frey, and tail pressure.
<figure><img src="image_11.png" /></figure>

**FIGURE 4** | *Scn10a*<sup>G1663S</sup> male mice show enhanced pain sensitivity to heat stimuli. (A) Male mutants show increased heat sensitivity in the Hargreaves test. Females *n* = 18–21/group; males *n* = 18–22/group. One-way ANOVA indicated a genotype difference in males and Sidak's multiple-comparison test showed a significant difference in wt vs. homozygous males. The bar and *p* value represent the genotype effect in wt vs. homozygous males. (B) Mutant females had increased latency in the tail flick test. Females *n* = 19–20/group; males *n* = 16–20/group. The one-way ANOVA showed an overall genotype difference in females. The bar and *p* value represent the global genotype effect in females. No genotype effect could be seen in wt vs. heterozygous or wt vs. heterozygous females. (C–E) Hot plate test. The latency to the first hindpaw reaction is shown in left panels and the number of coping reactions in right panels. (C) Mutant mice showed no phenotype on the 47°C hot plate. *n* = 17–22/group; males *n* = 14–19/group. Kruskal–Wallis analysis on sex-separated or grouped animals showed no genotype difference for first hindpaw sign latency and a trend for coping reactions. The bar and *p* value represent the global genotype effect for coping reaction in sex-grouped animals. (D) Mutant males or females showed no genotype difference on the 50°C hot plate. Females *n* = 19–21/group; males *n* = 16–20/group. One-way ANOVA on sex-separated or grouped animals. (E) The male mutants were more sensitive than their wt counterparts on the 54°C hot plate. Females *n* = 17–22/group; males *n* = 16–20/group. Kruskal–Wallis analysis showed a global genotype effect in males. The bar and *p* value represent the global genotype effect in males. Data are expressed as means ± SEM in all panels. See **Supplementary Table S7** for detailed statistical analysis.

Generalized linear model (GLM) (Figure 5A and Supplementary Figure S2.1A) and random forest (RF) (Figure 5A and Supplementary Figure S2.1B) classifiers identified von Frey as the most important variable discriminating between the three genotypes. The RF analysis also identified acetone, coping reactions, and latency on the 50°C hot plate and tail pressure

as major explanatory variables. For sex discrimination, GLM (Figure 5B and Supplementary Figure S2.2A) and RF (Figure 5B and Supplementary Figure S2.2B) identified tail pressure, paw lifts and jumps on the 5°C cold plate, and coping reactions on the 50°C hot plate as the most important variables. GLM also identified acetone as an important variable, and RF also
<figure><img src="image_12.png" /><figcaption>**FIGURE 5** | Gdaphen analyses to identify the variables contributing more to genotype and sex discrimination considering the three genotypes together, based on phenotypic data. (A) Identification of the power of each explanatory phenotypic variable to the genotype discrimination. The explanatory variables selected were the ones known to contribute more than a 30% to the genotype discrimination. Those variables were identified using a principal component analysis (PCA) after applying a multifactor analysis of mixed data (MFA) implemented using MFAmix function from the PCAmixdata R package (https://arxiv.org/abs/1411.4911). The relevance of those variables to the genotype discrimination was analyzed using two different statistical classifiers: Generalized linear models, noted as GLM taken from the caret R package, and Random forest, noted RF taken from the caret R package (Max Kuhn (2020). Caret: Classification and Regression Training. R package version 6.0-86. https://CRAN.R-project.org/package=caret)). All measures of importance are scaled to have a maximum value of 100 in the variable contributing the most to the discrimination. With GLM and RF, we can identify the importance of each variable to discriminate each genotype. (B) Identification of the power of each explanatory phenotypic variable known to contribute more than a 30% to the genotype discrimination, analyzed to see the existence of different responses based on the sex using the two different statistical classifiers. All measures of importance are scaled to have a maximum value of 100 in the variable contributing the most to the discrimination. (C) 3D-PCA plots showing the individual animals clustering on the 3D space based on the PCA analyses performed with all the phenotypic variables and colored based on genotype and sex. On the left panel showing the three genotypes and two-sex individuals, on the middle plot showing only male data, and on the right plot only female</figcaption></figure>

(Continued)
**FIGURE 5** | data. (D) Qualitative variable discrimination component map. Showing the distribution in 2D space of the qualitative variable coordinates calculated based on the PCA analysis performed after the MFA implemented using MFAmix function from the PCAmixdata R package. (E) Identification of correlation between each principal component of the separated analyses of each group of variables and the principal components (dimensions) calculated by the MFAmix function. By the arrow length, we could identify the contribution of each grouped variable to the discrimination on each dimension. Arrows that follow similar trajectories contribute to the separation of the individual data in the same dimensions.

identified coping reactions on the 47°C hot plate as an important one. Based on the qualitative variable discrimination 2D component map, sex had a stronger effect over genotype for discrimination between groups. Indeed, heterozygous and homozygous animals are more strongly separated from wt in dimensions 2 and 3, than in dimension 1; males and females are highly separated in dimension 1, but not in the other dimensions (Figure 5D). The multifactor analysis of mixed data (MFA) was used to assess the importance of each test after grouping the variables recorded in each test. By examining the correlation between the location of each grouped variable on the principal components, we show that von Frey and coping reactions on the 47°C hot plate contribute the most to genotype discrimination in dimensions 1 and 2 and dimensions 2 and 3. Tail pressure was also identified in dimensions 1 and 2 and paw lifts and jumps on the 5°C cold plate in dimensions 2 and 3. Coping reactions on the 50°C and 54°C hot plate and paw lifts and jumps on the 5°C cold plate were found to contribute the most to genotype discrimination in dimensions 1 and 3, with the strongest contribution in dimension 1 (Figure 5E).

DISCUSSION

The SCN10A<sup>G1662S</sup> mutation was described in three patients with painful SFN. We have generated a mouse model harboring this mutation using homologous recombination in ES cells and characterized this mouse line to understand the relationship between this mutation and the pain phenotype on the patients. The Scn10a mRNA is known to be expressed in DRGs that contain cell bodies of sensory neurons; thus, we studied the impact of the mutation on Scn10a transcript expression in DRGs. The Scn10a<sup>G1663S</sup> mutation did not alter the Scn10a transcript expression level in DRG in mutant mice of both sexes. IENFD was found unchanged in hindpaw skin of G1663S/G1663S mice. Finally, we investigated the sensitivity to mechanical, heat, and cold stimuli of the mutants in both sexes. Our results collectively indicate that the Scn10a<sup>G1663S</sup> mutant mice have an enhanced pain sensitivity toward mechanical, heat, and cold stimuli in some behavioral tests, which correlate with the pain symptoms of the SFN patients carrying the same mutation.

Normal IENFD in the Scn10a<sup>G1663S</sup> Mice

IENFD was normal in Scn10a<sup>G1663S</sup> in adult mutant mice, similarly to the SFN patients with SCN10A<sup>G1662S</sup> mutation that had normal IENFD values compared to age-matched and gender-matched normative values. However, patients carrying other point mutations in SCN10A have shown a loss of IENFs; thus, the loss of IENFs appears to be mutation dependent. A father-son pair harboring L554P mutation in the SCN10A gene has a complete depletion of IENFs and severe reduction in density of dermal nerve bundles, which appeared fragmented due to axonal degeneration (Faber et al.,

2012; Faber et al., 1944). Other SFN patients with I1706V and D1639N mutations in the SCN10A gene also showed a remarkably decreased IENFD (Huang et al., 2013; Dabby et al., 2016). Indeed, the study by Bianca de Greef et al. indicated that, out of 921 patients with SFN, 614 had normal IENFD with abnormal thermal thresholds (de Greef et al., 2018), indicating that abnormal IENFD is a secondary criterion for painful SFN. The normal IENFD found in Scn10a<sup>G1663S</sup> mutant mice therefore suggests that in mice also, Scn10a mutants can show normal IENFD and abnormal pain responses.

Enhanced Sensitivity to Mechanical Stimuli
in the Scn10a<sup>G1663S</sup> Mice

The enhanced sensitivity to mechanical stimuli in the Scn10a<sup>G1663S</sup> mice can be compared to other genetic mouse models for Scn10a. The Possum mice harboring a Scn10a<sup>T790A</sup> gain-of-function mutation were identified in an N-ethyl-N-nitrosourea-induced mutagenesis screen. Amanda Blasius and others reported abnormal mechanical behavior in Possum mice. Homozygous Possum mice showed no alteration in sensitivity to touch stimuli in the von Frey test while they showed increased sensitivity when a needle prick was applied to hindpaw (Blasius et al., 2011; Garrison et al., 2014). The Scn10a knockout (KO) mouse model generated by Akopian and others (Akopian et al., 1999) was studied for pain behaviors in many studies (see (Xue et al., 2021) for a recent review). The KO of the Scn10a gene did not impact on touch stimuli as assessed by the von Frey test (Akopian et al., 1999; Kerr et al., 2001; Laird et al., 2002; Nassar et al., 2005; Leo et al., 2010; Shields et al., 2012; Minett et al., 2014a) while it showed reduced pain sensitivity toward noxious mechanical stimuli as assessed by the tail pressure test (Akopian et al., 1999; Nassar et al., 2005; Minett et al., 2014b). Those findings suggested that the presence of Scn10a is important for sensitivity to noxious mechanical stimuli, but not to light touch. A third mouse model was developed where the Cre gene was knocked into Scn10a exon 1, the Scn10a<sup>Cre/+</sup> hemizygous mice. These mice had no alteration in response to von Frey filaments (Nassar et al., 2005; Stirling et al., 2005), which is in agreement with the results obtained with Scn10a KO mice. In addition, these mice behaved similarly to their wt littermates in the tail pressure test (Agarwal et al., 2004; Stirling et al., 2005), showing that inactivation of one Scn10a allele does not impair mechanical sensitivity. Our results show a reduced mechanical threshold of the Scn10a<sup>G1663S</sup> mice in the von Frey test. A phenotype in this test has not been previously reported in the other Scn10a mutant lines (Akopian et al., 1999; Kerr et al., 2001; Laird et al., 2002; Agarwal et al., 2004; Nassar et al., 2004; Nassar et al., 2005; Stirling et al., 2005; Leo et al., 2010; Blasius et al., 2011; Shields et al., 2012; Minett et al., 2014a; Garrison et al., 2014) which displayed altered sensitivities to noxious mechanical stimuli (Akopian et al., 1999; Minett et al., 2014b). With the help of Gdaphen analysis, we identified sensitivity to touch
(von Frey filaments) as the most important phenotype contributing to increased pain sensitivity in *Scn10a*<sup>G1663S</sup> mice. The increased sensitivity to touch in the mutant mice concurs with the *SCN10A*<sup>+/G1662S</sup> SFN patient reports on sheet intolerance. The increase in touch sensitivity in *Scn10a*<sup>G1663S</sup> mice is also an interesting finding for preclinical research on pain as the von Frey test is the most used test for assessing touch sensitivity in experimentally induced pain models in rodents. The higher touch sensitivity in mutant mice is consistent with the mechanical allodynia that is usually observed in preclinical models of neuropathic pain and could be interpreted to suggest an increase in the excitability of mechanical nociceptors.

## Sensitivity to Cooling and Cold Stimuli in the *Scn10a*<sup>G1663S</sup> Mice

Previous works reported that *Scn10a* mutant mice show alterations in the response to extreme cold as assayed by cold plate behaviors in *Scn10a* mutant mice. The Possum *Scn10a*<sup>T790A</sup> mice showed increased reactions on the cold plate at -1°C (Blasius et al., 2011). Collectively, previous studies with Cold Plate and acetone tests on *Scn10a* KO mice showed that *Scn10a* loss of function induces a loss of sensitivity to extreme cold but not to cool and mild cold (Zimmermann et al., 2007; Minett et al., 2012; Minett et al., 2014a; Luiz et al., 2019). More precisely, *Scn10a* KO mice were less sensitive than wt when nocifensive responses were measured on the cold plate at 0°C (Zimmermann et al., 2007). *Scn10a*<sup>Cre/Cre</sup> mice presented a reduced pain phenotype by jumping later than wt mice on the -5°C cold plate and no phenotype when forepaw lifting on the cold plate was assessed at 10°C and 5°C (Luiz et al., 2019). Our study shows no difference in the mean number of hindpaw lifts and total jumps on the cold plate at 5°C. In the future, further investigations using colder temperatures on the plate may reveal an influence of the *Scn10a*<sup>G1663S</sup> mutation on sensitivity to extreme cold as already shown in previous studies on other *Scn10a* mouse models (Zimmermann et al., 2007; Blasius et al., 2011; Luiz et al., 2019). Our findings in the acetone test suggest that the increase in *Scn10a* excitability induced by the G1663S mutation (Huang et al., 2013; Han et al., 2014) alters the sensitivity to cooling while the loss of function of *Scn10a*, investigated with the *Scn10a* KO mice, did not impact this sensitivity (Minett et al., 2012; Minett et al., 2014a; Luiz et al., 2019). Interestingly, G1663S/G1663S males did reactions to acetone similar to wt males while G1663S/G1663S females showed a tendency for longer duration of reactions. Thus, sex appears to influence the effect of the mutation on cooling sensitivity. Our results on increased sensitivity to cooling in the mutant females concur with the aggravated complaints of G1662S SFN patients when exposed to cold and with abnormal thresholds for cold sensation in the QST.

## Sensitivity to Heat Pain in the *Scn10a*<sup>G1663S</sup> Mice

In the Hargreaves test, the *Scn10a* KO animals displayed analgesia when a slow ramp for radiant heat was used (13-15 s latency in wt mice) (Akopian et al., 1999; Minett et al., 2014b) while no phenotype was found when using a faster heat ramp (7-8 s latency in wt mice)

(Nassar et al., 2005; Leo et al., 2010; Shields et al., 2012), indicating that the use of a slow temperature ramp allowed to detect analgesia in the *Scn10a* KO mice (Xue et al., 2021). In our study, the G1663S male animals were more sensitive in the Hargreaves test when using a heating ramp that led to a withdrawal latency of 7 s for wt mice similar to the latency described in the fast-ramp papers (Nassar et al., 2005; Leo et al., 2010; Shields et al., 2012). Therefore, in a fast-ramp condition previously reported to lead to no phenotype in the KO mice, we could evidence the increased sensitivity of G1663S male animals to radiant heat applied to the paw. In the tail flick assay, sex-grouped mice showed no phenotype, while mutant females had an overall higher reaction latency compared to control mice. An increased latency was found in *Scn10a* KO mice in this test under a heating condition leading to a longer latency in wt mice than in the present study, and with mice of unknown sex (Akopian et al., 1999). Also, *Scn10a* KO mice showed no phenotype in the tail immersion test (Akopian et al., 1999; Nassar et al., 2005; Leo et al., 2010). Altogether, the *Scn10a*<sup>G1663S</sup> line did not show increased sensitivity in the tail flick assay in the present experimental condition, and slower heating ramps may be further assayed. The sex difference for the mutation effect in the Hargreaves and tail flick tests could also be explained by sex differences in paw and tail skin innervation.

In the hot plate assay, a shorter latency to the first hindpaw reaction was found for G1663S male mice at 54°C, while no phenotype was found at the lower 47°C and 50°C temperatures. The Possum mice showed no alteration on the hot plate at 52°C (Blasius et al., 2011), and their sensitivity at other temperatures or for other heat tests was not reported. The *Scn10a* KO mice showed no phenotype in the hot plate assay when the latency to first hindpaw withdrawal reaction was recorded at 45°C, 50°C, 55°C, and 60°C (Akopian et al., 1999; Nassar et al., 2005; Leo et al., 2010). Also, *Scn10a*<sup>Cre/+</sup> mice had normal responses in the hot plate test at 50°C and 55°C (Agarwal et al., 2004; Nassar et al., 2005; Stirling et al., 2005).

Altogether, our findings on *Scn10a*<sup>G1663S</sup> mice in the three tests for heat sensitivity indicate that the *Scn10a*<sup>G1663S</sup> mutation leads to increased sensitivity to strong heating of the paw in males. This is in global agreement with previously published results on mutant mice in which *Scn10a* has been inactivated (Akopian et al., 1999; Nassar et al., 2005; Leo et al., 2010; Minett et al., 2014b; Xue et al., 2021), although the sex effect was not fully explored in this studies.

## Sex Effect on the *Scn10a*<sup>G1663S</sup> Mouse Behaviors

In our study, we assessed pain behaviors of *Scn10a*<sup>G1663S</sup> mice on both sexes and analyzed the data of each sex separately and sex-grouped. The results suggest that there are specific sex differences in *Scn10a*<sup>G1663S</sup> mutants. Mutant females showed increased reactions in the acetone test while mutant males had an increased sensitivity in the Hargreaves test and on the 54°C hot plate. The sex-grouped analysis allowed to detect an increased sensitivity to mechanical stimuli. This suggests that the effect of the mutation is influenced by the sex. The three SFN patients carrying the *SCN10A*<sup>G1662S</sup> mutation were women (Han
et al., 2014; Eijkenboom et al., 2019), whose symptoms globally correlate with the pain phenotype found in Scn10a<sup>G1663S</sup> mutant animals. Nevertheless, no SFN male patients with the SCN10A<sup>G1662S</sup> mutation have been reported until now. There are multifactorial and complex influences underlying the sex differences in pain in the clinical setting. Pain is reported more frequently by women than men. Healthy women have significantly lower pain thresholds for electrical, pressure, and thermal stimuli (Nasser and Afify, 2019). Sex hormones and sociocultural factors are known to impact men and women pain experience (Nasser and Afify, 2019). Women manifest their pain with stronger emotion and evoke pain distress, while men seem to withstand pain and less eager to report pain (Nasser and Afify, 2019). In addition, pain thresholds of both male and female rodents can be affected by the sex of the examiner (Kállai et al., 2004; Gijsbers and Nicholson, 2005). Also, sex differences in opioid-induced analgesia were reported (Niesters et al., 2010; Mogil, 2020). In future studies, it would be interesting to evaluate further mutant mice of both sexes for their response to inflammation-induced pain and for analgesia induced by opioids, gabapentinoids, and anti-depressants.

In recent years, more studies have used animals of both sexes and reported sex-dependent differences or no differences in nociceptive behaviors in rodent models (Mogil, 2020). In addition, the genetic background may impact on sex differences or no differences in pain. In the last decade, a neuroimmune mediation of pain has also been discovered, with the preferential involvement of spinal microglia in pain hypersensitivity in male rodents (Rosen et al., 2017; Sorge and Totsch, 2017; Mogil, 2020). The sex differences found in behavioral responses in Scn10a<sup>G1663S</sup> mutants contribute to complete the list of genes and mutations for which a sex effect has been detected.

## CONCLUSION AND PERSPECTIVES FOR
SCN10A<sup>G1663S</sup> MICE

Globally, the studies on Scn10a genetic rodent models have highlighted the important role of the NAV1.8 channel including the SCN10A protein in nociception (Xue et al., 2021). Also, the Nav beta4 subunit has been shown to contribute to sensory neuron excitability *in vitro* (Xiao et al., 2019). Still, rodent models of the human SCN10A mutations that were identified in SFN patients were needed to study those mutations *in vivo*. We have successfully created the mouse model for SCN10A<sup>G1662S</sup> mutation identified in SFN patients by using homologous recombination in ES cells. Scn10a<sup>G1663S</sup> mutant mice revealed a pain phenotype toward non-noxious and noxious mechanical stimuli, cooling stimuli, and some forms of heat stimuli. Moreover, by the help of the Gdaphen statistical analysis, we could identify which nociception tests allowed to discriminate the most mutants from wt control animals. Hereafter, other devices could also be used to assess thermal preference as complementary tests for heat sensitivity, such as the two-plate choice test (Bautista et al., 2007), the thermal gradient test (Bohic et al., 2020), or the operant plantar thermal assay (Reker et al., 2020). Our results show that Scn10a<sup>G1663S</sup> mutant animals concur in pain features with the SFN patients

carrying the SCN10A<sup>G1662S</sup> mutation, such as sensitive skin and abnormal thresholds for warm and cold stimuli as revealed by QST. In the future, spontaneous pain may also be assessed in the mutant mice. Different paradigms could be used to evaluate spontaneous pain, including weight bearing or conditioned place preference for the chamber associated with analgesic administration (Barrot, 2012; Tappe-Theodor et al., 2019).

The development of genetic rodent models recapitulating the SCN mutations found in SFN patients is a key step to explore pain mechanism involving sodium channels and the pillar to develop convenient novel analgesics. The Scn10a<sup>G1663S</sup> mouse model opens the door for further investigations in the aim of complementing the knowledge about the role of NAV1.8 channel in painful SFN.

## DATA AVAILABILITY STATEMENT

The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.

## ETHICS STATEMENT

The animal study was reviewed and approved by the ethical committee (“Com’Eth, Comité d’Ethique pour l’Expérimentation Animale” IGBMC-ICS, licence N° 17) with the agreement number 14691.

## AUTHOR CONTRIBUTIONS

All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.

## FUNDING

This work has been funded by the Molecule-to-Man Pain Network, a European Commission Multi-Center Collaborative Project through the European Union’s Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement No. 721841-PAIN-Net. This work has been supported by the National Centre for Scientific Research (CNRS), the French National Institute of health and medical research (INSERM) and the University of Strasbourg (Unistra). The French state funds through the “Agence Nationale de la Recherche” under the frame of Programme Investissements d’Avenir labelled ANR-10-IDEX-0002-02, ANR-10-LABX-0030-INRT, ANR-10-INBS-07 PHENOMIN.

## ACKNOWLEDGMENTS

We wish to thank Pr Giuseppe Lauria and all PAIN-Net members for fruitful discussions. We thank Laurence Schaeffer and Valerie
Erbs from the genetic engineering team, Sylvie Jacquot and Amelie Jeanblanc from the genotyping team, and Marie Wattenhofer-Donze from the ES cells team at ICS for their help to create and genotype the mutant mouse model. We thank Loïc Lindner and Pauline Cayrou for their help in ddPCR design and training. We thank the animal caretakers Sophie Brignon, Charley Pinault, and Dalila Ali-Hadji at PHENOMIN-ICS and the IGBMC animal facility. We wish to thank the Physiopathology of Aneuploidy, Gene Dosage Effect, and Down Syndrome Team for their help and in particular Michel

Roux, Maria-Victoria Hinckelman, Claire Chevalier, and Arnaud Duchon.

## SUPPLEMENTARY MATERIAL

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

## REFERENCES

Agarwal, N., Offermanns, S., and Kuner, R. (2004). Conditional Gene Deletion in Primary Nociceptive Neurons of Trigeminal Ganglia and Dorsal Root Ganglia. *Genesis* 38 (3), 122–129. doi:10.1002/gene.20010

Akopian, A. N., Souslova, V., England, S., Okuse, K., Ogata, N., Ure, J., et al. (1999). The Tetrodotoxin-Resistant Sodium Channel SNS Has a Specialized Function in Pain Pathways. *Nat. Neurosci.* 2 (6), 541–548. doi:10.1038/9195

Barrot, M. (2012). Tests and Models of Nociception and Pain in Rodents. *Neuroscience* 211, 39–50. doi:10.1016/j.neuroscience.2011.12.041

Bautista, D. M., Siemens, J., Glazer, J. M., Tsuruda, P. R., Basbaum, A. I., Stucky, C. L., et al. (2007). The Menthol Receptor TRPM8 Is the Principal Detector of Environmental Cold. *Nature* 448 (7150), 204–208. doi:10.1038/nature05910

Bennett, D. L., Clark, A. J., Huang, J., Waxman, S. G., and Dib-Hajj, S. D. (2019). The Role of Voltage-Gated Sodium Channels in Pain Signaling. *Physiol. Rev.* 99 (2), 1079–1151. doi:10.1152/physrev.00052.2017

Blasius, A. L., Dubin, A. E., Petrus, M. J., Lim, B. K., Narezkina, A., Criado, J. R., et al. (2011). Hypermorphogenic Mutation of the Voltage-Gated Sodium Channel Encoding Gene Scn10a Causes a Dramatic Stimulus-dependent Neurobehavioral Phenotype. *Proc. Natl. Acad. Sci. U S A.* 108 (48), 19413–19418. doi:10.1073/pnas.1117020108

Bohic, M., Marics, I., Santos, C., Malapert, P., Ben-Arie, N., Salio, C., et al. (2020). Loss of Bhla9 Impairs Thermotaxis and Formalin-Evoked Pain in a Sexually Dimorphic Manner. *Cell Rep* 30 (3), 602–e6. doi:10.1016/j.celrep.2019.12.041

Bouhassira, D. (2019). Neuropathic Pain: Definition, Assessment and Epidemiology. *Rev. Neurol. (Paris)* 175 (1-2), 16–25. doi:10.1016/j.neurol.2018.09.016

Brouwer, B. A., de Greef, B. T., Hoeijmakers, J. G., Geerts, M., van Kleef, M., Merkies, I. S., et al. (2015). Neuropathic Pain Due to Small Fiber Neuropathy in Aging: Current Management and Future Prospects. *Drugs Aging* 32 (8), 611–621. doi:10.1007/s40266-015-0283-8

Carter, R., Morton, A., and Dunnett, S. (2001). Motor Coordination and Balance in Rodents. *Curr. Protoc. Neurosci. Chapter 8, 8–12*. doi:10.1002/0471142301.ns0812s15

Colloca, L., Ludman, T., Bouhassira, D., Baron, R., Dickenson, A. H., Yarnitsky, D., et al. (2017). Neuropathic Pain. *Nat. Rev. Dis. Primers* 3, 17002. doi:10.1038/nrdp.2017.2

Dabby, R., Sadeh, M., Broitman, Y., Yosovich, K., Dickman, R., and Leshinsky-Silver, E. (2016). Painful Small Fiber Neuropathy with Gastroparesis: A New Phenotype with a Novel Mutation in the SCN10A Gene. *J. Clin. Neurosci.* 26, 84–88. doi:10.1016/j.jocn.2015.05.071

de Greef, B. T. A., Hoeijmakers, J. G. J., Gorissen-Brouwers, C. M. L., Geerts, M., Faber, C. G., and Merkies, I. S. J. (2018). Associated Conditions in Small Fiber Neuropathy - a Large Cohort Study and Review of the Literature. *Eur. J. Neurol.* 25 (2), 348–355. doi:10.1111/ene.13508

Eijkenboom, I., Sopacua, M., Hoeijmakers, J. G. J., de Greef, B. T. A., Lindsey, P., Almomani, R., et al. (2019). Yield of Peripheral Sodium Channels Gene Screening in Pure Small Fibre Neuropathy. *J. Neurol. Neurosurg. Psychiatry* 90 (3), 342–352. doi:10.1136/jnnp-2018-319042

Faber, C. G., Lauria, G., Merkies, I. S., Cheng, X., Han, C., Ahn, H. S., et al. (2012)19449. Gain-of-function Nav1.8 Mutations in Painful Neuropathy. *Proc. Natl. Acad. Sci. U S A.* 109 (47), 19444–19449. doi:10.1073/pnas.1216080109

Farrar, J. T., and Portenoy, R. K. (2001). Neuropathic Cancer Pain: the Role of Adjuvant Analgesics. *Oncology (Williston Park)* 15 (11), 1435–1450.

Garrison, S. R., Weyer, A. D., Barabas, M. E., Beutler, B. A., and Stucky, C. L. (2014). A Gain-Of-Function Voltage-Gated Sodium Channel 1.8 Mutation Drives Intense Hyperexcitability of A- and C-Fiber Neurons. *Pain* 155 (5), 896–905. doi:10.1016/j.pain.2014.01.012

Gijsbers, K., and Nicholson, F. (2005). Experimental Pain Thresholds Influenced by Sex of Experimenter. *Percept Mot. Skills* 101 (3), 803–807. doi:10.2466/pms.101.3.803-807

Han, C., Vasylyev, D., Macala, L. J., Gerrits, M. M., Hoeijmakers, J. G., Bekelaar, K. J., et al. (2014). The G1662S NaV1.8 Mutation in Small Fibre Neuropathy: Impaired Inactivation Underlying DRG Neuron Hyperexcitability. *J. Neurol. Neurosurg. Psychiatry* 85 (5), 499–505. doi:10.1136/jnnp-2013-306095

Huang, J., Yang, Y., Zhao, P., Gerrits, M. M., Hoeijmakers, J. G., Bekelaar, K., et al. (2013). Small-fiber Neuropathy Nav1.8 Mutation Shifts Activation to Hyperpolarized Potentials and Increases Excitability of Dorsal Root Ganglion Neurons. *J. Neurosci.* 33 (35), 14087–14097. doi:10.1523/JNEUROSCI.2710-13.2013

Huang, T., Lin, S. H., Malewicz, N. M., Zhang, Y., Zhang, Y., Goulding, M., et al. (2019). Identifying the Pathways Required for Coping Behaviours Associated with Sustained Pain. *Nature* 565 (7737), 86–90. doi:10.1038/s41586-018-0793-8

Huang, W., Liu, M., Yan, S. F., and Yan, N. (2017). Structure-based Assessment of Disease-Related Mutations in Human Voltage-Gated Sodium Channels. *Protein Cell* 8 (6), 401–438. doi:10.1007/s13238-017-0372-z

Källai, I., Barke, A., and Voss, U. (2004). The Effects of Experimenter Characteristics on Pain Reports in Women and Men. *Pain* 112 (1-2), 142–147. doi:10.1016/j.pain.2004.08.008

Kerr, B. J., Souslova, V., McMahon, S. B., and Wood, J. N. (2001). A Role for the TTX-Resistant Sodium Channel Nav 1.8 in NGF-Induced Hyperalgesia, but Not Neuropathic Pain. *Neuroreport* 12 (14), 3077–3080. doi:10.1097/00001756-200110080-00019

Laird, J. M., Souslova, V., Wood, J. N., and Cervero, F. (2002). Deficits in Visceral Pain and Referred Hyperalgesia in Nav1.8 (SNS/PN3)-null Mice. *J. Neurosci.* 22 (19), 8352–8356. doi:10.1523/jneurosci.22-19-08352.2002

Lauria, G., Cornblath, D. R., Johansson, O., McArthur, J. C., Mellgren, S. I., Nolano, M., et al. (2005). EFNS Guidelines on the Use of Skin Biopsy in the Diagnosis of Peripheral Neuropathy. *Eur. J. Neurol.* 12 (10), 747–758. doi:10.1111/j.1468-1331.2005.01260.x

Leo, S., D'Hooze, R., and Meert, T. (2010). Exploring the Role of Nociceptor-specific Sodium Channels in Pain Transmission Using Nav1.8 and Nav1.9 Knockout Mice. *Behav. Brain Res.* 208 (1), 149–157. doi:10.1016/j.bbr.2009.11.023

Luiz, A. P., MacDonald, D. I., Santana-Varela, S., Millet, Q., Sikandar, S., Wood, J. N., et al. (2019), Cold Sensing by NaV1.8-positive and NaV1.8-negative Sensory neurons. *Proc Natl Acad Sci U S A*, 116, 3811–3816.

Minett, M. S., Eijkelkamp, N., and Wood, J. N. (2014). Significant Determinants of Mouse Pain Behaviour. *PLoS One* 9 (8), e104458. doi:10.1371/journal.pone.0104458

Minett, M. S., Falk, S., Santana-Varela, S., Bogdanov, Y. D., Nassar, M. A., Heegaard, A. M., et al. (2014). Pain without Nociceptors? Nav1.7-independent Pain Mechanisms. *Cel Rep* 6 (2), 301–312. doi:10.1016/j.celrep.2013.12.033

Minett, M. S., Nassar, M. A., Clark, A. K., Passmore, G., Dickenson, A. H., Wang, F., et al. (2012). Distinct Nav1.7-dependent Pain Sensations Require Different Sets of Sensory and Sympathetic Neurons. *Nat. Commun.* 3, 791. doi:10.1038/ncomms1795
Mogil, J. S. (2020). Qualitative Sex Differences in Pain Processing: Emerging Evidence of a Biased Literature. *Nat. Rev. Neurosci.* 21 (7), 353–365. doi:10.1038/s41583-020-0310-6

Nassar, M. A., Levato, A., Stirling, L. C., and Wood, J. N. (2005). Neuropathic Pain Develops Normally in Mice Lacking Both Na(v)1.7 and Na(v)1.8. *Mol. Pain* 1, 24. doi:10.1186/1744-8069-1-24

Nassar, M. A., Stirling, L. C., Forlani, G., Baker, M. D., Matthews, E. A., Dickenson, A. H., et al. (2004). Nociceptor-specific Gene Deletion Reveals a Major Role for Nav1.7 (PN1) in Acute and Inflammatory Pain. *Proc. Natl. Acad. Sci. U S A.* 101 (34), 12706–12711. doi:10.1073/pnas.0404915101

Nasser, S. A., and Afify, E. A. (2019). Sex Differences in Pain and Opioid Mediated Antinociception: Modulatory Role of Gonadal Hormones. *Life Sci.* 237, 116926. doi:10.1016/j.lfs.2019.116926

Niesters, M., Dahan, A., Kest, B., Zacny, J., Stijnen, T., Aarts, L., et al. (2010). Do sex Differences Exist in Opioid Analgesia? A Systematic Review and Meta-Analysis of Human Experimental and Clinical Studies. *Pain* 151 (1), 61–68. doi:10.1016/j.pain.2010.06.012

Reker, A. N., Chen, S., Etter, K., Burger, T., Caudill, M., and Davidson, S. (2020). The Operant Plantar Thermal Assay: A Novel Device for Assessing Thermal Pain Tolerance in Mice. *eNeuro* 7 (2). doi:10.1523/ENEURO.0210-19.2020

Rosen, S., Ham, B., and Mogil, J. S. (2017). Sex Differences in Neuroimmunity and Pain. *J. Neurosci. Res.* 95 (1-2), 500–508. doi:10.1002/jnr.23831

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

Shields, S. D., Cheng, X., Uçeyler, N., Sommer, C., Dib-Hajj, S. D., and Waxman, S. G. (2012). Sodium Channel Na(v)1.7 Is Essential for Lowering Heat Pain Threshold after Burn Injury. *J. Neurosci.* 32 (32), 10819–10832. doi:10.1523/JNEUROSCI.0304-12.2012

Sorge, R. E., and Totsch, S. K. (2017). Sex Differences in Pain. *J. Neurosci. Res.* 95 (6), 1271–1281. doi:10.1002/jnr.23841

Stirling, L. C., Forlani, G., Baker, M. D., Wood, J. N., Matthews, E. A., Dickenson, A. H., et al. (2005). Nociceptor-specific Gene Deletion Using Heterozygous Nav1.8-Cre Recombinase Mice. *Pain* 113 (1-2), 27–36. doi:10.1016/j.pain.2004.08.015

Tappe-Theodor, A., King, T., and Morgan, M. M. (2019). Pros and Cons of Clinically Relevant Methods to Assess Pain in Rodents. *Neurosci. Biobehav Rev.* 100, 335–343. doi:10.1016/j.neubiorev.2019.03.009

Treede, R. D., Jensen, T. S., Campbell, J. N., Cruccu, G., Dostrovsky, J. O., Griffin, J. W., et al. (2008). Neuropathic Pain: Redefinition and a Grading System for Clinical and Research Purposes. *Neurology* 70 (18), 1630–1635. doi:10.1212/01.wnl.0000282763.29778.59

Weibel, R., Reiss, D., Karchewski, L., Gardon, O., Matifas, A., Filliol, D., et al. (2013). Mu Opioid Receptors on Primary Afferent nav1.8 Neurons Contribute to Opiate-Induced Analgesia: Insight from Conditional Knockout Mice. *PLoS One* 8 (9), e74706. doi:10.1371/journal.pone.0074706

Woolf, C. J. (2010). What Is This Thing Called Pain? *J. Clin. Invest.* 120 (11), 3742–3744. doi:10.1172/JCI45178

Xiao, Y., Barbosa, C., Pei, Z., Xie, W., Strong, J. A., Zhang, J. M., et al. (2019). Increased Resurgent Sodium Currents in Nav1.8 Contribute to Nociceptive Sensory Neuron Hyperexcitability Associated with Peripheral Neuropathies. *J. Neurosci.* 39 (8), 1539–1550. doi:10.1523/JNEUROSCI.0468-18.2018

Xue, Y., Chidiac, C., Herault, Y., and Gaveriaux-Ruff, C. (2021). Pain Behavior in SCN9A (Nav1.7) and SCN10A (Nav1.8) Mutant Rodent Models. *Neurosci. Lett.* 753, 135844. doi:10.1016/j.neulet.2021.135844

Zimmermann, K., Leffler, A., Babes, A., Cendan, C. M., Carr, R. W., Kobayashi, J., et al. (2007). Sensory Neuron Sodium Channel Nav1.8 Is Essential for Pain at Low Temperatures. *Nature* 447 (7146), 855–858. doi:10.1038/nature05880

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

**Publisher's Note:** All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Copyright © 2021 Chidiac, Xue, Muniz Moreno, Bakr Rasheed, Lorentz, Birling, Gaveriaux-Ruff and Herault. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.