|Home | About | Journals | Submit | Contact Us | Français|
siDirect (http://design.RNAi.jp/) is a web-based online software system for computing highly effective small interfering RNA (siRNA) sequences with maximum target-specificity for mammalian RNA interference (RNAi). Highly effective siRNA sequences are selected using novel guidelines that were established through an extensive study of the relationship between siRNA sequences and RNAi activity. Our efficient software avoids off-target gene silencing to enumerate potential cross-hybridization candidates that the widely used BLAST search may overlook. The website accepts an arbitrary sequence as input and quickly returns siRNA candidates, providing a wide scope of applications in mammalian RNAi, including systematic functional genomics and therapeutic gene silencing.
siRNA (small interfering RNA) 21–22 nt in length, the active agent in RNAi (RNA interference), serves as a guide for cognate mRNA degradation (1,2). In mammals, siRNA is expected to become a powerful tool, not only for the large-scale gene silencing essential for functional genomics, but also for therapeutic purposes, including antiviral treatments (3–5). siRNA-mediated RNAi appears to include two distinct siRNA-sequence-dependent steps: the formation of the RNA-induced silencing complex (RISC) and the recognition of target mRNA (3,6).
siRNA introduced into cells interacts with PIWI-family proteins to form RISC (7,8). Before, or at an early stage of, RISC formation, siRNA duplexes are unwound preferentially at the thermodynamically less stable end through the action of helicase. Accordingly, the two siRNA strands (sense and antisense strands) may not be incorporated into the RISC equally (9). An siRNA strand with a less stable 5′ end (i.e. a 5′ end situated in a thermodynamically less stable duplex end) is incorporated into the RISC more efficiently than the opposite strand (9–11). Consequently, in cells transfected with a highly effective siRNA, RISC with the antisense strand, may predominate. Recently, practical guidelines for selecting highly effective siRNA sequences for mammalian RNAi have been proposed (11). The rules indicate that highly effective RNAi occurs in mammalian cells and chick embryos by siRNA that satisfies the following four conditions at the same time: (i) A/U at the 5′ end of the antisense strand, (ii) G/C at the 5′ end of the sense strand, (iii) AU-richness in the 5′ terminal third of the antisense strand and (iv) the absence of any GC stretch over 9 bp in length (Figure (Figure1).1). These guidelines appear closely related to the molecular mechanism of RISC assembly as described previously (11). The first three guidelines may guarantee that the 5′ end of the antisense strand of siRNA is situated at or near the thermodynamically less stable siRNA duplex end.
The difficulty in designing siRNA sequences is the ‘off-target’ silencing effects of an siRNA duplex. Currently, we do not know how incomplete base-pairing between the antisense siRNA strand and non-targeted mRNA sequences affects silencing. Recent studies have shown that even a single mismatch in the center of an siRNA can abolish silencing, indicating the remarkable specificity of siRNA-based RNAi that permits allele-specific gene silencing (12,13). In contrast, Jackson et al. (14) showed that as few as 11 contiguous matches between siRNA and an unrelated mRNA might cause off-target silencing. These results motivated us to design siRNA sequences that contain several mismatches to all non-targeted mRNA sequences.
This paper presents a new web-based online software system, siDirect, for selecting highly effective siRNA sequences with maximal target-specificity for mammalian RNAi. The siDirect algorithm incorporates the guidelines mentioned above to favor efficient mammalian RNAi with a high success rate, and it investigates all potential cross-hybridization candidates to avoid off-target gene silencing effects.
Highly effective siRNA sequences are selected using an algorithm based on new guidelines developed by Ui-Tei et al. (11; Figure Figure1).1). Users can specify additional sequence conditions required for proper transcription initiation and termination (3,15) for designing short hairpin RNA, GC contents and custom rules.
The simple condition of requiring AA at the beginning of the target site is widely utilized to design siRNA sequences (2). This rule, however, frequently misses effective siRNA sequences that follow our four guidelines. Moreover, demanding AA at the beginning of the target site severely restricts the selection of possible siRNA sequence candidates that minimize off-target silencing effects. We therefore did not incorporate this rule into our system.
The minimization of off-target silencing effects calls for the selection of siRNA sequences that are guaranteed to have some mismatches to all unrelated sequences. We here use a rigorous specificity measure called the mismatch tolerance, the minimum number of mismatches between the siRNA sequence and any non-targeted sequence. For instance, an siRNA sequence of mismatch tolerance three does not match any off-target candidates with fewer than three mismatches. A higher mismatch tolerance of an siRNA sequence indicates its high specificity in the presence of some mismatches.
Exact mismatch-tolerance is costly to calculate, because it demands searching the entire sequence database to check whether individual siRNA oligos potentially cross-hybridize with irrelevant sequences. The Smith–Waterman local alignment algorithm (16) may return accurate answers but is very time-consuming to execute. In contrast, BLAST (17) is much faster than the Smith–Waterman algorithm, but it may overlook significant alignments.
The following alignment illustrates such an example where BLAST fails to identify the similarity between two 19 nt sequences that match with three mismatches at the 5th, 10th and 14th positions, because the two sequences do not share seven contiguous base matches, which is the shortest word (consecutive nucleotides) that BLAST requires tofind hits.
Moreover, BLAST with its default parameter values may fail to notice best alignments with minimum number of mismatches when it receives such short sequences of 19 bases as input. For instance, using BLAST, we searched the ‘nr’ database for ‘ACCGCAGTATATGGTTCTG’, a 19 nt siRNA sequence candidate for NM_000014, and we received the partial answer that the first 15 nt matched a substring of NM_002864 at 100% identity. In fact, BLAST search overlooked the following best alignment with 18/19 matches.
This search failure is due to the default parameter values of BLAST. Since ‘match reward’ and ‘mismatch penalty’ are respectively set to 1 and −3, the occurrence of one mismatch demands at least four additional matches to extend the running alignment. A partial solution to fully extend such alignments is reducing the penalty of one mismatch, which, however, is likely to output numerous, low homologous alignments with off-target candidates.
To date, most existing websites for designing siRNA sequences, such as siRNA Target Finder at the Ambion website, siDESIGN Center at Dharmacon, siRNA Target Finder at GenScript (18), and Gene specific siRNA selector (19) use BLAST to search for off-target candidates. We should bear in mind that the limitations of BLAST in seeking optimal alignments may not minimize off-target silencing effects, as illustrated in the above example. In contrast, Qiagen utilizes SSearch, a rigorous Smith–Waterman search which is computationally costly. It was the lack of efficient, accurate software for enumerating potential off-target candidates that motivated us to develop an efficient method for computing mismatch tolerance (20).
Another major issue to solve was the generation of a non-redundant sequence set of genes for checking the target specificity of siRNA sequences. Traditional non-redundant sequence datasets, such as UniGene (21) and RefSeq (21), are not suitable for this purpose, because alternative splice variants in these datasets share common exons that bring about duplication. Although searching for siRNA oligos on common exons is valuable for simultaneous silencing of all alternative splice variants, one siRNA oligo may hybridize to any of the redundant exons, calling for duplicate elimination to yield one representative exon so that siRNA oligos can be properly designed. In addition, it is also necessary to consider siRNA sequences that target the junction connecting two exons of a particular alternative splice variant. Thus, non-redundant sequences over exon–exon junctions together with duplicate-free exons ought to comprise the non-redundant sequence set of genes for checking target specificity (see Figure Figure22).
Since such a database was not available, we created one. First we aligned all the human RefSeq and Unique UniGene sequences onto the human genomic sequences (hg16). For each query sequence, we selected the best alignment that had >90% coverage ratio and >85% match ratio. We retrieved duplicate-free exons and sequences over exon–exon junctions. However, since some sequences were not totally aligned to the genomic sequence, due to sequence errors or the incompleteness of the genomic sequence, subsequences that failed to match were added to the non-redundant sequence set. The non-redundant sequence set of mouse genes was similarly generated.
The major benefit of using the non-redundant sequence set is that the mismatch tolerance of a ‘redundant sequence’, defined as a substring of more than one exon on the same locus of the genome, is likely to be higher in the non-redundant sequence set than in the original set of human RefSeq and Unique UniGene sequences. We will present the statistics below.
If an siRNA sequence is designed according to our four guidelines for effective sequences, the siRNA antisense strand is thought to be incorporated into the RISC more efficiently than the sense strand. This property may simply allow us to select effective siRNA sequences by considering only the sense target, the complement of the siRNA antisense strand, within the non-redundant sequence set. Thus, for siRNA sequences, we define, in particular, the plus-strand mismatch tolerance as mismatch tolerance calculated by using only the sense target sequence. However, we cannot disregard the possibility that the siRNA sense strand is also incorporated into the RISC and causes off-target effects. Thus, it is more reliable to take both strands of siRNA sequences into consideration. The both-strand mismatch tolerance is defined as the minimum number of mismatched bases that allow the siRNA antisense or sense strands to match a non-targeted sequence in the non-redundant sequence set. There remains the question of how much the mismatch tolerance of an siRNA sequence ought to be in order to treat the siRNA sequence target-specific.
The non-redundant sequence set of human genes was analyzed to obtain a comprehensive understanding of mismatch tolerance distribution of the 19 nt sequences that occur in the non-redundant set. The statistics for 19 nt sequences in Figure Figure3A3A shows that 9.5% are both-strand mismatch tolerance three or four but there exist no sequences of both-strand mismatch tolerance five or more. The fraction doubles if plus-strand mismatch tolerance is recalculated as illustrated in Figure Figure3B.3B. From these results, we anticipate that effective siRNA sequences of both/plus-strand mismatch tolerance three or four can be designed for most of mRNA sequences, and we define a sequence to be both-strand (plus-strand, respectively) specific if the both-strand (plus-strand) mismatch tolerance is three or more. In reality, Figure Figure3C3C shows that, for 96.3% of mRNA sequences in RefSeq, at least one effective both-strand specific siRNA sequence is designed. The fraction increases to 97.7% if the plus-strand specificity is considered instead.
Figure Figure3D3D verifies the usefulness of the non-redundant sequence set, because in the original sequence set of human RefSeq and Unique UniGene, most of redundant 19 nt sequences are both-strand mismatch tolerance zero, while, in the non-redundant sequence set, 10.7% of redundant 19 nt sequences are both-strand specific and are therefore mismatch tolerance three or more.
Figure Figure44 illustrates the flowchart of siRNA sequence selection by our system. First, our web server accepts an arbitrary sequence or an accession number to retrieve its sequence. Subsequently, the query is processed to calculate effective, gene-specific siRNA sequences by searching the non-redundant sequence set for individual 19-nt sequences. To accelerate the computational performance, we precompute all the both/plus-strand specific siRNA sequences. This precomputation makes it possible to take just a couple of seconds to return the complete list of both/plus-strand specific siRNA sequences for a typical mRNA sequence (Figure (Figure44B).
Care has to be taken to select an siRNA sequence since it may cross-react with off-target candidates. Our web server provides further information to examine the off-target silencing effects of a specific sequence. Clicking on the siRNA sequence asks the system to search the non-redundant sequence set for all the potential off-target candidates with which the siRNA sequence might cross-react. This complete search sounds to be computationally costly, but our algorithm (20) is capable of processing this request in less than a second. Subsequently, the server displays the alignment between each off-target candidate and the siRNA sequence in order to depict the locations of mismatches, which is useful in assessing the off-target silencing potential (Figure (Figure44C).
We plan to update our web server in response to major revisions to the human genome, the mouse genome, the RefSeq database and the UniGene database, though it is inevitable that such renewals will retract existing siRNA sequence candidates or add novel ones.
In reality, for ~2.3% of RefSeq sequences, we cannot design effective, plus-strand specific siRNA sequences. NM_010447, a mouse mRNA sequence, is a typical example. These cases are difficult to handle automatically. One may attempt to relax conditions on the target specificity by lowering the mismatch tolerance to two, which is likely to yield numerous effective siRNA sequence candidates, calling for some criteria for selecting more target-specific sequences from these numerous candidates.
Recall that Jackson et al. (14) stated that as few as 11 contiguous matches might cause off-target silencing. According to this observation, one promising measure would be the longest common factor of an siRNA sequence that is the length of the longest contiguous matches between the siRNA sequence and an unrelated sequence in the non-redundant set. A larger longest common factor is likely to indicate less possibility of cross-reaction with off-target candidates. Another criterion would examine the positions of mismatches. However, experimental confirmation of these criteria must be extensively performed to be utilized in practice, and hence we do not incorporate these measures into the current version of our website. Therefore, our recommendation is that users carefully select multiple siRNA sequences by investigating the list of off-target candidates (Figure (Figure4C)4C) and perform experiments to monitor off-target effects caused by individual siRNA sequences, even though the task is laborious.
We thank S. Zenno and F. Takahashi for helpful discussions and comments. The authors thank the anonymous referees for their suggestions, which have been valuable to improve our web server. This work was supported in part by funding from the Special Coordination Fund for Promoting Science and Technology to K.S., grants from the Ministry of Education, Culture, Sports, Science and Technology of Japan to K.S. and K.U.-T. and a Grant-in-Aid for Scientific Research on Priority Areas (Grant #12208003) to S.M.