# Muscle-selective RUNX3 dependence of sensorimotor circuit development

Yiqiao Wang, Haohao Wu, Pavel Zelenin, Paula Fontanet, Simone Wanderoy, Charles Petitpré, Glenda Comai, Carmelo Bellardita, Yongtao Xue-Franzén, Rosa-Eva Huettl, et al.

► To cite this version:

Yiqiao Wang, Haohao Wu, Pavel Zelenin, Paula Fontanet, Simone Wanderoy, et al. Muscle-selective RUNX3 dependence of sensorimotor circuit development. Development (Cambridge, England), 2019, 146 (20), pp.dev181750. 10.1242/dev.181750. hal-02367269

HAL Id: hal-02367269

https://hal.science/hal-02367269

Submitted on 18 Nov 2019

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
RESEARCH ARTICLE

# Muscle-selective RUNX3 dependence of sensorimotor circuit development

Yiqiao Wang<fnref n="1" />, Haohao Wu<fnref n="1" />, Pavel Zelenin<fnref n="1" />, Paula Fontanet<fnref n="1" />, Simone Wanderoy<fnref n="1" />, Charles Petitpré<fnref n="1" />, Glenda Comai<fnref n="2" />, Carmelo Bellardita³, Yongtao Xue-Franzén<fnref n="1" />, Rosa-Eva Huettl⁴,*, Andrea B. Huber⁴,‡, Shahragim Tajbakhsh<fnref n="2" />, Ole Kiehn<fnref n="1" />,³, Patrik Ernfors⁵, Tatiana G. Deliagina<fnref n="1" />, François Lallemand<fnref n="1" />,⁶,§,¶ and Saida Hadjab<fnref n="1" />,§,¶

ABSTRACT

The control of all our motor outputs requires constant monitoring by proprioceptive sensory neurons (PSNs) that convey continuous muscle sensory inputs to the spinal motor network. Yet the molecular programs that control the establishment of this sensorimotor circuit remain largely unknown. The transcription factor RUNX3 is essential for the early steps of PSNs differentiation, making it difficult to study its role during later aspects of PSNs specification. Here, we conditionally inactivate *Runx3* in PSNs after peripheral innervation and identify that RUNX3 is necessary for maintenance of cell identity of only a subgroup of PSNs, without discernable cell death. RUNX3 also controls the sensorimotor connection between PSNs and motor neurons at limb level, with muscle-by-muscle variable sensitivities to the loss of *Runx3* that correlate with levels of RUNX3 in PSNs. Finally, we find that muscles and neurotrophin 3 signaling are necessary for maintenance of RUNX3 expression in PSNs. Hence, a transcriptional regulator that is crucial for specifying a generic PSN type identity after neurogenesis is later regulated by target muscle-derived signals to contribute to the specialized aspects of the sensorimotor connection selectivity.

KEY WORDS: Sensory system, Sensorimotor circuit, Dorsal root ganglia, Neuronal specification, Neurotrophins

INTRODUCTION

The neuromuscular circuitry that controls all body movements relies on constant sensory feedback from the periphery to coordinate its

commands to hundreds of muscles. The sensory components of this feedback are the proprioceptive sensory neurons (PSNs) of the dorsal root ganglia (DRG), which convey information from individual muscles to specific neuron groups in the spinal cord. Previous studies have indicated that the basic design of this sensorimotor circuit is already established at birth (Mears and Frank, 1997), and that its construction is largely independent of patterned neuronal activity (Frank, 1990; Mendelson and Frank, 1991; Mendelsohn et al., 2015), implying differential recruitment of specific molecular pathways during the establishment of the sensorimotor connections.

PSNs are identified by their specific co-expression of the runt related transcription factor RUNX3, tropomyosin receptor kinase C (TRKC, receptor for neurotrophin 3, NT3; also known as Ntrk3), the ETS transcription factor ER81, parvalbumin (PV) and the vesicular glutamate transporter 1 (VGLUT1) (Oliveira et al., 2003; Lallemand and Ernfors, 2012). Peripherally, they terminate in muscles and innervate the Golgi tendon organ (GTO) (Ib PSNs) and muscle spindles (MSs) (Ia and II PSNs). Centrally, their axons connect with distinct classes of interneurons and α-motor neurons (MNs) in the deep dorsal and the ventral horn of the spinal cord (Lallemand and Ernfors, 2012). Only Ia proprioceptive afferents make monosynaptic connections with MNs, establishing the sensorimotor reflex arc. This circuit is very selective as Ia afferents connect to MNs supplying the same muscles and avoid making connections with MNs commanding antagonistic muscles (Eccles et al., 1957; Frank and Mendelson, 1990). Most studies aimed at understanding the molecular programs that control the development of the motor circuits in relation to the anatomical identity of the peripheral connections have focused on the spinal MN specification (Dasen, 2009; Dasen et al., 2005; Price et al., 2002). Spinal MNs innervating the limbs are organized into distinct anatomical columns, the identity of which is tightly controlled by intrinsic genetic programs during early stages of MNs development, before they project peripherally into the limb, and which are essential for the construction of the stereotypic connections with individual muscles (Dasen, 2009; Dasen et al., 2005; Price et al., 2002; De Marco Garcia and Jessell, 2008). Unlike MNs, PSNs innervating a specific muscle are scattered in a mosaic fashion throughout multiple DRGs (Honig, 1982; Honig et al., 1998), and do not seem to possess any specific identity other than a generic expression of genes common to all PSNs before innervating their muscle target (Wu et al., 2019). PSNs instead would acquire subclass identities through extrinsic, presumably target-derived, signals days after their neurogenesis (Wu et al., 2019; Poliak et al., 2016). Hence, genetic manipulation of signaling pathways or transcription factors affecting the afferent outgrowth and muscle targeting of PSNs often results in a lack of sensorimotor connections

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

<footnote n="1">¹Department of Neuroscience, Karolinska Institutet, Stockholm 17177, Sweden.</footnote>

<footnote n="2">²Department of Developmental and Stem Cell Biology, Institut Pasteur, CNRS UMR3738, Paris 75015, France. ³Department of Neuroscience, University of Copenhagen, Copenhagen 2200, Denmark. ⁴Helmholtz Zentrum München, German Research Center for Environmental Health, Institute of Developmental Genetics, Neuherberg 85764, Germany. ⁵Unit of Molecular Neurobiology, Department of Medical Biochemistry and Biophysics, Karolinska Institutet, Stockholm 17177, Sweden. ⁶Ming Wai Lau Centre for Reparative Medicine, Stockholm node, Karolinska Institutet, Stockholm 17177, Sweden.</footnote>

<footnote>*Present address: Department of Stress Neurobiology and Neurogenetics, Max-Planck-Institute of Psychiatry, Munich 80804, Germany. †Present address: ETH Zurich, Department of Biosystems Science and Engineering, Basel 4058, Switzerland.</footnote>

<footnote>§Co-senior authors</footnote>

<footnote>¶Authors for correspondence (francois.lallemand@ki.se, saida.hadjab@ki.se)</footnote>

<footnote>F.L., 0000-0001-5484-0011; S.H., 0000-0001-7953-8396</footnote>

<footnote>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.</footnote>
(Lallemand et al., 2012; Levanon et al., 2002; De Nooij et al., 2013; Patel et al., 2003). For example, deletion of RUNX3, which drives the specification of PSNs and is associated with a complete absence of muscle proprioceptive innervation (Lallemand et al., 2012), results in a large deficit of central innervation (Nakamura et al., 2008). Similarly, without NT3-TRKC signaling, PSNs fail to innervate their peripheral muscle targets and their central projections do not extend further into the ventral horn of the spinal cord (Patel et al., 2003). Earlier studies in chick and a more recent study in mice provide strong support for a role of the target in assigning muscle-specific identities of PSNs that are likely to play a role in the establishment of specific patterns of central sensorimotor connections (Poliak et al., 2016; Wenner and Frank, 1995). However, with the exception of a role for Sema3E-PlexinD1 signaling in gating one specialized aspect of sensorimotor connectivity (Fukuhara et al., 2013; Pecho-Vrieseling et al., 2009), little is known about the molecular mechanisms required for the establishment of sensorimotor connections between selective PSNs and the central motor network. Here, using an experimental strategy that depletes RUNX3 expression after peripheral innervation, we uncovered a transcriptional link between peripheral muscle target, likely involving NT3 levels, and selective sensorimotor connectivity during late fetal stages.

## RESULTS

### RUNX3 is essential for the maintenance of cell identity of a subgroup of PSNs

