EMX1 regulates NRP1-mediated wiring of the mouse anterior cingulate cortex

The anterior cingulate cortex in humans and mice is particularly important for the performance of executive functions and emotional processing (Gasquoine, 2013; Kim et al., 2014; Meechan et al., 2015). It comprises multiple subdivisions, each with distinct functions that are conferred by differing input and output projections (Vogt and Paxinos, 2014). However, little is known about the development and molecular regulation of the input and output projections of the anterior cingulate cortex.

Homotopic commissural projections from the anterior cingulate cortex are the first axons to cross the midline during corpus callosum development. These cells are primarily born between embryonic day (E) 13 and E14 in mice, and project their axons across the midline at E15.5 (Rash and Richards, 2001; Choe et al., 2012). As the first axons to navigate the complex midline environment during callosal development, these axons are thought to pioneer the initial trajectory that subsequently guides neocortical callosal axons across the midline from E16.5 (Ozaki and Wahlsten, 1998; Rash and Richards, 2001; Richards et al., 2004). Thus, defects in cingulate callosal axon pathfinding are likely to result in malformation of the corpus callosum, which is a relatively common neurodevelopmental disorder in humans (Richards et al., 2004; Paul et al., 2007; Edwards et al., 2014).

Neuropilin 1 (NRP1) and its ligands semaphorin (SEMA) 3A and SEMA3C are involved in the axon guidance of cingulate callosal neurons. NRP1 is expressed on cingulate callosal axons in both humans and mice (Ren et al., 2006; Piper et al., 2009). SEMA3A, which is expressed in a high lateral to low medial gradient in the neocortex during callosal development (Tamamaki et al., 2003; Zhao et al., 2011), is chemorepulsive to cingulate callosal neurons in mice (Piper et al., 2009). By contrast, SEMA3C is expressed by glutamatergic guidepost neurons at the midline, and is attractive to these axons (Niquille et al., 2009; Piper et al., 2009). We have also reported previously that cingulate callosal axons are misguided in mice with a mutated NRP1 receptor that lacks the semaphorin-binding domain (Nrp1^Sema− knock-in mice), demonstrating that these molecules are required for the correct guidance of cingulate axons across the midline (Piper et al., 2009).

The empty spiracles homeobox (EMX) transcription factors, EMX1 and EMX2, have been implicated in the development of the anterior cingulate cortex in mice (Pellegrini et al., 1996; Piper et al., 2009). EMX2 regulates the development of the caudal and medial regions of the telencephalon (Pellegrini et al., 1996; Yoshida et al., 1997; Bishop et al., 2000; Mallamaci et al., 2000), whereas much less is known about the function of EMX1 in forebrain development, as Emx1 knockout brains appear relatively normal (Qiu et al., 1996; Yoshida et al., 1997). However, as Emx1; Emx2 double knockout mice display a more severe phenotype than either single mutant alone, it is likely that the EMX transcription factors function cooperatively in some aspects of cortical development, and that EMX2 may act as a compensatory factor in the absence of EMX1 (Yoshida et al., 1997; Bishop et al., 2002,, 2003; Shinozaki et al., 2002; Shinozaki et al., 2004).

Originally, two Emx1 knockout mouse strains were generated in different laboratories, using either 129/SvJ or (F)1 CBA/C57Bl/6 embryonic stem cells, and both strains exhibit agenesis of the corpus callosum (Qiu et al., 1996; Yoshida et al., 1997). In the most severe cases, callosal axons that normally formed homotopic connections in the contralateral hemisphere did not cross the midline but projected rostro-caudally within the ipsilateral hemisphere to form aberrant longitudinal bundles called Probst bundles (Probst, 1901; Ren et al., 2007). Subsequently, Guo et al. (2000) generated a third Emx1 knockout strain using embryonic stem cells derived from the 129/Sv background and backcrossed this strain onto the C57Bl/6 background. They reported that, although the average brain weight of these mice was reduced, the corpus callosum appeared grossly normal, the exception being that its surface area in coronal sections was significantly smaller (Guo et al., 2000). Such phenotypic differences might be due to differences in the targeted insertion constructs used to generate each strain, or more likely the differing genetic backgrounds of these strains (Beck et al., 2000; Sellers, 2012). These conflicting results have led to uncertainty regarding the role of Emx1 in corpus callosum development.

Here, we report that EMX1 does indeed function to regulate corpus callosum development on the inbred C57Bl/6 background. By crossing the original Emx1 knockout mice generated by Qiu et al. (1996) onto the C57Bl/6 background, we find that EMX1 is required for the midline crossing of a subpopulation of callosal axons originating from the anterior cingulate cortex, and that these axons form aberrant projections in the absence of EMX1. We demonstrate that EMX1 functions by upregulating the expression of the gene encoding the axon guidance receptor NRP1, thereby establishing a novel regulatory mechanism of NRP1/SEMA signaling in corpus callosum development.

