Efficient microRNA capture and bar-coding via enzymatic oligonucleotide adenylation

Enhancing Transcriptome Mapping with Rapid PRO-seq Profiling of Nascent RNA

Precision nuclear run-on (PRO) sequencing (PRO-seq) is a powerful technique for mapping polymerase active sites with nucleotide resolution and measuring newly synthesized transcripts at both promoters and enhancer elements. The current PRO-seq protocol is time-intensive, technically challenging, and requires a large amount of starting material. To overcome these limitations, we developed rapid PRO-seq (rPRO-seq) which utilizes pre-adenylated single-stranded DNAs (AppDNA), a dimer blocking oligonucleotide (DBO), on-bead 5' RNA end repair, and column-based purification. These modifications enabled efficient transcriptome mapping within a single day (∼12 hours) increasing ligation efficiency, abolished adapter dimers, and reduced sample loss and RNA degradation. We demonstrate the reproducibility of rPRO-seq in measuring polymerases at promoters, gene bodies, and enhancers as compared to original PRO-seq protocols. Additionally, rPRO-seq is scalable, allowing for transcriptome mapping with as little as 25,000 cells. We apply rPRO-seq to study the role of Integrator in mouse hematopoietic stem and progenitor cell (mHSPC) homeostasis, identifying Ints11 as an essential component of transcriptional regulation and RNA processing in mHSPC homeostasis. Overall, rPRO-seq represents a significant advance in the field of nascent transcript analyses and will be a valuable tool for generating patient-specific genome-wide transcription profiles with minimal sample requirements.

Formation of an instantaneous nick for highly efficient adenylation of oligonucleotides by ligase without subsequent jointing

By forming a nick at the adenylation site instantaneously, nucleic acids are efficiently adenylated by T4 DNA ligase. The subsequent ligation is successfully suppressed in terms of rapid conversion of the instantaneous nick to a more stable gap. It is helpful to understand enzymatic ligation dynamics, and the adenylated products can be used for various practical applications.

Direct adenylation from 5'-OH-terminated oligonucleotides by a fusion enzyme containing Pfu RNA ligase and T4 polynucleotide kinase

5′-Adenylated oligonucleotides (AppOligos) are widely used for single-stranded DNA/RNA ligation in next-generation sequencing (NGS) applications such as microRNA (miRNA) profiling. The ligation between an AppOligo adapter and target molecules (such as miRNA) no longer requires ATP, thereby minimizing potential self-ligations and simplifying library preparation procedures. AppOligos can be produced by chemical synthesis or enzymatic modification. However, adenylation via chemical synthesis is inefficient and expensive, while enzymatic modification requires pre-phosphorylated substrate and additional purification. Here we cloned and characterized the Pfu RNA ligase encoded by the PF0353 gene in the hyperthermophilic archaea Pyrococcus furiosus. We further engineered fusion enzymes containing both Pfu RNA ligase and T4 polynucleotide kinase. One fusion enzyme, 8H-AP, was thermostable and can directly catalyze 5′-OH-terminated DNA substrates to adenylated products. The newly discovered Pfu RNA ligase and the engineered fusion enzyme may be useful tools for applications using AppOligos.

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

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

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.

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.

Capture, amplification, and global profiling of microRNAs from low quantities of whole cell lysate

MicroRNAs (miRNAs) are small non-coding RNAs that control gene expression at the post-transcriptional level via a complex regulatory network that requires genome-wide miRNA profiling to dissect. The patterns of miRNA expression at the genome scale are rich in diagnostic and prognostic information for human diseases such as cancers. This analysis, however, requires multi-step purification of RNAs from large quantities of cells, which is not only time consuming and costly but also challenging in situations where cell numbers are limited. In this study, we report direct capture, amplification, and library preparation of miRNAs from whole cell lysate without the need of pre-purification. As a result, it enables genome-wide miRNA profiling reproducibly with low quantity of cell samples (∼500 hematopoietic cells). Specifically, we conducted a systematic investigation of two key steps - cell lysis for miRNA release and 3' adaptor ligation required for direct miRNA capture and amplification. The obtained expression profile not only distinguishes cell types but also detects individual miRNA alterations in closely related isogenic cell lines. This approach, which is substantially simple as compared to the standard methods because of elimination of the need for RNA purification, is advantageous for the measurement of low quantity samples.

Holistic and Affordable Analyses of MicroRNA Expression Profiles Using Tagged cDNA Libraries and a Multiplex Sequencing Strategy

Small and long noncoding RNAs (ncRNAs) are key regulators of gene expression. Variations in ncRNA expression patterns can consequently affect the control of many cellular processes. Not just since 2006, when Andrew Z Fire and Craig C Mello were jointly awarded The Nobel Prize in Physiology or Medicine for their discovery of RNA interference, great efforts were undertaken to unleash the biomedical applicability of small noncoding RNAs, in particular microRNAs. With the technological evolution of massive parallel sequencing technologies over the last years, which now are available for an increasing number of scientists, there is a demand for comprehensible and efficient workflows reliable even for unique and valuable clinical specimens. Here we describe a highly reproducible low-cost protocol for analyses of miRNA expression patterns using tagged cDNA libraries and a multiplex sequencing strategy following an Illumina-like protocol. This protocol easily allows the identification of expression differences from samples of tissues of 1-2 mm³ and fluids of 50-200 μL. We further provide entry points into useful computational biology applications, whose target groups explicitly involve non-bioinformaticians.

Capture and sequencing of NAD-capped RNA sequences with NAD captureSeq

Here we describe a protocol for NAD captureSeq that allows for the identification of nicotinamide-adenine dinucleotide (NAD)-capped RNA sequences in total RNA samples from different organisms. NAD-capped RNA is first chemo-enzymatically biotinylated with high efficiency, permitting selective capture on streptavidin beads. Then, a highly efficient library preparation protocol tailored to immobilized, 5'-modified RNA is applied, with adaptor ligation to the RNA's 3' terminus and reverse transcription (RT) performed on-bead. Then, cDNA is released into solution, tailed, ligated to a second adaptor and PCR-amplified. After next-generation sequencing (NGS) of the DNA library, enriched sequences are identified by comparison with a control sample in which the first step of chemo-enzymatic biotinylation is omitted. Because the downstream protocol does not necessarily rely on NAD-modified but on 'clickable' or biotin-modified RNA, it can be applied to other RNA modifications or RNA-biomolecule interactions. The central part of this protocol can be completed in ∼7 d, excluding preparatory steps, sequencing and bioinformatic analysis.

Substantial and robust changes in microRNA transcriptome support postnatal development of the hypothalamus in rat

MicroRNAs (miRNAs) modulate gene expression in male germ cells and somatic tissues of mammals on a genome-wide scale. Hundreds of miRNAs are encoded by mammalian genomes, a large fraction of which is expressed in brain. Here we have investigated the complexity and dynamics of miRNA transcriptomes that associate with neuronal network maturation of hypothalamic arcuate nucleus and median eminence (ARC/ME) in rat by analysing more than 300 miRNAs from 3-7 biological replicates at 5 postnatal time-points. The network connecting ARC/ME to other hypothalamic and extra-hypothalamic regions maturates in an environment-dependent manner. We therefore analyzed miRNA transcriptomes of progeny of dams fed either a balanced or unbalanced diet during gestation and lactation. More than 30% of the miRNAs displayed significative changes of expression between stages P8 and P14, and P21 and P28; half of the changes were greater than 3-fold. Among those miRNAs were well-known and dozens of still poorly documented miRNAs. Progeny of dams fed an unbanced diet displayed a severe growth retardation phenotype, lower levels of plasma leptin but almost identical miRNA transcriptomes. Together these data demonstrate that two substantial and robust changes in miRNA transcriptome of ARC/ME occur at a period crucial for neuronal network functional organization.

Adenylylation of small RNA sequencing adapters using the TS2126 RNA ligase I

