Introduction

Tactile stimuli like brush, light pressure, or roughness engage highly specialized and diverse arrays of mechanoreceptors in both the hairy and glabrous skin (1–7). Somatosensory perception via these mechanoreceptors involves primary sensory neurons whose cell bodies reside within Dorsal Root Ganglia (DRG) and cranial sensory ganglia. These 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 (8–10).

Cutaneous mechanoreceptors or Low Threshold Mechanoreceptors (LMTRs) exhibit a variety of specialized terminal endings in the hairy and glabrous skin with strikingly unique morphologies (1–4, 6, 7, 11, 12). LTMRs projecting to the glabrous skin innervate Merkel cell complexes or Meissner corpuscles at the dermal-epidermal border. Those 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). On the other hand, 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, but not to static displacement (4). 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 (7, 13). A small number of LTMRs in the hairy skin are not activated by hair movement but show properties of rapidly-adapting mechanoreceptors (14). These were originally characterized as so-called field receptors (14, 15) and were recently shown to form circumferential endings around hair follicles (16). LTMRs tuned to high frequency vibration are called Aβ-RAII and innervate Pacinian corpuscles deep in the skin or on the bone (12). Ruffini endings that are thought to be innervated by stretch sensitive mechanoreceptors (Aβ-SAII) remain poorly characterized in mice (3). Recent advances in combining single cell transcriptomic, and deep RNA sequencing with genetic tracing have tremendously extended the classical subtypes repertoire and clustered at least 20 different subtypes of LTM neurons (5, 11, 16–18).

Cracking the transcriptional codes supporting sensory neurons identity and diversification has been the object of tremendous efforts in the last decades (5, 8, 9, 11, 17, 18). For instance, the functions of specification factors or terminal selectors, like Maf, Shox2, Runx3, Pea3 and ER81 have been functionally implicated in LTMR segregation (8, 9, 19–22). Whereas the specific function of adhesion molecules in shaping the assembly of touch circuitry is being unraveled (23), the transcriptional control of target cell innervation within the skin and of the establishment of specialized peripheral end-organ complexes is less understood. Meis2 is another TF expressed in LTMRs (11, 17, 18). Its mutation in human causes severe neurodevelopmental defects (24–26), an somatic mutations of its DNA consensus binding site are associated with neurodevelopmental defects (27). It belongs to a highly conserved homeodomain family containing three members in mammals, Meis1, Meis2 and Meis3 (28, 29), and Meis1 is necessary for target-field innervation of sympathetic peripheral neurons (30). We thus wondered if Meis2 could also be a pertinent regulator of late primary sensory neurons differentiation.

Here, we show that Meis2 regulates the innervation of specialized cutaneous end-organs important for LTMRs function. We confirmed that Meis2 expression is restricted to LTMR subclasses at late developmental stages compatible with functions in specification and/or target-field innervation. We generated mice carrying an inactive Meis2 gene in post-mitotic sensory neurons. While these animals are healthy and viable and do not exhibit any neuronal loss, they display tactile sensory defects in electro-physiological and behavioral assays. Consistent with these findings we identified major morphological alterations in LTMR end-organ structures in Meis2 null sensory neurons. Finally, transcriptomic analysis at late embryonic stages showed dysregulation of synapse and neuronal projection-related genes that underpin these functional and behavioral phenotypes.

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 and 2). 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 becomes restricted to a well-defined subpopulation in the ventro-lateral part of the DRG where LTMRs and proprioceptors are located(31) (Figure 1 Supplementary 2A). In both species, Meis2-positive cells also expressed the pan-neuronal marker Islet1, indicating that they are post-mitotic neurons. In chick, we estimated that Meis2-positive cells represented about 15% 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 Ntrk2- and Ntrk3-positive neurons confirming that Meis2-positive neurons belong to the LTMR and proprioceptive subpopulations. 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).

Figure 1:

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 co- expressed Ntrk1 at E11.5 and E18.5 (Figure 1 Supplementary 1C). In chick HH29 embryos, Meis2 expression was fully excluded from the Ntrk1 subpopulation (Figure 1 Supplementary 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 Nefh, 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 (11, 17, 18, 32). 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.

Target-derived signals are necessary to maintain Meis2 expression

The requirement for extrinsic signals provided by limb mesenchyme and muscles for proprioceptor and LTMR development has been documented (11, 21, 32–36). 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 HH29 and HH36, before and after ventro-lateral neurons are lost respectively (37–39) (Figure 1C and D; Figure 1 Supplementary 2D). In HH36 embryos, about 65% of Meis2-positive neurons were lost on the ablated side compared to the contralateral side (Figure 1C). This is consistent with a 30% loss of all sensory DRG neurons represented by the pan- neuronal marker Islet1, and the 50% and 65% loss of Ntrk2 and Ntrk3 positive VL-neurons respectively (Figure 1 Supplementary 2D). The number of Ntrk2-positive DL neurons were not significantly affected. In HH27 embryos, while 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 Meis2 inactivation 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). To validate the use of our strain, we first crossed the Meis2^LoxP/LoxP mice with the W n t^C1^re strain. This crossing efficiently inactivated Meis2 in the neural crest, and Wnt1^Cre::Meis2^LoxP/LoxPnew-born pups exhibited a cleft palate as previously reported in another conditional Meis2 mouse strain (40) (Figure 1 Supplementary 3B). They were however not viable, precluding functional and anatomical analyses at adult stages. To bypass this neural crest phenotype and to more specifically address Meis2 function in post-mitotic neurons, we crossed the Meis2^LoxP/LoxP mice with the Isl1^Cre/+ strain and focused our analysis on the Isl1^Cre/+::Meis2^LoxP/LoxP strain. Mutant pups were viable, appeared healthy and displayed a normal palate, allowing sensory behavior investigations.