The role of EMX1 in corpus callosum formation is unclear, as Emx1 knockout mice generated on different background strains display differing callosal phenotypes (Qiu et al., 1996; Yoshida et al., 1997; Guo et al., 2000). To resolve these differences, we crossed one of the original Emx1 targeted deletion mutants that displayed agenesis of the corpus callosum (derived with 129/SvJ embryonic stem cells; Qiu et al., 1996) onto a C57Bl/6 background for over ten generations (hereafter referred to as Emx1^−/−/C57Bl/6). The Emx1^−/−/C57Bl/6 mice did not display any gross anatomical defects (similar to their pre-backcrossed cousins; Qiu et al., 1996). Congenic analysis of the genomic DNA of the Emx1^−/−/C57Bl/6 mice demonstrated that they were 95.94% identical to wild-type C57Bl/6 mice.

To analyze whether forebrain commissures cross the midline in adult Emx1^−/−/C57Bl/6 mice, gold chloride staining was performed on brains bisected at the midline, identical to the analysis performed by Guo et al. (2000). We observed no gross commissural defects of the corpus callosum, hippocampal commissure or anterior commissure in the mid-sagittal plane in any of the genotypes examined (Fig. 1A-C). However, as Probst bundles do not traverse the mid-sagittal plane (Probst, 1901; Ren et al., 2007), we also stained coronal sections for an axonal marker to determine if such bundles were present. Strikingly, neurofilament labeling of axons in serial coronal adult brain sections demonstrated the presence of aberrant axonal bundles in the ventral portion of the corpus callosum (Fig. 1F,I, arrows) that were localized to rostral sections (Fig. 1D-O). Thus, in contrast to the findings of Guo et al. (2000), we observed that adult Emx1^−/−/C57Bl/6 mice display a subtle but evident callosal dysgenesis phenotype that may provide crucial insight into the function of EMX1 in this system.

As only a subset of callosal axons appeared to be misrouted in Emx1^−/−/C57Bl/6 mice, we sought to identify whether this subset was derived from a single origin, or if it was associated with multiple cortical areas. To investigate this, we used high angular resolution diffusion imaging (HARDI) coupled with tractography to characterize the origin and projection of the aberrant axonal bundles (Moldrich et al., 2010). HARDI is a non-invasive magnetic resonance imaging (MRI) technique that allows the three-dimensional mapping of brain connectivity based on the diffusion of water molecules within axon tracts (Tuch, 2004; Hess et al., 2006). Within each voxel (100 μm³), the least-restricted direction of water diffusion can be interpreted as the predominant direction of white matter tracts, and represented by different colors on a colored fractional anisotropy (FA) map. Tractography programs then allow regions of interest (ROIs) covering a few voxels to be selected and used as initiation points (or seeds). Streamlines representative of the predominant axonal tract that passes through these ROIs are then reconstructed.

Tractography was initiated from ROIs within the aberrant bundles in adult Emx1^−/−/C57Bl/6 brains, or in equivalent regions of the ventral corpus callosum in Emx1^+/+/C57Bl/6 and Emx1^+/−/C57Bl/6 brains (indicated by the boxed regions in Fig. 2A-C). From sagittal and horizontal views, it is evident that tractography-generated streamlines represent a strong homotopic connection between the anterior cingulate cortices in Emx1^+/+/C57Bl/6 and Emx1^+/−/C57Bl/6 brains (Fig. 2D,E,G,H, green streamlines), whereas streamlines generated from equivalent ROIs in Emx1^−/−/C57Bl/6 brains project rostro-caudally to form aberrant axonal bundles (Fig. 2F,I, purple and blue streamlines). We also generated streamlines passing through an inclusion ROI in the midsagittal corpus callosum, initiated from seed ROIs corresponding to the anterior cingulate cortex, and the primary motor and somatosensory cortices (Fig. 2J-P). This analysis suggests that homotopic connections between the anterior cingulate cortices, but not the primary motor or somatosensory cortices, are reduced in Emx1^−/−/C57Bl/6 animals compared with controls (Fig. 2J-P).

To verify that the aberrant axonal bundles originated from the anterior cingulate cortex, we performed stereotactic iontophoretic injections of an anterograde axonal tracer, biotinylated dextran amine, into adult Emx1^−/−/C57Bl/6 mice and Emx1^+/+/C57Bl/6 controls. In Emx1^+/+/C57Bl/6 controls, callosal axons from the anterior cingulate cortex projected through dorsal regions of the corpus callosum to form homotopic connections in the contralateral hemisphere (Fig. 3B-D). Although similar projections were present in the Emx1^−/−/C57Bl/6 mice (Fig. 3A), the labeled cingulate projections were de-fasciculated at the midline, with more cingulate axons projecting through the ventral region of the corpus callosum (Fig. 3E-G) compared with wild-type controls. A subpopulation of labeled axons also contributed to the aberrant axonal bundle, but only within the region ipsilateral to the injection site (Fig. 3F,G, arrowheads). Collectively, these data indicate that ipsilateral projections from the anterior cingulate cortex contribute to the aberrant axonal bundles in Emx1^−/−/C57Bl/6 mice.

