RNA-ligase-dependent biases in miRNA representation in deep-sequenced small RNA cDNA libraries

Guidelines for mitochondrial RNA analysis

Mitochondria are the energy-producing organelles of mammalian cells with critical involvement in metabolism and signaling. Studying their regulation in pathological conditions may lead to the discovery of novel drugs to treat, for instance, cardiovascular or neurological diseases, which affect high-energy-consuming cells such as cardiomyocytes, hepatocytes, or neurons. Mitochondria possess both protein-coding and noncoding RNAs, such as microRNAs, long noncoding RNAs, circular RNAs, and piwi-interacting RNAs, encoded by the mitochondria or the nuclear genome. Mitochondrial RNAs are involved in anterograde-retrograde communication between the nucleus and mitochondria and play an important role in physiological and pathological conditions. Despite accumulating evidence on the presence and biogenesis of mitochondrial RNAs, their study continues to pose significant challenges. Currently, there are no standardized protocols and guidelines to conduct deep functional characterization and expression profiling of mitochondrial RNAs. To overcome major obstacles in this emerging field, the EU-CardioRNA and AtheroNET COST Action networks summarize currently available techniques and emphasize critical points that may constitute sources of variability and explain discrepancies between published results. Standardized methods and adherence to guidelines to quantify and study mitochondrial RNAs in normal and disease states will improve research outputs, their reproducibility, and translation potential to clinical application.

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.

Small RNA structural biochemistry in a post-sequencing era

High-throughput sequencing has had an enormous impact on small RNA research during the past decade. However, sequencing only offers a one-dimensional view of the transcriptome and is often highly biased. Additionally, the 'sequence, map and annotate' approach, used widely in small RNA research, can lead to flawed interpretations of the data, lacking biological plausibility, due in part to database issues. Even in the absence of technical biases, the loss of three-dimensional information is a major limitation to understanding RNA stability, turnover and function. For example, noncoding RNA-derived fragments seem to exist mainly as dimers, tetramers or as nicked forms of their parental RNAs, contrary to widespread assumptions. In this perspective, we will discuss main sources of bias during small RNA-sequencing, present several useful bias-reducing strategies and provide guidance on the interpretation of small RNA-sequencing results, with emphasis on RNA fragmentomics. As sequencing offers a one-dimensional projection of a four-dimensional reality, prior structure-level knowledge is often needed to make sense of the data. Consequently, while less-biased sequencing methods are welcomed, integration of orthologous experimental techniques is also strongly recommended.

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.

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.

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.

Systematic analysis of cotranscriptional RNA folding using transcription elongation complex display

RNA can fold into structures that mediate diverse cellular functions. Understanding how RNA primary sequence directs the formation of functional structures requires methods that can comprehensively assess how changes in an RNA sequence affect its structure and function. Here we have developed a platform for performing high-throughput cotranscriptional RNA biochemical assays, called Transcription Elongation Complex display (TECdisplay). TECdisplay measures RNA function by fractionating a TEC library based on the activity of cotranscriptionally displayed nascent RNA. In this way, RNA function is measured as the distribution of template DNA molecules between fractions of the transcription reaction. This approach circumvents typical RNA sequencing library preparation steps that can cause technical bias. We used TECdisplay to characterize the transcription antitermination activity of 32,768 variants of the Clostridium beijerinckii pfl ZTP riboswitch designed to perturb steps within its cotranscriptional folding pathway. Our findings establish TECdisplay as an accessible platform for high-throughput RNA biochemical assays.

PRO-seq: Precise Mapping of Engaged RNA Pol II at Single-Nucleotide Resolution

Gene regulation is dependent on the production of mRNAs and a repertoire of non-coding RNAs by RNA polymerase II (RNAPII). Precision run-on sequencing (PRO-seq) maps the position of engaged RNAPII complexes at single-nucleotide resolution and can reveal direct targets of regulation, locations of enhancers, and transcription mechanisms that are difficult or impossible to measure by analysis of total cellular RNA. Briefly, this method first involves permeabilizing cells with mild detergents to remove intracellular NTPs and halt transcription. Transcription is then resumed in the presence of biotin-NTPs and sarkosyl to allow transcriptional incorporation of a single biotinylated NTP by RNAPII. The biotin moiety is then bound to streptavidin beads to stringently enrich for nascent RNAs. Sequencing libraries are then generated such that the first base read corresponds to the 3' end of the nascent transcript. Here, we describe our current protocol for generating PRO-seq libraries from metazoan cells, including adaptations of previously published protocols to incorporate unique molecular identifiers, reduce ligation bias, and improve library yields. Additional commentary describes quality control and processing of PRO-seq data and references for more advanced downstream analysis such as gene and enhancer identification. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Cell permeabilization for PRO-seq Basic Protocol 2: Construction of PRO-seq libraries Support Protocol: Adenylation of 3' adapter.

