Abstract Background High throughput DNA sequencing technology has enabled quantification of all the RNAs in a cell or tissue, a method widely known as RNA sequencing (RNA-Seq). However, non-coding RNAs such as rRNA are highly abundant and can consume >70% of sequencing reads. A common approach is to extract only polyadenylated mRNA; however, such approaches are blind to RNAs with short or no poly(A) tails, leading to an incomplete view of the transcriptome. Another challenge of preparing RNA-Seq libraries is to preserve the strand information of the RNAs. Design Here, we describe a procedure for preparing RNA-Seq libraries from 1 to 4 μg total RNA without poly(A) selection. Our method combines the deoxyuridine triphosphate (dUTP)/uracil-DNA glycosylase (UDG) strategy to achieve strand specificity with AMPure XP magnetic beads to perform size selection. Together, these steps eliminate gel purification, allowing a library to be made in less than two days. We barcode each library during the final PCR amplification step, allowing several samples to be sequenced in a single lane without sacrificing read length. Libraries prepared using this protocol are compatible with Illumina GAII, GAIIx and HiSeq 2000 platforms. Discussion The RNA-Seq protocol described here yields strand-specific transcriptome libraries without poly(A) selection, which provide approximately 90% mappable sequences. Typically, more than 85% of mapped reads correspond to protein-coding genes and only 6% derive from non-coding RNAs. The protocol has been used to measure RNA transcript identity and abundance in tissues from flies, mice, rats, chickens, and frogs, demonstrating its general applicability.

Long-read RNA sequencing has arisen as a counterpart to short-read sequencing, with the potential to capture full-length isoforms, albeit at the cost of lower depth. Yet this potential is not fully realized due to inherent limitations of current long-read assembly methods and underdeveloped approaches to integrate short-read data. Here, we critically compare the existing methods and develop a new integrative approach to characterize a particularly challenging pool of low-abundance long noncoding RNA (lncRNA) transcripts from short- and long-read sequencing in two distinct cell lines. Our analysis reveals severe limitations in each of the sequencing platforms. For short-read assemblies, coverage declines at transcript termini resulting in ambiguous ends, and uneven low coverage results in segmentation of a single transcript into multiple transcripts. Conversely, long-read sequencing libraries lack depth and strand-of-origin information in cDNA-based methods, culminating in erroneous assembly and quantitation of transcripts. We also discover a cDNA synthesis artifact in long-read datasets that markedly impacts the identity and quantitation of assembled transcripts. Towards remediating these problems, we develop a computational pipeline to "strand" long-read cDNA libraries that rectifies inaccurate mapping and assembly of long-read transcripts. Leveraging the strengths of each platform and our computational stranding, we also present and benchmark a hybrid assembly approach that drastically increases the sensitivity and accuracy of full-length transcript assembly on the correct strand and improves detection of biological features of the transcriptome. When applied to a challenging set of under-annotated and cell-type variable lncRNA, our method resolves the segmentation problem of short-read sequencing and the depth problem of long-read sequencing, resulting in the assembly of coherent transcripts with precise 5' and 3' ends. Our workflow can be applied to existing datasets for superior demarcation of transcript ends and refined isoform structure, which can enable better differential gene expression analyses and molecular manipulations of transcripts.

Tn5 transposon tagments double-stranded DNA and RNA/DNA hybrids to generate nucleic acids that are ready to be amplified for high-throughput sequencing. The nucleic acid substrates for the Tn5 transposon must be explored to increase the applications of Tn5. Here, we found that the Tn5 transposon can transpose oligos into the 5′ end of single-stranded DNA longer than 140 nucleotides. Based on this property of Tn5, we developed a tagmentation-based and ligation-enabled single-stranded DNA sequencing method called TABLE-seq. Through a series of reaction temperature, time, and enzyme concentration tests, we applied TABLE-seq to strand-specific RNA sequencing, starting with as little as 30 pg of total RNA. Moreover, compared with traditional dUTP-based strand-specific RNA sequencing, this method detects more genes, has a higher strand specificity, and shows more evenly distributed reads across genes. Together, our results provide insights into the properties of Tn5 transposons and expand the applications of Tn5 in cutting-edge sequencing techniques.

This unit presents protocols for construction of next-generation sequencing (NGS) directional RNA sequencing libraries for the Illumina HiSeq and MiSeq from a wide variety of input RNA sources. The protocols are based on the New England Biolabs (NEB) small RNA library preparation set for Illumina, although similar kits exist from different vendors. The protocol preserves the orientation of the original RNA in the final sequencing library, enabling strand-specific analysis of the resulting data. These libraries have been used for differential gene expression analysis and small RNA discovery and are currently being tested for de novo transcriptome assembly. The protocol is robust and applicable to a broad range of RNA input types and RNA quality, making it ideal for high-throughput laboratories.