RnaSeqSampleSize: real data based sample size estimation for RNA sequencing.

Area of Science:

• Genomics
• Bioinformatics
• Statistical Genetics

Background:

• Sample size estimation is critical but often overlooked in RNA sequencing (RNA-Seq) studies.
• Existing methods focus on single genes, neglecting the complexities of multi-gene analysis in RNA-Seq, such as false discovery rates and variable gene expression.
• RNA-Seq experiments involve simultaneous quantification and differential expression testing of thousands of genes, necessitating robust sample size determination.

Purpose of the Study:

• To develop a novel method for accurate sample size and power estimation in RNA sequencing experiments.
• To address the challenges posed by multiple testing and the wide distribution of read counts and dispersions across genes in RNA-Seq data.
• To provide a user-friendly tool for researchers to optimize their experimental design.

Main Methods:

• Developed RnaSeqSampleSize, a method based on the distributions of average read counts and dispersions from real RNA-Seq data.
• Utilized reference datasets, such as The Cancer Genome Atlas (TCGA), to estimate gene read counts and dispersions.
• Implemented the method in R and made it accessible via Bioconductor, with a web-based graphical interface.

Main Results:

• RnaSeqSampleSize enables estimation of power and sample size by leveraging distributions from existing RNA-Seq datasets.
• The tool provides a practical approach to address the complexities of sample size calculation in multi-gene RNA-Seq analyses.
• The R implementation and web interface facilitate easy adoption by researchers.

Conclusions:

• RnaSeqSampleSize offers a powerful and convenient solution for RNA-Seq power and sample size estimation.
• Unique features include estimation for specific genes or pathways, power curve visualization, and parameter optimization.
• The tool enhances the reliability and efficiency of RNA-Seq experimental design.