Identification and remediation of biases in the activity of RNA ligases in small-RNA deep sequencing

tRNA expression and modification landscapes, and their dynamics during zebrafish embryo development

tRNA genes exist in multiple copies in the genome of all organisms across the three domains of life. Besides the sequence differences across tRNA copies, extensive post-transcriptional modification adds a further layer to tRNA diversification. Whilst the crucial role of tRNAs as adapter molecules in protein translation is well established, whether all tRNAs are actually expressed, and whether the differences across isodecoders play any regulatory role is only recently being uncovered. Here we built upon recent developments in the use of NGS-based methods for RNA modification detection and developed tRAM-seq, an experimental protocol and in silico analysis pipeline to investigate tRNA expression and modification. Using tRAM-seq, we analysed the full ensemble of nucleo-cytoplasmic and mitochondrial tRNAs during embryonic development of the model vertebrate zebrafish. We show that the repertoire of tRNAs changes during development, with an apparent major switch in tRNA isodecoder expression and modification profile taking place around the start of gastrulation. Taken together, our findings suggest the existence of a general reprogramming of the expressed tRNA pool, possibly gearing the translational machinery for distinct stages of the delicate and crucial process of embryo development.

PT-seq: A method for metagenomic analysis of phosphorothioate epigenetics in complex microbial communities

Among dozens of known epigenetic marks, naturally occurring phosphorothioate (PT) DNA modifications are unique in replacing a non-bridging phosphate oxygen with redox-active sulfur and function in prokaryotic restriction-modification and transcriptional regulation. Interest in PTs has grown due to the widespread distribution of the dnd, ssp, and brx genes among bacteria and archaea, as well as the discovery of PTs in 5-10% of gut microbes. Efforts to map PTs in complex microbiomes using existing next-generation and direct sequencing technologies have failed due to poor sensitivity. Here we developed PT-seq as a high-sensitivity method to quantitatively map PTs across genomes and metagenomically identify PT-containing microbes in complex genomic mixtures. Like other methods for mapping PTs in individual genomes, PT-seq exploits targeted DNA strand cleavage at PTs by iodine, followed by sequencing library construction using ligation or template switching approaches. However, PT-specific sequencing reads are dramatically increased by adding steps to heat denature the DNA, block pre-existing 3'-ends, fragment DNA after T-tailing, and enrich iodine-induced breaks using biotin-labeling and streptavidin beads capture. Iterative optimization of the sensitivity and specificity of PT-seq is demonstrated with individual bacteria and human fecal DNA.

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.

A robust method for measuring aminoacylation through tRNA-Seq

Current methods to quantify the fraction of aminoacylated tRNAs, also known as the tRNA charge, are limited by issues with either low throughput, precision, and/or accuracy. Here, we present an optimized charge tRNA-Seq method that combines previous developments with newly described approaches to establish a protocol for precise and accurate tRNA charge measurements. We verify that this protocol provides robust quantification of tRNA aminoacylation and we provide an end-to-end method that scales to hundreds of samples including software for data processing. Additionally, we show that this method supports measurements of relative tRNA expression levels and can be used to infer tRNA modifications through reverse transcription misincorporations, thereby supporting multipurpose applications in tRNA biology.

Target-directed microRNA degradation: Mechanisms, significance, and functional implications

MicroRNAs (miRNAs) are small non-coding RNAs that play a fundamental role in enabling miRNA-mediated target repression, a post-transcriptional gene regulatory mechanism preserved across metazoans. Loss of certain animal miRNA genes can lead to developmental abnormalities, disease, and various degrees of embryonic lethality. These short RNAs normally guide Argonaute (AGO) proteins to target RNAs, which are in turn translationally repressed and destabilized, silencing the target to fine-tune gene expression and maintain cellular homeostasis. Delineating miRNA-mediated target decay has been thoroughly examined in thousands of studies, yet despite these exhaustive studies, comparatively less is known about how and why miRNAs are directed for decay. Several key observations over the years have noted instances of rapid miRNA turnover, suggesting endogenous means for animals to induce miRNA degradation. Recently, it was revealed that certain targets, so-called target-directed miRNA degradation (TDMD) triggers, can "trigger" miRNA decay through inducing proteolysis of AGO and thereby the bound miRNA. This process is mediated in animals via the ZSWIM8 ubiquitin ligase complex, which is recruited to AGO during engagement with triggers. Since its discovery, several studies have identified that ZSWIM8 and TDMD are indispensable for proper animal development. Given the rapid expansion of this field of study, here, we summarize the key findings that have led to and followed the discovery of ZSWIM8-dependent TDMD. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Turnover and Surveillance > Turnover/Surveillance Mechanisms RNA in Disease and Development > RNA in Development.

Evaluation of microRNA variant maturation prior to genome edition

Assessment of the functionality of individual microRNA/target sites is a crucial issue. Genome editing techniques should theoretically permit a fine functional exploration of such interactions, allowing the mutation of microRNAs or individual binding sites in a complete in vivo setting, therefore abrogating or restoring individual interactions on demand. A major limitation to this experimental strategy is the influence of microRNA sequence on its accumulation level, which introduces a confounding effect when assessing phenotypic rescue by compensatorily mutated microRNA and target site. Here we describe a simple assay to identify microRNA variants most likely to accumulate at wild-type levels even though their sequence has been mutated. In this assay, quantification of a reporter construct in cultured cells predicts the efficiency of an early biogenesis step, the Drosha-dependent cleavage of microRNA precursors, which appears to be a major determinant of microRNA accumulation in our variant collection. This system allowed the generation of a mutant Drosophila strain expressing a bantam microRNA variant at wild-type levels.

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.

Preparation and processing of dried blood spots for microRNA sequencing

