Neuroinformatics for genome-wide 3D gene expression mapping in the mouse brain

Unsupervised pattern identification in spatial gene expression atlas reveals mouse brain regions beyond established ontology

The rapid growth of large-scale spatial gene expression data demands efficient and reliable computational tools to extract major trends of gene expression in their native spatial context. Here, we used stability-driven unsupervised learning (i.e., staNMF) to identify principal patterns (PPs) of 3D gene expression profiles and understand spatial gene distribution and anatomical localization at the whole mouse brain level. Our subsequent spatial correlation analysis systematically compared the PPs to known anatomical regions and ontology from the Allen Mouse Brain Atlas using spatial neighborhoods. We demonstrate that our stable and spatially coherent PPs, whose linear combinations accurately approximate the spatial gene data, are highly correlated with combinations of expert-annotated brain regions. These PPs yield a brain ontology based purely on spatial gene expression. Our PP identification approach outperforms principal component analysis and typical clustering algorithms on the same task. Moreover, we show that the stable PPs reveal marked regional imbalance of brainwide genetic architecture, leading to region-specific marker genes and gene coexpression networks. Our findings highlight the advantages of stability-driven machine learning for plausible biological discovery from dense spatial gene expression data, streamlining tasks that are infeasible by conventional manual approaches.

Repairing the in situ hybridization missing data in the hippocampus region by using a 3D residual U-Net model

The hippocampus is a critical brain region. Transcriptome data provides valuable insights into the structure and function of the hippocampus at the gene level. However, transcriptome data is often incomplete. To address this issue, we use the convolutional neural network model to repair the missing voxels in the hippocampus region, based on Allen institute coronal slices in situ hybridization (ISH) dataset. Moreover, we analyze the gene expression correlation between coronal and sagittal dataset in the hippocampus region. The results demonstrated that the trend of gene expression correlation between the coronal and sagittal datasets remained consistent following the repair of missing data in the coronal ISH dataset. In the last, we use repaired ISH dataset to identify novel genes specific to hippocampal subregions. Our findings demonstrate the accuracy and effectiveness of using deep learning method to repair ISH missing data. After being repaired, ISH has the potential to improve our comprehension of the hippocampus's structure and function.

DeepSlice: rapid fully automatic registration of mouse brain imaging to a volumetric atlas

Registration of data to a common frame of reference is an essential step in the analysis and integration of diverse neuroscientific data. To this end, volumetric brain atlases enable histological datasets to be spatially registered and analyzed, yet accurate registration remains expertise-dependent and slow. In order to address this limitation, we have trained a neural network, DeepSlice, to register mouse brain histological images to the Allen Brain Common Coordinate Framework, retaining registration accuracy while improving speed by >1000 fold.

Endothelial EphrinB2 Regulates Sunitinib Therapy Response in Murine Glioma

Vascular guidance is critical in developmental vasculogenesis and pathological angiogenesis. Brain tumors are strongly vascularized, and antiangiogenic therapy was anticipated to exhibit a strong anti-tumor effect in this tumor type. However, vascular endothelial growth factor A (VEGFA) specific inhibition had no significant impact in clinical practice of gliomas. More research is needed to understand the failure of this therapeutic approach. EphrinB2 has been found to directly interact with vascular endothelial growth factor receptor 2 (VEGFR2) and regulate its activity. Here we analyzed the expression of ephrinB2 and EphB4 in human glioma, we observed vascular localization of ephrinB2 in physiology and pathology and found a significant survival reduction in patients with elevated ephrinB2 tumor expression. Induced endothelial specific depletion of ephrinB2 in the adult mouse (efnb2i∆EC) had no effect on the quiescent vascular system of the brain. However, we found glioma growth and perfusion altered in efnb2i∆EC animals similar to the effects observed with antiangiogenic therapy. No additional anti-tumor effect was observed in efnb2i∆EC animals treated with antiangiogenic therapy. Our data indicate that ephrinB2 and VEGFR2 converge on the same pathway and intervention with either molecule results in a reduction in angiogenesis.

Assessing the replicability of spatial gene expression using atlas data from the adult mouse brain

High-throughput, spatially resolved gene expression techniques are poised to be transformative across biology by overcoming a central limitation in single-cell biology: the lack of information on relationships that organize the cells into the functional groupings characteristic of tissues in complex multicellular organisms. Spatial expression is particularly interesting in the mammalian brain, which has a highly defined structure, strong spatial constraint in its organization, and detailed multimodal phenotypes for cells and ensembles of cells that can be linked to mesoscale properties such as projection patterns, and from there, to circuits generating behavior. However, as with any type of expression data, cross-dataset benchmarking of spatial data is a crucial first step. Here, we assess the replicability, with reference to canonical brain subdivisions, between the Allen Institute's in situ hybridization data from the adult mouse brain (Allen Brain Atlas (ABA)) and a similar dataset collected using spatial transcriptomics (ST). With the advent of tractable spatial techniques, for the first time, we are able to benchmark the Allen Institute's whole-brain, whole-transcriptome spatial expression dataset with a second independent dataset that similarly spans the whole brain and transcriptome. We use regularized linear regression (LASSO), linear regression, and correlation-based feature selection in a supervised learning framework to classify expression samples relative to their assayed location. We show that Allen Reference Atlas labels are classifiable using transcription in both data sets, but that performance is higher in the ABA than in ST. Furthermore, models trained in one dataset and tested in the opposite dataset do not reproduce classification performance bidirectionally. While an identifying expression profile can be found for a given brain area, it does not generalize to the opposite dataset. In general, we found that canonical brain area labels are classifiable in gene expression space within dataset and that our observed performance is not merely reflecting physical distance in the brain. However, we also show that cross-platform classification is not robust. Emerging spatial datasets from the mouse brain will allow further characterization of cross-dataset replicability ultimately providing a valuable reference set for understanding the cell biology of the brain.

Overcoming false-positive gene-category enrichment in the analysis of spatially resolved transcriptomic brain atlas data

