Thibaud Parpaite

Lucie Brosse

Nina Séjourné

Amandine Laur

Yasmine Mechioukhi

Patrick Delmas

Bertrand Coste
email: bertrand.coste@univ-amu.fr

Patch-seq of mouse DRG neurons reveals candidate genes for specific mechanosensory functions

Keywords: Somatosensation, mechanotransduction, nociception, ion channel, pain

à 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.

INTRODUCTION

Mechanosensory neurons are involved in multiple sensory modalities, including innocuous tickle, pleasant and discriminative touch, proprioception and kinesthesia as well as various mechanical pain-related sensations such as sharp or deep aching pain and visceral pain. A rich variety of dedicated dorsal root ganglion (DRG) mechanosensory neurons and structured endings innervate our skin, mucosa, muscles, tendons and joints to initiate these mechanosensory sub-modalities [START_REF] Abraira | The sensory neurons of touch[END_REF][START_REF] Proske | The proprioceptive senses: their roles in signaling body shape, body position and movement, and muscle force[END_REF][START_REF] Roudaut | Touch sense: functional organization and molecular determinants of mechanosensitive receptors[END_REF]. Physical cues are detected at the sensory nerve endings/auxiliary cells, where specialized mechanosensitive (MS) ion channels convert mechanical stimuli into electrochemical signals [START_REF] Basbaum | Cellular and Molecular Mechanisms of Pain[END_REF][START_REF] Kefauver | Discoveries in structure and physiology of mechanically activated ion channels[END_REF].

The discovery from low-throughput patch-clamp screening of PIEZO proteins as genuine forcesensing channels has provided better understanding of mechanical force transduction mechanisms in somatosensory neurons [START_REF] Coste | Piezo1 and Piezo2 are essential components of distinct mechanically activated cation channels[END_REF] as well as in multiple biological systems [START_REF] Murthy | Piezos thrive under pressure: mechanically activated ion channels in health and disease[END_REF][START_REF] Wu | Touch, Tension, and Transduction -The Function and Regulation of Piezo Ion Channels[END_REF]. PIEZO1 and PIEZO2 form a unique class of a conserved family of nonselective cation channels that are activated by physical stimuli, including pressure, stretch and shear stress [START_REF] Coste | Piezo proteins are pore-forming subunits of mechanically activated channels[END_REF]Ranade et al., 2014a). Both channels generate MS currents with rapid inactivation that convert force into cellular responses on a millisecond time scale. While PIEZO1 is broadly expressed in non-neuronal tissues exposed to pressure and fluid flow, e.g., kidneys, bladder, endothelial cells, and blood cells, PIEZO2 is predominantly expressed in sensory tissues, such as DRG, trigeminal ganglia and Merkel cells [START_REF] Murthy | Piezos thrive under pressure: mechanically activated ion channels in health and disease[END_REF][START_REF] Wu | Touch, Tension, and Transduction -The Function and Regulation of Piezo Ion Channels[END_REF]. PIEZOs are co-expressed in nodose and petrosal sensory ganglia where they are shown to be involved baroreceptor reflex [START_REF] Min | Arterial Baroreceptors Sense Blood Pressure through Decorated Aortic Claws[END_REF][START_REF] Zeng | PIEZOs mediate neuronal sensing of blood pressure and the baroreceptor reflex[END_REF]. Studies from knock-out (KO) mouse models have shown that Piezo2 plays a crucial role in innocuous touch sensation [START_REF] Ranade | Piezo2 is the major transducer of mechanical forces for touch sensation in mice[END_REF] and proprioception [START_REF] Florez-Paz | A critical role for Piezo2 channels in the mechanotransduction of mouse proprioceptive neurons[END_REF][START_REF] Woo | Piezo2 is the principal mechanotransduction channel for proprioception[END_REF]. Piezo2 also plays a more marginal role in mechanical nociception [START_REF] Murthy | The mechanosensitive ion channel Piezo2 mediates sensitivity to mechanical pain in mice[END_REF][START_REF] Szczot | PIEZO2 mediates injury-induced tactile pain in mice and humans[END_REF]. These functions appeared to be conserved in humans, as patients with loss-of-function variants in PIEZO2 show deficits in touch discrimination and joint proprioception, and failed to develop mechanical allodynia after skin inflammation [START_REF] Chesler | The Role of PIEZO2 in Human Mechanosensation[END_REF][START_REF] Szczot | PIEZO2 mediates injury-induced tactile pain in mice and humans[END_REF].

DRG neurons express a large repertoire of excitatory MS currents classified as rapidly-adapting (RA), intermediately-adapting (IA), and slowly/ultra slowly-adapting (SA/ultra SA) currents, according to their inactivation kinetics to sustained mechanical stimulation [START_REF] Coste | Piezo1 and Piezo2 are essential components of distinct mechanically activated cation channels[END_REF][START_REF] Drew | Acid-sensing ion channels ASIC2 and ASIC3 do not contribute to mechanically activated currents in mammalian sensory neurones[END_REF][START_REF] Drew | Distinct mechanosensitive properties of capsaicinsensitive and -insensitive sensory neurons[END_REF][START_REF] Francois | The Low-Threshold Calcium Channel Cav3.2 Determines Low-Threshold Mechanoreceptor Function[END_REF][START_REF] Hao | Multiple desensitization mechanisms of mechanotransducer channels shape firing of mechanosensory neurons[END_REF][START_REF] Hu | Mechanosensitive currents in the neurites of cultured mouse sensory neurones[END_REF][START_REF] Ranade | Piezo2 is the major transducer of mechanical forces for touch sensation in mice[END_REF][START_REF] Rugiero | Kinetic properties of mechanically activated currents in spinal sensory neurons[END_REF]. It is now well established that PIEZO2 channels sustain RA MS currents [START_REF] Coste | Piezo1 and Piezo2 are essential components of distinct mechanically activated cation channels[END_REF][START_REF] Ranade | Piezo2 is the major transducer of mechanical forces for touch sensation in mice[END_REF]. Importantly, IA and slow-type MS currents are largely unaffected in Piezo2 KO DRG neurons [START_REF] Ranade | Piezo2 is the major transducer of mechanical forces for touch sensation in mice[END_REF], demonstrating that as-yet-unknown MS channels must account for the other MS currents. It is worthwhile to identify these channels since many sensations related to innocuous and noxious mechanical stimuli are independent from PIEZO2 in humans [START_REF] Case | Innocuous pressure sensation requires A-type afferents but not functional RhoIotaEpsilonZetaOmicron2 channels in humans[END_REF][START_REF] Chesler | The Role of PIEZO2 in Human Mechanosensation[END_REF][START_REF] Szczot | PIEZO2 mediates injury-induced tactile pain in mice and humans[END_REF]. Recently, two other genes, namely Tmem120a (Tacan) and Tmem150c (Tentonin 3), have been proposed to encode ion channels sustaining slow MS currents in DRG neurons [START_REF] Beaulieu-Laroche | TACAN Is an Ion Channel Involved in Sensing Mechanical Pain[END_REF][START_REF] Hong | Tentonin 3/TMEM150c Confers Distinct Mechanosensitive Currents in Dorsal-Root Ganglion Neurons with Proprioceptive Function[END_REF], but some of these findings are controversial [START_REF] Anderson | TMEM150C/Tentonin3 Is a Regulator of Mechano-gated Ion Channels[END_REF][START_REF] Dubin | Endogenous Piezo1 Can Confound Mechanically Activated Channel Identification and Characterization[END_REF]. Hence, despite significant advances, the full panel of MS channels involved in mechanosensation has yet to be established.

Identification of MS channels is hampered by the large heterogeneity of DRG and trigeminal sensory neurons. This heterogeneity has been well documented by large scale single-cell transcriptomics enabling gene expression-based classification [START_REF] Li | Somatosensory neuron types identified by high-coverage single-cell RNA-sequencing and functional heterogeneity[END_REF][START_REF] Nguyen | Stereotyped transcriptomic transformation of somatosensory neurons in response to injury[END_REF][START_REF] Renthal | Transcriptional Reprogramming of Distinct Peripheral Sensory Neuron Subtypes after Axonal Injury[END_REF][START_REF] Sharma | The emergence of transcriptional identity in somatosensory neurons[END_REF][START_REF] Usoskin | Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencing[END_REF][START_REF] Zeisel | Molecular Architecture of the Mouse Nervous System[END_REF]. Eight main classes of DRG sensory neurons have been identified, including low-threshold mechanoreceptors, proprioceptors, thermoreceptors, nociceptors, itch sensitive neurons and type C low-threshold mechanoreceptors, which can be further subdivided into up to 17 clusters identified genetically [START_REF] Zeisel | Molecular Architecture of the Mouse Nervous System[END_REF]. A limitation of large-scale classification based on unsupervised grouping of neurons with similar expression profiles is the lack of direct functional correlates. Indeed, this would require a prior functional characterization of individual neurons by patch-clamp experiments, a laborious step incompatible with large samples of neurons. Nevertheless, the direct correlation of single cell transcriptomes with functional mechanosensory properties could serve as a basis for identifying the genes involved in mechanotransduction.

Here, we used Patch-seq methodology [START_REF] Cadwell | Electrophysiological, transcriptomic and morphologic profiling of single neurons using Patch-seq[END_REF][START_REF] Fuzik | Integration of electrophysiological recordings with single-cell RNA-seq data identifies neuronal subtypes[END_REF] to combine recordings of MS currents in cultured mouse DRG neurons with single cell RNA sequencing. In silico analysis of collected data with a previous large-scale dataset allowed to link MS current subtypes to molecular neuronal clusters. These data made it possible to precise the distribution of the four different types of MS currents in molecularly defined populations of DRG neurons. Moreover, differential comparison provided a list of candidate genes for mechanotransduction complexes. This dataset combining patch-clamp electrophysiology and quantitative single cell RNA sequencing constitutes an open-resource to explore further the cell-type-specific determinants of mechanosensory properties.

