NDLI: A flexible count data model to fit the wide diversity of expression profiles arising from extensively replicated RNA-seq experiments

Content Provider	Springer Nature : BioMed Central
Author	Esnaola, Mikel Puig, Pedro Gonzalez, David Castelo, Robert Gonzalez, Juan R
Abstract	Background High-throughput RNA sequencing (RNA-seq) offers unprecedented power to capture the real dynamics of gene expression. Experimental designs with extensive biological replication present a unique opportunity to exploit this feature and distinguish expression profiles with higher resolution. RNA-seq data analysis methods so far have been mostly applied to data sets with few replicates and their default settings try to provide the best performance under this constraint. These methods are based on two well-known count data distributions: the Poisson and the negative binomial. The way to properly calibrate them with large RNA-seq data sets is not trivial for the non-expert bioinformatics user. Results Here we show that expression profiles produced by extensively-replicated RNA-seq experiments lead to a rich diversity of count data distributions beyond the Poisson and the negative binomial, such as Poisson-Inverse Gaussian or Pólya-Aeppli, which can be captured by a more general family of count data distributions called the Poisson-Tweedie. The flexibility of the Poisson-Tweedie family enables a direct fitting of emerging features of large expression profiles, such as heavy-tails or zero-inflation, without the need to alter a single configuration parameter. We provide a software package for R called tweeDEseq implementing a new test for differential expression based on the Poisson-Tweedie family. Using simulations on synthetic and real RNA-seq data we show that tweeDEseq yields P-values that are equally or more accurate than competing methods under different configuration parameters. By surveying the tiny fraction of sex-specific gene expression changes in human lymphoblastoid cell lines, we also show that tweeDEseq accurately detects differentially expressed genes in a real large RNA-seq data set with improved performance and reproducibility over the previously compared methodologies. Finally, we compared the results with those obtained from microarrays in order to check for reproducibility. Conclusions RNA-seq data with many replicates leads to a handful of count data distributions which can be accurately estimated with the statistical model illustrated in this paper. This method provides a better fit to the underlying biological variability; this may be critical when comparing groups of RNA-seq samples with markedly different count data distributions. The tweeDEseq package forms part of the Bioconductor project and it is available for download at http://www.bioconductor.org .
Related Links	https://bmcbioinformatics.biomedcentral.com/counter/pdf/10.1186/1471-2105-14-254.pdf
Ending Page	22
Page Count	22
Starting Page	1
File Format	HTM / HTML
ISSN	14712105
DOI	10.1186/1471-2105-14-254
Journal	BMC Bioinformatics
Issue Number	1
Volume Number	14
Language	English
Publisher	BioMed Central
Publisher Date	2013-08-21
Access Restriction	Open
Subject Keyword	Bioinformatics Microarrays Computational Biology Computer Appl. in Life Sciences Algorithms Differential Expression Gene Negative Binomial Distribution Differential Expression Analysis Negative Binomial Negative Binomial Model Computational Biology/Bioinformatics
Content Type	Text
Resource Type	Article
Subject	Molecular Biology Biochemistry Computer Science Applications Applied Mathematics Structural Biology
Journal Impact Factor	2.9/2023
5-Year Journal Impact Factor	3.6/2023

Sl.	Authority	Responsibilities	Communication Details
1	Ministry of Education (GoI), Department of Higher Education	Sanctioning Authority	https://www.education.gov.in/ict-initiatives
2	Indian Institute of Technology Kharagpur	Host Institute of the Project: The host institute of the project is responsible for providing infrastructure support and hosting the project	https://www.iitkgp.ac.in
3	National Digital Library of India Office, Indian Institute of Technology Kharagpur	The administrative and infrastructural headquarters of the project	Dr. B. Sutradhar bsutra@ndl.gov.in
4	Project PI / Joint PI	Principal Investigator and Joint Principal Investigators of the project	Dr. B. Sutradhar bsutra@ndl.gov.in Prof. Saswat Chakrabarti will be added soon
5	Website/Portal (Helpdesk)	Queries regarding NDLI and its services	support@ndl.gov.in
6	Contents and Copyright Issues	Queries related to content curation and copyright issues	content@ndl.gov.in
7	National Digital Library of India Club (NDLI Club)	Queries related to NDLI Club formation, support, user awareness program, seminar/symposium, collaboration, social media, promotion, and outreach	clubsupport@ndl.gov.in
8	Digital Preservation Centre (DPC)	Assistance with digitizing and archiving copyright-free printed books	dpc@ndl.gov.in
9	IDR Setup or Support	Queries related to establishment and support of Institutional Digital Repository (IDR) and IDR workshops	idr@ndl.gov.in

A flexible count data model to fit the wide diversity of expression profiles arising from extensively replicated RNA-seq experiments

GC-Content Normalization for RNA-Seq Data

Negative binomial additive model for RNA-Seq data analysis

NBLDA: negative binomial linear discriminant analysis for RNA-Seq data

QNB: differential RNA methylation analysis for count-based small-sample sequencing data with a quad-negative binomial model

Sample size calculation based on exact test for assessing differential expression analysis in RNA-seq data

Power analysis for RNA-Seq differential expression studies using generalized linear mixed effects models

Comparative studies of differential gene calling using RNA-Seq data

Data-based RNA-seq simulations by binomial thinning

A flexible count data model to fit the wide diversity of expression profiles arising from extensively replicated RNA-seq experiments

Similar Documents

A flexible count data model to fit the wide diversity of expression profiles arising from extensively replicated RNA-seq experiments

GC-Content Normalization for RNA-Seq Data

Negative binomial additive model for RNA-Seq data analysis

NBLDA: negative binomial linear discriminant analysis for RNA-Seq data

QNB: differential RNA methylation analysis for count-based small-sample sequencing data with a quad-negative binomial model

Sample size calculation based on exact test for assessing differential expression analysis in RNA-seq data

Power analysis for RNA-Seq differential expression studies using generalized linear mixed effects models

Comparative studies of differential gene calling using RNA-Seq data

Data-based RNA-seq simulations by binomial thinning

A flexible count data model to fit the wide diversity of expression profiles arising from extensively replicated RNA-seq experiments