title: Modelling and simulating generic RNA-Seq experiments with the flux simulator.
abstract: High-throughput sequencing of cDNA libraries constructed from cellular RNA complements (RNA-Seq) naturally provides a digital quantitative measurement for every expressed RNA molecule. Nature, impact and mutual interference of biases in different experimental setups are, however, still poorly understood-mostly due to the lack of data from intermediate protocol steps. We analysed multiple RNA-Seq experiments, involving different sample preparation protocols and sequencing platforms: we broke them down into their common--and currently indispensable--technical components (reverse transcription, fragmentation, adapter ligation, PCR amplification, gel segregation and sequencing), investigating how such different steps influence abundance and distribution of the sequenced reads. For each of those steps, we developed universally applicable models, which can be parameterised by empirical attributes of any experimental protocol. Our models are implemented in a computer simulation pipeline called the Flux Simulator, and we show that read distributions generated by different combinations of these models reproduce well corresponding evidence obtained from the corresponding experimental setups. We further demonstrate that our in silico RNA-Seq provides insights about hidden precursors that determine the final configuration of reads along gene bodies; enhancing or compensatory effects that explain apparently controversial observations can be observed. Moreover, our simulations identify hitherto unreported sources of systematic bias from RNA hydrolysis, a fragmentation technique currently employed by most RNA-Seq protocols.
authors.affiliation: Bioinformatics and Genomics Program, Centre de Regulació Genòmica (CRG), 08003 Barcelona, Spain.
authors.forename: Thasso
authors.initials: T
authors.lastname: Griebel
authors.forename: Benedikt
authors.initials: B
authors.lastname: Zacher
authors.forename: Paolo
authors.initials: P
authors.lastname: Ribeca
authors.forename: Emanuele
authors.initials: E
authors.lastname: Raineri
authors.forename: Vincent
authors.initials: V
authors.lastname: Lacroix
authors.forename: Roderic
authors.initials: R
authors.lastname: Guigó
authors.forename: Michael
authors.initials: M
authors.lastname: Sammeth
mesh_terms: D003198:Computer Simulation; D020869:Gene Expression Profiling; D059014:High-Throughput Nucleotide Sequencing; D006868:Hydrolysis; D012313:RNA / Q000378:metabolism; D017423:Sequence Analysis, RNA
chemical_list: D012313:RNA