|Home | About | Journals | Submit | Contact Us | Français|
Respiration in mammals relies on the rhythmic firing of neurons within the Phrenic Motor Column (PMC), a motor neuron group that provides the sole source of diaphragm innervation. Despite their essential role in breathing, the specific determinants of PMC identity and patterns of connectivity are largely unknown. We show that two Hox genes, Hoxa5 and Hoxc5, control diverse aspects of PMC development including their clustering, intramuscular branching, and survival. In mice lacking Hox5 genes in motor neurons, axons extend to the diaphragm but fail to arborize, leading to respiratory failure. Genetic rescue of cell death fails to restore columnar organization and branching patterns, indicating these defects are independent of neuronal loss. Unexpectedly, late Hox5 removal preserves columnar organization but depletes PMC number and branches, demonstrating a continuous requirement for Hox function in motor neurons. These findings indicate that Hox5 genes orchestrate PMC development through deployment of temporally distinct wiring programs.
Breathing is a basic motor behavior essential to all terrestrial vertebrates. The frequency and amplitude of respiratory contractions are driven by neural networks residing in the brainstem that coordinate the activation of dedicated sets of spinal motor neurons. Respiratory rhythm generation occurs primarily in the Pre-Bötzinger complex and can be modified by other brain stem nuclei in response to stimuli such as pH changes 1. This rhythm is transmitted via descending pathways to motor nuclei that directly drive the activity of inspiratory and expiratory muscles. Despite the complexity of the networks that regulate respiratory rhythms, contraction of the diaphragm is controlled by a single input supplied by motor neurons within the PMC. Phrenic nerve lesions or spinal cord injuries at or above the fourth cervical segment (C4) result in diaphragm paralysis and respiratory failure, underscoring the vital role of PMC neurons within the respiratory system.
Motor neurons within the PMC are generated in the cervical spinal cord where they form a single clustered population spanning ~3 segments 2. Most PMC axons exit the spinal cord at the C4 level, initially projecting along a medioventral path before converging with other cervical axons at the brachial plexus. Following their separation from limb-innervating axons, PMC axons extend ventrally through the thoracic cavity towards the primordial diaphragm. Upon reaching their target, phrenic axons defasciculate from the main nerve and split into multiple finer branches, prior to forming synapses across the muscle length 3. Although PMC neurons have a central role in respiration, and their columnar organization has been recognized for over 100 years 4, 5, surprisingly little is known about their developmental origins.
All motor neuron subtype identities emerge from the intersection of transcription factor-based programs acting along the dorsoventral and rostrocaudal axes of the spinal cord 6. Motor neurons as a class are produced as an outcome of signaling pathways acting along the dorsoventral axis that specify features common to all subtypes, such as exit of axons from the spinal cord and neurotransmitter phenotype 7. These signaling pathways generate motor neurons that initially express a common set of transcription factors (Hb9, Isl1/2, and Lhx3) which distinguish them from other neuronal classes 8–10. While mutation of transcription factors required for core motor neuron programs results in phrenic nerve loss, largely due to conversion to interneuron fates 10, no selective determinants of PMC identity have been described.
Given their discrete position within the spinal cord, the specification of PMC neurons could involve the same programs contributing to motor neuron diversity along the rostrocaudal axis. Members of the Hox gene family are critical in generating segmentally-restricted motor neuron subtypes at limb and thoracic levels 11. At limb levels, the diversification of lateral motor column (LMC) neurons employs a network of ~20 Hox genes 12, while thoracic level motor neuron fates are determined by the single Hoxc9 gene 13. All Hox gene activities in spinal motor neurons are thought to require the transcription factor FoxP1, as limb-level and thoracic Hox-dependent subtypes are lost in Foxp1 mutants 14, 15. PMC neurons are however not depleted in Foxp1 mutants, but instead appear to increase in number 15. These observations raise the question of whether PMC neurons are specified through mechanisms independent of Hox activities, or whether certain Hox proteins contribute to motor neuron specification independent of Foxp1.
We show here that Hoxa5 and Hoxc5 have critical roles in phrenic motor neuron development. PMC neurons are defined through a broader network of Hox factors that constrain their position and number. Selective deletion of Hox5 genes from motor neurons leads to an extinction of PMC molecular determinants, cell body disorganization, and the progressive loss of PMC numbers. Hox5 genes are also essential for a diaphragm-specific pattern of intramuscular branching, independent of their roles in cell survival. Temporal analysis of Hox5 function in motor neurons indicates that survival and intramuscular branching programs are distinct from those controlling columnar organization. These results define a specific transcriptional program for PMC neurons and indicate Hox activities are required throughout motor neuron ontogeny.
Anatomical studies have identified a column of neurons in segments CIII–CV that projects along the phrenic nerve and innervates the diaphragm 16. To define the molecular identity of this motor neuron group we analyzed transcription factor profiles at cervical levels at embryonic (e) day 11.5 in mice (Fig. 1a–f). Motor neurons within this region expressed combinations of Isl1/2, Hb9, or Lhx3, a core set of transcription factors expressed by all spinal motor neurons (Fig. 1a,c,e). Limb-innervating LMC neurons are distinguished from other subtypes by expression of the transcription factor FoxP1 (Fig. 1b) 14, 15, while medial motor column (MMC) neurons targeting axial muscles coexpress Hb9 and Lhx3 (Fig. 1c,e) 17. We also characterized two additional motor neuron groups at this level; a lateral group that coexpressed Lhx3 and Sox5 (Fig. 1e) 18, and a medial group that expressed Scip (also known as Pou3f1), a POU–class transcription factor (Fig. 1a) 19. Scip+ motor neurons expressed high levels of Isl1/2 and Hb9 (Fig. 1a,c and Supplementary Fig. 1a) and excluded FoxP1 (Fig. 1b). In addition, Scip+ motor neurons expressed ALCAM (Fig. 1d), a cell adhesion molecule expressed by motor neurons at rostral cervical levels 20.
Scip has been suggested to be expressed by PMC neurons 15, and Scip mutant mice are not viable due to respiratory defects 19, although whether cervical Scip+ neurons correspond to the PMC is not known. To assess whether cervical Scip+ motor neurons target the diaphragm we used retrograde labeling assays. We injected rhodamine dextran (RhD) into the phrenic nerve and examined the transcriptional status of motor neurons that accumulated tracer by retrograde transport, using Hb9::GFP mice to aid in the identification of the phrenic nerve (Fig. 1g)8. After tracer injection, retrogradely labeled motor neurons expressed Scip (Fig. 1h and Supplementary Fig. 1c–d) and excluded FoxP1 (Fig. 1i), indicating that Scip expression at this level marks PMC neurons (Fig. 1f)13.
Hox protein activities are critical in motor neuron subtype differentiation along the rostrocaudal axis, and are thus potential determinants of PMC fate. At brachial levels (CIII– CVIII) several genes in the Hox4–Hox8 paralog groups are expressed by motor neurons targeting forelimb muscles12. To determine if a subset of these Hox genes selectively mark the PMC we examined Hox protein expression with respect to Scip+ neurons, and other cervical motor neuron subgroups (Fig. 1j–o). The majority of motor neurons generated at CII–CV expressed Hoxa5 and Hoxc5, including Scip+ PMC, rostral FoxP1+ LMC, and Sox5+ motor neurons (Fig. 1j–k and Supplementary Fig. 1b and data not shown). Unlike LMC neurons, however, Scip+ PMC neurons excluded Hoxc4, Hoxc6 and Hoxc8, although Hoxc6 and Hoxc4 were expressed by dorsal interneurons and other motor neuron subtypes at this level (Fig. 1l–n). PMC neurons can therefore be defined by the selective expression of Scip, ALCAM, Hoxa5 and Hoxc5, and the exclusion of FoxP1, Lhx3, Hoxc6, and Hoxc4.
The restriction of FoxP1 and specific Hox proteins from PMC neurons raises the question of whether their actions contribute to PMC position and distribution. The absence of Hoxc4 and Hoxc6 from PMC neurons suggests a possible repressive influence of these proteins. We therefore examined the number of PMC neurons in mice lacking either Hoxc4 or Hoxc6. While the PMC was grossly normal in Hoxc4 mutants (data not shown), in Hoxc6 mutants we observed a significant increase in the number of motor neurons expressing Scip and ALCAM at e12.5. Specifically we found a ~30% caudal extension of Scip+/ALCAM+ motor neurons and an overall 65% increase in total PMC number (516±38 neurons in wt vs 852 ±44 neurons in Hoxc6 mutants, n=8, P<10−4, Fig. 2a–f). Expression of Hox5 proteins was unchanged in Hoxc6 mutants (data not shown), indicating the expansion is not due to alteration in other Hox genes.
A primary function of Hoxc6 in motor neurons is to promote FoxP1 expression 14, and in Foxp1 mutants Scip expression expands throughout the spinal cord 15. Increased PMC numbers in Hoxc6 mutants might therefore be a consequence of FoxP1 loss. Consistent with this idea, we find that the number of FoxP1+ motor neurons is reduced in Hoxc6 mutants (Supplementary Fig. 2a–b). This observation raises the possibility that the absence of FoxP1 would promote PMC fate, independent of Hox genes. We therefore analyzed how the loss of FoxP1 influences PMC development. To circumvent the early lethality of the global Foxp1 mutation, we analyzed mice in which Foxp1 is selectively deleted from all motor neurons (Foxp1MNΔ) 21. Similar to Foxp1 global mutants 15, Scip expression expanded throughout motor neurons (Fig. 2g–i). Expression of ALCAM, however, was extended only within the Hox5+ domain (Fig. 2j–l), suggesting that FoxP1 exclusion and Scip expression alone are insufficient to specify PMC fate. Furthermore, tracer injection into the phrenic nerve of Foxp1MNΔ mice labeled Scip+ motor neurons that were confined to the Hox5+ domain with no labeling of caudal Hoxc8+ neurons (Fig. 2m–o and data not shown).
To further test the hypothesis that PMC specification relies on the regulation of Scip expression selectively within Hox5+ populations, we analyzed mice that express FoxP1 in all motor neurons using Hb9 regulatory sequences 14. While FoxP1 misexpression had no effect on motor neuron generation or Hox patterns, Scip was extinguished from PMC neurons (Fig. 2p–q). In addition to the PMC, Scip is expressed by pools of LMC neurons (Hoxc8+ Hox5− FoxP1+) targeting the median and ulnar nerves 12. In Hb9::Foxp1 embryos Scip expression is maintained in these pools (Fig. 2r–s). Thus the restriction of Scip to PMC neurons by FoxP1 is constrained to Hox5+ levels. Collectively, these results indicate that Hox5+, Scip+, FoxP1− motor neurons have the capacity to acquire the molecular features and projection characteristics of PMC neurons (Fig. 2t).
To determine whether Hox5 genes are required for PMC development we generated and analyzed mice lacking Hoxa5 and Hoxc5. Because Hoxa5−/– animals display non-neuronal defects in the respiratory system, including abnormal lung development 22, we generated mice in which Hoxa5 was selectively deleted from motor neurons by crossing a conditional Hoxa5 allele23 to Olig2::Cre mice (Hoxa5MNΔ) 24. We verified efficient Hoxa5 removal, finding that by e11.5 Hoxa5 is not expressed in motor neurons but is retained in neighboring cell types and cells of the lung (Supplementary Fig. 3a–f). Hoxa5 MNΔ mice were viable, and did not display discernable respiratory defects. Hoxc5−/− mice were also grossly normal and we therefore introduced the Hoxa5MNΔ mutant allele into a Hoxc5 mutant background. For simplicity we refer to Hoxa5 flox/flox; Hoxc5−/−; Olig2::Cre animals here as Hox5MNΔ mice, as Hoxb5 is normally excluded from motor neurons at this level of the neuraxis (Fig. 1o) and is not upregulated in PMC neurons of Hox5MNΔ mice (Supplementary Fig. 3g–j). At embryonic stages Hox5MNΔ mutants were recovered at the expected Mendelian ratios and showed no external anatomical differences when compared to controls (data not shown). However, all Hox5MNΔ neonates (from n>50 litters) failed to initiate breathing, were cyanotic, and died shortly after birth (Fig. 3a–b). Histological analysis revealed that the lungs of Hox5MNΔ mice were collapsed (Fig. 3c–d), and failed to surface when submerged in water (Supplementary Fig. 3k), indicating that they never inflated with air. Thus Hox5 genes are required in motor neurons for normal respiratory activity.
We next examined how Hox5 mutation affects the emergence of molecular features of PMC neurons. We first assessed the number and the distribution of Scip+ PMC neurons at multiple stages of development. Because Hoxc8 expression does not expand in Hox5MNΔ embryos (Supplementary Fig. 3l–m) 12, we used Hoxc8 as a landmark to define the caudal boundary of the former Hox5 domain and analyzed motor neurons within this region. At e11.5, when motor axons first exit the spinal cord, the number of Scip+ PMC neurons was similar between Hox5MNΔ and control animals (Fig. 3e–f). This number, however, began to decline at e12.5 in Hox5MNΔ mutants, and by e14.5 the number of Scip+ PMC neurons was reduced by 84% (Fig. 3g–l,o). By e17.5 Scip+ motor neurons were undetectable at cervical levels of the spinal cord (Fig. 3m–n). In contrast, the number of limb innervating LMC motor neurons was unchanged in Hox5MNΔ embryos, as assessed by FoxP1 expression at e13.5 (Supplementary Fig. 3n–o). To determine whether the attenuation of Scip expression reflects a selective depletion of PMC numbers, we also assessed the total number of non-LMC motor neurons in rostral cervical segments, by counting Isl1/2+, FoxP1− cells between e12.5–e14.5. This analysis revealed that the total number of non-LMC motor neurons also progressively declined, indicating that in Hox5MNΔ embryos PMC neurons selectively perish (Fig. 3p).
In addition to a reduction in Scip+ motor neurons in Hox5MNΔ mice, the remaining Scip+ neurons were disorganized at e12.5, and intercalated by other subtypes (Fig. 3q–r). However, they retained expression of Isl1/2 and Hb9 (Fig. 3q–r and Supplementary Fig. 3p–q). Expression of ALCAM was reduced at all developmental stages examined (Fig. 3s–t and Supplementary Fig. 3r–s), indicating PMC molecular features are lost prior to motor neuron death. This observation prompted us to examine whether there is an array of genes that are downregulated in PMC neurons in Hox5MNΔ mutants. We performed a microarray analysis of RNA purified from cervical motor neurons in control and Hox5MNΔ embryos, followed by in situ hybridization to verify candidates. Among the genes that showed reduced expression were lynx2, a modulator of nicotinic acetylcholine receptors (Fig. 3u–v) 25 and RTN4/NogoA, primarily known as an inhibitor of neuronal regeneration (Fig. 3w–x and Supplementary Table 1) 26.
Using a candidate approach, we also identified pleiotrophin (PTN), a known target of Hoxa5 that has been shown to have trophic effects in motor neurons 27, 28 as a PMC-restricted marker that is downregulated in Hox5MNΔ mutants (Fig. 3y–z). These observations indicate that Hox5 genes are required for multiple features of PMC fate.
To further assess the impact of Hox5 deletion on PMC development, we monitored the projection of motor axons along the phrenic nerve. We introduced the Hox5MNΔ allele into a Hb9::GFP background and compared projections of control and mutant mice at multiple stages. At e12.5 the phrenic nerve trajectory was similar between Hox5MNΔ and control mice (Fig. 4a–b). However, while all axons were directed to the diaphragm in control mice, in Hox5MNΔ mutants we observed branches straying from the main trunk that appeared to be directed towards the forelimb (Fig. 4b). These branches were transient, as they were not seen at later stages. At e13.5 and later stages we observed a thinning of the phrenic nerve, and by e14.5 the nerve diameter was reduced to 67% of wt (16.47 ±1.9 μm in wt vs 11.04 ±0.3 μm in mutants, P<0.05, Fig. 4c–f), reflecting the progressive loss of PMC neurons.
To assess the molecular identity of motor neurons projecting along the phrenic nerve we performed retrograde tracing experiments in Hox5MNΔ mice. Consistent with a reduction in Scip+ neurons, we observed a 40% decrease in the number of motor neurons that were retrogradely labeled with RhD at e12.5 in Hox5MNΔ mice (Fig. 4g–h,k) and a further reduction at later developmental stages (Supplementary Fig. 4a–d). At this stage labeled cells retained Scip expression and excluded FoxP1, although they were not tightly clustered (Fig. 4g–j). However, by e14.5, when Scip is further downregulated in Hox5MNΔ mice, 85% of retrogradely labeled neurons in Hox5MNΔ mice were Scip–, compared to 35% in control mice (Fig. 4l). These observations indicate a progressive loss of PMC identity in Hox5MNΔ mice, reflecting both a decline of PMC numbers and an altered molecular identity in the remaining motor neurons that have selected a phrenic trajectory.
Because axons extend to the diaphragm in Hox5MNΔ mice we next examined the behavior of PMC neurons at their target muscle. As axons from the phrenic nerve arrive at the diaphragm they split into three main branches, the two largest of which enter the muscle and project in opposing directions along its length. Subsequently, these main branches split further into finer arbors prior to establishing synapses. At e14.5 primary branches appeared to form normally in Hox5MNΔ mice, although there was noticeable absence of secondary and tertiary branches (Fig. 4e–f and Supplementary Fig. 4e–f). By e18.5, intramuscular branches were dramatically reduced in Hox5MNΔ mice, and in the most severe cases only a single branch was observed on one side of the muscle (Fig. 5a–d). Some synapses are formed between the phrenic nerves and the diaphragm as assessed by coincidence of synaptophysin and α-bungarotoxin staining (Fig. 5a–b), but they are localized to a small portion of the muscle. Therefore, the lack of sufficient synapses with the diaphragm is likely the cause of perinatal lethality in Hox5MNΔ mice.
The loss of diaphragm innervation in Hox5MNΔ mice could be a consequence of the absence of a specific molecular program initiated by Hox5 proteins, or simply a manifestation of the loss of PMC neurons. To address this question we introduced the Hox5MNΔ allele into a Bax mutant background to prevent programmed cell death. As reported previously, Bax−/− mice displayed a global increase in motor neuron number, marked by Isl1/2 expression (Fig. 6a–d,m)29. As a result, phrenic neuron number also increased, as evident by Scip and ALCAM expression at e15.5 (Fig. 6a–f,n). Deletion of Bax in a Hox5MNΔ background also dramatically increased Scip+ neurons (Fig. 6d,n). However these motor neurons were scattered (Fig. 6d) and ALCAM expression was not significantly recovered (Fig. 6f). We verified that motor neurons project towards the diaphragm, by DiI retrograde labeling in Hox5MNΔ Bax−/− mice at e18.5 (Fig. 6g–h), indicating rescued cells are not directed to other muscle targets.
An examination of the diaphragm in Bax−/− mice in a wildtype background revealed that the increase in Scip+ PMC neurons does not lead to diaphragm hyperinnervation, consistent with the observation that the overall pattern is not sensitive to elevated motor neuron number 30. Importantly in Hox5MNΔ Bax−/− we did not observe a recovery of the innervation defects seen in Hox5MNΔ mutants (Fig. 6i–l). These results indicate that the defects in Hox5MNΔ mice are not merely a consequence of decreased PMC number, but reflect a specific action of Hox5 genes in determining a molecular program required for diaphragm innervation.
Deletion of Hox5 genes from motor neuron progenitors results in a progressive decrease in the size of the PMC, cell body disorganization, and most strikingly, a loss of diaphragm innervation. These phenotypes become apparent at distinct stages of embryonic development, with cell loss and disorganization occurring at the time of PMC clustering (e11.5–12.5) and innervation defects manifesting during the onset of terminal arborization (e14.5–15.5). Hoxa5 expression persists in PMC neurons until late embryonic stages (e17.5, Supplementary Fig. 6a–b), suggesting a possible late role in PMC development.
To investigate the temporal role of Hox5 genes in PMC neurons we removed Hoxa5 function in a Hoxc5 mutant background by using a Cre recombinase controlled by ChAT regulatory sequences, a gene which is activated shortly after cervical motor neurons differentiate beginning at ~e9.5–e10.5. This strategy effectively removes Hoxa5 protein from motor neurons by e13.5, but preserves expression over a short window between e9.5 and e11.5 (Fig. 7a–d), during the peak of motor neuron generation and as axons select their initial trajectories. Unlike Hox5MNΔ mice, Hox5ChATMNΔ mice are viable, suggesting a less severe defect in PMC development. While the PMC was clustered at e13.5 in Hox5ChATMNΔ mice, there was a progressive reduction in PMC size which was reduced by 48% at e15.5 (Fig. 7a–l,q). Interestingly, while expression of ALCAM was normal at e14.5, PTN expression was attenuated (Supplementary Fig. 6c–f). These results indicate that late removal of Hox5 genes preserves some but not all molecular features of PMC neurons, and that Hox5 genes are required at later stages for maintaining PMC number.
Analysis of Hox5ChATMNΔ mice at e18.5 showed that the diaphragm was fully innervated, with the initial projection patterns similar to control mice. However there was a reduction in the number of terminal branches, resulting in a higher concentration of Ach receptors at individual axon terminals (Fig. 7m–p). While this branching phenotype does not affect viability, it nevertheless indicates that continuous Hox5 activity is required in PMC neurons for the stereotypic pattern of diaphragm innervation.
The neural networks that control breathing rely on the coordinate activation of a select set of muscle groups, which in mammals are supplied by motor neurons within the phrenic and hypaxial motor columns. While the steps that contribute to the specification of hypaxial projecting motor neurons are well documented 9, 14, 31, the programs that distinguish PMC neurons from other motor neuron subtypes are largely unknown. In this study we have found that Hox5 genes are expressed by PMC neurons and are required for multiple aspects of their development, including their organization, survival, and patterns of peripheral connectivity. Unexpectedly, sustained Hox5 activity is necessary for key aspects of PMC maturation subsequent to their initial differentiation. We discuss these findings in the context of Hox networks in the spinal cord and the role of Hox5 genes in respiratory motor neuron specification.
Innervation of the diaphragm is part of a motor circuit unique to mammals, and likely evolved to meet the increased metabolic demands of terrestrial life. How the PMC appeared during mammalian evolution is not known. Terrestrial tetrapods that lack a diaphragm muscle, and therefore a PMC, express Hox5 proteins in domains very similar to mammals 12. In non-mammalian vertebrates, such as chick, Hox5 proteins contribute to the diversity of limb-innervating motor neurons within the rostral half of the LMC 12, a program that requires the Hox cofactor FoxP1 14, 15. However, unlike other Hox-dependent motor neurons, we find that the development of the PMC does not require FoxP1, but rather that FoxP1 inhibits PMC differentiation. These findings indicate a novel FoxP1-independent strategy for motor neuron diversification.
At thoracic levels FoxP1 misexpression inhibits the differentiation of ventrally projecting hypaxial motor column (HMC) neurons, but not dorsally projecting medial motor column (MMC) neurons 14. Thus it is intriguing that the two columns that are selectively inhibited by FoxP1, the HMC and PMC, both project ventrally and are involved in respiratory motor function, suggesting a common origin. At most levels of the spinal cord motor neurons revert to a ground state of the HMC subtype in Foxp1 mutants 14, 15. Our studies indicate that at rostral cervical levels Hox5+ motor neurons acquire molecular features and projection characteristics of PMC neurons in the absence of Foxp1. The PMC therefore likely emerged from an HMC-like precursor that excluded LMC Hox determinants and acquired sensitivity to Hox5 activity.
Deletion of Hox5 genes in motor neuron progenitors leads to a variety of PMC defects, consistent with the idea that Hox genes control diverse facets of motor neuron identity. While Scip+ neurons are generated in Hox5 mutant mice, these motor neurons fail to cluster and expression of ALCAM is dramatically reduced. PMC neurons are progressively lost, and the motor axons that do reach the diaphragm fail to arborize. This innervation phenotype appears to be unique to the diaphragm, as in Foxp1 mutants the majority of limb muscles receive innervation, despite the loss of LMC molecular identity 14, 21. While this suggests a unique program of muscle-specific innervation by PMC neurons, similar phenotypes have been observed in genetic manipulations that affect motor pool specification within the LMC. Mutation in the transcription factor Pea3, for example, leads to motor pool disorganization and loss of intramuscular branches at target muscles 32. Pea3 expression requires peripheral signals from its target provided by glial-derived neurotrophic factor (GDNF) 33. The phenotypes in Hox5 mutants likely reflect a similar defect in nerve-muscle communication, due to an inability of PMC neurons to respond to trophic signals, leading to cell death and nerve degeneration.
Several gene mutations are known to affect diaphragm innervation patterns. For example, mutations in the netrin receptor Unc5c and the neuregulin receptor ErbB2 cause a loss of synapses at the diaphragm 34, 35 and mice lacking ephrins-A2/A5 exhibit impaired topographic innervation 36. Conversely mutations in robo/slit genes increase the number of phrenic branches as a consequence of premature defasciculation 37. Many of these genes are expressed broadly by motor neurons, and whether they contribute to Hox-dependent PMC programs is unclear. Given the variety of defects in Hox5 mutants, it is plausible that Hox5 proteins regulate both PMC-restricted target effectors, such as ALCAM and PTN, as well as determinants shared by many motor neuron subtypes, including RTN4 and lynx2. ALCAM and RTN4 are required for the branching of peripheral neurons 38, 39, while PTN appears to contribute to motor neuron survival28. The concurrent loss of these genes likely underlies the multiple PMC defects observed in Hox5 mutants. Finally, while the precise role of Scip in controlling PMC gene programs is unclear, it likely acts coordinately with Hox5 proteins, as Scip mutants also perish at birth due to respiratory deficiencies.
The actions of Hox factors in motor neuron differentiation are thought to be mediated through the induction of downstream transcription factors, which in turn control distinct aspects of their specification. An early target of Hox proteins is the Foxp1 gene, which is subsequently required for the expression of all LMC and PGC determinants 14, 15. Within LMC motor pools, Hox proteins control the transcription factors Nkx6.1 and Pea3, each contributing to different facets of pool identity. Mutation of Pea3 does not affect the trajectory of motor axons to their target but scrambles cell body positioning and intramuscular branching 32. In contrast, mutation of Nkx6.1 preserves cell body positioning but motor axons fail to project to their correct target 40. These findings raise the question of whether there are Hox-dependent actions in motor neurons that are independent of intermediate transcription factors, and whether Hox function is required subsequent to their induction.
Our findings indicate that Hox5 gene activity is required in PMC neurons subsequent to their differentiation. Motor neuron clustering and Scip expression are not affected under conditions of late depletion using ChAT::Cre, indicating early programs are preserved. However by e15.5 PMC number is decreased by half and terminal branches at the diaphragm are reduced (Supplementary Fig. 7a–b). While it is unclear whether PMC loss is a consequence of reduced branching, a branching defect could in principle limit the accessibility of terminals to trophic molecules leading to motor neuron death41. Alternatively, late Hox5-dependent genes like PTN could act in an autocrine fashion to promote PMC survival after differentiation. Thus, the observed defects could reflect the need for Hox5 genes to maintain expression of effectors required for terminal branching, and preserving communication between motor neurons and muscle.
Our data indicate that Hox5 genes control the branching patterns of PMC neurons, independent of their role in maintaining cell viability. Rescue of the cell death phenotype in the Hox5 mutant fails to restore the normal pattern of diaphragm innervation. While cell number may not directly influence nerve branching, the concurrence of both phenotypes could reflect a common upstream mechanism. For example, peripheral nerve growth factor (NGF) is required both for sympathetic neuron survival and peripheral target organ innervation 42. GDNF, a potent survival factor for some classes of motor neurons, is also required for the correct innervation of the cutaneous maximus and lattisimus dorsi muscles 32, 33. In the absence of Hox5 genes, PMC neurons may lose responsiveness to a muscle-derived signal that regulates both survival and nerve branching.
Some synapses are present at the diaphragm in Hox5 mutants although they are apparently insufficient to drive normal respiratory motor function, as all mutants perish shortly after birth. A contributing factor to the perinatal lethality of Hox5 mutants may be the disorganization of the PMC, which could affect the strength of inputs from hindbrain respiratory networks. Recent studies have suggested that cell position is a critical determinant in shaping the synaptic inputs from premotor sources 21, 43. Loss of Hox5 activity may additionally deplete the inputs from supraspinal networks, thereby exacerbating the defects caused by the loss of diaphragm innervation.
The establishment of connections between phrenic motor neurons and the diaphragm is a key step in the assembly of the neural networks that control breathing. In recent years significant progress has been made in defining the transcriptional codes that contribute to the assembly of the circuits that provide rhythm and modulation to respiratory motor networks 44-46. However there still remain significant gaps in our understanding of basic features of their wiring. For example, the sources of the neuronal subtypes that mediate the connections between respiratory networks to phrenic motor neurons are not clearly defined. Phrenic motor neurons are also known to fire in synchrony with HMC neurons 47, 48, and whether specificity of these connections relies on shared molecular determinants between PMC and HMC neurons is unknown. The definition of a selective molecular code for PMC neurons should aid in attempts to define the mechanisms that shape the specificity of connections between respiratory rhythm networks and the spinal cord.
The Hoxa5 floxed alleles 23, Hb9::GFP 8, ChAT::Cre 49, Olig2::Cre 24, Foxp1MNΔ 21, Hoxc5, Hoxc6 mutant strains 50, Hb9::Foxp1 14 and Bax −/− 51 lines were generated as described. Mouse colony maintenance and handling was performed in compliance with the protocols approved by the Institutional Animal Care and Use Committee of the New York University School of Medicine. Mice were housed in a 12hr light/dark cycle in cages containing no more than 5 animals at a time.
In situ hybridization and immunohistochemistry were performed as described 9. Whole-mount GFP staining was performed as described 40 and motor axons were visualized in projections of confocal Z-stacks (500–1000 μm). Whole mounts of diaphragm muscles from e18.5 mice were stained as described 52. Antibodies were generated as described 9, 12, 14. Other antibodies were used as follows: rabbit anti-GFP (1:1000; Invitrogen), rat anti-ALCAM (1:250; eBioscience), rabbits anti-Neurofilament (1:1000; Millipore), rabbit anti-synaptophysin (1:5; Invitrogen), anti-bungarotoxin Alexa594 conjugate (1:1000; Invitrogen). A Hoxb5 antibody was generated in rabbit using the peptide sequence: ESSRAFPASAQEPRFRQATSSC. For lung histology, lungs were dissected, fixed in 4% paraformaldehyde, paraffin-embedded, sectioned at 4μm and stained with hematoxylin and eosin. Images were obtained with a Zeiss (LSM 510 Meta) confocal microscope and analysed with LSM Image Browser.
Retrograde labeling of motor neurons was performed as described 14. Lysine-fixable dextran-tetramethylrhodamine (RhD, Molecular Probes) was injected into transected phrenic nerves of e12.5-e14.5 embryos. To aid in the identification of nerves, we used GFP fluorescence from Hb9::GFP transgenic mouse embryos, visualized using a MVX10 wide-field fluorescent microscope (Olympus). Nerves were severed using Vannas Spring Scissors (FST) and RhD was injected onto the cut terminal. Embryos were incubated for 5 to 6 hours in oxygenated F12/DMEM (50:50) solution at 32–34°C and subsequently fixed in 4% paraformaldehyde. For labeling of phrenic motor neurons at e18.5, crystals of carbocyanine dye, DiI were pressed onto the left phrenic nerve of eviscerated embryos and the embryos were incubated in 4% paraformaldehyde at 37°C in the dark for 2 weeks. Subsequently, spinal cords were dissected, embedded in 4% low melting point agarose (Invitrogen) and sectioned using a Leica VT1000S vibratome at 100μm.
Cervical spinal cords were dissected from 3 control and 3 Hox5MNΔ mutant embryos in a Hb9::GFP background at e12.5 and motor neurons were sorted by Fluorescence–activated cell sorting (FACS). RNA was extracted using the PicoPure RNA isolation system (Arcturus) and amplified using the Illumina TotalPrep RNA Amplification kit (Ambion). The cRNA from each sample was then hybridized to a MouseWG–6 v2.0 Expression BeadChip (Illumina).
For all experiments n always refers to the number of embryos used. Unless otherwise specified, a minimum of 3 embryos per genotype were used for all reported results. In experiments where motor neuron numbers are reported, cells were counted as motor neurons if they expressed Isl1/2 and as PMC motor neurons if they expressed both Scip and Isl1/2. Statistical analysis was performed by using Student’s t–test (two–tailed distribution, homoscedastic).
We thank Christopher Henderson, Kevin Kanning, and Ivo Lieberam for discussions and sharing unpublished observations, Myungin Baek, Heekyung Jung and Catarina Catela for comments on the manuscript, and Peter Hallock and Steve Burden for assistance with diaphragm preparations. L.J is supported by a grant from the Canadian Institutes of Health Research (MOP-15139). J.S.D is supported by grants from the McKnight Foundation, Alfred P. Sloan, Project ALS, NYSTEM, HHMI and NIH R01 NS062822.
Author Contributions P.P. performed all experiments described apart from the Hoxc6−/− and Hb9::Foxp1 analyses which were performed by J.S.D. C.W. provided technical assistance with in situ experiments and serial sectioning, J.A. performed lung histology, L.J. provided the Hoxa5 flox/flox mouse line, P.P. and J.S.D. designed the study, analyzed the data and wrote the manuscript.
Competing financial interests The authors declare no competing financial interests.