Small RNA expression profiling by high-throughput sequencing: implications of enzymatic manipulation

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.

Isolation and Library Preparation of Planarian piRNAs

In planarian flatworms, piRNAs and SMEDWI (Schmidtea mediterranea PIWI) proteins are both essential for the animals' impressive regenerative ability and for their survival. A knockdown of SMEDWI proteins disrupts the specification of the planarian germline and impairs stem cell differentiation, resulting in lethal phenotypes. As the molecular targets of PIWI proteins and thus their biological function are determined by PIWI-bound small RNAs, termed piRNAs (for PIWI-interacting RNAs), it is imperative to study the wealth of PIWI-bound piRNAs using next-generation sequencing-based techniques. Prior to sequencing, piRNAs bound to individual SMEDWI proteins must be isolated. To that end, we established an immunoprecipitation protocol that can be applied to all planarian SMEDWI proteins. Co-immunoprecipitated piRNAs are visualized by using qualitative radioactive 5'-end labeling, which detects even trace amounts of small RNAs. Next, isolated piRNAs are subjected to a library preparation protocol that has been optimized for the efficient capture of piRNAs, whose 3'-ends carry a 2'-O-methyl modification. Successfully prepared piRNA libraries are subjected to Illumina-based next-generation sequencing. Obtained data are analyzed as presented in the accompanying manuscript.

The role of noncoding RNAs in Parkinson's disease: biomarkers and associations with pathogenic pathways

The discovery of various noncoding RNAs (ncRNAs) and their biological implications is a growing area in cell biology. Increasing evidence has revealed canonical and noncanonical functions of long and small ncRNAs, including microRNAs, long ncRNAs (lncRNAs), circular RNAs, PIWI-interacting RNAs, and tRNA-derived fragments. These ncRNAs have the ability to regulate gene expression and modify metabolic pathways. Thus, they may have important roles as diagnostic biomarkers or therapeutic targets in various diseases, including neurodegenerative disorders, especially Parkinson's disease. Recently, through diverse sequencing technologies and a wide variety of bioinformatic analytical tools, such as reverse transcriptase quantitative PCR, microarrays, next-generation sequencing and long-read sequencing, numerous ncRNAs have been shown to be associated with neurodegenerative disorders, including Parkinson's disease. In this review article, we will first introduce the biogenesis of different ncRNAs, including microRNAs, PIWI-interacting RNAs, circular RNAs, long noncoding RNAs, and tRNA-derived fragments. The pros and cons of the detection platforms of ncRNAs and the reproducibility of bioinformatic analytical tools will be discussed in the second part. Finally, the recent discovery of numerous PD-associated ncRNAs and their association with the diagnosis and pathophysiology of PD are reviewed, and microRNAs and long ncRNAs that are transported by exosomes in biofluids are particularly emphasized.

cyPhyRNA-seq: a genome-scale RNA-seq method to detect active self-cleaving ribozymes by capturing RNAs with 2',3' cyclic phosphates and 5' hydroxyl ends

Self-cleaving ribozymes are catalytically active RNAs that cleave themselves into a 5′-fragment with a 2′,3′-cyclic phosphate and a 3′-fragment with a 5′-hydroxyl. They are widely applied for the construction of synthetic RNA devices and RNA-based therapeutics. However, the targeted discovery of self-cleaving ribozymes remains a major challenge. We developed a transcriptome-wide method, called cyPhyRNA-seq, to screen for ribozyme cleavage fragments in total RNA extract. This approach employs the specific ligation-based capture of ribozyme 5′-fragments using a variant of the Arabidopsis thaliana tRNA ligase we engineered. To capture ribozyme 3′-fragments, they are enriched from total RNA by enzymatic treatments. We optimized and enhanced the individual steps of cyPhyRNA-seq in vitro and in spike-in experiments. Then, we applied cyPhyRNA-seq to total RNA isolated from the bacterium Desulfovibrio vulgaris and detected self-cleavage of the three predicted type II hammerhead ribozymes, whose activity had not been examined to date. cyPhyRNA-seq can be used for the global analysis of active self-cleaving ribozymes with the advantage to capture both ribozyme cleavage fragments from total RNA. Especially in organisms harbouring many self-cleaving RNAs, cyPhyRNA-seq facilitates the investigation of cleavage activity. Moreover, this method has the potential to be used to discover novel self-cleaving ribozymes in different organisms.