Many high-throughput small RNA next-generation sequencing protocols use 5' preadenylylated DNA oligonucleotide adapters during cDNA library preparation. Preadenylylation of the DNA adapter's 5' end frees from ATP-dependence the ligation of the adapter to RNA collections, thereby avoiding ATP-dependent side reactions. However, preadenylylation of the DNA adapters can be costly and difficult. The currently available method for chemical adenylylation of DNA adapters is inefficient and uses techniques not typically practiced in laboratories profiling cellular RNA expression. An alternative enzymatic method using a commercial RNA ligase was recently introduced, but this enzyme works best as a stoichiometric adenylylating reagent rather than a catalyst and can therefore prove costly when several variant adapters are needed or during scale-up or high-throughput adenylylation procedures. Here, we describe a simple, scalable, and highly efficient method for the 5' adenylylation of DNA oligonucleotides using the thermostable RNA ligase 1 from bacteriophage TS2126. Adapters with 3' blocking groups are adenylylated at >95% yield at catalytic enzyme-to-adapter ratios and need not be gel purified before ligation to RNA acceptors. Experimental conditions are also reported that enable DNA adapters with free 3' ends to be 5' adenylylated at >90% efficiency.

Efficient synthesis of stably adenylated DNA and RNA adapters for microRNA capture using T4 RNA ligase 1

MicroRNA profiling methods have become increasingly important due to the rapid rise of microRNA in both basic and translational sciences. A critical step in many microRNA profiling assays is adapter ligation using pre-adenylated adapters. While pre-adenylated adapters can be chemically or enzymatically prepared, enzymatic adenylation is preferred due to its ease and high yield. However, previously reported enzymatic methods either require tedious purification steps or use thermostable ligases that can generate side products during the subsequent ligation step. We have developed a highly efficient, template- and purification-free, adapter adenylation method using T4 RNA ligase 1. This method is capable of adenylating large amounts of adapter at ~100% efficiency and can efficiently adenylate both DNA and RNA bases. We find that the adenylation reaction speed can differ between DNA and RNA and between terminal nucleotides, leading to bias if reactions are not allowed to run to completion. We further find that the addition of high PEG levels can effectively suppress these differences.

Blocking of targeted microRNAs from next-generation sequencing libraries

Highly abundant microRNAs (miRNAs) in small RNA sequencing libraries make it difficult to obtain efficient measurements of more lowly expressed species. We present a new method that allows for the selective blocking of specific, abundant miRNAs during preparation of sequencing libraries. This technique is specific with little off-target effects and has no impact on the reproducibility of the measurement of non-targeted species. In human plasma samples, we demonstrate that blocking of highly abundant hsa-miR-16–5p leads to improved detection of lowly expressed miRNAs and more precise measurement of differential expression overall. Furthermore, we establish the ability to target a second abundant miRNA and to multiplex the blocking of two miRNAs simultaneously. For small RNA sequencing, this technique could fill a similar role as do ribosomal or globin removal technologies in messenger RNA sequencing.

EnD-Seq and AppEnD: sequencing 3' ends to identify nontemplated tails and degradation intermediates

Existing methods for detecting RNA intermediates resulting from exonuclease degradation are low-throughput and laborious. In addition, mapping the 3' ends of RNA molecules to the genome after high-throughput sequencing is challenging, particularly if the 3' ends contain post-transcriptional modifications. To address these problems, we developed EnD-Seq, a high-throughput sequencing protocol that preserves the 3' end of RNA molecules, and AppEnD, a computational method for analyzing high-throughput sequencing data. Together these allow determination of the 3' ends of RNA molecules, including nontemplated additions. Applying EnD-Seq and AppEnD to histone mRNAs revealed that a significant fraction of cytoplasmic histone mRNAs end in one or two uridines, which have replaced the 1-2 nt at the 3' end of mature histone mRNA maintaining the length of the histone transcripts. Histone mRNAs in fly embryos and ovaries show the same pattern, but with different tail nucleotide compositions. We increase the sensitivity of EnD-Seq by using cDNA priming to specifically enrich low-abundance tails of known sequence composition allowing identification of degradation intermediates. In addition, we show the broad applicability of our computational approach by using AppEnD to gain insight into 3' additions from diverse types of sequencing data, including data from small capped RNA sequencing and some alternative polyadenylation protocols.

