{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,6,10]],"date-time":"2026-06-10T11:10:14Z","timestamp":1781089814337,"version":"3.54.1"},"reference-count":46,"publisher":"Oxford University Press (OUP)","issue":"9","content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2015,5,1]]},"abstract":"<jats:title>Abstract<\/jats:title><jats:p>Motivation: In recent years, gene expression studies have increasingly made use of high-throughput sequencing technology. In turn, research concerning the appropriate statistical methods for the analysis of digital gene expression (DGE) has flourished, primarily in the context of normalization and differential analysis.<\/jats:p><jats:p>Results: In this work, we focus on the question of clustering DGE profiles as a means to discover groups of co-expressed genes. We propose a Poisson mixture model using a rigorous framework for parameter estimation as well as the choice of the appropriate number of clusters. We illustrate co-expression analyses using our approach on two real RNA-seq datasets. A set of simulation studies also compares the performance of the proposed model with that of several related approaches developed to cluster RNA-seq or serial analysis of gene expression data.<\/jats:p><jats:p>Availability and and implementation: The proposed method is implemented in the open-source R package HTSCluster, available on CRAN.<\/jats:p><jats:p>Contact: \u00a0andrea.rau@jouy.inra.fr<\/jats:p><jats:p>Supplementary information: \u00a0Supplementary data are available at Bioinformatics online.<\/jats:p>","DOI":"10.1093\/bioinformatics\/btu845","type":"journal-article","created":{"date-parts":[[2015,1,7]],"date-time":"2015-01-07T03:18:32Z","timestamp":1420600712000},"page":"1420-1427","source":"Crossref","is-referenced-by-count":51,"title":["Co-expression analysis of high-throughput transcriptome sequencing data with Poisson mixture models"],"prefix":"10.1093","volume":"31","author":[{"given":"Andrea","family":"Rau","sequence":"first","affiliation":[{"name":"1 INRA, UMR1313 G\u00e9n\u00e9tique animale et biologie int\u00e9grative, Jouy-en-Josas, France, 2AgroParisTech, UMR1313 G\u00e9n\u00e9tique animale et biologie int\u00e9grative, Paris 05, France, 3Institut de Math\u00e9matiques de Toulouse, INSA de Toulouse, Universit\u00e9 de Toulouse, Toulouse, France, 4UMR AgroParisTech\/INRA MIA 518, Paris, France, 5INRA, UMR 1165 URGV, Saclay Plant Sciences, Evry, France, 6UEVE, UMR URGV, Saclay Plant Sciences, Evry, France, 7CNRS, ERL 8196, URGV, Saclay Plant Sciences, Evry, France and 8Inria Saclay - \u00cele-de-France, Orsay, France"},{"name":"1 INRA, UMR1313 G\u00e9n\u00e9tique animale et biologie int\u00e9grative, Jouy-en-Josas, France, 2AgroParisTech, UMR1313 G\u00e9n\u00e9tique animale et biologie int\u00e9grative, Paris 05, France, 3Institut de Math\u00e9matiques de Toulouse, INSA de Toulouse, Universit\u00e9 de Toulouse, Toulouse, France, 4UMR AgroParisTech\/INRA MIA 518, Paris, France, 5INRA, UMR 1165 URGV, Saclay Plant Sciences, Evry, France, 6UEVE, UMR URGV, Saclay Plant Sciences, Evry, France, 7CNRS, ERL 8196, URGV, Saclay Plant Sciences, Evry, France and 8Inria Saclay - \u00cele-de-France, Orsay, France"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Cathy","family":"Maugis-Rabusseau","sequence":"additional","affiliation":[{"name":"1 INRA, UMR1313 G\u00e9n\u00e9tique animale et biologie int\u00e9grative, Jouy-en-Josas, France, 2AgroParisTech, UMR1313 G\u00e9n\u00e9tique animale et biologie int\u00e9grative, Paris 05, France, 3Institut de Math\u00e9matiques de Toulouse, INSA de Toulouse, Universit\u00e9 de Toulouse, Toulouse, France, 4UMR AgroParisTech\/INRA MIA 518, Paris, France, 5INRA, UMR 1165 URGV, Saclay Plant Sciences, Evry, France, 6UEVE, UMR URGV, Saclay Plant Sciences, Evry, France, 7CNRS, ERL 8196, URGV, Saclay Plant Sciences, Evry, France and 8Inria Saclay - \u00cele-de-France, Orsay, France"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Marie-Laure","family":"Martin-Magniette","sequence":"additional","affiliation":[{"name":"1 INRA, UMR1313 G\u00e9n\u00e9tique animale et biologie int\u00e9grative, Jouy-en-Josas, France, 2AgroParisTech, UMR1313 G\u00e9n\u00e9tique animale et biologie int\u00e9grative, Paris 05, France, 3Institut de Math\u00e9matiques de Toulouse, INSA de Toulouse, Universit\u00e9 de Toulouse, Toulouse, France, 4UMR AgroParisTech\/INRA MIA 518, Paris, France, 5INRA, UMR 1165 URGV, Saclay Plant Sciences, Evry, France, 6UEVE, UMR URGV, Saclay Plant Sciences, Evry, France, 7CNRS, ERL 8196, URGV, Saclay Plant Sciences, Evry, France and 8Inria Saclay - \u00cele-de-France, Orsay, France"},{"name":"1 INRA, UMR1313 G\u00e9n\u00e9tique animale et biologie int\u00e9grative, Jouy-en-Josas, France, 2AgroParisTech, UMR1313 G\u00e9n\u00e9tique animale et biologie int\u00e9grative, Paris 05, France, 3Institut de Math\u00e9matiques de Toulouse, INSA de Toulouse, Universit\u00e9 de Toulouse, Toulouse, France, 4UMR AgroParisTech\/INRA MIA 518, Paris, France, 5INRA, UMR 1165 URGV, Saclay Plant Sciences, Evry, France, 6UEVE, UMR URGV, Saclay Plant Sciences, Evry, France, 7CNRS, ERL 8196, URGV, Saclay Plant Sciences, Evry, France and 8Inria Saclay - \u00cele-de-France, Orsay, France"},{"name":"1 INRA, UMR1313 G\u00e9n\u00e9tique animale et biologie int\u00e9grative, Jouy-en-Josas, France, 2AgroParisTech, UMR1313 G\u00e9n\u00e9tique animale et biologie int\u00e9grative, Paris 05, France, 3Institut de Math\u00e9matiques de Toulouse, INSA de Toulouse, Universit\u00e9 de Toulouse, Toulouse, France, 4UMR AgroParisTech\/INRA MIA 518, Paris, France, 5INRA, UMR 1165 URGV, Saclay Plant Sciences, Evry, France, 6UEVE, UMR URGV, Saclay Plant Sciences, Evry, France, 7CNRS, ERL 8196, URGV, Saclay Plant Sciences, Evry, France and 8Inria Saclay - \u00cele-de-France, Orsay, France"},{"name":"1 INRA, UMR1313 G\u00e9n\u00e9tique animale et biologie int\u00e9grative, Jouy-en-Josas, France, 2AgroParisTech, UMR1313 G\u00e9n\u00e9tique animale et biologie int\u00e9grative, Paris 05, France, 3Institut de Math\u00e9matiques de Toulouse, INSA de Toulouse, Universit\u00e9 de Toulouse, Toulouse, France, 4UMR AgroParisTech\/INRA MIA 518, Paris, France, 5INRA, UMR 1165 URGV, Saclay Plant Sciences, Evry, France, 6UEVE, UMR URGV, Saclay Plant Sciences, Evry, France, 7CNRS, ERL 8196, URGV, Saclay Plant Sciences, Evry, France and 8Inria Saclay - \u00cele-de-France, Orsay, France"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Gilles","family":"Celeux","sequence":"additional","affiliation":[{"name":"1 INRA, UMR1313 G\u00e9n\u00e9tique animale et biologie int\u00e9grative, Jouy-en-Josas, France, 2AgroParisTech, UMR1313 G\u00e9n\u00e9tique animale et biologie int\u00e9grative, Paris 05, France, 3Institut de Math\u00e9matiques de Toulouse, INSA de Toulouse, Universit\u00e9 de Toulouse, Toulouse, France, 4UMR AgroParisTech\/INRA MIA 518, Paris, France, 5INRA, UMR 1165 URGV, Saclay Plant Sciences, Evry, France, 6UEVE, UMR URGV, Saclay Plant Sciences, Evry, France, 7CNRS, ERL 8196, URGV, Saclay Plant Sciences, Evry, France and 8Inria Saclay - \u00cele-de-France, Orsay, France"}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"286","published-online":{"date-parts":[[2015,1,5]]},"reference":[{"key":"2023051308502889100_btu845-B1","doi-asserted-by":"crossref","first-page":"1","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":"2023051308502889100_btu845-B2","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1534\/genetics.110.114983","article-title":"Statistical design and analysis of RNA-Seq data","volume":"185","author":"Auer","year":"2010","journal-title":"Genetics"},{"key":"2023051308502889100_btu845-B3","doi-asserted-by":"crossref","first-page":"57","DOI":"10.1093\/bfgp\/elr041","article-title":"Differential expression\u2013the next generation and beyond","volume":"11","author":"Auer","year":"2012","journal-title":"Brief. Funct. Genomics"},{"key":"2023051308502889100_btu845-B4","doi-asserted-by":"crossref","first-page":"455","DOI":"10.1007\/s11222-011-9236-1","article-title":"Slope heuristics: overview and implementation","volume":"22","author":"Baudry","year":"2012","journal-title":"Stat. Comp."},{"key":"2023051308502889100_btu845-B5","doi-asserted-by":"crossref","first-page":"561","DOI":"10.1016\/S0167-9473(02)00163-9","article-title":"Choosing starting values for the EM algorithm for getting the highest likelihood in multivariate Gaussian mixture models","volume":"41","author":"Biernacki","year":"2003","journal-title":"Comp. Stat. Data Anal."},{"key":"2023051308502889100_btu845-B6","doi-asserted-by":"crossref","first-page":"203","DOI":"10.1007\/s100970100031","article-title":"Gaussian model selection","volume":"3","author":"Birg\u00e9","year":"2001","journal-title":"J. Eur. Math. Soc."},{"key":"2023051308502889100_btu845-B7","doi-asserted-by":"crossref","first-page":"33","DOI":"10.1007\/s00440-006-0011-8","article-title":"Minimal penalties for Gaussian model selection","volume":"138","author":"Birg\u00e9","year":"2006","journal-title":"Probab. Thoery Relat. Fields"},{"key":"2023051308502889100_btu845-B8","doi-asserted-by":"crossref","first-page":"R51","DOI":"10.1186\/gb-2004-5-7-r51","article-title":"Clustering analysis of SAGE data using a Poisson approach","volume":"5","author":"Cai","year":"2004","journal-title":"Genome Biol."},{"key":"2023051308502889100_btu845-B9","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1080\/03610927408827101","article-title":"A dendrite method for cluster analysis","volume":"3","author":"Cali\u0144ski","year":"1974","journal-title":"Commun. Stat. Theory Methods"},{"key":"2023051308502889100_btu845-B10","doi-asserted-by":"crossref","first-page":"315","DOI":"10.1016\/0167-9473(92)90042-E","article-title":"A classification EM algorithm for clustering and two stochastic versions","volume":"14","author":"Celeux","year":"1992","journal-title":"Comp. Stat. Data Anal."},{"key":"2023051308502889100_btu845-B11","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1111\/j.2517-6161.1977.tb01600.x","article-title":"Maximum likelihood from incomplete data via the EM algorithm","volume":"39","author":"Dempster","year":"1977","journal-title":"J. Royal Stat. Soc. Series B (Methodological)"},{"key":"2023051308502889100_btu845-B12","doi-asserted-by":"crossref","first-page":"671","DOI":"10.1093\/bib\/bbs046","article-title":"A comprehensive evaluation of normalization methods for Illumina high-throughput RNA sequencing data analysis","volume":"14","author":"Dillies","year":"2013","journal-title":"Brief. Bioinform."},{"key":"2023051308502889100_btu845-B13","doi-asserted-by":"crossref","first-page":"14863","DOI":"10.1073\/pnas.95.25.14863","article-title":"Cluster analysis and display of genome-wide expression patterns","volume":"95","author":"Eisen","year":"1998","journal-title":"PNAS"},{"key":"2023051308502889100_btu845-B14","doi-asserted-by":"crossref","DOI":"10.1186\/1471-2105-12-449","article-title":"ReCount: a multi-experiment resource of analysis-ready RNA-seq gene count datasets","volume":"12","author":"Frazee","year":"2011","journal-title":"BMC Bioinformatics"},{"key":"2023051308502889100_btu845-B15","doi-asserted-by":"crossref","first-page":"473","DOI":"10.1038\/nature09715","article-title":"The development transcriptome of Drosophila melanogaster","volume":"471","author":"Graveley","year":"2011","journal-title":"Nature"},{"key":"2023051308502889100_btu845-B16","doi-asserted-by":"crossref","first-page":"193","DOI":"10.1007\/BF01908075","article-title":"Comparing partitions","volume":"2","author":"Hubert","year":"1985","journal-title":"J. Classif."},{"key":"2023051308502889100_btu845-B17","doi-asserted-by":"crossref","first-page":"1370","DOI":"10.1109\/TKDE.2004.68","article-title":"Cluster analysis for gene expression data: a survey","volume":"16","author":"Jiang","year":"2004","journal-title":"IEEE Trans. Knowl. Data Eng."},{"key":"2023051308502889100_btu845-B18","doi-asserted-by":"crossref","first-page":"63","DOI":"10.1080\/0266476022000018510","article-title":"An EM algorithm for multivariate Poisson distribution and related models","volume":"30","author":"Karlis","year":"2003","journal-title":"J. Appl. Stat."},{"issue":"ISMB","key":"2023051308502889100_btu845-B19","doi-asserted-by":"crossref","first-page":"i383","DOI":"10.1093\/bioinformatics\/btr247","article-title":"Characterization and improvement of RNA-Seq precision in quantitative transcript expression profiling","volume":"27","author":"\u0141abaj","year":"2011","journal-title":"Bioinformatics"},{"key":"2023051308502889100_btu845-B20","doi-asserted-by":"crossref","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":"2023051308502889100_btu845-B21","doi-asserted-by":"crossref","first-page":"523","DOI":"10.1093\/biostatistics\/kxr031","article-title":"Normalization, testing, and false discovery rate estimation for RNA-sequencing data","volume":"13","author":"Li","year":"2012","journal-title":"Biostatistics"},{"key":"2023051308502889100_btu845-B22","first-page":"281","article-title":"Some methods for classification and analysis of multivariate observations","volume-title":"Proceedings of the 5th Berkeley Symposium on Mathematical Statistics and Probability","author":"MacQueen","year":"1967"},{"key":"2023051308502889100_btu845-B23","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":"2023051308502889100_btu845-B24","doi-asserted-by":"crossref","DOI":"10.4135\/9781412984713","volume-title":"Latent Class Analysis","author":"McCutcheon","year":"1987"},{"key":"2023051308502889100_btu845-B25","doi-asserted-by":"crossref","DOI":"10.1186\/1471-2164-12-293","article-title":"RNA-seq: technical variability and sampling","volume":"12","author":"McIntyre","year":"2011","journal-title":"BMC Genomics"},{"key":"2023051308502889100_btu845-B26","doi-asserted-by":"crossref","DOI":"10.1002\/047172842X","volume-title":"Analyzing Microarray Gene Expression Data","author":"McLachlan","year":"2004"},{"key":"2023051308502889100_btu845-B27","doi-asserted-by":"crossref","DOI":"10.1002\/0471721182","volume-title":"Finite Mixture Models","author":"McLachlan","year":"2000"},{"issue":"28","key":"2023051308502889100_btu845-B28","article-title":"A systematic comparison and evaluation of high density exon arrays and RNA-seq technology used to unravel the peripheral blood transcriptome of sickle cell disease","volume":"5","author":"Naghavachari","year":"2012","journal-title":"BMC Medical Genomics"},{"key":"2023051308502889100_btu845-B29","doi-asserted-by":"crossref","DOI":"10.1186\/1745-6150-4-14","article-title":"Transcript length bias in RNA-seq data confounds systems biology","volume":"4","author":"Oshlack","year":"2009","journal-title":"Biol. Direct"},{"key":"2023051308502889100_btu845-B30","article-title":"On the estimation of mixtures of Poisson regression models with large numbers of components","author":"Papastamoulis","year":"2014","journal-title":"Comp. Stat. Data Anal."},{"key":"2023051308502889100_btu845-B31","doi-asserted-by":"crossref","first-page":"2146","DOI":"10.1093\/bioinformatics\/btt350","article-title":"Data-based filtering for replicated high-throughput transcriptome sequencing experiments","volume":"29","author":"Rau","year":"2013","journal-title":"Bioinformatics"},{"key":"2023051308502889100_btu845-B32","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"},{"issue":"R25","key":"2023051308502889100_btu845-B33","article-title":"A scaling normalization method for differential expression analysis of RNA-seq data","volume":"11","author":"Robinson","year":"2010","journal-title":"Genome Biol."},{"key":"2023051308502889100_btu845-B34","doi-asserted-by":"crossref","first-page":"2881","DOI":"10.1093\/bioinformatics\/btm453","article-title":"Moderated statistical tests for assessing differences in tag abundance","volume":"23","author":"Robinson","year":"2007","journal-title":"Bioinformatics"},{"key":"2023051308502889100_btu845-B35","doi-asserted-by":"crossref","first-page":"461","DOI":"10.1214\/aos\/1176344136","article-title":"Estimating the dimension of a model","volume":"6","author":"Schwarz","year":"1978","journal-title":"Ann. Stat."},{"key":"2023051308502889100_btu845-B36","doi-asserted-by":"crossref","first-page":"903","DOI":"10.1038\/nbt.2957","article-title":"A comprehensive assessment of RNA-seq accuracy, reproducibility and information content by the Sequencing Quality Control Consortium","volume":"32","author":"SEQC\/MAQC-III Consortium","year":"2014","journal-title":"Nat. Biotechnol."},{"issue":"160","key":"2023051308502889100_btu845-B37","article-title":"RNA-Seq Atlas of Glycine max: a guide to the soybean transcriptome","volume":"10","author":"Severin","year":"2010","journal-title":"BMC Plant Biol."},{"key":"2023051308502889100_btu845-B38","doi-asserted-by":"crossref","first-page":"197","DOI":"10.1093\/bioinformatics\/btt632","article-title":"Model-based clustering for RNA-seq data","volume":"30","author":"Si","year":"2014","journal-title":"Bioinformatics"},{"issue":"91","key":"2023051308502889100_btu845-B39","article-title":"A comparison of methods for differential expression analysis of RNA-seq data","volume":"14","author":"Soneson","year":"2013","journal-title":"BMC Bioinformatics"},{"key":"2023051308502889100_btu845-B40","doi-asserted-by":"crossref","first-page":"956","DOI":"10.1126\/science.1160342","article-title":"A global view of gene activity and alternative splicing by deep sequencing of the human transcriptome","volume":"15","author":"Sultan","year":"2008","journal-title":"Science"},{"key":"2023051308502889100_btu845-B41","doi-asserted-by":"crossref","first-page":"511","DOI":"10.1038\/nbt.1621","article-title":"Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation","volume":"28","author":"Trapnell","year":"2010","journal-title":"Nat. Biotechnol."},{"key":"2023051308502889100_btu845-B42","doi-asserted-by":"crossref","first-page":"926","DOI":"10.1038\/nbt.3001","article-title":"The concordance between RNA-seq and microarray data depends on chemical treatment and transcript abundance","volume":"32","author":"Wang","year":"2014","journal-title":"Nat. Biotechnol."},{"key":"2023051308502889100_btu845-B43","doi-asserted-by":"crossref","first-page":"236","DOI":"10.1080\/01621459.1963.10500845","article-title":"Hierarchical grouping to optimize an objective function","volume":"58","author":"Ward","year":"1963","journal-title":"J. Am. Stat. Assoc."},{"key":"2023051308502889100_btu845-B44","doi-asserted-by":"crossref","first-page":"2493","DOI":"10.1214\/11-AOAS493","article-title":"Classification and clustering of sequencing data using a Poisson model","volume":"5","author":"Witten","year":"2011","journal-title":"Ann. Appl. Stat."},{"key":"2023051308502889100_btu845-B45","doi-asserted-by":"crossref","first-page":"977","DOI":"10.1093\/bioinformatics\/17.10.977","article-title":"Model-based clustering and data transformations for gene expression data","volume":"17","author":"Yeung","year":"2001","journal-title":"Bioinformatics"},{"key":"2023051308502889100_btu845-B46","doi-asserted-by":"crossref","first-page":"e91","DOI":"10.1093\/nar\/gku310","article-title":"Robustly detecting differential expression in RNA sequencing data using observation weights","volume":"42","author":"Zhou","year":"2014","journal-title":"Nucleic Acids Res."}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/31\/9\/1420\/50306125\/bioinformatics_31_9_1420.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/31\/9\/1420\/50306125\/bioinformatics_31_9_1420.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2024,6,6]],"date-time":"2024-06-06T11:01:30Z","timestamp":1717671690000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/31\/9\/1420\/200490"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2015,1,5]]},"references-count":46,"journal-issue":{"issue":"9","published-print":{"date-parts":[[2015,5,1]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btu845","relation":{},"ISSN":["1367-4811","1367-4803"],"issn-type":[{"value":"1367-4811","type":"electronic"},{"value":"1367-4803","type":"print"}],"subject":[],"published-other":{"date-parts":[[2015,5,1]]},"published":{"date-parts":[[2015,1,5]]}}}