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.     

Sorting of small RNAs into Arabidopsis argonaute complexes is directed by the 5' terminal nucleotide

Argonaute (AGO) proteins recruit small RNAs to form the core of RNAi effector complexes. Arabidopsis encodes ten AGO proteins and a large network of small RNAs. How these small RNAs are sorted into specific AGO complexes remains largely unknown. We have cataloged small RNAs resident in four AGO complexes. We found that AGO2 and AGO4 preferentially recruit small RNAs with a 5' terminal adenosine, whereas AGO1 harbors microRNAs (miRNAs) that favor a 5' terminal uridine. AGO5 predominantly binds small RNAs that initiate with cytosine. Changing the 5' terminal nucleotide of an miRNA predictably redirected it into a different AGO complex and alters its biological activity. These results reveal a role for small RNA sequences in assorting among AGO complexes. This suggests that specialization of AGO complexes might involve remodeling the 5' end-binding pocket to accept certain small RNA sequences, perhaps explaining the evolutionary drive for miRNAs to initiate with uridine.     

Ebbie: automated analysis and storage of small RNA cloning data using a dynamic web server

Background  

DNA sequencing is used ubiquitously: from deciphering genomes[1] to determining the primary sequence of small RNAs (smRNAs) [2–5]. The cloning of smRNAs is currently the most conventional method to determine the actual sequence of these important regulators of gene expression. Typical smRNA cloning projects involve the sequencing of hundreds to thousands of smRNA clones that are delimited at their 5' and 3' ends by fixed sequence regions. These primers result from the biochemical protocol used to isolate and convert the smRNA into clonable PCR products. Recently we completed a smRNA cloning project involving tobacco plants, where analysis was required for ~700 smRNA sequences[6]. Finding no easily accessible research tool to enter and analyze smRNA sequences we developed Ebbie to assist us with our study.     

Bombyx small RNAs: genomic defense system against transposons in the silkworm, Bombyx mori

Selfish genetic elements called transposons can insert themselves at new locations in host genomes to modify gene structure and alter gene expression. Expansion of transposons can occur when novel transposition events are transmitted to subsequent generations after germline hopping. Therefore, organisms seem likely to have evolved defense mechanisms to silence transposons in the germline. Recently, small RNAs interacting with Piwi proteins (piwi-interacting RNAs: piRNAs) have been demonstrated to be involved in genomic defense mechanism against transposons. Here, we show that piRNA-like small RNAs are present abundantly in the Bombyx ovary. We cloned 38,493 kinds of Bombyx small RNA from the ovary and performed functional characterization. Bombyx small RNAs showed a unimodal length distribution with a peak at 28nt and a strong bias for U at the 5' end. We found that 12,869 kinds of Bombyx small RNAs were associated with transposons or repetitive sequences. We classified them as repeat-associated small interfering RNAs (rasiRNAs), a subclass of piRNAs. Notably, antisense rasiRNAs have a strong bias toward U at 5' ends; in contrast, sense rasiRNAs have a strong bias toward A at nucleotide position 10, indicating that the piRNA amplification loop proposed in Drosophila is evolutionarily conserved in Bombyx. These results suggest that Bombyx small RNAs regulate transposon activity.     

Computational analysis of small RNA cloning data

Cloning and sequencing is the method of choice for small regulatory RNA identification. Using deep sequencing technologies one can now obtain up to a billion nucleotides--and tens of millions of small RNAs--from a single library. Careful computational analyses of such libraries enabled the discovery of miRNAs, rasiRNAs, piRNAs, and 21U RNAs. Given the large number of sequences that can be obtained from each individual sample, deep sequencing may soon become an alternative to oligonucleotide microarray technology for mRNA expression profiling. In this report we present the methods that we developed for the annotation and expression profiling of small RNAs obtained through large-scale sequencing. These include a fast algorithm for finding nearly perfect matches of small RNAs in sequence databases, a web-accessible software system for the annotation of small RNA libraries, and a Bayesian method for comparing small RNA expression across samples.     

The small nuclear RNAs of Drosophila

We have investigated the sequences of the major small nuclear RNAs of Drosophila cultured cells, with the objective of elucidating phylogenetically conserved primary and secondary structures by comparison of the data with previously determined sequences of these RNAs in vertebrate species. Our results reveal striking degrees of conservation between each Drosophila RNA and its vertebrate cognate, and also demonstrate blocks of homology among the Drosophila small nuclear RNAs, as previously described for vertebrates. The most conserved features include the 5' terminal region of U1 RNA, though to function in pre-mRNA splicing, most of the regions of U4 RNA recently implicated in 3' processing of pre-mRNA, and the major snRNP protein binding site ("domain A") that is also shared by vertebrate U1, U2, U4 and U5 RNAs. Several other conserved features have been revealed, suggesting additional regions of functional significance in these RNAs and also providing further insights into the evolutionary history of the small nuclear RNAs.     

