Automated pipeline for atlas-based annotation of gene expression patterns: Application to postnatal day 7 mouse brain

Massive amounts of image data have been collected and continue to be generated for representing cellular gene expression throughout the mouse brain. Critical to exploiting this key effort of the post-genomic era is the ability to place these data into a common spatial reference that enables rapid interactive queries, analysis, data sharing, and visualization. In this paper, we present a set of automated protocols for generating and annotating gene expression patterns suitable for the establishment of a database. The steps include imaging tissue slices, detecting cellular gene expression levels, spatial registration with an atlas, and textual annotation. Using high-throughput in situ hybridization to generate serial sets of tissues displaying gene expression, this process was applied toward the establishment of a database representing over 200 genes in the postnatal day 7 mouse brain. These data using this protocol are now well-suited for interactive comparisons, analysis, queries, and visualization.     

Selection and validation of microarray candidate genes from subregions and subnuclei of the brain

DNA array technology now allows an enormous amount of expression data to be obtained. For large-scale gene profiling enterprises, this is of course welcome. However, the scientist interested in follow-up studies of a handful of differentially expressed genes may find it hard to sift through the vast datasets to pinpoint genes with the most desirable and reliable behaviors. Here, we present the methodology we have employed to discover genes differentially expressed in the adult mouse brain. We first used Affymetrix microarrays to compare gene expression from five different brain regions: the amygdala, cerebellum, hippocampus, olfactory bulb, and periaqueductal gray. Second, we identified genes differentially expressed within three distinct amygdala subnuclei. In this case, the tissue was microdissected by laser-capture to minimize contamination from adjacent subnuclei, and extracted RNA was subjected to three rounds of linear amplification prior to hybridization to the microarrays. To select candidate genes, we developed a custom algorithm to identify those genes with the most robust changes in expression across different replicate samples. Confirmation of expression patterns with in situ hybridization uncovered further criteria to consider in the selection process.     

Microarray data integration for genome-wide analysis of human tissue-selective gene expression

Background

Microarray gene expression data are accumulating in public databases. The expression profiles contain valuable information for understanding human gene expression patterns. However, the effective use of public microarray data requires integrating the expression profiles from heterogeneous sources.

Results

In this study, we have compiled a compendium of microarray expression profiles of various human tissue samples. The microarray raw data generated in different research laboratories have been obtained and combined into a single dataset after data normalization and transformation. To demonstrate the usefulness of the integrated microarray data for studying human gene expression patterns, we have analyzed the dataset to identify potential tissue-selective genes. A new method has been proposed for genome-wide identification of tissue-selective gene targets using both microarray intensity values and detection calls. The candidate genes for brain, liver and testis-selective expression have been examined, and the results suggest that our approach can select some interesting gene targets for further experimental studies.

Conclusion

A computational approach has been developed in this study for combining microarray expression profiles from heterogeneous sources. The integrated microarray data can be used to investigate tissue-selective expression patterns of human genes.

     

Different coexpressions of arthritis-relevant genes between different body organs and different brain regions in the normal mouse population

Structural changes in different parts of the brain in rheumatoid arthritis (RA) patients have been reported. RA is not regarded as a brain disease. Body organs such as spleen and lung produce RA-relevant genes. We hypothesized that the structural changes in the brain are caused by changes of gene expression in body organs. Changes in different parts of the brain may be affected by altered gene expressions in different body organs. This study explored whether an association between gene expressions of an organ or a body part varies in different brain structures. By examining the association of the 10 most altered genes from a mouse model of spontaneous arthritis in a normal mouse population, we found two groups of gene expression patterns between five brain structures and spleen. The correlation patterns between the prefrontal cortex, nucleus accumbens, and spleen were similar, while the associations between the other three parts of the brain and spleen showed a different pattern. Among overall patterns of the associations between body organs and brain structures, spleen and lung had a similar pattern, and patterns for kidney and liver were similar. Analysis of the five additional known arthritis-relevant genes produced similar results. Analysis of 10 nonrelevant-arthritis genes did not result in a strong association of gene expression or clearly segregated patterns. Our data suggest that abnormal gene expressions in different diseased body organs may influence structural changes in different brain parts.     

Identifying differentially expressed genes in human acute leukemia and mouse brain microarray datasets utilizing QTModel

One of the essential issues in microarray data analysis is to identify differentially expressed genes (DEGs) under different experimental treatments. In this article, a statistical procedure was proposed to identify the DEGs for gene expression data with or without missing observations from microarray experiment with one- or two-treatment factors. An F statistic based on Henderson method III was constructed to test the significance of differential expression for each gene under different treatment(s) levels. The cutoff P value was adjusted to control the experimental-wise false discovery rate. A human acute leukemia dataset corrected from 38 leukemia patients was reanalyzed by the proposed method. In comparison to the results from significant analysis of microarray (SAM) and microarray analysis of variance (MAANOVA), it was indicated that the proposed method has similar performance with MAANOVA for data with one-treatment factor, but MAANOVA cannot directly handle missing data. In addition, a mouse brain dataset collected from six brain regions of two inbred strains (two-treatment factors) was reanalyzed to identify genes with distinct regional-specific expression patterns. The results showed that the proposed method could identify more distinct regional-specific expression patterns than the previous analysis of the same dataset. Moreover, a computer program was developed and incorporated in the software QTModel, which is freely available at .     

