Evaluation of protocols for rRNA depletion-based RNA sequencing of nanogram inputs of mammalian total RNA

A scalable and cost-efficient rRNA depletion approach to enrich RNAs for molecular biology investigations

Transcriptomics analyses play pivotal roles in understanding the complex regulatory networks that govern cellular processes. The abundance of rRNAs, which account for 80%-90% of total RNA in eukaryotes, limits the detection and investigation of other transcripts. While mRNAs and long noncoding RNAs have poly(A) tails that are often used for positive selection, investigations of poly(A)- RNAs, such as circular RNAs, histone mRNAs, and small RNAs, typically require the removal of the abundant rRNAs for enrichment. Current approaches to deplete rRNAs for downstream molecular biology investigations are hampered by restrictive RNA input masses and high costs. To address these challenges, we developed rRNA Removal by RNaseH (rRRR), a method to efficiently deplete rRNAs from a wide range of human, mouse, and rat RNA inputs and of varying qualities at a cost 10- to 20-fold cheaper than other approaches. We used probe-based hybridization and enzymatic digestion to selectively target and remove rRNA molecules while preserving the integrity of non-rRNA transcripts. Comparison of rRRR to two commercially available approaches showed similar rRNA depletion efficiencies and comparable off-target effects. Our developed method provides researchers with a valuable tool for investigating gene expression and regulatory mechanisms across a wide range of biological systems at an affordable price that increases the accessibility for researchers to enter the field, ultimately advancing our understanding of cellular processes.

Challenges and opportunities to computationally deconvolve heterogeneous tissue with varying cell sizes using single-cell RNA-sequencing datasets

Deconvolution of cell mixtures in "bulk" transcriptomic samples from homogenate human tissue is important for understanding disease pathologies. However, several experimental and computational challenges impede transcriptomics-based deconvolution approaches using single-cell/nucleus RNA-seq reference atlases. Cells from the brain and blood have substantially different sizes, total mRNA, and transcriptional activities, and existing approaches may quantify total mRNA instead of cell type proportions. Further, standards are lacking for the use of cell reference atlases and integrative analyses of single-cell and spatial transcriptomics data. We discuss how to approach these key challenges with orthogonal "gold standard" datasets for evaluating deconvolution methods.

RNA-seq RNAaccess identified as the preferred method for gene expression analysis of low quality FFPE samples

Clinical tumor tissues that are preserved as formalin-fixed paraffin-embedded (FFPE) samples result in extensive cross-linking, fragmentation, and chemical modification of RNA, posing significant challenges for RNA-seq-based gene expression profiling. This study sought to define an optimal RNA-seq protocol for FFPE samples. We employed a common RNA extraction method and then compared RNA-seq library preparation protocols including RNAaccess, RiboZero and PolyA in terms of sequencing quality and concordance of gene expression using FFPE and case-matched fresh-frozen (FF) triple-negative breast cancer (TNBC) tissues. We found that RNAaccess, a method based on exome capture, produced the most concordant results. Applying RNAaccess to FFPE gastric cancer tissues, we established a minimum RNA DV200 requirement of 10% and a RNA input amount of 10ng that generated highly reproducible gene expression data. Lastly, we demonstrated that RNAaccess and NanoString platforms produced highly concordant expression profiles from FFPE samples for shared genes; however, RNA-seq may be preferred for clinical biomarker discovery work because of the broader coverage of the transcriptome. Taken together, these results support the selection of RNA-seq RNAaccess method for gene expression profiling of FFPE samples. The minimum requirements for RNA quality and input established here may allow for inclusion of clinical FFPE samples of sub-optimal quality in gene expression analyses and ultimately increasing the statistical power of such analyses.

Performance of computational algorithms to deconvolve heterogeneous bulk ovarian tumor tissue depends on experimental factors

BACKGROUND

