A low-bias and sensitive small RNA library preparation method using randomized splint ligation

Insights into the target-directed miRNA degradation mechanism in Drosophila ovarian cell culture

Target-directed miRNA degradation (TDMD) is a process of post-transcriptional regulation of miRNA stability in animals induced by an extended pairing of Ago-bound miRNAs with specialized complementary RNA targets. As suggested by studies on human cell culture, Ago engaged with the extended duplex is recognized by the ZSWIM8 receptor of the Cullin-RING-ligase complex (CRL3), which also contains Cul3, EloB, and EloC proteins. The CRL activity is accelerated by the neddylation of Cul3 with the involvement of the E2 conjugating protein UbcE2M. The CRL ubiquitinates Ago, resulting in proteolysis of Ago and degradation of the released miRNAs. To date, the molecular mechanism of TDMD has not been studied in other species. To further characterize TDMD in animals, we investigated the protein Dora, the Drosophila ortholog of ZSWIM8, in the culture of Drosophila ovarian somatic cells (OSC). We showed that Dora in OSCs localizes in protein granules unrelated to P- and GW-bodies. The dora knockout resulted in the accumulation of multiple miRNAs, including miR-7-5p, and transcriptome-wide affected the mRNA targets of differentially expressed miRNAs. We also showed that Dora associates with proteins of the CRL3 complex, and the depletion of CRL3 components or inhibition of Cul3 neddylation upregulates miR-7-5p. We concluded that the molecular mechanism of TDMD is conserved in humans and Drosophila. Finally, we found that cells without Dora have an impaired Notch signaling pathway, indicating that TDMD in OSCs may contribute to the modulation of the Notch pathway.

Template switching enables chemical probing of native RNA structures

RNAs are often studied in nonnative sequence contexts to facilitate structural studies. However, seemingly innocuous changes to an RNA sequence may perturb the native structure and generate inaccurate or ambiguous structural models. To facilitate the investigation of native RNA secondary structure by selective 2' hydroxyl acylation analyzed by primer extension (SHAPE), we engineered an approach that couples minimal enzymatic steps to RNA chemical probing and mutational profiling (MaP) reverse transcription (RT) methods-a process we call template switching and mutational profiling (Switch-MaP). In Switch-MaP, RT templates and additional library sequences are added postprobing through ligation and template switching, capturing reactivities for every nucleotide. For a candidate SAM-I riboswitch, we compared RNA structure models generated by the Switch-MaP approach to those of traditional primer-based MaP, including RNAs with or without appended structure cassettes. Primer-based MaP masked reactivity data in the 5' and 3' ends of the RNA, producing ambiguous ensembles inconsistent with the conserved SAM-I riboswitch secondary structure. Structure cassettes enabled unambiguous modeling of an aptamer-only construct but introduced nonnative interactions in the full-length riboswitch. In contrast, Switch-MaP provided reactivity data for all nucleotides in each RNA and enabled unambiguous modeling of secondary structure, consistent with the conserved SAM-I fold. Switch-MaP is a straightforward alternative approach to primer-based and cassette-based chemical probing methods that precludes primer masking and the formation of alternative secondary structures due to nonnative sequence elements.

Untacking small RNA profiling and RNA fragment footprinting: Approaches and challenges in library construction

Small RNAs (sRNAs) with sizes ranging from 15 to 50 nucleotides (nt) are critical regulators of gene expression control. Prior studies have shown that sRNAs are involved in a broad range of biological processes, such as organ development, tumorigenesis, and epigenomic regulation; however, emerging evidence unveils a hidden layer of diversity and complexity of endogenously encoded sRNAs profile in eukaryotic organisms, including novel types of sRNAs and the previously unknown post-transcriptional RNA modifications. This underscores the importance for accurate, unbiased detection of sRNAs in various cellular contexts. A multitude of high-throughput methods based on next-generation sequencing (NGS) are developed to decipher the sRNA expression and their modifications. Nonetheless, distinct from mRNA sequencing, the data from sRNA sequencing suffer frequent inconsistencies and high variations emanating from the adapter contaminations and RNA modifications, which overall skew the sRNA libraries. Here, we summarize the sRNA-sequencing approaches, and discuss the considerations and challenges for the strategies and methods of sRNA library construction. The pros and cons of sRNA sequencing have significant implications for implementing RNA fragment footprinting approaches, including CLIP-seq and Ribo-seq. We envision that this review can inspire novel improvements in small RNA sequencing and RNA fragment footprinting in future. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA Processing > Processing of Small RNAs Regulatory RNAs/RNAi/Riboswitches > Biogenesis of Effector Small RNAs.

