Differential Expression Analysis in RNA-seq Data Using a Geometric Approach

Abstract

Although differential gene expression (DGE) profiling in RNA-seq is used by many researchers, new packages and pipelines are continuously being presented as a result of an ongoing investigation. In this work, a geometric approach based on Supervised Variational Relevance Learning (Suvrel) was compared with DEpackages (edgeR, DESEq, baySeq, PoissonSeq, and limma) in the DGE profiling. The Suvrel method seeks to determine the relevance of characteristics (e.g., gene or transcript) based on intraclass and interclass distances. The comparison was performed using technical and biological replicates. For technical replicates, we used receiver operating characteristic (ROC) analysis, while for the other ones, we used robustness analysis. From ROC analysis, we found that geometric approach had a better performance than the DEpackages. Particularly, for a reduced list of differentially expressed genes (DEG), we noticed that this method had a remarkable advantage in ranking of most DEG (with a specificity ranging from 1 to 0.8). From robustness analysis associated to biological replicates, we found that geometric approach has comparable performance to the DEpackages. We conclude that the geometric approach had a slight overall better performance than the other methods. Moreover, it is a simple method that does not make any assumption about the distribution associated with RNA-seq data set. From this perspective, the relevance of this study was to show that a simple method can provide as good performance as more complex methods.