Dried blood spots (DBS) are biological samples commonly collected from newborns and in geographic areas distanced from laboratory settings for the purposes of disease testing and identification. MicroRNAs (miRNAs)-small non-coding RNAs that regulate gene activity at the post-transcriptional level-are emerging as critical markers and mediators of disease, including cancer, infectious diseases, and mental disorders. This protocol describes optimized procedural steps for utilizing DBS as a reliable source of biological material for obtaining peripheral miRNA expression profiles. We outline key practices, such as the method of DBS rehydration that maximizes RNA extraction yield, and the use of degenerate oligonucleotide adapters to mitigate ligase-dependent biases that are associated with small RNA sequencing. The standardization of miRNA readout from DBS offers numerous benefits: cost-effectiveness in sample collection and processing, enhanced reliability and consistency of miRNA profiling, and minimal invasiveness that facilitates repeated testing and retention of participants. The use of DBS-based miRNA sequencing is a promising method to investigate disease mechanisms and to advance personalized medicine.

Reverse transcription-free digital-quantitative-PCR for microRNA analysis

MicroRNAs (miRNAs) are non-coding RNA sequences that regulate many biological processes and have become central targets of biomedical research. However, their naturally low abundances in biological samples necessitates the development of sensitive analytical techniques to conduct routine miRNA measurements in research laboratories. Digital PCR has the potential to meet this need because of its single-molecule detection capabilities, but PCR analyses of miRNAs are slowed by the ligation and reverse transcription steps first required to prepare samples. This report describes the development of a method to rapidly quantify miRNA in digital microwell arrays using base-stacking digital-quantitative-PCR (BS-dqPCR). BS-dqPCR expedites miRNA measurements by eliminating the need for ligation and reverse transcription steps, which reduces the time and cost compared to conventional miRNA PCR analyses. Under standard PCR thermocycling conditions, digital signals from miRNA samples were lower than expected, while signals from blanks were high. Therefore, a novel asymmetric thermocycling program was developed that maximized on-target signal from miRNA while minimizing non-specific amplification. The analytical response of BS-dqPCR was then evaluated over a range of miRNA concentrations. The digital PCR dimension increased in signal with increasing miRNA copy numbers. When the digital signal saturated, the quantitative PCR dimension readily discerned miRNA copy number differences. Overall, BS-dqPCR provides rapid, high-sensitivity measurements of miRNA over a wide dynamic range, which demonstrates its utility for routine miRNA analyses.

Spatial mapping of the total transcriptome by in situ polyadenylation

Spatial transcriptomics reveals the spatial context of gene expression, but current methods are limited to assaying polyadenylated (A-tailed) RNA transcripts. Here we demonstrate that enzymatic in situ polyadenylation of RNA enables detection of the full spectrum of RNAs, expanding the scope of sequencing-based spatial transcriptomics to the total transcriptome. We demonstrate that our spatial total RNA-sequencing (STRS) approach captures coding RNAs, noncoding RNAs and viral RNAs. We apply STRS to study skeletal muscle regeneration and viral-induced myocarditis. Our analyses reveal the spatial patterns of noncoding RNA expression with near-cellular resolution, identify spatially defined expression of noncoding transcripts in skeletal muscle regeneration and highlight host transcriptional responses associated with local viral RNA abundance. STRS requires adding only one step to the widely used Visium spatial total RNA-sequencing protocol from 10x Genomics, and thus could be easily adopted to enable new insights into spatial gene regulation and biology.

YBX1-interacting small RNAs and RUNX2 can be blocked in primary bone cancer using CADD522

Primary bone cancer (PBC) comprises several subtypes each underpinned by distinctive genetic drivers. This driver diversity produces novel morphological features and clinical behaviour that serendipitously makes PBC an excellent metastasis model. Here, we report that some transfer RNA-derived small RNAs termed tRNA fragments (tRFs) perform as a constitutive tumour suppressor mechanism by blunting a potential pro-metastatic protein-RNA interaction. This mechanism is reduced in PBC progression with a gradual loss of tRNAGly^TCC cleavage into 5' end tRF-Gly^TCC when comparing low-grade, intermediate-grade and high-grade patient tumours. We detected recurrent activation of miR-140 leading to upregulated RUNX2 expression in high-grade patient tumours. Both tRF-Gly^TCC and RUNX2 share a sequence motif in their 3' ends that matches the YBX1 recognition site known to stabilise pro-metastatic mRNAs. Investigating some aspects of this interaction network, gain- and loss-of-function experiments using small RNA mimics and antisense LNAs, respectively, showed that ectopic tRF-Gly^TCC reduced RUNX2 expression and dispersed 3D micromass architecture in vitro. iCLIP sequencing revealed YBX1 physical binding to the 3' UTR of RUNX2. The interaction between YBX1, tRF-Gly^TCC and RUNX2 led to the development of the RUNX2 inhibitor CADD522 as a PBC treatment. CADD522 assessment in vitro revealed significant effects on PBC cell behaviour. In xenograft mouse models, CADD522 as a single agent without surgery significantly reduced tumour volume, increased overall and metastasis-free survival and reduced cancer-induced bone disease. Our results provide insight into PBC molecular abnormalities that have led to the identification of new targets and a new therapeutic.

choros: correction of sequence-based biases for accurate quantification of ribosome profiling data

Ribosome profiling quantifies translation genome-wide by sequencing ribosome-protected fragments, or footprints. Its single-codon resolution allows identification of translation regulation, such as ribosome stalls or pauses, on individual genes. However, enzyme preferences during library preparation lead to pervasive sequence artifacts that obscure translation dynamics. Widespread over- and under-representation of ribosome footprints can dominate local footprint densities and skew estimates of elongation rates by up to five fold. To address these biases and uncover true patterns of translation, we present choros, a computational method that models ribosome footprint distributions to provide bias-corrected footprint counts. choros uses negative binomial regression to accurately estimate two sets of parameters: (i) biological contributions from codon-specific translation elongation rates; and (ii) technical contributions from nuclease digestion and ligation efficiencies. We use these parameter estimates to generate bias correction factors that eliminate sequence artifacts. Applying choros to multiple ribosome profiling datasets, we are able to accurately quantify and attenuate ligation biases to provide more faithful measurements of ribosome distribution. We show that a pattern interpreted as pervasive ribosome pausing near the beginning of coding regions is likely to arise from technical biases. Incorporating choros into standard analysis pipelines will improve biological discovery from measurements of translation.