RESULTS

Coupling patch-clamp recording of MS currents with single cell RNA sequencing

To identify molecular features associated with specific MS phenotypes, we combined single cell RNAsequencing (scRNAseq) profiling with electrophysiological characterization of MS currents in individual DRG neurons. A representation of the methodological procedure used for data collection is depicted in figure 1A. Briefly, sensory neurons from mouse DRG cultures were recorded using the whole-cell patch-clamp mode in order to characterize the properties of excitatory MS currents in response to cell poking with a mechanical probe [START_REF] Coste | Pharmacological Dissection and Distribution of NaN/Nav1.9, T-type Ca2+ Currents, and Mechanically Activated Cation Currents in Different Populations of DRG Neurons[END_REF][START_REF] Coste | Piezo1 and Piezo2 are essential components of distinct mechanically activated cation channels[END_REF][START_REF] Drew | Distinct mechanosensitive properties of capsaicinsensitive and -insensitive sensory neurons[END_REF][START_REF] Hao | Multiple desensitization mechanisms of mechanotransducer channels shape firing of mechanosensory neurons[END_REF][START_REF] Hu | Mechanosensitive currents in the neurites of cultured mouse sensory neurones[END_REF][START_REF] Mccarter | Mechanical transduction by rat dorsal root ganglion neurons in vitro[END_REF][START_REF] Michel | Maturational Changes in Mouse Cutaneous Touch and Piezo2-Mediated Mechanotransduction[END_REF][START_REF] Rugiero | Kinetic properties of mechanically activated currents in spinal sensory neurons[END_REF].

After functional characterization, the content of each cell was harvested through the recording pipette prior to the generation of cDNA libraries and RNA sequencing (see M&M for detailed procedure). Using this approach, we selected 62 DRG neurons covering the main types of MS currents (Fig. 1B) for cDNA library preparation, of which 53 cells proved suitable for scRNAseq analysis. A total of 1.6 billion pairedend reads were generated, among which 1.4 billion were correctly aligned to mouse genome, giving an average of 25.9 million mapped reads for a given neuron. The mean number of detectable genes per cell was 13,359 ± 1,478. Analysis of reference genes showed no differences between samples, which attested the absence of batch-effect between the sequencing runs (Fig. S1A). Quality metrics values from sequenced neurons were within the same range than large scale DRG scRNAseq studies [START_REF] Li | Somatosensory neuron types identified by high-coverage single-cell RNA-sequencing and functional heterogeneity[END_REF][START_REF] Usoskin | Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencing[END_REF][START_REF] Zeisel | Molecular Architecture of the Mouse Nervous System[END_REF], attesting that patch-clamping/poking the neurons prior to RNA sequencing analysis did not impact gene detection (Fig. S1B). Overall, the depth of our molecular analysis allowed us to quantitatively assay the expression of genes in each DRG neuron and to compare expression differences of any of these genes in light with mechanical properties.

Categorization of MS current-expressing neurons prior to transcriptome profiling

To tentatively classify mechanosensory neurons based on their mechanical properties, we examined the kinetics of MS currents over a large population of mouse DRG neurons (n = 312), not intended for scRNAseq. Previous studies have identified 3 to 4 kinetically distinct MS currents in rat and mouse DRG neurons, which may reflect the activation of non-uniform populations of ion channels [START_REF] Coste | Piezo1 and Piezo2 are essential components of distinct mechanically activated cation channels[END_REF][START_REF] Drew | Acid-sensing ion channels ASIC2 and ASIC3 do not contribute to mechanically activated currents in mammalian sensory neurones[END_REF][START_REF] Drew | Distinct mechanosensitive properties of capsaicinsensitive and -insensitive sensory neurons[END_REF][START_REF] Francois | The Low-Threshold Calcium Channel Cav3.2 Determines Low-Threshold Mechanoreceptor Function[END_REF][START_REF] Hao | Multiple desensitization mechanisms of mechanotransducer channels shape firing of mechanosensory neurons[END_REF][START_REF] Hu | Mechanosensitive currents in the neurites of cultured mouse sensory neurones[END_REF][START_REF] Ranade | Piezo2 is the major transducer of mechanical forces for touch sensation in mice[END_REF][START_REF] Rugiero | Kinetic properties of mechanically activated currents in spinal sensory neurons[END_REF]. Based on exponential fit of inactivation kinetics, we identified four kinetically distinct MS currents with gaussian mean value centered at 5.9 ± 1. 6, 15.1 ± 4.5, 75.4 ± 13.8 and 200.9 ± 66.8 ms (mean ± SD) (Fig. 1B). These currents were referred hereafter as rapidly adapting (RA; τinac ≤ 10 ms), intermediately adapting (IA; 10 ms < τinac ≤ 30 ms), slowly adapting (SA; 30 ms < τinac ≤ 110 ms) and ultra-SA (τinac > 110 ms). Prototypical recordings of these currents are depicted in figure 1C. It is worth noting that a substantial part of mechanosensory DRG neurons exhibited more than one single current component, potentially reflecting the activation of two or more populations of distinct MS channels. For analytical purpose, we only considered mixed currents best fitted by bi-exponential functions. The contribution of each component to the whole MS current was estimated from the amplitude of the respective exponential fit (Fig. S2, see M&M for detailed procedure). Representative examples of biphasic currents classified as RA/SA, RA/ultra-SA and IA/ultra-SA are illustrated in figure 1D. This kinetics-based categorization made it possible to classify our successfully sequenced DRG neurons.

Thus, from our sample of 53 scRNAseq neurons, 49% displayed a monophasic MS current, while the remaining 45.5% showed a mixed current, with variable MS current components. Only 3 scRNAseq neurons were unresponsive to mechanical stimulus, but this does not reflect physiological proportion [START_REF] Coste | Piezo1 and Piezo2 are essential components of distinct mechanically activated cation channels[END_REF][START_REF] Hu | Mechanosensitive currents in the neurites of cultured mouse sensory neurones[END_REF]. Average values of τinac were 5.9 ± 0.3, 18.8 ± 1.4, 82.2 ± 2.8 and 347.7 ± 26.5 ms for RA, IA, SA and ultra-SA current, respectively (n = 26, 15, 10 and 24; mean ± sem) (Fig. 1E). The detailed profiles of MS currents in individual neurons are described in Table S1.

Mapping Patch-seq transcriptomes of MS neurons to single-cell RNA-seq datasets

Previous scRNAseq studies performed on large samples of neurons have used gene expression profile analysis to sort out DRG neurons into molecularly defined neuronal subclasses [START_REF] Li | Somatosensory neuron types identified by high-coverage single-cell RNA-sequencing and functional heterogeneity[END_REF][START_REF] Sharma | The emergence of transcriptional identity in somatosensory neurons[END_REF][START_REF] Usoskin | Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencing[END_REF][START_REF] Zeisel | Molecular Architecture of the Mouse Nervous System[END_REF]. To achieve high-quality classification of our neuron sample, we merged our data with the dataset of 1580 scRNAseq DRG neurons from Zeisel and collaborators, which is available as open source [START_REF] Zeisel | Molecular Architecture of the Mouse Nervous System[END_REF]. Then, we performed a joint analysis using mutual nearest neighbors (MNN) correction by focusing on the expression magnitude of the 211 most highly enriched genes across DRG neuron clusters [START_REF] Zeisel | Molecular Architecture of the Mouse Nervous System[END_REF] (Fig. S1C). This procedure allowed to classify each neuron from our sample in transcriptional clusters identified previously [START_REF] Zeisel | Molecular Architecture of the Mouse Nervous System[END_REF] (Table S1). The merged datasets were mapped using t-distributed stochastic neighbor embedding (tSNE) projection (Fig. 2A). Neurons from our sample spread across seven out of eight functionally distinct DRG neuron populations, some of which are composed of multiple transcriptional clusters whose functional specificity remains to be determined. Eight neurons were classified into the NF group of neurons typically considered as low-threshold mechanoreceptors (LTMR) and proprioceptors. Six neurons belong to NP1 population related to polymodal nociceptors, whereas eleven neurons were assigned to NP2 population and a single neuron to NP3 population, which are related to itch-specific sensory neurons (Cavanaugh et al., 2009;Dong and Dong, 2018;[START_REF] Han | A subpopulation of nociceptors specifically linked to itch[END_REF]Wang and Zylka, 2009). Nine neurons were assigned to PEP1, a population of C-type thermo-nociceptors, and sixteen neurons to PEP2, a population of lightly myelinated Aδ nociceptors [START_REF] Usoskin | Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencing[END_REF]. Two neurons segregated into TH neuron population, representing C-low threshold mechanoreceptors (C-LTMR), which contribute to pleasant touch [START_REF] Li | The functional organization of cutaneous low-threshold mechanosensory neurons[END_REF][START_REF] Olausson | The neurophysiology of unmyelinated tactile afferents[END_REF] and mechanical allodynia [START_REF] Seal | Injury-induced mechanical hypersensitivity requires C-low threshold mechanoreceptors[END_REF]. Finally, none of our scRNAseq neurons belong to TrpM8 population. These neurons represent a labeled line for cold sensation, and their genetic ablation in mice has no effect on mechanical behaviors [START_REF] Knowlton | A sensory-labeled line for cold: TRPM8-expressing sensory neurons define the cellular basis for cold, cold pain, and cooling-mediated analgesia[END_REF]. Expression of marker genes and genes coding for thermosensitive and pruriceptive molecular sensors confirmed the appropriate mapping of scRNAseq neurons from our sample (Fig. 2B-C). NF LTMR and proprioceptors express the specific markers Nefh (neurofilament heavy chain), Pvalb (parvalbumin), Calb1 (calbindin1), Ntrk2 and Ntrk3 (TrkB and C), but thermosensitive and pruriceptive itch molecular sensors were not detected. NP1-3 sensory neurons express P2rx3 (P2X3), PlnxC1 (plexin C1) and Necab2 (N-terminal EF-hand calcium binding protein 2). The itch related genes Lpar3 and Lpar5 (lysophosphatidic acid receptors 3 and 5), and Mrgprd (MGRPRD) are expressed in NP1 neurons, whereas NP2 population expresses Mrgpra3 and Mrgprx1 (MRGPRA3 and X1). PEP neurons express Calca (CGRP) and Kit (c-Kit) but differ in the expression of Gal (galanin prepropeptide) for PEP1 C-type thermo-nociceptors and Nefh for PEP2 lightly myelinated Aδ nociceptors. The heat noxious sensors Trpa1, Trpv1 and Trpm3 [START_REF] Vandewauw | A TRP channel trio mediates acute noxious heat sensing[END_REF] were strongly detected in NP1, NP2 and PEP1 Ctype neurons compared with PEP2 Aδ-type neurons. This is consistent with previous observation in skin-nerve recording studies that showed a higher proportion of C-type than Aδ-type fibers that respond to noxious heat [START_REF] Cain | Response properties of mechanoreceptors and nociceptors in mouse glabrous skin: an in vivo study[END_REF][START_REF] Koltzenburg | Receptive properties of mouse sensory neurons innervating hairy skin[END_REF]. Finally, the TH C-LTMR characterized by expression of Th (tyrosine hydroxylase), Fam19a4 (TAFA4) and Slc17a8 (Vglut3), are devoid of thermosensitive or itch-related molecular sensors.