To understand more fully the mechanism underlying the formation of the aberrant axonal bundles in Emx1^−/−/C57Bl/6 mice, we first sought to determine when and how they arise. To investigate this, we performed immunohistochemical analysis for the axonal marker GAP43 on mice at either E17 or postnatal day (P) 0, when the majority of callosal axons would have just crossed the midline. Sections of developing brain from the pre-backcrossed Emx1^−/−/129/SvJ strain demonstrated agenesis of the corpus callosum and the presence of Probst bundles (Fig. 4A,B). In the post-backcrossed Emx1^−/−/C57Bl/6 strain, the aberrant axonal bundles were present in the ventral portion of the corpus callosum in rostral sections, similar to the findings in adult Emx1^−/−/C57Bl/6 brains (Fig. 4C-E). These results demonstrate that, during cortical development, EMX1 is required for the interhemispheric connectivity of a specific subpopulation of callosal axons.

De-fasciculation of cingulate callosal axons has also been reported in mice with specific deletion of Nrp1 in the primary motor cortex (Zhou et al., 2013), very similar to our observations in Emx1^−/−/C57Bl/6 mice. We therefore hypothesized that NRP1/SEMA signaling might be disrupted in the Emx1^−/−/C57Bl/6 mice. To investigate this, we first sought to determine how these proteins are expressed in vivo. We previously showed that NRP1 is expressed in cingulate callosal axons during initial callosal formation (E15-E17; Piper et al., 2009). EMX1 is expressed within the anterior cingulate cortex during a similar period (E14-E16; Fig. 5A-C). In a reporter mouse line expressing tdTomato in Emx1-positive cells (Emx1-iCre; ROSA26CAG-lox-STOP-lox-tdTomato; Kessaris et al., 2006; Madisen et al., 2010), Nrp1 mRNA was observed in tdTomato-positive cells in the anterior cingulate cortex at E15, suggesting that EMX1 and NRP1 are co-expressed within cingulate neurons (Fig. 5D-F). To determine if Nrp1 mRNA expression is altered in the anterior cingulate cortex of Emx1^−/−/C57Bl/6 mice, we dissected E16 brain tissue from the cingulate cortex for real-time quantitative PCR (qPCR) analyses. In the Emx1^−/−/C57Bl/6 mice, Nrp1 expression was significantly reduced in Emx1^−/−/C57Bl/6 mice to ∼40% of the levels observed in their wild-type littermates (Fig. 5G).

We next used immunohistochemistry to verify that NRP1 protein expression was similarly reduced, staining E16 coronal brain sections of Emx1^−/−/C57Bl/6 and Emx1^+/+/C57Bl/6 littermates. In the Emx1^−/−/C57Bl/6 mice, the expression of NRP1 was reduced in the cingulate cortex and at the midline where cingulate callosal axons were projecting into the contralateral hemisphere (Fig. 5H-M). This did not appear to be due to an overall decrease of NRP1 expression throughout the brain as NRP1 expression was retained caudally within the hippocampal commissure (Fig. 5M). Thus, Emx1^−/−/C57Bl/6 mice display reduced NRP1 expression in cingulate callosal axons that are forming the corpus callosum.

To determine whether reduced NRP1 expression in Emx1^−/−/C57Bl/6 cingulate callosal axons disrupts NRP1/SEMA3C signaling, we performed co-culture assays examining both neurite outgrowth and guidance in explants of cingulate cortex from E15.5 Emx1^−/−/ C57Bl/6 versus Emx1^+/+/C57Bl/6 mice challenged with HEK293T SEMA3C-overexpressing or control cell blocks. Axons from Emx1^+/+/C57Bl/6 cingulate cortex were significantly attracted to SEMA3C, as previously shown (Fig. 6A,B,E; Niquille et al., 2009; Piper et al., 2009), whereas this attraction was abolished in explants from Emx1^−/−/C57Bl/6 mice (Fig. 6C-E). No significant difference in outgrowth was observed between the two genotypes (Fig. 6F). These data demonstrate that the NRP1/SEMA3C-mediated guidance of cingulate axons is disrupted in Emx1^−/−/C57Bl/6 mice.

