GSEA-InContext: Identifying Novel and Common Patterns in Expression Experiments

Overview

Journal Bioinformatics

Publisher Oxford University Press

Specialty Biology

Date 2018 Jun 29

PMID 29950010

Citations 135

Authors

Rani K Powers

Andrew Goodspeed

Harrison Pielke-Lombardo

Aik-Choon Tan

James C Costello

Affiliations

Soon will be listed here.

Abstract

Motivation: Gene Set Enrichment Analysis (GSEA) is routinely used to analyze and interpret coordinate pathway-level changes in transcriptomics experiments. For an experiment where less than seven samples per condition are compared, GSEA employs a competitive null hypothesis to test significance. A gene set enrichment score is tested against a null distribution of enrichment scores generated from permuted gene sets, where genes are randomly selected from the input experiment. Looking across a variety of biological conditions, however, genes are not randomly distributed with many showing consistent patterns of up- or down-regulation. As a result, common patterns of positively and negatively enriched gene sets are observed across experiments. Placing a single experiment into the context of a relevant set of background experiments allows us to identify both the common and experiment-specific patterns of gene set enrichment.

Results: We compiled a compendium of 442 small molecule transcriptomic experiments and used GSEA to characterize common patterns of positively and negatively enriched gene sets. To identify experiment-specific gene set enrichment, we developed the GSEA-InContext method that accounts for gene expression patterns within a background set of experiments to identify statistically significantly enriched gene sets. We evaluated GSEA-InContext on experiments using small molecules with known targets to show that it successfully prioritizes gene sets that are specific to each experiment, thus providing valuable insights that complement standard GSEA analysis.

Availability And Implementation: GSEA-InContext implemented in Python, Supplementary results and the background expression compendium are available at: https://github.com/CostelloLab/GSEA-InContext.

Citing Articles

Integrative analysis of m7G methylation-associated genes prognostic signature with immunotherapy and identification of LARP1 as a key oncogene in head and neck squamous cell carcinoma.

Xu J, You Z, Zhu Z, Liu M, Zhang Z, Xu P Front Immunol. 2025; 16:1520070.

PMID: 40018039 PMC: 11864954. DOI: 10.3389/fimmu.2025.1520070.

Coix Seed Oil Alleviates DSS-Induced Ulcerative Colitis via Intestinal Barrier Repair and Ferroptosis Regulation.

Zeng Y, Li N, Deng B, Gu Y, Lu S, Liu Y J Inflamm Res. 2025; 18:2557-2581.

PMID: 39995826 PMC: 11849537. DOI: 10.2147/JIR.S501745.

Disulfidptosis-related gene signatures as prognostic biomarkers and predictors of immunotherapy response in HNSCC.

Qin H, Xu J, Yue Y, Chen M, Zhang Z, Xu P Front Immunol. 2025; 15:1456649.

PMID: 39896807 PMC: 11782277. DOI: 10.3389/fimmu.2024.1456649.

Screening of necroptosis-related genes and evaluating the prognostic capacity, clinical value, and the effect of their copy number variations in acute myeloid leukemia.

Wen D, Yan R, Zhang L, Zhang H, Chen X, Zhou J BMC Cancer. 2025; 25(1):71.

PMID: 39806277 PMC: 11727709. DOI: 10.1186/s12885-025-13439-y.

A comparative analysis of Marburg virus-infected bat and human models from public high-throughput sequencing data.

Xuan D, Yeh I, Liu H, Su C, Ko C, Ta H Int J Med Sci. 2025; 22(1):1-16.

PMID: 39744175 PMC: 11659840. DOI: 10.7150/ijms.100696.

References

Subramanian A, Tamayo P, Mootha V, Mukherjee S, Ebert B, Gillette M . Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005; 102(43):15545-50. PMC: 1239896. DOI: 10.1073/pnas.0506580102. View

Subramanian A, Narayan R, Corsello S, Peck D, Natoli T, Lu X . A Next Generation Connectivity Map: L1000 Platform and the First 1,000,000 Profiles. Cell. 2017; 171(6):1437-1452.e17. PMC: 5990023. DOI: 10.1016/j.cell.2017.10.049. View