Distribution of MS currents in genetically defined DRG subclasses

The characterization of MS current types expressed in scRNAseq neurons allowed us to examine their distribution in genetically defined DRG populations (Fig. 3A-B). Except for NF neurons, which exhibited almost exclusively the RA current, all other molecularly defined DRG subclasses displayed more than one type of MS current. However, the prevalence of each current component differed among DRG groups. The Figure 3C shows the proportion of each MS current type in genetically defined DRG neuron populations. The polymodal nociceptor NP1 subclass displayed all four types of MS currents, with individual contributions ranging from 10 to 33%. The itch specific NP2 subclass expressed almost no SA current and showed a higher amount of RA current (42%) compared with IA and ultra-SA currents (35 and 20%, respectively). The C-type thermo-nociceptor PEP1 subclass had prominent IA current (56%), a lower amount of RA current (33%), with no or little SA and ultra-SA current components, respectively. The Aδ nociceptor PEP2 subclass almost lacked IA current, and expressed prominent ultra-SA current (39%) and robust RA and SA currents (34 and 23%, respectively). Finally, C-LTMR TH neurons showed both RA and ultra-SA currents. Although our small sample size precluded definite conclusion about the nature of MS current in C-LTMR TH neurons, our data are consistent with previous work indicating that TH neurons mainly express RA and ultra-SA currents, or a mixture of both [START_REF] Delfini | TAFA4, a chemokine-like protein, modulates injury-induced mechanical and chemical pain hypersensitivity in mice[END_REF]. Altogether, these results revealed that the genetically defined DRG subclasses have different complements of MS currents, which may reflect some heterogeneity and functional differences within a given DRG subclass.

Expression analysis of genes hypothesized to encode MS channels

It has been remarkably difficult to identify candidate excitatory MS channels in mammals. With the notable exception of PIEZOs and of TMCs in auditory hair cells [START_REF] Corey | Function and Dysfunction of TMC Channels in Inner Ear Hair Cells[END_REF], other putative candidates have not been confirmed or might be modulators of mechanosensation rather than being directly involved in force transduction. A prerequisite for a candidate MS channel gene is to be found expressed in neurons displaying the purposed MS current. Therefore, we examined whether the expression of genes that have been proposed to encode excitatory MS channels in mammals correlates with the presence of a given MS current subtype across the scRNAseq neurons (Fig. 4A). As expected, Piezo2 was detected in all neurons exhibiting the RA current. However, Piezo2 was also detected in neurons for which no RA current could be extracted. Tmem120a, which has been proposed to sustain MS currents with slow kinetics in DRG neurons [START_REF] Beaulieu-Laroche | TACAN Is an Ion Channel Involved in Sensing Mechanical Pain[END_REF] was detected at high level in all neuron types, including neurons lacking MS current. Tmem150c (Tentonin 3), which has been proposed to encode MS channels with slow kinetics in DRG and nodose ganglion neurons [START_REF] Hong | Tentonin 3/TMEM150c Confers Distinct Mechanosensitive Currents in Dorsal-Root Ganglion Neurons with Proprioceptive Function[END_REF][START_REF] Lu | Tentonin 3/TMEM150C senses blood pressure changes in the aortic arch[END_REF], was detected in most neurons, regardless of their mechanical signature.

Tmem87a, a gene involved in melanoma adhesion/migration and sufficient to reconstitute MS currents in heterologous systems [START_REF] Patkunarajah | TMEM87a/Elkin1, a component of a novel mechanoelectrical transduction pathway, modulates melanoma adhesion and migration[END_REF], was detected in some neurons, irrespective of the MS current subtypes. All three Tmem63 genes which belong to an evolutionary conserved family of MS channel coding genes across eukaryotes (Murthy et al., 2018a) were differentially expressed.

Tmem63a was detected in few neurons but Tmem63b and Tmem63c were found in most MS currentexpressing neurons. Finally, Piezo1 was sporadically detected, while Tmc1 and Tmc2, which contribute to the pore of MS channels in auditory and vestibular hair cells [START_REF] Jia | TMC1 and TMC2 Proteins Are Pore-Forming Subunits of Mechanosensitive Ion Channels[END_REF][START_REF] Pan | TMC1 Forms the Pore of Mechanosensory Transduction Channels in Vertebrate Inner Ear Hair Cells[END_REF][START_REF] Pan | TMC1 and TMC2 are components of the mechanotransduction channel in hair cells of the mammalian inner ear[END_REF] were almost not expressed. Together, these results lend credence to the idea that PIEZO1 and TMC1 and 2 are unlikely to contribute to MS currents in DRG neurons. They also indicate that the other candidate genes are expressed in different MS current neuron families with no apparent clustering/specific expression.

Enrichment of transcripts as a mean to identify MS channel gene candidates

The above data failed to reveal specific expression of a given MS channel gene in a particular neuronal group. Therefore, we assessed whether some transcripts showed statistically over-representation (e.g. enrichment) in neurons clustered by current types. We determined the differential expression of a gene by comparing its expression in neurons displaying a specific current type against the neurons that do not express this current type. The plots illustrated in figure 4B provide graphical views of the enrichment of each MS channel gene in neurons with RA, IA, SA and ultra-SA components. Piezo2 is differentially expressed (P = 0.015) in neurons expressing RA currents, exhibiting a 2.3-fold enrichment. This result indicates that relative enrichment rather than specific expression may be a pertinent criterion when selecting putative MS current genes. Remarkably, none of the other MS channel gene candidates were found to be differentially expressed (P > 0.05) in IA, SA or ultra-SA MS currentexpressing neurons. This lack of enrichment raises doubt about their possible contribution to specific MS currents in DRG neurons. As this result is intriguing regarding the recent literature, we set up proofof-principle experiments to probe the functional relevance of Tmem120a and Tmem150c in MS currents.

Inhibition of Tmem120a expression does not alter DRG neuron MS currents

Our scRNAseq data indicated that Tmem120a is broadly expressed in all neuron subclasses. We further characterized Tmem120a expression in DRG slices using the RNAscope in situ hybridization (ISH) technology coupled with immunostaining for the broader nociceptive and non-nociceptive markers, Peripherin and NF200, respectively [START_REF] Goldstein | NF-L and peripherin immunoreactivities define distinct classes of rat sensory ganglion cells[END_REF][START_REF] Padilla | Expression and localization of the Nav1.9 sodium channel in enteric neurons and in trigeminal sensory endings: Implication for intestinal reflex function and orofacial pain[END_REF] (Fig. 5A-C, Fig. S3).

Tmem120a transcripts were detected in 75.0 % of NF200-positive neurons and 67.4 % of peripherinpositive DRG neurons (Fig. 5B). Moreover, analysis of the cross-sectional area distribution revealed that Tmem120a mRNA was detected in the majority of DRG neurons over a large range of cell soma diameters (Fig. 5C). Altogether, these results confirm that Tmem120a is widely expressed in the different subpopulations of sensory neurons.

We performed siRNA experiments in DRG neurons to determine whether Tmem120a could be responsible for a MS current. Since DRG neuron electroporation typically displays low transfection efficiency (≈10%), we first assessed Tmem120a siRNA efficacy by RT-qPCR in NIH/3T3 electroporated cells (transfection efficiency 60-80%), which constitutively express Tmem120a (Fig. S4A). In these cells, Tmem120a siRNA transfection induced a decrease of 72 ± 5 % of Tmem120a mRNA level compared to non-targeting siRNA transfection (Fig. S4D). Interestingly, knock-down of Piezo1 but not Tmem120a inhibited the MS current present in NIH/3T3 cells (Fig. S4F). The efficacy of siRNAs was confirmed in manually sorted electroporated DRG neurons by qPCR showing about 90% decrease of transcript levels (Fig. S4G).