Start codon-associated ribosomal frameshifting mediates nutrient stress adaptation

A translating ribosome is typically thought to follow the reading frame defined by the selected start codon. Using super-resolution ribosome profiling, here we report pervasive out-of-frame translation immediately from the start codon. Start codon-associated ribosomal frameshifting (SCARF) stems from the slippage of ribosomes during the transition from initiation to elongation. Using a massively paralleled reporter assay, we uncovered sequence elements acting as SCARF enhancers or repressors, implying that start codon recognition is coupled with reading frame fidelity. This finding explains thousands of mass spectrometry spectra that are unannotated in the human proteome. Mechanistically, we find that the eukaryotic initiation factor 5B (eIF5B) maintains the reading frame fidelity by stabilizing initiating ribosomes. Intriguingly, amino acid starvation induces SCARF by proteasomal degradation of eIF5B. The stress-induced SCARF protects cells from starvation by enabling amino acid recycling and selective mRNA translation. Our findings illustrate a beneficial effect of translational 'noise' in nutrient stress adaptation.

Optimized infrared photoactivatable ribonucleoside-enhanced crosslinking and immunoprecipitation (IR-PAR-CLIP) protocol identifies novel IGF2BP3-interacting RNAs in colon cancer cells

The conserved family of RNA-binding proteins (RBPs), IGF2BPs, plays an essential role in posttranscriptional regulation controlling mRNA stability, localization, and translation. Mammalian cells express three isoforms of IGF2BPs: IGF2BP1-3. IGF2BP3 is highly overexpressed in cancer cells, and its expression correlates with a poor prognosis in various tumors. Therefore, revealing its target RNAs with high specificity in healthy tissues and in cancer cells is of crucial importance. Photoactivatable-ribonucleoside-enhanced crosslinking and immunoprecipitation (PAR-CLIP) identifies the binding sites of RBPs on their target RNAs at nucleotide resolution in a transcriptome-wide manner. Here, we optimized the PAR-CLIP protocol to study RNA targets of endogenous IGF2BP3 in a human colorectal carcinoma cell line. To this end, we first established an immunoprecipitation protocol to obtain highly pure endogenous IGF2BP3-RNA complexes. Second, we modified the protocol to use highly sensitive infrared (IR) fluorescent dyes instead of radioactive probes to visualize IGF2BP3-crosslinked RNAs. We named the modified method "IR-PAR-CLIP." Third, we compared RNase cleavage conditions and found that sequence preferences of the RNases impact the number of the identified IGF2BP3 targets and introduce a systematic bias in the identified RNA motifs. Fourth, we adapted the single adapter circular ligation approach to increase the efficiency in library preparation. The optimized IR-PAR-CLIP protocol revealed novel RNA targets of IGF2BP3 in a human colorectal carcinoma cell line. We anticipate that our IR-PAR-CLIP approach provides a framework for studies of other RBPs.

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.

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.

Start codon-associated ribosomal frameshifting mediates nutrient stress adaptation

A translating ribosome is typically thought to follow the reading frame defined by the selected start codon. Using super-resolution ribosome profiling, here we report pervasive out-of-frame translation immediately from the start codon. The start codon-associated ribosome frameshifting (SCARF) stems from the slippage of ribosomes during the transition from initiation to elongation. Using a massively paralleled reporter assay, we uncovered sequence elements acting as SCARF enhancers or repressors, implying that start codon recognition is coupled with reading frame fidelity. This finding explains thousands of mass spectrometry spectra unannotated from human proteome. Mechanistically, we find that the eukaryotic initiation factor 5B (eIF5B) maintains the reading frame fidelity by stabilizing initiating ribosomes. Intriguingly, amino acid starvation induces SCARF by proteasomal degradation of eIF5B. The stress-induced SCARF protects cells from starvation by enabling amino acid recycling and selective mRNA translation. Our findings illustrate a beneficial effect of translational "noise" in nutrient stress adaptation.

APE1 controls DICER1 expression in NSCLC through miR-33a and miR-130b

Increasing evidence suggests different, not completely understood roles of microRNA biogenesis in the development and progression of lung cancer. The overexpression of the DNA repair protein apurinic/apyrimidinic endodeoxyribonuclease 1 (APE1) is an important cause of poor chemotherapeutic response in lung cancer and its involvement in onco-miRNAs biogenesis has been recently described. Whether APE1 regulates miRNAs acting as prognostic biomarkers of lung cancer has not been investigated, yet. In this study, we analyzed miRNAs differential expression upon APE1 depletion in the A549 lung cancer cell line using high-throughput methods. We defined a signature of 13 miRNAs that strongly correlate with APE1 expression in human lung cancer: miR-1246, miR-4488, miR-24, miR-183, miR-660, miR-130b, miR-543, miR-200c, miR-376c, miR-218, miR-146a, miR-92b and miR-33a. Functional enrichment analysis of this signature revealed its biological relevance in cancer cell proliferation and survival. We validated DICER1 as a direct functional target of the APE1-regulated miRNA-33a-5p and miR-130b-3p. Importantly, IHC analyses of different human tumors confirmed a negative correlation existing between APE1 and Dicer1 protein levels. DICER1 downregulation represents a prognostic marker of cancer development but the mechanisms at the basis of this phenomenon are still completely unknown. Our findings, suggesting that APE1 modulates DICER1 expression via miR-33a and miR-130b, reveal new mechanistic insights on DICER1 regulation, which are of relevance in lung cancer chemoresistance and cancer invasiveness.

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.

Topology preserving stratification of tissue neoplasticity using Deep Neural Maps and microRNA signatures

Background

Accurate cancer classification is essential for correct treatment selection and better prognostication. microRNAs (miRNAs) are small RNA molecules that negatively regulate gene expression, and their dyresgulation is a common disease mechanism in many cancers. Through a clearer understanding of miRNA dysregulation in cancer, improved mechanistic knowledge and better treatments can be sought.

Results

We present a topology-preserving deep learning framework to study miRNA dysregulation in cancer. Our study comprises miRNA expression profiles from 3685 cancer and non-cancer tissue samples and hierarchical annotations on organ and neoplasticity status. Using unsupervised learning, a two-dimensional topological map is trained to cluster similar tissue samples. Labelled samples are used after training to identify clustering accuracy in terms of tissue-of-origin and neoplasticity status. In addition, an approach using activation gradients is developed to determine the attention of the networks to miRNAs that drive the clustering. Using this deep learning framework, we classify the neoplasticity status of held-out test samples with an accuracy of 91.07%, the tissue-of-origin with 86.36%, and combined neoplasticity status and tissue-of-origin with an accuracy of 84.28%. The topological maps display the ability of miRNAs to recognize tissue types and neoplasticity status. Importantly, when our approach identifies samples that do not cluster well with their respective classes, activation gradients provide further insight in cancer subtypes or grades.

Conclusions

An unsupervised deep learning approach is developed for cancer classification and interpretation. This work provides an intuitive approach for understanding molecular properties of cancer and has significant potential for cancer classification and treatment selection.

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.

TGIRT-seq Protocol for the Comprehensive Profiling of Coding and Non-coding RNA Biotypes in Cellular, Extracellular Vesicle, and Plasma RNAs

High-throughput RNA sequencing (RNA-seq) has extraordinarily advanced our understanding of gene expression and disease etiology, and is a powerful tool for the identification of biomarkers in a wide range of organisms. However, most RNA-seq methods rely on retroviral reverse transcriptases (RTs), enzymes that have inherently low fidelity and processivity, to convert RNAs into cDNAs for sequencing. Here, we describe an RNA-seq protocol using Thermostable Group II Intron Reverse Transcriptases (TGIRTs), which have high fidelity, processivity, and strand-displacement activity, as well as a proficient template-switching activity that enables efficient and seamless RNA-seq adapter addition. By combining these activities, TGIRT-seq enables the simultaneous profiling of all RNA biotypes from small amounts of starting material, with superior RNA-seq metrics, and unprecedented ability to sequence structured RNAs. The TGIRT-seq protocol for Illumina sequencing consists of three steps: (i) addition of a 3' RNA-seq adapter, coupled to the initiation of cDNA synthesis at the 3' end of a target RNA, via template switching from a synthetic adapter RNA/DNA starter duplex; (ii) addition of a 5' RNA-seq adapter, by using thermostable 5' App DNA/RNA ligase to ligate an adapter oligonucleotide to the 3' end of the completed cDNA; (iii) minimal PCR amplification, to add capture sites and indices for Illumina sequencing. TGIRT-seq for the Illumina sequencing platform has been used for comprehensive profiling of coding and non-coding RNAs in ribodepleted, chemically fragmented cellular RNAs, and for the analysis of intact (non-chemically fragmented) cellular, extracellular vesicle (EV), and plasma RNAs, where it yields continuous full-length end-to-end sequences of structured small non-coding RNAs (sncRNAs), including tRNAs, snoRNAs, snRNAs, pre-miRNAs, and full-length excised linear intron (FLEXI) RNAs. Graphic abstract: Figure 1.Overview of the TGIRT-seq protocol for Illumina sequencing.Major steps are: (1) Template switching from a synthetic R2 RNA/R2R DNA starter duplex with a 1-nt 3' DNA overhang (a mixture of A, C, G, and T residues, denoted N) that base pairs to the 3' nucleotide of a target RNA, and upon initiating reverse transcription by adding dNTPs, seamlessly links an R2R adapter to the 5' end of the resulting cDNA; (2) Ligation of an R1R adapter to the 3' end of the completed cDNA; and (3) Minimal PCR amplification with primers that add Illumina capture sites (P5 and P7) and barcode sequences (indices 5 and 7). The index 7 barcode is required, while the index 5 barcode is optional, to provide unique dual indices (UDIs).

