Abstract:BackgroundRNA sequencing (RNA-Seq) is often used for transcriptome profiling as well as the identification of novel transcripts and alternative splicing events. Typically, RNA-Seq libraries are prepared from total RNA using poly(A) enrichment of the mRNA (mRNA-Seq) to remove ribosomal RNA (rRNA), however, this method fails to capture non-poly(A) transcripts or partially degraded mRNAs. Hence, a mRNA-Seq protocol will not be compatible for use with RNAs coming from Formalin-Fixed and Paraffin-Embedded (FFPE) samples.ResultsTo address the desire to perform RNA-Seq on FFPE materials, we evaluated two different library preparation protocols that could be compatible for use with small RNA fragments. We obtained paired Fresh Frozen (FF) and FFPE RNAs from multiple tumors and subjected these to different gene expression profiling methods. We tested 11 human breast tumor samples using: (a) FF RNAs by microarray, mRNA-Seq, Ribo-Zero-Seq and DSN-Seq (Duplex-Specific Nuclease) and (b) FFPE RNAs by Ribo-Zero-Seq and DSN-Seq. We also performed these different RNA-Seq protocols using 10 TCGA tumors as a validation set.The data from paired RNA samples showed high concordance in transcript quantification across all protocols and between FF and FFPE RNAs. In both FF and FFPE, Ribo-Zero-Seq removed rRNA with comparable efficiency as mRNA-Seq, and it provided an equivalent or less biased coverage on gene 3′ ends. Compared to mRNA-Seq where 69% of bases were mapped to the transcriptome, DSN-Seq and Ribo-Zero-Seq contained significantly fewer reads mapping to the transcriptome (20-30%); in these RNA-Seq protocols, many if not most reads mapped to intronic regions. Approximately 14 million reads in mRNA-Seq and 45–65 million reads in Ribo-Zero-Seq or DSN-Seq were required to achieve the same gene detection levels as a standard Agilent DNA microarray.ConclusionsOur results demonstrate that compared to mRNA-Seq and microarrays, Ribo-Zero-Seq provides equivalent rRNA removal efficiency, coverage uniformity, genome-based mapped reads, and consistently high quality quantification of transcripts. Moreover, Ribo-Zero-Seq and DSN-Seq have consistent transcript quantification using FFPE RNAs, suggesting that RNA-Seq can be used with FFPE-derived RNAs for gene expression profiling.

Gene expression profiling of non-polyadenylated RNA-seq across species.

Genomewide Characterization of Non-Polyadenylated RNAs.

Fractionation of non-polyadenylated and ribosomal-free RNAs from mammalian cells.

Novel Method of Full-Length RNA-seq That Expands the Identification of Non-Polyadenylated RNAs Using Nanopore Sequencing

Identification and analysis of mouse non-coding RNA using transcriptome data

Systematic identification and evolutionary features of rhesus monkey small nucleolar RNAs

RAG-seq: A NSR Primed and Transposase Tagmentation Mediated Strand-specific Total RNA Sequencing in Single Cell

Nano3P-seq: transcriptome-wide analysis of gene expression and tail dynamics using end-capture nanopore cDNA sequencing

Using Pan RNA-Seq Analysis to Reveal the Ubiquitous Existence of 5' and 3' End Small RNAs.

Accurate Profiling of Gene Expression and Alternative Polyadenylation with Whole Transcriptome Termini Site Sequencing (WTTS-Seq)

[Studies on the transcriptomes of non-model organisms]

Novel Rrna-Depletion Methods for Total RNA Sequencing and Ribosome Profiling Developed for Avian Species.

Poly(A)-seq: A Method for Direct Sequencing and Analysis of the Transcriptomic Poly(a)-Tails.

Gene Expression Analysis Using RNA-Seq from Organisms Lacking Substantial Genomic Resources

Digital Genome-Wide ncRNA Expression, Including SnoRNAs, across 11 Human Tissues Using PolyA-Neutral Amplification

Deep annotation of long noncoding RNAs by assembling RNA-seq and small RNA-seq data

Using pan RNA-seq analysis to reveal the ubiquitous existence of 5’ end and 3’ end small RNAs

Comparison of RNA-Seq by poly (A) capture, ribosomal RNA depletion, and DNA microarray for expression profiling

Nanopore-based native RNA sequencing of human transcriptomes reveals the complexity of mRNA modifications and crosstalk between RNA regulatory features

Polyadenylation-related Isoform Switching in Human Evolution Revealed by Full-Length Transcript Structure

An RNA tagging approach for system-wide RNA-binding proteome profiling and dynamics investigation upon transcription inhibition