The reduced expression of Nrp1 mRNA and protein in Emx1^−/−/ C57Bl/6 mice suggests that EMX1 acts as a transcriptional activator of Nrp1 expression. To test this, we identified potential EMX1 binding sites within the Nrp1 locus using the web-based tool FIMO (Grant et al., 2011) and consensus EMX1 binding sites represented as position weight matrices (Fig. S1A; Jolma et al., 2013). As EMX1 and NRP1 appear to be similarly required for callosal development in humans (Ren et al., 2006), we examined the Nrp1 promoter and 5′ untranslated regions (5′ UTR) shared between mouse and human, and identified four highly conserved sequences containing putative monomeric and dimeric EMX1 binding sites (Fig. 7A; Fig. S1B). To determine whether EMX1 is able to regulate Nrp1 transcription via this region, we cloned the mouse and human Nrp1 promoters encompassing these sequences to drive a firefly luciferase construct for dual-luciferase reporter assays in U251 MG cells. Overexpression of EMX1 resulted in a fivefold or 6.5-fold increase in Nrp1 promoter activity from the mouse or human Nrp1 promoter construct, respectively, compared with overexpression of a GFP control plasmid (Fig. 7B; Fig. S1C). This suggests that EMX1 overexpression is able to upregulate Nrp1 promoter activity in vitro.

In the full-length mouse Nrp1 promoter construct, three of the identified EMX1 binding sites were located within the 5′ UTR, with the remaining binding site located within the promoter region. To further investigate the role of individual EMX1 binding sites in inducing Nrp1 promoter activity, we designed truncated firefly luciferase constructs encompassing different regions of the full-length construct. When EMX1 was co-expressed in U251 MG cells in vitro, the luciferase construct containing the sole EMX1 binding site within the promoter (denoted ‘prom’ in Fig. 7B) upregulated Nrp1 promoter activity to levels resembling those of the full-length Nrp1 promoter, but this upregulation was abolished when the EMX1 binding site was removed (‘prom Δs1’ in Fig. 7B). Nrp1 promoter activity was also significantly upregulated in the truncated promoter and 5′ UTR construct containing all four EMX1 binding sites (denoted ‘5′ UTR+s1’), as well as the truncated 5′UTR containing three EMX1 binding sites (denoted ‘5′ UTR’). Together, these results suggest that, although the EMX1 binding sites may play an additive role in regulating Nrp1 promoter activity, the binding site within the Nrp1 promoter (‘s1’ in Fig. 7B) is most important for upregulating Nrp1 promoter activity in vitro.

Our findings suggest that EMX1 regulates the expression of Nrp1 that is required for midline crossing of cingulate axons. To confirm this in vivo, we specifically removed NRP1 in the cingulate cortex by electroporating Nrp1^flox/flox mice with a Cre-expressing plasmid at E14.5 (just prior to when cingulate callosal axons begin to project across the midline; Rash and Richards, 2001). Cingulate callosal axons electroporated with either EGFP or Cre-EGFP had projected across the midline and into the contralateral hemisphere by P5 in Nrp1^WT/WT mice (Fig. 8A,C). This was also evident in the cingulate axons of Nrp1^flox/flox mice electroporated with EGFP (Fig. 8B). However, the number of callosal axons crossing the midline and projecting into the contralateral hemisphere was reduced in Nrp1^flox/flox mice electroporated with Cre-EGFP (Fig. 8D,E). A proportion of these axons instead projected laterally out of the cingulate cortex and away from the midline (Fig. 8D, arrow). These data demonstrate that NRP1 is required for the appropriate guidance of cingulate axons across the corpus callosum in vivo, and that a similar reduction of NRP1 expression in the cingulate axons of Emx1^−/−/C57Bl/6 mice is likely to contribute to the misguidance of cingulate axons observed in these mice.