Cellular abundance shapes function in piRNA-guided genome defense

Defense against genome invaders universally relies on RNA-guided immunity. Prokaryotic CRISPR-Cas and eukaryotic RNA interference pathways recognize targets by complementary base-pairing, which places the sequences of their guide RNAs at the center of self/nonself discrimination. Here, we explore the sequence space of PIWI-interacting RNAs (piRNAs), the genome defense of animals, and establish functional priority among individual sequences. Our results reveal that only the topmost abundant piRNAs are commonly present in every cell, whereas rare sequences generate cell-to-cell diversity in flies and mice. We identify a skewed distribution of sequence abundance as a hallmark of piRNA populations and show that quantitative differences of more than a 1000-fold are established by conserved mechanisms of biogenesis. Finally, our genomics analyses and direct reporter assays reveal that abundance determines function in piRNA-guided genome defense. Taken together, we identify an effective sequence space and untangle two classes of piRNAs that differ in complexity and function. The first class represents the topmost abundant sequences and drives silencing of genomic parasites. The second class sparsely covers an enormous sequence space. These rare piRNAs cannot function in every cell, every individual, or every generation but create diversity with potential for adaptation in the ongoing arms race with genome invaders.

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.

Molecular markers of risk of subsequent invasive breast cancer in women with ductal carcinoma in situ: protocol for a population-based cohort study

INTRODUCTION

Ductal carcinoma in situ (DCIS) of the breast is a non-obligate precursor of invasive breast cancer (IBC). Many DCIS patients are either undertreated or overtreated. The overarching goal of the study described here is to facilitate detection of patients with DCIS at risk of IBC development. Here, we propose to use risk factor data and formalin-fixed paraffin-embedded (FFPE) DCIS tissue from a large, ethnically diverse, population-based cohort of 8175 women with a first diagnosis of DCIS and followed for subsequent IBC to: identify/validate miRNA expression changes in DCIS tissue associated with risk of subsequent IBC; evaluate ipsilateral IBC risk in association with two previously identified marker sets (triple immunopositivity for p16, COX-2, Ki67; Oncotype DX Breast DCIS score); examine the association of risk factor data with IBC risk.

METHODS AND ANALYSIS

We are conducting a series of case-control studies nested within the cohort. Cases are women with DCIS who developed subsequent IBC; controls (2/case) are matched to cases on calendar year of and age at DCIS diagnosis. We project 485 cases/970 controls in the aim focused on risk factors. We estimate obtaining FFPE tissue for 320 cases/640 controls for the aim focused on miRNAs; of these, 173 cases/346 controls will be included in the aim focused on p16, COX-2 and Ki67 immunopositivity, and of the latter, 156 case-control pairs will be included in the aim focused on the Oncotype DX Breast DCIS score®. Multivariate conditional logistic regression will be used for statistical analyses.

ETHICS AND DISSEMINATION

Ethics approval was obtained from the Institutional Review Boards of Albert Einstein College of Medicine (IRB 2014-3611), Kaiser Permanente Colorado, Kaiser Permanente Hawaii, Henry Ford Health System, Mayo Clinic, Marshfield Clinic Research Institute and Hackensack Meridian Health, and from Lifespan Research Protection Office. The study results will be presented at meetings and published in peer-reviewed journals.

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.

Statistical Assessment of Depth Normalization for Small RNA Sequencing

PURPOSE

Methods for depth normalization have been assessed primarily with simulated data or cell-line–mixture data. There is a pressing need for benchmark data enabling a more realistic and objective assessment, especially in the context of small RNA sequencing.

METHODS

We collected a unique pair of microRNA sequencing data sets for the same set of tumor samples; one data set was collected with and the other without uniform handling and balanced design. The former provided a benchmark for evaluating evidence of differential expression and the latter served as a test bed for normalization. Next, we developed a data perturbation algorithm to simulate additional data set pairs. Last, we assembled a set of computational tools to visualize and quantify the assessment.

RESULTS

We validated the quality of the benchmark data and showed the need for normalization of the test data. For illustration, we applied the data and tools to assess the performance of 9 existing normalization methods. Among them, trimmed mean of M-values was a better scaling method, whereas the median and the upper quartiles were consistently the worst performers; one variation of remove unwanted variation had the best chance of capturing true positives but at the cost of increased false positives. In general, these methods were, at best, moderately helpful when the level of differential expression was extensive and asymmetric.

