A Bayesian framework for ranking genes based on their statistical evidence for differential expression

Advances in sequencing technologies have revolutionised our ability to capture the complete RNA profile in tissue samples, allowing for comparative analyses of RNA levels between developmental processes, environmental responses, or treatments. However, quantifying changes in gene expression remains challenging, given inherent biological variability and technological limitations. To address this, we introduce a Bayesian framework for differential gene expression (DGE) analysis. Our framework unifies and streamlines a complex analysis, typically involving parameter estimations and multiple statistical tests, into a concise mathematical equation. This allows statistical evidence for differential expression to be computed rapidly and transparently. We show how this approach can be used to evaluate variabilty of individual genes between replicates. A comparison of our framework with existing tools revealed differences that can be explained by commonly employed thresh-olds in other packages. This motivated us to explore ranking genes based on their statistical evidence as opposed to a binary classification as DEGs. Our analysis leads us to advocate the use of Bayes factors within a rank-based approach. This framework offers enhanced computational efficiency and delivers a transparent way to analyse, interpret and communicate DGE results.