On the optimal design of metabolic RNA labeling experiments

Massively parallel RNA sequencing (RNA-seq) in combination with metabolic labeling has become the de facto standard approach to study alterations in RNA transcription, processing or decay. Regardless of advances in the experimental protocols and techniques, every experimentalist needs to specify the key aspects of experimental design: For example, which protocol should be used (biochemical separation vs. nucleotide conversion) and what is the optimal labeling time? In this work, we provide approximate answers to these questions using asymptotic theory of optimal design. Specifically, we derive the optimal labeling time for any given degradation rate and show that sub-optimal time points yield better rate estimates if they precede the optimal time point. Subsequently, we show that an increase in sample numbers should be preferred over an increase in sequencing depth. Lastly, we provide some guidance on use cases when laborious biochemical separation outcompetes recent nucleotide conversion based methods (such as SLAMseq).

[1]  M. Selbach,et al.  Global quantification of mammalian gene expression control , 2011, Nature.

[2]  Vladimir Benes,et al.  Translational control mediated by hnRNP K links NMHC IIA to erythroid enucleation , 2016, Journal of Cell Science.

[3]  J. Boothroyd,et al.  Biosynthetic labeling of RNA with uracil phosphoribosyltransferase allows cell-specific microarray analysis of mRNA synthesis and decay , 2005, Nature Biotechnology.

[4]  Meaghan C. Sullivan,et al.  TimeLapse-seq: Adding a temporal dimension to RNA sequencing through nucleoside recoding , 2018, Nature Methods.

[5]  Julien Gagneur,et al.  Measures of RNA metabolism rates: Toward a definition at the level of single bonds , 2016, Transcription.

[6]  L. Huang,et al.  Regulation of hypoxia-inducible factor 1α is mediated by an O2-dependent degradation domain via the ubiquitin-proteasome pathway , 1998 .

[7]  Florian Erhard,et al.  Dissecting newly transcribed and old RNA using GRAND-SLAM , 2018, Bioinform..

[8]  Matthias W. Hentze,et al.  Drosophila miR2 induces pseudo-polysomes and inhibits translation initiation , 2007, Nature.

[9]  Wolfgang Huber,et al.  Differential expression of RNA-Seq data at the gene level – the DESeq package , 2012 .

[10]  A. Bos.  Parameter Estimation for Scientists and Engineers , 2007 .

[11]  Dietmar Rieder,et al.  Osmium-Mediated Transformation of 4-Thiouridine to Cytidine as Key To Study RNA Dynamics by Sequencing. , 2017, Angewandte Chemie.

[12]  Christoph Dieterich,et al.  pulseR: Versatile computational analysis of RNA turnover from metabolic labeling experiments , 2017, Bioinform..

[13]  Susan R. Wilson,et al.  Efficient experimental design and analysis strategies for the detection of differential expression using RNA-Sequencing , 2012, BMC Genomics.

[14]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[15]  Johannes Zuber,et al.  Thiol-linked alkylation of RNA to assess expression dynamics , 2017, Nature Methods.

[16]  R. Zimmer,et al.  High-resolution gene expression profiling for simultaneous kinetic parameter analysis of RNA synthesis and decay. , 2008, RNA.

[17]  Y. Pawitan.  In all likelihood : statistical modelling and inference using likelihood , 2002 .

[18]  Gernot Marx,et al.  Translation control of TAK1 mRNA by hnRNP K modulates LPS-induced macrophage activation , 2014, RNA.

[19]  S. Salzberg,et al.  StringTie enables improved reconstruction of a transcriptome from RNA-seq reads , 2015, Nature Biotechnology.

[20]  J. Gagneur,et al.  TT-seq maps the human transient transcriptome , 2016, Science.

[21]  G. Fonarow,et al.  Epidemiology and aetiology of heart failure , 2016, Nature Reviews Cardiology.

[22]  Nicholas Denko,et al.  Hypoxia induces a novel signature of chromatin modifications and global repression of transcription. , 2008, Mutation research.

[23]  Mark B Gerstein,et al.  Tracking Distinct RNA Populations Using Efficient and Reversible Covalent Chemistry. , 2015, Molecular cell.

[24]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[25]  Achim Tresch,et al.  Dynamic transcriptome analysis measures rates of mRNA synthesis and decay in yeast , 2011, Molecular systems biology.

[26]  G. Semenza.  Targeting HIF-1 for cancer therapy , 2003, Nature Reviews Cancer.

[27]  Mircea Ivan,et al.  Post-Transcriptional Control of the Hypoxic Response by RNA-Binding Proteins and MicroRNAs , 2011, Front. Mol. Neurosci..

[28]  Julien Gagneur,et al.  Determinants of RNA metabolism in the Schizosaccharomyces pombe genome , 2015, bioRxiv.

[29]  Steven P Lund,et al.  Statistical Applications in Genetics and Molecular Biology Detecting Differential Expression in RNA-sequence Data Using Quasi-likelihood with Shrunken Dispersion Estimates , 2012 .

[30]  H. Chernoff.  Locally Optimal Designs for Estimating Parameters , 1953 .

[31]  W. Wilson,et al.  Targeting hypoxia in cancer therapy , 2011, Nature Reviews Cancer.

[32]  Lars Dölken,et al.  RNA dynamics revealed by metabolic RNA labeling and biochemical nucleoside conversions , 2018, Nature Methods.

[33]  Knut Reinert,et al.  Flexbar 3.0 ‐ SIMD and multicore parallelization , 2017, Bioinform..