CONCLUSION

Our study (1) provides the much-needed benchmark data and computational tools for assessing depth normalization, (2) shows the dependence of normalization performance on the underlying pattern of differential expression, and (3) calls for continued research efforts to develop more effective normalization methods.

Inducible and reversible inhibition of miRNA-mediated gene repression in vivo

Although virtually all gene networks are predicted to be controlled by miRNAs, the contribution of this important layer of gene regulation to tissue homeostasis in adult animals remains unclear. Gain and loss-of-function experiments have provided key insights into the specific function of individual miRNAs, but effective genetic tools to study the functional consequences of global inhibition of miRNA activity in vivo are lacking. Here we report the generation and characterization of a genetically engineered mouse strain in which miRNA-mediated gene repression can be reversibly inhibited without affecting miRNA biogenesis or abundance. We demonstrate the usefulness of this strategy by investigating the consequences of acute inhibition of miRNA function in adult animals. We find that different tissues and organs respond differently to global loss of miRNA function. While miRNA-mediated gene repression is essential for the homeostasis of the heart and the skeletal muscle, it is largely dispensable in the majority of other organs. Even in tissues where it is not required for homeostasis, such as the intestine and hematopoietic system, miRNA activity can become essential during regeneration following acute injury. These data support a model where many metazoan tissues primarily rely on miRNA function to respond to potentially pathogenic events.

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