Optimization of small RNA extraction and comparative study of NGS library preparation from low count sperm samples

Recent studies demonstrate that sperm epigenome is a vehicle that conveys paternal experiences to offspring phenotype. That evidence triggers interest of both experimental and epidemiological studies of epigenetic markers in sperm. Since samples are often unique in epidemiological studies, a careful and efficient use of the material is a critical requirement. The goal of this study was to provide optimization of methods for the isolation of small RNAs from spermatozoa and library preparation for sequencing. A total 67 fractionated sperm samples from the Russian Children's Study biobank prospectively collected at 18-20 years of age were used to isolate small RNAs with median (IQR) input total sperm count 17.0 (7.4-35.9) million. Twenty-four pairs of libraries were prepared using the NEBNext and NEXTFlex kits, 19 libraries using NEBNext and 6 using NEXTFlex. All libraries were sequenced on NextSeq 500, and the results were evaluated as a function of the number of small non-coding RNA (sncRNA) detected, quality parameters of sequencing libraries, as well as technical features of sample preparation. Although the same amount of miRNA input was used for NEBNext and NEXTFlex libraries, the concentration of DNA in NEBNext libraries was significantly higher in comparison with NEXTFlex libraries. In high input (sperm count >28 million and more than 25 ng miRNA in library) NEXTFlex Small RNA-Seq kit detected more microRNAs. In low input, the NEBNext proved more effective. The tricks and traps to protocol optimization are presented, including an efficient and effector gel-based system for the removal of sequencing library adaptors.

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.

Small RNAs Are Implicated in Regulation of Gene and Transposable Element Expression in the Protist Trichomonas vaginalis

Trichomoniasis, caused by the protozoan Trichomonas vaginalis, is the most common nonviral sexually transmitted infection in humans. The millions of cases each year have sequelae that may include complications during pregnancy and increased risk of HIV infection.

Trichomonas vaginalis is the causative agent of trichomoniasis, the most prevalent nonviral sexually transmitted infection worldwide. Repetitive elements, including transposable elements (TEs) and virally derived repeats, comprise more than half of the ∼160-Mb T. vaginalis genome. An intriguing question is how the parasite controls its potentially lethal complement of mobile elements, which can disrupt transcription of protein-coding genes and genome functions. In this study, we generated high-throughput RNA sequencing (RNA-Seq) and small RNA-Seq data sets in triplicate for the T. vaginalis G3 reference strain and characterized the mRNA and small RNA populations and their mapping patterns along all six chromosomes. Mapping the RNA-Seq transcripts to the genome revealed that the majority of genes predicted within repetitive elements are not expressed. Interestingly, we identified a novel species of small RNA that maps bidirectionally along the chromosomes and is correlated with reduced protein-coding gene expression and reduced RNA-Seq coverage in repetitive elements. This novel small RNA family may play a regulatory role in gene and repetitive element expression. Our results identify a possible small RNA pathway mechanism by which the parasite regulates expression of genes and TEs and raise intriguing questions as to the role repeats may play in shaping T. vaginalis genome evolution and the diversity of small RNA pathways in general.