We 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.008 to 1.4g to the hind paw and found the frequency of withdrawal responses to be significantly decreased in Isl1^+/Cre::Meis2^LoxP/LoxP mice compared to control WT and Isl1^+/Cre mice between 0.16 and 0.6g (Figure 1E), indicating that mutant mice are less responsive to light touch. No differences were observed between WT and Isl1^+/Cre mice for any of the stimuli applied. Behaviors evoked from stimulation of the glabrous skin were next assessed using the “cotton swab” dynamic touch assay (41). Here, responses were significantly decreased in Isl1^+/Cre::Meis2^LoxP/LoxP mice compared to control WT and Isl1^+/Cre littermates (Figure 1F). We also used the hot plate assay to assess noxious heat evoked behaviors and found no difference in response 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^LoxP/LoxP mice compared to WT and Isl1^Cre/+control mice (Figure 1H). Finally, although Meis2 and Isl1 are both expressed by spinal motor neurons and proprioceptors (42–44), 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).

Overall, these behavior analyses indicate that Meis2 gene inactivation specifically affects light touch sensation both in the glabrous and the hairy skin. The impaired behavioral response to light touch in Meis2 mutant suggests that Meis2 gene activity is necessary for the anatomical and functional maturation of LTMRs.

Meis2 is dispensable for LTMR specification and survival

To investigate whether Meis2 gene inactivation interfered with LTMR survival during embryonic development, we performed histological analysis of the Isl1^Cre/+::Meis2^LoxP/LoxP and Wnt1^Cre::Meis2^LoxP/LoxP strains (Figure 2). 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 2A) 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 2B). Finally, quantification of DRG neuron populations at P0 in Wnt1^Cre::Meis2^LoxP/LoxP mice showed similar results with no differences in the number of Ntrk2 and Ntrk3-positive neurons (Figure 2C). At this stage, phospho-Creb (pCreb) expression in Ntrk2 and Ntrk3-positive neurons was similar in WT and mutants (Figure 2 Supplementary 1), suggesting that Ntrk signaling is not affected. Altogether, these results show that Meis2 is dispensable for LTMR and proprioceptor survival and specification during embryogenesis.

Figure 2:

Meis2 is necessary for normal end-organ innervation

To better understand the molecular changes underlying tactile defects in Meis2 mutant mice, we performed RNAseq analysis on DRGs dissected from WT, Isl1^Cre/+ and Isl1^Cre/+::Meis2^LoxP/LoxP E18.5 embryos.