Current next-generation RNA-sequencing (RNA-seq) 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 RNA-sequencing (AQRNA-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 microRNA reference library, oligonucleotide standards of varying length, and RNA 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 transfer RNA 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.

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.

A non-radioactive, improved PAR-CLIP and small RNA cDNA library preparation protocol

Crosslinking and immunoprecipitation (CLIP) methods are powerful techniques to interrogate direct protein-RNA interactions and dissect posttranscriptional gene regulatory networks. One widely used CLIP variant is photoactivatable ribonucleoside enhanced CLIP (PAR-CLIP) that involves in vivo labeling of nascent RNAs with the photoreactive nucleosides 4-thiouridine (4SU) or 6-thioguanosine (6SG), which can efficiently crosslink to interacting proteins using UVA and UVB light. Crosslinking of 4SU or 6SG to interacting amino acids changes their base-pairing properties and results in characteristic mutations in cDNA libraries prepared for high-throughput sequencing, which can be computationally exploited to remove abundant background from non-crosslinked sequences and help pinpoint RNA binding protein binding sites at nucleotide resolution on a transcriptome-wide scale. Here we present a streamlined protocol for fluorescence-based PAR-CLIP (fPAR-CLIP) that eliminates the need to use radioactivity. It is based on direct ligation of a fluorescently labeled adapter to the 3'end of crosslinked RNA on immobilized ribonucleoproteins, followed by isolation of the adapter-ligated RNA and efficient conversion into cDNA without the previously needed size fractionation on denaturing polyacrylamide gels. These improvements cut the experimentation by half to 2 days and increases sensitivity by 10-100-fold.

Multisite Evaluation of Next-Generation Methods for Small RNA Quantification

Small RNAs (smRNAs) are important regulators of many biologic processes and are now most frequently characterized using Illumina sequencing. However, although standard RNA sequencing library preparation has become routine in most sequencing facilities, smRNA sequencing library preparation has historically been challenging because of high input requirements, laborious protocols involving gel purifications, inability to automate, and a lack of benchmarking standards. Additionally, studies have suggested that many of these methods are nonlinear and do not accurately reflect the amounts of smRNAs in vivo. Recently, a number of new kits have become available that permit lower input amounts and less laborious, gel-free protocol options. Several of these new kits claim to reduce RNA ligase-dependent sequence bias through novel adapter modifications and to lessen adapter-dimer contamination in the resulting libraries. With the increasing number of smRNA kits available, understanding the relative strengths of each method is crucial for appropriate experimental design. In this study, we systematically compared 9 commercially available smRNA library preparation kits as well as NanoString probe hybridization across multiple study sites. Although several of the new methodologies do reduce the amount of artificially over- and underrepresented microRNAs (miRNAs), we observed that none of the methods was able to remove all of the bias in the library preparation. Identical samples prepared with different methods show highly varied levels of different miRNAs. Even so, many methods excelled in ease of use, lower input requirement, fraction of usable reads, and reproducibility across sites. These differences may help users select the most appropriate methods for their specific question of interest.

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.

Commercially Available Blocking Oligonucleotides Effectively Suppress Unwanted Hemolysis-Related miRNAs in a Large Whole-Blood RNA Cohort

When sequencing small RNA libraries derived from whole blood, the most abundant microRNAs (miRs) detected are often miR-486-5p, miR-451a, and miR-92a-3p. These highly expressed erythropoietic miRs are released into the sample from red blood cell hemolysis. Next-generation sequencing of these unwanted miRs leads to a waste in sequencing cost and diminished detection of lowly expressed miRNAs, including many potential miRNA biomarkers. Previous work has developed a method to reduce targeted miRNAs using oligonucleotides that bind their target miRNA and prevent its ligation during library construction, although the extent to which oligonucleotides can be multiplexed and their effect on larger cohorts has not been thoroughly explored. We present a method for suppressing detection of three highly abundant heme miRs in a single multiplexed blocking oligonucleotide reaction. In a small paired-sample pilot (n = 8) and a large cohort of samples (n = 901), multiplexed oligos reduced detection of their target miRNAs by approximately 70%, allowing for an approximately 10-fold increase in reads mapping to nonheme miRs and increased detection of very lowly expressed miRs, with minimal off-target effects. By removing all three highly expressed erythropoietic miRNAs from next-generational sequencing libraries, this commercially available multiplexed blocking oligonucleotide method allows for greater detection of lowly expressed biomarkers, improving the efficacy, cost-efficiency, and sensitivity of biomarker studies and diagnostic tests.

Identification and Quantification of Small RNAs

RNA silencing plays a critical role in diverse biological processes in plants including growth, development, and responses to abiotic and biotic stresses. RNA silencing is guided by small non-coding RNAs (sRNAs) with the length of 21-24 nucleotides (nt) that are loaded into Argonaute (AGO) to repress expression of target loci and transcripts through transcriptional or posttranscriptional gene silencing mechanisms. Identification and quantitative characterization of sRNAs are crucial steps toward appreciation of their functions in biology. Here, we developed a step-by-step protocol to precisely illustrate the process of cloning of sRNA libraries and correspondingly computational analysis of the recovered sRNAs. This protocol can be used in all kinds of organisms, including Arabidopsis, and is compatible with various high-throughput sequence technologies such as Illumina Hiseq. Thus, we wish that this protocol represents an accurate way to identify and quantify sRNAs in vivo.

Human milk extracellular vesicle miRNA expression and associations with maternal characteristics in a population-based cohort from the Faroe Islands

Human milk plays a critical role in infant development and health, particularly in cognitive, immune, and cardiometabolic functions. Milk contains extracellular vesicles (EVs) that can transport biologically relevant cargo from mother to infant, including microRNAs (miRNAs). We aimed to characterize milk EV-miRNA profiles in a human population cohort, assess potential pathways and ontology, and investigate associations with maternal characteristics. We conducted the first study to describe the EV miRNA profile of human milk in 364 mothers from a population-based mother-infant cohort in the Faroe Islands using small RNA sequencing. We detected 1523 miRNAs with ≥ one read in 70% of samples. Using hierarchical clustering, we determined five EV-miRNA clusters, the top three consisting of 15, 27 and 67 miRNAs. Correlation coefficients indicated that the expression of many miRNAs within the top three clusters was highly correlated. Top-cluster human milk EV-miRNAs were involved in pathways enriched for the endocrine system, cellular community, neurodevelopment, and cancers. miRNA expression was associated with time to milk collection post-delivery, maternal body mass index, and maternal smoking, but not maternal parity. Future studies investigating determinants of human EV-miRNAs and associated health outcomes are needed to elucidate the role of human milk EV-miRNAs in health and disease.

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.

Charting Extracellular Transcriptomes in The Human Biofluid RNA Atlas

Extracellular RNAs present in biofluids have emerged as potential biomarkers for disease. Where most studies focus on blood-derived fluids, other biofluids may be more informative. We present an atlas of messenger, circular, and small RNA transcriptomes of a comprehensive collection of 20 human biofluids. By means of synthetic spike-in controls, we compare RNA content across biofluids, revealing a 10,000-fold difference in concentration. The circular RNA fraction is increased in most biofluids compared to tissues. Each biofluid transcriptome is enriched for RNA molecules derived from specific tissues and cell types. Our atlas enables an informed selection of the most relevant biofluid to monitor particular diseases. To verify the biomarker potential in these biofluids, four validation cohorts representing a broad spectrum of diseases were profiled, revealing numerous differential RNAs between case and control subjects. Spike-normalized data are publicly available in the R2 web portal for further exploration.

Applications of phage-derived RNA-based technologies in synthetic biology

As the most abundant biological entities with incredible diversity, bacteriophages (also known as phages) have been recognized as an important source of molecular machines for the development of genetic-engineering tools. At the same time, phages are crucial for establishing and improving basic theories of molecular biology. Studies on phages provide rich sources of essential elements for synthetic circuit design as well as powerful support for the improvement of directed evolution platforms. Therefore, phages play a vital role in the development of new technologies and central scientific concepts. After the RNA world hypothesis was proposed and developed, novel biological functions of RNA continue to be discovered. RNA and its related elements are widely used in many fields such as metabolic engineering and medical diagnosis, and their versatility led to a major role of RNA in synthetic biology. Further development of RNA-based technologies will advance synthetic biological tools as well as provide verification of the RNA world hypothesis. Most synthetic biology efforts are based on reconstructing existing biological systems, understanding fundamental biological processes, and developing new technologies. RNA-based technologies derived from phages will offer abundant sources for synthetic biological components. Moreover, phages as well as RNA have high impact on biological evolution, which is pivotal for understanding the origin of life, building artificial life-forms, and precisely reprogramming biological systems. This review discusses phage-derived RNA-based technologies terms of phage components, the phage lifecycle, and interactions between phages and bacteria. The significance of RNA-based technology derived from phages for synthetic biology and for understanding the earliest stages of biological evolution will be highlighted.

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.

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.

Classifying Lung Neuroendocrine Neoplasms through MicroRNA Sequence Data Mining

Lung neuroendocrine neoplasms (NENs) can be challenging to classify due to subtle histologic differences between pathological types. MicroRNAs (miRNAs) are small RNA molecules that are valuable markers in many neoplastic diseases. To evaluate miRNAs as classificatory markers for lung NENs, we generated comprehensive miRNA expression profiles from 14 typical carcinoid (TC), 15 atypical carcinoid (AC), 11 small cell lung carcinoma (SCLC), and 15 large cell neuroendocrine carcinoma (LCNEC) samples, through barcoded small RNA sequencing. Following sequence annotation and data preprocessing, we randomly assigned these profiles to discovery and validation sets. Through high expression analyses, we found that miR-21 and -375 are abundant in all lung NENs, and that miR-21/miR-375 expression ratios are significantly lower in carcinoids (TC and AC) than in neuroendocrine carcinomas (NECs; SCLC and LCNEC). Subsequently, we ranked and selected miRNAs for use in miRNA-based classification, to discriminate carcinoids from NECs. Using miR-18a and -155 expression, our classifier discriminated these groups in discovery and validation sets, with 93% and 100% accuracy. We also identified miR-17, -103, and -127, and miR-301a, -106b, and -25, as candidate markers for discriminating TC from AC, and SCLC from LCNEC, respectively. However, these promising findings require external validation due to sample size.

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.

cfRNAs as biomarkers in oncology - still experimental or applied tool for personalized medicine already?

Currently, the challenges of contemporary oncology are focused mainly on the development of personalized medicine and precise treatment, which could be achieved through the use of molecular biomarkers. One of the biological molecules with great potential are circulating free RNAs (cfRNAs) which are present in various types of body fluids, such as blood, serum, plasma, and saliva. Also, different types of cfRNA particles can be distinguished depending on their length and function: microRNA (miRNA), PIWI-interacting RNA (piRNA), tRNA-derived RNA fragments (tRFs), circular RNA (circRNA), long non-coding RNA (lncRNA), and messenger RNA (mRNA). Moreover, cfRNAs occur in various forms: as a free molecule alone, in membrane vesicles, such as exosomes, or in complexes with proteins and lipids. One of the modern approaches for monitoring patient's condition is a "liquid biopsy" that provides a non-invasive and easily available source of circulating RNAs. Both the presence of specific cfRNA types as well as their concentration are dependent on many factors including cancer type or even reaction to treatment. Despite the possibility of using circulating free RNAs as biomarkers, there is still a lack of validated diagnostic panels, defined protocols for sampling, storing as well as detection methods. In this work we examine different types of cfRNAs, evaluate them as possible biomarkers, and analyze methods of their detection. We believe that further research on cfRNA and defining diagnostic panels could lead to better and faster cancer identification and improve treatment monitoring.

Literature review and methodological considerations for understanding circulating risk biomarkers following trauma exposure

Exposure to traumatic events is common. While many individuals recover following trauma exposure, a substantial subset develop adverse posttraumatic neuropsychiatric sequelae (APNS) such as posttraumatic stress, major depression, and regional or widespread chronic musculoskeletal pain. APNS cause substantial burden to the individual and to society, causing functional impairment and physical disability, risk for suicide, lost workdays, and increased health care costs. Contemporary treatment is limited by an inability to identify individuals at high risk of APNS in the immediate aftermath of trauma, and an inability to identify optimal treatments for individual patients. Our purpose is to provide a comprehensive review describing candidate blood-based biomarkers that may help to identify those at high risk of APNS and/or guide individual intervention decision-making. Such blood-based biomarkers include circulating biological factors such as hormones, proteins, immune molecules, neuropeptides, neurotransmitters, mRNA, and noncoding RNA expression signatures, while we do not review genetic and epigenetic biomarkers due to other recent reviews of this topic. The current state of the literature on circulating risk biomarkers of APNS is summarized, and key considerations and challenges for their discovery and translation are discussed. We also describe the AURORA study, a specific example of current scientific efforts to identify such circulating risk biomarkers and the largest study to date focused on identifying risk and prognostic factors in the aftermath of trauma exposure.

Deep sequencing of non-enzymatic RNA primer extension

Life emerging in an RNA world is expected to propagate RNA as hereditary information, requiring some form of primitive replication without enzymes. Non-enzymatic template-directed RNA primer extension is a model of the copying step in this posited form of replication. The sequence space accessed by primer extension dictates potential pathways to self-replication and, eventually, ribozymes. Which sequences can be accessed? What is the fidelity of the reaction? Does the recently illuminated mechanism of primer extension affect the distribution of sequences that can be copied? How do sequence features respond to experimental conditions and prebiotically relevant contexts? To help answer these and related questions, we here introduce a deep-sequencing methodology for studying RNA primer extension. We have designed and vetted special RNA constructs for this purpose, honed a protocol for sample preparation and developed custom software that analyzes sequencing data. We apply this new methodology to proof-of-concept controls, and demonstrate that it works as expected and reports on key features of the sequences accessed by primer extension.

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.

A comparison between SOLiD 5500XLand Ion Torrent PGM-derived miRNA expression profiles in two breast cell lines

Next-generation sequencing (NGS) platforms allow the analysis of hundreds of millions of molecules in a single sequencing run, revolutionizing many research areas. NGS-based microRNA studies enable expression quantification in unprecedented scale without the limitations of closed-platforms. Yet, whereas a massive amount of data produced by these platforms is available, comparisons of quantification/discovery capabilities between platforms are still lacking. Here we compare two NGS-platforms: SOLiD and PGM, by evaluating their microRNA identification/quantification capabilities using two breast-derived cell-lines. A high expression correlation (R2 > 0.9) was achieved, encompassing 97% of the miRNAs, and the few discrepancies in miRNA counts were attributable to molecules that have very low expression. Quantification divergences indicative of artefactual representation were seen for 14 miRNAs (higher in SOLiD-reads) and another 10 miRNAs more abundant in PGM-data. An inspection of these revealed an increased and statistically significant count of uracyls and uracyl-stretches for PGM-enriched miRNAs, compared to SOLiD and to the miRBase. In parallel, adenines and adenine-stretches were enriched for SOLiDderived miRNA reads. We conclude that, whereas both platforms are overall consistent and can be used interchangeably for microRNA expression studies, particular sequence features appear to be indicative of specific platform bias, and their presence in microRNAs should be considered for database-analyses.

A whole-tissue RNA-seq toolkit for organism-wide studies of gene expression with PME-seq

The immune system operates at the scale of the whole organism in mammals. We currently lack experimental approaches to systematically track and study organism-wide molecular processes in mice. Here we describe an integrated toolkit for measuring gene expression in whole tissues, 3-prime mRNA extension sequencing, that is applicable to most mouse organs and any mouse model of interest. Further, the methods of RNA-seq described in this protocol are broadly applicable to other sample types beyond whole organs, such as tissue samples or isolated cell populations. We report procedures to collect, store and lyse a dozen organ types using conditions compatible with the extraction of high-quality RNA. In addition, we detail protocols to perform high-throughput and low-cost RNA extraction and sequencing, as well as downstream data analysis. The protocol takes 5 d to process 384 mouse organs from collecting tissues to obtaining raw sequencing data, with additional time required for data analysis and mining. The protocol is accessible to individuals with basic skills in (i) mouse perfusion and dissection for sample collection and (ii) computation using Unix and R for data analysis. Overall, the methods presented here fill a gap in our toolbox for studying organism-wide processes in immunology and physiology.

Depletion of erythropoietic miR-486-5p and miR-451a improves detectability of rare microRNAs in peripheral blood-derived small RNA sequencing libraries

Erythroid-specific miR-451a and miR-486-5p are two of the most dominant microRNAs (miRNAs) in human peripheral blood. In small RNA sequencing libraries, their overabundance reduces diversity as well as complexity and consequently causes negative effects such as missing detectability and inaccurate quantification of low abundant miRNAs. Here we present a simple, cost-effective and easy to implement hybridization-based method to deplete these two erythropoietic miRNAs from blood-derived RNA samples. By utilization of blocking oligonucleotides, this method provides a highly efficient and specific depletion of miR-486-5p and miR-451a, which leads to a considerable increase of measured expression as well as detectability of low abundant miRNA species. The blocking oligos are compatible with common 5′ ligation-dependent small RNA library preparation protocols, including commercially available kits, such as Illumina TruSeq and Perkin Elmer NEXTflex. Furthermore, the here described method and oligo design principle can be easily adapted to target many other miRNA molecules, depending on context and research question.