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.

SAM thresholding and false discovery rates for detecting differential gene expression in DNA microarrays

A statistical method for the conservative adjustment of false discovery rate (q-value)

Considerations when using the significance analysis of microarrays (SAM) algorithm

A two-step strategy for detecting differential gene expression in cDNA microarray data

Comparison of Methods Used for Detecting Differentially Expressed Genes

Detecting Differentially Expressed Genes by Relative Entropy.

Identifying differentially expressed genes in human acute leukemia and mouse brain microarray datasets utilizing QTModel

A Statistical Framework for the Design of Microarray Experiments and Effective Detection of Differential Gene Expression

Degps is a Powerful Tool for Detecting Differential Expression in RNA-sequencing Studies

A comparison of methods for differential expression analysis of RNA-seq data

Exploratory data analysis for large-scale multiple testing problems and its application in gene expression studies

Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences

A comparative review of statistical methods for discovering differentially expressed genes in replicated microarray experiments

Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation

SAM-dPCR: Real-Time and High-throughput Absolute Quantification of Biological Samples Using Zero-Shot Segment Anything Model

False discovery rate control: Moving beyond the Benjamini–Hochberg method

Improving data interpretability with new differential sample variance gene set tests

Calculation of reliable transcript levels of annotated genes on the basis of multiple probe-sets in Affymetrix microarrays.

edgeR: a Bioconductor package for differential expression analysis of digital gene expression data

Comparative analysis of differential gene expression analysis tools for single-cell RNA sequencing data

Investigation of reproducibility of differentially expressed genes in DNA microarrays through statistical simulation