Bayesian hierarchical modeling for time course microarray experiments

Time course microarray experiments designed to characterize the dynamic regulation of gene expression in biological systems are becoming increasingly important. One critical issue that arises when examining time course microarray data is the identification of genes that show different temporal expression patterns among biological conditions. Here we propose a Bayesian hierarchical model to incorporate important experimental factors and to account for correlated gene expression measurements over time and over different genes. A new gene selection algorithm is also presented with the model to simultaneously identify genes that show changes in expression among biological conditions, in response to time and other experimental factors of interest. The algorithm performs well in terms of the false positive and false negative rates in simulation studies. The methodology is applied to a mouse model time course experiment to correlate temporal changes in azoxymethane-induced gene expression profiles with colorectal cancer susceptibility.     

The mode-of-action by network identification (MNI) algorithm: a network biology approach for molecular target identification

This protocol details the use of the mode-of-action by network identification (MNI) algorithm to identify the gene targets of a drug treatment based on gene-expression data. Investigators might also use the MNI algorithm to identify the gene mediators of a disease or the physiological state of cells and tissues. The MNI algorithm uses a training data set of hundreds of expression profiles to construct a statistical model of gene-regulatory networks in a cell or tissue. The model describes combinatorial influences of genes on one another. The algorithm then uses the model to filter the expression profile of a particular experimental treatment and thereby distinguish the molecular targets or mediators of the treatment response from hundreds of additional genes that also exhibit expression changes. It takes approximately 1 h per run, although run time is significantly affected by the size of the genome and data set.     

Detection of perturbation phases and developmental stages in organisms from DNA microarray time series data

Available DNA microarray time series that record gene expression along the developmental stages of multicellular eukaryotes, or in unicellular organisms subject to external perturbations such as stress and diauxie, are analyzed. By pairwise comparison of the gene expression profiles on the basis of a translation-invariant and scale-invariant distance measure corresponding to least-rectangle regression, it is shown that peaks in the average distance values are noticeable and are localized around specific time points. These points systematically coincide with the transition points between developmental phases or just follow the external perturbations. This approach can thus be used to identify automatically, from microarray time series alone, the presence of external perturbations or the succession of developmental stages in arbitrary cell systems. Moreover, our results show that there is a striking similarity between the gene expression responses to these a priori very different phenomena. In contrast, the cell cycle does not involve a perturbation-like phase, but rather continuous gene expression remodeling. Similar analyses were conducted using three other standard distance measures, showing that the one we introduced was superior. Based on these findings, we set up an adapted clustering method that uses this distance measure and classifies the genes on the basis of their expression profiles within each developmental stage or between perturbation phases.     

maSigPro: a method to identify significantly differential expression profiles in time-course microarray experiments

MOTIVATION: Multi-series time-course microarray experiments are useful approaches for exploring biological processes. In this type of experiments, the researcher is frequently interested in studying gene expression changes along time and in evaluating trend differences between the various experimental groups. The large amount of data, multiplicity of experimental conditions and the dynamic nature of the experiments poses great challenges to data analysis. RESULTS: In this work, we propose a statistical procedure to identify genes that show different gene expression profiles across analytical groups in time-course experiments. The method is a two-regression step approach where the experimental groups are identified by dummy variables. The procedure first adjusts a global regression model with all the defined variables to identify differentially expressed genes, and in second a variable selection strategy is applied to study differences between groups and to find statistically significant different profiles. The methodology is illustrated on both a real and a simulated microarray dataset.     

Unbiased pattern detection in microarray data series

MOTIVATION: Following the advent of microarray technology in recent years, the challenge for biologists is to identify genes of interest from the thousands of genetic expression levels measured in each microarray experiment. In many cases the aim is to identify pattern in the data series generated by successive microarray measurements. RESULTS: Here we introduce a new method of detecting pattern in microarray data series which is independent of the nature of this pattern. Our approach provides a measure of the algorithmic compressibility of each data series. A series which is significantly compressible is much more likely to result from simple underlying mechanisms than series which are incompressible. Accordingly, the gene associated with a compressible series is more likely to be biologically significant. We test our method on microarray time series of yeast cell cycle and show that it blindly selects genes exhibiting the expected cyclic behaviour as well as detecting other forms of pattern. Our results successfully predict two independent non-microarray experimental studies.     

