Skip to main content

1Identification of genes differentially expressed in the embryonic pig cerebral cortex before and after appearance of gyration

Abstract

Background

Mammalian evolution is characterized by a progressive expansion of the surface area of the cerebral cortex, an increase that is accompanied by gyration of the cortical surface. The mechanisms controlling this gyration process are not well characterized but mutational analyses indicate that genes involved in neuronal migration play an important function. Due to the lack of gyration of the rodent brain it is important to establish alternative models to examine brain development during the gyration process. The pig brain is gyrated and accordingly is a candidate alternative model.

Findings

In this study we have identified genes differentially expressed in the pig cerebral cortex before and after appearance of gyration. Pig cortical tissue from two time points in development representing a non-folded, lissencephalic, brain (embryonic day 60) and primary-folded, gyrencephalic, brain (embryonic day 80) were examined by whole genome expression microarray studies. 91 differentially expressed transcripts (fold change >3) were identified. 84 transcripts were annotated and encoding proteins involved in for example neuronal migration, calcium binding, and cytoskeletal structuring. Quantitative real-time PCR was used to confirm the regulation of a subset of the identified genes.

Conclusion

This study provides identification of genes which are differentially expressed in the pig cerebral cortex before and after appearance of brain gyration. The identified genes include novel candidate genes which could have functional importance for brain development.

Findings

The complex architecture of the mammalian cerebral cortex is a consequence of the highly organized movement of neuronal cells. All neurons populating the six layered cerebral cortex undergo mitosis in distant compartments and achieve their final position following migration [1]. This migration is coordinated to obtain specific laminar position, orientation, and connections with other neurons. The number of neurons is dramatically increased throughout mammalian evolution but the expansion in the radial direction is relative limited, and the thickness of the cortex is relative constant [1, 2]. The main expansion is observed in the tangential direction which is the most variable and distinctive part of the central nervous system [3]. The proliferative zones lining the ventricular surface generate the neurons that migrate to their proper position guided by scaffolds formed by the transient extended fibers of the radial glial cells [4]. The radial glia cells define a radial unit consisting of a relative constant number of neurons in all mammalian species investigated [1, 57]. However, the number of radial units is increased throughout mammalian evolution. This increase in neuronal number is accomplished through extension of the time of neurogensis and decreasing the rate of neuronal death [3, 8]. Furthermore, in higher mammals the proliferative zones are expanded into structured subventricular zones (SVZs) enabling them to further amplify the neuronal number [6, 9, 10]. The neuronal output is increased through generation of intermediate progenitor cells (IPC) that divide symmetrically in the SVZ and generate either two neurons or two IPC and thereby amplify the number of neurons [11, 12].

The expansion of the mammalian brain during evolution is orchestrated with development of sulci and gyri, the convoluted folds of the cerebral cortex [13]. The process of cortical convolution is still poorly understood. The remarkable similarity of the gyration patterns among members within a species but with different patterns between species indicates that the convolution is a highly programmed process. Recent work suggested that the sites of gyral and sulcal formation can be predicted from the size of the SVZ [14]. The increase in primate SVZ complexity can, however, not solely explain the evolution of a gyrencephalic cortical surface since both lissencephalic and gyrencephalic brain structures are represented in diverse mammalian groups, including primates and rodents [6, 15]. Another hypothesis suggests that mechanical forces exerted by axonal fibers are generating the gyrencephalic cortical structure [16, 17]. The axonal and radial fibers attaching the growing cerebral cortex to the centre of the developing brain represent elastic elements in which the plasticity changes during elongation. Hence, the strongly interconnected cortical regions are pulled together resulting in gyri formation whereas weakly connected regions drift apart generating sulci [17, 18].