Characterisation and engineering of a thermophilic RNA ligase from Palaeococcus pacificus

RNA ligases are important enzymes in molecular biology and are highly useful for the manipulation and analysis of nucleic acids, including adapter ligation in next-generation sequencing of microRNAs. Thermophilic RNA ligases belonging to the RNA ligase 3 family are gaining attention for their use in molecular biology, for example a thermophilic RNA ligase from Methanobacterium thermoautotrophicum is commercially available for the adenylation of nucleic acids. Here we extensively characterise a newly identified RNA ligase from the thermophilic archaeon Palaeococcus pacificus (PpaRnl). PpaRnl exhibited significant substrate adenylation activity but low ligation activity across a range of oligonucleotide substrates. Mutation of Lys92 in motif I to alanine, resulted in an enzyme that lacked adenylation activity, but demonstrated improved ligation activity with pre-adenylated substrates (ATP-independent ligation). Subsequent structural characterisation revealed that in this mutant enzyme Lys238 was found in two alternate positions for coordination of the phosphate tail of ATP. In contrast mutation of Lys238 in motif V to glycine via structure-guided engineering enhanced ATP-dependent ligation activity via an arginine residue compensating for the absence of Lys238. Ligation activity for both mutations was higher than the wild-type, with activity observed across a range of oligonucleotide substrates with varying sequence and secondary structure.

An improved method for the highly specific detection of transcription start sites

Precise detection of the transcriptional start site (TSS) is a key for characterizing transcriptional regulation of genes and for annotation of newly sequenced genomes. Here, we describe the development of an improved method, designated 'TSS-seq2.' This method is an iterative improvement of TSS-seq, a previously published enzymatic cap-structure conversion method to detect TSSs in base sequences. By modifying the original procedure, including by introducing split ligation at the key cap-selection step, the yield and the accuracy of the reaction has been substantially improved. For example, TSS-seq2 can be conducted using as little as 5 ng of total RNA with an overall accuracy of 96%; this yield a less-biased and more precise detection of TSS. We then applied TSS-seq2 for TSS analysis of four plant species that had not yet been analyzed by any previous TSS method.

Terminal modifications independent cell-free RNA sequencing enables sensitive early cancer detection and classification

Cell-free RNAs (cfRNAs) offer an opportunity to detect diseases from a transcriptomic perspective, however, existing techniques have fallen short in generating a comprehensive cell-free transcriptome profile. We develop a sensitive library preparation method that is robust down to 100 µl input plasma to analyze cfRNAs independent of their 5'-end modifications. We show that it outperforms adapter ligation-based method in detecting a greater number of cfRNA species. We perform transcriptome-wide characterizations in 165 lung cancer, 30 breast cancer, 37 colorectal cancer, 55 gastric cancer, 15 liver cancer, and 133 cancer-free participants and demonstrate its ability to identify transcriptomic changes occurring in early-stage tumors. We also leverage machine learning analyses on the differentially expressed cfRNA signatures and reveal their robust performance in cancer detection and classification. Our work sets the stage for in-depth study of the cfRNA repertoire and highlights the value of cfRNAs as cancer biomarkers in clinical applications.

FindAdapt: A python package for fast and accurate adapter detection in small RNA sequencing

Adapter trimming is an essential step for analyzing small RNA sequencing data, where reads are generally longer than target RNAs ranging from 18 to 30 bp. Most adapter trimming tools require adapter information as input. However, adapter information is hard to access, specified incorrectly, or not provided with publicly available datasets, hampering their reproducibility and reusability. Manual identification of adapter patterns from raw reads is labor-intensive and error-prone. Moreover, the use of randomized adapters to reduce ligation biases during library preparation makes adapter detection even more challenging. Here, we present FindAdapt, a Python package for fast and accurate detection of adapter patterns without relying on prior information. We demonstrated that FindAdapt was far superior to existing approaches. It identified adapters successfully in 180 simulation datasets with diverse read structures and 3,184 real datasets covering a variety of commercial and customized small RNA library preparation kits. FindAdapt is stand-alone software that can be easily integrated into small RNA sequencing analysis pipelines.

