Gene expression analysis in Fmr1KO mice identifies an immunological signature in brain tissue and mGluR5-related signaling in primary neuronal cultures
Molecular Autism volume 6, Article number: 66 (2015)
Fragile X syndrome (FXS) is a neurodevelopmental disorder whose biochemical manifestations involve dysregulation of mGluR5-dependent pathways, which are widely modeled using cultured neurons. In vitro phenotypes in cultured neurons using standard morphological, functional, and chemical approaches have demonstrated considerable variability. Here, we study transcriptomes obtained in situ in the intact brain tissues of a murine model of FXS to see how they reflect the in vitro state.
We used genome-wide mRNA expression profiling as a robust characterization tool for studying differentially expressed pathways in fragile X mental retardation 1 (Fmr1) knockout (KO) and wild-type (WT) murine primary neuronal cultures and in embryonic hippocampal and cortical murine tissue. To study the developmental trajectory and to relate mouse model data to human data, we used an expression map of human development to plot murine differentially expressed genes in KO/WT cultures and brain.
We found that transcriptomes from cell cultures showed a stronger signature of Fmr1KO than whole tissue transcriptomes. We observed an over-representation of immunological signaling pathways in embryonic Fmr1KO cortical and hippocampal tissues and over-represented mGluR5-downstream signaling pathways in Fmr1KO cortical and hippocampal primary cultures. Genes whose expression was up-regulated in Fmr1KO murine cultures tended to peak early in human development, whereas differentially expressed genes in embryonic cortical and hippocampal tissues clustered with genes expressed later in human development.
The transcriptional profile in brain tissues primarily centered on immunological mechanisms, whereas the profiles from cell cultures showed defects in neuronal activity. We speculate that the isolation and culturing of neurons caused a shift in neurological transcriptome towards a “juvenile” or “de-differentiated” state. Moreover, cultured neurons lack the close coupling with glia that might be responsible for the immunological phenotype in the intact brain. Our results suggest that cultured cells may recapitulate an early phase of the disease, which is also less obscured with a consequent “immunological” phenotype and in vivo compensatory mechanisms observed in the embryonic brain. Together, these results suggest that the transcriptome of cultured primary neuronal cells, in comparison to whole brain tissue, more robustly demonstrated the difference between Fmr1KO and WT mice and might reveal a molecular phenotype, which is typically hidden by compensatory mechanisms present in vivo. Moreover, cultures might be useful for investigating the perturbed pathways in early human brain development and genes previously implicated in autism.
Fragile X syndrome (FXS) is the most common inherited form of mental retardation and the most common genetic cause of autism. It is caused by loss-of-function mutations in the gene fragile X mental retardation 1 (Fmr1) and a consequent loss of its product, fragile X mental retardation protein (FMRP) [1–3]. The pathophysiology of FXS involves dysregulation of numerous pathways. Recent studies have begun to provide insights into the biology of FMRP, and the results converged on metabotropic glutamate receptor (mGluR) signaling theory [1, 2]. Huber et al. discovered that a form of synaptic plasticity, mGluR-dependent long-term synaptic depression (mGluR-LTD) was exaggerated in Fmr1 knockout mice (Fmr1KO) . This discovery led to the mGluR theory of FXS , which suggests that many of its clinical features are due to exaggerated responses to activation of mGluR5. This theory was validated when multiple FXS phenotypes were rescued in Fmr1KO mice by reducing the production of mGluR5 protein [5–7].
Neurons from Fmr1KO mice and from patients with FXS consistently have increased spine densities, as well as longer spines, reminiscent of immature filopodia [5, 8–11]. Additional synaptic phenotypes, including hyper-connectivity and exaggerated responses to mGluR5 activation, are also found in Fmr1KO mice [4, 5]. These phenotypes can be recapitulated in cultured neurons [12, 13]. They can be corrected by treatment with mGluR5 antagonists in both cultured neurons and mice [5, 13]. Thus, in vitro models derived from Fmr1KO mice are reasonable platforms for modeling synaptic alterations occurring in FXS. Whether these in vitro models also mirror in vivo whole-brain transcriptional regulation has been less studied.
Animal models may be useful for studying the mechanism of disease in FXS, but current methods for characterizing neuronal phenotypes in these models may not be sufficiently robust for high-throughput pre-clinical screening of potential drugs. Furthermore, the molecular phenotypes are themselves heterogeneous. Therefore, identifying stable molecular phenotypes and finding a suitably robust characterization assay is a high priority.
Recent work has shown that RNA expression signatures can identify shared subsets of pathogenic pathways and produce short lists of affected marker genes in a variety of diseased tissues [14, 15]. Done correctly, this process can define a transcriptomic landscape of diseases and tissues and provide directionality with which to measure perturbations “to” and “from” the diseased status. Shared or partially overlapping mechanisms underlie complex phenotypes and may constitute a smaller set of pathways than the number of genetic variants or the genes that contain them. We hypothesized that genome-wide transcriptomic arrays would provide more accurate phenotypic readouts across various tissues and that subsequent comprehensive pathway-level analysis could pinpoint mechanisms relevant to a disease, whereas other morphological methods provide limited phenotypic information.
Our original motivation was to determine how well cell-culture transcriptomes recapitulate differences found in situ in the intact brain tissues of a murine model of FXS. Here, we studied in vitro cultured hippocampal and cortical neurons from Fmr1KO and wild-type (WT) mice and compared them to primary brain tissue. We also compared the location of each of these transcriptomes in the human transcriptomic developmental trajectories to determine how the Fmr1KO state affected the transcriptomic “age” [16, 17] in situ and in culture. We found that this examination of transcriptional changes is able to provide greater resolution on the differences between disease and wild-type conditions.
All animal experiments were conducted in accordance with the rules and regulations of the Institutional Animal Care and Use Committee at the Massachusetts Institute of Technology (MIT). Fmr1KO mice were originally obtained from The Jackson Laboratory (Bar Harbor, ME). The strain was maintained on a C57BL/6 background for at least six generations at MIT. All experiments were performed blind to genotype and were carried out on ex vivo brain tissue following euthanasia. Euthanasia methods are compatible with the recommendations of the Panel on Euthanasia of the American Veterinary Association.
Neuronal culture, tissue isolations, and RNA extraction
Hippocampal and cortical neuron cultures were prepared from C57BL/6 congenic WT and Fmr1KO embryos at 17–18 days post-conception (E17–E18) in parallel and allowed to mature for 14 days in vitro (DIV) as previously described . Five pairs of biological replicates for WT and Fmr1KO were processed on the same day. Cells were seeded at a concentration of 150,000/mL and maintained in Neurobasal media supplemented with B27, penicillin/streptomycin, and glutamax (Invitrogen) according to the manufacturer’s instructions. Cells were harvested on day 14 in PBS and snap frozen using dry ice. The estimated percentage of glial cells at 14 DIV is ~30 %. In parallel to the establishment of primary cultures, ~10 % of hippocampal and cortical tissues were dedicated to tissue analyses and stored in RNAlater for stabilization until further processing.
Four separate sample types were used per genotype (WT or KO): primary cortex (n = 5), primary hippocampus (n = 6), cortical culture (n = 5), and hippocampal culture (n = 5). Use of duplicates resulted in 42 samples. RNA was extracted using the miRNeasy Mini Kit (Qiagen). Eluted RNA was analyzed on a NanoDrop ND-1000 Spectrophotometer and an Agilent 2100 Bioanalyzer (OD260/280 ratio: 1.8–2.2; RNA Integrity Number >8).
Transcriptome profiling using microarrays
A total of 100 ng of RNA was processed using Affymetrix (Affymetrix Inc, Santa Clara, CA) protocols and kits, which were used to generate biotinylated amplified RNA (aRNA) and for hybridization, staining, and scanning of arrays (GeneChip 3’IVT Express Kit and GeneChip Expression Wash, Stain and Scan protocol). Total RNA was reverse-transcribed to synthesize double-stranded cDNA using T7 oligo(dT) primers and then transcribed in vitro into biotin-modified aRNA with IVT Labeling Master Mix. The aRNA was purified, quantified, and fragmented. Fragmented aRNA was hybridized onto Affymetrix Mouse Genome 430 2.0 arrays and scanned on an Affymetrix GeneChip scanner 3000 at 2.5 μm resolution. Microarray data are deposited at the Gene Expression Omnibus database (GSE71034).
Microarray analysis and pathway-level analysis
Expression values were extracted and normalized from .CEL files using the Affy package and the robust multi-array average (RMA) method in R/BioConductor (http://www.bioconductor.org), returning the measured gene expression signal of each microarray gene probe in a logarithmic base 2 scale. Differential gene expression analysis was performed using a linear regression model (lmFit) as implemented in the limma package in R/BioConductor, and significant differentially expressed probes (p < 0.05) were extracted. We used the Database for Annotation, Visualization and Integrated Discovery (DAVID, http://david.abcc.ncifcrf.gov) to identify enriched pathways in differentially expressed genes at a Fisher exact p value threshold (EASE score) less than 0.1, which were annotated with Entrez IDs. We annotated 45,101 probes with 21,141 unique Entrez Gene IDs by selecting the smallest Entrez Gene ID for any probe originally assigned to >1 ID. For pathway-level analysis, we focused on the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Additionally, we performed pathway-level analysis using MetaCore software (Gene Go Inc., St. Joseph, MI). Enrichment analysis consisted of mapping Entrez Gene IDs of differentially expressed genes in culture and in brain onto IDs in entities of built-in functional ontologies represented in MetaCore by process networks and diseases to identify biological processes that were over-represented.
Principal component analysis (PCA) was used to characterize the directions of maximal transcriptomic variance in the whole dataset [19–21]. PCA was performed on an RMA-normalized, ranked, and standardized matrix (mean zero and variance one) of 42 samples (45,101 probes represented on the array). The percentage variance captured by each of the first two principal components (PCs) were 61.3 % by PC1 and 11 % by PC2.
We examined the over-representation of differentially expressed genes in Fmr1KO hippocampal and cortical culture and embryonic brain tissues and three lists of genes associated with defects in cognitive development (Additional file 1: Table S1):
MGI genes associated with the behavioral/neurological phenotype in mice 
A search for the term “MP:0005386” in the JAX lab mouse database (http://www.informatics.jax.org/marker/) returned 3005 unique gene symbols, which are 2622 unique genes mapped to Entrez Gene IDs, out of which 2462 were measured on the Affymetrix Mouse Genome 430 2.0 microarray.
Human genes from SFARI : 705 genes, out of which 645 were mapped to mouse genes
FMRP binding targets as identified by Ascano et al. : 939 genes, out of which 898 were mapped to mouse genes
Human brain transcriptome data processing and integration with mouse data
A spatio-temporal transcriptome of the developing human brain (human brain transcriptome (HBT)) has been described previously , and these data are publicly available in NIH’s Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo) as GSE25219. These samples were profiled on Affymetrix Human Exon 1.0 ST Array, and we used the transcript (gene)-to-sample series matrix for the further analysis. In order to integrate mouse and human gene expression data, we converted 45,101 probes from the Mouse Affymetrix 430 Plus 2.0 array to 21,141 unique Entrez Gene IDs by selecting the smallest Entrez Gene ID for any probe originally assigned to >1 ID. For each differential comparison in mouse data, probe-to-gene unique-ing was done by picking the probe with smallest p value to represent the gene. After downloading human brain transcriptome data from GEO, we converted 17,565 probes to 16,492 unique min human Entrez Gene IDs, 14,653 unique corresponding min homologous mouse Entrez Gene IDs, and 13,830 unique mouse Entrez Gene IDs in common with the Affymetrix Mouse Genome 430 2.0 microarray. We restricted analysis of human data to hippocampus (HIP) and neocortex (NCX) at 15 developmental stages. NCX data included 11 areas collectively referred as the NCX region. In the NCX data for each probe, we computed the sum of coefficients of variance (coefvar) in stages 2–15. For each gene, we selected the probe with the minimal sum of coefvar for probe-to-gene unique-ing. We used the same probes for HIP probe-to-gene unique-ing. In all analyses comparing the two datasets above, we restricted data to 13,830 mouse genes in common between the datasets.
PCA of genes in the sample space was performed on 936 human NCX samples and 82 human HIP samples separately, and k-means clustering (k = 3) was used to generate temporal clusters. Before the PCA of human data, each individual sample, a vector of 13,830 genes, was first normalized to mean zero and variance one. Principal components were invariant under affine transformations.
Gene expression profiles in embryonic brain tissue capture fewer genotype-based differences than primary neuronal cultures
In order to find gene expression differences between the FXS and healthy states, we extracted embryonic hippocampal and cortical tissues from Fmr1KO and WT mice at the age of 17–18 days post-conception (E17–E18). One portion of the tissue was used to create primary neuronal cultures, and the other was used directly in the study of brain tissue.
Total RNA (100 ng) from each sample was extracted, amplified, labeled, and hybridized on Mouse Genome 430 2.0 Arrays (Affymetrix, Santa Clara, CA). Principal component analysis (PCA) was used to characterize the directions of maximal transcriptomic variance in the whole dataset (Fig. 1, each dot signifies a sample). The PCA captured 72.4 % of total variation in the data with the first two principal components (PCs) – 61.3 % captured by PC1 and 11.1 % by PC2. This analysis indicated that the greatest difference was observed between cultured cells and whole tissue, representing the effects of preparation, growth, and differentiation. The second greatest difference was between hippocampal and cortical tissue, representing different brain regions. In addition, Fig. 1 shows that transcriptomic scale differences between two genotypes Fmr1KO and WT are greater in the cultured neurons versus whole embryonic cortical and hippocampal tissue. Additional file 2: Table S2 shows the metrics used to represent the difference between Fmr1KO and WT mice. The PC1-axis centroid distance between KO and WT clusters in both cultures was at least three times greater than the distance between KO and WT in primary brain. In addition, the number of significant differentially expressed genes at an unadjusted p value of 0.05 was almost four times greater in vitro than in whole tissues. Our analysis identified differentially expressed genes (p value <0.05) in Fmr1KO versus WT as follows: 2648 in hippocampal culture, 3372 in cortical culture, 726 in primary hippocampus, and 866 in primary cortex (Additional file 3: Table S3). Fmr1 was the most significantly down-regulated gene in all systems studied (~twofold down-regulation with p value <3.75E−15).
Enriched biological pathways converge on mGluR5-downstream signaling pathways in cultured neurons, while the pathways in embryonic brain tissues center on genes associated with immunological signaling
Differential gene expression analysis was performed using a linear regression model as implemented in the limma library package in R/BioConductor, and significant differentially expressed probes (p < 0.05) were extracted. We used unadjusted p value < 0.05 for significance threshold. Next, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway-enrichment analysis of differentially expressed genes between KO and WT was performed using the Database for Annotation, Visualization and Integrated Discovery (DAVID; Table 1 and Additional file 4: Table S4).
The KEGG pathways that were most significantly enriched with the differentially expressed genes in Fmr1KO cultured cells and isolated embryonic hippocampal and cortical tissues were related to synaptic signaling, immunological response, and cell-cell interactions. Sixteen pathways were significantly enriched in Fmr1KO hippocampal culture compared to WT equivalent cultures (Table 1). Six have been extensively implicated in synaptic signaling. They include “insulin signaling pathway” (p value = 0.001), “long-term potentiation (LTP)” (p value = 0.003), “axon guidance” (p value = 0.009), “gap junction” (p value = 0.017), “neurotrophin signaling pathway” (p value = 0.017) and “mTOR signaling pathway” (p value = 0.08). Eight pathways were significantly enriched in Fmr1KO cortical cultures (Table 1), including “neurotrophin signaling pathway” (p value = 0.006), “MAPK signaling pathway” (p value = 0.016), and “long-term potentiation” (p value = 0.052). Five pathways were significantly enriched in Fmr1KO in situ hippocampus studies relating to immunological response (such as “cytokine-cytokine receptor interaction” (p value = 0.013)) and cell adhesion (“cell adhesion molecules” (p value = 0.017)). Two pathways were significantly enriched in Fmr1KO primary cortex. One of them was “arachidonic acid metabolism” (p value = 0.016), which has a role in inflammation and formation of an important group of inflammatory mediators. Several synaptic mGluR5-coupled signaling pathways perturbed in KO versus WT in cortical and hippocampal neuronal cultures include genes that may have roles in FMRP-regulated mRNA translation at the synapse (Additional file 5: Table S5). Importantly, the signature in embryonic brain tissue was primarily immunological, compared to a primarily synaptic signature in the cultured samples.
To validate and reproduce our findings, we also performed pathway-enrichment analysis on differentially expressed genes in culture and in brain using a separate tool, MetaCore, with similar outcomes. Networks associated with neuronal processes were over-represented in cultures (Development_Neurogenesis_Axonal Guidance, Additional file 6: Figure S1A). In contrast, networks associated with immune disorders were over-represented in brain samples (Autoimmune Diseases, Additional file 6: Figure S1B).
In addition to analyzing each system separately, we also examined the overlap of pathways and commonalities between embryonic brain and cultures and between two different brain regions (cortex and hippocampus; Additional file 7: Table S7). The commonalities between cortical and hippocampal cultures converged on pathways such as “ribosome,” “neurotrophin signaling pathway,” “long-term potentiation,” “endocytosis,” and “spliceosome,” which mostly relate to synaptic signaling. In brain, we observed convergence on the “cell adhesion molecules” pathway between hippocampus and cortex.
We examined the over-representation of differentially expressed genes in Fmr1KO cultures and brain tissues for three gene sets associated with cognitive development defects: (1) Mouse Gene Informatics (MGI) genes associated with behavioral/neurological phenotypes ; (2) human autism candidate genes from Simons Foundation Autism Research Initiative (SFARI) ; and (3) FMRP binding targets  (Additional file 8: Table S8). The genes down-regulated in murine Fmr1KO cultures were enriched for MGI behavioral/neurological phenotype genes, human SFARI-related genes, and FMRP targets. However, the overlap between these three lists and differentially expressed genes in embryonic brain tissue was not significant (Additional file 9: Table S9, Additional file 10: Table S10 and Additional file 11: Table S11). We also found a significant overlap between SFARI and FMRP (p value = 4.41E−7), SFARI and MGI behavioral (p value < 2.2E−16), FMRP and MGI behavioral (p value = 1.03E−8) gene lists.
Differentially expressed genes in Fmr1KO in cultured neurons are preferentially expressed in early human brain development, whereas genes in embryonic hippocampal and cortical tissues correspond more closely to later developmental stages
To identify the predominant gene clusters in mice that corresponded to developing human neocortex and hippocampus, we used PCA to compare 13,830 genes in sample space in the two datasets, respectively. Our reference model for development was the human neocortical and hippocampal transcriptome at 15 stages, profiled using Affymetrix Human Exon 1.0 ST Array as previously described . We integrated this reference data with our mouse data, which comprised 21,141 unique genes. Of the 16,492 unique minimal human Entrez IDs present on Affymetrix Exon ST 1.0 array, 14,653 genes correspond with homologous mouse Entrez Gene ID and 13,830 unique genes are shared between the Human Exon 1.0ST and Mouse 430 2.0 arrays (Additional file 12: Table S12).
Figure 2 shows the first two principal components for each gene (each dot signifies a gene). In the hippocampus, PC1 and PC2 captured 29 and 12 % of variation, respectively, and in the neocortex, they captured 27.7 and 11.3 % of variation, respectively. Visually, there appeared to be three large-scale patterns of co-expression in each of the developmental time series represented on PCA. Therefore, we used K-means clustering (k = 3) to clarify these clusters. The first cluster primarily had genes where the first PC was negative (PC1 < 0; H0 and C0 in the hippocampus and cortex respectively; magenta dots). Plotting these genes against developmental stage showed that they were typically up-regulated between the first and sixth stages of development and that their expression decreased subsequently (median value for each cluster based on RMA-normalized signal; right side of Fig. 2). Stages 1 to 7 are defined as fetal stages in Kang et al. . Genes with a positive first PC1 (PC1 > 0) were subdivided to two clusters, cluster 1 (green dots, H1/C1) and cluster 2 (blue dots, H2/C2). Genes whose expression was higher during developmental stages 7–11 are represented by H1 and C2 (infancy and childhood). H2 and C1 correspond primarily to genes up-regulated during young, middle, and late adulthood. The border between H0/H1 and C0/C2 may be defined as birth, with its associated broad changes in gene expression. Overall, we observed that during human hippocampal development, 42 % (5808), 31 % (4242), and 27 % (3780) of total genes resided in H0, H1, and H2 clusters, respectively. During cortical development, 40 % (5572), 30 % (4217), and 29 % (4041) of total genes resided in the C0, C1, and C2 clusters, respectively.
To investigate the expression profiles of murine homologues in the human data, we highlighted differentially expressed murine genes onto the human developmental transcriptome map. Figure 3 shows differentially expressed genes (blue: down-regulated, red: up-regulated, black: combined) for hippocampal and cortical KO/WT cultures and for embryonic hippocampus and cortex tissue.
The majority of genes that were up-regulated in murine Fmr1KO cortical and hippocampal cultures segregated into cluster 0. This finding is indicated by the increased density of red dots (genes) at the 9 o’clock periphery of the upper panels in Fig. 3 and the highest percentage in the corresponding cluster in Table 2. These findings indicate that genes up-regulated in KO mice tended to have increased expression in early development (than later) of the corresponding brain area.
Next, we highlighted differentially expressed genes in embryonic brain tissue in human development (Fig. 3c, d). We observed that the majority of differentially regulated genes in murine Fmr1KO primary hippocampus segregated to Cluster H1, which represents a later stage in human development (Table 2). The majority of up-regulated genes in murine Fmr1KO primary cortex segregated to cluster 2 (C2). This finding indicates that the majority of differentially expressed genes (mostly up-regulated genes which partition to H1 and C2 in the hippocampus and cortex) in primary brain tissue tended to be expressed at higher levels in later development. Overall, differentially expressed genes in murine culture were over-represented for early expression in humans (cortical culture: 95 % CI = 1.28–1.53, OR = 1.4; hippocampal culture: CI = 1.5–1.83, OR = 1.66), while the differentially expressed genes in tissue were under-represented at early stages (cortex: 95 % CI = 0.69–0.99, OR = 0.83; hippocampus: 95 % CI = 0.57–0.85, OR = 0.69) and over-represented in later stages (Additional file 13: Table S13).
Next, we determined if enrichment of up-regulated genes in KO mouse CNS in the early phases of development was a consequence of common non-cell-type-specific cell-cycle changes. We used the DAVID program to perform a KEGG pathway-level analysis of the 732 C0- and 821 H0-related up-regulated genes (red clusters in Fig. 3a, b). This analysis showed that in addition to cell-cycle-related processes, there was an over-representation of processes such as the neurotrophin signaling pathway and the T cell receptor signaling pathway in the C0-related subset of genes. In addition, the insulin signaling pathway, gap junction, endocytosis, ErbB signaling pathways and axon guidance were over-represented in the H0-related subset of genes (Additional file 14: Table S14).
The Fmr1KO mouse is a valuable model for studying FXS and a subset of patients with autism. We selected to work with FXS mouse model as it constitutes one of the more established and studied models of the more homogeneous monogenic forms of autism. Even so, we accept that it might not be fully representative for autism, which exhibits striking heterogeneity.
We investigated whether transcriptomic changes in Fmr1KO embryonic brain tissues were paralleled by those obtained in cultured Fmr1KO neurons. Although we found several commonalities, there were systematic differences that were surprising to us and, to our knowledge, are previously unreported. For example, in the in situ CNS of the Fmr1KO mouse, there was a significant enrichment of genes involved in immunological signaling; this pattern was not evident in cultured neurons. Conversely, in cultured neurons, we found a perturbation in several overlapping excitatory glutamatergic signaling cascades. Most are downstream of mGluR and neurotrophin signaling and may be coupled to FMRP-regulated mRNA translation.
As indicated by PCA, the distance between KO and WT clusters in cultures was higher than the distance between KO and WT in intact embryonic brain tissue. In addition, we observed a drastic difference in the numbers of differentially regulated genes in cultured cells and intact brain. Both findings suggest that the murine gene expression profile of cortical neuronal cultures and hippocampal neuronal cultures captures the difference between FXS and WT better than brain tissue. One possible explanation for this finding is that the effects of the Fmr1 mutation on gene expression are more directly observable in cultured cells because of the more simplistic nature of the model with fewer different cell types and absence of a buffer to perturbation (culture). Primary brain is a complex tissue that can develop compensatory mechanisms to disease and has many cell types, with different gene expression patterns, therefore diluting a particular phenotype. The phenotype is exaggerated in culture and subtle in the organism, with factors not present in vitro playing a role in disease pathogenesis and severity.
The findings in our cultured cells were congruent with the well-established classical mGluR theory of FXS [1, 12, 13]. Specifically, we found significant perturbation in LTP, neurotrophin signaling pathway, and mammalian target of rapamycin (mTOR) and mitogen-activated protein kinase (MAPK) signaling cascades in Fmr1KO cultures. Over-activation of these cascades could lead to abnormal synaptic function owing to exaggerated protein synthesis in FXS [2, 4]. Of the differentially expressed genes, several are key activators in these aforementioned pathways . Several have also been implicated in autism-related disorders , which have a high rate of comorbidity with FXS. However, as a large percentage (15–60 %) of children with FXS meet diagnostic criteria for autism, the findings might be relevant to a specific subset of cases with autism. Our findings support the hypothesis that the transcriptional variation in this Fragile X model is focused within mGluR-activated pathways related to synaptic plasticity.
Although the contrast between gene expression in brain tissue and in cultured neurons in Fmr1KO is novel, the finding of an immunological signal in the CNS is not. A module enriched for immune genes and glial markers has been observed in autism spectrum disorder (ASD) human post-mortem brain tissue [27, 28], and immunity response-related pathways were perturbed in Fmr1KO and Tsc2+/− murine cerebellum . Various studies have implicated a dysregulated immune system response in ASD [30–35], but the exact mechanisms by which this dysfunction relates to autism-related disorders are not well-established. Moreover, although the genes are labeled as immunological in function, they also function as morphogens in the developing CNS . Thus, their appearance in the intact CNS may mark a change in development rather than the inflammatory response they evince in the peripheral vasculature.
One possible explanation is that the end result of the FMR1-related disorders is an immunological phenotype that can be observed in the brain. However, culturing cells frees them from environmental factors that influence this phenotype, allowing them to manifest the synaptic-related dysfunctions of FXS. In addition, one of the reasons why synaptic and/or neuronal effects were more prominent in neuronal cultures may be because the immunological changes are less visible due to the lack of blood vessels or limited number and type of glia in cultured neurons. Another explanation is the presence of active glial cells in brain—especially microglia, which are the brain’s professional phagocytes. Transformation of microglia to reactive states in response to pathology has been known for decades as microglial activation . Neuroglial activation and neuroinflammation have been observed in the brains of autism patients [28, 35, 38–40] and reactive astrocytes in several brain regions of Fmr1KO mice were previously revealed . Moreover, the morphological phenotype observed in Fmr1KO neurons—the high density of dendritic spines [8–11]—might be linked to pruning defects by microglial cells [42–45]. The immunological phenotype therefore adds to prior evidence of a link between CNS pathology and glial activation [46, 47] in Fmr1KO mice .
We used Human Brain Transcriptome data to construct a map of genes that are active at different stages of human brain development. We defined three clusters of genes corresponding to prenatal, early, and late developmental time points in the cortex and hippocampus. We highlighted human genes homologous to the differentially expressed murine genes from culture and brain tissue, providing a temporal analysis of the expression of these genes. Up-regulated genes in KO cultures (as compared to WT cultures) tended to be genes whose physiological expression peaks early in healthy development. By contrast, differentially expressed genes in brain tissues were those whose expression peaks later in healthy development. Up-regulated genes from intact KO murine hippocampus and cortex were over-represented among genes with later physiological expression in human development. We hypothesize that the isolation and culture of neurons caused cortex and hippocampus cells to act “juvenile” or “de-differentiated” and therefore to express a “younger” transcriptomic profile. We have also found an enrichment of down-regulated genes in Fmr1KO cultures in literature databases on cognitive traits (MGI murine genes associated with behavioral/neurological phenotype, human SFARI genes, and FMRP binding targets). Together, we speculate that KO cultures bear a greater resemblance to earlier stages of development, and conversely, primary brain tissue bears a stronger resemblance to later stages of development.
The findings from this study are summarized in Table 3 and highlighted in the list below:
There is an over-representation of immune-related transcriptional activity in embryonic Fmr1KO cortex and hippocampus in comparison to neuronal signaling pathways in cortical and hippocampal primary cultures.
Importance: In addition to neuronal pathology, research on disease treatments should also consider the impact of immune dysregulation, which is comparatively less studied.
Up-regulated genes in murine KO cultures corresponded to genes whose peak expression occurs early in human development, while differentially expressed genes measured in murine KO embryonic brain tissue corresponded to genes whose peak expression occurs later in human development.
Importance: Cultured cells may recapitulate an early phase of the disease, which is also less obscured with “immunological” phenotype and in vivo compensatory mechanisms.
Down-regulated differentially expressed genes in KO cultures are enriched for (1) genes implicated in autism (SFARI), (2) genes associated with behavioral/neurological phenotypes of autism in mice (MGI), and (3) FMRP binding targets.
Importance: These findings suggest that there are potential shared disease-related mechanisms in fragile X and other autism phenotypes in mice.
Together, these results suggest that the Fmr1KO culture system:
Better captures the difference between genotypes in comparison to embryonic brain tissue.
Bears greater resemblance to earlier stages of human development
Bears the significant resemblance to the genes previously implicated with autism in human in comparison to brain tissue
To conclude, cultured hippocampal and cortical cells exhibited a more neurogenesis-related gene transcriptome, while primary brain tissue exhibited a more immunological transcriptome. These results are consistent with previous findings showing a role for the immunological response in neurodevelopmental disorders. Fmr1KO differentially expressed genes in culture coincided with genes active in early human sdevelopment, while the brain-related differentially expressed genes coincided with genes more active later in human development.
Availability of supporting data
The data set supporting the results of this article is available in the Gene Expression Omnibus repository with the accession identifier GSE71034 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE71034).
autism spectrum disorder
cell adhesion molecules
central nervous system
Database for Annotation, Visualization and Integrated Discovery
fragile X mental retardation 1
fragile X mental retardation protein
fragile X syndrome
human brain transcriptome
Kyoto Encyclopedia of Genes and Genomes
mitogen-activated protein kinase
Mouse Gene Informatics
metabotropic glutamate receptor
mammalian target of rapamycin
principal component analysis
robust multi-array average
Simons Foundation Autism Research Initiative
tuberous sclerosis protein 2
Bear MF, Huber KM, Warren ST. The mGluR theory of fragile X mental retardation. Trends Neurosci. 2004;27(7):370–7. doi:10.1016/j.tins.2004.04.009.
Bhakar AL, Dolen G, Bear MF. The pathophysiology of fragile X (and what it teaches us about synapses). Annu Rev Neurosci. 2012;35:417–43. doi:10.1146/annurev-neuro-060909-153138.
Koukoui SD, Chaudhuri A. Neuroanatomical, molecular genetic, and behavioral correlates of fragile X syndrome. Brain Res Rev. 2007;53(1):27–38. doi:10.1016/j.brainresrev.2006.06.001.
Huber KM, Gallagher SM, Warren ST, Bear MF. Altered synaptic plasticity in a mouse model of fragile X mental retardation. Proc Natl Acad Sci U S A. 2002;99(11):7746–50. doi:10.1073/pnas.122205699.
Dolen G, Osterweil E, Rao BS, Smith GB, Auerbach BD, Chattarji S, et al. Correction of fragile X syndrome in mice. Neuron. 2007;56(6):955–62. doi:10.1016/j.neuron.2007.12.001.
Dolen G, Carpenter RL, Ocain TD, Bear MF. Mechanism-based approaches to treating fragile X. Pharmacol Ther. 2010;127(1):78–93. doi:10.1016/j.pharmthera.2010.02.008.
Krueger DD, Bear MF. Toward fulfilling the promise of molecular medicine in fragile X syndrome. Annu Rev Med. 2011;62:411–29. doi:10.1146/annurev-med-061109-134644.
Antar LN, Li C, Zhang H, Carroll RC, Bassell GJ. Local functions for FMRP in axon growth cone motility and activity-dependent regulation of filopodia and spine synapses. Mol Cell Neurosci. 2006;32(1-2):37–48. doi:10.1016/j.mcn.2006.02.001.
Galvez R, Greenough WT. Sequence of abnormal dendritic spine development in primary somatosensory cortex of a mouse model of the fragile X mental retardation syndrome. Am J Med Genet A. 2005;135(2):155–60. doi:10.1002/ajmg.a.30709.
Irwin SA, Patel B, Idupulapati M, Harris JB, Crisostomo RA, Larsen BP, et al. Abnormal dendritic spine characteristics in the temporal and visual cortices of patients with fragile-X syndrome: a quantitative examination. Am J Med Genet. 2001;98(2):161–7.
Nimchinsky EA, Oberlander AM, Svoboda K. Abnormal development of dendritic spines in FMR1 knock-out mice. J Neurosci. 2001;21(14):5139–46.
Pfeiffer BE, Huber KM. Fragile X mental retardation protein induces synapse loss through acute postsynaptic translational regulation. J Neurosci. 2007;27(12):3120–30. doi:10.1523/JNEUROSCI.0054-07.2007.
Nakamoto M, Nalavadi V, Epstein MP, Narayanan U, Bassell GJ, Warren ST. Fragile X mental retardation protein deficiency leads to excessive mGluR5-dependent internalization of AMPA receptors. Proc Natl Acad Sci U S A. 2007;104(39):15537–42. doi:10.1073/pnas.0707484104.
Schmid PR, Palmer NP, Kohane IS, Berger B. Making sense out of massive data by going beyond differential expression. Proc Natl Acad Sci U S A. 2012;109(15):5594–9. doi:10.1073/pnas.1118792109.
Butte AJ, Kohane IS. Creation and implications of a phenome-genome network. Nat Biotechnol. 2006;24(1):55–62. doi:10.1038/nbt1150.
Kho AT, Zhao Q, Cai Z, Butte AJ, Kim JY, Pomeroy SL, et al. Conserved mechanisms across development and tumorigenesis revealed by a mouse development perspective of human cancers. Genes Dev. 2004;18(6):629–40. doi:10.1101/gad.1182504.
Naxerova K, Bult CJ, Peaston A, Fancher K, Knowles BB, Kasif S, et al. Analysis of gene expression in a developmental context emphasizes distinct biological leitmotifs in human cancers. Genome Biol. 2008;9(7):R108. doi:10.1186/gb-2008-9-7-r108.
Bhakar AL, Tannis LL, Zeindler C, Russo MP, Jobin C, Park DS, et al. Constitutive nuclear factor-kappa B activity is required for central neuron survival. J Neurosci. 2002;22(19):8466–75.
Melen E, Kho AT, Sharma S, Gaedigk R, Leeder JS, Mariani TJ, et al. Expression analysis of asthma candidate genes during human and murine lung development. Respir Res. 2011;12:86. doi:10.1186/1465-9921-12-86.
Alter O, Brown PO, Botstein D. Singular value decomposition for genome-wide expression data processing and modeling. Proc Natl Acad Sci U S A. 2000;97(18):10101–6.
Misra J, Schmitt W, Hwang D, Hsiao LL, Gullans S, Stephanopoulos G. Interactive exploration of microarray gene expression patterns in a reduced dimensional space. Genome Res. 2002;12(7):1112–20. doi:10.1101/gr.225302.
Mouse Gene Informatics. http://www.informatics.jax.org/searches/Phat.cgi?id=MP:0005386. Accessed May, 2014.
Simons Foundation Autism Research Initiative. https://gene.sfari.org. Accessed on May, 2015.
Ascano Jr M, Mukherjee N, Bandaru P, Miller JB, Nusbaum JD, Corcoran DL, et al. FMRP targets distinct mRNA sequence elements to regulate protein expression. Nature. 2012;492(7429):382–6. doi:10.1038/nature11737.
Kang HJ, Kawasawa YI, Cheng F, Zhu Y, Xu X, Li M, et al. Spatio-temporal transcriptome of the human brain. Nature. 2011;478(7370):483–9. doi:10.1038/nature10523.
Huguet G, Ey E, Bourgeron T. The genetic landscapes of autism spectrum disorders. Annu Rev Genomics Hum Genet. 2013;14:191–213. doi:10.1146/annurev-genom-091212-153431.
Voineagu I, Wang X, Johnston P, Lowe JK, Tian Y, Horvath S, et al. Transcriptomic analysis of autistic brain reveals convergent molecular pathology. Nature. 2011;474(7351):380–4. doi:10.1038/nature10110.
Gupta S, Ellis SE, Ashar FN, Moes A, Bader JS, Zhan J, et al. Transcriptome analysis reveals dysregulation of innate immune response genes and neuronal activity-dependent genes in autism. Nat Commun. 2014;5:5748. doi:10.1038/ncomms6748.
Kong SW, Sahin M, Collins CD, Wertz MH, Campbell MG, Leech JD, et al. Divergent dysregulation of gene expression in murine models of fragile X syndrome and tuberous sclerosis. Mol Autism. 2014;5(1):16. doi:10.1186/2040-2392-5-16.
Jung JY, Kohane IS, Wall DP. Identification of autoimmune gene signatures in autism. Transl Psychiatry. 2011;1:e63. doi:10.1038/tp.2011.62.
Kong SW, Collins CD, Shimizu-Motohashi Y, Holm IA, Campbell MG, Lee IH, et al. Characteristics and predictive value of blood transcriptome signature in males with autism spectrum disorders. PLoS One. 2012;7(12):e49475. doi:10.1371/journal.pone.0049475.
Garbett K, Ebert PJ, Mitchell A, Lintas C, Manzi B, Mirnics K, et al. Immune transcriptome alterations in the temporal cortex of subjects with autism. Neurobiol Dis. 2008;30(3):303–11. doi:10.1016/j.nbd.2008.01.012.
Korvatska E, Van de Water J, Anders TF, Gershwin ME. Genetic and immunologic considerations in autism. Neurobiol Dis. 2002;9(2):107–25. doi:10.1006/nbdi.2002.0479.
Cohly HH, Panja A. Immunological findings in autism. Int Rev Neurobiol. 2005;71:317–41.
Li X, Chauhan A, Sheikh AM, Patil S, Chauhan V, Li XM, et al. Elevated immune response in the brain of autistic patients. J Neuroimmunol. 2009;207(1-2):111–6. doi:10.1016/j.jneuroim.2008.12.002.
Garay PA, McAllister AK. Novel roles for immune molecules in neural development: implications for neurodevelopmental disorders. Front Synaptic Neurosci. 2010;2:136. doi:10.3389/fnsyn.2010.00136.
Hanisch UK, Kettenmann H. Microglia: active sensor and versatile effector cells in the normal and pathologic brain. Nat Neurosci. 2007;10(11):1387–94. doi:10.1038/nn1997.
Vargas DL, Nascimbene C, Krishnan C, Zimmerman AW, Pardo CA. Neuroglial activation and neuroinflammation in the brain of patients with autism. Ann Neurol. 2005;57(1):67–81. doi:10.1002/ana.20315.
Suzuki K, Sugihara G, Ouchi Y, Nakamura K, Futatsubashi M, Takebayashi K, et al. Microglial activation in young adults with autism spectrum disorder. JAMA Psychiatry. 2013;70(1):49–58. doi:10.1001/jamapsychiatry.2013.272.
Morgan JT, Chana G, Pardo CA, Achim C, Semendeferi K, Buckwalter J, et al. Microglial activation and increased microglial density observed in the dorsolateral prefrontal cortex in autism. Biol Psychiatry. 2010;68(4):368–76. doi:10.1016/j.biopsych.2010.05.024.
Yuskaitis CJ, Beurel E, Jope RS. Evidence of reactive astrocytes but not peripheral immune system activation in a mouse model of Fragile X syndrome. Biochim Biophys Acta. 2010;1802(11):1006–12. doi:10.1016/j.bbadis.2010.06.015.
Kim KH, Son SM, Mook-Jung I. Contributions of microglia to structural synaptic plasticity. J Exp Neurosci. 2013;7:85–91. doi:10.4137/JEN.S11269.
Paolicelli RC, Bolasco G, Pagani F, Maggi L, Scianni M, Panzanelli P, et al. Synaptic pruning by microglia is necessary for normal brain development. Science. 2011;333(6048):1456–8. doi:10.1126/science.1202529.
Schafer DP, Lehrman EK, Kautzman AG, Koyama R, Mardinly AR, Yamasaki R, et al. Microglia sculpt postnatal neural circuits in an activity and complement-dependent manner. Neuron. 2012;74(4):691–705. doi:10.1016/j.neuron.2012.03.026.
Schafer DP, Stevens B. Synapse elimination during development and disease: immune molecules take centre stage. Biochem Soc Trans. 2010;38(2):476–81. doi:10.1042/BST0380476.
Ashwood P, Wills S, Van de Water J. The immune response in autism: a new frontier for autism research. J Leukoc Biol. 2006;80(1):1–15. doi:10.1189/jlb.1205707.
Pardo CA, Vargas DL, Zimmerman AW. Immunity, neuroglia and neuroinflammation in autism. Int Rev Psychiatry. 2005;17(6):485–95. doi:10.1080/02646830500381930.
We thank the members of the laboratories of Dr. Isaac Kohane and Dr. Mark Bear for the helpful discussions and valuable advice. ISK, DP, and SWK acknowledge the support of the Conte Center for Computational Neuropsychiatric Genomics (NIH P50MH94267). ISK and NPP acknowledge the support of Neuropsychiatric Genome-Scale and RDOC Individualized Domains (NIMH/NHGRI P50MH106933). DP acknowledges the support of NIH/NHLBI 5U54HL117798. ATK acknowledges the support of NHLBI K25 HL091124. NPP acknowledges the support of NIH/NHGRI 1U54HG007963 and Aetna Foundation (A20357). SWK acknowledges the support of NIH/NIMH R01MH107205.
MFB holds patents on the use of mGluR5 inhibitors for treatment of fragile X and autism. DMM is a Chairman of the Board of Q-State Biosciences, a private company working in the space of disease modeling. The other authors declare no competing financial interests.
DP contributed to study design, carried out cellular and molecular studies, including RNA extraction and processing of microarrays, analyzed and interpreted the data, and wrote the manuscript. ATK carried out statistical analysis, contributed to data interpretation, and participated in the preparation of the manuscript. NPP contributed to study design, analyzed the data, and critically revised the manuscript. ALB contributed to study design and coordination, provided cultures and tissues from Fmr1 knockout and WT mice, contributed to cellular and molecular studies, including RNA extraction, and critically revised the manuscript. NSM contributed to study design, carried out cellular studies, and participated in the preparation of the manuscript. SWK contributed to data analysis and interpretation and critically revised the manuscript. DMM conceived the study, participated in its design and coordination, and critically revised the manuscript. MFB contributed to study concept and study design, contributed to data interpretation, and critically revised the manuscript. ISK conceived and supervised the study, participated in its design and coordination, and critically revised the manuscript. All authors have read and approved the final manuscript.
Genes associated with defects in cognitive development: MGI, SFARI, and FMRP targets. (XLSX 111 kb)
Centroid calculation for each cluster in Fig. 1. The distance between clusters is represented as a difference of absolute values on the PC1 and PC2 axes. (XLSX 9 kb)
List of significant differential gene expressed in four systems. (XLSX 1 mb)
Detailed KEGG pathway-enrichment analysis of differentially expressed genes between KO and WT mice. (XLSX 21 kb)
Selected over-represented KEGG pathways with core differentially expressed genes in Fmr1KO that made a pathway significant. (XLSX 11 kb)
Metacore enrichment analysis of differentially expressed genes in culture and in primary brain based on differentially affected functional ontologies such as (A) “Process Networks” and (B) “Diseases (by Biomarkers)”. (DOCX 124 kb)
Overlap of pathways and commonalities between two different types of preparations (brain and culture) and between two different brain regions (cortex and hippocampus). (XLSX 237 kb)
Unique and annotated mouse genes after differential gene expression analysis in four systems. Mapping to MGI genes, human SFARI genes, and FMRP targets. (XLSX 4 mb)
Significance of overlap for hippocampal or cortical cultures and embryonic brain tissues of KO/WT up/down/combine differentially expressed genes AND MGI genes with behavioral/neurological phenotype (http://www.informatics.jax.org/searches/Phat.cgi?id=MP:0005386). Fisher’s exact test was used to calculate the statistics. (XLSX 11 kb)
Significance of overlap for hippocampal or cortical cultures and embryonic brain tissues of KO/WT up/down/combine differentially expressed genes AND human SFARI genes. Fisher’s exact test was used to calculate the statistics. (XLSX 12 kb)
Significance of overlap for hippocampal or cortical cultures and embryonic brain tissues of KO/WT up/down/combine differentially expressed genes AND FMRP binding targets. Fisher’s exact test was used to calculate the statistics. (XLSX 12 kb)
Integration of mouse and human data to represent development. Final number of genes in common between the datasets: 13830. (XLSX 2 mb)
Significance of overlap for murine differentially expressed genes (up + down) and genes represented in a specific cluster in human development. Fisher’s exact test was used to calculate the statistics. (XLSX 16 kb)
Enriched pathways in murine differentially expressed up-regulated genes in cortical and hippocampal cultures highlighted on human development (C0 and H0 clusters respectively). (XLSX 11 kb)
About this article
Cite this article
Prilutsky, D., Kho, A.T., Palmer, N.P. et al. Gene expression analysis in Fmr1KO mice identifies an immunological signature in brain tissue and mGluR5-related signaling in primary neuronal cultures. Molecular Autism 6, 66 (2015). https://doi.org/10.1186/s13229-015-0061-9