For all analyses, consistent with the changes measured for Meis2 and Isl1 genes (Figure 3 Supplementary 1A and B), only DEGs with a minimal fold change of 20% and a p value lower than 0.05 were considered. Analyses of the dataset (n=3; p<0.05; Figure 3; Figure 3 Supplementary 1; Figure 3 Supplementary Table 1) identified 43 differentially expressed genes (DEGs) in the WT v s Isl1^+/Cre::Meis2^LoxP/LoxP comparison, 107 DEGs in the Isl1^+/Cre v s Isl1^+/Cre::Meis2^LoxP/LoxP comparison, and 109 DEGs in the WT v s Isl1^+/Cre comparison. Among them, only 10 DEGs were found in both WT v s Isl1^+/Cre::Meis2^LoxP/LoxP and Isl1^+/Cre vs Isl1^+/Cre::Meis2^LoxP/LoxP comparisons (Figure 3A). Half of them were down- or up-regulated (Figure 3 Supplementary 1C), and 8 were found to be expressed in sensory neurons expressing Meis2 (Figure 3 Supplementary 1D). These include 3 ncRNA (A230077H06Rik, Gm20163, Gm42418), the Adhesion G Protein-Coupled Receptor G3 (Adgrg3, also known as GPR97), the Cellular Repressor Of E1A Stimulated Genes 2 (Creg2), predicted to be located in Golgi apparatus and endoplasmic reticulum, the Tubulin Alpha 8 (Tuba8) mutated in Polymicrogyria, a developmental malformation of the cortex (45), the Hes Family BHLH Transcription Factor 5 (Hes5) activated downstream of the Notch pathway and largely involved in neuronal differentiation, the mitochondrial ribosomal protein s28 (Mrps28) which mutation severely impairs the development of the nervous system (46), the Phospholipase C Delta 1 (Plcd1) important for neuronal development and function of mature neurons, and the Pyridoxamine 5’-Phosphate Oxidase (Pnpo) involved in the synthesis of vitamin B6 and which mutation causes a form of neonatal epileptic encephalopathy and motor neuron disease. Gene ontology analysis for the 43 DEGs in the WT v s Isl1^+/Cre::Meis2^LoxP/LoxP comparison and for the 107 DEGs in the Isl1^+/Cre v s Isl1^+/Cre::Meis2^LoxP/LoxP comparison revealed significant relevant hits with many terms associated with neuronal projections and functions (Figure 3B and C; Figure 3 Supplementary Table 2; Figure 3 Supplementary 2 and 3). These include subsets for the GO term associated with synapse, dendrite and axons and more specifically with GABAergic synapses, dendritic shaft or postsynaptic membrane. None of these GO terms were significantly enriched in the WT v s Isl1^+/Cre comparison (Figure 3B and Figure 3 Supplementary 2 and 3) which overall showed lower enrichment scores and p values than in the two other datasets. It is important to note that many of the genes associated with neuron projection or synapse that were present in either WT vs Isl1^+/Cre::Meis2^LoxP/LoxP dataset or Isl1^+/Cre vs Isl1^+/Cre::Meis2^LoxP/LoxP dataset, such as Oprd1, Calb2, Whrn, Lrp2, Lypd6, Grid1 and Rps21, failed to enter the list of the 10 best DEGs, either because their fold changes were below the cutoff, either because their p-values were close to but higher than 0.05. Interestingly, a significant association with the GO term Cadherin in the Isl1^+/Cre v s Isl1^+/Cre::Meis2^LoxP/LoxP comparison points at the protocadherin family in which several members were down-regulated (Figure 3C). Finally, comparing these genes to scRNAseq analysis in adult DRG neurons (17) showed that most of them are expressed by Meis2-expressing DRG sensory neuron subtypes (Figure 3 Supplementary 4). These molecular analyses strongly support a role for Meis2 in regulating embryonic target-field innervation. We thus investigated this hypothesis, and 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 3E). In WT newborn mice, numerous Nefh⁺ sensory fibers surround all dermal papillae of the hairy skin and footpad of the glabrous skin, whereas in Wnt1^Cre::Meis2^LoxP/LoxP littermates, very few Nefh⁺ sensory fibers are present and they poorly innervate the dermal papillae and footpads.

Figure 3:

Meis2 gene is necessary for SA-LTMR morphology and function only in the glabrous skin

LTMRs form specialized sensory endings in a variety of end-organs specialized to shape the mechanoreceptor properties. We used the Isl1^+/Cre::Meis2^LoxP/LoxP mice to assess the effects of late loss of Meis2 on LTMR structure and function, and to investigate 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 (12, 13, 47).

We recorded 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 (12, 13). 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 rapidly-adapting or slowly-adapting mechanoreceptors (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 (Figure 4A). Such “tap” units have been found in several mice with deficits in sensory mechano-transduction (47, 48). 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/LoxPmice, but it was clear that many fibers in this sample had mechanical thresholds similar to those in the wild type (Figure 4B). The response of the same SA- LTMRs to a 25 Hz sinusoidal stimulus was unchanged in Isl 1^Cre/+::Meis2^LoxP/LoxP mice compared to controls (Figure 4B). 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 4C). 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 1).

Figure 4:

In both the glabrous and the hairy skin, Merkel cells are innervated by Slowly-Adapting Mechanoreceptor type I (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 Isl1^+/Cre (Figure 4D). 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 (Figure 4E). 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 (Figure 4F).

Altogether, these data indicate that Meis2 is necessary for Merkel cells innervation in the glabrous, but not in the hairy skin. In addition, electrophysiological recordings indicate that amongst SA- LTMRs, there was a light loss of sensitivity that could be associated with poor innervation of Merkel cells in the glabrous skin.

Meis2 is necessary for RA-LTMRs structure and function

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 5A, Figure 5 Supplementary video 1-2) with decreased complexity of the Nefh⁺ fibers within the corpuscle as shown by quantification of the average number of time that fibers cross the midline of the terminal structure. However recordings from RA- LTMRs innervating these structures in Isl1^Cre/+::Meis2^LoxP/LoxP animals showed largely normal physiological properties (Figure 5B). Thus RA-LTMRs recorded from Isl 1^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 5B).

Figure 5:

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 (4). Similar to Meissner corpuscles, they are tuned to frequencies between 10 and 50 Hz (12). 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 Isl1^Cre/+ (Figure 6A). Our analysis revealed significant decreases in both the number of plexus branch points and in the number of innervated hair follicles (Figure 6B).

Figure 6:

