S Desiderio

F Schwaller

G R Lewin

P Carroll

F Marmigère

Touch receptor end-organ innervation and function requires sensory expression of the transcription factor Meis2

whether they are published or not. The documents may come

Touch receptor end-organ innervation and function requires sensory expression of the transcription factor Meis2

Abstract:

Touch sensation is primarily encoded by mechanoreceptors, sometimes called Low-Threshold Mechanoreceptors (LTMRs), with their cell bodies in the Dorsal Root Ganglia (DRG). LTMRs make up no more that 20% of all sensory neurons and exhibit great diversity in terms of molecular signature, terminal ending morphology and electrophysiological properties, mirroring the complexity of tactile experience. LTMRs are an interesting model to study the molecular cues controlling neuronal diversification in terms of both molecular specification and target-field innervation. The morphological specialization of the sensory end-organ of LTMRs exhibits striking diversity between different mechanoreceptor types and whether it occurs in the glabrous or hairy skin. Much has been learnt about transcriptional codes that define different LTMR subtypes, but the identification of molecular players that participate in their late maturation has not been extensively addressed. Here we identified the TALE homeodomain transcription factor Meis2 as a key regulator of LTMR targetfield innervation. Meis2 is specifically expressed in cutaneous LTMRs and its expression depends on target-derived signals. Meis2 gene inactivation in mouse sensory neurons precursors or early postmitotic neurons allows normal survival and specification of LTMRs. However, LTMRs lacking Meis2

show pronounced defects in end-organ innervation which was accompanied by severely impaired receptor properties and behavioral responses. These results establish Meis2 as a major transcriptional regulator controlling the orderly formation of peripheral end-organs required for touch.

Key words:

Dorsal root ganglion, Meis2, sensory neurons, ASD, transcription factor, TALE homeodomain.

Introduction

Tactile stimuli like brush, light pressure, or roughness engage highly specialized, but also diverse arrays of mechanoreceptors in the hairy and glabrous skin [START_REF] Delmas | Molecular mechanisms of mechanotransduction in mammalian sensory neurons[END_REF][START_REF] Li | The functional organization of cutaneous low-threshold mechanosensory neurons[END_REF][START_REF] Abraira | The sensory neurons of touch[END_REF][START_REF] Lechner | Hairy sensation[END_REF][START_REF] Zimmerman | The gentle touch receptors of mammalian skin[END_REF][START_REF] Handler | The mechanosensory neurons of touch and their mechanisms of activation[END_REF][START_REF] Wu | Distinct subtypes of proprioceptive dorsal root ganglion neurons regulate adaptive proprioception in mice[END_REF]. Sensory neurons within the DRG can be broadly classified as nociceptors, mechanoreceptors or proprioceptors and each group is characterized by the expression of specific combination of genes, have distinctive physiological properties and projections within the spinal cord and periphery [START_REF] Marmigère | Specification and connectivity of neuronal subtypes in the sensory lineage[END_REF][START_REF] Lallemend | Molecular interactions underlying the specification of sensory neurons[END_REF][START_REF] Vermeiren | Vertebrate Sensory Ganglia: Common and Divergent Features of the Transcriptional Programs Generating Their Functional Specialization[END_REF]. Recent advances in combining single cell transcriptomic and genetic tracing have tremendously extended the classical subtypes repertoire [START_REF] Bai | Genetic Identification of an Expansive Mechanoreceptor Sensitive to Skin Stroking[END_REF][START_REF] Usoskin | Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencing[END_REF][START_REF] Zheng | Deep Sequencing of Somatosensory Neurons Reveals Molecular Determinants of Intrinsic Physiological Properties[END_REF][START_REF] Sharma | The emergence of transcriptional identity in somatosensory neurons[END_REF][START_REF] Wu | Distinct subtypes of proprioceptive dorsal root ganglion neurons regulate adaptive proprioception in mice[END_REF].

Cutaneous mechanoreceptors or Low Threshold Mechanoreceptors (LMTRs) exhibit a variety of specialized terminal endings in the hairy and glabrous skin with strikingly unique morphologies [START_REF] Delmas | Molecular mechanisms of mechanotransduction in mammalian sensory neurons[END_REF][START_REF] Li | The functional organization of cutaneous low-threshold mechanosensory neurons[END_REF][START_REF] Abraira | The sensory neurons of touch[END_REF][START_REF] Lechner | Hairy sensation[END_REF][START_REF] Zimmerman | The gentle touch receptors of mammalian skin[END_REF][START_REF] Sharma | The emergence of transcriptional identity in somatosensory neurons[END_REF][START_REF] Handler | The mechanosensory neurons of touch and their mechanisms of activation[END_REF][START_REF] Schwaller | USH2A is a Meissner's corpuscle protein necessary for normal vibration sensing in mice and humans[END_REF]. LTMRs projecting to the glabrous skin innervate Merkel cell complexes or Meissner corpuscles at the dermal-epidermal border. LTMRs innervating Merkel cells in the glabrous or hairy skin have large thickly myelinated axons (Aβ-fibers) and are characterized as slowly-adapting mechanoreceptors responding to skin movement and static displacement (also referred to as Aβ-SAIs). Meissner corpuscles are mechanoreceptors which are only sensitive to skin movement or vibration (rapidly-adapting mechanoreceptors) and are referred to as Aβ-RAs. LTMRs innervating hair follicles in the hairy skin can form lanceolate endings or circumferential endings. Virtually all mechanoreceptors innervating hairs show rapidly-adapting properties and respond only to hair movement and not static displacement [START_REF] Lechner | Hairy sensation[END_REF]. LTMRs with large myelinated axons innervating hairy skin are characterized as Aβ-RAs, and a specialized population of slowly conducting myelinated fibers called D-hair mechanoreceptors (or Aδ-RAs) also form lanceolate endings on small hairs. D-hair mechanoreceptors are most sensitive to low velocity stroking have large receptive fields and are directionally tuned [START_REF] Li | The functional organization of cutaneous low-threshold mechanosensory neurons[END_REF][START_REF] Walcher | Specialized mechanoreceptor systems in rodent glabrous skin[END_REF]. A small number of LTMRs in the hairy skin are not activated by hair movement but show properties of rapidly-adapting mechanoreceptors [START_REF] Lewin | Physiological properties of primary sensory neurons appropriately and inappropriately innervating skin in the adult rat[END_REF]. These were originally characterized as so-called field receptors [START_REF] Burgess | Specific regeneration of cutaneous fibers in the cat[END_REF][START_REF] Lewin | Physiological properties of primary sensory neurons appropriately and inappropriately innervating skin in the adult rat[END_REF] and were recently shown to form circumferential endings around hair follicles [START_REF] Bai | Genetic Identification of an Expansive Mechanoreceptor Sensitive to Skin Stroking[END_REF]. LTMRs tuned to high frequency vibration are called Aβ-RAII and innervate Pacinian corpuscles deep in the skin or on the bone [START_REF] Schwaller | USH2A is a Meissner's corpuscle protein necessary for normal vibration sensing in mice and humans[END_REF]. Ruffini endings that are thought to be innervated by stretch sensitive mechanoreceptors (Aβ-SAII) are poorly characterized in mice [START_REF] Handler | The mechanosensory neurons of touch and their mechanisms of activation[END_REF]. Whereas genetic tools allow to identify at least 5 LTM neurons populations [START_REF] Zimmerman | The gentle touch receptors of mammalian skin[END_REF][START_REF] Bai | Genetic Identification of an Expansive Mechanoreceptor Sensitive to Skin Stroking[END_REF], a recent update combining these tools with deep RNA sequencing clustered at least 20 different subtypes of LTM neurons [START_REF] Zheng | Deep Sequencing of Somatosensory Neurons Reveals Molecular Determinants of Intrinsic Physiological Properties[END_REF].

