Identification of Phox2b-regulated genes by expression profiling of cranial motoneuron precursors
© Pla et al.; licensee BioMed Central Ltd. 2008
Received: 27 November 2007
Accepted: 19 June 2008
Published: 19 June 2008
Skip to main content
© Pla et al.; licensee BioMed Central Ltd. 2008
Received: 27 November 2007
Accepted: 19 June 2008
Published: 19 June 2008
Branchiomotor neurons comprise an important class of cranial motor neurons that innervate the branchial-arch-derived muscles of the face, jaw and neck. They arise in the ventralmost progenitor domain of the rhombencephalon characterized by expression of the homeodomain transcription factors Nkx2.2 and Phox2b. Phox2b in particular plays a key role in the specification of branchiomotor neurons. In its absence, generic neuronal differentiation is defective in the progenitor domain and no branchiomotor neurons are produced. Conversely, ectopic expression of Phox2b in spinal regions of the neural tube promotes cell cycle exit and neuronal differentiation and, at the same time, induces genes and an axonal phenotype characteristic for branchiomotor neurons. How Phox2b exerts its pleiotropic functions, both as a proneural gene and a neuronal subtype determinant, has remained unknown.
To gain further insights into the genetic program downstream of Phox2b, we searched for novel Phox2b-regulated genes by cDNA microarray analysis of facial branchiomotor neuron precursors from heterozygous and homozygous Phox2b mutant embryos. We selected for functional studies the genes encoding the axonal growth promoter Gap43, the Wnt antagonist Sfrp1 and the transcriptional regulator Sox13, which were not previously suspected to play roles downstream of Phox2b and whose expression was affected by Phox2b misexpression in the spinal cord. While Gap43 did not produce an obvious phenotype when overexpressed in the neural tube, Sfrp1 induced the interneuron marker Lhx1,5 and Sox13 inhibited neuronal differentiation. We then tested whether Sfrp1 and Sox13, which are down-regulated by Phox2b in the facial neuron precursors, would antagonize some aspects of Phox2b activity. Co-expression of Sfrp1 prevented Phox2b from repressing Lhx1,5 and alleviated the commissural axonal phenotype. When expressed together with Sox13, Phox2b was still able to promote cell cycle exit and neuronal differentiation, but the cells failed to relocate to the mantle layer and to extinguish the neural stem cell marker Sox2.
Our results suggest novel roles for Sfrp1 and Sox13 in neuronal subtype specification and generic neuronal differentiation, respectively, and indicate that down-regulation of Sfrp1 and Sox13 are essential aspects of the genetic program controlled by Phox2b in cranial motoneurons.
Branchiomotor (bm) neurons comprise an important class of cranial motoneurons. They arise in the hindbrain and rostral cervical spinal cord from the ventralmost progenitor domain characterized by expression of the transcription factors Nkx2.2, Mash1 and Phox2b [1, 2]. Targeted mutation in the mouse has shown that their development crucially depends on Phox2b whereas the inactivation of Nkx2.2 or Mash1 results in only mild defects of bm development [2–5]. In the absence of Phox2b, the differentiation of all bm neurons is arrested at an early stage and their precursors either die by apoptosis or switch fate and become serotonergic neurons. Conversely, Phox2b induces neurons with bm characteristics when misexpressed in the neural tube of chicken embryos [2, 6]. Hence, Phox2b appears both necessary and sufficient for the implementation of a specific cranial motoneuronal phenotype. Loss and gain of function experiments also show that Phox2b has proneural activity  and promotes cell cycle exit and generic neuronal differentiation  (see for review ). The development of the facial bm (fbm) neurons that arise in rhombomere 4 (r4) has been particularly well studied in the Phox2b knockout embryos. A number of downstream target genes have been identified that depend on Phox2b for their expression in fbm neurons by examining expression of candidate genes and by a subtractive screening approach. They include generic markers for young post-mitotic neurons (Tubb3, Nfl, Ebf2, Math3) and genes more specifically associated with bm neuron differentiation (Phox2a, Islet1, Rgs4, Tbx20) ([4, 10] and unpublished results).
Despite this wealth of information, only part of the genetic program that depends on Phox2b in fbm neurons has been identified. Indeed, the subtractive screening approach proved far from exhaustive since most previously identified downstream genes were not detected and none were found to be up-regulated in the mutants. Here we have used microarrays to compare the patterns of gene expression in fbm precursors from homozygous (Phox2bLacZ/LacZ) mutants, in which fbm differentiation is arrested, with that of heterozygous Phox2b (Phox2bLacZ/+) embryos, in which fbm neurons are phenotypically normal. In addition to genes already identified in previous work, a number of new candidates turned up in the screen, among which many were up-regulated in the mutants and might thus be repressed by Phox2b. Most of the down-regulated genes were panneuronally expressed in young post-mitotic neurons. They were thus expected to pop up in the screen, since fewer post-mitotic neurons are produced in the mutants, and were not analyzed further. We selected a few other genes for further study by gain of function experiments in chicken embryos: Gap43, whose expression, although panneuronal, was strikingly enhanced in post-mitotic fbm precursors; and two that were up-regulated in the mutants, Sfrp1 and Sox13. Overexpression of Sox13, whose function in neural development has not been studied so far, inhibited neuronal differentiation and the relocation to the mantle layer (ML) caused by Phox2b misexpression in the embryonic spinal cord, while that of Sfrp1 prevented the repression of the interneuronal marker Lhx1,5.
Phox2b being an important regulator of the development of bm neurons, we wanted to further investigate its function in the specification of bm neurons by searching for novel genes regulated by Phox2b in fbm neuron precursors. We thus chose to compare gene expression in ventral r4 from heterozygous (Phox2bLacZ/+) and homozygous (Phox2bLacZ/LacZ) embryonic day 10.5 (E10.5) mutant embryos. At this stage, fbm differentiation is well underway in Phox2b+/+ and Phox2bLacZ/+animals while no signs of bm differentiation are seen in the homozygous mutants . Four pairs of Phox2bLacZ/+and Phox2bLacZ/LacZRNA samples were amplified from microdissected ventral r4, labeled with Cy3 or Cy5, mixed and hybridized to two mouse cDNA microarrays, a microarray of the 15 K NIA probe set representing approximately 12,000 different genes and a home-made NeuroDev microarray representing 2,000 genes chosen for their expression in developing nervous tissue. We adopted as criteria for putative Phox2b-regulated genes a significance analysis of microarray (SAM) score (see Materials and methods) below -3.5 or above 3.5 for genes up- or down-regulated in the mutants, respectively, and a false discovery rate of 5%. By these criteria, 51 genes were found to be down-regulated and 23 up-regulated in the homozygous mutants (Additional files 1 and 2). Surprisingly, only three genes showed a greater than twofold change in expression. However, out of the 12 genes we selected for validation by in situ hybridization (ISH), 10 were confirmed as being clearly differentially expressed, despite fold changes as low as 1.14 and 1.13 for two of them (see below). Part of the rather low signal ratios on microarrays may be due to the fact that the floor plate, which should not be affected in the mutants, was always included in the material used for RNA preparation.
We found Sfrp1 and Tcf3 to be normally expressed in and confined to the VZ, Sfrp1 only in the ventral part of the basal plate, and Tcf3 throughout the neural tube. In the mutant territory, their expression extended up to the neural tube border, thus including the post-mitotic neurons in the ML (Figure 2e–h). A similar expansion into the mutant ML has been found in previous work for Dll1 and Nkx2.2, which are normally confined to the VZ . In contrast to Tcf3, Sfrp1 appeared also to be up-regulated throughout the mutant VZ: the signal given by the Sfrp1 probe was stronger outside of the fbm progenitors in Phox2bLacZ/+embryos, but the reverse was true in Phox2bLacZ/LacZanimals. Tcfs are transcription factors that function in the Wnt signaling pathway [18, 19]. Secreted frizzled-related proteins (Sfrps) also function mainly in the Wnt pathway as antagonists of Wnt signaling, but Sfrp1 may also exert a Wnt-independent activity as a stimulator of neurite outgrowth and as a growth cone attractant for retinal axons [20–22].
In Hamburger and Hamilton stage 18 (HH18) chicken hindbrain, we also found Sox13 expression very weak or absent in the Phox2b+ domain that gives rise to bm/visceromotor neurons, while it extended into the ventral progenitor domain at HH24, when Phox2b has disappeared there. At this stage, a distinct Phox2b+ progenitor domain has formed in the basal alar plate. Strikingly, Sox13 expression in the VZ was interrupted at this location (Figure 3g–j). Together, these results suggest that Phox2b represses Sox13 expression in the hindbrain VZ. The effect may be indirect since Sox13 continues to be down-regulated in the progenitors when Mash1 is absent.
We then selected Gap43, Sfrp1 and Sox13 for gain of function experiments in the chicken embryo spinal cord to ask whether they were regulated by ectopically expressed Phox2b. We chose Gap43 as a panneuronal Phox2b-dependent gene because it was more strongly expressed in fbm neurons, Sfrp1 as an example of the VZ genes that expand into the mutant ML since it was at the same time up-regulated in the VZ, and Sox13 because Phox2b repressed it in bm progenitors.
Overexpression of Sfrp1 in the chicken embryo spinal cord also did not produce gross anatomical alterations in the neural tube detectable by green fluorescent protein (GFP) labeling (Figure 6a, compare with Figure 5a). Since Sfrp1 appeared to be repressed by Phox2b in fbm precursors and when misexpressed in the spinal cord, we asked whether Sfrp1 would counteract some of the effects produced by Phox2b misexpression. After co-electroporation of Sfrp1 together with Phox2b, virtually all transfected cells expressed both GFP from the Sfrp1 expression vector and Phox2b, showing that co-expression is highly efficient (Figure 6d). When expressed together with Phox2b in this way, Sfrp1 did not prevent Phox2b from either driving neuronal differentiation and relocation to the ML or promoting ectopic neural tube exit of axons (Figure 6c, d; data not shown).
Sfrp1 (SARP2) has been found to sensitize cells in culture to pro-apoptotic stimuli, whereas Sfrp2 had the opposite effect  and prevented programmed cell death in the neural crest . We thus tested whether Sfrp1 might influence cell death in our overexpression paradigm. In line with published data , we found very few apoptotic cells in the neural tube at the stages studied (one to four TUNEL+ cells per section), and their number differed neither between electroporated and non-electroporated sides of the neural tube nor among different conditions (electroporation of Srp1 or Phox2b or of Sfrp1 plus Phox2b) (not shown).
We next tested the idea that Sox13 overexpression may antagonize Phox2b in its capacity to drive neuronal differentiation. Phox2b misexpressed alone causes the transfected cells to accumulate in the ML and to express NeuN (Figure 8a), as expected from its propensity to promote neuronal differentiation and relocation to the ML. When Sox13 was misexpressed together with Phox2b, many transfected cells remained in the VZ where they either retained the morphology of neuroepithelial cells or aggregated near the lumen of the neural tube (Figure 8c, f, l; n = 3/3). Phox2b was still able to induce NeuN, but now many cells double-positive for NeuN and GFP were located in the VZ and had the morphology of neuroepithelial cells. Hence, Sox13, although unable to inhibit induction of NeuN by Phox2b, appears to prevent the transfected cells from losing their neuroepithelial morphology and from relocating to the ML. We then investigated further the phenotype of the cells constitutively co-expressing Phox2b and Sox13. We first established that, in the presence of Sox13, Phox2b was still able to promote cell cycle exit as shown by the fact that the transfected cells did not incorporate BrdU (Figure 8f; n = 3/3). However, while the neural stem cell marker Sox2  was never expressed in cells transfected with Phox2b alone, many doubly transfected cells still expressed Sox2 (Figure 8g, i; n = 3/3). Finally, we examined Islet1,2 as a marker for the bm phenotype induced by Phox2b, and Lhx1,5 as a marker for interneurons. Induction of Islet1,2 was completely repressed by co-transfection with Sox13 (Figure 8j, l; n = 3/3). By contrast, Sox13 did not prevent down-regulation of Lhx1,5 after Phox2b transfection, a result to be expected since Sox13 overexpression inhibits neuronal differentiation (Figure 8m; n = 3/3).
Phox2b plays a key role in the formation of two important classes of cranial motoneurons, the bm and the closely related visceromotor neurons. Its role in motoneuronal differentiation has been studied in most detail during differentiation of the fbm neurons. This work has yielded a few genetic targets whose expression is dependent on Phox2b in fbm precursors [4, 10]. To search for Phox2b-regulated genes by a more systematic approach, we have used gene expression profiling of ventral r4, in which the fbm precursors arise, from heterozygous and homozygous Phox2b mutant embryos. In this way we have identified a number of candidates that are either up- or down-regulated in the homozygous mutants. Two of these candidates regulated by Phox2b had been identified in previous work, but none of the others have been reported before. We restrict our discussion to those candidates we confirmed by ISH as being regulated by Phox2b in ventral r4.
Our analysis did not distinguish between genes that were under direct control of Phox2b and genes that depended on Phox2b indirectly, through the action of other transcription factors. However, the genes that we found induced or repressed in the chicken spinal cord already at 24 hae can be supposed to be proximal targets in the genetic program induced by Phox2b, and not to be dependent on a vast network of other factors. So far, only transcriptional activation by Phox2b has been demonstrated, but the possibility exists that it may also act as a repressor, as indicated by the observation that repressor forms of Phox2b reproduced some of the effects of Phox2b .
Since less post-mitotic neurons are produced in the mutant territory, it was to be expected that many genes that are normally expressed in young post-mitotic neurons would be underrepresented in Phox2bLacZ/LacZembryos. App, Mapk8, Nfl and Gap43 belong to this category. Although post-mitotic cells are produced by Phox2bLacZ/LacZprogenitors and migrate to the ML, albeit in reduced numbers, the ML cells failed to express these genes. The same pattern has been found in previous work for other post-mitotic markers . Some of the post-mitotic cells that arise in the mutants then go on to become serotonergic neurons [1, 7]. Apparently, in the absence of Phox2b and its proneural activity, these cells will acquire neuronal properties only at later stages. Gap43 is a case apart since its expression in post-mitotic neurons was distinctly stronger in fbm precursors than elsewhere in r4, suggesting that high Gap43 activity may be required for proper fbm differentiation. Gap43 was also markedly up-regulated by ectopic Phox2b expression in the spinal cord, an effect that could not be ascribed to an increase in generic neuronal differentiation since Ngn2 overexpression was without effect. Gap43 was originally described as a protein associated with growing axons, both in development and during regeneration, and is thought to function as a determinant of axonal growth competence [11, 12, 34]. Despite a wealth of in vitro data describing its growth-promoting properties, its in vivo function is still ill-defined. When constitutively expressed in adult mice, GAP43 induced branching of motor nerve terminals and promoted sprouting of nerves after lesion . In our study, Gap43 overexpression in the embryonic neural tube did not produce overt morphological alterations. Hence, its effects might be limited to terminal fields and situations of regeneration. A number of discrete developmental defects have been reported in Gap43 mutants, including pathfinding defects of the optic nerve , a reduction in serotonergic innervation of the forebrain  and a failure to form the anterior commissure . Such a subtle phenotype for a gene expressed in virtually all young neurons may be due to redundancy with the CAP23 gene .
Among the genes found to be up-regulated in mutant ventral r4, Hes1, at the stage studied, was expressed only in the dorsalmost part of the Phox2b+ VZ in wild-type or heterozygous animals. In Phox2bLacZ/LacZembryos, its expression invaded most of the normally Phox2b+ domain. The basic-helix-loop-helix transcription factor of the hairy and enhancer of split family, Hes1 is best known as an effector of Notch signaling and, together with its paralogues Hes3 and Hes5, plays a crucial role in maintaining the progenitor state of neuroepithelial cells [16, 17]. Since Phox2b drives neural progenitors to exit the cell cycle and to differentiate, it is perhaps not surprising that it shuts off a negative regulator of neural differentiation. This scheme, however, is complicated by the observation that its paralogue Hes5 has been found to be down-regulated in Phox2bLacZ/LacZembryos . It appears that the lack of proneural activity in Phox2b mutants and the consecutive down-regulation of the Notch ligand Dll1 in the VZ affect differently the two paralogous genes. The Lef/Tcf proteins are transcription factors that play pivotal roles in Wnt signaling, functioning mainly as repressors in the absence and as obligatory co-activators of beta-catenin in the presence of Wnt ligands . Tcf3 has remained the most enigmatic member of the family. Loss of Tcf3 results in early gastrulation defects and perhaps as a consequence to broad defects in neural tube patterning . Tcf3 has been said to become undetectable in the mouse embryo from E10.5 onward , but we found it strongly expressed in E10.5 r4 throughout the VZ. Normally excluded from the ML, its expression domain invaded the ML in the mutant territory. This has been observed in previous work also for Nkx2.2 and Dll1 . Apparently, the post-mitotic cells that accumulate in the mutant ML are not only unable to express generic neuronal markers but also do not switch off genes whose expression is normally confined to the VZ.
Sfrp1 shared with Tcf3 the tendency to be expressed in the mutant ML, from which it was normally excluded. However, the expression pattern of Sfrp1 in r4 was much more circumscribed than that of Tcf3. It was only expressed in ventral r4, including the Phox2b+ progenitor domain and extending up to the middle of the basal plate. In Phox2bLacZ/+embryos, its expression was stronger dorsal to the Phox2b+ progenitors, but the reverse was true in Phox2bLacZ/LacZembryos. In the embryonic chicken spinal cord, we found Sfrp1 to be strongly expressed in the VZ of the dorsal basal plate, in line with previous results , and to be massively repressed by Phox2b. Sfrp1 overexpression in the chicken spinal cord up-regulated the interneuron marker Lhx1,5. When co-expressed with Sfrp1, Phox2b was still capable of inducing Islet1,2, but became unable to repress expression of Lhx1,5. Hence, co-expression of Phox2b and Sfrp1 results in a mixed phenotype that is not normally encountered in the spinal cord, as also indicated by the observation that Sfrp1 partially restored the ability of Phox2b-transfected neurons to grow commissural axons. However, expression of Pax2, another interneuronal marker that was down-regulated by Phox2b, was already reduced by expression of Sfrp1 alone and continued to be down-regulated when Sfrp1 was co-expressed with Phox2b.
The Sfrps have been mostly described as secreted inhibitors of Wnt signaling, but there is evidence that they are able to promote rather than to inhibit Wnt signaling in some contexts and that they can act through Wnt-independent mechanisms [20–22]. In the Wnt1-/-; Wnt3a-/- embryonic spinal cord, the dI4 (Lhx1,5+; Pax2+) subclass of interneurons has been found to be expanded at the expense of the dI1 and dI2 interneuron populations . In line with these results, we found that Sfrp1 increased the population of Lhx1,5-expressing cells. Sfrp1, however, down-regulated Pax2, suggesting that the effect of Sfrp1 overexpression cannot be explained solely by inhibition of Wnt signaling. Gain of function experiments in the chicken retina have shown that Sfrp1 promotes neuronal differentiation and the generation of retinal ganglion cells and cones at the expense of other cell types by a mechanism that does not involve inhibition of Wnt signaling . In our study, Sfrp1 overexpression in the spinal cord did not cause an obvious increase in generic neuronal differentiation. Still, our results suggest that Sfrp1 plays a role in neuronal subtype specification as has been observed in the retina.
Transcription factors of the B-, C- and E-group Sox proteins have been found to play key roles during early neural development [33, 45–47]. By contrast, much less is known about the expression of the D-group Sox proteins in neural tissues and virtually nothing is known about their function during central nervous system development. In our screen, the D-group Sox13 gene was the only Sox gene that scored as a Phox2b-regulated gene. A previous study reported Sox13 to be expressed mainly by differentiating neurons in the ML . At the slightly earlier stages and more caudal regions studied here, we found, by both ISH and antibody staining, Sox13 to be expressed throughout and confined to the VZ in both the mouse and chicken nervous systems. We found Sox13 to be expressed in the VZ in a salt-and-pepper pattern, patches of high signal intensity alternating with cells expressing lower levels, suggesting a highly dynamic expression pattern perhaps in relation to the commitment to differentiation. VZ cells that expressed constitutively high levels of Sox13 after transfection were proliferating as indicated by BrdU incorporation and many expressed the neural cell stem marker Sox2, but none the neuronal marker NeuN, whose expression in the ML was inhibited. These results may be interpreted to mean that the cells normally expressing lower levels were those destined to become neurons. In E10.5 ventral r4, Sox13 mRNA and protein expression, which was strong in more dorsal regions, was much weaker in the fbm progenitor domain, but invaded the normally Phoxb+ region in the mutants. While Phox2b expression in the VZ and motor neuron production is prolonged in ventral r4, peak motor neuron production in r2 occurs at E9.5 and ceases around E10.5 together with Phox2b expression in the VZ, where the progenitors switch to the production of serotonergic neurons . Consequently, in E9.5 ventral r2, Sox13 expression was absent from ventral progenitors, but became detectable there one day later. At this stage, the basic-helix-loop-helix factor Mash1 has been found to be required in this domain as both the sole provider of proneural activity and a determinant of the serotonergic phenotype . As shown here, Sox13 expression did not expand into the ventral progenitor domain when Mash1 was absent, although Phox2b expression had ceased. One possibility is that Sox13 expression requires activation of the Notch pathway by its ligands, whose expression depends on proneural activity provided by Mash1, but more work is required to ascertain this. The idea that Sox13 may function as a Notch effector is also supported by the observation that its overexpression inhibits neurogenesis.
In contrast to what might have been expected from the pronounced down-regulation of Sox13 in Phox2b-expressing progenitors, Phox2b misexpression in spinal regions of the chicken neural tube did not result in downright repression, but areas of down-regulation alternated with regions with enhanced expression. Again, one reason may be that Sox13 expression depends on Notch activity, which may be repressed by Phox2b but enhanced in neighboring cells by the propensity of Phox2b to induce the Notch ligand Dll1 . Constitutively high expression of Sox13, however, counteracted Phox2b activity in our ectopic expression paradigm. A major effect was to block the relocation of the transfected cells to the ML. In the presence of forced expression of Sox13, Phox2b was still able to promote cell cycle exit and neuronal differentiation, but now the Phox2b+, NeuN+ and BrdU- cells normally located in the ML accumulated in the VZ. Many of the doubly transfected cells expressed Sox2, which was totally absent from cells transfected with Phox2b alone, and Islet1,2 induction by Phox2b was inhibited in them. Hence, Sox13 appears to counteract the capacity of Phox2b to silence Sox2 expression and to induce a bm phenotype. However, blocking migration to the ML may still be the primary effect, since persistent Sox2 expression and inhibition of Islet1,2 induction may be consequences of the former. So far, the only in vivo function of Sox13 that has been reported is the promotion of, and a requirement for, proper development of the gamma/delta T cell lineage . This has been attributed to the inhibition of Wnt signaling by direct interaction with the Wnt effector TCF1. Given that this also holds true for Sox13 function in the nervous system, it is interesting that two inhibitors of Wnt signaling, Sfrp1 and Sox13, appeared in our screen. However, Sfrp1 and Sox13 act on different processes: Sox13 inhibits generic neuronal differentiation while Sfrp1 affects cell type specification. Repression of both may thus be a prerequisite for proper development of bm neurons. Sox13 knockout mice die at birth from unknown reasons , and this may well be because of defective neural development. Our results showing that Sox13 inhibits neurogenesis and blocks migration to the ML give a first glance of the possible in vivo functions of Sox13 in neural development, opening the way for future mechanistic studies.
Phox2bLacZ/+ mice  were intercrossed and the hindbrains of E10.5 embryos were dissected out and incubated for 2 minutes at 37°C with 1 μg/μl fluorescein di-(beta-D-galactopyranoside) (Sigma-Aldrich, Lyon, France). The left and right beta-galactosidase expressing domains of ventral r4 with the floorplate in between were excised under a fluorescent binocular and stored at -20°C in RNALater solution (Ambion, Austin, TX, USA). The rest of the embryo was used for genotyping as described by Pattyn et al. . Mash1+/- mice (obtained from François Guillemot) were intercrossed and the embryos genotyped as described .
Four r4 samples of the same genotype (either Phox2bLacZ/+or Phox2bLacZ/LacZ) were pooled together and RNA was extracted and amplified using the RNAqueous-Micro Kit (Ambion) and the Message Amp aRNA amplification kit (Ambion) according to the instructions of the manufacturer. The quantity and the quality of the amplified RNA were verified using Agilent RNA 6000 Nano Assay and Bioanalyzer 2100 (Agilent Technologies, Massy, France).
Microarray hybridization was done using two cDNA collections, the NIA Mouse 15 k clone sets  and the NeuroDev clone sets (Ecole normale supérieure local collection), both spotted on Ultra GAPS II slides (Corning, Avon, France) using a BioRobotics MicroGrid microarrayer (Genomic Solutions, Huntingdon, UK) by the Transcriptome platform of ENS, Paris. Amplified RNA (2 μg) was retrotranscribed and labeled with dUTP-Cy5 or dUTP-Cy3 (Amersham, Piscataway, NJ, USA) using Superscript III Reverse Transcriptase (Invitrogen, Carlsbad, CA, USA) in the presence of 5 μg of random primers and appropriate buffer (Invitrogen). RNA was degraded by alkaline lysis and cDNA was purified using a QIAquick purification kit (Qiagen, Courtaboeuf, France). CDNA microarrays were pre-hybridized for 1 h at 42°C in 5 × SSC, 0.1% SDS, 1% bovine serum albumin, washed with distilled water and dried. They were hybridized overnight at 42°C in 25% formamide, 5 × SSC, 0.1% SDS with Cy5- and Cy3-labeled cDNA, washed in 1 × SSC, 0.2% SDS at 42°C, in 0.1 × SSC, 0.2% SDS and two times in 0.1 × SSC at room temperature and then spin-dried. Four biological replicas (independent dissections) for each microarray (15 k and NeuroDev) were done and each of these replicas was hybridized independently twice by inverting the labeling dyes (dye-swap). A total of 16 hybridizations where thus performed, 8 using NIA 15 k microarrays and 8 others using NeuroDev slides.
The hybridized microarrays were scanned using Genepix 4000B (Molecular Devices, Sunnyvale, CA, USA)) and the resulting image files analyzed by GenePix Pro 5.0 software (Axon). For each GenePix output file, two filters were applied, one to clear out spots and another one to discard saturating spots where the median foreground intensity was greater than 60,000 in one of the two channels. The resulting median foreground intensities were normalized, without background signal subtraction, for dye bias using a global Lowess correction . An expression matrix for each type of microarray (eight experiments each) was created gathering all normalized ratio data (M values). After normalization, we applied on each matrix a hierarchical clustering algorithm  using TMeV software  to detect a possible bias introduced by sample heterogeneity. The cluster obtained showed a clear separation according to dye-swap between hybridizations indicative of homogenous samples. The two matrices obtained were processed using the GEPAS web server  to remove genes with more than 30% missing values. A KNN impute step was performed to calculate missing values for the remaining genes using the 15 nearest neighbor profiles. A further scaling step was performed by dividing the ratios for each gene in one experiment by the experimental standard deviation and multiplying each value by the mean of all experimental standard deviations. The ratios were displayed such that a positive ratio describes a gene whose expression is higher in Phox2bLacZ/+than in Phox2bLacZ/LacZ embryos. SAM  was used to search for differentially expressed genes and to take advantage of our four biological replicates to calculate the SAM score equivalent to the T-statistic value and to estimate the associated false discovery rate. We fixed the false discovery rate at 5% using the 90th percentile of falsely called genes. We chose in addition a SAM score above 3.5 or below -3.5 (for genes down- or up-regulated in the mutants, respectively) as an arbitrary cut-off since smaller scores increased the number of genes whose expression changes could not be confirmed by ISH.
We used pCAGGS::Phox2b-IRES2-EGFP , pCAGGS::Gap43-IRES2-EGFP, pCAGGS::Sox13-IRES2-EGFP, pCAGGS::Sfrp1-IRES2-EGFP or pCAGGS::IRES2-EGFP to express mouse Phox2b, mouse Gap43, human SFRP1 (the latter two subcloned from IMAGE clones) and human Sox13 (subcloned from Ultimate ORF Clone IOH44934, Invitrogen) together with GFP or GFP alone in the chicken neural tube driven by a composite chicken beta actin-CMV promoter. The spinal region of HH12-14 embryonic chick neural tubes were electroporated in ovo with the appropriate expression vectors (2 mg/ml or 1 mg/ml) essentially as described . Correct expression of all constructs was verified by ISH with the cognate probes. After electroporation, embryos were allowed to develop at 38°C for 24–72 h.
Mouse embryos and well-electroporated chicken embryos, as assessed by GFP fluorescence, were dissected out, fixed in 4% paraformaldehyde overnight, embedded in Tissue-Tek® OCT™ (Sakura-Finetek, Zoeterwoude, The Netherlands) (mouse embryos) or gelatin (chicken embryos) and analyzed on 12 or 20 μm transverse cryosections, respectively. The methods used for ISH immunohistochemistry were as described previously [2, 8]. The following antibodies were used: rabbit anti-axonin-1 , mouse monoclonal anti-BrdU (1/100; Sigma), mouse monoclonal (Roche, Basel, Switzerland) and rabbit (Chemicon, Temecula, CA, USA) anti-GFP (1/400), mouse anti-Islet1,2 (1/100) , rabbit anti-Islet1 (1/500; Abcam, Paris, France), mouse anti-Lhx1,5 (Developmental Studies Hybridoma Bank), mouse anti-NeuN (1/500), rabbit anti-Sox2 (1/1,000; Abcam), rabbit anti-Sox13  and adequate fluorescent secondary antibodies. For NeuN and Sox2 staining, sections were pre-treated with sodium citrate pH 6, 0.05% Tween-20, and for BrdU staining, sections were pretreated in 2N HCl, 0.5% Triton X-100 for 20 minutes at room temperature and neutralized. Apoptotic cells were detected using the Apo Tag kit (Oncor, Gaithersburg, MD, USA)) according to the manufacturer's instructions.
green fluorescent protein
hours after electroporation
Hamburger and Hamilton stage
insulin-like growth factor
in situ hybridization
significance analysis of microarray
We like to thank F Guillemot for providing Mash1 mutant mice, ET Stoeckli for providing us with anti-axonin-1 antibodies, P Bovolenta for the chicken Sfrp1 probe and E Matsunaga for the chicken Pax2 probe. The monoclonal anti-Lhx1,5 antibodies developed by T Jessell were obtained from the Developmental Studies Hybridoma Bank developed under the auspices of the NICHD and maintained by the University of Iowa, Department of Biological Sciences. This work was supported by grants from Agence Nationale de la Recherche, Association Française contre les Myopathies, Association pour la Recherche sur le Cancer and the Fondation Franco-Norvegienne pour la Recherche Scientifique et Technique.
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.