Co-electroporation of siRNAs and GFP-plasmid was then performed in mouse DRG neurons, thereby allowing the fluorescent detection of transfected neurons. Non-targeting and Piezo2 siRNAs were used as negative and positive controls, respectively. We first compared the proportions of DRG neurons classified according to their main MS current component. As previously reported [START_REF] Coste | Piezo1 and Piezo2 are essential components of distinct mechanically activated cation channels[END_REF], inhibition of Piezo2 expression significantly (p < 0.0001, χ 2 test) altered the relative proportion of DRG neurons exhibiting the different MS currents, with notable decrease in the incidence of RA currentexpressing neurons from 38.1 to 10.0% and consequential increase of non-responsive neurons from 8.2 to 22.9%. By contrast, inhibition of Tmem120a expression had no effect on the relative proportion of MS currents compared to control conditions (p = 0.9, χ 2 test) (Fig. 5D). We next analyzed each MS current component quantified from mono-or bi-exponential fits of MS current recordings in siRNAtreated DRG neurons (Fig. 5E). Piezo2 siRNA treatment significantly reduced the proportion of neurons expressing the RA current component as well as the mean amplitude of the residual RA current. Once again, there were no significant changes in the proportion of any MS current components or in the mean amplitude of these components when Tmem120a expression was inhibited (Fig. 5E).

To restrain the neuronal population investigated, we took advantage of the observation that TMEM120A protein immunostaining was detected in 80% of IB4 + DRG neurons [START_REF] Beaulieu-Laroche | TACAN Is an Ion Channel Involved in Sensing Mechanical Pain[END_REF]. Therefore, we performed another set of experiments by selectively recording from IB4 + neurons identified using IB4-Alexa Fluor 568 conjugate staining. Consistent with our previous results, we found no differences regarding either the incidence or the amplitude of the different MS current components between Tmem120a and non-targeting siRNA-treated neurons (Fig. 5F-G). Furthermore, we cloned Tmem120a cDNA from mouse DRG into an expression vector. Cell-poking and pressureclamp experiments showed that Tmem120a overexpression did not induce MS currents in HEK-P1KO cells in which Piezo1 has been deleted [START_REF] Dubin | Endogenous Piezo1 Can Confound Mechanically Activated Channel Identification and Characterization[END_REF] (Fig. S5). Altogether, these data suggest that Tmem120a does not contribute to any MS currents in DRG neurons.

Inhibition of Tmem150c expression does not impair MS currents in DRG neurons

Tmem150c has been proposed to encode SA MS channels in proprioceptive DRG neurons contributing to motor coordination [START_REF] Hong | Tentonin 3/TMEM150c Confers Distinct Mechanosensitive Currents in Dorsal-Root Ganglion Neurons with Proprioceptive Function[END_REF]. This conclusion has been challenged by independent groups, because Tmem150c expression by itself does not appear to generate MS current but rather modulates the activity of mechano-gated channels including PIEZO1, PIEZO2 and the two-pore domain K + channel TREK1 [START_REF] Anderson | TMEM150C/Tentonin3 Is a Regulator of Mechano-gated Ion Channels[END_REF][START_REF] Dubin | Endogenous Piezo1 Can Confound Mechanically Activated Channel Identification and Characterization[END_REF]. Our scRNAseq data indicated that Tmem150c is not enriched in a particular subset of DRG neurons. Therefore, to gain more understanding about its role in MS current generation we recorded from Tmem150c siRNA electroporated DRG neurons. We found no differences whatsoever regarding either the incidence or the amplitude of the different MS current components in Tmem150c siRNA-treated neurons compared with non-targeting siRNA-treated neurons (Fig. 5D-E). These data argue against a role of Tmem150c in sustaining a particular type of MS current in DRG neurons.

Genes associated with DRG MS current types

Enrichment of Piezo2 in the subset of RA current expressing neurons provides a rationale to identify candidate genes coding for or regulating MS channels. Therefore, we looked for transcripts differentially expressed (P < 0.05) and enriched at least 2-fold in neurons grouped by MS current type (Fig. 6A-D). We provide a list of transcripts associated with RA, IA, SA and ultra-SA currents containing 212, 404, 125 and 84 protein coding genes, respectively (Table S2). As can be seen in figure 6E, Piezo2 emerged in the top 15% of most expressed genes amongst the enriched ones. Therefore, we reasoned that both enrichment and mRNA abundance in the corresponding neuronal population may be good criteria to pinpoint molecular determinants associated with MS current phenotype. Additionally, gene ontology (GO) analysis was used to annotate these genes with relevant knowledge such as detection of a mechanical stimulus, ion channel activity and integral component of the membrane. The top 40 of the most expressed genes associated with MS current type are listed in figure 6E. Among those, Piezo2 is the only gene known to encode a protein with demonstrated function in mechanical transduction. Moreover, the few genes known to encode ion channels, such as Scn10a encoding the voltage activated sodium channel Nav1.8 [START_REF] Han | Sodium channel Nav1.8: Emerging links to human disease[END_REF], are well characterized and unlikely to encode MS channels. Several genes are known to code for proteins inserted in the membrane, a feature common to ion channels. However, as GO resource reflects our current biological knowledge many of these candidate genes remain to be functionally characterized. By refining from 13,000 transcripts expressed per DRG neuron to only hundreds of candidates, this dataset represents a valuable resource for the identification of genes involved in the generation of DRG MS currents.

DISCUSSION

Patch-clamping of cultured DRG neurons is compatible with scRNAseq classification

Direct relationship between mechanosensory properties and molecular phenotypes of DRG neurons determined by whole transcriptome analysis has never been explored. To identify molecular features associated with specific mechanosensory functions and phenotypes, we have combined scRNA-seq profiling with electrophysiological recordings of individual DRG neurons. We provided full transcriptome data from single neurons in culture after characterization of MS currents, which were classified based on their inactivation kinetics. Coupling patch-clamp recording with mechanical stimulation required the use of adhering cultured DRG neurons. This approach is labor-intensive and low throughput and differs somehow from previous scRNAseq classification studies for which DRG neurons were collected shortly after ganglionic dissociation [START_REF] Li | Somatosensory neuron types identified by high-coverage single-cell RNA-sequencing and functional heterogeneity[END_REF][START_REF] Sharma | The emergence of transcriptional identity in somatosensory neurons[END_REF][START_REF] Usoskin | Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencing[END_REF][START_REF] Zeisel | Molecular Architecture of the Mouse Nervous System[END_REF]). Yet, 85% of the processed mRNA samples met quality standards for RNAseq. The molecular classification of our small and heterogeneous group of DRG neurons however was challenging because of the resulting low statistical power. Advances in microfluidics for single cell partitioning have enabled scRNAseq studies performed on samples of over one thousand cells, subdividing DRG neurons into a total of up to 17 transcriptional clusters [START_REF] Sharma | The emergence of transcriptional identity in somatosensory neurons[END_REF][START_REF] Zeisel | Molecular Architecture of the Mouse Nervous System[END_REF]. To achieve high-quality classification of our neuron sample, we combined our single-cell transcriptomes with the larger single-DRG neuron dataset from Zeisel and collaborators [START_REF] Zeisel | Molecular Architecture of the Mouse Nervous System[END_REF]. Therefore, despite its modest size, we succeeded to map all the neurons from our sample into DRG neuron transcriptional clusters. They are distributed among 11 of the 17 transcriptional clusters corresponding to 7 of the 8 main classes of neurons. The clusters which are not represented are PSNF2 (from NF neurons), PSNP2 (from NP1 neurons), PSPEP4 (from PEP1 neurons) and the three clusters PSPEP6-8, which together constitute the TrpM8 population. The absence of TrpM8 neurons in our sample is consistent with their lack of role in mechanosensation [START_REF] Knowlton | A sensory-labeled line for cold: TRPM8-expressing sensory neurons define the cellular basis for cold, cold pain, and cooling-mediated analgesia[END_REF]. It is possible that PSNP2 and PSPEP4 constitute subgroups of NP1 and PEP1 neurons that are not MS, but our relatively small sample of neurons does not allow a definite conclusion. However, the absence of PSNF2 neurons, thought to represent a cluster of proprioceptors [START_REF] Usoskin | Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencing[END_REF] likely stems from an undersized sample. Yet, our results attest that the culture of neurons, coupled to patch-clamp recording, does not prevent their classification and offers an extensive view of the main populations of DRG MS neurons. Thus, patch-clamp recording of MS currents is compatible with the production of a quantitative dataset that resolves mRNAs for all known genes in single mechanosensory neurons.

Distribution of MS current types in transcriptionally defined neurons