Small RNA Profiling by Next-Generation Sequencing Using High-Definition Adapters

Next-generation sequencing (NGS) of small RNA (sRNA) cDNA libraries permits the identification and characterization of sRNA species de novo. However, the method through which these libraries are constructed can often introduce artifacts such as over- or underrepresentation of specific sequences or adapter oligonucleotides due to sequence biases held by the enzymes used. In this chapter we describe a protocol for sRNA library construction making use of high-definition (HD) adapters for the Illumina sequencing platform, which reduce ligation bias. This protocol leads to drastically reduced direct 5'/3' adapter ligation products and can be used for the synthesis of sRNA libraries from total RNA or sRNA of various plant, animal, and fungal samples. This protocol also includes a method for total RNA extraction from plant leaf and cultured cells or body fluids.

Transfer RNAs-derived small RNAs and their application potential in multiple diseases

The role of tRNAs is best known as adapter components of translational machinery. According to the central dogma of molecular biology, DNA is transcribed to RNA and in turn is translated into proteins, in which tRNA outstands by its role of the cellular courier. Recent studies have led to the revision of the canonical function of transfer RNAs (tRNAs), which indicates that tRNAs also serve as a source for short non-coding RNAs called tRNA-derived small RNAs (tsRNAs). tsRNAs play key roles in cellular processes by modulating complicated regulatory networks beyond translation and are widely involved in multiple diseases. Herein, the biogenesis and classification of tsRNAs were firstly clarified. tsRNAs are generated from pre-tRNAs or mature tRNAs and are classified into tRNA-derived fragments (tRFs) and tRNA halves (tiRNA). The tRFs include five types according to the incision loci: tRF-1, tRF-2, tRF-3, tRF-5 and i-tRF which contain 3′ tiRNA and 5′ tiRNA. The functions of tsRNAs and their regulation mechanisms involved in disease processes are systematically summarized as well. The mechanisms can elaborate on the specific regulation of tsRNAs. In conclusion, the current research suggests that tsRNAs are promising targets for modulating pathological processes, such as breast cancer, ischemic stroke, respiratory syncytial virus, osteoporosis and so on, and maintain vital clinical implications in diagnosis and therapeutics of various diseases.

Duck sewage source coliphage P762 can lyse STEC and APEC

Multiple pathogenic types or serotypes restrict treatment for colibacillosis. In addition, rising antibiotic resistance has heightened public awareness to prevent and control pathogenic Escherichia coli. The bacteriophage is a viable technique to treat colibacillosis as an alternative to antibiotics. P762, a coliphage isolated from duck farm sewage, was demonstrated to cloud lyse Shiga toxin-producing Escherichia Coli serotypes O157 and non-O157 (17/39), Avian pathogenic E. coli covered serotype O78, O83, and O9 (5/19), and other pathogenic Escherichia coli (5/17). Additional fundamental biological characteristics analysis revealed that P762 is stable at pH 3 ~ 11 and temperature between 4 °C and 60 °C, and its optimum multiplicity of infection (MOI) is 0.1. The one-step curve of P762 exhibited three bursts of growth stage: two rapid and one slow stage. Furthermore, the first rapid burst size is 80 CFU/PFU, the burst size of the slow stage is 10 CFU/PFU, and the second rapid burst size is about 990 CFU/PFU. In addition, P762 can form a "halo" on a double agar plate, implying that the phage secretes depolymerase. With 95.14% identity and 90% query coverage, genome sequence analysis revealed that P762 is most closely related to Escherichia phage DY1, which belongs to the genus Kayfunavirus. After screening using RAST and VFDB, no virulence factors were discovered in P762. In vitro antibacterial tests revealed that P762 has high bactericidal activity in lettuce leaves contaminated with STEC. In conclusion, phage P762 might be employed in the future to prevent and control pathogenic Escherichia coli.

Next-generation sequencing: A new avenue to understand viral RNA-protein interactions

The genomes of RNA viruses present an astonishing source of both sequence and structural diversity. From intracellular viral RNA-host interfaces to interactions between the RNA genome and structural proteins in virus particles themselves, almost the entire viral lifecycle is accompanied by a myriad of RNA-protein interactions that are required to fulfill their replicative potential. It is therefore important to characterize such rich and dynamic collections of viral RNA-protein interactions to understand virus evolution and their adaptation to their hosts and environment. Recent advances in next-generation sequencing technologies have allowed the characterization of viral RNA-protein interactions, including both transient and conserved interactions, where molecular and structural approaches have fallen short. In this review, we will provide a methodological overview of the high-throughput techniques used to study viral RNA-protein interactions, their biochemical mechanisms, and how they evolved from classical methods as well as one another. We will discuss how different techniques have fueled virus research to characterize how viral RNA and proteins interact, both locally and on a global scale. Finally, we will present examples on how these techniques influence the studies of clinically important pathogens such as HIV-1 and SARS-CoV-2.

A Ribo-Seq Method to Study Genome-Wide Translational Regulation in Plants

Protein production from mRNA is one of the fundamental molecular processes in a cell. Accurate genome-wide information on the levels of translation and ribosome distribution on mRNA can be gathered by carrying out ribosome footprinting, aka Ribo-seq. Herein, we present a detailed protocol describing the construction of parallel Ribo-seq and RNA-seq libraries from Arabidopsis seedlings treated with the plant hormone auxin. The improved protocol for ribosome footprint library generation can be easily adapted to analyzing the effects on translation of genetic perturbations and various abiotic and biotic factors to shed the much-needed light on translational regulation in plants.

An evolutionarily conserved stop codon enrichment at the 5' ends of mammalian piRNAs

