Dataset Information

Bayesian hierarchical clustering for microarray time series data with replicates and outlier measurements.

ABSTRACT:

Background

Post-genomic molecular biology has resulted in an explosion of data, providing measurements for large numbers of genes, proteins and metabolites. Time series experiments have become increasingly common, necessitating the development of novel analysis tools that capture the resulting data structure. Outlier measurements at one or more time points present a significant challenge, while potentially valuable replicate information is often ignored by existing techniques.

Results

We present a generative model-based Bayesian hierarchical clustering algorithm for microarray time series that employs Gaussian process regression to capture the structure of the data. By using a mixture model likelihood, our method permits a small proportion of the data to be modelled as outlier measurements, and adopts an empirical Bayes approach which uses replicate observations to inform a prior distribution of the noise variance. The method automatically learns the optimum number of clusters and can incorporate non-uniformly sampled time points. Using a wide variety of experimental data sets, we show that our algorithm consistently yields higher quality and more biologically meaningful clusters than current state-of-the-art methodologies. We highlight the importance of modelling outlier values by demonstrating that noisy genes can be grouped with other genes of similar biological function. We demonstrate the importance of including replicate information, which we find enables the discrimination of additional distinct expression profiles.

Conclusions

By incorporating outlier measurements and replicate values, this clustering algorithm for time series microarray data provides a step towards a better treatment of the noise inherent in measurements from high-throughput genomic technologies. Timeseries BHC is available as part of the R package 'BHC' (version 1.5), which is available for download from Bioconductor (version 2.9 and above) via http://www.bioconductor.org/packages/release/bioc/html/BHC.html?pagewanted=all.

SUBMITTER: Cooke EJ

PROVIDER: S-EPMC3228548 | biostudies-literature | 2011 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Bayesian hierarchical clustering for microarray time series data with replicates and outlier measurements.

Cooke Emma J EJ Savage Richard S RS Kirk Paul D W PD Darkins Robert R Wild David L DL

BMC bioinformatics 20111013

<h4>Background</h4>Post-genomic molecular biology has resulted in an explosion of data, providing measurements for large numbers of genes, proteins and metabolites. Time series experiments have become increasingly common, necessitating the development of novel analysis tools that capture the resulting data structure. Outlier measurements at one or more time points present a significant challenge, while potentially valuable replicate information is often ignored by existing techniques.<h4>Results ...[more]

PMID: 21995452

Similar Datasets

Project description:BackgroundBarth syndrome (BTHS) is a rare genetic disease that is characterized by cardiomyopathy, skeletal myopathy, neutropenia, and growth abnormalities and often leads to death in childhood. Recently, elamipretide has been tested as a potential first disease-modifying drug. This study aimed to identify patients with BTHS who may respond to elamipretide, based on continuous physiological measurements acquired through wearable devices.ResultsData from a randomized, double-blind, placebo-controlled crossover trial of 12 patients with BTHS were used, including physiological time series data measured using a wearable device (heart rate, respiratory rate, activity, and posture) and functional scores. The latter included the 6-minute walk test (6MWT), Patient-Reported Outcomes Measurement Information System (PROMIS) fatigue score, SWAY Balance Mobile Application score (SWAY balance score), BTHS Symptom Assessment (BTHS-SA) Total Fatigue score, muscle strength by handheld dynamometry, 5 times sit-and-stand test (5XSST), and monolysocardiolipin to cardiolipin ratio (MLCL:CL). Groups were created through median split of the functional scores into "highest score" and "lowest score", and "best response to elamipretide" and "worst response to elamipretide". Agglomerative hierarchical clustering (AHC) models were implemented to assess whether physiological data could classify patients according to functional status and distinguish non-responders from responders to elamipretide. AHC models clustered patients according to their functional status with accuracies of 60-93%, with the greatest accuracies for 6MWT (93%), PROMIS (87%), and SWAY balance score (80%). Another set of AHC models clustered patients with respect to their response to treatment with elamipretide with perfect accuracy (all 100%).ConclusionsIn this proof-of-concept study, we demonstrated that continuously acquired physiological measurements from wearable devices can be used to predict functional status and response to treatment among patients with BTHS.

Dataset Information

Bayesian hierarchical clustering for microarray time series data with replicates and outlier measurements.

Background

Results

Conclusions

Publications

Bayesian hierarchical clustering for microarray time series data with replicates and outlier measurements.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets