# Touch receptor end-organ innervation and function require sensory neuron expression of the transcription factor Meis2

Simon Desiderio, Frederick Schwaller, Kevin Tartour, Kiran Padmanabhan, Gary Lewin, Patrick Carroll, Frederic Marmigere

► To cite this version:

Simon Desiderio, Frederick Schwaller, Kevin Tartour, Kiran Padmanabhan, Gary Lewin, et al.. Touch receptor end-organ innervation and function require sensory neuron expression of the transcription factor Meis2. eLife, 2024, 12, pp.RP89287. 10.7554/eLife.89287. hal-04480147

HAL Id: hal-04480147

https://hal.science/hal-04480147v1

Submitted on 6 Nov 2024

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

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

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

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

# Touch receptor end-organ innervation and function require sensory neuron expression of the transcription factor Meis2

Simon Desiderio¹†, Frederick Schwaller²†, Kevin Tartour³, Kiran Padmanabhan³, Gary R Lewin², Patrick Carroll¹, Frederic Marmigere³*

¹Institute for Neurosciences of Montpellier, University of Montpellier, INSERM U 1298, Montpellier, France; ²Department of Neuroscience, Max-Delbrück Centre for Molecular Medicine, Berlin-Buch, Germany; ³IGFL, UMR 5242 CNRS/ENS, Lyon, France

**Abstract** Touch sensation is primarily encoded by mechanoreceptors, called low-threshold mechanoreceptors (LTMRs), with their cell bodies in the dorsal root ganglia. Because of their great diversity in terms of molecular signature, terminal endings morphology, and electrophysiological properties, mirroring the complexity of tactile experience, LTMRs are a model of choice to study the molecular cues differentially controlling neuronal diversification. While the transcriptional codes that define different LTMR subtypes have been extensively studied, the molecular players that participate in their late maturation and in particular in the striking diversity of their end-organ morphological specialization are largely unknown. Here we identified the TALE homeodomain transcription factor Meis2 as a key regulator of LTMRs target-field innervation in mice. *Meis2* is specifically expressed in cutaneous LTMRs, and its expression depends on target-derived signals. While LTMRs lacking *Meis2* survived and are normally specified, their end-organ innervations, electrophysiological properties, and transcriptome are differentially and markedly affected, resulting in impaired sensory-evoked behavioral responses. These data establish *Meis2* as a major transcriptional regulator controlling the orderly formation of sensory neurons innervating peripheral end organs required for light touch.

## eLife assessment

This **fundamental** study identifies the homeodomain transcription factor Meis2 as a transcriptional regulator of maturation and end-organ innervation of low-threshold mechanoreceptors (LTMRs) in the dorsal root ganglia (DRG) of mice. The authors use histology, behavioral tests, RNA-sequencing, and electrophysiological recordings to provide evidence that conditional deletion of Meis2 in post-mitotic DRG neurons causes gene expression changes together with targeting errors and altered sensory neuron responses, ultimately resulting in reduced sensitivity to light touch in mutant animals. The data presented are **convincing**, the discussion comprehensive, and the conclusions drawn justified.

## Introduction

Tactile stimuli like brush, light pressure, or roughness engage highly specialized and diverse arrays of mechanoreceptors in both the hairy and glabrous skin (Abraira and Ginty, 2013; Delmas et al., 2011; Handler and Ginty, 2021; Lechner and Lewin, 2013; Wu et al., 2021; Zimmerman et al., 2014;

<footnote>*For correspondence: frederic.marmigere@ens-lyon.fr</footnote>

<footnote>†These authors contributed equally to this work</footnote>

<footnote>**Competing interest:** The authors declare that no competing interests exist.</footnote>

<footnote>**Funding:** See page 21</footnote>

<footnote>**Preprint posted**
20 May 2023</footnote>

<footnote>**Sent for Review**
06 July 2023</footnote>

<footnote>**Reviewed preprint posted**
02 October 2023</footnote>

<footnote>**Reviewed preprint revised**
05 February 2024</footnote>

<footnote>**Version of Record published**
22 February 2024</footnote>

<footnote>**Reviewing Editor:** K VijayRaghavan, National Centre for Biological Sciences, Tata Institute of Fundamental Research, India</footnote>

<footnote>© Copyright Desiderio, Schwaller et al. This article is distributed under the terms of the **Creative Commons Attribution License**, which permits unrestricted use and redistribution provided that the original author and source are credited.</footnote>
Li et al., 2011). 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 and have distinctive physiological properties and projections within the spinal cord and periphery (Lallemend and Ernfors, 2012; Marmigère and Ernfors, 2007; Vermeiren et al., 2020).

Cutaneous mechanoreceptors or low threshold mechanoreceptors (LMTRs) exhibit a variety of specialized terminal endings in the hairy and glabrous skin with strikingly unique morphologies (Abraira and Ginty, 2013; Delmas et al., 2011; Handler and Ginty, 2021; Lechner and Lewin, 2013; Zimmerman et al., 2014; Li et al., 2011; Sharma et al., 2020; Schwaller et al., 2021). 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 (Lechner and Lewin, 2013). 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 (Li et al., 2011; Walcher et al., 2018). A small number of LTMRs in the hairy skin are not activated by hair movement but show properties of rapidly adapting mechanoreceptors (Lewin and McMahon, 1991). These were originally characterized as so-called field receptors (Lewin and McMahon, 1991; Burgess and Horch, 1973) and were recently shown to form circumferential endings around hair follicles (Bai et al., 2015). LTMRs tuned to high-frequency vibration are called Aβ-RAII and innervate Pacinian corpuscles deep in the skin or on the bone (Schwaller et al., 2021). Ruffini endings that are thought to be innervated by stretch-sensitive mechanoreceptors (Aβ-SAll) remain poorly characterized in mice (Handler and Ginty, 2021). 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 (Wu et al., 2021; Sharma et al., 2020; Bai et al., 2015; Usoskin et al., 2015; Zheng et al., 2019).

Cracking the transcriptional codes supporting sensory neuron identity and diversification has been the object of tremendous efforts in the last decades (Wu et al., 2021; Lallemend and Ernfors, 2012; Marmigère and Ernfors, 2007; Sharma et al., 2020; Usoskin et al., 2015; Zheng et al., 2019). For instance, the functions of specification factors or terminal selectors, like Maf, Shox2, Runx3, Pea3, and ER81, have been functionally implicated in LTMR segregation (Lallemend and Ernfors, 2012; Marmigère and Ernfors, 2007; Scott et al., 2011; Abdo et al., 2011; Arber et al., 2000; Hu et al., 2012). Whereas the specific function of adhesion molecules in shaping the assembly of touch circuitry is being unraveled (Meltzer et al., 2023), 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 (Sharma et al., 2020; Usoskin et al., 2015; Zheng et al., 2019). Its mutation in humans causes severe neurodevelopmental defects (Gangfuß et al., 2021; Giliberti et al., 2020; Shimojima et al., 2017), and somatic mutations of its DNA consensus binding site are associated with neurodevelopmental defects (Bae et al., 2022). It belongs to a highly conserved homeodomain family containing three members in mammals, Meis1, Meis2, and Meis3 (Geerts et al., 2003; Longobardi et al., 2014), and Meis1 is necessary for target-field innervation of sympathetic peripheral neurons (Bouilloux et al., 2016). We thus wondered if Meis2 could also be a pertinent regulator of late primary sensory neuron differentiation.

