Transcriptome profiling in engrailed-2 mutant mice reveals common molecular pathways associated with autism spectrum disorders
Molecular Autism volume 4, Article number: 51 (2013)
Transcriptome analysis has been used in autism spectrum disorder (ASD) to unravel common pathogenic pathways based on the assumption that distinct rare genetic variants or epigenetic modifications affect common biological pathways. To unravel recurrent ASD-related neuropathological mechanisms, we took advantage of the En2-/- mouse model and performed transcriptome profiling on cerebellar and hippocampal adult tissues.
Cerebellar and hippocampal tissue samples from three En2-/- and wild type (WT) littermate mice were assessed for differential gene expression using microarray hybridization followed by RankProd analysis. To identify functional categories overrepresented in the differentially expressed genes, we used integrated gene-network analysis, gene ontology enrichment and mouse phenotype ontology analysis. Furthermore, we performed direct enrichment analysis of ASD-associated genes from the SFARI repository in our differentially expressed genes.
Given the limited number of animals used in the study, we used permissive criteria and identified 842 differentially expressed genes in En2-/- cerebellum and 862 in the En2-/- hippocampus. Our functional analysis revealed that the molecular signature of En2-/- cerebellum and hippocampus shares convergent pathological pathways with ASD, including abnormal synaptic transmission, altered developmental processes and increased immune response. Furthermore, when directly compared to the repository of the SFARI database, our differentially expressed genes in the hippocampus showed enrichment of ASD-associated genes significantly higher than previously reported. qPCR was performed for representative genes to confirm relative transcript levels compared to those detected in microarrays.
Despite the limited number of animals used in the study, our bioinformatic analysis indicates the En2-/- mouse is a valuable tool for investigating molecular alterations related to ASD.
Autism spectrum disorder (ASD) defines a complex group of neurodevelopmental disabilities characterized by a wide range of impairments in social and communicative skills, stereotyped behaviors, and restricted mental flexibility . The neurodevelopmental and neuroanatomical bases of ASD have been confirmed by a number of clinical, neuroimaging and neuropathological studies [1–3]. The most evident abnormality in ASD consists in an early (perinatal) brain overgrowth followed by an arrest of growth during the first year of age . Neuropathological studies on post-mortem samples from ASD patients also showed a number of cellular and cytoarchitectural abnormalities at the level of the cerebral cortex, cerebellum, amygdala and forebrain limbic structures. Anomalies in the cerebellum are the most reproducible neuropathological alterations in ASD patients [3, 5].
A large series of evidence clearly indicates that neuropathological features and behavioral deficits of ASD have a primarily genetic origin. However, the etiology of ASD remains essentially unknown [6, 7]. Transcriptome analysis has also been used to unravel common pathogenic pathways based on the assumption that distinct rare genetic variants or epigenetic modifications affect common biological pathways dysregulated in ASD . Several studies have analyzed genome-wide expression profiles of ASD patients using lymphoblastoid cell lines and blood samples, supporting upregulation of immune genes and downregulation of neurodevelopmental genes as key players in the pathogenesis of ASD (see  for a review). Recently, gene co-expression network analysis of autistic brain areas revealed defects in cortical patterning and an enrichment of differentially expressed genes associated with ASD .
The homeobox-containing transcription factor engrailed-2 (En2) is crucially involved in the regionalization, patterning and neuronal differentiation of the midbrain and hindbrain [10–15]. Human studies indicated association of two intronic single-nucleotide polymorphisms (SNPs) in the human engrailed-2 (EN2) gene with ASD [16, 17]. Furthermore, the ASD associated A-C haplotype markedly affected EN2 promoter activity when tested with a luciferase reporter assay in rat, mouse and human cell lines . A recent study of the epigenetic evaluation of EN2 in post-mortem cerebellar samples from autistic patients indicated a persistent upregulation of this homeobox gene induced by epigenetic abnormalities in histone methylation patterns that may contribute to Purkinje cell loss in some individuals with autism .
Mice lacking the homeobox domain of En2 (En2hd/hd mice; , referred to as En2-/-) have been proposed as a model for ASD, due to their complex neuroanatomical and behavioral phenotype. En2-/- mice display cerebellar hypoplasia, including a reduced number of Purkinje cells, and a defect in the antero-posterior pattern of cerebellar foliation [20–23]. The behavior of En2-/- mice is also reminiscent of some features of ASD individuals. Deficits in social behaviors were detected in En2-/- mice, including decreased play and reduced social interactions; locomotor impairment, as well as defective spatial learning and memory, was also reported in these mice [24–26]. Furthermore, we reported dysfunctions in GABAergic interneurons in adult En2-/- mice and demonstrated engrailed protein expression in specific subpopulations of adult hippocampal and cortical interneurons .
To unravel recurrent ASD-related neuropathological mechanisms, we took advantage of the En2-/- mouse model and performed genome-wide expression profiling on cerebellar and hippocampal adult tissues. Our transcriptome analysis of the cerebellum and hippocampus of En2-/- mice suggests convergent pathological pathways with ASD, including abnormal synaptic transmission and increased immune response. Furthermore, we provide evidence for a significant enrichment of differentially expressed genes associated to ASD in this mouse model of the disease.
Experiments were conducted in conformity with the European Communities Council Directive of 24 November 1986 (86/609/EEC) and were approved by the Italian Ministry of Health and Ethics Committee of the University of Trento. Animals were housed in a 12 hr light/dark cycle with food and water available ad libitum. All surgery was performed under chloral hydrate anesthesia, and all efforts were made to minimize suffering. The generation of En2-/- mice was previously described . The original En2 mutants (mixed 129Sv x C57BL/6 and outbred genetic background) were crossed at least five times into a C57BL/6 background. Heterozygous matings (En2+/- x En2+/-) were used to generate the En2+/+ (wild type, WT) and En2-/- littermates used in this study. PCR genotyping was performed according to the protocol available at the Jackson Laboratory website (http://www.jax.org; mouse strain En2tm1Alj). WT and En2-/- age-matched adult (3 to 5 months old; weight = 25 to 35 g) littermates mice of both sexes were used.
RNAs from dissected hippocampi and cerebella from three adult mice for each genotype were purified using standard column purification according to the manufacturer’s protocol (RNAeasy Mini Kit, Qiagen, USA). RNA quality was analyzed by microfluidic gel electrophoresis on RNA 6000 NanoChips using the Agilent 2100 Bioanalyzer. Only RNA with a high (>9) RNA integrity number was selected and used for subsequent retro-transcription, labeling and array hybridization according to Agilent protocols. Mouse gene expression arrays (Agilent 4X44K slides) were hybridized and scanned with the Agilent microarray station.
Intensity values were processed with Agi4x44PreProcess (http://bioconductor.org/packages/2.12/bioc/html/Agi4x44PreProcess.html) using default parameters to remove low-quality probes. Signals were then normalized by means of the quantile normalization method. To evaluate differential expression, we used RankProd (http://www.bioconductor.org/packages/2.11/bioc/html/RankProd.html) . RankProd utilizes the Rank Product (RP) non-parametric method  to identify up- or downregulated genes. The RP is equivalent to calculating the geometric mean rank with a statistical method (average rank) that is slightly more sensitive to outlier data and puts a higher premium on consistency between the ranks in various lists. To assess for functional categories overrepresented in the differentially expressed genes, we used DAVID (http://david.abcc.ncifcrf.gov) and Ingenuity Pathway Analysis (Ingenuity Systems, Inc., USA). To focus the functional analysis on brain expressed genes we used, as background for our functional analyses, a list of tissue specific ‘expressed genes’ for both the cerebellum and the hippocampus. Our ‘expressed genes’ lists were obtained by filtering the genes by the normalized expression values and excluding the ones with the lowest expression levels (<10th percentile), and include 13,652 genes for the cerebellum and 13,141 for the hippocampus. The hypergeometric test and the Student’s t-test were computed with R (http://www.r-project.org).
Total RNAs were extracted by Trizol™ reagent (Invitrogen Life Technologies, USA) from dissected hippocampi and cerebella from four WT and four En2-/- adult mice. RNAs were DNAse-treated and purified with RNeasy Mini Kit (Qiagen, USA). cDNA was synthesized from pooled RNAs (2 μg) using the SuperScript™ VILO™ (Invitrogen Life Technologies, USA) according to the manufacturer’s instructions. Individual PCR reactions were conducted in a volume of 20 μl using the KAPA FAST SYBR qPCR kit (KAPABiosystems, USA) according to manufacturer’s instructions. Mouse mitochondrial ribosomal protein L41 (Mrpl41) was used as a standard for quantification as previously shown . Primers (MWG, Germany) were designed on different exons to avoid amplification of genomic DNA. A list of primer sequences is reported in Additional file 1. Each PCR cycle consisted of denaturation for 10 s at 94°C, annealing for 20 s at 60°C and extension for 30 s at 72°C. The fluorescence intensity of SYBR green I was read and acquired at 72°C after completion of the extension step of each cycle. PCR conditions for individual primer sets were optimized by varying template cDNA and primer concentration in order to obtain a single PCR product and amplification efficiency >90%. Relative expression values were calculated using the Pfaffl method .
Differential gene expression in cerebellum and hippocampus of En2-/- mice
The cerebellum of En2-/- mice shows Purkinje cell loss and structural abnormalities resembling the neuropathological features observed in ASD patients [20, 21, 23]. To identify genes and pathways altered in En2-/- mice, we initially performed transcriptome profiling of the whole cerebellar tissue. En2-/- and WT cerebella from adult mice were assessed for differential gene expression by microarray and bioinformatical analysis (see Methods). We found 842 differentially expressed genes in the cerebellum of En2-/- mice compared to their control littermates. Among these, 407 and 435 were up- and downregulated, respectively. Alterations in limbic structures have also been shown in the temporal lobes of autistic patients. The main abnormalities were shown in the superior temporal sulcus and the ventral-basal temporal region, both of which are involved in decoding social stimuli and therefore are associated with the social deficits [32–34]. Most importantly, we previously showed anatomical defects in the En2-/- hippocampus that might contribute to the behavioral deficits displayed by these mutants . We therefore combined the cerebellar gene expression profile to that of the hippocampus. We found 862 differentially expressed genes in the hippocampus, among those 378 were upregulated and 484 were downregulated in En2-/- mice compared to their littermate controls. Additional file 2 shows the entire list of genes differentially expressed in the cerebellum and hippocampus of En2-/- mice with the differential expression P value and the percentage of false prediction (pfp) value calculated by RankProd. Given the profound structural and cytoarchitectural phenotype of the En2-/- cerebellum and the variability of the phenotype among individuals, we chose to be more permissive and include genes with smaller differential expression fold changes applying a moderate cutoff to the uncorrected P value [see Additional file 2]. Differentially expressed genes, which are common in the hippocampus and the cerebellum, are summarized in Figure 1 and listed in Additional file 3. Remarkably, En2 was not among the differentially expressed genes. The microarray probe for En2 (A_51_P397876) was designed on the 3′ untranslated region of the gene that is not deleted in the En2-/- locus . Our data confirm the previous reports indicating the presence of a residual transcript from the En2-/- locus . Furthermore, real-time quantitative PCR (qPCR) analysis of En2 expression using homeobox specific primers revealed no expression of the full-length gene in En2-/- mice [see Additional file 4].
To explore the biological processes most relevant to the phenotype of the En2-/- mice, the differentially expressed genes were analyzed through integrated gene-network analysis using the curated Ingenuity Pathway Analysis (IPA) database and the Database for Annotation, Visualization and Integrated Discovery (DAVID). Additional file 5 shows the most significantly enriched disease and cellular function categories obtained with IPA and gene ontology analysis. Functional categories included increased seizure and decreased neurotransmission release in the hippocampus, and decreased cancer-related diseases and development of lymphocytes as cellular function in the cerebellum. To identify statistically significant over-representation of key neurobiological processes, functional annotation analysis was performed with DAVID using a tissue specific list of ‘expressed genes’ as background (see Methods). To verify the tissue-expression pattern of the samples, we first classified the differentially expressed genes based on their tissue expression (P <0.05, calculated using Benjamini multiple testing correction). We only found significant tissue expression terms for En2-/- hippocampus showing significant enrichment in expression of genes related to ‘brain cortex’, ‘brain’, ‘hypothalamus’, ‘eye’ and ‘hippocampus’ [see Additional file 6]. We then analyzed the functional annotation using gene sets from the gene ontology (GO) public databases and our ‘expressed genes’ list as background. Figure 2 shows all GO terms enriched in the differentially expressed genes of En2-/- cerebellum with the gene counts and the relative P value (P <0.05, calculated using Benjamini multiple testing correction). Among the most represented functional categories were several terms related to the major histocompatibility complex (MHC)-mediated immunity and immune response. The GO terms related to the En2-/- hippocampus are shown in Figure 3 with the gene counts and the relative P value (P <0.05, calculated using Benjamini multiple testing correction). Of the enriched terms, many related to the cellular components synapse, synaptic vesicle and neuronal projection. Within the biological processes, the most represented terms were ‘neuronal activities’ and ‘calcium mediated signal’. Interestingly, the keyword ‘phosphoprotein’ was significantly enriched, suggesting a generic alteration in the protein phosphorylation state in the En2-/- hippocampus.
To assess the functional consequences of En2 ablation, we analyzed enrichment of the differentially expressed genes based on the mouse phenotype from the Mammalian Phenotype Ontology (MPO) project , using ToppGene . In the En2-/- cerebellum the only significantly enriched mouse phenotype was ‘loss of dopaminergic neurons’, a term associated with the cerebellar phenotype by the role of En2 in midbrain/hindbrain development and dopaminergic neuron survival [37, 38]. For the hippocampus we found significant enrichment for terms related to seizure and altered synaptic transmission (Figure 4). These data are in accordance with our previously described increase in seizure susceptibility and decrease in GABAergic neuron subpopulations in En2-/- mice [27, 30].
Remarkably, when analyzed specifically for ASD-associated genes, the differentially expressed genes showed significant over-representation of known ASD susceptibility genes when compared directly to the repository of the Simons Foundation and Autism Research Initiative (SFARI) database (gene.sfari.org). Since the SFARI database comprises genes related to neurodevelopmental disorders, therefore expressed in the brain, we filtered the SFARI database gene list with our list of ‘expressed genes’ in the cerebellum and in the hippocampus (see Methods for details). We used these SFARI gene lists to calculate enrichment with our differentially expressed genes. Table 1 shows the list of the ASD-associated differentially expressed genes for the cerebellum and the hippocampus. To compute significant enrichment between our differentially expressed genes and the SFARI genes, we employed the hypergeometric test. The statistical analysis indicated significant enrichment only for the hippocampus (P <0.05), whereas no significant enrichment was observed for the cerebellum.
To compare our findings with the three major genome-wide expression studies on ASD brain tissue [9, 39, 40], we matched the publicly available differentially expressed genes with the same ASD-related gene lists that we used for our analysis. To increase accuracy, we computed, for each study, the hypergeometric test and obtained an enrichment P value that we used for direct comparison. Table 2 summarizes the results of the enrichment analyses, separated in tissue-specific groups, and the comparison between the studies. For Voineagu et al. we re-analyzed the cerebellum data using the GEO2R tool with default parameters and used these results to evaluate the correspondence with our study. The results for cerebellum show a significant enrichment with ASD-associated genes only for the Voineagu et al. study (4.75% enrichment, P = 0.0343). Our study, however, was the only one to display significant enrichment of ASD-related genes (4.24% enrichment, p = 0.0265) in the limbic regions.
To validate microarray findings with qPCR analysis, we selected eight and eighteen representative genes from the cerebellum and the hippocampus differentially expressed gene lists, respectively. The selected genes reported differential expression values in microarray experiments ranging between 0.41 fold decrease to 2.5 fold increase. With qPCR, five of the eight selected genes showed statistically significant differential expression in the En2-/- cerebellum [glutamate receptor, metabotropic 5 (Grm5) P = 0.031; neurexin III (Nrxn3) P = 0.032; neuropilin 2 (Nrp2) P = 0.009; v-erb-a erythroblastic leukemia viral oncogene homolog 4 (Erbb4) P = 0.003; calcium channel, voltage-dependent, T type, alpha 1G subunit (Cacna1g) P = 0.003] (Figure 5a). Fifteen genes of the eighteen selected were instead significantly changed in the En2-/- hippocampus [netrin G1 (Ntng1) P = 0.008; Parkinson disease (autosomal recessive, juvenile) 2, parkin, (Park2) P = 0.004; gamma-aminobutyric acid (GABA) A receptor, subunit alpha 4 (Gabra4) P = 0.007; early growth response 2 (Egr2) P = 0.007; sodium channel, voltage-gated, type I, alpha (Scn1a) P = 0.001; neurotrophic tyrosine kinase, receptor, type 3 (Ntrk3) P = 0.0003; phospholipase C, beta 1 (Plcb1) P = 0.0001; cortactin (Cttn) P = 0.0005; synapsin I (Syn1) P = 0.007; fragile X mental retardation, autosomal homolog 2 (Fxr2) P = 0.0004; protein kinase C, beta (Prkcb) P = 0.001; gamma-aminobutyric acid (GABA) B receptor, 1, (Gabbr1) P = 0.03; neuroligin 2 (Nlgn2) P = 0.001; glycogen synthase kinase 3 beta (Gsk3b) P = 0.003; GNAS (guanine nucleotide binding protein, alpha stimulating) complex locus (Gnas) P = 0.01] (Figure 5b). Except for Egr2, in all tested genes the expression differences reported by qPCR correlated with the microarray data (Pearson r = 0.84, P < 0.00005).
To date, more than 500 autism-associated genes have been identified (SFARI Gene; gene.sfari.org; updated mar/2013); yet the etiology of ASD remains essentially unknown [6, 7]. The significance of animal models in ASD research has been widely recognized as important for unraveling the molecular, cellular, anatomical, electrophysiological and behavioral consequences of gene dysfunction in ASD. Here, we present a transcriptome analysis in a mouse model of ASD of two brain areas, the cerebellum and the hippocampus, areas that are profoundly affected in ASD patients. Despite the small number of samples used for the microarray analysis and the sample gender heterogeneity, the low genetic variance among individuals allowed a reasonable statistical power for our bioinformatic analysis. Our study revealed that the molecular signature of these two brain regions shares convergent pathological pathways with ASD, including abnormal synaptic transmission and increased immune response. Furthermore, when directly compared to the repository of the SFARI database (gene.sfari.org), our differentially expressed genes in the hippocampus show an enrichment of ASD-associated genes significantly higher than previously reported .
Transcriptome analysis has been employed to unravel common pathways based on the assumption that the core phenotypes of ASD may be caused by convergent molecular mechanisms . Several studies have analyzed genome-wide expression profiles of lymphoblastoid cell lines and blood samples from ASD patients, pointing to an upregulation of immune genes as key mechanisms in the pathogenesis of ASD . Despite the limited source of brain tissue samples from ASD cases and the technical restrictions, studies of ASD brain transcriptome are emerging as strategic for uncovering functionally relevant alterations in gene expression. A previous microarray study found alterations of glutamatergic neurotransmission in ASD cerebellum , and expression profiles from ASD patient temporal cortices showed upregulation of genes involved in innate immune response and downregulation of several neurodevelopmental genes . Moreover, the transcriptome profiles from three different brain regions (frontal cortex, temporal cortex and cerebellum) of nineteen autism cases and seventeen controls were investigated recently using classical differential expression analysis and a network-based approach . These analyses showed upregulation of genes involved in immune response and downregulation of genes involved in synaptic function and vesicular transport . Our results are in accordance with these findings. Using gene ontology enrichment, integrated gene-network analysis and mouse phenotypes analysis, we report significantly enriched functions and pathways that were previously associated to ASD . In detail, we found increased immune response and major histocompatibility complex-related immunity in the En2-/- cerebellum; decreased and abnormal neurotransmission and increased seizures in the En2-/- hippocampus [see Additional file 5 for details]. Moreover, by direct comparison with the SFARI repository of ASD-related genes, we show that the gene expression changes observed in the En2-/- hippocampus were significantly enriched in ASD-related genes. Furthermore, the proportion of ASD-associated genes enrichment in En2-/- hippocampus was significantly higher than previous studies (Table 2) when compared with Voineagu et al., likely the most comprehensive transcriptome study of ASD post-mortem brain to date. In the case of the cerebellum, in contrast to Voineagu et al. we did not find significant enrichment of ASD-associated genes in En2-/- mice. Such difference could be the result of the complex structural and cytoarchitectural abnormalities in En2-/- cerebellum [20, 21] and the consequent phenotypical variability, or could simply reflect differences between mouse and human phenotypes, as the incidence of cerebellar hypoplasia was not reported in the diagnostic criteria used in the study . Remarkably, EN2 was among the differentially expressed genes found in Voineagu et al., confirming our evidence about the role of En2 in the neuropathology of ASD, and in anterior brain structures .
Among the differentially expressed genes, Grm5, Nrxn3 and Scn1a are of particular interest for ASD. Grm5 encodes mGluR5, a G-protein coupled receptor for the neurotransmitter glutamate . In a recent study, mGluR5 has been shown to participate in the pathogenesis of fragile X syndrome (FXS) while genetic downregulation of Grm5 was able to compensate for some of the symptoms in a mouse model of FXS . Furthermore, Grm5 was shown to be downregulated in hippocampal neurons lacking Shank3, another ASD-associated gene . These data support a central role for Grm5 in neurobiological pathways related to ASD pathogenesis. Our results show an increased expression of Grm5 in the cerebellum of En2-/- mice, suggesting a role of Grm5 in the cerebellar phenotype of these mice. The contribution of Grm5 and its interaction with Fmr1 in the En2-/- hippocampus remains to be established and could open new perspective of pharmacological and genetic rescue of the ASD-related phenotype of these mice.
Nrxn3 encodes neuronal adhesion proteins of the Neurexin (NRXN) family. NRXNs are presynaptic cell adhesion proteins that form trans-synaptic complexes with their postsynaptic counterpart neuroligins (NLGNs) and have important roles in synapse development and function . Recently, a report of hemizygous and de novo deletions involving NRXN3 in ASD families provided strong support for a causative link between the loss of NRXN3 and the development of ASD . Our results of an increased expression of Nrxn3 in the cerebellum suggest alterations in Purkinje cell synaptic formation, where NRXNs have been shown to participate to the formation of glutamatergic synapses through interaction with Cerebellin 1 precursor protein (also downregulated in the En2-/- cerebellum) and GluR∂2 [48, 49].
Scn1a encodes the voltage-gated sodium channel alpha subunit. De novo null mutations in SCN1A result in severe myoclonic epilepsy of infancy . SCN1A mutations have been associated to a number of neurological disorders, including generalized epilepsy with febrile seizures plus, Dravet syndrome, borderline myoclonic epilepsy in infancy, intractable childhood epilepsy with generalized tonic-clonic seizures, familial hemiplegic migraine, and a number of cryptogenic focal and generalized epilepsies. Recently, de novo mutations in SCN1A have been associated with ASD , and a report of a recognized mutation in SCN1A suggests a wide phenotypic variation of the gene mutations causing a variety of neurologic disorders, including ASD . In mice, heterozygous loss-of-function mutation in Scn1a (Scn1a+/-), reproduces several of the symptoms associated to the human mutation, such as thermally induced and spontaneous seizures, premature death, ataxia and sleep disorder [53, 54]. Scn1a+/- mice show both cognitive deficits and autistic traits that are caused by impaired GABAergic neurotransmission and can be rescued by drug treatment. Scn1a down-regulation in the En2-/- hippocampus could contribute to the abnormal excitability and altered GABAergic neurotransmission shown in these mice by our previous studies [27, 30]. Pharmacological rescue of the hippocampal phenotype in the En2-/- with GABAergic drugs is currently under investigation.
Anomalies in the cerebellum are the most reproducible neuropathological alterations in ASD patients. Several cerebellar abnormalities have been observed in mouse models of both En2 gain- and loss-of-function. Ectopic overexpression of En2 in Purkinje cells during late embryonic and postnatal cerebellar development results in reduced cerebellar volume and loss of Purkinje cells and other cerebellar neurons [55, 56]. Interestingly, En2 knock-out causes defective cerebellar patterning, reduced Purkinje cell number and abnormal dendritic foliation [10, 57], indicating that alterations in En2 expression levels during development cause similar phenotypes. Furthermore, deficits in social behaviors as well as defective spatial learning and memory were also reported in En2-/- mice [24–26]. A recent epigenetic analysis of EN2 promoter methylation in the cerebellum of ASD individuals indicated hypermethylation of the promoter region and persistent upregulation of the gene. The authors report that promoter hypermethylation is normally associated with a decrease in gene expression and suggest the possibility of a developmental mechanism intended to support downregulation of EN2 during perinatal development . Taken together, this evidence suggests that an overall imbalance in EN2 expression may be relevant for ASD pathogenesis, as it could produce alterations in critical brain functions. Comparable evidence of a similar dosage effect has been reported in the case of mutations of other genes critically involved in gene expression regulation and maintenance of synaptic and neuronal homeostasis, such as MECP2 and ARX[58, 59]. It remains to be established whether En2 overexpressing mice display abnormal behaviors relevant to autism. Microarray data have been produced for En2 overexpressing Purkinje cells on a different platform; however, the results overlap only marginally with the herein reported study .
Using transcriptome analysis, we identified over 800 genes differentially expressed in the cerebellum and hippocampus of En2-/- mice. Despite the small number of samples used and the relatively small statistical power, our study is the first to analyze molecular changes occurring in two brain structures with neuroanatomical alterations relevant to ASD in a mouse model of this disease. Our bioinformatic analysis of the molecular signature of En2-/- cerebellum and hippocampus shows a significant convergence of neurobiological pathways previously linked to ASD pathology in brain samples from ASD patients. Overall, the present study points to a strong impact of transcriptome analysis on mouse models for identifying neurobiological pathways commonly altered when ASD genes are disrupted in a human patient and in a mouse model alike. Furthermore, together with the frequent association of cerebellar neuroanatomical alterations to the neuropathology of ASD, our molecular analysis suggests a contribution also for the hippocampus, where molecular changes relevant to ASD may occur also in human patients. This notion is supported by the consistent enrichment of ASD-related genes in the En2-/- hippocampus compared to the cerebellum and to other similar studies performed on ASD patient tissue samples . Finally, our results confirm the En2-/- mouse model of ASD as a valuable tool for investigating neuroanatomical, behavioral, as well as molecular alterations related to ASD.
Availability of supporting data
The data sets supporting the results of this article are available in the GEO repository, GSE51612.
Autism spectrum disorders
Percentage of false positives
Database for Annotation, Visualization and Integrated Discovery
Fragile X syndrome
Gamma aminobutyric acid
Ingenuity pathway analysis
Major histocompatibility complex
Mammalian phenotype ontology
Simons Foundation and Autism Research Initiative
Rapin I, Tuchman RF: Autism: definition, neurobiology, screening, diagnosis. Pediatr Clin North Am. 2008, 46: 1129-viii
Dicicco-Bloom E, Lord C, Zwaigenbaum L, Courchesne E, Dager SR, Schmitz C, Schultz RT, Crawley J, Young LJ: The developmental neurobiology of autism spectrum disorder. J Neurosci. 2006, 26: 6897-6906. 10.1523/JNEUROSCI.1712-06.2006.
Pardo CA, Eberhart CG: The neurobiology of autism. Brain Pathol. 2007, 17: 434-447. 10.1111/j.1750-3639.2007.00102.x.
Courchesne E: Brain development in autism: early overgrowth followed by premature arrest of growth. Ment Retard Dev Disabil Res Rev. 2004, 10: 106-111. 10.1002/mrdd.20020.
Bailey A, Luthert P, Dean A, Harding B, Janota I, Montgomery M, Rutter M, Lantos P: A clinicopathological study of autism. Brain. 1998, 121: 889-905. 10.1093/brain/121.5.889.
Voineagu I: Gene expression studies in autism: moving from the genome to the transcriptome and beyond. Neurobiol Dis. 2012, 45: 69-75. 10.1016/j.nbd.2011.07.017.
Huguet G, Ey E, Bourgeron T: The genetic landscapes of autism spectrum disorders. Annu Rev Genomics Hum Genet. 2013, 14: 191-213. 10.1146/annurev-genom-091212-153431.
Lintas C, Sacco R, Persico AM: Genome-wide expression studies in autism spectrum disorder, Rett syndrome, and Down syndrome. Neurobiol Dis. 2012, 45: 57-68. 10.1016/j.nbd.2010.11.010.
Voineagu I, Wang X, Johnston P, Lowe JK, Tian Y, Horvath S, Mill J, Cantor RM, Blencowe BJ, Geschwind DH: Transcriptomic analysis of autistic brain reveals convergent molecular pathology. Nature. 2011, 474: 380-384. 10.1038/nature10110.
Cheng Y, Sudarov A, Szulc KU, Sgaier SK, Stephen D, Turnbull DH, Joyner AL: The Engrailed homeobox genes determine the different foliation patterns in the vermis and hemispheres of the mammalian cerebellum. Development. 2010, 137: 519-529. 10.1242/dev.027045.
Gherbassi D, Simon HH: The engrailed transcription factors and the mesencephalic dopaminergic neurons. J Neural Transm Suppl. 2006, 70: 47-55.
Herrup K, Murcia C, Gulden F, Kuemerle B, Bilovocky N: The genetics of early cerebellar development: networks not pathways. Prog Brain Res. 2005, 148: 21-27.
Joyner AL: Engrailed, Wnt and Pax genes regulate midbrain–hindbrain development. Trends Genet. 1996, 12: 15-20. 10.1016/0168-9525(96)81383-7.
Sgaier SK, Lao Z, Villanueva MP, Berenshteyn F, Stephen D, Turnbull RK, Joyner AL: Genetic subdivision of the tectum and cerebellum into functionally related regions based on differential sensitivity to engrailed proteins. Development. 2007, 134: 2325-2335. 10.1242/dev.000620.
Orvis GD, Hartzell AL, Smith JB, Barraza LH, Wilson SL, Szulc KU, Turnbull DH, Joyner AL: The engrailed homeobox genes are required in multiple cell lineages to coordinate sequential formation of fissures and growth of the cerebellum. Dev Biol. 2012, 367: 25-39. 10.1016/j.ydbio.2012.04.018.
Gharani N, Benayed R, Mancuso V, Brzustowicz LM, Millonig JH: Association of the homeobox transcription factor, ENGRAILED 2, 3, with autism spectrum disorder. Mol Psychiatry. 2004, 9: 474-484. 10.1038/sj.mp.4001498.
Benayed R, Gharani N, Rossman I, Mancuso V, Lazar G, Kamdar S, Bruse SE, Tischfield S, Smith BJ, Zimmerman RA, Dicicco-Bloom E, Brzustowicz LM, Millonig JH: Support for the homeobox transcription factor gene ENGRAILED 2 as an autism spectrum disorder susceptibility locus. Am J Hum Genet. 2005, 77: 851-868. 10.1086/497705.
Benayed R, Choi J, Matteson PG, Gharani N, Kamdar S, Brzustowicz LM, Millonig JH: Autism-associated haplotype affects the regulation of the homeobox gene, ENGRAILED 2. Biol Psychiatry. 2009, 66: 911-917. 10.1016/j.biopsych.2009.05.027.
James SJ, Shpyleva S, Melnyk S, Pavliv O, Pogribny IP: Complex epigenetic regulation of Engrailed-2 (EN-2) homeobox gene in the autism cerebellum. Transl Psychiatry. 2013, 3: e232-10.1038/tp.2013.8.
Joyner AL, Herrup K, Auerbach BA, Davis CA, Rossant J: Subtle cerebellar phenotype in mice homozygous for a targeted deletion of the En-2 homeobox. Science. 1991, 251: 1239-1243. 10.1126/science.1672471.
Millen KJ, Wurst W, Herrup K, Joyner AL: Abnormal embryonic cerebellar development and patterning of postnatal foliation in two mouse Engrailed-2 mutants. Development. 1994, 120: 695-706.
Millen KJ, Hui CC, Joyner AL: A role for En-2 and other murine homologues of Drosophila segment polarity genes in regulating positional information in the developing cerebellum. Development. 1995, 121: 3935-3945.
Kuemerle B, Zanjani H, Joyner A, Herrup K: Pattern deformities and cell loss in Engrailed-2 mutant mice suggest two separate patterning events during cerebellar development. J Neurosci. 1997, 17: 7881-7889.
Gerlai R, Millen KJ, Herrup K, Fabien K, Joyner AL, Roder J: Impaired motor learning performance in cerebellar En-2 mutant mice. Behav Neurosci. 1996, 110: 126-133.
Cheh MA, Millonig JH, Roselli LM, Ming X, Jacobsen E, Kamdar S, Wagner GC: En2 knockout mice display neurobehavioral and neurochemical alterations relevant to autism spectrum disorder. Brain Res. 2006, 1116: 166-176. 10.1016/j.brainres.2006.07.086.
Brielmaier J, Matteson PG, Silverman JL, Senerth JM, Kelly S, Genestine M, Millonig JH, Dicicco-Bloom E, Crawley JN: Autism-relevant social abnormalities and cognitive deficits in engrailed-2 knockout mice. PLoS ONE. 2012, 7: e40914-10.1371/journal.pone.0040914.
Sgadò P, Genovesi S, Kalinovsky A, Zunino G, Macchi F, Allegra M, Murenu E, Provenzano G, Tripathi PP, Casarosa S, Joyner AL, Bozzi Y: Loss of GABAergic neurons in the hippocampus and cerebral cortex of Engrailed-2 null mutant mice: Implications for autism spectrum disorders. Exp Neurol. 2013, 247: 496-505.
Hong F, Breitling R, McEntee CW, Wittner BS, Nemhauser JL, Chory J: RankProd: a bioconductor package for detecting differentially expressed genes in meta-analysis. Bioinformatics. 2006, 22: 2825-2827. 10.1093/bioinformatics/btl476.
Breitling R, Armengaud P, Amtmann A, Herzyk P: Rank products: a simple, yet powerful, new method to detect differentially regulated genes in replicated microarray experiments. FEBS Lett. 2004, 573: 83-92. 10.1016/j.febslet.2004.07.055.
Tripathi PP, Sgado P, Scali M, Viaggi C, Casarosa S, Simon HH, Vaglini F, Corsini GU, Bozzi Y: Increased susceptibility to kainic acid-induced seizures in Engrailed-2 knockout mice. Neuroscience. 2009, 159: 842-849. 10.1016/j.neuroscience.2009.01.007.
Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29: e45-10.1093/nar/29.9.e45.
Lauvin M-A, Martineau J, Destrieux C, Andersson F, Bonnet-Brilhault F, Gomot M, El-Hage W, Cottier J-P: Functional morphological imaging of autism spectrum disorders: current position and theories proposed. Diagn Interv Imaging. 2012, 93: 139-147. 10.1016/j.diii.2012.01.007.
Boddaert N, Chabane N, Gervais H, Good CD, Bourgeois M, Plumet M-H, Barthélémy C, Mouren M-C, Artiges E, Samson Y, Brunelle F, Frackowiak RSJ, Zilbovicius M: Superior temporal sulcus anatomical abnormalities in childhood autism: a voxel-based morphometry MRI study. Neuroimage. 2004, 23: 364-369. 10.1016/j.neuroimage.2004.06.016.
Gendry Meresse I, Zilbovicius M, Boddaert N, Robel L, Philippe A, Sfaello I, Laurier L, Brunelle F, Samson Y, Mouren M-C, Chabane N: Autism severity and temporal lobe functional abnormalities. Ann Neurol. 2005, 58: 466-469. 10.1002/ana.20597.
Smith CL, Goldsmith C-AW, Eppig JT: The Mammalian Phenotype Ontology as a tool for annotating, analyzing and comparing phenotypic information. Genome Biol. 2005, 6: R7-10.1186/gb-2005-6-5-p7.
Chen J, Xu H, Aronow BJ, Jegga AG: Improved human disease candidate gene prioritization using mouse phenotype. BMC Bioinformatics. 2007, 8: 392-10.1186/1471-2105-8-392.
Alvarez-Fischer D, Fuchs J, Castagner F, Stettler O, Massiani-Beaudoin O, Moya KL, Bouillot C, Oertel WH, Lombès A, Faigle W, Joshi RL, Hartmann A, Prochiantz A: Engrailed protects mouse midbrain dopaminergic neurons against mitochondrial complex I insults. Nat Neurosci. 2011, 14: 1260-1266. 10.1038/nn.2916.
Sgadò P, Albéri L, Gherbassi D, Galasso SL, Ramakers GMJ, Alavian KN, Smidt MP, Dyck RH, Simon HH: Slow progressive degeneration of nigral dopaminergic neurons in postnatal Engrailed mutant mice. Proc Natl Acad Sci USA. 2006, 103: 15242-15247. 10.1073/pnas.0602116103.
Purcell AE, Jeon OH, Zimmerman AW, Blue ME, Pevsner J: Postmortem brain abnormalities of the glutamate neurotransmitter system in autism. Neurology. 2001, 57: 1618-1628. 10.1212/WNL.57.9.1618.
Garbett K, Ebert PJ, Mitchell A, Lintas C, Manzi B, Mirnics K, Persico AM: Immune transcriptome alterations in the temporal cortex of subjects with autism. Neurobiol Dis. 2008, 30: 303-311. 10.1016/j.nbd.2008.01.012.
Ziats MN, Rennert OM: Expression profiling of autism candidate genes during human brain development implicates central immune signaling pathways. PLoS ONE. 2011, 6: e24691-10.1371/journal.pone.0024691.
Buxbaum JD, Betancur C, Bozdagi O, Dorr NP, Elder GA, Hof PR: Optimizing the phenotyping of rodent ASD models: Enrichment analysis of mouse and human neurobiological phenotypes associated with high-risk autism genes identifies morphological, electrophysiological, neurological, and behavioral features. Mol Autism. 2012, 3: 1-10.1186/2040-2392-3-1.
Bear MF, Huber KM, Warren ST: The mGluR theory of fragile X mental retardation. Trends Neurosci. 2004, 27: 370-377. 10.1016/j.tins.2004.04.009.
Dölen G, Osterweil E, Rao BSS, Smith GB, Auerbach BD, Chattarji S, Bear MF: Correction of fragile X syndrome in mice. Neuron. 2007, 56: 955-962. 10.1016/j.neuron.2007.12.001.
Verpelli C, Dvoretskova E, Vicidomini C, Rossi F, Chiappalone M, Schoen M, Di Stefano B, Mantegazza R, Broccoli V, Böckers TM, Dityatev A, Sala C: Importance of Shank3 Protein in Regulating Metabotropic Glutamate Receptor 5 (mGluR5) Expression and Signaling at Synapses. J Biol Chem. 2011, 286: 34839-34850. 10.1074/jbc.M111.258384.
Siddiqui TJ, Craig AM: Synaptic organizing complexes. Curr Opin Neurobiol. 2011, 21: 132-143. 10.1016/j.conb.2010.08.016.
Vaags AK, Lionel AC, Sato D, Goodenberger M, Stein QP, Curran S, Ogilvie C, Ahn JW, Drmic I, Senman L, Chrysler C, Thompson A, Russell C, Prasad A, Walker S, Pinto D, Marshall CR, Stavropoulos DJ, Zwaigenbaum L, Fernandez BA, Fombonne E, Bolton PF, Collier DA, Hodge JC, Roberts W, Szatmari P, Scherer SW: Rare deletions at the neurexin 3 locus in autism spectrum disorder. Am J Hum Genet. 2012, 90: 133-141. 10.1016/j.ajhg.2011.11.025.
Uemura T, Lee S-J, Yasumura M, Takeuchi T, Yoshida T, Ra M, Taguchi R, Sakimura K, Mishina M: Trans-Synaptic Interaction of GluRδ2 and Neurexin through Cbln1 Mediates Synapse Formation in the Cerebellum. Cell. 2010, 141: 1068-1079. 10.1016/j.cell.2010.04.035.
Mishina M, Uemura T, Yasumura M, Yoshida T: Molecular mechanism of parallel fiber-Purkinje cell synapse formation. Front Neural Circuits. 2012, 6: 90-
Claes L, Del-Favero J, Ceulemans B, Lagae L, Van Broeckhoven C, De Jonghe P: De novo mutations in the sodium-channel gene SCN1A cause severe myoclonic epilepsy of infancy. Am J Hum Genet. 2001, 68: 1327-1332. 10.1086/320609.
O’Roak BJ, Vives L, Fu W, Egertson JD, Stanaway IB, Phelps IG, Carvill G, Kumar A, Lee C, Ankenman K, Munson J, Hiatt JB, Turner EH, Levy R, O’Day DR, Krumm N, Coe BP, Martin BK, Borenstein E, Nickerson DA, Mefford HC, Doherty D, Akey JM, Bernier R, Eichler EE, Shendure J: Multiplex targeted sequencing identifies recurrently mutated genes in autism spectrum disorders. Science. 2012, 338: 1619-1622. 10.1126/science.1227764.
Frosk P, Mhanni AA, Rafay MF: SCN1A Mutation associated with intractable Myoclonic Epilepsy and Migraine Headache. J Child Neurol. 2013, 28: 389-391. 10.1177/0883073812443309.
Yu FH, Mantegazza M, Westenbroek RE, Robbins CA, Kalume F, Burton KA, Spain WJ, McKnight GS, Scheuer T, Catterall WA: Reduced sodium current in GABAergic interneurons in a mouse model of severe myoclonic epilepsy in infancy. Nat Neurosci. 2006, 9: 1142-1149. 10.1038/nn1754.
Oakley JC, Kalume F, Yu FH, Scheuer T, Catterall WA: Temperature- and age-dependent seizures in a mouse model of severe myoclonic epilepsy in infancy. Proc Natl Acad Sci USA. 2009, 106: 3994-3999. 10.1073/pnas.0813330106.
Jankowski J, Holst MI, Liebig C, Oberdick J, Baader SL: Engrailed-2 negatively regulates the onset of perinatal Purkinje cell differentiation. J Comp Neurol. 2004, 472: 87-99. 10.1002/cne.20059.
Baader SL, Sanlioglu S, Berrebi AS, Parker-Thornburg J, Oberdick J: Ectopic overexpression of engrailed-2 in cerebellar Purkinje cells causes restricted cell loss and retarded external germinal layer development at lobule junctions. J Neurosci. 1998, 18: 1763-1773.
Sillitoe RV, Vogel MW, Joyner AL: Engrailed homeobox genes regulate establishment of the cerebellar afferent circuit map. J Neurosci. 2010, 30: 10015-10024. 10.1523/JNEUROSCI.0653-10.2010.
Olivetti PR, Noebels JL: Interneuron, interrupted: molecular pathogenesis of ARX mutations and X-linked infantile spasms. Curr Opin Neurobiol. 2012, 22: 859-865. 10.1016/j.conb.2012.04.006.
Na ES, Nelson ED, Kavalali ET, Monteggia LM: The impact of MeCP2 loss- or gain-of-function on synaptic plasticity. Neuropsychopharmacology. 2013, 38: 212-219. 10.1038/npp.2012.116.
Holst MI, Maercker C, Pintea B, Masseroli M, Liebig C, Jankowski J, Miething A, Martini J, Schwaller B, Oberdick J, Schilling K, Baader SL: Engrailed-2 regulates genes related to vesicle formation and transport in cerebellar Purkinje cells. Mol Cell Neurosci. 2008, 38: 495-504. 10.1016/j.mcn.2008.04.010.
P.S. is supported by Provincia Autonoma di Trento and the European Community’s FP7/2007-2013 under grant agreement Marie Curie FP7 - PCOFUND-GA-2008- 226070 ‘progetto Trentino’, project EnCort. This work was funded by the Italian Ministry of University and Research (PRIN 2008 grant # 200894SYW2_002 and PRIN 2010–2011 grant # 2010N8PBAA_002 to Y.B.) and the University of Trento (CIBIO start-up grant to S.C. and Y.B.). We thank Andrea Messina, Federico Vaggi and Tommaso Schiavinotto for helpful discussions, and Patrizia Paoli for administrative support.
The authors declare that they have no competing interests.
PS, GP and YB conceived the study. PS, GP, VA, GZ and SG performed microarray and qPCR experiments. PS, GP and ED performed bioinformatic analysis. PS, GP and YB analyzed data. PS and YB wrote the paper. PS, SC and YB provided funding. All authors read and approved the final manuscript.
Paola Sgadò, Giovanni Provenzano contributed equally to this work.
Electronic supplementary material
Additional file 2: Differentially expressed gene in the cerebellum and in the hippocampus of En2-/- mice. Tables showing differentially expressed genes in the En2-/- (a) cerebellum and (b) hippocampus with fold change expression compared to WT littermates, differential expression P value and percentage of false prediction (pfp) calculated with RankProd. P value cut-off P <0.01 for the cerebellum and P <0.005 for the hippocampus. (XLS 340 KB)
Additional file 3: Genes differentially expressed in both En2-/-cerebellum and hippocampus. Genes commonly regulated in the En2-/- cerebellum and hippocampus with differential expression P value and fold change. (XLS 74 KB)
Additional file 4: En2 full-length expression. Quantitative PCR analysis of En2 full-length expression in the WT and En2-/- cerebellum and hippocampus. Values are expressed as En2/L41 comparative quantitation ratios normalized on the expression of WT in the cerebellum (mean ± s.e.m of three replicates from pools of three animals per genotype; P <0.01, Student’s t-test, WT versus En2-/-). (PDF 53 KB)
Additional file 5: Integrated gene-network analysis and gene ontology enrichment. Tables a-f show the ten most significant disease and functional annotations obtained with Ingenuity Pathways Analysis (IPA) for En2-/- cerebellum (a-c) and En2-/- hippocampus (d-f) differentially expressed genes. Functional categories and annotations are shown with their P value, the predicted activation state, the activation z-score, the associated genes and the number of enriched genes. Activation z-score >2 or < -2 indicates significantly increased or decreased annotations. Tables g-j show gene ontology analysis details for En2-/- cerebellum (g, h) and En2-/- hippocampus (i, j) differentially expressed genes. Gene ontology annotations are shown with gene counts and corrected P values. Significant disease and functions for up- (c, f, h) and downregulated (b, e, j) genes are shown separately. (XLS 90 KB)
Additional file 6: Tissue expression of differentially expressed genes in En2-/- hippocampus. Table showing gene ontology tissue expression annotations for differentially expressed genes in the En2-/- hippocampus. No significant tissue expression annotations were found for differentially expressed genes in En2-/- cerebellum. (XLS 35 KB)
About this article
Cite this article
Sgadò, P., Provenzano, G., Dassi, E. et al. Transcriptome profiling in engrailed-2 mutant mice reveals common molecular pathways associated with autism spectrum disorders. Molecular Autism 4, 51 (2013). https://doi.org/10.1186/2040-2392-4-51