IMPORTANCE Trichomoniasis, caused by the protozoan Trichomonas vaginalis, is the most common nonviral sexually transmitted infection in humans. The millions of cases each year have sequelae that may include complications during pregnancy and increased risk of HIV infection. Given its evident success in this niche, it is paradoxical that T. vaginalis harbors in its genome thousands of transposable elements that have the potential to be extremely detrimental to normal genomic function. In many organisms, transposon expression is regulated by the activity of endogenously expressed short (∼21 to 35 nucleotides [nt]) small RNA molecules that effect gene silencing by targeting mRNAs for degradation or by recruiting epigenetic silencing machinery to locations in the genome. Our research has identified small RNA molecules correlated with reduced expression of T. vaginalis genes and transposons. This suggests that a small RNA pathway is a major contributor to gene expression patterns in the parasite and opens up new avenues for investigation into small RNA biogenesis, function, and diversity.

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.

DTA-SiST: de novo transcriptome assembly by using simplified suffix trees

Background

Alternative splicing allows the pre-mRNAs of a gene to be spliced into various mRNAs, which greatly increases the diversity of proteins. High-throughput sequencing of mRNAs has revolutionized our ability for transcripts reconstruction. However, the massive size of short reads makes de novo transcripts assembly an algorithmic challenge.

Results

We develop a novel radical framework, called DTA-SiST, for de novo transcriptome assembly based on suffix trees. DTA-SiST first extends contigs by reads that have the longest overlaps with the contigs’ terminuses. These reads can be found in linear time of the lengths of the reads through a well-designed suffix tree structure. Then, DTA-SiST constructs splicing graphs based on contigs for each gene locus. Finally, DTA-SiST proposes two strategies to extract transcript-representing paths: a depth-first enumeration strategy and a hybrid strategy based on length and coverage. We implemented the above two strategies and compared them with the state-of-the-art de novo assemblers on both simulated and real datasets. Experimental results showed that the depth-first enumeration strategy performs always better with recall and also better with precision for smaller datasets while the hybrid strategy leads with precision for big datasets.

Conclusions

DTA-SiST performs more competitive than the other compared de novo assemblers especially with precision measure, due to the read-based contig extension strategy and the elegant transcripts extraction rules.

Extracellular miRNAs: From Biomarkers to Mediators of Physiology and Disease

miRNAs can be found in serum and other body fluids and serve as biomarkers for disease. More importantly, secreted miRNAs, especially those in extracellular vesicles (EVs) such as exosomes, may mediate paracrine and endocrine communication between different tissues and thus modulate gene expression and the function of distal cells. When impaired, these processes can lead to tissue dysfunction, aging, and disease. Adipose tissue is an especially important contributor to the pool of circulating exosomal miRNAs. As a result, alterations in adipose tissue mass or function, which occur in many metabolic conditions, can lead to changes in circulating miRNAs, which then function systemically. Here we review the findings that led to these conclusions and discuss how this sets the stage for new lines of investigation in which extracellular miRNAs are recognized as important mediators of intercellular communication and potential candidates for therapy of disease.

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.

Stage-dependent piRNAs in chicken implicated roles in modulating male germ cell development

Background

The PIWI/piRNA pathway is a conserved machinery important for germ cell development and fertility. This piRNA-guided molecular machinery is best known for repressing derepressed transposable elements (TE) during epigenomic reprogramming. The extent to which piRNAs are involved in modulating transcripts beyond TEs still need to be clarified, and it may be a stage-dependent event. We chose chicken germline as a study model because of the significantly lower TE complexity in the chicken genome compared to mammalian species.

Results

We generated high-confidence piRNA candidates in various stages across chicken germline development by 3′-end-methylation-enriched small RNA sequencing and in-house bioinformatics analysis. We observed a significant developmental stage-dependent loss of TE association and a shifting of the ping-pong cycle signatures. Moreover, the stage-dependent reciprocal abundance of LINE retrotransposons, CR1-C, and its associated piRNAs implicated the developmental stage-dependent role of piRNA machinery. The stage dependency of piRNA expression and its potential functions can be better addressed by analyzing the piRNA precursors/clusters. Interestingly, the new piRNA clusters identified from embryonic chicken testes revealed evolutionary conservation between chickens and mammals, which was previously thought to not exist.

