GLORB: Robust Bayesian inference for differential expression underglobal expression shifts
- 1Cancer Early Detection Advanced Research Center, Oregon Health and Science University
Abstract
Estimating differential gene expression is a common task in RNA-seq. Current methods mostly rely on normalization to reduce variance and increase accuracy. These methods are widely used and provide invaluable information about transcriptomic changes between biological conditions. However, widely used normalization methods are known to distort estimates of transcript differential expression when a majority of genes are upregulated or downregulated, or when the total RNA content per cell changes. Despite the presence of global expression shifts in a number of contexts, few methods exist that can provide accurate normalization and estimate linear models under this context without spike-in controls. Here, we present \textbf{GLORB} (\textbf{G}\textbf{L}\textbf{O}bal-shift \textbf{R}obust \textbf{B}ayesian model), a method for estimating generalized linear models under global upregulation. We develop two models that are able to recover differentially expressed genes and linear model coefficients with lower distortion of results. We show that our model's method of accounting for library size variance is consistent with DESeq2's median of ratios, and edgeR's trimmed mean of M-values under conditions when a minority of genes are upregulated and outperforms them under circumstances when most genes are either increased or decreased between groups. Finally, because our method does not rely on calculating geometric means for each gene it is able to work in datasets with much higher sparsity.
Related
- NicheCorr: Identifying Niche-Specific Differential Gene Co-Expression in Spatial Transcriptomics
- Differential Gene Expression and Functional Enrichment Analysis of Breast Tumor and Normal Tissue RNA-seq Data: A Reanalysis of GEO Dataset GSE210787
- Differential Gene Expression in the Tropical House Cricket and Its Iridovirus in Healthy versus Diseased Specimens
- Behavioral Mapping of Hierarchical Dynamical Models Accelerates Bayesian Parameter Inference
- abcneuralnet: Simulation-based Inference with Bayesian Neural Networks and Explainable AI for Population Genomics Inference
BibTeX
@misc{callahan2026glorb,
title = {{GLORB: Robust Bayesian inference for differential expression underglobal expression shifts}},
author = {R. L. Callahan and S. D. Coleman and T. T. M. Ngo},
howpublished = {bioRxiv},
year = {2026},
month = aug,
doi = {10.64898/2026.08.28.747928},
url = {https://www.biorxiv.org/content/10.64898/2026.08.28.747928v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - GLORB: Robust Bayesian inference for differential expression underglobal expression shifts
AU - R. L. Callahan
AU - S. D. Coleman
AU - T. T. M. Ngo
PY - 2026
DA - 2026/08/28
PB - bioRxiv
DO - 10.64898/2026.08.28.747928
UR - https://www.biorxiv.org/content/10.64898/2026.08.28.747928v1
AB - Estimating differential gene expression is a common task in RNA-seq. Current methods mostly rely on normalization to reduce variance and increase accuracy. These methods are widely used and provide invaluable information about transcriptomic changes between biological conditions. However, widely used normalization methods are known to distort estimates of transcript differential expression when a majority of genes are upregulated or downregulated, or when the total RNA content per cell changes. Despite the presence of global expression shifts in a number of contexts, few methods exist that can provide accurate normalization and estimate linear models under this context without spike-in controls. Here, we present \textbf{GLORB} (\textbf{G}\textbf{L}\textbf{O}bal-shift \textbf{R}obust \textbf{B}ayesian model), a method for estimating generalized linear models under global upregulation. We develop two models that are able to recover differentially expressed genes and linear model coefficients with lower distortion of results. We show that our model's method of accounting for library size variance is consistent with DESeq2's median of ratios, and edgeR's trimmed mean of M-values under conditions when a minority of genes are upregulated and outperforms them under circumstances when most genes are either increased or decreased between groups. Finally, because our method does not rely on calculating geometric means for each gene it is able to work in datasets with much higher sparsity.
ER -