Mining co-regulated gene profiles for the detection of functional associations in gene expression data

MOTIVATION: Association pattern discovery (APD) methods have been successfully applied to gene expression data. They find groups of co-regulated genes in which the genes are either up- or down-regulated throughout the identified conditions. These methods, however, fail to identify similarly expressed genes whose expressions change between up- and down-regulation from one condition to another. In order to discover these hidden patterns, we propose the concept of mining co-regulated gene profiles. Co-regulated gene profiles contain two gene sets such that genes within the same set behave identically (up or down) while genes from different sets display contrary behavior. To reduce and group the large number of similar resulting patterns, we propose a new similarity measure that can be applied together with hierarchical clustering methods. RESULTS: We tested our proposed method on two well-known yeast microarray data sets. Our implementation mined the data effectively and discovered patterns of co-regulated genes that are hidden to traditional APD methods. The high content of biologically relevant information in these patterns is demonstrated by the significant enrichment of co-regulated genes with similar functions. Our experimental results show that the Mining Attribute Profile (MAP) method is an efficient tool for the analysis of gene expression data and competitive with bi-clustering techniques.     

Dynamic modeling of gene expression in prokaryotes: application to glucose-lactose diauxie in Escherichia coli

Coexpression of genes or, more generally, similarity in the expression profiles poses an unsurmountable obstacle to inferring the gene regulatory network (GRN) based solely on data from DNA microarray time series. Clustering of genes with similar expression profiles allows for a course-grained view of the GRN and a probabilistic determination of the connectivity among the clusters. We present a model for the temporal evolution of a gene cluster network which takes into account interactions of gene products with genes and, through a non-constant degradation rate, with other gene products. The number of model parameters is reduced by using polynomial functions to interpolate temporal data points. In this manner, the task of parameter estimation is reduced to a system of linear algebraic equations, thus making the computation time shorter by orders of magnitude. To eliminate irrelevant networks, we test each GRN for stability with respect to parameter variations, and impose restrictions on its behavior near the steady state. We apply our model and methods to DNA microarray time series' data collected on Escherichia coli during glucose-lactose diauxie and infer the most probable cluster network for different phases of the experiment. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11693-011-9079-2) contains supplementary material, which is available to authorized users.     

Mining gene expression profiles: an integrated implementation of kernel principal component analysis and singular value decomposition

The detection of genes that show similar profiles under different experimental conditions is often an initial step in inferring the biological significance of such genes. Visualization tools are used to identify genes with similar profiles in microarray studies. Given the large number of genes recorded in microarray experiments, gene expression data are generally displayed on a low dimensional plot, based on linear methods. However, microarray data show nonlinearity, due to high-order terms of interaction between genes, so alternative approaches, such as kernel methods, may be more appropriate. We introduce a technique that combines kernel principal component analysis (KPCA) and Biplot to visualize gene expression profiles. Our approach relies on the singular value decomposition of the input matrix and incorporates an additional step that involves KPCA. The main properties of our method are the extraction of nonlinear features and the preservation of the input variables (genes) in the output display. We apply this algorithm to colon tumor, leukemia and lymphoma datasets. Our approach reveals the underlying structure of the gene expression profiles and provides a more intuitive understanding of the gene and sample association.     

哺乳动物转录因子DNA结合谱

基因差异表达是生物发育和对刺激作出应答的分子基础,转录因子在这种基因差异表达中发挥着重要的调控作用。因此,要弄清楚转录因子调控基因差异表达的机理,就必须鉴定出它们全部的靶基因并构建其操纵的转录调控网络。对基因组DNA的序列特异性结合是转录因子调控基因转录的关键环节,因此,要鉴定转录因子的靶基因,就必须从它们与DNA相互作用的分子水平,鉴定它们能够识别并结合的全部DNA序列,即转录因子DNA结合谱。近年来随着DNA微阵列芯片和高通量DNA测序技术的产生和快速发展,出现了建立转录因子体内及体外DNA结合谱的一系列革命性的新技术,对该领域的研究带来重大影响。这些新技术主要包括建立转录因子体内DNA结合谱的染色质免疫沉淀-芯片技术(ChIP-chip)和染色质免疫沉淀-测序技术(ChIP-Seq),以及建立转录因子体外DNA结合谱的双链DNA微阵列芯片技术(dsDNA microarray)、指数富集配体系统进化-系列分析基因表达技术(SELEX-SAGE)、结合-n-测序技术(Bind-n-Seq)、多重大规模并行SELEX技术(MMP-SELEX)、凝胶迁移实验-测序技术(EMSA-Seq)和高通量测序-荧光配体互作图谱分析技术(HiTS-FLIP)。文章将对这些新技术做一综述。     