Consistent with the hypo-innervation of hair follicles in Isl1^Cre/+::Meis2^LoxP/LoxP, we observed robust deficits in the mechanosensitivity of RA-LTMRs in the hairy skin (Figure 6C and Figure 4 Supplementary 1B). Thus, we needed sinusoids of significantly larger amplitudes to evoke the first (threshold) spike in RA-LTMRs. We therefore measured the total number of spikes evoked by a sinusoid stimulus (25 Hz) of gradually increasing amplitude. Again, RA-LTMRs fired considerably less in Isl1^+/Cre::Meis2^LoxP/LoxP mutant than in control mice. This finding was confirmed using a series of vibration steps of increasing amplitudes again demonstrating decreased firing in response to 25 Hz vibration stimuli (Figure 6C). Thus, the functional deficits in RA-LTMRs correlate well with the defects in LTMR cutaneous projections we observed in Isl1^+/Cre::Meis2^LoxP/LoxP mutant mice.

Finally, D-Hair mechanoreceptors or Aδ−LTMRs are the most sensitive skin mechanoreceptors with very large receptive fields (7, 13, 49). They form lanceolate endings, are thinly myelinated and are activated by movement of the smaller zigzag hairs (4). Consistent with the lack of Meis2 expression in this population reported by single-cell RNA-seq databases, Aδ fibers D-hair in the hairy skin showed similar vibration responses in WT and Isl1^+/Cre::Meis2^LoxP/LoxP mice (Figure 5 Supplementary 1).

Discussion

The function of the Meis family of TFs in post-mitotic neurons has only been marginally addressed (30, 50, 51). Here, we showed that Meis2 is selectively expressed by subpopulations of early post- mitotic cutaneous LTMR and proprioceptive neurons during development of both mouse and chick, highlighting the conserved Meis2 expression across vertebrate species in those neurons. Our results on Meis2 expression are in agreement with previous combined single cell RNAseq (scRNAseq) analysis and genetic tracing reporting Meis2 in proprioceptive neurons, Aβ field-LTMR, Aβ-SA1-LTMR, and Aβ-RA-LTMR, but not in in C-LTMR, Aδ-LTMR, peptidergic and non-peptidergic nociceptive neurons (11, 17, 18). We unambiguously demonstrate that Meis2 differentially regulates target-field innervation and function of post-mitotic LTMR neurons. Meis2 inactivation in post-mitotic sensory neurons modified their embryonic transcriptomic profile and differentially impaired adult LTMR projections and functions without affecting their survival and molecular subtype identity.

The morphological and functional phenotypes we report following specific Meis2 gene inactivation in post-mitotic sensory neurons are consistent with its expression pattern, and ultimately, both defective morphological and electrophysiological responses result 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 responses to mechanical stimuli is consistent with Meis2 being expressed by Aβ-SA1-LTMR neurons. Interestingly, Meis2 gene inactivation compromises Merkel cells innervation and electrophysiological responses in the glabrous skin but not in touch domes of the hairy skin where innervation appeared unchanged. This difference supports previous work suggesting that the primary afferents innervating Merkel cells in the glabrous and the hairy skin maybe different (52, 53). 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, and another Ntrk3/Nefh-positive that do not depend on Ntrk1 signaling during development (52). 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 (54).

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 are 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 wider 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. Indeed, the electrophysiology methods used here can only identify sensory afferents that have a mechanosensitive receptive field. Primary afferents that have an axon in the skin but no mechanosensitvity can only be identified with a so- called electrical search protocol (47, 48) which was not used here. It is therefore quite likely that many primary afferents that failed to form endings would not be recorded in these experiments e.g. SA-LTMRs and RA-LTMRs that fail to innervate end-organs (Fig.4-6). 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 test performed within low forces and challenging light touch sensation could distinguish Merkel cells from Meissner corpuscles dysfunctions. Mice depleted of Merkel cells performed normally on this test while mice mutated at the Ntrk2 locus with Meissner corpuscle innervation deficits were less sensitive in response to filament within the 0.02-0.6 g range (55). Thus, our result in the Von Frey test likely reflects aberrant functioning of the RAM LTMR-Meissner corpuscle complex. Finally, the unaltered D-hair fibers electrophysiological responses and the normal noxious responses in the hot plate setting are consistent with the absence of Meis2 expression in Aδ-LTMR, peptidergic and non-peptidergic neurons. Surprisingly, although Meis2 is expressed in proprioceptive neurons (17, 18, 32), their function appeared not to be affected as seen by normal gait behavior in catwalk analysis. This is in agreement with studies in which HoxC8 inactivation, a classical Meis TF co-factor expressed by proprioceptive neurons from E11.5 to postnatal stages, affected neither proprioceptive neurons early molecular identity nor their survival (32). From our data, we could not conclude whether SA-LTMR electrophysiological responses are differentially affected in the glabrous versus hairy skin of Meis2 mutant as suggested by histological analysis. Similarly, the decreased sensitivity of Meis2 mutant mice in the cotton swab assay and the morphological defects of Meissner corpuscles evidenced in histological analysis do not correlate with RA-LTMR electrophysiological responses for which a tendency to decreased responses were however measured. The later might result from an insufficient number of fibers recording, whereas the first may be due of pooling SA-LTMR from both the hairy and glabrous skin.

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 (56, 57). 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 (58, 59). The proneural function of Ngn1 and Ngn2 genes in neural crest cells, the precursors of DRG sensory neurons is well demonstrated (9, 60), and several terminal selector genes shaping the different DRG sensory subpopulations have also been clearly identified including for cutaneous LTMRs (8, 9, 20, 53, 61–64). 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 (65).