To test directly whether the loss of NRP1 expression downstream of EMX1 was responsible for the pathfinding defects observed in the Emx1^−/−/C57Bl/6 mice, we performed a rescue experiment by over-expressing Nrp1 in the cingulate callosal neurons of these mice. To achieve this, cingulate neurons in Emx1^−/−/C57Bl/6 mice were electroporated at E14 with either Nrp1 or YFP control plasmid, and the number of sections displaying crossing cingulate axons were blindly quantified at E17 (Fig. 8G,I). In control YFP electroporations of Emx1^−/−/C57Bl/6 mice, we observed that 76% (62/82) of sections contained YFP-positive cingulate axons that crossed the midline by E17, substantiating our anterograde tracing experiments that demonstrated significant midline crossing in these mice. Despite this, in Emx1^−/−/C57Bl/6 mice electroporated with Nrp1, we observed a significant increase in the number of sections in which Nrp1-electroporated axons crossed the midline (88%; 72/82 sections) compared with that in control YFP-electroporated brains (P<0.05, one-sided Fisher's exact test, n=18 per condition). These data provide evidence that restoring Nrp1 expression to cingulate neurons at E14 is sufficient to partially rescue the cingulate axon guidance phenotype of Emx1^−/−/C57Bl/6 mice. We also found that Emx1^+/+/C57Bl/6 cingulate neurons electroporated with Nrp1 at E14 showed no differences in pathfinding by E17 when compared with YFP electroporations (Fig. 8F,H), suggesting that ectopic over-expression of Nrp1 does not disrupt midline guidance of cingulate callosal neurons.

The anterior cingulate cortex is part of a forebrain network that is required for processing emotion-related behaviors, including aggression, depression, self-esteem, pain and stress (Vogt, 2005; Weber et al., 2006; Toth et al., 2012; Yamashita et al., 2013; Agroskin et al., 2014). The molecular regulation of its development, and particularly its role in corpus callosum development, are only now becoming evident (Niquille et al., 2009; Piper et al., 2009). Here, we have provided evidence that Emx1 knockout mice display a subtle callosal dysgenesis phenotype that is associated with this brain region. Our findings demonstrate a novel role for EMX1 in regulating the pathfinding and fasciculation of callosal axons originating from the anterior cingulate cortex.

Our findings in the present study are in contrast to those of Guo et al. (2000), who analyzed a similar Emx1 null mutant on the C57Bl/6 background. These differences may be explained by discrepancies in the method of analysis. Specifically, Guo and colleagues used Nissl-stained coronal sections and gold chloride staining of mid-sagittally bisected brains for their analysis. Similar to their findings, we also observed no phenotype when we analyzed our Emx1^−/−/C57Bl/6 brains with gold chloride staining. However, upon further analysis with an axonal marker, we were able to demonstrate aberrant axonal bundles that originated from the anterior cingulate cortex. Similar analyses of the strain studied by Guo et al. (2000) would determine if it shares an identical phenotype.

Taken together with the previous studies that analyzed Emx1 null mutants on mixed backgrounds (Qiu et al., 1996; Yoshida et al., 1997), our findings demonstrate that the severity of callosal dysgenesis in these mice is dependent on genetic background. In humans, different callosal malformations have been observed in patients with similar genetic etiologies, suggesting that modifier genes are likely to influence the severity of this condition (Schaefer et al., 1997; Dastot-Le Moal et al., 2007; Edwards et al., 2014). The effects of genetic background have also been previously described for a number of other genes in the mouse central nervous system, including Scn1a (Miller et al., 2014) and Gnaz (van den Buuse et al., 2007). In these cases, quantitative trait locus mapping, coupled to gene expression and phenotypic analyses have identified the modifier loci and encompassing genes that might mediate these effects. Similar experiments using the Emx1^−/−/C57Bl/6 and Emx1^−/−/129/SvJ mouse models could identify the modifier loci influencing EMX1-mediated callosal dysgenesis. In turn, genes within the modifier loci themselves may potentially be demonstrated to influence callosal dysgenesis in humans and mice.

Axon targeting and de-fasciculation defects that are due to a loss of NRP1/SEMA signaling are well documented for neocortical callosal axons in vivo (Gu et al., 2003; Hatanaka et al., 2009; Zhao et al., 2011; Zylbersztejn et al., 2012; Zhou et al., 2013). Together with recent studies (Niquille et al., 2009; Piper et al., 2009), our findings demonstrate that NRP1/SEMA signaling regulates similar events during the development of cingulate callosal axons. Interestingly, defects in callosal axon pathfinding and fasciculation have been hypothesized to also affect callosal targeting within the contralateral cortex (Zhou et al., 2013; Fenlon and Richards, 2015). Thus, whether Emx1 knockout mice display defects in the final homotopic and heterotopic topography of callosal circuitry is of particular interest for further study.

The subtle phenotypes reported in Emx1 knockout mice as well as conflicting reports regarding whether EMX1 regulates callosal development have masked the importance of EMX1 in the developing mouse brain (Qiu et al., 1996; Yoshida et al., 1997; Guo et al., 2000). Specifically, how EMX1 functions independently of EMX2, particularly within postmitotic cortical neurons, which only express EMX1 and not EMX2 (Gulisano et al., 1996), is poorly understood. The present study demonstrates a novel role for EMX1 in regulating cingulate callosal guidance by the upregulation of Nrp1 expression. Even so, whether EMX1 binds to the Nrp1 promoter in vivo to regulate Nrp1 expression directly remains to be determined, as no anti-EMX1 antibodies have been reported for successful chromatin immunoprecipitation (including the antibodies we tested; see supplementary materials and methods). Our luciferase assays, however, suggest that EMX1 regulates Nrp1 expression via transcriptional activation of the Nrp1 promoter.

As EMX1 and NRP1 are expressed within both neocortical and cingulate callosal axons (Gulisano et al., 1996; Zhao et al., 2011; Zhou et al., 2013), the question that arises is why only a subpopulation of cingulate axons is misrouted in Emx1^−/−/C57Bl/6 mice. Further studies could investigate possible scenarios to explain this. Given that axonal growth cones integrate multiple cues to instruct pathfinding (Huber et al., 2003; O'Donnell et al., 2009), EMX1/Nrp1 signaling might simply increase the probability of cingulate axons crossing the midline, rather than being an absolute requirement for these axons to reach the adjacent hemisphere. Thus, defective EMX1/Nrp1 signaling in Emx1 knockout mice would increase the probability of misguidance, with the most severely affected axons failing to project across the midline.

Alternatively, or in concert with the previous scenario, callosal pathfinding in unaffected axons might be regulated through guidance mechanisms that function irrespective of NRP1 signaling. For instance, the Nrp1 homolog Nrp2 is expressed in the anterior cingulate cortex of E15 Emx1^+/+/C57Bl/6 and Emx1^−/−/C57Bl/6 mice (see Fig. S2), and NRP2 is similarly responsive to SEMA3C during cortical development (Giger et al., 2000). Whether NRP2 functions as a compensatory factor in Emx1^−/−/C57Bl/6 mice is of interest for future study. In addition, we recently demonstrated that cingulate callosal axons also express the axon guidance receptor DCC, and that these axons are attracted towards its ligand, NTN1, at the midline (Fothergill et al., 2014). Similarly, other short- and long-range axon guidance molecules have been demonstrated to function in regulating neocortical callosal pathfinding (Shu and Richards, 2001; Shu et al., 2003; Andrews et al., 2006; Islam et al., 2009; Nishikimi et al., 2011). As callosal axons from the dorso-lateral cortex, where Nrp1 expression is absent, are able to project their axons medially (Zhou et al., 2013), it appears likely that other axon guidance mechanisms can guide callosal pathfinding in the absence of Nrp1 expression.

Finally, EMX1 might be redundant in upregulating Nrp1 expression within the neocortex. Although EMX2 is unlikely to compensate for EMX1 in neocortical callosal neurons as it is not expressed in these cells (Gulisano et al., 1996), the nuclear receptor NR4A2 is reportedly expressed within the cortex during callosal development, and has been suggested to act as a transcriptional activator of Nrp1 (Zetterstrom et al., 1996; Hermanson et al., 2006). Whether NR4A2 compensates for the loss of EMX1 expression in Emx1^−/−/C57Bl/6 mice is therefore a question for future study.

The present study shows that EMX1 functions in corpus callosum development by specifically upregulating Nrp1 expression within the anterior cingulate cortex. Uncovering this mechanism has established several important regulatory features for further study. First, genetic modifiers that differ between background strains of mice can modulate EMX1 function. Second, EMX1 is required for cingulate axon pathfinding and fasciculation at the midline, with Emx1 mutants demonstrating the formation of aberrant axonal bundles. Finally, we have identified an animal model with specific anterior cingulate cortex wiring deficits that can be used in future studies to better understand the function of this circuit in behaviors such as depression and aggression.

Animals were maintained on a 12 h dark/light cycle with food and water provided ad libitum. All experiments were approved by the animal ethics committee of The University of Maryland School of Medicine, The University of Queensland or the Institute of Neuroscience, Shanghai Institutes of Biological Sciences, Chinese Academy of Sciences. The Emx1 knockout mice used were either on a mixed 129/SvJ/C57Bl/6 background (pre-backcrossed animals) or on a C57Bl/6 background (post-backcrossed). Post-backcrossed animals (Emx1/C57Bl/6) were generated by mating Emx1 knockout mice (B6;129X1-Emx1^tm1Jir/J strain; Qiu et al., 1996) with C57Bl/6 wild-type partners. This was repeated for over ten generations, with selection to retain the Emx1 mutation while alternating the C57Bl/6 wild type as the male or female between each generation. Heterozygotes from post-backcrossed animals were mated and the morning of the vaginal plug was designated E0. The Emx1-iCre; ROSA26CAG-lox-STOP-lox-tdTomato reporter mouse line was generated by crossing the Emx1-iCre mouse line (Kessaris et al., 2006) to the ROSA26CAG-lox-STOP-lox-tdTomato reporter mouse line obtained from the Jackson Laboratory (B6;129S6-Gt(ROSA)26Sor^{tm9(CAG-tdTomato)Hze}/J; Madisen et al., 2010). Nrp1^flox/flox mice (Gu et al., 2003) were a gift from Prof. Alex Kolodkin (Johns Hopkins University, MD, USA) and were maintained on a C57Bl/6 background. Genotyping of Emx1 knockout mice and Nrp1^flox/flox mice was performed using PCR primers as previously described (Qiu et al., 1996; Gu et al., 2003).

Adult mice were anesthetized with an intraperitoneal injection of sodium pentobarbitone (Lethabarb; Virbac; 185 mg/kg body weight) in PBS. Postnatal animals were anesthetized on ice. For embryo perfusion, the pregnant dam was anesthetized and embryos were sequentially removed from the uterus. Mice were transcardially perfused with 0.9% saline solution (0.9% w/v NaCl in Milli-Q water; Merck Millipore) (1 minute for embryos and postnatal animals; 5 minutes for adults) followed by 4% w/v paraformaldehyde (PFA; ProSciTech) in PBS (4 minutes for embryos and postnatal animals; 15 minutes for adults). All brains were post-fixed in 4% PFA and stored at 4°C until analysis.

Fixed adult Emx1/C57BL/6 brains were bisected along the midline and immersed in gold chloride solution as previously described (Guo et al., 2000; Wahlsten et al., 2003). For immunohistochemistry, fixed brains were set in 3% w/v Noble agar (Difco; Becton Dickinson) in PBS. Brains were sectioned coronally at 45 μm thickness or 200 μm thickness (data for Fig. 8A-E only) using a vibratome (Leica). Sections were then stained with the relevant antibodies as described in the supplementary materials and methods.

HARDI and tractography were performed as previously described (Moldrich et al., 2011). Briefly, Emx1 adult brains were incubated in 0.2% v/v Magnevist (Bayer HealthCare) over 4 days following fixation. Diffusion-weighted (DW) images were acquired with the brains immersed in Fomblin Y-LVAC fluid (Solvay Solexis), using a 16.4 Tesla vertical bore, small animal MRI scanner (Bruker Biospin; ParaVision v5.0) and a 15 mm linear surface acoustic wave coil (M2M Imaging). Three-dimensional 30 direction DW spin-echo sequences (b=5000 s/mm², 2×b=0) were acquired using the following parameters: repetition time=400 ms; echo time=22.8 ms; imaging resolution=0.1 mm isotropic. Reconstruction and tractography of the aberrant axonal bundles and equivalent callosal regions in controls were performed using Diffusion Toolkit and TrackVis (Ruopeng Wang and Van J. Wedeen, TrackVis.org, Martinos Center for Biomedical Imaging, Massachusetts General Hospital), according to the HARDI and Q-ball models (Hess et al., 2006). For characterization of homotopic callosal connectivity between the anterior cingulate, primary motor and primary somatosensory cortices, reconstruction and tractography was performed using MRtrix according to the constrained spherical deconvolution (CSD) model (Tournier et al., 2012). Cortical ROIs were defined by registering individual brains to a C57Bl/6J mouse neocortex atlas (Ullmann et al., 2013) using the FMRIB Software Library (Jenkinson et al., 2012). Tractography limits were set at FA values greater than 0.1 (fiber orientation density>0.1-0.25 for CSD tractography) and turning angle ≤45°.

Iontophoretic injection of the anterograde tracer was performed as previously described (Reiner et al., 2000), with an extended protocol delineated in the supplementary materials and methods. Briefly, an iontophoretic injector (Stoelting; 5 μA current pulses for 10 min; 7 s on, 7 s off) was used to deliver 15% w/v biotinylated dextran amine (10,000 MW; Life Technologies) in PBS into the anterior cingulate cortex through a glass micropipette (stereotactic coordinates: anterior posterior 1.42, lateral 0.2 and ventral −1.75).

In situ hybridization for Nrp1 was performed using an antisense riboprobe as previously described (Piper et al., 2009) with minor modifications. The fluorescein-labeled probe was detected with anti-Fluorescein-AP, Fab fragments (Roche) and amplified for visualization with the TSA Plus Fluorescein System (PerkinElmer). Sections were then stained for immunofluorescence as described to amplify the tdTomato fluorescence signal. In situ hybridization for Nrp2 was performed using a digoxygenin-labeled antisense riboprobe generated using the primers 5′-GCCTCTCTGTCAAAAATGGATATG-3′ and 5′-GCGTAGTCTGAGGTGAACTTG-3′, and visualized with BM Purple (Roche).

RNA was isolated from E16 cingulate cortex using TRIzol reagent (Invitrogen). Reverse transcription was performed using SuperScript III First-Strand Synthesis SuperMix (Invitrogen) with Oligo(dT) primers. Real-time qPCRs were performed using a Rotor-Gene 3000 (Corbett Life Science) and Platinum SYBR Green qPCR SuperMix-UDG (Invitrogen) with the following thermocycler conditions: 2 min at 50°C and 2 min at 95°C, followed by 50 cycles with 15 s denaturation at 95°C and 30 s annealing and extension at 60°C. Relative expression was determined using the ΔΔCt method with Actb as a housekeeping gene. The Nrp1 primers were 5′-GGTTTTCCATCCGCTATGAA-3′ and 5′-GGATTCGAGTCTTGCTCCAG-3′, and the Actb primers were 5′-CACGATGGAGGGGAATACAG-3′ and 5′-CAGCTTCTTTGCAGCTCCTT-3′.

The Sema3C expression construct in the pBK-CMV backbone was kindly provided by Prof. Andreas Püschel (Institut für Molekulare Zellbiologie, Westfälische Wilhelms-Universität Münster, Münster, Germany).

Plasmids of ubiquitin-EGFP and ubiquitin-Cre-2A-GFP were kindly provided by Prof. Zi-Long Qiu (Institute of Neuroscience, Chinese Academy of Sciences, Shanghai, China). Full-length Nrp1 in the pCAG-IRES-GFP backbone was kindly provided by Prof. Alex Kolodkin (Johns Hopkins University, MD, USA).

The in vitro axon guidance assay was performed and analyzed as previously described (Piper et al., 2009). At the completion of the culture period, explants were fixed in 4% paraformaldehyde in PBS and stained with anti-neuronal-specific beta III tubulin antibody (TuJ1 clone) (1:1000; R&D Systems Cat# MAB1195) and goat anti-rabbit Alexa Fluor 488 secondary antibody (1:500; Invitrogen Cat# A11008).

Ninety-six-well plates were individually seeded with 4500 U251 MG cells per well and maintained in DMEM (Life Technologies) supplemented with 10% v/v fetal bovine serum (SAFC Biosciences). After 24 h, a human or mouse Nrp1 promoter plasmid (or the empty firefly luciferase control plasmid) was co-transfected with either pCAG-Emx1-IRES-GFP or pCAG-IRES-GFP into these cells. The Renilla luciferase vector pRL-SV40 was co-transfected as an internal control in all experiments. Cells were lyzed after 40 h and luciferase activity was measured using the Dual-Luciferase Reporter Assay System (Promega) on a POLARstar OPTIMA plate reader (BMG Labtech). Cloning of luciferase DNA constructs are described in Table S1 and the supplementary materials and methods.

In utero electroporation was performed as previously described (Saito, 2006). Briefly, pregnant mice were anesthetized with 90 mg/kg (body weight) ketamine (Parnell Laboratories, Australia) and 10 mg/kg xylazine (Troy Laboratories, Australia). An abdominal incision was made to expose the uterus, and plasmid DNA diluted in 0.5 μl PBS was injected through the uterine wall into the medial region of the lateral ventricle of the brain with a glass pipette. Electrical pulses (five 35 V pulses of 50 ms applied at 1 second intervals) were then delivered to embryos by electrodes connected to a square-pulse generator [either a CUY21 (BEX) for experiments in Fig. 7 or a ECM 830 electroporator (BTX Harvard Apparatus) for experiments in Fig. 8]. Embryos were returned into the abdominal cavity and the abdominal cavity was sutured until embryo perfusion.

Immunohistochemical images were obtained with a Zeiss Z1 microscope (Carl Zeiss), using bright-field illumination at 2.5× magnification for adult sections and 5× magnification for embryonic sections. Immunofluorescence images were similarly obtained using 5× or 20× magnification. Gold chloride staining was imaged using a Zeiss Axiocam HrC camera (Carl Zeiss) attached to a stereodissecting microscope (Leica).

Sections were first imaged using a fluorescence microscope (Nikon E600FN) equipped with a Plan Apo 10× objective to obtain an image over a large area. Images were also acquired at a higher resolution using a confocal microscope (Nikon A1R equipped with a 40×, NA 1.30 lens) and were stitched over a large area for each brain section. The axon density ratio of callosal axon projections at the midline was quantitatively analyzed using ImageJ software as previously described (Zhou et al., 2013). The axon density ratio was obtained as the axon density at the midline normalized to that of the ipsilateral hemisphere.

This article is dedicated to the late Prof. Giorgio Corte (National Institute for Cancer Research and Department of Oncology, Biology and Genetics, University of Genova Medical School) who provided important reagents for this project. We are also grateful to Prof. Alex Kolodkin, Prof. Andreas Püschel, Prof. David Ginty and Prof. Zi-Long Qiu for providing important mouse lines and reagents; the Australian Genome Research Facility, Dr J Bertran-Gonzalez, Dr Oressia Zalucki, Chantelle Reid, Ilse Buttiens, John Baisden, Kim Valentino and Sarah Croft for technical assistance; Sarah Preston and Lynette Knowles for early work on the project; and Laura Morcom and Rowan Tweedale for reading the manuscript and providing critical comments.