Transcriptomic atlases have improved our understanding of the correlations between gene-expression patterns and spatially varying properties of brain structure and function. Gene-category enrichment analysis (GCEA) is a common method to identify functional gene categories that drive these associations, using gene-to-category annotation systems like the Gene Ontology (GO). Here, we show that applying standard GCEA methodology to spatial transcriptomic data is affected by substantial false-positive bias, with GO categories displaying an over 500-fold average inflation of false-positive associations with random neural phenotypes in mouse and human. The estimated false-positive rate of a GO category is associated with its rate of being reported as significantly enriched in the literature, suggesting that published reports are affected by this false-positive bias. We show that within-category gene-gene coexpression and spatial autocorrelation are key drivers of the false-positive bias and introduce flexible ensemble-based null models that can account for these effects, made available as a software toolbox.

Differential Spatiotemporal Expression of Type I and Type II Cadherins Associated With the Segmentation of the Central Nervous System and Formation of Brain Nuclei in the Developing Mouse

Type I and type II classical cadherins comprise a family of cell adhesion molecules that regulate cell sorting and tissue separation by forming specific homo and heterophilic bonds. Factors that affect cadherin-mediated cell-cell adhesion include cadherin binding affinity and expression level. This study examines the expression pattern of type I cadherins (Cdh1, Cdh2, Cdh3, and Cdh4), type II cadherins (Cdh6, Cdh7, Cdh8, Cdh9, Cdh10, Cdh11, Cdh12, Cdh18, Cdh20, and Cdh24), and the atypical cadherin 13 (Cdh13) during distinct morphogenetic events in the developing mouse central nervous system from embryonic day 11.5 to postnatal day 56. Cadherin mRNA expression levels obtained from in situ hybridization experiments carried out at the Allen Institute for Brain Science (https://alleninstitute.org/) were retrieved from the Allen Developing Mouse Brain Atlas. Cdh2 is the most abundantly expressed type I cadherin throughout development, while Cdh1, Cdh3, and Cdh4 are expressed at low levels. Type II cadherins show a dynamic pattern of expression that varies between neuroanatomical structures and developmental ages. Atypical Cdh13 expression pattern correlates with Cdh2 in abundancy and localization. Analyses of cadherin-mediated relative adhesion estimated from their expression level and binding affinity show substantial differences in adhesive properties between regions of the neural tube associated with the segmentation along the anterior–posterior axis. Differences in relative adhesion were also observed between brain nuclei in the developing subpallium (basal ganglia), suggesting that differential cell adhesion contributes to the segregation of neuronal pools. In the adult cerebral cortex, type II cadherins Cdh6, Cdh8, Cdh10, and Cdh12 are abundant in intermediate layers, while Cdh11 shows a gradated expression from the deeper layer 6 to the superficial layer 1, and Cdh9, Cdh18, and Cdh24 are more abundant in the deeper layers. Person’s correlation analyses of cadherins mRNA expression patterns between areas and layers of the cerebral cortex and the nuclei of the subpallium show significant correlations between certain cortical areas and the basal ganglia. The study shows that differential cadherin expression and cadherin-mediated adhesion are associated with a wide range of morphogenetic events in the developing central nervous system including the organization of neurons into layers, the segregation of neurons into nuclei, and the formation of neuronal circuits.

Constructing and optimizing 3D atlases from 2D data with application to the developing mouse brain

3D imaging data necessitate 3D reference atlases for accurate quantitative interpretation. Existing computational methods to generate 3D atlases from 2D-derived atlases result in extensive artifacts, while manual curation approaches are labor-intensive. We present a computational approach for 3D atlas construction that substantially reduces artifacts by identifying anatomical boundaries in the underlying imaging data and using these to guide 3D transformation. Anatomical boundaries also allow extension of atlases to complete edge regions. Applying these methods to the eight developmental stages in the Allen Developing Mouse Brain Atlas (ADMBA) led to more comprehensive and accurate atlases. We generated imaging data from 15 whole mouse brains to validate atlas performance and observed qualitative and quantitative improvement (37% greater alignment between atlas and anatomical boundaries). We provide the pipeline as the MagellanMapper software and the eight 3D reconstructed ADMBA atlases. These resources facilitate whole-organ quantitative analysis between samples and across development.

The research community needs precise, reliable 3D atlases of organs to pinpoint where biological structures and processes are located. For instance, these maps are essential to understand where specific genes are turned on or off, or the spatial organization of various groups of cells over time.

For centuries, atlases have been built by thinly ‘slicing up’ an organ, and then precisely representing each 2D layer. Yet this approach is imperfect: each layer may be accurate on its own, but inevitable mismatches appear between the slices when viewed in 3D or from another angle.

Advances in microscopy now allow entire organs to be imaged in 3D. Comparing these images with atlases could help to detect subtle differences that indicate or underlie disease. However, this is only possible if 3D maps are accurate and do not feature mismatches between layers. To create an atlas without such artifacts, one approach consists in starting from scratch and manually redrawing the maps in 3D, a labor-intensive method that discards a large body of well-established atlases.

Instead, Young et al. set out to create an automated method which could help to refine existing ‘layer-based’ atlases, releasing software that anyone can use to improve current maps. The package was created by harnessing eight atlases in the Allen Developing Mouse Brain Atlas, and then using the underlying anatomical images to resolve discrepancies between layers or fill out any missing areas. Known as MagellanMapper, the software was extensively tested to demonstrate the accuracy of the maps it creates, including comparison to whole-brain imaging data from 15 mouse brains. Armed with this new software, researchers can improve the accuracy of their atlases, helping them to understand the structure of organs at the level of the cell and giving them insight into a broad range of human disorders.

Automated identification of the mouse brain's spatial compartments from in situ sequencing data

Background

Neuroanatomical compartments of the mouse brain are identified and outlined mainly based on manual annotations of samples using features related to tissue and cellular morphology, taking advantage of publicly available reference atlases. However, this task is challenging since sliced tissue sections are rarely perfectly parallel or angled with respect to sections in the reference atlas and organs from different individuals may vary in size and shape and requires manual annotation. With the advent of in situ sequencing technologies and automated approaches, it is now possible to profile the gene expression of targeted genes inside preserved tissue samples and thus spatially map biological processes across anatomical compartments.

Results

Here, we show how in situ sequencing data combined with dimensionality reduction and clustering can be used to identify spatial compartments that correspond to known anatomical compartments of the brain. We also visualize gradients in gene expression and sharp as well as smooth transitions between different compartments. We apply our method on mouse brain sections and show that a fully unsupervised approach can computationally define anatomical compartments, which are highly reproducible across individuals, using as few as 18 gene markers. We also show that morphological variation does not always follow gene expression, and different spatial compartments can be defined by various cell types with common morphological features but distinct gene expression profiles.

Conclusion

We show that spatial gene expression data can be used for unsupervised and unbiased annotations of mouse brain spatial compartments based only on molecular markers, without the need of subjective manual annotations based on tissue and cell morphology or matching reference atlases.

Whole-Brain Image Analysis and Anatomical Atlas 3D Generation Using MagellanMapper

MagellanMapper is a software suite designed for visual inspection and end-to-end automated processing of large-volume, 3D brain imaging datasets in a memory-efficient manner. The rapidly growing number of large-volume, high-resolution datasets necessitates visualization of raw data at both macro- and microscopic levels to assess the quality of data, as well as automated processing to quantify data in an unbiased manner for comparison across a large number of samples. To facilitate these analyses, MagellanMapper provides both a graphical user interface for manual inspection and a command-line interface for automated image processing. At the macroscopic level, the graphical interface allows researchers to view full volumetric images simultaneously in each dimension and to annotate anatomical label placements. At the microscopic level, researchers can inspect regions of interest at high resolution to build ground truth data of cellular locations such as nuclei positions. Using the command-line interface, researchers can automate cell detection across volumetric images, refine anatomical atlas labels to fit underlying histology, register these atlases to sample images, and perform statistical analyses by anatomical region. MagellanMapper leverages established open-source computer vision libraries and is itself open source and freely available for download and extension. © 2020 Wiley Periodicals LLC. Basic Protocol 1: MagellanMapper installation Alternate Protocol: Alternative methods for MagellanMapper installation Basic Protocol 2: Import image files into MagellanMapper Basic Protocol 3: Region of interest visualization and annotation Basic Protocol 4: Explore an atlas along all three dimensions and register to a sample brain Basic Protocol 5: Automated 3D anatomical atlas construction Basic Protocol 6: Whole-tissue cell detection and quantification by anatomical label Support Protocol: Import a tiled microscopy image in proprietary format into MagellanMapper.

Whole-Brain Profiling of Cells and Circuits in Mammals by Tissue Clearing and Light-Sheet Microscopy

Tissue clearing and light-sheet microscopy have a 100-year-plus history, yet these fields have been combined only recently to facilitate novel experiments and measurements in neuroscience. Since tissue-clearing methods were first combined with modernized light-sheet microscopy a decade ago, the performance of both technologies has rapidly improved, broadening their applications. Here, we review the state of the art of tissue-clearing methods and light-sheet microscopy and discuss applications of these techniques in profiling cells and circuits in mice. We examine outstanding challenges and future opportunities for expanding these techniques to achieve brain-wide profiling of cells and circuits in primates and humans. Such integration will help provide a systems-level understanding of the physiology and pathology of our central nervous system.

An Application of Manifold Learning in Global Shape Descriptors

With the rapid expansion of applied 3D computational vision, shape descriptors have become increasingly important for a wide variety of applications and objects from molecules to planets. Appropriate shape descriptors are critical for accurate (and efficient) shape retrieval and 3D model classification. Several spectral-based shape descriptors have been introduced by solving various physical equations over a 3D surface model. In this paper, for the first time, we incorporate a specific manifold learning technique, introduced in statistics and machine learning, to develop a global, spectral-based shape descriptor in the computer graphics domain. The proposed descriptor utilizes the Laplacian Eigenmap technique in which the Laplacian eigenvalue problem is discretized using an exponential weighting scheme. As a result, our descriptor eliminates the limitations tied to the existing spectral descriptors, namely dependency on triangular mesh representation and high intra-class quality of 3D models. We also present a straightforward normalization method to obtain a scale-invariant and noise-resistant descriptor. The extensive experiments performed in this study using two standard 3D shape benchmarks—high-resolution TOSCA and McGill datasets—demonstrate that the present contribution provides a highly discriminative and robust shape descriptor under the presence of a high level of noise, random scale variations, and low sampling rate, in addition to the known isometric-invariance property of the Laplace–Beltrami operator. The proposed method significantly outperforms state-of-the-art spectral descriptors in shape retrieval and classification. The proposed descriptor is limited to closed manifolds due to its inherited inability to accurately handle manifolds with boundaries.

Discovering Conserved Properties of Brain Organization Through Multimodal Integration and Interspecies Comparison

The primate cerebral cortex is broadly organized along hierarchical processing streams underpinned by corresponding variation in the brain's microstructure and interareal connectivity patterns. Fulcher et al. recently demonstrated that a similar organization exists in the mouse cortex by combining independent datasets of cytoarchitecture, gene expression, cell densities, and long-range axonal connectivity. Using the T1w:T2w magnetic resonance imaging map as a common spatial reference for data-driven comparison of cortical gradients between mouse and human, we highlighted a common hierarchical expression pattern of numerous brain-related genes, providing new understanding of how systematic structural variation shapes functional specialization in mammalian brains. Reflecting on these findings, here we discuss how open neuroscience datasets, combined with advanced neuroinformatics approaches, will be crucial in the ongoing search for organization principles of brain structure. We explore the promises and challenges of integrative studies and argue that a tighter collaboration between experimental, statistical, and theoretical neuroscientists is needed to drive progress further.

Geometry processing of conventionally produced mouse brain slice images

BACKGROUND

Brain mapping research in most neuroanatomical laboratories relies on conventional processing techniques, which often introduce histological artifacts such as tissue tears and tissue loss.

NEW METHOD

In this paper, we present techniques and algorithms for automatic registration and 3D reconstruction of conventionally produced mouse brain slices in a standardized atlas space. This is achieved first by constructing a virtual 3D mouse brain model from annotated slices of Allen Reference Atlas (ARA). Virtual re-slicing of the reconstructed model generates ARA-based slice images corresponding to the microscopic images of histological brain sections. These image pairs are aligned using a geometric approach through contour images. Histological artifacts in the microscopic images are detected and removed using Constrained Delaunay Triangulation before performing global alignment. Finally, non-linear registration is performed by solving Laplace's equation with Dirichlet boundary conditions.

RESULTS

Our methods provide significant improvements over previously reported registration techniques for the tested slices in 3D space, especially on slices with significant histological artifacts. Further, as one of the application we count the number of neurons in various anatomical regions using a dataset of 51 microscopic slices from a single mouse brain.

COMPARISON WITH EXISTING METHOD(S)

To the best of our knowledge the presented work is the first that automatically registers both clean as well as highly damaged high-resolutions histological slices of mouse brain to a 3D annotated reference atlas space.

CONCLUSIONS

This work represents a significant contribution to this subfield of neuroscience as it provides tools to neuroanatomist for analyzing and processing histological data.

The BRAIN Initiative Cell Census Consortium: Lessons Learned toward Generating a Comprehensive Brain Cell Atlas

A comprehensive characterization of neuronal cell types, their distributions, and patterns of connectivity is critical for understanding the properties of neural circuits and how they generate behaviors. Here we review the experiences of the BRAIN Initiative Cell Census Consortium, ten pilot projects funded by the U.S. BRAIN Initiative, in developing, validating, and scaling up emerging genomic and anatomical mapping technologies for creating a complete inventory of neuronal cell types and their connections in multiple species and during development. These projects lay the foundation for a larger and longer-term effort to generate whole-brain cell atlases in species including mice and humans.

Molecular fingerprinting of principal neurons in the rodent hippocampus: A neuroinformatics approach

Neurons are often classified by their morphological and molecular properties. The online knowledge base Hippocampome.org primarily defines neuron types from the rodent hippocampal formation based on their main neurotransmitter (glutamate or GABA) and the spatial distributions of their axons and dendrites. For each neuron type, this open-access resource reports any and all published information regarding the presence or absence of known molecular markers, including calcium-binding proteins, neuropeptides, receptors, channels, transcription factors, and other molecules of biomedical relevance. The resulting chemical profile is relatively sparse: even for the best studied neuron types, the expression or lack thereof of fewer than 70 molecules has been firmly established to date. The mouse genome-wide in situ hybridization mapping of the Allen Brain Atlas provides a wealth of data that, when appropriately analyzed, can substantially augment the molecular marker knowledge in Hippocampome.org. Here we focus on the principal cell layers of dentate gyrus (DG), CA3, CA2, and CA1, which together contain approximately 90% of hippocampal neurons. These four anatomical parcels are densely packed with somata of mostly excitatory projection neurons. Thus, gene expression data for those layers can be justifiably linked to the respective principal neuron types: granule cells in DG and pyramidal cells in CA3, CA2, and CA1. In order to enable consistent interpretation across genes and regions, we screened the whole-genome dataset against known molecular markers of those neuron types. The resulting threshold values allow over 6000 very-high confidence (>99.5%) expressed/not-expressed assignments, expanding the biochemical information content of Hippocampome.org more than five-fold. Many of these newly identified molecular markers are potential pharmacological targets for major neurological and psychiatric conditions. Furthermore, our approach yields reasonable expression/non-expression estimates for every single gene in each of these four neuron types with >90% average confidence, providing a considerably complete genetic characterization of hippocampal principal neurons.

On Expression Patterns and Developmental Origin of Human Brain Regions

Anatomical substructures of the human brain have characteristic cell-types, connectivity and local circuitry, which are reflected in area-specific transcriptome signatures, but the principles governing area-specific transcription and their relation to brain development are still being studied. In adult rodents, areal transcriptome patterns agree with the embryonic origin of brain regions, but the processes and genes that preserve an embryonic signature in regional expression profiles were not quantified. Furthermore, it is not clear how embryonic-origin signatures of adult-brain expression interplay with changes in expression patterns during development. Here we first quantify which genes have regional expression-patterns related to the developmental origin of brain regions, using genome-wide mRNA expression from post-mortem adult human brains. We find that almost all human genes (92%) exhibit an expression pattern that agrees with developmental brain-region ontology, but that this agreement changes at multiple phases during development. Agreement is particularly strong in neuron-specific genes, but also in genes that are not spatially correlated with neuron-specific or glia-specific markers. Surprisingly, agreement is also stronger in early-evolved genes. We further find that pairs of similar genes having high agreement to developmental region ontology tend to be more strongly correlated or anti-correlated, and that the strength of spatial correlation changes more strongly in gene pairs with stronger embryonic signatures. These results suggest that transcription regulation of most genes in the adult human brain is spatially tuned in a way that changes through life, but in agreement with development-determined brain regions.

aMAP is a validated pipeline for registration and segmentation of high-resolution mouse brain data

The validation of automated image registration and segmentation is crucial for accurate and reliable mapping of brain connectivity and function in three-dimensional (3D) data sets. While validation standards are necessarily high and routinely met in the clinical arena, they have to date been lacking for high-resolution microscopy data sets obtained from the rodent brain. Here we present a tool for optimized automated mouse atlas propagation (aMAP) based on clinical registration software (NiftyReg) for anatomical segmentation of high-resolution 3D fluorescence images of the adult mouse brain. We empirically evaluate aMAP as a method for registration and subsequent segmentation by validating it against the performance of expert human raters. This study therefore establishes a benchmark standard for mapping the molecular function and cellular connectivity of the rodent brain.

Genome-wide coexpression of steroid receptors in the mouse brain: Identifying signaling pathways and functionally coordinated regions

Steroid receptors are pleiotropic transcription factors that coordinate adaptation to different physiological states. An important target organ is the brain, but even though their effects are well studied in specific regions, brain-wide steroid receptor targets and mediators remain largely unknown due to the complexity of the brain. Here, we tested the idea that novel aspects of steroid action can be identified through spatial correlation of steroid receptors with genome-wide mRNA expression across different regions in the mouse brain. First, we observed significant coexpression of six nuclear receptors (NRs) [androgen receptor (Ar), estrogen receptor alpha (Esr1), estrogen receptor beta (Esr2), glucocorticoid receptor (Gr), mineralocorticoid receptor (Mr), and progesterone receptor (Pgr)] with sets of steroid target genes that were identified in single brain regions. These coexpression relationships were also present in distinct other brain regions, suggestive of as yet unidentified coordinate regulation of brain regions by, for example, glucocorticoids and estrogens. Second, coexpression of a set of 62 known NR coregulators and the six steroid receptors in 12 nonoverlapping mouse brain regions revealed selective downstream pathways, such as Pak6 as a mediator for the effects of Ar and Gr on dopaminergic transmission. Third, Magel2 and Irs4 were identified and validated as strongly responsive targets to the estrogen diethylstilbestrol in the mouse hypothalamus. The brain- and genome-wide correlations of mRNA expression levels of six steroid receptors that we provide constitute a rich resource for further predictions and understanding of brain modulation by steroid hormones.

Cell-type-specific neuroanatomy of cliques of autism-related genes in the mouse brain

Two cliques of genes identified computationally for their high co-expression in the mouse brain according to the Allen Brain Atlas, and for their enrichment in genes related to autism spectrum disorder (ASD), have recently been shown to be highly co-expressed in the cerebellar cortex, compared to what could be expected by chance. Moreover, the expression of these cliques of genes is not homogeneous across the cerebellar cortex, and it has been noted that their expression pattern seems to highlight the granular layer. However, this observation was only made by eye, and recent advances in computational neuroanatomy allow to rank cell types in the mouse brain (characterized by their transcriptome profiles) according to the similarity between their spatial density profiles and the spatial expression profiles of the cliques. We establish by Monte Carlo simulation that with probability at least 99%, the expression profiles of the two cliques are more similar to the density profile of granule cells than 99% of the expression of cliques containing the same number of genes (Purkinje cells also score above 99% in one of the cliques). Thresholding the expression profiles shows that the signal is more intense in the granular layer. Finally, we work out pairs of cell types whose combined expression profiles are more similar to the expression profiles of the cliques than any single cell type. These pairs predominantly consist of one cortical pyramidal cell and one cerebellar cell (which can be either a granule cell or a Purkinje cell).

Cosplicing network analysis of mammalian brain RNA-Seq data utilizing WGCNA and Mantel correlations

Across species and tissues and especially in the mammalian brain, production of gene isoforms is widespread. While gene expression coordination has been previously described as a scale-free coexpression network, the properties of transcriptome-wide isoform production coordination have been less studied. Here we evaluate the system-level properties of cosplicing in mouse, macaque, and human brain gene expression data using a novel network inference procedure. Genes are represented as vectors/lists of exon counts and distance measures sensitive to exon inclusion rates quantifies differences across samples. For all gene pairs, distance matrices are correlated across samples, resulting in cosplicing or cotranscriptional network matrices. We show that networks including cosplicing information are scale-free and distinct from coexpression. In the networks capturing cosplicing we find a set of novel hubs with unique characteristics distinguishing them from coexpression hubs: heavy representation in neurobiological functional pathways, strong overlap with markers of neurons and neuroglia, long coding lengths, and high number of both exons and annotated transcripts. Further, the cosplicing hubs are enriched in genes associated with autism spectrum disorders. Cosplicing hub homologs across eukaryotes show dramatically increasing intronic lengths but stable coding region lengths. Shared transcription factor binding sites increase coexpression but not cosplicing; the reverse is true for splicing-factor binding sites. Genes with protein-protein interactions have strong coexpression and cosplicing. Additional factors affecting the networks include shared microRNA binding sites, spatial colocalization within the striatum, and sharing a chromosomal folding domain. Cosplicing network patterns remain relatively stable across species.

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.

Neuroinformatics of the Allen Mouse Brain Connectivity Atlas

The Allen Mouse Brain Connectivity Atlas is a mesoscale whole brain axonal projection atlas of the C57Bl/6J mouse brain. Anatomical trajectories throughout the brain were mapped into a common 3D space using a standardized platform to generate a comprehensive and quantitative database of inter-areal and cell-type-specific projections. This connectivity atlas has several desirable features, including brain-wide coverage, validated and versatile experimental techniques, a single standardized data format, a quantifiable and integrated neuroinformatics resource, and an open-access public online database (http://connectivity.brain-map.org/). Meaningful informatics data quantification and comparison is key to effective use and interpretation of connectome data. This relies on successful definition of a high fidelity atlas template and framework, mapping precision of raw data sets into the 3D reference framework, accurate signal detection and quantitative connection strength algorithms, and effective presentation in an integrated online application. Here we describe key informatics pipeline steps in the creation of the Allen Mouse Brain Connectivity Atlas and include basic application use cases.

PDM-ENLOR for segmentation of mouse brain gene expression images

Statistical shape models, such as Active Shape Models (ASMs), suffer from their inability to represent a large range of variations of a complex shape and to account for the large errors in detection of (point) landmarks. We propose a method, PDM-ENLOR (Point Distribution Model-based ENsemble of LOcal Regressors), that overcomes these limitations by locating each landmark individually using an ensemble of local regression models and appearance cues from selected landmarks. We first detect a set of reference landmarks which were selected based on their saliency during training. For each landmark, an ensemble of regressors is built. From the locations of the detected reference landmarks, each regressor infers a candidate location for that landmark using local geometric constraints, encoded by a point distribution model (PDM). The final location of that point is determined as a weighted linear combination, whose coefficients are learned from the training data, of candidates proposed by its ensemble's component regressors. We use multiple subsets of reference landmarks as explanatory variables for the component regressors to provide varying degrees of locality for the models in each ensemble. This helps our ensemble model to capture a larger range of shape variations as compared to a single PDM. We demonstrate the advantages of our method on the challenging problem of segmenting gene expression images of mouse brain. The overall mean and standard deviation of the Dice coefficient overlap over all 14 anatomical regions and all 100 test images were (88.1 ± 9.5)%.

Transcriptomic-anatomic analysis of the mouse habenula uncovers a high molecular heterogeneity among neurons in the lateral complex, while gene expression in the medial complex largely obeys subnuclear boundaries

The mammalian habenula with its medial and lateral complexes has gained much interest in recent years, while knowledge on the detailed biological functions of these nuclei is still scarce. Novel strategies to differentiate and identify habenular cell types are required. Such attempts have largely failed, most likely due to the lack of appropriate molecular markers. One important tool to approach this dilemma is available in form of the Allen Brain Atlas (ABA), which provides detailed expression patterns of many genes in the mouse brain. In the present report, ABA tools in combination with visual inspection of ISH images were used to detect transcripts, which are strongly expressed in medial (MHb) and lateral (LHb) habenular complexes. Against our expectations, most transcripts were differentially distributed throughout the LHb, disregarding boundaries of subnuclear areas. Nine distinct distribution patterns were recognized. Yet, several transcripts could not be attributed to one of these, suggesting a high diversity of neuron types in the LHb. In the MHb, in contrast, many transcripts tended to obey subnuclear boundaries. The differential distribution of others like Adcyap1, Chrna3, or Trp53i11 suggests the presence of a novel subfield adjacent to the region of the MHbVm, which now is termed intermediate field of the ventral MHb. In addition, the localizations of Amigo2, Adcyap1, and a couple of other transcripts suggest a lateral extension of the MHb, which is here, termed HbX area. Apparently, this area is composed of intermingled MHb and LHb neurons and may allow functional interaction between the both habenular complexes.

Combined neuroimaging and gene expression analysis of the genetic basis of brain plasticity indicates across species homology

Brain plasticity and memory formation depend on the expression of a large number of genes. This relationship had been studied using several experimental approaches and researchers have identified genes regulating plasticity through a variety of mechanisms. Despite this effort, a great deal remains unknown regarding the role of different genes in brain plasticity. Previous studies usually focused on specific brain structures and many of the genes influencing plasticity have yet to be identified. In this work, we integrate results of in vivo neuroimaging studies of plasticity with whole-brain gene expression data for the study of neuroplasticity. Brain regions, found in the imaging study to be involved in plasticity, are first spatially mapped to the anatomical framework of the genetic database. Feature ranking methods are then applied to identify genes that are differentially expressed in these regions. We find that many of our highly ranked genes are involved in synaptic transmission and that some of these genes have been previously associated with learning and memory. We show these results to be consistent when applying our method to gene expression data from four human subjects. Finally, by performing similar experiments in mice, we reveal significant cross species correlation in the ranking of genes. In addition to the identification of plasticity related candidate genes, our results also demonstrate the potential of data integration approaches as a tool to link high level phenomena such as learning and memory to underlying molecular mechanisms.

Anatomical characterization of Cre driver mice for neural circuit mapping and manipulation

Significant advances in circuit-level analyses of the brain require tools that allow for labeling, modulation of gene expression, and monitoring and manipulation of cellular activity in specific cell types and/or anatomical regions. Large-scale projects and individual laboratories have produced hundreds of gene-specific promoter-driven Cre mouse lines invaluable for enabling genetic access to subpopulations of cells in the brain. However, the potential utility of each line may not be fully realized without systematic whole brain characterization of transgene expression patterns. We established a high-throughput in situ hybridization (ISH), imaging and data processing pipeline to describe whole brain gene expression patterns in Cre driver mice. Currently, anatomical data from over 100 Cre driver lines are publicly available via the Allen Institute's Transgenic Characterization database, which can be used to assist researchers in choosing the appropriate Cre drivers for functional, molecular, or connectional studies of different regions and/or cell types in the brain.

Automated identification of cell-type-specific genes in the mouse brain by image computing of expression patterns

Background

Differential gene expression patterns in cells of the mammalian brain result in the morphological, connectional, and functional diversity of cells. A wide variety of studies have shown that certain genes are expressed only in specific cell-types. Analysis of cell-type-specific gene expression patterns can provide insights into the relationship between genes, connectivity, brain regions, and cell-types. However, automated methods for identifying cell-type-specific genes are lacking to date.

Results

Here, we describe a set of computational methods for identifying cell-type-specific genes in the mouse brain by automated image computing of in situ hybridization (ISH) expression patterns. We applied invariant image feature descriptors to capture local gene expression information from cellular-resolution ISH images. We then built image-level representations by applying vector quantization on the image descriptors. We employed regularized learning methods for classifying genes specifically expressed in different brain cell-types. These methods can also rank image features based on their discriminative power. We used a data set of 2,872 genes from the Allen Brain Atlas in the experiments. Results showed that our methods are predictive of cell-type-specificity of genes. Our classifiers achieved AUC values of approximately 87% when the enrichment level is set to 20. In addition, we showed that the highly-ranked image features captured the relationship between cell-types.

Conclusions

Overall, our results showed that automated image computing methods could potentially be used to identify cell-type-specific genes in the mouse brain.

A high-resolution spatiotemporal atlas of gene expression of the developing mouse brain

To provide a temporal framework for the genoarchitecture of brain development, we generated in situ hybridization data for embryonic and postnatal mouse brain at seven developmental stages for ∼2,100 genes, which were processed with an automated informatics pipeline and manually annotated. This resource comprises 434,946 images, seven reference atlases, an ontogenetic ontology, and tools to explore coexpression of genes across neurodevelopment. Gene sets coinciding with developmental phenomena were identified. A temporal shift in the principles governing the molecular organization of the brain was detected, with transient neuromeric, plate-based organization of the brain present at E11.5 and E13.5. Finally, these data provided a transcription factor code that discriminates brain structures and identifies the developmental age of a tissue, providing a foundation for eventual genetic manipulation or tracking of specific brain structures over development. The resource is available as the Allen Developing Mouse Brain Atlas (http://developingmouse.brain-map.org).

Allen Brain Atlas-Driven Visualizations: a web-based gene expression energy visualization tool

The Allen Brain Atlas-Driven Visualizations (ABADV) is a publicly accessible web-based tool created to retrieve and visualize expression energy data from the Allen Brain Atlas (ABA) across multiple genes and brain structures. Though the ABA offers their own search engine and software for researchers to view their growing collection of online public data sets, including extensive gene expression and neuroanatomical data from human and mouse brain, many of their tools limit the amount of genes and brain structures researchers can view at once. To complement their work, ABADV generates multiple pie charts, bar charts and heat maps of expression energy values for any given set of genes and brain structures. Such a suite of free and easy-to-understand visualizations allows for easy comparison of gene expression across multiple brain areas. In addition, each visualization links back to the ABA so researchers may view a summary of the experimental detail. ABADV is currently supported on modern web browsers and is compatible with expression energy data from the Allen Mouse Brain Atlas in situ hybridization data. By creating this web application, researchers can immediately obtain and survey numerous amounts of expression energy data from the ABA, which they can then use to supplement their work or perform meta-analysis. In the future, we hope to enable ABADV across multiple data resources.

Altered gene expression in the dorsolateral prefrontal cortex of individuals with schizophrenia

The underlying pathology of schizophrenia (SZ) is likely as heterogeneous as its symptomatology. A variety of cortical and subcortical regions, including the prefrontal cortex, have been implicated in its pathology, and a number of genes have been identified as risk factors for disease development. We used in situ hybridization (ISH) to examine the expression of 58 genes in the dorsolateral prefrontal cortex (DLPFC, comprised of Brodmann areas 9 and 46) from 19 individuals with a premorbid diagnosis of SZ and 33 control individuals. Genes were selected based on: (1) previous identification as risk factors for SZ; (2) cell type markers or (3) laminar markers. Cell density and staining intensity were compared in the DLPFC, as well as separately in Brodmann areas 9 and 46. The expression patterns of a variety of genes, many of which are associated with the GABAergic system, were altered in SZ when compared with controls. Additional genes, including C8orf79 and NR4A2, showed alterations in cell density or staining intensity between the groups, highlighting the need for additional studies. Alterations were, with only a few exceptions, limited to Brodmann area 9, suggesting regional specificity of pathology in the DLPFC. Our results agree with previous studies on the GABAergic involvement in SZ, and suggest that areas 9 and 46 may be differentially affected in the disease. This study also highlights additional genes that may be altered in SZ, and indicates that these potentially interesting genes can be identified by ISH and high-throughput image analysis techniques.

Cell-type-based model explaining coexpression patterns of genes in the brain

Spatial patterns of gene expression in the vertebrate brain are not independent, as pairs of genes can exhibit complex patterns of coexpression. Two genes may be similarly expressed in one region, but differentially expressed in other regions. These correlations have been studied quantitatively, particularly for the Allen Atlas of the adult mouse brain, but their biological meaning remains obscure. We propose a simple model of the coexpression patterns in terms of spatial distributions of underlying cell types and establish its plausibility using independently measured cell-type-specific transcriptomes. The model allows us to predict the spatial distribution of cell types in the mouse brain.

Coexpression and cosplicing network approaches for the study of mammalian brain transcriptomes

Next-generation sequencing experiments have demonstrated great potential for transcriptome profiling. While transcriptome sequencing greatly increases the level of biological detail, system-level analysis of these high-dimensional datasets is becoming essential. We illustrate gene network approaches to the analysis of transcriptional data, with particular focus on the advantage of RNA-Seq technology compared to microarray platforms. We introduce a novel methodology for constructing cosplicing networks, based on distance measures combined with matrix correlations. We find that the cosplicing network is distinct and complementary to the coexpression network, although it shares the scale-free properties. In the cosplicing network, we find a set of novel hubs that have unique characteristics distinguishing them from coexpression hubs: they are heavily represented in neurobiological functional pathways and have strong overlap with markers of neurons and neuroglia, long-coding lengths, and high number of both exons and annotated transcripts. We also find that gene networks are plastic in the face of genetic and environmental pressures.

Interactions between the neuromodulatory systems and the amygdala: exploratory survey using the Allen Mouse Brain Atlas

Neuromodulatory systems originate in nuclei localized in the subcortical region of the brain and control fundamental behaviors by interacting with many areas of the central nervous system. An exploratory survey of the cholinergic, dopaminergic, noradrenergic, and serotonergic receptor expression energy in the amygdala, and in the neuromodulatory areas themselves was undertaken using the Allen Mouse Brain Atlas. The amygdala was chosen because of its importance in cognitive behavior and its bidirectional interaction with the neuromodulatory systems. The gene expression data of 38 neuromodulatory receptor subtypes were examined across 13 brain regions. The substantia innominata of the basal forebrain and regions of the amygdala had the highest amount of receptor expression energy for all four neuromodulatory systems examined. The ventral tegmental area also displayed high receptor expression of all four neuromodulators. In contrast, the locus coeruleus displayed low receptor expression energy overall. In general, cholinergic receptor expression was an order of magnitude greater than other neuromodulatory receptors. Since the nuclei of these neuromodulatory systems are thought to be the source of specific neurotransmitters, the projections from these nuclei to target regions may be inferred by receptor expression energy. The comprehensive analysis revealed many connectivity relations and receptor localization that had not been previously reported. The methodology presented here may be applied to other neural systems with similar characteristics, and to other animal models as these brain atlases become available.

A probabilistic latent semantic analysis model for coclustering the mouse brain atlas

The mammalian brain contains cells of a large variety of types. The phenotypic properties of cells of different types are largely the results of distinct gene expression patterns. Therefore, it is of critical importance to characterize the gene expression patterns in the mammalian brain. The Allen Developing Mouse Brain Atlas provides spatiotemporal in situ hybridization gene expression data across multiple stages of mouse brain development. It provides a framework to explore spatiotemporal regulation of gene expression during development. We employ a graph approximation formulation to cocluster the genes and the brain voxels simultaneously for each time point. We show that this formulation can be expressed as a probabilistic latent semantic analysis (PLSA) model, thereby allowing us to use the expectation-maximization algorithm for PLSA to estimate the coclustering parameters. To provide a quantitative comparison with prior methods, we evaluate the coclustering method on a set of standard synthetic data sets. Results indicate that our method consistently outperforms prior methods. We apply our method to cocluster the Allen Developing Mouse Brain Atlas data. Results indicate that our clustering of voxels is more consistent with classical neuroanatomy than those of prior methods. Our analysis also yields sets of genes that are co-expressed in a subset of the brain voxels.

Neuron-Enriched Gene Expression Patterns are Regionally Anti-Correlated with Oligodendrocyte-Enriched Patterns in the Adult Mouse and Human Brain

An important goal in neuroscience is to understand gene expression patterns in the brain. The recent availability of comprehensive and detailed expression atlases for mouse and human creates opportunities to discover global patterns and perform cross-species comparisons. Recently we reported that the major source of variation in gene transcript expression in the adult normal mouse brain can be parsimoniously explained as reflecting regional variation in glia to neuron ratios, and is correlated with degree of connectivity and location in the brain along the anterior-posterior axis. Here we extend this investigation to two gene expression assays of adult normal human brains that consisted of over 300 brain region samples, and perform comparative analyses of brain-wide expression patterns to the mouse. We performed principal components analysis (PCA) on the regional gene expression of the adult human brain to identify the expression pattern that has the largest variance. As in the mouse, we observed that the first principal component is composed of two anti-correlated patterns enriched in oligodendrocyte and neuron markers respectively. However, we also observed interesting discordant patterns between the two species. For example, a few mouse neuron markers show expression patterns that are more correlated with the human oligodendrocyte-enriched pattern and vice-versa. In conclusion, our work provides insights into human brain function and evolution by probing global relationships between regional cell type marker expression patterns in the human and mouse brain.

A new sparse simplex model for brain anatomical and genetic network analysis

The Allen Brain Atlas (ABA) database provides comprehensive 3D atlas of gene expression in the adult mouse brain for studying the spatial expression patterns in the mammalian central nervous system. It is computationally challenging to construct the accurate anatomical and genetic networks using the ABA 4D data. In this paper, we propose a novel sparse simplex model to accurately construct the brain anatomical and genetic networks, which are important to reveal the brain spatial expression patterns. Our new approach addresses the shift-invariant and parameter tuning problems, which are notorious in the existing network analysis methods, such that the proposed model is more suitable for solving practical biomedical problems. We validate our new model using the 4D ABA data, and the network construction results show the superior performance of the proposed sparse simplex model.

Allen Brain Atlas: an integrated spatio-temporal portal for exploring the central nervous system

The Allen Brain Atlas (http://www.brain-map.org) provides a unique online public resource integrating extensive gene expression data, connectivity data and neuroanatomical information with powerful search and viewing tools for the adult and developing brain in mouse, human and non-human primate. Here, we review the resources available at the Allen Brain Atlas, describing each product and data type [such as in situ hybridization (ISH) and supporting histology, microarray, RNA sequencing, reference atlases, projection mapping and magnetic resonance imaging]. In addition, standardized and unique features in the web applications are described that enable users to search and mine the various data sets. Features include both simple and sophisticated methods for gene searches, colorimetric and fluorescent ISH image viewers, graphical displays of ISH, microarray and RNA sequencing data, Brain Explorer software for 3D navigation of anatomy and gene expression, and an interactive reference atlas viewer. In addition, cross data set searches enable users to query multiple Allen Brain Atlas data sets simultaneously. All of the Allen Brain Atlas resources can be accessed through the Allen Brain Atlas data portal.

neuroVIISAS: approaching multiscale simulation of the rat connectome

neuroVIISAS is a generic platform which allows the integration of neuroontologies, mapping functions for brain atlas development, and connectivity data administration; all of which are required for the analysis of structurally and neurobiologically realistic simulations of networks. What makes neuroVIISAS unique is the ability to integrate neuroontologies, image stacks, mappings, visualizations, analyzes and simulations to use them for modelling and simulations. Based on the analysis of over 2020 tracing studies, atlas terminologies and registered histological stacks of images, neuroVIISAS permits the definition of neurobiologically realistic networks that are transferred to the simulation engine NEST. The analysis on a local and global level, the visualization of connectivity data and the results of simulations offer new possibilities to study structural and functional relationships of neural networks. This paper describes the major components and techniques of how to analyse, visualize and simulate with neuroVIISAS shown on a model network at a coarse CNS level (106 regions, 1566 connections) out of 13681 regions and 134043 connections of the left and right part of the CNS. This network of major components of the left and right hemisphere has small-world properties of the Watts-Strogatz model. Furthermore, synchronized subpopulations, oscillations of rate distributions and a time shift of population activities of the left and right hemisphere were observed in the neurocomputational simulations. In summary, a generic platform has been developed that realizes data-analysis-visualization integration for the exploration of network dynamics on multiple levels.

Web-based digital gene expression atlases for the mouse

Over the past 15 years the publicly available mouse gene expression data determined by in situ hybridization have dramatically increased in scope and spatiotemporal resolution. As a consequence of resources and tools available in the post-genomic era, full transcriptomes in the mouse brain and in the mouse embryo can be studied. Here we introduce and discuss seven current databases (MAMEP, EMBRYS, GenePaint, EURExpress, EuReGene, BGEM, and GENSAT) that grant access to large collections of expression data in mouse. We review the experimental focus, coverage, data assessment, and annotation for each of these databases and the implementation of analytic tools and links to other relevant databases. We provide a user-oriented summary of how to interrogate each database.

Peripheral nervous system genes expressed in central neurons induce growth on inhibitory substrates

Trauma to the spinal cord and brain can result in irreparable loss of function. This failure of recovery is in part due to inhibition of axon regeneration by myelin and chondroitin sulfate proteoglycans (CSPGs). Peripheral nervous system (PNS) neurons exhibit increased regenerative ability compared to central nervous system neurons, even in the presence of inhibitory environments. Previously, we identified over a thousand genes differentially expressed in PNS neurons relative to CNS neurons. These genes represent intrinsic differences that may account for the PNS's enhanced regenerative ability. Cerebellar neurons were transfected with cDNAs for each of these PNS genes to assess their ability to enhance neurite growth on inhibitory (CSPG) or permissive (laminin) substrates. Using high content analysis, we evaluated the phenotypic profile of each neuron to extract meaningful data for over 1100 genes. Several known growth associated proteins potentiated neurite growth on laminin. Most interestingly, novel genes were identified that promoted neurite growth on CSPGs (GPX3, EIF2B5, RBMX). Bioinformatic approaches also uncovered a number of novel gene families that altered neurite growth of CNS neurons.

Transcripts with in silico predicted RNA structure are enriched everywhere in the mouse brain

Background

Post-transcriptional control of gene expression is mostly conducted by specific elements in untranslated regions (UTRs) of mRNAs, in collaboration with specific binding proteins and RNAs. In several well characterized cases, these RNA elements are known to form stable secondary structures. RNA secondary structures also may have major functional implications for long noncoding RNAs (lncRNAs). Recent transcriptional data has indicated the importance of lncRNAs in brain development and function. However, no methodical efforts to investigate this have been undertaken. Here, we aim to systematically analyze the potential for RNA structure in brain-expressed transcripts.

Results

By comprehensive spatial expression analysis of the adult mouse in situ hybridization data of the Allen Mouse Brain Atlas, we show that transcripts (coding as well as non-coding) associated with in silico predicted structured probes are highly and significantly enriched in almost all analyzed brain regions. Functional implications of these RNA structures and their role in the brain are discussed in detail along with specific examples. We observe that mRNAs with a structure prediction in their UTRs are enriched for binding, transport and localization gene ontology categories. In addition, after manual examination we observe agreement between RNA binding protein interaction sites near the 3’ UTR structures and correlated expression patterns.

Conclusions

Our results show a potential use for RNA structures in expressed coding as well as noncoding transcripts in the adult mouse brain, and describe the role of structured RNAs in the context of intracellular signaling pathways and regulatory networks. Based on this data we hypothesize that RNA structure is widely involved in transcriptional and translational regulatory mechanisms in the brain and ultimately plays a role in brain function.