Global gene expression profiling in Escherichia coli K12. The effects of integration host factor

We have used nylon membranes spotted in duplicate with full-length polymerase chain reaction-generated products of each of the 4,290 predicted Escherichia coli K12 open reading frames (ORFs) to measure the gene expression profiles in otherwise isogenic integration host factor IHF(+) and IHF(-) strains. Our results demonstrate that random hexamer rather than 3' ORF-specific priming of cDNA probe synthesis is required for accurate measurement of gene expression levels in bacteria. This is explained by the fact that the currently available set of 4,290 unique 3' ORF-specific primers do not hybridize to each ORF with equal efficiency and by the fact that widely differing degradation rates (steady-state levels) are observed for the 25-base pair region of each message complementary to each ORF-specific primer. To evaluate the DNA microarray data reported here, we used a linear analysis of variance (ANOVA) model appropriate for our experimental design. These statistical methods allowed us to identify and appropriately correct for experimental variables that affect the reproducibility and accuracy of DNA microarray measurements and allowed us to determine the statistical significance of gene expression differences between our IHF(+) and IHF(-) strains. Our results demonstrate that small differences in gene expression levels can be accurately measured and that the significance of differential gene expression measurements cannot be assessed simply by the magnitude of the fold difference. Our statistical criteria, supported by excellent agreement between previously determined effects of IHF on gene expression and the results reported here, have allowed us to identify new genes regulated by IHF with a high degree of confidence.     

GenXHC: a probabilistic generative model for cross-hybridization compensation in high-density genome-wide microarray data

MOTIVATION: Microarray designs containing millions to hundreds of millions of probes that tile entire genomes are currently being released. Within the next 2 months, our group will release a microarray data set containing over 12,000,000 microarray measurements taken from 37 mouse tissues. A problem that will become increasingly significant in the upcoming era of genome-wide exon-tiling microarray experiments is the removal of cross-hybridization noise. We present a probabilistic generative model for cross-hybridization in microarray data and a corresponding variational learning method for cross-hybridization compensation, GenXHC, that reduces cross-hybridization noise by taking into account multiple sources for each mRNA expression level measurement, as well as prior knowledge of hybridization similarities between the nucleotide sequences of microarray probes and their target cDNAs. RESULTS: The algorithm is applied to a subset of an exon-resolution genome-wide Agilent microarray data set for chromosome 16 of Mus musculus and is found to produce statistically significant reductions in cross-hybridization noise. The denoised data is found to produce enrichment in multiple gene ontology-biological process (GO-BP) functional groups. The algorithm is found to outperform robust multi-array analysis, another method for cross-hybridization compensation.     

Integration of Steady-State and Temporal Gene Expression Data for the Inference of Gene Regulatory Networks

We develop a new regression algorithm, cMIKANA, for inference of gene regulatory networks from combinations of steady-state and time-series gene expression data. Using simulated gene expression datasets to assess the accuracy of reconstructing gene regulatory networks, we show that steady-state and time-series data sets can successfully be combined to identify gene regulatory interactions using the new algorithm. Inferring gene networks from combined data sets was found to be advantageous when using noisy measurements collected with either lower sampling rates or a limited number of experimental replicates. We illustrate our method by applying it to a microarray gene expression dataset from human umbilical vein endothelial cells (HUVECs) which combines time series data from treatment with growth factor TNF and steady state data from siRNA knockdown treatments. Our results suggest that the combination of steady-state and time-series datasets may provide better prediction of RNA-to-RNA interactions, and may also reveal biological features that cannot be identified from dynamic or steady state information alone. Finally, we consider the experimental design of genomics experiments for gene regulatory network inference and show that network inference can be improved by incorporating steady-state measurements with time-series data.