Elimination of ligation dependent artifacts in T4 RNA ligase to achieve high efficiency and low bias microRNA capture

Adapter ligation is a critical first step in many microRNA analysis methods including microarray, qPCR, and sequencing. Previous studies have shown that ligation bias can have dramatic effects on both the fidelity of expression profiles and reproducibility across samples. We have developed a method for high efficiency and low bias microRNA capture by 3′ adapter ligation using T4 RNA ligase that does not require pooled adapters. Using a panel of 20 microRNA, we investigated the effects of ligase type, PEG concentration, ligase amount, adapter concentration, incubation time, incubation temperature, and adapter design on capture efficiency and bias. Of these factors, high PEG% was found to be critical in suppressing ligation bias. We obtained high average capture efficiency and low CV across the 20 microRNA panel, both in idealized buffer conditions (86%±10%) and total RNA spiking conditions (64%±17%). We demonstrate that this method is reliable across microRNA species that previous studies have had difficulty capturing and that our adapter design performs significantly better than the common adapter designs. Further, we demonstrate that the optimization methodology must be specifically designed for minimizing bias in order to obtain the ideal reaction parameters.

Library construction for next-generation sequencing: overviews and challenges

High-throughput sequencing, also known as next-generation sequencing (NGS), has revolutionized genomic research. In recent years, NGS technology has steadily improved, with costs dropping and the number and range of sequencing applications increasing exponentially. Here, we examine the critical role of sequencing library quality and consider important challenges when preparing NGS libraries from DNA and RNA sources. Factors such as the quantity and physical characteristics of the RNA or DNA source material as well as the desired application (i.e., genome sequencing, targeted sequencing, RNA-seq, ChIP-seq, RIP-seq, and methylation) are addressed in the context of preparing high quality sequencing libraries. In addition, the current methods for preparing NGS libraries from single cells are also discussed.

A high-throughput microRNA expression profiling system

As small noncoding RNAs, microRNAs (miRNAs) regulate diverse biological functions, including physiological and pathological processes. The expression and deregulation of miRNA levels contain rich information with diagnostic and prognostic relevance and can reflect pharmacological responses. The increasing interest in miRNA-related research demands global miRNA expression profiling on large numbers of samples. We describe here a robust protocol that supports high-throughput sample labeling and detection on hundreds of samples simultaneously. This method employs 96-well-based miRNA capturing from total RNA samples and on-site biochemical reactions, coupled with bead-based detection in 96-well format for hundreds of miRNAs per sample. With low-cost, high-throughput, high detection specificity, and flexibility to profile both small and large numbers of samples, this protocol can be adapted in a wide range of laboratory settings.

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.

Quantification of microRNA expression with next-generation sequencing

Rapid advancement of next-generation sequencing technologies has made it possible to study expression profiles of microRNAs (miRNAs) comprehensively and efficiently. Multiplexing miRNA libraries by barcoding can significantly reduce sequencing cost per sample without compromising library quality. This unit provides a step-by-step protocol for isolating miRNAs and constructing multiplexed miRNA libraries. Also described is a custom computational pipeline for analyzing the multiplexed miRNA library sequencing reads generated by Illumina-based technology.

MicroRNAs in human cancer

Mature microRNAs (miRNAs) are single-stranded RNA molecules of 20-23-nucleotide (nt) length that control gene expression in many cellular processes. These molecules typically reduce the translation and stability of mRNAs, including those of genes that mediate processes in tumorigenesis, such as inflammation, cell cycle regulation, stress response, differentiation, apoptosis, and invasion. miRNA targeting is initiated through specific base-pairing interactions between the 5' end ("seed" region) of the miRNA and sites within coding and untranslated regions (UTRs) of mRNAs; target sites in the 3' UTR lead to more effective mRNA destabilization. Since miRNAs frequently target hundreds of mRNAs, miRNA regulatory pathways are complex. To provide a critical overview of miRNA dysregulation in cancer, we first discuss the methods currently available for studying the role of miRNAs in cancer and then review miRNA genomic organization, biogenesis, and mechanism of target recognition, examining how these processes are altered in tumorigenesis. Given the critical role miRNAs play in tumorigenesis processes and their disease specific expression, they hold potential as therapeutic targets and novel biomarkers.

