Naught all zeros in sequence count data are the same

Justin D. Silverman, Kimberly Roche, Sayan Mukherjee, Lawrence A. David

    Research output: Contribution to journalReview articlepeer-review

    26 Scopus citations


    Genomic studies feature multivariate count data from high-throughput DNA sequencing experiments, which often contain many zero values. These zeros can cause artifacts for statistical analyses and multiple modeling approaches have been developed in response. Here, we apply different zero-handling models to gene-expression and microbiome datasets and show models can disagree substantially in terms of identifying the most differentially expressed sequences. Next, to rationally examine how different zero handling models behave, we developed a conceptual framework outlining four types of processes that may give rise to zero values in sequence count data. Last, we performed simulations to test how zero handling models behave in the presence of these different zero generating processes. Our simulations showed that simple count models are sufficient across multiple processes, even when the true underlying process is unknown. On the other hand, a common zero handling technique known as “zero-inflation” was only suitable under a zero generating process associated with an unlikely set of biological and experimental conditions. In concert, our work here suggests several specific guidelines for developing and choosing state-of-the-art models for analyzing sparse sequence count data.

    Original languageEnglish (US)
    Pages (from-to)2789-2798
    Number of pages10
    JournalComputational and Structural Biotechnology Journal
    StatePublished - Jan 2020

    All Science Journal Classification (ASJC) codes

    • Biotechnology
    • Biophysics
    • Structural Biology
    • Biochemistry
    • Genetics
    • Computer Science Applications


    Dive into the research topics of 'Naught all zeros in sequence count data are the same'. Together they form a unique fingerprint.

    Cite this