Single-cell gene expression profiling provides unique opportunities to understand tumor heterogeneity and the tumor microenvironment. Because of cost and feasibility, profiling bulk tumors remains the primary population-scale analytical strategy. Many algorithms can deconvolve these tumors using single-cell profiles to infer their composition. While experimental choices do not change the true underlying composition of the tumor, they can affect the measurements produced by the assay.

RESULTS

We generated a dataset of high-grade serous ovarian tumors with paired expression profiles from using multiple strategies to examine the extent to which experimental factors impact the results of downstream tumor deconvolution methods. We find that pooling samples for single-cell sequencing and subsequent demultiplexing has a minimal effect. We identify dissociation-induced differences that affect cell composition, leading to changes that may compromise the assumptions underlying some deconvolution algorithms. We also observe differences across mRNA enrichment methods that introduce additional discrepancies between the two data types. We also find that experimental factors change cell composition estimates and that the impact differs by method.

CONCLUSIONS

Previous benchmarks of deconvolution methods have largely ignored experimental factors. We find that methods vary in their robustness to experimental factors. We provide recommendations for methods developers seeking to produce the next generation of deconvolution approaches and for scientists designing experiments using deconvolution to study tumor heterogeneity.

Nuclei on the Rise: When Nuclei-Based Methods Meet Next-Generation Sequencing

In the last decade, we have witnessed an upsurge in nuclei-based studies, particularly coupled with next-generation sequencing. Such studies aim at understanding the molecular states that exist in heterogeneous cell populations by applying increasingly more affordable sequencing approaches, in addition to optimized methodologies developed to isolate and select nuclei. Although these powerful new methods promise unprecedented insights, it is important to understand and critically consider the associated challenges. Here, we provide a comprehensive overview of the rise of nuclei-based studies and elaborate on their advantages and disadvantages, with a specific focus on their utility for transcriptomic sequencing analyses. Improved designs and appropriate use of the various experimental strategies will result in acquiring biologically accurate and meaningful information.

Library adaptors with integrated reference controls improve the accuracy and reliability of nanopore sequencing

Library adaptors are short oligonucleotides that are attached to RNA and DNA samples in preparation for next-generation sequencing (NGS). Adaptors can also include additional functional elements, such as sample indexes and unique molecular identifiers, to improve library analysis. Here, we describe Control Library Adaptors, termed CAPTORs, that measure the accuracy and reliability of NGS. CAPTORs can be integrated within the library preparation of RNA and DNA samples, and their encoded information is retrieved during sequencing. We show how CAPTORs can measure the accuracy of nanopore sequencing, evaluate the quantitative performance of metagenomic and RNA sequencing, and improve normalisation between samples. CAPTORs can also be customised for clinical diagnoses, correcting systematic sequencing errors and improving the diagnosis of pathogenic BRCA1/2 variants in breast cancer. CAPTORs are a simple and effective method to increase the accuracy and reliability of NGS, enabling comparisons between samples, reagents and laboratories, and supporting the use of nanopore sequencing for clinical diagnosis.

A method for simultaneous detection of small and long RNA biotypes by ribodepleted RNA-Seq

RNA sequencing offers unprecedented access to the transcriptome. Key to this is the identification and quantification of many different species of RNA from the same sample at the same time. In this study we describe a novel protocol for simultaneous detection of coding and non-coding transcripts using modifications to the Ion Total RNA-Seq kit v2 protocol, with integration of QIASeq FastSelect rRNA removal kit. We report highly consistent sequencing libraries can be produced from both frozen high integrity mouse hippocampal tissue and the more challenging post-mortem human tissue. Removal of rRNA using FastSelect was extremely efficient, resulting in less than 1.5% rRNA content in the final library. We identified > 30,000 unique transcripts from all samples, including protein-coding genes and many species of non-coding RNA, in biologically-relevant proportions. Furthermore, the normalized sequencing read count for select genes significantly negatively correlated with Ct values from qRT-PCR analysis from the same samples. These results indicate that this protocol accurately and consistently identifies and quantifies a wide variety of transcripts simultaneously. The highly efficient rRNA depletion, coupled with minimized sample handling and without complicated and high-loss size selection protocols, makes this protocol useful to researchers wishing to investigate whole transcriptomes.