CLT-seq as a universal homopolymer-sequencing concept reveals poly(A)-tail-tuned ncRNA regulation

Dynamic tuning of the poly(A) tail is a crucial mechanism for controlling translation and stability of eukaryotic mRNA. Achieving a comprehensive understanding of how this regulation occurs requires unbiased abundance quantification of poly(A)-tail transcripts and simple poly(A)-length measurement using high-throughput sequencing platforms. Current methods have limitations due to complicated setups and elaborate library preparation plans. To address this, we introduce central limit theorem (CLT)-managed RNA-seq (CLT-seq), a simple and straightforward homopolymer-sequencing method. In CLT-seq, an anchor-free oligo(dT) primer rapidly binds to and unbinds from anywhere along the poly(A) tail string, leading to position-directed reverse transcription with equal probability. The CLT mechanism enables the synthesized poly(T) lengths, which correspond to the templated segment of the poly(A) tail, to distribute normally. Based on a well-fitted pseudogaussian-derived poly(A)-poly(T) conversion model, the actual poly(A)-tail profile is reconstructed from the acquired poly(T)-length profile through matrix operations. CLT-seq follows a simple procedure without requiring RNA-related pre-treatment, enrichment or selection, and the CLT-shortened poly(T) stretches are more compatible with existing sequencing platforms. This proof-of-concept approach facilitates direct homopolymer base-calling and features unbiased RNA-seq. Therefore, CLT-seq provides unbiased, robust and cost-efficient transcriptome-wide poly(A)-tail profiling. We demonstrate that CLT-seq on the most common Illumina platform delivers reliable poly(A)-tail profiling at a transcriptome-wide scale in human cellular contexts. We find that the poly(A)-tail-tuned ncRNA regulation undergoes a dynamic, complex process similar to mRNA regulation. Overall, CLT-seq offers a simplified, effective and economical approach to investigate poly(A)-tail regulation, with potential implications for understanding gene expression and identifying therapeutic targets.

Trans-species microRNA loci in the parasitic plant Cuscuta campestris have a U6-like snRNA promoter

Small regulatory RNAs can move between organisms and regulate gene expression in the recipient. Whether the trans-species small RNAs being exported are distinguished from the normal endogenous small RNAs of the source organism is not known. The parasitic plant Cuscuta campestris (dodder) produces many microRNAs that specifically accumulate at the host-parasite interface, several of which have trans-species activity. We found that induction of C. campestris interface-induced microRNAs is similar regardless of host species and occurs in C. campestris haustoria produced in the absence of any host. The loci-encoding C. campestris interface-induced microRNAs are distinguished by a common cis-regulatory element. This element is identical to a conserved upstream sequence element (USE) used by plant small nuclear RNA loci. The properties of the interface-induced microRNA primary transcripts strongly suggest that they are produced via U6-like transcription by RNA polymerase III. The USE promotes accumulation of interface-induced miRNAs (IIMs) in a heterologous system. This promoter element distinguishes C. campestris IIM loci from other plant small RNAs. Our data suggest that C. campestris IIMs are produced in a manner distinct from canonical miRNAs. All confirmed C. campestris microRNAs with documented trans-species activity are interface-induced and possess these features. We speculate that RNA polymerase III transcription of IIMs may allow these miRNAs to be exported to hosts.

Reassessment of miRNA variant (isomiRs) composition by small RNA sequencing

IsomiRs, sequence variants of mature microRNAs, are usually detected and quantified using high-throughput sequencing. Many examples of their biological relevance have been reported, but sequencing artifacts identified as artificial variants might bias biological inference and therefore need to be ideally avoided. We conducted a comprehensive evaluation of 10 different small RNA sequencing protocols, exploring both a theoretically isomiR-free pool of synthetic miRNAs and HEK293T cells. We calculated that, with the exception of two protocols, less than 5% of miRNA reads can be attributed to library preparation artifacts. Randomized-end adapter protocols showed superior accuracy, with 40% of true biological isomiRs. Nevertheless, we demonstrate concordance across protocols for selected miRNAs in non-templated uridyl additions. Notably, NTA-U calling and isomiR target prediction can be inaccurate when using protocols with poor single-nucleotide resolution. Our results highlight the relevance of protocol choice for biological isomiRs detection and annotation, which has key potential implications for biomedical applications.

A combined miRNA-piRNA signature in the serum and urine of rabbits infected with Toxoplasma gondii oocysts

BACKGROUND

Increasing evidence has shown that non-coding RNA (ncRNA) molecules play fundamental roles in cells, and many are stable in body fluids as circulating RNAs. Study on these ncRNAs will provide insights into toxoplasmosis pathophysiology and/or help reveal diagnostic biomarkers.

METHODS

We performed a high-throughput RNA-Seq study to comprehensively profile the microRNAs (miRNAs) and PIWI-interacting RNAs (piRNAs) in rabbit serum and urine after infection with Toxoplasma gondii oocysts during the whole infection process.

RESULTS

Total RNA extracted from serum and urine samples of acutely infected [8 days post-infection (DPI)], chronically infected (70 DPI) and uninfected rabbits were subjected to genome-wide small RNA sequencing. We identified 2089 miRNAs and 2224 novel piRNAs from the rabbit sera associated with T. gondii infection. Meanwhile, a total of 518 miRNAs and 4182 novel piRNAs were identified in the rabbit urine associated with T. gondii infection. Of these identified small ncRNAs, 1178 and 1317 serum miRNAs and 311 and 294 urine miRNAs were identified as differentially expressed (DE) miRNAs in the acute and chronic stages of infections, respectively. A total of 1748 and 1814 serum piRNAs and 597 and 708 urine piRNAs were found in the acute and chronic infection stages, respectively. Of these dysregulated ncRNAs, a total of 88 common DE miRNAs and 120 DE novel piRNAs were found in both serum and urine samples of infected rabbits.

CONCLUSIONS

These findings provide valuable data for revealing the physiology of herbivore toxoplasmosis caused by oocyst infection. Circulating ncRNAs identified in this study are potential novel diagnostic biomarkers for the detection/diagnosis of toxoplasmosis in herbivorous animals.

Enhanced transcriptome-wide RNA G-quadruplex sequencing for low RNA input samples with rG4-seq 2.0

Background

RNA G-quadruplexes (rG4s) are non-canonical structural motifs that have diverse functional and regulatory roles, for instance in transcription termination, alternative splicing, mRNA localization and stabilization, and translational process. We recently developed the RNA G-quadruplex structure sequencing (rG4-seq) technique and described rG4s in both eukaryotic and prokaryotic transcriptomes. However, rG4-seq suffers from a complicated gel purification step and limited PCR product yield, thus requiring a high amount of RNA input, which limits its applicability in more physiologically or clinically relevant studies often characterized by the limited availability of biological material and low RNA abundance. Here, we redesign and enhance the workflow of rG4-seq to address this issue.

Results

We developed rG4-seq 2.0 by introducing a new ssDNA adapter containing deoxyuridine during library preparation to enhance library quality with no gel purification step, less PCR amplification cycles and higher yield of PCR products. We demonstrate that rG4-seq 2.0 produces high-quality cDNA libraries that support reliable and reproducible rG4 identification at varying RNA inputs, including RNA mounts as low as 10 ng. rG4-seq 2.0 also improved the rG4-seq calling outcome and nucleotide bias in rG4 detection persistent in rG4-seq 1.0. We further provide in vitro mapping of rG4 in the HEK293T cell line, and recommendations for assessing RNA input and sequencing depth for individual rG4 studies based on transcript abundance.

Conclusions

rG4-seq 2.0 can improve the identification and study of rG4s in low abundance transcripts, and our findings can provide insights to optimize cDNA library preparation in other related methods.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01448-3.

Signature of circulating small non-coding RNAs during early fracture healing in mice

Fracture healing is a complex process with multiple overlapping metabolic and differentiation phases. Small non-coding RNAs are involved in the regulation of fracture healing and their presence in circulation is under current interest due to their obvious value as potential biomarkers. Circulating microRNAs (miRNAs) have been characterized to some extent but the current knowledge on tRNA-derived small RNA fragments (tsRNAs) is relatively scarce, especially in circulation.

In this study, the spectrum of circulating miRNAs and tsRNAs was analysed by next generation sequencing to show their differential expression during fracture healing in vivo. Analysed tsRNA fragments included stress-induced translation interfering tRNA fragments (tiRNAs or tRNA halves) and internal tRNA fragments (i-tRF), within the size range of 28–36 bp. To unveil the expression of these non-coding RNAs, genome-wide analysis was performed on two months old C57BL/6 mice on days 1, 5, 7, 10, and 14 (D1, D5, D7, D10, and D14) after a closed tibial fracture.

Valine isoacceptor tRNA-derived Val-AAC 5′end and Val-CAC 5′end fragments were the major types of 5′end tiRNAs in circulation, comprising about 65 % of the total counts. Their expression was not affected by fracture. After a fracture, the levels of two 5′end tiRNAs Lys-TTT 5′ and Lys-CTT 5′ were decreased and His-GTG 5′ was increased through D1-D14. The level of miR-451a was decreased on the first post-fracture day (D1), whereas miR-328-3p, miR-133a-3p, miR-375-3p, miR-423-5p, and miR-150-5p were increased post-fracture. These data provide evidence on how fracture healing could provoke systemic metabolic effects and further pinpoint the potential of small non-coding RNAs as biomarkers for tissue regeneration.

A comparative analysis of single cell small RNA sequencing data reveals heterogeneous isomiR expression and regulation

MicroRNAs (miRNAs) are non-coding small RNAs which play a critical role in the regulation of gene expression in cells. It is known that miRNAs are often expressed as multiple isoforms, called isomiRs, which may have alternative regulatory functions. Despite the recent development of several single cell small RNA sequencing protocols, these methods have not been leveraged to investigate isomiR expression and regulation to better understand their role on a single cell level. Here we integrate sequencing data from three independent studies and find substantial differences in isomiR composition that suggest that cell autonomous mechanisms may drive isomiR processing. We also find evidence of altered regulatory functions of different classes of isomiRs, when compared to their respective wild-type miRNA, which supports a biological role for many of the isomiRs that are expressed.

Solid-Phase-Supported Chemoenzymatic Synthesis of a Light-Activatable tRNA Derivative

Herein, we present a multi-cycle chemoenzymatic synthesis of modified RNA with simplified solid-phase handling to overcome size limitations of RNA synthesis. It combines the advantages of classical chemical solid-phase synthesis and enzymatic synthesis using magnetic streptavidin beads and biotinylated RNA. Successful introduction of light-controllable RNA nucleotides into the tRNAMet sequence was confirmed by gel electrophoresis and mass spectrometry. The methods tolerate modifications in the RNA phosphodiester backbone and allow introductions of photocaged and photoswitchable nucleotides as well as photocleavable strand breaks and fluorophores.

Rolling circle reverse transcription enables high fidelity nanopore sequencing of small RNA

Small RNAs (sRNAs) are an important group of non-coding RNAs that have great potential as diagnostic and prognostic biomarkers for treatment of a wide variety of diseases. The portability and affordability of nanopore sequencing technology makes it ideal for point of care and low resource settings. Currently sRNAs can't be reliably sequenced on the nanopore platform due to the short size of sRNAs and high error rate of the nanopore sequencer. Here, we developed a highly efficient nanopore-based sequencing strategy for sRNAs (SR-Cat-Seq) in which sRNAs are ligated to an adapter, circularized, and undergo rolling circle reverse transcription to generate concatemeric cDNA. After sequencing, the resulting tandem repeat sequences within the individual cDNA can be aligned to generate highly accurate consensus sequences. We compared our sequencing strategy with other sRNA sequencing methods on a short-read sequencing platform and demonstrated that SR-Cat-Seq can obtain low bias and highly accurate sRNA transcriptomes. Therefore, our method could enable nanopore sequencing for sRNA-based diagnostics and other applications.

