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

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.

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.

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

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

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.

QsRNA-seq: A protocol for generating libraries for high-throughput sequencing of small RNAs

Small RNAs (sRNAs) are 20-30 nt long non-coding RNA molecules that regulate essentially all cellular processes. Besides being an intensively studied topic in academic research, sRNAs also hold a promise as clinical biomarkers. While the need for expressional profiling of sRNAs is growing, preparation of sRNA libraries for high-throughput sequencing (HTS) remains technically challenging, due to their small size. The common PAGE-based protocol is time-consuming and inefficient due to material loss, while gel-free protocols generate libraries of insufficient quality. To overcome these shortcomings, we modified the conditions of size-selection by Solid Phase Reversible Immobilization (SPRI) in a way that allows separation of nucleic acids shorter than 100 nt and differing in length by only 20 nt. Implementing the method for preparation of small RNA libraries for HTS resulted in QsRNA-seq, a gel-free, fast and easy-to-perform protocol, amenable to automation, generating very clean libraries that result in high-depth expression data. The protocol also utilizes Unique Molecular Identifiers (UMI) for reduction of library preparation biases and to quantify expression levels. QsRNA-seq provides an excellent solution to the growing needs for small RNA expression profiling for research clinical use.

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.

Characterization of piRNAs across postnatal development in mouse brain

PIWI-interacting RNAs (piRNAs) are responsible for maintaining the genome stability by silencing retrotransposons in germline tissues– where piRNAs were first discovered and thought to be restricted. Recently, novel functions were reported for piRNAs in germline and somatic cells. Using deep sequencing of small RNAs and CAGE of postnatal development of mouse brain, we identified piRNAs only in adult mouse brain. These piRNAs have similar sequence length as those of MILI-bound piRNAs. In addition, we predicted novel candidate regulators and putative targets of adult brain piRNAs.

Addressing Bias in Small RNA Library Preparation for Sequencing: A New Protocol Recovers MicroRNAs that Evade Capture by Current Methods

Recent advances in sequencing technology have helped unveil the unexpected complexity and diversity of small RNAs. A critical step in small RNA library preparation for sequencing is the ligation of adapter sequences to both the 5' and 3' ends of small RNAs. Studies have shown that adapter ligation introduces a significant but widely unappreciated bias in the results of high-throughput small RNA sequencing. We show that due to this bias the two widely used Illumina library preparation protocols produce strikingly different microRNA (miRNA) expression profiles in the same batch of cells. There are 102 highly expressed miRNAs that are >5-fold differentially detected and some miRNAs, such as miR-24-3p, are over 30-fold differentially detected. While some level of bias in library preparation is not surprising, the apparent massive differential bias between these two widely used adapter sets is not well appreciated. In an attempt to mitigate this bias, the new Bioo Scientific NEXTflex V2 protocol utilizes a pool of adapters with random nucleotides at the ligation boundary. We show that this protocol is able to detect robustly several miRNAs that evade capture by the Illumina-based methods. While these analyses do not indicate a definitive gold standard for small RNA library preparation, the results of the NEXTflex protocol do correlate best with RT-qPCR. As increasingly more laboratories seek to study small RNAs, researchers should be aware of the extent to which the results may differ with different protocols, and should make an informed decision about the protocol that best fits their study.

Chromogen detection of microRNA in frozen clinical tissue samples using LNA™ probe technology

Specific chromogen- and fluorescence-based detection of microRNA by in situ hybridization (ISH) in formalin-fixed and paraffin-embedded (FFPE) tissue sections has been facilitated by locked nucleic acid (LNA)-based probe technology and can be performed within a single working day. In the current method paper, we present a similar simple 1-day ISH method developed for cryostat sections obtained from clinical cryo-embedded tissue samples. The presented chromogen-based ISH method does not involve proteolytic pretreatment, which is mandatory for FFPE sections, but still retains a sensitivity level similar to that obtained in FFPE sections. The LNA-based ISH method is not only applicable in situations where only access to cryo-embedded material is possible, but it also has a potential use if combining microRNA ISH with immunohistochemistry in double fluorescence staining with antibodies not being compatible with proteolytic predigestion.

An improved protocol for small RNA library construction using High Definition adapters