An Easy, Cost-Effective, and Scalable Method to Deplete Human Ribosomal RNA for RNA-seq

RNA sequencing (RNA-seq) is a powerful and increasingly prevalent method to characterize and quantify the transcriptome. Ribosomes are extremely abundant, however, and approximately 80% of total RNA is ribosomal RNA (rRNA). Therefore, to detect and quantify less abundant yet biologically important transcripts such as messenger RNA (mRNA) and long noncoding RNAs (lncRNA), it is essential to minimize the rRNA being sequenced. Although commercial methods exist to deplete rRNA from total RNA samples before sequencing, they are expensive and require specific amounts of input RNA, and the most commonly used kit is no longer available as a stand-alone product. Here, we present an optimized rRNA depletion protocol using RNase H and DNA oligonucleotides complementary to human rRNA transcripts. This protocol includes guidelines for DNA oligo preparation, RNA:DNA hybridization, RNase H cleavage and RNA cleanup, and benchmarking of rRNA depletion. The method is flexible because the user can include additional complementary DNA oligos directed against any abundant transcript in their particular system. Furthermore, the performance of this rRNA depletion approach is comparable to or better than that of commercial kits, at a fraction of the cost and across a wide range of input RNA amounts. © 2021 Wiley Periodicals LLC. Basic Protocol: Specific depletion of rRNA transcripts from human total RNA Support Protocol: Preparation of the rRNA depletion DNA oligonucleotide pool.

A Scalable Strand-Specific Protocol Enabling Full-Length Total RNA Sequencing From Single Cells

RNA sequencing (RNAseq) has been widely used to generate bulk gene expression measurements collected from pools of cells. Only relatively recently have single-cell RNAseq (scRNAseq) methods provided opportunities for gene expression analyses at the single-cell level, allowing researchers to study heterogeneous mixtures of cells at unprecedented resolution. Tumors tend to be composed of heterogeneous cellular mixtures and are frequently the subjects of such analyses. Extensive method developments have led to several protocols for scRNAseq but, owing to the small amounts of RNA in single cells, technical constraints have required compromises. For example, the majority of scRNAseq methods are limited to sequencing only the 3' or 5' termini of transcripts. Other protocols that facilitate full-length transcript profiling tend to capture only polyadenylated mRNAs and are generally limited to processing only 96 cells at a time. Here, we address these limitations and present a novel protocol that allows for the high-throughput sequencing of full-length, total RNA at single-cell resolution. We demonstrate that our method produced strand-specific sequencing data for both polyadenylated and non-polyadenylated transcripts, enabled the profiling of transcript regions beyond only transcript termini, and yielded data rich enough to allow identification of cell types from heterogeneous biological samples.

A smart polymer for sequence-selective binding, pulldown, and release of DNA targets

Selective isolation of DNA is crucial for applications in biology, bionanotechnology, clinical diagnostics and forensics. We herein report a smart methanol-responsive polymer (MeRPy) that can be programmed to bind and separate single- as well as double-stranded DNA targets. Captured targets are quickly isolated and released back into solution by denaturation (sequence-agnostic) or toehold-mediated strand displacement (sequence-selective). The latter mode allows 99.8% efficient removal of unwanted sequences and 79% recovery of highly pure target sequences. We applied MeRPy for the depletion of insulin, glucagon, and transthyretin cDNA from clinical next-generation sequencing (NGS) libraries. This step improved the data quality for low-abundance transcripts in expression profiles of pancreatic tissues. Its low cost, scalability, high stability and ease of use make MeRPy suitable for diverse applications in research and clinical laboratories, including enhancement of NGS libraries, extraction of DNA from biological samples, preparative-scale DNA isolations, and sorting of DNA-labeled non-nucleic acid targets.