Clustering gene expression time series data using an infinite Gaussian process mixture model.

Published

Journal Article

Transcriptome-wide time series expression profiling is used to characterize the cellular response to environmental perturbations. The first step to analyzing transcriptional response data is often to cluster genes with similar responses. Here, we present a nonparametric model-based method, Dirichlet process Gaussian process mixture model (DPGP), which jointly models data clusters with a Dirichlet process and temporal dependencies with Gaussian processes. We demonstrate the accuracy of DPGP in comparison to state-of-the-art approaches using hundreds of simulated data sets. To further test our method, we apply DPGP to published microarray data from a microbial model organism exposed to stress and to novel RNA-seq data from a human cell line exposed to the glucocorticoid dexamethasone. We validate our clusters by examining local transcription factor binding and histone modifications. Our results demonstrate that jointly modeling cluster number and temporal dependencies can reveal shared regulatory mechanisms. DPGP software is freely available online at https://github.com/PrincetonUniversity/DP_GP_cluster.

Full Text

Duke Authors

Cited Authors

  • McDowell, IC; Manandhar, D; Vockley, CM; Schmid, AK; Reddy, TE; Engelhardt, BE

Published Date

  • January 16, 2018

Published In

Volume / Issue

  • 14 / 1

Start / End Page

  • e1005896 -

PubMed ID

  • 29337990

Pubmed Central ID

  • 29337990

Electronic International Standard Serial Number (EISSN)

  • 1553-7358

International Standard Serial Number (ISSN)

  • 1553-734X

Digital Object Identifier (DOI)

  • 10.1371/journal.pcbi.1005896

Language

  • eng