5´XP sRNA-seq: efficient identification of transcripts with and without 5´ phosphorylation reveals evolutionary conserved small RNA

Small RNA (sRNA) sequencing has been critical for our understanding of many cellular processes, including gene regulation. Nonetheless, the varying biochemical properties of sRNA, such as 5´ nucleotide modifications, make many sRNA subspecies incompatible with common protocols for sRNA sequencing. Here we describe 5XP-seq that outlines a novel strategy that captures a more complete picture of sRNA. By tagging 5´P sRNA during library preparation, 5XP-seq combines an open approach that includes all types of 5'-terminal modifications (5´X), with a selective approach for 5-phosphorylated sRNA (5´P). We show that 5XP-seq not only enriches phosphorylated miRNA and piRNA but successfully discriminates these sRNA from all other sRNA species. We further demonstrate the importance of this strategy by successful inter-species validation of sRNAs that would have otherwise failed, including human to insect translation of several tRNA (tRFs) and rRNA (rRFs) fragments. By combining 5´ insensitive library strategies with 5´ sensitive tagging, we have successfully tackled an intrinsic bias in modern sRNA sequencing that will help us reveal the true complexity and the evolutionary significance of the sRNA world.

Chemical capping improves template switching and enhances sequencing of small RNAs

Template-switching reverse transcription is widely used in RNA sequencing for low-input and low-quality samples, including RNA from single cells or formalin-fixed paraffin-embedded (FFPE) tissues. Previously, we identified the native eukaryotic mRNA 5′ cap as a key structural element for enhancing template switching efficiency. Here, we introduce CapTS-seq, a new strategy for sequencing small RNAs that combines chemical capping and template switching. We probed a variety of non-native synthetic cap structures and found that an unmethylated guanosine triphosphate cap led to the lowest bias and highest efficiency for template switching. Through cross-examination of different nucleotides at the cap position, our data provided unequivocal evidence that the 5′ cap acts as a template for the first nucleotide in reverse transcriptase-mediated post-templated addition to the emerging cDNA—a key feature to propel template switching. We deployed CapTS-seq for sequencing synthetic miRNAs, human total brain and liver FFPE RNA, and demonstrated that it consistently improves library quality for miRNAs in comparison with a gold standard template switching-based small RNA-seq kit.

ggVennDiagram: An Intuitive, Easy-to-Use, and Highly Customizable R Package to Generate Venn Diagram