To investigate the role of RUNX3 on the development of PSNs and of the sensorimotor circuits independently of its function on neuronal survival and early aspects of peripheral innervation – which end at E12.5 and E14.0, respectively (Fig. 1A) (Fariñas et al., 1998; Tourtellotte et al., 2001) – we first generated *PV<sup>Cre</sup>;Runx3<sup>fl/fl</sup>* mice, with loxP-flanked *Runx3* alleles and Cre expressed in PV-expressing neurons. The *PV<sup>Cre</sup>;Runx3<sup>fl/fl</sup>* mouse line induced recombination in a large majority of PV<sup>+</sup> and RUNX3<sup>+</sup> PSNs at E16.5 when analyzing brachial DRG from *PV<sup>Cre</sup>;R26<sup>tdTOM</sup>* mice (Fig. S1A,B). However, in *PV<sup>Cre</sup>;Runx3<sup>fl/fl</sup>* mice, RUNX3 deletion was observed only from birth, with a modest 20% reduction in RUNX3 expression at P0 (Fig. S1C,D), when the sensorimotor circuit is already established (Mears and Frank, 1997). We thus decided to generate *Adv<sup>Cre</sup>;Runx3<sup>fl/fl</sup>* mice, with Cre expressed under the control of the advillin gene (Zhou et al., 2010). Compared with the *PV<sup>Cre</sup>* driver line, the *Adv<sup>Cre</sup>* line induced recombination in all PSNs between E13.5 and E15.5, and a cross with *Runx3<sup>fl/fl</sup>* mice resulted in a near-complete absence of RUNX3 expression in DRG at E15.5 (Fig. 1B,C; Fig. S1E,F). Together, this confirms the use of *Adv<sup>Cre</sup>* to delete *Runx3* in PSNs just after their peripheral innervation, and at the time they are growing their axons centrally within the spinal cord to reach their target. In P0 *Adv<sup>Cre</sup>;Runx3<sup>fl/fl</sup>* mice, focusing our analysis on the brachial segments, the number of NF200<sup>+</sup> DRG neurons (NF200 labels all myelinated DRG neurons, including PSNs) (Usoskin et al., 2015), and that of RUNX1<sup>+</sup> DRG neurons, which represent a large proportion of the unmyelinated DRG neurons at birth (Lallemand and Ernfors, 2012; Gascon et al., 2010), were similar between mutants and control animals (Fig. 1D,E). Hence, the absence of RUNX3 after peripheral innervation does not affect neuronal survival at this stage. However, the expression of factors necessary for the proper development and function of PSNs, such as ER81 and TRKC, and of the marker PV were all downregulated in P0 *Adv<sup>Cre</sup>;Runx3<sup>fl/fl</sup>* DRG compared with *Runx3<sup>fl/fl</sup>* control mice (Fig. 1F-K). These results show that

expression of RUNX3 in late embryonic DRG neurons is necessary for maintaining cell identity of a subgroup of PSNs.

### RUNX3 regulates development of central projection of subgroups of PSNs

From E13.5 to E15.5, the central afferents of PSNs enter the dorsal part of the spinal cord and send axons toward the MNs, which they make contact with at ~E17.5 (Mears and Frank, 1997; Ozaki and Snider, 1997). The development of the proprioceptive axonal inputs in the spinal cord has been suggested to involve NT3 and the transcription factors ER81 and RUNX3, with the deletion of any of them leading to severe deficits in the central projection pattern of PSNs (Patel et al., 2003; Nakamura et al., 2008; Arber et al., 2000). Yet those three factors are necessary for the peripheral outgrowth of PSNs and for their survival (De Nooij et al., 2013; Lallemand et al., 2012; Patel et al., 2003), preventing any study of their direct function on sensorimotor connectivity using null mutant mice. Here, we have explored whether the loss of RUNX3 in PSNs from E15 could perturb the formation of the sensorimotor connections between PSNs and MNs of the ventral spinal cord (phase iii, Fig. 1A). To this end, we analyzed the distribution of VGLUT1<sup>+</sup> sensory bouton contacts with CHAT<sup>+</sup> (choline acetyl transferase) MNs at P0, which reflects direct excitatory inputs from Ia PSNs on MNs (Oliveira et al., 2003; Pecho-Vrieseling et al., 2009). The innervation territory of PSNs in the spinal cord was divided into three distinct compartments (Fig. 2A): a ventro-medial (M) and a ventro-lateral (L) region corresponding presumably to axial and hypaxial muscle-derived Ia projections, respectively; and an intermediate zone (IZ), where most GTOs and type II MS afferents project (De Nooij et al., 2013; Lallemand and Ernfors, 2012). In *Adv<sup>Cre</sup>;Runx3<sup>fl/fl</sup>* and *Runx3<sup>fl/fl</sup>* mutant pups, VGLUT1 expression levels appeared unchanged in PSNs cell bodies (Fig. 2B,C). We observed a general decrease in the density of VGLUT1 labeling in these three regions at all brachial levels in the mutant (here shown for C5 and C8, Fig. 2D-F), with the largest defect found in the lateral and medial part of the ventral horn, which corresponds to the innervation territory of MNs (scheme in Fig. 1A; Fig. 2F). However, no defect was visible at the thoracic level (Fig. S2A). These results contrasted with the phenotype observed in the *Runx3<sup>-/-</sup>*; *Bax<sup>-/-</sup>* embryos in which, in spite of the absence of sensory neuron cell death due to null mutation of the pro-apototic gene *Bax* (Deckwerth et al., 1996), their central afferents did not grow beyond the dorsal aspect of the spinal cord, as revealed by immunostaining for peripherin at E15.5 and VGLUT1 at E18.5 (Fig. 2G,H). Taken together, our results indicate that PSNs at limb levels require RUNX3 expression to project correctly into the ventral spinal cord, independently of its role on the outgrowth dynamic of their peripheral projections. Nevertheless, some proprioceptive afferents still terminate both at the intermediate spinal cord and in the MN area in the absence of RUNX3. We thus hypothesized that RUNX3 might regulate the central connectivity of selective PSNs.

### Development of muscle-selective sensorimotor connections through RUNX3 activities

PSNs innervating specific skeletal muscles project centrally to and make contact with defined pools of MNs innervating the same muscles (Mears and Frank, 1997; Eccles et al., 1957; Dasen, 2009; Surmeli et al., 2011). This connection rule forms the structural basis of the muscle stretch reflex circuit. To explore whether the absence of RUNX3 affects preferentially the innervation of specific MN pools, we used Cholera toxin B subunit (CTB) injection into two major antagonistic groups of forelimb muscles: the biceps brachii
<figure><img src="image_3.png" /><figcaption>**Fig. 1.** Loss of cell identity in subgroups of PSNs following conditional targeting of RUNX3 after peripheral innervation. (A) Scheme representing the successive developmental steps of PSNs, which contribute to sensorimotor circuit. Early specification of PSNs (i) is followed by peripheral axonal growth and muscle targeting (circa. E14) (ii). After peripheral innervation, central afferents of PSNs project to the intermediate and then ventral regions of the spinal cord to contact interneurons and motor neurons (iii). (B) Ablation of *Runx3* from sensory neurons using *AdvCre* mice. At E13.5, RUNX3 expression is detectable in TRKC⁺ neurons with tdTomato (RFP) starting to be expressed in few neurons, while at E15.5, the recombination is fully efficient, all neurons expressing tdTomato and RUNX3 are strongly reduced in number. Scale bar: 50 µm. (C) Quantification of B, showing the recombination efficiency in TRKC⁺/RUNX3⁺ in *AdvCre;Runx3fl/fl* mice (n=3). ***P≤0.001; Student's t-test. Data are mean±s.e.m. (D) Immunostaining for NF200 and RUNX1 on DRG sections from P0 *AdvCre;Runx3fl/fl* and *Runx3fl/fl* animals identifies all myelinated sensory neurons (mechanoreceptive and proprioceptive neurons) and a large majority of nociceptive neurons (Lallemend and Ernfors, 2012; Gascon et al., 2010). Scale bar: 50 µm. (E) Quantification of D reveals absence of cell death in DRG neurons in the conditional *Runx3* mutants at P0. P>0.05, Student's t-test. Data are mean±s.e.m. (n=3). (F-K) Immunostaining for PSNs markers (F,H,J) and their quantification (G,I,K) in *AdvCre;Runx3fl/fl* and *Runx3fl/fl* P0 animals (n=3 per genotype). Scale bars: 100 µm in F; 50 µm in H,J. n=3 per genotype; *P*≤0.05; Student's t-test. (G) Data are mean±s.e.m.</figcaption></figure>

