Skip to main content

Advertisement

Download PDF

Hierarchical clustering of gene expression patterns in the Eomes + lineage of excitatory neurons during early neocortical development

BMC Neuroscience volume 13, Article number: 90 (2012) Cite this article

Abstract

Background

Cortical neurons display dynamic patterns of gene expression during the coincident processes of differentiation and migration through the developing cerebrum. To identify genes selectively expressed by the Eomes + (Tbr2) lineage of excitatory cortical neurons, GFP-expressing cells from Tg(Eomes::eGFP) Gsat embryos were isolated to > 99% purity and profiled.

Results

We report the identification, validation and spatial grouping of genes selectively expressed within the Eomes + cortical excitatory neuron lineage during early cortical development. In these neurons 475 genes were expressed ≥ 3-fold, and 534 genes ≤ 3-fold, compared to the reference population of neuronal precursors. Of the up-regulated genes, 328 were represented at the Genepaint in situ hybridization database and 317 (97%) were validated as having spatial expression patterns consistent with the lineage of differentiating excitatory neurons. A novel approach for quantifying in situ hybridization patterns (QISP) across the cerebral wall was developed that allowed the hierarchical clustering of genes into putative co-regulated groups. Forty four candidate genes were identified that show spatial expression with Intermediate Precursor Cells, 49 candidate genes show spatial expression with Multipolar Neurons, while the remaining 224 genes achieved peak expression in the developing cortical plate.

Conclusions

This analysis of differentiating excitatory neurons revealed the expression patterns of 37 transcription factors, many chemotropic signaling molecules (including the Semaphorin, Netrin and Slit signaling pathways), and unexpected evidence for non-canonical neurotransmitter signaling and changes in mechanisms of glucose metabolism. Over half of the 317 identified genes are associated with neuronal disease making these findings a valuable resource for studies of neurological development and disease.

Background

Excitatory, or projection, neurons comprise ~85% of the neuronal population in the mature cerebral cortex. These neurons share a stereotypical developmental trajectory, with most generated in the neocortical ventricular zone (VZ) and subsequently migrating outward, through complex and changing micro-environments, towards the cortical surface. During this radial migration neurons adopt multiple morphologies and display several distinct modes of migration [1].

The period of radial migration and early differentiation coincides with the sequential expression of a network of transcription factors. Core elements of the network include Pax6, Neurogenin, Neurod, Eomes (Tbr2), and Tbr1 [2–4]. Some of the network members are known to be essential for proper migration [5, 6], axogenesis [7, 8], dendritogenesis [9], and adoption of neuronal subtype [10–12]. The ability of neurons to simultaneously navigate through the developing cortex and acquire characteristic phenotypes is thus likely to depend on dynamic patterns of gene expression during the early postmitotic period.

Characterization of these cell-intrinsic, dynamic gene expression patterns remains incomplete. Transcriptional profiling of microdissected samples derived from developing brain, although revealing regional and temporal differences across the brain [13–15] is limited by the inclusion of blood vessel endothelial cells, inhibitory neurons, and glia within the dissected tissue: the derived transcriptional profile necessarily reflects contributions of multiple, and often non-neuronal, cell types. These interpretive obstacles can in principle be eliminated with a combined application of cell sorting from acutely dissociated brain [16] and subsequent transcriptional profiling [17, 18].

In this study fluorescence activated cell (FAC) sorting was used to purify neurons of the excitatory cortical lineage, from which RNA was isolated for transcriptional profiling. We focused on the neuronal lineage that expresses Eomes (Tbr2), a transcription factor that is expressed by intermediate precursor cells [19] and has been genetically linked to microcephaly in mouse [20] and human [21]. The target neurons were acutely dissociated at E14.5 from a reporter transgenic mouse (Eomes::eGFP)Gsat in which immature neurons of the excitatory lineage express eGFP [22]. The Genepaint in situ database [23, 24], which provides comprehensive coverage of gene expression patterns in mouse embryos at E14.5 was then utilized to validate the transcriptional profiling, from which dynamically expressed genes were identified. A simple but robust technique was developed to quantify the optical density of the in situ patterns (QISPs) across the cerebral wall. These patterns were in turn hierarchically clustered into putative functional groups.

Methods

Mice