The distribution of MS currents in DRG neurons has been addressed over the past by distinguishing neuronal subsets using morphological and electrophysiological parameters as well as IB4 staining [START_REF] Coste | Pharmacological Dissection and Distribution of NaN/Nav1.9, T-type Ca2+ Currents, and Mechanically Activated Cation Currents in Different Populations of DRG Neurons[END_REF][START_REF] Drew | Acid-sensing ion channels ASIC2 and ASIC3 do not contribute to mechanically activated currents in mammalian sensory neurones[END_REF][START_REF] Drew | Distinct mechanosensitive properties of capsaicinsensitive and -insensitive sensory neurons[END_REF][START_REF] Hu | Mechanosensitive currents in the neurites of cultured mouse sensory neurones[END_REF][START_REF] Lechner | Peripheral sensitisation of nociceptors via G-proteindependent potentiation of mechanotransduction currents[END_REF][START_REF] Patkunarajah | TMEM87a/Elkin1, a component of a novel mechanoelectrical transduction pathway, modulates melanoma adhesion and migration[END_REF]. These approaches led to the view that large diameter, non-nociceptive neurons express mainly RA MS currents whereas medium-and small-diameter nociceptors express MS currents with slower kinetics. Subsets of DRG neurons labeled in vivo using genetically modified mice confirmed the tendency of proprioceptors and non-nociceptive A-type LTMRs to express RA MS currents [START_REF] Woo | Piezo2 is the principal mechanotransduction channel for proprioception[END_REF][START_REF] Zheng | Deep Sequencing of Somatosensory Neurons Reveals Molecular Determinants of Intrinsic Physiological Properties[END_REF], and showed that C-LTMR express MS currents with rapid and/or slow kinetics [START_REF] Delfini | TAFA4, a chemokine-like protein, modulates injury-induced mechanical and chemical pain hypersensitivity in mice[END_REF][START_REF] Lou | Runx1 controls terminal morphology and mechanosensitivity of VGLUT3-expressing C-mechanoreceptors[END_REF]. The accurate characterization of four distinct types of DRG MS currents in classified neurons using scRNAseq enabled to study their distribution in DRG neuron subsets.

Our results showed that the RA MS current is the main subtype in NF neuronal population, in line with the remarkable phenotype of Piezo2 KO animals related to innocuous touch and proprioception [START_REF] Florez-Paz | A critical role for Piezo2 channels in the mechanotransduction of mouse proprioceptive neurons[END_REF][START_REF] Ranade | Piezo2 is the major transducer of mechanical forces for touch sensation in mice[END_REF][START_REF] Woo | Piezo2 is the principal mechanotransduction channel for proprioception[END_REF].

Consistent with Piezo2 modest contribution to mechanical nociception [START_REF] Murthy | The mechanosensitive ion channel Piezo2 mediates sensitivity to mechanical pain in mice[END_REF], RA current type contributes to a lesser extent to MS currents in NP1, PEP1 and PEP2 neuronal subpopulations. The PEP1 neurons express predominantly the IA MS current type. This suggests that IA MS current may be the main transducer of noxious mechanical stimuli in C-type peptidergic nociceptors. The contribution of SA and ultra-SA currents in these neurons seems minimal.

Interestingly, the opposite is true in PEP2 neurons, in which ultra-SA current type contributes the most, whereas IA current type is almost absent. This suggests that Aδ-nociceptors and C-type peptidergic nociceptors use distinct sets of MS ion channels for mechanical pain detection. These differences could contribute to the signaling of distinct mechanosensory modalities by PEP1 and PEP2 nociceptors.

Importantly, RA current type is the main contributor of MS currents in NP2 neurons. These neurons, which have been shown to respond to mechanical stimulus, heat, and a large array of pruritogens, mediate itch behavior but not pain [START_REF] Han | A subpopulation of nociceptors specifically linked to itch[END_REF]. Therefore, these results suggest that Piezo2 actively contributes to itch sensing, by opposition to its expression in Merkel cells in which Piezo2 inhibition produces alloknesis, a pathological itch sensation produced by innocuous touch [START_REF] Feng | Piezo2 channel-Merkel cell signaling modulates the conversion of touch to itch[END_REF]).

An elegant study by Von Buchholtz and colleagues mapped mechanosensory functions to transcriptomic classes by combining in vivo calcium imaging and multiplexed ISH [START_REF] Von Buchholtz | Decoding Cellular Mechanisms for Mechanosensory Discrimination[END_REF]. Notably, complete loss of Aβ LTMR responses and loss of C-LTMRs responses to gentle but not noxious stimulations are seen in Piezo2 deleted mice. Although this study was performed on trigeminal ganglion neurons, which differ notably from DRG neurons [START_REF] Lopes | The Molecular Fingerprint of Dorsal Root and Trigeminal Ganglion Neurons[END_REF][START_REF] Megat | Differences between Dorsal Root and Trigeminal Ganglion Nociceptors in Mice Revealed by Translational Profiling[END_REF][START_REF] Price | Critical evaluation of the colocalization between calcitonin generelated peptide, substance P, transient receptor potential vanilloid subfamily type 1 immunoreactivities, and isolectin B4 binding in primary afferent neurons of the rat and mouse[END_REF], our results are consistent with these observations. Indeed, we detected almost only RA currents in NF neurons related to Aβ LTMR, whereas both RA and ultra-SA currents are expressed in C-LTMR related TH neurons. This suggests that the detection of high-threshold mechanical stimulus in C-LTMRs is mediated by ultra-SA currents.

No evidence for Tmem120a and Tmem150c in sustaining slow MS currents

While a compelling body of evidence indicates that Piezo2 generates RA-type MS current in DRG neurons, the molecular identity of genes coding for slow MS currents has remained elusive [START_REF] Delmas | Mechano-gated ion channels in sensory systems[END_REF][START_REF] Ranade | Piezo1, a mechanically activated ion channel, is required for vascular development in mice[END_REF]. Recently, new gene candidates have been proposed to encode slow-type MS currents in DRG neurons. Tmem120a has been suggested to be expressed in small-to medium-diameter neurons, and to mediate the ultra-slowly adapting MS current (Beaulieu-Laroche et al., 2020). Our results do not support this conclusion. First, our single-cell transcriptome analysis indicates that Tmem120a expression is not restricted to nociceptors and is not enriched in slowly or ultra-slowly adapting MS current-expressing neurons. This is consistent with our analysis of the crosssectional area distribution of Tmem120a-positive DRG neurons in ISH experiments. When normalized to the overall population, Tmem120a was detected in all subtypes of DRG neuron cell body diameters, including the largest ones expressing NF200. Nevertheless, the possibility that Tmem120 contributes to MS currents without being selectively expressed/enriched in a particular population was still open. However, our Tmem120a siRNA experiments performed either on the whole DRG neuron population or restricted to IB4 + neurons showed no changes in the proportion of slowly or ultra-slowly MS currentexpressing neurons, nor in the mean amplitude of these current components. Our siRNA assays beneficiated from our extensive kinetics classification of MS currents and from the inclusion of a piezo2 siRNA effective target knockdown. Thus, the broad expression of Tmem120a in DRG neurons, along with the lack of appreciable effects of Tmem120a siRNA treatment on endogenous MS currents, cast doubt about its role in mediating SA/ultra-SA MS currents. Moreover, overexpression of mouse Tmem120a in heterologous system did not induce MS channel activity. However, our data do not question Tmem120a involvement in pain sensing, which has been shown recently in inflammatory, but not neuropathic, mechanical hyperalgesia in rats [START_REF] Bonet | Involvement of TACAN, a Mechanotransducing Ion Channel, in Inflammatory But Not Neuropathic Hyperalgesia in the Rat[END_REF]. The molecular function of Tmem120a needs to be clarified in future work.

Tmem150c has been proposed to encode slowly adapting MS channels in DRG proprioceptive neurons [START_REF] Hong | Tentonin 3/TMEM150c Confers Distinct Mechanosensitive Currents in Dorsal-Root Ganglion Neurons with Proprioceptive Function[END_REF] but this inference has been challenged. Indeed, we and others have shown that, contrary to what was initially reported, Tmem150c expression in heterologous system does not induce MS currents by itself [START_REF] Anderson | TMEM150C/Tentonin3 Is a Regulator of Mechano-gated Ion Channels[END_REF][START_REF] Dubin | Endogenous Piezo1 Can Confound Mechanically Activated Channel Identification and Characterization[END_REF], but instead modulates the properties of MS ion channels when co-expressed [START_REF] Anderson | TMEM150C/Tentonin3 Is a Regulator of Mechano-gated Ion Channels[END_REF]. In addition, no slowly adapting current type was detected in proprioceptive neurons labeled in vivo using genetically modified mice [START_REF] Woo | Piezo2 is the principal mechanotransduction channel for proprioception[END_REF]. Tmem150c was found to be broadly expressed in our single-cell transcriptome analysis, and our siRNA experiments showed that Tmem150c knock-down caused no changes in the proportion or amplitude of any MS currents, suggesting that Tmem150c is not a prime contributor to MS currents in DRG neurons. To note, Tmem150c has been proposed recently to contribute to slowly adapting MS currents in nodose ganglion baroreceptors [START_REF] Lu | Tentonin 3/TMEM150C senses blood pressure changes in the aortic arch[END_REF]. Future studies will help resolve these controversies.

Molecular bases of the kinetically distinct MS current subtypes in DRG

Intensive efforts have been deployed to identify mammalian MS channels for decades, since the earliest indications of their existence in auditory hair cells [START_REF] Corey | Response latency of vertebrate hair cells[END_REF]. Identifying MS channel genes represents a crucial challenge for the field of nociception and may have many clinical implications in acute as well as inflammatory and chronic pain. In addition, it would be potentially relevant beyond somatosensation as exemplified by PIEZO channel implication in many physiological functions and diseases in humans [START_REF] Ma | Common PIEZO1 Allele in African Populations Causes RBC Dehydration and Attenuates Plasmodium Infection[END_REF][START_REF] Ma | A role of PIEZO1 in iron metabolism in mice and humans[END_REF][START_REF] Marshall | PIEZO2 in sensory neurons and urothelial cells coordinates urination[END_REF][START_REF] Mcmillin | Mutations in PIEZO2 cause Gordon syndrome, Marden-Walker syndrome, and distal arthrogryposis type 5[END_REF]. Therefore, by correlating DRG MS current properties with co-expressed genes at the single cell level, our study provides a valuable resource for identifying new molecular entities involved in mechanotransduction.