Barcoded cDNA library preparation for small RNA profiling by next-generation sequencing

The characterization of post-transcriptional gene regulation by small regulatory (20-30 nt) RNAs, particularly miRNAs and piRNAs, has become a major focus of research in recent years. A prerequisite for characterizing small RNAs is their identification and quantification across different developmental stages, and in normal and disease tissues, as well as model cell lines. Here we present a step-by-step protocol for generating barcoded small RNA cDNA libraries compatible with Illumina HiSeq sequencing, thereby facilitating miRNA and other small RNA profiling of large sample collections.

A cost-effective method for Illumina small RNA-Seq library preparation using T4 RNA ligase 1 adenylated adapters

BACKGROUND

Deep sequencing is a powerful tool for novel small RNA discovery. Illumina small RNA sequencing library preparation requires a pre-adenylated 3' end adapter containing a 5',5'-adenyl pyrophosphoryl moiety. In the absence of ATP, this adapter can be ligated to the 3' hydroxyl group of small RNA, while RNA self-ligation and concatenation are repressed. Pre-adenylated adapters are one of the most essential and costly components required for library preparation, and few are commercially available.

RESULTS

We demonstrate that DNA oligo with 5' phosphate and 3' amine groups can be enzymatically adenylated by T4 RNA ligase 1 to generate customized pre-adenylated adapters. We have constructed and sequenced a small RNA library for tomato (Solanum lycopersicum) using the T4 RNA ligase 1 adenylated adapter.

CONCLUSION

We provide an efficient and low-cost method for small RNA sequencing library preparation, which takes two days to complete and costs around $20 per library. This protocol has been tested in several plant species for small RNA sequencing including sweet potato, pepper, watermelon, and cowpea, and could be readily applied to any RNA samples.

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

Eukaryotic regulatory small RNAs (sRNAs) play significant roles in many fundamental cellular processes. As such, they have emerged as useful biomarkers for diseases and cell differentiation states. sRNA-based biomarkers outperform traditional messenger RNA-based biomarkers by testing fewer targets with greater accuracy and providing earlier detection for disease states. Therefore, expression profiling of sRNAs is fundamentally important to further advance the understanding of biological processes, as well as diagnosis and treatment of diseases. High-throughput sequencing (HTS) is a powerful approach for both sRNA discovery and expression profiling. Here, we discuss the general considerations for sRNA-based HTS profiling methods from RNA preparation to sequencing library construction, with a focus on the causes of systematic error. By examining the enzymatic manipulation steps of sRNA expression profiling, this paper aims to demystify current HTS-based sRNA profiling approaches and to aid researchers in the informed design and interpretation of profiling experiments.

Multiple insert size paired-end sequencing for deconvolution of complex transcriptomes

Deep sequencing of transcriptomes allows quantitative and qualitative analysis of many RNA species in a sample, with parallel comparison of expression levels, splicing variants, natural antisense transcripts, RNA editing and transcriptional start and stop sites the ideal goal. By computational modeling, we show how libraries of multiple insert sizes combined with strand-specific, paired-end (SS-PE) sequencing can increase the information gained on alternative splicing, especially in higher eukaryotes. Despite the benefits of gaining SS-PE data with paired ends of varying distance, the standard Illumina protocol allows only non-strand-specific, paired-end sequencing with a single insert size. Here, we modify the Illumina RNA ligation protocol to allow SS-PE sequencing by using a custom pre-adenylated 3' adaptor. We generate parallel libraries with differing insert sizes to aid deconvolution of alternative splicing events and to characterize the extent and distribution of natural antisense transcription in C. elegans. Despite stringent requirements for detection of alternative splicing, our data increases the number of intron retention and exon skipping events annotated in the Wormbase genome annotations by 127% and 121%, respectively. We show that parallel libraries with a range of insert sizes increase transcriptomic information gained by sequencing and that by current established benchmarks our protocol gives competitive results with respect to library quality.

High-throughput multiplex sequencing of miRNA

