Gene expression analysis in Fmr1KO mice identifies an immunological signature in brain tissue and mGluR5-related signaling in primary neuronal cultures
© Prilutsky et al. 2015
Received: 20 August 2015
Accepted: 10 December 2015
Published: 21 December 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.
KeywordsFragile X syndrome Murine model Gene expression Neuronal cultures Brain
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.
- (1)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.
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
Enriched pathways in differentially expressed genes in brain and primary cultures of Fmr1KO mice
KEGG pathway (ENTREZ ID as input)
No. of genes
mmu04910:Insulin signaling pathway
mmu00900:Terpenoid backbone biosynthesis
mmu04120:Ubiquitin mediated proteolysis
mmu05211:Renal cell carcinoma
mmu04722:Neurotrophin signaling pathway
mmu04012:ErbB signaling pathway
mmu04914:Progesterone-mediated oocyte maturation
mmu04150:mTOR signaling pathway
mmu05220:Chronic myeloid leukemia
mmu04622:RIG-I-like receptor signaling pathway
mmu04722:Neurotrophin signaling pathway
mmu04660:T cell receptor signaling pathway
mmu04010:MAPK signaling pathway
mmu04120:Ubiquitin mediated proteolysis
mmu05220:Chronic myeloid leukemia
mmu00534:Heparan sulfate biosynthesis
mmu04062:Chemokine signaling pathway
Hippocampus in situ
mmu04060:Cytokine-cytokine receptor interaction
mmu05322:Systemic lupus erythematosus
mmu04514:Cell adhesion molecules (CAMs)
mmu05320:Autoimmune thyroid disease
mmu04672:Intestinal immune network for IgA production
mmu04610:Complement and coagulation cascades
Cortex in situ
mmu00590:Arachidonic acid metabolism
mmu04640:Hematopoietic cell lineage
mmu04060:Cytokine-cytokine receptor interaction
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).
Integration of human neocortex/hippocampus developmental genes and mouse differentially expressed genes in culture and brain
Differentially expressed genes in cultured cells
Hippo culture KO/WT down
Cortex culture KO/WT down
Hippo culture KO/WT up
Cortex culture KO/WT up
Hippo culture KO/WT combined
Cortex culture KO/WT combined
Differentially expressed genes in primary brain tissue
Hippo brain KO/WT down
Cortex brain KO/WT down
Hippo brain KO/WT up
Cortex brain KO/WT up
Hippo brain KO/WT combined
Cortex brain KO/WT combined
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.
Summary of findings
Differentially expressed genes in Fmr1KO vs WT (disease genes)
Number of differentially expressed genes
Distance between KO and WT clusters on PCA (PC1/PC2)
Early (H0 cluster)
Early (C0 cluster)
Later (H1 cluster)
Later (C2 cluster)
Down- and up- regulated
SFARI human genes
Significantly enriched in down-regulated genes
Significantly enriched in down-regulated genes
MGI mouse genes associated with behavioral/neurological phenotype
Significantly enriched in down-regulated genes
Significantly enriched in down-regulated genes
Significantly enriched in down-regulated genes
Significantly enriched in down- and up- regulated genes
- 1.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.
- 2.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.
- 3.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.
- 4.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
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.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Nimchinsky EA, Oberlander AM, Svoboda K. Abnormal development of dendritic spines in FMR1 knock-out mice. J Neurosci. 2001;21(14):5139–46.PubMedGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Butte AJ, Kohane IS. Creation and implications of a phenome-genome network. Nat Biotechnol. 2006;24(1):55–62. doi:10.1038/nbt1150.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Jung JY, Kohane IS, Wall DP. Identification of autoimmune gene signatures in autism. Transl Psychiatry. 2011;1:e63. doi:10.1038/tp.2011.62.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Cohly HH, Panja A. Immunological findings in autism. Int Rev Neurobiol. 2005;71:317–41.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Pardo CA, Vargas DL, Zimmerman AW. Immunity, neuroglia and neuroinflammation in autism. Int Rev Psychiatry. 2005;17(6):485–95. doi:10.1080/02646830500381930.PubMedView ArticleGoogle Scholar