{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"institution":[{"name":"bioRxiv"}],"indexed":{"date-parts":[[2026,4,11]],"date-time":"2026-04-11T20:44:09Z","timestamp":1775940249573,"version":"3.50.1"},"posted":{"date-parts":[[2025,3,11]]},"group-title":"Bioinformatics","reference-count":36,"publisher":"openRxiv","license":[{"start":{"date-parts":[[2025,3,11]],"date-time":"2025-03-11T00:00:00Z","timestamp":1741651200000},"content-version":"vor","delay-in-days":0,"URL":"http:\/\/creativecommons.org\/licenses\/by-nc-nd\/4.0\/"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"accepted":{"date-parts":[[2025,3,11]]},"abstract":"<jats:title>Abstract<\/jats:title>\n                <jats:p>\n                  Differential expression analysis is crucial in genomics, yet existing methods primarily focus on detecting mean shifts. Variance shifts in gene expression are well-documented in studies of cellular signaling pathways, and more recently they have characterized aging, thus motivating the need for flexible detection approaches that include tests of expression variance changes. In this work, we present QRscore (Quantile Rank Score), a general method for detecting distributional shifts in gene expression by extending the Mann-Whitney test into a flexible family of rank-based tests. Here, we focus on implementing QRscore to detect shifts in mean and variance in gene expression, using weights designed from negative binomial (NB) and zero-inflated negative binomial (ZINB) models to combine the strengths of parametric and non-parametric approaches. We show through simulations that QRscore not only achieves high statistical power while controlling the false discovery rate (FDR), but also outperforms existing methods in detecting variance shifts and mean shifts. Applying QRscore to bulk RNA-seq data from the Genotype-Tissue Expression (GTEx) project, we identified numerous differentially dispersed genes and differentially expressed genes across 33 tissues. Notably, many genes have significant variance shifts but non-significant mean shifts. QRscore augments the genome bioinformatics toolkit by offering a powerful and flexible approach for differential expression analysis. QRscore is available in R, at\n                  <jats:ext-link xmlns:xlink=\"http:\/\/www.w3.org\/1999\/xlink\" ext-link-type=\"uri\" xlink:href=\"https:\/\/github.com\/songlab-cal\/QRscore\">https:\/\/github.com\/songlab-cal\/QRscore<\/jats:ext-link>\n                  .\n                <\/jats:p>","DOI":"10.1101\/2025.03.06.641952","type":"posted-content","created":{"date-parts":[[2025,3,11]],"date-time":"2025-03-11T11:04:48Z","timestamp":1741691088000},"source":"Crossref","is-referenced-by-count":0,"title":["Robust and Adaptive Non-Parametric Tests for Detecting General Distributional Shifts in Gene Expression"],"prefix":"10.64898","author":[{"given":"Fanding","family":"Zhou","sequence":"first","affiliation":[{"name":"Biostatistics Division, University of California, Berkeley"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-9455-7878","authenticated-orcid":false,"given":"Alan J.","family":"Aw","sequence":"additional","affiliation":[{"name":"Department of Statistics, University of California, Berkeley"},{"name":"Department of Genetics, University of Pennsylvania"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Dan D.","family":"Erdmann-Pham","sequence":"additional","affiliation":[{"name":"Department of Statistics, Stanford University"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-1165-9930","authenticated-orcid":false,"given":"Jonathan","family":"Fischer","sequence":"additional","affiliation":[{"name":"Department of Biostatistics, University of Florida"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-0734-9868","authenticated-orcid":false,"given":"Yun S.","family":"Song","sequence":"additional","affiliation":[{"name":"Department of Statistics, University of California, Berkeley"},{"name":"Computer Science Division, University of California, Berkeley"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"54368","reference":[{"key":"2026041113011192000_2025.03.06.641952v1.1","doi-asserted-by":"publisher","DOI":"10.1186\/1752-0509-4-1"},{"key":"2026041113011192000_2025.03.06.641952v1.2","doi-asserted-by":"publisher","DOI":"10.1371\/journal.pgen.1002207"},{"key":"2026041113011192000_2025.03.06.641952v1.3","doi-asserted-by":"publisher","DOI":"10.1371\/journal.pgen.1005428"},{"key":"2026041113011192000_2025.03.06.641952v1.4","doi-asserted-by":"publisher","DOI":"10.1152\/physiolgenomics.00128.2018"},{"key":"2026041113011192000_2025.03.06.641952v1.5","doi-asserted-by":"publisher","DOI":"10.1038\/nature08781"},{"key":"2026041113011192000_2025.03.06.641952v1.6","doi-asserted-by":"publisher","DOI":"10.1093\/hmg\/ddx424"},{"key":"2026041113011192000_2025.03.06.641952v1.7","doi-asserted-by":"publisher","DOI":"10.1016\/j.semcdb.2017.06.012"},{"key":"2026041113011192000_2025.03.06.641952v1.8","doi-asserted-by":"publisher","DOI":"10.1038\/s41586-020-2499-y"},{"key":"2026041113011192000_2025.03.06.641952v1.9","doi-asserted-by":"publisher","DOI":"10.1038\/nature04844"},{"key":"2026041113011192000_2025.03.06.641952v1.10","doi-asserted-by":"publisher","DOI":"10.1093\/bioinformatics\/btn142"},{"key":"2026041113011192000_2025.03.06.641952v1.11","doi-asserted-by":"publisher","DOI":"10.1186\/gb-2014-15-1-r1"},{"key":"2026041113011192000_2025.03.06.641952v1.12","doi-asserted-by":"publisher","DOI":"10.1111\/j.1467-9876.2005.00510.x"},{"key":"2026041113011192000_2025.03.06.641952v1.13","doi-asserted-by":"publisher","DOI":"10.1093\/nar\/gkx456"},{"key":"2026041113011192000_2025.03.06.641952v1.14","doi-asserted-by":"publisher","DOI":"10.1093\/nargab\/lqab124"},{"key":"2026041113011192000_2025.03.06.641952v1.15","doi-asserted-by":"publisher","DOI":"10.1186\/s13059-020-02238-2"},{"key":"2026041113011192000_2025.03.06.641952v1.16","doi-asserted-by":"publisher","DOI":"10.1186\/s13059-022-02648-4"},{"key":"2026041113011192000_2025.03.06.641952v1.17","doi-asserted-by":"publisher","DOI":"10.1186\/s12859-018-2201-7"},{"key":"2026041113011192000_2025.03.06.641952v1.18","doi-asserted-by":"publisher","DOI":"10.1371\/journal.pone.0232271"},{"key":"2026041113011192000_2025.03.06.641952v1.19","doi-asserted-by":"publisher","DOI":"10.1038\/ng.2653"},{"key":"2026041113011192000_2025.03.06.641952v1.20","doi-asserted-by":"crossref","unstructured":"Alexander M Mood . On the asymptotic efficiency of certain nonparametric two-sample tests. The Annals of Mathematical Statistics, pages 514\u2013522, 1954.","DOI":"10.1214\/aoms\/1177728719"},{"key":"2026041113011192000_2025.03.06.641952v1.21","doi-asserted-by":"publisher","DOI":"10.18632\/aging.103885"},{"key":"2026041113011192000_2025.03.06.641952v1.22","unstructured":"Sevgi D\u00f6nd\u00fc \u00d6zen and Serkan Kir . Ectodysplasin A2 receptor signaling in skeletal muscle pathophysiology. Trends in Molecular Medicine, 2024."},{"key":"2026041113011192000_2025.03.06.641952v1.23","doi-asserted-by":"publisher","DOI":"10.1089\/thy.2017.0359"},{"key":"2026041113011192000_2025.03.06.641952v1.24","doi-asserted-by":"publisher","DOI":"10.1038\/s41467-022-33509-0"},{"key":"2026041113011192000_2025.03.06.641952v1.25","doi-asserted-by":"publisher","DOI":"10.3390\/ijms25168794"},{"key":"2026041113011192000_2025.03.06.641952v1.26","doi-asserted-by":"publisher","DOI":"10.1111\/acel.13583"},{"key":"2026041113011192000_2025.03.06.641952v1.27","doi-asserted-by":"publisher","DOI":"10.3390\/ijms22126351"},{"key":"2026041113011192000_2025.03.06.641952v1.28","doi-asserted-by":"publisher","DOI":"10.1016\/j.febslet.2004.07.006"},{"key":"2026041113011192000_2025.03.06.641952v1.29","doi-asserted-by":"publisher","DOI":"10.1016\/j.tig.2013.08.003"},{"key":"2026041113011192000_2025.03.06.641952v1.30","doi-asserted-by":"publisher","DOI":"10.1186\/s13059-020-1935-5"},{"key":"2026041113011192000_2025.03.06.641952v1.31","doi-asserted-by":"publisher","DOI":"10.1002\/jcsm.12099"},{"key":"2026041113011192000_2025.03.06.641952v1.32","doi-asserted-by":"publisher","DOI":"10.1186\/1471-2105-14-1"},{"key":"2026041113011192000_2025.03.06.641952v1.33","doi-asserted-by":"publisher","DOI":"10.3390\/cells12081163"},{"key":"2026041113011192000_2025.03.06.641952v1.34","doi-asserted-by":"publisher","DOI":"10.1007\/s13277-015-3584-y"},{"key":"2026041113011192000_2025.03.06.641952v1.35","article-title":"Exact and efficient multivariate two-sample tests through generalized linear rank statistics","year":"2022","journal-title":"arXiv"},{"issue":"8","key":"2026041113011192000_2025.03.06.641952v1.36","first-page":"1","article-title":"Regression models for count data in R","volume":"27","year":"2008","journal-title":"Journal of Statistical Software"}],"container-title":[],"original-title":[],"link":[{"URL":"https:\/\/syndication.highwire.org\/content\/doi\/10.1101\/2025.03.06.641952","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,4,11]],"date-time":"2026-04-11T20:01:19Z","timestamp":1775937679000},"score":1,"resource":{"primary":{"URL":"http:\/\/biorxiv.org\/lookup\/doi\/10.1101\/2025.03.06.641952"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2025,3,11]]},"references-count":36,"URL":"https:\/\/doi.org\/10.1101\/2025.03.06.641952","relation":{},"subject":[],"published":{"date-parts":[[2025,3,11]]},"subtype":"preprint"}}