Next-generation sequencing offers many advantages over other methods of microRNA (miRNA) expression profiling, such as sample throughput and the capability to discover novel miRNAs. As the sequencing depth of current sequencing platforms exceeds what is necessary to quantify miRNAs, multiplexing several samples in one sequencing run offers a significant cost advantage. Although previous studies have achieved this goal by adding bar codes to miRNA libraries at the ligation step, this was recently shown to introduce significant bias into the miRNA expression data. This bias can be avoided, however, by bar coding the miRNA libraries at the PCR step instead. Here, we describe a user-friendly PCR bar-coding method of preparing multiplexed microRNA libraries for Illumina-based sequencing. The method also prevents the production of adapter dimers and can be completed in one day.

Quantitative bias in Illumina TruSeq and a novel post amplification barcoding strategy for multiplexed DNA and small RNA deep sequencing

Here we demonstrate a method for unbiased multiplexed deep sequencing of RNA and DNA libraries using a novel, efficient and adaptable barcoding strategy called Post Amplification Ligation-Mediated (PALM). PALM barcoding is performed as the very last step of library preparation, eliminating a potential barcode-induced bias and allowing the flexibility to synthesize as many barcodes as needed. We sequenced PALM barcoded micro RNA (miRNA) and DNA reference samples and evaluated the quantitative barcode-induced bias in comparison to the same reference samples prepared using the Illumina TruSeq barcoding strategy. The Illumina TruSeq small RNA strategy introduces the barcode during the PCR step using differentially barcoded primers, while the TruSeq DNA strategy introduces the barcode before the PCR step by ligation of differentially barcoded adaptors. Results show virtually no bias between the differentially barcoded miRNA and DNA samples, both for the PALM and the TruSeq sample preparation methods. We also multiplexed miRNA reference samples using a pre-PCR barcode ligation. This barcoding strategy results in significant bias.

Deep sequencing of microRNA precursors reveals extensive 3' end modification

MicroRNAs (miRNAs) are small, noncoding RNAs that post-transcriptionally regulate gene expression. An emerging mechanism to control miRNA production is the addition of an oligo-uridine tail to the 3' end of the precursor miRNA. This has been demonstrated for the Let-7 family of miRNAs in embryonic cells. Additionally, nontemplated nucleotides have been found on mature miRNA species, though in most cases it is not known if nucleotide addition occurs at the precursor step or at the mature miRNA. To examine the diversity of nucleotide addition we have developed a high-throughput sequencing method specific for miRNA precursors. Here we report that nontemplated addition is a widespread phenomenon occurring in many miRNA families. As previously reported, Let-7 family members are oligo-uridylated in embryonic cells in a Lin28-dependent manner. However, we find that the fraction of uridylated precursors increases with differentiation, independent of Lin28, and is highest in adult mouse tissues, exceeding 30% of all sequence reads for some Let-7 family members. A similar fraction of sequence reads are modified for many other miRNA families. Mono-uridylation is most common, with cytidine and adenosine modification less frequent but occurring above the expected error rate for Illumina sequencing. Nucleotide addition in cell lines is associated with 3' end degradation, in contrast to adult tissues, where modification occurs predominantly on full-length precursors. This work provides an unprecedented view of the complexity of 3' modification and trimming of miRNA precursors.

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

Sequencing of small RNA cDNA libraries is an important tool for the discovery of new RNAs and the analysis of their mutational status as well as expression changes across samples. It requires multiple enzyme-catalyzed steps, including sequential oligonucleotide adapter ligations to the 3' and 5' ends of the small RNAs, reverse transcription (RT), and PCR. We assessed biases in representation of miRNAs relative to their input concentration, using a pool of 770 synthetic miRNAs and 45 calibrator oligoribonucleotides, and tested the influence of Rnl1 and two variants of Rnl2, Rnl2(1-249) and Rnl2(1-249)K227Q, for 3'-adapter ligation. The use of the Rnl2 variants for adapter ligations yielded substantially fewer side products compared with Rnl1; however, the benefits of using Rnl2 remained largely obscured by additional biases in the 5'-adapter ligation step; RT and PCR steps did not have a significant impact on read frequencies. Intramolecular secondary structures of miRNA and/or miRNA/3'-adapter products contributed to these biases, which were highly reproducible under defined experimental conditions. We used the synthetic miRNA cocktail to derive correction factors for approximation of the absolute levels of individual miRNAs in biological samples. Finally, we evaluated the influence of 5'-terminal 5-nt barcode extensions for a set of 20 barcoded 3' adapters and observed similar biases in miRNA read distribution, thereby enabling cost-saving multiplex analysis for large-scale miRNA profiling.