PIWI-interacting RNAs (piRNAs) are small RNAs required to recognize and silence transposable elements. The 5' ends of mature piRNAs are defined through cleavage of long precursor transcripts, primarily by Zucchini (Zuc). Zuc-dependent cleavage typically occurs immediately upstream of a uridine. However, Zuc lacks sequence preference in vitro, pointing towards additional unknown specificity factors. Here, we examine murine piRNAs and reveal a strong and specific enrichment of three sequences (UAA, UAG, UGA)-corresponding to stop codons-at piRNA 5' ends. Stop codon sequences are also enriched immediately after piRNA processing intermediates, reflecting their Zuc-dependent tail-to-head arrangement. Further analyses reveal that a Zuc in vivo cleavage preference at four sequences (UAA, UAG, UGA, UAC) promotes 5' end stop codons. This observation is conserved across mammals and possibly further. Our work provides new insights into Zuc-dependent cleavage and may point to a previously unrecognized connection between piRNA biogenesis and the translational machinery.

Exploring the expanding universe of small RNAs

The world of small noncoding RNAs (sncRNAs) is ever-expanding, from small interfering RNA, microRNA and Piwi-interacting RNA to the recently emerging non-canonical sncRNAs derived from longer structured RNAs (for example, transfer, ribosomal, Y, small nucleolar, small nuclear and vault RNAs), showing distinct biogenesis and functional principles. Here we discuss recent tools for sncRNA identification, caveats in sncRNA expression analysis and emerging methods for direct sequencing of sncRNAs and systematic mapping of RNA modifications that are integral to their function.

Small RNA-Sequencing for Analysis of Circulating miRNAs: Benchmark Study

Small RNA-sequencing (RNA-Seq) is being increasingly used for profiling of circulating microRNAs (miRNAs), a new group of promising biomarkers. Unfortunately, small RNA-Seq protocols are prone to biases limiting quantification accuracy, which motivated development of several novel methods. Here, we present comparison of all small RNA-Seq library preparation approaches that are commercially available for quantification of miRNAs in biofluids. Using synthetic and human plasma samples, we compared performance of traditional two-adaptor ligation protocols (Lexogen, Norgen), as well as methods using randomized adaptors (NEXTflex), polyadenylation (SMARTer), circularization (RealSeq), capture probes (EdgeSeq), or unique molecular identifiers (QIAseq). There was no single protocol outperforming others across all metrics. Limited overlap of measured miRNA profiles was documented between methods largely owing to protocol-specific biases. Methods designed to minimize bias largely differ in their performance, and contributing factors were identified. Usage of unique molecular identifiers has rather negligible effect and, if designed incorrectly, can even introduce spurious results. Together, these results identify strengths and weaknesses of all current methods and provide guidelines for applications of small RNA-Seq in biomarker research.

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.

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.

Global approaches for profiling transcription initiation

Transcription start site (TSS) selection influences transcript stability and translation as well as protein sequence. Alternative TSS usage is pervasive in organismal development, is a major contributor to transcript isoform diversity in humans, and is frequently observed in human diseases including cancer. In this review, we discuss the breadth of techniques that have been used to globally profile TSSs and the resulting insights into gene regulation, as well as future prospects in this area of inquiry.

Maternally inherited piRNAs direct transient heterochromatin formation at active transposons during early Drosophila embryogenesis

The PIWI-interacting RNA (piRNA) pathway controls transposon expression in animal germ cells, thereby ensuring genome stability over generations. In Drosophila, piRNAs are intergenerationally inherited through the maternal lineage, and this has demonstrated importance in the specification of piRNA source loci and in silencing of I- and P-elements in the germ cells of daughters. Maternally inherited Piwi protein enters somatic nuclei in early embryos prior to zygotic genome activation and persists therein for roughly half of the time required to complete embryonic development. To investigate the role of the piRNA pathway in the embryonic soma, we created a conditionally unstable Piwi protein. This enabled maternally deposited Piwi to be cleared from newly laid embryos within 30 min and well ahead of the activation of zygotic transcription. Examination of RNA and protein profiles over time, and correlation with patterns of H3K9me3 deposition, suggests a role for maternally deposited Piwi in attenuating zygotic transposon expression in somatic cells of the developing embryo. In particular, robust deposition of piRNAs targeting roo, an element whose expression is mainly restricted to embryonic development, results in the deposition of transient heterochromatic marks at active roo insertions. We hypothesize that roo, an extremely successful mobile element, may have adopted a lifestyle of expression in the embryonic soma to evade silencing in germ cells.

A microRNA expression and regulatory element activity atlas of the mouse immune system

To better define the control of immune system regulation, we generated an atlas of microRNA (miRNA) expression from 63 mouse immune cell populations and connected these signatures with assay for transposase-accessible chromatin using sequencing (ATAC-seq), chromatin immunoprecipitation followed by sequencing (ChIP-seq) and nascent RNA profiles to establish a map of miRNA promoter and enhancer usage in immune cells. miRNA complexity was relatively low, with >90% of the miRNA compartment of each population comprising <75 miRNAs; however, each cell type had a unique miRNA signature. Integration of miRNA expression with chromatin accessibility revealed putative regulatory elements for differentially expressed miRNAs, including miR-21a, miR-146a and miR-223. The integrated maps suggest that many miRNAs utilize multiple promoters to reach high abundance and identified dominant and divergent miRNA regulatory elements between lineages and during development that may be used by clustered miRNAs, such as miR-99a/let-7c/miR-125b, to achieve distinct expression. These studies, with web-accessible data, help delineate the cis-regulatory elements controlling miRNA signatures of the immune system.

Methyltransferase-directed orthogonal tagging and sequencing of miRNAs and bacterial small RNAs

Background

Targeted installation of designer chemical moieties on biopolymers provides an orthogonal means for their visualisation, manipulation and sequence analysis. Although high-throughput RNA sequencing is a widely used method for transcriptome analysis, certain steps, such as 3′ adapter ligation in strand-specific RNA sequencing, remain challenging due to structure- and sequence-related biases introduced by RNA ligases, leading to misrepresentation of particular RNA species. Here, we remedy this limitation by adapting two RNA 2′-O-methyltransferases from the Hen1 family for orthogonal chemo-enzymatic click tethering of a 3′ sequencing adapter that supports cDNA production by reverse transcription of the tagged RNA.

