Molecular networks involved in mouse cerebral corticogenesis and spatio-temporal regulation of Sox4 and Sox11 novel antisense transcripts revealed by transcriptome profiling
© Ling et al.; licensee BioMed Central Ltd. 2009
Received: 6 June 2009
Accepted: 2 October 2009
Published: 2 October 2009
Development of the cerebral cortex requires highly specific spatio-temporal regulation of gene expression. It is proposed that transcriptome profiling of the cerebral cortex at various developmental time points or regions will reveal candidate genes and associated molecular pathways involved in cerebral corticogenesis.
Serial analysis of gene expression (SAGE) libraries were constructed from C57BL/6 mouse cerebral cortices of age embryonic day (E) 15.5, E17.5, postnatal day (P) 1.5 and 4 to 6 months. Hierarchical clustering analysis of 561 differentially expressed transcripts showed regionalized, stage-specific and co-regulated expression profiles. SAGE expression profiles of 70 differentially expressed transcripts were validated using quantitative RT-PCR assays. Ingenuity pathway analyses of validated differentially expressed transcripts demonstrated that these transcripts possess distinctive functional properties related to various stages of cerebral corticogenesis and human neurological disorders. Genomic clustering analysis of the differentially expressed transcripts identified two highly transcribed genomic loci, Sox4 and Sox11, during embryonic cerebral corticogenesis. These loci feature unusual overlapping sense and antisense transcripts with alternative polyadenylation sites and differential expression. The Sox4 and Sox11 antisense transcripts were highly expressed in the brain compared to other mouse organs and are differentially expressed in both the proliferating and differentiating neural stem/progenitor cells and P19 (embryonal carcinoma) cells.
We report validated gene expression profiles that have implications for understanding the associations between differentially expressed transcripts, novel targets and related disorders pertaining to cerebral corticogenesis. The study reports, for the first time, spatio-temporally regulated Sox4 and Sox11 antisense transcripts in the brain, neural stem/progenitor cells and P19 cells, suggesting they have an important role in cerebral corticogenesis and neuronal/glial cell differentiation.
Complex behavioral tasks, from perception of sensory input and the control of motor output to cognitive functions such as learning and memory, are dependent on the precise development of innumerable interconnections of neuronal networks in the cerebral cortex. The development of the cerebral cortex (also known as cerebral corticogenesis) involves the specific influence of intrinsic and extrinsic mechanisms, which are triggered spatio-temporally [1–4]. Between embryonic day 11 (E11) and 18 (E18), the mouse cerebral cortex develops from a relatively homogenous band of mitotic multipotent progenitor cells into a complex laminated structure containing various classes of neuronal cells [2, 4–6]. Cerebral corticogenesis involves: proliferation of multipotent progenitors (E11 to E16.5); migration of postmitotic cells (E11 to E17); cell morphogenesis (E13 to E18); gliogenesis and synaptogenesis (E16 until early postnatal period); and reorganization, elimination and stabilization of neuronal networks (up to adulthood).
The mouse cerebral cortex develops in the latero-medial and rostro-caudal axes [7, 8]. At E11, the primordial plexiform layer begins to form in the most lateral part of the neural wall. Its growth continues in the latero-medial axis to the medial part of the telencephalon by E13. The primordial plexiform layer is also expanded in the rostro-caudal axis. The growth in this axis is always less than the growth in the latero-medial axis. The first wave of migratory neuronal cells form the cortical plate 2 days later after the development of the primordial plexiform layer. These events are followed by the development of the cortical plate into an organization of six distinct layers that forms the adult cerebral cortex. Generally, the rostral-most regions of the adult cerebral cortex consist of areas involved in executive functions and motor coordination, whilst the caudal-most regions consist of areas involved in sensory perception such as visual function. Although distinct functional arealization of the cerebral cortex do not fully apply to rodents, mounting evidence suggesting that regulated arealization exists has been shown in mice involving transcription factors such as empty spiracles homolog 2 (Emx2; Drosophila) [9, 10], paired box gene 6 (Pax6) , COUP transcription factor 1 (Coup-tf1) , Sp8 transcription factor (Sp8) [12, 13], distal-less homeobox 1/2 (Dlx1/2) and gastrulation brain homeobox 2 (Gbx2) . The extensive cyto-architectural and anatomical changes occurring in a spatio-temporal manner during the peak (E15) and at the end (E17) of embryonic cerebral neurogenesis as well as during early postnatal (P1) corticogenesis through to adulthood involves complex underlying molecular regulatory networks.
Complex molecular regulatory elements are important determinants in both spatial and temporal cerebral corticogenesis. These elements regulate gene expression at the chromatin, DNA or RNA, and protein levels through chromatin packaging or remodeling, histone acetylation and deacetylation, chromatin insulation, DNA methylation, post-transcriptional regulation and post-translational protein modification or degradation signaling pathways [15–17]. Other processes involved in such regulation include pre-mRNA processing and nuclear mRNA retention by nuclear-specific paraspeckle complexes [18, 19], microRNAs (miRNAs) that interfere with mRNA translation [20–22], and a new class of under-characterized non-coding RNA transcripts known as natural antisense transcripts (NATs) [23, 24]. These regulatory networks play a pivotal role in establishing when, where and how a multipotent progenitor cell should proliferate, migrate and settle at a designated position in the cortex. The information regarding regulatory networks during cerebral corticogenesis, however, remains incomplete and does not provide a comprehensive view of the underlying regulatory elements throughout this complex event.
In this study, we employed both short and long 3' serial analysis of gene expression (SAGE) technologies [25, 26] to identify differentially expressed regulatory elements by comparing transcriptomes of cerebral cortices generated from four selected developmental stages: E15.5, E17.5, P1.5 and 4 to 6 months old. We also compared rostral to caudal regions of E15.5 and left to right regions of adult cerebral cortices. We report temporally co-regulated gene clusters, novel molecular networks and associated pathways, novel candidates in regionalized development and genomic clustering of SRY-box containing gene 4 (Sox4) and SRY-box containing gene 11 (Sox11) sense and antisense transcripts. The genomic clustering analysis led us to the discovery of spatio-temporal regulation of novel Sox4 and Sox11 antisense transcripts as well as differential regulation of these transcripts in proliferating and differentiating neural stem/progenitor cells (NSPCs) and P19 (embryonal carcinoma) cells.
Generation and analysis of SAGE tags
SAGE library information
Generated SAGE tag (library)
Unique tags (after scaling to 100,000 tags/library)
GEO accession number
Rostral cerebral cortex
Caudal cerebral cortex
Left cerebral cortex
Left cerebral cortex
Left cerebral cortex
Left cerebral cortex
Right cerebral cortex
Right cerebral cortex
Analysis of differentially expressed transcripts/tags
Summary of tag classification into various categories and comparisons
DETs in various comparisons*
Number of DETs†
E15.5 R versus C
Adult L versus Ri
E15.5 versus E17.5
E15.5 versus P1.5
P1.5 versus Ad
E versus Ad
Approximately 69% of DETs have been assigned to known genes and 6% to expressed sequence tags (ESTs), while 10% of DETs have multiple matches (tag matching multiple gene identifiers) and 8% having ambiguous matches (tag matching the same gene identifier at multiple chromosome loci). Approximately 8% of DETs have no matches and it is most likely that these DETs belong to transcripts from less defined regions, such as centromeric and telomeric areas of a chromosome or assembly gaps of the mouse genome.
Hierarchical clustering of DETs
We performed a systematic gene ontology functional clustering using the standardized Gene Ontology term analysis tool DAVID (Database for Annotation, Visualization and Integrated Discovery)  (Additional data file 3). Functional analysis of these gene clusters reveals that they have distinctive roles during cerebral corticogenesis. Embryonic-specific gene clusters (1, 5 and 6) are dominated by genes associated with cellular protein and macromolecule metabolic processes or biosynthesis, and nervous system and neuron development. These results match with the expected functional ontologies during embryonic cerebral cortex development in which neuronal migration, differentiation and axonogenesis events are at their peak. In contrast, adult-specific gene clusters (4, 8 and 10) consist of genes related to biological processes in the mature cerebral cortex, such as ion homeostasis, synaptic transmission and regulation of neurotransmitter level. In addition, these gene clusters are also enriched for ribonucleotide/nucleotide binding activity and components of cytoplasmic membrane-bound vesicles. These functional ontologies are in accordance with adult cerebral cortex function, which features synaptogenesis and nerve impulse transmission at synapses. Gene-switching clusters 2 and 7 are enriched with gene ontologies that are similar to both the embryonic- and adult-specific gene clusters. In addition, these gene clusters are also enriched for microtubule cytoskeleton organization and biogenesis, nucleotide biosynthesis and regulation of mRNA translation processes.
Quantitative RT-PCR validation of DETs and gene clusters
RT-qPCR validation of SAGE profile for E versus adult comparison
calcium/calmodulin-dependent protein kinase II alpha, Camk2a
early growth response 1, Egr1
proteolipid protein (myelin) 1, Plp1
calcium/calmodulin-dependent protein kinase II inhibitor 1, Camk2n1
crystallin, alpha B, Cryab
EST from adult C57BL/6 subfornical organ and postrema tissues
septin 4, Sept4
EST from 27-32 days C57BL/6 hippocampus tissue
synuclein, beta, Sncb
neuronal pentraxin receptor, Nptxr
protein phosphatase 1, regulatory (inhibitor) subunit 1B, Ppp1r1b
glutathione S-transferase, mu 1, Gstm1
chromogranin B, Chgb
inositol 1,4,5-triphosphate receptor 1, Itpr1
regulator of G-protein signaling 4, Rgs4
Mus musculus 10 days neonate cortex cDNA
tetraspanin 7, Tspan7
signal-regulatory protein alpha, Sirpa
EST from mouse brain
protein phosphatase 3, catalytic subunit, alpha isoform, Ppp3ca
amyloid beta (A4) precursor protein, App
hypoxanthine guanine phosphoribosyl transferase 1, Hprt1
glutamate receptor, ionotropic, AMPA3 (alpha 3), Gria3
calmodulin 1, Calm1
cathepsin D, Ctsd
myelin basic protein, Mbp
phosphatidic acid phosphatase type 2B, Ppap2b
NADH dehydrogenase 6, mitochondrial, mt-ND6
Mus musculus adult male corpora quadrigemina cDNA
chloride channel 2, Clcn2
Mus musculus NOD-derived CD11c +ve dendritic cells cDNA
Kruppel-like factor 3 (basic), Klf3
eukaryotic translation initiation factor 2, subunit 3, structural gene X-linked, Eif2s3x
amyotrophic lateral sclerosis 2 (juvenile) chromosome region, candidate 2 (human), Als2cr2
tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, epsilon polypeptide, Ywhae
neural precursor cell expressed, developmentally down-regulated gene 8, Nedd8
ubiquitin-conjugating enzyme E2E 3, UBC4/5 homolog (yeast), Ube2e3
ribosomal protein S4, X-linked, Rps4x
proliferating cell nuclear antigen, Pcna
cyclin-dependent kinase inhibitor 1C (P57), Cdkn1c
ZW10 interactor, Zwint
RIKEN cDNA 5031439G07 gene, 5031439G07Rik
interleukin enhancer binding factor 3, Ilf3
B-cell CLL/lymphoma 11A (zinc finger protein), Bcl11a
H2A histone family, member Y, H2afy
neurogenic differentiation 1, Neurod1
myristoylated alanine rich protein kinase C substrate, Marcks
RNA binding motif protein 9, Rbm9
WD repeat and SOCS box-containing 1, Wsb1
B-cell translocation gene 1, anti-proliferative, Btg1
roundabout homolog 1 (Drosophila), Robo1
EST from M2 cells of skin melanoma
cyclin-dependent kinase 4, Cdk4
secreted frizzled-related protein 1, Sfrp1
enhancer of zeste homolog 2 (Drosophila), Ezh2
zinc finger protein 57, Zfp57
SRY-box containing gene 11, Sox11
RT-qPCR validation of SAGE profile for P1.5 versus adult stages
proteolipid protein (myelin) 1, Plp1
hypoxanthine guanine phosphoribosyl transferase 1, Hprt1
calmodulin 1, Calm1
myelin basic protein, Mbp
down-regulator of transcription 1, Dr1
zinc finger, SWIM domain containing 4, Zswim4
actin-like 6B, Actl6b
ribosomal protein S4, X-linked, Rps4x
ATPase, Cu++ transporting, alpha polypeptide, Atp7a
high mobility group box 3, Hmgb3
ZW10 interactor, Zwint
C57BL/6 whole brain E15.5 (or known as embryonic lethal, abnormal vision, Drosophila-like 3 (Hu antigen C), Elavl3)
actin, beta, cytoplasmic, Actb
myristoylated alanine rich protein kinase C substrate, Marcks
cullin-associated and neddylation-dissociated 2 (putative), Cand2
cyclin-dependent kinase 4, Cdk4
EST from M2 cells of skin melanoma
SRY-box containing gene 4, Sox4
cysteine and glycine-rich protein 2, Csrp2
zinc finger protein 57, Zfp57
RT-qPCR validation of SAGE profile for E15.5 versus P1.5 stages
microtubule-associated protein tau, Mapt
SRY-box containing gene 11, Sox11
RT-qPCR validation of SAGE profile for rostral E15.5 and caudal E15.5 regions
bladder cancer associated protein homolog (human), Blcap
EST sequence BC025816
ankyrin repeat and zinc finger domain containing 1, Ankzf1
actin, beta, cytoplasmic, Actb
thymosin, beta 4, X chromosome, Tmsb4x
According to the comparison between two developmental stages, the most abundantly expressed and validated DETs in the E versus Ad analysis are Camk2a, Egr1 and Plp1 (Table 3). In the adult cerebral cortex, the expression of these DETs is more than 100-fold greater than that at E15.5. Other validated adult-specific DETs with expression levels of approximately tenfold or greater than those in the embryonic brain are Camk2n1, Cryab, Nrgn, BQ176089, CD802535, Sncb and Nptxr. Conversely, Sox11, an embryonic specific DET, is expressed in the E15.5 and E17.5 cerebral cortices with an expression level of at least 100-fold greater than that in the adult. Other validated embryonic specific DETs are Dcx, Zfp57, Ezh2, Sfrp1 and Cdk4, which have expression levels approximately tenfold or greater than those in the adult.
In the P1.5 versus Ad analysis, Dcx is expressed at a level approximately 80-fold greater in P1.5 compared to the adult cerebral cortex (Table 4). Other P1.5 enriched and validated DETs are (in descending order of enrichment) Zfp57, Csrp2, AA122503, Cdk4, Sox4, Marcks, Actb, BQ177889, Hmgb3, Rps4x, Actl6b, Zswim4 and Dr1, whose expression ranges from 33- to 1.4-fold greater than in the adult. On the other hand, the Plp1 is expressed at a level 40 times greater in the adult cerebral cortex compared to P1.5. Other validated genes that are enriched in the adult cerebral cortex include (in descending order) Hprt1, Calm1 and Mbp, with a 2.3- to 1.8-fold greater expression than the P1.5 cerebral cortex. Comparison between E15.5 and P1.5 shows that Mapt has a 1.5-fold greater expression level in the P1.5 cerebral cortex while Sox11 expression is 3.3-fold lower (Table 5).
We were unable to validate all 17 DETs from L versus Ri regions, suggesting the left and right hemispheres of the adult mouse cerebral cortex are highly similar and indistinguishable at the molecular level. SAGE and RT-qPCR analyses for R versus C regions of E15.5 are discussed in a separate section below.
Functional analysis of validated gene clusters using Ingenuity Pathway Analysis
The validated DETs of embryonic, adult and gene-switching clusters were functionally characterized using proprietary software, Ingenuity Pathway Analysis (IPA) from Ingenuity Systems®, to identify enriched molecular networks and canonical pathways. Given a list of input genes (known as focus genes), IPA mapped these genes to a global molecular network developed from information contained in the Ingenuity knowledge base (a manually curated database of experimentally proven molecular interactions from published literature). Networks of these focus genes were then algorithmically generated based on their connectivity. IPA determined the most significantly enriched biological functions and/or related diseases by calculating the P-value using Fisher's exact test. Using similar methods, significantly represented canonical pathways in a set of focus genes were also determined using IPA (Section C in Additional data file 1).
In adult-specific gene clusters, two molecular networks (18 focus genes and 50 associated nodes; networks 3 and 4 in Figure 3) were identified and interconnected via a single gene, Mbp (Figures S6 and S7 in Additional data file 1). These molecular networks enrich for nine statistically significant canonical pathways (P < 0.05) such as synaptic long-term potentiation and depression, calcium signaling, B cell receptor signaling, cAMP-mediated signaling, GM-CSF signaling, amyotrophic lateral sclerosis signaling, G-protein-coupled receptor signaling and xenobiotic metabolism signaling pathways. Validated DETs such as Camk2a, Gria3, Itpr1, Egr1, Rgs4, Gstm1 and Ppp3ca are associated with cell-to-cell communication in the adult cerebral cortex. Biological functions such as cell proliferation, movement and death processes are also linked to these networks through the genes Camk2a, Cryab, Egr1, Gstm1, Itpr1, Mbp, Nrgn, Ppp1r1b, Ppp3ca, Rgs4, Sept4, Sirpa and Sncb. Known human neurological disorders associated with Network 3 include Pelizaeus-Merzbacher disease ([OMIM:312080]; PLP1) and dementia with Lewy bodies ([OMIM:127750]; SNCB), whereas X-linked mental retardation ([OMIM:300210]; TSPAN7 or TM4SF2) is associated with network 4. MBP, found in both network 3 and 4, is associated with genetic susceptibility to multiple sclerosis [OMIM:126200]. All the DETs implicated in these networks are expressed in layers I to VI of the cerebral cortex, with the exception of Gstm1 (layers II to V), Sirpa and Ppp3ca (layers II to VI), and Nptxr (no expression data available) (Table S5 in Additional data file 1) [32–34, 43–52].
Both networks 5 and 6 (Figure 3) are linked to gene-switching clusters. These two molecular networks (14 focus genes and 52 associated nodes) are associated with cellular morphogenesis, amino acid metabolism, cell death processes and developmental disorders (Figures S8 and S9 in Additional data file 1). These networks are connected by Ctsd and Ywhae and are implicated in amyloid processing (App and Mapt), aryl hydrocarbon receptor signaling (Ctsd and Nedd8) and cell cycle G2/M DNA damage checkpoint regulation pathways (Ywhae). App, Atp7a, Hprt1, Mapt, Ctsd and Clcn2 are involved in cell morphogenesis, assembly and organization. Among all the networks in the study, network 5 has the greatest number of DETs associated with known human neurological disorders with six out of ten focus genes being associated with known neurological disorders such as neuronal ceroid lipofuscinosis ([OMIM:610127]; CTSD), Alzheimer's disease ([OMIM:104300]; APP), Menkes disease ([OMIM:309400]; ATP7A), X-linked cutis laxa ([OMIM:304150]; ATP7A), epilepsy ([OMIM:606904, 607628 and 607631]; CLCN2), Lesch-Nyhan syndrome ([OMIM:300322]; HPRT1), frontotemporal dementia ([OMIM:600274]; MAPT), progressive supranuclear palsy ([OMIM:601104]; MAPT) and Pick disease of the brain ([OMIM:172700]; MAPT). All the DETs implicated in these networks are expressed in the cortical plate during embryonic development with the exception of Ube2e3. In the adult, all DETs are expressed in layers I to VI of the cerebral cortex except for Bcl11a, Clcn2, Hprt1 (layers II to VI), Atp7a and Ube2e3 (no expression was detected in the adult cerebral cortex) (Table S5 in Additional data file 1) [32–35, 46, 53–61].
To refine the functional analysis to a cellular level, we grouped all 40 RT-qPCR validated DETs from networks 3, 4, 5 and 6 (adult specific and gene-switching gene clusters) into three groups according to where they are expressed: only in cortical neurons (N group); only in cortical glia (G group); and in both cortical neurons and glia (B group). None of the DETs from networks 1 and 2 (embryonic-specific gene clusters) were analyzed because most of the cells within the cerebral cortex are committed to the neuronal lineage at E15.5. The methods used to tabulate and group validated DETs are detailed in Section D in Additional data file 1. DETs classified as part of the N group are Bcl11a, Calm1, Camk2a, Camk2n1, Hprt1, Itpr1, Mapt, Nedd8, Nrgn, Ppp1r1b, Ppp3ca, Rbm9, Rgs4, Sncb and Ywhae, whereas those in the G group include Cryab, Mbp, Nptxr, Plp1, Ppap2b, Sept4 and Tspan7. The B group consists of App, Atp7a, Chgb, Ctsd, Egr1, Gria3, Gstm1 and Sirpa. DETs without any cellular expression information, such as Actl6b, Clcn2, Hmgb3, Ube2e3, AK138272, AK139402, AK140219, AK154943, AU258168 and BQ176089, were placed into group B to facilitate downstream analysis. Of the 40 DETs, only 35 were mapped to the IPA knowledgebase and subjected to further analysis. Functional analysis of these 35 DETs show that the neuron-based DETs (N and B groups) are enriched for various human neurological disorders, such as schizophrenia (9 DETs; P = 4.36E-8), Huntington's disease (8 DETs; P = 1.03E-5), X-linked mental retardation (4 DETs; P = 3.78E-7), Parkinson's disease (3 DETs; P = 6.99E-3) and Alzheimer's disease (3 DETs; P = 1.52E-2) (Figure S10 in Additional data file 1). Thirteen DETs associated with these neurological disorders are also implicated in processes related to cell death. Of these, eight DETs are expressed in neurons only, two in glia only and three in both neurons and glia.
From the above IPA analysis, Cand2 (embryonic gene clusters), Camk2n1 (adult gene clusters), Hmgb3 (gene-switching clusters) and 10 ESTs (BQ176089, CD802535, AK139402, AU258168, AK138272, AK140219, AK154943, AA122503, NM_001033273 and BQ177886) are not currently connected to any networks.
Regionalized expression of DETs in the E15.5 cerebral cortex
Early regionalized development is an important event that could lay the foundations for adult arealization of cerebral functions. To identify genes with regionalized expression profiles, we compared SAGE libraries generated from rostral and caudal regions (equivalent to anterior and posterior regions of the human cerebral cortex) of the E15.5 cerebral cortex. We identified 44 DETs and selected 25 DETs (22 known genes, 1 EST and 2 ambiguous genes; Additional data file 2) and 2 genes of interest (bladder cancer associated protein, Blcap, and ankyrin repeat and zinc finger domain containing 1, Ankzf1) for RT-qPCR validation and further detailed region-based analysis. This was done using independent biological triplicates of clearly defined regions/quadrants, such as rostro-lateral (RL), rostro-medial (RM), caudo-lateral (CL) and caudo-medial (CM) of the E15.5 cerebral cortex (see Materials and methods). Two positive controls with known regionalised expression were included in the RT-qPCR: RAR-related orphan receptor beta (Rorb) and nuclear receptor subfamily 2, group F, member 1 (Nr2f1). Rorb is highly expressed in the rostral region whereas Nr2f1 is highly expressed in the caudal region of the cerebral cortex [1, 62].
ISH analysis showed that both Blcap and Ankzf1 are caudal specific. Serial coronal sections from rostral to caudal regions of the brain (Figure 5g-i) show Blcap is weakly expressed in the rostral cerebral cortex, particularly at the intermediate zone, subplate and the cortical plate, but is highly expressed in the hippocampus and thalamus (Figure 5i). Ankzf1 is expressed specifically in the ventricular zone as well as the cortical plate towards the caudal region of the cerebral cortex (Figure 5j-m). The distinctive expression of Ankzf1 in both the ventricular zone and cerebral cortex prompted us to extend our ISH analysis to various developmental stages starting from E11.5 to adulthood (Figure 5n-t). Ankzf1 is expressed in the primordial plexiform and the ventricular zone layers of the telencephalon at E11.5 (Figure 5n). By E13.5, Ankzf1 is weakly expressed in the ventricular zone, but is highly expressed in the cortical plate (Figure 5o). From E17.5 to P1.5, the expression of Ankzf1 is maintained in the cerebral cortex (Figure 5p, q). In the adult, Ankzf1 expression is obvious in the piriform cortex, hippocampus and cerebellum (Figure 5r-t).
Genomic clustering of sense-antisense SAGE tags at the Sox4 and Sox11loci
We performed genomic clustering analysis of the SAGE tags to determine any actively transcribed chromosomal loci throughout cerebral corticogenesis. Probabilities for chance occurrences of two, three, four, and five DETs being clustered within a window of ten adjacent tags present within each chromosomal location, irrespective of genetic distance, were calculated. This analysis was based on the DET lists described above (Additional data file 2). The analysis showed two overrepresented chromosome loci at Sox4 and Sox11, which derive from embryonic-specific gene clusters.
Genomic cluster at the Sox4gene locus
To further validate the expression profiles of the multiple Sox4 DETs, we performed 3' rapid amplification of cDNA ends (RACE)-Southern analysis using pooled adaptor-oligo-d [T]15 synthesized cDNAs from three mice at each developmental stage. Based on this method, we were able to semi-quantitatively and accurately measure the expression levels of individual SAGE tags at the locus. To show that the amplification was cDNA specific, we performed PCR by using the same primer sets on mouse genomic DNA under the same conditions. In all cases, no amplification was observed (data not shown). This analysis confirmed the presence of four out of seven alternative transcripts for Sox4 (Figure 7d). Corresponding tags were determined by estimating the amplicon sizes between the strand-specific primers used, the next downstream AAUAAA/AUUAAA polyadenylation signal (if any) and succeeding CATG sequence or SAGE tags. Figure 7d(1)-d(3) confirms the existence of sox4_tag10, sox4_tag12 and sox4_tag15. Of these tags, SAGE expression profiles of sox4_tag10 and sox4_tag15 were validated (embryonic-specific and reduced expression after P1.5) but not sox4_tag12 (E15.5 caudal region-specific). 3' RACE-Southern analysis using a sense probe detected bands in the rostral and caudal regions of E15.5, E15.5 and E17.5 cerebral cortices and, therefore, confirmed the existence of the Sox4 antisense transcripts (Figure 7d(4)). Even though none of these tags were differentially expressed in between these regions based on the SAGE analysis, our findings show distinctive regionalization for sox4_tag14 expression at the E15.5 rostral cerebral cortex. Proteasome (prosome, macropain) subunit, beta type 2 gene (Psmb2) and Hmbs were used as controls and no antisense or alternative transcripts were identified at these gene loci (Figure 7d(5)-d(8)).
Since 3' RACE-Southern analysis was dependent on oligo- [dT]15 priming, we could not rule out the possibility of amplicons that were generated by false priming on homopolymer-A stretches. Therefore, Northern analyses were performed on equally pooled total RNA extracted from the cerebral cortices of seven mice at E15.5, E17.5 and P150 (negative control). By using a double-stranded DNA probe at the 3' untranslated region (UTR) of Sox4 (Additional data file 8), we identified six bands ranging from approximately 2 kb to approximately 4.7 kb (Figure 7e). Sox4 sense transcripts are weakly expressed in the adult, but highly expressed in the embryonic cerebral cortex. The number of bands observed is similar and corresponds to overlapping mouse mRNAs as well as publicly available paired-end diTag (PET) sequences downloaded from Ensembl (Figure 7a; Table S6 in Additional data file 1) . Taken together, the analysis confirmed the existence of multiple overlapping variants of Sox4 sense and antisense transcripts at this gene locus.
To confirm the rostro-caudal expression of Sox4 sense transcripts, we performed ISH on sagittal sections of mouse brains using a Sox4 antisense riboprobe that spanned across the sox4_tag10, sox4_tag12, sox4_tag15 and sox4_tag16 SAGE tags. Sox4 showed regionalized expression at E15.5 and E17.5 (Figure 7f). At E15.5, Sox4 sense transcripts are expressed more in both the rostral- and caudal-end regions of the cortical plate compared to the intermediate region between them (Figure 7f(1)). By E17.5, expression of Sox4 sense transcripts is obvious in the rostral cortical plate (Figure 7f(2)). At both stages of development, Sox4 sense transcripts are uniformly expressed in the intermediate zone of the cerebral cortex. These findings correspond to the SAGE tag counts for E15.5 rostro-caudal regions of the cerebral cortex (Figure 7f). These observations explain the averaged total tag count per 100,000 tags for different Sox4 sense transcripts, which are predominantly expressed in both rostral and caudal regions of the cerebral cortex. The regionalized expression of Sox4 in the cortical plate is obvious only at E15.5 and E17.5, but not at other stages of development (Figure S11 in Additional data file 1).
In situ RNA hybridization of Sox4sense and antisense transcripts
Analysis of the Sox11genomic cluster
SAGE tags, which represent multiple overlapping sense and antisense transcript variants at the Sox11 genomic cluster, were validated using 3' RACE-Southern analysis as described above. See Section F in Additional data file 1 for a full description of the Sox11 results. ISH analysis did not confirm the expression of antisense transcripts of Sox11, but the presence of PETs spanning three out of five antisense tags confirmed the existence of Sox11 antisense transcripts (Table S8 and Figures S16, S17 and S18 in Additional data file 1). The discrepancy between ISH and RT-qPCR or 3' RACE-Southern analysis suggests that Sox11 antisense transcripts might be expressed at low levels or at specific locations of the cerebral cortex, and hence can be detectable only by using serial sections or whole mount ISH.
Screening of Sox4 and Sox11antisense transcripts in the adult mouse brain, organs, P19 cell line and neurospheres
Sox11 sense transcripts are expressed at the highest level in the olfactory bulb, approximately two- to seven-fold greater than those in other brain regions (Figure 10a). Sox11 antisense transcripts, on the other hand, are expressed in all brain regions screened and at a comparable level in the olfactory bulb, hippocampus, thalamus and cerebral cortex. In comparison to other adult mouse organs, Sox11 sense and antisense transcripts are highly expressed in the E15.5 whole brain, with Sox11 sense transcript levels at least 100-fold greater than those in other mouse organs (Figure 10b). On the other hand, Sox11 antisense expression is observed only in the E15.5 whole brain, skin and stomach. Notably, Sox11 sense transcripts are expressed more highly than antisense transcripts in the E15.5 whole brain and skin (23- and 4-fold, respectively).
In this study, SAGE was used to analyze global gene expression in the normal mouse cerebral cortex at various developmental stages. We report validated spatio-temporal regulation of genes involved in mouse cerebral cortex development from embryo to adulthood. The study highlights four main findings: association of DETs from different gene clusters with known functional processes or signaling pathways and disease-causative genes that are involved in cerebral corticogenesis; Ankzf1 and Sox4 sense transcripts are regionally expressed in the E15.5 cerebral cortex; multiple overlapping Sox4 and Sox11 sense and antisense transcripts are spatio-temporally regulated during cerebral corticogenesis; and Sox4 and Sox11 antisense transcripts are differentially regulated in both proliferating and differentiating embryonic-derived neurospheres and P19 cells.
We have shown that most tags generated in all libraries were singletons. The number of singletons could be reduced by increasing the number of tags sequenced. In mammalian cells, the number of additional unique transcripts identified approached zero when the number of SAGE tags sequenced reached approximately 600,000 . Increasing the number of tags sequenced could improve the sensitivity of the technique to identify weakly expressed or novel transcripts, and the application of massively parallel signature sequencing  using a next-generation sequencer would be an ideal solution to accomplish this. In addition, one of the benefits of SAGE is that datasets generated from different groups or in public repositories such as SAGE Genie  and GEO  are readily comparable and, hence, can increase the tag count and sensitivity of the technique in discovering DETs between SAGE libraries. However, any meta-analyses involving various SAGE datasets are affected by experimental and biological variation; thus, a careful selection of matching libraries is crucial to limit systematic error or biases.
Our SAGE analysis robustly detected DETs with a low false positive rate (for example, <0.001% for comparison between left and right hemispheres of the adult cerebral cortex). Of all the identified DETs, approximately 8% were not mapped to either a single locus in the mouse genome or any unique annotation. This problem could be overcome by generating additional information from the 5' end of the transcript through alternative techniques such as PET sequencing , cap analysis gene expression (CAGE)  and 5' LongSAGE .
We have identified functional ontologies, molecular interactions and enriched canonical pathways that are distinct to the stage-specific gene clusters of validated DETs. The IPA network analysis generated connections between validated DETs across various developmental stages in relation to well-established proteins or molecules and neurological disorders. In the study, members of the Wnt/β-catenin signaling pathway were enriched in networks 1 and 2 (embryonic-specific gene clusters). In neural development, Wnt/β-catenin signaling plays an important role in regulating regional specification of the cortex along the rostro-caudal and dorso-ventral axes, and proliferation of cortical progenitors . IPA highlighted three genes (Sox4, Sox11 and Sfrp1) associated with this pathway. In humans, the SFRP1 protein (secreted frizzled-related protein 1) stabilizes β-catenin and increases transcription from β-catenin-responsive promoters . In β-catenin-deficient mouse mutants, expression of both Sox4 and Sox11 is downregulated [75, 76]. Sox4 and Sox11 proteins play an important role in establishing neuronal properties, pan-neuronal gene expression and proper myelination of the central nervous system [65, 66]. This suggests that the role of the Wnt/β-catenin signaling pathway in regulating neuronal development could be mediated, at least in part, by the Sox4 and Sox11 proteins.
The role of Wnt/β-catenin signaling in regulating DETs (Btg1, Cdk4, Cdkn1c, Csrp2, Ezh2, Neurod1, Pcna, and Rps4x) involved in cell cycle and proliferation remains unclear. Gain- and loss-of-function studies have established that Wnt/β-catenin signaling is essential to maintain the pool of precursors for proper development of the cerebral cortex [77, 78]. To date, there is no direct evidence to show that Ezh2, Pcna, Rps4x and Btg1 are involved in cell cycle regulation during early embryonic neurogenesis. But their expression in the ventricular/subventricular zone of the E15.5 developing mouse cerebral cortex [33, 35, 36] suggests that they may be involved in regulating cell proliferation during neurogenesis. The role of these DETs and their association with the Wnt/β-catenin signaling pathway remains unclear and requires further experimentation.
Networks 3 and 4 from the adult gene cluster were associated with various canonical pathways, in particular, synaptic long term potentiation (LTP) and calcium signaling. LTP is a process of synapse enhancement, which is thought to underlie some forms of learning and memory . This process depends on Ca2+ and calmodulin, which are major components of the calcium signaling pathway. We identified four DETs, Gria3, Itpr1, Ppp3ca and Camk2a, in both canonical pathways. In particular, the Camk2a enzyme is highly expressed in the brain and regulates mainly glutamatergic synapses during LTP . DET products such as Nrgn and Camk2n1 can directly or indirectly regulate Camk2a or Ca2+/calmodulin and subsequently alter the outcome of the LTP pathway [80, 81]. Therefore, they may serve as important candidate genes in the future analysis of the synaptic LTP pathway involving neurodegenerative diseases that cause the loss of cognitive function and memory.
Networks 5 and 6 are enriched for genes in the amyloid processing signaling pathway. App and Mapt are associated with this pathway. Under normal circumstances, the App protein is required for proper migration of neuronal precursors into the cortical plate in early embryonic corticogenesis . The Mapt protein, on the other hand, plays an important role in maintaining the architecture of the neuronal cytoskeleton and intracellular trafficking. Overexpression of App protein and hyperphosphorylation of the Mapt protein have been implicated in the pathologies of Alzheimer's disease [83, 84]. Interestingly, Ctsd, Atp7a, Clcn2 and Hprt1, the genes responsible for other human neurological disorders such as neuronal ceroid lipofuscinosis (Ctsd), Menkes disease (Atp7a), epilepsy (Clcn2), and Lesch-Nyhan syndrome (Hprt1), are associated with App and Mapt. These candidate genes are also involved in cell morphogenesis, assembly and organization and could be linked to deterioration of neurons during the pathologic progression of these disorders.
Pathway analysis of DETs classified into the N, G and B groups showed DETs in the neuron (N and B groups) are associated with Huntington's disease and schizophrenia, which were not previously identified in networks 3 to 6. Our analysis showed that both disorders share three common DETs, namely Rgs4, Ppp1r1b and Chgb, whose expression is downregulated in humans with Huntington's disease or schizophrenia [85–88]. A proportion of patients with Huntington's disease also develop schizophrenia [89, 90]. Taken together, downregulation of Rgs4, Ppp1r1b and Chgb expression in neurons may contribute to the common symptoms in these disorders. Our findings imply that many DETs (including App, Hprt1 and Sncb) associated with both Huntington's disease and schizophrenia are also involved in neuronal/cell death processes [91–93]. Other DETs in the N group, not previously implicated in neuronal cell death, may serve as novel potential candidates during pathologic development in these disorders.
Regionalized development of the cerebral cortex involves the differential regulation of cell cycle exit, early migration and attainment of positional identity in neuronal fated cells. To date, only few genes have been associated with regionalized development of the cerebral cortex [3, 94]. In the regionalization analysis, we identified the highest number of DETs in the comparison of the RL and CM libraries, which signifies that these two regions of the cerebral cortex are the most different. This finding supports the notion that the cerebral cortex is developed in a latero-medial axis followed by a rostro-caudal axis [7, 8]. At E15.5, both Actb and Tmsb4x were expressed greater in the rostral cerebral cortex than in the caudal region. Both Actb and Tmsb4x proteins are involved in the actin cytoskeleton-signaling pathway . In particular, the Tmsb4x protein has been shown to promote cardiomyocyte migration  and axonal tract growth in zebrafish . Therefore, co-expression of Actb and Tmsb4x in the E15.5 mouse cortical plate suggests that they may have a synergistic role in early cortical cell development. Conversely, Blcap and Ankzf1 were expressed more highly in the caudal than in the rostral region of the E15.5 cerebral cortex. To date, the function of both Blcap and Ankzf1 in the cerebral cortex remains uncharacterized. This study provides the first comprehensive expression profile of Ankzf1 and suggests it could be an important transcription factor in cerebral corticogenesis.
At E15.5, Sox4 sense transcripts were expressed in a high-rostral and high-caudal manner with lesser expression within the intermediate region. By E17.5, Sox4 expression becomes obvious at the rostral cortical plate, which is similar to Rorb. But, we did not find that the regionalized expression of Sox4 sense transcripts resembles that of the restricted Rorb expression at E15.5 or in the postnatal brain [1, 62]. This finding could be caused by the combined expression profiles of different Sox4 sense transcripts that are present across the rostro-caudal axis of the cortical plate. The regionalized expression of Sox4 sense transcripts occurs only between E15.5 and E17.5. Because the thalamic axon innervates the cortical plate after E17.5 , the regionalization of Sox4 sense transcripts in early cortical development could be an outcome of an intrinsic instead of an extrinsic mechanism that regulates early patterning of the cerebral cortex.
Genomic clustering of DETs identified the differentially regulated Sox4 and Sox11 gene loci. These genomic clusters imply that there are multiple overlapping sense and antisense transcripts surrounding the same gene locus that are co-transcribed simultaneously during cerebral cortex development. Both Sox4 and Sox11 are single exon genes and these transcript variants are therefore likely to be generated due to alternative polyadenylation. The 3' UTRs of both Sox4 and Sox11 have tandem terminal polyadenylation signals on both sense and antisense strands (data not shown), which supports the occurrence of multiple transcript forms or SAGE tags. Multiple mRNA forms with different 3' UTRs can lead to cell-specific regulation, different nuclear or cytoplasmic mRNA stability and translation rates [99, 100]. The 3' UTR of Sox4 and Sox11 may contain AU-rich elements that play an important role in determining mRNA stability through deadenylation, decapping or 3' → 5' decay . Besides, different 3' UTR lengths may be targeted by different miRNAs, thus interfering with the translation process. Both Sox4 and Sox11 transcripts may be targeted by various miRNAs at different predicted positions across the 3' UTR (Tables S7 and S9 and Figures S15 and S18 in Additional data file 1). Therefore, 3' UTR lengths of Sox4 and Sox11 may be an important feature in the regulation of their protein expression during cerebral corticogenesis.
In the study, NATs were found at both the Sox4 and Sox11 gene loci overlapping the sense transcripts. Overlapping NATs may function as templates for the generation of pre-miRNA and mature miRNA with exceptional high sequence conservation that complement the overlapping sense protein-coding transcripts . To date, no mature or pre-miRNAs have been predicted on the Sox4 and Sox11 sense and antisense strands (data not shown) or have been reported in miRBase . In addition, NATs can self-complement to form double stranded RNA or pair with sense transcripts and function as templates for the generation of endogenous small interfering RNAs, which could subsequently interfere with translation or transcription of multiple protein-coding transcripts . Because both Sox4 and Sox11 proteins are highly expressed in the cerebral cortex, the overlapping NATs do not seem to be involved in the regulation of Sox4 and Sox11 through miRNA- or small interfering RNA-mediated translation repression mechanisms, but rather through antisense-regulated sense transcription within the nucleus.
Our ISH analysis showed complementary cellular expression profiles of Sox4 sense and antisense transcripts at the piriform cortex, olfactory bulb and dentate gyrus. This finding implies that Sox4 antisense transcripts may be essential in intracellular and interlocus negative feedback loop regulation of the Sox4 sense transcripts. Similar expression profiles of Sox4 sense and antisense transcripts in multiple mouse organs and brain regions, however, suggest that these transcripts may be co-expressed. This observation is also supported by the temporal co-expression of Sox4 sense and antisense transcripts in the cortical plate or layers I to III of the cerebral cortex. Taken together, the sense and antisense transcripts of Sox4 are co-expressed in some cells and expressed complementarily in other cells, suggesting crucial cell-type-specific regulation.
Sox4 and Sox11 have been shown to have redundant roles during mouse development , and Sox11 may play a compensatory role in the absence of Sox4 during brain development . We demonstrated that Sox4 and Sox11 sense and antisense transcripts have a similar expression in the brain, but not in other organs, suggesting a compensatory role for Sox11 only in the brain. Sox11 antisense transcripts were expressed in the brain, skin and stomach only, suggesting organ-specific regulation.
Our data show upregulation of Sox4 and Sox11 sense transcripts in differentiating P19 cells, consistent with the findings of others [107, 108], and demonstrate upregulation of antisense transcripts as well. We also find both Sox4 and Sox11 sense transcripts expressed in the NSPCs cultured as neurospheres, which is in agreement with Dy et al. . Furthermore, we identify upregulation of both Sox4 sense and antisense transcripts but downregulation of Sox11 sense transcripts in differentiating neurospheres. Taken together, our findings show that there are potentially common and distinct roles for Sox4 and Sox11 sense and antisense transcripts during neuronal and non-neuronal cell proliferation and differentiation. The underlying regulatory mechanism of these transcripts, particularly the antisense ones, remains unknown and requires further investigation.
This study provides avenues for future research focus in understanding the fundamental processes and development of neurological disorders related to the cerebral cortex. We confirm the regionalized expression of new candidate genes in the E15.5 cerebral cortex as well as differential regulation of multiple overlapping sense and novel antisense transcripts within Sox4 and Sox11 gene loci during cerebral corticogenesis. We also report for the first time the spatio-temporal regulation of Sox4 antisense transcripts in the brain as well as differential regulation of novel Sox4 and Sox11 antisense transcripts in various mouse organs and in proliferating and differentiating NSPCs and P19 cells. The finding provides an insight for future investigations into the role of antisense transcripts during cerebral corticogenesis and neuronal differentiation.
Materials and methods
Serial analysis of gene expression
Handling of animals and dissection of the cerebral cortex
All experiments that involved animal breeding and handling were performed according to protocols approved by the Melbourne Health Animal Ethics Committee (Project numbers 2001.045 and 2004.041). All animals involved in the study were C57BL/6 mice unless specified otherwise. All mice were kept under conditions of a 12-h light/12-h dark cycle with unlimited access to food and water. All mice were culled by cervical dislocation prior to dissection. Cortical tissue was procured in the following fashion. For adult samples, after removal of the meninges, coronal cuts were used to excise the olfactory bulb from the rostral region, and the superior colliculus from the caudal region. A sagittal cut to separate the two cortical hemispheres was performed. The cortical pallium was dissected from the subpallial striatum and the septum. The neocortex was then dissected away from the cingulate cortex and the entorhinal cortex. For embryonic samples, the cortical tissue was dissected free from the underlying ganglionic eminences at the pallial-subpallial border. An orthogonal cut was made to remove the presumptive striatum and the overlying piriform cortex. On the medial aspect, the medial limbic cortex was included for analysis, but the adjacent hippocampal primordium, including the cortical hem, was excluded. For the E15.5 cerebral cortex, the resulting hemispheres containing cortical tissue only were placed on the bottom of the Petri dish and, using a fine scalpel, divided into four equal quadrants per hemisphere, namely RL, RM, CL and CM. Rostral and caudal quadrants from both hemispheres were pooled for SAGE library construction but separately tracked for RT-qPCR analysis. Procurement of other adult brain tissues and related mouse organs for Sox4 and Sox11 antisense transcript screening was carried out according to the standard mouse necropsy protocol accessible at the National Institute of Allergy and Infectious Diseases (NIAID) website .
SAGE libraries and analysis of tags
Ten SAGE libraries were constructed from the cerebral cortex of E15.5, E17.5 and 4- to 6-month-old (Ad) mice according to either one of the two methods described previously [25, 26], using I-SAGE™ or I-SAGE™ Long Kits (Invitrogen, Mulgrave, Victoria, Australia). Additional libraries from E15.5 and P1.5 of the cerebral cortex described previously  were also included in the analysis. These libraries contain a total of 26,436 traces. SAGE tags were preprocessed - that is, TAGs were extracted and corrected for sequencing errors, and artifacts like SAGE linkers, ribosomal RNA and duplicated ditags were removed using the 'sagenhaft' package, which is available from the Bioconductor website [111, 112]. To compare libraries that contain long tags with those that contain short tags, all short tags were mapped to the existing long tags from the other libraries. A table for all libraries containing the unique long or short tags was generated and redundant tags were removed. Only tags with a total count >2 (across all libraries) were considered for subsequent comparisons. Each unique tag was mapped to the mouse genome using ESTgraph, which employs ESTs and their genomic position information. ESTgraph was created by Tim Beissbarth (unpublished) . Identity was assigned to these tags and they were further grouped into the following categories: matching to a gene, a genomic sequence, or an EST, or ambiguous matches or no alignment at all. All annotations were based on the latest mouse assembly (mm9 released in July 2007) accessible from the UCSC Genome Bioinformatics website .
Identification of differentially expressed tags
Library comparisons were performed using two methods. Fisher's exact test was used to compare two individual SAGE libraries. In the analysis, multiple testing correction  was carried out to control for false-discovery rate and adjusted P-value cutoffs (Q-values) were used to select DETs. In cases where several libraries were combined to focus on a specific biological comparison (for example, different stages of development), a Bayesian model, as described previously , was used to integrate multiple libraries in pairwise comparisons involving biological replicates of libraries. The model accounts for within-class variability by means of mixture distributions. The resulting E-values were used to select DETs. A table of all relevant comparisons, the comparison method and Q- or E-value cutoffs is provided in Table S1 in Additional data file 1.
Hierarchical clustering of SAGE tags
To identify co-regulated genes, the clustering of DETs was performed based on the log2 of normalized counts. Each library was normalized to 100,000 tags per library to account for size differences. A pseudocount of 0.5 was added before taking the log2 of the normalized tag counts. The tag-wise mean was subtracted from the log2 tag intensities before computing the Euclidean distance of the individual tag profiles. Hierarchical clustering was performed on the tags using the 'hclust' function and complete linkage, which was implemented using the statistical computing environment of R .
Genomic clustering of SAGE tags
To assess whether there was any genomic clustering of tags, a method previously described  was adopted. In brief, first gene lists (based on all DETs in both pairwise and multiple library comparisons as well as gene lists from the hierarchical clustering analysis) were selected. The genomic clustering of either of these selections was compared to the total unique tag list (all 25,165 unique tags). The tags were mapped to the mouse genome. The number of selected tags in ten consecutive tag positions for each window of the chromosome was calculated. One thousand permutations were used to compute the null distribution of maximum tag counts per window. The method was implemented using the statistical computing environment of R .
Functional classification and characterization of DETs
Gene Ontology enrichment analysis
The DET lists generated from various comparisons were subjected to systematic functional annotation using the standardized Gene Ontology term analysis tools at the DAVID . Functional clustering was performed using high stringency with a kappa similarity threshold of 0.85 and a minimum term overlap of 3. Classification was carried out using a multiple linkage threshold of 0.5 with both numbers of initial and final group members set to 3. A term was considered statistically significant when the computed P-value was < 0.05. All queries were performed in September 2009.
Molecular interactions and pathway analysis
Identification of molecular network interactions and pathway analysis of validated DETs or co-regulated genes was completed using the IPA  tools from Ingenuity Systems® (Redwood City, California, USA). Accession numbers for all genes with their corresponding fold changes or normalized counts were imported into the IPA software. No focus genes were set at the beginning of the analysis. To start building networks, the application queries the list of input genes and all other gene objects stored in the Ingenuity knowledge base. Networks with a maximum of 30 genes or proteins were constructed, and scores were computed based on the likelihood of the genes being connected together due to random chance. A score of 2 indicates that there is a 1/100 chance that these genes are connected in a network due to random chance. Therefore, any networks with a score of 2 or above are considered statistically significant (with >99% confidence). The most significant novel networks and their interactions with existing canonical pathways were investigated further.
Total RNA isolation and first strand cDNA synthesis
RT-qPCR was carried out to validate all selected candidate or co-regulated DETs. Biological triplicates from E15.5, E17.5, P1.5 and adult (5 to 6 months old) cerebral cortices were used. Genomic-free total RNA from independent mice was extracted using the RNeasy® Lipid Tissue Midi Kit (Qiagen, Doncaster, Victoria, Australia) according to the manufacturer's protocol (Section G in Additional data file 1). First strand cDNA was synthesized from high quality total RNA  using random hexamers (unless specified otherwise) and the SuperScript™ III Reverse Transcriptase Kit (Invitrogen, Mulgrave, Victoria, Australia), according to the manufacturer's protocol, under an RNAse-free environment.
Primer design and RT-qPCR
All primers were designed using ProbeFinder 2.35 at the Assay Design Center for Universal ProbeLibrary Assay provided by Roche Applied Science . All RT-qPCR reactions were prepared in 10 μl volumes in a 384-well plate format consisting of LC480 Master Probe Mix (Roche Diagnostics, Castle Hill, New South Wales, Australia), Universal ProbeLibrary (Roche Diagnostics, Castle Hill, New South Wales, Australia) and forward and reverse primers (GeneWorks, Hindmarsh, South Australia, Australia or Bioneer Corporation, Daedeok-gu, Daejeon, Korea) according to the manufacturers' protocols. RT-qPCRs were performed using the LightCycler® 480 System (Roche Diagnostics, Castle Hill, New South Wales, Australia) according to the manufacturer's protocol (Section G in Additional data file 1). A full list of primers and probes used in this study is included in Additional data file 8.
Relative quantification using a standard curve method
The crossing point (Cp) from each signal was calculated based on the Second Derivative Maximum method . A set of serially diluted cDNAs was used to construct a four-data point standard curve for every PCR system in each run. A total of three reference genes (from Hprt1, Psmb2, phosphoglycerate kinase 1 gene (Pgk1) or Hmbs) were used as endogenous controls. An estimated starting amount of each target gene was calculated and intra-samples multiple reference genes normalization was performed (Section G in Additional data file 1). A linear model was fitted to the time course of expression values for each gene. Genes differentially expressed between the various stages of development or regions were selected using empirical Bayesian moderated t-statistics, which borrow information between genes . Standard errors for the mean expression at various developmental stages were obtained from the linear model. For each comparison, P-values were adjusted using the Benjamini and Hochberg  method to control the false discovery rate. See Section H in Additional data file 1 for the R code used to execute the analysis.
Validation of sense-antisense and multiple overlapping transcripts in genomic clusters
Strand specific RT-PCR
All RNA was prepared as described above. Total RNA from all developmental stages (N = 3 per developmental stage) was equally pooled prior to cDNA synthesis. Four first strand cDNA synthesis reactions were prepared for each cluster as follows: with a primer complementary to the sense strand only; with a primer complementary to the antisense strand only; with oligo-d [T]15 as a positive control; and without any primers as a negative control. In all four reactions, both primers were added in subsequent PCRs (Section G in Additional data file 1). PCR amplifications were carried out using FastStart PCR High Fidelity System (Roche Diagnostics, Castle Hill, New South Wales, Australia) according to the manufacturer's protocol. More than one primer set was used in the sense-antisense strand specific RT-PCR (Additional data file 8).
First strand cDNA synthesis was carried out using pooled total RNA extracted from three biological replicates of rostral and caudal E15.5 and whole E15.5, E17.5, P1.5 and adult (5 to 6 months old) cerebral cortices. Oligo-d [T]15 with an adaptor sequence (5'-TACGACGTCTGCTAGGACTG-3') was used to prime the first strand cDNA synthesis. Second strand synthesis or PCR was then carried out using a strand-specific primer and the adaptor primers (Additional data file 8). All specific primers were designed to be complementary to the SAGE tags or their upstream sequences. PCR amplifications (Section G in Additional data file 1) were carried out using FastStart PCR High Fidelity System (Roche Diagnostics, Castle Hill, New South Wales, Australia) according to the manufacturer's protocol.
Southern blotting analysis
Amplified 3' RACE products were transferred to Hybond N+TM (GE Healthcare, Rydalmere, New South Wales, Australia) nylon membrane using the neutral transfer method. Prehybridization and hybridization steps were performed in Rapid-Hyb buffer (GE Healthcare, Rydalmere, New South Wales, Australia) according to the manufacturer's protocol. All oligonucleotides were designed to be complementary to sequence between the specific primer-priming site and the tag of interest. Synthetic oligonucleotides were 5' end-labeled using T4 polynucleotide kinase (Promega, Alexandria, New South Wales, Australia) and [γ-32P]ATP (GE Healthcare, Rydalmere, New South Wales, Australia) with modifications to the manufacturer's protocol. After the hybridization step, the membrane was washed with 5× sodium chloride sodium citrate solution (with 0.1% v/v sodium dodecyl sulphate (SDS)) and 1× sodium chloride sodium citrate solution (with 0.1% w/v SDS) (Section G in Additional data file 1). See Additional data file 8 for primer sequences and oligonucleotides used for detection.
Northern blotting analysis
Independent preparations of total RNA from the cerebral cortex of seven mice at E15.5 and E17.5, and three adult mice were equally pooled to a final concentration of 20 μg per developmental stage. These pooled total RNAs were electrophoresed overnight and capillary transferred onto Hybond N+TM (GE Healthcare, Rydalmere, New South Wales, Australia) nylon membrane. Double-stranded DNA probes were radioactive-labeled using the Amersham Megaprime DNA Labeling System (GE Healthcare, Rydalmere, New South Wales, Australia) and [β-32P]CTP, according to the manufacturer's protocol. Hybridization was carried out overnight (approximately 18 h) at 65°C in aqueous buffer (7% w/v SDS with 0.5 M phosphate). After hybridization, blots were washed using 1% w/v SDS at 65°C for 5 to 6 times until the background signal was low.
In situRNA hybridization
ISH was carried out using paraffin sections (5 μm) of embryonic, postnatal and adult brains (E11.5, E13.5, E15.5, E17.5, P1.5 and P150) and a related [35S]UTP-labeled complementary RNA probe (Additional data file 8) as described previously (Section G in Additional data file 1) .
Screening of Sox4 and Sox11sense and antisense transcript expression in the adult mouse brain, organs and both the proliferating and differentiating P19 cells and neurospheres
Total RNA was extracted from harvested organs using the TRIzol®'s reagent (Invitrogen, Mulgrave, Victoria, Australia) according to the manufacturer's protocol. To avoid genomic DNA contamination, all isolated total RNA was treated with the recombinant DNAse I enzyme provided by the DNA-free™ kit (Applied Biosystems, Scoresby, Victoria, Australia) according to the manufacturer's protocol. First strand cDNA synthesis was carried out using strand-specific primers followed by qPCR analysis as described above.
Embryonic neural stem cells grown as neurospheres
Mouse used for generation of neurospheres had a mixed genetic background including MF1, 129SvEv, C57BL/6 and CBA. Cerebral cortices from E14 embryos were dissected out into calcium-magnesium-free phosphate-buffered saline. The tissue was mechanically dissociated and centrifuged. The cells were plated in complete neuroculture medium (Section G in Additional data file 1) for 4 days followed by induction of neuronal differentiation. These cells were then plated on poly-D-lysine (catalogue number P6407, Sigma Aldrich, Castle Hill, New South Wales, Australia) and laminin (catalogue number 23017-015, Invitrogen, Mulgrave, Victoria, Australia) coated culture dishes in neuroculture medium with the presence of 2% (v/v) fetal bovine serum but not epidermal growth factor and basic fibroblast growth factor. The differentiation was allowed to proceed for 5 days. Total RNA was extracted from both proliferating and differentiating cells using TRIzol® reagent as described above.
P19 embryonal carcinoma cells
P19 mouse embryonal carcinoma cells were cultured and differentiated into neurons as described previously . Briefly, P19 cell cultures were maintained in P19GM complete medium (Section G in Additional data file 1). For induction of neuronal differentiation, 1 × 106 P19 cells were cultured in suspension form using bacteriological Petri dishes. The P19GM medium with additional supplementation of 5 × 10-7M all-trans retinoic acid (catalogue number R-2625; Sigma Aldrich, Castle Hill, New South Wales, Australia) was used for the induction. After 4 days, P19 cells formed embryoid body stages. Embryoid bodies were collected from suspension cultures and re-plated in adherent culture flasks in the P19GM medium with only 5% (v/v) fetal bovine serum and without retinoic acid supplementation. The cells were allowed to differentiate for 5 days. Total RNA was extracted from both proliferating and differentiating cells using TRIzol® reagent as described in above.
Additional data files
The following additional data are available with the online version of this paper: analysis of SAGE, DETs, IPA, Sox4 and Sox11 genomic cluster analysis, and R script for implementing empirical Bayesian moderated t-test on multiple groups (Additional data file 1); SAGE tag information for 561 DETs (Additional data file 2); functional annotations clustering analysis using DAVID (Additional data file 3); RT-qPCR validation of DETs based on multiple comparisons between two developmental stages (E versus Ad, PN1.5 versus Ad and E15.5 versus PN1.5) (Additional data file 4); RT-qPCR validation of gene clusters based on hierarchical clustering analysis (Additional data file 5); RT-qPCR validation of DETs based on the rostral versus caudal E15.5 cerebral cortex comparison (Additional data file 6); statistically significant over-represented genomic loci based on genomic clustering of tags (Additional data file 7); list of primers, probes, clones and assays designed for RT-qPCR, RACE, Southern, Northern and ISH analysis (Additional data file 8).
Database for Annotation, Visualization and Integrated Discovery
differentially expressed transcript/tag
expressed sequence tag
Gene Expression Omnibus
Ingenuity Pathway Analysis
in situ hybridization
long term potentiation
natural antisense transcript
neural stem/progenitor cell
Online Mendelian Inheritance in Man
rapid amplification of cDNA ends
serial analysis of gene expression
sodium dodecyl sulphate
University of California Santa Cruz
We thank Felix Ng, Sarah King and Caitlin Collins for technical assistance; Frank Weissenborn for assistance with tissue dissections; Ken Simpson and Keith Satterley for computational assistance. The P19 cell line was kindly provided by Dr K Jensen, School of Molecular and Biomedical Science, University of Adelaide, South Australia, Australia. This work was supported by a National Health and Medical Research Council fellowships (171601 and 461204 to HSS); National Health and Medical Research Council Grants 219176, 257501, 215201, and 257529 (to HSS, GKS, and S-ST), a fellowship from the Deutsche Forschungsgemeinschaft (to TB). PQT is a Pfizer Australia Research Fellow. K-HL was a recipient of the Melbourne International Fee Remission Scholarship (MIFRS) and Universiti Putra Malaysia Staff Training Scholarship (UPMSTS), and Adelaide International Fee Scholarship (AIFS) equivalent. The majority of DNA sequencing in production of SAGE libraries was performed by the Australian Genome Research Facility, which was established through the Commonwealth-funded Major National Research Facilities program.
- Miyashita-Lin EM, Hevner R, Wassarman KM, Martinez S, Rubenstein JL: Early neocortical regionalization in the absence of thalamic innervation. Science. 1999, 285: 906-909. 10.1126/science.285.5429.906.PubMedGoogle Scholar
- Molyneaux BJ, Arlotta P, Menezes JR, Macklis JD: Neuronal subtype specification in the cerebral cortex. Nat Rev Neurosci. 2007, 8: 427-437. 10.1038/nrn2151.PubMedGoogle Scholar
- O'Leary DD, Chou SJ, Sahara S: Area patterning of the mammalian cortex. Neuron. 2007, 56: 252-269. 10.1016/j.neuron.2007.10.010.PubMedGoogle Scholar
- Gupta A, Tsai LH, Wynshaw-Boris A: Life is a journey: a genetic look at neocortical development. Nat Rev Genet. 2002, 3: 342-355. 10.1038/nrg799.PubMedGoogle Scholar
- Marin O, Rubenstein JL: A long, remarkable journey: tangential migration in the telencephalon. Nat Rev Neurosci. 2001, 2: 780-790. 10.1038/35097509.PubMedGoogle Scholar
- Nadarajah B, Parnavelas JG: Modes of neuronal migration in the developing cerebral cortex. Nat Rev Neurosci. 2002, 3: 423-432.PubMedGoogle Scholar
- Smart IH: Histogenesis of the mesocortical area of the mouse telencephalon. J Anat. 1984, 138: 537-552.PubMedPubMed CentralGoogle Scholar
- Uylings HBM, Van Eden CG, Parnavelas JG, Kalsbeek A: The prenatal and postnatal development of rat cerebral cortex. The Cerebral Cortex of the Rat. Edited by: Kolb B, Tees RC. 1990, Cambridge, MA: MIT Press, 36-76.Google Scholar
- Bishop KM, Goudreau G, O'Leary DD: Regulation of area identity in the mammalian neocortex by Emx2 and Pax6. Science. 2000, 288: 344-349. 10.1126/science.288.5464.344.PubMedGoogle Scholar
- Mallamaci A, Muzio L, Chan CH, Parnavelas J, Boncinelli E: Area identity shifts in the early cerebral cortex of Emx2-/- mutant mice. Nat Neurosci. 2000, 3: 679-686. 10.1038/76630.PubMedGoogle Scholar
- Armentano M, Chou SJ, Tomassy GS, Leingartner A, O'Leary DD, Studer M: COUP-TFI regulates the balance of cortical patterning between frontal/motor and sensory areas. Nat Neurosci. 2007, 10: 1277-1286. 10.1038/nn1958.PubMedGoogle Scholar
- Sahara S, Kawakami Y, Izpisua Belmonte JC, O'Leary DD: Sp8 exhibits reciprocal induction with Fgf8 but has an opposing effect on anterior-posterior cortical area patterning. Neural Dev. 2007, 2: 10-10.1186/1749-8104-2-10.PubMedPubMed CentralGoogle Scholar
- Zembrzycki A, Griesel G, Stoykova A, Mansouri A: Genetic interplay between the transcription factors Sp8 and Emx2 in the patterning of the forebrain. Neural Dev. 2007, 2: 8-10.1186/1749-8104-2-8.PubMedPubMed CentralGoogle Scholar
- Bulfone A, Puelles L, Porteus MH, Frohman MA, Martin GR, Rubenstein JL: Spatially restricted expression of Dlx-1, Dlx-2 (Tes-1), Gbx-2, and Wnt-3 in the embryonic day 12.5 mouse forebrain defines potential transverse and longitudinal segmental boundaries. J Neurosci. 1993, 13: 3155-3172.PubMedGoogle Scholar
- Gaszner M, Felsenfeld G: Insulators: exploiting transcriptional and epigenetic mechanisms. Nat Rev Genet. 2006, 7: 703-713. 10.1038/nrg1925.PubMedGoogle Scholar
- Kadonaga JT: Eukaryotic transcription: an interlaced network of transcription factors and chromatin-modifying machines. Cell. 1998, 92: 307-313. 10.1016/S0092-8674(00)80924-1.PubMedGoogle Scholar
- Ladd-Acosta C, Pevsner J, Sabunciyan S, Yolken RH, Webster MJ, Dinkins T, Callinan PA, Fan JB, Potash JB, Feinberg AP: DNA methylation signatures within the human brain. Am J Hum Genet. 2007, 81: 1304-1315. 10.1086/524110.PubMedPubMed CentralGoogle Scholar
- Fox AH, Lam YW, Leung AK, Lyon CE, Andersen J, Mann M, Lamond AI: Paraspeckles: a novel nuclear domain. Curr Biol. 2002, 12: 13-25. 10.1016/S0960-9822(01)00632-7.PubMedGoogle Scholar
- Prasanth KV, Prasanth SG, Xuan Z, Hearn S, Freier SM, Bennett CF, Zhang MQ, Spector DL: Regulating gene expression through RNA nuclear retention. Cell. 2005, 123: 249-263. 10.1016/j.cell.2005.08.033.PubMedGoogle Scholar
- Schratt GM, Tuebing F, Nigh EA, Kane CG, Sabatini ME, Kiebler M, Greenberg ME: A brain-specific microRNA regulates dendritic spine development. Nature. 2006, 439: 283-289. 10.1038/nature04367.PubMedGoogle Scholar
- Stark KL, Xu B, Bagchi A, Lai WS, Liu H, Hsu R, Wan X, Pavlidis P, Mills AA, Karayiorgou M, Gogos JA: Altered brain microRNA biogenesis contributes to phenotypic deficits in a 22q11-deletion mouse model. Nat Genet. 2008, 40: 751-760. 10.1038/ng.138.PubMedGoogle Scholar
- Bak M, Silahtaroglu A, Moller M, Christensen M, Rath MF, Skryabin B, Tommerup N, Kauppinen S: MicroRNA expression in the adult mouse central nervous system. Rna. 2008, 14: 432-444. 10.1261/rna.783108.PubMedPubMed CentralGoogle Scholar
- Tochitani S, Hayashizaki Y: Nkx2.2 antisense RNA overexpression enhanced oligodendrocytic differentiation. Biochem Biophys Res Commun. 2008, 372: 691-696. 10.1016/j.bbrc.2008.05.127.PubMedGoogle Scholar
- Smalheiser NR, Lugli G, Torvik VI, Mise N, Ikeda R, Abe K: Natural antisense transcripts are co-expressed with sense mRNAs in synaptoneurosomes of adult mouse forebrain. Neurosci Res. 2008, 62: 236-239. 10.1016/j.neures.2008.08.010.PubMedPubMed CentralGoogle Scholar
- Saha S, Sparks AB, Rago C, Akmaev V, Wang CJ, Vogelstein B, Kinzler KW, Velculescu VE: Using the transcriptome to annotate the genome. Nat Biotechnol. 2002, 20: 508-512. 10.1038/nbt0502-508.PubMedGoogle Scholar
- Velculescu VE, Zhang L, Vogelstein B, Kinzler KW: Serial analysis of gene expression. Science. 1995, 270: 484-487. 10.1126/science.270.5235.484.PubMedGoogle Scholar
- Gunnersen JM, Augustine C, Spirkoska V, Kim M, Brown M, Tan SS: Global analysis of gene expression patterns in developing mouse neocortex using serial analysis of gene expression. Mol Cell Neurosci. 2002, 19: 560-573. 10.1006/mcne.2001.1098.PubMedGoogle Scholar
- NCBI Gene Expression Omnibus. [http://www.ncbi.nlm.nih.gov/geo/]
- Boon WM, Beissbarth T, Hyde L, Smyth G, Gunnersen J, Denton DA, Scott H, Tan SS: A comparative analysis of transcribed genes in the mouse hypothalamus and neocortex reveals chromosomal clustering. Proc Natl Acad Sci USA. 2004, 101: 14972-14977. 10.1073/pnas.0406296101.PubMedPubMed CentralGoogle Scholar
- Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4: P3-10.1186/gb-2003-4-5-p3.PubMedGoogle Scholar
- Yuan X, Chittajallu R, Belachew S, Anderson S, McBain CJ, Gallo V: Expression of the green fluorescent protein in the oligodendrocyte lineage: a transgenic mouse for developmental and physiological studies. J Neurosci Res. 2002, 70: 529-545. 10.1002/jnr.10368.PubMedGoogle Scholar
- Allen Institute for Brain Science. [http://mouse.brain-map.org/]
- Brain Gene Expression Map. [http://www.stjudebgem.org/]
- Gene Expression Nervous System Atlas. [http://www.gensat.org/]
- GenePaint.org. [http://www.genepaint.org/]
- Ino H, Chiba T: Expression of proliferating cell nuclear antigen (PCNA) in the adult and developing mouse nervous system. Brain Res Mol Brain Res. 2000, 78: 163-174. 10.1016/S0169-328X(00)00092-9.PubMedGoogle Scholar
- Lopez-Coviella I, Follettie MT, Mellott TJ, Kovacheva VP, Slack BE, Diesl V, Berse B, Thies RS, Blusztajn JK: Bone morphogenetic protein 9 induces the transcriptome of basal forebrain cholinergic neurons. Proc Natl Acad Sci USA. 2005, 102: 6984-6989. 10.1073/pnas.0502097102.PubMedPubMed CentralGoogle Scholar
- Bhat NR, Zhang P, Bhat AN: The expression of myristoylated alanine-rich C-kinase substrate in oligodendrocytes is developmentally regulated. Dev Neurosci. 1995, 17: 256-263. 10.1159/000111294.PubMedGoogle Scholar
- Calabrese B, Halpain S: Essential role for the PKC target MARCKS in maintaining dendritic spine morphology. Neuron. 2005, 48: 77-90. 10.1016/j.neuron.2005.08.027.PubMedGoogle Scholar
- 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: 3395-3407. 10.1523/JNEUROSCI.4605-06.2007.PubMedGoogle Scholar
- Sher F, Rossler R, Brouwer N, Balasubramaniyan V, Boddeke E, Copray S: Differentiation of neural stem cells into oligodendrocytes: involvement of the polycomb group protein Ezh2. Stem Cells. 2008, 26: 2875-2883. 10.1634/stemcells.2008-0121.PubMedGoogle Scholar
- Mermelstein F, Yeung K, Cao J, Inostroza JA, Erdjument-Bromage H, Eagelson K, Landsman D, Levitt P, Tempst P, Reinberg D: Requirement of a corepressor for Dr1-mediated repression of transcription. Genes Dev. 1996, 10: 1033-1048. 10.1101/gad.10.8.1033.PubMedGoogle Scholar
- Beiswanger CM, Diegmann MH, Novak RF, Philbert MA, Graessle TL, Reuhl KR, Lowndes HE: Developmental changes in the cellular distribution of glutathione and glutathione S-transferases in the murine nervous system. Neurotoxicology. 1995, 16: 425-440.PubMedGoogle Scholar
- Saha S, Datta K, Rangarajan P: Characterization of mouse neuronal Ca2+/calmodulin kinase II inhibitor alpha. Brain Res. 2007, 1148: 38-42. 10.1016/j.brainres.2007.02.018.PubMedGoogle Scholar
- Wilson CA, Murphy DD, Giasson BI, Zhang B, Trojanowski JQ, Lee VM: Degradative organelles containing mislocalized alpha-and beta-synuclein proliferate in presenilin-1 null neurons. J Cell Biol. 2004, 165: 335-346. 10.1083/jcb.200403061.PubMedPubMed CentralGoogle Scholar
- Bachoo RM, Kim RS, Ligon KL, Maher EA, Brennan C, Billings N, Chan S, Li C, Rowitch DH, Wong WH, DePinho RA: Molecular diversity of astrocytes with implications for neurological disorders. Proc Natl Acad Sci USA. 2004, 101: 8384-8389. 10.1073/pnas.0402140101.PubMedPubMed CentralGoogle Scholar
- Lencesova L, O'Neill A, Resneck WG, Bloch RJ, Blaustein MP: Plasma membrane-cytoskeleton-endoplasmic reticulum complexes in neurons and astrocytes. J Biol Chem. 2004, 279: 2885-2893. 10.1074/jbc.M310365200.PubMedGoogle Scholar
- Zeng H, Chattarji S, Barbarosie M, Rondi-Reig L, Philpot BD, Miyakawa T, Bear MF, Tonegawa S: Forebrain-specific calcineurin knockout selectively impairs bidirectional synaptic plasticity and working/episodic-like memory. Cell. 2001, 107: 617-629. 10.1016/S0092-8674(01)00585-2.PubMedGoogle Scholar
- Kim PM, Aizawa H, Kim PS, Huang AS, Wickramasinghe SR, Kashani AH, Barrow RK, Huganir RL, Ghosh A, Snyder SH: Serine racemase: activation by glutamate neurotransmission via glutamate receptor interacting protein and mediation of neuronal migration. Proc Natl Acad Sci USA. 2005, 102: 2105-2110. 10.1073/pnas.0409723102.PubMedPubMed CentralGoogle Scholar
- Mandir AS, Poitras MF, Berliner AR, Herring WJ, Guastella DB, Feldman A, Poirier GG, Wang ZQ, Dawson TM, Dawson VL: NMDA but not non-NMDA excitotoxicity is mediated by Poly(ADP-ribose) polymerase. J Neurosci. 2000, 20: 8005-8011.PubMedGoogle Scholar
- Yang X, Skoff RP: Proteolipid protein regulates the survival and differentiation of oligodendrocytes. J Neurosci. 1997, 17: 2056-2070.PubMedGoogle Scholar
- Southwood CM, Peppi M, Dryden S, Tainsky MA, Gow A: Microtubule deacetylases, SirT2 and HDAC6, in the nervous system. Neurochem Res. 2007, 32: 187-195. 10.1007/s11064-006-9127-6.PubMedGoogle Scholar
- Niciu MJ, Ma XM, El Meskini R, Ronnett GV, Mains RE, Eipper BA: Developmental changes in the expression of ATP7A during a critical period in postnatal neurodevelopment. Neuroscience. 2006, 139: 947-964. 10.1016/j.neuroscience.2006.01.044.PubMedGoogle Scholar
- Sanosaka T, Namihira M, Asano H, Kohyama J, Aisaki K, Igarashi K, Kanno J, Nakashima K: Identification of genes that restrict astrocyte differentiation of midgestational neural precursor cells. Neuroscience. 2008, 155: 780-788. 10.1016/j.neuroscience.2008.06.039.PubMedGoogle Scholar
- Mbebi C, See V, Mercken L, Pradier L, Muller U, Loeffler JP: Amyloid precursor protein family-induced neuronal death is mediated by impairment of the neuroprotective calcium/calmodulin protein kinase IV-dependent signaling pathway. J Biol Chem. 2002, 277: 20979-20990. 10.1074/jbc.M107948200.PubMedGoogle Scholar
- Jinnah HA, Hess EJ, Wilson MC, Gage FH, Friedmann T: Localization of hypoxanthine-guanine phosphoribosyltransferase messenger RNA in the mouse brain by in situ hybridization. Mol Cell Neurosci. 1992, 3: 64-78. 10.1016/1044-7431(92)90010-Y.PubMedGoogle Scholar
- Poet M, Kornak U, Schweizer M, Zdebik AA, Scheel O, Hoelter S, Wurst W, Schmitt A, Fuhrmann JC, Planells-Cases R, Mole SE, Hubner CA, Jentsch TJ: Lysosomal storage disease upon disruption of the neuronal chloride transport protein ClC-6. Proc Natl Acad Sci USA. 2006, 103: 13854-13859. 10.1073/pnas.0606137103.PubMedPubMed CentralGoogle Scholar
- Palfi A, Kortvely E, Fekete E, Kovacs B, Varszegi S, Gulya K: Differential calmodulin gene expression in the rodent brain. Life Sci. 2002, 70: 2829-2855. 10.1016/S0024-3205(02)01544-8.PubMedGoogle Scholar
- Toyo-oka K, Shionoya A, Gambello MJ, Cardoso C, Leventer R, Ward HL, Ayala R, Tsai LH, Dobyns W, Ledbetter D, Hirotsune S, Wynshaw-Boris A: 14-3-3epsilon is important for neuronal migration by binding to NUDEL: a molecular explanation for Miller-Dieker syndrome. Nat Genet. 2003, 34: 274-285. 10.1038/ng1169.PubMedGoogle Scholar
- Kuo TY, Hsueh YP: Expression of zinc finger transcription factor Bcl11A/Evi9/CTIP1 in rat brain. J Neurosci Res. 2007, 85: 1628-1636. 10.1002/jnr.21300.PubMedGoogle Scholar
- Li Q, Lee JA, Black DL: Neuronal regulation of alternative pre-mRNA splicing. Nat Rev Neurosci. 2007, 8: 819-831. 10.1038/nrn2237.PubMedGoogle Scholar
- Zhou C, Tsai SY, Tsai MJ: COUP-TFI: an intrinsic factor for early regionalization of the neocortex. Genes Dev. 2001, 15: 2054-2059. 10.1101/gad.913601.PubMedPubMed CentralGoogle Scholar
- UCSC Genome Bioinformatics. [http://genome.ucsc.edu/]
- Ensembl Genome Browser. [http://www.ensembl.org/]
- Bergsland M, Werme M, Malewicz M, Perlmann T, Muhr J: The establishment of neuronal properties is controlled by Sox4 and Sox11. Genes Dev. 2006, 20: 3475-3486. 10.1101/gad.403406.PubMedPubMed CentralGoogle Scholar
- Potzner MR, Griffel C, Lutjen-Drecoll E, Bosl MR, Wegner M, Sock E: Prolonged Sox4 expression in oligodendrocytes interferes with normal myelination in the central nervous system. Mol Cell Biol. 2007, 27: 5316-5326. 10.1128/MCB.00339-07.PubMedPubMed CentralGoogle Scholar
- Velculescu VE, Madden SL, Zhang L, Lash AE, Yu J, Rago C, Lal A, Wang CJ, Beaudry GA, Ciriello KM, Cook BP, Dufault MR, Ferguson AT, Gao Y, He TC, Hermeking H, Hiraldo SK, Hwang PM, Lopez MA, Luderer HF, Mathews B, Petroziello JM, Polyak K, Zawel L, Kinzler KW, et al: Analysis of human transcriptomes. Nat Genet. 1999, 23: 387-388. 10.1038/70487.PubMedGoogle Scholar
- Brenner S, Johnson M, Bridgham J, Golda G, Lloyd DH, Johnson D, Luo S, McCurdy S, Foy M, Ewan M, Roth R, George D, Eletr S, Albrecht G, Vermaas E, Williams SR, Moon K, Burcham T, Pallas M, DuBridge RB, Kirchner J, Fearon K, Mao J, Corcoran K: Gene expression analysis by massively parallel signature sequencing (MPSS) on microbead arrays. Nat Biotechnol. 2000, 18: 630-634. 10.1038/76469.PubMedGoogle Scholar
- The Cancer Genome Anatomy Project: SAGE Genie. [http://cgap.nci.nih.gov/SAGE]
- Ng P, Tan JJ, Ooi HS, Lee YL, Chiu KP, Fullwood MJ, Srinivasan KG, Perbost C, Du L, Sung WK, Wei CL, Ruan Y: Multiplex sequencing of paired-end ditags (MS-PET): a strategy for the ultra-high-throughput analysis of transcriptomes and genomes. Nucleic Acids Res. 2006, 34: e84-10.1093/nar/gkl444.PubMedPubMed CentralGoogle Scholar
- Shiraki T, Kondo S, Katayama S, Waki K, Kasukawa T, Kawaji H, Kodzius R, Watahiki A, Nakamura M, Arakawa T, Fukuda S, Sasaki D, Podhajska A, Harbers M, Kawai J, Carninci P, Hayashizaki Y: Cap analysis gene expression for high-throughput analysis of transcriptional starting point and identification of promoter usage. Proc Natl Acad Sci USA. 2003, 100: 15776-15781. 10.1073/pnas.2136655100.PubMedPubMed CentralGoogle Scholar
- Wei CL, Ng P, Chiu KP, Wong CH, Ang CC, Lipovich L, Liu ET, Ruan Y: 5' Long serial analysis of gene expression (LongSAGE) and 3' LongSAGE for transcriptome characterization and genome annotation. Proc Natl Acad Sci USA. 2004, 101: 11701-11706. 10.1073/pnas.0403514101.PubMedPubMed CentralGoogle Scholar
- Chenn A: Wnt/beta-catenin signaling in cerebral cortical development. Organogenesis. 2008, 4: 76-80.PubMedPubMed CentralGoogle Scholar
- Bafico A, Gazit A, Pramila T, Finch PW, Yaniv A, Aaronson SA: Interaction of frizzled related protein (FRP) with Wnt ligands and the frizzled receptor suggests alternative mechanisms for FRP inhibition of Wnt signaling. J Biol Chem. 1999, 274: 16180-16187. 10.1074/jbc.274.23.16180.PubMedGoogle Scholar
- Morkel M, Huelsken J, Wakamiya M, Ding J, Wetering van de M, Clevers H, Taketo MM, Behringer RR, Shen MM, Birchmeier W: Beta-catenin regulates Cripto- and Wnt3-dependent gene expression programs in mouse axis and mesoderm formation. Development. 2003, 130: 6283-6294. 10.1242/dev.00859.PubMedGoogle Scholar
- Fevr T, Robine S, Louvard D, Huelsken J: Wnt/beta-catenin is essential for intestinal homeostasis and maintenance of intestinal stem cells. Mol Cell Biol. 2007, 27: 7551-7559. 10.1128/MCB.01034-07.PubMedPubMed CentralGoogle Scholar
- Machon O, Backman M, Machonova O, Kozmik Z, Vacik T, Andersen L, Krauss S: A dynamic gradient of Wnt signaling controls initiation of neurogenesis in the mammalian cortex and cellular specification in the hippocampus. Dev Biol. 2007, 311: 223-237. 10.1016/j.ydbio.2007.08.038.PubMedGoogle Scholar
- Woodhead GJ, Mutch CA, Olson EC, Chenn A: Cell-autonomous beta-catenin signaling regulates cortical precursor proliferation. J Neurosci. 2006, 26: 12620-12630. 10.1523/JNEUROSCI.3180-06.2006.PubMedPubMed CentralGoogle Scholar
- Lisman J, Schulman H, Cline H: The molecular basis of CaMKII function in synaptic and behavioural memory. Nat Rev Neurosci. 2002, 3: 175-190. 10.1038/nrn753.PubMedGoogle Scholar
- Pak JH, Huang FL, Li J, Balschun D, Reymann KG, Chiang C, Westphal H, Huang KP: Involvement of neurogranin in the modulation of calcium/calmodulin-dependent protein kinase II, synaptic plasticity, and spatial learning: a study with knockout mice. Proc Natl Acad Sci USA. 2000, 97: 11232-11237. 10.1073/pnas.210184697.PubMedPubMed CentralGoogle Scholar
- Lepicard EM, Mizuno K, Antunes-Martins A, von Hertzen LS, Giese KP: An endogenous inhibitor of calcium/calmodulin-dependent kinase II is up-regulated during consolidation of fear memory. Eur J Neurosci. 2006, 23: 3063-3070. 10.1111/j.1460-9568.2006.04830.x.PubMedGoogle Scholar
- Young-Pearse TL, Bai J, Chang R, Zheng JB, LoTurco JJ, Selkoe DJ: A critical function for beta-amyloid precursor protein in neuronal migration revealed by in utero RNA interference. J Neurosci. 2007, 27: 14459-14469. 10.1523/JNEUROSCI.4701-07.2007.PubMedGoogle Scholar
- Alonso A, Zaidi T, Novak M, Grundke-Iqbal I, Iqbal K: Hyperphosphorylation induces self-assembly of tau into tangles of paired helical filaments/straight filaments. Proc Natl Acad Sci USA. 2001, 98: 6923-6928. 10.1073/pnas.121119298.PubMedPubMed CentralGoogle Scholar
- Rapoport M, Dawson HN, Binder LI, Vitek MP, Ferreira A: Tau is essential to beta-amyloid-induced neurotoxicity. Proc Natl Acad Sci USA. 2002, 99: 6364-6369. 10.1073/pnas.092136199.PubMedPubMed CentralGoogle Scholar
- Carter CJ: eIF2B and oligodendrocyte survival: where nature and nurture meet in bipolar disorder and schizophrenia?. Schizophr Bull. 2007, 33: 1343-1353. 10.1093/schbul/sbm007.PubMedPubMed CentralGoogle Scholar
- Craddock N, O'Donovan MC, Owen MJ: The genetics of schizophrenia and bipolar disorder: dissecting psychosis. J Med Genet. 2005, 42: 193-204. 10.1136/jmg.2005.030718.PubMedPubMed CentralGoogle Scholar
- Hodges A, Strand AD, Aragaki AK, Kuhn A, Sengstag T, Hughes G, Elliston LA, Hartog C, Goldstein DR, Thu D, Hollingsworth ZR, Collin F, Synek B, Holmans PA, Young AB, Wexler NS, Delorenzi M, Kooperberg C, Augood SJ, Faull RL, Olson JM, Jones L, Luthi-Carter R: Regional and cellular gene expression changes in human Huntington's disease brain. Hum Mol Genet. 2006, 15: 965-977. 10.1093/hmg/ddl013.PubMedGoogle Scholar
- Luthi-Carter R, Strand A, Peters NL, Solano SM, Hollingsworth ZR, Menon AS, Frey AS, Spektor BS, Penney EB, Schilling G, Ross CA, Borchelt DR, Tapscott SJ, Young AB, Cha JH, Olson JM: Decreased expression of striatal signaling genes in a mouse model of Huntington's disease. Hum Mol Genet. 2000, 9: 1259-1271. 10.1093/hmg/9.9.1259.PubMedGoogle Scholar
- Lovestone S, Hodgson S, Sham P, Differ AM, Levy R: Familial psychiatric presentation of Huntington's disease. J Med Genet. 1996, 33: 128-131. 10.1136/jmg.33.2.128.PubMedPubMed CentralGoogle Scholar
- Shiwach R: Psychopathology in Huntington's disease patients. Acta Psychiatr Scand. 1994, 90: 241-246. 10.1111/j.1600-0447.1994.tb01587.x.PubMedGoogle Scholar
- Chandra S, Gallardo G, Fernandez-Chacon R, Schluter OM, Sudhof TC: Alpha-synuclein cooperates with CSPalpha in preventing neurodegeneration. Cell. 2005, 123: 383-396. 10.1016/j.cell.2005.09.028.PubMedGoogle Scholar
- Evert BO, Wullner U, Klockgether T: Cell death in polyglutamine diseases. Cell Tissue Res. 2000, 301: 189-204. 10.1007/s004410000228.PubMedGoogle Scholar
- Morishima Y, Gotoh Y, Zieg J, Barrett T, Takano H, Flavell R, Davis RJ, Shirasaki Y, Greenberg ME: Beta-amyloid induces neuronal apoptosis via a mechanism that involves the c-Jun N-terminal kinase pathway and the induction of Fas ligand. J Neurosci. 2001, 21: 7551-7560.PubMedGoogle Scholar
- Mallamaci A, Stoykova A: Gene networks controlling early cerebral cortex arealization. Eur J Neurosci. 2006, 23: 847-856. 10.1111/j.1460-9568.2006.04634.x.PubMedGoogle Scholar
- Philp D, Huff T, Gho YS, Hannappel E, Kleinman HK: The actin binding site on thymosin beta4 promotes angiogenesis. Faseb J. 2003, 17: 2103-2105.PubMedGoogle Scholar
- Bock-Marquette I, Saxena A, White MD, Dimaio JM, Srivastava D: Thymosin beta4 activates integrin-linked kinase and promotes cardiac cell migration, survival and cardiac repair. Nature. 2004, 432: 466-472. 10.1038/nature03000.PubMedGoogle Scholar
- Roth LW, Bormann P, Bonnet A, Reinhard E: beta-thymosin is required for axonal tract formation in developing zebrafish brain. Development. 1999, 126: 1365-1374.PubMedGoogle Scholar
- Lopez-Bendito G, Molnar Z: Thalamocortical development: how are we going to get there?. Nat Rev Neurosci. 2003, 4: 276-289. 10.1038/nrn1075.PubMedGoogle Scholar
- Winter J, Kunath M, Roepcke S, Krause S, Schneider R, Schweiger S: Alternative polyadenylation signals and promoters act in concert to control tissue-specific expression of the Opitz Syndrome gene MID1. BMC Mol Biol. 2007, 8: 105-10.1186/1471-2199-8-105.PubMedPubMed CentralGoogle Scholar
- Abdel Wahab N, Gibbs J, Mason RM: Regulation of gene expression by alternative polyadenylation and mRNA instability in hyperglycaemic mesangial cells. Biochem J. 1998, 336: 405-411.PubMedPubMed CentralGoogle Scholar
- Wilusz CJ, Wilusz J: Bringing the role of mRNA decay in the control of gene expression into focus. Trends Genet. 2004, 20: 491-497. 10.1016/j.tig.2004.07.011.PubMedGoogle Scholar
- Lu C, Jeong D-H, Kulkarni K, Pillay M, Nobuta K, German R, Thatcher SR, Maher C, Zhang L, Ware D, Liu B, Cao X, Meyers BC, Green PJ: Genome-wide analysis for discovery of rice microRNAs reveals natural antisense microRNAs (nat-miRNAs). Proc Natl Acad Sci USA. 2008, 105: 4951-4956. 10.1073/pnas.0708743105.PubMedPubMed CentralGoogle Scholar
- miRBase. [http://microrna.sanger.ac.uk/]
- Watanabe T, Totoki Y, Toyoda A, Kaneda M, Kuramochi-Miyagawa S, Obata Y, Chiba H, Kohara Y, Kono T, Nakano T, Surani MA, Sakaki Y, Sasaki H: Endogenous siRNAs from naturally formed dsRNAs regulate transcripts in mouse oocytes. Nature. 2008, 453: 539-543. 10.1038/nature06908.PubMedGoogle Scholar
- Hoser M, Potzner MR, Koch JM, Bosl MR, Wegner M, Sock E: Sox12 deletion in the mouse reveals nonreciprocal redundancy with the related Sox4 and Sox11 transcription factors. Mol Cell Biol. 2008, 28: 4675-4687. 10.1128/MCB.00338-08.PubMedPubMed CentralGoogle Scholar
- Cheung M, Abu-Elmagd M, Clevers H, Scotting PJ: Roles of Sox4 in central nervous system development. Brain Res Mol Brain Res. 2000, 79: 180-191. 10.1016/S0169-328X(00)00109-1.PubMedGoogle Scholar
- Jin Z, Liu L, Bian W, Chen Y, Xu G, Cheng L, Jing N: Different transcription factors regulate nestin gene expression during P19 cell neural differentiation and central nervous system development. J Biol Chem. 2009, 284: 8160-8173. 10.1074/jbc.M805632200.PubMedPubMed CentralGoogle Scholar
- Lee MS, Jun DH, Hwang CI, Park SS, Kang JJ, Park HS, Kim J, Kim JH, Seo JS, Park WY: Selection of neural differentiation-specific genes by comparing profiles of random differentiation. Stem Cells. 2006, 24: 1946-1955. 10.1634/stemcells.2005-0325.PubMedGoogle Scholar
- Dy P, Penzo-Mendez A, Wang H, Pedraza CE, Macklin WB, Lefebvre V: The three SoxC proteins - Sox4, Sox11 and Sox12 - exhibit overlapping expression patterns and molecular properties. Nucleic Acids Res. 2008, 36: 3101-3117. 10.1093/nar/gkn162.PubMedPubMed CentralGoogle Scholar
- Virtual Mouse Necropsy. [http://www3.niaid.nih.gov/labs/aboutlabs/cmb/InfectiousDiseasePathogenesisSection/mouseNecropsy/]
- Beissbarth T, Hyde L, Smyth GK, Job C, Boon WM, Tan SS, Scott HS, Speed TP: Statistical modeling of sequencing errors in SAGE libraries. Bioinformatics. 2004, 20 (Suppl 1): I31-I39. 10.1093/bioinformatics/bth924.PubMedGoogle Scholar
- Bioconductor. [http://www.bioconductor.org/]
- Mapping SAGE Tags to Genes Using EST and Genomic Information. [http://www.mbgproject.org/tagmapping/]
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc. 1995, 57: 289-300.Google Scholar
- Vencio RZ, Brentani H, Patrao DF, Pereira CA: Bayesian model accounting for within-class biological variability in serial analysis of gene expression (SAGE). BMC Bioinformatics. 2004, 5: 119-10.1186/1471-2105-5-119.PubMedPubMed CentralGoogle Scholar
- The R Project for Statistical Computing. [http://www.r-project.org/]
- Gotter J, Brors B, Hergenhahn M, Kyewski B: Medullary epithelial cells of the human thymus express a highly diverse selection of tissue-specific genes colocalized in chromosomal clusters. J Exp Med. 2004, 199: 155-166. 10.1084/jem.20031677.PubMedPubMed CentralGoogle Scholar
- Ingenuity Pathway Analysis. [http://www.ingenuity.com]
- Schroeder A, Mueller O, Stocker S, Salowsky R, Leiber M, Gassmann M, Lightfoot S, Menzel W, Granzow M, Ragg T: The RIN: an RNA integrity number for assigning integrity values to RNA measurements. BMC Mol Biol. 2006, 7: 3-10.1186/1471-2199-7-3.PubMedPubMed CentralGoogle Scholar
- Assay Design Center for Universal ProbeLibrary Assay. [http://www.roche-applied-science.com/sis/rtpcr/upl/ezhome.html]
- Luu-The V, Paquet N, Calvo E, Cumps J: Improved real-time RT-PCR method for high-throughput measurements using second derivative calculation and double correction. Biotechniques. 2005, 38: 287-293. 10.2144/05382RR05.PubMedGoogle Scholar
- Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article3Google Scholar
- Thomas T, Voss AK, Chowdhury K, Gruss P: Querkopf, a MYST family histone acetyltransferase, is required for normal cerebral cortex development. Development. 2000, 127: 2537-2548.PubMedGoogle Scholar
- Rudnicki MA, McBurney MW: Cell culture methods and induction of differentiation of embryonal carcinoma cell lines. Teratocarcinomas and Embryonic Stem Cells: A Practical Approach. Edited by: Robertson EJ. 1987, Oxford: IRL Press, 19-49.Google Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.