In humans, at least 17 different mutations in the Meis2 gene have been associated with neurodevelopmental delay (24–26), emphasizing its essential function in neuronal differentiation. Meis2 function in late differentiation of post-mitotic peripheral sensory neurons adds 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 (50, 66–73) argues for diverse functions ranging from regulation of neuroblasts cell-cycle exit, to cell-fate decision, neurogenesis, neuronal specification, neurites outgrowth, synaptogenesis and maintenance of mature neurons. 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. However, our transcriptomic analyses strongly support that in LTMRs, Meis2 regulates other types 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 targets. Surprisingly, our RNAseq analysis only revealed 10 DEGs that could be unambiguously attributed to Meis2 activity. All of these genes are expressed in adult LTMRs (17), suggesting that they exert specific functions in the maintenance of these sensory subclasses, and that their up- or down- regulation might affect LTMRs maturation, but did not show any GEO term enrichment. Separate gene ontology analyses of our datasets however revealed alterations in pathways associated to synapses function and neurons projections. Therefore, from our results, we cannot exclude that dysregulation of those genes is secondary to the changed expression of one or more of the 10 above DEGs. Two GABA(A) Receptor Subunits (GABRA1 and GABRA4), the K-Cl cotransporters SLC12A5 associated to GABAergic neurotransmission, but also the glutamate receptor subunits GRID1 and GRIK3 are down or up-regulated in Meis2 mutants, questioning whether an imbalance Gabaergic and Glutamatergic transmission is responsible for Meis2 sensory phenotypes. Interestingly, Meis2 inactivation seems to interfere with the embryonic expression of many members of the protocadherin family, and the protocadherin γ cluster (Pcdhg) in particular has recently been highlighted as essential for building central and peripheral LTMRs innervation and synapses and establish proper peripheral target-field innervation and touch sensation (23). Finally, previous work on Islet1 conditional deletion in DRG sensory neurons report considerable changes in gene expression. Early homozygous Islet1 deletion results in increased sensory neurons apoptosis and a loss of Ntrk1- and Ntrk2-positive neurons, whereas late deletion seems to only affect the nociceptive subpopulations (74). Whereas the distal projection defects we report in Wnt^Cre::Meis2^LoxP/LoxP mutant can only be attributed to Meis2 inactivation, it is possible that among the DEGs we identified in Isl1^Cre/+::Meis2^LoxP/LoxP, some are epistatically regulated by the heterozygous Islet1 deletion in addition to Meis2 homozygous deletion. Such epistasis has been previously shown for Islet1 and the transcription factor Brn3a (74, 75).

In conclusion, this study reveals a novel function for the Meis2 transcription factor in selectively regulating target field innervation of LTMR neurons. More broadly, it opens new perspectives to molecularly understand how Meis2 is linked to neuronal development. Together with studies on Meis2 function in the SVZ where it is necessary to maintain the neurogenic effect of Pax6 in neural progenitors and is later expressed in their mature progenies (51), our results raise the possibility that this TF sets up a lineage specific platform on which various specific co-factors in turn participate in additional and/or subsequent steps of the neuronal differentiation program.

Acknowledgements

Stéphanie Ventéo for help with ChTx experiments, Anne-Laure Bonnefont for cat walk experiments, all staff at animal house and in particular Flora for great help at the animal facility. Yves Dusabyinema, Benjamin Gillet and Sandrine Hughes at the IGFL sequencing platform (PSI) for Illumina sequencing.

Contribution

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

Animals