Our extensive characterization of MS current properties from large neuronal samples showed that 4 distinct types of MS currents are present in mouse DRG neurons. It is now well established that Piezo2 encodes MS channels responsible for RA currents. Does Piezo2 contribute to other types of MS currents? It has been shown that overexpression of Tmem150c with Piezo2 in heterologous systems induces slowing of PIEZO2 current inactivation to values close to IA current type [START_REF] Anderson | TMEM150C/Tentonin3 Is a Regulator of Mechano-gated Ion Channels[END_REF]. However, Piezo2 deletion or knock-down studies in DRG neurons, including this work, have shown no alteration of IA currents [START_REF] Coste | Piezo1 and Piezo2 are essential components of distinct mechanically activated cation channels[END_REF][START_REF] Ranade | Piezo2 is the major transducer of mechanical forces for touch sensation in mice[END_REF], and a possible compensatory increase in IA current responses has been reported in proprioceptors [START_REF] Woo | Piezo2 is the principal mechanotransduction channel for proprioception[END_REF]. This suggests that Piezo2 does not encode MS channels responsible for IA currents in DRG neurons.

The situation differs in trigeminal corneal neurons where a non-specific decrease of currents with intermediate and slow kinetics has been reported in Piezo2 KO mice [START_REF] Fernandez-Trillo | Piezo2 Mediates Low-Threshold Mechanically Evoked Pain in the Cornea[END_REF]. We also observed a modest effect of Piezo2 inhibition on the amplitude of ultra-SA currents. It cannot be excluded that this decrease is due to a RA current component masked in neurons with large ultra-SA currents. However, this modest impairment is still observable when current amplitude was measured at the end of the mechanical stimulus, i.e. when RA currents are expected to be fully inactivated (data not shown). Inhibition of a Piezo2 splice variant could be responsible for this effect, as Piezo2 extensive alternative splicing affecting inactivation kinetics has been reported in DRG neurons, although Piezo2 splicing variants displayed inactivation kinetics still in the range of RA current subtype [START_REF] Szczot | Cell-Type-Specific Splicing of Piezo2 Regulates Mechanotransduction[END_REF]. Whatsoever, this effect is slight and does not alter the proportion of neurons expressing the ultra-SA current, arguing that the primary contributor of ultra-SA currents is distinct from PIEZO2.

Notwithstanding, the nature of the molecular entities involved in IA-, SA-and ultra SA-currents remains unknown. In addition, the question of whether these currents are generated by multiple or single molecular entities is not yet resolved. Indeed, proximate environment such as membrane composition, attachment to extracellular matrix or subcortical cytoskeleton organization can potentially modify membrane bilayer tension and inserted MS channels properties including kinetics [START_REF] Cox | Removal of the mechanoprotective influence of the cytoskeleton reveals PIEZO1 is gated by bilayer tension[END_REF][START_REF] Del Marmol | Piezo1 forms a slowlyinactivating mechanosensory channel in mouse embryonic stem cells[END_REF][START_REF] Syeda | Piezo1 Channels Are Inherently Mechanosensitive[END_REF][START_REF] Teng | The force-from-lipid (FFL) principle of mechanosensitivity, at large and in elements[END_REF]. In this scenario, accessory subunits conferring distinct inactivation properties to a unique MS channel could be encoded by genes found to be enriched in a subset of DRG neurons.

Evidence suggests that the different currents are sustained by molecularly distinct ion channels. Indeed, three types of cationic MS ion channels with sustained activity but distinguishable from their unitary conductance have been described in DRG neurons [START_REF] Cho | A novel mechanosensitive channel identified in sensory neurons[END_REF][START_REF] Cho | Mechanosensitive ion channels in cultured sensory neurons of neonatal rats[END_REF].

Moreover, ultra-SA currents are blocked selectively by the conopeptide NMB1 [START_REF] Drew | High-threshold mechanosensitive ion channels blocked by a novel conopeptide mediate pressure-evoked pain[END_REF], suggesting that this current is sustained by a specific pool of MS channels. Hence, several MS channel coding genes remain to be discovered in DRG neurons. By identifying the genes enriched in DRG neurons with distinct mechanical phenotypes, our study lays the foundations for the identification of their molecular entities. Testing these candidates using siRNA experiments coupled with patch-clamp recording in DRG neurons may identify genes involved in mechanosensation. Importantly, given the broad involvement of mechanotransduction in physiological functions and the dynamism of the mechanobiology field, illustrated by the development of novel mechanobiology assays [START_REF] Kurth | A new mechanobiological era: microfluidic pathways to apply and sense forces at the cellular level[END_REF][START_REF] Xu | GPR68 Senses Flow and Is Essential for Vascular Physiology[END_REF], our data could also be used to identify mechanosensors or associated proteins in other biological tissues.

Limitation of the study

Transcriptomics data were generated from DRG neurons cultured in vitro. This should be considered when comparing these data to those obtained in native state, as these conditions can potentially alter the gene expression profile of cells and reduce cellular heterogeneity of the donor tissue. Performing patch-seq experiments in DRG tissue slices could circumvent this limitation.

STAR METHODS

Detailed methods are provided in the online version of this paper.

Supplemental Excel table

Table S2. Set of transcripts enriched in neurons grouped by MS current types, related to Figure 6.

STAR METHODS

RESOURCE AVAILABILITY

Lead Contact

Further information and requests for reagents and resources should be directed to and will be fulfilled by the Lead Contact, Bertrand Coste (bertrand.coste@univ-amu.fr).

Materials Availability

This study did not generate new unique reagents.

Data and Code Availability

• The single-cell RNA-seq datasets generated and analyzed have been deposited in NCBI's Gene Expression Omnibus (Edgar et al., 2002) and are publicly available through the following GEO Series accession numbers: GSE168032.

• This paper does not report original code.

• Any additional information required to reanalyze the data reported in this work paper is available from the Lead Contact upon request.

EXPERIMENTAL MODEL AND SUBJECT DETAILS

Animals

All animals were used in accordance with the European Community guidelines in the care and use of animals (2010/63/UE). Male mice (C57Bl/6J background) were 10-12 weeks old. Animals were housed in a 12h light /12h dark cycle and fed ad libitum. This project was approved by the Institutional Review Board of the regional ethic committee.

Dorsal root ganglion neuron culture

Mice were anesthetized with isoflurane and killed by decapitation. DRG were isolated and incubated in enzyme solution containing 2 mg.ml -1 of collagenase 1A (Sigma-Aldrich) and 5 mg.ml -1 of dispase II (ThermoFisher) for 45 min at 37 °C. The tissue was washed several times and triturated in Hanks' balanced salt solution (HBSS, Gibco). The resulting suspension was centrifuged (300×g for 5 min) in 1 ml HBSS and 1 ml Bovine Serum Albumin (BSA, 15%, Sigma-Aldrich). The pellet was rinsed with Dulbecco's Phosphate Buffered Saline (DPBS, Gibco) and centrifuged (300×g for 2 min). Cells were then plated onto poly-D-lysine-coated 12-mm round glass coverslips (Corning) in 24-well plates and were grown in Dulbecco's modified Eagle's medium (DMEM, High glucose, pyruvate, ThermoFisher) supplemented with 10% heat-inactivated Fetal Bovine Serum (FBS, Gibco), 100 U.ml -1 (1%) penicillin-streptomycin (Gibco), 100 ng.ml -1 nerve growth factor (NGF, Merk Millipore), and 50 ng.ml -1 glialderived neurotrophic factor (GDNF, ThermoFisher). Cells were recorded at 3-5 days in vitro (DIV).

Cell lines

NIH/3T3 cells and PIEZO1-deficient HEK293T cells [START_REF] Dubin | Endogenous Piezo1 Can Confound Mechanically Activated Channel Identification and Characterization[END_REF] were grown in DMEM supplemented with 10% heat-inactivated FBS (Gibco) and 100 U.ml -1 (1%) penicillin-streptomycin (Gibco) and cultured in 37°C with 5% CO2. Cell passage was done when cells reached 80% confluency.

For recordings, cells were plated onto poly-D-lysine-coated 12-mm round glass coverslips (Corning) in 24-well plates.

METHOD DETAILS Electrophysiology

Patch-clamp experiments were performed under whole-cell configuration using an Axopatch 200B amplifier (Axon Instruments). Patch pipettes had resistances of 2-3 MΩ when filled with an internal solution consisting of (in mM) 140 CsCl, 10 Hepes, 5 EGTA, 1 CaCl2, 1 MgCl2, 4 MgATP and 0.4 Na2GTP (pH adjusted to 7.3 with CsOH). Internal solution was supplemented with 0.4 U/µl of RNaseOUT (Invitrogen) for cell harvesting experiments. The extracellular solution consisted of (in mM) 133 NaCl, 3 KCl, 1 MgCl2, 10 Hepes, 2.5 CaCl2, 10 glucose (pH adjusted to 7.3 with NaOH). All experiments were done at room temperature. Currents were sampled at 20 kHz and filtered at 2 kHz.

Mechanical stimulation

Mechanical stimulation was achieved using a firepolished glass pipette (tip diameter 3-4 μm) positioned at an angle of 80° and in contact with the cell being recorded. Downward movement of the probe toward the cell was driven by a Clampex controlled piezo-electric crystal microstage (E625 LVPZT Controller/Amplifier, Physik Instrumente). The probe had a velocity of 0.7 μm.ms -1 during the ramp segment of the command for forward motion, and the stimulus was applied for 150 ms. Inward MS currents were recorded at a holding potential of -80 mV.

