Repository logo
 

Comparison of sparse biclustering algorithms for gene expression datasets.

Accepted version
Peer-reviewed

Change log

Authors

Abstract

MOTIVATION: Gene clustering and sample clustering are commonly used to find patterns in gene expression datasets. However, genes may cluster differently in heterogeneous samples (e.g. different tissues or disease states), whilst traditional methods assume that clusters are consistent across samples. Biclustering algorithms aim to solve this issue by performing sample clustering and gene clustering simultaneously. Existing reviews of biclustering algorithms have yet to include a number of more recent algorithms and have based comparisons on simplistic simulated datasets without specific evaluation of biclusters in real datasets, using less robust metrics. RESULTS: We compared four classes of sparse biclustering algorithms on a range of simulated and real datasets. All algorithms generally struggled on simulated datasets with a large number of genes or implanted biclusters. We found that Bayesian algorithms with strict sparsity constraints had high accuracy on the simulated datasets and did not require any post-processing, but were considerably slower than other algorithm classes. We found that non-negative matrix factorisation algorithms performed poorly, but could be re-purposed for biclustering through a sparsity-inducing post-processing procedure we introduce; one such algorithm was one of the most highly ranked on real datasets. In a multi-tissue knockout mouse RNA-seq dataset, the algorithms rarely returned clusters containing samples from multiple different tissues, whilst such clusters were identified in a human dataset of more closely related cell types (sorted blood cell subsets). This highlights the need for further thought in the design and analysis of multi-tissue studies to avoid differences between tissues dominating the analysis. AVAILABILITY: Code to run the analysis is available at https://github.com/nichollskc/biclust_comp, including wrappers for each algorithm, implementations of evaluation metrics, and code to simulate datasets and perform pre- and post-processing. The full tables of results are available at https://doi.org/10.5281/zenodo.4581206.

Description

Keywords

biclustering, clustering, gene expression, multi-tissue, Algorithms, Databases, Genetic, Gene Expression Profiling, Gene Expression Regulation, Oligonucleotide Array Sequence Analysis

Journal Title

Brief Bioinform

Conference Name

Journal ISSN

1467-5463
1477-4054

Volume Title

Publisher

Oxford University Press (OUP)

Rights

All rights reserved
Sponsorship
Wellcome Trust (107881/Z/15/Z)
Medical Research Council (MC_UU_00002/4)
Engineering and Physical Sciences Research Council (EP/P020259/1)
National Institute for Health and Care Research (IS-BRC-1215-20014)
CW is supported by the Wellcome Trust (WT107881) and the Medical Research Council (MC UU 00002/4). KN is supported by the Wellcome Trust (220024/Z/19/Z). CW and KN are also supported by the NIHR Cambridge BRC.
Relationships
Is derived from: