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.     

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.     

DANCER: a program for digital anatomical reconstruction of gene expression data

A digital anatomy construction (DANCER) program was developed for gene expression data. DANCER can be used to reconstruct anatomical images from in situ hybridization images, microarray or other gene expression data. The program fills regions of a drawn figure with the corresponding values from a gene expression data set. The output of the program presents the expression levels of a particular gene in a particular region relative to other regions. The program was tested with values from experimental in situ hybridization autoradiographs and from a microarray experiment. Reconstruction of in situ hybridization data from adult rat brain made by DANCER corresponded well with the original autoradiograph. Reconstruction of microarray data from adult mouse brains provided images that reflect actual expression levels. This program should help to provide visualization and interpretation of data derived from gene expression experiments. DANCER may be freely downloaded.     

PEATmoss (Physcomitrella Expression Atlas Tool): a unified gene expression atlas for the model plant Physcomitrella patens

Physcomitrella patens is a bryophyte model plant that is often used to study plant evolution and development. Its resources are of great importance for comparative genomics and evo‐devo approaches. However, expression data from Physcomitrella patens were so far generated using different gene annotation versions and three different platforms: CombiMatrix and NimbleGen expression microarrays and RNA sequencing. The currently available P. patens expression data are distributed across three tools with different visualization methods to access the data. Here, we introduce an interactive expression atlas, Physcomitrella Expression Atlas Tool (PEATmoss), that unifies publicly available expression data for P. patens and provides multiple visualization methods to query the data in a single web‐based tool. Moreover, PEATmoss includes 35 expression experiments not previously available in any other expression atlas. To facilitate gene expression queries across different gene annotation versions, and to access P. patens annotations and related resources, a lookup database and web tool linked to PEATmoss was implemented. PEATmoss can be accessed at https://peatmoss.online.uni-marburg.de     

A Digital Framework to Build,Visualize and Analyze a Gene Expression Atlas with Cellular Resolution in Zebrafish Early Embryogenesis

A gene expression atlas is an essential resource to quantify and understand the multiscale processes of embryogenesis in time and space. The automated reconstruction of a prototypic 4D atlas for vertebrate early embryos, using multicolor fluorescence in situ hybridization with nuclear counterstain, requires dedicated computational strategies. To this goal, we designed an original methodological framework implemented in a software tool called Match-IT. With only minimal human supervision, our system is able to gather gene expression patterns observed in different analyzed embryos with phenotypic variability and map them onto a series of common 3D templates over time, creating a 4D atlas. This framework was used to construct an atlas composed of 6 gene expression templates from a cohort of zebrafish early embryos spanning 6 developmental stages from 4 to 6.3 hpf (hours post fertilization). They included 53 specimens, 181,415 detected cell nuclei and the segmentation of 98 gene expression patterns observed in 3D for 9 different genes. In addition, an interactive visualization software, Atlas-IT, was developed to inspect, supervise and analyze the atlas. Match-IT and Atlas-IT, including user manuals, representative datasets and video tutorials, are publicly and freely available online. We also propose computational methods and tools for the quantitative assessment of the gene expression templates at the cellular scale, with the identification, visualization and analysis of coexpression patterns, synexpression groups and their dynamics through developmental stages.     

Wholemount immunohistochemistry for revealing complex brain topography

The repeated and well-understood cellular architecture of the cerebellum make it an ideal model system for exploring brain topography. Underlying its relatively uniform cytoarchitecture is a complex array of parasagittal domains of gene and protein expression. The molecular compartmentalization of the cerebellum is mirrored by the anatomical and functional organization of afferent fibers. To fully appreciate the complexity of cerebellar organization we previously refined a wholemount staining approach for high throughput analysis of patterning defects in the mouse cerebellum. This protocol describes in detail the reagents, tools, and practical steps that are useful to successfully reveal protein expression patterns in the adult mouse cerebellum by using wholemount immunostaining. The steps highlighted here demonstrate the utility of this method using the expression of zebrinII/aldolaseC as an example of how the fine topography of the brain can be revealed in its native three-dimensional conformation. Also described are adaptations to the protocol that allow for the visualization of protein expression in afferent projections and large cerebella for comparative studies of molecular topography. To illustrate these applications, data from afferent staining of the rat cerebellum are included.     

Apolipoprotein D in the Niemann-Pick type C disease mouse brain: An ultrastructural immunocytochemical analysis

Enhanced apoD gene expression and abnormally high levels of apoD protein accumulation in the brain have been previously documented as features of the neurodegenerative disorder, Niemann-Pick Type C disease (NP-C). In the present study we have used immunocytochemistry and light and electron microscopy to elucidate the cellular and subcellular distribution of apoD in the Balb/c NIH npc1 ^?/? mouse brain. The normal mouse brain demonstrates low levels of apoD-expressing glia particularly in the cerebellar white matter. In contrast, abundant, strongly apoD-immunolabeled cells were observed in select grey matter nuclei, including the globus pallidus, thalamus, and substantia nigra, and in white matter tracts within the internal capsule and cerebellum of NP-C mouse brain. These brains regions have been previously shown to display the most significant neurodegenerative changes in the NP-C mouse. Ultrastructural analysis revealed dense apoD immunoreactivity on the nuclear envelopes of cells that have the morphological features of oligodendrocyte precursor-like cells and light staining on astrocytes. These results define the cellular and subcellular pattern of expression of apoD in NP-C mouse brain and suggest a possible role for this lipocalin in the pathophysiology of this disorder.     

EpoDB: a prototype database for the analysis of genes expressed during vertebrate erythropoiesis.

EpoDB is a database of genes expressed in vertebrate red blood cells. It is also a prototype for the creation of cell and tissue-specific databases from multiple external sources. The information in EpoDB obtained from GenBank, SWISS-PROT, Transfac, TRRD and GERD is curated to provide high quality data for sequence analysis aimed at understanding gene regulation during erythropoiesis. New protocols have been developed for data integration and updating entries. Using a BLAST-based algorithm, we have grouped GenBank entries representing the same gene together. This sequence similarity protocol was also used to identify new entries to be included in EpoDB. We have recently implemented our database in Sybase (relational tables) in addition to SICStus Prolog to provide us with greater flexibility in asking complex queries that utilize information from multiple sources. New additions to the public web site (http://www.cbil.upenn.edu/epodb) for accessing EpoDB are the ability to retrieve groups of entries representing different variants of the same gene and to retrieve gene expression data. The BLAST query has been enhanced by incorporating BLASTView, an interactive and graphical display of BLAST results. We have also enhanced the queries for retrieving sequence from specified genes by the addition of MEME, a motif discovery tool, to the integrated analysis tools which include CLUSTALW and TESS.     

AVA: visual analysis of gene expression microarray data

SUMMARY: AVA (Array Visual Analyzer) is a Java program that provides a graphical environment for visualization and analysis of gene expression microarray data. Together with its interactive visualization tools and a variety of built-in data analysis and filtration methods, AVA effectively integrates microarray data normalization, quality assessment, and data mining into one application. AVAILABILITY: The software is freely available for academic users on request from the authors.