Reducing ligation bias of small RNAs in libraries for next generation sequencing

SnoRNAs in cardiovascular development, function, and disease

Small nucleolar RNAs (snoRNAs) are emerging as important regulators of cardiovascular (patho)biology. Several roles of snoRNAs have recently been identified in heart development and congenital heart diseases, as well as their dynamic regulation in hypertrophic and dilated cardiomyopathies, coronary heart disease (CHD), myocardial infarction (MI), cardiac fibrosis, and heart failure. Furthermore, reports of changes in vesicular snoRNA expression and altered levels of circulating snoRNAs in response to cardiac stress suggest that snoRNAs also function in cardiac signaling and intercellular communication. In this review, we summarize and discuss key findings and outline the clinical potential of snoRNAs considering current challenges and gaps in the field of cardiovascular diseases (CVDs).

Isoform-specific RNA structure determination using Nano-DMS-MaP

RNA structure determination is essential to understand how RNA carries out its diverse biological functions. In cells, RNA isoforms are readily expressed with partial variations within their sequences due, for example, to alternative splicing, heterogeneity in the transcription start site, RNA processing or differential termination/polyadenylation. Nanopore dimethyl sulfate mutational profiling (Nano-DMS-MaP) is a method for in situ isoform-specific RNA structure determination. Unlike similar methods that rely on short sequencing reads, Nano-DMS-MaP employs nanopore sequencing to resolve the structures of long and highly similar RNA molecules to reveal their previously hidden structural differences. This Protocol describes the development and applications of Nano-DMS-MaP and outlines the main considerations for designing and implementing a successful experiment: from bench to data analysis. In cell probing experiments can be carried out by an experienced molecular biologist in 3-4 d. Data analysis requires good knowledge of command line tools and Python scripts and requires a further 3-5 d.

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

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

Characterisation and engineering of a thermophilic RNA ligase from Palaeococcus pacificus

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

NAP-seq reveals multiple classes of structured noncoding RNAs with regulatory functions

Up to 80% of the human genome produces "dark matter" RNAs, most of which are noncapped RNAs (napRNAs) that frequently act as noncoding RNAs (ncRNAs) to modulate gene expression. Here, by developing a method, NAP-seq, to globally profile the full-length sequences of napRNAs with various terminal modifications at single-nucleotide resolution, we reveal diverse classes of structured ncRNAs. We discover stably expressed linear intron RNAs (sliRNAs), a class of snoRNA-intron RNAs (snotrons), a class of RNAs embedded in miRNA spacers (misRNAs) and thousands of previously uncharacterized structured napRNAs in humans and mice. These napRNAs undergo dynamic changes in response to various stimuli and differentiation stages. Importantly, we show that a structured napRNA regulates myoblast differentiation and a napRNA DINAP interacts with dyskerin pseudouridine synthase 1 (DKC1) to promote cell proliferation by maintaining DKC1 protein stability. Our approach establishes a paradigm for discovering various classes of ncRNAs with regulatory functions.

Uncovering transcriptomic biomarkers for enhanced diagnosis of methamphetamine use disorder: a comprehensive review

Introduction

Methamphetamine use disorder (MUD) is a chronic relapsing disorder characterized by compulsive Methamphetamine (MA) use despite its detrimental effects on physical, psychological, and social well-being. The development of MUD is a complex process that involves the interplay of genetic, epigenetic, and environmental factors. The treatment of MUD remains a significant challenge, with no FDA-approved pharmacotherapies currently available. Current diagnostic criteria for MUD rely primarily on self-reporting and behavioral assessments, which have inherent limitations owing to their subjective nature. This lack of objective biomarkers and unidimensional approaches may not fully capture the unique features and consequences of MA addiction.

Methods

We performed a literature search for this review using the Boolean search in the PubMed database.

Results

This review explores existing technologies for identifying transcriptomic biomarkers for MUD diagnosis. We examined non-invasive tissues and scrutinized transcriptomic biomarkers relevant to MUD. Additionally, we investigated transcriptomic biomarkers identified for diagnosing, predicting, and monitoring MUD in non-invasive tissues.

Discussion

Developing and validating non-invasive MUD biomarkers could address these limitations, foster more precise and reliable diagnostic approaches, and ultimately enhance the quality of care for individuals with MA addiction.

The challenges and potential in developing microRNA associated with regeneration as biomarkers to improve prognostication for liver failure syndromes and hepatocellular carcinoma

INTRODUCTION

Determining the need for liver transplantation remains critical in the management of hepatocellular carcinoma (HCC) and liver failure syndromes (including acute liver failure and decompensated cirrhosis states). Conventional prognostic models utilize biomarkers of liver and non-liver failure and have limitations in their application. Novel biomarkers which predict regeneration may fulfil this niche. microRNA are implicated in health and disease and are present in abundance in the circulation. Despite this, they have not translated into mainstream clinical biomarkers.

AREAS COVERED

We will discuss current challenges in the prognostication of patients with liver failure syndromes as well as for patients with HCC. We will discuss biomarkers implicated with liver regeneration. We then provide an overview of the challenges in developing microRNA into clinically tractable biomarkers. Finally, we will provide a scoping review of microRNA which may have potential as prognostic biomarkers in liver failure syndromes and HCC.

EXPERT OPINION

Novel biomarkers are needed to improve prognostic models in liver failure syndromes and HCC. Biomarkers associated with liver regeneration are currently lacking and may fulfil this niche. microRNA have the potential to be developed into clinically tractable biomarkers but a consensus on standardizing methodology and reporting is required prior to large-scale studies.

Identification of ascorbate- and salicylate-responsive miRNAs and verification of the spectral control of miR395 in Arabidopsis