Results

We showed that the ssRNA-specific DmHen1 and dsRNA-specific AtHEN1 can be used to efficiently append an oligonucleotide adapter to the 3′ end of target RNA for sequencing library preparation. Using this new chemo-enzymatic approach, we identified miRNAs and prokaryotic small non-coding sRNAs in probiotic Lactobacillus casei BL23. We found that compared to a reference conventional RNA library preparation, methyltransferase-Directed Orthogonal Tagging and RNA sequencing, mDOT-seq, avoids misdetection of unspecific highly-structured RNA species, thus providing better accuracy in identifying the groups of transcripts analysed. Our results suggest that mDOT-seq has the potential to advance analysis of eukaryotic and prokaryotic ssRNAs.

Conclusions

Our findings provide a valuable resource for studies of the RNA-centred regulatory networks in Lactobacilli and pave the way to developing novel transcriptome and epitranscriptome profiling approaches in vitro and inside living cells. As RNA methyltransferases share the structure of the AdoMet-binding domain and several specific cofactor binding features, the basic principles of our approach could be easily translated to other AdoMet-dependent enzymes for the development of modification-specific RNA-seq techniques.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01053-w.

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.

CPA-seq reveals small ncRNAs with methylated nucleosides and diverse termini

High-throughput sequencing reveals the complex landscape of small noncoding RNAs (sRNAs). However, it is limited by requiring 5'-monophosphate and 3'-hydroxyl in RNAs for adapter ligation and hindered by methylated nucleosides that interfere with reverse transcription. Here we develop Cap-Clip acid pyrophosphatase (Cap-Clip), T4 polynucleotide kinase (PNK) and AlkB/AlkB(D135S)-facilitated small ncRNA sequencing (CPA-seq) to detect and quantify sRNAs with terminus multiplicities and nucleoside methylations. CPA-seq identified a large number of previously undetected sRNAs. Comparison of sRNAs with or without AlkB/AlkB(D135S) treatment reveals nucleoside methylations on sRNAs. Using CPA-seq, we profiled the sRNA transcriptomes (sRNomes) of nine mouse tissues and reported the extensive tissue-specific differences of sRNAs. We also observed the transition of sRNomes during hepatic reprogramming. Knockdown of mesenchymal stem cell-enriched U1-5' snsRNA promoted hepatic reprogramming. CPA-seq is a powerful tool with high sensitivity and specificity for profiling sRNAs with methylated nucleosides and diverse termini.

Quantitative mapping of the cellular small RNA landscape with AQRNA-seq

Current next-generation RNA sequencing methods do not provide accurate quantification of small RNAs within a sample due to sequence-dependent biases in capture, ligation, and amplification during library preparation. We present a method, Absolute Quantification (AQ) RNA-seq, that minimizes biases and provides a direct, linear correlation between sequencing read count and copy number for all small RNAs in a sample. Library preparation and data processing were optimized and validated using a 963-member miRNA reference library, oligonucleotide standards of varying lengths, and northern blots. Application of AQRNA-seq to a panel of human cancer cells revealed >800 detectable miRNAs that varied during cancer progression, while application to bacterial tRNA pools, with the challenges of secondary structure and abundant modifications, revealed 80-fold variation in tRNA isoacceptor levels, stress-induced site-specific tRNA fragmentation, quantitative modification maps, and evidence for stress-induced tRNA-driven codon-biased translation. AQRNA-seq thus provides a versatile means to quantitatively map the small RNA landscape in cells.

A novel platform for discovery of differentially expressed microRNAs in patients with repeated implantation failure

OBJECTIVE

To identify predictor microRNAs (miRNAs) from patients with repeated implantation failure (RIF).

DESIGN

Systemic analysis of miRNA profiles from the endometrium of patients undergoing in vitro fertilization (IVF).

SETTING

University research institute, private IVF center, and molecular testing laboratory.

PATIENT(S)

Twenty five infertile patients in the discovery cohort and 11 patients in the validation cohort.

INTERVENTIONS(S)

None.

MAIN OUTCOME MEASURE(S)

A signature set of miRNA associated with the risk of RIF.

RESULT(S)

We designed a reproductive disease-related PanelChip to access endometrium miRNA profiles in patients undergoing IVF. Three major miRNA signatures, including hsa-miR-20b-5p, hsa-miR-155-5p, and hsa-miR-718, were identified using infinite combination signature search algorithm analysis from 25 patients in the discovery cohort undergoing IVF. These miRNAs were used as biomarkers in the validation cohort of 11 patients. Finally, the 3-miRNA signature was capable of predicting patients with RIF with an accuracy >90%.

CONCLUSION(S)

Our findings indicated that specific endometrial miRNAs can be applied as diagnostic biomarkers to predict RIF. Such information will definitely help to increase the success rate of implantation practice.

A Small RNA-Seq Protocol with Less Bias and Improved Capture of 2'-O-Methyl RNAs

The study of small RNAs (sRNAs) by next-generation sequencing (NGS) is challenged by bias issues during library preparation. Several types of sRNAs such as plant microRNAs (miRNAs) carry a 2'-O-methyl (2'-OMe) modification at their 3' terminal nucleotide. This modification adds another level of difficulty as it inhibits 3' adapter ligation. We previously demonstrated that modified versions of the "TruSeq (TS)" protocol have less bias and an improved detection of 2'-OMe RNAs. Here we describe in detail protocol "TS5," which showed the best overall performance. We also provide guidelines for bioinformatics analysis of the sequencing data.

Structure-function analysis of microRNA 3'-end trimming by Nibbler