(flexor muscles) and the triceps brachii (extensor muscles, Fig. 3A). Biceps and triceps CHAT⁺ MNs (MNs<sup>Biceps</sup> and MNs<sup>Triceps</sup>, respectively) were retrogradely traced by injecting CTB in their respective muscles at postnatal day (P) 1, and their proprioceptive afferent terminals were examined at P2 by quantifying their number of contacts with VGLUT1⁺ boutons (Pecho-Vrieseling et al., 2009; Surmeli et al., 2011). The majority of MNs<sup>Biceps</sup> were found

essentially in C5-6 spinal segments, while MNs<sup>Triceps</sup> were found in C7-T1 spinal segments, as previously shown (Vrieseling and Arber, 2006). Gamma MNs, identified by their small size, were not considered for analysis. Quantification of the number of VGLUT1⁺ boutons making synapses with MNs<sup>Biceps</sup> revealed a 40% decrease in *AdvCre;Runx3fl/fl* mice compared with *Runx3fl/fl* control littermates (Fig. 3B,C). Strikingly, the majority of VGLUT1⁺
<figure><img src="image_4.png" /><figcaption>**Fig. 2.** Central afferentation deficit of PSNs after conditional targeting of RUNX3. (A) Representative regions of the spinal cord analyzed. To reveal central afferent terminations of PSNs, we used VGLUT1 immunostaining on cross-sections of spinal cord. The pattern of VGLUT1 reactivity was analyzed in three reference regions: the intermediate zone (IZ), the ventromedial (M) and the ventrolateral (L) regions. Scale bar: 100 µm. (B) VGLUT1 expression in *AdvCre; Runx3fl/fl* and *Runx3fl/fl* DRG sections at P0. Scale bar: 50 µm. (C) Quantification of the number of VGLUT1⁺ neurons per DRG section (left panel) and VGLUT1 intensity per cell (right panel) (from data in B) reveals absence of change in VGLUT1 expression in DRG from *AdvCre;Runx3fl/fl* mice (n=3). (D,E) Central innervation of PSNs in *AdvCre;Runx3fl/fl* and *Runx3fl/fl* mice at C5 (D) and C8 (E), as revealed by VGLUT1 immunostaining. Scale bar: 100 µm. (F) Quantification of the density of VGLUT1 staining in D and E in regions defined in A on one side of the spinal cord reveals deficits in central ingrowth of PSN afferents in conditional *Runx3* mutant mice (n=4 per genotype). A greater difference was observed in the lateral and medial (L and M) regions of the ventral spinal cord, which corresponds to the innervation of the MN pools (CHAT⁺ in A). **P**≤0.01, ***P**≤0.001; Student's *t*-test. Data are mean±s.e.m. (n=3). (G) Immunostaining for peripherin (PERI) on spinal cord sections shows complete absence of central PSN afferents in *Runx3⁻¹⁻;Bax⁻¹⁻* compared with *AdvCre;Runx3fl/fl* mice (see red arrowheads). Scale bar: 100 µm. (H) Similar to G, at P0 VGLUT1 immunostaining on spinal cord sections confirms the absence of central axon growth of PSNs in *Runx3⁻¹⁻;Bax⁻¹⁻* mice (see red arrowheads), a phenotype that differs from *AdvCre;Runx3fl/fl* mice (see D,E). Scale bar: 100 µm.</figcaption></figure>

boutons synapsing onto MNs<sup>Triceps</sup> were lost in the mutant (Fig. 3D,E), suggesting that RUNX3 expression in late embryonic stage has a selective role in the formation of sensorimotor connections between PSNs and defined pools of MNs.

## Selective loss of the peripheral innervation of PSNs in the absence of RUNX3

We next wondered whether the preferential sensorimotor connection deficits observed in *AdvCre;Runx3fl/fl* mice are reflected in the peripheral targeting of limb muscles by PSNs. MSs form after interaction between the sensory nerve endings of the PSNs and the targeted muscle fibers at ~E14 (Tourtellotte et al., 2001). Here, we quantified at P0 the VGLUT1⁺ sensory endings associated with MSs as readout of muscle innervation by PSNs (Fig. 4A) (Poliak et al., 2016). The GTO-associated nerve endings of the PSNs were not considered in the analysis. We focused on two pairs of antagonistic muscles spanning the proximo-distal axis of the forelimb: the extensors triceps and extensor carpi radialis (ECR); and the flexors biceps and flexor carpi radialis (FCR). Strikingly,

while the biceps and FCR muscles exhibited, respectively, a 30% decrease and normal incidence of VGLUT1⁺ PSN nerve endings in *AdvCre;Runx3fl/fl* mice, the triceps and ECR muscles displayed a large reduction in PSN nerve endings (Fig. 4B; Fig. S3A,B). In addition, the few remaining MSs observed in triceps and ECR muscles of *AdvCre;Runx3fl/fl* mutants were systematically of smaller size (Fig. 4A). These analyses at the forelimb level revealed a specific role for RUNX3 in PSN specification that affects MS differentiation in some defined muscle groups. They also show that the loss of intraspinal PSN axons targeting specific pools of MNs correlates well with the lack of sensory nerve endings in the limb muscles they innervate.

The observation that this innervation deficit is not associated with sensory neuron cell death in the *AdvCre;Runx3fl/fl* mutants at P0 raises the issue of what happens to those particular PSNs that have lost contact with their peripheral target. The peripheral projections of brachial PSNs make dorso-ventral choice within the forelimb before E12.5 (Kania et al., 2000; Huettl et al., 2011), and reach their peripheral target at E14 (Tourtellotte et al., 2001). By deleting
<figure><img src="image_5.png" /><figcaption>**Fig. 3.** Muscle-selective differential penetrance of central deficits of PSN connectivity in RUNX3 conditional mutants. (A) Experimental scheme describing muscle injection of the antagonistic muscles biceps and triceps, and retrograde labeling of specific MN pools at spinal segments C5-C6 for the biceps or segments C7-T1 for the triceps. (B) Immunostaining of spinal cord cross-sections representing MNs (CHAT+, in blue) traced by the CTB (red) from initial injection in the biceps. PSN synaptic contacts with the MNs are visualized by the VGLUT1+ synaptic bouton (green). Scale bars: 10 µm. (C) Quantification of 38 MNs (from B) in *Runx3fl/fl* and 39 MNs in *AdvCre;Runx3fl/fl*. ***P≤0.001; Student's t-test. Data are mean±s.e.m. (n=3). (D) Immunostaining of spinal cord cross-sections representing MNs (CHAT+, in blue) traced by the CTB (red) from initial injection in the triceps. PSN synaptic contacts with the MNs are visualized by the VGLUT1+ synaptic bouton (green). Scale bars: 10 µm. (E) Quantification of a total of 44 MNs (from D) in *Runx3fl/fl* and 46 MNs in *AdvCre;Runx3fl/fl*. ***P≤0.001; Student's t-test. Data are mean±s.e.m. (n=3).</figcaption></figure>

RUNX3 after this period using the *AdvCre* driver, it is thus unlikely that the PSNs losing their synaptic connection and identity will retract their projections to target another opposite muscle or different sensory organ such as the skin. Indeed, both the intrinsic signaling and environmental cues might not enable these cells to adopt a different sensory cell fate or trajectory. However, although future studies will be necessary to resolve this issue, the observed loss of identity features could eventually lead to cell death at later stage, as shown after the loss of many terminal selector genes that are needed to maintain a differentiated state (Deneris and Hobert, 2014).

## Role of muscle target and NT3 signaling in sustained RUNX3 expression

We next asked whether the phenotypic defects observed in PSNs in our *AdvCre;Runx3fl/fl* mutants could reflect a role for muscle-derived signals in PSN afferent specification. To assess this, we analyzed RUNX3 expression in *Lbx1* null mice, in which all extensor muscles are lacking at forelimb level although a large number of flexor muscles still develop (Gross et al., 2000; Brohmann et al., 2000). However, this phenotype is not accompanied by a loss in PSN number or by a defect in the peripheral outgrowth of PSN afferents within the limbs (Poliak et al., 2016). Examining brachial DRG neurons from E18.5 *Lbx1<sup>-/-</sup>* mice, we observed a 40% reduction in the number of neurons expressing RUNX3 (Fig. 5A,B).