Mutations causing lissencephaly, smooth brain, are identified in genes involved in neuronal migration suggesting that this step in brain development indeed constitutes a major determinant in establishment of the convoluted structure of the cerebral cortex. The first identified gene involved in human lissencephaly was encoding for the β-subunit of platelet-activating factor acetylhydrolase 1b, also known as lissencephaly 1 (LIS1) [19]. Doublecortin (DCX) was identified as the gene responsible for X-linked lissencephaly in males and subcortical heterotopia in females [20]. Mutations in the human Reelin gene are associated with recessive lissencephaly with cerebellar hypoplasia [21]. Altogether a large set of genetic data imply that cytoskeletal rearrangements and neuronal migration are crucial to development of gyrencephalic brains.

The pig brain as a model for mammalian brain gyration

To study the process of brain gyration the pig brain constitutes an attractive alternative model to more classical laboratory animals. The convolution of the pig brain is occurring in the period between embryonic day E60 and E80 of the 117 day gestation period and the relative developmental timing and the anatomical structure is comparable to the primate brain (Nielsen et al. unpublished results and Additional file 1). At the microscopic level the process of neuronal migration shows a pronounced similarity between primates and the pig and the neuronal proliferating layers and number of migrating neurons are more prominent at E60 compared to E80 (Nielsen et al. unpublished results).

We questioned which genes are differentially expressed between E60 and E80 and accordingly are candidate genes to be involved in brain development during this time frame. First we examined if a group of genes described to be involved in human brain convolution were differently expressed between E60 and E80 in pig cortical tissue. mRNA was isolated and by quantitative real-time PCR (qRT-PCR) the expression levels of such genes were measured. In the analysis we examined the expression of DCX, aristaless-related homeobox (ARX), G protein coupled receptor 56 isoform A (GPR56), filamin A gene (FLNA), Reelin, VLDLR, ApoER2, Dab1, FYN, LIS1, nuclear distribution element-like (NDEL1), and cyclin dependent kinase 5 (CDK5). For DCX, ARX, GPR56, FLNA, VLDLR, Dab1, and NDEL1 we observed a down-regulation of expression (Figure 1). LIS1 was the only of the genes up-regulated (Figure 1). Reelin, ApoER2, FYN, and CDK5 were not significantly up- or down-regulated (Figure 1). Examining also later cortical developmental time points supported the overall tendency in regulation Additional file 2. From the qRT-PCR expression analysis we conclude that some of the genes involved in the brain convolution process have an altered transcriptional activity during the timeframe in where convolution appears.

Figure 1
figure 1

Expression analysis of genes involved in mammalian brain convolution. mRNA was extracted from pig cortical tissue from E60 and E80. qRT-PCR analysis were performed on cDNA for the genes DCX, ARX, GPR56, FLNA, Reelin, VLDLR, ApoER2, Dab1, FYN, LIS1, NDEL1, and CDK5. The expression levels were normalized to GAPDH, Beta-actin and 18S rRNA expression using the geNorm program [30].

Expression microarray analysis

For a comprehensive analysis of the differences in gene expression profiles in the pig cerebral cortex between E60 and E80 we utilized the Affymetrix GeneChip® Pig Genome Array. This array contains 24,123 probes including 23,256 pig transcripts, which represents 20,201 pig genes [22]. Cortical RNA isolated from three pig E60 embryos and three E80 embryos were used for the microarray screening. Array data are available at Gene Expression Omnibus [GEO:GSE18467]. The inter-chip variance was small with a scale factor (SF) between 0.403 and 0.558 Additional file 3. The number of expressed genes present on each chip was very similar, 66%, in accordance with the number of genes hypothesised to be active during brain development Additional file 3[23]. Genes consistently up-regulated or down-regulated on the microarray triplets were further analysed. An arbitrary threshold was set at 3-fold before we classified a gene for differently regulated. This high threshold was selected to assure the minimization of background noise. 98 transcripts were found to be differential expressed by this stringency (Table 1 and Table 2). Six transcripts were represented more than once. Five transcripts were present two times; Glial fibrillary acidic protein (GFAP), the homeobox gene Meis1, megalencephalic leukoencephalopathy with subcortical cysts gene 1 (MLC1), Non-SMC element 1 homolog (NSE1), and proteolipid protein 1 (PLP) and one transcript was present three times; DEAD-box protein 17 (DBX17). Hence the total number of identified differential expressed genes was 91. Of the six genes represented multiple times on the microarray consistency in up and down regulation was observed.

