Ligation bias is a major contributor to nonstoichiometric abundances of secondary siRNAs and impacts analyses of microRNAs

ABSTRACTPlant small RNAs are a diverse and complex set of molecules, ranging in length from 21 to 24 nt, involved in a wide range of essential biological processes. High-throughput sequencing is used for the discovery and quantification of small RNAs. However, several biases can occur during the preparation of small RNA libraries, especially using low input RNA. We used two stages of maize anthers to evaluate the performance of seven commercially-available methods for small RNA library construction, using different RNA input amounts. We show that when working with plant material, library construction methods have differing capabilities to capture small RNAs, and that different library construction methods provide better results when applied to the detection of microRNAs, phasiRNAs, or tRNA-derived fragment. We also observed that ligation bias occurs at both ends of miRNAs and phasiRNAs, suggesting that the biased compositions observed in small RNA populations, including nonstoichiometric levels of phasiRNAs within a locus, may reflect a combination of biological and technical influences.

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.

Comparative Analysis of Biochemical Biases by Ligation- and Template-Switch-Based Small RNA Library Preparation Protocols

BACKGROUND Small RNAs are key players in the regulation of gene expression and differentiation. However, many different classes of small RNAs (sRNAs) have been described with distinct biogenesis pathways and, as a result, with different biochemical properties. To analyze sRNAs by deep sequencing, complementary DNA synthesis requires manipulation of the RNA molecule itself. Thus, enzymatic activities during library preparation bias the library content owing to biochemical criteria. METHODS We compared 4 different manipulations of RNA for library preparation: (a) a ligation-based procedure allowing only 5′-mono-phosphorylated RNA to enter the library, (b) a ligation-based procedure allowing additional 5′-triphosphates and Cap structures, (c) a ligation-independent, template-switch-based library preparation, and (d) a template-switch-based library preparation allowing 3′-phosphorylated RNAs to enter the library. RESULTS Our data show large differences between ligation-dependent and ligation-independent libraries in terms of their preference for individual sRNA classes such as microRNAs (miRNAs), Piwi-interacting RNAs (piRNAs), and transfer RNA fragments. Moreover, the miRNA composition is different between both procedures, and more microRNA isoforms (isomiRs) can be identified after pyrophosphatase treatment. piRNAs are enriched in template-switch libraries, and this procedure apparently includes more different RNA species. CONCLUSIONS Our data indicate that miRNAomics from both methods will hardly be comparable. Ligation-based libraries enrich for canonical miRNAs, which thus may be suitable methods for miRNAomics. Template-switch libraries contain increased numbers and different compositions of fragments and long RNAs. Following different interests for other small RNA species, ligation-independent libraries appear to show a more realistic sRNA landscape with lower bias against biochemical modifications.

An improved method of constructing degradome library suitable for sequencing using Illumina platform

Background Post-transcriptional gene regulation is one of the critical layers of overall gene expression programs and microRNAs (miRNAs) play an indispensable role in this process by guiding cleavage on the messenger RNA targets. The transcriptome-wide cleavages on the target transcripts can be identified by analyzing the degradome or PARE or GMUCT libraries. However, high-throughput sequencing of PARE or degradome libraries using Illumina platform, a widely used platform, is not so straightforward. Moreover, the currently used degradome or PARE methods utilize MmeI restriction site in the 5′ RNA adapter and the resulting fragments are only 20-nt long, which often poses difficulty in distinguishing between the members of the same target gene family or distinguishing miRNA biogenesis intermediates from the primary miRNA transcripts belonging to the same miRNA family. Consequently, developing a method which can generate longer fragments from the PARE or degradome libraries which can also be sequenced easily using Illumina platform is ideal. Results In this protocol, 3′ end of the 5′RNA adaptor of TruSeq small RNA library is modified by introducing EcoP15I recognition site. Correspondingly, the double-strand DNA (dsDNA) adaptor sequence is also modified to suit with the ends generated by the restriction enzyme EcoP15I. These modifications allow amplification of the degradome library by primer pairs used for small RNA library preparation, thus amenable for sequencing using Illumina platform, like small RNA library. Conclusions Degradome library generated using this improved protocol can be sequenced easily using Illumina platform, and the resulting tag length is ~ 27-nt, which is longer than the MmeI generated fragment (20-nt) that can facilitate better accuracy in validating target transcripts belonging to the same gene family or distinguishing miRNA biogenesis intermediates of the same miRNA family. Furthermore, this improved method allows pooling and sequencing degradome libraries and small RNA libraries simultaneously using Illumina platform.

Non-Coding RNA Sequencing of Equine Endometrium During Maternal Recognition of Pregnancy

Maternal recognition of pregnancy (MRP) in the mare is not well defined. In a non-pregnant mare, prostaglandin F2α (PGF) is released on day 14 post-ovulation (PO) to cause luteal regression, resulting in loss of progesterone production. Equine MRP occurs prior to day 14 to halt PGF production. Studies have failed to identify a gene candidate for MRP, so attention has turned to small, non-coding RNAs. The objective of this study was to evaluate small RNA (<200 nucleotides) content in endometrium during MRP. Mares were used in a cross-over design with each having a pregnant and non-mated cycle. Each mare was randomly assigned to collection day 11 or 13 PO (n = 3/day) and endometrial biopsies were obtained. Total RNA was isolated and sequencing libraries were prepared using a small RNA library preparation kit and sequenced on a HiSeq 2000. EquCab3 was used as the reference genome and DESeq2 was used for statistical analysis. On day 11, 419 ncRNAs, representing miRNA, snRNA, snoRNA, scaRNA, and vaultRNA, were different between pregnancy statuses, but none on day 13. Equine endometrial ncRNAs with unknown structure and function were also identified. This study is the first to describe ncRNA transcriptome in equine endometrium. Identifying targets of these ncRNAs could lead to determining MRP.

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.

Comprehensive analysis of small RNA profile by massive parallel sequencing in HTLV-1 asymptomatic subjects with monoclonal and polyclonal rearrangement of the T-cell antigen receptor γ-chain

IntroductionIn this study, we used a massive parallel sequencing technology to investigate the cellular small RNA (sRNA) operating in peripheral blood mononuclear cells (PBMCs) of the Human T-lymphotropic virus type I (HTLV-I) infected asymptomatic subjects with a monoclonal and polyclonal rearrangement of the T-cell antigen receptor γ-chain.Materials and MethodsBlood samples from 15 HTLV-1 asymptomatic carriers who were tested for clonalTCR-γ gene (seven and eight subjects presented monoclonal and polyclonal expansion of HTLV-1 infected cells, respectively), and were submitted to Illumina for small RNA library construction. sRNA libraries were prepared from cryoperserved PBMCs using TrueSeq Small RNA Library Preparation Kit (Illumina). The sRNA-Seq reads were aligned, annotated, and profiled by different bioinformatics tools.ResultsThrough bioinformatics analysis, we identified a total of 494 known sRNAs and 120 putative novel sRNAs. Twenty-two known and 15 novel sRNA showed a different expression (>2-fold) between the asymptomatic monoclonal (ASM) and asymptomatic polyclonal carriers (ASP). The hsa-mir-196a-5p was the most abundantly upregulated micro RNA (miRNA) and the hsa-mir-133a followed by hsa-mir-509-3p were significantly downregulated miRNAs with more than a three-fold difference in the ASM than ASP group. The target genes predicted to be regulated by the differentially expressed miRNAs play essential roles in diverse biological processes including cell proliferation, differentiation, and/or apoptosis.DiscussionOur results provide an opportunity for a further understanding of sRNA regulation and function in HTLV-1 infected subjects with monoclonality evidence.