Elimination of PCR duplicates in RNA-seq and small RNA-seq using unique molecular identifiers

Abstract

BackgroundRNA-seq and small RNA-seq are powerful, quantitative tools to study gene regulation and function. Common high-throughput sequencing methods rely on polymerase chain reaction (PCR) to expand the starting material, but not every molecule amplifies equally, causing some to be overrepresented. Unique molecular identifiers (UMIs) can be used to distinguish undesirable PCR duplicates derived from a single molecule and identical but biologically meaningful reads from different molecules.ResultsWe have incorporated UMIs into RNA-seq and small RNA-seq protocols and developed tools to analyze the resulting data. Our UMIs contain stretches of random nucleotides whose lengths sufficiently capture diverse molecule species in both RNA-seq and small RNA-seq libraries generated from mouse testis. Our approach yields high-quality data while allowing unique tagging of all molecules in high-depth libraries.ConclusionsUsing simulated and real datasets, we demonstrate that our methods increase the reproducibility of RNA-seq and small RNA-seq data. Notably, we find that the amount of starting material and sequencing depth, but not the number of PCR cycles, determine PCR duplicate frequency. Finally, we show that computational removal of PCR duplicates based only on their mapping coordinates introduces substantial bias into data analysis.Electronic supplementary materialThe online version of this article (10.1186/s12864-018-4933-1) contains supplementary material, which is available to authorized users.

Co-dependent Assembly of Drosophila piRNA Precursor Complexes and piRNA Cluster Heterochromatin

“…Strand specific RNA-seq libraries were constructed as described previously (Zhang et al, 2012b) with modification in the rRNA depletion procedure using enzymatic digestion of rRNA by Hybridase™ Thermostable RNase H (Epicenter) with a comprehensive mixture of antisense rRNA oligos (Fu et al, 2018). The small RNaseq library is constructed as detailed previously (Li et al, 2009) with 2S rRNA depletion as described in (Zhang et al, 2011).…”

Depletion of TRRAP Induces p53‐Independent Senescence in Liver Cancer by Down‐Regulating Mitotic Genes

“…Ribosomal RNA was depleted and RNA‐sequencing (RNA‐seq) libraries were prepared as previously described . The libraries were quantified using the KAPA Library Quantification Kit (Kapa Biosystems, Wilmington, MA) and sequenced using the Illumina NextSeq 500 system.…”

Maelstrom Represses Canonical Polymerase II Transcription within Bi-directional piRNA Clusters in Drosophila melanogaster

“…Construction and Analysis of RNA-Seq Libraries RNA-seq libraries were constructed as described (Zhang et al, 2012b) with several modifications, including the use of unique molecular identifiers to eliminate PCR duplicates (Fu et al, 2018a). For ribosomal RNA depletion, RNA was hybridized in 10 ml with a pool of 186 rRNA antisense oligos (0.05 mM/each; Morlan et al, 2012;Adiconis et al, 2013) in 10 mM Tris-HCl [pH 7.4] and 20 mM NaCl and heated to 95 C, then cooled at À0.1 C/sec to 22 C, and finally incubated at 22 C for 5 min.…”