Understanding how sensory neurons diversity is generated from an apparently generic progenitor and cracking the transcriptional codes supporting this diversification has been the object of tremendous efforts in the last decades [START_REF] Marmigère | Specification and connectivity of neuronal subtypes in the sensory lineage[END_REF][START_REF] Lallemend | Molecular interactions underlying the specification of sensory neurons[END_REF][START_REF] Usoskin | Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencing[END_REF][START_REF] Zheng | Deep Sequencing of Somatosensory Neurons Reveals Molecular Determinants of Intrinsic Physiological Properties[END_REF][START_REF] Sharma | The emergence of transcriptional identity in somatosensory neurons[END_REF][START_REF] Wu | Distinct subtypes of proprioceptive dorsal root ganglion neurons regulate adaptive proprioception in mice[END_REF]. Although the repertoire of transcription factor (TFs) of each subpopulation is largely established, their specific functions remain largely unexplored. The functions of specification factors or terminal selectors [START_REF] Hobert | Terminal Selectors of Neuronal Identity[END_REF][START_REF] Hobert | Neuronal identity control by terminal selectors in worms, flies, and chordates[END_REF], like Maf, Shox2, Runx3, Pea3 and ER81 have been functionally implicated in LTMR segregation [START_REF] Arber | ETS gene Er81 controls the formation of functional connections between group Ia sensory afferents and motor neurons[END_REF][START_REF] Marmigère | Specification and connectivity of neuronal subtypes in the sensory lineage[END_REF][START_REF] Abdo | Dependence on the transcription factor Shox2 for specification of sensory neurons conveying discriminative touch[END_REF][START_REF] Scott | Transcription factor short stature homeobox 2 is required for proper development of tropomyosin-related kinase B-expressing mechanosensory neurons[END_REF][START_REF] Hu | C-Maf is required for the development of dorsal horn laminae III/IV neurons and mechanoreceptive DRG axon projections[END_REF][START_REF] Lallemend | Molecular interactions underlying the specification of sensory neurons[END_REF]. Less is understood about the transcriptional control of target cell innervation within the skin and the establishment of specialized peripheral endorgan complexes. Interestingly, recent advances in the understanding of Autism Spectrum Disorder (ASD) suggest that centrally affected neurons in ASD and LTMRs share specific transcriptional programs regulating late neuronal differentiation [START_REF] Orefice | Peripheral Mechanosensory Neuron Dysfunction Underlies Tactile and Behavioral Deficits in Mouse Models of ASDs[END_REF][START_REF] Orefice | Targeting Peripheral Somatosensory Neurons to Improve Tactile-Related Phenotypes in ASD Models[END_REF]. Genes known to be involved in (ASD) are expressed by LTMRs, and mouse models for ASD exhibit sensory deficits [START_REF] Orefice | Peripheral Mechanosensory Neuron Dysfunction Underlies Tactile and Behavioral Deficits in Mouse Models of ASDs[END_REF][START_REF] Orefice | Targeting Peripheral Somatosensory Neurons to Improve Tactile-Related Phenotypes in ASD Models[END_REF]. Furthermore, DRG specific ablation of some of these genes causes ASDlike behaviors. Meis2, is another TF in which mutations are associated with severe ASD [START_REF] Gangfuß | Intellectual disability associated with craniofacial dysmorphism, cleft palate, and congenital heart defect due to a de novo MEIS2 mutation: A clinical longitudinal study[END_REF] in humans, and is also specifically expressed in LTMRs [START_REF] Usoskin | Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencing[END_REF][START_REF] Zheng | Deep Sequencing of Somatosensory Neurons Reveals Molecular Determinants of Intrinsic Physiological Properties[END_REF][START_REF] Sharma | The emergence of transcriptional identity in somatosensory neurons[END_REF]. Thus, because ASD is a highly genetically and phenotypically heterogeneous neurodevelopmental disorder in turn involving neurogenesis, neurite outgrowth, synaptogenesis and synaptic plasticity [START_REF] Mcfadden | Evidence for dysregulation of axonal growth and guidance in the etiology of ASD[END_REF][START_REF] Gilbert | Fundamental elements in autism: From neurogenesis and neurite growth to synaptic plasticity[END_REF]), Meis2 appears as a pertinent candidate to regulate late sensory neurons differentiation.

Meis2 belongs to a highly conserved homeodomain family of TFs containing three members in mammals, Meis1, Meis2 and Meis3 [START_REF] Geerts | The role of the MEIS homeobox genes in neuroblastoma[END_REF][START_REF] Longobardi | Biochemistry of the tale transcription factors PREP, MEIS, and PBX in vertebrates[END_REF]. We previously showed that Meis1 plays a role in target-field innervation of sympathetic peripheral neurons [START_REF] Bouilloux | Loss of the transcription factor Meis1 prevents sympathetic neurons target-field innervation and increases susceptibility to sudden cardiac death[END_REF]. Here, we have asked whether Meis2 could have a role in regulating the innervation of specialized cutaneous end-organs important for LTMR function. We confirmed that Meis2 expression is restricted to LTMR subclasses at late developmental stages compatible with a function in specification and/or target field innervation. Mice carrying Meis2 gene inactivation in post-mitotic sensory neurons are healthy and viable and do not exhibit any neuronal loss, but display severe tactile sensory defects as demonstrated physiologically and in behavioral assays. Consistent with these findings we found major morphological alterations in LTMR end-organ structures.

Results

:

Meis2 is expressed by cutaneous LTMRs.

We analyzed Meis2 expression using in situ hybridization (ISH) at various developmental stages in both mouse and chick lumbar DRG, combined with well-established molecular markers of sensory neuron subclasses (Figures 1A; Figures 1 Supplementary 1 and2). In mouse, Meis2 mRNA was first detected at embryonic day (E) 11.5 in a restricted group of large DRG neurons. This restricted expression pattern was maintained at E14.5, E18.5 and adult stages (Figure 1A). In chick, Meis2 was expressed in most DRG neurons at Hamburger-Hamilton stage (HH) 24, but later become restricted to a well-defined subpopulation in the ventro-lateral part of the DRG where LTMRs and proprioceptors are located [START_REF] Rifkin | Dynamic expression of neurotrophin receptors during sensory neuron genesis and differentiation[END_REF] (Figure 1 Supplementary 2A). In both species, Meis2positive cells also expressed the pan-neuronal marker Islet1, indicating that they are post-mitotic neurons. In chick, we estimated that Meis2-positive cells represented 16.4 ± 1.1 and 14.4 ± 1.2% (mean ± SEM; n = 3) of Islet1-positive DRG neurons at HH29 and HH36 respectively, suggesting a stable expression in given neuronal populations during embryonic development. Double ISH for Meis2 and Ntrk2, Ntrk3 or c-Ret mRNAs in E14.5 and E18.5 mouse embryonic DRG (Figure 1 Supplementary 1A and B) showed a large co-expression in Ntrk2and Ntrk3-positive neurons confirming that Meis2-positive neurons belong to the LTMR and proprioceptive subpopulations. We estimated that 16.0 ± 1.2, 28.7 ± 2.5 and 30.5 ± 3.5% (mean ± SEM; n = 3) of Meis2-positive neurons co-expressed Ntrk2 at E14.5, E18.5 and P8 respectively, and that 43.7 ± 2.1, 57.8 ± 5.1 and 39.5 ± 5.4% (mean ± SEM; n = 3) co-expressed Ntrk3 at these stages. Conversely, Meis2 was co-expressed in 55.9 ± 3.1, 52.1 ± 4.5 and 53.6 ± 9.4% (mean ± SEM; n = 3) of Ntrk2-positive neurons at E14.5, E18.5 and P8 respectively, and in 79.1 ± 3.4, 63.8 ± 5.1 and 78.5 ± 5.0% (mean ± SEM; n = 3) of Ntrk3positive neurons. Finally, double ISH for Meis2 and c-Ret in E14.5 mouse DRG showed that virtually all large c-Ret-positive neurons representing part of the LTMR pool co-expressed Meis2 at this stage before the emergence of the small nociceptive Ret-positive population. Similar results were found in chick at HH29 (Figure 1 Supplementary 2B) where we estimated that 38.4 ± 7.7 and 41.7 ± 3.0% (mean ± SEM; n = 3) of Meis2-positive neurons co-expressed Ntrk2 or Ntrk3 respectively, and inversely, that 62.1 ± 3.8 and 45.1 ± 4.7% (mean ± SEM; n = 3) of Ntrk2and Ntrk3-positive neurons co-expressed Meis2 mRNA respectively. In mouse, comparison of Meis2 mRNA expression to Ntrk1, a well-established marker for early nociceptive and thermo-sensitive neurons, showed that only few Meis2-positive neurons coexpressed Ntrk1 at E11.5 and E18.5 (Figure 1 Supplementary 1C). In chick HH29 embryos, Meis2 expression was fully excluded from the Ntrk1 [START_REF] Usoskin | Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencing[END_REF][START_REF] Zheng | Deep Sequencing of Somatosensory Neurons Reveals Molecular Determinants of Intrinsic Physiological Properties[END_REF][START_REF] Sharma | The emergence of transcriptional identity in somatosensory neurons[END_REF][START_REF] Shin | Intrinsic control of neuronal diversity and synaptic specificity in a proprioceptive circuit[END_REF].

The relatively lower co-incidence of Meis2 and Ntrk2 expressions compared to Ntrk3 is consistent with Meis2 being excluded from the Aδ-LTMRs (D-hair mechanoreceptors). The lack of co-expression with Ntrk1 and TrpV1 also confirmed Meis2 exclusion from peptidergic and non-peptidergic subpopulations.

The maintenance of Meis2 expression depends on target-derived signals.