Conclusions

In this report, we provided an original chicken RNA resource and proposed an analytical methodology that can be used to investigate stage-dependent changes in piRNA compositions and their potential roles in TE regulation and beyond, and also revealed possible conserved functions of piRNAs in developing germ cells.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4820-9) contains supplementary material, which is available to authorized users.

Capturing the 'ome': the expanding molecular toolbox for RNA and DNA library construction

All sequencing experiments and most functional genomics screens rely on the generation of libraries to comprehensively capture pools of targeted sequences. In the past decade especially, driven by the progress in the field of massively parallel sequencing, numerous studies have comprehensively assessed the impact of particular manipulations on library complexity and quality, and characterized the activities and specificities of several key enzymes used in library construction. Fortunately, careful protocol design and reagent choice can substantially mitigate many of these biases, and enable reliable representation of sequences in libraries. This review aims to guide the reader through the vast expanse of literature on the subject to promote informed library generation, independent of the application.

Systematic comparison of small RNA library preparation protocols for next-generation sequencing

Background

Next-generation sequencing technologies have revolutionized the study of small RNAs (sRNAs) on a genome-wide scale. However, classical sRNA library preparation methods introduce serious bias, mainly during adapter ligation steps. Several types of sRNA including plant microRNAs (miRNA), piwi-interacting RNAs (piRNA) in insects, nematodes and mammals, and small interfering RNAs (siRNA) in insects and plants contain a 2’-O-methyl (2’-OMe) modification at their 3′ terminal nucleotide. This inhibits 3′ adapter ligation and makes library preparation particularly challenging. To reduce bias, the NEBNext kit (New England Biolabs) uses polyethylene glycol (PEG), the NEXTflex V2 kit (BIOO Scientific) uses both randomised adapters and PEG, and the novel SMARTer (Clontech) and CATS (Diagenode) kits avoid ligation altogether. Here we compared these methods with Illumina’s classical TruSeq protocol regarding the detection of normal and 2’ OMe RNAs. In addition, we modified the TruSeq and NEXTflex protocols to identify conditions that improve performance.

Results

Among the five kits tested with their respective standard protocols, the SMARTer and CATS kits had the lowest levels of bias but also had a strong formation of side products, and as a result performed relatively poorly with biological samples; NEXTflex detected the largest numbers of different miRNAs. The use of a novel type of randomised adapters called MidRand-Like (MRL) adapters and PEG improved the detection of 2’ OMe RNAs both in the TruSeq as well as in the NEXTflex protocol.

Conclusions

While it is commonly accepted that biases in sRNA library preparation protocols are mainly due to adapter ligation steps, the ligation-free protocols were not the best performing methods. Our modified versions of the TruSeq and NEXTflex protocols provide an improved tool for the study of 2’ OMe RNAs.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4491-6) contains supplementary material, which is available to authorized users.

Coding and noncoding landscape of extracellular RNA released by human glioma stem cells

Tumor-released RNA may mediate intercellular communication and serve as biomarkers. Here we develop a protocol enabling quantitative, minimally biased analysis of extracellular RNAs (exRNAs) associated with microvesicles, exosomes (collectively called EVs), and ribonucleoproteins (RNPs). The exRNA complexes isolated from patient-derived glioma stem-like cultures exhibit distinct compositions, with microvesicles most closely reflecting cellular transcriptome. exRNA is enriched in small ncRNAs, such as miRNAs in exosomes, and precisely processed tRNA and Y RNA fragments in EVs and exRNPs. EV-enclosed mRNAs are mostly fragmented, and UTRs enriched; nevertheless, some full-length mRNAs are present. Overall, there is less than one copy of non-rRNA per EV. Our results suggest that massive EV/exRNA uptake would be required to ensure functional impact of transferred RNA on brain recipient cells and predict the most impactful miRNAs in such conditions. This study also provides a catalog of diverse exRNAs useful for biomarker discovery and validates its feasibility on cerebrospinal fluid.

