Rank–rank hypergeometric overlap: identification of statistically significant overlap between gene-expression signatures

Comparing independent high-throughput gene-expression experiments can generate hypotheses about which gene-expression programs are shared between particular biological processes. Current techniques to compare expression profiles typically involve choosing a fixed differential expression threshold to summarize results, potentially reducing sensitivity to small but concordant changes. We present a threshold-free algorithm called Rank–rank Hypergeometric Overlap (RRHO). This algorithm steps through two gene lists ranked by the degree of differential expression observed in two profiling experiments, successively measuring the statistical significance of the number of overlapping genes. The output is a graphical map that shows the strength, pattern and bounds of correlation between two expression profiles. To demonstrate RRHO sensitivity and dynamic range, we identified shared expression networks in cancer microarray profiles driving tumor progression, stem cell properties and response to targeted kinase inhibition. We demonstrate how RRHO can be used to determine which model system or drug treatment best reflects a particular biological or disease response. The threshold-free and graphical aspects of RRHO complement other rank-based approaches such as Gene Set Enrichment Analysis (GSEA), for which RRHO is a 2D analog. Rank–rank overlap analysis is a sensitive, robust and web-accessible method for detecting and visualizing overlap trends between two complete, continuous gene-expression profiles. A web-based implementation of RRHO can be accessed at http://systems.crump.ucla.edu/rankrank/

Determining PTEN Functional Status by Network Component Deduced Transcription Factor Activities

PTEN-controlled PI3K-AKT-mTOR pathway represents one of the most deregulated signaling pathways in human cancers. With many small molecule inhibitors that target PI3K-AKT-mTOR pathway being exploited clinically, sensitive and reliable ways of stratifying patients according to their PTEN functional status and determining treatment outcomes are urgently needed. Heterogeneous loss of PTEN is commonly associated with human cancers and yet PTEN can also be regulated on epigenetic, transcriptional or post-translational levels, which makes the use of simple protein or gene expression-based analyses in determining PTEN status less accurate. In this study, we used network component analysis to identify 20 transcription factors (TFs) whose activities deduced from their target gene expressions were immediately altered upon the re-expression of PTEN in a PTEN-inducible system. Interestingly, PTEN controls the activities (TFA) rather than the expression levels of majority of these TFs and these PTEN-controlled TFAs are substantially altered in prostate cancer mouse models. Importantly, the activities of these TFs can be used to predict PTEN status in human prostate, breast and brain tumor samples with enhanced reliability when compared to straightforward IHC-based or expression-based analysis. Furthermore, our analysis indicates that unique sets of PTEN-controlled TFAs significantly contribute to specific tumor types. Together, our findings reveal that TFAs may be used as “signatures” for predicting PTEN functional status and elucidate the transcriptional architectures underlying human cancers caused by PTEN loss

Rank-rank hypergeometric overlap: identification of statistically significant overlap between gene-expression signatures.

Comparing independent high-throughput gene-expression experiments can generate hypotheses about which gene-expression programs are shared between particular biological processes. Current techniques to compare expression profiles typically involve choosing a fixed differential expression threshold to summarize results, potentially reducing sensitivity to small but concordant changes. We present a threshold-free algorithm called Rank-rank Hypergeometric Overlap (RRHO). This algorithm steps through two gene lists ranked by the degree of differential expression observed in two profiling experiments, successively measuring the statistical significance of the number of overlapping genes. The output is a graphical map that shows the strength, pattern and bounds of correlation between two expression profiles. To demonstrate RRHO sensitivity and dynamic range, we identified shared expression networks in cancer microarray profiles driving tumor progression, stem cell properties and response to targeted kinase inhibition. We demonstrate how RRHO can be used to determine which model system or drug treatment best reflects a particular biological or disease response. The threshold-free and graphical aspects of RRHO complement other rank-based approaches such as Gene Set Enrichment Analysis (GSEA), for which RRHO is a 2D analog. Rank-rank overlap analysis is a sensitive, robust and web-accessible method for detecting and visualizing overlap trends between two complete, continuous gene-expression profiles. A web-based implementation of RRHO can be accessed at http://systems.crump.ucla.edu/rankrank/

Cross-platform comparison of independent datasets identifies an immune signature associated with improved survival in metastatic melanoma

Platform and study differences in prognostic signatures from metastatic melanoma (MM) gene expression reports often hinder consensus arrival. We performed survival/outcome-based pairwise comparisons of three independent MM gene expression profiles using the threshold-free algorithm rank-rank hypergeometric overlap analysis (RRHO). We found statistically significant overlap for genes overexpressed in favorable outcome (FO) groups, but no overlap for poor outcome (PO) groups. This "favorable outcome signature" (FOS) of 228 genes coinciding on all three overlapping gene lists showed immune function predominated in FO MM. Surprisingly, specific cell signature-enrichment analysis showed B cell-associated genes enriched in FO MM, along with T cell-associated genes. Higher levels of B and T cells (p<0.05) and their relative proximity (p<0.05) were detected in FO-to-PO tumor comparisons from an independent MM patients cohort. Finally, expression of FOS in two independent Stage III MM tumor datasets correctly predicted clinical outcome in 12/14 and 44/70 patients using a weighted gene voting classifier (area under the curve values 0.96 and 0.75, respectively). This RRHO-based, cross-study analysis emphasizes the RRHO approach power, confirms T cells relevance for prolonged MM survival, supports a favorable role for B cells in anti-melanoma immunity, and suggests B cells potential as means of intervention in melanoma treatment.Fil: Lardone, Ricardo Dante. The John Wayne Cancer Institute; Estados Unidos. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Plaisier, Seema. The John Wayne Cancer Institute; Estados UnidosFil: Navarrete, Marian S.. The John Wayne Cancer Institute; Estados UnidosFil: Shamonki, Jaime M.. California Cryobank; Estados UnidosFil: Jalas, John R.. Providence Saint John’s Health Center; Estados UnidosFil: Sieling, Peter A. The John Wayne Cancer Institute; Estados UnidosFil: Lee, Delphine J.. The John Wayne Cancer Institute; Estados Unido

Determining PTEN functional status by network component deduced transcription factor activities.

PTEN-controlled PI3K-AKT-mTOR pathway represents one of the most deregulated signaling pathways in human cancers. With many small molecule inhibitors that target PI3K-AKT-mTOR pathway being exploited clinically, sensitive and reliable ways of stratifying patients according to their PTEN functional status and determining treatment outcomes are urgently needed. Heterogeneous loss of PTEN is commonly associated with human cancers and yet PTEN can also be regulated on epigenetic, transcriptional or post-translational levels, which makes the use of simple protein or gene expression-based analyses in determining PTEN status less accurate. In this study, we used network component analysis to identify 20 transcription factors (TFs) whose activities deduced from their target gene expressions were immediately altered upon the re-expression of PTEN in a PTEN-inducible system. Interestingly, PTEN controls the activities (TFA) rather than the expression levels of majority of these TFs and these PTEN-controlled TFAs are substantially altered in prostate cancer mouse models. Importantly, the activities of these TFs can be used to predict PTEN status in human prostate, breast and brain tumor samples with enhanced reliability when compared to straightforward IHC-based or expression-based analysis. Furthermore, our analysis indicates that unique sets of PTEN-controlled TFAs significantly contribute to specific tumor types. Together, our findings reveal that TFAs may be used as "signatures" for predicting PTEN functional status and elucidate the transcriptional architectures underlying human cancers caused by PTEN loss

Enhanced robustness of TFA-based signatures in predicting PTEN status in human cancer.

<p>(A) The Kaplan-Meier survival curves of patients with brain tumors stratified according to PTEN-controlled TFA and IHC analyses. (B) Log₁₀ transformed p-values of the χ² test evaluating the association of PTEN status with the hierarchical clustering-determined groups of human tumors. Clustering results are based on PTEN-controlled TFAs (red; <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0031053#pone-0031053-g003" target="_blank">Figure 3</a>), prostate cancer-related TFAs (blue), TFA-based (gold) and gene expression-based (green) signatures derived from PTEN IHC data in breast tumors. When three major clusters are observed in prostate and breast cancers, the χ² tests are performed to associate different PTEN status in group 1 and groups 2 plus 3.</p