The requirement for extrinsic signals provided by limb mesenchyme and muscles for proprioceptor and LTMR development has been documented [START_REF] Arber | ETS gene Er81 controls the formation of functional connections between group Ia sensory afferents and motor neurons[END_REF][START_REF] Patel | Peripheral NT3 signaling is required for ETS protein expression and central patterning of proprioceptive sensory afferents[END_REF][START_REF] De Nooij | Etv1 Inactivation Reveals Proprioceptor Subclasses that Reflect the Level of NT3 Expression in Muscle Targets[END_REF][START_REF] Poliak | Muscle-type Identity of Proprioceptors Specified by Spatially Restricted Signals from Limb Mesenchyme[END_REF][START_REF] Wang | Muscle-selective RUNX3 dependence of sensorimotor circuit development[END_REF][START_REF] Sharma | The emergence of transcriptional identity in somatosensory neurons[END_REF][START_REF] Shin | Intrinsic control of neuronal diversity and synaptic specificity in a proprioceptive circuit[END_REF]. Whereas Hox gene expression in sensory neurons does not depend on limb derived signals [START_REF] Shin | Intrinsic control of neuronal diversity and synaptic specificity in a proprioceptive circuit[END_REF], Pea3, ER81 and Runx3 expressions are lost when limb-derived signals are removed [START_REF] Lin | Functionally related motor neuron pool and muscle sensory afferent subtypes defined by coordinate ETS gene expression[END_REF][START_REF] Lallemend | Molecular interactions underlying the specification of sensory neurons[END_REF]. ER81 and Runx3 expressions are regulated by Ntf3/Ntrk3 signaling and in this model, different levels of Ntf3 expression provided by each muscle seem to contribute to proprioceptive neurons diversity [START_REF] Patel | Peripheral NT3 signaling is required for ETS protein expression and central patterning of proprioceptive sensory afferents[END_REF][START_REF] De Nooij | Etv1 Inactivation Reveals Proprioceptor Subclasses that Reflect the Level of NT3 Expression in Muscle Targets[END_REF][START_REF] Wang | Muscle-selective RUNX3 dependence of sensorimotor circuit development[END_REF]. Early after limb ablation, sensory neuron differentiation and Ntrk receptors expression was unaffected, but later survival of ventro-lateral LTMRs and proprioceptors neurons is compromised and can be rescued by addition of BDNF and NT3 [START_REF] Oakley | Muscle sensory neurons require neurotrophin-3 from peripheral tissues during the period of normal cell death[END_REF][START_REF] Calderó | Peripheral target regulation of the development and survival of spinal sensory and motor neurons in the chick embryo[END_REF]. To test the influence of target-derived signals on Meis2 expression in sensory neurons, limb buds were unilaterally ablated in HH18 chick embryos. Embryos were harvested at HH27 and HH36, before and after ventro-lateral neurons are lost respectively [START_REF] Oakley | Muscle sensory neurons require neurotrophin-3 from peripheral tissues during the period of normal cell death[END_REF][START_REF] Oakley | Neurotrophin-3 promotes the differentiation of muscle spindle afferents in the absence of peripheral targets[END_REF][START_REF] Calderó | Peripheral target regulation of the development and survival of spinal sensory and motor neurons in the chick embryo[END_REF] (Figure 1C andD; Figure 1 Supplementary 2D). As previously reported, about 30% of all sensory DRG neurons represented by the pan-neuronal marker Islet1 were lost in HH36 embryos on the ablated side compared to the contralateral side (Figure 1C). About 50% of Ntrk2 and 65% of Ntrk3 positive VL-neurons were lost, but Ntrk2-positive DL neurons were not significantly affected (Figure 1 Supplementary 2D). Consistently, at this stage, about 65% of Meis2positive neurons were lost (Figure 1C). In HH27 embryos, whereas no significant loss of Islet1-positive neurons was detected following limb ablation, about 40% of Meis2-positive neurons were lost, and remaining Meis2-positive neurons expressed very low levels of Meis2 mRNAs (Figure 1D). These results indicate that target-derived signals are necessary for the maintenance, but not the induction, of Meis2 expression in sensory neurons.

Meis2 gene inactivation in post-mitotic sensory neurons induces severe behavioral defects.

We next asked whether inactivation of the ASD-associated Meis2 gene would induce changes in LTMR structure and function.

We generated a conditional mouse mutant strain for Meis2 (Meis2

LoxP/LoxP

) in which the first coding exon for the homeodomain was flanked by LoxP sites (Figure 1 Supplementary 3A). This strategy was previously successfully used to specifically disrupt Meis1 function in the developing sympathetic nervous system [START_REF] Bouilloux | Loss of the transcription factor Meis1 prevents sympathetic neurons target-field innervation and increases susceptibility to sudden cardiac death[END_REF]. This strain was independently crossed with the Wnt1 [START_REF] Ericson | Early stages of motor neuron differentiation revealed by expression of homeobox gene Islet-1[END_REF][START_REF] Dasen | A Hox regulatory network establishes motor neuron pool identity and target-muscle connectivity[END_REF][START_REF] Catela | Hox Proteins Coordinate Motor Neuron Differentiation and Connectivity Programs through Ret/Gfrα Genes[END_REF], we did not observe obvious motor deficits in Isl1 Cre/+ ::Meis2 LoxP/LoxP mice. Thus in a catwalk analysis we found no differences in any of the gait parameters measured between WT and mutant mice (Figure 1 Supplementary Table 1).

We next monitored tactile evoked behaviors in adult WT, Isl1 Cre/+ and Isl1 Cre/+ ::Meis2 LoxP/LoxP mice, using stimuli applied to both glabrous and hairy skin. We used von Frey filaments to apply a series of low forces ranging from 0.16 to 0.6g to the hind paw and found the frequency of withdrawal responses to be significantly decreased in Isl1 +/Cre littermates (Figure 1F). In the von Frey test, the thresholds for paw withdrawal were similar between all genotypes when using filaments exerting forces ranging from 1 to 1.4g, which likely reflects the activation of mechanical nociception suggesting that Meis2 gene inactivation did not affect nociceptor function. We also used hot plate assay to assess noxious heat evoked behaviors and found no difference in jump latencies between WT, Isl1 Cre/+ and Isl1 Cre/+ ::Meis2 LoxP/LoxP mice (Figure 1G). Finally, we compared the sensitivity of mice to stimuli applied to the hairy skin using the sticky tape test. Placing sticky tape on the back skin evoked attempts to remove the stimulus in a defined time window and we found that such bouts of behavior were significantly reduced in Isl1 Cre/+ ::Meis2 Meis2 inactivation does not compromise LTMR and proprioceptor survival or specification during embryogenesis. Similar to the Wnt1 Cre ::Meis2 LoxP/LoxP strain, there was no difference in the size of the DRGs between E16.5 WT and Isl1 Cre/+ ::Meis2 LoxP/LoxP embryos as well as in the number of Ntrk2 and Ntrk3-positive neurons (Figure 2B) suggesting no cell loss. In E18.5 embryonic DRGs, the number of LTMR and proprioceptors identified as positive for Ntrk2, Ntrk3, Ret and Maf were unchanged following Meis2 inactivation (Figure 2C). Consistent with the lack of Meis2 expression in nociceptors, the number of Ntrk1-positive neurons was also unaffected (Figure 2C). At this stage, phosphoCreb (pCreb) expression in Ntrk2 and Ntrk3-positive neurons was similar in WT and mutants (Figure 2 Supplementary 2), suggesting that Ntrk signaling is not affected. Altogether, these results show that Meis2 is dispensable for LTMR and proprioceptor survival and specification.

Our assays in chick limb ablation indicated that like other TFs involved in target field innervation, the maintenance of Meis2 expression also depends on target-derived signals. In the developing sympathetic nervous system, we previously demonstrated that another member of the Meis family, Meis1, is involved in distal target-field innervation of sympathetic neurons but was dispensable for their specification [START_REF] Bouilloux | Loss of the transcription factor Meis1 prevents sympathetic neurons target-field innervation and increases susceptibility to sudden cardiac death[END_REF]. In P0 Wnt1 Cre ::Meis2 LoxP/LoxP , Nefh staining in the hind paws showed strong innervation deficits as reflected by a paucity of neurofilament-positive myelinated branches in both the glabrous and hairy skin (Figure 2D). LTMRs form specialized sensory endings in a variety of end-organs specialized to shape the mechanoreceptor properties. We next used the Isl1 +/Cre ::Meis2 LoxP/LoxP mice to assess the effects of late loss of Meis2 on LTMR structure and function. Meis2 gene inactivation compromises SA-LTMR morphology and function in the glabrous, but not the hairy skin.

We next investigated if post mitotic Meis2 inactivation impacts terminal morphologies and physiological properties of LTMRs. We made recordings from single mechanoreceptors and probed their responses to defined mechanical stimuli in adult WT, Isl1 +/Cre and Isl1 +/Cre ::Meis2 LoxP/LoxP mice using ex vivo skin nerve preparations as previously described [START_REF] Wetzel | A stomatin-domain protein essential for touch sensation in the mouse[END_REF][START_REF] Walcher | Specialized mechanoreceptor systems in rodent glabrous skin[END_REF][START_REF] Schwaller | USH2A is a Meissner's corpuscle protein necessary for normal vibration sensing in mice and humans[END_REF].