Barcoding bias in high-throughput multiplex sequencing of miRNA

Second-generation sequencing is gradually becoming the method of choice for miRNA detection and expression profiling. Given the relatively small number of miRNAs and improvements in DNA sequencing technology, studying miRNA expression profiles of multiple samples in a single flow cell lane becomes feasible. Multiplexing strategies require marking each miRNA library with a DNA barcode. Here we report that barcodes introduced through adapter ligation confer significant bias on miRNA expression profiles. This bias is much higher than the expected Poisson noise and masks significant expression differences between miRNA libraries. This bias can be eliminated by adding barcodes during PCR amplification of libraries. The accuracy of miRNA expression measurement in multiplexed experiments becomes a function of sample number.

Simple and efficient synthesis of 5' pre-adenylated DNA using thermostable RNA ligase

We report a simple method of enzymatic synthesis of pre-adenylated DNA linkers/adapters for next-generation sequencing using thermostable RNA ligase from Methanobacterium thermoautotrophicum (MthRnl). Using RNA ligase for the reaction instead of the existing chemical or T4 DNA ligase-based methods allows quantitative conversion of 5′-phosphorylated single-stranded DNA (ssDNA) to the adenylated form. The MthRnl adenylation reaction is specific for ATP and either ssDNA or RNA. In the presence of Mg⁺², the reaction has a pH optimum of 6.0–6.5. Unlike reactions that use T4 DNA ligase, this protocol does not require synthesis of a template strand for adenylation. The high yield of the reaction simplifies isolation and purification of the adenylated product. Conducting the adenylation reaction at the elevated temperature (65°C) reduces structural constraints, while increased ATP concentrations allow quantitative adenylation of DNA with a 3′-unprotected end.

Reduction of non-insert sequence reads by dimer eliminator LNA oligonucleotide for small RNA deep sequencing

Here we describe a method for constructing small RNA libraries for high-throughput sequencing in which we have made a significant improvement to commonly available standard protocols. We added a locked nucleic acid (LNA) oligonucleotide--named dimer eliminator--that is complementary to the adapter-dimer ligation products during the reverse transcription reaction. It reduces adapter-dimers, which often contaminate standard libraries and increase the number of non-insert sequence reads. This simple technology can be used for simultaneous multiplex sequencing of various barcoded samples as well as nonbarcoded small RNA library sequencing. In this study we also evaluated the reproducibility and quantitative design of the eight barcoded tags by comparing the Pearson's correlation values in the expression analysis between each barcoded sample. This method improves the sequencing yield and efficiency, while simplifying library construction, and makes it easier to perform large-scale small RNA analysis under multiple conditions with next-generation sequencers.

Construction of small RNA cDNA libraries for deep sequencing

Since the phenomenon of small RNA-mediated gene silencing was first described over 15 years ago (Lee et al. Cell 75:843-854, 1993; Wightman et al. Cell 75:855-862, 1993), it has become evident that a variety of endogenous small RNAs play an important role in establishing and maintaining cell lineages. MicroRNAs (miRNAs), in particular, have been shown to exert regulatory control over the development and function of the many specialized cells that comprise the mammalian immune system (Baltimore et al. Nat Immunol 9:839-845, 2008; Kanellopoulous and Monticelli Semin Cancer Biol 18:79-88, 2008; Xiao and Rajewsky Cell 136:26-36, 2009). The advent of next generation sequencers provides an important tool for profiling the small RNA transcriptome of many diverse cell types. Compared to traditional Sanger sequencing, next generation sequencing machines can process millions of sequence reads in parallel, generating megabases of data within just a few days. The generation of small RNA libraries for sequencing is relatively straightforward and involves the ligation of platform-specific adapter sequences to small RNAs, followed by reverse transcription of the ligated species and PCR amplification. While other hybridization-based techniques are available for profiling well-characterized small RNAs, high-throughput sequencing remains the most powerful method for discovering novel small RNAs and posttranscriptional editing.

