During evolution, motor neurons became specialized to control movements of different body parts including head, trunk and limbs. Here we report that two enhancers of Isl1, E1 and E2, are active together with transcription factors in motor neurons. Surprisingly, E1 and its response to transcription factors has been conserved in evolution from the lamprey to man, whereas E2 is only found in animals with limbs. Our study provides an evolutionary example of how functional diversification of motor neurons is achieved by a dynamic interplay between enhancers and transcription factors.
During evolution, motor neurons became specialized to control movements of different body parts including head, trunk and limbs. Here we report that two enhancers of Isl1, E1 and E2, are active together with transcription factors in motor neurons. Surprisingly, E1 and its response to transcription factors has been conserved in evolution from the lamprey to man, whereas E2 is only found in animals with limbs. Our study provides an evolutionary example of how functional diversification of motor neurons is achieved by a dynamic interplay between enhancers and transcription factors.
Motor neurons are composed of multiple units called motor columns and motor pools specialized to innervate particular peripheral muscles . Motor neurons in the hindbrain innervate and control the movement of head and neck, while somatic motor (sm) neurons in the spinal cord control body muscles. Limbs and body walls are innervated by separate motor columns occupying different mediolateral positions in the ventral spinal cord; thus lateral motor column (LMC) neurons control limbs, preganglionic column (PGC) neurons control sympathetic ganglia, medial motor column (MMC) and hypaxial motor column (HMC, formerly known as MMCl) neurons control the body walls. The diversification of motor neurons is achieved by combinations of transcription factors restricted to particular motor columns: Phox2 factors specify cranial motor neurons and the LIM-homeodomain (HD) factors Isl1 and Lxh3 specify sm neurons [2–5]. Naturally, the interplay between cis-regulatory elements (i.e., promoters and enhancers) and trans-regulatory elements (i.e., transcription factors) is critical for selective gene expression in individual motor neuron subsets. Furthermore, although studies in Hox clusters revealed that signals of body patterning initiate motor neuron diversification, motor neuron-specific transcription control of downstream players in action is not still fully understood [6–8].
Isl1 is a member of the LIM-homeodomain (HD) transcription factor family present in somatic and visceral motor (vm) neurons once they are postmitotic . Genetic and biochemical studies have demonstrated that Isl1 is critical for assigning sm neuron identity, and forms a hexamer complex with Lhx3 [5,9]. Once their pan-motor neuronal identity is acquired via the Isl1-Lhx3 complex, motor neurons further diverge to create multiple motor columns. Motor neurons that retain the Isl1-Lhx3 complex become MMC neurons, while the expression of Foxp1 defines LMC and PGC neurons [10,11]. Isl1 continues to be expressed in most somatic and vm neurons, raising the possibility that dynamic transcriptional control of Isl1 is achieved by differences in the cellular environment of the individual motor neuron subsets. Searches for Isl1 enhancers by comparative functional genomics have revealed multiple cis-regulatory elements (CREs) specific for motor neurons, such as CREST1 and CREST2 identified in zebrafish [12,13]. However, the trans-regulating elements (TREs) that interact with them and the strategy used to achieve accurate spatiotemporal control of subtype-specific enhancer complexes remain unclear.
Interestingly, Isl1 is found in the motor neurons of many animal species, including primitive aquatic animals such as lampreys and ascidians [14,15]. This led us to reason that evolutionary diversification of motor neurons may have occurred along with the transcriptional control of Isl1 activity in newly-defined motor neuronal subsets. Indeed, chordate ascidians contain primitive vm neurons that share molecular characteristics of cranial motor neurons in the vertebrate CNS . Aquatic agnatha (jawless fish) vertebrates such as the lamprey only have sm neurons that contact the body wall, and display traits of MMC and HMC neurons [16,17]. The LMC and PGC neurons arose only later when paired appendages such as limbs (or lateral fins) and a sympathetic nervous system evolved in fish and amphibians [16,18,19]. Thus, motor neurons have constantly developed to expand the repertoire of motor neuron subsets and control novel body parts while the transcriptional control of Isl1 diversified in parallel.
In the present study, we asked whether transcription programs that diversify motor neurons are conserved or change during evolution and, if so, whether motor neurons build new programs when new paired appendages appear. We found that Isl1 expression in motor neurons was mainly controlled by two enhancers, CREST1 and CREST2 (herein called E1 and E2), with the help of the dedicated transcription factors Phox2, Isl1 and Lhx3, and onecut (OC) factor . Chromatin Immunoprecipitation Sequencing (ChIP-Seq) analysis and reporter assays demonstrated that Phox2, Isl1 and Lxh3 induce E1 activity in motor neurons in the hindbrain and the spinal cord, whereas OC-1 selectively induce E2 activity in limb-innervating motor neurons. Comparative genomic approaches showed that the core region of E1 was conserved from jawless fish to humans, whereas E2 was only found in animals with paired appendages. Together our findings demonstrate that motor neuron-specific expression of Isl1 has been conserved in evolution with the help of two major enhancers, and that new strategies were adopted to accommodate newly added paired appendages.
Characterization of Isl1 enhancers that label motor neuronal subsets
To understand the mechanism by which Isl1 becomes selectively expressed in certain neuronal subtypes, we chose to characterize two major Isl1 enhancers originally identified in the zebrafish . To examine their functions in more detail, we generated GFP reporters under the control of the enhancers, and electroporated them into the chick neural tube, aiming at either the hindbrain or spinal cord. There are three different types of motor neurons in the developing hindbrain: branchiomotor (bm), visceral motor (vm) and sm neurons . When the E1::GFP reporter was introduced into the hindbrain, whole mount immunostaining showed that GFP was expressed in the peripheral projections of all types of cranial motor neurons (Fig 1A). Transverse sectioning confirmed that the GFP signal co-distributed with Isl1 immunoreactivity in hindbrain motor neurons (Fig 1B–1D). The mini CMV promoter or minimal Isl1 promoter produced very low level activity by themselves, confirming the specificity of the E1 enhancer (S1A, S1B, S1D and S1E Fig). An E1 reporter with reverse orientation also showed motor neuron-specific activity, consistent with the orientation-independent character of enhancers (S1C Fig). Conversely, the activity of the E2 GFP reporter was not detectable in motor neurons, but only weakly found in sensory ganglia of the hindbrain (Fig 1E–1H).
In the spinal cord, the distribution of GFP-labeled axons in wholemount immunostained embryos and transverse sections showed that the E1 enhancer was active in all motor columns but not in sensory neurons (Fig 1I–1M). E2 enhancer activity was found in sensory neurons, and LMC and HMC neurons but not in MMC neurons (Fig 1N–1R). 3D reconstruction of z-slice images clearly demonstrated that the E1 but not the E2 reporter was active in MMC neurons (S1 and S2 Movies). Together these results show that E1 is active in motor neurons in both the hindbrain and spinal cord, while E2 activity is restricted to subsets of motor neurons in the spinal cord.
Spatiotemporal activation of the Isl1 enhancers occurs in specific motor neuron subtypes
To pinpoint the motor somata that were labeled by the E1 and E2 reporters, we constructed reporters for nuclear GFP (nGFP), which becomes localized to cell bodies. Expression of Isl1, Foxp1, Lhx3 and GFP was assessed by quadruple-immunostaining of individual sections to locate individual motor columns. At brachial levels, MMC neurons express Isl1 and Lhx3, while HMC neurons only express Isl1. LMC neurons are divided into medial LMC (LMCm) and lateral LMC (LMCl), which are Foxp1+Isl1+ and Foxp1+Isl1-, respectively [10,11]. At these levels, E1 was active in all motor neurons (42.0% in MMC neurons; 43.6% in HMC; 13.4% in LMCm; 9.0% in LMCl) (Fig 2A, 2B and 2Q). This pan-motor neuronal activity of E1, with expression even in the LMCl neurons, which do not express Isl1, led us to reason that stable expression of GFP may persist after the enhancer is no longer active. To monitor enhancer activity in situ, we constructed a reporter with destabilized nGFP (ndGFP), whose half-life is less than 4 hours . The majority of cells labeled by destabilized GFP under the control of the E1 enhancer were Lhx3+ MMC neurons (12.8%) rather than LMC neurons (0.2%) or HMC neurons (1.0%) (Fig 2E, 2F and 2R). ndGFP expression labeled a streak of cells next to the pMN domain; these were newborn migrating motor neurons (Fig 2E). In contrast, E2::nGFP expression was mostly found in LMC (12.1% in LMCm; 3.9% in LMCl) rather than MMC (0.3%) neurons, and expression of destabilized GFP was mostly found in LMCm neurons (10.4%) and in MMC neurons (1.9%) but not in LMCl neurons, in good agreement with endogenous Isl1 expression in LMCm neurons  (Fig 2I, 2J, 2M, 2N, 2S and 2T).
Next we examined the enhancer activities at thoracic levels where Isl1+Lhx3+ MMC, Isl1+Lhx3- HMC and Isl1+Foxp1+PGC neurons are present . The E1::nGFP reporter was detected in all motor columns (38.6% in MMC; 30.7% in PGC; 29.5% in HMC), whereas the destabilized E1::ndGFP reporter was restricted to MMCs and PGCs and was barely seen in HMCs (Fig 2C, 2D, 2G, 2H, 2Q and 2R). E2::nGFP reporter activity was found in MMCs (8.0%), PGCs (11.9%) and HMCs (1.3%), while the E2::ndGFP reporter was located in PGCs (1.0%) and HMCs (1.6%) (Fig 2K, 2L, 2O, 2P, 2S and 2T). In summary, E1 activity persists in MMC and PGC neurons and E2 is mostly active in LMCm neurons.
Phox2 homeodomain transcription factors regulate Isl1 expression via the E1 enhancer in bm/vm neurons
We next sought to identify transcription factors in motor neurons that bind to the Isl1 enhancer. When we looked for potential transcription factor binding sites using rVISTA (http://rvista.dcode.org/), we found multiple transcription factor binding motifs within the E1 enhancer, mostly motifs for homeodomain transcription factors. We then carried out luciferase assays using an E1::luciferase construct cotransduced with diverse homeodomain transcription factors and others that are present in motor neurons as follows: Phox2, Isl1, Lhx3, Barx2, Tbx20, Isl2, Nkx6.1, Hb9, Meis1, Pbx1, Otp, Shox2, Hmx3, Pax6, Sip and OC-1 [23–29]. Of the various homeodomain transcription factors, only Phox2a/b and Lhx3 induced expression of the E1 reporter gene (S2 Fig).
Phox2a and Phox2b are two paralogous homeodomain transcription factors present in branchio-and visceromotor (bm/vm) neurons in the hindbrain [2–4]. Misexpression of Phox2 factors induces ectopic expression of Isl1 when electroporated into the spinal cord, indicating that Phox2 factors may activate Isl1 enhancers [2,30]. We first explored ChIP-Seq data for Phox2 genomic binding sites in iNIP cells, an inducible embryonic stem cell line with bm/vm neuronal properties  (Fig 3A and S3 Fig). Remarkably, the highest Phox2 ChIP-Seq peak between 620 kb downstream and 540 kb upstream of the Isl1 locus was in E1, around 220 kb downstream of the Isl1 transcription start site. To further define Phox2 binding sites within E1, a series of deletions and point mutations was introduced into E1 (Fig 3B and S5A Fig). In luciferase assays, Phox2a induced E1 activity by 4.13-fold but did not induce E2, suggesting its specificity to E1 (S2B Fig). A fragment containing nt 320–479 of E1 proved to be sufficient for activation by Phox2a, whereas the mutE1-3 and mutE1-4 sites abolished activation (Fig 3B and S5A Fig).
We and others previously showed that overexpression of Phox2 factors induced ectopic expression of Isl1 within the chick spinal cord, and that the ectopic Isl1+ cells had bm/vm characteristics determined by their expression of the bm/vm marker Tbx20 (S4F–S4H and S4K–S4M Fig) [2,30]. The Isl1+ cells induced by Phox2a, however, did not express markers for other neuronal classes such as Brn3a (for dI3 neurons) and MNR2 (for sm neurons) (S4I, S4J, S4N and S4O Fig). Thus, we decided to use Isl1+ as a phenotypic readout for the ability to induce hindbrain motor neuron identity. When deletions and point mutations were introduced in E1, GFP expression in the spinal cord was not changed except in the cases of nt320-470 E1 and mutE1-3 (S4P–S4V and S4Y Fig). Nt320-470E1 was not active in spinal cord motor neurons because the Isl1-Lhx3 binding motif was partially deleted (S4Q and S5A Figs). Nevertheless it remained active in the hindbrain, implying that the fragment is involved specifically in bm/vm neuron expression (S4X Fig). When Phox2a was electroporated, dorsal expansion of GFP reporter activity was found in nt 320–479 E1, mutE1-1 and mutE1-2, but not in mutE1-3 and mutE1-4 (Fig 3G–3K and S5A Fig). MutE1-3 GFP expression was entirely lost in the neural tube, indicating that the mutE1-3 site is also essential in spinal cord sm neurons (S5A Fig). On the other hand, mutE1-4 only lost its activity in the dorsal spinal cord, in which bm marker Tbx20 but not Brn3a or MNR2 was induced, but its GFP expression was intact in sm neurons (Fig 3K and S4K–S4O Fig). This implies that Phox2 binding to the mutE1-4 motif turns on Isl1 selectively in hindbrain motor neurons (mostly branchiomotor/visceral) but not in spinal cord motor neurons (mostly somatic). Axonal projection in wholemount and transverse sections of chick embryos also confirmed the restricted activity of mutE1-4 in hindbrain sm neurons (S6A–S6H’ Fig). We also demonstrated specific binding of Phox2 at the E1-4 site by gel shift assays and chromatin IP (Fig 3M and 3N). We next tested whether Phox2 factor is necessary for inducing E1 activity in vivo. When we knocked down chick Phox2b by siRNA, E1 activity was diminished in bm neurons (Fig 3O–3Q and S7A Fig). Together, our results suggest that Phox2 factor is necessary and sufficient to induce E1 in hindbrain bm neurons.
Next we generated stable transgenic embryos carrying E1 and mutE1-4 GFP reporters. GFP-labeled peripheral projections were found in the hindbrain and spinal cord of wholemount E1::GFP embryos (Fig 4A–4B’). In transverse sections of these embryos GFP expression was found in Isl1+Phox2b+ facial motor neurons of rhombomere (r) 4, indicating that E1::GFP is active in bm/vm neurons (Fig 4C, 4C’, 4E and 4E’). Sm neurons in the caudal hindbrain also expressed E1::GFP (S6E, S6G and S6G’ Fig). Wholemount GFP staining of mutE1-4::GFP showed that GFP was expressed in the peripheral projections of the spinal cord but not in the hindbrain (Fig 4B and 4B’). Transverse sections of the hindbrains of mutE1-4::GFP embryos also showed that GFP expression was absent from facial motor neurons and other bm/vm neurons (Fig 4D, 4D’, 4F, 4F’, S6D and S6F Fig). However, the sm neurons in the caudal hindbrain expressed mutE1-4::GFP (S6F, S6H and S6H’ Fig). In the spinal cord, GFP expression by both E1 and mutE1-4 reporters was found in sm neurons at brachial and thoracic levels (Fig 4G–4J’). Thus, Phox2-E1 interaction is required for specific expression of Isl1 in bm/vm neurons.
Lhx3 and the Isl1-Lhx3 complex activate the E1 enhancer within somatic motor neurons
Since the E1 enhancer is active in spinal motor neurons and Phox2 is not present in the spinal cord, transcription factors other than Phox2 presumably control its activity in the spinal cord. Our luciferase assays with diverse transcription factors showed that Lhx3 also activated the E1 enhancer (See S2 Fig). Lhx3 is present in p2 and pMN domains and participates in producing V2a interneurons and motor neurons, respectively (Fig 5A–5C) . Isl1 appears when pMN progenitors just become postmitotic, which assigns pan-motor neuron identity by forming a hexamer complex with Lhx3 (Fig 5B and 5C) [5,32]. When motor neurons further diverge into multiple motor columns, only MMC neurons co-express Isl1 and Lhx3 [1,33]. In line with this, we detected E1::ndGFP reporter activity in cells that co-expressed Isl1 and Lhx3 (see Fig 2E, 2F and 2R). To test whether enhancer activity in the spinal cord was altered in the presence of Isl1 and Lhx3, we electroporated cells with the E1::GFP reporter together with Isl1 and Lhx3. When Isl1 alone was electroporated, E1 enhancer activity remained restricted to motor columns, as in the GFP controls (Fig 5D). In contrast, when Lhx3 was electroporated, E1 enhancer activity was expanded into the dorsal column where ectopic Chx10+ V2a interneurons arose, although to a lesser degree than the group electroporated with Isl1 and Lhx3 (Fig 5E and 5V). Co-electroporation of Isl1 and Lhx3, or introduction of DDI1L3, which mimics the Isl1-Lhx3 complex, also resulted in expansion of E1 activity in the dorsal spinal cord, in which ectopic MNR2+ motor neurons were induced (Fig 5F–5H’ and 5V) . We also tested whether misexpression of LIM factors induced transcription of the endogenous Isl1 gene in the dorsal spinal cord where the GFP reporter was induced. Since our Isl1 antibody does not distinguish between chick and mouse Isl1 (they are 99% similar), we used Isl2 instead of Isl1, a paralogue of Isl1 with similar biological activity [34,35]. No E1 reporter activity was found in the presence of Isl2 alone (Fig 5I and 5V). In contrast, when Isl2 and Lhx3 were introduced, ectopic Isl1 appeared in the dorsal spinal cord, consistent with the expansion of GFP expression driven by E1 (Fig 5J–5L and 5V). To determine whether the effect of Lhx3 required the LIM domain, we functionally blocked that domain by co-electroporating ΔL-Lhx3, LMO4 or the dimerized domain (DD) . The induction of the E1 enhancer was blocked in all three conditions, indicating that LIM domain-based complex formation is required for Lhx3 to activate the E1 enhancer (S8 Fig). Luciferase reporter assays also showed that E1 was induced by Lhx3, Isl1+Lhx3 and DDI1L3 but not by Isl1 (Fig 5W).
Repressors prevent abnormal activation of E1 in progenitors and V2 interneurons
Although Lhx3 is present in pMN progenitors and activates E1, we did not observe E1 activity in the pMN region (Fig 5A). Only a streak of migrating motor neuron progenitors that had just left the pMN domain expressed the E1::GFP reporter, and this coincided with the appearance of a low level of Isl1 protein (Fig 5B–5B”‘). Thus, some unknown repressor may suppress the initial expression of Isl1 in pMN progenitors and this repression may be released around the time when the progenitors become postmitotic motor neurons . To test this hypothesis, we tested the effect of Olig2 on Lhx3 expression since Olig2 is present in the pMN domain, is required for motor neuron specification and is extinguished in postmitotic motor neurons [36,37]. In HH stage 23 embryos, differentiating progenitors at lateral region of the pMN domain co-expressed Olig2 and Lhx3 but lacked E1::GFP activity (Fig 5A–5A”‘). As expected, co-electroporation of Olig2 with Lhx3 suppressed the activity of the latter as an inducer of E1 (Fig 5N, 5O and 5V). The induction of E1 luciferase reporter activity by Lhx3 was also inhibited by Olig2 (Fig 5W). Conversely, Olig2 was not effective in reducing E1::GFP expression in endogenous motor neurons or E1 GFP activity driven by Isl1 and Lhx3 (S5C Fig). Thus, Olig2 selectively blocks activity of Lhx3 complex. This was abolished when putative E-box sequence in E1 was mutated, indicating that Olig2 may bind to this site (S5A and S5K–S5M Fig).
Since Lhx3 is also present in the p2 domain where V2 interneurons arise, the E1 enhancer might also be expected to be active in the p2 domain and V2 interneurons. However, this was not observed in our experiments, raising the possibility that an unknown repressor suppresses E1 activity in the p2 domain. Therefore we tested the effect of two known repressors, Sox1 and Chx10, present in the p2 domain and postmitotic V2a interneurons, respectively [38–40]. When electroporated, Sox1 and Chx10 were effective in repressing the ectopic induction of GFP reporter activity driven by Lhx3 (Fig 5Q, 5R, 5T, 5U and 5V). There are two predicted Chx10-binding motifs in E1 but mutating them did not abolish the repressive activity of Chx10 (S5A and S5N–S5S Fig). Since Chx10 binds to AT-rich nucleotides and is known to block hexameric Isl1-Lhx3 at the Hb9 promoter, we decided to test whether the Isl1-Lhx3 binding motif E1-3 that we had identified was required for repression by Chx10 . Because mutating the E1-3 site abrogated E1 activity, we could not test the repressive activity of Chx10 using mutE1-3. Instead, we used tandem repeats of the E1-3 site (6xE1-3) for GFP reporter and found that its induction by Lhx3 or Isl1-Lhx3 was blocked by Chx10 (S5T–S5V, S5X and S5Y Fig). This appeared to require DNA binding since the Chx10 N51A point mutant defective in DNA binding failed to repress the induction of 6xE1-3 (S5W Fig) . In the case of Sox1, ΔNSox1 was still potent in its repression but C-Sox1 lacking the DNA binding domain was not (S5G–S5J Fig) . Thus, DNA binding is also required for Sox1-mediated repression. The E1 reporter was still active in motor neurons in the presence of Sox1 or Chx10, indicating that the repressors are only effective in non-motor neurons. Nevertheless, Chx10 but not Sox1 was able to inhibit the induction of E1::GFP reporter by exogenous Isl1 and Lhx3 (S5D and S5E Fig). We conclude that E1 activity is repressed by Sox1 in p2 progenitors and by Chx10 in V2a interneurons.
Isl1 and Lhx3 bind to E1 enhancer in somatic motor neurons
To search for Isl1-Lhx3-binding motifs in the E1 reporter, we examined the results of ChIP-Seq in NIL cells, inducible embryonic stem cells with sm neuronal traits . Strong binding of Isl1 and Lhx3 occurred in E1 (Fig 6A). Of the luciferase reporters with point mutations, all were induced by Isl1 and Lhx3 except mutE1-3, and the activity of all the GFP reporters except mutE1-3 appeared ectopically in the dorsal spinal cord when Isl1 and Lhx3 were co-electroporated (Fig 6B–6H and 6J–6L). E2 was not induced by Isl1 and Lhx3, suggesting the specificity of E1 (S2B Fig). GFP expression of mutE1-3 was completely lost even in the ventral spinal cord, indicating that E1-3 site activity relies on the endogenous Isl1 and Lhx3 present in motor neurons. MutE1-3 reporter activity was present in bm/vm neurons in the hindbrain when the reporter was introduced by electroporation, indicating that the E1-3 site is sm neuron-specific and its alteration by mutation does not affect basal transactivating ability (S6I and S6J Fig). To test whether the E1-3 site was sufficient to drive sm-specific gene expression, we generated a GFP reporter with six tandem repeats of the E1-3 site (6xE1-3). Expression of the 6xE1-3 reporter was selective in sm neurons by itself and expanded dorsally in the presence of Isl1 and Lhx3, indicating that Isl1-Lhx3 binding to the E1-3 site is sufficient to drive motor neuron-specific gene expression (Fig 6I, 6L and S4V Fig). Together our results suggest that Isl1-Lhx3 complex activates E1 in motor neurons.
OC-1 binds to E2 enhancers to drive LMCm-specific expression of Isl1 in tetrapods
E2 sequences are conserved in mouse, chicken, zebrafish and fugu, but no E2 sequences were identifiable in the non-tetrapod chordate, lamprey and the cephalochordate amphioxus, at least in any currently available database [42–44] (S9 Fig). This indicates that E2 might have appeared in genomes when limb structures arise in the early vertebrates. All the E2 sequences derived from mouse, chicken and zebrafish were active in motor columns of the ventral spinal cord (Fig 7A–7C). We also generated putative E2 reporter from fugu and found that it was specific for motor neurons (Fig 7D). To search for E2-binding transcription factors, we first examined ChIP-Seq data for Phox2, Isl1 and Lhx3 binding in embryonic stem cell lines but found no significant peaks in the E2 region (S9 Fig) . In a candidate approach, we decided to test the activity of OC factors, previously known to promote LMC identity and bind to E2 in ChIP assays . When the OC-1 factor was co-electroporated, E2::GFP reporters became active in the entire spinal cord (Fig 7H–7L and 7P). OC-1 did not induce E1::luciferase reporter as previously reported (S2B Fig) . Using the UCSC genome browser, we searched for two putative OC factor-binding sites in the E2 enhancer that were highly conserved in fugu, stickleback, zebrafish, coelacanth, chicken, mouse and human (Fig 7R). MutE2-1 and mutE2-2 retained their activity in motor neurons but only mutE2-1 responded to OC-1 (Fig 7E, 7F, 7H, 7M, 7N and 7P). When both sites were mutated in mutE2-3, the reporter was not induced by OC-1 (Fig 7G, 7O and 7Q). Similar responses were observed in luciferase assays in that the induction of E2 activity by OC-1 was significantly downregulated in mutE2-3 (Fig 7S). We also tested whether downregulation of OC factors affect E2 activity or not. When we knocked down OC-1 and OC-2 by siRNA, E2::GFP activity was significantly reduced (Fig 7T–7V and S7B–S7D Fig). Similarly, when we misexpressed Hoxc9 that switchs the LMC identity to thoracic motor neurons as determined by reduced expression of Foxp1 and Raldh2, E2 activity was also downregulated (S7E–S7L Fig) . Together, our results suggest that OC factors are necessary and sufficient to induce the LMC-specific activity of the E2 enhancer.
Evolutionary diversification of motor neuron-specific enhancers
It has been shown that the E1 enhancer sequence is strongly conserved from fugu to human. We confirmed this, and found that the Phox2-binding and Isl1-Lhx3-binding motifs we identified were also highly conserved (Fig 8A). To examine the vertebrate origin of motor neurons, we tested for the presence of the E1 enhancer in the lamprey, a living representative of the most ancient vertebrates [42,43]. Although there was generally little conservation of the entire E1 sequences, we found that the core region of E1 containing the Phox2 and Isl1-Lhx3 binding sites were relatively well conserved (Fig 8B and 8C). To test whether enhancer activity was conserved, we generated GFP enhancers from the E1 sequences of mouse, chick, zebrafish and lamprey and electroporated them into the chick spinal cord. All the reporters showed bm neuron-specific GFP expression in the chick hindbrain (Fig 8D, 8H, 8L and 8P). The mouse, chick and lamprey E1 enhancers were active in the motor neurons of the spinal cord, and expanded in the presence of Phox2 or Isl1-Lhx3 (Fig 8E–8G, 8I–8K and 8Q–8S). However, zebrafish E1 was not active in the spinal cord motor neurons but nevertheless was induced by mPhox2 and zIsl1-zLhx4 in the dorsal spinal cord (Fig 8L–8O). The same analysis was applied to the E2 enhancer sequence. In contrast to Uemura et al, we detected the E2 enhancer sequence with conserved OC binding motifs in fugu . We were unable to retrieve any putative E2 sequence from the lamprey genome when we conducted a BLAST search using E2 sequences from human, mouse, chicken, coelacanth, zebrafish, stickleback, and fugu . We also found no match in amphioxus genomic data . To be sure of the absence of E2 in the lamprey genome, we also analyzed vertebrate basewise conservation scores (phyloP) using the newest version, mmc10, which includes the lamprey genome (S9 Fig). Together, we conclude that motor neuron-specific E1 activity in the CNS is well-conserved from lamprey to man, whereas E2 appeared first in fugu along with the origin of limb/fin structures.
The interplay between the E1 enhancer and transcription factors occurs in parallel with the diversification of hindbrain and spinal cord motor neurons
Motor neurons consist of multiple motor columns that innervate distinct muscle targets. Each motor column expresses some combination of transcription factors that form the basis for diversification of motor neurons [1,33]. Naturally the interplay between cis-regulatory elements (i.e., promoters and enhancers) and trans-regulatory elements (i.e., transcription factors) is critical for selective gene expression in individual motor neuronal subsets. In the present study, we focused on transcriptional regulation of Islet1, one of the representative markers of motor neurons that define motor neuron identity. Using stable and destabilized GFP reporters introduced into chick embryos, we demonstrated that E1 and E2 are differentially expressed in motor neuron subsets: E1 is active in all motor neurons and later becomes restricted to cranial motor neurons and MMC and PGC neurons in the spinal cord; E2 is more selective, being active in limb-innervating LMCm neurons that express Isl1, in line with previous reports in the zebrafish  (Fig 7W).
In the present work we have further characterized the transcription factors that interact with enhancers and found that Phox2 binds to E1 and drives Isl1 expression in the hindbrain. Phox2a is necessary and sufficient to drive E1 activity since exogenous Phox2a induced E1 activity in embryos and cell lines, and knockdown of Phox2b by siRNA abrogated E1 activity. We designated as E1-4 the Phox2a binding site in E1 whose mutation abolished its activity in bm neurons but not in sm neurons. The mutE1-4 enhancer was not responsive to exogenous Phox2a but responded well to exogenous Isl1 and Lhx3, all of which supports the existence of bm neuron-specific gene transcription mechanisms. Similarly, we have demonstrated that the Isl1-Lhx3 complex induces E1 in sm neurons via the E1-3 site. Mutations at the E1-3 site abolished expression in sm neurons and prevented induction by Isl1-Lhx3. It is noteworthy that the mutE1-3 construct did not respond to exogenous Phox2a but was active in native bm neurons. Perhaps the cellular environment of the ectopic bm neurons was not as favorable as the native ones, and additional factors present in bona fide bm neurons may bind to this site and facilitate the action of Phox2a on E1. A likely candidate for such a factor is Isl1; the idea that it synergises with Phox2a is suggested by the surprisingly extensive overlap between binding sites for Phox2a and for Isl1 in a genomewide ChIP-Seq analysis . Alternatively, Phox2a may act together with bHLH factors, as has been suggested to occur in sympathetic neurons .
Previous studies by us and others have demonstrated that Isl1 switches the stoichiometry of the tetrameric Lhx3 complex to that of a hexameric Isl1-Lhx3 complex to induce motor neuron identity, and aberrant assembly of the Lhx3 complex in motor neurons is prevented by repressors such as LMO factors and Hb9 [5,32,35]. In this study, we unexpectedly found that the Lhx3 complex is capable of activating E1, which may contribute to the initiation of Isl1 transcription when motor neurons are about to become postmitotic. Nevertheless, ectopic expression of Isl1 in Lhx3-expressing motor neuron progenitors and V2 interneurons was blocked by the repressors Olig2, Sox1 and Chx10. These repressors were thus effective in blocking the activation of E1 by Lhx3. Interestingly, the ectopic formation of motor neurons and the E1 activity induced by Isl1-Lhx3 were also blocked by exogenous Chx10 but not by Olig2 and Sox1, suggesting that the potencies of these repressors or their mechanisms of action were different; repression may require DNA binding or protein-protein interaction. For instance, Olig2 may act as a repressor by binding to the E-box element, or it may squelch other bHLH factors such as Ngn2 by protein-protein interaction; either way it suppresses motor neuron formation . In our study, we mapped the E-box site in E1 and found that mutating it made E1 insensitive to Olig2. Thus, DNA binding of Olig2 is important for its repressive action on E1. In the case of Chx10, mutating two consensus Chx10 binding sites had no effect. Interestingly, it is reported that the binding sites for Chx10 and hexameric Isl1-Lhx3 are similar and therefore Chx10 may inhibit Isl1-Lhx3 binding to the Hb9 promoter . In agreement with this we showed that Chx10 inhibits the transcriptional activation of tandem repeats of the E1-3 site by Isl1 and Lhx3. The point mutant Chx10 N51A deficient in DNA binding failed to inhibit the activity of E1-3 site, implying that Chx10 may bind to this site. What then would be the explanation for selective inhibition of Lhx3 complex activity by repressors? Our results in chick embryos and cell lines showed that the Lhx3 complex is less potent than the Isl1-Lhx3 complex in activating E1. This could be due to the lower DNA binding affinity of the Lhx3 complex demonstrated previously . In addition, the binding elements for the Lhx3 and Isl1-Lhx3 complexes appear to be qualitatively different, and the different arrangements of multiple repeats and spacing between them should make the difference even greater, and may render the Lhx3 complex more susceptible to repressors .
It is noteworthy that Isl1 alone is not effective in inducing motor neuron identity or activating the E1, unlike Lhx3 or Isl1-Lhx3. Structure-function analysis suggests that NLI favors its binding to LIM domains of Lhx3 over Isl1, indicating that Lhx3 could be more efficient in forming a tetrameric complex in vivo . In addition, the Lhx3 binding domain (LBD) of Isl1 at the C-terminus binds to Lhx3 to form hexameric Isl1-Lhx3 complex, which makes tetrameric Lhx3 or Isl1 complexes less available . Recently, it is reported that the LBD domain also interacts with the LIM domains of Isl1 [50,51]. The intramolecular interaction is weak but specific, preventing unnecessary DNA binding and facilitating cofactor exchange in the cell. Indeed, Isl1 is known to act in synergy with other nuclear factors in neurons and other tissues [52–54]. Together, more diverse choice of cofactors and its tendency to suppress its own binding to DNA may explain weak biological activity of Isl1 by itself in motor neurons.
Considering the evolutionary emergence of transcription factors, it is noteworthy that Phox2 and Islet are found even in Ciona internalis in which Mnx-expressing somatic and Phox2-expressing visceral neurons are distinguishable [14,55]. Phox2 and Islet factors are also found in the lamprey and amphioxus, whose motor neurons are primitive, lacking motor columns and motor pools [56–58]. Thus, Phox2 and Isl1-Lhx3 were used in the motor neurons of ancient aquatic animals even before they diversified. In line with this, we identified Phox2 and Isl1-Lhx3 binding motifs in E1 enhancers from lamprey to man. Anatomical and genetic comparisons between species have suggested that the boundary between the hindbrain and spinal cord motor neurons is defined even in amphioxus (cephalochordate) and lamprey (chordate) [59–61]. Hence it appears that the motor neuron-specific activation of Isl1 with the help of Phox2 and Isl1-Lhx3 is evolutionarily conserved.
In this study, zebrafish E1 was not active in chick spinal cord motor neurons nor was it activated by exogenous Isl1-Lhx3 factors from other species. Nevertheless, we found that it was induced by zebrafish Isl1-Lhx4. There was only a one base difference between the Isl1-Lhx3 binding site in zE1 and the sites in mouse and chick. Substituting the binding site in mE1 for zE1 did not make it responsive to mIsl1 and mLhx3. It is possible that species differences in LIM-HD factors account for the absence of zE1 activity in chick spinal cords. However, the LIM-HD factors of zebrafish and chick differ by only one or two amino acids. In fact substantial cross-species variation has been reported in protein-DNA interactions that cannot be simply explained by changes in protein or DNA sequences [62,63]. Since an additional whole genome duplication event occurred in ray-finned fish including zebrafish, it is possible that extensive rearrangements of multiple cis- and trans-elements may have contributed to the zebrafish-specific transcription program.
The E2 enhancer and OC factor define limb-innervating motor neurons in tetrapod animals
We and others have shown that the E2 enhancer is selective for LMCm neurons in the limb-innervating motor columns that express Isl1 . Our detailed fate mapping analysis showed that E2 enhancer activity is strong in LMCm neurons. We and others also showed that E2 responded to OC-1 factors, which are important for LMC-specific gene expression . In mice deficient in OC factors there are fewer LMCm neurons and more LMCl and PGC neurons as a result of a reduced level of Isl1 protein . In line with this, we demonstrated that knockdown of OC factors or disruption of LMC identity by manipulating the Hox code coincided with reduced E2 activity. Taken together, these findings suggest that OC factors and E2 activity are important for maintaining Isl1 expression in LMCm neurons.
LMC neurons were added to motor columns when limb/fin structures appeared during evolution, and this coincided with the appearance of the Hox co-factor FoxP1 and loss of Lhx3 [16,19]. FoxP1 interacts with Hox proteins, which were altered in parallel with the emergence of paired fins/limbs [11,64]. Foxp1 also suppresses Lhx3, and this suppression is a prerequisite for LMC identity [10,11]. Likewise, in the absence of Foxp1, LMC neurons display HMC-like properties, reflecting their evolutionary origin . Nevertheless, LMCm neurons maintain Isl1 expression in the absence of Lhx3, and thus in the absence of the Isl1-Lhx3 complex. Hence it is likely that tetrapod animals adopted the alternative enhancer E2 together with OC factors to maintain Isl1 expression in LMCm neurons due to the absence of Lhx3.
Gradual specification of motor column identity
We used two reporters, one expressing stable GFP and the other, transient destabilized GFP, to trace cells that had activated the reporter at least once during their development and cells that were currently activating the enhancer, respectively. In this way, we expected to see whether motor column identity was acquired progressively or not. We found that the stable E1::GFP reporter marked all motor neurons, starting from the initial E1 activity in newborn motor neurons that had just begun to express Isl1. Thus, the transient assembly of Isl1 and Lhx3 in order to initially acquire motor neuron identity occurs in all motor neurons subsets, as previously suggested . On the other hand, the destablized E1 reporter representing current E1 activity only labeled MMC and PGC. This suggests that E1 continues to be active in MMC and PGC, but that the initial activity of E1 disappears in the other motor columns. LMC neurons appear to initially use the E1 enhancer to obtain pan-motor neuron identity but later shut off E1 activity when they lose Lhx3. In the case of E2 reporters, the stable and transient E2 GFP reporters were both expressed in motor neurons once they reached a lateral position, mostly in LMCm neurons. Thus, the latter may use E1 and E2 sequentially to drive Isl1 expression during their transition from pan-motor neurons to LMC neurons. LMCl neurons acquire their identity when they pass through the LMCm area, in which they are exposed to retinoic acid signaling [66,67]. This induces Lhx1, which suppresses expression of Isl1 and specifies the LMCl fate. Thus, E2 may be initially active in all LMCs but later cease to be active in LMCl neurons when they migrate more laterally, with the help of Lhx1 and retinoids. Thus, gradual changes in the interactions between transcription factors and enhancers lead motor neurons to acquire their columnar identity progressively.
We have seen that E1 has been evolutionarily conserved since the era of the lamprey, a jawless vertebrate which possesses primitive motor neuron . Recent comparative genomic studies have demonstrated that the segmentation of hindbrain and spinal cord is conserved in the lamprey, and is under the control of patterning signals such as Hox clusters [60,64,68]. The hindbrain and spinal cord motor neuron-specific expression of Islet1 could be achieved by E1 in the lamprey as part of the patterning program. However, the E2 enhancer, which is restricted to limb-innervating motor neurons, appeared when animals such as fugu and tetrapods developed limbs. E2 evolved to cooperate with OC factors and drive Isl1 expression in LMC neurons. Thus, our observations demonstrate that evolutionary conserved enhancers drive motor column-specific gene expression during motor neuron development.
Materials and Methods
Isl1 enhancers were amplified by PCR with primers using genomic DNA from mouse, chick, zebrafish, fugu and lamprey (S1 Table). PCR fragments were subcloned into the following reporter vectors: pCS2 mini CMV-GFP, pCS2 mini CMV-luciferase and tk-luciferase reporter vector . pCS2 mini CMV-GFP/luciferase contains a 60 bp TATA box and the transcription initiation site of the cytomegalovirus (CMV) promoter (GenBank accession no. X03922: nt 1090–1149). The mini CMV promoter and eGFP sequences were obtained by PCR amplification from pEGFP-N1 (Clontech). Mutations were introduced in E1 or Chx10 by PCR-based mutagenesis. Zebrafish and Fugu E2 DNAs were amplified from Dario rerio and Fugu rubripes, respectively.
In ovo electroporation was performed as described previously . In brief, about 1 μg/μl of DNA solution was electroporated into the chick spinal cord at Hamburger and Hamilton (HH) stages 12 and harvested at HH stages 23 to 29. In the case of OC-1, HH stage 14 embryos were used. For the hindbrain electroporation, DNA solution was injected into HH stage 8 to 10 embryonic neural tube anteriorward from the level of approximately the third somite as previously described . Electroporation was performed with a condition of 5 times of pulses, 20–24 volts, 50 msec, 1 sec intervals.
Immunohistochemistry and in situ hybridization
Immunohistochemistry was performed as described previously . Following antibodies were used: mouse anti-Isl1 (DSHB), mouse anti-Isl2 (DSHB), rabbit anti-Isl1/2 , mouse anti-MNR2 (DSHB), rabbit anti-GFP (Invitrogen), rabbit anti-Foxp1 (Abcam), rabbit anti-Sox1 (Cell signaling), guinea pig anti-Lhx3 , guinea pig anti-Chx10 , guinea pig anti-Olig2  and mouse anti-Neurofilament (DSHB). For wholemount immunostaining, Day 4.5 chick embryos electroporated at Day 2 were fixed and incubated with primary antibodies for 3 days and secondary antibodies for 1 day. For in situ hybridization, transverse sections were hybridized with digoxigenin-labeled probes specific for mouse Phox2a (full CDS), mouse OC-1 (partial CDS, 469–1263 bp). All images were captured with epifluorescent microscope or confocal microscope (Zeiss).
Cells in each quadrant of the ventral horn were counted on z-series of slice images using a confocal microscope and ZEN2009 imaging software (Zeiss). At least 3 embryos were quantified from each group, and 3–4 images were collected from each spinal cords. %GFP/column was calculated as the percentage of the number of GFP-expressing cells among the motor neurons in each column. To quantify GFP intensity induced by exogenous transcription factors in the chick spinal cords after electroporation, 12 μm-thick transverse sections were immunolabeled with GFP. The background-subtracted pixel intensities of GFP in 120 x 240 μm2 areas in the dorsal spinal cord were measured using ImageJ. At least 10 sections from 4 embryos were analyzed for each group. Statistical significance was analyzed by unpaired Student’s t-test and the Kruskal-Wallis test for multiple comparisons.
E1::GFP transgenic mice were generated in CL56BL6 background as described previously . MutE1-4::GFP transgenes were prepared and linearized with the SalI/EcoRI for microinjection as described . Transient transgenic embryos were generated by pronuclear injection into fertilized eggs. Animal experiments were performed under the guidelines of the Korea Ministry of Food and Drug Safety (Act No. 9025) with approval of procedures by the Institutional Animal Care and Use Committee of Gwangju Institute of Science and Technology (GIST-2010-12).
293T cells were transiently transfected with reporters and transcription factors using Lipofectamine 2000 reagent (Invitrogen). CMV-β-galactosidase plasmid was co-transfected to normalize transfection efficiency. After 2 days, cell extracts were assayed for luciferase assays and β-galactosidase assays. Data represent as means of triplicate values, repeated at least three times.
Enhancer sequences from human, mouse, chicken, zebrafish, and fugu were retrieved from UCSC genome browser. Comparative analyses of the E1 and E2 sequences were done with mVISTA (genome.lbl.gov/vista) using the LAGAN alignment tool [74,75]. All E1 and E2 sequences from five species were searched against the lamprey genome using BLAST with “somewhat similar sequences (blastn)” options and an E-value cut-off of 0.05. Potential transcription factor binding sites in these conserved sequences were predicted using rVISTA . A phylogenetic tree was constructed with MEGA version 6 maximum likelihood method with 1,000 bootstrap replications .
Electrophoretic mobility shift assays
The sequences of sense and antisense oligonucleotide used in our EMSA were as follows:
mutE1-4- AS ′-GTCTGCCACTCCACTCGCGAACCCGCGTTGCATTTTATATTGG-3′
Biotin-labeled probe was incubated with nuclear extract of HEK 293T cells transfected with mouse Phox2a in binding buffer (10 mM Tris, pH 7.5, 50 mM KCl, 5 mM MgCl2, 1 mM dithiothreitol, 0.05% Nonidet P-40, and 2.5% glycerol) with poly(dI-dC) at RT. Competition reactions were performed by adding a 200-fold excess of unlabeled double-stranded probe. The reactions were resolved on non-denaturing 6% polyacrylamide gels and visualized by chemiluminescence (Thermo Scientific).
Embryonic mouse hindbrains were micro-dissected at E13.5. Cells were lysed after cross-linking with 1% formaldehyde. Chromatin with a DNA fragment length of less than 500 bp was obtained by sonication and immunoprecipitated with rabbit anti-Phox2b , rabbit anti-Histone H3 (Millipore) and rabbit IgG (Vector Labs). The Phox2 binding motif in E1 was amplified by PCR in each sample. For ChIP-Seq analysis, previous results were retrieved and examined around the Isl1 locus .
1. Jessell TM (2000) Neuronal specification in the spinal cord: inductive signals and transcriptional codes. Nature Reviews Genetics 1: 20–29. 11262869
2. Dubreuil V, Hirsch M, Pattyn A, Brunet J, Goridis C (2000) The Phox2b transcription factor coordinately regulates neuronal cell cycle exit and identity. Development 127: 5191–5201. 11060244
3. Pattyn A, Hirsch M, Goridis C, Brunet J-F (2000) Control of hindbrain motor neuron differentiation by the homeobox gene Phox2b. Development 127: 1349–1358. 10704382
4. Brunet J-F, Pattyn A (2002) Phox2 genes-from patterning to connectivity. Current Opinion in Genetics & Development 12: 435–440.
5. Thaler JP, Lee S-K, Jurata LW, Gill GN, Pfaff SL (2002) LIM factor Lhx3 contributes to the specification of motor neuron and interneuron identity through cell-type-specific protein-protein interactions. Cell 110: 237–249. 12150931
6. Dasen J, Jessell T (2009) Hox networks and the origins of motor neuron diversity. Current topics in developmental biology 88: 169–200. doi: 10.1016/S0070-2153(09)88006-X 19651305
7. Alexander T, Nolte C, Krumlauf R (2009) Hox genes and segmentation of the hindbrain and axial skeleton. Annual Review of Cell and Developmental 25: 431–456.
9. Pfaff SL, Mendelsohn M, Stewart CL, Edlund T, Jessell TM (1996) Requirement for LIM Homeobox Gene Isl1 in Motor Neuron Generation Reveals a Motor Neuron–Dependent Step in Interneuron Differentiation. Cell 84: 309–320. 8565076
10. Rousso D, Gaber Z, Wellik D, Morrisey E, Novitch B (2008) Coordinated actions of the forkhead protein Foxp1 and Hox proteins in the columnar organization of spinal motor neurons. Neuron 59: 226–240. doi: 10.1016/j.neuron.2008.06.025 18667151
11. Dasen JS, De Camilli A, Wang B, Tucker PW, Jessell TM (2008) Hox repertoires for motor neuron diversity and connectivity gated by a single accessory factor, FoxP1. Cell 134: 304–316. doi: 10.1016/j.cell.2008.06.019 18662545
12. Uemura O, Okada Y, Ando H, Guedj M, Higashijima S, et al. (2005) Comparative functional genomics revealed conservation and diversification of three enhancers of the isl1 gene for motor and sensory neuron-specific expression. Developmental biology 278: 587–606. 15680372
13. Kappen C, Salbaum JM (2009) Identification of regulatory elements in the Isl1 gene locus. The International journal of developmental biology 53: 935–946. doi: 10.1387/ijdb.082819ck 19598113
14. Dufour HD, Chettouh Z, Deyts C, de Rosa R, Goridis C, et al. (2006) Precraniate origin of cranial motoneurons. Proceedings of the National Academy of Sciences 103: 8727–8732.
15. Sugahara F, Aota S-i, Kuraku S, Murakami Y, Takio-Ogawa Y, et al. (2011) Involvement of Hedgehog and FGF signalling in the lamprey telencephalon: evolution of regionalization and dorsoventral patterning of the vertebrate forebrain. Development 138: 1217–1226. doi: 10.1242/dev.059360 21343370
16. Fetcho JR (1992) The spinal motor system in early vertebrates and some of its evolutionary changes. Brain, behavior and evolution 40: 82–97. 1422809
17. Kusakabe R, Kuratani S (2005) Evolution and developmental patterning of the vertebrate skeletal muscles: perspectives from the lamprey. Developmental dynamics 234: 824–834. 16252276
18. Freitas R, Zhang G, Cohn MJ (2006) Evidence that mechanisms of fin development evolved in the midline of early vertebrates. Nature 442: 1033–1037. 16878142
19. Funakoshi K, Nakano M (2007) The sympathetic nervous system of anamniotes. Brain, behavior and evolution 69: 105–113. 17230018
20. Osumi N, Hirota A, Ohuchi H, Nakafuku M, Iimura T, et al. (1997) Pax-6 is involved in the specification of hindbrain motor neuron subtype. Development 124: 2961–2972. 9247338
21. Li X, Zhao X, Fang Y, Jiang X, Duong T, et al. (1998) Generation of destabilized green fluorescent protein as a transcription reporter. Journal of Biological Chemistry 273: 34970–34975. 9857028
22. Tsuchida T, Ensini M, Morton S, Baldassare M, Edlund T, et al. (1994) Topographic organization of embryonic motor neurons defined by expression of LIM homeobox genes. Cell 79: 957–970. 7528105
23. Stein S, Fritsch R, Lemaire L, Kessel M (1996) Checklist: vertebrate homeobox genes. Mechanisms of development 55: 91–108. 8734502
24. Arber S, Han B, Mendelsohn M, Smith M, Jessell T, et al. (1999) Requirement for the homeobox gene Hb9 in the consolidation of motor neuron identity. Neuron 23: 659–674. 10482234
25. Dasen J, Jessell T (2008) Hox networks and the origins of motor neuron diversity. Current topics in developmental biology 88: 169–200.
26. Hoffmann S, Berger IM, Glaser A, Bacon C, Li L, et al. (2013) Islet1 is a direct transcriptional target of the homeodomain transcription factor Shox2 and rescues the Shox2-mediated bradycardia. Basic research in cardiology 108: 1–11.
27. Francius C, Clotman F (2010) Dynamic expression of the Onecut transcription factors HNF-6, OC-2 and OC-3 during spinal motor neuron development. Neuroscience 165: 116–129. doi: 10.1016/j.neuroscience.2009.09.076 19800948
28. Sander M, Paydar S, Ericson J, Briscoe J, Berber E, et al. (2000) Ventral neural patterning by Nkx homeobox genes: Nkx6. 1 controls somatic motor neuron and ventral interneuron fates. Genes & development 14: 2134–2139.
29. Ericson J, Rashbass P, Schedl A, Brenner-Morton S, Kawakami A, et al. (1997) Pax6 controls progenitor cell identity and neuronal fate in response to graded Shh signaling. Cell 90: 169–180. 9230312
30. Song M-R, Shirasaki R, Cai C-L, Ruiz EC, Evans SM, et al. (2006) T-Box transcription factor Tbx20 regulates a genetic program for cranial motor neuron cell body migration. Development 133: 4945–4955. 17119020
31. Mazzoni EO, Mahony S, Closser M, Morrison CA, Nedelec S, et al. (2013) Synergistic binding of transcription factors to cell-specific enhancers programs motor neuron identity. Nature neuroscience 16: 1219–1227. doi: 10.1038/nn.3467 23872598
32. Lee S, Lee B, Joshi K, Pfaff SL, Lee JW, et al. (2008) A regulatory network to segregate the identity of neuronal subtypes. Developmental cell 14: 877–889. doi: 10.1016/j.devcel.2008.03.021 18539116
33. Shirasaki R, Pfaff SL (2002) Transcriptional codes and the control of neuronal identity. Annual review of neuroscience 25: 251–281. 12052910
34. Hutchinson SA, Eisen JS (2006) Islet1 and Islet2 have equivalent abilities to promote motoneuron formation and to specify motoneuron subtype identity. Development 133: 2137–2147. 16672347
35. Song M-R, Sun Y, Bryson A, Gill GN, Evans SM, et al. (2009) Islet-to-LMO stoichiometries control the function of transcription complexes that specify motor neuron and V2a interneuron identity. Development 136: 2923–2932. doi: 10.1242/dev.037986 19666821
36. Novitch BG, Chen AI, Jessell TM (2001) Coordinate regulation of motor neuron subtype identity and pan-neuronal properties by the bHLH repressor Olig2. Neuron 31: 773–789. 11567616
37. Zhou Q, Anderson DJ (2002) The bHLH transcription factors OLIG2 and OLIG1 couple neuronal and glial subtype specification. Cell 109: 61–73. 11955447
38. Genethliou N, Panayiotou E, Panayi H, Orford M, Mean R, et al. (2009) SOX1 links the function of neural patterning and Notch signalling in the ventral spinal cord during the neuron-glial fate switch. Biochemical and biophysical research communications 390: 1114–1120. doi: 10.1016/j.bbrc.2009.08.154 19723505
39. Karunaratne A, Hargrave M, Poh A, Yamada T (2002) GATA proteins identify a novel ventral interneuron subclass in the developing chick spinal cord. Developmental biology 249: 30–43. 12217316
40. Kan L, Israsena N, Zhang Z, Hu M, Zhao L, et al. (2004) Sox1 acts through multiple independent pathways to promote neurogenesis. Developmental biology 269: 580–594. 15110721
41. Dorval KM, Bobechko BP, Fujieda H, Chen S, Zack DJ, et al. (2006) CHX10 targets a subset of photoreceptor genes. Journal of Biological Chemistry 281: 744–751. 16236706
42. Smith JJ, Kuraku S, Holt C, Sauka-Spengler T, Jiang N, et al. (2013) Sequencing of the sea lamprey (Petromyzon marinus) genome provides insights into vertebrate evolution. Nature genetics 45: 415–421. doi: 10.1038/ng.2568 23435085
43. McEwen GK, Goode DK, Parker HJ, Woolfe A, Callaway H, et al. (2009) Early evolution of conserved regulatory sequences associated with development in vertebrates. PLoS genetics 5: e1000762. doi: 10.1371/journal.pgen.1000762 20011110
44. Putnam NH, Butts T, Ferrier DE, Furlong RF, Hellsten U, et al. (2008) The amphioxus genome and the evolution of the chordate karyotype. Nature 453: 1064–1071. doi: 10.1038/nature06967 18563158
45. Roy A, Francius C, Rousso DL, Seuntjens E, Debruyn J, et al. (2012) Onecut transcription factors act upstream of Isl1 to regulate spinal motoneuron diversification. Development 139: 3109–3119. doi: 10.1242/dev.078501 22833130
46. Jung H, Lacombe J, Mazzoni EO, Liem KF Jr, Grinstein J, et al. (2010) Global Control of Motor Neuron Topography Mediated by the Repressive Actions of a Single Hox Gene. Neuron 67: 781–796. doi: 10.1016/j.neuron.2010.08.008 20826310
47. Vincentz JW, VanDusen NJ, Fleming AB, Rubart M, Firulli BA, et al. (2012) A Phox2-and Hand2-dependent Hand1cis-regulatory element reveals a unique gene dosage requirement for Hand2 during sympathetic neurogenesis. The Journal of Neuroscience 32: 2110–2120. doi: 10.1523/JNEUROSCI.3584-11.2012 22323723
48. Lee S-K, Lee B, Ruiz EC, Pfaff SL (2005) Olig2 and Ngn2 function in opposition to modulate gene expression in motor neuron progenitor cells. Genes & Development 19: 282–294.
49. Gadd MS, Bhati M, Jeffries CM, Langley DB, Trewhella J, et al. (2011) Structural basis for partial redundancy in a class of transcription factors, the LIM homeodomain proteins, in neural cell type specification. Journal of Biological Chemistry 286: 42971–42980. doi: 10.1074/jbc.M111.248559 22025611
50. Gadd MS, Jacques DA, Nisevic I, Craig VJ, Kwan AH, et al. (2013) A structural basis for the regulation of the LIM-homeodomain protein Islet 1 (Isl1) by intra-and intermolecular interactions. Journal of Biological Chemistry 288: 21924–21935. doi: 10.1074/jbc.M113.478586 23750000
51. Sanchez-Garcia I, Osada H, Forster A, Rabbitts T (1993) The cysteine-rich LIM domains inhibit DNA binding by the associated homeodomain in Isl-1. The EMBO journal 12: 4243–4250. 7901000
52. Liu J, Hunter CS, Du A, Ediger B, Walp E, et al. (2011) Islet-1 Regulates Arx Transcription during Pancreatic Islet α-Cell Development. The Journal of Biological Chemistry 286: 15352–15360. doi: 10.1074/jbc.M111.231670 21388963
53. Zhang H, Wang W, Guo T, Yang J, Chen P, et al. (2009) The LIM-homeodomain protein ISL1 activates insulin gene promoter directly through synergy with BETA2. Journal of molecular biology 392: 566–577. doi: 10.1016/j.jmb.2009.07.036 19619559
54. Nasif S, de Souza F, González L, Yamashita M, Orquera D, et al. (2015) Islet 1 specifies the identity of hypothalamic melanocortin neurons and is critical for normal food intake and adiposity in adulthood. Proceedings of the National Academy of Sciences of the United States of America 112: E1861–1870. doi: 10.1073/pnas.1500672112 25825735
55. Giuliano P, Marino R, Pinto MR, De Santis R (1998) Identification and developmental expression of Ci-isl, a homologue of vertebrate islet genes, in the ascidian Ciona intestinalis. Mechanisms of development 78: 199–202. 9858732
56. Fetcho JR (1987) A review of the organization and evolution of motoneurons innervating the axial musculature of vertebrates. Brain Research Reviews 12: 243–280.
57. Fetcho J, Reich N (1992) Axial motor organization in postmetamorphic tiger salamanders (Ambystoma tigrinum): a segregation of epaxial and hypaxial motor pools is not necessarily associated with terrestrial locomotion. Brain, behavior and evolution 39: 219–228. 1633553
58. Jackman W, Langeland J, Kimmel C (2000) islet Reveals Segmentation in the Amphioxus Hindbrain Homolog. Developmental Biology 220: 16–26. 10720427
59. Murakami Y, Pasqualetti M, Takio Y, Hirano S, Rijli FM, et al. (2004) Segmental development of reticulospinal and branchiomotor neurons in lamprey: insights into the evolution of the vertebrate hindbrain. Development 131: 983–995. 14973269
60. Parker HJ, Bronner ME, Krumlauf R (2014) A Hox regulatory network of hindbrain segmentation is conserved to the base of vertebrates. Nature 514: 490–493. doi: 10.1038/nature13723 25219855
61. Holland LZ, Holland ND (1999) Chordate origins of the vertebrate central nervous system. Current opinion in neurobiology 9: 596–602. 10508734
62. Odom DT, Dowell RD, Jacobsen ES, Gordon W, Danford TW, et al. (2007) Tissue-specific transcriptional regulation has diverged significantly between human and mouse. Nature genetics 39: 730–732. 17529977
63. Ritter DI, Li Q, Kostka D, Pollard KS, Guo S, et al. (2010) The importance of being cis: evolution of orthologous fish and mammalian enhancer activity. Molecular biology and evolution 27: 2322–2332. doi: 10.1093/molbev/msq128 20494938
64. Jung H, Mazzoni EO, Soshnikova N, Hanley O, Venkatesh B, et al. (2014) Evolving Hox Activity Profiles Govern Diversity in Locomotor Systems. Developmental Cell 29: 171–187. doi: 10.1016/j.devcel.2014.03.008 24746670
65. Murakami Y, Tanaka M (2011) Evolution of motor innervation to vertebrate fins and limbs. Developmental biology 355: 164–172. doi: 10.1016/j.ydbio.2011.04.009 21540022
66. Eberhart J, Swartz M, Koblar SA, Pasquale E, Krull CE (2002) EphA4 constitutes a population-specific guidance cue for motor neurons. Developmental biology 247: 89–101. 12074554
67. Kania A, Jessell TM (2003) Topographic motor projections in the limb imposed by LIM homeodomain protein regulation of ephrin-A: EphA interactions. Neuron 38: 581–596. 12765610
68. Mehta TK, Ravi V, Yamasaki S, Lee AP, Lian MM, et al. (2013) Evidence for at least six Hox clusters in the Japanese lamprey (Lethenteron japonicum). Proceedings of the National Academy of Sciences 110: 16044–16049.
69. Lee S-K, Jurata LW, Funahashi J, Ruiz EC, Pfaff SL (2004) Analysis of embryonic motoneuron gene regulation: derepression of general activators function in concert with enhancer factors. Development 131: 3295–3306. 15201216
70. Yasugi S, Nakamura H (2000) Gene transfer into chicken embryos as an effective system of analysis in developmental biology. Development, growth & differentiation 42: 195–197.
71. Thaler JP, Koo SJ, Kania A, Lettieri K, Andrews S, et al. (2004) A postmitotic role for Isl-class LIM homeodomain proteins in the assignment of visceral spinal motor neuron identity. Neuron 41: 337–350. 14766174
72. Shirasaki R, Lewcock JW, Lettieri K, Pfaff SL (2006) FGF as a target-derived chemoattractant for developing motor axons genetically programmed by the LIM code. Neuron 50: 841–853. 16772167
73. Lee B, Rizzoti K, Kwon DS, Kim S-Y, Oh S, et al. (2012) Direct transcriptional regulation of Six6 is controlled by SoxB1 binding to a remote forebrain enhancer. Developmental biology 366: 393–403. doi: 10.1016/j.ydbio.2012.04.023 22561201
74. Brudno M, Do CB, Cooper GM, Kim MF, Davydov E, et al. (2003) LAGAN and Multi-LAGAN: efficient tools for large-scale multiple alignment of genomic DNA. Genome research 13: 721–731. 12654723
75. Frazer KA, Pachter L, Poliakov A, Rubin EM, Dubchak I (2004) VISTA: computational tools for comparative genomics. Nucleic acids research 32: W273–W279. 15215394
76. Loots GG, Ovcharenko I (2004) rVISTA 2.0: evolutionary analysis of transcription factor binding sites. Nucleic acids research 32: W217–W221. 15215384
77. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S (2013) MEGA6: molecular evolutionary genetics analysis version 6.0. Molecular biology and evolution 30: 2725–2729. doi: 10.1093/molbev/mst197 24132122
78. Pattyn A, Morin X, Cremer H, Goridis C, Brunet J-F (1997) Expression and interactions of the two closely related homeobox genes Phox2a and Phox2b during neurogenesis. Development 124: 4065–4075. 9374403