If I only had a brain: exploring mouse brain images in the Allen Brain Atlas

The combination of a powerful well-designed user interface with detailed high-quality data sets can create new possibilities for data exploration and analysis. The Allen Brain Atlas (http://www.brain-map.org) provides a collection of tools for examining a set of images that detail gene expression in the mouse brain. Powerful web-based viewers for individual images and parallel examination of related images interact with an external application for three-dimensional views. The underlying dataset, generated via high-throughput analysis of expression patterns of more than 21,000 genes in adult mouse brains, provides three-dimensional views of gene expression patterns displayed in the context of an anatomical ontology. Facilities for filtering views, saving views of interest, annotating images and sharing views via email support the ongoing process of analysis and provide a model for the future of integrated tools for analysing large image data sets.     

Surface-based mapping of gene expression and probabilistic expression maps in the mouse cortex

The Allen Brain Atlas (ABA, www.brain-map.org) is a genome wide, spatially registered collection of cellular resolution in situ hybridization gene expression image data of the C57Bl/6J mouse brain. Derived from the ABA, the Anatomic Gene Expression Atlas (AGEA, http://mouse.brain-map.org/agea) has demonstrated both laminar and areal spatial gene expression correlations in the mouse cortex. While the mouse cortex is lissencephalic, its curvature and substantial bending in boundary areas renders it difficult to visualize and analyze laminar versus areal effects in a rectilinear coordinate framework. In context of human and non-human primate cortex, surface-based representation has proven useful for understanding relative locations of laminar, columnar, and areal features. In this paper, we describe a methodology for constructing surface-based flatmaps of the mouse cortex that enables mapping of gene expression data from individual genes in the ABA, or probabilistic expression maps from the AGEA, to identify and visualize genetic relationships between layers and areas.     

Co-expression Profiling of Autism Genes in the Mouse Brain

Autism spectrum disorder (ASD) is one of the most prevalent and highly heritable neurodevelopmental disorders in humans. There is significant evidence that the onset and severity of ASD is governed in part by complex genetic mechanisms affecting the normal development of the brain. To date, a number of genes have been associated with ASD. However, the temporal and spatial co-expression of these genes in the brain remain unclear. To address this issue, we examined the co-expression network of 26 autism genes from AutDB (http://mindspec.org/autdb.html), in the framework of 3,041 genes whose expression energies have the highest correlation between the coronal and sagittal images from the Allen Mouse Brain Atlas database (http://mouse.brain-map.org). These data were derived from in situ hybridization experiments conducted on male, 56-day old C57BL/6J mice co-registered to the Allen Reference Atlas, and were used to generate a normalized co-expression matrix indicating the cosine similarity between expression vectors of genes in this database. The network formed by the autism-associated genes showed a higher degree of co-expression connectivity than seen for the other genes in this dataset (Kolmogorov–Smirnov P = 5×10⁻²⁸). Using Monte Carlo simulations, we identified two cliques of co-expressed genes that were significantly enriched with autism genes (A Bonferroni corrected P<0.05). Genes in both these cliques were significantly over-expressed in the cerebellar cortex (P = 1×10⁻⁵) suggesting possible implication of this brain region in autism. In conclusion, our study provides a detailed profiling of co-expression patterns of autism genes in the mouse brain, and suggests specific brain regions and new candidate genes that could be involved in autism etiology.     

Study of gene function based on spatial co-expression in a high-resolution mouse brain atlas

Background  

The Allen Brain Atlas (ABA) project systematically profiles three-dimensional high-resolution gene expression in postnatal mouse brains for thousands of genes. By unveiling gene behaviors at both the cellular and molecular levels, ABA is becoming a unique and comprehensive neuroscience data source for decoding enigmatic biological processes in the brain. Given the unprecedented volume and complexity of the in situ hybridization image data, data mining in this area is extremely challenging. Currently, the ABA database mainly serves as an online reference for visual inspection of individual genes; the underlying rich information of this large data set is yet to be explored by novel computational tools. In this proof-of-concept study, we studied the hypothesis that genes sharing similar three-dimensional expression profiles in the mouse brain are likely to share similar biological functions.     

Relationships between gene expression and brain wiring in the adult rodent brain

We studied the global relationship between gene expression and neuroanatomical connectivity in the adult rodent brain. We utilized a large data set of the rat brain "connectome" from the Brain Architecture Management System (942 brain regions and over 5000 connections) and used statistical approaches to relate the data to the gene expression signatures of 17,530 genes in 142 anatomical regions from the Allen Brain Atlas. Our analysis shows that adult gene expression signatures have a statistically significant relationship to connectivity. In particular, brain regions that have similar expression profiles tend to have similar connectivity profiles, and this effect is not entirely attributable to spatial correlations. In addition, brain regions which are connected have more similar expression patterns. Using a simple optimization approach, we identified a set of genes most correlated with neuroanatomical connectivity, and find that this set is enriched for genes involved in neuronal development and axon guidance. A number of the genes have been implicated in neurodevelopmental disorders such as autistic spectrum disorder. Our results have the potential to shed light on the role of gene expression patterns in influencing neuronal activity and connectivity, with potential applications to our understanding of brain disorders. Supplementary data are available at http://www.chibi.ubc.ca/ABAMS.     

Clustering of spatial gene expression patterns in the mouse brain and comparison with classical neuroanatomy

Spatial gene expression profiles provide a novel means of exploring the structural organization of the brain. Computational analysis of these patterns is made possible by genome-scale mapping of the C57BL/6J mouse brain in the Allen Brain Atlas. Here we describe methodology used to explore the spatial structure of gene expression patterns across a set of 3041 genes chosen on the basis of consistency across experimental observations (N = 2). The analysis was performed on smoothed, co-registered 3D expression volumes for each gene obtained by aggregating cellular resolution image data. Following dimensionality and noise reduction, voxels were clustered according to similarity of expression across the gene set. We illustrate the resulting parcellations of the mouse brain for different numbers of clusters (K) and quantitatively compare these parcellations with a classically-defined anatomical reference atlas at different levels of granularity, revealing a high degree of correspondence. These observations suggest that spatial localization of gene expression offers substantial promise in connecting knowledge at the molecular level with higher-level information about brain organization.     

Comparative analysis of selenocysteine machinery and selenoproteome gene expression in mouse brain identifies neurons as key functional sites of selenium in mammals

Although dietary selenium (Se) deficiency results in phenotypes associated with selenoprotein depletion in various organs, the brain is protected from Se loss. To address the basis for the critical role of Se in brain function, we carried out comparative gene expression analyses for the complete selenoproteome and associated biosynthetic factors. Using the Allen Brain Atlas, we evaluated 159 regions of adult mouse brain and provided experimental analyses of selected selenoproteins. All 24 selenoprotein mRNAs were expressed in the mouse brain. Most strikingly, neurons in olfactory bulb, hippocampus, cerebral cortex, and cerebellar cortex were exceptionally rich in selenoprotein gene expression, in particular in GPx4, SelK, SelM, SelW, and Sep15. Over half of the selenoprotein genes were also expressed in the choroid plexus. A unique expression pattern was observed for one of the highly expressed selenoprotein genes, SelP, which we suggest to provide neurons with Se. Cluster analysis of the expression data linked certain selenoproteins and selenocysteine machinery genes and suggested functional linkages among selenoproteins, such as that between SelM and Sep15. Overall, this study suggests that the main functions of selenium in mammals are confined to certain neurons in the brain.     

Harnessing Gene Expression Networks to Prioritize Candidate Epileptic Encephalopathy Genes

We apply a novel gene expression network analysis to a cohort of 182 recently reported candidate Epileptic Encephalopathy genes to identify those most likely to be true Epileptic Encephalopathy genes. These candidate genes were identified as having single variants of likely pathogenic significance discovered in a large-scale massively parallel sequencing study. Candidate Epileptic Encephalopathy genes were prioritized according to their co-expression with 29 known Epileptic Encephalopathy genes. We utilized developing brain and adult brain gene expression data from the Allen Human Brain Atlas (AHBA) and compared this to data from Celsius: a large, heterogeneous gene expression data warehouse. We show replicable prioritization results using these three independent gene expression resources, two of which are brain-specific, with small sample size, and the third derived from a heterogeneous collection of tissues with large sample size. Of the nineteen genes that we predicted with the highest likelihood to be true Epileptic Encephalopathy genes, two (GNAO1 and GRIN2B) have recently been independently reported and confirmed. We compare our results to those produced by an established in silico prioritization approach called Endeavour, and finally present gene expression networks for the known and candidate Epileptic Encephalopathy genes. This highlights sub-networks of gene expression, particularly in the network derived from the adult AHBA gene expression dataset. These networks give clues to the likely biological interactions between Epileptic Encephalopathy genes, potentially highlighting underlying mechanisms and avenues for therapeutic targets.     

The genomic determinants of alcohol preference in mice

Searches for the identity of genes that influence the levels of alcohol consumption by humans and other animals have often been driven by presupposition of the importance of particular gene products in determining positively or negatively reinforcing effects of ethanol. We have taken an unbiased approach and performed a meta-analysis across three types of mouse populations to correlate brain gene expression with levels of alcohol intake. Our studies, using filtering procedures based on QTL analysis, produced a list of eight candidate genes with highly heritable expression, which could explain a significant amount of the variance in alcohol preference in mice. Using the Allen Brain Atlas for gene expression, we noted that the candidate genes' expression was localized to the olfactory and limbic areas as well as to the orbitofrontal cortex. Informatics techniques and pathway analysis illustrated the role of the candidate genes in neuronal migration, differentiation, and synaptic remodeling. The importance of olfactory cues, learning and memory formation (Pavlovian conditioning), and cortical executive function, for regulating alcohol intake by animals (including humans), is discussed.     

Quantitative methods for genome-scale analysis of in situ hybridization and correlation with microarray data

With the emergence of genome-wide colorimetric in situ hybridization (ISH) data sets such as the Allen Brain Atlas, it is important to understand the relationship between this gene expression modality and those derived from more quantitative based technologies. This study introduces a novel method for standardized relative quantification of colorimetric ISH signal that enables a large-scale cross-platform expression level comparison of ISH with two publicly available microarray brain data sources.     

Modeling cholesterol metabolism by gene expression profiling in the hippocampus

An important part of the challenge of building models of biochemical reactions is determining reaction rate constants that transform substrates into products. We present a method to derive enzymatic kinetic values from mRNA expression levels for modeling biological networks without requiring further tuning. The core metabolic reactions of cholesterol in the brain, particularly in the hippocampus, were simulated. To build the model the baseline mRNA expression levels of genes involved in cholesterol metabolism were obtained from the Allen Mouse Brain Atlas. The model is capable of replicating the trends of relative cholesterol levels in Alzheimer's and Huntington's diseases; and reliably simulated SLOS, desmosterolosis, and Dhcr14/Lbr knockout studies. A sensitivity analysis correctly uncovers the Hmgcr, Idi2 and Fdft1 sites that regulate cholesterol homeostasis. Overall, our model and methodology can be used to pinpoint key reactions, which, upon manipulation, may predict altered cholesterol levels and reveal insights into potential drug therapy targets under diseased conditions.     

Validating Serum S100B and Neuron-Specific Enolase as Biomarkers for the Human Brain - A Combined Serum, Gene Expression and MRI Study

Introduction

Former studies have investigated the potential of serum biomarkers for diseases affecting the human brain. In particular the glial protein S100B, a neuro- and gliotrophin inducing plasticity, seems to be involved in the pathogenesis and treatment of psychiatric diseases such as major depression and schizophrenia. Neuron-specific enolase (NSE) is a specific serum marker for neuronal damage. However, the specificity of these biomarkers for cell type and brain region has not been investigated in vivo until now.

Methods

We acquired two magnetic resonance imaging parameters sensitive to changes in gray and white matter (T₁-weighted/diffusion tensor imaging) and obtained serum S100B and NSE levels of 41 healthy subjects. Additionally, we analyzed whole brain gene expressions of S100B in another male cohort of three subjects using the Allen Brain Atlas. Furthermore, a female post mortal brain was investigated using double immunofluorescence labelling with oligodendrocyte markers.

Results

We show that S100B is specifically related to white matter structures, namely the corpus callosum, anterior forceps and superior longitudinal fasciculus in female subjects. This effect was observed in fractional anisotropy and radial diffusivity – the latest an indicator of myelin changes. Histological data confirmed a co-localization of S100B with oligodendrocyte markers in the human corpus callosum. S100B was most abundantly expressed in the corpus callosum according to the whole genome Allen Human Brain Atlas. In addition, NSE was related to gray matter structures, namely the amygdala. This effect was detected across sexes.

Conclusion

Our data demonstrates a very high S100B expression in white matter tracts, in particular in human corpus callosum. Our study is the first in vivo study validating the specificity of the glial marker S100B for the human brain, and supporting the assumption that radial diffusivity represents a myelin marker. Our results open a new perspective for future studies investigating major neuropsychiatric disorders.     

Computational neuroanatomy and co-expression of genes in the adult mouse brain,analysis tools for the Allen Brain Atlas

We review quantitative methods and software developed to analyze genome-scale, brain-wide spatially-mapped gene-expression data. We expose new methods based on the underlying high-dimensional geometry of voxel space and gene space, and on simulations of the distribution of co-expression networks of a given size. We apply them to the Allen Atlas of the adult mouse brain, and to the co-expression network of a set of genes related to nicotine addiction retrieved from the NicSNP database. The computational methods are implemented in BrainGeneExpressionAnalysis (BGEA), a Matlab toolbox available for download.     

Laminar and Dorsoventral Molecular Organization of the Medial Entorhinal Cortex Revealed by Large-scale Anatomical Analysis of Gene Expression

Neural circuits in the medial entorhinal cortex (MEC) encode an animal’s position and orientation in space. Within the MEC spatial representations, including grid and directional firing fields, have a laminar and dorsoventral organization that corresponds to a similar topography of neuronal connectivity and cellular properties. Yet, in part due to the challenges of integrating anatomical data at the resolution of cortical layers and borders, we know little about the molecular components underlying this organization. To address this we develop a new computational pipeline for high-throughput analysis and comparison of in situ hybridization (ISH) images at laminar resolution. We apply this pipeline to ISH data for over 16,000 genes in the Allen Brain Atlas and validate our analysis with RNA sequencing of MEC tissue from adult mice. We find that differential gene expression delineates the borders of the MEC with neighboring brain structures and reveals its laminar and dorsoventral organization. We propose a new molecular basis for distinguishing the deep layers of the MEC and show that their similarity to corresponding layers of neocortex is greater than that of superficial layers. Our analysis identifies ion channel-, cell adhesion- and synapse-related genes as candidates for functional differentiation of MEC layers and for encoding of spatial information at different scales along the dorsoventral axis of the MEC. We also reveal laminar organization of genes related to disease pathology and suggest that a high metabolic demand predisposes layer II to neurodegenerative pathology. In principle, our computational pipeline can be applied to high-throughput analysis of many forms of neuroanatomical data. Our results support the hypothesis that differences in gene expression contribute to functional specialization of superficial layers of the MEC and dorsoventral organization of the scale of spatial representations.