Nibbler (Nbr) is a 3'-to-5' exoribonuclease whose catalytic 3'-end trimming activity impacts microRNA (miRNA) and PIWI-interacting RNA (piRNA) biogenesis. Here, we report on structural and functional studies to decipher the contributions of Nbr's N-terminal domain (NTD) and exonucleolytic domain (EXO) in miRNA 3'-end trimming. We have solved the crystal structures of the NTD core and EXO domains of Nbr, both in the apo-state. The NTD-core domain of Aedes aegypti Nbr adopts a HEAT-like repeat scaffold with basic patches constituting an RNA-binding surface exhibiting a preference for binding double-strand RNA (dsRNA) over single-strand RNA (ssRNA). Structure-guided functional assays in Drosophila S2 cells confirmed a principal role of the NTD in exonucleolytic miRNA trimming, which depends on basic surface patches. Gain-of-function experiments revealed a potential role of the NTD in recruiting Nbr to Argonaute-bound small RNA substrates. The EXO domain of A. aegypti and Drosophila melanogaster Nbr adopt a mixed α/β-scaffold with a deep pocket lined by a DEDDy catalytic cleavage motif. We demonstrate that Nbr's EXO domain exhibits Mn2+-dependent ssRNA-specific 3'-to-5' exoribonuclease activity. Modeling of a 3' terminal Uridine into the catalytic pocket of Nbr EXO indicates that 2'-O-methylation of the 3'-U would result in a steric clash with a tryptophan side chain, suggesting that 2'-O-methylation protects small RNAs from Nbr-mediated trimming. Overall, our data establish that Nbr requires its NTD as a substrate recruitment platform to execute exonucleolytic miRNA maturation, catalyzed by the ribonuclease EXO domain.

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

Small RNAs are important regulators of gene expression and are involved in human development and disease. Next generation sequencing (NGS) allows for scalable, genome-wide studies of small RNA; however, current methods are challenged by low sensitivity and high bias, limiting their ability to capture an accurate representation of the cellular small RNA population. Several studies have shown that this bias primarily arises during the ligation of single-strand adapters during library preparation, and that this ligation bias is magnified by 2′-O-methyl modifications (2′OMe) on the 3′ terminal nucleotide. In this study, we developed a novel library preparation process using randomized splint ligation with a cleavable adapter, a design which resolves previous challenges associated with this ligation strategy. We show that a randomized splint ligation based workflow can reduce bias and increase the sensitivity of small RNA sequencing for a wide variety of small RNAs, including microRNA (miRNA) and tRNA fragments as well as 2′OMe modified RNA, including Piwi-interacting RNA and plant miRNA. Finally, we demonstrate that this workflow detects more differentially expressed miRNA between tumorous and matched normal tissues. Overall, this library preparation process allows for highly accurate small RNA sequencing and will enable studies of 2′OMe modified RNA with new levels of detail.

A universal method for the rapid isolation of all known classes of functional silencing small RNAs

Diverse classes of silencing small (s)RNAs operate via ARGONAUTE-family proteins within RNA-induced-silencing-complexes (RISCs). Here, we have streamlined various embodiments of a Q-sepharose-based RISC-purification method that relies on conserved biochemical properties of all ARGONAUTEs. We show, in multiple benchmarking assays, that the resulting 15-min benchtop extraction procedure allows simultaneous purification of all known classes of RISC-associated sRNAs without prior knowledge of the samples-intrinsic ARGONAUTE repertoires. Optimized under a user-friendly format, the method – coined ‘TraPR’ for Trans-kingdom, rapid, affordable Purification of RISCs – operates irrespectively of the organism, tissue, cell type or bio-fluid of interest, and scales to minute amounts of input material. The method is highly suited for direct profiling of silencing sRNAs, with TraPR-generated sequencing libraries outperforming those obtained via gold-standard procedures that require immunoprecipitations and/or lengthy polyacrylamide gel-selection. TraPR considerably improves the quality and consistency of silencing sRNA sample preparation including from notoriously difficult-to-handle tissues/bio-fluids such as starchy storage roots or mammalian plasma, and regardless of RNA contaminants or RNA degradation status of samples.

Comparison of Reproducibility, Accuracy, Sensitivity, and Specificity of miRNA Quantification Platforms

Given the increasing interest in their use as disease biomarkers, the establishment of reproducible, accurate, sensitive, and specific platforms for microRNA (miRNA) quantification in biofluids is of high priority. We compare four platforms for these characteristics: small RNA sequencing (RNA-seq), FirePlex, EdgeSeq, and nCounter. For a pool of synthetic miRNAs, coefficients of variation for technical replicates are lower for EdgeSeq (6.9%) and RNA-seq (8.2%) than for FirePlex (22.4%); nCounter replicates are not performed. Receiver operating characteristic analysis for distinguishing present versus absent miRNAs shows small RNA-seq (area under curve 0.99) is superior to EdgeSeq (0.97), nCounter (0.94), and FirePlex (0.81). Expected differences in expression of placenta-associated miRNAs in plasma from pregnant and non-pregnant women are observed with RNA-seq and EdgeSeq, but not FirePlex or nCounter. These results indicate that differences in performance among miRNA profiling platforms impact ability to detect biological differences among samples and thus their relative utility for research and clinical use.

Using pools of synthetic RNA oligonucleotides and standardized extracellular RNA samples, Godoy et al. compare small RNA sequencing to three targeted miRNA quantification platforms to evaluate reproducibility, bias, specificity and sensitivity, and accuracy. Each platform has strengths and limitations important to consider for biomarker discovery, clinical validation, and broad clinical use.

Flåm S, Leithaug M, Dalland M, Rayner S, Lie BA, Gilfillan GD. An extension to: Systematic assessment of commercially available low-input miRNA library preparation kits

High-throughput sequencing has emerged as the favoured method to study microRNA (miRNA) expression, but biases introduced during library preparation have been reported. We recently compared the performance (sensitivity, reliability, titration response and differential expression) of six commercially-available kits on synthetic miRNAs and human RNA, where library preparation was performed by the vendors. We hereby supplement this study with data from two further commonly used kits (NEBNext, NEXTflex) whose manufacturers initially declined to participate. NEXTflex demonstrated the highest sensitivity, which may reflect its use of partially-randomized adapter sequences, but overall performance was lower than the QIAseq and TailorMix kits. NEBNext showed intermediate performance. We reaffirm that biases are kit specific, complicating the comparison of miRNA datasets generated using different kits.

