title: Synthetic spike-in standards for RNA-seq experiments.
abstract: High-throughput sequencing of cDNA (RNA-seq) is a widely deployed transcriptome profiling and annotation technique, but questions about the performance of different protocols and platforms remain. We used a newly developed pool of 96 synthetic RNAs with various lengths, and GC content covering a 2(20) concentration range as spike-in controls to measure sensitivity, accuracy, and biases in RNA-seq experiments as well as to derive standard curves for quantifying the abundance of transcripts. We observed linearity between read density and RNA input over the entire detection range and excellent agreement between replicates, but we observed significantly larger imprecision than expected under pure Poisson sampling errors. We use the control RNAs to directly measure reproducible protocol-dependent biases due to GC content and transcript length as well as stereotypic heterogeneity in coverage across transcripts correlated with position relative to RNA termini and priming sequence bias. These effects lead to biased quantification for short transcripts and individual exons, which is a serious problem for measurements of isoform abundances, but that can partially be corrected using appropriate models of bias. By using the control RNAs, we derive limits for the discovery and detection of rare transcripts in RNA-seq experiments. By using data collected as part of the model organism and human Encyclopedia of DNA Elements projects (ENCODE and modENCODE), we demonstrate that external RNA controls are a useful resource for evaluating sensitivity and accuracy of RNA-seq experiments for transcriptome discovery and quantification. These quality metrics facilitate comparable analysis across different samples, protocols, and platforms.
authors.affiliation: Section of Developmental Genomics, Laboratory of Cellular and Developmental Biology, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD 20892, USA.
authors.forename: Lichun
authors.initials: L
authors.lastname: Jiang
authors.forename: Felix
authors.initials: F
authors.lastname: Schlesinger
authors.forename: Carrie A
authors.initials: CA
authors.lastname: Davis
authors.forename: Yu
authors.initials: Y
authors.lastname: Zhang
authors.forename: Renhua
authors.initials: R
authors.lastname: Li
authors.forename: Marc
authors.initials: M
authors.lastname: Salit
authors.forename: Thomas R
authors.initials: TR
authors.lastname: Gingeras
authors.forename: Brian
authors.initials: B
authors.lastname: Oliver
mesh_terms: D000818:Animals; D015982:Bias; D020869:Gene Expression Profiling; D015723:Gene Library; D059014:High-Throughput Nucleotide Sequencing / Q000592:standards; D006801:Humans; D011786:Quality Control; D012313:RNA / Q000737:chemistry*; D015203:Reproducibility of Results; D012680:Sensitivity and Specificity; D017423:Sequence Analysis, RNA / Q000592:standards*
chemical_list: D012313:RNA