{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2022,4,3]],"date-time":"2022-04-03T02:35:45Z","timestamp":1648953345999},"reference-count":21,"publisher":"Springer Science and Business Media LLC","issue":"S16","license":[{"start":{"date-parts":[[2019,12,1]],"date-time":"2019-12-01T00:00:00Z","timestamp":1575158400000},"content-version":"tdm","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"},{"start":{"date-parts":[[2019,12,2]],"date-time":"2019-12-02T00:00:00Z","timestamp":1575244800000},"content-version":"vor","delay-in-days":1,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"}],"content-domain":{"domain":["link.springer.com"],"crossmark-restriction":false},"short-container-title":["BMC Bioinformatics"],"published-print":{"date-parts":[[2019,12]]},"abstract":"<jats:title>Abstract<\/jats:title><jats:sec>\n<jats:title>Background<\/jats:title>\n<jats:p>A fundamental problem in RNA-seq data analysis is to identify genes or exons that are differentially expressed with varying experimental conditions based on the read counts. The relativeness of RNA-seq measurements makes the between-sample normalization of read counts an essential step in differential expression (DE) analysis. In most existing methods, the normalization step is performed prior to the DE analysis. Recently, Jiang and Zhan proposed a statistical method which introduces sample-specific normalization parameters into a joint model, which allows for simultaneous normalization and differential expression analysis from log-transformed RNA-seq data. Furthermore, an <jats:italic>\u2113<\/jats:italic><jats:sub>0<\/jats:sub> penalty is used to yield a sparse solution which selects a subset of DE genes. The experimental conditions are restricted to be categorical in their work.<\/jats:p>\n<\/jats:sec><jats:sec>\n<jats:title>Results<\/jats:title>\n<jats:p>In this paper, we generalize Jiang and Zhan\u2019s method to handle experimental conditions that are measured in continuous variables. As a result, genes with expression levels associated with a single or multiple covariates can be detected. As the problem being high-dimensional, non-differentiable and non-convex, we develop an efficient algorithm for model fitting.<\/jats:p>\n<\/jats:sec><jats:sec>\n<jats:title>Conclusions<\/jats:title>\n<jats:p>Experiments on synthetic data demonstrate that the proposed method outperforms existing methods in terms of detection accuracy when a large fraction of genes are differentially expressed in an asymmetric manner, and the performance gain becomes more substantial for larger sample sizes. We also apply our method to a real prostate cancer RNA-seq dataset to identify genes associated with pre-operative prostate-specific antigen (PSA) levels in patients.<\/jats:p>\n<\/jats:sec>","DOI":"10.1186\/s12859-019-3070-4","type":"journal-article","created":{"date-parts":[[2019,12,2]],"date-time":"2019-12-02T07:00:35Z","timestamp":1575270035000},"update-policy":"http:\/\/dx.doi.org\/10.1007\/springer_crossmark_policy","source":"Crossref","is-referenced-by-count":0,"title":["Joint between-sample normalization and differential expression detection through \u21130-regularized regression"],"prefix":"10.1186","volume":"20","author":[{"given":"Kefei","family":"Liu","sequence":"first","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Li","family":"Shen","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Hui","family":"Jiang","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"297","published-online":{"date-parts":[[2019,12,2]]},"reference":[{"key":"3070_CR1","doi-asserted-by":"publisher","first-page":"139","DOI":"10.1093\/bioinformatics\/btp616","volume":"26","author":"MD Robinson","year":"2010","unstructured":"Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010; 26:139\u201340.","journal-title":"Bioinformatics"},{"issue":"11","key":"3070_CR2","doi-asserted-by":"publisher","first-page":"91","DOI":"10.1093\/nar\/gku310","volume":"42","author":"X Zhou","year":"2014","unstructured":"Zhou X, Lindsay H, Robinson MD. Robustly detecting differential expression in RNA sequencing data using observation weights. Nucleic Acids Res. 2014; 42(11):91.","journal-title":"Nucleic Acids Res"},{"issue":"12","key":"3070_CR3","doi-asserted-by":"publisher","first-page":"550","DOI":"10.1186\/s13059-014-0550-8","volume":"15","author":"MI Love","year":"2014","unstructured":"Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014; 15(12):550.","journal-title":"Genome Biol"},{"issue":"7","key":"3070_CR4","doi-asserted-by":"publisher","first-page":"47","DOI":"10.1093\/nar\/gkv007","volume":"43","author":"ME Ritchie","year":"2015","unstructured":"Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015; 43(7):47.","journal-title":"Nucleic Acids Res"},{"issue":"2","key":"3070_CR5","doi-asserted-by":"publisher","first-page":"29","DOI":"10.1186\/gb-2014-15-2-r29","volume":"15","author":"CW Law","year":"2014","unstructured":"Law CW, Chen Y, Shi W, Smyth GK. Voom: precision weights unlock linear model analysis tools for RNA-seq read counts. Genome Biol. 2014; 15(2):29.","journal-title":"Genome Biol"},{"issue":"6","key":"3070_CR6","doi-asserted-by":"publisher","first-page":"671","DOI":"10.1093\/bib\/bbs046","volume":"14","author":"M-A Dillies","year":"2013","unstructured":"Dillies M-A, Rau A, Aubert J, Hennequet-Antier C, Jeanmougin M, Servant N, Keime C, Marot G, Castel D, Estelle J, Guernec G, Jagla B, Jouneau L, Laloe D, Le Gall C, Schaeffer BLffer, Le Crom S, Guedj M, Jaffrezic F, FSC. A comprehensive evaluation of normalization methods for illumina high-throughput RNA sequencing data analysis. Brief Bioinform. 2013; 14(6):671\u201383.","journal-title":"Brief Bioinform"},{"issue":"3","key":"3070_CR7","doi-asserted-by":"publisher","first-page":"25","DOI":"10.1186\/gb-2010-11-3-r25","volume":"11","author":"MD Robinson","year":"2010","unstructured":"Robinson MD, Oshlack A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 2010; 11(3):25.","journal-title":"Genome Biol"},{"key":"3070_CR8","doi-asserted-by":"publisher","first-page":"94","DOI":"10.1186\/1471-2105-11-94","volume":"11","author":"JH Bullard","year":"2010","unstructured":"Bullard JH, Purdom E, Hansen KD, Dudoit S. Evaluation of statistical methods for normalization and differential expression in mRNA-Seq experiments. BMC Bioinf. 2010; 11:94.","journal-title":"BMC Bioinf"},{"issue":"10","key":"3070_CR9","doi-asserted-by":"publisher","first-page":"106","DOI":"10.1186\/gb-2010-11-10-r106","volume":"11","author":"S Anders","year":"2010","unstructured":"Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010; 11(10):106.","journal-title":"Genome Biol"},{"issue":"1","key":"3070_CR10","doi-asserted-by":"publisher","first-page":"178","DOI":"10.1007\/s12561-016-9182-8","volume":"9","author":"H Jiang","year":"2017","unstructured":"Jiang H, Zhan T. Unit-free and robust detection of differential expression from rna-seq data. Stat Biosci. 2017; 9(1):178\u201399.","journal-title":"Stat Biosci"},{"issue":"2","key":"3070_CR11","doi-asserted-by":"publisher","first-page":"442","DOI":"10.1109\/TCBB.2018.2790918","volume":"16","author":"K Liu","year":"2019","unstructured":"Liu K, Ye J, Yang Y, Shen L, Jiang H. A Unified Model for Joint Normalization and Differential Gene Expression Detection in RNA-Seq data. IEEE\/ACM Trans Comput Biol Bioinforma. 2019; 16(2):442\u201354. https:\/\/doi.org\/10.1109\/tcbb.2018.2790918.","journal-title":"IEEE\/ACM Trans Comput Biol Bioinforma"},{"issue":"7","key":"3070_CR12","doi-asserted-by":"publisher","first-page":"621","DOI":"10.1038\/nmeth.1226","volume":"5","author":"A Mortazavi","year":"2008","unstructured":"Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq,. Nat Methods. 2008; 5(7):621\u20138.","journal-title":"Nat Methods"},{"issue":"5","key":"3070_CR13","doi-asserted-by":"publisher","first-page":"511","DOI":"10.1038\/nbt.1621","volume":"28","author":"C Trapnell","year":"2010","unstructured":"Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010; 28(5):511\u20135.","journal-title":"Nat Biotechnol"},{"issue":"4","key":"3070_CR14","doi-asserted-by":"publisher","first-page":"493","DOI":"10.1093\/bioinformatics\/btp692","volume":"26","author":"B Li","year":"2010","unstructured":"Li B, Ruotti V, Stewart RM, Thomson JA, Dewey CN. RNA-Seq gene expression estimation with read mapping uncertainty. Bioinformatics. 2010; 26(4):493\u2013500.","journal-title":"Bioinformatics"},{"issue":"1","key":"3070_CR15","doi-asserted-by":"publisher","first-page":"1","DOI":"10.1080\/10618600.2017.1390471","volume":"28","author":"T-Y Liu","year":"2019","unstructured":"Liu T-Y, Jiang H. Minimizing sum of truncated convex functions and its applications. J Comput Graph Stat. 2019; 28(1):1\u201310. (just-accepted). https:\/\/doi.org\/10.1080\/10618600.2017.1390471.","journal-title":"J Comput Graph Stat"},{"issue":"18","key":"3070_CR16","doi-asserted-by":"publisher","first-page":"3629","DOI":"10.1093\/bioinformatics\/bti593","volume":"21","author":"H Ji","year":"2005","unstructured":"Ji H, Wong WH. TileMap: create chromosomal map of tiling array hybridizations. Bioinformatics. 2005; 21(18):3629\u201336.","journal-title":"Bioinformatics"},{"key":"3070_CR17","doi-asserted-by":"publisher","first-page":"1011","DOI":"10.1016\/j.cell.2015.10.025","volume":"163","author":"CGAR Network","year":"2015","unstructured":"Network CGAR. The molecular taxonomy of primary prostate cancer. Cell. 2015; 163:1011\u201325. https:\/\/doi.org\/10.1016\/j.cell.2015.10.025.","journal-title":"Cell"},{"issue":"1","key":"3070_CR18","doi-asserted-by":"publisher","first-page":"18","DOI":"10.1186\/s12885-015-1025-3","volume":"15","author":"S Li","year":"2015","unstructured":"Li S, Zhu Y, Ma C, Qiu Z, Zhang X, Kang Z, Wu Z, Wang H, Xu X, Zhang H, et al.Downregulation of EphA5 by promoter methylation in human prostate cancer. BMC Cancer. 2015; 15(1):18.","journal-title":"BMC Cancer"},{"issue":"3","key":"3070_CR19","doi-asserted-by":"publisher","first-page":"290","DOI":"10.1053\/hupa.2002.32218","volume":"33","author":"C Busch","year":"2002","unstructured":"Busch C, Hanssen TA, Wagener C, \u00d6brink B. Down-regulation of CEACAM1 in human prostate cancer: correlation with loss of cell polarity, increased proliferation rate, and gleason grade 3 to 4 transition. Hum Pathol. 2002; 33(3):290\u20138.","journal-title":"Hum Pathol"},{"key":"3070_CR20","first-page":"4497","volume":"58","author":"RK Ross","year":"1998","unstructured":"Ross RK, Pike MC, Coetzee GA, Reichardt JK, Yu MC, Feigelson H, Stanczyk FZ, Kolonel LN, Henderson BE. Androgen metabolism and prostate cancer: establishing a model of genetic susceptibility. Cancer Res. 1998; 58:4497\u2013504.","journal-title":"Cancer Res"},{"issue":"5","key":"3070_CR21","doi-asserted-by":"publisher","first-page":"0155901","DOI":"10.1371\/journal.pone.0155901","volume":"11","author":"S Alinezhad","year":"2016","unstructured":"Alinezhad S, V\u00e4\u00e4n\u00e4nen R-M, Mattsson J, Li Y, Tallgr\u00e9n T, Ochoa NT, Bjartell A, \u00c5kerfelt M, Taimen P, Bostr\u00f6m PJ, et al.Validation of novel biomarkers for prostate cancer progression by the combination of bioinformatics, clinical and functional studies. PLoS ONE. 2016; 11(5):0155901.","journal-title":"PLoS ONE"}],"container-title":["BMC Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"http:\/\/link.springer.com\/content\/pdf\/10.1186\/s12859-019-3070-4.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"text-mining"},{"URL":"http:\/\/link.springer.com\/article\/10.1186\/s12859-019-3070-4\/fulltext.html","content-type":"text\/html","content-version":"vor","intended-application":"text-mining"},{"URL":"http:\/\/link.springer.com\/content\/pdf\/10.1186\/s12859-019-3070-4.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2020,11,30]],"date-time":"2020-11-30T19:16:05Z","timestamp":1606763765000},"score":1,"resource":{"primary":{"URL":"https:\/\/bmcbioinformatics.biomedcentral.com\/articles\/10.1186\/s12859-019-3070-4"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2019,12]]},"references-count":21,"journal-issue":{"issue":"S16","published-print":{"date-parts":[[2019,12]]}},"alternative-id":["3070"],"URL":"https:\/\/doi.org\/10.1186\/s12859-019-3070-4","relation":{},"ISSN":["1471-2105"],"issn-type":[{"value":"1471-2105","type":"electronic"}],"subject":[],"published":{"date-parts":[[2019,12]]},"assertion":[{"value":"2 December 2019","order":1,"name":"first_online","label":"First Online","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"Not applicable.","order":1,"name":"Ethics","group":{"name":"EthicsHeading","label":"Ethics approval and consent to participate"}},{"value":"Not applicable.","order":2,"name":"Ethics","group":{"name":"EthicsHeading","label":"Consent for publication"}},{"value":"The authors declare that they have no competing interests.","order":3,"name":"Ethics","group":{"name":"EthicsHeading","label":"Competing interests"}}],"article-number":"593"}}