{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,2,27]],"date-time":"2026-02-27T15:04:04Z","timestamp":1772204644714,"version":"3.50.1"},"reference-count":37,"publisher":"Oxford University Press (OUP)","issue":"13","license":[{"start":{"date-parts":[[2018,6,27]],"date-time":"2018-06-27T00:00:00Z","timestamp":1530057600000},"content-version":"vor","delay-in-days":0,"URL":"http:\/\/creativecommons.org\/licenses\/by-nc\/4.0\/"}],"funder":[{"DOI":"10.13039\/100000001","name":"NSF","doi-asserted-by":"publisher","award":["CCF-1553281"],"award-info":[{"award-number":["CCF-1553281"]}],"id":[{"id":"10.13039\/100000001","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/100000199","name":"USDA","doi-asserted-by":"publisher","id":[{"id":"10.13039\/100000199","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/100005825","name":"NIFA","doi-asserted-by":"publisher","award":["06-505570-01006"],"award-info":[{"award-number":["06-505570-01006"]}],"id":[{"id":"10.13039\/100005825","id-type":"DOI","asserted-by":"publisher"}]},{"name":"TEES-AgriLife Center for Bioinformatics and Genomic Systems Engineering"},{"name":"CBGSE"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2018,7,1]]},"abstract":"<jats:title>Abstract<\/jats:title><jats:sec><jats:title>Motivation<\/jats:title><jats:p>High-throughput sequencing technologies, in particular RNA sequencing (RNA-seq), have become the basic practice for genomic studies in biomedical research. In addition to studying genes individually, for example, through differential expression analysis, investigating co-ordinated expression variations of genes may help reveal the underlying cellular mechanisms to derive better understanding and more effective prognosis and intervention strategies. Although there exists a variety of co-expression network based methods to analyze microarray data for this purpose, instead of blindly extending these methods for microarray data that may introduce unnecessary bias, it is crucial to develop methods well adapted to RNA-seq data to identify the functional modules of genes with similar expression patterns.<\/jats:p><\/jats:sec><jats:sec><jats:title>Results<\/jats:title><jats:p>We have developed a fully Bayesian covariate-dependent negative binomial factor analysis (dNBFA) method\u2014dNBFA\u2014for RNA-seq count data, to capture coordinated gene expression changes, while considering effects from covariates reflecting different influencing factors. Unlike existing co-expression network based methods, our proposed model does not require multiple ad-hoc choices on data processing, transformation, as well as co-expression measures and can be directly applied to RNA-seq data. Furthermore, being capable of incorporating covariate information, the proposed method can tackle setups with complex confounding factors in different experiment designs. Finally, the natural model parameterization removes the need for a normalization preprocessing step, as commonly adopted to compensate for the effect of sequencing-depth variations. Efficient Bayesian inference of model parameters is derived by exploiting conditional conjugacy via novel data augmentation techniques. Experimental results on several real-world RNA-seq datasets on complex diseases suggest dNBFA as a powerful tool for discovering the gene modules with significant differential expression and meaningful biological insight.<\/jats:p><\/jats:sec><jats:sec><jats:title>Availability and implementation<\/jats:title><jats:p>dNBFA is implemented in R language and is available at https:\/\/github.com\/siamakz\/dNBFA.<\/jats:p><\/jats:sec>","DOI":"10.1093\/bioinformatics\/bty237","type":"journal-article","created":{"date-parts":[[2018,4,12]],"date-time":"2018-04-12T19:32:51Z","timestamp":1523561571000},"page":"i61-i69","source":"Crossref","is-referenced-by-count":6,"title":["Covariate-dependent negative binomial factor analysis of RNA sequencing data"],"prefix":"10.1093","volume":"34","author":[{"given":"Siamak","family":"Zamani Dadaneh","sequence":"first","affiliation":[{"name":"Department of Electrical and Computer Engineering, TEES-AgriLife Center for Bioinformatics and Genomic Systems Engineering, Texas A&M University, College Station, TX, USA"}]},{"given":"Mingyuan","family":"Zhou","sequence":"additional","affiliation":[{"name":"Department of Information, Risk, and Operations Management, The University of Texas at Austin, Austin, TX, USA"}]},{"given":"Xiaoning","family":"Qian","sequence":"additional","affiliation":[{"name":"Department of Electrical and Computer Engineering, TEES-AgriLife Center for Bioinformatics and Genomic Systems Engineering, Texas A&M University, College Station, TX, USA"}]}],"member":"286","published-online":{"date-parts":[[2018,6,27]]},"reference":[{"key":"2023051604341170200_bty237-B1","doi-asserted-by":"crossref","first-page":"34.","DOI":"10.1186\/1743-7075-8-34","article-title":"Nutritional and metabolic status of children with autism vs. neurotypical children, and the association with autism severity","volume":"8","author":"Adams","year":"2011","journal-title":"Nutr. Metab"},{"key":"2023051604341170200_bty237-B2","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":"2023051604341170200_bty237-B3","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1189\/jlb.1205707","article-title":"The immune response in autism: a new frontier for autism research","volume":"80","author":"Ashwood","year":"2006","journal-title":"J. Leukocyte Biol"},{"key":"2023051604341170200_bty237-B4","doi-asserted-by":"crossref","first-page":"e500.","DOI":"10.1038\/tp.2014.136","article-title":"Association between SNAP-25 gene polymorphisms and cognition in autism: functional consequences and potential therapeutic strategies","volume":"5","author":"Braida","year":"2015","journal-title":"Trans. Psychiatry"},{"key":"2023051604341170200_bty237-B5","doi-asserted-by":"crossref","first-page":"1061","DOI":"10.1038\/nature07385","article-title":"Comprehensive genomic characterization defines human glioblastoma genes and core pathways","volume":"455","author":"Cancer Genome Atlas (TCGA) Research Network and Others","year":"2008","journal-title":"Nature"},{"key":"2023051604341170200_bty237-B6","doi-asserted-by":"crossref","first-page":"327","DOI":"10.1080\/00031305.1995.10476177","article-title":"Understanding the Metropolis-Hastings algorithm","volume":"49","author":"Chib","year":"1995","journal-title":"Am. Stat"},{"key":"2023051604341170200_bty237-B7","doi-asserted-by":"crossref","first-page":"109.","DOI":"10.1186\/1471-2105-10-109","article-title":"Identifying set-wise differential co-expression in gene expression microarray data","volume":"10","author":"Cho","year":"2009","journal-title":"BMC Bioinformatics"},{"key":"2023051604341170200_bty237-B8","doi-asserted-by":"crossref","first-page":"4348","DOI":"10.1093\/bioinformatics\/bti722","article-title":"Differential coexpression analysis using microarray data and its application to human cancer","volume":"21","author":"Choi","year":"2005","journal-title":"Bioinformatics"},{"key":"2023051604341170200_bty237-B9","doi-asserted-by":"crossref","first-page":"2780","DOI":"10.1093\/bioinformatics\/btp502","article-title":"Statistical methods for gene set co-expression analysis","volume":"25","author":"Choi","year":"2009","journal-title":"Bioinformatics"},{"key":"2023051604341170200_bty237-B10","article-title":"BNP-Seq: Bayesian nonparametric differential expression analysis of sequencing count data","author":"Dadaneh","year":"2017","journal-title":"J. Am. Stat. Assoc"},{"key":"2023051604341170200_bty237-B11","doi-asserted-by":"crossref","DOI":"10.1007\/978-3-319-07212-8","volume-title":"Statistical Analysis of Next Generation Sequencing Data","author":"Datta","year":"2014"},{"key":"2023051604341170200_bty237-B12","doi-asserted-by":"crossref","first-page":"7835","DOI":"10.1073\/pnas.0602460103","article-title":"Neurons produce type I interferon during viral encephalitis","volume":"103","author":"Delhaye","year":"2006","journal-title":"Proc. Natl. Acad. Sci"},{"key":"2023051604341170200_bty237-B13","doi-asserted-by":"crossref","first-page":"5748.","DOI":"10.1038\/ncomms6748","article-title":"Transcriptome analysis reveals dysregulation of innate immune response genes and neuronal activity-dependent genes in autism","volume":"5","author":"Gupta","year":"2014","journal-title":"Nat. Commun"},{"key":"2023051604341170200_bty237-B14","doi-asserted-by":"crossref","first-page":"4074.","DOI":"10.1038\/ncomms5074","article-title":"The impact of the metabotropic glutamate receptor and other gene family interaction networks on autism","volume":"5","author":"Hadley","year":"2014","journal-title":"Nat. Commun"},{"key":"2023051604341170200_bty237-B15","doi-asserted-by":"crossref","first-page":"275","DOI":"10.1002\/aur.42","article-title":"Reduced levels of immunoglobulin in children with autism correlates with behavioral symptoms","volume":"1","author":"Heuer","year":"2008","journal-title":"Autism Res"},{"key":"2023051604341170200_bty237-B16","doi-asserted-by":"crossref","first-page":"e1000117.","DOI":"10.1371\/journal.pcbi.1000117","article-title":"Geometric interpretation of gene coexpression network analysis","volume":"4","author":"Horvath","year":"2008","journal-title":"PLoS Comput. Biol"},{"key":"2023051604341170200_bty237-B17","doi-asserted-by":"crossref","DOI":"10.1002\/0471715816","volume-title":"Univariate Discrete Distributions","author":"Johnson","year":"2005"},{"key":"2023051604341170200_bty237-B18","doi-asserted-by":"crossref","first-page":"559.","DOI":"10.1186\/1471-2105-9-559","article-title":"WGCNA: an R package for weighted correlation network analysis","volume":"9","author":"Langfelder","year":"2008","journal-title":"BMC Bioinformatics"},{"key":"2023051604341170200_bty237-B19","doi-asserted-by":"crossref","first-page":"R29.","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":"2023051604341170200_bty237-B20","doi-asserted-by":"crossref","first-page":"2699","DOI":"10.1093\/bioinformatics\/btx214","article-title":"Network module-based model in the differential expression analysis for RNA-seq","volume":"33","author":"Lei","year":"2017","journal-title":"Bioinformatics"},{"key":"2023051604341170200_bty237-B21","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":"2023051604341170200_bty237-B22","doi-asserted-by":"crossref","first-page":"189","DOI":"10.1093\/bib\/bbn001","article-title":"Gene-set approach for expression pattern analysis","volume":"9","author":"Nam","year":"2008","journal-title":"Brief. Bioinformatics"},{"key":"2023051604341170200_bty237-B23","doi-asserted-by":"crossref","first-page":"389","DOI":"10.1016\/j.molmed.2011.03.001","article-title":"Maternal infection and immune involvement in autism","volume":"17","author":"Patterson","year":"2011","journal-title":"Trends Mol. Med"},{"key":"2023051604341170200_bty237-B24","doi-asserted-by":"crossref","first-page":"1339","DOI":"10.1080\/01621459.2013.829001","article-title":"Bayesian inference for logistic models using P\u00f3lya\u2013Gamma latent variables","volume":"108","author":"Polson","year":"2013","journal-title":"J. Am. Statist. Assoc"},{"key":"2023051604341170200_bty237-B25","doi-asserted-by":"crossref","first-page":"1551","DOI":"10.1126\/science.1073374","article-title":"Hierarchical organization of modularity in metabolic networks","volume":"297","author":"Ravasz","year":"2002","journal-title":"Science"},{"key":"2023051604341170200_bty237-B26","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":"2023051604341170200_bty237-B27","doi-asserted-by":"crossref","first-page":"497.","DOI":"10.1186\/1471-2105-11-497","article-title":"DiffCoEx: a simple and sensitive method to find differentially coexpressed gene modules","volume":"11","author":"Tesson","year":"2010","journal-title":"BMC Bioinformatics"},{"key":"2023051604341170200_bty237-B28","doi-asserted-by":"crossref","first-page":"311","DOI":"10.1016\/S0198-8859(02)00374-9","article-title":"The transmission disequilibrium test suggests that HLA-DR4 and DR13 are linked to autism spectrum disorder","volume":"63","author":"Torres","year":"2002","journal-title":"Human Immunol"},{"key":"2023051604341170200_bty237-B29","doi-asserted-by":"crossref","first-page":"1.","DOI":"10.1155\/2012\/959073","article-title":"HLA immune function genes in autism","volume":"2012","author":"Torres","year":"2012","journal-title":"Autism Res. Treatment"},{"key":"2023051604341170200_bty237-B30","doi-asserted-by":"crossref","first-page":"396","DOI":"10.1016\/j.cell.2008.10.015","article-title":"Autism and brain development","volume":"135","author":"Walsh","year":"2008","journal-title":"Cell"},{"key":"2023051604341170200_bty237-B31","doi-asserted-by":"crossref","first-page":"952","DOI":"10.1093\/bioinformatics\/btv677","article-title":"TCGA2STAT: simple TCGA data access for integrated statistical analysis in R","volume":"32","author":"Wan","year":"2015","journal-title":"Bioinformatics"},{"key":"2023051604341170200_bty237-B32","doi-asserted-by":"crossref","first-page":"57","DOI":"10.1038\/nrg2484","article-title":"RNA-Seq: a revolutionary tool for transcriptomics","volume":"10","author":"Wang","year":"2009","journal-title":"Nat. Rev. Genet"},{"key":"2023051604341170200_bty237-B33","doi-asserted-by":"crossref","first-page":"97","DOI":"10.1016\/0165-5728(96)00052-5","article-title":"Strong association of the third hypervariable region of HLA-DR\u03b21 with autism","volume":"67","author":"Warren","year":"1996","journal-title":"J. Neuroimmunol"},{"key":"2023051604341170200_bty237-B34","author":"Zhou","year":"2017"},{"key":"2023051604341170200_bty237-B35","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. Machine Intel"},{"key":"2023051604341170200_bty237-B36","author":"Zhou","year":"2012"},{"key":"2023051604341170200_bty237-B37","doi-asserted-by":"crossref","first-page":"14595","DOI":"10.1523\/JNEUROSCI.2257-10.2010","article-title":"Changes in prefrontal axons may disrupt the network in autism","volume":"30","author":"Zikopoulos","year":"2010","journal-title":"J. Neurosci"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/34\/13\/i61\/50316387\/bioinformatics_34_13_i61.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/34\/13\/i61\/50316387\/bioinformatics_34_13_i61.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2024,7,5]],"date-time":"2024-07-05T23:33:26Z","timestamp":1720222406000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/34\/13\/i61\/5045747"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2018,6,27]]},"references-count":37,"journal-issue":{"issue":"13","published-print":{"date-parts":[[2018,7,1]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/bty237","relation":{},"ISSN":["1367-4803","1367-4811"],"issn-type":[{"value":"1367-4803","type":"print"},{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2018,7,1]]},"published":{"date-parts":[[2018,6,27]]}}}