MS currents were characterized by applying series of 0.5 μm incremental steps every 10 s. This was done up to patch rupture, except for neurons dedicated to scRNAseq experiments. For this set of experiments, mechanical stimulation of neurons was stopped after MS currents with a peak amplitude of more than 100 pA were elicited in at least two consecutive incremental steps, preserving the integrity of the cellular content before the harvesting. The 3 non-responsive neurons included in our scRNAseq sample have been mechanically stimulated up to 10 µm probe displacement without disrupting the patch.

Pressure-clamp experiments

Pressure-clamp experiments were performed under cell-attached configuration. Patch pipettes had resistances of 2-3 MΩ when filled with a solution consisting of (in mM) 130 NaCl, 5 KCl, 10 HEPES, 1 CaCl2, 1 MgCl2, 10 TEA-Cl (pH 7.3 with NaOH). External solution used to zero the membrane potential consisted of (in mM) 145 KCl, 10 HEPES, 1 MgCl2, 10 glucose (pH 7.3 with KOH). Currents were sampled at 20 KHz and filtered at 2 KHz. Membrane patches were stimulated with negative pressure pulses through the recording electrode using a Clampex controlled pressure clamp HSPC-1 device (ALAscientific). Stretch-activated channels were recorded at a holding potential of -80 mV.

Single-cell library preparation

After electrophysiological recording, the cytoplasmic content of the neuron was harvested into the patch pipette, transferred into lysis solution (RNaseOUT 8 U.µl -1 , 2% Triton X-100) and flash frozen.

Reverse transcription was performed using a SMART-seq v4 Low Input Kit (Clontech) directly on the cell lysate according to the manufacturer's protocol. After cDNAds amplification, 50 μl of the sample were subjected to cDNA purification on AMPure XP beads. 0.5 ng of each purified cDNAds were used to construct the sequencing library using Nextera XT kit (Illumina) with fragments over 300-bp length from each neuron.

Processing, quality control and filtering of single-cell RNA-seq data cDNAds libraries were sequenced either in a NextSeq 500 (Illumina) using paired-end reads (76 pb) with an average depth of 29 million reads per cell. The raw reads (FastQ files) were cleaned by removing adaptor sequences, short sequences (length < 25 bp), low-quality bases (quality score < 30) and ambiguous sequences (i.e., reads with more than 30% unknown bases 'N') using CutAdapt 1.9.1 software [START_REF] Martin | Cutadapt removes adapter sequences from high-throughput sequencing reads[END_REF]. TopHat-2 v2.1.0 was used to map the cleaned RNA-seq reads to the mouse genome (mm10 / GRCm38.90) with two mismatches, two gaps and one multihit allowed. TopHat-2 splicing algorithm was used to map reads covering splice junctions, thereby improving the utilization of reads. After genome mapping, gene counts matrix was determined using GenomicAlignment (version 1.24.0) [START_REF] Lawrence | Software for computing and annotating genomic ranges[END_REF] using GTF annotation. A set of mitochondrial genes (mt-Atp6, mt-Atp8, mt-Co1, mt-Co2, mt-Co3, mt-Cytb, mt-Nd1, mt-Nd2, mt-Nd3, mt-Nd4, mt-Nd4l, mt-Nd5, mt-Nd6) were used to check viability of each sample: the sample is considered as non-apoptotic if the expression of these genes is less than 10% of total count in the same cell. Read counts were analyzed in the R/Bioconductor environment (version 3.12; www.bioconductor.org) with the R package edgeR.

The gene expression value was normalized between samples by TMM with a gene detection threshold of 20 counts in at least two cells. For the generation of heat maps shown in Figures 2B, z-transformed TMM count data was used. Relative expression levels were obtained by sample normalization using TPM (transcripts per kilobase million).

Dorsal root ganglion neuron clustering

The mapping of the scRNAseq neurons in DRG neuron transcriptional clusters was done by a crossdataset normalization approach. We used fastMNN, the new implementation of Mutual Nearest Neighbours (MNN) correction [START_REF] Haghverdi | Batch effects in single-cell RNAsequencing data are corrected by matching mutual nearest neighbors[END_REF] implemented in the R package batchelor to combine our dataset with Zeisel and collaborators scRNAseq sensory neuron dataset [START_REF] Zeisel | Molecular Architecture of the Mouse Nervous System[END_REF]. Briefly, a principal component analysis (PCA) was performed on the 211 most highly variable genes in DRG neurons (most highly enriched genes for each cluster) [START_REF] Zeisel | Molecular Architecture of the Mouse Nervous System[END_REF] to correct for batch effect. Identification of MNN was done in this reduced dimension space using the following parameters: dimension = 60 and k = 10. To assign each individual neuron to a specific cluster, we built a k nearest-neighbors (kNN) graph with k = 28 using Euclidean distance. Individual neurons were first assigned to one of the main cluster: NF, NP or PEP. Next, we repeated the procedure (PCA implemented in MNN correction, kNN clustering) for each main cluster. We were able to assign each cell to a specific sub-cluster: PSNF1, PSNF2, PSNF3, PSNP1, PSNP2, PSNP3, PSNP4, PSNP5, PSNP6, PSPEP1, PSPEP2, PSPEP3, PSPEP4, PSPEP5, PSPEP6, PSPEP7 or PSPEP8. Corrected expression values obtained during this procedure were only used for cluster assignment and visualization with tdistributed stochastic neighbor embedding (tSNE) projection (perplexity = 30, max iteration = 1000). Downstream analysis was performed using custom R scripts.

Electroporation

DRG cultures were electroporated using the Neon® Transfection System (Invitrogen) according to the manufacturer's protocol. After the centrifugation prior to cell plating, the cell pellet was resuspended in Neon® Transfection System resuspension Buffer R (Invitrogen) and gently mixed with 30 ng.µl -1 of pIRES2-AcGFP1 plasmid (Aequorea coerulescens GFP, Clontech; referred in results as GFP-plasmid) and 250 nM siRNA pool. Ten microliters of the cell suspension were then electroporated with the following program: 1200 V, 2 pulses, 20 ms. Electroporated cells were then plated onto poly-D-Lysine-coated 12mm round glass coverslips in 24-well plates in antibiotics-free media. Three hours later, the medium was replaced by DMEM supplemented with 10% heat-inactivated FBS, 1% penicillin-streptomycin, 100 ng.ml -1 NGF and 50 ng.ml -1 GDNF. Pools of 4 siRNA targeting Piezo2 or Tmem120a, or non-targeting siRNA for control experiments were purchased at Horizon Discovery (#L-163012-00-0005, #L-040281-01-0005 and #D-001810-10-05 for Piezo2, Tmem120a and control siRNA, respectively) and the pool of 4 siRNA targeting Tmem150c was purchased at Qiagen (#GS231503). Media were changed 48 hours later. Transfected DRG neurons were visualized by the expression of GFP. For experiments performed on IB4-positive neurons, neurons were incubated with 1 µg.ml -1 isolectin GS-IB4 from Griffonia simplicifolia, AlexaFluor 568 conjugate (Invitrogen) at 37°C for 10 minutes before recording. NIH/3T3 cells were electroporated using the Neon® Transfection System (Invitrogen) according to the manufacturer's protocol. Briefly, prior to cell plating, the cell pellet was resuspended in Neon® Transfection System Buffer R at a cell density of 1 ×10 7 cells.ml -1 . 250 nM siRNA was added to 10 µl of cell suspension and electroporation was done with the following program: 1400 V, 1 pulse, 30 ms. Electroporated cells were then transferred directly into 24-well tissue culture plates (1x10 5 cells per well) for RNA extraction, or onto poly-D-lysine-coated 12-mm round glass coverslips for patch-clamp recording experiments 72 hours later.

Molecular cloning of Tmem120a