We assumed that miRNAs might regulate the physiological and biochemical processes in plants through their effects on the redox system and phytohormones. To check this hypothesis, the transcriptome profile of wild-type Arabidopsis and lines with decreased ascorbate (Asc), glutathione (GSH), or salicylate (Sal) levels were compared. GSH deficiency did not influence the miRNA expression, whereas lower levels of Asc and Sal reduced the accumulation of 9 and 44 miRNAs, respectively, but only four miRNAs were upregulated. Bioinformatics analysis revealed that their over-represented target genes are associated with the synthesis of nitrogen-containing and aromatic compounds, nucleic acids, and sulphate assimilation. Among them, the sulphate reduction-related miR395 - ATP-sulfurylase couple was selected to check the assumed modulating role of the light spectrum. A greater induction of the Asc- and Sal-responsive miR395 was observed under sulphur starvation in far-red light compared to white and blue light in wild-type and GSH-deficient Arabidopsis lines. Sal deficiency inhibited the induction of miR395 by sulphur starvation in blue light, whereas Asc deficiency greatly reduced it independently of the spectrum. Interestingly, sulphur starvation decreased only the level of ATP sulfurylase 4 among the miR395 target genes in far-red light. The expression level of ATP sulfurylase 3 was higher in far-red light than in blue light in wild-type and Asc-deficient lines. The results indicate the coordinated control of miRNAs by the redox and hormonal system since 11 miRNAs were affected by both Asc and Sal deficiency. This process can be modulated by light spectrum, as shown for miR395.

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.

Intrinsically disordered proteins SAID1/2 condensate on SERRATE for dual inhibition of miRNA biogenesis in Arabidopsis

Intrinsically disordered proteins (IDPs) SAID1/2 are hypothetic dentin sialophosphoprotein-like proteins, but their true functions are unknown. Here, we identified SAID1/2 as negative regulators of SERRATE (SE), a core factor in miRNA biogenesis complex (microprocessor). Loss-of-function double mutants of said1; said2 caused pleiotropic developmental defects and thousands of differentially expressed genes that partially overlapped with those in se. said1; said2 also displayed increased assembly of microprocessor and elevated accumulation of microRNAs (miRNAs). Mechanistically, SAID1/2 promote pre-mRNA processing 4 kinase A-mediated phosphorylation of SE, causing its degradation in vivo. Unexpectedly, SAID1/2 have strong binding affinity to hairpin-structured pri-miRNAs and can sequester them from SE. Moreover, SAID1/2 directly inhibit pri-miRNA processing by microprocessor in vitro. Whereas SAID1/2 did not impact SE subcellular compartmentation, the proteins themselves exhibited liquid-liquid phase condensation that is nucleated on SE. Thus, we propose that SAID1/2 reduce miRNA production through hijacking pri-miRNAs to prevent microprocessor activity while promoting SE phosphorylation and its destabilization in Arabidopsis.

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.

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.

The Sum of Two Halves May Be Different from the Whole-Effects of Splitting Sequencing Samples Across Lanes

The advances in high-throughput sequencing (HTS) have enabled the characterisation of biological processes at an unprecedented level of detail; most hypotheses in molecular biology rely on analyses of HTS data. However, achieving increased robustness and reproducibility of results remains a main challenge. Although variability in results may be introduced at various stages, e.g., alignment, summarisation or detection of differential expression, one source of variability was systematically omitted: the sequencing design, which propagates through analyses and may introduce an additional layer of technical variation. We illustrate qualitative and quantitative differences arising from splitting samples across lanes on bulk and single-cell sequencing. For bulk mRNAseq data, we focus on differential expression and enrichment analyses; for bulk ChIPseq data, we investigate the effect on peak calling and the peaks' properties. At the single-cell level, we concentrate on identifying cell subpopulations. We rely on markers used for assigning cell identities; both smartSeq and 10× data are presented. The observed reduction in the number of unique sequenced fragments limits the level of detail on which the different prediction approaches depend. Furthermore, the sequencing stochasticity adds in a weighting bias corroborated with variable sequencing depths and (yet unexplained) sequencing bias. Subsequently, we observe an overall reduction in sequencing complexity and a distortion in the biological signal across technologies, experimental contexts, organisms and tissues.

Exosome RNA Sequencing as a Tool in the Search for Cancer Biomarkers

Numerous noninvasive methods are currently being used to determine biomarkers for diseases such as cancer. However, these methods are not always precise and reliable. Thus, there is an unmet need for better diagnostic and prognostic biomarkers that will be used to diagnose cancer in early, more treatable stages of the disease. Exosomes are extracellular vesicles of endocytic origin released by the majority of cells. Exosomes contain and transport nucleic acids, proteins, growth factors, and cytokines from their parent cells to surrounding or even distant cells via circulation in biofluids. Exosomes have attracted the interest of researchers, as recent data indicate that exosome content may be indicative of disease stages and may contribute to disease progression via exosome-mediated extracellular communication. Therefore, the contents of these vesicles are being investigated as possible biomarkers for disease diagnosis and prognosis. The functions of exosomes and their contents in disease development are becoming clearer as isolation and analytical methods, such as RNA sequencing, advance. In this review, we discuss current advances and challenges in exosomal content analyses with emphasis on information that can be generated using RNA sequencing. We also discuss how the RNA sequencing of exosomes may be used to discover novel biomarkers for the detection of different stages for various cancers using specific microRNAs that were found to be differentially expressed between healthy controls and cancer-diagnosed subjects.

Single-base resolution mapping of 2′-O-methylation sites by an exoribonuclease-enriched chemical method

