Abstract:Background: Gene set analysis methods have played a major role in generating biological interpretations from omics data such as gene expression datasets. However, most methods focus on detecting homogenous pattern changes in mean expression and methods detecting pattern changes in variance remain poorly explored. While a few studies attempted to use gene-level variance analysis, such approach remains under-utilized. When comparing two phenotypes, gene sets with distinct changes in subgroups under one phenotype are overlooked by available methods although they reflect meaningful biological differences between two phenotypes. Multivariate sample-level variance analysis methods are needed to detect such pattern changes. Results: We use ranking schemes based on minimum spanning tree to generalize the Cramer-Von Mises and Anderson-Darling univariate statistics into multivariate gene set analysis methods to detect differential sample variance or mean. We characterize these methods in addition to two methods developed earlier using simulation results with different parameters. We apply the developed methods to microarray gene expression dataset of prednisolone-resistant and prednisolone-sensitive children diagnosed with B-lineage acute lymphoblastic leukemia and bulk RNA-sequencing gene expression dataset of benign hyperplastic polyps and potentially malignant sessile serrated adenoma/polyps. One or both of the two compared phenotypes in each of these datasets have distinct molecular subtypes that contribute to heterogeneous differences. Our results show that methods designed to detect differential sample variance are able to detect specific hallmark signaling pathways associated with the two compared phenotypes as documented in available literature. Conclusions: The results in this study demonstrate the usefulness of methods designed to detect differential sample variance in providing biological interpretations when biologically relevant but heterogeneous changes between two phenotypes are prevalent in specific signaling pathways. Software implementation of the developed methods is available with detailed documentation from Bioconductor package GSAR. The available methods are applicable to gene expression datasets in a normalized matrix form and could be used with other omics datasets in a normalized matrix form with available collection of feature sets.

The Lorenz's side of Gene Set Enrichment Analysis: Similarity and divergences of the Gene Set Enrichment Analysis from the measure of wealth inequality

Gene Set linkage analysis: a tool for interpreting the overall functional impacts of observed transcriptomic changes

The Limitations of Simple Gene Set Enrichment Analysis Assuming Gene Independence

Analyzing Gene Expression Data in Terms of Gene Sets: Gene Set Enrichment Analysis

On testing the significance of sets of genes

Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles

Toward a gold standard for benchmarking gene set enrichment analysis

Improving GSEA for Analysis of Biologic Pathways for Differential Gene Expression Across a Binary Phenotype

Redundancy-aware unsupervised ranking based on game theory -- application to gene enrichment analysis

Identifying Differential Gene Sets using the Linear Combination of Genes with Maximum AUC

GAGE: generally applicable gene set enrichment for pathway analysis

Spectral gene set enrichment (SGSE)

Gene Set Enrichment Analysis (GSEA) for Interpreting Gene Expression Profiles

Toward Comprehensive Functional Analysis of Gene Lists Weighted by Gene Essentiality Scores.

Improving data interpretability with new differential sample variance gene set tests

GIGSEA: Genotype Imputed Gene Set Enrichment Analysis Using GWAS Summary Level Data

Moment based gene set tests

Investigating the Concordance of Gene Ontology Terms Reveals the Intra- and Inter-Platform Reproducibility of Enrichment Analysis

Assessment of Gene Set Enrichment Analysis using curated RNA-seq-based benchmarks

GSVA: gene set variation analysis for microarray and RNA-Seq data

Dual gene set enrichment analysis (dualGSEA); an R function that enables more robust biological discovery and pre-clinical model alignment from transcriptomics data