Abstract:It is crucial for researchers to optimize RNA-seq experimental designs for differential expression detection. Currently, the field lacks general methods to estimate power and sample size for RNA-Seq in complex experimental designs, under the assumption of the negative binomial distribution. We simulate RNA-Seq count data based on parameters estimated from six widely different public data sets (including cell line comparison, tissue comparison, and cancer data sets) and calculate the statistical power in paired and unpaired sample experiments. We comprehensively compare five differential expression analysis packages (DESeq, edgeR, DESeq2, sSeq, and EBSeq) and evaluate their performance by power, receiver operator characteristic (ROC) curves, and other metrics including areas under the curve (AUC), Matthews correlation coefficient (MCC), and F -measures. DESeq2 and edgeR tend to give the best performance in general. Increasing sample size or sequencing depth increases power; however, increasing sample size is more potent than sequencing depth to increase power, especially when the sequencing depth reaches 20 million reads. Long intergenic noncoding RNAs (lincRNA) yields lower power relative to the protein coding mRNAs, given their lower expression level in the same RNA-Seq experiment. On the other hand, paired-sample RNA-Seq significantly enhances the statistical power, confirming the importance of considering the multifactor experimental design. Finally, a local optimal power is achievable for a given budget constraint, and the dominant contributing factor is sample size rather than the sequencing depth. In conclusion, we provide a power analysis tool ( http://www2.hawaii.edu/~lgarmire/RNASeqPowerCalculator.htm ) that captures the dispersion in the data and can serve as a practical reference under the budget constraint of RNA-Seq experiments.

Power analysis and sample size estimation for sequence-based association studies

Statistical power and significance testing in large-scale genetic studies

Impact of genotyping errors on statistical power of association tests in genomic analyses: A case study.

Biases and errors on allele frequency estimation and disease association tests of next-generation sequencing of pooled samples.

A Robust and Powerful Set-Valued Approach to Rare Variant Association Analyses of Secondary Traits in Case-Control Sequencing Studies.

scPower accelerates and optimizes the design of multi-sample single cell transcriptomic studies

Powerful Rare‐Variant Association Analysis of Secondary Phenotypes

Power Analysis of Principal Components Regression in Genetic Association Studies.

Meta-SAIGE: Scalable and Accurate Meta-Analysis for Rare Variants

A Robust Model-free Approach for Rare Variants Association Studies Incorporating Gene-Gene and Gene-Environmental Interactions

SAIGE-GENE+ improves the efficiency and accuracy of set-based rare variant association tests

Power analysis and sample size estimation for RNA-Seq differential expression

Robust Reference Powered Association Test of Genome-Wide Association Studies

Computationally efficient meta-analysis of gene-based tests using summary statistics in large-scale genetic studies

Enhancing power of rare variant association test by Zoom-Focus Algorithm (ZFA) to locate optimal testing region

Extreme Sampling Design in Genetic Association Mapping of Quantitative Trait Loci Using Balanced and Unbalanced Case-Control Samples

Pooled Association Tests for Rare Genetic Variants: A Review and Some New Results

Alternative Methods for H1 Simulations in Genome Wide Association Studies

Power and sample size calculations for testing the ratio of reproductive values in phylogenetic samples

Linear Models for Analysis of Multiple Single Nucleotide Polymorphisms with Quantitative Traits in Unrelated Individuals

RnaSeqSampleSize: real data based sample size estimation for RNA sequencing