We made recordings from single myelinated afferents in the saphenous nerve which innervates the hairy skin of the foot or from the tibial nerve that innervates the glabrous skin of the foot [START_REF] Walcher | Specialized mechanoreceptor systems in rodent glabrous skin[END_REF][START_REF] Schwaller | USH2A is a Meissner's corpuscle protein necessary for normal vibration sensing in mice and humans[END_REF]. In control nerves all the single units (n=78) with conduction velocities in the Aβ-fiber range (>10 m/s) could be easily classified as either a rapidly-adapting mechanoreceptor or a slowly-adapting mechanoreceptor (RA-LTMR or SA-LTMRs respectively), using a set of standard quantitative mechanical stimuli. However, in the Isl1 +/Cre ::Meis2 LoxP/LoxP mice about 10 and 18% of Aβ fibers in the hairy and glabrous skin respectively could not be reliably activated by any of the quantitative mechanical stimuli used. Sensory neurons that could not be activated by our standard array of mechanical stimuli, but could still be activated by rapid manual application of force with a glass rod were classified as so called "Tap" units. Such "tap" units have been found in several mice with deficits in sensory mechanotransduction [START_REF] Wetzel | A stomatin-domain protein essential for touch sensation in the mouse[END_REF][START_REF] Ranade | Piezo2 is the major transducer of mechanical forces for touch sensation in mice[END_REF]. We made recordings from SA-LTMRs from both glabrous and hairy skin, but decided to pool the data as there was an insufficient sample size from either skin area alone. We reasoned that electrophysiological recordings would pick up primarily receptors that had successfully innervated Merkel cells and miss those fibers that had failed to innervate end-organs and would likely not be activated by mechanical stimuli. In this mixed sample of SA-LTMRs the mean vibration threshold was significantly elevated in Isl1 Cre/+ ::Meis2 LoxP/LoxP mice, but it was clear that many fibers in this sample had mechanical thresholds similar to those in the wild type (Figure 3F). The response of the same SA-LTMRs to a 25 Hz sinusoidal stimulus was unchanged in Isl1 Cre/+ ::Meis2 LoxP/LoxP mice compared to controls (Figure 3F). The response of these fibers to ramp stimuli of increasing velocities or to increasing amplitudes of ramp and hold stimuli were also not significantly different in mutant mice compared to controls (Figure 3F). Finally, consistent with the lack of neuronal loss in Isl1 Cre/+ ::Meis2 LoxP/LoxP , the number of recorded fibers were identical in WT and Isl1 Cre/+ ::Meis2 LoxP/LoxP (Figure 4 Supplementary 1E). Altogether, these data indicate that Meis2 is necessary for Merkel cells innervation in the glabrous, but not in the hairy skin. There was some indication that amongst SA-LTMRs with mechanosensitivity there was a light loss of sensitivity that could be associated with poor innervation of Merkel cells in the glabrous skin.

Meis2 gene inactivation compromises structure and function of RA-LTMRs.

In the glabrous skin, Meissner corpuscles are located in the dermal papillae and are innervated by Rapidly-Adapting type LTMR (RA-LTMR) that detect small-amplitude skin vibrations <80 Hz.

Histological analysis of the glabrous skin showed that Nefh-positive innervation of the Meissner corpuscles was strongly disorganized (Figure 4A, Figure 4 Supplementary video 1-2). However recordings from RA-LTMRs innervating these structures in Isl1 Cre/+ ::Meis2 LoxP/LoxP animals showed largely normal physiological properties (Figure 4B). Thus RA-LTMRs recorded from

Isl1 Cre/+ ::Meis2 LoxP/LoxP displayed normal vibration sensitivity in terms of absolute threshold and their ability to follow 25 sinusoids. There was a tendency for RA-LTMRs in Isl1 Cre/+ ::Meis2 LoxP/LoxP mutant mice to fire fewer action potentials to sinusoids and to the ramp phase of a series 2 second duration ramp and hold stimuli, but these differences were not statistically significant (Figure 4B).

In the hairy skin, RA-LTMRs form longitudinal lanceolate endings parallel to the hair shaft of guard and awl/auchene hairs and respond to hair deflection only during hair movement, but not during maintained displacement [START_REF] Lechner | Hairy sensation[END_REF]. Similar to Meissner corpuscles, they are tuned to frequencies between 10 and 50 Hz [START_REF] Schwaller | USH2A is a Meissner's corpuscle protein necessary for normal vibration sensing in mice and humans[END_REF]. Whole mount analysis of Nefh-positive fibers in the adult back skin showed an overall decrease in the innervation density of hairs in Isl1 Cre/+ ::Meis2 LoxP/LoxP animals compared to WT (Figure 4C). Our analysis revealed significant decreases in both the number of plexus branch points and in the number of innervated hair follicles (Figure 4D).

Consistent with the hypo innervation of hair follicle in Isl1 Cre/+ ::Meis2 LoxP/LoxP , we observed robust deficits in the mechanosensitivity of RA-LTMRs in the hairy skin (Figure 4E). Thus we needed sinusoids of significantly larger amplitudes to evoke the first (threshold) spike in RA-LTMRs. We measured the total number of spikes evoked by a sinusoid stimulus (25 Hz) of gradually increasing amplitude and here too RA-LTMRs in the in Isl1 +/Cre ::Meis2 LoxP/LoxP mutant mice fired considerably less than in controls. This finding was confirmed using a series of vibration steps of increasing amplitudes again demonstrating decreased firing to 25 Hz vibration stimuli (Figure 4E). and another Ntrk3/Nefh-positive that do not depend on Ntrk1 signaling during development [START_REF] Niu | Dual Innervation of neonatal Merkel cells in mouse touch domes[END_REF]. However, the functional significance of these different innervations is unknown. Denervation in rat also pointed at differences between Merkel cells of the glabrous and the hairy skin. Following denervation, Merkel cells of the touch dome almost fully disappear, whereas in the footpad, Merkel cells developed normally [START_REF] Mills | The neural dependency of Merkel cell development in the rat: The touch domes and foot pads contrasted[END_REF]. Thus to expression in Aδ-LTMR, peptidergic and non-peptidergic neurons. Surprisingly, although Meis2 is expressed in proprioceptive neurons [START_REF] Usoskin | Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencing[END_REF][START_REF] Zheng | Deep Sequencing of Somatosensory Neurons Reveals Molecular Determinants of Intrinsic Physiological Properties[END_REF][START_REF] Shin | Intrinsic control of neuronal diversity and synaptic specificity in a proprioceptive circuit[END_REF], their function appeared not to be dramatically affected as seen by normal gait behavior in catwalk analysis. This is in agreement with studies in which HoxC8 inactivation, a classical Meis TF cofactor expressed by proprioceptive neurons from E11.5 to postnatal stages, did not affect proprioceptive neurons early molecular identity or survival [START_REF] Shin | Intrinsic control of neuronal diversity and synaptic specificity in a proprioceptive circuit[END_REF]. In the HoxC8 mutant, proprioceptive neurons normally innervate their peripheral target and extend axons in the central nervous system, but instead exhibit ectopic synapses.

D-Hair mechanoreceptors or

our
Understanding the transcriptional programs controlling each step in the generation of a given fully differentiated and specified neuron is an extensive research field in developmental neurobiology.

Basic studies in model organisms led to a functional classification of TFs. Proneural TFs such as

Neurogenins control the expression of generic pan-neuronal genes and are able to reprogram nearly any cell type into immature neurons [START_REF] Guillemot | Beyond proneural: emerging functions and regulations of proneural proteins[END_REF][START_REF] Baker | All in the family: Proneural bHLH genes and neuronal diversity[END_REF]. Terminal selectors are TFs mastering the initiation and maintenance of terminal identity programs through direct regulation of neuron type-specific effector genes critical for neuronal identity and function such as genes involved in neurotransmitters synthesis and transport, ion channels, receptors, synaptic connectivity or neuropeptide content [START_REF] Hobert | Terminal Selectors of Neuronal Identity[END_REF][START_REF] Hobert | Neuronal identity control by terminal selectors in worms, flies, and chordates[END_REF]. The proneural function of Ngn1 and Ngn2 genes in neural crest cells, the precursors of DRG sensory neurons is well demonstrated [START_REF] Ma | NEUROGENIN1 and NEUROGENIN2 control two distinct waves of neurogenesis in developing dorsal root ganglia[END_REF][START_REF] Marmigère | Specification and connectivity of neuronal subtypes in the sensory lineage[END_REF], and several terminal selector genes shaping the different DRG sensory subpopulations have also been clearly identified including for cutaneous LTMRs [START_REF] Marmigère | Specification and connectivity of neuronal subtypes in the sensory lineage[END_REF][START_REF] Lallemend | Molecular interactions underlying the specification of sensory neurons[END_REF][START_REF] Olson | The specification and wiring of mammalian cutaneous low-threshold mechanoreceptors[END_REF]. Maf, Runx3, Shox2, ER81 and Pea3 are part of this regulatory transcriptional network regulating cutaneous LTMR neurons diversification through intermingled crossed activation and/or repression of subclasses specific effector genes.

In Shox2 mutant mice, lanceolate endings in the hairy skin and Nefh-positive innervation of Merkel cells both in the glabrous skin and in touch domes are defective [START_REF] Abdo | Dependence on the transcription factor Shox2 for specification of sensory neurons conveying discriminative touch[END_REF]. Meissner corpuscles are also poorly innervated in these mice [START_REF] Abdo | Dependence on the transcription factor Shox2 for specification of sensory neurons conveying discriminative touch[END_REF]. Maf is important for Aβ RA- LTMR both in the hairy and the glabrous skin, and in Maf mutant mice, the number of Meissner corpuscles are greatly reduced and circumferential and lanceolate endings are affected while Merkel cells neurites complexes remain intact although Aβ-SAI LTMR have higher mechanical thresholds [START_REF] Wende | The transcription factor c-Maf controls touch receptor development and function[END_REF]. In Runx3 mutants, Ntrk3-positive proprioceptive neurons failed to extend peripherally and die [START_REF] Inoue | Runx3 controls the axonal projection of proprioceptive dorsal root ganglion neurons[END_REF], but Runx3 expression by some LTMR suggests that it play additional functions in these populations [START_REF] Yoshikawa | Coexpression of Runx1 and Runx3 in mechanoreceptive dorsal root ganglion neurons[END_REF]. In addition to their functions as terminal selectors for the proprioceptive subclasses [START_REF] Lin | Functionally related motor neuron pool and muscle sensory afferent subtypes defined by coordinate ETS gene expression[END_REF] ranging from regulation of neuroblasts cell-cycle exit, to cell-fate decision, neurogenesis, neuronal specification, neurites outgrowth and mature neurons maintenance. In DRG LTMRs, Meis2 fulfils some but not all of the criteria defining terminal selector genes. Its inactivation in neural crest cells does not affect sensory neurons generation nor pan-neuronal features, clearly excluding it from a proneural TF function. Although Meis2 expression is continuously maintained in defined sensory neurons subtypes starting from early post-mitotic neurons throughout life, its expression is not restricted to a unique neuronal identity and its early or late inactivation in either sensory neurons progenitors or post-mitotic neurons does not influence neuronal subtypes identity nor survival as seen by the normal numbers of Ntrk2, Ntrk3 or c-Ret positive neurons. It is however possible that in LTMRs, Meis2 regulate other type of terminal effector genes such as genes participating in neurotransmitters machinery specification and/or recognition, establishment and maintenance of physical interactions between LTMRs and their peripheral target.

To conclude, together with studies on Meis2 function in the SVZ where Meis2 is necessary to maintain the neurogenic effect of Pax6 in neural progenitors and is later expressed in their mature progenies [START_REF] Agoston | Meis2 is a Pax6 co-factor in neurogenesis and dopaminergic periglomerular fate specification in the adult olfactory bulb[END_REF], our results raise the possibility that this TF sets up a lineage specific platform on which various specific co-factors in turn participate in different steps of the neuronal differentiation program. Sensory neurons subtypes specification is the consequence of intermingled repressive and activator activities of identified terminal selector TFs that activate subclasses specific effector genes while repressing those specific to other subclasses. Thus, in this scenario, Meis2 could be a central hub dynamically interacting with different TFs belonging to distinct functional classes including terminal selector genes. Figure 1 Supplementary Table 1: Isl1 +/Cre ::Meis2 LoxP/LoxP adult mice exhibit normal locomotion. Table recapitulating different Catwalk two-paw analysis parameters in 3 months old female mice. Several recordings were performed for each mice. Only sequence when mice showed a constant and straight locomotion with an average speed between 25 and 55 cm/sec were selected for analysis. Student-ttest analysis showed no significant differences for any of the parameters. For retro-tracing experiments, newborn pups were anesthetized on ice and cholera toxin B coupled to Cholera Toxin Subunit B conjugated with Alexa Fluo 488 (Thermofisher) was injected using a glass micropipette in several points of the glabrous and hairy forepaw. Mice were sacrificed 7 days after injection and L4 to L6 DRGs were collected for analysis.

For behavioral assays, skin-nerve preparation and electrophysiological recording, sex-matched 12weeks old mutant and WT littermates mice were used.

Chick: Fertilized eggs were incubated at 37°C in a humidified incubator. For limb ablation experiments, eggs were opened on the third day of incubation (embryonic day 3, stage 17/18) [START_REF] Hamburger | A series of normal stages in the development of the chick embryo[END_REF] and the right hind limb bud was surgically removed as previously reported [START_REF] Oakley | Muscle sensory neurons require neurotrophin-3 from peripheral tissues during the period of normal cell death[END_REF]. Eggs were closed with tape and further grown in the incubator for 4 (HH27) or 7 (HH36) additional days before collection.

Tissue preparation.

Mouse and chick embryos were collected at different stages, fixed in 4% paraformaldehyde/PBS overnight at 4°C and incubated overnight at 4°C for cryopreservation in increasing sucrose/PBS solutions (10 to 30% sucrose). After snap freezing in TissueTek, embryos were sectioned at 14-µm thickness and stored at -20°C until use.

Cloning of mouse and chick Meis2 and probes preparation.

For preparation of digoxigenin-and fluorescein-labeled probes, RNA from whole mouse or chick embryos was extracted using Absolutely RNATM Nanoprep kit (Stratagene) following manufacturer's instruction. Reverse transcription (RT) was carried out 10 minutes at 65°C followed by 1 hour at 42°C

and 15 minutes at 70°C in 20 µl reactions containing 0.5 mM dNTP each, 10 mM DTT, 0.5 μg oligod(T)15 (Promega) and 200 U of Super Script II RT (Gibco BRL Life Technologies). A 1206 bp long and a 1201 bp long Meis2 fragments were amplified from mouse and gallus cDNA respectively using the following primers: mMeis2 forward: 5'-ATGGCGCAAAGGTACGATGAGCT-3'; mMeis2 reverse: 5'-TTACATATAGTGCCACTGCCCATC-3'; gMeis2 forward: 5'-ATGGCGCAAAGGTACGATGAG-3'; gMeis2 reverse: 5'-TTACATGTAGTGCCATTGCCCAT-3'. PCR were conducted in 50 µl reactions containing 10%

RT product, 200 μM each dNTP, 10 pmol of each primer (MWG-Biotech AG), 3 mM MgCl2, 6% DMSO and 2.5 U of Herculase hotstart DNA polymerase (Stratagene). cDNA was denatured 10 minutes at 98°C and amplified for 35 cycles in a three steps program as following: 1 minute denaturation at 98°C, 1 minute annealing at annealing temperature and then 1.5 minutes polymerization at 72°C.

PCR products were separated on 2% agarose gels containing ethidium bromide. Bands at the expected size were excised, DNA was extracted, and the fragment was cloned into pCR4Blunt-TOPO vector (Invitrogen) and confirmed by sequencing. Other probes used for ISH have been described

elsewhere [START_REF] Bouilloux | Loss of the transcription factor Meis1 prevents sympathetic neurons target-field innervation and increases susceptibility to sudden cardiac death[END_REF].

In situ hybridization (ISH). Immunochemistry.

Immunochemistry was performed as previously described [START_REF] Bouilloux | Loss of the transcription factor Meis1 prevents sympathetic neurons target-field innervation and increases susceptibility to sudden cardiac death[END_REF]. In situ hybridized sections or new sections were washed 3x10 min with PBS, blocked with 4% normal goat serum, 1% bovine serum albumin and 0.1% Triton X100 in PBS and incubated overnight at 4°C with primary antibodies. After washing 3x10 min with PBS incubation occurred for 2-4 hours with secondary species and isotype-specific fluorescent antibodies (Alexa Fluor Secondary Antibodies, Molecular Probes). After repeated washing with PBS, slides were mounted in Glycerol/PBS (9/1) or Mowiol.

Picture were taken using a confocal microscope (Leica SP5-SMD, Germany). Confocal images are presented as maximal projections.

The following antibodies were used for immunochemistry: mouse-anti-islet1 39.4D used for mouse Whole-mount immunohistochemistry.

Whole-mount immunohistochemistry of adult mice back hairy skin was performed a described elsewhere [START_REF] Chang | Flat mount imaging of mouse skin and its application to the analysis of hair follicle patterning and sensory axon morphology[END_REF]. Briefly, mice were euthanized by CO 2 asphyxiation. The back skin was shaved and cleaned with commercial hair remover. The back skin was removed, carefully handled with curved forceps and fixated in 4% PFA at 4°C for 20'. The tissue was then washed with PBS containing 0.3% Triton X-100 (PBST) every 30 minutes for 3 to 5 hours and kept overnight in the washing solution. The next day, the skin was incubated for 5 days with primary antibodies diluted in PBST containing 5% donkey serum and 20% DMSO. The skin was washed the following day 8 to 10 times over a day before being incubated with secondary antibodies diluted in PBST containing 5% donkey serum and 20% DMSO. The skin was then washed every 30 minutes for 6 to 8 hours before being dehydrated in successive baths of 25%, 50%, 75% and 100% methanol. They were then incubated overnight in a 1:2 mixture of benzyl alcohol and benzyl benzoate before being mounted and sealed into chambers filled with the same medium.

Hematoxylin-Eosin staining.

As previously described [START_REF] Bouilloux | Loss of the transcription factor Meis1 prevents sympathetic neurons target-field innervation and increases susceptibility to sudden cardiac death[END_REF], air dried frozen sections were washed in water then stained with hematoxylin for 1^min at room temperature and washed extensively with water. After dehydration in PBS/alcohol (70%), slides were stained with eosin for 30 sec at room temperature.

After serial wash in water, sections were dehydrated in PBS solutions with increasing alcohol concentration (50%, 75%, 95%, and 100%), mounted and observed with a microscope (Leica DMRB, Germany).

Mouse skin-nerve preparation and sensory afferent recordings

Cutaneous sensory fiber recordings were performed using the ex vivo skin-nerve preparation as previously described [START_REF] Schwaller | USH2A is a Meissner's corpuscle protein necessary for normal vibration sensing in mice and humans[END_REF]. Mice were euthanized by CO2 inhalation for 2-4^min followed by cervical dislocation. Three different preparations were performed in separate experiments using different paw regions: the saphenous nerve innervating the hairy hind paw skin;

the tibial nerve innervating the glabrous hind paw skin; and the medial and ulnar nerves innervating the forepaw glabrous skin. In all preparations, the hairy skin of the limb was shaved and the skin and nerve were dissected free and transferred to the recording chamber, where muscle, bone and tendon tissues were removed from the skin to improve recording quality. The recording chamber was perfused with a 32^°C synthetic interstitial fluid: 123^mM NaCl, 3.5^mM KCl, 0.7^mM MgSO4, 1.7^mM NaH2PO4, 2.0^mM CaCl2, 9.5^mM sodium gluconate, 5.5^mM glucose, 7.5^mM sucrose and 10^mM 4-(2-hydroxyethyl)-1-piperazine-ethanesulfonic acid (Hepes), pH 7.4. The skin was pinned out and stretched, such that the outside of the skin could be stimulated using stimulator probes. The peripheral nerve was fed through to an adjacent chamber in mineral oil, where fine filaments were teased from the nerve and placed on a silver-wire recording electrode.

The receptive fields of individual mechanoreceptors were identified by mechanically probing the surface of the skin with a blunt glass rod or blunt forceps. Analog output from a Neurolog amplifier was filtered and digitized using the Powerlab 4/30 system and Labchart 7.1 software (AD instruments). Spike-histogram extension for Labchart 7.1 was used to sort spikes of individual units.

Electrical stimuli (1^Hz, square pulses of 50-500^ms) were delivered to single-unit receptive fields to measure conduction velocity and enable classification as C-fibers (velocity <1.2^m^s-1), Aδ-fibers

(1.2-10^m^s-1) or Aβ-fibers (>10^m^s-1). Mechanical stimulation of the receptive fields of neurons were performed using a piezo actuator (Physik Instrumente, catalog no. P-841.60) and a doubleended Nanomotor (Kleindiek Nanotechnik, catalog no. MM-NM3108) connected to a force measurement device (Kleindiek Nanotechnik, catalog no. PL-FMS-LS). Calibrated force measurements were acquired simultaneously using the Powerlab system and Labchart software during the experiment.

As different fiber types have different stimulus-tuning properties, different mechanical stimuli protocols were used based on the unit type. Low-threshold Aβ-fibers (RAMs and SAMs) and Aδ-fiber D-hairs were stimulated with the piezo actuator with three vibration stimuli (5^Hz, 25^Hz and 50^Hz, distortions introduced by the in-series force sensor precluded using frequencies >50^Hz) with increasing amplitude over six steps (peak-to-peak amplitudes of ~6-65^mN; 20 cycles per step), and a dynamic stimulus sequence with four ramp-and-hold waveforms with varying probe deflection velocities (3^s duration; 0.075, 0.15, 0.45 and 1.5^mm^s-1; average amplitude 100^mN). Aβ-fiber SAMs and RAMs were classified by the presence or absence of firing during the static phase of a ramp-and-hold stimulus, respectively, as previously described. Single units were additionally stimulated with a series of five static mechanical stimuli with ramp-and-hold waveforms of increasing amplitude (3^s duration; ranging from ~10^mN to 260^mN). Low-threshold SAMs, high-threshold

Aδ-fibers and C-fibers were also stimulated using the nanomotor with five ramp-and-hold stimuli with increasing amplitudes 20.

Behavioral assays.

Von Frey paw withdrawal test: Mice were placed on an elevated wire mesh grid into PVC chambers.

Before the test, mice were habituated to the device for one hour for two consecutive days. On the day test, mice were placed in the chamber one hour before Von Frey filaments application. The test was performed as previously described [START_REF] Neubarth | Meissner corpuscles and their spatially intermingled afferents underlie gentle touch perception[END_REF]. During the test, withdrawal response following Von Frey filament application on the palm of the left hind paw was measured. Starting with the lowest force, each filament ranging from 0.008 g to 1.4 g was applied ten times in a row with a break of 30 seconds following the fifth application. The number of paw withdrawals for each filament was counted.

Hot Plate test: Before starting the test, mice were habituated to the experimentation room for at least five minutes. Mice were individually placed on the hot plate set up at 53°C and removed at the first signs of aversive behavior (paw licking or shaking). The time to this first stimulus was recorded. A 30 sec cut off was applied to avoid skin damages. After 5 min recovery in their home cage, the test was repeated three times for each mouse and averaged. Data are shown as the average of these three measurements.

Sticky tape test: A two cm 2

of laboratory tape was placed on the upper back skin of mice just before they were placed on an elevated wire mesh grid into PVC chambers. The number of tape-directed reactions were then counted during 5 min. Considered responses were body shaking like a "wetdog", hindlimb scratching directed to the tape, trying to reach the tape with the snout and grooming of the neck with forepaws.

Dynamic touch test: Mice were placed in the same conditions as described above for the Von Frey paw withdrawal test. Sensitivity to dynamic touch was performed by stroking hind paws with a tapered cotton-swab in a heel to toe direction. The stimulation was repeated 10 times by alternating left and right hind paws and the number of paw withdrawals was counted.

Gait analysis: Gait was analyzed using the Catwalk™ system (Noldus Information Technology, Netherlands) in a dark room with minimized light emission from the computer screen. Mice were allowed to voluntarily cross a 100-cm-long, 5-cm-wide walkway with a glass platform illuminated by green fluorescent light. An illuminated image is produced when a mouse paw touches the glass floor through dispersion of the green light, and footprints were captured by a high-speed camera placed under the glass floor. Data were analyzed using the CatWalk™ XT 10.1 software. For each mouse, several recordings were performed until at least 3 runs met criteria defined by a minimum of three consecutive complete step cycles of all four paws without stopping or hesitation and within the range of 25 to 50 cm.s -1

. Data are reported as the average of at least three runs per mouse.

  (SAI-LTMR) neurons responding to both static skin indentation and moving stimuli such as vibration. In the glabrous skin, Merkel cells form clusters in the basal layer of the epidermis, and in the hairy skin, similar clusters of Merkel cells called touch-domes are located at the bulge region of guard hairs. Histological analysis indicated that in the forepaw glabrous skin of Isl1 +/Cre ::Meis2 LoxP/LoxP adult mice, the number of Merkel cells contacted by Nefh-positive fibers was strongly decreased compared to WT (Figure3B and E). However, in contrast to the glabrous skin, Merkel cell innervation by Nefh-positive fibers appeared largely unaffected in the hairy skin of Isl1 Cre/+ ::Meis2 LoxP/LoxP mice (Figure3D and E). Whole mount analysis of CK8-positive Merkel cells in the hairy back skin of E18.5 embryos showed that the overall number of touch domes and of Merkel cells per touch dome were unchanged in mutant animals compared to WT (Figure3C).

  best knowledge, Meis2 may represent the first TF which inactivation differentially affect Merkel cells innervation in the glabrous skin and hairy touch domes.Because touch domes innervation and Aδ fibers D-hair vibration responses were unaffected in Isl1 +/Cre ::Meis2 LoxP/LoxP mice, we postulate that the innervation defects we observed in the hairy skin is supported by defects in lanceolate endings with RA-LTMR electrophysiological properties. However, the increased number of "tap" units both in the hairy and glabrous skin is compatible with wilder deficits also including Aβ-field LTMRs peripheral projections. Similarly, although the severely disorganized Meissner corpuscle architecture did not result in significant consequences on RAM fibers electrophysiological responses in the glabrous skin, it is possible that the large increase in the number of "tap units" within Aβ fibers of the glabrous skin represent Meissner corpuscles whose normal electrophysiological responses are abolished. In agreement, challenging sensory responses in the glabrous skin with either von Frey filament application or cotton swab stroking clearly showed a dramatic loss of mechanical sensitivity specifically within the range of gentle touch neurons. Recent work reported that the von Frey filaments test performed within low forces and challenging light touch sensation could distinguish Merkel cells from Meissner corpuscles dysfunctions: mice depleted from Merkel cells performed normally on this test while mice mutated at the Ntrk2 locus that display Meissner corpuscle innervation deficits showed decreased hind paw withdrawal threshold in response to filament within the 0.02-0.6 g range[START_REF] Neubarth | Meissner corpuscles and their spatially intermingled afferents underlie gentle touch perception[END_REF]. Thus, our result in the von Frey filament test possibly involves Meissner corpuscles functioning. Finally, the unaltered D-hair fibers electrophysiological responses, the normal noxious responses in the hot plate setting and the normal responses to high forces Von Frey filaments are consistent with the absence of Meis2

Figure 1 :

 1 Figure legends:

Figure 2 :

 2 Figure 2: Meis2 is necessary for LTM neurons peripheral projections but dispensable for their survival and specification. A) Representatives images showing IF for Ntrk2 or Ntrk3 (green) with Pvalb or Maf (red) in P0 WT and Wnt1 Cre ::Meis2 LoxP/LoxP DRGs. Box plots showing that the number of

Figure 3 :

 3 Figure 3: Meis2 gene inactivation compromised Merkel cells innervation in the glabrous skin and increased SAM vibration threshold. A) Percentage of tap units among Aβ fibers in the hairy and glabrous skin detected by electrophysiological recording in the nerve-skin preparation. Note that Tap-units are only present in both the hairy and glabrous skin of adult Isl1 +/Cre ::Meis2 LoxP/LoxP mice but