Venn diagrams are widely used diagrams to show the set relationships in biomedical studies. In this study, we developed ggVennDiagram, an R package that could automatically generate high-quality Venn diagrams with two to seven sets. The ggVennDiagram is built based on ggplot2, and it integrates the advantages of existing packages, such as venn, RVenn, VennDiagram, and sf. Satisfactory results can be obtained with minimal configurations. Furthermore, we designed comprehensive objects to store the entire data of the Venn diagram, which allowed free access to both intersection values and Venn plot sub-elements, such as set label/edge and region label/filling. Therefore, high customization of every Venn plot sub-element can be fulfilled without increasing the cost of learning when the user is familiar with ggplot2 methods. To date, ggVennDiagram has been cited in more than 10 publications, and its source code repository has been starred by more than 140 GitHub users, suggesting a great potential in applications. The package is an open-source software released under the GPL-3 license, and it is freely available through CRAN (https://cran.r-project.org/package=ggVennDiagram).

A New Specific and Sensitive RT-qPCR Method Based on Splinted 5' Ligation for the Quantitative Detection of RNA Species Shorter than microRNAs

Recently, we discovered a new family of unusually short RNAs mapping to 5.8S ribosomal RNA (rRNA) and which we named dodecaRNAs (doRNAs), according to the number of core nucleotides (12 nt) their members contain. To confirm these small RNA-sequencing (RNA-Seq) data, validate the existence of the two overly abundant doRNAs-the minimal core 12-nt doRNA sequence and its + 1-nt variant bearing a 5' Cytosine, C-doRNA-and streamline their analysis, we developed a new specific and sensitive splinted 5' ligation reverse transcription (RT)-quantitative polymerase chain reaction (qPCR) method. This method is based on a splint-assisted ligation of an adapter to the 5' end of doRNAs, followed by RT-qPCR amplification and quantitation. Our optimized protocol, which may discriminate between doRNA, C-doRNA, mutated and precursor sequences, can accurately detect as low as 240 copies and is quantitatively linear over a range of 7 logs. This method provides a unique tool to expand and facilitate studies exploring the molecular and cellular biology of RNA species shorter than microRNAs.

Single-cell microRNA sequencing method comparison and application to cell lines and circulating lung tumor cells

Molecular single cell analyses provide insights into physiological and pathological processes. Here, in a stepwise approach, we first evaluate 19 protocols for single cell small RNA sequencing on MCF7 cells spiked with 1 pg of 1,006 miRNAs. Second, we analyze MCF7 single cell equivalents of the eight best protocols. Third, we sequence single cells from eight different cell lines and 67 circulating tumor cells (CTCs) from seven SCLC patients. Altogether, we analyze 244 different samples. We observe high reproducibility within protocols and reads covered a broad spectrum of RNAs. For the 67 CTCs, we detect a median of 68 miRNAs, with 10 miRNAs being expressed in 90% of tested cells. Enrichment analysis suggested the lung as the most likely organ of origin and enrichment of cancer-related categories. Even the identification of non-annotated candidate miRNAs was feasible, underlining the potential of single cell small RNA sequencing.

Technologies for small non-coding RNA sequencing at the single-cell level are less mature than for sequencing mRNAs. Here the authors evaluate available protocols for analysis of circulating lung cancer tumour cells.

Small RNA-Sequencing: Approaches and Considerations for miRNA Analysis

MicroRNAs (miRNAs) are a class of small RNA molecules that have an important regulatory role in multiple physiological and pathological processes. Their disease-specific profiles and presence in biofluids are properties that enable miRNAs to be employed as non-invasive biomarkers. In the past decades, several methods have been developed for miRNA analysis, including small RNA sequencing (RNA-seq). Small RNA-seq enables genome-wide profiling and analysis of known, as well as novel, miRNA variants. Moreover, its high sensitivity allows for profiling of low input samples such as liquid biopsies, which have now found applications in diagnostics and prognostics. Still, due to technical bias and the limited ability to capture the true miRNA representation, its potential remains unfulfilled. The introduction of many new small RNA-seq approaches that tried to minimize this bias, has led to the existence of the many small RNA-seq protocols seen today. Here, we review all current approaches to cDNA library construction used during the small RNA-seq workflow, with particular focus on their implementation in commercially available protocols. We provide an overview of each protocol and discuss their applicability. We also review recent benchmarking studies comparing each protocol's performance and summarize the major conclusions that can be gathered from their usage. The result documents variable performance of the protocols and highlights their different applications in miRNA research. Taken together, our review provides a comprehensive overview of all the current small RNA-seq approaches, summarizes their strengths and weaknesses, and provides guidelines for their applications in miRNA research.

Transfer RNA-Derived Small RNAs: Another Layer of Gene Regulation and Novel Targets for Disease Therapeutics

Decades after identification as essential for protein synthesis, transfer RNAs (tRNAs) have been implicated in various cellular processes beyond translation. tRNA-derived small RNAs (tsRNAs), referred to as tRNA-derived fragments (tRFs) or tRNA-derived, stress-induced RNAs (tiRNAs), are produced by cleavage at different sites from mature or pre-tRNAs. They are classified into six major types representing potentially thousands of unique sequences and have been implicated to play a wide variety of regulatory roles in maintaining normal homeostasis, cancer cell viability, tumorigenesis, ribosome biogenesis, chromatin remodeling, translational regulation, intergenerational inheritance, retrotransposon regulation, and viral replication. However, the detailed mechanisms governing these processes remain unknown. Aberrant expression of tsRNAs is found in various human disease conditions, suggesting that a further understanding of the regulatory role of tsRNAs will assist in identifying novel biomarkers, potential therapeutic targets, and gene-regulatory tools. Here, we highlight the classification, biogenesis, and biological role of tsRNAs in regulatory mechanisms of normal and disease states.