All animal procedures were approved by the IACUC of SUNY Upstate Medical University. Tg(Eomes::eGFP)Gsat mice (The Gene Expression Nervous System Atlas Project, NINDS Contract # N01NS02331 to The Rockefeller University, New York, NY) were mated to C57Bl6 (Taconic Labs, Germantown, NY) to generate embryos. The day of plug discovery was designated embryonic day 0.5 (E0.5). Eomesodermin (Eomes) – also known as Tbr2 – is a transcription factor of the T-box family expressed by intermediate precursor cells in the excitatory neuron lineage of the cortex [3, 25]. Eomes::eGFP mice contain a BAC transgene with 185 kb of 5’ and 18 kb of 3’ sequence surrounding the Eomes locus. In these mice selective GFP expression is observed in developing excitatory neurons indicating the fidelity of the large cis-regulatory elements in the BAC transgene [26].

Tissue dissociation and cell sorting

To generate RNA from enhanced Green Fluorescent Protein (eGFP) expressing neurons, dorsal neocortices from Tg(Eomes::eGFP)Gsat embryos were rapidly dissected at E13.5 and E14.5. The meninges were removed and the cortical tissue was dissociated in 0.25% Trypsin-EDTA solution (Invitrogen) for subsequent fluorescence activated cell (FAC) sorting (Figure 1). Cortices from 2 embryos were used per sorted sample and 3 samples were generated per developmental time point. Dissociated cells were kept in ice-cold buffer (PBS + 2% BSA) until FAC sorting, with the total elapsed time between dissection and sorting < 1 hr. Cell suspensions were sorted into GFP + and GFP- populations on a Becton Dickinson FACS Vantage Flow Cytometer Cell Sorter (SUNY Upstate Medical Flow Cytometry Unit). Sorting continued until a minimum of 100,000 GFP + cells were generated per sample. The total time required for sorting the 3 samples was less than 30 min. Cells were sorted directly into collection tubes containing RNALater™(Qiagen Inc., Valencia, CA) to minimize post-sort RNA degradation.

Figure 1
figure1

Isolation, profiling, and QISP analysis of GFP + (Eomes+) cortical neurons. a, Neocortex (NCX) containing GFP + neurons is dissected, and cells isolated with flourescence activated cell sorting (FACS). Following transcriptional profiling, genes are sorted by RMA and their level of fold up- or down-regulation compared to the FACS GFP- population. Verification of up-regulated genes is accomplished with quantified in situ hybridization (QISP), and the verified genes hierarchically clustered (see Methods). b, Fluorescence microscopic images of a section from an E14.5 Eomes::eGFP brain, displaying GFP + (green) and GFP- cells, that are blue due to Hoechst nuclear labeling. Subregions of Interest (sROIs) 1-20 are superimposed in white on the image. Scale bar, 25 μm. MZ, marginal zone; CP, cortical plate; IZ, intermediate zone; SVZ, subventricular zone; VZ, ventricular zone; LGE, lateral ganglionic eminence. c, Representative outcome of FAC sort of GFP + and GFP- cells. A GFP + (P3) and GFP- (P4) were isolated based on fluorescence (GFP-A) and side scatter (SSC-A). The boxes denote the relevant sorted populations. All in situ hybridizations were obtained from Genepaint.

Full size image

RNA isolation and quality assessment

High quality total RNA was isolated using the RNAeasy™ (Qiagen) kit and flash frozen. From a total volume of 35 μl, 1 μl of the RNAeasy fluid was run on an Agilent 2100 Bioanalyzer using the RNA PicoChip to assess RNA structural quality and quantity. Peak areas for 28 S and 18 S rRNA were determined. Degraded samples have reduced 28 S/18 S ratios and contain additional small RNA fragments that appear as smears in the gel photo and small humps in the electropherogram making them easily detectable. Samples were discarded if the total RNA profile displayed evidence of significant degradation, as indicated by a ribosomal RNA peak ratio (28 S/18 S) below 1.0 (Applied Biosystems/Ambion TechNote 11:2). No samples were lost due to RNA degradation provided the acutely dissociated cells were sorted into collection tubes that contained RNALater™. The typical yield of total RNA was ~700 ng per 100,000 cells. Approximately 100 ng of total RNA was used for transcriptional profiling.

Transcriptional profiling and quality assessment

The Affymetrix Gene 1.0 ST Array (Affymetrix, Santa Clara, CA) provides whole-transcript coverage of 28,853 mouse genes. Each gene is represented on the array by 27 probes that interrogate the full length of the gene, rather than 3’-biased interrogations used in prior technologies (e.g., the Murine U74v2 set). Purification of RNA, probe generation, chip hybridization and array scanning was performed at the SUNY Microarray Core facility. Amplified and terminal-labeled cDNA were generated by the WT Sense Target Labeling Protocol (Affymetrix). cDNA was hybridized to the array with the Affymetrix Fluidics Station 450, and scanned with the Affymetrix GeneChip® Scanner3000. Spike controls were used to determine the quality of the hybridization. A 3’/5’ probe signal ratio significantly greater or less than 1 was interpreted as indicating RNA degradation, in which case the experiment was repeated.

Quantified in situ patterns (QISPs) and network analysis

QISPs were generated by measuring the optical density of in situ images derived from Genepaint.org, a public database of in situ hybridizations performed on E14.5 mouse embryos [23]. The sampled regions of interest (ROI) spanned the cerebral wall from the marginal zone (MZ) to the VZ. A single ROI was first placed in rostral cortex from an in situ section that included olfactory bulb as an orienting landmark. This original ROI was then subdivided into 20 sub-regions (sROI1-20), each sub-region corresponding to 5% of the cortical wall, and the optical density of each sROI was measured. The measured sROI optical density was corrected by subtracting background optical density, and each sROI was then averaged with adjacent sROIs to smooth the QISP profile (Figure 1).

Hierarchical clustering of QISPS using Pearson’s correlation

The 20 optical density measurements corresponding to each QISP (sROI1-20) were imported into Multiple Experiment Viewer [27, 28] and normalized for grouping. Multiple Experiment Viewer can cluster data by Euclidean distance or Pearson’s correlation. Pearson’s correlation quantifies the covariance between samples producing values from -1 to 1. As such it is well suited to quantify the relationship between two QISPs. Importantly two QISP with identical profiles across sROIs1-20 but different magnitudes will have a correlation value of 1 and be tightly clustered, thereby facilitating the analysis of co-regulated genes independent of the absolute expression levels of each.

IPA analysis of gene groups

The IPA network analysis tool (Ingenuity Systems, Redwood City, CA) was used to analyze clustered genes. IPA is a literature-based analysis tool that uses relationships derived from reviewed journal articles as well as additional sources (e.g, KEGG, Gene Ontology) to generate networks of interactions within a data set. Clustered genes were imported and network relationships were ascertained using the “Direct and Indirect interactions” and included “Endogeneous Chemicals” (e.g., retinoic acid).

Results

General analysis of the E14.5 excitatory neuron transcriptome

Previous work from our lab indicated that although Eomes (Tbr2) protein is transiently expressed in the subventricular zone (SVZ), GFP from the transgene endures in post-mitotic neurons for several days after the downregulation of endogenous Eomes [22]. Thus the GFP + neuron population in this system includes not only IPCs from the SVZ, but also migrating and differentiating excitatory cortical neurons. FAC sorting purified this GFP + population to near homogeneity; analytic re-sorting of FAC sorted GFP + cells indicated a purity of > 99%. mRNA was successfully isolated from both the GFP + and GFP- cell populations, processed and profiled on the Affymetrix Gene 1.0 ST Array (Figure 1). To further improve the analytic stringency we restricted analysis to genes with an expression level of RMA (robust multichip average) ≥ 7 in the GFP + population. This threshold was based on an internal standard; for 10 eye-specific genes ([29]; Additional file 1) we calculated from the GFP + population a RMA value of 5.5 ± 0.7 (mean ± SD). All analyzed genes in this report are thus expressed at a minimum of 2.8 fold over this empirically-defined “not expressed” RMA level.

The total transcriptome of the E14.5 GFP + and E13.5 GFP + population was examined. We identified transcript that showed ≥ 3-fold difference (up or down) between E14.5 GFP + and E13.5 GFP + (Figure 2). The gene expression was similar at the two time points: we found 9 up-regulated and 35 down-regulated transcripts (Figure 2a, 2a’ and Additional file 2 and Additional file 3). We chose to focus our further analysis on the E14.5 GFP + data, as this is the same developmental time point used by the Genepaint mouse in situ database (http://www.genepaint.org). For the E14.5 excitatory neurons more than 13,000 unique mRNAs, out of ~25,000 represented by the 1.0 ST Array, were above the defined threshold (RMA = 7.0). Although numerous (i.e., 53% of all unique mRNAs are expressed above threshold) it is similar to recent, independent estimates of cortical neuron transcriptome size (66%, [30]; 60-70%, [31]). The difference between estimates may arise from distinct threshold assignments as well as differences in sampling technique; that is, purified FAC sorted embryonic neurons versus dissections of postnatal cortical tissue.

Figure 2
figure2

Comparison of RMA values between samples. a, a’, Comparison of RMA values between E14.5 GFP + and E13.5 GFP + cells. b, b’) Comparison of RMA values between E14.5 GFP + and E13.5 GFP- cells.

Full size image

In an effort to focus analysis upon genes most robustly expressed by the E14.5 GFP + excitatory neuron population, we identified 475 transcripts up-regulated ≥ 3-fold compared to the E13.5 GFP- reference sample (Figure 2b and 2b’ and Additional file 4). An additional 534 genes were identified as being down-regulated by ≥ 3-fold (Additional file 5). To validate the up-regulated genes, Genepaint.org, a public database of in situ hybridizations performed on E14.5 mouse embryos, was queried [23]. Of the original 475 up-regulated transcripts, 352 were represented in the Genepaint in situ database. Fort two of the up regulated transcripts not represented by Genepaint in situs encoded small nucleolar RNAs, 6 encoded RIKEN cDNA clones, 5 encoded predicted genes while the remainder of non-represented genes were not readily grouped and presumably reflected Genepaint program priorities. Twenty-four of the 352 (6.8%) in situ patterns either displayed no signal in the cerebral cortex or were represented by a single in situ section with incomplete sampling of the cortical wall. The remaining 328 in situ patterns were then processed for additional validation and grouping.

Validation of the analytic grouping technique

A method to quantify the optical density of the in situ patterns (QISPs) across the cerebral wall was developed (see Methods). The method of QISP generation is potentially subject to error due to variability of ROI positioning and histological heterogeneity between cortical areas and sections. To assess the precision of the QISP method, multiple measurements along the rostral caudal axis were performed using a set of E14.5 in situ hybridizations for Dab1 (Additional file 6), which at E14.5 is strongly expressed in the VZ and cortical plate (CP) [32], providing two distinct regions of peak signal. The Dab1 in situ profile was measured across the cerebral wall at 5 different points within the rostral cortex on a single in situ image (Additional file 6a). The Pearson’s correlation of these five distinct quantifications ranged from 0.94 to 0.99 with an average value of 0.97 (n = 5). This indicates that values obtained with the ROI approach were consistent and largely insensitive to the precise position of the ROI. Similarly, we compared the Pearson’s correlation of Dab1 signal from three different, sagittal in situ images derived from sections of rostral cortex along the lateral-to-medial axis (Additional file 6b). These correlations ranged from 0.87 to 0.97 with an average correlation of 0.93 (n = 4). Importantly, for both the rostral-caudal and lateral-medial planes the peak signals within the VZ and CP were always restricted to within 10% of the cerebral wall (i.e., ± 2 sROI). These tests validated the QISP analysis by indicating robust concordance of optical density measurements between ROIs from different positions in the rostral cortex and along the lateral-to- medial axis.

To objectively evaluate genes of interest we quantified the in situ pattern for Eomes itself, reasoning that up-regulated genes would display a QISP equivalent to, or localized more superficially than, the QISP derived for Eomes (Tbr2). This assumption is based on the accepted observation that, at this developmental stage, the most mature neurons are located more superficially in the cortex than immature neurons [33]. Eomes displayed little expression in ventricular sROIs1-2 and higher expression in the more superficial sROIs3-11. We therefore compared the average in situ expression in sROIs1-2 versus the average expression in sROIs3-20. Transcripts where this ratio exceeded 1.0 (i.e., where expression was strongest in the deep VZ) were considered invalid. This procedure identified 11 invalid genes which were then removed from the initial, post-Genepaint list of 328 genes. Thus 97% (n = 317) of the ≥ 3-fold up-regulated, RMA ≥ 7 transcripts displayed a Genepaint-derived in situ hybridization pattern that was spatially consistent with the Eomes + excitatory neuron lineage.

Analytic grouping of the E14.5 excitatory neuron transcriptome

To explore in greater detail groups of spatially co-expressed genes, the 317 QISP-validated genes were imported into the Multiple Experiment Viewer [27, 28], the signal (expression) values normalized by gene, and then hierarchically clustered. This procedure groups together similar QISPs, and can reveal groups of genes that are spatially co-expressed and thus potentially co-regulated [34].

The clustering of similar QISPs revealed several distinct spatial distributions of dynamically regulated genes in the radial dimension of the developing cerebral cortex (Figure 3). Approximately 70% of up-regulated QISPs within the Eomes + lineage display peak expression in the CP and upper intermediate zone (IZ). The majority of these QISPs displayed similar (i.e., tightly clustered) expression patterns, characterized by gradually increasing expression in the IZ and peak expression in the CP. The most heterogeneous expression patterns were observed in those QISPs with peak expression in the VZ and IZ, with a number of QISPs displaying dual expression peaks (e.g., Gdap1, Bcl6) or expression restricted to either the IZ or VZ (e.g., Neurod1 or Sstr2).

Figure 3
figure3

Identification and hierarchical clustering of QISPs. a, Hierarchical clustering of 317 QISPs representing transcripts expressed ≥ 3 fold higher in Eomes + (GFP+) neurons, compared to GFP- cells. From left to right: gene clusters denoted numerically (1-15); expression levels in the MZ + CP, IZ, and VZ (red indicates high expression, green indicates low expression) and gene number (1-317). b, Gene lists of the clustered QISPs from panel a. QISPs were imported into Multiple Experiment Viewer (MeV) and clustered using Pearson’s correlation (see Methods). White and shaded areas demarcate the different gene clusters.

Full size image

Four groups of genes – that is, four groups of clustered QISPs – were examined using a network analysis tool (IPA, Ingenuity Systems). The four groups were: (1) all 317 genes; (2) 44 “deep group” genes that achieve highest in situ expression in the VZ; (3) 49 “middle group” genes that achieve highest expression in the IZ and lower CP; and (4) 224 “superficial group” genes with greatest expression in the CP. This approach was adopted because neurons in the deep, middle, and superficial regions of the cortex have distinct morphologies and modes of migration [1] and thus can be expected to manifest dynamic regulation of genes related to morphology, navigation and migration.

Network analysis of all up-regulated genes

Network analysis of all 317 genes revealed 164 associated with neurological disease and 85 associated with psychological disorders; i.e., ≥ 52% of the genes have previously been implicated in human health (Table 1). One hundred thirteen of the genes were categorized as Nervous System Development and Function, indicating a substantial number (36%) of genes in this group having been previously identified and functionally characterized. The top “Molecular and Cellular Functions” (using the IPA terminology) include Cell to Cell Signaling and Interactions, Cellular Assembly and Organization, and Cellular Movement, with the top “Canonical Pathways” being Axonal Guidance Signaling and CDK5 signaling. Both of these signaling pathways are known to be active during cortical neuron migration and differentiation [35–37]. This network analysis showed that the entire ≥ 3-fold up-regulated gene category contains a large number of annotated genes associated with neuronal development, migration, differentiation, and neurological disease.

Table 1 Summary of IPA Ingenuity analysis for the upper, middle, and deep cortical zones
Full size table

The “deep group” of 44 up-regulated genes (14% of the total), which is spatially co-expressed with Eomes (2.2 fold up), achieved highest expression in the VZ (sROIs 1-10). This group includes clusters 1, 5, 7, 8 and 11-15 (Figure 3). Along with Eomes, the deep group is characterized by the expression of the transcriptional regulators Neurod4, Bcl6, Tanc2 and Myt1. Earlier studies focusing upon this deep area of cortex revealed (a) neurons with a long leading process that translocate out of the VZ [19, 38, 39] and (b) compact intermediate precursor cells in the SVZ [40]. Thus, at this point in development the GFP + cells in the VZ and SVZ are predicted to represent a mixture of intermediate precursors and postmitotic neurons [41]. Consistent with this prediction, IPA analysis identified DNA Replication, Recombination and Repair as a functional category for these group (p = 1.37E-04; Table 1, Additional file 7). Surprisingly, Carbohydrate Metabolism was also identified (p = 3.00E-05): the deep group includes genes involved in the synthesis (C1ql3 and Lpl) as well as level regulation (Chga, Gipr, Pcsk2 and Sstr2) of D-glucose. The top “Canonical Pathways” include Protein Kinase A Signaling and Cardiac β-Adrenergic Signaling. Additionally, a number of genes involved in intracellular calcium signaling (Ryr3, Camk2b, Cacna1h) showed elevated expression in this group. Network analysis of deep group genes thus suggests that the phenotypic transition from radial glial precursor to intermediate precursor/immature cortical neuron coincides with significant changes in mechanisms of glucose metabolism and second messenger signaling.

The “middle group” comprises 49 genes, or 16% of the ≥ 3 fold up data set. The middle group QISPs, comprising clusters 9 and 10, are maximal in the IZ, where developing excitatory neurons display a multipolar morphology. Multipolar neurons (MPNs) possess multiple fine processes that lack definitive orientation, and migrate relatively slowly in the radial direction [1]. The top “Molecular and Cellular Functions” in this group are Cellular Development and Cell Morphology (Table 1, Additional file 7). The middle group is dominated by the basic helix loop helix (bHLH) transcription regulators Neurod6 (71.5 fold), Neurod1 (30.9 fold), Nhlh1 (11.6 fold) and Neurod2 (7.3 fold), which represent some of the most highly up-regulated genes in the entire data set. These bHLH transcription factors have well-described functions in excitatory cortical neuron differentiation and development [6, 42]. The top “Canonical Pathways” include Axonal Guidance Signaling and Semaphorin Signaling in Neurons, as well as Inositol Phosphate Metabolism. Of particular interest within the middle group genes is cluster 9, consisting of 26 genes that includes seven with known chemotropic signaling functions in axon guidance. The expression of this set of axonal guidance genes could promote the generation and extension of axons that are observed to trail the migrating cortical neuron in the upper IZ and CP [39, 43]. In addition to the highly up-regulated transcription factors, the top up-regulated genes include Vglut1/Slc17A6 (19.4 fold) a glutamate transporter that identifies excitatory cortical neurons [44]. Network analysis of “middle group” genes thus indicated that developing excitatory neurons in the IZ are molecularly defined by a high expression of canonical bHLH neuronal transcription factors and axonal guidance molecules, particularly those of the Semaphorin signaling pathway.

The “upper group” of genes (n = 224) were observed in clusters 2, 3, 4 and 6. These genes, comprising 71% of the up-regulated gene population, achieved highest expression in the developing cortical plate (CP), an area that contains migrating neurons with an elongated or bipolar shape, and post-migratory differentiating neurons. Within the upper group the most prominent transcriptional regulators were Tbr1 (43.1 fold) [45], Bhlhb5 (19 fold) [46], and Hivep2/Schnurri2 (18 fold). Additional transcription factors in the upper group included Prdm8 [47], Satb2 [48], Sox5 [11], and Bcl11b/Ctip2 [49]. All of these transcriptions factors, with the exception of Schnurri2, have demonstrated roles in corticogenesis. Schnurri2 is, however, a BMP signaling effector, with a neuronal function suggested by the observations that mice deficient in Schnurri2 are hyperactive and exhibit hypersensitivity to stress [50].

Neurons within the forming CP migrate along radial glia fibers, from which they later detach and extend dendrites. It is therefore not surprising that the top “Molecular and Cellular Functions” in the upper group are Cell-To-Cell Signaling and Interaction, Cellular Assembly and Organization, and Cellular Movement (Table 1, Additional file 7). The top “Canonical Pathways” include Corticotropin Releasing Hormone Signaling and CDK5 Signaling. A role for CDK5 signaling mechanisms is established for many neurodevelopmental functions, including the migration of cortical neurons along radial glial fibers and subsequent cortical layer formation and dendritogenesis [51, 52]. In contrast the role(s) of Corticotropin Releasing Hormone Signaling in these neurons is less clear. The eight genes from the ≥ 3 fold up-regulated category that populate the Corticotropin Releasing Hormone Signaling category include four with known functions in brain development, including the cannabinoid receptor Cnr1 [53], the transcription factor Mef2C [54], and nitric oxide synthase Nos1 [55]. The remaining genes in this group include the protein kinases Prkar1β, Prkcε and Prkcζ Network analysis of “upper group” genes thus suggest maturing excitatory neurons in the developing CP express a number of regulatory and signaling factors, including the transcription factors Tbr1, Bhlhb5, Hivep2/Schnurri2, Prdm8, Satb2, Sox5, and Bcl11b/Ctip2.

Analysis of highly up- and down-regulated transcription factors

To explore the spatial relationships of key transcription factors (TFs) in the cortical lineage, we identified all ≥ 3 fold up- and down-regulated genes associated with the Gene Ontology descriptor “transcription, nuclear” and clustered their respective QISPs (Figure 4). These transcription/nuclear QISPs span the cerebral wall at E14.5, with subsets of QISPs showing peak expression in the VZ, IZ and CP (Figure 4a). We were initially surprised that Eomes itself was upregulated only 2.2 fold, a value below our threshold for analysis. There are several reasons for this low fold up value. The first is that GFP- cells express Eomes mRNA (RMA = 10.7). The presence of Eomes mRNA in the absence of GFP protein could imply post-transcriptional regulation of the 5’ UTR of the Eomes message that is included in the Eomes: :eGFP transgene [26]. Similarly, the presence of Eomes mRNA in GFP- cells might reflect the rapid up regulation of Eomes mRNA in a subpopulation of cells that has yet to express sufficient GFP protein for FAC detection and sorting into the GFP + group. Finally, Eomes mRNA is sharply down regulated as neurons leave the SVZ [3], while GFP protein expression persists and is expressed in all cortical zones superficial to the SVZ, at this developmental time point. Thus the GFP + cell population is likely larger than the Eomes mRNA + cell population, effectively reducing the average Eomes mRNA expression level in the GFP + population. Thus post-transcriptional regulation and/or the difference between sharp Eomes mRNA regulation and the delayed and extended GFP protein expression may account for this relatively low fold up level for this critical TF.

Figure 4
figure4

Transcripts encoding transcription factors. Functional groupings of QISPs based upon Gene Ontology (GO) molecular function and biological process annotations. a, Clustered QISPs of transcripts expressed ≥ 3 fold higher in Eomes + neurons compared to GFP- cells. The left-most dendrogram reveals the hierarchical arrangement of QISPs. b, Clustered QISPs of transcripts expressed ≥ 3 fold lower in Eomes + neurons compared to GFP- cells. a’ and b’.

Full size image

We first compared these genes to an earlier study of bHLH-family TFs denoted as the “Dorsal Telencephalon Network” [4]. Comparison to this independent data set of 13 TFs revealed 7 (54%) to be members of our ≥ 3 fold up-regulated category: Neurod1 (31 fold), Neurod2 (7.3 fold), Neurod6 (72 fold), Tbr1 (43 fold), Id2 (4.1 fold), Nhlh1 (12 fold) and Satb2 (12 fold) (Figure 4). Other dorsal telencephalon TFs outside the ≥ 3-fold up-regulated category include the early genes Pax6 (7.6 fold down-regulated), Neurog 1 (1.3 fold down-regulated), Neurog 2 (1.1 fold down-regulated) and Eomes itself (2.2 fold up-regulated). Although Pax6 expression in the GFP + population was below threshold (RMA = 6.2) Neurog1, Neurog2, and Eomes were not (RMA values of 7.6, 10.1, and 11.9, respectively). This suggests that the GFP + neurons in our study substantially down-regulate Pax6 relative to neuronal precursors, but do express the neuronal specifying TFs Neurog1, Neurog2 and Eomes. Additionally, two genes specific to cortical layer V-VI, Etv1 and Otx1 [56], were down-regulated in our data (5 fold and 2 fold, respectively), an observation supported by strong Genepaint in situ hybridization signals at E14.5 in the ventricular zone for both genes. Because postnatal expression of Etv1 and Otx1 is specific to cortical layer V-VI, this result suggests temporally differential (i.e., developmentally early and late) requirements for these TFs during cortical development. In contrast to the “Dorsal Telencephalon Network” we identified the layer II-IV marker Satb2 within our ≥ 3 fold up-regulated data set (12 fold), an observation confirmed by Genepaint-derived in situ hybridization. Consistent with the dorsal telencephalic identity of the Eomes + population, none of the nine genes associated with the “Ventral Telencephalon Network” [4] were observed in the ≥ 3 fold up-regulated population. By comparison to published studies of excitatory neuron differentiation [2–4, 57], this TF expression analysis developmentally positions the “maturation spectrum” of the isolated E14.5 GFP + neurons to a period after the expression of early dorsal specification genes (e.g., Pax6), coincident with neural specification genes (e.g., Neurog1, Neurog2), and prior to the expression of late layer-specific markers (e.g., Etv1 and Otx1).

We hypothesized that the highly expressed “middle” and “deep” group TFs are critical molecular regulators, perhaps functioning to simultaneously repress precursor-specific TFs and promote the expression of TFs required for proper differentiation of cortical excitatory neurons. To test this idea, a network analysis tool (IPA; see Methods) was used to explore within our dataset possible regulatory interactions among the dynamically expressed TFs (Figure 5). All TFs from the ≥ 3 fold up- and down-regulated populations were imported into the pathway analysis and the network of direct and indirect interactions was identified. All interacting genes that were identified, but absent from the populations of ≥ 3 fold up- and down-regulated TFs, were assigned to a group denoted “IPA Interacting Genes.”

Figure 5
figure5

Network analyses of dynamically regulated transcription factors (TFs) in differentiating cortical excitatory neurons, derived from the IPA-Ingenuity Network Tool (see Methods). a, ≥ 3 fold up- and ≥ 3 fold down-regulated TFs were imported along with endogenous chemicals (e.g., retinoic acid) to populate the network. Grey boxes collect the up-regulated genes expressed in the “upper,” “mid” and “deep” cortical regions (CP + MZ, IZ, VZ alone) as well as the down-regulated genes (VZ) and the group of IPA Interacting Genes. Gray lines that are numbered indicate interactions between cortical zones (upper, middle, and deep); light blue lines indicate direct interactions of upper, middle, and deep up-regulated genes, and down-regulated genes, with IPA interacting genes. Dark blue lines indicate interactions between IPA interacting genes (see Results). b, Tabular summary of interactions between cortical zones. 1, Interaction identified by IPA Ingenuity analysis.

Full size image

Approximately 430 connections were identified for TFs in the up- and down-regulated groups.

Only 4 direct interactions were reported between the “deep group” or “middle group” genes and either the “upper group” genes or the down-regulated genes (Figure 5). Numerous direct interactions, however (50 of 54), were detected that connected “middle” and “deep” genes to the “IPA Interacting Genes” group. Included in this group are signaling molecules that function non-cell autonomously in neuronal development: Notch signaling [58], Wnt-βcatenin signaling [59], VEGF [60], TGF β [37] and retinoic acid signaling [61]. Each of these signaling pathways has been shown to regulate either the production and/or differentiation of excitatory cortical neurons. Surprisingly, only 14, or ~3%, of all detected interactions were direct between any group of dynamically regulated TF genes (i.e., between the 4 lower boxes in Figure 5). The results of this network analysis support models of neurogenesis that predict complex interactions between specific cell autonomous (e.g. transcription regulatory networks) and non-cell autonomous signaling mechanisms (e.g., Notch, TGFβ signaling) during the maturation of cortical excitatory neurons.

Analysis of highly up- and down-regulated receptor genes

To more fully explore the role(s) of dynamically regulated receptors in the differentiation and migration of cortical neurons, we identified within the ≥ 3 fold up- and down-regulated gene populations those with the Gene Ontology descriptor “receptor activity,” and clustered their respective QISPs (Figure 6). The down-regulated receptors were enriched for signaling systems with established functions in neuronal precursors (Figure 6b, 6b’), including the Wnt-βcatenin signaling family members Sfrp2, Fzd8, Lrp4, and Wnt7a. Wnt-βcatenin signaling controls the production of neurons by neural precursors [62] and regulates intermediate precursor proliferation [59, 63]. The FGF receptors Fgfr2 and Fgfr3, and Notch family members Notchr1, Notchr2 and Notchr3 were also down-regulated. Notch receptors have multiple roles in neurogenesis and differentiation, including the regulation of neuronal precursor differentiation [58, 64], the promotion of glial fate [65], and dendritic growth [66, 67]. Fgfr2 is required for excitatory neuron development [68], with Fgfr3 required for attainment of normal cortical and hippocampal volume and caudal cortex development [69]. These results suggest that, along with a dynamic enhancement of pro-differentiation genes, genes responsible for the maintenance of precursor/progenitor/glial phenotypes are dyna mically suppressed at E14.5 in the lineage of excitatory cortical neurons upon exiting the VZ.

Figure 6
figure6

Transcripts encoding receptors. Functional groupings of QISPs based on Gene Ontology molecular function and biological process annotations. Graphical representation as in Figure 4. Clustered QISPs of transcripts expressed (a) ≥ 3 fold higher and (b) ≥ 3 fold lower in Eomes:: GFP + cells compared to eGFP- cells.

Full size image

Thirty four receptors were identified in the ≥ 3 fold up category (Figure 6a, 6a’). The top category of up-regulated receptors includes those with established roles in axonogenesis (8.14E-11 p-value): members of the Semaphorin signaling family (Nrp1, Sema6d, Plxna2, Sema7a, Plxna4, Sema4g, Plxnd1); Netrin receptors (Unc5d, Dcc); the Slit signaling receptor Robo1. Because they are observed in clusters 9 and 10 (Figure 3), many of these canonical chemotropic signaling receptors are co-expressed in the IZ. Specifically, Semaphorin signaling molecules are expressed by excitatory neurons in the IZ soon after exit from the VZ, with Netrin receptors and Robo1 expressed in more superficial domains. These chemotropic signaling systems have multiple roles and may be enabled for axonal growth [70–72], neuronal migration [73], dendritic growth [74–76] and cell survival [77].

Also included among the most highly up-regulated receptors were eight GTP-binding protein-coupled receptors (Gprs) (2.87 E-5 p-value), including Mcr4 (melanocortin receptor 4), Cnr1 (cannabinoid receptor1), Sstr2 (somatostatin receptor 2), Gipr (gastric inhibitory polypeptide receptor), and the muscarinic acetylcholine receptor Chrm3. Cnr1 has multiple developmental roles including axonal guidance [78]. Specific developmental functions of Mcr4, Gipr, and Chrm3 are unknown, although Mcr4 interacts with Cnr1 to regulate food intake in adults [79], Gipr regulates insulin secretion after feeding [80], and Chrm3 is required for pancreatic insulin and glucagon release [81] and neuronal endocannabinoid release [82]. The up-regulated Gprs also include the orphan receptor Gpr21, as well as a highly conserved receptor, Gpr85 (also known as Sreb2), known to control brain size [83]. The Flamingo homolog Celsr3, which is required for neuronal projections to the anterior commissure as well as subcerebral targets, is also up-regulated [84].

Analysis of highly up- and down-regulated neurotransmitter receptors and synaptic genes

Neurotransmitters and their receptors have important functions in the developing nervous system [85]. Although our dataset contained no neurotransmitter receptors down-regulated ≥ 3 fold, the up-regulated category contained a variety of neurotransmitter receptor subunits. Cholinergic receptor subunit β2 (Chrnb2) and Chrna7 are up-regulated in the IZ and CP by 4.5 fold and 2.8 fold, respectively. These receptor subunits form heteromeric receptor complexes in the developing brain, suggesting an early role for cholinergic signaling as excitatory neurons migrate out of the VZ. Although deletion of Chrna7 in mice does not cause overt structural abnormalities [86], recent studies indicate that haploinsufficiency in CHRNA7 may cause a range of neurodevelopmental phenotypes associated with the 15q13.3 syndrome, including developmental delay, mental retardation and seizure [87]. Although no glutamate receptors were upregulated ≥ 3 fold, N-methyl-D-aspartate (NMDA) receptor subunits Grina and Grin1 were mildly up-regulated (2.6 fold and 2.0 fold, respectively; data not shown). In addition to mediating adult forms of synaptic plasticity, NMDA receptors may mediate forms of spontaneous activity in the developing cortex [88] as well as structural plasticity in the late embryonic and early postnatal period [89, 90]. mRNA encoding GABA (gamma-amino-butyric-acid) receptor subunits were also up-regulated, including Gabarg2 (4.6 fold) and Gabbrb2 (3.9 fold). Gabarb3 was upregulated 2.8 fold. Deficiency in Gabarg2 may be causal in Dravet syndrome, a form of childhood epilepsy [91]. GABA functions as an excitatory neurotransmitter in early development, xand deficiency in GABAA receptor subunits alter Inhibitory Post Synaptic Potentials (IPSPs) but not synapse number [92]. During early periods GABA can function as a depolarizing neurotransmitter and may relieve Mg2+ blockade of NMDA receptors [93, 94]. One glycine receptor subunit, Glra2 (10.7 fold up-regulated) was identified in this group, although mice lacking Glra2 display no overt abnormalities of cortical structure [95].

Synaptogenesis begins quite early in cortical development, at ~ E15 in subplate neurons in the rodent temporal cortex [96], and continues well into the postnatal period [97]. At the developmental point of our analysis, E14.5, few if any synapses are present in the developing CP. Despite this ultrastructural feature, we detected not only multiple mRNAs encoding neurotransmitter receptors, but also message for 22 other synaptic genes identified by the Gene Ontology descriptor “synapse” or “synaptic” (Figure 7). Included in this group are synaptotagmin IV (SytIV), a calcium sensing protein involved in neurotransmitter release [98], and synaptic vesicle glycoproteins Sv2a and Sv2b that may also be involved in regulated secretion [99]. Additionally detected were three members of the Leucine Rich Repeat Transmembrane family (Lrrtm1-3), which function as synaptic organizers [100]. These results admit the possibilities that neurotransmitter receptors and synaptic mRNAs expressed by these immature excitatory neurons either (a) do not encode functional proteins, or (b) encode proteins involved in non-canonical functions. In either case, the results are consistent with important roles for classical neurotransmitter signaling molecules in shaping the developmental trajectory of differentiating excitatory cortical neurons.

Figure 7
figure7

Transcripts encoding synaptic proteins. Functional groupings of QISPs based on Gene Ontology molecular function and biological process annotations. Graphical representation as in Figure 4. Clustered QISPs of transcripts expressed (a) ≥ 3 fold higher in Eomes-GFP + cells compared to GFP- cells. Note: there were no downregulated synaptic transcripts.

Full size image

Analysis of highly up- and down-regulated voltage-dependent channels

Early patterns of electrical activity and calcium transients are known to shape the developing nervous system [101–104]. For this reason we queried the grouping of QISPs associated with known voltage-dependent channels (Figure 8). One down-regulated, and ten up-regulated, transcripts were identified and grouped. Of the up-regulated transcripts, four encode subunits of voltage-dependent calcium channels (Cacna1e, Cacna1h, Cacnb3, Cacna2d1), four encode voltage-dependent potassium channels (Kcna1, Kcnj3, Kcnq3, Kcna6), and two encode voltage-dependent sodium channels (Scn3b, Scn2a1). Hierarchical clustering of these transcripts’ in situ hybridization profiles revealed peak expression within the upper-most ten sROIs, corresponding to the IZ, CP and MZ. This expression pattern is spatially consistent with the acquisition of voltage-dependent ionic currents in multipolar and migrating neurons [103]. Only one gene (Kcnj10), encoding the weakly inwardly rectifying potassium channel subunit Kir4.1, was down-regulated (3.3 fold), and this observation was consistent with the corresponding in situ hybridization data. Although robustly down-regulated at this developmental time point by differentiating excitatory cortical neurons, we note that deficiency in KCNJ10 produces seizures, sensorineural deafness, ataxia, mental retardation, and electrolyte imbalance (SeSAME syndrome) [105, 106]. This analysis is consistent with an important contribution of electrical activity toward shaping the development of cortical neurons (and presumably their synaptic/functional networks), but also implicates specific channel types and sub-units in the manifestation of these early activity-dependent processes by cortical excitatory neurons.

Figure 8
figure8

Transcripts encoding ion channels and solute carriers. Functional groupings of QISPs based on Gene Ontology molecular function and biological process annotations. Graphical representation as in Figure 4. Clustered QISPs of transcripts expressed (a) ≥ 3 fold higher and (b) ≥ 3 fold lower in Eomes-GFP + cells compared to GFP- cells.

Full size image

Analysis of highly up- and down-regulated adhesion molecules

Spatial distributions of adhesion molecules, and changes in such distributions, can in principal convey substantial morphologic and migratory changes upon developing neurons. To group expression patterns of known adhesion molecules within the Eomes + neuronal lineage, we sorted within the ≥ 3 fold up-regulated category for the Gene Ontology descriptor “adhesion” (Figure 9). Twenty one genes were identified that were verified by screening the Genepaint database. These genes include several Immunoglobulin Super Family Members (IGSFMs), such as L1cam, Dscam, Nfasc, Cadm3, Cntn1 and Cntn2. All of these IGSFM genes have confirmed neurodevelopmental functions: L1cam is involved in neuronal migration and neurite outgrowth [107]; Dscam is a co-receptor for Netrin signaling and is required for commissural axon pathfinding [108]; Cadm3, also known as Nectin-like 1 or Syncam3, is a cell adhesion molecule involved in axonal myelination [109]; Nfasc and Cntn1 interact to promote the formation of paranodal axo-glial junctions during axonal myelination [110]; Cntn2, also known as Tag-1, is highly expressed on axons and can serve as a substrate for interneuron migration [111] (although Tag-1 knockout mice display no overt cortical phenotype [112], there is a disruption of paranodal myelin junctions [113]). Several Cadherins (Cadh8, Cadh11 and Cadh13) as well as Protocadherin 20 were observed to be up-regulated, although to date none of these Cadherin genes have identified roles in cortical development. Surprisingly, no members of the integrin family were observed in the ≥ 3 fold, or the ≥ 2 fold, up-regulated groups, suggesting this important family of adhesion proteins is not dynamically regulated at the transcriptional level in the Eomes + lineage at E14.5.

Figure 9
figure9

Transcripts encoding adhesion molecules. Functional groupings of QISPs based on Gene Ontology molecular function and biological process annotations. Graphical representation as in Figure 4. Clustered QISPs of transcripts expressed (a) ≥ 3 fold higher and (b) ≥ 3 fold lower in Eomes-GFP + cells compared to GFP- cells.

Full size image

In contrast other adhesion molecules, such as Fibronectin (Fn1), were down-regulated in the developing Eomes + lineage. Two integrins that via heterodimerization form a fibronectin receptor (and an activating integrin for TGFβ signaling), Itgv and Itgβ8, were also down-regulated. Conditional deletion of Itgαv in radial glia and neurons causes cerebral hemorrhage [114] while reduced expression of αvβ8 integrin is associated with arteriorvenous malformations [115]. Other down-regulated genes encoding adhesion molecules were: Jam3, a tight-junction protein essential for maintaining the integrity of the cerebrovascular endothelium [116]; Cdh6, which is expressed in the VZ where it may serve to assign the cortico-striatal boundary [117]; two additional IGSFM genes, Boc and Cdon, that are evident in the VZ, where they function in Sonic Hedgehog signaling and are implicated in forms of holoprosencephaly [118, 119]. Taken together these groupings suggest that expression (and suppressed expression) of numerous adhesion molecules, particularly specific members of the cadherin and IGSF families, are required elements for migration, neurite formation, and subsequent myelination of developing excitatory cortical neurons.

Discussion

At E14.5 the cerebral neocortex is composed primarily of neuronal precursors and immature neurons, with a small number of differentiated preplate neurons. To better identify specific patterns of gene expression within the excitatory neuron lineage during early cortical development, we combined FAC sorting and transcriptional profiling with a new, comprehensive in situ hybridization validation procedure, termed QISP. This approach revealed 317 validated genes in the ≥ 3 fold up category of transcripts. We then clustered the QISPs into groups of spatially co-expressed genes that can contribute to the combined processes of migration and differentiation. This approach allowed us to identify 44 genes expressed in the deep VZ that are likely to be IPC specific. These deep genes would be difficult to isolate by a microdissection based profiling and in situ approach [13–15] because of the mixture of cell types and maturity. Similarly, we identified a group of 49 genes that are expressed in the IZ and may contribute to the recently identified multipolar neuron stage cortical development [1], as well identifying 224 genes expressed in the developing cortical plate that in part contribute to the early morphologic differentiation of cortical excitatory neurons (Figure 10). This approach allows inferences about the temporal sequence of neuronal development: While the in situ analysis was performed at one time point (E14.5) and is formally interpreted as a spatial pattern of gene expression at that time point, an approximate developmental timeline of gene expression within the lineage can be inferred, based on the well-established observation that more superficial neurons in the developing mammalian cortex are more mature [33, 120]. This inference is particularly valid at early developmental time points, such as E14.5, when the CP is starting to form and few neurons have “settled” or been displaced to deeper positions by subsequent, later-born neurons.

Figure 10
figure10

Schematic of excitatory neuron migration, differentiation, and patterns of gene expression. a, Illustration of radial migration, and general morphological characteristics, of excitatory neurons. RGC, radial glial cell; IPC, intermediate precursor cell; MPN, multipolar neuron; BPN, bipolar neuron; PMN, postmigratory neuron; VZ, ventricular zone; IZ, intermediate zone; CP, cortical plate; MZ, marginal zone. b, Pie charts representing gene clusters 1-14 (see Figure 3) as a function of their relative contribution to the population of ≥ 3 fold up-regulated genes as a function of cortical location (upper, middle, deep). The accompanying table (bottom) indicates the Gene Ontology (GO) classification for each of the gene clusters, and the percentage of genes within each GO classification in each zone (D, M, U). c, The individual, greatest up-regulated genes in the Eomes lineage for each cortical zone.

Full size image

One of the initially surprising features of the data set is the large number of genes that exceeded our empirical expression threshold (RMA 7.0). We found that 53% of all unique mRNAs are expressed above this threshold. Although initially surprising, this result is consistent with prior studies of more mature cortical tissues [30, 31] that contain a mixture of fully differentiated neurons with differing molecular properties. Although mRNA expression is not equivalent to protein expression, the fact that newly-generated excitatory cortical neurons display a complex transcriptome points to the potential plasticity of these cells from their earliest developmental origins. The broad spectrum of expressed receptors indicates chemosensitivity to major chemotropic signaling mechanisms (Semaphorin, Ephrin, Netrin, and Slit) and classical neurotransmitters including GABA, Glutamate, Acetylcholine and Glycine. Immature excitatory neurons are therefore potentially broadly sensitive to environmental cues, and could harness these cues to guide, and ultimately define, their early differentiation and migration across the cerebral wall. While the large number of signaling pathways underscores the potential plasticity and adaptability of immature neurons, this observation also highlights their potential susceptibility to toxic insult and drugs of abuse.

The Eomes::eGFP transgene causes GFP expression for a defined interval (approximately 3 days) in the early postmitotic life of excitatory neurons. Our analyses showed that the gene expressed by neurons during this ~ 3 day period of postmitotic life, is similar regardless of whether the GFP + cells were isolated from E13.5 or E14.5 embryos (Figure 2a and 2a’). This presumably reflects the shared initial developmental trajectory of these early born neurons, regardless of birthdate. This finding supports the observation of a core excitatory neuron developmental program centered on the transcription factors Pax6, Eomes and Tbr1 [57]. Thus, while our analyses focused on the comparison of E14.5 GFP + cells to GFP- precursors, we expect that many (but not all) of the genes we identified as upregulated in the transition from precursor to immature deep layer neuron will also be upregulated during the differentiation of later born, upper layer cortical neurons.

The complex interplay between the developing neuron and its environment was highlighted in our analysis of the network of transcription factors (TFs) that drive excitatory neuron differentiation. Although we were able to identify a total of 74 dynamically regulated TFs in the excitatory lineage, only 3% of the IPA literature-based interactions connected these TFs between different spatial zones (i.e., “deep,” middle,” and “upper” cortical plate). The remaining ~97% of interactions connected dynamically regulated TFs to a group of “IPA Interacting Genes,” a collection that included a number of molecules involved in non-cell autonomous signaling, such as Tgfβ, Vegf, and Wnt- βcatenin signaling systems. While this analysis only reports known interactions, and is therefore likely to under represent the true number of functional TF-TF interactions, it is presumably also likely to under represent indirect molecular interactions. So instead of a network composed exclusively of sequentially active (or activated) TFs, our data is consistent with an emerging model of early cortical neuron differentiation that predicts interaction between cell autonomous and non-cell autonomous mechanisms [102, 121–124]. This model helps account for one of the paradoxical features of neuronal development: the stochastic cell type composition of clonally related neurons observed in, for example, the neural retina [125–127] and cortex [128, 129]. In contrast to the invariant lineages that would be expected from a strictly deterministic developmental program, the variable mixture of neuronal subtypes within each clone may indicate a “noisy” deterministic program and/or a developmental program that relies on non-cell autonomous cues (e.g., growth factors or neurotransmitters) that may change in functional abundance and/or localization during development.

Another important result emerging from recent work is that specific members of TF families control restricted sets of phenotypic attributes. For example, Neurog2 enacts a program of motility in migrating neurons in the IZ through activation of the Rho family GTPase Rnd2 [130]. Later in development a second bHLH TF, Ascl1, regulates the transition to locomotory migration on radial glial fibers via Rnd3 activation [131]. Additionally, the subcortical projection patterns of cortical neurons depends on the expression of Sox5 [11] and may involve other TFs as well [7, 132, 133]. This stepwise control of neuronal phenotype may enable the appropriate maturation of neurons, independent of the size of the developing cerebral cortex or the time required for neurons to migrate across the cerebral wall. During early corticogenesis, mouse layer 6 neurons can transition from cell cycle exit to a postmigratory differentiating neuron in ~1.5 days [22, 134], whereas layer 2/3 cortical neurons require 4-5 days for this transition [135]. Nevertheless, a similar, if not identical, TF program is thought to underlie the maturation of both classes of neurons [3, 57]. Developing neurons’ interactions with external cues that promote (or suppress) the expression of stage-specific TFs would thus permit the coordinated transitions of neuronal phenotype (e.g., multipolar neuron to bipolar migrating neuron) in a manner independent of a cell-intrinsic “timer” of differentiation.

In this regard the cluster of genes expressed in the IZ may be important for coordinating the multipolar neuron stage, an early stage of migration, and later axon initiation. In clusters 9 and 10 we observed a significant enrichment of the Gene Ontology categories Cellular Movement Signaling (p = 4.59E-05) as well as Axonal Guidance Signaling Molecules (p = 1.28E-07) and bHLH TFs including Neurod6 (72 fold) and Neurod1 (31 fold). Neurod6 in particular may have a role in regulating the mitochondrial mass increase and cytoskeletal rearrangements prior to axonal growth [136], while other bHLHs, including Neurod1, promote the expression of genes critical for cortical migration, including the microtubule associated protein Dcx, the Cdk5 regulatory subunit p35, and the small GTPase RhoA [137]. The multipolar neuron phenotype could thus be explained, in part, by cellular specification provided by Neurod1 and Neurod6, two of the most highly upregulated TFs in our data set.

Conclusion

In summary, this study provides a rigorously validated list of transcripts expressed by immature excitatory neurons that have been grouped by in situ expression pattern. We identified the transcriptional complexity of these cells, and highlighted the potential functionality of these neurons with regards to the expression of specific chemotropic receptors, neurotransmitter receptors, ion channels, adhesion molecules, and synaptic proteins. The established association of many of the validated genes with human neurological disorders revealed new molecular bases for understanding how specific genes, and functional networks of such genes, control neuronal development and contribute to pathology. Future challenges include relating identified TFs and TF networks to specific developmental events, and identifying the exogenous cues that coordinate these events.

Abbreviations

BAC:

Bacterial Artificial Chromosome

CP:

Cortical Plate

eGFP:

enhanced Green Fluorescent Protein

Eomes:

Eomesodermin

Gsat:

Gensat

IACUC:

Institutional Animal Care and Use Committee

IPC:

Intermediate Precursor Cell

IZ:

Intermediate Zone

MPN:

Multipolar Neuron

MZ:

Marginal Zone

QISP:

Quantified In Situ Pattern

ROI:

Region of Interest

sROI:

sub –Region of Interest

sVZ:

sub-Ventricular Zone

Tg:

Transgene

VZ:

Ventricular Zone.

References

  1. 1.

    Tabata H, Nakajima K: Multipolar migration: the third mode of radial neuronal migration in the developing cerebral cortex. J Neurosci. 2003, 23 (31): 9996-10001.

  2. 2.

    Schuurmans C, Armant O, Nieto M, Stenman JM, Britz O, Klenin N, Brown C, Langevin LM, Seibt J, Tang H, et al: Sequential phases of cortical specification involve neurogenin-dependent and -independent pathways. Embo J. 2004, 23 (14): 2892-2902. 10.1038/sj.emboj.7600278.

  3. 3.

    Englund C, Fink A, Lau C, Pham D, Daza RA, Bulfone A, Kowalczyk T, Hevner RF: Pax6, Tbr2, and Tbr1 are expressed sequentially by radial glia, intermediate progenitor cells, and postmitotic neurons in developing neocortex. J Neurosci. 2005, 25 (1): 247-251. 10.1523/JNEUROSCI.2899-04.2005.

  4. 4.

    Gohlke JM, Armant O, Parham FM, Smith MV, Zimmer C, Castro DS, Nguyen L, Parker JS, Gradwohl G, Portier CJ, et al: Characterization of the proneural gene regulatory network during mouse telencephalon development. BMC Biol. 2008, 6: 15-10.1186/1741-7007-6-15.

  5. 5.

    Castro DS, Martynoga B, Parras C, Ramesh V, Pacary E, Johnston C, Drechsel D, Lebel-Potter M, Garcia LG, Hunt C, et al: A novel function of the proneural factor Ascl1 in progenitor proliferation identified by genome-wide characterization of its targets. Genes & development. 2011, 25 (9): 930-945. 10.1101/gad.627811.

  6. 6.

    Ge W, He F, Kim KJ, Blanchi B, Coskun V, Nguyen L, Wu X, Zhao J, Heng JI, Martinowich K, et al: Coupling of cell migration with neurogenesis by proneural bHLH factors. Proc Natl Acad Sci U S A. 2006, 103 (5): 1319-1324. 10.1073/pnas.0510419103.

  7. 7.

    Weimann JM, Zhang YA, Levin ME, Devine WP, Brulet P, McConnell SK: Cortical neurons require Otx1 for the refinement of exuberant axonal projections to subcortical targets. Neuron. 1999, 24 (4): 819-831. 10.1016/S0896-6273(00)81030-2.

  8. 8.

    Mattar P, Britz O, Johannes C, Nieto M, Ma L, Rebeyka A, Klenin N, Polleux F, Guillemot F, Schuurmans C: A screen for downstream effectors of Neurogenin2 in the embryonic neocortex. Dev Biol. 2004, 273 (2): 373-389. 10.1016/j.ydbio.2004.06.013.

  9. 9.

    Aizawa H, Hu SC, Bobb K, Balakrishnan K, Ince G, Gurevich I, Cowan M, Ghosh A: Dendrite development regulated by CREST, a calcium-regulated transcriptional activator. Science. 2004, 303 (5655): 197-202. 10.1126/science.1089845.

  10. 10.

    Azim E, Shnider SJ, Cederquist GY, Sohur US, Macklis JD: Lmo4 and Clim1 progressively delineate cortical projection neuron subtypes during development. Cerebral cortex. 2009, 19 (Suppl 1): i62-69.

  11. 11.

    Lai T, Jabaudon D, Molyneaux BJ, Azim E, Arlotta P, Menezes JR, Macklis JD: SOX5 controls the sequential generation of distinct corticofugal neuron subtypes. Neuron. 2008, 57 (2): 232-247. 10.1016/j.neuron.2007.12.023.

  12. 12.

    Molyneaux BJ, Arlotta P, Menezes JR, Macklis JD: Neuronal subtype specification in the cerebral cortex. Nat Rev Neurosci. 2007, 8 (6): 427-437.

  13. 13.

    Sun T, Collura RV, Ruvolo M, Walsh CA: Genomic and evolutionary analyses of asymmetrically expressed genes in human fetal left and right cerebral cortex. Cerebral cortex. 2006, 16 (Suppl 1): i18-25.

  14. 14.

    Ayoub AE, Oh S, Xie Y, Leng J, Cotney J, Dominguez MH, Noonan JP, Rakic P: Transcriptional programs in transient embryonic zones of the cerebral cortex defined by high-resolution mRNA sequencing. Proceedings of the National Academy of Sciences of the United States of America. 2011, 108 (36): 14950-14955. 10.1073/pnas.1112213108.

  15. 15.

    Kang HJ, Kawasawa YI, Cheng F, Zhu Y, Xu X, Li M, Sousa AM, Pletikos M, Meyer KA, Sedmak G, et al: Spatio-temporal transcriptome of the human brain. Nature. 2011, 478 (7370): 483-489. 10.1038/nature10523.

  16. 16.

    Hempel CM, Sugino K, Nelson SB: A manual method for the purification of fluorescently labeled neurons from the mammalian brain. Nature protocols. 2007, 2 (11): 2924-2929. 10.1038/nprot.2007.416.

  17. 17.

    Okaty BW, Miller MN, Sugino K, Hempel CM, Nelson SB: Transcriptional and electrophysiological maturation of neocortical fast-spiking GABAergic interneurons. J Neurosci. 2009, 29 (21): 7040-7052. 10.1523/JNEUROSCI.0105-09.2009.

  18. 18.

    Sugino K, Hempel CM, Miller MN, Hattox AM, Shapiro P, Wu C, Huang ZJ, Nelson SB: Molecular taxonomy of major neuronal classes in the adult mouse forebrain. Nature neuroscience. 2006, 9 (1): 99-107. 10.1038/nn1618.

  19. 19.

    Nadarajah B, Brunstrom JE, Grutzendler J, Wong RO, Pearlman AL: Two modes of radial migration in early development of the cerebral cortex. Nat Neurosci. 2001, 4 (2): 143-150. 10.1038/83967.

  20. 20.

    Arnold SJ, Huang GJ, Cheung AF, Era T, Nishikawa S, Bikoff EK, Molnar Z, Robertson EJ, Groszer M: The T-box transcription factor Eomes/Tbr2 regulates neurogenesis in the cortical subventricular zone. Genes & development. 2008, 22 (18): 2479-2484. 10.1101/gad.475408.

  21. 21.

    Baala L, Briault S, Etchevers HC, Laumonnier F, Natiq A, Amiel J, Boddaert N, Picard C, Sbiti A, Asermouh A, et al: Homozygous silencing of T-box transcription factor EOMES leads to microcephaly with polymicrogyria and corpus callosum agenesis. Nature genetics. 2007, 39 (4): 454-456. 10.1038/ng1993.

  22. 22.

    Nichols AJ, Olson EC: Reelin promotes neuronal orientation and dendritogenesis during preplate splitting. Cerebral cortex. 2010, 20 (9): 2213-2223. 10.1093/cercor/bhp303.

  23. 23.

    Visel A, Thaller C, Eichele G: GenePaint.org: an atlas of gene expression patterns in the mouse embryo. Nucleic acids research. 2004, 32 (Database issue): 552-556.

  24. 24.

    Carson JP, Ju T, Thaller C, Warren J, Bello M, Kakadiaris I, Chiu W, Eichele G: Automated characterization of gene expression patterns with an atlas of the mouse brain. Conference proceedings: Annual International Conference of the IEEE Engineering in Medicine and Biology Society IEEE Engineering in Medicine and Biology Society Conference. 2004, 4: 2917-2920.

  25. 25.

    Bulfone A, Martinez S, Marigo V, Campanella M, Basile A, Quaderi N, Gattuso C, Rubenstein JL, Ballabio A: Expression pattern of the Tbr2 (Eomesodermin) gene during mouse and chick brain development. Mech Dev. 1999, 84 (1–2): 133-138.

  26. 26.

    Kwon GS, Hadjantonakis AK: Eomes:: GFP-a tool for live imaging cells of the trophoblast, primitive streak, and telencephalon in the mouse embryo. Genesis. 2007, 45 (4): 208-217. 10.1002/dvg.20293.

  27. 27.

    Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, Braisted J, Klapa M, Currier T, Thiagarajan M, et al: TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 2003, 34 (2): 374-378.

  28. 28.

    Saeed AI, Bhagabati NK, Braisted JC, Liang W, Sharov V, Howe EA, Li J, Thiagarajan M, White JA, Quackenbush J: TM4 microarray software suite. Methods in enzymology. 2006, 411: 134-193.

  29. 29.

    Cameron DA, Gentile KL, Middleton FA, Yurco P: Gene expression profiles of intact and regenerating zebrafish retina. Molecular vision. 2005, 11: 775-791.

  30. 30.

    Nielsen KB, Kruhoffer M, Holm IE, Jorgensen AL, Nielsen AL: 1Identification of genes differentially expressed in the embryonic pig cerebral cortex before and after appearance of gyration. BMC Res Notes. 2010, 3: 127-10.1186/1756-0500-3-127.

  31. 31.

    Ramskold D, Wang ET, Burge CB, Sandberg R: An abundance of ubiquitously expressed genes revealed by tissue transcriptome sequence data. PLoS Comput Biol. 2009, 5 (12): e1000598-10.1371/journal.pcbi.1000598.

  32. 32.

    Rice DS, Sheldon M, D'Arcangelo G, Nakajima K, Goldowitz D, Curran T: Disabled-1 acts downstream of Reelin in a signaling pathway that controls laminar organization in the mammalian brain. Development. 1998, 125 (18): 3719-3729.

  33. 33.

    Angevine JB SR: Autoradiographic study of cell migration during histogenesis of cerebral cortex in the mouse. Nature. 1961, 192: 766-768.

  34. 34.

    Alvarez-Bolado G, Eichele G: Analysing the developing brain transcriptome with the GenePaint platform. The Journal of physiology. 2006, 575 (Pt 2): 347-352.

  35. 35.

    Kwon YT, Tsai LH: The role of the p35/cdk5 kinase in cortical development. Results Probl Cell Differ. 2000, 30: 241-253. 10.1007/978-3-540-48002-0_10.

  36. 36.

    Takei Y, Teng J, Harada A, Hirokawa N: Defects in axonal elongation and neuronal migration in mice with disrupted tau and map1b genes. J Cell Biol. 2000, 150 (5): 989-1000. 10.1083/jcb.150.5.989.

  37. 37.

    Yi JJ, Barnes AP, Hand R, Polleux F, Ehlers MD: TGF-beta signaling specifies axons during brain development. Cell. 2010, 142 (1): 144-157. 10.1016/j.cell.2010.06.010.

  38. 38.

    Chenn A, McConnell SK: Cleavage orientation and the asymmetric inheritance of Notch1 immunoreactivity in mammalian neurogenesis. Cell. 1995, 82 (4): 631-641. 10.1016/0092-8674(95)90035-7.

  39. 39.

    Noctor SC, Martinez-Cerdeno V, Ivic L, Kriegstein AR: Cortical neurons arise in symmetric and asymmetric division zones and migrate through specific phases. Nat Neurosci. 2004, 7 (2): 136-144. 10.1038/nn1172.

  40. 40.

    Pontious A, Kowalczyk T, Englund C, Hevner RF: Role of intermediate progenitor cells in cerebral cortex development. Developmental neuroscience. 2008, 30 (1–3): 24-32.

  41. 41.

    Wrobel CN, Mutch CA, Swaminathan S, Taketo MM, Chenn A: Persistent expression of stabilized beta-catenin delays maturation of radial glial cells into intermediate progenitors. Developmental Biology. 2007, 309 (2): 285-297. 10.1016/j.ydbio.2007.07.013.

  42. 42.

    Nieto M, Schuurmans C, Britz O, Guillemot F: Neural bHLH genes control the neuronal versus glial fate decision in cortical progenitors. Neuron. 2001, 29 (2): 401-413. 10.1016/S0896-6273(01)00214-8.

  43. 43.

    de Anda FC, Meletis K, Ge X, Rei D, Tsai LH: Centrosome motility is essential for initial axon formation in the neocortex. J Neurosci . 2010, 30 (31): 10391-10406. 10.1523/JNEUROSCI.0381-10.2010.

  44. 44.

    Micheva KD, Busse B, Weiler NC, O'Rourke N, Smith SJ: Single-synapse analysis of a diverse synapse population: proteomic imaging methods and markers. Neuron. 2010, 68 (4): 639-653. 10.1016/j.neuron.2010.09.024.

  45. 45.

    Hevner RF, Shi L, Justice N, Hsueh Y, Sheng M, Smiga S, Bulfone A, Goffinet AM, Campagnoni AT, Rubenstein JL: Tbr1 regulates differentiation of the preplate and layer 6. Neuron. 2001, 29 (2): 353-366. 10.1016/S0896-6273(01)00211-2.

  46. 46.

    Joshi PS, Molyneaux BJ, Feng L, Xie X, Macklis JD, Gan L: Bhlhb5 regulates the postmitotic acquisition of area identities in layers II-V of the developing neocortex. Neuron. 2008, 60 (2): 258-272. 10.1016/j.neuron.2008.08.006.

  47. 47.

    Kinameri E, Inoue T, Aruga J, Imayoshi I, Kageyama R, Shimogori T, Moore AW: Prdm proto-oncogene transcription factor family expression and interaction with the Notch-Hes pathway in mouse neurogenesis. PLoS ONE. 2008, 3 (12): e3859-10.1371/journal.pone.0003859.

  48. 48.

    Britanova O, de Juan Romero C, Cheung A, Kwan KY, Schwark M, Gyorgy A, Vogel T, Akopov S, Mitkovski M, Agoston D, et al: Satb2 is a postmitotic determinant for upper-layer neuron specification in the neocortex. Neuron. 2008, 57 (3): 378-392. 10.1016/j.neuron.2007.12.028.

  49. 49.

    Chen B, Wang SS, Hattox AM, Rayburn H, Nelson SB, McConnell SK: The Fezf2-Ctip2 genetic pathway regulates the fate choice of subcortical projection neurons in the developing cerebral cortex. Proceedings of the National Academy of Sciences of the United States of America. 2008, 105 (32): 11382-11387. 10.1073/pnas.0804918105.

  50. 50.

    Takagi T, Jin W, Taya K, Watanabe G, Mori K, Ishii S: Schnurri-2 mutant mice are hypersensitive to stress and hyperactive. Brain research. 2006, 1108 (1): 88-97. 10.1016/j.brainres.2006.06.018.

  51. 51.

    Gupta A, Sanada K, Miyamoto DT, Rovelstad S, Nadarajah B, Pearlman AL, Brunstrom J, Tsai LH: Layering defect in p35 deficiency is linked to improper neuronal-glial interaction in radial migration. Nat Neurosci. 2003, 6 (12): 1284-1291. 10.1038/nn1151.

  52. 52.

    Ohshima T, Hirasawa M, Tabata H, Mutoh T, Adachi T, Suzuki H, Saruta K, Iwasato T, Itohara S, Hashimoto M, et al: Cdk5 is required for multipolar-to-bipolar transition during radial neuronal migration and proper dendrite development of pyramidal neurons in the cerebral cortex. Development. 2007, 134 (12): 2273-2282. 10.1242/dev.02854.

  53. 53.

    Zorina Y, Iyengar R, Bromberg KD: Cannabinoid 1 receptor and interleukin-6 receptor together induce integration of protein kinase and transcription factor signaling to trigger neurite outgrowth. The Journal of biological chemistry. 2010, 285 (2): 1358-1370. 10.1074/jbc.M109.049841.

  54. 54.

    Nowakowska BA, Obersztyn E, Szymanska K, Bekiesinska-Figatowska M, Xia Z, Ricks CB, Bocian E, Stockton DW, Szczaluba K, Nawara M, et al: Severe mental retardation, seizures, and hypotonia due to deletions of MEF2C. Am J Med Genet B Neuropsychiatr Genet . 2010, 153B (5): 1042-1051.

  55. 55.

    Zoubovsky SP, Pogorelov VM, Taniguchi Y, Kim SH, Yoon P, Nwulia E, Sawa A, Pletnikov MV, Kamiya A: Working memory deficits in neuronal nitric oxide synthase knockout mice: potential impairments in prefrontal cortex mediated cognitive function. Biochemical and biophysical research communications. 2011, 408 (4): 707-712. 10.1016/j.bbrc.2011.04.097.

  56. 56.

    Hadjivassiliou G, Martinian L, Squier W, Blumcke I, Aronica E, Sisodiya SM, Thom M: The application of cortical layer markers in the evaluation of cortical dysplasias in epilepsy. Acta neuropathologica. 2010, 120 (4): 517-528. 10.1007/s00401-010-0686-x.

  57. 57.

    Kowalczyk T, Pontious A, Englund C, Daza RA, Bedogni F, Hodge R, Attardo A, Bell C, Huttner WB, Hevner RF: Intermediate neuronal progenitors (Basal Progenitors) produce pyramidal-projection neurons for all layers of cerebral cortex. Cereb Cortex. 2009, 19 (10): 2439-2450. 10.1093/cercor/bhn260.

  58. 58.

    Pierfelice T, Alberi L, Gaiano N: Notch in the vertebrate nervous system: an old dog with new tricks. Neuron. 2011, 69 (5): 840-855. 10.1016/j.neuron.2011.02.031.

  59. 59.

    Munji RN, Choe Y, Li G, Siegenthaler JA, Pleasure SJ: Wnt signaling regulates neuronal differentiation of cortical intermediate progenitors. J Neurosci. 2011, 31 (5): 1676-1687. 10.1523/JNEUROSCI.5404-10.2011.

  60. 60.

    de Almodovar R, Fabre PJ, Knevels E, Coulon C, Segura I, Haddick PC, Aerts L, Delattin N, Strasser G, Oh WJ, et al: VEGF mediates commissural axon chemoattraction through its receptor Flk1. Neuron. 2011, 70 (5): 966-978. 10.1016/j.neuron.2011.04.014.

  61. 61.

    Siegenthaler JA, Ashique AM, Zarbalis K, Patterson KP, Hecht JH, Kane MA, Folias AE, Choe Y, May SR, Kume T, et al: Retinoic acid from the meninges regulates cortical neuron generation. Cell. 2009, 139 (3): 597-609. 10.1016/j.cell.2009.10.004.

  62. 62.

    Woodhead GJ, Mutch CA, Olson EC, Chenn A: Cell-autonomous beta-catenin signaling regulates cortical precursor proliferation. J Neurosci. 2006, 26 (48): 12620-12630. 10.1523/JNEUROSCI.3180-06.2006.

  63. 63.

    Mutch CA, Schulte JD, Olson E, Chenn A: Beta-catenin signaling negatively regulates intermediate progenitor population numbers in the developing cortex. PLoS ONE. 2010, 5 (8): e12376-10.1371/journal.pone.0012376.

  64. 64.

    Gaiano N, Fishell G: The role of notch in promoting glial and neural stem cell fates. Annu Rev Neurosci. 2002, 25: 471-490. 10.1146/annurev.neuro.25.030702.130823.

  65. 65.

    Gaiano N, Nye JS, Fishell G: Radial glial identity is promoted by Notch1 signaling in the murine forebrain. Neuron. 2000, 26 (2): 395-404. 10.1016/S0896-6273(00)81172-1.

  66. 66.

    Berezovska O, Frosch M, McLean P, Knowles R, Koo E, Kang D, Shen J, Lu FM, Lux SE, Tonegawa S, et al: The Alzheimer-related gene presenilin 1 facilitates notch 1 in primary mammalian neurons. Brain research Molecular brain research. 1999, 69 (2): 273-280. 10.1016/S0169-328X(99)00119-9.

  67. 67.

    Sestan N, Artavanis-Tsakonas S, Rakic P: Contact-dependent inhibition of cortical neurite growth mediated by notch signaling. Science. 1999, 286 (5440): 741-746. 10.1126/science.286.5440.741.

  68. 68.

    Stevens HE, Smith KM, Maragnoli ME, Fagel D, Borok E, Shanabrough M, Horvath TL, Vaccarino FM: Fgfr2 is required for the development of the medial prefrontal cortex and its connections with limbic circuits. J Neurosci. 2010, 30 (16): 5590-5602. 10.1523/JNEUROSCI.5837-09.2010.

  69. 69.

    Moldrich RX, Mezzera C, Holmes WM, Goda S, Brookfield SJ, Rankin AJ, Barr E, Kurniawan N, Dewar D, Richards LJ, et al: Fgfr3 regulates development of the caudal telencephalon. Dev Dyn. 2011, 240 (6): 1586-1599. 10.1002/dvdy.22636.

  70. 70.

    Polleux F, Giger RJ, Ginty DD, Kolodkin AL, Ghosh A: Patterning of cortical efferent projections by semaphorin-neuropilin interactions. Science. 1998, 282 (5395): 1904-1906.

  71. 71.

    Metin C, Deleglise D, Serafini T, Kennedy TE, Tessier-Lavigne M: A role for netrin-1 in the guidance of cortical efferents. Development. 1997, 124 (24): 5063-5074.

  72. 72.

    Lopez-Bendito G, Flames N, Ma L, Fouquet C, Di Meglio T, Chedotal A, Tessier-Lavigne M, Marin O: Robo1 and Robo2 cooperate to control the guidance of major axonal tracts in the mammalian forebrain. J Neurosci. 2007, 27 (13): 3395-3407. 10.1523/JNEUROSCI.4605-06.2007.

  73. 73.

    Chen G, Sima J, Jin M, Wang KY, Xue XJ, Zheng W, Ding YQ, Yuan XB: Semaphorin-3A guides radial migration of cortical neurons during development. Nat Neurosci. 2008, 11 (1): 36-44. 10.1038/nn2018.

  74. 74.

    Polleux F, Morrow T, Ghosh A: Semaphorin 3A is a chemoattractant for cortical apical dendrites. Nature. 2000, 404 (6778): 567-573. 10.1038/35007001.

  75. 75.

    Fenstermaker V, Chen Y, Ghosh A, Yuste R: Regulation of dendritic length and branching by semaphorin 3A. J Neurobiol. 2004, 58 (3): 403-412. 10.1002/neu.10304.

  76. 76.

    Shelly M, Cancedda L, Lim BK, Popescu AT, Cheng PL, Gao H, Poo MM: Semaphorin3A regulates neuronal polarization by suppressing axon formation and promoting dendrite growth. Neuron. 2011, 71 (3): 433-446. 10.1016/j.neuron.2011.06.041.

  77. 77.

    Takemoto M, Hattori Y, Zhao H, Sato H, Tamada A, Sasaki S, Nakajima K, Yamamoto N: Laminar and areal expression of unc5d and its role in cortical cell survival. Cerebral cortex. 2011, 21 (8): 1925-1934. 10.1093/cercor/bhq265.

  78. 78.

    Argaw A, Duff G, Zabouri N, Cecyre B, Chaine N, Cherif H, Tea N, Lutz B, Ptito M, Bouchard JF: Concerted action of CB1 cannabinoid receptor and deleted in colorectal cancer in axon guidance. J Neurosci. 2011, 31 (4): 1489-1499. 10.1523/JNEUROSCI.4134-09.2011.

  79. 79.

    Verty AN, McFarlane JR, McGregor IS, Mallet PE: Evidence for an interaction between CB1 cannabinoid and melanocortin MCR-4 receptors in regulating food intake. Endocrinology. 2004, 145 (7): 3224-3231. 10.1210/en.2004-0059.

  80. 80.

    Kieffer TJ: GIP or not GIP? That is the question. Trends in pharmacological sciences. 2003, 24 (3): 110-112. 10.1016/S0165-6147(03)00031-2.

  81. 81.

    Duttaroy A, Zimliki CL, Gautam D, Cui Y, Mears D, Wess J: Muscarinic stimulation of pancreatic insulin and glucagon release is abolished in m3 muscarinic acetylcholine receptor-deficient mice. Diabetes. 2004, 53 (7): 1714-1720. 10.2337/diabetes.53.7.1714.

  82. 82.

    Ohno-Shosaku T, Matsui M, Fukudome Y, Shosaku J, Tsubokawa H, Taketo MM, Manabe T, Kano M: Postsynaptic M1 and M3 receptors are responsible for the muscarinic enhancement of retrograde endocannabinoid signalling in the hippocampus. The European journal of neuroscience. 2003, 18 (1): 109-116. 10.1046/j.1460-9568.2003.02732.x.

  83. 83.

    Matsumoto M, Straub RE, Marenco S, Nicodemus KK, Matsumoto S, Fujikawa A, Miyoshi S, Shobo M, Takahashi S, Yarimizu J, et al: The evolutionarily conserved G protein-coupled receptor SREB2/GPR85 influences brain size, behavior, and vulnerability to schizophrenia. Proceedings of the National Academy of Sciences of the United States of America. 2008, 105 (16): 6133-6138. 10.1073/pnas.0710717105.

  84. 84.

    Zhou L, Tissir F, Goffinet AM: The atypical cadherin Celsr3 regulates the development of the axonal blueprint. Novartis Foundation symposium. 2007, 288: 130-134. discussion 134-140, 276-181

  85. 85.

    Conhaim J, Easton CR, Becker MI, Barahimi M, Cedarbaum ER, Moore JG, Mather LF, Dabagh S, Minter DJ, Moen SP, et al: Developmental changes in propagation patterns and transmitter dependence of waves of spontaneous activity in the mouse cerebral cortex. The Journal of physiology. 2011, 589 (Pt 10): 2529-2541.

  86. 86.

    Orr-Urtreger A, Goldner FM, Saeki M, Lorenzo I, Goldberg L, De Biasi M, Dani JA, Patrick JW, Beaudet AL: Mice deficient in the alpha7 neuronal nicotinic acetylcholine receptor lack alpha-bungarotoxin binding sites and hippocampal fast nicotinic currents. J Neurosci. 1997, 17 (23): 9165-9171.

  87. 87.

    Shinawi M, Schaaf CP, Bhatt SS, Xia Z, Patel A, Cheung SW, Lanpher B, Nagl S, Herding HS, Nevinny-Stickel C, et al: A small recurrent deletion within 15q13.3 isassociated with a range of neurodevelopmental phenotypes. Nature genetics. 2009, 41 (12): 1269-1271. 10.1038/ng.481.

  88. 88.

    McCabe AK, Easton CR, Lischalk JW, Moody WJ: Roles of glutamate and GABA receptors in setting the developmental timing of spontaneous synchronized activity in the developing mouse cortex. Dev Neurobiol. 2007, 67 (12): 1574-1588. 10.1002/dneu.20533.

  89. 89.

    Katz LC, Shatz CJ: Synaptic activity and the construction of cortical circuits. Science. 1996, 274 (5290): 1133-1138. 10.1126/science.274.5290.1133.

  90. 90.

    Nase G, Weishaupt J, Stern P, Singer W, Monyer H: Genetic and epigenetic regulation of NMDA receptor expression in the rat visual cortex. The European journal of neuroscience. 1999, 11 (12): 4320-4326. 10.1046/j.1460-9568.1999.00859.x.

  91. 91.

    Marini C, Scheffer IE, Nabbout R, Suls A, De Jonghe P, Zara F, Guerrini R: The genetics of Dravet syndrome. Epilepsia. 2011, 52 (Suppl 2): 24-29.

  92. 92.

    Bosman LW, Heinen K, Spijker S, Brussaard AB: Mice lacking the major adult GABAA receptor subtype have normal number of synapses, but retain juvenile IPSC kinetics until adulthood. Journal of neurophysiology. 2005, 94 (1): 338-346. 10.1152/jn.00084.2005.

  93. 93.

    Ben-Ari Y, Khazipov R, Leinekugel X, Caillard O, Gaiarsa JL: GABAA, NMDA and AMPA receptors: a developmentally regulated 'menage a trois'. Trends in neurosciences. 1997, 20 (11): 523-529. 10.1016/S0166-2236(97)01147-8.

  94. 94.

    Rheims S, Minlebaev M, Ivanov A, Represa A, Khazipov R, Holmes GL, Ben-Ari Y, Zilberter Y: Excitatory GABA in rodent developing neocortex in vitro. Journal of neurophysiology. 2008, 100 (2): 609-619. 10.1152/jn.90402.2008.

  95. 95.

    Young-Pearse TL, Ivic L, Kriegstein AR, Cepko CL: Characterization of mice with targeted deletion of glycine receptor alpha 2. Molecular and cellular biology. 2006, 26 (15): 5728-5734. 10.1128/MCB.00237-06.

  96. 96.

    Konig N, Roch G, Marty R: The onset of synaptogenesis in rat temporal cortex. Anatomy and embryology. 1975, 148 (1): 73-87. 10.1007/BF00315564.

  97. 97.

    Welker C, Woolsey TA: Structure of layer IV in the somatosensory neocortex of the rat: description and comparison with the mouse. The Journal of comparative neurology. 1974, 158 (4): 437-453. 10.1002/cne.901580405.

  98. 98.

    Johnson SL, Franz C, Kuhn S, Furness DN, Ruttiger L, Munkner S, Rivolta MN, Seward EP, Herschman HR, Engel J, et al: Synaptotagmin IV determines the linear Ca2+ dependence of vesicle fusion at auditory ribbon synapses. Nature neuroscience. 2010, 13 (1): 45-52. 10.1038/nn.2456.

  99. 99.

    Nowack A, Yao J, Custer KL, Bajjalieh SM: SV2 regulates neurotransmitter release via multiple mechanisms. American journal of physiology Cell physiology. 2010, 299 (5): C960-967. 10.1152/ajpcell.00259.2010.

  100. 100.

    Linhoff MW, Lauren J, Cassidy RM, Dobie FA, Takahashi H, Nygaard HB, Airaksinen MS, Strittmatter SM, Craig AM: An unbiased expression screen for synaptogenic proteins identifies the LRRTM protein family as synaptic organizers. Neuron. 2009, 61 (5): 734-749. 10.1016/j.neuron.2009.01.017.

  101. 101.

    Spitzer NC: Spontaneous activity: functions of calcium transients in neuronal differentiation. Perspectives on developmental neurobiology. 1995, 2 (4): 379-386.

  102. 102.

    Spitzer NC: Electrical activity in early neuronal development. Nature. 2006, 444 (7120): 707-712. 10.1038/nature05300.

  103. 103.

    Picken Bahrey HL, Moody WJ: Early development of voltage-gated ion currents and firing properties in neurons of the mouse cerebral cortex. Journal of neurophysiology. 2003, 89 (4): 1761-1773.

  104. 104.

    Corlew R, Bosma MM, Moody WJ: Spontaneous, synchronous electrical activity in neonatal mouse cortical neurones. The Journal of physiology. 2004, 560 (Pt 2): 377-390.

  105. 105.

    Scholl UI, Choi M, Liu T, Ramaekers VT, Hausler MG, Grimmer J, Tobe SW, Farhi A, Nelson-Williams C, Lifton RP: Seizures, sensorineural deafness, ataxia, mental retardation, and electrolyte imbalance (SeSAME syndrome) caused by mutations in KCNJ10. Proceedings of the National Academy of Sciences of the United States of America. 2009, 106 (14): 5842-5847. 10.1073/pnas.0901749106.

  106. 106.

    Neusch C, Rozengurt N, Jacobs RE, Lester HA, Kofuji P: Kir4.1 Potassium channel subunit is crucial for oligodendrocyte development and in vivo myelination. J Neurosci. 2001, 21 (15): 5429-5438.

  107. 107.

    Schmid RS, Maness PF: L1 and NCAM adhesion molecules as signaling coreceptors in neuronal migration and process outgrowth. Current opinion in neurobiology. 2008, 18 (3): 245-250. 10.1016/j.conb.2008.07.015.

  108. 108.

    Liu G, Li W, Wang L, Kar A, Guan KL, Rao Y, Wu JY: DSCAM functions as a netrin receptor in commissural axon pathfinding. Proceedings of the National Academy of Sciences of the United States of America. 2009, 106 (8): 2951-2956. 10.1073/pnas.0811083106.

  109. 109.

    Park J, Liu B, Chen T, Li H, Hu X, Gao J, Zhu Y, Zhu Q, Qiang B, Yuan J, et al: Disruption of Nectin-like 1 cell adhesion molecule leads to delayed axonal myelination in the CNS. J Neurosci. 2008, 28 (48): 12815-12819. 10.1523/JNEUROSCI.2665-08.2008.

  110. 110.

    Thaxton C, Pillai AM, Pribisko AL, Labasque M, Dupree JL, Faivre-Sarrailh C, Bhat MA: In vivo deletion of immunoglobulin domains 5 and 6 in neurofascin (Nfasc) reveals domain-specific requirements in myelinated axons. J Neurosci. 2010, 30 (14): 4868-4876. 10.1523/JNEUROSCI.5951-09.2010.

  111. 111.

    Denaxa M, Chan CH, Schachner M, Parnavelas JG, Karagogeos D: The adhesion molecule TAG-1 mediates the migration of cortical interneurons from the ganglionic eminence along the corticofugal fiber system. Development. 2001, 128 (22): 4635-4644.

  112. 112.

    Denaxa M, Kyriakopoulou K, Theodorakis K, Trichas G, Vidaki M, Takeda Y, Watanabe K, Karagogeos D: The adhesion molecule TAG-1 is required for proper migration of the superficial migratory stream in the medulla but not of cortical interneurons. Dev Biol. 2005, 288 (1): 87-99. 10.1016/j.ydbio.2005.09.021.

  113. 113.

    Savvaki M, Panagiotaropoulos T, Stamatakis A, Sargiannidou I, Karatzioula P, Watanabe K, Stylianopoulou F, Karagogeos D, Kleopa KA: Impairment of learning and memory in TAG-1 deficient mice associated with shorter CNS internodes and disrupted juxtaparanodes. Molecular and cellular neurosciences. 2008, 39 (3): 478-490. 10.1016/j.mcn.2008.07.025.

  114. 114.

    McCarty JH, Lacy-Hulbert A, Charest A, Bronson RT, Crowley D, Housman D, Savill J, Roes J, Hynes RO: Selective ablation of alphav integrins in the central nervous system leads to cerebral hemorrhage, seizures, axonal degeneration and premature death. Development. 2005, 132 (1): 165-176.

  115. 115.

    Su H, Kim H, Pawlikowska L, Kitamura H, Shen F, Cambier S, Markovics J, Lawton MT, Sidney S, Bollen AW, et al: Reduced expression of integrin alphavbeta8 is associated with brain arteriovenous malformation pathogenesis. The American journal of pathology. 2010, 176 (2): 1018-1027. 10.2353/ajpath.2010.090453.

  116. 116.

    Mochida GH, Ganesh VS, Felie JM, Gleason D, Hill RS, Clapham KR, Rakiec D, Tan WH, Akawi N, Al-Saffar M, et al: A homozygous mutation in the tight-junction protein JAM3 causes hemorrhagic destruction of the brain, subependymal calcification, and congenital cataracts. American journal of human genetics. 2010, 87 (6): 882-889. 10.1016/j.ajhg.2010.10.026.

  117. 117.

    Inoue T, Tanaka T, Takeichi M, Chisaka O, Nakamura S, Osumi N: Role of cadherins in maintaining the compartment boundary between the cortex and striatum during development. Development. 2001, 128 (4): 561-569.

  118. 118.

    Bae GU, Domene S, Roessler E, Schachter K, Kang JS, Muenke M, Krauss RS: Mutations in CDON, encoding a hedgehog receptor, result in holoprosencephaly and defective interactions with other hedgehog receptors. American journal of human genetics. 2011, 89 (2): 231-240. 10.1016/j.ajhg.2011.07.001.

  119. 119.

    Zhang W, Hong M, Bae GU, Kang JS, Krauss RS: Boc modifies the holoprosencephaly spectrum of Cdo mutant mice. Dis Model Mech. 2011, 4 (3): 368-380. 10.1242/dmm.005744.

  120. 120.

    Takahashi T, Goto T, Miyama S, Nowakowski RS, Caviness VS: Sequence of neuron origin and neocortical laminar fate: relation to cell cycle of origin in the developing murine cerebral wall. J Neurosci. 1999, 19 (23): 10357-10371.

  121. 121.

    Agathocleous M, Harris WA: From progenitors to differentiated cells in the vertebrate retina. Annual review of cell and developmental biology. 2009, 25: 45-69. 10.1146/annurev.cellbio.042308.113259.

  122. 122.

    Agathocleous M, Iordanova I, Willardsen MI, Xue XY, Vetter ML, Harris WA, Moore KB: A directional Wnt/beta-catenin-Sox2-proneural pathway regulates the transition from proliferation to differentiation in the Xenopus retina. Development. 2009, 136 (19): 3289-3299. 10.1242/dev.040451.

  123. 123.

    Marek KW, Kurtz LM, Spitzer NC: cJun integrates calcium activity and tlx3 expression to regulate neurotransmitter specification. Nature neuroscience. 2010, 13 (8): 944-950. 10.1038/nn.2582.

  124. 124.

    Demarque M, Spitzer NC: Neurotransmitter phenotype plasticity: an unexpected mechanism in the toolbox of network activity homeostasis. Dev Neurobiol. 2012, 72 (1): 22-32. 10.1002/dneu.20909.

  125. 125.

    Holt CE, Bertsch TW, Ellis HM, Harris WA: Cellular determination in the Xenopus retina is independent of lineage and birth date. Neuron. 1988, 1 (1): 15-26. 10.1016/0896-6273(88)90205-X.

  126. 126.

    Turner DL, Snyder EY, Cepko CL: Lineage-independent determination of cell type in the embryonic mouse retina. Neuron. 1990, 4 (6): 833-845. 10.1016/0896-6273(90)90136-4.

  127. 127.

    Wetts R, Fraser SE: Multipotent precursors can give rise to all major cell types of the frog retina. Science. 1988, 239 (4844): 1142-1145. 10.1126/science.2449732.

  128. 128.

    Mione MC, Danevic C, Boardman P, Harris B, Parnavelas JG: Lineage analysis reveals neurotransmitter (GABA or glutamate) but not calcium-binding protein homogeneity in clonally related cortical neurons. J Neurosci. 1994, 14 (1): 107-123.

  129. 129.

    Walsh C, Reid C: Cell lineage and patterns of migration in the developing cortex. Ciba Found Symp. 1995, 193: 21-40. discussion 59-70

  130. 130.

    Heng JI, Nguyen L, Castro DS, Zimmer C, Wildner H, Armant O, Skowronska-Krawczyk D, Bedogni F, Matter JM, Hevner R, et al: Neurogenin 2 controls cortical neuron migration through regulation of Rnd2. Nature. 2008, 455 (7209): 114-118. 10.1038/nature07198.

  131. 131.

    Pacary E, Heng J, Azzarelli R, Riou P, Castro D, Lebel-Potter M, Parras C, Bell DM, Ridley AJ, Parsons M, et al: Proneural transcription factors regulate different steps of cortical neuron migration through Rnd-mediated inhibition of RhoA signaling. Neuron. 2011, 69 (6): 1069-1084. 10.1016/j.neuron.2011.02.018.

  132. 132.

    Molyneaux BJ, Arlotta P, Hirata T, Hibi M, Macklis JD: Fezl is required for the birth and specification of corticospinal motor neurons. Neuron. 2005, 47 (6): 817-831. 10.1016/j.neuron.2005.08.030.

  133. 133.

    Fame RM, MacDonald JL, Macklis JD: Development, specification, and diversity of callosal projection neurons. Trends in neurosciences. 2011, 34 (1): 41-50. 10.1016/j.tins.2010.10.002.

  134. 134.

    Bayer SA, Altman J, Russo RJ, Dai XF, Simmons JA: Cell migration in the rat embryonic neocortex. J Comp Neurol. 1991, 307 (3): 499-516. 10.1002/cne.903070312.

  135. 135.

    Olson EC, Kim S, Walsh CA: Impaired neuronal positioning and dendritogenesis in the neocortex after cell-autonomous Dab1 suppression. J Neurosci. 2006, 26 (6): 1767-1775. 10.1523/JNEUROSCI.3000-05.2006.

  136. 136.

    Kathleen Baxter K, Uittenbogaard M, Yoon J, Chiaramello A: The neurogenic basic helix-loop-helix transcription factor NeuroD6 concomitantly increases mitochondrial mass and regulates cytoskeletal organization in the early stages of neuronal differentiation. ASN Neuro. 2009, 1 (4): e00016.

  137. 137.

    Ackman JB, LoTurco JJ: The potential of endogenous neuronal replacement in developing cerebral cortex following hypoxic injury. Exp Neurol. 2006, 199 (1): 5-9. 10.1016/j.expneurol.2006.03.007.

Download references

Acknowledgements

The authors thank the following: Dr. Nick Gonchoroff, Director of Flow Cytometery at SUNY Upstate for FAC sorting, Karen Gentile, Technical and Adminstrative Director of the SUNY Microarray Core Facility for Affymetrix profiling and Dr. Robert Quinn and the staff in the Department of Laboratory Animal Resources for animal care.

Funding

This work was supported by the Department of Neuroscience and Physiology, the Hendricks Foundation, and a grant from the NINDS (NS066071) to ECO.

Author information

Affiliations

  1. Department of Neuroscience and Physiology, SUNY Upstate Medical University, Syracuse, NY, 13210, USA
    • David A Cameron
    • , Frank A Middleton
    •  & Eric C Olson
  2. Department of Electrical Engineering and Computer Science, Syracuse University, Syracuse, NY, 13210, USA
    • Eric C Olson
  3. Department of Pathology, Feinberg School of Medicine, Northwestern University, Chicago, IL, 60611, USA
    • Anjen Chenn
  4. Department of Neuroscience and Physiology, SUNY Upstate Medical University, 3295 Weiskotten Hall, 750 E, Adams St., Syracuse, NY, 13210, USA
    • Eric C Olson
Authors
  1. David A Cameron
    View author publications
    You can also search for this author in
  2. Frank A Middleton
    View author publications
    You can also search for this author in
  3. Anjen Chenn
    View author publications
    You can also search for this author in
  4. Eric C Olson
    View author publications
    You can also search for this author in

Corresponding author

Correspondence to Eric C Olson.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

EO and FM designed the study. EO performed the study and co-authored the manuscript with DC and AC. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Comparison of QISPs from repeated measurement of Dab1 in situ expression in rostral neocortex along ( a, b ) the rostral-caudal axis and (c, d) the lateral-medial axis. Panels b and d plot strength of the in situ hybridization label (“intensity”) as a function of cortical position. All in situ hybridizations are derived from Genepaint. (PDF )

Additional file 2: Table S1. RMA values for eye-specific genes that were used to establish not-expressed threshold for this study. (XLSX 55 KB)

Additional file 3: Table S2. Complete list of all ≥ 3 fold up-regulated genes in E14.5 GFP + Eomes lineage compared to E13.5 GFP + Eomes lineage. (XLSX 58 KB)

Additional file 4: Table S3. Complete list of all ≥ 3 fold down-regulated genes in E14.5 GFP + Eomes lineage compared to E13.5 GFP + Eomes lineage. (XLSX 89 KB)

Additional file 5: Table S4. Complete list of all ≥ 3 fold up-regulated genes in the E14.5 GFP + Eomes lineage compared to E13.5 GFP- precursors. Blue cells indicate genes that were unavailable or could not be analyzed by Genepaint in situ. (XLSX 66 KB)

Additional file 6: Table S5. Complete list of all ≥ 3 fold down-regulated genes in the E14.5 GFP + Eomes lineage compared to GFP- precursors. (XLSX )

Additional file 7: Functional groupings of QISPs based on Gene Ontology (GO) molecular function and biological process annotations. Genes are grouped according to zone of expression: deep (D, dark gray), middle (M, light gray), upper (U, white). Listed genes are ranked by fold up-regulation. (PDF 194 KB)

Authors’ original submitted files for images

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Cameron, D.A., Middleton, F.A., Chenn, A. et al. Hierarchical clustering of gene expression patterns in the Eomes + lineage of excitatory neurons during early neocortical development. BMC Neurosci 13, 90 (2012). https://doi.org/10.1186/1471-2202-13-90

Download citation

Keywords

  • Transcriptome
  • Migration
  • Cortex
  • Profiling
  • Excitatory neuron

BMC Neuroscience

ISSN: 1471-2202

Contact us