{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,2,26]],"date-time":"2026-02-26T20:34:35Z","timestamp":1772138075083,"version":"3.50.1"},"reference-count":24,"publisher":"Oxford University Press (OUP)","issue":"6","license":[{"start":{"date-parts":[[2025,6,11]],"date-time":"2025-06-11T00:00:00Z","timestamp":1749600000000},"content-version":"vor","delay-in-days":10,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/100004440","name":"Wellcome","doi-asserted-by":"publisher","award":["206328\/Z\/17\/Z"],"award-info":[{"award-number":["206328\/Z\/17\/Z"]}],"id":[{"id":"10.13039\/100004440","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/100004440","name":"Wellcome","doi-asserted-by":"publisher","award":["203151\/Z\/16\/Z"],"award-info":[{"award-number":["203151\/Z\/16\/Z"]}],"id":[{"id":"10.13039\/100004440","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2025,6,2]]},"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:sec>\n                    <jats:title>Motivation<\/jats:title>\n                    <jats:p>A defining characteristic of all metazoan organisms is the existence of different cell states or cell types, driven by changes in gene expression kinetics, principally transcription, splicing and degradation rates. The RNA velocity framework utilizes both spliced and unspliced reads in single cell mRNA preparations to predict future cellular states and estimate transcriptional kinetics. However, current models assume either constant kinetic rates, rates equal for all genes, or rates completely independent of progression through differentiation. Consequently, current models for rate estimation are either underparametrized or overparametrized.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Results<\/jats:title>\n                    <jats:p>Here, we developed a new method (diffGEK) which overcomes this issue, and allows comparison of transcriptional rates across different biological conditions. diffGEK assumes that rates can vary over a trajectory, but are smooth functions of the differentiation process. Analysing Jak2 V617F mutant versus wild type mice for erythropoiesis, and Ezh2 KO versus wild type mice in myelopoiesis, revealed which genes show altered transcription, splicing or degradation rates between different conditions. Moreover, we observed that, for some genes, compensatory changes between different rates can result in comparable overall mRNA levels, thereby masking highly dynamic changes in gene expression kinetics in conventional expression analysis. Collectively, we report a robust pipeline for comparative expression analysis based on altered transcriptional kinetics to discover mechanistic differences missed by conventional approaches, with broad applicability across any biomedical research question where single cell expression data are available for both wild type and treatment\/mutant conditions.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Availability and implementation<\/jats:title>\n                    <jats:p>This study does not include new data. All the codes are available on github: https:\/\/github.com\/mebarile\/transcriptional_kinetics.<\/jats:p>\n                  <\/jats:sec>","DOI":"10.1093\/bioinformatics\/btaf316","type":"journal-article","created":{"date-parts":[[2025,5,21]],"date-time":"2025-05-21T09:10:19Z","timestamp":1747818619000},"source":"Crossref","is-referenced-by-count":0,"title":["diffGEK: differential gene expression kinetics"],"prefix":"10.1093","volume":"41","author":[{"ORCID":"https:\/\/orcid.org\/0009-0004-7740-473X","authenticated-orcid":false,"given":"Melania","family":"Barile","sequence":"first","affiliation":[{"name":"Department of Haematology, Wellcome\u2013MRC Cambridge Stem Cell Institute, University of Cambridge, Jeffrey Cheah Biomedical Centre , Cambridge CB2 0AW,","place":["United Kingdom"]},{"name":"Centre for Translational Stem Cell Biology, HKSTP , Hong Kong SAR,","place":["China"]}]},{"given":"Shirom","family":"Chabra","sequence":"additional","affiliation":[{"name":"Department of Haematology, Wellcome\u2013MRC Cambridge Stem Cell Institute, University of Cambridge, Jeffrey Cheah Biomedical Centre , Cambridge CB2 0AW,","place":["United Kingdom"]}]},{"given":"Tomoya","family":"Isobe","sequence":"additional","affiliation":[{"name":"Department of Haematology, Wellcome\u2013MRC Cambridge Stem Cell Institute, University of Cambridge, Jeffrey Cheah Biomedical Centre , Cambridge CB2 0AW,","place":["United Kingdom"]}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-6302-5705","authenticated-orcid":false,"given":"Berthold","family":"Gottgens","sequence":"additional","affiliation":[{"name":"Department of Haematology, Wellcome\u2013MRC Cambridge Stem Cell Institute, University of Cambridge, Jeffrey Cheah Biomedical Centre , Cambridge CB2 0AW,","place":["United Kingdom"]}]}],"member":"286","published-online":{"date-parts":[[2025,6,10]]},"reference":[{"key":"2025070408274002300_btaf316-B1","doi-asserted-by":"crossref","first-page":"197","DOI":"10.1186\/s13059-021-02414-y","article-title":"Coordinated changes in gene expression kinetics underlie both mouse and human erythroid maturation","volume":"22","author":"Barile","year":"2021","journal-title":"Genome Biol"},{"key":"2025070408274002300_btaf316-B2","doi-asserted-by":"crossref","first-page":"1408","DOI":"10.1038\/s41587-020-0591-3","article-title":"Generalizing RNA velocity to transient cell states through dynamical modeling","volume":"38","author":"Bergen","year":"2020","journal-title":"Nat Biotechnol"},{"key":"2025070408274002300_btaf316-B3","doi-asserted-by":"crossref","DOI":"10.3390\/jcm9092934","article-title":"Targeting PI3K\/Akt\/mTOR in AML: rationale and clinical evidence","volume":"9","author":"Darici","year":"2020","journal-title":"J Clin Med"},{"key":"2025070408274002300_btaf316-B4","doi-asserted-by":"crossref","first-page":"461","DOI":"10.1038\/s41587-019-0088-0","article-title":"Inferring population dynamics from single-cell RNA-sequencing time series data","volume":"37","author":"Fischer","year":"2019","journal-title":"Nat Biotechnol"},{"key":"2025070408274002300_btaf316-B5","doi-asserted-by":"crossref","first-page":"6586","DOI":"10.1038\/s41467-022-34188-7","article-title":"UniTVelo: temporally unified RNA velocity reinforces single-cell trajectory inference","volume":"13","author":"Gao","year":"2022","journal-title":"Nat Commun"},{"key":"2025070408274002300_btaf316-B6","doi-asserted-by":"crossref","first-page":"100426","DOI":"10.1016\/j.xgen.2023.100426","article-title":"Preleukemic single-cell landscapes reveal mutation-specific mechanisms and gene programs predictive of aml patient outcomes","volume":"3","author":"Isobe","year":"2023","journal-title":"Cell Genom"},{"key":"2025070408274002300_btaf316-B7","doi-asserted-by":"crossref","first-page":"2540","DOI":"10.1182\/blood-2013-07-517847","article-title":"The il-2\/cd25 axis maintains distinct subsets of chronic myeloid leukemia-initiating cells","volume":"123","author":"Kobayashi","year":"2014","journal-title":"Blood"},{"key":"2025070408274002300_btaf316-B8","doi-asserted-by":"crossref","first-page":"494","DOI":"10.1038\/s41586-018-0414-6","article-title":"RNA velocity of single cells","volume":"560","author":"La Manno","year":"2018","journal-title":"Nature"},{"key":"2025070408274002300_btaf316-B9","doi-asserted-by":"crossref","first-page":"159","DOI":"10.1038\/s41592-021-01346-6","article-title":"CellRank for directed single-cell fate mapping","volume":"19","author":"Lange","year":"2022","journal-title":"Nat Methods"},{"key":"2025070408274002300_btaf316-B10","doi-asserted-by":"crossref","first-page":"89","DOI":"10.1182\/blood.2019003439","article-title":"p53 activation during ribosome biogenesis regulates normal erythroid differentiation","volume":"137","author":"Le Goff","year":"2021","journal-title":"Blood"},{"key":"2025070408274002300_btaf316-B11","doi-asserted-by":"publisher","first-page":"99","DOI":"10.1038\/s41587-023-01728-5","article-title":"A relay velocity model infers cell-dependent RNA velocity","volume":"42","author":"Li","year":"2024","journal-title":"Nat Biotechnol"},{"key":"2025070408274002300_btaf316-B12","doi-asserted-by":"crossref","first-page":"6323","DOI":"10.1038\/s41467-021-26582-4","article-title":"Convergence of oncogenic cooperation at single-cell and single-gene levels drives leukemic transformation","volume":"12","author":"Liu","year":"2021","journal-title":"Nat Commun"},{"key":"2025070408274002300_btaf316-B13","doi-asserted-by":"crossref","DOI":"10.1126\/sciimmunol.aal2192","article-title":"Single-cell RNA-seq and computational analysis using temporal mixture modelling resolves Th1\/tfh fate bifurcation in malaria","volume":"2","author":"L\u00f6nnberg","year":"2017","journal-title":"Sci Immunol"},{"key":"2025070408274002300_btaf316-B14","doi-asserted-by":"crossref","first-page":"550","DOI":"10.1186\/s13059-014-0550-8","article-title":"Moderated estimation of fold change and dispersion for rna-seq data with deseq2","volume":"15","author":"Love","year":"2014","journal-title":"Genome Biol"},{"key":"2025070408274002300_btaf316-B15","doi-asserted-by":"crossref","first-page":"e2102669","DOI":"10.1002\/advs.202102669","article-title":"Fine-tuning of cholesterol homeostasis controls erythroid differentiation","volume":"9","author":"Lu","year":"2022","journal-title":"Adv Sci"},{"key":"2025070408274002300_btaf316-B16","doi-asserted-by":"crossref","first-page":"109113","DOI":"10.1016\/j.biopha.2019.109113","article-title":"Selective inhibition of aurora a and b kinases effectively induces cell cycle arrest in t(8; 21) acute myeloid leukemia","volume":"117","author":"Qi","year":"2019","journal-title":"Biomed Pharmacother"},{"key":"2025070408274002300_btaf316-B17","doi-asserted-by":"crossref","first-page":"e2105859118","DOI":"10.1073\/pnas.2105859118","article-title":"Representation learning of RNA velocity reveals robust cell transitions","volume":"118","author":"Qiao","year":"2021","journal-title":"Proc Natl Acad Sci USA"},{"key":"2025070408274002300_btaf316-B18","doi-asserted-by":"crossref","first-page":"979","DOI":"10.1038\/nmeth.4402","article-title":"Reversed graph embedding resolves complex single-cell trajectories","volume":"14","author":"Qiu","year":"2017","journal-title":"Nat Methods"},{"key":"2025070408274002300_btaf316-B19","doi-asserted-by":"crossref","first-page":"690","DOI":"10.1016\/j.cell.2021.12.045","article-title":"Mapping transcriptomic vector fields of single cells","volume":"185","author":"Qiu","year":"2022","journal-title":"Cell"},{"key":"2025070408274002300_btaf316-B20","doi-asserted-by":"crossref","first-page":"139","DOI":"10.1093\/bioinformatics\/btp616","article-title":"edgeR: a bioconductor package for differential expression analysis of digital gene expression data","volume":"26","author":"Robinson","year":"2010","journal-title":"Bioinformatics"},{"key":"2025070408274002300_btaf316-B21","doi-asserted-by":"crossref","first-page":"163","DOI":"10.1007\/s12185-022-03406-9","article-title":"The role of specialized cell cycles during erythroid lineage development: insights from single-cell rna sequencing","volume":"116","author":"Socolovsky","year":"2022","journal-title":"Int J Hematol"},{"key":"2025070408274002300_btaf316-B22","doi-asserted-by":"crossref","first-page":"1079","DOI":"10.1093\/biostatistics\/kxae017","article-title":"DifferentialRegulation: a bayesian hierarchical approach to identify differentially regulated genes","volume":"25","author":"Tiberi","year":"2024","journal-title":"Biostatistics"},{"key":"2025070408274002300_btaf316-B23","doi-asserted-by":"crossref","first-page":"1201","DOI":"10.1038\/s41467-020-14766-3","article-title":"Trajectory-based differential expression analysis for single-cell sequencing data","volume":"11","author":"Van den Berge","year":"2020","journal-title":"Nat Commun"},{"key":"2025070408274002300_btaf316-B24","doi-asserted-by":"crossref","first-page":"694","DOI":"10.1016\/j.cels.2024.07.002","article-title":"Global transcription regulation revealed from dynamical correlations in time-resolved single-cell RNA sequencing","volume":"15","author":"Volteras","year":"2024","journal-title":"Cell Syst"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/advance-article-pdf\/doi\/10.1093\/bioinformatics\/btaf316\/63465472\/btaf316.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/41\/6\/btaf316\/63465472\/btaf316.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/41\/6\/btaf316\/63465472\/btaf316.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,7,4]],"date-time":"2025-07-04T08:27:50Z","timestamp":1751617670000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/doi\/10.1093\/bioinformatics\/btaf316\/8160081"}},"subtitle":[],"editor":[{"given":"Laura","family":"Cantini","sequence":"additional","affiliation":[]}],"short-title":[],"issued":{"date-parts":[[2025,6]]},"references-count":24,"journal-issue":{"issue":"6","published-print":{"date-parts":[[2025,6,2]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btaf316","relation":{"has-preprint":[{"id-type":"doi","id":"10.1101\/2024.08.21.608952","asserted-by":"object"}]},"ISSN":["1367-4811"],"issn-type":[{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2025,6]]},"published":{"date-parts":[[2025,6]]},"article-number":"btaf316"}}