Here, we show that Meis2 regulates 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 functions in specification and/or target-field innervation. We generated mice carrying an inactive Meis2 gene in postmitotic sensory neurons. While these animals are healthy and viable and do not exhibit any neuronal loss, they display tactile sensory defects in
electrophysiological and behavioral assays. Consistent with these findings, we identified major morpho-
logical 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 (Figure 1A, Figure 1—figure supplements 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 ventrolateral part of the DRG where LTMRs and proprioceptors are located (Rifkin et al., 2000; Figure 1—figure supplement 2A). In both species, *Meis2*-positive cells also expressed the pan-neuronal marker Islet1, indicating that they are postmitotic 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 *Ret* mRNAs in E14.5 and E18.5 mouse embryonic DRG (Figure 1—figure supplement 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 *Ret* in E14.5 mouse DRG showed that virtually all large *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—figure supplement 2B).

In mouse, comparison of *Meis2* mRNA expression to *Ntrk1*, a well-established marker for early nociceptive and thermo-sensitive neurons, showed that only a few *Meis2*-positive neurons co-expressed *Ntrk1* at E11.5 and E18.5 (Figure 1—figure supplement 1C). In chick HH29 embryos, *Meis2* expression was fully excluded from the *Ntrk1* subpopulation (Figure 1—figure supplement 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 thermosensitive populations of DRG neurons. Instead, a large proportion of *Meis2*-positive neurons co-expressed *Nehf*, a marker for large myelinated neurons including LTMR and proprioceptors, and *Pvalb*, a specific marker for proprioceptors (Figure 1—figure supplement 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 3 d later at P8 showed that many retrogradely labeled sensory neurons were also immunopositive for *Meis2*, *Maf*, *Ntrk2*, and *Ntrk3* (Figure 1B).

Altogether, our results on Meis2 co-localization with Nehf, Ntrk2, Ntrk3, Ret, Pvalb, and Maf at different embryonic and postnatal stages are consistent with previous reports on restricted Meis2 expression to the Aβ-field, Aβ-SA1 and Aβ-RA subclasses of LTMR neurons and proprioceptive neurons (Sharma et al., 2020; Usoskin et al., 2015; Zheng et al., 2019; Shin et al., 2020). The relatively lower coincidence 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 (Sharma et al., 2020; Arber et al., 2000; Shin et al.,
2020; Patel et al., 2003; de Nooij et al., 2013; Poliak et al., 2016; Wang et al., 2019). 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 (Oakley et al., 1995; Oakley et al., 1997; Calderó et al., 1998;
Figure 1C and D, Figure 1—figure supplement 2D). In HH36 embryos, about 65% of Meis2-positive
A

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

C

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

B

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

D

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

E

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

F

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

G

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

H

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

**Figure 1.** *Meis2* is expressed in subclasses of dorsal root ganglia (DRG) cutaneous mechanoreceptive neurons in mouse embryos. (A) In situ hybridization (ISH) for *Meis2* mRNA showed expression in a subpopulation of DRG sensory neurons at embryonic stages E11.5, E14.5, and E18.5, at P0, and at adult stages. Dashed lines delineate the DRG. Scale bar = 50 µm. (B) IF for *Meis2* (red) and c-Maf, Ntrk2, or Ntrk3 (blue) at P7 following injection of cholera toxin B subunit (CTB in green) in the skin of newborn pups. Note that *Meis2*/CTB⁺ retro-traced sensory neurons co-expressed c-Maf, Ntrk2, or Ntrk3 (arrows). Scale bar = 50 µm. We estimated that 30.5 ± 3.5% (mean ± SEM; n = 3) of *Meis2*-positive neurons co-expressed Ntrk2, and that 39.5 ± 5.4% co-expressed Ntrk3. Conversely, *Meis2* was co-expressed in 53.6 ± 9.4% of Ntrk2-positive neurons, and in 78.5 ± 5.0% of Ntrk3-positive neurons. *Meis2* expression depends on target-derived signals (C, D). (C) Representative images of *Meis2* mRNA expression (blue or pseudo-colored in red) and islet1 (green) in DRGs of Hamburger-Hamilton stage (HH) 36 chick embryos on the ablated and contralateral sides. Box plots showing the number of Islet1⁺/Meis2⁺ DRG neurons per section at stage HH36 following limb bud ablation. For Islet1-positive neurons, the contralateral side was considered as 100% of neurons per section. For *Meis2*-positive neurons, values represent the percentage of *Meis2*⁺ over Islet1⁺ neurons. (D) Representative images of *Meis2* mRNA expression (blue or pseudo-colored in red) and islet1 (green) in DRGs of HH29 chick embryos on the ablated and contralateral sides. Box plots showing the quantification of Islet1⁺/Meis2⁺ neurons number per section at stage HH29 on the contralateral and ablated sides. Arrowheads point to

*Figure 1 continued on next page*
**Figure 1 continued**

at remaining Meis2-positive VL neurons. Dashed lines encircle the DRGs. **p**≤0.005; ***p**≤0.0005; ns = not significant following Student's t-test. n = 3 chick embryos. Scale bar = 100 µm. Altered touch perception in *Meis2* mutant mice (E–H). (E) Box plots showing the responses following application of Von Frey filaments of different forces. *IsI1<sup>+/Cre</sup>::Meis2<sup>LoxP/LoxP</sup> mice exhibited a significantly reduced sensitivity to the 0.16, 0.4, and 0.6 g Von Frey filaments but not to higher forces filaments compared to WT and *IsI1<sup>+/Cre</sup>* littermates. *p*≤0.05; **p**≤0.005; ***p**≤0.001 following Kruskal–Wallis statistical analysis. (F) Box plots showing the dynamic touch responses when the hind paw palms of individual mice were stroked with a tapered cotton swab. Analysis showed that *IsI1<sup>+/Cre</sup>::Meis2<sup>LoxP/LoxP</sup> mice were less responsive to the stimulus than WT and *IsI1<sup>+/Cre</sup>* littermates. ***p**≤0.0001 following a one-way ANOVA statistical analysis. (G) Box plots indicating that the latency to the first signs of aversive behavior in the hot plate test is similar in all groups of mice. WT, n = 19; *IsI1<sup>+/Cre</sup>*, n = 16; *IsI1<sup>+/Cre</sup>::Meis2<sup>LoxP/LoxP</sup>, n = 9. (H) Box plots showing the number of bouts when a sticky paper tape was applied on the back skin of mice. Analysis indicated a significant decrease in the number of bouts in *IsI1<sup>+/Cre</sup>::Meis2<sup>LoxP/LoxP</sup> mice compared to WT and *IsI1<sup>+/Cre</sup>* littermates. *p*≤0.05 following a one-way ANOVA statistical analysis.

The online version of this article includes the following source data and figure supplement(s) for figure 1:

**Source data 1.** *IsI1<sup>+/Cre</sup>::Meis2<sup>LoxP/LoxP</sup>* adult mice exhibit normal locomotion.

**Figure supplement 1.** *Meis2* mRNA expression in low-threshold mechanoreceptor (LTMR) neurons of mouse dorsal root ganglia (DRG).

**Figure supplement 2.** *Meis2* is expressed in a subset of chick ventrolateral dorsal root ganglia (DRG) sensory neurons during embryogenesis.

**Figure supplement 3.** Mice with a conditional deletion of *Meis2* gene in neural crest derivatives (*Wnt1<sup>Cre</sup>*) exhibited cleft palate and died at birth.

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—figure supplement 2D). The number of *Ntrk2*-positive DL neurons was 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 postmitotic 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<sup>LoxP/LoxP</sup>*) in which the first coding exon for the homeodomain was flanked by *LoxP* sites (Figure 1—figure supplement 3A). To validate the use of our strain, we first crossed the *Meis2<sup>LoxP/LoxP</sup>* mice with the *Wnt1<sup>Cre</sup>* strain. This crossing efficiently inactivated *Meis2* in the neural crest, and *Wnt1<sup>Cre</sup>::Meis2<sup>LoxP/LoxP</sup>* newborn pups exhibited a cleft palate as previously reported in another conditional *Meis2* mouse strain (Machon et al., 2015; Figure 1—figure supplement 3B). They were, however, not viable, precluding functional and anatomical analyses at adult stages. To bypass this neural crest phenotype and more specifically address *Meis2* function in postmitotic neurons, we crossed the *Meis2<sup>LoxP/LoxP</sup>* mice with the *IsI1<sup>Cre/+</sup>* strain and focused our analysis on the *IsI1<sup>Cre/+</sup>::Meis2<sup>LoxP/LoxP</sup>* strain. Mutant pups were viable, appeared healthy, and displayed a normal palate, allowing sensory behavior investigations.

We monitored tactile-evoked behaviors in adult WT, *IsI1<sup>Cre/+</sup>* and *IsI1<sup>Cre/+</sup>::Meis2<sup>LoxP/LoxP</sup> 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.4 g to the hind paw and found the frequency of withdrawal responses to be significantly decreased in *IsI1<sup>+/Cre</sup>::Meis2<sup>LoxP/LoxP</sup>* mice compared to control WT and *IsI1<sup>+/Cre</sup>* mice between 0.16 and 0.6 g (Figure 1E), indicating that mutant mice are less responsive to light touch. No differences were observed between WT and *IsI1<sup>+/Cre</sup>* 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 (Bourane et al., 2015). Here, responses were significantly decreased in *IsI1<sup>+/Cre</sup>::Meis2<sup>LoxP/LoxP</sup>* mice compared to control WT and *IsI1<sup>+/Cre</sup>* 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, *IsI1<sup>Cre/+</sup>* and *IsI1<sup>Cre/+</sup>::Meis2<sup>LoxP/LoxP</sup>* 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 *IsI1<sup>Cre/+</sup>::Meis2<sup>LoxP/LoxP</sup>* mice compared to WT and *IsI1<sup>Cre/+</sup>* control mice
(Figure 1H). Finally, although *Meis2* and *Isl1* are both expressed by spinal motor neurons and proprioceptors (*Catela et al., 2016; Dasen et al., 2005; Ericson et al., 1992), we did not observe obvious motor deficits in *Isl1*<sup>Cre/+</sup>::*Meis2*<sup>LoxP/LoxP</sup> mice. Thus, in a catwalk analysis we found no differences in any of the gait parameters measured between WT and mutant mice (Figure 1—source data 1).

Overall, these behavioral 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*<sup>Cre/+</sup>::*Meis2*<sup>LoxP/LoxP</sup> and *Wnt1*<sup>Cre</sup>::*Meis2*<sup>LoxP/LoxP</sup> strains (Figure 2). There was no difference in the size of the DRGs between E18.5 WT and *Isl1*<sup>Cre/+</sup>::*Meis2*<sup>LoxP/LoxP</sup> 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 was 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*<sup>Cre</sup>::*Meis2*<sup>LoxP/LoxP</sup> 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—figure supplement 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.

## *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*<sup>Cre/+</sup>, and *Isl1*<sup>Cre/+</sup>::*Meis2*<sup>LoxP/LoxP</sup> E18.5 embryos.

For all analyses, consistent with the changes measured for *Meis2* and *Isl1* genes (Figure 3—figure supplement 1A and B), only DEGs with a minimal fold change of 20% and a p-value <0.05 were considered. Analyses of the dataset (n = 3; p<0.05; Figure 3, Figure 3—figure supplement 1, Figure 3—source data 1) identified 43 differentially expressed genes (DEGs) in the WT vs *Isl1*<sup>+/Cre</sup>::*Meis2*<sup>LoxP/LoxP</sup> comparison, 107 DEGs in the *Isl1*<sup>+/Cre</sup> vs *Isl1*<sup>+/Cre</sup>::*Meis2*<sup>LoxP/LoxP</sup> comparison, and 109 DEGs in the WT vs *Isl1*<sup>+/Cre</sup> comparison. Among them, only 10 DEGs were found in both WT vs *Isl1*<sup>+/Cre</sup>::*Meis2*<sup>LoxP/LoxP</sup> and *Isl1*<sup>+/Cre</sup> vs *Isl1*<sup>+/Cre</sup>::*Meis2*<sup>LoxP/LoxP</sup> comparisons (Figure 3A). Half of them were down- or upregulated (Figure 3—figure supplement 1C), and eight were found to be expressed in sensory neurons expressing *Meis2* (Figure 3—figure supplement 1D). These include three 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 (Abdollahi et al., 2009), 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) whose mutation severely impairs the development of the nervous system (Pulman et al., 2019), 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 whose mutation causes a form of neonatal epileptic encephalopathy and motor neuron disease. Gene Ontology (GO) analysis for the 43 DEGs in the WT vs *Isl1*<sup>+/Cre</sup>::*Meis2*<sup>LoxP/LoxP</sup> comparison and for the 107 DEGs in the *Isl1*<sup>+/Cre</sup> vs *Isl1*<sup>+/Cre</sup>::*Meis2*<sup>LoxP/LoxP</sup> comparison revealed significant relevant hits with many terms associated with neuronal projections and functions (Figure 3B and C, Figure 3—source data 2, Figure 3—figure supplements 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 vs *Isl1*<sup>+/Cre</sup> comparison (Figure 3B, Figure 3—figure supplements 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
<figure><img src="image_11.png" /><figcaption>**Figure 2.** *Meis2* is dispensable for low-threshold mechanoreceptor (LTMR) neuron survival and specification. (A) Box plots showing that the dorsal root ganglia (DRG) volumes along the rostrocaudal axis are similar in embryonic day (E) 16.5 WT and *Is1<sup>Cre/+</sup>::Meis2<sup>LoxP/LoxP</sup>* embryos. IF for Ntrk2 or Ntrk3 (red) and Islet1 (blue) and box plots analysis indicating that the percentage of Ntrk2+ and Ntrk3+ neurons is not affected in E16.5 *Is1<sup>Cre/+</sup>::Meis2<sup>LoxP/LoxP</sup>* Figure 2 continued on next page</figcaption></figure>
Figure 2 continued

LoxP. Dashed lines encircle the DRGs. n = 4; n.s. = not significant. Scale bar = 20 µm. (B) Representative images showing IF for Ntrk1, Ret, Ntrk2, Ntrk3, and Maf in E18.5 WT and *Isi1<sup>Cre/+</sup>::Meis2<sup>LoxP/LoxP</sup> DRGs. Box plots showing that the number of Ret+, Ntrk2+, Ntrk3+, and Maf+ LTMR neurons and of Ntrk1+ nociceptive neurons are similar in E18.5 WT and *Isi1<sup>Cre/+</sup>::Meis2<sup>LoxP/LoxP</sup> DRGs. n = 3; n.s. = not significant. Scale bar = 100 µm. (C) Representative images showing IF for Ntrk2 or Ntrk3 (green) with Pvalb or Maf (red) in P0 WT and Wnt1<sup>Cre/+</sup>::Meis2<sup>LoxP/LoxP</sup> DRGs. Box plots showing that the number of Ntrk2+ and Ntrk3+ neurons is unchanged in P0 WT and Wnt1<sup>Cre/+</sup>::Meis2<sup>LoxP/LoxP</sup> DRGs. n = 3, n.s. = not significant. Scale bar = 20 µm.

The online version of this article includes the following figure supplement(s) for figure 2:

**Figure supplement 1.** *Meis2* gene inactivation does not affect phospho-Creb expression.

synapse that were present in either WT vs *Isi1<sup>+/Cre</sup>::Meis2<sup>LoxP/LoxP</sup> dataset or *Isi1<sup>+/Cre</sup> vs *Isi1<sup>+/Cre</sup>::Meis2<sup>LoxP/LoxP</sup>* 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 or their p-values were close to but higher than 0.05. Interestingly, a significant association with the GO term Cadherin in the *Isi1<sup>+/Cre</sup>* vs *Isi1<sup>+/Cre</sup>::Meis2<sup>LoxP/LoxP</sup>* comparison points at the protocadherin family in which several members were downregulated (**Figure 3C**). Finally, comparing these genes to single-cell RNAseq (scRNAseq) analysis in adult DRG neurons (*Usoskin et al., 2015*) showed that most of them are expressed by Meis2-expressing DRG sensory neuron subtypes (**Figure 3—figure supplement 4**). These molecular analyses strongly support the role of *Meis2* in regulating embryonic target-field innervation. We thus investigated this hypothesis, and in P0 Wnt1<sup>Cre/+</sup>::Meis2<sup>LoxP/LoxP</sup>, 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<sup>Cre/+</sup>::Meis2<sup>LoxP/LoxP</sup> littermates, very few Nefh+ sensory fibers are present and they poorly innervate the dermal papillae and footpads.

**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 mech-
anoreceptor properties. We used the *Isi1<sup>+/Cre</sup>::Meis2<sup>LoxP/LoxP</sup>* mice to assess the effects of late loss of
*Meis2* on LTMR structure and function and investigate whether postmitotic *Meis2* inactivation impacts
terminal morphologies and physiological properties of LTMRs. We made recordings from single mech-
anoreceptors and probed their responses to defined mechanical stimuli in adult WT, *Isi1<sup>+/Cre</sup>*; and *Isi1<sup>+/Cre</sup>\*::Meis2<sup>LoxP/LoxP</sup>* mice using ex vivo skin nerve preparations as previously described (*Schwaller et al.,* 2021; *Walcher et al., 2018; Wetzel et al., 2007*).

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 (Schwaller et al., 2021; Walcher et al., 2018). 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 *Isi1<sup>+/Cre</sup>::Meis2<sup>LoxP/LoxP</sup>* 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 (Wetzel et al., 2007; Ranade et al., 2014). 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 *Isi1<sup>Cre/+</sup>::Meis2<sup>LoxP/LoxP</sup>* mice, 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 *Isi1<sup>Cre/+</sup>::Meis2<sup>LoxP/LoxP</sup>* mice compared to controls (Figure 4B). The response of these
<figure><img src="image_12.png" /><figcaption>**Figure 3.** *Meis2* inactivation dysregulates genes linked to neuronal projections and synaptogenesis. (A) Venn diagram comparing the number of differentially expressed genes (DEGs) between each genotype (n = 3; *p*<0.05). This comparison identified 10 DEGs that were differentially expressed compared to both control genotypes (WT or *Isl1*/**Cre** embryos), and a total of 140 genes that were differentially expressed in *Meis2* mutant compared to either WT or *Isl1*/**Cre** embryos. (B) Gene Ontology (GO) analysis for the three paired-analysis (WT vs *Isl1*/**Cre**; *Meis2**^LOXP/LOXP^; *Isl1*/**Cre** vs *Isl1*/**Cre**; *Meis2**^LOXP/LOXP^, and WT vs *Isl1*/**Cre**) datasets using DAVID and the full RNAseq gene list as background. Graphs show the comparison of the fold enrichment and the -log10(p value) of selected significant (*p*<0.05) GO or KEGG_PATHWAY terms associated to synapse and neuron projections whatever the number of genes. Blue dotted line indicates a *p*-value of 0.05. Note that following DAVID analysis GO terms associated to synapse and neuron projections were overrepresented in the WT vs *Isl1*/**Cre**; *Meis2**^LOXP/LOXP^ and the *Isl1*/**Cre** vs *Isl1*/**Cre**; *Meis2**^LOXP/LOXP^ datasets compared to the WT vs *Isl1*/**Cre** dataset. (C)</figcaption></figure>

*Figure 3 continued on next page*
Figure 3 continued

Heat maps showing the DEGs related to the GO terms synapse, neuron projection including dendrite, and protocadherin. (D) Representative images showing a strong overall deficit of Nefh⁺ (red) sensory projections innervating the dermal papillae in the hairy and the foot pads in the glabrous skin of P0 Wnt1<sup>Cre/::Meis2<sup>LoxP/LoxP</sup></sup> neonates forepaw compared to WT littermates. Dashed lines delineate the hair follicle and the epidermis. Scale bar = 50 µm.

The online version of this article includes the following source data and figure supplement(s) for figure 3:

**Source data 1.** Table showing the results of the bulk RNAseq analysis.

**Source data 2.** Table showing the results of the GO terms analysis performed with AVID.

**Figure supplement 1.** Best dysregulated genes in Meis2 mutant adult DRG.

**Figure supplement 2.** Best DAVID GO terms for each dataset in the RNAseq paired analysis.

**Figure supplement 3.** Best DAVID terms other than GO for each dataset in the RNAseq paired analysis.

**Figure supplement 4.** Comparison of differentially expressed genes (DEGs) in Meis2 mutant embryonic dorsal root ganglia (DRG) neurons with their expression in the different adult DRG sensory neuron subtypes.

fibers to ramp stimuli of increasing velocities or to increasing amplitudes of ramp and hold stimuli was also not significantly different in mutant mice compared to controls (Figure 4C).

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<sup>+/+</sup>::Meis2<sup>LoxP/LoxP</sup>* adult mice, the number of Merkel cells contacted by Nefh-positive fibers was strongly decreased compared to *Isl1<sup>+/Cre</sup>* (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<sup>Cre/+</sup>::Meis2<sup>LoxP/LoxP</sup>* 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 Merkel cells per touch dome was unchanged in mutant animals compared to WT (Figure 4F).

Altogether, these data indicate that *Meis2* is necessary for Merkel cell innervation in the glabrous, but not in the hairy skin. In addition, electrophysiological recordings indicate that among 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-LTMR 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—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<sup>Cre/+</sup>::Meis2<sup>LoxP/LoxP</sup>* animals showed largely normal physiological properties (Figure 5B). Thus RA-LTMRs recorded from *Isl1<sup>Cre/+</sup>::Meis2<sup>LoxP/LoxP</sup>* 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<sup>Cre/+</sup>::Meis2<sup>LoxP/LoxP</sup>* mutant mice to fire fewer action potentials to sinusoids and the ramp phase of a series 2 s duration ramp and hold stimuli, but these differences were not statistically significant (Figure 5B).

In the hairy skin, RA-LTMRs form longitudinal lanceolate endings parallel to the hair shaft of guard and awl/auuchene hairs and respond to hair deflection only during hair movement, but not during maintained displacement (Lechner and Lewin, 2013). Similar to Meissner corpuscles, they are tuned to frequencies between 10 and 50 Hz (Schwaller et al., 2021). Whole-mount analysis of Nefh-positive fibers in the adult back skin showed an overall decrease in the innervation density of hairs in *Isl1<sup>Cre/+</sup>::Meis2<sup>LoxP/LoxP</sup>* animals compared to *Isl1<sup>Cre/+</sup>* (Figure 6A). Our analysis revealed significant decreases in both the number of plexus branch points and the number of innervated hair follicles (Figure 6B).

Consistent with the hypo-innervation of hair follicles in *Isl1<sup>Cre/+</sup>::Meis2<sup>LoxP/LoxP</sup>* we observed robust deficits in the mechanosensitivity of RA-LTMRs in the hairy skin (Figure 6C). Thus, we needed sinusoids of significantly larger amplitudes to evoke the first (threshold) spike in RA-LTMRs. We therefore
<figure><img src="image_13.png" /><figcaption>A</figcaption></figure>

Recording

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

Hairy skin

Glabrous skin

<figure><img src="image_15.png" /><figcaption>B</figcaption></figure>

SA1-LTMR

<figure><img src="image_16.png" /><figcaption>C</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>D</figcaption></figure>

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

<figure><img src="image_19.png" /><figcaption>E</figcaption></figure>

<figure><img src="image_20.png" /><figcaption>CK8/Nefh</figcaption></figure>

<figure><img src="image_21.png" /><figcaption>F</figcaption></figure>

<figure><img src="image_22.png" /><figcaption>**Figure 4.** *Meis2* gene inactivation compromised Merkel cell innervation in the glabrous skin and increased slowly adapting mechanoreceptor (SAM) vibration threshold. (A) Graph showing the percentage of tap units among all recorded Aβ fibers in the nerve skin preparation both in the hairy and glabrous skins. The number of tap units over the number of recorded fibers is indicated. Note that Tap units are only present in both the hairy and glabrous skin of adult *Isl1+/Cre::Meis2LoxP/LoxP* mice but not in WT littermates. (B) In the hairy and glabrous skins, SAMs in *Isl1+/Cre::Meis2LoxP/LoxP* mice (n*Figure 4 continued on next page*)</figcaption></figure>
Figure 4 continued

= 22 from six mice) had significantly increased vibration threshold compared to WT mice (n = 29 from six mice), but normal firing activity to a 25 Hz vibration. Trace shows the stimulation applied to the skin, and red squares indicate the time frame during which the number of spikes was calculated. (C) SAM responses to a ramp of 50 Hz vibration with increasing amplitude are similar in WT, *Isl1+/Cre*, and *Isl1+/Cre*::Meis2*LoxP/LoxP* mice. SAM responses to ramp stimuli and their static force responses were also identical in the different genotypes. Fibers from WT and *Isl1+/Cre* mice (n = 5) displayed similar responses. \* p≤0.05; \*\* p≤0.005. Traces show the applied stimulus and red squares the time frame during which the parameters below were measured. (D) Confocal images of Nefh⁺ innervation (green) of CK8⁺ Merkel cells (red) in the forepaw glabrous skin of *Isl1Cre/+* and *Isl1Cre/+*::Meis2*LoxP/LoxP* adult mice. Dotted white squares indicate the close-up of CK8⁺ Merkel cells. Note the lack of Nefh⁺ fibers innervating Merkel cells in mutant mice. White arrows point at contact between NF200⁺ fibers and CK8⁺ Merkel cells. Scale bar = 10 µm. The box plot indicates the percentage of Merkel cells in contact with Nefh⁺ fibers. n = 4. \* p≤0.05 in Mann–Whitney test. (E) Confocal images of Nefh⁺ innervation (green) and CK8⁺ Merkel cells (red) of guard hairs in the hairy back skin of *Isl1Cre/+* and *Isl1Cre/+*::Meis2*LoxP/LoxP* adult mice. Dotted white squares indicate the close-up of CK8⁺ Merkel cells with apparently normal Nefh⁺ innervation. White arrows point at contacts between Nefh⁺ fibers and CK8⁺ Merkel cells. Scale bar = 10 µm. The box plot indicates the percentage of Merkel cells contacted by Nefh⁺ fibers. n = 4. (F) Representative images of whole-mount staining for CK8 in the hairy back skin of WT and *Isl1Cre/+*::Meis2*LoxP/LoxP* E18.5 embryos showing no difference in the number of touch dome between genotypes. Box plots show the number of touch domes per surface area and the number of Merkel cells per touch dome. No significant differences were found between both genotypes in Mann–Whitney test. n = 5 (WT) and 4 (*Isl1+/Cre*::Meis2*LoxP/LoxP*).

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 (Li et al., 2011; Walcher et al., 2018; Shin et al., 2003). They form lanceolate endings, are thinly myelinated, and are activated by movement of the smaller zigzag hairs (Lechner and Lewin, 2013). Consistent with the lack of Meis2 expression in this population reported by scRNAseq databases, Aδ fibers D-hair in the hairy skin showed similar vibration responses in WT and *Isl1+/Cre*::Meis2*LoxP/LoxP* mice (Figure 5—figure supplement 1).

Discussion

The function of the Meis family of TFs in postmitotic neurons has only been marginally addressed (Bouilloux et al., 2016; Jakovcevski et al., 2015; Agoston et al., 2014). Here, we showed that *Meis2* is selectively expressed by subpopulations of early postmitotic cutaneous LTMR and proprioceptive neurons during the 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 scRNAseq analysis and genetic tracing reporting Meis2 in proprioceptive neurons, Aβ field-LTMR, Aβ-SA1-LTMR, and Aβ-RA-LTMR, but not in C-LTMR, Aδ-LTMR, peptidergic, and non-peptidergic nociceptive neurons (Sharma et al., 2020; Usoskin et al., 2015; Zheng et al., 2019). We unambiguously demonstrate that *Meis2* differentially regulates target-field innervation and function of postmitotic LTMR neurons. *Meis2* inactivation in postmitotic 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 postmitotic 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 are consistent with *Meis2* being expressed by Aβ-SA1-LTMR neurons. Interestingly, *Meis2* gene inactivation compromises Merkel cell 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 (Niu et al., 2014; Olson et al., 2016). 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
<figure><img src="image_23.png" /><figcaption>**Figure 5.** *Meis2* gene inactivation affects Meissner corpuscles morphology. (A) Representative images showing S100β⁺ Meissner corpuscles (red) and their innervation by Nefh⁺ fibers (green) in the glabrous skin of WT and *Isl1⁺/Cre::Meis2LoxP/LoxP* adult mice. Scale bar = 10 µm. The box plot shows the average number of times Nefh⁺ fibers cross the midline of the Meissner corpuscles. Dashed blue lines indicate the Meissner corpuscle midline. Blue arrowheads indicate sites where Nefh⁺ fibers cross this midline. (B) RAMs of the glabrous skin exhibited similar vibration threshold and firing activity to a 25 Hz vibration in WT (n = 16 from four mice) and *Isl1⁺/Cre::Meis2LoxP/LoxP* mice (n = 21 from six mice). Glabrous RAMs showed a nonsignificant decrease in firing activity to a ramp of 50 Hz vibration with increasing amplitude in *Isl1⁺/Cre::Meis2LoxP/LoxP* compared to WT littermates, but their response to ramp stimuli was similar in both genotypes. Traces indicate the type of stimulation and red squares the time frame during which the number of spikes was calculated. *** p≤0.001; Student's t-test.</figcaption></figure>

The online version of this article includes the following video and figure supplement(s) for figure 5:

**Figure supplement 1.** Normal electrophysiological responses of D-hair mechanoreceptors following *Meis2* gene inactivation.

**Figure 5—video 1.** Meissner corpuscle in WT.

https://elifesciences.org/articles/89287/figures#fig5video1

**Figure 5—video 2.** Meissner corpuscles in *Isl1⁺/Cre::Meis2LoxP/LoxP*.

https://elifesciences.org/articles/89287/figures#fig5video2

depends on Ntrk1 for survival and innervation, and another Ntrk3/Nefh-positive that does not depend on Ntrk1 signaling during development (Niu et al., 2014). However, the functional significance of these different innervations is unknown. Denervation in rat also pointed to the 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 (Mills et al., 1989).

Because touch domes innervation and Aδ fibers D-hair vibration responses were unaffected in *Isl1⁺/Cre::Meis2LoxP/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
<figure><img src="image_24.png" /><figcaption>**Figure 6.** *Meis2* gene inactivation affects hair follicle innervation and RAM fibers electrophysiological responses in the hairy skin. (A) Representative images of whole-mount immunostaining for Nefh⁺ sensory projections (green) in the hairy skin of adult WT and *Isl1⁺/Cre::Meis2<sup>LoxP/LoxP</sup>* embryos counterstained with S100β (red) to highlight terminal Schwann cells decorating the periphery of hair follicles. Scale bar = 100 µm. (B) Box plots showing the quantification for the number of branch points in the innervation network and the number of innervated hair follicles. n = 3; \*p≤0.05. (C) RAMs in the hairy skin of *Isl1⁺/Cre::Meis2<sup>LoxP/LoxP</sup>* mice (n = 24 from three mice) exhibited significantly increased vibration threshold and reduced firing activity to a 25 Hz vibration compared to WT mice (n = 20 from three mice). RAMs in the hairy skin of *Isl1⁺/Cre::Meis2<sup>LoxP/LoxP</sup>* mice also showed a reduced firing activity in response to a ramp of 50 Hz vibration with increasing amplitude compared to WT and *Isl1⁺/Cre* animals. Fibers recorded from *Isl1⁺/Cre* mice (n = 5) showed similar responses than those recorded from WT mice. \* p≤0.05; \*\*p≤0.005.</figcaption></figure>

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 represents 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 mechanosensitivity can only be identified with a so-called electrical search protocol (Wetzel et al., 2007; Ranade et al., 2014) 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, for example, SA-LTMRs and RA-LTMRs that fail to innervate end organs (Figures 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 (*Neubarth et al., 2020*). 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 A5-LTMR, peptidergic, and non-peptidergic neurons. Surprisingly, although *Meis2* is expressed in proprioceptive neurons (*Usoskin et al., 2015; Zheng et al., 2019; Shin et al., 2020*), 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 (*Shin et al., 2020*). From our data, we could not conclude whether SA-LTMR electrophysiological responses are differentially affected in the glabrous vs 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 was, however, measured. The latter might result from an insufficient number of fibers recording, whereas the first may be due to 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 types into immature neurons (*Guillemot and Hassan, 2017; Baker and Brown, 2018*). 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 (*Hobert, 2016; Hobert and Kratsios, 2019*). The proneural function of *Ngn1* and *Ngn2* genes in neural crest cells, the precursors of DRG sensory neurons, is well demonstrated (*Marmigère and Ernfors, 2007; Ma et al., 1999*), and several terminal selector genes shaping the different DRG sensory subpopulations have also been clearly identified including for cutaneous LTMRs (*Lallemand and Ernfors, 2012; Marmigère and Ernfors, 2007; Abdo et al., 2011; Olson et al., 2016; Wende et al., 2012; Yoshikawa et al., 2013; Lin et al., 1998; Inoue et al., 2002*). *Maf*, *Runx3*, *Shox2*, *ER81*, and *Pea3* are part of this regulatory transcriptional network regulating cutaneous LTMR neuron diversification through intermingled crossed activation and/or repression of subclasses-specific effector genes (*Marmigère and Carroll, 2014*).

In humans, at least 17 different mutations in the *Meis2* gene have been associated with neurodevelopmental delay (*Gangfuß et al., 2021; Giliberti et al., 2020; Shimojima et al., 2017*), emphasizing its essential function in neuronal differentiation. *Meis2* function in late differentiation of postmitotic peripheral sensory neurons adds to the wide actions of this TF in the developing and adult nervous system in the number of regions of the mouse nervous system. Its expression both in dividing neural progenitors, immature neurons, and discrete populations of mature neurons (*Jakovcevski et al., 2015; Toresson et al., 1999; Toresson et al., 2000; Chang and Parrilla, 2016; Frazer et al., 2017; Bumsted-O'Brien et al., 2007; Yan et al., 2020; Yang et al., 2021; Allen et al., 2007*) argues for diverse functions ranging from regulation of neuroblasts cell-cycle exit, to cell-fate decision, neurogenesis, neuronal specification, neurite outgrowth, synaptogenesis, and maintenance of mature neurons. In DRG LTMRs, *Meis2* fulfills some but not all of the criteria defining terminal selector genes. Its inactivation in neural crest cells does not affect sensory neuron generation nor pan-neuronal features, clearly excluding it from a proneural TF function. Although *Meis2* expression is continuously maintained in defined sensory neuron subtypes starting from early postmitotic neurons throughout life, its expression is not restricted to a unique neuronal identity and its early or late inactivation in either sensory neuron progenitors or postmitotic neurons does not influence neuronal subtypes identity nor survival as seen by the normal numbers of *Ntrk2*, *Ntrk3*, or *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 neurotransmitter 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 (Usoskin et al., 2015), suggesting that they exert specific functions in the maintenance of these sensory subclasses, and that their up- or downregulation might affect LTMRs maturation, but did not show any GO term enrichment. Separate GO analyses of our datasets, however, revealed alterations in pathways associated to synapse function and neuron 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-CI cotransporter SLC12A5 associated to GABAergic neurotransmission, but also the glutamate receptor subunits GRID1 and GRIK3 are down- or upregulated 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 (Meltzer et al., 2023). Finally, previous work on *Islet1* conditional deletion in DRG sensory neurons reports considerable changes in gene expression. Early homozygous *Islet1* deletion results in increased sensory neuron apoptosis and a loss of Ntrk1- and Ntrk2-positive neurons, whereas late deletion seems to only affect the nociceptive subpopulations (Sun et al., 2008). Whereas the distal projection defects we report in *Wnt*<sup>Cre</sup>::*Meis2*<sup>LoxP/LoxP</sup> mutant can only be attributed to Meis2 inactivation, it is possible that among the DEGs we identified in *Isl1*<sup>Cre/+</sup>::*Meis2*<sup>LoxP/LoxP</sup>, 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 (Sun et al., 2008; Dykes et al., 2011).

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 (Agoston et al., 2014), 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.

# Materials and methods

## Animals

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

## Mice strains

*Wnt1*<sup>Cre</sup> and *Islet1*<sup>+/Cre</sup> mice were previously described (Yang et al., 2006; Lewis et al., 2013). We previously reported the conditional mutant strain for Meis2 (*Meis2*<sup>LoxP/LoxP</sup>) used in the present study (Roussel et al., 2022). 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) using standard homologous recombination technology in mouse embryonic stem cells. FLP-FRT recombination was used to remove the neomycin selection cassette, and the *Meis2*<sup>LoxP/LoxP</sup> mutant mice were backcrossed for at least eight generations onto the C57BL/6 background before use. The 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 CTB coupled to Cholera Toxin Subunit B conjugated with Alexa Fluor 488 (Thermo Fisher) was injected using a glass
micropipette in several points of the glabrous and hairy forepaw. Mice were sacrificed 7 d after injec-
tion and L4 to L6 DRGs were collected for analysis.

For behavioral assays, skin-nerve preparation, and electrophysiological recording, sex-matched
12-week-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) (Hamburger and Hamilton, 1992) and the right hind limb bud was surgically removed as previously reported (Oakley et al., 1995). Eggs were closed with tape and further grown in the incubator for four (HH27) or seven (HH36) additional days before collection.

Tissue preparation

Mouse and chick embryos were collected at different stages, fixed in 4% paraformaldehyde/PBS over-
night at 4°C, and incubated overnight at 4°C for cryopreservation in increasing sucrose/PBS solutions
(10–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 RNA Nanoprep kit (Stratagene) following the manufacturer's
instruction. Reverse transcription (RT) was carried out for 10 min at 65°C followed by 1 hr at 42°C
and 15 min 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 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'-TTACATAT
AGTGCCACTGCCCATC-3'; gMeis2 forward: 5'-ATGGCGCAAAGGTACGATGAG-3'; gMeis2 reverse:
5'-TTACATGTAGTGCCATTGCCCAT-3'. PCRs were conducted in 50 µl reactions containing 10% RT
product, 200 µM each dNTP, 10 pmol of each primer (MWG-Biotech AG), 3 mM MgCl₂, 6% DMSO,
and 2.5 U of Herculase hotstart DNA polymerase (Stratagene). cDNA was denatured for 10 min at
98°C and amplified for 35 cycles in a three-step program as follows: 1 min denaturation at 98°C,
1 min annealing at annealing temperature and then 1.5 min polymerization at 72°C. PCR products
were separated into 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 (Bouilloux
et al., 2016).

In situ hybridization (ISH)

Before hybridization, slides were air-dried for 2-3 hr at room temperature. Plasmids containing probes
were used to synthesize digoxigenin-labeled or fluorescein-labeled antisense riboprobes according to
the supplier's instructions (Roche) and purified by LiCl precipitation. Sections were hybridized over-
night at 70°C with a solution containing 0.19 M NaCl, 10 mM Tris (pH 7.2), 5 mM NaH₂PO₄*2H₂O/
Na₂HPO₄ (pH 6.8), 50 mM EDTA, 50% formamide, 10% dextran sulfate, 1 mg/ml yeast tRNA, 1x
Denhardt solution, and 100-200 ng/ml of probe. Sections were then washed four times for 20 min
at 65°C in 0.4x SSC pH 7.5, 50% formamide, 0.1% Tween 20, and three times for 20 min at room
temperature in 0.1 M maleic acid, 0.15 M NaCl, and 0.1% Tween 20 (pH 7.5). Sections were blocked
for 1 hr at room temperature in the 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 an NBT/BCIP reaction in 0.1 M Tris-HCl
pH 9.5, 100 mM NaCl, 50 mM MgCl₂, and 0.1% Tween 20.

For double ISH, 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 MgCl₂, 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 alkaline phosphatase was inhibited by incubation in a solution of 0.1 M glycin pH 2.2, 0.2% Tween 20 for 30 min 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 were taken. Wide-field microscopy (Leica DMRB, Germany) was only used for ISH and ISH combined with immunochemistry.

## Immunochemistry

Immunochemistry was performed as previously described (Bouilloux et al., 2016). 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 hr 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. Pictures 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&D Systems, Cat# AF1494); goat anti-TrkC (1/1000; R&D Systems, Cat# AF1404); rabbit anti-TrkA (1/500; Millipore, Cat# 06-574); rabbit anti-parvalbumin antibody (1:500, Swant, Cat# PV 25); guinea pig anti-calcitonin gene-related peptide (CGRP) antibody (1:500, BMA Biomedicals, Cat# T-5053); rabbit anti-Nefh (Sigma, rabbit 1:1000, Cat# N4142), guinea pig anti-c-maf (generous gift of C. Birchmeier, MDC, 1/10,000), rabbit anti-TrpV1 (Sigma, Cat# V2764, 1:1000), mouse-anti-S100β (Sigma, Cat# S2532, 1:1000), rabbit anti-phospho-CREB (Cell Signaling, Cat# 87G3, 1/200, Germany); goat anti-Ret (R&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 (Chang et al., 2014). 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 min for 3–5 hr and kept overnight in the washing solution. The next day, the skin was incubated for 5 d with primary antibodies diluted in PBST containing 5% donkey serum and 20% DMSO. The skin was washed the following day 8–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 min for 6–8 hr 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 (Bouilloux et al., 2016), 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 s at room temperature. After serial washing 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 the 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, Thermo Fisher) 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 follows: denatured libraries were loaded on a HighOutput flowcell kit v2 and sequenced in single-end 84 pb 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) (Smith et al., 2017) default parameters followed by STAR alignment (version 2.7.10) (Dobin et al., 2013) 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) (Anders et al., 2015) was used to aggregate read count per gene followed by differential gene expression analysis with Limma voom on Galaxy (version 3.50.1) (Law et al., 2014). 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 GO 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 (Schwaller et al., 2021). Mice were euthanized by CO₂ 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 MgSO₄, 1.7 mM NaH₂PO₄, 2.0 mM CaCl₂, 9.5 mM sodium gluconate, 5.5 mM glucose, 7.5 mM sucrose, and 10 mM 4-(2-hydroxyethyl)-1-piperazine-etha nesulfonic 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 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 ms⁻¹), Aδ-fibers (1.2–10 ms⁻¹), or Aβ-fibers (>10 ms⁻¹). Mechanical stimulation of the receptive fields of neurons was performed using a piezo actuator (Physik Instrumente, Cat# P-841.60) and a double-ended Nanomotor (Kleindiek Nanotechnik, Cat# MM-NM3108) connected to a force measurement device (Kleindiek Nanotechnik, Cat# 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 mms⁻¹; 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.

## 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 1 hr for two consecutive days. On the testing day, mice were placed in the chamber 1 hr before Von Frey filaments application. The test was performed as previously described (Neubarth et al., 2020). 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 10 times in a row with a break of 30 s following the fifth application. During each application, bend filament was maintained for 4–5 s. 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 5 min. 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 s cutoff was applied to avoid skin damage. 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 2 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 was 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, the 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 three runs met the 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–50 cm s⁻¹. Data are reported as the average of at least three runs per mouse.

### Materials availability

Probes generated for ISH are available upon request to the corresponding author. Meis²LoxP/LoxP mice strain is available upon request to the corresponding author by signing a material transfer agreement (MTA). Sequencing data are available under the following accession code: GSE223788.

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

## Additional information

### Funding

No external funding was received for this work.

### Author contributions

Simon Desiderio, Frederick Schwaller, Formal analysis, Investigation, Methodology; Kevin Tartour, Formal analysis; Kiran Padmanabhan, Funding acquisition, Project administration, Writing – review and editing; Gary R Lewin, Data curation, Formal analysis, Methodology, Writing – review and editing; Patrick Carroll, Funding acquisition, Writing – review and editing; Frederic Marmigere, Conceptualization, Data curation, Formal analysis, Supervision, Funding acquisition, Investigation, Methodology, Writing – original draft, Project administration, Writing – review and editing

### Author ORCIDs

Frederick Schwaller https://orcid.org/0000-0002-7685-5537

Kiran Padmanabhan https://orcid.org/0000-0001-7020-1682

Gary R Lewin https://orcid.org/0000-0002-2890-6352

Patrick Carroll https://orcid.org/0000-0003-3287-8195

Frederic Marmigere https://orcid.org/0000-0002-0515-7483

### Ethics

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 4th of 2020).

### Peer review material

Joint Public Review: https://doi.org/10.7554/eLife.89287.3.sa1

Author Response https://doi.org/10.7554/eLife.89287.3.sa2

## Additional files

### Supplementary files

* MDAR checklist

### Data availability

Sequencing data have been deposited in GEO under accession codes GSE223788.

The following dataset was generated:

<table><thead><tr><td>Author(s)</td><td>Year</td><td>Dataset title</td><td>Dataset URL</td><td>Database and Identifier</td></tr></thead><tbody><tr><td>Tartour K, Carroll P, Padmanabhan K, Marmigère F</td><td>2024</td><td>Touch receptor end-organ innervation and function requires sensory expression of the transcription factor Meis2</td><td>http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE223788</td><td>NCBI Gene Expression Omnibus, GSE223788</td></tr></tbody></table>

## References

Abdo H, Li L, Lallemand F, Bachy I, Xu XJ, Rice FL, Ernfors P. 2011. Dependence on the transcription factor Shox2 for specification of sensory neurons conveying discriminative touch. *The European Journal of Neuroscience* **34**:1529–1541. DOI: https://doi.org/10.1111/j.1460-9568.2011.07883.x, PMID: 22103411
Abdollahi MR, Morrison E, Sirey T, Molnar Z, Hayward BE, Carr IM, Springell K, Woods CG, Ahmed M, Hattingh L, Corry P, Pilz DT, Stoodley N, Crow Y, Taylor GR, Bonthron DT, Sheridan E. 2009. Mutation of the Variant α-Tubulin TUBA8 results in polymicrogyria with optic nerve hypoplasia. *The American Journal of Human Genetics* **85**:737–744. DOI: https://doi.org/10.1016/j.ajhg.2009.10.007

Abraira VE, Ginty DD. 2013. The sensory neurons of touch. *Neuron* **79**:618–639. DOI: https://doi.org/10.1016/j.neuron.2013.07.051, PMID: 23972592

Agoston Z, Heine P, Brill MS, Grebbin BM, Hau AC, Kallenborn-Gerhardt W, Schramm J, Götz M, Schulte D. 2014. Meis2 is a Pax6 co-factor in neurogenesis and dopaminergic periglomerular fate specification in the adult olfactory bulb. *Development* **141**:28–38. DOI: https://doi.org/10.1242/dev.097295

Allen ZJ, Waclaw RR, Colbert MC, Campbell K. 2007. Molecular identity of olfactory bulb interneurons: transcriptional codes of periglomerular neuron subtypes. *Journal of Molecular Histology* **38**:517–525. DOI: https://doi.org/10.1007/s10735-007-9115-4, PMID: 17624499

Anders S, Pyl PT, Huber W. 2015. HTSeq--a Python framework to work with high-throughput sequencing data. *Bioinformatics* **31**:166–169. DOI: https://doi.org/10.1093/bioinformatics/btu638, PMID: 25260700

Arber S, Ladle DR, Lin JH, Frank E, Jessell TM. 2000. ETS gene Er81 controls the formation of functional connections between group Ia sensory afferents and motor neurons. *Cell* **101**:485–498. DOI: https://doi.org/10.1016/s0092-8674(00)80859-4, PMID: 10850491

Bae T, Fasching L, Wang Y, Shin JH, Suvakov M, Jang Y, Norton S, Dias C, Mariani J, Jourdon A, Wu F, Panda A, Pattni R, Chahine Y, Yeh R, Roberts RC, Huttner A, Kleinman JE, Hyde TM, Straub RE, et al. 2022. Analysis of somatic mutations in 131 human brains reveals aging-associated hypermutability. *Science* **377**:511–517. DOI: https://doi.org/10.1126/science.abm6222

Bai L, Lehnert BP, Liu J, Neubarth NL, Dickendesher TL, Nwe PH, Cassidy C, Woodbury CJ, Ginty DD. 2015. Genetic identification of an expansive mechanoreceptor sensitive to skin stroking. *Cell* **163**:1783–1795. DOI: https://doi.org/10.1016/j.cell.2015.11.060

Baker NE, Brown NL. 2018. All in the family: proneural bHLH genes and neuronal diversity. *Development* **145**:dev159426. DOI: https://doi.org/10.1242/dev.159426, PMID: 29720483

Bouilloux F, Thireau J, Ventéo S, Farah C, Karam S, Dauvilliers Y, Valmier J, Copeland NG, Jenkins NA, Richard S, Marmigère F. 2016. Loss of the transcription factor Meis1 prevents sympathetic neurons target-field innervation and increases susceptibility to sudden cardiac death. *eLife* **5**:11627. DOI: https://doi.org/10.7554/eLife.11627

Bourane S, Grossmann KS, Britz O, Dalet A, Del Barrio MG, Stam FJ, Garcia-Campmany L, Koch S, Goulding M. 2015. Identification of a spinal circuit for light touch and fine motor control. *Cell* **160**:503–515. DOI: https://doi.org/10.1016/j.cell.2015.01.011

Bumsted-O'Brien KM, Hendrickson A, Haverkamp S, Ashery-Padan R, Schulte D. 2007. Expression of the homeodomain transcription factor Meis2 in the embryonic and postnatal retina. *The Journal of Comparative Neurology* **505**:58–72. DOI: https://doi.org/10.1002/cne.21458, PMID: 17729288

Burgess PR, Horch KW. 1973. Specific regeneration of cutaneous fibers in the cat. *Journal of Neurophysiology* **36**:101–114. DOI: https://doi.org/10.1152/jn.1973.36.1.101

Calderó J, Prevette D, Mei X, Oakley RA, Li L, Milligan C, Houenou L, Burek M, Oppenheim RW. 1998. Peripheral target regulation of the development and survival of spinal sensory and motor neurons in the chick embryo. *The Journal of Neuroscience* **18**:356–370. DOI: https://doi.org/10.1523/JNEUROSCI.18-01-00356.1998, PMID: 9412513

Catela C, Shin MM, Lee DH, Liu JP, Dasen JS. 2016. Hox proteins coordinate motor neuron differentiation and connectivity programs through Ret/Gfrα Genes. *Cell Reports* **14**:1901–1915. DOI: https://doi.org/10.1016/j.celrep.2016.01.067, PMID: 26904955

Chang H, Wang Y, Wu H, Nathans J. 2014. Flat mount imaging of mouse skin and its application to the analysis of hair follicle patterning and sensory axon morphology. *Journal of Visualized Experiments* **1**:e51749. DOI: https://doi.org/10.3791/51749, PMID: 24999071

Chang I, Parrilla M. 2016. Expression patterns of homeobox genes in the mouse vomeronasal organ at postnatal stages. *Gene Expression Patterns* **21**:69–80. DOI: https://doi.org/10.1016/j.gep.2016.08.001, PMID: 27521061

Dasen JS, Tice BC, Brenner-Morton S, Jessell TM. 2005. A Hox regulatory network establishes motor neuron pool identity and target-muscle connectivity. *Cell* **123**:477–491. DOI: https://doi.org/10.1016/j.cell.2005.09.009, PMID: 16269338

Delmas P, Hao J, Rodat-Despoix L. 2011. Molecular mechanisms of mechanotransduction in mammalian sensory neurons. *Nature Reviews. Neuroscience* **12**:139–153. DOI: https://doi.org/10.1038/nrn2993, PMID: 21304548

de Nooij JC, Doobar S, Jessell TM. 2013. Etv1 inactivation reveals proprioceptor subclasses that reflect the level of NT3 Expression in muscle targets. *Neuron* **77**:1055–1068. DOI: https://doi.org/10.1016/j.neuron.2013.01.015

Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, Batut P, Chaisson M, Gingeras TR. 2013. STAR: ultrafast universal RNA-seq aligner. *Bioinformatics* **29**:15–21. DOI: https://doi.org/10.1093/bioinformatics/bts635, PMID: 23104886

Dykes IM, Tempest L, Lee SI, Turner EE. 2011. Brn3a and Islet1 act epistatically to regulate the gene expression program of sensory differentiation. *Journal of Neuroscience* **31**:9789–9799. DOI: https://doi.org/10.1523/JNEUROSCI.0901-11.2011

Ericson J, Thor S, Edlund T, Jessell TM, Yamada T. 1992. Early stages of motor neuron differentiation revealed by expression of homeobox gene Islet-1. *Science* **256**:1555–1560. DOI: https://doi.org/10.1126/science.1350865, PMID: 1350865
Frazer S, Prados J, Niquille M, Cadilhac C, Markopoulos F, Gomez L, Tomasello U, Telley L, Holtmaat A, Jabaudon D, Dayer A. 2017. Transcriptomic and anatomic parcellation of 5-HT<sub>3A</sub>R expressing cortical interneuron subtypes revealed by single-cell RNA sequencing. *Nature Communications* **8**:14219. DOI: https://doi.org/10.1038/ncomms14219, PMID: 28134272

Gangfuß A, Yigit G, Altmüller J, Nürnberg P, Czeschik JC, Wollnik B, Bögershausen N, Burfeind P, Wieczorek D, Kaiser F, Roos A, Kölbel H, Schara-Schmidt U, Kuechler A. 2021. Intellectual disability associated with craniofacial dysmorphism, cleft palate, and congenital heart defect due to A de novo MEIS2 mutation: A clinical longitudinal study. *American Journal of Medical Genetics. Part A* **185**:1216–1221. DOI: https://doi.org/10.1002/ajmg.a.62070, PMID: 33427397

Geerts D, Schilderink N, Jorritsma G, Versteeg R. 2003. The role of the MEIS homeobox genes in neuroblastoma. *Cancer Letters* **197**:87–92. DOI: https://doi.org/10.1016/s0304-3835(03)00087-9, PMID: 12880965

Giliberti A, Currò A, Papa FT, Frullanti E, Ariani F, Coriolani G, Grosso S, Renieri A, Mari F. 2020. MEIS2 gene is responsible for intellectual disability, cardiac defects and a distinct facial phenotype. *European Journal of Medical Genetics* **63**:103627. DOI: https://doi.org/10.1016/j.ejmg.2019.01.017, PMID: 30735726

Guillemot F, Hassan BA. 2017. Beyond proneural: emerging functions and regulations of proneural proteins. *Current Opinion in Neurobiology* **42**:93–101. DOI: https://doi.org/10.1016/j.conb.2016.11.011, PMID: 28025176

Hamburger V, Hamilton HL. 1992. A series of normal stages in the development of the chick embryo. 1951. *Developmental Dynamics* **195**:231–272. DOI: https://doi.org/10.1002/aja.1001950404, PMID: 1304821

Handler A, Ginty DD. 2021. The mechanosensory neurons of touch and their mechanisms of activation. *Nature Reviews. Neuroscience* **22**:521–537. DOI: https://doi.org/10.1038/s41583-021-00489-x, PMID: 34312536

Hobert O. 2016. Terminal selectors of neuronal identity. Rougvie AE (Ed). *Current Topics in Developmental Biology*. Academic Press Inc. p. 455–475.

Hobert O, Kratsios P. 2019. Neuronal identity control by terminal selectors in worms, flies, and chordates. *Current Opinion in Neurobiology* **56**:97–105. DOI: https://doi.org/10.1016/j.conb.2018.12.006, PMID: 30665084

Hu J, Huang T, Li T, Guo Z, Cheng L. 2012. C-Maf is required for the development of dorsal horn laminae III/IV neurons and mechanoreceptive DRG axon projections. *The Journal of Neuroscience* **32**:5362–5373. DOI: https://doi.org/10.1523/JNEUROSCI.6239-11.2012, PMID: 22514301

Inoue K, Ozaki S, Shiga T, Ito K, Masuda T, Okado N, Iseda T, Kawaguchi S, Ogawa M, Bae S-C, Yamashita N, Itohara S, Kudo N, Ito Y. 2002. Runx3 controls the axonal projection of proprioceptive dorsal root ganglion neurons. *Nature Neuroscience* **5**:946–954. DOI: https://doi.org/10.1038/nn925, PMID: 12352981

Jakovcevski M, Ruan H, Shen EY, Dincer A, Javidfar B, Ma Q, Peter CJ, Cheung I, Mitchell AC, Jiang Y, Lin CL, Pothula V, Stewart AF, Ernst P, Yao WD, Akbarian S. 2015. Neuronal Kmt2a/Mll1 histone methyltransferase is essential for prefrontal synaptic plasticity and working memory. *The Journal of Neuroscience* **35**:5097–5108. DOI: https://doi.org/10.1523/JNEUROSCI.3004-14.2015, PMID: 25834037

Lallemend F, Ernfors P. 2012. Molecular interactions underlying the specification of sensory neurons. *Trends in Neurosciences* **35**:373–381. DOI: https://doi.org/10.1016/j.tins.2012.03.006, PMID: 22516617

Law CW, Chen Y, Shi W, Smyth GK. 2014. voom: Precision weights unlock linear model analysis tools for RNA-seq read counts. *Genome Biology* **15**:R29. DOI: https://doi.org/10.1186/gb-2014-15-2-r29, PMID: 24485249

Lechner SG, Lewin GR. 2013. Hairy sensation. *Physiology* **28**:142–150. DOI: https://doi.org/10.1152/physiol.00059.2012, PMID: 23636260

Lewin GR, McMahon SB. 1991. Physiological properties of primary sensory neurons appropriately and inappropriately innervating skin in the adult rat. *Journal of Neurophysiology* **66**:1205–1217. DOI: https://doi.org/10.1152/jn.1991.66.4.1205

Lewis AE, Vasudevan HN, O'Neill AK, Soriano P, Bush JO. 2013. The widely used Wnt1-Cre transgene causes developmental phenotypes by ectopic activation of Wnt signaling. *Developmental Biology* **379**:229–234. DOI: https://doi.org/10.1016/j.ydbio.2013.04.026

Li L, Rutlin M, Abraira VE, Cassidy C, Kus L, Gong S, Jankowski MP, Luo W, Heintz N, Koerber HR, Woodbury CJ, Ginty DD. 2011. The functional organization of cutaneous low-threshold mechanosensory neurons. *Cell* **147**:1615–1627. DOI: https://doi.org/10.1016/j.cell.2011.11.027, PMID: 22196735

Lin JH, Saito T, Anderson DJ, Lance-Jones C, Jessell TM, Arber S. 1998. Functionally related motor neuron pool and muscle sensory afferent subtypes defined by coordinate ETS gene expression. *Cell* **95**:393–407. DOI: https://doi.org/10.1016/s0092-8674(00)81770-5, PMID: 9814709

Longobardi E, Penkov D, Mateos D, De Florian G, Torres M, Blasi F. 2014. Biochemistry of the tale transcription factors PREP, MEIS, and PBX in vertebrates. *Developmental Dynamics* **243**:59–75. DOI: https://doi.org/10.1002/dvdy.24016, PMID: 23873833

Ma Q, Fode C, Guillemot F, Anderson DJ. 1999. NEUROGENIN1 and NEUROGENIN2 control two distinct waves of neurogenesis in developing dorsal root ganglia. *Genes & Development* **13**:1717–1728. DOI: https://doi.org/10.1101/gad.13.13.1717, PMID: 10398684

Machon O, Masek J, Machonova O, Krauss S, Kozmik Z. 2015. Meis2 is essential for cranial and cardiac neural crest development. *BMC Developmental Biology* **15**:40. DOI: https://doi.org/10.1186/s12861-015-0093-6, PMID: 26545946

Marmigère F, Ernfors P. 2007. Specification and connectivity of neuronal subtypes in the sensory lineage. *Nature Reviews. Neuroscience* **8**:114–127. DOI: https://doi.org/10.1038/nrn2057, PMID: 17237804
Marmigère F, Carroll P. 2014. Neurotrophin signalling and transcription programmes interactions in the development of Somatosensory neurons. Kulkarni SK (Ed). *Handbook of Experimental Pharmacology*. Springer. p. 329–353. DOI: https://doi.org/10.1007/978-3-642-45106-5_13, PMID: 24668479

Meltzer S, Boulanger KC, Chirila AM, Osei-Asante E, DeLisle M, Zhang Q, Kalish BT, Tasnim A, Huey EL, Fuller LC, Flaherty EK, Maniatis T, Garrett AM, Weiner JA, Ginty DD. 2023. γ-Protocadherins control synapse formation and peripheral branching of touch sensory neurons. *Neuron* **111**:1776–1794. DOI: https://doi.org/10.1016/j.neuron.2023.03.012, PMID: 37028432

Mills LR, Nurse CA, Diamond J. 1989. The neural dependency of Merkel cell development in the rat: the touch domes and foot pads contrasted. *Developmental Biology* **136**:61–74. DOI: https://doi.org/10.1016/0012-1606(89)90130-9, PMID: 2806725

Neubarth NL, Emanuel AJ, Liu Y, Springel MW, Handler A, Zhang Q, Lehnert BP, Guo C, Orefice LL, Abdelaziz A, DeLisle MM, Iskols M, Rhyins J, Kim SJ, Cattel SJ, Regehr W, Harvey CD, Drugowitsch J, Ginty DD. 2020. Meissner corpuscles and their spatially intermingled afferents underlie gentle touch perception. *Science* **368**:eabb2751. DOI: https://doi.org/10.1126/science.abb2751, PMID: 32554568

Niu J, Vysochan A, Luo W. 2014. Dual innervation of neonatal merkel cells in mouse touch domes. *PLOS ONE* **9**:e92027. DOI: https://doi.org/10.1371/journal.pone.0092027

Oakley RA, Garner AS, Large TH, Frank E. 1995. Muscle sensory neurons require neurotrophin-3 from peripheral tissues during the period of normal cell death. *Development* **121**:1341–1350. DOI: https://doi.org/10.1242/dev.121.5.1341, PMID: 7789265

Oakley RA, Lefcort FB, Clary DO, Reichardt LF, Prevette D, Oppenheim RW, Frank E. 1997. Neurotrophin-3 promotes the differentiation of muscle spindle afferents in the absence of peripheral targets. *The Journal of Neuroscience* **17**:4262–4274. DOI: https://doi.org/10.1523/JNEUROSCI.17-11-04262.1997, PMID: 9151743

Olson W, Dong P, Fleming M, Luo W. 2016. The specification and wiring of mammalian cutaneous low-threshold mechanoreceptors. *Wiley Interdisciplinary Reviews. Developmental Biology* **5**:389–404. DOI: https://doi.org/10.1002/wdev.229, PMID: 26992078

Patel TD, Kramer I, Kucera J, Niederkofler V, Jessell TM, Arber S, Snider WD. 2003. Peripheral NT3 signaling is required for ETS protein expression and central patterning of proprioceptive sensory afferents. *Neuron* **38**:403–416. DOI: https://doi.org/10.1016/S0896-6273(03)00261-7

Poliak S, Norovich AL, Yamagata M, Sanes JR, Jessell TM. 2016. Muscle-type identity of proprioceptors specified by spatially restricted signals from limb mesenchyme. *Cell* **164**:512–525. DOI: https://doi.org/10.1016/j.cell.2015.12.049, PMID: 26824659

Pulman J, Ruzzenente B, Bianchi L, Rio M, Boddaert N, Munnich A, Rötig A, Metodiev MD. 2019. Mutations in the MRPS28 gene encoding the small mitoribosomal subunit protein bS1m in a patient with intrauterine growth retardation, craniofacial dysmorphism and multisystemic involvement. *Human Molecular Genetics* **28**:1445–1462. DOI: https://doi.org/10.1093/hmg/ddy441, PMID: 30566640

Ranade SS, Woo SH, Dubin AE, Moshourab RA, Wetzel C, Petrus M, Mathur J, Bégay V, Coste B, Mainquist J, Wilson AJ, Francisco AG, Reddy K, Qiu Z, Wood JN, Lewin GR, Patapoutian A. 2014. Piezo2 is the major transducer of mechanical forces for touch sensation in mice. *Nature* **516**:121–125. DOI: https://doi.org/10.1038/nature13980

Rifkin JT, Todd VJ, Anderson LW, Lefcort F. 2000. Dynamic expression of neurotrophin receptors during sensory neuron genesis and differentiation. *Developmental Biology* **227**:465–480. DOI: https://doi.org/10.1006/dbio.2000.9841, PMID: 11071767

Roussel J, Larcher R, Sicard P, Bideaux P, Richard S, Marmigère F, Thireau J. 2022. The autism-associated Meis2 gene is necessary for cardiac baroreflex regulation in mice. *Scientific Reports* **12**:20150. DOI: https://doi.org/10.1038/s41598-022-24616-5, PMID: 36418415

Schwaller F, Bégay V, García-García G, Taberner FJ, Moshourab R, McDonald B, Docter T, Kühnemund J, Ojeda-Alonso J, Paricio-Montesinos R, Lechner SG, Poulet JFA, Millan JM, Lewin GR. 2021. USH2A is a Meissner's corpuscle protein necessary for normal vibration sensing in mice and humans. *Nature Neuroscience* **24**:74–81. DOI: https://doi.org/10.1038/s41593-020-00751-y, PMID: 33288907

Scott A, Hasegawa H, Sakurai K, Yaron A, Cobb J, Wang F. 2011. Transcription factor short stature homeobox 2 is required for proper development of tropomyosin-related kinase B-expressing mechanosensory neurons. *The Journal of Neuroscience* **31**:6741–6749. DOI: https://doi.org/10.1523/JNEUROSCI.5883-10.2011, PMID: 21543603

Sharma N, Flaherty K, Lezgiyeva K, Wagner DE, Klein AM, Ginty DD. 2020. The emergence of transcriptional identity in somatosensory neurons. *Nature* **577**:392–398. DOI: https://doi.org/10.1038/s41586-019-1900-1, PMID: 31915380

Shimojima K, Ondo Y, Okamoto N, Yamamoto T. 2017. A 15q14 microdeletion involving MEIS2 identified in A patient with autism spectrum disorder. *Human Genome Variation* **4**:17029. DOI: https://doi.org/10.1038/hgv.2017.29, PMID: 28736618

Shin JB, Martinez-Salgado C, Heppenstall PA, Lewin GR. 2003. A T-type calcium channel required for normal function of A mammalian mechanoreceptor. *Nature Neuroscience* **6**:724–730. DOI: https://doi.org/10.1038/nn1076, PMID: 12808460

Shin MM, Catela C, Dasen J. 2020. Intrinsic control of neuronal diversity and synaptic specificity in a proprioceptive circuit. *eLife* **9**:e56374. DOI: https://doi.org/10.7554/eLife.56374, PMID: 32808924

Smith T, Heger A, Sudbery I. 2017. UMI-tools: modeling sequencing errors in Unique Molecular Identifiers to improve quantification accuracy. *Genome Research* **27**:491–499. DOI: https://doi.org/10.1101/gr.209601.116, PMID: 28100584
Sun Y, Dykes IM, Liang X, Eng SR, Evans SM, Turner EE. 2008. A central role for Islet1 in sensory neuron development linking sensory and spinal gene regulatory programs. *Nature Neuroscience* **11**:1283-1293. DOI: https://doi.org/10.1038/nn.2209, PMID: 18849985

Toresson H, Mata de Urquiza A, Fagerström C, Perlmann T, Campbell K. 1999. Retinoids are produced by glia in the lateral ganglionic eminence and regulate striatal neuron differentiation. *Development* **126**:1317-1326. DOI: https://doi.org/10.1242/dev.126.6.1317, PMID: 10021349

Toresson H, Parmar M, Campbell K. 2000. Expression of Meis and Pbx genes and their protein products in the developing telencephalon: implications for regional differentiation. *Mechanisms of Development* **94**:183-187. DOI: https://doi.org/10.1016/s0925-4773(00)00324-5, PMID: 10842069

Usoskin D, Furlan A, Islam S, Abdo H, Lönnerberg P, Lou D, Hjerling-Leffler J, Haeggström J, Kharchenko O, Kharchenko PV, Linnarsson S, Ernfors P. 2015. Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencing. *Nature Neuroscience* **18**:145-153. DOI: https://doi.org/10.1038/nn.3881, PMID: 25420068

Vermeiren S, Bellefroid EJ, Desiderio S. 2020. Vertebrate sensory ganglia: common and divergent features of the transcriptional programs generating their functional specialization. *Frontiers in Cell and Developmental Biology* **8**:587699. DOI: https://doi.org/10.3389/fcell.2020.587699, PMID: 33195244

Walcher J, Ojeda-Alonso J, Haseleu J, Oosthuizen MK, Rowe AH, Bennett NC, Lewin GR. 2018. Specialized mechanoreceptor systems in rodent glabrous skin. *The Journal of Physiology* **596**:4995-5016. DOI: https://doi.org/10.1113/JP276608, PMID: 30132906

Wang Y, Wu H, Zelenin P, Fontanet P, Wanderoy S, Petitpré C, Comai G, Bellardita C, Xue-Franzén Y, Huettl R-E, Huber AB, Tajbakhsh S, Kiehn O, Ernfors P, Deliagina TG, Lallemand F, Hadjab S. 2019. Muscle-selective RUNX3 dependence of sensorimotor circuit development. *Development* **146**:dev181750. DOI: https://doi.org/10.1242/dev.181750, PMID: 31575648

Wende H, Lechner SG, Cheret C, Bourane S, Kolanczyk ME, Pattyn A, Reuter K, Munier FL, Carroll P, Lewin GR, Birchmeier C. 2012. The transcription factor c-Maf controls touch receptor development and function. *Science* **335**:1373-1376. DOI: https://doi.org/10.1126/science.1214314, PMID: 22345400

Wetzel C, Hu J, Riethmacher D, Benckendorff A, Harder L, Eilers A, Moshourab R, Kozlenkov A, Labuz D, Caspani O, Erdmann B, Machelska H, Heppenstall PA, Lewin GR. 2007. A stomatin-domain protein essential for touch sensation in the mouse. *Nature* **445**:206-209. DOI: https://doi.org/10.1038/nature05394, PMID: 17167420

Wu H, Petitpré C, Fontanet P, Sharma A, Bellardita C, Quadros RM, Jannig PR, Wang Y, Heimel JA, Cheung KKY, Wanderoy S, Xuan Y, Meletis K, Ruas J, Gurumurthy CB, Kiehn O, Hadjab S, Lallemand F. 2021. Distinct subtypes of proprioceptive dorsal root ganglion neurons regulate adaptive proprioception in mice. *Nature Communications* **12**:1026. DOI: https://doi.org/10.1038/s41467-021-21173-9, PMID: 33589589

Yan W, Laboulaye MA, Tran NM, Whitney IE, Benhar I, Sanes JR. 2020. Mouse retinal cell atlas: molecular identification of over sixty amacrine cell types. *The Journal of Neuroscience* **40**:5177-5195. DOI: https://doi.org/10.1523/JNEUROSCI.0471-20.2020, PMID: 32457074

Yang L, Cai CL, Lin L, Qyang Y, Chung C, Monteiro RM, Mummery CL, Fishman GI, Cogen A, Evans S. 2006. Isl1Cre reveals a common Bmp pathway in heart and limb development. *Development* **133**:1575-1585. DOI: https://doi.org/10.1242/dev.02322, PMID: 16556916

Yang L, Su Z, Wang Z, Li Z, Shang Z, Du H, Liu G, Qi D, Yang Z, Xu Z, Zhang Z. 2021. Transcriptional profiling reveals the transcription factor networks regulating the survival of striatal neurons. *Cell Death & Disease* **12**:262. DOI: https://doi.org/10.1038/s41419-021-03552-8, PMID: 33712552

Yoshikawa M, Murakami Y, Senzaki K, Masuda T, Ozaki S, Ito Y, Shiga T. 2013. Coexpression of Runx1 and Runx3 in mechanoreceptive dorsal root ganglion neurons. *Developmental Neurobiology* **73**:469-479. DOI: https://doi.org/10.1002/dneu.22073, PMID: 23378040

Zheng Y, Liu P, Bai L, Trimmer JS, Bean BP, Ginty DD. 2019. Deep sequencing of somatosensory neurons reveals molecular determinants of intrinsic physiological properties. *Neuron* **103**:598-616. DOI: https://doi.org/10.1016/j.neuron.2019.05.039, PMID: 31248728

Zimmerman A, Bai L, Ginty DD. 2014. The gentle touch receptors of mammalian skin. *Science* **346**:950-954. DOI: https://doi.org/10.1126/science.1254229, PMID: 25414303