Identification of a potential non-coding RNA biomarker signature for amyotrophic lateral sclerosis

Objective biomarkers for the clinically heterogeneous adult-onset neurodegenerative disorder amyotrophic lateral sclerosis are crucial to facilitate assessing emerging therapeutics and improve the diagnostic pathway in what is a clinically heterogeneous syndrome. With non-coding RNA transcripts including microRNA, piwi-RNA and transfer RNA present in human biofluids, we sought to identify whether non-coding RNA in serum could be biomarkers for amyotrophic lateral sclerosis. Serum samples from our Oxford Study for Biomarkers in motor neurone disease/amyotrophic lateral sclerosis discovery cohort of amyotrophic lateral sclerosis patients (n = 48), disease mimics (n = 16) and age- and sex-matched healthy controls (n = 24) were profiled for non-coding RNA expression using RNA-sequencing, which showed a wide range of non-coding RNA to be dysregulated. We confirmed significant alterations with reverse transcription-quantitative PCR in the expression of hsa-miR-16-5p, hsa-miR-21-5p, hsa-miR-92a-3p, hsa-piR-33151, TRV-AAC4-1.1 and TRA-AGC6-1.1. Furthermore, hsa-miR-206, a previously identified amyotrophic lateral sclerosis biomarker, showed a binary-like pattern of expression in our samples. Using the expression of these non-coding RNA, we were able to discriminate amyotrophic lateral sclerosis samples from healthy controls in our discovery cohort using a random forest analysis with 93.7% accuracy with promise in predicting progression rate of patients. Importantly, cross-validation of this novel signature using a new geographically distinct cohort of samples from the United Kingdom and Germany with both amyotrophic lateral sclerosis and control samples (n = 156) yielded an accuracy of 73.9%. The high prediction accuracy of this non-coding RNA-based biomarker signature, even across heterogeneous cohorts, demonstrates the strength of our approach as a novel platform to identify and stratify amyotrophic lateral sclerosis patients.

High-Throughput Identification of Adapters in Single-Read Sequencing Data

Sequencing datasets available in public repositories are already high in number, and their growth is exponential. Raw sequencing data files constitute a substantial portion of these data, and they need to be pre-processed for any downstream analyses. The removal of adapter sequences is the first essential step. Tools available for the automated detection of adapters in single-read sequencing protocol datasets have certain limitations. To explore these datasets, one needs to retrieve the information on adapter sequences from the methods sections of appropriate research articles. This can be time-consuming in metadata analyses. Moreover, not all research articles provide the information on adapter sequences. We have developed adapt_find, a tool that automates the process of adapter sequences identification in raw single-read sequencing datasets. We have verified adapt_find through testing a number of publicly available datasets. adapt_find secures a robust, reliable and high-throughput process across different sequencing technologies and various adapter designs. It does not need prior knowledge of the adapter sequences. We also produced associated tools: random_mer, for the detection of random N bases either on one or both termini of the reads, and fastqc_parser, for consolidating the results from FASTQC outputs. Together, this is a valuable tool set for metadata analyses on multiple sequencing datasets.

Simple and efficient profiling of transcription initiation and transcript levels with STRIPE-seq

Accurate mapping of transcription start sites (TSSs) is key for understanding transcriptional regulation. However, current protocols for genome-wide TSS profiling are laborious and/or expensive. We present Survey of TRanscription Initiation at Promoter Elements with high-throughput sequencing (STRIPE-seq), a simple, rapid, and cost-effective protocol for sequencing capped RNA 5' ends from as little as 50 ng total RNA. Including depletion of uncapped RNA and reaction cleanups, a STRIPE-seq library can be constructed in about 5 h. We show application of STRIPE-seq to TSS profiling in yeast and human cells and show that it can also be effectively used for quantification of transcript levels and analysis of differential gene expression. In conjunction with our ready-to-use computational workflows, STRIPE-seq is a straightforward, efficient means by which to probe the landscape of transcriptional initiation.

Increased Hematopoietic Extracellular RNAs and Vesicles in the Lung during Allergic Airway Responses

Extracellular RNAs (exRNAs) can be released by numerous cell types in vitro, are often protected within vesicles, and can modify recipient cell function. To determine how the composition and cellular sources of exRNAs and the extracellular vesicles (EVs) that carry them change in vivo during tissue inflammation, we analyzed bronchoalveolar lavage fluid (BALF) from mice before and after lung allergen challenge. In the lung, extracellular microRNAs (ex-miRNAs) had a composition that was highly correlated with airway-lining epithelium. Using cell type-specific membrane tagging and single vesicle flow, we also found that 80% of detected vesicles were of epithelial origin. After the induction of allergic airway inflammation, miRNAs selectively expressed by immune cells, including miR-223 and miR-142a, increased and hematopoietic-cell-derived EVs also increased >2-fold. These data demonstrate that infiltrating immune cells release ex-miRNAs and EVs in inflamed tissues to alter the local extracellular environment.

Small-RNA-mediated transgenerational silencing of histone genes impairs fertility in piRNA mutants

PIWI-interacting RNAs (piRNAs) promote fertility in many animals. However, whether this is due to their conserved role in repressing repetitive elements (REs) remains unclear. Here, we show that the progressive loss of fertility in Caenorhabditis elegans lacking piRNAs is not caused by derepression of REs or other piRNA targets but, rather, is mediated by epigenetic silencing of all of the replicative histone genes. In the absence of piRNAs, downstream components of the piRNA pathway relocalize from germ granules and piRNA targets to histone mRNAs to synthesize antisense small RNAs (sRNAs) and induce transgenerational silencing. Removal of the downstream components of the piRNA pathway restores histone mRNA expression and fertility in piRNA mutants, and the inheritance of histone sRNAs in wild-type worms adversely affects their fertility for multiple generations. We conclude that sRNA-mediated silencing of histone genes impairs the fertility of piRNA mutants and may serve to maintain piRNAs across evolution.

Highly efficient single-stranded DNA ligation technique improves low-input whole-genome bisulfite sequencing by post-bisulfite adaptor tagging