Table 1 Whole pig genome microarray results.
Table 2 Unannotated genes.

10 of the genes differentially expressed between the two gestational time points were selected for qRT-PCR verification. These genes were GFAP, apolipoprotein E (ApoE), Calbindin 2, Neurofilament heavy chain (200 kDa), S100 calcium binding protein A1 (S100A1), Tubulin-alpha1 (TUBA1), Neurogranin, Actinin-alpha2 (ACTN2), N-chimaerin (CHN1), and DBX17 (Figure 2). For all the examined genes we observed the same tendency in regulation by the qRT-PCR analysis and the microarray analysis but the differences in fold regulation were not equal illustrating the use of two fundamental different detection approaches (see also Additional file 4). However, we conclude that an overall consistency exists between the microarray and qRT-PCR data.

Figure 2
figure 2

Verification of microarray data with qRT-PCR. Nine genes differently expressed from the microarray analysis were examined for the expression level using mRNA extracted from pig cortical tissue from E60 and E80. qRT-PCR analysis were performed on cDNA for the genes GFAP, ApoE, calbindin-2, Neurofilament (200 kDa), S100A1, TUBA1, Neurogranin, ACTN2, CHN1, and Dbx17. The expression levels were normalized to GAPDH, Beta-actin and 18S rRNA expression using the geNorm program [30].

The degree of annotation made available by Affymetrix covers only about 10% of the genes represented on the chip. Further annotations were accomplished through comparing the Affymetrix pig target sequences by BLAST against the Ensembl human cDNA sequence database or retrieved from the annotation list by Tsai et al., 2006 [24]. Of the identified 91 differently expressed transcripts 84 could be annotated by this method (Table 1). We note that among the annotated genes were for example TUBA1 and MLC1, both directly correlated with brain abnormalities in humans suffering of lissencephaly and macrocephaly, respectively [2529]. Seven transcripts could not be annotated [GenBank: CO986932; BI402064; BF712758; CN155998; AU060035; CF367810; BF712467] (Table 2). Gene ontology (GO) http://www.geneontology.org annotations were determined for each individual transcript. GO terms covers the consistent descriptions of gene represented in different databases including biological processes, cellular components in which they exist, and the molecular functions they perform. The majority of the identified differentially regulated genes are highly expressed in the brain. Using biological process annotations, the genes differential expressed between E60 and E80 could be classified into groups according to calcium binding, cytoskeleton organization and biogenesis, transcription activation, signal transduction, and CNS development. A subset of the genes could not be placed in any of these functional groups.

Conclusions

The aim of this study was to identify genes which are differentially expressed during the time of gyration of the pig cerebral cortex. It is important to notice that by screening these two time points we are not screening for genes specifically involved in brain gyration but for genes those expression is changed during the developmental and differentiation processes occurring in the time frame before and after appearance of brain gyration. We have identified for us a surprisingly low number of genes to be differently expressed between the examined embryonic time points supporting that the majority of the cortical cells have not undergone specific differentiation processes at E80 compared to E60. In this line it is important to state that neuronal migration is evident both at E60 and E80, but decreased at the latter time point. We have identified several differentially expressed genes that are described to be functional involved in neuronal migration, apoptosis, angiogenesis, myelination, and brain gyration but also a number of genes not characterised for such functions and which accordingly could be interesting new candidate genes. Further analysis will be required to determine the function of these genes during brain development.

Abbreviations

ACTN2:

Actinin-alpha2

ApoE:

apolipoprotein E

ARX:

aristaless-related homeobox

CDK5:

cyclin dependent kinase 5

CHN1:

N-chimaerin

DBX17:

DEAD-box protein 17

DCX:

Doublecortin

E60:

embryonic day 60

E80:

embryonic day 80

FLNA:

filamin A gene

GFAP:

Glial fibrillary acidic protein

GO:

Gene ontology

GPR56:

G protein coupled receptor 56 isoform A