2'-O-methylation (Nm) is one of the most abundant RNA epigenetic modifications and plays a vital role in the post-transcriptional regulation of gene expression. Current Nm mapping approaches are normally limited to highly abundant RNAs and have significant technical hurdles in mRNAs or relatively rare non-coding RNAs (ncRNAs). Here, we developed a new method for enriching Nm sites by using RNA exoribonuclease and periodate oxidation reactivity to eliminate 2'-hydroxylated (2'-OH) nucleosides, coupled with sequencing (Nm-REP-seq). We revealed several novel classes of Nm-containing ncRNAs as well as mRNAs in humans, mice, and drosophila. We found that some novel Nm sites are present at fixed positions in different tRNAs and are potential substrates of fibrillarin (FBL) methyltransferase mediated by snoRNAs. Importantly, we discovered, for the first time, that Nm located at the 3'-end of various types of ncRNAs and fragments derived from them. Our approach precisely redefines the genome-wide distribution of Nm and provides new technologies for functional studies of Nm-mediated gene regulation.

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.

Mechanistic insights into non-coding Y RNA processing

Y RNAs (84–112 nt) are non-coding RNAs transcribed by RNA polymerase III and are characterized by a distinctive secondary structure. Human Y RNAs interact with the autoimmune proteins SSB and RO60 that together form a ribonucleoprotein (RNP) complex termed RoRNP and Y RNAs also perform regulatory roles in DNA and RNA replication and stability, which has major implications for diseases including cancer. During cellular stress and apoptosis, Y RNAs are cleaved into 3’ and 5’ end fragments termed Y RNA-derived small RNAs (ysRNAs). Although some ysRNA functions in stress, apoptosis and cancer have been reported, their fundamental biogenesis has not been described. Here we report that 3’ end RNY5 cleavage is structure dependent. In high throughput mutagenesis experiments, cleavage occurred between the 2^nd and 3^rd nt above a double stranded stem comprising high GC content. We demonstrate that an internal loop above stem S3 is critical for producing 3’ end ysRNAs (31 nt) with mutants resulting in longer or no ysRNAs. We show a UGGGU sequence motif at position 22 of RNY5 is critical for producing 5’ end ysRNAs (22–25 nt). We show that intact RO60 is critical for ysRNA biogenesis. We conclude that ribonuclease L (RNASEL) contributes to Y RNA cleavage in mouse embryonic fibroblasts but is not the only endoribonuclease important in human cells.

Technological Approaches in the Analysis of Extracellular Vesicle Nucleotide Sequences

Together with metabolites, proteins, and lipid components, the EV cargo consists of DNA and RNA nucleotide sequence species, which are part of the intracellular communication network regulating specific cellular processes and provoking distinct target cell responses. The extracellular vesicle (EV) nucleotide sequence cargo molecules are often investigated in association with a particular pathology and may provide an insight into the physiological and pathological processes in hard-to-access organs and tissues. The diversity and biological function of EV nucleotide sequences are distinct regarding EV subgroups and differ in tissue- and cell-released EVs. EV DNA is present mainly in apoptotic bodies, while there are different species of EV RNAs in all subgroups of EVs. A limited sample volume of unique human liquid biopsy provides a small amount of EVs with limited isolated DNA and RNA, which can be a challenging factor for EV nucleotide sequence analysis, while the additional difficulty is technical variability of molecular nucleotide detection. Every EV study is challenged with its first step of the EV isolation procedure, which determines the EV's purity, yield, and diameter range and has an impact on the EV's downstream analysis with a significant impact on the final result. The gold standard EV isolation procedure with ultracentrifugation provides a low output and not highly pure isolated EVs, while modern techniques increase EV's yield and purity. Different EV DNA and RNA detection techniques include the PCR procedure for nucleotide sequence replication of the molecules of interest, which can undergo a small-input EV DNA or RNA material. The nucleotide sequence detection approaches with their advantages and disadvantages should be considered to appropriately address the study problem and to extract specific EV nucleotide sequence information with the detection using qPCR or next-generation sequencing. Advanced next-generation sequencing techniques allow the detection of total EV genomic or transcriptomic data even at the single-molecule resolution and thus, offering a sensitive and accurate EV DNA or RNA biomarker detection. Additionally, with the processes where the EV genomic or transcriptomic data profiles are compared to identify characteristic EV differences in specific conditions, novel biomarkers could be discovered. Therefore, a suitable differential expression analysis is crucial to define the EV DNA or RNA differences between conditions under investigation. Further bioinformatics analysis can predict molecular cell targets and identify targeted and affected cellular pathways. The prediction target tools with functional studies are essential to help specify the role of the investigated EV-targeted nucleotide sequences in health and disease and support further development of EV-related therapeutics. This review will discuss the biological diversity of human liquid biopsy-obtained EV nucleotide sequences DNA and RNA species reported as potential biomarkers in health and disease and methodological principles of their detection, from human liquid biopsy EV isolation, EV nucleotide sequence extraction, techniques for their detection, and their cell target prediction.

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.

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.

miR-7b-3p Exerts a Dual Role After Spinal Cord Injury, by Supporting Plasticity and Neuroprotection at Cortical Level

Spinal cord injury (SCI) affects 6 million people worldwide with no available treatment. Despite research advances, the inherent poor regeneration potential of the central nervous system remains a major hurdle. Small RNAs (sRNAs) 19-33 nucleotides in length are a set of non-coding RNA molecules that regulate gene expression and have emerged as key players in regulating cellular events occurring after SCI. Here we profiled a class of sRNA known as microRNAs (miRNAs) following SCI in the cortex where the cell bodies of corticospinal motor neurons are located. We identified miR-7b-3p as a candidate target given its significant upregulation after SCI in vivo and we screened by miRWalk PTM the genes predicted to be targets of miR-7b-3p (among which we identified Wipf2, a gene regulating neurite extension). Moreover, 16 genes, involved in neural regeneration and potential miR-7b-3p targets, were found to be downregulated in the cortex following SCI. We also analysed miR-7b-3p function during cortical neuron development in vitro: we observed that the overexpression of miR-7b-3p was important (1) to maintain neurons in a more immature and, likely, plastic neuronal developmental phase and (2) to contrast the apoptotic pathway; however, in normal conditions it did not affect the Wipf2 expression. On the contrary, the overexpression of miR-7b-3p upon in vitro oxidative stress condition (mimicking the SCI environment) significantly reduced the expression level of Wipf2, as observed in vivo, confirming it as a direct miR-7b-3p target. Overall, these data suggest a dual role of miR-7b-3p: (i) the induction of a more plastic neuronal condition/phase, possibly at the expense of the axon growth, (ii) the neuroprotective role exerted through the inhibition of the apoptotic cascade. Increasing the miR-7b-3p levels in case of SCI could reactivate in adult neurons silenced developmental programmes, supporting at the same time the survival of the axotomised neurons.

Endogenous activated small interfering RNAs in virus-infected Brassicaceae crops show a common host gene-silencing pattern affecting photosynthesis and stress response

Viral infections are accompanied by a massive production of small interfering RNAs (siRNAs) of plant origin, such as virus-activated (va)siRNAs, which drive the widespread silencing of host gene expression, and whose effects in plant pathogen interactions remain unknown. By combining phenotyping and molecular analyses, we characterized vasiRNAs that are associated with typical mosaic symptoms of cauliflower mosaic virus infection in two crops, turnip (Brassica rapa) and oilseed rape (Brassica napus), and the reference plant Arabidopsis thaliana. We identified 15 loci in the three infected plant species, whose transcripts originate vasiRNAs. These loci appear to be generally affected by virus infections in Brassicaceae and encode factors that are centrally involved in photosynthesis and stress response, such as Rubisco activase (RCA), senescence-associated protein, heat shock protein HSP70, light harvesting complex, and membrane-related protein CP5. During infection, the expression of these factors is significantly downregulated, suggesting that their silencing is a central component of the plant's response to virus infections. Further findings indicate an important role for 22 nt long vasiRNAs in the plant's endogenous RNA silencing response. Our study considerably enhances knowledge about the new class of vasiRNAs that are triggered in virus-infected plants and will help to advance strategies for the engineering of gene clusters involved in the development of crop diseases.

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.

microRNA-seq of cartilage reveals an overabundance of miR-140-3p which contains functional isomiRs

miR-140 is selectively expressed in cartilage. Deletion of the entire Mir140 locus in mice results in growth retardation and early-onset osteoarthritis-like pathology; however, the relative contribution of miR-140-5p or miR-140-3p to the phenotype remains to be determined. An unbiased small RNA sequencing approach identified miR-140-3p as significantly more abundant (>10-fold) than miR-140-5p in human cartilage. Analysis of these data identified multiple miR-140-3p isomiRs differing from the miRBase annotation at both the 5' and 3' end, with >99% having one of two seed sequences (5' bases 2-8). Canonical (miR-140-3p.2) and shifted (miR-140-3p.1) seed isomiRs were overexpressed in chondrocytes and transcriptomics performed to identify targets. miR-140-3p.1 and miR-140-3p.2 significantly down-regulated 694 and 238 genes, respectively, of which only 162 genes were commonly down-regulated. IsomiR targets were validated using 3'UTR luciferase assays. miR-140-3p.1 targets were enriched within up-regulated genes in rib chondrocytes of Mir140-null mice and within down-regulated genes during human chondrogenesis. Finally, through imputing the expression of miR-140 from the expression of the host gene WWP2 in 124 previously published data sets, an inverse correlation with miR-140-3p.1 predicted targets was identified. Together these data suggest the novel seed containing isomiR miR-140-3p.1 is more functional than original consensus miR-140-3p seed containing isomiR.

mirnaQC: a webserver for comparative quality control of miRNA-seq data

Although miRNA-seq is extensively used in many different fields, its quality control is frequently restricted to a PhredScore-based filter. Other important quality related aspects like microRNA yield, the fraction of putative degradation products (such as rRNA fragments) or the percentage of adapter-dimers are hard to assess using absolute thresholds. Here we present mirnaQC, a webserver that relies on 34 quality parameters to assist in miRNA-seq quality control. To improve their interpretability, quality attributes are ranked using a reference distribution obtained from over 36 000 publicly available miRNA-seq datasets. Accepted input formats include FASTQ and SRA accessions. The results page contains several sections that deal with putative technical artefacts related to library preparation, sequencing, contamination or yield. Different visualisations, including PCA and heatmaps, are available to help users identify underlying issues. Finally, we show the usefulness of this approach by analysing two publicly available datasets and discussing the different quality issues that can be detected using mirnaQC.

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.

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.

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.

PAREameters: a tool for computational inference of plant miRNA-mRNA targeting rules using small RNA and degradome sequencing data

MicroRNAs (miRNAs) are short, non-coding RNAs that modulate the translation-rate of messenger RNAs (mRNAs) by directing the RNA-induced silencing complex to sequence-specific targets. In plants, this typically results in cleavage and subsequent degradation of the mRNA. Degradome sequencing is a high-throughput technique developed to capture cleaved mRNA fragments and thus can be used to support miRNA target prediction. The current criteria used for miRNA target prediction were inferred on a limited number of experimentally validated A. thaliana interactions and were adapted to fit these specific interactions; thus, these fixed criteria may not be optimal across all datasets (organisms, tissues or treatments). We present a new tool, PAREameters, for inferring targeting criteria from small RNA and degradome sequencing datasets. We evaluate its performance using a more extensive set of experimentally validated interactions in multiple A. thaliana datasets. We also perform comprehensive analyses to highlight and quantify the differences between subsets of miRNA-mRNA interactions in model and non-model organisms. Our results show increased sensitivity in A. thaliana when using the PAREameters inferred criteria and that using data-driven criteria enables the identification of additional interactions that further our understanding of the RNA silencing pathway in both model and non-model organisms.

Artificially induced phased siRNAs promote virus resistance in transgenic plants

We previously developed transgenic tobacco plants that were resistant to two geminiviruses. We generated resistance using RNAi constructs that produced trans-acting siRNA (tasiRNA) like secondary siRNAs known as phased siRNA (phasiRNA) that targeted several regions of Tomato Leaf Curl New Delhi Virus (ToLCNDV) and Tomato Leaf Curl Gujarat Virus (ToLCGV) transcripts encoding the RNA silencing suppressor proteins AC2 and AC4. Here, we performed degradome analysis to determine the precise cleavage sites of RNA-RNA interaction between phasiRNA and viral transcripts. We then applied our RNAi technology in tomato, which is the natural host for ToLCNDV and ToLCGV. The relative ease of developing and using phasiRNA constructs represents a significant technical advance in imparting virus resistance in crops and/or important model systems.

Accurate Adapter Information Is Crucial for Reproducibility and Reusability in Small RNA Seq Studies

A necessary pre-processing data analysis step is the removal of adapter sequences from the raw reads. While most adapter trimming tools require adapter sequence as an essential input, adapter information is often incomplete or missing. This can impact quantification of features, reproducibility of the study and might even lead to erroneous conclusions. Here, we provide examples to highlight the importance of specifying the adapter sequence by demonstrating the effect of using similar but different adapter sequences and identify additional potential sources of errors in the adapter trimming step. Finally, we propose solutions by which users can ensure their small RNA-seq data is fully annotated with adapter information.

Mapping RNA-chromatin interactions by sequencing with iMARGI

RNA-chromatin interactions represent an important aspect of the transcriptional regulation of genes and transposable elements. However, analyses of chromatin-associated RNAs (caRNAs) are often limited to one caRNA at a time. Here, we describe the iMARGI (in situ mapping of RNA-genome interactome) technique, which is used to discover caRNAs and reveal their respective genomic interaction loci. iMARGI starts with in situ crosslinking and genome fragmentation, followed by converting each proximal RNA-DNA pair into an RNA-linker-DNA chimeric sequence. These chimeric sequences are subsequently converted into a sequencing library suitable for paired-end sequencing. A standardized bioinformatic software package, iMARGI-Docker, is provided to decode the paired-end sequencing data into caRNA-DNA interactions. Compared to its predecessor MARGI (mapping RNA-genome interactions), the number of input cells for iMARGI is 3-5 million (a 100-fold reduction), experimental time is reduced, and clear checkpoints have been established. It takes a few hours a day and a total of 8 d to complete the construction of an iMARGI sequencing library and 1 d to carry out data processing with iMARGI-Docker.

Bias-minimized quantification of microRNA reveals widespread alternative processing and 3' end modification

MicroRNAs (miRNAs) modulate diverse biological and pathological processes via post-transcriptional gene silencing. High-throughput small RNA sequencing (sRNA-seq) has been widely adopted to investigate the functions and regulatory mechanisms of miRNAs. However, accurate quantification of miRNAs has been limited owing to the severe ligation bias in conventional sRNA-seq methods. Here, we quantify miRNAs and their variants (known as isomiRs) by an improved sRNA-seq protocol, termed AQ-seq (accurate quantification by sequencing), that utilizes adapters with terminal degenerate sequences and a high concentration of polyethylene glycol (PEG), which minimize the ligation bias during library preparation. Measurement using AQ-seq allows us to correct the previously misannotated 5' end usage and strand preference in public databases. Importantly, the analysis of 5' terminal heterogeneity reveals widespread alternative processing events which have been underestimated. We also identify highly uridylated miRNAs originating from the 3p strands, indicating regulations mediated by terminal uridylyl transferases at the pre-miRNA stage. Taken together, our study reveals the complexity of the miRNA isoform landscape, allowing us to refine miRNA annotation and to advance our understanding of miRNA regulation. Furthermore, AQ-seq can be adopted to improve other ligation-based sequencing methods including crosslinking-immunoprecipitation-sequencing (CLIP-seq) and ribosome profiling (Ribo-seq).

Systematic assessment of commercially available low-input miRNA library preparation kits

High-throughput sequencing is increasingly favoured to assay the presence and abundance of microRNAs (miRNAs) in biological samples, even from low RNA amounts, and a number of commercial vendors now offer kits that allow miRNA sequencing from sub-nanogram (ng) inputs. Although biases introduced during library preparation have been documented, the relative performance of current reagent kits has not been investigated in detail. Here, six commercial kits capable of handling <100ng total RNA input were used for library preparation, performed by kit manufactures, on synthetic miRNAs of known quantities and human total RNA samples. We compared the performance of miRNA detection sensitivity, reliability, titration response and the ability to detect differentially expressed miRNAs. In addition, we assessed the use of unique molecular identifiers (UMI) sequence tags in one kit. We observed differences in detection sensitivity and ability to identify differentially expressed miRNAs between the kits, but none were able to detect the full repertoire of synthetic miRNAs. The reliability within the replicates of all kits was good, while larger differences were observed between the kits, although none could accurately quantify the relative levels of the majority of miRNAs. UMI tags, at least within the input ranges tested, offered little advantage to improve data utility. In conclusion, biases in miRNA abundance are heavily influenced by the kit used for library preparation, suggesting that comparisons of datasets prepared by different procedures should be made with caution. This article is intended to assist researchers select the most appropriate kit for their experimental conditions.

CRISPR/Cas9-targeted removal of unwanted sequences from small-RNA sequencing libraries

In small RNA (smRNA) sequencing studies, highly abundant molecules such as adapter dimer products and tissue-specific microRNAs (miRNAs) inhibit accurate quantification of lowly expressed species. We previously developed a method to selectively deplete highly abundant miRNAs. However, this method does not deplete adapter dimer ligation products that, unless removed by gel-separation, comprise most of the library. Here, we have adapted and modified recently described methods for CRISPR/Cas9-based Depletion of Abundant Species by Hybridization ('DASH') to smRNA-seq, which we have termed miRNA and Adapter Dimer-DASH (MAD-DASH). In MAD-DASH, Cas9 is complexed with single guide RNAs (sgRNAs) targeting adapter dimer ligation products, alongside highly expressed tissue-specific smRNAs, for cleavage in vitro. This process dramatically reduces adapter dimer and targeted smRNA sequences, can be multiplexed, shows minimal off-target effects, improves the quantification of lowly expressed miRNAs from human plasma and tissue derived RNA, and obviates the need for gel-separation, greatly increasing sample throughput. Additionally, the method is fully customizable to other smRNA-seq preparation methods. Like depletion of ribosomal RNA for mRNA-seq and mitochondrial DNA for ATAC-seq, our method allows for greater proportional read-depth of non-targeted sequences.

Distinct roles of Argonaute in the green alga Chlamydomonas reveal evolutionary conserved mode of miRNA-mediated gene expression

The unicellular green alga Chlamydomonas reinhardtii is evolutionarily divergent from higher plants, but has a fully functional silencing machinery including microRNA (miRNA)-mediated translation repression and mRNA turnover. However, distinct from the metazoan machinery, repression of gene expression is primarily associated with target sites within coding sequences instead of 3′UTRs. This feature indicates that the miRNA-Argonaute (AGO) machinery is ancient and the primary function is for post transcriptional gene repression and intermediate between the mechanisms in the rest of the plant and animal kingdoms. Here, we characterize AGO2 and 3 in Chlamydomonas, and show that cytoplasmically enriched Cr-AGO3 is responsible for endogenous miRNA-mediated gene repression. Under steady state, mid-log phase conditions, Cr-AGO3 binds predominantly miR-C89, which we previously identified as the predominant miRNA with effects on both translation repression and mRNA turnover. In contrast, the paralogue Cr-AGO2 is nuclear enriched and exclusively binds to 21-nt siRNAs. Further analysis of the highly similar Cr-AGO2 and Cr-AGO 3 sequences (90% amino acid identity) revealed a glycine-arginine rich N-terminal extension of ~100 amino acids that, given previous work on unicellular protists, may associate AGO with the translation machinery. Phylogenetic analysis revealed that this glycine-arginine rich N-terminal extension is present outside the animal kingdom and is highly conserved, consistent with our previous proposal that miRNA-mediated CDS-targeting operates in this green alga.

Comprehensive assessment of multiple biases in small RNA sequencing reveals significant differences in the performance of widely used methods

BACKGROUND

RNA sequencing offers advantages over other quantification methods for microRNA (miRNA), yet numerous biases make reliable quantification challenging. Previous evaluations of these biases have focused on adapter ligation bias with limited evaluation of reverse transcription bias or amplification bias. Furthermore, evaluations of the quantification of isomiRs (miRNA isoforms) or the influence of starting amount on performance have been very limited. No study had yet evaluated the quantification of isomiRs of altered length or compared the consistency of results derived from multiple moderate starting inputs. We therefore evaluated quantifications of miRNA and isomiRs using four library preparation kits, with various starting amounts, as well as quantifications following removal of duplicate reads using unique molecular identifiers (UMIs) to mitigate reverse transcription and amplification biases.

RESULTS

All methods resulted in false isomiR detection; however, the adapter-free method tested was especially prone to false isomiR detection. We demonstrate that using UMIs improves accuracy and we provide a guide for input amounts to improve consistency.

CONCLUSIONS

Our data show differences and limitations of current methods, thus raising concerns about the validity of quantification of miRNA and isomiRs across studies. We advocate for the use of UMIs to improve accuracy and reliability of miRNA quantifications.

A comparison of RNA extraction and sequencing protocols for detection of small RNAs in plasma

BACKGROUND

Circulating microRNAs (miRNAs) are attractive non-invasive biomarkers for a variety of conditions due to their stability and altered pathophysiological expression levels. Reliable detection of global expression profiles is required to maximise miRNA biomarker discovery. Although developments in small RNA-Seq technology have improved detection of plasma-based miRNAs, the low RNA content and sequencing bias introduced during library preparation remain challenging. In this study we compare commercially available RNA extraction methods using MagnaZol (Bioo Scientific) or miRNeasy (QIAGEN) and three library preparation methods - CleanTag (TriLink), NEXTflex (Bioo Scientific) and QIAseq (QIAGEN) - which aim to address one or both of these issues.

RESULTS

Different RNA extractions and library preparation protocols result in differential detection of miRNAs. A greater proportion of reads mapped to miRNAs in libraries prepared with MagnaZol RNA than with miRNeasy RNA. Libraries prepared using QIAseq demonstrated the greatest miRNA diversity with many more very low abundance miRNAs detected (~ 2-3 fold more with < 10 reads), whilst CleanTag detected the fewest individual miRNAs and considerably over-represented miR-486-5p. Libraries prepared with QIAseq had the strongest correlation with RT-qPCR quantification. Analysis of unique molecular indices (UMIs) incorporated in the QIAseq protocol indicate that little PCR bias is introduced during small RNA library preparation.

CONCLUSIONS

Small RNAs were consistently detected using all RNA extraction and library preparation protocols tested, but with some miRNAs at significantly different levels. Choice of the most suitable protocol should be informed by the relative importance of minimising the total sequencing required, detection of rare miRNAs or absolute quantification.

Molecular insights into an ancient form of Paget's disease of bone

We identify an ancient and atypical form of Paget’s disease of bone (PDB) in a collection of medieval skeletons exhibiting unusually extensive pathological changes, high disease prevalence, and low age-at-death estimations. Proteomic analysis of ancient bone-preserved proteins combined with analysis of small RNAs supports a retrospective diagnosis of PDB. Remains affected by other skeletal disorders may therefore hold a chemical memory amenable to similar molecular interrogation. Abnormalities in a contemporary PDB-linked protein detected in ancient tooth samples indicate that dentition may represent an unexplored storehouse for the study of skeletal disorders. Our work provides insights into the natural history of PDB and prompts a similar revaluation of other archaeological collections.

Paget’s disease of bone (PDB) is a chronic skeletal disorder that can affect one or several bones in individuals older than 55 y of age. PDB-like changes have been reported in archaeological remains as old as Roman, although accurate diagnosis and natural history of the disease is lacking. Six skeletons from a collection of 130 excavated at Norton Priory in the North West of England, which dates to medieval times, show atypical and extensive pathological changes resembling contemporary PDB affecting as many as 75% of individual skeletons. Disease prevalence in the remaining collection is high, at least 16% of adults, with age at death estimations as low as 35 y. Despite these atypical features, paleoproteomic analysis identified sequestosome 1 (SQSTM1) or p62, a protein central to the pathological milieu of PDB, as one of the few noncollagenous human sequences preserved in skeletal samples. Targeted proteomic analysis detected >60% of the ancient p62 primary sequence, with Western blotting indicating p62 abnormalities, including in dentition. Direct sequencing of ancient DNA excluded contemporary PDB-associated SQSTM1 mutations. Our observations indicate that the ancient p62 protein is likely modified within its C-terminal ubiquitin-associated domain. Ancient miRNAs were remarkably preserved in an osteosarcoma from a skeleton with extensive disease, with miR-16 expression consistent with that reported in contemporary PDB-associated bone tumors. Our work displays the use of proteomics to inform diagnosis of ancient diseases such as atypical PDB, which has unusual features presumably potentiated by yet-unidentified environmental or genetic factors.

Evaluating and Correcting Inherent Bias of microRNA Expression in Illumina Sequencing Analysis

microRNA (miRNA) expression profiles based on the highly powerful Illumina sequencing technology rely on the construction of cDNA libraries in which adaptor ligation is known to deeply favor some miRNAs over others. This introduces erroneous measurements of the miRNA abundances and relative miRNA quantities in biological samples. Here, by using the commercial miRXplore Universal Reference that contains an equimolar mixture of 963 animal miRNAs and TruSeq or bulged adaptors, we describe a method for correcting ligation biases in expression profiles obtained with standard protocols of cDNA library construction and provide data for quantifying the true miRNA abundances in biological samples. Ligation biases were evaluated at three ratios of miRNA to 3′-adaptor and four numbers of polymerase chain reaction amplification cycles by calculating efficiency captures/correcting factors for each miRNA. We show that ligation biases lead to over- or under-expression covering a 10⁵ amplitude range. We also show that, at each miRNA:3′-adaptor ratio, coefficients of variation (CVs) of efficiency captures calculated over the four number of amplification cycles using sliding windows of 10 values ranged from 0.1 for the miRNAs of high expression to 0.6 for the miRNAs of low expression. Efficiency captures of miRNAs of high and low expression in profiles are therefore differently impacted by the number of amplification cycles. Importantly, we observed that at a given number of amplification cycles, CVs of efficiency captures calculated over the three miRNA:3′-adaptor ratios displayed a steady value of 0.3 +/− 0.05 STD for miRNAs of high and low expression. This allows, at a given number of amplification cycles, accurate comparison of miRNA expression between biological samples over a substantial expression range. Finally we provide tables of correcting factors that allow to measure the abundances of 963 miRNAs in biological samples from TruSeq-based expression profiles and, an example of their use by characterizing miRNAs of the let-7, miR-26, miR-29, and miR-30 families as the more abundant miRNAs of the rat adult cerebellum.

Functional assays for transcription mechanisms in high-throughput

Dramatic increases in the scale of programmed synthesis of nucleic acid libraries coupled with deep sequencing have powered advances in understanding nucleic acid and protein biology. Biological systems centering on nucleic acids or encoded proteins greatly benefit from such high-throughput studies, given that large DNA variant pools can be synthesized and DNA, or RNA products of transcription, can be easily analyzed by deep sequencing. Here we review the scope of various high-throughput functional assays for studies of nucleic acids and proteins in general, followed by discussion of how these types of study have yielded insights into the RNA Polymerase II (Pol II) active site as an example. We discuss methodological considerations in the design and execution of these experiments that should be valuable to studies in any system.

Small RNA-seq: The RNA 5'-end adapter ligation problem and how to circumvent it

The preparation of small RNA cDNA sequencing libraries depends on the unbiased ligation of adapters to the RNA ends. Small RNA with 5' recessed ends are poor substrates for enzymatic adapter ligation, but this 5' adapter ligation problem can go undetected if the library preparation steps are not monitored. Here we illustrate the severity of the 5' RNA end ligation problem using several pre-miRNA-like hairpins that allow us to expand the definition of the problem to include 5' ends close to a hairpin stem, whether recessed or in a short extension. The ribosome profiling method can avoid a difficult 5' adapter ligation, but the enzyme typically used to circularize the cDNA has been reported to be biased, calling into question the benefit of this workaround. Using the TS2126 RNA ligase 1 (a.k.a. CircLigase) as the circularizing enzyme, we devised a bias test for the circularization of first strand cDNA. All possible dinucleotides were circle-ligated with similar efficiency. To re-linearize the first strand cDNA in the ribosome profiling approach, we introduce an improved method wherein a single ribonucleotide is placed between the sequencing primer binding sites in the reverse transcriptase primer, which later serves as the point of re-linearization by RNase A. We incorporate this step into the ribosomal profiling method and describe a complete improved library preparation method, Coligo-seq, for the sequencing of small RNA with secondary structure close to the 5' end. This method accepts a variety of 5' modified RNA, including 5' monophosphorylated RNA, as demonstrated by the construction of a HeLa cell microRNA cDNA library.

High-throughput, Efficient, and Unbiased Capture of Small RNAs from Low-input Samples for Sequencing

MicroRNAs hold great promise as biomarkers of disease. However, there are few efficient and robust methods for measuring microRNAs from low input samples. Here, we develop a high-throughput sequencing protocol that efficiently captures small RNAs while minimizing inherent biases associated with library production. The protocol is based on early barcoding such that all downstream manipulations can be performed on a pool of many samples thereby reducing reagent usage and workload. We show that the optimization of adapter concentrations along with the addition of nucleotide modifications and random nucleotides increases the efficiency of small RNA capture. We further show, using unique molecular identifiers, that stochastic capture of low input RNA rather than PCR amplification influences the biased quantitation of intermediately and lowly expressed microRNAs. Our improved method allows the processing of tens to hundreds of samples simultaneously while retaining high efficiency quantitation of microRNAs in low input samples from tissues or bodily fluids.

The evolution of RNA structural probing methods: From gels to next-generation sequencing

RNA molecules are important players in all domains of life and the study of the relationship between their multiple flexible states and the associated biological roles has increased in recent years. For several decades, chemical and enzymatic structural probing experiments have been used to determine RNA structure. During this time, there has been a steady improvement in probing reagents and experimental methods, and today the structural biologist community has a large range of tools at its disposal to probe the secondary structure of RNAs in vitro and in cells. Early experiments used radioactive labeling and polyacrylamide gel electrophoresis as read-out methods. This was superseded by capillary electrophoresis, and more recently by next-generation sequencing. Today, powerful structural probing methods can characterize RNA structure on a genome-wide scale. In this review, we will provide an overview of RNA structural probing methodologies from a historical and technical perspective. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry RNA Methods > RNA Analyses in vitro and In Silico RNA Methods > RNA Analyses in Cells.

RATEmiRs: the rat atlas of tissue-specific and enriched miRNAs database

Background

MicroRNAs (miRNAs) regulate gene expression and have been targeted as indicators of environmental/toxicologic stressors. Using the data from our deep sequencing of miRNAs in an extensive sampling of rat tissues, we developed a database called RATEmiRs for the Rat Atlas of Tissue-specific and Enriched miRNAs to allow users to dynamically determine mature-, iso- and pre-miR expression abundance, enrichment and specificity in rat tissues and organs.

Results

Illumina sequencing count data from mapped reads and meta data from the miRNA body atlas consisting of 21 and 23 tissues (14 organs) of toxicologic interest from 12 to 13 week old male and female Sprague Dawley rats respectively, were managed in a relational database with a user-friendly query interface. Data-driven pipelines are available to tailor the identification of tissue-enriched (TE) and tissue-specific (TS) miRNAs. Data-driven organ-specific (OS) pipelines reveal miRNAs that are expressed predominately in a given organ. A user-driven approach is also available to assess the tissue expression of user-specified miRNAs. Using one tissue vs other tissues and tissue(s) of an organ vs other organs, we illustrate the utility of RATEmiRs to facilitate the identification of candidate miRNAs. As a use case example, RATEmiRs revealed two TS miRNAs in the liver: rno-miR-122-3p and rno-miR-122-5p. When liver is compared to just the brain tissues for example, rno-miR-192-5p, rno-miR-193-3p, rno-miR-203b-3p, rno-miR-3559-5p, rno-miR-802-3p and rno-miR-802-5p are also detected as abundantly expressed in liver. As another example, 55 miRNAs from the RATEmiRs query of ileum vs brain tissues overlapped with miRNAs identified from the same comparison of tissues in an independent, publicly available dataset of 10 week old male rat microarray data suggesting that these miRNAs are likely not age-specific, platform-specific nor pipeline-dependent. Lastly, we identified 10 miRNAs that have conserved tissue/organ-specific expression between the rat and human species.

Conclusions

RATEmiRs provides a new platform for identification of TE, TS and OS miRNAs in a broad array of rat tissues. RATEmiRs is available at: https://www.niehs.nih.gov/ratemirs

Evaluation of methodologies for microRNA biomarker detection by next generation sequencing

In recent years, microRNAs (miRNAs) in tissues and biofluids have emerged as a new class of promising biomarkers for numerous diseases. Blood-based biomarkers are particularly desirable since serum or plasma is easily accessible and can be sampled repeatedly. To comprehensively explore the biomarker potential of miRNAs, sensitive, accurate and cost-efficient miRNA profiling techniques are required. Next generation sequencing (NGS) is emerging as the preferred method for miRNA profiling; offering high sensitivity, single-nucleotide resolution and the possibility to profile a considerable number of samples in parallel. Despite the excitement about miRNA biomarkers, challenges associated with insufficient characterization of the sequencing library preparation efficacy, precision and method-related quantification bias have not been addressed in detail and are generally underappreciated in the wider research community. Here, we have tested in parallel four commercially available small RNA sequencing kits against a cohort of samples comprised of human plasma, human serum, murine brain tissue and a reference library containing ~ 950 synthetic miRNAs. We discuss the advantages and limits of these methodologies for massive parallel microRNAs profiling. This work can serve as guideline for choosing an adequate library preparation method, based on sensitivity, specificity and accuracy of miRNA quantification, workflow convenience and potential for automation.

Decreasing miRNA sequencing bias using a single adapter and circularization approach

The ability to accurately quantify all the microRNAs (miRNAs) in a sample is important for understanding miRNA biology and for development of new biomarkers and therapeutic targets. We develop a new method for preparing miRNA sequencing libraries, RealSeq®-AC, that involves ligating the miRNAs with a single adapter and circularizing the ligation products. When compared to other methods, RealSeq®-AC provides greatly reduced miRNA sequencing bias and allows the identification of the largest variety of miRNAs in biological samples. This reduced bias also allows robust quantification of miRNAs present in samples across a wide range of RNA input levels.

Comprehensive RNA-Sequencing Analysis in Serum and Muscle Reveals Novel Small RNA Signatures with Biomarker Potential for DMD

Extracellular small RNAs (sRNAs), including microRNAs (miRNAs), are promising biomarkers for diseases such as Duchenne muscular dystrophy (DMD), although their biological relevance is largely unknown. To investigate the relationship between intracellular and extracellular sRNA levels on a global scale, we performed sRNA sequencing in four muscle types and serum from wild-type, dystrophic mdx, and mdx mice in which dystrophin protein expression was restored by exon skipping. Differentially abundant sRNAs were identified in serum (mapping to miRNA, small nuclear RNA [snRNA], and PIWI-interacting RNA [piRNA] loci). One novel candidate biomarker, miR-483, was increased in both mdx serum and muscle, and also elevated in DMD patient sera. Dystrophin restoration induced global shifts in miRNA (including miR-483) and snRNA-fragment abundance toward wild-type levels. Specific serum piRNA-like sRNAs also responded to exon skipping therapy. Absolute miRNA expression in muscle was positively correlated with abundance in the circulation, although multiple highly expressed miRNAs in muscle were not elevated in mdx serum, suggesting that both passive and selective release mechanisms contribute to serum miRNA levels. In conclusion, this study has revealed new insights into the sRNA biology of dystrophin deficiency and identified novel DMD biomarkers.

Extracellular vesicles and their nucleic acids for biomarker discovery

Extracellular vesicles (EVs) are a heterogenous population of vesicles originate from cells. EVs are found in different biofluids and carry different macromolecules, including proteins, lipids, and nucleic acids, providing a snap shot of the parental cells at the time of release. EVs have the ability to transfer molecular cargoes to other cells and can initiate different physiological and pathological processes. Mounting lines of evidence demonstrated that EVs' cargo and machinery is affected in disease states, positioning EVs as potential sources for the discovery of novel biomarkers. In this review, we demonstrate a conceptual overview of the EV field with particular focus on their nucleic acid cargoes. Current knowledge of EV subtypes, nucleic acid cargo and pathophysiological roles are outlined, with emphasis placed on advantages against competing analytes. We review the utility of EVs and their nucleic acid cargoes as biomarkers and critically assess the newly available advances in the field of EV biomarkers and high throughput technologies. Challenges to achieving the diagnostic potential of EVs, including sample handling, EV isolation, methodological considerations, and bioassay reproducibility are discussed. Future implementation of 'omics-based technologies and integration of systems biology approaches for the development of EV-based biomarkers and personalized medicine are also considered.