Interestingly, the remaining RUNX3+ cells in *Lbx1* mutants exhibited reduced levels of expression similar to low RUNX3-expressing cells in wild-type animals [*Lbx1<sup>+/+</sup>*; 1±0.077; *Lbx1<sup>-/-</sup>*; 0.401±0.069 (expressed as mean±s.d.); *P*<0.01, normalized to *Lbx1<sup>+/+</sup>*], indicating that RUNX3 levels in PSNs might correlate with the specific target muscle groups they innervate. To test this, we analyzed backfilled PSNs following retrograde tracing using rhodamine dextran (Rh.dex.) injection in triceps and biceps from E16.5 wild-type embryos (Fig. 5C). Biceps- and triceps-innervating PSNs were primarily found in C5 and C8 segments, and expressed (on average) low and high levels of RUNX3, respectively, when compared with the whole PSN population (Fig. 5C-E; Fig. S4A). Together, these data suggest that RUNX3 expression and its level in PSNs depend on limb muscles at late embryonic stages. To support this, we analyzed *Bax<sup>-/-</sup>* mice devoid of NT3 signaling and in which PSNs completely fail to project into limb muscles but survive even in reduced neurotrophic factor conditions (Patel et al., 2003). RUNX3 expression was dramatically decreased in E15.5 *TrkC<sup>-/-</sup>*; *Bax<sup>-/-</sup>* DRG and virtually absent at birth (P0; Fig. 5F,G; Fig. S4B). Its expression in E13.5 *TrkC<sup>-/-</sup>*;*Bax<sup>-/-</sup>* DRG, i.e. prior to target innervation, however, was unchanged (Fig. S4B), indicating a late function of muscle-derived signals, following sensory afferent-muscle fiber contact, in the regulation of RUNX3 expression in PSNs. NT3 itself was shown to be a plausible candidate for controlling later aspects of PSNs specification (Wang et al., 2007;

<figure><img src="image_6.png" /><figcaption>**Fig. 4.** Muscle-specific MSs deficits in conditional *Runx3* mutants at birth. (A) Immunostaining for VGLUT1 (for MSs) and myosin on cross-sections from biceps and triceps. (B) Quantification reveals a muscle-selective MS deficiency in *AdvCre;Runx3fl/fl* with a significant decrease in the numbers of MSs in triceps (n=3 animals). **P**≤0.01, ****P≤0.0001; Student's t-test. Data are mean±s.e.m. Scale bar: 100 µm.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Fig. 5. RUNX3 dependence of PSN connectivity correlates with RUNX3 levels and muscle target NT3 levels.** (A) Immunostaining for ISL1 and RUNX3 on DRG sections from *Lbx1*<sup>-/-</sup> and *Lbx1*<sup>+/+</sup> mice. Scale bar: 50 µm. (B) Quantification of RUNX3<sup>+</sup> neurons (from A) in *Lbx1*<sup>-/-</sup> (*n*=2) and their control littermates (*n*=4) at brachial levels shows significant reduction in the full mutants compared with control animals. **P**≤0.01. Data are mean±s.e.m. (C) Experimental design. (Left) Rhodamine dextran (Rh. dex.) injection in specific muscle will retrogradely trace their innervating PSNs in the DRG from E16.5 wild-type animals. The dextran is injected in biceps in one of the forelimbs and in triceps in the contralateral limb. (Right) Distribution within DRG of the PSNs innervating biceps and triceps, showing biceps- and triceps-innervating PSNs located mostly in DRG C5 and C8, respectively (data from 5 animals for triceps and 11 animals for biceps). (D,E) Quantification of RUNX3 expression per cell in retrogradely traced PSNs (as in C) versus all RUNX3<sup>+</sup> PSNs after injection of Rh. dex. in triceps (analysis at C8 level, D) or in biceps (analysis at C5 level, E). **P**≤0.01; Student's *t*-test. Data are mean±s.e.m. (single neurons analyzed from five embryos). (F) RUNX3 expression is largely reduced in *TrkC*<sup>-/-</sup>:*Bax*<sup>-/-</sup> P0 mice: a mouse model of peripheral outgrowth deficits. Scale bar: 100 µm. (G) Quantification of F reveals an almost complete absence of RUNX3 expression in *TrkC*<sup>-/-</sup>:*Bax*<sup>-/-</sup> P0 mice compared with their control littermates (*n*=2). **P**≤0.01; Student's *t*-test. Data are mean±s.e.m. (H) X-Gal reaction on E15.5 in *Ntf3LacZ/+* forelimb embryos show heterogeneous muscle-specific expression of NT3. There is a large difference in NT3 levels between biceps and triceps. (I) Quantification of NT3 mRNA (*Ntf3*) in the biceps (Bic) and triceps (Tri) shows a twofold increase in triceps compared with biceps in control animals (Ctr). *P*=0.052; Student's *t*-test (*n*=2 samples from four animals). Data are mean±s.e.m. (J) *Bax*<sup>-/-</sup> mice DRG explants (E15.5) in culture with or without NT3 for 48 h *in vitro* (HIV) reveal a decreased expression of RUNX3 in the absence of NT3. Scale bar: 100 µm. (K) Quantification of J shows a significant increase of RUNX3 intensity per cell, while the number of positive neurons remains unchanged. **P**≤0.01; Student's *t*-test. Data are mean±s.e.m. (*n*=3). (L) Immunostaining for ISL1, RET and RUNX3 on spinal cord (SC) sections from *Hb9Cre;Islet2DTA* mice shows complete absence of MNs at E16.5. Scale bar: 50 µm. (M) Immunostaining for TRKC, RUNX3 and ISL1 on DRG sections from *Hb9Cre;Islet2DTA* and *Islet2DTA* E16.5 mice shows no deficits in RUNX3 and TRKC expression in the absence of MNs, confirmed by the quantification of the number of PSNs in *Hb9Cre;Islet2DTA* and *Islet2DTA* (*n*=2, right panel). Scale bar: 20 µm.</figcaption></figure>

De Nooij et al., 2013). Notably, elevated levels of NT3 in muscles late during embryonic development has been shown to disrupt the selective pattern of synaptic connections between sensory afferents and MNs (Wang et al., 2007). In addition, various levels of NT3 in hindlimb muscles have been correlated with subclass-specific sensitivity of PSNs to cell death following loss of ER81 (De Nooij et al., 2013). This heterogeneity in NT3 levels among muscles was also observed here in forelimb muscles using β-gal activity in

whole-mounted *Ntf3LacZ/+* E15.5 embryos and confirmed using qPCR (Fig. 5H,I). To directly test the role of NT3 on RUNX3 expression at late embryonic stage, we cultured E15.5 DRG explants from *Bax*<sup>-/-</sup> embryos for 2 days in the presence or absence of NT3 (Fig. 5J). Sensory neurons from *Bax*<sup>-/-</sup> DRG could survive several days *in vitro* without neurotrophin signaling (Lentz et al., 1999; Zhong et al., 2007). We found higher levels of RUNX3 expression in the presence of NT3 compared with the control condition
(Fig. 5K). Supporting these data, addition of NT3 to DRG explants from HH27 chicken embryos resulted in marked upregulation of *Runx3* mRNA levels compared with untreated DRG explants (Fig. S4C,D). Thus, after peripheral innervation, distinct levels of *RUNX3* expression in PSN subgroups are controlled by muscle-derived signals and involve NT3 signaling.

We also considered whether MNs, through direct axon-axon contact or indirect signaling, could act on RUNX3 expression in PSNs. To test this, we used *Hb9cre;Isl2DTA* mice, in which diphtheria toxin (DTA) is selectively expressed in HB9⁺/ISL2⁺ MNs, ablating MNs as they exit the cell cycle, before peripheral innervation (Yang et al., 2001) (Fig. 5L). In the absence of MNs, PSNs density in DRG is unchanged, and their nerve endings can still be found in association with MSs (Poliak et al., 2016). In *Hb9cre; Isl2DTA* mice examined at E16.5, RUNX3 expression in PSNs was found unchanged compared with *Isl2DTA* mice (Fig. 5M). These data thus argue against the idea that MNs could provide signals necessary for regulating or maintaining RUNX3 in PSNs at late embryonic stages.

## DISCUSSION

Muscle sensory feedback is essential for controlling coordination of locomotor behavior. It is conveyed by different subgroups of PSNs that form specific connections with second order neurons and MNs in the spinal cord. Despite the importance of this feedback in motor behavior, our knowledge of the molecular programs that control the assembly of the sensorimotor circuits remains limited. This is essentially due to the fact that the major genetic determinants of PSNs identity, which continue to be expressed in the mature neurons, are also necessary for their early peripheral outgrowth and survival, limiting the study of their direct function in stage-specific maintenance of PSNs identity and connectivity. Here, we have used conditional Cre/loxP gene-targeting approach to delete the major regulatory factor *Runx3* in PSNs after the period of natural cell death and peripheral innervation. We find that RUNX3 is continuously required throughout development to maintain specific identity features of particular subgroups of PSNs innervating the limb, independent of cell death. We also show that the maintained expression of RUNX3 in PSNs is regulated by muscle target-derived signals and NT3 signaling which is differentially required in subgroups of PSNs to setup specific sensorimotor connections.

## Dedicated maintenance factors for PSN subgroups

The generation of proprioceptive afferents during embryonic development follows a systematic sequence of events, from neurogenesis, early cell fate specification and target innervation to finer aspects of PSN subgroup-specific identity and the formation of the sensorimotor circuits (Lallemand and Ernfors, 2012). Hence, soon after neurogenesis, TRKC⁺ PSNs acquire RUNX3 expression, which is necessary to consolidate a proprioceptive fate and direct the extent of peripheral outgrowth of PSNs projections (Lallemand et al., 2012; Levanon et al., 2002; Kramer et al., 2006). As their axons grow peripherally, PSNs, through target-derived NT3 signaling, will initiate expression of ER81, which later will define the survival competence of particular subgroups of PSNs in correlation with their target muscles (Patel et al., 2003; De Nooij et al., 2013). At the same time, different subtypes of PSNs emerge on the basis of their preferential expression of recognition molecules defined by interactions with signals notably from the limb (Poliak et al., 2016; Wu et al., 2019). In this context, it has been unclear whether early specification markers continue to be crucial during

late embryogenesis for preservation and maturation of the sensory neuron type identity and connectivity. Our data indicate that RUNX3 is not essential for the survival of PSNs following peripheral innervation. This suggests that the loss of PSNs observed in *Runx3*-null mice at E12.5 is likely a consequence of the severe reduction in TRKC expression and survival signaling within PSNs during the developmental cell death period (Lallemand et al., 2012; Levanon et al., 2002; Kramer et al., 2006). This is reminiscent of a recent study exploring the role of BRN3A in the diversification of sensory neurons of the cochlea, where conditional deletion of *Pou4f1* (encoding BRN3A) from E14.5 (while it is first expressed at E10.5) did not affect the survival of cochlear neurons as seen in the full mutant where TRKC is found to be downregulated (Huang et al., 2001; Sherrill et al., 2019). In contrast, classic key molecular features of PSN identity, such as ER81, TRKC and PV expression, showed a deficit in our conditional *Runx3* mutant. Interestingly, only a subgroup of PSNs, representing about half of its population at the forelimb level, showed a loss of identity marker expression. This implies the existence of two distinct subgroups of PSNs at limb levels based on their dependence on RUNX3 expression for the later aspect of their differentiation. In support of this, key molecular markers of PSNs are found unchanged in thoracic DRG in our conditional *Runx3* mutant mice. One possible mechanism that could participate in the development of these two distinct subgroups of PSNs is the graded level of RUNX3 activity itself. A previous study in chicken embryos using *in ovo* electroporation suggested that the status of RUNX3 expression could participate in the segregation of the different subtypes of PSNs (innervating GTOs or MSs), notably by regulating their central projection trajectory (Chen et al., 2006). Our data in mice indicate instead that the levels of RUNX3 expression and its role in PSN diversification might correlate with the identity of the target muscles being innervated, with mean RUNX3 levels decreasing in remaining PSNs projecting to forelimb flexor muscles in *Lbx1* null mice, which are devoid of extensor muscles (Gross et al., 2000; Brohmann et al., 2000). These contrasting data could indicate a species difference. In support, siRNA-based extinction of RUNX3 expression in embryonic chicken PSNs did not lead to loss of TRKC (Chen et al., 2006), which would be expected from the data found in mice (Lallemand et al., 2012; Kramer et al., 2006). Our data remind us, however, of the differential activity of the NT3-ER81 signaling pathway in PSNs and its role in delineating two classes of PSNs based on their dependence on ER81 for their survival (De Nooij et al., 2013). Although results from our work would need further investigation of a direct link between RUNX3 levels and PSN diversification, graded levels of transcription factor activity may play a crucial role in the development of neuronal identity. In this context, the regional differences in the expression at the target of signaling cues might serve at the nerve endings as key factors eliciting various strengths of retrograde signaling operating in the neurons to drive their identity. A very good example is the development of the trigeminal ganglion, in which the coincident expression of both BDNF and BMP4 in discrete regions of the peripheral target field plays an important role in controlling the development of specific neuronal pools within the ganglion (Ji and Jaffrey, 2012). Distinct signals derived from limb mesenchyme have also been hypothesized to drive the subclass identity of PSNs. In a similar way, we show here that the presence of muscles is necessary for maintaining the expression of RUNX3 in PSNs. Moreover, NT3, which is expressed at various levels in developing skeletal muscles (Fariñas et al., 1996; De Nooij et al., 2013), can directly regulate RUNX3 expression *ex vivo*. Hence, it is conceivable that various extents of NT3-TRKC
signaling activities between PSNs subgroups innervating distinct muscles might regulate their later cell identity aspects through graded RUNX3 expression. Following peripheral innervation, such differential activities among distinct PSNs subgroups could eventually serve to link the targeted muscles with the central neurons in the spinal cord the PSNs must connect.

## Selectivity of sensorimotor connectivity

Certainly the most studied circuit in sensorimotor network connecting PSNs and spinal cord neurons is the reflex arc in which group Ia PSN afferents make specifically strong connections with MNs supplying the same muscle and weaker connections with motor neurons supplying synergistic muscles, while no connection is established with MNs of antagonistic muscles (Eccles et al., 1957; Frank and Mendelson, 1990). During the formation of this circuit, an elegant study has shown that the topography of the PSN endings along the proximodistal axis of the limb could predict the dorsoventral location and identity of their target MNs, suggesting an area-specific targeting of PSNs within the ventral spinal cord (Surmeli et al., 2011). Others emphasized the importance of a direct MN-dependent mechanism that involves cell-to-cell recognition with PSNs and the use of Sema3E-plexin D1 signaling for the formation of sensorimotor connectivity patterns (Fukuhara et al., 2013; Peko-Vrieseling et al., 2009). This implies a stepwise targeting approach making use of different signaling and strategies that PSN must use to eventually synapse to the appropriate MNs. Our data here on the selective function of RUNX3 in the spinal cord ingrowth of particular subgroups of PSNs indicate another level of regulation by which only a subset of neurons critically depends on the activity of a specific transcriptional program to project within the spinal cord and eventually develop a proper connectivity. Interestingly, this RUNX3 dependence of PSN-MN connections exhibits a preference for neurons innervating forelimb extensors, when compared with flexors and axial muscles at limb level, as well as thoracic muscles. In adult animals, such a deficit of muscle spindles in limb extensors, in acute situation, would most likely cause a substantial decrease in muscle tone of extensors, resulting in a change in the basic body configuration (a configuration with more flexed limbs as compared with control) and a decrease in efficacy of feedback postural corrections. In our mutant mice, due to adaptive plastic changes that could take place during early life, these deficits could be compensated for, e.g. by an increase in the activity of neurons of the vestibulospinal tract leading to an increase in extensor muscle tone and/or by a substitution of proprioceptive information (which in control subjects plays a crucial role in the generation of postural corrections) by visual and vestibular signals. However, the most likely specific aspects of postural corrections might still differ from those observed in control animals.

The RUNX3-dependent bias to neurons innervating forelimb extensors, compared with flexors muscles, evokes the ER81 dependence of PSNs for their survival, where PSNs innervating thoraco-hypaxial and axial muscles exhibit an almost complete dependence on ER81, whereas those innervating hindlimb muscles exhibit a mosaic, muscle-by-muscle, sensitivity to ER81 inactivation (De Nooij et al., 2013). Strikingly, the levels of NT3 expression in individual muscles strongly correlated with this dependence on ER81 of the survival of particular PSNs subgroups. Although our study only considers selected sensorimotor connections, it suggests a role for peripheral gradation of NT3 in setting up the sensitivity of PSNs for RUNX3-dependent sensorimotor connectivity, as described above. Even if further programs are certainly required to create the fine topography of

PSNs to MNs connectivity, this earlier step already distinguishes populations of neurons as they need to enter the spinal cord, as shown in chick and mouse for the selectivity of the sensorimotor circuits by regionally restricted limb signals (Wenner and Frank, 1995; Poliak et al., 2016). This could serve to coordinate spatially and temporally separated developmental events, as observed during the development of olfactory receptor neurons and their selective connectivity with their glomerular target (Chou et al., 2010). The implication of these results is also that, within proprioceptors, higher levels of RUNX3 might engage a distinct set of transcriptional targets that are necessary for their proper maturation. The lower RUNX3-expressing PSNs would then require the activation of another set of transcriptional regulators. Such a mechanism, which most likely involves differences in target-derived signals, including neurotrophins, would thus participate in a spatial patterning of transcription factor activities among sensory neurons.

## MATERIAL AND METHODS

### Animals

Wild-type C57BL6 mice were used unless specified otherwise. *Runx3*<sup>−/−</sup>, *Advillin*<sup>Cre</sup> (*Adv*<sup>Cre</sup>) *Lbx1*<sup>−/−</sup>, *Hb9*<sup>Cre</sup> and *Isl2*<sup>DTA</sup> mouse strains have been described elsewhere (Levanon et al., 2002; Yang et al., 2001; Gross et al., 2000; Zhou et al., 2010). *Runx3*<sup>Im1Itan</sup> (*Runx3*<sup>fl</sup>), *Bax*<sup>−/−</sup>, *TrkC*<sup>−/−</sup> and *Ai14* mice were purchased from Jackson Laboratories, and *NT3*<sup>−/−</sup> from MMRRC. Both male and female animals were included in this study, except for the *Adv*<sup>Cre</sup>; *Runx3*<sup>fl/fl</sup>, where only males were analyzed because of weak expression of Cre in oocytes and possible leakage or generation of 'deleter' mice. Animals were group-housed, with food and water available *ad libitum*, under 12 h light-dark cycle conditions. Fertile white Leghorn eggs were incubated at 38°C and embryos were staged according to Hamburger-Hamilton (HH) tables. All animal work was performed in accordance with the national guidelines and approved by the local ethics committee of Stockholm, Stockholms Norra djurförsöksetiska nämd.

### In vitro cultures of whole DRG

Whole DRG were cultured on coated plates (5% matrigel in ice-cold PBS) with N2 medium (DMEM-F12/glutamax medium with N2 supplement; Gibco) supplemented with pen/strep, gentamicin, the pan-caspase inhibitor Q-VD-Oph (2 µM, Sigma) and with NT3 (Peprotech), retinoic acid, bFGF and IGF when specified, as previously described (Hadjab et al., 2013). DRG analyzed for mRNA expression were cultured directly in eppendorf tubes.

### qPCR

Tissue was freshly dissected and directly placed into lysis buffer. RNA was extracted using the Qiagen RNAEasy Mini Kit according to the manufacturer's instructions, including DNase I to degrade potential genomic DNA contamination. The amount of RNA was quantified using the Qubit RNA BR Assay Kit, and 50 ng of each biological replicate used for reverse transcription using Biorad iScript in a 20 µl reaction according to the manufacturer's instructions. Produced cDNAs were used to analyze transcript levels using real-time PCR in a Biorad CFX96 Real-time System and using Biorad iTaq Universal SYBR Green Supermix in 20 µl reactions for 40 cycles (95°C for 10 s, 60°C for 20 s, 72°C for 30 s followed by fluorescence measurement at 75°C to reduce the influence of primer dimer on quantification). Primers for *Nft3* and *Runx3* were designed using Primer3 Blast (NCBI). Primers for *Nft3*, *Ntrk3* (TRKC), *Isl1* and *Runx3* were designed using Primer3 Blast (NCBI). *Gapdh* was used as a reference gene. A full list of primers used is given below. Negative controls did not amplify. Data were analyzed using Bio-Rad CFX manager. *Gapdh*, *AACTCCCACT-TCTCCACCTTC* and *GATAGGGCCTCTCTTGCTCAG*; *Nft3*, *ATAAAA-TTCGTGTGCTTGCCTTCC* and *GAGAGCCCAATCACAAAACAAGG*; *Runx3*, *ACCAAGTGGCGAGATTTAACGA* and *ACGGTGACTTTAATGGCTCGG*; *Isl1*, *AAAAGAAGCATTATGATGAAGCAA* and *CATGTC-TCTCCGGACTAGCAG*; *Ntrk3*, *TGATCCTCGTGGATGGGCAG* and *CTTCACCAGCAGGTTGGCTCC*.
## Immunostaining

Mouse embryos were collected and fixed for 1 to 6 h at +4°C (4% PFA in PBS) depending on the stage, washed in PBS three times (30 min each), equilibrated in 20% and 30% sucrose in PBS, embedded in OCT (Tissue-Tek) and cryosectioned at 14 µm. Sections were incubated for one or two nights with primary antibodies diluted in blocking solution (2% donkey serum, 0.0125% NaN₃, 0.5% Triton X-100 in PBS). Primary antibodies used were: rabbit anti-RUNX3 (a gift from T. M. Jessell and D. Levanon, Columbia University and the Weizmann Institute of Science), goat anti-TRKA (R&D Systems, AF1056, 1/500), goat anti-TRKB (R&D Systems, AF1494, 1/500), goat anti-TRKC (R&D Systems, AF1404, 1/500), goat anti-RET (R&D Systems, AF482, 1/500), mouse anti-ISL1 (Developmental Studies Hybridoma Bank, 39.4D5, 1/200), rabbit anti-VGLUT1 (SYSY 135303, 1/1000), chicken anti-RFP (Rockland, 600-401-379, 1/500), goat anti-ChAT (Millipore, AB144p, 1/1000), mouse anti-NF200 (Sigma, N0142), rabbit anti-peripherin (Millipore, AB1530, 1/500), rabbit anti-PEA3 (a gift from S. Arber, Friedrich Miescher Institute, Basel, Switzerland) and mouse anti-myosin (DSHB, F59-s, 1/50). After washing with PBS, Alexa Fluor secondary antibodies (Life Technology; 1:500 in blocking solution) were applied overnight (at +4°C). Samples were then washed in PBS and mounted in DAKO fluorescent mounting medium. Staining was imaged by confocal microscopy (Zeiss LSM700 or LSM800) using identical settings between control and experimental groups. Optical pinholes were 2 µm under 20× magnification unless specified otherwise.

## Quantification of neurons

For cell type counting quantifications, ImageJ software was used. Only neurons with a visible nucleus were considered for analysis. Quantification of molecular markers in the DRG was carried out on five DRG sections/animal, selected from the most equatorial region of each DRG and covering the segments C5-T1 (see figure legends for *n* values and genotypes).

## Quantification of muscle spindles

Limb muscle VGLUT1⁺ muscle spindles were counted on longitudinal sections of whole muscles.

## Quantification of PSN collateral density

Quantitative analysis of PSN axon terminals in the ventrolateral horn of the spinal cord was performed using ImageJ analysis software. For each segment level (C5 and C8) the total VGLUT1⁺ collateral surface area (expressed as squared microns) was measured within a confined intermediate zone, ventromedial area and ventrolateral area of the ventral horn (adjusted for each segmental level). The surface areas analyzed were as follows: C5 ventrolateral, 210,000 µm²; C5 ventromedial, 112,000 µm²; C5 intermediate zone, 37,400 µm²; C8 ventrolateral, 160,000 µm²; C8 ventromedial, 63,000 µm²; and C8 intermediate zone, 30,000 µm². For each genotype, measurements for all segmental levels were performed from at least two individual experiments and three z-stacks per segment.

## Quantification of sensory synaptic contacts with motor neurons

Quantification of VGLUT1⁺ sensory bouton contacts between PSNs and motor neurons somata was performed using 0.5 µm confocal Z-scans of 10 to 20 µm thick sections. Motor neuron surface area was determined using Image J. Synapses were counted manually and filtered for sizes larger than 1 µm².

## β-Galactosidase staining

LacZ expression was detected by staining using X-Gal (5-bromo-4-chloro-3-indolyl β-d-galactoside) for β-galactosidase activity as described previously (Lallemend et al., 2012).

## Retrograde labeling of motor and sensory neurons

Newborn mice were anesthetized on ice by hypothermia and muscles of interest were injected with 1% cholera toxin B subunit/Alexa555 (CTB555). The amount of CTB555 was enough to fill the whole muscle. The following day, spinal columns were dissected and fixed in 4% PFA for 6 h then processed for immunostaining and analysis. Sensory neuron tracing was