A Penalized Likelihood Approach for Bivariate Conditional Normal Models for Dynamic Co‐expression Analysis

Summary Gene co‐expressions have been widely used in the analysis of microarray gene expression data. However, the co‐expression patterns between two genes can be mediated by cellular states, as reflected by expression of other genes, single nucleotide polymorphisms, and activity of protein kinases. In this article, we introduce a bivariate conditional normal model for identifying the variables that can mediate the co‐expression patterns between two genes. Based on this model, we introduce a likelihood ratio (LR) test and a penalized likelihood procedure for identifying the mediators that affect gene co‐expression patterns. We propose an efficient computational algorithm based on iterative reweighted least squares and cyclic coordinate descent and have shown that when the tuning parameter in the penalized likelihood is appropriately selected, such a procedure has the oracle property in selecting the variables. We present simulation results to compare with existing methods and show that the LR‐based approach can perform similarly or better than the existing method of liquid association and the penalized likelihood procedure can be quite effective in selecting the mediators. We apply the proposed method to yeast gene expression data in order to identify the kinases or single nucleotide polymorphisms that mediate the co‐expression patterns between genes.     

MS2 phage ribonucleoproteins as exogenous internal control for RT-qPCR data normalization in gene expression study of developing rat brain

The most popular strategy for normalization of RT-qPCR data involves presenting them in comparison with expression of “housekeeping” genes. However, the required stable expression of the control genes is not always achievable. As an alternative, we used ribonucleoprotein phage particles as an exogenous internal control and demonstrated that this type of normalization provides a simple and reliable method for quantification in RT-qPCR experiments. Using phage-based normalization, we analyzed mRNA levels of three popular housekeeping genes coding β-actin, glyceraldehyde-3-phosphate dehydrogenase, and ribosomal protein L30 and showed high variability in their expression patterns during rat brain development, indicating that they should not be used as controls in gene expression studies of the developing brain either individually or in combination. Using phage-based controls, we showed interstrain differences and age-related changes in the expression of genes involved in proteoglycan biosynthesis and degradation in developing brain of senescenceaccelerated OXYS rats and control Wistar rats.     

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.     

Gene expression in the rodent brain is associated with its regional connectivity

The putative link between gene expression of brain regions and their neural connectivity patterns is a fundamental question in neuroscience. Here this question is addressed in the first large scale study of a prototypical mammalian rodent brain, using a combination of rat brain regional connectivity data with gene expression of the mouse brain. Remarkably, even though this study uses data from two different rodent species (due to the data limitations), we still find that the connectivity of the majority of brain regions is highly predictable from their gene expression levels-the outgoing (incoming) connectivity is successfully predicted for 73% (56%) of brain regions, with an overall fairly marked accuracy level of 0.79 (0.83). Many genes are found to play a part in predicting both the incoming and outgoing connectivity (241 out of the 500 top selected genes, p-value<1e-5). Reassuringly, the genes previously known from the literature to be involved in axon guidance do carry significant information about regional brain connectivity. Surveying the genes known to be associated with the pathogenesis of several brain disorders, we find that those associated with schizophrenia, autism and attention deficit disorder are the most highly enriched in the connectivity-related genes identified here. Finally, we find that the profile of functional annotation groups that are associated with regional connectivity in the rodent is significantly correlated with the annotation profile of genes previously found to determine neural connectivity in C. elegans (Pearson correlation of 0.24, p<1e-6 for the outgoing connections and 0.27, p<1e-5 for the incoming). Overall, the association between connectivity and gene expression in a specific extant rodent species' brain is likely to be even stronger than found here, given the limitations of current data.     

Localizing Genes to Cerebellar Layers by Classifying ISH Images

Gene expression controls how the brain develops and functions. Understanding control processes in the brain is particularly hard since they involve numerous types of neurons and glia, and very little is known about which genes are expressed in which cells and brain layers. Here we describe an approach to detect genes whose expression is primarily localized to a specific brain layer and apply it to the mouse cerebellum. We learn typical spatial patterns of expression from a few markers that are known to be localized to specific layers, and use these patterns to predict localization for new genes. We analyze images of in-situ hybridization (ISH) experiments, which we represent using histograms of local binary patterns (LBP) and train image classifiers and gene classifiers for four layers of the cerebellum: the Purkinje, granular, molecular and white matter layer. On held-out data, the layer classifiers achieve accuracy above 94% (AUC) by representing each image at multiple scales and by combining multiple image scores into a single gene-level decision. When applied to the full mouse genome, the classifiers predict specific layer localization for hundreds of new genes in the Purkinje and granular layers. Many genes localized to the Purkinje layer are likely to be expressed in astrocytes, and many others are involved in lipid metabolism, possibly due to the unusual size of Purkinje cells.     