While circulating DNA has been extensively explored as a potential cancer biomarker, RNA potential has been overlooked so far. Here the authors present a comprehensive analysis of extracellular RNA secreted by glioblastoma cells that could prove a valuable resource for biomarker discovery and a means of intercellular communication.

Frequent hypermethylation of orphan CpG islands with enhancer activity in cancer

BACKGROUND

CpG islands (CGIs) are interspersed DNA sequences that have unusually high CpG ratios and GC contents. CGIs are typically located in the promoter of protein-coding genes. They normally lack DNA methylation but become hypermethylated and induce repression of associated genes in cancer. However, the biological functions of non-promoter CGIs (orphan CGIs) largely remain unclear.

RESULTS

Here, we identify orphan CGIs that do not map to the promoter of any protein-coding or non-coding transcripts but possess chromatin and transcriptional marks that reflect enhancer activity (termed eCGIs). They exhibit three-dimensional chromatin looping toward multiple target genes with high affinity. Intriguingly, transcription regulators were frequently associated with such CGI-containing enhancers. Remarkably, our analyses in cell lines and clinical tissues showed that eCGIs have more dynamic DNA methylation changes in cancer relative to promoter CGIs. The observed eCGI hypermethylation was accompanied by a loss of enhancer marks and transcriptional inactivation of the target genes.

CONCLUSION

Our results suggest that eCGIs may constitute a distinct class of enhancers and perform a more instrumental role in tumorigenesis than typical CGIs in gene promoters.

Sources and Functions of Extracellular Small RNAs in Human Circulation

Various biotypes of endogenous small RNAs (sRNAs) have been detected in human circulation, including microRNAs, transfer RNAs, ribosomal RNA, and yRNA fragments. These extracellular sRNAs (ex-sRNAs) are packaged and secreted by many different cell types. Ex-sRNAs exhibit differences in abundance in several disease states and have, therefore, been proposed for use as effective biomarkers. Furthermore, exosome-borne ex-sRNAs have been reported to elicit physiological responses in acceptor cells. Exogenous ex-sRNAs derived from diet (most prominently from plants) and microorganisms have also been reported in human blood. Essential issues that remain to be conclusively addressed concern the (a) presence and sources of exogenous ex-sRNAs in human bodily fluids, (b) detection and measurement of ex-sRNAs in human circulation, (c) selectivity of ex-sRNA export and import, (d) sensitivity and specificity of ex-sRNA delivery to cellular targets, and (e) cell-, tissue-, organ-, and organism-wide impacts of ex-sRNA-mediated cell-to-cell communication. We survey the present state of knowledge of most of these issues in this review.

Toward reliable biomarker signatures in the age of liquid biopsies - how to standardize the small RNA-Seq workflow

Small RNA-Seq has emerged as a powerful tool in transcriptomics, gene expression profiling and biomarker discovery. Sequencing cell-free nucleic acids, particularly microRNA (miRNA), from liquid biopsies additionally provides exciting possibilities for molecular diagnostics, and might help establish disease-specific biomarker signatures. The complexity of the small RNA-Seq workflow, however, bears challenges and biases that researchers need to be aware of in order to generate high-quality data. Rigorous standardization and extensive validation are required to guarantee reliability, reproducibility and comparability of research findings. Hypotheses based on flawed experimental conditions can be inconsistent and even misleading. Comparable to the well-established MIQE guidelines for qPCR experiments, this work aims at establishing guidelines for experimental design and pre-analytical sample processing, standardization of library preparation and sequencing reactions, as well as facilitating data analysis. We highlight bottlenecks in small RNA-Seq experiments, point out the importance of stringent quality control and validation, and provide a primer for differential expression analysis and biomarker discovery. Following our recommendations will encourage better sequencing practice, increase experimental transparency and lead to more reproducible small RNA-Seq results. This will ultimately enhance the validity of biomarker signatures, and allow reliable and robust clinical predictions.