Figure 4 :

 4 Figure 4: Meis2 gene inactivation affects Meissner corpuscles morphology, hair follicle innervation and RAM fibers electrophysiological responses in the hairy but not the glabrous skin. A) Representative images showing S100β + Meissner corpuscles (red) and their innervation by Nefh +

Figure 1

 1 Figure 1 supplementary 1: Meis2 mRNA expression in LTM neurons of mouse DRG. A) Double ISH for Meis2 (blue) and Ntrk2, Ntrk3 or c-Ret (red) showed that Meis2 mRNA partly colocalises with mRNA for Ntrk2, Ntrk3 and c-Ret in E14.5 mouse embryos. Arrows point at double positive neurons. Arrowheads point at only Meis2 + neurons. Stars indicate Meis2 -/Ntrk2 + or Ntrk3 + neurons. Note that all large c-Ret-positive neurons (pseudo-color in green) express Meis2 mRNA. Scale bar=25µm. B) Double ISH for Meis2 (blue) and Ntrk2 or Ntrk3 (red) on E18.5 embryos. Arrows point at double positive neurons. Arrowheads point at only Meis2 + neurons. Stars indicate Meis2 -/Ntrk2 + or Ntrk3 + neurons. Scale bar=25µm. C) Combined ISH for Meis2 mRNA with IF for Ntrk1 (red) and Islet1 (green) showed that Meis2 is expressed by Islet1-positive post-mitotic neurons and is mostly excluded from the Ntrk1-positive subpopulation of DRG sensory neurons at E11.5 and E18.5. Arrowheads point at Meis2 + /Ntrk1 neurons, arrows point at Meis2 + /Ntrk1 + neurons. Note that the level of Meis2 mRNA

Figure 1

 1 Figure 1 Supplementary 2: Meis2 is expressed in a subset of chick ventro-lateral DRG sensory neurons during embryogenesis. A) Developmental expression of Meis2 visualized by ISH in chick DRG at HH24, HH29 and HH36. Dashed lines delineate the DRGs split into the ventro-lateral (VL) and dorso-medial (DM) parts. At early stages after DRG condensation, Meis2 showed a broad expression in chick DRGs neurons. As differentiation progresses (HH29 and HH36), Meis2 expression became progressively restricted to the VL population of sensory neurons which represents the Ntrk2 + and

Figure 1

 1 Figure 1 Supplementary 3: Mice with a conditional deletion of Meis2 gene in neural crest derivatives (Wnt1 Cre ) exhibited cleft-palate and died at birth. A) Targeting vector used for the generation of a conditional knockout mouse strain for Meis2 (Meis2 LoxP/LoxP ). The Exon 8 (first coding exon of the DNA binding Homeodomain) is flanked by LoxP sites allowing deletion of Meis2 DNA binding domain. B) Representative images showing that new-born Wnt1 Cre ::Meis2 LoxP/LoxP pups mutant pups are smaller at birth compared to WT littermates, unable to stand on their legs and exhibit a cleft palate phenotype. Representatives images showing coronal sections of new-born WT

Figure 2

 2 Figure 2 Supplementary 1: Meis2 gene inactivation does not affect phospho-Creb expression. Representative images showing phosphor-Creb (pCreb in green), Ntrk2 or Ntrk3 (red) and Isl1 (blue) expression in WT and Isl 1+/Cre ::Meis2 LoxP/LoxP E16.5 DRG embryos.

Figure 3

 3 Figure 3 Supplementary 1: WT and Isl1 +/Cre mice exhibit similar electrophysiological responses. Graphs show the same results than in Figure 3 and 4 but including Isl1 +/Cre mice (n=5 from 2 mice) a s

Figure 4

 4 Figure 4 Supplementary video 1: Meissner corpuscle in WT. 3D visualization of Meissner corpuscles in adult WT glabrous skin visualized by IF against Nefh and S100β.

Figure 4

 4 Figure 4 Supplementary video 2: Meissner corpuscles in Isl1 +/Cre ::Meis2 LoxP/LoxP . 3D visualization of

Figure 4

 4 Figure 4 Supplementary 1: Normal electrophysiological responses of D-Hair following Meis2 gene inactivation. A) D-Hair of the hairy skin exhibited similar vibration threshold and firing activity to a 25-Hz vibration in WT (n=8 from 3 mice) and Isl1 +/Cre ::Meis2 LoxP/LoxP mice (n = 10 from 3 mice). B) D-Hair showed similar responses to a ramp of 50 Hz vibration with increasing amplitude and to a ramp stimuli in WT and Isl1 +/Cre ::Meis2 LoxP/LoxP animals. Traces showed the type of stimulation and red squares indicate the time frame during when the number of spikes was calculated.

  Before hybridization, slides were air dried for 2-3 hours at room temperature. Plasmids containing probes were used to synthesize digoxigenin-labeled or fluorescein-labeled antisense riboprobes according to supplier's instructions (Roche) and purified by LiCl precipitation. Sections were hybridized overnight at 70°C with a solution containing 0.19 M NaCl, 10 mM Tris (pH 7.2), 5 mM NaH2PO4*2H2O/Na2HPO4 (pH 6.8), 50 mM EDTA, 50% formamide, 10% dextran sulphate, 1 mg/ml yeast tRNA, 1XDenhardt solution and 100 to 200 ng/ml of probe. Sections were then washed four times for 20 minutes at 65°C in 0.4X SSC pH 7.5, 50% formamide, 0.1% Tween 20 and three times for 20 minutes at room temperature in 0.1 M maleic acid, 0.15 M NaCl and 0.1% Tween 20 (pH 7.5).Sections were blocked 1 hour at room temperature in presence of 20% goat serum and 2% blocking agent (Roche) prior to incubation overnight with AP-conjugated anti-DIG-Fab-fragments (Roche, 1:2000). After extensive washing, hybridized riboprobes were revealed by performing a NBT/BCIP reaction in 0.1 M Tris HCl pH 9.5, 100 mM NaCl, 50 mM MgCl2 and 0.1% Tween 20.For double in situ hybridization, the procedure was the same except that hybridization was conducted by incubation with 100-200 ng/ml of one digoxigenin-labeled probe and 100-200 ng/ml of -conjugated anti-fluorescein-Fab-fragment (Roche, 1:2000) and further incubation with Fast Red tablets in 0.1 M Tris HCl pH 8.5, 100 mM NaCl, 50 mM MgCl2 and 0.1% Tween 20.Pictures of fluorescein alone were taken after mounting in glycerol/PBS (1:9). To reveal the digoxigenin-labeled probe, sections were unmounted, washed extensively in PBS and Alkalin Phosphatase was inhibited by incubation in a solution of 0.1M glycin pH2.2, 0.2% Tween 20 for 30 minutes at room temperature. After extensive washing in PBS, digoxigenin-labeled-probe wa s revealed as described using the AP-conjugated anti-DIG-Fab-fragments (Roche, 1:2000) and the NBT/BCIP reaction. Sections were mounted again in glycerol/PBS (1:9) and pictures of both fluorescein and digoxigenin were taken. For removing the Fast Red staining, sections were unmounted again, washed extensively in PBS and incubated in increasing solutions of ethanol/PBS solutions (20-100% ethanol). After extensive washing in PBS, sections were mounted in glycerol/PBS (1:9) and pictures of the digoxigenin staining alone. Wide field microscopy (Leica DMRB, Germany) was only used for ISH and ISH combined with immunochemistry.

  and chick (diluted 1:100, Developmental Studies Hybridoma Bank); TrkB (1/2000; R and D Systems); TrkC (1/1000; R and D Systems); CGRP (1/500; Sigma-Aldrich); TrkA (1/500; Millipore); rabbit antiparvalbumin antibody (1:500, Swant); guinea pig anti-calcitonin gene-related peptide (CGRP) antibody (1:500, Peninsula Laboratories); Nfh (Sigma N4142, rabbit 1:1000), c-maf (generous gift of C. Birchmeier, MDC, 1/1000), TrpV1 (Sigma V2764, rabbit 1:1000), S100β (Sigma S2532, mouse 1:1000), Phospho-CREB (Cell Signaling, 87G3, 1/200, Germany); Goat anti-c Ret (R and D Systems, Cat#AF482, 1/100) and Meis2 (Sigma Aldrich, WH0004212M1 or Abcam ab244267, 1/500). The chick TrkB and C antibodies were a generous gift from LF Reichardt, UCSF and have been previously reported.

  Figure 1

  subpopulation (Figure1Supplementary 2C). In adult mouse DRG, comparison of Meis2 mRNA to Ntrk1, Calca and TrpV1 immunostaining confirmed that Meis2-expressing neurons are largely excluded from the nociceptive and thermo-sensitive populations of DRG neurons. Instead, a large proportion of Meis2-positive neurons co-expressed Nefh, a marker for large myelinated neurons including LTMR and proprioceptors, and Pvalb, a specific

	marker for proprioceptors (Figure 1 Supplementary 1D). Finally, Meis2 expression in LTMRs
	projecting to the skin was confirmed by retrograde-tracing experiments using Cholera toxin B subunit
	(CTB) coupled with a fluorochrome injected into hind paw pads of P5 newborn mice. Analyses of CTB
	expression in lumbar DRG three days later at P8 showed that many retrogradely labelled sensory
	neurons were also immuno-positive for Meis2, Maf, Ntrk2 and Ntrk3 (Figure 1B).
	Altogether, our results on Meis2 co-localization with Nfeh, NtrK2, Ntrk3, Ret, Pvalb and Maf at

different embryonic and postnatal stages are consistent with previous report on restricted Meis2 expression to the Aβ-field, Aβ-SA1 and Aβ-RA subclasses of LTMR neurons and proprioceptive neurons

  The function of the Meis family of TFs in post-mitotic neurons has only been marginally addressed[START_REF] Agoston | Meis2 is a Pax6 co-factor in neurogenesis and dopaminergic periglomerular fate specification in the adult olfactory bulb[END_REF][START_REF] Jakovcevski | Neuronal Kmt2a/Mll1 histone methyltransferase is essential for prefrontal synaptic plasticity and working memory[END_REF][START_REF] Bouilloux | Loss of the transcription factor Meis1 prevents sympathetic neurons target-field innervation and increases susceptibility to sudden cardiac death[END_REF]. In agreement with previous reports on Meis2 expression in DRG sensory neurons[START_REF] Usoskin | Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencing[END_REF][START_REF] Zheng | Deep Sequencing of Somatosensory Neurons Reveals Molecular Determinants of Intrinsic Physiological Properties[END_REF][START_REF] Sharma | The emergence of transcriptional identity in somatosensory neurons[END_REF][START_REF] Shin | Intrinsic control of neuronal diversity and synaptic specificity in a proprioceptive circuit[END_REF], we showed that Meis2 is selectively expressed by This difference support previous work suggesting that the primary afferents innervating Merkel cells in the glabrous and the hairy skin maybe different[START_REF] Niu | Dual Innervation of neonatal Merkel cells in mouse touch domes[END_REF][START_REF] Olson | The specification and wiring of mammalian cutaneous low-threshold mechanoreceptors[END_REF]. Whereas Merkel cells of the glabrous skin are exclusively contacted by large Ntrk3/Nefh-positive Aβ afferents, neonatal mouse touch domes receive innervation of two types of neuronal populations, a Ret/Ntrk1-positive one that depend on Ntrk1 for survival and innervation,

	Discussion:
	subpopulations	of	early	post-mitotic	cutaneous	LTMR	and	proprioceptive	neurons	during
	development of both mouse and chick, highlighting the conserved Meis2 expression across
	vertebrates species in those neurons. We unambiguously demonstrate that Meis2 differentially
	regulates target-field innervation and function of post-mitotic LTMR neurons by using different
	recombination strategies: Meis2 gene inactivation prior to birth in migratory and dividing neural
	crest cells only (Wnt1 Cre ) or in early Isl1-positive post-mitotic sensory neurons. In addition, Meis2
	and Isl1 +/Cre ::Meis2 LoxP/LoxP The morphological and functional phenotypes we report following specific Meis2 gene inactivation in mice (Supplementary 1).
	The functional deficits in RA-LTMRs correlate well with the defects in LTMR cutaneous projections we all post-mitotic sensory neurons are consistent with Meis2 expression pattern, and ultimately, both
	mutant mice. electrophysiological observed in Isl1 +/Cre ::Meis2 LoxP/LoxP defective morphological and	responses	resumed	in	specifically	impaired
	behavioral responses to light touch mechanical stimuli. In these mutants, the decreased innervation
	of Merkel cells in the glabrous skin and the decreased sensitivity in SA-LTMR electrophysiological

Aδ-LTMRs are the most sensitive skin mechanoreceptors with very large receptive fields

[START_REF] Shin | A T-type calcium channel required for normal function of a mammalian mechanoreceptor[END_REF][START_REF] Li | The functional organization of cutaneous low-threshold mechanosensory neurons[END_REF][START_REF] Walcher | Specialized mechanoreceptor systems in rodent glabrous skin[END_REF]

. They form lanceolate endings, are thinly myelinated and are activated by movement of the smaller zigzag hairs

[START_REF] Lechner | Hairy sensation[END_REF]

. Consistent with the lack of Meis2 expression in this population reported by singlecell RNA-seq databases, Aδ fibers D-hair in the hairy skin showed similar vibration responses in WT inactivation in post-mitotic sensory neurons differentially impaired LTMRs projections and functions without affecting their survival and molecular identity.

Combined single cell RNAseq (scRNAseq) analysis and genetic tracing have accurately resolved Meis2 expression

[START_REF] Usoskin | Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencing[END_REF][START_REF] Zheng | Deep Sequencing of Somatosensory Neurons Reveals Molecular Determinants of Intrinsic Physiological Properties[END_REF][START_REF] Sharma | The emergence of transcriptional identity in somatosensory neurons[END_REF]

. A scRNAseq unbiased classification of adult mouse sensory neurons identified more than 11 molecularly distinct subtypes of neurons and Meis2 is expressed in four of the five clusters of Nefh-and Pvalb-positive myelinated neurons representing LTMR and proprioceptive neurons. Meis2 expression is mostly excluded from the C-LTMR population and the five other peptidergic and non-peptidergic nociceptive populations

[START_REF] Usoskin | Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencing[END_REF]

. In a combined subclasses-specific genetic labelling and RNA bulk deepsequencing strategy

[START_REF] Zheng | Deep Sequencing of Somatosensory Neurons Reveals Molecular Determinants of Intrinsic Physiological Properties[END_REF]

, the transcriptomes of five functionally defined LTMR subtypes (C-LTMRs, Aδ-LTMRs, Aβ RALTMRs, Aβ SA1-LTMRs, Aβ field-LTMRs) and proprioceptors were established. These analyses confirmed Meis2 expression in proprioceptive neurons, Aβ field- LTMR and Aβ-SA1-LTMR, and a weaker expression in Aβ-RA-LTMR. This study also confirmed the lack of Meis2 expression in C-LTMR, Aδ-LTMR, peptidergic and non-peptidergic neurons. responses to mechanical stimuli is consistent with Meis2 being expressed by Aβ-SA1-LTMR neurons. Interestingly, Meis2 gene inactivation compromise Merkel cells innervation and electrophysiological responses in the glabrous skin but not in touch domes of the hairy skin where innervation appeared unchanged.

  , ER81 and Pea3, two members of the Ets family of TF are also involved in Pacinian corpuscles formation. In ER81 mutant mice, Pacinian corpuscles are absent and in Pea3 mutant as well but to a lesser extent.Despite the strong implication of Meis TFs in many human neurological disorders[START_REF] Hammerschlag | Genome-wide association analysis of insomnia complaints identifies risk genes and genetic overlap with psychiatric and metabolic traits[END_REF][START_REF] Lane | Genome-wide association analyses of sleep disturbance traits identify new loci and highlight shared genetics with neuropsychiatric and metabolic traits[END_REF][START_REF] Jiménez-Jiménez | Genetics of restless legs syndrome: An update[END_REF][START_REF] Sarayloo | Meis1 and restless legs syndrome: A comprehensive review[END_REF], their mechanisms of action in particular in post-mitotic neurons are poorly understood. In humans, at least 17 different mutations in the Meis2 gene have been associated to neurodevelopmental delay and ASD[START_REF] Giliberti | MEIS2 gene is responsible for intellectual disability, cardiac defects and a distinct facial phenotype[END_REF][START_REF] Gangfuß | Intellectual disability associated with craniofacial dysmorphism, cleft palate, and congenital heart defect due to a de novo MEIS2 mutation: A clinical longitudinal study[END_REF], emphasizing its essential function in neuronal differentiation. The Meis2 function in late differentiation of post-mitotic peripheral sensory neurons add to the wide actions of this TF in the developing and adult nervous system in number of regions of the mouse nervous system. Its expression both in dividing neural progenitors, in immature neurons and in discrete populations of mature neurons[START_REF] Toresson | Retinoids are produced by glia in the lateral ganglionic eminence and regulate striatal neuron differentiation[END_REF][START_REF] Toresson | Expression of Meis and Pbx genes and their protein products in the developing telencephalon: Implications for regional differentiation[END_REF][START_REF] Allen | Molecular identity of olfactory bulb interneurons: Transcriptional codes of periglomerular neuron subtypes[END_REF][START_REF] Bumsted-O'brien | Expression of the homeodomain transcription factor Meis2 in the embryonic and postnatal retina[END_REF][START_REF] Jakovcevski | Neuronal Kmt2a/Mll1 histone methyltransferase is essential for prefrontal synaptic plasticity and working memory[END_REF][START_REF] Chang | Expression patterns of homeobox genes in the mouse vomeronasal organ at postnatal stages[END_REF][START_REF] Frazer | Transcriptomic and anatomic parcellation of 5-HT3A R expressing cortical interneuron subtypes revealed by single-cell RNA sequencing[END_REF][START_REF] Yan | Mouse Retinal Cell Atlas: Molecular Identification of over Sixty Amacrine Cell Types[END_REF][START_REF] Yang | Transcriptional profiling reveals the transcription factor networks regulating the survival of striatal neurons[END_REF] argues for diverse functions

Acknowledgements:

S.V. for help with ChTx experiments, A.L.B for cat walk experiments, al staff at animal house and in particular Flora for great help at the animal facility.

Contribution:

Design of experiments, data analysis and writing: GL, PC, FM Electrophysiology experiments: GL and FS Histology and behavior: SD and FM