BACKGROUND
the transcriptome catalogues the complete set of transcripts in a cell.
transcriptomic regulation is critical to all physiological, developmental and pathological processes  <cit> , and mrna expression profiles can represent the characteristics of a cell at a specific state and help to govern its present and future activities  <cit> .
the profiles of a transcriptome in terms of alterations in response to specific biological stimuli provides valuable insights for interpreting functional elements of the genome, revealing the molecular constituents of cells, and also understanding developmental and disease processes.
different types of technologies have been developed to interrogate transcript abundance, including hybridization-based and sequencing-based approaches.
hybridization-based microarrays have been the primary transcriptomic high-throughput tool for almost two decades, which has accelerated the study of transcriptome analysis by profiling thousands of genes simultaneously  <cit> .
however, microarray technology has several limitations including: indirect quantification by hybridization-signal intensities  <cit> , background and cross-hybridization problems  <cit>  and reproducibility issues  <cit> .
the development of next generation sequencing with improved qualitative and quantitative measurements holds great promise in transcriptome analysis.
rna-seq is a recently developed approach to map and quantify transcriptomes by digitally recording how frequently each transcript is represented in a sequence sample.
after poly  selection, rna is fragmented to small fragments and converted into a cdna library, which provides a simple and more comprehensive way to measure transcriptome composition and to discover new genes by high-throughput sequencing without bacterial cloning of cdna input  <cit> .
studies using this technology have already altered our views regarding the extent and complexity of transcriptomes in an organism and dramatically improved our understanding of transcriptome.
rna-seq has several advantages over micorarrays including: 1) rna-seq is not dependent on prior knowledge about the target sequence; 2) it has a large dynamic range and sensitivity due to its digital nature, which is especially important for highly abundant and extremely low abundant genes; 3) the survey of a transcriptome is more accurate because the quantification of each transcript is directly based on digital counts of the transcript.
therefore, rna-seq offers both single-base resolution for annotation and “digital” quantification at the rna level, which allows the entire transcriptome to be analyzed in a high-throughput and quantitative manner  <cit> .
however, the expense per sample for rna- seq is still a limiting factor in preventing researchers from sequencing multiple biological replicates per group, which are needed for statistically-significant analysis.
it is common to adopt a pooling strategy to reduce the cost for rna-seq studies  <cit> .
with the continued enhancement of sequencing output and the development of multiplex labelling techniques, the cost per sample could be significantly reduced if several samples are multiplexed and sequenced in the same lane, given sufficient transcriptome coverage per sample.
therefore, it is imperative to address the trade-off between the depth of rna-seq and the coverage of the transcriptome in an organism.
the objective of this study was to evaluate what coverage or sequencing depth of transcriptome would be sufficient to interrogate gene expression profiling in the chicken by rna-seq.
methods
rna preparation
total rna was isolated from four chicken lungs from two genetic chicken lines  by trizol  according to the manufacturer’s protocol.
two rna samples from the same line were pooled to generate totals of two pooled rna samples .
dnase i  digestion was carried out after rna isolation and the rna concentration and purity were determined by measuring absorbance at  <dig> nm and a260/a <dig> ratio using a nanodrop nd- <dig> spectrophotometer .
rna samples were stored at - <dig> °c until further use.
cdna library preparation and sequencing by rna-seq
total rna  was subjected to two rounds of hybridization to oligo  beads  to enrich mrna.
ribosomal rna contamination was evaluated by rna pico chip using a bioanalyzer .
the resulting mrna was then used to prepare cdna libraries using the rna sequencing sample preparation kit .
sample <dig> and sample <dig> were sequenced by illumina genome analyzer and then genome analyzer ii, which generated four datasets: s1-r <dig>  s2-r <dig>  and s1-r <dig>  s2-r <dig>  respectively.
data filtering, mapping reads and identifying transcriptome contents
the sequences generated went through a filtering process first.
any reads that contained numerous interspersed ns in their sequences, or had relatively short reads , were removed for the following analysis.
sequence reads obtained after quality control with filtering were analyzed using clc genomics workbench  <dig> .
after mapping, the gene expression level was quantified by simply dividing the number of reads mapped to each gene by the size of its transcripts, commonly known as the number of reads per kilobase of exon per million mapped reads   <cit> , for all  <dig>  annotated chicken genes in the database.
the gene expression level was then log <dig> transformed.
random sampling of s1-r <dig> and s2-r2
we have obtained rna-seq data in three different levels of depth:  <dig>  m,  <dig>  m, and about  <dig> m reads.
clearly, there was a big gap between  <dig>  m and  <dig> m reads.
in order to identify the appropriate depth of transcriptome per sample that is sufficient for whole genome transcriptome profiling, it is important to generate additional datasets at different levels of depth.
it would be very costly to re-sequence each sample to generate rna-seq data at different levels of sequencing depth.
random sampling from current dataset might provide a cost-effective approach for this purpose.
this procedure synthetically created samples from the originally sequenced samples.
thus, for samples s1-r <dig> and s2-r <dig> data sets, by drawing without replacement a fixed number of reads from the overall data set, we randomly selected  <dig> m,  <dig> m and finally  <dig> m reads.
these random selections are repeated  <dig> times, resulting in total of  <dig> technical replicates with different transcriptome depth.
each one of the reads in the fastq format of the input files, which were used in the sampling, was selected equally likely.
a program in perl was written to serve this purpose.
then, the resulting replicate datasets were uploaded into clc genomics workbench in the fastq format for the analysis individually.
the correlation coefficients of gene expression  between replicates of each sequencing depth were calculated by jmp .
the average rpkms of transcripts identified by each sequencing depth  were calculated to represent the gene expression values for further analysis.
correlation coefficients between different sequencing depths from the same sample
in order to evaluate how reliable the sequence data is at each level of sequencing depth, correlation coefficients between lower and high depth sequence data for each sample were calculated by jmp .
any genes with no gene expression at either sequence dataset were excluded from the correlation coefficient computation.
the depth of sequencing is highly correlated with the abundance of gene expression, so we divided genes into the four quartile groups based on expression levels for each dataset, from the bottom 25%  to the top 25% .
subsequently, correlation coefficients between lower-depth sequence data and the data with  <dig> -  <dig>  m reads within each of the four quartile groups were calculated.
RESULTS
rna-seq for cdna libraries
the two chicken cdna libraries  were sequenced by the illumina genome analyzer, which generated  <dig>  m  reads  and  <dig>  m  reads , respectively.
then, two technical replicate cdna libraries from the same rna samples were re-sequenced using the genome analyzer ii, which generated  <dig>  m  reads  and  <dig>  m  reads .
random sampling of the s1-r <dig> and s2-r2
the datasets of s1-r <dig>  and s2-r <dig>  were each randomly re-sampled into  <dig> m,  <dig> m, and  <dig> m reads with four replicates each.
the correlation coefficients between every two replicates for each re-sampled level  within each sample  were all greater than  <dig> , which demonstrated that the sampling procedure is consistent and accurate.
averages of gene expression from the four replicates at each re-sampled level for each sample were used for further analysis.
effects of sequence depth on the distributions of transcripts
the distributions of transcript abundance at different levels of sequence depth from sample <dig> and sample <dig> are presented in figures  <dig> and  <dig>  respectively.
in general, the median and 75% percentile were similar across five different levels of depth, while the 95% percentile showed a slight increase, especially from  <dig> m to  <dig> m. significant decreases at the 25% and 5% percentile were observed, especially from  <dig> m to  <dig> m. in addition, a significant decrease was also observed from  <dig>  m to  <dig> m in sample <dig>
coverage of annotated chicken genes
there are about  <dig>  annotated chicken genes in the ncbi database.
number of detected chicken genes at different levels of sequencing depth from samples  <dig> and  <dig> are presented in table  <dig>  there were  <dig>  genes detected in s1-r <dig>  and  <dig>  genes in s2-r <dig> , which accounted for  <dig> % and  <dig> % of all  <dig>  annotated chicken genes in the database, respectively.
with the reduction of sequencing depth, the number of detectable genes also decreased from 91% to 78% in sample <dig> , and from 90% to 68% in sample <dig> .
two significant drops of transcriptome coverage were observed: from  <dig> m to  <dig> m, and  <dig> m to  <dig>  m.
sample1
sample2
numbers of genes
numbers of genes
correlation coefficients between different sequencing depths
to evaluate the appropriate depth of sequence that is needed for transcriptome profiling, correlation coefficients between different levels of sequencing depth and the most abundant reads for each sample were calculated.
for sample <dig>  overall correlation coefficients at four different levels of depth were greater than  <dig> .
when we examined the four quartile groups based on expression level , correlation coefficients ranged from  <dig>  to  <dig>  for the 1st quartile,  <dig>  to  <dig>  for the 2nd quartile,  <dig>  to  <dig>  for the 3rd quartile, and  <dig>  to  <dig>  for the 4th quartile.
a similar pattern in terms of correlation coefficient change was observed for the 1st, 2nd, and 3rd quartiles; a significantly increased correlation coefficient from  <dig>  m to  <dig> m, and kept relative flat from  <dig> m to  <dig> m. for the 4th quartile, correlation coefficients at four different levels of depth were greater than  <dig> .
from the 1st to the 4th quartiles, there were significant increases for correlation coefficients between every two quartile groups .
for sample <dig>  overall correlation coefficients at four different levels of depth were greater than  <dig> , except for  <dig>  m at  <dig> .
correlation coefficients ranged from  <dig>  to  <dig>  for the 1st quartile,  <dig>  to  <dig>  for the 2nd quartile,  <dig>  to  <dig>  for the 3rd quartile, and  <dig>  to  <dig>  for the 4th quartile .
the same pattern in terms of correlation coefficients changes at different levels of depth between sample <dig> and sample <dig> was also observed.
discussion
in the current study, rna-seq was performed twice using two chicken cdna samples.
the first run of rna-seq had fewer number of reads and larger variation in terms of total number of reads between the two samples, while the second run had greater number of reads and very small variation between the two samples.
the first run was performed at the very early stage of the sequencing technology when it was still in the testing phase.
the lower reads and larger variation in the first run may be coming from two major sources of technical error: the purification of cdna templates during the library preparation, and the loading of libraries onto flow cells .
these potential sources of errors were corrected during the second rna-seq analysis, which provided very good sequencing depth with greater number of reads.
the first rna-seq datasets were directly derived from actual experiment, which made the results more informative than replicating datasets by random sampling.
therefore, we chose to include these two early datasets in the analysis in the current study.
furthermore, all of the reads from each sample were normalized by the rpkm and the datasets can serve as a reference for random sampling at different sequencing depths from the exact same samples.
the capacity of sequencing length of  <dig> bp for the first run was increased to  <dig> bp at the second rna-seq analysis.
longer reads should reduce estimation error and mapping uncertainty, and read lengths have been consistently increasing with improving illumina massively parallel sequencing technology.
however, people have noted that the number of reads is more important than the read length once reaching a minimum read length of  <dig> bp  <cit> .
the read lengths  in the current study were larger than  <dig> bp; therefore, the read length will not affect the overall conclusions drawn.
as a powerful new technology for transcriptome analysis, rna-seq provides a more comprehensive view of the transcriptome than earlier technologies.
besides its ability to detect splicing variation, rna editing and discovery of new transcripts  <cit> , rna-seq can also function in the role of a conventional microarray in measuring gene expression due to its accurate measurements.
in order to detect less abundant transcripts, appropriate sequencing depth is needed.
the transcriptome coverage or sequencing depth needed for a given study can be affected by several factors such as genome size, transcriptome complexity and objectives of the study.
in general, the more complex the transcriptome, the more sequencing depth is required for adequate coverage  <cit> .
depending on the purpose of transcriptome analysis, the requirement of sequencing depth varies.
in most transcriptome studies, quantifying mrna abundance is one of the major objectives.
there is a certain sequencing depth that is sufficient in simple transcriptomes.
for example, in the yeast genome, a  <dig>  m  reads dataset was generated by rna-seq which was able to get 100% transcriptome coverage  <cit> .
the number of transcripts detected by rna-seq in the yeast dataset was able to reach 80% transcriptome coverage at 4m mapped reads, and even though the sequencing depth doubled as 8m reads, the transcriptome coverage only increased 10%  <cit> .
these results suggest that the improvement of sequencing depth or transcriptome coverage after reaching a certain sequencing depth had relatively less impact on detecting low abundant genes  <cit> .
in addition, the cost per sample per lane by rna-seq is still not affordable for most laboratories.
recent development in multiplex labelling using bar-coded libraries by illumina and continued increase in sequence output have made it possible to sequence multiple samples per lane without extra cost or running time  <cit> .
therefore, it is imperative to examine the correlation between sequencing depth and transcriptome coverage; in other words, what sequencing depth might be sufficient in reaching a certain level of transcriptome coverage and reliable measurement for rna-seq.
in order to accomplish this objective, two approaches could be applied: experimental or simulation methods.
both methods have been applied in this study.
high correlation among replicates within each sequencing depth, gradual increase in correlation coefficients from  <dig> m to  <dig> m, and consistent patterns observed between samples <dig> and  <dig>  have demonstrated that random sampling was an effective and reliable method in reaching the goals of this study.
transcriptome coverage is one of the most important parameters in profiling global gene expression.
the number and level of transcript isoforms is not always known and transcription activity varies across the genome  <cit> .
this was confirmed in a study by using the number of unique transcription start sites as a measure of coverage in mouse cells  <cit> .
in the current study, we took a more practical approach using all annotated genes in the chicken genome.
because the chicken genome is far under-annotated, we assume that the  <dig>  annotated chicken genes in the database would well represent different levels of expression abundance in the chicken genome, which is essential for the analysis of transcriptome coverage in this study.
since gene expression depends on tissue and time of biological process  <cit> , it is impossible for any single tissue to have all genes in the genome expressed.
ninety percent of all annotated genes  detected at about  <dig> m  reads might represent a saturated detection of the whole genome.
the analysis results indicated significant improvements of transcriptome coverage occurred from  <dig>  m to  <dig>  m and from  <dig> m to  <dig> m. depending on the purpose of transcriptome analysis, the current study suggested that  <dig> m  reads could have 80% of transcriptome coverage, while  <dig> m  reads could reach 90% of transcriptome coverage.
when we analyzed overall correlation coefficients at different levels of sequencing depth regardless of gene expression level, we observed very high correlation coefficients between each level of sequencing depth compared with  <dig> m, except for  <dig>  m. one might draw a conclusion that there is no significant difference among different levels of sequencing depth.
but as we see in figure  <dig>  this might be true in the case of highly abundant genes , but not in the case of the 1st to 3rd quartile groups, especially the first two quartile groups .
high abundant genes will be less affected by sequencing depth than low abundant genes, because high abundant genes are more likely to be captured even when the sequencing depth varies  <cit> .
this is also confirmed by our analysis.
collectively, the following points can be inferred: 1) sequencing depth below  <dig> m  reads had a significant effect on detecting transcript levels of medium and low abundant transcripts; 2) sequencing depth at both  <dig>  m and  <dig>  m would result in unreliable mrna expression on all genes except highly abundant transcripts; 3) there was no significant improvement in correlation coefficients when sequencing depth doubled from  <dig> m to  <dig> m. based on these analysis, the results suggested: 1)  <dig> m reads might be sufficient in obtaining reliable mrna expression measurement on highly abundant transcripts; 2) when sequencing depth is beyond  <dig> m reads, a relatively reliable measurement of mrna expression is expected, especially for abundant transcripts; 3) it seems that  <dig> m of reads is needed to achieve reliable measurement of mrna expression across all genes in the chicken genome.
to our knowledge, this is the first study evaluating the appropriate sequencing depth using rna-seq in farm animals and will provide the first reference for similar studies.
the knowledge generated from this study has laid a solid foundation for applying this analysis to other species.
CONCLUSIONS
in summary, the analysis from this study demonstrated that  <dig> m  reads is sufficient to detect all annotated genes and provide a reliable measurement of mrna abundance in chicken lungs using rna-seq.
as we expected, the depth of sequencing had significant impact on low abundance transcripts, but not on high abundance transcripts.
in practice, if rna-seq serves as a replacement of microarray technology,  <dig> m  reads would allow detection of about 80% of annotated chicken genes.
finally, increasing the depth of sequencing from  <dig> m to  <dig> m reads did not have a significant effect on transcriptome coverage and reliability of mrna measurements, whereas  <dig> m reads was needed to achieve reliable measurement of mrna expression across all genes in the genome.
competing interests
the authors declare that there are no competing interests.
authors' contributions
yw analyzed the data and drafted the manuscript, ng wrote the code for sampling the rna-seq data, cj and ub contributed ideas for the data analysis and revised the manuscript, hw and rc ran rna-seq samples, and hz designed the experiment, provided the concept of the analysis.
all authors submitted comments, read and approved the final manuscript.

 the transcriptome catalogues the complete set of transcripts in a cell. transcriptomic regulation is critical to all physiological, developmental and pathological processes, and mrna expression profiles can represent the characteristics of a cell at a specific state and help to govern its present and future activities. the profiles of a transcriptome in terms of alterations in response to specific biological stimuli provides valuable insights for interpreting functional elements of the genome, revealing the molecular constituents of cells, and also understanding developmental and disease processes. different types of technologies have been developed to interrogate transcript abundance, including hybridization-based and sequencing-based approaches. hybridization-based microarrays have been the primary transcriptomic high-throughput tool for almost two decades, which has accelerated the study of transcriptome analysis by profiling thousands of genes simultaneously. however, microarray technology has several limitations including: indirect quantification by hybridization-signal intensities, background and cross-hybridization problems and reproducibility issues. the development of next generation sequencing with improved qualitative and quantitative measurements holds great promise in transcriptome analysis. rna-seq is a recently developed approach to map and quantify transcriptomes by digitally recording how frequently each transcript is represented in a sequence sample. after poly selection, rna is fragmented to small fragments and converted into a cdna library, which provides a simple and more comprehensive way to measure transcriptome composition and to discover new genes by high-throughput sequencing without bacterial cloning of cdna input. studies using this technology have already altered our views regarding the extent and complexity of transcriptomes in an organism and dramatically improved our understanding of transcriptome. rna-seq has several advantages over micorarrays including: 1) rna-seq is not dependent on prior knowledge about the target sequence; 2) it has a large dynamic range and sensitivity due to its digital nature, which is especially important for highly abundant and extremely low abundant genes; 3) the survey of a transcriptome is more accurate because the quantification of each transcript is directly based on digital counts of the transcript. therefore, rna-seq offers both single-base resolution for annotation and “digital” quantification at the rna level, which allows the entire transcriptome to be analyzed in a high-throughput and quantitative manner. however, the expense per sample for rna- seq is still a limiting factor in preventing researchers from sequencing multiple biological replicates per group, which are needed for statistically-significant analysis. it is common to adopt a pooling strategy to reduce the cost for rna-seq studies. with the continued enhancement of sequencing output and the development of multiplex labelling techniques, the cost per sample could be significantly reduced if several samples are multiplexed and sequenced in the same lane, given sufficient transcriptome coverage per sample. therefore, it is imperative to address the trade-off between the depth of rna-seq and the coverage of the transcriptome in an organism. the objective of this study was to evaluate what coverage or sequencing depth of transcriptome would be sufficient to interrogate gene expression profiling in the chicken by rna-seq. methods rna preparation total rna was isolated from four chicken lungs from two genetic chicken lines by trizol according to the manufacturer’s protocol. two rna samples from the same line were pooled to generate totals of two pooled rna samples. dnase i digestion was carried out after rna isolation and the rna concentration and purity were determined by measuring absorbance at nm and a260/a ratio using a nanodrop nd- spectrophotometer. rna samples were stored at - °c until further use. cdna library preparation and sequencing by rna-seq total rna was subjected to two rounds of hybridization to oligo beads to enrich mrna. ribosomal rna contamination was evaluated by rna pico chip using a bioanalyzer. the resulting mrna was then used to prepare cdna libraries using the rna sequencing sample preparation kit. sample and sample were sequenced by illumina genome analyzer and then genome analyzer ii, which generated four datasets: s1-r s2-r and s1-r s2-r respectively. data filtering, mapping reads and identifying transcriptome contents the sequences generated went through a filtering process first. any reads that contained numerous interspersed ns in their sequences, or had relatively short reads, were removed for the following analysis. sequence reads obtained after quality control with filtering were analyzed using clc genomics workbench. after mapping, the gene expression level was quantified by simply dividing the number of reads mapped to each gene by the size of its transcripts, commonly known as the number of reads per kilobase of exon per million mapped reads, for all annotated chicken genes in the database. the gene expression level was then log transformed. random sampling of s1-r and s2-r2 we have obtained rna-seq data in three different levels of depth: m, m, and about m reads. clearly, there was a big gap between m and m reads. in order to identify the appropriate depth of transcriptome per sample that is sufficient for whole genome transcriptome profiling, it is important to generate additional datasets at different levels of depth. it would be very costly to re-sequence each sample to generate rna-seq data at different levels of sequencing depth. random sampling from current dataset might provide a cost-effective approach for this purpose. this procedure synthetically created samples from the originally sequenced samples. thus, for samples s1-r and s2-r data sets, by drawing without replacement a fixed number of reads from the overall data set, we randomly selected m, m and finally m reads. these random selections are repeated times, resulting in total of technical replicates with different transcriptome depth. each one of the reads in the fastq format of the input files, which were used in the sampling, was selected equally likely. a program in perl was written to serve this purpose. then, the resulting replicate datasets were uploaded into clc genomics workbench in the fastq format for the analysis individually. the correlation coefficients of gene expression between replicates of each sequencing depth were calculated by jmp. the average rpkms of transcripts identified by each sequencing depth were calculated to represent the gene expression values for further analysis. correlation coefficients between different sequencing depths from the same sample in order to evaluate how reliable the sequence data is at each level of sequencing depth, correlation coefficients between lower and high depth sequence data for each sample were calculated by jmp. any genes with no gene expression at either sequence dataset were excluded from the correlation coefficient computation. the depth of sequencing is highly correlated with the abundance of gene expression, so we divided genes into the four quartile groups based on expression levels for each dataset, from the bottom 25% to the top 25%. subsequently, correlation coefficients between lower-depth sequence data and the data with - m reads within each of the four quartile groups were calculated.  rna-seq for cdna libraries the two chicken cdna libraries were sequenced by the illumina genome analyzer, which generated m reads and m reads, respectively. then, two technical replicate cdna libraries from the same rna samples were re-sequenced using the genome analyzer ii, which generated m reads and m reads. random sampling of the s1-r and s2-r2 the datasets of s1-r and s2-r were each randomly re-sampled into m, m, and m reads with four replicates each. the correlation coefficients between every two replicates for each re-sampled level within each sample were all greater than, which demonstrated that the sampling procedure is consistent and accurate. averages of gene expression from the four replicates at each re-sampled level for each sample were used for further analysis. effects of sequence depth on the distributions of transcripts the distributions of transcript abundance at different levels of sequence depth from sample and sample are presented in figures and respectively. in general, the median and 75% percentile were similar across five different levels of depth, while the 95% percentile showed a slight increase, especially from m to m. significant decreases at the 25% and 5% percentile were observed, especially from m to m. in addition, a significant decrease was also observed from m to m in sample coverage of annotated chicken genes there are about annotated chicken genes in the ncbi database. number of detected chicken genes at different levels of sequencing depth from samples and are presented in table there were genes detected in s1-r and genes in s2-r, which accounted for % and % of all annotated chicken genes in the database, respectively. with the reduction of sequencing depth, the number of detectable genes also decreased from 91% to 78% in sample, and from 90% to 68% in sample. two significant drops of transcriptome coverage were observed: from m to m, and m to m. sample1 sample2 numbers of genes numbers of genes correlation coefficients between different sequencing depths to evaluate the appropriate depth of sequence that is needed for transcriptome profiling, correlation coefficients between different levels of sequencing depth and the most abundant reads for each sample were calculated. for sample overall correlation coefficients at four different levels of depth were greater than. when we examined the four quartile groups based on expression level, correlation coefficients ranged from to for the 1st quartile, to for the 2nd quartile, to for the 3rd quartile, and to for the 4th quartile. a similar pattern in terms of correlation coefficient change was observed for the 1st, 2nd, and 3rd quartiles; a significantly increased correlation coefficient from m to m, and kept relative flat from m to m. for the 4th quartile, correlation coefficients at four different levels of depth were greater than. from the 1st to the 4th quartiles, there were significant increases for correlation coefficients between every two quartile groups. for sample overall correlation coefficients at four different levels of depth were greater than, except for m at. correlation coefficients ranged from to for the 1st quartile, to for the 2nd quartile, to for the 3rd quartile, and to for the 4th quartile. the same pattern in terms of correlation coefficients changes at different levels of depth between sample and sample was also observed. discussion in the current study, rna-seq was performed twice using two chicken cdna samples. the first run of rna-seq had fewer number of reads and larger variation in terms of total number of reads between the two samples, while the second run had greater number of reads and very small variation between the two samples. the first run was performed at the very early stage of the sequencing technology when it was still in the testing phase. the lower reads and larger variation in the first run may be coming from two major sources of technical error: the purification of cdna templates during the library preparation, and the loading of libraries onto flow cells. these potential sources of errors were corrected during the second rna-seq analysis, which provided very good sequencing depth with greater number of reads. the first rna-seq datasets were directly derived from actual experiment, which made the results more informative than replicating datasets by random sampling. therefore, we chose to include these two early datasets in the analysis in the current study. furthermore, all of the reads from each sample were normalized by the rpkm and the datasets can serve as a reference for random sampling at different sequencing depths from the exact same samples. the capacity of sequencing length of bp for the first run was increased to bp at the second rna-seq analysis. longer reads should reduce estimation error and mapping uncertainty, and read lengths have been consistently increasing with improving illumina massively parallel sequencing technology. however, people have noted that the number of reads is more important than the read length once reaching a minimum read length of bp. the read lengths in the current study were larger than bp; therefore, the read length will not affect the overall conclusions drawn. as a powerful new technology for transcriptome analysis, rna-seq provides a more comprehensive view of the transcriptome than earlier technologies. besides its ability to detect splicing variation, rna editing and discovery of new transcripts, rna-seq can also function in the role of a conventional microarray in measuring gene expression due to its accurate measurements. in order to detect less abundant transcripts, appropriate sequencing depth is needed. the transcriptome coverage or sequencing depth needed for a given study can be affected by several factors such as genome size, transcriptome complexity and objectives of the study. in general, the more complex the transcriptome, the more sequencing depth is required for adequate coverage. depending on the purpose of transcriptome analysis, the requirement of sequencing depth varies. in most transcriptome studies, quantifying mrna abundance is one of the major objectives. there is a certain sequencing depth that is sufficient in simple transcriptomes. for example, in the yeast genome, a m reads dataset was generated by rna-seq which was able to get 100% transcriptome coverage. the number of transcripts detected by rna-seq in the yeast dataset was able to reach 80% transcriptome coverage at 4m mapped reads, and even though the sequencing depth doubled as 8m reads, the transcriptome coverage only increased 10%. these results suggest that the improvement of sequencing depth or transcriptome coverage after reaching a certain sequencing depth had relatively less impact on detecting low abundant genes. in addition, the cost per sample per lane by rna-seq is still not affordable for most laboratories. recent development in multiplex labelling using bar-coded libraries by illumina and continued increase in sequence output have made it possible to sequence multiple samples per lane without extra cost or running time. therefore, it is imperative to examine the correlation between sequencing depth and transcriptome coverage; in other words, what sequencing depth might be sufficient in reaching a certain level of transcriptome coverage and reliable measurement for rna-seq. in order to accomplish this objective, two approaches could be applied: experimental or simulation methods. both methods have been applied in this study. high correlation among replicates within each sequencing depth, gradual increase in correlation coefficients from m to m, and consistent patterns observed between samples and have demonstrated that random sampling was an effective and reliable method in reaching the goals of this study. transcriptome coverage is one of the most important parameters in profiling global gene expression. the number and level of transcript isoforms is not always known and transcription activity varies across the genome. this was confirmed in a study by using the number of unique transcription start sites as a measure of coverage in mouse cells. in the current study, we took a more practical approach using all annotated genes in the chicken genome. because the chicken genome is far under-annotated, we assume that the annotated chicken genes in the database would well represent different levels of expression abundance in the chicken genome, which is essential for the analysis of transcriptome coverage in this study. since gene expression depends on tissue and time of biological process, it is impossible for any single tissue to have all genes in the genome expressed. ninety percent of all annotated genes detected at about m reads might represent a saturated detection of the whole genome. the analysis results indicated significant improvements of transcriptome coverage occurred from m to m and from m to m. depending on the purpose of transcriptome analysis, the current study suggested that m reads could have 80% of transcriptome coverage, while m reads could reach 90% of transcriptome coverage. when we analyzed overall correlation coefficients at different levels of sequencing depth regardless of gene expression level, we observed very high correlation coefficients between each level of sequencing depth compared with m, except for m. one might draw a conclusion that there is no significant difference among different levels of sequencing depth. but as we see in figure this might be true in the case of highly abundant genes, but not in the case of the 1st to 3rd quartile groups, especially the first two quartile groups. high abundant genes will be less affected by sequencing depth than low abundant genes, because high abundant genes are more likely to be captured even when the sequencing depth varies. this is also confirmed by our analysis. collectively, the following points can be inferred: 1) sequencing depth below m reads had a significant effect on detecting transcript levels of medium and low abundant transcripts; 2) sequencing depth at both m and m would result in unreliable mrna expression on all genes except highly abundant transcripts; 3) there was no significant improvement in correlation coefficients when sequencing depth doubled from m to m. based on these analysis, the results suggested: 1) m reads might be sufficient in obtaining reliable mrna expression measurement on highly abundant transcripts; 2) when sequencing depth is beyond m reads, a relatively reliable measurement of mrna expression is expected, especially for abundant transcripts; 3) it seems that m of reads is needed to achieve reliable measurement of mrna expression across all genes in the chicken genome. to our knowledge, this is the first study evaluating the appropriate sequencing depth using rna-seq in farm animals and will provide the first reference for similar studies. the knowledge generated from this study has laid a solid foundation for applying this analysis to other species.  in summary, the analysis from this study demonstrated that m reads is sufficient to detect all annotated genes and provide a reliable measurement of mrna abundance in chicken lungs using rna-seq. as we expected, the depth of sequencing had significant impact on low abundance transcripts, but not on high abundance transcripts. in practice, if rna-seq serves as a replacement of microarray technology, m reads would allow detection of about 80% of annotated chicken genes. finally, increasing the depth of sequencing from m to m reads did not have a significant effect on transcriptome coverage and reliability of mrna measurements, whereas m reads was needed to achieve reliable measurement of mrna expression across all genes in the genome. competing interests the authors declare that there are no competing interests. authors' contributions yw analyzed the data and drafted the manuscript, ng wrote the code for sampling the rna-seq data, cj and ub contributed ideas for the data analysis and revised the manuscript, hw and rc ran rna-seq samples, and hz designed the experiment, provided the concept of the analysis. all authors submitted comments, read and approved the final manuscript. 
