gene expression is the appearance of a characteristic or effect in the phenotype that can be attributed to a certain gene . for over a decade , microarray technology was the dominating technology for high - throughput gene expression profiling until the introduction of rnaseq technology .
rnaseq technology is a form of next - generation sequencing ( ngs ) technology , sometimes referred to as high - throughput sequencing technology . since its introduction
rnaseq is the use of ngs technology to sequence cdna ( reversed transcribed from rna ) in order to obtain information about rna .
first , rnaseq allows the detection of all isoforms of a gene , even novel ones .
microarray , on the other hand , relies purely on previous knowledge regarding genes to design probes for detection , thus it can not be used for novel detection .
second , the resolution of microarray usually stays at the gene and exon level , but the resolution of rnaseq can reach the level of a single nucleotide , which allows detection of single nucleotide variance and structural variants , such as small insertions , deletions , alternative splicing , and gene fusion . although the cost of rnaseq has become comparable to microarray , researchers have unanimously agreed that rnaseq has replaced microarray as the go to technology of high - throughput gene expression profiling.14 although rnaseq technology has introduced exciting opportunities to biomedical researchers , it has also created several challenges for analysts . for example , the data storage cost is significantly higher compared to microarray . and
the lack of consensus on the best statistical method for detecting differentially expressed genes5 is really noticeable .
another easily overlooked area is sample size and power calculation for rnaseq - based experiment designs .
experiment design is a pivotal step and the first step of any successful study , and one of the most important questions to address during experiment design is what sample size is needed to achieve the desired statistical power within the financial budget .
sample size and power analyses are required for most national institutes of health ( nih ) funding applications . for microarray - based studies ,
sample size and power calculation has been well established.68 the sample size and power for microarray - based studies are relatively easy to compute , because gene expression of microarray data follows a normal distribution .
however , rnaseq data are count based , which is usually modeled as a poisson9,10 or negative binomial distribution.11,12 several studies have attempted to produce methods to estimate sample size and power based on poisson and negative binomial distributions for rnaseq - based experiment design .
for example , fang and cui13 derived a sample size calculation formula based on the wald test for a single - gene differential expression test .
busby et al.14 introduced scotty , a web tool for computed sample size using a t - test for rnaseq data based on a poisson - lognormal distribution .
hart et al.15 introduced rnaseqpower , a method for sample size calculation using the score test based on a negative binomial distribution .
r code and the excel worksheet for this method were provided by the authors . in reality ,
rnaseq experiment is capable of detecting thousands of genes , and those genes are tested for the significance of differential expression simultaneously . in such cases , the correction of error rates for multiple comparisons is required .
however , both scotty and rnaseqpower fail to discuss how to calculate sample size under multiple comparison testing . to address this issue ,
later , li et al.17 proposed a sample size calculation method based on the exact test.18 of all the developed methods , only scotty has a web - based graphical user interface .
a convenient , user - friendly rnaseq sample size and power calculation tool is highly desirable . here
, we present rnaseqps , a web - based power and sample size calculation tool , based on both the poisson and negative binomial distributions .
rnaseqps distinguishes itself from other rnaseq sample size and power calculation methods in several aspects .
also , it does not require the skill of running r code , a person without any previous experience in r can operate it with ease .
there are other features of rnaseqps which are discussed in detail in the method section .
for the poisson distribution , we considered six different methods : ( 1 ) wald test ; ( 2 ) score test ; ( 3 ) likelihood ratio test ; ( 4 ) log transformation of wald statistic ; ( 5 ) log transformation of score statistic ; and ( 6 ) transformation of poisson for calculating sample size under both single- and multiple - gene comparison scenarios .
the detail of the implementation of these tests can be found.16 under the poisson distribution , we assume that mean equals variance .
when biological replicates are used , rnaseq data could exhibit variation significantly greater than the mean ( over - dispersion ) . in such scenario ,
the negative binomial distribution can be used as a natural extension of the poisson model . for the negative binomial distribution ,
we provide a method to calculate power for the exact test proposed.18 based on a given p value under the single - gene comparison scenario .
the detail of the implementation of this test can be found.17 to address the multiple comparison issue , we further incorporated false discovery rate ( fdr)19 controlling into our methods described.16,17 rnaseqps implements sample size and power calculation from both poisson and negative binomial distributions based on the methods proposed.16,17
rnaseqps is written using the shiny package in r. shiny is a powerful web framework for building interactive web applications using r. it combines the computation power of r with the interactivity of modern web .
applications built using shiny can automatically react with the input and output parameters of an application .
thus , when new parameters are input into rnaseqps , sample size or power is automatically recomputed and displayed on the screen .
users without any experiences in statistics and programing languages can navigate through it without trouble .
for example , a user can specify sample size as the x - axis and selects a reasonable range for sample size .
rnaseqps can plot relationship of the sample size with power while other parameters stay unchanged .
such plots help users visualize the intrinsic relationship between the parameters and power and can be easily exported as high - resolution figures that can be used for nih grant funding applications .
each computation of sample size or power takes up to 30 seconds to 1 minute of time . to reduce the user wait time
the database lookup function works as follows . to estimate the power of the rnaseq experiment under the poisson distribution , the following eight parameters are required as input : ( 1 ) sample size ; ( 2 ) expected fold change ; ( 3 ) average reads per gene ; ( 4 ) read count ratio between two conditions ; ( 5 ) fdr threshold ; ( 6 ) total number of genes ; ( 7 ) expected number of significant genes and ( 8) dispersion ( only required when using negative binomial distribution ) .
the details of these parameters are explained in table 1 . for each of these parameters ,
for example , for expected fold change parameter , the range is usually between absolute value 1 to 10 because the absolute value of the fold change can not be less than 1 and fold change 10 is a reasonable assumption of the high ceiling . to produce the database
, we pre - computed the power and sample size using the parameter values within these most likely ranges for all parameters and stored them in a my sql database . for each round of input parameters , rnaseqps first queries the database to see
if samples size or power for such set of parameters has already been entered into the database . if the input parameters match the existing parameters in the database , the power or sample size is automatically returned , otherwise power or sample size will be computed on the fly .
when using slide bar , the input is guaranteed to match an input in the database .
in such cases , the sample size and power will be computed on the fly . by using the look up function
the dominance of rnaseq over microarray in gene expression research has been documented by multiple studies.14 even though the price of rnaseq has matched that of microarray , large studies can still be costly .
sample size estimation is a critical issue for any biological studies , including rnaseq - based studies .
better management of the tradeoff between cost and sample size is key to a successful study .
sample size and power analysis for rnaseq is an underdeveloped area when compared with microarray . in comparison , sample size and power
is easier to calculate for microarray , because microarray follows a normal distribution , while rnaseq data are count based and follow a poisson or negative binomial distribution .
these methods adequately address the sample size and power calculation needs for an rnaseq experiment design but still have two issues .
rnaseqps addresses the multiple comparisons problem by implementing the methods introduced by li et al . , which incorporate fdr controls . also , rnaseqps implements methods based on both the poisson and negative binomial distributions .
this graphical user interface helps users input the parameters with ease and greatly enhances users abilities to understand the sample size and power analysis with the help of visualization .
overall , rnaseqps addresses major issues currently lacking in sample size and power analysis for rnaseq experiment design .
it greatly benefits researchers who plan to use rnaseq as the primary technology for their study .
 sample size and power determination is the first step in the experimental design of a successful study . 
 sample size and power calculation is required for applications for national institutes of health ( nih ) funding . 
 sample size and power calculation is well established for traditional biological studies such as mouse model , genome wide association study ( gwas ) , and microarray studies . 
 recent developments in high - throughput sequencing technology have allowed rnaseq to replace microarray as the technology of choice for high - throughput gene expression profiling . 
 however , the sample size and power analysis of rnaseq technology is an underdeveloped area . here , we present rnaseqps , an advanced online rnaseq power and sample size calculation tool based on the poisson and negative binomial distributions . 
 rnaseqps was built using the shiny package in r. it provides an interactive graphical user interface that allows the users to easily conduct sample size and power analysis for rnaseq experimental design . 
 rnaseqps can be accessed directly at http://cqs.mc.vanderbilt.edu / shiny / rnaseqps/. 
Introduction
Method
Results
Discussion