Whole-genome bisulfite sequencing (WGBS) is the current gold standard of methylome analysis. Post-bisulfite adaptor tagging (PBAT) is an increasingly popular WGBS protocol because of high sensitivity and low bias. PBAT originally relied on two rounds of random priming for adaptor-tagging of single-stranded DNA (ssDNA) to attain high efficiency but at a cost of library insert length. To overcome this limitation, we developed terminal deoxyribonucleotidyl transferase (TdT)-assisted adenylate connector-mediated ssDNA (TACS) ligation as an alternative to random priming. In this method, TdT attaches adenylates to the 3′-end of input ssDNA, which are then utilized by RNA ligase as an efficient connector to the ssDNA adaptor. A protocol that uses TACS ligation instead of the second random priming step substantially increased the lengths of PBAT library fragments. Moreover, we devised a dual-library strategy that splits the input DNA to prepare two libraries with reciprocal adaptor polarity, combining them prior to sequencing. This strategy ensured an ideal base–color balance to eliminate the need for DNA spike-in for color compensation, further improving the throughput and quality of WGBS. Adopting the above strategies to the HiSeq X Ten and NovaSeq 6000 platforms, we established a cost-effective, high-quality WGBS, which should accelerate various methylome analyses.

Assessment of methods for serum extracellular vesicle small RNA sequencing to support biomarker development

Extracellular vesicles (EVs) have great potential as a source for clinically relevant biomarkers since they can be readily isolated from biofluids and carry microRNA (miRNA), mRNA, and proteins that can reflect disease status. However, the biological and technical variability of EV content is unknown making comparisons between healthy subjects and patients difficult to interpret. In this study, we sought to establish a laboratory and bioinformatics analysis pipeline to analyse the small RNA content within EVs from patient serum that could serve as biomarkers and to assess the biological and technical variability of EV RNA content in healthy individuals. We sequenced EV small RNA from multiple individuals (biological replicates) and sequenced multiple replicates per individual (technical replicates) using the Illumina Truseq protocol. We observed that the replicates of samples clustered by subject indicating that the biological variability (~95%) was greater than the technical variability (~0.50%). We observed that ~30% of the sequencing reads were miRNAs. We evaluated the technical parameters of sequencing by spiking the EV RNA preparation with a mix of synthetic small RNA and demonstrated a disconnect between input concentration of the spike-in RNA and sequencing read frequencies indicating that bias was introduced during library preparation. To determine whether there are differences between library preparation platforms, we compared the Truseq with the Nextflex protocol that had been designed to reduce bias in library preparation. While both methods were technically robust, the Nextflex protocol reduced the bias and exhibited a linear range across input concentrations of the synthetic spike-ins. Altogether, our results indicate that technical variability is much smaller than biological variability supporting the use of EV small RNAs as potential biomarkers. Our findings also indicate that the choice of library preparation method leads to artificial differences in the datasets generated invalidating the comparability of sequencing data across library preparation platforms.

Comparative Analysis of Biochemical Biases by Ligation- and Template-Switch-Based Small RNA Library Preparation Protocols

BACKGROUND

Small RNAs are key players in the regulation of gene expression and differentiation. However, many different classes of small RNAs (sRNAs) have been described with distinct biogenesis pathways and, as a result, with different biochemical properties. To analyze sRNAs by deep sequencing, complementary DNA synthesis requires manipulation of the RNA molecule itself. Thus, enzymatic activities during library preparation bias the library content owing to biochemical criteria.

METHODS

We compared 4 different manipulations of RNA for library preparation: (a) a ligation-based procedure allowing only 5'-mono-phosphorylated RNA to enter the library, (b) a ligation-based procedure allowing additional 5'-triphosphates and Cap structures, (c) a ligation-independent, template-switch-based library preparation, and (d) a template-switch-based library preparation allowing 3'-phosphorylated RNAs to enter the library.

RESULTS

Our data show large differences between ligation-dependent and ligation-independent libraries in terms of their preference for individual sRNA classes such as microRNAs (miRNAs), Piwi-interacting RNAs (piRNAs), and transfer RNA fragments. Moreover, the miRNA composition is different between both procedures, and more microRNA isoforms (isomiRs) can be identified after pyrophosphatase treatment. piRNAs are enriched in template-switch libraries, and this procedure apparently includes more different RNA species.

CONCLUSIONS

Our data indicate that miRNAomics from both methods will hardly be comparable. Ligation-based libraries enrich for canonical miRNAs, which thus may be suitable methods for miRNAomics. Template-switch libraries contain increased numbers and different compositions of fragments and long RNAs. Following different interests for other small RNA species, ligation-independent libraries appear to show a more realistic sRNA landscape with lower bias against biochemical modifications.

Non-templated addition and template switching by Moloney murine leukemia virus (MMLV)-based reverse transcriptases co-occur and compete with each other

Single-cell RNA-Seq (scRNA-Seq) has led to an unprecedented understanding of gene expression and regulation in individual cells. Many scRNA-Seq approaches rely upon the template switching property of Moloney murine leukemia virus (MMLV)-type reverse transcriptases. Template switching is believed to happen in a sequential process involving nontemplated addition of three protruding nucleotides (+CCC) to the 3′-end of the nascent cDNA, which can then anneal to the matching rGrGrG 3′-end of the template-switching oligo (TSO), allowing the reverse transcriptase (RT) to switch templates and continue copying the TSO sequence. In this study, we present a detailed analysis of template switching biases with respect to the RNA template, specifically of the role of the sequence and nature of its 5′-end (capped versus noncapped) in these biases. Our findings confirmed that the presence of a 5′-m⁷G cap enhances template switching efficiency. We also profiled the composition of the nontemplated addition in the absence of TSO and observed that the 5′-end of RNA template influences the terminal transferase activity of the RT. Furthermore, we found that designing new TSOs that pair with the most common nontemplated additions did little to improve template switching efficiency. Our results provide evidence suggesting that, in contrast to the current understanding of the template switching process, nontemplated addition and template switching are concurrent and competing processes.