Using normalization to resolve RNA-seq biases caused by amplification from minimal input

Physiological Genomics

Volumn 46, Issue 21, 2014, Pages 808-820

  Ager Wick, Eirill ^a   Henkel, Christiaan V ^a,b   Haug, Trude M ^c   Weltzien, Finn Arne ^a,c  

^a NORWEGIAN UNIVERSITY OF LIFE SCIENCES   (Norway)

^b LEIDEN UNIVERSITY   (Netherlands)

^c UNIVERSITY OF OSLO   (Norway)

Author keywords

low RNA input; Medaka; Normalization; Pituitary; RNA Seq

Indexed keywords

RNA; TRANSCRIPTOME; GREEN FLUORESCENT PROTEIN;

ANIMAL EXPERIMENT; ARTICLE; CELL POPULATION; CELL SELECTION; DNA SEQUENCE; FEMALE; GENE AMPLIFICATION; GENE EXPRESSION; HYPOPHYSIS CELL; MALE; NONHUMAN; ORYZIAS; RNA ISOLATION; RNA SEQUENCE; ANIMAL; CELL CULTURE; CYTOLOGY; GENETICS; HYPOPHYSIS; ISOLATION AND PURIFICATION; PROCEDURES; REPRODUCIBILITY; SEQUENCE ANALYSIS; TRANSGENIC ANIMAL;

ANIMALS; ANIMALS, GENETICALLY MODIFIED; CELLS, CULTURED; FEMALE; GENE EXPRESSION; GREEN FLUORESCENT PROTEINS; ORYZIAS; PITUITARY GLAND; REPRODUCIBILITY OF RESULTS; RNA; SEQUENCE ANALYSIS, RNA;

EID: 84908459005     PISSN: 10948341     EISSN: 15312267     Source Type: Journal    
DOI: 10.1152/physiolgenomics.00196.2013     Document Type: Article

References (38)