Where does transcription start? 5'-RACE adapted to next-generation sequencing

The variability and complexity of the transcription initiation process was examined by adapting RNA ligase-mediated rapid amplification of 5' cDNA ends (5'-RACE) to Next-Generation Sequencing (NGS). We oligo-labelled 5'-m(7)G-capped mRNA from two genes, the simple mono-exonic Beta-2-Adrenoceptor (ADRB2R)and the complex multi-exonic Glucocorticoid Receptor (GR, NR3C1), and detected a variability in TSS location that has received little attention up to now. Transcription was not initiated at a fixed TSS, but from loci of 4 to 10 adjacent nucleotides. Individual TSSs had frequencies from <0.001% to 38.5% of the total gene-specific 5' m(7)G-capped transcripts. ADRB2R used a single locus consisting of 4 adjacent TSSs. Unstimulated, the GR used a total of 358 TSSs distributed throughout 38 loci, that were principally in the 5' UTRs and were spliced using established donor and acceptor sites. Complete demethylation of the epigenetically sensitive GR promoter with 5-azacytidine induced one new locus and 127 TSSs, 12 of which were unique. We induced GR transcription with dexamethasone and Interferon-γ, adding one new locus and 185 additional TSSs distributed throughout the promoter region. In-vitro the TSS microvariability regulated mRNA translation efficiency and the relative abundance of the different GRN-terminal protein isoform levels.

Bias in ligation-based small RNA sequencing library construction is determined by adaptor and RNA structure

High-throughput sequencing (HTS) has become a powerful tool for the detection of and sequence characterization of microRNAs (miRNA) and other small RNAs (sRNA). Unfortunately, the use of HTS data to determine the relative quantity of different miRNAs in a sample has been shown to be inconsistent with quantitative PCR and Northern Blot results. Several recent studies have concluded that the major contributor to this inconsistency is bias introduced during the construction of sRNA libraries for HTS and that the bias is primarily derived from the adaptor ligation steps, specifically where single stranded adaptors are sequentially ligated to the 3’ and 5’-end of sRNAs using T4 RNA ligases. In this study we investigated the effects of ligation bias by using a pool of randomized ligation substrates, defined mixtures of miRNA sequences and several combinations of adaptors in HTS library construction. We show that like the 3’ adaptor ligation step, the 5’ adaptor ligation is also biased, not because of primary sequence, but instead due to secondary structures of the two ligation substrates. We find that multiple secondary structural factors influence final representation in HTS results. Our results provide insight about the nature of ligation bias and allowed us to design adaptors that reduce ligation bias and produce HTS results that more accurately reflect the actual concentrations of miRNAs in the defined starting material.

A non-radioactive method for small RNA detection by northern blotting

BACKGROUND

Small non-coding RNAs are essential regulators of gene expression at the transcriptional and posttranscriptional levels. High-throughput sequencing has revealed thousands of predicted small RNAs; however, only a few of these have been well characterized. Northern blotting is the most convincing method for small RNA validation.

FINDINGS

In this study, we improved the Northern blot method by using biotin-labeled probes. miRNAs and siRNAs derived from both Arabidopsis thaliana and Oryza sativa were investigated. The results suggest that this improved method is sensitive and efficient, with approximately 5 μg of total RNA being sufficient for detection. Furthermore, long-term storage of probes labeled in this manner is more convenient, less contaminative and degradative compared with traditional probes.

CONCLUSIONS

This protocol is an alternative strategy for small RNA detection and represents an efficient means of researching small RNAs.

Evaluating bias-reducing protocols for RNA sequencing library preparation

Background