All procedures involving animals and their care were conducted according to European Parliament Directive 2010/63/EU and the 22 September 2010 Council on the protection of animals, and were approved by the French Ministry of research (APAFIS#17869-2018112914501928 v2, June the 4^th of 2020).

Mice strains

Wnt1^Cre and Islet1^+/Cre mice were previously described(76, 77). We previously reported on the conditional mutant strain for Meis2 (Meis2^LoxP/LoxP) used in the present study(78). To generate this strain, exon 8 of the Meis2 gene was flanked by the LoxP recognition elements for the Cre recombinase at the ITL (Ingenious Targeting Laboratory, NY, USA) using standard homologous recombination technology in mouse embryonic stem cells. FLP-FRT recombination was used to remove the neomycin selection cassette and the Meis2^LoxP/LoxP mutant mice were backcrossed for at least 8 generations onto the C57BL/6 background before use. Primers used to genotype the different strains were: Meis2 sense 5’-TGT TGG GAT CTG GTG ACT TG-3’; Meis2 antisense 5’-ACT TCA TGG GCT CCT CAC AG-3’; Cre sense 5’-TGC CAG GAT CAG GGT TAA AG-3’; Cre antisense 5’-GCT TGC ATG ATC TCC GGT AT-3’. Mice were kept in an animal facility and gestational stages were determined according to the date of the vaginal plug.

For retro-tracing experiments, newborn pups were anesthetized on ice and cholera toxin B coupled to Cholera Toxin Subunit B conjugated with Alexa Fluor 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 12- weeks 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) (79) and the right hind limb bud was surgically removed as previously reported (37). 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 (30).

In situ hybridization (ISH)

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 one fluorescein-labeled probe. Fluorescein-labeled probe was first revealed after overnight incubation with AP-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 was 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.

Immunochemistry

Immunochemistry was performed as previously described (30). In situ hybridized sections or new sections were washed 3×10 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 3×10 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 and chick (diluted 1:100, Developmental Studies Hybridoma Bank); goat anti-TrkB (1/2000; R and D Systems); goat anti-TrkC (1/1000; R and D Systems); rabbit anti-TrkA (1/500; Millipore); rabbit anti- parvalbumin antibody (1:500, Swant); guinea pig anti-calcitonin gene-related peptide (CGRP) antibody (1:500, Peninsula Laboratories); rabbit anti-Nefh (Sigma N4142, rabbit 1:1000), guinea pig anti-c-maf (generous gift of C. Birchmeier, MDC, 1/10000), rabbit anti-TrpV1 (Sigma V2764, 1:1000), mouse-anti-S100β (Sigma S2532, 1:1000), rabbit anti-Phospho-CREB (Cell Signaling, 87G3, 1/200, Germany); Goat anti-c Ret (R and D Systems, Cat#AF482, 1/100) and mouse or rabbit anti-Meis2 (Sigma Aldrich, WH0004212M1 or Abcam ab244267, 1/500) and rat anti-CK8 (DSHB, TROMA-I, 1/100). The chick rabbit anti-TrkB and C antibodies were a generous gift from LF Reichardt, UCSF and have been previously reported.

Whole-mount immunohistochemistry

Whole-mount immunohistochemistry of adult mice back hairy skin was performed as described elsewhere (80). Briefly, mice were euthanized by CO₂ 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 min. 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(30), air dried frozen sections were washed in water then stained with hematoxylin for 1Lmin 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).

RNA-sequencing and analysis