Next generation sequencing of small RNA (sRNA) libraries is widely used for studying sRNAs in various biological systems. However, cDNA libraries of sRNAs are biased for molecules that are ligated to adapters more or less efficiently than other molecules. One approach to reduce this ligation bias is to use a pool of adapters instead of a single adapter sequence, which allows many sRNAs to be ligated efficiently. We previously developed High Definition (HD) adapters for the Illumina sequencing platform, which contain degenerate nucleotides at the ligating ends of the adapters. However, the current commercial kits produced a large amount of 5’ adapter – 3’ adapter ligation product without the cDNA insert when HD adapters were used to replace the kit adapters. Here, we report a protocol to generate sRNA libraries using HD adapters with greatly reduced proportion of adapter-adapter products due to the degradation of nonligated 3’ adapters. The libraries can be completed within two days and can be used for various biological and clinical samples. As examples for using this protocol, we constructed sRNA libraries using total RNA extracted from cultured mammalian cells and plant leaf tissue.

Potential small guide RNAs for tRNase ZL from human plasma, peripheral blood mononuclear cells, and cultured cell lines

Several pieces of evidence suggest that small RNA degradation products together with tRNase ZL appear to form another layer of the whole gene regulatory network. The degraded RNA such as a 5'-half-tRNA and an rRNA fragment function as small guide RNA (sgRNA) to guide the enzyme to target RNA. We were curious whether there exist RNAs in plasma that can function as sgRNAs for tRNase ZL, whether these RNAs are working as signaling molecules between cells to fulfill physiological roles, and whether there are any differences in plasma sgRNA species and levels between normal and pathological conditions. Here, we analyzed small plasma RNAs from three healthy persons and three multiple myeloma patients for potential sgRNAs by deep sequencing. We also examined small RNAs from peripheral blood mononuclear cells (PBMC) of three healthy persons and three myeloma patients and from various cultured human cell lines for sgRNAs. We found that read-number distribution patterns of plasma and PBMC RNAs differ between persons in the range of 5-40 nt and that there are many RNA species that exist significantly more or less abundantly in the plasma or PBMC of the myeloma patients than those of the healthy persons. Furthermore, we found that there are many potential sgRNAs in the 5-40-nt RNAs and that, among them, a 31-nt RNA fragment derived from 94-nt Y4-RNA, which can function as a 5'-half-tRNA-type sgRNA, is overwhelmingly abundant in the plasma of 2/3 of the examinees. These observations suggest that the gene regulatory network via tRNase ZL and sgRNA may be extended intercellularly.

Methods of quantifying microRNAs for hypoxia research: classic and next generation

SIGNIFICANCE

Recent evidence is uncovering the ever-increasing importance of microribonucleic acids (miRNAs) in the hypoxia response. In order to investigate the important roles that these small RNAs play, methods of quantification whether using classic single-gene methods or genome-wide technologies are necessary to obtain a global picture of the differential expression of miRNAs in hypoxia and their interplay with protein coding genes.

RECENT ADVANCES

Building on the groundwork of classic quantitative polymerase chain reaction (qPCR) and microarrays, the advent of next-generation sequencing technology has revolutionized how small RNAs can be detected and quantified on a genome-wide scale and without a priori knowledge of the small RNA sequence. This method delivers accurate and comprehensive data on the expression and sequence of all expressed small RNAs, and the data can be further combined with other sequencing modalities to better understand miRNAs via integrated genomic analyses.

CRITICAL ISSUES

Advancing technology has increased the need for better methods of sample and library preparation and for bioinformatics tools. Speed, cost, sample input, and analysis expertise remain the mainstay critical issues of small RNA sequencing.

FUTURE DIRECTIONS

Future hypoxia research will benefit from the application of genome-wide sequencing technologies. Analyses that combine genomic, transcriptomic, chromosome conformation, DNA/RNA-protein binding, and proteomics will help greatly advance hypoxia miRNA research.

Multiplicity of 5' cap structures present on short RNAs

Most RNA molecules are co- or post-transcriptionally modified to alter their chemical and functional properties to assist in their ultimate biological function. Among these modifications, the addition of 5′ cap structure has been found to regulate turnover and localization. Here we report a study of the cap structure of human short (<200 nt) RNAs (sRNAs), using sequencing of cDNA libraries prepared by enzymatic pretreatment of the sRNAs with cap sensitive-specificity, thin layer chromatographic (TLC) analyses of isolated cap structures and mass spectrometric analyses for validation of TLC analyses. Processed versions of snoRNAs and tRNAs sequences of less than 50 nt were observed in capped sRNA libraries, indicating additional processing and recapping of these annotated sRNAs biotypes. We report for the first time 2,7 dimethylguanosine in human sRNAs cap structures and surprisingly we find multiple type 0 cap structures (mGpppC, 7mGpppG, GpppG, GpppA, and 7mGpppA) in RNA length fractions shorter than 50 nt. Finally, we find the presence of additional uncharacterized cap structures that wait determination by the creation of needed reference compounds to be used in TLC analyses. These studies suggest the existence of novel biochemical pathways leading to the processing of primary and sRNAs and the modifications of their RNA 5′ ends with a spectrum of chemical modifications.

