» Articles » PMID: 39698741

To Tweak or Not to Tweak. How Exploiting Flexibilities in Gene Set Analysis Leads to Overoptimism

Overview
Journal Biom J
Specialty Public Health
Date 2024 Dec 19
PMID 39698741
Authors
Affiliations
Soon will be listed here.
Abstract

Gene set analysis, a popular approach for analyzing high-throughput gene expression data, aims to identify sets of genes that show enriched expression patterns between two conditions. In addition to the multitude of methods available for this task, users are typically left with many options when creating the required input and specifying the internal parameters of the chosen method. This flexibility can lead to uncertainty about the "right" choice, further reinforced by a lack of evidence-based guidance. Especially when their statistical experience is scarce, this uncertainty might entice users to produce preferable results using a "trial-and-error" approach. While it may seem unproblematic at first glance, this practice can be viewed as a form of "cherry-picking" and cause an optimistic bias, rendering the results nonreplicable on independent data. After this problem has attracted a lot of attention in the context of classical hypothesis testing, we now aim to raise awareness of such overoptimism in the different and more complex context of gene set analyses. We mimic a hypothetical researcher who systematically selects the analysis variants yielding their preferred results, thereby considering three distinct goals they might pursue. Using a selection of popular gene set analysis methods, we tweak the results in this way for two frequently used benchmark gene expression data sets. Our study indicates that the potential for overoptimism is particularly high for a group of methods frequently used despite being commonly criticized. We conclude by providing practical recommendations to counter overoptimism in research findings in gene set analysis and beyond.

References
1.
Love M, Huber W, Anders S . Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014; 15(12):550. PMC: 4302049. DOI: 10.1186/s13059-014-0550-8. View

2.
Reimand J, Isserlin R, Voisin V, Kucera M, Tannus-Lopes C, Rostamianfar A . Pathway enrichment analysis and visualization of omics data using g:Profiler, GSEA, Cytoscape and EnrichmentMap. Nat Protoc. 2019; 14(2):482-517. PMC: 6607905. DOI: 10.1038/s41596-018-0103-9. View

3.
Huang D, Sherman B, Lempicki R . Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2008; 37(1):1-13. PMC: 2615629. DOI: 10.1093/nar/gkn923. View

4.
Ballouz S, Pavlidis P, Gillis J . Using predictive specificity to determine when gene set analysis is biologically meaningful. Nucleic Acids Res. 2017; 45(4):e20. PMC: 5389513. DOI: 10.1093/nar/gkw957. View

5.
Maleki F, Ovens K, Hogan D, Kusalik A . Gene Set Analysis: Challenges, Opportunities, and Future Research. Front Genet. 2020; 11:654. PMC: 7339292. DOI: 10.3389/fgene.2020.00654. View