Primers were designed from cDNA sequence of Tmem120a from the NCBI database (NM_172541). An 1.096 kb fragment was amplified from cDNA libraries generated from adult C57BL/6J DRG total RNA using primers fwd (5' GATTATGCATGCCGTGGACAAAGACATGCAGT 3') and rev (5' GATTATGGATCCTCAGTCCTTCTTGTTCCCGT 3') and cloned into pIRES2-AcGFP1 vector (Clontech, #632435) with NheI and BamHI restriction sites. The protein sequence of Tmem120a that we cloned is identical to NCBI Reference Sequence NP_766129.

Transfection of HEK-P1KO cells

PIEZO1-deficient HEK293T cells [START_REF] Dubin | Endogenous Piezo1 Can Confound Mechanically Activated Channel Identification and Characterization[END_REF] were transfected using lipofectamine 2000 (Invitrogen) according to the manufacturer's instruction. Each well was transfected with 600 ng.ml -1 of specified constructs. GFP positive cells were recorded 48 hours later.

qPCR experiments

Electroporated NIH/3T3 cells were harvested after 48 hours of incubation at 37°C in antibiotics-free growth medium. Total RNA was extracted using NucleoSpin RNA Mini kit (Macherey-Nagel) according to the manufacturer's instructions. DNA-free DNase Treatment (Ambion) was additionally used to remove genomic DNA contamination and iScript TM Reverse Transcription Supermix for RT-qPCR (Bio-Rad) was employed for cDNA synthesis. Samples generated with the iScript No-RT Control Supermix included in the kit were also synthesized for each sample as an additional verification step for the absence of genomic DNA contamination. DRG neurons were electroporated and cultured as described above. Forty-eight hours later, GFPexpressing DRG neurons were harvested individually and transferred into lysis solution (RNaseOUT 8 U.µl -1 , 2% Triton X-100). For each condition, 4 pools of 5 neurons were collected and flash frozen.

Reverse transcription was performed using a SMART-seq v4 Low Input Kit (Clontech) directly on the cell lysates according to the manufacturer's protocol.

Quantitative PCR (qPCR) was performed on an Applied Biosystems 7500 Fast Real-Time PCR System (Applied Biosystems) using KAPA SYBR® FAST qPCR Kit Master Mix (2X) Universal (Kapa Biosystems). A 12 µl reaction mix was made using 4 µl of diluted cDNA template (1 ng/µl for NIH/3T3 cells and 10 ng/µl for DRG neurons), specific primers (200 nM each), ROX Reference Dye Low (0.24 µl) and SYBR Green I Master Mix (6 µL) as recommended by the manufacturer. Cycling conditions were as follows:

95°C for 3 seconds, then 60°C for 30 seconds, 40 cycles. Reactions were performed in duplicate and melting-curve analysis was performed to assess the specificity of each amplification. Actb and Gapdh were selected as reference genes.

RT-qPCR amplicons produced in the presence or absence of RT enzyme were size-separated by electrophoresis on a 3.5% agarose gel stained with GelRed (Biotum), in comparison with ExactGene 100 bp ladder (Fisher BioReagents TM ) and photographed under UV light. The expected amplicon sizes are 144 bp for Tmem120a, 100 bp for Piezo1, 99 bp for Piezo2 and 85 bp for Tmem150c.

In Situ Hybridization experiments

Mice were anesthetized with isoflurane and killed by decapitation. Thoraco-lumbar vertebral spines were freshly collected, cut transversely, and transferred to a 4% sucrose solution (in 1X PBS) during 1h30 at 4°C, then to a 20% sucrose solution (in 1X PBS) overnight at 4°C. Then, DRG were isolated and embedded in Optimum Cutting Temperature medium (OCT; CellPath Ltd; #KMA-0100-00A) and cut into 14 µm sections with a CM350 S cryostat (Leica). In situ hybridization (ISH) was carried out using the RNAscope fluorescent multiplex reagent kit (Advanced Cell Diagnostics, a Bio-Techne brand, United States; #323133), according to the manufacturer's instructions. Briefly, tissue sections were post-fixed with 4% paraformaldehyde during 30 min at 4°C and rinsed with 1X PBS. After ethanol dehydration and a Protease IV pretreatment for 30 minutes at room temperature, sections were incubated with a mouse Tmem120a probe (ACD; #513211, accession number #NM_172541.2) for 2h at 40°C. The negative control, designed by ACD Bio, contains probes targeting the DapB gene from the Bacillus subtilis strain SMY (ACD; #320871, accession number #EF_191515). The signal was amplified with a succession of four pre-amplifiers supplied in the kit and detected with Alexa 488.

Immunofluorescence staining

NF-200 and peripherin immunostainings were done following ISH. After washing for 10 min at room temperature with PHEM buffer (consisting of in mM: 60 Pipes, 25 Hepes, 10 EGTA, 2 MgCl2) supplemented with 0.1% Triton X-100, non-specific sites were blocked with a reagent containing 3% Gelatin from cold water fish skin (Sigma-Aldrich; #G7765). Samples were incubated with primary antibodies against peripherin (1:200; rabbit polyclonal; Sigma-Aldrich; #AB1530) and NF-200 (1:1000; chicken polyclonal; AVES Labs; #NFH) overnight at 4°C. Primary antibodies were detected by using secondary antibodies conjugated to Alexa Fluor 647 (1 :1000 ; Donkey anti-rabbit polyclonal IgG ; Jackson ImmunoResearch, #711-605-152) and Alexa Fluor 568 (1:1000; Goat anti-chicken polyclonal IgY; Invitrogen, #A-11041), respectively, incubated for 1 h at room temperature. Sections were rinsed with PHEM buffer, dried and mounted in Prolong Gold Antifade Mountant (ThermoFisher; #P36930).

Gene Ontology (GO)

Gene Ontology (GO) annotations [START_REF] Ashburner | Gene ontology: tool for the unification of biology. The Gene Ontology Consortium[END_REF][START_REF] Ontology | The Gene Ontology resource: enriching a GOld mine[END_REF] were determined using the R package limma (v.3.44.3) and custom R scripts. GO terms were taken from the Bioconductor annotation package org. Mm.eg.db (v.3.11.4) released on 2020-05-02. The following GO terms were selected: detection of mechanical stimulus (GO:0050982), ion channel activity (GO:0005216) and integral component of membrane (GO:0016021). 𝐶. C was set to the baseline value of current once the mechanical stimulus is removed. Inactivation time constant τ was used to identify the type of MS current components. To note, τ measurements of ultra-SA currents are approximave due to the relative short lasting (150 ms) of the mechanical stimulus. When fitted with bi-exponential equation, the contribution (in percentage) of each component to the whole current is determined from the ratio of Ai to the amplitude of the peak current (see Fig. S2B). Recordings of MS currents displaying no clear deactivation once the mechanical stimulus is removed were discarded. We arbitrarily set a component detection threshold at 10% of the peak current, i.e. so that only detected components contributing to at least 10% of MS currents were included into analyzes, to avoid mis-identification of current components.

QUANTIFICATION AND STATISTICAL ANALYSIS

Mechanosensitive current analysis

Differential expression analysis

Analysis of differential gene expression was made between the cells expressing a given type of MS current (RA, IA, SA or ultra-SA current) versus cells that do not express this specific current (control group). LogP value is the Pvalue of Fold Change determined with exact test for negative binomial distribution implanted in the R package edgeR between two groups of count libraries. We set up a 2fold enrichment threshold (with Pvalue < 0.05) to consider genes displaying expression enrichment in the group of interest compared to the control group. For graphical representation, fold change is displayed as Log2 values (corresponding to Log2FC) and Pvalue as -log10 values (corresponding to -Log10 Pvalue).

Abundance of transcripts is expressed by the Transcript Per kilobase Million value (TPM). Mean TPM values in a group of neurons expressing a specific current type can be used to rank genes by transcript abundance in this group, which differs from fold change values that do not consider absolute level of expression but reflect expression difference between two groups. In figure 6, we used both parameters to present gene candidates, i.e. we listed among genes that are significantly enriched (FC > 2, Pvalue < 0.05) in the group of interest those that are the most abundantly expressed in this group (ranked by mean TPM).

Image Acquisition and Analysis

Images were acquired using a LSM 780 laser-scanning confocal microscope (Zeiss) with Ar and He-Ne lasers using a 63X oil-immersion objective. Images were treated and analysed with ImageJ software (v 1.50b, NIH, United States). The "analyze particle" function was used to quantify in situ hybridization signal. Neurons containing at least two dots were considered positive.

Statistical analysis

Unless mentioned otherwise in Method Details above, data were analyzed using GraphPad Prism 8. 1 0 0 -2 0 0 2 0 0 -3 0 0 3 0 0 -4 0 0 4 0 0 -5 0 0 5 0 0 -6 0 0 6 0 0 -7 0 0 7 0 0 -8 0 0 8 0 0 -9 0 0 9 0 0 -1 0 0 0 > 1 0 0 0 Cross-sectional area (µm²)

FIGURE LEGENDS Figure 1 .

 LEGENDS1 FIGURE LEGENDS Figure 1. Single cell RNAseq coupled to MS current characterization

Figure 2 .

 2 Figure 2. Cluster analysis of scRNAseq neurons

Figure 3 .

 3 Figure 3. Distribution of MS currents among transcriptional clusters of DRG neurons

Figure 4 .

 4 Figure 4. Expression analysis of genes proposed to encode MS channels

Figure 5 .

 5 Figure 5. Tmem120a and Tmem150c do not contribute to DRG MS currents

Figure 6 .

 6 Figure 6. Potential molecular determinants of DRG MS channels

  Clampfit 10.7 (Molecular Devices) software was used to analyze recordings and biophysical parameters. The decay of inactivation from the peak of current was fitted with mono or bi-exponential function of the form 𝑓(𝑡) = ∑ 𝐴 𝑖 𝑒 -

  , Ct values obtained for gene of interest (GOI) were normalized in relation to the average Ct values of Gapdh and Actb genes (RG) (ΔCt = CtGOI -Mean CtRG). The difference between the ΔCt of the sample of interest and the control sample (ΔΔCt = ΔCtsample of interest -ΔCtcontrol sample) was then calculated, and finally the 2-ΔΔCt value to express the results as relative expression levels of gene of interest in samples of interest versus the control sample (non-targeting siRNA).

2. 0 (

 0 figure legend and are two-sided. In all figure legends the exact value and definition of n is indicated. In all panels: not significant, ns p > 0.05; * p < 0.05; ** p < 0.01; *** p < 0.001.

  Figure 1

  Figure 5

Acknowledgements

This project has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (Grant agreement No. 678610). RNA sequencing was done by Profilexpert, "Génomique & Microgénomique" service, Lyon, FRANCE. We thank Virginie Penalba and Axel Fernandez for expert technical assistance, and Dr. Nancy Osorio, Dr.

Françoise Padilla and Dr. Nicolas Grillet for critical reading of the manuscript.

Declaration of interests

The authors declare no competing interests. 3 RNA sequenced DRG neurons are labeled according to their genetically identified neuronal type and their MS currents characterized by patch clamp experiments. Mixed currents are MS currents for which inactivation kinetics could be fitted with bi-exponential function, revealing the presence of two distinct MS components, with the smaller (current type 2) contributing to at least 10% of the peak current amplitude.