Abstract
In biological systems that undergo processes such as differentiation, a clear concept of progression exists. We present a novel computational approach, called Sample Progression Discovery (SPD), to discover patterns of biological progression underlying microarray gene expression data. SPD assumes that individual samples of a microarray dataset are related by an unknown biological process (i.e., differentiation, development, cell cycle, disease progression), and that each sample represents one unknown point along the progression of that process. SPD aims to organize the samples in a manner that reveals the underlying progression and to simultaneously identify subsets of genes that are responsible for that progression. We demonstrate the performance of SPD on a variety of microarray datasets that were generated by sampling a biological process at different points along its progression, without providing SPD any information of the underlying process. When applied to a cell cycle time series microarray dataset, SPD was not provided any prior knowledge of samples' time order or of which genes are cell-cycle regulated, yet SPD recovered the correct time order and identified many genes that have been associated with the cell cycle. When applied to B-cell differentiation data, SPD recovered the correct order of stages of normal B-cell differentiation and the linkage between preB-ALL tumor cells with their cell origin preB. When applied to mouse embryonic stem cell differentiation data, SPD uncovered a landscape of ESC differentiation into various lineages and genes that represent both generic and lineage specific processes. When applied to a prostate cancer microarray dataset, SPD identified gene modules that reflect a progression consistent with disease stages. SPD may be best viewed as a novel tool for synthesizing biological hypotheses because it provides a likely biological progression underlying a microarray dataset and, perhaps more importantly, the candidate genes that regulate that progression.
MeSH Terms
Algorithms
Animals
B-Lymphocytes/cytology
Cell Differentiation
Cell Lineage
Cluster Analysis
Computational Biology/methods
Databases, Factual
Embryonic Stem Cells/cytology
Gene Expression Profiling/methods
Humans
Male
Mice
Oligonucleotide Array Sequence Analysis/methods
Prostatic Neoplasms/metabolism
Software
Authors & Affiliations
3 authors, click to expand affiliations / ORCID
Qiu Peng
Department of Radiology, Stanford University, Stanford, California, USA.
[email protected]
Gentles Andrew J
Plevritis Sylvia K
References (27)
27 references, click to expand
-
Clustering gene expression data using a graph-theoretic approach: an application of minimum spanning trees.
Bioinformatics. 2002 Apr;18(4):536-45
PMID: 12016051
-
Efficient and accurate construction of genetic linkage maps from the minimum spanning tree of a graph.
PLoS Genet. 2008 Oct;4(10):e1000212
PMID: 18846212
-
Gene network interconnectedness and the generalized topological overlap measure.
BMC Bioinformatics. 2007 Jan 24;8:22
PMID: 17250769
-
Fast calculation of pairwise mutual information for gene regulatory network reconstruction.
Comput Methods Programs Biomed. 2009 May;94(2):177-80
PMID: 19167129
-
Ensemble dependence model for classification and prediction of cancer and normal gene expression data.
Bioinformatics. 2005 Jul 15;21(14):3114-21
PMID: 15879455
-
Significance analysis of time course microarray experiments.
Proc Natl Acad Sci U S A. 2005 Sep 6;102(36):12837-42
PMID: 16141318
-
Cluster analysis and display of genome-wide expression patterns.
Proc Natl Acad Sci U S A. 1998 Dec 8;95(25):14863-8
PMID: 9843981
-
Significance analysis of microarrays applied to the ionizing radiation response.
Proc Natl Acad Sci U S A. 2001 Apr 24;98(9):5116-21
PMID: 11309499
-
Transcription control of early B cell differentiation.
Curr Opin Immunol. 2010 Apr;22(2):161-7
PMID: 20144854
-
Tumor classification using phylogenetic methods on expression data.
J Theor Biol. 2004 Jun 21;228(4):477-96
PMID: 15178197
-
Gene expression profiles of prostate cancer reveal involvement of multiple molecular pathways in the metastatic process.
BMC Cancer. 2007 Apr 12;7:64
PMID: 17430594
-
Extracting dynamics from static cancer expression data.
IEEE/ACM Trans Comput Biol Bioinform. 2008 Apr-Jun;5(2):172-82
PMID: 18451427
-
Simultaneous class discovery and classification of microarray data using spectral analysis.
J Comput Biol. 2009 Jul;16(7):935-44
PMID: 19580522
-
ARACNE: an algorithm for the reconstruction of gene regulatory networks in a mammalian cellular context.
BMC Bioinformatics. 2006 Mar 20;7 Suppl 1:S7
PMID: 16723010
-
Neural network analysis of lymphoma microarray data: prognosis and diagnosis near-perfect.
BMC Bioinformatics. 2003 Apr 10;4:13
PMID: 12697066
-
Support vector machine classification and validation of cancer tissue samples using microarray expression data.
Bioinformatics. 2000 Oct;16(10):906-14
PMID: 11120680
-
Network-based inference of cancer progression from microarray data.
IEEE/ACM Trans Comput Biol Bioinform. 2009 Apr-Jun;6(2):200-12
PMID: 19407345
-
Analysis techniques for microarray time-series data.
J Comput Biol. 2002;9(2):317-30
PMID: 12015884
-
Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles.
Proc Natl Acad Sci U S A. 2005 Oct 25;102(43):15545-50
PMID: 16199517
-
Identification of genes periodically expressed in the human cell cycle and their expression in tumors.
Mol Biol Cell. 2002 Jun;13(6):1977-2000
PMID: 12058064
-
Network constrained clustering for gene microarray data.
Bioinformatics. 2005 Nov 1;21(21):4014-20
PMID: 16141248
-
Defining developmental potency and cell lineage trajectories by expression profiling of differentiating mouse embryonic stem cells.
DNA Res. 2009 Feb;16(1):73-80
PMID: 19112179
-
Characterization of early stages of human B cell development by gene expression profiling.
J Immunol. 2007 Sep 15;179(6):3662-71
PMID: 17785802
-
A general framework for weighted gene co-expression network analysis.
Stat Appl Genet Mol Biol. 2005;4:Article17
PMID: 16646834
-
Systematic determination of genetic network architecture.
Nat Genet. 1999 Jul;22(3):281-5
PMID: 10391217
-
Precedence Temporal Networks to represent temporal relationships in gene expression data.
J Biomed Inform. 2007 Dec;40(6):761-74
PMID: 17644489
-
Reconstructing the temporal ordering of biological samples using microarray data.
Bioinformatics. 2003 May 1;19(7):842-50
PMID: 12724294