Goeman J, Buhlmann P . Analyzing gene expression data in terms of gene sets: methodological issues. Bioinformatics. 2007; 23(8):980-7. DOI: 10.1093/bioinformatics/btm051. View

Shoemaker R . The NCI60 human tumour cell line anticancer drug screen. Nat Rev Cancer. 2006; 6(10):813-23. DOI: 10.1038/nrc1951. View

Chen L, Fink T, Ebbesen P, Zachar V . Temporal transcriptome of mouse ATDC5 chondroprogenitors differentiating under hypoxic conditions. Exp Cell Res. 2006; 312(10):1727-44. DOI: 10.1016/j.yexcr.2006.02.013. View

Liberzon A, Birger C, Thorvaldsdottir H, Ghandi M, Mesirov J, Tamayo P . The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Syst. 2016; 1(6):417-425. PMC: 4707969. DOI: 10.1016/j.cels.2015.12.004. View

Gustafsson M, Zheng X, Pereira T, Gradin K, Jin S, Lundkvist J . Hypoxia requires notch signaling to maintain the undifferentiated cell state. Dev Cell. 2005; 9(5):617-28. DOI: 10.1016/j.devcel.2005.09.010. View

Khatri P, Sirota M, Butte A . Ten years of pathway analysis: current approaches and outstanding challenges. PLoS Comput Biol. 2012; 8(2):e1002375. PMC: 3285573. DOI: 10.1371/journal.pcbi.1002375. View

Sahlgren C, Gustafsson M, Jin S, Poellinger L, Lendahl U . Notch signaling mediates hypoxia-induced tumor cell migration and invasion. Proc Natl Acad Sci U S A. 2008; 105(17):6392-7. PMC: 2359811. DOI: 10.1073/pnas.0802047105. View

10.

Baker M . Gene data to hit milestone. Nature. 2012; 487(7407):282-3. DOI: 10.1038/487282a. View

11.

Jiang Z, Gentleman R . Extensions to gene set enrichment. Bioinformatics. 2006; 23(3):306-13. DOI: 10.1093/bioinformatics/btl599. View

12.

Eklund A, Szallasi Z . Correction of technical bias in clinical microarray data improves concordance with known biological information. Genome Biol. 2008; 9(2):R26. PMC: 2374720. DOI: 10.1186/gb-2008-9-2-r26. View

13.

Moellering R, Cornejo M, Davis T, Del Bianco C, Aster J, Blacklow S . Direct inhibition of the NOTCH transcription factor complex. Nature. 2009; 462(7270):182-8. PMC: 2951323. DOI: 10.1038/nature08543. View

14.

Barbie D, Tamayo P, Boehm J, Kim S, Moody S, Dunn I . Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1. Nature. 2009; 462(7269):108-12. PMC: 2783335. DOI: 10.1038/nature08460. View

15.

Bolstad B, Irizarry R, Astrand M, Speed T . A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003; 19(2):185-93. DOI: 10.1093/bioinformatics/19.2.185. View

16.

Yoo M, Shin J, Kim J, Ryall K, Lee K, Lee S . DSigDB: drug signatures database for gene set analysis. Bioinformatics. 2015; 31(18):3069-71. PMC: 4668778. DOI: 10.1093/bioinformatics/btv313. View

17.

Davis S, Meltzer P . GEOquery: a bridge between the Gene Expression Omnibus (GEO) and BioConductor. Bioinformatics. 2007; 23(14):1846-7. DOI: 10.1093/bioinformatics/btm254. View

18.

Culhane A, Schroder M, Sultana R, Picard S, Martinelli E, Kelly C . GeneSigDB: a manually curated database and resource for analysis of gene expression signatures. Nucleic Acids Res. 2011; 40(Database issue):D1060-6. PMC: 3245038. DOI: 10.1093/nar/gkr901. View

19.

Arora V, Schenkein E, Murali R, Subudhi S, Wongvipat J, Balbas M . Glucocorticoid receptor confers resistance to antiandrogens by bypassing androgen receptor blockade. Cell. 2013; 155(6):1309-22. PMC: 3932525. DOI: 10.1016/j.cell.2013.11.012. View

20.

Mootha V, Lindgren C, Eriksson K, Subramanian A, Sihag S, Lehar J . PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003; 34(3):267-73. DOI: 10.1038/ng1180. View