Abstract:Background The increasing availability of time-series expression data opens up new possibilities to study functional linkages of genes. Present methods used to infer functional linkages between genes from expression data are mainly based on a point-to-point comparison. Change trends between consecutive time points in time-series data have been so far not well explored. Results In this work we present a new method based on extracting main features of the change trend and level of gene expression between consecutive time points. The method, termed as trend correlation (TC), includes two major steps: 1, calculating a maximal local alignment of change trend score by dynamic programming and a change trend correlation coefficient between the maximal matched change levels of each gene pair; 2, inferring relationships of gene pairs based on two statistical extraction procedures. The new method considers time shifts and inverted relationships in a similar way as the local clustering (LC) method but the latter is merely based on a point-to-point comparison. The TC method is demonstrated with data from yeast cell cycle and compared with the LC method and the widely used Pearson correlation coefficient (PCC) based clustering method. The biological significance of the gene pairs is examined with several large-scale yeast databases. Although the TC method predicts an overall lower number of gene pairs than the other two methods at a same p-value threshold, the additional number of gene pairs inferred by the TC method is considerable: e.g. 20.5% compared with the LC method and 49.6% with the PCC method for a p-value threshold of 2.7E-3. Moreover, the percentage of the inferred gene pairs consistent with databases by our method is generally higher than the LC method and similar to the PCC method. A significant number of the gene pairs only inferred by the TC method are process-identity or function-similarity pairs or have well-documented biological interactions, including 443 known protein interactions and some known cell cycle related regulatory interactions. It should be emphasized that the overlapping of gene pairs detected by the three methods is normally not very high, indicating a necessity of combining the different methods in search of functional association of genes from time-series data. For a p-value threshold of 1E-5 the percentage of process-identity and function-similarity gene pairs among the shared part of the three methods reaches 60.2% and 55.6% respectively, building a good basis for further experimental and functional study. Furthermore, the combined use of methods is important to infer more complete regulatory circuits and network as exemplified in this study. Conclusion The TC method can significantly augment the current major methods to infer functional linkages and biological network and is well suitable for exploring temporal relationships of gene expression in time-series data.

A hidden spatial-temporal Markov random field model for network-based analysis of time course gene expression data

Using Matrix of Thresholding Partial Correlation Coefficients to Infer Regulatory Network

Constructing Gene Co-Expression Networks and Predicting Functions of Unknown Genes by Random Matrix Theory

Mining gene expression data using a novel approach based on hidden Markov models.

A bayesian multivariate mixture model for high throughput spatial transcriptomics

Analyzing time-course gene expression data using profile-state hidden Markov model

Discovery of Time-Delayed Gene Regulatory Networks Based on Temporal Gene Expression Profiling

Bayesian hidden mark interaction model for detecting spatially variable genes in imaging-based spatially resolved transcriptomics data

Quantitative inference of dynamic regulatory pathways via microarray data

A statistical analysis of memory CD8 T cell differentiation: An application of a hierarchical state space model to a short time course microarray experiment

In Search of Functional Association from Time-Series Microarray Data Based on the Change Trend and Level of Gene Expression

Co-regulated Gene Module Detection for Time Series Gene Expression Data

A Markov random field-based approach to characterizing human brain development using spatial-temporal transcriptome data

Differential Markov Random Field Analysis with an Application to Detecting Differential Microbial Community Networks.

TSMiner: a novel framework for generating time-specific gene regulatory networks from time-series expression profiles

A robust statistical approach for finding informative spatially associated pathways

Predicting viral exposure response from modeling the changes of co-expression networks using time series gene expression data

A New Approach to Identify Functional Modules Using Random Matrix Theory

Gaussian Graphical Model for Identifying Significantly Responsive Regulatory Networks from Time Course High-Throughput Data

STW-MD: a novel spatio-temporal weighting and multi-step decision tree method for considering spatial heterogeneity in brain gene expression data

Bayesian Inference of Genetic Regulatory Networks from Time Series Microarray Data Using Dynamic Bayesian Networks