Detection of patient subgroups with differential expression in omics data: a comprehensive comparison of univariate measures

PLoS One. 2013 Nov 22;8(11):e79380. doi: 10.1371/journal.pone.0079380. eCollection 2013.

Abstract

Detection of yet unknown subgroups showing differential gene or protein expression is a frequent goal in the analysis of modern molecular data. Applications range from cancer biology over developmental biology to toxicology. Often a control and an experimental group are compared, and subgroups can be characterized by differential expression for only a subgroup-specific set of genes or proteins. Finding such genes and corresponding patient subgroups can help in understanding pathological pathways, diagnosis and defining drug targets. The size of the subgroup and the type of differential expression determine the optimal strategy for subgroup identification. To date, commonly used software packages hardly provide statistical tests and methods for the detection of such subgroups. Different univariate methods for subgroup detection are characterized and compared, both on simulated and on real data. We present an advanced design for simulation studies: Data is simulated under different distributional assumptions for the expression of the subgroup, and performance results are compared against theoretical upper bounds. For each distribution, different degrees of deviation from the majority of observations are considered for the subgroup. We evaluate classical approaches as well as various new suggestions in the context of omics data, including outlier sum, PADGE, and kurtosis. We also propose the new FisherSum score. ROC curve analysis and AUC values are used to quantify the ability of the methods to distinguish between genes or proteins with and without certain subgroup patterns. In general, FisherSum for small subgroups and t-test for large subgroups achieve best results. We apply each method to a case-control study on Parkinson's disease and underline the biological benefit of the new method.

Publication types

  • Research Support, Non-U.S. Gov't

MeSH terms

  • Computer Simulation*
  • Gene Expression Profiling
  • Humans
  • Models, Theoretical
  • Oligonucleotide Array Sequence Analysis

Grants and funding

This work was supported by the European Regional Development Fond (ERDF) of the European Union; the Ministerium für Innovation, Wissenschaft und Forschung des Landes Nordrhein-Westfalen (ParkCHIP, FZ 280381102); P.U.R.E. (Protein Unit for Research in Europe), a project of Nordrhein-Westfalen, a federal state of Germany, and from Collaborative Research Center SFB 876 Providing Information by Resource-Constrained Analysis, project A3. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.