DRGs were dissected from E18.5 mouse embryos, collected in lysis buffer and stored at −80°C until RNA extraction with RNeasy extraction kit (Qiagen). After mRNA purification using the NEBNext® Poly(A) mRNA Magnetic Isolation (NEB), libraries were prepared with the CORALL mRNA-Seq Library Prep Kits with UDIs (Lexogen) following manufacturer’s recommendations. After a qPCR assay to determine the optimal PCR cycle number for endpoint PCR, 14 PCR cycles were completed to finalize the library preparation. Quantitation and quality assessment of each library were performed using Qubit 4.0 (HS DNA kit, Thermofisher) and 4150 Tapestation analyzer (D5000 ScreenTape kit, Agilent). Indexed libraries were sequenced in an equimolar manner on NextSeq 500 Illumina sequencer. Sequencing conditions were as follow: denatured libraries were loaded on a HighOutput flowcell kit v2 and sequenced in single-end 84pb reads. Data were extracted and processed following Illumina recommendations. After a quality check of the fastq files with FastQC, UMI sequences were extracted with UMI tools (version 1.1.2) (81) default parameters followed by STAR alignement (version 2.7.10) (82) on mm10 genome and removal of PCR duplicate with UMI tools, default parameters. Uniquely mapped sequences from the STAR output files (bam format) were then used for further analysis. HT-seq count (version 0.6) (83) was used to aggregate read count per gene followed by differential gene expression analysis with Limma voom on Galaxy (version 3.50.1) (84). Data are available under the following accession codes: GSE223788. Only genes that exhibited more than 100 reads in any of the samples were kept in the analysis. Genes with more than 1.2-fold differential expression and p-value <0.05 were used for gene ontology analysis (https://david.ncifcrf.gov) using the list of expressed genes in our experiment as background.

Mouse skin-nerve preparation and sensory afferent recordings

Cutaneous sensory fiber recordings were performed using the ex vivo skin-nerve preparation as previously described (12). Mice were euthanized by CO2 inhalation for 2–4Lmin 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 32L°C synthetic interstitial fluid: 123LmM NaCl, 3.5LmM KCl, 0.7LmM MgSO4, 1.7LmM NaH2PO4, 2.0LmM CaCl2, 9.5LmM sodium gluconate, 5.5LmM glucose, 7.5LmM sucrose and 10LmM 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 (1LHz, square pulses of 50–500Lms) were delivered to single-unit receptive fields to measure conduction velocity and enable classification as C-fibers (velocity <1.2LmLs−1), Aδ-fibers (1.2–10LmLs−1) or Aβ-fibers (>10LmLs−1). Mechanical stimulation of the receptive fields of neurons were performed using a piezo actuator (Physik Instrumente, catalog no. P-841.60) and a double-ended 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 (5LHz, 25LHz and 50LHz, distortions introduced by the in-series force sensor precluded using frequencies >50LHz) with increasing amplitude over six steps (peak-to-peak amplitudes of ∼6–65LmN; 20 cycles per step), and a dynamic stimulus sequence with four ramp-and-hold waveforms with varying probe deflection velocities (3Ls duration; 0.075, 0.15, 0.45 and 1.5LmmLs−1; average amplitude 100LmN). 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 (3Ls duration; ranging from ∼10LmN to 260LmN). 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.

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 testing day, mice were placed in the chamber one hour before Von Frey filaments application. The test was performed as previously described (55). 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. During each application, bend filament was maintained for four to five seconds. 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² 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 “wet- dog”, 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.

Figure 1 supplementary 1: Meis2 mRNA expression in LTMR 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 Meis2⁺ only 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. We estimated that 16.0 ± 1.2% (mean ± SEM; n = 3) of Meis2-positive neurons co-expressed Ntrk2 and that 43.7 ± 2.1% co-expressed Ntrk3. Conversely, Meis2 was co-expressed in 55.9 ± 3.1% of Ntrk2- positive neurons and in 79.1 ± 3.4% of Ntrk3-positive neurons. 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. We estimated that 28.7 ± 2.5 of Meis2-positive neurons co-expressed Ntrk2 and that 57.8 ± 5.1 co- expressed Ntrk3. Conversely, Meis2 was co-expressed in 52.1 ± 4.5 % of Ntrk2-positive neurons and 63.8 ± 5.1 ± 3.4% of Ntrk3-positive neurons. 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 expression in Ntrk1⁺/Meis2⁺ neurons is very low at the limit of detection. Dashed lines delineate the DRG. Scale bar=50µm. D) ISH for Meis2 (blue) combined with IF against Ntrk1, Calca, Trpv1, Pvalb (red) and Islet1 or Nefh (green) in adult mouse lumbar DRG. Arrows point at Meis2⁺ neurons. Arrowheads point at Ntrk1⁺, Calca⁺, Trpv1⁺ or Pvalb⁺ neurons. Stars indicate neurons that are both positive for Meis2 and Ntrk1, Calca, Trpv1 or Pvalb. Graphs showing the percentage of Meis2⁺ neurons showing immunoreactivity for Ntrk1, Trpv1, Calca, Pvalb and Nefh, and the percentage of Isl1⁺, Ntrk1⁺, Calca⁺, Trpv1⁺, Pvalb⁺ or Nefh⁺ neurons co-expressing Meis2. Scale bar=50µm.

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 becomes progressively restricted to the VL population of sensory neurons which represents the Ntrk2⁺ and Ntrk3⁺ populations of mechano- and proprioceptive neurons. 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, B) Combined ISH for Meis2 (blue) with IF against Islet1 (green) and Ntrk2 or Ntrk3 (red) in chick lumbar DRGs at HH29 showed that Meis2 expression is shared between the Ntrk2⁺ and Ntrk3⁺ subpopulations of sensory neurons. Arrowheads point at Meis2⁺/Islet1⁺/Ntrk2^- and Meis2⁺/Islet1⁺/Ntrk3^- neurons; arrows point at Meis2⁺/Islet1⁺/Ntrk2⁺ and Meis2⁺/Islet1⁺/Ntrk3⁺ neurons; stars indicate Meis2^-/Islet1⁺/Ntrk2⁺ and Meis2^-/Islet1⁺/Ntrk3⁺ neurons. 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 Ntrk2- and Ntrk3-positive neurons co-expressed Meis2 mRNA respectively. C) Combined ISH for Meis2 (blue) with IF against Islet1 (green) and Ntrk1 (red) in chick lumbar DRGs at HH29 showed that Meis2 is excluded from the Ntrk1⁺ population of sensory neurons. Arrowheads point at Meis2⁺/Islet1⁺/Ntrk1^- neurons. Enlargement is indicated by a dashed square. D) Representative images showing immunostaining for Ntrk2 or Ntrk3 (red) and Islet1 (green) on HH36 chick embryos DRG sections following limb ablation at HH18. Dashed lines delineate the DRGs. Arrows point at Ntrk2-positive VL neurons. Arrow heads point at Ntrk2-positive DL neurons. C) Box plots showing the percentage of Ntrk2 and Ntrk3 VL neurons and of Ntrk2 DL neurons. *** p≤0.0005; ns= not significant following Student t-test. n=3 chick embryos.

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) Upper right, representative images showing that new-born Wnt1^Cre::Meis2^LoxP/LoxP pups mutant pups were smaller in size at birth compared to WT littermates and were unable to stand on their legs. Upper left, representative images showing that new-born Wnt1^Cre::Meis2^LoxP/LoxP pups mutant pups exhibited a cleft palate phenotype (white arrowheads). Bottom, representatives images showing coronal sections of new-born WT and Wnt1^Cre::Meis2^LoxP/LoxP heads stained by eosin/hematoxylin treatment. Black arrows indicate the cleft palate.