Comparison of CAGE and RNA-seq transcriptome profiling using clonally amplified and single-molecule next-generation sequencing

CAGE (cap analysis gene expression) and RNA-seq are two major technologies used to identify transcript abundances as well as structures. They measure expression by sequencing from either the 5' end of capped molecules (CAGE) or tags randomly distributed along the length of a transcript (RNA-seq). Library protocols for clonally amplified (Illumina, SOLiD, 454 Life Sciences [Roche], Ion Torrent), second-generation sequencing platforms typically employ PCR preamplification prior to clonal amplification, while third-generation, single-molecule sequencers can sequence unamplified libraries. Although these transcriptome profiling platforms have been demonstrated to be individually reproducible, no systematic comparison has been carried out between them. Here we compare CAGE, using both second- and third-generation sequencers, and RNA-seq, using a second-generation sequencer based on a panel of RNA mixtures from two human cell lines to examine power in the discrimination of biological states, detection of differentially expressed genes, linearity of measurements, and quantification reproducibility. We found that the quantified levels of gene expression are largely comparable across platforms and conclude that CAGE and RNA-seq are complementary technologies that can be used to improve incomplete gene models. We also found systematic bias in the second- and third-generation platforms, which is likely due to steps such as linker ligation, cleavage by restriction enzymes, and PCR amplification. This study provides a perspective on the performance of these platforms, which will be a baseline in the design of further experiments to tackle complex transcriptomes uncovered in a wide range of cell types.

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.

Suppression of artifacts and barcode bias in high-throughput transcriptome analyses utilizing template switching

Template switching (TS) has been an inherent mechanism of reverse transcriptase, which has been exploited in several transcriptome analysis methods, such as CAGE, RNA-Seq and short RNA sequencing. TS is an attractive option, given the simplicity of the protocol, which does not require an adaptor mediated step and thus minimizes sample loss. As such, it has been used in several studies that deal with limited amounts of RNA, such as in single cell studies. Additionally, TS has also been used to introduce DNA barcodes or indexes into different samples, cells or molecules. This labeling allows one to pool several samples into one sequencing flow cell, increasing the data throughput of sequencing and takes advantage of the increasing throughput of current sequences. Here, we report TS artifacts that form owing to a process called strand invasion. Due to the way in which barcodes/indexes are introduced by TS, strand invasion becomes more problematic by introducing unsystematic biases. We describe a strategy that eliminates these artifacts in silico and propose an experimental solution that suppresses biases from TS.

Site-specific DICER and DROSHA RNA products control the DNA-damage response

Non-coding RNAs (ncRNAs) are involved in an increasing number of cellular events¹. Some ncRNAs are processed by DICER and DROSHA ribonucleases to give rise to small double-stranded RNAs involved in RNA interference (RNAi)². The DNA-damage response (DDR) is a signaling pathway that originates from the DNA lesion and arrests cell proliferation³. So far, DICER or DROSHA RNA products have not been reported to control DDR activation. Here we show that DICER and DROSHA, but not downstream elements of the RNAi pathway, are necessary to activate DDR upon oncogene-induced genotoxic stress and exogenous DNA damage, as studied also by DDR foci formation in mammalian cells and zebrafish and by checkpoint assays. DDR foci are sensitive to RNase A treatment, and DICER- and DROSHA-dependent RNA products are required to restore DDR foci in treated cells. Through RNA deep sequencing and studies of DDR activation at an inducible unique DNA double-strand break (DSB), we demonstrate that DDR foci formation requires site-specific DICER- and DROSHA-dependent small RNAs, named DDRNAs, which act in a MRE11-RAD50-NBS1 (MRN) complex-dependent manner. Chemically synthesized or in vitro-generated by DICER cleavage, DDRNAs are sufficient to restore DDR in RNase A-treated cells, also in the absence of other cellular RNAs. Our results describe an unanticipated direct role of a novel class of ncRNAs in the control of DDR activation at sites of DNA damage.

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.