performed *ex vivo* by injection of fluorescently labeled dextrans (3000MW, Invitrogen) in limb muscles using tightly fitting glass capillaries. Following injection, the *ex vivo* preparation (brachial region with associated limbs) was incubated for 8 h in oxygenated ACSF before fixation.

## Statistics

Data were analyzed using GraphPad Prism 6 and are expressed as mean±s.e.m. The statistical test performed is reported in the figure legend. *t*-tests were two-sided (*P*≤0.05, ***P*≤0.01, ***P*≤0.001). No animals or data points were excluded from the analyses. No statistical methods were used to pre-determine sample size but our sample sizes are similar to those generally employed in the field.

## Acknowledgements

We thank Prof. Yoram Groner and Ditsa Levanon for the *Runx3⁻/⁻* mouse line. We thank the CLICK imaging Facility supported by the Knut and Alice Wallenberg Foundation.

## Competing interests

The authors declare no competing or financial interests.

## Author contributions

Conceptualization: S.H., F.L.; Methodology: Y.W., H.W., P.F., S.W., C.P., R.-E.H., A.H.B., Y.X.-F.; Validation: H.W., P.F., S.W.; Formal analysis: Y.W., H.W., P.Z., P.F., S.W., C.P., C.B., T.G.D.; Investigation: S.H., F.L.; Resources: P.Z., G.C., C.B., S.T., O.K., P.E., T.G.D.; Data curation: S.H., F.L.; Writing - original draft: S.H., F.L.; Writing - review & editing: Y.W., H.W., P.F., S.W., C.P., G.C., C.B., Y.X.-F.; Visualization: S.H., Y.W.; Supervision: S.H., F.L.; Project administration: S.H., F.L.; Funding acquisition: S.H., F.L.

## Funding

This work was supported by StratNeuro (2013-0237), the Vetenskapsrådet (2012-04708), a Karolinska Institute doctoral grant the Knut och Alice Wallenberg Stiftelse (Wallenberg Academy Fellow), the Hjärnfonden (FO2014-0048), the Karolinska Institutet (Faculty Funded Career Position) and the Ragnar Söderbergs stiftelse (M48/12) (Ragnar Söderberg Fellow in Medicine), and by a Ming Wai Lau research grant. F.L. is a Ragnar Söderberg fellow in Medicine, a Wallenberg Academy Fellow in Medicine and a Ming Wai Lau Center investigator. Deposited in PMC for immediate release.

## Supplementary information

Supplementary information available online at http://dev.biologists.org/lookup/doi/10.1242/dev.181750 supple mental

## References

* Arber, S., Ladle, D. R., Lin, J. H., Frank, E. and Jessell, T. M. (2000). ETS gene Er81 controls the formation of functional connections between group Ia sensory afferents and motor neurons. *Cell* **101**, 485-498. doi:10.1016/S0092-8674(00)80859-4

* Brohmann, H., Jagla, K. and Birchmeier, C. (2000). The role of Lbx1 in migration of muscle precursor cells. *Development* **127**, 437-445.

* Chen, A. I., de Nooij, J. C. and Jessell, T. M. (2006). Graded activity of transcription factor Runx3 specifies the laminar termination pattern of sensory axons in the developing spinal cord. *Neuron* **49**, 395-408. doi:10.1016/j.neuron.2005.12.028

* Chou, Y.-H., Zheng, X., Beachy, P. A. and Luo, L. (2010). Patterning axon targeting of olfactory receptor neurons by coupled hedgehog signaling at two distinct steps. *Cell* **142**, 954-966. doi:10.1016/j.cell.2010.08.015

* Dasen, J. S. (2009). Transcriptional networks in the early development of sensory-motor circuits. *Curr. Top. Dev. Biol.* **87**, 119-148. doi:10.1016/S0070-2153(09)01204-6

* Dasen, J. S., Tice, B. C., Brenner-Morton, S. and Jessell, T. M. (2005). A Hox regulatory network establishes motor neuron pool identity and target-muscle connectivity. *Cell* **123**, 477-491. doi:10.1016/j.cell.2005.09.009

* Deckwerth, T. L., Elliott, J. L., Knudson, C. M., Johnson, E. M., Jr, Snider, W. D. and Korsmeyer, S. J. (1996). BAX is required for neuronal death after trophic factor deprivation and during development. *Neuron* **17**, 401-411. doi:10.1016/S0896-6273(00)80173-7

* De Marco Garcia, N. V. and Jessell, T. M. (2008). Early motor neuron pool identity and muscle nerve trajectory defined by postmitotic restrictions in Nkx6.1 activity. *Neuron* **57**, 217-231. doi:10.1016/j.neuron.2007.11.033

* Deneris, E. S. and Hobert, O. (2014). Maintenance of postmitotic neuronal cell identity. *Nat. Neurosci.* **17**, 899-907. doi:10.1038/nn.3731

* De Nooij, J. C., Doobar, S. and Jessell, T. M. (2013). Etv1 inactivation reveals proprioceptor subclasses that reflect the level of NT3 expression in muscle targets. *Neuron* **77**, 1055-1068. doi:10.1016/j.neuron.2013.01.015
Eccles, J. C., Eccles, R. M. and Lundberg, A. (1957). The convergence of monosynaptic excitatory afferents on to many different species of alpha motoneurons. *J. Physiol.* **137**, 22-50. doi:10.1113/jphysiol.1957.sp005794

Fariñas, I., Yoshida, C. K., Backus, C. and Reichardt, L. F. (1996). Lack of neurotrophin-3 results in death of spinal sensory neurons and premature differentiation of their precursors. *Neuron* **17**, 1065-1078. doi:10.1016/S0896-6273(00)80240-8

Fariñas, I., Wilkinson, G. A., Backus, C., Reichardt, L. F. and Patapoutian, A. (1998). Characterization of neurotrophin and Trk receptor functions in developing sensory ganglia: direct NT-3 activation of TrkB neurons in vivo. *Neuron* **21**, 325-334. doi:10.1016/S0896-6273(00)80542-5

Frank, E. (1990). The formation of specific synaptic connections between muscle sensory and motor neurons in the absence of coordinated patterns of muscle activity. *J. Neurosci.* **10**, 2250-2260. doi:10.1523/JNEUROSCI.10-07-02250. 1990

Frank, E. and Mendelson, B. (1990). Specification of synaptic connections mediating the simple stretch reflex. *J. Exp. Biol.* **153**, 71-84.

Fukuhara, K., Imai, F., Ladle, D. R., Katayama, K., Leslie, J. R., Arber, S., Jessell, T. M. and Yoshida, Y. (2013). Specificity of monosynaptic sensory-motor connections imposed by repellent Sema3E-PlexinD1 signaling. *Cell Rep.* **5**, 748-758. doi:10.1016/j.celrep.2013.10.005

Gascon, E., Gaillard, S., Malapert, P., Liu, Y., Rodat-Despoix, L., Samokhvalov, I. M., Delmas, P., Helmbacher, F., Maina, F. and Moqrich, A. (2010). Hepatocyte growth factor-Met signaling is required for Runx1 extinction and peptidergic differentiation in primary nociceptive neurons. *J. Neurosci.* **30**, 12414-12423. doi:10.1523/JNEUROSCI.3135-10.2010

Gross, M. K., Moran-Rivard, L., Velasquez, T., Nakatsu, M. N., Jagla, K. and Goulding, M. (2000). Lbx1 is required for muscle precursor migration along a lateral pathway into the limb. *Development* **127**, 413-424.

Hadjab, S., Franck, M. C. M., Wang, Y., Sterzenbach, U., Sharma, A., Ernfors, P. and Lallemand, F. (2013). A local source of FGF initiates development of the unmyelinated lineage of sensory neurons. *J. Neurosci.* **33**, 17656-17666. doi:10.1523/JNEUROSCI.1090-13.2013

Honig, M. G. (1982). The development of sensory projection patterns in embryonic chick hind limb. *J. Physiol.* **330**, 175-202. doi:10.1113/jphysiol.1982.sp014336

Honig, M. G., Frase, P. A. and Camilli, S. J. (1998). The spatial relationships among cutaneous, muscle sensory and motoneuron axons during development of the chick hindlimb. *Development* **125**, 995-1004.

Huang, E. J., Liu, W., Fritzsch, B., Bianchi, L. M., Reichardt, L. F. and Xiang, M. (2001). Brn3a is a transcriptional regulator of soma size, target field innervation and axon pathfinding of inner ear sensory neurons. *Development* **128**, 2421-2432.