1. Adiconis X, Borges-Rivera D, Satija R, DeLuca DS, Busby MA, Berlin AM, Sivachenko A, Thompson DA, Wysoker A, Fennell T, Gnirke A, Pochet N, Regev A, Levin JZ. Comparative analysis of RNA sequencing methods for degraded or low-input samples. Nat Methods 10: 623-629, 2013.
2. Ager-Wick E, Dirks RP, Burgerhout E, Nourizadeh-Lillabadi R, de Wijze DL, Spaink HP, van den Thillart GEEJM, Tsukamoto K, Dufour S, Weltzien FA, Henkel C.V. The pituitary gland of the European eel reveals massive expression of genes involved in the melanocotin system. PLoS One 8: e77396, 2013.
3. Bhargava V, Head SR, Ordoukhanian P, Mercola M, Subramaniam S. Technical variations in low-input RNA-Seq methodologies. Sci Rep 4: 3678, 2014.
4. Bolstad BM, Irizarry RA, Astrand M, Speed TP. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics 19: 185-193, 2003.
5. Bullard JH, Purdom E, Hansen KD, Dudoit S. Evaluation of statistical methods for normalization and differential expression in mRNA-Seq experiments. BMC Bioinformatics 11: 94, 2010.
6. Cloonan N, Forrest AR, Kolle G, Gardiner BB, Faulkner GJ, Brown MK, Taylor DF, Steptoe AL, Wani S, Bethel G, Robertson AJ, Perkins AC, Bruce SJ, Lee CC, Ranade SS, Peckham HE, Manning JM, McKernan KJ, Grimmond SM. Stem cell transcriptome profiling via massive-scale mRNA sequencing. Nat Methods 5: 613-619, 2008.
7. Efron B, Tibshirani R, Storey JD, Tusher V. Empirical Bayes analysis of microarray experiment. J Am Stat Assoc 96: 1151-1160, 2001.
8. Hansen KD, Irizarry RA, Wu Z. Removing technical variability in RNA-Seq data using conditional quantile normalization. Biostatistics 13: 204-216, 2012.
9. Hildahl J, Sandvik GK, Lifjeld R, Hodne K, Nagahama Y, Haug TM, Okubo K, Weltzien FA. Developmental tracing of luteinizing hormone beta-subunit gene expression using green fluorescent protein transgenic medaka (Oryzias latipes) reveals a putative novel developmental function. Dev Dyn 241: 1665-1677, 2012.
10. Hodne K, von Krogh K, Weltzien FA, Sand O, Haug TM. Optimized conditions for primary culture of pituitary cells from the Atlantic cod (Gadus morhua). The importance of osmolality, pCO2, and pH. Gen Comp Endocrinol 178: 206-215, 2012.
11. Kanda S, Okubo K, Oka Y. Differential regulation of the luteinizing hormone genes in teleosts and tetrapods due to their distinct genomic environments-Insights into gonadotropin beta subunit evolution. Gen Comp Endocrinol 173: 253-258, 2011.
12. Kasahara M, Naruse K, Sasaki S, Nakatani Y, Qu W, Ahsan B, Yamada T, Nagayasu Y, Doi K, Kasai Y, Jindo T, Kobayashi D, Shimada A, Toyoda A, Kuroki Y, Fujiyama A, Sasaki T, Shimizu A, Asakawa S, Shimizu N, Hashimoto S, Yang J, Lee Y, Matsushima K, Sugano S, Sakaizumi M, Narita T, Ohishi K, Haga S, Ohta F, Nomoto H, Nogata K, Morishita T, Endo T, Shin I, Takeda H, Morishita S, Kohara Y. The medaka draft genome and insights into vertebrate genome evolution. Nature 447: 714-719, 2007.
13. Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol 14: R36, 2013.
14. Kinoshita M, Murata K, Naruse K, Tanaka M. Medaka: Biology, Management, and Experimental Protocols. Wiley-Blackwell, 2009.
15. Kurn N, Chen PC, Heath JD, Kopf-Sill A, Stephens KM, Wang SL. Novel isothermal, linear nucleic acid amplification systems for highly multiplexed applications. Clin Chem 51: 1973-1981, 2005.
16. Li B, Dewey CN. RSEM: accurate transcript quantification from RNASeq data with or without a reference genome. BMC Bioinformatics 12: 323, 2011.
17. Li B, Ruotti V, Stewart R., Thomson JA, Dewey CN. RNA-Seq gene expression estimation with read mapping uncertainty. Bioinformatics 26: 493-500, 2010.
18. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25: 2078-2079, 2009.
19. Malboeuf CM, Yang X, Charlebois P, Qu J, Berlin AM, Casali M, Pesko KN, Boutwell CL, DeVincenzo JP, Ebel GD, Allen TM, Zody MC, Henn MR, Levin JZ. Complete viral RNA genome sequencing of ultra-low copy samples by sequence-independent amplification. Nucleic Acids Res 41: e13, 2013.
20. Marioni JC, Mason CE, Mane SM, Stephens M, Gilad Y. RNA-Seq: an assessment of technical reproducibility and comparison with gene expression arrays. Genome Res 18: 1509-1517, 2008.
21. Matsuda M, Nagahama Y, Shinomiya A, Sato T, Matsuda C, Kobayashi T, Morrey CE, Shibata N, Asakawa S, Shimizu N, Hori H, Hamaguchi S, Sakaizumi M. DMY is a Y-specific DM-domain gene required for male development in the medaka fish. Nature 417: 559-563, 2002.
22. Mortazavi A, Williams BA, Mccue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods 5: 621-628, 2008.
23. Nagalakshmi U, Wang Z, Waern K, Shou C, Raha D, Gerstein M, Snyder M. The transcriptional landscape of the yeast genome defined by RNA sequencing. Science 320: 1344-1349, 2008.
24. Nanda I, Kondo M, Hornung U, Asakawa S, Winkler C, Shimizu A, Shan ZH, Haaf T, Shimizu N, Shima A, Schmid M, Schartl M. A duplicated copy of DMRT1 in the sex-determining region of the Y chromosome of the medaka, Oryzias latipes. Proc Natl Acad Sci USA 99: 11778-11783, 2002.
25. Oshlack A, Robinson MD, Young MD. From RNA-Seq reads to differential expression results. Genome Biol 11: 220, 2010.
26. Pan XH, Durrett RE, Zhu HY, Tanaka Y, Li YM, Zi XY, Marjani SL, Euskirchen G, Ma C, LaMotte RH, Park IH, Snyder M., Mason CE, Weissman SM. Two methods for full-length RNA sequencing for low quantities of cells and single cells. Proc Natl Acad Sci USA 110: 594-599, 2013.
27. Robinson JT, Thorvaldsdottir H, Winckler W, Guttman M, Lander ES, Getz G, Mesirov JP. Integrative genomics viewer. Nat Biotechnol 29: 24-26, 2011.
28. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26: 139-140, 2010.
29. Roche Diagnostics GmbH. Transcriptome profiling of low-input RNA samples. http://454.com/downloads/NuGENTranscriptomeAppNote_Finalv3.pdf (02.05.2014).
30. Sengupta S, Ruotti V, Bolin J, Elwell A, Hernandez A, Thomson JA, Stewart R. Highly consistent, fully representative mRNA-Seq libraries from ten nanograms of total RNA. Biotechniques 49: 898-904, 2010.
31. Soneson C, Delorenzi M. A comparison of methods for differential expression analysis of RNA-Seq data. BMC Bioinformatics 14: 91, 2013.
32. Strandabø RAU, Hodne K, Ager-Wick E, Sand O, Weltzien FA Haug TM. Signal transduction involved in GnRH2-stimulation of identified LH-producing gonadotropes from lhb-GFP transgenic medaka (Oryzias latipes). Mol Cell Endocrinol 372: 128-139, 2013.
33. Tarazona S, Furio P, Ferrer A, Conesa A. NOISeq: an RNA-Seq package for differential expression in RNA-Seq using biological replicates. Institute of Computational Genomics, Príncipe Felipe Research Center, Spain. http://bioinfo.cipf.es/noiseq/lib/exe/fetch.php?medianoiseqbio_techreport.pdf (02.05.2014).
34. Tarazona S, Garcia-Alcalde F, Dopazo J, Ferrer A, Conesa A. Differential expression in RNA-Seq: a matter of depth. Genome Res 21: 2213-2223, 2011.
35. Tariq MA, Kim HJ, Jejelowo O, Pourmand N. Whole-transcriptome RNAseq analysis from minute amount of total RNA. Nucleic Acids Res 39: e120, 2011.
36. Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics 25: 1105-1111, 2009.
37. Vallandingham J, Fleharty B, Peak A, Staehling K, Perera A, Li H. Evaluation of whole transcriptome amplification methods by RNA-Seq. J Biomol Tech 24, Suppl: S57-S58, 2013.
38. Wickham H. ggplot2: Elegant Graphics for Data Analysis. Springer, 2009.