Next-generation sequencing does not yield fully unbiased estimates for read abundance, which may impact on the conclusions that can be drawn from sequencing data. The ligation step in RNA sequencing library generation is a known source of bias, motivating developments in enzyme technology and library construction protocols. We present the first comparison of the standard duplex adaptor protocol supplied by Life Technologies for use on the Ion Torrent PGM with an alternate single adaptor approach involving CircLigase (CircLig protocol).

A correlation between over-representation in sequenced libraries and degree of secondary structure has been reported previously, therefore we also investigated whether bias could be reduced by ligation with an enzyme that functions at a temperature not permissive for such structure.

Results

A pool of small RNA fragments of known composition was converted into a sequencing library using one of three protocols and sequenced on an Ion Torrent PGM. The CircLig protocol resulted in less over-representation of specific sequences than the standard protocol. Over-represented sequences are more likely to be predicted to have secondary structure and to co-fold with adaptor sequences. However, use of the thermostable ligase Methanobacterium thermoautotrophicum RNA ligase K97A (Mth K97A) was not sufficient to reduce bias.

Conclusions

The single adaptor CircLigase-based approach significantly reduces, but does not eliminate, bias in Ion Torrent data. Ligases that function at temperatures to remove the possible influence of secondary structure on library generation may be of value, although Mth K97A is not effective in this case.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-569) contains supplementary material, which is available to authorized users.

Terminator oligo blocking efficiently eliminates rRNA from Drosophila small RNA sequencing libraries

A large number of methods are available to deplete ribosomal RNA reads from high-throughput RNA sequencing experiments. Such methods are critical for sequencing Drosophila small RNAs between 20 and 30 nucleotides because size selection is not typically sufficient to exclude the highly abundant class of 30 nucleotide 2S rRNA. Here we demonstrate that pre-annealing terminator oligos complimentary to Drosophila 2S rRNA prior to 5' adapter ligation and reverse transcription efficiently depletes 2S rRNA sequences from the sequencing reaction in a simple and inexpensive way. This depletion is highly specific and is achieved with minimal perturbation of miRNA and piRNA profiles.

Small RNA cloning and sequencing strategy affects host and viral microRNA expression signatures

The establishment of the microRNA (miRNA) expression signatures is the basic element to investigate the role played by these regulatory molecules in the biology of an organism. Marek's disease virus 1 (MDV-1) is an avian herpesvirus that naturally infects chicken and induces T cells lymphomas. During latency, MDV-1, like other herpesviruses, expresses a limited subset of transcripts. These include three miRNA clusters. Several studies identified the expression of virus and host encoded miRNAs from MDV-1 infected cell cultures and chickens. But a high discrepancy was observed when miRNA cloning frequencies obtained from different cloning and sequencing protocols were compared. Thus, we analyzed the effect of small RNA library preparation and sequencing on the miRNA frequencies obtained from the same RNA samples collected during MDV-1 infection of chicken at different steps of the oncoviral pathogenesis. Qualitative and quantitative variations were found in the data, depending on the strategy used. One of the mature miRNA derived from the latency-associated-transcript (LAT), mdv1-miR-M7-5p, showed the highest variation. Its cloning frequency was 50% of the viral miRNA counts when a small scale sequencing approach was used. Its frequency was 100 times less abundant when determined through the deep sequencing approach. Northern blot analysis showed a better correlation with the miRNA frequencies found by the small scale sequencing approach. By analyzing the cellular miRNA repertoire, we also found a gap between the two sequencing approaches. Collectively, our study indicates that next-generation sequencing data considered alone are limited for assessing the absolute copy number of transcripts. Thus, the quantification of small RNA should be addressed by compiling data obtained by using different techniques such as microarrays, qRT-PCR and NB analysis in support of high throughput sequencing data. These observations should be considered when miRNA variations are studied prior addressing functional studies.

Library preparation methods for next-generation sequencing: tone down the bias