IPCs:

intermediate progenitor cells

LIS1:

lissencephaly 1

MLC1:

megalencephalic leukoencephalopathy with subcortical cysts gene 1

NDEL1:

nuclear distribution element-like

NSE1:

Non-SMC element 1 homolog

PLP:

proteolipid protein 1

qRT-PCR:

quantitative real-time PCR

S100A1:

S100 calcium binding protein A1

SF:

scale factor

SVZs:

subventricular zones

TUBA1:

Tubulin-alpha1

References

  1. Rakic P: A small step for the cell, a giant leap for mankind: a hypothesis of neocortical expansion during evolution. Trends Neurosci. 1995, 18 (9): 383-388. 10.1016/0166-2236(95)93934-P.

    Article  CAS  PubMed  Google Scholar 

  2. Fish JL, Dehay C, Kennedy H, Huttner WB: Making bigger brains-the evolution of neural-progenitor-cell division. J Cell Sci. 2008, 121 (Pt 17): 2783-2793. 10.1242/jcs.023465.

    Article  CAS  PubMed  Google Scholar 

  3. Caviness VS, Takahashi T, Nowakowski RS: Numbers, time and neocortical neuronogenesis: a general developmental and evolutionary model. Trends in Neurosciences. 1995, 18 (9): 379-383. 10.1016/0166-2236(95)93933-O.

    Article  CAS  PubMed  Google Scholar 

  4. Rakic P: Mode of cell migration to the superficial layers of fetal monkey neocortex. J Comp Neurol. 1972, 145 (1): 61-83. 10.1002/cne.901450105.

    Article  CAS  PubMed  Google Scholar 

  5. Rakic P: Specification of cerebral cortical areas. Science. 1988, 241 (4862): 170-176. 10.1126/science.3291116.

    Article  CAS  PubMed  Google Scholar 

  6. Cheung AFP, Pollen AA, Tavare A, DeProto J, Molnar Z: Comparative aspects of cortical neurogenesis in vertebrates. Journal of Anatomy. 2007, 211 (2): 164-176. 10.1111/j.1469-7580.2007.00769.x.

    Article  PubMed Central  PubMed  Google Scholar 

  7. Herculano-Houzel S, Collins CE, Wong P, Kaas JH, Lent R: The basic nonuniformity of the cerebral cortex. Proc Natl Acad Sci USA. 2008, 105 (34): 12593-12598. 10.1073/pnas.0805417105.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  8. Kornack DR, Rakic P: Changes in cell-cycle kinetics during the development and evolution of primate neocortex. Proc Natl Acad Sci USA. 1998, 95 (3): 1242-1246. 10.1073/pnas.95.3.1242.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  9. Lukaszewicz A, Savatier P, Cortay V, Giroud P, Huissoud C, Berland M, Kennedy H, Dehay C: G1 phase regulation, area-specific cell cycle control, and cytoarchitectonics in the primate cortex. Neuron. 2005, 47 (3): 353-364. 10.1016/j.neuron.2005.06.032.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  10. Smart IH, Dehay C, Giroud P, Berland M, Kennedy H: Unique morphological features of the proliferative zones and postmitotic compartments of the neural epithelium giving rise to striate and extrastriate cortex in the monkey. Cereb Cortex. 2002, 12 (1): 37-53. 10.1093/cercor/12.1.37.

    Article  PubMed Central  PubMed  Google Scholar 

  11. Haubensak W, Attardo A, Denk W, Huttner WB: Neurons arise in the basal neuroepithelium of the early mammalian telencephalon: a major site of neurogenesis. Proc Natl Acad Sci USA. 2004, 101 (9): 3196-3201. 10.1073/pnas.0308600100.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  12. Kriegstein AR, Noctor SC: Patterns of neuronal migration in the embryonic cortex. Trends Neurosci. 2004, 27 (7): 392-399. 10.1016/j.tins.2004.05.001.

    Article  CAS  PubMed  Google Scholar 

  13. Martinez-Cerdeno V, Noctor SC, Kriegstein AR: The Role of Intermediate Progenitor Cells in the Evolutionary Expansion of the Cerebral Cortex. Cereb Cortex. 2006, 16 (suppl_1): i152-161. 10.1093/cercor/bhk017.

    Article  PubMed  Google Scholar 

  14. Kriegstein A, Noctor S, Martinez-Cerdeno V: Patterns of neural stem and progenitor cell division may underlie evolutionary cortical expansion. Nat Rev Neurosci. 2006, 7 (11): 883-890. 10.1038/nrn2008.

    Article  CAS  PubMed  Google Scholar 

  15. Pontious A, Kowalczyk T, Englund C, Hevner RF: Role of intermediate progenitor cells in cerebral cortex development. Dev Neurosci. 2008, 30 (1-3): 24-32. 10.1159/000109848.

    Article  CAS  PubMed  Google Scholar 

  16. Hilgetag CC, Barbas H: Role of mechanical factors in the morphology of the primate cerebral cortex. PLoS Comput Biol. 2006, 2 (3): e22-10.1371/journal.pcbi.0020022.

    Article  PubMed Central  PubMed  Google Scholar 

  17. Van Essen DC: A tension-based theory of morphogenesis and compact wiring in the central nervous system. Nature. 1997, 385 (6614): 313-318. 10.1038/385313a0.

    Article  CAS  PubMed  Google Scholar 

  18. Toro R, Burnod Y: A Morphogenetic Model for the Development of Cortical Convolutions. Cereb Cortex. 2005, 15 (12): 1900-1913. 10.1093/cercor/bhi068.

    Article  PubMed  Google Scholar 

  19. Reiner O, Carrozzo R, Shen Y, Wehnert M, Faustinella F, Dobyns WB, Caskey CT, Ledbetter DH: Isolation of a Miller-Dicker lissencephaly gene containing G protein [beta]-subunit-like repeats. Nature. 1993, 364 (6439): 717-721. 10.1038/364717a0.

    Article  CAS  PubMed  Google Scholar 

  20. Ayala R, Shu T, Tsai LH: Trekking across the brain: the journey of neuronal migration. Cell. 2007, 128 (1): 29-43. 10.1016/j.cell.2006.12.021.

    Article  CAS  PubMed  Google Scholar 

  21. Hong SE, Shugart YY, Huang DT, Shahwan SA, Grant PE, Hourihane JOB, Martin NDT, Walsh CA: Autosomal recessive lissencephaly with cerebellar hypoplasia is associated with human RELN mutations. Nat Genet. 2000, 26 (1): 93-96. 10.1038/79246.

    Article  CAS  PubMed  Google Scholar 

  22. GeneChip® Porcine Genome Array. [http://www.affymetrix.com/browse/products.jsp?productId=131488&navMode=34000&navAction=jump&aId=productsNav#1_3]

  23. O'Connor A: Public Atlas Of the Brain In the Works. In The New York Times Volume September 16, 2003. 2003, New York The New York Times Company

    Google Scholar 

  24. Tsai S, Cassady JP, Freking BA, Nonneman DJ, Rohrer GA, Piedrahita JA: Annotation of the Affymetrix1 porcine genome microarray. Animal Genetics. 2006, 37 (4): 423-424. 10.1111/j.1365-2052.2006.01460.x.

    Article  CAS  PubMed  Google Scholar 

  25. Keays DA, Tian G, Poirier K, Huang GJ, Siebold C, Cleak J, Oliver PL, Fray M, Harvey RJ, Molnar Z, et al: Mutations in alpha-tubulin cause abnormal neuronal migration in mice and lissencephaly in humans. Cell. 2007, 128 (1): 45-57. 10.1016/j.cell.2006.12.017.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  26. Poirier K, Keays DA, Francis F, Saillour Y, Bahi N, Manouvrier S, Fallet-Bianco C, Pasquier L, Toutain A, Phan Dinh Tuy F, et al: Large spectrum of lissencephaly and pachygyria phenotypes resulting from de novo missense mutations in tubulin alpha 1A (TUBA1A). Human Mutation. 2007, 28 (11): 1055-1064. 10.1002/humu.20572.

    Article  CAS  PubMed  Google Scholar 

  27. Tian G, Kong X-P, Jaglin XH, Chelly J, Keays D, Cowan NJ: A Pachygyria-causing {alpha}-Tubulin Mutation Results in Inefficient Cycling with CCT and a Deficient Interaction with TBCB. Mol Biol Cell. 2008, 19 (3): 1152-1161. 10.1091/mbc.E07-09-0861.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  28. Ben-Zeev B, Levy-Nissenbaum E, Lahat H, Anikster Y, Shinar Y, Brand N, Gross-Tzur V, MacGregor D, Sidi R, Kleta R, et al: Megalencephalic leukoencephalopathy with subcortical cysts; a founder effect in Israeli patients and a higher than expected carrier rate among Libyan Jews. Human Genetics. 2002, 111 (2): 214-218. 10.1007/s00439-002-0770-y.

    Article  CAS  PubMed  Google Scholar 

  29. Leegwater PAJ, Boor PKI, Yuan BQ, Steen van der J, Visser A, Könst AAM, Oudejans CBM, Schutgens RBM, Pronk JC, van der K, M S: Identification of novel mutations in MLC1 responsible for megalencephalic leukoencephalopathy with subcortical cysts. Human Genetics. 2002, 110 (3): 279-283. 10.1007/s00439-002-0682-x.

    Article  CAS  PubMed  Google Scholar 

  30. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3 (7): 10.1186/gb-2002-3-7-research0034.

Download references

Acknowledgements

The study was supported by The Danish Research Council FSS (A.L.N.), the Danish Trans-Disciplinary Research Programme: "Comparative Genomics Focusing on Man and Pig" (A.L.J.), The Lundbeck Foundation (A.L.N.), a Hallas Møller fellowship from the NovoNordisk Foundation to (A.L.N.), and the Faculty of Health Science, Aarhus University Graduate Scholarship (K.B.N.).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Karsten B Nielsen.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

KBN carried out the molecular genetic and bioinformatics studies and drafted the manuscript. MK carried out the array experiments and participated in the array data analysis. IEH, ALJ, and ALN participated in the design of the study, in the evaluation of the results, and in drafting and finalizing the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

13104_2010_531_MOESM1_ESM.PDF

Additional file 1:Lissencephalic and gyrencephalic cortical structure during pig embryogenesis. Gross appearance of the pig brain during development at E60 and E80. (PDF 42 KB)

13104_2010_531_MOESM2_ESM.PDF

Additional file 2:qRT-PCR analysis of selected genes representing all three trimesters of pig cortical brain development. qRT-PCR analysis of a group of selected genes at cortical developmental times E40, E60, E80, E100, and E115. The examined time points are indicated. The expression levels were normalized to GAPDH, Beta-actin and 18S rRNA expression using the geNorm program. (PDF 520 KB)

13104_2010_531_MOESM3_ESM.PDF

Additional file 3:Inter-Chip Comparison. The Scaling factor (SF) measures the differences in signal intensity between the individual chips. Affymetrix recommends that the SF value should be within 3-fold of one another. The number of genes expressed in cortex tissue is ~66% of the total gene number represented on the chip. (PDF 49 KB)

13104_2010_531_MOESM4_ESM.PDF

Additional file 4:Verification of chip data by qRT-PCR. Comparison of fold change in expression between E60 and E80 obtained from microarray and qRT-PCR analysis. (PDF 43 KB)

Authors’ original submitted files for images

Rights and permissions

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.

Reprints and permissions

About this article

Cite this article

Nielsen, K.B., Kruhøffer, M., Holm, I.E. et al. 1Identification of genes differentially expressed in the embryonic pig cerebral cortex before and after appearance of gyration. BMC Res Notes 3, 127 (2010). https://0-doi-org.brum.beds.ac.uk/10.1186/1756-0500-3-127

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/1756-0500-3-127

Keywords