{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,10,1]],"date-time":"2026-10-01T11:46:31Z","timestamp":1790855191047,"version":"4.1.0"},"reference-count":27,"publisher":"Oxford University Press (OUP)","issue":"4","license":[{"start":{"date-parts":[[2023,3,22]],"date-time":"2023-03-22T00:00:00Z","timestamp":1679443200000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2023,4,3]]},"abstract":"<jats:title>Abstract<\/jats:title><jats:sec><jats:title>Motivation<\/jats:title><jats:p>Finding outliers in RNA-sequencing (RNA-Seq) gene expression (GE) can help in identifying genes that are aberrant and cause Mendelian disorders. Recently developed models for this task rely on modeling RNA-Seq GE data using the negative binomial distribution (NBD). However, some of those models either rely on procedures for inferring NBD\u2019s parameters in a nonbiased way that are computationally demanding and thus make confounder control challenging, while others rely on less computationally demanding but biased procedures and convoluted confounder control approaches that hinder interpretability.<\/jats:p><\/jats:sec><jats:sec><jats:title>Results<\/jats:title><jats:p>In this article, we present OutSingle (Outlier detection using Singular Value Decomposition), an almost instantaneous way of detecting outliers in RNA-Seq GE data. It uses a simple log-normal approach for count modeling. For confounder control, it uses the recently discovered optimal hard threshold (OHT) method for noise detection, which itself is based on singular value decomposition (SVD). Due to its SVD\/OHT utilization, OutSingle\u2019s model is straightforward to understand and interpret. We then show that our novel method, when used on RNA-Seq GE data with real biological outliers masked by confounders, outcompetes the previous state-of-the-art model based on an ad hoc denoising autoencoder. Additionally, OutSingle can be used to inject artificial outliers masked by confounders, which is difficult to achieve with previous approaches. We describe a way of using OutSingle for outlier injection and proceed to show how OutSingle outperforms its competition on 16 out of 18 datasets that were generated from three real datasets using OutSingle\u2019s injection procedure with different outlier types and magnitudes. Our methods are applicable to other types of similar problems involving finding outliers in matrices under the presence of confounders.<\/jats:p><\/jats:sec><jats:sec><jats:title>Availability and implementation<\/jats:title><jats:p>The code for OutSingle is available at https:\/\/github.com\/esalkovic\/outsingle.<\/jats:p><\/jats:sec>","DOI":"10.1093\/bioinformatics\/btad142","type":"journal-article","created":{"date-parts":[[2023,3,22]],"date-time":"2023-03-22T08:21:57Z","timestamp":1679473317000},"source":"Crossref","is-referenced-by-count":13,"title":["OutSingle: a novel method of detecting and injecting outliers in RNA-Seq count data using the optimal hard threshold for singular values"],"prefix":"10.1093","volume":"39","author":[{"given":"Edin","family":"Salkovic","sequence":"first","affiliation":[{"name":"College of Science Engineering, Hamad Bin Khalifa University , Doha, Qatar"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Mohammad Amin","family":"Sadeghi","sequence":"additional","affiliation":[{"name":"Qatar Computing Research Institute, Hamad Bin Khalifa University , Doha, Qatar"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-8742-5519","authenticated-orcid":false,"given":"Abdelkader","family":"Baggag","sequence":"additional","affiliation":[{"name":"Qatar Computing Research Institute, Hamad Bin Khalifa University , Doha, Qatar"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Ahmed Gamal Rashed","family":"Salem","sequence":"additional","affiliation":[{"name":"Department of Computer Sciences, College of Engineering, Qatar University , P.O. Box: 2713 , Doha, Qatar"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-6700-5752","authenticated-orcid":false,"given":"Halima","family":"Bensmail","sequence":"additional","affiliation":[{"name":"Qatar Computing Research Institute, Hamad Bin Khalifa University , Doha, Qatar"}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"286","published-online":{"date-parts":[[2023,3,22]]},"reference":[{"key":"2023041120321980400_","doi-asserted-by":"crossref","first-page":"R106","DOI":"10.1186\/gb-2010-11-10-r106","article-title":"Differential expression analysis for sequence count data","volume":"11","author":"Anders","year":"2010","journal-title":"Genome Biol"},{"key":"2023041120321980400_","doi-asserted-by":"crossref","first-page":"448","DOI":"10.1016\/j.ajhg.2019.07.011","article-title":"Mendelian gene discovery: fast and furious with no end in sight","volume":"105","author":"Bamshad","year":"2019","journal-title":"Am J Hum Genet"},{"key":"2023041120321980400_","first-page":"1165","article-title":"The control of the false discovery rate in multiple testing under dependency","author":"Benjamini","year":"2001","journal-title":"Ann Stat"},{"key":"2023041120321980400_","doi-asserted-by":"crossref","first-page":"907","DOI":"10.1016\/j.ajhg.2018.10.025","article-title":"OUTRIDER: a statistical method for detecting aberrantly expressed genes in RNA sequencing data","volume":"103","author":"Brechtmann","year":"2018","journal-title":"Am J Hum Genet"},{"key":"2023041120321980400_","doi-asserted-by":"crossref","DOI":"10.1017\/9781108380690","volume-title":"Data-driven Science and Engineering: Machine Learning, Dynamical Systems, and Control","author":"Brunton","year":"2019"},{"issue":"386","key":"2023041120321980400_","doi-asserted-by":"crossref","DOI":"10.1126\/scitranslmed.aal5209","article-title":"Improving genetic diagnosis in Mendelian disease with transcriptome sequencing","volume":"9","author":"Cummings","year":"2017","journal-title":"Sci Transl Med"},{"key":"2023041120321980400_","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1186\/s12864-020-06938-8","article-title":"Bayesian gamma-negative binomial modeling of single-cell RNA sequencing data","volume":"21","author":"Dadaneh","year":"2020","journal-title":"BMC Genomics"},{"key":"2023041120321980400_","doi-asserted-by":"crossref","DOI":"10.1145\/1143844.1143874","volume-title":"The Relationship between Precision-Recall and ROC Curves","author":"Davis","year":"2006"},{"key":"2023041120321980400_","doi-asserted-by":"crossref","first-page":"5040","DOI":"10.1109\/TIT.2014.2323359","article-title":"The optimal hard threshold for singular values is","volume":"60","author":"Gavish","year":"2014","journal-title":"IEEE Trans Inf Theory"},{"key":"2023041120321980400_","doi-asserted-by":"crossref","first-page":"15824","DOI":"10.1038\/ncomms15824","article-title":"Genetic diagnosis of Mendelian disorders via RNA sequencing","volume":"8","author":"Kremer","year":"2017","journal-title":"Nat Commun"},{"key":"2023041120321980400_","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1186\/gb-2014-15-2-r29","article-title":"voom: precision weights unlock linear model analysis tools for RNA-seq read counts","volume":"15","author":"Law","year":"2014","journal-title":"Genome Biol"},{"key":"2023041120321980400_","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":"2023041120321980400_","doi-asserted-by":"crossref","first-page":"4288","DOI":"10.1093\/nar\/gks042","article-title":"Differential expression analysis of multifactor RNA-seq experiments with respect to biological variation","volume":"40","author":"McCarthy","year":"2012","journal-title":"Nucleic Acids Res"},{"key":"2023041120321980400_","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1038\/s41467-020-20573-7","article-title":"Detection of aberrant splicing events in RNA-seq data using FRASER","volume":"12","author":"Mertes","year":"2021","journal-title":"Nat Commun"},{"key":"2023041120321980400_","doi-asserted-by":"crossref","DOI":"10.1172\/JCI141500","article-title":"Transcriptome-directed analysis for mendelian disease diagnosis overcomes limitations of conventional genomic testing","volume":"131","author":"Murdock","year":"2021","journal-title":"J Clin Investig"},{"key":"2023041120321980400_","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":"2023041120321980400_","doi-asserted-by":"crossref","first-page":"75789","DOI":"10.1109\/ACCESS.2021.3082311","article-title":"A novel Bayesian outlier score based on the negative binomial distribution for detecting aberrantly expressed genes in RNA-Seq gene expression count data","volume":"9","author":"Salkovic","year":"2021","journal-title":"IEEE Access"},{"key":"2023041120321980400_","doi-asserted-by":"crossref","first-page":"101245","DOI":"10.1016\/j.jocs.2020.101245","article-title":"OutPyR: Bayesian inference for RNA-Seq outlier detection","volume":"47","author":"Salkovic","year":"2020","journal-title":"J Comput Sci"},{"key":"2023041120321980400_","doi-asserted-by":"crossref","first-page":"500","DOI":"10.1038\/nprot.2011.457","article-title":"Using probabilistic estimation of expression residuals (PEER) to obtain increased power and interpretability of gene expression analyses","volume":"7","author":"Stegle","year":"2012","journal-title":"Nat Protoc"},{"key":"2023041120321980400_","author":"Stewart","year":"1998"},{"key":"2023041120321980400_","doi-asserted-by":"crossref","first-page":"96","DOI":"10.1038\/s41586-020-2434-2","article-title":"Whole-genome sequencing of patients with rare diseases in a national health system","volume":"583","author":"Turro","year":"2020","journal-title":"Nature"},{"key":"2023041120321980400_","doi-asserted-by":"crossref","first-page":"261","DOI":"10.1038\/s41592-019-0686-2","article-title":"SciPy 1.0: fundamental algorithms for scientific computing in python","volume":"17","author":"Virtanen","year":"2020","journal-title":"Nat Methods"},{"key":"2023041120321980400_","doi-asserted-by":"crossref","first-page":"437","DOI":"10.1007\/s10545-015-9823-y","article-title":"Whole exome sequencing of suspected mitochondrial patients in clinical practice","volume":"38","author":"Wortmann","year":"2015","journal-title":"J Inherit Metab Dis"},{"key":"2023041120321980400_","doi-asserted-by":"crossref","first-page":"232","DOI":"10.1093\/biostatistics\/kxs033","article-title":"A new shrinkage estimator for dispersion improves differential expression detection in RNA-seq data","volume":"14","author":"Wu","year":"2013","journal-title":"Biostatistics"},{"key":"2023041120321980400_","doi-asserted-by":"crossref","first-page":"1276","DOI":"10.1038\/s41596-020-00462-5","article-title":"Detection of aberrant gene expression events in RNA sequencing data","volume":"16","author":"Y\u00e9pez","year":"2021","journal-title":"Nat Protoc"},{"key":"2023041120321980400_","doi-asserted-by":"crossref","first-page":"307","DOI":"10.1109\/TPAMI.2013.211","article-title":"Negative binomial process count and mixture modeling","volume":"37","author":"Zhou","year":"2015","journal-title":"IEEE Trans Pattern Anal Mach Intell"},{"key":"2023041120321980400_","doi-asserted-by":"crossref","first-page":"2084","DOI":"10.1093\/bioinformatics\/bty895","article-title":"Heavy-tailed prior distributions for sequence count data: removing the noise and preserving large differences","volume":"35","author":"Zhu","year":"2019","journal-title":"Bioinformatics"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/advance-article-pdf\/doi\/10.1093\/bioinformatics\/btad142\/49589872\/btad142.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/39\/4\/btad142\/49852162\/btad142.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/39\/4\/btad142\/49852162\/btad142.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,12,9]],"date-time":"2023-12-09T07:16:36Z","timestamp":1702106196000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/doi\/10.1093\/bioinformatics\/btad142\/7083276"}},"subtitle":[],"editor":[{"given":"Christina","family":"Kendziorski","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"editor"}]}],"short-title":[],"issued":{"date-parts":[[2023,3,22]]},"references-count":27,"journal-issue":{"issue":"4","published-print":{"date-parts":[[2023,4,3]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btad142","relation":{},"ISSN":["1367-4811"],"issn-type":[{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2023,4,1]]},"published":{"date-parts":[[2023,3,22]]},"article-number":"btad142"}}