to evaluate the appropriate depth needed for mrna profiling, a random sampling method was used to generate different number of reads from each sample. the objective of this analysis is to evaluate what sequencing depth might be sufficient to interrogate gene expression profiling in the chicken by rna-seq. with continued improvement in sequence output, it would be cost-effective if multiple samples are multiplexed and sequenced in a single lane with sufficient transcriptome coverage. the correlation coefficients of gene expression between technical replicates within the same sample were and. there was a significant increase in correlation coefficients from a sequencing depth of m to m for all genes except highly abundant genes. rna-seq is the recently developed high-throughput sequencing technology for profiling the entire transcriptome in any organism. with significant improvements in sequencing technology, two technical replicate cdna libraries were re-sequenced. no significant improvement was observed from the depth of m to m reads. furthermore, the depth of sequencing had a significant impact on measuring gene expression of low abundant genes. two cdna libraries from chicken lungs were sequenced initially, and million and m reads were generated, respectively. the analysis from the current study demonstrated that m reads is sufficient to detect all annotated genes in chicken lungs. ten million reads could detect about 80% of annotated chicken genes, and rna-seq at this depth can serve as a replacement of microarray technology. finally, the combination of experimental and simulation approaches is a powerful approach to address the relationship between the depth of sequencing and transcriptome coverage.
 rna-seq is the recently developed high-throughput sequencing technology for profiling the entire transcriptome in any organism. it has several major advantages over current hybridization-based approach such as microarrays. however, the cost per sample by rna-seq is still prohibitive for most laboratories. with continued improvement in sequence output, it would be cost-effective if multiple samples are multiplexed and sequenced in a single lane with sufficient transcriptome coverage. the objective of this analysis is to evaluate what sequencing depth might be sufficient to interrogate gene expression profiling in the chicken by rna-seq.  two cdna libraries from chicken lungs were sequenced initially, and million and m reads were generated, respectively. with significant improvements in sequencing technology, two technical replicate cdna libraries were re-sequenced. totals of m and m reads were obtained with the two samples. more than 90% of annotated genes were detected in the data sets with - m reads, while only 68% of genes were detected in the data set with m reads. the correlation coefficients of gene expression between technical replicates within the same sample were and. to evaluate the appropriate depth needed for mrna profiling, a random sampling method was used to generate different number of reads from each sample. there was a significant increase in correlation coefficients from a sequencing depth of m to m for all genes except highly abundant genes. no significant improvement was observed from the depth of m to m reads.  the analysis from the current study demonstrated that m reads is sufficient to detect all annotated genes in chicken lungs. ten million reads could detect about 80% of annotated chicken genes, and rna-seq at this depth can serve as a replacement of microarray technology. furthermore, the depth of sequencing had a significant impact on measuring gene expression of low abundant genes. finally, the combination of experimental and simulation approaches is a powerful approach to address the relationship between the depth of sequencing and transcriptome coverage. 1- april eighth annual mcbios conference. computational biology and bioinformatics for a new decade college station, tx, usa 