|
|
||||||||
|
First published online August 18, 2006; 10.1104/pp.106.083295 Plant Physiology 142:429-440 (2006) © 2006 American Society of Plant Biologists OPEN ACCESS ARTICLE
Computational Estimation and Experimental Verification of Off-Target Silencing during Posttranscriptional Gene Silencing in Plants1,[W],[OA]Plant Biology Division, The Samuel Roberts Noble Foundation, Ardmore, Oklahoma 73401
Successful application of posttranscriptional gene silencing (PTGS) for gene function study in both plants and animals depends on high target specificity and silencing efficiency. By computational analysis with genome and/or transcriptome sequences of 25 plant species, we predicted that about 50% to 70% of gene transcripts in plants have potential off-targets when used for PTGS that could obscure experimental results. We have developed a publicly available Web-based computational tool called siRNA Scan to identify potential off-targets during PTGS. Some of the potential off-targets obtained from this tool were tested by measuring the amount of off-target transcripts using quantitative reverse transcription-PCR. Up to 50% of the predicted off-target genes tested in plants were actually silenced when tested experimentally. Our results suggest that a high risk of off-target gene silencing exists during PTGS in plants. Our siRNA Scan tool is useful to design better constructs for PTGS by minimizing off-target gene silencing in both plants and animals.
Posttranscriptional gene silencing (PTGS), also known as RNA interference (RNAi) in animals, cosuppression in plants, and RNA quelling in fungi, is an epigenetic phenomenon that results in sequence-specific degradation of endogenous mRNAs (Cogoni and Macino, 2000
In plants, PTGS can be induced by antisense and sense transgenic technology and it is also achieved by expressing dsRNA through stable or transient transformation with RNAi constructs (binary vectors) to knock down the expression of target genes (Wesley et al., 2001
Theoretically, PTGS functions in a siRNA-specific rather than a target-specific manner. However, analyses of mammalian cells transfected with different siRNAs against a target gene by two different research groups led to contradictory conclusions about silencing unintended genes (Chi et al., 2003 Investigating off-target gene silencing is crucial for accurate interpretation of gene function by PTGS and for use of PTGS application in agriculture. In this study, we estimate potential off-targets based on sequence identity for 25 plant species whose genomic or expressed sequence tag (EST) sequences are publicly available. Experimentally, we investigated the expression level of several potential off-target genes in an RNAi transgenic Arabidopsis (Arabidopsis thaliana) line and in gene-silenced (by VIGS) Nicotiana benthamiana plants. To assist the design of PTGS constructs to minimize off-target gene silencing or to identify potential off-targets from a particular PTGS construct, we have developed a Web-based computational tool. The tool provides an integrated sequence similarity search environment for plant and animal species, identifies potential off-targets, and predicts putatively effective siRNAs from the target query sequences.
Computational Analysis of Potential Off-Targets during PTGS
PTGS in plants is mediated by siRNAs derived from the cleavage of dsRNA produced through binary vector-based RNAi or VIGS constructs. The siRNAs recognize their targets by base pairing. Therefore, the potential off-targets are identified by a direct sequence identity and reverse complimentary sequence identity search. In this study, we define the cDNA producing dsRNA for silencing of the intended target gene (from which the cDNA comes) as the trigger and the unintended genes that share a contiguous
Using the full-length cDNA sequence of every transcript in the Arabidopsis genome mRNA dataset (ATH1) as a trigger, we estimated the patterns of potential off-targets by searching for identical or reverse complementary regions of
On average, 68.7% of Arabidopsis transcripts can potentially silence 3.9 off-target genes (Table I). Each candidate off-target has an average of three off-target regions with an average length of over 40 nt (Table I). However, this could be an overestimation because not all siRNAs derived from the cleavage of dsRNAs are efficient at silencing (Khvorova et al., 2003
In addition to Arabidopsis, we also estimated the potential off-targets in 24 other plant species, using Gene Indices assembled from ESTs and other expressed transcripts provided by The Institute for Genomic Research (TIGR; Quackenbush et al., 2001
Off-target silencing, however, may be desirable in PTGS applications to generate a loss-of-function phenotype if the target and off-targets are functionally redundant. For example, the same family members or genes encoding proteins that share conserved functional protein domains may be simultaneously silenced. Arabidopsis gene family information at The Arabidopsis Information Resource (TAIR) Web site (http://www.arabidopsis.org) was used to evaluate the chance that potential off-targets belong to the same gene family as the target gene. The data from the TAIR ftp site, after removing ambiguous entries and single-member families, contain 5,842 genes from 674 families (see Supplemental Table S1), with family size ranging from two to 307 members. Our analysis indicates that 4,677 gene sequences (79.9% of 5,842 genes) of these families may trigger off-target silencing of other members of the same family when a full-length sequence is used as a trigger. Among the families with different numbers of gene members, the number of off-targets was different but generally fewer than six (Fig. 2
; Supplemental Table S1). Moreover, only 418 genes from 192 families containing fewer than 10 members share an identical region of
An Integrated Search Environment for PTGS
An online tool (available at http://bioinfo2.noble.org/RNAiScan.htm) has been developed to evaluate and minimize the risk of off-target effects during PTGS. The tool is a search environment with several integrated components, including a sequence similarity search to identify potential off-targets, efficiency estimation of siRNAs, and functional analysis of off-target genes. In the sequence similarity search, the query sequence is compared with a gene transcript dataset using the BLASTn program (Altschul et al., 1997
Expression of at Least Three Off-Target Genes Was Knocked Down in BTI1 RNAi Transgenic Lines
To rule out the possibility that the lowered expression of some of these selected genes might be due to the loss of function of the BTI1 protein instead of off-target silencing, we analyzed the expression levels of all the above investigated genes in an Arabidopsis bti1-2 null mutant line (T-DNA knockout line; Salk-032220). As expected, the expression of BTI1 was below the detectable levels in the mutant plants, but the expression levels of all the other genes tested were similar to those in nontransgenic Arabidopsis plants (Fig. 3). Expression of two additional housekeeping genes, eIF4 and actin II, was also investigated for the confirmation of equal RNA amounts used for real-time qRT-PCR. Expression levels were similar among the samples from the three groups, further confirming the equal loading of total RNA in the reactions (Fig. 3). Thus, it is off-target silencing caused by BTI1-derived dsRNA that down-regulated the expression of three of the five selected potential off-targets in BTI1 RNAi transgenic plants. In addition, there is a striking phenotypic difference between BTI1 RNAi transgenic lines and the bti1-2 mutant line. Plants from four independent BTI1 RNAi transgenic lines flowered earlier than the wild-type and bti1-2 mutant plants (Fig. 4 ) and this is probably due to off-target gene silencing in the RNAi plants.
Expression of Some Potential Off-Targets Was Down-Regulated in N. benthamiana Plants Infiltrated with VIGS Constructs Computational analysis showed that more than one-half of the EST contigs (TCs) in the N. benthamiana database have multiple potential off-targets during PTGS (Table I). We randomly chose two genes for experimental analyses. They were TC381 and TC1146 encoding the U2 small nuclear ribonucleoprotein A (snRNAP A) and pyruvate decarboxylase, respectively. Partial fragments from both genes were separately cloned into the tobacco rattle virus (TRV) RNA2 vector (Liu et al., 2002 gene was used as the endogenous reaction control to normalize the relative quantity and expression of the -tubulin gene was investigated to confirm the equal loading of total RNA. The results showed a greater silencing of TC10748 than the target gene TC381 and very little silencing (not statistically significant) for CK286172 and CK288691 genes (Fig. 5A
). The analysis was repeated twice with two different groups of plants and similar results were observed.
A similar scenario was found in TRV2::TC1146-infiltrated plants. Five unique EST entries share at least 21 nt of contiguous reverse complementary sequences to the cloned TC1146 fragment (Table III). They are CK282591 encoding a protein similar to the VPI/ABI3 family regulatory protein, CK292351 encoding a Ser carboxypeptidase II-2 precursor, and CK287535, CK296810, and TC8666 encoding three different members of the auxin efflux protein family. The expression levels of these genes and TC1146 in silenced plants were compared between the silenced plants and TRV-RNA1 + TRV2:00-infiltrated plants. The expression of target gene TC1146 was decreased to about 50% in TC1146-silenced plants, whereas among the five analyzed potential off-targets, CK287535, CK296810, and TC8666 were down-regulated to a greater extent than TC1146 (Fig. 5B). These three genes all belong to the auxin efflux protein family and have the same 22-nt reverse complementary sequences to TC1146. The expression level of CK282591 was slightly reduced, whereas there was no significant change in the expression of CK292351 in TC1146-silenced plants. These experiments were repeated twice with two individual silenced plants and the results were similar. In conclusion, some of the potential off-target genes were silenced to different degrees when VIGS was used to silence target genes of interest in N. benthamiana. Because we do not have null mutations of these target genes in N. benthamiana, we cannot rule out the possibility that the reduced transcripts of the off-target genes result from the reduction in the amount of target protein.
The use of RNA silencing/interference for suppressing gene expression has become a powerful and promising approach in gene function exploration and disease treatment in both plants and animals. Its successful application relies on specific and efficient silencing of particular genes or gene families. Exquisite specificity of RNAi through siRNA in animal cells has been supported by several studies (Tuschl et al., 1999
Direct introduction of siRNA into plant protoplasts for PTGS was reported, but it is difficult to apply in intact plant tissues (Vanitharani et al., 2003
Statistically, longer siRNAs should be less likely to silence unintended genes by chance. Therefore, predicted off-target risk might be overestimated for siRNAs longer than 21 nt. However, because on average each predicted off-target for all plant species analyzed here has multiple off-target regions and each average off-target region is longer than the longest siRNAs reported in plants (Table I), the overall trend of off-target risk predicted from this study should be valid for longer siRNAs such as 22 to 26 nt. More importantly, because siRNAs of various lengths simultaneously coexist within the plant cell (Hamilton et al., 2002
Although our computational analysis showed a high risk of off-target gene silencing during PTGS in plants, the efficiency of off-target gene silencing should be the main factor that affects the functional analysis of a particular target gene. Results from our experimental analyses showed a varied reduction of expression levels of the potential off-targets in both Arabidopsis and N. benthamiana. The reduction ranges from none to a greater reduction than that of the target gene expression. The underlying mechanisms for this variation in expression levels of off-targets are not yet clear. Gene silencing efficiency is correlated with siRNA sequence-specific features (Khvorova et al., 2003
Short reverse complementary sequences (over 21 nt) were found in some completely unrelated genes in our computational analyses. In animals, siRNAs can simultaneously induce sequence-specific degradation of two endogenous mammalian transcripts oriented in opposite directions (Hu et al., 2004
The understanding of off-target silencing is crucial for accurate interpretation of gene function by PTGS. Our computational analyses with the genome and transcriptome sequences from 25 plant species showed a high risk of off-target gene silencing when a full-length sequence of each transcript entry in the datasets was used as an RNAi trigger. This off-target gene silencing risk was confirmed by our experimental analysis with both RNAi transgenic Arabidopsis lines and N. benthamiana plants infected by VIGS constructs. So far, the contribution of silenced off-targets to the silencing phenotype is not known and needs to be analyzed in exploring target gene function. For example, although the BTI1 RNAi lines had off-target gene silencing, the transformation recalcitrant phenotype of these lines is not in question because it was confirmed by null mutant lines of the target gene (Hwang and Gelvin, 2004
Datasets Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa subsp. japonica) genome mRNAs were downloaded from TIGR (http://www.tigr.org). Gene Indices for 25 plant organisms were downloaded from TIGR. The Gene Ontology (GO) database was provided by the GO Consortium (http://geneontology.org). Arabidopsis gene family information and GO annotation were downloaded from TAIR (http://www.arabidopsis.org). Protein domain data (http://pfam.wustl.edu) was also downloaded for off-target function analysis.
The gene transcript sequence dataset for each organism was searched against itself for contiguous
Arabidopsis (ecotype Columbia) and BTI1 RNAi transgenic lines (Hwang and Gelvin, 2004
The sequences of TC381 and TC1146 were obtained from the TIGR database of N. benthamiana. cDNA fragments, including the sequences with >20-nt continuous identity or complementarity to their potential off-target genes, were amplified by RT-PCR. Primers for the amplification of specific TC381 and TC1146 fragments are shown in Supplemental Table S2. About 5 µg total RNA from N. benthamiana leaves were used for RT at 42°C with NNpoly(dT)20 as the primer for 2 h. The RT product was used for PCR amplification for TC381 and TC1146 fragments, respectively, in a PTC-100 Peltier thermal cycler (M.J. Research). The PCR products were cloned into the pGEM-T-easy vector (Promega). Clones with inserts that are identical to the area in the target genes were amplified by PCR with primers that contain attb recombination sequences adapted to the previous primer pairs. PCR products were purified, sequenced, and cloned into the pTRV2 VIGS vector using the GATEWAY cloning system with the protocol from the manufacturer (Invitrogen). The recombinant plasmids were named TRV2::TC381 and TRV2::TC1146.
Agrobacterium strain GV2260 containing pTRV1, TRV2::00, TRV2::TC381, and TRV2::TC1146 were grown in an incubator at 28°C on Luria-Bertani broth with 10 µg L1 rifampicin and 50 µg L1 kanamycin for 2 d. Inoculum was prepared with the protocol published previously (Ryu et al., 2004
For Arabidopsis, total RNA was isolated from leaf tissues using TRIzol reagent (Invitrogen), followed by RNase-free DNase treatment (Promega). First-strand cDNA was synthesized with the Omniscript RT kit (Qiagen) using oligo(dT)15 according to the manufacturer's instructions. For qRT-PCR, real-time experiments were conducted in an ABI PRISM 7000 sequence detection system (Applied Biosystems), using the intercalation dye SYBR Green I as a fluorescent reporter. Quantification of PCR products was performed via a calibration curve procedure using EF-1 For N. benthamiana, total RNA from each plant at 15 d postinfiltration was extracted using TRIzol reagent with the protocol provided by the manufacturer (Molecular Research Center). First-strand cDNA was synthesized with 2 µg of total RNA using primer NNpoly(dT)20 as described above. Semiquantitative RT-PCR was performed with the modified program from the one described in "Plasmid Construction" with annealing temperature at 60°C and 20 cycles for TC381, 33 cycles for TC1146. The primer sequences are listed in Supplemental Table S2. Further real-time qRT-PCR analysis for each comparison was done with the total RNA of one TRV2::00-infiltrated and one recombinant clone-infiltrated plant, as described above for Arabidopsis samples. The analysis was repeated twice with two different groups of plants for each combination. The primer sequences used for real-time qRT-PCR are listed in Supplemental Table S3.
Data were subjected to ANOVA using Student's t test software of Excel 2003 with 98% confidence.
The following materials are available in the online version of this article.
We thank Dr. Stan Gelvin and Dr. Hau-Hsuan Hwang for sending us the Arabidopsis seeds of BTI1 RNAi lines and the bti1-2 T-DNA mutant; Dr. S.P. Dinesh-Kumar for providing GATEWAY-ready TRV-VIGS vectors; Dr. Choong-Min Ryu for help with infiltrations and discussions; and Dr. Elison Blancaflor, Dr. Stan Gelvin, and Dr. Jianzhong Liu for reviewing this manuscript. Received May 8, 2006; accepted August 7, 2006; published August 18, 2006.
1 This work was supported by The Samuel Roberts Noble Foundation.
2 These authors contributed equally to the paper. The author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (www.plantphysiol.org) is: Kirankumar S. Mysore (ksmysore{at}noble.org).
[W] The online version of this article contains Web-only data.
[OA] Open Access articles can be viewed online without a subscription. www.plantphysiol.org/cgi/doi/10.1104/pp.106.083295 * Corresponding author; e-mail ksmysore{at}noble.org; fax 5802246692.
Altschul SF, Madden TL, Schäffer A, Zhang J, Zhang Z, Miller W, Lipman DJ (1997) Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25: 33893402 Amarzguioui M, Holen T, Babaie E, Prydz H (2003) Tolerance for mutations and chemical modifications in a siRNA. Nucleic Acids Res 31: 589595 Amarzquioui M, Prydz H (2004) An algorithm for selection of functional siRNA sequences. Biochem Biophys Res Commun 316: 10501058[CrossRef][ISI][Medline] Baulcombe D (1999) Fast forward genetics based on virus-induced gene silencing. Curr Opin Plant Biol 2: 109130[CrossRef][ISI][Medline] Baulcombe D (2004) RNA silencing in plants. Nature 431: 356363[CrossRef][Medline] Birmingham A, Anderson EM, Reynolds A, Ilsley-Tyree D, Leake D, Fedorov Y, Baskerville S, Maksimova E, Robinson K, Karpilow J, et al (2006) 3' UTR seed matches, but not overall identity, are associated with RNAi off-targets. Nat Methods 3: 199204[CrossRef][ISI][Medline] Bridge AJ, Pebernard S, Ducraux A, Nicoulaz A-L, Iggo R (2003) Induction of an interferon response by RNAi vectors in mammalian cells. Nat Genet 34: 263264[CrossRef][ISI][Medline] Brown KM, Chu CY, Rana TM (2005) Target accessibility dictates the potency of human RISC. Nat Struct Mol Biol 12: 469470[CrossRef][ISI][Medline] Burch-Smith TM, Anderson JC, Martin GB, Dinesh-Kumar SP (2004) Applications and advantages of virus-induced gene silencing for gene function studies in plants. Plant J 39: 734746[CrossRef][ISI][Medline] Chi J-T, Chang HY, Wang NN, Chang DS, Dunphy N, Brown PO (2003) Genomewide view of gene silencing by small interfering RNAs. Proc Natl Acad Sci USA 100: 63436346 Chuang C-F, Meyerowitz EM (2000) Specific and heritable genetic interference by double-stranded RNA in Arabidopsis thaliana. Proc Natl Acad Sci USA 97: 49854990 Cogoni C, Macino G (2000) Post-transcriptional gene silencing across kingdoms. Curr Opin Genet Dev 10: 638643[CrossRef][ISI][Medline] Cutter AD, Payseur BA, Salcedo T, Estes AM, Good JM, Wood E, Hartl T, Maughan H, Strempel J, Wang B, et al (2003) Molecular correlates of genes exhibiting RNAi phenotypes in Caenorhabditis elegans. Genome Res 13: 26512657 Ekengren SK, Liu Y, Schiff M, Dinesh-Kumar SP, Martin GB (2003) Two MAPK cascades, NPR1, and TGA transcription factors play a role in Pto-mediated disease resistance in tomato. Plant J 36: 905917[CrossRef][ISI][Medline] Elbashir SM, Garborth J, Lendeckel W, Yalcin A, Weber K, Tuschl T (2001) Duplexes of 21-nucleotide RNAs mediate RNA interference in cultured mammalian cells. Nature 411: 494498[CrossRef][Medline] Hamilton A, Voinnet O, Chappell L, Baulcombe D (2002) Two classes of short interfering RNA in RNA silencing. EMBO J 21: 46714679[CrossRef][ISI][Medline] Hannon GJ, Rossi JJ (2004) Unlocking the potential of the human genome with RNA interference. Nature 431: 371378[CrossRef][Medline] Holen T, Amarzguioui M, Wiiger MT, Babaie E, Prydz H (2002) Positional effects of short interfering RNAs targeting the human coagulation trigger tissue factor. Nucleic Acids Res 30: 17571766 Hu X, Hipolito S, Lynn R, Abraham V, Ramos S, Wong-Staal F (2004) Relative gene-silencing efficiencies of small interfering RNAs targeting sense and antisense transcripts from the same genetic locus. Nucleic Acids Res 32: 46094617 Hwang H-H, Gelvin SB (2004) Plant proteins that interact with VirB2, the Agrobacterium tumefaciens pilin protein, mediate plant transformation. Plant Cell 16: 31483167 Jackson AL, Bartz SR, Schelter J, Kobayashi SV, Burchard J, Mao M, Li B, Cavet G, Linsley PS (2003) Expression profiling reveals off-target gene regulation by RNAi. Nat Biotechnol 21: 635637[CrossRef][ISI][Medline] Jackson AL, Linsley PS (2004) Noise amidst the silence: off-target effects of siRNAs? Trends Genet 20: 521524[CrossRef][ISI][Medline] Kamath RS, Fraser AG, Dong Y, Poulin R, Gotta M, Kanapin A, Le Bot N, Moreno S, Sohrmann M, Welchman DP, et al (2003) Systematic functional analysis of the Caenorhabditis elegans genome using RNAi. Nature 421: 231237[CrossRef][Medline] Kerschen A, Napoli CA, Jorgensen RA, Muller AE (2004) Effectiveness of RNA interference in transgenic plants. FEBS Lett 466: 223228 Khvorova A, Reynolds A, Jayasena SD (2003) Functional siRNAs and miRNAs exhibit strand bias. Cell 115: 209216[CrossRef][ISI][Medline] Kuttenkeuler D, Boutros M (2004) Genome-wide RNAi as a route to gene function in Drosophila. Brief Func Genomics Proteomics 3: 168176 Liu Y, Schiff M, Dinesh-Kumar SP (2002) Virus-induced gene silencing in tomato. Plant J 31: 777786[CrossRef][ISI][Medline] Liu Y, Schiff M, Dinesh-Kumar SP (2004) Involvement of MEK1 MAPKK, NTF6 MAPK, WRKY/MYB transcription factors, COI1 and CTR1 in N-mediated resistance to tobacco mosaic virus. Plant J 38: 800809[CrossRef][ISI][Medline] Llave C, Kasschau KD, Rector MA, Carrington JC (2002) Endogenous and silencing-associated small RNAs in plants. Plant Cell 14: 16051619 Luo KQ, Chang DC (2004) The gene-silencing efficiency of siRNA is strongly dependent on the local structure of mRNA at the targeted region. Biochem Biophys Res Commun 318: 303310[CrossRef][ISI][Medline] McManus MT, Petersen CP, Haines BB, Chen J, Sharp PA (2002) Gene silencing using micro-RNA designed hairpins. RNA 8: 842850[Abstract] Miki D, Shimamoto K (2004) Simple RNAi vectors for stable and transient suppression of gene function in rice. Plant Cell Physiol 45: 490495 Pancoska P, Moravek Z, Moll UM (2004) Efficient RNA interference depends on global context of the target sequence: quantitative analysis of silencing efficiency using Eulerian graph representation of siRNA. Nucleic Acids Res 32: 14691479 Qi Y, Denli AM, Hannon GJ (2005) Biochemical specialization within Arabidopsis RNA silencing pathways. Mol Cell 19: 421428[CrossRef][ISI][Medline] Qiu S, Adema CM, Lane T (2005) A computational study of off-target effects of RNA interference. Nucleic Acids Res 33: 18341847 Quackenbush J, Cho J, Lee D, Liang F, Holt I, Karamycheva S, Parvizi B, Pertea G, Sultana R, White J (2001) The TIGR gene indices: analysis of gene transcript sequences in highly sampled eukaryotic species. Nucleic Acids Res 29: 159164 Reynolds A, Leake D, Boese Q, Scaringe S, Marshall WS, Khvorova A (2004) Rational siRNA design for RNA interference. Nat Biotechnol 22: 326330[CrossRef][ISI][Medline] Ryu C-M, Anand A, Kang L, Mysore KS (2004) Agrodrench: a novel and effective agroinoculation method for virus-induced gene silencing in roots and diverse Solanaceous species. Plant J 40: 322331[CrossRef][ISI][Medline] Saetrom P, Snove O Jr (2004) A comparison of siRNA efficacy predictors. Biochem Biophys Res Commun 321: 247253[CrossRef][ISI][Medline] Saxena S, Jónsson SO, Dutta Z (2003) Small RNAs with imperfect match to endogenous mRNA repress translation. J Biol Chem 278: 4431244319 Scacheri PC, Rozenblatt-Rosen O, Caplen NJ, Wolfsberg TG, Umayam L, Lee JC, Hughes CM, Shanmugam KS, Bhattacharjee A, Meyerson M, et al (2004) Short interfering RNAs can induce unexpected and divergent changes in the levels of untargeted proteins in mammalian cells. Proc Natl Acad Sci USA 101: 18921897 Schwarz DS, Hutvágner G, Du T, Xu Z, Aronin N, Zamore PD (2003) Asymmetry in the assembly of the RNAi enzyme complex. Cell 115: 199208[CrossRef][ISI][Medline] Semizarov D, Frost L, Sarthy A, Kroeger P, Halbert DN, Fesik SW (2003) Specificity of short interfering RNA determined through gene expression signatures. Proc Natl Acad Sci USA 100: 63476352 Shirane D, Sugao K, Namiki S, Tanabe M, Iino M, Hirose K (2004) Enzymatic production of RNAi libraries from cDNAs. Nat Genet 36: 190196[CrossRef][ISI][Medline] Sledz CA, Holko M, de Veer MJ, Silverman RH, Williams BRG (2003) Activation of the interferon system by short-interfering RNAs. Nat Cell Biol 5: 834839[CrossRef][ISI][Medline] Sonnhammer EL, Eddy SR, Birney E, Bateman A, Durbin R (1998) Pfam: multiple sequence alignments and HMM-profiles of protein domains. Nucleic Acids Res 26: 320322 Tang G, Reinhart BJ, Bartel DP, Zamore PD (2003) A biochemical framework for RNA silencing in plants. Genes Dev 17: 4963 Thomas CL, Jones L, Baulcombe DC, Maule AJ (2001) Size constraints for targeting post-transcriptional gene silencing and for RNA-directed methylation in Nicotiana benthamiana using a potato virus X vector. Plant J 25: 417425[CrossRef][ISI][Medline] Tuschl T, Zamore PD, Lehmann R, Bartel DP, Sharp PA (1999) Targeted mRNA degradation by double-stranded RNA in vitro. Genes Dev 13: 31913197 Ui-Tei K, Naito Y, Takahashi F, Haraguchi T, Ohki-Hamazaki H, Juni A, Ueda R, Saigo K (2004) Guidelines for the selection of highly effective siRNA sequences for mammalian and chick RNA interference. Nucleic Acids Res 32: 936948 Vanitharani R, Chellappan P, Fauquet CM (2003) Short interfering RNA-mediated interference of gene expression and viral DNA accumulation in cultured plant cells. Proc Natl Acad Sci USA 100: 96329636 Wesley SV, Helliwell CA, Smith NA, Wang M, Rouse DT, Liu Q, Gooding PS, Singh SP, Abbott D, Stoutjesdijk PA, et al (2001) Construct design for efficient, effective and high-throughput gene silencing in plants. Plant J 27: 581590[CrossRef][ISI][Medline] Yokota T, Sakamoto N, Enomoto N, Tanabe Y, Miyagishi M, Maekawa S, Yi L, Kurosaki M, Taira K, Watanabe M, et al (2003) Inhibition of intracellular hepatitis C virus replication by synthetic and vector-derived small interfering RNAs. EMBO Rep 4: 602608[CrossRef][ISI][Medline] Yu JY, DeRuiter SL, Turner DL (2002) RNA interference by expression of short-interfering RNAs and hairpin | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||