Figure 1 Supplementary Table: 1Isl1^+/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-t- test analysis showed no significant differences for any of the parameters.

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 Supplementary Table 1: Table showing the results of the bulk RNAseq analysis.

Figure 3 Supplementary Table 2: Table showing the results of the GO-terms analysis performed with David.

Figure 3 Supplementary 1: A) Volcano plots showing differentially expressed genes (DEGs) in red including genes with of minimal fold change of 20% (n=3; p<0.05). Volcano plots show the comparison between the 3 genotypes (WT, Isl1^Cre/+and Isl1^Cre/+::Meis2^LoxP/LoxP). Plots reporting Meis2 and Isl1 mRNA expression are in blue and green respectively. B) Graphs showing the individual number of reads for Meis2 and Isl1 genes in each genotype. n=3; p values are indicated. C) Heat map showing the 10 best DEGs (significantly and differentially expressed in both WT vs Isl1^Cre/+::Meis2^LoxP/LoxP and Isl1^Cre/+ vs Isl1^Cre/+::Meis2^LoxP/LoxP). D) Expression of the 10 best DEGs in the different adult DRG sensory neurons subtypes according to Usoskin et al. (17). Red frames indicate Meis2-expressing subpopulations of LTMR and proprioceptive neurons. Note that all DEG are expressed in Meis2-expressing neurons.

Figure 3 Supplementary 2: Best David GEO terms for each dataset in the RNAseq paired-analysis. Graphs show the fold enrichment and –log10(p value) of the first 12 best GEO terms sorted according to their p value and fold enrichment for A) WT vs Isl1^+/Cre::Meis2^LoxP/LoxP, B) Isl1^+/Cre vs Isl1^+/Cre::Meis2^LoxP/LoxP and C) WT vs Isl1^+/Cre. Only terms associated with more than 5 genes were considered. Terms associated with more than 15 genes were excluded as they referred to very generic keywords. Red bars indicate terms associated to the nervous system. Note that none of the best GO terms in the WT vs Isl1^+/Cre dataset associates to the nervous system, and that most of the fold enrichment and p values are lower in the WT vs Isl^1+/Cre dataset compared to the other datasets.

Figure 3 Supplementary 3: Best David terms other than GEO for each dataset in the RNAseq paired- analysis. Graphs show the fold enrichment and –log10(p value) of all the terms other than GEO sorted according to their p value and fold enrichment for A) WT vs Isl1^+/Cre::Meis2^LoxP/LoxP, B) Isl1^+/Cre vs Isl1^+/Cre::Meis2^LoxP/LoxP and C) WT vs Isl1^+/Cre. Only terms associated with more than 5 genes were considered. Terms associated with more than 15 genes were excluded as they referred to very generic keywords. Red bars indicate terms associated to the nervous system. Note that none of these terms in the WT vs Isl1^+/Cre dataset associates to the nervous system, and that most of the fold enrichment and p values are lower in the WT vs Isl1^+/Cre dataset compared to the 2 others datasets.

Figure 3 Supplementary 4: Comparison of DEGs in Meis2 mutant embryonic DRG neurons with their expression in the different adult DRG sensory neurons subtypes. A) Heat map showing the expression of DEGs of the WT vs Isl1^+/Cre::Meis2^LoxP/LoxP - (A) and of the Isl1^+/Cre vs Isl1^+/Cre::Meis2^LoxP/LoxP (B) datasets in the different adult DRG sensory neurons subtypes according to Usoskin et al. (17). Red frames indicate Meis2-expressing subpopulations.

Figure 4 Supplementary 1: Graphs showing the number of recorded fibers in the nerve-skin preparation in WT and Isl1^+/Cre::Meis2^LoxP/LoxP.

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

Figure 5 Supplementary video 2: Meissner corpuscles in Isl1^+/Cre::Meis2^LoxP/LoxP. 3D visualization of Meissner corpuscles in adult Isl1^+/Cre::Meis2^LoxP/LoxP glabrous skin visualized by IF against Nefh and S100β.

Figure 5 Supplementary 1: Normal electrophysiological responses of D-Hair mechanoreceptors following Meis2 gene inactivation. A) D-Hair mechanoreceptors 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 mechanoreceptors 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. Trace shows the type of stimulation and red square indicates the time frame during when the number of spikes was calculated.