How well do RNA-Seq differential gene expression tools perform in a complex eukaryote? A case study in Arabidopsis thaliana

Kimon Froussios, Nicholas Schurch, Katarzyna Mackinnon, Marek Gierlinski, Celine Duc, Gordon G. Simpson (Lead / Corresponding author), Geoffrey J. Barton (Lead / Corresponding author)

Research output: Contribution to journalArticle

94 Downloads (Pure)

Abstract

RNA-seq experiments are usually carried out in three or fewer replicates. In order to work well with so few samples, differential gene expression (DGE) tools typically assume the form of the underlying gene expression distribution. In this paper, the statistical properties of gene expression from RNA-seq are investigated in the complex eukaryote, Arabidopsis thaliana, extending and generalizing the results of previous work in the simple eukaryote Saccharomyces cerevisiae. Results: We show that, consistent with the results in S.cerevisiae, more gene expression measurements in A.thaliana are consistent with being drawn from an underlying negative binomial distribution than either a log-normal distribution or a normal distribution, and that the size and complexity of the A.thaliana transcriptome does not influence the false positive rate performance of nine widely used DGE tools tested here. We therefore recommend the use of DGE tools that are based on the negative binomial distribution. Availability and implementation: The raw data for the 17 WT Arabidopsis thaliana datasets is available from the European Nucleotide Archive (E-MTAB-5446). The processed and aligned data can be visualized in context using IGB (Freese et al., 2016), or downloaded directly, using our publicly available IGB quickload server at https://compbio.lifesci.dundee.ac.uk/arabidopsisQuickload/public-quickload/ under 'RNAseq>Froussios2019'. All scripts and commands are available from github at https://github.com/bartongroup/KF-arabidopsis-GRNA. Supplementary information: Supplementary data are available at Bioinformatics online.

Original languageEnglish
Pages (from-to)3372-3377
Number of pages6
JournalBioinformatics
Volume35
Issue number18
Early online date6 Feb 2019
DOIs
Publication statusPublished - 15 Sep 2019

Fingerprint Dive into the research topics of 'How well do RNA-Seq differential gene expression tools perform in a complex eukaryote? A case study in Arabidopsis thaliana'. Together they form a unique fingerprint.

  • Projects

    Cite this