Integrated analysis of gene expression and copy number data on gene shaving using independent component analysis

DNA microarray gene expression and microarray-based comparative genomic hybridization (aCGH) have been widely used for biomedical discovery. Because of the large number of genes and the complex nature of biological networks, various analysis methods have been proposed. One such method is "gene shaving," a procedure which identifies subsets of the genes with coherent expression patterns and large variation across samples. Since combining genomic information from multiple sources can improve classification and prediction of diseases, in this paper we proposed a new method, "ICA gene shaving" (ICA, independent component analysis), for jointly analyzing gene expression and copy number data. First we used ICA to analyze joint measurements, gene expression and copy number, of a biological system and project the data onto statistically independent biological processes. Next, we used these results to identify patterns of variation in the data and then applied an iterative shaving method. We investigated the properties of our proposed method by analyzing both simulated and real data. We demonstrated that the robustness of our method to noise using simulated data. Using breast cancer data, we showed that our method is superior to the Generalized Singular Value Decomposition (GSVD) gene shaving method for identifying genes associated with breast cancer.     

A statistical method for identifying differential gene-gene co-expression patterns

MOTIVATION: To understand cancer etiology, it is important to explore molecular changes in cellular processes from normal state to cancerous state. Because genes interact with each other during cellular processes, carcinogenesis related genes may form differential co-expression patterns with other genes in different cell states. In this study, we develop a statistical method for identifying differential gene-gene co-expression patterns in different cell states. RESULTS: For efficient pattern recognition, we extend the traditional F-statistic and obtain an Expected Conditional F-statistic (ECF-statistic), which incorporates statistical information of location and correlation. We also propose a statistical method for data transformation. Our approach is applied to a microarray gene expression dataset for prostate cancer study. For a gene of interest, our method can select other genes that have differential gene-gene co-expression patterns with this gene in different cell states. The 10 most frequently selected genes, include hepsin, GSTP1 and AMACR, which have recently been proposed to be associated with prostate carcinogenesis. However, genes GSTP1 and AMACR cannot be identified by studying differential gene expression alone. By using tumor suppressor genes TP53, PTEN and RB1, we identify seven genes that also include hepsin, GSTP1 and AMACR. We show that genes associated with cancer may have differential gene-gene expression patterns with many other genes in different cell states. By discovering such patterns, we may be able to identify carcinogenesis related genes.     

Registration of the expression patterns of Drosophila segmentation genes by two independent methods

MOTIVATION: To construct an integrated map of Drosophila segmentation gene expression from partial data taken from individual embryos. RESULTS: Spline and wavelet based registration techniques were developed to register Drosophila segmentation gene expression data. As ground control points for registration we used the locations of extrema on gene expression patterns, represented in 1D. The registration method was characterized by unprecedented high accuracy. A method for constructing the integrated pattern of gene expression at cellular resolution was designed. These patterns were constructed for 9 segmentation genes belonging to gap and pair-rule classes.     

Temporal classification of Drosophila segmentation gene expression patterns by the multi-valued neural recognition method

In order to reconstruct the establishment of the body pattern over time in Drosophila embryos, we have developed automated methods for detecting the age of an embryo on the basis of knowledge about its gene expression patterns. In this paper we perform temporal classification of confocal images of expression patterns of genes controlling segmentation by means of a neural network based on multi-valued neurons (MVN). MVN are artificial neural processing elements with complex-valued weights and high functionality, which proved to be efficient for solving the image recognition problems. The results obtained by this method confirm its efficiency for image recognition and indicate that the method can detect characteristic features of expression patterns which mark their development over time.     

Predicting and analyzing early wake-up associated gene expressions by integrating GWAS and eQTL studies

Circadian rhythms are endogenous 24-hour rhythmic oscillations affecting human behaviors, such as sleep, blood pressure and other biological processes, the disturbance of which lead to circadian rhythm sleep disorders (CRSDs). In this study, based on the data from genome-wide association studies (GWASs) and expression quantitative trait loci (eQTLs), we tried to identify novel gene expression patterns in brain tissues that were associated with early wake-up. First, the maximum-relevance-minimum-redundancy (mRMR) method was adopted to analyze the involved gene expression patterns, yielding a feature list. Second, the incremental feature selection (IFS) method and the Dagging algorithm were applied to extract important gene expression patterns, which yield the best performance for Dagging. As a result, 4374 gene expression patterns were obtained, and they were further used to build an optimal classifier with a good performance of a Matthews's correlation coefficient of 0.933. Furthermore, the most important 49 gene expression patterns were extensively analyzed. Four genes were found to be related to circadian rhythm, as reported in previous studies. As a first attempt in identifying the target genes whose expression levels are associated with sleep-wake rhythms through integrating GWAS and eQTL results, this study can motivate more investigations in this regard.This article is part of a Special Issue entitled: Accelerating Precision Medicine through Genetic and Genomic Big Data Analysis edited by Yudong Cai & Tao Huang.