miRNAs in human cancer

Mature microRNAs (miRNAs) are single-stranded RNA molecules of 20-23 nucleotide (nt) length that control gene expression in many cellular processes. These molecules typically reduce the stability of mRNAs, including those of genes that mediate processes in tumorigenesis, such as inflammation, cell cycle regulation, stress response, differentiation, apoptosis and invasion. miRNA targeting is mostly achieved through specific base-pairing interactions between the 5' end ('seed' region) of the miRNA and sites within coding and untranslated regions (UTRs) of mRNAs; target sites in the 3' UTR lead to more effective mRNA destabilization. Since miRNAs frequently target hundreds of mRNAs, miRNA regulatory pathways are complex. To provide a critical overview of miRNA dysregulation in cancer, we first discuss the methods currently available for studying the role of miRNAs in cancer and then review miRNA genomic organization, biogenesis and mechanism of target recognition, examining how these processes are altered in tumorigenesis. Given the critical role miRNAs play in tumorigenesis processes and their disease-specific expression, they hold potential as therapeutic targets and novel biomarkers.

Recognition of 5' triphosphate by RIG-I helicase requires short blunt double-stranded RNA as contained in panhandle of negative-strand virus

Antiviral immunity is triggered by immunorecognition of viral nucleic acids. The cytosolic helicase RIG-I is a key sensor of viral infections and is activated by RNA containing a triphosphate at the 5' end. The exact structure of RNA activating RIG-I remains controversial. Here, we established a chemical approach for 5' triphosphate oligoribonucleotide synthesis and found that synthetic single-stranded 5' triphosphate oligoribonucleotides were unable to bind and activate RIG-I. Conversely, the addition of the synthetic complementary strand resulted in optimal binding and activation of RIG-I. Short double-strand conformation with base pairing of the nucleoside carrying the 5' triphosphate was required. RIG-I activation was impaired by a 3' overhang at the 5' triphosphate end. These results define the structure of RNA for full RIG-I activation and explain how RIG-I detects negative-strand RNA viruses that lack long double-stranded RNA but do contain blunt short double-stranded 5' triphosphate RNA in the panhandle region of their single-stranded genome.

Suppression subtractive hybridisation allows selective sampling of metagenomic subsets of interest

Metagenomic studies bypass the requirement of a pure culture for analysis, focusing instead on the genetic information present in a given sample. Metagenomics have been applied to various studies, with objectives ranging from genome reconstruction, gene prospecting and ecology. However, the use of metagenomics in comparative studies has been constrained by sequencing costs and computational limitations. Efforts are underway to improve current sequencing methods and reduce the expense involved. We suggest an alternative approach - pretreatment of the sample of interest to enrich for desired subsets prior to deep sequencing. In this study, we tested the use of suppression subtractive hybridisation (SSH) for in vitro separation of metagenomic samples based on temporal variance. Faecal samples were taken from pigs at different timepoints and extracted DNA was whole genome-amplified using multiple displacement amplification (MDA). A sample collected at 31 days of age was designated the tester while a 24 day sample was denoted the driver. Following hybridisation and subtraction, tester-specific sequences are expected to be enriched in the final sample while driver-specific and common sequences are removed. Using denaturing gel gradient electrophoresis (DGGE), we found that driver-specific bands were completely removed from all final profiles while an average of 70% of common bands were successfully subtracted. Final profiles retained an average of 70% of tester-specific sequences and new sequences contributed an average of 36% of the band mobilities found in the final profiles. Tester-unique sequences were inferred to make up 78% of the final profile after SSH. We expect that using subtractive hybridisation for separation of metagenomic samples into desired subsets will provide a more effective and targeted approach to comparative studies.