Huetti, R.-E., Soellner, H., Bianchi, E., Novitch, B. G. and Huber, A. B. (2011). Npn-1 contributes to axon-axon interactions that differentially control sensory and motor innervation of the limb. *PLoS Biol.* **9**, e1001020. doi:10.1371/journal.pbio.1001020

Ji, S.-J. and Jaffrey, S. R. (2012). Intra-axonal translation of SMAD1/5/8 mediates retrograde regulation of trigeminal ganglia subtype specification. *Neuron* **74**, 95-107. doi:10.1016/j.neuron.2012.02.022

Kania, A., Johnson, R. L. and Jessell, T. M. (2000). Coordinate roles for LIM homeobox genes in directing the dorsoventral trajectory of motor axons in the vertebrate limb. *Cell* **102**, 161-173. doi:10.1016/S0092-8674(00)00022-2

Kramer, I., Sigrist, M., De Nooij, J. C., Taniuchi, I., Jessell, T. M. and Arber, S. (2006). A role for Runx transcription factor signaling in dorsal root ganglion sensory neuron diversification. *Neuron* **49**, 379-393. doi:10.1016/j.neuron.2006.01.008

Lallemand, F. and Ernfors, P. (2012). Molecular interactions underlying the specification of sensory neurons. *Trends Neurosci.* **35**, 373-381. doi:10.1016/j.tins.2012.03.006

Lallemand, F., Sterzenbach, U., Hadjab-Lallemand, S., Aquino, J. B., Castelo-Branco, G., Sinha, I., Villaescusa, J. C., Levanon, D., Wang, Y., Franck, M. C. et al. (2012). Positional differences of axon growth rates between sensory neurons encoded by runx3. *EMBO J.* **31**, 3718-3729. doi:10.1038/emboj.2012.228

Lentz, S. I., Knudson, C. M., Korsmeyer, S. J. and Snider, W. D. (1999). Neurotrophins support the development of diverse sensory axon morphologies. *J. Neurosci.* **19**, 1038-1048. doi:10.1523/JNEUROSCI.19-03-01038.1999

Levanon, D., Bettoun, D., Harris-Cerruti, C., Woolf, E., Negreanu, V., Eilam, R., Bernstein, Y., Goldenberg, D., Xiao, C., Fliegauf, M. et al. (2002). The Runx3 transcription factor regulates development and survival of TrkC dorsal root ganglia neurons. *EMBO J.* **21**, 3454-3463. doi:10.1093/emboj/cdf370

Mears, S. C. and Frank, E. (1997). Formation of specific monosynaptic connections between muscle spindle afferents and motoneurons in the mouse. *J. Neurosci.* **17**, 3128-3135. doi:10.1523/JNEUROSCI.17-09-03128.1997

Mendelsohn, A. I., Simon, C. M., Abbott, L. F., Mentis, G. Z. and Jessell, T. M. (2015). Activity regulates the incidence of heteronymous sensory-motor connections. *Neuron* **87**, 111-123. doi:10.1016/j.neuron.2015.05.045

Mendelson, B. and Frank, E. (1991). Specific monosynaptic sensory-motor connections form in the absence of patterned neural activity and motoneuronal cell death. *J. Neurosci.* **11**, 1390-1403. doi:10.1523/JNEUROSCI.11-05-01390. 1991

Nakamura, S., Senzaki, K., Yoshikawa, M., Nishimura, M., Inoue, K., Ito, Y., Ozaki, S. and Shiga, T. (2008). Dynamic regulation of the expression of neurotrophin receptors by Runx3. *Development* **135**, 1703-1711. doi:10.1242/ dev.015248

Oliveira, A. L. R., Hydling, F., Olsson, E., Shi, T., Edwards, R. H., Fujiyama, F., Kaneko, T., Hökfelt, T., Cullheim, S. and Meister, B. (2003). Cellular localization of three vesicular glutamate transporter mRNAs and proteins in rat spinal cord and dorsal root ganglia. *Synapse* **50**, 117-129. doi:10.1002/syn.10249

Ozaki, S. and Snider, W. D. (1997). Initial trajectories of sensory axons toward laminar targets in the developing mouse spinal cord. *J. Comp. Neurol.* **380**, 215-229. doi:10.1002/(SICI)1096-9861(19970407)380:2<125::AID-CNE5>3.0.CO;2-6

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

Pecho-Vrieseing, E., Sigrist, M., Yoshida, Y., Jessell, T. M. and Arber, S. (2009). Specificity of sensory-motor connections encoded by Sema3e-Plnd1 recognition. *Nature* **459**, 842-846. doi:10.1038/nature08000

Poliak, S., Norovich, A. L., Yamagata, M., Sanes, J. R. and Jessell, T. M. (2016). Muscle-type identity of proprioceptors specified by spatially restricted signals from limb mesenchyme. *Cell* **164**, 512-525. doi:10.1016/j.cell.2015.12.049

Price, S. R., De Marco Garcia, N. V., Ranscht, B. and Jessell, T. M. (2002). Regulation of motor neuron pool sorting by differential expression of type II cadherins. *Cell* **109**, 205-216. doi:10.1016/S0092-8674(02)00695-5

Sherrill, H. E., Jean, P., Driver, E. C., Sanders, T. R., Fitzgerald, T. S., Moser, T. and Kelley, M. W. (2019). Pou4f1 defines a subgroup of Type I spiral ganglion neurons and is necessary for normal inner hair cell presynaptic Ca(2+) signaling. *J. Neurosci.* **39**, 5284-5298. doi:10.1523/JNEUROSCI.2728-18.2019

Surmeli, G., Akay, T., Ippolito, G. C., Tucker, P. W. and Jessell, T. M. (2011). Patterns of spinal sensory-motor connectivity prescribed by a dorsoventral positional template. *Cell* **147**, 653-665. doi:10.1016/j.cell.2011.10.012

Tourtellotte, W. G., Keller-Peck, C., Milbrandt, J. and Kucera, J. (2001). The transcription factor Egr3 modulates sensory axon-myotube interactions during muscle spindle morphogenesis. *Dev. Biol.* **232**, 388-399. doi:10.1006/dbio.2001.0202

Usoskin, D., Furlan, A., Islam, S., Abdo, H., Lönnerberg, P., Lou, D., Hjerling-Leffler, J., Haeggstrom, J., Kharchenko, O., Kharchenko, P. V. et al. (2015). Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencing. *Nat. Neurosci.* **18**, 145-153. doi:10.1038/nn.3881

Vrieseling, E. and Arber, S. (2006). Target-induced transcriptional control of dendritic patterning and connectivity in motor neurons by the ETS gene Pea3. *Cell* **127**, 1439-1452. doi:10.1016/j.cell.2006.10.042

Wang, Z., Li, L. Y., Taylor, M. D., Wright, D. E. and Frank, E. (2007). Prenatal exposure to elevated NT3 disrupts synaptic selectivity in the spinal cord. *J. Neurosci.* **27**, 3686-3694. doi:10.1523/JNEUROSCI.0197-07.2007

Wenner, P. and Frank, E. (1995). Peripheral target specification of synaptic connectivity of muscle spindle sensory neurons with spinal motoneurons. *J. Neurosci.* **15**, 8191-8198. doi:10.1523/JNEUROSCI.15-12-08191.1995

Wu, D., Schieren, I., Qian, Y., Zhang, C., Jessell, T. M. and de Nooij, J. C. (2019). A role for sensory end organ-derived signals in regulating muscle spindle proprioceptor phenotype. *J. Neurosci.* **39**, 4252-4267. doi:10.1523/JNEUROSCI. 2671-18.2019

Yang, X., Arber, S., William, C., Li, L., Tanabe, Y., Jessell, T. M., Birchmeier, C. and Burden, S. J. (2001). Patterning of muscle acetylcholine receptor gene expression in the absence of motor innervation. *Neuron* **30**, 399-410. doi:10.1016/ S0896-6273(01)00287-2

Zhong, J., Li, X., McNamee, C., Chen, A. P., Baccarini, M. and Snider, W. D. (2007). Raf kinase signaling functions in sensory neuron differentiation and axon growth in vivo. *Nat. Neurosci.* **10**, 598-607. doi:10.1038/nn1898

Zhou, X., Wang, L., Hasegawa, H., Amin, P., Han, B.-X., Kaneko, S., He, Y. and Wang, F. (2010). Deletion of PIK3C3/Vps34 in sensory neurons causes rapid neurodegeneration by disrupting the endosomal but not the autophagic pathway. *Proc. Natl. Acad. Sci. USA* **107**, 9424-9429. doi:10.1073/pnas.0914725107