Next-generation sequencing (NGS) has caused a revolution in biology. NGS requires the preparation of libraries in which (fragments of) DNA or RNA molecules are fused with adapters followed by PCR amplification and sequencing. It is evident that robust library preparation methods that produce a representative, non-biased source of nucleic acid material from the genome under investigation are of crucial importance. Nevertheless, it has become clear that NGS libraries for all types of applications contain biases that compromise the quality of NGS datasets and can lead to their erroneous interpretation. A detailed knowledge of the nature of these biases will be essential for a careful interpretation of NGS data on the one hand and will help to find ways to improve library quality or to develop bioinformatics tools to compensate for the bias on the other hand. In this review we discuss the literature on bias in the most common NGS library preparation protocols, both for DNA sequencing (DNA-seq) as well as for RNA sequencing (RNA-seq). Strikingly, almost all steps of the various protocols have been reported to introduce bias, especially in the case of RNA-seq, which is technically more challenging than DNA-seq. For each type of bias we discuss methods for improvement with a view to providing some useful advice to the researcher who wishes to convert any kind of raw nucleic acid into an NGS library.

Small RNA library construction for high-throughput sequencing

Since their discovery about 20 years ago, small RNAs have been shown to play a critical role in a myriad of biological processes. The greater availability of high-throughput sequencing has been invaluable to furthering our understanding of small RNAs as regulatory molecules. In particular, these sequencing technologies have been crucial in understanding the role of small RNAs in reproductive tissues, where millions of individual sequences are generated. In this context, high-throughput sequencing provides the requisite level of resolution that other procedures, like northern blotting, would not be able to achieve. Here, we describe a protocol for the preparation of small RNA libraries for sequencing using the Solexa/Illumina technology.

Novel small RNA expression libraries uncover hsa-miR-30b and hsa-miR-30c as important factors in anoikis resistance

MicroRNAs (miRNAs) have been widely studied in order to elucidate their biological functions. MicroRNA microarrays or miRNA overexpression libraries generated by synthesis and cloning of individual miRNAs have been used to study their different roles. In this work, we have developed a novel methodology to express mature miRNAs and other small RNAs from a double convergent RNA polymerase III promoter. We show that the generated miRNAs function similarly to those processed from primary transcripts or pri-miRNAs. This system allowed us to produce a lentiviral library expressing the whole population of small RNAs present in a metastatic cell line. A functional screening using this library led to the identification of hsa-miR-30b and hsa-miR-30c as negative regulators of cell death induced by loss of attachment (anoikis). Importantly, we demonstrated that the acquisition of anoikis resistance via these miRNAs is achieved through down-regulation of caspase 3 expression. Moreover, overexpression of these miRNAs resulted in a decrease of other types of caspase 3-dependent cell death and enhanced the survival of MCF10A acinar cells in morphogenesis assays, suggesting a putative role as oncomirs. In summary, this novel methodology provides a powerful and effective way for identifying novel small RNAs involved in a particular biological process.

Biases in small RNA deep sequencing data

High-throughput RNA sequencing (RNA-seq) is considered a powerful tool for novel gene discovery and fine-tuned transcriptional profiling. The digital nature of RNA-seq is also believed to simplify meta-analysis and to reduce background noise associated with hybridization-based approaches. The development of multiplex sequencing enables efficient and economic parallel analysis of gene expression. In addition, RNA-seq is of particular value when low RNA expression or modest changes between samples are monitored. However, recent data uncovered severe bias in the sequencing of small non-protein coding RNA (small RNA-seq or sRNA-seq), such that the expression levels of some RNAs appeared to be artificially enhanced and others diminished or even undetectable. The use of different adapters and barcodes during ligation as well as complex RNA structures and modifications drastically influence cDNA synthesis efficacies and exemplify sources of bias in deep sequencing. In addition, variable specific RNA G/C-content is associated with unequal polymerase chain reaction amplification efficiencies. Given the central importance of RNA-seq to molecular biology and personalized medicine, we review recent findings that challenge small non-protein coding RNA-seq data and suggest approaches and precautions to overcome or minimize bias.