Identification and remediation of biases in the activity of RNA ligases in small-RNA deep sequencing

Deep sequencing of small RNAs (sRNA-seq) is now the gold standard for small RNA profiling and discovery. Biases in sRNA-seq have been reported, but their etiology remains unidentified. Through a comprehensive series of sRNA-seq experiments, we establish that the predominant cause of the bias is the RNA ligases. We further demonstrate that RNA ligases have strong sequence-specific biases which distort the small RNA profiles considerably. We have devised a pooled adapter strategy to overcome this bias, and validated the method through data derived from microarray and qPCR. In light of our findings, published small RNA profiles, as well as barcoding strategies using adapter-end modifications, may need to be revisited. Importantly, by providing a wide spectrum of substrate for the ligase, the pooled-adapter strategy developed here provides a means to overcome issues of bias, and generate more accurate small RNA profiles.     

Intersection of Small RNA Pathways in Arabidopsis thaliana Sub-Nuclear Domains

In Arabidopsis thaliana, functionally diverse small RNA (smRNA) pathways bring about decreased RNA accumulation of target genes via several different mechanisms. Cytological experiments have suggested that the processing of microRNAs (miRNAs) and heterochromatic small interfering RNAs (hc-siRNAs) occurs within a specific nuclear domain that can present Cajal Body (CB) characteristics. It is unclear whether single or multiple smRNA-related domains are found within the same CB and how specialization of the smRNA pathways is determined within this specific sub-compartment. To ascertain whether nuclear smRNA centers are spatially related, we localized key proteins required for siRNA or miRNA biogenesis by immunofluorescence analysis. The intranuclear distribution of the proteins revealed that hc-siRNA, miRNA and trans-acting siRNA (ta-siRNA) pathway proteins accumulate and colocalize within a sub-nuclear structure in the nucleolar periphery. Furthermore, colocalization of miRNA- and siRNA-pathway members with CB markers, and reduced wild-type localization patterns in CB mutants indicates that proper nuclear localization of these proteins requires CB integrity. We hypothesize that these nuclear domains could be important for RNA silencing and may partially explain the functional redundancies and interactions among components of the same protein family. The CB may be the place in the nucleus where Dicer-generated smRNA precursors are processed and assigned to a specific pathway, and where storage, recycling or assembly of RNA interference components takes place.     

The virulent satellite RNA of turnip crinkle virus has a major domain homologous to the 3' end of the helper virus genome

RNA C (355 bases), RNA D (194 bases) and RNA F (230 bases) are small, linear satellite RNAs of turnip crinkle virus (TCV) which have been cloned as cDNAs and sequenced in this study. These RNAs produce dramatically different disease symptoms in infected plants. RNA C is a virulent satellite that intensifies virus symptoms when co-inoculated with its helper virus in turnip plants, while RNA D and RNA F are avirulent. RNA D and RNA F, the avirulent satellites, are closely related to each other except that RNA F has a 36-base insert near its 3' end, not found in RNA D. The 189 bases at the 5' end of RNA C, the virulent satellite, are homologous to the entire sequence of RNA D. However, the 3' half of RNA C, is composed of 166 bases which are nearly identical to two regions at the 3' end of the TCV helper virus genome. Hence, the virulent satellite is a composite molecule with one domain at its 5' end homologous to the other avirulent satellites and another domain at its 3' end homologous to the helper virus genome. All four TCV RNAs, RNAs C, D and F and the helper virus genome have identical 7 bases at their 3' ends. The secondary structure of RNA C deduced from the sequence can be folded into two separate domains — the domain of helper virus genome homology and the domain homologous to other TCV satellite RNAs. Comparative sequences of several different RNA C clones reveal that this satellite is a population of molecules with sequence and length heterogeneity.     

戊型肝炎病毒亚基因组RNA的研究

本文利用同位素代谢标记在HEV感染85～10.5,6.5～7.5h分别检测到1及2个亚基因组RNA,而感染21h后及在成熟的病毒颗粒内未能检测到亚基因组RNA。通过杂交实验,发现HEV的亚基因组RNA具有典型的共3′端的半套式结构,且基因组RNA与亚基因组RNA的5′端不存在共同的引导序列。通过紫外转录图谱发现HEV的亚基因组RNA是通过独立转录的方式产生的。利用引物延伸反应发现两种亚基因组RNA的转录起始位点分别位于RNA聚合酶区及非结构区、结构区的基因间序列。     

Sorting out small RNAs

Small RNAs carry out their functions by guiding Argonaute (AGO) proteins to their targets. Diverse types of small RNAs and multiple AGO proteins exist in most eukaryotic species, but how small RNAs are sorted into specific AGO complexes remains unclear. Two papers in this issue (Mi et al., 2008; Montgomery et al., 2008) now reveal the importance of the 5' terminal nucleotide of the small RNA in the sorting process in Arabidopsis.