Chromatin-associated RNA interference components contribute to transcriptional regulation in Drosophila

RNAi pathways have evolved as important modulators of gene expression that act in the cytoplasm by degrading RNA target molecules via the activity of short (21-30nt) RNAs^1-6 RNAi components have been reported to play a role in the nucleus as they are involved in epigenetic regulation and heterochromatin formation^7-10. However, although RNAi-mediated post-transcriptional silencing (PTGS) is well documented, mechanisms of RNAi-mediated transcriptional gene silencing (TGS) and in particular the role of RNAi components in chromatin, especially in higher eukaryotes, are still elusive. Here we show that key RNAi components Dicer-2 (Dcr2) and and Argonaute-2 (AGO2) AGO2 associate with chromatin, with strong preference for euchromatic, transcriptionally active loci and interact with core transcription machinery. Notably Dcr2 and AGO2 loss of function show that transcriptional defects are accompanied by perturbation of Pol II positioning on promoters. Further, both Dcr2 and Ago2 null mutations as well as missense mutations compromising the RNAi activity impair global Pol II dynamics upon heat shock. Finally, AGO2 RIP-seq experiments reveal that, AGO2 is strongly enriched in small-RNAs encompassing promoter as well as other parts of heat shock and other gene loci on both sense and antisense, with a strong bias for antisense, particularly after heat shock. Taken together our results reveal a new scenario in which Dcr2 and AGO2 are globally associated with transcriptionally active loci and may play a pivotal role in shaping the transcriptome by controlling RNA Pol II processivity.

pre-miRNA profiles obtained through application of locked nucleic acids and deep sequencing reveals complex 5'/3' arm variation including concomitant cleavage and polyuridylation patterns

Recent research hints at an underappreciated complexity in pre-miRNA processing and regulation. Global profiling of pre-miRNA and its potential to increase understanding of the pre-miRNA landscape is impeded by overlap with highly expressed classes of other non coding (nc) RNA. Here, we present a data set excluding these RNA before sequencing through locked nucleic acids (LNA), greatly increasing pre-miRNA sequence counts with no discernable effect on pre-miRNA or mature miRNA sequencing. Analysis of profiles generated in total, nuclear and cytoplasmic cell fractions reveals that pre-miRNAs are subject to a wide range of regulatory processes involving loci-specific 3′- and 5′-end variation entailing complex cleavage patterns with co-occurring polyuridylation. Additionally, examination of nuclear-enriched flanking sequences of pre-miRNA, particularly those derived from polycistronic miRNA transcripts, provides insight into miRNA and miRNA-offset (moRNA) production, specifically identifying novel classes of RNA potentially functioning as moRNA precursors. Our findings point to particularly intricate regulation of the let-7 family in many ways reminiscent of DICER1-independent, pre-mir-451-like processing, introduce novel and unify known forms of pre-miRNA regulation and processing, and shed new light on overlooked products of miRNA processing pathways.

miRNA profiling of bilateral rat hippocampal CA3 by deep sequencing

MicroRNAs (miRNAs) have been demonstrated to be potent post-trascriptional modulators of protein expression. miRNA expression was profiled in the left and right dorsal hippocampal CA3 of mature rats by high-throughput deep sequencing. Among the sequenced and cross-mapped small RNAs, 88% belonged to the miRNAs annotated in the miRBase 15 database. Nearly half of the small RNAs belonged to the let-7 family miRNA. Seven percent of the sequenced small RNAs were not annotated in miRBase 15. Bioinformatic analysis of the unannotated small RNA sequences suggested seventeen novel miRNA candidates with relatively high expression levels (>100 tags per million). The left:right expression ratios were similar for all highly expressed miRNAs with less than 10% differences. These results provide a basic idea of the relative expression strengths of known and unknown miRNAs in the dorsal hippocampal CA3.

Whole transcriptome analysis: what are we still missing?

New technologies such as tag-based sequencing and tiling arrays have provided unique insights into the transcriptional output of cells. Many new RNA classes have been uncovered in the past decade, despite limitations in current technologies. Even as the repertoire of known functional elements of the transcriptome increases and contemporary technologies become mainstream, inadequacies in conventional protocols for library preparation, sequencing and mapping continue to hamper revelation of the entire transcriptome of cells. In this article, we review current protocols and outline their deficiencies. We also provide our view on what we may be overlooking in the transcriptome, despite exhaustive investigations, and indicate future areas of technological development and research.